Search

Domestication in the cotton genus is remarkable in that it has occurred independently four different times at two different ploidy levels. Relatively little is known about genome evolution and domestication in the cultivated diploid species Gossypium herbaceum and Gossypium arboreum, due to the absence of wild representatives for the latter species, their ancient domestication, and their joint history of human-mediated dispersal and interspecific gene flow. Using in-depth resequencing of a broad sampling from both species, we provide support for their independent domestication, as opposed to a progenitor–derivative relationship, showing that diversity (mean π = 6 × 10−3) within species is similar, and that divergence between species is modest (FST = 0.413). Individual accessions were homozygous for ancestral single-nucleotide polymorphisms at over half of variable sites, while fixed, derived sites were at modest frequencies. Notably, two chromosomes with a paucity of fixed, derived sites (i.e., chromosomes 7 and 10) were also strongly implicated as having experienced high levels of introgression. Collectively, these data demonstrate variable permeability to introgression among chromosomes, which we propose is due to divergent selection under domestication and/or the phenomenon of F2 breakdown in interspecific crosses. Our analyses provide insight into the evolutionary forces that shape diversity and divergence in the diploid cultivated species and establish a foundation for understanding the contribution of introgression and/or strong parallel selection to the extensive morphological similarities shared between species.

Rice (Oryza sativa), a major staple throughout the world and a model system for plant genomics and breeding, was the first crop genome sequenced almost two decades ago. However, reference genomes for all higher organisms to date contain gaps and missing sequences. Here, we report the assembly and analysis of gap-free reference genome sequences for two elite O. sativa xian/indica rice varieties, Zhenshan 97 and Minghui 63, which are being used as a model system for studying heterosis and yield. Gap-free reference genomes provide the opportunity for a global view of the structure and function of centromeres. We show that all rice centromeric regions share conserved centromere-specific satellite motifs with different copy numbers and structures. In addition, the similarity of CentO repeats in the same chromosome is higher than across chromosomes, supporting a model of local expansion and homogenization. Both genomes have over 395 non-TE genes located in centromere regions, of which ∼41% are actively transcribed. Two large structural variants at the end of chromosome 11 affect the copy number of resistance genes between the two genomes. The availability of the two gap-free genomes lays a solid foundation for further understanding genome structure and function in plants and breeding climate-resilient varieties.

The evolution and domestication of cotton is of great interest from both economic and evolutionary standpoints. Although many genetic and genomic resources have been generated for cotton, the genetic underpinnings of the transition from wild to domesticated cotton remain poorly known. Here we generated an intraspecific QTL mapping population specifically targeting domesticated cotton phenotypes. We used 466 F2 individuals derived from an intraspecific cross between the wild Gossypium hirsutum var. yucatanense (TX2094) and the elite cultivar G. hirsutum cv. Acala Maxxa, in two environments, to identify 120 QTL associated with phenotypic changes under domestication. While the number of QTL recovered in each subpopulation was similar, only 22 QTL were considered coincident (i.e., shared) between the two locations, eight of which shared peak markers. Although approximately half of QTL were located in the A-subgenome, many key fiber QTL were detected in the D-subgenome, which was derived from a species with unspinnable fiber. We found that many QTL are environment-specific, with few shared between the two environments, indicating that QTL associated with G. hirsutum domestication are genomically clustered but environmentally labile. Possible candidate genes were recovered and are discussed in the context of the phenotype. We conclude that the evolutionary forces that shape intraspecific divergence and domestication in cotton are complex, and that phenotypic transformations likely involved multiple interacting and environmentally responsive factors.

