Investigation into Host Selection of the Cecal Acetogen Population in Rabbits after Weaning

Abstract

Homoacetogenic bacteria have received attention as a hydrogenotrophic population that offers a significant energetic advantage to the host animal. Reductive acetogenesis is likely an important hydrogen disposal mechanism in the cecum of rabbits. However, molecular ecology information about cecal acetogen candidates has rarely been reported. To better understand the effect of host selection in the rabbit cecal acetogen community with respect to growth, rabbits at four different age stages (30, 60, 120 and 180 days) with the same diet were studied. Although the abundance of potential acetogens and methanogens was high in the cecum of rabbits undergoing growth, many novel potential acetogen populations were observed in the cecum of rabbits across all age groups. Young and adult rabbits had their own distinct acetogen community although they received the same diet, which suggests that as the rabbit ages, acetogens in the cecum undergo developmental changes because of host selection that are independent of diet, and perhaps the different acetogen communities result in different hydrogenotrophic characteristics. The within-group similarity increased with age, indicating that the acetogen community converges to a more homogeneous and stable arrangement with aging.

About the Authors:

Chunlei Yang

Affiliation: Institute of Dairy Science, MoE Key Laboratory of Molecular Animal Nutrition, College of Animal Sciences, Zhejiang University, Hangzhou 310058, China

Lan Mi

Affiliation: Institute of Dairy Science, MoE Key Laboratory of Molecular Animal Nutrition, College of Animal Sciences, Zhejiang University, Hangzhou 310058, China

Xialu Hu

Affiliation: Institute of Dairy Science, MoE Key Laboratory of Molecular Animal Nutrition, College of Animal Sciences, Zhejiang University, Hangzhou 310058, China

Jianxin Liu

Affiliation: Institute of Dairy Science, MoE Key Laboratory of Molecular Animal Nutrition, College of Animal Sciences, Zhejiang University, Hangzhou 310058, China

Jiakun Wang

* E-mail: [email protected]

Affiliation: Institute of Dairy Science, MoE Key Laboratory of Molecular Animal Nutrition, College of Animal Sciences, Zhejiang University, Hangzhou 310058, China

ORCID http://orcid.org/0000-0002-7213-3721

Introduction

Homoacetogenic bacteria (acetogens) are a group of obligate anaerobic bacteria that utilize the acetyl coenzyme-A (CoA) pathway to synthesize acetate from syngas [1]. Acetogenesis is of great importance to the global carbon cycle. Approximately 10¹³ kg of acetate is formed annually in anaerobic habitats [2]. Acetogens are phylogenetically quite diverse and metabolically versatile [3, 4], but only some strains, such as Moorella thermoacetica [5, 6], Clostridium ljungdahlii [7, 8], Clostridium ragsdalei [8, 9] and their genetic modifications, are used in syngas fermentation for biofuels. Furthermore, homoacetogenesis is a promising pathway to compete with methanogens in the rumen, because methane producers cause 23% of the global anthropogenic methane emissions [10]. A better understanding of the molecular ecology of the acetogen population will help develop new acetogen products and enhance its catalyst function.

Methanogens have a greater thermodynamic advantage than acetogens when competing for hydrogen in anaerobic habitats [11–13]. Rabbits, as herbivores, have less methane production [14] and lower energy loss from methane production per unit of body mass than ruminants [15]. The estimated hydrogen recoveries for methane are 24.7% and 85.4% in the rabbit cecum [16] and goat rumen [17], respectively, and the ratio of acetate to propionate production is much higher in the rabbit cecum compared to rumens (17.8 vs. 2.81) [18]. The reason for this may be because reductive acetogenesis is the dominant hydrogen disposal pathway in the cecum of rabbits [19, 20]. Therefore, acetogens that are, more efficient at syngas fermentation may exist in the cecum of rabbits.

The structure and composition of the gut microbiota is driven by lifestyle strategies, such as the growth rate, substrate utilization patterns, and host selection for specific bacteria with emergent collective behavior that is beneficial to the host [21]. Studies on the bacterial 16S rRNA gene revealed that the fecal microbiota of mammals is specific and rather stable for their particular host species, to a large extent [22–24], suggesting that mechanisms exist to recruit and maintain selected bacterial populations.

The cellulolytic capability of the gut microbiota is enhanced with an increase in the growth stage of the host [25, 26], and hence, more hydrogen is produced. However, there is a lack of knowledge on the stability or fluctuations in the cecum of rabbits maintained on a uniform diet and the selection of the acetogen population under this stability or fluctuation. Knowledge of the acetogen distribution with growth stage will help us understand the microbial host selection process and develop new acetogen products in the future. Therefore, host selection of the cecal acetogen population was studied in the cecum of rabbits at different growth stages maintained under constant conditions, including a uniform diet.

Materials and Methods

Experimental design and sampling

