Genetic Diversity of the Two Commercial Tetraploid Cotton Species in the Gossypium Diversity Reference Set

Regulation of a single inositol 1-phosphate synthase homeologue by HSFA6B contributes to fibre yield maintenance under drought conditions in upland cotton

Drought stress substantially impacts crop physiology resulting in alteration of growth and productivity. Understanding the genetic and molecular crosstalk between stress responses and agronomically important traits such as fibre yield is particularly complicated in the allopolyploid species, upland cotton (Gossypium hirsutum), due to reduced sequence variability between A and D subgenomes. To better understand how drought stress impacts yield, the transcriptomes of 22 genetically and phenotypically diverse upland cotton accessions grown under well-watered and water-limited conditions in the Arizona low desert were sequenced. Gene co-expression analyses were performed, uncovering a group of stress response genes, in particular transcription factors GhDREB2A-A and GhHSFA6B-D, associated with improved yield under water-limited conditions in an ABA-independent manner. DNA affinity purification sequencing (DAP-seq), as well as public cistrome data from Arabidopsis, were used to identify targets of these two TFs. Among these targets were two lint yield-associated genes previously identified through genome-wide association studies (GWAS)-based approaches, GhABP-D and GhIPS1-A. Biochemical and phylogenetic approaches were used to determine that GhIPS1-A is positively regulated by GhHSFA6B-D, and that this regulatory mechanism is specific to Gossypium spp. containing the A (old world) genome. Finally, an SNP was identified within the GhHSFA6B-D binding site in GhIPS1-A that is positively associated with yield under water-limiting conditions. These data lay out a regulatory connection between abiotic stress and fibre yield in cotton that appears conserved in other systems such as Arabidopsis.

First Report on the Genetic Diversity of Populations of Gossypium barbadense L. and Gossypium hirsutun L. in the Amazonian Native Communities, Cusco-Peru

The genus Gossypium has important ethnobotanical and economic value for Amazonian Native Communities (A.N.C.). However, little research has been undertaken on the distribution and genetic diversity of cotton populations maintained in the Peruvian rainforest. This work aims to present the first report on the genetic diversity of Gossypium spp. populations in the A.N.C. of the province of La Convención, Cusco-Peru. The methodology was based on exploring, collecting, identifying, and characterizing the Gossypium populations present in the A.N.C. Twenty-six descriptors were evaluated (9 quantitative and 17 qualitative), and with this information, distribution, correlation, and principal component (PC) analyses were carried out. As a result, plants of two species [G. barbadense L. (44 samples) and G. hirsutum L. (19 samples)], one variety [G. barbadense var. brasiliensis (75 samples)], and three previously unidentified variations (9 samples) were identified. Altogether, 147 samples were collected. G. barbadense var. brasiliensis, which was always found in association with other economic crops within an altitude range of 338 to 1086 m, was the most predominant (51%), distributed in eleven A.N.C. and always in small plots (up to 2 ha). G. barbadense L. was cultivated between 397 and 1137 m of altitude in eight A.N.C. in plots of up to 3 ha in marginal lands. G. hirsutum L., with a smaller distribution (13%), was found between 334 and 497 m of altitude in only three communities; this species is cultivated in marginal areas throughout the year. The variability found for the first two PCs when considering the quantitative and qualitative descriptors was high (74.7%) and moderate (48.2%), respectively. When combining all the descriptors, the analysis showed that the first two PCs accounted for 51.8% of the total variability of the data. The PCs of the two types of data and their combination confirmed that the three populations found were grouped. The nine undefined samples were close to or intermediate between the described ones, showing that these samples may be the result of spontaneous crosses; as such, these samples need to be better evaluated with other tools for further definition. The information obtained shows that in the A.N.C. of Cusco-Peru, there is variability conserved by the inhabitants, who have been able to maintain and use these genotypes, even from their Amazonian indigenous ancestry, and the environment has been able to generate variability among the species, as will be highlighted in future works.

Genetic polymorphism detection in brazilian perennial cottons (Gossypium spp.) using an ISSR marker system and its application for molecular interspecific differentiation

