Search

Thermophily is a trait scattered across the fungal tree of life, with its highest prevalence within three fungal families (Chaetomiaceae, Thermoascaceae, and Trichocomaceae), as well as some members of the phylum Mucoromycota. We examined 37 thermophilic and thermotolerant species and 42 mesophilic species for this study and identified thermophily as the ancestral state of all three prominent families of thermophilic fungi. Thermophilic fungal genomes were found to encode various thermostable enzymes, including carbohydrate-active enzymes such as endoxylanases, which are useful for many industrial applications. At the same time, the overall gene counts, especially in gene families responsible for microbial defense such as secondary metabolism, are reduced in thermophiles compared to mesophiles. We also found a reduction in the core genome size of thermophiles in both the Chaetomiaceae family and the Eurotiomycetes class. The Gene Ontology terms lost in thermophilic fungi include primary metabolism, transporters, UV response, and O-methyltransferases. Comparative genomics analysis also revealed higher GC content in the third base of codons (GC3) and a lower effective number of codons in fungal thermophiles than in both thermotolerant and mesophilic fungi. Furthermore, using the Support Vector Machine classifier, we identified several Pfam domains capable of discriminating between genomes of thermophiles and mesophiles with 94% accuracy. Using AlphaFold2 to predict protein structures of endoxylanases (GH10), we built a similarity network based on the structures. We found that the number of disulfide bonds appears important for protein structure, and the network clusters based on protein structures correlate with the optimal activity temperature. Thus, comparative genomics offers new insights into the biology, adaptation, and evolutionary history of thermophilic fungi while providing a parts list for bioengineering applications.

BACKGROUND: The colonization of land and the diversification of terrestrial plants is intimately linked to the evolutionary history of their symbiotic fungal partners. Extant representatives of these fungal lineages include mutualistic plant symbionts, the arbuscular mycorrhizal (AM) fungi in Glomeromycota and fine root endophytes in Endogonales (Mucoromycota), as well as fungi with saprotrophic, pathogenic and endophytic lifestyles. These fungal groups separate into three monophyletic lineages but their evolutionary relationships remain enigmatic confounding ancestral reconstructions. Their taxonomic ranks are currently fluid. RESULTS: In this study, we recognize these three monophyletic linages as phyla, and use a balanced taxon sampling and broad taxonomic representation for phylogenomic analysis that rejects a hard polytomy and resolves Glomeromycota as sister to a clade composed of Mucoromycota and Mortierellomycota. Low copy numbers of genes associated with plant cell wall degradation could not be assigned to the transition to a plant symbiotic lifestyle but appears to be an ancestral phylogenetic signal. Both plant symbiotic lineages, Glomeromycota and Endogonales, lack numerous thiamine metabolism genes but the lack of fatty acid synthesis genes is specific to AM fungi. Many genes previously thought to be missing specifically in Glomeromycota are either missing in all analyzed phyla, or in some cases, are actually present in some of the analyzed AM fungal lineages, e.g. the high affinity phosphorus transporter Pho89. CONCLUSION: Based on a broad taxon sampling of fungal genomes we present a well-supported phylogeny for AM fungi and their sister lineages. We show that among these lineages, two independent evolutionary transitions to mutualistic plant symbiosis happened in a genomic background profoundly different from that known from the emergence of ectomycorrhizal fungi in Dikarya. These results call for further reevaluation of genomic signatures associated with plant symbiosis.

The model moss species Physcomitrium patens has long been used for studying divergence of land plants spanning from bryophytes to angiosperms. In addition to its phylogenetic relationships, the limited number of differential tissues, and comparable morphology to the earliest embryophytes provide a system to represent basic plant architecture. Based on plant-fungal interactions today, it is hypothesized these kingdoms have a long-standing relationship, predating plant terrestrialization. Mortierellaceae have origins diverging from other land fungi paralleling bryophyte divergence, are related to arbuscular mycorrhizal fungi but are free-living, observed to interact with plants, and can be found in moss microbiomes globally. Due to their parallel origins, we assess here how two Mortierellaceae species, Linnemannia elongata and Benniella erionia, interact with P. patens in coculture. We also assess how Mollicute-related or Burkholderia-related endobacterial symbionts (MRE or BRE) of these fungi impact plant response. Coculture interactions are investigated through high-throughput phenomics, microscopy, RNA-sequencing, differential expression profiling, gene ontology enrichment, and comparisons among 99 other P. patens transcriptomic studies. Here we present new high-throughput approaches for measuring P. patens growth, identify novel expression of over 800 genes that are not expressed on traditional agar media, identify subtle interactions between P. patens and Mortierellaceae, and observe changes to plant-fungal interactions dependent on whether MRE or BRE are present. Our study provides insights into how plants and fungal partners may have interacted based on their communications observed today as well as identifying L. elongata and B. erionia as modern fungal endophytes with P. patens.

Aspergillus terreus has attracted interest due to its application in industrial biotechnology, particularly for the production of itaconic acid and bioactive secondary metabolites. As related species also seem to possess a prosperous secondary metabolism, they are of high interest for genome mining and exploitation. Here, we present draft genome sequences for six species from Aspergillus section Terrei and one species from Aspergillus section Nidulantes. Whole-genome phylogeny confirmed that section Terrei is monophyletic. Genome analyses identified between 70 and 108 key secondary metabolism genes in each of the genomes of section Terrei, the highest rate found in the genus Aspergillus so far. The respective enzymes fall into 167 distinct families with most of them corresponding to potentially unique compounds or compound families. Moreover, 53% of the families were only found in a single species, which supports the suitability of species from section Terrei for further genome mining. Intriguingly, this analysis, combined with heterologous gene expression and metabolite identification, suggested that species from section Terrei use a strategy for UV protection different to other species from the genus Aspergillus. Section Terrei contains a complete plant polysaccharide degrading potential and an even higher cellulolytic potential than other Aspergilli, possibly facilitating additional applications for these species in biotechnology.

DNA N6-adenine methylation (6mA) has recently gained importance as an epigenetic modification in eukaryotes. Its function in lineages with high levels, such as early-diverging fungi (EDF), is of particular interest. Here, we investigated the biological significance and evolutionary implications of 6mA in EDF, which exhibit divergent evolutionary patterns in 6mA usage. The analysis of two Mucorales species displaying extreme 6mA usage reveals that species with high 6mA levels show symmetric methylation enriched in highly expressed genes. In contrast, species with low 6mA levels show mostly asymmetric 6mA. Interestingly, transcriptomic regulation throughout development and in response to environmental cues is associated with changes in the 6mA landscape. Furthermore, we identify an EDF-specific methyltransferase, likely originated from endosymbiotic bacteria, as responsible for asymmetric methylation, while an MTA-70 methylation complex performs symmetric methylation. The distinct phenotypes observed in the corresponding mutants reinforced the critical role of both types of 6mA in EDF.