Domestication in the cotton genus is remarkable in that it has occurred independently four different times at two different ploidy levels. Relatively little is known about genome evolution and domestication in the cultivated diploid species Gossypium herbaceum and G. arboreum, because of the absence of wild representatives for the latter species, their ancient domestication, and their joint history of human-mediated dispersal and interspecific gene flow. Using in-depth resequencing of a broad sampling from both species, we confirm their independent domestication, as opposed to a progenitor-derivative relationship, showing that diversity (mean π = 2.3×10-3) within species is similar, and that divergence between species is modest (weighted FST=0.4430). Individual accessions were homozygous for ancestral SNPs at over half of variable sites, while fixed, derived sites were at modest frequencies. Notably, two chromosomes with a paucity of fixed, derived sites (i.e., chromosomes 7 and 10) were also strongly implicated in introgression analyses. Collectively, these data demonstrate variable permeability to introgression among chromosomes, which we propose is due to divergent selection under domestication and/or the phenomenon of F2 breakdown in interspecific crosses. Our analyses provide insight into the evolutionary forces influencing diversity and divergence in the diploid cultivated species, and establish a foundation for understanding the contribution of introgression and/or strong parallel selection to the extensive morphological similarities shared between species.Significance statementThe cotton genus (Gossypium) contains four different species that were independently domesticated at least 4,000 years ago. Relatively little is understood about diversity and evolution in the two diploid African-Asian sister-species G. herbaceum and G. arboreum, despite their historical importance in the region and contemporary cultivation, largely in the Indian subcontinent. Here we address questions regarding the relationship between the two species, their contemporary levels of diversity, and their patterns of interspecific gene flow accompanying their several millennia history of human-mediated dispersal and contact. We validate independent domestication of the two species and document the genomic distribution of interspecific genetic exchange.

Sap-sucking insects cause severe damage to cotton production. Long non-coding RNAs (lncRNAs) play vital regulatory roles in various development processes and stress response, however, the function of lncRNAs during sap-sucking insect infection in cotton is largely unknown. In this study, the transcriptome profiles between resistant (HR) and susceptible (ZS) cotton cultivars under whitefly infestation at different time points (0, 4, 12, 24, and 48 h) were compared. A total of 6,651 lncRNAs transcript and 606 differentially expressed lncRNAs were identified from the RNA-seq data. A co-expression network indicated that lncA07 and lncD09 were potential hub genes that play a regulatory role in cotton defense against aphid infestation. Furthermore, CRISPR/Cas9 knock-out mutant of lncD09 and lncA07 showed a decrease of jasmonic acid (JA) content, which potentially lead to increased susceptibility toward insect infestation. Differentially expressed genes between wild type and lncRNA knock-out plants are enriched in modulating development and resistance to stimulus. Additionally, some candidate genes such as Ghir_A01G022270, Ghir_D04G014430, and Ghir_A01G022270 are involved in the regulation of the JA-mediated signaling pathway. This result provides a novel insight of the lncRNA role in the cotton defense system against pests.

Cotton is an important crop that has been the beneficiary of multiple genome sequencing efforts, including diverse representatives of wild species for germplasm development. Gossypium anomalum is a wild African diploid species that harbors stress-resistance and fiber-related traits with potential application to modern breeding efforts. In addition, this species is a natural source of cytoplasmic male sterility and a resource for understanding hybrid lethality in the genus. Here we report a high-quality de novo genome assembly for G. anomalum and characterize this genome relative to existing genome sequences in cotton. In addition, we use the synthetic allopolyploids 2(A2D1) and 2(A2D3) to discover regions in the G. anomalum genome potentially involved in hybrid lethality, a possibility enabled by introgression of regions homologous to the D3 (G. davidsonii) lethality loci into the synthetic 2(A2D3) allopolyploid.

