Telomeres in the rice blast fungus Magnaporthe oryzae: the world of the end as we know it

First telomere-to-telomere gapless assembly of the rice blast fungus Pyricularia oryzae

Rice blast caused by Pyricularia oryzae (syn., Magnaporthe oryzae) was one of the most destructive diseases of rice throughout the world. Genome assembly was fundamental to genetic variation identification and critically impacted the understanding of its ability to overcome host resistance. Here, we report a gapless genome assembly of rice blast fungus P. oryzae strain P131 using PacBio, Illumina and high throughput chromatin conformation capture (Hi-C) sequencing data. This assembly contained seven complete chromosomes (43,237,743 bp) and a circular mitochondrial genome (34,866 bp). Approximately 14.31% of this assembly carried repeat sequences, significantly greater than its previous assembled version. This assembly had a 99.9% complement in BUSCO evaluation. A total of 14,982 genes protein-coding genes were predicted. In summary, we assembled the first telomere-to-telomere gapless genome of P. oryzae, which would be a valuable genome resource for future research on the genome evolution and host adaptation.

Gapless reference genome assembly of Didymella glomerata, a new fungal pathogen of maize causing Didymella leaf blight

Didymella leaf blight (DLB) caused by Didymella glomerata is a new fungal disease of maize (Zea mays), first detected in 2021 in Panjin, Liaoning province of China. Here we report the reference genome assembly of D. glomerata to unravel how the fungal pathogen controls its virulence on maize at the molecular level. A maize-infecting strain Pj-2 of the pathogen was sequenced on the Illumina NovaSeq 6000 and PacBio Sequel II platforms at a 575-fold genomic coverage. The 33.17 Mb gapless genome assembly comprises 32 scaffolds with L/N₅₀ of 11/1.36 Mb, four of which represent full-length chromosomes. The Pj-2 genome is predicted to contain 10,334 protein-coding genes, of which 211, 12 and 134 encode effector candidates, secondary metabolite backbone-forming enzymes and CAZymes, respectively. Some of these genes are potentially implicated in niche adaptation and expansion, such as colonizing new hosts like maize. Phylogenomic analysis of eight strains of six Didymella spp., including three sequenced strains of D. glomerata, reveals that the maize (Pj-2)- and Chrysanthemum (CBS 528.66)-infecting strains of D. glomerata are genetically similar (sharing 92.37% genome with 98.89% identity), whereas Pj-2 shows truncated collinearity with extensive chromosomal rearrangements with the Malus-infecting strain M27-16 of D. glomerata (sharing only 55.01% genome with 88.20% identity). Pj-2 and CBS 528.66 carry four major reciprocal translocations in their genomes, which may enable them to colonize the different hosts. Furthermore, germplasm screening against Pj-2 led to the identification of three sources of DLB resistance in maize, including a tropical inbred line CML496. DLB resistance in the line is attributed to the accumulation of ROS H₂O₂ in the apoplastic space of the infected cells, which likely restricts the fungal growth and proliferation.

How to Completely Squeeze a Fungus-Advanced Genome Mining Tools for Novel Bioactive Substances

Fungal species have the capability of producing an overwhelming diversity of bioactive substances that can have beneficial but also detrimental effects on human health. These so-called secondary metabolites naturally serve as antimicrobial "weapon systems", signaling molecules or developmental effectors for fungi and hence are produced only under very specific environmental conditions or stages in their life cycle. However, as these complex conditions are difficult or even impossible to mimic in laboratory settings, only a small fraction of the true chemical diversity of fungi is known so far. This also implies that a large space for potentially new pharmaceuticals remains unexplored. We here present an overview on current developments in advanced methods that can be used to explore this chemical space. We focus on genetic and genomic methods, how to detect genes that harbor the blueprints for the production of these compounds (i.e., biosynthetic gene clusters, BGCs), and ways to activate these silent chromosomal regions. We provide an in-depth view of the chromatin-level regulation of BGCs and of the potential to use the CRISPR/Cas technology as an activation tool.

Telomere Roles in Fungal Genome Evolution and Adaptation

Telomeres form the ends of linear chromosomes and usually comprise protein complexes that bind to simple repeated sequence motifs that are added to the 3′ ends of DNA by the telomerase reverse transcriptase (TERT). One of the primary functions attributed to telomeres is to solve the “end-replication problem” which, if left unaddressed, would cause gradual, inexorable attrition of sequences from the chromosome ends and, eventually, loss of viability. Telomere-binding proteins also protect the chromosome from 5′ to 3′ exonuclease action, and disguise the chromosome ends from the double-strand break repair machinery whose illegitimate action potentially generates catastrophic chromosome aberrations. Telomeres are of special interest in the blast fungus, Pyricularia, because the adjacent regions are enriched in genes controlling interactions with host plants, and the chromosome ends show enhanced polymorphism and genetic instability. Previously, we showed that telomere instability in some P. oryzae strains is caused by novel retrotransposons (MoTeRs) that insert in telomere repeats, generating interstitial telomere sequences that drive frequent, break-induced rearrangements. Here, we sought to gain further insight on telomeric involvement in shaping Pyricularia genome architecture by characterizing sequence polymorphisms at chromosome ends, and surrounding internalized MoTeR loci (relics) and interstitial telomere repeats. This provided evidence that telomere dynamics have played historical, and likely ongoing, roles in shaping the Pyricularia genome. We further demonstrate that even telomeres lacking MoTeR insertions are poorly preserved, such that the telomere-adjacent sequences exhibit frequent presence/absence polymorphism, as well as exchanges with the genome interior. Using TERT knockout experiments, we characterized chromosomal responses to failed telomere maintenance which suggested that much of the MoTeR relic-/interstitial telomere-associated polymorphism could be driven by compromised telomere function. Finally, we describe three possible examples of a phenomenon known as “Adaptive Telomere Failure,” where spontaneous losses of telomere maintenance drive rapid accumulation of sequence polymorphism with possible adaptive advantages. Together, our data suggest that telomere maintenance is frequently compromised in Pyricularia but the chromosome alterations resulting from telomere failure are not as catastrophic as prior research would predict, and may, in fact, be potent drivers of adaptive polymorphism.

Use of Telomere Fingerprinting to Identify Clonal Lineages of Colletotrichum fioriniae in Kentucky Mixed-Fruit Orchards

Multiple species in the fungal genus Colletotrichum cause anthracnose fruit rot diseases that are responsible for major yield losses of as much as 100%. Individual species of Colletotrichum typically have broad host ranges and can infect multiple fruit species. Colletotrichum fioriniae causes anthracnose fruit rots of apples, blueberries, and strawberries in Kentucky orchards where these fruits grow in close proximity. This raises the possibility of cross-infection, which may have significant management implications. The potential occurrence of cross-infection was investigated by using telomere fingerprinting to identify C. fioriniae clones in several mixed-fruit orchards. Telomere fingerprints were highly polymorphic among a test group of C. fioriniae strains and effectively defined clonal lineages. Fingerprints were compared among apple, blueberry, and strawberry isolates of C. fioriniae from three different orchards and similarity matrices were calculated to build phylograms for each orchard group. Multiple clonal lineages of C. fioriniae were identified within each orchard on the same fruit host. Related lineages were found among isolates from different hosts, but the results did not provide direct evidence for cross-infection of different fruit species by the same clones. Recovery of the same clonal lineages within orchards across multiple years suggested that local dispersal was important in pathogen population structure and that C. fioriniae strains persisted within orchards over time. Isolates from blueberry were less diverse than isolates from apple, perhaps related to more intensive anthracnose management protocols on apple versus blueberry. Telomere fingerprinting is a valuable tool for understanding population dynamics of Colletotrichum fruit rot fungi.

Reference genome assembly for Australian Ascochyta lentis isolate Al4

