Diversifying selection in the wheat stem rust fungus acts predominantly on pathogen-associated gene families and reveals candidate effectors

Genome comparisons reveal accessory genes crucial for the evolution of apple Glomerella leaf spot pathogenicity in Colletotrichum fungi

Apple Glomerella leaf spot (GLS) is an emerging fungal disease caused by Colletotrichum fructicola and other Colletotrichum species. These species are polyphyletic and it is currently unknown how these pathogens convergently evolved to infect apple. We generated chromosome-level genome assemblies of a GLS-adapted isolate and a non-adapted isolate in C. fructicola using long-read sequencing. Additionally, we resequenced 17 C. fructicola and C. aenigma isolates varying in GLS pathogenicity using short-read sequencing. Genome comparisons revealed a conserved bipartite genome architecture involving minichromosomes (accessory chromosomes) shared by C. fructicola and other closely related species within the C. gloeosporioides species complex. Moreover, two repeat-rich genomic regions (1.61 Mb in total) were specifically conserved among GLS-pathogenic isolates in C. fructicola and C. aenigma. Single-gene deletion of 10 accessory genes within the GLS-specific regions of C. fructicola identified three that were essential for GLS pathogenicity. These genes encoded a putative non-ribosomal peptide synthetase, a flavin-binding monooxygenase and a small protein with unknown function. These results highlight the crucial role accessory genes play in the evolution of Colletotrichum pathogenicity and imply the significance of an unidentified secondary metabolite in GLS pathogenesis.

Soybean-Phakopsora pachyrhizi interactions: towards the development of next-generation disease-resistant plants

Soybean rust (SBR), caused by the obligate biotrophic fungus Phakopsora pachyrhizi, is a devastating foliar disease threatening soybean production. To date, no commercial cultivars conferring durable resistance to SBR are available. The development of long-lasting SBR resistance has been hindered by the lack of understanding of this complex pathosystem, encompassing challenges posed by intricate genetic structures in both the host and pathogen, leading to a gap in the knowledge of gene-for-gene interactions between soybean and P. pachyrhizi. In this review, we focus on recent advancements and emerging technologies that can be used to improve our understanding of the P. pachyrhizi-soybean molecular interactions. We further explore approaches used to combat SBR, including conventional breeding, transgenic approaches and RNA interference, and how advances in our understanding of plant immune networks, the availability of new molecular tools, and the recent sequencing of the P. pachyrhizi genome could be used to aid in the development of better genetic resistance against SBR. Lastly, we discuss the research gaps of this pathosystem and how new technologies can be used to shed light on these questions and to develop durable next-generation SBR-resistant soybean plants.

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.

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.

Novel Genetic Dysregulations and Oxidative Damage in Fusarium graminearum Induced by Plant Defense Eliciting Psychrophilic Bacillus atrophaeus TS1

This study elaborates inter-kingdom signaling mechanisms, presenting a sustainable and eco-friendly approach to combat biotic as well as abiotic stress in wheat. Fusarium graminearum is a devastating pathogen causing head and seedling blight in wheat, leading to huge yield and economic losses. Psychrophilic Bacillus atrophaeus strain TS1 was found as a potential biocontrol agent for suppression of F. graminearum under low temperature by carrying out extensive biochemical and molecular studies in comparison with a temperate biocontrol model strain Bacillus amyloliquefaciens FZB42 at 15 and 25 °C. TS1 was able to produce hydrolytic extracellular enzymes as well as antimicrobial lipopeptides, i.e., surfactin, bacillomycin, and fengycin, efficiently at low temperatures. The Bacillus strain-induced oxidative cellular damage, ultrastructural deformities, and novel genetic dysregulations in the fungal pathogen as the bacterial treatment at low temperature were able to downregulate the expression of newly predicted novel fungal genes potentially belonging to necrosis inducing protein families (fgHCE and fgNPP1). The wheat pot experiments conducted at 15 and 25 °C revealed the potential of TS1 to elicit sudden induction of plant defense, namely, H₂O₂ and callose enhanced activity of plant defense-related enzymes and induced over-expression of defense-related genes which accumulatively lead to the suppression of F. graminearum and decreased diseased leaf area.

Paediatric Antimicrobial Stewardship for Respiratory Infections in the Emergency Setting: A Systematic Review

