The goal of the Hungate1000 project is to generate a reference set of rumen microbial genome sequences. Toward this goal we have carried out a meta-analysis using information from culture collections, scientific literature, and the NCBI and RDP databases and linked this with a comparative study of several rumen 16S rRNA gene-based surveys. In this way we have attempted to capture a snapshot of rumen bacterial diversity to examine the culturable fraction of the rumen bacterial microbiome. Our analyses have revealed that for cultured rumen bacteria, there are many genera without a reference genome sequence. Our examination of culture-independent studies highlights that there are few novel but many uncultured taxa within the rumen bacterial microbiome. Taken together these results have allowed us to compile a list of cultured rumen isolates that are representative of abundant, novel and core bacterial species in the rumen. In addition, we have identified taxa, particularly within the phylum Bacteroidetes, where further cultivation efforts are clearly required. This information is being used to guide the isolation efforts and selection of bacteria from the rumen microbiota for sequencing through the Hungate1000.

The rumen is known to harbor dense populations of bacteriophages (phages) predicted to be capable of infecting a diverse range of rumen bacteria. While bacterial genome sequencing projects are revealing the presence of phages which can integrate their DNA into the genome of their host to form stable, lysogenic associations, little is known of the genetics of phages which utilize lytic replication. These phages infect and replicate within the host, culminating in host lysis, and the release of progeny phage particles. While lytic phages for rumen bacteria have been previously isolated, their genomes have remained largely uncharacterized. Here we report the first complete genome sequences of lytic phage isolates specifically infecting three genera of rumen bacteria: Bacteroides, Ruminococcus, and Streptococcus. All phages were classified within the viral order Caudovirales and include two phage morphotypes, representative of the Siphoviridae and Podoviridae families. The phage genomes displayed modular organization and conserved viral genes were identified which enabled further classification and determination of closest phage relatives. Co-examination of bacterial host genomes led to the identification of several genes responsible for modulating phage:host interactions, including CRISPR/Cas elements and restriction-modification phage defense systems. These findings provide new genetic information and insights into how lytic phages may interact with bacteria of the rumen microbiome.

Bacterial species belonging to the genus Pseudobutyrivibrio are important members of the rumen microbiome contributing to the degradation of complex plant polysaccharides. Pseudobutyrivibrio xylanivorans MA3014 was selected for genome sequencing to examine its ability to breakdown and utilize plant polysaccharides. The complete genome sequence of MA3014 is 3.58 Mb, consists of three replicons (a chromosome, chromid, and plasmid), has an overall G + C content of 39.6%, and encodes 3,265 putative protein-coding genes (CDS). Comparative pan-genomic analysis of all cultivated and currently available P. xylanivorans genomes has revealed a strong correlation of orthologous genes within this rumen bacterial species. MA3014 is metabolically versatile and capable of growing on a range of simple mono- or oligosaccharides derived from complex plant polysaccharides such as pectins, mannans, starch, and hemicelluloses, with lactate, butyrate, and formate as the principal fermentation end products. The genes encoding these metabolic pathways have been identified and MA3014 is predicted to encode an extensive range of Carbohydrate-Active enZYmes with 78 glycoside hydrolases, 13 carbohydrate esterases, and 54 glycosyl transferases, suggesting an important role in solubilization of plant matter in the rumen.

Plant polysaccharide breakdown by microbes in the rumen is fundamental to digestion in ruminant livestock. Bacterial species belonging to the rumen genera Butyrivibrio and Pseudobutyrivibrio are important degraders and utilizers of lignocellulosic plant material. These bacteria degrade polysaccharides and ferment the released monosaccharides to yield short-chain fatty acids that are used by the ruminant for growth and the production of meat, milk, and fiber products. Although rumen Butyrivibrio and Pseudobutyrivibrio species are regarded as common rumen inhabitants, their polysaccharide-degrading and carbohydrate-utilizing enzymes are not well understood. In this study, we analyzed the genomes of 40 Butyrivibrio and 6 Pseudobutyrivibrio strains isolated from the plant-adherent fraction of New Zealand dairy cows to explore the polysaccharide-degrading potential of these important rumen bacteria. Comparative genome analyses combined with phylogenetic analysis of their 16S rRNA genes and short-chain fatty acid production patterns provide insight into the genomic diversity and physiology of these bacteria and divide Butyrivibrio into 3 species clusters. Rumen Butyrivibrio bacteria were found to encode a large and diverse spectrum of degradative carbohydrate-active enzymes (CAZymes) and binding proteins. In total, 4,421 glycoside hydrolases (GHs), 1,283 carbohydrate esterases (CEs), 110 polysaccharide lyases (PLs), 3,605 glycosyltransferases (GTs), and 1,706 carbohydrate-binding protein modules (CBM) with predicted activities involved in the depolymerization and transport of the insoluble plant polysaccharides were identified. Butyrivibrio genomes had similar patterns of CAZyme families but varied greatly in the number of genes within each category in the Carbohydrate-Active Enzymes database (CAZy), suggesting some level of functional redundancy. These results suggest that rumen Butyrivibrio species occupy similar niches but apply different degradation strategies to be able to coexist in the rumen.IMPORTANCE Feeding a global population of 8 billion people and climate change are the primary challenges facing agriculture today. Ruminant livestock are important food-producing animals, and maximizing their productivity requires an understanding of their digestive systems and the roles played by rumen microbes in plant polysaccharide degradation. Members of the genera Butyrivibrio and Pseudobutyrivibrio are a phylogenetically diverse group of bacteria and are commonly found in the rumen, where they are a substantial source of polysaccharide-degrading enzymes for the depolymerization of lignocellulosic material. Our findings have highlighted the immense enzymatic machinery of Butyrivibrio and Pseudobutyrivibrio species for the degradation of plant fiber, suggesting that these bacteria occupy similar niches but apply different degradation strategies in order to coexist in the competitive rumen environment.

Farmed ruminants are the largest source of anthropogenic methane emissions globally. The methanogenic archaea responsible for these emissions use molecular hydrogen (H2), produced during bacterial and eukaryotic carbohydrate fermentation, as their primary energy source. In this work, we used comparative genomic, metatranscriptomic and co-culture-based approaches to gain a system-wide understanding of the organisms and pathways responsible for ruminal H2 metabolism. Two-thirds of sequenced rumen bacterial and archaeal genomes encode enzymes that catalyse H2 production or consumption, including 26 distinct hydrogenase subgroups. Metatranscriptomic analysis confirmed that these hydrogenases are differentially expressed in sheep rumen. Electron-bifurcating [FeFe]-hydrogenases from carbohydrate-fermenting Clostridia (e.g., Ruminococcus) accounted for half of all hydrogenase transcripts. Various H2 uptake pathways were also expressed, including methanogenesis (Methanobrevibacter), fumarate and nitrite reduction (Selenomonas), and acetogenesis (Blautia). Whereas methanogenesis-related transcripts predominated in high methane yield sheep, alternative uptake pathways were significantly upregulated in low methane yield sheep. Complementing these findings, we observed significant differential expression and activity of the hydrogenases of the hydrogenogenic cellulose fermenter Ruminococcus albus and the hydrogenotrophic fumarate reducer Wolinella succinogenes in co-culture compared with pure culture. We conclude that H2 metabolism is a more complex and widespread trait among rumen microorganisms than previously recognised. There is evidence that alternative hydrogenotrophs, including acetogenic and respiratory bacteria, can prosper in the rumen and effectively compete with methanogens for H2. These findings may help to inform ongoing strategies to mitigate methane emissions by increasing flux through alternative H2 uptake pathways, including through animal selection, dietary supplementation and methanogenesis inhibitors.

Lactococcus lactis subsp. cremoris HP(T) has been widely used in studies of the metabolism of lactococcal dairy starter cultures. A comparison of the draft HP(T) genome with those from other strains of L. lactis subsp. cremoris will aid our understanding of the domestication and evolution of these important industrial cultures.

Ruminant animals contribute significantly to the global value of agriculture and rely on a complex microbial community for efficient digestion. However, little is known of how this microbial-host relationship develops and is maintained. To begin to address this, we have determined the ability of three Bifidobacterium species isolated from the faeces of newborn calves to grow on carbohydrates typical of a newborn ruminant diet. Genome sequences have been determined for these bacteria with analysis of the genomes providing insights into the host association and identification of several genes that may mediate interactions with the ruminant gastrointestinal tract. The present study provides a starting point from which we can define the role of potential beneficial microbes in the nutrition of young ruminants and begin to influence the interactions between the microbiota and the host. The differences observed in genomic content hint at niche partitioning among the bifidobacterial species analysed and the different strategies they employ to successfully adapt to this habitat.

Digestive processes in the rumen lead to the release of methyl-compounds, mainly methanol and methylamines, which are used by methyltrophic methanogens to form methane, an important agricultural greenhouse gas. Methylamines are produced from plant phosphatidylcholine degradation, by choline trimethylamine lyase, while methanol comes from demethoxylation of dietary pectins via pectin methylesterase activity. We have screened rumen metagenomic and metatranscriptomic datasets, metagenome assembled genomes, and the Hungate1000 genomes to identify organisms capable of producing methyl-compounds. We also describe the enrichment of pectin-degrading and methane-forming microbes from sheep rumen contents and the analysis of their genomes via metagenomic assembly.

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.

Methane (CH(4)) is a potent greenhouse gas (GHG), having a global warming potential 21 times that of carbon dioxide (CO(2)). Methane emissions from agriculture represent around 40% of the emissions produced by human-related activities, the single largest source being enteric fermentation, mainly in ruminant livestock. Technologies to reduce these emissions are lacking. Ruminant methane is formed by the action of methanogenic archaea typified by Methanobrevibacter ruminantium, which is present in ruminants fed a wide variety of diets worldwide. To gain more insight into the lifestyle of a rumen methanogen, and to identify genes and proteins that can be targeted to reduce methane production, we have sequenced the 2.93 Mb genome of M. ruminantium M1, the first rumen methanogen genome to be completed.

A large collection of Lactococcus lactis strains, including wild-type isolates and dairy starter cultures, were screened on the basis of their phenotype and the macrorestriction patterns produced from pulsed-field gel electrophoresis (PFGE) analysis of SmaI digests of genomic DNA. Three groups of dairy starter cultures, used for different purposes in the dairy industry, and a fourth group made up of strains isolated from the environment were selected for analysis of their chromosomal diversity using the endonuclease I-CeuI. Chromosome architecture was largely conserved with each strain having six copies of the rRNA genes, and the chromosome size of individual strains ranged between 2,240 and 2,688 kb. The origin of L. lactis strains showed the greatest correlation with chromosome size, and dairy strains, particularly those with the cremoris phenotype, had smaller chromosomes than wild-type strains. Overall, this study, coupled with analysis of the sequenced L. lactis genomes, provides evidence that defined strain dairy starter cultures have arisen from plant L. lactis strains. Adaptation of these strains to the dairy environment has involved loss of functions resulting in smaller chromosomes and acquisition of genes (usually plasmid associated) that facilitate growth in milk. We conclude that dairy starter cultures generally and the industrially used cremoris and diacetylactis phenotype strains in particular comprise a specialized group of L. lactis strains that have been selected to become an essential component of industrial processes and have evolved accordingly, so that they are no longer fit to survive outside the dairy environment.

Productivity of ruminant livestock depends on the rumen microbiota, which ferment indigestible plant polysaccharides into nutrients used for growth. Understanding the functions carried out by the rumen microbiota is important for reducing greenhouse gas production by ruminants and for developing biofuels from lignocellulose. We present 410 cultured bacteria and archaea, together with their reference genomes, representing every cultivated rumen-associated archaeal and bacterial family. We evaluate polysaccharide degradation, short-chain fatty acid production and methanogenesis pathways, and assign specific taxa to functions. A total of 336 organisms were present in available rumen metagenomic data sets, and 134 were present in human gut microbiome data sets. Comparison with the human microbiome revealed rumen-specific enrichment for genes encoding de novo synthesis of vitamin B12, ongoing evolution by gene loss and potential vertical inheritance of the rumen microbiome based on underrepresentation of markers of environmental stress. We estimate that our Hungate genome resource represents ∼75% of the genus-level bacterial and archaeal taxa present in the rumen.

The taxonomy and associated nomenclature of many taxa of rumen bacteria are poorly defined within databases of 16S rRNA genes. This lack of resolution results in inadequate definition of microbial community structures, with large parts of the community designated as incertae sedis, unclassified, or uncultured within families, orders, or even classes. We have begun resolving these poorly-defined groups of rumen bacteria, based on our desire to name these for use in microbial community profiling. We used the previously-reported global rumen census (GRC) dataset consisting of >4.5 million partial bacterial 16S rRNA gene sequences amplified from 684 rumen samples and representing a wide range of animal hosts and diets. Representative sequences from the 8,985 largest operational units (groups of sequence sharing >97% sequence similarity, and covering 97.8% of all sequences in the GRC dataset) were used to identify 241 pre-defined clusters (mainly at genus or family level) of abundant rumen bacteria in the ARB SILVA 119 framework. A total of 99 of these clusters (containing 63.8% of all GRC sequences) had no unique or had inadequate taxonomic identifiers, and each was given a unique nomenclature. We assessed this improved framework by comparing taxonomic assignments of bacterial 16S rRNA gene sequence data in the GRC dataset with those made using the original SILVA 119 framework, and three other frameworks. The two SILVA frameworks performed best at assigning sequences to genus-level taxa. The SILVA 119 framework allowed 55.4% of the sequence data to be assigned to 751 uniquely identifiable genus-level groups. The improved framework increased this to 87.1% of all sequences being assigned to one of 871 uniquely identifiable genus-level groups. The new designations were included in the SILVA 123 release (https://www.arb-silva.de/documentation/release-123/) and will be perpetuated in future releases.

Ruminants are important contributors to global methane emissions via microbial fermentation in their reticulo-rumens. This study is part of a larger program, characterising the rumen microbiomes of sheep which vary naturally in methane yield (g CH4/kg DM/day) and aims to define differences in microbial communities, and in gene and transcript abundances that can explain the animal methane phenotype.

In silico, secretome proteins can be predicted from completely sequenced genomes using various available algorithms that identify membrane-targeting sequences. For metasecretome (collection of surface, secreted and transmembrane proteins from environmental microbial communities) this approach is impractical, considering that the metasecretome open reading frames (ORFs) comprise only 10% to 30% of total metagenome, and are poorly represented in the dataset due to overall low coverage of metagenomic gene pool, even in large-scale projects.

The research was aimed at finding which membrane proteins of the rumen bacterium Butyrivibrio proteoclasticus are involved in the uptake of carbohydrates resulting from extracellular enzymatic degradation of hemicellulose and fructan. The proteomic analysis of cells grown with fructose or xylan as the sole substrate identified 13 membrane proteins predicted to function as carbohydrate transporters. One protein detected was the membrane component of a fructose-specific phosphoenolpyruvate:sugar phosphotransferase system believed to be involved in the fructose uptake following extracellular fructan breakdown. The other 12 proteins were all ABC transport system substrate-binding proteins, nine of which belong to functional category COG1653 that includes proteins predicted to transport oligosaccharides. Four of the SBPs were significantly upregulated in xylan grown cells, and three of these were found in polysaccharide utilisation loci where they are clustered with other genes involved in hemicellulose breakdown and metabolism. It is possible that the carbon source available regulates a wider network of genes. The information on the mechanisms used by rumen bacteria to take up carbohydrates from their environment may improve our understanding of the ruminant digestion and facilitate strategies for improved pasture and stored feed utilisation.

Determining the role of rumen microbes and their enzymes in plant polysaccharide breakdown is fundamental to understanding digestion and maximising productivity in ruminant animals. Butyrivibrio proteoclasticus B316(T) is a gram-positive, butyrate-forming rumen bacterium with a key role in plant polysaccharide degradation. The 4.4 Mb genome consists of 4 replicons; a chromosome, a chromid and two megaplasmids. The chromid is the smallest reported for all bacteria, and the first identified from the phylum Firmicutes. B316 devotes a large proportion of its genome to the breakdown and reassembly of complex polysaccharides and has a highly developed glycobiome when compared to other sequenced bacteria. The secretion of a range of polysaccharide-degrading enzymes which initiate the breakdown of pectin, starch and xylan, a subtilisin family protease active against plant proteins, and diverse intracellular enzymes to break down oligosaccharides constitute the degradative capability of this organism. A prominent feature of the genome is the presence of multiple gene clusters predicted to be involved in polysaccharide biosynthesis. Metabolic reconstruction reveals the absence of an identifiable gene for enolase, a conserved enzyme of the glycolytic pathway. To our knowledge this is the first report of an organism lacking an enolase. Our analysis of the B316 genome shows how one organism can contribute to the multi-organism complex that rapidly breaks down plant material in the rumen. It can be concluded that B316, and similar organisms with broad polysaccharide-degrading capability, are well suited to being early colonizers and degraders of plant polysaccharides in the rumen environment.

Many bacteria are facultative anaerobes, and can proliferate in both anoxic and oxic environments. Under anaerobic conditions, fermentation is the primary means of energy generation in contrast to respiration. Furthermore, the rates and spectra of spontaneous mutations that arise during anaerobic growth differ to those under aerobic growth. A long-term selection experiment was undertaken to investigate the genetic changes that underpin how the facultative anaerobe, Escherichia coli, adapts to anaerobic environments.

Oxidative stress is a major cause of mutation but little is known about how growth in the absence of oxygen impacts the rate and spectrum of mutations. We employed long-term mutation accumulation experiments to directly measure the rates and spectra of spontaneous mutation events in Escherichia coli populations propagated under aerobic and anaerobic conditions. To detect mutations, whole genome sequencing was coupled with methods of analysis sufficient to identify a broad range of mutational classes, including structural variants (SVs) generated by movement of repetitive elements. The anaerobically grown populations displayed a mutation rate nearly twice that of the aerobic populations, showed distinct asymmetric mutational strand biases, and greater insertion element activity. Consistent with mutation rate and spectra observations, genes for transposition and recombination repair associated with SVs were up-regulated during anaerobic growth. Together, these results define differences in mutational spectra affecting the evolution of facultative anaerobes.

The growth and productivity of ruminants depends on a complex microbial community found in their fore-stomach (rumen), which is able to breakdown plant polysaccharides and ferment the released sugars. Butyrivibrio proteoclasticus B316T is a Gram-positive polysaccharide-degrading, butyrate-producing bacterium that is present at high numbers in the rumen of animals consuming pasture or grass silage based diets. B316T is one of a small number of rumen fibrolytic microbes capable of efficiently degrading and utilizing xylan, as well as being capable of utilizing arabinose, xylose, pectin and starch. We have therefore carried out a proteomic analysis of B316T to identify intracellular enzymes that are implicated in the metabolism of internalized xylan. Three hundred and ninety four proteins were identified including enzymes that have potential to metabolize assimilated products of extracellular xylan digestion. Identified enzymes included arabinosidases, esterases, an endoxylanase, and β-xylosidase. The presence of intracellular debranching enzymes indicated that some hemicellulosic side-chains may not be removed until oligosaccharides liberated by extracellular digestion have been assimilated by the cells. The results support a model of extracellular digestion of hemicellulose to oligosaccharides that are then transported to the cytoplasm for further digestion by intracellular enzymes.