Genome-Wide Survey of Gut Fungi (Harpellales) Reveals the First Horizontally Transferred Ubiquitin Gene from a Mosquito Host

Arctic cyanobacterial mat community diversity decreases with latitude across the Canadian Arctic

Cyanobacterial mats are commonly reported as hotspots of microbial diversity across polar environments. These thick, multilayered microbial communities provide a refuge from extreme environmental conditions, with many species able to grow and coexist despite the low allochthonous nutrient inputs. The visibly dominant phototrophic biomass is dependent on internal nutrient recycling by heterotrophic organisms within the mats; however, the specific contribution of heterotrophic protists remains little explored. In this study, mat community diversity was examined along a latitudinal gradient (55-83°N), spanning subarctic taiga, tundra, polar desert, and the High Arctic ice shelves. The prokaryotic and eukaryotic communities were targeted, respectively, by V4 16S ribosomal RNA (rRNA) and V9 18S rRNA gene amplicon high-throughput sequencing. Prokaryotic and eukaryotic richness decreased, in tandem with decreasing temperatures and shorter seasons of light availability, from the subarctic to the High Arctic. Taxonomy-based annotation of the protist community revealed diverse phototrophic, mixotrophic, and heterotrophic genera in all mat communities, with fewer parasitic taxa in High Arctic communities. Co-occurrence network analysis identified greater heterogeneity in eukaryotic than prokaryotic community structure among cyanobacterial mats across the Canadian Arctic. Our findings highlight the sensitivity of microbial eukaryotes to environmental gradients across northern high latitudes.

Conceptual risk assessment of mosquito population modification gene-drive systems to control malaria transmission: preliminary hazards list workshops

The field-testing and eventual adoption of genetically-engineered mosquitoes (GEMs) to control vector-borne pathogen transmission will require them meeting safety criteria specified by regulatory authorities in regions where the technology is being considered for use and other locales that might be impacted. Preliminary risk considerations by researchers and developers may be useful for planning the baseline data collection and field research used to address the anticipated safety concerns. Part of this process is to identify potential hazards (defined as the inherent ability of an entity to cause harm) and their harms, and then chart the pathways to harm and evaluate their probability as part of a risk assessment. The University of California Malaria Initiative (UCMI) participated in a series of workshops held to identify potential hazards specific to mosquito population modification strains carrying gene-drive systems coupled to anti-parasite effector genes and their use in a hypothetical island field trial. The hazards identified were placed within the broader context of previous efforts discussed in the scientific literature. Five risk areas were considered i) pathogens, infections and diseases, and the impacts of GEMs on human and animal health, ii) invasiveness and persistence of GEMs, and interactions of GEMs with target organisms, iii) interactions of GEMs with non-target organisms including horizontal gene transfer, iv) impacts of techniques used for the management of GEMs and v) evolutionary and stability considerations. A preliminary hazards list (PHL) was developed and is made available here. This PHL is useful for internal project risk evaluation and is available to regulators at prospective field sites. UCMI project scientists affirm that the subsequent processes associated with the comprehensive risk assessment for the application of this technology should be driven by the stakeholders at the proposed field site and areas that could be affected by this intervention strategy.

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.

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.

Comparative Genomics of Mortierellaceae Provides Insights into Lipid Metabolism: Two Novel Types of Fatty Acid Synthase

