Genome sequence of Gossypium anomalum facilitates interspecific introgression breeding

A GDSL esterase/lipase gene GbGELP identified from a fiber micronaire QTL qMIC-A11 modulates cell elongation and fiber development

KEY MESSAGE

A fiber micronaire QTL qMIC-A11 was fine-mapped, and the GDSL esterase/lipase gene GbGELP was identified as the causal gene of the QTL. GbGELP modulates cell elongation and cotton fiber development. Fine mapping and map-based cloning of fiber micronaire (MIC)-related quantitative trait loci (QTL) have not been reported to date. Here, we utilized a G. hirsutum (Gh) acc. TM-1-G. barbadense (Gb) acc. Hai7124 introgression line CSSL47, which exhibits a significant decrease in MIC compared to TM-1, to cross with TM-1 and develop the F2 and F2:3 secondary segregating populations. Further, a stable MIC QTL qMIC-A11 was simultaneously detected in the F2 and F2:3 populations and anchored within a 407 kb region. Among them, GB_A11G1593 encoding a GDSL esterase/lipase, exhibited substantially higher expression levels at fiber elongation period in CSSL47 compared to TM-1, which was temporally identified as the causal gene for qMIC-A11 and named as GbGELP. The heterologous expression of GbGELP in Arabidopsis showed increased root length, root cell length, rosette leaf growth, and trichome density. However, knockdown of GbGELP homologs in CSSL47 significantly decreased the fiber length. Further investigation found that there was an A/T single-nucleotide polymorphism variation (SNP) in the promoter of GELP orthologs between CSSL47 and TM-1, which results in a differential CATTAAATT/CATTTAATT HAHR1-box cis-acting element, a binding site for the homeodomain-leucine zipper IV (HD-ZIP IV). GbGELP was regulated by a HD-ZIP IV transcription factor GhHDG2 via binding to the CATTAAATT element in the GbGELP promoter, while GhGELP could be activated due to GhHDG2 unable to bind the CATTTAATT element in the GhGELP promoter. The fine-mapped MIC QTL qMIC-A11, along with the causal gene GbGELP, will be utilized to improve the fiber quality in cotton breeding.

Haplotype-resolved genome assembly provides insights into the genetic basis of green peach aphid resistance in peach

Green peach aphid (GPA) is one of the most destructive pests of peach, threatening both growth and fruit quality. However, the mechanism underlying GPA resistance remains unclear. Here, we performed haplotype-resolved genome assembly of a GPA-resistant cultivar and identified an allele-specific expressed gene, PpNLR1, responsible for the GPA-resistant trait. A genome-wide association study (GWAS) revealed a functional 20-bp insertion or deletion (indel) in the PpNLR1 promoter, which co-segregated with the GPA-resistant trait and directly influenced promoter activity. Furthermore, jasmonate (JA) signaling, activated during GPA infestation, induced the transcription of PpERF109. This transcription factor specifically bound to the "CAAGT" motif within the GWAS-identified 20-bp insertion of the PpNLR1 promoter, resulting in allele-specific expression (ASE). Functional validation of the two alleles (PpNLR1-Hap1 and PpNLR1-Hap2) in both peach and Arabidopsis demonstrated their role in aphid resistance. Additionally, two GPA salivary proteins were identified as effectors, triggering reactive oxygen species (ROS) and activating the peach immune system in conjunction with the PpNLR1 protein. Comparative genomics and phylogenetic analysis indicated that an ∼53.6-kb genomic variation surrounding PpNLR1 underwent negative selection during peach evolution. In conclusion, the JA-mediated PpERF109-PpNLR1 module and GPA effector proteins significantly contribute to GPA resistance in peach. The novel haplotype-resolved genome assembly and identified key genes provide valuable resources for future genomic research and GPA resistance breeding in peach.

Monkeys at Rigged Typewriters: A Population and Network View of Plant Immune System Incompatibility

Immune system incompatibilities between naturally occurring genomic variants underlie many hybrid defects in plants and present a barrier for crop improvement. In this review, we approach immune system incompatibilities from pan-genomic and network perspectives. Pan-genomes offer insights into how natural variation shapes the evolutionary landscape of immune system incompatibilities, and through it, selection, polymorphisms, and recombination resistance emerge as common features that synergistically drive these incompatibilities. By contextualizing incompatibilities within the immune network, immune receptor promiscuity, complex dysregulation, and single-point failure appear to be recurrent themes of immune system defects. As geneticists break genes to investigate their function, so can we investigate broken immune systems to enrich our understanding of plant immune systems and work toward improving them.

