Diploid-dominant life cycles characterize the early evolution of Fungi

Phylogenomics including new sequence data of phytoplankton-infecting chytrids reveals multiple independent lifestyle transitions across the phylum

Parasitism is the most common lifestyle on Earth and has emerged many times independently across the eukaryotic tree of life. It is frequently found among chytrids (Chytridiomycota), which are early-branching unicellular fungi that feed osmotrophically via rhizoids as saprotrophs or parasites. Chytrids are abundant in most aquatic and terrestrial environments and fulfil important ecosystem functions. As parasites, they can have significant impacts on host populations. They cause global amphibian declines and influence the Earth's carbon cycle by terminating algal blooms. To date, the evolution of parasitism within the chytrid phylum remains unclear due to the low phylogenetic resolution of rRNA genes for the early diversification of fungi, and because few parasitic lineages have been cultured and genomic data for parasites is scarce. Here, we combine transcriptomics, culture-independent single-cell genomics and a phylogenomic approach to overcome these limitations. We newly sequenced 29 parasitic taxa and combined these with existing data to provide a robust backbone topology for the diversification of Chytridiomycota. Our analyses reveal multiple independent lifestyle transitions between parasitism and saprotrophy among chytrids and multiple host shifts by parasites. Based on these results and the parasitic lifestyle of other early-branching holomycotan lineages, we hypothesise that the chytrid last common ancestor was a parasite of phytoplankton.

Alternative ergosterol biosynthetic pathways confer antifungal drug resistance in the human pathogens within the Mucor species complex

Mucormycoses are emerging fungal infections caused by a variety of heterogeneous species within the Mucorales order. Among the Mucor species complex, Mucor circinelloides is the most frequently isolated pathogen in mucormycosis patients and despite its clinical significance, there is an absence of established genome manipulation techniques to conduct molecular pathogenesis studies. In this study, we generated a spontaneous uracil auxotrophic strain and developed a genetic transformation procedure to analyze molecular mechanisms conferring antifungal drug resistance. With this new model, phenotypic analyses of gene deletion mutants were conducted to define Erg3 and Erg6a as key biosynthetic enzymes in the M. circinelloides ergosterol pathway. Erg3 is a C-5 sterol desaturase involved in growth, sporulation, virulence, and azole susceptibility. In other fungal pathogens, erg3 mutations confer azole resistance because Erg3 catalyzes the production of a toxic diol upon azole exposure. Surprisingly, M. circinelloides produces only trace amounts of this toxic diol and yet, it is still susceptible to posaconazole and isavuconazole due to alterations in membrane sterol composition. These alterations are severely aggravated by erg3Δ mutations, resulting in ergosterol depletion and consequently, hypersusceptibility to azoles. We also identified Erg6a as the main C-24 sterol methyltransferase, whose activity may be partially rescued by the paralogs Erg6b and Erg6c. Loss of Erg6a function diverts ergosterol synthesis to the production of cholesta-type sterols, resulting in resistance to amphotericin B. Our findings suggest that mutations or epimutations causing loss of Erg6 function may arise during human infections, resulting in antifungal drug resistance to first-line treatments against mucormycosis.

Signatures of transposon-mediated genome inflation, host specialization, and photoentrainment in Entomophthora muscae and allied entomophthoralean fungi

Despite over a century of observations, the obligate insect parasites within the order Entomophthorales remain poorly characterized at the genetic level. In this manuscript, we present a genome for a laboratory-tractable Entomophthora muscae isolate that infects fruit flies. Our E. muscae assembly is 1.03 Gb, consists of 7810 contigs and contains 81.3% complete fungal BUSCOs. Using a comparative approach with recent datasets from entomophthoralean fungi, we show that giant genomes are the norm within Entomophthoraceae owing to extensive, but not recent, Ty3 retrotransposon activity. In addition, we find that E. muscae and its closest allies possess genes that are likely homologs to the blue-light sensor white-collar 1, a Neurospora crassa gene that has a well-established role in maintaining circadian rhythms. We uncover evidence that E. muscae diverged from other entomophthoralean fungi by expansion of existing families, rather than loss of particular domains, and possesses a potentially unique suite of secreted catabolic enzymes, consistent with E. muscae's species-specific, biotrophic lifestyle. Finally, we offer a head-to-head comparison of morphological and molecular data for species within the E. muscae species complex that support the need for taxonomic revision within this group. Altogether, we provide a genetic and molecular foundation that we hope will provide a platform for the continued study of the unique biology of entomophthoralean fungi.

