GhRCD1 promotes cotton tolerance to cadmium by regulating the GhbHLH12-GhMYB44-GhHMA1 transcriptional cascade.

Heavy metal pollution poses a significant risk to human health and wreaks havoc on agricultural productivity. Phytoremediation, a plant-based, environmentally benign, and cost-effective method, is employed to remove heavy metals from contaminated soil, particularly in agricultural or heavy metal-sensitive lands. However, the phytoremediation capacity of various plant species and germplasm resources display significant genetic diversity, and the mechanisms underlying these differences remain hitherto obscure. Given its potential benefits, genetic improvement of plants is essential for enhancing their uptake of heavy metals, tolerance to harmful levels, as well as overall growth and development in contaminated soil. In this study, we uncover a molecular cascade that regulates cadmium (Cd2+ ) tolerance in cotton, involving GhRCD1, GhbHLH12, GhMYB44, and GhHMA1. We identified a Cd2+ -sensitive cotton T-DNA insertion mutant with disrupted GhRCD1 expression. Genetic knockout of GhRCD1 by CRISPR/Cas9 technology resulted in reduced Cd2+ tolerance in cotton seedlings, while GhRCD1 overexpression enhanced Cd2+ tolerance. Through molecular interaction studies, we demonstrated that, in response to Cd2+ presence, GhRCD1 directly interacts with GhbHLH12. This interaction activates GhMYB44, which subsequently activates a heavy metal transporter, GhHMA1, by directly binding to a G-box cis-element in its promoter. These findings provide critical insights into a novel GhRCD1-GhbHLH12-GhMYB44-GhHMA1 regulatory module responsible for Cd2+ tolerance in cotton. Furthermore, our study paves the way for the development of elite Cd2+ -tolerant cultivars by elucidating the molecular mechanisms governing the genetic control of Cd2+ tolerance in cotton.

GhRCD1 regulates cotton somatic embryogenesis by modulating the GhMYC3-GhMYB44-GhLBD18 transcriptional cascade.

Plant somatic embryogenesis (SE) is a multifactorial developmental process where embryos that can develop into whole plants are produced from somatic cells rather than through the fusion of gametes. The molecular regulation of plant SE, which involves the fate transition of somatic cells into embryogenic cells, is intriguing yet remains elusive. We deciphered the molecular mechanisms by which GhRCD1 interacts with GhMYC3 to regulate cell fate transitions during SE in cotton. While silencing of GhMYC3 had no discernible effect on SE, its overexpression accelerated callus formation, and proliferation. We identified two of GhMYC3 downstream SE regulators, GhMYB44 and GhLBD18. GhMYB44 overexpression was unconducive to callus growth but bolstered EC differentiation. However, GhLBD18 can be triggered by GhMYC3 but inhibited by GhMYB44, which positively regulates callus growth. On top of the regulatory cascade, GhRCD1 antagonistically interacts with GhMYC3 to inhibit the transcriptional function of GhMYC3 on GhMYB44 and GhLBD18, whereby a CRISPR-mediated rcd1 mutation expedites cell fate transition, resembling the effects of GhMYC3 overexpression. Furthermore, we showed that reactive oxygen species (ROS) are involved in SE regulation. Our findings elucidated that SE homeostasis is maintained by the tetrapartite module, GhRCD1-GhMYC3-GhMYB44-GhLBD18, which acts to modulate intracellular ROS in a temporal manner.

High-quality genome assembly of Verticillium dahliae VD991 allows for screening and validation of pathogenic genes.