Centromere-size reduction and chromatin state dynamics following intergenomic hybridization in cotton

Centromeres are pivotal for accurate chromosome segregation, yet their regulation and evolutionary dynamics remain poorly understood. Here, we investigate centromeres of the diploid species Gossypium anomalum (Ga, B-genome) that were transferred into tetraploid cotton G. hirsutum (Gh, AD-genome) as either an additional or integrated chromosome, as well as in synthetic allohexaploid (AABBDD) lines. We demonstrate consistent size reduction for all Ga centromeres in the Gh background. Histone modification profiling across 10 marks revealed heightened levels of both active and repressive chromatin marks within the Ga centromeres when transferred into the Gh background, particularly for H3K36me2. The centromeric histone modification perturbation extended into pericentromeric regions, with variable CENH3-binding domains consistently exhibiting a more pronounced increase in histone modification levels compared to stable centromere regions, highlighting the role of histone modification elevation in centromere dynamics. In addition, we observed enhanced chromatin accessibility and the presence of non-B-form DNA motifs, such as A-phased DNA repeats within stable centromere domains that are correlated with centromere stability. Hi-C analysis reveals a reorganized 3D chromatin architecture within the introgression line centromeres, including the formation of new topologically associating domains linked to H3K36me2 dynamics, emphasizing the importance of H3K36me2 in centromere organization. Together, these findings elucidate epigenetic mechanisms underlying centromere composition following intergenomic hybridization and allopolyploid formation, offering insights into centromere evolution in plants and its myriad epigenetic and potentially functional dimensions.

Cotton2035: From genomics research to optimized breeding

Cotton is the world's most important natural fiber crop and serves as an ideal model for studying plant genome evolution, cell differentiation, elongation, and cell wall biosynthesis. The first draft genome assembly for Gossypium raimondii, completed in 2012, marked the beginning of global efforts in studying cotton genomics. Over the past decade, the cotton research community has continued to assemble and refine the genomes for both wild and cultivated Gossypium species. With the accumulation of de novo genome assemblies and resequencing data across virous cotton populations, significant progress has been made in uncovering the genetic basis of key agronomic traits. Achieving the goal of cotton genomics-to-breeding (G2B) will require a deeper understanding of the spatiotemporal regulatory mechanisms involved in genome information storage and expression. We advocate for a cotton ENCODE project to systematically decode the functional elements and regulatory networks within the cotton genome. Technological advances, particularly on single-cell sequencing and high-resolution spatiotemporal omics, will be essential for elucidating these regulatory mechanisms. By integrating multi-omics data, genome editing tools, and artificial intelligence, these efforts will empower the genomics-driven strategies needed for future cotton G2B breeding.

Transcriptome Expression Profiling Reveals the Molecular Response to Salt Stress in Gossypium anomalum Seedlings

Some wild cotton species are remarkably tolerant to salt stress, and hence represent valuable resources for improving salt tolerance of the domesticated allotetraploid species Gossypium hirsutum L. Here, we first detected salt-induced stress changes in physiological and biochemical indexes of G. anomalum, a wild African diploid cotton species. Under 350 mmol/L NaCl treatment, the photosynthetic parameters declined significantly, whereas hydrogen peroxide (H2O2) and malondialdehyde (MDA) contents increased. Catalase (CAT), superoxide dismutase (SOD), and peroxidase (POD) activity and proline (PRO) content also significantly increased, reaching peak values at different stages of salt stress. We used RNA-Seq to characterize 15,476 differentially expressed genes in G. anomalum roots after 6, 12, 24, 72, and 144 h of salt stress. Gene Ontology enrichment analysis revealed these genes to be related to sequence-specific DNA and iron ion binding and oxidoreductase, peroxidase, antioxidant, and transferase activity; meanwhile, the top enriched pathways from the Kyoto Encyclopedia of Genes and Genomes database were plant hormone signal transduction, phenylpropanoid biosynthesis, fatty acid degradation, carotenoid biosynthesis, zeatin biosynthesis, starch and sucrose metabolism, and MAPK signaling. A total of 1231 transcription factors were found to be expressed in response to salt stress, representing ERF, MYB, WRKY, NAC, C2H2, bZIP, and HD-ZIP families. Nine candidate genes were validated by quantitative real-time PCR and their expression patterns were found to be consistent with the RNA-Seq data. These data promise to significantly advance our understanding of the molecular response to salt stress in Gossypium spp., with potential value for breeding applications.

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.

