Skip to main page content
U.S. flag

An official website of the United States government

Dot gov

The .gov means it’s official.
Federal government websites often end in .gov or .mil. Before sharing sensitive information, make sure you’re on a federal government site.

Https

The site is secure.
The https:// ensures that you are connecting to the official website and that any information you provide is encrypted and transmitted securely.

Access keys NCBI Homepage MyNCBI Homepage Main Content Main Navigation
. 2023 Dec 20;25(1):65.
doi: 10.3390/ijms25010065.

Integrated Bulk Segregant Analysis, Fine Mapping, and Transcriptome Revealed QTLs and Candidate Genes Associated with Drought Adaptation in Wild Watermelon

Affiliations

Integrated Bulk Segregant Analysis, Fine Mapping, and Transcriptome Revealed QTLs and Candidate Genes Associated with Drought Adaptation in Wild Watermelon

Ahmed Mahmoud et al. Int J Mol Sci. .

Abstract

Drought stress has detrimental effects on crop productivity worldwide. A strong root system is crucial for maintaining water and nutrients uptake under drought stress. Wild watermelons possess resilient roots with excellent drought adaptability. However, the genetic factors controlling this trait remain uninvestigated. In this study, we conducted a bulk segregant analysis (BSA) on an F2 population consisting of two watermelon genotypes, wild and domesticated, which differ in their lateral root development under drought conditions. We identified two quantitative trait loci (qNLR_Dr. Chr01 and qNLR_Dr. Chr02) associated with the lateral root response to drought. Furthermore, we determined that a small region (0.93 Mb in qNLR_Dr. Chr01) is closely linked to drought adaptation through quantitative trait loci (QTL) validation and fine mapping. Transcriptome analysis of the parent roots under drought stress revealed unique effects on numerous genes in the sensitive genotype but not in the tolerant genotype. By integrating BSA, fine mapping, and the transcriptome, we identified six genes, namely L-Ascorbate Oxidase (AO), Cellulose Synthase-Interactive Protein 1 (CSI1), Late Embryogenesis Abundant Protein (LEA), Zinc-Finger Homeodomain Protein 2 (ZHD2), Pericycle Factor Type-A 5 (PFA5), and bZIP transcription factor 53-like (bZIP53-like), that might be involved in the drought adaptation. Our findings provide valuable QTLs and genes for marker-assisted selection in improving water-use efficiency and drought tolerance in watermelon. They also lay the groundwork for the genetic manipulation of drought-adapting genes in watermelon and other Cucurbitacea species.

Keywords: drought tolerance; genetic resources; genomics-assisted breeding; quantitative trait loci (QTL) mapping; transcriptome; wild watermelon.

PubMed Disclaimer

Conflict of interest statement

The authors declare no conflict of interest.

Figures

