<figure><img src="image_1.png" /></figure>

Check for
updates

# Multi-Omics Analysis Reveals a Correlation between the Host Phylogeny, Gut Microbiota and Metabolite Profiles in Cyprinid Fishes

Tongtong Li<fnref n="1" />,²†, Meng Long<fnref n="1" />†, Huan Li², François-Joël Gatesoupe³, Xujie Zhang⁴, Qianqian Zhang<fnref n="1" />,<fnref n="5" />, Dongyue Feng⁶* and Aihua Li<fnref n="1" />,<fnref n="5" />*

Gut microbiota play key roles in host nutrition and metabolism. However, little is known about the relationship between host genetics, gut microbiota and metabolic profiles. Here, we used high-throughput sequencing and gas chromatography/mass spectrometry approaches to characterize the microbiota composition and the metabolite profiles in the gut of five cyprinid fish species with three different feeding habits raised under identical husbandry conditions. Our results showed that host species and feeding habits significantly affect not only gut microbiota composition but also metabolite profiles (ANOSIM, *p* ≤ 0.05). Mantel test demonstrated that host phylogeny, gut microbiota, and metabolite profiles were significantly related to each other (*p* ≤ 0.05). Additionally, the carps with the same feeding habits had more similarity in gut microbiota composition and metabolite profiles. Various metabolites were correlated positively with bacterial taxa involved in food degradation. Our results shed new light on the microbiome and metabolite profiles in the gut content of cyprinid fishes, and highlighted the correlations between host genotype, fish gut microbiome and putative functions, and gut metabolite profiles.

**Keywords:** correlation, host phylogeny, gut microbiota, metabolite profiles, cyprinid fishes

## INTRODUCTION

Vertebrates harbor vast and complex microbial communities that colonize their gastrointestinal tracts (Walter et al., 2011). As a result of this intimate relationship, the gut microbiome has become an integral part of the digestive system. The gut microbiota strongly influences fish health by stimulating the development of the intestinal epithelium and the immune system, and impeding pathogenic microorganisms to colonize the intestinal tract (Sugita et al., 1991; Flint et al., 2008; Ley et al., 2008b; Roeselers et al., 2011; Li T. et al., 2016). In the gastrointestinal tract, commensal bacteria can synthesize essential amino acids, vitamins and short-chain fatty acids (SCFAs), and contribute to feed efficiency, especially by degrading indigestible plant polysaccharides

**OPEN ACCESS**

**Edited by:**
Malka Halpern
University of Haifa, Israel

**Reviewed by:**
David William Waite,
University of Queensland, Australia
Renee Maxine Petri,
Veterinärmedizinische Universität,
Austria
Frank O'Neill Aylward,
Center for Microbial Oceanography:
Research and Education (C-MORE),
USA

**\*Correspondence:**

Dongyue Feng
fengdy76@sina.com
Aihua Li
liaihua@ihb.ac.cn

**Specialty section:**
This article was submitted to
Microbial Symbioses,
a section of the journal
Frontiers in Microbiology

Received: 16 December 2016

Accepted: 03 March 2017

Published: 17 March 2017

**Citation:**
Li T, Long M, Li H, Gatesoupe F-J, Zhang X, Zhang Q, Feng D and Li A (2017) Multi-Omics Analysis Reveals a Correlation between the Host Phylogeny, Gut Microbiota and Metabolite Profiles in Cyprinid Fishes. *Front. Microbiol.* 8:454. doi: 10.3389/fmicb.2017.00454

<footnote n="1">¹ State Key Laboratory of Freshwater Ecology and Biotechnology, Institute of Hydrobiology, Chinese Academy of Sciences, Wuhan, China, ² Key Laboratory of Environmental and Applied Microbiology, and Environmental Microbiology Key Laboratory of Sichuan Province, Chengdu Institute of Biology, Chinese Academy of Sciences, Chengdu, China, ³ Nutrition, Métabolisme et Aquaculture, Institut National de la Recherche Agronomique, University of Pau and Pays de l'Adour, Saint-Pée-sur-Nivelle, France, ⁴ College of Fisheries and Life Science, Shanghai Ocean University, Shanghai, China,</footnote>

<footnote n="5">⁵ Freshwater Aquaculture Collaborative Innovation Center of Hubei Province, Huazhong Agricultural University, Wuhan, China, ⁶ National Fisheries Technical Extension Center, Ministry of Agriculture, Beijing, China</footnote>

<footnote>† These authors have contributed
equally to this work.</footnote>
(Gill et al., 2006). Therefore, gut content and feces contain small molecules that are considered to result from co-metabolism or metabolic exchanges between microbes and host cells (Chen et al., 2012).

The composition of gut microbial communities is shaped by various internal and external factors, such as host genotype, diet, lifestyle, and surrounding environment (e.g., water temperature, salinity) (Nayak, 2010; Sullam et al., 2012). Previous studies of mammalian species revealed that their gut microbiota clustered according to diet rather than host phylogeny (Ley et al., 2008a; Muegge et al., 2011). It was also shown that identical twins had still significant differences in their gut microbiota, although they shared much higher similarity between gut microbiota structures than genetically unrelated married couples (Zoetendal et al., 2001). In addition, some scientists demonstrated that the phylogenetic relationships of hominids were completely consistent with those of gut microbiota (Ochman et al., 2010; Moeller et al., 2014). During evolution, changes in the composition of gut microbiota may lead to shifts in its functions, which may finally influence host nutrition and environmental adaptability (Amato, 2013). Comparative analysis among various hosts and their microbiota revealed that both diet and host phylogeny have driven the evolution of gut microbiota (Ley et al., 2008a). However, the main driving factor remains controversial. Thus, identifying shifts in gut microbiota composition and diversity over evolutionary timescales will be crucial to understanding how gut microbiota of cyprinid fishes is involved in evolution and environmental adaptation.

The concept of metabonomics was first defined as “the quantitative measurement of the dynamic multiparametric metabolic response of living systems to pathophysiological stimuli or genetic modification” (Nicholson et al., 1999; Nicholson and Wilson, 2003). Nowadays, metabonomics provides a systematic approach to characterize the metabolic phenotype, which results from a coordinated physiological response to various intrinsic and extrinsic parameters including environment, drugs, dietary, lifestyle, genetics, and microbiome (Rezzi et al., 2007; Li et al., 2008; Zhou et al., 2014). Recently, complementary metabonomic approaches have been employed for the biochemical characterization of metabolic changes triggered by gut microbiota, dietary variation, and stress interactions (Martin et al., 2008; Wikoffa et al., 2009; Vitali et al., 2010). Solid phase micro-extraction followed by gas chromatography and mass spectrometry (GC/MS) represents a novel method for studying metabolic profiles of biological samples, which was considered as the gold standard in metabonomics (Harrigan and Goodacre, 2012). With this method, the volatile compounds, or those that can be made volatile, or stabilized by derivatization, are separated by gas chromatography and then detected by mass spectrometry (Roessner et al., 2000; Vernocchi et al., 2012, 2016). This approach has been increasingly applied to humans and other terrestrial vertebrates (Garner et al., 2007; Vitali et al., 2010; Antharam et al., 2016; Faber et al., 2016). However, studies on fish gut metabolome are scarce compared with those on terrestrial vertebrates using GC/MS. Global metabolite profiling

performed here in on fish gut samples provided insights into the relationship between microbial populations and metabolites.