Antimicrobial resistance occurs due to the propensity of microbial pathogens to develop resistance to antibiotics over time. Antimicrobial stewardship programs (ASPs) have been developed in response to this growing crisis, to limit unnecessary antibiotic prescription through initiatives such as education-based seminars, prescribing guidelines, and rapid respiratory pathogen (RRP) testing. Paediatric patients who present to the emergency setting with respiratory symptoms are a particularly high-risk population susceptible to inappropriate antibiotic prescribing behaviours and are therefore an ideal cohort for focused ASPs. The purpose of this systematic review was to assess the efficacy and safety of ASPs in this clinical context. A systematic search of PubMed, Medline, EMBASE and the Cochrane Database of Systematic Reviews was conducted to review the current evidence. Thirteen studies were included in the review and these studies assessed a range of stewardship interventions and outcome measures. Overall, ASPs reduced the rates of antibiotic prescription, increased the prescription of narrow-spectrum antibiotics, and shortened the duration of antibiotic therapy. Multimodal interventions that were education-based and those that used RRP testing were found to be the most effective. Whilst we found strong evidence that ASPs are effective in reducing antibiotic prescribing, further studies are required to assess whether they translate to equivalent clinical outcomes.

Genome-Wide Analysis of Four Pathotypes of Wheat Rust Pathogen (Puccinia graminis) Reveals Structural Variations and Diversifying Selection

Diseases caused by Puccinia graminis are some of the most devastating diseases of wheat. Extensive genomic understanding of the pathogen has proven helpful not only in understanding host- pathogen interaction but also in finding appropriate control measures. In the present study, whole-genome sequencing of four diverse P. graminis pathotypes was performed to understand the genetic variation and evolution. An average of 63.5 Gb of data per pathotype with about 100× average genomic coverage was achieved with 100-base paired-end sequencing performed with Illumina Hiseq 1000. Genome structural annotations collectively predicted 9273 functional proteins including ~583 extracellular secreted proteins. Approximately 7.4% of the genes showed similarity with the PHI database which is suggestive of their significance in pathogenesis. Genome-wide analysis demonstrated pathotype 117-6 as likely distinct and descended through a different lineage. The 3-6% more SNPs in the regulatory regions and 154 genes under positive selection with their orthologs and under negative selection in the other three pathotypes further supported pathotype 117-6 to be highly diverse in nature. The genomic information generated in the present study could serve as an important source for comparative genomic studies across the genus Puccinia and lead to better rust management in wheat.

Positive Selection of Transcription Factors Is a Prominent Feature of the Evolution of a Plant Pathogenic Genus Originating in the Miocene

Tests based on the d_N/d_S statistic are used to identify positive selection of nonsynonymous polymorphisms. Using these tests on alignments of all orthologs from related species can provide insights into which gene categories have been most frequently positively selected. However, longer alignments have more power to detect positive selection, creating a detection bias that could create misleading results from functional enrichment tests. Most studies of positive selection in plant pathogens focus on genes with specific virulence functions, with little emphasis on broader molecular processes. Furthermore, no studies in plant pathogens have accounted for detection bias due to alignment length when performing functional enrichment tests. To address these research gaps, we analyze 12 genomes of the phytopathogenic fungal genus Botrytis, including two sequenced in this study. To establish a temporal context, we estimated fossil-calibrated divergence times for the genus. We find that Botrytis likely originated 16–18 Ma in the Miocene and underwent continuous radiation ending in the Pliocene. An untargeted scan of Botrytis single-copy orthologs for positive selection with three different statistical tests uncovered evidence for positive selection among proteases, signaling proteins, CAZymes, and secreted proteins. There was also a strong overrepresentation of transcription factors among positively selected genes. This overrepresentation was still apparent after two complementary controls for detection bias due to sequence length. Positively selected sites were depleted within DNA-binding domains, suggesting changes in transcriptional responses to internal and external cues or protein–protein interactions have undergone positive selection more frequently than changes in promoter fidelity.

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.

Plant-pathogen interactions: MicroRNA-mediated trans-kingdom gene regulation in fungi and their host plants

MicroRNAs (miRNAs) have been prevalently studied in plants, animals, and viruses. However, recent studies show evidences of miRNA-like RNAs (milRNAs) in fungi as well. It is known that after successful infection, pathogens hijack the host machinery and use it for their own growth and multiplication. Alternatively, resistant plants can overcome the pathogen attack by a variety of mechanisms. Based on this prior knowledge, we computationally predicted milRNAs from 13 fungi, and identified their targets in transcriptomes of the respective fungi as well as their host plants. The expressions of the milRNAs and targets were confirmed using qRT-PCR. We found that plant miRNAs targeted fungal virulence genes, while fungal milRNAs targeted plant resistance genes; corroborating miRNA-mediated trans-kingdom gene regulation and the roles of miRNAs in plant-pathogen interactions. Transgenic plants with miRNAs targeting fungal virulence genes, or anti-sense of fungal milRNAs, would be expected to be highly resistant to the fungal pathogens.

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.

