Complete genome sequence of Methanosphaera sp. ISO3-F5, a rumen methylotrophic methanogen

Methanosphaera spp. are methylotrophic methanogenic archaea and members of the order Methanobacteriales with few cultured representatives. Methanosphaera sp. ISO3-F5 was isolated from sheep rumen contents in New Zealand. Here, we report its complete genome, consisting of a large chromosome and a megaplasmid (GenBank accession numbers CP118753 and CP118754, respectively).

Methanobrevibacter boviskoreani JH1T growth on alcohols allows development of a high throughput bioassay to detect methanogen inhibition

Rumen methanogenic archaea use by-products of fermentation to carry out methanogenesis for energy generation. A key fermentation by-product is hydrogen (H₂), which acts as the source of reducing potential for methane (CH₄) formation in hydrogenotrophic methanogens. The in vitro cultivation of hydrogenotrophic rumen methanogens requires pressurised H₂ which limits the ability to conduct high-throughput screening experiments with these organisms. The genome of the hydrogenotrophic methanogen Methanobrevibacter boviskoreani JH1^T harbors genes encoding an NADP-dependent alcohol dehydrogenase and a F₄₂₀-dependent NADP reductase, which may facilitate the transfer of reducing potential from ethanol to F₄₂₀ via NADP. The aim of this study was to explore the anaerobic culturing of JH1^T without pressurised H₂, using a variety of short chain alcohols. The results demonstrate that in the absence of H₂, JHI^T can use ethanol, 1-propanol, and 1-butanol but not methanol, as a source of reducing potential for methanogenesis. The ability to use ethanol to drive CH₄ formation in JH1^T makes it possible to develop a high throughput culture-based bioassay enabling screening of potential anti-methanogen compounds. The development of this resource will help researchers globally to accelerate the search for methane mitigation technologies for ruminant animals. Global emissions pathways that are consistent with the temperature goal of the Paris Agreement, rely on substantial reductions of agricultural greenhouse gasses, particularly from ruminant animals.

Hydrogen and formate production and utilisation in the rumen and the human colon

Molecular hydrogen (H₂) and formate (HCOO⁻) are metabolic end products of many primary fermenters in the mammalian gut. Both play a vital role in fermentation where they are electron sinks for individual microbes in an anaerobic environment that lacks external electron acceptors. If H₂ and/or formate accumulate within the gut ecosystem, the ability of primary fermenters to regenerate electron carriers may be inhibited and microbial metabolism and growth disrupted. Consequently, H₂- and/or formate-consuming microbes such as methanogens and homoacetogens play a key role in maintaining the metabolic efficiency of primary fermenters. There is increasing interest in identifying approaches to manipulate mammalian gut environments for the benefit of the host and the environment. As H₂ and formate are important mediators of interspecies interactions, an understanding of their production and utilisation could be a significant entry point for the development of successful interventions. Ruminant methane mitigation approaches are discussed as a model to help understand the fate of H₂ and formate in gut systems.

Electron flow: key to mitigating ruminant methanogenesis

Disposal of electrons generated during the fermentation of ingested feed is a fundamental feature of anaerobic microbial gut ecosystems. Here, we focus on the well-studied rumen environment to highlight how electrons are transferred through anaerobic fermentation pathways and how manipulating this electron flow is important to reducing methane emissions from ruminants. Priorities for research that can accelerate understanding in this area are highlighted.

Complete Genome Sequence of the Polysaccharide-Degrading Rumen Bacterium Pseudobutyrivibrio xylanivorans MA3014 Reveals an Incomplete Glycolytic Pathway

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.

Occurrence and expression of genes encoding methyl-compound production in rumen bacteria

Background

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.

Results

Screens of metagenomic data using the protein domains of choline trimethylamine lyase (CutC), and activator protein (CutD) found good matches only to Olsenella umbonata and to Caecibacter, while the Hungate1000 genomes and metagenome assembled genomes from the cattle rumen found bacteria within the phyla Actinobacteria, Firmicutes and Proteobacteria. The cutC and cutD genes clustered with genes that encode structural components of bacterial microcompartment proteins. Prevotella was the dominant genus encoding pectin methyl esterases, with smaller numbers of sequences identified from other fibre-degrading rumen bacteria. Some large pectin methyl esterases (> 2100 aa) were found to be encoded in Butyrivibrio genomes. The pectin-utilising, methane-producing consortium was composed of (i) a putative pectin-degrading bacterium (phylum Tenericutes, class Mollicutes), (ii) a galacturonate-using Sphaerochaeta sp. predicted to produce acetate, lactate, and ethanol, and (iii) a methylotrophic methanogen, Methanosphaera sp., with the ability to form methane via a primary ethanol-dependent, hydrogen-independent, methanogenesis pathway.

Conclusions

The main bacteria that produce methyl-compounds have been identified in ruminants. Their enzymatic activities can now be targeted with the aim of finding ways to reduce the supply of methyl-compound substrates to methanogens, and thereby limit methylotrophic methanogenesis in the rumen.

Applications of the Soil, Plant and Rumen Microbiomes in Pastoral Agriculture

The production of dairy, meat, and fiber by ruminant animals relies on the biological processes occurring in soils, forage plants, and the animals' rumens. Each of these components has an associated microbiome, and these have traditionally been viewed as distinct ecosystems. However, these microbiomes operate under similar ecological principles and are connected via water, energy flows, and the carbon and nitrogen nutrient cycles. Here, we summarize the microbiome research that has been done in each of these three environments (soils, forage plants, animals' rumen) and investigate what additional benefits may be possible through understanding the interactions between the various microbiomes. The challenge for future research is to enhance microbiome function by appropriate matching of plant and animal genotypes with the environment to improve the output and environmental sustainability of pastoral agriculture.

Diverse hydrogen production and consumption pathways influence methane production in ruminants

Farmed ruminants are the largest source of anthropogenic methane emissions globally. The methanogenic archaea responsible for these emissions use molecular hydrogen (H₂), 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 H₂ metabolism. Two-thirds of sequenced rumen bacterial and archaeal genomes encode enzymes that catalyse H₂ 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 H₂ 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 H₂ 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 H₂. These findings may help to inform ongoing strategies to mitigate methane emissions by increasing flux through alternative H₂ uptake pathways, including through animal selection, dietary supplementation and methanogenesis inhibitors.

Improved taxonomic assignment of rumen bacterial 16S rRNA sequences using a revised SILVA taxonomic framework

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.

Toward Understanding Phage:Host Interactions in the Rumen; Complete Genome Sequences of Lytic Phages Infecting Rumen Bacteria

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.

Gene and transcript abundances of bacterial type III secretion systems from the rumen microbiome are correlated with methane yield in sheep

Background

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 CH₄/kg DM/day) and aims to define differences in microbial communities, and in gene and transcript abundances that can explain the animal methane phenotype.

Methods

Rumen microbiome metagenomic and metatranscriptomic data were analysed by Gene Set Enrichment, sparse partial least squares regression and the Wilcoxon Rank Sum test to estimate correlations between specific KEGG bacterial pathways/genes and high methane yield in sheep. KEGG genes enriched in high methane yield sheep were reassembled from raw reads and existing contigs and analysed by MEGAN to predict their phylogenetic origin. Protein coding sequences from Succinivibrio dextrinosolvens strains were analysed using Effective DB to predict bacterial type III secreted proteins. The effect of S. dextrinosolvens strain H5 growth on methane formation by rumen methanogens was explored using co-cultures.

Results

Detailed analysis of the rumen microbiomes of high methane yield sheep shows that gene and transcript abundances of bacterial type III secretion system genes are positively correlated with methane yield in sheep. Most of the bacterial type III secretion system genes could not be assigned to a particular bacterial group, but several genes were affiliated with the genus Succinivibrio, and searches of bacterial genome sequences found that strains of S. dextrinosolvens were part of a small group of rumen bacteria that encode this type of secretion system. In co-culture experiments, S. dextrinosolvens strain H5 showed a growth-enhancing effect on a methanogen belonging to the order Methanomassiliicoccales, and inhibition of a representative of the Methanobrevibacter gottschalkii clade.

Conclusions

This is the first report of bacterial type III secretion system genes being associated with high methane emissions in ruminants, and identifies these secretions systems as potential new targets for methane mitigation research. The effects of S. dextrinosolvens on the growth of rumen methanogens in co-cultures indicate that bacteria-methanogen interactions are important modulators of methane production in ruminant animals.

Electronic supplementary material

The online version of this article (doi:10.1186/s13104-017-2671-0) contains supplementary material, which is available to authorized users.

Dynamics and genetic diversification of Escherichia coli during experimental adaptation to an anaerobic environment

Background

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.

Methods

Twenty-one populations of E. coli REL4536, an aerobically evolved 10,000th generation descendent of the E. coli B strain, REL606, were established from a clonal ancestral culture. These were serially sub-cultured for 2,000 generations in a defined minimal glucose medium in strict aerobic and strict anaerobic environments, as well as in a treatment that fluctuated between the two environments. The competitive fitness of the evolving lineages was assessed at approximately 0, 1,000 and 2,000 generations, in both the environment of selection and the alternative environment. Whole genome re-sequencing was performed on random colonies from all lineages after 2,000-generations. Mutations were identified relative to the ancestral genome, and based on the extent of parallelism, traits that were likely to have contributed towards adaptation were inferred.

Results

There were increases in fitness relative to the ancestor among anaerobically evolved lineages when tested in the anaerobic environment, but no increases were found in the aerobic environment. For lineages that had evolved under the fluctuating regime, relative fitness increased significantly in the anaerobic environment, but did not increase in the aerobic environment. The aerobically-evolved lineages did not increase in fitness when tested in either the aerobic or anaerobic environments. The strictly anaerobic lineages adapted more rapidly to the anaerobic environment than did the fluctuating lineages. Two main strategies appeared to predominate during adaptation to the anaerobic environment: modification of energy generation pathways, and inactivation of non-essential functions. Fermentation pathways appeared to alter through selection for mutations in genes such as nadR, adhE, dcuS/R, and pflB. Mutations were frequently identified in genes for presumably dispensable functions such as toxin-antitoxin systems, prophages, virulence and amino acid transport. Adaptation of the fluctuating lineages to the anaerobic environments involved mutations affecting traits similar to those observed in the anaerobically evolved lineages.

Discussion

There appeared to be strong selective pressure for activities that conferred cell yield advantages during anaerobic growth, which include restoring activities that had previously been inactivated under long-term continuous aerobic evolution of the ancestor.

Anaerobically Grown Escherichia coli Has an Enhanced Mutation Rate and Distinct Mutational Spectra

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 complete genome sequence of the methanogenic archaeon ISO4-H5 provides insights into the methylotrophic lifestyle of a ruminal representative of the Methanomassiliicoccales

Methane emissions from agriculture represent around 9 % of global anthropogenic greenhouse emissions. The single largest source of this methane is animal enteric fermentation, predominantly from ruminant livestock where it is produced mainly in their fermentative forestomach (or reticulo-rumen) by a group of archaea known as methanogens. In order to reduce methane emissions from ruminants, it is necessary to understand the role of methanogenic archaea in the rumen, and to identify their distinguishing characteristics that can be used to develop methane mitigation technologies. To gain insights into the role of methylotrophic methanogens in the rumen environment, the genome of a methanogenic archaeon has been sequenced. This isolate, strain ISO4-H5, was isolated from the ovine rumen and belongs to the order Methanomassiliicoccales. Genomic analysis suggests ISO4-H5 is an obligate hydrogen-dependent methylotrophic methanogen, able to use methanol and methylamines as substrates for methanogenesis. Like other organisms within this order, ISO4-H5 does not possess genes required for the first six steps of hydrogenotrophic methanogenesis. Comparison between the genomes of different members of the order Methanomassiliicoccales revealed strong conservation in energy metabolism, particularly in genes of the methylotrophic methanogenesis pathway, as well as in the biosynthesis and use of pyrrolysine. Unlike members of Methanomassiliicoccales from human sources, ISO4-H5 does not contain the genes required for production of coenzyme M, and so likely requires external coenzyme M to survive.

The complete genome sequence of the rumen methanogen Methanobrevibacter millerae SM9

Methanobrevibacter millerae SM9 was isolated from the rumen of a sheep maintained on a fresh forage diet, and its genome has been sequenced to provide information on the phylogenetic diversity of rumen methanogens with a view to developing technologies for methane mitigation. It is the first rumen isolate from the Methanobrevibacter gottschalkii clade to have its genome sequence completed. The 2.54 Mb SM9 chromosome has an average G + C content of 31.8 %, encodes 2269 protein-coding genes, and harbors a single prophage. The overall gene content is comparable to that of Methanobrevibacter ruminantium M1 and the type strain of M. millerae (ZA-10^T) suggesting that the basic metabolism of these two hydrogenotrophic rumen methanogen species is similar. However, M. millerae has a larger complement of genes involved in methanogenesis including genes for methyl coenzyme M reductase II (mrtAGDB) which are not found in M1. Unusual features of the M. millerae genomes include the presence of a tannase gene which shows high sequence similarity with the tannase from Lactobacillus plantarum, and large non-ribosomal peptide synthase genes. The M. millerae sequences indicate that methane mitigation strategies based on the M. ruminantium M1 genome sequence are also likely to be applicable to members of the M. gottschalkii clade.

The complete genome sequence of Eubacterium limosum SA11, a metabolically versatile rumen acetogen

Acetogens are a specialized group of anaerobic bacteria able to produce acetate from CO2 and H2 via the Wood-Ljungdahl pathway. In some gut environments acetogens can compete with methanogens for H2, and as a result rumen acetogens are of interest in the development of microbial approaches for methane mitigation. The acetogen Eubacterium limosum SA11 was isolated from the rumen of a New Zealand sheep and its genome has been sequenced to examine its potential application in methane mitigation strategies, particularly in situations where hydrogenotrophic methanogens are inhibited resulting in increased H2 levels in the rumen. The 4.15 Mb chromosome of SA11 has an average G + C content of 47 %, and encodes 3805 protein-coding genes. There is a single prophage inserted in the chromosome, and several other gene clusters appear to have been acquired by horizontal transfer. These include genes for cell wall glycopolymers, a type VII secretion system, cell surface proteins and chemotaxis. SA11 is able to use a variety of organic substrates in addition to H2/CO2, with acetate and butyrate as the principal fermentation end-products, and genes involved in these metabolic pathways have been identified. An unusual feature is the presence of 39 genes encoding trimethylamine methyltransferase family proteins, more than any other bacterial genome. Overall, SA11 is a metabolically versatile organism, but its ability to grow on such a wide range of substrates suggests it may not be a suitable candidate to take the place of hydrogen-utilizing methanogens in the rumen.

The Cytosolic Oligosaccharide-Degrading Proteome of Butyrivibrio Proteoclasticus

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.

The complete genome sequence of the rumen methanogen Methanosarcina barkeri CM1

Methanosarcina species are the most metabolically versatile of the methanogenic Archaea and can obtain energy for growth by producing methane via the hydrogenotrophic, acetoclastic or methylotrophic pathways. Methanosarcina barkeri CM1 was isolated from the rumen of a New Zealand Friesian cow grazing a ryegrass/clover pasture, and its genome has been sequenced to provide information on the phylogenetic diversity of rumen methanogens with a view to developing technologies for methane mitigation. The 4.5 Mb chromosome has an average G + C content of 39 %, and encodes 3523 protein-coding genes, but has no plasmid or prophage sequences. The gene content is very similar to that of M. barkeri Fusaro which was isolated from freshwater sediment. CM1 has a full complement of genes for all three methanogenesis pathways, but its genome shows many differences from those of other sequenced rumen methanogens. Consequently strategies to mitigate ruminant methane need to include information on the different methanogens that occur in the rumen.

The complete genome sequence of the rumen methanogen Methanobacterium formicicum BRM9

Methanobacterium formicicum BRM9 was isolated from the rumen of a New Zealand Friesan cow grazing a ryegrass/clover pasture, and its genome has been sequenced to provide information on the phylogenetic diversity of rumen methanogens with a view to developing technologies for methane mitigation. The 2.45 Mb BRM9 chromosome has an average G + C content of 41%, and encodes 2,352 protein-coding genes. The genes involved in methanogenesis are comparable to those found in other members of the Methanobacteriaceae with the exception that there is no [Fe]-hydrogenase dehydrogenase (Hmd) which links the methenyl-H4MPT reduction directly with the oxidation of H2. Compared to the rumen Methanobrevibacter strains, BRM9 has a much larger complement of genes involved in determining oxidative stress response, signal transduction and nitrogen fixation. BRM9 also has genes for the biosynthesis of the compatible solute ectoine that has not been reported to be produced by methanogens. The BRM9 genome has a prophage and two CRISPR repeat regions. Comparison to the genomes of other Methanobacterium strains shows a core genome of ~1,350 coding sequences and 190 strain-specific genes in BRM9, most of which are hypothetical proteins or prophage related.

Determining the culturability of the rumen bacterial microbiome

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.

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 CH₄ mitigation at the levels of microbiota composition and transcriptional regulation.

Metasecretome-selective phage display approach for mining the functional potential of a rumen microbial community

Background

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.

Results

By combining secretome-selective phage display and next-generation sequencing, we focused the sequence analysis of complex rumen microbial community on the metasecretome component of the metagenome. This approach achieved high enrichment (29 fold) of secreted fibrolytic enzymes from the plant-adherent microbial community of the bovine rumen. In particular, we identified hundreds of heretofore rare modules belonging to cellulosomes, cell-surface complexes specialised for recognition and degradation of the plant fibre.

Conclusions

As a method, metasecretome phage display combined with next-generation sequencing has a power to sample the diversity of low-abundance surface and secreted proteins that would otherwise require exceptionally large metagenomic sequencing projects. As a resource, metasecretome display library backed by the dataset obtained by next-generation sequencing is ready for i) affinity selection by standard phage display methodology and ii) easy purification of displayed proteins as part of the virion for individual functional analysis.

Electronic supplementary material

The online version of this article (doi:10.1186/1471-2164-15-356) contains supplementary material, which is available to authorized users.

Interaction between the genomes of Lactococcus lactis and phages of the P335 species

Phages of the P335 species infect Lactococcus lactis and have been particularly studied because of their association with strains of L. lactis subsp. cremoris used as dairy starter cultures. Unlike other lactococcal phages, those of the P335 species may have a temperate or lytic lifestyle, and are believed to originate from the starter cultures themselves. We have sequenced the genome of L. lactis subsp. cremoris KW2 isolated from fermented corn and found that it contains an integrated P335 species prophage. This 41 kb prophage (Φ KW2) has a mosaic structure with functional modules that are highly similar to several other phages of the P335 species associated with dairy starter cultures. Comparison of the genomes of 26 phages of the P335 species, with either a lytic or temperate lifestyle, shows that they can be divided into three groups and that the morphogenesis gene region is the most conserved. Analysis of these phage genomes in conjunction with the genomes of several L. lactis strains shows that prophage insertion is site specific and occurs at seven different chromosomal locations. Exactly how induced or lytic phages of the P335 species interact with carbohydrate cell surface receptors in the host cell envelope remains to be determined. Genes for the biosynthesis of a variable cell surface polysaccharide and for lipoteichoic acids (LTAs) are found in L. lactis and are the main candidates for phage receptors, as the genes for other cell surface carbohydrates have been lost from dairy starter strains. Overall, phages of the P335 species appear to have had only a minor role in the adaptation of L. lactis subsp. cremoris strains to the dairy environment, and instead they appear to be an integral part of the L. lactis chromosome. There remains a great deal to be discovered about their role, and their contribution to the evolution of the bacterial genome.

The Complete Genome Sequence of Methanobrevibacter sp. AbM4

Methanobrevibacter sp. AbM4 was originally isolated from the abomasal contents of a sheep and was chosen as a representative of the Methanobrevibacter wolinii clade for genome sequencing. The AbM4 genome is smaller than that of the rumen methanogen M. ruminantium M1 (2.0 Mb versus 2.93 Mb), encodes fewer open reading frames (ORFs) (1,671 versus 2,217) and has a lower G+C percentage (29% versus 33%). Overall, the composition of the AbM4 genome is very similar to that of M1 suggesting that the methanogenesis pathway and central metabolism of these strains are highly similar, and both organisms are likely to be amenable to inhibition by small molecule inhibitors and vaccine-based methane mitigation technologies targeting these conserved features. The main differences compared to M1 are that AbM4 has a complete coenzyme M biosynthesis pathway and does not contain a prophage or non-ribosomal peptide synthase genes. However, AbM4 has a large CRISPR region and several type I and type II restriction-modification system components. Unusually, DNA-directed RNA polymerase B' and B'' subunits of AbM4 are joined, a feature only previously observed in some thermophilic archaea. AbM4 has a much reduced complement of genes encoding adhesin-like proteins which suggests it occupies a ruminal niche different from that of M1.

Carbohydrate transporting membrane proteins of the rumen bacterium, Butyrivibrio proteoclasticus

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.

Extracellular polysaccharide-degrading proteome of Butyrivibrio proteoclasticus

Plant polysaccharide-degrading rumen microbes are fundamental to the health and productivity of ruminant animals. Butyrivibrio proteoclasticus B316(T) is a gram-positive, butyrate-producing anaerobic bacterium with a key role in hemicellulose degradation in the rumen. Gel-based proteomics was used to examine the growth-phase-dependent abundance patterns of secreted proteins recovered from cells grown in vitro with xylan or xylose provided as the sole supplementary carbon source. Five polysaccharidases and two carbohydrate-binding proteins (CBPs) were among 30 identified secreted proteins. The endo-1,4-β-xylanase Xyn10B was 17.5-fold more abundant in the culture medium of xylan-grown cells, which suggests it plays an important role in hemicellulose degradation. The secretion of three nonxylanolytic enzymes and two CBPs implies they augment hemicellulose degradation by hydrolysis or disruption of associated structural polysaccharides. Sixteen ATP-binding cassette (ABC) transporter substrate-binding proteins were identified, several of which had altered relative abundance levels between growth conditions, which suggests they are important for oligosaccharide uptake. This study demonstrates that B. proteoclasticus modulates the secretion of hemicellulose-degrading enzymes and ATP-dependent sugar uptake systems in response to growth substrate and supports the notion that this organism makes an important contribution to polysaccharide degradation in the rumen.

Functional genome analysis of Bifidobacterium breve UCC2003 reveals type IVb tight adherence (Tad) pili as an essential and conserved host-colonization factor

Development of the human gut microbiota commences at birth, with bifidobacteria being among the first colonizers of the sterile newborn gastrointestinal tract. To date, the genetic basis of Bifidobacterium colonization and persistence remains poorly understood. Transcriptome analysis of the Bifidobacterium breve UCC2003 2.42-Mb genome in a murine colonization model revealed differential expression of a type IVb tight adherence (Tad) pilus-encoding gene cluster designated "tad(2003)." Mutational analysis demonstrated that the tad(2003) gene cluster is essential for efficient in vivo murine gut colonization, and immunogold transmission electron microscopy confirmed the presence of Tad pili at the poles of B. breve UCC2003 cells. Conservation of the Tad pilus-encoding locus among other B. breve strains and among sequenced Bifidobacterium genomes supports the notion of a ubiquitous pili-mediated host colonization and persistence mechanism for bifidobacteria.

Transposition of Tn916 in the four replicons of the Butyrivibrio proteoclasticus B316(T) genome

The rumen bacterium Butyrivibrio proteoclasticus B316(T) has a 4.4-Mb genome composed of four replicons (approximately 3.55 Mb, 361, 302 and 186 kb). Mutagenesis of B316(T) was performed with the broad host-range conjugative transposon Tn916 to screen for functionally important characteristics. The insertion sites of 123 mutants containing a single copy of Tn916 were identified and corresponded to 53 different insertion points, of which 18 (34.0%), representing 39 mutants (31.7%), were in ORFs and 12 were where transposition occurred in both directions (top and bottom DNA strand). Up to eight mutants from several independent conjugation experiments were found to have the same integration site. Although transposition occurred in all four replicons, the number of specific insertion sites, transposition frequency and the average intertransposon distance between insertions varied between the four replicons. In silico analysis of the 53 insertion sites was used to model a target consensus sequence for Tn916 integration into B316(T) . A search of the B316(T) genome using the modelled target consensus sequence (up to two mismatches) identified 39 theoretical Tn916 insertion sites (19 coding, 20 noncoding), of which nine corresponded to Tn916 insertions identified in B316(T) mutants during our conjugation experiments.

Chromosomal diversity in Lactococcus lactis and the origin of dairy starter cultures

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.

The glycobiome of the rumen bacterium Butyrivibrio proteoclasticus B316(T) highlights adaptation to a polysaccharide-rich environment

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.4Mb 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.

The genome sequence of the rumen methanogen Methanobrevibacter ruminantium reveals new possibilities for controlling ruminant methane emissions

BACKGROUND

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.

METHODOLOGY/PRINCIPAL FINDINGS

The M1 genome was sequenced, annotated and subjected to comparative genomic and metabolic pathway analyses. Conserved and methanogen-specific gene sets suitable as targets for vaccine development or chemogenomic-based inhibition of rumen methanogens were identified. The feasibility of using a synthetic peptide-directed vaccinology approach to target epitopes of methanogen surface proteins was demonstrated. A prophage genome was described and its lytic enzyme, endoisopeptidase PeiR, was shown to lyse M1 cells in pure culture. A predicted stimulation of M1 growth by alcohols was demonstrated and microarray analyses indicated up-regulation of methanogenesis genes during co-culture with a hydrogen (H(2)) producing rumen bacterium. We also report the discovery of non-ribosomal peptide synthetases in M. ruminantium M1, the first reported in archaeal species.

CONCLUSIONS/SIGNIFICANCE

The M1 genome sequence provides new insights into the lifestyle and cellular processes of this important rumen methanogen. It also defines vaccine and chemogenomic targets for broad inhibition of rumen methanogens and represents a significant contribution to worldwide efforts to mitigate ruminant methane emissions and reduce production of anthropogenic greenhouse gases.