The Lipomyces clade contains oleaginous yeast species with advantageous metabolic features for biochemical and biofuel production. Limited knowledge about the metabolic networks of the species and limited tools for genetic engineering have led to a relatively small amount of research on the microbes. Here, a genome-scale metabolic model (GSM) of Lipomyces starkeyi NRRL Y-11557 was built using orthologous protein mappings to model yeast species. Phenotypic growth assays were used to validate the GSM (66% accuracy) and indicated that NRRL Y-11557 utilized diverse carbohydrates but had more limited catabolism of organic acids. The final GSM contained 2,193 reactions, 1,909 metabolites, and 996 genes and was thus named iLst996. The model contained 96 of the annotated carbohydrate-active enzymes. iLst996 predicted a flux distribution in line with oleaginous yeast measurements and was utilized to predict theoretical lipid yields. Twenty-five other yeasts in the Lipomyces clade were then genome sequenced and annotated. Sixteen of the Lipomyces species had orthologs for more than 97% of the iLst996 genes, demonstrating the usefulness of iLst996 as a broad GSM for Lipomyces metabolism. Pathways that diverged from iLst996 mainly revolved around alternate carbon metabolism, with ortholog groups excluding NRRL Y-11557 annotated to be involved in transport, glycerolipid, and starch metabolism, among others. Overall, this study provides a useful modeling tool and data for analyzing and understanding Lipomyces species metabolism and will assist further engineering efforts in Lipomyces.

BackgroundCost-effective production of biofuels from lignocellulose requires the fermentation of D-xylose. Many yeast species within and closely related to the genera Spathaspora and Scheffersomyces (both of the order Serinales) natively assimilate and ferment xylose. Other species consume xylose inefficiently, leading to extracellular accumulation of xylitol. Xylitol excretion is thought to be due to the different cofactor requirements of the first two steps of xylose metabolism. Xylose reductase (XR) generally uses NADPH to reduce xylose to xylitol, while xylitol dehydrogenase (XDH) generally uses NAD+ to oxidize xylitol to xylulose, creating an imbalanced redox pathway. This imbalance is thought to be particularly consequential in hypoxic or anoxic environments.ResultsWe screened the growth of xylose-fermenting yeast species in high and moderate aeration and identified both ethanol producers and xylitol producers. Selected species were further characterized for their XR and XDH cofactor preferences by enzyme assays and gene expression patterns by RNA-Seq. Our data revealed that xylose metabolism is more redox balanced in some species, but it is strongly affected by oxygen levels. Under high aeration, most species switched from ethanol production to xylitol accumulation, despite the availability of ample oxygen to accept electrons from NADH. This switch was followed by decreases in enzyme activity and the expression of genes related to xylose metabolism, suggesting that bottlenecks in xylose fermentation are not always due to cofactor preferences. Finally, we expressed XYL genes from multiple Scheffersomyces species in a strain of Saccharomyces cerevisiae. Recombinant S. cerevisiae expressing XYL1 from Scheffersomyces xylosifermentans, which encodes an XR without a cofactor preference, showed improved anaerobic growth on xylose as the primary carbon source compared to S. cerevisiae strain expressing XYL genes from Scheffersomyces stipitis.ConclusionCollectively, our data do not support the hypothesis that xylitol accumulation occurs primarily due to differences in cofactor preferences between xylose reductase and xylitol dehydrogenase; instead, gene expression plays a major role in response to oxygen levels. We have also identified the yeast Sc. xylosifermentans as a potential source for genes that can be engineered into S. cerevisiae to improve xylose fermentation and biofuel production.

Intimate associations between fungi and intracellular bacterial endosymbionts are becoming increasingly well understood. Phylogenetic analyses demonstrate that bacterial endosymbionts of Mucoromycota fungi are related either to free-living Burkholderia or Mollicutes species. The so-called Burkholderia-related endosymbionts or BRE comprise Mycoavidus, Mycetohabitans and Candidatus Glomeribacter gigasporarum. These endosymbionts are marked by genome contraction thought to be associated with intracellular selection. However, the conclusions drawn thus far are based on a very small subset of endosymbiont genomes, and the mechanisms leading to genome streamlining are not well understood. The purpose of this study was to better understand how intracellular existence shapes Mycoavidus and BRE functionally at the genome level. To this end we generated and analyzed 14 novel draft genomes for Mycoavidus living within the hyphae of Mortierellomycotina fungi. We found that our novel Mycoavidus genomes were significantly reduced compared to free-living Burkholderiales relatives. Using a genome-scale phylogenetic approach including the novel and available existing genomes of Mycoavidus, we show that the genus is an assemblage composed of two independently derived lineages including three well supported clades of Mycoavidus. Using a comparative genomic approach, we shed light on the functional implications of genome reduction, documenting shared and unique gene loss patterns between the three Mycoavidus clades. We found that many endosymbiont isolates demonstrate patterns of vertical transmission and host-specificity, but others are present in phylogenetically disparate hosts. We discuss how reductive evolution and host specificity reflect convergent adaptation to the intrahyphal selective landscape, and commonalities of eukaryotic endosymbiont genome evolution.

The order Sordariales is taxonomically diverse, and harbours many species with different lifestyles and large economic importance. Despite its importance, a robust genome-scale phylogeny, and associated comparative genomic analysis of the order is lacking. In this study, we examined whole-genome data from 99 Sordariales, including 52 newly sequenced genomes, and seven outgroup taxa. We inferred a comprehensive phylogeny that resolved several contentious relationships amongst families in the order, and cleared-up intrafamily relationships within the Podosporaceae. Extensive comparative genomics showed that genomes from the three largest families in the dataset (Chaetomiaceae, Podosporaceae and Sordariaceae) differ greatly in GC content, genome size, gene number, repeat percentage, evolutionary rate, and genome content affected by repeat-induced point mutations (RIP). All genomic traits showed phylogenetic signal, and ancestral state reconstruction revealed that the variation of the properties stems primarily from within-family evolution. Together, the results provide a thorough framework for understanding genome evolution in this important group of fungi.

Diverse members of early-diverging Mucoromycota, including mycorrhizal taxa and soil-associated Mortierellaceae, are known to harbor Mollicutes-related endobacteria (MRE). It has been hypothesized that MRE were acquired by a common ancestor and transmitted vertically. Alternatively, MRE endosymbionts could have invaded after the divergence of Mucoromycota lineages and subsequently spread to new hosts horizontally. To better understand the evolutionary history of MRE symbionts, we generated and analyzed four complete MRE genomes from two Mortierellaceae genera: Linnemannia (MRE-L) and Benniella (MRE-B). These genomes include the smallest known of fungal endosymbionts and showed signals of a tight relationship with hosts including a reduced functional capacity and genes transferred from fungal hosts to MRE. Phylogenetic reconstruction including nine MRE from mycorrhizal fungi revealed that MRE-B genomes are more closely related to MRE from Glomeromycotina than MRE-L from the same host family. We posit that reductions in genome size, GC content, pseudogene content, and repeat content in MRE-L may reflect a longer-term relationship with their fungal hosts. These data indicate Linnemannia and Benniella MRE were likely acquired independently after their fungal hosts diverged from a common ancestor. This work expands upon foundational knowledge on minimal genomes and provides insights into the evolution of bacterial endosymbionts.

Fungi have evolved over millions of years and their species diversity is predicted to be the second largest on the earth. Fungi have cross-kingdom interactions with many organisms that have mutually shaped their evolutionary trajectories. Zygomycete fungi hold a pivotal position in the fungal tree of life and provide important perspectives on the early evolution of fungi from aquatic to terrestrial environments. Phylogenomic analyses have found that zygomycete fungi diversified into two separate clades, the Mucoromycota which are frequently associated with plants and Zoopagomycota that are commonly animal-associated fungi. Genetic elements that contributed to the fitness and divergence of these lineages may have been shaped by the varied interactions these fungi have had with plants, animals, bacteria, and other microbes. To investigate this, we performed comparative genomic analyses of the two clades of zygomycetes in the context of Kingdom Fungi, benefiting from our generation of a new collection of zygomycete genomes, including nine produced for this study. We identified lineage-specific genomic content that may contribute to the disparate biology observed in these zygomycetes. Our findings include the discovery of undescribed diversity in CotH, a Mucormycosis pathogenicity factor, which was found in a broad set of zygomycetes. Reconciliation analysis identified multiple duplication events and an expansion of CotH copies throughout the Mucoromycotina, Mortierellomycotina, Neocallimastigomycota, and Basidiobolus lineages. A kingdom-level phylogenomic analysis also identified new evolutionary relationships within the subphyla of Mucoromycota and Zoopagomycota, including supporting the sister-clade relationship between Glomeromycotina and Mortierellomycotina and the placement of Basidiobolus as sister to other Zoopagomycota lineages.

Improved sequencing technologies have profoundly altered global views of fungal diversity and evolution. High-throughput sequencing methods are critical for studying fungi due to the cryptic, symbiotic nature of many species, particularly those that are difficult to culture. However, the low coverage genome sequencing (LCGS) approach to phylogenomic inference has not been widely applied to fungi. Here we analyzed 171 Kickxellomycotina fungi using LCGS methods to obtain hundreds of marker genes for robust phylogenomic reconstruction. Additionally, we mined our LCGS data for a set of nine rDNA and protein coding genes to enable analyses across species for which no LCGS data were obtained. The main goals of this study were to: 1) evaluate the quality and utility of LCGS data for both phylogenetic reconstruction and functional annotation, 2) test relationships among clades of Kickxellomycotina, and 3) perform comparative functional analyses between clades to gain insight into putative trophic modes. In opposition to previous studies, our nine-gene analyses support two clades of arthropod gut dwelling species and suggest a possible single evolutionary event leading to this symbiotic lifestyle. Furthermore, we resolve the mycoparasitic Dimargaritales as the earliest diverging clade in the subphylum and find four major clades of Coemansia species. Finally, functional analyses illustrate clear variation in predicted carbohydrate active enzymes and secondary metabolites (SM) based on ecology, that is biotroph versus saprotroph. Saprotrophic Kickxellales broadly lack many known pectinase families compared with saprotrophic Mucoromycota and are depauperate for SM but have similar numbers of predicted chitinases as mycoparasitic.

The first genome sequenced of a eukaryotic organism was for Saccharomyces cerevisiae, as reported in 1996, but it was more than 10 years before any of the zygomycete fungi, which are the early-diverging terrestrial fungi currently placed in the phyla Mucoromycota and Zoopagomycota, were sequenced. The genome for Rhizopus delemar was completed in 2008; currently, more than 1000 zygomycete genomes have been sequenced. Genomic data from these early-diverging terrestrial fungi revealed deep phylogenetic separation of the two major clades-primarily plant-associated saprotrophic and mycorrhizal Mucoromycota versus the primarily mycoparasitic or animal-associated parasites and commensals in the Zoopagomycota. Genomic studies provide many valuable insights into how these fungi evolved in response to the challenges of living on land, including adaptations to sensing light and gravity, development of hyphal growth, and co-existence with the first terrestrial plants. Genome sequence data have facilitated studies of genome architecture, including a history of genome duplications and horizontal gene transfer events, distribution and organization of mating type loci, rDNA genes and transposable elements, methylation processes, and genes useful for various industrial applications. Pathogenicity genes and specialized secondary metabolites have also been detected in soil saprobes and pathogenic fungi. Novel endosymbiotic bacteria and viruses have been discovered during several zygomycete genome projects. Overall, genomic information has helped to resolve a plethora of research questions, from the placement of zygomycetes on the evolutionary tree of life and in natural ecosystems, to the applied biotechnological and medical questions.

The global regulator LaeA controls secondary metabolism in diverse Aspergillus species. Here we explored its role in regulation of itaconic acid production in Aspergillus pseudoterreus. To understand its role in regulating metabolism, we deleted and overexpressed laeA, and assessed the transcriptome, proteome, and secreted metabolome prior to and during initiation of phosphate limitation induced itaconic acid production. We found that secondary metabolite clusters, including the itaconic acid biosynthetic gene cluster, are regulated by laeA and that laeA is required for high yield production of itaconic acid. Overexpression of LaeA improves itaconic acid yield at the expense of biomass by increasing the expression of key biosynthetic pathway enzymes and attenuating the expression of genes involved in phosphate acquisition and scavenging. Increased yield was observed in optimized conditions as well as conditions containing excess nutrients that may be present in inexpensive sugar containing feedstocks such as excess phosphate or complex nutrient sources. This suggests that global regulators of metabolism may be useful targets for engineering metabolic flux that is robust to environmental heterogeneity.

Most of the described species in kingdom Fungi are contained in two phyla, the Ascomycota and the Basidiomycota (subkingdom Dikarya). As a result, our understanding of the biology of the kingdom is heavily influenced by traits observed in Dikarya, such as aerial spore dispersal and life cycles dominated by mitosis of haploid nuclei. We now appreciate that Fungi comprises numerous phylum-level lineages in addition to those of Dikarya, but the phylogeny and genetic characteristics of most of these lineages are poorly understood due to limited genome sampling. Here, we addressed major evolutionary trends in the non-Dikarya fungi by phylogenomic analysis of 69 newly generated draft genome sequences of the zoosporic (flagellated) lineages of true fungi. Our phylogeny indicated five lineages of zoosporic fungi and placed Blastocladiomycota, which has an alternation of haploid and diploid generations, as branching closer to the Dikarya than to the Chytridiomyceta. Our estimates of heterozygosity based on genome sequence data indicate that the zoosporic lineages plus the Zoopagomycota are frequently characterized by diploid-dominant life cycles. We mapped additional traits, such as ancestral cell-cycle regulators, cell-membrane- and cell-wall-associated genes, and the use of the amino acid selenocysteine on the phylogeny and found that these ancestral traits that are shared with Metazoa have been subject to extensive parallel loss across zoosporic lineages. Together, our results indicate a gradual transition in the genetics and cell biology of fungi from their ancestor and caution against assuming that traits measured in Dikarya are typical of other fungal lineages.

Forest soil microbiomes have crucial roles in carbon storage, biogeochemical cycling and rhizosphere processes. Wildfire season length, and the frequency and size of severe fires have increased owing to climate change. Fires affect ecosystem recovery and modify soil microbiomes and microbially mediated biogeochemical processes. To study wildfire-dependent changes in soil microbiomes, we characterized functional shifts in the soil microbiota (bacteria, fungi and viruses) across burn severity gradients (low, moderate and high severity) 1 yr post fire in coniferous forests in Colorado and Wyoming, USA. We found severity-dependent increases of Actinobacteria encoding genes for heat resistance, fast growth, and pyrogenic carbon utilization that might enhance post-fire survival. We report that increased burn severity led to the loss of ectomycorrhizal fungi and less tolerant microbial taxa. Viruses remained active in post-fire soils and probably influenced carbon cycling and biogeochemistry via turnover of biomass and ecosystem-relevant auxiliary metabolic genes. Our genome-resolved analyses link post-fire soil microbial taxonomy to functions and reveal the complexity of post-fire soil microbiome activity.

