Allotetraploid cotton is an economically important natural-fiber-producing crop worldwide. After polyploidization, Gossypium hirsutum L. evolved to produce a higher fiber yield and to better survive harsh environments than Gossypium barbadense, which produces superior-quality fibers. The global genetic and molecular bases for these interspecies divergences were unknown. Here we report high-quality de novo-assembled genomes for these two cultivated allotetraploid species with pronounced improvement in repetitive-DNA-enriched centromeric regions. Whole-genome comparative analyses revealed that speciesspecific alterations in gene expression, structural variations and expanded gene families were responsible for speciation and the evolutionary history of these species. These findings help to elucidate the evolution of cotton genomes and their domestication history. The information generated not only should enable breeders to improve fiber quality and resilience to ever-changing environmental conditions but also can be translated to other crops for better understanding of their domestication history and use in improvement.

To increase the numbers of microsatellites available for use in constructing a genetic map, and facilitate the use of functional genomics to elucidate fiber development and breeding in cotton, we sampled microsatellite sequences from expressed sequence tags (ESTs) transcribed during fiber elongation in the A-genome species Gossypium arboreum to evaluate their frequency of occurrence, level of polymorphism and distribution in the At and Dt subgenomes of tetraploid cotton. From among ESTs derived from G. arboreum fibers at 7-10 days post anthesis (dpa), 931 ESTs were found to contain simple sequence repeats (SSRs); 544 (58.4%) EST-SSR primer pairs were developed, and 468 (86%) amplified PCR products from allotetraploid cotton (G. hirsutumcv. TM-1 and G. barbadense cv. Hai7124). However, only 99 (18.2%) of these were found to be polymorphic and segregating in our interspecific BC1 mapping population [(TM-1xHai7124)xTM-1]. In these amplified and informative EST-SSRs, hexa- and tri-nucleotide repeat motifs were the most frequent, representing 40.1 and 30%, respectively, of the total. A total of 111 loci detected with these 99 EST-SSRs were integrated into our backbone map including 511 SSR loci. The distribution of the EST-SSRs appeared to be non-random, since 72 loci were anchored to the At and 37 to the Dt subgenome of allotetraploid cotton based on linkage tests. Interestingly, out of the 10 pairs of duplicate loci amplified, seven were mapped to the corresponding homologous linkage groups and/or chromosomes. BLASTX analysis revealed that 69 of the 99 ESTs showed significant similarities to known genes. Some genes important for fiber development, such as sucrose synthase, were mapped to corresponding chromosomes. These EST-SSRs provide structural and functional genomic information that will be useful for understanding cotton fiber development.

We characterize GoSP genes underlying the development of cotton plants with short branches and clustered bolls, a phenotype that allows higher planting density and promotes increased fiber yield per acre.

Interspecific genomic variation can provide a genetic basis for local adaptation and domestication. A series of studies have presented its role of interspecific haplotypes and introgressions in adaptive traits, but few studies have addressed their role in improving agronomic character. Two allotetraploid Gossypium species, Gossypium barbadense (Gb) and G. hirsutum (Gh) originating from the Americas, are cultivated independently. Here, through sequencing and the comparison of one GWAS panel in 229 Gb accessions and two GWAS panels in 491 Gh accessions, we found that most associated loci or functional haplotypes for agronomic traits were highly divergent, representing the strong divergent improvement between Gb and Gh. Using a comprehensive interspecific haplotype map, we revealed that six interspecific introgressions from Gh to Gb were significantly associated with the phenotypic performance of Gb, which could explain 5%-40% of phenotypic variation in yield and fibre qualities. In addition, three introgressions overlapped with six associated loci in Gb, indicating that these introgression regions were under further selection and stabilized during improvement. A single interspecific introgression often possessed yield-increasing potential but decreased fibre qualities, or the opposite, making it difficult to simultaneously improve yield and fibre qualities. Our study not only has proved the importance of interspecific functional haplotypes or introgressions in the divergent improvement of Gb and Gh, but also supports their potential value in further humanmediated hybridization or precision breeding.