Genome sequencing of the Australian wild diploid species Gossypium australe highlights disease resistance and delayed gland morphogenesis

Phytol-induced interplant signaling in maize facilitates EXP-A20-driven resistance through ACO31-dependent ethylene accumulation against Ostrinia furnacalis

Plants have evolved sophisticated defense mechanisms against insect herbivores, including cell wall fortification through lignin biosynthesis. Insect attack primes systemic acquired resistance in plants, preparing them to respond more swiftly and vigorously to subsequent insect assaults. Here, we found that Beauveria bassiana-exposed maize plants can emit phytol upon infestation by Spodoptera frugiperda, inducing plant-to-plant (PTP) communication of alert signals for neighboring plants, and revealed the expansin protein EXP-A20 as a pivotal node mediating maize defense responses in neighboring plants against the destructive pest Ostrinia furnacalis via stimulation of ethylene (ET) synthesis and lignin production. Through virus-induced gene silencing, we showed that EXP-A20 is essential for maize resistance, while downregulating ET and lignin pathways. Critically, protein-protein interactions determined via luciferase complementation and yeast two-hybrid assays demonstrated that EXP-A20 binds to and likely activates the ET-forming enzyme gene ACO31 to initiate defense signaling cascades, representing a novel signaling modality for expansins. Treatment with the plant volatile phytol has known insecticidal/priming activity, but we found that its effectiveness requires EXP-A20. This finding highlights the importance of EXP-A20 upstream of hormone-cell wall crosstalk in defense activation by volatiles. Overall, our multifaceted dissection of EXP-A20 revealed key molecular intersections underlying inducible maize immunity against herbivores. Furthermore, we provide functional evidence that extensive cell growth processes directly stimulate defense programs in plants. Our work opens new avenues for enhancing durable, broad-spectrum pest resistance in maize through the use of volatile organic compounds and PTP interactions.

Evaluation of resistance to Verticillium wilt in Gossypium hirsutum-Gossypium arboreum introgression lines and identification of putative resistance genes using RNA-seq

Verticillium wilt (VW), a fungal disease caused by Verticillium dahliae (Vd), is one of the most destructive threats to cotton production. Moreover, widely cultivated upland cotton (Gossypium hirsutum, 2n = 4x = AADD = 52) often demonstrates low resistance to Vd. In contrast, G. arboreum (2n = 2x = AA = 26) shows high resistant to VW, making it a valuable source for breeding, despite the challenges posed by hybridization incompatibility between the two species. Here, a population of introgression lines derived from G. hirsutum and G. arboreum was evaluated for resistance to VW through both glasshouse and field tests. Among these lines, DM11039 demonstrated high resistance to VW. Both DM11039 and the recipient TM-1 underwent transcriptome sequencing during Vd infection at 0, 4, 12, 24, 48, and 96 h post inoculation. The analysis identified differentially expressed genes (DEGs), which were predominantly associated with resistance mechanisms. Based on the results from transcriptome sequencing and weighted correlation network analysis, three DEGs from each parent-G. arboreum and G. hirsutum- in DM11039 were subjected to virus-induced gene silencing in cotton seedlings. The findings revealed that silencing of GaPP2A1, GaPDH-E1, or GaLRK10L-1.2, which are located within the introgression segments from G. arboreum, significantly impaired disease resistance in cotton. This suggests that these genes are potentially linked to the disease phenotype. In contrast, silencing of GHA13G1263, GhZIP1 or GHA10G2498 from G. hirsutum did not result in any changes in disease resistance in DM11039. The results indicate G. arboreum harbors resistance genes to VW. Furthermore, the introgression population presents a valuable resource for future cotton breeding.

Histone deacetylase GhHDA5 negatively regulates Verticillium wilt resistance in cotton

Verticillium wilt (VW) caused by Verticillium dahliae (V. dahliae) is one of the most destructive diseases in cotton (Gossypium spp.). Histone acetylation plays critical roles in plant development and adaptive responses to biotic and abiotic stresses. However, the relevance of histone acetylation in cotton VW resistance remains largely unclear. Here, we identified histone deacetylase 5 (GhHDA5) from upland cotton (Gossypium hirsutum L.), as a negative regulator of VW resistance. GhHDA5 expression was responsive to V. dahliae infection. Silencing GhHDA5 in upland cotton led to improved resistance to V. dahliae, while heterologous expression of GhHDA5 in Arabidopsis (Arabidopsis thaliana) compromised V. dahliae tolerance. GhHDA5 repressed the expression of several lignin biosynthesis-related genes, such as 4-coumarate:CoA ligase gene Gh4CL3 and ferulate 5-hydroxylase gene GhF5H, through reducing the acetylation level of histone H3 lysine 9 and 14 (H3K9K14ac) at their promoter regions, thereby resulting in an increased deposition of lignin, especially S monomers, in the GhHDA5-silenced cotton plants. The silencing of GhF5H impaired cotton VW tolerance. Additionally, the silencing of GhHDA5 also promoted the production of reactive oxygen species (ROS), elevated the expression of several pathogenesis-related genes (PRs), and altered the content and signaling of the phytohormones salicylic acid (SA), jasmonic acid (JA), and strigolactones (SLs) after V. dahliae infection. Taken together, our findings suggest that GhHDA5 negatively regulates cotton VW resistance through modulating disease-induced lignification and the ROS- and phytohormone-mediated defense response.

GhHAM regulates GoPGF-dependent gland development and contributes to broad-spectrum pest resistance in cotton

Cotton is a globally cultivated crop, producing 87% of the natural fiber used in the global textile industry. The pigment glands, unique to cotton and its relatives, serve as a defense structure against pests and pathogens. However, the molecular mechanism underlying gland formation and the specific role of pigment glands in cotton's pest defense are still not well understood. In this study, we cloned a gland-related transcription factor GhHAM and generated the GhHAM knockout mutant using CRISPR/Cas9. Phenotypic observations, transcriptome analysis, and promoter-binding experiments revealed that GhHAM binds to the promoter of GoPGF, regulating pigment gland formation in cotton's multiple organs via the GoPGF-GhJUB1 module. The knockout of GhHAM significantly reduced gossypol production and increased cotton's susceptibility to pests in the field. Feeding assays demonstrated that more than 80% of the cotton bollworm larvae preferred ghham over the wild type. Furthermore, the ghham mutants displayed shorter cell length and decreased gibberellins (GA) production in the stem. Exogenous application of GA3 restored stem cell elongation but not gland formation, thereby indicating that GhHAM controls gland morphogenesis independently of GA. Our study sheds light on the functional differentiation of HAM proteins among plant species, highlights the significant role of pigment glands in influencing pest feeding preference, and provides a theoretical basis for breeding pest-resistant cotton varieties to address the challenges posed by frequent outbreaks of pests.

Endophytic Fungi-Mediated Defense Signaling in Maize: Unraveling the Role of WRKY36 in Regulating Immunity against Spodoptera frugiperda

Seed priming with beneficial endophytic fungi is an emerging sustainable strategy for enhancing plant resistance against insect pests. This study examined the effects of Beauvaria bassiana Bb20091317 and Metarhizium rileyi MrCDTLJ1 fungal colonization on maize growth, defence signalling, benzoxazinoid levels and gene expression. The colonization did not adversely affect plant growth but reduced larval weights of Spodoptera frugiperda. Maize leaves treated with M. rileyi exhibited higher levels of jasmonic acid, jasmonoyl-Isoleucine, salicylic acid, and indole acetic acid compared to control. B. bassiana and M. rileyi accelerated phytohormone increase upon S. frugiperda herbivory. Gene expression analysis revealed modulation of benzoxazinoid biosynthesis genes. We further elucidated the immune regulatory role of the transcription factor zmWRKY36 using virus-induced gene silencing (VIGS) in maize. zmWRKY36 positively regulates maize immunity against S. frugiperda, likely by interacting with defense-related proteins. Transient overexpression of zmWRKY36 in tobacco-induced cell death, while silencing in maize reduced chitin-triggered reactive oxygen species burst, confirming its immune function. Overall, B. bassiana and M. rileyi successfully colonized maize, impacting larval growth, defense signalling, and zmWRKY36-mediated resistance. This sheds light on maize-endophyte-insect interactions for sustainable plant protection.

