Tracing the Origin and Evolutionary History of Pyricularia oryzae Infecting Maize and Barnyard Grass

Distinct genomic contexts predict gene presence-absence variation in different pathotypes of Magnaporthe oryzae

Fungi use the accessory gene content of their pangenomes to adapt to their environments. While gene presence-absence variation contributes to shaping accessory gene reservoirs, the genomic contexts that shape these events remain unclear. Since pangenome studies are typically species-wide and do not analyze different populations separately, it is yet to be uncovered whether presence-absence variation patterns and mechanisms are consistent across populations. Fungal plant pathogens are useful models for studying presence-absence variation because they rely on it to adapt to their hosts, and members of a species often infect distinct hosts. We analyzed gene presence-absence variation in the blast fungus, Magnaporthe oryzae (syn. Pyricularia oryzae), and found that presence-absence variation genes involved in host-pathogen and microbe-microbe interactions may drive the adaptation of the fungus to its environment. We then analyzed genomic and epigenomic features of presence-absence variation and observed that proximity to transposable elements, gene GC content, gene length, expression level in the host, and histone H3K27me3 marks were different between presence-absence variation genes and conserved genes. We used these features to construct a model that was able to predict whether a gene is likely to experience presence-absence variation with high precision (86.06%) and recall (92.88%) in M. oryzae. Finally, we found that presence-absence variation genes in the rice and wheat pathotypes of M. oryzae differed in their number and their genomic context. Our results suggest that genomic and epigenomic features of gene presence-absence variation can be used to better understand and predict fungal pangenome evolution. We also show that substantial intra-species variation can exist in these features.

Pyricularia oryzae: Lab star and field scourge

Pyricularia oryzae (syn. Magnaporthe oryzae), is a filamentous ascomycete that causes a major disease called blast on cereal crops, as well as on a wide variety of wild and cultivated grasses. Blast diseases have a tremendous impact worldwide particularly on rice and on wheat, where the disease emerged in South America in the 1980s, before spreading to Asia and Africa. Its economic importance, coupled with its amenability to molecular and genetic manipulation, have inspired extensive research efforts aiming at understanding its biology and evolution. In the past 40 years, this plant-pathogenic fungus has emerged as a major model in molecular plant-microbe interactions. In this review, we focus on the clarification of the taxonomy and genetic structure of the species and its host range determinants. We also discuss recent molecular studies deciphering its lifecycle.

TAXONOMY

Kingdom: Fungi, phylum: Ascomycota, sub-phylum: Pezizomycotina, class: Sordariomycetes, order: Magnaporthales, family: Pyriculariaceae, genus: Pyricularia.

HOST RANGE

P. oryzae has the ability to infect a wide range of Poaceae. It is structured into different host-specialized lineages that are each associated with a few host plant genera. The fungus is best known to cause tremendous damage to rice crops, but it can also attack other economically important crops such as wheat, maize, barley, and finger millet.

DISEASE SYMPTOMS

P. oryzae can cause necrotic lesions or bleaching on all aerial parts of its host plants, including leaf blades, sheaths, and inflorescences (panicles, spikes, and seeds). Characteristic symptoms on leaves are diamond-shaped silver lesions that often have a brown margin and whose appearance is influenced by numerous factors such as the plant genotype and environmental conditions. USEFUL WEBSITES Resources URL Genomic data repositories http://genome.jouy.inra.fr/gemo/ Genomic data repositories http://openriceblast.org/ Genomic data repositories http://openwheatblast.net/ Genome browser for fungi (including P. oryzae) http://fungi.ensembl.org/index.html Comparative genomics database https://mycocosm.jgi.doe.gov/mycocosm/home T-DNA mutant database http://atmt.snu.kr/ T-DNA mutant database http://www.phi-base.org/ SNP and expression data https://fungidb.org/fungidb/app/.

The complete mitochondrial genome of the rice blast fungus Pyricularia oryzae Cavara 1892 strain Guy11 and phylogenetic analysis

