kb_DRAM: annotation and metabolic profiling of genomes with DRAM in KBase.

Microbial genome annotation is the process of identifying structural and functional elements in DNA sequences and subsequently attaching biological information to those elements. DRAM is a tool developed to annotate bacterial, archaeal, and viral genomes derived from pure cultures or metagenomes. DRAM goes beyond traditional annotation tools by distilling multiple gene annotations to genome level summaries of functional potential. Despite these benefits, a downside of DRAM is the requirement of large computational resources, which limits its accessibility. Further, it did not integrate with downstream metabolic modeling tools that require genome annotation. To alleviate these constraints, DRAM and the viral counterpart, DRAM-v, are now available and integrated with the freely accessible KBase cyberinfrastructure. With kb_DRAM users can generate DRAM annotations and functional summaries from microbial or viral genomes in a point-and-click interface, as well as generate genome-scale metabolic models from DRAM annotations. AVAILABILITY AND IMPLEMENTATION: For kb_DRAM users, the kb_DRAM apps on KBase can be found in the catalog at https://narrative.kbase.us/#catalog/modules/kb_DRAM. For kb_DRAM users, a tutorial workflow with all documentation is available at https://narrative.kbase.us/narrative/129480. For kb_DRAM developers, software is available at https://github.com/shafferm/kb_DRAM.

High-resolution sequencing reveals unexplored archaeal diversity in freshwater wetland soils.

Despite being key contributors to biogeochemical processes, archaea are frequently outnumbered by bacteria, and consequently are underrepresented in combined molecular surveys. Here, we demonstrate an approach to concurrently survey the archaea alongside the bacteria with high-resolution 16S rRNA gene sequencing, linking these community data to geochemical parameters. We applied this integrated analysis to hydric soils sampled across a model methane-emitting freshwater wetland. Geochemical profiles, archaeal communities, and bacterial communities were independently correlated with soil depth and water cover. Centimeters of soil depth and corresponding geochemical shifts consistently affected microbial community structure more than hundreds of meters of lateral distance. Methanogens with diverse metabolisms were detected across the wetland, but displayed surprising OTU-level partitioning by depth. Candidatus Methanoperedens spp. archaea thought to perform anaerobic oxidation of methane linked to iron reduction were abundant. Domain-specific sequencing also revealed unexpectedly diverse non-methane-cycling archaeal members. OTUs within the underexplored Woesearchaeota and Bathyarchaeota were prevalent across the wetland, with subgroups and individual OTUs exhibiting distinct occupancy and abundance distributions aligned with environmental gradients. This study adds to our understanding of ecological range for key archaeal taxa in a model freshwater wetland, and links these taxa and individual OTUs to hypotheses about processes governing biogeochemical cycling.

Disturbed subsurface microbial communities follow equivalent trajectories despite different structural starting points.

Microbial community structure, and niche and neutral processes can all influence response to disturbance. Here, we provide experimental evidence for niche versus neutral and founding community effects during a bioremediation-related organic carbon disturbance. Subsurface sediment, partitioned into 22 flow-through columns, was stimulated in situ by the addition of acetate as a carbon and electron donor source. This drove the system into a new transient biogeochemical state characterized by iron reduction and enriched Desulfuromonadales, Comamonadaceae and Bacteroidetes lineages. After approximately 1 month conditions favoured sulfate reduction, and were accompanied by a substantial increase in the relative abundance of Desulfobulbus, Desulfosporosinus, Desulfitobacterium and Desulfotomaculum. Two subsets of four to five columns each were switched from acetate to lactate amendment during either iron (earlier) or sulfate (later) reduction. Hence, subsets had significantly different founding communities. All lactate treatments exhibited lower relative abundances of Desulfotomaculum and Bacteroidetes, enrichments of Clostridiales and Psychrosinus species, and a temporal succession from highly abundant Clostridium sensu stricto to Psychrosinus. Regardless of starting point, lactate-switch communities followed comparable structural trajectories, whereby convergence was evident 9 to 16 days after each switch, and significant after 29 to 34 days of lactate addition. Results imply that neither the founding community nor neutral processes influenced succession following perturbation.

