Search

Endofungal Mycetohabitans (formerly Burkholderia) spp. rely on a type III secretion system to deliver mostly unidentified effector proteins when colonizing their host fungus, Rhizopus microsporus. The one known secreted effector family from Mycetohabitans consists of homologues of transcription activator-like (TAL) effectors, which are used by plant pathogenic Xanthomonas and Ralstonia spp. to activate host genes that promote disease. These 'Burkholderia TAL-like (Btl)' proteins bind corresponding specific DNA sequences in a predictable manner, but their genomic target(s) and impact on transcription in the fungus are unknown. Recent phenotyping of Btl mutants of two Mycetohabitans strains revealed that the single Btl in one Mycetohabitans endofungorum strain enhances fungal membrane stress tolerance, while others in a Mycetohabitans rhizoxinica strain promote bacterial colonization of the fungus. The phenotypic diversity underscores the need to assess the sequence diversity and, given that sequence diversity translates to DNA targeting specificity, the functional diversity of Btl proteins. Using a dual approach to maximize capture of Btl protein sequences for our analysis, we sequenced and assembled nine Mycetohabitans spp. genomes using long-read PacBio technology and also mined available short-read Illumina fungal-bacterial metagenomes. We show that btl genes are present across diverse Mycetohabitans strains from Mucoromycota fungal hosts yet vary in sequences and predicted DNA binding specificity. Phylogenetic analysis revealed distinct clades of Btl proteins and suggested that Mycetohabitans might contain more species than previously recognized. Within our data set, Btl proteins were more conserved across M. rhizoxinica strains than across M. endofungorum, but there was also evidence of greater overall strain diversity within the latter clade. Overall, the results suggest that Btl proteins contribute to bacterial-fungal symbioses in myriad ways.

Every fungal cell is encapsulated in a cell wall, essential for cell viability, morphogenesis, and pathogenesis. Most knowledge of the cell wall composition in fungi has focused on ascomycetes, especially human pathogens, but considerably less is known about early divergent fungal groups, such as species in the Zoopagomycota and Mucoromycota phyla. To shed light on evolutionary changes in the fungal cell wall, we studied the monosaccharide composition of the cell wall of 18 species including early diverging fungi and species in the Basidiomycota and Ascomycota phyla with a focus on those with pathogenic lifestyles and interactions with plants. Our data revealed that chitin is the most characteristic component of the fungal cell wall, and was found to be in a higher proportion in the early divergent groups. The Mucoromycota species possess few glucans, but instead have other monosaccharides such as fucose and glucuronic acid that are almost exclusively found in their cell walls. Additionally, we observed that hexoses (glucose, mannose and galactose) accumulate in much higher proportions in species belonging to Dikarya. Our data demonstrate a clear relationship between phylogenetic position and fungal cell wall carbohydrate composition and lay the foundation for a better understanding of their evolution and their role in plant interactions.

The global panzootic lineage (GPL) of the pathogenic fungus Batrachochytrium dendrobatidis (Bd) has caused severe amphibian population declines, yet the drivers underlying the high frequency of GPL in regions of amphibian decline are unclear. Using publicly available Bd genome sequences, we identified multiple non-GPL Bd isolates that contain a circular Rep-encoding single-stranded (CRESS)-like DNA virus, which we named Bd DNA virus 1 (BdDV-1). We further sequenced and constructed genome assemblies with long read sequences to find that the virus is integrated into the nuclear genome in some strains. Attempts to cure virus-positive isolates were unsuccessful; however, phenotypic differences between naturally virus-positive and virus-negative Bd isolates suggested that BdDV-1 decreases the growth of its host in vitro but increases the virulence of its host in vivo. BdDV-1 is the first-described CRESS DNA mycovirus of zoosporic true fungi, with a distribution inversely associated with the emergence of the panzootic lineage.

Reptiles and amphibians (herptiles) are some of the most endangered and threatened species on the planet and numerous conservation strategies are being implemented with the goal of ensuring species recovery. Little is known, however, about the gut microbiome of wild herptiles and how it relates to the health of these populations. Here, we report results from the gut microbiome characterization of both a broad survey of herptiles, and the correlation between the fungus Basidiobolus, and the bacterial community supported by a deeper, more intensive sampling of Plethodon glutinosus, known as slimy salamanders. We demonstrate that bacterial communities sampled from frogs, lizards, and salamanders are structured by the host taxonomy and that Basidiobolus is a common and natural component of these wild gut microbiomes. Intensive sampling of multiple hosts across the ecoregions of Tennessee revealed that geography and host:geography interactions are strong predictors of distinct Basidiobolus operational taxonomic units present within a given host. Co-occurrence analyses of Basidiobolus and bacterial community diversity support a correlation and interaction between Basidiobolus and bacteria, suggesting that Basidiobolus may play a role in structuring the bacterial community. We further the hypothesis that this interaction is advanced by unique specialized metabolism originating from horizontal gene transfer from bacteria to Basidiobolus and demonstrate that Basidiobolus is capable of producing a diversity of specialized metabolites including small cyclic peptides.IMPORTANCEThis work significantly advances our understanding of biodiversity and microbial interactions in herptile microbiomes, the role that fungi play as a structural and functional members of herptile gut microbiomes, and the chemical functions that structure microbiome phenotypes. We also provide an important observational system of how the gut microbiome represents a unique environment that selects for novel metabolic functions through horizontal gene transfer between fungi and bacteria. Such studies are needed to better understand the complexity of gut microbiomes in nature and will inform conservation strategies for threatened species of herpetofauna.

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.

