Next Article in Journal
The Structural, Functional and Evolutionary Impact of Transposable Elements in Eukaryotes
Previous Article in Journal
Sex-Differences in Traumatic Brain Injury in the Absence of Tau in Drosophila
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Non-Specific Immunity Associated Gut Microbiome in Aristichthys nobilis under Different Rearing Strategies

1
Key Laboratory of Freshwater Fish Resources and Reproductive Development, Ministry of Education, Southwest University, Chongqing 400715, China
2
Chongqing Fisheries Technology Promotion Station, Chongqing 400020, China
3
Chengdu Life Baseline Technology, Chengdu 610041, China
4
Australian Centre for Blood Diseases, Central Clinical School, Monash University, Melbourne, VIC 3004, Australia
*
Author to whom correspondence should be addressed.
Genes 2021, 12(6), 916; https://doi.org/10.3390/genes12060916
Submission received: 16 May 2021 / Revised: 31 May 2021 / Accepted: 7 June 2021 / Published: 14 June 2021
(This article belongs to the Section Microbial Genetics and Genomics)

Abstract

:
To understand the intestinal microbial diversity and community structure of bighead carp (Aristichthys nobilis) under different feeding strategies, 39 fish from three groups (A: 9 fish, natural live food only; B: 15 fish, natural live food + fish formulated feeds; C: 15 fish, natural live food + fish formulated feed + lactic acid bacteria) were obtained for the high throughput 16S rRNA gene sequencing. We first examined five non-specific immunity indications of the carp—lysozyme (LZM), catalase (CAT), glutathione reductase (GR), glutathione peroxidase (GSH-PX), and superoxide dismutase (SOD). Interestingly, the composition of gut microbiota and related non-specific immune indices were affected by the feeding treatment of the bighead carp. Notably, all enzyme activity indexes were significantly different (p < 0.01) in the spleen and three enzyme activity indexes (LZM, GSH-PX, and SOD) had significant differences in the hepatopancreas (p < 0.001) of the carp from the three groups. The 16S rRNA gene sequencing showed higher diversity in groups B and C. Compared to group A, the relative abundance of Actinobacteria increased significantly and the relative abundance of Proteobacteria and Firmicutes decreased significantly in groups B and C at the phylum level. Functional analysis revealed the association between non-specific immune indicators and import genera in the hepatopancreas and spleen of bighead carp. This study provides new insights into the gut microbiomes and non-specific immune of bighead carp.

1. Introduction

The bighead carp (Aristichthys nobilis) is a cyprinid fish that is native to South China and mainly feeds on zooplankton and phytoplankton in natural waters [1,2]. As a well-known economic fish, the global annual production of bighead carp is about ~10 million tons per year. In recent years, with increasing demand for bighead carp consumption, there has been wide concern about high-density pond farming of bighead carp [3]. In addition, feeding bighead carp with various formulation feeds to increase yield has become popular and acceptable [4]. However, due to the large-scale use of fertilizers and various fish feeds, the fish breeding environment has gradually deteriorated, resulting in increased stress factors and decreased autoimmune capacity, which can cause chronic diseases of the fish, and is one of the most pressing concerns in modern fish farming [5].
Lactic acid bacteria are natural active microorganisms that have regulatory effects on the internal environment of the gastrointestinal tract of animals [6]. They were reported to have the ability to improve the ecological environment of intestinal microbes through biological antagonisms, strengthen the barrier layer of intestinal epithelial cells, and improve the immune function of animals [7]. Previous studies on lactic acid bacteria in feed research showed that lactic acid bacteria can improve the intestinal environment of livestock and poultry, inhibit the growth and reproduction of harmful bacteria, and enhance the disease resistance in pigs and chickens [8]. Some studies were conducted to investigate the association between lactic acid bacteria and fish farming. For example, Son et al. showed improved weight gain rate and feed efficiency of groupers after adding lactobacillus into their feeds for 4 weeks [9]. Aly et al. found that the weight and survival rate of tilapia greatly improved after lactic acid bacteria were applied to their feeds [10]. Suzer et al. applied lactobacillus to black snapper at the development stage of larva and found that the growth performance and digestive enzyme activity of larva improved, especially when lactobacillus was added to live feeds or water [11]. However, little is known about the effects of lactic acid bacteria in bighead carp.
Innate immunity usually plays a very important role in combating microbial infection for all animals [12], and in general, the lysozyme (LZM), catalase (CAT), glutathione reductase (GR), glutathione peroxidase (GSH-PX), and superoxide dismutase (SOD) are key indices representing the non-specific immunity of animals in vivo [13]. LZM is a widely distributed cationic enzyme, which can destroy and eliminate foreign matter invading the body [14]. CAT and SOD are the main antioxidant enzymes of the oxidation resistance defense system in fish [14]. GSH-PX is an enzyme that prevents damage to cellular membranes and participates in the immune response of animals, which can directly reflect the body’s antioxidant capacity and immune level [15]. Thus, these indicators can be used to show the immune ability of fish and increase of the immune enzyme activity can enhance the immune system of fish.
Intestinal microbiota, as an important part of the host body, also plays an important role in fish health [16]. The dominant microorganisms in the fish gut are usually bacteria [17], which are mainly composed of aerobic bacteria and facultative anaerobic bacteria, among them, three core microorganisms were dominant: Firmicutes, Proteobacteria, and Bacteroidetes [18,19]. However, our knowledge about the intestinal flora of filter-feeding fish fed with formulated feeds is still poor. Traditional microbial research methods mainly rely on microbial culture technology, which can only detect less than 1% of microorganisms in nature [20]. Compared to traditional microbial identification and analysis methods, in the last ten years, the 16S rRNA high-throughput sequencing technology, a new generation sequencing method that possesses high accuracy and sensitivity [21,22], was widely used in the intestinal microbial studies, can greatly improve the depth and breadth of the analysis of microbial diversity in samples, and can improve our understanding of the structure and function of intestinal microflora. Considering the influence of lactobacillus on the fish gut bacterial community [8,9,10,11], using 16S rRNA sequencing technology to explore how formulated feeds change the intestinal microbial composition of bighead carp would be of certain interest. The aims of the present study include the investigation of the intestinal microbiome differences under natural feeding and fodder feeding, specially focusing on the effect of lactic acid bacteria on the intestinal flora and non-specific immune enzyme activity of bighead fish. The results will improve our understanding of the gut microbiome of bighead carp and will provide new ideas to study the fish intestinal flora.

2. Materials and Methods

2.1. Experimental Designs and Fish Management

