Repeat-Induced Point Mutation and Other Genome Defense Mechanisms in Fungi

Genomic Signatures of Domestication in a Fungus Obligately Farmed by Leafcutter Ants

The naturally selected fungal crop (Leucoagaricus gongylophorus) farmed by leafcutter ants shows striking parallels with artificially selected plant crops domesticated by humans (e.g. polyploidy, engorged nutritional rewards, and dependence on cultivation). To date, poorly resolved L. gongylophorus genome assemblies based on short-read sequencing have constrained hypotheses about how millions of years under cultivation by ants shaped the fungal crop genome and potentially drove domestication. We use PacBio HiFi sequencing of L. gongylophorus from the leafcutter ant Atta colombica to identify 18 putatively novel biosynthetic gene clusters that likely cemented life as a cultivar (e.g. plant fragment degradation, ant-farmer communication, and antimicrobial defense). Comparative analyses with cultivated and free-living fungi showed genomic signatures of stepwise domestication transitions: (i) free-living to ant-cultivated: loss of genes conferring stress response and detoxification; (ii) hyphal food to engorged nutritional rewards: expansions of genes governing cellular homeostasis, carbohydrate metabolism, and siderophore biosynthesis; and (iii) detrital provisioning to freshly cut plant fragments: gene expansions promoting cell wall biosynthesis, fatty acid metabolism, and DNA repair. Comparisons across L. gongylophorus fungi farmed by 3 leafcutter ant species highlight genomic signatures of exclusively vertical clonal propagation and widespread transposable element activity. These results show how natural selection can shape domesticated cultivar genomes toward long-term ecological resilience of farming systems that have thrived across millennia.

Starships: a new frontier for fungal biology

Transposable elements (TEs) are semiautonomous genetic entities that proliferate in genomes. We recently discovered the Starships, a previously hidden superfamily of giant TEs found in a diverse subphylum of filamentous fungi, the Pezizomycotina. Starships are unlike other eukaryotic TEs because they have evolved mechanisms for both mobilizing entire genes, including those encoding conditionally beneficial phenotypes, and for horizontally transferring between individuals. We argue that Starships have unrivaled capacity to engage their fungal hosts as genetic parasites and mutualists, revealing unexplored terrain for investigating the ecoevolutionary dynamics of TE-eukaryote interactions. We build on existing models of fungal genome evolution by conceptualizing Starships as a distinct genomic compartment whose dynamics profoundly shape fungal biology.

DNA cytosine methyltransferases differentially regulate genome-wide hypermutation and interhomolog recombination in Trichoderma reesei meiosis

Trichoderma reesei is an economically important enzyme producer with several unique meiotic features. spo11, the initiator of meiotic double-strand breaks (DSBs) in most sexual eukaryotes, is dispensable for T. reesei meiosis. T. reesei lacks the meiosis-specific recombinase Dmc1. Rad51 and Sae2, the activator of the Mre11 endonuclease complex, promote DSB repair and chromosome synapsis in wild-type and spo11Δ meiosis. DNA methyltransferases (DNMTs) perform multiple tasks in meiosis. Three DNMT genes (rid1, dim2 and dimX) differentially regulate genome-wide cytosine methylation and C:G-to-T:A hypermutations in different chromosomal regions. We have identified two types of DSBs: type I DSBs require spo11 or rid1 for initiation, whereas type II DSBs do not rely on spo11 and rid1 for initiation. rid1 (but not dim2) is essential for Rad51-mediated DSB repair and normal meiosis. rid1 and rad51 exhibit a locus heterogeneity (LH) relationship, in which LH-associated proteins often regulate interconnectivity in protein interaction networks. This LH relationship can be suppressed by deleting dim2 in a haploid rid1Δ (but not rad51Δ) parental strain, indicating that dim2 and rid1 share a redundant function that acts earlier than rad51 during early meiosis. In conclusion, our studies provide the first evidence of the involvement of DNMTs during meiotic initiation and recombination.

Comparative genomics of the extremophile Cryomyces antarcticus and other psychrophilic Dothideomycetes

Over a billion years of fungal evolution has enabled representatives of this kingdom to populate almost all parts of planet Earth and to adapt to some of its most uninhabitable environments including extremes of temperature, salinity, pH, water, light, or other sources of radiation. Cryomyces antarcticus is an endolithic fungus that inhabits rock outcrops in Antarctica. It survives extremes of cold, humidity and solar radiation in one of the least habitable environments on Earth. This fungus is unusual because it produces heavily melanized, meristematic growth and is thought to be haploid and asexual. Due to its growth in the most extreme environment, it has been suggested as an organism that could survive on Mars. However, the mechanisms it uses to achieve its extremophilic nature are not known. Comparative genomics can provide clues to the processes underlying biological diversity, evolution, and adaptation. This effort has been greatly facilitated by the 1000 Fungal Genomes project and the JGI MycoCosm portal where sequenced genomes have been assembled into phylogenetic and ecological groups representing different projects, lifestyles, ecologies, and evolutionary histories. Comparative genomics within and between these groups provides insights into fungal adaptations, for example to extreme environmental conditions. Here, we analyze two Cryomyces genomes in the context of additional psychrophilic fungi, as well as non-psychrophilic fungi with diverse lifestyles selected from the MycoCosm database. This analysis identifies families of genes that are expanded and contracted in Cryomyces and other psychrophiles and may explain their extremophilic lifestyle. Higher GC contents of genes and of bases in the third positions of codons may help to stabilize DNA under extreme conditions. Numerous smaller contigs in C. antarcticus suggest the presence of an alternative haplotype that could indicate the sequenced isolate is diploid or dikaryotic. These analyses provide a first step to unraveling the secrets of the extreme lifestyle of C. antarcticus.

Remodeling of perturbed chromatin can initiate de novo transcriptional and post-transcriptional silencing

In eukaryotes, repetitive DNA can become silenced de novo, either transcriptionally or post-transcriptionally, by processes independent of strong sequence-specific cues. The mechanistic nature of such processes remains poorly understood. We found that in the fungus Neurospora crassa, de novo initiation of both transcriptional and post-transcriptional silencing was linked to perturbed chromatin, which was produced experimentally by the aberrant activity of transcription factors at the tetO operator array. Transcriptional silencing was mediated by canonical constitutive heterochromatin. On the other hand, post-transcriptional silencing resembled repeat-induced quelling but occurred normally when homologous recombination was inactivated. All silencing of the tetO array was dependent on SAD-6, fungal ortholog of the SWI/SNF chromatin remodeler ATRX (Alpha Thalassemia/Mental Retardation Syndrome X-Linked), which was required to maintain nucleosome occupancy at the perturbed locus. In addition, we found that two other types of sequences (the lacO array and native AT-rich DNA) could also undergo recombination-independent quelling associated with perturbed chromatin. These results suggested a model in which the de novo initiation of transcriptional and post-transcriptional silencing is coupled to the remodeling of perturbed chromatin.

Functions and mechanisms of A-to-I RNA editing in filamentous ascomycetes

Although lack of ADAR (adenosine deaminase acting on RNA) orthologs, genome-wide A-to-I editing occurs specifically during sexual reproduction in a number of filamentous ascomycetes, including Fusarium graminearum and Neurospora crassa. Unlike ADAR-mediated editing in animals, fungal A-to-I editing has a strong preference for hairpin loops and U at -1 position, which leads to frequent editing of UAG and UAA stop codons. Majority of RNA editing events in fungi are in the coding region and cause amino acid changes. Some of these editing events have been experimentally characterized for providing heterozygote and adaptive advantages in F. graminearum. Recent studies showed that FgTad2 and FgTad3, 2 ADAT (adenosine deaminase acting on tRNA) enzymes that normally catalyze the editing of A34 in the anticodon of tRNA during vegetative growth mediate A-to-I mRNA editing during sexual reproduction. Stage specificity of RNA editing is conferred by stage-specific expression of short transcript isoforms of FgTAD2 and FgTAD3 as well as cofactors such as AME1 and FIP5 that facilitate the editing of mRNA in perithecia. Taken together, fungal A-to-I RNA editing during sexual reproduction is catalyzed by ADATs and it has the same sequence and structural preferences with editing of A34 in tRNA.

