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.

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

[This corrects the article DOI: 10.1371/journal.ppat.1010687.].

Population structure of Pyricularia oryzae on rice in Vietnam reveals diversified populations with four pandemic and two endemic clusters.

We characterized the genetic structure of 609 strains of Pyricularia oryzae, the fungal pathogen causing rice blast disease, in three main regions in Vietnam using microsatellites (SSR) markers. From the 447 distinct multilocus genotypes identified, six genetic clusters were defined, all of them showing elevated genetic and genotypic diversities. Four of these clusters were related to rice-attacking lineages already described at the worldwide scale, whereas the two remaining clusters were endemic to Vietnam. Strains were unevenly distributed into the six clusters depending on their groups of rice variety (indica / japonica) or type of varieties (traditional / modern) of origin, but none of the clusters was specifically related to these two factors. The highest diversity of blast population was found in Northern mountainous area, and the lowest in Red River Delta in both terms of genetic diversity and gene diversity. Hierarchical AMOVAs confirmed that all three factors considered (rice variety group, type of variety origin and geography) significantly contributed to the population structure of P. oryzae in Vietnam, with highest contribution from rice variety group. Mating types were unevenly distributed among clusters. Combined with results of female fertility and linkage disequilibirum, we hypothesized that clonal reproduction probably occurred in all clusters, but that sexual reproduction likely took place at least in some restricted areas in the Northern mountainous area for strains belonging to the cluster related to the previously described recombinant lineage (worldwide lineage 1). Our study pictures the genetic diversity, population structure and reproductive mode of the blast fungus in central and north Vietnam, and shows that the observed population structure is explained by several factors, the most important one being the variability of rice variety. All these new information might help for elaborating appropriate strategies to controlling the blast disease.

Evolution of the rice blast pathogen on spatially structured rice landraces maintains multiple generalist fungal lineages.

Traditional agrosystems, where humans, crops and microbes have coevolved over long periods, can serve as models to understand the ecoevolutionary determinants of disease dynamics and help the engineering of durably resistant agrosystems. Here, we investigated the genetic and phenotypic relationship between rice (Oryza sativa) landraces and their rice blast pathogen (Pyricularia oryzae) in the traditional Yuanyang terraces of flooded rice paddies in China, where rice landraces have been grown and bred over centuries without significant disease outbreaks. Analyses of genetic subdivision revealed that indica rice plants clustered according to landrace names. Three new diverse lineages of rice blast specific to the Yuanyang terraces coexisted with lineages previously detected at the worldwide scale. Population subdivision in the pathogen population did not mirror pattern of population subdivision in the host. Measuring the pathogenicity of rice blast isolates on landraces revealed generalist life history traits. Our results suggest that the implementation of disease control strategies based on the emergence or maintenance of a generalist lifestyle in pathogens may sustainably reduce the burden of disease in crops.

Male fertility in Pyricularia oryzae: Microconidia are spermatia.

Sexual reproduction in Ascomycetes is well described in several model organisms such as Neurospora crassa or Podospora anserina. Deciphering the biological process of sexual reproduction (from the recognition between compatible partners to the formation of zygote) can be a major advantage to better control sexually reproducing pathogenic fungi. In Pyricularia oryzae, the fungal pathogen causing blast diseases on several Poaceae species, the biology of sexual reproduction remains poorly documented. Besides the well-documented production of asexual macroconidia, the production of microconidia was seldom reported in P. oryzae, and their role as male gamete (i.e., spermatia) and in male fertility has never been explored. Here, we characterised the morphological features of microconidia and demonstrated that they are bona fide spermatia. Contrary to macroconidia, microconidia are not able to germinate and seem to be the only male gametes in P. oryzae. We show that fruiting body (perithecium) formation requires microconidia to get in contact with mycelium of strains of opposite mating type, to presumably fertilise the female gametes.

Congruent population genetic structures and divergence histories in anther-smut fungi and their host plants Silene italica and the Silene nutans species complex.

