Multistage temporal transcriptomic atlas unveils major contributor to reproductive phase in upland cotton.

Flowering is a major developmental transition in plants, but asynchronous flowering hinders the utilization of wild cotton relatives in breeding programs. We performed comparative transcriptomic profiling of early- and late-flowering Gossypium hirsutum genotypes to elucidate genetic factors influencing reproductive timing. Shoot apices were sampled from the photoperiod-sensitive landrace G. hirsutum purpurascens (GhP) and early-maturing variety ZhongMianSuo (ZMS) at five time points following the emergence of sympodial nodes. RNA-sequencing revealed extensive transcriptional differences during floral transition. Numerous flowering-associated genes exhibited genotype-specific expression, including FLOWERING LOCUS T (FT) homologs upregulated in ZMS. FT-interacting factors like SOC1 and CO-like also showed higher expression in ZMS, implicating florigen pathways in early flowering. Additionally, circadian clock and light signalling components were misregulated between varieties, suggesting altered photoperiod responses in GhP. Weighted co-expression network analysis specifically linked a module enriched for circadian-related genes to GhP's late flowering. Through an integrated transcriptome analysis, we defined a regulatory landscape of reproductive phase change in cotton. Differentially expressed genes related to photoperiod, circadian clock, and light signalling likely contribute to delayed flowering in wild cottons. Characterization of upstream flowering regulators will enable modifying photoperiod sensitivity and expand germplasm use for cotton improvement. This study provides candidate targets for elucidating interactive mechanisms that control cotton flowering time across diverse genotypes.

Exploring the regulatory role of non-coding RNAs in fiber development and direct regulation of GhKCR2 in the fatty acid metabolic pathway in upland cotton.

Cotton fiber holds immense importance as the primary raw material for the textile industry. Consequently, comprehending the regulatory mechanisms governing fiber development is pivotal for enhancing fiber quality. Our study aimed to construct a regulatory network of competing endogenous RNAs (ceRNAs) and assess the impact of non-coding RNAs on gene expression throughout fiber development. Through whole transcriptome data analysis, we identified differentially expressed genes (DEGs) regulated by non-coding RNA (ncRNA) that were predominantly enriched in phenylpropanoid biosynthesis and the fatty acid elongation pathway. This analysis involved two contrasting phenotypic materials (J02-508 and ZRI015) at five stages of fiber development. Additionally, we conducted a detailed analysis of genes involved in fatty acid elongation, including KCS, KCR, HACD, ECR, and ACOT, to unveil the factors contributing to the variation in fatty acid elongation between J02-508 and ZRI015. Through the integration of histochemical GUS staining, dual luciferase assay experiments, and correlation analysis of expression levels during fiber development stages for lncRNA MSTRG.44818.23 (MST23) and GhKCR2, we elucidated that MST23 positively regulates GhKCR2 expression in the fatty acid elongation pathway. This identification provides valuable insights into the molecular mechanisms underlying fiber development, emphasizing the intricate interplay between non-coding RNAs and protein-coding genes.

Identification and characterization of candidate genes for primary root length in Asiatic cotton (Gossypium arboreum L.).

KEY MESSAGE: One major gene controlling primary root length (PRL) in Gossypium arboreum is identified and this research provides a theoretical basis for root development for cotton. Primary root elongation is an essential process in plant root system structure. Here, we investigated the primary root length (PRL) of 215 diploid cotton (G. arboreum) accessions at 5, 8, 10, 15 days after sowing. A Genome-wide association study was performed for the PRL, resulting in 49 significant SNPs associated with 32 putative candidate genes. The SNP with the strongest signal (Chr07_8047530) could clearly distinguish the PRLs between accessions with two haplotypes. GamurG is the only gene that showed higher relative expression in the long PRL genotypes than the short PRL genotypes, which indicated it was the most likely candidate gene for regulating PRL. Moreover, the GamurG-silenced cotton seedlings showed a shorter PRL, while the GamurG-overexpressed Arabidopsis exhibited a significantly longer PRL. Our findings provide insight into the regulation mechanism of cotton root growth and will facilitate future breeding programs to optimize the root system structure in cotton.

Transcriptomic analysis reveals the beneficial effects of salt priming on enhancing defense responses in upland cotton under successive salt stress.