Ascochyta lentis causes ascochyta blight in lentil (Lens culinaris Medik.) and yield loss can be as high as 50%. With careful agronomic management practices, fungicide use, and advances in breeding resistant lentil varieties, disease severity and impact to farmers have been largely controlled. However, evidence from major lentil producing countries, Canada and Australia, suggests that A. lentis isolates can change their virulence profile and level of aggressiveness over time and under different selection pressures. In this paper, we describe the first genome assembly for A. lentis for the Australian isolate Al4, through the integration of data from Illumina and PacBio SMRT sequencing. The Al4 reference genome assembly is almost 42 Mb in size and encodes 11,638 predicted genes. The Al4 genome comprises 21 full-length and gapless chromosomal contigs and two partial chromosome contigs each with one telomere. We predicted 31 secondary metabolite clusters, and 38 putative protein effectors, many of which were classified as having an unknown function. Comparison of A. lentis genome features with the recently published reference assembly for closely related A. rabiei show that genome synteny between these species is highly conserved. However, there are several translocations and inversions of genome sequence. The location of secondary metabolite clusters near transposable element and repeat-rich genomic regions was common for A. lentis as has been reported for other fungal plant pathogens.

Keeping Balance Between Genetic Stability and Plasticity at the Telomere and Subtelomere of Trypanosoma brucei

Telomeres, the nucleoprotein complexes at chromosome ends, are well-known for their essential roles in genome integrity and chromosome stability. Yet, telomeres and subtelomeres are frequently less stable than chromosome internal regions. Many subtelomeric genes are important for responding to environmental cues, and subtelomeric instability can facilitate organismal adaptation to extracellular changes, which is a common theme in a number of microbial pathogens. In this review, I will focus on the delicate and important balance between stability and plasticity at telomeres and subtelomeres of a kinetoplastid parasite, Trypanosoma brucei, which causes human African trypanosomiasis and undergoes antigenic variation to evade the host immune response. I will summarize the current understanding about T. brucei telomere protein complex, the telomeric transcript, and telomeric R-loops, focusing on their roles in maintaining telomere and subtelomere stability and integrity. The similarities and differences in functions and underlying mechanisms of T. brucei telomere factors will be compared with those in human and yeast cells.

Chromosome-level genome assembly and manually-curated proteome of model necrotroph Parastagonospora nodorum Sn15 reveals a genome-wide trove of candidate effector homologs, and redundancy of virulence-related functions within an accessory chromosome

Background

The fungus Parastagonospora nodorum causes septoria nodorum blotch (SNB) of wheat (Triticum aestivum) and is a model species for necrotrophic plant pathogens. The genome assembly of reference isolate Sn15 was first reported in 2007. P. nodorum infection is promoted by its production of proteinaceous necrotrophic effectors, three of which are characterised – ToxA, Tox1 and Tox3.

Results

A chromosome-scale genome assembly of P. nodorum Australian reference isolate Sn15, which combined long read sequencing, optical mapping and manual curation, produced 23 chromosomes with 21 chromosomes possessing both telomeres. New transcriptome data were combined with fungal-specific gene prediction techniques and manual curation to produce a high-quality predicted gene annotation dataset, which comprises 13,869 high confidence genes, and an additional 2534 lower confidence genes retained to assist pathogenicity effector discovery. Comparison to a panel of 31 internationally-sourced isolates identified multiple hotspots within the Sn15 genome for mutation or presence-absence variation, which was used to enhance subsequent effector prediction. Effector prediction resulted in 257 candidates, of which 98 higher-ranked candidates were selected for in-depth analysis and revealed a wealth of functions related to pathogenicity. Additionally, 11 out of the 98 candidates also exhibited orthology conservation patterns that suggested lateral gene transfer with other cereal-pathogenic fungal species. Analysis of the pan-genome indicated the smallest chromosome of 0.4 Mbp length to be an accessory chromosome (AC23). AC23 was notably absent from an avirulent isolate and is predominated by mutation hotspots with an increase in non-synonymous mutations relative to other chromosomes. Surprisingly, AC23 was deficient in effector candidates, but contained several predicted genes with redundant pathogenicity-related functions.

Conclusions

We present an updated series of genomic resources for P. nodorum Sn15 – an important reference isolate and model necrotroph – with a comprehensive survey of its predicted pathogenicity content.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-021-07699-8.

The Impact of Blast Disease: Past, Present, and Future

Rice blast disease is both the most explosive and potentially damaging disease of the world's rice (Oryza sativa) crop and a model system for research on the molecular mechanisms that fungi use to cause plant disease. The blast fungus, Magnaporthe oryzae, is highly evolved to sense when it is on a leaf surface; to develop a pressurized cell, the appressorium, to punch through the leaf cuticle; and then to hijack living rice cells to assist it in causing disease. Host specificity, determining which plants particular fungal strains can infect, is also an important topic for research. The blast fungus is a moving target, quickly overcoming rice resistance genes we deploy to control it, and recently emerging to cause devastating disease on an entirely new cereal crop, wheat. M. oryzae is highly adaptable, with multiple examples of genetic instability at certain gene loci and in certain genomic regions. Understanding the biology of the fungus in the field, and its potential for genetic and genome variability, is key to keep it from adapting to life in the research laboratory and losing relevance to the significant impact it has on global food security.

Transposon-mediated telomere destabilization: a driver of genome evolution in the blast fungus

The fungus Magnaporthe oryzae causes devastating diseases of crops, including rice and wheat, and in various grasses. Strains from ryegrasses have highly unstable chromosome ends that undergo frequent rearrangements, and this has been associated with the presence of retrotransposons (Magnaporthe oryzae Telomeric Retrotransposons-MoTeRs) inserted in the telomeres. The objective of the present study was to determine the mechanisms by which MoTeRs promote telomere instability. Targeted cloning, mapping, and sequencing of parental and novel telomeric restriction fragments (TRFs), along with MinION sequencing of genomic DNA allowed us to document the precise molecular alterations underlying 109 newly-formed TRFs. These included truncations of subterminal rDNA sequences; acquisition of MoTeR insertions by 'plain' telomeres; insertion of the MAGGY retrotransposons into MoTeR arrays; MoTeR-independent expansion and contraction of subtelomeric tandem repeats; and a variety of rearrangements initiated through breaks in interstitial telomere tracts that are generated during MoTeR integration. Overall, we estimate that alterations occurred in approximately sixty percent of chromosomes (one in three telomeres) analyzed. Most importantly, we describe an entirely new mechanism by which transposons can promote genomic alterations at exceptionally high frequencies, and in a manner that can promote genome evolution while minimizing collateral damage to overall chromosome architecture and function.

Genome comparisons suggest an association between Ceratocystis host adaptations and effector clusters in unique transposable element families

Ceratocystis fimbriata is a host specific fungal pathogen of sweet potato (Ipomoea batatas). The closely related species, C. manginecans, is an important pathogen of trees (e.g. Acacia mangium and Mangifera indica) but has never been isolated from tuber crops. The genetic factors that determine the host range and host specificity of these species have not been determined. The aim of this study was to compare the genomes of C. fimbriata and C. manginecans in order to identify species-specific genetic differences that could be associated with host specificity. This included whole-genome alignments as well as comparisons of gene content and transposable elements (TEs). The genomes of the two species were found to be very similar, sharing similar catalogues of CAZymes, peptidases and lipases. However, the genomes of the two species also varied, harbouring species-specific genes (e.g. small secreted effectors, nutrient processing proteins and stress response proteins). A portion of the TEs identified (17%) had a unique distribution in each species. Transposable elements appeared to have played a prominent role in the divergence of the two species because they were strongly associated with chromosomal translocations and inversions as well as with unique genomic regions containing species-specific genes. Two large effector clusters, with unique TEs in each species, were identified. These effectors displayed non-synonymous mutations and deletions, conserved within a species, and could serve as mutational hot-spots for the development of host specificity in the two species.

Reference Genome Assembly for Australian Ascochyta rabiei Isolate ArME14