Population-level transposable element expression dynamics influence trait evolution in a fungal crop pathogen

The rapid adaptive evolution of microbes is driven by strong selection pressure acting on genetic variation. How adaptive genetic variation is generated within species and how such variation influences phenotypic trait expression is often not well understood though. We focused on the recent activity of transposable elements (TEs) using deep population genomics and transcriptomics analyses of a fungal plant pathogen with a highly active content of TEs in the genome. Zymoseptoria tritici causes one of the most damaging diseases on wheat, with recent adaptation to the host and environment being facilitated by TE-associated mutations. We obtained genomic and RNA-sequencing data from 146 isolates collected from a single wheat field. We established a genome-wide map of TE insertion polymorphisms in the population by analyzing recent TE insertions among individuals. We quantified the locus-specific transcription of individual TE copies and found considerable population variation at individual TE loci in the population. About 20% of all TE copies show transcription in the genome suggesting that genomic defenses such as repressive epigenetic marks and repeat-induced polymorphisms are at least partially ineffective at preventing the proliferation of TEs in the genome. A quarter of recent TE insertions are associated with expression variation of neighboring genes providing broad potential to influence trait expression. We indeed found that TE insertions are likely responsible for variation in virulence on the host and potentially diverse components of secondary metabolite production. Our large-scale transcriptomics study emphasizes how TE-derived polymorphisms segregate even in individual microbial populations and can broadly underpin trait variation in pathogens.IMPORTANCEPathogens can rapidly adapt to new hosts, antimicrobials, or changes in the environment. Adaptation arises often from mutations in the genome; however, how such variation is generated remains poorly understood. We investigated the most dynamic regions of the genome of Zymoseptoria tritici, a major fungal pathogen of wheat. We focused on the transcription of transposable elements. A large proportion of the transposable elements not only show signatures of potential activity but are also variable within a single population of the pathogen. We find that this variation in activity is likely influencing many important traits of the pathogen. Hence, our work provides insights into how a microbial species can adapt over the shortest time periods based on the activity of transposable elements.

A systematic screen for co-option of transposable elements across the fungal kingdom

How novel protein functions are acquired is a central question in molecular biology. Key paths to novelty include gene duplications, recombination or horizontal acquisition. Transposable elements (TEs) are increasingly recognized as a major source of novel domain-encoding sequences. However, the impact of TE coding sequences on the evolution of the proteome remains understudied. Here, we analyzed 1237 genomes spanning the phylogenetic breadth of the fungal kingdom. We scanned proteomes for evidence of co-occurrence of TE-derived domains along with other conventional protein functional domains. We detected more than 13,000 predicted proteins containing potentially TE-derived domain, of which 825 were identified in more than five genomes, indicating that many host-TE fusions may have persisted over long evolutionary time scales. We used the phylogenetic context to identify the origin and retention of individual TE-derived domains. The most common TE-derived domains are helicases derived from Academ, Kolobok or Helitron. We found putative TE co-options at a higher rate in genomes of the Saccharomycotina, providing an unexpected source of protein novelty in these generally TE depleted genomes. We investigated in detail a candidate host-TE fusion with a heterochromatic transcriptional silencing function that may play a role in TE and gene regulation in ascomycetes. The affected gene underwent multiple full or partial losses within the phylum. Overall, our work establishes a kingdom-wide view of putative host-TE fusions and facilitates systematic investigations of candidate fusion proteins.

Evolutionary dynamics of the LTR-retrotransposon crapaud in the Podospora anserina species complex and the interaction with repeat-induced point mutations

BACKGROUND

The genome of the filamentous ascomycete Podospora anserina shows a relatively high abundance of retrotransposons compared to other interspersed repeats. The LTR-retrotransposon family crapaud is particularly abundant in the genome, and consists of multiple diverged sequence variations specifically localized in the 5' half of both long terminal repeats (LTRs). P. anserina is part of a recently diverged species-complex, which makes the system ideal to classify the crapaud family based on the observed LTR variation and to study the evolutionary dynamics, such as the diversification and bursts of the elements over recent evolutionary time.

RESULTS

We developed a sequence similarity network approach to classify the crapaud repeats of seven genomes representing the P. anserina species complex into 14 subfamilies. This method does not utilize a consensus sequence, but instead it connects any copies that share enough sequence similarity over a set sequence coverage. Based on phylogenetic analyses, we found that the crapaud repeats likely diversified in the ancestor of the complex and have had activity at different time points for different subfamilies. Furthermore, while we hypothesized that the evolution into multiple subfamilies could have been a direct effect of escaping the genome defense system of repeat induced point mutations, we found this not to be the case.

CONCLUSIONS

Our study contributes to the development of methods to classify transposable elements in fungi, and also highlights the intricate patterns of retrotransposon evolution over short timescales and under high mutational load caused by nucleotide-altering genome defense.

Dynamics of transposable element accumulation in the non-recombining regions of mating-type chromosomes in anther-smut fungi

In the absence of recombination, the number of transposable elements (TEs) increases due to less efficient selection, but the dynamics of such TE accumulations are not well characterized. Leveraging a dataset of 21 independent events of recombination cessation of different ages in mating-type chromosomes of Microbotryum fungi, we show that TEs rapidly accumulated in regions lacking recombination, but that TE content reached a plateau at ca. 50% of occupied base pairs by 1.5 million years following recombination suppression. The same TE superfamilies have expanded in independently evolved non-recombining regions, in particular rolling-circle replication elements (Helitrons). Long-terminal repeat (LTR) retrotransposons of the Copia and Ty3 superfamilies also expanded, through transposition bursts (distinguished from gene conversion based on LTR divergence), with both non-recombining regions and autosomes affected, suggesting that non-recombining regions constitute TE reservoirs. This study improves our knowledge of genome evolution by showing that TEs can accumulate through bursts, following non-linear decelerating dynamics.

A DEAD-box RNA helicase mediates meiotic silencing by unpaired DNA

During the sexual phase of Neurospora crassa, unpaired genes are subject to a silencing mechanism known as meiotic silencing by unpaired DNA (MSUD). MSUD targets the transcripts of an unpaired gene and utilizes typical RNA interference factors for its process. Using a reverse genetic screen, we have identified a meiotic silencing gene called sad-9, which encodes a DEAD-box RNA helicase. While not essential for vegetative growth, SAD-9 plays a crucial role in both sexual development and MSUD. Our results suggest that SAD-9, with the help of the SAD-2 scaffold protein, recruits the SMS-2 Argonaute to the perinuclear region, the center of MSUD activity.

Polymorphism of Avirulence Genes and Adaptation to Brassica Resistance Genes Is Gene-Dependent in the Phytopathogenic Fungus Leptosphaeria maculans

The fungal phytopathogen Leptosphaeria maculans, which causes stem canker (blackleg) of rapeseed (Brassica napus), is mainly controlled worldwide by genetic resistance, which includes major resistance genes (Rlm). This model is one of those for which the highest number of avirulence genes (AvrLm) has been cloned. In many systems, including the L. maculans-B. napus interaction, intense use of resistance genes exerts strong selection pressure on the corresponding avirulent isolates, and the fungi may rapidly escape resistance through various molecular events which modify the avirulence genes. In the literature, the study of polymorphism at avirulence loci is often focused on single genes under selection pressure. In this study, we investigate allelic polymorphism at 11 avirulence loci in a French population of 89 L. maculans isolates collected on a trap cultivar in four geographic locations in the 2017-2018 cropping season. The corresponding Rlm genes have been (i) used for a long time, (ii) recently used, or (iii) unused in agricultural practice. The sequence data generated indicate an extreme diversity of situations. For example, genes submitted to an ancient selection may have either been deleted in populations (AvrLm1) or replaced by a single-nucleotide mutated virulent version (AvrLm2, AvrLm5-9). Genes that have never been under selection may either be nearly invariant (AvrLm6, AvrLm10A, AvrLm10B), exhibit rare deletions (AvrLm11, AvrLm14), or display a high diversity of alleles and isoforms (AvrLmS-Lep2). These data suggest that the evolutionary trajectory of avirulence/virulence alleles is gene-dependent and independent of selection pressure in L. maculans. [Formula: see text] Copyright © 2023 The Author(s). This is an open access article distributed under the CC BY-NC-ND 4.0 International license.

