A comprehensive analysis of genes encoding small secreted proteins identifies candidate effectors in Melampsora larici-populina (poplar leaf rust)

Transposon-Associated Small RNAs Involved in Plant Defense in Poplar

Utilizing high-throughput Illumina sequencing, we examined how small RNA (sRNA) profiles vary in Chinese white poplar (Populus tomentosa) across two pivotal infection stages by the rust fungus Melampsora larici-populina: the biotrophic growth phase (T02; 48 h post infection) and the urediniospore development and dispersal phase (T03; 168 h), both essential for plant colonization and prolonged biotrophic engagement. Far exceeding random expectations, siRNA clusters predominantly arose from transposon regions, with pseudogenes also contributing significantly, and infection-stage-specific variations were notably tied to these transposon-derived siRNAs. As the infection advanced, clusters of 24 nt siRNAs in transposon and intergenic regions exhibited pronounced abundance shifts. An analysis of targets indicated that Populus sRNAs potentially regulate 95% of Melampsora larici-populina genes, with pathogen effector genes showing heightened targeting by sRNAs during the biotrophic and urediniospore phases compared to controls, pointing to selective sRNA-target interactions. In contrast to conserved miRNAs across plant species, Populus-specific miRNAs displayed a markedly greater tendency to target NB-LRR genes. These observations collectively highlight the innovative roles of sRNAs in plant defense, their evolutionary roots, and their dynamic interplay with pathogen coevolution.

Host-Specific and Homologous Pairs of Melampsora larici-populina Effectors Unveil Novel Nicotiana benthamiana Stromule Induction Factors

The poplar rust fungus Melampsora larici-populina is part of one of the most devastating group of fungi (Pucciniales) and causes important economic losses to the poplar industry. Because M. larici-populina is a heteroecious obligate biotroph, its spread depends on its ability to carry out its reproductive cycle through larch and then poplar parasitism. Genomic approaches have identified more than 1,000 candidate secreted effector proteins (CSEPs) from the predicted secretome of M. larici-populina that are potentially implicated in the infection process. In this study, we selected CSEP pairs (and one triplet) among CSEP gene families that share high sequence homology but display specific gene expression profiles among the two distinct hosts. We determined their subcellular localization by confocal microscopy through expression in the heterologous plant system Nicotiana benthamiana. Five out of nine showed partial or complete chloroplastic localization. We also screened for potential protein interactors from larch and poplar by yeast two-hybrid assays. One pair of CSEPs and the triplet shared common interactors, whereas the members of the two other pairs did not have common targets from either host. Finally, stromule induction quantification revealed that two pairs and the triplet of CSEPs induced stromules when transiently expressed in N. benthamiana. The use of N. benthamiana eds1 and nrg1 knockout lines showed that CSEPs can induce stromules through an eds1-independent mechanism. However, CSEP homologs shared the same impact on stromule induction and contributed to discovering a new stromule induction cascade that can be partially and/or fully independent of eds1. [Formula: see text] Copyright © 2024 The Author(s). This is an open access article distributed under the CC BY-NC-ND 4.0 International license.

Suppression of Chitin-Triggered Immunity by Plant Fungal Pathogens: A Case Study of the Cucurbit Powdery Mildew Fungus Podosphaera xanthii

Fungal pathogens are significant plant-destroying microorganisms that present an increasing threat to the world's crop production. Chitin is a crucial component of fungal cell walls and a conserved MAMP (microbe-associated molecular pattern) that can be recognized by specific plant receptors, activating chitin-triggered immunity. The molecular mechanisms underlying the perception of chitin by specific receptors are well known in plants such as rice and Arabidopsis thaliana and are believed to function similarly in many other plants. To become a plant pathogen, fungi have to suppress the activation of chitin-triggered immunity. Therefore, fungal pathogens have evolved various strategies, such as prevention of chitin digestion or interference with plant chitin receptors or chitin signaling, which involve the secretion of fungal proteins in most cases. Since chitin immunity is a very effective defensive response, these fungal mechanisms are believed to work in close coordination. In this review, we first provide an overview of the current understanding of chitin-triggered immune signaling and the fungal proteins developed for its suppression. Second, as an example, we discuss the mechanisms operating in fungal biotrophs such as powdery mildew fungi, particularly in the model species Podosphaera xanthii, the main causal agent of powdery mildew in cucurbits. The key role of fungal effector proteins involved in the modification, degradation, or sequestration of immunogenic chitin oligomers is discussed in the context of fungal pathogenesis and the promotion of powdery mildew disease. Finally, the use of this fundamental knowledge for the development of intervention strategies against powdery mildew fungi is also discussed.

The secreted FolAsp aspartic protease facilitates the virulence of Fusarium oxysporum f. sp. lycopersici

Pathogens utilize secretory effectors to manipulate plant defense. Fusarium oxysporum f. sp. lycopersici (Fol) is the causal agent of Fusarium wilt disease in tomatoes. We previously identified 32 secreted effector candidates by LC-MS analysis. In this study, we functionally identified one of the secreted proteins, FolAsp, which belongs to the aspartic proteases (Asp) family. The FolAsp was upregulated with host root specifically induction. Its N-terminal 1-19 amino acids performed the secretion activity in the yeast system, which supported its secretion in Fol. Phenotypically, the growth and conidia production of the FolAsp deletion mutants were not changed; however, the mutants displayed significantly reduced virulence to the host tomato. Further study revealed the FolAsp was localized at the apoplast and inhibited INF1-induced cell death in planta. Meanwhile, FolAsp could inhibit flg22-mediated ROS burst. Furthermore, FolAsp displayed protease activity on host protein, and overexpression of FolAsp in Fol enhanced pathogen virulence. These results considerably extend our understanding of pathogens utilizing secreted protease to inhibit plant defense and promote its virulence, which provides potential applications for tomato improvement against disease as the new drug target.

The secreted FoAPY1 peptidase promotes Fusarium oxysporum invasion

The secretion of peptidases from several pathogens has been reported, but the biological function of these proteins in plant-pathogen interactions is poorly understood. Fusarium oxysporum, a soil-borne plant pathogenic fungus that causes Fusarium wilt in its host, can secrete proteins into host plant cells during the infection process to interfere with the host plant defense response and promote disease occurrence. In this study, we identified a peptidase, FoAPY1, that could be secreted from F. oxysporum depending on the N-terminal signal peptide of the protein. FoAPY1 belongs to the peptidase M28 family and exerts peptidase activity in vitro. Furthermore, the FoAYP1 gene knockout strain (∆FoAYP1) presented reduced virulence to tomato plants, but its mycelial growth and conidiation were unchanged. Moreover, FoAYP1 overexpression tomato seedlings exhibited enhanced susceptibility to F. oxysporum and Botrytis cinerea strains. These data demonstrated that FoAYP1 contributes to the virulence of F. oxysporum may through peptidase activity against host plant proteins.

Suppression of Chitin-Triggered Immunity by a New Fungal Chitin-Binding Effector Resulting from Alternative Splicing of a Chitin Deacetylase Gene

Phytopathogenic fungi have evolved mechanisms to manipulate plant defences, such as chitin-triggered immunity, a plant defensive response based on the recognition of chitin oligomers by plant-specific receptors. To cope with chitin resistance, fungal pathogens have developed different strategies to prevent chitin recognition, such as binding, breaking, or modifying immunogenic oligomers. In powdery mildew fungi, the activity of chitin deacetylase (CDA) is crucial for this purpose, since silencing of the CDA gene leads to a rapid activation of chitin signalling and the subsequent suppression of fungal growth. In this work, we have identified an unusually short CDA transcript in Podosphaera xanthii, the cucurbit powdery mildew pathogen. This transcript, designated PxCDA3, appears to encode a truncated version of CDA resulting from an alternative splicing of the PxCDA gene, which lacked most of the chitin deacetylase activity domain but retained the carbohydrate-binding module. Experiments with the recombinant protein showed its ability to bind to chitin oligomers and prevent the activation of chitin signalling. Furthermore, the use of fluorescent fusion proteins allowed its localization in plant papillae at pathogen penetration sites. Our results suggest the occurrence of a new fungal chitin-binding effector, designated CHBE, involved in the manipulation of chitin-triggered immunity in powdery mildew fungi.

A decade after the first Pucciniales genomes: A bibliometric snapshot of (post) genomics studies in three model rust fungi