Priming-mediated stress tolerance in plants stimulates defense mechanisms and enables plants to cope with future stresses. Seed priming has been proven effective for tolerance against abiotic stresses; however, underlying genetic mechanisms are still unknown. We aimed to assess upland cotton genotypes and their transcriptional behaviors under salt priming and successive induced salt stress. We pre-selected 16 genotypes based on previous studies and performed morpho-physiological characterization, from which we selected three genotypes, representing different tolerance levels, for transcriptomic analysis. We subjected these genotypes to four different treatments: salt priming (P0), salt priming with salinity dose at 3-true-leaf stage (PD), salinity dose at 3-true-leaf stage without salt priming (0D), and control (CK). Although the three genotypes displayed distinct expression patterns, we identified common differentially expressed genes (DEGs) under PD enriched in pathways related to transferase activity, terpene synthase activity, lipid biosynthesis, and regulation of acquired resistance, indicating the beneficial role of salt priming in enhancing salt stress resistance. Moreover, the number of unique DEGs associated with G. hirsutum purpurascens was significantly higher compared to other genotypes. Coexpression network analysis identified 16 hub genes involved in cell wall biogenesis, glucan metabolic processes, and ribosomal RNA binding. Functional characterization of XTH6 (XYLOGLUCAN ENDOTRANSGLUCOSYLASE/HYDROLASE) using virus-induced gene silencing revealed that suppressing its expression improves plant growth under salt stress. Overall, findings provide insights into the regulation of candidate genes in response to salt stress and the beneficial effects of salt priming on enhancing defense responses in upland cotton.

Genome-wide analysis of the serine carboxypeptidase-like protein family reveals Ga09G1039 is involved in fiber elongation in cotton.

The Gossypium is a model genus for understanding polyploidy and the evolutionary pattern of inheritance. This study aimed to investigate the characteristics of SCPLs in different cotton species and their role in fiber development. A total of 891 genes from one typical monocot and ten dicot species were naturally divided into three classes based on phylogenetic analysis. The SCPL gene family in cotton has undergone intense purifying selection with some functional variation. Segmental duplication and whole genome duplication were shown to be the two main reasons for the increase in the number of genes during cotton evolution. The identification of Gh_SCPL genes exhibiting differential expression in particular tissues or response to environmental stimuli provides a new measure for the in-depth characterization of selected genes of importance. Ga09G1039 was involved in the developmental process of fibers and ovules, and it is significantly different from proteins from other cotton species in terms of phylogenetic, gene structure, conserved protein motifs and tertiary structure. Overexpression of Ga09G1039 significantly increased the length of stem trichomes. Ga09G1039 may be a serine carboxypeptidase protein with hydrolase activity, according to functional region, prokaryotic expression, and western blotting analysis. The results provide a comprehensive overview of the genetic basis of SCPLs in Gossypium and further our knowledge in understanding the key aspects of SCPLs in cotton with their potential role in fiber development and stress resistance.

Silencing of GhORP_A02 enhances drought tolerance in Gossypium hirsutum.

BACKGROUND: ORP (Oxysterol-binding protein-related proteins) genes play a role in lipid metabolism, vesicular transferring and signaling, and non-vesicular sterol transport. However, no systematic identification and analysis of ORP genes have been reported in cotton. RESULT: In this study, we identified 14, 14, 7, and 7 ORP genes in G. hirsutum, G. barbadense, G. arboreum, and G. raimondii, respectively. Phylogenetic analysis showed that all ORP genes could be classified into four groups. Gene structure and conserved motif analysis suggest that the function of this gene family was conserved. The Ka/Ks analysis showed that this gene family was exposed to purifying selection during evolution. Transcriptome data showed that four ORP genes, especially GhORP_A02, were induced by abiotic stress treatment. The cis-acting elements in the ORP promoters were responsive to phytohormones and various abiotic stresses. The silenced plants of GhORP_A02 were more sensitive to drought stress when compared to control. CONCLUSION: The major finding of this study shed light on the potential role of ORP genes in abiotic stress and provided a fundamental resource for further analysis in cotton.

Transcriptomic and metabolomic profiling of flavonoid biosynthesis provides novel insights into petals coloration in Asian cotton (Gossypium arboreum L.).

