Comparative genome analyses highlight transposon-mediated genome expansion and the evolutionary architecture of 3D genomic folding in cotton

Chromosome-level genome assembly and annotation of a potential model organism Gossypium arboreum ZB-1

Recent advancements in plant regeneration and synthetic polyploid creation have been documented in Gossypium arboreum ZB-1. These developments make ZB-1 a potential model within the Gossypium genus for investigating gene function and polyploidy. This work generated the sequence and annotation of the ZB-1 genome. The contig-level genome was constructed using the PacBio high-fidelity reads, encompassing 81 contigs with an N50 length of 112.12 Mb. The Hi-C data assisted the construction of the chromosome-level genome, which consists of 13 pseudo-chromosomes and 39 un-anchored contigs, with a total length of about 1.67 Gb. Repetitive sequences accounted for about 69.7% of the genome in length. Based on ab initio and evidence-based prediction, we have identified 48,021 protein-coding genes in the ZB-1 genome. Comparative genomics analysis revealed conserved gene content and arrangement between ZB-1 and G. arboreum SXY1. The single nucleotide polymorphism occurrence rate between ZB-1 and SXY1 was about 0.54 per 1,000 nucleotides. This study enriched the genomic resources for further exploration into cotton regeneration and polyploidy mechanisms.

Evolutionary insights into the organization of chromatin structure and landscape of transcriptional regulation in plants

Development of complex traits necessitates the functioning and coordination of intricate regulatory networks involving multiple genes. Understanding 3D chromatin structure can facilitate insight into the regulation of gene expression by regulatory elements. This potential, of visualizing the role of chromatin organization in the evolution and function of regulatory elements, remains largely unexplored. Here, we describe new perspectives that arise from the dual considerations of sequence variation of regulatory elements and chromatin structure, with a special focus on whole-genome doubling or polyploidy. We underscore the significance of hierarchical chromatin organization in gene regulation during evolution. In addition, we describe strategies for exploring chromatin organization in future investigations of regulatory evolution in plants, enabling insights into the evolutionary influence of regulatory elements on gene expression and, hence, phenotypes.

Cotton transposon-related variome reveals roles of transposon-related variations in modern cotton cultivation

INTRODUCTION

Transposon plays a vital role in cotton genome evolution, contributing to the expansion and divergence of genomes within the Gossypium genus. However, knowledge of transposon activity in modern cotton cultivation is limited.

OBJECTIVES

In this study, we aimed to construct transposon-related variome within Gossypium genus and reveal role of transposon-related variations during cotton cultivation. In addition, we try to identify valuable transposon-related variations for cotton breeding.

METHODS

We utilized graphical genome construction to build up the graphical transposon-related variome. Based on the graphical variome, we integrated t-test, eQTL analysis and Mendelian Randomization (MR) to identify valuable transposon activities and elite genes. In addition, a convolutional neural network (CNN) model was constructed to evaluate epigenomic effects of transposon-related variations.

RESULTS

We identified 35,980 transposon activities among 10 cotton genomes, and the diversity of genomic and epigenomic features was observed among 21 transposon categories. The graphical cotton transposon-related variome was constructed, and 9,614 transposon-related variations with plasticity in the modern cotton cohort were used for eQTL, phenotypic t-test and Mendelian Randomization. 128 genes were identified as gene resources improving fiber length and strength simultaneously. 4 genes were selected from 128 genes to construct the elite gene panel whose utility has been validated in a natural cotton cohort and 2 accessions with phenotypic divergence. Based on the eQTL analysis results, we identified transposon-related variations involved in cotton's environmental adaption and human domestication, providing evidence of their role in cotton's adaption-domestication cooperation.

CONCLUSIONS

The cotton transposon-related variome revealed the role of transposon-related variations in modern cotton cultivation, providing genomic resources for cotton molecular breeding.

High-quality genomes of Bombax ceiba and Ceiba pentandra provide insights into the evolution of Malvaceae species and differences in their natural fiber development

Members of the Malvaceae family, including Corchorus spp., Gossypium spp., Bombax spp., and Ceiba spp., are important sources of natural fibers. In the past decade, the genomes of several Malvaceae species have been assembled; however, the evolutionary history of Malvaceae species and the differences in their fiber development remain to be clarified. Here, we report the genome assembly and annotation of two natural fiber plants from the Malvaceae, Bombax ceiba and Ceiba pentandra, whose assembled genome sizes are 783.56 Mb and 1575.47 Mb, respectively. Comparative analysis revealed that whole-genome duplication and Gypsy long terminal repeat retroelements have been the major causes of differences in chromosome number (2n = 14 to 2n = 96) and genome size (234 Mb to 2676 Mb) among Malvaceae species. We also used comparative genomic analyses to reconstruct the ancestral Malvaceae karyotype with 11 proto-chromosomes, providing new insights into the evolutionary trajectories of Malvaceae species. MYB-MIXTA-like 3 is relatively conserved among the Malvaceae and functions in fiber cell-fate determination in the epidermis. It appears to perform this function in any tissue where it is expressed, i.e. in fibers on the endocarp of B. ceiba and in ovule fibers of cotton. We identified a structural variation in a cellulose synthase gene and a higher copy number of cellulose synthase-like genes as possible causes of the finer, less spinnable, weaker fibers of B. ceiba. Our study provides two high-quality genomes of natural fiber plants and offers insights into the evolution of Malvaceae species and differences in their natural fiber formation and development through multi-omics analysis.