All experimental bighead carp species were cultivated from the same batch of larvae from the Aquatic Products Research and Development Center of Jingzhi (Chongqing) Smart Agriculture Company. The experiment was carried out in eight outer pond cages (2 m × 2 m × 1.5 m) with 2.5 m of average water depth. Three groups of bighead carp were fed by different materials, including group A (2 parallel)—the control group fed with the ingested plankton in the water; group B (3 parallel) mainly fed with the artificial diets; and group C (3 parallel) fed with artificial feed plus lactic acid bacteria. Thus, a total of eight cages were used for these three groups and they were transferred from Jingzhi Smart Agriculture Company to the Beibei District Fish Base of the Southwest University (Chongqing). The artificial feed for this research uses the special feed for bighead carp from Aolong Kexiong Feed Co., Ltd. Chongqing, China. The biochemical composition of the artificial feed consisted of moisture (8.96%), crude protein (31.82%), crude fat (6.86%), crude ash (9.30%), crude fiber (4.91%), calcium (0.86%), and total phosphorus (1.39%). The culture experiment lasted for 8 weeks (56 days). Fish from groups B and C were fed twice a day at 8:30 a.m. and 5:30 p.m. The pond was regularly disinfected, bottomed, and oxygenated, according to normal aquaculture production management, to ensure water quality was under optimal conditions.

2.2. Sample Collection and Pre-Processing

At the end of the trial, fish were carefully captured by the nets in order to avoid external physiological stimulation. Five bighead carp with similar weight were randomly selected from each cage and anesthetized with MS-222 (Syndel, Ferndale, WA, USA) on ice. A total of 39 fish (group A: 9; group B: 15; group C: 15) were dissected. After the intestines were cut and taken out, the outer wall of the intestinal tube was wiped with 75% ethanol and then rinsed with sterile saline. The intestinal contents were evenly collected by 5 mL EP tubes and immediately placed in liquid nitrogen for 6 h. All samples were transported to the laboratory and stored at −80 °C refrigerator in time until DNA extraction.
Fish hepatopancreas and spleen were separated in order to determine related enzyme activity. The hepatopancreas and spleen were blotted with a filter paper to absorb the surface water and placed in sterile 1.5 mL centrifuge tubes. After the samples were centrifuged at 10,000 rpm for 15 min, the supernatant was collected and then stored at −20 °C to determine the enzyme activity. The experiments were performed in accordance with the Regulations for the Administration of Affairs Concerning Experimental Animals of China and the protocols applied in the present study were approved by the Ethics Committee of Southwest University (approved IACUC NO: 20200720-1).

2.3. Non-Specific Immunity Related Enzyme Index Determination

Lysozyme (LZM), catalase (CAT), glutathione reductase (GR), glutathione peroxidase (GSH-PX), superoxide dismutase (SOD) activities were measured by using commercial kits according to the manufacturers’ instructions (Nanjing JianCheng Bioengineering Institute, Nanjing, China).

2.4. DNA Extraction, PCR, and Sequencing

For bacterial 16S ribosomal RNA gene sequencing, microbial DNA was extracted from all 39 fish gut contents using the Qiagen Stool DNA Kit (Norcross, GA, USA) following the manufacturer’s protocol. The V4–V5 region of the bacterial 16S ribosomal RNA gene was amplified using the primers 515F (5′-GTGCCAGCMGCCGCGG-3′) and 907R (5′-CCGTCAATTCMTTTRAGTTT-3′), where the barcode was an eight-base sequence unique to each sample. The PCR experiment was carried out as described [23]. All samples were carried out in accordance with standard experimental condition in triplicates. The PCR products were then mixed and run on a 2% agarose gel. After the targeting regions were cut and recovered using the AxyPrepDNA Gel Recovery Kit (AXYGEN, Union City, CA, USA), they were quantified by Qubit® 3.0 (Invitrogen, Waltham, MA, USA) and mixed equally with barcodes. Then, the mixture was used to prepare the DNA pair-end library following the Illumina library preparation protocol and sequenced on the Illumina MiSeq platform in Shanghai BIOZERON (Shanghai, China) with paired-end 250 strategy. All raw sequencing data generated by this study can be accessed from the NCBI short read archive (SRA) platform under the accession number PRJNA695396.

2.5. Statistical and Bioinformatics Analysis

Raw FASTQ files were demultiplexed using in-house Perl scripts, according to the barcode sequences of all samples, and cleaned with the following criteria: (i) the 250 bp reads were truncated at any site receiving an average quality score <20 over a 10 bp sliding window and truncated reads shorter than 50 bp were discarded; (ii) exact barcode match, reads containing ambiguous characters were removed. We removed the reads that were not assembled [24]. After data cleaning, we identified operation taxonomic units (OTUs) for each sample using UPARSE (V7.1) with 97% similarity [25]. UCHIME was applied to identify and remove chimeric sequences based on both reference database and de novo assembly analysis [26]. The phylogenetic affiliation analysis of each 16S rRNA gene sequences was introduced by RDP Classifier against the SILVA (SSU115) 16S rRNA database with a confidence threshold of 70% [27,28].
The rarefaction analysis based on mothur (v.1.21.1) was conducted to identify the diversity indices [29], including the Chao, Shannon, and Simpson diversity indices. The unweighted pair-group method based on Bray–Curtis dissimilarity was used to calculate the distance between two samples, subsequent β-diversity analysis, and visual statistical analysis were performed with the data. Principal co-ordinate analysis (PCoA) was performed using the community ecology package and R-forge, and the results were visualized using the Vegan 2.0 package. Linear discriminant analysis effect size (LEfSe) analysis used to identify differentially abundant bacteria taxa between groups (A vs. B and A vs. C) according to the relative abundance and the p-value less than 0.05 was considered as significant difference [30]. Kruskal–Wallis sum-rank test was applied to examine the changes and differences between groups, followed by the LDA analysis to determine the size effect of each distinctively abundant taxa [31]. The enzyme activity data were presented in mean ± SD. We used the ‘boxplert’ function in the R for data statistics. One-way ANOVA analysis was used to assess the differences in bighead carp enzyme activity among the different groups. Once the significant difference was found, Dunn’s Multiple Test was used to perform the multiple comparisons. For correlation analysis, the relationships between genera OUT abundance and enzymatic activity were investigated separately through Pearson correlation calculation to explore the potential metabolic functions; the six genus abundance data were log2 transformed before Pearson analysis to reduce the effect of extremely high values, and the p-value less than 0.05 was considered a significant pairwise relationship.

3. Result