A conserved electrostatic membrane-binding surface in synaptotagmin-like proteins revealed using molecular phylogenetic analysis and homology modeling.

Protein structure prediction has emerged as a core technology for understanding biomolecules and their interactions. Here, we combine homology-based structure prediction with molecular phylogenetic analysis to study the evolution of electrostatic membrane binding among the vertebrate synaptotagmin-like protein (Slp) family. Slp family proteins play key roles in the membrane trafficking of large dense-core secretory vesicles. Our previous experimental and computational study found that the C2A domain of Slp-4 (also called granuphilin) binds with high affinity to anionic phospholipids in the cytoplasmic leaflet of the plasma membrane through a large positively charged protein surface centered on a cluster of phosphoinositide-binding lysine residues. Because this surface contributes greatly to Slp-4 C2A domain membrane binding, we hypothesized that the net charge on the surface might be evolutionarily conserved. To test this hypothesis, the known C2A sequences of Slp-4 among vertebrates were organized by class (from mammalia to pisces) using molecular phylogenetic analysis. Consensus sequences for each class were then identified and used to generate homology structures, from which Poisson-Boltzmann electrostatic potentials were calculated. For comparison, homology structures and electrostatic potentials were also calculated for the five human Slp protein family members. The results demonstrate that the charge on the membrane-binding surface is highly conserved throughout the evolution of Slp-4, and more highly conserved than many individual residues among the human Slp family paralogs. Such molecular phylogenetic-driven computational analysis can help to describe the evolution of electrostatic interactions between proteins and membranes which are crucial for their function.

A novel three-unit tRNA splicing endonuclease found in ultrasmall Archaea possesses broad substrate specificity.

tRNA splicing endonucleases, essential enzymes found in Archaea and Eukaryotes, are involved in the processing of pre-tRNA molecules. In Archaea, three types of splicing endonuclease [homotetrameric: α(4), homodimeric: α(2), and heterotetrameric: (αß)(2)] have been identified, each representing different substrate specificity during the tRNA intron cleavage. Here, we discovered a fourth type of archaeal tRNA splicing endonuclease (ε(2)) in the genome of the acidophilic archaeon Candidatus Micrarchaeum acidiphilum, referred to as ARMAN-2 and its closely related species, ARMAN-1. The enzyme consists of two duplicated catalytic units and one structural unit encoded on a single gene, representing a novel three-unit architecture. Homodimeric formation was confirmed by cross-linking assay, and site-directed mutagenesis determined that the conserved L10-pocket interaction between catalytic and structural unit is necessary for the assembly. A tRNA splicing assay reveal that ε(2) endonuclease cleaves both canonical and non-canonical bulge-helix-bulge motifs, similar to that of (αß)(2) endonuclease. Unlike other ARMAN and Euryarchaeota, tRNAs found in ARMAN-2 are highly disrupted by introns at various positions, which again resemble the properties of archaeal species with (αß)(2) endonuclease. Thus, the discovery of ε(2) endonuclease in an archaeon deeply branched within Euryarchaeota represents a new example of the coevolution of tRNA and their processing enzymes.

A Conserved Electrostatic Membrane-Binding Surface in Synaptotagmin-Like Proteins Revealed Using Molecular Phylogenetic Analysis and Homology Modeling.