Pangenome analysis reveals transposon-driven genome evolution in cotton

BACKGROUND

Transposable elements (TEs) have a profound influence on the trajectory of plant evolution, driving genome expansion and catalyzing phenotypic diversification. The pangenome, a comprehensive genetic pool encompassing all variations within a species, serves as an invaluable tool, unaffected by the confounding factors of intraspecific diversity. This allows for a more nuanced exploration of plant TE evolution.

RESULTS

Here, we constructed a pangenome for diploid A-genome cotton using 344 accessions from representative geographical regions, including 223 from China as the main component. We found 511 Mb of non-reference sequences (NRSs) and revealed the presence of 5479 previously undiscovered protein-coding genes. Our comprehensive approach enabled us to decipher the genetic underpinnings of the distinct geographic distributions of cotton. Notably, we identified 3301 presence-absence variations (PAVs) that are closely tied to gene expression patterns within the pangenome, among which 2342 novel expression quantitative trait loci (eQTLs) were found residing in NRSs. Our investigation also unveiled contrasting patterns of transposon proliferation between diploid and tetraploid cotton, with long terminal repeat (LTR) retrotransposons exhibiting a synchronized surge in polyploids. Furthermore, the invasion of LTR retrotransposons from the A subgenome to the D subgenome triggered a substantial expansion of the latter following polyploidization. In addition, we found that TE insertions were responsible for the loss of 36.2% of species-specific genes, as well as the generation of entirely new species-specific genes.

CONCLUSIONS

Our pangenome analyses provide new insights into cotton genomics and subgenome dynamics after polyploidization and demonstrate the power of pangenome approaches for elucidating transposon impacts and genome evolution.

Implications of the three-dimensional chromatin organization for genome evolution in a fungal plant pathogen

The spatial organization of eukaryotic genomes is linked to their biological functions, although it is not clear how this impacts the overall evolution of a genome. Here, we uncover the three-dimensional (3D) genome organization of the phytopathogen Verticillium dahliae, known to possess distinct genomic regions, designated adaptive genomic regions (AGRs), enriched in transposable elements and genes that mediate host infection. Short-range DNA interactions form clear topologically associating domains (TADs) with gene-rich boundaries that show reduced levels of gene expression and reduced genomic variation. Intriguingly, TADs are less clearly insulated in AGRs than in the core genome. At a global scale, the genome contains bipartite long-range interactions, particularly enriched for AGRs and more generally containing segmental duplications. Notably, the patterns observed for V. dahliae are also present in other Verticillium species. Thus, our analysis links 3D genome organization to evolutionary features conserved throughout the Verticillium genus.

High-quality Gossypium hirsutum and Gossypium barbadense genome assemblies reveal the landscape and evolution of centromeres

Centromere positioning and organization are crucial for genome evolution; however, research on centromere biology is largely influenced by the quality of available genome assemblies. Here, we combined Oxford Nanopore and Pacific Biosciences technologies to de novo assemble two high-quality reference genomes for Gossypium hirsutum (TM-1) and Gossypium barbadense (3-79). Compared with previously published reference genomes, our assemblies show substantial improvements, with the contig N50 improved by 4.6-fold and 5.6-fold, respectively, and thus represent the most complete cotton genomes to date. These high-quality reference genomes enable us to characterize 14 and 5 complete centromeric regions for G. hirsutum and G. barbadense, respectively. Our data revealed that the centromeres of allotetraploid cotton are occupied by members of the centromeric repeat for maize (CRM) and Tekay long terminal repeat families, and the CRM family reshapes the centromere structure of the At subgenome after polyploidization. These two intertwined families have driven the convergent evolution of centromeres between the two subgenomes, ensuring centromere function and genome stability. In addition, the repositioning and high sequence divergence of centromeres between G. hirsutum and G. barbadense have contributed to speciation and centromere diversity. This study sheds light on centromere evolution in a significant crop and provides an alternative approach for exploring the evolution of polyploid plants.

Evolutionary insights into 3D genome organization and epigenetic landscape of Vigna mungo