Nuclear genome organization in fungi: from gene folding to Rabl chromosomes

Comparative genomics has recently provided unprecedented insights into the biology and evolution of the fungal lineage. In the postgenomics era, a major research interest focuses now on detailing the functions of fungal genomes, i.e. how genomic information manifests into complex phenotypes. Emerging evidence across diverse eukaryotes has revealed that the organization of DNA within the nucleus is critically important. Here, we discuss the current knowledge on the fungal genome organization, from the association of chromosomes within the nucleus to topological structures at individual genes and the genetic factors required for this hierarchical organization. Chromosome conformation capture followed by high-throughput sequencing (Hi-C) has elucidated how fungal genomes are globally organized in Rabl configuration, in which centromere or telomere bundles are associated with opposite faces of the nuclear envelope. Further, fungal genomes are regionally organized into topologically associated domain-like (TAD-like) chromatin structures. We discuss how chromatin organization impacts the proper function of DNA-templated processes across the fungal genome. Nevertheless, this view is limited to a few fungal taxa given the paucity of fungal Hi-C experiments. We advocate for exploring genome organization across diverse fungal lineages to ensure the future understanding of the impact of nuclear organization on fungal genome function.

Current genetic strategies to investigate gene functions in Trichoderma reesei

The filamentous fungus Trichoderma reesei (teleomorph Hypocrea jecorina, Ascomycota) is a well-known lignocellulolytic enzymes-producing strain in industry. To increase the fermentation titer of lignocellulolytic enzymes, random mutagenesis and rational genetic engineering in T. reesei were carried out since it was initially found in the Solomon Islands during the Second World War. Especially the continuous exploration of the underlying regulatory network during (hemi)cellulase gene expression in the post-genome era provided various strategies to develop an efficient fungal cell factory for these enzymes' production. Meanwhile, T. reesei emerges competitiveness potential as a filamentous fungal chassis to produce proteins from other species (e.g., human albumin and interferon α-2b, SARS-CoV-2 N antigen) in virtue of the excellent expression and secretion system acquired during the studies about (hemi)cellulase production. However, all the achievements in high yield of (hemi)cellulases are impossible to finish without high-efficiency genetic strategies to analyze the proper functions of those genes involved in (hemi)cellulase gene expression or secretion. Here, we in detail summarize the current strategies employed to investigate gene functions in T. reesei. These strategies are supposed to be beneficial for extending the potential of T. reesei in prospective strain engineering.

Repeat-Induced Point Mutation and Gene Conversion Coinciding with Heterochromatin Shape the Genome of a Plant-Pathogenic Fungus

Meiosis is associated with genetic changes in the genome-via recombination, gene conversion, and mutations. The occurrence of gene conversion and mutations during meiosis may further be influenced by the chromatin conformation, similar to the effect of the chromatin conformation on the mitotic mutation rate. To date, however, the exact distribution and type of meiosis-associated changes and the role of the chromatin conformation in this context are largely unexplored. Here, we determine recombination, gene conversion, and de novo mutations using whole-genome sequencing of all meiotic products of 23 individual meioses in Zymoseptoria tritici, an important pathogen of wheat. We confirm a high genome-wide recombination rate of 65 centimorgan (cM)/Mb and see higher recombination rates on the accessory compared to core chromosomes. A substantial fraction of 0.16% of all polymorphic markers was affected by gene conversions, showing a weak GC-bias and occurring at higher frequency in regions of constitutive heterochromatin, indicated by the histone modification H3K9me3. The de novo mutation rate associated with meiosis was approximately three orders of magnitude higher than the corresponding mitotic mutation rate. Importantly, repeat-induced point mutation (RIP), a fungal defense mechanism against duplicated sequences, is active in Z. tritici and responsible for the majority of these de novo meiotic mutations. Our results indicate that the genetic changes associated with meiosis are a major source of variability in the genome of an important plant pathogen and shape its evolutionary trajectory. IMPORTANCE The impact of meiosis on the genome composition via gene conversion and mutations is mostly poorly understood, in particular, for non-model species. Here, we sequenced all four meiotic products for 23 individual meioses and determined the genetic changes caused by meiosis for the important fungal wheat pathogen Zymoseptoria tritici. We found a high rate of gene conversions and an effect of the chromatin conformation on gene conversion rates. Higher conversion rates were found in regions enriched with the H3K9me3-a mark for constitutive heterochromatin. Most importantly, meiosis was associated with a much higher frequency of de novo mutations than mitosis; 78% of the meiotic mutations were caused by repeat-induced point mutations-a fungal defense mechanism against duplicated sequences. In conclusion, the genetic changes associated with meiosis are therefore a major factor shaping the genome of this fungal pathogen.

Starships are active eukaryotic transposable elements mobilized by a new family of tyrosine recombinases

Transposable elements in eukaryotic organisms have historically been considered "selfish," at best conferring indirect benefits to their host organisms. The Starships are a recently discovered feature in fungal genomes that are, in some cases, predicted to confer beneficial traits to their hosts and also have hallmarks of being transposable elements. Here, we provide experimental evidence that Starships are indeed autonomous transposons, using the model Paecilomyces variotii, and identify the HhpA "Captain" tyrosine recombinase as essential for their mobilization into genomic sites with a specific target site consensus sequence. Furthermore, we identify multiple recent horizontal gene transfers of Starships, implying that they jump between species. Fungal genomes have mechanisms to defend against mobile elements, which are frequently detrimental to the host. We discover that Starships are also vulnerable to repeat-induced point mutation defense, thereby having implications on the evolutionary stability of such elements.

A population genetics theory for piRNA-regulated transposable elements

Transposable elements (TEs) are self-reproducing selfish DNA sequences that can invade the genome of virtually all living species. Population genetics models have shown that TE copy numbers generally reach a limit, either because the transposition rate decreases with the number of copies (transposition regulation) or because TE copies are deleterious, and thus purged by natural selection. Yet, recent empirical discoveries suggest that TE regulation may mostly rely on piRNAs, which require a specific mutational event (the insertion of a TE copy in a piRNA cluster) to be activated - the so-called TE regulation "trap model". We derived new population genetics models accounting for this trap mechanism, and showed that the resulting equilibria differ substantially from previous expectations based on a transposition-selection equilibrium. We proposed three sub-models, depending on whether or not genomic TE copies and piRNA cluster TE copies are selectively neutral or deleterious, and we provide analytical expressions for maximum and equilibrium copy numbers, as well as cluster frequencies for all of them. In the full neutral model, the equilibrium is achieved when transposition is completely silenced, and this equilibrium does not depend on the transposition rate. When genomic TE copies are deleterious but not cluster TE copies, no long-term equilibrium is possible, and active TEs are eventually eliminated after an active incomplete invasion stage. When all TE copies are deleterious, a transposition-selection equilibrium exists, but the invasion dynamics is not monotonic, and the copy number peaks before decreasing. Mathematical predictions were in good agreement with numerical simulations, except when genetic drift and/or linkage disequilibrium dominates. Overall, the trap-model dynamics appeared to be substantially more stochastic and less repeatable than traditional regulation models.

The model diatom Phaeodactylum tricornutum provides insights into the diversity and function of microeukaryotic DNA methyltransferases