The dual role of GoPGF reveals that the pigment glands are synthetic sites of gossypol in aerial parts of cotton

Gossypol and the related terpenoids are stored in the pigment gland to protect cotton plants from biotic stresses, but little is known about the synthetic sites of these metabolites. Here, we showed that GoPGF, a key gene regulating gland formation, was expressed in gland cells and roots. The chromatin immunoprecipitation sequencing (ChIP-seq) analysis demonstrated that GoPGF targets GhJUB1 to regulate gland morphogenesis. RNA-sequencing (RNA-seq) showed high accumulation of gossypol biosynthetic genes in gland cells. Moreover, integrated analysis of the ChIP-seq and RNA-seq data revealed that GoPGF binds to the promoter of several gossypol biosynthetic genes. The cotton callus overexpressing GoPGF had dramatically increased the gossypol levels, indicating that GoPGF can directly activate the biosynthesis of gossypol. In addition, the gopgf mutant analysis revealed the existence of both GoPGF-dependent and -independent regulation of gossypol production in cotton roots. Our study revealed that the pigment glands are synthetic sites of gossypol in aerial parts of cotton and that GoPGF plays a dual role in regulating gland morphogenesis and gossypol biosynthesis. The study provides new insights for exploring the complex relationship between glands and the metabolites they store in cotton and other plant species.

General Regulatory Factor7 regulates innate immune signalling to enhance Verticillium wilt resistance in cotton

Sessile growing plants are always vulnerable to microbial pathogen attacks throughout their lives. To fend off pathogen invasion, plants have evolved a sophisticated innate immune system that consists of cell surface receptors and intracellular receptors. Somatic embryogenesis receptor kinases (SERKs) belong to a small group of leucine-rich repeat receptor-like kinases (LRR-RLKs) that function as co-receptors regulating diverse physiological processes. GENRAL REGULATORY FACTOR (GRF) proteins play an important role in physiological signalling transduction. However, the function of GRF proteins in plant innate immune signalling remains elusive. Here, we identified a GRF gene, GauGRF7, that is expressed both constitutively and in response to fungal pathogen infection. Intriguingly, silencing of GRF7 compromised plant innate immunity, resulting in susceptibility to Verticillium dahliae infection. Both transgenic GauGRF7 cotton and transgenic GauGRF7 Arabidopsis lines enhanced the innate immune response to V. dahliae infection, leading to high expression of two helper NLRs (hNLR) genes (ADR1 and NRG1) and pathogenesis-related genes, and increased ROS production and salicylic acid level. Moreover, GauGRF7 interacted with GhSERK1, which positively regulated GRF7-mediated innate immune response in cotton and Arabidopsis. Our findings revealed the molecular mechanism of the GRF protein in plant immune signaling and offer potential opportunities for improving plant resistance to V. dahliae infection.

RVE2, a new regulatory factor in jasmonic acid pathway, orchestrates resistance to Verticillium wilt

Verticillium dahliae, one of the most destructive fungal pathogens of several crops, challenges the sustainability of cotton productivity worldwide because very few widely-cultivated Upland cotton varieties are resistant to Verticillium wilt (VW). Here, we report that REVEILLE2 (RVE2), the Myb-like transcription factor, confers the novel function in resistance to VW by regulating the jasmonic acid (JA) pathway in cotton. RVE2 expression was essentially required for the activation of JA-mediated disease-resistance response. RVE2 physically interacted with TPL/TPRs and disturbed JAZ proteins to recruit TPL and TPR1 in NINJA-dependent manner, which regulated JA response by relieving inhibited-MYC2 activity. The MYC2 then bound to RVE2 promoter for the activation of its transcription, forming feedback loop. Interestingly, a unique truncated RVE2 widely existing in D-subgenome (GhRVE2D) of natural Upland cotton represses the ability of the MYC2 to activate GhRVE2A promoter but not GausRVE2 or GbRVE2. The result could partially explain why Gossypium barbadense popularly shows higher resistance than Gossypium hirsutum. Furthermore, disturbing the JA-signalling pathway resulted into the loss of RVE2-mediated disease-resistance in various plants (Arabidopsis, tobacco and cotton). RVE2 overexpression significantly enhanced the resistance to VW. Collectively, we conclude that RVE2, a new regulatory factor, plays a pivotal role in fine-tuning JA-signalling, which would improve our understanding the mechanisms underlying the resistance to VW.

Decoding the guardians of cotton resilience: A comprehensive exploration of the βCA genes and its role in Verticillium dahliae resistance

Plant Carbonic anhydrases (Cas) have been shown to be stress-responsive enzymes that may play a role in adapting to adverse conditions. Cotton is a significant economic crop in China, with upland cotton (Gossypium hirsutum) being the most widely cultivated species. We conducted genome-wide identification of the βCA gene in six cotton species and preliminary analysis of the βCA gene in upland cotton. In total, 73 βCA genes from six cotton species were identified, with phylogenetic analysis dividing them into five subgroups. GHβCA proteins were predominantly localized in the chloroplast and cytoplasm. The genes exhibited conserved motifs, with motifs 1, 2, and 3 being prominent. GHβCA genes were unevenly distributed across chromosomes and were associated with stress-responsive cis-regulatory elements, including those responding to light, MeJA, salicylic acid, abscisic acid, cell cycle regulation, and defence/stress. Expression analysis indicated that GHβCA6, GHβCA7, GHβCA10, GHβCA15, and GHβCA16 were highly expressed under various abiotic stress conditions, whereas GHβCA3, GHβCA9, GHβCA10, and GHβCA18 had higher expression patterns under Verticillium dahliae infection at different time intervals. In Gossypium thurberi, GthβCA1, GthβCA2, and GthβCA4 showed elevated expression across stress conditions and tissues. Silencing GHβCA10 through VIGS increased Verticillium wilt severity and reduced lignin deposition compared to non-silenced plants. GHβCA10 is crucial for cotton's defense against Verticillium dahliae. Further research is needed to understand the underlying mechanisms and develop strategies to enhance resistance against Verticillium wilt.

Dynamic characteristics and functional analysis provide new insights into the role of GauERF105 for resistance against Verticillium dahliae in cotton

BACKGROUND

The cotton industry suffers significant yield losses annually due to Verticillium wilt, which is considered the most destructive disease affecting the crop. However, the precise mechanisms behind this disease in cotton remain largely unexplored.

METHODS

Our approach involved utilizing transcriptome data from G. australe which was exposed to Verticillium dahliae infection. From this data, we identified ethylene-responsive factors and further investigated their potential role in resistance through functional validations via Virus-induced gene silencing (VIGS) in cotton and overexpression in Arabidopsis.

RESULTS

A total of 23 ethylene response factors (ERFs) were identified and their expression was analyzed at different time intervals (24 h, 48 h, and 72 h post-inoculation). Among them, GauERF105 was selected based on qRT-PCR expression analysis for further investigation. To demonstrate the significance of GauERF105, VIGS was utilized, revealing that suppressing GauERF105 leads to more severe infections in cotton plants compared to the wild-type. Additionally, the silenced plants exhibited reduced lignin deposition in the stems compared to the WT plants, indicating that the silencing of GauERF105 also impacts lignin content. The overexpression of GauERF105 in Arabidopsis confirmed its pivotal role in conferring resistance against Verticillium dahliae infection. Our results suggest that WT possesses higher levels of the oxidative stress markers MDA and H2O2 as compared to the overexpressed lines. In contrast, the activities of the antioxidant enzymes SOD and POD were higher in the overexpressed lines compared to the WT. Furthermore, DAB and trypan staining of the overexpressed lines suggested a greater impact of the disease in the wild-type compared to the transgenic lines.

