Redefining genera of cereal pathogens: Oculimacula, Rhynchosporium and Spermospora.

The taxonomy of Oculimacula, Rhynchosporium and Spermospora is re-evaluated, along with that of phylogenetically related genera. Isolates are identified using comparisons of DNA sequences of the internal transcribed spacer ribosomal RNA locus (ITS), partial translation elongation factor 1-alpha (tef1), actin (act), DNA-directed RNA polymerase II largest (rpb1) and second largest subunit (rpb2) genes, and the nuclear ribosomal large subunit (LSU), combined with their morphological characteristics. Oculimacula is restricted to two species, O. acuformis and O. yallundae, with O. aestiva placed in Cyphellophora, and O. anguioides accommodated in a new genus, Helgardiomyces. Rhynchosporium s. str. is restricted to species with 1-septate conidia and hooked apical beaks, while Rhynchobrunnera is introduced for species with 1-3-septate, straight conidia, lacking any apical beak. Rhynchosporium graminicola is proposed to replace the name R. commune applied to the barley scald pathogen based on nomenclatural priority. Spermospora is shown to be paraphyletic, representing Spermospora (type: S. subulata), with three new species, S. arrhenatheri, S. loliiphila and S. zeae, and Neospermospora gen. nov. (type: N. avenae). Ypsilina (type: Y. graminea), is shown to be monophyletic, but appears to be of minor importance on cereals. Finally, Vanderaaea gen. nov. (type: V. ammophilae), is introduced as a new coelomycetous fungus occurring on dead leaves of Ammophila arenaria. Citation: Crous PW, Braun U, McDonald BA, Lennox CL, Edwards J, Mann RC, Zaveri A, Linde CC, Dyer PS, Groenewald JZ (2020). Redefining genera of cereal pathogens: Oculimacula, Rhynchosporium and Spermospora. Fungal Systematics and Evolution 7: 67-98. doi: 10.3114/fuse.2021.07.04.

Genome-scale phylogenies reveal relationships among Parastagonospora species infecting domesticated and wild grasses.

Several plant pathogenic Parastagonospora species have been identified infecting wheat and other cereals over the past 50 years. As new lineages were discovered, naming conventions grew unwieldy and the relationships with previously recognized species remained unclear. We used genome sequencing to clarify relationships among these species and provided new names for most of these species. Six of the nine described Parastagonospora species were recovered from wheat, with five of these species coming from Iran. Genome sequences revealed that three strains thought to be hybrids between P. nodorum and P. pseudonodorum were not actually hybrids, but rather represented rare gene introgressions between those species. Our data are consistent with the hypothesis that P. nodorum originated as a pathogen of wild grasses in the Fertile Crescent, then emerged as a wheat pathogen via host-tracking during the domestication of wheat in the same region. The discovery of a diverse array of Parastagonospora species infecting wheat in Iran suggests that new wheat pathogens could emerge from this region in the future. Citation: Croll D, Crous PW, Pereira D, et al. 2021. Genome-scale phylogenies reveal relationships among Parastagonospora species infecting domesticated and wild grasses. Persoonia 46: 116-128. https://doi.org/10.3767/persoonia.2021.46.04.

Kupffer Cells Undergo Fundamental Changes during the Development of Experimental NASH and Are Critical in Initiating Liver Damage and Inflammation.

Non-alcoholic fatty liver disease has become the leading liver disease in North America and is associated with the progressive inflammatory liver disease non-alcoholic steatohepatitis (NASH). Considerable effort has been made to understand the role of resident and recruited macrophage populations in NASH however numerous questions remain. Our goal was to characterize the dynamic changes in liver macrophages during the initiation of NASH in a murine model. Using the methionine-choline deficient diet we found that liver-resident macrophages, Kupffer cells were lost early in disease onset followed by a robust infiltration of Ly-6C+ monocyte-derived macrophages that retained a dynamic phenotype. Genetic profiling revealed distinct patterns of inflammatory gene expression between macrophage subsets. Only early depletion of liver macrophages using liposomal clodronate prevented the development of NASH in mice suggesting that Kupffer cells are critical for the orchestration of inflammation during experimental NASH. Increased understanding of these dynamics may allow us to target potentially harmful populations whilst promoting anti-inflammatory or restorative populations to ultimately guide the development of effective treatment strategies.