Eukaryotic genomes show an intricate three-dimensional (3D) organization within the nucleus that regulates multiple biological processes including gene expression. Contrary to animals, understanding of 3D genome organization in plants remains at a nascent stage. Here, we investigate the evolution of 3D chromatin architecture in legumes. By using cutting-edge PacBio, Illumina, and Hi-C contact reads, we report a gap-free, chromosome-scale reference genome assembly of Vigna mungo, an important minor legume cultivated in Southeast Asia. We spatially resolved V. mungo chromosomes into euchromatic, transcriptionally active A compartment and heterochromatic, transcriptionally-dormant B compartment. We report the presence of TAD-like-regions throughout the diagonal of the HiC matrix that resembled transcriptional quiescent centers based on their genomic and epigenomic features. We observed high syntenic breakpoints but also high coverage of syntenic sequences and conserved blocks in boundary regions than in the TAD-like region domains. Our findings present unprecedented evolutionary insights into spatial 3D genome organization and epigenetic patterns and their interaction within the V. mungo genome. This will aid future genomics and epigenomics research and breeding programs of V. mungo.

Cotton pedigree genome reveals restriction of cultivar-driven strategy in cotton breeding

BACKGROUND

Many elite genes have been identified from the available cotton genomic data, providing various genetic resources for gene-driven breeding. However, backbone cultivar-driven breeding is the most widely applied strategy. Revealing the genetic basis of cultivar-driven strategy's restriction is crucial for transition of cotton breeding strategy.

RESULT

CRI12 is a backbone cultivar in cultivar-driven breeding. Here we sequence the pedigree of CRI12 using Nanopore long-read sequencing. We construct a graphical pedigree genome using the high-quality CRI12 genome and 13,138 structural variations within 20 different pedigree members. We find that low hereditary stability of elite segments in backbone cultivars is a drawback of cultivar-driven strategy. We also identify 623 functional segments in CRI12 for multiple agronomic traits in presence and absence variation-based genome-wide association study on three cohorts. We demonstrate that 25 deleterious segments are responsible for the geographical divergence of cotton in pathogen resistance. We also characterize an elite pathogen-resistant gene (GhKHCP) utilized in modern cotton breeding. In addition, we identify 386 pedigree fingerprint segments by comparing the segments of the CRI12 pedigree with those of a large cotton population.

CONCLUSION

We characterize the genetic patterns of functional segments in the pedigree of CRI12 using graphical genome method, revealing restrictions of cultivar-driven strategies in cotton breeding. These findings provide theoretical support for transitioning from cultivar-driven to gene-driven strategy in cotton breeding.

Domains Rearranged Methylase 2 maintains DNA methylation at large DNA hypomethylated shores and long-range chromatin interactions in rice

DNA methylation plays an important role in gene regulation and genomic stability. However, large DNA hypomethylated regions known as DNA methylation valleys (DMVs) or canyons have also been suggested to serve unique regulatory functions, largely unknown in rice (Oryza sativa). Here, we describe the DMVs in rice seedlings, which were highly enriched with developmental and transcription regulatory genes. Further detailed analysis indicated that grand DMVs (gDMVs) might be derived from nuclear integrants of organelle DNA (NORGs). Furthermore, Domains Rearranged Methylase 2 (OsDRM2) maintained DNA methylation at short DMV (sDMV) shores. Epigenetic maps indicated that sDMVs were marked with H3K4me3 and/or H3K27me3, although the loss of DNA methylation had a negligible effect on histone modification within these regions. In addition, we constructed H3K27me3-associated interaction maps for homozygous T-DNA insertion mutant of the gene (osdrm2) and wild type (WT). From a global perspective, most (90%) compartments were stable between osdrm2 and WT plants. At a high resolution, we observed a dramatic loss of long-range chromatin loops in osdrm2, which suffered an extensive loss of non-CG (CHG and CHH, H = A, T, or C) methylation. From another viewpoint, the loss of non-CG methylation at sDMV shores in osdrm2 could disrupt H3K27me3-mediated chromatin interaction networks. Overall, our results demonstrated that DMVs are a key genomic feature in rice and are precisely regulated by epigenetic modifications, including DNA methylation and histone modifications. OsDRM2 maintained DNA methylation at sDMV shores, while OsDRM2 deficiency strongly affected three-dimensional (3D) genome architectures.

Phylogenomic and syntenic data demonstrate complex evolutionary processes in early radiation of the rosids

