Duplication,divergence and persistence in the Phytochrome photoreceptor gene family of cottons (<Emphasis Type="Italic">Gossypium</Emphasis> spp.)

Background  

Phytochromes are a family of red/far-red photoreceptors that regulate a number of important developmental traits in cotton (Gossypium spp.), including plant architecture, fiber development, and photoperiodic flowering. Little is known about the composition and evolution of the phytochrome gene family in diploid (G. herbaceum, G. raimondii) or allotetraploid (G. hirsutum, G. barbadense) cotton species. The objective of this study was to obtain a preliminary inventory and molecular-evolutionary characterization of the phytochrome gene family in cotton.     

Sampling nucleotide diversity in cotton

Background  

Cultivated cotton is an annual fiber crop derived mainly from two perennial species, Gossypium hirsutum L. or upland cotton, and G. barbadense L., extra long-staple fiber Pima or Egyptian cotton. These two cultivated species are among five allotetraploid species presumably derived monophyletically between G. arboreum and G. raimondii. Genomic-based approaches have been hindered by the limited variation within species. Yet, population-based methods are being used for genome-wide introgression of novel alleles from G. mustelinum and G. tomentosum into G. hirsutum using combinations of backcrossing, selfing, and inter-mating. Recombinant inbred line populations between genetics standards TM-1, (G. hirsutum) × 3-79 (G. barbadense) have been developed to allow high-density genetic mapping of traits.     

Genetic diversity and relationships of diploid and tetraploid cottons revealed using AFLP

Gossypium species (± 49) represent a vast resource of genetic diversity for the improvement of cultivated cotton. To determine intra- and inter-specific genetic relationships within a diverse collection of Gossypium taxa, we employed 16 AFLP primer combinations on three diploid species, Gossypium herbaceum L. (A1), Gossypium arboreum L. (A2) and Gossypium raimondii Ulbrich (D5), and 26 AD allotetraploid accessions (Gossypium barbadense L. and Gossypium hirsutum L.). A total of 1180 major AFLP bands were observed; 368 of these (31%) were polymorphic. Genetic similarities among all taxa ranged from 0.21 (between the diploid species G. arboreum and G. raimondii) up to 0.89 (within G. barbadense). Phenetic trees based on genetic similarities (UPGMA, N-J) were consistent with known taxonomic relationships. In some cases, well-supported phylogenetic relationships, as well as evidence of genetic reticulation, could also be inferred. UPGMA trees and principal coordinate analysis based on genetic similarity matrices were used to identify genetically distinct cultivars that are potentially important sources of germplasm for cotton improvement, particularly of fiber quality traits. We show that AFLP is useful for estimating genetic relationships across a wide range of taxonomic levels, and for analyzing the evolutionary and historical development of cotton cultivars at the genomic level. Received: 17 January 2000 / Accepted: 4 May 2000     

Expansion of MIR482/2118 by a class‐II transposable element in cotton

Some plant microRNA (miRNA) families contain multiple members generating identical or highly similar mature miRNA variants. Mechanisms underlying the expansion of miRNA families remain elusive, although tandem and/or segmental duplications have been proposed. In this study of two tetraploid cottons, Gossypium hirsutum and Gossypium barbadense, and their extant diploid progenitors, Gossypium arboreum and Gossypium raimondii, we investigated the gain and loss of members of the miR482/2118 superfamily, which modulates the expression of nucleotide‐binding site leucine‐rich repeat (NBS‐LRR) disease resistance genes. We found significant expansion of MIR482/2118d in G. barbadense, G. hirsutum and G. raimondii, but not in G. arboreum. Several newly expanded MIR482/2118d loci have mutated to produce different miR482/2118 variants with altered target‐gene specificity. Based on detailed analysis of sequences flanking these MIR482/2118 loci, we found that this expansion of MIR482/2118d and its derivatives resulted from an initial capture of an MIR482/2118d by a class‐II DNA transposable element (TE) in G. raimondii prior to the tetraploidization event, followed by transposition to new genomic locations in G. barbadense, G. hirsutum and G. raimondii. The ‘GosTE’ involved in the capture and proliferation of MIR482/2118d and its derivatives belongs to the PIF/Harbinger superfamily, generating a 3‐bp target site duplication upon insertion at new locations. All orthologous MIR482/2118 loci in the two diploids were retained in the two tetraploids, but mutation(s) in miR482/2118 were observed across all four species as well as in different cultivars of both G. barbadense and G. hirsutum, suggesting a dynamic co‐evolution of miR482/2118 and its NBS‐LRR targets. Our results provide fresh insights into the mechanisms contributing to MIRNA proliferation and enrich our knowledge on TEs.     

