Rumen microbial community composition varies with diet and host, but a core microbiome is found across a wide geographical range

<jats:title>Abstract</jats:title><jats:p>Ruminant livestock are important sources of human food and global greenhouse gas emissions. Feed degradation and methane formation by ruminants rely on metabolic interactions between rumen microbes and affect ruminant productivity. Rumen and camelid foregut microbial community composition was determined in 742 samples from 32 animal species and 35 countries, to estimate if this was influenced by diet, host species, or geography. Similar bacteria and archaea dominated in nearly all samples, while protozoal communities were more variable. The dominant bacteria are poorly characterised, but the methanogenic archaea are better known and highly conserved across the world. This universality and limited diversity could make it possible to mitigate methane emissions by developing strategies that target the few dominant methanogens. Differences in microbial community compositions were predominantly attributable to diet, with the host being less influential. There were few strong co-occurrence patterns between microbes, suggesting that major metabolic interactions are non-selective rather than specific.</jats:p>

section_content: Rumen microbes can be assigned to different functional groups, such as cellulolytics, amylolytics, proteolytics, etc., which degrade the wide variety of feed components or further metabolize some of the products formed by other microbes 3. For example, methanogens, the methane-forming archaea, are among those that metabolize hydrogen formed by some fermentative microbes to form methane. The methane generated during this fermentation contributes to global anthropogenic greenhouse gas emissions 4 and represents a 2-12% loss of feed energy for the animal 5. Differences in rumen microbial communities underlie variations in methane formation 6 and the conversion of feed to animal products 7, 8. Therefore, understanding these communities is key to understanding ruminal transformations of plant material to both undesirable and useful ruminant products. 

The aim of this study was to determine the composition of the microbiota in rumen and foregut samples from 742 individual animals from around the world. The resulting dataset allowed us to determine that dietary factors dominate over host species in determining microbial community composition, identify the dominant microbes and their potential associations, and describe the degree of similarity of rumen microbial communities worldwide.

section_content: This is the largest single study to examine microbial communities across a range of ruminant and camelid species, diets, and geographical regions. A standardised pipeline was used to process samples in order to minimise variation introduced by processing steps such as DNA extraction or PCR amplification. This is important for detecting authentic patterns rather than ones introduced by methodological differences between different studies 9. The primers chosen amplify, to the best of our current knowledge, the target gene regions from nearly all known bacteria, archaea, and rumen ciliates.
section_name: Results and Discussion

section_content: Despite the range of ruminants with different feeding strategies and diets, similar rumen bacteria were abundant around the world (Fig. 1 ). There was some variation in bacterial community compositions in animals from different regions, likely to be caused by differences in diet, climate, and farming practices. The 30 most abundant bacterial groups (Greengenes 10 taxonomy summarised at the genus-level) were all found in over 90% of samples, and together comprised 89. 4% of all sequence data (Supplementary Table 1 ) and were similar to those described in an earlier meta-analysis of rumen microbial communities 11. All 30 are known rumen-inhabiting bacteria. Because the samples came from a wide range of ruminant species, diets, and geographical locations, these data suggest that new dominant bacteria are not likely to be found in future studies. The seven most abundant bacterial PowerPoint and Adobe Illustrator. https://creativecommons. org/licenses/by-sa/3. 0/deed. en groups comprised 67. 1% of all bacterial sequence data, were detected in all samples (Supplementary Fig. 1 ), and can be considered the "dominant" rumen bacteria. They were Prevotella, Butyrivibrio, and Ruminococcus, as well as unclassified Lachnospiraceae, Ruminococcaceae, Bacteroidales, and Clostridiales. These might be considered a "core bacterial microbiome" at the genus level or higher, because they are present in a large selection of ruminants, so confirming the suggestion that there is a core rumen microbiome 9. However, these bacterial groups were not equally abundant in all animal species (P ≤ 0. 005; Supplementary Table 2 ). With the exception of Butyrivibrio 12, these groups are not adequately represented by characterised cultures 13, and their functions are not well understood. 

