Microbiome and Human Health: Current Understanding, Engineering, and Enabling Technologies.

The human microbiome is composed of a collection of dynamic microbial communities that inhabit various anatomical locations in the body. Accordingly, the coevolution of the microbiome with the host has resulted in these communities playing a profound role in promoting human health. Consequently, perturbations in the human microbiome can cause or exacerbate several diseases. In this Review, we present our current understanding of the relationship between human health and disease development, focusing on the microbiomes found across the digestive, respiratory, urinary, and reproductive systems as well as the skin. We further discuss various strategies by which the composition and function of the human microbiome can be modulated to exert a therapeutic effect on the host. Finally, we examine technologies such as multiomics approaches and cellular reprogramming of microbes that can enable significant advancements in microbiome research and engineering.

Proposal for reclassification of the species Hungatella xylanolytica as Lacrimispora xylanisolvens nom. nov. and transfer of the genus Hungatella to the family Lachnospiraceae.

Hungatella xylanolytica X5-1T is an anaerobic, xylan-fermenting bacterium first isolated from methane-producing cattle manure. Initially identified as Bacteroides xylanolyticus, this species was later reclassified as H. xylanolytica in 2019. Although this reclassification found support through Genome blast Distance Phylogeny analysis which placed H. xylanolytica X5-1T into the same clade as Hungatella effluvii DSM 24995T, it was contradicted by 16S rRNA gene phylogenetic analysis, which associated it with a set of misnamed Clostridium species later reassigned into the genus Lacrimispora. To ascertain its taxonomic position, comparative analyses were performed to re-examine the relationship between H. xylanolytica X5-1T and all species of the genera Hungatella and Lacrimispora. The ranges of 16S rRNA gene sequence similarity, average amino acid identity, and percentage of conserved protein prediction values were higher between H. xylanolytica X5-1T and species of the genus Lacrimispora than Hungatella. In addition, H. xylanolytica X5-1T was found to harbour genes and pathways conserved and exclusive to species within the genus Lacrimispora but not Hungatella. Essentially, in both the 16S rRNA gene phylogenetic tree and the core-genome phylogenomic tree, H. xylanolytica X5-1T clustered into the same clade as species of the genus Lacrimispora, distinct from species of the genus Hungatella. It is thus clear that H. xylanolytica X5-1T represents a species within the genus Lacrimispora, which we propose to reclassify as Lacrimispora xylanisolvens nom. nov. Finally, based on the results from the phylogenetic and comparative analyses, the genus Hungatella was transferred to the family Lachnospiraceae.

Lactiplantibacillus brownii sp. nov., a novel psychrotolerant species isolated from sauerkraut.

A Gram-stain-positive, rod-shaped, facultatively anaerobic and homofermentative strain, named WILCCON 0030T, was isolated from sauerkraut (fermented cabbage) collected from a local market in the Moscow region of Russia. Comparative analyses based on 16S rRNA gene sequence similarity and whole genome relatedness indicated that strain WILCCON 0030T was most closely related to the type strains Lactiplantibacillus nangangensis NCIMB 15186T, Lactiplantibacillus daoliensis LMG 31171T and Lactiplantibacillus pingfangensis LMG 31176T. However, the average nucleotide identity and digital DNA-DNA hybridization prediction values with these closest relatives only ranged from 84.6 to 84.9 % and from 24.1 to 24.7 %, respectively, and were below the 95.0 and 70.0% thresholds for species delineation. Substantiated by further physiological and biochemical analyses, strain WILCCON 0030T represents a novel species within the genus Lactiplantibacillus for which we propose the name Lactiplantibacillus brownii sp. nov. (type strain WILCCON 0030T=DSM 116485T=LMG 33211T).

Ligilactobacillus ubinensis sp. nov., a novel species isolated from the wild ferment of a durian fruit (Durio zibethinus).