The Animal Care Committee of Zhejiang University (Hangzhou, China) approved all experiments, and the experimental procedures used in this study were in accordance with the university’s guidelines for animal research. Rabbits at four different age stages were fed the same diets and were used to study the change in the acetogen population in the cecum. Six male New Zealand White rabbits at ages 30, 60, 120 and 180 days were purchased simultaneously from the Zhejiang Academy of Agricultural Science. They were housed in indoor three-layer cages (60×50×35 cm) with natural lighting and raised with the same commercial pellet diets consisting of 12% corn, 18% bran, 8% soya bean cake, 31% meal, 10% malt root, 16% chaff, 0.8% powder, 0.5% salt and 4% commercial additives. Only one rabbit was housed in each individual cage. All rabbits were sacrificed on the same day. The performance of euthanasia and removal of cecum were undertaken by trained animal technician (license number: 15128, issued by Zhejiang University Laboratory Animal Center). Body weight was measured for each rabbit to ensure intravenous injection of 100 mg/kg (body weight) phenobarbital sodium (Sigma, Saint Louis, MO, USA) [27]. Once the rabbits have lost toe pinch, leg withdrawal reflex, ear intravenous injection of 20 ml air [28] was then performed to accelerate the death of rabbit. After sacrifice, the cecum was immediately removed and isolated by tying off the extremities with a nylon string to prevent movement of the digesta. The contents of the proximal, medial and distal parts of the cecum were squeezed into 50-ml sterilized eppendorf tubes and mixed for pH measurement. The samples were stored on dry ice before shipment to laboratory where they were then stored at -80°C.

Measurement of cecal fermentation parameters

The samples for the measurement of volatile fatty acids (VFA) and ammonia nitrate were prepared as described by Yang et al. [29]. Approximately two grams of each stored sample was vortexed with 9 ml sterile PBS (pH 7.0). Then the supernatant was obtained after centrifugation at 13,000 × g at 4°C for 15 min. Colorimetry [30] and gas liquid chromatography (GC-2010, Shimadzu, Kyoto, Japan) [31] were used to determine the ammonia nitrate and VFA concentrations, respectively.

DNA extraction and real-time PCR quantification

Total DNA was extracted from 0.1 gram of each of the stored samples using the cetyltrimethylammonium bromide (CTAB) method as described by Gagen et al. [32]. Primers (Table 1) specific for the formyltetrahydrofolate synthetase gene (fhs) [33] and the methyl coenzyme-M reductase A gene (mcrA) [34] were used to amplify potential acetogens and methanogens, respectively. The 16S rRNA gene from total bacteria was amplified with primers as reported by Denman and McSweeney [35]. Real-time PCR was performed with SYBR green in an ABI 7500 (Life Technologies, Singapore) using the following program: one cycle of initial denaturation at 95°C for 10 s, followed by 40 cycles of denaturation at 95°C for 15 s and annealing at 60°C for 1 min. The standards were prepared from serial dilutions of the linearized plasmid containing between 10² and 10⁹ target gene copies calculated from the concentration of plasmids. The copy numbers of the standards and the copy numbers of the target genes per gram of cecal sample were calculated according to the method of Ritalahti et al. [36].

[Figure omitted. See PDF.]

Table 1. Primers used in this study.

https://doi.org/10.1371/journal.pone.0158768.t001

Analysis of the acetogen community

Amplification of subunit B of the acetyl-CoA synthasegene (acsB) from each extracted DNA sample was performed using the TaKaRa PCR Kit (TaKaRa, Dalian, China) with 1.0 ng template DNA per reaction. The specific primers for acsB [32] (Table 1) were used with the addition of 6-base barcodes and an Illumina adapter sequence [37]. Amplification was performed on the ABI 9700 PCR System using the conditions described by Gagen et al. [32]. The PCR products for each sample were separated from the primers and primer dimers by electrophoresis on a 1.5% agarose gel and the QIAquick gel extraction kit (QIAGEN, Shanghai, China). Then, the products for each sample were mixed in equal ng quantities using Qubit 2.0 (Invitrogen, USA) and sent to a commercial company (Novogene, Beijing, China) for Illumina MiSeq sequencing.

After sequencing, the size filtering, quality control (Q>20) and chimera removal of the raw data were performed using the pipeline in QIIME [38]. The sequences that passed quality control were aligned in Mega version 4.0 [39] to strictly select the sequences that were approximately 216 bp long and contained the forward and reverse acsB primers. Then, these nucleotide sequences were translated into amino acid sequences in a specific open reading frame using the translate tool in ExPASy [40]. The obtained ACS amino acid sequences were aligned again and length-based filtering was performed. The final obtained ACS amino acid sequences were clustered into operational taxonomic units (OTUs) at 0.05 distance limit using the CD-HIT program [41]. Additional analysis of OTUs was performed in the R packages phyloseq, DEseq, ggplot2 and gplots [42–45]. RAxML version 7.0.3 [46] and ARB [47] were used to construct the maximum-likelihood trees of the deduced ACS amino acids as described by Gagen et al. [32]. The sequences obtained in this study were deposited in the European Nucleotide Archive (ENA) under accession number PRJEB12535.

Statistical analysis

The abundance of the 16S rRNA and fhs and mcrA genes, and the concentration or molar percentage of the total and individual VFA were analyzed by one-way analysis of variance in R 3.2 with individual rabbits as the experimental unit and age as the main effect. Multiple comparisons of means among different age groups were conducted using Turkey multiple range tests. Significance was declared at P<0.05.

Results and Discussion

Change in the cecum environment with growth stage