Development and mapping of SNP assays in allotetraploid cotton

A narrow germplasm base and a complex allotetraploid genome have made the discovery of single nucleotide polymorphism (SNP) markers difficult in cotton (Gossypium hirsutum). To generate sequence for SNP discovery, we conducted a genome reduction experiment (EcoRI, BafI double digest, followed by adapter ligation, biotin–streptavidin purification, and agarose gel separation) on two accessions of G. hirsutum and two accessions of G. barbadense. From the genome reduction experiment, a total of 2.04 million genomic sequence reads were assembled into contigs with an N₅₀ of 508 bp and analyzed for SNPs. A previously generated assembly of expressed sequence tags (ESTs) provided an additional source for SNP discovery. Using highly conservative parameters (minimum coverage of 8× at each SNP and 20% minor allele frequency), a total of 11,834 and 1,679 non-genic SNPs were identified between accessions of G. hirsutum and G. barbadense in genome reduction assemblies, respectively. An additional 4,327 genic SNPs were also identified between accessions of G. hirsutum in the EST assembly. KBioscience KASPar assays were designed for a portion of the intra-specific G. hirsutum SNPs. From 704 non-genic and 348 genic markers developed, a total of 367 (267 non-genic, 100 genic) mapped in a segregating F₂ population (Acala Maxxa × TX2094) using the Fluidigm EP1 system. A G. hirsutum genetic linkage map of 1,688 cM was constructed based entirely on these new SNP markers. Of the genic-based SNPs, we were able to identify within which genome (‘A’ or ‘D’) each SNP resided using diploid species sequence data. Genetic maps generated by these newly identified markers are being used to locate quantitative, economically important regions within the cotton genome.     

Spatial and temporal expression of the two sucrose synthase genes in maize: immunohistological evidence

Summary The DNAs of two diploid species of Gossypium, G. herbaceum var. africanum (A₁ genome) and G. raimondii (D₅ genome), and the allotetraploid species, G. hirsutum (A_h and D_h genomes), were characterized by kinetic analyses of single copy and repetitive sequences. Estimated haploid genome sizes of A₁ and D₅ were 1.04 pg and 0.68 pg, respectively, in approximate agreement with cytological observations that A genome chromosomes are about twice the size of D genome chromosomes. This differences in genome size was accounted for entirely by differences in the major repetitive fraction (0.56 pg versus 0.20 pg), as single copy fractions of the two genomes were essentially identical (0.41 pg for A₁ and 0.43 pg for D₅). Kinetic analyses and thermal denaturation measurements of single copy duplexes from reciprocal intergenomic hybridizations showed considerable sequence similarity between A₁ and D₅ genomes (77% duplex formation with an average thermal depression of 6 °C). Moreover, little sequence divergence was detectable between diploid single copy sequences and their corresponding genomes in the allotetraploid, consistent with previous chromosome pairing observations in interspecific F₁ hybrids.Journal paper No. 4461 of the Arizona Agricultural Experiment Station     

Rapid proliferation and nucleolar organizer targeting centromeric retrotransposons in cotton

Centromeric chromatin in most eukaryotes is composed of highly repetitive centromeric retrotransposons and satellite repeats that are highly variable even among closely related species. The evolutionary mechanisms that underlie the rapid evolution of centromeric repeats remain unknown. To obtain insight into the evolution of centromeric repeats following polyploidy, we studied a model diploid progenitor (Gossypium raimondii, D‐genome) of the allopolyploid (AD‐genome) cottons, G. hirsutum and G. barbadense. Sequence analysis of chromatin‐immunoprecipitated DNA showed that the G. raimondii centromeric repeats originated from retrotransposon‐related sequences. Comparative analysis showed that nine of the 10 analyzed centromeric repeats were absent from the centromeres in the A‐genome and related diploid species (B‐, F‐ and G‐genomes), indicating that they colonized the centromeres of D‐genome lineage after the divergence of the A‐ and D‐ ancestral species or that they were ancestrally retained prior to the origin of Gossypium. Notably, six of the nine repeats were present in both the A‐ and D‐subgenomes in tetraploid G. hirsutum, and increased in abundance in both subgenomes. This finding suggests that centromeric repeats may spread and proliferate between genomes subsequent to polyploidization. Two repeats, Gr334 and Gr359 occurred in both the centromeres and nucleolar organizer regions (NORs) in D‐ and AD‐genome species, yet localized to just the NORs in A‐, B‐, F‐, and G‐genome species. Contained within is a story of an established centromeric repeat that is eliminated and allopolyploidization provides an opportunity for reinvasion and reestablishment, which broadens our evolutionary understanding behind the cycles of centromeric repeat establishment and targeting.     