QTL mapping of temperature sensitivity reveals candidate genes for thermal adaptation and growth morphology in the plant pathogenic fungus Zymoseptoria tritici.

Different thermal environments impose strong, differential selection on populations, leading to local adaptation, but the genetic basis of thermal adaptation is poorly understood. We used quantitative trait locus (QTL) mapping in the fungal wheat pathogen Zymoseptoria tritici to study the genetic architecture of thermal adaptation and identify candidate genes. Four wild-type strains originating from the same thermal environment were crossed to generate two mapping populations with 263 (cross 1) and 261 (cross 2) progeny. Restriction site-associated DNA sequencing was used to genotype 9745 (cross 1) and 7333 (cross 2) single-nucleotide polymorphism markers segregating within the mapping population. Temperature sensitivity was assessed using digital image analysis of colonies growing at two different temperatures. We identified four QTLs for temperature sensitivity, with unique QTLs found in each cross. One QTL had a logarithm of odds score >11 and contained only six candidate genes, including PBS2, encoding a mitogen-activated protein kinase kinase associated with low temperature tolerance in Saccharomyces cerevisiae. This and other QTLs showed evidence for pleiotropy among growth rate, melanization and growth morphology, suggesting that many traits can be correlated with thermal adaptation in fungi. Higher temperatures were highly correlated with a shift to filamentous growth among the progeny in both crosses. We show that thermal adaptation has a complex genetic architecture, with natural populations of Z. tritici harboring significant genetic variation for this trait. We conclude that Z. tritici populations have the potential to adapt rapidly to climate change and expand into new climatic zones.

Pyricularia graminis-tritici, a new Pyricularia species causing wheat blast.

Pyricularia oryzae is a species complex that causes blast disease on more than 50 species of poaceous plants. Pyricularia oryzae has a worldwide distribution as a rice pathogen and in the last 30 years emerged as an important wheat pathogen in southern Brazil. We conducted phylogenetic analyses using 10 housekeeping loci for 128 isolates of P. oryzae sampled from sympatric populations of wheat, rice, and grasses growing in or near wheat fields. Phylogenetic analyses grouped the isolates into three major clades. Clade 1 comprised isolates associated only with rice and corresponds to the previously described rice blast pathogen P. oryzae pathotype Oryza (PoO). Clade 2 comprised isolates associated almost exclusively with wheat and corresponds to the previously described wheat blast pathogen P. oryzae pathotype Triticum (PoT). Clade 3 contained isolates obtained from wheat as well as other Poaceae hosts. We found that Clade 3 is distinct from P. oryzae and represents a new species, Pyricularia graminis-tritici (Pgt). No morphological differences were observed among these species, but a distinctive pathogenicity spectrum was observed. Pgt and PoT were pathogenic and highly aggressive on Triticum aestivum (wheat), Hordeum vulgare (barley), Urochloa brizantha (signal grass), and Avena sativa (oats). PoO was highly virulent on the original rice host (Oryza sativa), and also on wheat, barley, and oats, but not on signal grass. We conclude that blast disease on wheat and its associated Poaceae hosts in Brazil is caused by multiple Pyricularia species. Pyricularia graminis-tritici was recently found causing wheat blast in Bangladesh. This indicates that P. graminis-tritici represents a serious threat to wheat cultivation globally.

First Report of Resistance to QoI Fungicides in North American Populations of Zymoseptoria tritici, Causal Agent of Septoria Tritici Blotch of Wheat.

