Promotion of apoplastic oxidative burst by artificially selected GhCBSX3A enhances Verticillium dahliae resistance in upland cotton.

Verticillium wilt (VW) is a devasting disease affecting various plants, including upland cotton, a crucial fiber crop. Despite its impact, the genetic basis underlying cotton's susceptibility or defense against VW remains unclear. Here, we conducted a genome-wide association study on VW phenotyping in upland cotton and identified a locus on A13 that is significantly associated with VW resistance. We then identified a cystathionine ß-synthase domain gene at A13 locus, GhCBSX3A, which was induced by Verticillium dahliae. Functional analysis, including expression silencing in cotton and overexpression in Arabidopsis thaliana, confirmed that GhCBSX3A is a causal gene at the A13 locus, enhancing SAR-RBOHs-mediated apoplastic oxidative burst. We found allelic variation on the TATA-box of GhCBSX3A promoter attenuated its expression in upland cotton, thereby weakening VW resistance. Interestingly, we discovered that altered artificial selection of GhCBSX3A_R (an elite allele for VW) under different VW pressures during domestication and other improved processes allows specific human needs to be met. Our findings underscore the importance of GhCBSX3A in response to VW, and we propose a model for defense-associated genes being selected depending on the pathogen's pressure. The identified locus and gene serve as promising targets for VW resistance enhancement in cotton through genetic engineering.

Transcriptional and translational landscape fine-tune genome annotation and explores translation control in cotton.

INTRODUCTION: The unavailability of intergenic region annotation in whole genome sequencing and pan-genomics hinders efforts to enhance crop improvement. OBJECTIVES: Despite advances in research, the impact of post-transcriptional regulation on fiber development and translatome profiling at different stages of fiber growth in cotton (G. hirsutum) remains unexplored. METHODS: We utilized a combination of reference-guided de novo transcriptome assembly and ribosome profiling techniques to uncover the hidden mechanisms of translational control in eight distinct tissues of upland cotton. RESULTS: Our study identified P-site distribution at three-nucleotide periodicity and dominant ribosome footprint at 27 nucleotides. Specifically, we have detected 1,589 small open reading frames (sORFs), including 1,376 upstream ORFs (uORFs) and 213 downstream ORFs (dORFs), as well as 552 long non-coding RNAs (lncRNAs) with potential coding functions, which fine-tune the annotation of the cotton genome. Further, we have identified novel genes and lncRNAs with strong translation efficiency (TE), while sORFs were found to affect mRNA transcription levels during fiber elongation. The reliability of these findings was confirmed by the high consistency in correlation and synergetic fold change between RNA-sequencing (RNA-seq) and Ribosome-sequencing (Ribo-seq) analyses. Additionally, integrated omics analysis of the normal fiber ZM24 and short fiber pag1 cotton mutant revealed several differentially expressed genes (DEGs), and fiber-specific expressed (high/low) genes associated with sORFs (uORFs and dORFs). These findings were further supported by the overexpression and knockdown of GhKCS6, a gene associated with sORFs in cotton, and demonstrated the potential regulation of the mechanism governing fiber elongation on both the transcriptional and post-transcriptional levels. CONCLUSION: Reference-guided transcriptome assembly and the identification of novel transcripts fine-tune the annotation of the cotton genome and predicted the landscape of fiber development. Our approach provided a high-throughput method, based on multi-omics, for discovering unannotated ORFs, hidden translational control, and complex regulatory mechanisms in crop plants.

Oil candidate genes in seeds of cotton (Gossypium hirsutum L.) and functional validation of GhPXN1.

BACKGROUND: Cottonseed oil is a promising edible plant oil with abundant unsaturated fatty acids. However, few studies have been conducted to explore the characteristics of cottonseed oil. The molecular mechanism of cottonseed oil accumulation remains unclear. RESULTS: In the present study, we conducted comparative transcriptome and weighted gene co-expression network (WGCNA) analysis for two G. hirsutum materials with significant difference in cottonseed oil content. Results showed that, between the high oil genotype 6053 (H6053) and the low oil genotype 2052 (L2052), a total of 412, 507, 1,121, 1,953, and 2,019 differentially expressed genes (DEGs) were detected at 10, 15, 20, 25, and 30 DPA, respectively. Remarkably, a large number of the down-regulated DEGs were enriched in the phenylalanine metabolic processes. Investigation into the dynamic changes of expression profiling of genes associated with both phenylalanine metabolism and oil biosynthesis has shed light on a significant competitive relationship in substrate allocation during cottonseed development. Additionally, the WGCNA analysis of all DEGs identified eight distinct modules, one of which includes GhPXN1, a gene closely associated with oil accumulation. Through phylogenetic analysis, we hypothesized that GhPXN1 in G. hirsutum might have been introgressed from G. arboreum. Overexpression of the GhPXN1 gene in tobacco leaf suggested a significant reduction in oil content compared to the empty-vector transformants. Furthermore, ten other crucial oil candidate genes identified in this study were also validated using quantitative real-time PCR (qRT-PCR). CONCLUSIONS: Overall, this study enhances our comprehension of the molecular mechanisms underlying cottonseed oil accumulation.