The pH of the cecum increased with age. However, the total VFA and ammonia nitrate decreased with age (Table 2). The ratio of acetate to propionate did not change significantly but was higher than 9.81:1. The primary activity of the bacterial population that resides in the hindgut is fermentation, and the major products are VFA. These can be absorbed across the gut mucosa and utilized by the rabbit as a significant energy resource [48]. With increased body weight, the energy requirement should also increase; therefore, the decrease in the total VFA in the cecum of older rabbits may be a result of increased absorption. The high ratio of acetate to propionate further confirms that acetogenesis is an important metabolic pathway present in the rabbit cecum.

[Figure omitted. See PDF.]

Table 2. Fermentation parameters of cecal content.

https://doi.org/10.1371/journal.pone.0158768.t002

Comparison of the acetogen community

The deduced amino acids for acsB from the cecum of rabbit clustered into 6101 OTUs, suggesting that many acetogen species exist in the cecum of rabbits. However, most of the detected sequences lack representatives from previously described sequences, which is in consistent with previous studies [32, 49], suggesting that the novelty of the rabbit cecal flora is likely because the rabbit is a cecotrophic species. In addition, these novel acsB sequences amplified in our study expand the database of ACS amino acid sequences for acetogens. Alpha diversity, which is described in terms of richness or the number of different species present, and evenness were initially investigated to observe changes in acetogen diversity in the cecum of rabbits during growth and development (Fig 1). Our results showed that both the richness and evenness of acetogen populations in the cecum of rabbit at each growth stage did not change significantly (P>0.05), with no significant changes in the Chao 1 and Shannon indices (Fig 1). However, young rabbits showed a larger animal variance compared to the older rabbits (Fig 1), which indicates that during the first days of life, the acetogen population in the cecum has more heterogeneity and becomes more homogeneous and stable with growth.

[Figure omitted. See PDF.]

Fig 1. Alpha diversity measures of the deduced ACS amino acid sequences recovered from the cecum of rabbits across different age groups.

Observed = Observed OTU numbers; Chao1 = Chao estimator; Shannon = Shannon diversity index.

https://doi.org/10.1371/journal.pone.0158768.g001

The structure of the acetogen population was also analyzed with beta diversity measures using principal coordinates analysis (PCoA) (Fig 2). Clear alterations of the acetogen community in the cecum of rabbit were observed with growth. Animals at different growth stages were divergent in their acetogen composition, suggesting that rabbits at each growth stage have their own distinct acetogen community, although they received the same diet. Variations in the cecal acetogen community between animals at the same growth stage were also observed, suggesting that the effect of rabbit-specific host selection determines the distinct acetogen population in each animal.

[Figure omitted. See PDF.]

Fig 2. Covariation of the community structure of potential acetogens in the cecum of rabbits at different age stages using principal coordinate analysis (PCoA).

The PCoA plot shows that the samples from different age groups span their first two principal components, principal coordinate 1 (29% of variance) and principal coordinate 2 (13% of variance).

https://doi.org/10.1371/journal.pone.0158768.g002

The between-group similarity was further analyzed using a heat map (Fig 3), and the number of presented OTUs was significantly changed between groups at each of the two growth stages (Fig 4). Obvious differences were found in the presented ACS OTUs between the growth stage and host animals (Fig 3) by observing the clusters. The groups aged 30 and 60 days were clearly separated from the groups aged 120 and 180 days, with the exception of one rabbit at 60 days and two rabbits at 120 days. Furthermore, a significant individual effect was observed because rabbits of the same age showed different sub-clusters. From analysis of the number of OTUs that were significantly changed (fold change ≥2) (P<0.05) in every two age groups (Fig 4), the between-group similarity increased in a growth-dependent manner. There were 131 ACS OTUs that showed a greater than 2-fold change between groups aged 30 and 180 days, whereas the changed OTUs were reduced to 35 and 0 when the 180-day group was compared to 60 days and 120 days groups, respectively, suggesting that the cecum of rabbits becomes a more stable arrangement of the acetogen community from youth to adult. Similar results were observed from previous studies focusing on the bacterial community. The pairwise similarity of human bacterial diversity increases with age [50], and the gut microbiome also undergoes a convergence to a mature bacterial arrangement with age [51].

[Figure omitted. See PDF.]

Fig 3. Assessment of differentially expressed ACS genes in samples across all age groups in heat map format.

A heat map showing the rlog transformed expression data of the first 50 ACS amino acid sequences in the top ranked ACS sequences detected in the cecum of rabbits of different ages. The color scale ranges from blue (high expression) to white (low expression).

https://doi.org/10.1371/journal.pone.0158768.g003

[Figure omitted. See PDF.]

Fig 4. Analysis of the within-group similarity of two different age groups.

The higher number of significantly changed expressed deduced ACS amino acid sequences between two age groups with a lower similarity between them.

https://doi.org/10.1371/journal.pone.0158768.g004

Phylogenetic analysis of the acetogen community