The chloroplast is one of two organelles containing a separate genome that codes for essential and distinct cellular functions such as photosynthesis. Given the importance of chloroplasts in plant metabolism, the genomic architecture and gene content have been strongly conserved through long periods of time and as such are useful molecular tools for evolutionary inferences. At present, complete chloroplast genomes from over 4000 species have been deposited into publicly accessible databases. Despite the large number of complete chloroplast genomes, comprehensive analyses regarding genome architecture and gene content have not been conducted for many lineages with complete species sampling. In this study, we employed the genus Populus to assess how more comprehensively sampled chloroplast genome analyses can be used in understanding chloroplast evolution in a broadly studied lineage of angiosperms. We conducted comparative analyses across Populus in order to elucidate variation in key genome features such as genome size, gene number, gene content, repeat type and number, SSR (Simple Sequence Repeat) abundance, and boundary positioning between the four main units of the genome. We found that some genome annotations were variable across the genus owing in part from errors in assembly or data checking and from this provided corrected annotations. We also employed complete chloroplast genomes for phylogenetic analyses including the dating of divergence times throughout the genus. Lastly, we utilized re-sequencing data to describe the variations of pan-chloroplast genomes at the population level for P. euphratica. The analyses used in this paper provide a blueprint for the types of analyses that can be conducted with publicly available chloroplast genomes as well as methods for building upon existing datasets to improve evolutionary inference.

Cotton (Gossypium hirsutum) is an allotetraploid species and a typical thermophilic crop that can survive and grow well under temperatures up to 45°C. CRISPR/LbCpf1 (LbCas12a) is a temperature‐sensitive system for plant genome editing (Malzahn et al., 2019) and has been successfully applied in species such as rice, soybean, tobacco, maize and cotton (Lee et al., 2019; Li et al., 2018; Tang et al., 2017; Xu et al., 2019). However, the temperature sensitivity of LbCpf1 has not been tested in cotton yet, a high temperature‐resistant crop. In order to improve LbCpf1 efficiency and determine the optimum temperature for cotton genome editing, we investigated the effects of different temperatures on LbCpf1 activity and genome editing efficiency. Moreover, we created nontransgenic and glandless cotton plants with seeds free of gossypol, representing a valuable germplasm resource for cotton breeding.

BackgroundMillennia 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.ResultsWe 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 fromGossypium hirsutumandGossypium barbadensereference 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.ConclusionsThis 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.

Cotton is an important textile crop whose gains in production over the last century have been challenged by various diseases. Because many modern cultivars are susceptible to several pests and pathogens, breeding efforts have included attempts to introgress wild, naturally resistant germplasm into elite lines. Gossypium stocksii is a wild cotton species native to Africa, which is part of a clade of vastly understudied species. Most of what is known about this species comes from pest resistance surveys and/or breeding efforts, which suggests that G. stocksii could be a valuable reservoir of natural pest resistance. Here we present a high-quality de novo genome sequence for G. stocksii. We compare the G. stocksii genome with resequencing data from a closely related, understudied species (G. somalense) to generate insight into the relatedness of these cotton species. Finally, we discuss the utility of the G. stocksii genome for understanding pest resistance in cotton, particularly resistance to cotton leaf curl virus.

Salicylic acid (SA) and brassinosteroids (BRs) are well known to regulate diverse processes of plant development and stress responses, but the mechanisms by which these phytohormones mediate the growth and defense trade-off are largely unclear. In addition, little is known about the roles of DEHYDRATION RESPONSIVE ELEMENT BINDING transcription factors, especially in biotic stress and plant growth. Here, we identified a cotton (Gossypium hirsutum) APETALA2/ETHYLENE RESPONSIVE FACTOR gene GhTINY2 that is strongly induced by Verticillium dahliae. Overexpression of GhTINY2 in cotton and Arabidopsis enhanced tolerance to V. dahliae, while knockdown of expression increased the susceptibility of cotton to the pathogen. GhTINY2 was found to promote SA accumulation and SA signaling transduction by directly activating expression of WRKY51. Moreover, GhTINY2-overexpressing cotton and Arabidopsis showed retardation of growth, increased sensitivity to inhibitors of BR biosynthesis, down-regulation of several BR-induced genes, and up-regulation of BR-repressed genes, while GhTINY2-RNAi cotton showed the opposite effects. We further determined that GhTINY2 negatively regulates BR signaling by interacting with BRASSINAZOLE-RESISTANT 1 (BZR1) and restraining its transcriptional activation of the expression of INDOLE-3-ACETIC ACID INDUCIBLE 19 (IAA19). These findings indicate that GhTINY2 fine-tunes the trade-off between immunity and growth via indirect crosstalk between WRKY51-mediated SA biosynthesis and BZR1-IAA19-regulated BR signaling.