A Gram-stain-positive, rod-shaped, non-spore-forming, catalase-negative, urease-negative, homofermentative and facultatively anaerobic strain, named WILCCON 0076T, was isolated from a wild ferment of pieces of a 'Kampung' durian fruit collected on the island of Ubin (Pulau Ubin), Singapore. The durian had fallen to the ground from a durian tree (Durio zibethinus), on which a group of long-tailed macaques had been observed picking and eating the fruits. Comparative analyses of 16S rRNA gene sequences indicated that WILCCON 0076T potentially represented a novel species within the genus Ligilactobacillus, with the most closely related type strain being Ligilactobacillus agilis DSM 20509T (16S rRNA gene sequence similarity of 97.2 %). Average nucleotide identity and digital DNA-DNA hybridization prediction values were only 86.0% and 18.9 %, respectively. On the basis of the results of a polyphasic approach that included phylogenomic, chemotaxonomic and morphological analyses, we propose a novel species with the name Ligilactobacillus ubinensis sp. nov. (type strain WILCCON 0076T=DSM 114293T=LMG 32698T).

Genome-wide high-throughput signal peptide screening via plasmid pUC256E improves protease secretion in Lactiplantibacillus plantarum and Pediococcus acidilactici.

BACKGROUND: Proteases catalyze the hydrolysis of peptide bonds of proteins, thereby improving dietary protein digestibility, nutrient availability, as well as flavor and texture of fermented food and feed products. The lactobacilli Lactiplantibacillus plantarum (formerly Lactobacillus plantarum) and Pediococcus acidilactici are widely used in food and feed fermentations due to their broad metabolic capabilities and safe use. However, extracellular protease activity in these two species is low. Here, we optimized protease expression and secretion in L. plantarum and P. acidilactici via a genetic engineering strategy. RESULTS: To this end, we first developed a versatile and stable plasmid, pUC256E, which can propagate in both L. plantarum and P. acidilactici. We then confirmed expression and secretion of protease PepG1 as a functional enzyme in both strains with the aid of the previously described L. plantarum-derived signal peptide LP_0373. To further increase secretion of PepG1, we carried out a genome-wide experimental screening of signal peptide functionality. A total of 155 predicted signal peptides originating from L. plantarum and 110 predicted signal peptides from P. acidilactici were expressed and screened for extracellular proteolytic activity in the two different strains, respectively. We identified 12 L. plantarum signal peptides and eight P. acidilactici signal peptides that resulted in improved yield of secreted PepG1. No significant correlation was found between signal peptide sequence properties and its performance with PepG1. CONCLUSION: The vector developed here provides a powerful tool for rapid experimental screening of signal peptides in both L. plantarum and P. acidilactici. Moreover, the set of novel signal peptides identified was widely distributed across strains of the same species and even across some closely related species. This indicates their potential applicability also for the secretion of other proteins of interest in other L. plantarum or P. acidilactici host strains. Our findings demonstrate that screening a library of homologous signal peptides is an attractive strategy to identify the optimal signal peptide for the target protein, resulting in improved protein export.

Characterization and rank assignment criteria for the anaerobic fungi (Neocallimastigomycota).

Establishing a solid taxonomic framework is crucial for enabling discovery and documentation efforts. This ensures effective communication between scientists as well as reproducibility of results between laboratories, and facilitates the exchange and preservation of biological material. Such framework can only be achieved by establishing clear criteria for taxa characterization and rank assignment. Within the anaerobic fungi (phylum Neocallimastigomycota), the need for such criteria is especially vital. Difficulties associated with their isolation, maintenance and long-term storage often result in limited availability and loss of previously described taxa. To this end, we provide here a list of morphological, microscopic, phylogenetic and phenotypic criteria for assessment and documentation when characterizing newly obtained Neocallimastigomycota isolates. We also recommend a polyphasic rank-assignment scheme for novel genus-, species- and strain-level designations for newly obtained Neocallimastigomycota isolates.

Phylotype-Level Characterization of Complex Communities of Lactobacilli Using a High-Throughput, High-Resolution Phenylalanyl-tRNA Synthetase (pheS) Gene Amplicon Sequencing Approach.