The G143A mutation in cytb (cytochrome b gene) is associated with high levels of resistance to quinone outside inhibitor (QoI or strobilurin) fungicides that disrupt electron transport during cellular respiration (1). The G143A mutation in Zymoseptoria tritici (synonyms: Mycosphaerella graminicola and Septoria tritici), the causal agent of septoria tritici blotch of wheat (Triticum aestivum), was first reported in Europe in 2001 (1). Although Z. tritici has a global distribution (3), G143A mutants of Z. tritici have not been reported outside of Europe. We used PCR-RFLP (4) to estimate the frequencies of G143A mutants in Z. tritici populations at two locations in the Willamette Valley of western Oregon: the Hyslop Crop Science Field Research Laboratory (Hyslop Farm, HF), Benton County (44°37'52.85″ N, 123°11'55.19″ W) and research plots planted in a commercial wheat field in Washington County (45°33'58.53″ N, 123°00'11.78″ W) (North Valley Farm, NVF). Isolates originated from flag leaf collections from two cultivars ('Bobtail' and 'Tubbs 06') made in April and June of 2012 from plants in a replicated fungicide-treatment experiment, with isolates collected from both sprayed and unsprayed plots. Sixteen of the 169 isolates (9.5%) from HF possessed the G143A mutation (7 of 132 isolates from plots not receiving a QoI fungicide and 9 of 37 isolates collected from plots receiving two applications of the QoI azoxystrobin). One hundred forty six of the 175 isolates (83.4%) from NVF were G143A mutants (101 of 129 isolates from plots receiving no QoI fungicide and 45 of 46 isolates from plots receiving two applications of azoxystrobin). Results of phenotypic assays of a subset of 10 isolates from each location (5 mutants, 5 wild types from each location; 20 isolates altogether) supported a high level of resistance to azoxystrobin only in the G143A mutants. All 10 G143A mutants developed colonies after 8 days of growth on YMA plates amended with SHAM (2) and 1 ppm or 10 ppm azoxystrobin, with nine and eight G143A mutant isolates developing colonies on plates amended with 1 ppm and 10 ppm azoxystrobin, respectively. None of the wild-type isolates developed colonies on plates amended with SHAM and 1 ppm azoxystrobin, nor on plates amended with SHAM and 10 ppm azoxystrobin. All 20 isolates developed colonies on YMA plates lacking azoxystrobin, and treatments produced identical results across three replicates. These results are consistent with findings of higher levels of azoxystrobin resistance in G143A mutants compared to wild types in European populations (1). Isolates from HF and NVF differ in their previous exposure to QoI fungicides. The majority of the wheat area at HF is planted to breeding plots that are not sprayed with fungicide. Plots at NVF were planted in a commercial wheat field in a county where most wheat fields were treated with two to three applications of strobilurins each year over the past 4 years. Future monitoring for G143A mutants of Z. tritici throughout its range in North America will be necessary to assess whether strobilurin resistance will spread via wind-dispersal of ascospores or emerge de novo in treated fields. In Europe, stobilurins were first applied to wheat in 1996. G143A mutants of Z. tritici emerged de novo several times (4) and were widespread by 2007. References: (1) B. A. Fraaje et al. Phytopathology 95:933, 2005. (2) J. A. LaMondia. Tob. Sci. 49:1, 2012. (3) E. S. Orton et al. Mol. Plant Pathol. 12:413, 2011. (4) S. F. F. Torriani et al. Pest Manag. Sci. 65:155, 2008.

Divergence between sympatric rice- and maize-infecting populations of Rhizoctonia solani AG-1 IA from Latin America.

ABSTRACT The basidiomycetous fungus Rhizoctonia solani anastomosis group (AG)-1 IA is a major pathogen in Latin America causing sheath blight (SB) of rice. Particularly in Venezuela, the fungus also causes banded leaf and sheath blight (BLSB) on maize, which is considered an emerging disease problem where maize replaced traditional rice-cropping areas or is now planted in adjacent fields. Our goals in this study were to elucidate (i) the effects of host specialization on gene flow between sympatric and allopatric rice and maize-infecting fungal populations and (ii) the reproductive mode of the fungus, looking for evidence of recombination. In total, 375 isolates of R. solani AG1 IA sampled from three sympatric rice and maize fields in Venezuela (Portuguesa State) and two allopatric rice fields from Colombia (Meta State) and Panama (Chiriquí State) were genotyped using 10 microsatellite loci. Allopatric populations from Venezuela, Colombia, and Panama were significantly differentiated (Phi(ST) of 0.16 to 0.34). Partitioning of the genetic diversity indicated differentiation between sympatric populations from different host species, with 17% of the total genetic variation distributed between hosts while only 3 to 6% was distributed geographically among the sympatric Venezuelan fields. We detected symmetrical historical migration between the rice- and the maize-infecting populations from Venezuela. Rice- and maize-derived isolates were able to infect both rice and maize but were more aggressive on their original hosts, consistent with host specialization. Because the maize- and rice-infecting populations are still cross-pathogenic, we postulate that the genetic differentiation was relatively recent and mediated via a host shift. An isolation with migration analysis indicated that the maize-infecting population diverged from the rice-infecting population between 40 and 240 years ago. Our findings also suggest that maize-infecting populations have a mainly recombining reproductive system whereas the rice-infecting populations have a mixed reproductive system in Latin America.

