Evidence for Natural Selection in the Mitochondrial Genome of Mycosphaerella graminicola

Multiple rearrangements and low inter- and intra-species mitogenome sequence variation in the Heterobasidion annosum s.l. species complex

Introduction

Mitochondria are essential organelles in the eukaryotic cells and responsible for the energy production but are also involved in many other functions including virulence of some fungal species. Although the evolution of fungal mitogenomes have been studied at some taxonomic levels there are still many things to be learned from studies of closely related species.

Methods

In this study, we have analyzed 60 mitogenomes in the five species of the Heterobasidion annosum sensu lato complex that all are necrotrophic pathogens on conifers.

Results and Discussion

Compared to other fungal genera the genomic and genetic variation between and within species in the complex was low except for multiple rearrangements. Several translocations of large blocks with core genes have occurred between the five species and rearrangements were frequent in intergenic areas. Mitogenome lengths ranged between 108 878 to 116 176 bp, mostly as a result of intron variation. There was a high degree of homology of introns, homing endonuclease genes, and intergenic ORFs among the five Heterobasidion species. Three intergenic ORFs with unknown function (uORF6, uORF8 and uORF9) were found in all five species and was located in conserved synteny blocks. A 13 bp long GC-containing self-complementary palindrome was discovered in many places in the five species that were optional in presence/absence. The within species variation is very low, among 48 H. parviporum mitogenomes, there was only one single intron exchange, and SNP frequency was 0.28% and indel frequency 0.043%. The overall low variation in the Heterobasidion annosum sensu lato complex suggests a slow evolution of the mitogenome.

Mitochondrial Genome Contributes to the Thermal Adaptation of the Oomycete Phytophthora infestans

As a vital element of climate change, elevated temperatures resulting from global warming present new challenges to natural and agricultural sustainability, such as ecological disease management. Mitochondria regulate the energy production of cells in responding to environmental fluctuation, but studying their contribution to the thermal adaptation of species is limited. This knowledge is needed to predict future disease epidemiology for ecology conservation and food security. Spatial distributions of the mitochondrial genome (mtDNA) in 405 Phytophthora infestans isolates originating from 15 locations were characterized. The contribution of MtDNA to thermal adaptation was evaluated by comparative analysis of mtDNA frequency and intrinsic growth rate, relative population differentiation in nuclear and mtDNA, and associations of mtDNA distribution with local geography climate conditions. Significant variation in frequency, intrinsic growth rate, and spatial distribution was detected in mtDNA. Population differentiation in mtDNA was significantly higher than that in the nuclear genome, and spatial distribution of mtDNA was strongly associated with local climatic conditions and geographic parameters, particularly air temperature, suggesting natural selection caused by a local temperature is the main driver of the adaptation. Dominant mtDNA grew faster than the less frequent mtDNA. Our results provide useful insights into the evolution of pathogens under global warming. Given its important role in biological functions and adaptation to local air temperature, mtDNA intervention has become an increasing necessity for future disease management. To secure ecological integrity and food production under global warming, a synergistic study on the interactive effect of changing temperature on various components of biological and ecological functions of mitochondria in an evolutionary frame is urgently needed.

Mitonuclear Genetic Interactions in the Basidiomycete Heterobasidion parviporum Involve a Non-conserved Mitochondrial Open Reading Frame

The mitochondrial and nuclear genomes of Eukaryotes are inherited separately and consequently follow distinct evolutionary paths. Nevertheless, the encoding of many mitochondrial proteins by the nuclear genome shows the high level of integration they have reached, which makes mitonuclear genetic interactions all the more conceivable. For each species, natural selection has fostered the evolution of coadapted alleles in both genomes, but a population-wise divergence of such alleles could lead to important phenotypic variation, and, ultimately, to speciation. In this study in the Basidiomycete Heterobasidion parviporum, we have investigated the genetic basis of phenotypic variation among laboratory-designed heterokaryons carrying the same pair of haploid nuclei, but a different mitochondrial genome. Radial growth rate data of thirteen unrelated homokaryotic parents and of their heterokaryotic offspring were combined with SNP data extracted from parental genome sequences to identify nuclear and mitochondrial loci involved in mitonuclear interactions. Two nuclear loci encoding mitochondrial proteins appeared as best candidates to engage in a genetic interaction affecting radial growth rate with a non-conserved mitochondrial open reading frame of unknown function and not reported apart from the Russulales order of Basidiomycete fungi. We believe our approach could be useful to investigate several important traits of fungal biology where mitonuclear interactions play a role, including virulence of fungal pathogens.

Fungal Mitogenomes: Relevant Features to Planning Plant Disease Management

Mitochondrial genomes (mt-genomes) are characterized by a distinct codon usage and their autonomous replication. Mt-genomes encode highly conserved genes (mt-genes), like proteins involved in electron transport and oxidative phosphorylation but they also carry highly variable regions that are in part responsible for their high plasticity. The degree of conservation of their genes is such that they allow the establishment of phylogenetic relationships even across distantly related species. Here, we describe the mechanisms that generate changes along mt-genomes, which play key roles at enlarging the ability of fungi to adapt to changing environments. Within mt-genomes of fungal pathogens, there are dispensable as well as indispensable genes for survival, virulence and/or pathogenicity. We also describe the different complexes or mechanisms targeted by fungicides, thus addressing a relevant issue regarding disease management. Despite the controversial origin and evolution of fungal mt-genomes, the intrinsic mechanisms and molecular biology involved in their evolution will help to understand, at the molecular level, the strategies for fungal disease management.

Overview of genomic and bioinformatic resources for Zymoseptoria tritici

Zymoseptoria tritici (syn. Mycosphaerella graminicola, Septoria tritici) is a haploid fungus belonging to the class Dothideomycetes. It is the causal agent of septoria leaf blotch - one of the world's most significant diseases of wheat. Here we review the genomic and bioinformatic resources that have been generated for Z. tritici. These include the whole-genome reference assembly for isolate IPO323, genome resequencing of alternate isolates, mitochondrial genome sequences, transcriptome sequences and expression data, and annotations of gene structure and function. We also highlight important advances in our fundamental knowledge of genome evolution and its effects on adaptation and pathogenicity in Z. tritici that have been facilitated by these resources.

Genetics of resistance to Zymoseptoria tritici and applications to wheat breeding

This paper reviews current knowledge about genes for resistance to Septoria tritici blotch (STB) of wheat, caused by Zymoseptoria tritici (formerly Mycosphaerella graminicola). These genes can be placed into two classes, although a few may have characteristics of both classes. Qualitative resistance is controlled by genes which control large fractions of genetic variation, 21 of which have been discovered and mapped so far. Most of them have been shown to be genotype-specific, being effective against the minority of Z. tritici isolates which are avirulent, and Stb6 has been shown to control a gene-for-gene relationship. Most qualitative resistances are unlikely to be durable and some formerly effective genes have been overcome by the evolution of pathogen virulence. Quantitative resistance is generally controlled by genes with small-to-moderate effects on STB. They have generally weaker specificity than qualitative genes and have provided more durable resistance. 89 genome regions carrying quantitative trait loci (QTL) or meta-QTL have been identified to date. Some QTL have been mapped at or near loci of qualitative genes, especially Stb6, which is present in several sources of resistance. Another gene of particular interest is Stb16q, which has been effective against all Z. tritici isolates tested so far. In addition to resistance, the susceptibility of wheat cultivars to STB can also be reduced by disease escape traits, some of which may be undesirable in breeding. The fundamental requirements for breeding for STB-resistance are genetic diversity for resistance in wheat germplasm and a field trial site at which STB epidemics occur regularly and effective selection can be conducted for resistance combined with other desirable traits. If these are in place, knowledge of resistance genes can be applied to improving control of STB.

QoI Resistance and Mitochondrial Genetic Structure of Zymoseptoria tritici in Morocco

In total, 230 single-conidial isolates of the fungal wheat pathogen Zymoseptoria tritici (formerly Septoria tritici, teleomorph: Mycosphaerella graminicola) were sampled in Morocco in 2008 and 2010 to assess resistance against quinone outside inhibitors (QoIs), a widely used group of fungicides in wheat pest management. All 134 isolates sampled in 2008 were QoI sensitive. In contrast, 9 of the 96 isolates from the 2010 collection were resistant, suggesting a recent emergence of the resistance. Mitochondrial (mt)DNA-sequence analyses identified four haplotypes among the resistant isolates. Wright's F statistics (F_ST) analyses from mtDNA sequences revealed a shallow population structure of Z. tritici within Morocco and a substantial asymmetric gene flow from Europe into Morocco. A phylogenetic reconstruction including Moroccan and European isolates clustered the haplotypes regardless of their geographic origin. The four Moroccan QoI-resistant mitochondrial haplotypes clustered in two distinct clades in the tree topology, suggesting at least two independent origins of the resistance. This study reported, for the first time, the occurrence of QoI-resistant genotypes of Z. tritici in Morocco. Our findings are consistent with the hypothesis that QoI resistance emerged very recently through parallel genetic adaptation in Morocco, although gene flow from Europe cannot be excluded.

Population biology of fungal plant pathogens

Studies of the population genetics of fungal and oomycetous phytopathogens are essential to clarifying the disease epidemiology and devising management strategies. Factors commonly associated with higher organisms such as migration, natural selection, or recombination, are critical for the building of a clearer picture of the pathogen in the landscape. In this chapter, we focus on a limited number of experimental and analytical methods that are commonly applied in population genetics. At first, we present different types of qualitative and quantitative traits that could be identified morphologically (phenotype). Subsequently, we describe several molecular methods based on dominant and codominant markers, and we provide our assessment of the advantages and shortfalls of these methods. Third, we discuss various analytical methods, which include phylogenies, summary statistics as well as coalescent-based methods, and we elaborate on the benefits associated with each approach. Last, we develop a case study in which we investigate the population structure of the fungal phytopathogen Verticillium dahliae in coastal California, and assess the hypotheses of transcontinental gene flow and recombination in a fungus that is described as asexual.

Evolutionary history of the mitochondrial genome in Mycosphaerella populations infecting bread wheat, durum wheat and wild grasses

Plant pathogens emerge in agro-ecosystems following different evolutionary mechanisms over different time scales. Previous analyses based on sequence variation at six nuclear loci indicated that Mycosphaerella graminicola diverged from an ancestral population adapted to wild grasses during the process of wheat domestication approximately 10,500 years ago. We tested this hypothesis by conducting coalescence analyses based on four mitochondrial loci using 143 isolates that included four closely related pathogen species originating from four continents. Pathogen isolates from bread and durum wheat were included to evaluate the emergence of specificity towards these hosts in M. graminicola. Although mitochondrial and nuclear genomes differed greatly in degree of genetic variability, their coalescence was remarkably congruent, supporting the proposed origin of M. graminicola through host tracking. The coalescence analysis was unable to trace M. graminicola host specificity through recent evolutionary time, indicating that the specificity towards durum or bread wheat emerged following the domestication of the pathogen on wheat.

Population genetics of fungal and oomycete effectors involved in gene-for-gene interactions

Antagonistic coevolution between plants and pathogens has generated a broad array of attack and defense mechanisms. In the classical avirulence (Avr) gene-for-gene model, the pathogen gene evolves to escape host recognition while the host resistance (R) gene evolves to track the evolving pathogen elicitor. In the case of host-specific toxins (HST), the evolutionary arms race may be inverted, with the gene encoding the pathogen toxin evolving to maintain recognition of the host sensitivity target while the host sensitivity gene evolves to escape binding with the toxin. Pathogen effector genes, including those encoding Avr elicitors and HST, often show elevated levels of polymorphism reflecting the coevolutionary arms race between host and pathogen. However, selection can also eliminate variation in the coevolved gene and its neighboring regions when advantageous alleles are swept to fixation. The distribution and diversity of corresponding host genes will have a major impact on the distribution and diversity of effectors in the pathogen population. Population genetic analyses including both hosts and their pathogens provide an essential tool to understand the diversity and dynamics of effector genes. Here, we summarize current knowledge about the population genetics of fungal and oomycete effector genes, focusing on recent studies that have used both spatial and temporal collections to assess the diversity and distribution of alleles and to monitor changes in allele frequencies over time. These studies illustrate that effector genes exhibit a significant degree of diversity at both small and large sampling scales, suggesting that local selection plays an important role in their evolution. They also illustrate that Avr elicitors and HST may be recognizing the same R genes in plants, leading to evolutionary outcomes that differ for necrotrophs and biotrophs while affecting the evolution of the corresponding R genes. Under this scenario, the optimal number of R genes in a plant genome may be determined by the relative abundance of necrotrophic and biotrophic pathogens in the plant's environment.

Concordant evolution of mitochondrial and nuclear genomes in the wheat pathogen Phaeosphaeria nodorum

We compared patterns of mitochondrial restriction fragment length polymorphism (RFLP) diversity with patterns of nuclear RFLP diversity to investigate the effects of selection, gene flow, and sexual reproduction on the population genetic structure and evolutionary history of the wheat pathogen Phaeosphaeria nodorum. A total of 315 fungal isolates from Texas, Oregon, and Switzerland were analyzed using seven nuclear RFLP probes that hybridized to discrete loci and purified mitochondrial DNA that hybridized to the entire mtDNA genome. Forty-two different mitochondrial haplotypes and 298 different nuclear haplotypes were detected. The two most frequent mtDNA haplotypes were present in every population and represented 32% of all isolates. High levels of gene flow, low levels of population subdivision, no evidence for either host specificity or cyto-nuclear disequilibrium were inferred from the analysis of both genomes. The concordance in estimates of these population genetic parameters from both genomes suggests that the two genomes experienced similar degrees of migration, genetic drift and selection.

Significant difference in pathogenicity between MAT1-1 and MAT1-2 isolates in the wheat pathogen Mycosphaerella graminicola

Five Mycosphaerella graminicola populations from four geographic regions (Australia, Israel, Switzerland, and the USA) were assayed for neutral RFLP markers and mating type idiomorphs. On average, 25-30 genetically distinct isolates were selected from each population and their pathogenicity was measured on two wheat cultivars in a common garden experiment conducted in a greenhouse. A significant difference in pathogenicity was found between MAT1-1 and MAT1-2 isolates. On average, MAT1-1 isolates had 14-22% greater pathogenicity than MAT1-2 isolates. The pattern of higher pathogenicity in MAT1-1 isolates was consistent across four geographical populations and on two wheat cultivars. A uniform and continuous variation in pathogenicity was found among isolates within each mating type, but no genetic differentiation in selectively neutral RFLP loci was found between mating types, consistent with the hypothesis that differences in pathogenicity were not due to the effects of specific pathogenicity genes or non-random genetic backgrounds.

Differential Selection on Rhynchosporium secalis During Parasitic and Saprophytic Phases in the Barley Scald Disease Cycle

ABSTRACT Competition among eight Rhynchosporium secalis isolates was assessed during parasitic and saprophytic phases of the disease cycle in field experiments conducted at two locations and over two growing seasons. The eight isolates were inoculated onto six barley populations exhibiting varying degrees of resistance. Microsatellite analysis of 2,866 isolates recovered from the field experiments showed significant, and sometimes opposite, changes in the frequencies of R. secalis genotypes during the growing season (parasitic phase) and between growing seasons (saprophytic phase). Isolates that showed the most complex virulence in greenhouse seedling assays had the lowest fitness in the field experiment. Significant differences in isolate fitness were found on different host populations and in different environments. Selection coefficients were large, indicating that evolution can occur rapidly in field populations. Although inoculated isolates had the lowest overall fitness on the moderately resistant landrace cv. Arabi Aswad, some isolates were more virulent and consistently increased in frequency on this landrace, suggesting a risk of directional selection and possible erosion of the resistance following its widespread deployment in monoculture. These results provide the first direct evidence that R. secalis pathogen genotypes differ in their saprophytic ability and parasitic fitness under field conditions.

MgSlt2, a cellular integrity MAP kinase gene of the fungal wheat pathogen Mycosphaerella graminicola, is dispensable for penetration but essential for invasive growth

Among expressed sequence tag libraries of Mycosphaerella graminicola isolate IPO323, we identified a full-length cDNA clone with high homology to the mitogen-activated protein (MAP) kinase Slt2 in Saccharomyces cerevisiae. This MAP kinase consists of a 1242-bp open reading frame, and encodes a 414-amino-acid protein. We designated this homolog MgSlt2, generated MgSlt2 knockout strains in M. graminicola isolate IPO323, and found several altered phenotypes in vitro as well as in planta. In yeast glucose broth, MgSlt2 disruptants showed a defective polarized growth in the tip cells upon aging, causing substantial local enlargements culminating in large swollen cells containing two to four nuclei. The MgSlt2 disruptants showed a significantly increased sensitivity to several fungicides, including miconazole (2x), bifonazole (>4x), imazalil (5x), and cyproconazole (10x), and were hypersensitive to glucanase. Unlike the wild type, MgSlt2 disruptants did not produce aerial mycelia and did not melanize on potato dextrose agar. Although cytological analysis in planta showed normal penetration of wheat stomata by the germ tubes of the MgSlt2 disruptants, subsequently formed hyphal filaments frequently were unable to branch out and establish invasive growth resulting in highly reduced virulence, and prevented pycnidia formation. Therefore, we conclude that MgSlt2 is a new pathogenicity factor in M. graminicola.

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.