CONCLUSIONS

Our findings provide confirmation that GauERF105 is a crucial candidate in the defense mechanism of cotton against Verticillium dahliae invasion, and plays a pivotal role in this process. These results have the potential to facilitate the development of germplasm resistance in cotton.

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.

An Integrated Multi-Omics and Artificial Intelligence Framework for Advance Plant Phenotyping in Horticulture

This review discusses the transformative potential of integrating multi-omics data and artificial intelligence (AI) in advancing horticultural research, specifically plant phenotyping. The traditional methods of plant phenotyping, while valuable, are limited in their ability to capture the complexity of plant biology. The advent of (meta-)genomics, (meta-)transcriptomics, proteomics, and metabolomics has provided an opportunity for a more comprehensive analysis. AI and machine learning (ML) techniques can effectively handle the complexity and volume of multi-omics data, providing meaningful interpretations and predictions. Reflecting the multidisciplinary nature of this area of research, in this review, readers will find a collection of state-of-the-art solutions that are key to the integration of multi-omics data and AI for phenotyping experiments in horticulture, including experimental design considerations with several technical and non-technical challenges, which are discussed along with potential solutions. The future prospects of this integration include precision horticulture, predictive breeding, improved disease and stress response management, sustainable crop management, and exploration of plant biodiversity. The integration of multi-omics and AI holds immense promise for revolutionizing horticultural research and applications, heralding a new era in plant phenotyping.

Evolutionary roles of polyploidization-derived structural variations in the phenotypic diversification of Panax species

Genomic structural variations (SVs) are widespread in plant and animal genomes and play important roles in phenotypic novelty and species adaptation. Frequent whole genome duplications followed by (re)diploidizations have resulted in high diversity of genome architecture among extant species. In this study, we identified abundant genomic SVs in the Panax genus that are hypothesized to have occurred through during the repeated polyploidizations/(re)diploidizations. Our genome-wide comparisons demonstrated that although these polyploidization-derived SVs have evolved at distinct evolutionary stages, a large number of SV-intersecting genes showed enrichment in functionally important pathways related to secondary metabolites, photosynthesis and basic cellular activities. In line with these observations, our metabolic analyses of these Panax species revealed high diversity of primary and secondary metabolites both at the tissue and interspecific levels. In particular, genomic SVs identified at ginsenoside biosynthesis genes, including copy number variation and large fragment deletion, appear to have played important roles in the evolution and diversification of ginsenosides. A further herbivore deterrence experiment demonstrated that, as major triterpenoidal saponins found exclusively in Panax, ginsenosides provide protection against insect herbivores. Our study provides new insights on how polyploidization-derived SVs have contributed to phenotypic novelty and plant adaptation.

Establishing a Virus-Induced Gene Silencing System in Lycoris chinensis

Lycoris is an important plant with both medicinal and ornamental values. However, it does not have an efficient genetic transformation system, which makes it difficult to study gene function of the genus. Virus-induced gene silencing (VIGS) is an effective technique for studying gene functions in plants. In this study, we develop an efficient virus-induced gene-silencing (VIGS) system using the leaf tip needle injection method. The widely used TRV vector is constructed, and the Cloroplastos Alterados 1 (CLA1) and Phytoene Desaturase (PDS) genes are selected as visual indicators in the VIGS system. As a result, it is observed that leaves infected with TRV-LcCLA1 and TRV-LcPDS both show a yellowing phenotype (loss of green), and the chlorosis range of TRV-LcCLA1 was larger and deeper than that of TRV-LcPDS. qRT-PCR results show that the expression levels of LcCLA1 and LcPDS are significantly reduced, and the silencing efficiency of LcCLA1 is higher than that of LcPDS. These results indicate that the VIGS system of L. chinensis was preliminarily established, and LcCLA1 is more suitable as a gene-silencing indicator. For the monocotyledonous plant leaves with a waxy surface, the leaf tip injection method greatly improves the infiltration efficiency. The newly established VIGS system will contribute to gene functional research in Lycoris species.

Genome-wide identification and functional analysis of ICE genes reveal that Gossypium thurberi "GthICE2" is responsible for cold and drought stress tolerance

Cold stress has been found to have a negative impact on cotton growth and annual production. To address this issue, the utilization of cold-tolerant gene resources from wild species of Gossypium is crucial for genetic improvements in cultivated cotton. ICE (inducer of CBF expression) are the key regulators of cold tolerance in plants, however, there is relatively little information on ICE genes in cotton. Herein, we performed comprehensive bioinformatics analyses of the ICE gene family in eight cotton species. Phylogenetic analysis showed that 52 ICE genes were clustered into four subgroups. Cis-regulatory elements analysis suggests that the expression of ICE genes might be regulated by light, plant hormones, and various environment stresses. Higher expression of GthICE2 was observed in leaves as compared to roots and stems, in response to cold, drought, and exogenous hormone ABA. Furthermore, overexpression of GthICE2 in A. thaliana led to higher germination and survival rates, longer root length, lower ion leakage, and induction under cold and drought stress. Histochemical staining showed that oxidative damage in transgenic lines was much lower compared to wild-type plants. Lower MDA contents and higher SOD and POD activities were observed in overexpressed plants. Y1H and LUC assays revealed that GthICE2 might activate the expression of GthCBF4, a cold-responsive gene, by connecting with the MYC cis-element present in the promoter of GthCBF4. GthICE2 confers cold and drought stress tolerance in cotton. Our findings add significantly to the existing knowledge regarding cold stress tolerance and helps to elucidate cold response mechanisms in cotton.

Strigolactones positively regulate Verticillium wilt resistance in cotton via crosstalk with other hormones

Verticillium wilt caused by Verticillium dahliae is a serious vascular disease in cotton (Gossypium spp.). V. dahliae induces the expression of the CAROTENOID CLEAVAGE DIOXYGENASE 7 (GauCCD7) gene involved in strigolactone (SL) biosynthesis in Gossypium australe, suggesting a role for SLs in Verticillium wilt resistance. We found that the SL analog rac-GR24 enhanced while the SL biosynthesis inhibitor TIS108 decreased cotton resistance to Verticillium wilt. Knock-down of GbCCD7 and GbCCD8b genes in island cotton (Gossypium barbadense) decreased resistance, whereas overexpression of GbCCD8b in upland cotton (Gossypium hirsutum) increased resistance to Verticillium wilt. Additionally, Arabidopsis (Arabidopsis thaliana) SL mutants defective in CCD7 and CCD8 putative orthologs were susceptible, whereas both Arabidopsis GbCCD7- and GbCCD8b-overexpressing plants were more resistant to Verticillium wilt than wild-type (WT) plants. Transcriptome analyses showed that several genes related to the jasmonic acid (JA)- and abscisic acid (ABA)-signaling pathways, such as MYELOCYTOMATOSIS 2 (GbMYC2) and ABA-INSENSITIVE 5, respectively, were upregulated in the roots of WT cotton plants in responses to rac-GR24 and V. dahliae infection but downregulated in the roots of both GbCCD7- and GbCCD8b-silenced cotton plants. Furthermore, GbMYC2 suppressed the expression of GbCCD7 and GbCCD8b by binding to their promoters, which might regulate the homeostasis of SLs in cotton through a negative feedback loop. We also found that GbCCD7- and GbCCD8b-silenced cotton plants were impaired in V. dahliae-induced reactive oxygen species (ROS) accumulation. Taken together, our results suggest that SLs positively regulate cotton resistance to Verticillium wilt through crosstalk with the JA- and ABA-signaling pathways and by inducing ROS accumulation.

Comparative Transcriptome Analysis Revealed Key Genes Regulating Gossypol Synthesis in Tetraploid Cultivated Cotton