Inspection of the most abundant and prevalent bacterial operational taxonomic units (OTUs) in the dataset showed that only 14% fell within a named species, and 70% were not even within a formally recognised genus (Fig. 2a ). When cultured isolates from as-yet unnamed species were included in the analysis, the dominant OTUs were better (35%) but still poorly represented by cultures that belonged to potentially the same species (Fig. 2b ). This study shows that, while we appear to recognize the dominant rumen bacteria, considerable microbiological effort is still required to understand them. Some efforts have been made to isolate more cultures and gather more information about these bacteria 13, 14. For example, the genomes of Prevotella aff. ruminicola Tc2-24, rumen bacterium R-7, and other isolates whose 16S rRNA gene sequences are similar to those of dominant rumen bacterial OTUs (Fig. 2b ), have been sequenced as part of the Hungate1000 project 15. 

Because there is a flux of both liquids and solids through the rumen 16, microbes must actively metabolize to gain energy and multiply to counteract washout and so maintain populations in the rumen 17. The dominant bacteria found in this study are therefore likely to be responsible for the majority of the 8 and 9 8 and 9 ). Bacterial OTU abundances were multiplied by a factor of 15 relative to archaeal OTUs. Mbb. Methanobrevibacter. transformation of ingested feed in the rumen and camelid foregut, especially of cellulose, hemicellulose, pectin, starch, fructan, organic acids, and protein, as these are the major energy-yielding substrates used for microbial growth 17. There is also a convergence of bacterial community structure in the rumen and in the crop of the hoatzin, a bird that relies on a foregut fermentation of ingested leaves 18. Thus microbial community structure seems to be driven by the similarity of organ function extending across the rumen, the camelid foregut, and the crop of this unusual bird. More efforts should go into characterizing the metabolism and roles of these bacteria that are the responsible for the majority of feed fermentation, with the aim of enhancing animal productivity and reducing methane emissions. 

Nearly all the archaea were identified as methanogens known to be residents of the rumen (Supplementary Table 3, Supplementary Text 1), and their relative abundances were comparable to previous studies 19. The dominant archaeal groups were remarkably similar in all regions of the world (Fig. 1 ). This universality and limited diversity was also recently noted in survey of archaea in New Zealand ruminants 20 and could make it possible to successfully mitigate methane emissions by developing strategies, such as vaccines or small-molecule inhibitors, that target the few dominant methanogens. 

Members of the Methanobrevibacter gottschalkii and Methanobrevibacter ruminantium clades were found in almost all samples, and were the two largest groups, accounting for 74% of all archaea. Together with a Methanosphaera sp. and two Methanomassiliicoccaceae-affiliated groups, the five dominant methanogen groups comprised 89. 2% of the archaeal communities (Supplementary Fig. 1 ), showing that rumen archaea are much less diverse than rumen bacteria. This likely reflects the narrow range of substrates they use. Methanomicrobium has previously been reported as abundant in ruminants in Asia 19. In our study, they were found to comprise > 5% of the archaeal community of some Australian, Brazilian, Chinese, North American, and South African cattle, as well as South African sheep, showing them to be widely distributed, but not universally prevalent. The five dominant methanogen groups were not equally abundant in all animal species groups (P ≤ 0. 005; Supplemental Table 4 ). In contrast to bacteria, the rumen archaea are better represented by cultures, with 58% of the most abundant and prevalent OTUs falling within a named species, and all but 22% within named genera (Fig. 2c ). All of the latter were members of Methanomassiliicoccales, which is an order of relatively poorly-characterised methanogens 21 for which representative cultures of as-yet unnamed species and genera are available (Fig. 2d ) 22. The 50 most abundant OTUs accounted for 74. 5% of the archaeal sequence data, again indicating a much lower diversity than in the bacteria, where the 50 most abundant OTUs made up only 11. 0%. 