BACKGROUND: Asian cotton (Gossypium arboreum L.), as a precious germplasm resource of cotton with insect resistance and stress tolerance, possesses a broad spectrum of phenotypic variation related to pigmentation. Flower color affects insect pollination and the ornamental value of plants. Studying flower color of Asian cotton varieties improves the rate of hybridization and thus enriches the diversity of germplasm resources. Meanwhile, it also impacts the development of the horticultural industry. Unfortunately, there is a clear lack of studies concerning intricate mechanisms of cotton flower-color differentiation. Hereby, we report an integrative approach utilizing transcriptome and metabolome concerning flower color variation in three Gossypium arboreum cultivars. RESULTS: A total of 215 differentially accumulated metabolites (DAMs) were identified, including 83 differentially accumulated flavonoids (DAFs). Colorless kaempferol was more abundant in white flowers, while gossypetin-fer showed specificity in white flowers. Quercetin and gossypetin were the main contributors to yellow petal formation. Pelargonidin 3-O-beta-D-glucoside and cyanidin-3-O-(6''-Malonylglucoside) showed high accumulation levels in purple petals. Quercetin and gossypetin pigments also promoted purple flower coloration. Moreover, 8178 differentially expressed genes (DEGs) were identified by RNA sequencing. The correlation results between total anthocyanins and DEGs were explored, indicating that 10 key structural genes and 29 transcription factors promoted anthocyanin biosynthesis and could be candidates for anthocyanin accumulation. Ultimately, we constructed co-expression networks of key DAFs and DEGs and demonstrated the interactions between specific metabolites and transcripts in different color flowers. CONCLUSION: This study provides new insights into elucidating the regulatory mechanisms of cotton flower color and lays a potential foundation for generate cotton varieties with highly attractive flowers for pollinators.

Differential seedling growth and tolerance indices reflect drought tolerance in cotton.

BACKGROUND: Cotton production is adversely effected by drought stress. It is exposed to drought stress at various critical growth stages grown under a water scarcity environment. Roots are the sensors of plants; they detect osmotic stress under drought stress and play an important role in plant drought tolerance mechanisms. The seedling stage is very sensitive to drought stress, and it needed to explore the methods and plant characteristics that contribute to drought tolerance in cotton. RESULTS: Initially, seedlings of 18 genotypes from three Gossypium species: G. hirsutum, G. barbadense, and G. arboreum, were evaluated for various seedling traits under control (NS) and drought stress (DS). Afterward, six genotypes, including two of each species, one tolerant and one susceptible, were identified based on the cumulative drought sensitivity response index (CDSRI). Finally, growth rates (GR) were examined for shoot and root growth parameters under control and DS in experimental hydroponic conditions. A significant variation of drought stress responses was observed across tested genotypes and species. CDSRI allowed here to identify the drought-sensitive and drought-resistant cultivar of each investigated species. Association among root and shoots growth traits disclosed influential effects of enduring the growth under DS. The traits including root length, volume, and root number were the best indicators with significantly higher differential responses in the tolerant genotypes. These root growth traits, coupled with the accumulation of photosynthates and proline, were also the key indicators of the resistance to drought stress. CONCLUSION: Tolerant genotypes have advanced growth rates and the capacity to cop with drought stress by encouraging characteristics, including root differential growth traits coupled with physiological traits such as chlorophyll and proline contents. Tolerant and elite genotypes of G. hirsutum were more tolerant of drought stress than obsolete genotypes of G. barbadense and G. arboreum. Identified genotypes have a strong genetic basis of drought tolerance, which can be used in cotton breeding programs.

Integrating Genome-wide association and whole transcriptome analysis to reveal genetic control of leaf traits in Gossypium arboreum L.

Leaves are important organs for crop photosynthesis and transpiration, and their morphological characteristics can directly reflect the growth state of plants. Accurate measurement of leaf traits and mining molecular markers are of great significance to the study of cotton growth. Here, we performed a Genome-wide association study on 7 leaf traits in 213 Asian cotton accessions. 32 significant SNPs and 44 genes were identified. A field experiment showed significant difference in leaf hair and leaf area between DPL971 and its natural mutant DPL972. We also compared the leaf transcriptome difference between DPL971 and DPL972, and found a batch of differentially expressed genes and non-coding RNAs (including lncRNAs, microRNAs, and circRNAs). After integrating the GWAS and transcriptome results, we finally selected two coding genes (Ga03G2383 and Ga05G3412) and two microRNAs (hbr-miR156, unconservative_Chr03_contig343_2364) as the candidate for leaf traits. Those findings will provide important genomic resources for cotton leaf improvement breeding.