Antarctica's hidden mycoviral treasures in fungi isolated from mosses: A first genomic approach

This study investigates the presence of mycoviruses in Antarctic fungi and elucidates their evolutionary relationships. To achieve this, we aligned mycoviral gene sequences with genomes of previously sequenced Antarctic endophytic fungi, made available by our research group and accessible via Joint Genome Institute. Our findings reveal that the most prevalent genetic regions in all endophytic fungi are homologous to Partitiviruses, Baculoviridae, and Phycodnaviridae. These regions display evidence of positive selection pressure, suggesting genetic diversity and the accumulation of nonsynonymous mutations. This phenomenon implies a crucial role for these regions in the adaptation and survival of these fungi in the challenging Antarctic ecosystems. The presence of mycoviruses in Antarctic endophytic fungi may indicate shared survival strategies between the virus and its host, shedding light on their evolutionary dynamics. This study underscores the significance of exploring mycoviruses within endophytic fungi and their contributions to genetic diversity. Future research avenues could delve into the functional implications of these conserved mycoviral genetic regions in Antarctic endophytic fungi, providing a comprehensive understanding of this intriguing association and genomic retention of viral region in fungi.

Coevolution of a generalist pathogen with many hosts: the case of the amphibian chytrid Batrachochytrium dendrobatidis

Generalist pathogens maintain infectivity in numerous hosts; how this broad ecological niche impacts host-pathogen coevolution remains to be widely explored. Batrachochytrium dendrobatidis (Bd) is a highly generalist pathogenic fungus that has caused devastating declines in hundreds of amphibian species worldwide. This review examines amphibian chytridiomycosis host-pathogen interactions and available evidence for coevolution between Bd and its numerous hosts. We summarize recent evidence showing that Bd genotypes vary in geographic distribution and virulence, and that amphibian species also vary in Bd susceptibility according to their geographic distribution. How much variation can be explained by phenotypic plasticity or genetic differences remains uncertain. Recent research suggests that Bd genotypes display preferences for specific hosts and that some hosts are undergoing evolution as populations rebound from Bd outbreaks. Taken together, these findings suggest the potential for coevolution to occur and illuminate a path for addressing open questions through integrating historical and contemporary genetic data.

Signatures of transposon-mediated genome inflation, host specialization, and photoentrainment in Entomophthora muscae and allied entomophthoralean fungi

Despite over a century of observations, the obligate insect parasites within the order Entomophthorales remain poorly characterized at the genetic level. This is in part due to their large genome sizes and difficulty in obtaining sequenceable material. In this manuscript, we leveraged a recently-isolated, laboratory-tractable Entomophthora muscae isolate and improved long-read sequencing to obtain a largely-complete entomophthoralean genome. Our E. muscae assembly is 1.03 Gb, consists of 7,810 contigs and contains 81.3% complete fungal BUSCOs. Using a comparative approach with other available (transcriptomic and genomic) datasets from entomophthoralean fungi, we provide new insight into the biology of these understudied pathogens. We offer a head-to-head comparison of morphological and molecular data for species within the E. muscae species complex. Our findings suggest that substantial taxonomic revision is needed to define species within this group and we provide recommendations for differentiating strains and species in the context of the existing body of E. muscae scientific literature. We show that giant genomes are the norm within Entomophthoraceae owing to extensive, but not recent, Ty3 retrotransposon activity, despite the presence of machinery to defend against transposable elements(RNAi). In addition, we find that E. muscae and its closest allies are enriched for M16A peptidases and possess genes that are likely homologs to the blue-light sensor white-collar 1, a Neurospora crassa gene that has a well-established role in maintaining circadian rhythms. We find that E. muscae has an expanded group of acid-trehalases, consistent with trehalose being the primary sugar component of fly (and insect) hemolymph. We uncover evidence that E. muscae diverged from other entomophthoralean fungi by expansion of existing families, rather than loss of particular domains, and possesses a potentially unique suite of secreted catabolic enzymes, consistent with E. muscae's species-specific, biotrophic lifestyle. Altogether, we provide a genetic and molecular foundation that we hope will provide a platform for the continued study of the unique biology of entomophthoralean fungi.

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.

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.

