Whole genome comparative analysis of transposable elements provides new insight into mechanisms of their inactivation in fungal genomes

Section-level genome sequencing and comparative genomics of Aspergillus sections Cavernicolus and Usti

The genus Aspergillus is diverse, including species of industrial importance, human pathogens, plant pests, and model organisms. Aspergillus includes species from sections Usti and Cavernicolus, which until recently were joined in section Usti, but have now been proposed to be non-monophyletic and were split by section Nidulantes, Aenei and Raperi. To learn more about these sections, we have sequenced the genomes of 13 Aspergillus species from section Cavernicolus (A. cavernicola, A. californicus, and A. egyptiacus), section Usti (A. carlsbadensis, A. germanicus, A. granulosus, A. heterothallicus, A. insuetus, A. keveii, A. lucknowensis, A. pseudodeflectus and A. pseudoustus), and section Nidulantes (A. quadrilineatus, previously A. tetrazonus). We compared these genomes with 16 additional species from Aspergillus to explore their genetic diversity, based on their genome content, repeat-induced point mutations (RIPs), transposable elements, carbohydrate-active enzyme (CAZyme) profile, growth on plant polysaccharides, and secondary metabolite gene clusters (SMGCs). All analyses support the split of section Usti and provide additional insights: Analyses of genes found only in single species show that these constitute genes which appear to be involved in adaptation to new carbon sources, regulation to fit new niches, and bioactive compounds for competitive advantages, suggesting that these support species differentiation in Aspergillus species. Sections Usti and Cavernicolus have mainly unique SMGCs. Section Usti contains very large and information-rich genomes, an expansion partially driven by CAZymes, as section Usti contains the most CAZyme-rich species seen in genus Aspergillus. Section Usti is clearly an underutilized source of plant biomass degraders and shows great potential as industrial enzyme producers. Citation: Nybo JL, Vesth TC, Theobald S, Frisvad JC, Larsen TO, Kjaerboelling I, Rothschild-Mancinelli K, Lyhne EK, Barry K, Clum A, Yoshinaga Y, Ledsgaard L, Daum C, Lipzen A, Kuo A, Riley R, Mondo S, LaButti K, Haridas S, Pangalinan J, Salamov AA, Simmons BA, Magnuson JK, Chen J, Drula E, Henrissat B, Wiebenga A, Lubbers RJM, Müller A, dos Santos Gomes AC, Mäkelä MR, Stajich JE, Grigoriev IV, Mortensen UH, de Vries RP, Baker SE, Andersen MR (2025). Section-level genome sequencing and comparative genomics of Aspergillus sections Cavernicolus and Usti. Studies in Mycology 111: 101-114. doi: 10.3114/sim.2025.111.03.

Evolutionary History of the DD41D Family of Tc1/Mariner Transposons in Two Mayetiola Species

Tc1/mariner elements are ubiquitous in eukaryotic genomes including insects. They are diverse and divided into families and sub-families. The DD34D family including mauritiana and irritans subfamilies have already been identified in two closely related species of Cecidomyiids M. destructor and M. hordei. In the current study the de novo and similarity-based methods allowed the identification for the first time of seven consensuses in M. destructor and two consensuses in M. hordei belonging to DD41D family whereas the in vitro method allowed the amplification of two and three elements in these two species respectively. Most of identified elements accumulated different mutations and long deletions spanning the N-terminal region of the transposase. Phylogenetic analyses showed that the DD41D elements were clustered in two groups belonging to rosa and Long-TIR subfamilies. The age estimation of the last transposition events of the identified Tc1/mariner elements in M. destructor showed different evolutionary histories. Indeed, irritans elements have oscillated between periods of silencing and reappearance while rosa and mauritiana elements have shown regular activity with large recent bursts. The study of insertion sites showed that they are mostly intronic and that some recently transposed elements occurred in genes linked to putative DNA-binding domains and enzymes involved in metabolic chains. Thus, this study gave evidence of the existence of DD41D family in two Mayetiola species and an insight on their evolutionary history.

Complete genome of the Medicago anthracnose fungus, Colletotrichum destructivum, reveals a mini-chromosome-like region within a core chromosome

Colletotrichum destructivum (Cd) is a phytopathogenic fungus causing significant economic losses on forage legume crops (Medicago and Trifolium species) worldwide. To gain insights into the genetic basis of fungal virulence and host specificity, we sequenced the genome of an isolate from Medicago sativa using long-read (PacBio) technology. The resulting genome assembly has a total length of 51.7 Mb and comprises ten core chromosomes and two accessory chromosomes, all of which were sequenced from telomere to telomere. A total of 15, 631 gene models were predicted, including genes encoding potentially pathogenicity-related proteins such as candidate-secreted effectors (484), secondary metabolism key enzymes (110) and carbohydrate-active enzymes (619). Synteny analysis revealed extensive structural rearrangements in the genome of Cd relative to the closely related Brassicaceae pathogen, Colletotrichum higginsianum. In addition, a 1.2 Mb species-specific region was detected within the largest core chromosome of Cd that has all the characteristics of fungal accessory chromosomes (transposon-rich, gene-poor, distinct codon usage), providing evidence for exchange between these two genomic compartments. This region was also unique in having undergone extensive intra-chromosomal segmental duplications. Our findings provide insights into the evolution of accessory regions and possible mechanisms for generating genetic diversity in this asexual fungal pathogen.

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 chromosome-scale genome assembly of the grape powdery mildew pathogen Erysiphe necator reveals its genomic architecture and previously unknown features of its biology

Erysiphe necator is an obligate fungal pathogen that causes grape powdery mildew, globally the most important disease on grapevines. Previous attempts to obtain a quality genome assembly for this pathogen were hindered by its high repetitive DNA content. Here, chromatin conformation capture (Hi-C) with long-read PacBio sequencing was combined to obtain a chromosome-scale assembly and a high-quality annotation for E. necator isolate EnFRAME01. The resulting 81.1 Mb genome assembly is 98% complete and consists of 34 scaffolds, 11 of which represent complete chromosomes. All chromosomes contain large centromeric-like regions and lack synteny to the 11 chromosomes of the cereal PM pathogen Blumeria graminis. Further analysis of their composition showed that repeats and transposable elements (TEs) occupy 62.7% of their content. TEs were almost evenly interspersed outside centromeric and telomeric regions and massively overlapped with regions of annotated genes, suggesting that they could have a significant functional impact. Abundant gene duplicates were observed as well, particularly in genes encoding candidate secreted effector proteins. Moreover, younger in age gene duplicates exhibited more relaxed selection pressure and were more likely to be located physically close in the genome than older duplicates. A total of 122 genes with copy number variations among six isolates of E. necator were also identified and were enriched in genes that were duplicated in EnFRAME01, indicating they may reflect an adaptive variation. Taken together, our study illuminates higher-order genomic architectural features of E. necator and provides a valuable resource for studying genomic structural variations in this pathogen. IMPORTANCE Grape powdery mildew caused by the ascomycete fungus Erysiphe necator is economically the most important and recurrent disease in vineyards across the world. The obligate biotrophic nature of E. necator hinders the use of typical genetic methods to elucidate its pathogenicity and adaptation to adverse conditions, and thus comparative genomics has been a major method to study its genome biology. However, the current reference genome of E. necator isolate C-strain is highly fragmented with many non-coding regions left unassembled. This incompleteness prohibits in-depth comparative genomic analyses and the study of genomic structural variations (SVs) that are known to affect several aspects of microbial life, including fitness, virulence, and host adaptation. By obtaining a chromosome-scale genome assembly and a high-quality gene annotation for E. necator, we reveal the organization of its chromosomal content, unearth previously unknown features of its biology, and provide a reference for studying genomic SVs in this pathogen.