Study of the congruence of population genetic structure between hosts and pathogens gives important insights into their shared phylogeographical and coevolutionary histories. We studied the population genetic structure of castrating anther-smut fungi (genus Microbotryum) and of their host plants, the Silene nutans species complex, and the morphologically and genetically closely related Silene italica, which can be found in sympatry. Phylogeographical population genetic structure related to persistence in separate glacial refugia has been recently revealed in the S. nutans plant species complex across Western Europe, identifying several distinct lineages. We genotyped 171 associated plant-pathogen pairs of anther-smut fungi and their host plant individuals using microsatellite markers and plant chloroplastic single nucleotide polymorphisms. We found clear differentiation between fungal populations parasitizing S. nutans and S. italica plants. The population genetic structure of fungal strains parasitizing the S. nutans plant species complex mirrored the host plant genetic structure, suggesting that the pathogen was isolated in glacial refugia together with its host and/or that it has specialized on the plant genetic lineages. Using random forest approximate Bayesian computation (ABC-RF), we found that the divergence history of the fungal lineages on S. nutans was congruent with that previously inferred for the host plant and probably occurred with ancient but no recent gene flow. Genome sequences confirmed the genetic structure and the absence of recent gene flow between fungal genetic lineages. Our analyses of individual host-pathogen pairs contribute to a better understanding of co-evolutionary histories between hosts and pathogens in natural ecosystems, in which such studies remain scarce.

A Genomic Approach to Develop a New qPCR Test Enabling Detection of the Pyricularia oryzae Lineage Causing Wheat Blast.

Rapid detection is key to managing emerging diseases because it allows their spread around the world to be monitored and limited. The first major wheat blast epidemics were reported in 1985 in the Brazilian state of Paraná. Following this outbreak, the disease quickly spread to neighboring regions and countries and, in 2016, the first report of wheat blast disease outside South America was released. This Asian outbreak was due to the trade of infected South American seed, demonstrating the importance of detection tests in order to avoid importing contaminated biological material into regions free from the pathogen. Genomic analysis has revealed that one particular lineage within the fungal species Pyricularia oryzae is associated with this disease: the Triticum lineage. A comparison of 81 Pyricularia genomes highlighted polymorphisms specific to the Triticum lineage, and this study developed a real-time PCR test targeting one of these polymorphisms. The test's performance was then evaluated in order to measure its analytical specificity, analytical sensitivity, and robustness. The C17 quantitative PCR test detected isolates belonging to the Triticum lineage with high sensitivity, down to 13 plasmid copies or 1 pg of genomic DNA per reaction tube. The blast-based approach developed here to study P. oryzae can be transposed to other emerging diseases.

Structure Analysis Uncovers a Highly Diverse but Structurally Conserved Effector Family in Phytopathogenic Fungi.

Phytopathogenic ascomycete fungi possess huge effector repertoires that are dominated by hundreds of sequence-unrelated small secreted proteins. The molecular function of these effectors and the evolutionary mechanisms that generate this tremendous number of singleton genes are largely unknown. To get a deeper understanding of fungal effectors, we determined by NMR spectroscopy the 3-dimensional structures of the Magnaporthe oryzae effectors AVR1-CO39 and AVR-Pia. Despite a lack of sequence similarity, both proteins have very similar 6 ß-sandwich structures that are stabilized in both cases by a disulfide bridge between 2 conserved cysteins located in similar positions of the proteins. Structural similarity searches revealed that AvrPiz-t, another effector from M. oryzae, and ToxB, an effector of the wheat tan spot pathogen Pyrenophora tritici-repentis have the same structures suggesting the existence of a family of sequence-unrelated but structurally conserved fungal effectors that we named MAX-effectors (Magnaporthe Avrs and ToxB like). Structure-informed pattern searches strengthened this hypothesis by identifying MAX-effector candidates in a broad range of ascomycete phytopathogens. Strong expansion of the MAX-effector family was detected in M. oryzae and M. grisea where they seem to be particularly important since they account for 5-10% of the effector repertoire and 50% of the cloned avirulence effectors. Expression analysis indicated that the majority of M. oryzae MAX-effectors are expressed specifically during early infection suggesting important functions during biotrophic host colonization. We hypothesize that the scenario observed for MAX-effectors can serve as a paradigm for ascomycete effector diversity and that the enormous number of sequence-unrelated ascomycete effectors may in fact belong to a restricted set of structurally conserved effector families.