3.1. Biochemical Parameters Related to Enzyme Activity

To understand the effect of lactic acid bacteria on the intestinal flora and non-specific immune enzyme activity of bighead carp, we bred a total of 400 fish with similar weight in eight cages for three groups—group A (two parallel, fed with the ingested plankton in the water); group B (three parallel, fed with the artificial diets); and group C (three parallel, fed with the artificial diet plus lactic acid bacteria). The fish were maintained in the pond for two months and the characteristics can be found in Table 1. It is notable that compared to group A, bighead carp from groups B and C had significant increase of weight and visceral body ratio after two-month breeding. Then, the intestines were cut and taken out for subsequence analyses.
We first examined the enzyme activities of bighead carp form the eight cages (Table S1) and statistical analysis showed that feeding regime can significantly affect the enzymatic activities of the fish visceral tissue. Three indicators of enzymatic activity (LZM, GSH-PX, and SOD) showed significant differences (p < 0.001) in hepatopancreas (Table 1). The LZM activity of bighead carp was significantly lower (p < 0.001) in groups B and C, compared to group A. The GSH-PX activity was significantly higher in groups A and B, compared to group C, indicating that lactic acid bacteria may affect the level of enzyme activity in hepatopancreas. In addition, we found that bighead carp of group B had the highest SOD activity, and that group A had the lowest SOD activity (p < 0.001). Overall, artificial aquaculture practices can significantly increase the SOD activity. No significant (p > 0.05) differences in CAT and GR activities were observed in these groups. Interestingly, we found these five indicators with significant difference in the spleen (Table 1). Contrary to the results of hepatopancreas, the activity of LZM was significantly higher (p < 0.01) in groups B and C, compared to group A in the spleen. Higher CAT activity was found in group B and lower GR activity was found in group A. In addition, the GSH-PX activity was significant higher (p < 0.001) in group A compared to B and C.

3.2. Gut Bacterial Diversity

After data cleaning of the sequencing data, a total of 2,360,939 high-quality bacterial reads were obtained for the 39 samples, with an average of 60,654 (min: 32,680, max: 74,751) reads per sample (Table S2). Then, the OTUs were taxonomically classified into different levels using the RDP classifier at a 97% confidence. Based on four different indicators (the numbers of OTU observed, Chao, Shannon, and Simpson indexes) we observed that bighead carp had the highest biodiversity (the highest OUTS, Chao, and Shannon, the lowest Simpson indexes) in group C among the three groups (Table 2). Notably, compared with group B all the four indicators changed in group C, indicating that lactic acid bacteria may affect the diversity of intestinal microbiome in bighead carp. Further, the Shannon and Simpson indices showed significant difference in group B vs. C and group A vs. C (Figure 1a,b). Interestingly, no significant difference of Shannon and Simpson indices was observed between groups A and B (p > 0.05). We next performed the principal co-ordinate analysis (PCoA) based on Bray–Curtis distance and clearly separated the bighead carp of the three groups (Figure 1c). These results suggest that the microbiota profiles of carp might be modified by the artificial diet on the OUT level, but not significant (Table 2). In addition, by analyzing the intra-group and inter-group distances of all samples, we found that group B had a larger intra-group distance compare to groups A and C, at the same time, there was a larger intergroup distance between group B and group C (Figure 1d).

3.3. Gut Bacterial Composition

We analyzed the bacteria community composition at three different taxonomic levels. First, we found that at phylum level phyla Firmicutes, Actinobacteria, Proteobacteria, and Fusobacteria were the dominant bacteria in all three groups (Table S3). Compared to group A, Firmicutes and Proteobacteria were significantly decreased but Actinobacteria and Fusobacteria were significantly increased in both groups B and C (Table S4). Notably, the minor bacteria phylum Planctomycetes showed higher expression (Table S4, p < 0.05) in group C. At the order level, we found the predominant orders of microorganism species of the bighead carp were Erysipelotrichales, Propionibacteriales, Clostridiales, Fusobacteriales, and Betaproteobacteriales (Figure 2b). Interestingly, the relative abundance of order Erysipelotrichales significantly decreased after the bighead carp were fed by the artificial diet while the Propionibacteriales and Pirellulales increased (Table S4). At the genus level (Figure 2c), we found that genus Dielma and Burkholderiaceae_uncultured dominated the intestinal microbiome in both groups A and B. Compared with group A, group B had one more dominant bacteria species Cetobacterium. However, the genus Dielma OUT was significant lower in group C, compared to groups A and B (Table S4). The dominant bacteria species in group C included Propionibacteriaceae_uncultured (phylum Bacteroidetes, order Propionibacteriales), Chloroplast norank (phylum Cyanobacteria, order Chloroplast) and Saccharimonadales_norank (phylum Bacteroidetes, order Saccharimonadales). In all three groups, we found that some sequences (approximately 30%) were assigned to unknown bacteria.

3.4. Differential Bacterial Composition

We next identified the significant differences in bacteria communities between the three different groups using LEfSe. Compared to group A, we identified 47 bacteria species differentially expressed (p < 0.05) in group B, including 21 downregulated and 26 upregulated (Table S4). The overexpressed bacteria in group B included Fusobacteria, Actinobacteria, Planctomycetes, Patescibacteria, and Cyanobacteria while the downregulated bacteria included Bacteroidetes and Proteobacteria at phylum level (Figure 3a). After the lactic acid bacteria was added to the feeds for bighead carp, we identified 64 bacteria species differentially expressed (p < 0.05) in group C (Figure 3b), including 24 downregulated (e.g., Proteobacteria, Firmicutes, Bacteroidetes) and 40 upregulated (e.g., Chloroflexi, Cyanobacteria, Actinobacteria, Planctomycetes, Patescibacteria, and Fusobacteria). Notably, the phylum Chloroflexi was the only bacteria significantly increased in group C compared to group B, which indicated that lactic acid bacteria may have increase the abundance of Chloroflexi in bighead carp gut flora. Compared with group B, Firmicutes was found to be decreased in group C.

3.5. Functional Analysis