Label free proteomics and systematic analysis of secretome reveals effector candidates regulated by SGE1 and FTF1 in the plant pathogen Fusarium oxysporum f. sp. cubense tropical race 4

BACKGROUND

Phytopathogens secreted effectors during host colonization to suppress or trigger plant immunity. Identification of new effectors is one of the research focuses in recent years. There is only a limited knowledge about effectors of Fusarium oxysporum f. sp. Cubense tropical race 4 (Foc TR4), the causal agent of wilt disease in Cavendish banana.

RESULTS

Two transcription factors, SGE1 and FTF1, were constitutively over-expressed in Foc TR4 to partially mimic the in-planta state. Secreted proteins with high purity were prepared through a two-round extraction method. Then the secretome were analyzed via label free proteomics method. A total of 919 non-redundant proteins were detected, of which 74 proteins were predicted to be effector candidates. Among these candidates, 29 were up-regulated and 13 down-regulated in the strain over-expressing SGE1 and FTF1, 8 were up-regulated and 4 down-regulated in either SGE1 or FTF1 over expression strain.

CONCLUSIONS

Through label free proteomics analysis, a series of effector candidates were identified in secretome of Foc TR4. Our work put a foundation for functional research of these effectors.

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.

DsEcp2-1 is a polymorphic effector that restricts growth of Dothistroma septosporum in pine

The detrimental effect of fungal pathogens on forest trees is an increasingly important problem that has implications for the health of our planet. Despite this, the study of molecular plant-microbe interactions in forest trees is in its infancy, and very little is known about the roles of effector molecules from forest pathogens. Dothistroma septosporum causes a devastating needle blight disease of pines, and intriguingly, is closely related to Cladosporium fulvum, a tomato pathogen in which pioneering effector biology studies have been carried out. Here, we studied D. septosporum effectors that are shared with C. fulvum, by comparing gene sequences from global isolates of D. septosporum and assessing effector function in both host and non-host plants. Many of the effectors were predicted to be non-functional in D. septosporum due to their pseudogenization or low expression in planta, suggesting adaptation to lifestyle and host. Effector sequences were polymorphic among a global collection of D. septosporum isolates, but there was no evidence for positive selection. The DsEcp2-1 effector elicited cell death in the non-host plant Nicotiana tabacum, whilst D. septosporum DsEcp2-1 mutants showed increased colonization of pine needles. Together these results suggest that DsEcp2-1 might be recognized by an immune receptor in both angiosperm and gymnosperm plants. This work may lead to the identification of plant targets for DsEcp2-1 that will provide much needed information on the molecular basis of gymnosperm-pathogen interactions in forests, and may also lead to novel methods of disease control.

Identification and Functional Characterization of an Effector Secreted by Cronartium ribicola

Cri-9402 was identified as a protein effector from Cronartium ribicola, based on the effect of its expression on growth of Pseudomonas syringae Psm ES4326 introduced into transiently transformed tobacco leaves and stably transformed Arabidopsis seedlings. In tobacco leaves transiently expressing its coding sequence, growth of P. syringae Psm ES4326 was inhibited. Expression of pathogenesis-related (PR) protein 2 (PR2), PR4a, endochitinase B, hypersensitive-related 201 (HSR201), HSR203J, and proteinase inhibitor 1 was upregulated but expression of PR1, coronatine insensitive 1, and abscisic acid 1 was significantly suppressed. In transformed Arabidopsis seedlings, the effector stimulated growth of P. syringae Psm ES4326; significantly suppressed expression of PR1, PR2, nonexpresser of pathogenesis-related genes 1 (NPR1), NPR3, NPR4, phytoalexin deficient 4, and salicylic acid induction deficient 2; and enhanced expression of plant defensin 1.2 (PDF1.2). The above results showed that the majority of responses to this effector in tobacco leaves were converse to those in transformed Arabidopsis. We could conclude that Cri-9402 promoted disease resistance in tobacco leaves and disease susceptibility in Arabidopsis seedlings. Its transcript was mainly expressed in aeciospores of C. ribicola and was probably involved in production or germination of aeciospores, and it was an effector potentially functioning in white pine-blister rust interactions.

The emergence of the multi-species NIP1 effector in Rhynchosporium was accompanied by high rates of gene duplications and losses