Tetraploid cultivated cotton (Gossypium spp.) produces cottonseeds rich in protein and oil. Gossypol and related terpenoids, stored in the pigment glands of cottonseeds, are toxic to human beings and monogastric animals. However, a comprehensive understanding of the genetic basis of gossypol and gland formation is still lacking. We performed a comprehensive transcriptome analysis of four glanded versus two glandless tetraploid cultivars distributed in Gossypium hirsutum and Gossypium barbadense. A weighted gene co-expression network analysis (WGCNA) based on 431 common differentially expressed genes (DEGs) uncovered a candidate module that was strongly associated with the reduction in or disappearance of gossypol and pigment glands. Further, the co-expression network helped us to focus on 29 hub genes, which played key roles in the regulation of related genes in the candidate module. The present study contributes to our understanding of the genetic basis of gossypol and gland formation and serves as a rich potential source for breeding cotton cultivars with gossypol-rich plants and gossypol-free cottonseed, which is beneficial for improving food safety, environmental protection, and economic gains of tetraploid cultivated cotton.

Insights to Gossypium defense response against Verticillium dahliae: the Cotton Cancer

The soil-borne pathogen Verticillium dahliae, also referred as "The Cotton Cancer," is responsible for causing Verticillium wilt in cotton crops, a destructive disease with a global impact. To infect cotton plants, the pathogen employs multiple virulence mechanisms such as releasing enzymes that degrade cell walls, activating genes that contribute to virulence, and using protein effectors. Conversely, cotton plants have developed numerous defense mechanisms to combat the impact of V. dahliae. These include strengthening the cell wall by producing lignin and depositing callose, discharging reactive oxygen species, and amassing hormones related to defense. Despite the efforts to develop resistant cultivars, there is still no permanent solution to Verticillium wilt due to a limited understanding of the underlying molecular mechanisms that drive both resistance and pathogenesis is currently prevalent. To address this challenge, cutting-edge technologies such as clustered regularly interspaced short palindromic repeats (CRISPR)/CRISPR-associated protein 9 (Cas9), host-induced gene silencing (HIGS), and gene delivery via nano-carriers could be employed as effective alternatives to control the disease. This article intends to present an overview of V. dahliae virulence mechanisms and discuss the different cotton defense mechanisms against Verticillium wilt, including morphophysiological and biochemical responses and signaling pathways including jasmonic acid (JA), salicylic acid (SA), ethylene (ET), and strigolactones (SLs). Additionally, the article highlights the significance of microRNAs (miRNAs), circular RNAs (circRNAs), and long non-coding RNAs (lncRNAs) in gene expression regulation, as well as the different methods employed to identify and functionally validate genes to achieve resistance against this disease. Gaining a more profound understanding of these mechanisms could potentially result in the creation of more efficient strategies for combating Verticillium wilt in cotton crops.

Single-cell transcriptomic analysis reveals the developmental trajectory and transcriptional regulatory networks of pigment glands in Gossypium bickii

Comprehensive utilization of cottonseeds is limited by the presence of pigment glands and its inclusion gossypol. The ideal cotton has glandless seeds but a glanded plant, a trait found in only a few Australian wild cotton species, including Gossypium bickii. Introgression of this trait into cultivated species has proved to be difficult. Understanding the biological processes toward pigment gland morphogenesis and the associated underlying molecular mechanisms will facilitate breeding of cultivated cotton varieties with the trait of glandless seeds and glanded plant. In this study, single-cell RNA sequencing (scRNA-seq) was performed on 12 222 protoplasts isolated from cotyledons of germinating G. bickii seeds 48 h after imbibition. Clustered into 14 distinct clusters unsupervisedly, these cells could be grouped into eight cell populations with the assistance of known cell marker genes. The pigment gland cells were well separated from others and could be separated into pigment gland parenchyma cells, secretory cells, and apoptotic cells. By integrating the pigment gland cell developmental trajectory, transcription factor regulatory networks, and core transcription factor functional validation, we established a model for pigment gland formation. In this model, light and gibberellin were verified to promote the formation of pigment glands. In addition, three novel genes, GbiERF114 (ETHYLENE RESPONSE FACTOR 114), GbiZAT11 (ZINC FINGER OF ARABIDOPSIS THALIANA 11), and GbiNTL9 (NAC TRANSCRIPTION FACTOR-LIKE 9), were found to affect pigment gland formation. Collectively, these findings provide new insights into pigment gland morphogenesis and lay the cornerstone for future cotton scRNA-seq investigations.

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.

Comparative Genomics and Functional Studies of Putative m6A Methyltransferase (METTL) Genes in Cotton

N6-methyladenosine (m6A) RNA modification plays important regulatory roles in plant development and adapting to the environment, which requires methyltransferases to achieve the methylation process. However, there has been no research regarding m6A RNA methyltransferases in cotton. Here, a systematic analysis of the m6A methyltransferase (METTL) gene family was performed on twelve cotton species, resulting in six METTLs identified in five allotetraploid cottons, respectively, and three to four METTLs in the seven diploid species. Phylogenetic analysis of protein-coding sequences revealed that METTL genes from cottons, Arabidopsis thaliana, and Homo sapiens could be classified into three clades (METTL3, METTL14, and METTL-like clades). Cis-element analysis predicated the possible functions of METTL genes in G. hirsutum. RNA-seq data revealed that GhMETTL14 (GH_A07G0817/GH_D07G0819) and GhMETTL3 (GH_A12G2586/GH_D12G2605) had high expressions in root, stem, leaf, torus, petal, stamen, pistil, and calycle tissues. GhMETTL14 also had the highest expression in 20 and 25 dpa fiber cells, implying a potential role at the cell wall thickening stage. Suppressing GhMETTL3 and GhMETTL14 by VIGS caused growth arrest and even death in G. hirsutum, along with decreased m6A abundance from the leaf tissues of VIGS plants. Overexpression of GhMETTL3 and GhMETTL14 produced distinct differentially expressed genes (DEGs) in A. thaliana, indicating their possible divergent functions after gene duplication. Overall, GhMETTLs play indispensable but divergent roles during the growth of cotton plants, which provides the basis for the systematic investigation of m6A in subsequent studies to improve the agronomic traits in cotton.

Evolutionary divergence of duplicated genomes in newly described allotetraploid cottons

Wild relatives of domesticated plants provide a rich resource for crop improvement and a valuable comparative perspective for understanding genomic, physiological, and agricultural traits. Here, we provide high-quality reference genomes of one early domesticated form of the economically most important cotton species, Gossypium hirsutum, and two other wild species, to clarify evolutionary relationships and understand the genomic changes that characterize these species and their close relatives. We document abundant gene resources involved in adaptation to environmental challenges, highlighting the potential for introgression of favorable genes into domesticated cotton and for increasing resilience to climate variability. Our study complements other recent genomic analyses in the cotton genus and provides a valuable foundation for breeding improved cotton varieties.

Allotetraploid cotton (Gossypium) species represents a model system for the study of plant polyploidy, molecular evolution, and domestication. Here, chromosome-scale genome sequences were obtained and assembled for two recently described wild species of tetraploid cotton, Gossypium ekmanianum [(AD)₆, Ge] and Gossypium stephensii [(AD)₇, Gs], and one early form of domesticated Gossypium hirsutum, race punctatum [(AD)₁, Ghp]. Based on phylogenomic analysis, we provide a dated whole-genome level perspective for the evolution of the tetraploid Gossypium clade and resolved the evolutionary relationships of Gs, Ge, and domesticated G. hirsutum. We describe genomic structural variation that arose during Gossypium evolution and describe its correlates—including phenotypic differentiation, genetic isolation, and genetic convergence—that contributed to cotton biodiversity and cotton domestication. Presence/absence variation is prominent in causing cotton genomic structural variations. A presence/absence variation-derived gene encoding a phosphopeptide-binding protein is implicated in increasing fiber length during cotton domestication. The relatively unimproved Ghp offers the potential for gene discovery related to adaptation to environmental challenges. Expanded gene families enoyl-CoA δ isomerase 3 and RAP2-7 may have contributed to abiotic stress tolerance, possibly by targeting plant hormone-associated biochemical pathways. Our results generate a genomic context for a better understanding of cotton evolution and for agriculture.