Cytosine methylation is an important epigenetic mark involved in the transcriptional control of transposable elements in mammals, plants and fungi. The Stramenopiles-Alveolate-Rhizaria (SAR) lineages are a major group of ecologically important marine microeukaryotes, including the phytoplankton groups diatoms and dinoflagellates. However, little is known about their DNA methyltransferase diversity. Here, we performed an in-silico analysis of DNA methyltransferases found in marine microeukaryotes and showed that they encode divergent DNMT3, DNMT4, DNMT5 and DNMT6 enzymes. Furthermore, we found three classes of enzymes within the DNMT5 family. Using a CRISPR/Cas9 strategy we demonstrated that the loss of the DNMT5a gene correlates with a global depletion of DNA methylation and overexpression of young transposable elements in the model diatom Phaeodactylum tricornutum. The study provides a view of the structure and function of a DNMT family in the SAR supergroup using an attractive model species.

A thousand-genome panel retraces the global spread and adaptation of a major fungal crop pathogen

Human activity impacts the evolutionary trajectories of many species worldwide. Global trade of agricultural goods contributes to the dispersal of pathogens reshaping their genetic makeup and providing opportunities for virulence gains. Understanding how pathogens surmount control strategies and cope with new climates is crucial to predicting the future impact of crop pathogens. Here, we address this by assembling a global thousand-genome panel of Zymoseptoria tritici, a major fungal pathogen of wheat reported in all production areas worldwide. We identify the global invasion routes and ongoing genetic exchange of the pathogen among wheat-growing regions. We find that the global expansion was accompanied by increased activity of transposable elements and weakened genomic defenses. Finally, we find significant standing variation for adaptation to new climates encountered during the global spread. Our work shows how large population genomic panels enable deep insights into the evolutionary trajectory of a major crop pathogen.

Investigating the origin of subtelomeric and centromeric AT-rich elements in Aspergillus flavus

An in silico study of Aspergillus flavus genome stability uncovered significant variations in both coding and non-coding regions. The non-coding insertions uniformly consisted of AT-rich sequences that are evolutionarily maintained, albeit distributed at widely different sites in an array of A. flavus strains. A survey of ≥ 2kb AT-rich elements (AT ≥ 70%; ATEs) in non-centromeric regions uncovered two major categories of ATEs. The first category is composed of homologous insertions at ectopic, non-allelic sites that contain homology to transposable elements (TEs; Classes B, C, D, and E). Strains differed significantly in frequency, position, and TE type, but displayed a common enrichment in subtelomeric regions. The TEs were heavily mutated, with patterns consistent with the ancestral activity of repeat-induced point mutations (RIP). The second category consists of a conserved set of novel subtelomeric ATE repeats (Classes A, G, G, H, I and J) which lack discernible TEs and, unlike TEs, display a constant polarity relative to the telomere. Members of one of these classes are derivatives of a progenitor ATE that is predicted to have undergone extensive homologous recombination during evolution. A third category of ATEs consists of ~100 kb regions at each centromere. Centromeric ATEs and TE clusters within these centromeres display a high level of sequence identity between strains. These studies suggest that transposition and RIP are forces in the evolution of subtelomeric and centromeric structure and function.

Two-speed genome evolution drives pathogenicity in fungal pathogens of animals

The origins and evolution of virulence in amphibian-infecting chytrids Batrachochytrium dendrobatidis (Bd) and Batrachochytrium salamandrivorans (Bsal) are largely unknown. Here, we use deep nanopore sequencing of Bsal and comparative genomics against 21 high-quality genome assemblies that span the fungal Chytridiomycota. We discover that Bsal has the most repeat-rich genome of the Chytridiomycota, comprising 40.9% repetitive elements; this genome has expanded to more than 3× the length of its conspecific Bd, with autonomous and fully functional LTR/Gypsy elements contributing significantly to the expansion. The M36 metalloprotease virulence factors are highly expanded (n = 177) in Bsal, most of which (53%) are flanked by transposable elements, suggesting they have a repeat-associated expansion. We find enrichment upstream of M36 metalloprotease genes of three novel repeat families belonging to the repeat superfamily of LINEs that are implicated with gene copy number variations. Additionally, Bsal has a highly compartmentalized genome architecture, with virulence factors enriched in gene-sparse/repeat-rich compartments, while core conserved genes are enriched in gene-rich/repeat-poor compartments. Genes upregulated during infection are primarily found in the gene-sparse/repeat-rich compartment in both Bd and Bsal. Furthermore, genes with signatures of positive selection in Bd are enriched in repeat-rich regions, suggesting these regions are a cradle for the evolution of chytrid pathogenicity. These are the hallmarks of two-speed genome evolution, and this study provides evidence of two-speed genomes in an animal pathogen, shedding light on the evolution of fungal pathogens of vertebrates driving global declines and extinctions.

Heterochromatin in the fungal plant pathogen, Zymoseptoria tritici: Control of transposable elements, genome plasticity and virulence

Heterochromatin is a repressive chromatin state that plays key roles in the functional organisation of eukaryotic genomes. In fungal plant pathogens, effector genes that are required for host colonization tend to be associated with heterochromatic regions of the genome that are enriched with transposable elements. It has been proposed that the heterochromatin environment silences effector genes in the absence of host and dynamic chromatin remodelling facilitates their expression during infection. Here we discuss this model in the context of the key wheat pathogen, Zymoseptoria tritici. We cover progress in understanding the deposition and recognition of heterochromatic histone post translational modifications in Z. tritici and the role that heterochromatin plays in control of genome plasticity and virulence.

Encyclopaedia of eukaryotic DNA methylation: from patterns to mechanisms and functions

DNA methylation is an epigenetic modification with a very long evolutionary history. However, DNA methylation evolves surprisingly rapidly across eukaryotes. The genome-wide distribution of methylation diversifies rapidly in different lineages, and DNA methylation is lost altogether surprisingly frequently. The growing availability of genomic and epigenomic sequencing across organisms highlights this diversity but also illuminates potential factors that could explain why both the DNA methylation machinery and its genome-wide distribution evolve so rapidly. Key to this are new discoveries about the fitness costs associated with DNA methylation, and new theories about how the fundamental biochemical mechanisms of DNA methylation introduction and maintenance could explain how new genome-wide patterns of methylation evolve.

Telomere-to-Telomere Genome Sequences across a Single Genus Reveal Highly Variable Chromosome Rearrangement Rates but Absolute Stasis of Chromosome Number

Genome rearrangements in filamentous fungi are prevalent but little is known about the modalities of their evolution, in part because few complete genomes are available within a single genus. To address this, we have generated and compared 15 complete telomere-to-telomere genomes across the phylogeny of a single genus of filamentous fungi, Epichloë. We find that the striking distinction between gene-rich and repeat-rich regions previously reported for isolated species is ubiquitous across the Epichloë genus. We built a species phylogeny from single-copy gene orthologs to provide a comparative framing to study chromosome composition and structural change through evolutionary time. All Epichloë genomes have exactly seven nuclear chromosomes, but despite this conserved ploidy, analyses reveal low synteny and substantial rearrangement of gene content across the genus. These rearrangements are highly lineage-dependent, with most occurring over short evolutionary distances, with long periods of structural stasis. Quantification of chromosomal rearrangements shows they are uncorrelated with numbers of substitutions and evolutionary distances, suggesting that different modes of evolution are acting to create nucleotide and chromosome-scale changes.

Intra-Species Genomic Variation in the Pine Pathogen Fusarium circinatum

Fusarium circinatum is an important global pathogen of pine trees. Genome plasticity has been observed in different isolates of the fungus, but no genome comparisons are available. To address this gap, we sequenced and assembled to chromosome level five isolates of F. circinatum. These genomes were analysed together with previously published genomes of F. circinatum isolates, FSP34 and KS17. Multi-sample variant calling identified a total of 461,683 micro variants (SNPs and small indels) and a total of 1828 macro structural variants of which 1717 were copy number variants and 111 were inversions. The variant density was higher on the sub-telomeric regions of chromosomes. Variant annotation revealed that genes involved in transcription, transport, metabolism and transmembrane proteins were overrepresented in gene sets that were affected by high impact variants. A core genome representing genomic elements that were conserved in all the isolates and a non-redundant pangenome representing all genomic elements is presented. Whole genome alignments showed that an average of 93% of the genomic elements were present in all isolates. The results of this study reveal that some genomic elements are not conserved within the isolates and some variants are high impact. The described genome-scale variations will help to inform novel disease management strategies against the pathogen.