Plant pathogens secrete effector proteins to manipulate the host and facilitate infection. Cognate hosts trigger strong defence responses upon detection of these effectors. Consequently, pathogens and hosts undergo rapid coevolutionary arms races driven by adaptive evolution of effectors and receptors. Because of their high rate of turnover, most effectors are thought to be species-specific and the evolutionary trajectories are poorly understood. Here, we investigate the necrosis-inducing protein 1 (NIP1) effector in the multihost pathogen genus Rhynchosporium. We retraced the evolutionary history of the NIP1 locus using whole-genome assemblies of 146 strains covering four closely related species. NIP1 orthologues were present in all species but the locus consistently segregated presence-absence polymorphisms suggesting long-term balancing selection. We also identified previously unknown paralogues of NIP1 that were shared among multiple species and showed substantial copy-number variation within R. commune. The NIP1A paralogue was under significant positive selection suggesting that NIP1A is the dominant effector variant coevolving with host immune receptors. Consistent with this prediction, we found that copy number variation at NIP1A had a stronger effect on virulence than NIP1B. Our analyses unravelled the origins and diversification mechanisms of a pathogen effector family shedding light on how pathogens gain adaptive genetic variation.

Mining the effector repertoire of the biotrophic fungal pathogen Ustilago hordei during host and non-host infection

The success of plant-pathogenic fungi mostly relies on their arsenal of virulence factors which are expressed and delivered into the host tissue during colonization. The biotrophic fungal pathogen Ustilago hordei causes covered smut disease on both barley and oat. In this study, we combined cytological, genomics and molecular biological methods to achieve a better understanding of the molecular interactions in the U. hordei-barley pathosystem. Microscopic analysis revealed that U. hordei densely colonizes barley leaves on penetration, in particular the vascular system. Transcriptome analysis of U. hordei at different stages of host infection revealed differential expression of the transcript levels of 273 effector gene candidates. Furthermore, U. hordei transcriptionally activates core effector genes which may suppress even non-host early defence responses. Based on expression profiles and novelty of sequences, knockout studies of 14 effector candidates were performed in U. hordei, which resulted in the identification of four virulence factors required for host colonization. Yeast two-hybrid screening identified potential barley targets for two of the effectors. Overall, this study provides a first systematic analysis of the effector repertoire of U. hordei and identifies four effectors (Uvi1-Uvi4) as virulence factors for the infection of barley.

Evolutionarily conserved partial gene duplication in the Triticeae tribe of grasses confers pathogen resistance

Background

The large and highly repetitive genomes of the cultivated species Hordeum vulgare (barley), Triticum aestivum (wheat), and Secale cereale (rye) belonging to the Triticeae tribe of grasses appear to be particularly rich in gene-like sequences including partial duplicates. Most of them have been classified as putative pseudogenes. In this study we employ transient and stable gene silencing- and over-expression systems in barley to study the function of HvARM1 (for H. vulgare Armadillo 1), a partial gene duplicate of the U-box/armadillo-repeat E3 ligase HvPUB15 (for H. vulgare Plant U-Box 15).

Results

The partial ARM1 gene is derived from a gene-duplication event in a common ancestor of the Triticeae and contributes to quantitative host as well as nonhost resistance to the biotrophic powdery mildew fungus Blumeria graminis. In barley, allelic variants of HvARM1 but not of HvPUB15 are significantly associated with levels of powdery mildew infection. Both HvPUB15 and HvARM1 proteins interact in yeast and plant cells with the susceptibility-related, plastid-localized barley homologs of THF1 (for Thylakoid formation 1) and of ClpS1 (for Clp-protease adaptor S1) of Arabidopsis thaliana. A genome-wide scan for partial gene duplicates reveals further events in barley resulting in stress-regulated, potentially neo-functionalized, genes.

Conclusion

The results suggest neo-functionalization of the partial gene copy HvARM1 increases resistance against powdery mildew infection. It further links plastid function with susceptibility to biotrophic pathogen attack. These findings shed new light on a novel mechanism to employ partial duplication of protein-protein interaction domains to facilitate the expansion of immune signaling networks.

Electronic supplementary material

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

Restriction associated DNA-genotyping at multiple spatial scales in Arabidopsis lyrata reveals signatures of pathogen-mediated selection

Background

Genome scans based on outlier analyses have revolutionized detection of genes involved in adaptive processes, but reports of some forms of selection, such as balancing selection, are still limited. It is unclear whether high throughput genotyping approaches for identification of single nucleotide polymorphisms have sufficient power to detect modes of selection expected to result in reduced genetic differentiation among populations. In this study, we used Arabidopsis lyrata to investigate whether signatures of balancing selection can be detected based on genomic smoothing of Restriction Associated DNA sequencing (RAD-seq) data. We compared how different sampling approaches (both within and between subspecies) and different background levels of polymorphism (inbreeding or outcrossing populations) affected the ability to detect genomic regions showing key signatures of balancing selection, specifically elevated polymorphism, reduced differentiation and shifts towards intermediate allele frequencies. We then tested whether candidate genes associated with disease resistance (R-gene analogs) were detected more frequently in these regions compared to other regions of the genome.