The entomopathogenic fungus Conoideocrella luteorostrata has recently been implicated in natural epizootics among exotic elongate hemlock scale (EHS) insects in Christmas tree farms in the eastern United States. Since 1913, C. luteorostrata has been reported from various plant feeding Hemiptera in the southeastern United States, but comprehensive morphological and phylogenetic studies of U.S. populations are lacking. The recovery of multiple strains of C. luteorostrata from mycosed EHS in North Carolina provided an opportunity to conduct pathogenicity assays and morphological and phylogenetic studies to investigate genus- and species-level boundaries among the Clavicipitaceae. Pathogenicity assays confirmed that C. luteorostrata causes mortality of EHS crawlers, an essential first step in developing this fungus as a biocontrol. Morphological studies revealed that conidia aligned with previous measurements of the Paecilomyces-like asexual state of C. luteorostrata, with conidiophore morphology consistent with historical observations. Additionally, a Hirsutella-like synanamorph was observed in select C. luteorostrata strains. In both a four-locus, 54-taxon Clavicipitaceae-wide phylogenetic analysis including D1-D2 domains of the nuclear 28S rRNA region (28S), elongation factor 1 alpha (EF1-α), DNA-directed RNA polymerase II subunit 1 (RPB1), and DNA-directed RNA polymerase II subunit 2 (RPB2) and a two-locus, 38-taxon (28S and EF1-α) phylogenetic analysis, all three Conoideocrella species were resolved as strongly supported monophyletic lineages across all loci and both methods (maximum likelihood and Bayesian inference) of phylogenetic inference except for 28S for C. tenuis. Despite the strong support for individual Conoideocrella species, none of the analyses supported the monophyly of Conoideocrella with the inclusion of Dussiella. Due to the paucity of RPB1 and RPB2 sequence data, EF1-α provided superior delimitation of intraspecies groupings for Conoideocrella and should be used in future studies. Further development of C. luteorostrata as a biocontrol against EHS will require additional surveys across diverse Hemiptera and expanded pathogenicity testing to clarify host range and efficacy of this fungus.

Survival and growth strategies of Antarctic endolithic microbes residing in Earth's driest and coldest desert remain virtually unknown. From 109 endolithic microbiomes, 4539 metagenome-assembled genomes were generated, 49.3 % of which were novel candidate bacterial species. We present evidence that trace gas oxidation and atmospheric chemosynthesis may be the prevalent strategies supporting metabolic activity and persistence of these ecosystems at the fringe of life and the limits of habitability.

A 30.28 Mb draft genome sequence was assembled and annotated for the melanized ascomycetous fungus Exophiala xenobiotica NRRL_64630 (Pezizomycotina; Chaetothyriales) isolated from La Brea Tar Pits, Los Angeles, California. Species identification was made by phylogenetic assessment of the Internal Transcribed Spacer. This is the first isolated fungal species from this historic space.

We present six whole community shotgun metagenomic sequencing data sets of two types of biological soil crusts sampled at the ecotone of the Mojave Desert and Colorado Desert in California. These data will help us understand the diversity and function of biocrust microbial communities, which are essential for desert ecosystems.

We present eight metatranscriptomic datasets of light algal and cyanolichen biological soil crusts from the Mojave Desert in response to wetting. These data will help us understand gene expression patterns in desert biocrust microbial communities after they have been reactivated by the addition of water.

Members of the family Trichomeriaceae, belonging to the Chaetothyriales order and the Ascomycota phylum, are known for their capability to inhabit hostile environments characterized by extreme temperatures, oligotrophic conditions, drought, or presence of toxic compounds. The genus Knufia encompasses many polyextremophilic species. In this report, the genomic and morphological features of the strain FJI-L2-BK-P2 presented, which was isolated from the Mars 2020 mission spacecraft assembly facility located at the Jet Propulsion Laboratory in Pasadena, California. The identification is based on sequence alignment for marker genes, multi-locus sequence analysis, and whole genome sequence phylogeny. The morphological features were studied using a diverse range of microscopic techniques (bright field, phase contrast, differential interference contrast and scanning electron microscopy). The phylogenetic marker genes of the strain FJI-L2-BK-P2 exhibited highest similarities with type strain of Knufia obscura (CBS 148926T) that was isolated from the gas tank of a car in Italy. To validate the species identity, whole genomes of both strains (FJI-L2-BK-P2 and CBS 148926T) were sequenced, annotated, and strain FJI-L2-BK-P2 was confirmed as K. obscura. The morphological analysis and description of the genomic characteristics of K. obscura FJI-L2-BK-P2 may contribute to refining the taxonomy of Knufia species. Key morphological features are reported in this K. obscura strain, resembling microsclerotia and chlamydospore-like propagules. These features known to be characteristic features in black fungi which could potentially facilitate their adaptation to harsh environments.

Giant viruses of the phylum Nucleocytoviricota have captured researchers’ attention due to their increasingly recognized impacts on eukaryotic genome evolution. Their origins are hypothesized to predate or coincide with the diversification of eukaryotes and their natural hosts span the eukaryotic tree of life. But surprisingly, giant viruses have not been found in the fungal kingdom, though recent metagenomic work suggests a putative association. Here we show that early-diverging fungal lineages maintain both “viral fossils” and host active infections by giant viruses, which form a monophyletic clade that we name Mycodnaviridae . Viral genomes described here span up to 350kb and encode over 300 genes, including many genes characteristic of giant viruses. Interestingly, we observed variation in infection status among the isolates including apparent active infection and transcriptionally-suppressed states, suggestive that viral activation may be constrained to certain life stages of the host. Our experimental findings add to the scant few natural virus-host systems available in culture for the study of giant viruses. These viruses have likely shaped the early evolution of these fungal lineages and should prove a useful model for their study.

Acarospora socialis, the bright cobblestone lichen, is commonly found in southwestern North America. This charismatic yellow lichen is a species of key ecological significance as it is often a pioneer species in new environments. Despite their ecological importance virtually no research has been conducted on the genomics of A. socialis. To address this, we used long-read sequencing to generate the first high-quality draft genome of A. socialis. Lichen thallus tissue was collected from Pinkham Canyon in Joshua Tree National Park, California and deposited in the UC Riverside herbarium under accession #295874. The de novo assembly of the mycobiont partner of the lichen was generated from Pacific Biosciences HiFi long reads and Dovetail Omni-C chromatin capture data. After removing algal and bacterial contigs, the fungal genome was approximately 31.2 Mb consisting of 38 scaffolds with contig and scaffold N50 of 2.4 Mb. The BUSCO completeness score of the assembled genome was 97.5% using the Ascomycota gene set. Information on the genome of A. socialis is important for California conservation purposes given that this lichen is threatened in some places locally by wildfires due to climate change. This reference genome will be used for understanding the genetic diversity, population genomics, and comparative genomics of A. socialis species. Genomic resources for this species will support population and landscape genomics investigations, exploring the use of A. socialis as a bioindicator species for climate change, and in studies of adaptation by comparing populations that occur across aridity gradients in California.