By assigning physiologies (Supplementary Table 5 ) to the sequence abundance information (Supplementary Table 3 ), it can be concluded that 77. 7% of archaea were hydrogenotrophic methanogens, while 22. 1% had the ability to grow with hydrogen plus methyl groups derived from methanol or methylamines. Methanogens able to form methane from acetate (Methanosarcina spp. and Methanosaeta spp. ) were extremely rare (< 0. 015%; Supplemental Data 1), as expected based on their general slow growth rates that would not allow them to be maintained in the rumen under normal conditions. 

Almost all protozoal sequence data (> 99. 9%) were assigned to 12 genus-equivalent protozoal groups (Supplementary Table 6 ). It was apparent that the variability of protozoa between and within cohorts of co-located animals was much greater than that of bacteria and archaea (Supplementary Fig. 2 ). It has been reported that there is strong host individuality of rumen protozoal community structure 9, and this is evident in our study. The genera Entodinium and Epidinium dominated, occurring in more than 90% of samples and representing 54. 7% of protozoal sequence data (Supplementary Fig. 1 ). Many of the protozoal genera were present in greater than 70% of the samples, indicating a wide prevalence. Genera such as Enoploplastron and Ophryoscolex had a wider than expected host distribution. They are considered to be mainly present in sheep and cattle, respectively 23, but we also found Enoploplastron in cattle, deer, and reindeer samples from twelve countries, and Ophryoscolex in buffalo, goats, deer, sheep, and giraffe samples from 18 countries. Although different rumen protozoa are reported to have limited host and geographical distributions, host specificity has been questioned 24. It seems likely that further investigation will demonstrate greater ubiquity of the rumen protozoa.
section_name: Dominant rumen microbes.

section_content: Because the abundance of microbial groups varied between animal species groups and cohorts (Supplementary Tables 2 and 4 ; Supplementary Fig. 2 ), we looked for factors that might underlie this. Rumen and camelid foregut microbial community structure could be expected to be shaped by morphological, physiological, and even behavioural characteristics that evolved along with the varied feeding strategies in the various ruminant lineages 2. Indeed, adaptation has resulted in a diversity of rumen sizes and passage rates of rumen contents, allowing ruminant species to exploit a range of feed types. In addition to feed composition effects 25, these host adaptations might also play a role in regulating rumen microbial community structure. Because our dataset was from ruminants and camelids from different lineages consuming a range of diets, host and diet effects on rumen microbial community structure could be separated. 

To look at diet and host effects, we classified the diets based on forage and browse or concentrate content (Supplementary Table 7 ) and grouped the animals according to their lineage (Supplementary Data 1). Microbial communities could clearly be discriminated by both host and diet (Fig. 3a ), with bacteria being the main drivers behind the observed differences (Fig. 3b ). This probably reflects their more diverse metabolic capabilities compared with the less versatile archaea and protozoa. We investigated the patterns of microbial abundances across hosts and diets (Fig. 3c, Supplementary Fig. 3 4 5 6 ). Ruminococcus, one of the dominant bacteria, was relatively evenly distributed, but this was an exception. For many bacteria, diet was the major factor determining relative abundance. Bacterial communities from forage-fed animals were similar to each other, those from concentrate-fed animals were similar to each other, but distinct from those in forage-fed animals, and those from animals fed mixed diets were intermediate between these. Unclassified Bacteroidales and Ruminococcaceae were more abundant in all animals fed forages. Some as-yet poorly characterised Bacteroidales are postulated to be able to degrade cellulose, and their genomes encode a broad range of plant polysaccharide degrading capabilities 26, 27, which could explain their pattern of distribution. In contrast, members of Prevotella and unclassified Succinivibrionaceae were more abundant in animals fed diets containing concentrate. Based on the physiologies of cultured relatives 28, 29, these are probably major producers of propionate and the propionate-precursor succinate, and so are responsible for the greater levels of propionate formed from concentrate-rich diets 25. The abundance of only a few other major bacterial groups was associated with host lineage (Fig. 3c ). For example, unclassified Veillonellaceae were proportionally more abundant in sheep, deer, and camelids (Fig. 3d ). This may be related to differences in rumen and camelid foregut sizes, anatomy, and feeding frequencies compared to bovines 2. 