Molecular evidence for recent founder populations and human-mediated migration in the barley scald pathogen Rhynchosporium secalis.

Rhynchosporium secalis is an important pathogen of barley globally. Fourteen polymorphic microsatellites were analyzed for 1664 R. secalis isolates sampled from 37 field populations to infer their demographic history. The results falsified the hypothesis that R. secalis co-evolved with its barley host in the Middle East. Populations from Scandinavia had significantly higher allelic diversities, the greatest number of private alleles and the highest genotypic diversities. All but three of the analyzed populations had an excess of gene diversity compared to the number of alleles, consistent with a recent population bottleneck. The remaining populations had a gene diversity deficit consistent with a population expansion following a recent population bottleneck in the last +/-100 years. A coalescent analysis revealed that the effective population sizes based on theta, of the analyzed populations were small relative to their ancestral population sizes, indicating that only a fraction of the diversity present in the ancestral populations was transmitted into current populations. These findings are consistent with the hypothesis that the pathogen population on barley experienced a selection bottleneck imposed by the host and/or are founder populations. The mean estimate of migration rates was 2.2 (avg 90% confidence interval=1.3-3.1). Major migration routes were identified among populations separated by long distances, eg between South Africa and Australia, as well as among North Africa, the Middle East and California, suggesting contemporary exchange of infected barley seed. In contrast with earlier findings, most populations exhibited significant gametic disequilibrium, probably as a result of genetic drift. We conclude that the majority of R. secalis populations have experienced human-mediated migration that led to numerous and relatively recent founder events around the world.

Genetic structure of populations of Rhizoctonia solani anastomosis group-1 IA from soybean in Brazil.

The Basidiomycete fungus Rhizoctonia solani anastomosis group (AG)-1 IA is a major pathogen of soybean in Brazil, where the average yield losses have reached 30 to 60% in some states in Northern Brazil. No information is currently available concerning levels of genetic diversity and population structure for this pathogen in Brazil. A total of 232 isolates of R. solani AG1 IA were collected from five soybean fields in the most important soybean production areas in central-western, northern, and northeastern Brazil. These isolates were genotyped using 10 microsatellite loci. Most of the multilocus genotypes (MLGTs) were site-specific, with few MLGTs shared among populations. Significant population subdivision was evident. High levels of admixture were observed for populations from Mato Grosso and Tocantins. After removing admixed genotypes, three out of five field populations (Maranhao, Mato Grosso, and Tocantins), were in Hardy-Weinberg (HW) equilibrium, consistent with sexual recombination. HW and gametic disequilibrium were found for the remaining soybean-infecting populations. The findings of low genotypic diversity, departures from HW equilibrium, gametic disequilibrium, and high degree of population subdivision in these R. solani AG-1 IA populations from Brazil are consistent with predominantly asexual reproduction, short-distance dispersal of vegetative propagules (mycelium or sclerotia), and limited long-distance dispersal, possibly via contaminated seed. None of the soybean-infecting populations showed a reduction in population size (bottleneck effect). We detected asymmetric historical migration among the soybean-infecting populations, which could explain the observed levels of subdivision.

Evolution of the CYP51 gene in Mycosphaerella graminicola: evidence for intragenic recombination and selective replacement.