Ascochyta rabiei is the causal organism of ascochyta blight of chickpea and is present in chickpea crops worldwide. Here we report the release of a high-quality PacBio genome assembly for the Australian A. rabiei isolate ArME14. We compare the ArME14 genome assembly with an Illumina assembly for Indian A. rabiei isolate, ArD2. The ArME14 assembly has gapless sequences for nine chromosomes with telomere sequences at both ends and 13 large contig sequences that extend to one telomere. The total length of the ArME14 assembly was 40,927,385 bp, which was 6.26 Mb longer than the ArD2 assembly. Division of the genome by OcculterCut into GC-balanced and AT-dominant segments reveals 21% of the genome contains gene-sparse, AT-rich isochores. Transposable elements and repetitive DNA sequences in the ArME14 assembly made up 15% of the genome. A total of 11,257 protein-coding genes were predicted compared with 10,596 for ArD2. Many of the predicted genes missing from the ArD2 assembly were in genomic regions adjacent to AT-rich sequence. We compared the complement of predicted transcription factors and secreted proteins for the two A. rabiei genome assemblies and found that the isolates contain almost the same set of proteins. The small number of differences could represent real differences in the gene complement between isolates or possibly result from the different sequencing methods used. Prediction pipelines were applied for carbohydrate-active enzymes, secondary metabolite clusters and putative protein effectors. We predict that ArME14 contains between 450 and 650 CAZymes, 39 putative protein effectors and 26 secondary metabolite clusters.

Telomere and Subtelomere R-loops and Antigenic Variation in Trypanosomes

Trypanosoma brucei is a kinetoplastid parasite that causes African trypanosomiasis, which is fatal if left untreated. T. brucei regularly switches its major surface antigen, VSG, to evade the host immune responses. VSGs are exclusively expressed from subtelomeric expression sites (ESs) where VSG genes are flanked by upstream 70 bp repeats and downstream telomeric repeats. The telomere downstream of the active VSG is transcribed into a long-noncoding RNA (TERRA), which forms RNA:DNA hybrids (R-loops) with the telomeric DNA. At an elevated level, telomere R-loops cause more telomeric and subtelomeric double-strand breaks (DSBs) and increase VSG switching rate. In addition, stabilized R-loops are observed at the 70 bp repeats and immediately downstream of ES-linked VSGs in RNase H defective cells, which also have an increased amount of subtelomeric DSBs and more frequent VSG switching. Although subtelomere plasticity is expected to be beneficial to antigenic variation, severe defects in subtelomere integrity and stability increase cell lethality. Therefore, regulation of the telomere and 70 bp repeat R-loop levels is important for the balance between antigenic variation and cell fitness in T. brucei. In addition, the high level of the active ES transcription favors accumulation of R-loops at the telomere and 70 bp repeats, providing an intrinsic mechanism for local DSB formation, which is a strong inducer of VSG switching.

In silico prediction and characterisation of secondary metabolite clusters in the plant pathogenic fungus Verticillium dahliae

Fungi are renowned producers of natural compounds, also known as secondary metabolites (SMs) that display a wide array of biological activities. Typically, the genes that are involved in the biosynthesis of SMs are located in close proximity to each other in so-called secondary metabolite clusters. Many plant-pathogenic fungi secrete SMs during infection in order to promote disease establishment, for instance as cytocoxic compounds. Verticillium dahliae is a notorious plant pathogen that can infect over 200 host plants worldwide. However, the SM repertoire of this vascular pathogen remains mostly uncharted. To unravel the potential of V. dahliae to produce SMs, we performed in silico predictions and in-depth analyses of its secondary metabolite clusters. Using distinctive traits of gene clusters and the conserved signatures of core genes 25 potential SM gene clusters were identified. Subsequently, phylogenetic and comparative genomics analyses were performed, revealing that two putative siderophores, ferricrocin and TAFC, DHN-melanin and fujikurin may belong to the SM repertoire of V. dahliae.

Opposite Polarity Monospore Genome De Novo Sequencing and Comparative Analysis Reveal the Possible Heterothallic Life Cycle of Morchella importuna

Morchella is a popular edible fungus worldwide due to its rich nutrition and unique flavor. Many research efforts were made on the domestication and cultivation of Morchella all over the world. In recent years, the cultivation of Morchella was successfully commercialized in China. However, the biology is not well understood, which restricts the further development of the morel fungus cultivation industry. In this paper, we performed de novo sequencing and assembly of the genomes of two monospores with a different mating type (M04M24 and M04M26) isolated from the commercially cultivated strain M04. Gene annotation and comparative genome analysis were performed to study differences in CAZyme (Carbohydrate-active enzyme) enzyme content, transcription factors, duplicated sequences, structure of mating type sites, and differences at the gene and functional levels between the two monospore strains of M. importuna. Results showed that the de novo assembled haploid M04M24 and M04M26 genomes were 48.98 and 51.07 Mb, respectively. A complete fine physical map of M. importuna was obtained from genome coverage and gene completeness evaluation. A total of 10,852 and 10,902 common genes and 667 and 868 endemic genes were identified from the two monospore strains, respectively. The Gene Ontology (GO) and KAAS (KEGG Automatic Annotation Serve) enrichment analyses showed that the endemic genes performed different functions. The two monospore strains had 99.22% collinearity with each other, accompanied with certain position and rearrangement events. Analysis of complete mating-type loci revealed that the two monospore M. importuna strains contained an independent mating-type structure and remained conserved in sequence and location. The phylogenetic and divergence time of M. importuna was analyzed at the whole-genome level for the first time. The bifurcation time of morel and tuber was estimated to be 201.14 million years ago (Mya); the two monospore strains with a different mating type represented the evolution of different nuclei, and the single copy homologous genes between them were also different due to a genetic differentiation distance about 0.65 Mya. Compared with truffles, M. importuna had an extension of 28 clusters of orthologous genes (COGs) and a contraction of two COGs. The two different polar nuclei with different degrees of contraction and expansion suggested that they might have undergone different evolutionary processes. The different mating-type structures, together with the functional clustering and enrichment analysis results of the endemic genes of the two different polar nuclei, imply that M. importuna might be a heterothallic fungus and the interaction between the endemic genes may be necessary for its complete life history. Studies on the genome of M. importuna facilitate a better understanding of morel biology and evolution.

Telomeres and genomic evolution

The terminal regions of eukaryotic chromosomes, composed of telomere repeat sequences and sub-telomeric sequences, represent some of the most variable and rapidly evolving regions of the genome. The sub-telomeric regions are characterized by segmentally duplicated repetitive DNA elements, interstitial telomere repeat sequences and families of variable genes. Sub-telomeric repeat sequence families are shared among multiple chromosome ends, often rendering detailed sequence characterization difficult. These regions are composed of constitutive heterochromatin and are subjected to high levels of meiotic recombination. Dysfunction within telomere repeat arrays, either due to disruption in the chromatin structure or because of telomere shortening, can lead to chromosomal fusion and the generation of large-scale genomic rearrangements across the genome. The dynamic nature of telomeric regions, therefore, provides functionally useful variation to create genetic diversity, but also provides a mechanism for rapid genomic evolution that can lead to reproductive isolation and speciation. This article is part of the theme issue 'Understanding diversity in telomere dynamics'.This article is part of the theme issue 'Understanding diversity in telomere dynamics'.

In silico prediction and characterization of secondary metabolite biosynthetic gene clusters in the wheat pathogen Zymoseptoria tritici

BACKGROUND

Fungal pathogens of plants produce diverse repertoires of secondary metabolites, which have functions ranging from iron acquisition, defense against immune perturbation, to toxic assaults on the host. The wheat pathogen Zymoseptoria tritici causes Septoria tritici blotch, a foliar disease which is a significant threat to global food security. Currently, there is limited knowledge of the secondary metabolite arsenal produced by Z. tritici, which significantly restricts mechanistic understanding of infection. In this study, we analyzed the genome of Z. tritici isolate IP0323 to identify putative secondary metabolite biosynthetic gene clusters, and used comparative genomics to predict their encoded products.

RESULTS

We identified 32 putative secondary metabolite clusters. These were physically enriched at subtelomeric regions, which may facilitate diversification of cognate products by rapid gene rearrangement or mutations. Comparative genomics revealed a four gene cluster with significant similarity to the ferrichrome-A biosynthetic locus of the maize pathogen Ustilago maydis, suggesting this siderophore is deployed by Z. tritici to acquire iron. The Z. tritici genome also contains several isoprenoid biosynthetic gene clusters, including one with high similarity to a carotenoid/opsin producing locus in several fungi. Furthermore, we identify putative phytotoxin biosynthetic clusters, suggesting Z. tritici can produce an epipolythiodioxopiperazine, and a polyketide and non-ribosomal peptide with predicted structural similarities to fumonisin and the Alternaria alternata AM-toxin, respectively. Interrogation of an existing transcriptional dataset suggests stage specific deployment of numerous predicted loci during infection, indicating an important role of these secondary metabolites in Z. tritici disease.

CONCLUSIONS

We were able to assign putative biosynthetic products to numerous clusters based on conservation amongst other fungi. However, analysis of the majority of secondary metabolite loci did not enable prediction of a cluster product, and consequently the capacity of these loci to play as yet undetermined roles in disease or other stages of the Z. tritici lifecycle is significant. These data will drive future experimentation for determining the role of these clusters and cognate secondary metabolite products in Z. tritici virulence, and may lead to discovery of novel bioactive molecules.

Small RNA Functions Are Required for Growth and Development of Magnaporthe oryzae

RNA interference (RNAi) is conserved in eukaryotic organisms, and it has been well studied in many animal and plant species and some fungal species, yet it is not well studied in fungal plant pathogens. In the rice blast fungus Magnaporthe oryzae, we examined small RNA (sRNA) and their biogenesis in the context of growth and pathogenicity. Through genetic and genomic analyses, we demonstrate that loss of a single gene encoding Dicer, RNA-dependent RNA polymerase, or Argonaute reduces sRNA levels. These three proteins are required for the biogenesis of sRNA-matching genome-wide regions (coding regions, repeats, and intergenic regions). The loss of one Argonaute reduced both sRNA and fungal virulence on barley leaves. Transcriptome analysis of multiple mutants revealed that sRNA play an important role in transcriptional regulation of repeats and intergenic regions in M. oryzae. Together, these data support that M. oryzae sRNA regulate developmental processes including, fungal growth and virulence.

Genome analysis and avirulence gene cloning using a high-density RADseq linkage map of the flax rust fungus, Melampsora lini

BACKGROUND

Rust fungi are an important group of plant pathogens that cause devastating losses in agricultural, silvicultural and natural ecosystems. Plants can be protected from rust disease by resistance genes encoding receptors that trigger a highly effective defence response upon recognition of specific pathogen avirulence proteins. Identifying avirulence genes is crucial for understanding how virulence evolves in the field.

RESULTS

To facilitate avirulence gene cloning in the flax rust fungus, Melampsora lini, we constructed a high-density genetic linkage map using single nucleotide polymorphisms detected in restriction site-associated DNA sequencing (RADseq) data. The map comprises 13,412 RADseq markers in 27 linkage groups that together span 5860 cM and contain 2756 recombination bins. The marker sequences were used to anchor 68.9 % of the M. lini genome assembly onto the genetic map. The map and anchored assembly were then used to: 1) show that M. lini has a high overall meiotic recombination rate, but recombination distribution is uneven and large coldspots exist; 2) show that substantial genome rearrangements have occurred in spontaneous loss-of-avirulence mutants; and 3) identify the AvrL2 and AvrM14 avirulence genes by map-based cloning. AvrM14 is a dual-specificity avirulence gene that encodes a predicted nudix hydrolase. AvrL2 is located in the region of the M. lini genome with the lowest recombination rate and encodes a small, highly-charged proline-rich protein.

CONCLUSIONS

The M. lini high-density linkage map has greatly advanced our understanding of virulence mechanisms in this pathogen by providing novel insights into genome variability and enabling identification of two new avirulence genes.

OcculterCut: A Comprehensive Survey of AT-Rich Regions in Fungal Genomes

We present a novel method to measure the local GC-content bias in genomes and a survey of published fungal species. The method, enacted as "OcculterCut" (https://sourceforge.net/projects/occultercut, last accessed April 30, 2016), identified species containing distinct AT-rich regions. In most fungal taxa, AT-rich regions are a signature of repeat-induced point mutation (RIP), which targets repetitive DNA and decreases GC-content though the conversion of cytosine to thymine bases. RIP has in turn been identified as a driver of fungal genome evolution, as RIP mutations can also occur in single-copy genes neighboring repeat-rich regions. Over time RIP perpetuates "two speeds" of gene evolution in the GC-equilibrated and AT-rich regions of fungal genomes. In this study, genomes showing evidence of this process are found to be common, particularly among the Pezizomycotina. Further analysis highlighted differences in amino acid composition and putative functions of genes from these regions, supporting the hypothesis that these regions play an important role in fungal evolution. OcculterCut can also be used to identify genes undergoing RIP-assisted diversifying selection, such as small, secreted effector proteins that mediate host-microbe disease interactions.

A game of hide and seek between avirulence genes AvrLm4-7 and AvrLm3 in Leptosphaeria maculans

Extending the durability of plant resistance genes towards fungal pathogens is a major challenge. We identified and investigated the relationship between two avirulence genes of Leptosphaeria maculans, AvrLm3 and AvrLm4-7. When an isolate possesses both genes, the Rlm3-mediated resistance of oilseed rape (Brassica napus) is not expressed due to the presence of AvrLm4-7 but virulent isolates toward Rlm7 recover the AvrLm3 phenotype. Combining genetic and genomic approaches (genetic mapping, RNA-seq, BAC (bacterial artificial chromosome) clone sequencing and de novo assembly) we cloned AvrLm3, a telomeric avirulence gene of L. maculans. AvrLm3 is located in a gap of the L. maculans reference genome assembly, is surrounded by repeated elements, encodes for a small secreted cysteine-rich protein and is highly expressed at early infection stages. Complementation and silencing assays validated the masking effect of AvrLm4-7 on AvrLm3 recognition by Rlm3 and we showed that the presence of AvrLm4-7 does not impede AvrLm3 expression in planta. Y2H assays suggest the absence of physical interaction between the two avirulence proteins. This unusual interaction is the basis for field experiments aiming to evaluate strategies that increase Rlm7 durability.

A complete annotation of the chromosomes of the cellulase producer Trichoderma reesei provides insights in gene clusters, their expression and reveals genes required for fitness

BACKGROUND

Investigations on a few eukaryotic model organisms showed that many genes are non-randomly distributed on chromosomes. In addition, chromosome ends frequently possess genes that are important for the fitness of the organisms. Trichoderma reesei is an industrial producer of enzymes for food, feed and biorefinery production. Its seven chromosomes have recently been assembled, thus making an investigation of its chromosome architecture possible.

RESULTS

We manually annotated and mapped 9194 ORFs on their respective chromosomes and investigated the clustering of the major gene categories and of genes encoding carbohydrate-active enzymes (CAZymes), and the relationship between clustering and expression. Genes responsible for RNA processing and modification, amino acid metabolism, transcription, translation and ribosomal structure and biogenesis indeed showed loose clustering, but this had no impact on their expression. A third of the genes encoding CAZymes also occurred in loose clusters that also contained a high number of genes encoding small secreted cysteine-rich proteins. Five CAZyme clusters were located less than 50 kb apart from the chromosome ends. These genes exhibited the lowest basal (but not induced) expression level, which correlated with an enrichment of H3K9 methylation in the terminal 50 kb areas indicating gene silencing. No differences were found in the expression of CAZyme genes present in other parts of the chromosomes. The putative subtelomeric areas were also enriched in genes encoding secreted proteases, amino acid permeases, enzyme clusters for polyketide synthases (PKS)-non-ribosomal peptide synthase (NRPS) fusion proteins (PKS-NRPS) and proteins involved in iron scavenging. They were strongly upregulated during conidiation and interaction with other fungi.

CONCLUSIONS

Our findings suggest that gene clustering on the T. reesei chromosomes occurs but generally has no impact on their expression. CAZyme genes, located in subtelomers, however, exhibited a much lower basal expression level. The gene inventory of the subtelomers suggests a major role of competition for nitrogen and iron supported by antibiosis for the fitness of T. reesei. The availability of fully annotated chromosomes will facilitate the use of genetic crossings in identifying still unknown genes responsible for specific traits of T. reesei.

Deciphering Genome Content and Evolutionary Relationships of Isolates from the Fungus Magnaporthe oryzae Attacking Different Host Plants

Deciphering the genetic bases of pathogen adaptation to its host is a key question in ecology and evolution. To understand how the fungus Magnaporthe oryzae adapts to different plants, we sequenced eight M. oryzae isolates differing in host specificity (rice, foxtail millet, wheat, and goosegrass), and one Magnaporthe grisea isolate specific of crabgrass. Analysis of Magnaporthe genomes revealed small variation in genome sizes (39–43 Mb) and gene content (12,283–14,781 genes) between isolates. The whole set of Magnaporthe genes comprised 14,966 shared families, 63% of which included genes present in all the nine M. oryzae genomes. The evolutionary relationships among Magnaporthe isolates were inferred using 6,878 single-copy orthologs. The resulting genealogy was mostly bifurcating among the different host-specific lineages, but was reticulate inside the rice lineage. We detected traces of introgression from a nonrice genome in the rice reference 70-15 genome. Among M. oryzae isolates and host-specific lineages, the genome composition in terms of frequencies of genes putatively involved in pathogenicity (effectors, secondary metabolism, cazome) was conserved. However, 529 shared families were found only in nonrice lineages, whereas the rice lineage possessed 86 specific families absent from the nonrice genomes. Our results confirmed that the host specificity of M. oryzae isolates was associated with a divergence between lineages without major gene flow and that, despite the strong conservation of gene families between lineages, adaptation to different hosts, especially to rice, was associated with the presence of a small number of specific gene families. All information was gathered in a public database (http://genome.jouy.inra.fr/gemo).

Chromatin-based control of effector gene expression in plant-associated fungi

Plant-associated fungi often present in their genome areas enriched in repeat sequences and effector genes, the latter being specifically induced in planta. The location of effector genes in regions enriched in repeats has been shown to have an impact on adaptability of fungi but could also provide for tight control of effector gene expression through chromatin-based regulation. The distribution of two repressive histone marks was shown to be an important regulatory layer in two fungal species with different lifestyles. Chromatin-based control of effector gene expression is likely to provide an evolutionary advantage by preventing the expression of genes not needed during vegetative growth and allow for a massive concerted expression at particular time-points of plant infection.

Complete Genome Sequence of Sporisorium scitamineum and Biotrophic Interaction Transcriptome with Sugarcane

Sporisorium scitamineum is a biotrophic fungus responsible for the sugarcane smut, a worldwide spread disease. This study provides the complete sequence of individual chromosomes of S. scitamineum from telomere to telomere achieved by a combination of PacBio long reads and Illumina short reads sequence data, as well as a draft sequence of a second fungal strain. Comparative analysis to previous available sequences of another strain detected few polymorphisms among the three genomes. The novel complete sequence described herein allowed us to identify and annotate extended subtelomeric regions, repetitive elements and the mitochondrial DNA sequence. The genome comprises 19,979,571 bases, 6,677 genes encoding proteins, 111 tRNAs and 3 assembled copies of rDNA, out of our estimated number of copies as 130. Chromosomal reorganizations were detected when comparing to sequences of S. reilianum, the closest smut relative, potentially influenced by repeats of transposable elements. Repetitive elements may have also directed the linkage of the two mating-type loci. The fungal transcriptome profiling from in vitro and from interaction with sugarcane at two time points (early infection and whip emergence) revealed that 13.5% of the genes were differentially expressed in planta and particular to each developmental stage. Among them are plant cell wall degrading enzymes, proteases, lipases, chitin modification and lignin degradation enzymes, sugar transporters and transcriptional factors. The fungus also modulates transcription of genes related to surviving against reactive oxygen species and other toxic metabolites produced by the plant. Previously described effectors in smut/plant interactions were detected but some new candidates are proposed. Ten genomic islands harboring some of the candidate genes unique to S. scitamineum were expressed only in planta. RNAseq data was also used to reassure gene predictions.

Genetics, genomics and evolution of ergot alkaloid diversity

The ergot alkaloid biosynthesis system has become an excellent model to study evolutionary diversification of specialized (secondary) metabolites. This is a very diverse class of alkaloids with various neurotropic activities, produced by fungi in several orders of the phylum Ascomycota, including plant pathogens and protective plant symbionts in the family Clavicipitaceae. Results of comparative genomics and phylogenomic analyses reveal multiple examples of three evolutionary processes that have generated ergot-alkaloid diversity: gene gains, gene losses, and gene sequence changes that have led to altered substrates or product specificities of the enzymes that they encode (neofunctionalization). The chromosome ends appear to be particularly effective engines for gene gains, losses and rearrangements, but not necessarily for neofunctionalization. Changes in gene expression could lead to accumulation of various pathway intermediates and affect levels of different ergot alkaloids. Genetic alterations associated with interspecific hybrids of Epichloë species suggest that such variation is also selectively favored. The huge structural diversity of ergot alkaloids probably represents adaptations to a wide variety of ecological situations by affecting the biological spectra and mechanisms of defense against herbivores, as evidenced by the diverse pharmacological effects of ergot alkaloids used in medicine.

Loss of a 1.6 Mb chromosome in Pyricularia oryzae harboring two alleles of AvrPik leads to acquisition of virulence to rice cultivars containing resistance alleles at the Pik locus

A small and extra chromosome of 1.6 Mb was previously identified in a Pyricularia oryzae strain, 84R-62B. To understand a role of the 1.6 Mb chromosome in the pathogenic changeability of P. oryzae, we performed experiments designed to characterize the 1.6 Mb chromosome in the present study. A gene family encoding secreted protein Pex31s in P. oryzae consists of five homologs, Pex31-A to -E. Among them, Pex31-A and -D are known to be recognized by Pik-m and Pik/Pik-m/Pik-p, respectively. In the present study, we identified Pex31-A and -D in the genome of 84R-62B. Segregation analyses using an F1 population between 84R-62B and another rice blast strain, Y93-245c-2, revealed a strong linkage between the two homologs and the 1.6 Mb chromosome of 84R-62B. A CHEF-Southern analysis revealed an association between the 1.6 Mb chromosome and the homologs, indicating that both homologs are located on the 1.6 Mb chromosome of 84R-62B. The loss of the 1.6 Mb chromosome was observed in subcultures of a F1 progeny, F1-327. These subcultures concomitantly acquired virulence on Pik, Pik-m, and Pik-p. The present study is the first report showing that loss of a small and extra chromosome leads to pathogenic mutation of P. oryzae and may provide a new insight into the mechanisms generating pathogenic variation of this fungus.

[Transposable elements reshaping genomes and favouring the evolutionary and adaptive potential of fungal phytopathogens]

Phytopathogenic fungi are a major threat for global food security and show an extreme plasticity in pathogenicity behaviours. They often have a high adaptive potential allowing them to rapidly counteract the control methods used by men in agrosystems. In this paper, we evaluate the link between genome plasticity and adaptive potential using genomics and comparative genomics approaches. Our model is the evolutionary series Leptosphaeria maculans-Leptosphaeria biglobosa, encompassing five distinct entities, whose conspecificity or heterospecificity status is unclear, and which all are pathogens of cruciferous plants. They however differ by their host range and pathogenicity. Compared to other species of the species complex, the species best adapted to oilseed rape, L. maculans "brassicae", causing important losses in the crop, has a genome that was submitted to a recent and massive burst of transposition by a few families of transposable elements (TEs). Whether the genome invasion contributed to speciation is still unclear to-date but there is a coincidence between this burst of TEs and divergence between two species. This TE burst contributed to diversification of effector proteins and thus to generation of novel pathogenic specificities. In addition, the location of effector genes within genome regions enriched in TEs has direct consequences on adaptation to plant resistance and favours a multiplicity of mutation events allowing "breakdown" of resistance. These data are substantiated by other examples in the literature showing that fungi tend to have a "two-speed" genome, in which a plastic compartment enriched in TE host genes is involved in pathogenicity and adaptation to host.

The Fusarium graminearum histone H3 K27 methyltransferase KMT6 regulates development and expression of secondary metabolite gene clusters

The cereal pathogen Fusarium graminearum produces secondary metabolites toxic to humans and animals, yet coordinated transcriptional regulation of gene clusters remains largely a mystery. By chromatin immunoprecipitation and high-throughput DNA sequencing (ChIP-seq) we found that regions with secondary metabolite clusters are enriched for trimethylated histone H3 lysine 27 (H3K27me3), a histone modification associated with gene silencing. H3K27me3 was found predominantly in regions that lack synteny with other Fusarium species, generally subtelomeric regions. Di- or trimethylated H3K4 (H3K4me2/3), two modifications associated with gene activity, and H3K27me3 are predominantly found in mutually exclusive regions of the genome. To find functions for H3K27me3, we deleted the gene for the putative H3K27 methyltransferase, KMT6, a homolog of Drosophila Enhancer of zeste, E(z). The kmt6 mutant lacks H3K27me3, as shown by western blot and ChIP-seq, displays growth defects, is sterile, and constitutively expresses genes for mycotoxins, pigments and other secondary metabolites. Transcriptome analyses showed that 75% of 4,449 silent genes are enriched for H3K27me3. A subset of genes that were enriched for H3K27me3 in WT gained H3K4me2/3 in kmt6. A largely overlapping set of genes showed increased expression in kmt6. Almost 95% of the remaining 2,720 annotated silent genes showed no enrichment for either H3K27me3 or H3K4me2/3 in kmt6. In these cases mere absence of H3K27me3 was insufficient for expression, which suggests that additional changes are required to activate genes. Taken together, we show that absence of H3K27me3 allowed expression of an additional 14% of the genome, resulting in derepression of genes predominantly involved in secondary metabolite pathways and other species-specific functions, including putative secreted pathogenicity factors. Results from this study provide the framework for novel targeted strategies to control the “cryptic genome”, specifically secondary metabolite expression.

Changes in chromatin structure are required for time- and tissue-specific gene regulation. How exactly these changes are mediated is under intense scrutiny. The interplay between activating histone modifications, e.g. H3K4me, and the silencing H3K27me3 mark has been recognized as critical to orchestrate differentiation and development in plants and animals. Here we show that filamentous fungi, exemplified by the cereal pathogen Fusarium graminearum, can use H3K27 methylation to generate silenced, facultative heterochromatin, covering more than a third of the genome, much more than the 5–8% of Neurospora or metazoan genomes. Removal of the silencing mark by mutation of the methyltransferase subunit of the PRC2 silencing complex resulted in activation of more than 1,500 genes, 14% of the genome. We show that generation of facultative heterochromatin by H3K27 methylation is an ancestral process that has been lost in certain lineages (e.g. at least some hemiascomycetes, the genus Aspergillus and some basidiomycetes). Our studies will open the door to future precise “epigenetic engineering” of gene clusters that generate bioactive compounds, e.g. putative mycotoxins, antibiotics and industrial feedstocks. Availability of tractable fungal model systems for studies of the control and function of H3K27 methylation may accelerate mechanistic research.

Cross-kingdom gene transfer facilitates the evolution of virulence in fungal pathogens

The constant interaction between plants and their pathogens has resulted in the evolution of a diverse array of microbial infection strategies. It is increasingly evident that horizontal acquisition of new virulence functions in fungi is one of the evolutionary processes that maintain pathogens' competitive edge over host plants. Genome analyses of fungi are pointing towards this phenomenon being particularly prevalent in the subphylum Pezizomycota. While the extent of cross-kingdom gene transfer can be determined with existing genomic tools and databases, so far very few horizontally transmitted genes have been functionally characterised, and an understanding of their physiological roles in virulence has been determined for even fewer genes. Understanding the evolutionary selection pressures that drive the retention of acquired genes in particular fungal lineages is important, as it will undoubtedly reveal new insights into both fungal virulence mechanisms and corresponding plant defence processes in the future.

Comparative genome structure, secondary metabolite, and effector coding capacity across Cochliobolus pathogens

The genomes of five Cochliobolus heterostrophus strains, two Cochliobolus sativus strains, three additional Cochliobolus species (Cochliobolus victoriae, Cochliobolus carbonum, Cochliobolus miyabeanus), and closely related Setosphaeria turcica were sequenced at the Joint Genome Institute (JGI). The datasets were used to identify SNPs between strains and species, unique genomic regions, core secondary metabolism genes, and small secreted protein (SSP) candidate effector encoding genes with a view towards pinpointing structural elements and gene content associated with specificity of these closely related fungi to different cereal hosts. Whole-genome alignment shows that three to five percent of each genome differs between strains of the same species, while a quarter of each genome differs between species. On average, SNP counts among field isolates of the same C. heterostrophus species are more than 25× higher than those between inbred lines and 50× lower than SNPs between Cochliobolus species. The suites of nonribosomal peptide synthetase (NRPS), polyketide synthase (PKS), and SSP-encoding genes are astoundingly diverse among species but remarkably conserved among isolates of the same species, whether inbred or field strains, except for defining examples that map to unique genomic regions. Functional analysis of several strain-unique PKSs and NRPSs reveal a strong correlation with a role in virulence.

Regulation of fungal secondary metabolism

Fungi produce a multitude of low-molecular-mass compounds known as secondary metabolites, which have roles in a range of cellular processes such as transcription, development and intercellular communication. In addition, many of these compounds now have important applications, for instance, as antibiotics or immunosuppressants. Genome mining efforts indicate that the capability of fungi to produce secondary metabolites has been substantially underestimated because many of the fungal secondary metabolite biosynthesis gene clusters are silent under standard cultivation conditions. In this Review, I describe our current understanding of the regulatory elements that modulate the transcription of genes involved in secondary metabolism. I also discuss how an improved knowledge of these regulatory elements will ultimately lead to a better understanding of the physiological and ecological functions of these important compounds and will pave the way for a novel avenue to drug discovery through targeted activation of silent gene clusters.

Genome structure and reproductive behaviour influence the evolutionary potential of a fungal phytopathogen

Modern agriculture favours the selection and spread of novel plant diseases. Furthermore, crop genetic resistance against pathogens is often rendered ineffective within a few years of its commercial deployment. Leptosphaeria maculans, the cause of phoma stem canker of oilseed rape, develops gene-for-gene interactions with its host plant, and has a high evolutionary potential to render ineffective novel sources of resistance in crops. Here, we established a four-year field experiment to monitor the evolution of populations confronted with the newly released Rlm7 resistance and to investigate the nature of the mutations responsible for virulence against Rlm7. A total of 2551 fungal isolates were collected from experimental crops of a Rlm7 cultivar or a cultivar without Rlm7. All isolates were phenotyped for virulence and a subset was genotyped with neutral genetic markers. Virulent isolates were investigated for molecular events at the AvrLm4-7 locus. Whilst virulent isolates were not found in neighbouring crops, their frequency had reached 36% in the experimental field after four years. An extreme diversity of independent molecular events leading to virulence was identified in populations, with large-scale Repeat Induced Point mutations or complete deletion of AvrLm4-7 being the most frequent. Our data suggest that increased mutability of fungal genes involved in the interactions with plants is directly related to their genomic environment and reproductive system. Thus, rapid allelic diversification of avirulence genes can be generated in L. maculans populations in a single field provided that large population sizes and sexual reproduction are favoured by agricultural practices.

Lifestyle transitions in plant pathogenic Colletotrichum fungi deciphered by genome and transcriptome analyses

Colletotrichum species are fungal pathogens that devastate crop plants worldwide. Host infection involves the differentiation of specialized cell types that are associated with penetration, growth inside living host cells (biotrophy) and tissue destruction (necrotrophy). We report here genome and transcriptome analyses of Colletotrichum higginsianum infecting Arabidopsis thaliana and Colletotrichum graminicola infecting maize. Comparative genomics showed that both fungi have large sets of pathogenicity-related genes, but families of genes encoding secreted effectors, pectin-degrading enzymes, secondary metabolism enzymes, transporters and peptidases are expanded in C. higginsianum. Genome-wide expression profiling revealed that these genes are transcribed in successive waves that are linked to pathogenic transitions: effectors and secondary metabolism enzymes are induced before penetration and during biotrophy, whereas most hydrolases and transporters are upregulated later, at the switch to necrotrophy. Our findings show that preinvasion perception of plant-derived signals substantially reprograms fungal gene expression and indicate previously unknown functions for particular fungal cell types.

Telomere-targeted retrotransposons in the rice blast fungus Magnaporthe oryzae: agents of telomere instability

The fungus Magnaporthe oryzae is a serious pathogen of rice and other grasses. Telomeric restriction fragments in Magnaporthe isolates that infect perennial ryegrass (prg) are hotspots for genomic rearrangement and undergo frequent, spontaneous alterations during fungal culture. The telomeres of rice-infecting isolates are very stable by comparison. Sequencing of chromosome ends from a number of prg-infecting isolates revealed two related non-LTR retrotransposons (M. oryzae Telomeric Retrotransposons or MoTeRs) inserted in the telomere repeats. This contrasts with rice pathogen telomeres that are uninterrupted by other sequences. Genetic evidence indicates that the MoTeR elements are responsible for the observed instability. MoTeRs represent a new family of telomere-targeted transposons whose members are found exclusively in fungi.

Genomic advances will herald new insights into the Brassica: Leptosphaeria maculans pathosystem

The study of the relationship between plants and phytopathogenic fungi is one of the most rapidly moving fields in the plant sciences, the findings of which have contributed to the development of new strategies and technologies to protect crops. Plants employ sophisticated mechanisms to perceive and appropriately defend themselves against pathogens. A good example of plant and pathogen evolution is the gene-for-gene interaction between the fungal pathogen Leptosphaeria maculans, the causal agent of blackleg disease, and Brassica crops. This interaction has been studied at the genetic and physiological level due to its agro-economic importance. The newly available genome sequence for Brassica spp. and L. maculans will provide the resources to study the co-evolution of this plant and pathogen. Particularly, an understanding of the co-evolution of genes responsible for virulence and resistance will lead to improved plant protection strategies for Brassica canola and provide a model to understand plant-pathogen interactions in other major crops. This review summarises the research-to-date in the study of the Brassica-L. maculans gene-for-gene interaction, with a focus on the genetics of resistance in Brassica and the wealth of information to be gained from genome sequencing efforts.

Microevolution of Cryptococcus neoformans driven by massive tandem gene amplification

The subtelomeric regions of organisms ranging from protists to fungi undergo a much higher rate of rearrangement than is observed in the rest of the genome. While characterizing these ~40-kb regions of the human fungal pathogen Cryptococcus neoformans, we have identified a recent gene amplification event near the right telomere of chromosome 3 that involves a gene encoding an arsenite efflux transporter (ARR3). The 3,177-bp amplicon exists in a tandem array of 2-15 copies and is present exclusively in strains with the C. neoformans var. grubii subclade VNI A5 MLST profile. Strains bearing the amplification display dramatically enhanced resistance to arsenite that correlates with the copy number of the repeat; the origin of increased resistance was verified as transport-related by functional complementation of an arsenite transporter mutant of Saccharomyces cerevisiae. Subsequent experimental evolution in the presence of increasing concentrations of arsenite yielded highly resistant strains with the ARR3 amplicon further amplified to over 50 copies, accounting for up to ~1% of the whole genome and making the copy number of this repeat as high as that seen for the ribosomal DNA. The example described here therefore represents a rare evolutionary intermediate-an array that is currently in a state of dynamic flux, in dramatic contrast to relatively common, static relics of past tandem duplications that are unable to further amplify due to nucleotide divergence. Beyond identifying and engineering fungal isolates that are highly resistant to arsenite and describing the first reported instance of microevolution via massive gene amplification in C. neoformans, these results suggest that adaptation through gene amplification may be an important mechanism that C. neoformans employs in response to environmental stresses, perhaps including those encountered during infection. More importantly, the ARR3 array will serve as an ideal model for further molecular genetic analyses of how tandem gene duplications arise and expand.

Genome evolution in the eremothecium clade of the Saccharomyces complex revealed by comparative genomics

We used comparative genomics to elucidate the genome evolution within the pre-whole-genome duplication genus Eremothecium. To this end, we sequenced and assembled the complete genome of Eremothecium cymbalariae, a filamentous ascomycete representing the Eremothecium type strain. Genome annotation indicated 4712 gene models and 143 tRNAs. We compared the E. cymbalariae genome with that of its relative, the riboflavin overproducer Ashbya (Eremothecium) gossypii, and the reconstructed yeast ancestor. Decisive changes in the Eremothecium lineage leading to the evolution of the A. gossypii genome include the reduction from eight to seven chromosomes, the downsizing of the genome by removal of 10% or 900 kb of DNA, mostly in intergenic regions, the loss of a TY3-Gypsy-type transposable element, the re-arrangement of mating-type loci, and a massive increase of its GC content. Key species-specific events are the loss of MNN1-family of mannosyltransferases required to add the terminal fourth and fifth α-1,3-linked mannose residue to O-linked glycans and genes of the Ehrlich pathway in E. cymbalariae and the loss of ZMM-family of meiosis-specific proteins and acquisition of riboflavin overproduction in A. gossypii. This reveals that within the Saccharomyces complex genome, evolution is not only based on genome duplication with subsequent gene deletions and chromosomal rearrangements but also on fungi associated with specific environments (e.g. involving fungal-insect interactions as in Eremothecium), which have encountered challenges that may be reflected both in genome streamlining and their biosynthetic potential.

Multiple translocation of the AVR-Pita effector gene among chromosomes of the rice blast fungus Magnaporthe oryzae and related species

Magnaporthe oryzae is the causal agent of rice blast disease, a devastating problem worldwide. This fungus has caused breakdown of resistance conferred by newly developed commercial cultivars. To address how the rice blast fungus adapts itself to new resistance genes so quickly, we examined chromosomal locations of AVR-Pita, a subtelomeric gene family corresponding to the Pita resistance gene, in various isolates of M. oryzae (including wheat and millet pathogens) and its related species. We found that AVR-Pita (AVR-Pita1 and AVR-Pita2) is highly variable in its genome location, occurring in chromosomes 1, 3, 4, 5, 6, 7, and supernumerary chromosomes, particularly in rice-infecting isolates. When expressed in M. oryzae, most of the AVR-Pita homologs could elicit Pita-mediated resistance, even those from non-rice isolates. AVR-Pita was flanked by a retrotransposon, which presumably contributed to its multiple translocation across the genome. On the other hand, family member AVR-Pita3, which lacks avirulence activity, was stably located on chromosome 7 in a vast majority of isolates. These results suggest that the diversification in genome location of AVR-Pita in the rice isolates is a consequence of recognition by Pita in rice. We propose a model that the multiple translocation of AVR-Pita may be associated with its frequent loss and recovery mediated by its transfer among individuals in asexual populations. This model implies that the high mobility of AVR-Pita is a key mechanism accounting for the rapid adaptation toward Pita. Dynamic adaptation of some fungal plant pathogens may be achieved by deletion and recovery of avirulence genes using a population as a unit of adaptation.

Effector diversification within compartments of the Leptosphaeria maculans genome affected by Repeat-Induced Point mutations

Fungi are of primary ecological, biotechnological and economic importance. Many fundamental biological processes that are shared by animals and fungi are studied in fungi due to their experimental tractability. Many fungi are pathogens or mutualists and are model systems to analyse effector genes and their mechanisms of diversification. In this study, we report the genome sequence of the phytopathogenic ascomycete Leptosphaeria maculans and characterize its repertoire of protein effectors. The L. maculans genome has an unusual bipartite structure with alternating distinct guanine and cytosine-equilibrated and adenine and thymine (AT)-rich blocks of homogenous nucleotide composition. The AT-rich blocks comprise one-third of the genome and contain effector genes and families of transposable elements, both of which are affected by repeat-induced point mutation, a fungal-specific genome defence mechanism. This genomic environment for effectors promotes rapid sequence diversification and underpins the evolutionary potential of the fungus to adapt rapidly to novel host-derived constraints.

In silico reversal of repeat-induced point mutation (RIP) identifies the origins of repeat families and uncovers obscured duplicated genes

Background

Repeat-induced point mutation (RIP) is a fungal genome defence mechanism guarding against transposon invasion. RIP mutates the sequence of repeated DNA and over time renders the affected regions unrecognisable by similarity search tools such as BLAST.

Results

DeRIP is a new software tool developed to predict the original sequence of a RIP-mutated region prior to the occurrence of RIP. In this study, we apply deRIP to the genome of the wheat pathogen Stagonospora nodorum SN15 and predict the origin of several previously uncharacterised classes of repetitive DNA.

Conclusions

Five new classes of transposon repeats and four classes of endogenous gene repeats were identified after deRIP. The deRIP process is a new tool for fungal genomics that facilitates the identification and understanding of the role and origin of fungal repetitive DNA. DeRIP is open-source and is available as part of the RIPCAL suite at http://www.sourceforge.net/projects/ripcal.

Evolution of linked avirulence effectors in Leptosphaeria maculans is affected by genomic environment and exposure to resistance genes in host plants

Brassica napus (canola) cultivars and isolates of the blackleg fungus, Leptosphaeria maculans interact in a 'gene for gene' manner whereby plant resistance (R) genes are complementary to pathogen avirulence (Avr) genes. Avirulence genes encode proteins that belong to a class of pathogen molecules known as effectors, which includes small secreted proteins that play a role in disease. In Australia in 2003 canola cultivars with the Rlm1 resistance gene suffered a breakdown of disease resistance, resulting in severe yield losses. This was associated with a large increase in the frequency of virulence alleles of the complementary avirulence gene, AvrLm1, in fungal populations. Surprisingly, the frequency of virulence alleles of AvrLm6 (complementary to Rlm6) also increased dramatically, even though the cultivars did not contain Rlm6. In the L. maculans genome, AvrLm1 and AvrLm6 are linked along with five other genes in a region interspersed with transposable elements that have been degenerated by Repeat-Induced Point (RIP) mutations. Analyses of 295 Australian isolates showed deletions, RIP mutations and/or non-RIP derived amino acid substitutions in the predicted proteins encoded by these seven genes. The degree of RIP mutations within single copy sequences in this region was proportional to their proximity to the degenerated transposable elements. The RIP alleles were monophyletic and were present only in isolates collected after resistance conferred by Rlm1 broke down, whereas deletion alleles belonged to several polyphyletic lineages and were present before and after the resistance breakdown. Thus, genomic environment and exposure to resistance genes in B. napus has affected the evolution of these linked avirulence genes in L. maculans.

The choanoflagellate Monosiga brevicollis karyotype revealed by the genome sequence: telomere-linked helicase genes resemble those of some fungi

The approximately 42 Mbp assembled genome sequence for the choanoflagellate Monosiga brevicollis reveals that most of the large scaffolds of 300-2,600 kb represent entire chromosomes or chromosome arms. Telomeres are partially assembled at the termini of 37 scaffolds, while another 43 scaffolds end in telomere-associated regions containing distinctive gene sets. Potential centromeric regions were identified on 39 scaffolds. Together, these observations suggest a karyotype of approximately 40 metacentric and submetacentric chromosomes averaging 1 Mbp in size. Genes encoding RecQ family DNA helicases, along with ankyrin-domain proteins and serine/threonine kinases, are associated with most telomeres, a feature shared with some fungi. This telomere-linked helicase gene arrangement might be ancestral to both fungi and choanoflagellates in the super-kingdom Opisthokonta; however, the great lability of telomere architecture suggests that it could also be a convergent feature.

The genome of Nectria haematococca: contribution of supernumerary chromosomes to gene expansion

The ascomycetous fungus Nectria haematococca, (asexual name Fusarium solani), is a member of a group of >50 species known as the "Fusarium solani species complex". Members of this complex have diverse biological properties including the ability to cause disease on >100 genera of plants and opportunistic infections in humans. The current research analyzed the most extensively studied member of this complex, N. haematococca mating population VI (MPVI). Several genes controlling the ability of individual isolates of this species to colonize specific habitats are located on supernumerary chromosomes. Optical mapping revealed that the sequenced isolate has 17 chromosomes ranging from 530 kb to 6.52 Mb and that the physical size of the genome, 54.43 Mb, and the number of predicted genes, 15,707, are among the largest reported for ascomycetes. Two classes of genes have contributed to gene expansion: specific genes that are not found in other fungi including its closest sequenced relative, Fusarium graminearum; and genes that commonly occur as single copies in other fungi but are present as multiple copies in N. haematococca MPVI. Some of these additional genes appear to have resulted from gene duplication events, while others may have been acquired through horizontal gene transfer. The supernumerary nature of three chromosomes, 14, 15, and 17, was confirmed by their absence in pulsed field gel electrophoresis experiments of some isolates and by demonstrating that these isolates lacked chromosome-specific sequences found on the ends of these chromosomes. These supernumerary chromosomes contain more repeat sequences, are enriched in unique and duplicated genes, and have a lower G+C content in comparison to the other chromosomes. Although the origin(s) of the extra genes and the supernumerary chromosomes is not known, the gene expansion and its large genome size are consistent with this species' diverse range of habitats. Furthermore, the presence of unique genes on supernumerary chromosomes might account for individual isolates having different environmental niches.

Molecular cloning and characterization of the AVR-Pia locus from a Japanese field isolate of Magnaporthe oryzae

In order to clone and analyse the avirulence gene AVR-Pia from Japanese field isolates of Magnaporthe oryzae, a mutant of the M. oryzae strain Ina168 was isolated. This mutant, which was named Ina168m95-1, gained virulence towards the rice cultivar Aichi-asahi, which contains the resistance gene Pia. A DNA fragment (named PM01) that was deleted in the mutant and that co-segregated with avirulence towards Aichi-asahi was isolated. Three cosmid clones that included the regions that flanked PM01 were isolated from a genomic DNA library. One of these clones (46F3) complemented the mutant phenotype, which indicated clearly that this clone contained the avirulence gene AVR-Pia. Clone 46F3 contained insertions of transposable elements. The 46F3 insert was divided into fragments I-VI, and these were cloned individually into a hygromycin-resistant vector for the transformation of the mutant Ina168m95-1. An inoculation assay of the transformants revealed that fragment V (3.5 kb) contained AVR-Pia. By deletion analysis of fragment V, AVR-Pia was localized to an 1199-bp DNA fragment, which included a 255-bp open reading frame with weak homology to a bacterial cytochrome-c-like protein. Restriction fragment length polymorphism analysis of this region revealed that this DNA sequence co-segregated with the AVR-Pia locus in a genetic map that was constructed using Chinese isolates.