Potential for Use of Species in the Subfamily Erynioideae for Biological Control and Biotechnology

The fungal order Entomophthorales in the Zoopagomycota includes many fungal pathogens of arthropods. This review explores six genera in the subfamily Erynioideae within the family Entomophthoraceae, namely, Erynia, Furia, Orthomyces, Pandora, Strongwellsea, and Zoophthora. This is the largest subfamily in the Entomophthorales, including 126 described species. The species diversity, global distribution, and host range of this subfamily are summarized. Relatively few taxa are geographically widespread, and few have broad host ranges, which contrasts with many species with single reports from one location and one host species. The insect orders infected by the greatest numbers of species are the Diptera and Hemiptera. Across the subfamily, relatively few species have been cultivated in vitro, and those that have require more specialized media than many other fungi. Given their potential to attack arthropods and their position in the fungal evolutionary tree, we discuss which species might be adopted for biological control purposes or biotechnological innovations. Current challenges in the implementation of these species in biotechnology include the limited ability or difficulty in culturing many in vitro, a correlated paucity of genomic resources, and considerations regarding the host ranges of different species.

Single-cell genomics reveals new rozellid lineages and supports their sister relationship to Microsporidia

The phylum Rozellomycota has been proposed for a group of early-branching holomycotan lineages representing obligate parasites and hyperparasites of zoosporic fungi, oomycotes or phytoplankton. Given their predominantly intracellular lifestyle, rozellids are typically known from environmental ribosomal DNA data, except for the well-studied Rozella species. To date, the phylogenetic relationship between rozellids and microsporidians (Microsporidia) is not fully understood and most reliable hypotheses are based on phylogenomic analyses that incorporate the only publicly available rozellid genome of Rozella allomycis. Here, we provide genomic data of three new rozellid lineages obtained by single-cell sequencing from environmental samples and show with a phylogenomic approach that rozellids form a monophyletic group that is sister to microsporidians, corroborating the previously proposed phylum Rozellomycota. Whereas no mitochondrial genes coding for the respiratory Complex I could be found, we discovered a gene coding for a nucleotide phosphate transporter in one of the three draft genomes. The scattered absence of Complex I genes and scattered presence of nucleotide transporter genes across diverse microsporidian and rozellid lineages suggest that these adaptations to a parasitic lifestyle, which reduce the parasite's capability to synthesize ATP but enables it to steal ATP from its host, evolved independently in microsporidians and rozellids.

What can be lost? Genomic perspective on the lipid metabolism of Mucoromycota

Mucoromycota is a phylum of early diverging fungal (EDF) lineages, of mostly plant-associated terrestrial fungi. Some strains have been selected as promising biotechnological organisms due to their ability to produce polyunsaturated fatty acids and efficient conversion of nutrients into lipids. Others get their lipids from the host plant and are unable to produce even the essential ones on their own. Following the advancement in EDF genome sequencing, we carried out a systematic survey of lipid metabolism protein families across different EDF lineages. This enabled us to explore the genomic basis of the previously documented ability to produce several types of lipids within the fungal tree of life. The core lipid metabolism genes showed no significant diversity in distribution, however specialized lipid metabolic pathways differed in this regard among different fungal lineages. In total 165 out of 202 genes involved in lipid metabolism were present in all tested fungal lineages, while remaining 37 genes were found to be absent in some of fungal lineages. Duplications were observed for 69 genes. For the first time we demonstrate that ergosterol is not being produced by several independent groups of plant-associated fungi due to the losses of different ERG genes. Instead, they possess an ancestral pathway leading to the synthesis of cholesterol, which is absent in other fungal lineages. The lack of diacylglycerol kinase in both Mortierellomycotina and Blastocladiomycota opens the question on sterol equilibrium regulation in these organisms. Early diverging fungi retained most of beta oxidation components common with animals including Nudt7, Nudt12 and Nudt19 pointing at peroxisome divergence in Dikarya. Finally, Glomeromycotina and Mortierellomycotina representatives have a similar set of desaturases and elongases related to the synthesis of complex, polyunsaturated fatty acids pointing at an ancient expansion of fatty acid metabolism currently being explored by biotechnological studies.

