Environment shapes the fecal microbiome of invasive carp species
© The Author(s). 2016
Received: 8 February 2016
Accepted: 2 August 2016
Published: 12 August 2016
Although the common, silver, and bighead carps are native and sparsely distributed in Eurasia, these fish have become abundant and invasive in North America. An understanding of the biology of these species may provide insights into sustainable control methods. The animal-associated microbiome plays an important role in host health. Characterization of the carp microbiome and the factors that affect its composition is an important step toward understanding the biology and interrelationships between these species and their environments.
We compared the fecal microbiomes of common, silver, and bighead carps from wild and laboratory environments using Illumina sequencing of bacterial 16S ribosomal RNA (rRNA). The fecal bacterial communities of fish were diverse, with Shannon indices ranging from 2.3 to 4.5. The phyla Proteobacteria, Firmicutes, and Fusobacteria dominated carp guts, comprising 76.7 % of total reads. Environment played a large role in shaping fecal microbial community composition, and microbiomes among captive fishes were more similar than among wild fishes. Although differences among wild fishes could be attributed to feeding preferences, diet did not strongly affect microbial community structure in laboratory-housed fishes. Comparison of wild- and lab-invasive carps revealed five shared OTUs that comprised approximately 40 % of the core fecal microbiome.
The environment is a dominant factor shaping the fecal bacterial communities of invasive carps. Captivity alters the microbiome community structure relative to wild fish, while species differences are pronounced within habitats. Despite the absence of a true stomach, invasive carp species exhibited a core microbiota that warrants future study.
KeywordsFecal microbiome Invasive carp Microbiota Community structure 16SrRNA gene
Aquatic invasive species are among the greatest threats to aquatic ecosystems. In particular, species of carp, such as the common, bighead, and silver carp, can consume large quantities of food and disrupt food chains, while potentially out-competing native species and reaching great densities in invaded ranges [1–4]. The common carp is the most widespread invasive fish in the world. It was introduced to the USA over a century ago, and it has gradually spread throughout lakes and rivers. Once established, common carp act as ecosystem engineers, uprooting aquatic vegetation and increasing phosphorus availability, resulting in eutrophication and ecosystem degradation [5, 6]. In contrast, bighead and silver carps, members of a group known as the “Asian carps,” were recently introduced into the USA in the 1970s. Their range, however, is rapidly expanding and now stretches from the lower Mississippi River to its northern reaches and tributaries.
Although the potential ecosystem impacts of silver and bighead carp are not well known, as large filter feeders they could potentially affect native fish species directly by inducing changes in zooplankton communities [7, 8]. While the management of common carp has occasionally been successful through poisoning entire lakes or manipulating weaknesses in their life history characteristics (e.g., manipulation of predators in spawning habitats) , no control strategies have been successful in reducing the population of silver and bighead carps. Novel options for control must therefore come through a more thorough understanding of the biology of these invasive fishes.
The microbiome, the collection of microorganisms associated with an animal, is essential for optimal growth and survival of the host species [10–12]. In particular, the digestive tract microbiota plays an integral role in the breakdown of food, provision of energy, vitamin production, and shaping innate immunity. In humans and other mammals, the gut microbiota has been shown to have broad effects on health and behavior . In fish, the gut microbiota has largely been studied in the context of aquaculture in order to identify or examine the effect of probiotics to enhance growth or health [14–16]. A better understanding of the gut microbiome of fishes, however, might reveal potential for control of these invasive species, since dysbiosis of the gut microbiome has been found to contribute to disease manifestation in humans and other vertebrates . Manipulation of gut microbiota to influence health has recently received greater attention in humans through procedures such as fecal microbiota transplantation [18, 19]. Thus, characterizing the gut microbiome of invasive fishes is an important step toward understanding the community that comprises this “hidden organ,” which might be eventually exploited for species control purposes.
Fish possess a gut microbiome that is distinct from other animals and the microbial communities of water and soil . Previous studies have shown that the fish gut microbiome is dominated by members of the phyla Proteobacteria, Firmicutes, Bacteroidetes, Actinobacteria, and Fusobacteria , and both trophic level and salinity predominantly influence the fish gut microbial community [20–22]. While diet can also affect the gut microbiome, the significance and magnitude of the effect are variable [23–25]. The microbiota of prey items has been shown to influence the gut microbiome in three-spined stickleback; however, host genotype exhibited a larger effect . Gut microbiome diversity was inversely related with dietary diversity in two species of freshwater fishes , whereas the effect of diet on Trinidadian guppies was negligible . The gut microbiome can also reflect relative preference for cyanobacteria as a food source . In silver carp, the gut microbiome has also been shown to be geographically and temporally variable .
Like other vertebrates, fish likely harbor a core microbiome. Roeselers et al.  identified a core microbiome of zebrafish through comparison of lab-raised and wild stocks. Further support of this concept was demonstrated in a reciprocal transplant of microbiota between zebrafish and mice . After transplantation, the microbial community gradually shifted to resemble the typical structure of its new host. However, habitat changes, such as the transition from wild to captive environments can lead to dramatic changes in the gut microbiome of fishes, including decreased gut microbiome diversity [25, 28, 32].
Although our understanding of the structure of the fish microbiome has increased in recent years, there are still important gaps in our current knowledge regarding the factors that shape the fish gut microbiome. The advent of metagenomics and high-throughput amplicon sequencing technologies has demonstrated that culture-based studies of the fish microbiome are inherently biased and do not reflect total community diversity [14, 16]. In the first study of carp using high-throughput sequencing, van Kessel et al.  found that nearly half of the sequences in captive carp belonged to the phylum Fusobacteria, which were conspicuously absent in culture-dependent studies. However, of the studies that have utilized metagenomic approaches to examine the gut microbiome of carp species, only two have included wild fish [22, 29]. In addition, only one of these studies was done in an invaded range. Ye et al.  compared the microbiome of silver carp to gizzard shad, a native fish species that is planktivorous, but does not share the same taxonomic order. This study found differences in microbiota between fish species that could be explained by gut morphology and feeding preferences. Additional research is needed to compare the gut microbiome among closely related species of carps in their invaded range and to characterize the core microbiome.
In this study, we characterized the fecal microbiomes of three species of invasive carps (silver, bighead, and common carps) and determined the relationships between microbial community structure and environment, diet, and fish species. The aims of this study were to determine (1) the effect of environment on the microbiome of invasive carps, (2) how microbial communities of fish species differ within environments and the extent to which diet plays a role, (3) to what extent a core microbiome exists among invasive carps, and (4) how differences in microbiome structure among environments and fish species affect inferred bacterial community function.
Our results shed new light on the understanding of the microbiome of invasive carps and highlight the dominant role of the environment in shaping the fish microbiome.
Description of fishes used in this study
Algal feed mixture, see 
Algal feed mixture, see 
2.5-mm pellet feed (Oncor Fry, Skretting USA, Tooele, UT)
Frozen brine shrimp (San Francisco Bay Brand, Newark, CA)
Flake food (Color Tropical Marine, Pentair Aquatic Ecosystems, Apopka, FL)
Asian carps were field-identified as either bighead or silver carp, and fin clips were taken for SNP analysis to determine species and exclude F1 hybrids . An additional 13 common carp were caught from Lotus (44° 52′ 28′′ N, 93° 31′ 48′′ W), Riley (44° 50′ 10′′ N, 93° 31′ 17′′ W), and Susan (44° 51′ 5′′ N, 93° 32′ 27′′ W) lakes in MN, USA, in August and September 2012. Four fish each were collected from Lotus and Riley lakes, while five were collected from Lake Susan.
Laboratory-housed bighead carp (H. nobilis), silver carp (H. molitrix), common carp (C. carpio), and goldfish (Carassius auratus) were also used in this study (Table 1). Common carp and goldfish were obtained from a commercial fish hatchery (Osage Catfisheries Inc., Osage Beach, MO; Hunting Creek Fisheries, PA, respectively). Bighead and silver carps were obtained from an experimental research facility (US Geological Survey, Columbia, MO, USA). Goldfish were included for comparison because they are closely related and hypothesized to have originated from a wild population of Prussian carp (Carrassius gibelio) and a model organism. All fish were juvenile development stage and had been housed in the laboratory for at least 6 months.
All laboratory fishes were kept at 18–20 °C, within flow-through tanks with constant aeration. Fish were fed ad libitum once daily, and food source was consistent for at least 3 weeks prior to sampling (Table 1). Fish were sampled from at a minimum of three different tanks. All fish were held in accordance with the University of Minnesota’s Institutional Animal Care and Use Committee (IACUC) (Protocol: 1407-31659A).
Fecal specimens from wild fish were collected within 1 h of fish capture, whereas laboratory fishes were anesthetized in an aerated anesthetic bath (0.01 % MS-222; Syndel, CO, USA) prior to handling. Fecal specimens were collected by manual stripping of live or recently deceased fish. Fish were stroked firmly from pelvic fins to anus, and fecal material was collected in a sterile microcentrifuge tube and stored frozen at −20 °C.
DNA extraction, PCR, and sequencing
DNA was extracted from frozen fecal samples using the QIAamp DNA Stool Mini Kit (Qiagen, Hilden, Germany) according to the manufacturer’s protocol. The V6 hypervariable region of the 16S rRNA gene was amplified as previously described using a mixture of five forward primers and a barcoded reverse primer to amplify triplicate 50-μL reactions containing 25 ng of fecal DNA each . PCR products were visualized on a 2 % agarose gel and purified using the QIAquick Gel Extraction Kit (Qiagen, Hilden, Germany). Replicate purified reactions were pooled, and DNA was quantified using the QuantiFluor-ST and the dsDNA System (Promega, Madison, WI, USA). Purified amplicons were pooled in equal concentrations.
DNA was sequenced using an Illumina HiSeq 2000 platform at the University of Minnesota Genomics Center (Saint Paul, MN, USA) with up to 20 pooled samples per lane for a total of six runs. Paired-end sequences (100-bp read length) were joined as previously described . Sequencing results were submitted to NCBI under project number SRP071816.
Sequence reads were processed using the mothur software package version 1.36.1, as previously described [35, 36]. Sequence reads that had ambiguous bases, more than one mismatch to primer sequences, homopolymers > 8 nt, and quality scores < 35 in a 50-nt window were removed. UCHIME  was used to identify possible chimeric sequences in mothur using the default parameters and the SILVA database release 102 of bacterial reference sequences . Sequences identified as possible chimeras were subsequently removed from the dataset. Sequences were aligned using the SILVA database . The threshold for aligning the reverse compliment sequence was set to 0.75, and all other settings were set to default parameters. OTUs were clustered by furthest neighbor at a 97 % similarity cutoff. OTUs were classified using a naïve Bayesian classifier and the Ribosomal Database Project (RDP) taxonomic database release 9 and mother training set 9, with a probability cutoff of 60 % . Chloroplast sequences and sequences that were unclassified at the kingdom level were removed, which comprised 0.001 % of total reads. Sequences that had a frequency ≤ 10 were removed from the dataset , removing a total of 1.7 % of reads. The number of reads from each sample was normalized to 150,000 from a maximum of 1,841,796 by randomly subsampling.
The Shannon index of diversity (alpha diversity) was calculated for each experimental group (Table 1) using mothur. One-way ANOVA and Tukey post hoc tests were used to examine differences in OTUs observed (Sobs) and Shannon index of diversity using JMP, Version 10 (SAS Institute Inc., Cary, NC). The distances among environment-species groups were calculated based on Bray-Curtis dissimilarity . The resulting distances were used to perform a hierarchical clustering of experimental groups using Unweighted Pair Group Method (UPGMA). Results were graphed using TreeView v 1.6.6. Bray-Curtis dissimilarity distances were calculated among individual samples and used to perform ordination and statistical tests. Ordination was done using non-metric multidimensional scaling (NMDS) in two dimensions . NMDS was done separately for lab and wild fishes. Ten iterations were performed, and the iteration resulting in the lowest stress was plotted. Analysis of molecular variance (AMOVA) was used to test for significance of observed groupings by testing whether the genetic diversity within each group was different from the pooled genetic diversity.
Analysis of similarity (ANOSIM) and Kruskal-Wallis test were done in mothur. ANOSIM was also used to compare the community composition (beta diversity) among sampling groups. Kruskal-Wallis test was used to examine significant differences in abundance of individual OTUs among groups. Results are discussed at the taxonomic level of order, due to short sequence reads . All analyses were done at α = 0.05.
Putative microbiota functions were predicted using PICRUSt . OTUs were mapped to the Greengenes database version 13.5, and 85 % of genes were classified to a Tier 1 KO function. The weighted nearest sequenced taxon index (NSTI) scores averaged 0.056 ± 0.018, indicating a relatively good match to reference genomes (ideal NSTI ≤ 0.03). Functional predictions were assigned up to KO tier 2. Tier 1 KO were compared among fishes using one-way ANOVA and Tukey post hoc tests for wild and lab fishes separately. Student’s t test was used to compare KO between wild and lab fishes and between lake and river environments for common carp. Due to numerous significant differences among groups in tier 2 KO, data were visualized using PCA. Functional classifications of chitinases and vitamin B12 synthesis enzymes were compared between wild and laboratory-housed bighead carp using Student’s t test. All statistical analysis of functional data was done using JMP, Version 10 (SAS Institute Inc., Cary, NC).
Diversity and richness
Distinct and shared OTUs
There were many OTUs that differed among fish species within and between environments. Of the 13,793 OTUs identified in wild fish, 76.6 % differed in abundance across species (p ≤ 0.05). However, OTUs that differed significantly among species made up 94 to 99 % of total reads per species. Among the laboratory fishes examined, only 20.5 % of the 7262 OTUs identified varied among species, making up 78 to 86 % of reads.
Mean relative abundance of OTUs in invasive carp species across environments
Description, similarity and accession number, and source of closest NCBI Blast N match
Uncultured prokaryote, 98 % identity to KC601630, Asian seabass intestine
Psychrobacter sp., 100 % identity to EU753148, marine intertidal flat
Uncultured clone, 98 % identity to KM312603, earthworm gut
Cetobacterium sp., 98 % identity to KM85610, Zebrafish intestine
Uncultured prokaryote, 98 % identity to KC601623, Asian seabass intestine
Differences in community composition
ANOSIM showed significant differences among groups in community composition and abundance (beta diversity). All wild fish species had distinct microbial communities (p < 0.0001), and the gut microbiota between wild and lab fish of the same species were different (p < 0.0001). For lab fishes, however, fewer differences were observed. Bighead carp were not significantly different from common carp fed pellet feed (p = 0.08) or brine shrimp (p = 0.18) diet. Diet did not change the bacterial community composition of common carp (p = 0.30). All other comparisons among lab fish were significant at α = 0.05.
Axis loadings for principal components analysis of tier 2 KO functions
o Xenobiotics biodegradation and metabolism
o Metabolism of terpenoids and polyketides
o Lipid metabolism
o Signal transduction
o Cellular processes and signaling
o Metabolism of other amino acids
o Carbohydrate metabolism
o Nucleotide metabolism
o Replication and repair
o Lipid metabolism
o Metabolism of terpenoids and polyketides
o Xenobiotics biodegradation and metabolism
o Membrane transport
o Carbohydrate metabolism
o Cellular processes and signaling
o Amino acid metabolism
o Enzyme families
In a comparison of wild and laboratory-housed fish of the same species, there were no differences in tier 1 KO functions that were consistent across the three species. Comparisons of the proportion of proposed function for five vitamin B12 biosynthesis proteins between wild and lab fishes showed some differences. For instance, there was a higher proportion of vitamin B12 biosynthesis for two enzymes in wild silver carp and all five in common carp (p ≤ 0.05). Chitinase functional classification was not different (p = 0.20) between wild and laboratory bighead carp.
The role of diet and environment in shaping the fish microbiome
Bacterial community structure and function in the guts of carps was strongly affected by environment. There were large differences between laboratory-housed and wild fish of the same species and between common carp and lake and river environments. Many of the differences between wild and lab fishes can likely be partially explained in the context of diet, and the three taxonomic orders that varied most between wild and lab fish were Clostridiales, Fusobacteriales, Vibrionales, and Aeromonadales.
Clostridiales, which were more abundant in all wild-invasive carp species compared to those reared in the laboratory habitat, are associated with the degradation and metabolism of carbohydrates, specifically sugars [46, 47]. Clostridiaceae also aid in glucose fermentation as has been shown in earthworm guts  and are responsible for producing short chain fatty acids in vertebrates. Prebiotic arabinoxylan oligosaccharides have been shown to increase Clostridium spp. in sturgeon [23, 49]. Therefore, greater plant matter intake might increase the proportion of Clostridiales. Despite this supposition, however, both silver and bighead carps in the lab were fed a diet composed primarily of cyanobacteria and green algae. Moreover, feeding preference for cyanobacteria was similarly confirmed for wild silver carp ; thus, differences in plant matter intake cannot solely explain patterns of abundance of the Clostridiales.
An alternative explanation is that feed timing or availability in the wild carps resulted in changes in their gut microbiota. Transgenic common carp exhibited higher food intake and growth rate, which was associated with a greater proportion of Clostridiales in the gut microbiome . In animals and humans, this increase is thought to be due to enhanced production of the short chain fatty acids acetate, butyrate, and propionate . Although lab fish were fed to satiation, they ate only once daily during the daytime. Wild-invasive carps, in contrast, are able to consume food over a larger time frame, and they feed primarily at night [7, 51]. There is limited information on how food availability or timing affects fish gut microbiota; however, a study on Asian seabass showed starvation-induced changes in proportion of Bacteroides . While members of the genus Bacteroides were not abundant in invasive carps, it is unclear what effect food limitation might have on invasive carp species, and this topic warrants further study.
Laboratory-housed carp had a much higher proportion of Fusobacteriales compared to wild fish. Over 95 % of Fusobacteriales were classified to the genus Cetobacterium. Cetobacterium somerae, a bacterium within the order Fusobacteriales, is a common and widely distributed species within the guts of freshwater fishes, and its prevalence is negatively correlated with dietary availability of vitamin B12 (cobalamin) [53, 54]. Hence, the main role of Cetobacterium somerae in the fish gut is assumed to be synthesis of vitamin B12 [53, 54]. The relative proportion of vitamin B12 within each food type in this study is not known, so we are unable to make comparisons in the relative proportion of Fusobacteriales across diets. Bighead and silver carp can obtain B12 from algal food sources. Cyanobacteria are capable of synthesizing vitamin B12, and although many eukaryotic algae do not directly synthesize B12, they have symbioses with B12-producing bacteria [55, 56]. As wild common carp eat a significant amount of detritus, they may also satisfy their B12 requirements through consumption of bacteria . However, we found that wild common carp and silver carp had a higher proportion of vitamin B12 functional classifications than captive fish with a higher abundance of total Fusobacteriales.
For bighead carp, the bacterial order Vibrionales exhibited the greatest difference in abundance between wild and laboratory environments. These bacteria are may be dietary in origin. Vibrionales are associated with exoskeletons of zooplankton , a major food source for bighead carps. Although silver carp consume zooplankton, such as cladocerans and copepods, these organisms make up a smaller proportion of the silver carp diet compared to bighead carp . Dietary chitin has not been conclusively shown to increase the proportion of chitinase-producing bacteria in fish guts in aquaculture applications ; however, chitinase-producing bacteria have been isolated from many fishes, including common carp . We found no difference in the proportion of chitinase function between wild and captive bighead carp, further indicating that Vibrionales were not contributing significantly to chitinase activity in the gut.
The fecal microbiota of common carp from lake and river habitats were different with respect to bacterial community structure, richness, and diversity. Generally, all these parameters were greater in river fish than those dwelling in lakes. Common carp exhibit similar diets in both lake and river environments. The primary component of the stomach contents of both lake and river carp is detritus, while seeds and invertebrates make up the majority of the remaining contents [1, 2, 60]. Therefore, differences are likely due to environment-specific factors rather than diet alone.
In wild fishes, patterns in relative abundance of bacterial phyla reflected feeding preferences. While silver carp and bighead carp are both filter feeders, phytoplankton makes up a larger proportion of the silver carp diet . Common carp are omnivores, and their diet includes detritus, invertebrates, and plant matter , whereas freshwater drum are carnivorous [61, 62]
We found no consistent trends in phyla across trophic level, possibly due to the large variation in common carp between lake and river habitats.
We found a much greater proportion of unclassified bacteria than in previous studies [22, 29]. When considering the cumulative frequency of unclassified reads, most were present in high abundance (Additional file 2: Figure S1). Moreover, abundant sequences that were unclassified at the phylum level were distributed across species. For example, the three most abundant sequences are very prevalent across species (Additional file 3: Figure S2), and they are present in 88 to 99 % of individual fish sampled. While removal of sequence reads that were ≤ 10 in abundance decreased overall diversity, results of statistical analyses and beta diversity indices remained the same. Taken together, the abundance of bacteria unclassified at the phylum level and their distribution across fish species indicates that there is considerable diversity in the fish gut that is uncharacterized.
There are differences in the microbial community structure observed in the present study to previous studies on wild carps. For example, we did not observe a significant proportion of Bacteroidetes or Actinobacteria in the silver carp microbiome, as was observed by Ye et al. . The differences may be due to differences in the geographic location of the fish sampled or gut samples were taken. Generally, fish in the present study were sampled further north, and fecal material was collected by stripping rather than dissecting. Li et al.  found a predominance of Proteobacteria and comparatively few Firmicutes in bighead carp in their native range. However, samples were collected in the winter, when fish were under starvation conditions, which is known to alter gut microbial composition .
The distribution of phyla in captive carps was similar to that found in previous studies, with a few exceptions. In common carp, we found nearly 50 % of phyla classified to Fusobacteria, similar to van Kessel et al. . The proportion of Fusobacteria and Firmicutes in bighead carps was similar to Li et al.’s . However, we found a smaller proportion of Bacteroidetes than both studies.
Effect of captivity on the fish microbiome
Although patterns in microbiota in wild fish appeared to be linked with feeding preferences, we found that diet had little effect on fecal microbiome in a lab setting. For example, there was no difference in microbial community of common carp fed pellet or brine shrimp by any statistical measure of community composition. Previous studies have shown that diet can influence the fish gut microbiota by introduction of prey-associated microbes . However, we do not have any data in this study on the microbes present in the foods ingested by fishes or the microbiome of the surrounding habitat. Intraspecific changes in the gut microbiome may also result from alteration in food metabolism, as shown in Eastern African cichlids and surgeon fishes [25, 63]. However, changes in response to diet can be slight  and may depend on species  or ecotype . We also found that bighead and silver carp had different microbiota, despite being fed a similar algal feed mixture. Our results support those of Li et al.  who found that paddlefish and bighead carp reared in the same pond had distinct intestinal microbiota.
In lab fishes, species differences were not as pronounced as for wild fishes. For example, fewer OTUs differed in their relative abundance among lab fish species, and fewer differences among bacterial communities were observed. Temperature is a driving force in the biology and behavior of poikilotherms, such as fishes . In laboratory or captive environments, fishes are not exposed to diurnal cycles in water temperature, nor is there the level of habitat complexity that is found in the natural environment. Differences in habitat use between closely related species of cichlids within a lake were associated with differences in microbiota . Thus, differences in microbiota among species in the wild are likely due to a combination of the effects of diet, habitat usage, temperature, physiology, and taxonomy. Thus, the homogeneous and homeostatic environment of the laboratory might modulate behavioral effects on microbiota, such as habitat usage.
Only common carp exhibited a difference in alpha diversity between wild and captive fish, with lower diversity in captive fish. The effect of captivity on gut microbiome diversity may be species-dependent. Captivity did not reduce the diversity of zebrafish gut microbiome . However, captivity dramatically reduced in cichlid fish  and mummichog .
Our results indicate a striking effect of environment on the fecal microbiome of invasive carps and, in particular, a dramatic effect of captivity. In a previous study, few differences between the gut microbiota of wild and lab populations of zebrafish were observed . But this may be due to rearing practices as the wild zebrafish were held for approximately 1 month under laboratory conditions prior to sampling, and this time frame was previously shown to alter the gut microbial community of silver carp . Several studies have observed differences in microbiota of fishes in captive versus wild environments that mirror our findings. In a study of pond-reared and wild grass carp, Fusobacteria was more prevalent in pond fish . A predominance of Fusobacteria was also found for grass carp, crucian carp, and bighead carp held in a rearing pond and fed a commercial feed . Thus, we support the contention that environment shapes the fish gut microbiome and that a true understanding of their gut microbiota needs to come from wild-caught fish.
The core microbiome of invasive carps
Ringø and Birkbeck  enumerated five criteria required to be considered core gut microbiota in fishes: (1) they must be present in healthy individuals, (2) they colonize the gut at early life stages and persist throughout the lifespan of the fish, (3) they are found in both wild and cultured fish populations, (4) they are able to grow anaerobically, and (5) they are associated with the stomach, foregut, or hindgut. We identified five OTUs that comprise a large proportion of the core microbiome of silver, bighead, and common carps. These OTUs satisfied criteria 1, 3, and 5, but additional experiments would be needed to evaluate whether these OTUs satisfy all criteria.
Previous studies have shown that the proportion of bacteria that make up the core microbiome is variable . In zebrafish, 21 OTUs comprise the core microbiome . In rainbow trout, the core microbiota makes up over 80 % of the total community, and it is resistant to environmental factors . Future work is needed to further identify and characterize the core microbiota of carps. In addition, studies should assess when these bacteria initially colonize the gut and their degree of persistence over time. Moreover, due to the disproportionate abundance of some OTUs in bighead and silver carp, the potential for utilizing core microbiota to identify the presence of these invasive fishes in water bodies through the identification of their associated microbes be explored [70, 71].
Our results indicate that environment is an important factor controlling invasive carp fecal microbiota. We draw this conclusion from the difference between lab and wild fishes of the same species and the difference between lake and river habitats for common carp. Diet may partially explain some patterns in phyla abundance for wild fishes. However, in laboratory-housed fishes, diet did not exert a strong effect on fish gut microbiota, rather, fish species was the factor controlling differences among lab fishes. Future studies are needed to tease apart the multitude of factors which potentially control fish microbiota in wild populations. The mechanisms underlying differences in microbiota between lake- and river-dwelling common carp, for example, are not known. Due to the role of environment in shaping the microbiome, source-tracking markers for specific fish species should be developed from large representative samples of individuals from different geographical areas and habitat types. Studies on laboratory populations of fishes should be interpreted with caution, as lab fish have distinct microbiome structure from those of wild fishes, a pattern which is apparent, but not well characterized, in other studies. Invasive carps have a core fecal microbiome comprised primarily of five bacterial species, which make up approximately 40 % of the total fish gut microbiome. Future research is needed on the specific functional role of these organisms within the carp fecal microbiome. In addition, potential for augmenting these microbes for aquaculture applications or species detection or control should be evaluated.
AMOVA, analysis of molecular variance; ANOVA, analysis of variance; IACUC, Institutional Animal Care and Use Committee; KO, KEGG Orthology; NMDS, non-metric multidimensional scaling; NSTI, nearest sequenced taxon index; OTU, operational taxonomic unit; PCA, principal components analysis; RDP, Ribosomal Database Project; UPGMA, Unweighted Pair Group Method
We thank Kevin Irons, David Wyffels, and their crew from the Illinois Department of Natural Resources for the help in obtaining fish from the Illinois River. We also thank Mary Headrick, Brett Miller, and Reid Swanson for obtaining samples from Minnesota lakes, Ratna Ghosal and Aaron Claus for providing assistance with the collection of fecal samples from fishes, James Lamer for analyzing Asian carp DNA to determine hybridization, Chanlan Chun for assistance with data analysis, and Sendréa Best for help with the preparation of samples for sequencing. The authors acknowledge the Minnesota Supercomputing Institute (MSI) at the University of Minnesota for providing resources that contributed to the research results reported within this paper.
Funding for this project was provided by the Minnesota Environment and Natural Resources Trust Fund as recommended by the Legislative‐Citizen Commission on Minnesota Resources (LCCMR).
Availability of data and material
Sequencing results are available in the Sequence Read Archive (SRA) database at NCBI under BioProject ID PRJNA311219, accession number SRP071816.
PWS and MJS conceived of the study, participated in its design, and helped to draft the manuscript. MJH and CS helped to analyze the data and draft the manuscript. JJE collected and analyzed the data and drafted the manuscript. All authors read and approved the final manuscript.
The authors declare that they have no competing interests.
Consent for publication
Ethics approval and consent to participate
All animals used in this study were treated in accordance with the approved University of Minnesota’s Institutional Animal Care and Use Committee (IACUC) Protocol: 1407-31659A.
Open AccessThis article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. 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.
- Lubinsky KS, Van Vooren A, Farabee G, Janecek J, Jackson SD. Common carp in the Upper Mississippi River. Hydrobiologia. 2014;154:141–53.Google Scholar
- Crivelli AJ. The biology of the common carp, Cyprinus carpio L. in the Camargue, southern France. J Fish Biol. 1981;18:271–90.View ArticleGoogle Scholar
- Williamson CJ, Garvey JE. Growth, fecundity, and diets of newly established silver carp in the middle Mississippi River. Trans Am Fish Soc. 2005;134:1423–30.View ArticleGoogle Scholar
- Koehn JD. Carp (Cyprinus carpio) as a powerful invader in Australian waterways. Freshw Biol. 2004;49:882–94.View ArticleGoogle Scholar
- Bajer PG, Sorensen PW. Effects of common carp on phosphorus concentrations, water clarity, and vegetation density: a whole system experiment in a thermally stratified lake. Hydrobiologia. 2015;746:303–11.View ArticleGoogle Scholar
- Breukelaar AW, Lammens E, Breteler J, Tátrai I. Effects of benthivorous bream (Abramis brama) and carp (Cyprinus carpio) on sediment resuspension and concentrations of nutrients and chlorophyll a. Freshw Biol. 1994;32:113–21.View ArticleGoogle Scholar
- Kolar CS, Chapman DC, Courtenay WR, Housel CM, Williams JD, Jennings DP. Asian carps of the genus Hypophthalmichthys: A Biological Synopsis and Environmental Risk Assessment. Wahington, DC: US Fish and Wildlife Service; 2005.Google Scholar
- Sass GG, Hinz C, Erickson AC, McClelland NN, McClelland MA, Epifanio JM. Invasive bighead and silver carp effects on zooplankton communities in the Illinois River, Illinois, USA. J Great Lakes Res. 2014;40:911–21.View ArticleGoogle Scholar
- Bajer PG, Chizinski CJ, Silbernagel JJ, Sorensen PW. Variation in native micro-predator abundance explains recruitment of a mobile invasive fish, the common carp, in a naturally unstable environment. Biol Invasions. 2012;14:1919–29.View ArticleGoogle Scholar
- Human Microbiome Project Consortium. Structure, function, and diversity of the healthy human microbiome. Nature. 2012;486:207–14.View ArticleGoogle Scholar
- O’Hara AM, Shanahan F. The gut flora as a forgotten organ. EMBO Rep. 2006;7:688–93.View ArticlePubMedPubMed CentralGoogle Scholar
- Ray AK, Ghosh K, Ringø E. Enzyme-producing bacteria isolated from fish gut: a review. Aquac Nutr. 2012;18:465–92.View ArticleGoogle Scholar
- Sampson TR, Mazmanian SK. Control of brain development, function, and behavior by the microbiome. Cell Host Microbe. 2015;17:565–76.View ArticlePubMedPubMed CentralGoogle Scholar
- Ghanbari M, Kneifel W, Domig KJ. A new view of the fish gut microbiome: advances from next-generation sequencing. Aquaculture. 2015;448:464–75.View ArticleGoogle Scholar
- Nayak SK. Role of gastrointestinal microbiota in fish. Aquac Res. 2010;41:1553–73.View ArticleGoogle 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;5:1–17.View ArticleGoogle Scholar
- de Vos WM, de Vos EA. Role of the intestinal microbiome in health and disease: from correlation to causation. Nutr Rev. 2012;70:S45–56.View ArticlePubMedGoogle Scholar
- Khoruts A, Dicksved J, Jansson JK, Sadowsky MJ. Changes in the composition of the human fecal microbiome after bacteriotherapy for recurrent Clostridium difficile-associated diarrhea. J Clin Gastroenterol. 2010;44:354–60.PubMedGoogle Scholar
- Hamilton MJ, Weingarden AR, Unno T, Khoruts A, Sadowsky MJ. High-throughput DNA sequence analysis reveals stable engraftment of gut microbiota following transplantation of previously frozen fecal bacteria. Gut Microbes. 2013;4:125–35.View ArticlePubMedPubMed CentralGoogle Scholar
- Sullam KE, Essinger SD, Lozupone CA, O’Connor MP, Rosen GL, Knight R, Kilham SS, Russell JA. Environmental and ecological factors that shape the gut bacterial communities of fish: a meta-analysis. Mol Ecol. 2012;21:3363–78.View ArticlePubMedGoogle Scholar
- Wong SD, Rawls JF. Intestinal microbiota composition in fishes is influenced by host ecology and environment. Mol Ecol. 2012;21:3100–2.View ArticlePubMedPubMed CentralGoogle Scholar
- Li T, Long M, Gatesoupe FJ, Zhang Q, Li A, Gong X. Comparative analysis of the intestinal bacterial communities in different species of carp by pyrosequencing. Microb Ecol. 2014;69:25–36.View ArticlePubMedGoogle Scholar
- Geraylou Z, Souffreau C, Rurangwa E, D’Hondt S, Callewaert L, Courtin CM, Delcour JA, Buyse J, Ollevier F. Effects of arabinoxylan-oligosaccharides (AXOS) on juvenile Siberian sturgeon (Acipenser baerii) performance, immune responses, and gastrointestinal microbial community. Fish Shellfish Immunol. 2012;33:718–24.View ArticlePubMedGoogle Scholar
- Wong S, Waldrop T, Summerfelt S, Davidson J, Barrows F, Kenney PB, Welch T, Wiens GD, Snekvik K, Rawls JF, Good C. 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.View ArticlePubMed CentralGoogle Scholar
- Baldo L, Riera JL, Tooming-Klunderud A, Albà MM, Salzburger W. Gut microbiota dynamics during dietary shift in Eastern African cichlid fishes. PLoS One. 2015;10(5), e0127462.View ArticlePubMedPubMed CentralGoogle Scholar
- Smith CR, Snowberg LK, Caporaso G, Knight R, Bolnick DI. Dietary input of microbes and host genetic variation shape among-population differences in stickleback gut microbiota. ISME J. 2015;9:2515–26.View ArticlePubMedGoogle Scholar
- Bolnick DI, Snowberg LK, Hirsch PE, Lauber CL, Knight R, Caporaso JG, Svanbäck R. Individuals’ diet diversity influences gut microbial diversity in two freshwater fish (threespine stickleback and Eurasian perch). Ecol Lett. 2014;17:979–87.View ArticlePubMedPubMed CentralGoogle Scholar
- Sullam KE, Rubin BER, Dalton CM, Kilham SS, Flecker AS, Russell JA. Divergence across diet, time and populations rules out parallel evolution in the gut microbiomes of Trinidadian guppies. ISME J. 2015;9:1508–22.View ArticlePubMedPubMed CentralGoogle Scholar
- Ye L, Amberg J, Chapman D, Gaikowski M, Liu WT. Fish gut microbiota analysis differentiates physiology and behavior of invasive Asian carp and indigenous American fish. ISME J. 2014;8:541–51.View ArticlePubMedGoogle Scholar
- Roeselers G, Mittge EK, Stephens WZ, Parichy DM, Cavanaugh CM, Guillemin K, Rawls JF. Evidence for a core gut microbiota in the zebrafish. ISME J. 2011;5:1595–608.View ArticlePubMedPubMed CentralGoogle Scholar
- Rawls JF, Mahowald MA, Ley RE, Gordon JI. Reciprocal gut microbiota transplants from zebrafish and mice to germ-free recipients reveal host habitat selection. Cell. 2006;127:423–33.View ArticlePubMedPubMed CentralGoogle Scholar
- Clements KD, Angert ER, Montgmery L, Choat JH. Intestinal microbiota in fishes: what’s known and what’s not. Mol Ecol. 2014;23:1891–8.View ArticlePubMedGoogle Scholar
- van Kessel MAHJ, Dutilh BE, Neveling K, Kwint MP, Veltman JA, Flik G, Jetten MSM, Klaren PHM, Op Den Camp HJM. Pyrosequencing of 16S rRNA gene amplicons to study the microbiota in the gastrointestinal tract of carp (Cyprinus carpio L.). AMB Express. 2011;1:41.View ArticlePubMedPubMed CentralGoogle Scholar
- Lamer JT, Ruebush BC, Arbieva ZH, McClelland MA, Epifanio JM, Sass GG. Diagnostic SNPs reveal widespread introgressive hybridization between introduced bighead and silver carp in the Mississippi River Basin. Mol Ecol. 2015;24:3931–43.View ArticlePubMedGoogle Scholar
- Staley C, Unno T, Gould TJ, Jarvis B, Phillips J, Cotner JB, Sadowsky MJ. Application of Illumina next-generation sequencing to characterize the bacterial community of the Upper Mississippi River. J Appl Microbiol. 2013;115:1147–58.View ArticlePubMedGoogle Scholar
- Schloss PD, Westcott SL, Ryabin T, Hall JR, Hartmann M, Hollister EB, Lesniewski RA, Oakley BB, Parks DH, Robinson CJ, Sahl JW, Stres B, Thallinger GG, Van Horn DJ, Weber CF. Introducing mothur: open-source, platform-independent, community-supported software for describing and comparing microbial communities. Appl Environ Microbiol. 2009;75:7537–41.View ArticlePubMedPubMed CentralGoogle Scholar
- Edgar RC, Haas BJ, Clemente JC, Quince C, Knight R. UCHIME improves sensitivity and speed of chimera detection. Bioinformatics. 2011;27:2194–200.View ArticlePubMedPubMed CentralGoogle Scholar
- Pruesse E, Quast C, Knittel K, Fuchs BM, Ludwig W, Peplies J, Glockner FO. SILVA: a comprehensive online resource for quality checked and aligned ribosomal RNA sequence data compatible with ARB. Nucleic Acids Res. 2007;35:7188–96.View ArticlePubMedPubMed CentralGoogle 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.View ArticlePubMedPubMed CentralGoogle Scholar
- Brown SP, Veach AM, Rigdon-Huss AR, Grond K, Lickteig SK, Lothamer K, Oliver AK, Jumpponen A. Scraping the bottom of the barrel: are rare high throughput sequences artifacts? Fungal Ecology. 2015;13:221–225.Google Scholar
- Bray JR, Curtis JT. An ordination of the upland forest communities of southern Wisconsin. Ecol Monogr. 1957;27:325–349.Google Scholar
- Kenkel N, Orloci L. Applying metric and nonmetric multidimensional scaling to ecological studies: some new results. Ecology. 1986;67:919–28.View ArticleGoogle Scholar
- Mizrahi-Man O, Davenport ER, Gilad Y. Taxonomic classification of bacterial 16S rRNA genes using short sequencing reads: evaluation of effective study designs. PLoS One. 2013;8:18–23.View ArticleGoogle Scholar
- Langille MGI, Zaneveld J, Caporaso JG, McDonald D, Knights D, Reyes JA, Clemente JC, Burkepile DE, Vega Thurber RL, Knight R, Beiko RG, Huttenhower C. Predictive functional profiling of microbial communities using 16S rRNA marker gene sequences. Nat Biotechnol. 2013;31:814–21.View ArticlePubMedPubMed CentralGoogle Scholar
- Ringø E, Birkbeck TH. Intestinal microflora of fish larvae and fry. Aquac Res. 1999;30:73–93.View ArticleGoogle Scholar
- Clements KD, Pasch IBY, Moran D, Turner SJ. Clostridia dominate 16S rRNA gene libraries prepared from the hindgut of temperate marine herbivorous fishes. Mar Biol. 2007;150:1431–40.View ArticleGoogle Scholar
- Zhu L, Wu Q, Dai J, Zhang S, Wei F. Evidence of cellulose metabolism by the giant panda gut microbiome. Proc Natl Acad Sci. 2011;108:17714–9.View ArticlePubMedPubMed CentralGoogle Scholar
- Wüst PK, Horn MA, Drake HL. Clostridiaceae and Enterobacteriaceae as active fermenters in earthworm gut content. ISME J. 2011;5:92–106.View ArticlePubMedGoogle Scholar
- Geraylou Z, Souffreau C, Rurangwa E, Maes GE, Spanier KI, Courtin CM, Delcour JA, Buyse J, Ollevier F. Prebiotic effects of arabinoxylan oligosaccharides on juvenile Siberian sturgeon (Acipenser baerii) with emphasis on the modulation of the gut microbiota using 454 pyrosequencing. FEMS Microbiol Ecol. 2013;86:357–71.View ArticlePubMedGoogle 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.View ArticlePubMedPubMed CentralGoogle Scholar
- Garcia LM, Adelman IR. An in situ estimate of daily food consumption and alimentary canal evacuation rates of common carp, Cyprinus carpio L. J Fish Biol. 1985;27:487–93.View ArticleGoogle Scholar
- Xia J, Lin G, Fu G, Wan Z, Lee M, Wang L, Liu X, Yue G. The intestinal microbiome of fish under starvation. BMC Genomics. 2014;15:266.View ArticlePubMedPubMed CentralGoogle Scholar
- Sugita H, Miyajima C. The vitamin B12-producing ability of the intestinal microflora of freshwater fish. Aquaculture. 1991;92:267–76.View ArticleGoogle Scholar
- Tsuchiya C, Sakata T, Sugita H. Novel ecological niche of Cetobacterium somerae, an anaerobic bacterium in the intestinal tracts of freshwater fish. Lett Appl Microbiol. 2008;46:43–8.PubMedGoogle Scholar
- Croft MT, Lawrence AD, Raux-Deery E, Warren MJ, Smith AG. Algae acquire vitamin B12 through a symbiotic relationship with bacteria. Nature. 2005;438:90–3.View ArticlePubMedGoogle Scholar
- Rodionov DA. Comparative genomics of the vitamin B12 metabolism and regulation in prokaryotes. J Biol Chem. 2003;278:41148–59.View ArticlePubMedGoogle Scholar
- Huq A, Small EB, West PA, Huq MI. Rahman, Colwell RR. Ecological relationships between Vibrio cholerae and planktonic crustacean copepods. Appl Environ Microbiol. 1983;45:275–83.PubMedPubMed CentralGoogle Scholar
- Sampson SJ, Chick JH, Pegg MA. Diet overlap among two Asian carp and three native fishes in backwater lakes on the Illinois and Mississippi rivers. Biol Invasions. 2008;11:483–96.View ArticleGoogle Scholar
- Ringø E, Zhou Z, Olsen RE, Song SK. Use of chitin and krill in aquaculture - the effect on gut microbiota and the immune system: a review. Aquac Nutr. 2012;18:117–31.View ArticleGoogle Scholar
- Crivelli AJ. The destruction of aquatic vegetation by carp. Hydrobiologia. 1981;106:37–41.View ArticleGoogle Scholar
- Jacquemin SJ, Pyron M, Allen M, Etchison L. Wabash River freshwater drum Aplodinotus grunniens diet: effects of body size, sex, and river gradient. J Fish Wildl Manag. 2014;5:133–40.View ArticleGoogle Scholar
- Dettmers JM, Wahl DH, Soluk DA, Gutreuter S. Life in the fast lane: fish and foodweb structure in the main channel of large rivers. J N Am Benthol Soc. 2014;20:255–65.View ArticleGoogle Scholar
- Miyake S, Ngugi DK, Stingl U. Diet strongly influences the gut microbiota of surgeonfishes. Mol Ecol. 2015;24:656–72.View ArticlePubMedGoogle Scholar
- Givens CE. A fish tale: comparison of the gut microbiome of 15 fish species and the influence of diet and temperature on its composition, Doctoral thesis. Athens: University of Georgia; 2007.Google Scholar
- Li XM, Zhu YJ, Yan QY, Ringo E, Yang DG. Do the microbiota differ between paddlefish (Polyodon spathala) and bighead carp (Aristichthy nobilis) reared in the same pond? Appl Microbiol. 2014;117:1245–52.View ArticleGoogle Scholar
- Beitinger TL, Bennet WA, McCauley RW. Temperature tolerances of North American freshwater fishes exposed to dynamic changes in temperature. Environ Biol Fishes. 2000;58:237–75.View ArticleGoogle Scholar
- Franchini P, Fruciano C, Frickey T, Jones JC, Meyer A. The gut microbial community of midas cichlid fish in repeatedly evolved limnetic-benthic species pairs. PLoS One. 2014;9, e95027.View ArticlePubMedPubMed CentralGoogle Scholar
- Givens CE, Ransom B, Bano N, Hollibaugh JT. Comparison of the gut microbiomes of 12 bony fish and 3 shark species. Marine Ecol Prog Ser. 2015;518:209–23.View ArticleGoogle Scholar
- Ni J, Yan Q, Yu Y, Zhang T. Factors influencing the grass carp gut microbiome and its effect on metabolism. FEMS Microbiol Ecol. 2014;87:704–14.View ArticlePubMedGoogle Scholar
- Yan T, Sadowsky MJ. Determining sources of fecal bacteria in waterways. Environ Monit Assess. 2007;129:97–106.View ArticlePubMedGoogle Scholar
- Roslev P, Bukh AS. State of the aret molecular markers for fecal pollution source tracking in water. Appl Microbiol Biotechnol. 2011;89:1341–55.View ArticlePubMedGoogle Scholar
- Hansen A, Ghosal R, Caprio J, Claus AW, Sorensen PW. Anatomical and physiological studies of bigheaded carps demonstrate that the epibranchial organ functions as a pharyngeal taste organ. J Exp Biol. 2014;217:3945–54.View ArticlePubMedPubMed CentralGoogle Scholar