A recessive LRR-RLK gene causes hybrid breakdown in cotton

KEY MESSAGE

An LRR-RLK gene causing interspecific hybrid breakdown between Gossypium. anomalum and G. hirsutum was identified by deploying a map-based cloning strategy. The self-destructing symptoms of hybrid incompatibility in most cases are attributed to autoimmunity. The cloning of genes responsible for hybrid incompatibility in cotton is helpful to clarify the mechanisms underlying hybrid incompatibility and can break the barriers in distant hybridization. In this study, a temperature-dependent lethality was identified in CSSL11-9 (chromosome segment substitution line) with Gossypium anomalum chromosome segment on chromosome A11. Transcriptome analysis showed the differentially expressed genes related to autoimmune responses were highly enriched, suggesting that expression of CSSL11-9 plant lethal gene activated autoimmunity in the absence of any pathogen or external stimulus, inducing programmed cell death (PCD) and causing a lethal phenotype. The lethal phenotype was controlled by a pair of recessive genes and then fine mapped between JAAS3191-JAAS3050 interval, which covered 63.87 kb in G. hirsutum genome and 98.66 kb in G. anomalum. We demonstrated that an LRR-RLK gene designated as hybrid breakdown 1 (GoanoHBD1) was the causal gene underlying this locus for interspecific hybrid incompatibility between G. anomalum and G. hirsutum. Silencing this LRR-RLK gene could restore CSSL11-9 plants from a lethal to a normal phenotype. Our findings provide new insights into reproductive isolation and may benefit cotton breeding.

Deciphering the complex cotton genome for improving fiber traits and abiotic stress resilience in sustainable agriculture

BACKGROUND

Understanding the complex cotton genome is of paramount importance in devising a strategy for sustainable agriculture. Cotton is probably the most economically important cash crop known for its cellulose-rich fiber content. The cotton genome has become an ideal model for deciphering polyploidization due to its polyploidy, setting it apart from other major crops. However, the main challenge in understanding the functional and regulatory functions of many genes in cotton is still the complex cotton polyploidy genome, which is not limited to a single role. Cotton production is vulnerable to the sensitive effects of climate change, which can alter or aggravate soil, pests, and diseases. Thus, conventional plant breeding coupled with advanced technologies has led to substantial progress being made in cotton production.

GENOMICS APPROACHES IN COTTON

In the frontier areas of genomics research, cotton genomics has gained momentum accomplished by robust high-throughput sequencing platforms combined with novel computational tools to make the cotton genome more tractable. Advances in long-read sequencing have allowed for the generation of the complete set of cotton gene transcripts giving incisive scientific knowledge in cotton improvement. In contrast, the integration of the latest sequencing platforms has been used to generate multiple high-quality reference genomes in diploid and tetraploid cotton. While pan-genome and 3D genomic studies are still in the early stages in cotton, it is anticipated that rapid advances in sequencing, assembly algorithms, and analysis pipelines will have a greater impact on advanced cotton research.

CONCLUSIONS

This review article briefly compiles substantial contributions in different areas of the cotton genome, which include genome sequencing, genes, and their molecular regulatory networks in fiber development and stress tolerance mechanism. This will greatly help us in understanding the robust genomic organization which in turn will help unearth candidate genes for functionally important agronomic traits.

Diamonds in the not-so-rough: Wild relative diversity hidden in crop genomes

Crop production is becoming an increasing challenge as the global population grows and the climate changes. Modern cultivated crop species are selected for productivity under optimal growth environments and have often lost genetic variants that could allow them to adapt to diverse, and now rapidly changing, environments. These genetic variants are often present in their closest wild relatives, but so are less desirable traits. How to preserve and effectively utilize the rich genetic resources that crop wild relatives offer while avoiding detrimental variants and maladaptive genetic contributions is a central challenge for ongoing crop improvement. This Essay explores this challenge and potential paths that could lead to a solution.

Transcriptome Revealed GhPP2C43-A Negatively Regulates Salinity Tolerance in an Introgression Line from a Semi-wild Upland Cotton