Results

We found that genomic regions showing elevated polymorphism contained a significantly higher density of R-gene analogs predicted to be under pathogen-mediated selection than regions of non-elevated polymorphism, and that many of these also showed evidence for an intermediate site-frequency spectrum based on Tajima’s D. However, we found few genomic regions that showed both elevated polymorphism and reduced F_ST among populations, despite strong background levels of genetic differentiation among populations. This suggests either insufficient power to detect the reduced population structure predicted for genes under balancing selection using sparsely distributed RAD markers, or that other forms of diversifying selection are more common for the R-gene analogs tested.

Conclusions

Genome scans based on a small number of individuals sampled from a wide range of populations were sufficient to confirm the relative scarcity of signatures of balancing selection across the genome, but also identified new potential disease resistance candidates within genomic regions showing signatures of balancing selection that would be strong candidates for further sequencing efforts.

Electronic supplementary material

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

ApoplastP: prediction of effectors and plant proteins in the apoplast using machine learning

The plant apoplast is integral to intercellular signalling, transport and plant-pathogen interactions. Plant pathogens deliver effectors both into the apoplast and inside host cells, but no computational method currently exists to discriminate between these localizations. We present ApoplastP, the first method for predicting whether an effector or plant protein localizes to the apoplast. ApoplastP uncovers features of apoplastic localization common to both effectors and plant proteins, namely depletion in glutamic acid, acidic amino acids and charged amino acids and enrichment in small amino acids. ApoplastP predicts apoplastic localization in effectors with a sensitivity of 75% and a false positive rate of 5%, improving the accuracy of cysteine-rich classifiers by > 13%. ApoplastP does not depend on the presence of a signal peptide and correctly predicts the localization of unconventionally secreted proteins. The secretomes of fungal saprophytes as well as necrotrophic, hemibiotrophic and extracellular fungal pathogens are enriched for predicted apoplastic proteins. Rust pathogens have low proportions of predicted apoplastic proteins, but these are highly enriched for predicted effectors. ApoplastP pioneers apoplastic localization prediction using machine learning. It will facilitate functional studies and will be valuable for predicting if an effector localizes to the apoplast or if it enters plant cells.

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.

Effector Biology in Focus: A Primer for Computational Prediction and Functional Characterization

Plant-pathogen interactions are controlled by a multilayered immune system, which is activated by pathogen recognition in the host. Pathogens secrete effector molecules to interfere with the immune recognition or signaling network and reprogram cell structure or metabolism. Understanding the effector repertoires of diverse pathogens will contribute to unraveling the molecular mechanism of virulence and developing sustainable disease-control strategies for crops and natural ecosystems. Effector functionality has been investigated extensively in only a small number of pathogen species. However, many more pathogen genomes are becoming available, and much can be learned from a broader view of effector biology in diverse pathosystems. The purpose of this review is to summarize methodology for computational prediction of protein effectors, functional characterization of effector proteins and their targets, and the use of effectors as probes to screen for new sources of host resistance. Although these techniques were generally developed in model pathosystems, many of the approaches are directly applicable for exploration and exploitation of effector biology in pathosystems that are less well studied. We hope to facilitate such exploration, which will broaden understanding of the mechanisms that underpin the biological diversity of plant-pathogen interactions, and maximize the impact of new approaches that leverage effector biology for disease control.

Comparative secretome analysis of Rhizoctonia solani isolates with different host ranges reveals unique secretomes and cell death inducing effectors

Rhizoctonia solani is a fungal pathogen causing substantial damage to many of the worlds’ largest food crops including wheat, rice, maize and soybean. Despite impacting global food security, little is known about the pathogenicity mechanisms employed by R. solani. To enable prediction of effectors possessing either broad efficacy or host specificity, a combined secretome was constructed from a monocot specific isolate, a dicot specific isolate and broad host range isolate infecting both monocot and dicot hosts. Secretome analysis suggested R. solani employs largely different virulence mechanisms to well-studied pathogens, despite in many instances infecting the same host plants. Furthermore, the secretome of the broad host range AG8 isolate may be shaped by maintaining functions for saprophytic life stages while minimising opportunities for host plant recognition. Analysis of possible co-evolution with host plants and in-planta up-regulation in particular, aided identification of effectors including xylanase and inhibitor I9 domain containing proteins able to induce cell death in-planta. The inhibitor I9 domain was more abundant in the secretomes of a wide range of necrotising fungi relative to biotrophs. These findings provide novel targets for further dissection of the virulence mechanisms and potential avenues to control this under-characterised but important pathogen.