Verticillium dahliae (V. dahliae) is a notorious soil-borne pathogen causing Verticillium wilt in more than 400 dicotyledonous plants, including a wide range of economically important crops, such as cotton, tomato, lettuce, potato, and romaine lettuce, which can result in extensive economic losses. In the last decade, several studies have been conducted on the physiological and molecular mechanisms of plant resistance to V. dahliae. However, the lack of a complete genome sequence with a high-quality assembly and complete genomic annotations for V. dahliae has limited these studies. In this study, we produced a full genomic assembly for V. dahliae VD991 using Nanopore sequencing technology, consisting of 35.77 Mb across eight pseudochromosomes and with a GC content of 53.41%. Analysis of the genome completeness assessment (BUSCO alignment: 98.62%; Illumina reads alignment: 99.17%) indicated that our efforts resulted in a nearly complete and high-quality genomic assembly. We selected 25 species closely related to V. dahliae for evolutionary analysis, confirming the evolutionary relationship between V. dahliae and related species, and the identification of a possible whole genome duplication event in V. dahliae. The interaction between cotton and V. dahliae was investigated by transcriptome sequencing resulting in the identification of many genes and pathways associated with cotton disease resistance and V. dahliae pathogenesis. These results will provide new insights into the pathogenic mechanisms of V. dahliae and contribute to the cultivation of cotton varieties resistant to Verticillium wilt.

Interactions between Verticillium dahliae and cotton: pathogenic mechanism and cotton resistance mechanism to Verticillium wilt.

Cotton is widely grown in many countries around the world due to the huge economic value of the total natural fiber. Verticillium wilt, caused by the soil-borne pathogen Verticillium dahliae, is the most devastating disease that led to extensive yield losses and fiber quality reduction in cotton crops. Developing resistant cotton varieties through genetic engineering is an effective, economical, and durable strategy to control Verticillium wilt. However, there are few resistance gene resources in the currently planted cotton varieties, which has brought great challenges and difficulties for breeding through genetic engineering. Further revealing the molecular mechanism between V. dahliae and cotton interaction is crucial to discovering genes related to disease resistance. In this review, we elaborated on the pathogenic mechanism of V. dahliae and the resistance mechanism of cotton to Verticillium wilt. V. dahliae has evolved complex mechanisms to achieve pathogenicity in cotton, mainly including five aspects: (1) germination and growth of microsclerotia; (2) infection and successful colonization; (3) adaptation to the nutrient-deficient environment and competition of nutrients; (4) suppression and manipulation of cotton immune responses; (5) rapid reproduction and secretion of toxins. Cotton has evolved multiple physiological and biochemical responses to cope with V. dahliae infection, including modification of tissue structures, accumulation of antifungal substances, homeostasis of reactive oxygen species (ROS), induction of Ca2+ signaling, the mitogen-activated protein kinase (MAPK) cascades, hormone signaling, and PAMPs/effectors-triggered immune response (PTI/ETI). This review will provide an important reference for the breeding of new cotton germplasm resistant to Verticillium wilt through genetic engineering.

N6 -Methyladenosine mRNA modification regulates transcripts stability associated with cotton fiber elongation.

N6 -Methyladenosine (m6 A) is the most abundant methylation modification in eukaryotic mRNA. The discovery of the dynamic and reversible regulatory mechanism of m6 A has greatly promoted the development of m6 A-led epitranscriptomics. However, the characterization of m6 A in cotton fiber is still unknown. Here, we reveal the potential link between m6 A modification and cotton fiber elongation by parallel m6 A-immunoprecipitation-sequencing (m6 A-seq) and RNA-seq analysis of fibers from the short fiber mutants Ligonliness-2 (Li2 ) and wild-type (WT). This study demonstrated a higher level of m6 A in the Li2 mutant, with the enrichment of m6 A modifications in the stop codon, 3'-untranslated region and coding sequence regions than in WT cotton. In the correlation analysis between genes containing differential m6 A modifications and differentially expressed genes, we identified several genes that could potentially regulate fiber elongation, including cytoskeleton, microtubule binding, cell wall and transcription factors (TFs). We further confirmed that the methylation of m6 A affected the mRNA stability of these fiber elongation-related genes including the TF GhMYB44, which showed the highest expression level in the RNA-seq data and m6 A methylation in the m6 A-seq data. Next, the overexpression of GhMYB44 reduces fiber elongation, whereas the silencing of GhMYB44 produces longer fibers. In summary, these results uncover that m6 A methylation regulated the expression of genes related to fiber development by affecting mRNA's stability, ultimately affecting cotton fiber elongation.