Some of the most vexing problems of deep level relationship that remain in angiosperms involve the superrosids. The superrosid clade contains a quarter of all angiosperm species, with 18 orders in three subclades (Vitales, Saxifragales and core rosids) exhibiting remarkable morphological and ecological diversity. To help resolve deep-level relationships, we constructed a high-quality chromosome-level genome assembly for Tiarella polyphylla (Saxifragaceae) thus providing broader genomic representation of Saxifragales. Whole genome microsynteny analysis of superrosids showed that Saxifragales shared more synteny clusters with core rosids than Vitales, further supporting Saxifragales as more closely related with core rosids. To resolve the ordinal phylogeny of superrosids, we screened 122 single copy nuclear genes from genomes of 36 species, representing all 18 superrosid orders. Vitales were recovered as sister to all other superrosids (Saxifragales + core rosids). Our data suggest dramatic differences in relationships compared to earlier studies within core rosids. Fabids should be restricted to the nitrogen-fixing clade, while Picramniales, the Celastrales-Malpighiales (CM) clade, Huerteales, Oxalidales, Sapindales, Malvales and Brassicales formed an "expanded" malvid clade. The Celastrales-Oxalidales-Malpighiales (COM) clade (sensu APG IV) was not monophyletic. Crossosomatales, Geraniales, Myrtales and Zygophyllales did not belong to either of our well-supported malvids or fabids. There is strong discordance between nuclear and plastid phylogenetic hypotheses for superrosid relationships; we show that this is best explained by a combination of incomplete lineage sorting and ancient reticulation.

A comprehensive overview of cotton genomics, biotechnology and molecular biological studies

Cotton is an irreplaceable economic crop currently domesticated in the human world for its extremely elongated fiber cells specialized in seed epidermis, which makes it of high research and application value. To date, numerous research on cotton has navigated various aspects, from multi-genome assembly, genome editing, mechanism of fiber development, metabolite biosynthesis, and analysis to genetic breeding. Genomic and 3D genomic studies reveal the origin of cotton species and the spatiotemporal asymmetric chromatin structure in fibers. Mature multiple genome editing systems, such as CRISPR/Cas9, Cas12 (Cpf1) and cytidine base editing (CBE), have been widely used in the study of candidate genes affecting fiber development. Based on this, the cotton fiber cell development network has been preliminarily drawn. Among them, the MYB-bHLH-WDR (MBW) transcription factor complex and IAA and BR signaling pathway regulate the initiation; various plant hormones, including ethylene, mediated regulatory network and membrane protein overlap fine-regulate elongation. Multistage transcription factors targeting CesA 4, 7, and 8 specifically dominate the whole process of secondary cell wall thickening. And fluorescently labeled cytoskeletal proteins can observe real-time dynamic changes in fiber development. Furthermore, research on the synthesis of cotton secondary metabolite gossypol, resistance to diseases and insect pests, plant architecture regulation, and seed oil utilization are all conducive to finding more high-quality breeding-related genes and subsequently facilitating the cultivation of better cotton varieties. This review summarizes the paramount research achievements in cotton molecular biology over the last few decades from the above aspects, thereby enabling us to conduct a status review on the current studies of cotton and provide strong theoretical support for the future direction.

Identification and characterization of the LDAP family revealed GhLDAP2_Dt enhances drought tolerance in cotton

Lipid droplet-associated proteins (LDAPs) play essential roles in tissue growth and development and in drought stress responses in plants. Cotton is an important fiber and cash crop; however, the LDAP family has not been characterized in cotton. In this study, a total of 14, six, seven, and seven genes were confirmed as LDAP family members in Gossypium hirsutum, Gossypium raimondii, Gossypium arboreum, and Gossypium stocksii, respectively. Additionally, expansion in the LDAP family occurred with the formation of Gossypium, which is mirrored in the number of LDAPs found in five Malvaceae species (Gossypioides kirkii, Bombax ceiba, Durio zibethinus, Theobroma cacao, and Corchorus capsularis), Arabidopsis thaliana, and Carica papaya. The phylogenetic tree showed that the LDAP genes in cotton can be divided into three groups (I, II, and III). The analysis of gene structure and conserved domains showed that LDAPs derived from group I (LDAP1/2/3) are highly conserved during evolution, while members from groups II and III had large variations in both domains and gene structures. The gene expression pattern analysis of LDAP genes showed that they are expressed not only in the reproductive organs (ovule) but also in vegetative organs (root, stem, and leaves). The expression level of two genes in group III, GhLDAP6_At/Dt, were significantly higher in fiber development than in other tissues, indicating that it may be an important regulator of cotton fiber development. In group III, GhLDAP2_At/Dt, especially GhLDAP2_Dt was strongly induced by various abiotic stresses. Decreasing the expression of GhLDAP2_Dt in cotton via virus-induced gene silencing increased the drought sensitivity, and the over-expression of GhLDAP2_Dt led to increased tolerance to mannitol-simulated osmotic stress at the germination stage. Thus, we conclude that GhLDAP2_Dt plays a positive role in drought tolerance.