Divergent and convergent modes of interaction between wheat and Puccinia graminis f. sp. tritici isolates revealed by the comparative gene co-expression network and genome analyses

BACKGROUND

Two opposing evolutionary constraints exert pressure on plant pathogens: one to diversify virulence factors in order to evade plant defenses, and the other to retain virulence factors critical for maintaining a compatible interaction with the plant host. To better understand how the diversified arsenals of fungal genes promote interaction with the same compatible wheat line, we performed a comparative genomic analysis of two North American isolates of Puccinia graminis f. sp. tritici (Pgt).

RESULTS

The patterns of inter-isolate divergence in the secreted candidate effector genes were compared with the levels of conservation and divergence of plant-pathogen gene co-expression networks (GCN) developed for each isolate. Comprative genomic analyses revealed substantial level of interisolate divergence in effector gene complement and sequence divergence. Gene Ontology (GO) analyses of the conserved and unique parts of the isolate-specific GCNs identified a number of conserved host pathways targeted by both isolates. Interestingly, the degree of inter-isolate sub-network conservation varied widely for the different host pathways and was positively associated with the proportion of conserved effector candidates associated with each sub-network. While different Pgt isolates tended to exploit similar wheat pathways for infection, the mode of plant-pathogen interaction varied for different pathways with some pathways being associated with the conserved set of effectors and others being linked with the diverged or isolate-specific effectors.

CONCLUSIONS

Our data suggest that at the intra-species level pathogen populations likely maintain divergent sets of effectors capable of targeting the same plant host pathways. This functional redundancy may play an important role in the dynamic of the "arms-race" between host and pathogen serving as the basis for diverse virulence strategies and creating conditions where mutations in certain effector groups will not have a major effect on the pathogen's ability to infect the host.

Using Population and Comparative Genomics to Understand the Genetic Basis of Effector-Driven Fungal Pathogen Evolution

Epidemics caused by fungal plant pathogens pose a major threat to agro-ecosystems and impact global food security. High-throughput sequencing enabled major advances in understanding how pathogens cause disease on crops. Hundreds of fungal genomes are now available and analyzing these genomes highlighted the key role of effector genes in disease. Effectors are small secreted proteins that enhance infection by manipulating host metabolism. Fungal genomes carry 100s of putative effector genes, but the lack of homology among effector genes, even for closely related species, challenges evolutionary and functional analyses. Furthermore, effector genes are often found in rapidly evolving chromosome compartments which are difficult to assemble. We review how population and comparative genomics toolsets can be combined to address these challenges. We highlight studies that associated genome-scale polymorphisms with pathogen lifestyles and adaptation to different environments. We show how genome-wide association studies can be used to identify effectors and other pathogenicity-related genes underlying rapid adaptation. We also discuss how the compartmentalization of fungal genomes into core and accessory regions shapes the evolution of effector genes. We argue that an understanding of genome evolution provides important insight into the trajectory of host-pathogen co-evolution.

Computational Methods for Predicting Effectors in Rust Pathogens

Lower costs and improved sequencing technologies have led to a large number of high-quality rust pathogen genomes and deeper characterization of gene expression profiles during early and late infection stages. However, the set of secreted proteins expressed during infection is too large for experimental investigations and contains not only effectors but also proteins that play a role in niche colonization or in fighting off competing microbes. Therefore, accurate computational prediction is essential for identifying high-priority rust effector candidates from secretomes.

Addressing the Challenges of Pathogen Evolution on the World's Arable Crops

Advances in genomic and molecular technologies coupled with an increasing understanding of the fine structure of many resistance and infectivity genes, have opened up a new era of hope in controlling the many plant pathogens that continue to be a major source of loss in arable crops. Some new approaches are under consideration including the use of nonhost resistance and the targeting of critical developmental constraints. However, the major thrust of these genomic and molecular approaches is to enhance the identification of resistance genes, to increase their ease of manipulation through marker and gene editing technologies and to lock a range of resistance genes together in simply manipulable resistance gene cassettes. All these approaches essentially continue a strategy that assumes the ability to construct genetic-based resistance barriers that are insurmountable to target pathogens. Here we show how the recent advances in knowledge and marker technologies can be used to generate more durable disease resistance strategies that are based on broad evolutionary principles aimed at presenting pathogens with a shifting, landscape of fluctuating directional selection.