Genome-Wide Investigation and Co-Expression Network Analysis of SBT Family Gene in Gossypium.

Subtilases (SBTs), which belong to the serine peptidases, control plant development by regulating cell wall properties and the activity of extracellular signaling molecules, and affect all stages of the life cycle, such as seed development and germination, and responses to biotic and abiotic environments. In this study, 146 Gossypium hirsutum, 138 Gossypium barbadense, 89 Gossypium arboreum and 84 Gossypium raimondii SBTs were identified and divided into six subfamilies. Cotton SBTs are unevenly distributed on chromosomes. Synteny analysis showed that the members of SBT1 and SBT4 were expanded in cotton compared to Arabidopsis thaliana. Co-expression network analysis showed that six Gossypium arboreum SBT gene family members were in a network, among which five SBT1 genes and their Gossypium hirsutum and Arabidopsis thaliana direct homologues were down-regulated by salt treatment, indicating that the co-expression network might share conserved functions. Through co-expression network and annotation analysis, these SBTs may be involved in the biological processes of auxin transport, ABA signal transduction, cell wall repair and root tissue development. In summary, this study provides valuable information for the study of SBT genes in cotton and excavates SBT genes in response to salt stress, which provides ideas for cotton breeding for salinity resistance.

Integration of eQTL Analysis and GWAS Highlights Regulation Networks in Cotton under Stress Condition.

The genus Gossypium is one of the most economically important crops in the world. Here, we used RNA-seq to quantify gene expression in a collection of G. arboreum seedlings and performed eGWAS on 28,382 expressed genes. We identified a total of 30,089 eQTLs in 10,485 genes, of which >90% were trans-regulate target genes. Using luciferase assays, we confirmed that different cis-eQTL haplotypes could affect promoter activity. We found ~6600 genes associated with ~1300 eQTL hotspots. Moreover, hotspot 309 regulates the expression of 325 genes with roles in stem length, fresh weight, seed germination rate, and genes related to cell wall biosynthesis and salt stress. Transcriptome-wide association study (TWAS) identified 19 candidate genes associated with the cotton growth and salt stress response. The variation in gene expression across the population played an essential role in population differentiation. Only a small number of the differentially expressed genes between South China, the Yangtze River region, and the Yellow River region sites were located in different chromosomal regions. The eQTLs found across the duplicated gene pairs showed conservative cis- or trans- regulation and that the expression levels of gene pairs were correlated. This study provides new insights into the evolution of gene expression regulation in cotton, and identifies eQTLs in stress-related genes for use in breeding improved cotton varieties.

Verticillium dahliae secreted protein Vd424Y is required for full virulence, targets the nucleus of plant cells, and induces cell death.

Fungal pathogens secrete effector proteins that regulate host immunity and can suppress basal defence mechanisms against colonization in plants. Verticillium dahliae is a widespread and destructive soilborne fungus that can cause vascular wilt disease and reduces plant yields. However, little is currently known about how the effectors secreted by V. dahliae function. In this study, we analysed and identified 34 candidate effectors in the V. dahliae secretome and found that Vd424Y, a glycoside hydrolase family 11 protein, was highly upregulated during the early stages of V. dahliae infection in cotton plants. This protein was located in the nucleus and its deletion compromised the virulence of the fungus. The transient expression of Vd424Y in Nicotiana benthamiana induced BAK1- and SOBIR1-dependent cell death and activated both salicylic acid and jasmonic acid signalling. This enhanced its resistance to the oomycetes Phytophthora capsici in a way that depended on its nuclear localization signal and signal peptides. Our results demonstrate that Vd424Y is an important effector protein targeting the host nucleus to regulate and activate effector-triggered immunity 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.

Cotton D genome assemblies built with long-read data unveil mechanisms of centromere evolution and stress tolerance divergence.