Nutritional deprivation triggers a switch from a saprotrophic to predatory lifestyle in soil-dwelling nematode-trapping fungi (NTF). In particular, the NTF Arthrobotrys oligospora secretes food and sex cues to lure nematodes to its mycelium and is triggered to develop specialized trapping devices. Captured nematodes are then invaded and digested by the fungus, thus serving as a food source. In this study, we examined the transcriptomic response of A. oligospora across the stages of sensing, trap development, and digestion upon exposure to the model nematode Caenorhabditis elegans. A. oligospora enacts a dynamic transcriptomic response, especially of protein secretion-related genes, in the presence of prey. Two-thirds of the predicted secretome of A. oligospora was up-regulated in the presence of C. elegans at all time points examined, and among these secreted proteins, 38.5% are predicted to be effector proteins. Furthermore, functional studies disrupting the t-SNARE protein Sso2 resulted in impaired ability to capture nematodes. Additionally, genes of the DUF3129 family, which are expanded in the genomes of several NTF, were highly up-regulated upon nematode exposure. We observed the accumulation of highly expressed DUF3129 proteins in trap cells, leading us to name members of this gene family as Trap Enriched Proteins (TEPs). Gene deletion of the most highly expressed TEP gene, TEP1, impairs the function of traps and prevents the fungus from capturing prey efficiently. In late stages of predation, we observed up-regulation of a variety of proteases, including metalloproteases. Following penetration of nematodes, these metalloproteases facilitate hyphal growth required for colonization of prey. These findings provide insights into the biology of the predatory lifestyle switch in a carnivorous fungus and provide frameworks for other fungal-nematode predator-prey systems.

The glassy-winged sharpshooter, Homalodisca vitripennis Germar, is an invasive xylem-feeding leafhopper with a devastating economic impact on California agriculture through transmission of the plant pathogen, Xylella fastidiosa. While studies have focused on X. fastidiosa or known symbionts of H. vitripennis, little work has been done at the scale of the microbiome (the bacterial community) or mycobiome (the fungal community). Here, we characterize the mycobiome and the microbiome of H. vitripennis across Southern California and explore correlations with captivity and host insecticide resistance status. Using high-throughput sequencing of the ribosomal internal transcribed spacer 1 region and the 16S rRNA gene to profile the mycobiome and microbiome, respectively, we found that while the H. vitripennis mycobiome significantly varied across Southern California, the microbiome did not. We also observed a significant difference in both the mycobiome and microbiome between captive and wild H. vitripennis. Finally, we found that the mycobiome, but not the microbiome, was correlated with insecticide resistance status in wild H. vitripennis. This study serves as a foundational look at the H. vitripennis mycobiome and microbiome across Southern California. Future work should explore the putative link between microbes and insecticide resistance status and investigate whether microbial communities should be considered in H. vitripennis management practices. IMPORTANCE The glassy-winged sharpshooter is an invasive leafhopper that feeds on the xylem of plants and transmits the devastating pathogen, Xylella fastidiosa, resulting in significant economic damage to California's agricultural system. While studies have focused on this pathogen or obligate symbionts of the glassy-winged sharpshooter, there is limited knowledge of the bacterial and fungal communities that make up its microbiome and mycobiome. To address this knowledge gap, we explored the composition of the mycobiome and the microbiome of the glassy-winged sharpshooter across Southern California and identified differences associated with geography, captivity, and host insecticide resistance status. Understanding sources of variation in the microbial communities associated with the glassy-winged sharpshooter is an important consideration for developing management strategies to control this invasive insect. This study is a first step toward understanding the role microbes may play in the glassy-winged sharpshooter's resistance to insecticides.

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.

The genetic disease cystic fibrosis (CF) frequently leads to chronic lung infections by bacteria and fungi. We identified three individuals with CF with persistent lung infections dominated by Clavispora (Candida) lusitaniae. Whole-genome sequencing analysis of multiple isolates from each infection found evidence for selection for mutants in the gene MRS4 in all three distinct lung-associated populations. In each population, we found one or two unfixed, non-synonymous mutations in MRS4 relative to the reference allele found in multiple environmental and clinical isolates including the type strain. Genetic and phenotypic analyses found that all evolved alleles led to loss of function (LOF) of Mrs4, a mitochondrial iron transporter. RNA-seq analyses found that Mrs4 variants with decreased activity led to increased expression of genes involved in iron acquisition mechanisms in both low iron and replete iron conditions. Furthermore, surface iron reductase activity and intracellular iron were much higher in strains with Mrs4 LOF variants. Parallel studies found that a subpopulation of a CF-associated Exophiala dermatitidis infection also had a non-synonymous LOF mutation in MRS4. Together, these data suggest that MRS4 mutations may be beneficial during chronic CF lung infections in diverse fungi, perhaps, for the purposes of adaptation to an iron-restricted environment with chronic infections. IMPORTANCE The identification of MRS4 mutations in Clavispora (Candida) lusitaniae and Exophiala dermatitidis in individuals with cystic fibrosis (CF) highlights a possible adaptive mechanism for fungi during chronic CF lung infections. The findings of this study suggest that loss of function of the mitochondrial iron transporter Mrs4 can lead to increased activity of iron acquisition mechanisms, which may be advantageous for fungi in iron-restricted environments during chronic infections. This study provides valuable information for researchers working toward a better understanding of the pathogenesis of chronic lung infections and more effective therapies to treat them.