Salt damage is one of the major threats to sustainable cotton production owing to the limited arable land in China mainly occupied by the production of staple food crops. Salt-stress tolerant cotton varieties are lacking in production and, the mechanisms underpinning salt-stress tolerance in cotton remain enigmatic. Here, DM37, an intraspecific introgression line from G. hirsutum race yucatanense acc TX-1046 into the G. hirsutum acc TM-1 background, was found to be highly tolerant to salt stress. Its seed germination rate and germination potential were significantly higher than the recipient TM-1 under salt stress. Physiological analysis showed DM37 had higher proline content and Peroxidase activity, as well as lower Na+/K+ ratios at the seedling stage, consistent with higher seedling survival rate after durable salt stress. Furthermore, comparative transcriptome analysis revealed that responsive patterns to salt stress in DM37 were different from TM-1. Weighted Correlation Network Analysis (WGCNA) demonstrated that co-expression modules associated with salt stress in DM37 also differed from TM-1. Out of them, GhPP2C43-A, a phosphatase gene, exhibited negative regulation of salt-stress tolerance verified by VIGS and transgenic Arabidopsis. Gene expression showed GhPP2C43-A in TM-1 was induced by durable salt stress but not in DM37 probably attributing to the variation of cis-element in its promoter, thereby being conferred different salt-stress tolerance. Our result would provide new genes/germplasms from semi-wild cotton in salt-stress tolerant cotton breeding. This study would give us new insights into the mechanisms underpinning the salt-stress tolerance in cotton.

Integrated metabolomic and transcriptomic strategies to reveal alkali-resistance mechanisms in wild soybean during post-germination growth stage

The keys to alkali-stress resistance of barren-tolerant wild soybean lay in enhanced reutilization of reserves in cotyledons as well as improved antioxidant protection and organic acid accumulation in young roots. Soil alkalization of farmlands is increasingly serious, adversely restricting crop growth and endangering food security. Here, based on integrated analysis of transcriptomics and metabolomics, we systematically investigated changes in cotyledon weight and young root growth in response to alkali stress in two ecotypes of wild soybean after germination to reveal alkali-resistance mechanisms in barren-tolerant wild soybean. Compared with barren-tolerant wild soybean, the dry weight of common wild soybean cotyledons under alkali stress decreased slowly and the length of young roots shortened. In barren-tolerant wild soybean, nitrogen-transport amino acids asparagine and glutamate decreased in cotyledons but increased in young roots, and nitrogen-compound transporter genes and genes involved in asparagine metabolism were significantly up-regulated in both cotyledons and young roots. Moreover, isocitric, succinic, and L-malic acids involved in the glyoxylate cycle significantly accumulated and the malate synthetase gene was up-regulated in barren-tolerant wild soybean cotyledons. In barren-tolerant wild soybean young roots, glutamate and glycine related to glutathione metabolism increased significantly and the glutathione reductase gene was up-regulated. Pyruvic acid and citric acid involved in pyruvate-citrate metabolism increased distinctly and genes encoding pyruvate decarboxylase and citrate synthetase were up-regulated. Integrated analysis showed that the keys to alkali-stress resistance of barren-tolerant wild soybean lay in enhanced protein decomposition, amino acid transport, and lipolysis in cotyledons as well as improved antioxidant protection and organic acid accumulation in young roots. This study provides new ideas for the exploitation and utilization of wild soybean resources.

Characterization of centromeric DNA of Gossypium anomalum reveals sequence-independent enrichment dynamics of centromeric repeats

Centromeres in eukaryotes are composed of highly repetitive DNAs, which evolve rapidly and are thought to achieve a favorable structure in mature centromeres. However, how the centromeric repeat evolves into an adaptive structure is largely unknown. We characterized the centromeric sequences of Gossypium anomalum through chromatin immunoprecipitation against CENH3 antibodies. We revealed that the G. anomalum centromeres contained only retrotransposon-like repeats but were depleted in long arrays of satellites. These retrotransposon-like centromeric repeats were present in the African-Asian and Australian lineage species, suggesting that they might have arisen in the common ancestor of these diploid species. Intriguingly, we observed a substantial increase and decrease in copy numbers among African-Asian and Australian lineages, respectively, for the retrotransposon-derived centromeric repeats without apparent structure or sequence variation in cotton. This result indicates that the sequence content is not a decisive aspect of the adaptive evolution of centromeric repeats or at least retrotransposon-like centromeric repeats. In addition, two active genes with potential roles in gametogenesis or flowering were identified in CENH3 nucleosome-binding regions. Our results provide new insights into the constitution of centromeric repetitive DNA and the adaptive evolution of centromeric repeats in plants.

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.

