Divergent Evolution of Early Terrestrial Fungi Reveals the Evolution of Mucormycosis Pathogenicity Factors.

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.

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.

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.

Phylogenetic and morphological analyses of the mycoparasitic genus Piptocephalis.

The Piptocephalidaceae (Zoopagales, Zoopagomycota) contains three genera of mycoparasitic, haustoria-forming fungi: Kuzuhaea, Piptocephalis, and Syncephalis. Although the species in this family are diverse and ubiquitous in soil and dung, they are among the least studied fungi. Co-cultures of Piptocephalis and their hosts are relatively easy to isolate from soil and dung samples across the globe, making them a good model taxon for the order Zoopagales. This study focuses on the systematics of the genus Piptocephalis. Despite the fact that there are approximately 40 described Piptocephalis species, there are no modern taxonomic or molecular phylogenetic treatments of this group. Minimal sequence data are available, and relatively little is known about the true diversity or biogeography of the genus. Our study addresses two aspects: Piptocephalis systematics and analyses of the length and inter- and infraspecific variation of the nuc rDNA internal transcribed spacer (ITS1-5.8S-ITS2 = ITS) region. First, we generated a large subunit (28S) nuc rDNA phylogeny and evaluated several morphological characters by testing their correlation with the phylogeny using Bayesian Tip-association Significance testing (BaTS). We found monophyly of Piptocephalis species identified based on morphological traits, but morphological character states were not conserved across clades, suggesting that there have been multiple gains and losses of morphological characters. We also found that Kuzhuaea is nested within Piptocephalis. Second, we amplified the ITS from many Piptocephalis isolates, created a sequence alignment, and measured the lengths using the software ITSx. Piptocephalis species had ITS regions that were longer than the average for most Dikarya but were similar in length to those of the related genus Syncephalis.

Leveraging single-cell genomics to expand the fungal tree of life.

Environmental DNA surveys reveal that most fungal diversity represents uncultured species. We sequenced the genomes of eight uncultured species across the fungal tree of life using a new single-cell genomics pipeline. We show that, despite a large variation in genome and gene space recovery from each single amplified genome (SAG), ≥90% can be recovered by combining multiple SAGs. SAGs provide robust placement for early-diverging lineages and infer a diploid ancestor of fungi. Early-diverging fungi share metabolic deficiencies and show unique gene expansions correlated with parasitism and unculturability. Single-cell genomics holds great promise in exploring fungal diversity, life cycles and metabolic potential.

Resolving relationships at the animal-fungal divergence: A molecular phylogenetic study of the protist trichomycetes (Ichthyosporea, Eccrinida).

Trichomycetes is a group of microorganisms that was considered a class of fungi comprising four orders of commensal, gut-dwelling endosymbionts obligately associated with arthropods. Since molecular phylogenies revealed two of those orders (Amoebidiales and Eccrinales="protist trichos") to be closely related to members of the protist class Ichthyosporea (=Mesomycetozoea), trichomycetes have been considered an ecological association of both early-diverging fungi and protists. Understanding of the taxonomy, evolution, and diversity of the protist trichos is lacking largely due to the difficulties inherent in species collection that have contributed to undersampling and understudy. The most recent classification divides the protist trichos between two families, Amoebidiidae and Eccrinidae (suborder Trichomycina, order Eccrinida). However, there is no comprehensive molecular phylogeny available for this group and major questions about the systematics of protist trichos remain unanswered. Therefore, we generated 18S and 28S rDNA sequences for 106 protist tricho samples and combined them with publicly available Eccrinida sequences for phylogenetic analyses. We also sequenced a conserved protein-coding gene (heat-shock 70 protein) to obtain a multigene data set. We conducted ancestral state reconstruction (ASR) and Bayesian tip-association significance test (BaTS) analyses by mapping six morphological and ecological characters onto the resulting phylogenetic trees. Our results demonstrate: (1) several ecological and morphological character states (habitat, host type, host stage at time of infestation, location within host, spore production, and growth form) are significantly correlated with the phylogeny, and (2) two additional protist tricho families should be incorporated into the taxonomy to reflect phylogenetic relationships. Our data suggest that an integrated strategy that combines morphological, ecological, and molecular characters is needed to further resolve and clarify the systematics of the Eccrinida.

Mutations that rescue the paralysis of Caenorhabditis elegans ric-8 (synembryn) mutants activate the G alpha(s) pathway and define a third major branch of the synaptic signaling network.

To identify hypothesized missing components of the synaptic G alpha(o)-G alpha(q) signaling network, which tightly regulates neurotransmitter release, we undertook two large forward genetic screens in the model organism C. elegans and focused first on mutations that strongly rescue the paralysis of ric-8(md303) reduction-of-function mutants, previously shown to be defective in G alpha(q) pathway activation. Through high-resolution mapping followed by sequence analysis, we show that these mutations affect four genes. Two activate the G alpha(q) pathway through gain-of-function mutations in G alpha(q); however, all of the remaining mutations activate components of the G alpha(s) pathway, including G alpha(s), adenylyl cyclase, and protein kinase A. Pharmacological assays suggest that the G alpha(s) pathway-activating mutations increase steady-state neurotransmitter release, and the strongly impaired neurotransmitter release of ric-8(md303) mutants is rescued to greater than wild-type levels by the strongest G alpha(s) pathway activating mutations. Using transgene induction studies, we show that activating the G alpha(s) pathway in adult animals rapidly induces hyperactive locomotion and rapidly rescues the paralysis of the ric-8 mutant. Using cell-specific promoters we show that neuronal, but not muscle, G alpha(s) pathway activation is sufficient to rescue ric-8(md303)'s paralysis. Our results appear to link RIC-8 (synembryn) and a third major G alpha pathway, the G alpha(s) pathway, with the previously discovered G alpha(o) and G alpha(q) pathways of the synaptic signaling network.

Convergent, RIC-8-dependent Galpha signaling pathways in the Caenorhabditis elegans synaptic signaling network.

We used gain-of-function and null synaptic signaling network mutants to investigate the relationship of the G alpha(q) and G alpha(s) pathways to synaptic vesicle priming and to each other. Genetic epistasis studies using G alpha(q) gain-of-function and null mutations, along with a mutation that blocks synaptic vesicle priming and the synaptic vesicle priming stimulator phorbol ester, suggest that the G alpha(q) pathway generates the core, obligatory signals for synaptic vesicle priming. In contrast, the G alpha(s) pathway is not required for the core priming function, because steady-state levels of neurotransmitter release are not significantly altered in animals lacking a neuronal G alpha(s) pathway, even though these animals are strongly paralyzed as a result of functional (nondevelopmental) defects. However, our genetic analysis indicates that these two functionally distinct pathways converge and that they do so downstream of DAG production. Further linking the two pathways, our epistasis analysis of a ric-8 null mutant suggests that RIC-8 (a receptor-independent G alpha guanine nucleotide exchange factor) is required to maintain both the G alpha(q) vesicle priming pathway and the neuronal G alpha(s) pathway in a functional state. We propose that the neuronal G alpha(s) pathway transduces critical positional information onto the core G alpha(q) pathway to stabilize the priming of selected synapses that are optimal for locomotion.