BACKGROUND: Many of genome features which could help unravel the often complex post-speciation evolution of closely related species are obscured because of their location in chromosomal regions difficult to accurately characterize using standard genome analysis methods, including centromeres and repeat regions. RESULTS: Here, we analyze the genome evolution and diversification of two recently diverged sister cotton species based on nanopore long-read sequence assemblies and Hi-C 3D genome data. Although D genomes are conserved in gene content, they have diversified in gene order, gene structure, gene family diversification, 3D chromatin structure, long-range regulation, and stress-related traits. Inversions predominate among D genome rearrangements. Our results support roles for 5mC and 6mA in gene activation, and 3D chromatin analysis showed that diversification in proximal-vs-distal regulatory-region interactions shape the regulation of defense-related-gene expression. Using a newly developed method, we accurately positioned cotton centromeres and found that these regions have undergone obviously more rapid evolution relative to chromosome arms. We also discovered a cotton-specific LTR class that clarifies evolutionary trajectories among diverse cotton species and identified genetic networks underlying the Verticillium tolerance of Gossypium thurberi (e.g., SA signaling) and salt-stress tolerance of Gossypium davidsonii (e.g., ethylene biosynthesis). Finally, overexpression of G. thurberi genes in upland cotton demonstrated how wild cottons can be exploited for crop improvement. CONCLUSIONS: Our study substantially deepens understanding about how centromeres have developed and evolutionarily impacted the divergence among closely related cotton species and reveals genes and 3D genome structures which can guide basic investigations and applied efforts to improve crops.

The miR164-GhCUC2-GhBRC1 module regulates plant architecture through abscisic acid in cotton.

Branching determines cotton architecture and production, but the underlying regulatory mechanisms remain unclear. Here, we report that the miR164-GhCUC2 (CUP-SHAPED COTYLEDON2) module regulates lateral shoot development in cotton and Arabidopsis. We generated OE-GhCUC2m (overexpression GhCUC2m) and STTM164 (short tandem target mimic RNA of miR164) lines in cotton and heterologous expression lines for gh-miR164, GhCUC2 and GhCUC2m in Arabidopsis to study the mechanisms controlling lateral branching. GhCUC2m overexpression resulted in a short-branch phenotype similar to STTM164. In addition, heterologous expression of GhCUC2m led to decreased number and length of branches compared with wild type, opposite to the effects of the OE-gh-pre164 line in Arabidopsis. GhCUC2 interacted with GhBRC1 and exhibited similar negative regulation of branching. Overexpression of GhBRC1 in the brc1-2 mutant partially rescued the mutant phenotype and decreased branch number. GhBRC1 directly bound to the NCED1 promoter and activated its transcription, leading to local abscisic acid (ABA) accumulation and response. Mutation of the NCED1 promoter disrupted activation by GhBRC1. This finding demonstrates a direct relationship between BRC1 and ABA signalling and places ABA downstream of BRC1 in the control of branching development. The miR164-GhCUC2-GhBRC1-GhNCED1 module provides a clear regulatory axis for ABA signalling to control plant architecture.

Transcriptome Analysis Revealed GhWOX4 Intercedes Myriad Regulatory Pathways to Modulate Drought Tolerance and Vascular Growth in Cotton.

Cotton is a paramount cash crop around the globe. Among all abiotic stresses, drought is a leading cause of cotton growth and yield loss. However, the molecular link between drought stress and vascular growth and development is relatively uncharted. Here, we validated a crucial role of GhWOX4, a transcription factor, modulating drought stress with that of vasculature growth in cotton. Knock-down of GhWOX4 decreased the stem width and severely compromised vascular growth and drought tolerance. Conversely, ectopic expression of GhWOX4 in Arabidopsis enhanced the tolerance to drought stress. Comparative RNAseq analysis revealed auxin responsive protein (AUX/IAA), abscisic acid (ABA), and ethylene were significantly induced. Additionally, MYC-bHLH, WRKY, MYB, homeodomain, and heat-shock transcription factors (HSF) were differentially expressed in control plants as compared to GhWOX4-silenced plants. The promotor zone of GhWOX4 was found congested with plant growth, light, and stress response related cis-elements. differentially expressed genes (DEGs) related to stress, water deprivation, and desiccation response were repressed in drought treated GhWOX4-virus-induced gene silencing (VIGS) plants as compared to control. Gene ontology (GO) functions related to cell proliferation, light response, fluid transport, and flavonoid biosynthesis were over-induced in TRV: 156-0 h/TRV: 156-1 h (control) in comparison to TRV: VIGS-0 h/TRV: VIGS-1 h (GhWOX4-silenced) plants. This study improves our context for elucidating the pivotal role of GhWOX4 transcription factors (TF), which mediates drought tolerance, plays a decisive role in plant growth and development, and is likely involved in different regulatory pathways in cotton.

