First Report of Fecal Microflora of Wild Bar-Headed Goose in Tibet Plateau

Abstract

The bar-headed goose (Anser indicus) has two black spots on its head. It is considered an important bird in China. It breeds in plateau lakes, especially saltwater lakes, and swamp areas. However, the intestinal flora of wild bar-headed geese in the Tibet Autonomous Region is currently not known. In this study, 16S rDNA sequencing was performed on the intestinal microbes of wild bar-headed geese. A total of 513,505 reads of raw data were obtained, and the results analyzed the average number of 128,376 ± 2,392 reads per sample. The microbiota of all samples consists of 10 main bacterial phyla, including Firmicutes, Proteobacteria, Bacteroidetes, Actinobacteria, Cyanobacteria, Patescibacteria, Deferribacteres, Planctomy-cetes, Fusobacteria, and Tenericutes. The results indicated that Firmicutes (67.34%) was the predominant phylum, followed by Proteobacteria (29.03%) and Cyanobacteria (1.97%). In our research, we identified the intestinal flora of the wild bar-headed goose, which provides valuable information for further research on the gene function of the bar-headed goose and the intestinal flora of wild animals. These findings are also useful and valuable for genetic and high-altitude research in the Tibet Autonomous Region.

Introduction

Research on the development and progress of biotechnology (1) has shown that gut microbes not only affect the health of the host but also play an important role in the occurrence of diseases including obesity (2), cancer (3), diabetes (4), and increased cardiovascular risks through different metabolic processes (5). Among vertebrates, birds are considered to be the most common and diverse on earth, with more than 10,000 living species (6). It is reported that a variety of birds in the world constantly migrate from one place to another every year. They can migrate over vast geographical areas and cross biological and geographical boundaries (7). Birds can act as a source of infections for public health and animals through direct contact or as carriers of pathogens including avian and zoonotic pathogens and antibiotic-resistant bacteria (8). For these reasons, birds play an important role in the spread of microorganisms from one place to another, thereby affecting their dynamics ecology and the evolution of various viruses and bacteria. In the published literature, research on the gut microbiota of birds is very limited and most of them are researches on some artificially farmed economic species.

The bar-headed geese are a species unique to Asia, and most bar-headed geese breed in different places during the winter season including Tibet, south central Tibet, and India (9, 10). During migration, the bar-headed geese avoid areas with high altitudes, harsh climates, and barren vegetation and find suitable habitats that are convenient for harvesting food. Birds may carry and spread different infectious agents during long-distance migration including avian influenza virus and other pathogens that pose a threat to wildlife and public health.