Due to the limited acsB gene database, most of the detected novel ACS OTUs in the cecum of rabbits could not be classified at the species or genus level. Therefore, a maximum likelihood tree (Fig 5) was constructed to determine the phylogenetic information of acetogens at the family level. In a previous analysis, the groups of 30 and 180 days showed the largest between-group difference. To identify the potential link between acetogen species and young and adult growth stages, the presence of OTUs that changed significantly and the dominant OTUs that were not different in the 30- and 180-day groups were included in the maximum likelihood tree (Fig 5). The most abundant OTU, OTU 1, was affiliated with the Lachnospiraceae family but was phylogenetically placed far away from any of the known isolates present across all growth stages. OTU 1 accounted for 57.66% of all species detected in our study, showing an unchanged abundance between the 30- and 180-day groups; this was followed by OTU 6 with a proportion of 10.73%, and it was phylogenetically affiliated with the Blautia group (Fig 5). The presence of 85 ACS OTUs (a total of 61 in groups 1, 2, 3, 4, 6, 7 and 10; 12 in group 8; 2 in group 9; as well as OTUs 41, 43, 79, 134, 162, 249, 292, 456, 734, and 1732 rabbit) significantly increased (fold change ≥2) (P<0.05), and 46 ACS OTUs (3 in group 5; 39 in group 8; 1 in group 9; and OTUs 7, 23, and 541 rabbit) significantly decreased (fold change ≥2) (P<0.05) in the adult rabbit aged 180 days (Fig 5). These significantly changed (fold change ≥2) ACS sequences were phylogenetically placed at 23 different phylogenetic places in the maximum-likelihood trees of ACS amino acid sequences (Fig 5). Comparison of the tree topology between the 16S rRNA gene and the ACS amino acid sequences published previously show [32] the presence of a total of 34 ACS OTUs (group 1), including a dominant OTU, OTU 4, which clustered together in the ACS tree placed phylogenetically between the Blautia group and Ruminococcaceae, as well as one OTU in group 2 and OTU 179, which was phylogenetically affiliated with the Blautia group, that significantly (P<0.05) increased in adult rabbits aged 180 days compared to the young rabbits aged 30 days. The presence of OTU 43, which was phylogenetically affiliated with Ruminococcaceae, significantly (P<0.05) increased in older rabbits aged 180 days. In addition, another 13 ACS OTUs (group 10) and OTUs 41 and 456 were phylogenetically affiliated with Clostridiaceae, showing a significant increase (P<0.05) in rabbits aged 180 days. Another 8 clusters (groups 3, 4, 6, and 7 and OTUs 41, 134, 162, 249, 292, 456, 734, and 1732) of ACS sequences containing the dominant OTU 9 were phylogenetically affiliated with the Lachnospiraceae family and were more highly present (P<0.05) in rabbits aged 180 days but were far away from any known cultured acetogens. The presence of OTUs 7, 23, and 541 and group 5 in the trees was affiliated with Lachnospiraceae but were far away from any known cultured acetogens and was higher (P<0.05) in young rabbits aged 30 days. A previous study indicated that over 90% of the bacterial sequences in the cecum of rabbit belong to the Firmicutes phylum [52]. Ruminococcaceae and Lachnospiraceae were the dominant families, accounting for 45% and 35% of all sequences, respectively [53]. In our study, the great diversity of ACS sequences was broadly placed phylogenetically into the family groups Lachnospiraceae, Ruminococcaceae, Blautia, and Clostridiaceae, which were shown from the maximum-likelihood trees of the 16S rRNA gene and the ACS amino acid sequences for comparison to the tree topology reported by Gagen et al. [32]. Adult rabbits at 180 days showed a more diverse phylogenetic acetogen community compared to young rabbits aged 30 days. These phylogenetic differences in the potential acetogen populations in the cecum of young and adult rabbits may contribute to their different roles in hydrogenotrophy and energy utilization in the cecal ecosystem. Previous studies showed that a novel sequence, closely related to the known reductive acetogen Blautia coccoides, was associated with carbon dioxide and hydrogen metabolism [54]. Mixotrophic acetogens, which can simultaneously use hydrogen and organic substrates for growth, may be more energetically competitive with methanogens in some situations [55].

[Figure omitted. See PDF.]

Fig 5. Phylogenetic analysis of significantly changed (fold change ≥2) acetyl-CoA synthase amino acid sequences in the cecum of rabbits aged 30 and 180 days with reference ACS sequences (GenBank accession numbers are shown after the species names).

Bootstrap values that ≥75% are shown at the nodes. The scale bar represents a sequence divergence of 10%. The number of ACS OTUs more highly present in the 30-day group or the 180-day group, as well as the dominant OTUs across all age groups, is indicated in brackets.

https://doi.org/10.1371/journal.pone.0158768.g005

Comparison of bacteria and archaea populations

From analysis of the gene copy number of the 16S rRNA targeting total bacteria, large amounts of bacteria were observed in the cecum of rabbits across all growth stages (P>0.05) (Fig 6), which is consistent with previous studies [56, 57]. Real-time PCR was performed to specifically quantify the gene copy number of fhs and mcrA to determine the potential abundance of acetogens, which is currently the best target for quantifying acetogens [33, 58], and with the mcrA gene to determine the methane producing hydrogenotrophy, which is likely the dominant hydrogen disposal method in most anaerobic environments. The abundance of the mcrA gene was at 10⁶ to 10⁸ gene copy numbers per gram of cecal content in most of the samples, which is consistent with a previous study (approximately 7.4 log₁₀ copy number per gram of cecal content at day 35) [52] and was significantly higher (P<0.05) in samples aged 60 days compared with the other age groups (Fig 6). However, the abundance of the fhs gene was at 10⁵ to 10⁶ gene copy numbers per gram of cecal content in most of the samples and showed no significant differences (P>0.05) with growth stage (Fig 6).