The relative abundances of several major bacterial groups were affected by both host and diet (Fig. 3c ). Unclassified Clostridiales were most abundant in bovines fed forage and least abundant in bovines fed high concentrate diets, while in caprids, cervids, and camelids these diet differences were far less pronounced. Butyrivibrio was most abundant in rumen samples from bovines fed mixes of forage and concentrates. Fibrobacter was most abundant in bovines fed forage. When concentrate was included in bovine diets, the relative abundance of Fibrobacter was decreased, but it was still more abundant than in other animals. To examine its distribution in more detail, we compared Fibrobacter abundances across different ruminant species and found significantly higher levels in bovines compared to deer, sheep, or camelids (Fig. 3e ). These data suggest that Fibrobacter is favoured in the bovine rumen and, given that it is cellulose degrader 30, may play an essential role in the degradation of plant fibre in cattle. 

Overall, diet was a major determinant of bacterial community structure. This may be because physical and chemical characteristics of the feed determine the different microbial niches available. In contrast to the post-gastric mammalian digestive tract 31, and due to the sheer volume of digesta and feed input, there is probably less shaping of the rumen microbial community by local host biological factors such as the immune system, secreted antimicrobial peptides, host-cell glycosylation, and host-derived nutrients. 

Associations between rumen microbes. The abundance patterns within bacterial, archaeal, and protozoal communities in different hosts fed different diets showed that certain microbes exhibited parallel patterns of relative abundance (Fig. 3c, Supplementary Figs 2 3 4 5 6 ). We therefore looked for correlations within and between bacteria, archaea, and protozoa (Fig. 4 and Supplementary Fig. 7 ), reasoning that specific associations should be seen across diets, hosts, and geography. Negative correlations of abundances of groups were observed within the bacteria, archaea, and protozoa, including replacement effects between dominant groups within each of these (Supplementary Text 1 and Supplementary Fig. 7 ). Few strong positive correlations were found within bacteria, archaea, and protozoa. For example, there was a strong correlation between Veillonellaceae and the TG5 group, driven by their co-occurrence within cervids and caprids. These microbes may cooperate in the rumen, or they may share similar requirements and so certain hosts and diets would offer better opportunities for their growth. This explanation could also underlie the strong positive correlations observed between different groups of methylotrophic methanogens (Supplementary Fig. 7 ). They may be responding to diets rich in methyl groups, such as feeds with high levels of pectins or osmolytes such as betaine. The strongest correlation within protozoa was a positive one between Dasytricha and Isotricha. These two genera of holotrichous protozoa display very similar spectra of substrate use, including use of plant soluble sugars and storage carbohydrates 24, again suggesting that co-occurrence may be due to exploitation of similar opportunities. 

We also investigated associations between bacteria, archaea, and protozoa. Strikingly, no strong correlations were detected between archaea and protozoa (Supplementary Fig. 7 ). Methanogens are known to colonize protozoa, and this mutualistic relationship is believed to enhance methane formation in the rumen 32. The occurrence of specific symbioses between methanogens and rumen protozoa has been speculated on, but not convincingly demonstrated 33. The lack of strong co-occurrence patterns within this study indicates that these undoubtedly important associations are probably non-specific, or occur at a strain level. Further investigation is required to corroborate this interesting finding, as mechanisms that mediate the colonization of protozoa by archaea remain to be elucidated. These could have interesting evolutionary aspects if they allow non-specific interactions to form or are mediated by strain-specific mechanisms that confer different partner specificities within archaeal or protozoal species. In contrast, there were some positive associations between bacterial and protozoal groups. Most noticeable were the associations of Isotricha and Dasytricha with Fibrobacter. Fibrobacter were reported to decrease in abundance in animals where protozoa were eliminated 34, indicating that there may be a mutually beneficial relationship between these protozoa and Fibrobacter, which are surface colonizers of plant material 23. 

