Early-diverging Fungi Research Articles

Unveiling species diversity within early-diverging fungi from China I: three new species of Backusella (Backusellaceae, Mucoromycota).

The species diversity of early-diverging fungi has long lagged behind that of higher fungi, posing a significant obstacle to our comprehensive understanding of the fungal kingdom. Our ongoing research endeavors aim to address this gap by exploring the species diversity of early-diverging fungi in China. In this study, we describe three novel species within the Backusella, namely B.elliptica sp. nov., B.fujianensis sp. nov., and B.variispora sp. nov., based on phylogenetic and morphological analyses. In the phylogenetic analysis of the ITS (internal transcribed spacer), LSU (large subunit of ribosomal RNA gene), and RPB1 (RNA polymerase II largest subunit gene) regions, the B.elliptica and B.fujianensis cluster closely with B.gigacellularis, B.ovalispora, and B.solicola, and the B.variispora is closely related to B.locustae and B.pernambucensis. Morphologically, B.elliptica is distinguished by elliptical sporangiospores, as well as cylindrical and hemispherical columellae. The B.fujianensis is characterized by elliptical sporangiospores, and various types of columellae such as hemispherical, subglobose, depressed globose and conical. The B.variispora is characterized by subglobose to globose sporangiospores, as well as hemispherical, subglobose to globose columellae. Additionally, the sporangiophores are long and monopodially branched in B.elliptica and B.fujianensis, while short and simple or sympodially branched in B.variispora. Physiologically, the maximum growth temperatures of B.elliptica (32 °C), B.fujianensis (35 °C), and B.variispora were (35 °C) were determined. With the inclusion of these newly described taxa, the total number of Backusella species known from China now stands at 12. Finally, we provide a key to facilitate the morphological identification of Backusella species from Asia.

Symmetric and asymmetric DNA N6-adenine methylation regulates different biological responses in Mucorales

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.

A minimal Fanconi Anemia complex in early diverging fungi

Fanconi Anemia (FA) pathway resolves DNA interstrand cross links (ICL). The FA pathway was initially recognized in vertebrates, but was later confirmed in other animals and speculated in fungi. FA proteins FANCM, FANCL and FANCJ are present in Saccharomyces cerevisiae but, their mechanism of interaction to resolve ICL is still unclear. Unlike Dikarya, early diverging fungi (EDF) possess more traits shared with animals. We traced the evolutionary history of the FA pathway across Opisthokonta. We scanned complete proteomes for FA-related homologs to establish their taxonomic distribution and analyzed their phylogenetic trees. We checked transcription profiles of FA genes to test if they respond to environmental conditions and their genomic localizations for potential co-localization. We identified fungal homologs of the activation and ID complexes, 5 out of 8 core proteins, all of the endonucleases, and deubiquitination proteins. All fungi lack FANCC, FANCF and FANCG proteins responsible for post-replication repair and chromosome stability in animals. The observed taxonomic distribution can be attributed to a gradual degradation of the FA pathway from EDF to Dikarya. One of the key differences is that EDF have the ID complex recruiting endonucleases to the site of ICL. Moreover, 21 out of 32 identified FA genes are upregulated in response to different growth conditions. Several FA genes are co-localized in fungal genomes which also could facilitate co-expression. Our results indicate that a minimal FA pathway might still be functional in Mucoromycota with a gradual loss of components in Dikarya ancestors.

Microbial hitchhikers on microplastics: The exchange of aquatic microbes across distinct aquatic habitats.

Microplastics (MPs) have the potential to modify aquatic microbial communities and distribute microorganisms, including pathogens. This poses a potential risk to aquatic life and human health. Despite this, the fate of 'hitchhiking' microbes on MPs that traverse different aquatic habitats remains largely unknown. To address this, we conducted a 50-day microcosm experiment, manipulating estuarine conditions to study the exchange of bacteria and microeukaryotes between river, sea and plastisphere using a long-read metabarcoding approach. Our findings revealed a significant increase in bacteria on the plastisphere, including Pseudomonas, Sphingomonas, Hyphomonas, Brevundimonas, Aquabacterium and Thalassolituus, all of which are known for their pollutant degradation capabilities, specifically polycyclic aromatic hydrocarbons. We also observed a strong association of plastic-degrading fungi (i.e., Cladosporium and Plectosphaerella) and early-diverging fungi (Cryptomycota, also known as Rozellomycota) with the plastisphere. Sea MPs were primarily colonised by fungi (70%), with a small proportion of river-transported microbes (1%-4%). The mere presence of MPs in seawater increased the relative abundance of planktonic fungi from 2% to 25%, suggesting significant exchanges between planktonic and plastisphere communities. Using microbial source tracking, we discovered that MPs only dispersed 3.5% and 5.5% of river bacterial and microeukaryotic communities into the sea, respectively. Hence, although MPs select and facilitate the dispersal of ecologically significant microorganisms, drastic compositional changes across distinct aquatic habitats are unlikely.