Fungal species in the family Mortierellaceae are important for their remarkable capability to synthesize large amounts of polyunsaturated fatty acids, especially arachidonic acid (ARA). Although many genomes have been published, the quality of these data is not satisfactory, resulting in an incomplete understanding of the lipid pathway in Mortierellaceae. We provide herein two novel and high-quality genomes with 55.32% of syntenic gene pairs for Mortierella alpina CGMCC 20262 and M. schmuckeri CGMCC 20261, spanning 28 scaffolds of 40.22 Mb and 25 scaffolds of 49.24 Mb, respectively. The relative smaller genome for the former is due to fewer protein-coding gene models (11,761 vs. 13,051). The former yields 45.57% of ARA in total fatty acids, while the latter 6.95%. The accumulation of ARA is speculated to be associated with delta-5 desaturase (Delta5) and elongation of very long chain fatty acids protein 3 (ELOVL3). A further genomic comparison of 19 strains in 10 species in three genera in the Mortierellaceae reveals three types of fatty acid synthase (FAS), two of which are new to science. The most common type I exists in 16 strains of eight species of three genera, and was discovered previously and consists of a single unit with eight active sites. The newly revealed type II exists only in M. antarctica KOD 1030 where the unit is separated into two subunits α and β comprised of three and five active sites, respectively. Another newly revealed type III exists in M. alpina AD071 and Dissophora globulifera REB-010B, similar to type II but different in having one more acyl carrier protein domain in the α subunit. This study provides novel insights into the enzymes related to the lipid metabolism, especially the ARA-related Delta5, ELOVL3, and FAS, laying a foundation for genetic engineering of Mortierellaceae to modulate yield in polyunsaturated fatty acids.

Distribution of methionine sulfoxide reductases in fungi and conservation of the free-methionine-R-sulfoxide reductase in multicellular eukaryotes

Methionine, either as a free amino acid or included in proteins, can be oxidized into methionine sulfoxide (MetO), which exists as R and S diastereomers. Almost all characterized organisms possess thiol-oxidoreductases named methionine sulfoxide reductase (Msr) enzymes to reduce MetO back to Met. MsrA and MsrB reduce the S and R diastereomers of MetO, respectively, with strict stereospecificity and are found in almost all organisms. Another type of thiol-oxidoreductase, the free-methionine-R-sulfoxide reductase (fRMsr), identified so far in prokaryotes and a few unicellular eukaryotes, reduces the R MetO diastereomer of the free amino acid. Moreover, some bacteria possess molybdenum-containing enzymes that reduce MetO, either in the free or protein-bound forms. All these Msrs play important roles in the protection of organisms against oxidative stress. Fungi are heterotrophic eukaryotes that colonize all niches on Earth and play fundamental functions, in organic matter recycling, as symbionts, or as pathogens of numerous organisms. However, our knowledge on fungal Msrs is still limited. Here, we performed a survey of msr genes in almost 700 genomes across the fungal kingdom. We show that most fungi possess one gene coding for each type of methionine sulfoxide reductase: MsrA, MsrB, and fRMsr. However, several fungi living in anaerobic environments or as obligate intracellular parasites were devoid of msr genes. Sequence inspection and phylogenetic analyses allowed us to identify non-canonical sequences with potentially novel enzymatic properties. Finaly, we identified several ocurences of msr horizontal gene transfer from bacteria to fungi.

The genome of Geosiphon pyriformis reveals ancestral traits linked to the emergence of the arbuscular mycorrhizal symbiosis

Arbuscular mycorrhizal fungi (AMF) (subphylum Glomeromycotina)¹ are among the most prominent symbionts and form the Arbuscular Mycorrhizal symbiosis (AMS) with over 70% of known land plants.^2,3 AMS allows plants to efficiently acquire poorly soluble soil nutrients⁴ and AMF to receive photosynthetically fixed carbohydrates. This plant-fungus symbiosis dates back more than 400 million years⁵ and is thought to be one of the key innovations that allowed the colonization of lands by plants.⁶ Genomic and genetic analyses of diverse plant species started to reveal the molecular mechanisms that allowed the evolution of this symbiosis on the host side, but how and when AMS abilities emerged in AMF remain elusive. Comparative phylogenomics could be used to understand the evolution of AMS.^7,8 However, the availability of genome data covering basal AMF phylogenetic nodes (Archaeosporales, Paraglomerales) is presently based on fragmentary protein coding datasets.⁹Geosiphon pyriformis (Archaeosporales) is the only fungus known to produce endosymbiosis with nitrogen-fixing cyanobacteria (Nostoc punctiforme) presumably representing the ancestral AMF state.^10-12 Unlike other AMF, it forms long fungal cells ("bladders") that enclose cyanobacteria. Once in the bladder, the cyanobacteria are photosynthetically active and fix nitrogen, receiving inorganic nutrients and water from the fungus. Arguably, G. pyriformis represents an ideal candidate to investigate the origin of AMS and the emergence of a unique endosymbiosis. Here, we aimed to advance knowledge in these questions by sequencing the genome of G. pyriformis, using a re-discovered isolate.