Recent findings are consistent with a slow but constant shift towards reduced sensitivity of Mycosphaerella graminicola to azole fungicides, which target the CYP51 gene. The goal of this study was to elucidate the evolutionary mechanisms through which CYP51-based mutations associated with altered sensitivity have evolved in M. graminicola over space and time. To accomplish this, we sequenced and compared a portion of the CYP51 gene encompassing the main mutations associated with altered sensitivity towards demethylation inhibitor fungicides. The CYP51 gene showed an extraordinary dynamic shift consistent with a selective haplotype replacement both in space and in time. No mutations associated with increased resistance to azoles were found in non-European populations. These mutations were also absent in the oldest collections from Europe, whereas they dominated in the recent European populations. Intragenic recombination was identified as an important evolutionary process in populations affected by high fungicide selection, suggesting the creation of novel alleles among existing mutations as a potential source of novel resistance alleles. We propose that CYP51 mutations giving resistance in M. graminicola arose only locally (perhaps in Denmark or the UK) and were then spread eastward across Europe through wind-dispersed ascospores. We conclude that recurring cycles of recombination coupled with selection due to the widespread use of azole fungicides will increase the frequency of novel mutants or recombinants with higher resistance. Long-distance gene flow due to wind dispersal of ascospores will move the resulting new alleles to new areas following the prevailing wind directions. A selective replacement favouring haplotypes with various coding mutations at the target site for azole fungicides during the last 5-10 years is the most likely cause of the decrease in sensitivity reported for many azole fungicides in the same period.

Highly polymorphic microsatellite loci in the rice- and maize-infecting fungal pathogen Rhizoctonia solani anastomosis group 1 IA.

Ten polymorphic microsatellite loci were isolated and characterized from the rice- and maize-infecting Basidiomycete fungus Rhizoctonia solani anastomosis group AG-1 IA. All loci were polymorphic in two populations from Louisiana in USA and Venezuela. The total number of alleles per locus ranged from four to eight. All 10 loci were also useful for genotyping soybean-infecting R. solani AG-1 isolates from Brazil and USA. One locus, TC06, amplified across two other AG groups representing different species, showing species-specific repeat length polymorphism. This marker suite will be used to determine the global population structure of this important pathogenic fungus.

The origin and colonization history of the barley scald pathogen Rhynchosporium secalis.

The origins of pathogens and their past and present migration patterns are often unknown. We used phylogenetic haplotype clustering in conjunction with model-based coalescent approaches to reconstruct the genetic history of the barley leaf pathogen Rhynchosporium secalis using the avirulence gene NIP1 and its flanking regions. Our results falsify the hypothesis that R. secalis emerged in association with its host during the domestication of barley 10,000 to 15,000 years ago in the Fertile Crescent and was introduced into Europe through the migration of Neolithic farmers. Estimates of time since most recent common ancestor (2500-5000 BP) placed the emergence of R. secalis clearly after the domestication of barley. We propose that modern populations of R. secalis originated in northern Europe following a host switch, most probably from a wild grass onto cultivated barley shortly after barley was introduced into northern Europe. R. secalis subsequently spread southwards into already established European barley-growing areas.

Sexual reproduction facilitates the adaptation of parasites to antagonistic host environments: Evidence from empirical study in the wheat-Mycosphaerella graminicola system.

Most eukaryotes use sexual reproduction to transmit genetic information from generation to generation despite the advantages offered by asexual reproduction. One theory to explain the origin and maintenance of sexual reproduction hypothesises that sexual recombination generates genetic variation that allows faster adaptation to fluctuating and/or stressful environments. We used a combination of ecological, molecular genetic, statistical and experimental evolution approaches to test this hypothesis in an agricultural plant-pathogen system. We inoculated wheat hosts with 10 strains of the fungal pathogen Mycosphaerella graminicola in a field experiment and estimated the contributions of sexual reproduction, asexual reproduction and immigration to the genetic composition of fungal populations sampled from moderately resistant and susceptible hosts through the course of an epidemic cycle. We found that a significant proportion of the M. graminicola population in the late phase of the epidemic originated from sexual reproduction among isolates that had been introduced into the field plots at the beginning of the epidemic. Recombinants were recovered at a higher frequency on the moderately resistant plant host Madsen than on the susceptible host Stephens. By the end of the growing season, we estimated that approximately 13% of the strains sampled from the resistant host were recombinants, compared with 9% in the samples collected from the susceptible host. We also found that pathogen strains originating from the resistant cultivar displayed higher levels of fitness, virulence and fungicide tolerance than those originating from the susceptible cultivar. Our results provide empirical support for the hypothesis that sexual reproduction facilitates the evolution of parasites to overcome host resistance.

Selection for increased cyproconazole tolerance in Mycosphaerella graminicola through local adaptation and in response to host resistance.

SUMMARY Sterol demethylation inhibitors (DMIs) represent one of the largest groups of systemic fungicides that have been used to control agriculturally important fungal pathogens. Knowledge regarding the evolution of fungicide resistance in agricultural ecosystems is fragmentary and a better understanding of the processes driving the development of DMI resistance in populations of fungal pathogens is needed by plant pathologists and the agrochemical industry. We considered some of these processes using approaches based on molecular population and quantitative genetics. Five Mycosphaerella graminicola populations sampled from unsprayed wheat fields on four continents were assayed for eight restriction fragment length polymorphism (RFLP) markers and their level of tolerance to cyproconazole. DMI fungicides such as cyproconazole inhibit the enzyme eburicol 14-alpha-demethylase. The gene encoding this target, CYP51, was sequenced for all isolates. We found unimodal, continuous variations in cyproconazole tolerance among the M. graminicola isolates sampled from individual fields, consistent with a polygenic mode of inheritance. We also found that population differentiation for cyproconazole tolerance (Q(ST)) among the five M. graminicola populations was significantly higher than the corresponding population differentiation for neutral RFLP markers (G(ST)), suggesting that selection for cyproconazole tolerance in the Swiss population has already led to local adaptation that can be seen even in an unsprayed population. The Swiss population displayed the highest level of tolerance to cyproconazole, in addition to a lower than expected quantitative variation in fungicide tolerance and a skewed distribution, indicating that selection had increased the overall tolerance of this population. Further analysis with DNA sequencing showed that the population from Switzerland was dominated by isolates with several point mutations and a 6-bp deletion in CYP51. This deletion and one of the point mutations were previously related to increased resistance in field isolates. The fungal population from Oregon sampled from an unsprayed resistant host cultivar displayed the same gene diversity in RFLP loci but higher cyproconazole tolerance and quantitative variation in tolerance than the fungal population from the same field sampled from an unsprayed susceptible host cultivar.

Variation for neutral markers is correlated with variation for quantitative traits in the plant pathogenic fungus Mycosphaerella graminicola.

We compared genetic variation and population differentiation at RFLP marker loci with seven quantitative characters including fungicide resistance, temperature sensitivity, pycnidial size, pycnidial density, colony size, percentage of leaves covered by pycnidia (PLACP) and percentage of leaves covered by lesions (PLACL) in Mycosphaerella graminicola populations sampled from four regions. Wide variation in population differentiation was found across the quantitative traits assayed. Fungicide resistance, temperature sensitivity, and PLACP displayed a significantly higher Q(ST) than G(ST), consistent with selection for local adaptation, while pycnidial size, pycnidial density and colony size displayed a lower or significantly lower Q(ST) than G(ST), consistent with constraining selection. There was not a statistical difference between Q(ST) and G(ST) in PLACL. We also found a positive and significant correlation between genetic variation in molecular marker loci and quantitative traits at the multitrait scale, suggesting that estimates of overall genetic variation for quantitative traits in M. graminicola could be derived from analysis of the molecular genetic markers.

Migration patterns among global populations of the pathogenic fungus Mycosphaerella graminicola.

DNA sequences from five nuclear loci and data from three microsatellites were collected from 360 isolates representing 14 globally distributed populations of the plant pathogenic fungus Mycosphaerella graminicola. Haplotype networks were constructed for the five sequence loci and population subdivision was assessed using Hudson's permutation test. Migration estimates were calculated using six regional populations for both the sequence and microsatellite loci. While subdivision was detected among the six regional populations, significant gene flow was indicated among some of the populations. The European and Israeli populations contributed the majority of historical immigrants to the New World. Migration estimates for microsatellite loci were used to infer more recent migration events among specific New World populations. We conclude that gene flow was an important factor in determining the demographic history of Mycosphaerella graminicola.

The interaction among evolutionary forces in the pathogenic fungus Mycosphaerella graminicola.

The population genetic dynamic of a species is driven by interactions among mutation, migration, drift, mating system, and selection, but it is rare to have sufficient empirical data to estimate values for all of these forces and to allow comparison of the relative magnitudes of these evolutionary forces. We combined data from a mark-release-recapture experiment, extensive population surveys, and computer simulations to evaluate interactions among these evolutionary forces in the pathogenic fungus Mycosphaerella graminicola. The results from these studies showed that, on average, the immigration rate was 0.027, the fraction of outcrossing individuals was 0.035, and the selection coefficient associated with immigrants was 0.106 each generation. We also estimated that effective population sizes for this fungus were larger than 24,000 and the mutation rate for the RFLP markers used in surveys and field experiments was approximately 4 x 10(-5). Computer simulations based on these estimates indicate that, on average, the global population of M. graminicola has reached equilibrium. Population genetic parameters including number of alleles, gene diversity, and population subdivision estimated from the computer simulations were surprisingly close to empirical estimates. Simulations also revealed that random drift is the major evolutionary force decreasing genetic variation in this fungus, followed by natural selection. The major force adding to genetic variation was mutation, followed by gene flow and sexual recombination. Gene flow played the leading role in decreasing population subdivision while natural selection was the major factor increasing population subdivision.

Phylogenetic analysis of globally distributed Mycosphaerella graminicola populations based on three DNA sequence loci.

DNA sequence data from three nuclear loci were collected from 384 isolates representing fourteen globally distributed populations of the plant pathogenic fungus Mycosphaerella graminicola. Gene genealogies were constructed for the actin and beta-tubulin loci as well as for the previously characterized RFLP locus STS2. The STS2 and beta-tubulin loci showed greater potential for phylogenetic studies than the actin locus. Greater sequence diversity was found in the "Old World" populations (Middle East and Europe) than in the "New World" populations (North and South America and Australia). The gene trees were rooted using homologous DNA sequences of Septoria passerinii, the closest known relative to M. graminicola, as well as coalescent rooting. Based on the rooted trees, a tentative phylogenetic history of these populations was inferred. The Middle East appears to be the most likely center of origin, while European populations are more ancient than New World populations. A test for neutrality indicated that the intron in the actin locus could be under selection, while the other two sequence loci were neutral.

Evidence for Natural Selection in the Mitochondrial Genome of Mycosphaerella graminicola.

ABSTRACT Pathogenicity assays were combined with restriction fragment length polymorphism (RFLP) markers in the mitochondrial and nuclear genomes to compare Mycosphaerella graminicola populations adapted to bread wheat (Triticum aestivum) and durum wheat (T. turgidum) in the Mediterranean Basin. The majority of isolates had unique nuclear DNA fingerprints and multilocus haplotypes. Only six mitochondrial DNA (mtDNA) haplotypes were identified among 108 isolates assayed. There were minor differences in frequencies of alleles at nuclear RFLP loci between the two host-adapted populations, but differences in the frequencies of mtDNA haplotypes were highly significant (P < 0.0001). mtDNA haplotype 1 dominated on the isolates adapted to bread wheat, and its frequency was twice as high as for the isolates adapted to durum wheat. mtDNA haplotype 4, which contained a unique approximately 3-kb insertion, was detected only in isolates showing specificity toward durum wheat and was the dominant haplotype on this species. We propose that the low mitochondrial diversity in this pathogenic fungus is due to a selective sweep and that differences in the frequencies of mtDNA haplotypes between the two host-adapted populations were due to natural selection according to host species.

The global genetic structure of the wheat pathogen Mycosphaerella graminicola is characterized by high nuclear diversity, low mitochondrial diversity, regular recombination, and gene flow.

A total of 1673 Mycosphaerella graminicola strains were assayed for DNA fingerprints and restriction fragment length polymorphism (RFLP) markers in the nuclear and mitochondrial genomes. The isolates were collected from 17 wheat fields located in 11 countries on five continents over a six year period (1989-1995). Our results indicate that genetic diversity in the nuclear genome of this fungus was high for all but three of the field populations surveyed and that populations sampled from different continents had similar frequencies for the most common RFLP alleles. Hierarchical analysis revealed that more than 90% of global gene diversity was distributed within a wheat field, while approximately 5% of gene diversity was distributed among fields within regions and approximately 3% was distributed among regions on different continents. These findings suggest that gene flow has occurred on a global scale. On average, each leaf was colonized by a different nuclear genotype. In contrast, only seven mtDNA haplotypes were detected among the 1673 isolates and the two most common mtDNA haplotypes represented approximately 93% of the world population, consistent with a selective sweep. Analysis of multilocus associations indicated that all field populations were in gametic equilibrium, suggesting that sexual recombination is a regular occurrence globally.