The interkingdom horizontal gene transfer in 44 early diverging fungi boosted their metabolic, adaptive, and immune capabilities.

Numerous studies have been devoted to individual cases of horizontally acquired genes in fungi. It has been shown that such genes expand the hosts' metabolic capabilities and contribute to their adaptations as parasites or symbionts. Some studies have provided an extensive characterization of the horizontal gene transfer (HGT) in Dikarya. However, in the early diverging fungi (EDF), a similar characterization is still missing. In order to fill this gap, we have designed a computational pipeline to obtain a statistical sample of reliable HGT events with a low false discovery rate. We have analyzed 44 EDF proteomes and identified 829 xenologs in fungi ranging from Chytridiomycota to Mucoromycota. We have identified several patterns and statistical properties of EDF HGT. We show that HGT is driven by bursts of gene exchange and duplication, resulting in highly divergent numbers and molecular properties of xenologs between fungal lineages. Ancestrally aquatic fungi are generally more likely to acquire foreign genetic material than terrestrial ones. Endosymbiotic bacteria can be a source of useful xenologs, as exemplified by NOD-like receptors transferred to Mortierellomycota. Closely related fungi have similar rates of intronization of xenologs. Posttransfer gene fusions and losses of protein domains are common and may influence the encoded proteins' functions. We argue that there is no universal approach for HGT identification and inter- and intra-kingdom transfers require tailored identification methods. Our results help to better understand how and to what extent HGT has shaped the metabolic, adaptive, and immune capabilities of fungi.

Genetic transformation of the frog-killing chytrid fungus Batrachochytrium dendrobatidis

Batrachochytrium dendrobatidis (Bd), a causative agent of chytridiomycosis, is decimating amphibian populations around the world. Bd belongs to the chytrid lineage, a group of early-diverging fungi that are widely used to study fungal evolution. Like all chytrids, Bd develops from a motile form into a sessile, growth form, a transition that involves drastic changes in its cytoskeletal architecture. Efforts to study Bd cell biology, development, and pathogenicity have been limited by the lack of genetic tools with which to test hypotheses about underlying molecular mechanisms. Here, we report the development of a transient genetic transformation system for Bd. We used electroporation to deliver exogenous DNA into Bd cells and detected transgene expression for up to three generations under both heterologous and native promoters. We also adapted the transformation protocol for selection using an antibiotic resistance marker. Finally, we used this system to express fluorescent protein fusions and, as a proof of concept, expressed a genetically encoded probe for the actin cytoskeleton. Using live-cell imaging, we visualized the distribution and dynamics of polymerized actin at each stage of the Bd life cycle, as well as during key developmental transitions. This transformation system enables direct testing of key hypotheses regarding mechanisms of Bd pathogenesis. This technology also paves the way for answering fundamental questions of chytrid cell, developmental, and evolutionary biology.

Phylogenetic and functional analyses of N6-methyladenosine RNA methylation factors in the wheat scab fungus Fusarium graminearum.