Genome characterization of two novel deep-sea sediment fungi, Penicillium pacificagyrus sp. nov. and Penicillium pacificasedimenti sp. nov., from South Pacific Gyre subseafloor sediments, highlights survivability

BACKGROUND

Marine deep subsurface sediments were once thought to be devoid of eukaryotic life, but advances in molecular technology have unlocked the presence and activity of well-known closely related terrestrial and marine fungi. Commonly detected fungi in deep marine sediment environments includes Penicillium, Aspergillus, Cladosporium, Fusarium, and Schizophyllum, which could have important implications in carbon and nitrogen cycling in this isolated environment. In order to determine the diversity and unknown metabolic capabilities of fungi in deep-sea sediments, their genomes need to be fully analyzed. In this study, two Penicillium species were isolated from South Pacific Gyre sediment enrichments during Integrated Ocean Drilling Program Expedition 329. The inner gyre has very limited productivity, organic carbon, and nutrients.

RESULTS

Here, we present high-quality genomes of two proposed novel Penicillium species using Illumina HiSeq and PacBio sequencing technologies. Single-copy homologues within the genomes were compared to other closely related genomes using OrthoMCL and maximum-likelihood estimation, which showed that these genomes were novel species within the genus Penicillium. We propose to name isolate SPG-F1 as Penicillium pacificasedimenti sp. nov. and SPG-F15 as Penicillium pacificagyrus sp. nov. The resulting genome sizes were 32.6 Mbp and 36.4 Mbp, respectively, and both genomes were greater than 98% complete as determined by the presence of complete single-copy orthologs. The transposable elements for each genome were 4.87% for P. pacificasedimenti and 10.68% for P. pacificagyrus. A total of 12,271 genes were predicted in the P. pacificasedimenti genome and 12,568 genes in P. pacificagyrus. Both isolates contained genes known to be involved in the degradation of recalcitrant carbon, amino acids, and lignin-derived carbon.

CONCLUSIONS

Our results provide the first constructed genomes of novel Penicillium isolates from deep marine sediments, which will be useful for future studies of marine subsurface fungal diversity and function. Furthermore, these genomes shed light on the potential impact fungi in marine sediments and the subseafloor could have on global carbon and nitrogen biogeochemical cycles and how they may be persisting in the most energy-limited sedimentary biosphere.

Complete Genome Sequences of Septoria linicola: A Resource for Studying a Damaging Flax Pathogen

Fungal genus Septoria causes diseases in a wide range of plants. Here, we report the first genome sequences of two strains of Septoria linicola, the causal agent of the pasmo disease of flax (Linum usitatissimum). The genome of the first strain, SE15195, was fully assembled in 16 chromosomes, while 35 unitigs were obtained for a second strain, SE14017. Structural annotations predicted 13,096 and 13,085 protein-encoding genes and transposable elements content of 19.0 and 18.1% of the genome for SE15195 and SE14017, respectively. The four smaller chromosomes 13 to 16 show genomics features of potential accessory chromosomes. The assembly of these two genomes is a new resource for studying S. linicola and improving management of pasmo. [Formula: see text] Copyright © 2022 The Author(s). This is an open access article distributed under the CC BY-NC-ND 4.0 International license.

Evolution of pathogenicity in obligate fungal pathogens and allied genera

Obligate fungal pathogens (ascomycetes and basidiomycetes) and oomycetes are known to cause diseases in cereal crop plants. They feed on living cells and most of them have learned to bypass the host immune machinery. This paper discusses some of the factors that are associated with pathogenicity drawing examples from ascomycetes, basidiomycetes and oomycetes, with respect to their manifestation in crop plants. The comparisons have revealed a striking similarity in the three groups suggesting convergent pathways that have arisen from three lineages independently leading to an obligate lifestyle. This review has been written with the intent, that new information on adaptation strategies of biotrophs, modifications in pathogenicity strategies and population dynamics will improve current strategies for breeding with stable resistance.

Comparative Genome Analyses of Plant Rust Pathogen Genomes Reveal a Confluence of Pathogenicity Factors to Quell Host Plant Defense Responses

Switchgrass rust caused by Puccinia novopanici (P. novopanici) has the ability to significantly affect the biomass yield of switchgrass, an important biofuel crop in the United States. A comparative genome analysis of P. novopanici with rust pathogen genomes infecting monocot cereal crops wheat, barley, oats, maize and sorghum revealed the presence of larger structural variations contributing to their genome sizes. A comparative alignment of the rust pathogen genomes resulted in the identification of collinear and syntenic relationships between P. novopanici and P. sorghi; P. graminis tritici 21–0 (Pgt 21) and P. graminis tritici Ug99 (Pgt Ug99) and between Pgt 21 and P. triticina (Pt). Repeat element analysis indicated a strong presence of retro elements among different Puccinia genomes, contributing to the genome size variation between ~1 and 3%. A comparative look at the enriched protein families of Puccinia spp. revealed a predominant role of restriction of telomere capping proteins (RTC), disulfide isomerases, polysaccharide deacetylases, glycoside hydrolases, superoxide dismutases and multi-copper oxidases (MCOs). All the proteomes of Puccinia spp. share in common a repertoire of 75 secretory and 24 effector proteins, including glycoside hydrolases cellobiohydrolases, peptidyl-propyl isomerases, polysaccharide deacetylases and protein disulfide-isomerases, that remain central to their pathogenicity. Comparison of the predicted effector proteins from Puccinia spp. genomes to the validated proteins from the Pathogen–Host Interactions database (PHI-base) resulted in the identification of validated effector proteins PgtSR1 (PGTG_09586) from P. graminis and Mlp124478 from Melampsora laricis across all the rust pathogen genomes.

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.

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.

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.

Dynamics of transposable elements in recently diverged fungal pathogens: lineage-specific transposable element content and efficiency of genome defenses