Additional file 2. Figures S1-S30.

Additional file 3. Supplementary Notes.

Additional file 4. Review history.

We present an overview of Gossypium genome evolution and the implications of this understanding for targeted breeding objectives. The cotton genus (Gossypium) contains more than 50 species distributed in arid to semiarid regions of the tropic and subtropics. Following the genus origin approximately 10–15 million years ago, a rapid global radiation leads to eight major genome groups (A through G and K) of diploids (n = 13). Allopolyploid cottons appeared within the last 1–2 million years, as a consequence of transoceanic dispersal of an A-genome taxon to the New World and subsequent hybridization with an indigenous D-genome diploid. The nascent allopolyploid radiated into three modern lineages of seven described species, including the agronomically important species G. hirsutum L. and G. barbadense L. These two allopolyploids, together with two A-genome diploids from Africa-Asia, G. arboreum L. and G. herbaceum L., were independently domesticated for their seed fiber, representing a remarkable case of human-driven parallel evolution. Recent investigations have clarified many aspects of this evolutionary history of Gossypium with genomic insights, including the paleopolyploid history of diploid species, a surprisingly high frequency of natural interspecific hybridization within and among genome groups, myriad interactions of molecular mechanisms underlying allopolyploid genome evolution, and a much-refined evolutionary relationship among gene pools of each of the four cultivated species. The extraordinary natural diversity in Gossypium in fiber morphology, stress tolerance, and other agronomic characteristics provides ample resources for breeders to develop new cotton varieties.

Rice (Oryza sativa), a major staple throughout the world and a model system for plant genomics and breeding, was the first crop genome completed almost two decades ago. However, all sequenced genomes to date contain gaps and missing sequences. Here, we report, for the first time, the assembly and analyses of two gap-free reference genome sequences of the eliteO. sativa xian/indicarice varieties ‘Zhenshan 97 (ZS97)’ and ‘Minghui 63 (MH63)’ that are being used as a model system to study heterosis. Gap-free reference genomes also provide global insights into the structure and function of centromeres. All rice centromeric regions share conserved centromere-specific satellite motifs but with different copy numbers and structures. Importantly, we demonstrate that >1,500 genes are located in centromere regions, of which ~15.6% are actively transcribed. The generation and release of both the ZS97 and MH63 gap-free genomes lays a solid foundation for the comprehensive study of genome structure and function in plants and breed climate resilient varieties for the 21stcentury.

Background Heterosis has been exploited for decades in different crops for the dramatic increase in yield but less documented molecular evidences reported in cotton. Results An elite cotton hybrid variety ‘Huaza Mian H318’ (H318) and its parental lines was used to explore the source of its yield heterosis. A four years’ investigation of yield related traits showed that the boll number of the H318 was stable higher compare with its two parents whether in the climate suitable or unsuitable years. Besides, the hybrid H318 grows faster and showed higher fresh and dry weight than its parental lines at the seedling stage. Transcriptome analysis of seedlings identified 17,308 differential expressed genes (DEGs) between H318 and its parental lines, and 945 extremely changed DEGs were screened out for later analysis. In which, most DEGs (931/945) were gathered between H318 and its paternal line (4–5), but only 78 DEGs were found between H318 and its maternal line (B0011), which implied that H318 display more similar gene action with its maternal parent. GO and KEGG analysis showed that these DEGs highly enriched in photosynthesis and carbohydrate metabolic processes, the expression level of these DEGs was significantly higher in the H318 relative to its parental lines, which implies that sucrose and starch biosynthesis were enhanced in H318. Conclusion The enhanced photosynthesis and carbohydrate metabolic capability contributes to the heterosis of H318 at seedling stage, and establish the material foundation for later higher boll setting rates in the complex field environments.