Tracking of Diversity and Evolution in the Brown Rot Fungi Monilinia fructicola, Monilinia fructigena, and Monilinia laxa

Monilinia species are among the most devastating fungi worldwide as they cause brown rot and blossom blight on fruit trees. To understand the molecular bases of their pathogenic lifestyles, we compared the newly assembled genomes of single strains of Monilinia fructicola, M. fructigena and M. laxa, with those of Botrytis cinerea and Sclerotinia sclerotiorum, as the closest species within Sclerotiniaceae. Phylogenomic analysis of orthologous proteins and syntenic investigation suggest that M. laxa is closer to M. fructigena than M. fructicola, and is closest to the other investigated Sclerotiniaceae species. This indicates that M. laxa was the earliest result of the speciation process. Distinct evolutionary profiles were observed for transposable elements (TEs). M. fructicola and M. laxa showed older bursts of TE insertions, which were affected (mainly in M. fructicola) by repeat-induced point (RIP) mutation gene silencing mechanisms. These suggested frequent occurrence of the sexual process in M. fructicola. More recent TE expansion linked with low RIP action was observed in M. fructigena, with very little in S. sclerotiorum and B. cinerea. The detection of active non-syntenic TEs is indicative of horizontal gene transfer and has resulted in alterations in specific gene functions. Analysis of candidate effectors, biosynthetic gene clusters for secondary metabolites and carbohydrate-active enzymes, indicated that Monilinia genus has multiple virulence mechanisms to infect host plants, including toxins, cell-death elicitor, putative virulence factors and cell-wall-degrading enzymes. Some species-specific pathogenic factors might explain differences in terms of host plant and organ preferences between M. fructigena and the other two Monilinia species.

Genetic Characteristics and Metabolic Interactions between Pseudocercospora fijiensis and Banana: Progress toward Controlling Black Sigatoka

The international importance of banana and severity of black Sigatoka disease have led to extensive investigations into the genetic characteristics and metabolic interactions between the Dothideomycete Pseudocercospora fijiensis and its banana host. P. fijiensis was shown to have a greatly expanded genome compared to other Dothideomycetes, due to the proliferation of retrotransposons. Genome analysis suggests the presence of dispensable chromosomes that may aid in fungal adaptation as well as pathogenicity. Genomic research has led to the characterization of genes and metabolic pathways involved in pathogenicity, including: secondary metabolism genes such as PKS10-2, genes for mitogen-activated protein kinases such as Fus3 and Slt2, and genes for cell wall proteins such as glucosyl phosphatidylinositol (GPI) and glycophospholipid surface (Gas) proteins. Studies conducted on resistance mechanisms in banana have documented the role of jasmonic acid and ethylene pathways. With the development of banana transformation protocols, strategies for engineering resistance include transgenes expressing antimicrobial peptides or hydrolytic enzymes as well as host-induced gene silencing (HIGS) targeting pathogenicity genes. Pseudocercospora fijiensis has been identified as having high evolutionary potential, given its large genome size, ability to reproduce both sexually and asexually, and long-distance spore dispersal. Thus, multiple control measures are needed for the sustainable control of black Sigatoka disease.

Speciation of Genes and Genomes: Conservation of DNA Polymorphism by Barriers to Recombination Raised by Mismatch Repair System

Some basic aspects of human and animal biology and evolution involve the establishment of biological uniqueness of species and individuals within their huge variety. The discrimination among closely related species occurs in their offspring at the level of chromosomal DNA sequence homology, which is required for fertility as the hallmark of species. Biological identification of individuals, i.e., of their biological “self”, occurs at the level of protein sequences presented by the MHC/HLA complex as part of the immune system that discriminates non-self from self. Here, a mechanistic molecular model is presented that can explain how DNA sequence divergence and the activity of key mismatch repair proteins, MutS and MutL, lead to 1) genetic separation of closely related species (sympatric speciation) (Fitch and Ayala, Proceedings of the National Academy of Sciences, 1994, 91, 6717–6720), 2) the stability of genomes riddled by diverged repeated sequences, and 3) conservation of highly polymorphic DNA sequence blocks that constitute the immunological self. All three phenomena involve suppression of recombination between diverged homologies, resulting in prevention of gene sharing between closely related genomes (evolution of new species) as well as sequence sharing between closely related genes within a genome (e.g., evolution of immunoglobulin, MHC, and other gene families bearing conserved polymorphisms).

A-to-I mRNA editing controls spore death induced by a fungal meiotic drive gene in homologous and heterologous expression systems

Spore killers are meiotic drive elements that can block development of sexual spores in fungi. In the maize ear rot and mycotoxin-producing fungus Fusarium verticillioides, a spore killer called SkK has been mapped to a 102-kb interval of chromosome V. Here, we show that a gene within this interval, SKC1, is required for SkK-mediated spore killing and meiotic drive. We also demonstrate that SKC1 is associated with at least four transcripts, two sense (sense-SKC1a and sense-SKC1b) and two antisense (antisense-SKC1a and antisense-SKC1b). Both antisense SKC1 transcripts lack obvious protein-coding sequences and thus appear to be non-coding RNAs. In contrast, sense-SKC1a is a protein-coding transcript that undergoes A-to-I editing to sense-SKC1b in sexual tissue. Translation of sense-SKC1a produces a 70 amino acid protein (Skc1a), whereas translation of sense-SKC1b produces an 84 amino acid protein (Skc1b). Heterologous expression analysis of SKC1 transcripts shows that sense-SKC1a also undergoes A-to-I editing to sense-SKC1b during the Neurospora crassa sexual cycle. Site directed mutagenesis studies indicate that Skc1b is responsible for spore killing in F. verticillioides and that it induces most meiotic cells to die in N. crassa. Finally, we report that SKC1 homologs are present in over 20 Fusarium species. Overall, our results demonstrate that fungal meiotic drive elements like SKC1 can influence the outcome of meiosis by hijacking a cell's A-to-I editing machinery and that the involvement of A-to-I editing in a fungal meiotic drive system does not preclude its horizontal transfer to a distantly related species.

Chromatin profiling reveals heterogeneity in clinical isolates of the human pathogen Aspergillus fumigatus

Invasive Pulmonary Aspergillosis, which is caused by the filamentous fungus Aspergillus fumigatus, is a life-threatening infection for immunosuppressed patients. Chromatin structure regulation is important for genome stability maintenance and has the potential to drive genome rearrangements and affect virulence and pathogenesis of pathogens. Here, we performed the first A. fumigatus global chromatin profiling of two histone modifications, H3K4me3 and H3K9me3, focusing on the two most investigated A. fumigatus clinical isolates, Af293 and CEA17. In eukaryotes, H3K4me3 is associated with active transcription, while H3K9me3 often marks silent genes, DNA repeats, and transposons. We found that H3K4me3 deposition is similar between the two isolates, while H3K9me3 is more variable and does not always represent transcriptional silencing. Our work uncovered striking differences in the number, locations, and expression of transposable elements between Af293 and CEA17, and the differences are correlated with H3K9me3 modifications and higher genomic variations among strains of Af293 background. Moreover, we further showed that the Af293 strains from different laboratories actually differ in their genome contents and found a frequently lost region in chromosome VIII. For one such Af293 variant, we identified the chromosomal changes and demonstrated their impacts on its secondary metabolites production, growth and virulence. Overall, our findings not only emphasize the influence of genome heterogeneity on A. fumigatus fitness, but also caution about unnoticed chromosomal variations among common laboratory strains.

Complete Genome Sequences and Genome-Wide Characterization of Trichoderma Biocontrol Agents Provide New Insights into their Evolution and Variation in Genome Organization, Sexual Development, and Fungal-Plant Interactions