Individuals with cystic fibrosis (CF) are susceptible to chronic lung infections that lead to inflammation and irreversible lung damage. While most respiratory infections that occur in CF are caused by bacteria, some are dominated by fungi such as the slow-growing black yeast Exophiala dermatitidis. Here, we analyze isolates of E. dermatitidis cultured from two samples, collected from a single subject 2 years apart. One isolate genome was sequenced using long-read Nanopore technology as an in-population reference to use in comparative single nucleotide polymorphism and insertion-deletion variant analyses of 23 isolates. We then used population genomics and phylo-genomics to compare the isolates to each other as well as the reference genome strain E. dermatitidis NIH/UT8656. Within the CF lung population, three E. dermatitidis clades were detected, each with varying mutation rates. Overall, the isolates were highly similar suggesting that they were recently diverged. All isolates were MAT 1-1, which was consistent with their high relatedness and the absence of evidence for mating or recombination between isolates. Phylogenetic analysis grouped sets of isolates into clades that contained isolates from both early and late time points indicating there are multiple persistent lineages. Functional assessment of variants unique to each clade identified alleles in genes that encode transporters, cytochrome P450 oxidoreductases, iron acquisition, and DNA repair processes. Consistent with the genomic heterogeneity, isolates showed some stable phenotype heterogeneity in melanin production, subtle differences in antifungal minimum inhibitory concentrations, and growth on different substrates. The persistent population heterogeneity identified in lung-derived isolates is an important factor to consider in the study of chronic fungal infections, and the analysis of changes in fungal pathogens over time may provide important insights into the physiology of black yeasts and other slow-growing fungi in vivo.

During the construction and assembly of the Mars 2020 mission components at two different NASA cleanrooms, several fungal strains were isolated. Based on their colony morphology, two strains that showed yeast-like appearance were further characterized for their phylogenetic position. The species-level classification of these two novel strains, using traditional colony and cell morphology methods combined with the phylogenetic reconstructions using multi-locus sequence analysis (MLSA) based on several gene loci (ITS, LSU, SSU, RPB1, RPB2, CYTB and TEF1), and whole genome sequencing (WGS) was carried out. This polyphasic taxonomic approach supported the conclusion that the two basidiomycetous yeasts belong to hitherto undescribed species. The strain FJI-L2-BK-P3T, isolated from the Jet Propulsion Laboratory Spacecraft Assembly Facility, was placed in the Naganishia albida clade (Filobasidiales, Tremellomycetes), but is genetically and physiologically different from other members of the clade. Another yeast strain FKI-L6-BK-PAB1T, isolated from the Kennedy Space Center Payload Hazardous and Servicing Facility, was placed in the genus Cystobasidium (Cystobasidiales, Cystobasidiomycetes) and is distantly related to C. benthicum. Here we propose two novel species with the type strains, Naganishia kalamii sp. nov. (FJI-L2-BK-P3T = NRRL 64466 = DSM 115730) and Cystobasidium onofrii sp. nov. (FKI-L6-BK-PAB1T = NRRL 64426 = DSM 114625). The phylogenetic analyses revealed that single gene phylogenies (ITS or LSU) were not conclusive, and MLSA and WGS-based phylogenies were more advantageous for species discrimination in the two genera. The genomic analysis predicted proteins associated with dehydration and desiccation stress-response and the presence of genes that are directly related to osmotolerance and psychrotolerance in both novel yeasts described. Cells of these two newly-described yeasts were exposed to UV-C radiation and compared with N. onofrii, an extremophilic UV-C resistant cold-adapted Alpine yeast. Both novel species were UV resistant, emphasizing the need for collecting and characterizing extremotolerant microbes, including yeasts, to improve microbial reduction techniques used in NASA planetary protection programs.

Lichtheimia ornata is an emerging opportunistic Mucorales pathogen that is associated with fatal infections in immunocompromised individuals. While these environmentally acquired infections have rarely been reported to date, cases were noted in a recent analysis of coronavirus disease 2019 (COVID-19)-associated mucormycosis in India. Here, we report the annotated genome sequence of the environmental isolate CBS 291.66.

There are few places on Earth that are truly aseptic. Even environments that we may consider 'extreme', such as glaciers, deserts, or hypersaline bodies of water (Figure 1), can harbour life. The organisms that thrive in such environments - mostly microbes - are often referred to as 'extremophiles'. However, what constitutes extreme is in the eye of the beholder. Extremophilic organisms are so adapted to their environment that they perceive extreme conditions as optimal for their growth and can sometimes even be stressed by what we perceive as moderate. Stress is therefore not an optimal criterion for defining what is extreme. Instead, extreme conditions can be seen as those in which the majority of species cannot grow or even survive.

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.

The genomes of eighteen Fusarium isolates cultured from diseased and healthy citrus trees were sequenced, assembled, and annotated. Isolate species identification was confirmed using single marker (TEF1-alpha) phylogenetic assessment. Studies of the traits and genotypes of plant-associated isolates are important to understanding the fungal contribution to phytobiomes of citrus.

Filamentous fungi are an essential source of bioactive mycotoxins. Recent efforts have focused on developing antifungal agents that are effective against invasive yeasts, such as Candida spp. By screening fungal strains isolated from regions surrounding the Chernobyl nuclear power plant disaster for antifungal activity against Candida albicans, we found that Aspergillus melleus IMV 01140 produced compounds that inhibited the growth of the yeast. The active compound produced by A. melleus was isolated and found to be neoaspergillic acid, a compound that is closely related to aspergillic acid. While aspergillic acid and its derivatives have been characterized and were found to have antibacterial and antifungal properties, neoaspergillic acid has been much less studied. Even though neoaspergillic acid and related compounds were found to have antibacterial and antitumoral effects, further investigation into this group of compounds is limited by challenges associated with large-scale production, isolation, and purification. The production of neoaspergillic acid has been shown to require co-cultivation methods or special growth conditions. In this work, neoaspergillic acid and related compounds were found to be produced by A. melleus under laboratory growth conditions. The biosynthetic gene cluster of neoaspergillic acid was predicted using the aspergillic acid gene cluster as a model. The biosynthetic pathway for neoaspergillic acid was then confirmed by establishing an in vitro CRISPR-ribonucleoprotein system to individually delete genes within the cluster. A negative transcriptional factor, mcrA, was also eliminated to further improve the production of neoaspergillic acid and the related compounds for future studies.