The polyculture of carps with different feeding habits is a traditional method to optimize the use of trophic resources in ponds (Li et al., 2015). In the present study, the intestinal microbial community structure and global metabolite profiles of five cyprinid fish species cohabitated in the same environment were investigated: herbivorous grass carp (Ctenopharyngodon idellus, HG) and blunt snout bream (Megalobrama amblycephala, HB); omnivorous crucian carp (Carassius auratus, OC); filter-feeding silver carp (Hypophthalmichthys molitrix, FS) and bighead carp (Hypophthalmichthys nobilis, FB). These five freshwater fishes are the major carps in Chinese aquaculture and widely cultivated for food. In 2012, the production of these species reached 14.48 million tons in China, accounting for about 62.02% of the total freshwater-cultured fish annual output (MoA, 2012). Furthermore, understanding the microbial community and metabolite profile in the gut of these fish species can provide useful information to improve health and productivity of these commercially valuable freshwater species.

To date, many studies focused on fish gastrointestinal microbiota, and such studies only concentrated on some factors (e.g., diet, feed habits, and genotype of host) that may affect fish gastrointestinal microbiota (Wu et al., 2012; Li et al., 2013, 2015; Ye et al., 2013). Little is known about the relationship among host evolutionary distance, gut microbiota and metabolic profiles in cyprinid fishes. Here, we investigated two key questions. First, do host genotype, gut microbiota and gut metabolic profiles correlate between them? Secondly, what is the relationship between dominant microbes and key metabolic products in the fish guts?

# MATERIALS AND METHODS

## Sample Collection and Pyrosequencing

All fish samples were harvested from Dongxihu Fish Farm, Wuhan City, Hubei Province, China, on April 10th, 2012, when the water temperature in the pond was 18°C. The adult cyprinid fishes were harvested from the same earth pond of approximately 0.5 ha. The fish were fed with commercial diet from JiuZhou Shennong Pharmaceutical Co., Ltd., Wuhan, China (crude protein ≥ 30.0%, crude fiber ≥ 12.0%, crude ash ≤ 15.0%, calcium = 0.4-2.5%, phosphate ≥ 0.7%, salt = 0.3-1.2%, moisture ≤ 12.5%, and lysine ≥ 1.2%). The five fish species with different feeding habits were as follows: herbivorous grass carp (Ctenopharyngodon idellus, HG) and blunt snout bream (Megalobrama amblycephala, HB); omnivorous crucian carp (Carassius auratus, OC); filter-feeding silver carp (Hypophthalmichthys molitrix, FS) and bighead carp (Hypophthalmichthys nobilis, FB). Three healthy individuals of each species were randomly harvested using nets. The average body weights were c. 1.2, 0.5, 0.3, 1.0, and 1.0 kg for grass carp, blunt snout bream, crucian carp, silver carp, and bighead carp, respectively. The fish were carried on ice to the laboratory within 2 h. Prior to dissection, all fish were euthanized with an overdose of MS 222 (Sigma, Germany). All procedures for
handling and euthanasia of fish were conducted as described by Wu et al. (2012). The intestines were aseptically removed from the abdominal cavity under sterile environments and the contents of hindgut (lower one-third of the full intestine) from the three individuals of each species were gently squeezed out, placed in a sterile tube and stored at -80°C. The methods used in this study were reviewed and approved by the ethics committee of the Institute of Hydrobiology, Chinese Academy of Sciences, and were carried out in accordance with the relevant guidelines, including any relevant details (Approval ID: keshuizhuan 08529).

DNA preparation, PCR amplification, and pyrosequencing were performed as described by Li et al. (2015). Briefly, total DNA was extracted from 250 mg of contents by using E.Z.N.A. Stool DNA kit (OMEGA, Bio-Tek, USA). A primer set (515F/926R) was used to amplify the V4 region of the 16S rRNA gene for analyzing gut microbiota. Each sample was amplified in triplicate in a reaction volume of 20 µL containing 1 × Ex Taq PCR buffer, 10 pM of each primer, 1.25 U Takara Ex Taq (all TaKaRa Biotechnology Co., Ltd., Dalian, China), and 5 ng genomic DNA using the following program: 5 min at 95°C, followed by 30 cycles of 94°C for 30 s, 56°C for 30 s, and 72°C for 30 s, and finally, 10 min at 72°C. PCR products were purified using the QIAEX II Gel Extraction Kit (QIAGEN). Library preparation and pyrosequencing was performed at the Chinese National Human Genome Center in Shanghai with 454 GS FLX sequencing platform (Roche). The raw sequences are available through the NCBI/EBI/DDBJ Sequence Read Archive (Accession No. DRA002627 and DRA001264¹).

## Phylogenetic Analysis of Fishes

Phylogenetic analysis of five cyprinid species (grass carp, blunt snout bream, crucian carp, silver carp, and bighead carp) were constructed from two mitochondrial genes (COI and cytb) using MEGA program (version 6.0) (Tamura et al., 2013). The corresponding gene sequences were downloaded from GenBank (accession numbers in Supplementary Table S1). Common rabbits (*Oryctolagus cuniculus*) was used as the outgroup taxon. The concatenated sequences of the two mitochondrial genes were aligned using ClustalW in MEGA (pairwise and multiple alignment parameters: gap opening penalty 15, gap extension penalty 6.66, delay divergent sequences 30%, DNA transition weight 0.5, and no use of a negative matrix) (Thompson et al., 1994; Tamura et al., 2013). The phylogenetic tree was constructed based on the aligned DNA sequences by the neighbor-joining method using Kimura 2-parameter model in MEGA (Test of phylogeny options: Bootstrap 999 replicates; Substitutions to include: Transitions + Transversions; Rates among sites: Uniform rates; Pattern among lineages: Same (Homogeneous); Gaps/Missing Data: Complete Deletion; Codon Positions: 1st + 2nd + 3rd + Non-coding) (Kimura, 1980; Tamura et al., 2013). The identical triplicates of the concatenated sequences were used to generate the phylogenetic distance matrix for Mantel tests in MEGA using the same methods as described above.

## Metabolite Extraction and Metabolite Profiling Analysis

Extraction of metabolites from gut content samples was according to Lisec et al. (2006) with some modification. Gas chromatography was performed on a HP-5MS capillary column [5% phenyl/95% methylpolysiloxane (30 m × 250 µm i.d., 0.25 µm film thickness, Agilent J&W Scientific, Folsom, CA, USA)] to separate the derivatives at a constant flow of 1 mL/min helium. Raw gas chromatography/mass spectrometry data (GC/MS) data were converted into CDF format (NetCDF) files by Agilent GC/MS 5975. Data Analysis was processed by the XCMS<fnref n="2" />. Detailed protocols for metabolite extraction and GC/MS analyses were provided in Supporting Information Methods.

## Pyrosequencing Data Processing and Bioinformatics Analysis

Quality filtering, denoising, and chimera checking of the sequences obtained from pyrosequencing were performed with the mothur software (Schloss et al., 2009), as described previously (Li et al., 2015). The average read length was 250 bp. The resulting sequences were rarefied separately by random subsampling of 5,935 sequences from each sample in mothur, thereby equalizing the sampling effort across samples (Schloss et al., 2009). To predict the functional profiles of microbial communities, OTUs were picked using a closed reference (Greengenes ver. 13.5) at 97% sequence similarity, with normalization to control for differences in 16S rRNA copy number among OTUs, and their metagenomic contributions were predicted using PICRUSt, from the Kyoto Encyclopedia of Genes and Genomes (KEGG) pathways (DeSantis et al., 2006; Kanehisa et al., 2012; Langille et al., 2013). The pertinence of the metagenome predictions was assessed by computing NSTI (Nearest Sequenced Taxon Index).