BACKGROUND

The semi-domesticated Brazilian perennial cotton (Gossypium spp.) germplasm is considered a source of variability for creating modern upland cotton varieties. Here we used Inter-simple Sequence Repeat (ISSR) markers to detect intra and interspecific genetic polymorphism in Gossypium hirsutum L. r. marie-galante and Gossypium barbadense L. and to use molecular data to assessing genetic diversity and molecular discrimination of these species.

METHODS AND RESULTS

The sets contained 12 G. barbadense genotypes and 16 G. hirsutum genotypes from a Brazilian collection. The 11 ISSR primers were used for genotyping yielded 101 bands (polymorphism = 47.5%) and were classified as moderately informative (PIC = 0.304). The ISSR markers exposed a greater diversity in G. hirsutum (P = 24.72%; H_E =0.071 and I = 0.111) as compared to G. barbadense (P = 17.98%, H_E = 0.043 and I = 0.070). The AMOVA analysis showed that 89.47% of the genetic variation was partitioned within species which is supported by Nei's genetic differentiation (G_st = 0.598) and gene flow (N_m = 0.338), suggesting that strong reproductive barriers between species. The UPGMA Cluster Analysis, Principal Coordinate Analysis and Bayesian Model-Based Structural Analysis divided the 28 genotypes into two main clades consistent with the taxonomical delimitation.

CONCLUSION

The ISSR marker system offers a new approach to determining molecular differences between two cotton species (G. hirsutum L. r. marie-galante and G. barbadense L.). This study can expand the molecular marker resources for the identification and improvement of our knowledge about the genetic diversity and relationships between perennial cotton genotypes.

The final piece of the Triangle of U: Evolution of the tetraploid Brassica carinata genome

Ethiopian mustard (Brassica carinata) is an ancient crop with remarkable stress resilience and a desirable seed fatty acid profile for biofuel uses. Brassica carinata is one of six Brassica species that share three major genomes from three diploid species (AA, BB, and CC) that spontaneously hybridized in a pairwise manner to form three allotetraploid species (AABB, AACC, and BBCC). Of the genomes of these species, that of B. carinata is the least understood. Here, we report a chromosome scale 1.31-Gbp genome assembly with 156.9-fold sequencing coverage for B. carinata, completing the reference genomes comprising the classic Triangle of U, a classical theory of the evolutionary relationships among these six species. Our assembly provides insights into the hybridization event that led to the current B. carinata genome and the genomic features that gave rise to the superior agronomic traits of B. carinata. Notably, we identified an expansion of transcription factor networks and agronomically important gene families. Completion of the Triangle of U comparative genomics platform has allowed us to examine the dynamics of polyploid evolution and the role of subgenome dominance in the domestication and continuing agronomic improvement of B. carinata and other Brassica species.

Genome-Wide Dissection of the Genetic Basis for Drought Tolerance in Gossypium hirsutum L. Races

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.

Cotton pan-genome retrieves the lost sequences and genes during domestication and selection

BACKGROUND

Millennia of directional human selection has reshaped the genomic architecture of cultivated cotton relative to wild counterparts, but we have limited understanding of the selective retention and fractionation of genomic components.

RESULTS

We construct a comprehensive genomic variome based on 1961 cottons and identify 456 Mb and 357 Mb of sequence with domestication and improvement selection signals and 162 loci, 84 of which are novel, including 47 loci associated with 16 agronomic traits. Using pan-genome analyses, we identify 32,569 and 8851 non-reference genes lost from Gossypium hirsutum and Gossypium barbadense reference genomes respectively, of which 38.2% (39,278) and 14.2% (11,359) of genes exhibit presence/absence variation (PAV). We document the landscape of PAV selection accompanied by asymmetric gene gain and loss and identify 124 PAVs linked to favorable fiber quality and yield loci.

CONCLUSIONS

This variation repertoire points to genomic divergence during cotton domestication and improvement, which informs the characterization of favorable gene alleles for improved breeding practice using a pan-genome-based approach.