Pucciniales (rust fungi) are one of the largest fungal order of plant pathogens. They collectively infect key crops such as wheat and soybean, and threaten global food security. In the early 2010s, the genome sequences of three rust fungi were released: Melampsora larici-populina (the poplar leaf rust fungus), Puccinia graminis f. sp. tritici (the wheat stem rust fungus), and Puccinia striiformis f. sp. triciti (the wheat stripe rust or wheat yellow rust fungus). The availability of those genomes has forwarded rust biology into the post-genomic era, sparking a series of genomics, transcriptomics, in silico, and functional studies. Here, we snapshot the last 10 years of post-genomics studies addressing M. larici-populina, P. graminis f. sp. tritici, and/or P. striiformis f. sp. tritici. This mini-review notably reveals the model species-centered structure of the research community, and highlights the drastic increase of the number of functional studies focused on effectors since 2014, which notably revealed chloroplasts as a central host compartment targeted by rust fungi. This mini-review also discusses genomics-facilitated studies in other rust species, and emerging post-genomic research trends related to fully-phased rust genomes.

Reference Genes Selection of Gymnosporangium yamadae during the Interaction with Apple Leaves

Apple rust disease caused by Gymnosporangium yamadae is the one of the major threats to the development of the apple industry in China, but the pathogenic molecular mechanism of the disease remains unclear. It is imperative to screen out appropriate reference genes during the interaction between G. yamadae and apple leaves to analyze the gene expression patterns during the pathogenesis of G. yamadae. ACT, EF1, EF2, GAPDH, 40S, 60S, α-TUB, β-TUB and UBCE3 were selected as candidate reference genes based on the transcriptomic dataset of G. yamadae. The expression levels were tested by real-time quantitative PCR during time-course infection of apple leaves and the expression stabilities were evaluated by △Ct method as well as by three software (NormFinder, geNorm and BestKeeper) and one web-based analysis software (RefFinder). The expression stability of the candidate reference genes was further validated by using the effector candidate gene Cluster-3395.48660 as the target gene in RT-qPCR. According to the results by △Ct and BestKeeper, 40S, EF2 and EF1 were the most stable reference genes, while EF1, EF2 and GAPDH were the most stable reference genes based on the NormFinder analysis result. The geNorm recommended the most stable genes EF1, EF2 and α-TUB as reference genes. Comprehensive analysis results of the RefFinder indicated EF1, EF2 and α-TUB were the most suitable genes. Based on these results, EF1, EF2 and α-TUB were considered as reference genes for analyzing the gene expression profiles of Cluster-3395.48660 in different infection stages, and the results were consistent with the transcriptome data. All the results suggest that the combination of EF1, EF2 and α-TUB proved to be acceptable reference genes during the interaction between G. yamadae and apple leaves.

Comparative Secretome Analysis of Magnaporthe oryzae Identified Proteins Involved in Virulence and Cell Wall Integrity

Plant fungal pathogens secrete numerous proteins into the apoplast at the plant-fungus contact sites to facilitate colonization. However, only a few secretory proteins were functionally characterized in Magnaporthe oryzae, the fungal pathogen causing rice blast disease worldwide. Asparagine-linked glycosylation 3 (Alg3) is an α-1,3-mannosyltransferase functioning in the N-glycan synthesis of N-glycosylated secretory proteins. Fungal pathogenicity and cell wall integrity are impaired in Δalg3 mutants, but the secreted proteins affected in Δalg3 mutants are largely unknown. In this study, we compared the secretomes of the wild-type strain and the Δalg3 mutant and identified 51 proteins that require Alg3 for proper secretion. These proteins were predicted to be involved in metabolic processes, interspecies interactions, cell wall organization, and response to chemicals. Nine proteins were selected for further validation. We found that these proteins were localized at the apoplastic region surrounding the fungal infection hyphae. Moreover, the N-glycosylation of these proteins was significantly changed in the Δalg3 mutant, leading to the decreased protein secretion and abnormal protein localization. Furthermore, we tested the biological functions of two genes, INV1 (encoding invertase 1, a secreted invertase) and AMCase (encoding acid mammalian chinitase, a secreted chitinase). The fungal virulence was significantly reduced, and the cell wall integrity was altered in the Δinv1 and Δamcase mutant strains. Moreover, the N-glycosylation was essential for the function and secretion of AMCase. Taken together, our study provides new insight into the role of N-glycosylated secretory proteins in fungal virulence and cell wall integrity.

Elevated Ozone Concentration and Nitrogen Addition Increase Poplar Rust Severity by Shifting the Phyllosphere Microbial Community

Tropospheric ozone and nitrogen deposition are two major environmental pollutants. A great deal of research has focused on the negative impacts of elevated O₃ and the complementary effect of soil N addition on the physiological properties of trees. However, it has been overlooked how elevated O₃ and N addition affect tree immunity in face of pathogen infection, as well as of the important roles of phyllosphere microbiome community in host–pathogen–environment interplay. Here, we examined the effects of elevated O₃ and soil N addition on poplar leaf rust [Melampsora larici-populina] severity of two susceptible hybrid poplars [clone ‘107’: Populus euramericana cv. ‘74/76’; clone ‘546’: P. deltoides Í P. cathayana] in Free-Air-Controlled-Environment plots, in addition, the link between Mlp-susceptibility and changes in microbial community was determined using Miseq amplicon sequencing. Rust severity of clone ‘107’ significantly increased under elevated O₃ or N addition only; however, the negative impact of elevated O₃ could be significantly mitigated when accompanied by N addition, likewise, this trade-off was reflected in its phyllosphere microbial α-diversity responding to elevated O₃ and N addition. However, rust severity of clone ‘546’ did not differ significantly in the cases of elevated O₃ and N addition. Mlp infection altered microbial community composition and increased its sensitivity to elevated O₃, as determined by the markedly different abundance of taxa. Elevated O₃ and N addition reduced the complexity of microbial community, which may explain the increased severity of poplar rust. These findings suggest that poplars require a changing phyllosphere microbial associations to optimize plant immunity in response to environmental changes.

Annotation survey and life-cycle transcriptomics of transcription factors in rust fungi (Pucciniales) identify a possible role for cold shock proteins in dormancy exit

Fungi of the order Pucciniales are obligate plant biotrophs causing rust diseases. They exhibit a complex life cycle with the production of up to five spore types, infection of two unrelated hosts and an overwintering stage. Transcription factors (TFs) are key regulators of gene expression in eukaryote cells. In order to better understand genetic programs expressed during major transitions of the rust life cycle, we surveyed the complement of TFs in fungal genomes with an emphasis on Pucciniales. We found that despite their large gene numbers, rust genomes have a reduced repertoire of TFs compared to other fungi. The proportions of C2H2 and Zinc cluster -two of the most represented TF families in fungi- indicate differences in their evolutionary relationships in Pucciniales and other fungal taxa. The regulatory gene family encoding cold shock protein (CSP) showed a striking expansion in Pucciniomycotina with specific duplications in the order Pucciniales. The survey of expression profiles collected by transcriptomics along the life cycle of the poplar rust fungus revealed TF genes related to major biological transitions, e.g. response to environmental cues and host infection. Particularly, poplar rust CSPs were strongly expressed in basidia produced after the overwintering stage suggesting a possible role in dormancy exit. Expression during transition from dormant telia to basidia confirmed the specific expression of the three poplar rust CSP genes. Their heterologous expression in yeast improved cell growth after cold stress exposure, suggesting a probable regulatory function when the poplar rust fungus exits dormancy. This study addresses for the first time TF and regulatory genes involved in developmental transition in the rust life cycle opening perspectives to further explore molecular regulation in the biology of the Pucciniales.

Convergence Analysis of Rust Fungi and Anther Smuts Reveals Their Common Molecular Adaptation to a Phytoparasitic Lifestyle