Diversity of rhodopsin cyclases in zoospore-forming fungi

Light perception for orientation in zoospore-forming fungi is linked to homo- or heterodimeric rhodopsin-guanylyl cyclases (RGCs). Heterodimeric RGCs, first identified in the chytrid Rhizoclosmatium globosum, consist of an unusual near-infrared absorbing highly fluorescent sensitizer neorhodopsin (NeoR) that is paired with a visual light-absorbing rhodopsin responsible for enzyme activation. Here, we present a comprehensive analysis of the distribution of RGC genes in early-branching fungi using currently available genetic data. Among the characterized RGCs, we identified red-sensitive homodimeric RGC variants with maximal light activation close to 600 nm, which allow for red-light control of GTP to cGMP conversion in mammalian cells. Heterodimeric RGC complexes have evolved due to a single gene duplication within the branching of Chytridiales and show a spectral range for maximal light activation between 480 to 600 nm. In contrast, the spectral sensitivity of NeoRs is reaching into the near-infrared range with maximal absorption between 641 and 721 nm, setting the low energy spectral edge of rhodopsins so far. Based on natural NeoR variants and mutational studies, we reevaluated the role of the counterion-triad proposed to cause the extreme redshift. With the help of chimera constructs, we disclose that the cyclase domain is crucial for functioning as homo- or heterodimers, which enables the adaptation of the spectral sensitivity by modular exchange of the photosensor. The extreme spectral plasticity of retinal chromophores in native photoreceptors provides broad perspectives on the achievable spectral adaptation for rhodopsin-based molecular tools ranging from UVB into the near-infrared.

Genomes and transcriptomes help unravel the complex life cycle of the blastoclad fungus, Coelomomyces lativittatus, an obligate parasite of mosquitoes and microcrustaceans

Species of the phylum Blastocladiomycota, early-diverging zoosporic (flagellated) lineages of fungi, are vastly understudied. This phylum includes the genus Coelomomyces, which consists of more than 80 fungal species that are obligate parasites of arthropods. Known Coelomomyces species lack a complete asexual life cycle, instead surviving through an obligate heterecious alternation of generations life cycle. Despite their global distribution and interesting life cycle, little is known about the genomics of any Coelomomyces species. To address this, we generated three draft-level genomes and annotations for C. lativittatus representing its haploid meiospore, orange gamete, and amber gamete life stages. These draft genome assemblies ranged in size from 5002 to 5799 contigs, with a total length of 19.8-22.8 Mb and a mean of 7416 protein-coding genes. We then demonstrated the utility of these genomes by combining the draft annotations as a reference for analysis of C. lativittatus transcriptomes. We analyzed transcriptomes from across host-associated life stages, including infected larvae and excised mature sporangia from the mosquito Anopheles quadrimaculatus. We identified differentially expressed genes and enriched GO terms both across and within life stages and used these to make hypotheses about C. lativittatus biology. Generally, we found the C. lativittatus transcriptome to be a complex and dynamic expression landscape; GO terms related to metabolism and transport processes were enriched during infection and terms related to dispersal were enriched during sporulation. We further identified five high mobility group (HMG)-box genes in C. lativittatus, three belonging to clades with mating type (MAT) loci from other fungi, as well as four ortholog expansions in C. lativittatus compared with other fungi. The C. lativittatus genomes and transcriptomes reported here are a valuable resource and may be leveraged toward furthering understanding of the biology of these and other early-diverging fungal lineages.

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.

Diploid Nuclei Occur throughout the Life Cycles of Pucciniales Fungi