Identification and characterization of genes related to salt stress tolerance within segregation distortion regions of genetic map in F2 population of upland cotton

Segregation distortion (SD) is a genetic mechanism commonly found in segregating or stable populations. The principle behind this puzzles many researchers. The F₂ generation developed from wild Gossypium darwinii and G. hirsutum CCRI12 species was used to investigate the possible transcription factors within the segregation distortion regions (SDRs). The 384 out of 2763 markers were distorted in 29 SDRs on 18 chromosomes. Good collinearity was observed among genetic and physical maps of G. hirsutum and G. barbadense syntenic blocks. Total 568 genes were identified from SDRs of 18 chromosomes. Out of these genes, 128 belonged to three top-ranked salt-tolerant gene families. The DUF597 contained 8 uncharacterized genes linked to Pkinase (PF00069) gene family in the phylogenetic tree, while 15 uncharacterized genes clustered with the zinc finger gene family. Two hundred thirty four miRNAs targeted numerous genes, including ghr-miR156, ghr-miR399 and ghr-miR482, while others targeted top-ranked stress-responsive transcription factors. Moreover, these genes were involved in the regulation of numerous stress-responsive cis-regulatory elements. The RNA sequence data of fifteen upregulated genes were verified through the RT-qPCR. The expression profiles of two highly upregulated genes (Gh_D01G2015 and Gh_A01G1773) in salt-tolerant G. darwinii showed antagonistic expression in G. hirsutum. The results indicated that salt-tolerant genes have been possibly transferred from the wild G. darwinii species. A detailed functional analysis of these genes can be carried out which might be helpful in the future for gene cloning, transformation, gene editing and the development of salt-resistant cotton varieties.

Networks of Physiological Adjustments and Defenses, and Their Synergy With Sodium (Na+) Homeostasis Explain the Hidden Variation for Salinity Tolerance Across the Cultivated Gossypium hirsutum Germplasm

The abilities to mobilize and/or sequester excess ions within and outside the plant cell are important components of salt-tolerance mechanisms. Mobilization and sequestration of Na+ involves three transport systems facilitated by the plasma membrane H+/Na+ antiporter (SOS1), vacuolar H+/Na+ antiporter (NHX1), and Na+/K+ transporter in vascular tissues (HKT1). Many of these mechanisms are conserved across the plant kingdom. While Gossypium hirsutum (upland cotton) is significantly more salt-tolerant relative to other crops, the critical factors contributing to the phenotypic variation hidden across the germplasm have not been fully unraveled. In this study, the spatio-temporal patterns of Na+ accumulation along with other physiological and biochemical interactions were investigated at different severities of salinity across a meaningful genetic diversity panel across cultivated upland Gossypium. The aim was to define the importance of holistic or integrated effects relative to the direct effects of Na+ homeostasis mechanisms mediated by GhHKT1, GhSOS1, and GhNHX1. Multi-dimensional physio-morphometric attributes were investigated in a systems-level context using univariate and multivariate statistics, randomForest, and path analysis. Results showed that mobilized or sequestered Na+ contributes significantly to the baseline tolerance mechanisms. However, the observed variance in overall tolerance potential across a meaningful diversity panel were more significantly attributed to antioxidant capacity, maintenance of stomatal conductance, chlorophyll content, and divalent cation (Mg2+) contents other than Ca2+ through a complex interaction with Na+ homeostasis. The multi-tier macro-physiological, biochemical and molecular data generated in this study, and the networks of interactions uncovered strongly suggest that a complex physiological and biochemical synergy beyond the first-line-of defense (Na+ sequestration and mobilization) accounts for the total phenotypic variance across the primary germplasm of Gossypium hirsutum. These findings are consistent with the recently proposed Omnigenic Theory for quantitative traits and should contribute to a modern look at phenotypic selection for salt tolerance in cotton breeding.

Construction of Gossypium barbadense Mutant Library Provides Genetic Resources for Cotton Germplasm Improvement