Non-metric Multidimensional Scaling (NMDS) plots based on Bray-Curtis distance were used to visualize the separation of microbiota structure, metabolic profiles, and presumptive functions across different groups. Statistical testing among variation in microbial community composition was carried out using the analysis of similarity (ANOSIM). Mantel tests based on Bray-Curtis distance and Spearman's correlation analysis were applied to evaluate the correlations between gut microbial dissimilarity, metabolic profiles and host evolutionary distance. Differences between groups were evaluated by One-way Analysis of Variance (One-way ANOVA) followed by Tukey's post hoc test. The heatmap was constructed using the heatmap 2 function of the R gplots package (Gentleman et al., 2004). The machine learning algorithm Random Forests (Breiman, 2001) was used to identify the discriminatory OTUs and presumptive functions between species using the package R randomForest with 5000 trees and all default settings. The Boruta algorithm was used to select the features with predictive power (Kursa and Rudnicki, 2010). In order to understand the relationship between some specific gut microbes and certain metabolites, redundancy analysis (RDA) was performed. All statistical

<footnote n="1">¹http://trace.ddbj.nig.ac.jp/DRASearch/</footnote>

<footnote n="2">²www.bioconductor.org</footnote>
analyses were performed using R version 3.1.0. All permutation tests (i.e., ANOSIM and Mantel test) were conducted using 999 permutations.

# RESULTS

## Differences in Bacterial Community Diversity and Metabolite Profiles

The completeness of sampling was estimated with Good's coverage, by calculating the probability that a randomly selected amplicon sequence was already detected in the same sample (Table 1). The coverage ranged from 93.29 to 99.63% (97.04 ± 2.64 %), indicating that between 15 and 270 [1/(1-'Good's coverage')] additional reads would need to be sequenced before detecting a new OTU. This level of coverage suggested that the majority of bacterial OTUs present in the samples were identified in this study.

The gut microbiota of filter-feeding silver carp and bighead carp samples yielded the significantly higher alpha-diversity indices (ANOVA, *p* < 0.05), followed by those of herbivorous grass carp, blunt snout bream, and omnivorous crucian carp (Table 1 and Supplementary Figure S1).

Similarities of the bacterial communities and metabolite profiles between samples were compared by ANOSIM and NMDS based on Bray-Curtis (Singh et al., 2015). ANOSIM revealed significant differences in the structure (Figure 1A; ANOSIM, *r* = 0.59, *p* = 0.001) of gut microbiota among different fish species. There was a tendency for the bacterial profiles of fish to separate by feeding habits (ANOSIM, *r* = 0.60, *p* = 0.001). In addition to differences in bacterial communities, the five species of fish had also marked differences in metabolite profiles (Figure 1B; ANOSIM, *r* = 0.60, *p* = 0.001). There was a

consistent effect of feeding habits on the metabolite profiles of fish (ANOSIM, *r* = 0.74, *p* = 0.001).

Moreover, the distances of the bacterial communities among the individuals between different feeding habits were calculated (Supplementary Figure S2). We found that distances between filter-feeding and herbivorous carps (fil vs. her), as well as between filter-feeding and omnivorous carps (fil vs. omn) were significantly higher as compared with those between herbivorous and omnivorous species (her vs. omn) (Supplementary Figure S2A; ANOVA, *p* < 0.05). For the herbivorous and omnivorous species, the distances among the individuals within the same feeding habits (HG vs. HB) were significantly lower than those in different feeding habits (HG vs. OC and HB vs. OC) (Supplementary Figure S2B; ANOVA, *p* < 0.05). Similar results were observed for the profiles of metabolites (Supplementary Figures S2C,D). The filter-feeding carps displayed greater within-group variability in gut microbiota and metabolic profiles, compared to the herbivorous species (Supplementary Figures S3A,B).

The NMDS plot indicated that the gut bacterial communities in each species clustered firstly together, and then they tend to cluster with those fish species within the same feeding habits (Figures 1A,B). Similar clusters were observed for metabolite profiling (Figures 1C,D).

## Variation in Gut Microbiota Composition, Gut Metabolites and Evolutionary Distance

The abundances of 10 major phyla were commonly observed across all the samples (Figure 2). Those sequences that were not assigned to known microbial phyla were designated as "Unclassified Bacteria," which represented 1.08% of the entire data set. A few phyla occurred at low abundance and sporadically in some samples were referred to as "Others" (≈0.24% of the total sequences). Fusobacteria was the most dominant phylum in the gut bacterial communities, representing 59.96% of the total sequences. This phylum was the highly dominant, especially in the gut samples of omnivorous crucian carp (OC1-OC3), followed by herbivorous grass carp (HG1HG3) and blunt snout bream (HB1-HB3), and then filter-feeding bighead carp (FB1-FB3) and silver carp (FS1-FS3). Proteobacteria and Planctomycetes were two major phyla of the gut content samples in filter-feeding silver carp and bighead carp. Bacteroidetes and Firmicutes were also the dominant phyla, except in crucian carp.

The two most predominant bacteria in the gut microbiota of crucian carp were two OTUs assigned to genera *Cetobacterium* (97.09% ± 0.36%) and *Aeromonas* (1.45% ± 0.83%) in the phyla Fusobacteria and Proteobacteria, respectively (Figure 3A). The genus *Cetobacterium* was also abundant in the gut content of the other fishes: in herbivorous blunt snout bream (89.87% ± 0.66%) and grass carp (63.80% ± 12.02%); in filter-feeding bighead carp (40.68% ± 17.63%) and silver carp (7.64% ± 4.90%). In grass carp, a large proportion of OTUs belonging to the genus *Bacteroides* (19.11% ± 5.96%), followed by the genera *Aeromonas* (5.13% ± 3.72%), *Clostridium* XVIII (1.92% ± 1.18%), *Clostridium* XIVb (0.80% ± 0.63%),

<table><caption>TABLE 1 | Number of sequences analyzed, observed diversity richness (OTUs), and estimated sample coverage for 16S rRNA gene libraries of the different samples.</caption><thead><tr><th>Samples</th><th>No. sequences</th><th>OTUs</th><th>Coverage (%)</th></tr></thead><tbody><tr><td>FB1</td><td>5935</td><td>835</td><td>93.36</td></tr><tr><td>FB2</td><td>5935</td><td>427</td><td>96.41</td></tr><tr><td>FB3</td><td>5935</td><td>761</td><td>94.00</td></tr><tr><td>FS1</td><td>5935</td><td>913</td><td>93.55</td></tr><tr><td>FS2</td><td>5935</td><td>793</td><td>93.68</td></tr><tr><td>FS3</td><td>5935</td><td>790</td><td>93.29</td></tr><tr><td>HB1</td><td>5935</td><td>113</td><td>99.21</td></tr><tr><td>HB2</td><td>5935</td><td>159</td><td>98.57</td></tr><tr><td>HB3</td><td>5935</td><td>201</td><td>98.23</td></tr><tr><td>HG1</td><td>5935</td><td>143</td><td>98.75</td></tr><tr><td>HG2</td><td>5935</td><td>153</td><td>98.92</td></tr><tr><td>HG3</td><td>5935</td><td>124</td><td>99.26</td></tr><tr><td>OC1</td><td>5935</td><td>90</td><td>99.22</td></tr><tr><td>OC2</td><td>5935</td><td>60</td><td>99.63</td></tr><tr><td>OC3</td><td>5935</td><td>55</td><td>99.51</td></tr></tbody></table>

FB1-FB3 for bighead carp, FS1-FS3 for silver carp, HB1-HB3 for blunt snout bream, HG1-HG3 for grass carp, and OC1-OC3 for crucian carp.
<figure><img src="image_2.png" /></figure>

**FIGURE 1** | Differences in the composition of gut microbial and metabolite profiles. Non-metric Multidimensional Scaling (NMDS) plot showing variation in the composition (Bray–Curtis distance) of gut bacterial communities (A) and metabolite profiles (B) among different fish species. Hierarchical clustering of sample groups based on the mean relative abundance of each OTU (C) and metabolite (D) in each fish species. Abbreviations: FS, silver carp; FB, bighead carp; HG, grass carp; HB, blunt snout bream; OC, crucian carp.

<figure><img src="image_3.png" /></figure>