Protein structure prediction has emerged as a core technology for understanding biomolecules and their interactions. Here, we combine homology-based structure prediction with molecular phylogenetic analysis to study the evolution of electrostatic membrane binding among vertebrate synaptotagmin-like proteins (Slps). Slp family proteins play key roles in the membrane trafficking of large dense-core secretory vesicles. Our previous experimental and computational study found that the C2A domain of Slp-4 (also called granuphilin) binds with high affinity to anionic phospholipids in the cytoplasmic leaflet of the plasma membrane through a large positively charged protein surface centered on a cluster of phosphoinositide-binding lysine residues. Because this surface contributes greatly to Slp-4 C2A domain membrane binding, we hypothesized that the net charge on the surface might be evolutionarily conserved. To test this hypothesis, the known C2A sequences of Slp-4 among vertebrates were organized by class (from mammalia to pisces) using molecular phylogenetic analysis. Consensus sequences for each class were then identified and used to generate homology structures, from which Poisson-Boltzmann electrostatic potentials were calculated. For comparison, homology structures and electrostatic potentials were also calculated for the five human Slp protein family members. The results demonstrate that the charge on the membrane-binding surface is highly conserved throughout the evolution of Slp-4, and more highly conserved than many individual residues among the human Slp family paralogs. Such molecular phylogenetic-driven computational analysis can help to describe the evolution of electrostatic interactions between proteins and membranes which are crucial for their function.

A functional microbiome catalog crowdsourced from North American rivers.

Predicting elemental cycles and maintaining water quality under increasing anthropogenic influence requires understanding the spatial drivers of river microbiomes. However, the unifying microbial processes governing river biogeochemistry are hindered by a lack of genome-resolved functional insights and sampling across multiple rivers. Here we employed a community science effort to accelerate the sampling, sequencing, and genome-resolved analyses of river microbiomes to create the Genome Resolved Open Watersheds database (GROWdb). This resource profiled the identity, distribution, function, and expression of thousands of microbial genomes across rivers covering 90% of United States watersheds. Specifically, GROWdb encompasses 1,469 microbial species from 27 phyla, including novel lineages from 10 families and 128 genera, and defines the core river microbiome for the first time at genome level. GROWdb analyses coupled to extensive geospatial information revealed local and regional drivers of microbial community structuring, while also presenting a myriad of foundational hypotheses about ecosystem function. Building upon the previously conceived River Continuum Concept 1 , we layer on microbial functional trait expression, which suggests the structure and function of river microbiomes is predictable. We make GROWdb available through various collaborative cyberinfrastructures 2, 3 so that it can be widely accessed across disciplines for watershed predictive modeling and microbiome-based management practices.

Heterotrophic archaea contribute to carbon cycling in low-pH, suboxic biofilm communities.

Archaea are widely distributed and yet are most often not the most abundant members of microbial communities. Here, we document a transition from Bacteria- to Archaea-dominated communities in microbial biofilms sampled from the Richmond Mine acid mine drainage (AMD) system (∼pH 1.0, ∼38°C) and in laboratory-cultivated biofilms. This transition occurs when chemoautotrophic microbial communities that develop at the air-solution interface sink to the sediment-solution interface and degrade under microaerobic and anaerobic conditions. The archaea identified in these sunken biofilms are from the class Thermoplasmata, and in some cases, the highly divergent ARMAN nanoarchaeal lineage. In several of the sunken biofilms, nanoarchaea comprise 10 to 25% of the community, based on fluorescent in situ hybridization and metagenomic analyses. Comparative community proteomic analyses show a persistence of bacterial proteins in sunken biofilms, but there is clear evidence for amino acid modifications due to acid hydrolysis. Given the low representation of bacterial cells in sunken biofilms based on microscopy, we infer that hydrolysis reflects proteins derived from lysed cells. For archaea, we detected ∼2,400 distinct proteins, including a subset involved in proteolysis and peptide uptake. Laboratory cultivation experiments using complex carbon substrates demonstrated anaerobic enrichment of Ferroplasma and Aplasma coupled to the reduction of ferric iron. These findings indicate dominance of acidophilic archaea in degrading biofilms and suggest that they play roles in anaerobic nutrient cycling at low pH.

Host-derived oxidized phospholipids and HDL regulate innate immunity in human leprosy.