Allotetraploid cotton (Gossypium hirsutum and Gossypium barbadense) are cultivated worldwide for its white fiber. For centuries, conventional breeding approaches increase cotton yield at the cost of extensive erosion of natural genetic variability. Sea Island cotton (G. barbadense) is known for its superior fiber quality, but show poor adaptability as compared to Upland cotton. Here, in this study, we use ethylmethanesulfonate (EMS) as a mutagenic agent to induce genome-wide point mutations to improve the current germplasm resources of Sea Island cotton and develop diverse breeding lines with improved adaptability and excellent economic traits. We determined the optimal EMS experimental procedure suitable for construction of cotton mutant library. At M₆ generation, mutant library comprised of lines with distinguished phenotypes of the plant architecture, leaf, flower, boll, and fiber. Genome-wide analysis of SNP distribution and density in yellow leaf mutant reflected the better quality of mutant library. Reduced photosynthetic efficiency and transmission electron microscopy of yellow leaf mutants revealed the effect of induced mutations at physiological and cellular level. Our mutant collection will serve as the valuable resource for basic research on cotton functional genomics, as well as cotton breeding.

Genetic Basis of Fiber Improvement and Decreased Stress Tolerance in Cultivated Versus Semi-Domesticated Upland Cotton

Crop domestication from wild ancestors has resulted in the wide adaptation coupled with improved yield and quality traits. However, the genetic basis of many domesticated characteristics remains to be explored. Upland cotton (Gossypium hirsutum) is the most important tetraploid cotton species, accounting for about 90% of world cotton commerce. Here, we reveal the effects of domestication on fiber and stress traits through comprehensive analyses of semi-domesticated races and cultivated cotton accessions. A total of 416 cotton accessions were genotyped, and a decrease in genetic diversity from races to landraces and modern cultivars was detected. Furthermore, 71 domestication selective sweeps (DSS) and 14 improvement selective sweeps (ISS) were identified, with the Dt sub-genome experiencing stronger selection than the At sub-genome during the both selection types. The more expressed genes and a delay in the expression peak of genes related to secondary cell wall (SCW) development in modern cultivars compared to semi-domesticated cotton races, may have contributed to long fibers in these plants. However, down-regulation of genes related to stress response was responsible for decreasing stress tolerance in modern cultivars. We further experimentally confirmed that silencing of PR1 and WRKY20, genes that showed higher expression in the semi-domesticated races, drastically compromised cotton resistance to V. dahliae. Our results reveal fiber improvement and decreased stress tolerance as a result of the domestication of modern upland cotton cultivars.

SSR-Linkage map of interspecific populations derived from Gossypium trilobum and Gossypium thurberi and determination of genes harbored within the segregating distortion regions

Wild cotton species have significant agronomic traits that can be introgressed into elite cultivated varieties. The use of a genetic map is important in exploring, identification and mining genes which carry significant traits. In this study, 188 F₂mapping individuals were developed from Gossypium thurberi (female) and Gossypium trilobum (male), and were genotyped by using simple sequence repeat (SSR) markers. A total of 12,560 simple sequence repeat (SSR) markers, developed by Southwest University, thus coded SWU were screened out of which only 994 were found to be polymorphic, and 849 markers were linked in all the 13 chromosomes. The map had a length of 1,012.458 cM with an average marker distance of 1.193 cM. Segregation distortion regions (SDRs) were observed on Chr01, Chr02, Chr06, Chr07 Chr09, Chr10 and Chr11 with a large proportion of the SDR regions segregating towards the heterozygous allele. There was good syntenic block formation that revealed good collinearity between the genetic and physical map of G. raimondii, compared to the Dt_sub genome of the G. hirsutum and G. barbadense. A total of 2,496 genes were mined within the SSR related regions. The proteins encoding the mined genes within the SDR had varied physiochemical properties; their molecular weights ranged from 6.586 to 252.737 kDa, charge range of -39.5 to 52, grand hydropathy value (GRAVY) of -1.177 to 0.936 and isoelectric (pI) value of 4.087 to 12.206. The low GRAVY values detected showed that the proteins encoding these genes were hydrophilic in nature, a property common among the stress responsive genes. The RNA sequence analysis revealed more of the genes were highly upregulated in various stages of fiber development for instance; Gorai.002G241300 was highly up regulated at 5, 10, 20 and 25 day post anthesis (DPA). Validation through RT-qPCR further revealed that these genes mined within the SDR regions might be playing a significant role under fiber development stages, therefore we infer that Gorai.007G347600 (TFCA), Gorai.012G141600 (FOLB1), Gorai.006G024500 (NMD3), Gorai.002G229900 (LST8) and Gorai.002G235200 (NSA2) are significantly important in fiber development and in turn the quality, and further researches needed to be done to elucidate their exact roles in the fiber development process. The construction of the genetic map between the two wild species paves away for the mapping of quantitative trait loci (QTLs) since the average distance between the markers is small, and mining of genes on the SSR regions will provide an insight in identifying key genes that can be introgressed into the cultivated cotton cultivars.