Testing Differences Between Pathogen Compositions with Small Samples and Sparse Data.

The structure of pathogen populations is an important driver of epidemics affecting crops and natural plant communities. Comparing the composition of two pathogen populations consisting of assemblages of genotypes or phenotypes is a crucial, recurrent question encountered in many studies in plant disease epidemiology. Determining whether there is a significant difference between two sets of proportions is also a generic question for numerous biological fields. When samples are small and data are sparse, it is not straightforward to provide an accurate answer to this simple question because routine statistical tests may not be exactly calibrated. To tackle this issue, we built a computationally intensive testing procedure, the generalized Monte Carlo plug-in test with calibration test, which is implemented in an R package available at https://doi.org/10.5281/zenodo.635791 . A simulation study was carried out to assess the performance of the proposed methodology and to make a comparison with standard statistical tests. This study allows us to give advice on how to apply the proposed method, depending on the sample sizes. The proposed methodology was then applied to real datasets and the results of the analyses were discussed from an epidemiological perspective. The applications to real data sets deal with three topics in plant pathology: the reproduction of Magnaporthe oryzae, the spatial structure of Pseudomonas syringae, and the temporal recurrence of Puccinia triticina.

Emergence of wheat blast in Bangladesh was caused by a South American lineage of Magnaporthe oryzae.

BACKGROUND: In February 2016, a new fungal disease was spotted in wheat fields across eight districts in Bangladesh. The epidemic spread to an estimated 15,000 hectares, about 16 % of the cultivated wheat area in Bangladesh, with yield losses reaching up to 100 %. Within weeks of the onset of the epidemic, we performed transcriptome sequencing of symptomatic leaf samples collected directly from Bangladeshi fields. RESULTS: Reinoculation of seedlings with strains isolated from infected wheat grains showed wheat blast symptoms on leaves of wheat but not rice. Our phylogenomic and population genomic analyses revealed that the wheat blast outbreak in Bangladesh was most likely caused by a wheat-infecting South American lineage of the blast fungus Magnaporthe oryzae. CONCLUSION: Our findings suggest that genomic surveillance can be rapidly applied to monitor plant disease outbreaks and provide valuable information regarding the identity and origin of the infectious agent.

The rice resistance protein pair RGA4/RGA5 recognizes the Magnaporthe oryzae effectors AVR-Pia and AVR1-CO39 by direct binding.

Resistance (R) proteins recognize pathogen avirulence (Avr) proteins by direct or indirect binding and are multidomain proteins generally carrying a nucleotide binding (NB) and a leucine-rich repeat (LRR) domain. Two NB-LRR protein-coding genes from rice (Oryza sativa), RGA4 and RGA5, were found to be required for the recognition of the Magnaporthe oryzae effector AVR1-CO39. RGA4 and RGA5 also mediate recognition of the unrelated M. oryzae effector AVR-Pia, indicating that the corresponding R proteins possess dual recognition specificity. For RGA5, two alternative transcripts, RGA5-A and RGA5-B, were identified. Genetic analysis showed that only RGA5-A confers resistance, while RGA5-B is inactive. Yeast two-hybrid, coimmunoprecipitation, and fluorescence resonance energy transfer-fluorescence lifetime imaging experiments revealed direct binding of AVR-Pia and AVR1-CO39 to RGA5-A, providing evidence for the recognition of multiple Avr proteins by direct binding to a single R protein. Direct binding seems to be required for resistance as an inactive AVR-Pia allele did not bind RGA5-A. A small Avr interaction domain with homology to the Avr recognition domain in the rice R protein Pik-1 was identified in the C terminus of RGA5-A. This reveals a mode of Avr protein recognition through direct binding to a novel, non-LRR interaction domain.