Trichoderma spp. represent one of the most important fungal genera to mankind and in natural environments. The genus harbors prolific producers of wood-decaying enzymes, biocontrol agents against plant pathogens, plant-growth-promoting biofertilizers, as well as model organisms for studying fungal-plant-plant pathogen interactions. Pursuing highly accurate, contiguous, and chromosome-level reference genomes has become a primary goal of fungal research communities. Here, we report the chromosome-level genomic sequences and whole-genome annotation data sets of four strains used as biocontrol agents or biofertilizers (Trichoderma virens Gv29-8, Trichoderma virens FT-333, Trichoderma asperellum FT-101, and Trichoderma atroviride P1). Our results provide comprehensive categorization, correct positioning, and evolutionary detail of both nuclear and mitochondrial genomes, including telomeres, AT-rich blocks, centromeres, transposons, mating-type loci, nuclear-encoded mitochondrial sequences, as well as many new secondary metabolic and carbohydrate-active enzyme gene clusters. We have also identified evolutionarily conserved core genes contributing to plant-fungal interactions, as well as variations potentially linked to key behavioral traits such as sex, genome defense, secondary metabolism, and mycoparasitism. The genomic resources we provide herein significantly extend our knowledge not only of this economically important fungal genus, but also fungal evolution and basic biology in general. IMPORTANCE Telomere-to-telomere and gapless reference genome assemblies are necessary to ensure that all genomic variants are studied and discovered, including centromeres, telomeres, AT-rich blocks, mating type loci, biosynthetic, and metabolic gene clusters. Here, we applied long-range sequencing technologies to determine the near-completed genome sequences of four widely used biocontrol agents or biofertilizers: Trichoderma virens Gv29-8 and FT-333, Trichoderma asperellum FT-101, and Trichoderma atroviride P1. Like those of three Trichoderma reesei wild isolates [QM6a, CBS999.97(MAT1-1) and CBS999.97(MAT1-2)] we reported previously, these four biocontrol agent genomes each contain seven nuclear chromosomes and a circular mitochondrial genome. Substantial intraspecies and intragenus diversities are also discovered, including single nucleotide polymorphisms, chromosome shuffling, as well as genomic relics derived from historical transposition events and repeat-induced point (RIP) mutations.

Transcriptional activity and epigenetic regulation of transposable elements in the symbiotic fungus Rhizophagus irregularis

Arbuscular mycorrhizal (AM) fungi form mutualistic relationships with most land plant species. AM fungi have long been considered as ancient asexuals. Long-term clonal evolution would be remarkable for a eukaryotic lineage and suggests the importance of alternative mechanisms to promote genetic variability facilitating adaptation. Here, we assessed the potential of transposable elements for generating such genomic diversity. The dynamic expression of TEs during Rhizophagus irregularis spore development suggests ongoing TE activity. We find Mutator-like elements located near genes belonging to highly expanded gene families. Whole-genome epigenomic profiling of R. irregularis provides direct evidence of DNA methylation and small RNA production occurring at TE loci. Our results support a model in which TE activity shapes the genome, while DNA methylation and small RNA-mediated silencing keep their overproliferation in check. We propose that a well-controlled TE activity directly contributes to genome evolution in AM fungi.

Modulation of C-to-T mutation by recombination-independent pairing of closely positioned DNA repeats

Repeat-induced point mutation is a genetic process that creates cytosine-to-thymine (C-to-T) transitions in duplicated genomic sequences in fungi. Repeat-induced point mutation detects duplications (irrespective of their origin, specific sequence, coding capacity, and genomic positions) by a recombination-independent mechanism that likely matches intact DNA double helices directly, without relying on the annealing of complementary single strands. In the fungus Neurospora crassa, closely positioned repeats can induce mutation of the adjoining nonrepetitive regions. This process is related to heterochromatin assembly and requires the cytosine methyltransferase DIM-2. Using DIM-2-dependent mutation as a readout of homologous pairing, we find that GC-rich repeats produce a much stronger response than AT-rich repeats, independently of their intrinsic propensity to become mutated. We also report that direct repeats trigger much stronger DIM-2-dependent mutation than inverted repeats. These results can be rationalized in the light of a recently proposed model of homologous DNA pairing, in which DNA double helices associate by forming sequence-specific quadruplex-based contacts with a concomitant release of supercoiling. A similar process featuring pairing-induced supercoiling may initiate epigenetic silencing of repetitive DNA in other organisms, including humans.

Involvement of RNA granule proteins in meiotic silencing by unpaired DNA

In Neurospora crassa, expression from an unpaired gene is suppressed by a mechanism known as meiotic silencing by unpaired DNA (MSUD). MSUD utilizes common RNA interference (RNAi) factors to silence target mRNAs. Here, we report that Neurospora CAR-1 and CGH-1, homologs of two Caenorhabditis elegans RNA granule components, are involved in MSUD. These fungal proteins are found in the perinuclear region and P-bodies, much like their worm counterparts. They interact with components of the meiotic silencing complex (MSC), including the SMS-2 Argonaute. This is the first time MSUD has been linked to RNA granule proteins.

A devil's bargain with transposable elements in plant pathogens

Transposable elements (TEs) spread in genomes through self-copying mechanisms and are a major cause of genome expansions. Plant pathogens have finely tuned the expression of virulence factors to rely on epigenetic control targeted at nearby TEs. Stress experienced during the plant infection process leads to derepression of TEs and concurrently allows the expression of virulence factors. We argue that the derepression of TEs elements causes an evolutionary conflict by favoring TEs that can be reactivated. Active TEs and recent genome size expansions indicate that plant pathogens could face long-term consequences from the short-term benefit of fine-tuning the infection process. Hence, encoding key virulence factors close to TEs under epigenetic control constitutes a devil's bargain for pathogens.

Retrotransposons as pathogenicity factors of the plant pathogenic fungus Botrytis cinerea

BACKGROUND

Retrotransposons are genetic elements inducing mutations in all domains of life. Despite their detrimental effect, retrotransposons can become temporarily active during epigenetic reprogramming and cellular stress response, which may accelerate host genome evolution. In fungal pathogens, a positive role has been attributed to retrotransposons when shaping genome architecture and expression of genes encoding pathogenicity factors; thus, retrotransposons are known to influence pathogenicity.

RESULTS

We uncover a hitherto unknown role of fungal retrotransposons as being pathogenicity factors, themselves. The aggressive fungal plant pathogen, Botrytis cinerea, is known to deliver some long-terminal repeat (LTR) deriving regulatory trans-species small RNAs (BcsRNAs) into plant cells to suppress host gene expression for infection. We find that naturally occurring, less aggressive B. cinerea strains possess considerably lower copy numbers of LTR retrotransposons and had lost retrotransposon BcsRNA production. Using a transgenic proof-of-concept approach, we reconstitute retrotransposon expression in a BcsRNA-lacking B. cinerea strain, which results in enhanced aggressiveness in a retrotransposon and BcsRNA expression-dependent manner. Moreover, retrotransposon expression in B. cinerea leads to suppression of plant defence-related genes during infection.

CONCLUSIONS

We propose that retrotransposons are pathogenicity factors that manipulate host plant gene expression by encoding trans-species BcsRNAs. Taken together, the novelty that retrotransposons are pathogenicity factors will have a broad impact on studies of host-microbe interactions and pathology.

Recombination-independent recognition of DNA homology for meiotic silencing in Neurospora crassa