Distribution of 5S and 18S–28S rDNA loci in a tetraploid cotton (Gossypium hirsutum L.) and its putative diploid ancestors

The most widely cultivated species of cotton,Gossypium hirsutum, is a disomic tetraploid (2n=4x=52). It has been proposed previously that extant A- and D-genome species are most closely related to the diploid progenitors of the tetraploid. We used fluorescent in situ hybridization (FISH) to determine the distribution of 5S and 18S-28S rDNA loci in the A-genome speciesG. herbaceum andG. arboreum, the D-genome speciesG. raimondii andG. thurberi, and the AD tetraploidG. hirsutum. High signal-to-noise, single-label FISH was used to enumerate rDNA loci, and simultaneous, dual-label FISH was used to determine the syntenic relationships of 5S rDNA loci relative to 18S–28S rDNA loci. These techniques provided greater sensitivity than our previous methods and permitted detection of six newG. hirsutum 18S–28S rDNA loci, bringing the total number of observed loci to 11. Differences in the intensity of the hybrizization signal at these loci allowed us to designate them as major, intermediate, or minor 18–28S loci. Using genomic painting with labeled A-genome DNA, five 18S–28S loci were localized to theG. hirsutum A-subgenome and six to the D-subgenome. Four of the 11 18S–28S rDNA loci inG. hirsutum could not be accounted for in its presumed diploid progenitors, as both A-genome species has three loci and both D-genome species had four.G. hirsutum has two 5S rDNA loci, both of which are syntenic to major 18S–28S rDNA loci. All four of the diploid genomes wer examined contained a single 5S locus. InG. herbaceum (A₁) andG. thurberi (D₁), the 5S locus is syntenic to a major 18S–28S locus, but inG. arboreum (A₂) andG. raimondii (D₅), the proposed D-genome progenitor ofG. hirsutum, the 5S loci are syntenic tominor and intermediate 18S–28S loci, respecitively. The multiplicity, variation in size and site number, and lack of additivity between the tetraploid species and its putative diploid ancestors indicate that the behavior of rDNA loci in cotton is nondogmatic, and considerably more complex and dynamic than previously envisioned. The relative variability of 18S–28S rDNA loci versus 5S rDNA loci suggests that the behavior of tandem repearts can differ widely. Edited by: R. Appels     

Evaluation of Gossypium species for resistance to cotton leaf curl Burewala virus

Cotton leaf curl disease (CLCuD), caused by cotton leaf curl Burewala virus (CLCuBV), has emerged as a major threat to cotton production in Pakistan. Resistance to CLCuBV was evaluated in cultivated and wild cotton genotypes representing six Gossypium species by visual symptom scoring and virus assessment using PCR tests. Considerable variation in responses was observed when using whitefly and graft transmission to inoculate Gossypium genotypes with CLCuBV in field and greenhouse experiments. Under field evaluation, all cultivated genotypes of Gossypium hirsutum and three genotypes of G. barbadense were susceptible. Eleven genotypes that represented six wild and cultivated Gossypium species were considered to be highly resistant as they were free from infection. Similar results were obtained when these genotypes were tested using whitefly transmission. To verify these findings, 132 cultivated and wild genotypes were tested by graft inoculation. All G. hirsutum genotypes (116 cultivated, 1 wild, 1 transgenic Coker-312 and 1 non-transgenic Coker-312), three G. barbadense genotypes and one G. thurberi genotype were highly susceptible and exhibited symptoms 9–12 days after grafting. Four genotypes of G. arboreum and one genotype of G. anomalum did not express symptoms but had a detectable level of virus. One genotype of G. herbaceum and three wild genotypes of G. hirsutum showed mild symptoms (severity indexes of 1–2) and exhibited delayed disease development. These genotypes were classified as moderately resistant to resistant. Resistant genotypes that were identified in this study will be useful sources for exploitation of breeding programmes aimed at developing CLCuBV-resistant varieties and increasing genetic diversity.     

Somatic embryogenesis and plant regeneration from different wild diploid cotton (Gossypium) species