Within Eukaryotes, fungi are the typical representatives of haplontic life cycles. Basidiomycota fungi are dikaryotic in extensive parts of their life cycle, but diploid nuclei are known to form only in basidia. Among Basidiomycota, the Pucciniales are notorious for presenting the most complex life cycles, with high host specialization, and for their expanded genomes. Using cytogenomic (flow cytometry and cell sorting on propidium iodide-stained nuclei) and cytogenetic (FISH with rDNA probe) approaches, we report the widespread occurrence of replicating haploid and diploid nuclei (i.e., 1C, 2C and a small proportion of 4C nuclei) in diverse life cycle stages (pycnial, aecial, uredinial, and telial) of all 35 Pucciniales species analyzed, but not in sister taxa. These results suggest that the Pucciniales life cycle is distinct from any cycle known, i.e., neither haplontic, diplontic nor haplodiplontic, corroborating patchy and disregarded previous evidence. However, the biological basis and significance of this phenomenon remain undisclosed. IMPORTANCE Within Eukaryotes, fungi are the typical representatives of haplontic life cycles, contrasting with plants and animals. As such, fungi thus contain haploid nuclei throughout their life cycles, with sexual reproduction generating a single diploid cell upon karyogamy that immediately undergoes meiosis, thus resuming the haploid cycle. In this work, using cytogenetic and cytogenomic tools, we demonstrate that a vast group of fungi presents diploid nuclei throughout their life cycles, along with haploid nuclei, and that both types of nuclei replicate. Moreover, haploid nuclei are absent from urediniospores. The phenomenon appears to be transversal to the organisms in the order Pucciniales (rust fungi) and it does not occur in neighboring taxa, but a biological explanation or function for it remains elusive.

Genomes of fungi and relatives reveal delayed loss of ancestral gene families and evolution of key fungal traits

Fungi are ecologically important heterotrophs that have radiated into most niches on Earth and fulfil key ecological services. Despite intense interest in their origins, major genomic trends of their evolutionary route from a unicellular opisthokont ancestor to derived multicellular fungi remain poorly known. Here we provide a highly resolved genome-wide catalogue of gene family changes across fungal evolution inferred from the genomes of 123 fungi and relatives. We show that a dominant trend in early fungal evolution has been the gradual shedding of protist genes and the punctuated emergence of innovation by two main gene duplication events. We find that the gene content of non-Dikarya fungi resembles that of unicellular opisthokonts in many respects, owing to the conservation of protist genes in their genomes. The most rapidly duplicating gene groups included extracellular proteins and transcription factors, as well as ones linked to the coordination of nutrient uptake with growth, highlighting the transition to a sessile osmotrophic feeding strategy and subsequent lifestyle evolution as important elements of early fungal history. These results suggest that the genomes of pre-fungal ancestors evolved into the typical filamentous fungal genome by a combination of gradual gene loss, turnover and several large duplication events rather than by abrupt changes. Consequently, the taxonomically defined Fungi represents a genomically non-uniform assemblage of species.

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.

Tubulin Polymerization Promoting Proteins (TPPPs) of Aphelidiomycota: Correlation between the Incidence of p25alpha Domain and the Eukaryotic Flagellum

The seven most early diverging lineages of the 18 phyla of fungi are the non-terrestrial fungi, which reproduce through motile flagellated zoospores. There are genes/proteins that are present only in organisms with flagellum or cilium. It was suggested that TPPP-like proteins (proteins containing at least one complete or partial p25alpha domain) are among them, and a correlation between the incidence of the p25alpha domain and the eukaryotic flagellum was hypothesized. Of the seven phyla of flagellated fungi, six have been known to contain TPPP-like proteins. Aphelidiomycota, one of the early-branching phyla, has some species (e.g., Paraphelidium tribonematis) that retain the flagellum, whereas the Amoeboaphelidium genus has lost the flagellum. The first two Aphelidiomycota genomes (Amoeboaphelidium protococcorum and Amoeboaphelidium occidentale) were sequenced and published last year. A BLASTP search revealed that A. occidentale does not have a TPPP, but A. protococcorum, which possesses pseudocilium, does have a TPPP. This TPPP is the ‘long-type’ which occurs mostly in animals as well as other Opisthokonta. P. tribonematis has a ‘fungal-type’ TPPP, which is found only in some flagellated fungi. These data on Aphelidiomycota TPPP proteins strengthen the correlation between the incidence of p25alpha domain-containing proteins and that of the eukaryotic flagellum/cilium.