[Figure omitted. See PDF.]

Fig 6. Abundance of the fhs, mcrA, and 16S rRNA genes in total bacteria in samples from different age groups.

Groups with different superscripts indicate that they differ significantly (P<0.05).

https://doi.org/10.1371/journal.pone.0158768.g006

General Discussion

Compared to the published acetogen strains used in microbial production platforms, many new acetogen strains were found in the cecum of rabbits in our study, with the dominant acetogen affiliated with the Blautia family, in which the known strains Blautia coccoides, Blautia hydrogenotrophica and Blautia schinkii can use carbon dioxide and hydrogen to synthesize acetate as a main product [59, 60]. In addition, the minimum doubling time of Blautia producta with carbon monoxide or carbon dioxide and hydrogen as a substrate is only 1.5 or 5 h [61, 62], suggesting that the cecum of the rabbit is a promising resource for harvesting novel acetogen strains with potential industrial applications.

The true extent of acetogen biodiversity in the cecum of rabbits is not known because most organisms cannot be cultured ex vivo. However, there was a clear rabbit-specific effect on the cecal acetogen community in individuals used in this study. Only rabbits of the same species (New Zealand White rabbit) and with the same type of diet were used, which minimized both the genetic variability and dietary variation in the study to examine the effect of growth. However, because of the limited database of known acetogens, the acetogen strain changes with growth could not be classified accurately. Some of the acetogen strains were more present in adult rabbits (180 days) and were affiliated with the Blautia group, Ruminococcaceae and Clostridiaceae. Blautia schinkii and Blautia hydrogenotrophica are known to use carbon dioxide and hydrogen as substrates for growth, and Blautia spp. may play a key role in the kangaroo foregut and are associated with carbon dioxide and hydrogen metabolism [54]. The Clostridiaceae family contains the Clostridium genus, which is now an important candidate in the biofuel industry. Therefore, because of the growth effect from host selection, the cecum of adult rabbits may be a good source for enriching acetogen candidates for biofuel production.

The inter-individual variation of the cecal bacteria in rabbits was previously observed, although acetogen was not specifically tested [52]. Comparable results were also described regarding the bacterial community in the feces of ponies. A significantly different population was found in each pony during 11 weeks of observation with the same diet [63]. In addition, from the analysis of the bacterial community in cattle and human, the diversity and within-group similarity increases with growth, suggesting a convergence to a more diverse, but homogeneous and specific, mature bacterial arrangement with growth [25, 51].

As a hindgut fermenter, the greatest numbers of bacteria reside in the cecum of rabbits. Relationships between the gut bacterial communities and their mammalian hosts contribute to the host’s well-being and proper function [22, 64]. In the cecum, archaea and acetogens play important roles in the final stage of organic matter decay, reducing carbon dioxide into methane or acetate, and the latter plays an important role in host animal nutrition [65]. Although high levels of gene copy number both for mcrA and fhs were observed in our study, a previous study indicated that the cecum of rabbits harbors a methanogenic community with low complexity of the genus Methanobrevibacter, and it is possibly one-species dominated because all sequences detected presented similarity values of 99% or higher [66]. However, many novel and host-specific acetogen populations were observed in our study, suggesting an important role for reductive acetogenesis in the cecum of rabbits.

To better understand the function of acetogens in the cecum of rabbit, further studies, such as tracking of the biological fate of carbon dioxide in the cecum of rabbits using RNA stable isotope probing (RNA-SIP) techniques, should be performed to identify the dominant organisms and pathways involved in hydrogenotrophy of rabbit cecum. This was previously performed in the kangaroo forestomach to explain the dominance of reductive acetogenesis as a hydrogen removal pathway [54].

Conclusion

The cecum of the rabbit is a good resource to identify and develop new acetogen products. Hosts select preferable microbiota for their growth, and the microbiota in the cecum are more homogeneous and stable with growth and development. To better understand the contribution of acetogen populations to the host, stable isotope tracing can be performed to identify the dominant organisms and pathways involved in hydrogenotrophy in the cecum of rabbit.

Acknowledgments

The authors are grateful to Dr. Stuart Denman at CSIRO, Australia, for his assistance in the data analysis of this manuscript. The authors acknowledge the financial support from Natural Science Foundation of China (31372337).

Compliance with ethics guidelines

All management and experimental procedures were approved by the Animal Care Committee of Zhejiang University, China, and were performed in accordance with the University’s guidelines for animal research.

Author Contributions

Conceived and designed the experiments: JW. Performed the experiments: CY LM XH. Analyzed the data: CY. Contributed reagents/materials/analysis tools: JW. Wrote the paper: CY JL JW.

