Draft Genome Sequence of Neurospora crassa Strain FGSC 73

Whole-genome sequencing and evolutionary analysis of the wild edible mushroom, Morchella eohespera

Morels (Morchella, Ascomycota) are an extremely desired group of edible mushrooms with worldwide distribution. Morchella eohespera is a typical black morel species, belonging to the Elata clade of Morchella species. The biological and genetic studies of this mushroom are rare, largely hindering the studies of molecular breeding and evolutionary aspects. In this study, we performed de novo sequencing and assembly of the M. eohespera strain m200 genome using the third-generation nanopore sequencing platform. The whole-genome size of M. eohespera was 53.81 Mb with a contig N50 of 1.93 Mb, and the GC content was 47.70%. A total of 9,189 protein-coding genes were annotated. Molecular dating showed that M. eohespera differentiated from its relative M. conica at ~19.03 Mya (million years ago) in Burdigalian. Evolutionary analysis showed that 657 gene families were contracted and 244 gene families expanded in M. eohespera versus the related morel species. The non-coding RNA prediction results showed that there were 336 tRNAs, 76 rRNAs, and 45 snRNAs in the M. eohespera genome. Interestingly, there was a high degree of repetition (20.93%) in the M. eohespera genome, and the sizes of long interspersed nuclear elements, short interspersed nuclear elements, and long terminal repeats were 0.83 Mb, 0.009 Mb, and 4.56 Mb, respectively. Additionally, selection pressure analysis identified that a total of 492 genes in the M. eohespera genome have undergone signatures of positive selection. The results of this study provide new insights into the genome evolution of M. eohespera and lay the foundation for in-depth research into the molecular biology of the genus Morchella in the future.

Terabase-Scale Coassembly of a Tropical Soil Microbiome

Petabases of environmental metagenomic data are publicly available, presenting an opportunity to characterize complex environments and discover novel lineages of life. Metagenome coassembly, in which many metagenomic samples from an environment are simultaneously analyzed to infer the underlying genomes' sequences, is an essential tool for achieving this goal. We applied MetaHipMer2, a distributed metagenome assembler that runs on supercomputing clusters, to coassemble 3.4 terabases (Tbp) of metagenome data from a tropical soil in the Luquillo Experimental Forest (LEF), Puerto Rico. The resulting coassembly yielded 39 high-quality (>90% complete, <5% contaminated, with predicted 23S, 16S, and 5S rRNA genes and ≥18 tRNAs) metagenome-assembled genomes (MAGs), including two from the candidate phylum Eremiobacterota. Another 268 medium-quality (≥50% complete, <10% contaminated) MAGs were extracted, including the candidate phyla Dependentiae, Dormibacterota, and Methylomirabilota. In total, 307 medium- or higher-quality MAGs were assigned to 23 phyla, compared to 294 MAGs assigned to nine phyla in the same samples individually assembled. The low-quality (<50% complete, <10% contaminated) MAGs from the coassembly revealed a 49% complete rare biosphere microbe from the candidate phylum FCPU426 among other low-abundance microbes, an 81% complete fungal genome from the phylum Ascomycota, and 30 partial eukaryotic MAGs with ≥10% completeness, possibly representing protist lineages. A total of 22,254 viruses, many of them low abundance, were identified. Estimation of metagenome coverage and diversity indicates that we may have characterized ≥87.5% of the sequence diversity in this humid tropical soil and indicates the value of future terabase-scale sequencing and coassembly of complex environments. IMPORTANCE Petabases of reads are being produced by environmental metagenome sequencing. An essential step in analyzing these data is metagenome assembly, the computational reconstruction of genome sequences from microbial communities. "Coassembly" of metagenomic sequence data, in which multiple samples are assembled together, enables more complete detection of microbial genomes in an environment than "multiassembly," in which samples are assembled individually. To demonstrate the potential for coassembling terabases of metagenome data to drive biological discovery, we applied MetaHipMer2, a distributed metagenome assembler that runs on supercomputing clusters, to coassemble 3.4 Tbp of reads from a humid tropical soil environment. The resulting coassembly, its functional annotation, and analysis are presented here. The coassembly yielded more, and phylogenetically more diverse, microbial, eukaryotic, and viral genomes than the multiassembly of the same data. Our resource may facilitate the discovery of novel microbial biology in tropical soils and demonstrates the value of terabase-scale metagenome sequencing.