In order to explore the potential non-specific immune function related gut microbes in bighead carp, based on the LEfSe results we manually screened six significant differential expressed genera in group C or group B compared to group A, including Cetobacterium, ZOR0006, Dielma, Acetobacteroides, Mycobacterium, and Pirellula. Pearson correlation analysis revealed that the genera Acetobacteroides, Cetobacterium, Dielma, Mycobacterium, and Pirellula were significantly correlated (p < 0.05) with the five enzyme activities in hepatopancreas and spleen, including 14 positive correlations and 13 negative correlations (Table S5). It is interesting that LZM was negatively correlated with the genus Pirellula in hepatopancreas (Figure 4a) but positively correlated with the genus Pirellula and Mycobacterium in the spleen (Figure 4b,c). GR was positively correlated with the genus Mycobacterium and Pirellula in the spleen (Figure 4d,e). In order to investigate the different effects of intestinal flora in hepatopancreas and spleen on the immune function of bighead carp, we compared all 14 positively correlated and 13 negatively correlated pairs. It is notable that LZM had opposite correlations with the enzyme activities of hepatopancreas and spleen (Table S5). Three pairs were found with consistent correlations in hepatopancreas and spleen, including Mycobacterium ~ GSH-PX (negative), Pirellula ~ GSH-PX (negative) and Cetobacterium ~ SOD (positive). The different non-specific immune functions of LZM in different tissues require further experiments to be explored.

4. Discussion

Fish intestinal microbes play an important role in inhibiting pathogenic microorganisms, the dominant flora in the intestine can protect the host from infection and damage of the complex microbes in the environment, the presence of certain intestinal microbial metabolites can stimulate the proliferation of intestinal epithelial cell and immune system response [32]. By analyzing the intestinal microbiota of healthy and diseased largemouth bronze gudgeon, Li et al. reported that the bacteria diversity was lower in diseased fish compared to healthy fish and that genus Aeromonas may be a key factor of causing infectious diseases of farmed fish [33]. Stagaman et al. analyzed the intestinal microorganisms of 68 wild type and 61 zebrafish lacking B or T cell receptor function, and showed that the immune could filter the composition of intestinal microorganisms, which might reflect the host immune status [34].
The mechanical, immune and biological barriers formed by the combination of animal intestinal flora and intestinal mucosa not only play an important role in maintaining the stability of the internal environment, but also effectively prevent the invasion of pathogenic substances and the displacement of bacterial endotoxin [35]. In this study, five indicators were considered to present the fish nonspecific immunity in the spleen and HP. Based on our analyses, LZM activity of bighead carp increased in group C compared to group B in the spleen after adding lactic acid bacteria, which is consistent with the findings by Son et al. in grouper blood [9], while in HP LZM, activity in groups B and C significantly decreased compared with group A (p < 0.001, Table 1). This intriguing finding might be consequences of the weakly acidic water after the addition of artificial feed or lactic acid bacteria and the HP is highly sensitive to water changes and may be suppressed to some extent. In addition, CAT and GR activities were changed in the spleen of bighead carp when after adding lactic acid bacteria alone (group C vs. group B), but there was no significant change in the HP, the results showed that lactic acid bacteria had different effects on CAT and GR activities in the spleen and HP of bighead carp. In a previous study, Lin et al. reported that lactic acid bacteria can significantly increase GSH-PX and SOD activities in the liver of mice compared with the control group [36]; however, very little research has been done in this area on fish. Our results showed that the GSH-PX and SOD activities decreased in HP but no change in the spleen of bighead carp when lactic acid bacteria were added to the feed, which can be explained by that there are certain differences between HP and spleen in the non-specific immune physiological functions of fish. Overall, although some enzyme activities were not observed with the same trends in HP and spleen of bighead carp when supplemented with lactic acid bacteria. However, according to the analysis of related enzyme activity index in the spleen, except CAT, all other enzyme activity indexes were increased or unchanged, these results roughly support that lactic acid bacteria are beneficial to improve the non-specific immune of bighead carp. Meanwhile, the effects and mechanisms of lactic acid bacteria in different fish organs may require further experiments to be explored.
It was reported that the gut specific microbial community of fish could be affected by the dietary manipulations [37,38,39]. However, there are few studies about the effects on gut microbiota of filter-feeding fish fed formulated feed. Earlier studies indicated the dominant phyla in the intestinal microbiota of healthy fish were Proteobacteria and Firmicutes [40,41], this is consistent with the overall results of our study. In addition, Actinobacteria was also found as a dominant phylum in the current study (Figure 2a), which is in accordance with that revealed for rainbow trout (Oncorhynchus mykiss) [42]. Actinobacteria, a phylum of Gram-positive bacteria, is an important secondary metabolite producer and plays a key role in animal intestines [43]. Actinobacteria have the ability to biosynthesize secondary metabolites as antibiotics against invasive pathogens [44]. Interestingly, we found that the relative abundance of Actinobacteria increased when the bighead carp were fed with formulated feed and/or lactic acid bacteria (Figure 2a). Likewise, Maria et al. found the phylum Actinobacteria was significantly increased when the colitis associated colon cancer model (C57BL/6 mice) was supplemented with probiotics [45]. While Ma et al. reported that the Actinobacteria as a dominant phylum in intestines of diseased fish compared to healthy controls [39], indicating that Actinobacteria may be useful biomarkers in immunocompromised transfused fish patients. Therefore, we believe that whether the phyla Actinobacteria is related to fish immunity needs to be further studied.
In this study, the genus Cetobacterium has the largest number of significant pair relationships in functional analysis (7 pairs, Table S5). It has been reported that genus Cetobacterium is widely distributed within the guts of freshwater fishes [46] and its prevalence is mainly contributing to the production of vitamin B12 in human [47]. Hence, Cetobacterium has been speculated to have a role in the synthesis of vitamin B12 in the fish gut [48]. However, in this study, vitamin B12 within the three groups was not investigated; further studies on this topic are needed. Because the abundance of Cetobacterium is positively correlated with the content of vitamin B12, supplementation of vitamin B12 in formulated fish feed may increase the abundance of Cetobacterium in group B and group C. In addition, as the dominant genera in LEfSe analyzed in group B and group C compared to group A, the Cetobacterium was revealed to be positively associated with SOD activity in HP and the spleen, and simultaneously suggested that the abundance of genus Cetobacterium may have a positive effect on the non-specific immune of fish. However, the genus Cetobacterium also showed a significant negative correlation with GSH-Px in the spleen; more microbial sequencing and functional activity study of intestinal microbiome in bighead carp are needed in the future.
In summary, the addition of lactic acid bacteria can significantly altered the composition of bighead carp intestinal flora, but the effects on non-specific immune indexes in HP and spleen were complex. It should be noted that this study did not carry out a variety of lactic acid bacteria experiments in bighead carp breeding, and the excess lactic acid bacteria may decrease the pH value of the water, thereby causing a certain stress to the HP. Therefore, more experiments are needed to verify this hypothesis. Currently, there are insufficient studies on the relationship between microorganisms and non-specific immunity in fish, thus, the research on the function and effect of fish intestinal microorganism needs to be strengthened.