Citation: Yang C, Mi L, Hu X, Liu J, Wang J (2016) Investigation into Host Selection of the Cecal Acetogen Population in Rabbits after Weaning. PLoS ONE 11(7): e0158768. https://doi.org/10.1371/journal.pone.0158768

References

1. Bengelsdorf FR, Straub M, Dürre P. Bacterial synthesis gas (syngas) fermentation. Environ Technol. 2013;34(13–14):1639–51.

2. Wood H. Life with CO or CO2 and H2 as a source of carbon and energy. FASEB J. 1991;5(2):156–63. pmid:1900793

3. Drake HL. Acetogenesis, acetogenic bacteria, and the acetyl-CoA “Wood/Ljungdahl” pathway: past and current perspectives. Acetogenesis; Springer US; 1994. pp. 3–60.

4. Schink B. Diversity, ecology, and isolation of acetogenic bacteria. Acetogenesis; Springer US; 1994. pp. 197–235.

5. Sakai S, Nakashimada Y, Inokuma K, Kita M, Okada H, Nishio N. Acetate and ethanol production from H2 and CO2 by Moorella sp. using a repeated batch culture. J Biosci Bioeng. 2005;99(3):252–8. pmid:16233785

6. Kita A, Iwasaki Y, Sakai S, Okuto S, Takaoka K, Suzuki T, et al. Development of genetic transformation and heterologous expression system in carboxydotrophic thermophilic acetogen Moorella thermoacetica. J Biosci Bioeng. 2013;115(4):347–52. pmid:23177215

7. Köpke M, Held C, Hujer S, Liesegang H, Wiezer A, Wollherr A, et al. Clostridium ljungdahlii represents a microbial production platform based on syngas. Proc Natl Acad Sci USA. 2010;107(29):13087–92. pmid:20616070

8. Köpke M, Mihalcea C, Liew F, Tizard JH, Ali MS, Conolly JJ, et al. 2, 3-Butanediol production by acetogenic bacteria, an alternative route to chemical synthesis, using industrial waste gas. Appl Environ Microbiol. 2011;77(15):5467–75. pmid:21685168

9. Kundiyana DK, Huhnke RL, Wilkins MR. Effect of nutrient limitation and two-stage continuous fermentor design on productivities during “Clostridium ragsdalei” syngas fermentation. Bioresour Technol. 2011;102(10):6058–64. pmid:21470855

10. Pachauri RK, Allen M, Barros V, Broome J, Cramer W, Christ R, et al. Climate Change 2014: Synthesis Report. Contribution of Working Groups I, II and III to the Fifth Assessment Report of the Intergovernmental Panel on Climate Change. 2014. Available: http://epic.awi.de/37530/.

11. Cord-Ruwisch R, Seitz H-J, Conrad R. The capacity of hydrogenotrophic anaerobic bacteria to compete for traces of hydrogen depends on the redox potential of the terminal electron acceptor. Arch Microbiol. 1988;149(4):350–7.

12. Gibson G, Cummings J, Macfarlane G, Allison C, Segal I, Vorster H, et al. Alternative pathways for hydrogen disposal during fermentation in the human colon. Gut. 1990;31(6):679–83. pmid:2379871

13. Joblin K. Ruminal acetogens and their potential to lower ruminant methane emissions. Crop Pasture Sci. 1999;50(8):1307–14.

14. Dong H, Mangino J, McAllister T, Have D. Emissions from livestock and manure management. Climate Change 2006: IPCC Guidelines for National Greenhouse Gas Inventories. 2006. Available: http://www.citeulike.org/group/10326/article/4425199.

15. Franz R, Soliva CR, Kreuzer M, Hummel J, Clauss M. Methane output of rabbits (Oryctolagus cuniculus) and guinea pigs (Cavia porcellus) fed a hay-only diet: Implications for the scaling of methane production with body mass in non-ruminant mammalian herbivores. Comp Biochem Physiol A Mol Integr Physiol. 2011;158(1):177–81. pmid:20971203

16. Belenguer A, Fondevila M, Balcells J, Abecia L, Lachica M, Carro M. Methanogenesis in rabbit caecum as affected by the fermentation pattern: in vitro and in vivo measurements. World Rabbit Sci. 2011;19(2):75–83.

17. Martinez G, Abecia L, Martin-Garcia A, Ramos-Morales E, Molina-Alcaide E, Ranilla M, et al. In vitro evaluation of some plant extracts as methane-inhibiting agents with diets with different degradability using rumen liquor from goats. Adv Anim Biosci. 2011;2:488.

18. Abecia L, Fondevila M, Rodríguez-Romero N, Martínez G, Yáñez-Ruiz D. Comparative study of fermentation and methanogen community structure in the digestive tract of goats and rabbits. J Anim Physiol Anim Nutr (Berl). 2013;97(s1):80–8.

19. Marounek M, Fievez V, Mbanzamihigo L, Demeyer D, Maertens L. Age and incubation time effects on in vitro caecal fermentation pattern in rabbits before and after weaning. Arch Anim Nutr. 1999;52(2):195–201.

20. Piattoni F, Demeyer D, Maertens L. In vitro study of the age-dependent caecal fermentation pattern and methanogenesis in young rabbits. Reprod Nutr Dev. 1996;36(3):253–61. pmid:8766730