Phylogenomic Analyses of Non-Dikarya Fungi Supports Horizontal Gene Transfer Driving Diversification of Secondary Metabolism in the Amphibian Gastrointestinal Symbiont, Basidiobolus

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

Influences of a Prolific Gut Fungus (Zancudomyces culisetae) on Larval and Adult Mosquito (Aedes aegypti)-Associated Microbiota

Adult mosquitoes inherit a bacterial community from larvae via transstadial transmission, an understudied process that may influence host-microbe interactions. Microbes contribute to important host life history traits, and analyzing transmitted microbial communities, the interrelationship between larval and adult-associated microbiota, and factors influencing host-microbe relationships provides targets for research. During its larval stage, the yellow fever mosquito (Aedes aegypti) hosts the trichomycete gut fungus Zancudomyces culisetae, and fungal colonization coincides with environmental perturbations in the digestive tract microecosystem. Natural populations are differentially exposed to fungi, thereby potentially harboring distinct microbiota and experiencing disparate host-microbe interactions. This study's objectives were to characterize larval and initial adult microbiomes, investigate variation in diversity and distribution of microbial communities across individuals, and assess whether larval fungal colonization impacted microbiomes at these developmental stages. Laboratory-based fungal infestation assays, sequencing of 16S rRNA gene amplicons, and bacterial load quantification protocols revealed that initial adult microbiomes varied in diversity and distribution. Larval fungal colonization had downstream effects on initial adult microbiomes, significantly reducing microbial community variation, shifting relative abundances of certain bacterial families, and influencing transstadial transmission outcomes of particular genera. Further, abundances of several families consistently decreased in adults relative to levels in larvae, possibly reflecting impacts of host development on specific bacterial taxa. These findings demonstrated that a prolific gut fungus impacted mosquito-associated microbiota at two developmental stages in an insect connected with global human health.IMPORTANCE Mosquitoes are widespread vectors of numerous human pathogens and harbor microbiota known to affect host phenotypic traits. However, little research has directly investigated how bacterial communities associated with larvae and adults are connected. We characterized whole-body bacterial communities in mosquito larvae preceding pupation and in newly emerged adults, and investigated whether a significant biotic factor, fungal colonization of the larval hindgut, impacted these microbiomes. Results showed that fungal colonization reduced microbial community variation across individuals and differentially impacted the outcomes of transstadial transmission for certain bacterial genera, revealing downstream effects of the fungus on initial adult microbiomes. The importance of our research is in providing a thorough comparative analysis of whole-body microbiota harbored in larvae and adults of the yellow fever mosquito (Aedes aegypti) and in demonstrating the important role a widespread gut fungus played in a host-associated microbiome.

Advances and Applications in the Quest for Orthologs

Gene families evolve by the processes of speciation (creating orthologs), gene duplication (paralogs), and horizontal gene transfer (xenologs), in addition to sequence divergence and gene loss. Orthologs in particular play an essential role in comparative genomics and phylogenomic analyses. With the continued sequencing of organisms across the tree of life, the data are available to reconstruct the unique evolutionary histories of tens of thousands of gene families. Accurate reconstruction of these histories, however, is a challenging computational problem, and the focus of the Quest for Orthologs Consortium. We review the recent advances and outstanding challenges in this field, as revealed at a symposium and meeting held at the University of Southern California in 2017. Key advances have been made both at the level of orthology algorithm development and with respect to coordination across the community of algorithm developers and orthology end-users. Applications spanned a broad range, including gene function prediction, phylostratigraphy, genome evolution, and phylogenomics. The meetings highlighted the increasing use of meta-analyses integrating results from multiple different algorithms, and discussed ongoing challenges in orthology inference as well as the next steps toward improvement and integration of orthology resources.

