Integrative transcriptomic and gene co-expression network analysis of host responses upon Verticillium dahliae infection in Gossypium hirsutum

SR45a plays a key role in enhancing cotton resistance to Verticillium dahliae by alternative splicing of immunity genes

Alternative splicing (AS) of pre-mRNAs increases the diversity of transcriptome and proteome and plays fundamental roles in plant development and stress responses. However, the prevalent changes in AS events and the regulating mechanisms of plants in response to pathogens remain largely unknown. Here, we show that AS changes are an important mechanism conferring cotton immunity to Verticillium dahliae (Vd). GauSR45a, encoding a serine/arginine-rich RNA binding protein, was upregulated expression and underwent AS in response to Vd infection in Gossypium australe, a wild diploid cotton species highly resistant to Vd. Silencing GauSR45a substantially reduced the splicing ratio of Vd-induced immune-associated genes, including GauBAK1 (BRI1-associated kinase 1) and GauCERK1 (chitin elicitor receptor kinase 1). GauSR45a binds to the GAAGA motif that is commonly found in the pre-mRNA of genes essential for PTI, ETI, and defense. The binding between GauSR45a and the GAAGA motif in the pre-mRNA of BAK1 was enhanced by two splicing factors of GauU2AF35B and GauU1-70 K, thereby facilitating exon splicing; silencing either AtU2AF35B or AtU1-70 K decreased the resistance to Vd in transgenic GauSR45a Arabidopsis. Overexpressing the short splicing variant of BAK1GauBAK1.1 resulted in enhanced Verticillium wilt resistance rather than the long one GauBAK1.2. Vd-induced far more AS events were in G. barbadense (resistant tetraploid cotton) than those in G. hirsutum (susceptible tetraploid cotton) during Vd infection, indicating resistance divergence in immune responses at a genome-wide scale. We provided evidence showing a fundamental mechanism by which GauSR45a enhances cotton resistance to Vd through global regulation of AS of immunity genes.

Genome-Wide and Expression Pattern Analysis of the HIT4 Gene Family Uncovers the Involvement of GHHIT4_4 in Response to Verticillium Wilt in Gossypium hirsutum

Chromatin remodelers are essential for regulating plant growth, development, and responses to environmental stresses. HIT4 (HEAT-INTOLERANT 4) is a novel stress-induced chromatin remodeling factor that has been less studied in abiotic stress and stress resistance, particularly in cotton. In this study, we conducted a comprehensive analysis of the members of the HIT4 gene family in Gossypium hirsutum using bioinformatics methods, including phylogenetic relationships, gene organization, transcription profiles, phylogenetic connections, selection pressure, and stress response. A total of 18 HIT4 genes were identified in four cotton species, with six HIT4 gene members in upland cotton. Based on the evolutionary relationships shown in the phylogenetic tree, the 18 HIT4 protein sequences were classified into four distinct subgroups. Furthermore, we conducted chromosome mapping to determine the genomic locations of these genes and visually represented the structural characteristics of HIT4 in G. hirsutum. In addition, we predicted the regulatory elements in HIT4 in G. hirsutum and conducted an analysis of repetitive sequences and gene collinearity among HIT4 in four cotton species. Moreover, we calculated the Ka/Ks ratio for homologous genes to assess the selection pressure acting on HIT4. Using RNA-seq, we explored the expression patterns of HIT4 genes in G. hirsutum and Gossypium barbadense. Through weighted gene co-expression network analysis (WGCNA), we found that GHHIT4_4 belonged to the MEblue module, which was mainly enriched in pathways such as DNA replication, phagosome, pentose and glucuronate interconversions, steroid biosynthesis, and starch and sucrose metabolism. This module may regulate the mechanism of upland cotton resistance to Verticillium wilt through DNA replication, phagosome, and various metabolic pathways. In addition, we performed heterologous overexpression of GH_D11G0591 (GHHIT4_4) in tobacco, and the results showed a significant reduction in disease index compared to the wild type, with higher expression levels of disease resistance genes in the transgenic tobacco. After conducting a VIGS (virus-induced gene silencing) experiment in cotton, the results indicated that silencing GHHIT4_4 had a significant impact, the resistance to Verticillium wilt weakened, and the internode length of the plants significantly decreased by 30.7% while the number of true leaves increased by 41.5%. qRT-PCR analysis indicated that GHHIT4_4 mainly enhanced cotton resistance to Verticillium wilt by indirectly regulating the PAL, 4CL, and CHI genes. The subcellular localization results revealed that GHHIT4_4 was predominantly distributed in the mitochondria and nucleus. This study offers preliminary evidence for the involvement of the GHHIT4_4 in cotton resistance to Verticillium wilt and lays the foundation for further research on the disease resistance mechanism of this gene in cotton.

Genome-wide identification of GhRLCK-VII subfamily genes in Gossypium hirsutum and investigation of their functions in resistance to Verticillium wilt

BACKGROUND

