Gut microbiota of the scimitar-horned oryx, Oryx dammah.
There are three oryx species (Bovidae) in northern Africa: scimitar-horned oryx (Oryx dammah), Cape oryx (O. gazella) and Arabian oryx (O. leucoryx). Amongst these, the scimitar-horned oryx is listed as Extinct in the Wild in the International Union for Conservation of Nature Red List because there has been no confirmed evidence of the survival of this species since the early 1990s despite extensive surveys dedicated to the detection of Sahelo-Saharan antelopes conducted in Chad and Niger. This species is well adapted for survival in dry areas and is a social bovid.
Herbivores prefer to consume fresh shoots because these usually contain low concentrations of secondary plant compounds and are protein-rich, making them more palatable, easier to digest and more nourishing (Villalba et al. 2002, 2004). Wild scimitar-horned oryx inhabited steppe and desert ecosystems, where they consumed foliage, grass, herbs, shrubs, succulent plants, legumes, juicy roots, buds and fruits. However, replicating their natural diet in captivity is challenging as food items vary greatly in terms of availability, palatability and nutrient content throughout the year. Nonetheless, captive diets are often designed to mimic the natural diet of animals.
Gut microbiota plays a vital role in nutrient uptake, food digestion, energy harvesting and vitamin synthesis in humans and other animals (Hooper et al. 2002, Backhed et al. 2005). Previous studies have revealed that Firmicutes and Bacteroidetes bacteria are predominant in the rumen and faeces of Bovidae (Shanks et al. 2011, Singh et al. 2012, Thoetkiattikul et al. 2013, Xue et al. 2016). The documentation of concurrent shifts in gut microbiota will aid the establishment of baseline rates of changes in the composition of gut microbial communities across different species and diets as well as the assessment of the relationship of such gut microbial shifts with species and diet. The composition, structure and role of the gut microbiota of the scimitar-horned oryx remain largely unknown. In the present study, next-generation sequencing targeting the V3-V4 region of the bacterial 16S rRNA gene was employed to characterise the gut microbiota in faecal samples of captive scimitar-horned oryx.
Material and Methods
The study was approved by the Binzhou University Animal Care and Use Committee. Fresh faecal samples (~150 g) of five healthy scimitar-horned oryx individuals at the Dongying Wild Animal Park, China, were collected in the late morning during March 2018 (Table 1). No animal had received anti-infammatory or anti-microbial drugs within the past four months, and none suffered from gastrointestinal diseases. The samples were collected from the ground immediately after defecation and instantly transferred into separate plastic tubes. The samples were stored at 4 [degrees]C, with long-term storage at -80 [degrees]C, for further analysis. DNA extraction, 16S rRNA gene PCR and sequencing were performed. Briefly, total genomic DNA was isolated from the faecal samples using a QIAamp DNA Stool Mini Kit (Qiagen, Germany) according to the manufacturer's protocol. A UV-vis spectrophotometer (NanoDrop 2000c, U.S.A.) was used to determine DNA concentration and quality. Only samples with DNA concentrations of > 20 ng/[micro]l were selected for PCR amplification of the V3-V4 hypervariable region of the 16S rRNA gene using the primers 341F (CCTAYGGGRBGCASCAG) and 806R (GGACTACNNGGGTATCTAAT) (Barfod et al. 2013). Amplification was performed using a final volume of 30 [micro]l containing 15 [micro]l of Phusion Master Mix (2x), 1.5 [micro]l of each primer and 10 [micro]l of microbial genomic DNA. PCR conditions were as follows: denaturation for 1 min at 98 [degrees]C, followed by denaturation with 35 cycles of 10 s at 98 [degrees]C, annealing for 30 s at 55 [degrees]C and elongation for 30 s at 72 [degrees]C and final extension for 5 min at 72 [degrees]C. The PCR products were stored at 4 [degrees]C until further analysis (within two hours) using electrophoresis in 2 % agarose gel and purification using the Qiagen Gel Extraction Kit (Qiagen, Germany). Sequence libraries were constructed using TruSeq[R] DNA PCR-Free Sample Preparation Kit (Illumina, U.S.A.) according to the manufacturer's protocol, and index codes were added. Library quality was assessed using Qubit[R] 2.0 Fluorometer (Thermo Scientific) and Agilent Bioanalyzer 2100. Finally, the libraries were sequenced using an Illumina HiSeq 2500 platform, and 250-bp paired-end reads were generated. The data are publicly available in the NCBI Sequence Read Archive (accession number PRJNA507332).
For bioinformatic analyses, paired-end reads were assigned to samples based on their unique barcodes and truncated by cutting off the barcode and primer sequence. The paired-end reads were merged using FLASH (v1.2.7) (Magoc & Salzberg 2011) to obtain the raw tags. Quality filtering of the raw tags was performed under specific filtering conditions to obtain high-quality clean tags (Bokulich et al. 2013) using QIIME (v1.7.0) (Caporaso et al. 2010). The tags were compared with the reference database (Gold database) using the UCHIME algorithm (UCHIME algorithm) (Edgar et al. 2011) to detect chimaeras, which were removed (Haas et al. 2011) to finally obtain the effective tags. Sequence analyses were performed using UPARSE (v7.0.1001) and multiple sequence alignment was conducted using MUSCLE (version 3.8.31) (Edgar 2004). For each representative sequence, the Greengene database (DeSantis et al. 2006) was used with the RDP Classifier (v2.2) algorithm (Wang et al. 2007) to annotate taxonomic information. Subsequent analyses of alpha diversities were performed using the normalised output data. Alpha diversity was applied to analyse the complexity of species diversity through six indices, including observed-species, Chao1, Shannon, Simpson, ACE and Good's coverage. All indices were calculated using QIIME (v1.7.0) and analysed using R (v2.15.3).
Results and Discussion
Similar to other herbivores, the scimitar-horned oryx relies upon symbiosis with microorganisms in the digestive system to utilise cellulose and hemicellulose. In the present study, for the first time, the gut microbiota of the scimitar-horned oryx was characterised using faecal samples from five captive individuals. After removing low-quality reads and chimaeras, 367172 high-quality reads were obtained (mean 73434 reads per sample). The rarefaction curves tended to approach the saturation plateau (Fig. 1). The average Good's coverage was 99.8 %. For alpha diversity analyses, the sequence number was normalised to 61567 by random sub-sampling to standardise the sampling effort. Species diversity was measured using the Shannon and Simpson indices, and species richness was measured using Chao1 and ACE indices (Table 2).
A total of 15 prokaryotic phyla were identified in the five samples (Fig. 2), including Firmicutes (53.40-72.01 %), Bacteroidetes (12.94-23.72 %) and Proteobacteria (1.03-31.74 %), which accounted for > 96 % of all the sequences. However, microbial composition varied across samples. For instance, samples BC5 and BC3, which were both females, showed higher relative abundances of Proteobacteria (31.74 % and 14.57 %, respectively) than the remaining three samples (< 4.64 %), which were all males. In contrast, the abundances of Bacteroidetes in samples BC5 and BC3 (12.94 % and 18.44 %, respectively) were lower than those in the remaining three samples (20.67-23.72 %). There are many studies regarding faecal microbes in herbivores (Costa et al. 2012, Bian et al. 2013, Guan et al. 2017). According to Guan et al. (2017), Firmicutes (77.62 %), Bacteroidetes (18.28 %) and Tenericutes (1.34 %) were the most predominant phyla in wild sika deer, while Firmicutes (50.71 %) was the dominant phylum, followed by Bacteroidetes (31.99 %) and Proteobacteria (4.80 %), in captive sika deer. In a previous study, the bacterial communities in cape oryx were dominated by Firmicutes (42.81-55.29 %), Bacteroidetes (21.26-27.82 %) and Proteobacteria (3.05-7.14 %), represented by the families Ruminococcaceae, Lachnospiraceae, Prevotellaceae, Porphyromonadaceae, Succinivibrionaceae and Rikenellaceae (Chen et al. 2018).
At the family level, there were < 3.5 % unclassified sequences. The ten most common families are shown in Fig. 3. In samples BC5 and BC3, the predominant family was Ruminococcaceae, followed by Succinivibrionaceae, Lachnospiraceae, Rikenellaceae, Christensenellaceae and Bacteroidales_RF16_ group. Similarly, in the remaining three samples, Ruminococcaceae was the predominant family, followed by Rikenellaceae, Lachnospiraceae, Bacteroidales_RF16_group and Christensenellaceae. These comprised 84.69 %, 81.31 %, 85.47 %, 76.74 % and 89.12 % of the sequences in samples BC1, BC2, BC3, BC4 and BC5, respectively.
There were 159 genera detected in the five samples. Unclassified bacteria accounted respectively for 23.94 %, 21.12 %, 20.55 %, 19.69 % and 16.87 % of the sequences in samples BC1, BC2, BC3, BC4 and BC5, and the ten most common genera are shown in Fig. 4. Although the most common genera were identical in all the samples, the relative abundances of Succinivibrio in samples BC5 and BC3 were > 30 times higher than those in the remaining three samples. Previous studies have demonstrated that the most common genera were the bacteria responsible for cellulose digestion; thus, their higher abundances are associated with the herbivorous diet of the host (Daly et al. 2012, Jami & Mizrahi 2012).
Previous studies have shown that Proteobacteria are mainly associated with energy storage in humans, monkeys and mice (Bryant & Small 1956, Koren et al. 2012, Amato et al. 2014, Chevalier et al. 2015). Similarly, Sun et al. (2016) found that the relative abundances of Proteobacteria in the gut of Tibetan macaques were significantly higher in winter than in spring; they suspected that the increase in relative abundance during winter helps these macaques cope with the cold weather (Sun et al. 2016). Additionally, Wu et al. (2017) discovered that Succinivibrio spp. require carbon dioxide for anaerobic growth and ferment organic matter produced by the Krebs cycle to generate acetic and succinic acids. In this study, the high abundance of these bacteria might also be related to the energy produced by the species.
This work was financially supported by the National Natural Science Foundation of China (#41476091, #U1806213 and #ZR2018PC010) and the Doctoral Scientific Fund of Binzhou University (20118Y17).
Amato K.R., Leigh S.R., Kent A. et al. 2014: The role of gut microbes in satisfying the nutritional demands of adult and juvenile wild, black howler monkeys (Alouatta pigra). Am. J. Phys. Anthropol. 155: 652-664.
Backhed F., Ley R.E., Sonnenburg J.L. et al. 2005: Host-bacterial mutualism in the human intestine. Science 307: 1915-1920.
Barfod K.K., Roggenbuck M., Hansen L.H. et al. 2013: The murine lung microbiome in relation to the intestinal and vaginal bacterial communities. BMC Microbiol. 13: 303.
Bian G., Ma L., Su Y. & Zhu W. 2013: The microbial community in the feces of the white rhinoceros (Ceratotherium simum) as determined by barcoded pyrosequencing analysis. PLOS ONE 8: e70103.
Bokulich N.A., Subramanian S., Faith J.J. et al. 2013: Quality-filtering vastly improves diversity estimates from Illumina amplicon sequencing. Nat. Methods 10: 57-59.
Bryant M.P. & Small N. 1956: Characteristics of two new genera of anaerobic curved rods isolated from the rumen of cattle. J. Bacteriol. 72: 22-26.
Caporaso J.G., Kuczynski J., Stombaugh J. et al. 2010: QIIME allows analysis of high-throughput community sequencing data. Nat. Methods 7: 335-336.
Chen J., Shang S., Wu X.Y. et al. 2018: The microbial community in the feces of cape oryx (Oryx gazella) as determined by high-throughput Illumina sequencing technology. Pak. J. Zool. 50: 1601-1610.
Chevalier C., Stojanovic O., Colin D.J. et al. 2015: Gut microbiota orchestrates energy homeostasis during cold. Cell 163: 1360-1374.
Costa M.C., Arroyo L.G., Allen-Vercoe E. et al. 2012: 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 7: e41484.
Daly K., Proudman C.J., Duncan S.H. et al. 2012: Alterations in microbiota and fermentation products in equine large intestine in response to dietary variation and intestinal disease. Br. J. Nutr. 107: 989-995.
DeSantis T.Z., Hugenholtz P., Larsen N. et al. 2006: Greengenes, a chimera-checked 16S rRNA gene database and workbench compatible with ARB. Appl. Environ. Microbiol. 72: 5069-5072.
Edgar R.C. 2004: MUSCLE: multiple sequence alignment with high accuracy and high throughput. Nucleic Acids Res. 32: 1792-1797.
Edgar R.C., Haas B.J., Clemente J.C. et al. 2011: UCHIME improves sensitivity and speed of chimera detection. Bioinformatics 27: 2194-2200.
Guan Y., Yang H., Han S. et al. 2017: Comparison of the gut microbiota composition between wild and captive sika deer (Cervus nippon hortulorum) from feces by high-throughput sequencing. AMB Express 7: 212.
Haas B.J., Gevers D., Earl A.M. et al. 2011: Chimeric 16S rRNA sequence formation and detection in Sanger and 454-pyrosequenced PCR amplicons. Genome Res. 21: 494-504.
Hooper L.V., Midtvedt T. & Gordon J.I. 2002: How host-microbial interactions shape the nutrient environment of the mammalian intestine. Annu. Rev. Nutr. 22: 283-307.
Jami E. & Mizrahi I. 2012: Composition and similarity of bovine rumen microbiota across individual animals. PLOS ONE 7: e33306.
Koren O., Goodrich J.K., Cullender T.C. et al. 2012: Host remodeling of the gut microbiome and metabolic changes during pregnancy. Cell 150: 470-480.
Magoc T. & Salzberg S.L. 2011: FLASH: fast length adjustment of short reads to improve genome assemblies. Bioinformatics 27: 2957-2963.
Shanks O.C., Kelty C.A., Archibeque S. et al. 2011: Community structures of fecal bacteria in cattle from different animal feeding operations. Appl. Environ. Microbiol. 77: 2992-3001.
Singh K.M., Ahir V.B., Tripathi A.K. et al. 2012: Metagenomic analysis of Surti buffalo (Bubalus bubalis) rumen: a preliminary study. Mol. Biol. Rep. 39: 4841-4848.
Sun B., Wang X., Bernstein S. et al. 2016: Marked variation between winter and spring gut microbiota in free-ranging Tibetan macaques (Macaca thibetana). Sci. Rep. 6: 26035.
Thoetkiattikul H., Mhuantong W., Laothanachareon T. et al. 2013: Comparative analysis of microbial profiles in cow rumen fed with different dietary fiber by tagged 16S rRNA gene pyrosequencing. Curr. Microbiol. 67: 130-137.
Villalba J.J., Provenza F.D. & Bryant J. 2002: Consequences of the interaction between nutrients and plant secondary metabolites on herbivore selectivity: benefits or detriments for plants? Oikos 97: 282-292.
Villalba J.J., Provenza F.D. & Han G. 2004: Experience influences diet mixing by herbivores: implications for plant biochemical diversity. Oikos 107: 100-109.
Wang Q., Garrity G.M., Tiedje J.M. et al. 2007: Naive Bayesian classifier for rapid assignment of rRNA sequences into the new bacterial taxonomy. Appl. Environ. Microbiol. 73: 5261-5267.
Wu X., Zhang H., Chen J. et al. 2017: Analysis and comparison of the wolf microbiome under different environmental factors using three different data of Next Generation Sequencing. Sci. Rep. 7: 11332.
Xue D., Chen H., Chen F. et al. 2016: Analysis of the rumen bacteria and methanogenic archaea of yak (Bos grunniens) steers grazing on the Qinghai-Tibetan Plateau. Livest. Sci. 188: 61-71.
Shuai SHANG (1,2,3*), Longxiang LIU (2*), Jiayun REN (1,2), Hanjie ZHANG (2), Wangfeng ZHAO (2), Jikun XU (2), Qiang MAN (2), Liping ZHAO (2), Bin DONG (2), Jiangbao XIA (1,2) and Xuexi TANG (3)
(1) Shandong Key Laboratory of Eco-Environmental Science for the Yellow River Delta, Huanghe Road, Binzhou City, 256600 Shandong, China; e-mail: firstname.lastname@example.org, email@example.com
(2) College of Biological and Environmental Engineering, Binzhou University, Huanghe Road, Binzhou City, 256600 Shandong, China
(3) College of Marine Life Sciences, Ocean University of China, Qingdao, 266000 Shandong, China
Received 2 April 2019; Accepted 7 September 2019
(*) Corresponding Author
Table1. Basic information on scimitar-horned oryx used in this study. Scimitar-horned Health oryx Sex Age (years) condition B1 male 5 health B2 male 5 health B3 female 6 health B4 male 5 health B5 female 6 health Table 2. Bacterial community diversity of five samples measured using different indexes. Sample name Observed species Shannon Simpson Chao1 ACE BC1 1143 7.809 0.989 1521.477 1264.965 BC2 1099 7.758 0.988 1143.706 1145.836 BC3 1113 7.455 0.978 1149.093 1159.498 BC4 1123 8.123 0.992 1169.146 1171.415 BC5 1035 6.317 0.924 1070 1076.664 Sample name Goods coverage BC1 0.997 BC2 0.998 BC3 0.998 BC4 0.998 BC5 0.999
|Printer friendly Cite/link Email Feedback|
|Author:||Shang, Shuai; Liu, Longxiang; Ren, Jiayun; Zhang, Hanjie; Zhao, Wangfeng; Xu, Jikun; Man, Qiang; Zha|
|Date:||Nov 1, 2019|
|Previous Article:||Influence of seasonality, environmental and anthropic factors on crop damage by wild boar Sus scrofa.|
|Next Article:||Using video surveillance to monitor feeding behaviour and kleptoparasitism at Eurasian lynx kill sites.|