Integrated single-molecule real-time sequencing and RNA sequencing reveal the molecular mechanisms of salt tolerance in a novel synthesized polyploid genetic bridge between maize and its wild relatives

BACKGROUND

Tripsacum dactyloides (2n = 4x = 72) and Zea perennis (2n = 4x = 40) are tertiary gene pools of Zea mays L. and exhibit many abiotic adaptations absent in modern maize, especially salt tolerance. A previously reported allopolyploid (hereafter referred to as MTP, 2n = 74) synthesized using Zea mays, Tripsacum dactyloides, and Zea perennis has even stronger salt tolerance than Z. perennis and T. dactyloides. This allopolyploid will be a powerful genetic bridge for the genetic improvement of maize. However, the molecular mechanisms underlying its salt tolerance, as well as the key genes involved in regulating its salt tolerance, remain unclear.

RESULTS

Single-molecule real-time sequencing and RNA sequencing were used to identify the genes involved in salt tolerance and reveal the underlying molecular mechanisms. Based on the SMRT-seq results, we obtained 227,375 reference unigenes with an average length of 2300 bp; most of the unigenes were annotated to Z. mays sequences (76.5%) in the NR database. Moreover, a total of 484 and 1053 differentially expressed genes (DEGs) were identified in the leaves and roots, respectively. Functional enrichment analysis of DEGs revealed that multiple pathways responded to salt stress, including "Flavonoid biosynthesis," "Oxidoreductase activity," and "Plant hormone signal transduction" in the leaves and roots, and "Iron ion binding," "Acetyl-CoA carboxylase activity," and "Serine-type carboxypeptidase activity" in the roots. Transcription factors, such as those in the WRKY, B3-ARF, and bHLH families, and cytokinin negatively regulators negatively regulated the salt stress response. According to the results of the short time series-expression miner analysis, proteins involved in "Spliceosome" and "MAPK signal pathway" dynamically responded to salt stress as salinity changed. Protein-protein interaction analysis revealed that heat shock proteins play a role in the large interaction network regulating salt tolerance.

CONCLUSIONS

Our results reveal the molecular mechanism underlying the regulation of MTP in the response to salt stress and abundant salt-tolerance-related unigenes. These findings will aid the retrieval of lost alleles in modern maize and provide a new approach for using T. dactyloides and Z. perennis to improve maize.

Transcriptome Profiling of Gossypium anomalum Seedlings Reveals Key Regulators and Metabolic Pathways in Response to Drought Stress

Drought stress is a key limiting factor for cotton (Gossypium spp.) growth, production, development, and production worldwide. Some wild diploid cotton species are remarkably tolerant of water deficit and constitute an important reservoir for understanding the molecular mechanisms of Gossypium spp. drought tolerance and improving cultivated upland cotton. Here, we utilized RNA-Seq technology to characterize the leaf transcriptomes of a wild African diploid cotton species, Gossypium anomalum, under drought stress. A total of 12,322 differentially expressed genes (DEGs) were identified after mapping valid clean reads to the reference genome of G. anomalum, of which 1243 were commonly differentially expressed at all stages of drought stress. These genes were significantly enriched for molecular functions Gene Ontology terms related to cytoskeleton, hydrolase activity, cellular redox, and binding. Additionally, a substantial proportion of enriched biological process terms concerned cell or subcellular processes, while most in the cellular components category concerned membrane function and photosynthesis. An enrichment analysis against the Kyoto Encyclopedia of Genes and Genomes showed the top significantly enriched pathways to be photosynthesis-antenna proteins, amino sugar and nucleotide sugar metabolism, starch and sucrose metabolism, MAPK signaling pathway, glutathione metabolism, and plant hormone signal transduction. The DEGs also exhibited interestingly significant enrichments for drought stress-induced tandemly repeated genes involved in iron ion binding, oxidoreductase activity, heme binding, and other biological processes. A large number of genes encoding transcription factors, such as MYB, bHLH, ERF, NAC, WRKY, and bZIP, were identified as playing key roles in acclimatizing to drought stress. These results will provide deeper insights into the molecular mechanisms of drought stress adaptation in Gossypium spp.