EffectorP: predicting fungal effector proteins from secretomes using machine learning

Eukaryotic filamentous plant pathogens secrete effector proteins that modulate the host cell to facilitate infection. Computational effector candidate identification and subsequent functional characterization delivers valuable insights into plant-pathogen interactions. However, effector prediction in fungi has been challenging due to a lack of unifying sequence features such as conserved N-terminal sequence motifs. Fungal effectors are commonly predicted from secretomes based on criteria such as small size and cysteine-rich, which suffers from poor accuracy. We present EffectorP which pioneers the application of machine learning to fungal effector prediction. EffectorP improves fungal effector prediction from secretomes based on a robust signal of sequence-derived properties, achieving sensitivity and specificity of over 80%. Features that discriminate fungal effectors from secreted noneffectors are predominantly sequence length, molecular weight and protein net charge, as well as cysteine, serine and tryptophan content. We demonstrate that EffectorP is powerful when combined with in planta expression data for predicting high-priority effector candidates. EffectorP is the first prediction program for fungal effectors based on machine learning. Our findings will facilitate functional fungal effector studies and improve our understanding of effectors in plant-pathogen interactions. EffectorP is available at http://effectorp.csiro.au.

Changing the Game: Using Integrative Genomics to Probe Virulence Mechanisms of the Stem Rust Pathogen Puccinia graminis f. sp. tritici

The recent resurgence of wheat stem rust caused by new virulent races of Puccinia graminis f. sp. tritici (Pgt) poses a threat to food security. These concerns have catalyzed an extensive global effort toward controlling this disease. Substantial research and breeding programs target the identification and introduction of new stem rust resistance (Sr) genes in cultivars for genetic protection against the disease. Such resistance genes typically encode immune receptor proteins that recognize specific components of the pathogen, known as avirulence (Avr) proteins. A significant drawback to deploying cultivars with single Sr genes is that they are often overcome by evolution of the pathogen to escape recognition through alterations in Avr genes. Thus, a key element in achieving durable rust control is the deployment of multiple effective Sr genes in combination, either through conventional breeding or transgenic approaches, to minimize the risk of resistance breakdown. In this situation, evolution of pathogen virulence would require changes in multiple Avr genes in order to bypass recognition. However, choosing the optimal Sr gene combinations to deploy is a challenge that requires detailed knowledge of the pathogen Avr genes with which they interact and the virulence phenotypes of Pgt existing in nature. Identifying specific Avr genes from Pgt will provide screening tools to enhance pathogen virulence monitoring, assess heterozygosity and propensity for mutation in pathogen populations, and confirm individual Sr gene functions in crop varieties carrying multiple effective resistance genes. Toward this goal, much progress has been made in assembling a high quality reference genome sequence for Pgt, as well as a Pan-genome encompassing variation between multiple field isolates with diverse virulence spectra. In turn this has allowed prediction of Pgt effector gene candidates based on known features of Avr genes in other plant pathogens, including the related flax rust fungus. Upregulation of gene expression in haustoria and evidence for diversifying selection are two useful parameters to identify candidate Avr genes. Recently, we have also applied machine learning approaches to agnostically predict candidate effectors. Here, we review progress in stem rust pathogenomics and approaches currently underway to identify Avr genes recognized by wheat Sr genes.

Genome-wide gene expression dynamics of the fungal pathogen Dothistroma septosporum throughout its infection cycle of the gymnosperm host Pinus radiata

We present genome-wide gene expression patterns as a time series through the infection cycle of the fungal pine needle blight pathogen, Dothistroma septosporum, as it invades its gymnosperm host, Pinus radiata. We determined the molecular changes at three stages of the disease cycle: epiphytic/biotrophic (early), initial necrosis (mid) and mature sporulating lesion (late). Over 1.7 billion combined plant and fungal reads were sequenced to obtain 3.2 million fungal-specific reads, which comprised as little as 0.1% of the sample reads early in infection. This enriched dataset shows that the initial biotrophic stage is characterized by the up-regulation of genes encoding fungal cell wall-modifying enzymes and signalling proteins. Later necrotrophic stages show the up-regulation of genes for secondary metabolism, putative effectors, oxidoreductases, transporters and starch degradation. This in-depth through-time transcriptomic study provides our first snapshot of the gene expression dynamics that characterize infection by this fungal pathogen in its gymnosperm host.

Common protein sequence signatures associate with Sclerotinia borealis lifestyle and secretion in fungal pathogens of the Sclerotiniaceae