Convergent evolution between distantly related taxa often mirrors adaptation to similar environments. Rust fungi and anther smuts, which belong to different classes in Pucciniomycotina, have independently evolved a phytoparasitic lifestyle, representing an example of convergent evolution in the fungal kingdom. To investigate their adaptations and the genetic bases underlying their phytoparasitic lifestyles, we performed genome-wide convergence analysis of amino acid substitutions, evolutionary rates, and gene gains and losses. Convergent substitutions were detected in ATPeV0D and RP-S27Ae, two genes important for the generation of turgor pressure and ribosomal biosynthesis, respectively. A total of 51 positively selected genes were identified, including eight genes associated with translation and three genes related to the secretion pathway. In addition, rust fungi and anther smuts contained more proteins associated with oligopeptide transporters and vacuolar proteases than did other fungi. For rust fungi and anther smuts, these forms of convergence suggest four adaptive mechanisms for a phytoparasitic lifestyle: 1) reducing the metabolic demand for hyphal growth and penetration at the pre-penetration stage, 2) maintaining the efficiency of protein synthesis during colonization, 3) ensuring the normal secretion of rapidly evolving secreted proteins, and 4) improving the capacity for oligopeptide metabolism. Our results are the first to shed light on the genetic convergence mechanisms and molecular adaptation underlying phytoparasitic lifestyles in fungi.

Comparative transcriptomics and eQTL mapping of response to Melampsora americana in selected Salix purpurea F2 progeny

Background

Melampsora spp. rusts are the greatest pathogen threat to shrub willow (Salix spp.) bioenergy crops. Genetic resistance is key to limit the effects of these foliar diseases on host response and biomass yield, however, the genetic basis of host resistance has not been characterized. The addition of new genomic resources for Salix provides greater power to investigate the interaction between S. purpurea and M. americana, species commonly found in the Northeast US. Here, we utilize 3′ RNA-seq to investigate host-pathogen interactions following controlled inoculations of M. americana on resistant and susceptible F₂S. purpurea genotypes identified in a recent QTL mapping study. Differential gene expression, network analysis, and eQTL mapping were used to contrast the response to inoculation and to identify associated candidate genes.

Results

Controlled inoculation in a replicated greenhouse study identified 19 and 105 differentially expressed genes between resistant and susceptible genotypes at 42 and 66 HPI, respectively. Defense response gene networks were activated in both resistant and susceptible genotypes and enriched for many of the same defense response genes, yet the hub genes of these common response modules showed greater mean expression among the resistant plants. Further, eight and six eQTL hotspots were identified at 42 and 66 HPI, respectively. The combined results of three analyses highlight 124 candidate genes in the host for further analysis while analysis of pathogen RNA showed differential expression of 22 genes, two of which are candidate pathogen effectors.

Conclusions

We identified two differentially expressed M. americana transcripts and 124 S. purpurea genes that are good candidates for future studies to confirm their role in conferring resistance.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-021-08254-1.

Genomic signatures of a major adaptive event in the pathogenic fungus Melampsora larici-populina

The recent availability of genome-wide sequencing techniques has allowed systematic screening for molecular signatures of adaptation, including in nonmodel organisms. Host–pathogen interactions constitute good models due to the strong selective pressures that they entail. We focused on an adaptive event which affected the poplar rust fungus Melampsora larici-populina when it overcame a resistance gene borne by its host, cultivated poplar. Based on 76 virulent and avirulent isolates framing narrowly the estimated date of the adaptive event, we examined the molecular signatures of selection. Using an array of genome scan methods based on different features of nucleotide diversity, we detected a single locus exhibiting a consistent pattern suggestive of a selective sweep in virulent individuals (excess of differentiation between virulent and avirulent samples, linkage disequilibrium, genotype–phenotype statistical association, and long-range haplotypes). Our study pinpoints a single gene and further a single amino acid replacement which may have allowed the adaptive event. Although our samples are nearly contemporary to the selective sweep, it does not seem to have affected genome diversity further than the immediate vicinity of the causal locus, which can be explained by a soft selective sweep (where selection acts on standing variation) and by the impact of recombination in mitigating the impact of selection. Therefore, it seems that properties of the life cycle of M. larici-populina, which entails both high genetic diversity and outbreeding, has facilitated its adaptation.

Host Adaptation and Virulence in Heteroecious Rust Fungi

Rust fungi (Pucciniales, Basidiomycota) are obligate biotrophic pathogens that cause rust diseases in plants, inflicting severe damage to agricultural crops. Pucciniales possess the most complex life cycles known in fungi. These include an alternation of generations, the development of up to five different sporulating stages, and, for many species, the requirement of infecting two unrelated host plants during different parts of their life cycle, termed heteroecism. These fungi have been extensively studied in the past century through microscopy and inoculation studies, providing precise descriptions of their infection processes, although the molecular mechanisms underlying their unique biology are poorly understood. In this review, we cover recent genomic and life cycle transcriptomic studies in several heteroecious rust species, which provide insights into the genetic tool kits associated with host adaptation and virulence, opening new avenues for unraveling their unique evolution.

Multinuclei Occurred Under Cryopreservation and Enhanced the Pathogenicity of Melampsora larici-populina

Melampsora larici-populina is a macrocyclic rust, and the haploid stage with two nuclei and the diploid of mononuclear sequentially occur annually. During the preservation of dry urediniospores at −80°C, we found that one isolate, ΔTs₀₆, was different from the usual wild-type isolate Ts₀₆ at −20°C because it has mixed polykaryotic urediniospores. However, the other spores, including the 0, I, III, and IV stages of a life cycle, were the same as Ts₀₆. After five generations of successive inoculation and harvest of urediniospores from the compatible host Populus purdomii, the isolate ΔTs₀₆ steadily maintained more than 20% multiple nucleus spores. To test the pathogenesis variation of ΔTs₀₆, an assay of host poplars was applied to evaluate the differences between ΔTs₀₆ and Ts₀₆. After ΔTs₀₆ and Ts₀₆ inoculation, leaves of P. purdomii were used to detect the expression of small secreted proteins (SSPs) and fungal biomasses using quantitative real-time PCR (qRT-PCR) and trypan blue staining. ΔTs₀₆ displayed stronger expression of five SSPs and had a shorter latent period, a higher density of uredinia, and higher DNA mass. A transcriptomic comparison between ΔTs₀₆ and Ts₀₆ revealed that 3,224 were differentially expressed genes (DEGs), 55 of which were related to reactive oxygen species metabolism, the Mitogen-activated protein kinase (MAPK) signaling pathway, and the meiosis pathway. Ten genes in the mitotic and meiotic pathways and another two genes associated with the “response to DNA damage stimulus” all had an upward expression, which were detected by qRT-PCR in ΔTs₀₆ during cryopreservation. Gas chromatography–mass spectrometry (GC-MS) confirmed that the amounts of hexadecanoic acid and octadecadienoic acid were much more in ΔTs₀₆ than in Ts₀₆. In addition, using spectrophotometry, hydrogen peroxide (H₂O₂) was also present in greater quantities in ΔTs₀₆ compared with those found in Ts₀₆. Increased fatty acids metabolism could prevent damage to urediniospores in super-low temperatures, but oxidant species that involved H₂O₂ may destroy tube proteins of mitosis and meiosis, which could cause abnormal nuclear division and lead to multinucleation, which has a different genotype. Therefore, the multinuclear isolate is different from the wild-type isolate in terms of phenotype and genotype; this multinucleation phenomenon in urediniospores improves the pathogenesis and environmental fitness of M. larici-populina.

The Fungal Effector Mlp37347 Alters Plasmodesmata Fluxes and Enhances Susceptibility to Pathogen

Melampsora larici-populina (Mlp) is a devastating pathogen of poplar trees, causing the defoliating poplar leaf rust disease. Genomic studies have revealed that Mlp possesses a repertoire of 1184 small secreted proteins (SSPs), some of them being characterized as candidate effectors. However, how they promote virulence is still unclear. This study investigates the candidate effector Mlp37347’s role during infection. We developed a stable Arabidopsis transgenic line expressing Mlp37347 tagged with the green fluorescent protein (GFP). We found that the effector accumulated exclusively at plasmodesmata (PD). Moreover, the presence of the effector at plasmodesmata favors enhanced plasmodesmatal flux and reduced callose deposition. Transcriptome profiling and a gene ontology (GO) analysis of transgenic Arabidopsis plants expressing the effector revealed that the genes involved in glucan catabolic processes are up-regulated. This effector has previously been shown to interact with glutamate decarboxylase 1 (GAD1), and in silico docking analysis supported the strong binding between Mlp37347 and GAD1 in this study. In infection assays, the effector promoted Hyalonoperospora arabidopsidis growth but not bacterial growth. Our investigation suggests that the effector Mlp37347 targets PD in host cells and promotes parasitic growth.