Figure 1
Figure 1
Drought tolerance index (DTI) based on the number of lateral roots (NLR) in 38 watermelon accessions exposed to 15% Polyethylene Glycol (PEG, drought) and distilled water (CK, control) for 4 days in pouches. The DTI was calculated as the ratio between NLR with drought and CK treatments. The tested genotypes were categorized into five groups based on their DTI. The plotted values are averages of three independent biological replications (nine plants/replication/treatment).
Figure 2
Figure 2
Root phenotyping and validation of drought tolerance in the parents, F1 and F2 individuals. (a) Root phenotypes of ZJU076 (tolerant), ZJU196 (sensitive), and F1 after 4 days of growth with 15% PEG6000 (drought, Dr) and distilled water (control, CK). (b) Drought tolerance index (DTI) of the parents and F1. Values are means ± S.D. (n = 18). (c) After 15 days of water withholding (the final soil water content ≈ is 2%), (d) Frequency distribution of the number of lateral roots (as drought tolerance indicator) among 484 F2 individuals exposed to 15% Polyethylene Glycol (PEG) for four days.
Figure 3
Figure 3
Bulk segregant analysis (BSA) results of the drought tolerance in watermelon. (a) The single nucleotide polymorphism (SNP) index of the high pool, (b) the SNP index of the low pool, and (c) the delta SNP index values used for the association analysis. The x and y axes show the 11 chromosomes of watermelon and the SNP index, respectively. The curved line indicates the fitted SNP index or delta SNP index. The horizontal line indicates the association threshold with FFN of a 95% confidence interval. (d) Major QTLs for drought tolerance in watermelon detected by the G prime (G’) method. The QTL1, and QTL2 represent qNLR_Dr. Chr01 and qNLR_Dr. Chr02, respectively.
Figure 4
Figure 4
Validation of the detected QTLs through haplotype analysis using six Kompetitive allele specific PCR (KASP) markers for each QTL in 305 F2 individuals. Green and red indicate homozygous segments from sensitive and tolerant parents, respectively. Yellow indicates heterozygous segments. The sidebars (on the left) represent the root traits of F2 individuals: primary root length (PRL), total root system (TRS), lateral root system (LRS), and the number of lateral roots (NLR). The single nucleotide polymorphism (SNP) data are arranged in ascending order from top to bottom based on the NLR values. Zooming in will make the sample IDs and root trait values much more visible.
Figure 5
Figure 5
Fine mapping of the detected QTLs with 12 Kompetitive allele specific PCR (KASP) markers. (a) Parent genotypes and two F2 recombinants in the target regions. The left part represents qNLR_Dr. Chr01, while the right part shows qNLR_Dr. Chr02. (b) The recombinant offspring individuals of F2 are sorted based on qNLR_Dr. Chr02. (c) The recombinant offspring individuals of F2 are sorted based on qNLR_Dr. Chr01. The F2 offspring were classified into two groups based on the segment origins. Red indicates the homozygous ZJU076 (tolerant) segment, green indicates the homozygous ZJU196 (sensitive) segment, yellow shows the heterozygous region, and gray represents mixed (involved individuals similar to the parents or heterozygous). The average number of lateral roots (NLR) of each family was calculated from ten individuals.
Figure 6
Figure 6
The significant differentially expressed genes (DEGs) within the delimited region of qNLR_Dr. Chr01. Log2 (fold change) of the 13 DEGs (p < 0.05) in CK_76 vs. DR_76 and CK_ 196 vs. D.R._ 196 within the delimited region of qNLR_Dr. Chr01.

Similar articles

Cited by

References

    1. Regan P.M., Kim H., Maiden E. Climate change, adaptation, and agricultural output. Reg. Environ. Chang. 2018;19:113–123. doi: 10.1007/s10113-018-1364-0. - DOI
    1. FAOSTAT Food and Agriculture Organization of the United Nations. [(accessed on 16 August 2022)]. Available online: http://faostat.fao.org.
    1. Mall M., Kumar R., Akhtar M.Q. Chapter 1—Horticultural crops and abiotic stress challenges. In: Chandra Rai A., Rai A., Kumar Rai K., Rai V.P., Kumar A., editors. Stress Tolerance in Horticultural Crops. Woodhead Publishing; Sawston, UK: 2021. pp. 1–19.
    1. Li H., Mo Y., Cui Q., Yang X., Guo Y., Wei C., Yang J., Zhang Y., Ma J., Zhang X. Transcriptomic and physiological analyses reveal drought adaptation strategies in drought-tolerant and -susceptible watermelon genotypes. Plant Sci. 2019;278:32–43. doi: 10.1016/j.plantsci.2018.10.016. - DOI - PubMed
    1. Fahad S., Bajwa A.A., Nazir U., Anjum S.A., Farooq A., Zohaib A., Sadia S., Nasim W., Adkins S., Saud S., et al. Crop Production under Drought and Heat Stress: Plant Responses and Management Options. Front. Plant Sci. 2017;8:1147. doi: 10.3389/fpls.2017.01147. - DOI - PMC - PubMed