Assessment of Genetic Diversity, Population Structure, and Evolutionary Relationship of Uncharacterized Genes in a Novel Germplasm Collection of Diploid and Allotetraploid Gossypium Accessions Using EST and Genomic SSR Markers

This study evaluated the genetic diversity and population structures in a novel cotton germplasm collection comprising 132 diploids, including Glossypium klotzschianum and allotetraploid cotton accessions, including Glossypium barbadense, Glossypium darwinii, Glossypium tomentosum, Glossypium ekmanianum, and Glossypium stephensii, from Santa Cruz, Isabella, San Cristobal, Hawaiian, Dominican Republic, and Wake Atoll islands. A total of 111 expressed sequence tag (EST) and genomic simple sequence repeat (gSSR) markers produced 382 polymorphic loci with an average of 3.44 polymorphic alleles per SSR marker. Polymorphism information content values counted 0.08 to 0.82 with an average of 0.56. Analysis of a genetic distance matrix revealed values of 0.003 to 0.53 with an average of 0.33 in the wild cotton collection. Phylogenetic analysis supported the subgroups identified by STRUCTURE and corresponds well with the results of principal coordinate analysis with a cumulative variation of 45.65%. A total of 123 unique alleles were observed among all accessions and 31 identified only in G. ekmanianum. Analysis of molecular variance revealed highly significant variation between the six groups identified by structure analysis with 49% of the total variation and 51% of the variation was due to diversity within the groups. The highest genetic differentiation among tetraploid populations was observed between accessions from the Hawaiian and Santa Cruz regions with a pairwise FST of 0.752 (p < 0.001). DUF819 containing an uncharacterized gene named yjcL linked to genomic markers has been found to be highly related to tryptophan-aspartic acid (W-D) repeats in a superfamily of genes. The RNA sequence expression data of the yjcL-linked gene Gh_A09G2500 was found to be upregulated under drought and salt stress conditions. The existence of genetic diversity, characterization of genes and variation in novel germplasm collection will be a landmark addition to the genetic study of cotton germplasm.

Genetic diversity, linkage disequilibrium, and association mapping analyses of Gossypium barbadense L. germplasm

Limited polymorphism and narrow genetic base, due to genetic bottleneck through historic domestication, highlight a need for comprehensive characterization and utilization of existing genetic diversity in cotton germplasm collections. In this study, 288 worldwide Gossypium barbadense L. cotton germplasm accessions were evaluated in two diverse environments (Uzbekistan and USA). These accessions were assessed for genetic diversity, population structure, linkage disequilibrium (LD), and LD-based association mapping (AM) of fiber quality traits using 108 genome-wide simple sequence repeat (SSR) markers. Analyses revealed structured population characteristics and a high level of intra-variability (67.2%) and moderate interpopulation differentiation (32.8%). Eight percent and 4.3% of markers revealed LD in the genome of the G. barbadense at critical values of r2 ≥ 0.1 and r2 ≥ 0.2, respectively. The LD decay was on average 24.8 cM at the threshold of r2 ≥ 0.05. LD retained on average distance of 3.36 cM at the threshold of r2 ≥ 0.1. Based on the phenotypic evaluations in the two diverse environments, 100 marker loci revealed a strong association with major fiber quality traits using mixed linear model (MLM) based association mapping approach. Fourteen marker loci were found to be consistent with previously identified quantitative trait loci (QTLs), and 86 were found to be new unreported marker loci. Our results provide insights into the breeding history and genetic relationship of G. barbadense germplasm and should be helpful for the improvement of cotton cultivars using molecular breeding and omics-based technologies.