The lactobacilli identified to date encompass more than 270 closely related species that were recently reclassified into 26 genera. Because of their relevance to industry, there is a need to distinguish between closely related and yet metabolically and regulatory distinct species, e.g., during monitoring of biotechnological processes or screening of samples of unknown composition. Current available methods, such as shotgun metagenomics or rRNA gene-based amplicon sequencing, have significant limitations (high cost, low resolution, etc.). Here, we generated a phylogeny of lactobacilli based on phenylalanyl-tRNA synthetase (pheS) genes and, from it, developed a high-resolution taxonomic framework which allows for comprehensive and confident characterization of the community diversity and structure of lactobacilli at the species level. This framework is based on a total of 445 pheS gene sequences, including sequences of 276 validly described species and subspecies (of a total of 282, including the proposed L. timonensis species and the reproposed L. zeae species; coverage of 98%), and allows differentiation between 265 species-level clades of lactobacilli and the subspecies of L. sakei The methodology was validated through next-generation sequencing of mock communities. At a sequencing depth of ∼30,000 sequences, the minimum level of detection was approximately 0.02 pg per µl DNA (equaling approximately 10 genome copies per µl template DNA). The pheS approach, along with parallel sequencing of partial 16S rRNA genes, revealed considerable diversity of lactobacilli and distinct community structures across a broad range of samples from different environmental niches. This novel complementary approach may be applicable to industry and academia alike.IMPORTANCE Species formerly classified within the genera Lactobacillus and Pediococcus have been studied extensively at the genomic level. To accommodate their exceptional functional diversity, the over 270 species were recently reclassified into 26 distinct genera. Despite their relevance to both academia and industry, methods that allow detailed exploration of their ecology are still limited by low resolution, high cost, or copy number variations. The approach described here makes use of a single-copy marker gene which outperforms other markers with regard to species-level resolution and availability of reference sequences (98% coverage). The tool was validated against a mock community and used to address diversity of lactobacilli and community structure in various environmental matrices. Such analyses can now be performed at a broader scale to assess and monitor the assembly, structure, and function of communities of lactobacilli at the species level (and, in some cases, even at the subspecies level) across a wide range of academic and commercial applications.

Methane yield phenotypes linked to differential gene expression in the sheep rumen microbiome.

Ruminant livestock represent the single largest anthropogenic source of the potent greenhouse gas methane, which is generated by methanogenic archaea residing in ruminant digestive tracts. While differences between individual animals of the same breed in the amount of methane produced have been observed, the basis for this variation remains to be elucidated. To explore the mechanistic basis of this methane production, we measured methane yields from 22 sheep, which revealed that methane yields are a reproducible, quantitative trait. Deep metagenomic and metatranscriptomic sequencing demonstrated a similar abundance of methanogens and methanogenesis pathway genes in high and low methane emitters. However, transcription of methanogenesis pathway genes was substantially increased in sheep with high methane yields. These results identify a discrete set of rumen methanogens whose methanogenesis pathway transcription profiles correlate with methane yields and provide new targets for CH4 mitigation at the levels of microbiota composition and transcriptional regulation.

An adhesin from hydrogen-utilizing rumen methanogen Methanobrevibacter ruminantium†M1 binds a broad range of hydrogen-producing microorganisms.

Symbiotic associations are ubiquitous in the microbial world and have a major role in shaping the evolution of both partners. One of the most interesting mutualistic relationships exists between protozoa and methanogenic archaea in the fermentative forestomach (rumen) of ruminant animals. Methanogens reside within and on the surface of protozoa as symbionts, and interspecies hydrogen transfer is speculated to be the main driver for physical associations observed between the two groups. In silico analyses of several rumen methanogen genomes have previously shown that up to 5% of genes encode adhesin-like proteins, which may be central to rumen interspecies attachment. We hypothesized that adhesin-like proteins on methanogen cell surfaces facilitate attachment to protozoal hosts. Using phage display technology, we have identified a protein (Mru_1499) from Methanobrevibacter ruminantium M1 as an adhesin that binds to a broad range of rumen protozoa (including the genera Epidinium and Entodinium). This unique adhesin also binds the cell surface of the bacterium Butyrivibrio proteoclasticus, suggesting a broad adhesion spectrum for this protein.

Natural variation in methane emission of sheep fed on a lucerne pellet diet is unrelated to rumen ciliate community type.