5. Conclusions

In conclusion, we observed that the enzyme activities of bighead carp were modulated under different feeding strategies, especially when lactic acid was added to the feed. As a result of enzymatic analysis, compared with group B, the activities of LZM and GR in group C mainly increased in the spleen, while the activities of LZM, GSH-Px, and SOD decreased in the hepatopancreas. After the 16S rRNA gene sequencing data were analyzed for bighead carp with different strategies, we found higher bacterial α-diversities of the intestinal flora in the fish fed by formulated food with/without lactic acid bacteria, compared to group A. Moreover, we identified that three phyla—Proteobacteria, Actinobacteria and Firmicutes—were dominant in intestinal flora of bighead carp among the three groups. Importantly, when lactic acid bacteria were added to the feed, the relative abundance of Actinobacteria increased while the relative abundance of Proteobacteria and Firmicutes decreased. Correlation analysis identified 27 significantly correlated pairs (14 positive and 13 negative correlations) between the intestinal microbiome and the metabolic functions. In particular, LZM was found to affect the non-specific immune function of bighead carp in the spleen, and the genus Cetobacterium may be beneficial for non-specific immunity in fish. This study reveals the complex relationship between the gut microbes and non-specific immunity of bighead carp; the result will improve our understanding of the composition of bighead carp intestinal microbes under new rearing strategies.

Supplementary Materials

The following are available online at https://www.mdpi.com/article/10.3390/genes12060916/s1. Table S1. Enzyme activities in the hepatopancreas and spleen of bighead carp. Table S2. The status of sequencing data for each sample. Table S3. The bacteria community composition identified in the bighead with different feeding strategies on different taxonomic levels. Table S4. Differentially expressed bacteria identified in the intestinal flora of bighead carp with different feed strategies. Table S5. Pearson correlation analysis of enzyme activities and genera in the bighead carp from group C.

Author Contributions

Conceptualization, J.Y., Z.W., B.W., and Y.Z.; methodology, J.Y., H.M., and X.Z.; formal analysis, J.Y., Z.Z., and M.C.; investigation, J.Y., Z.W., and B.W.; resources, J.Y. and Z.W.; data curation, Z.Z.; writing–original draft preparation, J.Y. and M.C.; writing—review and editing, M.C. and Y.Z.; supervision, Y.Z.; project administration, Y.Z.; funding acquisition, Y.Z. All authors have read and agreed to the published version of the manuscript.

Funding

This research was founded by the Chongqing Science and Technology Bureau (cstc2017shms-zdyf0343).

Institutional Review Board Statement

This study was approved by Institutional Animal Care and Use Committee (IACUC) of Southwest University on 20 July 2020 (IACUC No. Approved: 20200720-01).

Informed Consent Statement

Not applicable.

Data Availability Statement

The 16S rRNA sequencing datasets of bighead carp intestinal bacteria was submitted to NCBI SRA (BioProject ID: PRJNA695396).

Acknowledgments