Labyrinthulomycetes are a group of ubiquitous and diverse unicellular Stramenopiles and have long been known for their vital role in ocean carbon cycling. However, their ecological function from the perspective of organic matter degradation remains poorly understood. This study reports high-quality genomes of two newly isolated Labyrinthulomycetes strains, namely, Botryochytrium sp. strain S-28 and Oblongichytrium sp. strain S-429, and provides molecular analysis of their ecological functions using comparative genomics and a biochemical assay. Our results suggest that Labyrinthulomycetes may occupy multiple ecological niches in marine ecosystems because of the significant differences in gene function among different genera. Certain strains could degrade wheat bran independently by secreting cellulase. The key glycoside hydrolase families (GH1, GH5, and GH9) related to cellulase and the functional domains of carbohydrate-active enzymes (CAZymes) were more enriched in their genomes. This group can actively participate in marine biochemical cycles as decomposers. In contrast, other strains that could not produce cellulase may thrive as "leftover scavengers" and act as a source of nutrients to the higher-trophic-level plankton. In addition, our findings emphasize the dual roles of endoglucanase, acting as both exo- and endoglucanases, in the process of cellulose degradation. Using genomic, biochemical, and phylogenetic analyses, our study provides a broader insight into the nutritional patterns and ecological functions of Labyrinthulomycetes. IMPORTANCE Unicellular heterotrophic eukaryotes are an important component of marine ecosystems. However, their ecological functions and modes of nutrition remain largely unknown. Our current understanding of marine microbial ecology is incomplete without integrating these heterotrophic microeukaryotes into the food web models. This study focuses on the unicellular fungus-like protists Labyrinthulomycetes and provides two high-quality genomes of cellulase-producing Labyrinthulomycetes. Our study uncovers the basis of their cellulase production by deciphering the results of genomic, biochemical, and phylogenetic analyses. This study instigates a further investigation of the molecular mechanism of organic matter utilization by Labyrinthulomycetes in the world's oceans.

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.

The anaerobic gut fungi (AGF) represent a coherent phylogenetic clade within the Mycota. Twenty genera have been described so far. Currently, the phylogenetic and evolutionary relationships between AGF genera remain poorly understood. Here, we utilized 52 transcriptomic datasets from 14 genera to resolve AGF inter-genus relationships using phylogenomics, and to provide a quantitative estimate (amino acid identity, AAI) for intermediate rank assignments. We identify four distinct supra-genus clades, encompassing all genera producing polyflagellated zoospores, bulbous rhizoids, the broadly circumscribed genus Piromyces, and the Anaeromyces and affiliated genera. We also identify the genus Khoyollomyces as the earliest evolving AGF genus. Concordance between phylogenomic outputs and RPB1 and D1/D2 LSU, but not RPB2, MCM7, EF1α or ITS1, phylogenies was observed. We combine phylogenomic analysis and AAI outputs with informative phenotypic traits to propose accommodating 14/20 AGF genera into four families: Caecomycetaceae fam. nov. (encompassing the genera Caecomyces and Cyllamyces), Piromycetaceae fam. nov. (encompassing the genus Piromyces), emend the description of the family Neocallimastigaceae to encompass the genera Neocallimastix, Orpinomyces, Pecoramyces, Feramyces, Ghazallomyces, Aestipascuomyces and Paucimyces, as well as the family Anaeromycetaceae to include the genera Oontomyces, Liebetanzomyces and Capellomyces in addition to Anaeromyces. We refrain from proposing families for the deeply branching genus Khoyollomyces and for genera with uncertain position (Buwchfawromyces, Joblinomyces, Tahromyces, Agriosomyces and Aklioshbomyces) pending availability of additional isolates and sequence data; and these genera are designated as 'genera incertae sedis' in the order Neocallimastigales. Our results establish an evolutionary-grounded Linnaean taxonomic framework for the AGF, provide quantitative estimates for rank assignments, and demonstrate the utility of RPB1 as an additional informative marker in Neocallimastigomycota taxonomy.

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 Global Panzootic Lineage (GPL) of the pathogenic fungus Batrachochytrium dendrobatidis ( Bd ) has caused severe amphibian population declines, yet the drivers underlying the high frequency of GPL in regions of amphibian decline are unclear. Using publicly available Bd genome sequences, we identified multiple non-GPL Bd isolates that contain a circular Rep-encoding single stranded DNA-like virus which we named BdDV-1. We further sequenced and constructed genome assemblies with long read sequences to find that the virus is likely integrated into the nuclear genome. Attempts to cure virus positive isolates were unsuccessful, however, phenotypic differences between naturally virus positive and virus negative Bd isolates suggested that BdDV-1 decreases the growth of its host in vitro but increases the virulence of its host in vivo . BdDV-1 is the first described CRESS DNA mycovirus to increase host virulence with a distribution inversely associated with the emergence of the panzootic lineage.

Endofungal Mycetohabitans (formerly Burkholderia) spp. rely on a type III secretion system to deliver mostly unidentified effector proteins when colonizing their host fungus, Rhizopus microsporus. The one known secreted effector family from Mycetohabitans consists of homologs of transcription activator-like (TAL) effectors, which are used by plant pathogenic Xanthomonas and Ralstonia spp. to activate host genes that promote disease. These Burkholderia TAL-like (Btl) proteins bind corresponding specific DNA sequences in a predictable manner, but their genomic target(s) and impact on transcription in the fungus are unknown. Recent phenotyping of Btl mutants of two Mycetohabitans strains revealed that the single Btl in one M. endofungorum strain enhances fungal membrane stress tolerance, while others in a M. rhizoxinica strain promote bacterial colonization of the fungus. The phenotypic diversity underscores the need to assess the sequence diversity and, given that sequence diversity translates to DNA targeting specificity, the functional diversity of Btl proteins. Using a dual approach to maximize capture of Btl protein sequences for our analysis, we sequenced and assembled nine Mycetohabitans spp. genomes using long-read PacBio technology and also mined available short-read Illumina fungal-bacterial metagenomes. We show that btl genes are present across diverse Mycetohabitans strains from Mucoromycota fungal hosts yet vary in sequences and predicted DNA binding specificity. Phylogenetic analysis revealed distinct clades of Btl proteins and suggested that Mycetohabitans might contain more species than previously recognized. Within our data set, Btl proteins were more conserved across Mycetohabitans rhizoxinica strains than across Mycetohabitans endofungorum, but there was also evidence of greater overall strain diversity within the latter clade. Overall, the results suggest that Btl proteins contribute to bacterial-fungal symbioses in myriad ways.