**FIGURE 2** | Bacterial composition of the different communities (relative read abundance of bacterial phyla within each community). Lanes FB1-FB3, FS1-FS3, HB1-HB3, HG1-HG3, and OC1-OC3 correspond to the three individuals of bighead carp, silver carp, blunt snout bream, grass carp, and crucian carp, respectively.
<figure><img src="image_4.png" /></figure>

**FIGURE 3** | Summary of the variation in gut microbiota, gut metabolite, and evolutionary distance among different fish species. (A) A heat map of the mean abundances of the prominent OTUs (average abundance > 0.1%) assigned to genus level among different fish species. The rows represent the 63 predominant bacterial OTUs, and the values in the heatmap represent the Z-transformed relative percentage of each OTU. Phylum and genus level classifications of OTUs (OTU ID in parentheses) are noted also. (B) A heat map of the mean percent abundances of top 50 metabolite among different fish species. Both heat map based on the same hierarchical clustering solution (Bray–Curtis distance metric and complete clustering method). (C) Phylogenetic tree of the five cyprinid fish species based on the concatenated sequences of the two mitochondrial genes. The numbers on the nodes are neighbor-joining bootstrap values (values > 50 are shown), and bootstrap was replicated 999 times. The bar shows the relative branch length distance computed using the Kimura 2-parameter method. HG, grass carp; HB, blunt snout bream; OC, crucian carp; FS, silver carp; and FB, bighead carp.

and *Proteocatella* (0.71% ± 0.39%). OTUs belonging to *Aeromonas* (2.86% ± 2.06%), *Clostridium* XI (0.34% ± 0.25%), *Proteocatella* (0.28% ± 0.38%), and *Pirellula* (0.23% ± 0.07%) were also common in blunt snout bream. Other genera were relatively abundant in the gut microbiota of silver carp and bighead carp, namely: *Pirellula*, *Gp6*, *Anaerorhabdus*, *Steroidobacter*, *Proteocatella*, *Dechloromonas*, *GpIIa*, *Clostridium* XI and *Clostridium sensu stricto* (cluster I). In the gut content of the two filter-feeding species, *Acidobacteria* clade *Gp6* and *Cyanobacteria* clade *GpIIa* were also abundant (**Figure 3A**).

The prominent OTUs (mentioned in **Figure 3A**) used as inputs in the Random-Forest model paired with Boruta feature selection. An importance score, the mean decrease accuracy (MDA) was assigned to each OTU, based on the error increase corresponding to the removal of this feature from the predictors. A total of 24 significant OTUs were selected based on Boruta algorithm. The importance score (MDA) of these OTUs was illustrated in **Figure 4A**. Interestingly, the relative abundance

of OTUs belonging to *Pirellula*, *Bacteroides*, *Clostridium* (cluster I, XI, and XVIII), and *Cetobacterium* were among the most discriminative features between fish species.

In this investigation, 99 different metabolites were detected in the gut content samples by means of GC/MS analysis. The metabolite profiles of gut contents varied also among different fish species, mainly due to lactic acid and some amino acids. The gut samples of cyprinid fishes, especially crucian carp and silver carp, had a high percentage of lactic acid (**Figure 3B**). The gut samples of crucian carp had higher percentages of phosphoric acid, whereas the gut samples of filter-feeding carps had higher percentages of leucine, followed by herbivorous carps. Other dominant metabolites corresponded to alanine, isoleucine, and valine, which were more abundant in the filter-feeding species.

The phylogenetic structure inferred here was consistent with that of previous studies (Figure 3C) (Kong et al., 2007; Tao et al., 2010; Li et al., 2012). The results showed that silver carp and bighead carp were the most closely related species
<figure><img src="image_5.png" /><figcaption>**FIGURE 4** | The correlation between some specific microbes and metabolite. (A) Top 24 OTUs that differentiate among different fish species as revealed by Random Forest. The phylum and genus classifications are also provided. The OTU ID is given after the genus. (B) Redundancy analysis (RDA) of some specific gut microbes responding to metabolite. Only those metabolite with relative abundances > 3% were shown.</figcaption></figure>

(Figure 3C), and the two filter-feeding species shared more diverse gut bacteria, such as OTUs of *Proteocatella*, *Clostridium* XI and *Clostridium sensu stricto*, and metabolites, such as leucine, alanine, isoleucine, and valine (Figures 3A,B). Among herbivorous and omnivorous carps, more closely related species tended to harbor more similar gut microbiota composition and metabolite profiles (Figures 3A–C). OTUs belonging to *Aeromonas*, *Clostridium* XI, *Proteocatella*, *Pirellula* and relatively high percentages of leucine were common in both grass carp and blunt snout bream, while, *Cetobacterium*, lactic acid and phosphoric acid were more abundant in crucian carp.

## The Relationship between Evolutionary Distance, Gut Microbiota and Metabolite Profiles

Mantel tests were performed to detect the correlations between the fish evolutionary distance, gut microbiota and metabolite profiles. In particular, we found that the evolutionary distance was positively associated with the dissimilarity of gut microbiota (Mantel test, *r* = 0.34, *p* = 0.009) and that of metabolite profiles (Mantel test, *r* = 0.36, *p* = 0.003). These correlations were

relatively low, but significant. The gut bacterial and metabolite profiles were also significantly related (Mantel test, *r* = 0.61, *p* = 0.001).

The correlation between some specific gut microbes and metabolites was detected using RDA (Figure 4B). For example, the OTUs assigned to *Pirellula* were positively correlated with isoleucine, valine, leucine, and alanine, while OTUs belonging to *Clostridium* (cluster I and XI) showed positive correlations with lactic acid, glucose, and phosphoric acid. In addition, OTUs of *Bacteroides* and *Clostridium* XVIII correlated positively with proline, serine, and threonine. *Cetobacterium* OTU1079 and OTU766 correlated positively with methionine and N,N-Di-(2-Hydroxyethyl)-methanamine, while OTU2177 in the same genus was positively correlated with phenylalanine.

The presumptive functions of the gut microbial communities were also examined using PICRUSt. The fish gut samples had NSTI values of 0.11 ± 0.05. For comparison, Langille et al. (2013) found that Human Microbiome Project samples had the lowest (best) NSTI values (0.03 ± 0.2). Other mammalian guts had a higher mean NSTI value (0.14 ± 0.06), and diverse communities such as soil also had a much higher NSTI value (0.17 ± 0.02). Thus, the fish gut samples provide a reasonable data set to examine predictions from PICRUSt. We found there were
significant differences in predicted microbial functions across species.

A total of 29 significant level 3 KEGG Orthology (KO) functions were selected based on Boruta algorithm. The majority of the most discriminative putative functions were in the category of metabolism (Supplementary Figure S4A). In this context, the most noteworthy functional genes were involved in carbohydrate metabolism (butanoate metabolism, citrate cycle, glycolysis/gluconeogenesis, and pyruvate metabolism), amino acid metabolism (arginine and proline metabolism, tyrosine metabolism and selenocompound metabolism), energy metabolism (nitrogen metabolism), lipid metabolism (glycerolipid metabolism), metabolism of cofactors and vitamins (porphyrin and chlorophyll metabolism), and xenobiotics biodegradation and metabolism (chloroalkane and chloroalkene degradation). The result indicated that these categories were particularly important in differentiating the putative functions of gut microbiota among fish species. The distribution of points representing the inferred function of gut microbiota suggested that fishes within the same feeding habits had more similarity in the predicted functions of the microbiome than fishes with different feeding habits (Supplementary Figure S4B).

## DISCUSSION