Evolution of Compatibility Range in the Rice-Magnaporthe oryzae System: An Uneven Distribution of R Genes Between Rice Subspecies.

Efficient strategies for limiting the impact of pathogens on crops require a good understanding of the factors underlying the evolution of compatibility range for the pathogens and host plants, i.e., the set of host genotypes that a particular pathogen genotype can infect and the set of pathogen genotypes that can infect a particular host genotype. Until now, little is known about the evolutionary and ecological factors driving compatibility ranges in systems implicating crop plants. We studied the evolution of host and pathogen compatibility ranges for rice blast disease, which is caused by the ascomycete Magnaporthe oryzae. We challenged 61 rice varieties from three rice subspecies with 31 strains of M. oryzae collected worldwide from all major known genetic groups. We determined the compatibility range of each plant variety and pathogen genotype and the severity of each plant-pathogen interaction. Compatibility ranges differed between rice subspecies, with the most resistant subspecies selecting for pathogens with broader compatibility ranges and the least resistant subspecies selecting for pathogens with narrower compatibility ranges. These results are consistent with a nested distribution of R genes between rice subspecies.

Population structure and temporal maintenance of the multihost fungal pathogen Botrytis cinerea: causes and implications for disease management.

Understanding the causes of population subdivision is of fundamental importance, as studying barriers to gene flow between populations may reveal key aspects of the process of adaptive divergence and, for pathogens, may help forecasting disease emergence and implementing sound management strategies. Here, we investigated population subdivision in the multihost fungus Botrytis cinerea based on comprehensive multiyear sampling on different hosts in three French regions. Analyses revealed a weak association between population structure and geography, but a clear differentiation according to the host plant of origin. This was consistent with adaptation to hosts, but the distribution of inferred genetic clusters and the frequency of admixed individuals indicated a lack of strict host specificity. Differentiation between individuals collected in the greenhouse (on Solanum) and outdoor (on Vitis and Rubus) was stronger than that observed between individuals from the two outdoor hosts, probably reflecting an additional isolating effect associated with the cropping system. Three genetic clusters coexisted on Vitis but did not persist over time. Linkage disequilibrium analysis indicated that outdoor populations were regularly recombining, whereas clonality was predominant in the greenhouse. Our findings open up new perspectives for disease control by managing plant debris in outdoor conditions and reinforcing prophylactic measures indoor.

South-East Asia is the center of origin, diversity and dispersion of the rice blast fungus, Magnaporthe oryzae.

• Inferring invasion routes and identifying reservoirs of diversity of plant pathogens are essential in proposing new strategies for their control. Magnaporthe oryzae, the fungus responsible for rice blast disease, has invaded all rice growing areas. Virulent genotypes regularly (re)emerge, causing rapid resistance breakdowns. However, the world-wide genetic subdivision of M. oryzae populations on rice and its past history of invasion have never been elucidated. • In order to investigate the centers of diversity, origin and migration of M. oryzae on rice, we analyzed the genetic diversity of 55 populations from 15 countries. • Three genetic clusters were identified world-wide. Asia was the center of diversity and the origin of most migrations to other continents. In Asia, two centers of diversity were revealed in the Himalayan foothills: South China-Laos-North Thailand, and western Nepal. Sexual reproduction persisted only in the South China-Laos-North Thailand region, which was identified as the putative center of origin of all M. oryzae populations on rice. • Our results suggest a scenario of early evolution of M. oryzae on rice that matches the past history of rice domestication. This study confirms that crop domestication may have considerable influence on the pestification process of natural enemies.

Genomic analysis of the necrotrophic fungal pathogens Sclerotinia sclerotiorum and Botrytis cinerea.