Only limited information is available on the roles of different rumen ciliate community types, first described by Eadie in 1962, in enteric methane (CH4) formation by their ruminant hosts. If the different types were differentially associated with CH4 formation, then ciliate community typing could be used to identify naturally high and low CH4-emitting animals. Here we measured the CH4 yields [g CH4 (kg feed dry matter intake, DMI)(-1)] of 118 sheep fed a standard pelleted lucerne diet at two different times, at least 2 weeks apart. There were significant differences (P < 2.2 × 10(-16), Wilcoxon rank sum test) in the CH4 yields (± sd) from sheep selected as high [16.7 ± 1.5 g CH4 (kg DMI)(-1)] and low emitters [13.3 ± 1.5 g CH4 (kg DMI)(-1)]. A rumen sample was collected after each of the two measurements, and ciliate composition was analysed using barcoded 454 Titanium pyrosequencing of 18S rRNA genes. The genera found, in order of mean relative abundance, were Epidinium, Entodinium, Dasytricha, Eudiplodinium, Polyplastron, Isotricha and Anoplodinium-Diplodinium, none of which was significantly correlated with the CH4 emissions ranking associated with the rumen sample. Ciliate communities naturally assembled into four types (A, AB, B and O), characterized by the presence and absence of key genera. There was no difference in CH4 yield between sheep that harboured different ciliate community types, suggesting that these did not underlie the natural variation in CH4 yields. Further research is needed to unravel the nature of interactions between ciliate protozoa and other rumen micro-organisms, which may ultimately lead to contrasting CH4 emission phenotypes.

Few highly abundant operational taxonomic units dominate within rumen methanogenic archaeal species in New Zealand sheep and cattle.

Sequencing and analyses of 16S rRNA gene amplicons were performed to estimate the composition of the rumen methanogen community in 252 samples from eight cohorts of sheep and cattle, separated into 16 different sample groups by diet, and to determine which methanogens are most prominent in the rumens of farmed New Zealand ruminants. Methanobacteriales (relative abundance ± standard deviation, 89.6% ± 9.8%) and Methanomassiliicoccales (10.4% ± 9.8%) were the two major orders and contributed 99.98% (±0.1%) to the rumen methanogen communities in the samples. Sequences from Methanobacteriales were almost entirely from only four different species (or clades of very closely related species). Each was detectable in at least 89% of the samples. These four species or clades were the Methanobrevibacter gottschalkii clade and Methanobrevibacter ruminantium clade with a mean abundance of 42.4% (±19.5% standard deviation) and 32.9% (±18.8%), respectively, and Methanosphaera sp. ISO3-F5 (8.2% ± 6.7%) and Methanosphaera sp. group5 (5.6% ± 5.7%). These four species or clades appeared to be primarily represented by only one or, in one case, two dominant sequence types per species or clade when the sequences were grouped into operational taxonomic units (OTUs) at 99% sequence identity. The mean relative abundance of Methanomassiliicoccales in the samples was relatively low but exceeded 40% in some of the treatment groups. Animal feed affected the apparent methanogen community structure of both orders, as evident from differences in relative abundances of the major OTUs in animals under different feeding regimens.

Buccal swabbing as a noninvasive method to determine bacterial, archaeal, and eukaryotic microbial community structures in the rumen.

Analysis of rumen microbial community structure based on small-subunit rRNA marker genes in metagenomic DNA samples provides important insights into the dominant taxa present in the rumen and allows assessment of community differences between individuals or in response to treatments applied to ruminants. However, natural animal-to-animal variation in rumen microbial community composition can limit the power of a study considerably, especially when only subtle differences are expected between treatment groups. Thus, trials with large numbers of animals may be necessary to overcome this variation. Because ruminants pass large amounts of rumen material to their oral cavities when they chew their cud, oral samples may contain good representations of the rumen microbiota and be useful in lieu of rumen samples to study rumen microbial communities. We compared bacterial, archaeal, and eukaryotic community structures in DNAs extracted from buccal swabs to those in DNAs from samples collected directly from the rumen by use of a stomach tube for sheep on four different diets. After bioinformatic depletion of potential oral taxa from libraries of samples collected via buccal swabs, bacterial communities showed significant clustering by diet (R = 0.37; analysis of similarity [ANOSIM]) rather than by sampling method (R = 0.07). Archaeal, ciliate protozoal, and anaerobic fungal communities also showed significant clustering by diet rather than by sampling method, even without adjustment for potentially orally associated microorganisms. These findings indicate that buccal swabs may in future allow quick and noninvasive sampling for analysis of rumen microbial communities in large numbers of ruminants.