In recent years, advancements in sequencing technology have rapidly improved (researchers' ability in) understanding of the intestinal microbiota of different animals including humans, mice, and other mammals. However, the application of sequencing technology in the study of bird intestinal flora is still limited (11). Compared with mammals, the gastrointestinal tract of birds is relatively short, and the time to digest food is also short, which enables them to form highly selective and adaptable microbiota. Studies have shown that the diversity of an intestinal bacterial community of white-headed cranes in different regions has been investigated using 16S rDNA sequencing techniques. Previous studies have found that there are significant differences in the structure and diversity of the intestinal bacterial community of hooded cranes in different locations in winter and different climatic conditions that affect the composition of the intestinal bacterial community (12). Few studies have reported sequencing techniques for feeding habits (insectivorous or omnivorous) to determine the differences in the digestive tract microbiota of passerine birds like those in New Guinea (13). At present, there are few reports on the use of 16S rDNA sequencing technology to study the intestinal microbes of the wild bar-headed geese in Tibet. Therefore, this study used 16S rDNA sequencing technology to analyze the intestinal microbial diversity of wild bar-headed geese in Tibet.

Materials and Methods

Sample Collection

Fecal samples were collected from four different bar-headed geese living in the southern part of Naidong District, Shannan City (29°10′46″N, 91°46′15″E, elevation 3,552 m), China. Fresh fecal samples were collected and stored in sterile test tubes. All the freshly excreted fecal samples by bar-headed geese were collected, and the middle part of the sample was processed (Figure 1). All samples were transported to the laboratory using a −20°C portable refrigerator and stored at −80°C until further processing.

Figure 1

DNA Extraction, PCR Amplification, and High-Throughput Sequencing

Microbial DNA was extracted using the HiPure Stool DNA Kits (Magen, Guangzhou, China) according to the manufacturer's guidelines. The extracted DNA was quantified and evaluated for purity using NanoDrop 2000 UV-vis spectrophotometer and 1% agarose gel electrophoresis. The hypervariable region V3–V4 of the 16S rDNA genes was amplified using primers 341F(5′-CCTACGGGNGGCWGCAG-3′) and 806R(5′-GGACTACHVGGGTATCTAAT-3′) (14). The PCR conditions were as follows: initial denaturation at 94°C for 2 min, denaturation at 98°C for 10 s, 62°C annealing for 30 s, 68°C extension for 30 s using 30 cycles, and a final extension at 68°C for 5 min. PCR reactions were performed in triplicate 50 μl mixture containing 5 μl of 10 × KOD buffer, 5 μl of 2 mM dNTPs, 3 μl of 25 mM MgSO₄, 1.5 μl of each primer (10 μM), 1 μl of KOD polymerase, and 100 ng of template DNA. Amplicons were extracted on 2% agarose gel and purified using the AxyPrep DNA Gel Extraction Kit (Axygen Biosciences, Union City, CA, USA) according to the manufacturer's instructions. All the amplicons were quantified using ABI Step One Plus Real-Time PCR System (Life Technologies, Foster City, USA). The purified amplicons were pooled in equimolar and paired-end sequenced (2 × 250) on an Illumina platform according to the standard protocols.

Bioinformatic and Statistical Analyses

Raw data containing adapters or low quality reads may affect the following assembly and analysis. In order to get high quality clean reads, the raw reads were further filtered following the guidelines using FASTP to remove the reads containing more than 10% of unknown nucleotides and to remove the reads containing less than 80% of bases with quality (Q-value) > 20. After that the paired end clean reads were merged as raw tags using FLSAH (15) with a minimum overlap of 10 bp and mismatch error rates of 2%. Noisy sequences of raw tags were filtered by QIIME (16) pipeline under specific filtering conditions (17) to obtain the high-quality clean tags. Clean tags were searched against the reference database (http://drive5.com/uchime/uchime_download.html) to perform reference-based chimera checking using UCHIME algorithm (http://www.drive5.com/usearch/manual/uchime_algo.html). After that, all the chimeric tags were removed, and effective tags were finally obtained and used for further analysis. The effective tags were clustered into operational taxonomic units (OTUs) of ≥ 97% using UPARSE (18) pipeline. The tag sequence with highest abundance was selected as representative sequence within each cluster. The representative sequences were classified into organisms by a naive Bayesian model using RDP classifier (19) based on SILVA (20) Database (https://www.arb-silva.de/), with the confidence threshold values ranging from 0.8 to 1.

The abundance statistics of each taxonomy were visualized using Krona (21). The stacked bar plot of the community composition was visualized in R project ggplot2 package (version 2.2.1). Chao1, Simpson, and all other alpha diversity indexes were calculated in QIIME.FAPROTAX database (Functional Annotation of Prokaryotic Taxa), and associated software (22) (version 1.0) were used for generating the ecological functional profiles of bacteria. The Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analysis of the OTUs was inferred using Tax4Fun (version 1.0) or PICRUSt (version 2.1.4).

Statistical Analysis

The abundance statistics of each taxonomy were visualized using Krona (21). Chao1, Simpson and all other alpha diversity indexes were calculated in QIIME. Alpha index comparison between groups was calculated by Welch's t-test and Wilcoxon rank test in R project. FAPROTAX database and associated software (22) (version 1.0) were used for generating the ecological functional profiles of bacteria. Analysis of function difference between groups was calculated by Welch's t-test, Wilcoxon rank test, Kruskal–Wallis H-test, and Tukey's honestly significant difference (HSD) test in R Project Vegan package (version 2.5.3).

Results

Effective Sequence Quality Assessment

Results on the microbiome analysis indicated that totals of 118,238, 120,100, 123,894, and 121,612 original sequences were acquired from B1, B2, B3, and B4 groups, respectively (Table 1). After eliminating the unqualified data, a total number of 474,918 high-quality tags were achieved from all the samples. The length distribution of all samples ranged from 200 to 475 bp.

All the optimized sequences are aligned to OTU representative sequences by UPARSE software, and sequences with more than 97% similarity with representative sequences are selected to generate OTUs. After classification matching, a total of 1,106 OTUs (B1 = 229, B2 = 262, B3 = 282, B4 = 333) were obtained (Figure 2A).

Figure 2

The Diversity of the Intestinal Microbial

Generally, the alpha diversity of the gut microbial community can be analyzed by mainstream alpha diversity indexes including species (Sob), Chao1, ACE, Shannon, Simpson, Good's Coverage, pielou, and PD-whole tree. Our results showed that both the Shannon curve and the PD-whole tree curve have reached the plateau, indicating that the sequencing results can reflect the diversity of the present samples. Furthermore, the rank abundance curve is wide and the downward trend is flat, showing excellent abundance and evenness (Figures 2B,C and Table 2). The repetition between samples is well, and it can effectively reflect the bar-headed geese that live in Tibet.

Microbial Community Structure Analysis

This study analyzed the composition of intestinal microbes at the phylum level. The results on 16S rDNA gene analysis indicated the presence of high proportions of Firmicutes, Proteobacteria, and Cyanobacteria (Figure 3). As for the fecal samples of four groups (B1, B2, B3, and B4) of bar-headed geese, the main bacterial phylum was Firmicutes (69.60, 54.92, 91.76, and 53.09%), Proteobacteria (26.89, 42.64, 5.12, and 41.47%), and Cyanobacteria (2.94, 1.46, 2.43, and 1.03%), accounting for approximately 97% of the taxonomic groups identified (Figure 3A). The results indicated that the predominant phylum of the B3 group was Firmicutes (91.76%), and Proteobacteria (5.12%) was the secondary phylum (5.12%), which was significantly different from other groups. At the genus level, Bacillus (41.5%) was the most dominant genus, followed by Solibacillus (16.35%), Exiguobacterium (3.6%), Acinetobacter (3.34%), Lysinibacillus (3%), and Pseudomonas (2.1%) (Figures 3B–D).

Figure 3

Function Prediction

The horizontal coordinates of the stack diagram in samples are different, and the columns with different colors represent the relative abundance of different ecological functions. Among them, the results showed that the abundance of the top 11 functions and the abundance of other functions were combined into other categories (Figure 4).

Figure 4

Pathway results showed that bacteria were mainly related to metabolism, environmental information processing, genetic information processing, cellular processes, multiple levels of human diseases, and organismal systems (Table 3).

Discussion

Bar-headed geese not only survive on the plains but can also adjust to the hostile environment of the plateau. As an adaptive species of plateaus and plains, the bar-headed goose is one of the birds that efficiently pay attention (9), and some researchers have conducted special studies on the changes in oxygen in their bodies. Some researchers have also studied the temperature regulation (17) and childbirth (23) of bar-headed geese during their high-altitude migration. In this study, the 16S rDNA amplicons were sequenced for the first time in the feces of bar-headed geese in Tibet, and the diversity of intestinal microbes has also been investigated. Comparing this with previous studies on the microbial composition of birds at the phylum level, Firmicutes, Proteobacteria, Bacteroidetes, Actinobacteria, and Fusobacteria were dominant (24). The findings of this study regarding 16s rRNA sequencing are similar to an earlier study on the intestinal microbial diversity of cultured bar-headed goose (25). Previously, different studies also sequenced the diversity of the gut microbiota of bar-headed geese (wild, semi-captive, and captive), and the results indicated that the dominant phyla included Firmicutes, Proteobacteria, Actinobacteria, Bacteroidetes, and Fusobacteria (26) and in this study 23.33% Turicibacter in Firmicutes in the semi-artificial feeding group. A total of 77.67% Lactococcus in Firmicutes in the wild feeding group and 51% SMB53 in the feeding group were recorded. However, at the genus level in this study, Bacillus (41.5%) is the most dominant genus. Among Firmicutes (67.34%), Proteobacteria (29.03%) and Cyanobacteria (1.97%), the most dominant phylum is Firmicutes (67.34%). This may be due to the different living environment of the bar-headed geese. The samples in this study were collected from bar-headed geese living in Tibet more than 3,000 m above sea level where the growth of plants varies at different altitudes. Therefore, the variations in the results of our study from different other reports may be associated to this reason.

Over the past few years, frequent and extensive investigations have been conducted on the fecal genome of birds (27). The diversity of the fecal genome of birds is an important part of the gut microbiota. In this study, the wild bar-headed goose living in Tibet is an endangered protected animal. Its intestinal flora is more scientific and logical than that of the captive-bred bar-headed goose. Moreover, due to the living environment in Tibet, its intestinal flora also has unique characteristics of Tibet. Firmicutes and Bacteroidetes were the most common bacteria in all animal feces in previous studies (28), and Firmicutes are the most important microflora to promote the decomposition of cellulose by host gastrointestinal microorganisms (29). The results showed that Firmicutes and Bacteroidetes were the main bacteria encoding active enzymes of carbohydrate (30). The ratio of Firmicutes to Bacteroidetes in the gut (F/B) affects the host's ability to obtain energy from food (31). The relative abundance of Firmicutes in poultry cecum is positively correlated with body weight gain (32) and egg production performance (33). Another study showed that early use of antibiotics in chickens reduced the relative abundance of Firmicutes and increased the relative abundance of Proteobacteria. In the later stage, the immune activity of individual T cells decreases with low relative abundance of Firmicutes (32), indicating the influence of Firmicutes and Proteobacteria on the immune functions of the host. The third dominant phylum in feces is Cyanobacteria. Previously, lactic acid bacteria were added to chicken feed to improve the growth performance and immune response of chickens. The same bacteria, Cyanobacteria, Proteobacteria, Bacteroides, and Actinomycetes, were also detected (34). In terms of genera, the top two genera sequenced are Bacillus (41.5%) and Solibacillus (16.35%). Although Bacillus accounted for 41%, the study indicated no significant presence of Bacillus. Bacillus belongs to bacillus family, and the genus is Bacillus, which is a kind of gram-positive bacteria that can produce endophytic spores and a large amount of calcium pyridine dicarboxylic acid. The antibacterial substances produced by bacillus generally have a wide antibacterial spectrum, which can kill bacteria including drug-resistant strains, some fungi, parasites, certain viruses, and tumor cells and can bind lipopolysaccharides and neutralize endotoxins. Therefore, the researchers focused their studies on the diversity of intestinal microbes. The probiotics prepared by bacillus play an important role in the treatment of disorders in intestinal microflora, candida infection, and wound infections.

In this study, compared with the 16S sequencing results of previous reports, the focus is on the inclusion of chloroplast-containing fecal Yanobacteria at the phylum level. This may be due to the lack of food supply in Tibet in winter and the need to consume chloroplast-containing plants for survival. In terms of genera, Bacillus is the most dominant, which reduces the incidence of disease in poultry. These are the self-regulation mechanisms of bar-headed geese to adapt to life in Tibet. For wild animals raised in captivity, although wild animal resources are protected, it has caused the destruction of microbial diversity. While protecting animal germplasm resources, we should also make appropriate ecological improvements to wild animals to protect the diversity of microbial communities.

Conclusion

In our study, we investigated the intestinal microbiome in the feces of wild bar-headed geese which provides valuable resources for further research on the gene functions of different bar-headed geese and the intestinal microbiome of wild animals. These are also valuable for genetic and high-altitude research in the Tibet Autonomous Region.

Funding

This study was supported by the major science and technology projects of the Tibet autonomous region (XZ202101ZD0005N), Key R&D plan of Bayi District, Nyingchi City (2021-GX-SY-01), and the Basic Research Funds of China Agricultural University (2021TC002).

Publisher's Note

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.

References

• 1.

  SchroederBOBäCkhedF. Signals from the gut microbiota to distant organs in physiology and disease. Nat Med. (2016) 22:1079. 10.1038/nm.4185

  • Pubmed Abstract
  • CrossRef
  • Google Scholar
  • View reference in article

• 2.

  SanmiguelCGuptaAMayerEA. Gut microbiome and obesity: a plausible explanation for obesity. Curr Obes Rep. (2015) 4:250–61. 10.1007/s13679-015-0152-0

  • Pubmed Abstract
  • CrossRef
  • Google Scholar
  • View reference in article

• 3.

  Pevsner-FischerMTuganbaevTMeijerMZhangSHZengZRChenMHet al. Role of the microbiome in non-gastrointestinal cancers. World J Clin Oncol. (2016) 7:14. 10.5306/wjco.v7.i2.200

  • Pubmed Abstract
  • CrossRef
  • Google Scholar
  • View reference in article

• 4.

  ForslundKHildebrandFNielsenTFalonyGChatelierELSunagawaSet al. Disentangling type 2 diabetes and metformin treatment signatures in the human gut microbiota. Nature. (2017) 545:116. 10.1038/nature15766

  • Pubmed Abstract
  • CrossRef
  • Google Scholar
  • View reference in article

• 5.

  VijayAAstburySPanayiotisLMarquesFZSpectorTDMenniCet al. Dietary interventions reduce traditional and novel cardiovascular risk markers by altering the gut microbiome and their metabolites. Front Cardiovasc Med. (2021) 8:691564. 10.3389/fcvm.2021.691564

  • Pubmed Abstract
  • CrossRef
  • Google Scholar

• 6.

  PrumROBervJSDornburgAFieldDJTownsendJPLemmonEMet al. A comprehensive phylogeny of birds (aves) using targeted next-generation dna sequencing. Nature. (2015) 526:569–73. 10.1038/nature19417

  • Pubmed Abstract
  • CrossRef
  • Google Scholar
  • View reference in article

• 7.

  OlivaresMCalbetASaizE. Effects of multigenerational rearing, ontogeny and predation threat on copepod feeding rhythms. Aquat Ecol. (2020) 2020:2194–201. 10.1007/s10452-020-09768-8

  • CrossRef
  • Google Scholar
  • View reference in article

• 8.

  ElmbergJBergCLernerHWaldenstrmJHesselR. Potential disease transmission from wild geese and swans to livestock, poultry and humans: a review of the scientific literature from a One Health perspective. Infect Ecol Epidemiol. (2017) 7:1300450. 10.1080/20008686.2017.1300450

  • Pubmed Abstract
  • CrossRef
  • Google Scholar
  • View reference in article

• 9.

  JessicaUMJuliaMYBevACWilhelminaJLucyAHWilliamKM. Reduced metabolism supports hypoxic flight in the high-flying bar-headed goose (Anser indicus). Elife. (2019) 8:e44986. 10.7554/eLife.44986

  • Pubmed Abstract
  • CrossRef
  • Google Scholar
  • View reference in article

• 10.

  ZhangYHaoMTakekawaJYLeiFYanBProsserDJet al. Tracking the autumn migration of the bar-headed goose (anser indicus) with satellite telemetry and relationship to environmental conditions. Int J Zool. (2011) 2011-10-16:620–2. 10.1007/s10336-010-0492-1

  • CrossRef
  • Google Scholar
  • View reference in article

• 11.

  BennettDCTunHMKimJELeungFCChengKM. Characterization of cecal microbiota of the emu (dromaius novaehollandiae). Vet Microbiol. (2013) 166:304–10. 10.1016/j.vetmic.2013.05.018

  • Pubmed Abstract
  • CrossRef
  • Google Scholar
  • View reference in article

• 12.

  DongYXiangXZhaoGSongYZhouL. Variations in gut bacterial communities of hooded crane (grus monacha) over spatial-temporal scales. PeerJ. (2019) 7:e7045. 10.7717/peerj.7045

  • Pubmed Abstract
  • CrossRef
  • Google Scholar
  • View reference in article

• 13.

  BodawattaKHKaterinaSJøNssonKAMichaelP. Comparative analyses of the digestive tract microbiota of new guinean passerine birds. Front Microbiol. (2018) 9:1830. 10.3389/fmicb.2018.01830

  • Pubmed Abstract
  • CrossRef
  • Google Scholar
  • View reference in article

• 14.

  GuoMFWu HaoGQiQRongLLiNet al. Bacillus subtilis improves immunity and disease resistance in rabbits. Front Immunol. (2017) 8:354. 10.1038/nature22318

  • Pubmed Abstract
  • CrossRef
  • Google Scholar
  • View reference in article

• 15.

  MagočTSalzbergSL. Flash: fast length adjustment of short reads to improve genome assemblies. Bioinformatics. (2011) 27:2957–63. 10.1093/bioinformatics/btr507

  • Pubmed Abstract
  • CrossRef
  • Google Scholar
  • View reference in article

• 16.

  CaporasoJKuczynskiJStombaughJBittingerKBushmanF. QIIME allows integration and analysis of high-throughput community sequencing data. Nat Methods. (2010) 7:335–6. 10.1038/nmeth.f.303

  • Pubmed Abstract
  • CrossRef
  • Google Scholar
  • View reference in article

• 17.

  ParrNBishopCMBatbayarNPatrickJBBeverlyCWilliamKMet al. Tackling the Tibetan Plateau in a down suit: insights into thermoregulation by bar-headed geese during migration. J Exp Biol. (2020) 223:jeb226639. 10.1242/jeb.203695

  • Pubmed Abstract
  • CrossRef
  • Google Scholar
  • View reference in article

• 18.

  EdgarRC. UPARSE: highly accurate OTU sequences from microbial amplicon reads. Nat Methods. (2013) 10:996. 10.1038/nmeth.2604

  • Pubmed Abstract
  • CrossRef
  • Google Scholar
  • View reference in article

• 19.

  WangQGarrityGMTiedjeJMColeJR. Naive Bayesian classifier for rapid assignment of rRNA sequences into the new bacterial taxonomy. Appl Environ Microbiol. (2007) 73:5261–7. 10.1128/aem.00062-07

  • Pubmed Abstract
  • CrossRef
  • Google Scholar
  • View reference in article

• 20.

  ElmarPChristianQKatrinKFuchsBMWolfgangLJrgPet al. Silva: a comprehensive online resource for quality checked and aligned ribosomal rna sequence data compatible with arb. Nucleic Acids Res. (2007) 21:7188–96. 10.1093/nar/gkm864

  • Pubmed Abstract
  • CrossRef
  • Google Scholar
  • View reference in article

• 21.

  OndovBDBergmanNHPhillippyAM. Interactive metagenomic visualization in a web browser. BMC Bioinformatics. (2011) 12:385. 10.1186/1471-2105-12-385

  • Pubmed Abstract
  • CrossRef
  • Google Scholar
  • View reference in article

• 22.

  LoucaSParfreyLWDoebeliM. Decoupling function and taxonomy in the global ocean microbiome. Science. (2016) 353:1272. 10.1126/science.aaf4507

  • Pubmed Abstract
  • CrossRef
  • Google Scholar
  • View reference in article

• 23.

  LagueSLChuaBFarrellAPWangYMilsomWK. Altitude matters: differences in cardiovascular and respiratory responses to hypoxia in bar-headed geese reared at high and low altitudes. J Exp Biol. (2016) 219:1974–84. 10.1242/jeb.132431

  • Pubmed Abstract
  • CrossRef
  • Google Scholar
  • View reference in article

• 24.

  DavidWWaiteMTaylorW. Exploring the avian gut microbiota: current trends and future directions. Front Microbiol. (2015) 6:673. 10.3389/fmicb.2015.00673

  • Pubmed Abstract
  • CrossRef
  • Google Scholar
  • View reference in article

• 25.

  WangWWangFWangASharshovKWangS. Characterization of the microbiome along the gastrointestinal tracts of semi-artificially reared bar-headed geese (anser indicus). Folia Microbiol. (2019) 65:533–43. 10.1007/s12223-019-00758-4

  • Pubmed Abstract
  • CrossRef
  • Google Scholar
  • View reference in article

• 26.

  WangWCaoJLiJRYangFLiZLiLX. Comparative analysis of the gastrointestinal microbial communities of bar-headed goose (Anser indicus) in different breeding patterns by high-throughput sequencing. Microbiol Res. (2016) 182:59–67. 10.1016/j.micres.2015.10.003

  • Pubmed Abstract
  • CrossRef
  • Google Scholar
  • View reference in article

• 27.

  ZhangGRahbekCGravesGLeiFJarvisEGilbertM. Genomics: bird sequencing project takes off. Nature. (2015) 522:34. 10.1038/522034d

  • Pubmed Abstract
  • CrossRef
  • Google Scholar
  • View reference in article

• 28.

  LeyREHamadyMLozuponeCTurnbaughPJRameyRRBircherJSet al. Evolution of mammals and their gut microbes. Science. (2008) 320:1647–51. 10.1126/science.1155725

  • Pubmed Abstract
  • CrossRef
  • Google Scholar
  • View reference in article

• 29.

  PatelVPatelAKParmarNRPatelABReddyBJoshiCG. Characterization of the rumen microbiome of indian kankrej cattle (bos indicus) adapted to different forage diet. Appl Microbiol Biotechnol. (2014) 98:9749–61. 10.1007/s00253-014-6153-1

  • Pubmed Abstract
  • CrossRef
  • Google Scholar
  • View reference in article

• 30.

  KaoutariAEArmougomFGordonJIRaoultDHenrissatB. The abundance and variety of carbohydrate-active enzymes in the human gut microbiota.Nat. Rev. Microbiol. (2013) 11:497–504. 10.1038/nrmicro3050

  • Pubmed Abstract
  • CrossRef
  • Google Scholar
  • View reference in article

• 31.

  TurnbaughPJLeyREMahowaldMAMagriniVMardisERGordonJI. An obesity-associated gut microbiome with increased capacity for energy harvest. Nature. (2006) 444:1027–31. 10.1038/nature05414

  • Pubmed Abstract
  • CrossRef
  • Google Scholar
  • View reference in article

• 32.

  ZhangLLiJYunTTQiWTLiAK. Effects of pre-encapsulated and pro-encapsulated enterococcus faecalis on growth performance, blood characteristics, and cecal microflora in broiler chickens. Poult Sci. (2015) 94:2821–30. 10.3382/ps/pev262

  • Pubmed Abstract
  • CrossRef
  • Google Scholar
  • View reference in article

• 33.

  YangZZhangCWangJCeliPShanZ. Characterization of the intestinal microbiota of broiler breeders with different egg laying rate. Frontiers in Veterinary Science. (2020) 7:599337. 10.3389/fvets.2020.599337

  • Pubmed Abstract
  • CrossRef
  • Google Scholar
  • View reference in article

• 34.

  KumarSLeeHCSunKJKimDByunSJ. Inclusion of lactobacillus salivarius strain revealed a positive effect on improving growth performance, fecal microbiota and immunological responses in chicken. Arch Microbiol. (2020) 203:847–53. 10.1007/s00203-020-02088-3

  • Pubmed Abstract
  • CrossRef
  • Google Scholar
  • View reference in article