Transposable elements (TEs) impact genome plasticity, architecture, and evolution in fungal plant pathogens. The wide range of TE content observed in fungal genomes reflects diverse efficacy of host-genome defense mechanisms that can counter-balance TE expansion and spread. Closely related species can harbor drastically different TE repertoires. The evolution of fungal effectors, which are crucial determinants of pathogenicity, has been linked to the activity of TEs in pathogen genomes. Here, we describe how TEs have shaped genome evolution of the fungal wheat pathogen Zymoseptoria tritici and four closely related species. We compared de novo TE annotations and repeat-induced point mutation signatures in 26 genomes from the Zymoseptoria species-complex. Then, we assessed the relative insertion ages of TEs using a comparative genomics approach. Finally, we explored the impact of TE insertions on genome architecture and plasticity. The 26 genomes of Zymoseptoria species reflect different TE dynamics with a majority of recent insertions. TEs associate with accessory genome compartments, with chromosomal rearrangements, with gene presence/absence variation, and with effectors in all Zymoseptoria species. We find that the extent of RIP-like signatures varies among Z. tritici genomes compared to genomes of the sister species. The detection of a reduction of RIP-like signatures and TE recent insertions in Z. tritici reflects ongoing but still moderate TE mobility.

The Assembled and Annotated Genome of the Fairy-Ring Fungus Marasmius oreades

Marasmius oreades is a basidiomycete fungus that grows in so called "fairy rings," which are circular, underground mycelia common in lawns across temperate areas of the world. Fairy rings can be thought of as natural, long-term evolutionary experiments. As each ring has a common origin and expands radially outwards over many years, different sectors will independently accumulate mutations during growth. The genotype can be followed to the next generation, as mushrooms producing the sexual spores are formed seasonally at the edge of the ring. Here, we present new genomic data from 95 single-spore isolates of the species, which we used to construct a genetic linkage map and an updated version of the genome assembly. The 44-Mb assembly was anchored to 11 linkage groups, producing chromosome-length scaffolds. Gene annotation revealed 13,891 genes, 55% of which contained a pfam domain. The repetitive fraction of the genome was 22%, and dominated by retrotransposons and DNA elements of the KDZ and Plavaka groups. The level of assembly contiguity we present is so far rare in mushroom-forming fungi, and we expect studies of genomics, transposons, phylogenetics, and evolution to be facilitated by the data we present here of the iconic fairy-ring mushroom.

Transposable Elements Contribute to Genome Dynamics and Gene Expression Variation in the Fungal Plant Pathogen Verticillium dahliae

Transposable elements (TEs) are a major source of genetic and regulatory variation in their host genome and are consequently thought to play important roles in evolution. Many fungal and oomycete plant pathogens have evolved dynamic and TE-rich genomic regions containing genes that are implicated in host colonization and adaptation. TEs embedded in these regions have typically been thought to accelerate the evolution of these genomic compartments, but little is known about their dynamics in strains that harbor them. Here, we used whole-genome sequencing data of 42 strains of the fungal plant pathogen Verticillium dahliae to systematically identify polymorphic TEs that may be implicated in genomic as well as in gene expression variation. We identified 2,523 TE polymorphisms and characterize a subset of 8% of the TEs as polymorphic elements that are evolutionary younger, less methylated, and more highly expressed when compared with the remaining 92% of the total TE complement. As expected, the polyrmorphic TEs are enriched in the adaptive genomic regions. Besides, we observed an association of polymorphic TEs with pathogenicity-related genes that localize nearby and that display high expression levels. Collectively, our analyses demonstrate that TE dynamics in V. dahliae contributes to genomic variation, correlates with expression of pathogenicity-related genes, and potentially impacts the evolution of adaptive genomic regions.

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.

A chromosome-scale assembly of the smallest Dothideomycete genome reveals a unique genome compaction mechanism in filamentous fungi

Background

The wide variation in the size of fungal genomes is well known, but the reasons for this size variation are less certain. Here, we present a chromosome-scale assembly of ectophytic Peltaster fructicola, a surface-dwelling extremophile, based on long-read DNA sequencing technology, to assess possible mechanisms associated with genome compaction.

Results

At 18.99 million bases (Mb), P. fructicola possesses one of the smallest known genomes sequence among filamentous fungi. The genome is highly compact relative to other fungi, with substantial reductions in repeat content, ribosomal DNA copies, tRNA gene quantity, and intron sizes, as well as intergenic lengths and the size of gene families. Transposons take up just 0.05% of the entire genome, and no full-length transposon was found. We concluded that reduced genome sizes in filamentous fungi such as P. fructicola, Taphrina deformans and Pneumocystis jirovecii occurred through reduction in ribosomal DNA copy number and reduced intron sizes. These dual mechanisms contrast with genome reduction in the yeast fungus Saccharomyces cerevisiae, whose small and compact genome is associated solely with intron loss.

Conclusions

Our results reveal a unique genomic compaction architecture of filamentous fungi inhabiting plant surfaces, and broaden the understanding of the mechanisms associated with compaction of fungal genomes.

Unusual genome expansion and transcription suppression in ectomycorrhizal Tricholoma matsutake by insertions of transposable elements

Genome sequencing of Tricholoma matsutake revealed its unusually large size as 189.0 Mbp, which is a consequence of extraordinarily high transposable element (TE) content. We identified that 702 genes were surrounded by TEs, and 83.2% of these genes were not transcribed at any developmental stage. This observation indicated that the insertion of TEs alters the transcription of the genes neighboring these TEs. Repeat-induced point mutation, such as C to T hypermutation with a bias over "CpG" dinucleotides, was also recognized in this genome, representing a typical defense mechanism against TEs during evolution. Many transcription factor genes were activated in both the primordia and fruiting body stages, which indicates that many regulatory processes are shared during the developmental stages. Small secreted protein genes (<300 aa) were dominantly transcribed in the hyphae, where symbiotic interactions occur with the hosts. Comparative analysis with 37 Agaricomycetes genomes revealed that IstB-like domains (PF01695) were conserved across taxonomically diverse mycorrhizal genomes, where the T. matsutake genome contained four copies of this domain. Three of the IstB-like genes were overexpressed in the hyphae. Similar to other ectomycorrhizal genomes, the CAZyme gene set was reduced in T. matsutake, including losses in the glycoside hydrolase genes. The T. matsutake genome sequence provides insight into the causes and consequences of genome size inflation.

Repeat-Induced Point Mutations Drive Divergence between Fusarium circinatum and Its Close Relatives