The pairing of homologous chromosomes represents a critical step of meiosis in nearly all sexually reproducing species. In many organisms, pairing involves chromosomes that remain apparently intact. The mechanistic nature of homology recognition at the basis of such pairing is unknown. Using "meiotic silencing by unpaired DNA" (MSUD) as a model process, we demonstrate the existence of a cardinally different approach to DNA homology recognition in meiosis. The main advantage of MSUD over other experimental systems lies in its ability to identify any relatively short DNA fragment lacking a homologous allelic partner. Here, we show that MSUD does not rely on the canonical mechanism of meiotic recombination, yet it is promoted by REC8, a conserved component of the meiotic cohesion complex. We also show that certain patterns of interspersed homology are recognized as pairable during MSUD. Such patterns need to be colinear and must contain short tracts of sequence identity spaced apart at 21 or 22 base pairs. By using these periodicity values as a guiding parameter in all-atom molecular modeling, we discover that homologous DNA molecules can pair by forming quadruplex-based contacts with an interval of 2.5 helical turns. This process requires right-handed plectonemic coiling and additional conformational changes in the intervening double-helical segments. Our results 1) reconcile genetic and biophysical evidence for the existence of direct homologous double-stranded DNA (dsDNA)-dsDNA pairing, 2) identify a role for this process in initiating RNA interference, and 3) suggest that chromosomes can be cross-matched by a precise mechanism that operates on intact dsDNA molecules.

Population genomics of transposable element activation in the highly repressive genome of an agricultural pathogen

The activity of transposable elements (TEs) can be an important driver of genetic diversity with TE-mediated mutations having a wide range of fitness consequences. To avoid deleterious effects of TE activity, some fungi have evolved highly sophisticated genomic defences to reduce TE proliferation across the genome. Repeat-induced point mutation (RIP) is a fungal-specific TE defence mechanism efficiently targeting duplicated sequences. The rapid accumulation of RIPs is expected to deactivate TEs over the course of a few generations. The evolutionary dynamics of TEs at the population level in a species with highly repressive genome defences is poorly understood. Here, we analyse 366 whole-genome sequences of Parastagonospora nodorum, a fungal pathogen of wheat with efficient RIP. A global population genomics analysis revealed high levels of genetic diversity and signs of frequent sexual recombination. Contrary to expectations for a species with RIP, we identified recent TE activity in multiple populations. The TE composition and copy numbers showed little divergence among global populations regardless of the demographic history. Miniature inverted-repeat transposable elements (MITEs) and terminal repeat retrotransposons in miniature (TRIMs) were largely underlying recent intra-species TE expansions. We inferred RIP footprints in individual TE families and found that recently active, high-copy TEs have possibly evaded genomic defences. We find no evidence that recent positive selection acted on TE-mediated mutations rather that purifying selection maintained new TE insertions at low insertion frequencies in populations. Our findings highlight the complex evolutionary equilibria established by the joint action of TE activity, selection and genomic repression.

Loss of EZH2-like or SU(VAR)3-9-like proteins causes simultaneous perturbations in H3K27 and H3K9 tri-methylation and associated developmental defects in the fungus Podospora anserina

BACKGROUND

Selective gene silencing is key to development. It is generally accepted that H3K27me3-enriched heterochromatin maintains transcriptional repression established during early development and regulates cell fate. Conversely, H3K9me3-enriched heterochromatin prevents differentiation but constitutes protection against transposable elements. We exploited the fungus Podospora anserina, a valuable alternative to higher eukaryote models, to question the biological relevance and functional interplay of these two distinct heterochromatin conformations.

RESULTS

We established genome-wide patterns of H3K27me3 and H3K9me3 modifications, and found these marks mutually exclusive within gene-rich regions but not within repeats. We generated the corresponding histone methyltransferase null mutants and showed an interdependence of H3K9me3 and H3K27me3 marks. Indeed, removal of the PaKmt6 EZH2-like enzyme resulted not only in loss of H3K27me3 but also in significant H3K9me3 reduction. Similarly, removal of PaKmt1 SU(VAR)3-9-like enzyme caused loss of H3K9me3 and substantial decrease of H3K27me3. Removal of the H3K9me binding protein PaHP1 provided further support to the notion that each type of heterochromatin requires the presence of the other. We also established that P. anserina developmental programs require H3K27me3-mediated silencing, since loss of the PaKmt6 EZH2-like enzyme caused severe defects in most aspects of the life cycle including growth, differentiation processes and sexual reproduction, whereas loss of the PaKmt1 SU(VAR)3-9-like enzyme resulted only in marginal defects, similar to loss of PaHP1.

CONCLUSIONS

Our findings support a conserved function of the PRC2 complex in fungal development. However, we uncovered an intriguing evolutionary fluidity in the repressive histone deposition machinery, which challenges canonical definitions of constitutive and facultative heterochromatin.

Austropuccinia psidii, causing myrtle rust, has a gigabase-sized genome shaped by transposable elements

Austropuccinia psidii, originating in South America, is a globally invasive fungal plant pathogen that causes rust disease on Myrtaceae. Several biotypes are recognized, with the most widely distributed pandemic biotype spreading throughout the Asia-Pacific and Oceania regions over the last decade. Austropuccinia psidii has a broad host range with more than 480 myrtaceous species. Since first detected in Australia in 2010, the pathogen has caused the near extinction of at least three species and negatively affected commercial production of several Myrtaceae. To enable molecular and evolutionary studies into A. psidii pathogenicity, we assembled a highly contiguous genome for the pandemic biotype. With an estimated haploid genome size of just over 1 Gb (gigabases), it is the largest assembled fungal genome to date. The genome has undergone massive expansion via distinct transposable element (TE) bursts. Over 90% of the genome is covered by TEs predominantly belonging to the Gypsy superfamily. These TE bursts have likely been followed by deamination events of methylated cytosines to silence the repetitive elements. This in turn led to the depletion of CpG sites in TEs and a very low overall GC content of 33.8%. Compared to other Pucciniales, the intergenic distances are increased by an order of magnitude indicating a general insertion of TEs between genes. Overall, we show how TEs shaped the genome evolution of A. psidii and provide a greatly needed resource for strategic approaches to combat disease spread.

TSETA: A Third-Generation Sequencing-Based Computational Tool for Mapping and Visualization of SNPs, Meiotic Recombination Products, and RIP Mutations

TSETA (Third-generation Sequencing to Enable Tetrad Analysis) is a fungus-centric software pipeline that utilizes chromosome-level sequence assembly for genome-wide and single-nucleotide-resolution mapping of single-nucleotide polymorphisms (SNPs), meiotic recombination products, illegitimate mutations (IMs) and repeat-induced point (RIP) mutations. It utilizes a newly invented algorithm (i.e., BLASTN-guided sectional MAFFT) to perform fast, accurate, and low-cost multiple genome sequence alignments. This new algorithm outcompetes next-generation sequencing (NGS)-based variant-calling approaches for accurate and comprehensive identification of single-nucleotide variants (SNVs) and insertion/deletion mutations (Indels) among the near-complete genome sequences of any two or more intraspecific strains, as well as sequences before and after meiosis, with single-nucleotide precision. TSETA also has a powerful tool for the visualization of the results from the scale of the chromosomal landscape to individual nucleotides. The data output files are user-friendly for researchers and students lacking computational expertise to analyze and reason about data and evidence.

The challenges of predicting transposable element activity in hybrids

Transposable elements (TEs) are ubiquitous mobile genetic elements that hold both disruptive and adaptive potential for species. It has long been postulated that their activity may be triggered by hybridization, a hypothesis that received mixed support from studies in various species. While host defense mechanisms against TEs are being elucidated, the increasing volume of genomic data and bioinformatic tools specialized in TE detection enable in-depth characterization of TEs at the levels of species and populations. Here, I borrow elements from the genome ecology theory to illustrate how knowledge of the diversity of TEs and host defense mechanisms may help predict the activity of TEs in the face of hybridization, and how current limitations make this task especially challenging.

Trichoderma reesei Rad51 tolerates mismatches in hybrid meiosis with diverse genome sequences

Sexual eukaryotes fall into two groups with respect to their RecA-like recombinases. The first group possesses Rad51 (ubiquitous) and Dmc1 (meiosis-specific), which cooperate to conduct interhomolog recombination in zygotes with high sequence heterogeneity. Interestingly, Dmc1 was lost from the second group of eukaryotic organisms. Here we used the industrial workhorse fungus Trichoderma reesei to address if and how Rad51-only eukaryotes carry out hybrid meiosis. We show that T. reesei Rad51 (TrRad51) is indispensable for interhomolog recombination during meiosis and that TrRad51, like Saccharomyces cerevisiae Dmc1, possesses a better mismatch tolerability than S. cerevisiae Rad51. Our results indicate that the ancestral TrRad51 evolved to acquire Dmc1-like properties by adopting multiple structural variations in the L1 and L2 DNA-binding loops.