The receptor-like cytoplasmic kinases subfamily VII (RLCK-VII) is critical in regulating plant growth, development, and pattern-triggered immunity. However, a comprehensive exploration of these genes in the allotetraploid Gossypium hirsutum is still lacking. This study aimed to identify RLCK-VII genes in G. hirsutum and investigate their evolutionary history, structural features, expression patterns, and role in plant defense.

RESULTS

Seventy-two RLCK-VII genes in the G. hirsutum genome were unveiled and classified into nine groups following their phylogenetic analysis with Arabidopsis thaliana. Group VII-1 was the largest, accounting for 28%, while Groups VII-2 and VII-3 had only one member each. The analysis using MCScanX revealed that these 72 genes formed 166 collinear gene pairs and were resided on 26 chromosomes of G. hirsutum, suggesting that they were derived from whole genome segmental duplication events. Their calculated Ka/Ks values were below one, implying the occurrence of purification selection during the evolution and inhibition of gene function differentiation/loss. All members of the RLCK-VII subfamily possessed two conserved domains, PKinase-Tyr and PKinase, and several conserved PBS1 kinase subdomains, individually included in one of the ten motifs identified using MEME. The RNA-Seq results showed that RLCK-VII genes exhibited different spatiotemporal expression, indicating their involvement in cotton growth, development, and defense responses to Verticillium dahliae. The transcription patterns of RLCK-VII genes found by RNA-Seq were further validated using qRT-PCR assays after inoculating "20B12" (cotton cultivar) with "V991" (V. dahliae). The virus-induced gene silencing (VIGS) assays uncovered that two RLCK-VII genes (Gohir.A13G227248 and Gohir.A10G219900) were essential to G. hirsutum resistance to Verticillium wilt.

CONCLUSIONS

These observations offer valuable insight into the attributes and roles of RLCK-VII genes in G. hirsutum, potentially enable the breeding of new cotton cultivars with enhanced resistance to Verticillium wilt.

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.

Genome-wide identification and analysis of a cotton secretome reveals its role in resistance against Verticillium dahliae

BACKGROUND

The extracellular space between the cell wall and plasma membrane is a battlefield in plant-pathogen interactions. Within this space, the pathogen employs its secretome to attack the host in a variety of ways, including immunity manipulation. However, the role of the plant secretome is rarely studied for its role in disease resistance.

RESULTS

Here, we examined the secretome of Verticillium wilt-resistant Gossypium hirsutum cultivar Zhongzhimian No.2 (ZZM2, encoding 95,327 predicted coding sequences) to determine its role in disease resistance against the wilt causal agent, Verticillium dahliae. Bioinformatics-driven analyses showed that the ZZM2 genome encodes 2085 secreted proteins and that these display disequilibrium in their distribution among the chromosomes. The cotton secretome displayed differences in the abundance of certain amino acid residues as compared to the remaining encoded proteins due to the localization of these putative proteins in the extracellular space. The secretome analysis revealed conservation for an allotetraploid genome, which nevertheless exhibited variation among orthologs and comparable unique genes between the two sub-genomes. Secretome annotation strongly suggested its involvement in extracellular stress responses (hydrolase activity, oxidoreductase activity, and extracellular region, etc.), thus contributing to resistance against the V. dahliae infection. Furthermore, the defense response genes (immunity marker NbHIN1, salicylic acid marker NbPR1, and jasmonic acid marker NbLOX4) were activated to varying degrees when Nicotina benthamiana leaves were agro-infiltrated with 28 randomly selected members, suggesting that the secretome plays an important role in the immunity response. Finally, gene silencing assays of 11 members from 13 selected candidates in ZZM2 displayed higher susceptibility to V. dahliae, suggesting that the secretome members confer the Verticillium wilt resistance in cotton.

CONCLUSIONS

Our data demonstrate that the cotton secretome plays an important role in Verticillium wilt resistance, facilitating the development of the resistance gene markers and increasing the understanding of the mechanisms regulating disease resistance.

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.

Genome Resource for the Verticillium Wilt Resistant Gossypium hirsutum Cultivar Zhongzhimian No. 2

Verticillium wilt, caused by the fungal pathogen Verticillium dahliae, is the major cause of disease-related yield losses in cotton (Gossypium hirsutum). Despite these losses, the major cultivars of G. hirsutum remain highly susceptible to Verticillium wilt. The lack of understanding on the genetic basis for Verticillium wilt resistance may further hinder progress in deploying elite cultivars with proven resistance, such as the wilt resistant G. hirsutum cultivar Zhongzhimian No. 2. To help remedy this knowledge gap, we sequenced the whole genome of Zhongzhimian No. 2 and assembled it from a combination of PacBio long reads, Illumina short reads, and high-throughput chromosome conformation capture technologies. The final assembly of the genome was 2.33 Gb, encoding 95,327 predicted coding sequences. The GC content was 34.39% with 99.2% of the bases anchored to 26 pseudo-chromosomes that ranged from 53.8 to 127.7 Mb. This resource will help gain a detailed understanding of the genomic features governing high yield and Verticillium wilt resistance in this cultivar. Comparative genomics will be particularly helpful, since there are several published genomes of other Gossypium species. [Formula: see text] Copyright © 2022 The Author(s). This is an open access article distributed under the CC BY 4.0 International license.