The Repeat-Induced Point (RIP) mutation pathway is a fungal-specific genome defense mechanism that counteracts the deleterious effects of transposable elements. This pathway permanently mutates its target sequences by introducing cytosine to thymine transitions. We investigated the genome-wide occurrence of RIP in the pitch canker pathogen, Fusarium circinatum, and its close relatives in the Fusarium fujikuroi species complex (FFSC). Our results showed that the examined fungi all exhibited hallmarks of RIP, but that they differed in terms of the extent to which their genomes were affected by this pathway. RIP mutations constituted a large proportion of all the FFSC genomes, including both core and dispensable chromosomes, although the latter were generally more extensively affected by RIP. Large RIP-affected genomic regions were also much more gene sparse than the rest of the genome. Our data further showed that RIP-directed sequence diversification increased the variability between homologous regions of related species, and that RIP-affected regions can interfere with homologous recombination during meiosis, thereby contributing to post-mating segregation distortion. Taken together, these findings suggest that RIP can drive the independent divergence of chromosomes, alter chromosome architecture, and contribute to the divergence among F. circinatum and other members of this economically important group of fungi.

Direct Homologous dsDNA-dsDNA Pairing: How, Where, and Why?

The ability of homologous chromosomes (or selected chromosomal loci) to pair specifically in the apparent absence of DNA breakage and recombination represents a prominent feature of eukaryotic biology. The mechanism of homology recognition at the basis of such recombination-independent pairing has remained elusive. A number of studies have supported the idea that sequence homology can be sensed between intact DNA double helices in vivo. In particular, recent analyses of the two silencing phenomena in fungi, known as "repeat-induced point mutation" (RIP) and "meiotic silencing by unpaired DNA" (MSUD), have provided genetic evidence for the existence of the direct homologous dsDNA-dsDNA pairing. Both RIP and MSUD likely rely on the same search strategy, by which dsDNA segments are matched as arrays of interspersed base-pair triplets. This process is general and very efficient, yet it proceeds normally without the RecA/Rad51/Dmc1 proteins. Further studies of RIP and MSUD may yield surprising insights into the function of DNA in the cell.

At the nexus of three kingdoms: the genome of the mycorrhizal fungus Gigaspora margarita provides insights into plant, endobacterial and fungal interactions

As members of the plant microbiota, arbuscular mycorrhizal fungi (AMF, Glomeromycotina) symbiotically colonize plant roots. AMF also possess their own microbiota, hosting some uncultivable endobacteria. Ongoing research has revealed the genetics underlying plant responses to colonization by AMF, but the fungal side of the relationship remains in the dark. Here, we sequenced the genome of Gigaspora margarita, a member of the Gigasporaceae in an early diverging group of the Glomeromycotina. In contrast to other AMF, G. margarita may host distinct endobacterial populations and possesses the largest fungal genome so far annotated (773.104 Mbp), with more than 64% transposable elements. Other unique traits of the G. margarita genome include the expansion of genes for inorganic phosphate metabolism, the presence of genes for production of secondary metabolites and a considerable number of potential horizontal gene transfer events. The sequencing of G. margarita genome reveals the importance of its immune system, shedding light on the evolutionary pathways that allowed early diverging fungi to interact with both plants and bacteria.

The Parauncinula polyspora Draft Genome Provides Insights into Patterns of Gene Erosion and Genome Expansion in Powdery Mildew Fungi

Powdery mildew fungi are widespread and agronomically relevant phytopathogens causing major yield losses. Their genomes have disproportionately large numbers of mobile genetic elements, and they have experienced a significant loss of highly conserved fungal genes. In order to learn more about the evolutionary history of this fungal group, we explored the genome of an Asian oak tree pathogen, Parauncinula polyspora, a species that diverged early during evolution from the remaining powdery mildew fungi. We found that the P. polyspora draft genome is comparatively compact, has a low number of protein-coding genes, and, despite the absence of a dedicated genome defense system, lacks the massive proliferation of repetitive sequences. Based on these findings, we infer an evolutionary trajectory that shaped the genomes of powdery mildew fungi.

Due to their comparatively small genome size and short generation time, fungi are exquisite model systems to study eukaryotic genome evolution. Powdery mildew fungi present an exceptional case because of their strict host dependency (termed obligate biotrophy) and the atypical size of their genomes (>100 Mb). This size expansion is largely due to the pervasiveness of transposable elements on 70% of the genome and is associated with the loss of multiple conserved ascomycete genes required for a free-living lifestyle. To date, little is known about the mechanisms that drove these changes, and information on ancestral powdery mildew genomes is lacking. We report genome analysis of the early-diverged and exclusively sexually reproducing powdery mildew fungus Parauncinula polyspora, which we performed on the basis of a natural leaf epiphytic metapopulation sample. In contrast to other sequenced species of this taxonomic group, the assembled P. polyspora draft genome is surprisingly small (<30 Mb), has a higher content of conserved ascomycete genes, and is sparsely equipped with transposons (<10%), despite the conserved absence of a common defense mechanism involved in constraining repetitive elements. We speculate that transposable element spread might have been limited by this pathogen’s unique reproduction strategy and host features and further hypothesize that the loss of conserved ascomycete genes may promote the evolutionary isolation and host niche specialization of powdery mildew fungi. Limitations associated with this evolutionary trajectory might have been in part counteracted by the evolution of plastic, transposon-rich genomes and/or the expansion of gene families encoding secreted virulence proteins.

A RID-like putative cytosine methyltransferase homologue controls sexual development in the fungus Podospora anserina

DNA methyltransferases are ubiquitous enzymes conserved in bacteria, plants and opisthokonta. These enzymes, which methylate cytosines, are involved in numerous biological processes, notably development. In mammals and higher plants, methylation patterns established and maintained by the cytosine DNA methyltransferases (DMTs) are essential to zygotic development. In fungi, some members of an extensively conserved fungal-specific DNA methyltransferase class are both mediators of the Repeat Induced Point mutation (RIP) genome defense system and key players of sexual reproduction. Yet, no DNA methyltransferase activity of these purified RID (RIP deficient) proteins could be detected in vitro. These observations led us to explore how RID-like DNA methyltransferase encoding genes would play a role during sexual development of fungi showing very little genomic DNA methylation, if any. To do so, we used the model ascomycete fungus Podospora anserina. We identified the PaRid gene, encoding a RID-like DNA methyltransferase and constructed knocked-out ΔPaRid defective mutants. Crosses involving P. anserina ΔPaRid mutants are sterile. Our results show that, although gametes are readily formed and fertilization occurs in a ΔPaRid background, sexual development is blocked just before the individualization of the dikaryotic cells leading to meiocytes. Complementation of ΔPaRid mutants with ectopic alleles of PaRid, including GFP-tagged, point-mutated and chimeric alleles, demonstrated that the catalytic motif of the putative PaRid methyltransferase is essential to ensure proper sexual development and that the expression of PaRid is spatially and temporally restricted. A transcriptomic analysis performed on mutant crosses revealed an overlap of the PaRid-controlled genetic network with the well-known mating-types gene developmental pathway common to an important group of fungi, the Pezizomycotina.

Genomic analysis of the aggressive tree pathogen Ceratocystis albifundus