To our knowledge, the present study was the first one that addressed the relationship between host phylogeny, gut microbiota and metabolic profiles in cyprinid fishes raised under identical husbandry conditions. Our results showed (1) the presence of specific fermentative bacteria populations in different cyprinid species; (2) the correlations between host evolutionary distance, gut microbiota and metabolic profiles of cyprinid fishes; (3) the influence of feeding habits and fish species on gut microbial and metabolic profiles. Li et al. (2014) explored the gastrointestinal microbiota of eight fish species with different feeding habits, but these wild fishes were not raised under identical husbandry conditions. In particular, the differences of gut microbiota composition between closely related fish species were not clearly defined, and the specificity as an adaptive measure remains undetermined. Our results shed new light on the microbiota of cyprinid fishes and highlighted the correlations between host genotype, gut microbiota and presumptive functions, and gut metabolite profile. Furthermore, understanding the microbial community and metabolite profile in the guts of these fish species can provide useful information to improve the health and productivity of these commercially valuable freshwater species.

## The Composition of Gut Microbial and Metabolite Profiles Differed

The filter-feeding carps, which have the smallest phylogenetic distance between them, displayed higher alpha-diversity and greater intragroup variability, compared with the other dietary groups (Supplementary Figures S1, S3). An explanation may be ascribed to feeding behavior, although feed pellets were supplied in the pond. In general, diet-associated microbes have a wider

diversity when they originate from diverse feeds, and the animals that consume more diverse feeds may be exposed to carry more diverse microbes (Laparra and Sanz, 2010). Compared with other species, the filter-feeding carps actively swim all around and expand a variety of feed items by filtering large volumes of water (Fang et al., 2013; Zhang and Kang, 2013). Such behavior can increase the variability in the composition of feed intake. In fact, it has been shown that diet diversity was associated with the variation of gut microbiota composition (Li H. et al., 2016). This filter-feeding activity resulted likely in more diverse microbes and greater variation of microbiota composition in filter-feeding carps (Verschuere et al., 2000; Tetlock et al., 2012; Li H. et al., 2016).

The environment affects the gut microbiota of fish and mammals (Ley et al., 2008a; Wu et al., 2012), but the same environment in the present study did not result in similar gut microbiota and metabolite profiles among the five carp species. Our results demonstrated that host species could affect not only gut microbiota, but also metabolite profiles (Figure 1). Recent studies showed that some particular genes may influence host immunity or physiology, and control the size of individual microbial populations (Benson et al., 2010; Bolnick et al., 2014). However, diet may also influence the phylogenetic relationship of gut microbiota among hosts. In our study, the natural feeding habits among the five carp species differed, and their feed intake was likely different, even though commercial food was supplied in the pond. As most animals lack the ability to degrade and to digest cellulose, however, certain species are capable of digesting cellulose because of their gut microbiota (Poulsen et al., 2014; Wang et al., 2015). Plant polysaccharide-degrading bacteria are particularly important for food degradation in the gut of herbivorous carps, namely grass carp and blunt snout bream (Wu et al., 2012, 2015). Contrastingly, filter-feeding activity resulted in a much higher bacterial diversity in the gut of silver carp and bighead carp, compared to those of the other species (Verschuere et al., 2000; Tetlock et al., 2012). Thus, the diet may also explain the differences of phylogenetic relationship in host-gut microbiota.

As the feeding habits similarity was related to the phylogenetic distance between fishes (Figure 3C), it was not possible to differentiate the relative influence of the diet from that of genetic factors, when studying the gut microbial composition and metabolite profiles in this study. Future interspecific studies with well-defined diets are needed for addressing the issue.

## Correlation between Host Genetics, Gut Microbiota and Metabolic Profiles of Cyprinid Fishes

Additionally, the carps with the same feeding habits had more similarity in gut microbiota composition and metabolite profiles (Figures 3A,B). This similarity among closely related host species suggested some degree of correlations between hosts and microbes, and further analysis showed the positive correlations between evolutionary distance and gut microbiota or metabolite profiles. In vertebrates, closely related host lineages harbor more similar gut microbiota than distantly related
lineages (Ochman et al., 2010; Sanders et al., 2014). *Bacteroides*, *Aeromonas*, *Clostridium* XVIII, *Clostridium* XIVb, *Clostridium* XI, and *Pirellula* were the dominant genera identified in the herbivorous carps (**Figure 3A**). Some strains of these genera are associated with plant-rich diets, and widely known as cellulose-degrading bacteria, which are particularly important for food degradation in the gut of herbivorous carps (Hayashi et al., 2002; Flint et al., 2008; King et al., 2012; Wu et al., 2012; Ziemer, 2013; Li et al., 2015). The two filter-feeding species shared more diverse dominant gut bacteria (**Figure 3A**). Some members of these bacterial communities included not only species with saccharolytic and fiber-fermenting activities, but also proteolytic species, such as *Proteocatella*, *Clostridium* XI and *Clostridium sensu stricto* (Lubbs et al., 2009; Pikuta et al., 2009; Schwab et al., 2011). Our results are consistent with another convincing study, suggesting that differences in the gut microbiota of closely related stickleback populations were correlated with host genetic divergence, even after controlling for food and water microbes (Smith et al., 2015). This host genetic control may explain the conserved composition of the gut microbiota with closely related species. By exerting top-town selection pressure, host genetic control may overcome microbial competition within the gut ecosystem and promote microbes that benefit the host (Benson et al., 2010). However, a correlation with host phylogeny by itself does not allow to discern whether genetic factors are acting through the direct physiological control of gut microbiota by the host (Benson et al., 2010; Bolnick et al., 2014; Thaiss et al., 2014), or whether the correlation was due to the dietary strategy of the species (Ley et al., 2008a; Muegge et al., 2011).

The metabolic importance of gut microbiota was illustrated by the fact that genetically homogeneous fishes may have diverse metabolic profiles when they have structurally different gut microbiota. The gut bacteria in carps belong to clades of fermentative bacteria, whose metabolism may yield varying amounts of products (Ni et al., 2014; Wu et al., 2015). These differences suggest that variation in the structure of gut bacterial communities should result in variation in the structure of gut metabolic profiles. Here we found that the gut bacterial and metabolic profiles of cyprinid fishes was significantly related (Mantel test, *r* = 0.61, *p* = 0.001; **Figures 3A,B**). The discriminative OTUs between different fish species (**Figure 4B**) included members of *Pirellula*, *Bacteroides*, *Clostridium* (cluster I, XI, and XVIII), and *Cetobacterium*. Most of these genera are involved in food degradation, as described above (Lubbs et al., 2009; Schwab et al., 2011; King et al., 2012; Ziemer, 2013). A difference in abundance of metabolites suggested a difference in the ability of food utilization by gut bacteria among cyprinid fishes. Due to the complex composition of gut microbiota, various metabolites were detected which correlated positively with various bacterial taxa (**Figure 4B**). Our results showed that fermentative subsets of gut bacterial taxa might be selectively stimulated by food (carbohydrates, proteins or lipids) and might therefore contribute to digestion, especially by degrading plant-derived polysaccharides in fish gut. However, this study disclosed the correlations between metabolism and gut microbiota, but not the causal relationships. In order to investigate the functional roles of gut microbiota, especially those

associated with food degradation and digestive metabolism, a strict control of experimental conditions in laboratories should be applied, and multiple analytical approaches are necessary.

Lactic acid and SCFAs are among the dominant metabolites produced by commensal microbiota. Aldunate et al. (2015) suggested their potential use as biomarkers of disease and/or disease susceptibility, due to antimicrobial and immune modulatory properties. Interestingly, amino acids can serve as precursors for the synthesis of SCFA by bacteria, suggesting an interplay between microbial activity and host amino acid and SCFA homeostasis (Neis et al., 2015). The PICRUSt results suggested also that gut microbiota may exert physiological functions linked to amino acid metabolism and glycolysis (Supplementary Figure S4), which is consistent with another study about the gut microbiome of grass carp (Ni et al., 2014). However, these results still need further confirmation.