In eukaryotes, N6-methyladenosine (m6A) RNA modification plays a crucial role in governing the fate of RNA molecules and has been linked to various developmental processes. However, the phyletic distribution and functions of genetic factors responsible for m6A modification remain largely unexplored in fungi. To get insights into the evolution of m6A machineries, we reconstructed global phylogenies of potential m6A writers, readers, and erasers in fungi. Substantial copy number variations were observed, ranging from up to five m6A writers in early-diverging fungi to a single copy in the subphylum Pezizomycotina, which primarily comprises filamentous fungi. To characterize m6A factors in a phytopathogenic fungus Fusarium graminearum, we generated knockout mutants lacking potential m6A factors including the sole m6A writer MTA1. However, the resulting knockouts did not exhibit any noticeable phenotypic changes during vegetative and sexual growth stages. As obtaining a homozygous knockout lacking MTA1 was likely hindered by its essential role, we generated MTA1-overexpressing strains (MTA1-OE). The MTA1-OE5 strain showed delayed conidial germination and reduced hyphal branching, suggesting its involvement during vegetative growth. Consistent with these findings, the expression levels of MTA1 and a potential m6A reader YTH1 were dramatically induced in germinating conidia, followed by the expression of potential m6A erasers at later vegetative stages. Several genes including transcription factors, transporters, and various enzymes were found to be significantly upregulated and downregulated in the MTA1-OE5 strain. Overall, our study highlights the functional importance of the m6A methylation during conidial germination in F. graminearum and provides a foundation for future investigations into m6A modification sites in filamentous fungi.IMPORTANCEN6-methyladenosine (m6A) RNA methylation is a reversible posttranscriptional modification that regulates RNA function and plays a crucial role in diverse developmental processes. This study addresses the knowledge gap regarding phyletic distribution and functions of m6A factors in fungi. The identification of copy number variations among fungal groups enriches our knowledge regarding the evolution of m6A machinery in fungi. Functional characterization of m6A factors in a phytopathogenic filamentous fungus Fusarium graminearum provides insights into the essential role of the m6A writer MTA1 in conidial germination and hyphal branching. The observed effects of overexpressing MTA1 on fungal growth and gene expression patterns of m6A factors throughout the life cycle of F. graminearum further underscore the importance of m6A modification in conidial germination. Overall, this study significantly advances our understanding of m6A modification in fungi, paving the way for future research into its roles in filamentous growth and potential applications in disease control.

Non-canonical two-step biosynthesis of anti-oomycete indole alkaloids in Kickxellales

BackgroundFungi are prolific producers of bioactive small molecules of pharmaceutical or agricultural interest. The secondary metabolism of higher fungi (Dikarya) has been well-investigated which led to > 39,000 described compounds. However, natural product researchers scarcely drew attention to early-diverging fungi (Mucoro- and Zoopagomycota) as they are considered to rarely produce secondary metabolites. Indeed, only 15 compounds have as yet been isolated from the entire phylum of the Zoopagomycota.ResultsHere, we showcase eight species of the order Kickxellales (phylum Zoopagomycota) as potent producers of the indole-3-acetic acid (IAA)-derived compounds lindolins A and B. The compounds are produced both under laboratory conditions and in the natural soil habitat suggesting a specialized ecological function. Indeed, lindolin A is a selective agent against plant-pathogenic oomycetes such as Phytophthora sp. Lindolin biosynthesis was reconstituted in vitro and relies on the activity of two enzymes of dissimilar evolutionary origin: Whilst the IAA–CoA ligase LinA has evolved from fungal 4-coumaryl-CoA synthetases, the subsequently acting IAA-CoA:anthranilate N-indole-3-acetyltransferase LinB is a unique enzyme across all kingdoms of life.ConclusionsThis is the first report on bioactive secondary metabolites in the subphylum Kickxellomycotina and the first evidence for a non-clustered, two-step biosynthetic route of secondary metabolites in early-diverging fungi. Thus, the generally accepted “gene cluster hypothesis” for natural products needs to be reconsidered for early diverging fungi.

A combined microscopy and single-cell sequencing approach reveals the ecology, morphology, and phylogeny of uncultured lineages of zoosporic fungi.