On the evolution of variation in sexual reproduction through the prism of eukaryotic microbes

Almost all eukaryotes undergo sexual reproduction to generate diversity and select for fitness in their population pools. Interestingly, the systems by which sex is defined are highly diverse and can even differ between evolutionarily closely related species. While the most commonly known form of sex determination involves males and females in animals, eukaryotic microbes can have as many as thousands of different mating types for the same species. Furthermore, some species have found alternatives to sexual reproduction and prefer to grow clonally and yet undergo infrequent facultative sexual reproduction. These organisms are mainly invertebrates and microbes, but several examples are also present among vertebrates suggesting that alternative modes of sexual reproduction evolved multiple times throughout evolution. In this review, we summarize the sex-determination modes and variants of sexual reproduction found across the eukaryotic tree of life and suggest that eukaryotic microbes provide unique opportunities to study these processes in detail. We propose that understanding variations in modes of sexual reproduction can serve as a foundation to study the evolution of sex and why and how it evolved in the first place.

Mycoparasites, Gut Dwellers, and Saprotrophs: Phylogenomic Reconstructions and Comparative Analyses of Kickxellomycotina Fungi

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.

New insights into the genome and transmission of the microsporidian pathogen Nosema muscidifuracis

Introduction

Nosema is a diverse genus of unicellular microsporidian parasites of insects and other arthropods. Nosema muscidifuracis infects parasitoid wasp species of Muscidifurax zaraptor and M. raptor (Hymenoptera: Pteromalidae), causing ~50% reduction in longevity and ~90% reduction in fecundity.

Methods and Results

Here, we report the first assembly of the N. muscidifuracis genome (14,397,169 bp in 28 contigs) of high continuity (contig N50 544.3 Kb) and completeness (BUSCO score 97.0%). A total of 2,782 protein-coding genes were annotated, with 66.2% of the genes having two copies and 24.0% of genes having three copies. These duplicated genes are highly similar, with a sequence identity of 99.3%. The complex pattern suggests extensive gene duplications and rearrangements across the genome. We annotated 57 rDNA loci, which are highly GC-rich (37%) in a GC-poor genome (25% genome average). Nosema-specific qPCR primer sets were designed based on 18S rDNA annotation as a diagnostic tool to determine its titer in host samples. We discovered high Nosema titers in Nosema-cured M. raptor and M. zaraptor using heat treatment in 2017 and 2019, suggesting that the remedy did not completely eliminate the Nosema infection. Cytogenetic analyses revealed heavy infections of N. muscidifuracis within the ovaries of M. raptor and M. zaraptor, consistent with the titer determined by qPCR and suggesting a heritable component of infection and per ovum vertical transmission.

Discussion

The parasitoids-Nosema system is laboratory tractable and, therefore, can serve as a model to inform future genome manipulations of Nosema-host system for investigations of Nosemosis.

A Genome Sequence Assembly of the Phototactic and Optogenetic Model Fungus Blastocladiella emersonii Reveals a Diversified Nucleotide-Cyclase Repertoire

The chytrid fungus Blastocladiella emersonii produces spores with swimming tails (zoospores); these cells can sense and swim toward light. Interest in this species stems from ongoing efforts to develop B. emersonii as a model for understanding the evolution of phototaxis and the molecular cell biology of the associated optogenetic circuits. Here, we report a highly contiguous genome assembly and gene annotation of the B. emersonii American Type Culture Collection 22665 strain. We integrate a PacBio long-read library with an Illumina paired-end genomic sequence survey leading to an assembly of 21 contigs totaling 34.27 Mb. Using these data, we assess the diversity of sensory system encoding genes. These analyses identify a rich complement of G-protein-coupled receptors, ion transporters, and nucleotide cyclases, all of which have been diversified by domain recombination and tandem duplication. In many cases, these domain combinations have led to the fusion of a protein domain to a transmembrane domain, tying a putative signaling function to the cell membrane. This pattern is consistent with the diversification of the B. emersonii sensory-signaling systems, which likely plays a varied role in the complex life cycle of this fungus.