Antifungal effects of volatile organic compounds produced by Trichoderma koningiopsis T2 against Verticillium dahliae

Volatile organic compounds (VOCs) produced by microorganisms are considered promising environmental-safety fumigants for controlling soil-borne diseases. Verticillium dahliae, a notorious fungal pathogen, causes economically important wilt diseases in agriculture and forestry industries. Here, we determined the antifungal activity of VOCs produced by Trichoderma koningiopsis T2. The VOCs from T. koningiopsis T2 were trapped by solid-phase microextraction (SPME) and tentatively identified through gas chromatography–mass spectrometry (GC/MS). The microsclerotia formation, cell wall-degrading enzymes and melanin synthesis of V. dahliae exposed to the VOC mixtures and selected single standards were examined. The results showed that the VOCs produced by strain T2 significantly inhibited the growth of V. dahliae mycelium and reduced the severity of Verticillium wilt in tobacco and cotton. Six individual compounds were identified in the volatilome of T. koningiopsis T2, and the dominant compounds were 3-octanone, 3-methyl-1-butanol, butanoic acid ethyl ester and 2-hexyl-furan. The VOCs of strain T2 exert a significant inhibitory effect on microsclerotia formation and decreased the activities of pectin lyase and endo-β-1,4-glucanase in V. dahliae. VOCs also downregulated the VdT3HR, VdT4HR, and VdSCD genes related to melanin synthesis by 29. 41-, 10. 49-, and 3.11-fold, respectively. Therefore, T. koningiopsis T2 has potential as a promising biofumigant for the biocontrol of Verticillium wilt disease.

Revealing Genetic Differences in Fiber Elongation between the Offspring of Sea Island Cotton and Upland Cotton Backcross Populations Based on Transcriptome and Weighted Gene Coexpression Networks

Fiber length is an important indicator of cotton fiber quality, and the time and rate of cotton fiber cell elongation are key factors in determining the fiber length of mature cotton. To gain insight into the differences in fiber elongation mechanisms in the offspring of backcross populations of Sea Island cotton Xinhai 16 and land cotton Line 9, we selected two groups with significant differences in fiber length (long-fiber group L and short-fiber group S) at different fiber development stages 0, 5, 10 and 15 days post-anthesis (DPA) for transcriptome comparison. A total of 171.74 Gb of clean data was obtained by RNA-seq, and eight genes were randomly selected for qPCR validation. Data analysis identified 6055 differentially expressed genes (DEGs) between two groups of fibers, L and S, in four developmental periods, and gene ontology (GO) term analysis revealed that these DEGs were associated mainly with microtubule driving, reactive oxygen species, plant cell wall biosynthesis, and glycosyl compound hydrolase activity. Kyoto encyclopedia of genes and genomes (KEGG) pathway analysis indicated that plant hormone signaling, mitogen-activated protein kinase (MAPK) signaling, and starch and sucrose metabolism pathways were associated with fiber elongation. Subsequently, a sustained upregulation expression pattern, profile 19, was identified and analyzed using short time-series expression miner (STEM). An analysis of the weighted gene coexpression network module uncovered 21 genes closely related to fiber development, mainly involved in functions such as cell wall relaxation, microtubule formation, and cytoskeletal structure of the cell wall. This study helps to enhance the understanding of the Sea Island–Upland backcross population and identifies key genes for cotton fiber development, and these findings will provide a basis for future research on the molecular mechanisms of fiber length formation in cotton populations.

Transcriptome Analysis of a Cotton Cultivar Provides Insights into the Differentially Expressed Genes Underlying Heightened Resistance to the Devastating Verticillium Wilt

Cotton is an important economic crop worldwide. Verticillium wilt (VW) caused by Verticillium dahliae (V. dahliae) is a serious disease in cotton, resulting in massive yield losses and decline of fiber quality. Breeding resistant cotton cultivars is an efficient but elaborate method to improve the resistance of cotton against V. dahliae infection. However, the functional mechanism of several excellent VW resistant cotton cultivars is poorly understood at present. In our current study, we carried out RNA-seq to discover the differentially expressed genes (DEGs) in the roots of susceptible cotton Gossypium hirsutum cultivar Junmian 1 (J1) and resistant cotton G.hirsutum cultivar Liaomian 38 (L38) upon Vd991 inoculation at two time points compared with the mock inoculated control plants. The potential function of DEGs uniquely expressed in J1 and L38 was also analyzed by GO enrichment and KEGG pathway associations. Most DEGs were assigned to resistance-related functions. In addition, resistance gene analogues (RGAs) were identified and analyzed for their role in the heightened resistance of the L38 cultivar against the devastating Vd991. In summary, we analyzed the regulatory network of genes in the resistant cotton cultivar L38 during V. dahliae infection, providing a novel and comprehensive insight into VW resistance in cotton.