Fungi survive in diverse ecological niches by secreting proteins and other molecules into the environment to acquire food and interact with various biotic and abiotic stressors. Fungal secretome content is, therefore, believed to be tightly linked to fungal ecologies. We sampled 132 genomes from the early-diverging terrestrial fungal lineage zygomycetes (Mucoromycota and Zoopagomycota) and characterized their secretome composition. Our analyses revealed that phylogeny played an important role in shaping the secretome composition of zygomycete fungi with trophic mode contributing a smaller amount. Reconstruction of the evolution of secreted digestive enzymes revealed lineage-specific expansions, indicating that Mucoromycota and Zoopagomycota followed different trajectories early in their evolutionary history. We identified the presence of multiple pathogenicity-related proteins in the lineages known as saprotrophs, suggesting that either the ecologies of these fungi are incompletely known, and/or that these pathogenicity-related proteins have important functions associated with saprotrophic ecologies, both of which invite further investigation.

Profiling the host-mycobiota interactions in healthy vs. diseased forest ecosystems helps understand the dynamics of understudied yet increasingly important threats to forest health that are emerging due to climate change. We analyzed the structural and functional changes of the mycobiota and the responses of Pinus contorta in the Lophodermella needle cast pathosystem through metabarcoding and metatranscriptomics. When needles transitioned from asymptomatic to symptomatic, dysbiosis of the mycobiota occurred, but with an enrichment of Lophodermella pathogens. Many pathogenicity-related genes were highly expressed by the mycobiota at the necrotrophic phase, showing an active pathogen response that are absent in asymptomatic needles. This study also revealed that Lophodermella spp. are members of a healthy needle mycobiota that have latent lifestyles suggesting that other pine needle pathogens may have similar biology. Interestingly, Pinus contorta upregulated defense genes in healthy needles, indicating response to fungal recognition, while a variety of biotic and abiotic stresses genes were activated in diseased needles. Further investigation to elucidate the possible antagonistic interplay of other biotic members leading to disease progression and/or suppression is warranted. This study provides insights into microbial interactions in non-model pathosystems and contributes to the development of new forest management strategies against emerging latent pathogens.

Although secondary metabolites are typically associated with competitive or pathogenic interactions, the high bioactivity of endophytic fungi in the Xylariales, coupled with their abundance and broad host ranges spanning all lineages of land plants and lichens, suggests that enhanced secondary metabolism might facilitate symbioses with phylogenetically diverse hosts. Here, we examined secondary metabolite gene clusters (SMGCs) across 96 Xylariales genomes in two clades (Xylariaceae s.l. and Hypoxylaceae), including 88 newly sequenced genomes of endophytes and closely related saprotrophs and pathogens. We paired genomic data with extensive metadata on endophyte hosts and substrates, enabling us to examine genomic factors related to the breadth of symbiotic interactions and ecological roles. All genomes contain hyperabundant SMGCs; however, Xylariaceae have increased numbers of gene duplications, horizontal gene transfers (HGTs) and SMGCs. Enhanced metabolic diversity of endophytes is associated with a greater diversity of hosts and increased capacity for lignocellulose decomposition. Our results suggest that, as host and substrate generalists, Xylariaceae endophytes experience greater selection to diversify SMGCs compared with more ecologically specialised Hypoxylaceae species. Overall, our results provide new evidence that SMGCs may facilitate symbiosis with phylogenetically diverse hosts, highlighting the importance of microbial symbioses to drive fungal metabolic diversity.

Pangenome analyses are increasingly being utilized to study the evolution of eukaryotic organisms. While pangenomes can provide insight into polymorphic gene content, inferences about the ecological and adaptive potential of such organisms also need to be accompanied by additional supportive genomic analyses. In this study we constructed a pangenome of Claviceps purpurea from 24 genomes and examined the positive selection and recombination landscape of an economically important fungal organism for pharmacology and agricultural research. Together, these analyses revealed that C. purpurea has a relatively large accessory genome (~ 38%), high recombination rates (ρ = 0.044), and transposon mediated gene duplication. However, due to observations of relatively low transposable element (TE) content (8.8%) and a lack of variability in genome sizes, prolific TE expansion may be controlled by frequent recombination. We additionally identified that within the ergoline biosynthetic cluster the lpsA1 and lpsA2 were the result of a recombination event. However, the high recombination rates observed in C. purpurea may be influencing an overall trend of purifying selection across the genome. These results showcase the use of selection and recombination landscapes to identify mechanisms contributing to pangenome structure and primary factors influencing the evolution of an organism.

The ascomycetous yeast Candida membranifaciens has been isolated from diverse habitats, including humans, insects, and environmental sources, exhibiting a remarkable ability to use different carbon sources that include pentoses, melibiose, and inulin. In this study, we isolated four C. membranifaciens strains from soil and investigated their potential to overproduce riboflavin. C. membranifaciens IST 626 was found to produce the highest concentrations of riboflavin. The volumetric production of this vitamin was higher when C. membranifaciens IST 626 cells were cultured in a commercial medium without iron and when xylose was the available carbon source compared to the same basal medium with glucose. Supplementation of the growth medium with 2 g/L glycine favored the metabolization of xylose, leading to biomass increase and consequent enhancement of riboflavin volumetric production that reached 120 mg/L after 216 h of cultivation. To gain new insights into the molecular basis of riboflavin production and carbon source utilization in this species, the first annotated genome sequence of C. membranifaciens is reported in this article, as well as the result of a comparative genomic analysis with other relevant yeast species. A total of 5619 genes were predicted to be present in C. membranifaciens IST 626 genome sequence (11.5 Mbp). Among them are genes involved in riboflavin biosynthesis, iron homeostasis, and sugar uptake and metabolism. This work put forward C. membranifaciens IST 626 as a riboflavin overproducer and provides valuable molecular data for future development of superior producing strains capable of using the wide range of carbon sources, which is a characteristic trait of the species.

Aspergillus flavus is an opportunistic pathogen of crops, including peanuts and maize, and is the second leading cause of aspergillosis in immunocompromised patients. A. flavus is also a major producer of the mycotoxin, aflatoxin, a potent carcinogen, which results in significant crop losses annually. The A. flavus isolate NRRL 3357 was originally isolated from peanut and has been used as a model organism for understanding the regulation and production of secondary metabolites, such as aflatoxin. A draft genome of NRRL 3357 was previously constructed, enabling the development of molecular tools and for understanding population biology of this particular species. Here, we describe an updated, near complete, telomere-to-telomere assembly and re-annotation of the eight chromosomes of A. flavus NRRL 3357 genome, accomplished via long-read PacBio and Oxford Nanopore technologies combined with Illumina short-read sequencing. A total of 13,715 protein-coding genes were predicted. Using RNA-seq data, a significant improvement was achieved in predicted 5' and 3' untranslated regions, which were incorporated into the new gene models.

Anaerobic fungi (class Neocallimastigomycetes) thrive as low-abundance members of the herbivore digestive tract. The genomes of anaerobic gut fungi are poorly characterized and have not been extensively mined for the biosynthetic enzymes of natural products such as antibiotics. Here, we investigate the potential of anaerobic gut fungi to synthesize natural products that could regulate membership within the gut microbiome. Complementary 'omics' approaches were combined to catalog the natural products of anaerobic gut fungi from four different representative species: Anaeromyces robustus (Arobustus), Caecomyces churrovis (Cchurrovis), Neocallimastix californiae (Ncaliforniae), and Piromyces finnis (Pfinnis). In total, 146 genes were identified that encode biosynthetic enzymes for diverse types of natural products, including nonribosomal peptide synthetases and polyketide synthases. In addition, N. californiae and C. churrovis genomes encoded seven putative bacteriocins, a class of antimicrobial peptides typically produced by bacteria. During standard laboratory growth on plant biomass or soluble substrates, 26% of total core biosynthetic genes in all four strains were transcribed. Across all four fungal strains, 30% of total biosynthetic gene products were detected via proteomics when grown on cellobiose. Liquid chromatography-tandem mass spectrometry (LC-MS/MS) characterization of fungal supernatants detected 72 likely natural products from A. robustus alone. A compound produced by all four strains of anaerobic fungi was putatively identified as the polyketide-related styrylpyrone baumin. Molecular networking quantified similarities between tandem mass spectrometry (MS/MS) spectra among these fungi, enabling three groups of natural products to be identified that are unique to anaerobic fungi. Overall, these results support the finding that anaerobic gut fungi synthesize natural products, which could be harnessed as a source of antimicrobials, therapeutics, and other bioactive compounds.