Cotton is an important crop that has made significant gains in production over the last century. Emerging pests such as the reniform nematode have threatened cotton production. The rare African diploid species Gossypium longicalyx is a wild species that has been used as an important source of reniform nematode immunity. While mapping and breeding efforts have made some strides in transferring this immunity to the cultivated polyploid species, the complexities of interploidal transfer combined with substantial linkage drag have inhibited progress in this area. Moreover, this species shares its most recent common ancestor with the cultivated A-genome diploid cottons, thereby providing insight into the evolution of long, spinnable fiber. Here we report a newly generated de novo genome assembly of G. longicalyx. This high-quality genome leveraged a combination of PacBio long-read technology, Hi-C chromatin conformation capture, and BioNano optical mapping to achieve a chromosome level assembly. The utility of the G. longicalyx genome for understanding reniform immunity and fiber evolution is discussed.

Cis and trans regulatory divergence underlies phenotypic and evolutionary diversification. Relatively little is understood about the complexity of regulatory evolution accompanying crop domestication, particularly for polyploid plants. Here, we compare the fiber transcriptomes between wild and domesticated cotton (Gossypium hirsutum) and their reciprocal F1 hybrids, revealing genome-wide (~15%) and often compensatory cis and trans regulatory changes under divergence and domestication. The high level of trans evolution (54%–64%) observed is likely enabled by genomic redundancy following polyploidy. Our results reveal that regulatory variation is significantly associated with sequence evolution, inheritance of parental expression patterns, co-expression gene network properties, and genomic loci responsible for domestication traits. With respect to regulatory evolution, the two subgenomes of allotetraploid cotton are often uncoupled. Overall, our work underscores the complexity of regulatory evolution during fiber domestication and may facilitate new approaches for improving cotton and other polyploid plants., Relatively little is known about the complexity of regulatory evolution accompanying polyploid crop domestication. Here, using reciprocal hybrids between wild and domesticated allotetraploid cotton lines, the authors catalog cis and trans regulatory variants and show their equivalent effects on cotton fiber domestication.

Comparative population genomics offers an excellent opportunity for unraveling the genetic history of crop domestication. Upland cotton (Gossypium hirsutum) has long been an important economic crop, but a genome-wide and evolutionary understanding of the effects of human selection is lacking. Here, we describe a variation map for 352 wild and domesticated cotton accessions. We scanned 93 domestication sweeps occupying 74 Mb of the A subgenome and 104 Mb of the D subgenome, and identified 19 candidate loci for fiber-quality-related traits through a genome-wide association study. We provide evidence showing asymmetric subgenome domestication for directional selection of long fibers. Global analyses of DNase I-hypersensitive sites and 3D genome architecture, linking functional variants to gene transcription, demonstrate the effects of domestication on cis-regulatory divergence. This study provides new insights into the evolution of gene organization, regulation and adaptation in a major crop, and should serve as a rich resource for genome-based cotton improvement.

Black pepper (Piper nigrum), dubbed the ‘King of Spices’ and ‘Black Gold’, is one of the most widely used spices. Here, we present its reference genome assembly by integrating PacBio, 10x Chromium, BioNano DLS optical mapping, and Hi-C mapping technologies. The 761.2 Mb sequences (45 scaffolds with an N50 of 29.8 Mb) are assembled into 26 pseudochromosomes. A phylogenomic analysis of representative plant genomes places magnoliids as sister to the monocots-eudicots clade and indicates that black pepper has diverged from the shared Laurales-Magnoliales lineage approximately 180 million years ago. Comparative genomic analyses reveal specific gene expansions in the glycosyltransferase, cytochrome P450, shikimate hydroxycinnamoyl transferase, lysine decarboxylase, and acyltransferase gene families. Comparative transcriptomic analyses disclose berry-specific upregulated expression in representative genes in each of these gene families. These data provide an evolutionary perspective and shed light on the metabolic processes relevant to the molecular basis of species-specific piperine biosynthesis., Black pepper (Piper nigrum) belongs to the long-isolated lineage of basal angiosperm and its fruit has been used for food spice and phytomedicines for thousands of years. Here, the authors assemble the reference genome of this species and analyze gene families associated with piperine biosynthesis.