Identification of fertility restoration candidate genes from a restorer line R186 for Gossypium harknessii cytoplasmic male sterile cotton

BACKGROUND

The utilization of heterosis based on three-line system is an effective strategy in crop breeding. However, cloning and mechanism elucidation of restorer genes for cytoplasmic male sterility (CMS) in upland cotton have yet been realized.

RESULTS

This research is based on CMS line 2074A with the cytoplasm from Gossypium harknessii (D_2-2) and restorer line R186. The offspring of 2074A × R186 were used to conduct genetic analysis. The fertility mechanism of 2074A can be speculated to be governed by multiple genes, since neither the single gene model nor the double genes model could be used. The bulked segregant analysis (BSA) for (2074A × R186) F₂ determined the genetic interval of restorer genes on a region of 4.30 Mb on chromosome D05 that contains 77 annotated genes. Four genes were identified as candidates for fertility restoration using the RNA-seq data of 2074A, 2074B, and R186. There are a number of large effect variants in the four genes between 2074A and R186 that could cause amino acid changes. Evolutionary analysis and identity analysis revealed that GH_D05G3183, GH_D05G3384, and GH_D05G3490 have high identity with their homologs in D_2-2, respectively. Tissue differential expression analysis revealed that the genes GH_D05G3183, GH_D05G3384, and GH_D05G3490 were highly expressed in the buds of the line R186. The predicted results demonstrated that GH_D05G3183, GH_D05G3384 and GH_D05G3490 might interact with GH_A02G1295 to regulate orf610a in mitochondria.

CONCLUSION

Our study uncovered candidate genes for fertility restoration in the restorer line R186 and predicted the possible mechanism for restoring the male fertility in 2074A. This research provided valuable insight into the nucleoplasmic interactions.

Genomic rearrangements and evolutionary changes in 3D chromatin topologies in the cotton tribe (Gossypieae)

BACKGROUND

Analysis of the relationship between chromosomal structural variation (synteny breaks) and 3D-chromatin architectural changes among closely related species has the potential to reveal causes and correlates between chromosomal change and chromatin remodeling. Of note, contrary to extensive studies in animal species, the pace and pattern of chromatin architectural changes following the speciation of plants remain unexplored; moreover, there is little exploration of the occurrence of synteny breaks in the context of multiple genome topological hierarchies within the same model species.

RESULTS

Here we used Hi-C and epigenomic analyses to characterize and compare the profiles of hierarchical chromatin architectural features in representative species of the cotton tribe (Gossypieae), including Gossypium arboreum, Gossypium raimondii, and Gossypioides kirkii, which differ with respect to chromosome rearrangements. We found that (i) overall chromatin architectural territories were preserved in Gossypioides and Gossypium, which was reflected in their similar intra-chromosomal contact patterns and spatial chromosomal distributions; (ii) the non-random preferential occurrence of synteny breaks in A compartment significantly associate with the B-to-A compartment switch in syntenic blocks flanking synteny breaks; (iii) synteny changes co-localize with open-chromatin boundaries of topologically associating domains, while TAD stabilization has a greater influence on regulating orthologous expression divergence than do rearrangements; and (iv) rearranged chromosome segments largely maintain ancestral in-cis interactions.

CONCLUSIONS

Our findings provide insights into the non-random occurrence of epigenomic remodeling relative to the genomic landscape and its evolutionary and functional connections to alterations of hierarchical chromatin architecture, on a known evolutionary timescale.

Resequencing of durian genomes reveals large genetic variations among different cultivars

Durian (Durio zibethinus), which yields the fruit known as the "King of Fruits," is an important economic crop in Southeast Asia. Several durian cultivars have been developed in this region. In this study, we resequenced the genomes of three popular durian cultivars in Thailand, including Kradumthong (KD), Monthong (MT), and Puangmanee (PM) to investigate genetic diversities of cultivated durians. KD, MT, and PM genome assemblies were 832.7, 762.6, and 821.6 Mb, and their annotations covered 95.7, 92.4, and 92.7% of the embryophyta core proteins, respectively. We constructed the draft durian pangenome and analyzed comparative genomes with related species in Malvales. Long terminal repeat (LTR) sequences and protein families in durian genomes had slower evolution rates than that in cotton genomes. However, protein families with transcriptional regulation function and protein phosphorylation function involved in abiotic and biotic stress responses appeared to evolve faster in durians. The analyses of phylogenetic relationships, copy number variations (CNVs), and presence/absence variations (PAVs) suggested that the genome evolution of Thai durians was different from that of the Malaysian durian, Musang King (MK). Among the three newly sequenced genomes, the PAV and CNV profiles of disease resistance genes and the expressions of methylesterase inhibitor domain containing genes involved in flowering and fruit maturation in MT were different from those in KD and PM. These genome assemblies and their analyses provide valuable resources to gain a better understanding of the genetic diversity of cultivated durians, which may be useful for the future development of new durian cultivars.