No strong associations were found between the most abundant bacteria and archaea (Fig. 4 ). This was surprising, since rumen bacteria degrade feed and produce the substrates that methanogens use for growth, mainly hydrogen and methyl groups. In contrast, there were distinct positive associations between some less abundant bacteria and archaea. The strongest association was between bacteria such as the succinate-producing Succinivibrionaceae, the succinate-using Dialister, and the amino-acid-fermenting Acidaminococcus, and methanogens belonging to the Methanomassiliicoccaceae, Methanosphaera sp. A4, and Methanobrevibacter boviskoreani. Succinivibrio spp. degrade pectin 28, and methanol is required for growth of Methanomassiliicoccaceae 35 and Methanosphaera 36, explaining part of this pattern. Other associations were between the methylotrophic methanogen Methanosphaera sp. ISO3-F5 and different bacteria, including members of Lachnospiraceae. These associations may be based on the ability of Lachnospiraceae to degrade pectin and so provide methanol as a substrate for the methylotrophs 37. The associations between other Methanomassiliicoccaceae groups and various unclassified members of Bacteroidales suggest the possibility of yet further methanol-dependent metabolic interactions. In contrast to archaeal-protozoal interactions, these findings suggest that some archaeal-bacterial interactions are specific, inferring specialised mechanisms for partner recognition or very similar requirements for growth. The basis for these associations remains to be determined. However, the general lack of strong association patterns between protozoa and the major bacteria on the one hand, and the major methanogen groups on the other, suggests that conserved mechanisms may mediate the interactions between hydrogen producing and hydrogen consuming microbes, allowing flexible interactions. This may aid methane mitigation research, since interfering with these potentially universal mechanisms could slow the rate of hydrogen transfer and so slow methane formation 38. It may also be that the interactions mainly occur via pools of common metabolites, especially where the end products of one group form the substrates of another. 

The results of this survey showed that the rumen microbial ecosystem is dominated by a core community composed of poorly-characterised microbes, especially amongst the bacteria. Diet had more influence than animal species on rumen or camelid foregut microbial community composition. Rumen ecosystems are typified by strong metabolic interactions between microbes that facilitate the fermentation of plant material to products useful for both the host and other rumen microbes 3, 17, 25, 32. The relatively few co-occurrence patterns seen in this study suggest that these microbial interactions do not rely on exclusive associations, and could indicate considerable promiscuity between members of interacting functional groups. Analysis at metagenomic and metatranscriptomic levels could in future uncover whether common functional elements that facilitate interactions are shared among multiple species. It seems plausible that functional redundancy among the microbes 9 means that multiple microbial species can fulfil the same function, with different combinations of microbes being co-selected depending on 1 and 3 ). Mbb. Methanobrevibacter, Mmc. Methanomassiliicoccales, *indicates unclassified bacteria within a family. 

the diet. This flexibility of rumen microbial community structure would confer on the ruminant host the ability to exploit a variety of different plant feeds.
section_name: Effects of diet and host on microbial community composition.

section_content: Geographical distribution and diversity of gastrointestinal tract content samples. A total of 742 samples from 32 species or sub-species of ruminants and other foregut fermenters in 35 countries and seven global regions were selected for sequencing of microbial marker genes (Fig. 1, Supplementary Data 1). The samples were from cattle, bison, and buffalo (bovines), sheep and goats (caprids), deer (cervids), and alpacas, llamas, and guanacos (camelids), including diverse breeds of domestic cattle, sheep, and goats, and were largely made up of small cohorts of four or more co-located individuals consuming the same diet. We included foregut samples of camelids in this study, recognizing that these organs have a common function but evolved separately 39. The use of animals, including welfare, husbandry, experimental procedures, and the collection of samples used for this study, was, where applicable, approved by named institutional and/or licensing committees and performed in accordance with approved institutional and regulatory guidelines (please refer to Supplementary Data 1 for details of these).
section_name: Methods