Intracellular pathogens survive by evading the host immune system and accessing host metabolic pathways to obtain nutrients for their growth. Mycobacterium leprae, the causative agent of leprosy, is thought to be the mycobacterium most dependent on host metabolic pathways, including host-derived lipids. Although fatty acids and phospholipids accumulate in the lesions of individuals with the lepromatous (also known as disseminated) form of human leprosy (L-lep), the origin and significance of these lipids remains unclear. Here we show that in human L-lep lesions, there was preferential expression of host lipid metabolism genes, including a group of phospholipases, and that these genes were virtually absent from the mycobacterial genome. Host-derived oxidized phospholipids were detected in macrophages within L-lep lesions, and 1 specific oxidized phospholipid, 1-palmitoyl-2-(5,6-epoxyisoprostane E2)-sn-glycero-3-phosphorylcholine (PEIPC), accumulated in macrophages infected with live mycobacteria. Mycobacterial infection and host-derived oxidized phospholipids both inhibited innate immune responses, and this inhibition was reversed by the addition of normal HDL, a scavenger of oxidized phospholipids, but not by HDL from patients with L-lep. The accumulation of host-derived oxidized phospholipids in L-lep lesions is strikingly similar to observations in atherosclerosis, which suggests that the link between host lipid metabolism and innate immunity contributes to the pathogenesis of both microbial infection and metabolic disease.

Glycoside hydrolase activities of thermophilic bacterial consortia adapted to switchgrass.

Industrial-scale biofuel production requires robust enzymatic cocktails to produce fermentable sugars from lignocellulosic biomass. Thermophilic bacterial consortia are a potential source of cellulases and hemicellulases adapted to harsher reaction conditions than commercial fungal enzymes. Compost-derived microbial consortia were adapted to switchgrass at 60°C to develop thermophilic biomass-degrading consortia for detailed studies. Microbial community analysis using small-subunit rRNA gene amplicon pyrosequencing and short-read metagenomic sequencing demonstrated that thermophilic adaptation to switchgrass resulted in low-diversity bacterial consortia with a high abundance of bacteria related to thermophilic paenibacilli, Rhodothermus marinus, and Thermus thermophilus. At lower abundance, thermophilic Chloroflexi and an uncultivated lineage of the Gemmatimonadetes phylum were observed. Supernatants isolated from these consortia had high levels of xylanase and endoglucanase activities. Compared to commercial enzyme preparations, the endoglucanase enzymes had a higher thermotolerance and were more stable in the presence of 1-ethyl-3-methylimidazolium acetate ([C2mim][OAc]), an ionic liquid used for biomass pretreatment. The supernatants were used to saccharify [C2mim][OAc]-pretreated switchgrass at elevated temperatures (up to 80°C), demonstrating that these consortia are an excellent source of enzymes for the development of enzymatic cocktails tailored to more extreme reaction conditions.

The Sorcerer II Global Ocean Sampling expedition: expanding the universe of protein families.

Metagenomics projects based on shotgun sequencing of populations of micro-organisms yield insight into protein families. We used sequence similarity clustering to explore proteins with a comprehensive dataset consisting of sequences from available databases together with 6.12 million proteins predicted from an assembly of 7.7 million Global Ocean Sampling (GOS) sequences. The GOS dataset covers nearly all known prokaryotic protein families. A total of 3,995 medium- and large-sized clusters consisting of only GOS sequences are identified, out of which 1,700 have no detectable homology to known families. The GOS-only clusters contain a higher than expected proportion of sequences of viral origin, thus reflecting a poor sampling of viral diversity until now. Protein domain distributions in the GOS dataset and current protein databases show distinct biases. Several protein domains that were previously categorized as kingdom specific are shown to have GOS examples in other kingdoms. About 6,000 sequences (ORFans) from the literature that heretofore lacked similarity to known proteins have matches in the GOS data. The GOS dataset is also used to improve remote homology detection. Overall, besides nearly doubling the number of current proteins, the predicted GOS proteins also add a great deal of diversity to known protein families and shed light on their evolution. These observations are illustrated using several protein families, including phosphatases, proteases, ultraviolet-irradiation DNA damage repair enzymes, glutamine synthetase, and RuBisCO. The diversity added by GOS data has implications for choosing targets for experimental structure characterization as part of structural genomics efforts. Our analysis indicates that new families are being discovered at a rate that is linear or almost linear with the addition of new sequences, implying that we are still far from discovering all protein families in nature.