Recent progression and future perspectives in cotton genomic breeding

Upland cotton is an important global cash crop for its long seed fibers and high edible oil and protein content. Progress in cotton genomics promotes the advancement of cotton genetics, evolutionary studies, functional genetics, and breeding, and has ushered cotton research and breeding into a new era. Here, we summarize high-impact genomics studies for cotton from the last 10 years. The diploid Gossypium arboreum and allotetraploid Gossypium hirsutum are the main focus of most genetic and genomic studies. We next review recent progress in cotton molecular biology and genetics, which builds on cotton genome sequencing efforts, population studies, and functional genomics, to provide insights into the mechanisms shaping abiotic and biotic stress tolerance, plant architecture, seed oil content, and fiber development. We also suggest the application of novel technologies and strategies to facilitate genome-based crop breeding. Explosive growth in the amount of novel genomic data, identified genes, gene modules, and pathways is now enabling researchers to utilize multidisciplinary genomics-enabled breeding strategies to cultivate "super cotton", synergistically improving multiple traits. These strategies must rise to meet urgent demands for a sustainable cotton industry.

The Gossypium herbaceum L. Wagad genome as a resource for understanding cotton domestication

Gossypium herbaceum is a species of cotton native to Africa and Asia that is one of the 2 domesticated diploids. Together with its sister-species G. arboreum, these A-genome taxa represent models of the extinct A-genome donor of modern polyploid cotton, which provide about 95% of cotton grown worldwide. As part of a larger effort to characterize variation and improve resources among diverse diploid and polyploid cotton genomes, we sequenced and assembled the genome of G. herbaceum cultivar (cv.) Wagad, representing the first domesticated accession for this species. This chromosome-level genome was generated using a combination of PacBio long-read technology, HiC, and Bionano optical mapping and compared to existing genome sequences in cotton. We compare the genome of this cultivar to the existing genome of wild G. herbaceum subspecies africanum to elucidate changes in the G. herbaceum genome concomitant with domestication and extend these analyses to gene expression using available RNA-seq. Our results demonstrate the utility of the G. herbaceum cv. Wagad genome in understanding domestication in the diploid species, which could inform modern breeding programs.

Genomic innovation and regulatory rewiring during evolution of the cotton genus Gossypium

Phenotypic diversity and evolutionary innovation ultimately trace to variation in genomic sequence and rewiring of regulatory networks. Here, we constructed a pan-genome of the Gossypium genus using ten representative diploid genomes. We document the genomic evolutionary history and the impact of lineage-specific transposon amplification on differential genome composition. The pan-3D genome reveals evolutionary connections between transposon-driven genome size variation and both higher-order chromatin structure reorganization and the rewiring of chromatin interactome. We linked changes in chromatin structures to phenotypic differences in cotton fiber and identified regulatory variations that decode the genetic basis of fiber length, the latter enabled by sequencing 1,005 transcriptomes during fiber development. We showcase how pan-genomic, pan-3D genomic and genetic regulatory data serve as a resource for delineating the evolutionary basis of spinnable cotton fiber. Our work provides insights into the evolution of genome organization and regulation and will inform cotton improvement by enabling regulome-based approaches.

Genome-wide identification and functional analysis of class III peroxidases in Gossypium hirsutum

Class III peroxidase (PRX) genes play essential roles in various processes, such as auxin catabolism, removal of H₂O₂, crosslinking cell wall components, and response to biotic and abiotic stresses. In this study, we identified 166, 78 and 89 PRX genes from G. hirsutum, G. arboretum and G. raimondii, respectively. These PRX genes were classified into seven subfamilies based on phylogenetic tree analysis and the classification of PRX genes in Arabidopsis. Segmental duplication and purifying selection were the major factors driving the evolution of GhPRXs. GO and KEGG enrichment analysis revealed that GhPRX genes were mainly associated with responding to oxidative stresses, peroxidase activities and phenylpropanoid biosynthesis pathways. Transcriptome data analysis showed that GhPRX genes expression were significantly different in microspore development between the sterility line-JinA and the maintainer line MB177. We confirmed the up-regulation of GhPRX107 and down-regulation of GhPRX128 in the sterile line compared to its maintainer line using qRT-PCR, suggesting their roles in pollen fertility. In addition, silencing GhPRX107 in cotton showed a significant decrease of the reactive oxygen species (ROS) levels of microsporocyte stage anthers compared to control. Overexpressing GhPRX107 in Arabidopsis significantly increased the ROS levels of anthers compared to wild type. In conclusion, we identified GhPRX107 as a determinant of ROS levels in anther. This work sets a foundation for PRX studies in pollen development.