Calli were successfully induced from hypocotyls of eight wild diploid cotton species (Gossypium) on MSB (MS salts and B₅ vitamins) medium supplemented with 0.09 μM 2,4-D (2,4-dichlorophenoxyacetic acid) and 2.32 μM KT (kinetin). Plant growth regulator (PGR) combinations, adding GA₃ (Gibberellic acid), high inorganic salt stress, and PGR-free media were used to induce embryogenic calli from nonembryogenic calli. Embryogenic cultures were induced from G. aridum S. (D₄ genome), G. davidsonii K. (D₃-d genome), G. klotzschianum A. (D₃-k genome), G. raimondii U. (D₅ genome), and G. stocksii M. (E₁ genome). We then observed somatic embryogenesis in the five species while calli of G. africanum V. (A₁-2 genome), G. anomalum W. (B₁ genome), and G. bickii P. (G genome) remained nonembryogenic. Somatic embryogenesis was adjusted by changing sugar sources, regulating combinations of PGRs, and using cell suspension culture. Embryos at various developmental stages produced mature and germinating embryos when cultured on filter paper placed on the media containing different sugar sources. The utility of different sugar sources promoted globular embryos developing into cotyledonary stage and increased the frequency of cotyledonary embryos developing into normal plants. Normal plantlets were regenerated from G. davidsonii, G. klotzschianum, G. raimondii, and G. stocksii. Only abnormal plantlets were obtained in G. aridum. This work will contribute to broadening the number of regenerable cotton species and provide foundations for somatic hybridization in cotton to create new germplasm.     

Phylogenetic diversity and relationship among Gossypium germplasm using SSRs markers

Gossypium species represent a vast resource of genetic multiplicity for the improvement of cultivated cotton. To determine genetic diversity and relationships within a diverse collection of Gossypium, we employed 120 SSR primers on 20 diploid species representing seven basic genome groups of the genus Gossypium, five AD allotetraploid cotton accessions while T. populnea served as an outgroup species. Out of 120 SSR primers, 49 pairs are polymorphic, which produced a total of 99 distinct alleles with an average of 2.0 alleles per primer pair. A total of 1139 major SSR bands were observed. Genetic similarities among all the diploid species ranged from 0.582 (between G. herbaceum and G. trilobum) up to 0.969 (between G. arboreum and G. herbaceum). Phylogenetic trees based on genetic similarities were consistent with known taxonomic relationships. The results also indicated that G. raimondii is the closest living relative of the ancestral D-genome donor of tetraploid species and the A-genome donor is much similar to the present-day G. herbaceum and G. arboreum. Ancient tetraploid cotton species were formed by hybridizing and chromosome doubling between them, then different tetraploid cotton species appeared by further geographical and genetic isolation and separating differentiation. The results showed that SSRs could be an ideal means for the identification of the genetic diversity and relationship of cotton resources at the genomic level.     

Genetic diversity and geographic pattern in early South American cotton domestication

Amplified fragment length polymorphism fingerprinting was applied to survey the genetic diversity of primitive South American Gossypium barbadense cotton for establishing a possible link to its pre-Columbian expansion. New germplasm was collected along coastal Peru and over an Andean transect in areas where most of the archaeological evidence relating to cotton domestication has been recorded. Gene bank material of three diploid (G. raimondii, G. arboreum, and G. herbaceum) and four allotetraploid cotton species (G. hirsutum, G. mustelinum, G. tomentosum and additional G. barbadense) was added for inter- and intra-specific comparison. Eight primer combinations yielded 340 polymorphic bands among the 131 accessions. The obtained neighbor joining and unweighted pair-group method with arithmetic means are in full agreement with the known cytogenetics of the tetraploid cottons and their diploid genome donors. The four tetraploid species are clearly distinct based on taxonomic classification. The genetic diversity within G. barbadense reveals geographic patterns. The locally maintained cottons from coastal Peru display a distinct genetic diversity that mirrors their primitive agro-morphological traits. Accessions from the northernmost coast of Peru and from southwestern (SW) Ecuador cluster basal to the east-of-Andes accessions. The remaining accessions from Bolivia, Brazil, Columbia, Venezuela, and the Caribbean and Pacific islands cluster with the east-of-Andes accessions. Northwestern Peru/SW Ecuador (the area flanking the Guayaquil gulf) appears to be the center of the primitive domesticated G. barbadense cotton from where it spread over the Andes and expanded into its pre-Columbian range.This publication is dedicated to Prof. Dr. Drs.h.c. Gerhard Röbbelen on the occasion of his 75th birthday     