Sclerotinia sclerotiorum and Botrytis cinerea are closely related necrotrophic plant pathogenic fungi notable for their wide host ranges and environmental persistence. These attributes have made these species models for understanding the complexity of necrotrophic, broad host-range pathogenicity. Despite their similarities, the two species differ in mating behaviour and the ability to produce asexual spores. We have sequenced the genomes of one strain of S. sclerotiorum and two strains of B. cinerea. The comparative analysis of these genomes relative to one another and to other sequenced fungal genomes is provided here. Their 38-39 Mb genomes include 11,860-14,270 predicted genes, which share 83% amino acid identity on average between the two species. We have mapped the S. sclerotiorum assembly to 16 chromosomes and found large-scale co-linearity with the B. cinerea genomes. Seven percent of the S. sclerotiorum genome comprises transposable elements compared to <1% of B. cinerea. The arsenal of genes associated with necrotrophic processes is similar between the species, including genes involved in plant cell wall degradation and oxalic acid production. Analysis of secondary metabolism gene clusters revealed an expansion in number and diversity of B. cinerea-specific secondary metabolites relative to S. sclerotiorum. The potential diversity in secondary metabolism might be involved in adaptation to specific ecological niches. Comparative genome analysis revealed the basis of differing sexual mating compatibility systems between S. sclerotiorum and B. cinerea. The organization of the mating-type loci differs, and their structures provide evidence for the evolution of heterothallism from homothallism. These data shed light on the evolutionary and mechanistic bases of the genetically complex traits of necrotrophic pathogenicity and sexual mating. This resource should facilitate the functional studies designed to better understand what makes these fungi such successful and persistent pathogens of agronomic crops.

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 unconventional resistance protein PTR recognizes the Magnaporthe oryzae effector AVR-Pita in an allele-specific manner.

Blast disease caused by the fungus Magnaporthe oryzae is one of the most devastating rice diseases. Disease resistance genes such as Pi-ta or Pi-ta2 are critical in protecting rice production from blast. Published work reports that Pi-ta codes for a nucleotide-binding and leucine-rich repeat domain protein (NLR) that recognizes the fungal protease-like effector AVR-Pita by direct binding. However, this model was challenged by the recent discovery that Pi-ta2 resistance, which also relies on AVR-Pita detection, is conferred by the unconventional resistance gene Ptr, which codes for a membrane protein with a cytoplasmic armadillo repeat domain. Here, using NLR Pi-ta and Ptr RNAi knockdown and CRISPR/Cas9 knockout mutant rice lines, we found that AVR-Pita recognition relies solely on Ptr and that the NLR Pi-ta has no role in it, indicating that it is not the Pi-ta resistance gene. Different alleles of Ptr confer different recognition specificities. The A allele of Ptr (PtrA) detects all natural sequence variants of the effector and confers Pi-ta2 resistance, while the B allele of Ptr (PtrB) recognizes a restricted set of AVR-Pita alleles and, thereby, confers Pi-ta resistance. Analysis of the natural diversity in AVR-Pita and of mutant and transgenic strains identified one specific polymorphism in the effector sequence that controls escape from PtrB-mediated resistance. Taken together, our work establishes that the M. oryzae effector AVR-Pita is detected in an allele-specific manner by the unconventional rice resistance protein Ptr and that the NLR Pi-ta has no function in Pi-ta resistance and the recognition of AVR-Pita.

Deep population structure linked to host vernalization requirement in the barley net blotch fungal pathogen.