The overall goal of this study was to determine whether the genome of an important plant pathogen in Africa, Ceratocystis albifundus, is structured into subgenomic compartments, and if so, to establish how these compartments are distributed across the genome. For this purpose, the publicly available genome of C. albifundus was complemented with the genome sequences for four additional isolates using the Illumina HiSeq platform. In addition, a reference genome for one of the individuals was assembled using both PacBio and Illumina HiSeq technologies. Our results showed a high degree of synteny between the five genomes, although several regions lacked detectable long-range synteny. These regions were associated with the presence of accessory genes, lower genetic similarity, variation in read-map depth, as well as transposable elements and genes associated with host-pathogen interactions (e.g. effectors and CAZymes). Such patterns are regarded as hallmarks of accelerated evolution, particularly of accessory subgenomic compartments in fungal pathogens. Our findings thus showed that the genome of C. albifundus is made-up of core and accessory subgenomic compartments, which is an important step towards characterizing its pangenome. This study also highlights the value of comparative genomics for understanding mechanisms that may underly and influence the biology and evolution of pathogens.

RepetDB: a unified resource for transposable element references

Background

Thanks to their ability to move around and replicate within genomes, transposable elements (TEs) are perhaps the most important contributors to genome plasticity and evolution. Their detection and annotation are considered essential in any genome sequencing project. The number of fully sequenced genomes is rapidly increasing with improvements in high-throughput sequencing technologies. A fully automated de novo annotation process for TEs is therefore required to cope with the deluge of sequence data.

However, all automated procedures are error-prone, and an automated procedure for TE identification and classification would be no exception. It is therefore crucial to provide not only the TE reference sequences, but also evidence justifying their classification, at the scale of the whole genome. A few TE databases already exist, but none provides evidence to justify TE classification. Moreover, biological information about the sequences remains globally poor.

Results

We present here the RepetDB database developed in the framework of GnpIS, a genetic and genomic information system. RepetDB is designed to store and retrieve detected, classified and annotated TEs in a standardized manner. RepetDB is an implementation with extensions of InterMine, an open-source data warehouse framework used here to store, search, browse, analyze and compare all the data recorded for each TE reference sequence. InterMine can display diverse information for each sequence and allows simple to very complex queries. Finally, TE data are displayed via a worldwide data discovery portal. RepetDB is accessible at urgi.versailles.inra.fr/repetdb.

Conclusions

RepetDB is designed to be a TE knowledge base populated with full de novo TE annotations of complete (or near-complete) genome sequences. Indeed, the description and classification of TEs facilitates the exploration of specific TE families, superfamilies or orders across a large range of species. It also makes possible cross-species searches and comparisons of TE family content between genomes.

Analysis of Repeat Induced Point (RIP) Mutations in Leptosphaeria maculans Indicates Variability in the RIP Process Between Fungal Species

Gene duplication contributes to evolutionary potential, yet many duplications in a genome arise from the activity of "selfish" genetic elements such as transposable elements. Fungi have a number of mechanisms by which they limit the expansion of transposons, including Repeat Induced Point mutation (RIP). RIP has been best characterized in the Sordariomycete Neurospora crassa, wherein duplicated DNA regions are recognized after cell fusion, but before nuclear fusion during the sexual cycle, and then mutated. While "signatures" of RIP appear in the genome sequences of many fungi, the species most distant from N. crassa in which the process has been experimentally demonstrated to occur is the Dothideomycete Leptosphaeria maculans In the current study, we show that similar to N. crassa, nonlinked duplications can trigger RIP; however, the frequency of the generated RIP mutations is extremely low in L maculans (< 0.1%) and requires a large duplication to initiate RIP, and that multiple premeiotic mitoses are involved in the RIP process. However, a single sexual cycle leads to the generation of progeny with unique haplotypes, despite progeny pairs being generated from mitosis. We hypothesize that these different haplotypes may be the result of the deamination process occurring post karyogamy, leading to unique mutations within each of the progeny pairs. These findings indicate that the RIP process, while common to many fungi, differs between fungi and that this impacts on the fate of duplicated DNA.

Genome-wide DNA methylation and transcriptomic profiles in the lifestyle strategies and asexual development of the forest fungal pathogen Heterobasidion parviporum

Heterobasidion parviporum is the most devastating fungal pathogen of conifer forests in Northern Europe. The fungus has dual life strategies, necrotrophy on living trees and saprotrophy on dead woods. DNA cytosine methylation is an important epigenetic modification in eukaryotic organisms. Our presumption is that the lifestyle transition and asexual development in H. parviporum could be driven by epigenetic effects. Involvements of DNA methylation in the regulation of aforementioned processes have never been studied thus far. RNA-seq identified lists of highly induced genes enriched in carbohydrate-active enzymes during necrotrophic interaction with host trees and saprotrophic sawdust growth. It also highlighted signaling- and transcription factor-related genes potentially associated with the transition of saprotrophic to necrotrophic lifestyle and groups of primary cellular activities throughout asexual development. Whole-genome bisulfite sequencing revealed that DNA methylation displayed pronounced preference in CpG dinucleotide context across the genome and mostly targeted transposable element (TE)-rich regions. TE methylation level demonstrated a strong negative correlation with TE expression, reinforcing the protective function of DNA methylation in fungal genome stability. Small groups of genes putatively subject to methylation transcriptional regulation in response to saprotrophic and necrotrophic growth in comparison with free-living mycelia were also explored. Our study reported on the first methylome map of a forest pathogen. Analysis of transcriptome and methylome variations associated with asexual development and different lifestyle strategies provided further understanding of basic biological processes in H. parviporum. More importantly, our work raised additional potential roles of DNA methylation in fungi apart from controlling the proliferation of TEs.

Genomic overview of closely related fungi with different Protea host ranges

Genome comparisons of species with distinctive ecological traits can elucidate genetic divergence that influenced their differentiation. The interaction of a microorganism with its biotic environment is largely regulated by secreted compounds, and these can be predicted from genome sequences. In this study, we considered Knoxdaviesia capensis and Knoxdaviesia proteae, two closely related saprotrophic fungi found exclusively in Protea plants. We investigated their genome structure to compare their potential inter-specific interactions based on gene content. Their genomes displayed macrosynteny and were approximately 10 % repetitive. Both species had fewer secreted proteins than pathogens and other saprotrophs, reflecting their specialized habitat. The bulk of the predicted species-specific and secreted proteins coded for carbohydrate metabolism, with a slightly higher number of unique carbohydrate-degrading proteins in the broad host-range K. capensis. These fungi have few secondary metabolite gene clusters, suggesting minimal competition with other microbes and symbiosis with antibiotic-producing bacteria common in this niche. Secreted proteins associated with detoxification and iron sequestration likely enable these Knoxdaviesia species to tolerate antifungal compounds and compete for resources, facilitating their unusual dominance. This study confirms the genetic cohesion between Protea-associated Knoxdaviesia species and reveals aspects of their ecology that have likely evolved in response to their specialist niche.