Fungal plant pathogens produce secreted proteins adapted to function outside fungal cells to facilitate colonization of their hosts. In many cases such as for fungi from the Sclerotiniaceae family the repertoire and function of secreted proteins remains elusive. In the Sclerotiniaceae, whereas Sclerotinia sclerotiorum and Botrytis cinerea are cosmopolitan broad host-range plant pathogens, Sclerotinia borealis has a psychrophilic lifestyle with a low optimal growth temperature, a narrow host range and geographic distribution. To spread successfully, S. borealis must synthesize proteins adapted to function in its specific environment. The search for signatures of adaptation to S. borealis lifestyle may therefore help revealing proteins critical for colonization of the environment by Sclerotiniaceae fungi. Here, we analyzed amino acids usage and intrinsic protein disorder in alignments of groups of orthologous proteins from the three Sclerotiniaceae species. We found that enrichment in Thr, depletion in Glu and Lys, and low disorder frequency in hot loops are significantly associated with S. borealis proteins. We designed an index to report bias in these properties and found that high index proteins were enriched among secreted proteins in the three Sclerotiniaceae fungi. High index proteins were also enriched in function associated with plant colonization in S. borealis, and in in planta-induced genes in S. sclerotiorum. We highlight a novel putative antifreeze protein and a novel putative lytic polysaccharide monooxygenase identified through our pipeline as candidate proteins involved in colonization of the environment. Our findings suggest that similar protein signatures associate with S. borealis lifestyle and with secretion in the Sclerotiniaceae. These signatures may be useful for identifying proteins of interest as targets for the management of plant diseases.

Genome-Wide Analysis in Three Fusarium Pathogens Identifies Rapidly Evolving Chromosomes and Genes Associated with Pathogenicity

Pathogens and hosts are in an ongoing arms race and genes involved in host–pathogen interactions are likely to undergo diversifying selection. Fusarium plant pathogens have evolved diverse infection strategies, but how they interact with their hosts in the biotrophic infection stage remains puzzling. To address this, we analyzed the genomes of three Fusarium plant pathogens for genes that are under diversifying selection. We found a two-speed genome structure both on the chromosome and gene group level. Diversifying selection acts strongly on the dispensable chromosomes in Fusarium oxysporum f. sp. lycopersici and on distinct core chromosome regions in Fusarium graminearum, all of which have associations with virulence. Members of two gene groups evolve rapidly, namely those that encode proteins with an N-terminal [SG]-P-C-[KR]-P sequence motif and proteins that are conserved predominantly in pathogens. Specifically, 29 F. graminearum genes are rapidly evolving, in planta induced and encode secreted proteins, strongly pointing toward effector function. In summary, diversifying selection in Fusarium is strongly reflected as genomic footprints and can be used to predict a small gene set likely to be involved in host–pathogen interactions for experimental verification.

Common protein sequence signatures associate with Sclerotinia borealis lifestyle and secretion in fungal pathogens of the Sclerotiniaceae

Fungal plant pathogens produce secreted proteins adapted to function outside fungal cells to facilitate colonization of their hosts. In many cases such as for fungi from the Sclerotiniaceae family the repertoire and function of secreted proteins remains elusive. In the Sclerotiniaceae, whereas Sclerotinia sclerotiorum and Botrytis cinerea are cosmopolitan broad host-range plant pathogens, Sclerotinia borealis has a psychrophilic lifestyle with a low optimal growth temperature, a narrow host range and geographic distribution. To spread successfully, S. borealis must synthesize proteins adapted to function in its specific environment. The search for signatures of adaptation to S. borealis lifestyle may therefore help revealing proteins critical for colonization of the environment by Sclerotiniaceae fungi. Here, we analyzed amino acids usage and intrinsic protein disorder in alignments of groups of orthologous proteins from the three Sclerotiniaceae species. We found that enrichment in Thr, depletion in Glu and Lys, and low disorder frequency in hot loops are significantly associated with S. borealis proteins. We designed an index to report bias in these properties and found that high index proteins were enriched among secreted proteins in the three Sclerotiniaceae fungi. High index proteins were also enriched in function associated with plant colonization in S. borealis, and in in planta-induced genes in S. sclerotiorum. We highlight a novel putative antifreeze protein and a novel putative lytic polysaccharide monooxygenase identified through our pipeline as candidate proteins involved in colonization of the environment. Our findings suggest that similar protein signatures associate with S. borealis lifestyle and with secretion in the Sclerotiniaceae. These signatures may be useful for identifying proteins of interest as targets for the management of plant diseases.