Phylogeny of intestinal ciliates, including Charonina ventriculi, and comparison of microscopy and 18S rRNA gene pyrosequencing for rumen ciliate community structure analysis.

The development of high-throughput methods, such as the construction of 18S rRNA gene clone or pyrosequencing libraries, has allowed evaluation of ciliate community composition in hundreds of samples from the rumen and other intestinal habitats. However, several genera of mammalian intestinal ciliates have been described based only on morphological features and, to date, have not been identified using molecular methods. Here, we isolated single cells of one of the smallest but widely distributed intestinal ciliates, Charonina ventriculi, and sequenced its 18S rRNA gene. We verified the sequence in a full-cycle rRNA approach using fluorescence in situ hybridization and thereby assigned an 18S rRNA gene sequence to this species previously known only by its morphology. Based on its full-length 18S rRNA gene sequence, Charonina ventriculi was positioned within the phylogeny of intestinal ciliates in the subclass Trichostomatia. The taxonomic framework derived from this phylogeny was used for taxonomic assignment of trichostome ciliate 18S rRNA gene sequence data stemming from high-throughput amplicon pyrosequencing of rumen-derived DNA samples. The 18S rRNA gene-based ciliate community structure was compared to that obtained from microscopic counts using the same samples. Both methods allowed identification of dominant members of the ciliate communities and classification of the rumen ciliate community into one of the types first described by Eadie in 1962. Notably, each method is associated with advantages and disadvantages. Microscopy is a highly accurate method for evaluation of total numbers or relative abundances of different ciliate genera in a sample, while 18S rRNA gene pyrosequencing represents a valuable alternative for comparison of ciliate community structure in a large number of samples from different animals or treatment groups.

Whole-Genome Sequence of Ligilactobacillus faecis WILCCON 0062, Isolated from Feces of a Wild Boar (Sus scrofa).

We report the whole genome of a strain of Ligilactobacillus faecis. The complete circular chromosome and plasmid of strain WILCCON 0062 were obtained through a combination of short- and long-read sequencing and may be used to derive unprecedented insights into the genome-level phylogeny and functional capacities of Ligilactobacillus faecis.

Genomic insights into the physiology of Quinella, an iconic uncultured rumen bacterium.

Quinella is a genus of iconic rumen bacteria first reported in 1913. There are no cultures of these bacteria, and information on their physiology is scarce and contradictory. Increased abundance of Quinella was previously found in the rumens of some sheep that emit low amounts of methane (CH4) relative to their feed intake, but whether Quinella contributes to low CH4 emissions is not known. Here, we concentrate Quinella cells from sheep rumen contents, extract and sequence DNA, and reconstruct Quinella genomes that are >90% complete with as little as 0.20% contamination. Bioinformatic analyses of the encoded proteins indicate that lactate and propionate formation are major fermentation pathways. The presence of a gene encoding a potential uptake hydrogenase suggests that Quinella might be able to use free hydrogen (H2). None of the inferred metabolic pathways is predicted to produce H2, a major precursor of CH4, which is consistent with the lower CH4 emissions from those sheep with high abundances of this bacterium.

Five-week dietary exposure to dry diets alters the faecal bacterial populations in the domestic cat (Felis catus).