Genome sequence of Gossypium anomalum facilitates interspecific introgression breeding

Crop wild relatives are an important reservoir of natural biodiversity. However, incorporating wild genetic diversity into breeding programs is often hampered by reproductive barriers and a lack of accurate genomic information. We assembled a high-quality, accurately centromere-anchored genome of Gossypium anomalum, a stress-tolerant wild cotton species. We provided a strategy to discover and transfer agronomically valuable genes from wild diploid species to tetraploid cotton cultivars. With a (Gossypium hirsutum × G. anomalum)2 hexaploid as a bridge parent, we developed a set of 74 diploid chromosome segment substitution lines (CSSLs) of the wild cotton species G. anomalum in the G. hirsutum background. This set of CSSLs included 70 homozygous substitutions and four heterozygous substitutions, and it collectively contained about 72.22% of the G. anomalum genome. Twenty-four quantitative trait loci associated with plant height, yield, and fiber qualities were detected on 15 substitution segments. Integrating the reference genome with agronomic trait evaluation of the CSSLs enabled location and cloning of two G. anomalum genes that encode peroxiredoxin and putative callose synthase 8, respectively, conferring drought tolerance and improving fiber strength. We have demonstrated the power of a high-quality wild-species reference genome for identifying agronomically valuable alleles to facilitate interspecific introgression breeding in crops.

The Conservation of Long Intergenic Non-Coding RNAs and Their Response to Verticillium dahliae Infection in Cotton

Long intergenic non-coding RNAs (lincRNAs) have been demonstrated to be vital regulators of diverse biological processes in both animals and plants. While many lincRNAs have been identified in cotton, we still know little about the repositories and conservativeness of lincRNAs in different cotton species or about their role in responding to biotic stresses. Here, by using publicly available RNA-seq datasets from diverse sources, including experiments of Verticillium dahliae (Vd) infection, we identified 24,425 and 17,713 lincRNAs, respectively, in Gossypium hirsutum (Ghr) and G. barbadense (Gba), the two cultivated allotetraploid cotton species, and 6933 and 5911 lincRNAs, respectively, in G. arboreum (Gar) and G. raimondii (Gra), the two extant diploid progenitors of the allotetraploid cotton. While closely related subgenomes, such as Ghr_At and Gba_At, tend to have more conserved lincRNAs, most lincRNAs are species-specific. The majority of the synthetic and transcribed lincRNAs (78.2%) have a one-to-one orthologous relationship between different (sub)genomes, although a few of them (0.7%) are retained in all (sub)genomes of the four species. The Vd responsiveness of lincRNAs seems to be positively associated with their conservation level. The major functionalities of the Vd-responsive lincRNAs seem to be largely conserved amongst Gra, Ghr, and Gba. Many Vd-responsive Ghr-lincRNAs overlap with Vd-responsive QTL, and several lincRNAs were predicted to be endogenous target mimicries of miR482/2118, with a pair being highly conserved between Ghr and Gba. On top of the confirmation of the feature characteristics of the lincRNAs previously reported in cotton and other species, our study provided new insights into the conservativeness and divergence of lincRNAs during cotton evolution and into the relationship between the conservativeness and Vd responsiveness of lincRNAs. The study also identified candidate lincRNAs with a potential role in disease response for functional characterization.

Conservation and Divergence of Phosphoenolpyruvate Carboxylase Gene Family in Cotton

Phosphoenolpyruvate carboxylase (PEPC) is an important enzyme in plants, which regulates carbon flow through the TCA cycle and controls protein and oil biosynthesis. Although it is important, there is little research on PEPC in cotton, the most important fiber crop in the world. In this study, a total of 125 PEPCs were identified in 15 Gossypium genomes. All PEPC genes in cotton are divided into six groups and each group generally contains one PEPC member in each diploid cotton and two in each tetraploid cotton. This suggests that PEPC genes already existed in cotton before their divergence. There are additional PEPC sub-groups in other plant species, suggesting the different evolution and natural selection during different plant evolution. PEPC genes were independently evolved in each cotton sub-genome. During cotton domestication and evolution, certain PEPC genes were lost and new ones were born to face the new environmental changes and human being needs. The comprehensive analysis of collinearity events and selection pressure shows that genome-wide duplication and fragment duplication are the main methods for the expansion of the PEPC family, and they continue to undergo purification selection during the evolutionary process. PEPC genes were widely expressed with temporal and spatial patterns. The expression patterns of PEPC genes were similar in G. hirsutum and G. barbadense with a slight difference. PEPC2A and 2D were highly expressed in cotton reproductive tissues, including ovule and fiber at all tested developmental stages in both cultivated cottons. However, PEPC1A and 1D were dominantly expressed in vegetative tissues. Abiotic stress also induced the aberrant expression of PEPC genes, in which PEPC1 was induced by both chilling and salinity stresses while PEPC5 was induced by chilling and drought stresses. Each pair (A and D) of PEPC genes showed the similar response to cotton development and different abiotic stress, suggesting the similar function of these PEPCs no matter their origination from A or D sub-genome. However, some divergence was also observed among their origination, such as PEPC5D was induced but PEPC5A was inhibited in G. barbadense during drought treatment, suggesting that a different organized PEPC gene may evolve different functions during cotton evolution. During cotton polyploidization, the homologues genes may refunction and play different roles in different situations.

Comprehensive identification and expression analysis of CRY gene family in Gossypium

BACKGROUND

The cryptochromes (CRY) are specific blue light receptors of plants and animals, which play crucial roles in physiological processes of plant growth, development, and stress tolerance.

RESULTS

In the present work, a systematic analysis of the CRY gene family was performed on twelve cotton species, resulting in 18, 17, 17, 17, and 17 CRYs identified in five alloteraploid cottons (Gossypium hirsutum, G. barbadense, G. tomentosum, G. mustelinum and G. darwinii), respectively, and five to nine CRY genes in the seven diploid species. Phylogenetic analysis of protein-coding sequences revealed that CRY genes from cottons and Arabidopsis thaliana could be classified into seven clades. Synteny analysis suggested that the homoeolog of G. hirsutum Gh_A02G0384 has undergone an evolutionary loss event in the other four allotetraploid cotton species. Cis-element analysis predicated the possible functions of CRY genes in G. hirsutum. RNA-seq data revealed that Gh_D09G2225, Gh_A09G2012 and Gh_A11G1040 had high expressions in fiber cells of different developmental states. In addition, the expression levels of one (Gh_A03G0120), 15 and nine GhCRY genes were down-regulated following the PEG, NaCl and high-temperature treatments, respectively. For the low-temperature treatment, five GhCRY genes were induced, and five were repressed. These results indicated that most GhCRY genes negatively regulate the abiotic stress treatments.

CONCLUSION

We report the structures, domains, divergence, synteny, and cis-elements analyses systematically of G. hirsutum CRY genes. Possible biological functions of GhCRY genes in differential tissues as well as in response to abiotic stress during the cotton plant life cycle were predicted.

GRAND: An Integrated Genome, Transcriptome Resources, and Gene Network Database for Gossypium

With the increasing amount of cotton omics data, breeding scientists are confronted with the question of how to use massive cotton data to mine effective breeding information. Here, we construct a Gossypium Resource And Network Database (GRAND), which integrates 18 cotton genome sequences, genome annotations, two cotton genome variations information, and also four transcriptomes for Gossypium species. GRAND allows to explore and mine this data with the help of a toolbox that comprises a flexible search system, BLAST and BLAT suite, orthologous gene ID, networks of co-expressed genes, primer design, Gbrowse and Jbrowse, and drawing instruments. GRAND provides important information regarding Gossypium resources and hopefully can accelerate the progress of cultivating cotton varieties.

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.

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.