In summary, this is the first report that addressed the relationship between host phylogeny, gut microbiota and metabolic profiles in cyprinid fishes using next-generation sequencing and GC/MS. The results showed that host genotype, gut microbiota, and gut metabolite profile have concordant phylogenic relationship. However, only five fish species with a small sample size of three individuals per species were used to construct the relationships between host phylogeny, gut microbiota and metabolic profiles. The correlation scores presented in our study are between 0.3 and 0.4. This relatively low correlation coefficient may be due to insufficient sampling. More diverse host species and larger populations should be surveyed in the future work to confirm these results and extend the knowledge about the relationship among host evolutionary distance, gut microbiota and metabolic profiles in cyprinid fishes. In addition, it remains to be determined whether evolutionary distance is associated with gut microbial diversity or metabolic profiles in other host species.

## AUTHOR CONTRIBUTIONS

TL, ML, DF, and AL conceived the research. TL, ML, HL, and XZ performed the experiments. TL wrote the manuscript. TL, HL, F-JG, and AL edited the manuscript. QZ, XZ, and ML contributed sampling, reagents or data analysis pipeline. All authors reviewed and accepted the manuscript.

## ACKNOWLEDGMENTS

This work was supported by grants from FEBL project (2016FBZ04), Jiangsu Marine and Fisheries Bureau project (D2015-11), and the Natural Science Foundation of China (No. 30670112 and 31070112).

## SUPPLEMENTARY MATERIAL

The Supplementary Material for this article can be found online at: http://journal.frontiersin.org/article/10.3389/fmicb.2017.00454/full#supplementary-material
REFERENCES

Aldunate, M., Srinovski, D., Hearps, A. C., Latham, C. F., Ramsland, P. A., Gugasyan, R., et al. (2015). Antimicrobial and immune modulatory effects of lactic acid and short chain fatty acids produced by vaginal microbiota associated with eubiosis and bacterial vaginosis. *Front. Physiol.* 6:164. doi: 10.3389/fphys. 2015.00164

Amato, K. R. (2013). Co-evolution in context: the importance of studying gut microbiomes in wild animals. *Microbiome Sci. Med.* 1, 10–29. doi: 10.2478/ micsm-2013-0002

Antharam, V. C., McEwen, D. C., Garrett, T. J., Dossey, A. T., Li, E. C., Kozlov, A. N., et al. (2016). An integrated metabolomic and microbiome analysis identified specific gut microbiota associated with fecal cholesterol and coprostanol in clostridium difficile infection. *PLoS ONE* 11:e0148824. doi: 10.1371/journal.pone.0148824

Benson, A. K., Kelly, S. A., Legge, R., Ma, F., Low, S. J., Kim, J., et al. (2010). Individuality in gut microbiota composition is a complex polygenic trait shaped by multiple environmental and host genetic factors. *Proc. Natl. Acad. Sci. U.S.A.* 107, 18933–18938. doi: 10.1073/pnas.1007028107

Bolnick, D. I., Snowberg, L. K., Caporaso, J. G., Lauber, C., Knight, R., and Stutz, W. E. (2014). Major histocompatibility complex class IIb polymorphism influences gut microbiota composition and diversity. *Mol. Ecol.* 23, 4831–4845. doi: 10.1111/mec.12846

Breiman, L. (2001). Random forests. *Mach. Learn.* 45, 5–32. doi: 10.1023/A: 1010933404324

Chen, W., Liu, F., Ling, Z., Tong, X., and Xiang, C. (2012). Human intestinal lumen and mucosa-associated microbiota in patients with colorectal cancer. *PLoS ONE* 7:e39743. doi: 10.1371/journal.pone.0039743

DeSantis, T. Z., Hugenholtz, P., Larsen, N., Rojas, M., Brodie, E. L., Keller, K., et al. (2006). Greengenes, a chimera-checked 16S rRNA gene database and workbench compatible with ARB. *Appl. Environ. Microbiol.* 72, 5069–5072. doi: 10.1128/aem.03006-05

Faber, F., Tran, L., Byndloss, M. X., Lopez, C. A., Velazquez, E. M., Kerrinnes, T., et al. (2016). Host-mediated sugar oxidation promotes post-antibiotic pathogen expansion. *Nature* 534, 697–699. doi: 10.1038/nature18597

Fang, M., Cai, L., Gao, Y., Tu, Z.-Y., Wang, C.-F., and Huang, Y.-P. (2013). Effect of temperature on swimming capability and oxygen consumption of juvenile *Hypophthalmichthys molitrix*. *J. Hydroecol.* 34, 49–53. doi: 10.1007/s00360-012- 0690-7

Flint, H. J., Bayer, E. A., Rincon, M. T., Lamed, R., and White, B. A. (2008). Polysaccharide utilization by gut bacteria: potential for new insights from genomic analysis. *Nat. Rev. Microbiol.* 6, 121–131. doi: 10.1038/nrmicro1817

Garner, C. E., Smith, S., de Lacy Costello, B., White, P., Spencer, R., Probert, C. S., et al. (2007). Volatile organic compounds from feces and their potential for diagnosis of gastrointestinal disease. *FASEB J.* 21, 1675–1688. doi: 10.1096/fj. 06-6927.com

Gentleman, R. C., Carey, V. J., Bates, D. M., Bolstad, B., Dettling, M., Dudoit, S., et al. (2004). Bioconductor: open software development for computational biology and bioinformatics. *Genome Biol.* 5:R80. doi: 10.1186/gb-2004-5- 10-r80

Gill, S. R., Pop, M., Deboy, R. T., Eckburg, P. B., Turnbaugh, P. J., Samuel, B. S., et al. (2006). Metagenomic analysis of the human distal gut microbiome. *Science* 312, 1355–1359. doi: 10.1126/science.1124234

Harrigan, G. G., and Goodacre, R. (2012). *Metabolic Profiling: Its Role in Biomarker Discovery and Gene Function Analysis*. Berlin: Springer Science & Business Media.

Hayashi, H., Sakamoto, M., and Benno, Y. (2002). Fecal microbial diversity in a strict vegetarian as determined by molecular analysis and cultivation. *Microbiol. Immunol.* 46, 819–831. doi: 10.1111/j.1348-0421.2002.tb02769.x

Kanehisa, M., Goto, S., Sato, Y., Furumichi, M., and Tanabe, M. (2012). KEGG for integration and interpretation of large-scale molecular data sets. *Nucleic Acids Res.* 40, D109–D114. doi: 10.1093/nar/gkr988

Kimura, M. (1980). A simple method for estimating evolutionary rates of base substitutions through comparative studies of nucleotide sequences. *J. Mol. Evol.* 16, 111–120. doi: 10.1007/BF01731581

King, G. M., Judd, C., Kuske, C. R., and Smith, C. (2012). Analysis of stomach and gut microbiomes of the eastern oyster (*Crassostrea virginica*) from coastal Louisiana, USA. *PLoS ONE* 7:e51475. doi: 10.1371/journal.pone.0051475

Kong, X., Wang, X., Gan, X., Li, J., and He, S. (2007). Phylogenetic relationships of Cyprinidae (Teleostei: Cypriniformes) inferred from the partial S6K1 gene sequences and implication of indel sites in intron 1. *Sci. China C Life Sci.* 50, 780–788. doi: 10.1007/s11427-007-0076-3

Kursa, M. B., and Rudnicki, W. R. (2010). Feature selection with the boruta package. *J. Stat. Softw.* 36, 1–13. doi: 10.18637/jss.v036.i11

Langille, M. G., Zaneveld, J., Caporaso, J. G., McDonald, D., Knights, D., Reyes, J. A., et al. (2013). Predictive functional profiling of microbial communities using 16S rRNA marker gene sequences. *Nat. Biotechnol.* 31, 814–821. doi: 10.1038/nbt.2676