The effects of wet (canned) or dry (kibbled) diets on faecal bacterial populations in the cat were investigated in eight domestic short-haired cats (four males and four females; averaging 6 years of age and 3.4 kg) in a nested design. The cats were fed ad libitum a commercially available wet diet (moisture 82.0 %, crude protein 51.7 %, fat 28.9 %, carbohydrate (CHO) 8.9 % and ash 10.6 % DM) for 5 weeks. On the fifth week, individual feed intakes and faecal outputs were determined. Fresh faecal samples were collected twice daily, mixed for homogeneity, subsampled and stored at - 85 °C until analysis. The cats were then switched to a commercially available dry diet (moisture 8.5 %, crude protein 33.0 %, fat 11.0 %, CHO 49.4 % and ash 6.6 % DM) for 5 weeks, and fresh faeces were sampled as described previously. Energy intake tended to be higher in cats fed dry diets (P < 0.10), but body weight was similar between the two feeding periods (P>0.05). Denaturing gradient gel electrophoresis (DGGE) of bacterial 16S rRNA genes amplified from DNA extracted from faeces was performed. The unweighted pair group method with arithmetic mean cluster analysis of bacterial community profiles using Pearson's correlation revealed diet-specific clustering when the same cats were fed on either a dry or a wet diet (dissimilarity between the groups, 88.6 %; P < 0.001). Subsequent cloning and sequencing of five selected distinct DGGE bands indicated that members of the Pelomonas and Fusobacteriaceae were influenced by a short-term change in diet format. This suggests that 5-week dietary exposure is sufficient to alter gastrointestinal microflora.

Effects of long-acting, broad spectra anthelmintic treatments on the rumen microbial community compositions of grazing sheep.

Anthelmintic treatment of adult ewes is widely practiced to remove parasite burdens in the expectation of increased ruminant productivity. However, the broad activity spectra of many anthelmintic compounds raises the possibility of impacts on the rumen microbiota. To investigate this, 300 grazing ewes were allocated to treatment groups that included a 100-day controlled release capsule (CRC) containing albendazole and abamectin, a long-acting moxidectin injection (LAI), and a non-treated control group (CON). Rumen bacterial, archaeal and protozoal communities at day 0 were analysed to identify 36 sheep per treatment with similar starting compositions. Microbiota profiles, including those for the rumen fungi, were then generated for the selected sheep at days 0, 35 and 77. The CRC treatment significantly impacted the archaeal community, and was associated with increased relative abundances of Methanobrevibacter ruminantium, Methanosphaera sp. ISO3-F5, and Methanomassiliicoccaceae Group 12 sp. ISO4-H5 compared to the control group. In contrast, the LAI treatment increased the relative abundances of members of the Veillonellaceae and resulted in minor changes to the bacterial and fungal communities by day 77. Overall, the anthelmintic treatments resulted in few, but highly significant, changes to the rumen microbiota composition.

Resilience of Faecal Microbiota in Stabled Thoroughbred Horses Following Abrupt Dietary Transition between Freshly Cut Pasture and Three Forage-Based Diets.

The management of competition horses in New Zealand often involves rotations of short periods of stall confinement and concentrate feeding, with periods of time at pasture. Under these systems, horses may undergo abrupt dietary changes, with the incorporation of grains or concentrate feeds to the diet to meet performance needs, or sudden changes in the type of forage fed in response to a lack of fresh or conserved forage. Abrupt changes in dietary management are a risk factor for gastrointestinal (GI) disturbances, potentially due to the negative effects observed on the population of GI microbiota. In the present study, the faecal microbiota of horses was investigated to determine how quickly the bacterial communities; (1) responded to dietary change, and (2) stabilised following abrupt dietary transition. Six Thoroughbred mares were stabled for six weeks, consuming freshly cut pasture (weeks 1, 3 and 5), before being abruptly transitioned to conserved forage-based diets, both offered ad libitum. Intestinal markers were administered to measure digesta transit time immediately before each diet change. The conserved forage-based diets were fed according to a 3 × 3 Latin square design (weeks 2, 4 and 6), and comprised a chopped ensiled forage fed exclusively (Diet FE) or with whole oats (Diet FE + O), and perennial ryegrass hay fed with whole oats (Diet H + O). Faecal samples were collected at regular intervals from each horse following the diet changes. High throughput 16S rRNA gene sequencing was used to evaluate the faecal microbiota. There were significant differences in alpha diversity across diets (p < 0.001), and a significant effect of diet on the beta diversity (ANOSIM, p = 0.001), with clustering of samples observed by diet group. There were differences in the bacterial phyla across diets (p < 0.003), with the highest relative abundances observed for Firmicutes (62-64%) in the two diets containing chopped ensiled forage, Bacteroidetes (32-38%) in the pasture diets, and Spirochaetes (17%) in the diet containing hay. Major changes in relative abundances of faecal bacteria appeared to correspond with the cumulative percentage of intestinal markers retrieved in the faeces as the increasing amounts of digesta from each new diet transited the animals. A stable faecal microbiota profile was observed in the samples from 96 h after abrupt transition to the treatment diets containing ensiled chopped forage. The present study confirmed that the diversity and community structure of the faecal bacteria in horses is diet-specific and resilient following dietary transition and emphasised the need to have modern horse feeding management that reflects the ecological niche, particularly by incorporating large proportions of forage into equine diets.

