Comparing microbiotas in the upper aerodigestive and lower respiratory tracts of lambs
© The Author(s). 2017
Received: 19 July 2016
Accepted: 19 October 2017
Published: 27 October 2017
Recently, the importance of the lung microbiota during health and disease has been examined in humans and in small animal models. Whilst sheep have been proposed as an appropriate large animal model for studying the pathophysiology of a number of important human respiratory diseases, it is clearly important to continually define the limits of agreement between these systems as new concepts emerge. In humans, it has recently been established that the lung microbiota is seeded by microbes from the oral cavity. We sought to determine whether the same was true in sheep.
We took lung fluid and upper aerodigestive tract (oropharyngeal) swab samples from 40 lambs (7 weeks old). DNA extraction was performed, and the V2-V3 region of the 16S rRNA gene was amplified by PCR then sequenced via Illumina Miseq. Oropharyngeal swabs were either dominated by bacteria commonly associated with the rumen or by bacteria commonly associated with the upper aerodigestive tract. Lung microbiota samples did not resemble either the upper aerodigestive tract samples or reagent-only controls. Some rumen-associated bacteria were found in lung fluids, indicating that inhalation of ruminal bacteria does occur. We also identified several bacteria which were significantly more abundant in lung fluids than in the upper aerodigestive tract swabs, the most predominant of which was classified as Staphylococcus equorum.
In contrast to humans, we found that the lung microbiota of lambs is dissimilar to that of the upper aerodigestive tract, and we suggest that this may be related to physiological and anatomical differences between sheep and humans. Understanding the comparative physiology and anatomy underlying differences in lung microbiota between species will provide a foundation upon which to interpret changes associated with disease and/or environment.
The use of 16S ribosomal RNA (rRNA) gene sequencing has facilitated the study of difficult to culture, low biomass microbial communities present in the lower respiratory tract. The impact of the lung microbiota on human health is a rapidly growing area of research. In order to understand this impact, it is important to also understand the lung microbiota dynamics during health and to include healthy controls in disease studies. To achieve this, the majority of previous studies have relied on human volunteers.
However, many individuals are hesitant to participate in research bronchoscopy due to the perceived inconvenience and a fear of complications , despite the low risk involved. Mice and rats have been used to explore the relationship between the lung microbiota and airway inflammation [2–4], microbiota at different body sites , the environment , acute lung injury  and antibiotic  and corticosteroid exposure . However, rodents are of limited use when exploring spatial or longitudinal lung microbiota dynamics due to their small lung size. Recognising the utility of large animal models in this regard, and the anatomical and immunological relevance of sheep as models [10–13], our group has previously used this species to explore the changes in the lung microbiota upon Pseudomonas aeruginosa infection  and to explore the spatial variability present within the healthy lung .
Subclinical microaspiration of pharyngeal secretions is a feature of health and this can contribute to the lung microbiota composition , and the microbiome of the human lungs more closely resembles that of the mouth than that of the nose or the lower gastrointestinal tract . It is not yet known whether the same relationship holds for species other than humans. In ruminating sheep, where the oropharynx is exposed to ruminal contents on a frequent basis, one would anticipate that lung microbiota would similarly reflect this influence. In this paper, we find that the presence of rumen-like bacteria in the upper aerodigestive tract is correlated with changes in the lung microbiota and rumen-type bacteria are present in lamb lungs. We also identify bacteria which are more indicative of the lungs than the oropharynx, indicating that the presence of the sheep lung microbiota is not merely due to passive diffusion of microbes from the upper aerodigestive tract.
Animals and sampling
Scottish Mule X Suffolk lambs (20 males and 20 females, unweaned), raised on pasture from 48 h after birth, were used in this study. These lambs were part of a study on the animal welfare implications of prenatal stress which was approved by Scotland’s Rural College’s (SRUC) Animal Experiments Committee and was conducted under Home Office licence. All lambs were raised by their dams, and prior to euthanasia, their only food sources were ewe’s milk and pasture. At 7 weeks old (mean age = 48.8 days ± 0.8 standard deviation (SD); mean weight ± SD = 20.6 kg ± 2.6 kg), the lambs were euthanised by barbiturate overdose then the cadavers were transported from the farm to the dissecting room (~ 5 min). Oropharyngeal swabs were taken using cotton-tipped swabs (Swab Plain Wood Cotton Tip Sterile (710-0181), Copan, Italy). To prevent oral contamination, the swabs were stored in protective plastic sheaths from which the swab could be advanced and retracted once it was positioned at the sampling site. Swabs were then transferred into a new plastic sheath and stored on ice.
The ventral aspect of the neck was shaved, and a sterile scalpel used to incise through the skin and subcutaneous tissues to expose the ventral surface of the trachea. A sampling site was identified on the exposed ventral surface, and the trachea cranial to this site was completely closed off by both string ligature and clamp placement. The selected sampling site was then heat seared, and 50 ml of sterile phosphate-buffered saline (PBS) was injected through the seared section into the tracheal lumen. The head and neck were oriented such that the PBS would flow caudally down the thorax. A second clamp was immediately placed caudal to the site of injection to prevent backflow, leakage and potential contamination. The lamb cadavers were then tipped so that the PBS would run caudally into their lungs and then tipped back again so that the fluid would collect in the tracheal lumen immediately caudal to the position of the second clamp. A sampling site identified on the ventral surface of the trachea was seared, and a needle and syringe were used to collect the pooled fluid. On average, 4 ± 1.7 ml (mean ± SD) of lung fluid was collected per animal. Lung fluid was stored on ice until further processing. Oropharyngeal swabs were sterilely cut into 500 μl PBS. Lung fluids were centrifuged at 13,000g for 5 min. The supernatant was removed, and the pellets were resuspended in 500 μl PBS. The oropharyngeal swabs and lung fluids were stored at −80 °C until DNA extraction.
DNA extraction, amplification and sequencing
DNA extractions using the PowerSoil DNA isolation Kit (Mo Bio, Carlsbad, USA) and quantitative PCR (qPCR) using the 16S rRNA gene qPCR primers UniF340 (5-ACTCCTACGGGAGGCAGCAGT-3) and UniR514 (5-ATTACCGCGGCTGCTGGC-3) were performed as described previously . Extraction kit reagent controls, consisting of reagent-only extractions, were produced for every day DNA extractions were performed. PBS controls were created by extracting DNA from 500 μl of the PBS which had also been used to wash out the lamb lungs. A mock community control was included which has been described previously .
A nested PCR reaction was used to produce amplicons for sequencing; this technique was chosen to reduce PCR bias caused by barcoded primers . The first round of PCR amplified the V1-V4 16S hypervariable regions using the primers 28F (5-GAGTTTGATCNTGGCTCAG-3) and 805R (5-GACTACCAGGGTATCTAATC-3). The conditions were 94 °C for 2 min followed by 20 cycles of 94 °C for 1 min, 55 °C for 45 s and 72 °C for 1.5 min followed by a final extension step of 72 °C for 20 min. Clean-up was performed using the AMPure XP PCR purification system (Beckman Coulter, Brea, USA).
In a previous study, we found that PCR bias in high template samples could be reduced by diluting amplicons from the first round of PCR to a similar concentration to those of lung fluid samples . Therefore, in this study, we used our qPCR values to calculate the dilutions needed to achieve this. The second round of PCR used the barcoded V2-V3 primers 104F (5-GGCGVACGGGTGAGTAA-3) and 519R (5-GTNTTACNGCGGCKGCTG-3). The dilutions and barcoded primers used for each sample can be found in Additional file 1. The PCR conditions were 98 °C for 30 s followed by 20 cycles of 98 °C for 10 s, 67 °C for 30 s and 72 °C for 10 s followed by a final extension step of 72 °C for 2 min. The amplicons were again purified using the AMPure XP PCR purification system.
Bioinformatic and statistical analysis
Samples were sequenced via either Illumina Miseq or Hiseq runs (Illumina, San Diego, USA) (Additional file 1) producing 250 base pair paired-end reads. Cutadapt  was used to remove primers. Quality control and taxonomic assignment of sequences was carried out within mothur  following a protocol created by the mothur developers , adjusted to suit our dataset . Sequences were subsampled before statistical analysis. The sequencing error rate, principal coordinate analysis (PCOA) graphs, analyses of molecular variance (AMOVAs), Good’s coverage analyses , richness (Chao 1 index) and diversity (inverse Simpson index) calculations and indicator analyses  were all calculated within mothur. Clustering of microbial communities into metacommunities was also carried out within mothur using a probabilistic modelling technique based upon work by Holmes et al. . The significance of differences between the diversity and richness of groups was calculated using either the two-sample t test (normal data) or the Mann-Whitney U test (non-normal data) within Minitab 16 for Windows (Minitab, Coventry, UK). Heatmaps were constructed in R Version 3.2.2  using the Vegan , RColorBrewer , gplots  and heatplus  packages. Clustering within heatmaps was performed using the Bray-Curtis dissimilarity . Sequences can be accessed via the Bioproject accession number PRJNA317719.
Quality assurance of methodology
Bacterial OTUs found to be > 5% abundant in extraction kit reagent controls
Date of DNA extraction
17 July 2014
25 March 2015
26 March 2015
Lamb oropharyngeal swabs cluster into two distinct community types
Oropharyngeal swabs were taken from 40 lambs. Using the Laplace approximation, it was found that the swabs could be partitioned into two separate groups based upon the types of bacteria present. These appeared to correspond to either oropharyngeal-type (partition 1) or rumen-type (partition 2) bacteria (Additional file 4). The oropharyngeal-type communities were dominated by the OTUs Pasteurellaceae (22%), Mannheimia (14%), Fusobacterium (11%), Bibersteinia trehalosi (8%), Neisseriaceae (7%), Moraxella (6%) and Bibersteinia (5%). The rumen-type communities were dominated by the OTUs Prevotella (36%), Clostridiales (11%), Ruminococcaceae (7%), Lachnospiraceae (6%) and Butyrivibrio (6%).
The richness (Chao 1: non-normal data) and diversity (inverse Simpsons: normal data) of the partitions were compared. There was no significant difference in richness or diversity between the rumen-type partition and the oropharyngeal-type partition.
Dichotomous oropharyngeal microbiota are associated with different lung community structures
We compared the proportions of the dominant OTUs in rumen-type swabs with their corresponding proportions in lung samples. On average, these OTUs were found in the following proportions in lung samples: Prevotella (5%), Clostridiales (2%), Ruminococcaceae (3%), Lachnospiraceae (1%) and Butyrivibrio (1%).
The presence of a lung-specific microbiota
OTUs significantly more indicative of lung fluids than oropharyngeal swabs
Average proportion in lung fluids (mean ± SD) (%)
Average proportion in oropharyngeal swabs (mean ± SD) (%)
Highest proportion in PBS controls (%)
1.0 ± 1.7
0.035 ± 0.10
1.2 ± 1.4
0.064 ± 0.24
1.9 ± 2.3
0.065 ± 0.19
0.80 ± 1.7
0 ± 0
0.65 ± 2.6
0.0063 ± 0.029
0.79 ± 1.4
0.077 ± 0.31
0.57 ± 1.4
0.0023 ± 0.0068
1.6 ± 2.1
0.040 ± 0.10
1.2 ± 2.8
0.0045 ± 0.013
0.77 ± 1.6
0.0080 ± 0.029
0.96 ± 1.5
0.012 ± 0.043
0.65 ± 1.1
0.00057 ± 0.0036
1.8 ± 2.2
0.050 ± 0.11
0.84 ± 2.3
0.0040 ± 0.020
0.51 ± 1.1
0 ± 0
0.58 ± 1.2
0.0057 ± 0.021
0.52 ± 1.2
0.0011 ± 0.0071
0.71 ± 1.3
0.0045 ± 0.018
0.53 ± 0.56
0 ± 0
3.9 ± 5.4
0.060 ± 0.18
13.3 ± 9.6
0.32 ± 0.97
6.4 ± 5.3
0.18 ± 0.59
2.0 ± 3.8
0.025 ± 0.096
1.0 ± 1.8
0.0074 ± 0.021
1.2 ± 1.2
0.024 ± 0.060
Sheep are commonly used as large animal models of the respiratory system due to the physiological and immunological similarities of their lungs to those of humans [10, 11, 31, 32]. We have previously used sheep to study both the extent of variation in the lung microbiota  and the direct and remote changes in the lung microbiota caused by localised P. aeruginosa infection and antibiotic treatment . As the sheep is an important agricultural animal, studies of its respiratory microbial communities may also be of interest from an animal health perspective.
It has previously been demonstrated that microaspiration of microbes from the upper aerodigestive tract is common in humans and can lead to an inflammatory phenotype . When microbial communities from healthy human lungs are characterised, they are often found to contain microbes associated with the upper aerodigestive tract . The healthy human lung microbiota is thought to be formed predominantly from the neutral dispersal of these upper aerodigestive tract microbes into the lungs rather than by the differential growth of lung-adapted microbial communities . We sought to identify whether this was also the case in sheep.
Sheep oropharyngeal swabs could be partitioned into two separate groups which were predominantly composed of OTUs identified as bacteria which are well known members of either the rumen (Prevotella, Clostridiales, Ruminococcaceae, Lachnospiraceae and Butyrivibrio [36–38]) or respiratory tract microbiotas (Pasteurellaceae, Mannheimia, Fusobacterium, Bibersteinia trehalosi, Neisseriaceae, Moraxella and Bibersteinia [39–41]). These bacteria were also detected in a previous study examining sheep buccal swabs .
It is not possible to identify whether this dichotomy reflected recent rumination or some stochastic post-mortem leakage of rumen fluid into the oropharynx in some individuals. The lambs during this study were not weaned but were at an age when it is expected that they all would be regularly supplementing their diet with grass and would be ruminating.
Regardless of the drivers of this oropharyngeal dichotomy, the microbial communities found in the lungs were very different to those found in both the rumen- and oropharyngeal-type swabs. A large number of bacterial OTUs were found to be significantly more abundant in lung fluids in comparison to oropharyngeal swabs, including Staphylococcus equorum which was by far the most common bacterial OTU found. Several OTUs which are commonly associated with the rumen were also identified in lung fluids. Our lung fluid samples will have been more affected by reagent contamination than the oropharyngeal swabs due to the lower quantity of bacterial DNA present . However, the microbial communities found in lung fluids did not reflect the bacteria found in reagent-only controls processed on the same day, so the presence of bacteria in the lamb lung cannot be attributed purely to sample contamination nor can it be attributed to disease as no lambs showed clinically overt signs of respiratory illness during the study.
There are several reasons why the microbes found in lamb lungs might not reflect those found in the upper aerodigestive tract to the same extent as is found in humans. Sheep have evolved to cope with rumination and thereby may have more efficient anatomical barriers to microaspiration . Physiological and anatomical differences such as the horizontal disposition of the lungs, increased nasal breathing and increased saliva production [45, 46] may also contribute.
In this study, we examined oropharyngeal swab and lung fluid samples taken from healthy lambs to characterise the bacterial communities present and to assess the impact of rumination on these communities. We found that the oropharyngeal swabs were dominated by either rumen-type or oropharyngeal-type microbial communities. We also found that lung bacteria did not greatly resemble either rumen- or oropharyngeal-type swabs and identified several bacterial OTUs which were more indicative of lung fluids. The lungs did contain several rumen-associated bacteria which may indicate that there is a certain degree of microaspiration of ruminal contents in lambs.
Sheep are not human, but the opportunities that they, and other large animals, present offer valuable insights into the dynamic relationship of the upper aerodigestive and lower respiratory tract microbiota in health. In the future, their value may extend to developing an understanding of the factors that predispose the upper aerodigestive tract microbiota towards disease in the lower respiratory tract.
The authors would like to thank the farm staff for the care of the animals used in this study and the staff and students from SRUC and the University of Edinburgh who aided during sample collection. The following reagent was obtained through BEI Resources, NIAID, NIH as part of the Human Microbiome Project: Genomic DNA from Microbial Mock Community B (Even, Low Concentration), v5.1L, for 16S rRNA Gene Sequencing, HM-782D.
This study was funded by the Scottish Government’s Rural and Environment Science and Analytical Services Division (RESAS) as part of the ‘Strategic Partnership for Animal Science Excellence’ programme. Laura Glendinning was funded by the BBSRC (Biotechnology and Biological Sciences Research Council) under the EASTBIO doctoral training programme (grant number: BB/J01446X/1). The funders had no role in the design of this microbiota study, the analysis and interpretation of the data or the writing of this manuscript.
Availability of data and materials
The datasets generated during and/or analysed during the current study are available in the Short Read Archive and can be accessed at Bioproject accession number PRJNA317719 (http://www.ncbi.nlm.nih.gov/bioproject/PRJNA317719).
LG contributed to study design, drafting of the paper and sample collection, processing and analysis. DC contributed to study design, drafting of the paper and sample collection. SW contributed to drafting of the paper and sample collection. KR contributed to study design and drafting of the paper. GM contributed to study design, drafting of the paper and sample collection. All authors read and approved the final manuscript.
Animals were used which were part of a study examining the effects of prenatal stress. This study was approved by Scotland’s Rural College’s (SRUC) Animal Experiments Committee and was conducted under Home Office licence.
Consent for publication
The authors declare that they have no competing interests.
Open AccessThis article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.
- Martinsen EMH, Leiten EO, Bakke PS, Eagan TML, Gronseth R. Participation in research bronchoscopy: a literature review. Eur Respir J. 2016;3:29511. https://doi.org/10.3402/ecrj.v3.29511.Google Scholar
- Yadava K, Pattaroni C, Sichelstiel AK, Trompette A, Gollwitzer ES, Salami O, et al. Microbiota promotes chronic pulmonary inflammation by enhancing IL-17A and autoantibodies. Am J Respir Crit Care Med. 2015;115(1):63–8. https://doi.org/10.1164/rccm.201504-0779OC.Google Scholar
- Yu W, Yuan X, Xu X, Ding R, Pang L, Liu Y, et al. Reduced airway microbiota diversity is associated with elevated allergic respiratory inflammation. Ann Allergy Asthma Immunol. 2015;115(1):63–8. https://doi.org/10.1016/j.anai.2015.04.025.View ArticlePubMedGoogle Scholar
- Stokes K, Wood A, August A. The microbiome regulates the eosinophil requirement for development of allergic lung inflammation. J Immunol. 2015;194(1 Supplement):53.18.Google Scholar
- Barfod KK, Roggenbuck M, Hansen LH, Schjørring S, Larsen ST, Sørensen SJ, et al. The murine lung microbiome in relation to the intestinal and vaginal bacterial communities. BMC Microbiol. 2013;13(1):1–12. https://doi.org/10.1186/1471-2180-13-303.View ArticleGoogle Scholar
- Yun Y, Srinivas G, Kuenzel S, Linnenbrink M, Alnahas S, Bruce KD, et al. Environmentally determined differences in the murine lung microbiota and their relation to alveolar architecture. PLoS One. 2014;9(12):e113466. https://doi.org/10.1371/journal.pone.0113466.View ArticlePubMedPubMed CentralGoogle Scholar
- Poroyko V, Meng F, Meliton A, Afonyushkin T, Ulanov A, Semenyuk E, et al. Alterations of lung microbiota in a mouse model of LPS-induced lung injury. Am J Physiol Lung Cell Mol Physiol. 2015;309(1):L76–83. https://doi.org/10.1152/ajplung.00061.2014.View ArticlePubMedPubMed CentralGoogle Scholar
- Barfod KK, Vrankx K, Mirsepasi-Lauridsen HC, Hansen JS, Hougaard KS, Larsen ST, et al. The murine lung microbiome changes during lung inflammation and intranasal vancomycin treatment. Open Microbiol J. 2015;9:167–79. https://doi.org/10.2174/1874285801509010167.View ArticlePubMedPubMed CentralGoogle Scholar
- Singanayagam A, Glanville N, Pearson R, James P, Cuthbertson L, Cox M, et al. T1 Fluticasone propionate alters the resident airway microbiota and impairs anti-viral and anti-bacterial immune responses in the airways. Thorax. 2015;70(Suppl 3):A1. https://doi.org/10.1136/thoraxjnl-2015-207770.1.View ArticleGoogle Scholar
- Bouljihad M, Leipold HW. An ultrastructural study of pulmonary bronchiolar and alveolar epithelium in sheep. J Vet Med A. 1994;41(8):573–86.View ArticleGoogle Scholar
- Griebel PJ, Entrican G, Rocchi M, Beskorwayne T, Davis WC. Cross-reactivity of mAbs to human CD antigens with sheep leukocytes. Vet Immunol Immunopathol. 2007;119(1–2):115–22. https://doi.org/10.1016/j.vetimm.2007.06.015.View ArticlePubMedGoogle Scholar
- Scheerlinck JPY, Snibson KJ, Bowles VM, Sutton P. Biomedical applications of sheep models: from asthma to vaccines. Trends Biotechnol. 2008;26(5):259–66. https://doi.org/10.1016/j.tibtech.2008.02.002.View ArticlePubMedGoogle Scholar
- Meeusen EN, Snibson KJ, Hirst SJ, Bischof RJ. Sheep as a model species for the study and treatment of human asthma and other respiratory diseases. Drug Discov Today Dis Models. 2009;6(4):101–6. https://doi.org/10.1016/j.ddmod.2009.12.002.View ArticleGoogle Scholar
- Collie D, Glendinning L, Govan J, Wright S, Thornton E, Tennant P, et al. Lung microbiota changes associated with chronic Pseudomonas aeruginosa lung infection and the impact of intravenous colistimethate sodium. PLoS One. 2015;10(11):e0142097. https://doi.org/10.1371/journal.pone.0142097.View ArticlePubMedPubMed CentralGoogle Scholar
- Glendinning L, Wright S, Pollock J, Tennant P, Collie D, McLachlan G. Variability of the sheep lung microbiota. Appl Environ Microbiol. 2016;82(11):3225–38. https://doi.org/10.1128/aem.00540-16.View ArticlePubMedPubMed CentralGoogle Scholar
- Dickson RP, Huffnagle GB. The lung microbiome: new principles for respiratory bacteriology in health and disease. PLoS Pathog. 2015;11(7):e1004923. https://doi.org/10.1371/journal.ppat.1004923.View ArticlePubMedPubMed CentralGoogle Scholar
- Bassis CM, Erb-Downward JR, Dickson RP, Freeman CM, Schmidt TM, Young VB, et al. Analysis of the upper respiratory tract microbiotas as the source of the lung and gastric microbiotas in healthy individuals. MBio. 2015;6(2):e00037–15. https://doi.org/10.1128/mBio.00037-15.View ArticlePubMedPubMed CentralGoogle Scholar
- Berry D, Ben Mahfoudh K, Wagner M, Loy A. Barcoded primers used in multiplex amplicon pyrosequencing bias amplification. Appl Environ Microbiol. 2011;77(21):7846–9. https://doi.org/10.1128/AEM.05220-11.View ArticlePubMedPubMed CentralGoogle Scholar
- Martin M. Cutadapt removes adapter sequences from high-throughput sequencing reads. Embnetjournal. 2011;17(1):10–2. 10.14806/ej.17.1.200.Google Scholar
- Schloss PD, Westcott SL, Ryabin T, Hall JR, Hartmann M, Hollister EB, et al. Introducing mothur: open-source, platform-independent, community-supported software for describing and comparing microbial communities. Appl Environ Microbiol. 2009;75(23):7537–41. https://doi.org/10.1128/AEM.01541-09.View ArticlePubMedPubMed CentralGoogle Scholar
- Kozich JJ, Westcott SL, Baxter NT, Highlander SK, Schloss PD. Development of a dual-index sequencing strategy and curation pipeline for analyzing amplicon sequence data on the MiSeq illumina sequencing platform. Appl Environ Microbiol. 2013;79(17):5112–20. https://doi.org/10.1128/AEM.01043-13.View ArticlePubMedPubMed CentralGoogle Scholar
- Esty WW. The efficiency of Good’s nonparametric coverage estimator. Ann Stat. 1986;14(3):1257–60.View ArticleGoogle Scholar
- Dufrêne M, Legendre P. Species assemblages and indicator species:the need for a flexible assymetrical approach. Ecol Monogr. 1997;67(3):345–66.Google Scholar
- Holmes I, Harris K, Quince C. Dirichlet multinomial mixtures: generative models for microbial metagenomics. PLoS One. 2012;7(2):e30126. https://doi.org/10.1371/journal.pone.0030126.View ArticlePubMedPubMed CentralGoogle Scholar
- R Core Team. R: a language and environment for statistical computing. https://www.R-project.org/: R Foundation for Statistical Computing; 2015.
- Oksanen J, Blanchet F, Kindt R, Legendre P, Minchin P, O'Hara R et al. Vegan: community ecology package. R package version 2.3-4/r2988 ed. http://R-Forge.R-project.org/projects/vegan/}2016.
- Neuwirth E. RColorBrewer: ColorBrewer palettes. R package version 1.1-2 ed. https://cran.r-project.org/web/packages/RColorBrewer/RColorBrewer.pdf.
- Warnes G, Bolker B, Bonebakker L, Gentleman R, Liaw W, Lumley T et al. gplots: various R programming tools for plotting data. R package version 2.17.0 ed. https://cran.r-project.org/web/packages/gplots/index.html.
- Ploner A. Heatplus: heatmaps with row and/or column covariates and colored clusters. R package version 2.16.0 ed. http://bioconductor.org/packages/release/bioc/html/Heatplus.html.
- Bray JR, Curtis JT. An ordination of the upland forest communities of southern Wisconsin. Ecol Monogr. 1957;27(4):326–49. https://doi.org/10.2307/1942268.View ArticleGoogle Scholar
- Entrican G, Wattegedera SR, Griffiths DJ. Exploiting ovine immunology to improve the relevance of biomedical models. Mol Immunol. 2015;66(1):68–77. https://doi.org/10.1016/j.molimm.2014.09.002.View ArticlePubMedPubMed CentralGoogle Scholar
- Enkhbaatar P, Nelson C, Salsbury JR, Carmical JR, Torres KEO, Herndon D, et al. Comparison of gene expression by sheep and human blood stimulated with the TLR4 agonists lipopolysaccharide and monophosphoryl lipid A. PLoS One. 2015;10(12):e0144345. https://doi.org/10.1371/journal.pone.0144345.View ArticlePubMedPubMed CentralGoogle Scholar
- Segal LN, Clemente JC, Tsay J-CJ, Koralov SB, Keller BC, Wu BG, et al. Enrichment of the lung microbiome with oral taxa is associated with lung inflammation of a Th17 phenotype. Nat Microbiol. 2016;1:16031. https://doi.org/10.1038/nmicrobiol.2016.31.View ArticlePubMedPubMed CentralGoogle Scholar
- Segal L, Alekseyenko A, Clemente J, Kulkarni R, Wu B, Chen H, et al. Enrichment of lung microbiome with supraglottic taxa is associated with increased pulmonary inflammation. Microbiome. 2013;1(1):19. https://doi.org/10.1186/2049-2618-1-19.View ArticlePubMedPubMed CentralGoogle Scholar
- Venkataraman A, Bassis CM, Beck JM, Young VB, Curtis JL, Huffnagle GB, et al. Application of a neutral community model to assess structuring of the human lung microbiome. MBio. 2015;6(1):e0084–14. https://doi.org/10.1128/mBio.02284-14.View ArticleGoogle Scholar
- Huws SA, Kim EJ, Lee MRF, Scott MB, Tweed JKS, Pinloche E, et al. As yet uncultured bacteria phylogenetically classified as Prevotella, Lachnospiraceae incertae sedis and unclassified Bacteroidales, Clostridiales and Ruminococcaceae may play a predominant role in ruminal biohydrogenation. Environ Microbiol. 2011;13(6):1500–12. https://doi.org/10.1111/j.1462-2920.2011.02452.x.View ArticlePubMedGoogle Scholar
- Zened A, Combes S, Cauquil L, Mariette J, Klopp C, Bouchez O, et al. Microbial ecology of the rumen evaluated by 454 GS FLX pyrosequencing is affected by starch and oil supplementation of diets. FEMS Microbiol Ecol. 2013;83(2):504–14. https://doi.org/10.1111/1574-6941.12011.View ArticlePubMedGoogle Scholar
- Omoniyi LA, Jewell KA, Isah OA, Neumann AP, Onwuka CFI, Onagbesan OM, et al. An analysis of the ruminal bacterial microbiota in West African Dwarf sheep fed grass- and tree-based diets. J Appl Microbiol. 2014;116(5):1094–105. https://doi.org/10.1111/jam.12450.View ArticlePubMedGoogle Scholar
- Christensen H, Kuhnert P, Nørskov-Lauritsen N, Planet PJ, Bisgaard M. The family Pasteurellaceae. In: Rosenberg E, EF DL, Lory S, Stackebrandt E, Thompson F, editors. The prokaryotes: Gammaproteobacteria. Berlin, Heidelberg: Springer Berlin Heidelberg; 2014. p. 535–64.Google Scholar
- Ling Z, Liu X, Wang Y, Li L, Xiang C. Pyrosequencing analysis of the salivary microbiota of healthy Chinese children and adults. Microb Ecol. 2012;65(2):487–95. https://doi.org/10.1007/s00248-012-0123-x.View ArticlePubMedGoogle Scholar
- Cardenas PA, Cooper PJ, Cox MJ, Chico M, Arias C, Moffatt MF, et al. Upper airways microbiota in antibiotic-naïve wheezing and healthy infants from the tropics of rural Ecuador. PLoS One. 2012;7(10):e46803. https://doi.org/10.1371/journal.pone.0046803.View ArticlePubMedPubMed CentralGoogle Scholar
- Kittelmann S, Kirk MR, Jonker A, McCulloch A, Janssen PH. Buccal swabbing as a noninvasive method to determine bacterial, archaeal, and eukaryotic microbial community structures in the rumen. Appl Environ Microbiol. 2015;81(21):7470–83. https://doi.org/10.1128/AEM.02385-15.View ArticlePubMedPubMed CentralGoogle Scholar
- Salter SJ, Cox MJ, Turek EM, Calus ST, Cookson WO, Moffatt MF, et al. Reagent and laboratory contamination can critically impact sequence-based microbiome analyses. BMC Biol. 2014;12:87. https://doi.org/10.1186/s12915-014-0087-z.View ArticlePubMedPubMed CentralGoogle Scholar
- May NDS. The anatomy of the sheep. St. Lucia: University of Queensland Press; 1964.Google Scholar
- Hungate RE. Chapter V The rumen as a continuous fermentation system. The rumen and its microbes. NY, USA: Academic Press; 1966. p. 206–44.Google Scholar
- Sarkar A, Singh H. Oral behaviour of food emulsions. Food oral processing. Wiley-Blackwell; 2012. p. 111-137. doi: 10.1002/9781444360943.