The use of ribosome-nascent chain complex-seq to reveal the translated mRNA profile and the role of ASN1 in resistance to Verticillium wilt in cotton

We combined traditional mRNA-seq and RNC-seq together to reveal post-transcriptional regulation events impacting gene expression and interactions between the serious fungal pathogen Verticillium dahliae and a susceptible host, Gossypium hirsutum TM-1. After screening the differentially expressed and translated genes, V. dahliae infection was observed to influence gene transcription and translation in its host. Interestingly, the asparagine synthase (ASN1) gene transcripts increased significantly with the increase of infection time, while the rate of ASN1 protein accumulation in host TM-1 was distinctly lower than that in resistant hosts. We knocked down the ASN1 gene in resistant plants (ZZM2), and found that Verticillium-resistance was significantly reduced upon knockdown of ASN1. Our study revealed both transcriptional and post-transcriptional regulation of gene expression in TM-1 cotton plants infected by V. dahliae, and showed that ASN1 functions in the V. dahliae resistance process. These insights support breeding of disease resistance in cotton.

Comparative Genome Analyses Highlight Transposon-Mediated Genome Expansion and the Evolutionary Architecture of 3D Genomic Folding in Cotton

Transposable element (TE) amplification has been recognized as a driving force mediating genome size expansion and evolution, but the consequences for shaping 3D genomic architecture remains largely unknown in plants. Here, we report reference-grade genome assemblies for three species of cotton ranging 3-fold in genome size, namely Gossypium rotundifolium (K2), G. arboreum (A2), and G. raimondii (D5), using Oxford Nanopore Technologies. Comparative genome analyses document the details of lineage-specific TE amplification contributing to the large genome size differences (K2, 2.44 Gb; A2, 1.62 Gb; D5, 750.19 Mb) and indicate relatively conserved gene content and synteny relationships among genomes. We found that approximately 17% of syntenic genes exhibit chromatin status change between active ("A") and inactive ("B") compartments, and TE amplification was associated with the increase of the proportion of A compartment in gene regions (∼7,000 genes) in K2 and A2 relative to D5. Only 42% of topologically associating domain (TAD) boundaries were conserved among the three genomes. Our data implicate recent amplification of TEs following the formation of lineage-specific TAD boundaries. This study sheds light on the role of transposon-mediated genome expansion in the evolution of higher-order chromatin structure in plants.

Recent Advances and Future Perspectives in Cotton Research

Cotton is not only the world's most important natural fiber crop, but it is also an ideal system in which to study genome evolution, polyploidization, and cell elongation. With the assembly of five different cotton genomes, a cotton-specific whole-genome duplication with an allopolyploidization process that combined the A- and D-genomes became evident. All existing A-genomes seemed to originate from the A₀-genome as a common ancestor, and several transposable element bursts contributed to A-genome size expansion and speciation. The ethylene production pathway is shown to regulate fiber elongation. A tip-biased diffuse growth mode and several regulatory mechanisms, including plant hormones, transcription factors, and epigenetic modifications, are involved in fiber development. Finally, we describe the involvement of the gossypol biosynthetic pathway in the manipulation of herbivorous insects, the role of GoPGF in gland formation, and host-induced gene silencing for pest and disease control. These new genes, modules, and pathways will accelerate the genetic improvement of cotton.

Genome-Wide Comparative Analysis of Heat Shock Transcription Factors Provides Novel Insights for Evolutionary History and Expression Characterization in Cotton Diploid and Tetraploid Genomes

Heat shock transcription factors (HSFs) are involved in environmental stress response and plant development, such as heat stress and flowering development. According to the structural characteristics of the HSF gene family, HSF genes were classified into three major types (HSFA, HSFB, and HSFC) in plants. Using conserved domains of HSF genes, we identified 621 HSF genes among 13 cotton genomes, consisting of eight diploid and five tetraploid genomes. Phylogenetic analysis indicated that HSF genes among 13 cotton genomes were grouped into two different clusters: one cluster contained all HSF genes of HSFA and HSFC, and the other cluster contained all HSF genes of HSFB. Comparative analysis of HSF genes in Arabidopsis thaliana, Gossypium herbaceum (A1), Gossypium arboreum (A2), Gossypium raimondii (D5), and Gossypium hirsutum (AD1) genomes demonstrated that four HSF genes were inherited from a common ancestor, A0, of all existing cotton A genomes. Members of the HSF gene family in G. herbaceum (A1) genome indicated a significant loss compared with those in G. arboretum (A2) and G. hirsutum (AD1) A genomes. However, HSF genes in G. raimondii (D5) showed relative loss compared with those in G. hirsutum (AD1) D genome. Analysis of tandem duplication (TD) events of HSF genes revealed that protein-coding genes among different cotton genomes have experienced TD events, but only the two-gene tandem array was detected in Gossypium thurberi (D1) genome. The expression analysis of HSF genes in G. hirsutum (AD1) and Gossypium barbadense (AD2) genomes indicated that the expressed HSF genes were divided into two different groups, respectively, and the expressed HSF orthologous genes between the two genomes showed totally different expression patterns despite the implementation of the same abiotic stresses. This work will provide novel insights for the study of evolutionary history and expression characterization of HSF genes in different cotton genomes and a widespread application model for the study of HSF gene families in plants.

Genome-Wide Analysis of DA1-Like Genes in Gossypium and Functional Characterization of GhDA1-1A Controlling Seed Size

Cotton (Gossypium spp.) is an economically important crop grown for natural fiber and seed oil production. DA1 is a ubiquitin receptor that determines final seed and organ size by restricting the period of cell proliferation. In the present study, we identified 7 DA1-like genes each in cultivated tetraploid (AADD) G. hirsutum and G. barbadense, and 4 and 3 DA1-like genes in their ancestral diploid G. arboreum (A2A2) and G. raimondii (D5D5), respectively. The 7 GhDA1 genes were confirmed to be distributed on four At and three Dt subgenome chromosomes in G. hirsutum. GhDA1-1A showed a high sequence similarity to AtDA1 in Arabidopsis, and they possessed the same functional domains, suggesting conserved functions. The overexpression of GhDA1-1A ^R301K in Arabidopsis significantly increased seed size and seed weight, indicating that GhDA1-1A is a promising target for cotton improvement. This study provides information on the molecular evolutionary properties of DA1-like genes in cotton, which will be useful for the genetic improvement of cotton.

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

Transposable element (TE) amplification has been recognized as a driving force mediating genome size expansion and evolution, but the consequences for shaping 3D genomic architecture remains largely unknown in plants. Here, we report reference-grade genome assemblies for three species of cotton ranging 3-fold in genome size, namely Gossypium rotundifolium (K₂), G. arboreum (A₂), and G. raimondii (D₅), using Oxford Nanopore Technologies. Comparative genome analyses document the details of lineage-specific TE amplification contributing to the large genome size differences (K₂, 2.44 Gb; A₂, 1.62 Gb; D₅, 750.19 Mb) and indicate relatively conserved gene content and synteny relationships among genomes. We found that approximately 17% of syntenic genes exhibit chromatin status change between active (“A”) and inactive (“B”) compartments, and TE amplification was associated with the increase of the proportion of A compartment in gene regions (∼7,000 genes) in K₂ and A₂ relative to D₅. Only 42% of topologically associating domain (TAD) boundaries were conserved among the three genomes. Our data implicate recent amplification of TEs following the formation of lineage-specific TAD boundaries. This study sheds light on the role of transposon-mediated genome expansion in the evolution of higher-order chromatin structure in plants.