GhTIE1 Regulates Branching Through Modulating the Transcriptional Activity of TCPs in Cotton and Arabidopsis.

Transcription factors (TFs) and transcriptional regulators are important switches in transcriptional networks. In recent years, the transcriptional regulator TIE1 (TCP interactor containing EAR motif protein 1) was identified as a nuclear transcriptional repressor which regulates leaf development and controls branch development. However, the function and regulatory network of GhTIE1 has not been studied in cotton. Here, we demonstrated that GhTIE1 is functionally conserved in controlling shoot branching in cotton and Arabidopsis. Overexpression of GhTIE1 in Arabidopsis leads to higher bud vigor and more branches, while silencing GhTIE1 in cotton reduced bud activity and increased branching inhibition. Yeast two-hybrid (Y2H) and bimolecular fluorescence complementation (BiFC) assays showed that GhTIE1 directly interacted with subclass II TCPs (GhBRC1, GhBRC2, and GhTCP13) in vivo and in vitro. Overexpression of GhBRC1, GhBRC2, and GhTCP13 in mutant brc1-2 partially rescued the mutant phenotype and decreased the number of branches, showing that these TCPs are functionally redundant in controlling branching. A transient dual-luciferase reporter assay indicated that GhTIE1 repressed the protein activity of GhBRC1 and GhTCP13, and thereby decreased the expression of their target gene GhHB21. Gene expression level analysis in GhTIE1-overexpressed and silenced plants also proved that GhTIE1 regulated shoot branching via repressing the activity of BRC1, HB21, HB40, and HB53. Our data reveals that shoot branching can be controlled via modulation of the activity of the TIE1 and TCP proteins and provides a theoretical basis for cultivating cotton varieties with ideal plant types.

Over-expression of an R2R3 MYB Gene, GhMYB73, increases tolerance to salt stress in transgenic Arabidopsis.

MYB family genes act as important regulators modulating the response to abiotic stress in plants. However, much less is known about MYB proteins in cotton. Here, we found that a cotton MYB gene, GhMYB73, was induced by NaCl and abscisic acid (ABA). Silencing GhMYB73 expression in cotton increased sensitivity to salt stress. The cotyledon greening rate of Arabidopsis thaliana over-expressing GhMYB73 under NaCl or mannitol treatment was significantly enhanced during the seedling germination stage. What's more, several osmotic stress-induced genes, such as AtNHX1, AtSOS3 and AtP5CS1, were more highly induced in the over-expression lines than in wild type under salt treatment, supporting the hypothesis that GhMYB73 contributes to salinity tolerance by improving osmotic stress resistance. Arabidopsis lines over-expressing GhMYB73 had superior germination and cotyledon greening under ABA treatment, and some abiotic stress-induced genes involved in ABA pathways (AtPYL8, AtABF3, AtRD29B and AtABI5), had increased transcription levels under salt-stress conditions in these lines. Furthermore, we found that GhMYB73 physically interacts with GhPYL8 and AtPYL8, suggesting that GhMYB73 regulates ABA signaling during salinity stress response. Taken together, over-expression of GhMYB73 significantly increases tolerance to salt and ABA stress, indicating that it can potentially be used in transgenic technology approaches to improve cotton salt tolerance.