Wheat Domestication Accelerated Evolution and Triggered Positive Selection in the β-Xylosidase Enzyme of Mycosphaerella graminicola

Plant cell wall degrading enzymes (PCWDEs) of plant pathogens are receiving increasing interest for their potential to trigger plant defense reactions. In an antagonistic co-evolutionary arms race between host and pathogen, PCWDEs could be under strong selection. Here, we tested the hypothesis that PCWDEs in the fungal wheat pathogen Mycosphaerella graminicola have been positively selected by analyzing ratios of non-synonymous and synonymous nucleotide changes in the genes encoding these enzymes. Analyses of five PCWDEs demonstrated that one (β-xylosidase) has been under strong positive selection and experienced an accelerated rate of evolution. In contrast, PCWDEs in the closest relatives of M. graminicola collected from wild grasses did not show evidence for selection or deviation from a molecular clock. Since the genealogical divergence of M. graminicola from these latter species coincided with the onset of agriculture, we hypothesize that the recent domestication of the host plant and/or agricultural practices triggered positive selection in β-xylosidase and that this enzyme played a key role in the emergence of a host-specialized pathogen

Processing and analyzing multiple genomes alignments with MafFilter

As the number of available genome sequences from both closely related species and individuals withinspecies increased, theoretical and methodological convergences between the fields of phylogenomics andpopulation genomics emerged. Population genomics typically focuses on the analysis of variants, whilephylogenomics heavily relies on genome alignments. However, these are playing an increasingly importantrole in studies at the population level. Multiple genome alignments of individuals are used when structuralvariation is of primary interest and when genome architecture permits to assemblede novogenomesequences. Here I describe MafFilter, a command-line-driven program allowing to process genome align-ments in the Multiple Alignment Format (MAF). Using concrete examples based on publicly availabledatasets, I demonstrate how MafFilter can be used to develop efficient and reproducible pipelines withquality assurance for downstream analyses. I further show how MafFilter can be used to perform both basicand advanced population genomic analyses in order to infer the patterns of nucleotide diversity alonggenomes

Evolutionary relationships between Rhynchosporium lolii sp. nov. and other Rhynchosporium species on grass.

Copyright: 2013 King et al. This is an open-access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are creditedThe fungal genus Rhynchosporium (causative agent of leaf blotch) contains several host-specialised species, including R. commune (colonising barley and brome-grass), R. agropyri (couch-grass), R. secalis (rye and triticale) and the more distantly related R. orthosporum (cocksfoot). This study used molecular fingerprinting, multilocus DNA sequence data, conidial morphology, host range tests and scanning electron microscopy to investigate the relationship between Rhynchosporium species on ryegrasses, both economically important forage grasses and common wild grasses in many cereal growing areas, and other plant species. Two different types of Rhynchosporium were found on ryegrasses in the UK. Firstly, there were isolates of R. commune that were pathogenic to both barley and Italian ryegrass. Secondly, there were isolates of a new species, here named R. lolii, that were pathogenic only to ryegrass species. R. lolii was most closely related to R. orthosporum, but exhibited clear molecular, morphological and host range differences. The species was estimated to have diverged from R. orthosporum ca. 5735 years before the present. The colonisation strategy of all of the different Rhynchosporium species involved extensive hyphal growth in the sub-cuticular regions of the leaves. Finally, new species-specific PCR diagnostic tests were developed that could distinguish between these five closely related Rhynchosporium species.Peer reviewedFinal Published versio

Effect of hosts on competition among clones and evidence of differential selection between pathogenic and saprophytic phases in experimental populations of the wheat pathogen Phaeosphaeria nodorum

<p>Abstract</p> <p>Background</p> <p>Monoculture, multi-cropping and wider use of highly resistant cultivars have been proposed as mechanisms to explain the elevated rate of evolution of plant pathogens in agricultural ecosystems. We used a mark-release-recapture experiment with the wheat pathogen <it>Phaeosphaeria nodorum </it>to evaluate the impact of two of these mechanisms on the evolution of a pathogen population. Nine <it>P. nodorum </it>isolates marked with ten microsatellite markers and one minisatellite were released onto five replicated host populations to initiate epidemics of Stagonospora nodorum leaf blotch. The experiment was carried out over two consecutive host growing seasons and two pathogen collections were made during each season.</p> <p>Results</p> <p>A total of 637 pathogen isolates matching the marked inoculants were recovered from inoculated plots over two years. Genetic diversity in the host populations affected the evolution of the corresponding <it>P. nodorum </it>populations. In the cultivar mixture the relative frequencies of inoculants did not change over the course of the experiment and the pathogen exhibited a low variation in selection coefficients.</p> <p>Conclusions</p> <p>Our results support the hypothesis that increasing genetic heterogeneity in host populations may retard the rate of evolution in associated pathogen populations. Our experiment also provides indirect evidence of fitness costs associated with host specialization in <it>P. nodorum </it>as indicated by differential selection during the pathogenic and saprophytic phases.</p

Effect of Biodiversity Changes in Disease Risk: Exploring Disease Emergence in a Plant-Virus System

The effect of biodiversity on the ability of parasites to infect their host and cause disease (i.e. disease risk) is a major question in pathology, which is central to understand the emergence of infectious diseases, and to develop strategies for their management. Two hypotheses, which can be considered as extremes of a continuum, relate biodiversity to disease risk: One states that biodiversity is positively correlated with disease risk (Amplification Effect), and the second predicts a negative correlation between biodiversity and disease risk (Dilution Effect). Which of them applies better to different host-parasite systems is still a source of debate, due to limited experimental or empirical data. This is especially the case for viral diseases of plants. To address this subject, we have monitored for three years the prevalence of several viruses, and virus-associated symptoms, in populations of wild pepper (chiltepin) under different levels of human management. For each population, we also measured the habitat species diversity, host plant genetic diversity and host plant density. Results indicate that disease and infection risk increased with the level of human management, which was associated with decreased species diversity and host genetic diversity, and with increased host plant density. Importantly, species diversity of the habitat was the primary predictor of disease risk for wild chiltepin populations. This changed in managed populations where host genetic diversity was the primary predictor. Host density was generally a poorer predictor of disease and infection risk. These results support the dilution effect hypothesis, and underline the relevance of different ecological factors in determining disease/infection risk in host plant populations under different levels of anthropic influence. These results are relevant for managing plant diseases and for establishing conservation policies for endangered plant species

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

Tracing the role of human civilization in the globalization of plant pathogens

Acknowledgments. We apologize to all those colleagues whose work was not cited because of space restrictions.Peer reviewedPostprin

Whole-Genome and Chromosome Evolution Associated with Host Adaptation and Speciation of the Wheat Pathogen Mycosphaerella graminicola

The fungus Mycosphaerella graminicola has been a pathogen of wheat since host domestication 10,000–12,000 years ago in the Fertile Crescent. The wheat-infecting lineage emerged from closely related Mycosphaerella pathogens infecting wild grasses. We use a comparative genomics approach to assess how the process of host specialization affected the genome structure of M. graminicola since divergence from the closest known progenitor species named M. graminicola S1. The genome of S1 was obtained by Illumina sequencing resulting in a 35 Mb draft genome sequence of 32X. Assembled contigs were aligned to the previously sequenced M. graminicola genome. The alignment covered >90% of the non-repetitive portion of the M. graminicola genome with an average divergence of 7%. The sequenced M. graminicola strain is known to harbor thirteen essential chromosomes plus eight dispensable chromosomes. We found evidence that structural rearrangements significantly affected the dispensable chromosomes while the essential chromosomes were syntenic. At the nucleotide level, the essential and dispensable chromosomes have evolved differently. The average synonymous substitution rate in dispensable chromosomes is considerably lower than in essential chromosomes, whereas the average non-synonymous substitution rate is three times higher. Differences in molecular evolution can be related to different transmission and recombination patterns, as well as to differences in effective population sizes of essential and dispensable chromosomes. In order to identify genes potentially involved in host specialization or speciation, we calculated ratios of synonymous and non-synonymous substitution rates in the >9,500 aligned protein coding genes. The genes are generally under strong purifying selection. We identified 43 candidate genes showing evidence of positive selection, one encoding a potential pathogen effector protein. We conclude that divergence of these pathogens was accompanied by structural rearrangements in the small dispensable chromosomes, while footprints of positive selection were present in only a small number of protein coding genes

Evolution of Linked Avirulence Effectors in Leptosphaeria maculans Is Affected by Genomic Environment and Exposure to Resistance Genes in Host Plants

Brassica napus (canola) cultivars and isolates of the blackleg fungus, Leptosphaeria maculans interact in a ‘gene for gene’ manner whereby plant resistance (R) genes are complementary to pathogen avirulence (Avr) genes. Avirulence genes encode proteins that belong to a class of pathogen molecules known as effectors, which includes small secreted proteins that play a role in disease. In Australia in 2003 canola cultivars with the Rlm1 resistance gene suffered a breakdown of disease resistance, resulting in severe yield losses. This was associated with a large increase in the frequency of virulence alleles of the complementary avirulence gene, AvrLm1, in fungal populations. Surprisingly, the frequency of virulence alleles of AvrLm6 (complementary to Rlm6) also increased dramatically, even though the cultivars did not contain Rlm6. In the L. maculans genome, AvrLm1 and AvrLm6 are linked along with five other genes in a region interspersed with transposable elements that have been degenerated by Repeat-Induced Point (RIP) mutations. Analyses of 295 Australian isolates showed deletions, RIP mutations and/or non-RIP derived amino acid substitutions in the predicted proteins encoded by these seven genes. The degree of RIP mutations within single copy sequences in this region was proportional to their proximity to the degenerated transposable elements. The RIP alleles were monophyletic and were present only in isolates collected after resistance conferred by Rlm1 broke down, whereas deletion alleles belonged to several polyphyletic lineages and were present before and after the resistance breakdown. Thus, genomic environment and exposure to resistance genes in B. napus has affected the evolution of these linked avirulence genes in L. maculans

Gene expression polymorphism underpins evasion of host immunity in an asexual lineage of the Irish potato famine pathogen

BACKGROUND: Outbreaks caused by asexual lineages of fungal and oomycete pathogens are a continuing threat to crops, wild animals and natural ecosystems (Fisher MC, Henk DA, Briggs CJ, Brownstein JS, Madoff LC, McCraw SL, Gurr SJ, Nature 484:186-194, 2012; Kupferschmidt K, Science 337:636-638, 2012). However, the mechanisms underlying genome evolution and phenotypic plasticity in asexual eukaryotic microbes remain poorly understood (Seidl MF, Thomma BP, BioEssays 36:335-345, 2014). Ever since the 19th century Irish famine, the oomycete Phytophthora infestans has caused recurrent outbreaks on potato and tomato crops that have been primarily caused by the successive rise and migration of pandemic asexual lineages (Goodwin SB, Cohen BA, Fry WE, Proc Natl Acad Sci USA 91:11591-11595, 1994; Yoshida K, Burbano HA, Krause J, Thines M, Weigel D, Kamoun S, PLoS Pathog 10:e1004028, 2014; Yoshida K, Schuenemann VJ, Cano LM, Pais M, Mishra B, Sharma R, Lanz C, Martin FN, Kamoun S, Krause J, et al. eLife 2:e00731, 2013; Cooke DEL, Cano LM, Raffaele S, Bain RA, Cooke LR, Etherington GJ, Deahl KL, Farrer RA, Gilroy EM, Goss EM, et al. PLoS Pathog 8:e1002940, 2012). However, the dynamics of genome evolution within these clonal lineages have not been determined. The objective of this study was to use a comparative genomics and transcriptomics approach to determine the molecular mechanisms that underpin phenotypic variation within a clonal lineage of P. infestans. RESULTS: Here, we reveal patterns of genomic and gene expression variation within a P. infestans asexual lineage by comparing strains belonging to the South American EC-1 clone that has dominated Andean populations since the 1990s (Yoshida K, Burbano HA, Krause J, Thines M, Weigel D, Kamoun S, PLoS Pathog 10e1004028, 2014; Yoshida K, Schuenemann VJ, Cano LM, Pais M, Mishra B, Sharma R, Lanz C, Martin FN, Kamoun S, Krause J, et al. eLife 2:e00731, 2013; Delgado RA, Monteros-Altamirano AR, Li Y, Visser RGF, van der Lee TAJ, Vosman B, Plant Pathol 62:1081-1088, 2013; Forbes GA, Escobar XC, Ayala CC, Revelo J, Ordonez ME, Fry BA, Doucett K, Fry WE, Phytopathology 87:375-380, 1997; Oyarzun PJ, Pozo A, Ordonez ME, Doucett K, Forbes GA, Phytopathology 88:265-271, 1998). We detected numerous examples of structural variation, nucleotide polymorphisms and loss of heterozygosity within the EC-1 clone. Remarkably, 17 genes are not expressed in one of the two EC-1 isolates despite apparent absence of sequence polymorphisms. Among these, silencing of an effector gene was associated with evasion of disease resistance conferred by a potato immune receptor. CONCLUSIONS: Our findings highlight the molecular changes underpinning the exceptional genetic and phenotypic plasticity associated with host adaptation in a pandemic clonal lineage of a eukaryotic plant pathogen. We observed that the asexual P. infestans lineage EC-1 can exhibit phenotypic plasticity in the absence of apparent genetic mutations resulting in virulence on a potato carrying the Rpi-vnt1.1 gene. Such variant alleles may be epialleles that arose through epigenetic changes in the underlying genes