Most eukaryotes possess two RecA-like recombinases (ubiquitous Rad51 and meiosis-specific Dmc1) to promote interhomolog recombination during meiosis. However, some eukaryotes have lost Dmc1. Given that mammalian and yeast Saccharomyces cerevisiae (Sc) Dmc1 have been shown to stabilize recombination intermediates containing mismatches better than Rad51, we used the Pezizomycotina filamentous fungus Trichoderma reesei to address if and how Rad51-only eukaryotes conduct interhomolog recombination in zygotes with high sequence heterogeneity. We applied multidisciplinary approaches (next- and third-generation sequencing technology, genetics, cytology, bioinformatics, biochemistry, and single-molecule biophysics) to show that T. reesei Rad51 (TrRad51) is indispensable for interhomolog recombination during meiosis and, like ScDmc1, TrRad51 possesses better mismatch tolerance than ScRad51 during homologous recombination. Our results also indicate that the ancestral TrRad51 evolved to acquire ScDmc1-like properties by creating multiple structural variations, including via amino acid residues in the L1 and L2 DNA-binding loops.

Diversity of Fungal DNA Methyltransferases and Their Association With DNA Methylation Patterns

DNA methyltransferases (DNMTs) are a group of proteins that catalyze DNA methylation by transferring a methyl group to DNA. The genetic variation in DNMTs results in differential DNA methylation patterns associated with various biological processes. In fungal species, DNMTs and their DNA methylation profiles were found to be very diverse and have gained many research interests. We reviewed fungal DNMTs in terms of their biological functions, protein domain structures, and their associated epigenetic regulations compared to those known in plant and animal systems. In addition, we summarized recent reports on potential RNA-directed DNA methylation (RdDM) related to DNMT5 in fungi. We surveyed up to 40 fungal species with published genome-wide DNA methylation profiles (methylomes) and presented the associations between the specific patterns of fungal DNA methylation and their DNMTs based on a phylogenetic tree of protein domain structures. For example, the main DNMTs in Basidiomycota, DNMT1 with RFD domain + DNMT5, contributing to CG methylation preference, were distinct from RID + Dim-2 in Ascomycota, resulting in a non-CG methylation preference. Lastly, we revealed that the dynamic methylation involved in fungal life stage changes was particularly low in mycelium and DNA methylation was preferentially located in transposable elements (TEs). This review comprehensively discussed fungal DNMTs and methylomes and their connection with fungal development and taxonomy to present the diverse usages of DNA methylation in fungal genomes.

Population Genomics of Fungal Plant Pathogens and the Analyses of Rapidly Evolving Genome Compartments

Genome sequencing of fungal pathogens have documented extensive variation in genome structure and composition between species and in many cases between individuals of the same species. This type of genomic variation can be adaptive for pathogens to rapidly evolve new virulence phenotypes. Analyses of genome-wide variation in fungal pathogen genomes rely on high quality assemblies and methods to detect and quantify structural variation. Population genomic studies in fungi have addressed the underlying mechanisms whereby structural variation can be rapidly generated. Transposable elements, high mutation and recombination rates as well as incorrect chromosome segregation during mitosis and meiosis contribute to extensive variation observed in many species. We here summarize key findings in the field of fungal pathogen genomics and we discuss methods to detect and characterize structural variants including an alignment-based pipeline to study variation in population genomic data.

Metarhizium: jack of all trades, master of many

The genus Metarhizium and Pochonia chlamydosporia comprise a monophyletic clade of highly abundant globally distributed fungi that can transition between long-term beneficial associations with plants to transitory pathogenic associations with frequently encountered protozoans, nematodes or insects. Some very common 'specialist generalist' species are adapted to particular soil and plant ecologies, but can overpower a wide spectrum of insects with numerous enzymes and toxins that result from extensive gene duplications made possible by loss of meiosis and associated genome defence mechanisms. These species use parasexuality instead of sex to combine beneficial mutations from separate clonal individuals into one genome (Vicar of Bray dynamics). More weakly endophytic species which kill a narrow range of insects retain sexuality to facilitate host-pathogen coevolution (Red Queen dynamics). Metarhizium species can fit into numerous environments because they are very flexible at the genetic, physiological and ecological levels, providing tractable models to address how new mechanisms for econutritional heterogeneity, host switching and virulence are acquired and relate to diverse sexual life histories and speciation. Many new molecules and functions have been discovered that underpin Metarhizium associations, and have furthered our understanding of the crucial ecology of these fungi in multiple habitats.

The insertion of a mitochondrial selfish element into the nuclear genome and its consequences

Homing endonucleases (HE) are enzymes capable of cutting DNA at highly specific target sequences, the repair of the generated double-strand break resulting in the insertion of the HE-encoding gene ("homing" mechanism). HEs are present in all three domains of life and viruses; in eukaryotes, they are mostly found in the genomes of mitochondria and chloroplasts, as well as nuclear ribosomal RNAs. We here report the case of a HE that accidentally integrated into a telomeric region of the nuclear genome of the fungal maize pathogen Ustilago maydis. We show that the gene has a mitochondrial origin, but its original copy is absent from the U. maydis mitochondrial genome, suggesting a subsequent loss or a horizontal transfer from a different species. The telomeric HE underwent mutations in its active site and lost its original start codon. A potential other start codon was retained downstream, but we did not detect any significant transcription of the newly created open reading frame, suggesting that the inserted gene is not functional. Besides, the insertion site is located in a putative RecQ helicase gene, truncating the C-terminal domain of the protein. The truncated helicase is expressed during infection of the host, together with other homologous telomeric helicases. This unusual mutational event altered two genes: The integrated HE gene subsequently lost its homing activity, while its insertion created a truncated version of an existing gene, possibly altering its function. As the insertion is absent in other field isolates, suggesting that it is recent, the U. maydis 521 reference strain offers a snapshot of this singular mutational event.

The effect of hybridization on transposable element accumulation in an undomesticated fungal species

Transposable elements (TEs) are mobile genetic elements that can profoundly impact the evolution of genomes and species. A long-standing hypothesis suggests that hybridization could deregulate TEs and trigger their accumulation, although it received mixed support from studies mostly in plants and animals. Here, we tested this hypothesis in fungi using incipient species of the undomesticated yeast Saccharomyces paradoxus. Population genomic data revealed no signature of higher transposition in natural hybrids. As we could not rule out the elimination of past transposition increase signatures by natural selection, we performed a laboratory evolution experiment on a panel of artificial hybrids to measure TE accumulation in the near absence of selection. Changes in TE copy numbers were not predicted by the level of evolutionary divergence between the parents of a hybrid genotype. Rather, they were highly dependent on the individual hybrid genotypes, showing that strong genotype-specific deterministic factors govern TE accumulation in yeast hybrids.

Accessory Chromosomes in Fusarium oxysporum

Most genomes within the species complex of Fusarium oxysporum are organized into two compartments: the core chromosomes (CCs) and accessory chromosomes (ACs). As opposed to CCs, which are conserved and vertically transmitted to carry out essential housekeeping functions, lineage- or strain-specific ACs are believed to be initially horizontally acquired through unclear mechanisms. These two genomic compartments are different in terms of gene density, the distribution of transposable elements, and epigenetic markers. Although common in eukaryotes, the functional importance of ACs is uniquely emphasized among fungal species, specifically in relationship to fungal pathogenicity and their adaptation to diverse hosts. With a focus on the cross-kingdom fungal pathogen F. oxysporum, this review provides a summary of the differences between CCs and ACs based on current knowledge of gene functions, genome structures, and epigenetic signatures, and explores the transcriptional crosstalk between the core and accessory genomes.