21. Bäckhed F, Ley RE, Sonnenburg JL, Peterson DA, Gordon JI. Host-bacterial mutualism in the human intestine. Science. 2005;307(5717):1915–20. pmid:15790844

22. Ley RE, Hamady M, Lozupone C, Turnbaugh PJ, Ramey RR, Bircher JS, et al. Evolution of mammals and their gut microbes. Science. 2008;320(5883):1647–51. pmid:18497261

23. Ochman H, Worobey M, Kuo C-H, Ndjango J-BN, Peeters M, Hahn BH, et al. Evolutionary relationships of wild hominids recapitulated by gut microbial communities. PLoS Biol. 2010;8(11) e1000546. pmid:21103409

24. Schloss PD, Schubert AM, Zackular JP, Iverson KD, Young VB, Petrosino JF. Stabilization of the murine gut microbiome following weaning. Gut Microbes. 2012;3(4):383–93. pmid:22688727

25. Jami E, Israel A, Kotser A, Mizrahi I. Exploring the bovine rumen bacterial community from birth to adulthood. ISME J. 2013;7(6):1069–79. pmid:23426008

26. Gallois M, Le Huërou-Luron I, Fortun-Lamothe L, Lalles J, Gidenne T. Adaptability of the digestive function according to age at weaning in the rabbit: I. Effect on feed intake and digestive functionality. Animal. 2008;2(4):525–35. pmid:22443566

27. Association APoEAVM. 2000 Report of the AVMA Panel on Euthanasia. J Am Vet Med Assoc. 2001;218(5):669. pmid:11280396

28. Ernster L, Herlin L, Zetterström R. Experimental studies on the pathogenesis of kernicterus. Pediatrics. 1957;20(4):647–52. pmid:13484304

29. Yang B, He B, Wang SS, Liu JX, Wang JK. Early supplementation of starter pellets with alfalfa improves the performance of pre-and postweaning Hu lambs. J Anim Sci. 2015;93(10):4984–4994. pmid:26523591

30. Weatherburn M. Phenol-hypochlorite reaction for determination of ammonia. Anal Chem. 1967;39(8):971–974.

31. Hu WL, Liu JX, Ye JA, Wu YM, Guo YQ. Effect of tea saponin on rumen fermentation in vitro. Anim Feed Sci Technol. 2005;120(3):333–339.

32. Gagen EJ, Denman SE, Padmanabha J, Zadbuke S, Al Jassim R, Morrison M, et al. Functional gene analysis suggests different acetogen populations in the bovine rumen and tammar wallaby forestomach. Appl Environ Microbiol. 2010;76(23):7785–95. pmid:20889794

33. Xu K, Liu H, Du G, Chen J. Real-time PCR assays targeting formyltetrahydrofolate synthetase gene to enumerate acetogens in natural and engineered environments. Anaerobe. 2009;15(5):204–13. pmid:19328859

34. Denman SE, Tomkins NW, McSweeney CS. Quantitation and diversity analysis of ruminal methanogenic populations in response to the antimethanogenic compound bromochloromethane. FEMS Microbiol Ecol. 2007;62(3):313–22. pmid:17949432

35. Denman SE, McSweeney CS. Development of a real-time PCR assay for monitoring anaerobic fungal and cellulolytic bacterial populations within the rumen. FEMS Microbiol Ecol. 2006;58(3):572–82. pmid:17117998

36. Ritalahti KM, Amos BK, Sung Y, Wu Q, Koenigsberg SS, Löffler FE. Quantitative PCR targeting 16S rRNA and reductive dehalogenase genes simultaneously monitors multiple Dehalococcoides strains. Appl Environ Microbiol. 2006;72(4):2765–74. pmid:16597981

37. Bartram AK, Lynch MD, Stearns JC, Moreno-Hagelsieb G, Neufeld JD. Generation of multimillion-sequence 16S rRNA gene libraries from complex microbial communities by assembling paired-end Illumina reads. Appl Environ Microbiol. 2011;77(11):3846–52. pmid:21460107

38. Caporaso JG, Kuczynski J, Stombaugh J, Bittinger K, Bushman FD, Costello EK, et al. QIIME allows analysis of high-throughput community sequencing data. Nat Methods. 2010;7(5):335–6. pmid:20383131

39. Tamura K, Dudley J, Nei M, Kumar S. MEGA4: molecular evolutionary genetics analysis (MEGA) software version 4.0. Mol Biol Evol. 2007;24(8):1596–9. Epub 2007 May 7. pmid:17488738

40. Gasteiger E, Gattiker A, Hoogland C, Ivanyi I, Appel RD, Bairoch A. ExPASy: the proteomics server for in-depth protein knowledge and analysis. Nucleic Acids Res. 2003;31(13):3784–8. pmid:12824418

41. Li W, Godzik A. Cd-hit: a fast program for clustering and comparing large sets of protein or nucleotide sequences. Bioinformatics. 2006;22(13):1658–9. pmid:16731699

42. McMurdie PJ, Holmes S. Phyloseq: an R package for reproducible interactive analysis and graphics of microbiome census data. PLoS One. 2013;22:8(4):e61217. pmid:23630581