Laparra, J. M., and Sanz, Y. (2010). Interactions of gut microbiota with functional food components and nutraceuticals. *Pharmacol. Res.* 61, 219–225. doi: 10.1016/j.phrs.2009.11.001

Ley, R. E., Hamady, M., Lozupone, C., Turnbaugh, P. J., Ramey, R. R., Bircher, J. S., et al. (2008a). Evolution of mammals and their gut microbes. *Science* 320, 1647–1651. doi: 10.1126/science.1155725

Ley, R. E., Lozupone, C. A., Hamady, M., Knight, R., and Gordon, J. I. (2008b). Worlds within worlds: evolution of the vertebrate gut microbiota. *Nat. Rev. Microbiol.* 6, 776–788. doi: 10.1038/nrmicro1978

Li, H., Li, T., Beasley, D. E., Hedenec, P., Xiao, Z., Zhang, S., et al. (2016). Diet diversity is associated with beta but not alpha diversity of pika gut microbiota. *Front. Microbiol.* 7:1169. doi: 10.3389/fmicb.2016.01169

Li, J., Ni, J., Li, J., Wang, C., Li, X., Wu, S., et al. (2014). Comparative study on gastrointestinal microbiota of eight fish species with different feeding habits. *J. Appl. Microbiol.* 117, 1750–1760. doi: 10.1111/jam.12663

Li, M., Wang, B., Zhang, M., Rantalainen, M., Wang, S., Zhou, H., et al. (2008). Symbiotic gut microbes modulate human metabolic phenotypes. *Proc. Natl. Acad. Sci. U.S.A.* 105, 2117–2122. doi: 10.1073/pnas.0712038105

Li, T., Long, M., Gatesoupe, F. J., Zhang, Q., Li, A., and Gong, X. (2015). Comparative analysis of the intestinal bacterial communities in different species of carp by pyrosequencing. *Microb. Ecol.* 69, 25–36. doi: 10.1007/s00248-014- 0480-8

Li, T., Long, M., Ji, C., Shen, Z., Gatesoupe, F. J., Zhang, X., et al. (2016). Alterations of the gut microbiome of largemouth bronze gudgeon (*Coreius guichenoti*) suffering from furunculosis. *Sci. Rep.* 6:30606. doi: 10.1038/srep30606

Li, X., Yan, Q., Xie, S., Hu, W., Yu, Y., and Hu, Z. (2013). Gut microbiota contributes to the growth of fast-growing transgenic common carp (*Cyprinus carpio* L.). *PLoS ONE* 8:e64577. doi: 10.1371/journal.pone.0064577

Li, X., Yu, Y., Feng, W., Yan, Q., and Gong, Y. (2012). Host species as a strong determinant of the intestinal microbiota of fish larvae. *J. Microbiol.* 50, 29–37. doi: 10.1007/s12275-012-1340-1

Liscic, J., Schauer, N., Kopka, J., Willmitzer, L., and Fernie, A. R. (2006). Gas chromatography mass spectrometry-based metabolite profiling in plants. *Nat. Protoc.* 1, 387–396. doi: 10.1038/nprot.2006.59

Lubbs, D., Vester, B., Fastinger, N., and Swanson, K. (2009). Dietary protein concentration affects intestinal microbiota of adult cats: a study using DGGE and qPCR to evaluate differences in microbial populations in the feline gastrointestinal tract. *J. Anim. Physiol. Anim. Nutr.* 93, 113–121. doi: 10.1111/j. 1439-0396.2007.00788.x

Martin, F. P. J., Wang, Y., Sprenger, N., Yap, I. K. S., Lundstedt, T., Lek, P., et al. (2008). Probiotic modulation of symbiotic gut microbial–host metabolic interactions in a humanized microbiome mouse model. *Mol. Syst. Biol.* 4, 157. doi: 10.1038/msb4100190

MoA, B. O. F. (2012). *China Fishery Statistics Yearbook 2012*. Beijing: China Agriculture Press.

Moeller, A. H., Li, Y., Ngole, E. M., Ahuka-Mundeke, S., Lonsdorf, E. V., Pusey, A. E., et al. (2014). Rapid changes in the gut microbiome during human evolution. *Proc. Natl. Acad. Sci. U.S.A.* 111, 16431–16435. doi: 10.1073/pnas. 1419136111

Muegge, B. D., Kuczynski, J., Knights, D., Clemente, J. C., González, A., Fontana, L., et al. (2011). Diet drives convergence in gut microbiome functions across mammalian phylogeny and within humans. *Science* 332, 970–974. doi: 10.1126/ science.1198719

Nayak, S. K. (2010). Role of gastrointestinal microbiota in fish. *Aquac. Res.* 41, 1553–1573. doi: 10.1111/j.1365-2109.2010.02546.x

Neis, E. P., Dejong, C. H., and Rensen, S. S. (2015). The role of microbial amino acid metabolism in host metabolism. *Nutrients* 7, 2930–2946. doi: 10.3390/ nu7042930
Ni, J., Yan, Q., Yu, Y., and Zhang, T. (2014). Factors influencing the grass carp gut microbiome and its effect on metabolism. *FEMS Microbiol. Ecol.* 87, 704–714. doi: 10.1111/1574-6941.12256

Nicholson, J. K., Lindon, J. C., and Holmes, E. (1999). 'Metabonomics': understanding the metabolic responses of living systems to pathophysiological stimuli via multivariate statistical analysis of biological NMR spectroscopic data. *Xenobiotica* 29, 1181–1189. doi: 10.1080/004982599238047

Nicholson, J. K., and Wilson, I. D. (2003). Opinion: understanding 'global' systems biology: metabonomics and the continuum of metabolism. *Nat. Rev. Drug Discov.* 2, 668–676. doi: 10.1038/nrd1157

Ochman, H., Worobey, M., Kuo, C.-H., Ndjango, J.-B. N., Peeters, M., Hahn, B. H., et al. (2010). Evolutionary relationships of wild hominids recapitulated by gut microbial communities. *PLoS Biol.* 8:e1000546. doi: 10.1371/journal.pbio.1000546

Pikuta, E. V., Hoover, R. B., Marsic, D., Whitman, W. B., Lupa, B., Tang, J., et al. (2009). *Proteocatella sphenisci* gen. nov., sp. nov., a psychrotolerant, spore-forming anaerobe isolated from penguin guano. *Int. J. Syst. Evol. Microbiol.* 59, 2302–2307. doi: 10.1099/ijs.0.002816-0

Poulsen, M., Hu, H., Li, C., Chen, Z., Xu, L., Otani, S., et al. (2014). Complementary symbiont contributions to plant decomposition in a fungus-farming termite. *Proc. Natl. Acad. Sci. U.S.A.* 111, 14500–14505. doi: 10.1073/pnas.131971811

Rezzi, S., Ramadan, Z., Fay, L. B., and Kochhar, S. (2007). Nutritional metabonomics: applications and perspectives. *J. Proteome Res.* 6, 513–525. doi: 10.1021/pr060522z

Roeselers, G., Mittge, E. K., Stephens, W. Z., Parichy, D. M., Cavanaugh, C. M., Guillemin, K., et al. (2011). Evidence for a core gut microbiota in the zebrafish. *ISME J.* 5, 1595–1608. doi: 10.1038/ismej.2011.38

Roessner, U., Wagner, C., Kopka, J., Trethewey, R. N., and Willmitzer, L. (2000). Technical advance: simultaneous analysis of metabolites in potato tuber by gas chromatography-mass spectrometry. *Plant J.* 23, 131–142. doi: 10.1046/j.1365-313x.2000.00774.x

Sanders, J. G., Powell, S., Kronauer, D. J., Vasconcelos, H. L., Frederickson, M. E., and Pierce, N. E. (2014). Stability and phylogenetic correlation in gut microbiota: lessons from ants and apes. *Mol. Ecol.* 23, 1268–1283. doi: 10.1111/mec.12611