Population Genomics of the Maize Pathogen Ustilago maydis: Demographic History and Role of Virulence Clusters in Adaptation

The tight interaction between pathogens and their hosts results in reciprocal selective forces that impact the genetic diversity of the interacting species. The footprints of this selection differ between pathosystems because of distinct life-history traits, demographic histories, or genome architectures. Here, we studied the genome-wide patterns of genetic diversity of 22 isolates of the causative agent of the corn smut disease, Ustilago maydis, originating from five locations in Mexico, the presumed center of origin of this species. In this species, many genes encoding secreted effector proteins reside in so-called virulence clusters in the genome, an arrangement that is so far not found in other filamentous plant pathogens. Using a combination of population genomic statistical analyses, we assessed the geographical, historical, and genome-wide variation of genetic diversity in this fungal pathogen. We report evidence of two partially admixed subpopulations that are only loosely associated with geographic origin. Using the multiple sequentially Markov coalescent model, we inferred the demographic history of the two pathogen subpopulations over the last 0.5 Myr. We show that both populations experienced a recent strong bottleneck starting around 10,000 years ago, coinciding with the assumed time of maize domestication. Although the genome average genetic diversity is low compared with other fungal pathogens, we estimated that the rate of nonsynonymous adaptive substitutions is three times higher in genes located within virulence clusters compared with nonclustered genes, including nonclustered effector genes. These results highlight the role that these singular genomic regions play in the evolution of this pathogen.

Unrelated Fungal Rust Candidate Effectors Act on Overlapping Plant Functions

Rust fungi cause epidemics that threaten the production of important plant species, such as wheat and soy. Melampsora larici-populina (Mlp) causes the poplar rust and encodes at least 1184 candidate effectors (CEs) whose functions are poorly known. In this study, we sequenced the transcriptome and used mass spectrometry to analyze the metabolome of Arabidopsis plants constitutively expressing 14 Mlp CEs and of a control line to discover alterations leading to plant susceptibility. We found 2299 deregulated genes across the experiment. Genes involved in pattern-triggered immunity, such as FRK1, PR1, RBOHD, and WRKY33, as well as AUX/IAA genes were down-regulated. We further observed that 680 metabolites were deregulated in at least one CE-expressing transgenic line, with “highly unsaturated and phenolic compounds” and “peptides” enriched among down- and up-regulated metabolites. Interestingly, transgenic lines expressing unrelated CEs had correlated patterns of gene and metabolite deregulation, while expression of CEs belonging to the same family deregulated different genes and metabolites. Thus, our results uncouple effector sequence similarity and function. This supports that effector functional investigation in the context of their virulence activity and effect on plant susceptibility requires the investigation of the individual effector and precludes generalization based on sequence similarity.

The Ustilago hordei-Barley Interaction Is a Versatile System for Characterization of Fungal Effectors

Obligate biotrophic fungal pathogens, such as Blumeria graminis and Puccinia graminis, are amongst the most devastating plant pathogens, causing dramatic yield losses in many economically important crops worldwide. However, a lack of reliable tools for the efficient genetic transformation has hampered studies into the molecular basis of their virulence or pathogenicity. In this study, we present the Ustilago hordei-barley pathosystem as a model to characterize effectors from different plant pathogenic fungi. We generate U. hordei solopathogenic strains, which form infectious filaments without the presence of a compatible mating partner. Solopathogenic strains are suitable for heterologous expression system for fungal virulence factors. A highly efficient Crispr/Cas9 gene editing system is made available for U. hordei. In addition, U. hordei infection structures during barley colonization are analyzed using transmission electron microscopy, showing that U. hordei forms intracellular infection structures sharing high similarity to haustoria formed by obligate rust and powdery mildew fungi. Thus, U. hordei has high potential as a fungal expression platform for functional studies of heterologous effector proteins in barley.

In-depth secretome analysis of Puccinia striiformis f. sp. tritici in infected wheat uncovers effector functions

The importance of wheat yellow rust disease, caused by Puccinia striiformis f. sp. tritici (Pst), has increased substantially due to the emergence of aggressive new Pst races in the last couple of decades. In an era of escalating human populations and climate change, it is vital to understand the infection mechanism of Pst in order to develop better strategies to combat wheat yellow disease. The present study focuses on the identification of small secreted proteins (SSPs) and candidate-secreted effector proteins (CSEPs) that are used by the pathogen to support infection and control disease development. We generated de novo assembled transcriptomes of Pst collected from wheat fields in central Anatolia. We inoculated both susceptible and resistant seedlings with Pst and analyzed haustoria formation. At 10 days post-inoculation (dpi), we analyzed the transcriptomes and identified 10550 Differentially Expressed Unigenes (DEGs), of which 6072 were Pst-mapped. Among those Pst-related genes, 227 were predicted as PstSSPs. In silico characterization was performed using an approach combining the transcriptomic data and data mining results to provide a reliable list to narrow down the ever-expanding repertoire of predicted effectorome. The comprehensive analysis detected 14 Differentially Expressed Small-Secreted Proteins (DESSPs) that overlapped with the genes in available literature data to serve as the best CSEPs for experimental validation. One of the CSEPs was cloned and studied to test the reliability of the presented data. Biological assays show that the randomly selected CSEP, Unigene17495 (PSTG_10917), localizes in the chloroplast and is able to suppress cell death induced by INF1 in a Nicotiana benthamiana heterologous expression system.

A Poplar Rust Effector Protein Associates with Protein Disulfide Isomerase and Enhances Plant Susceptibility

Melampsora larici-populina (Mlp), the causal agent of Populus leaf rust, secretes an array of effectors into the host through the haustorium to gain nutrients and suppress immunity. The precise mechanisms by which these effectors promote virulence remain unclear. To address this question, we developed a transgenic Arabidopsis line expressing a candidate effector, Mlp124357. Constitutive expression of the effector increased plant susceptibility to pathogens. A GxxxG motif present in Mlp124357 is required for its subcellular localization at the vacuolar membrane of the plant cell, as replacement of the glycine residues with alanines led to the delocalization of Mlp124357 to the nucleus and cytoplasm. We used immunoprecipitation and mass spectrometry (MS) to identify Mlp124357 interaction partners. Only one of the putative interaction partners knock-out line caused delocalization of the effector, indicating that Arabidopsis protein disulfide isomerase-11 (AtPDI-11) is required for the effector localization. This interaction was further confirmed by a complementation test, a yeast-two hybrid assay and a molecular modeling experiment. Moreover, localization results and infection assays suggest that AtPDI-11 act as a helper for Mlp124357. In summary, our findings established that one of Mlp effectors resides at the vacuole surface and modulates plant susceptibility.

SHR4z, a novel decoy effector from the haustorium of the parasitic weed Striga gesnerioides, suppresses host plant immunity

Cowpea (Vigna unguiculata) cultivar B301 is resistant to races SG4 and SG3 of the root parasitic weed Striga gesnerioides, developing a hypersensitive response (HR) at the site of parasite attachment. By contrast, race SG4z overcomes B301 resistance and successfully parasitises the plant. Comparative transcriptomics and in silico analysis identified a small secreted effector protein dubbed Suppressor of Host Resistance 4z (SHR4z) in the SG4z haustorium that upon transfer to the host roots causes a loss of host immunity (i.e. decreased HR and increased parasite growth). SHR4z has significant homology to the short leucine-rich repeat (LRR) domain of SOMATIC EMBRYOGENESIS RECEPTOR-LIKE KINASE (SERK) family proteins and functions by binding to VuPOB1, a host BTB-BACK domain-containing ubiquitin E3 ligase homologue, leading to its rapid turnover. VuPOB1 is shown to be a positive regulator of HR since silencing of VuPOB1 expression in transgenic B301 roots lowers the frequency of HR and increases the levels of successful SG4 parasitism and overexpression decreases parasitism by SG4z. These findings provide new insights into how parasitic weeds overcome host defences and could potentially contribute to the development of novel strategies for controlling Striga and other parasitic weeds thereby enhancing crop productivity and food security globally.

Transcriptome Analysis Identifies Plasmodiophora brassicae Secondary Infection Effector Candidates