Using inferred residue contacts to distinguish between correct and incorrect protein models.

MOTIVATION: The de novo prediction of 3D protein structure is enjoying a period of dramatic improvements. Often, a remaining difficulty is to select the model closest to the true structure from a group of low-energy candidates. To what extent can inter-residue contact predictions from multiple sequence alignments, information which is orthogonal to that used in most structure prediction algorithms, be used to identify those models most similar to the native protein structure? RESULTS: We present a Bayesian inference procedure to identify residue pairs that are spatially proximal in a protein structure. The method takes as input a multiple sequence alignment, and outputs an accurate posterior probability of proximity for each residue pair. We exploit a recent metagenomic sequencing project to create large, diverse and informative multiple sequence alignments for a test set of 1656 known protein structures. The method infers spatially proximal residue pairs in this test set with good accuracy: top-ranked predictions achieve an average accuracy of 38% (for an average 21-fold improvement over random predictions) in cross-validation tests. Notably, the accuracy of predicted 3D models generated by a range of structure prediction algorithms strongly correlates with how well the models satisfy probable residue contacts inferred via our method. This correlation allows for confident rejection of incorrect structural models. AVAILABILITY: An implementation of the method is freely available at http://www.doe-mbi.ucla.edu/services.

Assessing the reliability of scores produced by the substance abuse subtle screening inventory.

The Substance Abuse Subtle Screening Inventory (SASSI) is a 10 scale indirect screening instrument used to detect substance use disorders. The current meta-analytic study described reliability reporting practices across 48 studies involving the SASSI. Reliability generalization methods were then employed to evaluate typical score reliability for the screening measure. Results showed approximately 73% of studies did not report reliability estimates. Analysis of data from the remaining studies revealed adequate reliability for the total scale (alpha = .87) and face valid scales (FVA alpha = .88 and FVOD alpha = .92), but substantially lower reliability estimates for the indirect scales (range of alpha = .23-.65). The study's findings underscore the need for improved reliability reporting for the SASSI and suggest cautious use of the measure, especially its indirect scales, as an indicator of problematic substance use/abuse in clinical settings.

Uncovering the Diversity and Activity of Methylotrophic Methanogens in Freshwater Wetland Soils.

Wetland soils are one of the largest natural contributors to the emission of methane, a potent greenhouse gas. Currently, microbial contributions to methane emissions from these systems emphasize the roles of acetoclastic and hydrogenotrophic methanogens, while less frequently considering methyl-group substrates (e.g., methanol and methylamines). Here, we integrated laboratory and field experiments to explore the potential for methylotrophic methanogenesis in Old Woman Creek (OWC), a temperate freshwater wetland located in Ohio, USA. We first demonstrated the capacity for methylotrophic methanogenesis in these soils using laboratory soil microcosms amended with trimethylamine. However, subsequent field porewater nuclear magnetic resonance (NMR) analyses to identify methanogenic substrates failed to detect evidence for methylamine compounds in soil porewaters, instead noting the presence of the methylotrophic substrate methanol. Accordingly, our wetland soil-derived metatranscriptomic data indicated that methanol utilization by the Methanomassiliicoccaceae was the likely source of methylotrophic methanogenesis. Methanomassiliicoccaceae relative contributions to mcrA transcripts nearly doubled with depth, accounting for up to 8% of the mcrA transcripts in 25-cm-deep soils. Longitudinal 16S rRNA amplicon and mcrA gene surveys demonstrated that Methanomassiliicoccaceae were stably present over 2 years across lateral and depth gradients in this wetland. Meta-analysis of 16S rRNA sequences similar (>99%) to OWC Methanomassiliicoccaceae in public databases revealed a global distribution, with a high representation in terrestrial soils and sediments. Together, our results demonstrate that methylotrophic methanogenesis likely contributes to methane flux from climatically relevant wetland soils.IMPORTANCE Understanding the sources and controls on microbial methane production from wetland soils is critical to global methane emission predictions, particularly in light of changing climatic conditions. Current biogeochemical models of methanogenesis consider only acetoclastic and hydrogenotrophic sources and exclude methylotrophic methanogenesis, potentially underestimating microbial contributions to methane flux. Our multi-omic results demonstrated that methylotrophic methanogens of the family Methanomassiliicoccaceae were present and active in a freshwater wetland, with metatranscripts indicating that methanol, not methylamines, was the likely substrate under the conditions measured here. However, laboratory experiments indicated the potential for other methanogens to become enriched in response to trimethylamine, revealing the reservoir of methylotrophic methanogenesis potential residing in these soils. Collectively, our approach used coupled field and laboratory investigations to illuminate metabolisms influencing the terrestrial microbial methane cycle, thereby offering direction for increased realism in predictive process-oriented models of methane flux in wetland soils.