Uncovering genomic and transcriptional variations facilitates utilization of wild resources in cotton disease resistance improvement.

BACKGROUND: Upland cotton wild/landraces represent a valuable resource for disease resistance alleles. Genetic differentiation between genotypes, as well as variation in Verticillium wilt (VW) resistance, has been poorly characterized for upland cotton accessions on the domestication spectrum (from wild/landraces to elite lines). RESULTS: To illustrate the effects of modern breeding on VW resistance in upland cotton, 37 wild/landraces were resequenced and phenotyped for VW resistance. Genomic patterns of differentiation were identified between wild/landraces and improved upland cotton, and a significant decline in VW resistance was observed in association with improvement. Four genotypes representing different degrees of improvement were used in a full-length transcriptome analysis to study the genetic basis of VW resistance. ROS signaling was highly conserved at the transcriptional level, likely providing the basis for VW resistance in upland cotton. ASN biosynthesis and HSP90-mediated resistance moderated the response to VW in wild/landraces, and loss of induction activity of these genes resulted in VW susceptibility. The observed genomic differentiation contributed to the loss of induction of some important VW resistance genes such as HSP90.4 and PR16. CONCLUSIONS: Besides providing new insights into the evolution of upland cotton VW resistance, this study also identifies important resistance pathways and genes for both fundamental research and cotton breeding.

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.

Synergistic optimization of crops by combining early maturation with other agronomic traits.

Many newly created early maturing varieties exhibit poor stress resistance and low yield, whereas stress-resistant varieties are typically late maturing. For this reason, the polymerization of early maturity and other desired agronomic qualities requires overcoming the negative connection between early maturity, multi-resistance, and yield, which presents a formidable challenge in current breeding techniques. We review the most salient constraints of early maturity breeding in current crop planting practices and the molecular mechanisms of different maturation timeframes in diverse crops from their origin center to production areas. We explore current breeding tactics and the future direction of crop breeding and the issues that must be resolved to accomplish the polymerization of desirable traits in light of the current obstacles and limitations.

Genome-wide association analysis reveals a novel pathway mediated by a dual-TIR domain protein for pathogen resistance in cotton.

BACKGROUND: Verticillium wilt is one of the most devasting diseases for many plants, leading to global economic loss. Cotton is known to be vulnerable to its fungal pathogen, Verticillium dahliae, yet the related genetic mechanism remains unknown. RESULTS: By genome-wide association studies of 419 accessions of the upland cotton, Gossypium hirsutum, we identify ten loci that are associated with resistance against Verticillium wilt. Among these loci, SHZDI1/SHZDP2/AYDP1 from chromosome A10 is located on a fragment introgressed from Gossypium arboreum. We characterize a large cluster of Toll/interleukin 1 (TIR) nucleotide-binding leucine-rich repeat receptors in this fragment. We then identify a dual-TIR domain gene from this cluster, GhRVD1, which triggers an effector-independent cell death and is induced by Verticillium dahliae. We confirm that GhRVD1 is one of the causal gene for SHZDI1. Allelic variation in the TIR domain attenuates GhRVD1-mediated resistance against Verticillium dahliae. Homodimerization between TIR1-TIR2 mediates rapid immune response, while disruption of its αD- and αE-helices interface eliminates the autoactivity and self-association of TIR1-TIR2. We further demonstrate that GhTIRP1 inhibits the autoactivity and self-association of TIR1-TIR2 by competing for binding to them, thereby preventing the resistance to Verticillium dahliae. CONCLUSIONS: We propose the first working model for TIRP1 involved self-association and autoactivity of dual-TIR domain proteins that confer compromised pathogen resistance of dual-TIR domain proteins in plants. The findings reveal a novel mechanism on Verticillium dahliae resistance and provide genetic basis for breeding in future.

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.

Efficient genotype-independent cotton genetic transformation and genome editing.

Cotton (Gossypium spp.) is one of the most important fiber crops worldwide. In the last two decades, transgenesis and genome editing have played important roles in cotton improvement. However, genotype dependence is one of the key bottlenecks in generating transgenic and gene-edited cotton plants through either particle bombardment or Agrobacterium-mediated transformation. Here, we developed a shoot apical meristem (SAM) cell-mediated transformation system (SAMT) that allowed the transformation of recalcitrant cotton genotypes including widely grown upland cotton (Gossypium hirsutum), Sea island cotton (Gossypium barbadense), and Asiatic cotton (Gossypium arboreum). Through SAMT, we successfully introduced two foreign genes, GFP and RUBY, into SAM cells of some recalcitrant cotton genotypes. Within 2-3 months, transgenic adventitious shoots generated from the axillary meristem zone could be recovered and grown into whole cotton plants. The GFP fluorescent signal and betalain accumulation could be observed in various tissues in GFP- and RUBY-positive plants, as well as in their progenies, indicating that the transgenes were stably integrated into the genome and transmitted to the next generation. Furthermore, using SAMT, we successfully generated edited cotton plants with inheritable targeted mutagenesis in the GhPGF and GhRCD1 genes through CRISPR/Cas9-mediated genome editing. In summary, the established SAMT transformation system here in this study bypasses the embryogenesis process during tissue culture in a conventional transformation procedure and significantly accelerates the generation of transgenic and gene-edited plants for genetic improvement of recalcitrant cotton varieties.

CottonMD: a multi-omics database for cotton biological study.