Gene Expression Correlation Analysis Reveals MYC-NAC Regulatory Network in Cotton Pigment Gland Development

Plant NAC (NAM, ATAF1/2, and CUC2) family is involved in various development processes including Programmed Cell Death (PCD) associated development. However, the relationship between NAC family and PCD-associated cotton pigment gland development is largely unknown. In this study, we identified 150, 153 and 299 NAC genes in newly updated genome sequences of G. arboreum, G. raimondii and G. hirsutum, respectively. All NAC genes were divided into 8 groups by the phylogenetic analysis and most of them were conserved during cotton evolution. Using the vital regulator of gland formation GhMYC2-like as bait, expression correlation analysis screened out 6 NAC genes which were low-expressed in glandless cotton and high-expressed in glanded cotton. These 6 NAC genes acted downstream of GhMYC2-like and were induced by MeJA. Silencing CGF1(Cotton Gland Formation1), another MYC-coding gene, caused almost glandless phenotype and down-regulated expression of GhMYC2-like and the 6 NAC genes, indicating a MYC-NAC regulatory network in gland development. In addition, predicted regulatory mechanism showed that the 6 NAC genes were possibly regulated by light, various phytohormones and transcription factors as well as miRNAs. The interaction network and DNA binding sites of the 6 NAC transcription factors were also predicted. These results laid the foundation for further study of gland-related genes and gland development regulatory network.

The Gossypium stocksii genome as a novel resource for cotton improvement

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 (Gossypium 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.

Molecular cloning and characterization of GhERF105, a gene contributing to the regulation of gland formation in upland cotton (Gossypium hirsutum L.)

BACKGROUND

Gossypium hirsutum L. (cotton) is one of the most economically important crops in the world due to its significant source of fiber, feed, foodstuff, oil and biofuel products. However, the utilization of cottonseed was limited due to the presence of small and darkly pigmented glands that contain large amounts of gossypol, which is toxic to human beings and non-ruminant animals. To date, some progress has been made in the pigment gland formation, but the underlying molecular mechanism of its formation was still unclear.

RESULTS

In this study, we identified an AP2/ERF transcription factor named GhERF105 (GH_A12G2166), which was involved in the regulation of gland pigmentation by the comparative transcriptome analysis of the leaf of glanded and glandless plants. It encoded an ERF protein containing a converved AP2 domain which was localized in the nucleus with transcriptional activity, and showed the high expression in glanded cotton accessions that contained much gossypol. Virus-induced gene silencing (VIGS) against GhERF105 caused the dramatic reduction in the number of glands and significantly lowered levels of gossypol in cotton leaves. GhERF105 showed the patterns of spatiotemporal and inducible expression in the glanded plants.

CONCLUSIONS

These results suggest that GhERF105 contributes to the pigment gland formation and gossypol biosynthesis in partial organs of glanded plant. It also provides a potential molecular basis to generate 'glandless-seed' and 'glanded-plant' cotton cultivar.

Transcriptomic and proteomic analyses of a new cytoplasmic male sterile line with a wild Gossypium bickii genetic background

Background

Cotton is an important fiber crop but has serious heterosis effects, and cytoplasmic male sterility (CMS) is the major cause of heterosis in plants. However, to the best of our knowledge,

no studies have investigated CMS Yamian A in cotton with the genetic background of Australian wild Gossypium bickii. Conjoint transcriptomic and proteomic analysis was first performed between Yamian A and its maintainer Yamian B.

Results

We detected 550 differentially expressed transcript-derived fragments (TDFs) and at least 1013 proteins in anthers at various developmental stages. Forty-two TDFs and 11 differentially expressed proteins (DEPs) were annotated by analysis in the genomic databases of G. austral, G. arboreum and G. hirsutum. Gene Ontology and Kyoto Encyclopedia of Genes and Genomes pathway analyses were performed to better understand the functions of these TDFs and DEPs. Transcriptomic and proteomic results showed that UDP-glucuronosyl/UDP-glucosyltransferase, 60S ribosomal protein L13a-4-like, and glutathione S-transferase were upregulated; while heat shock protein Hsp20, ATPase, F0 complex, and subunit D were downregulated at the microspore abortion stage of Yamian A. In addition, several TDFs from the transcriptome and several DEPs from the proteome were detected and confirmed by quantitative real-time PCR as being expressed in the buds of seven different periods of development. We established the databases of differentially expressed genes and proteins between Yamian A and its maintainer Yamian B in the anthers at various developmental stages and constructed an interaction network based on the databases for a comprehensive understanding of the mechanism underlying CMS with a wild cotton genetic background.

Conclusion

We first analyzed the molecular mechanism of CMS Yamian A from the perspective of omics, thereby providing an experimental basis and theoretical foundation for future research attempting to analyze the abortion mechanism of new CMS with a wild Gossypium bickii background and to realize three-line matching.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-020-07261-y.

Genomics and lipidomics analysis of the biotechnologically important oleaginous red yeast Rhodotorula glutinis ZHK provides new insights into its lipid and carotenoid metabolism

Background

Rhodotorula glutinis is recognized as a biotechnologically important oleaginous red yeast, which synthesizes numerous meritorious compounds with wide industrial usages. One of the most notable properties of R. glutinis is the formation of intracellular lipid droplets full of carotenoids. However, the basic genomic features that underlie the biosynthesis of these valuable compounds in R. glutinis have not been fully documented. To reveal the biotechnological potential of R. glutinis, the genomics and lipidomics analysis was performed through the Next-Generation Sequencing and HPLC-MS-based metabolomics technologies.

Results

Here, we firstly assemble the genome of R. glutinis ZHK into 21.8 Mb, containing 30 scaffolds and 6774 predicted genes with a N50 length of 14, 66,672 bp and GC content of 67.8%. Genome completeness assessment (BUSCO alignment: 95.3%) indicated the genome assembly with a high-quality features. According to the functional annotation of the genome, we predicted several key genes involved in lipids and carotenoids metabolism as well as certain industrial enzymes biosynthesis. Comparative genomics results suggested that most of orthologous genes have underwent the strong purifying selection within the five Rhodotorula species, especially genes responsible for carotenoids biosynthesis. Furthermore, a total of 982 lipids were identified using the lipidomics approaches, mainly including triacylglycerols, diacylglyceryltrimethylhomo-ser and phosphatidylethanolamine.

Conclusion

Using whole genome shotgun sequencing, we comprehensively analyzed the genome of R. glutinis and predicted several key genes involved in lipids and carotenoids metabolism. By performing comparative genomic analysis, we show that most of the ortholog genes have undergone strong purifying selection within the five Rhodotorula species. Furthermore, we identified 982 lipid species using lipidomic approaches. These results provided valuable resources to further advance biotechnological applications of R .glutinis.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-020-07244-z.

Transcriptome analysis reveals underlying immune response mechanism of fungal (Penicillium oxalicum) disease in Gastrodia elata Bl. f. glauca S. chow (Orchidaceae)

BACKGROUND

Gastrodia elata Bl. f. glauca S. Chow is a medicinal plant. G. elata f. glauca is unavoidably infected by pathogens in their growth process. In previous work, we have successfully isolated and identified Penicillium oxalicum from fungal diseased tubers of G. elata f. glauca. As a widespread epidemic, this fungal disease seriously affected the yield and quality of G. elata f. glauca. We speculate that the healthy G. elata F. glauca might carry resistance genes, which can resist against fungal disease. In this study, healthy and fungal diseased mature tubers of G. elata f. glauca from Changbai Mountain area were used as experimental materials to help us find potential resistance genes against the fungal disease.

RESULTS

A total of 7540 differentially expressed Unigenes (DEGs) were identified (FDR < 0.01, log2FC > 2). The current study screened 10 potential resistance genes. They were attached to transcription factors (TFs) in plant hormone signal transduction pathway and plant pathogen interaction pathway, including WRKY22, GH3, TIFY/JAZ, ERF1, WRKY33, TGA. In addition, four of these genes were closely related to jasmonic acid signaling pathway.

CONCLUSIONS

The immune response mechanism of fungal disease in G. elata f. glauca is a complex biological process, involving plant hormones such as ethylene, jasmonic acid, salicylic acid and disease-resistant transcription factors such as WRKY, TGA.

Comprehensive Analysis of Respiratory Burst Oxidase Homologs (Rboh) Gene Family and Function of GbRboh5/18 on Verticillium Wilt Resistance in Gossypium barbadense

Respiratory burst oxidase homologs (Rbohs) play a predominant role in reactive oxygen species (ROS) production, which is crucial in plant growth, differentiation, as well as their responses to biotic and abiotic stresses. To date, however, there is little knowledge about the function of cotton Rboh genes. Here, we identified a total of 87 Rbohs from five sequenced Gossypium species (the diploids Gossypium arboreum, Gossypium raimondii, and Gossypium australe, and the allotetraploids Gossypium hirsutum and Gossypium barbadense) via BLAST searching their genomes. Phylogenetic analysis of the putative 87 cotton Rbohs revealed that they were divided into seven clades. All members within the same clade are generally similar to each other in terms of gene structure and conserved domain arrangement. In G. barbadense, the expression levels of GbRbohs in the CladeD were induced in response to a fungal pathogen and to hormones (i.e., jasmonic acid and abscisic acid), based upon which the main functional member in CladeD was discerned to be GbRboh5/18. Further functional and physiological analyses showed that the knock-down of GbRboh5/18 expression attenuates plant resistance to Verticillium dahliae infection. Combined with the molecular and biochemical analyses, we found less ROS accumulation in GbRboh5/18-VIGS plants than in control plants after V. dahliae infection. Overexpression of GbRboh5/18 in G. barbadense resulted in more ROS accumulation than in control plants. These results suggest that GbRboh5/18 enhances the cotton plants' resistance against V. dahliae by elevating the levels of ROS accumulation. By integrating phylogenetic, molecular, and biochemical approaches, this comprehensive study provides a detailed overview of the number, phylogeny, and evolution of the Rboh gene family from five sequenced Gossypium species, as well as elucidating the function of GbRboh5/18 for plant resistance against V. dahliae. This study sheds fresh light on the molecular evolutionary properties and function of Rboh genes in cotton, and provides a reference for improving cotton's responses to the pathogen V. dahliae.

The genome of kenaf (Hibiscus cannabinus L.) provides insights into bast fibre and leaf shape biogenesis

Kenaf is an annual crop that is widely cultivated as a source of bast (phloem) fibres, the phytoremediation of heavy metal-contaminated farmlands and textile-relevant compounds. Leaf shape played a unique role in kenaf improvement, due to the inheritance as a single locus and the association with fibre development in typical lobed-leaf varieties. Here we report a high-quality genome assembly and annotation for var. 'Fuhong 952' with 1078 Mbp genome and 66 004 protein-coding genes integrating single-molecule real-time sequencing, a high-density genetic map and high-throughput chromosome conformation capture techniques. Gene mapping assists the identification of a homeobox transcription factor LATE MERISTEM IDENTITY 1 (HcLMI1) gene controlling lobed-leaf. Virus-induced gene silencing (VIGS) of HcLMI1 in a lobed-leaf variety was critical to induce round (entire)-like leaf formation. Candidate genes involved in cell wall formation were found in quantitative trait loci (QTL) for fibre yield and quality-related traits. Comparative genomic and transcriptome analyses revealed key genes involved in bast fibre formation, among which there are twice as many cellulose synthase A (CesA) genes due to a recent whole-genome duplication after divergence from Gossypium. Population genomic analysis showed two recent population bottlenecks in kenaf, suggesting domestication and improvement process have led to an increase in fibre biogenesis and yield. This chromosome-scale genome provides an important framework and toolkit for sequence-directed genetic improvement of fibre crops.

Genome-wide analysis and comparison of the DNA-binding one zinc finger gene family in diploid and tetraploid cotton (Gossypium)

The Dof (DNA-binding one zinc finger) transcription factor family is a representative of plant-specific classes of transcription factors. In this study, we performed a genome-wide screening and characterization of the Dof gene family within two tetraploid species Gossypium barbadense, Gossypium hirsutum, and two diploid species Gossypium arboreum, Gossypium raimondii. 115, 116, 55 and 56 Dof genes were identified respectively and all of the genes contain a sequence encoding the Dof DNA-binding domain. Those genes were unevenly distributed across 13/26 chromosomes of the cotton. Genome comparison revealed that segmental duplication may have played crucial roles in the expansion of the cotton Dof gene family, and tandem duplication also played a minor role. Analysis of RNA-Seq data indicated that cotton Dof gene expression levels varied across different tissues and in response to different abiotic stress. Overall, our results could provide valuable information for better understanding the evolution of cotton Dof genes, and lays a foundation for future investigation in cotton.

Comparative Transcriptomic Analysis to Identify the Genes Related to Delayed Gland Morphogenesis in Gossypium bickii

Cotton is one of the major industrial crops that supply natural fibers and oil for industries. This study was conducted to understand the mechanism of delayed gland morphogenesis in seeds of Gossypium bickii. In this study, we compared glandless seeds of G. bickii with glanded seeds of Gossypium arboreum. High-throughput sequencing technology was used to explore and classify the expression patterns of gland-related genes in seeds and seedlings of cotton plants. Approximately 131.33 Gigabases of raw data from 12 RNA sequencing samples with three biological replicates were generated. A total of 7196 differentially-expressed genes (DEGs) were identified in all transcriptome data. Among them, 3396 genes were found up-regulated and 3480 genes were down-regulated. Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) annotations were performed to identify different functions between genes unique to glandless imbibed seeds and glanded seedlings. Co-expression network analysis revealed four modules that were identified as highly associated with the development of glandless seeds. Here the hub genes in each module were identified by weighted gene co-expression network analysis (WGCNA). In total, we have selected 13 genes involved in transcription factors, protein and MYB-related functions, that were differentially expressed in transcriptomic data and validated by quantitative reverse-transcription polymerase chain reaction (qRT-PCR). These selected genes may play an important role for delayed gland morphogenesis. Our study provides comprehensive insight into the key genes related to glandless traits of seeds and plants, and can be further exploited by functional and molecular studies.

