Stefano F.F. Torriani

Finished Genome of the Fungal Wheat Pathogen Mycosphaerella graminicola Reveals Dispensome Structure, Chromosome Plasticity, and Stealth Pathogenesis

The plant-pathogenic fungus Mycosphaerella graminicola (asexual stage: Septoria tritici) causes septoria tritici blotch, a disease that greatly reduces the yield and quality of wheat. This disease is economically important in most wheat-growing areas worldwide and threatens global food production. Control of the disease has been hampered by a limited understanding of the genetic and biochemical bases of pathogenicity, including mechanisms of infection and of resistance in the host. Unlike most other plant pathogens, M. graminicola has a long latent period during which it evades host defenses. Although this type of stealth pathogenicity occurs commonly in Mycosphaerella and other Dothideomycetes, the largest class of plant-pathogenic fungi, its genetic basis is not known. To address this problem, the genome of M. graminicola was sequenced completely. The finished genome contains 21 chromosomes, eight of which could be lost with no visible effect on the fungus and thus are dispensable. This eight-chromosome dispensome is dynamic in field and progeny isolates, is different from the core genome in gene and repeat content, and appears to have originated by ancient horizontal transfer from an unknown donor. Synteny plots of the M. graminicola chromosomes versus those of the only other sequenced Dothideomycete, Stagonospora nodorum, revealed conservation of gene content but not order or orientation, suggesting a high rate of intra-chromosomal rearrangement in one or both species. This observed “mesosynteny” is very different from synteny seen between other organisms. A surprising feature of the M. graminicola genome compared to other sequenced plant pathogens was that it contained very few genes for enzymes that break down plant cell walls, which was more similar to endophytes than to pathogens. The stealth pathogenesis of M. graminicola probably involves degradation of proteins rather than carbohydrates to evade host defenses during the biotrophic stage of infection and may have evolved from endophytic ancestors.

Dothideomycete–Plant Interactions Illuminated by Genome Sequencing and EST Analysis of the Wheat Pathogen Stagonospora nodorum

Stagonospora nodorum is a major necrotrophic fungal pathogen of wheat (Triticum aestivum) and a member of the Dothideomycetes, a large fungal taxon that includes many important plant pathogens affecting all major crop plant families. Here, we report the acquisition and initial analysis of a draft genome sequence for this fungus. The assembly comprises 37,164,227 bp of nuclear DNA contained in 107 scaffolds. The circular mitochondrial genome comprises 49,761 bp encoding 46 genes, including four that are intron encoded. The nuclear genome assembly contains 26 classes of repetitive DNA, comprising 4.5% of the genome. Some of the repeats show evidence of repeat-induced point mutations consistent with a frequent sexual cycle. ESTs and gene prediction models support a minimum of 10,762 nuclear genes. Extensive orthology was found between the polyketide synthase family in S. nodorum and Cochliobolus heterostrophus, suggesting an ancient origin and conserved functions for these genes. A striking feature of the gene catalog was the large number of genes predicted to encode secreted proteins; the majority has no meaningful similarity to any other known genes. It is likely that genes for host-specific toxins, in addition to ToxA, will be found among this group. ESTs obtained from axenic mycelium grown on oleate (chosen to mimic early infection) and late-stage lesions sporulating on wheat leaves were obtained. Statistical analysis shows that transcripts encoding proteins involved in protein synthesis and in the production of extracellular proteases, cellulases, and xylanases predominate in the infection library. This suggests that the fungus is dependant on the degradation of wheat macromolecular constituents to provide the carbon skeletons and energy for the synthesis of proteins and other components destined for the developing pycnidiospores.

QoI resistance emerged independently at least 4 times in European populations of Mycosphaerella graminicola