Environmental DNA analyses of fungal communities typically reveal a much larger diversity than can be ascribed to known species. Much of this hidden diversity lies within undescribed fungal lineages, especially the early diverging fungi (EDF). Although these EDF often represent new lineages even at the phylum level, they have never been cultured, making their morphology and ecology uncertain. One of the methods to characterize these uncultured fungi is a single-cell DNA sequencing approach. In this study, we established a large data set of single-cell sequences of EDF by manually isolating and photographing parasitic fungi on various hosts such as algae, protists, and micro-invertebrates, combined with subsequent long-read sequencing of the ribosomal DNA locus (rDNA). We successfully obtained rDNA sequences of 127 parasitic fungal cells, which clustered into 71 phylogenetic lineages belonging to seven phylum-level clades of EDF: Blastocladiomycota, Chytridiomycota, Aphelidiomycota, Rozellomycota, and three unknown phylum-level clades. Most of our single cells yielded novel sequences distinguished from both described taxa and existing metabarcoding data, indicating an expansive and hidden diversity of parasitic taxa of EDF. We also revealed an unexpected diversity of endobiotic Olpidium-like chytrids and hyper-parasitic lineages. Overall, by combining photographs of parasitic fungi with phylogenetic analyses, we were able to better understand the ecological function and morphology of many of the branches on the fungal tree of life known only from DNA sequences. IMPORTANCE Much of the diversity of microbes from natural habitats, such as soil and freshwater, comprise species and lineages that have never been isolated into pure culture. In part, this stems from a bias of culturing in favor of saprotrophic microbes over the myriad symbiotic ones that include parasitic and mutualistic relationships with other taxa. In the present study, we aimed to shed light on the ecological function and morphology of the many undescribed lineages of aquatic fungi by individually isolating and sequencing molecular barcodes from 127 cells of host-associated fungi using single-cell sequencing. By adding these sequences and their photographs into the fungal tree, we were able to understand the morphology of reproductive and vegetative structures of these novel fungi and to provide a hypothesized ecological function for them. These individual host-fungal cells revealed themselves to be complex environments despite their small size; numerous samples were hyper-parasitized with other zoosporic fungal lineages such as Rozellomycota.

Sequencing the Genomes of the First Terrestrial Fungal Lineages: What Have We Learned?

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.

Chytrids in Soil Environments: Unique Adaptations and Distributions

Chytridiomycota (zoosporic true fungi) have a consistent presence in soils and have been frequently identified within many diverse terrestrial environments. However, Chytridiomycota and other early-diverging fungi have low representation in whole-genome sequencing databases compared to Dikarya. New molecular techniques have provided insights into the diversity and abundance of chytrids in soils and the changes in their populations both spatially and temporally. Chytrids complete their life cycle within rapidly changing soil environments where they may be more common within micropores due to protection from predation, desiccation, and extreme temperatures. Reproductive and morphological changes occur in response to environmental changes including pH, fluctuating nutrient concentrations, and metals at levels above toxic thresholds. Rhizoids share some features of hyphae, including the spatial regulation of branching and the ability to attach, adapt to, and proliferate in different substrates, albeit on a microscale. Soil chytrids provide a pool of novel enzymes and proteins which enable a range of lifestyles as saprotrophs or parasites, but also can be utilised as alternative tools with some biotechnological applications. Thus, 3D live-cell imaging and micromodels such as MicroCT may provide insight into zoospore functions and rhizoid plasticity, respectively, in response to various conditions. A combination of classical techniques of soil chytrid baiting with simultaneous molecular and ecological data will provide insights into temporal population changes in response to environmental change. The authors emphasise the need to review and improve DNA-based methodologies for identifying and quantifying chytrids within the soil microbiome to expand our knowledge of their taxonomy, abundance, diversity, and functionality within soil environments.

Heterochromatin and RNAi act independently to ensure genome stability in Mucorales human fungal pathogens

Chromatin modifications play a fundamental role in controlling transcription and genome stability and yet despite their importance, are poorly understood in early-diverging fungi. We present a comprehensive study of histone lysine and DNA methyltransferases across the Mucoromycota, emphasizing heterochromatin formation pathways that rely on the Clr4 complex involved in H3K9-methylation, the Polycomb-repressive complex 2 driving H3K27-methylation, or DNMT1-like methyltransferases that catalyze 5mC DNA methylation. Our analysis uncovered H3K9-methylated heterochromatin as the major chromatin modification repressing transcription in these fungi, which lack both Polycomb silencing and cytosine methylation. Although small RNAs generated by RNA interference (RNAi) pathways facilitate the formation of heterochromatin in many eukaryotic organisms, we show that RNAi is not required to maintain either genomic or centromeric heterochromatin in Mucor. H3K9-methylation and RNAi act independently to control centromeric regions, suggesting a functional subspecialization. Whereas the H3K9 methyltransferase Clr4 and heterochromatin formation are essential for cell viability, RNAi is dispensable for viability yet acts as the main epigenetic, regulatory force repressing transposition of centromeric GremLINE1 elements. Mutations inactivating canonical RNAi lead to rampant transposition and insertional inactivation of targets resulting in antimicrobial drug resistance. This fine-tuned, Rdrp2-dependent RNAi activity is critical for genome stability, restricting GremLINE1 retroelements to the centromeres where they occupy long heterochromatic islands. Taken together, our results suggest that RNAi and heterochromatin formation are independent genome defense and regulatory mechanisms in the Mucorales, contributing to a paradigm shift from the cotranscriptional gene silencing observed in fission yeasts to models in which heterochromatin and RNAi operate independently in early-diverging fungi.

A genetic tool to express long fungal biosynthetic genes

BackgroundSecondary metabolites (SMs) from mushroom-forming fungi (Basidiomycota) and early diverging fungi (EDF) such as Mucoromycota are scarcely investigated. In many cases, production of SMs is induced by unknown stress factors or is accompanied by seasonable developmental changes on fungal morphology. Moreover, many of these fungi are considered as non-culturable under laboratory conditions which impedes investigation into SM. In the post-genomic era, numerous novel SM genes have been identified especially from EDF. As most of them encode multi-module enzymes, these genes are usually long which limits cloning and heterologous expression in traditional hosts.ResultsAn expression system in Aspergillus niger is presented that is suitable for the production of SMs from both Basidiomycota and EDF. The akuB gene was deleted in the expression host A. niger ATNT∆pyrG, resulting in a deficient nonhomologous end-joining repair mechanism which in turn facilitates the targeted gene deletion via homologous recombination. The ∆akuB mutant tLK01 served as a platform to integrate overlapping DNA fragments of long SM genes into the fwnA locus required for the black pigmentation of conidia. This enables an easy discrimination of correct transformants by screening the transformation plates for fawn-colored colonies. Expression of the gene of interest (GOI) is induced dose-dependently by addition of doxycycline and is enhanced by the dual TetON/terrein synthase promoter system (ATNT) from Aspergillus terreus. We show that the 8 kb polyketide synthase gene lpaA from the basidiomycete Laetiporus sulphureus is correctly assembled from five overlapping DNA fragments and laetiporic acids are produced. In a second approach, we expressed the yet uncharacterized > 20 kb nonribosomal peptide synthetase gene calA from the EDF Mortierella alpina. Gene expression and subsequent LC–MS/MS analysis of mycelial extracts revealed the production of the antimycobacterial compound calpinactam. This is the first report on the heterologous production of a full-length SM multidomain enzyme from EDF.ConclusionsThe system allows the assembly, targeted integration and expression of genes of > 20 kb size in A. niger in one single step. The system is suitable for evolutionary distantly related SM genes from both Basidiomycota and EDF. This uncovers new SM resources including genetically intractable or non-culturable fungi.

Co-occurrences enhance our understanding of aquatic fungal metacommunity assembly and reveal potential host-parasite interactions.

Our knowledge of aquatic fungal communities, their assembly, distributions and ecological roles in marine ecosystems is scarce. Hence, we aimed to investigate fungal metacommunities of coastal habitats in a subarctic zone (northern Baltic Sea, Sweden). Using a novel joint species distribution model and network approach, we quantified the importance of biotic associations contributing to the assembly of mycoplankton, further, detected potential biotic interactions between fungi-algae pairs, respectively. Our long-read metabarcoding approach identified 493 fungal taxa, of which a dominant fraction (44.4%) was assigned as early-diverging fungi (i.e. Cryptomycota and Chytridiomycota). Alpha diversity of mycoplankton declined and community compositions changed along inlet-bay-offshore transects. The distributions of most fungi were rather influenced by environmental factors than by spatial drivers, and the influence of biotic associations was pronounced when environmental filtering was weak. We found great number of co-occurrences (120) among the dominant fungal groups, and the 25 associations between fungal and algal OTUs suggested potential host-parasite and/or saprotroph links, supporting a Cryptomycota-based mycoloop pathway. We emphasize that the contribution of biotic associations to mycoplankton assembly are important to consider in future studies as it helps to improve predictions of species distributions in aquatic ecosystems.

Analysis of Cyp51 protein sequences shows 4 major Cyp51 gene family groups across fungi.

Azole drugs target fungal sterol biosynthesis and are used to treat millions of human fungal infections each year. Resistance to azole drugs has emerged in multiple fungal pathogens including Candida albicans, Cryptococcus neoformans, Histoplasma capsulatum, and Aspergillus fumigatus. The most well-studied resistance mechanism in A. fumigatus arises from missense mutations in the coding sequence combined with a tandem repeat in the promoter of cyp51A, which encodes a cytochrome P450 enzyme in the fungal sterol biosynthesis pathway. Filamentous members of Ascomycota such as A. fumigatus have either 1 or 2 of 3 Cyp51 paralogs (Cyp51A, Cyp51B, and Cyp51C). Most previous research in A. fumigatus has focused on Cyp51A due to its role in azole resistance. We used the A. fumigatus Cyp51A protein sequence as the query in database searches to identify Cyp51 proteins across fungi. We found 435 Cyp51 proteins in 295 species spanning from early-diverging fungi (Blastocladiomycota, Chytridiomycota, Zoopagomycota, and Mucormycota) to late-diverging fungi (Ascomycota and Basidiomycota). We found these sequences formed 4 major Cyp51 groups: Cyp51, Cyp51A, Cyp51B, and Cyp51C. Surprisingly, we found all filamentous Ascomycota had a Cyp51B paralog, while only 50% had a Cyp51A paralog. We created maximum likelihood trees to investigate the evolution of Cyp51 in fungi. Our results suggest Cyp51 is present in all fungi with 3 paralogs emerging in Pezizomycotina, including Cyp51C which appears to have diverged from the progenitor of the Cyp51A and Cyp51B groups.

In Silico Predictions of Ecological Plasticity Mediated by Protein Family Expansions in Early-Diverging Fungi.

Early-diverging fungi (EDF) are ubiquitous and versatile. Their diversity is reflected in their genome sizes and complexity. For instance, multiple protein families have been reported to expand or disappear either in particular genomes or even whole lineages. The most commonly mentioned are CAZymes (carbohydrate-active enzymes), peptidases and transporters that serve multiple biological roles connected to, e.g., metabolism and nutrients intake. In order to study the link between ecology and its genomic underpinnings in a more comprehensive manner, we carried out a systematic in silico survey of protein family expansions and losses among EDF with diverse lifestyles. We found that 86 protein families are represented differently according to EDF ecological features (assessed by median count differences). Among these there are 19 families of proteases, 43 CAZymes and 24 transporters. Some of these protein families have been recognized before as serine and metallopeptidases, cellulases and other nutrition-related enzymes. Other clearly pronounced differences refer to cell wall remodelling and glycosylation. We hypothesize that these protein families altogether define the preliminary fungal adaptasome. However, our findings need experimental validation. Many of the protein families have never been characterized in fungi and are discussed in the light of fungal ecology for the first time.

Understudied, underrepresented, and unknown: Methodological biases that limit detection of early diverging fungi from environmental samples.

Metabarcoding is an important tool for understanding fungal communities. The internal transcribed spacer (ITS) rDNA is the accepted fungal barcode but has known problems. The large subunit (LSU) rDNA has also been used to investigate fungal communities but available LSU metabarcoding primers were mostly designed to target Dikarya (Ascomycota+Basidiomycota) with little attention to early diverging fungi (EDF). However, evidence from multiple studies suggests that EDF comprise a large portion of unknown diversity in community sampling. Here, we investigate how DNA marker choice and methodological biases impact recovery of EDF from environmental samples. We focused on one EDF lineage, Zoopagomycota, as an example. We evaluated three primer sets (ITS1F/ITS2, LROR/LR3, and LR3 paired with new primer LR22F) to amplify and sequence a Zoopagomycota mock community and a set of 146 environmental samples with Illumina MiSeq. We compared two taxonomy assignment methods and created an LSU reference database compatible with AMPtk software. The two taxonomy assignment methods recovered strikingly different communities of fungi and EDF. Target fragment length variation exacerbated PCR amplification biases and influenced downstream taxonomic assignments, but this effect was greater for EDF than Dikarya. To improve identification of LSU amplicons we performed phylogenetic reconstruction and illustrate the advantages of this critical tool for investigating identified and unidentified sequences. Our results suggest much of the EDF community may be missed or misidentified with "standard" metabarcoding approaches and modified techniques are needed to understand the role of these taxa in a broader ecological context.