section_content: To minimise variation introduced by differing methodologies, such as choice of sampling or DNA extraction method 40 and primer-driven gene amplification biases 41, we used a standardised pipeline to process samples (unless indicated otherwise in Supplementary Data 1). Briefly, approximately 20 g of whole (i. e., solid and liquid) mid-rumen or camelid foregut contents were collected via stomach tube, cannula, or post mortem as previously described 35. Samples were immediately frozen, freeze-dried, and then couriered to AgResearch. Freeze-dried samples were homogenised in a coffee blender and DNA was extracted from a representative 30 mg subsample using the PCQI method 40, 42. We assessed the structure of microbial communities by sequencing regions of bacterial and archaeal 16S rRNA genes and ciliate protozoal 18S rRNA genes in triplicate as described previously 35, 37 using primers comprised of (5′ to 3′ ) a sequencing adapter (A or B), a sample-unique 12-base error-correcting Golay barcode on one of each primer pair, a two-base linker, and a group-specific sequence targeting the marker gene. For bacteria, the primers were Ba515Rmod1 (adapter A-barcode-GT-CCGCGGCKGCTGGCAC) and Ba9F (adapter B-AC-GAGTTTGATCMTGGCTCAG). For archaea, the primers were Ar915aF (adapter A-barcode-GT-AGGAATTGGCGGGGGAGCAC) and Ar1386R (adapter B-CA-GCGGTGTGTGCAAGGAGC). For protozoa, the primers were Reg1320R (adapter A-barcode-TC-AATTGCAAAGATCTATCCC) and RP841F (adapter B-AA-GACTAGGGATTGGARTGG). Linker A was CCATCTCATCCCTGCGTGTCTCCGACTCAG and linker B was CCTATCCCCTGTGTGCCTTGGCAGTCTCAG. Amplicons were sequenced using 454 GS FLX Titanium chemistry at Eurofins MWG Operon (Ebersberg, Germany). Sample processing and pipeline reproducibility controls were performed to identify variation introduced during sample processing (Supplementary Text 1). Sequence data are available from GenBank [accession numbers PRJNA272135, PRJNA272136, and PRJNA273417]. 

Phylogenetic analysis of sequencing data. Pyrosequence data were processed and analysed using the QIIME software package version 1. 8 43. Sequences over 400 bp in length with an average quality score over 25 were assigned to a specific sample via the barcodes. The number of bacterial, archaeal, and ciliate protozoal sequencing reads available for analysis are summarised in Supplementary Data 1. Sequence data were grouped into operational taxonomic units (OTUs) sharing over 97% (bacteria -UCLUST 44 ), 99% (archaea -UCLUST) or 100% (ciliate protozoa -prefix_suffix option in QIIME) sequence similarity. Sequences were assigned to phylogenetic groups by BLAST 45. Bacterial 16S rRNA genes were assigned using the Greengenes database version 13_5 10, archaeal 16S rRNA genes using RIM-DB version 13_11_13 22 and ciliate protozoal 18S rRNA genes against an in-house database 46. Bacterial and ciliate protozoal data were summarised at the genus level. Archaea were summarised at the species level. Samples for which low read numbers were obtained or that contained high proportions of sequences from "exogenous" bacteria (i. e., likely environmental contaminants such as Stenotrophomonas) were excluded from further analyses (Supplementary Text 1). 

The identity of the most abundant and prevalent OTUs was determined using BLAST 45 against sequences from type material and against all sequences (excluding sequences from model organisms or environmental samples) in the nt database 47. Bellerophon (version 3, 200 bp window, Huber-Hugenholtz correction 48 ) was used to identify chimeric OTU sequences. Sequence similarities greater than 97% and 93% were used as cut-offs to classify OTUs at species-and genus level, respectively. The rationale for these cut-offs was discussed by Kenters et al. 49.
section_name: Sample collection, DNA extraction, amplification and processing of samples for high-throughput sequencing.