Molecular Dating of the Emergence of Anaerobic Rumen Fungi and the Impact of Laterally Acquired Genes

The anaerobic gut fungi (AGF), or Neocallimastigomycota, inhabit the rumen and alimentary tract of herbivorous mammals, where they play important roles in the degradation of plant fiber. Comparative genomic and phylogenomic analyses of the AGF have long been hampered by their fastidious growth condition, as well as their large (up to 200 Mb) and AT-biased (78 to 84%) genomes. We sequenced 21 AGF transcriptomes and combined them with 5 available AGF genome sequences to explore their evolutionary relationships, time their divergence, and characterize gene gain/loss patterns associated with their evolution. We estimate that the most recent common ancestor of the AGF diverged 66 (±10) million years ago, a time frame that coincides with the evolution of grasses (Poaceae), as well as the mammalian transition from insectivory to herbivory. The concordance of independent estimations suggests that AGF have been important in shaping the success of mammalian herbivory transition by improving the efficiency of energy acquisition from recalcitrant plant materials. Comparative genomics identified multiple lineage-specific genes in the AGF, two of which were acquired from rumen gut bacteria and animal hosts via horizontal gene transfer (HGT). A third AGF domain, plant-like polysaccharide lyase, represents a novel gene in fungi that potentially aids AGF to degrade pectin. Analysis of genomic and transcriptomic sequences confirmed both the presence and expression of these lineage-specific genes in nearly all AGF clades. These genetic elements may contribute to the exceptional abilities of AGF to degrade plant biomass and enable metabolism of the rumen microbes and animal hosts.IMPORTANCE Anaerobic fungi living in the rumen of herbivorous mammals possess an extraordinary ability to degrade plant biomass. We examined the origin and genomic composition of these poorly characterized anaerobic gut fungi using both transcriptome and genomic data. Phylogenomics and molecular dating analyses found remarkable concurrence of the divergence times of the rumen fungi, the forage grasses, and the dietary shift of ancestral mammals from primarily insectivory to herbivory. Comparative genomics identified unique machinery in these fungi to utilize plant polysaccharides. The rumen fungi were also identified with the ability to code for three protein domains with putative functions in plant pectin degradation and microbial defense, which were absent from all other fungal organisms (examined over 1,000 fungal genomes). Two of these domains were likely acquired from rumen gut bacteria and animal hosts separately via horizontal gene transfer. The third one is a plant-like polysaccharide lyase, representing a unique fungal enzyme with potential pectin breakdown abilities.

Antibiotic resistant bacteria and commensal fungi are common and conserved in the mosquito microbiome

The emerging and increasing prevalence of bacterial antibiotic resistance is a significant public health challenge. To begin to tackle this problem, it will be critical to not only understand the origins of this resistance but also document environmental reservoirs of antibiotic resistance. In this study we investigated the possibility that both colony and field caught mosquitoes could harbor antibiotic resistant bacteria. Specifically, we characterized the antibiotic resistant bacterial populations from colony-reared Aedes aegypti larvae and adults and two field caught mosquito species Coquillettidia perturbans and Ochlerotatus canadensis. The cultured bacterial populations were dominated by isolates belonging to the class Gammaproteobacteria. Among the antibiotic resistant populations, we found bacteria resistant to carbenicillin, kanamycin, and tetracycline, including bacteria resistant to a cocktail of all three antibiotics in combination. The antibiotic resistant bacteria were numerically rare, at most 5% of total cell counts. Isolates were characterized by 16S rRNA gene sequencing, and clustering into Operational Taxonomic Units (OTUs; 99% sequence identity). 27 antibiotic resistant OTUs were identified, although members of an OTU did not always share the same resistance profile. This suggests the clustering was either not sensitive enough to distinguish different bacteria taxa or different antibiotic resistant sub-populations exist within an OTU. Finally, the antibiotic selection opened up a niche to culture mosquito-associated fungi, and 10 fungal OTUs (28S rRNA gene sequencing) were identified. Two fungal OTUs both classified to the class Microbotryomycetes were commonly identified in the field-caught mosquitoes. Thus, in this study we demonstrate that antibiotic resistant bacteria and certain fungi are common and conserved mosquito microbiome members. These observations highlight the potential of invertebrates to serve as vehicles for the spread of antibiotic resistance throughout the environment.