Gene prediction is required to obtain optimal biologically meaningful information from genomic sequences, but automated gene prediction software is imperfect. In this study, we compare the original annotation of the Coccidioides immitis RS genome (the reference strain of C. immitis) to annotations using the Funannotate and Augustus genome prediction pipelines. A total of 25% of the originally predicted genes (denoted CIMG) were not found in either the Funannotate or Augustus predictions. A comparison of Funannotate and Augustus predictions also found overlapping but not identical sets of genes. The predicted genes found only in the original annotation (referred to as CIMG-unique) were less likely to have a meaningful functional annotation and a lower number of orthologs and homologs in other fungi than all CIMG genes predicted by the original annotation. The CIMG-unique genes were also more likely to be lineage-specific and poorly expressed. In addition, the CIMG-unique genes were found in clusters and tended to be more frequently associated with transposable elements than all CIMG-predicted genes. The CIMG-unique genes were more likely to have experimentally determined transcription start sites that were further away from the originally predicted transcription start sites, and experimentally determined initial transcription was less likely to result in stable CIMG-unique transcripts. A sample of CIMG-unique genes that were relatively well expressed and differentially expressed in mycelia and spherules was inspected in a genome browser, and the structure of only about half of them was found to be supported by RNA-seq data. These data suggest that some of the CIMG-unique genes are not authentic gene predictions. Genes that were predicted only by the Funannotate pipeline were also less likely to have a meaningful functional annotation, be shorter, and express less well than all the genes predicted by Funannotate. C. immitis genes predicted by more than one annotation are more likely to have predicted functions, many orthologs and homologs, and be well expressed. Lineage-specific genes are relatively uncommon in this group. These data emphasize the importance and limitations of gene prediction software and suggest that improvements to the annotation of the C. immitis genome should be considered.

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.

Despite their role in host nutrition, the anaerobic gut fungal (AGF) component of the herbivorous gut microbiome remains poorly characterized. Here, to examine global patterns and determinants of AGF diversity, we generate and analyze an amplicon dataset from 661 fecal samples from 34 mammalian species, 9 families, and 6 continents. We identify 56 novel genera, greatly expanding AGF diversity beyond current estimates (31 genera and candidate genera). Community structure analysis indicates that host phylogenetic affiliation, not domestication status and biogeography, shapes the community rather than. Fungal-host associations are stronger and more specific in hindgut fermenters than in foregut fermenters. Transcriptomics-enabled phylogenomic and molecular clock analyses of 52 strains from 14 genera indicate that most genera with preferences for hindgut hosts evolved earlier (44-58 Mya) than those with preferences for foregut hosts (22-32 Mya). Our results greatly expand the documented scope of AGF diversity and provide an ecologically and evolutionary-grounded model to explain the observed patterns of AGF diversity in extant animal hosts.

The citrus root and rhizosphere microbiomes have been relatively well described in the literature, especially in the context of Huanglonbing disease. Yet questions addressing the assembly of root microbial endophytes have remained unanswered. In the above ground tree tissues, leaves and stems have been the research focus point, while flush and flower microbiomes, two important tissues in the vegetative and reproductive cycles of the tree, are not well described. In this study, the fungal and bacterial taxa in five biocompartments (bulk soil, rhizosphere, root endosphere, flower and flush) of citrus trees grown in a single California orchard were profiled using an amplicon-based metagenomic Illumina sequencing approach. Trees with no observable signs of abiotic or biotic stresses were sampled for two consecutive years during the floral development phase. The rhizosphere was the most biodiverse compartment compared to bulk soil, root endosphere, flower and flush microbiomes. In addition, the belowground bacteriome was more diverse than the mycobiome. Microbial richness decreased significantly from the root exosphere to the endosphere and was overall low in the above ground tissues. Root endophytic microbial community composition shared strong similarities to the rhizosphere but also contained few taxa from above ground tissues. Our data indicated compartmentalization of the microbiome with distinct profiles between above and below ground microbial communities. However, several taxa were present across all compartments suggesting the existence of a core citrus microbiota. These findings highlight key microbial taxa that could be engineered as biopesticides and biofertilizers for citriculture.

BackgroundRock-dwelling microorganisms are key players in ecosystem functioning of Antarctic ice free-areas. Yet, little is known about their diversity and ecology, and further still, viruses in these communities have been largely unexplored despite important roles related to host metabolism and nutrient cycling. To begin to address this, we present a large-scale viral catalog from Antarctic rock microbial communities.ResultsWe performed metagenomic analyses on rocks from across Antarctica representing a broad range of environmental and spatial conditions, and which resulted in a predicted viral catalog comprising > 75,000 viral operational taxonomic units (vOTUS). We found largely undescribed, highly diverse and spatially structured virus communities which had predicted auxiliary metabolic genes (AMGs) with functions indicating that they may be potentially influencing bacterial adaptation and biogeochemistry.ConclusionThis catalog lays the foundation for expanding knowledge of virosphere diversity, function, spatial ecology, and dynamics in extreme environments. This work serves as a step towards exploring adaptability of microbial communities in the face of a changing climate. Video Abstract.

