Research article | Open | Published:
Molecular analysis of the microbiota in hard feces from healthy rabbits (Oryctolagus cuniculus) medicated with long term oral meloxicam
BMC Veterinary Researchvolume 10, Article number: 62 (2014)
Analgesia is often indicated in rabbits undergoing surgical procedures or suffering from various painful conditions and the most common adverse effects associated with NSAIDs occur in the gastrointestinal tract (GIT). The objective of this study was to determine the potential effect of long-term (21 days) meloxicam administration on the fecal bacterial microbiota in healthy rabbits.
Samples of hard feces were collected from six rabbits treated with meloxicam (1 mg/kg orally once every 24 h) on days 0,6,14 and 21. Next generation sequencing of V4 16S rRNA gene products was performed.
A total of 2589912 V4 rRNA gene sequences passed all quality control filters. Firmicutes predominated (82.0 ± 6.2%). Sixteen other phyla were also identified but other than Verrucomicrobia (4.4 ± 4.9%), all accounted for less than 1% of the identified sequences. Within Firmicutes, Clostridia was the dominant class, accounting for 76% of operational taxon units (OTUs). In general, there were only few differences observed between time points and different rabbits at the phylum level. A significant change was observed in the relative abundance of Proteobacteria over the 4 time points (P = 0.02).
The gastrointestinal tract of rabbits harbors dense and diverse microbiota. Significant alteration of the hard fecal microbiota does not appear to be a considerable adverse effect expected in rabbits treated for 21 days with oral meloxicam at a dose of 1 mg/kg.
Analgesia is a basic aspect of veterinary medicine, indicated in both laboratory and pet-owned animals undergoing surgical procedures or suffering from various painful conditions . One of the most common classes of drugs used by veterinarians for analgesia is non-steroidal anti-inflammatory drugs (NSAIDs) . Nonsteroidal anti-inflammatory drugs reduce inflammation by inhibiting the action of Cyclooxygenase (COX) enzymes, which convert arachidonic acid into prostanoids and analgesia associated with NSAIDs is mainly due to inhibition of COX-2 . The most common adverse effects associated with NSAIDs occur in the gastrointestinal tract (GIT), and to a lesser degree, impairment of renal blood flow [3, 4]. Reported NSAIDs-related GIT pathologies include perforation, ulceration, bleeding and local irritation [3, 4]. Because prostaglandins also play a role in the rabbit’s GIT motility , NSAIDs have the potential to also impair GIT function. To date, several studies have been performed to determine pharmacokinetics and safety of meloxicam in rabbits but none had looked into the potential drug’s effect on GIT microbiota [6–8].
The fecal or intestinal microbiota of domestic rabbits has been extensively studied through culture dependent methods, as well as earlier generation culture-independent methods, with a focus on the cecal microbiota [9–13]. However, limitations of these methods for understanding the overall microbial composition are well known. In particular, culture-dependent studies are only able to provide a very superficial (and sometimes misleading) assessment of this complex polymicrobial environment [14–16]. While earlier studies have provided some important information, the ability of next generation sequencing based studies such as this to evaluate millions of bacterial sequences from diverse (and often unculturable) microbial populations is now possible.
This study explores the stability over time of bacterial communities in hard feces in healthy rabbits medicated with oral meloxicam, using high throughput sequencing of the V4 region of the 16S rRNA gene.
Ethical approval was obtained from the Institutional Animal Care and Use Committee of Kansas State University. Rabbits enrolled in this study were part of a broader study evaluating long term (29 days) administration of oral meloxicam. Six clinically normal 3-month-old intact female New Zealand White rabbits (body weight range, 2.52 to 2.71 kg) were included in this study. The rabbits were obtained from a commercial source and were specific pathogen (Pasteurella sp.) free. The rabbits were housed individually at the research facilities of the Kansas State University College of Veterinary Medicine at a constant temperature (21°C) and humidity (60%), and were exposed to cycles of 16 hours of light and 8 hours of dark per day. Rabbits were fed free choice of timothy haya, timothy-based pelleted dietb and water was available ad libitum. Rabbits were acclimated to the facility for 5 days after their arrival and were habituated to handling prior to initiation of the study. Immediately prior to the start of the study, each rabbit was determined to be behaviorally normal and healthy based on a thorough physical examination and clinical pathology testing (complete blood count [CBC], serum biochemistry profile and urinalysis), with results compared to published reference ranges . Physical examination, weight and blood testing were repeated weekly in all rabbits enrolled in the study. Fecal production was subjectively monitored for shape, color and consistency, and objectively by weight and fecal pellet size. At the end of the study, all rabbits were humanely euthanatized using pentobarbitalc (100 mg/kg Intravenous), and necropsy was performed within 15 minutes of euthanasia.
Experimental design and sample collection
Meloxicamd (1 mg/kg) was orally administered daily for 29 days to each rabbit via a 3-mL syringe. The rabbits’ behaviour, attitude, mentation, level of activity, eating and drinking, and defecation were subjectively assessed daily by the researchers. Hard fecal samples were collected prior to treatment (day 0) and on days 6, 14 and 21 from the start of the meloxicam treatment. At each time point morning (08:00), 10 fecal pellets from a fresh elimination were placed into an individual collection bag. Therefore, four fecal samples were collected from each rabbit during the study period for determination of fecal microbiota. Each sample (10 fecal pellets) was weighed using a digital gram scale and each fecal pellet was measured using hand calipers. All bagged samples were stored at −70°C until analysis.
DNA extraction and quality control
DNA extraction was performed using a commercial kite following the manufacturer’s “stool DNA protocol for pathogen detection”, which includes an initial bead beating step and proteinase K treatment. DNA quantity and quality were accessed by spectrophotometry using the NanoDropf.
16 s rRNA gene amplification and sequencing
PCR amplification of the V4 region of the 16S rRNA gene was designed based on Klindworth et al.  using the primers forward S-D-Bact-0564-a-S-15 (5′-AYTGGGYDTAAAGNG-3′) and reverse S-D-Bact-0785-b-A-18 (5′-TACNVGGGTATCTAATCC-3′). The forward and reverse primers were designed containing an overlapping region of the forward and reverse Illumina sequencing primers (TCGTCGGCAGCGTCAGATGTGTATAAGAGACAG and GTCTCGTGGGCTCGGAGATGTGTATAAGAGACAG, respectively) in order to anneal them to primers containing the Illumina adaptors plus the 8 bp identifier indices (forward: AATGATACGGCGACCACCGAGATCTACAC-index-TCGTCGGCAGCGTC; reverse: CAAGCAGAAGACGGCATACGAGAT-index-GTCTCGTGGGCTCGG). For a final volume of 50 μL, 2 μL of each DNA sample were added to a solution containing 18.7 μL of water, 25 μL of ReadMisg, 1.3 μL of BSAh, 0.5 μL of each 16S primer (1000 pmol/μL) and 1 μL of each Illumina primer (1000 pmol/μL). The mixture was subjected to the following PCR conditions, 3 minutes at 94˚C for denaturing, and 35 cycles of 45 seconds at 94˚C for denaturing, 60 seconds at 50˚C for annealing and 90 seconds at 72˚C for elongation followed by a final period of 10 minutes at 72˚C and kept at 4˚C until purification.
PCR products were evaluated by electrophoresis in 2% agarose gel and purified with the Agencourt AMPure XPi by mixing 22 μL of amplicon with 72 μL of AMPure on a 96 well plate. After 5 minutes at room temperature, beads were separated and washed twice with 80% ethanol and eluted in 30 μL of water. After purification samples were quantified by spectrophotometry using the NanoDropj and normalized to a final concentration of 2 nM. Sequencing of the library pool was performed at the University of Guelph’s Advanced Analysis Centre using an Illumina MiSeqk and 2x250 chemistry.
The MOTHUR package of algorithms (v1.30) was used for analysis . Paired end reads were aligned. Sequences >275 bp and those containing ambiguous base calls or long runs (>8 bp) of holopolymers were removed, as were sequences that did not align with the correct region. Chimeras were detected using uchime  and removed. Sequences from chloroplasts, mitochondria, Archaea and eukaryotes were also removed.
Sequences were binned into OTUs at a 3% (0.03) dissimilarity level, and taxonomy was assigned using the Silva 16S rRNA reference database .
Coverage was assessed using Good’s coverage. Population diversity was described using the inverse Simpson’s index. Relative abundances were compared between animals and between timepoints using ANOVA.
Dissimilarity between the two groups was assessed through creation of dendrograms using the Yue & Clayton measure of dissimilarity (a measure of community structure, which considers shared OTUs and their relative abundances) and traditional Jaccard index (a measure of community membership, which considers the number of shared OTUs, not their abundance). Parsimony and unweighted unifrac  tests were applied to evaluate the impact of day of treatment on population structure. Principal coordinate analysis (PCoA) was also performed, based on both Jaccard and Yue and Clayton indices. Analysis of molecular variance (AMOVA), a non-parametric evaluation of genetic diversity, was used to evaluate the impact of day.
Meloxicam was easily administered daily to each rabbit, and all of the rabbits seemed to remain healthy during the study. Subjectively, none of the rabbits had changes in behavior, attitude, mentation, level of activity, amount of food or water consumed, or fecal production. Objectively, all participating rabbits continued to gain weight (final body weight range 2.61-2.91 Kg) when compared to the beginning of the study period. Also, all values in the weekly serum biochemistry analysis were within normal limits both before and during the study. No changes were observed in hard feces’ shape, color and consistency (no diarrhea or ileus). The mean weight of a 10-pellet fecal sample was 6.75 g (range 6–9 g) and mean measured individual fecal pellet size was 6.2 mm (range 5–8 mm). Post-study gross necropsy and histopathologic evaluation of all participating rabbits did not reveal any major abnormalities.
A total of 3453182 V4 16S RNA gene sequences were obtained; 2589912 of which passed all quality control filters. The number of sequences per sample ranged from 35178 to 250720 (mean 107913, SD 42703). Further analysis was performed on random subsampling of 35178 sequences per sample. Excellent sample coverage was obtained with this subsampled population, as demonstrated by a mean Good’s coverage value of 0.972 (SD 0.002, range 0.967-0.975) and rarefaction curves (Figure 1). The number of OTUs ranged from 1814–2525 (mean 2118, SD 170). There was no impact of rabbit (P = 0.16) or day (P = 0.16) on the number of different OTUs that were identified. The overall mean inverse Simpson’s value was 73.6 (SD 23.6, range 19.5-119), with no impact of day (P = 0.54) or rabbit (P = 0.07).
The microbiota of all rabbits was dominated at the phylum level by Firmicutes, which accounted for 66.6 ± 6.7% (mean ± SD) of sequences (Figure 2). Sixteen other phyla were identified; however, only Verrucomicrobia (4.4 ± 5.1%) and accounted for greater than 1% of sequences. These other phyla consisted of Actinobacteria (0.18 ± 0.13%), Bacteroidetes (0.38 ± 0.26%), Chlamydiae (0.001 ± 0.00%), Cyanobacteria (0.003 ± 0.003), Deinococcus-Thermus (0.0001 ± 0.0005%), Fibrobacteres (0.021 ± 0.014%), Fusobacteria (0.003 ± 0.002%), Lentisphaerae (0.002 ± 0.002%, Plantomycetes (0.0002 ± 0.0008), Proteobacteria (0.57 ± 0.33%), Spirochaetes (0.043 ± 0.24%), Synergistetes (0.005 ± 0.001%), Tenericutes (0.021 ± 0.017%), SR1 (0.0006 ± 0.001%) and TM7 (0.06 ± 0.07%). Sequences that were unclassified at the phylum level accounted for 27.7 ± 5.5% of sequences.
Within Firmicutes, Clostridia was the dominant class, accounting for 76% of OTUs from Firmicutes, with Ruminococcaceae and Lachnospiraceae being the main families (order, Clostridiales), representing 48 and 20% of sequences from Clostridia, respectively. Unclassified Clostridia were common (24%), followed by smaller numbers of Bacilli (2.9%).
Within Verrucomicrobia, Verrucomicrobiaceae was the predominant (98%) family. Within that family, Akkermansia and Persicirhabdus were the two identified genera, but they only accounted for 9.9 and 0.2% of sequences from Verrucomicrobia, respectively. The remaining members of that family were unclassified at the genus level.
There were few differences between different timepoints and different rabbits at the phylum level. Rabbit 4 had a significantly higher relative abundance of Verrucomicrobia (12%) compared to the others (2.1-3.5%) (P = 0.03). No other differences between rabbits were evident at the phylum level. There was a significant change in the relative abundance of Proteobacteria over the four timepoints (P = 0.02), with day 21 having a significant higher relative abundance than the other three time points. Otherwise, no significant differences between timepoints were identified.
Dendrograms based on Jaccard and Yue & Clayton indices are presented in Figures 3 and 4. There were few differences in microbial population structure. No effect of sampling day was evident by unweighted unifrac (Table 1). A significant difference was noted with parsimony test but only between days 0 and 21 and using the Jaccard index.
Principle coordinate analysis results are presented in Figures 5 and 6. Using AMOVA, there was no difference between days using the Yue & Clayton index (P = 0.082), but there was a difference using the Jaccard index (P = 0.011), with significant differences between days 0 and 21 (P = 0.002) and days 14 and 21 (P = 0.030).
Rabbits enrolled in this study did not show any abnormal clinical signs, hematological values, postmortem evaluation, or physical changes of the hard feces throughout the observed study period, thus suggesting that the tested meloxicam dose can be considered safe for clinical use.
This study identified marked (but perhaps unsurprising) richness and diversity of the fecal microbiome that is incomparable with culture-dependent studies, as an average of over 2000 OTUs were identified. This and the corresponding diversity indices highlight the complexity of the rabbit fecal microbiota, even in rabbits housed under controlled conditions with little or no variation in diet and the need to use broad approaches such as this to understand this environment.
The predominance of Firmicutes was not particularly surprising, as this has been reported in other culture-independent studies of the rabbit  and can be found in other hindgut fermenters such as the horse [23, 24]. It was interesting to observe that Ruminococcaceae and Lachnospiraceae were leading members of this phylum. Those two families are receiving increasing attention in various species as potentially key components of the ‘beneficial’ microbiota [23, 25, 26]. There is limited information about these bacteria in rabbits, although they were also reported to be dominant families in a study of 9-week old rabbits . These may be critical components of the rabbit microbiome and groups to study further in terms of gut health and the loosely defined ‘intestinal dysbiosis’ that is often diagnosed in rabbits with gastrointestinal disease .
While the predominance of Firmicutes was not unexpected, the utter dominance of this phylum was remarkable. In studies of the fecal microbiome of most species, while one phylum may be most common, there are typically multiple phyla with relative abundances greater than 5% (or higher) [24, 29–31]. Here, only one phylum beyond Firmicutes accounted for greater than 1% of sequences. That was Verrucomicrobia, which accounted for an average of 4.4% of sequences. Verrucomicrobia is a relatively newly described and poorly characterized phylum. It has been reported in similar relative abundances in other species, including horses [24, 32, 33], and has been shown to increase in abundance in mice fed type 2 resistant starch , perhaps indicating an important role in hindgut fermentation. Akkermansia (a member of this phylum) has been shown to be decreased in obese versus lean children and mice [35, 36], to be associated with improved glucose homeostasis in metformin treated mice , to be associated with reversal of high fat diet induced metabolic disorders in mice  and to be reduced in humans with ulcerative colitis . These suggest that this phylum may play an important role in gastrointestinal health in rabbits.
The number of sequences that were unclassified at the phylum level was large, although within the range that has been reported in other fecal microbiome studies and lower than that reported in earlier cloning-based studies of the rabbit cecal microbiota [10, 39]. Even within sequences classified as members of the main phyla, Firmicutes and Verrucomicrobia, there were large percentages of sequences that were not identified to lower taxonomic orders. The relevance of these unclassified organisms is unclear. Unclassified sequences are not necessarily unknown phyla. Unclassified could indicate a sequence that does not have any closely related sequences in the reference database, but it could also indicate an inability of the analysis to assign a sequence to known phyla (or lower taxonomic levels) because of similarities in sequences or because of limitations in discriminatory power. Yet, unknown phyla could be present, particularly because of the paucity of knowledge regarding the fecal microbiota. The relevance of these in health and disease is completely unknown but given the high relative abundance, it must be assumed that this unclassified component of the microbiota could be of importance. This underscores our limited understanding of this complex polymicrobial ecosystem and the need for further study to define the microbiota and its individual constituents.
When assessing the microbial composition of these rabbits, there was surprisingly little clustering by individual, as is evident from the dendrograms and PCoA plot. With studies of other species, it is typically reported that there is less intra-individual variation than inter-individual variation [40, 41], and if that was the case here, clustering by animal would have been expected. There are a few potential explanations. It is possible that there was enough of an effect of treatment to alter the microbiota of each rabbit, thereby disrupting intra-individual consistency, but that the alternation was not in a consistent or strong enough manner to result in a statistically significant effect of sampling day. Another more likely explanation might be the nature of this rabbit population, with similar genetics, housing and diet that resulted in a very homogenous microbiota between individuals. Study of variation in the microbiota of rabbits outside of a research colony would be enlightening. It is also possible that there is more inherent variation in the rabbit microbiota, since a study of rabbits reported similar intra- and inter-rabbit variation . However, that study also involved a homogenously managed group of research rabbits so care must be taken when interpreting the degree of inter- and intra- individual variation from these studies.
Little impact of treatment was identified and there were inconsistent results between different tests. Differences were noted between days 0 and 21 using parsimony test and AMOVA, as well as days 14 and 21 using AMOVA. However, these were only when analysis was based on the Jaccard index. It is important to consider differences between these indices. The Yue & Clayton index is a measure of population dissimilarity based on both microbial presence and abundance, while the Jaccard index compares community membership without assessment of the relative abundance of different members of the community. Thus, differences in Jaccard but not Yue & Clayton indices suggest that changes in the community membership (addition or loss of members) occurred, but that this may have mainly involved low abundance members of community based on a lack of difference when relative abundance is considered. Other than the potential long-term microbiota adjustment to a new diet  the accumulated effect of meloxicam should also be considered, and in the absence of clinical or histologic abnormalities, only advanced molecular testing methods had the sensitivity to detect subtle changes in the rabbit’s GIT microbiota.
Limitations of the study must be considered. Since there was no untreated control group raised in parallel with these animals and they were all treated at the same time, it is possible that there was some unidentified factor such as a subtle change in diet or stress that accounted for the microbiota change. Additionally, it is known that the fecal microbiome of laboratory raised rabbits is different from rabbits raised under other conditions (e.g. farms, households) [42–44]. As with similar studies, there is always the question of how well the fecal microbiome represents the microbiome in more proximal sections of the gastrointestinal tract. That is an inherent limitation given the inability to easily and non-invasively obtain samples from different areas of the gastrointestinal tract over time. It has been reported that hard feces (which were studied here) are less representative of the cecal microbiota than soft feces  but conflicting data have also been reported [15, 16]. These data should therefore be interpreted as a description and comparison of the hard fecal microbiota. Optimal approaches for studying the proximal intestinal tract in rabbits are not known. The sample size may also have limited statistical power; however, given the depth of data that are obtained and the similar sample sizes that have been used in other fecal microbiome studies, the sample size was deemed adequate to identify significant changes in the microbiota.
The fecal microbiome of the rabbit is a complex and diverse polymicrobial ecosystem, composed of a wide range of both well-understood and completely unknown bacteria. While there were some noted differences over time, there was no clear and consistent impact of meloxicam treatment on the fecal microbiome. Given the importance of the intestinal microbiome in health and disease, particularly amongst hindgut fermenting species such as the rabbit, it is clear that there is inadequate understanding of the ‘normal’ microbiota, how it varies within and between rabbits, and how this complex polymicrobial ecosystem influences (and is influenced by) disease, nutrition, management and myriad other factors.
In the absence of adverse clinical, histologic and significant microbial changes observed in rabbits enrolled in this study, it is suggested that alteration of the hard fecal microbiota is not a significant adverse reaction expected with daily oral meloxicam administrated at the tested dose.
Availability of supporting data
The data set supporting the results of this article are available at the MG-RAST metagenomics analysis server (project 8102, http://metagenomics.anl.gov).
aWestern Timothy Hay, Oxbow Animal Health, Murdock, Neb.
bEssentials – Young Rabbit, Oxbow Animal Health, Murdock, Neb.
cVortech Pharmaceuticals, Dearborne, MO.
dMetacam 1.5 mg/ml oral suspension, Boehringer Ingelheim Vetmedica, St Joseph, MO.
eE.Z.N.A. Stool DNA Kit, Omega Bio-Tek Inc., Doraville, GA.
fRoche, Mississauga, ON, Canada.
gInvitrogen, San Diego, CA.
hInvitrogen, San Diego, CA.
iBeckman Coulter Inc, Mississauga, ON, Canada.
jRoche, Mississauga, ON, Canada.
kIllumina Inc, San Diego, CA.
Wenger S: Anesthesia and analgesia in rabbits and rodents. J Exotic Pet Med. 2012, 21: 7-16. 10.1053/j.jepm.2011.11.010.
Streppa HK, Jones CJ, Budsberg SC: Cyclooxygenase selectivity of nonsteroidal anti-inflammatory drugs in canine blood. Am J Vet Res. 2002, 63: 91-94. 10.2460/AJVR.2002.63.91.
Forsyth SF, Guilford WG, Haslett SJ, Hill PB: Endoscopy of the gastroduodenal mucosa after carprofen, meloxicam and ketoprofen administration in dogs. J Small Anim Pract. 1998, 39: 421-424. 10.1111/j.1748-5827.1998.tb03748.x.
Luna SPL, Basilio AC, Steagall PVM, Machado LP, Moutinho FQ, Takahira RK, Brandao CV: Evaluation of adverse effects of long-term oral administration of carprofen, etodolac, flunixin meglumine, ketoprofen, and meloxicam in dogs. Am J Vet Res. 2007, 68: 258-264. 10.2460/ajvr.68.3.258.
Pairet M, Bouyssou T, Ruckebusch Y: Colonic formation of soft feces in rabbits: a role for endogenous prostaglandins. Am J Physiol-Gastrointestinal Liver Physiol. 1986, 250: G302-G308.
Turner PV, Chen HC, Taylor MW: Pharmacokinetics of meloxicam in rabbits after single and repeat oral dosing. Comp Med. 2006, 56: 63-67.
Carpenter JW, Pollock CG, Koch DE, Hunter RP: Single and multiple-dose pharmacokinetics of meloxicam after oral administration to the rabbit (Oryctolagus cuniculus). J Zoo Wildl Med. 2009, 40: 601-606. 10.1638/2007-0115.1.
Fredholm DV, Carpenter JW, KuKanich B, Kohles M: Pharmacokinetics of meloxicam in rabbits after oral administration of single and multiple doses. Am J Vet Res. 2013, 74: 636-641. 10.2460/ajvr.74.4.636.
Straw TE: Bacteria of the rabbit gut and their role in the health of the rabbit. J Appl Rabbit Res. 1988, 11: 142-150.
Abecia L, Fondevila M, Balcells J, Edwards JE, Newbold CJ, McEwan RN: Molecular profiling of bacterial species in the rabbit caecum. FEMS Microbiol Lett. 2005, 244: 111-115. 10.1016/j.femsle.2005.01.028.
Furet JP, Firmesse O, Gourmelon M, Bridonneau C, Tap J, Mondot S, Dore J, Corthier G: Comparative assessment of human and farm animal faecal microbiota using real‒time quantitative PCR. FEMS Microbiol Ecol. 2009, 68: 351-362. 10.1111/j.1574-6941.2009.00671.x.
Kušar D, Avguštin G: Molecular profiling and identification of methanogenic archaeal species from rabbit caecum. FEMS Microbiol Ecol. 2010, 74: 623-630. 10.1111/j.1574-6941.2010.00980.x.
Combes S, Michelland RJ, Monteils V, Cauquil L, Soulie V, Tran NU, Gidenne T, Fortun-Lamothe L: Postnatal development of the rabbit caecal microbiota composition and activity. FEMS Microbiol Ecol. 2011, 77: 680-689. 10.1111/j.1574-6941.2011.01148.x.
Suau A, Bonnet R, Sutren M, Godon JJ, Gibson GR, Collins MD, Dore J: Direct analysis of genes encoding 16S rRNA from complex communities reveals many novel molecular species within the human gut. Appl Environ Microbiol. 1999, 65: 4799-4807.
Michelland RJ, Combes S, Monteils V, Cauquil L, Gidenne T, Fortun-Lamothe L: Molecular analysis of the bacterial community in digestive tract of rabbit. Anaerobe. 2010, 16: 61-65. 10.1016/j.anaerobe.2009.05.002.
Michelland RJ, Monteils V, Combes S, Cauquil L, Gidenne T, Fortun-Lamothe L: Comparison of the archaeal community in the fermentative compartment and faeces of the cow and the rabbit. Anaerobe. 2010, 16: 396-401. 10.1016/j.anaerobe.2010.04.004.
Fiorello CV, Divers SJ: Rabbits. Exotic Animal Formulary. 4th edition. Edited by: Carpenter JW. St Louis: Elsevier; 2012:517-559.
Klindworth A, Pruesse E, Schweer T, Peplies J, Quast C, Horn M, Glockner FO: Evaluation of general 16S ribosomal RNA gene PCR primers for classical and next-generation sequencing-based diversity studies. Nucleic Acids Res. 2013, 41: e1-10.1093/nar/gks808.
Schloss PD, Westcott SL, Ryabin T, Hall JR, Hartmann M, Hollister EB, Lesniewski RA, Oakley BB, Parks DH, Robinson CJ, Sahl JW, Stres B, Thalinger GG, Van Horn DJ, Weber CF: Introducing mothur: open-source, platform-independent, community-supported software for describing and comparing microbial communities. Appl Environ Microbiol. 2009, 75: 7537-7541. 10.1128/AEM.01541-09.
Edgar RC, Haas BJ, Clemente JC, Quince C, Knight R: UCHIME improves sensitivity and speed of chimera detection. Bioinformatics. 2011, 27: 2194-2200. 10.1093/bioinformatics/btr381.
Quast C, Pruesse E, Yilma P, Gerken J, Schweer T, Yarza P, Peplies J, Glockner FO: The SILVA ribosomal RNA gene database project: improved data processing and web-based tools. Nucleic Acids Res. 2013, 41: D590-D596. 10.1093/nar/gks1219.
Lozupone C, Hamady M, Knight R: UniFrac–an online tool for comparing microbial community diversity in a phylogenetic context. BMC Bioinforma. 2006, 7: 371-10.1186/1471-2105-7-371.
Costa MC, Arroyo LG, Allen-Vercoe E, Staempfli HR, Kim PT, Sturgeon A, Weese JS: Comparison of the fecal microbiota of healthy horses and horses with colitis by high throughput sequencing of the V3-V5 region of the 16S rRNA gene. PloS One. 2012, 7: e41484-10.1371/journal.pone.0041484.
O’Donnell MM, Harris H, Jeffery IB, Claesson MJ, Younge B, O’Toole PW, Ross RP: The core faecal bacterial microbiome of Irish Thoroughbred racehorses. Lett Appl Microbiol. 2013.
Suchodolski JS: Companion animal’s symposium: microbes and gastrointestinal health of dogs and cats. J Anim Sci. 2011, 89: 1520-1530. 10.2527/jas.2010-3377.
Petrof EO, Gloor GB, Vanner SJ, Weese JS, Carter D, Daigneault MC, Brown EM, Schroeter K, Allen-Vercoe E: Stool substitute transplant therapy for the eradication of Clostridium difficile infection: ‘RePOOPulating’ the gut. Microbiome. 2013, 1: 1-12. 10.1186/2049-2618-1-1.
Massip K, Combes S, Cauquil L, Zemb O, Gidenne T: High throughput 16S-DNA sequencing for phylogenetic affiliation of the caecal bacterial community in the rabbit - Impact of the hygiene of housing and of the intake level. Proceeding of Symposium on Gut Microbiology. Gut Microbiota: Friend or Foe?: 17–20 June 2012. Clermont-Ferrand, France, 2012;56.
Harrenstien L: Gastrointestinal diseases of pet rabbits. Sem. Avian Exotic Pet Med. 1999, 8: 83-89. 10.1016/S1055-937X(99)80040-1.
Turnbaugh PJ, Hamady M, Yatsunenko T, Cantarel BL, Duncan A, Ley RE, Sogin ML, Jones WJ, Roe BA, Affourtit JP, Egholm M, Henrissat B, Heath AC, Knight R, Gordon JI: A core gut microbiome in obese and lean twins. Nature. 2008, 457: 480-484.
Barry KA, Middelbos IS, Vester Boler BM, Dowd SE, Suchodolski JS, Henrissat B, Coutinho PM, White BA, Fahey GC, Swanson KS: Effects of dietary fiber on the feline gastrointestinal metagenome. J Proteome Res. 2012, 11: 5924-5933.
Suchodolski JS, Markel ME, Garcia-Mazcorro JF, Unterer S, Heilmann RM, Dowd SE, Kachroo P, Ivanov I, Minamoto Y, Dillman EM, Steiner JM, Cook AK, Toresson L: The fecal microbiome in dogs with acute diarrhea and idiopathic inflammatory bowel disease. PloS One. 2012, 7: e51907-10.1371/journal.pone.0051907.
Steelman SM, Chowdhary BP, Dowd S, Suchodolski JS, Janecka JE: Pyrosequencing of 16S rRNA genes in fecal samples reveals high diversity of hindgut microflora in horses and potential links to chronic laminitis. BMC Vet Res. 2012, 8: 231-10.1186/1746-6148-8-231.
Shepherd ML, Swecker WS, Jensen RV, Ponder MA: Characterization of the fecal bacteria communities of forage‒fed horses by pyrosequencing of 16S rRNA V4 gene amplicons. FEMS Microbiol Lett. 2012, 326: 62-68. 10.1111/j.1574-6968.2011.02434.x.
Tachon S, Zhou J, Keenan M, Martin R, Marco ML: The intestinal microbiota in aged mice is modulated by dietary resistant starch and correlated with improvements in host responses. FEMS Microbiol Ecol. 2013, 83: 299-309. 10.1111/j.1574-6941.2012.01475.x.
Karlsson CL, Önnerfält J, Xu J, Molin G, Ahrne S, Thorngren-Jerneck K: The microbiota of the gut in preschool children with normal and excessive body weight. Obesity. 2012, 20: 2257-2261. 10.1038/oby.2012.110.
Everard A, Belzer C, Geurts L, Ouwerkerk JP, Druart C, Bindels LB, Guiot Y, Derrien M, Muccioloi GG, Deizenne NM, de Vos WM, Cani PD: Cross-talk between Akkermansia muciniphila and intestinal epithelium controls diet-induced obesity. Proc Natl Acad Sci. 2013, 110: 9066-9071. 10.1073/pnas.1219451110.
Shin NR, Lee JC, Lee HY, Kim MS, Whon TW, Lee MS, Bae JW: An increase in the Akkermansia spp. population induced by metformin treatment improves glucose homeostasis in diet-induced obese mice. Gut Online First. 2013, 0: 1-9.
Rajilic-Stojanovic M, Shanahan F, Guarner F, de Vos WM: Phylogenetic analysis of dysbiosis in ulcerative colitis during remission. Inflamm Bowel Dis. 2013, 19: 481-488. 10.1097/MIB.0b013e31827fec6d.
Monteils V, Cauquil L, Combes S, Godon JJ, Gidenne T: Potential core species and satellite species in the bacterial community within the rabbit caecum. FEMS Microbiol Ecol. 2008, 66: 620-629. 10.1111/j.1574-6941.2008.00611.x.
Lazarevic V, Whiteson K, Hernandez D, Francois P, Schrenzel J: Study of inter-and intra-individual variations in the salivary microbiota. BMC Genomics. 2010, 11: 523-534. 10.1186/1471-2164-11-523.
Jalanka-Tuovinen J, Salonen A, Nikkilä J, Immonen O, Kekkonen R, Lahti L, Palva A, de Vos WM: Intestinal microbiota in healthy adults: temporal analysis reveals individual and common core and relation to intestinal symptoms. PloS One. 2011, 6: e23035-10.1371/journal.pone.0023035.
Rodríguez-Romero N, Abecia L, Fondevila M: Microbial ecosystem and fermentation traits in the caecum of growing rabbits given diets varying in neutral detergent soluble and insoluble fibre levels. Anaerobe. 2013, 20: 50-57.
Bennegadi N, Fonty G, Millet L, Gidenne T, Licois D: Effects of age and dietary fibre level on caecal microbial communities of conventional and specific pathogen-free rabbits. Microb Ecol Health Dis. 2003, 15: 23-32. 10.1080/08910600310015574.
Abecia L, Rodríguez-Romero N, Yañez-Ruiz DR, Fondevilla M: Biodiversity and fermentative activity of caecal microbial communities in wild and farm rabbits from Spain. Anaerobe. 2012, 18: 344-349. 10.1016/j.anaerobe.2012.04.004.
Michelland RJ, Combes S, Monteils V, Cauquil L, Gidenne T, Fortun-Lamothe L: Rapid adaptation of the bacterial community in the growing rabbit caecum after a change in dietary fibre supply. Animal. 2011, 5: 1761-1768. 10.1017/S1751731111001005.
The authors would like to thank Dr. James W. Carpenter for allowing the fecal samples collection from his research rabbits.
The authors declare that they have no competing interests.
DE participated in fecal collection, data generation and interpretation, and manuscript preparation. JSW participated in data analysis and interpretation, and manuscript preparation. Both authors read and approved the final manuscript.