Cotton is an important economic crop, and many loci for important traits have been identified, but it remains challenging and time-consuming to identify candidate or causal genes/variants and clarify their roles in phenotype formation and regulation. Here, we first collected and integrated the multi-omics datasets including 25 genomes, transcriptomes in 76 tissue samples, epigenome data of five species and metabolome data of 768 metabolites from four tissues, and genetic variation, trait and transcriptome datasets from 4180 cotton accessions. Then, a cotton multi-omics database (CottonMD, http://yanglab.hzau.edu.cn/CottonMD/) was constructed. In CottonMD, multiple statistical methods were applied to identify the associations between variations and phenotypes, and many easy-to-use analysis tools were provided to help researchers quickly acquire the related omics information and perform multi-omics data analysis. Two case studies demonstrated the power of CottonMD for identifying and analyzing the candidate genes, as well as the great potential of integrating multi-omics data for cotton genetic breeding and functional genomics research.

Functional divergence of GLP genes between G. barbadense and G. hirsutum in response to Verticillium dahliae infection.

Germin-like proteins (GLPs) play important roles in plant disease resistance but are rarely reported in cotton. We compared the expression of GLPs in Verticillium dahliae inoculate G. hirsutum (susceptible) and G. barbadense (resistant) and enriched 11 differentially expressed GLPs. 2741 GLP proteins identified from 53 species determined that GLP probably originated from algae and could be classified into 7 clades according to phylogenetic analysis, among which Clade I is likely the most ancient. Cotton GLP (two allopolyploids and two diploids) genes within a shared clade were highly conserved. Intriguingly, clade VII genes were mainly located in gene clusters that derived from the expansion of LTR transposons. Clade VII members expressed mainly in root which is the first battle against Verticillium dahlia and could be induced more intensely in G. barbadense than G. hirsutum. The GLP genes are resistant to Verticillium dahliae, which can be further investigated against Verticillium wilt.

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.

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.

Regulatory Network of Cotton Genes in Response to Salt, Drought and Wilt Diseases (Verticillium and Fusarium): Progress and Perspective.

In environmental conditions, crop plants are extremely affected by multiple abiotic stresses including salinity, drought, heat, and cold, as well as several biotic stresses such as pests and pathogens. However, salinity, drought, and wilt diseases (e.g., Fusarium and Verticillium) are considered the most destructive environmental stresses to cotton plants. These cause severe growth interruption and yield loss of cotton. Since cotton crops are central contributors to total worldwide fiber production, and also important for oilseed crops, it is essential to improve stress tolerant cultivars to secure future sustainable crop production under adverse environments. Plants have evolved complex mechanisms to respond and acclimate to adverse stress conditions at both physiological and molecular levels. Recent progresses in molecular genetics have delivered new insights into the regulatory network system of plant genes, which generally includes defense of cell membranes and proteins, signaling cascades and transcriptional control, and ion uptake and transport and their relevant biochemical pathways and signal factors. In this review, we mainly summarize recent progress concerning several resistance-related genes of cotton plants in response to abiotic (salt and drought) and biotic (Fusarium and Verticillium wilt) stresses and classify them according to their molecular functions to better understand the genetic network. Moreover, this review proposes that studies of stress related genes will advance the security of cotton yield and production under a changing climate and that these genes should be incorporated in the development of cotton tolerant to salt, drought, and fungal wilt diseases (Verticillium and Fusarium).

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.

BIN2 negatively regulates plant defence against Verticillium dahliae in Arabidopsis and cotton.

Verticillium wilt is caused by the soil-borne vascular pathogen Verticillium dahliae, and affects a wide range of economically important crops, including upland cotton (Gossypium hirsutum). Previous studies showed that expression levels of BIN2 were significantly down-regulated during infestation with V. dahliae. However, the underlying molecular mechanism of BIN2 in plant regulation against V. dahliae remains enigmatic. Here, we characterized a protein kinase GhBIN2 from Gossypium hirsutum, and identified GhBIN2 as a negative regulator of resistance to V. dahliae. The Verticillium wilt resistance of Arabidopsis and cotton were significantly enhanced when BIN2 was knocked down. Constitutive expression of BIN2 attenuated plant resistance to V. dahliae. We found that BIN2 regulated plant endogenous JA content and influenced the expression of JA-responsive marker genes. Further analysis revealed that BIN2 interacted with and phosphorylated JAZ family proteins, key repressors of the JA signalling pathway in both Arabidopsis and cotton. Spectrometric analysis and site-directed mutagenesis showed that BIN2 phosphorylated AtJAZ1 at T196, resulting in the degradation of JAZ proteins. Collectively, these results show that BIN2 interacts with JAZ proteins and plays a negative role in plant resistance to V. dahliae. Thus, BIN2 may be a potential target gene for genetic engineering against Verticillium wilt in crops.

Freeze substitution Hi-C, a convenient and cost-effective method for capturing the natural 3D chromatin conformation from frozen samples.

Chromatin interactions functionally affect genome architecture and gene regulation, but to date, only fresh samples must be used in High-through chromosome conformation capture (Hi-C) to keep natural chromatin conformation intact. This requirement has impeded the advancement of 3D genome research by limiting sample collection and storage options for researchers and severely limiting the number of samples that can be processed in a short time. Here, we develop a freeze substitution Hi-C (FS-Hi-C) technique that overcomes the need for fresh samples. FS-Hi-C can be used with samples stored in liquid nitrogen (LN2): the water in a vitreous form in the sample cells is replaced with ethanol via automated freeze substitution. After confirming that the FS step preserves the natural chromosome conformation during sample thawing, we tested the performance of FS-Hi-C with Drosophila melanogaster and Gossypium hirsutum. Beyond allowing the use of frozen samples and confirming that FS-Hi-C delivers robust data for generating contact heat maps and delineating A/B compartments and topologically associating domains, we found that FS-Hi-C outperforms the in situ Hi-C in terms of library quality, reproducibility, and valid interactions. Thus, FS-Hi-C will probably extend the application of 3D genome structure analysis to the vast number of experimental contexts in biological and medical research for which Hi-C methods have been unfeasible to date.

Genetic and evolution analysis of extrafloral nectary in cotton.

Extrafloral nectaries are a defence trait that plays important roles in plant-animal interactions. Gossypium species are characterized by cellular grooves in leaf midribs that secret large amounts of nectar. Here, with a panel of 215 G. arboreum accessions, we compared extrafloral nectaries to nectariless accessions to identify a region of Chr12 that showed strong differentiation and overlapped with signals from GWAS of nectaries. Fine mapping of an F2 population identified GaNEC1, encoding a PB1 domain-containing protein, as a positive regulator of nectary formation. An InDel, encoding a five amino acid deletion, together with a nonsynonymous substitution, was predicted to cause 3D structural changes in GaNEC1 protein that could confer the nectariless phenotype. mRNA-Seq analysis showed that JA-related genes are up-regulated and cell wall-related genes are down-regulated in the nectary. Silencing of GaNEC1 led to a smaller size of foliar nectary phenotype. Metabolomics analysis identified more than 400 metabolites in nectar, including expected saccharides and amino acids. The identification of GaNEC1 helps establish the network regulating nectary formation and nectar secretion, and has implications for understanding the production of secondary metabolites in nectar. Our results will deepen our understanding of plant-mutualism co-evolution and interactions, and will enable utilization of a plant defence trait in cotton breeding efforts.

Gossypium Genomics: Trends, Scope, and Utilization for Cotton Improvement.

Cotton (Gossypium spp.) is the most important natural fiber crop worldwide. The diversity of Gossypium species also provides an ideal model for investigating evolution and domestication of polyploids. However, the huge and complex cotton genome hinders genomic research. Technical advances in high-throughput sequencing and bioinformatics analysis have now largely overcome these obstacles, bringing about a new era of cotton genomics. Here, we review recent progress in Gossypium genomics based on whole genome sequencing, resequencing, and comparative genomics, which have provided insights about the genomic basis of fiber biogenesis and the landscape of cotton functional genomics. We address current challenges and present multidisciplinary genomics-enabled breeding strategies covering the breadth of high fiber yield, quality, and environmental resilience for future cotton breeding programs.