The Landscape of Repetitive Elements in the Refined Genome of Chilli Anthracnose Fungus Colletotrichum truncatum

The ascomycete fungus Colletotrichum truncatum is a major phytopathogen with a broad host range which causes anthracnose disease of chilli. The genome sequencing of this fungus led to the discovery of functional categories of genes that may play important roles in fungal pathogenicity. However, the presence of gaps in C. truncatum draft assembly prevented the accurate prediction of repetitive elements, which are the key players to determine the genome architecture and drive evolution and host adaptation. We re-sequenced its genome using single-molecule real-time (SMRT) sequencing technology to obtain a refined assembly with lesser and smaller gaps and ambiguities. This enabled us to study its genome architecture by characterising the repetitive sequences like transposable elements (TEs) and simple sequence repeats (SSRs), which constituted 4.9 and 0.38% of the assembled genome, respectively. The comparative analysis among different Colletotrichum species revealed the extensive repeat rich regions, dominated by Gypsy superfamily of long terminal repeats (LTRs), and the differential composition of SSRs in their genomes. Our study revealed a recent burst of LTR amplification in C. truncatum, C. higginsianum, and C. scovillei. TEs in C. truncatum were significantly associated with secretome, effectors and genes in secondary metabolism clusters. Some of the TE families in C. truncatum showed cytosine to thymine transitions indicative of repeat-induced point mutation (RIP). C. orbiculare and C. graminicola showed strong signatures of RIP across their genomes and "two-speed" genomes with extensive AT-rich and gene-sparse regions. Comparative genomic analyses of Colletotrichum species provided an insight into the species-specific SSR profiles. The SSRs in the coding and non-coding regions of the genome revealed the composition of trinucleotide repeat motifs in exons with potential to alter the translated protein structure through amino acid repeats. This is the first genome-wide study of TEs and SSRs in C. truncatum and their comparative analysis with six other Colletotrichum species, which would serve as a useful resource for future research to get insights into the potential role of TEs in genome expansion and evolution of Colletotrichum fungi and for development of SSR-based molecular markers for population genomic studies.

Small RNAs from cereal powdery mildew pathogens may target host plant genes

Small RNAs (sRNAs) play a key role in eukaryotic gene regulation, for example by gene silencing via RNA interference (RNAi). The biogenesis of sRNAs depends on proteins that are generally conserved in all eukaryotic lineages, yet some species that lack part or all the components of the mechanism exist. Here we explored the presence of the RNAi machinery and its expression as well as the occurrence of sRNA candidates and their putative endogenous as well as host targets in phytopathogenic powdery mildew fungi. We focused on the species Blumeria graminis, which occurs in various specialized forms (formae speciales) that each have a strictly limited host range. B. graminis f. sp. hordei and B. graminis f. sp. tritici, colonizing barley and wheat, respectively, have genomes that are characterized by extensive gene loss. Nonetheless, we find that the RNAi machinery appears to be largely complete and expressed during infection. sRNA sequencing data enabled the identification of putative sRNAs in both pathogens. While a considerable part of the sRNA candidates have predicted target sites in endogenous genes and transposable elements, a small proportion appears to have targets in planta, suggesting potential cross-kingdom RNA transfer between powdery mildew fungi and their respective plant hosts.

Biosynthesis of abscisic acid in fungi: identification of a sesquiterpene cyclase as the key enzyme in Botrytis cinerea

While abscisic acid (ABA) is known as a hormone produced by plants through the carotenoid pathway, a small number of phytopathogenic fungi are also able to produce this sesquiterpene but they use a distinct pathway that starts with the cyclization of farnesyl diphosphate (FPP) into 2Z,4E-α-ionylideneethane which is then subjected to several oxidation steps. To identify the sesquiterpene cyclase (STC) responsible for the biosynthesis of ABA in fungi, we conducted a genomic approach in Botrytis cinerea. The genome of the ABA-overproducing strain ATCC58025 was fully sequenced and five STC-coding genes were identified. Among them, Bcstc5 exhibits an expression profile concomitant with ABA production. Gene inactivation, complementation and chemical analysis demonstrated that BcStc5/BcAba5 is the key enzyme responsible for the key step of ABA biosynthesis in fungi. Unlike what is observed for most of the fungal secondary metabolism genes, the key enzyme-coding gene Bcstc5/Bcaba5 is not clustered with the other biosynthetic genes, i.e., Bcaba1 to Bcaba4 that are responsible for the oxidative transformation of 2Z,4E-α-ionylideneethane. Finally, our study revealed that the presence of the Bcaba genes among Botrytis species is rare and that the majority of them do not possess the ability to produce ABA.

The genome assembly of the fungal pathogen Pyrenochaeta lycopersici from Single-Molecule Real-Time sequencing sheds new light on its biological complexity

The first draft genome sequencing of the non-model fungal pathogen Pyrenochaeta lycopersici showed an expansion of gene families associated with heterokaryon incompatibility and lacking of mating-type genes, providing insights into the genetic basis of this “imperfect” fungus which lost the ability to produce the sexual stage. However, due to the Illumina short-read technology, the draft genome was too fragmented to allow a comprehensive characterization of the genome, especially of the repetitive sequence fraction. In this work, the sequencing of another P. lycopersici isolate using long-read Single Molecule Real-Time sequencing technology was performed with the aim of obtaining a gapless genome. Indeed, a gapless genome assembly of 62.7 Mb was obtained, with a fraction of repetitive sequences representing 30% of the total bases. The gene content of the two P. lycopersici isolates was very similar, and the large difference in genome size (about 8 Mb) might be attributable to the high fraction of repetitive sequences detected for the new sequenced isolate. The role of repetitive elements, including transposable elements, in modulating virulence effectors is well established in fungal plant pathogens. Moreover, transposable elements are of fundamental importance in creating and re-modelling genes, especially in imperfect fungi. Their abundance in P. lycopersici, together with the large expansion of heterokaryon incompatibility genes in both sequenced isolates, suggest the presence of possible mechanisms alternative to gene re-assorting mediated by sexual recombination. A quite large fraction (~9%) of repetitive elements in P. lycopersici, has no homology with known classes, strengthening this hypothesis. The availability of a gapless genome of P. lycopersici allowed the in-depth analysis of its genome content, by annotating functional genes and TEs. This goal will be an important resource for shedding light on the evolution of the reproductive and pathogenic behaviour of this soilborne pathogen and the onset of a possible speciation within this species.

Novel jack-in-the-box effector of the barley powdery mildew pathogen?

This article comments on: Nottensteiner M, Zechmann B, McCollum C, Hückelhoven R. 2018. A barley powdery mildew fungus non-autonomous retrotransposon encodes a peptide that supports penetration success on barley. Journal of Experimental Botany 69, 3745–3758.

Intraspecific comparative genomics of isolates of the Norway spruce pathogen (Heterobasidion parviporum) and identification of its potential virulence factors

BACKGROUND

Heterobasidion parviporum is an economically most important fungal forest pathogen in northern Europe, causing root and butt rot disease of Norway spruce (Picea abies (L.) Karst.). The mechanisms underlying the pathogenesis and virulence of this species remain elusive. No reference genome to facilitate functional analysis is available for this species.

RESULTS

To better understand the virulence factor at both phenotypic and genomic level, we characterized 15 H. parviporum isolates originating from different locations across Finland for virulence, vegetative growth, sporulation and saprotrophic wood decay. Wood decay capability and latitude of fungal origins exerted interactive effects on their virulence and appeared important for H. parviporum virulence. We sequenced the most virulent isolate, the first full genome sequences of H. parviporum as a reference genome, and re-sequenced the remaining 14 H. parviporum isolates. Genome-wide alignments and intrinsic polymorphism analysis showed that these isolates exhibited overall high genomic similarity with an average of at least 96% nucleotide identity when compared to the reference, yet had remarkable intra-specific level of polymorphism with a bias for CpG to TpG mutations. Reads mapping coverage analysis enabled the classification of all predicted genes into five groups and uncovered two genomic regions exclusively present in the reference with putative contribution to its higher virulence. Genes enriched for copy number variations (deletions and duplications) and nucleotide polymorphism were involved in oxidation-reduction processes and encoding domains relevant to transcription factors. Some secreted protein coding genes based on the genome-wide selection pressure, or the presence of variants were proposed as potential virulence candidates.

CONCLUSION

Our study reported on the first reference genome sequence for this Norway spruce pathogen (H. parviporum). Comparative genomics analysis gave insight into the overall genomic variation among this fungal species and also facilitated the identification of several secreted protein coding genes as putative virulence factors for the further functional analysis. We also analyzed and identified phenotypic traits potentially linked to its virulence.

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

Transposable elements have colonized the genomes of nearly all organisms, including fungi. Although transposable elements may sometimes provide beneficial functions to their hosts their overall impact is considered deleterious. As a result, the activity of transposable elements needs to be counterbalanced by the host genome defenses. In fungi, the primary genome defense mechanisms include repeat-induced point mutation (RIP) and methylation induced premeiotically, meiotic silencing by unpaired DNA, sex-induced silencing, cosuppression (also known as somatic quelling), and cotranscriptional RNA surveillance. Recent studies of the filamentous fungus Neurospora crassa have shown that the process of repeat recognition for RIP apparently involves interactions between coaligned double-stranded segments of chromosomal DNA. These studies have also shown that RIP can be mediated by the conserved pathway that establishes transcriptional (heterochromatic) silencing of repetitive DNA. In light of these new findings, RIP emerges as a specialized case of the general phenomenon of heterochromatic silencing of repetitive DNA.

Global DNA Methylation in the Chestnut Blight Fungus Cryphonectria parasitica and Genome-Wide Changes in DNA Methylation Accompanied with Sectorization

Mutation in CpBck1, an ortholog of the cell wall integrity mitogen-activated protein kinase kinase kinase (MAPKKK) of Saccharomyces cerevisiae, in the chestnut blight fungus Cryphonectria parasitica resulted in a sporadic sectorization as culture proceeded. The progeny from the sectored area maintained the characteristics of the sector, showing a massive morphogenetic change, including robust mycelial growth without differentiation. Epigenetic changes were investigated as the genetic mechanism underlying this sectorization. Quantification of DNA methylation and whole-genome bisulfite sequencing revealed genome-wide DNA methylation of the wild-type at each nucleotide level and changes in DNA methylation of the sectored progeny. Compared to the wild-type, the sectored progeny exhibited marked genome-wide DNA hypomethylation but increased methylation sites. Expression analysis of two DNA methyltransferases, including two representative types of DNA methyltransferase (DNMTase), demonstrated that both were significantly down-regulated in the sectored progeny. However, functional analysis using mutant phenotypes of corresponding DNMTases demonstrated that a mutant of CpDmt1, an ortholog of RID of Neurospora crassa, resulted in the sectored phenotype but the CpDmt2 mutant did not, suggesting that the genetic basis of fungal sectorization is more complex. The present study revealed that a mutation in a signaling pathway component resulted in sectorization accompanied with changes in genome-wide DNA methylation, which suggests that this signal transduction pathway is important for epigenetic control of sectorization via regulation of genes involved in DNA methylation.

Ceratocystis cacaofunesta genome analysis reveals a large expansion of extracellular phosphatidylinositol-specific phospholipase-C genes (PI-PLC)

BACKGROUND

The Ceratocystis genus harbors a large number of phytopathogenic fungi that cause xylem parenchyma degradation and vascular destruction on a broad range of economically important plants. Ceratocystis cacaofunesta is a necrotrophic fungus responsible for lethal wilt disease in cacao. The aim of this work is to analyze the genome of C. cacaofunesta through a comparative approach with genomes of other Sordariomycetes in order to better understand the molecular basis of pathogenicity in the Ceratocystis genus.

RESULTS

We present an analysis of the C. cacaofunesta genome focusing on secreted proteins that might constitute pathogenicity factors. Comparative genome analyses among five Ceratocystidaceae species and 23 other Sordariomycetes fungi showed a strong reduction in gene content of the Ceratocystis genus. However, some gene families displayed a remarkable expansion, in particular, the Phosphatidylinositol specific phospholipases-C (PI-PLC) family. Also, evolutionary rate calculations suggest that the evolution process of this family was guided by positive selection. Interestingly, among the 82 PI-PLCs genes identified in the C. cacaofunesta genome, 70 genes encoding extracellular PI-PLCs are grouped in eight small scaffolds surrounded by transposon fragments and scars that could be involved in the rapid evolution of the PI-PLC family. Experimental secretome using LC-MS/MS validated 24% (86 proteins) of the total predicted secretome (342 proteins), including four PI-PLCs and other important pathogenicity factors.

CONCLUSION

Analysis of the Ceratocystis cacaofunesta genome provides evidence that PI-PLCs may play a role in pathogenicity. Subsequent functional studies will be aimed at evaluating this hypothesis. The observed genetic arsenals, together with the analysis of the PI-PLC family shown in this work, reveal significant differences in the Ceratocystis genome compared to the classical vascular fungi, Verticillium and Fusarium. Altogether, our analyses provide new insights into the evolution and the molecular basis of plant pathogenicity.

Describing Genomic and Epigenomic Traits Underpinning Emerging Fungal Pathogens