The complete mitochondrial genome of Pyricularia oryzae Cavara 1892 strain Guy11 is 34,865 bp in length (GenBank accession number OP095391), containing 29 tRNA genes, 2 rRNA genes, and 15 protein-coding genes (PCGs). The gene order and orientation are novel compared to other Sordariomycetes species with sequenced mitogenomes in the GenBank database. Phylogenetic analysis suggests that P. oryzae Guy11 and 19 other Sordariomycetes species form a monophyletic group. The complete mitochondrial sequence of P. oryzae Guy11 will be a valuable resource for species identification, population genetics, phylogenetics, and comparative genomics studies in Sordariomycetes and Magnaporthales.

Adaptive evolution in virulence effectors of the rice blast fungus Pyricularia oryzae

Plant pathogens secrete proteins called effectors that target host cellular processes to promote disease. Recently, structural genomics has identified several families of fungal effectors that share a similar three-dimensional structure despite remarkably variable amino-acid sequences and surface properties. To explore the selective forces that underlie the sequence variability of structurally-analogous effectors, we focused on MAX effectors, a structural family of effectors that are major determinants of virulence in the rice blast fungus Pyricularia oryzae. Using structure-informed gene annotation, we identified 58 to 78 MAX effector genes per genome in a set of 120 isolates representing seven host-associated lineages. The expression of MAX effector genes was primarily restricted to the early biotrophic phase of infection and strongly influenced by the host plant. Pangenome analyses of MAX effectors demonstrated extensive presence/absence polymorphism and identified gene loss events possibly involved in host range adaptation. However, gene knock-in experiments did not reveal a strong effect on virulence phenotypes suggesting that other evolutionary mechanisms are the main drivers of MAX effector losses. MAX effectors displayed high levels of standing variation and high rates of non-synonymous substitutions, pointing to widespread positive selection shaping the molecular diversity of MAX effectors. The combination of these analyses with structural data revealed that positive selection acts mostly on residues located in particular structural elements and at specific positions. By providing a comprehensive catalog of amino acid polymorphism, and by identifying the structural determinants of the sequence diversity, our work will inform future studies aimed at elucidating the function and mode of action of MAX effectors.

Maintenance of divergent lineages of the Rice Blast Fungus Pyricularia oryzae through niche separation, loss of sex and post-mating genetic incompatibilities

Many species of fungal plant pathogens coexist as multiple lineages on the same host, but the factors underlying the origin and maintenance of population structure remain largely unknown. The rice blast fungus Pyricularia oryzae is a widespread model plant pathogen displaying population subdivision. However, most studies of natural variation in P. oryzae have been limited in genomic or geographic resolution, and host adaptation is the only factor that has been investigated extensively as a contributor to population subdivision. In an effort to complement previous studies, we analyzed genetic and phenotypic diversity in isolates of the rice blast fungus covering a broad geographical range. Using single-nucleotide polymorphism genotyping data for 886 isolates sampled from 152 sites in 51 countries, we showed that population subdivision of P. oryzae in one recombining and three clonal lineages with broad distributions persisted with deeper sampling. We also extended previous findings by showing further population subdivision of the recombining lineage into one international and three Asian clusters, and by providing evidence that the three clonal lineages of P. oryzae were found in areas with different prevailing environmental conditions, indicating niche separation. Pathogenicity tests and bioinformatic analyses using an extended set of isolates and rice varieties indicated that partial specialization to rice subgroups contributed to niche separation between lineages, and differences in repertoires of putative virulence effectors were consistent with differences in host range. Experimental crosses revealed that female sterility and early post-mating genetic incompatibilities acted as strong additional barriers to gene flow between clonal lineages. Our results demonstrate that the spread of a fungal pathogen across heterogeneous habitats and divergent populations of a crop species can lead to niche separation and reproductive isolation between distinct, widely distributed, lineages.

Population Genomic- and Phylogenomic-Enabled Advances to Increase Insight Into Pathogen Biology and Epidemiology

Population genetics has been a key discipline in phytopathology for many years. The recent rise in cost-effective, high-throughput DNA sequencing technologies, allows sequencing of dozens, if not hundreds of specimens, turning population genetics into population genomics and opening up new, exciting opportunities as described in this Focus Issue. Without the limitations of genetic markers and the availability of whole or near whole-genome data, population genomics can give new insights into the biology, evolution and adaptation, and dissemination patterns of plant-associated microbes.