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 thermophilic Campylobacter jejuni and Campylobacter coli are considered weakly clonal populations where incongruences between genetic markers are assumed to be due to random horizontal transfer of genomic DNA. In order to investigate the population genetics structure we extracted a set of 1180 core gene families (CGF) from 27 sequenced genomes of C. jejuni and C. coli. We adopted a principal component analysis (PCA) on the normalized evolutionary distances in order to reveal any patterns in the evolutionary signals contained within the various CGFs.

Antibody-mediated immunity is highly protective against disease. The majority of current vaccines confer protection through humoral immunity, but there is high variability in responsiveness across populations. Identifying immune mechanisms that mediate low antibody responsiveness may provide potential strategies to boost vaccine efficacy. Here, we report diverse antibody responsiveness to unadjuvanted as well as adjuvanted immunization in substrains of BALB/c mice, resulting in high and low antibody response phenotypes. Furthermore, these antibody phenotypes were not affected by changes in environmental factors such as the gut microbiota composition. Antigen-specific B cells following immunization had a marked difference in capability to class switch, resulting in perturbed IgG isotype antibody production. In vitro, a B-cell intrinsic defect in the regulation of class-switch recombination was identified in mice with low IgG antibody production. Whole genome sequencing identified polymorphisms associated with the magnitude of antibody produced, and we propose candidate genes that may regulate isotype class-switching capability. This study highlights that mice sourced from different vendors can have significantly altered humoral immune response profiles, and provides a resource to interrogate genetic regulators of antibody responsiveness. Together these results further our understanding of immune heterogeneity and suggest additional research on the genetic influences of adjuvanted vaccine strategies is warranted for enhancing vaccine efficacy.

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.

Hydrogen cross-feeding microbes form a functionally important subset of the human colonic microbiota. The three major hydrogenotrophic functional groups of the colon: sulphate-reducing bacteria (SRB), methanogens and reductive acetogens, have been linked to wide ranging impacts on host physiology, health and wellbeing.

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.

The bacterial production of acetate via reductive acetogenesis along the Wood-Ljungdahl metabolic pathway is an important source of this molecule in several environments, ranging from industrial bioreactors to the human gastrointestinal tract. Here, we contributed to the study of reductive acetogens by considering mathematical modelling techniques for the prediction of bacterial growth and acetate production. We found that the incorporation of a hydrogen uptake concentration threshold into the models improves their predictions and we calculated this threshold as 86.2 mM (95% confidence interval 6.1-132.6 mM). Monod kinetics and first-order kinetics models, with the inclusion of two candidate threshold terms or reversible Michaelis-Menten kinetics, were compared to experimental data and the optimal formulation for predicting both growth and metabolism was found. The models were then used to compare the efficacy of two growth media for acetogens. We found that the recently described general acetogen medium was superior to the DSMZ medium in terms of unbiased estimation of acetogen growth and investigated the contribution of yeast extract concentration to acetate production and bacterial growth in culture. The models and their predictions will be useful to those studying both industrially and environmentally relevant reductive acetogenesis and allow for straightforward adaptation to similar cases with different organisms.

Nonribosomal peptide synthetases (NRPSs) are multimodular enzymes that produce a wide range of bioactive peptides, such as siderophores, toxins, and antibacterial and insecticidal agents. NRPSs are dynamic proteins characterized by extensive interdomain communications as a consequence of their assembly-line mode of synthesis. Hence, crystal structures of multidomain fragments of NRPSs have aided in elucidating crucial interdomain interactions that occur during different steps of the NRPS catalytic cycle. One crucial yet unexplored interaction is that between the reductase (R) domain and the peptide carrier protein (PCP) domain. R domains are members of the short-chain dehydrogenase/reductase family and function as termination domains that catalyze the reductive release of the final peptide product from the terminal PCP domain of the NRPS. Here, we report the crystal structure of an archaeal NRPS PCP-R didomain construct. This is the first NRPS R domain structure to be determined together with the upstream PCP domain and is also the first structure of an archaeal NRPS to be reported. The structure reveals that a novel helix-turn-helix motif, found in NRPS R domains but not in other short-chain dehydrogenase/reductase family members, plays a major role in the interface between the PCP and R domains. The information derived from the described PCP-R interface will aid in gaining further mechanistic insights into the peptide termination reaction catalyzed by the R domain and may have implications in engineering NRPSs to synthesize novel peptide products.

Eight Faecalibacterium sp. strains were isolated from feces of healthy human volunteers. Here, we describe their genome sequences. The genome sizes ranged from 2.78 Mbp to 3.23 Mbp, with an average GC content of 56.6% and encoding 2,795 protein-coding genes on average.

Lactobacillus fermentum is commonly found in food products, and some strains are known to have beneficial effects on human health. However, our previous research indicated that L. fermentum AGR1487 decreases in vitro intestinal barrier integrity. The hypothesis was that cell surface structures of AGR1487 are responsible for the observed in vitro effect. AGR1487 was compared to another human oral L. fermentum strain, AGR1485, which does not cause the same effect. The examination of phenotypic traits associated with the composition of cell surface structures showed that compared to AGR1485, AGR1487 had a smaller genome, utilized different sugars, and had greater tolerance to acid and bile. The effect of the two strains on intestinal barrier integrity was determined using two independent measures of paracellular permeability of the intestinal epithelial Caco-2 cell line. The transepithelial electrical resistance (TEER) assay specifically measures ion permeability, whereas the mannitol flux assay measures the passage of uncharged molecules. Both live and UV-inactivated AGR1487 decreased TEER across Caco-2 cells implicating the cell surfaces structures in the effect. However, only live AGR1487, and not UV-inactivated AGR1487, increased the rate of passage of mannitol, implying that a secreted component(s) is responsible for this effect. These differences in barrier integrity results are likely due to the TEER and mannitol flux assays measuring different characteristics of the epithelial barrier, and therefore imply that there are multiple mechanisms involved in the effect of AGR1487 on barrier integrity.

Lateral gene transfer (LGT) appears to promote genotypic and phenotypic variation in microbial communities in a range of environments, including the mammalian intestine. However, the extent and mechanisms of LGT in intestinal microbial communities of non-mammalian hosts remains poorly understood.

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.

Equally important and challenging as genome annotation, is the subsequent classification of predicted genes into their respective pathways. The Kyoto Encyclopedia of Genes and Genomes (KEGG) represents a database consisting of known genes and their respective biochemical functionalities. Although accessible online, analyses of multiple genes are time consuming and are not suitable for analyzing data sets that are proprietary.

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.

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.

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.