The evolution and domestication of cotton is of great interest from both economic and evolutionary standpoints. Although many genetic and genomic resources have been generated for cotton, the genetic underpinnings of the transition from wild to domesticated cotton remain poorly known. Here we generated an intraspecific QTL mapping population specifically targeting domesticated cotton phenotypes. We used 466 F2individuals derived from an intraspecific cross between the wildGossypium hirsutumvar.yucatanense(TX2094) and the elite cultivarG. hirsutumcv. Acala Maxxa, in two environments, to identify 120 QTL associated with phenotypic changes under domestication. While the number of QTL recovered in each subpopulation was similar, only 22 QTL were considered coincident (i.e., shared) between the two locations, eight of which shared peak markers. Although approximately half of QTL were located in the A-subgenome, many key fiber QTL were detected in the D-subgenome, which was derived from a species with unspinnable fiber. We found that many QTL are environment-specific, with few shared between the two environments, indicating that QTL associated withG. hirsutumdomestication are genomically clustered but environmentally labile. Possible candidate genes were recovered and are discussed in the context of the phenotype. We conclude that the evolutionary forces that shape intraspecific divergence and domestication in cotton are complex, and that phenotypic transformations likely involved multiple interacting and environmentally responsive factors.SummaryAn F2population between wild and domesticated cotton was used to identify QTL associated with selection under domestication. Multiple traits characterizing domesticated cotton were evaluated, and candidate genes underlying QTL are described for all traits. QTL are unevenly distributed between subgenomes of the domesticated polyploid, with many fiber QTL located on the genome derived from the D parent, which does not have spinnable fiber, but a majority of QTL overall located on the A subgenome. QTL are many (120) and environmentally labile. These data, together with candidate gene analyses, suggest recruitment of many environmentally responsive factors during cotton domestication.

Epigenetic modifications are highlighted for their great importance in regulating plant development, but their function associated with single-cell differentiation remains undetermined. Here, we used the cotton fibre, which is the epidermal hair on the cotton ovule, as a model to investigate the regulatory role of DNA methylation in cell differentiation. The level of CHH (H = A, T, or C) DNA methylation level was found to increase during fibre development, accompanied by a decrease in RNA-directed DNA methylation (RdDM). Examination of nucleosome positioning revealed a gradual transition from euchromatin to heterochromatin for chromatin dynamics in developing fibres, which could shape the DNA methylation landscape. The observed increase in DNA methylation in fibres, compared with other ovule tissue, was demonstrated to be mediated predominantly by an active H3K9me2-dependent pathway rather than the RdDM pathway, which was inactive. Furthermore, integrated multi-omics analyses revealed that dynamic DNA methylation played a role in the regulation of lipid biosynthesis and spatio-temporal modulation of reactive oxygen species during fibre differentiation. Our study illustrates two divergent pathways mediating a continuous increase of DNA methylation and also sheds further light on the epigenetic basis for single-cell differentiation in plants. These data and analyses are made available to the wider research community through a comprehensive web portal.

Allotetraploid cotton species (Gossypium hirsutum and Gossypium barbadense) have long been cultivated worldwide for natural renewable textile fibers. The draft genome sequences of both species are available but they are highly fragmented and incomplete1-4. Here we report reference-grade genome assemblies and annotations for G. hirsutum accession Texas Marker-1 (TM-1) and G. barbadense accession 3-79 by integrating single-molecule real-time sequencing, BioNano optical mapping and high-throughput chromosome conformation capture techniques. Compared with previous assembled draft genomes1,3, these genome sequences show considerable improvements in contiguity and completeness for regions with high content of repeats such as centromeres. Comparative genomics analyses identify extensive structural variations that probably occurred after polyploidization, highlighted by large paracentric/pericentric inversions in 14 chromosomes. We constructed an introgression line population to introduce favorable chromosome segments from G. barbadense to G. hirsutum, allowing us to identify 13 quantitative trait loci associated with superior fiber quality. These resources will accelerate evolutionary and functional genomic studies in cotton and inform future breeding programs for fiber improvement.