IntroductionSoil microbial communities, including biological soil crust microbiomes, play key roles in water, carbon and nitrogen cycling, biological weathering, and other nutrient releasing processes of desert ecosystems. However, our knowledge of microbial distribution patterns and ecological drivers is still poor, especially so for the Chihuahuan Desert.MethodsThis project investigated the effects of trampling disturbance on surface soil microbiomes, explored community composition and structure, and related patterns to abiotic and biotic landscape characteristics within the Chihuahuan Desert biome. Composite soil samples were collected in disturbed and undisturbed areas of 15 long-term ecological research plots in the Jornada Basin, New Mexico. Microbial diversity of cross-domain microbial groups (total Bacteria, Cyanobacteria, Archaea, and Fungi) was obtained via DNA amplicon metabarcode sequencing. Sequence data were related to landscape characteristics including vegetation type, landforms, ecological site and state as well as soil properties including gravel content, soil texture, pH, and electrical conductivity.ResultsFilamentous Cyanobacteria dominated the photoautotrophic community while Proteobacteria and Actinobacteria dominated among the heterotrophic bacteria. Thaumarchaeota were the most abundant Archaea and drought adapted taxa in Dothideomycetes and Agaricomycetes were most abundant fungi in the soil surface microbiomes. Apart from richness within Archaea (p = 0.0124), disturbed samples did not differ from undisturbed samples with respect to alpha diversity and community composition (p ≥ 0.05), possibly due to a lack of frequent or impactful disturbance. Vegetation type and landform showed differences in richness of Bacteria, Archaea, and Cyanobacteria but not in Fungi. Richness lacked strong relationships with soil variables. Landscape features including parent material, vegetation type, landform type, and ecological sites and states, exhibited stronger influence on relative abundances and microbial community composition than on alpha diversity, especially for Cyanobacteria and Fungi. Soil texture, moisture, pH, electrical conductivity, lichen cover, and perennial plant biomass correlated strongly with microbial community gradients detected in NMDS ordinations.DiscussionOur study provides first comprehensive insights into the relationships between landscape characteristics, associated soil properties, and cross-domain soil microbiomes in the Chihuahuan Desert. Our findings will inform land management and restoration efforts and aid in the understanding of processes such as desertification and state transitioning, which represent urgent ecological and economical challenges in drylands around the world.

National Aeronautics and Space Administration’s (NASA) spacecraft assembly facilities are monitored for the presence of any bacteria or fungi that might conceivably survive a transfer to an extraterrestrial environment. Fungi present a broad and diverse range of phenotypic and functional traits to adapt to extreme conditions, hence the detection of fungi and subsequent eradication of them are needed to prevent forward contamination for future NASA missions. During the construction and assembly for the Mars 2020 mission, three fungal strains with unique morphological and phylogenetic properties were isolated from spacecraft assembly facilities. The reconstruction of phylogenetic trees based on several gene loci (ITS, LSU, SSU, RPB, TUB, TEF1) using multi-locus sequence typing (MLST) and whole genome sequencing (WGS) analyses supported the hypothesis that these were novel species. Here we report the genus or species-level classification of these three novel strains via a polyphasic approach using phylogenetic analysis, colony and cell morphology, and comparative analysis of WGS. The strain FJI-L9-BK-P1 isolated from the Jet Propulsion Laboratory Spacecraft Assembly Facility (JPL-SAF) exhibited a putative phylogenetic relationship with the strain Aaosphaeria arxii CBS175.79 but showed distinct morphology and microscopic features. Another JPL-SAF strain, FJII-L3-CM-DR1, was phylogenetically distinct from members of the family Trichomeriaceae and exhibited morphologically different features from the genera Lithohypha and Strelitziana. The strain FKI-L1-BK-DR1 isolated from the Kennedy Space Center facility was identified as a member of Dothideomycetes incertae sedis and is closely related to the family Kirschsteiniotheliaceae according to a phylogenetic analysis. The polyphasic taxonomic approach supported the recommendation for establishing two novel genera and one novel species. The names Aaosphaeria pasadenensis (FJI-L9-BK-P1 = NRRL 64424 = DSM 114621), Pasadenomyces melaninifex (FJII-L3-CM-DR1 = NRRL 64433 = DSM 114623), and Floridaphiala radiotolerans (FKI-L1-BK-DR1 = NRRL 64434 = DSM 114624) are proposed as type species. Furthermore, resistance to ultraviolet-C and presence of specific biosynthetic gene cluster(s) coding for metabolically active compounds are unique to these strains.

Kākāpō are a critically endangered species of parrots restricted to a few islands off the coast of New Zealand. Kākāpō are very closely monitored, especially during nesting seasons. In 2019, during a highly successful nesting season, an outbreak of aspergillosis affected 21 individuals and led to the deaths of 9, leaving a population of only 211 kākāpō. In monitoring this outbreak, cultures of aspergillus were grown, and genome sequenced. These sequences demonstrate that, very unusually for an aspergillus outbreak, a single strain of aspergillus caused the outbreak. This strain was found on two islands, but only one had an outbreak of aspergillosis; indicating that the strain was necessary, but not sufficient, to cause disease. Our analysis provides an understanding of the 2019 outbreak and provides potential ways to manage such events in the future.

The fungal kingdom represents an extraordinary diversity of organisms with profound impacts across animal, plant, and ecosystem health. Fungi simultaneously support life, by forming beneficial symbioses with plants and producing life-saving medicines, and bring death, by causing devastating diseases in humans, plants, and animals. With climate change, increased antimicrobial resistance, global trade, environmental degradation, and novel viruses altering the impact of fungi on health and disease, developing new approaches is now more crucial than ever to combat the threats posed by fungi and to harness their extraordinary potential for applications in human health, food supply, and environmental remediation. To address this aim, the Canadian Institute for Advanced Research (CIFAR) and the Burroughs Wellcome Fund convened a workshop to unite leading experts on fungal biology from academia and industry to strategize innovative solutions to global challenges and fungal threats. This report provides recommendations to accelerate fungal research and highlights the major research advances and ideas discussed at the meeting pertaining to 5 major topics: (1) Connections between fungi and climate change and ways to avert climate catastrophe; (2) Fungal threats to humans and ways to mitigate them; (3) Fungal threats to agriculture and food security and approaches to ensure a robust global food supply; (4) Fungal threats to animals and approaches to avoid species collapse and extinction; and (5) Opportunities presented by the fungal kingdom, including novel medicines and enzymes.

Feline-transmitted sporotrichosis has garnered attention due to the recent high incidence and the lack of efficient control in the epicenter of the epidemic, Rio de Janeiro, Brazil. Sporothrix brasiliensis is the major pathogen involved in feline-to-human sporotrichosis in Brazil and displays more virulent genotypes than the closely related species S. schenckii. Over the last two decades, several reports of antifungal-resistant strains have emerged. Sequencing and comparison analysis of the outbreak strains allowed us to observe that the azole non-wild-type S. brasiliensis strain CFP 1054 had significant chromosomal variations compared to wild-type strains. One of these variants includes a region of 231 Kb containing 75 duplicated genes, which were overrepresented for lipid and isoprenoid metabolism. We also identified an additional strain (CFP 1055) that was resistant to itraconazole and amphotericin B, which had a single nucleotide polymorphism in the tac1 gene. The patients infected with these two strains showed protracted clinical course and sequelae. Even though our sample size is modest, these results suggest the possibility of identifying specific point mutations and large chromosomal duplications potentially associated with antifungal resistance and clinical outcomes of sporotrichosis.

