Population genetics of fungal and oomycete effectors involved in gene-for-gene interactions

Quantitative pathogenicity and host adaptation in a fungal plant pathogen revealed by whole-genome sequencing

Knowledge of genetic determinism and evolutionary dynamics mediating host-pathogen interactions is essential to manage fungal plant diseases. Studies on the genetic architecture of fungal pathogenicity often focus on large-effect effector genes triggering strong, qualitative resistance. It is not clear how this translates to predominately quantitative interactions. Here, we use the Zymoseptoria tritici-wheat model to elucidate the genetic architecture of quantitative pathogenicity and mechanisms mediating host adaptation. With a multi-host genome-wide association study, we identify 19 high-confidence candidate genes associated with quantitative pathogenicity. Analysis of genetic diversity reveals that sequence polymorphism is the main evolutionary process mediating differences in quantitative pathogenicity, a process that is likely facilitated by genetic recombination and transposable element dynamics. Finally, we use functional approaches to confirm the role of an effector-like gene and a methyltransferase in phenotypic variation. This study highlights the complex genetic architecture of quantitative pathogenicity, extensive diversifying selection and plausible mechanisms facilitating pathogen adaptation.

Analysis of five near-complete genome assemblies of the tomato pathogen Cladosporium fulvum uncovers additional accessory chromosomes and structural variations induced by transposable elements effecting the loss of avirulence genes

BACKGROUND

Fungal plant pathogens have dynamic genomes that allow them to rapidly adapt to adverse conditions and overcome host resistance. One way by which this dynamic genome plasticity is expressed is through effector gene loss, which enables plant pathogens to overcome recognition by cognate resistance genes in the host. However, the exact nature of these loses remains elusive in many fungi. This includes the tomato pathogen Cladosporium fulvum, which is the first fungal plant pathogen from which avirulence (Avr) genes were ever cloned and in which loss of Avr genes is often reported as a means of overcoming recognition by cognate tomato Cf resistance genes. A recent near-complete reference genome assembly of C. fulvum isolate Race 5 revealed a compartmentalized genome architecture and the presence of an accessory chromosome, thereby creating a basis for studying genome plasticity in fungal plant pathogens and its impact on avirulence genes.

RESULTS

Here, we obtained near-complete genome assemblies of four additional C. fulvum isolates. The genome assemblies had similar sizes (66.96 to 67.78 Mb), number of predicted genes (14,895 to 14,981), and estimated completeness (98.8 to 98.9%). Comparative analysis that included the genome of isolate Race 5 revealed high levels of synteny and colinearity, which extended to the density and distribution of repetitive elements and of repeat-induced point (RIP) mutations across homologous chromosomes. Nonetheless, structural variations, likely mediated by transposable elements and effecting the deletion of the avirulence genes Avr4E, Avr5, and Avr9, were also identified. The isolates further shared a core set of 13 chromosomes, but two accessory chromosomes were identified as well. Accessory chromosomes were significantly smaller in size, and one carried pseudogenized copies of two effector genes. Whole-genome alignments further revealed genomic islands of near-zero nucleotide diversity interspersed with islands of high nucleotide diversity that co-localized with repeat-rich regions. These regions were likely generated by RIP, which generally asymmetrically affected the genome of C. fulvum.

CONCLUSIONS

Our results reveal new evolutionary aspects of the C. fulvum genome and provide new insights on the importance of genomic structural variations in overcoming host resistance in fungal plant pathogens.

Virulence Associations and Global Context of AvrStb6 Genetic Diversity in Iranian Populations of Zymoseptoria tritici

Managing pathogen damage in wheat production is important for sustaining yields. Fungal plant pathogen genomes encode many small secreted proteins acting as effectors that play key roles in the successful colonization of host tissue and triggering host defenses. AvrStb6 is the first described Zymoseptoria tritici avirulence effector, which triggers Stb6-mediated immunity in the wheat host in a gene-for-gene manner. Evasion of major resistance factors such as Stb6 challenges deployment decisions on wheat cultivars. In this study, we analyzed the evolution of the AvrStb6 effector in Iranian isolates of Z. tritici. In total, 78 isolates were isolated and purified from 30 infected wheat specimens collected from the East Azerbaijan and Ardabil provinces of Iran. The pathogenicity of all isolates was evaluated on the susceptible wheat cultivar 'Tajan'. A subset of 40 isolates were also tested for pathogenicity on the resistant cultivar 'Shafir' carrying Stb6. Genetic diversity at the AvrStb6 locus was analyzed for 14 isolates covering the breadth of the observed disease severity. The AvrStb6 sequence variation was high, with virulent isolates carrying highly diverse AvrStb6 haplotypes. In an analysis including more than 1,000 additional AvrStb6 sequences from a global set of isolates, we found that virulent isolates carried AvrStb6 haplotypes either clustering with known virulent haplotypes on different continents or constituting previously unknown haplotypes. Furthermore, we found that AvrStb6 variants from avirulent isolates clustered with known avirulent genotypes from Europe. Our study highlights the relevance of AvrStb6 for Z. tritici virulence and the exceptional global diversity patterns of this effector.

A global-temporal analysis on Phytophthora sojae resistance-gene efficacy

Plant disease resistance genes are widely used in agriculture to reduce disease outbreaks and epidemics and ensure global food security. In soybean, Rps (Resistance to Phytophthora sojae) genes are used to manage Phytophthora sojae, a major oomycete pathogen that causes Phytophthora stem and root rot (PRR) worldwide. This study aims to identify temporal changes in P. sojae pathotype complexity, diversity, and Rps gene efficacy. Pathotype data was collected from 5121 isolates of P. sojae, derived from 29 surveys conducted between 1990 and 2019 across the United States, Argentina, Canada, and China. This systematic review shows a loss of efficacy of specific Rps genes utilized for disease management and a significant increase in the pathotype diversity of isolates over time. This study finds that the most widely deployed Rps genes used to manage PRR globally, Rps1a, Rps1c and Rps1k, are no longer effective for PRR management in the United States, Argentina, and Canada. This systematic review emphasizes the need to widely introduce new sources of resistance to P. sojae, such as Rps3a, Rps6, or Rps11, into commercial cultivars to effectively manage PRR going forward.

Variable Tandem Glycine-Rich Repeats Contribute to Cell Death-Inducing Activity of a Glycosylphosphatidylinositol-Anchored Cell Wall Protein That Is Associated with the Pathogenicity of Sclerotinia sclerotiorum

Glycosylphosphatidylinositol (GPI) anchoring of proteins is a conserved posttranslational modification in eukaryotes. GPI-anchored proteins are widely distributed in fungal plant pathogens, but the specific roles of the GPI-anchored proteins in the pathogenicity of Sclerotinia sclerotiorum, a devastating necrotrophic plant pathogen with a worldwide distribution, remain largely unknown. This research addresses SsGSR1, which encodes an S. sclerotiorum glycine- and serine-rich protein named SsGsr1 with an N-terminal secretory signal and a C-terminal GPI-anchor signal. SsGsr1 is located at the cell wall of hyphae, and deletion of SsGSR1 leads to abnormal cell wall architecture and impaired cell wall integrity of hyphae. The transcription levels of SsGSR1 were maximal in the initial stage of infection, and SsGSR1-deletion strains showed impaired virulence in multiple hosts, indicating that SsGSR1 is critical for the pathogenicity. Interestingly, SsGsr1 targeted the apoplast of host plants to induce cell death that relies on the glycine-rich 11-amino-acid repeats arranged in tandem. The homologs of SsGsr1 in Sclerotinia, Botrytis, and Monilinia species contain fewer repeat units and have lost their cell death activity. Moreover, allelic variants of SsGSR1 exist in field isolates of S. sclerotiorum from rapeseed, and one of the variants lacking one repeat unit results in a protein that exhibits loss of function relative to the cell death-inducing activity and the virulence of S. sclerotiorum. Taken together, our results demonstrate that a variation in tandem repeats provides the functional diversity of GPI-anchored cell wall protein that, in S. sclerotiorum and other necrotrophic pathogens, allows successful colonization of the host plants. IMPORTANCE Sclerotinia sclerotiorum is an economically important necrotrophic plant pathogen and mainly applies cell wall-degrading enzymes and oxalic acid to kill plant cells before colonization. In this research, we characterized a glycosylphosphatidylinositol (GPI)-anchored cell wall protein named SsGsr1, which is critical for the cell wall architecture and the pathogenicity of S. sclerotiorum. Additionally, SsGsr1 induces rapid cell death of host plants that is dependent on glycine-rich tandem repeats. Interestingly, the number of repeat units varies among homologs and alleles of SsGsr1, and such a variation creates alterations in the cell death-inducing activity and the role in pathogenicity. This work advances our understanding of the variation of tandem repeats in accelerating the evolution of a GPI-anchored cell wall protein associated with the pathogenicity of necrotrophic fungal pathogens and prepares the way toward a fuller understanding of the interaction between S. sclerotiorum and host plants.