Seasonal Variation in the Faecal Microbiota of Mature Adult Horses Maintained on Pasture in New Zealand.

Seasonal variation in the faecal microbiota of forage-fed horses was investigated over a 12-month period to determine whether the bacterial diversity fluctuated over time. Horses (n = 10) were maintained on pasture for one year, with hay supplemented from June to October. At monthly intervals, data were recorded on pasture availability and climate (collected continuously and averaged on monthly basis), pasture and hay samples were collected for nutrient analysis, and faecal samples were collected from all horses to investigate the diversity of faecal microbiota using next-generation sequencing on the Illumina MiSeq platform. The alpha diversity of bacterial genera was high in all samples (n = 118), with significantly higher Simpson's (p < 0.001) and Shannon-Wiener (p < 0.001) diversity indices observed during the months when horses were kept exclusively on pasture compared to the months when pasture was supplemented with hay. There were significant effects of diet, season, and month (ANOSIM, p < 0.01 for each comparison) on the beta diversity of bacterial genera identified in the faeces. While there was some inter-horse variation, hierarchical clustering of beta diversity indices showed separate clades originating for samples obtained during May, June, and July (late-autumn to winter period), and January, February, and March (a period of drought), with a strong association between bacterial taxa and specific nutrients (dry matter, protein, and structural carbohydrates) and climate variables (rainfall and temperature). Our study supports the hypothesis that the diversity and community structure of the faecal microbiota of horses kept on pasture varied over a 12-month period, and this variation reflects changes in the nutrient composition of the pasture, which in turn is influenced by climatic conditions. The findings of this study may have implications for grazing management and the preparation of conserved forages for those horses susceptible to perturbations of the hindgut microbiota.

A restriction enzyme reduced representation sequencing approach for low-cost, high-throughput metagenome profiling.

Microbial community profiles have been associated with a variety of traits, including methane emissions in livestock. These profiles can be difficult and expensive to obtain for thousands of samples (e.g. for accurate association of microbial profiles with traits), therefore the objective of this work was to develop a low-cost, high-throughput approach to capture the diversity of the rumen microbiome. Restriction enzyme reduced representation sequencing (RE-RRS) using ApeKI or PstI, and two bioinformatic pipelines (reference-based and reference-free) were compared to bacterial 16S rRNA gene sequencing using repeated samples collected two weeks apart from 118 sheep that were phenotypically extreme (60 high and 58 low) for methane emitted per kg dry matter intake (n = 236). DNA was extracted from freeze-dried rumen samples using a phenol chloroform and bead-beating protocol prior to RE-RRS. The resulting sequences were used to investigate the repeatability of the rumen microbial community profiles, the effect of laboratory and analytical method, and the relationship with methane production. The results suggested that the best method was PstI RE-RRS analyzed with the reference-free approach, which accounted for 53.3±5.9% of reads, and had repeatabilities of 0.49±0.07 and 0.50±0.07 for the first two principal components (PC1 and PC2), phenotypic correlations with methane yield of 0.43±0.06 and 0.46±0.06 for PC1 and PC2, and explained 41±8% of the variation in methane yield. These results were significantly better than for bacterial 16S rRNA gene sequencing of the same samples (p<0.05) except for the correlation between PC2 and methane yield. A Sensitivity study suggested approximately 2000 samples could be sequenced in a single lane on an Illumina HiSeq 2500, meaning the current work using 118 samples/lane and future proposed 384 samples/lane are well within that threshold. With minor adaptations, our approach could be used to obtain microbial profiles from other metagenomic samples.