A Genetic Engineering Toolbox for the Lignocellulolytic Anaerobic Gut Fungus Neocallimastix frontalis

Anaerobic fungi are powerful platforms for biotechnology that remain unexploited due to a lack of genetic tools. These gut fungi encode the largest number of lignocellulolytic carbohydrate active enzymes (CAZymes) in the fungal kingdom, making them attractive for applications in renewable energy and sustainability. However, efforts to genetically modify anaerobic fungi have remained limited due to inefficient methods for DNA uptake and a lack of characterized genetic parts. We demonstrate that anaerobic fungi are naturally competent for DNA and leverage this to develop a nascent genetic toolbox informed by recently acquired genomes for transient transformation of anaerobic fungi. We validate multiple selectable markers (HygR and Neo), an anaerobic reporter protein (iRFP702), enolase and TEF1A promoters, TEF1A terminator, and a nuclear localization tag for protein compartmentalization. This work establishes novel methods to reliably transform the anaerobic fungus Neocallimastix frontalis, thereby paving the way for strain development and various synthetic biology applications.

First genome-scale insights into the virulence of the snow mold causal fungus Microdochium nivale

Pink snow mold, caused by a phytopathogenic and psychrotolerant fungus, Microdochium nivale, is a severe disease of winter cereals and grasses that predominantly occurs under snow cover or shortly after its melt. Snow mold has significantly progressed during the past decade, often reaching epiphytotic levels in northern countries and resulting in dramatic yield losses. In addition, M. nivale gradually adapts to a warmer climate, spreading to less snowy territories and causing different types of plant diseases throughout the growing period. Despite its great economic importance, M. nivale is poorly investigated; its genome has not been sequenced and its crucial virulence determinants have not been identified or even predicted. In our study, we applied a hybrid assembly based on Oxford Nanopore and Illumina reads to obtain the first genome sequence of M. nivale. 11,973 genes (including 11,789 protein-encoding genes) have been revealed in the genome assembly. To better understand the genetic potential of M. nivale and to obtain a convenient reference for transcriptomic studies on this species, the identified genes were annotated and split into hierarchical three-level functional categories. A file with functionally classified M. nivale genes is presented in our study for general use. M. nivale gene products that best meet the criteria for virulence factors have been identified. The genetic potential to synthesize human-dangerous mycotoxins (fumonisin, ochratoxin B, aflatoxin, and gliotoxin) has been revealed for M. nivale. The transcriptome analysis combined with the assays for extracellular enzymatic activities (conventional virulence factors of many phytopathogens) was carried out to assess the effect of host plant (rye) metabolites on the M. nivale phenotype. In addition to disclosing plant-metabolite-upregulated M. nivale functional gene groups (including those related to host plant protein destruction and amino acid metabolism, xenobiotic detoxication (including phytoalexins benzoxazinoids), cellulose destruction (cellulose monooxygenases), iron transport, etc.), the performed analysis pointed to a crucial role of host plant lipid destruction and fungal lipid metabolism modulation in plant-M. nivale interactions.

Cytospora paraplurivora sp. nov. isolated from orchards with fruit tree decline syndrome in Ontario, Canada

A new species of Cytospora was isolated from cankered wood of Prunus spp. during a survey of orchards exhibiting symptoms of fruit tree decline syndrome in southern Ontario, Canada. We found isolates that are morphologically similar to species in the Cytosporaceae family, which is characterized by single or labyrinthine locules, filamentous conidiophores or clavate to elongate obovoid asci and allantoid, hyaline conidia. Multi-gene phylogenetic analysis of ITS, LSU, act and tef1- α showed that the isolates form a distinct clade, sister to Cytospora plurivora. Morphologically, our isolates showed differences in the length of conidia and culture characteristics compared to C. plurivora, suggesting the establishment of a new species. The species is described as Cytospora paraplurivora sp. nov. and placed in the family Cytosporaceae of Diaporthales. Additionally, we sequenced, assembled and characterized the genome of the representative isolate for this new species. The phylogenomic analysis confirms the species order and family level classification. C. paraplurivora sp. nov. has the potential to severely affect stone fruits production, causing cankers and dieback in stressed trees, and eventually leads to tree decline. Pathogenicity tests show that the species is pathogenic to Prunus persica var. persica.

FunOrder: A robust and semi-automated method for the identification of essential biosynthetic genes through computational molecular co-evolution

Secondary metabolites (SMs) are a vast group of compounds with different structures and properties that have been utilized as drugs, food additives, dyes, and as monomers for novel plastics. In many cases, the biosynthesis of SMs is catalysed by enzymes whose corresponding genes are co-localized in the genome in biosynthetic gene clusters (BGCs). Notably, BGCs may contain so-called gap genes, that are not involved in the biosynthesis of the SM. Current genome mining tools can identify BGCs, but they have problems with distinguishing essential genes from gap genes. This can and must be done by expensive, laborious, and time-consuming comparative genomic approaches or transcriptome analyses. In this study, we developed a method that allows semi-automated identification of essential genes in a BGC based on co-evolution analysis. To this end, the protein sequences of a BGC are blasted against a suitable proteome database. For each protein, a phylogenetic tree is created. The trees are compared by treeKO to detect co-evolution. The results of this comparison are visualized in different output formats, which are compared visually. Our results suggest that co-evolution is commonly occurring within BGCs, albeit not all, and that especially those genes that encode for enzymes of the biosynthetic pathway are co-evolutionary linked and can be identified with FunOrder. In light of the growing number of genomic data available, this will contribute to the studies of BGCs in native hosts and facilitate heterologous expression in other organisms with the aim of the discovery of novel SMs.

The discovery and description of novel fungal secondary metabolites promises novel antibiotics, pharmaceuticals, and other useful compounds. A way to identify novel secondary metabolites is to express the corresponding genes in a suitable expression host. Consequently, a detailed knowledge or an accurate prediction of these genes is necessary. In fungi, the genes are co-localized in so-called biosynthetic gene clusters. Notably, the clusters may also contain genes that are not necessary for the biosynthesis of the secondary metabolites, so-called gap genes. We developed a method to detect co-evolved genes within the clusters and demonstrated that essential genes are co-evolving and can thus be differentiated from the gap genes. This adds an additional layer of information, which can support researchers with their decisions on which genes to study and express for the discovery of novel secondary metabolites.

Comparative Transcriptome and Endophytic Bacterial Community Analysis of Morchella conica SH

The precious rare edible fungus Morchella conica is popular worldwide for its rich nutrition, savory flavor, and varieties of bioactive components. Due to its high commercial, nutritional, and medicinal value, it has always been a hot spot. However, the molecular mechanism and endophytic bacterial communities in M. conica were poorly understood. In this study, we sequenced, assembled, and analyzed the genome of M. conica SH. Transcriptome analysis reveals significant differences between the mycelia and fruiting body. As shown in this study, 1,329 and 2,796 genes were specifically expressed in the mycelia and fruiting body, respectively. The Gene Ontology (GO) enrichment showed that RNA polymerase II transcription activity-related genes were enriched in the mycelium-specific gene cluster, and nucleotide binding-related genes were enriched in the fruiting body-specific gene cluster. Further analysis of differentially expressed genes in different development stages resulted in finding two groups with distinct expression patterns. Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment displays that glycan degradation and ABC transporters were enriched in the group 1 with low expressed level in the mycelia, while taurine and hypotaurine metabolismand tyrosine metabolism-related genes were significantly enriched in the group 2 with high expressed level in mycelia. Moreover, a dynamic shift of bacterial communities in the developing fruiting body was detected by 16S rRNA sequencing, and co-expression analysis suggested that bacterial communities might play an important role in regulating gene expression. Taken together, our study provided a better understanding of the molecular biology of M. conica SH and direction for future research on artificial cultivation.

Distribution of methionine sulfoxide reductases in fungi and conservation of the free-methionine-R-sulfoxide reductase in multicellular eukaryotes

Methionine, either as a free amino acid or included in proteins, can be oxidized into methionine sulfoxide (MetO), which exists as R and S diastereomers. Almost all characterized organisms possess thiol-oxidoreductases named methionine sulfoxide reductase (Msr) enzymes to reduce MetO back to Met. MsrA and MsrB reduce the S and R diastereomers of MetO, respectively, with strict stereospecificity and are found in almost all organisms. Another type of thiol-oxidoreductase, the free-methionine-R-sulfoxide reductase (fRMsr), identified so far in prokaryotes and a few unicellular eukaryotes, reduces the R MetO diastereomer of the free amino acid. Moreover, some bacteria possess molybdenum-containing enzymes that reduce MetO, either in the free or protein-bound forms. All these Msrs play important roles in the protection of organisms against oxidative stress. Fungi are heterotrophic eukaryotes that colonize all niches on Earth and play fundamental functions, in organic matter recycling, as symbionts, or as pathogens of numerous organisms. However, our knowledge on fungal Msrs is still limited. Here, we performed a survey of msr genes in almost 700 genomes across the fungal kingdom. We show that most fungi possess one gene coding for each type of methionine sulfoxide reductase: MsrA, MsrB, and fRMsr. However, several fungi living in anaerobic environments or as obligate intracellular parasites were devoid of msr genes. Sequence inspection and phylogenetic analyses allowed us to identify non-canonical sequences with potentially novel enzymatic properties. Finaly, we identified several ocurences of msr horizontal gene transfer from bacteria to fungi.

Evolution and functional characterization of pectate lyase PEL12, a member of a highly expanded Clonostachys rosea polysaccharide lyase 1 family

Background

Pectin is one of the major and most complex plant cell wall components that needs to be overcome by microorganisms as part of their strategies for plant invasion or nutrition. Microbial pectinolytic enzymes therefore play a significant role for plant-associated microorganisms and for the decomposition and recycling of plant organic matter. Recently, comparative studies revealed significant gene copy number expansion of the polysaccharide lyase 1 (PL1) pectin/pectate lyase gene family in the Clonostachys rosea genome, while only low numbers were found in Trichoderma species. Both of these fungal genera are widely known for their ability to parasitize and kill other fungi (mycoparasitism) and certain species are thus used for biocontrol of plant pathogenic fungi.

Results

In order to understand the role of the high number of pectin degrading enzymes in Clonostachys, we studied diversity and evolution of the PL1 gene family in C. rosea compared with other Sordariomycetes with varying nutritional life styles. Out of 17 members of C. rosea PL1, we could only detect two to be secreted at acidic pH. One of them, the pectate lyase pel12 gene was found to be strongly induced by pectin and, to a lower degree, by polygalacturonic acid. Heterologous expression of the PEL12 in a PL1-free background of T. reesei revealed direct enzymatic involvement of this protein in utilization of pectin at pH 5 without a requirement for Ca²⁺. The mutants showed increased utilization of pectin compounds, but did not increase biocontrol ability in detached leaf assay against the plant pathogen Botrytis cinerea compared to the wild type.

Conclusions

In this study, we aimed to gain insight into diversity and evolution of the PL1 gene family in C. rosea and other Sordariomycete species in relation to their nutritional modes. We show that C. rosea PL1 expansion does not correlate with its mycoparasitic nutritional mode and resembles those of strong plant pathogenic fungi. We further investigated regulation, specificity and function of the C. rosea PEL12 and show that this enzyme is directly involved in degradation of pectin and pectin-related compounds, but not in C. rosea biocontrol.

Electronic supplementary material

The online version of this article (10.1186/s12866-018-1310-9) contains supplementary material, which is available to authorized users.

Genome analysis of the ubiquitous boxwood pathogen Pseudonectria foliicola

Boxwood (Buxus spp.) are broad-leaved, evergreen landscape plants valued for their longevity and ornamental qualities. Volutella leaf and stem blight, caused by the ascomycete fungi Pseudonectria foliicola and P. buxi, is one of the major diseases affecting the health and ornamental qualities of boxwood. Although this disease is less severe than boxwood blight caused by Calonectria pseudonaviculata and C. henricotiae, its widespread occurrence and disfiguring symptoms have caused substantial economic losses to the ornamental industry. In this study, we sequenced the genome of P. foliicola isolate ATCC13545 using Illumina technology and compared it to other publicly available fungal pathogen genomes to better understand the biology of this organism. A de novo assembly estimated the genome size of P. foliicola at 28.7 Mb (425 contigs; N50 = 184,987 bp; avg. coverage 188×), with just 9,272 protein-coding genes. To our knowledge, P. foliicola has the smallest known genome within the Nectriaceae. Consistent with the small size of the genome, the secretome, CAzyme and secondary metabolite profiles of this fungus are reduced relative to two other surveyed Nectriaceae fungal genomes: Dactylonectria macrodidyma JAC15-245 and Fusarium graminearum Ph-1. Interestingly, a large cohort of genes associated with reduced virulence and loss of pathogenicity was identified from the P. foliicola dataset. These data are consistent with the latest observations by plant pathologists that P. buxi and most likely P. foliicola, are opportunistic, latent pathogens that prey upon weak and stressed boxwood plants.

Fungal Genomes and Genotyping