Quantitative and qualitative plant-pathogen interactions call upon similar pathogenicity genes with a spectrum of effects

Septoria leaf blotch is a foliar wheat disease controlled by a combination of plant genetic resistances and fungicides use. R-gene-based qualitative resistance durability is limited due to gene-for-gene interactions with fungal avirulence (Avr) genes. Quantitative resistance is considered more durable but the mechanisms involved are not well documented. We hypothesize that genes involved in quantitative and qualitative plant-pathogen interactions are similar. A bi-parental population of Zymoseptoria tritici was inoculated on wheat cultivar 'Renan' and a linkage analysis performed to map QTL. Three pathogenicity QTL, Qzt-I05-1, Qzt-I05-6 and Qzt-I07-13, were mapped on chromosomes 1, 6 and 13 in Z. tritici, and a candidate pathogenicity gene on chromosome 6 was selected based on its effector-like characteristics. The candidate gene was cloned by Agrobacterium tumefaciens-mediated transformation, and a pathology test assessed the effect of the mutant strains on 'Renan'. This gene was demonstrated to be involved in quantitative pathogenicity. By cloning a newly annotated quantitative-effect gene in Z. tritici that is effector-like, we demonstrated that genes underlying pathogenicity QTL can be similar to Avr genes. This opens up the previously probed possibility that 'gene-for-gene' underlies not only qualitative but also quantitative plant-pathogen interactions in this pathosystem.

Transgenic Improvement for Biotic Resistance of Crops

Biotic constraints, including pathogenic fungi, viruses and bacteria, herbivory insects, as well as parasitic nematodes, cause significant yield loss and quality deterioration of crops. The effect of conventional management of these biotic constraints is limited. The advances in transgenic technologies provide a direct and directional approach to improve crops for biotic resistance. More than a hundred transgenic events and hundreds of cultivars resistant to herbivory insects, pathogenic viruses, and fungi have been developed by the heterologous expression of exogenous genes and RNAi, authorized for cultivation and market, and resulted in a significant reduction in yield loss and quality deterioration. However, the exploration of transgenic improvement for resistance to bacteria and nematodes by overexpression of endogenous genes and RNAi remains at the testing stage. Recent advances in RNAi and CRISPR/Cas technologies open up possibilities to improve the resistance of crops to pathogenic bacteria and plant parasitic nematodes, as well as other biotic constraints.

Improving sustainable crop protection using population genetics concepts

Growing genetically resistant plants allows pathogen populations to be controlled and reduces the use of pesticides. However, pathogens can quickly overcome such resistance. In this context, how can we achieve sustainable crop protection? This crucial question has remained largely unanswered despite decades of intense debate and research effort. In this study, we used a bibliographic analysis to show that the research field of resistance durability has evolved into three subfields: (i) 'plant breeding' (generating new genetic material), (ii) 'molecular interactions' (exploring the molecular dialogue governing plant-pathogen interactions) and (iii) 'epidemiology and evolution' (explaining and forecasting of pathogen population dynamics resulting from selection pressure(s) exerted by resistant plants). We argue that this triple split of the field impedes integrated research progress and ultimately compromises the sustainable management of genetic resistance. After identifying a gap among the three subfields, we argue that the theoretical framework of population genetics could bridge this gap. Indeed, population genetics formally explains the evolution of all heritable traits, and allows genetic changes to be tracked along with variation in population dynamics. This provides an integrated view of pathogen adaptation, in particular via evolutionary-epidemiological feedbacks. In this Opinion Note, we detail examples illustrating how such a framework can better inform best practices for developing and managing genetically resistant cultivars.

Biological Control of Plant Diseases: An Evolutionary and Eco-Economic Consideration

Biological control is considered as a promising alternative to pesticide and plant resistance to manage plant diseases, but a better understanding of the interaction of its natural and societal functions is necessary for its endorsement. The introduction of biological control agents (BCAs) alters the interaction among plants, pathogens, and environments, leading to biological and physical cascades that influence pathogen fitness, plant health, and ecological function. These interrelationships generate a landscape of tradeoffs among natural and social functions of biological control, and a comprehensive evaluation of its benefits and costs across social and farmer perspectives is required to ensure the sustainable development and deployment of the approach. Consequently, there should be a shift of disease control philosophy from a single concept that only concerns crop productivity to a multifaceted concept concerning crop productivity, ecological function, social acceptability, and economical accessibility. To achieve these goals, attempts should make to develop “green” BCAs used dynamically and synthetically with other disease control approaches in an integrated disease management scheme, and evolutionary biologists should play an increasing role in formulating the strategies. Governments and the public should also play a role in the development and implementation of biological control strategies supporting positive externality.

Multiple Mechanisms Drive the Evolutionary Adaptation of Phytophthora infestans Effector Avr1 to Host Resistance

Effectors, a group of small proteins secreted by pathogens, play a central role in antagonistic interactions between plant hosts and pathogens. The evolution of effector genes threatens plant disease management and sustainable food production, but population genetic analyses to understand evolutionary mechanisms of effector genes are limited compared to molecular and functional studies. Here we investigated the evolution of the Avr1 effector gene from 111 Phytophthora infestans isolates collected from six areas covering three potato cropping regions in China using a population genetic approach. High genetic variation of the effector gene resulted from diverse mechanisms including base substitution, pre-termination, intragenic recombination and diversifying selection. Nearly 80% of the 111 sequences had a point mutation in the 512th nucleotide (T512G), which generated a pre-termination stop codon truncating 38 amino acids in the C-terminal, suggesting that the C-terminal may not be essential to ecological and biological functions of P. infestans. A significant correlation between the frequency of Avr1 sequences with the pre-termination and annual mean temperature in the collection sites suggests that thermal heterogeneity might be one of contributors to the diversifying selection, although biological and biochemical mechanisms of the likely thermal adaptation are not known currently. Our results highlight the risk of rapid adaptation of P. infestans and possibly other pathogens as well to host resistance, and the application of eco-evolutionary principles is necessary for sustainable disease management in agricultural ecosystems.

How Do Pathogens Evolve Novel Virulence Activities?