43. Anders S, Huber W. Differential expression analysis for sequence count data. Genome Biol. 2010;11(10):R106. pmid:20979621

44. Wickham H. ggplot2: elegant graphics for data analysis: Springer Science & Business Media. 2009.

45. Warnes GR, Bolker B, Bonebakker L, Gentleman R, Huber W, Liaw A, et al. gplots: Various R programming tools for plotting data. R package version 2(4). 2009.

46. Stamatakis A. RAxML-VI-HPC: maximum likelihood-based phylogenetic analyses with thousands of taxa and mixed models. Bioinformatics. 2006;22(21):2688–90. pmid:16928733

47. Ludwig W, Strunk O, Westram R, Richter L, Meier H, Buchner A, et al. ARB: a software environment for sequence data. Nucleic Acids Res. 2004;32(4):1363–71. pmid:14985472

48. Bovera F, Calabrò S, Cutrignelli MI, Infascelli F, Piccolo G, Nizza S, et al. Prediction of rabbit caecal fermentation characteristics from faeces by in vitro gas production technique: roughages. J Anim Physiol Anim Nutr (Berl). 2008;92(3):260–71.

49. Abecia L, Fondevila M, Balcells J, Edwards JE, Newbold CJ, McEwan NR. Molecular profiling of bacterial species in the rabbit caecum. FEMS Microbiol Lett. 2005;244(1):111–5. pmid:15727829

50. Palmer C, Bik EM, DiGiulio DB, Relman DA, Brown PO. Development of the human infant intestinal microbiota. PLoS Biol. 2007;5(7):e177. pmid:17594176

51. Yatsunenko T, Rey FE, Manary MJ, Trehan I, Dominguez-Bello MG, Contreras M, et al. Human gut microbiome viewed across age and geography. Nature. 2012;486(7402):222–7. pmid:22699611

52. Combes S, Michelland RJ, Monteils V, Cauquil L, Soulié V, Tran NU, et al. Postnatal development of the rabbit caecal microbiota composition and activity. FEMS Microbiol Ecol. 2011;77(3):680–9. pmid:21658088

53. Eshar D, Weese JS. Molecular analysis of the microbiota in hard feces from healthy rabbits (Oryctolagus cuniculus) medicated with long term oral meloxicam. BMC Vet Res. 2014;10:62 pmid:24618207

54. Godwin S, Kang A, Gulino L-M, Manefield M, Gutierrez-Zamora M-L, Kienzle M, et al. Investigation of the microbial metabolism of carbon dioxide and hydrogen in the kangaroo foregut by stable isotope probing. ISME J. 2014;8(9):1855–65. pmid:24621520

55. Breznak JA, Blum JS. Mixotrophy in the termite gut acetogen, Sporomusa termitida. Arch Microbiol. 1991;156(2):105–10.

56. Forsythe S, Parker D. Nitrogen metabolism by the microbial flora of the rabbit caecum. J Appl Bacteriol. 1985;58(4):363–9. pmid:3997689

57. Gouet P, Fonty G, editors. Evolution de la microflore digestive du lapin holoxénique de la naissance au sevrage. Ann Biol Anim Biochim Biophys. 1973;13(4):733–735.

58. Leaphart AB, Lovell CR. Recovery and analysis of formyltetrahydrofolate synthetase gene sequences from natural populations of acetogenic bacteria. Appl Environ Microbiol. 2001;67(3):1392–5. pmid:11229939

59. Bernalier A, Willems A, Leclerc M, Rochet V, Collins MD. Ruminococcus hydrogenotrophicus sp. nov., a new H2/CO2-utilizing acetogenic bacterium isolated from human feces. Arch Microbiol. 1996;166(3):176–83. pmid:8703194

60. Rieu-Lesme F, Morvan B, Collins M, Fonty G, Willems A. A new H2/CO2-using acetogenic bacterium from the rumen: Description of Ruminococcus schinkii sp. nov. FEMS Microbiol Lett. 1996;140(23):281–6.

61. Lorowitz WH, Bryant MP. Peptostreptococcus productus strain that grows rapidly with CO as the energy source. Appl Environ Microbiol. 1984;47(5):961–4. pmid:6430231

62. Geerligs G, Aldrich H, Harder W, Diekert G. Isolation and characterization of a carbon monoxide utilizing strain of the acetogen Peptostreptococcus productus. Arch Microbiol. 1987;148(4):305–313.

63. Blackmore TM, Dugdale A, Argo CM, Curtis G, Pinloche E, Harris PA, et al. Strong stability and host specific bacterial community in faeces of ponies. PLoS One. 2013;8(9):e75079. pmid:24040388

64. Zilber-Rosenberg I, Rosenberg E. Role of microorganisms in the evolution of animals and plants: the hologenome theory of evolution. FEMS Microbiol Rev. 2008;32(5):723–35. pmid:18549407

65. Breznak JA, Kane MD. Microbial H2/CO2 acetogenesis in animal guts: nature and nutritional significance. FEMS Microbiol Rev. 1990;7(3–4):309–13. pmid:2128799

66. Kušar D, Avguštin G. Molecular profiling and identification of methanogenic archaeal species from rabbit caecum. FEMS Microbiol Ecol. 2010;74(3):623–30. pmid:20950344