Breeding for ”low-gossypol seed and high-gossypol plants” in upland cotton. Analysis of tri-species hybrids and backcross progenies using AFLPs and mapped RFLPs

This work aims at breeding upland cotton [Gossypium hirsutum L., 2(AD)₁ genome] with a reduced level of gossypol in the seeds for optimal food and feed uses, and a high gossypol level in the remaining organs for resistance to pests. Two tri-species Gossypium hybrids, (G. thurberi–G. sturtianum–G. hirsutum and G. hirsutum–G. raimondii–G. sturtianum) including G. sturtianum (2C₁) as a donor, G. thurberi (2D₁) and G. raimondii (2D₅) as a bridge species, were created. Recurrent selection initiated with these tri-species hybrids produced backcross (BC) progenies expressing the ”low-gossypol seed and high-gossypol plant” trait at different levels. We used AFLP markers to assess the genetic similarity among the germplasm and RFLP probes to tag the introgression of specific chomosome segments from the parental species. Five pairs of AFLP primers generated 477 fragments, among which 417 (87.4%) were polymorphic. The genetic similarity between the upland cotton and the wild species ranged from 29.5 to 43.2%, while similarity reached 80% between upland cotton and BC3 plants. Introgression of species-specific AFLPs was evident from all the parental species and confirmed the hybrid origin of the analyzed progenies. Southern-blot analysis based on 49 RFLP probes allowed us to trace the introgression of parental DNA segments in the tri-species hybrids and in three generations of backcross. Introgression was evident from 11, 8 and 7 linkage groups of G. sturtianum, G. raimondii and G. thurberi respectively. The types of introgression revealed by RFLP probes are discussed, and breeding schemes to enhance recombination are proposed. The ability to trace DNA segments of known chromosomal locations from the donor G. sturtianum through segregating generations is a starting point to map the ”low-gossypol seed and high-gossypol plant” traits. Received: 5 January 1999 / Accepted: 17 June 1999     

Comparative Genome Analyses Highlight Transposon-Mediated Genome Expansion and the Evolutionary Architecture of 3D Genomic Folding in Cotton

Transposable element (TE) amplification has been recognized as a driving force mediating genome size expansion and evolution, but the consequences for shaping 3D genomic architecture remains largely unknown in plants. Here, we report reference-grade genome assemblies for three species of cotton ranging 3-fold in genome size, namely Gossypium rotundifolium (K₂), G. arboreum (A₂), and G. raimondii (D₅), using Oxford Nanopore Technologies. Comparative genome analyses document the details of lineage-specific TE amplification contributing to the large genome size differences (K₂, 2.44 Gb; A₂, 1.62 Gb; D₅, 750.19 Mb) and indicate relatively conserved gene content and synteny relationships among genomes. We found that approximately 17% of syntenic genes exhibit chromatin status change between active (“A”) and inactive (“B”) compartments, and TE amplification was associated with the increase of the proportion of A compartment in gene regions (∼7,000 genes) in K₂ and A₂ relative to D₅. Only 42% of topologically associating domain (TAD) boundaries were conserved among the three genomes. Our data implicate recent amplification of TEs following the formation of lineage-specific TAD boundaries. This study sheds light on the role of transposon-mediated genome expansion in the evolution of higher-order chromatin structure in plants.     

Ancient origin of the vacuolar H+-ATPase 69-kilodalton catalytic subunit superfamily

Recently, two distinct cDNA clones encoding the catalytic subunit of the vacuolar H⁺-ATPase (V-ATPase) were isolated from the allotetraploid cotton species Gossypium hirsutum L. cv Acala SJ-2 (Wilkins 1992, 1993). Differences in the nucleotide sequence of these clones were used as molecular markers to explore the organization and structure of the V-ATPase catalytic subunit genes in the A and D genomes of diploid and allotetraploid cotton species. Nucleotide sequencing of polymerase chain reaction (PCR) products amplified from G. arboreum (A₂, 2n=26), G. raimondii (D₅, 2n=26), and G. hirsutum cv Acala SJ-2 [(AD)₁, 2n=4x=52] revealed a V-ATPase catalytic subunit organization more complex than indicated hitherto in any species, including higher plants. In the genus Gossypium, the V-ATPase catalytic subunit genes are organized as a superfamily comprising two diverse but closely related multigene families, designated as vat69A and vat69B, present in both diploid and allotetraploid species. As expected, each vat69 subfamily is correspondingly more complex in the allotetraploid species due to the presence of both A and D alloalleles. Because of this, about one-half of the complex organization of V-ATPase catalytic subunit genes predates polyploidization and speciation of New World tetraploid species. Comparison of plant and fungal V-ATPase catalytic subunit gene structure indicates that introns accrued in the plant homologs following the bifurcation of plant and fungi but prior to the gene duplication event that gave rise to the vat69A and vat69B genes approximately 45 million years ago. The structural complexity of plant V-ATPase catalytic subunit genes is highly conserved, indicating the presence of at least ten introns dispersed throughout the coding region.     