This article is part of the Top 10 Unanswered Questions in MPMI invited review series.We consider the state of knowledge on pathogen evolution of novel virulence activities, broadly defined as anything that increases pathogen fitness with the consequence of causing disease in either the qualitative or quantitative senses, including adaptation of pathogens to host immunity and physiology, host species, genotypes, or tissues, or the environment. The evolution of novel virulence activities as an adaptive trait is based on the selection exerted by hosts on variants that have been generated de novo or arrived from elsewhere. In addition, the biotic and abiotic environment a pathogen experiences beyond the host may influence pathogen virulence activities. We consider host-pathogen evolution, host range expansion, and external factors that can mediate pathogen evolution. We then discuss the mechanisms by which pathogens generate and recombine the genetic variation that leads to novel virulence activities, including DNA point mutation, transposable element activity, gene duplication and neofunctionalization, and genetic exchange. In summary, if there is an (epi)genetic mechanism that can create variation in the genome, it will be used by pathogens to evolve virulence factors. Our knowledge of virulence evolution has been biased by pathogen evolution in response to major gene resistance, leaving other virulence activities underexplored. Understanding the key driving forces that give rise to novel virulence activities and the integration of evolutionary concepts and methods with mechanistic research on plant-microbe interactions can help inform crop protection.[Formula: see text] Copyright © 2021 The Author(s). This is an open access article distributed under the CC BY-NC-ND 4.0 International license.

Effector Avr4 in Phytophthora infestans Escapes Host Immunity Mainly Through Early Termination

Effector genes play critical roles in the antagonistic interactions between plants and pathogens. However, knowledge of mutation mechanisms and evolutionary processes in effector genes and the contribution of climatic factors to the evolution of effector genes are fragmented but important in sustainable management of plant diseases and securing food supply under changing climates. Here, we used a population genetic approach to explore the evolution of the Avr4 gene in Phytophthora infestans, the causal agent of potato blight. We found that the Avr4 gene exhibited a high genetic diversity generated by point mutation and sequence deletion. Frameshifts caused by a single base-pair deletion at the 194th nucleotide position generate two stop codons, truncating almost the entire C-terminal, which is important for effector function and R4 recognition in all sequences. The effector is under natural selection for adaptation supported by comparative analyses of population differentiation (FST ) and isolation-by-distance between Avr4 sequences and simple sequence repeat marker loci. Furthermore, we found that local air temperature was positively associated with pairwise FST in the Avr4 sequences. These results suggest that the evolution of the effector gene is influenced by local air temperature, and the C-terminal truncation is one of the main mutation mechanisms in the P. infestans effector gene to circumvent the immune response of potato plants. The implication of these results to agricultural and natural sustainability in future climate conditions is discussed.

Multilocus Sequence Analysis of Selected Housekeeping- and Pathogenicity-Related Genes in Venturia inaequalis

The relationship between housekeeping and pathogenicity-related genes and virulence or avirulence towards the primary Malus resistance genes (R) has not been previously studied for Venturia inaequalis fungus, the causal agent of apple scab. In this study, the sequences of two housekeeping genes encoding elongation factor alpha (EF-1α) and β-tubulin and two previously unstudied effector genes of V. inaequalis from mannosidase and glucosidase families of 100 strains collected from apple cultivars with Rvi6, Rvi1, and Rvi17 and without known scab resistance genes were submitted to the analyses. Based on the phylogenetic and diversity data, as well as recombination analyses of the sequenced regions, we assessed the phylogenetic relationships and genetic structure of the pathogen within the species and the evolutionary forces that are currently acting upon this microorganism. The topology of the obtained phylograms demonstrates the lack of a relationship between the phylogenetic position of the strain and the host cultivar and the geographical origin or race of the strain. The isolates from different hosts were differentiated but did not form diagnosable, distinct phylogenetic groups. These results suggest that the analyzed genes may be too conserved to reflect the adaptation of pathogens to apple genotypes with different R genes; thus, they do not adequately reflect race discrimination. In contrast, based on variation and gene flow estimation, genetic divergence was observed among strains virulent to apple trees containing Rvi6. The results of this study confirmed a lack of free recombination between strains and demonstrated that the analyzed regions are in linkage disequilibrium and contain non-random polymorphisms associated with the strain.

Flourishing in water: the early evolution and diversification of plant receptor-like kinases

Receptor-like kinases (RLKs) play significant roles in mediating innate immunity and development of plants. The evolution of plant RLKs has been characterized by extensive variation in copy numbers and domain configurations. However, much remains unknown about the origin, evolution, and early diversification of plant RLKs. Here, we perform phylogenomic analyses of RLKs across plants (Archaeplastida), including embryophytes, charophytes, chlorophytes, prasinodermophytes, glaucophytes, and rhodophytes. We identify the presence of RLKs in all the streptophytes (land plants and charophytes), nine out of 18 chlorophytes, one prasinodermophyte, and one glaucophyte, but not in rhodophytes. Interestingly, the copy number of RLKs increased drastically in streptophytes after the split of the clade of Mesostigmatophyceae and Chlorokybophyceae and other streptophytes. Moreover, phylogenetic analyses suggest RLKs from charophytes form diverse distinct clusters, and are dispersed along the diversity of land plant RLKs, indicating that RLKs have extensively diversified in charophytes and charophyte RLKs seeded the major diversity of land plant RLKs. We identify at least 81 and 76 different kinase-associated domains for charophyte and land plant RLKs, 23 of which are shared, suggesting that RLKs might have evolved in a modular fashion through frequent domain gains or losses. We also detect signatures of positive selection for many charophyte RLK groups, indicating potential functions in host-microbe interaction. Taken together, our findings provide significant insights into the early evolution and diversification of plant RLKs and the ancient evolution of plant-microbe symbiosis.

Deciphering the role of microRNAs during Pi54 gene mediated Magnaporthe oryzae resistance response in rice

The broad-spectrum resistance gene Pi54 confers resistance to multiple isolates of Magnaporthe oryzae in rice. In order to decipher the molecular mechanism underlying the Pi54 mediated resistance in rice line Taipei309 ^Pi54 (carrying Pi54), miRNAome study was performed at 24 h post-inoculation (hpi) with M. oryzae. A total of 222 known miRNAs representing 101 miRNA families were found in this study. Of these, 29 and 24 miRNAs were respectively up- and down-regulated in the resistant Taipei309 ^Pi54 . Defence response (DR) genes, like, NBSGO35, and OsWAK129b, and genes related to transcription factors were up-regulated in Taipei309 ^Pi54 line. The vast array of miRNA candidates identified here are miR159c, miR167c, miR2100, miR2118o, miR2118l, miR319a, miR393, miR395l, miR397a, miR397b, miR398, miR439g, miR531b, miR812f, and miR815c, and they manifest their role in balancing the interplay between various DR genes during Pi54 mediated resistance. We also validated miRNA/target gene pairs involved in hormone signalling, and cross-talk among hormone pathways regulating the rice immunity. This study suggests that the Pi54 gene mediated blast resistance is influenced by several microRNAs through PTI and ETI components in the rice line Taipei309 ^Pi54 , leading to incompatible host-pathogen interaction.

Novel Variation and Evolution of AvrPiz-t of Magnaporthe oryzae in Field Isolates

The product of the avirulence (Avr) gene of Magnaporthe oryzae can be detected by the product of the corresponding resistance (R) gene of rice and activates immunity to rice mediated by the R gene. The high degree of variability of M. oryzae isolates in pathogenicity makes the control of rice blast difficult. That resistance of the R gene in rice has been lost has been ascribed to the instability of the Avr gene in M. oryzae. Further study on the variation of the Avr genes in M. oryze field isolates may yield valuable information on the durable and effective deployment of R genes in rice production areas. AvrPiz-t and Piz-t are a pair of valuable genes in the Rice-Magnaporthe pathosystem. AvrPiz-t is detectable by Piz-t and determines the effectiveness of Piz-t. To effectively deploy the R gene Piz-t, the distribution, variation, and evolution of the corresponding Avr gene AvrPiz-t were found among 312 M. oryzae isolates collected from Yunnan rice production areas of China. PCR amplification and pathogenicity assays of AvrPiz-t showed that 202 isolates (64.7%) held AvrPiz-t alleles and were avirulent to IRBLzt-T (holding Piz-t). There were 42.3–83.3% avirulent isolates containing AvrPiz-t among seven regions in Yunnan Province. Meanwhile, 11 haplotypes of AvrPiz-t encoding three novel AvrPiz-t variants were identified among 100 isolates. A 198 bps insertion homologous to solo-LTR of the retrotransposon inago2 in the promoter region of AvrPiz-t in one isolate and a frameshift mutation of CDS in another isolate were identified among 100 isolates, and those two isolates had evolved to virulent from avirulent. Synonymous mutation and non-AUG-initiated N-terminal extensions keeps the AvrPiz-t gene avirulence function in M. oryzae field isolates in Yunnan. A haplotype network showed that H3 was an ancestral haplotype. Structure variance for absence (28.2%) or partial fragment loss (71.8%) of AvrPiz-t was found among 39 virulent isolates and may cause the AvrPiz-t avirulence function to be lost. Overall, AvrPiz-t evolved to virulent from avirulent forms via point mutation, retrotransposon, shift mutation, and structure variance under field conditions.

Recent insights into barley and Rhynchosporium commune interactions

Rhynchosporium commune is the causal pathogen of scald in barley (Hordeum vulgare), a foliar disease that can reduce yield by up to 40% in susceptible cultivars. R. commune is found worldwide in all temperate growing regions and is regarded as one of the most economically important barley pathogens. It is a polycyclic pathogen with the ability to rapidly evolve new virulent strains in response to resistance genes deployed in commercial cultivars. Hence, introgression and pyramiding of different loci for resistance (qualitative or quantitative) through marker-assisted selection is an effective way to improve scald resistance in barley. This review summarizes all 148 resistance quantitative trait loci reported at the date of submission of this review and projects them onto the barley physical map, where it is clear many loci co-locate on chromosomes 3H and 7H. We have summarized the major named resistance loci and reiterated the renaming of Rrs15 (CI8288) to Rrs17. This review provides a comprehensive resource for future discovery and breeding efforts of qualitative and quantitative scald resistance loci.

The rise and fall of genes: origins and functions of plant pathogen pangenomes

Plant pathogens can rapidly overcome resistance of their hosts by mutating key pathogenicity genes encoding for effectors. Pathogen adaptation is fuelled by extensive genetic variability in populations and different strains may not share the same set of genes. Recently, such an intra-specific variation in gene content became formalized as pangenomes distinguishing core genes (i.e. shared) and accessory genes (i.e. lineage or strain-specific). Across pathogens species, key effectors tend to be part of the rapidly evolving accessory genome. Here, we show how the construction and analysis of pathogen pangenomes provide deep insights into the dynamic host adaptation process. We also discuss how pangenomes should ideally be built and how geography, niche and lifestyle likely determine pangenome sizes.

Advances and Challenges in Bacterial Spot Resistance Breeding in Tomato (Solanum lycopersicum L.)

Bacterial spot is a serious disease of tomato caused by at least four species of Xanthomonas. These include X. euvesicatoria (race T1), X. vesicatoria (race T2), X. perforans (races T3 and T4), and X. gardneri, with the distinct geographical distribution of each group. Currently, X. gardneri and X. perforans are two major bacterial pathogens of tomato in North America, with X. perforans (race T4) dominating in east-coast while X. gardneri dominating in the Midwest. The disease causes up to 66% yield loss. Management of this disease is challenging due to the lack of useful chemical control measures and commercial resistant cultivars. Although major genes for resistance (R) and quantitative resistance have been identified, breeding tomato for resistance to bacterial spot has been impeded by multiple factors including the emergence of new races of the pathogen that overcome the resistance, multigenic control of the resistance, linkage drag, non-additive components of the resistance and a low correlation between seedling assays and field resistance. Transgenic tomato with Bs2 and EFR genes was effective against multiple races of Xanthomonas. However, it has not been commercialized because of public concerns and complex regulatory processes. The genomics-assisted breeding, effectors-based genomics breeding, and genome editing technology could be novel approaches to achieve durable resistance to bacterial spot in tomato. The main goal of this paper is to understand the current status of bacterial spot of tomato including its distribution and pathogen diversity, challenges in disease management, disease resistance sources, resistance genetics and breeding, and future prospectives with novel breeding approaches.

Local adaptation drives the diversification of effectors in the fungal wheat pathogen Parastagonospora nodorum in the United States

Filamentous fungi rapidly evolve in response to environmental selection pressures in part due to their genomic plasticity. Parastagonospora nodorum, a fungal pathogen of wheat and causal agent of septoria nodorum blotch, responds to selection pressure exerted by its host, influencing the gain, loss, or functional diversification of virulence determinants, known as effector genes. Whole genome resequencing of 197 P. nodorum isolates collected from spring, durum, and winter wheat production regions of the United States enabled the examination of effector diversity and genomic regions under selection specific to geographically discrete populations. 1,026,859 SNPs/InDels were used to identify novel loci, as well as SnToxA and SnTox3 as factors in disease. Genes displaying presence/absence variation, predicted effector genes, and genes localized on an accessory chromosome had significantly higher pN/pS ratios, indicating a higher rate of sequence evolution. Population structure analyses indicated two P. nodorum populations corresponding to the Upper Midwest (Population 1) and Southern/Eastern United States (Population 2). Prevalence of SnToxA varied greatly between the two populations which correlated with presence of the host sensitivity gene Tsn1 in the most prevalent cultivars in the corresponding regions. Additionally, 12 and 5 candidate effector genes were observed to be under diversifying selection among isolates from Population 1 and 2, respectively, but under purifying selection or neutrally evolving in the opposite population. Selective sweep analysis revealed 10 and 19 regions that had recently undergone positive selection in Population 1 and 2, respectively, involving 92 genes in total. When comparing genes with and without presence/absence variation, those genes exhibiting this variation were significantly closer to transposable elements. Taken together, these results indicate that P. nodorum is rapidly adapting to distinct selection pressures unique to spring and winter wheat production regions by rapid adaptive evolution and various routes of genomic diversification, potentially facilitated through transposable element activity.

A genome-wide scan for genes under balancing selection in the plant pathogen Ralstonia solanacearum

Background

Plant pathogens are under significant selective pressure by the plant host. Consequently, they are expected to have adapted to this condition or contribute to evading plant defenses. In order to acquire long-term fitness, plant bacterial pathogens are usually forced to maintain advantageous genetic diversity in populations. This strategy ensures that different alleles in the pathogen’s gene pool are maintained in a population at frequencies larger than expected under neutral evolution. This selective process, known as balancing selection, is the subject of this work in the context of a common bacterial phytopathogen. We performed a genome-wide scan of Ralstonia solanacearum species complex, an aggressive plant bacterial pathogen that shows broad host range and causes a devastating disease called ‘bacterial wilt’.

Results

Using a sliding window approach, we analyzed 57 genomes from three phylotypes of the R. solanacearum species complex to detect signatures of balancing selection. A total of 161 windows showed extreme values in three summary statistics of population genetics: Tajima’s D, θ_w and Fu & Li’s D*. We discarded any confounding effects due to demographic events by means of coalescent simulations of genetic data. The prospective windows correspond to 78 genes with known function that map in any of the two main replicons (1.7% of total number of genes). The candidate genes under balancing selection are related to primary metabolism and other basal activities (51.3%) or directly associated to virulence (48.7%), the latter being involved in key functions targeted to dismantle plant defenses or to participate in critical stages in the pathogenic process.

Conclusions

We identified various genes under balancing selection that play a significant role in basic metabolism as well as in virulence of the R. solanacearum species complex. These genes are useful to understand and monitor the evolution of bacterial pathogen populations and emerge as potential candidates for future treatments to induce specific plant immune responses.

Electronic supplementary material

The online version of this article (10.1186/s12862-019-1456-6) contains supplementary material, which is available to authorized users.

Phylogenetically Diverse Fusarium Species Associated with Sorghum (Sorghum Bicolor L. Moench) and Finger Millet (Eleusine Coracana L. Garten) Grains from Ethiopia

Fusarium is one of the most diverse fungal genera affecting several crops around the world. This study describes the phylogeny of Fusarium species associated with grains of sorghum and finger millet from different parts of Ethiopia. Forty-two sorghum and 34 finger millet grain samples were mycologically analysed. All of the sorghum and more than 40% of the finger millet grain samples were contaminated by the Fusarium species. The Fusarium load was higher in sorghum grains than that in finger millet grains. In addition, 67 test isolates were phylogenetically analysed using EF-1α and β-tubulin gene primers. Results revealed the presence of eight phylogenetic placements within the genus Fusarium, where 22 of the isolates showed a close phylogenetic relation to the F. incarnatum–equiseti species complex. Nevertheless, they possess a distinct shape of apical cells of macroconidia, justifying the presence of new species within the Fusarium genus. The new species was the most dominant, represented by 33% of the test isolates. The current work can be seen as an important addition to the knowledge of the biodiversity of fungal species that exists within the Fusarium genus. It also reports a previously unknown Fusarium species that needs to be investigated further for toxin production potential.

Novel haplotypes and networks of AVR-Pik alleles in Magnaporthe oryzae

BACKGROUND

Rice blast disease is one of the most destructive fungal disease of rice worldwide. The avirulence (AVR) genes of Magnaporthe oryzae are recognized by the cognate resistance (R) genes of rice and trigger race-specific resistance. The variation in AVR is one of the major drivers of new races. Detecting the variation in the AVR gene in isolates from a population of Magnaporthe oryzae collected from rice production fields will aid in evaluating the effectiveness of R genes in rice production areas. The Pik gene contains 5 R alleles (Pik, Pikh, Pikp, Pikm and Piks) corresponding to the AVR alleles (AVR-Pik/kh/kp/km/ks) of M. oryzae. The Pik gene specifically recognizes and prevents infections by isolates of M. oryzae that contain AVR-Pik. The molecular variation in AVR-Pik alleles of M. oryzae and Pik alleles of rice remains unclear.

RESULTS

We studied the possible evolutionary pathways of AVR-Pik alleles by analyzing their DNA sequence variation and assaying their avirulence to the cognate Pik alleles of resistance genes under field conditions in China. The results of PCR products from genomic DNA showed that 278 of the 366 isolates of M. oryzae collected from Yunnan Province, China, carried AVR-Pik alleles. Among the isolates from six regions of Yunnan, 66.7-90.3% carried AVR-Pik alleles. Moreover, 10 AVR-Pik haplotypes encoding five novel AVR-Pik variants were identified among 201 isolates. The AVR-Pik alleles evolved to virulent from avirulent forms via stepwise base substitution. These findings demonstrate that AVR-Pik alleles are under positive selection and that mutations are responsible for defeating race-specific resistant Pik alleles in nature.

CONCLUSIONS

We demonstrated the polymorphism and distribution of AVR-Pik alleles in Yunnan Province, China. By pathogenicity assays used to detect the function of the different haplotypes of AVR-Pik, for the first time, we showed the avoidance and stepwise evolution of AVR-Pik alleles in rice production areas of Yunnan. The functional AVR-Pik possesses diversified sequence structures and is under positive selection in nature.

Geographic Distribution of Avirulence Genes of the Rice Blast Fungus Magnaporthe oryzae in the Philippines

A total of 131 contemporary and 33 reference isolates representing a number of multi-locus genotypes of Magnaporthe oryzae were subjected to a PCR test to detect the presence/absence of avirulence (Avr) genes. Results revealed that the more frequently occurring genes were Avr-Pik (81.50%), Avr-Pita (64.16%) and Avr-Pii (47.98%), whereas the less frequently occurring genes were Avr-Pizt (19.08%) and Avr-Pia (5.20%). It was also laid out that the presence of Avr genes in M. oryzae is strongly associated with agroecosystems where the complementary resistant (R) genes exist. No significant association, however, was noted on the functional Avr genes and the major geographic locations. Furthermore, it was identified that the upland varieties locally known as "Milagrosa" and "Waray" contained all the R genes complementary to the Avr genes tested.

Determining Whether Geographic Origin and Potato Genotypes Shape the Population Structure of Phytophthora infestans in the Central Region of Colombia

Pathogen variation plays an important role in the dynamics of infectious diseases. In this study, the genetic variation of 279 Phytophthora infestans isolates was assessed using a combination of 12 microsatellite simple-sequence repeat markers. Isolates were collected from 11 different potato cultivars in 11 different geographic localities of the central region of Colombia. The objective of this study was to determine whether populations were differentiated by host genotype or geographic origin. Within a single clonal lineage, EC-1, 76 genotypes were detected. An analysis of molecular variance attributed most of the variation to differences within host genotypes rather than among the host genotypes, suggesting that host cultivars do not structure the populations of the pathogen. Furthermore, the lack of a genetic population structure according to the host cultivar was confirmed by all of the analyses, including the Bayesian clustering analysis and the minimum spanning network that used the Bruvo genetic distance, which suggested that there are no significant barriers to gene flow for P. infestans among potato cultivars. According to the geographic origin, the populations of P. infestans were also not structured, and most of the variation among the isolates was attributed to differences within localities. Only some but not all localities in the north and west of the central region of Colombia showed some genetic differentiation from the other regions. The absence of sexual reproduction of this pathogen in Colombia was also demonstrated. Important insights are discussed regarding the genetic population dynamics of the P. infestans populations of the central region of Colombia that were provided by the results. In Colombia, there is a high genetic variation within the EC-1 clonal lineage with closely related genotypes, none dominant, that coexist in a wide geographic area and on several potato cultivars.

Genetic Analysis of Victorin Sensitivity and Identification of a Causal Nucleotide-Binding Site Leucine-Rich Repeat Gene in Phaseolus vulgaris

Cochliobolus victoria, the causal agent of Victoria blight, is pathogenic due to its production of a toxin called victorin. Victorin sensitivity in oats, barley, Brachypodium spp., and Arabidopsis has been associated with nucleotide-binding site leucine-rich repeat (NLR) genes, a class of genes known for conferring disease resistance. In this work, we investigated the sensitivity of Phaseolus vulgaris to victorin. We found that victorin sensivity in Phaseolus vulgaris is a developmentally regulated, quantitative trait. A single quantitative trait locus (QTL) accounted for 34% of the phenotypic variability in victorin sensitivity among Stampede × Red Hawk (S×R) recombinant inbred lines. We cloned two NLR-encoding genes within this QTL and showed one, Phvul05G031200 (PvLOV), confers victorin-dependent cell death when overexpressed in Nicotiana benthamiana. Protein sequences of PvLOV from victorin-sensitive and the victorin-resistant bean parents differ by two amino acids in the leucine-rich repeat region, but both proteins confer victorin-dependent cell death when overexpressed in N. benthamiana.

Comparative whole-genome analysis reveals artificial selection effects on Ustilago esculenta genome

Ustilago esculenta, infects Zizania latifolia, and induced host stem swollen to be a popular vegetable called Jiaobai in China. It is the long-standing artificial selection that maximizes the occurrence of favourable Jiaobai, and thus maintaining the plant-fungi interaction and modulating the fungus evolving from plant pathogen to entophyte. In this study, whole genome of U. esculenta was sequenced and transcriptomes of the fungi and its host were analysed. The 20.2 Mb U. esculenta draft genome of 6,654 predicted genes including mating, primary metabolism, secreted proteins, shared a high similarity to related Smut fungi. But U. esculenta prefers RNA silencing not repeat-induced point in defence and has more introns per gene, indicating relatively slow evolution rate. The fungus also lacks some genes in amino acid biosynthesis pathway which were filled by up-regulated host genes and developed distinct amino acid response mechanism to balance the infection-resistance interaction. Besides, U. esculenta lost some surface sensors, important virulence factors and host range-related effectors to maintain the economic endophytic life. The elucidation of the U. esculenta genomic information as well as expression profiles can not only contribute to more comprehensive insights into the molecular mechanism underlying artificial selection but also into smut fungi-host interactions.

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.

On plant defense signaling networks and early land plant evolution

All land plants must cope with phytopathogens. Algae face pathogens, too, and it is reasonable to assume that some of the strategies for dealing with pathogens evolved prior to the origin of embryophytes – plant terrestrialization simply changed the nature of the plant-pathogen interactions. Here we highlight that many potential components of the angiosperm defense toolkit are i) found in streptophyte algae and non-flowering embryophytes and ii) might be used in non-flowering plant defense as inferred from published experimental data. Nonetheless, the common signaling networks governing these defense responses appear to have become more intricate during embryophyte evolution. This includes the evolution of the antagonistic signaling pathways of jasmonic and salicylic acid, multiple independent expansions of resistance genes, and the evolution of resistance gene-regulating microRNAs. Future comparative studies will illuminate which modules of the streptophyte defense signaling network constitute the core and which constitute lineage- and/or environment-specific (peripheral) signaling circuits.

Genome-wide evidence for divergent selection between populations of a major agricultural pathogen

The genetic and environmental homogeneity in agricultural ecosystems is thought to impose strong and uniform selection pressures. However, the impact of this selection on plant pathogen genomes remains largely unknown. We aimed to identify the proportion of the genome and the specific gene functions under positive selection in populations of the fungal wheat pathogen Zymoseptoria tritici. First, we performed genome scans in four field populations that were sampled from different continents and on distinct wheat cultivars to test which genomic regions are under recent selection. Based on extended haplotype homozygosity and composite likelihood ratio tests, we identified 384 and 81 selective sweeps affecting 4% and 0.5% of the 35 Mb core genome, respectively. We found differences both in the number and the position of selective sweeps across the genome between populations. Using a XtX‐based outlier detection approach, we identified 51 extremely divergent genomic regions between the allopatric populations, suggesting that divergent selection led to locally adapted pathogen populations. We performed an outlier detection analysis between two sympatric populations infecting two different wheat cultivars to identify evidence for host‐driven selection. Selective sweep regions harboured genes that are likely to play a role in successfully establishing host infections. We also identified secondary metabolite gene clusters and an enrichment in genes encoding transporter and protein localization functions. The latter gene functions mediate responses to environmental stress, including interactions with the host. The distinct gene functions under selection indicate that both local host genotypes and abiotic factors contributed to local adaptation.

The ash dieback invasion of Europe was founded by two genetically divergent individuals

Accelerating international trade and climate change make pathogen spread an increasing concern. Hymenoscyphus fraxineus, the causal agent of ash dieback, is a fungal pathogen that has been moving across continents and hosts from Asian to European ash. Most European common ash trees (Fraxinus excelsior) are highly susceptible to H. fraxineus, although a minority (~5%) have partial resistance to dieback. Here, we assemble and annotate a H. fraxineus draft genome which approaches chromosome scale. Pathogen genetic diversity across Europe and in Japan, reveals a strong bottleneck in Europe, though a signal of adaptive diversity remains in key host interaction genes. We find that the European population was founded by two divergent haploid individuals. Divergence between these haplotypes represents the ancestral polymorphism within a large source population. Subsequent introduction from this source would greatly increase adaptive potential of the pathogen. Thus, further introgression of H. fraxineus into Europe represents a potential threat and Europe-wide biological security measures are needed to manage this disease.

Boxwood blight: an ongoing threat to ornamental and native boxwood

Boxwood blight is an emerging disease of ornamental and native boxwood plants in the family Buxaceae. First documented in the 1990s at a single location in England, the disease is now reported throughout Europe, Asia, New Zealand, and North America. To address the growing concern over boxwood blight, ongoing research focuses on multiple biological and genetic aspects of the causal pathogens and susceptible host plants. Characterization of genetic variation among the Calonectria fungi that cause boxwood blight shows that two unique sister species with different geographic distributions incite the disease. Studies of the pathogen life cycle show the formation of long-lived survival structures and that host infection is dependent on inoculum density, temperature, and humidity. Host range investigations detail high levels of susceptibility among boxwood as well as the potential for asymptomatic boxwood infection and for other plants in the family Buxaceae to serve as additional hosts. Multiple DNA-based diagnostic assays are available, ranging from probe-based quantitative PCR assays to the use of comparative genomics to develop robust diagnostic markers or provide whole genome-scale identifications. Though many questions remain, the research that continues to address boxwood blight demonstrates the importance of applying a multidisciplinary approach to understand and control emerging plant diseases.

Distinct Trajectories of Massive Recent Gene Gains and Losses in Populations of a Microbial Eukaryotic Pathogen

Differences in gene content are a significant source of variability within species and have an impact on phenotypic traits. However, little is known about the mechanisms responsible for the most recent gene gains and losses. We screened the genomes of 123 worldwide isolates of the major pathogen of wheat Zymoseptoria tritici for robust evidence of gene copy number variation. Based on orthology relationships in three closely related fungi, we identified 599 gene gains and 1,024 gene losses that have not yet reached fixation within the focal species. Our analyses of gene gains and losses segregating in populations showed that gene copy number variation arose preferentially in subtelomeres and in proximity to transposable elements. Recently lost genes were enriched in virulence factors and secondary metabolite gene clusters. In contrast, recently gained genes encoded mostly secreted protein lacking a conserved domain. We analyzed the frequency spectrum at loci segregating a gene presence–absence polymorphism in four worldwide populations. Recent gene losses showed a significant excess in low-frequency variants compared with genome-wide single nucleotide polymorphism, which is indicative of strong negative selection against gene losses. Recent gene gains were either under weak negative selection or neutral. We found evidence for strong divergent selection among populations at individual loci segregating a gene presence–absence polymorphism. Hence, gene gains and losses likely contributed to local adaptation. Our study shows that microbial eukaryotes harbor extensive copy number variation within populations and that functional differences among recently gained and lost genes led to distinct evolutionary trajectories.

Fine-Scale Recombination Maps of Fungal Plant Pathogens Reveal Dynamic Recombination Landscapes and Intragenic Hotspots

Meiotic recombination is an important driver of evolution. Variability in the intensity of recombination across chromosomes can affect sequence composition, nucleotide variation, and rates of adaptation. In many organisms, recombination events are concentrated within short segments termed recombination hotspots. The variation in recombination rate and positions of recombination hotspot can be studied using population genomics data and statistical methods. In this study, we conducted population genomics analyses to address the evolution of recombination in two closely related fungal plant pathogens: the prominent wheat pathogen Zymoseptoria tritici and a sister species infecting wild grasses Z. ardabiliae We specifically addressed whether recombination landscapes, including hotspot positions, are conserved in the two recently diverged species and if recombination contributes to rapid evolution of pathogenicity traits. We conducted a detailed simulation analysis to assess the performance of methods of recombination rate estimation based on patterns of linkage disequilibrium, in particular in the context of high nucleotide diversity. Our analyses reveal overall high recombination rates, a lack of suppressed recombination in centromeres, and significantly lower recombination rates on chromosomes that are known to be accessory. The comparison of the recombination landscapes of the two species reveals a strong correlation of recombination rate at the megabase scale, but little correlation at smaller scales. The recombination landscapes in both pathogen species are dominated by frequent recombination hotspots across the genome including coding regions, suggesting a strong impact of recombination on gene evolution. A significant but small fraction of these hotspots colocalize between the two species, suggesting that hotspot dynamics contribute to the overall pattern of fast evolving recombination in these species.

Widespread signatures of selection for secreted peptidases in a fungal plant pathogen

BACKGROUND

Fungal plant pathogens secrete a large arsenal of hydrolytic enzymes during the course of infection, including peptidases. Secreted peptidases have been extensively studied for their role as effectors. In this study, we combined transcriptomics, comparative genomics and evolutionary analyses to investigate all 39 secreted peptidases in the fungal wheat pathogen Zymoseptoria tritici and its close relatives Z. pseudotritici and Z. ardabiliae.

RESULTS

RNA-seq data revealed that a majority of the secreted peptidases displayed differential transcription during the course of Z. tritici infection, indicative of specialization for different stages in the life cycle. Evolutionary analyses detected widespread evidence of adaptive evolution acting on at least 28 of the peptidases. A few peptidases displayed lineage-specific rates of molecular evolution, suggesting altered selection pressure in Z. tritici following host specialization on domesticated wheat. The peptidases belonging to MEROPS families A1 and G1 emerged as a particularly interesting group that may play key roles in host-pathogen co-evolution, host adaptation and pathogenicity. Sister genes in the A1 and G1 families showed accelerated substitution rates after gene duplications.

CONCLUSIONS

These results suggest widespread evolution of secreted peptidases leading to novel gene functions, consistent with predicted models of "escape from adaptive conflict" and "neo-functionalization". Our analyses identified candidate genes worthy of functional analyses that may encode effector functions, for example by suppressing plant defenses during the biotrophic phase of infection.

Evolutionary analyses of the avirulence effector AvrStb6 in global populations of Zymoseptoria tritici identify candidate amino acids involved in recognition

We analysed the population genetic diversity of AvrStb6, the first avirulence gene cloned from the wheat pathogen Zymoseptoria tritici, using 142 Z. tritici strains sampled from four wheat fields growing on three continents. Although AvrStb6 was located in a recombination hotspot, it was found in every strain, with 71 polymorphic sites that produced 41 distinct DNA haplotypes encoding 30 AvrStb6 protein isoforms. An AvrStb6 homologue was found in the closest known relative, Z. pseudotritici, but not in three other closely related Zymoseptoria species, indicating that this gene has emerged in Zymoseptoria quite recently. Two AvrStb6 homologues with nucleotide similarities greater than 70% were identified on chromosome 10 in all Z. tritici isolates, suggesting that AvrStb6 belongs to a multigene family of candidate effectors that has expanded recently through gene duplication. The AvrStb6 sequences exhibited strong evidence for non-neutral evolution, including a large number of non-synonymous mutations, with significant positive diversifying selection operating on nine of the 82 codons. It appears that balancing selection is operating across the entire gene in natural field populations. There was also evidence for co-evolving codons within the gene that may reflect compensatory mutations associated with the evasion of recognition by Stb6. Intragenic recombination also appears to have affected the diversity of AvrStb6.

Quantitative trait locus mapping reveals complex genetic architecture of quantitative virulence in the wheat pathogen Zymoseptoria tritici

We conducted a comprehensive analysis of virulence in the fungal wheat pathogen Zymoseptoria tritici using quantitative trait locus (QTL) mapping. High-throughput phenotyping based on automated image analysis allowed the measurement of pathogen virulence on a scale and with a precision that was not previously possible. Across two mapping populations encompassing more than 520 progeny, 540 710 pycnidia were counted and their sizes and grey values were measured. A significant correlation was found between pycnidia size and both spore size and number. Precise measurements of percentage leaf area covered by lesions provided a quantitative measure of host damage. Combining these large and accurate phenotypic datasets with a dense panel of restriction site-associated DNA sequencing (RADseq) genetic markers enabled us to genetically dissect pathogen virulence into components related to host damage and those related to pathogen reproduction. We showed that different components of virulence can be under separate genetic control. Large- and small-effect QTLs were identified for all traits, with some QTLs specific to mapping populations, cultivars and traits and other QTLs shared among traits within the same mapping population. We associated the presence of four accessory chromosomes with small, but significant, increases in several virulence traits, providing the first evidence for a meaningful function associated with accessory chromosomes in this organism. A large-effect QTL involved in host specialization was identified on chromosome 7, leading to the identification of candidate genes having a large effect on virulence.

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.

A fungal wheat pathogen evolved host specialization by extensive chromosomal rearrangements

Fungal pathogens can rapidly evolve virulence towards resistant crops in agricultural ecosystems. Gains in virulence are often mediated by the mutation or deletion of a gene encoding a protein recognized by the plant immune system. However, the loci and the mechanisms of genome evolution enabling rapid virulence evolution are poorly understood. We performed genome-wide association mapping on a global collection of 106 strains of Zymoseptoria tritici, the most damaging pathogen of wheat in Europe, to identify polymorphisms linked to virulence on two wheat varieties. We found 25 distinct genomic loci associated with reproductive success of the pathogen. However, no locus was shared between the host genotypes, suggesting host specialization. The main locus associated with virulence encoded a highly expressed, small secreted protein. Population genomic analyses showed that the gain in virulence was explained by a segregating gene deletion polymorphism. The deletion was likely adaptive by preventing detection of the encoded protein. Comparative genomics of closely related species showed that the locus emerged de novo since speciation. A large cluster of transposable elements in direct proximity to the locus generated extensive rearrangements leading to multiple independent gene losses. Our study demonstrates that rapid turnover in the chromosomal structure of a pathogen can drive host specialization.

Unusual evolutionary mechanisms to escape effector-triggered immunity in the fungal phytopathogen Leptosphaeria maculans

Leptosphaeria maculans is the fungus responsible for the stem canker disease of oilseed rape (Brassica napus). AvrLm3 and AvrLm4-7, two avirulence effector genes of L. maculans, are involved in an unusual relationship: AvrLm4-7 suppresses the Rlm3-mediated resistance. Here, we assessed AvrLm3 polymorphism in a collection of 235 L. maculans isolates. No field isolates exhibited deletion or inactivating mutations in AvrLm3, as observed for other L. maculans avirulence genes. Eleven isoforms of the AvrLm3 protein were found. In isolates virulent towards both Rlm3 and Rlm7 (a3a7), the loss of the Rlm3-mediated resistance response was due to two distinct mechanisms. First, when AvrLm4-7 was inactivated (deletion or inactivating mutations), amino acid substitutions in AvrLm3 generated virulent isoforms of the protein. Second, when only point mutations were observed in AvrLm4-7, a3a7 isolates still contained an avirulent allele of AvrLm3. Directed mutagenesis confirmed that some point mutations in AvrLm4-7 were sufficient for the fungus to escape Rlm7-mediated resistance while maintaining the suppression of the AvrLm3 phenotype. Signatures of positive selection were also identified in AvrLm3. The complex evolutionary mechanisms enabling L. maculans to escape Rlm3-mediated resistance while preserving AvrLm3 integrity, along with observed reduced aggressiveness of isolates silenced for AvrLm3, serves to emphasize the importance of this effector in pathogenicity towards B. napus. While the common response to resistance gene pressure is local selection of isolates depleted in the cognate avirulence gene, this example contributes to complexify the gene-for-gene concept of plant-pathogen evolution with a 'camouflaged' model allowing retention of nondispensable avirulence effectors.

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.

Population Genomics of Fungal and Oomycete Pathogens

We are entering a new era in plant pathology in which whole-genome sequences of many individuals of a pathogen species are becoming readily available. Population genomics aims to discover genetic mechanisms underlying phenotypes associated with adaptive traits such as pathogenicity, virulence, fungicide resistance, and host specialization, as genome sequences or large numbers of single nucleotide polymorphisms become readily available from multiple individuals of the same species. This emerging field encompasses detailed genetic analyses of natural populations, comparative genomic analyses of closely related species, identification of genes under selection, and linkage analyses involving association studies in natural populations or segregating populations resulting from crosses. The era of pathogen population genomics will provide new opportunities and challenges, requiring new computational and analytical tools. This review focuses on conceptual and methodological issues as well as the approaches to answering questions in population genomics. The major steps start with defining relevant biological and evolutionary questions, followed by sampling, genotyping, and phenotyping, and ending in analytical methods and interpretations. We provide examples of recent applications of population genomics to fungal and oomycete plant pathogens.

Weeds, as ancillary hosts, pose disproportionate risk for virulent pathogen transfer to crops

Background

The outcome of the arms race between hosts and pathogens depends heavily on the interactions between their genetic diversity, population size and transmission ability. Theory predicts that genetically diverse hosts will select for higher virulence and more diverse pathogens than hosts with low genetic diversity. Cultivated hosts typically have lower genetic diversity and thus small effective population sizes, but can potentially harbour large pathogen population sizes. On the other hand, hosts, such as weeds, which are genetically more diverse and thus have larger effective population sizes, usually harbour smaller pathogen population sizes. Large pathogen population sizes may lead to more opportunities for mutation and hence more diverse pathogens. Here we test the predictions that pathogen neutral genetic diversity will increase with large pathogen population sizes and host diversity, whereas diversity under selection will increase with host diversity. We assessed and compared the diversity of a fungal pathogen, Rhynchosporium commune, on weedy barley grass (which have a large effective population size) and cultivated barley (low genetic diversity) using microsatellites, effector locus nip1 diversity and pathogen aggressiveness in order to assess the importance of weeds in the evolution of the neutral and selected diversity of pathogens.

Results

The findings indicated that the large barley acreage and low host diversity maintains higher pathogen neutral genetic diversity and lower linkage disequilibrium, while the weed maintains more pathotypes and higher virulence diversity at nip1. Strong evidence for more pathogen migration from barley grass to barley suggests transmission of virulence from barley grass to barley is common.

Conclusions

Pathogen census population size is a better predictor for neutral genetic diversity than host diversity. Despite maintaining a smaller pathogen census population size, barley grass acts as an important ancillary host to R. commune, harbouring highly virulent pathogen types capable of transmission to barley. Management of disease on crops must therefore include management of weedy ancillary hosts, which may harbour disproportionate supplies of virulent pathogen strains.

Electronic supplementary material

The online version of this article (doi:10.1186/s12862-016-0680-6) contains supplementary material, which is available to authorized users.

A Tale of Genome Compartmentalization: The Evolution of Virulence Clusters in Smut Fungi

Smut fungi are plant pathogens mostly parasitizing wild species of grasses as well as domesticated cereal crops. Genome analysis of several smut fungi including Ustilago maydis revealed a singular clustered organization of genes encoding secreted effectors. In U. maydis, many of these clusters have a role in virulence. Reconstructing the evolutionary history of clusters of effector genes is difficult because of their intrinsically fast evolution, which erodes the phylogenetic signal and homology relationships. Here, we describe the use of comparative evolutionary analyses of quality draft assemblies of genomes to study the mechanisms of this evolution. We report the genome sequence of a South African isolate of Sporisorium scitamineum, a smut fungus parasitizing sugar cane with a phylogenetic position intermediate to the two previously sequenced species U. maydis and Sporisorium reilianum. We show that the genome of S. scitamineum contains more and larger gene clusters encoding secreted effectors than any previously described species in this group. We trace back the origin of the clusters and find that their evolution is mainly driven by tandem gene duplication. In addition, transposable elements play a major role in the evolution of the clustered genes. Transposable elements are significantly associated with clusters of genes encoding fast evolving secreted effectors. This suggests that such clusters represent a case of genome compartmentalization that restrains the activity of transposable elements on genes under diversifying selection for which this activity is potentially beneficial, while protecting the rest of the genome from its deleterious effect.

Comparative Transcriptomic Analysis of Virulence Factors in Leptosphaeria maculans during Compatible and Incompatible Interactions with Canola

Leptosphaeria maculans is a hemibiotrophic fungus that causes blackleg of canola (Brassica napus), one of the most devastating diseases of this crop. In the present study, transcriptome profiling of L. maculans was performed in an effort to understand and define the pathogenicity genes that govern both the biotrophic and the necrotrophic phase of the fungus, as well as those that separate a compatible from an incompatible interaction. For this purpose, comparative RNA-seq analyses were performed on L. maculans isolate D5 at four different time points following inoculation on susceptible cultivar Topas-DH16516 or resistant introgression line Topas-Rlm2. Analysis of 1.6 billion Illumina reads readily identified differentially expressed genes that were over represented by candidate secretory effector proteins, CAZymes, and other pathogenicity genes. Comparisons between the compatible and incompatible interactions led to the identification of 28 effector proteins whose chronology and level of expression suggested a role in the establishment and maintenance of biotrophy with the plant. These included all known Avr genes of isolate D5 along with eight newly characterized effectors. In addition, another 15 effector proteins were found to be exclusively expressed during the necrotrophic phase of the fungus, which supports the concept that L. maculans has a separate and distinct arsenal contributing to each phase. As for CAZymes, they were often highly expressed at 3 dpi but with no difference in expression between the compatible and incompatible interactions, indicating that other factors were necessary to determine the outcome of the interaction. However, their significantly higher expression at 11 dpi in the compatible interaction confirmed that they contributed to the necrotrophic phase of the fungus. A notable exception was LysM genes whose high expression was singularly observed on the susceptible host at 7 dpi. In the case of TFs, their higher expression at 7 and 11 dpi on susceptible Topas support an important role in regulating the genes involved in the different pathogenic phases of L. maculans. In conclusion, comparison of the transcriptome of L. maculans during compatible and incompatible interactions has led to the identification of key pathogenicity genes that regulate not only the fate of the interaction but also lifestyle transitions of the fungus.

Molecular recognition of avirulence protein (avrxa5) by eukaryotic transcription factor xa5 of rice (Oryza sativa L.): insights from molecular dynamics simulations

The avirulence gene avrxa5 of bacterial blight pathogen Xanthomonas oryzae pv. oryzae (Xoo) recognized by the resistant rice lines having corresponding resistance (xa5) gene in a gene-for-gene manner. We used a combinatorial approach involving protein-protein docking, molecular dynamics (MD) simulations and binding free energy calculations to gain novel insights into the gene-for-gene mechanism that governs the direct interaction of R-Avr protein. From the best three binding poses predicted by molecular docking, MD simulations were performed to explore the dynamic binding mechanism of xa5 and avrxa5. Molecular Mechanics/Poisson Boltzmann Surface Area (MM/PBSA) techniques were employed to calculate the binding free energy and to uncover the thriving force behind the molecular recognition of avrxa5 by eukaryotic transcription factor xa5. Binding free energy analysis revealed van der Waals term as the most constructive component that favors the xa5 and avrxa5 interaction. In addition, hydrogen bonds (H-bonds) and essential electrostatic interactions analysis highlighted amino acid residues Lys54/Asp870, Lys56/Ala868, Lys56/Ala866, Lys56/Glu871, Ile59/His862, Gly61/Phe858, His62/Arg841, His62/Leu856, Ser101/Ala872 and Ser105/Asp870 plays pivotal role for the energetically stability of the R-Avr complex. Insights gained from the present study are expected to unveil the molecular mechanisms that define the transcriptional activator mediated transcriptome modification in host plants.