Plasmodiophora brassicae (Wor.) is an obligate intracellular plant pathogen affecting Brassicas worldwide. Identification of effector proteins is key to understanding the interaction between P. brassicae and its susceptible host plants. To date, there is very little information available on putative effector proteins secreted by P. brassicae during a secondary infection of susceptible host plants, resulting in root gall production. A bioinformatics pipeline approach to RNA-Seq data from Arabidopsis thaliana (L.) Heynh. root tissues at 17, 20, and 24 d postinoculation (dpi) identified 32 small secreted P. brassicae proteins (SSPbPs) that were highly expressed over this secondary infection time frame. Functional signal peptides were confirmed for 31 of the SSPbPs, supporting the accuracy of the pipeline designed to identify secreted proteins. Expression profiles at 0, 2, 5, 7, 14, 21, and 28 dpi verified the involvement of some of the SSPbPs in secondary infection. For seven of the SSPbPs, a functional domain was identified using Blast2GO and 3D structure analysis and domain functionality was confirmed for SSPbP22, a kinase localized to the cytoplasm and nucleus.

Bioinformatic Detection of Positive Selection Pressure in Plant Pathogens: The Neutral Theory of Molecular Sequence Evolution in Action

The genomes of plant pathogenic fungi and oomycetes are often exposed to strong positive selection pressure. During speciation, shifts in host range and preference can lead to major adaptive changes. Furthermore, evolution of total host resistance to most isolates can force rapid evolutionary changes in host-specific pathogens. Crop pathogens are subjected to particularly intense selective pressures from monocultures and fungicides. Detection of the footprints of positive selection in plant pathogen genomes is a worthwhile endeavor as it aids understanding of the fundamental biology of these important organisms. There are two main classes of test for detection of positively selected alleles. Tests based on the ratio of non-synonymous to synonymous substitutions per site detect the footprints of multiple fixation events between divergent lineages. Thus, they are well-suited to the study of ancient adaptation events spanning speciations. On the other hand, tests that scan genomes for local fluctuations in allelic diversity within populations are suitable for detection of recent positive selection in populations. In this review, I briefly describe some of the more widely used tests of positive selection and the theory underlying them. I then discuss various examples of their application to plant pathogen genomes, emphasizing the types of genes that are associated with signatures of positive selection. I conclude with a discussion of the practicality of such tests for identification of pathogen genes of interest and the important features of pathogen ecology that must be taken into account for accurate interpretation.

Conservation and expansion of a necrosis-inducing small secreted protein family from host-variable phytopathogens of the Sclerotiniaceae

Fungal effector proteins facilitate host-plant colonization and have generally been characterized as small secreted proteins (SSPs). We classified and functionally tested SSPs from the secretomes of three closely related necrotrophic phytopathogens: Ciborinia camelliae, Botrytis cinerea, and Sclerotinia sclerotiorum. Alignment of predicted SSPs identified a large protein family that share greater than 41% amino acid identity and that have key characteristics of previously described microbe-associated molecular patterns (MAMPs). Strikingly, 73 of the 75 SSP family members were predicted within the secretome of the host-specialist C. camelliae with single-copy homologs identified in the secretomes of the host generalists S. sclerotiorum and B. cinerea. To explore the potential function of this family of SSPs, 10 of the 73 C. camelliae proteins, together with the single-copy homologs from S. sclerotiorum (SsSSP3) and B. cinerea (BcSSP2), were cloned and expressed as recombinant proteins. Infiltration of SsSSP3 and BcSSP2 into host tissue induced rapid necrosis. In contrast, only one of the 10 tested C. camelliae SSPs was able to induce a limited amount of necrosis. Analysis of chimeric proteins consisting of domains from both a necrosis-inducing and a non-necrosis-inducing SSP demonstrated that the C-terminus of the S. sclerotiorum SSP is essential for necrosis-inducing function. Deletion of the BcSSP2 homolog from B. cinerea did not affect growth or pathogenesis. Thus, this research uncovered a family of highly conserved SSPs present in diverse ascomycetes that exhibit contrasting necrosis-inducing functions.

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.

The Conifer Root and Stem Rot Pathogen (Heterobasidion parviporum): Effectome Analysis and Roles in Interspecific Fungal Interactions

Heterobasidion parviporum Niemelä & Korhonen is an economically important basidiomycete, causing root and stem rot disease of Norway spruce (Picea abies (L.) Karst) in Northern Europe. The H. parviporum genome encodes numerous small secreted proteins, which might be of importance for interacting with mycorrhiza symbionts, endophytes, and other saprotrophs. We hypothesized that small secreted proteins from H. parviporum (HpSSPs) are involved in interspecific fungal interaction. To identify HpSSP-coding genes potentially involved, we screened the H. parviporum effectome and compared their transcriptomic profiles during fungal development and in planta tree infection. We further conducted phylogenetic analysis, and identified a subset of hypothetical proteins with nonpredicted domain or unknown function as HpSSPs candidates for further characterization. The HpSSPs candidates were selected based on high-quality sequence, cysteine residue frequency, protein size, and in planta expression. We subsequently explored their roles during in vitro interaction in paired cultures of H. parviporum with ectomycorrhizal Cortinarius gentilis, endophytic Phialocephala sphaeroides, saprotrophs (Mycena sp., Phlebiopsis gigantea, and Phanerochaete chrysosporium), respectively. The transcriptomic profile revealed that a large proportion of effector candidates was either barely expressed or highly expressed under all growth conditions. In vitro dual-culture test showed that P. sphaeroides and C. gentilis were overgrown by H. parviporum. The barrage zone formation or no physical contact observed in paired cultures with the saprotrophs suggest they had either combative interaction or antibiosis effect with H. parviporum. Several HpSSPs individuals were up- or downregulated during the nonself interactions. The results of HpSSPs gene expression patterns provide additional insights into the diverse roles of SSPs in tree infection and interspecific fungal interactions.

Structural genomics applied to the rust fungus Melampsora larici-populina reveals two candidate effector proteins adopting cystine knot and NTF2-like protein folds

Rust fungi are plant pathogens that secrete an arsenal of effector proteins interfering with plant functions and promoting parasitic infection. Effectors are often species-specific, evolve rapidly, and display low sequence similarities with known proteins. How rust fungal effectors function in host cells remains elusive, and biochemical and structural approaches have been scarcely used to tackle this question. In this study, we produced recombinant proteins of eleven candidate effectors of the leaf rust fungus Melampsora larici-populina in Escherichia coli. We successfully purified and solved the three-dimensional structure of two proteins, MLP124266 and MLP124017, using NMR spectroscopy. Although both MLP124266 and MLP124017 show no sequence similarity with known proteins, they exhibit structural similarities to knottins, which are disulfide-rich small proteins characterized by intricate disulfide bridges, and to nuclear transport factor 2-like proteins, which are molecular containers involved in a wide range of functions, respectively. Interestingly, such structural folds have not been reported so far in pathogen effectors, indicating that MLP124266 and MLP124017 may bear novel functions related to pathogenicity. Our findings show that sequence-unrelated effectors can adopt folds similar to known proteins, and encourage the use of biochemical and structural approaches to functionally characterize effector candidates.

Comparative transcriptomics of Gymnosporangium spp. teliospores reveals a conserved genetic program at this specific stage of the rust fungal life cycle

BACKGROUND

Gymnosporangium spp. are fungal plant pathogens causing rust disease and most of them are known to infect two different host plants (heteroecious) with four spore stages (demicyclic). In the present study, we sequenced the transcriptome of G. japonicum teliospores on its host plant Juniperus chinensis and we performed comparison to the transcriptomes of G. yamadae and G. asiaticum at the same life stage, that happens in the same host but on different organs.

RESULTS

Functional annotation for the three Gymnosporangium species showed the expression of a conserved genetic program with the top abundant cellular categories corresponding to energy, translation and signal transduction processes, indicating that this life stage is particularly active. Moreover, the survey of predicted secretomes in the three Gymnosporangium transcriptomes revealed shared and specific genes encoding carbohydrate active enzymes and secreted proteins of unknown function that could represent candidate pathogenesis effectors. A transcript encoding a hemicellulase of the glycoside hydrolase 26 family, previously identified in other rust fungi, was particularly highly expressed suggesting a general role in rust fungi. The comparison between the transcriptomes of the three Gymnosporangium spp. and selected Pucciniales species in different taxonomical families allowed to identify lineage-specific protein families that may relate to the biology of teliospores in rust fungi. Among clustered gene families, 205, 200 and 152 proteins were specifically identified in G. japonicum, G. yamadae and G. asiaticum, respectively, including candidate effectors expressed in teliospores.