The 3D architecture of the pepper genome and its relationship to function and evolution

The organization of chromatin into self-interacting domains is universal among eukaryotic genomes, though how and why they form varies considerably. Here we report a chromosome-scale reference genome assembly of pepper (Capsicum annuum) and explore its 3D organization through integrating high-resolution Hi-C maps with epigenomic, transcriptomic, and genetic variation data. Chromatin folding domains in pepper are as prominent as TADs in mammals but exhibit unique characteristics. They tend to coincide with heterochromatic regions enriched with retrotransposons and are frequently embedded in loops, which may correlate with transcription factories. Their boundaries are hotspots for chromosome rearrangements but are otherwise depleted for genetic variation. While chromatin conformation broadly affects transcription variance, it does not predict differential gene expression between tissues. Our results suggest that pepper genome organization is explained by a model of heterochromatin-driven folding promoted by transcription factories and that such spatial architecture is under structural and functional constraints.

The organization of chromatin into self-interacting domains is universal among eukaryotic genomes. Here, the authors report a reference-grade pepper genome assembly and use this reference to help describe the relationship among 3D chromatin conformation, chromatin function, and gene expression.

Spatial Features and Functional Implications of Plant 3D Genome Organization

The advent of high-throughput sequencing-based methods for chromatin conformation, accessibility, and immunoprecipitation assays has been a turning point in 3D genomics. Altogether, these new tools have been pushing upward the interpretation of pioneer cytogenetic evidence for a higher order in chromatin packing. Here, we review the latest development in our understanding of plant spatial genome structures and different levels of organization and discuss their functional implications. Then, we spotlight the complexity of organellar (i.e., mitochondria and plastids) genomes and discuss their 3D packing into nucleoids. Finally, we propose unaddressed research axes to investigate functional links between chromatin-like dynamics and transcriptional regulation within organellar nucleoids.

Inheritance, QTLs, and Candidate Genes of Lint Percentage in Upland Cotton

Cotton (Gossypium spp.) is an important natural fiber plant. Lint percentage (LP) is one of the most important determinants of cotton yield and is a typical quantitative trait with high variation and heritability. Many cotton LP genetic linkages and association maps have been reported. This work summarizes the inheritance, quantitative trait loci (QTLs), and candidate genes of LP to facilitate LP genetic study and molecular breeding. More than 1439 QTLs controlling LP have been reported. Excluding replicate QTLs, 417 unique QTLs have been identified on 26 chromosomes, including 243 QTLs identified at LOD >3. More than 60 are stable, major effective QTLs that can be used in marker-assisted selection (MAS). More than 90 candidate genes for LP have been reported. These genes encode MYB, HOX, NET, and other proteins, and most are preferentially expressed during fiber initiation and elongation. A putative molecular regulatory model of LP was constructed and provides the foundation for the genetic study and molecular breeding of LP.

Dynamic 3D genome architecture of cotton fiber reveals subgenome-coordinated chromatin topology for 4-staged single-cell differentiation

Background

Despite remarkable advances in our knowledge of epigenetically mediated transcriptional programming of cell differentiation in plants, little is known about chromatin topology and its functional implications in this process.

Results

To interrogate its significance, we establish the dynamic three-dimensional (3D) genome architecture of the allotetraploid cotton fiber, representing a typical single cell undergoing staged development in plants. We show that the subgenome-relayed switching of the chromatin compartment from active to inactive is coupled with the silencing of developmentally repressed genes, pinpointing subgenome-coordinated contribution to fiber development. We identify 10,571 topologically associating domain-like (TAD-like) structures, of which 25.6% are specifically organized in different stages and 75.23% are subject to partition or fusion between two subgenomes. Notably, dissolution of intricate TAD-like structure cliques showing long-range interactions represents a prominent characteristic at the later developmental stage. Dynamic chromatin loops are found to mediate the rewiring of gene regulatory networks that exhibit a significant difference between the two subgenomes, implicating expression bias of homologous genes.

Conclusions

This study sheds light on the spatial-temporal asymmetric chromatin structures of two subgenomes in the cotton fiber and offers a new insight into the regulatory orchestration of cell differentiation in plants.

Supplementary Information

The online version contains supplementary material available at 10.1186/s13059-022-02616-y.

The hexokinase Gene Family in Cotton: Genome-Wide Characterization and Bioinformatics Analysis