Investigation of the Influence of Leaf Thickness on Canopy Reflectance and Physiological Traits in Upland and Pima Cotton Populations

Many systems for field-based, high-throughput phenotyping (FB-HTP) quantify and characterize the reflected radiation from the crop canopy to derive phenotypes, as well as infer plant function and health status. However, given the technology's nascent status, it remains unknown how biophysical and physiological properties of the plant canopy impact downstream interpretation and application of canopy reflectance data. In that light, we assessed relationships between leaf thickness and several canopy-associated traits, including normalized difference vegetation index (NDVI), which was collected via active reflectance sensors carried on a mobile FB-HTP system, carbon isotope discrimination (CID), and chlorophyll content. To investigate the relationships among traits, two distinct cotton populations, an upland (Gossypium hirsutum L.) recombinant inbred line (RIL) population of 95 lines and a Pima (G. barbadense L.) population composed of 25 diverse cultivars, were evaluated under contrasting irrigation regimes, water-limited (WL) and well-watered (WW) conditions, across 3 years. We detected four quantitative trait loci (QTL) and significant variation in both populations for leaf thickness among genotypes as well as high estimates of broad-sense heritability (on average, above 0.7 for both populations), indicating a strong genetic basis for leaf thickness. Strong phenotypic correlations (maximum r = -0.73) were observed between leaf thickness and NDVI in the Pima population, but not the RIL population. Additionally, estimated genotypic correlations within the RIL population for leaf thickness with CID, chlorophyll content, and nitrogen discrimination ([Formula: see text] = -0.32, 0.48, and 0.40, respectively) were all significant under WW but not WL conditions. Economically important fiber quality traits did not exhibit significant phenotypic or genotypic correlations with canopy traits. Overall, our results support considering variation in leaf thickness as a potential contributing factor to variation in NDVI or other canopy traits measured via proximal sensing, and as a trait that impacts fundamental physiological responses of plants.

Diversity analysis of cotton (Gossypium hirsutum L.) germplasm using the CottonSNP63K Array

BACKGROUND

Cotton germplasm resources contain beneficial alleles that can be exploited to develop germplasm adapted to emerging environmental and climate conditions. Accessions and lines have traditionally been characterized based on phenotypes, but phenotypic profiles are limited by the cost, time, and space required to make visual observations and measurements. With advances in molecular genetic methods, genotypic profiles are increasingly able to identify differences among accessions due to the larger number of genetic markers that can be measured. A combination of both methods would greatly enhance our ability to characterize germplasm resources. Recent efforts have culminated in the identification of sufficient SNP markers to establish high-throughput genotyping systems, such as the CottonSNP63K array, which enables a researcher to efficiently analyze large numbers of SNP markers and obtain highly repeatable results. In the current investigation, we have utilized the SNP array for analyzing genetic diversity primarily among cotton cultivars, making comparisons to SSR-based phylogenetic analyses, and identifying loci associated with seed nutritional traits.

RESULTS

The SNP markers distinctly separated G. hirsutum from other Gossypium species and distinguished the wild from cultivated types of G. hirsutum. The markers also efficiently discerned differences among cultivars, which was the primary goal when designing the CottonSNP63K array. Population structure within the genus compared favorably with previous results obtained using SSR markers, and an association study identified loci linked to factors that affect cottonseed protein content.

CONCLUSIONS

Our results provide a large genome-wide variation data set for primarily cultivated cotton. Thousands of SNPs in representative cotton genotypes provide an opportunity to finely discriminate among cultivated cotton from around the world. The SNPs will be relevant as dense markers of genome variation for association mapping approaches aimed at correlating molecular polymorphisms with variation in phenotypic traits, as well as for molecular breeding approaches in cotton.