Deltapine 15 contributes to the genomic architecture of modern upland cotton cultivars.

Foundation parents play a critical role in the genetic constituents of the derived genotypes. Deltapine-15 (DLP-15), introduced to China in 1950, is one of the most commonly used parents for early breeding programs in China. However, the formation and inheritance patterns of genomic constituents have not been studied. Therefore, this study aimed at understanding and exploring the genomic architecture of 146 DLP-15 derived cultivars with a common foundation parent DLP-15. Population structure based on sequencing data clustered genotypes into two groups (G1 and G2) supported by principal component analysis. Further exploration led to the identification of Chr-A08 with significantly differentiated regions between two groups. Moreover, we identified genome-wide identity by descent (IBD) segments (840 segments) to understand the genomic inheritance pattern in DLP-15 derived cultivars, spanning the 20-95 Mb region on Chr-A08. Interestingly, Chr-A08 depicted a unique inheritance pattern from DLP-15 to its derived cultivars. IBD-segment-based haplotype analysis suggested significant differences among the two groups. Phenotypic trait association with DLP-derived haplotypes concerning Chr-A08 suggested a significant increase in yield and fiber quality. Furthermore, distinguished IBD segments overlapped with previously reported QTLs concerning fiber yield and quality. Our results systematically identified genomic signatures transmitted from the foundation parent DLP-15 to its derived cultivars and provided a basis for further exploiting excellent haplotypes associated with DLP-15.

Genome-wide association study for seedling biomass-related traits in Gossypium arboreum L.

BACKGROUND: Seedling stage plant biomass is usually used as an auxiliary trait to study plant growth and development or stress adversities. However, few molecular markers and candidate genes of seedling biomass-related traits were found in cotton. RESULT: Here, we collected 215 Gossypium arboreum accessions, and investigated 11 seedling biomass-related traits including the fresh weight, dry weight, water content, and root shoot ratio. A genome-wide association study (GWAS) utilizing 142,5003 high-quality SNPs identified 83 significant associations and 69 putative candidate genes. Furthermore, the transcriptome profile of the candidate genes emphasized higher expression of Ga03G1298, Ga09G2054, Ga10G1342, Ga11G0096, and Ga11G2490 in four representative cotton accessions. The relative expression levels of those five genes were further verified by qRT-PCR. CONCLUSIONS: The significant SNPs, candidate genes identified in this study are expected to lay a foundation for studying the molecular mechanism for early biomass development and related traits in Asian cotton.

Genetic Factors Underlying Single Fiber Quality in A-Genome Donor Asian Cotton (Gossypium arboreum).

The study of A-genome Asian cotton as a potential fiber donor in Gossypium species may offer an enhanced understanding of complex genetics and novel players related to fiber quality traits. Assessment of individual fibers providing classified fiber quality information to the textile industry is Advanced Fiber Information System (AFIS) in the recent technological era. Keeping the scenario, a diverse collection of 215 Asiatic cotton accessions were evaluated across three agro-ecological zones of China. Genome-Wide Association Studies (GWAS) was performed to detect association signals related to 17 AFIS fiber quality traits grouped into four categories viz: NEPs, fiber length, maturity, and fineness. Significant correlations were found within as well as among different categories of various traits related to fiber quality. Fiber fineness has shown a strong correlation to all other categories, whereas these categories are shown interrelationships via fiber-fineness. A total of 7,429 SNPs were found in association with 17 investigated traits, of which 177 were selected as lead SNPs. In the vicinity of these lead SNPs, 56 differentially expressed genes in various tissues/development stages were identified as candidate genes. This compendium connecting trait-SNP-genes may allow further prioritization of genes in GWAS loci to enable mechanistic studies. These identified quantitative trait nucleotides (QTNs) may prove helpful in fiber quality improvement in Asian cotton through marker-assisted breeding as well as in reviving eroded genetic factors of G. hirsutum via introgression breeding.

GhGASA10-1 promotes the cell elongation in fiber development through the phytohormones IAA-induced.

