Gene Flow between Divergent Cereal- and Grass-Specific Lineages of the Rice Blast Fungus Magnaporthe oryzae

ABSTRACT Delineating species and epidemic lineages in fungal plant pathogens is critical to our understanding of disease emergence and the structure of fungal biodiversity and also informs international regulatory decisions. Pyricularia oryzae (syn. Magnaporthe oryzae) is a multihost pathogen that infects multiple grasses and cereals, is responsible for the most damaging rice disease (rice blast), and is of growing concern due to the recent introduction of wheat blast to Bangladesh from South America. However, the genetic structure and evolutionary history of M. oryzae, including the possible existence of cryptic phylogenetic species, remain poorly defined. Here, we use whole-genome sequence information for 76 M. oryzae isolates sampled from 12 grass and cereal genera to infer the population structure of M. oryzae and to reassess the species status of wheat-infecting populations of the fungus. Species recognition based on genealogical concordance, using published data or extracting previously used loci from genome assemblies, failed to confirm a prior assignment of wheat blast isolates to a new species (Pyricularia graminis-tritici). Inference of population subdivisions revealed multiple divergent lineages within M. oryzae, each preferentially associated with one host genus, suggesting incipient speciation following host shift or host range expansion. Analyses of gene flow, taking into account the possibility of incomplete lineage sorting, revealed that genetic exchanges have contributed to the makeup of multiple lineages within M. oryzae. These findings provide greater understanding of the ecoevolutionary factors that underlie the diversification of M. oryzae and highlight the practicality of genomic data for epidemiological surveillance in this important multihost pathogen. IMPORTANCE Infection of novel hosts is a major route for disease emergence by pathogenic microorganisms. Understanding the evolutionary history of multihost pathogens is therefore important to better predict the likely spread and emergence of new diseases. Magnaporthe oryzae is a multihost fungus that causes serious cereal diseases, including the devastating rice blast disease and wheat blast, a cause of growing concern due to its recent spread from South America to Asia. Using whole-genome analysis of 76 fungal strains from different hosts, we have documented the divergence of M. oryzae into numerous lineages, each infecting a limited number of host species. Our analyses provide evidence that interlineage gene flow has contributed to the genetic makeup of multiple M. oryzae lineages within the same species. Plant health surveillance is therefore warranted to safeguard against disease emergence in regions where multiple lineages of the fungus are in contact with one another. IMPORTANCE Infection of novel hosts is a major route for disease emergence by pathogenic microorganisms. Understanding the evolutionary history of multihost pathogens is therefore important to better predict the likely spread and emergence of new diseases. Magnaporthe oryzae is a multihost fungus that causes serious cereal diseases, including the devastating rice blast disease and wheat blast, a cause of growing concern due to its recent spread from South America to Asia. Using whole-genome analysis of 76 fungal strains from different hosts, we have documented the divergence of M. oryzae into numerous lineages, each infecting a limited number of host species. Our analyses provide evidence that interlineage gene flow has contributed to the genetic makeup of multiple M. oryzae lineages within the same species. Plant health surveillance is therefore warranted to safeguard against disease emergence in regions where multiple lineages of the fungus are in contact with one another.

[1]  D. J. Funk,et al.  REPRODUCTIVE ISOLATION CAUSED BY NATURAL SELECTION AGAINST IMMIGRANTS FROM DIVERGENT HABITATS , 2005 .

[2]  A. Narechania,et al.  Architecture of a Species: Phylogenomics of Staphylococcus aureus. , 2017, Trends in microbiology.

[3]  R. D. de Vries,et al.  Resolving the polyphyletic nature of Pyricularia (Pyriculariaceae) , 2014, Studies in mycology.

[4]  N. Talbot,et al.  Emergence of wheat blast in Bangladesh was caused by a South American lineage of Magnaporthe oryzae , 2016, BMC Biology.

[5]  Liyuan Liu,et al.  Comparative Analysis of the Genomes of Two Field Isolates of the Rice Blast Fungus Magnaporthe oryzae , 2012, PLoS genetics.

[6]  M. Farman Pyricularia grisea Isolates Causing Gray Leaf Spot on Perennial Ryegrass (Lolium perenne) in the United States: Relationship to P. grisea Isolates from Other Host Plants. , 2002, Phytopathology.

[7]  P. Smouse,et al.  genalex 6: genetic analysis in Excel. Population genetic software for teaching and research , 2006 .

[8]  Thomas L. Madden,et al.  Gapped BLAST and PSI-BLAST: a new generation of protein database search programs. , 1997, Nucleic acids research.

[9]  S. Gavrilets,et al.  Linking the emergence of fungal plant diseases with ecological speciation. , 2010, Trends in Ecology & Evolution.

[10]  A. Auton,et al.  Recombination rate estimation in the presence of hotspots. , 2007, Genome research.

[11]  Frédéric Delsuc,et al.  MACSE: Multiple Alignment of Coding SEquences Accounting for Frameshifts and Stop Codons , 2011, PloS one.

[12]  Åke Brännström,et al.  Genomics and the origin of species , 2014, Nature Reviews Genetics.

[13]  Tandy J. Warnow,et al.  ASTRAL-II: coalescent-based species tree estimation with many hundreds of taxa and thousands of genes , 2015, Bioinform..

[14]  Daniel H. Huson,et al.  SplitsTree: analyzing and visualizing evolutionary data , 1998, Bioinform..

[15]  B. McDonald,et al.  Pyricularia graminis-tritici, a new Pyricularia species causing wheat blast , 2016, Persoonia.

[16]  Eric Bazin,et al.  Sex at the origin: an Asian population of the rice blast fungus Magnaporthe oryzae reproduces sexually , 2012, Molecular ecology.

[17]  M. Siol,et al.  EggLib: processing, analysis and simulation tools for population genetics and genomics , 2012, BMC Genetics.

[18]  M. Stephens,et al.  Inference of population structure using multilocus genotype data: dominant markers and null alleles , 2007, Molecular ecology notes.

[19]  Siavash Mirarab,et al.  Fast Coalescent-Based Computation of Local Branch Support from Quartet Frequencies , 2016, Molecular biology and evolution.

[20]  K. Tanaka,et al.  Genetic Analysis of Host Species Specificity of Magnaporthe oryzae Isolates from Rice and Wheat. , 2006, Phytopathology.

[21]  B. McDonald,et al.  Wheat blast disease caused by Pyricularia graminis-tritici sp. nov. , 2016 .

[22]  M. Rep,et al.  Mitochondrial genomes reveal recombination in the presumed asexual Fusarium oxysporum species complex , 2017, BMC Genomics.

[23]  Korbinian Strimmer,et al.  APE: Analyses of Phylogenetics and Evolution in R language , 2004, Bioinform..

[24]  Simon H. Martin,et al.  Genome-wide evidence for speciation with gene flow in Heliconius butterflies , 2013, Genome research.

[25]  Scott V Edwards,et al.  Implementing and testing the multispecies coalescent model: A valuable paradigm for phylogenomics. , 2016, Molecular phylogenetics and evolution.

[26]  D. Falush,et al.  Inference of Population Structure using Dense Haplotype Data , 2012, PLoS genetics.

[27]  F. Chumley,et al.  Electrophoretic karyotypes of Magnaporthe grisea pathogens of diverse grasses , 1996 .

[28]  R. Terauchi,et al.  Host specialization of the blast fungus Magnaporthe oryzae is associated with dynamic gain and loss of genes linked to transposable elements , 2016, BMC Genomics.

[29]  A. S. Urashima,et al.  Host range, mating type, and fertility of Pyricularia grisea from wheat in Brazil , 1993 .

[30]  M. Kohli,et al.  Pyricularia blast - a threat to wheat cultivation. , 2018 .

[31]  N. Nga,et al.  Genetic analysis of host–pathogen incompatibility between Lolium isolates of Pyriculariaoryzae and wheat , 2013, Journal of General Plant Pathology.

[32]  B. Valent,et al.  Wheat blast disease: danger on the move , 2017, Tropical Plant Pathology.

[33]  J. Nottéghem,et al.  Production of perithecia of Magnaporthe grisea on rice plants , 1990 .

[34]  T. Giraud,et al.  PHYLOGENETIC EVIDENCE OF HOST-SPECIFIC CRYPTIC SPECIES IN THE ANTHER SMUT FUNGUS , 2007, Evolution; international journal of organic evolution.

[35]  P. Gladieux,et al.  Emergence of new virulent populations of apple scab from nonagricultural disease reservoirs. , 2016, The New phytologist.

[36]  T. P. Tiwari,et al.  First Report of Wheat Blast Caused by Magnaporthe oryzae Pathotype triticum in Bangladesh , 2016 .

[37]  D. Hibbett,et al.  Phylogenetic species recognition and species concepts in fungi. , 2000, Fungal genetics and biology : FG & B.

[38]  H. Nakayashiki,et al.  Pathogenicity, Mating Ability and DNA Restriction Fragment Length Polymorphisms of Pyricularia Populations Isolated from Gramineae, Bambusideae and Zingiberaceae Plants , 2000, Journal of General Plant Pathology.

[39]  A MULTILOCUS GENEALOGICAL APPROACH TO PHYLOGENETIC SPECIES RECOGNITION IN THE MODEL EUKARYOTE NEUROSPORA , 2003, Evolution; international journal of organic evolution.

[40]  Tatiana Giraud,et al.  Partition of the Botrytis cinerea complex in France using multiple gene genealogies , 2005, Mycologia.

[41]  David S. Hibbett,et al.  Fungal systematics: is a new age of enlightenment at hand? , 2013, Nature Reviews Microbiology.

[42]  David Reich,et al.  Testing for ancient admixture between closely related populations. , 2011, Molecular biology and evolution.

[43]  Marc-Henri Lebrun,et al.  Origins of Host-Specific Populations of the Blast Pathogen Magnaporthe oryzae in Crop Domestication With Subsequent Expansion of Pandemic Clones on Rice and Weeds of Rice , 2005, Genetics.

[44]  Sonja J. Prohaska,et al.  Proteinortho: Detection of (Co-)orthologs in large-scale analysis , 2011, BMC Bioinformatics.

[45]  T. Giraud,et al.  Importance of the life cycle in sympatric host race formation and speciation of pathogens. , 2006, Phytopathology.

[46]  E. Birney,et al.  Velvet: algorithms for de novo short read assembly using de Bruijn graphs. , 2008, Genome research.

[47]  Alexandros Stamatakis,et al.  RAxML version 8: a tool for phylogenetic analysis and post-analysis of large phylogenies , 2014, Bioinform..

[48]  Haibao Tang,et al.  Directional Selection from Host Plants Is a Major Force Driving Host Specificity in Magnaporthe Species , 2016, Scientific Reports.

[49]  R. Terauchi,et al.  Evolution of the wheat blast fungus through functional losses in a host specificity determinant , 2017, Science.

[50]  B. Henrissat,et al.  Deciphering Genome Content and Evolutionary Relationships of Isolates from the Fungus Magnaporthe oryzae Attacking Different Host Plants , 2015, Genome biology and evolution.

[51]  M. Stephens,et al.  Inference of population structure using multilocus genotype data: linked loci and correlated allele frequencies. , 2003, Genetics.

[52]  Cathryn J. Rehmeyer,et al.  The genome sequence of the rice blast fungus Magnaporthe grisea , 2005, Nature.

[53]  H. Matsumura,et al.  Association Genetics Reveals Three Novel Avirulence Genes from the Rice Blast Fungal Pathogen Magnaporthe oryzae[W][OA] , 2009, The Plant Cell Online.

[54]  Martin Kollmar,et al.  A novel hybrid gene prediction method employing protein multiple sequence alignments , 2011, Bioinform..

[55]  D. J. Funk,et al.  REPRODUCTIVE ISOLATION CAUSED BY NATURAL SELECTION AGAINST IMMIGRANTS FROM DIVERGENT HABITATS , 2005 .

[56]  Tandy J. Warnow,et al.  ASTRAL: genome-scale coalescent-based species tree estimation , 2014, Bioinform..

[57]  M. Stephens,et al.  Inferring weak population structure with the assistance of sample group information , 2009, Molecular ecology resources.

[58]  Simon H. Martin,et al.  Evaluating the Use of ABBA–BABA Statistics to Locate Introgressed Loci , 2014, bioRxiv.

[59]  F. Balloux,et al.  Discriminant analysis of principal components: a new method for the analysis of genetically structured populations , 2010, BMC Genetics.

[60]  B. McDonald,et al.  Population structure and pathotype diversity of the wheat blast pathogen Magnaporthe oryzae 25 years after its emergence in Brazil. , 2014, Phytopathology.

[61]  D. Hershman,et al.  The Lolium Pathotype of Magnaporthe oryzae Recovered from a Single Blasted Wheat Plant in the United States. , 2017, Plant disease.

[62]  M. Fisher,et al.  Fungal multilocus sequence typing--it's not just for bacteria. , 2003, Current opinion in microbiology.

[63]  K. Pedley,et al.  Genomics-Based Marker Discovery and Diagnostic Assay Development for Wheat Blast. , 2017, Plant disease.

[64]  M. Fisher,et al.  The rise and rise of emerging infectious fungi challenges food security and ecosystem health , 2011 .

[65]  D. Huson,et al.  Application of phylogenetic networks in evolutionary studies. , 2006, Molecular biology and evolution.

[66]  Y. Tosa,et al.  A Gene-for-Gene Relationship Underlying the Species-Specific Parasitism of Avena/Triticum Isolates of Magnaporthe grisea on Wheat Cultivars. , 2002, Phytopathology.