BACKGROUND QoI fungicides or quinone outside inhibitors (also called strobilurins) have been widely used to control agriculturally important fungal pathogens since their introduction in 1996. Strobilurins block the respiration pathway by inhibiting the cytochrome bc1 complex in mitochondria. Several plant pathogenic fungi have developed field resistance. The first QoI resistance in Mycosphaerella graminicola (Fuckel) Schroter was detected retrospectively in UK in 2001 at a low frequency in QoI-treated plots. During the following seasons, resistance reached high frequencies across northern Europe. The aim of this study was to identify the main evolutionary forces driving the rapid emergence and spread of QoI resistance in M. graminicola populations. RESULTS The G143A mutation causing QoI resistance was first detected during 2002 in all tested populations and in eight distinct mtDNA sequence haplotypes. By 2004, 24 different mtDNA haplotypes contained the G143A mutation. Phylogenetic analysis showed that strobilurin resistance was acquired independently through at least four recurrent mutations at the same site of cytochrome b. Estimates of directional migration rates showed that the majority of gene flow in Europe had occurred in a west-to-east direction. CONCLUSION This study demonstrated that recurring mutations independently introduced the QoI resistance allele into different genetic and geographic backgrounds. The resistant haplotypes then increased in frequency owing to the strong fungicide selection and spread eastward through wind dispersal of ascospores.

Coevolution and Life Cycle Specialization of Plant Cell Wall Degrading Enzymes in a Hemibiotrophic Pathogen

Zymoseptoria tritici is an important fungal pathogen on wheat that originated in the Fertile Crescent. Its closely related sister species Z. pseudotritici and Z. ardabiliae infect wild grasses in the same region. This recently emerged host–pathogen system provides a rare opportunity to investigate the evolutionary processes shaping the genome of an emerging pathogen. Here, we investigate genetic signatures in plant cell wall degrading enzymes (PCWDEs) that are likely affected by or driving coevolution in plant-pathogen systems. We hypothesize four main evolutionary scenarios and combine comparative genomics, transcriptomics, and selection analyses to assign the majority of PCWDEs in Z. tritici to one of these scenarios. We found widespread differential transcription among different members of the same gene family, challenging the idea of functional redundancy and suggesting instead that specialized enzymatic activity occurs during different stages of the pathogen life cycle. We also find that natural selection has significantly affected at least 19 of the 48 identified PCWDEs. The majority of genes showed signatures of purifying selection, typical for the scenario of conserved substrate optimization. However, six genes showed diversifying selection that could be attributed to either host adaptation or host evasion. This study provides a powerful framework to better understand the roles played by different members of multigene families and to determine which genes are the most appropriate targets for wet laboratory experimentation, for example, to elucidate enzymatic function during relevant phases of a pathogen’s life cycle.

Evidence for Extensive Recent Intron Transposition in Closely Related Fungi

Though spliceosomal introns are a major structural component of most eukaryotic genes and intron density varies by more than three orders of magnitude among eukaryotes [1-3], the origins of introns are poorly understood, and only a few cases of unambiguous intron gain are known [4-8]. We utilized population genomic comparisons of three closely related fungi to identify crucial transitory phases of intron gain and loss. We found 74 intron positions showing intraspecific presence-absence polymorphisms (PAPs) for the entire intron. Population genetic analyses identified intron PAPs at different stages of fixation and showed that intron gain or loss was very recent. We found direct support for extensive intron transposition among unrelated genes. A substantial proportion of highly similar introns in the genome either were recently gained or showed a transient phase of intron PAP. We also identified an intron transfer among paralogous genes that created a new intron. Intron loss was due mainly to homologous recombination involving reverse-transcribed mRNA. The large number of intron positions in transient phases of either intron gain or loss shows that intron evolution is much faster than previously thought and provides an excellent model to study molecular mechanisms of intron gain.

Intraspecific comparison and annotation of two complete mitochondrial genome sequences from the plant pathogenic fungus Mycosphaerella graminicola

The mitochondrial genomes of two isolates of the wheat pathogen Mycosphaerella graminicola were sequenced completely and compared to identify polymorphic regions. This organism is of interest because it is phylogenetically distant from other fungi with sequenced mitochondrial genomes and it has shown discordant patterns of nuclear and mitochondrial diversity. The mitochondrial genome of M. graminicola is a circular molecule of approximately 43,960bp containing the typical genes coding for 14 proteins related to oxidative phosphorylation, one RNA polymerase, two rRNA genes and a set of 27 tRNAs. The mitochondrial DNA of M. graminicola lacks the gene encoding the putative ribosomal protein (rps5-like), commonly found in fungal mitochondrial genomes. Most of the tRNA genes were clustered with a gene order conserved with many other ascomycetes. A sample of 35 additional strains representing the known global mt diversity was partially sequenced to measure overall mitochondrial variability within the species. Little variation was found, confirming previous RFLP-based findings of low mitochondrial diversity. The mitochondrial sequence of M. graminicola is the first reported from the family Mycosphaerellaceae or the order Capnodiales. The sequence also provides a tool to better understand the development of fungicide resistance and the conflicting pattern of high nuclear and low mitochondrial diversity in global populations of this fungus.

Comparative analysis of mitochondrial genomes from closely related Rhynchosporium species reveals extensive intron invasion.

We sequenced and annotated the complete mitochondrial (mt) genomes of four closely related Rhynchosporium species that diverged ∼14,000-35,000years ago. During this time frame, three of the mt genomes expanded significantly due to an invasion of introns into three genes (cox1, cox2, and nad5). The enlarged mt genomes contained ∼40% introns compared to 8.1% in uninvaded relatives. Many intron gains were accompanied by co-conversion of flanking exonic regions. The comparative analysis revealed a highly variable set of non-intronic, free-standing ORFs of unknown function (uORFs). This is consistent with a rapidly evolving accessory compartment in the mt genome of these closely related species. Only one free-standing uORF was shared among all mt genomes analyzed. This uORF had a mutation rate similar to the core mt protein-encoding genes, suggesting conservation of function among the species. The nucleotide composition of the core protein-encoding genes significantly differed from those of introns and uORFs. The mt mutation rate was 77 times higher than the nuclear mutation rate, indicating that the phylogeny inferred from mt genes may better resolve the phylogenetic relationships among closely related Rhynchosporium species than phylogenies inferred from nuclear genes.

Microbial mediation of complex subterranean mineral structures

Helictites—an enigmatic type of mineral structure occurring in some caves—differ from classical speleothems as they develop with orientations that defy gravity. While theories for helictite formation have been forwarded, their genesis remains equivocal. Here, we show that a remarkable suite of helictites occurring in Asperge Cave (France) are formed by biologically-mediated processes, rather than abiotic processes as had hitherto been proposed. Morphological and petro-physical properties are inconsistent with mineral precipitation under purely physico-chemical control. Instead, microanalysis and molecular-biological investigation reveals the presence of a prokaryotic biofilm intimately associated with the mineral structures. We propose that microbially-influenced mineralization proceeds within a gliding biofilm which serves as a nucleation site for CaCO3, and where chemotaxis influences the trajectory of mineral growth, determining the macroscopic morphology of the speleothems. The influence of biofilms may explain the occurrence of similar speleothems in other caves worldwide, and sheds light on novel biomineralization processes.

Comparative Transcriptome Analyses in Zymoseptoria tritici Reveal Significant Differences in Gene Expression Among Strains During Plant Infection

Zymoseptoria tritici is an ascomycete fungus that causes Septoria tritici blotch, a globally distributed foliar disease on wheat. Z. tritici populations are highly polymorphic and exhibit significant quantitative variation for virulence. Despite its importance, the genes responsible for quantitative virulence in this pathogen remain largely unknown. We investigated the expression profiles of four Z. tritici strains differing in virulence in an experiment conducted under uniform environmental conditions. Transcriptomes were compared at four different infection stages to characterize the regulation of gene families thought to be involved in virulence and to identify new virulence factors. The major components of the fungal infection transcriptome showed consistent expression profiles across strains. However, strain-specific regulation was observed for many genes, including some encoding putative virulence factors. We postulate that strain-specific regulation of virulence factors can determine the outcome of Z. tritici infections. We show that differences in gene expression may be major determinants of virulence variation among Z. tritici strains, adding to the already known contributions to virulence variation based on differences in gene sequence and gene presence/absence polymorphisms.

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. Wrights F statistics (FST) 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.