BACKGROUND: Cotton is an important cash crop. The fiber length has always been a hot spot, but multi-factor control of fiber quality makes it complex to understand its genetic basis. Previous reports suggested that OsGASR9 promotes germination, width, and thickness by GAs in rice, while the overexpression of AtGASA10 leads to reduced silique length, which is likely to reduce cell wall expansion. Therefore, this study aimed to explore the function of GhGASA10 in cotton fibers development. RESULTS: To explore the molecular mechanisms underlying fiber elongation regulation concerning GhGASA10-1, we revealed an evolutionary basis, gene structure, and expression. Our results emphasized the conservative nature of GASA family with its origin in lower fern plants S. moellendorffii. GhGASA10-1 was localized in the cell membrane, which may synthesize and transport secreted proteins to the cell wall. Besides, GhGASA10-1 promoted seedling germination and root extension in transgenic Arabidopsis, indicating that GhGASA10-1 promotes cell elongation. Interestingly, GhGASA10-1 was upregulated by IAA at fiber elongation stages. CONCLUSION: We propose that GhGASA10-1 may promote fiber elongation by regulating the synthesis of cellulose induced by IAA, to lay the foundation for future research on the regulation networks of GASA10-1 in cotton fiber development.

Genomic insight into the divergence and adaptive potential of a forgotten landrace G. hirsutum L. purpurascens.

Wild progenitors are an excellent source for strengthening the genetic basis and accumulation of desirable variation lost because of directional selection and adaptation in modern cultivars. Here, we re-evaluate a landrace of Gossypium hirsutum, formerly known as Gossypium purpurascens. Our study seeks to understand the genomic structure, variation, and breeding potential of this landrace, providing potential insights into the biogeographic history and genomic changes likely associated with domestication. A core set of accessions, including current varieties, obsolete accessions, G. purpurascens, and other geographical landraces, are subjected to genotyping along with multilocation phenotyping. Population fixation statistics suggests a marked differentiation between G. purpurascens and three other groups, emphasizing the divergent genomic behavior of G. purpurascens. Phylogenetic analysis establishes the primitive nature of G. purpurascens, identifying it as a vital source of functional variation, the inclusion of which in the upland cotton (cultivated G. hirsutum) gene pool may broaden the genetic basis of modern cultivars. Genome-wide association results indicate multiple loci associated with domestication regions corresponding to flowering and fiber quality. Moreover, the conserved nature of G. purpurascens can also provide insights into the evolutionary process of G. hirsutum.

Dissecting Source-Sink Relationship of Subtending Leaf for Yield and Fiber Quality Attributes in Upland Cotton (Gossypium hirsutum L.).

Photosynthesis as a source is a significant contributor to the reproductive sink affecting cotton yield and fiber quality. Moreover, carbon assimilation from subtending leaves adds up a significant proportion to the reproductive sink. Therefore, this study aimed to address the source-sink relationship of boll subtending leaf with fiber quality and yield related traits in upland cotton. A core collection of 355 upland cotton accessions was subjected to subtending leaf removal treatment effects across 2 years. The analysis of variance suggested a significant effect range in the source-sink relationship under subtending leaf removal effects at different growth stages. Further insight into the variation was provided by the correlation analysis and principal component analysis. A significant positive correlation between different traits was observed and the multivariate analysis including hierarchical clustering and principal component analysis (PCA) categorised germplasm accessions into three groups on the basis of four subtending leaf removal treatment effects across 2 years. A set of genotypes with the lowest and highest treatment effects has been identified. Selected accessions and the outcome of the current study may provide a basis for a further study to explore the molecular mechanism of source-sink relationship of boll subtending leaf and utilization of breeding programs focused on cotton improvement.

GWAS Mediated Elucidation of Heterosis for Metric Traits in Cotton (Gossypium hirsutum L.) Across Multiple Environments.