The authors would like to thank and acknowledge Tong Xu and Chengdu Life Baseline Company for the assistance with bioinformatics analysis.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Dong, S.; Li, D. Comparative studies on the feeding selectivity of silver carp Hypophthalmichthys molitrix and bighead carp Aristichthys nobilis. J. Fish Biol. 1994, 44, 621–626. [Google Scholar] [CrossRef]
  2. Anderson, K.R.; Chapman, D.C.; Hayer, C.A. Assessment of dreissenid biodeposits as a potential food resource for invasive Asian carp. Bioinvasions Rec. 2016, 5, 251–257. [Google Scholar] [CrossRef]
  3. Mi, H.; Wen, Y.; Ge, X.; Sui, R.; Liao, L.; Li, B.; Ren, M.; Zeng, B. The aquaculture situation and development trend of bighead carp in the Pearl River Delta. Sci. Fish Farming 2016, 32, 82–84. [Google Scholar]
  4. Huang, A. The technology of breeding bighead with compound feed. Fish. Guide Rich 2008, 000, 25–26. [Google Scholar]
  5. Jorgensen, L.V.G. Zebrafish as a Model for Fish Diseases in Aquaculture. Pathogens 2020, 9, 609. [Google Scholar] [CrossRef]
  6. Vine, N.G.; Leukes, W.D.; Horst, K. Probiotics in marine larviculture. FEMS Microbiol. Rev. 2006, 30, 404–427. [Google Scholar] [CrossRef] [Green Version]
  7. Rengpipat, S.; Rukpratanporn, S.; Piyatiratitivorakul, S.; Menasaveta, P. Immunity enhancement in black tiger shrimp (Penaeus monodon) by a probiont bacterium (Bacillus S11). Aquaculture 2000, 191, 271–288. [Google Scholar] [CrossRef]
  8. Huang, Q.; Wang, J. The mechanism and application of feed lactic acid bacteria probiotics. Chin. J. Anim. Nutr. 2002, 4, 13–18. [Google Scholar]
  9. Son, V.M.; Chang, C.C.; Wu, M.C.; Guu, Y.K.; Chiu, C.H.; Cheng, W. Dietary administration of the probiotic, Lactobacillus plantarum, enhanced the growth, innate immune responses, and disease resistance of the grouper Epinephelus coioides. Fish Shellfish Immunol. 2009, 26, 691–698. [Google Scholar] [CrossRef] [PubMed]
  10. Aly, S.M.; Abdel-Galil Ahmed, Y.; Abdel-Aziz Ghareeb, A.; Mohamed, M.F. Studies on Bacillus subtilis and Lactobacillus acidophilus, as potential probiotics, on the immune response and resistance of Tilapia nilotica (Oreochromis niloticus) to challenge infections. Fish Shellfish Immunol. 2008, 25, 128–136. [Google Scholar] [CrossRef] [PubMed]
  11. Suzer, C.; Çoban, D.; Kamaci, H.O.; Saka, Ş.; Firat, K.; Otgucuoğlu, Ö.; Küçüksari, H. Lactobacillus spp. bacteria as probiotics in gilthead sea bream (Sparus aurata, L.) larvae: Effects on growth performance and digestive enzyme activities. Aquaculture 2008, 280, 140–145. [Google Scholar] [CrossRef]
  12. De Gregorio, E.; Spellman, P.T.; Tzou, P.; Rubin, G.M.; Lemaitre, B. The Toll and Imd pathways are the major regulators of the immune response in Drosophila. EMBO J. 2002, 21, 2568–2579. [Google Scholar] [CrossRef] [Green Version]
  13. Perez, T.; Balcazar, J.L.; Ruiz-Zarzuela, I.; Halaihel, N.; Vendrell, D.; de Blas, I.; Muzquiz, J.L. Host-microbiota interactions within the fish intestinal ecosystem. Mucosal Immunol. 2010, 3, 355–360. [Google Scholar] [CrossRef] [PubMed]
  14. Radovanović, T.; Borković-Mitić, S.; Perendija, B.R.; Despotović, S.; Saičić, Z. Superoxide dismutase and catalase activities in the liver and muscle of barbel (Barbus barbus) and its intestinal parasite (Pomphoryinchus laevis) from the Danube river, Serbia. Arch. Biol. Sci. 2014, 62, 97–105. [Google Scholar] [CrossRef]
  15. Paolicchi, F.; Perea, J.; Cseh, S.; Morsella, C. Relationship between Paratuberculosis and the microelements Copper, Zinc, Iron, Selenium and Molybdenum in Beef Cattle. Braz. J. Microbiol. 2013, 44, 153–160. [Google Scholar] [CrossRef]
  16. Saurabh, S.; Sahoo, P.K. Lysozyme: An important defence molecule of fish innate immune system. Aquac. Res. 2010, 39, 223–239. [Google Scholar] [CrossRef]
  17. Gajardo, K.; Rodiles, A.; Kortner, T.M.; Krogdahl, A.; Bakke, A.M.; Merrifield, D.L.; Sorum, H. A high-resolution map of the gut microbiota in Atlantic salmon (Salmo salar): A basis for comparative gut microbial research. Sci. Rep. 2016, 6, 30893. [Google Scholar] [CrossRef] [Green Version]
  18. Gatesoupe, F.J. Live yeasts in the gut: Natural occurrence, dietary introduction, and their effects on fish health and development. Aquaculture 2007, 267, 20–30. [Google Scholar] [CrossRef] [Green Version]
  19. Saha, S.; Roy, R.N.; Sen, S.K.; Ray, A.K. Characterization of cellulase-producing bacteria from the digestive tract of tilapia, Oreochromis mossambica (Peters) and grass carp, Ctenopharyngodon idella (Valenciennes). Aquac. Res. 2010, 37, 380–388. [Google Scholar] [CrossRef]
  20. Torsvik, V.; Ovreas, L. Microbial diversity and function in soil: From genes to ecosystems. Curr. Opin. Microbiol. 2002, 5, 240–245. [Google Scholar] [CrossRef]
  21. Fjellheim, A.J.; Playfoot, K.J.; Skjermo, J.; Vadstein, O. Vibrionaceae dominates the microflora antagonistic towards Listonella anguillarum in the intestine of cultured Atlantic cod (Gadus morhua L.) larvae. Aquaculture 2007, 269, 98–106. [Google Scholar] [CrossRef]
  22. Caporaso, J.G.; Lauber, C.L.; Walters, W.A.; Berg-Lyons, D.; Lozupone, C.A.; Turnbaugh, P.J.; Fierer, N.; Knight, R. Global patterns of 16S rRNA diversity at a depth of millions of sequences per sample. Proc. Natl. Acad. Sci. USA 2011, 108, 4516–4522. [Google Scholar] [CrossRef] [Green Version]
  23. Gu, C.; Yang, Y.; Xiang, H.; Li, S.; Liang, L.; Sui, H.; Zhan, L.; Lu, X. Deciphering bacterial community changes in zucker diabetic fatty rats based on 16S rRNA gene sequences analysis. Oncotarget 2016, 7, 48941–48952. [Google Scholar] [CrossRef] [Green Version]
  24. Sun, Q.; Li, A.; Li, M.; Hou, B. Effect of pH on biodiesel production and the microbial structure of glucose-fed activated sludge. Int. Biodeterior. Biodegrad. 2015, 104, 224–230. [Google Scholar] [CrossRef]
  25. Edgar, R.C. UPARSE: Highly accurate OTU sequences from microbial amplicon reads. Nat. Methods 2013, 10, 996–998. [Google Scholar] [CrossRef]
  26. Edgar, R.C.; Haas, B.J.; Clemente, J.C.; Quince, C.; Knight, R. UCHIME improves sensitivity and speed of chimera detection. Bioinformatics 2011, 27, 2194–2200. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  27. Schloss, P.D.; Westcott, S.L. Assessing and improving methods used in operational taxonomic unit-based approaches for 16S rRNA gene sequence analysis. Appl. Environ. Microbiol. 2011, 77, 3219–3226. [Google Scholar] [CrossRef] [Green Version]
  28. Westcott, S.L.; Schloss, P.D. De novo clustering methods outperform reference-based methods for assigning 16S rRNA gene sequences to operational taxonomic units. Peer J. 2015, 3, e1487. [Google Scholar] [CrossRef]
  29. Schloss, P.D.; Westcott, S.L.; Ryabin, T.; Hall, J.R.; Hartmann, M.; Hollister, E.B.; Lesniewski, R.A.; Oakley, B.B.; Parks, D.H.; Robinson, C.J.; et al. Introducing mothur: Open-source, platform-independent, community-supported software for describing and comparing microbial communities. Appl. Environ. Microbiol. 2009, 75, 7537–7541. [Google Scholar] [CrossRef] [Green Version]
  30. Segata, N.; Izard, J.; Waldron, L.; Gevers, D.; Miropolsky, L.; Garrett, W.S.; Huttenhower, C. Metagenomic biomarker discovery and explanation. Genome Biol. 2011, 12, R60. [Google Scholar] [CrossRef] [Green Version]
  31. Ijaz, M.U.; Ahmed, M.I.; Zou, X.; Hussain, M.; Zhang, M.; Zhao, F.; Xu, X.; Zhou, G.; Li, C. Beef, Casein, and Soy Proteins Differentially Affect Lipid Metabolism, Triglycerides Accumulation and Gut Microbiota of High-Fat Diet-Fed C57BL/6J Mice. Front. Microbiol. 2018, 9, 2200. [Google Scholar] [CrossRef]
  32. Rawls, J.F.; Samuel, B.S.; Gordon, J.I. Gnotobiotic zebrafish reveal evolutionarily conserved responses to the gut microbiota. Proc. Natl. Acad. Sci. USA 2004, 101, 4596–4601. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  33. Li, T.; Long, M.; Ji, C.; Shen, Z.; Gatesoupe, F.J.; Zhang, X.; Zhang, Q.; Zhang, L.; Zhao, Y.; Liu, X.; et al. Alterations of the gut microbiome of largemouth bronze gudgeon (Coreius guichenoti) suffering from furunculosis. Sci. Rep. 2016, 6, 30606. [Google Scholar] [CrossRef]
  34. Stagaman, K.; Burns, A.R.; Guillemin, K.; Bohannan, B.J. The role of adaptive immunity as an ecological filter on the gut microbiota in zebrafish. ISME J. 2017, 11, 1630–1639. [Google Scholar] [CrossRef]
  35. Balcázar, J.L.; De, B.I.; Ruiz-Zarzuela, I.; Cunningham, D.; Vendrell, D.; Múzquiz, J.L. The role of probiotics in aquaculture. Vet. Microbiol. 2006, 114, 173–186. [Google Scholar] [CrossRef]
  36. Lin, X.; Xia, Y.; Wang, G.; Yang, Y.; Xiong, Z.; Lv, F.; Zhou, W.; Ai, L. Lactic Acid Bacteria With Antioxidant Activities Alleviating Oxidized Oil Induced Hepatic Injury in Mice. Front. Microbiol. 2018, 9, 2684. [Google Scholar] [CrossRef] [PubMed]
  37. Muegge, B.D.; Kuczynski, J.; Knights, D.; Clemente, J.C.; Gonzalez, A.; Fontana, L.; Henrissat, B.; Knight, R.; Gordon, J.I. Diet drives convergence in gut microbiome functions across mammalian phylogeny and within humans. Science 2011, 332, 970–974. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  38. Wu, G.D.; Chen, J.; Hoffmann, C.; Bittinger, K.; Chen, Y.Y.; Keilbaugh, S.A.; Bewtra, M.; Knights, D.; Walters, W.A.; Knight, R.; et al. Linking long-term dietary patterns with gut microbial enterotypes. Science 2011, 334, 105–108. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  39. Ma, C.; Chen, C.; Jia, L.; He, X.; Zhang, B. Comparison of the intestinal microbiota composition and function in healthy and diseased Yunlong Grouper. AMB Express 2019, 9, 187. [Google Scholar] [CrossRef] [Green Version]
  40. Sullam, K.E.; Essinger, S.D.; Lozupone, C.A.; O’Connor, M.P.; Rosen, G.L.; Knight, R.; Kilham, S.S.; Russell, J.A. Environmental and ecological factors that shape the gut bacterial communities of fish: A meta-analysis. Mol. Ecol. 2012, 21, 3363–3378. [Google Scholar] [CrossRef] [Green Version]
  41. Estruch, G.; Collado, M.C.; Penaranda, D.S.; Tomas Vidal, A.; Jover Cerda, M.; Perez Martinez, G.; Martinez-Llorens, S. Impact of Fishmeal Replacement in Diets for Gilthead Sea Bream (Sparus aurata) on the Gastrointestinal Microbiota Determined by Pyrosequencing the 16S rRNA Gene. PLoS ONE 2015, 10, e0136389. [Google Scholar] [CrossRef]
  42. Ricaud, K.; Rey, M.; Plagnes-Juan, E.; Larroquet, L.; Even, M.; Quillet, E.; Skiba-Cassy, S.; Panserat, S. Composition of Intestinal Microbiota in Two Lines of Rainbow Trout (Oncorhynchus Mykiss) Divergently Selected for Muscle Fat Content. Open Microbiol. J. 2018, 12, 308–320. [Google Scholar] [CrossRef] [PubMed]
  43. Guarner, F.; Malagelada, J.R. Gut flora in health and disease. Lancet 2003, 361, 512–519. [Google Scholar] [CrossRef]
  44. Penn, K.; Jenkins, C.; Nett, M.; Udwary, D.W.; Gontang, E.A.; McGlinchey, R.P.; Foster, B.; Lapidus, A.; Podell, S.; Allen, E.E.; et al. Genomic islands link secondary metabolism to functional adaptation in marine Actinobacteria. ISME J. 2009, 3, 1193–1203. [Google Scholar] [CrossRef]
  45. Mendes, M.C.S.; Paulino, D.S.; Brambilla, S.R.; Camargo, J.A.; Persinoti, G.F.; Carvalheira, J.B.C. Microbiota modification by probiotic supplementation reduces colitis associated colon cancer in mice. World J. Gastroenterol. 2018, 24, 1995–2008. [Google Scholar] [CrossRef]
  46. Eichmiller, J.J.; Hamilton, M.J.; Staley, C.; Sadowsky, M.J.; Sorensen, P.W. Environment shapes the fecal microbiome of invasive carp species. Microbiome 2016, 4, 44. [Google Scholar] [CrossRef] [Green Version]
  47. Finegold, S.M.; Vaisanen, M.L.; Molitoris, D.R.; Tomzynski, T.J.; Song, Y.; Liu, C.; Collins, M.D.; Lawson, P.A. Cetobacterium somerae sp. nov. from human feces and emended description of the genus Cetobacterium. Syst. Appl. Microbiol. 2003, 26, 177–181. [Google Scholar] [CrossRef]
  48. Sugita, H.; Miyajima, C.; Deguchi, Y. The vitamin B12-producing ability of the intestinal microflora of freshwater fish. Aquaculture 1991, 92, 267–276. [Google Scholar] [CrossRef]