CONCLUSIONS

This comprehensive comparative transcriptomics study of three Gymnosporangium spp. identified gene functions and metabolic pathways particularly expressed in teliospores, a stage of the life cycle that is mostly overlooked in rust fungi. Secreted protein encoding transcripts expressed in teliospores may reveal new candidate effectors related to pathogenesis. Although this spore stage is not involved in host plant infection but in the production of basidiospores infecting plants in the Amygdaloideae, we speculate that candidate effectors may be expressed as early as the teliospore stage for preparing further infection by basidiospores.

The haustorial transcriptome of the cucurbit pathogen Podosphaera xanthii reveals new insights into the biotrophy and pathogenesis of powdery mildew fungi

Background

Podosphaera xanthii is the main causal agent of powdery mildew disease in cucurbits and is responsible for important yield losses in these crops worldwide. Powdery mildew fungi are obligate biotrophs. In these parasites, biotrophy is determined by the presence of haustoria, which are specialized structures of parasitism developed by these fungi for the acquisition of nutrients and the delivery of effectors. Detailed molecular studies of powdery mildew haustoria are scarce due mainly to difficulties in their isolation. Therefore, their analysis is considered an important challenge for powdery mildew research. The aim of this work was to gain insights into powdery mildew biology by analysing the haustorial transcriptome of P. xanthii.

Results

Prior to RNA isolation and massive-scale mRNA sequencing, a flow cytometric approach was developed to isolate P. xanthii haustoria free of visible contaminants. Next, several commercial kits were used to isolate total RNA and to construct the cDNA and Illumina libraries that were finally sequenced by the Illumina NextSeq system. Using this approach, the maximum amount of information from low-quality RNA that could be obtained was used to accomplish the de novo assembly of the P. xanthii haustorial transcriptome. The subsequent analysis of this transcriptome and comparison with the epiphytic transcriptome allowed us to identify the importance of several biological processes for haustorial cells such as protection against reactive oxygen species, the acquisition of different nutrients and genetic regulation mediated by non-coding RNAs. In addition, we could also identify several secreted proteins expressed exclusively in haustoria such as cell adhesion proteins that have not been related to powdery mildew biology to date.

Conclusions

This work provides a novel approach to study the molecular aspects of powdery mildew haustoria. In addition, the results of this study have also allowed us to identify certain previously unknown processes and proteins involved in the biology of powdery mildews that could be essential for their biotrophy and pathogenesis.

Electronic supplementary material

The online version of this article (10.1186/s12864-019-5938-0) contains supplementary material, which is available to authorized users.

Advances in understanding obligate biotrophy in rust fungi

Contents Summary 1190 I. Introduction 1190 II. Rust fungi: a diverse and serious threat to agriculture 1191 III. The different facets of rust life cycles and unresolved questions about their evolution 1191 IV. The biology of rust infection 1192 V. Rusts in the genomics era: the ever-expanding list of candidate effector genes 1195 VI. Functional characterization of rust effectors 1197 VII. Putting rusts to sleep: Pucciniales research outlooks 1201 Acknowledgements 1202 References 1202 SUMMARY: Rust fungi (Pucciniales) are the largest group of plant pathogens and represent one of the most devastating threats to agricultural crops worldwide. Despite the economic importance of these highly specialized pathogens, many aspects of their biology remain obscure, largely because rust fungi are obligate biotrophs. The rise of genomics and advances in high-throughput sequencing technology have presented new options for identifying candidate effector genes involved in pathogenicity mechanisms of rust fungi. Transcriptome analysis and integrated bioinformatics tools have led to the identification of key genetic determinants of host susceptibility to infection by rusts. Thousands of genes encoding secreted proteins highly expressed during host infection have been reported for different rust species, which represents significant potential towards understanding rust effector function. Recent high-throughput in planta expression screen approaches (effectoromics) have pushed the field ahead even further towards predicting high-priority effectors and identifying avirulence genes. These new insights into rust effector biology promise to inform future research and spur the development of effective and sustainable strategies for managing rust diseases.

Prediction of the Diplocarpon rosae secretome reveals candidate genes for effectors and virulence factors

Rose black spot is one of the most severe diseases of field-grown roses. Though R-genes have been characterised, little information is known about the molecular details of the interaction between pathogen and host. Based on the recently published genome sequence of the black spot fungus, we analysed gene models with various bioinformatic tools utilising the expression data of infected host tissues, which led to the prediction of 827 secreted proteins. A significant proportion of the predicted secretome comprises enzymes for the degradation of cell wall components, several of which were highly expressed during the first infection stages. As the secretome comprises major factors determining the ability of the fungus to colonise its host, we focused our further analyses on predicted effector candidates. In total, 52 sequences of 251 effector candidates matched several bioinformatic criteria of effectors, contained a Y/F/WxC motif, and did not match annotated proteins from other fungi. Additional sequences were identified based on their high expression levels during the penetration/haustorium formation phase and/or by matching known effectors from other fungi. Several host genotypes that are resistant to the sequenced isolate but differ in the R-genes responsible for this resistance are available. The combination of these genotypes with functional studies of the identified candidate effectors will allow the mechanisms of the rose black spot interaction to be dissected.

Detection of Race-Specific Resistance Against Puccinia coronata f. sp. avenae in Brachypodium Species

Oat crown rust caused by Puccinia coronata f. sp. avenae is the most destructive foliar disease of cultivated oat. Characterization of genetic factors controlling resistance responses to Puccinia coronata f. sp. avenae in nonhost species could provide new resources for developing disease protection strategies in oat. We examined symptom development and fungal colonization levels of a collection of Brachypodium distachyon and B. hybridum accessions infected with three North American P. coronata f. sp. avenae isolates. Our results demonstrated that colonization phenotypes are dependent on both host and pathogen genotypes, indicating a role for race-specific responses in these interactions. These responses were independent of the accumulation of reactive oxygen species. Expression analysis of several defense-related genes suggested that salicylic acid and ethylene-mediated signaling but not jasmonic acid are components of resistance reaction to P. coronata f. sp. avenae. Our findings provide the basis to conduct a genetic inheritance study to examine whether effector-triggered immunity contributes to nonhost resistance to P. coronata f. sp. avenae in Brachypodium spp.

A rust fungal effector binds plant DNA and modulates transcription

The basidiomycete Melampsora larici-populina causes poplar rust disease by invading leaf tissues and secreting effector proteins through specialized infection structures known as haustoria. The mechanisms by which rust effectors promote pathogen virulence are poorly understood. The present study characterized Mlp124478, a candidate effector of M. larici-populina. We used the models Arabidopsis thaliana and Nicotiana benthamiana to investigate the function of Mlp124478 in plant cells. We established that Mlp124478 accumulates in the nucleus and nucleolus, however its nucleolar accumulation is not required to promote growth of the oomycete pathogen Hyaloperonospora arabidopsidis. Stable constitutive expression of Mlp124478 in A. thaliana repressed the expression of genes involved in immune responses, and also altered leaf morphology by increasing the waviness of rosette leaves. Chip-PCR experiments showed that Mlp124478 associats'e with the TGA1a-binding DNA sequence. Our results suggest that Mlp124478 exerts a virulence activity and binds the TGA1a promoter to suppress genes induced in response to pathogen infection.

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.

Whole-Genome Resequencing and Pan-Transcriptome Reconstruction Highlight the Impact of Genomic Structural Variation on Secondary Metabolite Gene Clusters in the Grapevine Esca Pathogen Phaeoacremonium minimum

The Ascomycete fungus Phaeoacremonium minimum is one of the primary causal agents of Esca, a widespread and damaging grapevine trunk disease. Variation in virulence among Pm. minimum isolates has been reported, but the underlying genetic basis of the phenotypic variability remains unknown. The goal of this study was to characterize intraspecific genetic diversity and explore its potential impact on virulence functions associated with secondary metabolism, cellular transport, and cell wall decomposition. We generated a chromosome-scale genome assembly, using single molecule real-time sequencing, and resequenced the genomes and transcriptomes of multiple isolates to identify sequence and structural polymorphisms. Numerous insertion and deletion events were found for a total of about 1 Mbp in each isolate. Structural variation in this extremely gene dense genome frequently caused presence/absence polymorphisms of multiple adjacent genes, mostly belonging to biosynthetic clusters associated with secondary metabolism. Because of the observed intraspecific diversity in gene content due to structural variation we concluded that a transcriptome reference developed from a single isolate is insufficient to represent the virulence factor repertoire of the species. We therefore compiled a pan-transcriptome reference of Pm. minimum comprising a non-redundant set of 15,245 protein-coding sequences. Using naturally infected field samples expressing Esca symptoms, we demonstrated that mapping of meta-transcriptomics data on a multi-species reference that included the Pm. minimum pan-transcriptome allows the profiling of an expanded set of virulence factors, including variable genes associated with secondary metabolism and cellular transport.

Genome sequence and effectorome of Moniliophthora perniciosa and Moniliophthora roreri subpopulations

Background

The hemibiotrophic pathogens Moniliophthora perniciosa (witches’ broom disease) and Moniliophthora roreri (frosty pod rot disease) are among the most important pathogens of cacao. Moniliophthora perniciosa has a broad host range and infects a variety of meristematic tissues in cacao plants, whereas M. roreri infects only pods of Theobroma and Herrania genera. Comparative pathogenomics of these fungi is essential to understand Moniliophthora infection strategies, therefore the detection and in silico functional characterization of effector candidates are important steps to gain insight on their pathogenicity.

Results

Candidate secreted effector proteins repertoire were predicted using the genomes of five representative isolates of M. perniciosa subpopulations (three from cacao and two from solanaceous hosts), and one representative isolate of M. roreri from Peru. Many putative effectors candidates were identified in M. perniciosa: 157 and 134 in cacao isolates from Bahia, Brazil; 109 in cacao isolate from Ecuador, 92 and 80 in wild solanaceous isolates from Minas Gerais (Lobeira) and Bahia (Caiçara), Brazil; respectively. Moniliophthora roreri showed the highest number of effector candidates, a total of 243. A set of eight core effectors were shared among all Moniliophthora isolates, while others were shared either between the wild solanaceous isolates or among cacao isolates. Mostly, candidate effectors of M. perniciosa were shared among the isolates, whereas in M. roreri nearly 50% were exclusive to the specie. In addition, a large number of cell wall-degrading enzymes characteristic of hemibiotrophic fungi were found. From these, we highlighted the proteins involved in cell wall modification, an enzymatic arsenal that allows the plant pathogens to inhabit environments with oxidative stress, which promotes degradation of plant compounds and facilitates infection.

Conclusions

The present work reports six genomes and provides a database of the putative effectorome of Moniliophthora, a first step towards the understanding of the functional basis of fungal pathogenicity.

Electronic supplementary material

The online version of this article (10.1186/s12864-018-4875-7) contains supplementary material, which is available to authorized users.

The Rust Fungus Melampsora larici-populina Expresses a Conserved Genetic Program and Distinct Sets of Secreted Protein Genes During Infection of Its Two Host Plants, Larch and Poplar

Mechanisms required for broad-spectrum or specific host colonization of plant parasites are poorly understood. As a perfect illustration, heteroecious rust fungi require two alternate host plants to complete their life cycles. Melampsora larici-populina infects two taxonomically unrelated plants, larch, on which sexual reproduction is achieved, and poplar, on which clonal multiplication occurs, leading to severe epidemics in plantations. We applied deep RNA sequencing to three key developmental stages of M. larici-populina infection on larch: basidia, pycnia, and aecia, and we performed comparative transcriptomics of infection on poplar and larch hosts, using available expression data. Secreted protein was the only significantly overrepresented category among differentially expressed M. larici-populina genes between the basidial, the pycnial, and the aecial stages, highlighting their probable involvement in the infection process. Comparison of fungal transcriptomes in larch and poplar revealed a majority of rust genes were commonly expressed on the two hosts and a fraction exhibited host-specific expression. More particularly, gene families encoding small secreted proteins presented striking expression profiles that highlight probable candidate effectors specialized on each host. Our results bring valuable new information about the biological cycle of rust fungi and identify genes that may contribute to host specificity.

Crystal structure of the Melampsora lini effector AvrP reveals insights into a possible nuclear function and recognition by the flax disease resistance protein P

The effector protein AvrP is secreted by the flax rust fungal pathogen (Melampsora lini) and recognized specifically by the flax (Linum usitatissimum) P disease resistance protein, leading to effector-triggered immunity. To investigate the biological function of this effector and the mechanisms of specific recognition by the P resistance protein, we determined the crystal structure of AvrP. The structure reveals an elongated zinc-finger-like structure with a novel interleaved zinc-binding topology. The residues responsible for zinc binding are conserved in AvrP effector variants and mutations of these motifs result in a loss of P-mediated recognition. The first zinc-coordinating region of the structure displays a positively charged surface and shows some limited similarities to nucleic acid-binding and chromatin-associated proteins. We show that the majority of the AvrP protein accumulates in the plant nucleus when transiently expressed in Nicotiana benthamiana cells, suggesting a nuclear pathogenic function. Polymorphic residues in AvrP and its allelic variants map to the protein surface and could be associated with differences in recognition specificity. Several point mutations of residues on the non-conserved surface patch result in a loss of recognition by P, suggesting that these residues are required for recognition.

The Hydrophobin-Like OmSSP1 May Be an Effector in the Ericoid Mycorrhizal Symbiosis

Mutualistic and pathogenic plant-colonizing fungi use effector molecules to manipulate the host cell metabolism to allow plant tissue invasion. Some small secreted proteins (SSPs) have been identified as fungal effectors in both ectomycorrhizal and arbuscular mycorrhizal fungi, but it is currently unknown whether SSPs also play a role as effectors in other mycorrhizal associations. Ericoid mycorrhiza is a specific endomycorrhizal type that involves symbiotic fungi mostly belonging to the Leotiomycetes (Ascomycetes) and plants in the family Ericaceae. Genomic and RNASeq data from the ericoid mycorrhizal fungus Oidiodendron maius led to the identification of several symbiosis-upregulated genes encoding putative SSPs. OmSSP1, the most highly symbiosis up-regulated SSP, was found to share some features with fungal hydrophobins, even though it lacks the Pfam hydrophobin domain. Sequence alignment with other hydrophobins and hydrophobin-like fungal proteins placed OmSSP1 within Class I hydrophobins. However, the predicted features of OmSSP1 may suggest a distinct type of hydrophobin-like proteins. The presence of a predicted signal peptide and a yeast-based signal sequence trap assay demonstrate that OmSSP1 is secreted. OmSSP1 null-mutants showed a reduced capacity to form ericoid mycorrhiza with Vaccinium myrtillus roots, suggesting a role as effectors in the ericoid mycorrhizal interaction.

Comparative systems analysis of the secretome of the opportunistic pathogen Aspergillus fumigatus and other Aspergillus species

Aspergillus fumigatus and multiple other Aspergillus species cause a wide range of lung infections, collectively termed aspergillosis. Aspergilli are ubiquitous in environment with healthy immune systems routinely eliminating inhaled conidia, however, Aspergilli can become an opportunistic pathogen in immune-compromised patients. The aspergillosis mortality rate and emergence of drug-resistance reveals an urgent need to identify novel targets. Secreted and cell membrane proteins play a critical role in fungal-host interactions and pathogenesis. Using a computational pipeline integrating data from high-throughput experiments and bioinformatic predictions, we have identified secreted and cell membrane proteins in ten Aspergillus species known to cause aspergillosis. Small secreted and effector-like proteins similar to agents of fungal-plant pathogenesis were also identified within each secretome. A comparison with humans revealed that at least 70% of Aspergillus secretomes have no sequence similarity with the human proteome. An analysis of antigenic qualities of Aspergillus proteins revealed that the secretome is significantly more antigenic than cell membrane proteins or the complete proteome. Finally, overlaying an expression dataset, four A. fumigatus proteins upregulated during infection and with available structures, were found to be structurally similar to known drug target proteins in other organisms, and were able to dock in silico with the respective drug.

Show me the way: rust effector targets in heterologous plant systems

For years, the study of rust fungal effectors has been impeded by the lack of molecular genetic tools in rust pathosystems. The recent use of heterologous plants to perform effector screens (effectoromics)-including effector localisation (cellular targets) and protein interactors (molecular targets) in plant cells-has changed the game. These screens revealed that many candidate effectors from various rust fungi target specific plant cell compartments, including chloroplasts, and associate with specific plant protein complexes. Such information represents unparalleled opportunities to understand how effectors sustain extreme parasitic interactions and obligate biotrophy. Despite their limitations, we here portray how the use of heterologous expression systems has been essential for gaining new insight into rust effectors.

Positively Selected Effector Genes and Their Contribution to Virulence in the Smut Fungus Sporisorium reilianum

Plants and fungi display a broad range of interactions in natural and agricultural ecosystems ranging from symbiosis to parasitism. These ecological interactions result in coevolution between genes belonging to different partners. A well-understood example is secreted fungal effector proteins and their host targets, which play an important role in pathogenic interactions. Biotrophic smut fungi (Basidiomycota) are well-suited to investigate the evolution of plant pathogens, because several reference genomes and genetic tools are available for these species. Here, we used the genomes of Sporisorium reilianum f. sp. zeae and S. reilianum f. sp. reilianum, two closely related formae speciales infecting maize and sorghum, respectively, together with the genomes of Ustilago hordei, Ustilago maydis, and Sporisorium scitamineum to identify and characterize genes displaying signatures of positive selection. We identified 154 gene families having undergone positive selection during species divergence in at least one lineage, among which 77% were identified in the two investigated formae speciales of S. reilianum. Remarkably, only 29% of positively selected genes encode predicted secreted proteins. We assessed the contribution to virulence of nine of these candidate effector genes in S. reilianum f. sp. zeae by deleting individual genes, including a homologue of the effector gene pit2 previously characterized in U. maydis. Only the pit2 deletion mutant was found to be strongly reduced in virulence. Additional experiments are required to understand the molecular mechanisms underlying the selection forces acting on the other candidate effector genes, as well as the large fraction of positively selected genes encoding predicted cytoplasmic proteins.

A root-knot nematode small glycine and cysteine-rich secreted effector, MiSGCR1, is involved in plant parasitism

Root-knot nematodes, Meloidogyne spp., are obligate endoparasites that maintain a biotrophic relationship with their hosts. They infect roots as microscopic vermiform second-stage juveniles, and establish specialized feeding structures called 'giant-cells', from which they withdraw water and nutrients. The nematode effector proteins secreted in planta are key elements in the molecular dialogue of parasitism. Here, we compared Illumina RNA-seq transcriptomes for M. incognita obtained at various points in the lifecycle, and identified 31 genes more strongly expressed in parasitic stages than in preparasitic juveniles. We then selected candidate effectors for functional characterization. Quantitative real-time PCR and in situ hybridizations showed that the validated differentially expressed genes are predominantly specifically expressed in oesophageal glands of the nematode. We also soaked the nematodes in siRNA to silence these genes and to determine their role in pathogenicity. The silencing of the dorsal gland specific-Minc18876 and its paralogues resulted in a significant, reproducible decrease in the number of mature females with egg masses, demonstrating a potentially important role for the small glycine- and cysteine-rich effector MiSGCR1 in early stages of plant-nematode interaction. Finally, we report that MiSGCR1 suppresses plant cell death induced by bacterial or oomycete triggers of plant defense.

Identification of Plasmodiophora brassicae effectors - A challenging goal

Clubroot is an economically important disease affecting Brassica plants worldwide. Plasmodiophora brassicae is the protist pathogen associated with the disease, and its soil-borne obligate parasitic nature has impeded studies related to its biology and the mechanisms involved in its infection of the plant host. The identification of effector proteins is key to understanding how the pathogen manipulates the plant's immune response and the genes involved in resistance. After more than 140 years studying clubroot and P. brassicae, very little is known about the effectors playing key roles in the infection process and subsequent disease progression. Here we analyze the information available for identified effectors and suggest several features of effector genes that can be used in the search for others. Based on the information presented in this review, we propose a comprehensive bioinformatics pipeline for effector identification and provide a list of the bioinformatics tools available for such.

Non-parent of Origin Expression of Numerous Effector Genes Indicates a Role of Gene Regulation in Host Adaption of the Hybrid Triticale Powdery Mildew Pathogen

Powdery mildew is an important disease of cereals. It is caused by one species, Blumeria graminis, which is divided into formae speciales each of which is highly specialized to one host. Recently, a new form capable of growing on triticale (B.g. triticale) has emerged through hybridization between wheat and rye mildews (B.g. tritici and B.g. secalis, respectively). In this work, we used RNA sequencing to study the molecular basis of host adaptation in B.g. triticale. We analyzed gene expression in three B.g. tritici isolates, two B.g. secalis isolates and two B.g. triticale isolates and identified a core set of putative effector genes that are highly expressed in all formae speciales. We also found that the genes differentially expressed between isolates of the same form as well as between different formae speciales were enriched in putative effectors. Their coding genes belong to several families including some which contain known members of mildew avirulence (Avr) and suppressor (Svr) genes. Based on these findings we propose that effectors play an important role in host adaptation that is mechanistically based on Avr-Resistance gene-Svr interactions. We also found that gene expression in the B.g. triticale hybrid is mostly conserved with the parent-of-origin, but some genes inherited from B.g. tritici showed a B.g. secalis-like expression. Finally, we identified 11 unambiguous cases of putative effector genes with hybrid-specific, non-parent of origin gene expression, and we propose that they are possible determinants of host specialization in triticale mildew. These data suggest that altered expression of multiple effector genes, in particular Avr and Svr related factors, might play a role in mildew host adaptation based on hybridization.

Infection assays in Arabidopsis reveal candidate effectors from the poplar rust fungus that promote susceptibility to bacteria and oomycete pathogens

Fungi of the Pucciniales order cause rust diseases which, altogether, affect thousands of plant species worldwide and pose a major threat to several crops. How rust effectors-virulence proteins delivered into infected tissues to modulate host functions-contribute to pathogen virulence remains poorly understood. Melampsora larici-populina is a devastating and widespread rust pathogen of poplar, and its genome encodes 1184 identified small secreted proteins that could potentially act as effectors. Here, following specific criteria, we selected 16 candidate effector proteins and characterized their virulence activities and subcellular localizations in the leaf cells of Arabidopsis thaliana. Infection assays using bacterial (Pseudomonas syringae) and oomycete (Hyaloperonospora arabidopsidis) pathogens revealed subsets of candidate effectors that enhanced or decreased pathogen leaf colonization. Confocal imaging of green fluorescent protein-tagged candidate effectors constitutively expressed in stable transgenic plants revealed that some protein fusions specifically accumulate in nuclei, chloroplasts, plasmodesmata and punctate cytosolic structures. Altogether, our analysis suggests that rust fungal candidate effectors target distinct cellular components in host cells to promote parasitic growth.

Different waves of effector genes with contrasted genomic location are expressed by Leptosphaeria maculans during cotyledon and stem colonization of oilseed rape

Leptosphaeria maculans, the causal agent of stem canker disease, colonizes oilseed rape (Brassica napus) in two stages: a short and early colonization stage corresponding to cotyledon or leaf colonization, and a late colonization stage during which the fungus colonizes systemically and symptomlessly the plant during several months before stem canker appears. To date, the determinants of the late colonization stage are poorly understood; L. maculans may either successfully escape plant defences, leading to stem canker development, or the plant may develop an 'adult-stage' resistance reducing canker incidence. To obtain an insight into these determinants, we performed an RNA-sequencing (RNA-seq) pilot project comparing fungal gene expression in infected cotyledons and in symptomless or necrotic stems. Despite the low fraction of fungal material in infected stems, sufficient fungal transcripts were detected and a large number of fungal genes were expressed, thus validating the feasibility of the approach. Our analysis showed that all avirulence genes previously identified are under-expressed during stem colonization compared with cotyledon colonization. A validation RNA-seq experiment was then performed to investigate the expression of candidate effector genes during systemic colonization. Three hundred and seven 'late' effector candidates, under-expressed in the early colonization stage and over-expressed in the infected stems, were identified. Finally, our analysis revealed a link between the regulation of expression of effectors and their genomic location: the 'late' effector candidates, putatively involved in systemic colonization, are located in gene-rich genomic regions, whereas the 'early' effector genes, over-expressed in the early colonization stage, are located in gene-poor regions of the genome.