Comparative genomics reveals cotton-specific virulence factors in flexible genomic regions in Verticillium dahliae and evidence of horizontal gene transfer from Fusarium

Comparative Genomics and Transcriptomics Analyses Reveal a Unique Environmental Adaptability of Vibrio fujianensis

The genus Vibrio is ubiquitous in marine environments and uses numerous evolutionary characteristics and survival strategies in order to occupy its niche. Here, a newly identified species, Vibrio fujianensis, was deeply explored to reveal a unique environmental adaptability. V. fujianensis type strain FJ201301^T shared 817 core genes with the Vibrio species in the population genomic analysis, but possessed unique genes of its own. In addition, V. fujianensis FJ201301^T was predicated to carry 106 virulence-related factors, several of which were mostly found in other pathogenic Vibrio species. Moreover, a comparative transcriptome analysis between the low-salt (1% NaCl) and high-salt (8% NaCl) condition was conducted to identify the genes involved in salt tolerance. A total of 913 unigenes were found to be differentially expressed. In a high-salt condition, 577 genes were significantly upregulated, whereas 336 unigenes were significantly downregulated. Notably, differentially expressed genes have a significant association with ribosome structural component and ribosome metabolism, which may play a role in salt tolerance. Transcriptional changes in ribosome genes indicate that V. fujianensis may have gained a predominant advantage in order to adapt to the changing environment. In conclusion, to survive in adversity, V. fujianensis has enhanced its environmental adaptability and developed various strategies to fill its niche.

Genome Sequences of Verticillium dahliae Defoliating Strain XJ592 and Nondefoliating Strain XJ511

Verticillium dahliae is a widely distributed soilborne pathogen that causes vascular wilt in more than 200 plant species. Defoliating and nondefoliating symptoms caused by the disease that result in either the loss or retention of leaves in infected plants, respectively, in hosts such as cotton, olive, and okra, divide the causal agent into defoliating and nondefoliating pathotypes. Our goal in this current work was to generate genome resources for the defoliating strain XJ592 and the nondefoliating strain XJ511 of V. dahliae isolated from cotton in China.

GhWRKY70D13 Regulates Resistance to Verticillium dahliae in Cotton Through the Ethylene and Jasmonic Acid Signaling Pathways

Verticillium wilt caused by Verticillium dahliae is a destructive cotton disease causing severe yield and quality losses worldwide. WRKY transcription factors play important roles in plant defense against pathogen infection. However, little has been reported on the functions of WRKYs in cotton's resistance to V. dahliae. Here, we identified 5, 5, and 10 WRKY70 genes in Gossypium arboreum, Gossypium raimondii, and Gossypium hirsutum, respectively, and investigated the expression profiles of all GhWRKY70 genes in various cotton tissues and in response to hormone treatment or V. dahliae infection. Reverse transcription-quantitative PCR analysis showed that GhWRKY70D13 was expressed higher in roots and stems than in other tissues, and up-regulated after V. dahliae inoculation. Knock-down of GhWRKY70D13 improved resistance to V. dahliae in both resistant and susceptible cotton cultivars. Comparative analysis of transcriptomes generated from wild-type and stable RNAi (RNA interference) plant with down-regulated GhWRKY70D13 showed that genes involved in ethylene (ET) and jasmonic acid (JA) biosynthesis and signaling were significantly upregulated in the GhWRKY70D13 RNAi plants. Consistently, the contents of 1-aminocyclopropane-1-carboxylic (ACC), JA, and JA-isoleucine levels were significantly higher in the GhWRKY70D13 RNAi plants than in wild-type. Following V. dahliae infection, the levels of ACC and JA decreased in the GhWRKY70D13 RNAi plants but still significantly higher (for ACC) than that in wild-type or at the same level (for JA) as in non-infected wild-type plants. Collectively, our results suggested that GhWRKY70D13 negatively regulates cotton's resistance to V. dahliae mainly through its effect on ET and JA biosynthesis and signaling pathways.

The Verticillium dahliae Sho1-MAPK pathway regulates melanin biosynthesis and is required for cotton infection

Verticillium dahliae is a soil‐borne fungus that causes vascular wilt on numerous plants worldwide. The fungus survives in the soil for up to 14 years by producing melanized microsclerotia. The protective function of melanin in abiotic stresses is well documented. Here, we found that the V. dahliae tetraspan transmembrane protein VdSho1, a homolog of the Saccharomyces cerevisiae Sho1, acts as an osmosensor, and is required for plant penetration and melanin biosynthesis. The deletion mutant ΔSho1 was incubated on a cellophane membrane substrate that mimics the plant epidermis, revealing that the penetration of ΔSho1 strain was reduced compared to the wild‐type strain. Furthermore, VdSho1 regulates melanin biosynthesis by a signalling mechanism requiring a kinase‐kinase signalling module of Vst50‐Vst11‐Vst7. Strains, ΔVst50, ΔVst7 and ΔVst11 also displayed defective penetration and melanin production like the ΔSho1 strain. Defects in penetration and melanin production in ΔSho1 were restored by overexpression of Vst50, suggesting that Vst50 lies downstream of VdSho1 in the regulatory pathway governing penetration and melanin biosynthesis. Data analyses revealed that the transmembrane portion of VdSho1 was essential for both membrane penetration and melanin production. This study demonstrates that Vst50‐Vst11‐Vst7 module regulates VdSho1‐mediated plant penetration and melanin production in V. dahliae, contributing to virulence.

Transcriptome sequencing of Verticillium dahliae from a cotton farm reveals positive correlation between virulence and tolerance of sugar-induced hyperosmosis

Verticillium dahliae causes disease symptoms in its host plants; however, due to its rapid variability, V. dahliae is difficult to control. To analyze the reason for this pathogenic differentiation, 22 V. dahliae strains with different virulence were isolated from a cotton farm. The genetic diversity of cotton varieties make cotton cultivars have different Verticillium wilt resistance, so the Xinluzao 7 (susceptible to V. dahliae), Zhongmian 35 (tolerant), and Xinluzao 33 (resistant) were used to investigate the pathogenicity of the strains in a green house. Vegetative compatibility groups (VCGs) assays, Internal Transcribed Spacer (ITS) PCR, and pathogenicity analysis showed that SHZ-4, SHZ-5, and SHZ-9 had close kinship and significantly different pathogenicity. Transcriptome sequencing of the three strains identified 19 of 146 unigenes in SHZ-4_vs_ SHZ-5, SHZ-5_vs_ SHZ-9, and SHZ-4_vs_ SHZ-9. In these unigenes, three proteinase and four polysaccharide degrading hydrolases were found to be associated with the pathogenicity. However, due to a number of differentially expressed genes in the transport, these unigenes not only played a role in nutrition absorption but might also contribute to the resistance of sugar-induced hyperosmosis. Moreover, the tolerance ability was positively related to the pathogenicity of V. dahliae. This resistance to sugar-induced hyperosmosis might help V. dahliae to access the nutrition of the host. The pathogenicity of V. dahliae correlated with the resistance of sugar-induced-hyperosmosis, which provides clues for the cultivation of V. dahliae resistant varieties.

Trans-kingdom RNA interactions drive the evolutionary arms race between hosts and pathogens

Trans-kingdom RNA plays a key role in host-parasite interactions. Hosts export specific endogenous microRNAs (miRNAs) into pathogens to target pathogen virulence genes and inhibit their invasion. In addition, trans-kingdom sRNAs produced by parasites may function as RNA effectors to suppress host immunity. Here, we summarize recent, important findings regarding trans-kingdom RNA and focus on the roles of trans-kingdom RNA in driving an evolutionary arms race between host and pathogen. We suggest that trans-kingdom RNA is a new platform for such arms races. Furthermore, we conjecture that trans-kingdom RNA contributes to horizontal gene transfer (HGT) involved in host-pathogen interactions. In addition, we propose that trans-kingdom RNA exchange and RNA driven HGT can have a great impact on the evolutionary ecology of interacting species.

VdOGDH is involved in energy metabolism and required for virulence of Verticillium dahliae

Verticillium dahliae, a soil-borne fungus, can invade plant vascular tissue and cause Verticillium wilt. The enzyme α-oxoglutarate dehydrogenase (OGDH), catalyzing the oxidation of α-oxoglutarate in the tricarboxylic acid cycle (TCA), is vital for energy metabolism in the fungi. Here, we identified the OGDH gene in V. dahliae (VdOGDH, VDAG_10018) and investigated its function in virulence by generating gene deletion mutants (ΔVdOGDH) and complementary mutants (ΔVdOGDH-C). When the ΔVdOGDH mutants were supplemented with different carbon sources, vegetative growth on Czapek Dox medium was significantly impaired, suggesting that VdOGDH is crucial for vegetative growth and carbon utilization. Conidia of the ΔVdOGDH mutants were atypically rounded or spherical, and hyphae were irregularly branched and lacked typical whorled branches. Mutants ΔVdOGDH-1 and ΔVdOGDH-2 were highly sensitive to H₂O₂ in the medium plates and had higher intracellular ROS levels. ΔVdOGDH mutants also had elevated expression of oxidative response-related genes, indicating that VdOGDH is involved in response to oxidative stress. In addition, the disruption of VdOGDH caused a significant increase in the expression of energy metabolism-related genes VdICL, VdICDH, VdMDH, and VdPDH and melanin-related genes Vayg1, VdSCD, VdLAC, VT4HR, and VaflM in the ΔVdOGDH mutants; thus, VdOGDH is also important for energy metabolism and melanin accumulation. Cotton plants inoculated with ΔVdOGDH mutants exhibited mild leaf chlorosis and the disease index was lower compared with wild type and ΔVdOGDH-C strains. These results together show that VdOGDH involved in energy metabolism of V. dahliae, is also essential for full virulence by regulating multiple fungal developmental factors.

The Gossypium hirsutum TIR-NBS-LRR gene GhDSC1 mediates resistance against Verticillium wilt

Improving genetic resistance is a preferred method to manage Verticillium wilt of cotton and other hosts. Identifying host resistance is difficult because of the dearth of resistance genes against this pathogen. Previously, a novel candidate gene involved in Verticillium wilt resistance was identified by a genome-wide association study using a panel of Gossypium hirsutum accessions. In this study, we cloned the candidate resistance gene from cotton that encodes a protein sharing homology with the TIR-NBS-LRR receptor-like defence protein DSC1 in Arabidopsis thaliana (hereafter named GhDSC1). GhDSC1 expressed at higher levels in response to Verticillium wilt and jasmonic acid (JA) treatment in resistant cotton cultivars as compared to susceptible cultivars and its product was localized to nucleus. The transfer of GhDSC1 to Arabidopsis conferred Verticillium resistance in an A. thaliana dsc1 mutant. This resistance response was associated with reactive oxygen species (ROS) accumulation and increased expression of JA-signalling-related genes. Furthermore, the expression of GhDSC1 in response to Verticillium wilt and JA signalling in A. thaliana displayed expression patterns similar to GhCAMTA3 in cotton under identical conditions, suggesting a coordinated DSC1 and CAMTA3 response in A. thaliana to Verticillium wilt. Analyses of GhDSC1 sequence polymorphism revealed a single nucleotide polymorphism (SNP) difference between resistant and susceptible cotton accessions, within the P-loop motif encoded by GhDSC1. This SNP difference causes ineffective activation of defence response in susceptible cultivars. These results demonstrated that GhDSC1 confers Verticillium resistance in the model plant system of A. thaliana, and therefore represents a suitable candidate for the genetic engineering of Verticillium wilt resistance in cotton.

The Gossypium hirsutum TIR-NBS-LRR gene GhDSC1 mediates resistance against Verticillium wilt

Improving genetic resistance is a preferred method to manage Verticillium wilt of cotton and other hosts. Identifying host resistance is difficult because of the dearth of resistance genes against this pathogen. Previously, a novel candidate gene involved in Verticillium wilt resistance was identified by a genome-wide association study using a panel of Gossypium hirsutum accessions. In this study, we cloned the candidate resistance gene from cotton that encodes a protein sharing homology with the TIR-NBS-LRR receptor-like defence protein DSC1 in Arabidopsis thaliana (hereafter named GhDSC1). GhDSC1 expressed at higher levels in response to Verticillium wilt and jasmonic acid (JA) treatment in resistant cotton cultivars as compared to susceptible cultivars and its product was localized to nucleus. The transfer of GhDSC1 to Arabidopsis conferred Verticillium resistance in an A. thaliana dsc1 mutant. This resistance response was associated with reactive oxygen species (ROS) accumulation and increased expression of JA-signalling-related genes. Furthermore, the expression of GhDSC1 in response to Verticillium wilt and JA signalling in A. thaliana displayed expression patterns similar to GhCAMTA3 in cotton under identical conditions, suggesting a coordinated DSC1 and CAMTA3 response in A. thaliana to Verticillium wilt. Analyses of GhDSC1 sequence polymorphism revealed a single nucleotide polymorphism (SNP) difference between resistant and susceptible cotton accessions, within the P-loop motif encoded by GhDSC1. This SNP difference causes ineffective activation of defence response in susceptible cultivars. These results demonstrated that GhDSC1 confers Verticillium resistance in the model plant system of A. thaliana, and therefore represents a suitable candidate for the genetic engineering of Verticillium wilt resistance in cotton.

Population genomics demystifies the defoliation phenotype in the plant pathogen Verticillium dahliae

Verticillium dahliae is a broad host-range pathogen that causes vascular wilts in plants. Interactions between three hosts and specific V. dahliae genotypes result in severe defoliation. The underlying mechanisms of defoliation are unresolved. Genome resequencing, gene deletion and complementation, gene expression analysis, sequence divergence, defoliating phenotype identification, virulence analysis, and quantification of V. dahliae secondary metabolites were performed. Population genomics previously revealed that G-LSR2 was horizontally transferred from the fungus Fusarium oxysporum f. sp. vasinfectum to V. dahliae and is exclusively found in the genomes of defoliating (D) strains. Deletion of seven genes within G-LSR2, designated as VdDf genes, produced the nondefoliation phenotype on cotton, olive, and okra but complementation of two genes restored the defoliation phenotype. Genes VdDf5 and VdDf6 associated with defoliation shared homology with polyketide synthases involved in secondary metabolism, whereas VdDf7 shared homology with proteins involved in the biosynthesis of N-lauroylethanolamine (N-acylethanolamine (NAE) 12:0), a compound that induces defoliation. NAE overbiosynthesis by D strains also appears to disrupt NAE metabolism in cotton by inducing overexpression of fatty acid amide hydrolase. The VdDfs modulate the synthesis and overproduction of secondary metabolites, such as NAE 12:0, that cause defoliation either by altering abscisic acid sensitivity, hormone disruption, or sensitivity to the pathogen.

Functional Characterization of Target of Rapamycin Signaling in Verticillium dahliae

More than 200 plants have been suffering from Verticillium wilt caused by Verticillium dahliae (V. dahliae) across the world. The target of rapamycin (TOR) is a lethal gene and controls cell growth and development in various eukaryotes, but little is known about TOR signaling in V. dahliae. Here, we found that V. dahliae strain is hypersensitive to rapamycin in the presence of rapamycin binding protein VdFKBP12 while the deletion mutant aaavdfkbp12 is insensitive to rapamycin. Heterologous expressing VdFKBP12 in Arabidopsis conferred rapamycin sensitivity, indicating that VdFKBP12 can bridge the interaction between rapamycin and TOR across species. The key across species of TOR complex 1 (TORC1) and TORC2 have been identified in V. dahliae, suggesting that TOR signaling pathway is evolutionarily conserved in eukaryotic species. Furthermore, the RNA-seq analysis showed that ribosomal biogenesis, RNA polymerase II transcription factors and many metabolic processes were significantly suppressed in rapamycin treated cells of V. dahliae. Importantly, transcript levels of genes associated with cell wall degrading enzymes (CWEDs) were dramatically down-regulated in TOR-inhibited cells. Further infection assay showed that the pathogenicity of V. dahliae and occurrence of Verticillium wilt can be blocked in the presence of rapamycin. These observations suggested that VdTOR is a key target of V. dahliae for controlling and preventing Verticillium wilt in plants.

The Genome Biology of Effector Gene Evolution in Filamentous Plant Pathogens

Filamentous pathogens, including fungi and oomycetes, pose major threats to global food security. Crop pathogens cause damage by secreting effectors that manipulate the host to the pathogen's advantage. Genes encoding such effectors are among the most rapidly evolving genes in pathogen genomes. Here, we review how the major characteristics of the emergence, function, and regulation of effector genes are tightly linked to the genomic compartments where these genes are located in pathogen genomes. The presence of repetitive elements in these compartments is associated with elevated rates of point mutations and sequence rearrangements with a major impact on effector diversification. The expression of many effectors converges on an epigenetic control mediated by the presence of repetitive elements. Population genomics analyses showed that rapidly evolving pathogens show high rates of turnover at effector loci and display a mosaic in effector presence-absence polymorphism among strains. We conclude that effective pathogen containment strategies require a thorough understanding of the effector genome biology and the pathogen's potential for rapid adaptation.

SNARE-Encoding Genes VdSec22 and VdSso1 Mediate Protein Secretion Required for Full Virulence in Verticillium dahliae

Proteins that mediate cellular and subcellular membrane fusion are key factors in vesicular trafficking in all eukaryotic cells, including the secretion and transport of plant pathogen virulence factors. In this study, we identified vesicle-fusion components that included 22 soluble N-ethylmaleimide-sensitive factor attachment protein receptors (SNAREs), four Sec1/Munc18 (SM) family proteins, and 10 Rab GTPases encoded in the genome of the vascular wilt pathogen Verticillium dahliae Vd991. Targeted deletion of two SNARE-encoding genes in V. dahliae, VdSec22 and VdSso1, significantly reduced virulence of both mutants on cotton, relative to the wild-type Vd991 strain. Comparative analyses of the secreted protein content (exoproteome) revealed that many enzymes involved in carbohydrate hydrolysis were regulated by VdSec22 or VdSso1. Consistent with a role of these enzymes in plant cell-wall degradation, pectin, cellulose, and xylan utilization were reduced in the VdSec22 or VdSso1 mutant strains along with a loss of exoproteome cytotoxic activity on cotton leaves. Comparisons with a pathogenicity-related exoproteome revealed that several known virulence factors were not regulated by VdSec22 or VdSso1, but some of the proteins regulated by VdSec22 or VdSso1 displayed different characteristics, including the lack of a typical signal peptide, suggesting that V. dahliae employs more than one secretory route to transport proteins to extracellular sites during infection.

The island cotton NBS-LRR gene GbaNA1 confers resistance to the non-race 1 Verticillium dahliae isolate Vd991

Wilt caused by Verticillium dahliae significantly reduces cotton yields, as host resistance in commercially cultivated Gossypium species is lacking. Understanding the molecular basis of disease resistance in non-commercial Gossypium species could galvanize the development of Verticillium wilt resistance in cultivated species. Nucleotide-binding site leucine-rich repeat (NBS-LRR) proteins play a central role in plant defence against pathogens. In this study, we focused on the relationship between a locus enriched with eight NBS-LRR genes and Verticillium wilt resistance in G. barbadense. Independent virus-induced gene silencing of each of the eight NBS-LRR genes in G. barbadense cultivar Hai 7124 revealed that silencing of GbaNA1 alone compromised the resistance of G. barbadense to V. dahliae isolate Vd991. In cultivar Hai 7124, GbaNA1 could be induced by V. dahliae isolate Vd991 and by ethylene, jasmonic acid and salicylic acid. Nuclear protein localization of GbaNA1 was demonstrated by transient expression. Sequencing of the GbaNA1 orthologue in nine G. hirsutum accessions revealed that all carried a non-functional allele, caused by a premature peptide truncation. In addition, all 10 G. barbadense and nine G. hirsutum accessions tested carried a full-length (∼1140 amino acids) homologue of the V. dahliae race 1 resistance gene Gbve1, although some sequence polymorphisms were observed. Verticillium dahliae Vd991 is a non-race 1 isolate that lacks the Ave1 gene. Thus, the resistance imparted by GbaNA1 appears to be mediated by a mechanism distinct from recognition of the fungal effector Ave1.

A Verticillium dahliae Extracellular Cutinase Modulates Plant Immune Responses

Cutinases have been implicated as important enzymes during the process of fungal infection of aerial plant organs. The function of cutinases in the disease cycle of fungal pathogens that invade plants through the roots has been less studied. Here, functional analysis of 13 cutinase (carbohydrate esterase family 5 domain-containing) genes (VdCUTs) in the highly virulent vascular wilt pathogen Verticillium dahliae Vd991 was performed. Significant sequence divergence in cutinase family members was observed in the genome of V. dahliae Vd991. Functional analyses demonstrated that only VdCUT11, as purified protein, induced cell death and triggered defense responses in Nicotiana benthamiana, cotton, and tomato plants. Virus-induced gene silencing showed that VdCUT11 induces plant defense responses in Nicotiana benthamania in a BAK1 and SOBIR-dependent manner. Furthermore, coinfiltration assays revealed that the carbohydrate-binding module family 1 protein (VdCBM1) suppressed VdCUT11-induced cell death and other defense responses in N. benthamiana. Targeted deletion of VdCUT11 in V. dahliae significantly compromised virulence on cotton plants. The cutinase VdCUT11 is an important secreted enzyme and virulence factor that elicits plant defense responses in the absence of VdCBM1.