Complex Evolutionary History of Translation Elongation Factor 2 and Diphthamide Biosynthesis in Archaea and Parabasalids.

Diphthamide is a modified histidine residue which is uniquely present in archaeal and eukaryotic elongation factor 2 (EF-2), an essential GTPase responsible for catalyzing the coordinated translocation of tRNA and mRNA through the ribosome. In part due to the role of diphthamide in maintaining translational fidelity, it was previously assumed that diphthamide biosynthesis genes (dph) are conserved across all eukaryotes and archaea. Here, comparative analysis of new and existing genomes reveals that some archaea (i.e., members of the Asgard superphylum, Geoarchaea, and Korarchaeota) and eukaryotes (i.e., parabasalids) lack dph. In addition, while EF-2 was thought to exist as a single copy in archaea, many of these dph-lacking archaeal genomes encode a second EF-2 paralog missing key residues required for diphthamide modification and for normal translocase function, perhaps suggesting functional divergence linked to loss of diphthamide biosynthesis. Interestingly, some Heimdallarchaeota previously suggested to be most closely related to the eukaryotic ancestor maintain dph genes and a single gene encoding canonical EF-2. Our findings reveal that the ability to produce diphthamide, once thought to be a universal feature in archaea and eukaryotes, has been lost multiple times during evolution, and suggest that anticipated compensatory mechanisms evolved independently.

Microbial Biodiversity: Straight from the Dolphin's Mouth.

Advances in metagenomic sequencing and bioinformatics have vastly expanded our knowledge of microbial phylogenetic and functional diversity. In this issue, Dudek et al. show that shotgun metagenomic sequencing of a less-well-studied environment - dolphin gums - uncovers surprising novelty in the bacterial tree of life, underscoring the promise of future discovery.

Methanogenesis in oxygenated soils is a substantial fraction of wetland methane emissions.

The current paradigm, widely incorporated in soil biogeochemical models, is that microbial methanogenesis can only occur in anoxic habitats. In contrast, here we show clear geochemical and biological evidence for methane production in well-oxygenated soils of a freshwater wetland. A comparison of oxic to anoxic soils reveal up to ten times greater methane production and nine times more methanogenesis activity in oxygenated soils. Metagenomic and metatranscriptomic sequencing recover the first near-complete genomes for a novel methanogen species, and show acetoclastic production from this organism was the dominant methanogenesis pathway in oxygenated soils. This organism, Candidatus Methanothrix paradoxum, is prevalent across methane emitting ecosystems, suggesting a global significance. Moreover, in this wetland, we estimate that up to 80% of methane fluxes could be attributed to methanogenesis in oxygenated soils. Together, our findings challenge a widely held assumption about methanogenesis, with significant ramifications for global methane estimates and Earth system modeling.