Hexokinase (HXK) is involved in hexose phosphorylation, sugar sensing, and signal transduction, all of which regulate plant growth and adaptation to stresses. Gossypium hirsutum L. is one of the most important fiber crops in the world, however, little is known about the HXKs gene family in G. hirsutum L. We identified 17 GhHXKs from the allotetraploid G. hirsutum L. genome (AADD). G. raimondii (DD) and G. arboreum (AA) are the diploid progenitors of G. hirsutum L. and contributed equally to the At_genome and Dt_genome GhHXKs genes. The chromosomal locations and exon-intron structures of GhHXK genes among cotton species are conservative. Phylogenetic analysis grouped the HXK proteins into four and three groups based on whether they were monocotyledons and dicotyledons, respectively. Duplication event analysis demonstrated that HXKs in G. hirsutum L. primarily originated from segmental duplication, which prior to diploid hybridization. Experiments of qRT-PCR, transcriptome and promoter cis-elements demonstrated that GhHXKs' promoters have auxin and GA responsive elements that are highly expressed in the fiber initiation and elongation stages, while the promoters contain ABA-, MeJA-, and SA-responsive elements that are highly expressed during the synthesis of the secondary cell wall. We performed a comprehensive analysis of the GhHXK gene family is a vital fiber crop, which lays the foundation for future studies assessing its role in fiber development.

CottonGen: The Community Database for Cotton Genomics, Genetics, and Breeding Research

Over the last eight years, the volume of whole genome, gene expression, SNP genotyping, and phenotype data generated by the cotton research community has exponentially increased. The efficient utilization/re-utilization of these complex and large datasets for knowledge discovery, translation, and application in crop improvement requires them to be curated, integrated with other types of data, and made available for access and analysis through efficient online search tools. Initiated in 2012, CottonGen is an online community database providing access to integrated peer-reviewed cotton genomic, genetic, and breeding data, and analysis tools. Used by cotton researchers worldwide, and managed by experts with crop-specific knowledge, it continuous to be the logical choice to integrate new data and provide necessary interfaces for information retrieval. The repository in CottonGen contains colleague, gene, genome, genotype, germplasm, map, marker, metabolite, phenotype, publication, QTL, species, transcriptome, and trait data curated by the CottonGen team. The number of data entries housed in CottonGen has increased dramatically, for example, since 2014 there has been an 18-fold increase in genes/mRNAs, a 23-fold increase in whole genomes, and a 372-fold increase in genotype data. New tools include a genetic map viewer, a genome browser, a synteny viewer, a metabolite pathways browser, sequence retrieval, BLAST, and a breeding information management system (BIMS), as well as various search pages for new data types. CottonGen serves as the home to the International Cotton Genome Initiative, managing its elections and serving as a communication and coordination hub for the community. With its extensive curation and integration of data and online tools, CottonGen will continue to facilitate utilization of its critical resources to empower research for cotton crop improvement.

Programmed DNA elimination: silencing genes and repetitive sequences in somatic cells

In a multicellular organism, the genomes of all cells are in general the same. Programmed DNA elimination is a notable exception to this genome constancy rule. DNA elimination removes genes and repetitive elements in the germline genome to form a reduced somatic genome in various organisms. The process of DNA elimination within an organism is highly accurate and reproducible; it typically occurs during early embryogenesis, coincident with germline-soma differentiation. DNA elimination provides a mechanism to silence selected genes and repeats in somatic cells. Recent studies in nematodes suggest that DNA elimination removes all chromosome ends, resolves sex chromosome fusions, and may also promote the birth of novel genes. Programmed DNA elimination processes are diverse among species, suggesting DNA elimination likely has evolved multiple times in different taxa. The growing list of organisms that undergo DNA elimination indicates that DNA elimination may be more widespread than previously appreciated. These various organisms will serve as complementary and comparative models to study the function, mechanism, and evolution of programmed DNA elimination in metazoans.

The Gossypium anomalum genome as a resource for cotton improvement and evolutionary analysis of hybrid incompatibility

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 (Gossypium davidsonii) lethality loci into the synthetic 2(A2D3) allopolyploid.

Disruption of topologically associating domains by structural variations in tetraploid cottons

Structural variations (SVs) are recognized to have an important role in transcriptional regulation, especially in the light of resolved 3D genome structure using high-throughput chromosome conformation capture (Hi-C) technology in mammals. However, the effect of SVs on 3D genome organization in plants remains rarely understood. In this study, we identified 295,496 SVs and 5251 topologically associating domains (TADs) in two diploid and two tetraploid cottons. We observed that approximately 16% of SVs occurred in TAD boundary regions that were called boundary affecting-structural variations (BA-SVs), and had a large effect on disrupting TAD organization. Nevertheless, SVs preferred occurring in TAD interior instead of TAD boundary, probably associated with the relaxed evolutionary selection pressure. We noticed the biased evolution of the At and Dt subgenomes of tetraploid cottons, in terms of SV-mediated disruption of 3D genome structure relative to diploids. In addition, we provide evidence showing that both SVs and TAD disruption could lead to expression difference of orthologous genes. This study advances our understanding of the effect of SVs on 3D genome organization and gene expression regulation in plants.