- 1Henan Institute of Science and Technology, Collaborative Innovation Center of Modern Biological Breeding of Henan Province, Henan Key Laboratory Molecular Ecology and Germplasm Innovation of Cotton and Wheat, Xinxiang, China
- 2Institute of Vegetables and Flowers, Chinese Academy of Agricultural Sciences, Beijing, China
- 3State Key Laboratory of Cotton Biology, Institute of Cotton Research, Chinese Academy of Agricultural Sciences, Anyang, China
Drought seriously threats the growth and development of Gossypium hirsutum L. To dissect the genetic basis for drought tolerance in the G. hirsutum L. germplasm, a population, consisting of 188 accessions of G. hirsutum races and a cultivar (TM-1), was genotyped using the Cotton80KSNP biochip, and 51,268 high-quality single-nucleotide polymorphisms (SNPs) were obtained. Based on the phenotypic data of eight drought relative traits from four environments, we carried out association mapping with five models using GAPIT software. In total, thirty-six SNPs were detected significantly associated at least in two environments or two models. Among these SNPs, 8 and 28 (including 24 SNPs in 5 peak regions) were distributed in the A and D subgenome, respectively; eight SNPs were found to be distributed within separate genes. An SNP, TM73079, located on chromosome D10, was simultaneously associated with leaf fresh weight, leaf wilted weight, and leaf dry weight. Another nine SNPs, TM47696, TM33865, TM40383, TM10267, TM59672, TM59675, TM59677, TM72359, and TM72361, on chromosomes A13, A10, A12, A5, D6, and D9, were localized within or near previously reported quantitative trait loci for drought tolerance. Moreover, 520 genes located 200 kb up- and down-stream of 36 SNPs were obtained and analyzed based on gene annotation and transcriptome sequencing. The results showed that three candidate genes, Gh_D08G2462, Gh_A03G0043, and Gh_A12G0369, may play important roles in drought tolerance. The current GWAS represents the first investigation into mapping QTL for drought tolerance in G. hirsutum races and provides important information for improving cotton cultivars.
Introduction
Cotton is the most important fiber crop and an important oilseed crop. The genus Gossypium contains more than 50 species, of which four are cultivated. G. hirsutum L. (AD1), also known as Upland cotton, is the most widely cultivated cotton plant worldwide; it is grown on approximately 95% of cotton acreage and accounts for more than 95% of cotton production because of its high yield potential, attractive fiber quality, and wide adaptability. Even though G. hirsutum L. is adapted to subtropical climates where drought is prevalent, it is not recognized as a drought-tolerant crop and is not very efficient in terms of water use (Dabbert and Gore, 2014; Abdelraheem et al., 2019). Drought has a deleterious effect on cotton at all growth stages and overall, productivity and quality are seriously diminished. Moreover, with the increasing global population, the increasing water demand for agriculture, and the decreasing availability of fresh water, the impact of drought will be further exacerbated (Gupta et al., 2020). Therefore, determining the mechanisms underlying the response of G. hirsutum to drought stress and making full use of natural genetic variation for improving drought tolerance (DT) in cotton cultivars are urgent needs for cotton breeding programs.
Drought relative traits are complex quantitative traits that are controlled by quantitative trait loci (QTL) and environmental factors. Moreover, there is a complicated relationship between drought-tolerance, quality, and yield. Therefore, it is difficult to improve these three traits simultaneously using traditional breeding techniques. Using molecular marker techniques to identify QTL for drought relative traits and applying the molecular markers that are closely linked to, or significantly associated with, target QTL for marker-assisted selection (MAS) could rapidly produce drought-tolerant cotton varieties with high yield and excellent quality. Linkage mapping is a traditional method for gene/QTL analysis. To dissect the genetic variation underlying drought relative traits, ample QTL have been identified based on intra- and interspecific genetic linkage maps (Saeed et al., 2011; Abdelraheem et al., 2015a,b; Amjid et al., 2015; Pauli et al., 2016; Zheng et al., 2016; Pawlowicz et al., 2017; Ulloa et al., 2020). These QTL provide ample marker resources for MAS in drought-tolerant cotton breeding programs.
However, limited recombination events and low genetic diversity in the limited parental mapping population are major obstacles for conventional linkage mapping. Genome-wide association study (GWAS), as an alternative strategy, has been found to be an effective method with unique advantages for identifying the genetic variants controlling traits. Several researchers have reported markers associated with drought tolerance in G. hirsutum L. (Jia et al., 2014; Hou et al., 2018; Li B. Q. et al., 2020; Abdelraheem et al., 2021b). Jia et al. (2014) performed association mapping for DT based on 323 G. hirsutum accessions, and identified 15 significantly associated markers. Hou et al. (2018) performed association mapping with 319 Upland cotton accessions and detected 20 quantitative trait nucleotides (QTNs) significantly associated with 6 drought tolerance-related traits. Li B. Q. et al. (2020) used a natural population of 517 Upland cotton accessions to conduct association mapping and discovered 39 QTL related to drought tolerance. Abdelraheem et al. (2021a) carried out association mapping for DT and identified 53 QTL based on an association mapping panel of 376 Upland cotton accessions. Using a multiparent advanced generation intercross population of 550 recombinant inbred lines (RILs) together with their 11 Upland cotton parents, Abdelraheem et al. (2021b) carried out association mapping for drought tolerance at the seedling stage and detected 20 QTL including 13 and 7 QTL for plant height and dry shoot weight (DSW), respectively. These studies have supplemented lots of marker information for molecular marker-assisted breeding.
Although the studies above have been performed to improve the drought tolerance of cotton, the accessions they used were cultivars which were domesticated and selected repeatedly for their high yield traits. The genetic base of cotton is so narrow that it is necessary to broaden our germplasm resources for a high level of genetic diversity. G. hirsutum races, named semi-wild cotton, are the progenitors of cultivated G. hirsutum and are classified into seven races according to their geographical distribution. G. hirsutum races have wide phenotypic and high genetic diversity levels (Liu et al., 2000; Lacape et al., 2006; Abdurakhmonov et al., 2008; Hinze et al., 2015; Zhang et al., 2019). Moreover, there is no reproductive isolation between G. hirsutum races and cultivars. Many researchers have used G. hirsutum races to broaden the genetic base for improving yield (Zhang et al., 2016; Liu X. et al., 2021), fiber quality (Feng et al., 2019, 2020), cotton fiber color (Liu D. et al., 2021; Liu X. et al., 2021), and abiotic and/or biotic stress tolerance (Hinze et al., 2015; Wang et al., 2016; Xu et al., 2020). In this study, we used a natural population consisting of 188 accessions of G. hirsutum races and TM-1 to carry out a GWAS of drought tolerance based on the phenotypes in 4 different environments and the genotypes performed by the Cotton80KSNP biochip. This study aims to uncover the genetic variation in drought tolerance, to select candidate genes, and to predict the breeding potential of lines for the genetic improvement of drought-tolerance traits.
Materials and Methods
Plant Materials and Growth Conditions
The association mapping panel consisted of 189 accessions (Supplementary Table 1), which represented all seven races according to geographical distribution. Among them, 188 accessions were selected from G. hirsutum races (including 20 punctatum, 113 latifolium, 18 marie-galante, 7 palmeri, 20 morrilli, 2 yucatanense and 8 richmondi). A representative of cultivated Upland cotton, the genetic standard G. hirsutum L. acc. Texas Marker-1 (TM-1), was drought-tolerant, and was used as the control. All the accessions were perennially preserved and purified in the National Wild Cotton Nursery, Sanya, Hainan, China, which is supervised by the Institute of Cotton Research of Chinese Academy of Agricultural Sciences (ICR-CAAS), Anyang, Henan, China. The 188 accessions of G. hirsutum races were originally introduced from the USDA-ARS Southern Agricultural Research Center in College Station, Texas, USA. All the accessions were legally planted in the field/greenhouse with a single row plot during 2015–2016 and 2016–2017 in Damao, Yacheng, and Baogang, Hainan, China. In each environment, 20–25 plants were arranged in each row, with a row length of 5.0 m and a row interval of 1.0 m. Normal watering experiments were performed in the field, and drought stress experiments were carried out in the greenhouse. For all the other activities, standard local management practices were performed. For description purposes, the four environments, 2015–2016 Damao in the field, 2015–2016 Damao in the greenhouse, 2016–2017 Yacheng in the greenhouse, and 2016–2017 Baogang in the field, were designated E1, E2, E3, and E4, respectively.
Evaluation of Drought-Tolerance Traits
Eight traits contributing to drought tolerance, i.e., chlorophyll content (CHL, unit; measured using a SPAD 502 Plus (Minolta Instruments, Osaka, Japan)), leaf area (LA, cm2), plant height (PH, cm), stem diameter (SD, cm), number of leaves (LEAVES), leaf fresh weight (LFW, g plant−1), leaf wilted weight (LWW, g plant−1), and leaf dry weight (LDW, g plant−1), were investigated (Supplementary Table 2). Ten consecutive plants in the middle of each row were tagged for trait phenotyping. To reduce environmental error, we estimated the best linear unbiased prediction (BLUP) for each trait per genotype using the lme4 package in the R program. Both the BLUP and single environment values were imported for association mapping. We used R to carry out the descriptive statistical analyses, to draw the histograms for the phenotypic change trends in each environment, to perform ANOVA, and to realize the correlation analysis.
DNA Extraction and SNP Genotyping
The seeds of each accession were sown in soil with added nutrients in the greenhouse, and young leaves of each individual were collected at the seedling stage. High-quality genomic DNA was extracted by the modified CTAB method, and the concentration was diluted to 50 ng/μL. The genomic DNAs of all samples were hybridized to the CottonSNP80K biochip according to the Illumina protocols (Illumina Inc., San Diego, USA). The specifically hybridized loci were processed for single base extension, stained, and imaged on an Illumina iScan Reader. GenomeStudio Genotyping software (v1.9.4, Illumina) was used to cluster the row data for SNP genotype calls (Cai et al., 2017). The raw genotyping data were filtered, and polymorphic SNPs with a minor allele frequency (MAF) <0.05, integrity < 50%, or call rate <90% were removed.
Population Structure and Kinship Analysis
Only SNPs with a MAF ≥ 0.05, integrity ≥ 50%, and call rate ≥ 90% were used for population structure and kinship analyses. Population structure was determined using ADMIXTURE software (Alexander et al., 2009). The number of genetic clusters (K) was predefined as 1–20 to explore the population structure. This analysis provided cross-validation (CV) error estimates derived from each of the K populations, and the K with the minimum number of CV errors was taken as the number of subgroups. The kinship coefficient of each pair of individuals was calculated using Tassel 5.0 software. Linkage disequilibrium (LD) analysis was performed to determine the mapping resolution for GWAS. Pairwise LD between markers was calculated as the squared correlation coefficient (R2) of alleles using LDdecay software.
GWAS and Candidate Genes Selection
The GWAS was performed using the general linear model (GLM), mixed linear model (MLM), compressed linear mixed model (CMLM), settlement of kinship under progressively exclusive relationship (SUPER), and fixed and random model circuitous probability unification (FarmCPU), provided by the genome associated prediction integrated tool (GAPIT) program. For all of the models, the “PCA.total” was set as two. Manhattan plots were drawn using the R package qqman (Turner, 2014). The significance levels of the associations between SNPs and traits were estimated based on the threshold of the Bonferroni correction for multiple tests (1/n), where n was the total number of SNPs used in the association mapping.
To identify potential candidate genes for drought tolerance, we collected genes located within 400 kb (200 kb upstream and downstream) of significantly trait-associated SNPs based on TM-1 genome sequencing data (NBI v1.1) (Zhang et al., 2015). Then, we implemented the gene annotations from Gene Ontology (GO) analysis (http://www.geneontology.org/) (Ashburner et al., 2000).
Expression Patterns Analysis of Candidate Genes
To analyze the expression patterns of candidate genes, G. hirsutum race marie-galante85 (Mar85) and G. hirsutum race latifolium40 (Lat40) were used. The seeds of the two accessions were germinated in filter paper rolls at 28°C under 16 h light/8 h dark conditions. The seedlings were grown in Hoagland nutrient solution for 3 weeks and then transferred to 17% polyethylene glycol 6000 (PEG6000) solution. The root and leaf samples were collected at 0, 24, and 48 h after treatment. The sampling was carried out with three repeats. High-quality total RNA from each sample was extracted using TRIzol Reagent (Life Technologies, California, USA) following the manufacturer's instructions. After constructing the cDNA libraries, performing sequencing and quality control, mapping the reads to the TM-1 genome (Zhang et al., 2015), and quantifying, we calculated the fragments per kilobase of transcript per million fragments mapped reads (FPKM) values. For each accession, the normalized values represented the expression level of each candidate gene, and heatmaps of the candidate genes expression patterns were created using Mev 4.9 (Saeed et al., 2003).
Results
Phenotypic Variation and Correlation Analysis
We evaluated eight drought-tolerance traits of 189 accessions across four environments during 2015–2016 and 2016–2017 (Table 1). All eight traits displayed broad variation when grown under individual environmental condition. The average phenotypic values of the traits were as follows: CHL, 50.18 (range 35.43–66.72); LA, 130.68 cm2 (range 23.70–249.64 cm2); PH, 69.34 cm (range 26.20–146.10 cm); SD, 0.95 cm (range 0.37−2.12 cm); LEAVES, 15.115 (range 5.1–31.6); LFW, 16.02 g plant−1 (range 1.99–43.74 g); LWW, 11.35 g plant−1 (range 0.82–50.83 g); and LDW, 3.78 g plant−1 (range 0.28–12.58 g). For the eight traits, the coefficient of variation (CV) ranged from 6.87% (CHL in E2) to 43.17% (LDW in E3). As indicated by the descriptive statistics, the traits exhibited an approximately normal distribution in all environments, and the histogram of the traits investigated in each environment showed that drought-tolerance traits exhibited the genetic characteristics of quantitative traits, with continuous distributions across different environments (Supplementary Figure 1).
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/876095/fpls-13-876095-HTML/image_m/fpls-13-876095-t001.jpg)
Table 1. Descriptive statistics for the eight drought-tolerance related traits across four environments.
The correlation analysis showed that LA, LFW, LWW and LDW were significantly, highly, and positively correlated, with a correlation coefficient between 0.82 and 0.97; LEAVES showed a significant and moderate negative correlation with those four traits, with a correlation coefficient between 0.47 and 0.50, but showed a significant and moderate positive correlation with SD and PH, with a correlation coefficient between 0.50 and 0.59. Furthermore, CHL and LA, CHL and SD showed weakly negative correlations (Figure 1). The broad-sense heritability (H2) of the eight traits were estimated based on ANOVA. The H2 for PH (0.3061), SD (0.3364) and CHL (0.3771) were low (<0.5), and the H2 for LDW (0.5066), LWW (0.5418), LFW (0.5813), LEAVES (0.5901), and LA (0.5990) were moderate (> 0.5 but < 0.7) (Abdelraheem et al., 2021a). The results indicate that 30–60% of the phenotypic variation of these traits was determined by genetic factors.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/876095/fpls-13-876095-HTML/image_m/fpls-13-876095-g001.jpg)
Figure 1. Correlation analysis among the eight drought-tolerance traits of Gossypium hirsutum races. CHL, chlorophyll content; LA, leaf area; PH, plant height; SD, stem diameter; LEAVES, number of leaves; LFW, leaf fresh weight; LWW, leaf wilted weight; and LDW, leaf dry weight. *, significant at the α = 0.05 level; **, significant at the α = 0.01 level.
Distribution of Polymorphic SNPs
The 189 accessions were genotyped by CotttonSNP80K arrays, with the Illumina GenomeStudio software. In total, 51,268 filtered SNPs were detected based on a MAF ≥ 0.05, integrity ≥ 50%, and call rate ≥ 90%, and they were unevenly distributed throughout the G. hirsutum genome (Table 2; Figure 2), with 27,673 (54%) and 23,596 (46%) SNPs in the A and D subgenomes, respectively. The average marker density was approximately one SNP per 38.66 kb. In A subgenome, chromosome A8 contained the largest number of markers (4,644), with a marker density of one SNP per 22.26 kb, and chromosome A4 had the fewest markers (918), with a marker density of one SNP per 68.32 kb. For the D subgenome, chromosome D6 had the most markers (2,954), with a marker density of one SNP per 21.76 kb, and D4 had the fewest markers (939), with a marker density of one SNP per 54.63 kb. The polymorphism information content (PIC) values ranged from 0.011 (TM50634) to 0.597 (TM29059) among the SNPs, and from 0.295 (chromosome D8) to 0.341 (chromosome A8) among the chromosomes. The mean PIC values of the A and D subgenome were 0.3271 and 0.3275, respectively, and the mean PIC value of the whole genome was 0.3273.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/876095/fpls-13-876095-HTML/image_m/fpls-13-876095-g002.jpg)
Figure 2. Single-nucleotide polymorphism (SNP) distributions on the 26 chromosomes of G. hirsutum races. A1–A13 and D1–D13 on the vertical axis represent the serial numbers of the 26 chromosomes; the horizontal axis shows chromosome length; the window size is 1 Mb.
Population Structure and Relative Kinship Analysis
Population structure analysis of the 189 accessions was performed by ADMIXTURE with all 51,268 SNPs, and the results indicated that CV error was minimal at K = 2. Therefore, the optimum K was 2, and the tested population could be separated into two subgroups (Figure 3A). Subgroups 1 and 2 included 90 and 99 accessions, respectively. Most of the accessions in each subgroup had mixed ancestry, and an obvious geographic subpopulation structure was not observed, indicating that these accessions might have experienced introgression or gene flow in the complex natural environment. Moreover, 51,268 high-quality SNPs were used to evaluate the kinship coefficient. The kinship coefficient of each pair of individuals was calculated using Tassel software. A pairwise relative kinship value of 0 accounted for 62.31% of all the kinship coefficients. In addition, kinship values from 0 to 0.05 accounted for more than 70.91% of all the pairwise relative kinship coefficients (Figure 3B). Only 2.12% of the pairwise relative kinship coefficients were > 0.5. The kinship analysis indicated that most accessions had no, or a weak, relationship with each other in the population. Moreover, LD decay was measured using all 51,268 SNPs. The LD analysis estimated the average LD decay distance of our population across the whole genome to be approximately 550 kb, with an r2 = 0.5 (Figure 3C). These results indicated that the association population was suitable for association mapping.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/876095/fpls-13-876095-HTML/image_m/fpls-13-876095-g003.jpg)
Figure 3. The results of population structure, kinship coefficients, and linkage disequilibrium (LD) decay of 189 G. hirsutum accessions. (A) The total panel was divided into two subgroups, and each color represents one subgroup. (B) Histogram of the frequency distribution of pairwise relative kinship coefficients. (C) Genome-wide average LD decay estimated for the whole genome.
Genome-Wide Association Mapping
Genome-wide association mapping for the eight drought traits was performed using the GLM, MLM, CMLM, SUPER, and FarmCPU models, which were provided by the GAPIT program. SNPs with -log10P values >4.71 (1/51, 268) were selected as significantly trait-associated SNPs (Sun Z. et al., 2017). Here, based on the 5 models, 30, 43, 26, 34, and 190 SNPs randomly distributed on 26 chromosomes were associated with the 8 drought-tolerant traits based on the values for each environment, including BLUPs (Supplementary Table 3). For CHL, 12, 16, 16, 16, and 18 SNPs were identified using the GLM, MLM, CMLM, SUPER, and FarmCPU models, respectively; 4, 3, 1, 3, and 28 were associated with PH; 6, 5, 1, 3, and 14 were associated with SD; 1, 12, 3, 4, and 25 were associated with LWW; and 3, 4, 3, 4, and 41 were associated with LDW. For LA, 1, 1, and 21 SNPs were identified by GLM, SUPER, and FarmCPU; for LFW, 3, 2, 2 and 27 were identified by GLM, MLM, SUPER and FarmCPU, and for LEAVES, 1, 2, 1, and 16 SNPs were identified by MLM, cMLM, SUPER, and FarmCPU. Of the SNPs described above, 36 were significantly associated with drought-tolerant traits repeatedly detected at least in two environments or two models (Table 3; Supplementary Figure 2).
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/876095/fpls-13-876095-HTML/image_m/fpls-13-876095-t003.jpg)
Table 3. Significant SNPs associated with drought-tolerance traits that were repeatedly detected in at least two environments or two models.
Among the 36 SNPs, 24 SNPs, which were respectively associated with 5 traits, formed 5 association peaks. For CHL, an association peak, including 14 significantly associated SNPs (TM62605, TM62606, TM62607, TM62609, TM62610, TM62611, TM62612, TM62613, TM62614, TM62617, TM62621, TM62623, TM62625, and TM62626), was detected in a 165.60 kb region on chromosome D6 and was significantly associated with the E1 values. For PH, two SNPs (TM76072 and TM76073) formed a GWAS peak, which was detected at a distance of 3.4 kb on chromosome D11. These two SNPs were significantly associated with the E1 values, simultaneously detected by GLM, MLM, SUPER, and FarmCPU. For SD, three SNPs (TM74311, TM74312, and TM74314) formed an association peak, which was detected at a distance of 24.56 kb on chromosome D10. They were significantly associated with the E4 values detected by GLM, MLM, SUPER, and FarmCPU. For LWW, an association peak, including 3 SNPs (TM59672, TM59675, and TM59677), was identified at a distance of 62.75 kb on chromosome D6. These SNPs were significantly associated with the E2 values detected by the MLM, CMLM, and FarmCPU models. For LDW, an association peak including 2 SNPs (TM72359 and TM72361), was detected at a distance of 17.44 kb on chromosome D9. These SNPs were significantly associated with the E4 values detected by the GLM, MLM, CMLM, SUPER, and/or FarmCPU models.
In addition, one SNP, TM73079, located on chromosome D10, was simultaneously associated with LFW, LWW, and LDW, which could result from the pleiotropy of a single causal gene or the tight linkage of multiple causal genes.
Comparison of the Significantly Trait-Associated SNPs with those from Previous Reports
In the current study, the 36 trait-associated SNPs detected repeatedly were located on chromosomes A3, A5, A7, A9, A10, A12, A13, D5, D6, D8, D9, D10, and D11. To validate these SNPs, we investigated previously mapped markers/QTL from 23 publications on drought tolerance in cotton, and obtained 426 markers located on these 11 chromosomes. The physical locations of these markers were searched by e-PCR. Finally, nine SNPs, TM47696, TM33865, TM40383, TM10267, TM59672, TM59675, TM59677, TM72359, and TM72361, which were identified in the current study and located on chromosomes A13, A10, A12, A5, D6, and D9, were localized within or near previously reported markers/QTL for drought tolerance (Figure 4). Specifically, the marker TM47696 on chromosome A13 associated with LDW was near qDSW-A13-5 (A13_79227754—A13_79519328), which was detected by Abdelraheem et al. (2021b), with the distance between them being 90 kb. The marker TM33865 on chromosome A10 associated with LDW was near qDSW.PEG3NM.c10 (SNP0445–SNP0002), which was identified by Abdelraheem et al. (2017), with the distance between them being 810 kb. The marker TM40383 on chromosome A12 associated with LFW was near qtl_SLW_A12_E2-1 (A12_5033411 – A12_5294757), which was detected in our previous study (Magwanga et al., 2020), with the distance between them being 1.81 Mb. The marker TM10267 on chromosome A5 associated with the number of leaves was near a QTL (SNP0029–SNP0316) for LAI and qtl_SLW_A05_CA-2 (A05_7387677–A05_91982182), which were detected by Pauli et al. (2016) and/or detected in our previous research (Magwanga et al., 2020), with the distances between them being 300 kb and 2.06 Mb, respectively. The three markers, TM59672, TM59675 and TM59677, which were located on chromosome D6 and associated with LWW, were on the QTL (SNP0132–SNP0028) of LAI overlap detected by Pauli et al. (2016), and were near the two SNPs (D06_13569102 and D06_14173430) associated with DSW, which were detected by Abdelraheem et al. (2021a), with the distances between them being 1.66 Mb and 2.26 Mb, respectively. The two markers, TM72359 and TM72361 on chromosome D9, were associated with LDW, located near SNP33396, which was associated with DSW by Mahmood et al. (2019), with the distance between them being approximately 4 Mb.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/876095/fpls-13-876095-HTML/image_m/fpls-13-876095-g004.jpg)
Figure 4. Meta-analysis of markers related to drought-tolerance traits. The markers shown in red represent the SNPs detected in this study.
Distribution and Annotation of Candidate Genes
According to previous studies (Sun Z. et al., 2017; Li et al., 2018; Song et al., 2019), the genes located within ± 200 kb of the 36 significantly associated SNPs were collected based on their position in the G. hirsutum genome in this study. According to this rule, a total of 520 genes were collected (Supplementary Table 4). We then compared the distributions of these genes between A and D subgenome. The results showed that the number of genes detected within the D subgenome (275) was not significantly different from that within the A subgenome (245). On the A subgenome, chromosome A13 contained the largest numbers of genes (61), whereas chromosome A9 had only five genes. On the D subgenome, chromosome D10 had the largest number of genes (86), whereas chromosome D11 had only 27 genes. In addition, we observed the position of these candidate genes and found that eight genes included one significantly associated SNP, respectively (Table 4).
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/876095/fpls-13-876095-HTML/image_m/fpls-13-876095-t004.jpg)
Table 4. Candidate genes that include significantly associated SNPs for traits involved in drought tolerance.
To predict the functions of these 520 genes, each was annotated with GO terms. Three hundred and twenty-seven genes had annotation information, and were classified into three main types including molecular function, biological process, and cellular components (Supplementary Table 4; Supplementary Figure 3). In the molecular function category, most genes were enriched in binding, catalytic activity, transferase activity, kinase activity, oxidoreductase activity, transporter activity, and phosphatase activity. In the biological process category, 119 genes were related to stimulus response, such as water deprivation, desiccation, cadmium ion, salt stress, cold, light, wounding, bacterium, auxin, abscisic acid; 30 genes were involved in multiple signaling pathway, including ethylene, jasmonic acid, abscisic acid, auxin, and gibberellic acid mediated signaling pathway; 62, 55, and 39 genes were involved in transport, growth and developmental process, and regulation, respectively. In the cellular component category, 104 genes were mainly related to membrane and the integral component of the membrane.
Transcriptome Sequencing Analysis
To explore the response patterns of these 520 genes under drought, we carried out transcriptome sequencing using the leaf and root tissues of Mar85 and Lat40 at 0 h, 24 h and 48 h after 17% PEG6000 treatment. Among these genes, 278 and 256 genes, barely or weakly expressed (FPKM values <5) in leaf or root tissue, were removed. Moreover, the expression changes of another 43 and 49 genes in both tissues were <2-fold. The remaining 199 and 215 genes were divided into two categories in both tissues (Figures 5A–D; Supplementary Figure 4). After 24 h of drought treatment, 142 and 132 genes showed upregulation in leaf and root tissues, respectively (Figures 5A,B); while the expression of the other 57 and 83 genes displayed downregulated expression in these two tissues (Figures 5C,D).
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/876095/fpls-13-876095-HTML/image_m/fpls-13-876095-g005.jpg)
Figure 5. The expression change trends of selected genes. (A) Genes down regulated in leaf tissue. (B) Genes down regulated in root tissue. (C) Genes up regulated in leaf tissue. (D) Genes up regulated in root tissue. (E) expression pattern of drought-tolerance candidate genes in leaf tissue. (F) expression pattern of drought-tolerance candidate genes in root tissue.
Combined with the gene annotations of their homologs in Arabidopsis, 21 genes may be involved in the drought response (Figures 5E,F). Among these genes, the expression pattern of nine genes showed the same change trend in both leaf and roots, seven of which (Gh_A03G0043, Gh_A05G0473, Gh_A05G0479, Gh_A05G0494, Gh_A12G0364, Gh_A12G0371, and Gh_D10G0066, which are homologous to AT4G17615.3, AT2G25180.1, AT2G24270.3, AT5G57910.1, AT2G21180.1, AT4G38810.2, and AT1G26420.1, respectively) were upregulated, and two of which (Gh_D08G2462 and Gh_D10G0046, which are homologous to AT5G11270.1 and AT1G05850.2, respectively) were downregulated. Moreover, two genes (Gh_A10G0290 and Gh_A13G1950, which are homologous to AT1G29050.1 and AT3G09390.2, respectively) displayed opposite change trends, with decreased expression in root tissue and increased expression in leaf tissue. Of the remaining ten genes, three (Gh_A05G0490, Gh_A10G0303, and Gh_A10G0316, which are homologous to AT4G31830.1, AT1G29640.1, and AT1G29330.1) were increased in leaf tissue, one (Gh_D09G1312, which is homologous to AT3G52710.1) was increased in root tissue, two (Gh_A12G0369 and Gh_A13G1960, which are homologous to AT2G45960.1 and AT3G46970.1) was decreased in leaf tissue, and four (Gh_A13G1977, Gh_D05G1180, Gh_D06G1648, and Gh_D08G2457, which are homologous to AT2G29130.1, AT3G15810.1, AT1G62280.1, and AT3G12830.1, respectively) were decreased in the root tissue.
Discussion
Abundant Genetic Variation Among G. hirsutum Race Accessions is the Basis for Performing a GWAS
SNPs have become the preferred and widely used marker system for analyzing the genetic diversity of biological species because of the potential to identify a substantial number of SNPs throughout the whole genome. PIC is an index for evaluating the allele frequencies at single loci or their sum over multiple loci, which reflects the genetic diversity of the population. Several studies have reported the PIC of G. hirsutum based on SNPs. Based on the Cotton 63K Illumina Infinium SNP array, Sun Z. et al. (2017) genotyped a natural population of 719 Upland cotton accessions, identified 10,511 polymorphic SNPs, and calculated the mean PIC of the A and D subgenome as 0.287 and 0.283, respectively; Song et al. (2019) genotyped 279 Upland cotton accessions and calculated the mean PIC of the whole genome as 0.250 with 10,660 high-quality polymorphic SNPs. Based on the Cotton80KSNP array, Xing et al. (2019) evaluated the mean PIC of 196 Upland cotton accessions as 0.267 with 41,815 polymorphic SNPs. Using 169 cultivated Upland cotton accessions, Li et al. (2018) identified the mean PIC values of the A and D subgenome as 0.282 and 0.281, respectively, with 49,650 polymorphic SNPs. Here, we identified the mean PIC of 189 accessions of G. hirsutum as 0.3273 using 51,268 polymorphic SNPs. The results above suggest that the diversity of G. hirsutum races is greater than that of cultivated Upland cotton, which was consistent with previous reports (Liu et al., 2000; Lacape et al., 2006; Abdurakhmonov et al., 2008; Hinze et al., 2015, 2016, 2017; Wang et al., 2016; Okubazghi et al., 2017). The complex natural habitat and changeable climate where G. hirsutum races originated, which lacked strong selection pressure, played a decisive role in the accumulation of genetic variation. In contrast, Hou et al. (2018) genotyped 319 cultivars using the Cotton80KSNP array, produced 55,060 polymorphic SNPs, and calculated the mean PIC of the A and D subgenome as 0.338 and 0.334, respectively. Su et al. (2018) genotyped 355 cultivated G. hirsutum accessions through specific-locus amplified fragment sequencing, detected 93,250 polymorphic SNPs, and calculated the mean PIC as 0.34. It is undeniable that domestication can create some diversity. However, the diversity analyses are characteristically limited by the number of markers and the number of materials used. Therefore, many more accessions of G. hirsutum races will be used to evaluate the genetic diversity efficiently through genotyping by sequencing, which produces many more polymorphic SNPs, and that genetic diversity can be fully utilized in current breeding programs.
The Character of Drought-Tolerant Genetic Variation
The drought tolerance of cotton is a complex and comprehensive quantitative trait. In this study, the H2 of eight drought relative traits meant that 30% to 60% of the phenotypic variation of these traits is determined by genetic factors. Similarly, Pauli et al. (2016), Zheng et al. (2016), Abdelraheem et al. (2017, 2021a), and Li B. Q. et al. (2020) estimated the H2 for DT were 0.12–0.88, 0.40–0.666, 0.23–0.60, 0.39–0.65, and 0.08–0.76. On the contrary, Abdelraheem et al. (2015a, 2021b) estimated the H2 for DT were 0.60–0.75, and 0.681–0.691 under drought stress. The results above indicated that drought tolerance was heavily influenced by environments, and that QTL identification was urgently needed for excavating the accumulation of minor genes and major genes on DT.
In this study, based on a G. hirsutum races population with abundant genetic variation, we carried out a GWAS of eight drought-tolerant traits and identified eight and twenty-eight significantly associated SNPs in the A subgenome and D subgenome, respectively. Similarly, several previous publications indicated that more QTL for drought stresses were detected in the D subgenome than those in the A subgenome (Jia et al., 2014; Abdelraheem et al., 2015a,b, 2017, 2019, 2021a; Hou et al., 2018). Similarly, several researchers reported that the number of QTL for salt stress in the D subgenome was greater than that in the A subgenome (Du et al., 2016; Oluoch et al., 2016; Yasir et al., 2019; Yuan et al., 2019). The results above supported that the D subgenome contributed to adaptation during the long-term domestication of allotetraploid G. hirsutum (Zhang et al., 2015). However, in our previous research, the QTL for drought stress were mapped asymmetrically within the two subgenomes, including 17 and 13 QTL located in the A and D subgenomes, respectively (Magwanga et al., 2020). Similarly, Abdelraheem et al. (2021b) detected 13 and 8 QTL for DT in the A and D subgenome, respectively. The contribution of the D subgenome to drought tolerance needs to be more widely investigated based on more indicators and populations with greater genetic diversity.
Among all 36 significantly associated SNPs identified in the present study, TM73079 on chromosome D10 was simultaneously associated with LFW, LWW, and LDW. Moreover, phenotypic analysis showed that these three traits were significantly highly correlated with each other. These results may be due to the pleiotropy of a single causal gene or the tight linkage of multiple causal genes. Furthermore, another nine SNPs, TM47696, TM33865, TM40383, TM10267, TM59672, TM59675, TM59677, TM72359, and TM72361, located on chromosomes A13, A10, A12, A5, D6, and D9, were localized within or near previously reported markers/QTL for drought tolerance (Figure 4). Four SNPs, TM47696, TM33865, TM72359, and TM72361, associated with LDW were localized within QTL detected for DSW (Abdelraheem et al., 2017, 2021b); three SNPs, TM59672, TM59675, and TM59677, associated with LWW, overlapped with a QTL for LAI (Pauli et al., 2016), and were near D06_13569102 and D06_14173430, which were detected associated with DSW (Abdelraheem et al., 2021b); TM40383, associated with LFW, was near a QTL for SLW (Magwanga et al., 2020); and TM10267, associated with LEAVES, was near QTL for LAI and SLW (Pauli et al., 2016; Magwanga et al., 2020). Similarly, Li et al. (2019) detected a QTL on A11 at 86.28 Mb that associated with the plant wilting score under drought stress; this QTL overlapped with the QTL for PH at 74.64–89.93 Mb detected by (Abdelraheem et al., 2021a). The colocalized QTL found for different traits suggest that there may be a common response to drought stress and these regions may be involved in drought tolerance. Moreover, moderately or highly significant correlations (p <0.01) among LDW, LWW, LFW, LA, and LEAVES were observed in the present study, indicating that these traits have a close relationship with each other, and suggesting that they are better indicators for analyzing the genetic basis of drought tolerance in G. hirsutum.
Moreover, we compared the genomic positions of the 36 significantly associated SNPs with the relative candidate genes and subsequently identified that only eight significant SNPs were distributed within separate putative candidate genes. However, transcriptome sequencing could provide expression level of each gene and the variations in the gene coding regions. Using transcriptome-wide association studies (TWAS) to associate the target trait and the expression level of genes in the whole genome is available for predicting candidate genes in studying complex traits. Till now, TWAS has been successfully implemented in cotton and rapeseed to identify casual genes for important traits (Li Z. H. et al., 2020; Ma et al., 2021; Tang et al., 2021). For drought tolerance of G. hirsutum races, a transcriptome-wide association study (TWAS) should be confirmed to identify pivotal expression-trait association and to detect functional SNPs directly from transcriptionally active regions of the genome and potential genes.
Potential Candidate Genes of Drought Tolerance
Using the G. hirsutum population containing 189 accessions, we identified 36 significant SNPs associated with eight drought-tolerance traits, analyzed the annotation and expression pattern of 520 candidate genes, and predicted that some genes may be involved in the drought response. For example, Gh_D08G2462, containing one significantly associated SNP, TM69863 (Table 4), is homologous to AT5G11270.1 (OCP3, a member of the homeobox transcription factor family). The loss-of-function Atocp3 mutant exhibits drought tolerance through the modulation of ABA-mediated stomatal closure (Ramirez et al., 2009). In addition, Gh_A03G0043 is homologous to AT4G17615.3 (CBL1), which codes for the calcineurin B-like (CBL) protein; this protein is a primary calcium sensor and plays crucial roles in the response to drought in plants (Cheong et al., 2003). Here the gene expression of Gh_A03G0043 increased in leaf and root tissues after drought treatment, which is similar to our previous report on wild type Arabidopsis (Lu et al., 2019; Magwanga et al., 2019), and the expression of its homologous genes are similar in several plants, such as Arabidopsis (Cheong et al., 2003), P. betulifolia Bunge (Xu et al., 2015), and C. sinensis (Shu et al., 2020). Moreover, several reports indicated that overexpression of CBL1 may confer tolerance to drought (Cheong et al., 2003). Another candidate gene, Gh_A12G0369, is an ortholog of AT2G45960.1 (plasma membrane intrinsic protein 1, PIP1B), which encodes a plasma membrane intrinsic protein and functions as a water channel in plants. The gene expression of Gh_A12G0369 decreased in both tissues after drought treatment in this study, which is consistent with the results in A. thaliana (Jang et al., 2004), F. arundinacea (Pawlowicz et al., 2017), and S. bigelovii (Sun X. et al., 2017). Many studies have indicated that PIP1B was involved in drought stress. For example, overexpression of SbPIP1 in tobacco plants improved drought tolerance (Sun X. et al., 2017); overexpression of NnPIP1-2, ScPIP1, and SpPIP1, enhanced the tolerance of Arabidopsis to stresses (Chen et al., 2018; Wang et al., 2019; Liu et al., 2020). Overexpression of AtPIP1;2 in transgenic tobacco plants increased plant vigor under favorable growth conditions but had no beneficial effects under drought stress (Aharon et al., 2003). However, the decrease of PIP1b gene expression in Arabidopsis resulted in reduced water permeability in leaf mesophyll protoplasts and an increased root:shoot ratio (Kaldenhoff et al., 1998). Meanwhile, Gh_A05G0479 is homologous to AT2G24270.3 (aldehyde dehydrogenase, ALDH11A3), Gh_A13G1950 is homologous to AT3G09390.2 (metallothionein 2a, MT2A), and Gh_A05G0473 is homologous to AT2G25180.1 (Response Regulators 12, RR12). These homologous genes in Arabidopsis may be involved in the drought response. These results suggest that Gh_D08G2462, Gh_A03G0043, and Gh_A12G0369 may be candidate genes involved in the drought response, and future research will be directed to understanding their regulatory mechanisms.
Data Availability Statement
The original contributions presented in the study are included in the article/Supplementary Materials, further inquiries can be directed to the corresponding author/s.
Author Contributions
XC and YX designed the experiment. XG, YuaW, YH, RS, TQ, and ZH analyzed the data. XG, YH, ZZ, KW, FL, YuhW, XC, and YX carried out the experiments. XG, YuaW, and YH prepared the manuscript. All authors have read, discussed, and approved the current version of the manuscript.
Funding
This work was supported by the Science and Technology Development Project of Henan Province (212102110054 and 222102110146), the Key Scientific Research Projects of Higher Education of Henan Province (22A210014), State Key Laboratory of Cotton Biology Open Fund (CB2021A02), and the Postgraduate Education Reform and Quality Improvement Project of Henan Province (Yu degree [2018] No. 23).
Conflict of Interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Publisher's Note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
Supplementary Material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fpls.2022.876095/full#supplementary-material
Supplementary Figure 1. Histogram of drought-tolerance traits investigated in four environments. E1, E2, E3, and E4 indicate the four environments 2015–2016 Damao in the field, 2015–2016 Damao in the greenhouse, 2016–2017 Yacheng in the greenhouse, and 2016–17 Baogang in the field, respectively; CHL, chlorophyll content; LA, leaf area; PH, plant height; SD, stem diameter; LEAVES, number of leaves; LFW, leaf fresh weight; LWW, leaf wilted weight; LDW, leaf dry weight.
Supplementary Figure 2. Manhattan plot of genome-wide association (GWA) mapping results for drought-tolerance traits in G. hirsutum races. 1–26 on the horizontal axis represent chromosome A1-A13, D2, D1, D7, D3, D13, D5, D10, D11, D4, D9, D8, D6, and D12, respectively; (A) Chlorophyll content (CHL); (B) Leaf area (LA); (C) Plant height (PH); (D) Stem diameter (SD); (E) Number of leaves (LEAVES); (F) Leaf fresh weight (LFW); (G) Leaf wilted weight (LWW); (H) Leaf dry weight (LDW); the dashed line represents the significance threshold (-log10P > 4.71).
Supplementary Figure 3. GO annotation of the genes within 200 kb upstream and downstream of the significantly associated SNPs.
Supplementary Figure 4. Heatmap of the expression levels of 199 and 215 genes in leaf and root tissues of MAR85 and LAT40 after drought treatment. (A) The 142 genes showed up-regulation in leaf tissue; (B) The 132 genes showed up-regulation in root tissue; (C) The 57 genes showed down-regulation in leaf tissue; (D) The 83 genes showed down-regulation in root tissue. Red indicates high expression, and green indicates low expression.
Supplementary Table 1. List of G. hirsutum accessions used in this study.
Supplementary Table 2. Data acquisition for eight traits in four environments. “Y” represented that data were collected; “-“ indicated that data were missing.
Supplementary Table 3. SNPs significantly associated with drought-tolerance that were detected using five models with GAPIT software. E1, E2, E3, and E4 indicate the four environments 2015–2016 Damao in the field, 2015–2016 Damao in the greenhouse, 2016–2017 Yacheng in the greenhouse, and 2016–2017 Baogang in the field, respectively. CHL, chlorophyll content; LA, leaf area; PH, plant height; SD, stem diameter; LEAVES, number of leaves; LFW, leaf fresh weight; LWW, leaf wilted weight; and LDW, leaf dry weight. GLM, general linear model; MLM, mixed linear model; CMLM, compressed linear mixed model; SUPER, settlement of kinship under progressively exclusive relationship; FarmCPU, fixed and random model circuitous probability unification.
Supplementary Table 4. Annotation of 520 genes identified based on 36 significantly associated SNPs.
References
Abdelraheem, A., Esmaeili, N., O'Connell, M., and Zhang, J. F. (2019). Progress and perspective on drought and salt stress tolerance in cotton. Ind. Crops Prod. 130, 118–129. doi: 10.1016/j.indcrop.2018.12.070
Abdelraheem, A., Fang, D. D., and Zhang, J. F. (2017). Quantitative trait locus mapping of drought and salt tolerance in an introgressed recombinant inbred line population of upland cotton under the greenhouse and field conditions. Euphytica. 17, 214. doi: 10.1007/s10681-017-2095-x
Abdelraheem, A., Hughs, S. E., Jones, D. C., Zhang, J. F., and Jenkins, J. (2015a). Genetic analysis and quantitative trait locus mapping of PEG-induced osmotic stress tolerance in cotton. Plant Breed. 134, 111–120. doi: 10.1111/pbr.12228
Abdelraheem, A., Kuraparthy, V., Hinze, L., Stelly, D., Wedegaertner, T., and Zhang, J. F. (2021a). Genome-wide association study for tolerance to drought and salt tolerance and resistance to thrips at the seedling growth stage in US Upland upland cotton. Ind. Crops Prod. 169. doi: 10.1016/j.indcrop.2021.113645
Abdelraheem, A., Mahdy, E., and Zhang, J. F. (2015b). The first linkage map for a recombinant inbred line population in cotton (Gossypium barbadense) and its use in studies of PEG-induced dehydration tolerance. Euphytica. 205, 941–958. doi: 10.1007/s10681-015-1448-6
Abdelraheem, A., Thyssen, G. N., Fang, D. D., Jenkins, J. N., McCarty, J. C., Wedegaertner, T., et al. (2021b). GWAS reveals consistent QTL for drought and salt tolerance in a MAGIC population of 550 lines derived from intermating of 11 Upland cotton (Gossypium hirsutum) parents. Mol. Genet. Genomics. 296, 119–129. doi: 10.1007/s00438-020-01733-2
Abdurakhmonov, I. Y., Kohel, R. J., Yu, J. Z., Pepper, A. E., Abdullaev, A. A., Kushanov, F. N., et al. (2008). Molecular diversity and association mapping of fiber quality traits in exotic G. hirsutum L. germplasm. Genomics. 92, 478–487. doi: 10.1016/j.ygeno.2008.07.013
Aharon, R., Shahak, Y., Wininger, S., Bendov, R., Kapulnik, Y., and Galili, G. (2003). Overexpression of a plasma membrane aquaporin in transgenic tobacco improves plant vigor under favorable growth conditions but not under drought or salt stress. Plant Cell. 15, 439–447. doi: 10.1105/tpc.009225
Alexander, D. H., Novembre, J., and Lange, K. (2009). Fast model-based estimation of ancestry in unrelated individuals. Genome Res. 19, 1655–1664. doi: 10.1101/gr.094052.109
Amjid, M. W., Malik, T. A., Shakeel, A., and Wahid, A. (2015). QTL mapping for relative leaf water contents, cell Membrane stability and excised leaf water loss under drought by using EST-SSR markers in Gossypium hirsutum. Int. J. Agric. Biol. 17, 779–784. doi: 10.17957/IJAB/14.0011
Ashburner, M., Ball, C. A., Blake, J. A., Botstein, D., Butler, H., Cherry, J. M., et al. (2000). Gene ontology: tool for the unification of biology. the gene ontology consortium. Nat Genet. 25, 25–29. doi: 10.1038/75556
Cai, C. P., Zhu, G. Z., Zhang, T. Z., and Guo, W. Z. (2017). High-density 80 K SNP array is a powerful tool for genotyping G. hirsutum accessions and genome analysis. BMC Genomics. 18, 654. doi: 10.1186/s12864-017-4062-2
Chen, Q., Yang, S. H., Kong, X. X., Wang, C. T., Xiang, N., Yang, Y. Q., et al. (2018). Molecular cloning of a plasma membrane aquaporin in Stipa purpurea, and exploration of its role in drought stress tolerance. Gene. 665, 41–48. doi: 10.1016/j.gene.2018.04.056
Cheong, Y. H., Kim, K. N., Pandey, G. K., Gupta, R., Grant, J. J., and Luan, S. (2003). CBL1, a calcium sensor that differentially regulates salt, drought, and cold responses in Arabidopsis. Plant Cell. 15, 1833–1845. doi: 10.1105/tpc.012393
Dabbert, T. A., and Gore, M. A. (2014). Challenges and perspectives on improving heat and drought stress resilience in cotton. J. Cotton Sci. 18, 393–409.
Du, L., Cai, C. P., Wu, S., Zhang, F., Hou, S., and Guo, W. Z. (2016). Evaluation and exploration of favorable QTL alleles for salt stress related traits in cotton cultivars (G. hirsutum L.). PLoS One. 11, e0151076. doi: 10.1371/journal.pone.0151076
Feng, L. C., Zhang, S. W., Xing, L. T., Yang, B., Gao, X., Xie, X. M., et al. (2019). QTL analysis for yield and fibre quality traits using three sets of introgression lines developed from three Gossypium hirsutum race stocks. Mol. Genet. Genom. 294, 789–810. doi: 10.1007/s00438-019-01548-w
Feng, L. C., Zhou, C. H., Su, Q., Xu, M., Yue, H. R., Zhang, S. W., et al. (2020). Fine-mapping and candidate gene analysis of qFS-Chr. D02, a QTL for fibre strength introgressed from a semi-wild cotton into Gossypium hirsutum. Plant Sci. 297, 110524. doi: 10.1016/j.plantsci.2020.110524
Gupta, A., Rico-Medina, A., and Cano-Delgado, A. I. (2020). The physiology of plant responses to drought. Science. 368, 266–269. doi: 10.1126/science.aaz7614
Hinze, L. L., Fang, D. D., Gore, M. A., Scheffler, B. E., Yu, J. Z., Frelichowski, J., et al. (2015). Molecular characterization of the Gossypium diversity reference set of the US national cotton germplasm collection. Theor. Appl. Genet. 128, 313–327. doi: 10.1007/s00122-014-2431-7
Hinze, L. L., Gazave, E., Gore, M. A., Fang, D. D., Scheffler, B. E., Yu, J. Z., et al. (2016). Genetic diversity of the two commercial tetraploid cotton species in the Gossypium diversity reference set. J. Hered. 107, 274–286. doi: 10.1093/jhered/esw004
Hinze, L. L., Hulse-Kemp, A. M., Wilson, I. W., Zhu, Q. H., Llewellyn, D. J., Taylor, J. M., et al. (2017). Diversity analysis of cotton (Gossypium hirsutum L.) germplasm using the CottonSNP63K Array. BMC Plant. Biol. 17, 37. doi: 10.1186/s12870-017-0981-y
Hou, S., Zhu, G. Z., Li, Y., Li, W. X., Fu, J., Niu, E. L., et al. (2018). Genome-wide association studies reveal genetic variation and candidate genes of drought stress related traits in cotton (Gossypium hirsutum L.). Front Plant Sci. 9, 1276. doi: 10.3389/fpls.2018.01276
Jang, J. Y., Kim, D. G., Kim, Y. O., Kim, J. S., and Kang, H. (2004). An expression analysis of a gene family encoding plasma membrane aquaporins in response to abiotic stresses in Arabidopsis thaliana. Plant Mol. Biol. 54, 713–725. doi: 10.1023/B:PLAN.0000040900.61345.a6
Jia, Y. H., Sun, J. L., Wang, X. W., Zhou, Z. L., Pan, Z. E., He, S. P., et al. (2014). Molecular diversity and association analysis of drought and salt tolerance in Gossypium hirsutum L. Germplasm. J. Integr. Agric. 13, 1845–1853. doi: 10.1016/S2095-3119(13)60668-1
Kaldenhoff, R., Grote, K., Zhu, J. J., and Zimmermann, U. (1998). Significance of plasmalemma aquaporins for water-transport in Arabidopsis thaliana. Plant J. 14, 121–128. doi: 10.1046/j.1365-313X.1998.00111.x
Lacape, J. M., Dessauw, D., Rajab, M., Noyer, J. L., and Hau, B. (2006). Microsatellite diversity in tetraploid Gossypium germplasm: assembling a highly informative genotyping set of cotton SSRs. Mol. Breed. 19, 45–58. doi: 10.1007/s11032-006-9042-1
Li, B. Q., Tian, Q., Wang, X. W., Han, B., Liu, L., ong, X. H., et al. (2020). Phenotypic plasticity and genetic variation of cotton yield and its related traits under water-limited conditions. Crop J. 8, 966–976. doi: 10.1016/j.cj.2020.02.003
Li, C. Q., Wang, Y. Y., Ai, N. J., Li, Y., and Song, J. F. (2018). A genome-wide association study of early-maturation traits in upland cotton based on the CottonSNP80K array. J. Integr. Plant Biol. 60, 970–985. doi: 10.1111/jipb.12673
Li, H. M., Liu, S. D., Ge, C. W., Zhang, X. M., Zhang, S. P., Chen, J., et al. (2019). Analysis of drought tolerance and associated traits in Upland cotton at the seedling stage. Int. J. Mol. Sci. 20. doi: 10.3390/ijms20163888
Li, Z. H., Wang, P. C., You, C. Y., Yu, J. W., Zhang, X. N., Yan, F. L., et al. (2020). Combined GWAS and eQTL analysis uncovers a genetic regulatory network orchestrating the initiation of secondary cell wall development in cotton. New Phytol. 226, 1738–1752. doi: 10.1111/nph.16468
Liu, D., Liu, X., Su, Y., Zhang, X., Guo, K., Teng, Z., et al. (2021). Genetic mapping and identification of Lg loci controlling green fuzz in Upland cotton (Gossypium hirsutum L.). Crop J. 9, 777–784. doi: 10.1016/j.cj.2020.09.009
Liu, S., Cantrell, R. G., McCarty, J. C., and Stewart, J. M. (2000). Simple sequence repeat-based assessment of genetic diversity in cotton race stock accessions. Crop Sci. 40, 1459–1469. doi: 10.2135/cropsci2000.4051459x
Liu, X., Yang, L., Wang, J., Wang, Y., Guo, Z., Li, Q., et al. (2021). Analyzing quantitative trait loci for fiber quality and yield-related traits from a recombinant inbred line population with Gossypium hirsutum race palmeri as one parent. Front. Plant Sci. 12, 817748. doi: 10.3389/fpls.2021.817748
Liu, Z. Y., Wang, C. F., Yang, J. J., Liu, X., Li, L. J., Cheng, L. B., et al. (2020). Molecular cloning and functional analysis of lotus salt-induced NnDREB2C, NnPIP1-2 and NnPIP2-1 in Arabidopsis thaliana. Mol. Biol. Rep. 47, 497–506. doi: 10.1007/s11033-019-05156-0
Lu, P., Magwanga, R. O., Kirungu, J. N., Hu, Y., Dong, Q., Cai, X., et al. (2019). Overexpression of cotton a DTX/MATE gene enhances drought, salt, and cold stress tolerance in transgenic Arabidopsis. Front. Plant Sci. 10, 299. doi: 10.3389/fpls.2019.00299
Ma, Y. Z., Min, L., Wang, J. D., Li, Y. Y., Wu, Y. L., Hu, Q., et al. (2021). A combination of genome-wide and transcriptome-wide association studies reveals genetic elements leading to male sterility during high temperature stress in cotton. New Phytol. 231, 165–181. doi: 10.1111/nph.17325
Magwanga, R. O., Kirungu, J. N., Lu, P., Yang, X., Dong, Q., Cai, X., et al. (2019). Genome wide identification of the trihelix transcription factors and overexpression of Gh_A05G2067 (GT-2), a novel gene contributing to increased drought and salt stresses tolerance in cotton. Physiol Plant. 167, 447–464. doi: 10.1111/ppl.12920
Magwanga, R. O., Lu, P., Kirungu, J. N., Cai, X., Zhou, Z., Agong, S. G., et al. (2020). Identification of QTLs and candidate genes for physiological traits associated with drought tolerance in cotton. J. Cotton Res. 3, 43. doi: 10.1186/s42397-020-0043-0
Mahmood, T., Khalid, S., Abdullah, M., Ahmed, Z., Shah, M. K. N., Ghafoor, A., et al. (2019). Insights into drought stress signaling in plants and the molecular genetic basis of cotton drought tolerance. Cells. 9, 105. doi: 10.3390/cells9010105
Okubazghi, K. W., Li, X. N., Cai, X. Y., Wang, X. X., Chen, H. D., Zhou, Z. L., et al. (2017). Genome-wide assessment of genetic diversity and fiber quality traits characterization in Gossypium hirsutum races. J. Integr. Agric. 16, 2402–2412. doi: 10.1016/S2095-3119(17)61671-X
Oluoch, G., Zheng, J., Wang, X., Khan, M. K. R., Zhou, Z., Cai, X., et al. (2016). QTL mapping for salt tolerance at seedling stage in the interspecific cross of Gossypium tomentosum with Gossypium hirsutum. Euphytica. 209, 223–235. doi: 10.1007/s10681-016-1674-6
Pauli, D., Andrade-Sanchez, P., Carmo-Silva, A. E., Gazave, E., French, A. N., Heun, J., et al. (2016). Field-based high-throughput plant phenotyping reveals the temporal patterns of quantitative trait loci associated with stress-responsive traits in cotton. G3 (Bethesda). 6, 865–879. doi: 10.1534/g3.115.023515
Pawlowicz, I., Rapacz, M., Perlikowski, D., Gondek, K., and Kosmala, A. (2017). Abiotic stresses influence the transcript abundance of PIP and TIP aquaporins in Festuca species. J. Appl. Genet. 58, 421–435. doi: 10.1007/s13353-017-0403-8
Ramirez, V., Coego, A., Lopez, A., Agorio, A., Flors, V., and Vera, P. (2009). Drought tolerance in Arabidopsis is controlled by the OCP3 disease resistance regulator. Plant J. 58, 578–591. doi: 10.1111/j.1365-313X.2009.03804.x
Saeed, A. I., Sharov, V., White, J., Li, J., Liang, W., Bhagabati, N., et al. (2003). TM4: a free, open-source system for microarray data management and analysis. Biotechniques. 34, 374–378. doi: 10.2144/03342mt01
Saeed, M., Guo, W. Z., Ullah, I., Tabbasam, N., Zafar, Y., Ur-Rahman, M., et al. (2011). QTL mapping for physiology, yield and plant architecture traits in cotton (Gossypium hirsutum L.) grown under well-watered versus water-stress conditions. Electron J. Biotechnol. 14, 3. doi: 10.2225/vol14-issue3-fulltext-3
Shu, B., Cai, D., Zhang, F., Zhang, D. J., Liu, C. Y., Wu, Q. S., et al. (2020). Identifying citrus CBL and CIPK gene families and their expressions in response to drought and arbuscular mycorrhizal fungi colonization. Biol plant. 64, 773–783. doi: 10.32615/bp.2020.123
Song, C. X., Li, W., Pei, X. Y., Liu, Y. G., Ren, Z. Y., He, K. L., et al. (2019). Dissection of the genetic variation and candidate genes of lint percentage by a genome-wide association study in upland cotton. Theor. Appl. Genet. 132, 1991–2002. doi: 10.1007/s00122-019-03333-0
Su, J. J., Li, L. B., Zhang, C., Wang, C. X., Gu, L. J., Wang, H. T., et al. (2018). Genome-wide association study identified genetic variations and candidate genes for plant architecture component traits in Chinese upland cotton. Theor. Appl. Genet. 131, 1299–1314. doi: 10.1007/s00122-018-3079-5
Sun, X., Deng, Y., Liang, L., Jia, X., Xiao, Z., and Su, J. (2017). Overexpression of a PIP1 gene from Salicornia bigelovii in tobacco plants improves their drought tolerance. J. Am. Soc. Hortic. Sci. 142, 235–245. doi: 10.21273/JASHS04098-17
Sun, Z., Wang, X., Liu, Z., Gu, Q., Zhang, Y., Li, Z., et al. (2017). Genome-wide association study discovered genetic variation and candidate genes of fibre quality traits in Gossypium hirsutum L. Plant Biotechnol. J. 15, 982–996. doi: 10.1111/pbi.12693
Tang, S., Zhao, H., Lu, S. P., Yu, L. Q., Zhang, G. F., Zhang, Y. T., et al. (2021). Genome- and transcriptome-wide association studies provide insights into the genetic basis of natural variation of seed oil content in Brassica napus. Mol. Plant. 14, 470–487. doi: 10.1016/j.molp.2020.12.003
Turner, S. D. (2014). qqman: an R package for visualizing GWAS results using Q-Q and manhattan plots. bioRxiv preprint. doi: 10.1101/005165
Ulloa, M., De Santiago, L. M., Hulse-Kemp, A. M., Stelly, D. M., and Burke, J. J. (2020). Enhancing Upland cotton for drought resilience, productivity, and fiber quality: comparative evaluation and genetic dissection. Mol. Genet. Genom. 295, 155–176. doi: 10.1007/s00438-019-01611-6
Wang, X., Gao, F., Bing, J., Sun, W., Feng, X., Ma, X., et al. (2019). Overexpression of the jojoba aquaporin gene, ScPIP1, enhances drought and salt tolerance in transgenic Arabidopsis. Int. J. Mol. Sci. 20, 153. doi: 10.3390/ijms20010153
Wang, Y. Y., Zhou, Z. L., Wang, X. X., Cai, X. Y., Li, X. N., Wang, C. Y., et al. (2016). Genome-wide association mapping of glyphosate-resistance in Gossypium hirsutum races. Euphytica. 209, 209–221. doi: 10.1007/s10681-016-1663-9
Xing, H. X., Yuan, Y. C., Zhang, H. J., Wang, L. Y., Mao, L. L., Tao, J. C., et al. (2019). Multi-environments and multi-models association mapping identified candidate genes of lint percentage and seed index in Gossypium hirsutum L. Mol Breed. 39, 7. doi: 10.1007/s11032-019-1063-7
Xu, Y., Magwanga, R. O., Yang, X., Jin, D., Cai, X., Hou, Y., et al. (2020). Genetic regulatory networks for salt-alkali stress in Gossypium hirsutum with differing morphological characteristics. BMC Genomics. 21, 15. doi: 10.1186/s12864-019-6375-9
Xu, Y. Y., Li, H., Lin, J., Li, X. G., Chang, Y. H. (2015). Isolation and characterization of Calcineurin B-like gene (PbCBL1) and its promoter in birch-leaf pear (Pyrus betulifolia Bunge). Genet. Mol. Res. 14, 16756–16770. doi: 10.4238/2015.December.14.3
Yasir, M., He, S., Sun, G., Geng, X., Pan, Z., Gong, W., et al. (2019). A genome-wide association study revealed key SNPs/Genes associated with salinity stress tolerance in Upland cotton. Genes (Basel). 10, 829. doi: 10.3390/genes10100829
Yuan, Y., Xing, H., Zeng, W., Xu, J., Mao, L., Wang, L., et al. (2019). Genome-wide association and differential expression analysis of salt tolerance in Gossypium hirsutum L at the germination stage. BMC Plant Biol. 19, 394. doi: 10.1186/s12870-019-1989-2
Zhang, S., Cai, Y., Guo, J., Li, K., Peng, R., Liu, F., et al. (2019). Genotyping-by-sequencing of Gossypium hirsutum races and cultivars uncovers novel patterns of genetic relationships and domestication footprints. Evol Bioinform Online. 15, 48. doi: 10.1177/1176934319889948
Zhang, S. W., Feng, L. C., Xing, L. T., Yang, B., Gao, X., Zhu, X. F., et al. (2016). New QTLs for lint percentage and boll weight mined in introgression lines from two feral landraces into Gossypium hirsutum acc TM-1. Plant Breed. 135, 90–101. doi: 10.1111/pbr.12337
Zhang, T., Hu, Y., Jiang, W., Fang, L., Guan, X., Chen, J., et al. (2015). Sequencing of allotetraploid cotton (Gossypium hirsutum L. acc. TM-1) provides a resource for fiber improvement. Nat. Biotechnol. 33, 531–537. doi: 10.1038/nbt.3207
Keywords: cotton, G. hirsutum race, drought, genome-wide association study, Cotton80KSNP, GAPIT, candidate gene prediction
Citation: Guo X, Wang Y, Hou Y, Zhou Z, Sun R, Qin T, Wang K, Liu F, Wang Y, Huang Z, Xu Y and Cai X (2022) Genome-Wide Dissection of the Genetic Basis for Drought Tolerance in Gossypium hirsutum L. Races. Front. Plant Sci. 13:876095. doi: 10.3389/fpls.2022.876095
Received: 15 February 2022; Accepted: 31 May 2022;
Published: 28 June 2022.
Edited by:
Ting Peng, Henan Agricultural University, ChinaReviewed by:
Chengqi Li, Yuncheng University, ChinaAwais Rasheed, Quaid-i-Azam University, Pakistan
Copyright © 2022 Guo, Wang, Hou, Zhou, Sun, Qin, Wang, Liu, Wang, Huang, Xu and Cai. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Yanchao Xu, xuyanchao@foxmail.com; Xiaoyan Cai, cxycri@163.com
† These authors have contributed equally to this work and share first authorship