The Database of Interacting Proteins: 2004 update.

The Database of Interacting Proteins (http://dip.doe-mbi.ucla.edu) aims to integrate the diverse body of experimental evidence on protein-protein interactions into a single, easily accessible online database. Because the reliability of experimental evidence varies widely, methods of quality assessment have been developed and utilized to identify the most reliable subset of the interactions. This CORE set can be used as a reference when evaluating the reliability of high-throughput protein-protein interaction data sets, for development of prediction methods, as well as in the studies of the properties of protein interaction networks.

Fungi Contribute Critical but Spatially Varying Roles in Nitrogen and Carbon Cycling in Acid Mine Drainage.

The ecosystem roles of fungi have been extensively studied by targeting one organism and/or biological process at a time, but the full metabolic potential of fungi has rarely been captured in an environmental context. We hypothesized that fungal genome sequences could be assembled directly from the environment using metagenomics and that transcriptomics and proteomics could simultaneously reveal metabolic differentiation across habitats. We reconstructed the near-complete 27 Mbp genome of a filamentous fungus, Acidomyces richmondensis, and evaluated transcript and protein expression in floating and streamer biofilms from an acid mine drainage (AMD) system. A. richmondensis transcripts involved in denitrification and in the degradation of complex carbon sources (including cellulose) were up-regulated in floating biofilms, whereas central carbon metabolism and stress-related transcripts were significantly up-regulated in streamer biofilms. These findings suggest that the biofilm niches are distinguished by distinct carbon and nitrogen resource utilization, oxygen availability, and environmental challenges. An isolated A. richmondensis strain from this environment was used to validate the metagenomics-derived genome and confirm nitrous oxide production at pH 1. Overall, our analyses defined mechanisms of fungal adaptation and identified a functional shift related to different roles in carbon and nitrogen turnover for the same species of fungi growing in closely located but distinct biofilm niches.

Perturbation and restoration of the fathead minnow gut microbiome after low-level triclosan exposure.

BACKGROUND: Triclosan is a widely used antimicrobial compound and emerging environmental contaminant. Although the role of the gut microbiome in health and disease is increasingly well established, the interaction between environmental contaminants and host microbiome is largely unexplored, with unknown consequences for host health. This study examined the effects of low, environmentally relevant levels of triclosan exposure on the fish gut microbiome. Developing fathead minnows (Pimephales promelas) were exposed to two low levels of triclosan over a 7-day exposure. Fish gastrointestinal tracts from exposed and control fish were harvested at four time points: immediately preceding and following the 7-day exposure and after 1 and 2 weeks of depuration. RESULTS: A total of 103 fish gut bacterial communities were characterized by high-throughput sequencing and analysis of the V3-V4 region of the 16S rRNA gene. By measures of both alpha and beta diversity, gut microbial communities were significantly differentiated by exposure history immediately following triclosan exposure. After 2 weeks of depuration, these differences disappear. Independent of exposure history, communities were also significantly structured by time. This first detailed census of the fathead minnow gut microbiome shows a bacterial community that is similar in composition to those of zebrafish and other freshwater fish. Among the triclosan-resilient members of this host-associated community are taxa associated with denitrification in wastewater treatment, taxa potentially able to degrade triclosan, and taxa from an unstudied host-associated candidate division. CONCLUSIONS: The fathead minnow gut microbiome is rapidly and significantly altered by exposure to low, environmentally relevant levels of triclosan, yet largely recovers from this short-term perturbation over an equivalently brief time span. These results suggest that even low-level environmental exposure to a common antimicrobial compound can induce significant short-term changes to the gut microbiome, followed by restoration, demonstrating both the sensitivity and resilience of the gut flora to challenges by environmental toxicants. This short-term disruption in a developing organism may have important long-term consequences for host health. The identification of multiple taxa not often reported in the fish gut suggests that microbial nitrogen metabolism in the fish gut may be more complex than previously appreciated.