Perturbation and restoration of the fathead minnow gut microbiome after low-level triclosan exposure
© Narrowe et al.; licensee BioMed Central. 2015
Received: 9 August 2014
Accepted: 29 January 2015
Published: 3 March 2015
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.
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.
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.
KeywordsTriclosan Fathead minnow 16S Gut microbiome Denitrification
The importance of host-associated gut microbiota to the normal development and overall health of the host organism is well-established and increasingly appreciated [1-3]. The gut microbiome has been shown to be both stable over the long-term  and vulnerable to disruption  which may have long-term implications for host health [6,7]. Studies on the chemical disruption of host-associated microbiota have generally been concerned with the clinical use of antibiotics [6,8,9], the intentional exposure to personal care products , or have focused on the effect of antimicrobials on specific taxa . Less well-characterized are the challenges to the healthy host-associated microbiome from common contaminants, including antimicrobial compounds, particularly at low, but environmentally relevant levels. Even short-duration, low-concentration exposures may alter the gut flora during developmentally important windows.
One such emerging contaminant is triclosan (5-chloro-2-(2,4-dichlorophenoxy)phenol), a chlorinated aromatic compound that has been used as an antimicrobial since the 1960s . Triclosan has limited clinical application, but is frequently found as a component of personal care products and household products . Triclosan is ubiquitous, enters the wastewater system, persists through the wastewater treatment process , and has been detected in surface waters at concentrations up to 2.3 μg L−1 . Triclosan has been associated with adverse physiological and developmental outcomes [16-18], is photo-degraded to produce dioxins , and is potentially an endocrine disruptor in fish . Recently, the United States Food and Drug Administration has reopened discussion of the regulation of triclosan , and the state of Minnesota has banned the sale of consumer products containing this compound . Despite this increased interest in triclosan, nothing is known about indirect effects on host health, either due to the triclosan-mediated alteration of the microbiome, or due to the microbially mediated transformation of triclosan.
We hypothesized that low, environmentally relevant levels of triclosan exposure are sufficient to disrupt the fish gut microbiome. To test this hypothesis, we exposed larval fathead minnows (Pimephales promelas) to two levels of triclosan (100 and 1,000 ng L−1) in a controlled laboratory experiment and used high-throughput 16S rDNA sequencing to profile gut microbial communities before exposure, immediately following an acute exposure window, and after depuration. We also sought to provide a census of the gut bacterial community of the untreated juvenile fathead minnow, an important environmental toxicology model organism .
16S V3-V4 variable region sequencing of fish gut microbiomes
Composition and accuracy of identification of mock community 16S rDNA sequences
Pipeline assigned taxonomy
Mock community ID (ATCC or GenBank)
k__Bacteria; p__Proteobacteria; c__Gammaproteobacteria; o__Alteromonadales; f__Shewanellaceae; g__Shewanella; s__
Shewanella amazonenesis SB2B ATCC: BAA-1098
k__Bacteria; p__Firmicutes; c__Bacilli; o__Lactobacillales; f__Lactobacillaceae; g__Pediococcus; s__
Pediococcus pentosaceus ATCC: 25745
k__Bacteria; p__Firmicutes; c__Bacilli; o__Bacillales; f__Bacillaceae; g__Bacillus; s__
Uncharacterized Bacillus isolate GenBank ID: KP025972
k__Bacteria; p__Firmicutes; c__Bacilli; o__Lactobacillales; f__Lactobacillaceae; g__Lactobacillus; s__zeae
Lactobacillus casei ATCC: 334
k__Bacteria; p__Firmicutes; c__Bacilli; o__Lactobacillales; f__Streptococcaceae; g__Lactococcus; s__
Lactococcus lactis SK11 ATCC: BAA-493 Lactococcus lactis IL403
k__Bacteria; p__Firmicutes; c__Bacilli; o__Lactobacillales; f__Lactobacillaceae; g__Lactobacillus; s__brevis
Lactobacillus brevis ATCC: 367
k__Bacteria; p__Proteobacteria; c__Gammaproteobacteria; o__Enterobacteriales; f__Enterobacteriaceae; g__; s__
Escherichia coli K12 GenBank ID: U00096
k__Bacteria; p__Actinobacteria; c__Actinobacteria; o__Actinomycetales; f__; g__; s__
Acidothermus cellulolyticus ATCC: 43068
k__Bacteria; p__Proteobacteria; c__Gammaproteobacteria; o__Enterobacteriales; f__Enterobacteriaceae; g__; s__
Escherichia coli K12 GenBank ID: U00096
k__Archaea; p__Euryarchaeota; c__Halobacteria; o__Halobacteriales; f__Halobacteriaceae; g__Halobacterium; s__
Halobacterium salinarum NRC-1 ATCC: 700922
V3-V4 sequencing of fish samples is accurate and technically reproducible
In order to assess the accuracy and reproducibility of V3-V4 hypervariable region library preparation and sequencing, we sequenced one mock community sample, three mock community spike-in control samples, and one replicate of sample D0C9. Community membership and abundance were highly correlated across replicate samples (Pearson r = 0.99 between D0C9 replicates; Additional file 4: Figure S4). We used the mock community sample to guide quality control and filtering of low-level contaminants and spurious OTUs. Our pipeline identified all ten expected mock community V3-V4 sequences at 100% identity (Table 1; Additional file 4: Figure S4) with the following exceptions: the two subspecies of Lactococcus lactis were 100% identical over the V3-V4 region surveyed and were collapsed together into OTU20. Escherichia coli K12 was represented by two OTUs. The E. coli K12 genome contains seven copies of the 16S rRNA gene, which are identical over the V3-V4 region. The second OTU assigned to E. coli was substantially less abundant (OTU294: mock community relative abundance 0.005) than the expected OTU (OTU30: mock community relative abundance 0.063) and may represent a contaminant or variant in the culture from which the DNA was acquired.
Our pipeline assigned 37 OTUs to our mock community sample. Of the reads from this sample, 99.89% were assigned to one of the ten expected mock community sequences, which were also the most abundant OTUs. The remaining 27 OTUs accounted for a total of only 103 reads (0.11%; Additional file 4: Figure S4). These 27 OTUs were all found within the fish samples. Although we cannot rule out the unlikely possibility of miscalled multiplexing indices falsely placing reads into the wrong sample, these 27 OTUs are presumed to be low-level PCR cross-contamination occurring during sequencing library preparation, reflecting the sensitivity inherent to deep DNA sequencing. The most abundant non-mock community OTU found in the mock community sample was OTU1, a member of the family Aeromonadaceae. OTU1 represented 32% of all reads in the unfiltered OTU table and was the most abundant OTU in the study. Within the mock community, this OTU accounted for only 0.03% of reads; this abundance level was used as the basis for study-wide OTU filtering (see the ‘Methods’ section and Additional file 4: Figure S4). If present in the experimental samples, the mock community OTUs were successfully filtered out as contaminants. In the spike-in samples, the few remaining non-filtered reads (mean n = 263) closely reconstructed the expected D0C9 community (Additional file 4: Figure S4; mean pairwise Pearson r = 0.85). Thus, we proceeded to estimate fish gut microbial community composition under the assumption that any observed inter-individual differences were not primarily due to technical artifacts.
The fathead minnow gut microbiome resembles that of other freshwater fish
Twenty OTUs form a study-wide core fathead minnow gut microbiome
Mean relative abundance
0.053 ± 0.057
0.017 ± 0.02
Acinetobacter johnsonii (species)
0.004 ± 0.006
0.37 ± 0.18
0.08 ± 0.11
0.02 ± 0.09
0.007 ± 0.012
0.008 ± 0.013
Cetobacterium somerae (species)
0.03 ± 0.04
CK-1C4-19 (candidate division)
0.03 ± 0.06
0.02 ± 0.04
0.01 ± 0.01
0.005 ± 0.007
0.03 ± 0.05
0.02 ± 0.03
Pseudomonas alcaligenes (species)
0.002 ± 0.003
0.03 ± 0.06
0.1 ± 0.06
0.001 ± 0.001
0.03 ± 0.04
Gut community composition changes over time
Environmental exposure to triclosan alters fish gut microbial community composition
Differential response to time and triclosan exposure can be seen for individual taxa
Fathead minnow GI tracts contain highly abundant members of the poorly characterized candidate division CK-1C4-19
Despite the low but environmentally relevant  levels (100 to 1,000 ng L−1) used here, acute triclosan exposure is sufficient to disrupt fish gut bacterial community structure (Figure 4). In addition, fish microbiomes from the two exposure levels are visually separated in the PCoA, which may indicate a dosage effect. Triclosan targets the final step of the bacterial fatty acid synthesis pathway, inhibiting the active site of the enoyl-acyl reductase FabI . Insensitivity to triclosan results from mutations in the fabI gene [35,36], FabI overexpression , the expression of alternate Fab proteins such as FabK, FabL, and FabV [37-39], or through efflux . Initially described in Vibrio sp., FabV is reported to confer greater resistance to triclosan than mutant forms of FabI . FabV has also been identified in other Gram-negative bacteria including members of the Aeromonas, Shewanella, and Pseudomonas genera, organisms highly abundant in freshwater fish gut microbiomes. The prevalence of fish gut bacterial species maintaining variant Fab genes may explain the stable presence of certain taxa across our study, such as Aeromonadaceae and Shewanella. In addition, triclosan may be either bactericidal or bacteriostatic depending on concentration . Finally, overall microbial community structure may contribute to differences in triclosan bioavailability and effect, if the community contains members that are resistant to or are able to metabolize triclosan.
Although order-level (Figure 5) or genus-level (Additional file 7: Figure S1) patterns of changes in abundance among exposure groups were clearly visible, many interesting differential responses to triclosan at the OTU-level were obscured by summarizing at higher taxonomic levels. For example, while as an aggregate, all members of the order Pseudomonadales are highly abundant over the course of this study (Figure 5); the abundance patterns of the five Acinetobacter OTUs in this order suggest differing OTU-level responses to environmental factors, including triclosan (Additional file 12: Figure S8). In addition to Acinetobacter, the 26 OTUs whose abundance increased at day 7 in association with triclosan exposure included Flavobacterium, Chryseobacterium, and Shewanella; all four genera contain a subset of species that are documented pathogens of fish and humans [42-45]. Another differentially abundant OTU was classified as candidate division CK-1C4-19, a rarely identified taxon previously documented within the zebrafish gut microbiota . Repeated detection of CK-1C4-19 within the guts of closely related hosts despite geographical distance and its apparent insensitivity to triclosan (Figure 8) warrant further study examining the potential functional contribution to or dependence on its host.
Also included among the OTUs significantly increased in abundance with triclosan exposure were species of Methylobacterium, Hydrogenophaga, and Acidovorax (Figure 6). Methylotrophic proteobacteria such as these are commonly observed as members of denitrifying consortia in wastewater treatment plants [46,47]. These consortia can include Thauera and Methylotenera [47-49], other genera identified in this study. The microbial denitrification process in wastewater treatment is often facilitated by amendment with carbon substrates including methanol  which can be utilized by methylotrophic denitrifiers such as those identified here . Over the course of this study, which used methanol for triclosan solubilization, OTUs identified as Hydrogenophaga, Methylophilaceae, Methylotenera, Thauera, and Acidovorax displayed a similar pattern of changes in relative abundance (Figure 7 and Additional file 11: Table S4). These OTUs were low in overall abundance at day 7 and became more prominent at later time points, particularly in triclosan-exposed and solvent-exposed samples. This pattern of increase in abundance suggests the possibility that methanol-enabled denitrification can occur within the fish gastrointestinal tract. Alternatively, because we did not measure changes in available nitrogen species in the water or characterize water microbial community changes over time, it is also possible that the combination of methanol (during the 7-day exposure) and the nitrogenous waste from the fish favored the growth of these organisms and that the measured increase in relative abundance of these taxa in the gut may in fact be due to a relative increase in transient bacteria that were sourced from the housing water. While these bacterial taxa were measured from the gut, this study did not differentiate between adherent or transient members of the gut microbial communities.
Hydrogenophaga spp. have been identified in association with amphibian hosts , but, to our knowledge, Thauera sp. have not been reported as members of vertebrate host-associated bacterial communities. Acidovorax OTUs have been noted in the trout gut  and also within the gut of the marine sea bream in a recent study that speculated on the possibilities of microbial denitrification occurring within the fish gut . The presence of annamox bacterial taxa within fish guts has also been shown within the carp gut . Host-associated denitrification and nitrous oxide production has been studied in freshwater invertebrates  and in earthworms , which also host Acidovorax species , but is less explored within vertebrate gut microbial communities, where removal of nitrate via reduction to ammonia rather than via denitrification is noted . Thus, microbial nitrogen cycling pathways within the fish gut may be more diverse than have been reported within mammalian guts. Because of the obvious limitations to attempting to infer function from taxonomic information , functional studies are needed, including direct measurement of denitrification-associated gene expression within fish gut microbiomes and housing or habitat waters.
While most studies performed to identify triclosan-degrading bacteria have focused on free-living rather than host-associated organisms, this catabolic ability has been demonstrated in bacteria across a wide variety of taxa , and it would not be surprising to find that members of the fish gut microbiome also possess the ability to metabolize or co-metabolize triclosan. Bacterial species related to those identified here have been associated with biodegradation of triclosan or other aromatic or halogenated compounds [47,61-66]. For example, OTU38 is a 100% BLAST match to a fully sequenced strain of Acidovorax that degrades poly-chlorinated biphenolic compounds , and Thauera has been specifically identified as a denitrifying organism with the potential to degrade aromatic compounds . Thauera and Hydrogenophaga were differentially abundant at day 14 and/or day 21 in high samples with respect to solvent samples (Additional file 10: Table S3), suggesting that these organisms may derive some competitive benefit from even a short exposure to triclosan, independent of the presence of the methanol solvent. Microbial degradation of halogenated aromatics typically requires an aerobic environment . Oxygen availability is likely to be variable over the length of the fish gastrointestinal tract, so sufficient oxygen may be present to support this process. If occurring in the gut, such microbially mediated processes could result in the direct intra-lumenal exposure to triclosan degradation products, including the lipophilic end product methyl-triclosan, which has been shown to bioaccumulate within fish .
Broadly, the differences in microbial community structure seen immediately following triclosan exposure (day 7) do not persist during depuration. By the conclusion of 2 weeks of recovery (day 21), as a whole, the communities cannot be distinguished on the basis of prior triclosan exposure, suggesting that most of the gut communities have returned to the same developmental path as those of the unexposed fish. Three samples from the triclosan-exposed groups remained distant from the main group of samples (Figure 4c), and in some cases, restoration of the initial community may take longer than our experimental window permits us to observe, or may never be complete, as has been shown for human microbiomes perturbed with antibiotic exposure . The apparent increase in alpha diversity associated with triclosan exposure (Figure 2) is difficult to interpret, as short-term antibiotic exposure has been shown to decrease [5,6] or increase  diversity in other vertebrate guts.
In the case of low concentrations of triclosan, as used here, a bacteriostatic effect on most taxa, rather than a bactericidal one, could result in a near-complete restoration of community structure after depuration; however, even short-term disruption to the gut ecology, as demonstrated here, may be harmful to a developing host, with the potential for both immediate and long-ranging effects. Typical fish commensals such as Aeromonas, Deefgea, and Flavobacterium sp. can also be pathogens [42,71,72] and may be held in check by less abundant, triclosan-sensitive members of the fish gut microbiome [24,73]. Thus, even a brief imbalance from acute exposure may precipitate opportunistic infections . The fish gut microbiome is implicated in nutrient absorption and growth [75,76], so juvenile dysbiosis as a result of environmental toxicants may impact long-term fitness at the individual or population level in contaminated habitats. In surface water ecosystems, increased time (constant rather than acute) and intensity (concentration) of contaminant exposure could easily result in permanent alteration of host microbiome, with potential for ecosystem-scale consequences.
Microbial communities of the fish gut are underexplored relative to the contribution of fish species to overall vertebrate diversity . The fathead minnow (P. promelas) is an important model organism for aquatic environmental toxicology  and is widely distributed in North America, and its developmental and reproductive response to environmental contaminants is well-characterized . The gut microbiome described here adds to the limited catalog of fish microbiomes characterized by high-throughput sequencing . Despite the temporal and triclosan exposure-associated differences in the microbial community structure, 20 of the 94 OTUs reported here constitute ‘core’ organisms present in 95% of the 103 samples sequenced (Table 2). These study-wide core OTUs include members of Aeromonadaceae, Bacteroidaceae, Shewanella, Pseudomonas, Deefgea, Acinetobacter, Flavobacterium, Cetobacterium, and others. Both technical differences (cultured vs. sequenced organisms, clones vs. short amplicons, choice of variable region sequenced, primer bias, sequencing depth, analysis) as well as differences related to the fish (age, diet, fresh caught vs. domesticated, water quality, husbandry, habitat) limit the ability to make direct comparisons of relative abundances across studies of other freshwater omnivorous fish gut microbiomes. Despite these caveats, these fathead minnow core microbiome OTUs are generally similar to the most common genera reported previously for zebrafish , the common carp  (both also family Cyprinidae), other freshwater omnivorous fish [24,25,80], and, in a limited fashion, for the fathead minnow . There are also some differences between the fathead minnow core microbiome and that of other related fishes. For example, the phylum Firmicutes occurs with notable relative abundance in some (but not all) gut communities reported for zebrafish and guppy , but is not present in the fathead minnow core microbiome or abundant at any point in the full experiment. Similarity of commensal bacterial communities among phylogenetically related and anatomically similar hosts has been shown for non-fish vertebrate species  and suggested for fish . Our study adds to the body of evidence that gut microbial community structure may also be conserved among phylogenetically related fishes.
The gut microbial community of individual fish is likely to be most strongly structured by the interacting effects of 1) host environment and diet, 2) host developmental stage, and 3) triclosan exposure history. The baseline communities were markedly different than days 7 to 21, and this difference made a strong contribution to the overall beta-diversity temporal trajectory (Figure 3a). While almost all OTUs present at days 7 to 21 are also present in the baseline (day 0) samples, because we did not sequence samples of water (autoclaved) or food (commercially prepared), we cannot definitively rule out water or food as a source of the few new OTUs we detected. All exposure groups received food and water from identical stock at each feeding or water change, so any new OTUs that might have arisen from these sources are expected to be evenly distributed across cohorts. The baseline samples were collected immediately upon arrival in our facility, and differing community structure more likely reflects differing conditions in the rearing facility from which the fish were acquired, as diet, water chemistry, and stress have been suggested to be predictive of fish gut microbiome structure [24,25,77,80]. However, even when the highly divergent baseline samples were removed from consideration, communities were in large part structured by time point/developmental stage (Figure 3b). Our study was designed to examine the effect of short-term triclosan exposure on the endogenous gut microbiota of developing, larval fish. Gut microbiomes of developing vertebrates, including fish, are dynamic, showing complex successional processes [82-85]. Thus, the temporally dynamic communities observed here, in addition to effects of triclosan exposure, likely respond to the combined effects of host developmental processes and initial changes to environment.
Most research on the impact of environmental contaminants has focused on aquatic animals and invertebrates  or on water, sediment, and soil microbial communities [51,87]. The effects of environmental contaminants on the host-associated microbiome are largely unexplored. This study demonstrates a shift in the fish gut bacterial community following a 7-day exposure to low, environmentally relevant levels of triclosan. Taxa whose relative abundances change with triclosan exposure include those potentially involved in nitrogen cycling and triclosan metabolism. Even short-duration disruption to the host microbiome such as that shown here may induce long-term effects on the host organism and larger ecosystem.
The effects of acute early-life stage exposure to environmentally relevant concentrations of triclosan on the composition of the gastrointestinal tract microbiome in larval fathead minnow (P. promelas) were evaluated in an immersion exposure experiment. Fathead minnow larvae (approximately 8 weeks post-hatch) were obtained from Aquatic Biosystems (Fort Collins, CO). At this developmental stage, fish were sexually undifferentiated or were undergoing differentiation and gender was not determined. The rearing facility classifies fish by hatch age, and the fish in this study were all hatched within a 4-day period. The fish are not isogenic, rather are the progeny of multiple breeding groups. Thus, sibship and age are untested variables in our study.
Eight randomly selected fish were sampled upon arrival as initial controls (baseline; day 0) as described below. The remaining fish were randomly assigned to experimental and control groups. For each group, 15 larvae were placed into each of two 4-L glass jars containing 1 L of solution. All glassware and aeration tubing was autoclaved prior to use, and food was introduced using sterile, single-use serological pipettes. A stock solution of sterile hatched brine shrimp (Hikari Bio-Pure Baby Brine Shrimp, Hayward CA) was prepared daily, and larvae were presented with 2 mL of this stock solution and allowed to feed ad libitum. We did not re-verify food sterility in-house or sequence a food-only sample, and thus cannot rule out the possibility that food contributed OTUs to our study; however, at each feeding, all fish were fed from the same stock solution. Our design and methods do not differentiate between transient (including any food-associated) bacteria and adherent bacteria. Test solutions for the 7-day acute exposure were a) control: moderately hard reconstituted water (MHRW) as defined by US Environmental Protection Agency (EPA) protocol #EPA-821-R-02-013 , b) solvent control: 0.0001 mg L−1 methanol in MHRW, c) low triclosan: 100 ng L−1 triclosan in MHRW, and d) high triclosan: 1,000 ng L−1 in MHRW. Triclosan (Sigma-Aldrich) was solubilized in methanol, and methanol concentrations were identical across solvent control and triclosan solutions. Methanol was chosen as the solvent for consistency with a recent study focusing on the direct physiological effects of triclosan on the fathead minnow . MHRW involves autoclaved (121°C, 15 psi, 30 min), sterilized water in addition to the additives described in EPA protocol #EPA-821-R-02-013. Our measured pH and alkalinity align with those of similar husbandry schemes and experimental designs, including those of the breeding facility. To utilize environmentally relevant concentrations of triclosan, we chose 100 and 1,000 ng/L which allowed us to work within the Kolpin et al. study median of 140 and max of 2,300 ng/L , while still working below the reported LC50 for P. promelas (260 μg/L at 96-h duration) .
Over a 7-day acute exposure, >90% of the test solution was replaced daily. Daily static renewal exchanges consisted of 100 μL aliquots of the pertinent spike aliquot dissolved into 1 L of MHRW. At each daily renewal, every cohort was moved into a freshly autoclaved jar. Following the 7-day acute exposure, larvae from all exposure groups were maintained in control conditions with daily renewals into MHRW and fresh jars until the end of the experiment on day 21 (Figure 1). Photoperiod (14 h light: 10 h dark), temperature (22 ± 1°C), and dissolved oxygen (>85% saturation) did not differ among exposure groups during the 21-day experiment. Animal care and handling was in accordance with the Institutional Animal Care and Use Committee of the University of Colorado Denver, #92514(05)1E.
Multiple intact gastrointestinal (GI) tracts were collected at four time points for each experimental group: prior to exposure (baseline), immediately after a 7-day exposure (day 7), after 1 week of post-exposure depuration (day 14), and finally at the end of the 21-day experiment (day 21). Fish were anesthetized in ice water prior to euthanasia by rapid decapitation. Body length and mass measurements of anesthetized fish were collected and are included in Additional file 1: Table S1. Freshly dissected GI tracts were placed into filter-sterilized PBS and frozen at −20°C until DNA extraction. Since the gastrointestinal organ is less developed in larval fish compared to adult fish, the transition between esophagus, stomach, and intestines was not distinguishable during dissection. Manipulation instruments were autoclaved, and instruments and surfaces were cleaned with ethanol after each specimen. All work was performed in a biosafety cabinet.
Total bacterial and host DNA was extracted using the PowerSoil DNA Isolation Kit (MO BIO Laboratories, Carlsbad CA), with modifications to the standard protocol noted below. Samples were thawed and transferred with 50 μL of PBS storage buffer to the bead tube, which was then vortexed at maximum speed for 10 min to disrupt the intact fish GI tracts. After the addition of solution C1, the tube was heated at 65°C for 10 min, followed by 10 min of vortexing at maximum speed. Following the C3 incubation step, the centrifugation time was increased from 1 to 2 min. Extracted DNA was stored at −20°C. Because these extractions contained varying, unknown ratios of host to microbial gDNA, we used the spectroscopic DNA quantification as rough estimates of extraction success, but were unable to use these measurements for normalization of microbial gDNA mass for downstream processing.
Marker gene (16S rDNA) amplification and sequencing
Amplicon sequencing of 108 samples targeting the V3-V4 hypervariable regions of the 16S rRNA gene was performed following the dual-indexing strategy of Kozich et al. . The 103 experimental samples included eight baseline samples and eight samples per triclosan exposure group for each of the three subsequent time-points (Figure 1), with the exception of ‘day 14 low’ which had only seven samples due to the failure of one sample in DNA extraction. In addition to these 103 samples, we sequenced one technical replicate of sample D0C9, one sample of a mock bacterial community consisting of equal mass of gDNA from ten different bacterial and archaeal species  (Table 1), and three mock/spike technical replicates containing a 1:1 combination (by gDNA mass) of the mock community and sample D0C9. D0C9 was chosen because sufficient gDNA was available, and a pilot study suggested that this sample was representative of other baseline samples. PCR reactions included 10 μL Q5 2× Hot Start Master Mix (New England Biolabs, Ipswich MA), 0.5 μL each primer (from 10 μM stock), 1 μL template DNA, and 8 μL nuclease-free water under the following conditions: 30 s at 98°C; 25 cycles of (10 s at 98°C, 15 s at 55°C, 20 s at 72°C); 2 min at 72°C; 4°C hold. Triplicate reactions per sample were combined and cleaned up using the Zymo Clean and Concentrate-5 kit (Zymo Research Corporation, Irvine CA), eluted into 17 μL nuclease-free water and quantified using the Qubit BR dsDNA kit (Life Technologies, Grand Island NY). Additional PCR reactions were performed as needed to generate the 10 ng of cleaned amplicon from each sample included in the pooled sequencing library. Sequencing was performed at the University of Colorado Denver Genomics and Microarray Core with a single lane of Illumina MiSeq using 2 × 251 bp paired end reads and V2 chemistry, with 8% PhiX added to the library. All reads and metadata are deposited in the Sequence Read Archive (SRA; http://www.ncbi.nlm.nih.gov/sra) under BioProject PRJNA257816 (SRA accession SRP045371).
Read preprocessing, OTU assignment, and OTU filtering
Demultiplexing was performed with CASAVA v. 1.8, and reads representing the PhiX or reads not matching indices were removed. The remaining reads were assigned to OTUs following the UPARSE pipeline (usearch v7.0.1090_i86linux32)  with the following non-default parameters: -fastq_mergepairs (-fastq_truncqual 3, -fastq_minmergelen 250); -fastq_filter (-fastq_maxee 1.0, -fastq_truncqual 10, -fastq_minlen 300); -usearch_global (-strand plus, -id 0.97) and custom scripting to accommodate large file sizes. The uchime_ref step was omitted, and the final OTUs were checked for chimeras using the DECIPHER web tool and the short sequences option . Taxonomy was assigned to OTUs in QIIME v 1.8.0 , using the RDP classifier v2.2  which was retrained against the Greengenes 13_8 rep set  that had been trimmed using PrimerProspector v. 1.0.1  to the V3-V4 region . A phylogenetic tree of the OTU sequences was constructed using FastTree .
We used the mock community sample to determine a threshold for filtering our OTU table for spurious or contaminant OTUs. Of the mock community reads, 99.89% were assigned to OTUs representing mock community 16S rRNA sequences. All OTUs in the remaining 0.11% coincided with sequences found in the fish samples, and the most abundant non-mock sequence seen in the mock community sample constituted 0.03% of the reads. We used this value to filter the whole-study OTU table according to the following criteria. To be retained for further analysis, an OTU should be found at greater than 0.03% relative abundance in at least N − 3 samples where N is the size of a time + exposure group (e.g., day 7, solvent: N = 8). This procedure should eliminate or greatly reduce instances of PCR contamination between samples, as we do not expect contamination to occur in a biologically meaningful pattern consistent with our experimental design. Note that OTUs can still occur at <0.03% relative abundance in any given sample. This conservative filtering procedure runs the risk of excluding naturally occurring, low abundance OTUs near the detection limits of the filtering protocol; however, this is not expected to impact our ability to answer the diversity-based question posed in this study regarding the effects of triclosan exposure. We explicitly tested this by using a variety of beta diversity distance metrics on the filtered and unfiltered OTU tables and by calculating alpha diversity on the filtered and unfiltered OTU tables. All analyses presented here are conducted on the filtered OTU table unless otherwise noted, with no additional normalizations or rarefying procedures . Both filtered and unfiltered OTU tables (Additional files 13 and 14) and OTU FASTA files (Additional files 15 and 16) are available.
Publicly available sequences classified as belonging to the taxon CK-1C4-19 were downloaded from the Silva rRNA database SSU Ref web release 117 . The V3-V4 region was extracted in silico from these 124 sequences using PrimerProspector v. 1.0.1 . Of the 124 sequences, 62 contained priming sites deemed likely to successfully amplify. Examination of sequences predicted not to amplify revealed many misclassified eukaryotes or sequences with low pintail scores, which may be chimeras. The remaining 62 sequences were combined with CK-1C4-19 sequences from this study and several Aeromonas sequences as an outgroup, aligned using MUSCLE v3.8.31 with default parameters , and an approximately-maximum-likelihood phylogenetic tree generated using FastTree 2.1.5 SSE3 .
Data visualization and statistical analyses were conducting using QIIME, R (http://www.r-project.org), phyloseq , vegan , and ggplot2 . Generation of taxa summary bar charts and core microbiome calculation was performed using QIIME v. 1.8.0. All distance measures and ordinations were calculated in R v. 3.1.0 using phyloseq v. 1.8.2, and vegan v. 2.0-10. Shannon’s diversity index was calculated using phyloseq, with ANOVA and Tukey’s HSD tests performed in R. We used PERMANOVA (vegan::adonis)  with the weighted UniFrac distance matrix in order to test the ability of multiple variables (time, triclosan exposure) to account for observed variance in inter-sample distances. To calculate the significance of clusters observed in principle coordinates analysis, we used MRPP (vegan::mrpp), a univariate analysis that compares mean within-group distance against the within-group distance of randomly permuted groups . MRPP was performed with 10,000 permutations on a pairwise weighted UniFrac distance matrix for all samples within each time point, with groups defined by triclosan exposure category. To identify differentially abundant OTUs, we used the DESeq2  package for R which has been extended to the analysis of microbial community data via phyloseq  with parameters: test = Ward, fit = local, P ≤ 0.005.
Availability of supporting data
The data sets supporting the results of this article are available in the NCBI Sequence Read Archive, BioProject PRJNA257816, SRA accession SRP045371 (http://www.ncbi.nlm.nih.gov/sra).
This work was funded by the University of Colorado Denver, College of Liberal Arts and Sciences. We thank Aaron N. Johnson for the microdissection assistance and advice. We thank Michael Wunder for the advice on statistical reporting. We thank Angelina Baroffio for the assistance with dissections and fish husbandry. We thank Katrina Diener, Ted Shade, and the University of Colorado Denver Genomics and Microarray Core for the sequencing and informatics assistance. MAL is supported by funding from the Saudi Arabian Cultural Mission.
- Relman DA. The human microbiome: ecosystem resilience and health. Nutr Rev. 2012;70 Suppl 1:S2–9.PubMed CentralPubMedGoogle Scholar
- Huttenhower C, Gevers D, Knight R, Abubucker S, Badger JH, Chinwalla AT, et al. Structure, function and diversity of the healthy human microbiome. Nature. 2012;486:207–14.Google Scholar
- Kovatcheva-Datchary P, Tremaroli V, Bäckhed F. The gut microbiota. In: Rosenberg E, DeLong EF, Lory S, Stackebrandt E, Thompson F, editors. The prokaryotes. Berlin, Heidelberg: Springer Berlin Heidelberg; 2013. p. 3–24.Google Scholar
- Faith JJ, Guruge JL, Charbonneau M, Subramanian S, Seedorf H, Goodman AL, et al. The long-term stability of the human gut microbiota. Science. 2013;341:1237439.PubMed CentralPubMedGoogle Scholar
- Dethlefsen L, Relman DA. Incomplete recovery and individualized responses of the human distal gut microbiota to repeated antibiotic perturbation. Proc Natl Acad Sci U S A. 2011;108(Suppl):4554–61.PubMed CentralPubMedGoogle Scholar
- Cho I, Yamanishi S, Cox L, Methé BA, Zavadil J, Li K, et al. Antibiotics in early life alter the murine colonic microbiome and adiposity. Nature. 2012;488:621–6.PubMed CentralPubMedGoogle Scholar
- Lozupone CA, Stombaugh JI, Gordon JI, Jansson JK, Knight R. Diversity, stability and resilience of the human gut microbiota. Nature. 2012;489:220–30.PubMed CentralPubMedGoogle Scholar
- Jernberg C, Löfmark S, Edlund C, Jansson JK. Long-term impacts of antibiotic exposure on the human intestinal microbiota. Microbiology. 2010;156(Pt 11):3216–23.PubMedGoogle Scholar
- Pérez-Cobas AE, Gosalbes MJ, Friedrichs A, Knecht H, Artacho A, Eismann K, et al. Gut microbiota disturbance during antibiotic therapy: a multi-omic approach. Gut. 2013;62:1591–601.PubMed CentralPubMedGoogle Scholar
- Kraglund F. Triclosan produces statistically significant reduction in plaque, gingivitis and caries but not clinically important benefit. Evid Based Dent. 2014;15:6–7.PubMedGoogle Scholar
- Syed AK, Ghosh S, Love NG, Boles BR. Triclosan promotes Staphylococcus aureus nasal colonization. MBio. 2014;5:e01015.PubMed CentralPubMedGoogle Scholar
- Halden RU. On the need and speed of regulating triclosan and triclocarban in the United States. Environ Sci Technol. 2014;48:3603–11.PubMed CentralPubMedGoogle Scholar
- Dann AB, Hontela A. Triclosan: environmental exposure, toxicity and mechanisms of action. J Appl Toxicol. 2011;31:285–311.PubMedGoogle Scholar
- Heidler J, Halden RU. Mass balance assessment of triclosan removal during conventional sewage treatment. Chemosphere. 2007;66:362–9.PubMedGoogle Scholar
- Kolpin DW, Furlong ET, Meyer MT, Thurman EM, Zaugg SD, Barber LB, et al. Pharmaceuticals, hormones, and other organic wastewater contaminants in U.S. streams, 1999–2000: a national reconnaissance. Environ Sci Technol. 2002;36(6):1202–11.PubMedGoogle Scholar
- Lankester J, Patel C, Cullen MR, Ley C, Parsonnet J. Urinary triclosan is associated with elevated body mass index in NHANES. PLoS One. 2013;8:e80057.PubMed CentralPubMedGoogle Scholar
- Fritsch EB, Connon RE, Werner I, Davies RE, Beggel S, Feng W, et al. Triclosan impairs swimming behavior and alters expression of excitation-contraction coupling proteins in fathead minnow (Pimephales promelas). Environ Sci Technol. 2013;47:2008–17.PubMed CentralPubMedGoogle Scholar
- Cherednichenko G, Zhang R, Bannister RA, Timofeyev V, Li N, Fritsch EB, et al. Triclosan impairs excitation-contraction coupling and Ca2+ dynamics in striated muscle. Proc Natl Acad Sci U S A. 2012;109:14158–63.PubMed CentralPubMedGoogle Scholar
- Raut SA, Angus RA. Triclosan has endocrine-disrupting effects in male western mosquitofish, Gambusia affinis. Environ Toxicol Chem. 2010;29:1287–91.PubMedGoogle Scholar
- Kuehn BM. FDA pushes makers of antimicrobial soap to prove safety and effectiveness. JAMA. 2014;311:234.PubMedGoogle Scholar
- Marty J, Hawj F, Hoffman JA, Eaton CA, Scalze B, Hortman M, et al. SF2192/HF2452. Minnesota State Legislature. 2014;277:10841.Google Scholar
- Ankley GT, Villeneuve DL. The fathead minnow in aquatic toxicology: past, present and future. Aquat Toxicol. 2006;78:91–102.PubMedGoogle Scholar
- Roeselers G, Mittge EK, Stephens WZ, Parichy DM, Cavanaugh CM, Guillemin K, et al. Evidence for a core gut microbiota in the zebrafish. ISME J. 2011;5:1595–608.PubMed CentralPubMedGoogle Scholar
- Nayak SK. Role of gastrointestinal microbiota in fish. Aquac Res. 2010;41:1553–73.Google Scholar
- Sullam KE, Essinger SD, Lozupone CA, O’Connor MP, Rosen GL, Knight R, et al. Environmental and ecological factors that shape the gut bacterial communities of fish: a meta-analysis. Mol Ecol. 2012;21:3363–78.PubMedGoogle Scholar
- Wong S, Waldrop T, Summerfelt S, Davidson J, Barrows F, Kenney PB, et al. Aquacultured rainbow trout (Oncorhynchus mykiss) possess a large core intestinal microbiota that is resistant to variation in diet and rearing density. Appl Environ Microbiol. 2013;79:4974–84.PubMed CentralPubMedGoogle Scholar
- Lochmann R, Phillips H, Xie L. Effects of a dairy-yeast prebiotic and water hardness on the growth performance, mineral composition and gut microflora of fathead minnow (Pimephales promelas) in recirculating systems. Aquaculture. 2011;320:76–81.Google Scholar
- Lozupone C, Lladser ME, Knights D, Stombaugh J, Knight R. UniFrac: an effective distance metric for microbial community comparison. ISME J. 2011;5:169–72.PubMed CentralPubMedGoogle Scholar
- Anderson MJ. A new method for non-parametric multivariate analysis of variance. Austral Ecol. 2001;26:32–46.Google Scholar
- Kuczynski J, Liu Z, Lozupone C, McDonald D, Fierer N, Knight R. Microbial community resemblance methods differ in their ability to detect biologically relevant patterns. Nat Methods. 2010;7:813–9.PubMed CentralPubMedGoogle Scholar
- Mielke PW, Berry KJ, Johnson ES. Multi-response permutation procedures for a priori classifications. Commun Stat - Theor Meth. 1976;A5:1409–24.Google Scholar
- Love MI, Huber W, Anders S. Moderated estimation of fold change and dispersion for RNA-seq data with DESeq2. 2014.Google Scholar
- DeSantis TZ, Hugenholtz P, Larsen N, Rojas M, Brodie EL, Keller K, et al. Greengenes, a chimera-checked 16S rRNA gene database and workbench compatible with ARB. Appl Environ Microbiol. 2006;72:5069–72.PubMed CentralPubMedGoogle Scholar
- Quast C, Pruesse E, Yilmaz P, Gerken J, Schweer T, Yarza P, et al. The SILVA ribosomal RNA gene database project: Improved data processing and web-based tools. Nucleic Acids Res. 2013;41:D590–6.PubMed CentralPubMedGoogle Scholar
- McMurry LM, Oethinger M, Levy SB. Triclosan targets lipid synthesis. Nature. 1998;394:531–2.PubMedGoogle Scholar
- Heath RJ, Rubin JR, Holland DR, Zhang E, Snow ME, Rock CO. Mechanism of triclosan inhibition of bacterial fatty acid synthesis. J Biol Chem. 1999;274:11110–4.PubMedGoogle Scholar
- Heath RJ, Rock CO. A triclosan-resistant bacterial enzyme. Nature. 2000;406:145–6.PubMedGoogle Scholar
- Heath RJ, Su N, Murphy CK, Rock CO. The enoyl-[acyl-carrier-protein] reductases FabI and FabL from Bacillus subtilis. J Biol Chem. 2000;275:40128–33.PubMedGoogle Scholar
- Massengo-Tiassé RP, Cronan JE. Vibrio cholerae FabV defines a new class of enoyl-acyl carrier protein reductase. J Biol Chem. 2008;283:1308–16.PubMedGoogle Scholar
- Chuanchuen R, Karkhoff-Schweizer RR, Schweizer HP. High-level triclosan resistance in Pseudomonas aeruginosa is solely a result of efflux. Am J Infect Control. 2003;31:124–7.PubMedGoogle Scholar
- Russell AD. Whither triclosan? J Antimicrob Chemother. 2004;53:693–5.PubMedGoogle Scholar
- Decostere A. Flavobacterium columnare infections in fish: the agent and its adhesion to the gill tissue. Verh K Acad Geneeskd Belg. 2002;64:421–30.PubMedGoogle Scholar
- Zamora L, Vela AI, Palacios MA, Domínguez L, Fernández-Garayzábal JF. First isolation and characterization of Chryseobacterium shigense from rainbow trout. BMC Vet Res. 2012;8:77.PubMed CentralPubMedGoogle Scholar
- Janda JM, Abbott SL. The genus Shewanella: from the briny depths below to human pathogen. Crit Rev Microbiol. 2014;40:293–312.PubMedGoogle Scholar
- Vignier N, Barreau M, Olive C, Baubion E, Théodose R, Hochedez P, et al. Human infection with Shewanella putrefaciens and S. algae: report of 16 cases in Martinique and review of the literature. Am J Trop Med Hyg. 2013;89:151–6.PubMed CentralPubMedGoogle Scholar
- Hoshino T, Terahara T, Tsuneda S, Hirata A, Inamori Y. Molecular analysis of microbial population transition associated with the start of denitrification in a wastewater treatment process. J Appl Microbiol. 2005;99:1165–75.PubMedGoogle Scholar
- Etchebehere C, Errazquin I, Barrandeguy E, Dabert P, Moletta R, Muxí L. Evaluation of the denitrifying microbiota of anoxic reactors. FEMS Microbiol Ecol. 2001;35:259–65.PubMedGoogle Scholar
- Lu H, Chandran K, Stensel D. Microbial ecology of denitrification in biological wastewater treatment. Water Res. 2014;64:237–54.PubMedGoogle Scholar
- Liu B, Zhang F, Feng X, Liu Y, Yan X, Zhang X, et al. Thauera and Azoarcus as functionally important genera in a denitrifying quinoline-removal bioreactor as revealed by microbial community structure comparison. FEMS Microbiol Ecol. 2006;55:274–86.PubMedGoogle Scholar
- Hallin S, Throbäck IN, Dicksved J, Pell M. Metabolic profiles and genetic diversity of denitrifying communities in activated sludge after addition of methanol or ethanol. Appl Environ Microbiol. 2006;72:5445–52.PubMed CentralPubMedGoogle Scholar
- Drury B, Scott J, Rosi-Marshall EJ, Kelly JJ. Triclosan exposure increases triclosan resistance and influences taxonomic composition of benthic bacterial communities. Environ Sci Technol. 2013;47:8923–30.PubMedGoogle Scholar
- Walke JB, Becker MH, Loftus SC, House LL, Cormier G, Jensen RV, et al. Amphibian skin may select for rare environmental microbes. ISME J. 2014;8:2207–17.PubMedGoogle Scholar
- Kormas KA, Meziti A, Mente E, Frentzos A. Dietary differences are reflected on the gut prokaryotic community structure of wild and commercially reared sea bream (Sparus aurata). Microbiologyopen. 2014;3:718–28.PubMed CentralPubMedGoogle Scholar
- Van Kessel MA, Dutilh BE, Neveling K, Kwint MP, Veltman JA, Flik G, et al. Pyrosequencing of 16S rRNA gene amplicons to study the microbiota in the gastrointestinal tract of carp (Cyprinus carpio L.). AMB Express. 2011;1:41.PubMed CentralPubMedGoogle Scholar
- Stief P, Poulsen M, Nielsen LP, Brix H, Schramm A. Nitrous oxide emission by aquatic macrofauna. Proc Natl Acad Sci U S A. 2009;106:4296–300.PubMed CentralPubMedGoogle Scholar
- Wüst PK, Horn MA, Henderson G, Janssen PH, Rehm BHA, Drake HL. Gut-associated denitrification and in vivo emission of nitrous oxide by the earthworm families Megascolecidae and Lumbricidae in New Zealand. Appl Environ Microbiol. 2009;75:3430–6.PubMed CentralPubMedGoogle Scholar
- Davidson SK, Stahl DA. Selective recruitment of bacteria during embryogenesis of an earthworm. ISME J. 2008;2:510–8.PubMedGoogle Scholar
- Lundberg JO, Weitzberg E, Cole JA, Benjamin N. Nitrate, bacteria and human health. Nat Rev Microbiol. 2004;2:593–602.PubMedGoogle Scholar
- Meyer F, Trimble W, Chang E, Handley K. Functional predictions from inference and observation in sequence-based inflammatory bowel disease research. Genome Biol. 2012;13:169.PubMed CentralPubMedGoogle Scholar
- Mechichi T, Stackebrandt E, Gad’on N, Fuchs G. Phylogenetic and metabolic diversity of bacteria degrading aromatic compounds under denitrifying conditions, and description of Thauera phenylacetica sp. nov., Thauera aminoaromaticasp. nov., and Azoarcus buckelii sp. nov. Arch Microbiol. 2002;178:26–35.PubMedGoogle Scholar
- Monferrán MV, Echenique JR, Wunderlin DA. Degradation of chlorobenzenes by a strain of Acidovorax avenae isolated from a polluted aquifer. Chemosphere. 2005;61:98–106.PubMedGoogle Scholar
- Lambo AJ, Patel TR. Biodegradation of polychlorinated biphenyls in Aroclor 1232 and production of metabolites from 2,4,4′-trichlorobiphenyl at low temperature by psychrotolerant Hydrogenophaga sp. strain IA3-A. J Appl Microbiol. 2007;102:1318–29.PubMedGoogle Scholar
- Sueoka K, Satoh H, Onuki M, Mino T. Microorganisms involved in anaerobic phenol degradation in the treatment of synthetic coke-oven wastewater detected by RNA stable-isotope probing. FEMS Microbiol Lett. 2009;291:169–74.PubMedGoogle Scholar
- Lolas IB, Chen X, Bester K, Nielsen JL. Identification of triclosan-degrading bacteria using stable isotope probing, fluorescence in situ hybridization and microautoradiography. Microbiology. 2012;158(Pt 11):2796–804.PubMedGoogle Scholar
- Lee DG, Chu K-H. Effects of growth substrate on triclosan biodegradation potential of oxygenase-expressing bacteria. Chemosphere. 2013;93:1904–11.Google Scholar
- Lee DG, Cho K-C, Chu K-H. Identification of triclosan-degrading bacteria in a triclosan enrichment culture using stable isotope probing. Biodegradation. 2014;25:55–65.Google Scholar
- Ohtsubo Y, Maruyama F, Mitsui H, Nagata Y, Tsuda M. Complete genome sequence of Acidovorax sp. strain KKS102, a polychlorinated-biphenyl degrader. J Bacteriol. 2012;194:6970–1.PubMed CentralPubMedGoogle Scholar
- Shapleigh JP. The denitrifying prokaryotes. In: Rosenberg E, DeLong EF, Lory S, Stackebrandt E, Thompson F, editors. The prokaryotes. Berlin, Heidelberg: Springer Berlin Heidelberg; 2013. p. 405–25.Google Scholar
- Leahy JG, Colwell RR. Microbial degradation of hydrocarbons in the environment. Microbiol Rev. 1990;54:305–15.PubMed CentralPubMedGoogle Scholar
- Balmer ME, Poiger T, Droz C, Romanin K, Bergqvist P-A, Müller MD, et al. Occurrence of methyl triclosan, a transformation product of the bactericide triclosan, in fish from various lakes in Switzerland. Environ Sci Technol. 2004;38:390–5.PubMedGoogle Scholar
- Janda JM, Abbott SL. The genus Aeromonas: taxonomy, pathogenicity, and infection. Clin Microbiol Rev. 2010;23:35–73.PubMed CentralPubMedGoogle Scholar
- Jung A, Jung-Schroers V. Detection of Deefgea chitinilytica in freshwater ornamental fish. Lett Appl Microbiol. 2011;52:497–500.PubMedGoogle Scholar
- Chu W, Zhou S, Zhu W, Zhuang X. Quorum quenching bacteria Bacillus sp. QSI-1 protect zebrafish (Danio rerio) from Aeromonas hydrophila infection. Sci Rep. 2014;4:5446.PubMed CentralPubMedGoogle Scholar
- Boutin S, Bernatchez L, Audet C, Derôme N. Network analysis highlights complex interactions between pathogen, host and commensal microbiota. PLoS One. 2013;8:e84772.PubMed CentralPubMedGoogle Scholar
- Li X, Yan Q, Xie S, Hu W, Yu Y, Hu Z. Gut microbiota contributes to the growth of fast-growing transgenic common carp (Cyprinus carpio L.). PLoS One. 2013;8:e64577.PubMed CentralPubMedGoogle Scholar
- Semova I, Carten JD, Stombaugh J, Mackey LC, Knight R, Farber SA, et al. Microbiota regulate intestinal absorption and metabolism of fatty acids in the zebrafish. Cell Host Microbe. 2012;12:277–88.PubMedGoogle Scholar
- Clements KD, Angert ER, Montgomery WL, Choat JH. Intestinal microbiota in fishes: what’s known and what’s not. Mol Ecol. 2014;23:1891–8.PubMedGoogle Scholar
- Vajda AM, Barber LB, Gray JL, Lopez EM, Bolden AM, Schoenfuss HL, et al. Demasculinization of male fish by wastewater treatment plant effluent. Aquat Toxicol. 2011;103:213–21.PubMedGoogle Scholar
- Llewellyn MS, Boutin S, Hoseinifar SH, Derome N. Teleost microbiomes: the state of the art in their characterization, manipulation and importance in aquaculture and fisheries. Front Microbiol. 2014;1:1.Google Scholar
- Wong S, Rawls JF. Intestinal microbiota composition in fishes is influenced by host ecology and environment. Mol Ecol. 2012;21:3100–2.PubMedGoogle Scholar
- Ley RE, Lozupone CA, Hamady M, Knight R, Gordon JI. Worlds within worlds: evolution of the vertebrate gut microbiota. Nat Rev Microbiol. 2008;6:776–88.PubMed CentralPubMedGoogle Scholar
- Palmer C, Bik EM, DiGiulio DB, Relman DA, Brown PO. Development of the human infant intestinal microbiota. PLoS Biol. 2007;5:1556–73.Google Scholar
- Koenig JE, Spor A, Scalfone N, Fricker AD, Stombaugh J, Knight R, et al. Succession of microbial consortia in the developing infant gut microbiome. Proc Natl Acad Sci U S A. 2011;108(Suppl):4578–85.PubMed CentralPubMedGoogle Scholar
- Kohl KD, Cary TL, Karasov WH, Dearing MD. Restructuring of the amphibian gut microbiota through metamorphosis. Environ Microbiol Rep. 2013;5:899–903.PubMedGoogle Scholar
- Vadstein O, Bergh Ø, Gatesoupe F-J, Galindo-Villegas J, Mulero V, Picchietti S, et al. Microbiology and immunology of fish larvae. Rev Aquac. 2013;5:S1–25.Google Scholar
- Chalew TE, Halden RU. Environmental exposure of aquatic and terrestrial biota to triclosan and triclocarban. J Am Water Works Assoc. 2009;45:4–13.PubMed CentralPubMedGoogle Scholar
- Butler E, Whelan MJ, Ritz K, Sakrabani R, van Egmond R. Effects of triclosan on soil microbial respiration. Environ Toxicol Chem. 2011;30:360–6.PubMedGoogle Scholar
- United States Environmental Protection Agency. Short-term methods for estimating the chronic toxicity of effluents and receiving water to freshwater organisms. EPA Publication No. 821-R-02–013. 4th ed. Washington, D.C: US EPA; 2002.Google Scholar
- Orvos DR, Versteeg DJ, Inauen J, Capdevielle M, Rothenstein A, Cunningham V. Aquatic toxicity of triclosan. Environ Toxicol Chem. 2002;21:1338–49.PubMedGoogle Scholar
- Kozich JJ, Westcott SL, Baxter NT, Highlander SK, Schloss PD. Development of a dual-index sequencing strategy and curation pipeline for analyzing amplicon sequence data on the MiSeq Illumina sequencing platform. Appl Environ Microbiol. 2013;79:5112–20.PubMed CentralPubMedGoogle Scholar
- Morgan JL, Darling AE, Eisen JA. Metagenomic sequencing of an in vitro-simulated microbial community. PLoS One. 2010;5:e10209.PubMed CentralPubMedGoogle Scholar
- Edgar RC. UPARSE: highly accurate OTU sequences from microbial amplicon reads. Nat Methods. 2013;10:996–8.PubMedGoogle Scholar
- Wright ES, Yilmaz LS, Noguera DR. DECIPHER, a search-based approach to chimera identification for 16S rRNA sequences. Appl Environ Microbiol. 2012;78:717–25.PubMed CentralPubMedGoogle Scholar
- Caporaso JG, Kuczynski J, Stombaugh J, Bittinger K, Bushman FD, Costello EK, et al. QIIME allows analysis of high-throughput community sequencing data. Nat Methods. 2010;7:335–6.PubMed CentralPubMedGoogle Scholar
- Wang Q, Garrity GM, Tiedje JM, Cole JR. Naive Bayesian classifier for rapid assignment of rRNA sequences into the new bacterial taxonomy. Appl Environ Microbiol. 2007;73:5261–7.PubMed CentralPubMedGoogle Scholar
- Walters WA, Caporaso JG, Lauber CL, Berg-Lyons D, Fierer N, Knight R. PrimerProspector: de novo design and taxonomic analysis of barcoded polymerase chain reaction primers. Bioinformatics. 2011;27:1159–61.PubMed CentralPubMedGoogle Scholar
- Werner JJ, Koren O, Hugenholtz P, DeSantis TZ, Walters WA, Caporaso JG, et al. Impact of training sets on classification of high-throughput bacterial 16 s rRNA gene surveys. ISME J. 2012;6:94–103.PubMed CentralPubMedGoogle Scholar
- Price MN, Dehal PS, Arkin AP. FastTree 2 - approximately maximum-likelihood trees for large alignments. PLoS One. 2010;5:e9490.PubMed CentralPubMedGoogle Scholar
- McMurdie PJ, Holmes S. Waste not, want not: why rarefying microbiome data is inadmissible. PLoS Comput Biol. 2014;10:e1003531.PubMed CentralPubMedGoogle Scholar
- Edgar RC. MUSCLE: multiple sequence alignment with high accuracy and high throughput. Nucleic Acids Res. 2004;32:1792–7.PubMed CentralPubMedGoogle Scholar
- McMurdie PJ, Holmes S. phyloseq: an R package for reproducible interactive analysis and graphics of microbiome census data. PLoS One. 2013;8:e61217.PubMed CentralPubMedGoogle Scholar
- Oksanen J, Blanchet FG, Kindt R, Legendre P, Minchin PR, O’Hara RB, et al. Package vegan. R Packag ver. 2013;254:20–8.Google Scholar
- Wickham H. Elegant graphics for data analysis, vol. 35. New York: Springer; 2009. p. 211.Google Scholar
This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/4.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly credited. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.