Genome sequencing of filamentous fungi has demonstrated that most secondary metabolite biosynthetic gene clusters (BGCs) are silent under standard laboratory conditions. In this work, we have established an in vitro CRISPR-Cas9 system in Aspergillus wentii. To activate otherwise silent BGCs, we deleted the negative transcriptional regulator mcrA. Deletion of mcrA (mcrAΔ) resulted in differential production of 17 SMs in total when the strain was cultivated on potato dextrose media (PDA). Nine out of fifteen of these SMs were fully characterized, including emodin (1), physcion (2), sulochrin (3), physcion bianthrone (4), 14-O-demethylsulochrin (5), (trans/cis)-emodin bianthrone (6 and 7), and (trans/cis)-emodin physcion bianthrone (8 and 9). These compounds were all found to be produced by the same polyketide synthase (PKS) BGC. We then performed a secondary knockout targeting this PKS cluster in the mcrAΔ background. The metabolite profile of the dual-knockout strain revealed new metabolites that were not previously detected in the mcrAΔ parent strain. Two additional SMs were purified from the dual-knockout strain and were characterized as aspergillus acid B (16) and a structurally related but previously unidentified compound (17). For the first time, this work presents a facile genetic system capable of targeted gene editing in A. wentii. This work also illustrates the utility of performing a dual knockout to eliminate major metabolic products, enabling additional SM discovery.

We report six metagenome-assembled genomes (MAGs) associated with Massospora cicadina strain MCPNR19 (ARSEF 14555), an obligate entomopathogenic fungus of periodical cicadas. The MAGs include representatives of Pantoea, Pseudomonas, Lactococcus, and one potential new Chryseobacterium species. Future research is needed to resolve the ecology of these MAGs and determine whether they represent symbionts or contaminants.

A 1.488-Gb draft genome sequence was assembled for the fungus Massospora cicadina, an obligate parasite of periodical cicadas. The M. cicadina genome has experienced massive expansion via transposable elements (TEs), which account for 92% of the genome.

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.

Cancer-microbe associations have been explored for centuries, but cancer-associated fungi have rarely been examined. Here, we comprehensively characterize the cancer mycobiome within 17,401 patient tissue, blood, and plasma samples across 35 cancer types in four independent cohorts. We report fungal DNA and cells at low abundances across many major human cancers, with differences in community compositions that differ among cancer types, even when accounting for technical background. Fungal histological staining of tissue microarrays supported intratumoral presence and frequent spatial association with cancer cells and macrophages. Comparing intratumoral fungal communities with matched bacteriomes and immunomes revealed co-occurring bi-domain ecologies, often with permissive, rather than competitive, microenvironments and distinct immune responses. Clinically focused assessments suggested prognostic and diagnostic capacities of the tissue and plasma mycobiomes, even in stage I cancers, and synergistic predictive performance with bacteriomes.

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.

The role of microbes in insecticide resistance is an emerging question. Here, we describe six metagenome-assembled genomes (MAGs) associated with the glassy-winged sharpshooter (Homalodisca vitripennis [Germar, 1821]) (Hemiptera, Cicadellidae). MAGs representing the obligate symbionts Candidatus Sulcia muelleri and Candidatus Baumannia cicadellinicola and the facultative symbiont Wolbachia were obtained from imidacloprid-resistant and imidacloprid-susceptible sharpshooters.

Abstract While the study of the tumor microbiome and its effects on cancer biology has expanded considerably over the last few years, most of this research focused on bacteria and viruses, leaving behind the fungal kingdom. Recently, a few studies have demonstrated that specific fungi may promote tumor progression, stressing the importance of comprehensively studying the tumor mycobiome and its effects. To address this, we have characterized the mycobiome in 1,183 human tumors and their adjacent tissues, originating from eight major solid tumor types. Staining and imaging demonstrated the presence of fungi in both cancer and immune cells, with tumor-type specific distribution patterns. Quantitative PCR of the fungal 5.8s rDNA revealed the presence of fungal DNA in all tumor types. To characterize the tumor mycobiome and address potential contamination during tissue handling and processing, we subjected all samples, as well as 295 negative controls of different types, to sequencing of the ITS2 region that is situated between fungal rRNA genes. We found cancer-type specific mycobial signatures with relatively high similarity between tumors and their adjacent tissues. While the fungal mycobiome had a lower species richness as compared to the bacterial microbiome of the same tumors, fungi showed significant co-occurrences with specific bacteria, suggesting the existence of ecological niches within the tumors. We also found significant correlations with clinical parameters such as patient’s age, tumor stage, progression-free survival, overall survival, and response to immune checkpoint blockade therapy. Characterization of the tumor mycobiome may add a biologically relevant, previously overlooked, component to be considered in the study of cancer, including its effects on tumor initiation, progression, diagnosis, and response to therapy. Citation Format: Lian Narunsky Haziza, Gregory D. Sepich-Poore, Ilana Livyatan, Omer Asraf, Cameron Martino, Deborah Nejman, Nancy Gavert, Jason E. Stajich, Guy Amit, Antonio González, Stephen Wandro, Gili Perry, Ruthie Ariel, Arnon Meltser, Justin P. Shaffer, Qiyun Zhu, Nora Balint-Lahat, Iris Barshack, Maya Dadian, Einav N. Gal-Yam, Sandip P. Pate, Amir Bashan, Austin D. Swafford, Yitzhak Pilpel, Rob Knight, Ravid Straussman. Pan-cancer characterization of the tumor mycobiome and its clinical effects [abstract]. In: Proceedings of the American Association for Cancer Research Annual Meeting 2022; 2022 Apr 8-13. Philadelphia (PA): AACR; Cancer Res 2022;82(12_Suppl):Abstract nr 3054.