Comparative Genomics Reveals the Core Gene Toolbox for the Fungus-Insect Symbiosis

Modern genomics has shed light on many entomopathogenic fungi and expanded our knowledge widely; however, little is known about the genomic features of the insect-commensal fungi. Harpellales are obligate commensals living in the digestive tracts of disease-bearing insects (black flies, midges, and mosquitoes). In this study, we produced and annotated whole-genome sequences of nine Harpellales taxa and conducted the first comparative analyses to infer the genomic diversity within the members of the Harpellales. The genomes of the insect gut fungi feature low (26% to 37%) GC content and large genome size variations (25 to 102 Mb). Further comparisons with insect-pathogenic fungi (from both Ascomycota and Zoopagomycota), as well as with free-living relatives (as negative controls), helped to identify a gene toolbox that is essential to the fungus-insect symbiosis. The results not only narrow the genomic scope of fungus-insect interactions from several thousands to eight core players but also distinguish host invasion strategies employed by insect pathogens and commensals. The genomic content suggests that insect commensal fungi rely mostly on adhesion protein anchors that target digestive system, while entomopathogenic fungi have higher numbers of transmembrane helices, signal peptides, and pathogen-host interaction (PHI) genes across the whole genome and enrich genes as well as functional domains to inactivate the host inflammation system and suppress the host defense. Phylogenomic analyses have revealed that genome sizes of Harpellales fungi vary among lineages with an integer-multiple pattern, which implies that ancient genome duplications may have occurred within the gut of insects.

Insect guts harbor various microbes that are important for host digestion, immune response, and disease dispersal in certain cases. Bacteria, which are among the primary endosymbionts, have been studied extensively. However, fungi, which are also frequently encountered, are poorly known with respect to their biology within the insect guts. To understand the genomic features and related biology, we produced the whole-genome sequences of nine gut commensal fungi from disease-bearing insects (black flies, midges, and mosquitoes). The results show that insect gut fungi tend to have low GC content across their genomes. By comparing these commensals with entomopathogenic and free-living fungi that have available genome sequences, we found a universal core gene toolbox that is unique and thus potentially important for the insect-fungus symbiosis. This comparative work also uncovered different host invasion strategies employed by insect pathogens and commensals, as well as a model system to study ancient fungal genome duplication within the gut of insects.

Sources of Fungal Genetic Variation and Associating It with Phenotypic Diversity

The first eukaryotic genome to be sequenced was fungal, and there continue to be more sequenced genomes in the kingdom Fungi than in any other eukaryotic kingdom. Comparison of these genomes reveals many sources of genetic variation, from single nucleotide polymorphisms to horizontal gene transfer and on to changes in the arrangement and number of chromosomes, not to mention endofungal bacteria and viruses. Population genomics shows that all sources generate variation all the time and implicate natural selection as the force maintaining genome stability. Variation in wild populations is a rich resource for associating genetic variation with phenotypic variation, whether through quantitative trait locus mapping, genome-wide association studies, or reverse ecology. Subjects of studies associating genetic and phenotypic variation include model fungi, e.g., Saccharomyces and Neurospora, but pioneering studies have also been made with fungi pathogenic to plants, e.g., Pyricularia (= Magnaporthe), Zymoseptoria, and Fusarium, and to humans, e.g., Coccidioides, Cryptococcus, and Candida.

Draft Genome Sequence of Capniomyces stellatus, the Obligate Gut Fungal Symbiont of Stonefly

Capniomyces stellatus is a host-specific endosymbiotic fungus, living in the hindgut of stoneflies (especially in Allocapnia). Here, we present the first draft genome sequence of the fungus, as well as the ab initio gene prediction and function analyses, which will facilitate the study and comparative analyses of insect-associated fungi.