For about a century, plant breeding has widely exploited the heterosis phenomenon-often considered as hybrid vigor-to increase agricultural productivity. The ensuing F1 hybrids can substantially outperform their progenitors due to heterozygous combinations that mitigate deleterious mutations occurring in each genome. However, only fragmented knowledge is available concerning the underlying genes and processes that foster heterosis. Although cotton is among the highly valued crops, its improvement programs that involve the exploitation of heterosis are still limited in terms of significant accomplishments to make it broadly applicable in different agro-ecological zones. Here, F1 hybrids were derived from mating a diverse Upland Cotton germplasm with commercially valuable cultivars in the Line × Tester fashion and evaluated across multiple environments for 10 measurable traits. These traits were dissected into five different heterosis types and specific combining ability (SCA). Subsequent genome-wide predictions along-with association analyses uncovered a set of 298 highly significant key single nucleotide polymorphisms (SNPs)/Quantitative Trait Nucleotides (QTNs) and 271 heterotic Quantitative Trait Nucleotides (hQTNs) related to agronomic and fiber quality traits. The integration of a genome wide association study with RNA-sequence analysis yielded 275 candidate genes in the vicinity of key SNPs/QTNs. Fiber micronaire (MIC) and lint percentage (LP) had the maximum number of associated genes, i.e., each with 45 related to QTNs/hQTNs. A total of 54 putative candidate genes were identified in association with HETEROSIS of quoted traits. The novel players in the heterosis mechanism highlighted in this study may prove to be scientifically and biologically important for cotton biologists, and for those breeders engaged in cotton fiber and yield improvement programs.

Identification and functional analysis of 9-cis-epoxy carotenoid dioxygenase (NCED) homologs in G. hirsutum.

9-cis-epoxy carotenoid dioxygenase (NCED) is a fundamental enzyme, which plays an essential role in the process of organ development and stress resistance by regulating abscisic acid (ABA) synthesis in plant. In this study, a total of 7, 7, 14 and 14 NCED genes were identified from the genomes of G. arboreum, G. raimondii, G. barbadense and G. hirsutum, respectively. Phylogenetic tree showed that all forty-two NCED genes could be classified into three groups in cotton genus. Collinear analysis revealed that the NCED genes in G. hirsutum were not amplified by tandem repeats after polyploidy events. The function of NCED genes was evaluated between two accessions with contrasting plant height. The results showed that expression of the NCED genes in dwarf accession was higher than that in taller ones. GhNCED1-silenced cotton plants confirmed that suppression of NCED genes could increase the plant height, but reduce the resistance abilities to drought and salt stress. Our study systematically identified the homologs of NCED genes and their functions in cotton, which could provide new genetic resources for improving plant height and stress in future cotton breeding.

Introgression Leads to Genomic Divergence and Responsible for Important Traits in Upland Cotton.

Understanding the genetic diversity and population structure of germplasms is essential when selecting parents for crop breeding. The genomic changes that occurred during the domestication and improvement of Upland cotton (Gossypium hirsutum) remains poorly understood. Besides, the available genetic resources from cotton cultivars are limited. By applying restriction site-associated DNA marker sequencing (RAD-seq) technology to 582 tetraploid cotton accessions, we confirmed distinct genomic regions on chromosomes A06 and A08 in Upland cotton cultivar subgroups. Based on the pedigree, reported QTLs, introgression analyses, and genome-wide association study (GWAS), we suggest that these divergent regions might have resulted from the introgression of exotic lineages of G. hirsutum landraces and their wild relatives. These regions were the typical genomic signatures that might be responsible for maturity and fiber quality on chromosome A06 and chromosome A08, respectively. Moreover, these genomic regions are located in the putative pericentromeric regions, implying that their application will be challenging. In the study, based on high-density SNP markers, we reported two genomic signatures on chromosomes A06 and A08, which might originate from the introgression events in the Upland cotton population. Our study provides new insights for understanding the impact of historic introgressions on population divergence and important agronomic traits of modern Upland cotton cultivars.

Genomic Insight into Differentiation and Selection Sweeps in the Improvement of Upland Cotton.

Upland cotton is the most economically important fibre crop. The human-mediated selection has resulted in modern upland cultivars with higher yield and better fibre quality. However, changes in genome structure resulted from human-mediated selection are poorly understood. Comparative population genomics offers us tools to dissect the genetic history of domestication and helps to understand the genome-wide effects of human-mediated selection. Hereby, we report a comprehensive assessment of Gossypium hirsutum landraces, obsolete cultivars and modern cultivars based on high throughput genome-wide sequencing of the core set of genotypes. As a result of the genome-wide scan, we identified 93 differential regions and 311 selection sweeps associated with domestication and improvement. Furthermore, we performed genome-wide association studies to identify traits associated with the differential regions and selection sweeps. Our study provides a genetic basis to understand the domestication process in Chinese cotton cultivars. It also provides a comprehensive insight into changes in genome structure due to selection and improvement during the last century. We also identified multiple genome-wide associations (GWAS associations) for fibre yield, quality and other morphological characteristics.