Anaerobic gut fungi in the phylum Neocallimastigomycota typically inhabit the digestive tracts of large mammalian herbivores, where they play an integral role in the decomposition of raw lignocellulose into its constitutive sugar monomers. However, quantitative tools to study their physiology are lacking, partially due to their complex and unresolved metabolism that includes the largely uncharacterized fungal hydrogenosome. Modern omics approaches combined with metabolic modeling can be used to establish an understanding of gut fungal metabolism and develop targeted engineering strategies to harness their degradation capabilities for lignocellulosic bioprocessing. Here, we introduce a high-quality genome of the anaerobic fungus Neocallimastix lanati from which we constructed the first genome-scale metabolic model of an anaerobic fungus. Relative to its size (200 Mbp, sequenced at 62× depth), it is the least fragmented publicly available gut fungal genome to date. Of the 1,788 lignocellulolytic enzymes annotated in the genome, 585 are associated with the fungal cellulosome, underscoring the powerful lignocellulolytic potential of N. lanati The genome-scale metabolic model captures the primary metabolism of N. lanati and accurately predicts experimentally validated substrate utilization requirements. Additionally, metabolic flux predictions are verified by 13C metabolic flux analysis, demonstrating that the model faithfully describes the underlying fungal metabolism. Furthermore, the model clarifies key aspects of the hydrogenosomal metabolism and can be used as a platform to quantitatively study these biotechnologically important yet poorly understood early-branching fungi.IMPORTANCE Recent genomic analyses have revealed that anaerobic gut fungi possess both the largest number and highest diversity of lignocellulolytic enzymes of all sequenced fungi, explaining their ability to decompose lignocellulosic substrates, e.g., agricultural waste, into fermentable sugars. Despite their potential, the development of engineering methods for these organisms has been slow due to their complex life cycle, understudied metabolism, and challenging anaerobic culture requirements. Currently, there is no framework that can be used to combine multi-omic data sets to understand their physiology. Here, we introduce a high-quality PacBio-sequenced genome of the anaerobic gut fungus Neocallimastix lanati Beyond identifying a trove of lignocellulolytic enzymes, we use this genome to construct the first genome-scale metabolic model of an anaerobic gut fungus. The model is experimentally validated and sheds light on unresolved metabolic features common to gut fungi. Model-guided analysis will pave the way for deepening our understanding of anaerobic gut fungi and provides a systematic framework to guide strain engineering efforts of these organisms for biotechnological use.

The genus Claviceps has been known for centuries as an economically important fungal genus for pharmacology and agricultural research. Only recently have researchers begun to unravel the evolutionary history of the genus, with origins in South America and classification of four distinct sections through ecological, morphological, and metabolic features (Claviceps sects. Citrinae, Paspalorum, Pusillae, and Claviceps). The first three sections are additionally characterized by narrow host range, whereas section Claviceps is considered evolutionarily more successful and adaptable as it has the largest host range and biogeographical distribution. However, the reasons for this success and adaptability remain unclear. Our study elucidates factors influencing adaptability by sequencing and annotating 50 Claviceps genomes, representing 21 species, for a comprehensive comparison of genome architecture and plasticity in relation to host range potential. Our results show the trajectory from specialized genomes (sects. Citrinae and Paspalorum) toward adaptive genomes (sects. Pusillae and Claviceps) through colocalization of transposable elements around predicted effectors and a putative loss of repeat-induced point mutation resulting in unconstrained tandem gene duplication coinciding with increased host range potential and speciation. Alterations of genomic architecture and plasticity can substantially influence and shape the evolutionary trajectory of fungal pathogens and their adaptability. Furthermore, our study provides a large increase in available genomic resources to propel future studies of Claviceps in pharmacology and agricultural research, as well as, research into deeper understanding of the evolution of adaptable plant pathogens.