section_content: The range of diets consumed by the animals from which the samples came was highly diverse and complex. For this reason, and where the information was available, diets were categorised in terms of forage type, forage plant, and forage to concentrate ratio (Supplementary Table 7 ). Diets likely to contain > 5% starch (e. g., whole or grain crops of maize, barley, wheat, rice, as well as pea, potato, sorghum, etc. ) or > 5% pectin (e. g., beets or legumes such as alfalfa and clover) were also identified. Animals that had been fed their respective diets for less than a two-week period were noted in Supplementary Data 1. Factors such as gender, age, modifications (e. g., cannulation), treatments (e. g., antibiotics, drench, surgery), farming conditions, season, contact with other animals, and sample processing steps that may affect apparent microbial community compositions (e. g., DNA extraction method, sample fraction used, sample storage, etc. ) were also recorded (Supplementary Data 1). Where details were not provided, latitude, longitude, and elevation were estimated using http://www. mapcoordinates. net/en. Climate zones were designated according to the Köppen-Geiger climate classification scheme 50.
section_name: Simplified classification of dietary information and other factors.

section_content: The resulting dataset allowed us to establish whether animal or dietary factors relate to rumen and camelid foregut microbial community composition, identify the dominant microbes and their potential associations, and describe the degree of similarity of rumen and camelid foregut microbial communities worldwide. Statistical analyses of microbial data were performed using GenStat for Windows 51, R software 52, and QIIME 43. Principal coordinate analysis of Bray-Curtis dissimilarity matrices, analysis of variance, sparse partial least squares discriminant analysis (sPLS-DA, using a sPLS regression approach), and canonical discriminant analyses (CDA) of microbial community composition data in context of the metadata (Supplementary Data 1) were used to identify impacts of factors such as host lineage, diet, etc. on rumen and camelid foregut microbial communities and to identify the groups associated with these factors. Pearson, Spearman, SparCC 53, and regularised canonical correlation analyses (CCA) were used to identify associations within and between archaeal, bacterial, and protozoal groups. Association scores were visualised as relevance networks and clustered image maps (CIM, heatmaps) representing the first two dimensions. González et al. provides a comprehensive overview of sPLS-DA, CCA and the corresponding 'pairwise associations', network and CIM techniques and their application 54. 

AgResearch Limited, Grasslands Research Centre, Palmerston North 4442, New Zealand In addition, there was a typographical error in Affiliation 39 which was incorrectly listed as 'Department of Animal Science & Technology, Sunchon National University, Suncheon, Jeonnam 540-743, Korea'. The correct affiliation is listed below: Department of Animal Science & Technology, Sunchon National University, Suncheon, Jeonnam 540-742, Korea Lastly, there was a typographical error in the project number in the Acknowledgements section. This work is licensed under a Creative Commons Attribution 4. 0 International License. The images or other third party material in this article are included in the article's Creative Commons license, unless indicated otherwise in the credit line; if the material is not included under the Creative Commons license, users will need to obtain permission from the license holder to reproduce the material. To view a copy of this license, visit http://creativecommons. org/licenses/by/4. 0/
section_name: Statistical analyses.

section_name: Acknowledgements