Invasive fungal pathogens pose a substantial threat to widely cultivated crop species, owing to their capacity to adapt to new hosts and new environmental conditions. Gaining insights into the demographic history of these pathogens and unravelling the mechanisms driving coevolutionary processes are crucial for developing durably effective disease management programmes. Pyrenophora teres is a significant fungal pathogen of barley, consisting of two lineages, Ptt and Ptm, with global distributions and demographic histories reflecting barley domestication and spread. However, the factors influencing the population structure of P. teres remain poorly understood, despite the varietal and environmental heterogeneity of barley agrosystems. Here, we report on the population genomic structure of P. teres in France and globally. We used genotyping-by-sequencing to show that Ptt and Ptm can coexist in the same area in France, with Ptt predominating. Furthermore, we showed that differences in the vernalization requirement of barley varieties were associated with population differentiation within Ptt in France and at a global scale, with one population cluster found on spring barley and another population cluster found on winter barley. Our results demonstrate how cultivation conditions, possibly associated with genetic differences between host populations, can be associated with the maintenance of divergent invasive pathogen populations coexisting over large geographic areas. This study not only advances our understanding of the coevolutionary dynamics of the Pt-barley pathosystem but also prompts further research on the relative contributions of adaptation to the host versus adaptation to abiotic conditions in shaping Ptt populations.

Arms race co-evolution of Magnaporthe oryzae AVR-Pik and rice Pik genes driven by their physical interactions.

Attack and counter-attack impose strong reciprocal selection on pathogens and hosts, leading to development of arms race evolutionary dynamics. Here we show that Magnaporthe oryzae avirulence gene AVR-Pik and the cognate rice resistance (R) gene Pik are highly variable, with multiple alleles in which DNA replacements cause amino acid changes. There is tight recognition specificity of the AVR-Pik alleles by the various Pik alleles. We found that AVR-Pik physically binds the N-terminal coiled-coil domain of Pik in a yeast two-hybrid assay as well as in an in planta co-immunoprecipitation assay. This binding specificity correlates with the recognition specificity between AVR and R genes. We propose that AVR-Pik and Pik are locked into arms race co-evolution driven by their direct physical interactions.

Asexual reproduction induces a rapid and permanent loss of sexual reproduction capacity in the rice fungal pathogen Magnaporthe oryzae: results of in vitro experimental evolution assays.

BACKGROUND: Sexual reproduction is common in eukaryotic microorganisms, with few species reproducing exclusively asexually. However, in some organisms, such as fungi, asexual reproduction alternates with episodic sexual reproduction events. Fungi are thus appropriate organisms for studies of the reasons for the selection of sexuality or clonality and of the mechanisms underlying this selection. Magnaporthe oryzae, an Ascomycete causing blast disease on rice, reproduces mostly asexually in natura. Sexual reproduction is possible in vitro and requires (i) two strains of opposite mating types including (ii) at least one female-fertile strain (i.e. a strain able to produce perithecia, the female organs in which meiosis occurs). Female-fertile strains are found only in limited areas of Asia, in which evidence for contemporary recombination has recently been obtained. We induced the forced evolution of four Chinese female-fertile strains in vitro by the weekly transfer of asexual spores (conidia) between Petri dishes. We aimed to determine whether female fertility was rapidly lost in the absence of sexual reproduction and whether this loss was controlled genetically or epigenetically. RESULTS: All the strains became female-sterile after 10 to 19 rounds of selection under asexual conditions. As no single-spore isolation was carried out, the observed decrease in the production of perithecia reflected the emergence and the invasion of female-sterile mutants. The female-sterile phenotype segregated in the offspring of crosses between female-sterile evolved strains and female-fertile wild-type strains. This segregation was maintained in the second generation in backcrosses. Female-sterile evolved strains were subjected to several stresses, but none induced the restoration of female fertility. This loss of fertility was therefore probably due to genetic rather than epigenetic mechanisms. In competition experiments, female-sterile mutants produced similar numbers of viable conidia to wild-type strains, but released them more efficiently. This advantage may account for the invasion of our populations by female-sterile mutants. CONCLUSIONS: We show for the first time that, in the absence of sexual reproduction, female-sterile mutants of M. oryzae rice strains can arise and increase in abundance in asexual generations. This change in phenotype was frequent and probably caused by mutation. These results suggest that female fertility may have been lost rapidly during the dispersion of the fungus from Asia to the rest of the world.