Algae are a diverse, polyphyletic group of photosynthetic eukaryotes spanning nearly all eukaryotic lineages of life and collectively responsible for ∼50% of photosynthesis on Earth. Sequenced algal genomes, critical to understanding their complex biology, are growing in number and require efficient tools for analysis. PhycoCosm (https://phycocosm.jgi.doe.gov) is an algal multi-omics portal, developed by the US Department of Energy Joint Genome Institute to support analysis and distribution of algal genome sequences and other 'omics' data. PhycoCosm provides integration of genome sequence and annotation for >100 algal genomes with available multi-omics data and interactive web-based tools to enable algal research in bioenergy and the environment, encouraging community engagement and data exchange, and fostering new sequencing projects that will further these research goals.

Anaerobic fungi are a potential biotechnology platform to produce biomass-degrading enzymes. Unlike model fungi such as yeasts, stress responses that are relevant during bioprocessing have not yet been established for anaerobic fungi. In this work, we characterize both the heat shock and unfolded protein responses of four strains of anaerobic fungi (Anaeromyces robustus, Caecomyces churrovis, Neocallimastix californiae, and Piromyces finnis). The inositol-requiring 1 (Ire1) stress sensor, which typically initiates the fungal UPR, was conserved in all four genomes. However, these genomes also encode putative transmembrane kinases with catalytic domains that are similar to the metazoan stress-sensing enzyme PKR-like endoplasmic reticulum kinase (PERK), although whether they function in the UPR of anaerobic fungi remains unclear. Furthermore, we characterized the global transcriptional responses of Anaeromyces robustus and Neocallimastix californiae to a transient heat shock. Both fungi exhibited the hallmarks of ER stress, including upregulation of genes with functions in protein folding, ER-associated degradation, and intracellular protein trafficking. Relative to other fungi, the genomes of Neocallimastigomycetes contained the greatest gene percentage of HSP20 and HSP70 chaperones, which may serve to stabilize their asparagine-rich genomes. Taken together, these results delineate the unique stress response of anaerobic fungi, which is an important step toward their development as a biotechnology platform to produce enzymes and valuable biomolecules.

Anaerobic fungi and methanogenic archaea are two classes of microorganisms found in the rumen microbiome that metabolically interact during lignocellulose breakdown. Here, stable synthetic co-cultures of the anaerobic fungus Caecomyces churrovis and the methanogen Methanobacterium bryantii (not native to the rumen) were formed, demonstrating that microbes from different environments can be paired based on metabolic ties. Transcriptional and metabolic changes induced by methanogen co-culture were evaluated in C. churrovis across a variety of substrates to identify mechanisms that impact biomass breakdown and sugar uptake. A high-quality genome of C. churrovis was obtained and annotated, which is the first sequenced genome of a non-rhizoid-forming anaerobic fungus. C. churrovis possess an abundance of CAZymes and carbohydrate binding modules and, in agreement with previous studies of early-diverging fungal lineages, N6-methyldeoxyadenine (6mA) was associated with transcriptionally active genes. Co-culture with the methanogen increased overall transcription of CAZymes, carbohydrate binding modules, and dockerin domains in co-cultures grown on both lignocellulose and cellulose and caused upregulation of genes coding associated enzymatic machinery including carbohydrate binding modules in family 18 and dockerin domains across multiple growth substrates relative to C. churrovis monoculture. Two other fungal strains grown on a reed canary grass substrate in co-culture with the same methanogen also exhibited high log2-fold change values for upregulation of genes encoding carbohydrate binding modules in families 1 and 18. Transcriptional upregulation indicated that co-culture of the C. churrovis strain with a methanogen may enhance pyruvate formate lyase (PFL) function for growth on xylan and fructose and production of bottleneck enzymes in sugar utilization pathways, further supporting the hypothesis that co-culture with a methanogen may enhance certain fungal metabolic functions. Upregulation of CBM18 may play a role in fungal-methanogen physical associations and fungal cell wall development and remodeling.

Research into secondary metabolism (SM) production by fungi has resulted in the discovery of diverse, biologically active compounds with significant medicinal applications. The fungi rich in SM production are taxonomically concentrated in the subkingdom Dikarya, which comprises the phyla Ascomycota and Basidiomycota. Here, we explore the potential for SM production in Mucoromycota and Zoopagomycota, two phyla of nonflagellated fungi that are not members of Dikarya, by predicting and identifying core genes and gene clusters involved in SM. The majority of non-Dikarya have few genes and gene clusters involved in SM production except for the amphibian gut symbionts in the genus Basidiobolus Basidiobolus genomes exhibit an enrichment of SM genes involved in siderophore, surfactin-like, and terpene cyclase production, all these with evidence of constitutive gene expression. Gene expression and chemical assays also confirm that Basidiobolus has significant siderophore activity. The expansion of SMs in Basidiobolus are partially due to horizontal gene transfer from bacteria, likely as a consequence of its ecology as an amphibian gut endosymbiont.

Symbioses of bacteria with fungi have only recently been described and are poorly understood. In the symbiosis of Mycetohabitans (formerly Burkholderia) rhizoxinica with the fungus Rhizopus microsporus, bacterial type III (T3) secretion is known to be essential. Proteins resembling T3-secreted transcription activator-like (TAL) effectors of plant pathogenic bacteria are encoded in the three sequenced Mycetohabitans spp. genomes. TAL effectors nuclear-localize in plants, where they bind and activate genes important in disease. The Burkholderia TAL-like (Btl) proteins bind DNA but lack the N- and C-terminal regions, in which TAL effectors harbor their T3 and nuclear localization signals, and activation domain. We characterized a Btl protein, Btl19-13, and found that, despite the structural differences, it can be T3-secreted and can nuclear-localize. A btl19 -13 gene knockout did not prevent the bacterium from infecting the fungus, but the fungus became less tolerant to cell membrane stress. Btl19-13 did not alter transcription in a plant-based reporter assay, but 15 R. microsporus genes were differentially expressed in comparisons both of the fungus infected with the wild-type bacterium vs. the mutant and with the mutant vs. a complemented strain. Southern blotting revealed btl genes in 14 diverse Mycetohabitans isolates. However, banding patterns and available sequences suggest variation, and the btl19-13 phenotype could not be rescued by a btl gene from a different strain. Our findings support the conclusion that Btl proteins are effectors that act on host DNA and play important but varied or possibly host genotype-specific roles in the M. rhizoxinica-R. microsporus symbiosis.

Rhizopus microsporus is an early-diverging fungal species with importance in ecology, agriculture, food production, and public health. Pathogenic strains of R. microsporus harbor an intracellular bacterial symbiont, Mycetohabitans (formerly named Burkholderia). This vertically transmitted bacterial symbiont is responsible for the production of toxins crucial to the pathogenicity of Rhizopus and remarkably also for fungal reproduction. Here we show that R. microsporus can live not only in symbiosis with bacteria but also with two viral members of the genus Narnavirus. Our experiments revealed that both viruses replicated similarly in the growth conditions we tested. Viral copies were affected by the developmental stage of the fungus, the substrate, and the presence or absence of Mycetohabitans. Absolute quantification of narnaviruses in isolated asexual sporangiospores and sexual zygospores indicates their vertical transmission. By curing R. microsporus of its viral and bacterial symbionts and reinfecting bacteria to reestablish symbiosis, we demonstrate that these viruses affect fungal biology. Narnaviruses decrease asexual reproduction, but together with Mycetohabitans, are required for sexual reproductive success. This fungal-bacterial-viral system represents an outstanding model to investigate three-way microbial symbioses and their evolution.

Research into secondary metabolism (SM) production by fungi has resulted in the discovery of diverse, biologically active compounds with significant medicinal applications. However, the fungi rich in SM production are taxonomically restricted to Dikarya, two phyla of Kingdom Fungi, Ascomycota and Basidiomycota. Here, we explore the potential for SM production in Mucoromycota and Zoopagomycota, two phyla of nonflagellated fungi that are not members of Dikarya, by predicting and identifying core genes and gene clusters involved in SM. The majority of non-Dikarya have few genes and gene clusters involved in SM production except for the amphibian gut symbionts in the genus Basidiobolus . Basidiobolus genomes exhibit an enrichment of SM genes involved in siderophore, surfactin-like, and terpene cyclase production, all these with evidence of constitutive gene expression. Gene expression and chemical assays confirm that Basidiobolus has significant siderophore activity. The expansion of SMs in Basidiobolus are partially due to horizontal gene transfer from bacteria, likely as a consequence of its ecology as an amphibian gut endosymbiont.

Section Flavi encompasses both harmful and beneficial Aspergillus species, such as Aspergillus oryzae, used in food fermentation and enzyme production, and Aspergillus flavus, food spoiler and mycotoxin producer. Here, we sequence 19 genomes spanning section Flavi and compare 31 fungal genomes including 23 Flavi species. We reassess their phylogenetic relationships and show that the closest relative of A. oryzae is not A. flavus, but A. minisclerotigenes or A. aflatoxiformans and identify high genome diversity, especially in sub-telomeric regions. We predict abundant CAZymes (598 per species) and prolific secondary metabolite gene clusters (73 per species) in section Flavi. However, the observed phenotypes (growth characteristics, polysaccharide degradation) do not necessarily correlate with inferences made from the predicted CAZyme content. Our work, including genomic analyses, phenotypic assays, and identification of secondary metabolites, highlights the genetic and metabolic diversity within section Flavi.

Coniochaeta species are versatile ascomycetes that have great capacity to deconstruct lignocellulose. Here, we explore the transcriptome of Coniochaeta sp. strain 2T2.1 from wheat straw-driven cultures with the fungus growing alone or as a member of a synthetic microbial consortium with Sphingobacterium multivorum w15 and Citrobacter freundii so4. The differential expression profiles of carbohydrate-active enzymes indicated an onset of (hemi)cellulose degradation by 2T2.1 during the initial 24 hours of incubation. Within the tripartite consortium, 63 transcripts of strain 2T2.1 were differentially expressed at this time point. The presence of the two bacteria significantly upregulated the expression of one galactose oxidase, one GH79-like enzyme, one multidrug transporter, one laccase-like protein (AA1 family) and two bilirubin oxidases, suggesting that inter-kingdom interactions (e.g. amensalism) take place within this microbial consortium. Overexpression of multicopper oxidases indicated that strain 2T2.1 may be involved in lignin depolymerization (a trait of enzymatic synergism), while S. multivorum and C. freundii have the metabolic potential to deconstruct arabinoxylan. Under the conditions applied, 2T2.1 appears to be a better degrader of wheat straw when the two bacteria are absent. This conclusion is supported by the observed suppression of its (hemi)cellulolytic arsenal and lower degradation percentages within the microbial consortium.

BackgroundParticular species of the genus Coniochaeta (Sordariomycetes) exhibit great potential for bioabatement of furanic compounds and have been identified as an underexplored source of novel lignocellulolytic enzymes, especially Coniochaeta ligniaria. However, there is a lack of information about their genomic features and metabolic capabilities. Here, we report the first in-depth genome/transcriptome survey of a Coniochaeta species (strain 2T2.1).ResultsThe genome of Coniochaeta sp. strain 2T2.1 has a size of 74.53 Mbp and contains 24,735 protein-encoding genes. Interestingly, we detected a genome expansion event, resulting ~ 98% of the assembly being duplicated with 91.9% average nucleotide identity between the duplicated regions. The lack of gene loss, as well as the high divergence and strong genome-wide signatures of purifying selection between copies indicates that this is likely a recent duplication, which arose through hybridization between two related Coniochaeta-like species (allopolyploidization). Phylogenomic analysis revealed that 2T2.1 is related Coniochaeta sp. PMI546 and Lecythophora sp. AK0013, which both occur endophytically. Based on carbohydrate-active enzyme (CAZy) annotation, we observed that even after in silico removal of its duplicated content, the 2T2.1 genome contains exceptional lignocellulolytic machinery. Moreover, transcriptomic data reveal the overexpression of proteins affiliated to CAZy families GH11, GH10 (endoxylanases), CE5, CE1 (xylan esterases), GH62, GH51 (α-l-arabinofuranosidases), GH12, GH7 (cellulases), and AA9 (lytic polysaccharide monoxygenases) when the fungus was grown on wheat straw compared with glucose as the sole carbon source.ConclusionsWe provide data that suggest that a recent hybridization between the genomes of related species may have given rise to Coniochaeta sp. 2T2.1. Moreover, our results reveal that the degradation of arabinoxylan, xyloglucan and cellulose are key metabolic processes in strain 2T2.1 growing on wheat straw. Different genes for key lignocellulolytic enzymes were identified, which can be starting points for production, characterization and/or supplementation of enzyme cocktails used in saccharification of agricultural residues. Our findings represent first steps that enable a better understanding of the reticulate evolution and "eco-enzymology" of lignocellulolytic Coniochaeta species.

Centromeres are rapidly evolving across eukaryotes, despite performing a conserved function to ensure high-fidelity chromosome segregation. CENP-A chromatin is a hallmark of a functional centromere in most organisms. Due to its critical role in kinetochore architecture, the loss of CENP-A is tolerated in only a few organisms, many of which possess holocentric chromosomes. Here, we characterize the consequence of the loss of CENP-A in the fungal kingdom. Mucor circinelloides, an opportunistic human pathogen, lacks CENP-A along with the evolutionarily conserved CENP-C but assembles a monocentric chromosome with a localized kinetochore complex throughout the cell cycle. Mis12 and Dsn1, two conserved kinetochore proteins, were found to co-localize to a short region, one in each of nine large scaffolds, composed of an ∼200-bp AT-rich sequence followed by a centromere-specific conserved motif that echoes the structure of budding yeast point centromeres. Resembling fungal regional centromeres, these core centromere regions are embedded in large genomic expanses devoid of genes yet marked by Grem-LINE1s, a novel retrotransposable element silenced by the Dicer-dependent RNAi pathway. Our results suggest that these hybrid features of point and regional centromeres arose from the absence of CENP-A, thus defining novel mosaic centromeres in this early-diverging fungus.

The evolution of complex multicellularity has been one of the major transitions in the history of life. In contrast to simple multicellular aggregates of cells, it has evolved only in a handful of lineages, including animals, embryophytes, red and brown algae, and fungi. Despite being a key step toward the evolution of complex organisms, the evolutionary origins and the genetic underpinnings of complex multicellularity are incompletely known. The development of fungal fruiting bodies from a hyphal thallus represents a transition from simple to complex multicellularity that is inducible under laboratory conditions. We constructed a reference atlas of mushroom formation based on developmental transcriptome data of six species and comparisons of >200 whole genomes, to elucidate the core genetic program of complex multicellularity and fruiting body development in mushroom-forming fungi (Agaricomycetes). Nearly 300 conserved gene families and >70 functional groups contained developmentally regulated genes from five to six species, covering functions related to fungal cell wall remodeling, targeted protein degradation, signal transduction, adhesion, and small secreted proteins (including effector-like orphan genes). Several of these families, including F-box proteins, expansin-like proteins, protein kinases, and transcription factors, showed expansions in Agaricomycetes, many of which convergently expanded in multicellular plants and/or animals too, reflecting convergent solutions to genetic hurdles imposed by complex multicellularity among independently evolved lineages. This study provides an entry point to studying mushroom development and complex multicellularity in one of the largest clades of complex eukaryotic organisms.

The generation of thousands of fungal genomes is leading to a better understanding of genes and genomic organization within the kingdom. However, the epigenome, which includes DNA and chromatin modifications, remains poorly investigated in fungi. Large comparative studies in animals and plants have deepened our understanding of epigenomic variation, particularly of the modified base 5-methylcytosine (5mC), but taxonomic sampling of disparate groups is needed to develop unifying explanations for 5mC variation. Here, we utilize the largest phylogenetic resolution of 5mC methyltransferases (5mC MTases) and genome evolution to better understand levels and patterns of 5mC across fungi. We show that extant 5mC MTase genotypes are descendent from ancestral maintenance and de novo genotypes, whereas the 5mC MTases DIM-2 and RID are more recently derived, and that 5mC levels are correlated with 5mC MTase genotype and transposon content. Our survey also revealed that fungi lack canonical gene-body methylation, which distinguishes fungal epigenomes from certain insect and plant species. However, some fungal species possess independently derived clusters of contiguous 5mC encompassing many genes. In some cases, DNA repair pathways and the N6-methyladenine DNA modification negatively coevolved with 5mC pathways, which additionally contributed to interspecific epigenomic variation across fungi.

Arbuscular mycorrhizal fungi (AMF, Glomeromycotina), in addition to forming symbioses with the majority of land plants, harbor vertically transmitted endosymbiotic bacteria 'Candidatus Glomeribacter gigasporarum' (CaGg) and 'Candidatus Moeniiplasma glomeromycotorum' (CaMg). CaGg is a nonessential mutualist of AMF, whereas the lifestyle of CaMg is unknown. To start unraveling the interactions between AMF and their endosymbionts in nature, we examined diversity and distribution of AMF-associated endobacteria in North Atlantic dunes at Cape Cod. Of nearly 500 foredune AMF isolates successfully genotyped during a systematic study, 94% were classified as Gigasporaceae. Two percent of all AMF spores harbored CaGg, and 88% contained CaMg. CaGg was found only in the Gigasporaceae, whereas CaMg was present in Gigasporaceae, Acaulosporaceae, and Diversisporaceae. Incidence of CaGg across AMF was not affected by any of the environmental parameters measured, whereas distribution of CaMg in one of the fungal hosts was impacted by plant density. CaMg populations associated with AMF individuals displayed high levels of genetic diversity but no evidence of gene flow, suggesting that host physical proximity is not sufficient to facilitate horizontal transmission of CaMg. Finally, in addition to a novel lineage of CaGg, we discovered that AMF likely harbor Burkholderia-related bacteria with close phylogenetic affinity to free-living Burkholderia and endobacteria of other Mucoromycota fungi.

Ustilaginomycotina is home to a broad array of fungi including important plant pathogens collectively called smut fungi. Smuts are biotrophs that produce characteristic perennating propagules called teliospores, one of which, Ustilago maydis, is a model genetic organism. Broad exploration of smut biology has been hampered by limited phylogenetic resolution of Ustilaginiomycotina as well as an overall lack of genomic data for members of this subphylum. In this study, we sequenced eight Ustilaginomycotina genomes from previously unrepresented lineages, deciphered ordinal-level phylogenetic relationships for the subphylum, and performed comparative analyses. Unlike other Basidiomycota subphyla, all sampled Ustilaginomycotina genomes are relatively small and compact. Ancestral state reconstruction analyses indicate that teliospore formation was present at the origin of the subphylum. Divergence time estimation dates the divergence of most extant smut fungi after that of grasses (Poaceae). However, we found limited conservation of well-characterized genes related to smut pathogenesis from U. maydis, indicating dissimilar pathogenic mechanisms exist across other smut lineages. The genomes of Malasseziomycetes are highly diverged from the other sampled Ustilaginomycotina, likely due to their unique history as mammal-associated lipophilic yeasts. Despite extensive genomic data, the phylogenetic placement of this class remains ambiguous. Although the sampled Ustilaginomycotina members lack many core enzymes for plant cell wall decomposition and starch catabolism, we identified several novel carbohydrate active enzymes potentially related to pectin breakdown. Finally, ∼50% of Ustilaginomycotina species-specific genes are present in previously undersampled and rare lineages, highlighting the importance of exploring fungal diversity as a resource for novel gene discovery.

Phylogenomic approaches have the potential to improve confidence about the inter-relationships of species in the order Mucorales within the fungal tree of life. Rhizopus species are especially important as plant and animal pathogens and bioindustrial fermenters for food and metabolite production. A dataset of 192 orthologous genes was used to construct a phylogenetic tree of 21 Rhizopus strains, classified into four species isolated from habitats of industrial, medical and environmental importance. The phylogeny indicates that the genus Rhizopus consists of three major clades, with R. microsporus as the basal species and the sister lineage to R. stolonifer and two closely related species R. arrhizus and R. delemar A comparative analysis of the mating type locus across Rhizopus reveals that its structure is flexible even between different species in the same genus, but shows similarities between Rhizopus and other mucoralean fungi. The topology of single-gene phylogenies built for two genes involved in mating is similar to the phylogenomic tree. Comparison of the total length of the genome assemblies showed that genome size varies by as much as threefold within a species and is driven by changes in transposable element copy numbers and genome duplications.

The basidiomycete yeast Rhodosporidium toruloides (also known as Rhodotorula toruloides) accumulates high concentrations of lipids and carotenoids from diverse carbon sources. It has great potential as a model for the cellular biology of lipid droplets and for sustainable chemical production. We developed a method for high-throughput genetics (RB-TDNAseq), using sequence-barcoded Agrobacterium tumefaciens T-DNA insertions. We identified 1,337 putative essential genes with low T-DNA insertion rates. We functionally profiled genes required for fatty acid catabolism and lipid accumulation, validating results with 35 targeted deletion strains. We identified a high-confidence set of 150 genes affecting lipid accumulation, including genes with predicted function in signaling cascades, gene expression, protein modification and vesicular trafficking, autophagy, amino acid synthesis and tRNA modification, and genes of unknown function. These results greatly advance our understanding of lipid metabolism in this oleaginous species and demonstrate a general approach for barcoded mutagenesis that should enable functional genomics in diverse fungi.

The fungal genus of Aspergillus is highly interesting, containing everything from industrial cell factories, model organisms, and human pathogens. In particular, this group has a prolific production of bioactive secondary metabolites (SMs). In this work, four diverse Aspergillus species (A. campestris, A. novofumigatus, A. ochraceoroseus, and A. steynii) have been whole-genome PacBio sequenced to provide genetic references in three Aspergillus sections. A. taichungensis and A. candidus also were sequenced for SM elucidation. Thirteen Aspergillus genomes were analyzed with comparative genomics to determine phylogeny and genetic diversity, showing that each presented genome contains 15-27% genes not found in other sequenced Aspergilli. In particular, A. novofumigatus was compared with the pathogenic species A. fumigatus This suggests that A. novofumigatus can produce most of the same allergens, virulence, and pathogenicity factors as A. fumigatus, suggesting that A. novofumigatus could be as pathogenic as A. fumigatus Furthermore, SMs were linked to gene clusters based on biological and chemical knowledge and analysis, genome sequences, and predictive algorithms. We thus identify putative SM clusters for aflatoxin, chlorflavonin, and ochrindol in A. ochraceoroseus, A. campestris, and A. steynii, respectively, and novofumigatonin, ent-cycloechinulin, and epi-aszonalenins in A. novofumigatus Our study delivers six fungal genomes, showing the large diversity found in the Aspergillus genus; highlights the potential for discovery of beneficial or harmful SMs; and supports reports of A. novofumigatus pathogenicity. It also shows how biological, biochemical, and genomic information can be combined to identify genes involved in the biosynthesis of specific SMs.

The draft genome, morphological description, and phylogenetic placement of Coniella lustricola sp. nov. (Schizoparmeaceae) are provided. The species was isolated from submerged detritus in a fen at Black Moshannon State Park, Pennsylvania, USA and differs from all other Coniella species by having ellipsoid to fusoid, inequilateral conidia that are rounded on one end and truncate or obtuse on the other end, with a length to width ratio of 2.8. The draft genome is 36.56 Mbp and consists of 870 contigs on 634 scaffolds (L50 = 0.14 Mb, N50 = 76 scaffolds), with 0.5% of the scaffold length in gaps. It contains 11,317 predicted gene models, including predicted genes for cellulose, hemicellulose, and xylan degradation, as well as predicted regions encoding for amylase, laccase, and tannase enzymes. Many members of the Schizoparmeaceae are plant pathogens of agricultural crops. This draft genome represents the first sequenced Coniella genome and will be a valuable tool for comparisons among pathogenic Coniella species.

Species in the genus Paecilomyces, a member of the fungal order Eurotiales, are ubiquitous in nature and impact a variety of human endeavors. Here, the biology of one common species, Paecilomyces variotii, was explored using genomics and functional genetics. Sequencing the genome of two isolates revealed key genome and gene features in this species. A striking feature of the genome was the two-part nature, featuring large stretches of DNA with normal GC content separated by AT-rich regions, a hallmark of many plant-pathogenic fungal genomes. These AT-rich regions appeared to have been mutated by repeat-induced point (RIP) mutations. We developed methods for genetic transformation of P. variotii, including forward and reverse genetics as well as crossing techniques. Using transformation and crossing, RIP activity was identified, demonstrating for the first time that RIP is an active process within the order Eurotiales. A consequence of RIP is likely reflected by a reduction in numbers of genes within gene families, such as in cell wall degradation, and reflected by growth limitations on P. variotii on diverse carbon sources. Furthermore, using these transformation tools we characterized a conserved protein containing a domain of unknown function (DUF1212) and discovered it is involved in pigmentation.