Schloss, P. D., Westcott, S. L., Ryabin, T., Hall, J. R., Hartmann, M., Hollister, E. B., et al. (2009). Introducing mothur: open-source, platform-independent, community-supported software for describing and comparing microbial communities. *Appl. Environ. Microbiol.* 75, 7537–7541. doi: 10.1128/AEM.01541-09

Schwab, C., Cristescu, B., Northrup, J. M., Stenhouse, G. B., and Gänzle, M. (2011). Diet and environment shape fecal bacterial microbiota composition and enteric pathogen load of grizzly bears. *PLoS ONE* 6:e27905. doi: 10.1371/journal.pone.0027905

Singh, P., Teal, T. K., Marsh, T. L., Tiedje, J. M., Mosci, R., Jernigan, K., et al. (2015). Intestinal microbial communities associated with acute enteric infections and disease recovery. *Microbiome* 3, 45. doi: 10.1186/s40168-015-0109-2

Smith, C. C., Snowberg, L. K., Gregory Caporaso, J., Knight, R., and Bolnick, D. I. (2015). Dietary input of microbes and host genetic variation shape among-population differences in stickleback gut microbiota. *ISME J.* 9, 2515–2526. doi: 10.1038/ismej.2015.64

Sugita, H., Miyajima, C., and Deguchi, Y. (1991). The vitamin B12-producing ability of the intestinal microflora of freshwater fish. *Aquaculture* 92, 267–276. doi: 10.1016/0044-8486(91)90028-6

Sullam, K. E., Essinger, S. D., Lozupone, C. A., O’CONNOR, M. P., Rosen, G. L., Knight, R., et al. (2012). Environmental and ecological factors that shape the gut bacterial communities of fish: a meta-analysis. *Mol. Ecol.* 21, 3363–3378. doi: 10.1111/j.1365-294X.2012.05552.x

Tamura, K., Stecher, G., Peterson, D., Filipski, A., and Kumar, S. (2013). MEGA6: molecular evolutionary genetics analysis version 6.0. *Mol. Biol. Evol.* 30, 2725–2729. doi: 10.1093/molbev/mst197

Tao, W., Zou, M., Wang, X., Gan, X., Mayden, R. L., and He, S. (2010). Phylogenomic analysis resolves the formerly intractable adaptive diversification of the endemic clade of east Asian Cyprinidae (Cypriniformes). *PLoS ONE* 5:e13508. doi: 10.1371/journal.pone.0013508

Tetlock, A., Yost, C. K., Stavrinides, J., and Manzon, R. G. (2012). Changes in the gut microbiome of the sea lamprey during metamorphosis. *Appl. Environ. Microbiol.* 78, 7638–7644. doi: 10.1128/AEM.01640-12

Thaiss, C. A., Zeevi, D., Levy, M., Zilberman-Schapira, G., Suez, J., Tengeler, A. C., et al. (2014). Transkingdom control of microbiota diurnal oscillations promotes metabolic homeostasis. *Cell* 159, 514–529. doi: 10.1016/j.cell.2014.09.048

Thompson, J. D., Higgins, D. G., and Gibson, T. J. (1994). CLUSTAL W: improving the sensitivity of progressive multiple sequence alignment through sequence weighting, position-specific gap penalties and weight matrix choice. *Nucleic Acids Res.* 22, 4673–4680. doi: 10.1128/MMBR.64.4.655-671.2000

Vernocchi, P., Del Chierico, F., and Putignani, L. (2016). Gut microbiota profiling: metabolomics based approach to unravel compounds affecting human health. *Front. Microbiol.* 7:1144. doi: 10.3389/fmicb.2016.01144

Vernocchi, P., Vannini, L., Gottardi, D., Del Chierico, F., Serrazanetti, D. I., Ndagijimana, M., et al. (2012). Integration of datasets from different analytical techniques to assess the impact of nutrition on human metabolome. *Front. Cell. Infect. Microbiol.* 2:156. doi: 10.3389/fcimb.2012.00156

Verschuere, L., Rombaut, G., Sorgeloos, P., and Verstraete, W. (2000). Probiotic bacteria as biological control agents in aquaculture. *Microbiol. Mol. Biol. Rev.* 64, 655–671. doi: 10.1128/MMBR.64.4.655-671.2000

Vitali, B., Ndagijimana, M., Cruciani, F., Carnevali, P., Candela, M., Guerzoni, M. E., et al. (2010). Impact of a synbiotic food on the gut microbial ecology and metabolic profiles. *BMC Microbiol.* 10:4. doi: 10.1186/1471-2180-10-4

Walter, J., Britton, R. A., and Roos, S. (2011). Host-microbial symbiosis in the vertebrate gastrointestinal tract and the *Lactobacillus reuteri* paradigm. *Proc. Natl. Acad. Sci. U.S.A.* 108(Suppl. 1), 4645–4652. doi: 10.1073/pnas.100099107

Wang, Y., Lu, Y., Zhang, Y., Ning, Z., Li, Y., Zhao, Q., et al. (2015). The draft genome of the grass carp (*Ctenopharyngodon idellus*) provides insights into its evolution and vegetarian adaptation. *Nat. Genet.* 47, 625–631. doi: 10.1038/ng.3280

Wikoffa, W. R., Anforab, A. T., Liub, J., Schultzb, P. G., Lesleyb, S. A., Petersb, E. C., et al. (2009). Metabolomics analysis reveals large effects of gut microflora on mammalian blood metabolites. *Proc. Natl. Acad. Sci. U.S.A.* 106, 3698–3703. doi: 10.1073/pnas.0812874106

Wu, S., Ren, Y., Peng, C., Hao, Y., Xiong, F., Wang, G., et al. (2015). Metatranscriptomic discovery of plant biomass-degrading capacity from grass carp intestinal microbiomes. *FEMS Microbiol. Ecol.* 91:fiv107. doi: 10.1093/femsec/fiv107

Wu, S., Wang, G., Angert, E. R., Wang, W., Li, W., and Zou, H. (2012). Composition, diversity, and origin of the bacterial community in grass carp intestine. *PLoS ONE* 7:e03440. doi: 10.1371/journal.pone.0030440

Ye, L., Amberg, J., Chapman, D., Gaikowski, M., and Liu, W.-T. (2013). Fish gut microbiota analysis differentiates physiology and behavior of invasive Asian carp and indigenous American fish. *ISME J.* 8, 541–551. doi: 10.1038/ismej.2013.181

Zhang, Q., and Kang, B. (2013). Effect of water velocity on swimming behavior of juvenile *Megalobrama amblycephala*. *Zool. Res.* 34, 429–436. doi: 10.11813/j.issn.0254-5853.2013.4.0429

Zhou, Y., Liao, Q., Lin, M., Deng, X., Zhang, P., Yao, M., et al. (2014). Combination of (1)H NMR- and GC-MS-based metabonomics to study on the toxicity of *Coptidis Rhizome* in rats. *PLoS ONE* 9:e88281. doi: 10.1371/journal.pone.0088281

Ziemer, C. J. (2013). Newly cultured bacteria with broad diversity isolated from 8 week continuous culture enrichments of cow feces on complex polysaccharides. *Appl. Environ. Microbiol.* 80, 574–585. doi: 10.1128/AEM.03016-13

Zoetendal, E. G., Akkermans, A. D., Akkermans-van Vliet, W. M., de Visser, J. A. G., and de Vos, W. M. (2001). The host genotype affects the bacterial community in the human gastrinointestinal tract. *Microb. Ecol. Health Dis.* 13, 129–134. doi: 10.1080/08910600175046269

**Conflict of Interest Statement:** The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Copyright © 2017 Li, Long, Li, Gatesoupe, Zhang, Zhang, Feng and Li. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) or licensor are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.