An unprecedented number of pathogenic fungi are emerging and causing disease in animals and plants, putting the resilience of wild and managed ecosystems in jeopardy. While the past decades have seen an increase in the number of pathogenic fungi, they have also seen the birth of new big data technologies and analytical approaches to tackle these emerging pathogens. We review how the linked fields of genomics and epigenomics are transforming our ability to address the challenge of emerging fungal pathogens. We explore the methodologies and bioinformatic toolkits that currently exist to rapidly analyze the genomes of unknown fungi, then discuss how these data can be used to address key questions that shed light on their epidemiology. We show how genomic approaches are leading a revolution into our understanding of emerging fungal diseases and speculate on future approaches that will transform our ability to tackle this increasingly important class of emerging pathogens.

Plasticity of the MFS1 Promoter Leads to Multidrug Resistance in the Wheat Pathogen Zymoseptoria tritici

The ascomycete Zymoseptoria tritici is the causal agent of Septoria leaf blotch on wheat. Disease control relies mainly on resistant wheat cultivars and on fungicide applications. The fungus displays a high potential to circumvent both methods. Resistance against all unisite fungicides has been observed over decades. A different type of resistance has emerged among wild populations with multidrug-resistant (MDR) strains. Active fungicide efflux through overexpression of the major facilitator gene MFS1 explains this emerging resistance mechanism. Applying a bulk-progeny sequencing approach, we identified in this study a 519-bp long terminal repeat (LTR) insert in the MFS1 promoter, a relic of a retrotransposon cosegregating with the MDR phenotype. Through gene replacement, we show the insert as a mutation responsible for MFS1 overexpression and the MDR phenotype. Besides this type I insert, we found two different types of promoter inserts in more recent MDR strains. Type I and type II inserts harbor potential transcription factor binding sites, but not the type III insert. Interestingly, all three inserts correspond to repeated elements present at different genomic locations in either IPO323 or other Z. tritici strains. These results underline the plasticity of repeated elements leading to fungicide resistance in Z. tritici and which contribute to its adaptive potential. IMPORTANCE Disease control through fungicides remains an important means to protect crops from fungal diseases and to secure the harvest. Plant-pathogenic fungi, especially Zymoseptoria tritici, have developed resistance against most currently used active ingredients, reducing or abolishing their efficacy. While target site modification is the most common resistance mechanism against single modes of action, active efflux of multiple drugs is an emerging phenomenon in fungal populations reducing additionally fungicides' efficacy in multidrug-resistant strains. We have investigated the mutations responsible for increased drug efflux in Z. tritici field strains. Our study reveals that three different insertions of repeated elements in the same promoter lead to multidrug resistance in Z. tritici. The target gene encodes the membrane transporter MFS1 responsible for drug efflux, with the promoter inserts inducing its overexpression. These results underline the plasticity of repeated elements leading to fungicide resistance in Z. tritici.

Gapless genome assembly of Colletotrichum higginsianum reveals chromosome structure and association of transposable elements with secondary metabolite gene clusters

BACKGROUND

The ascomycete fungus Colletotrichum higginsianum causes anthracnose disease of brassica crops and the model plant Arabidopsis thaliana. Previous versions of the genome sequence were highly fragmented, causing errors in the prediction of protein-coding genes and preventing the analysis of repetitive sequences and genome architecture.

RESULTS

Here, we re-sequenced the genome using single-molecule real-time (SMRT) sequencing technology and, in combination with optical map data, this provided a gapless assembly of all twelve chromosomes except for the ribosomal DNA repeat cluster on chromosome 7. The more accurate gene annotation made possible by this new assembly revealed a large repertoire of secondary metabolism (SM) key genes (89) and putative biosynthetic pathways (77 SM gene clusters). The two mini-chromosomes differed from the ten core chromosomes in being repeat- and AT-rich and gene-poor but were significantly enriched with genes encoding putative secreted effector proteins. Transposable elements (TEs) were found to occupy 7% of the genome by length. Certain TE families showed a statistically significant association with effector genes and SM cluster genes and were transcriptionally active at particular stages of fungal development. All 24 subtelomeres were found to contain one of three highly-conserved repeat elements which, by providing sites for homologous recombination, were probably instrumental in four segmental duplications.

CONCLUSION

The gapless genome of C. higginsianum provides access to repeat-rich regions that were previously poorly assembled, notably the mini-chromosomes and subtelomeres, and allowed prediction of the complete SM gene repertoire. It also provides insights into the potential role of TEs in gene and genome evolution and host adaptation in this asexual pathogen.

Filamentous ascomycete genomes provide insights into Copia retrotransposon diversity in fungi

Background

The relative scarcity of Copia retrotransposons has been recently characterized in metazoans in comparison with the other superfamilies of LTR elements. Furthermore, Copia retrotransposons have often a particular dynamics that results in a highly predominant single clade of elements within a large host taxon, such as the GalEa-like retrotransposons in crustaceans. Taking advantage of the skyrocketing amount of genomic data available for fungi, we carried out the first large-scale comparative genomic analysis of the Copia clades in filamentous ascomycetes.

Results

Screening 30 completely sequenced genomes allowed us to identify more than 2500 Copia copies with conserved LTR, which are distributed in 138 families. Their characterization revealed that fungal Copia diversity is much broader than previously thought with at least 27 clades, 23 of which likely correspond to new ones. While the Copia copy number is low in most species, the two clades GalEa and FunCo1 are widely distributed and highly dominate Copia content as they both account for 80% of the detected sequences.

Conclusions

In Fungi, GalEa retrotransposons are restricted to Pezizomycotina in which they can make up an outstandingly high proportion of the genome (up to 10% in Cenococcum geophilum). At last, we revealed that fungal GalEa elements structurally differ from all other Copia elements with an absence of Primer Binding Site. These elements however harbor a Conserved Hairpin Site which is probably essential for their transposition.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-017-3795-2) contains supplementary material, which is available to authorized users.

DNA sequence homology induces cytosine-to-thymine mutation by a heterochromatin-related pathway in Neurospora

Eukaryotic genomes contain substantial amounts of repetitive DNA organized in the form of constitutive heterochromatin and associated with repressive epigenetic modifications, such as H3K9me3 and C5-cytosine methylation (5mC). In the fungus Neurospora crassa, H3K9me3 and 5mC are catalyzed, respectively, by a conserved SUV39 histone methyltransferase DIM-5 and a DNMT1-like cytosine methyltransferase DIM-2. Here we show that DIM-2 can also mediate Repeat-Induced Point mutation (RIP) of repetitive DNA in N. crassa. We further show that DIM-2-dependent RIP requires DIM-5, HP1, and other known heterochromatin factors, implying the role of a repeat-induced heterochromatin-related process. Our previous findings suggest that the mechanism of repeat recognition for RIP involves direct interactions between homologous double-stranded (ds) DNA segments. We thus now propose that, in somatic cells, homologous dsDNA/dsDNA interactions between a small number of repeat copies can nucleate a transient heterochromatic state, which, on longer repeat arrays, may lead to the formation of constitutive heterochromatin.