DNA Methylation on N6-Adenine Regulates the Hyphal Development during Dimorphism in the Early-Diverging Fungus Mucor lusitanicus.

The epigenetic modifications control the pathogenicity of human pathogenic fungi, which have been poorly studied in Mucorales, causative agents of mucormycosis. This order belongs to a group referred to as early-diverging fungi that are characterized by high levels of N6-methyldeoxy adenine (6mA) in their genome with dense 6mA clusters associated with actively expressed genes. AlkB enzymes can act as demethylases of 6mA in DNA, with the most remarkable eukaryotic examples being mammalian ALKBH1 and Caenorhabditis elegans NMAD-1. The Mucor lusitanicus (formerly M. circinelloides f. lusitanicus) genome contains one gene, dmt1, and two genes, dmt2 and dmt3, encoding proteins similar to C. elegans NMAD-1 and ALKBH1, respectively. The function of these three genes was analyzed by the generation of single and double deletion mutants for each gene. Multiple processes were studied in the mutants, but defects were only found in single and double deletion mutants for dmt1. In contrast to the wild-type strain, dmt1 mutants showed an increase in 6mA levels during the dimorphic transition, suggesting that 6mA is associated with dimorphism in M. lusitanicus. Furthermore, the spores of dmt1 mutants challenged with macrophages underwent a reduction in polar growth, suggesting that 6mA also has a role during the spore–macrophage interaction that could be important in the infection process.

The methylome of the model arbuscular mycorrhizal fungus, Rhizophagus irregularis, shares characteristics with early diverging fungi and Dikarya

Early-diverging fungi (EDF) are distinct from Dikarya and other eukaryotes, exhibiting high N6-methyldeoxyadenine (6mA) contents, rather than 5-methylcytosine (5mC). As plants transitioned to land the EDF sub-phylum, arbuscular mycorrhizal fungi (AMF; Glomeromycotina) evolved a symbiotic lifestyle with 80% of plant species worldwide. Here we show that these fungi exhibit 5mC and 6mA methylation characteristics that jointly set them apart from other fungi. The model AMF, R. irregularis, evolved very high levels of 5mC and greatly reduced levels of 6mA. However, unlike the Dikarya, 6mA in AMF occurs at symmetrical ApT motifs in genes and is associated with their transcription. 6mA is heterogeneously distributed among nuclei in these coenocytic fungi suggesting functional differences among nuclei. While far fewer genes are regulated by 6mA in the AMF genome than in EDF, most strikingly, 6mA methylation has been specifically retained in genes implicated in components of phosphate regulation; the quintessential hallmark defining this globally important symbiosis.

Terpenoid Biosynthesis Dominates among Secondary Metabolite Clusters in Mucoromycotina Genomes.

Early-diverging fungi harbour unprecedented diversity in terms of living forms, biological traits and genome architecture. Before the sequencing era, non-Dikarya fungi were considered unable to produce secondary metabolites (SM); however, this perspective is changing. The main classes of secondary metabolites in fungi include polyketides, nonribosomal peptides, terpenoids and siderophores that serve different biological roles, including iron chelation and plant growth promotion. The same classes of SM are reported for representatives of early-diverging fungal lineages. Encouraged by the advancement in the field, we carried out a systematic survey of SM in Mucoromycotina and corroborated the presence of various SM clusters (SMCs) within the phylum. Among the core findings, considerable representation of terpene and nonribosomal peptide synthetase (NRPS)-like candidate SMCs was found. Terpene clusters with diverse domain composition and potentially highly variable products dominated the landscape of candidate SMCs. A uniform low-copy distribution of siderophore clusters was observed among most assemblies. Mortierellomycotina are highlighted as the most potent SMC producers among the Mucoromycota and as a source of novel peptide products. SMC identification is dependent on gene model quality and can be successfully performed on a batch scale with genomes of different quality and completeness.