section_content: We thank Ron Ronimus, Paul Newton, and Christina Moon for reading and commenting on the manuscript. We thank all who provided assistance that allowed Global Rumen Census collaborators to supply samples and metadata (Supplemental Text 1). AgResearch was funded by the New Zealand Government as part of its support for the Global Research Alliance on Agricultural Greenhouse Gases. The following funding sources allowed Global Rumen Census collaborators to supply samples and metadata, listed with the primary contact(s) for each funding source: Agencia Nacional de Investigación e Innovación, Martín Fraga ; Alberta Livestock and Meat Agency, Canada, Tim A. McAllister ; Area de Ciencia y Técnica, Universidad Juan A Maza (Resolución Proy. N° 508/2012 ), Diego Javier Grilli ; Canada British Columbia Ranching Task Force Funding Initiative, John Church ; CNPq, Hilário Cuquetto Mantovani, Luiz Gustavo Ribeiro Pereira ; FAPEMIG, Hilário Cuquetto Mantovani ; FAPEMIG, PECUS RumenGases, Luiz Gustavo Ribeiro Pereira ; Cooperative Research Program for Agriculture Science & Technology Development (project number PJ010906 ), Rural Development Administration, Republic of Korea, Sang-Suk Lee ; Dutch Dairy Board & Product Board Animal Feed, André Bannink, Kasper Dieho, Jan Dijkstra ; Ferdowsi University of Mashhad, Vahideh Heidarian Miri ; Finnish Ministry of Agriculture and Forestry, Ilma Tapio ; Instituto Nacional de Tecnología Agropecuaria, Argentina (Project PNBIO1431044 ), Silvio Cravero, María Cerón Cucchi ; Irish Department of Agriculture, Fisheries and Food, Alexandre B. De Menezes ; Meat & Livestock Australia ; and Department of Agriculture, Fisheries & Forestry (Australian Government ), Chris McSweeney ; Ministerio de Agricultura y desarrollo sostenible (Colombia), Olga Lucía Mayorga ; Montana Agricultural Experiment Station project ( MONB00113 ), Carl Yeoman ; Multistate project W-3177 Enhancing the competitiveness of US beef ( MONB00195 ), Carl Yeoman ; NSW Stud Merino Breeders' Association, Alexandre Vieira Chaves; Queensland Enteric Methane Hub, Diane Ouwerkerk; RuminOmics, Jan Kopecny, Ilma Tapio ; Rural and Environment Science and Analytical Services Division (RESAS) of the Scottish Government and the Technology Strategy Board, UK, R. John Wallace ; Science Foundation Ireland ( 09/RFP/GEN2447 ), Sinead Waters ; Secretaría de Agricultura, Ganadería, Desarrollo Rural, Pesca y Alimentación, Mario A. Cobos-Peralta ; Slovenian Research Agency (project number J1-6732 and P4-0097 ), Blaz Stres ; Strategic Priority Research
section_name: Acknowledgements

section_name: Author Contributions

section_content: Program, Climate Change: Carbon Budget and Relevant Issues (Grant No. XDA05020700), ZhiLiang Tan; The European Research Commission Starting Grant Fellowship (336355-MicroDE), Phil B. Pope; The Independent Danish Research Council (project number 4002-00036), Torsten Nygaard Kristensen; and The Independent Danish Research Council (Technology and Production, project number 11-105913), Jan Lassen. These funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.

section_content: The study was designed by GH and PHJ. Sample processing was carried out by FC and GH. Data analysis was done by GH, PHJ, SG, AJ, and WY. The manuscript was written by GH and PHJ, assisted by all coauthors. Global Rumen Census Collaborators performed animal trials, collected and provided samples and associated metadata used in this study, critically revised the manuscript, contributed intellectual content, and approved the manuscript for publication.
section_name: Author Contributions

section_content: Supplementary information accompanies this paper at http://www. nature. com/srep
section_name: Additional Information

section_content: The AgResearch component of this study was funded by the New Zealand Government via the Ministry for Primary Industries (MPI) as part of MPI's support for the Global Research Alliance on Agricultural Greenhouse Gases. The publication of the data reported here is at the discretion of MPI. MPI did not control which data were presented or how these data were interpreted within this paper. This does not alter the authors' adherence to all the Scientific Reports policies on sharing data and materials. Material transfer agreements, limiting the use of samples to this study, are in place between AgResearch and Global Rumen Census Collaborators from The University of Alberta (Canada), The Department of Agriculture, Fisheries and Forestry (Queensland, Australia), The University of Aberdeen (Scotland), and The National Institute of Livestock and Grassland Science (Japan). There are no patents, products in development or marketed products to declare. No competing interests were declared by Global Rumen Census Collaborators.
section_name: Competing financial interests: