Comprehensive analysis of Verticillium nonalfalfae in silico secretome uncovers putative effector proteins expressed during hop invasion

Novel Identification of the Collection of Pathogenic Fungal Species Verticillium with the Development of Species-Specific SSR Markers

The genus Verticillium is a group of ascomycete fungi that includes several pathogenic plant species. In 2011, a new taxonomic classification, proposed by Inderbitzin and coworkers (2011), re-defined the genus as Verticillium sensu stricto. The objective of our study was the re-classification of the fungal species held in the culture collection in the Slovenian Institute of Hop Research and Brewing in accordance with the newly established taxonomy. With the PCR marker system proposed by Inderbitzin and coworkers in 2011, we re-classified 88 Verticillium isolates out of the 105 samples that are held in the institute's bank, which were obtained from different geographic locations in Europe, North America, and Japan, and from different host plants, including alfalfa, cotton, hop, olive, potato, and tomato. However, the PCR marker for the V. dahliae identification proved to be less specific, and it resulted in the positive amplification of Gibellulopsis nigrescens, V. isaacii, and V. longisporum. To enable the accurate distinction of the fungi, the SSR and LAMP markers were added to the analyses. The 12 newly identified SSR markers, which were used in simplex PCR reactions or in combination, enabled the accurate identification of all included Verticillium isolates and could potentially be used as biomarkers for rapid and easy species identification.

Discovery of microRNA-like Small RNAs in Pathogenic Plant Fungus Verticillium nonalfalfae Using High-Throughput Sequencing and qPCR and RLM-RACE Validation

Verticillium nonalfalfae (V. nonalfalfae) is one of the most problematic hop (Humulus lupulus L.) pathogens, as the highly virulent fungal pathotypes cause severe annual yield losses due to infections of entire hop fields. In recent years, the RNA interference (RNAi) mechanism has become one of the main areas of focus in plant-fungal pathogen interaction studies and has been implicated as one of the major contributors to fungal pathogenicity. MicroRNA-like RNAs (milRNAs) have been identified in several important plant pathogenic fungi; however, to date, no milRNA has been reported in the V. nonalfalfae species. In the present study, using a high-throughput sequencing approach and extensive bioinformatics analysis, a total of 156 milRNA precursors were identified in the annotated V. nonalfalfae genome, and 27 of these milRNA precursors were selected as true milRNA candidates, with appropriate microRNA hairpin secondary structures. The stem-loop RT-qPCR assay was used for milRNA validation; a total of nine V. nonalfalfae milRNAs were detected, and their expression was confirmed. The milRNA expression patterns, determined by the absolute quantification approach, imply that milRNAs play an important role in the pathogenicity of highly virulent V. nonalfalfae pathotypes. Computational analysis predicted milRNA targets in the V. nonalfalfae genome and in the host hop transcriptome, and the activity of milRNA-mediated RNAi target cleavage was subsequently confirmed for two selected endogenous fungal target gene models using the 5' RLM-RACE approach.

Infection Process and Genome Assembly Provide Insights into the Pathogenic Mechanism of Destructive Mycoparasite Calcarisporium cordycipiticola with Host Specificity

Calcarisporium cordycipiticola is the pathogen in the white mildew disease of Cordyceps militaris, one of the popular mushrooms. This disease frequently occurs and there is no effective method for disease prevention and control. In the present study, C. militaris is found to be the only host of C. cordycipiticola, indicating strict host specificity. The infection process was monitored by fluorescent labeling and scanning and transmission electron microscopes. C. cordycipiticola can invade into the gaps among hyphae of the fruiting bodies of the host and fill them gradually. It can degrade the hyphae of the host by both direct contact and noncontact. The parasitism is initially biotrophic, and then necrotrophic as mycoparasitic interaction progresses. The approximate chromosome-level genome assembly of C. cordycipiticola yielded an N50 length of 5.45 Mbp and a total size of 34.51 Mbp, encoding 10,443 proteins. Phylogenomic analysis revealed that C. cordycipiticola is phylogenetically close to its specific host, C. militaris. A comparative genomic analysis showed that the number of CAZymes of C. cordycipiticola was much less than in other mycoparasites, which might be attributed to its host specificity. Secondary metabolite cluster analysis disclosed the great biosynthetic capabilities and potential mycotoxin production capability. This study provides insights into the potential pathogenesis and interaction between mycoparasite and its host.

In silico approach to predict pancreatic β-cells classically secreted proteins

Pancreatic β-cells, residents of the islets of Langerhans, are the unique insulin-producers in the body. Their physiology is a topic of intensive studies aiming to understand the biology of insulin production and its role in diabetes pathology. However, investigations about these cells' subset of secreted proteins, the secretome, are surprisingly scarce and a list describing islet/β-cell secretome upon glucose-stimulation is not yet available. In silico predictions of secretomes are an interesting approach that can be employed to forecast proteins likely to be secreted. In this context, using the rationale behind classical secretion of proteins through the secretory pathway, a Python tool capable of predicting classically secreted proteins was developed. This tool was applied to different available proteomic data (human and rodent islets, isolated β-cells, β-cell secretory granules, and β-cells supernatant), filtering them in order to selectively list only classically secreted proteins. The method presented here can retrieve, organize, search and filter proteomic lists using UniProtKB as a central database. It provides analysis by overlaying different sets of information, filtering out potential contaminants and clustering the identified proteins into functional groups. A range of 70-92% of the original proteomes analyzed was reduced generating predicted secretomes. Islet and β-cell signal peptide-containing proteins, and endoplasmic reticulum-resident proteins were identified and quantified. From the predicted secretomes, exemplary conservational patterns were inferred, as well as the signaling pathways enriched within them. Such a technique proves to be an effective approach to reduce the horizon of plausible targets for drug development or biomarkers identification.

Unraveling the Roles of Vascular Proteins Using Proteomics

Vascular bundles play important roles in transporting nutrients, growth signals, amino acids, and proteins between aerial and underground tissues. In order to understand these sophisticated processes, a comprehensive analysis of the roles of the components located in the vascular tissues is required. A great deal of data has been obtained from proteomic analyses of vascular tissues in plants, which mainly aim to identify the proteins moving through the vascular tissues. Here, different aspects of the phloem and xylem proteins are reviewed, including their collection methods, and their main biological roles in growth, and biotic and abiotic stress responses. The study of vascular proteomics shows great potential to contribute to our understanding of the biological mechanisms related to development and defense in plants.

Hop Polyphenols in Relation to Verticillium Wilt Resistance and Their Antifungal Activity

(1) Background: Verticillium wilt (VW) of hop is a devastating disease caused by the soil-borne fungi Verticillium nonalfalfae and Verticillium dahliae. As suggested by quantitative trait locus (QTL) mapping and RNA-Seq analyses, the underlying molecular mechanisms of resistance in hop are complex, consisting of preformed and induced defense responses, including the synthesis of various phenolic compounds. (2) Methods: We determined the total polyphenolic content at two phenological stages in roots and stems of 14 hop varieties differing in VW resistance, examined the changes in the total polyphenols of VW resistant variety Wye Target (WT) and susceptible Celeia (CE) on infection with V. nonalfalfae, and assessed the antifungal activity of six commercial phenolic compounds and total polyphenolic extracts from roots and stems of VW resistant WT and susceptible CE on the growth of two different V. nonalfalfae hop pathotypes. (3) Results: Generally, total polyphenols were higher in roots than stems and increased with maturation of the hop. Before flowering, the majority of VW resistant varieties had a significantly higher content of total polyphenols in stems than susceptible varieties. At the symptomatic stage of VW disease, total polyphenols decreased in VW resistant WT and susceptible CE plants in both roots and stems. The antifungal activity of total polyphenolic extracts against V. nonalfalfae was higher in hop extracts from stems than those from roots. Among the tested phenolic compounds, only p-coumaric acid and tyrosol markedly restricted fungal growth. (4) Conclusions: Although the correlation between VW resistance and total polyphenols content is not straightforward, higher levels of total polyphenols in the stems of the majority of VW resistant hop varieties at early phenological stages probably contribute to fast and efficient activation of signaling pathways, leading to successful defense against V. nonalfalfae infection.

Verticillium longisporum Elicits Media-Dependent Secretome Responses With Capacity to Distinguish Between Plant-Related Environments

Verticillia cause a vascular wilt disease affecting a broad range of economically valuable crops. The fungus enters its host plants through the roots and colonizes the vascular system. It requires extracellular proteins for a successful plant colonization. The exoproteomes of the allodiploid Verticillium longisporum upon cultivation in different media or xylem sap extracted from its host plant Brassica napus were compared. Secreted fungal proteins were identified by label free liquid chromatography-tandem mass spectrometry screening. V. longisporum induced two main secretion patterns. One response pattern was elicited in various non-plant related environments. The second pattern includes the exoprotein responses to the plant-related media, pectin-rich simulated xylem medium and pure xylem sap, which exhibited similar but additional distinct features. These exoproteomes include a shared core set of 221 secreted and similarly enriched fungal proteins. The pectin-rich medium significantly induced the secretion of 143 proteins including a number of pectin degrading enzymes, whereas xylem sap triggered a smaller but unique fungal exoproteome pattern with 32 enriched proteins. The latter pattern included proteins with domains of known pathogenicity factors, metallopeptidases and carbohydrate-active enzymes. The most abundant proteins of these different groups are the necrosis and ethylene inducing-like proteins Nlp2 and Nlp3, the cerato-platanin proteins Cp1 and Cp2, the metallopeptidases Mep1 and Mep2 and the carbohydrate-active enzymes Gla1, Amy1 and Cbd1. Their pathogenicity contribution was analyzed in the haploid parental strain V. dahliae. Deletion of the majority of the corresponding genes caused no phenotypic changes during ex planta growth or invasion and colonization of tomato plants. However, we discovered that the MEP1, NLP2, and NLP3 deletion strains were compromised in plant infections. Overall, our exoproteome approach revealed that the fungus induces specific secretion responses in different environments. The fungus has a general response to non-plant related media whereas it is able to fine-tune its exoproteome in the presence of plant material. Importantly, the xylem sap-specific exoproteome pinpointed Nlp2 and Nlp3 as single effectors required for successful V. dahliae colonization.

EffHunter: A Tool for Prediction of Effector Protein Candidates in Fungal Proteomic Databases

Pathogens are able to deliver small-secreted, cysteine-rich proteins into plant cells to enable infection. The computational prediction of effector proteins remains one of the most challenging areas in the study of plant fungi interactions. At present, there are several bioinformatic programs that can help in the identification of these proteins; however, in most cases, these programs are managed independently. Here, we present EffHunter, an easy and fast bioinformatics tool for the identification of effectors. This predictor was used to identify putative effectors in 88 proteomes using characteristics such as size, cysteine residue content, secretion signal and transmembrane domains.

Time-Course RNAseq Reveals Exserohilum turcicum Effectors and Pathogenicity Determinants

Exserohilum turcicum (sexual stage Setosphaeria turcica) is the hemibiotrophic causal agent of northern leaf blight of maize and sorghum. This study aimed to identify the genes involved in host colonization during the biotrophic and necrotrophic phases of infection. It also aimed to identify race-specific differences in gene expression. RNAseq of maize seedlings inoculated with a race 13N or 23N E. turcicum isolate was conducted before inoculation and at 2, 5, 7, and 13 days post-inoculation (dpi). Biological replicates were pooled per time point for each race and sequenced. A bioinformatics pipeline was used to identify candidate effectors, and expression was validated for selected candidates. Fungal biomass was positively correlated with the percentages of E. turcicum reads mapped, which were low at early time points (2-7 dpi) with a significant increase at 13 dpi, indicating a lifestyle switch from biotrophy to necrotrophy between 7 and 13 dpi. AVRHt1 is the putative E. turcicum effector recognized by the maize resistance gene Ht1. Consistent with this, AVRHt1 was expressed in planta by race 23N, but transcripts were absent in race 13N. In addition, specific transposable elements were expressed in 23N only. Genes encoding the virulence-associated peptidases leupeptin-inhibiting protein 1 and fungalysin were expressed in planta. Transcriptional profiles of genes involved in secondary metabolite synthesis or cell wall degradation revealed the importance of these genes during late stages of infection (13 dpi). A total of 346 expressed candidate effectors were identified, including Ecp6 and proteins similar to the secreted in xylem (SIX) effectors common to formae speciales of Fusarium oxysporum, SIX13 and SIX5. Expression profiling of Ecp6 and SIX13-like indicated a peak in expression at 5 and 7 dpi compared to 2 and 13 dpi. Sequencing of SIX13-like from diverse isolates of E. turcicum revealed host-specific polymorphisms that were mostly non-synonymous, resulting in two groups of SIX13-like proteins that corresponded to the maize or sorghum origin of each isolate. This study suggests putative mechanisms whereby E. turcicum causes disease. Identification of the candidate effector SIX13-like is consistent with the infection mode of E. turcicum through the xylem of susceptible hosts.

Flax rust infection transcriptomics reveals a transcriptional profile that may be indicative for rust Avr genes

Secreted effectors of fungal pathogens are essential elements for disease development. However, lack of sequence conservation among identified effectors has long been a problem for predicting effector complements in fungi. Here we have explored the expression characteristics of avirulence (Avr) genes and candidate effectors of the flax rust fungus, Melampsora lini. We performed transcriptome sequencing and real-time quantitative PCR (qPCR) on RNA extracted from ungerminated spores, germinated spores, isolated haustoria and flax seedlings inoculated with M. lini isolate CH5 during plant infection. Genes encoding two categories of M. lini proteins, namely Avr proteins and plant cell wall degrading enzymes (CWDEs), were investigated in detail. Analysis of the expression profiles of 623 genes encoding predicted secreted proteins in the M. lini transcriptome shows that the six known Avr genes (i.e. AvrM (avrM), AvrM14, AvrL2, AvrL567, AvrP123 (AvrP) and AvrP4) fall within a group of 64 similarly expressed genes that are induced in planta and show a peak of expression early in infection with a subsequent decline towards sporulation. Other genes within this group include two paralogues of AvrL2, an AvrL567 virulence allele, and a number of genes encoding putative effector proteins. By contrast, M. lini genes encoding CWDEs fall into different expression clusters with their distribution often unrelated to their catalytic activity or substrate targets. These results suggest that synthesis of M. lini Avr proteins may be regulated in a coordinated fashion and that the expression profiling-based analysis has significant predictive power for the identification of candidate Avr genes.

Chitin-Binding Protein of Verticillium nonalfalfae Disguises Fungus from Plant Chitinases and Suppresses Chitin-Triggered Host Immunity

During fungal infections, plant cells secrete chitinases, which digest chitin in the fungal cell walls. The recognition of released chitin oligomers via lysin motif (LysM)-containing immune host receptors results in the activation of defense signaling pathways. We report here that Verticillium nonalfalfae, a hemibiotrophic xylem-invading fungus, prevents these digestion and recognition processes by secreting a carbohydrate-binding motif 18 (CBM18)-chitin-binding protein, VnaChtBP, which is transcriptionally activated specifically during the parasitic life stages. VnaChtBP is encoded by the Vna8.213 gene, which is highly conserved within the species, suggesting high evolutionary stability and importance for the fungal lifestyle. In a pathogenicity assay, however, Vna8.213 knockout mutants exhibited wilting symptoms similar to the wild-type fungus, suggesting that Vna8.213 activity is functionally redundant during fungal infection of hop. In a binding assay, recombinant VnaChtBP bound chitin and chitin oligomers in vitro with submicromolar affinity and protected fungal hyphae from degradation by plant chitinases. Moreover, the chitin-triggered production of reactive oxygen species from hop suspension cells was abolished in the presence of VnaChtBP, indicating that VnaChtBP also acts as a suppressor of chitin-triggered immunity. Using a yeast-two-hybrid assay, circular dichroism, homology modeling, and molecular docking, we demonstrated that VnaChtBP forms dimers in the absence of ligands and that this interaction is stabilized by the binding of chitin hexamers with a similar preference in the two binding sites. Our data suggest that, in addition to chitin-binding LysM (CBM50) and Avr4 (CBM14) fungal effectors, structurally unrelated CBM18 effectors have convergently evolved to prevent hydrolysis of the fungal cell wall against plant chitinases and to interfere with chitin-triggered host immunity.

RNA interference core components identified and characterised in Verticillium nonalfalfae, a vascular wilt pathogenic plant fungi of hops

The conserved RNA interference mechanism (RNAi) in the fungal kingdom has become a focus of intense scientific investigation. The three catalytic core components, Dicer-like (DCL), Argonaute (AGO), and RNA-dependent RNA polymerase (RdRP), and their associated small interfering RNA molecules (siRNAs) have been identified and characterised in several fungal species. Recent studies have proposed that RNAi is a major contributor to the virulence of fungal pathogens as a result of so-called trans-kingdom RNA silencing. In the present study, we report on the existence of three core RNAi proteins in the pathogenic plant fungus Verticillium nonalfalfae, which is a soilborne plant pathogen that causes severe wilting disease in hops (Humulus lupulus L.). Two DCL proteins, two AGO proteins, and two RdRP proteins were identified, and their conserved RNAi domains were characterised. Our phylogeny results confirm the existing taxonomic relationships in the Ascomycete fungal phylum and show that the fungi of the Hypocreomycetidae subclass of the Sordariomycetes class have high amino acid sequence similarity. The expression analysis revealed a potential role of RNAi in the pathogenicity of the fungi, since all the RNAi genes were highly upregulated in the highly virulent isolate T2 and were also differentially expressed in the V. nonalfalfae-susceptible Celeia and V. nonalfalfae-resistant Wye Target cultivars.