GbEXPATR,a species‐specific expansin,enhances cotton fibre elongation through cell wall restructuring

Cotton provides us the most important natural fibre. High fibre quality is the major goal of cotton breeding, and introducing genes conferring longer, finer and stronger fibre from Gossypium barbadense to Gossypium hirsutum is an important breeding strategy. We previously analysed the G. barbadense fibre development mechanism by gene expression profiling and found two homoeologous fibre‐specific α‐expansins from G. barbadense, GbEXPA2 and GbEXPATR. GbEXPA2 (from the D_T genome) is a classical α‐expansin, while its homoeolog, GbEXPATR (A_T genome), encodes a truncated protein lacking the normal C‐terminal polysaccharide‐binding domain of other α‐expansins and is specifically expressed in G. barbadense. Silencing EXPA in G. hirsutum induced shorter fibres with thicker cell walls. GbEXPA2 overexpression in G. hirsutum had no effect on mature fibre length, but produced fibres with a slightly thicker wall and increased crystalline cellulose content. Interestingly, GbEXPATR overexpression resulted in longer, finer and stronger fibres coupled with significantly thinner cell walls. The longer and thinner fibre was associated with lower expression of a number of secondary wall‐associated genes, especially chitinase‐like genes, and walls with lower cellulose levels but higher noncellulosic polysaccharides which advocated that a delay in the transition to secondary wall synthesis might be responsible for better fibre. In conclusion, we propose that α‐expansins play a critical role in fibre development by loosening the cell wall; furthermore, a truncated form, GbEXPATR, has a more dramatic effect through reorganizing secondary wall synthesis and metabolism and should be a candidate gene for developing G. hirsutum cultivars with superior fibre quality.     

Pros and cons of using genomic SSRs and EST-SSRs for resolving phylogeny of the genus Gossypium

The genus Gossypium is comprised of 50 diverse cotton species representing eight different genomes (A through G and K), however, phylogenetic relationship using various DNA marker types such as RAPD and SSRs was determined on limited number of cotton species. In this report, we have demonstrated the application of genomic SSRs (gSSRs) and EST-SSRs, and after combining both the data sets, for resolving the phylogenies of 36 cotton species including seven races. Out of the 100 primer pairs surveyed (50 for gSSRs and 50 for EST-SSRs), 75 produced scorable amplification products in all species. Out of these, 73 were found to be polymorphic and amplified 135 alleles ranging from 1 to 5 alleles per SSR marker (average 2.87 alleles per marker). The gSSRs amplified higher number of alleles (72) compared to the EST-SSRs (63). In total 22 highly informative SSRs with PIC values ≥0.5 were identified. Genomic SSRs containing di-while EST-SSRs containing tri-nucleotide repeats exhibited high polymorphism compared to the other nucleotide repeats containing gSSRs/EST-SSRs. Number of tandem repeats and polymorphism were positively correlated. Neither the type of chromosome nor the location of the SSRs showed association with the polymorphism. Gossypium herbaceum var. africanum (Watt) Hutch. ex and Ghose and Gossypium robinsonii F. Muell. were found the most genetically diverse, while among races of Gossypium hirsutum L. “yucatanense” and G. hirsutum “punctatum” were found genetically diverse. Of the three data sets, clustering analysis based on EST-SSRs and combined data sets, revealed parallel results reported in earlier studies. This study further confirmed that Gossypium darwinii Watt has close relationship with Gossypium barbadense L. Moreover, Gossypium raimondii Ulbr. and G. herbaceum/Gossypium arboreum L. are close living relatives of the ancestor allotetraploid species. Our studies suggest that for resolving phylogenetic relationship among the various plant species EST-SSRs could be a better choice. This information can be instrumental in transferring novel alleles or loci from the wild species into the cultivated cotton species which would set a stage for cultivating genetically diverse cultivars—a way to achieve sustainable cotton production in changing climate.