The availability of complete fungal genomes is expanding rapidly and is offering an extensive and accurate view of this "kingdom." The scientific milestone of free access to more than 1000 fungal genomes of different species was reached, and new and stimulating projects have meanwhile been released. The "1000 Fungal Genomes Project" represents one of the largest sequencing initiative regarding fungal organisms trying to fill some gaps on fungal genomics. Presently, there are 329 fungal families with at least one representative genome sequenced, but there is still a large number of fungal families without a single sequenced genome. In addition, additional sequencing projects helped to understand the genetic diversity within some fungal species. The availability of multiple genomes per species allows to support taxonomic organization, brings new insights for fungal evolution in short-time scales, clarifies geographical and dispersion patterns, elucidates outbreaks and transmission routes, among other objectives. Genotyping methodologies analyze only a small fraction of an individual's genome but facilitate the comparison of hundreds or thousands of isolates in a small fraction of the time and at low cost. The integration of whole genome strategies and improved genotyping panels targeting specific and relevant SNPs and/or repeated regions can represent fast and practical strategies for studying local, regional, and global epidemiology of fungi.

Diverse data supports the transition of filamentous fungal model organisms into the post-genomics era

Filamentous fungi have been important as model organisms since the beginning of modern biological inquiry and have benefitted from open data since the earliest genetic maps were shared. From early origins in simple Mendelian genetics of mating types, parasexual genetics of colony colour, and the foundational demonstration of the segregation of a nutritional requirement, the contribution of research systems utilising filamentous fungi has spanned the biochemical genetics era, through the molecular genetics era, and now are at the very foundation of diverse omics approaches to research and development. Fungal model organisms have come from most major taxonomic groups although Ascomycete filamentous fungi have seen the most major sustained effort. In addition to the published material about filamentous fungi, shared molecular tools have found application in every area of fungal biology. Similarly, shared data has contributed to the success of model systems. The scale of data supporting research with filamentous fungi has grown by 10 to 12 orders of magnitude. From genetic to molecular maps, expression databases, and finally genome resources, the open and collaborative nature of the research communities has assured that the rising tide of data has lifted all of the research systems together.

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

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

Genetics of mating in members of the Chaetomiaceae as revealed by experimental and genomic characterization of reproduction in Myceliophthora heterothallica

Members of the Chaetomiaceae are among the most studied fungi in industry and among the most reported in investigations of biomass degradation in both natural and laboratory settings. The family is recognized for production of carbohydrate-active enzymes and antibiotics. Thermophilic species are of special interest for their abilities to produce thermally stable enzymes and to be grown under conditions that are unsuitable for potential contaminant microorganisms. Such interests led to the recent acquisition of genome sequences from several members of the family, including thermophilic species, several of which are reported here for the first time. To date, however, thermophilic fungi in industry have served primarily as parts reservoirs and there has been no good genetic model for species in the family Chaetomiaceae or for thermophiles in general. We report here on the reproductive biology of the thermophile Myceliophthora heterothallica, which is heterothallic, unlike most described species in the family. We confirmed heterothallism genetically by following the segregation of mating type idiomorphs and other markers. We have expanded the number of known sexually-compatible individuals from the original isolates from Indiana and Germany to include several isolates from New Mexico. An interesting aspect of development in M. heterothallica is that ascocarp formation is optimal at approximately 30 °C, whereas vegetative growth is optimal at 45 °C. Genome sequences obtained from several strains, including isolates of each mating type, revealed mating-type regions whose genes are organized similarly to those of other members of the Sordariales, except for the presence of a truncated version of the mat A-1 (MAT1-1-1) gene in mating-type a (MAT1-2) strains. In M. heterothallica and other Chaetomiaceae, mating-type A (MAT1-1) strains have the full-length version of mat A-1 that is typical of mating-type A strains of diverse Ascomycota, whereas a strains have only the truncated version. This truncated mat A-1 has an intact open reading frame and a derived start codon that is not present in mat A-1 from A strains. The predicted protein contains a region that is conserved across diverse mat A-1 genes, but it lacks the major alpha1 domain, which characterizes proteins in this family and is known to be required for fertility in A strains from other Ascomycota. Finally, we have used genes from M. heterothallica to probe for mating genes in other homothallic and heterothallic members of the Chaetomiaceae. The majority of homothallic species examined have a typical mat A-1,2,3 (MAT1-1-1,2,3) region in addition to an unlinked mat a-1 (MAT1-2-1) gene, reflecting one type of homothallism commonly observed in diverse Ascomycota.