Differential Expression of Circular RNAs in Polytocous and Monotocous Uterus during the Reproductive Cycle of Sheep

Simple Summary

The uterus is an important reproductive organ that provides nutrition and place for embryonic development. In this study, we identified circular RNAs by deep sequencing and analyzed their expression in the uteri of polytocous and monotocous sheep (FecB++) during follicular and luteal phases. Gene Ontology (GO) and KEGG enrichment analyses revealed that the source genes of these differential circular RNAs (circRNAs) were mainly enriched in reproductive hormone- and energy metabolism-related pathways. These results provide information on the molecular mechanisms of sheep prolificacy.

Abstract

CircRNA plays important roles in cell proliferation, differentiation, autophagy and apoptosis during development. However, there are few reports on circRNAs related to livestock reproduction. In this study, we identified circRNAs by deep sequencing and analyzed their expression in the uteri of polytocous and monotocous sheep (FecB++) during follicular and luteal phases. There were 147 and 364 circRNAs with differential expression in the follicular and luteal phases, respectively. GO and KEGG enrichment analysis was performed for the host genes of the circRNAs to predict the functions of differentially expressed circRNAs. These source genes were mainly involved in the estrogen signaling pathway, TGFβ signaling pathway, GnRH signaling pathway, oxytocin signaling pathway, pentose phosphate pathway, and starch and sucrose metabolism related to reproduction and energy metabolism. CircRNA expression patterns were validated by RT-qPCR. Our findings provide a solid foundation for the identification and characterization of key important circRNAs involved in reproduction.