Figure 1. α-diversity Shannon (a) and Simpson (b) indexes based on average OTUs of bighead carp from the three groups. The values on the horizontal line represent the p value of the Wilcoxon signed-rank test. Error bars indicate the SD values. (c) Principal co-ordinate analysis (Bray–Curtis distance) of β-diversity. PCoA1, PCoA2, and PcoA3 represent the top three principal coordinates that, together, captured most of the diversity. (d) Boxplot for comparisons of intra-group and inter-group based on group distances. A, nature live food only; B, nature live food + fish feed; C, nature live food + fish feed + lactic acid bacteria.
Figure 1. α-diversity Shannon (a) and Simpson (b) indexes based on average OTUs of bighead carp from the three groups. The values on the horizontal line represent the p value of the Wilcoxon signed-rank test. Error bars indicate the SD values. (c) Principal co-ordinate analysis (Bray–Curtis distance) of β-diversity. PCoA1, PCoA2, and PcoA3 represent the top three principal coordinates that, together, captured most of the diversity. (d) Boxplot for comparisons of intra-group and inter-group based on group distances. A, nature live food only; B, nature live food + fish feed; C, nature live food + fish feed + lactic acid bacteria.
Genes 12 00916 g001
Figure 2. Global composition of bacteria at different taxonomic levels. (a) Phyla level (b) order level (c) genus level. A, nature live food only; B, nature live food + fish feed; C, nature live food + fish feed + lactic acid bacteria. Only the top 15 organisms are presented.
Figure 2. Global composition of bacteria at different taxonomic levels. (a) Phyla level (b) order level (c) genus level. A, nature live food only; B, nature live food + fish feed; C, nature live food + fish feed + lactic acid bacteria. Only the top 15 organisms are presented.
Genes 12 00916 g002
Figure 3. Bacteria taxa differentially expressed in fish intestinal microbiota of the three groups of bighead carp. (a) Differentially expressed bacteria identified in group B (nature live food + fish feed, n = 15) compared to group A (nature live food only, n = 9). (b) Differentially expressed bacteria identified in group C (nature live food + fish feed + lactic acid bacteria, n = 15). The histograms represent the linear discriminant analysis (LDA) score computed for taxa differentially abundant between the different groups. *** p < 0.001; ** p < 0.01; * p < 0.05.
Figure 3. Bacteria taxa differentially expressed in fish intestinal microbiota of the three groups of bighead carp. (a) Differentially expressed bacteria identified in group B (nature live food + fish feed, n = 15) compared to group A (nature live food only, n = 9). (b) Differentially expressed bacteria identified in group C (nature live food + fish feed + lactic acid bacteria, n = 15). The histograms represent the linear discriminant analysis (LDA) score computed for taxa differentially abundant between the different groups. *** p < 0.001; ** p < 0.01; * p < 0.05.
Genes 12 00916 g003
Figure 4. Scatter plot (ae) showing the top 5 Pearson correlation coefficient between genus taxonomy and enzyme activities. a, in hepatopancreas; b–e, in the spleen. X-axis, enzyme activity indicator; Y-axis, log10 value of OUT abundance. LZM, lysozyme; GR, glutathione reductase.
Figure 4. Scatter plot (ae) showing the top 5 Pearson correlation coefficient between genus taxonomy and enzyme activities. a, in hepatopancreas; b–e, in the spleen. X-axis, enzyme activity indicator; Y-axis, log10 value of OUT abundance. LZM, lysozyme; GR, glutathione reductase.
Genes 12 00916 g004
Table 1. Characteristics and non-specific immunity related enzyme index of bighead carp.
Table 1. Characteristics and non-specific immunity related enzyme index of bighead carp.
Characteristics/TissueIndexGroup A (n = 9)Group B (n = 15)Group C (n = 15)p-Value
CharacteristicsInitial weight (g)47.57 ± 6.0852.40 ± 1.5551.87 ± 2.13NS
Final weight (g)77.85 ± 11.86121.80 ± 2.90120.54 ± 5.83<0.01
Weight gain rate (%)63.40 ± 4.04132.46 ± 1.75132.34 ± 1.67<0.001
Growth rate (%)50.47 ± 9.63115.66 ± 2.37114.44 ± 6.16<0.001
Survival rate (%)100100100NS
Fatness (%)1.97 ± 0.092.00 ± 0.092.00 ± 0.08NS
Visceral body ratio (%)3.77 ± 0.546.17 ± 0.746.14 ± 0.64<0.01
HepatopancreasLZM (ug/mL)3.06 ± 0.39 a1.92 ± 0.28 b0.57 ± 0.02 c<0.001
CAT (U/mgprot)13.75 ± 2.3711.00 ± 0.6511.86 ± 0.34NS
GR (U/gprot)18.36 ± 3.8415.35 ± 3.7218.46 ± 1.01NS
GSH-PX (U/mgprot)1182.86 ± 65.99 a1162.77 ± 35.16 a760.89 ± 49.84 b<0.001
SOD (U/mgprot)34.85 ± 2.58 c130.29 ± 8.76 a86.57 ± 2.98 b<0.001
SpleenLZM (ug/mL)14.66 ± 0.31 b16.82 ± 1.14 b22.27 ± 2.22 a<0.01
CAT (U/mgprot)1.11 ± 0.06 b2.84 ± 0.30 a1.07 ± 0.29 b<0.001
GR (U/gprot)7.61 ± 0.75 c14.83 ± 1.43 b20.44 ± 0.78 a<0.001
GSH-PX (U/mgprot)670.89 ± 27.66 a449.05 ± 38.08 b452.16 ± 15.25 b<0.001
SOD (U/mgprot)381.52 ± 0.99 b449.27 ± 14.30 a450.89 ± 18.71 a<0.01
Note: values are expressed in mean ± SD. LZM, lysozyme; CAT, catalase; GR, glutathione reductase; GSH-PX, glutathione peroxidase; SOD, superoxide dismutase; HP, hepatopancreas; NS, not significant; a, b and c in the cells denote significant difference of group A, B and C, respectively, among the three groups (p < 0.05).
Table 2. Single sample diversity α-diversity analysis (mean +/− SD).
Table 2. Single sample diversity α-diversity analysis (mean +/− SD).
GroupOTUsChaoShannonSimpsonCoverage
A1806 +/− 2673111.11 +/− 543.804.20 +/− 0.470.1084 +/− 0.06590.9862 +/− 0.0029
B1764 +/− 3032844.00 +/− 562.984.52+/− 0.600.075 +/− 0.04990.9871 +/− 0.0027
C2298 +/− 1963739.67 +/− 315.295.16 +/− 0.2910.0393 +/− 0.01370.9837 +/− 0.0026
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Yuan, J.; Wang, Z.; Wang, B.; Mei, H.; Zhai, X.; Zhuang, Z.; Chen, M.; Zhang, Y. Non-Specific Immunity Associated Gut Microbiome in Aristichthys nobilis under Different Rearing Strategies. Genes 2021, 12, 916. https://doi.org/10.3390/genes12060916

AMA Style

Yuan J, Wang Z, Wang B, Mei H, Zhai X, Zhuang Z, Chen M, Zhang Y. Non-Specific Immunity Associated Gut Microbiome in Aristichthys nobilis under Different Rearing Strategies. Genes. 2021; 12(6):916. https://doi.org/10.3390/genes12060916

Chicago/Turabian Style

Yuan, Jianming, Zhijian Wang, Bo Wang, Huiqing Mei, Xuliang Zhai, Zhenhua Zhuang, Maoshan Chen, and Yaoguang Zhang. 2021. "Non-Specific Immunity Associated Gut Microbiome in Aristichthys nobilis under Different Rearing Strategies" Genes 12, no. 6: 916. https://doi.org/10.3390/genes12060916

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop