- Research article
Gene and transcript abundances of bacterial type III secretion systems from the rumen microbiome are correlated with methane yield in sheep
BMC Research Notesvolume 10, Article number: 367 (2017)
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.
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.
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.
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.
Methane emissions from enteric fermentation in livestock are important contributors to global warming [1, 2] and most of these emissions are from the action of methanogenic archaea in the reticulo-rumen of ruminant animals. In the reticulo-rumen, methanogens use the products of bacterial fermentation, (hydrogen, carbon dioxide, formate and methyl compounds) to produce methane [3, 4]. Considerable amounts of methane are produced each day, a typical sheep produces ~40 L of methane per day, while a dairy cow belches around 450 L each day [5,6,7]. Methane is a particularly strong greenhouse gas with a global warming potential of 34× that of CO2 , and many researchers around the world are investigating ways to reduce methane emissions from ruminants. Animal selection is being investigated as a means of breeding low methane emitting animals [8,9,10,11,12,13,14], and sheep with consistently high methane yields (HMY; defined in g methane/kg dry matter intake DMI/day) or low methane yield (LMY) have been identified. The methane yield trait is heritable , and the animal characteristics controlling feed particle retention time [15,16,17,18] and rumen volume  are likely to contribute to the trait. Microorganisms in the rumen also contribute to the methane yield phenotype, and differences between HMY and LMY sheep in rumen microbiome composition  and methanogen gene expression  have been reported. We have recently found that methanogen changes in LMY animals are paralleled by an enrichment in Sharpea spp. in the bacterial community and higher abundance of the genes and transcripts encoding lactate formation and metabolism pathways that lead to lower hydrogen production and therefore lower methane formation in the rumen . Here, we report the analysis of the rumen microbiomes of HMY animals and report the finding that gene and transcript abundances of bacterial type III secretion system (T3SS) genes are positively correlated with HMY in sheep.
Study aim and design
This report is part of a larger study , the aim of which was to characterise the rumen microbiomes of sheep which differed in their methane yield (g CH4/kg DMI/day) to investigate whether there were corresponding differences in microbial communities, and in gene and transcript abundances that could explain the animal methane phenotype. Rams (11 HMY and 11 LMY) were selected based on their methane yields and breeding values  and their methane yields were re-measured twice (measurements were separated by 2 weeks) in respiration chambers at the New Zealand Ruminant Methane Measurement Centre, AgResearch Grasslands, Palmerston North, New Zealand, on a pelleted lucerne diet . Samples of rumen contents were collected at the end of each methane measurement period, by stomach intubation 4 h after the morning feeding. The pH of rumen contents was measured and the samples were immediately snap-frozen as pellets in liquid N2, and stored at −85 °C prior to DNA and RNA extractions. Based on the re-measured methane yields, frozen rumen samples were selected from 4 HMY (group mean 15.85 g CH4/kg DMI/day) and 4 LMY (group mean 11.44 g CH4/kg DMI/day) animals (one sample each per methane measurement, 16 samples), and from 2 animals which had intermediate methane yield (IMY, group mean 13.77 g CH4/kg DMI/day) (one each at the two methane measurement points, 4 samples) and used for DNA and RNA extractions (20 samples for each extraction). Nucleic acids extractions, purifications and library construction followed protocols described in the parent study . Briefly, DNA was extracted using the “Repeated Bead Beating and Column (RBB + C) purification” method . For large paired-end insert libraries, high-molecular-weight DNA was extracted using the method of Rosewarne et al. . RNA extraction was via a hot lysis-acid phenol extraction method . For cDNA library construction, total RNA was enriched for mRNA using the Ribo-Zero TM rRNA Removal Kit (Meta-Bacteria, Epicenter Biotechnologies, Madison, WI, USA), and fragmented using mRNA Fragmentation Reagents (Ambion, Foster City, CA, USA). Double-stranded cDNA (ds cDNA) was synthesised using SuperScript II reverse transcriptase (Invitrogen, Carlsbad, CA, USA), using random hexamers (MBI Fermentas, NY, USA). The cDNA sequencing libraries were generated and amplified using the Illumina TruSeqTM genomic sample prep kit (Illumina, San Diego, CA, USA) following the manufacturer’s instructions. The amplified libraries were purified and size-selected and the pooled library was sequenced using the Illumina HiSeq 2000 platform. The metagenomic sequencing produced, on average, 216 M reads per sample (51 Gb sequence), while the metatranscriptome sequencing produced 35 M reads (5.3 Gb non-rRNA sequence) per sample. All sequencing data used in this study has been previously published by our group and detailed description of data processing can be found in these publications [20, 21]. Briefly, quality-checked 16S rRNA amplicon sequencing data was processed using the QIME pipeline  and phylogenetically assigned using the Greengenes database (version gg_13_5) . Metagenomic and metatranscriptome sequencing datasets were quality and artifact filtered and paired end reads merged using FLASH  and used for read count based comparative analysis by screening against the KEGG database , (version 58.1) using USEARCH 6.0  at an E-value cutoff of 1 × 10−5 and read count matrices were constructed and normalised to reads per million (RPM).
Analysis of fermentation acids from rumen samples was performed on a Shimadzu 5050a GC–MS (Shimadzu, Kyoto, Japan) equipped with a ZB-5 MS (30 m × 0.25 mm ID × 0.25 µm film thickness, Phenomenex, USA) capillary column. Acids were extracted from the acidified samples with ether and derivatised to their t-butyl-dimethylsilyl esters, which detects both volatile and non-volatile fatty acids.
Statistical analyses were conducted using reads per million (RPM)-normalised read count matrices as previously described . Gene set enrichment analysis (GSEA, ) was used to estimate differential gene and transcript abundances between HMY or LMY animals (n = 16) at the pathway level, based on KEGG pathways. The GSEA-P application  was used to pre-rank genes based on signal-to-noise metric scores and to estimate normalised enrichment scores (NES), nominal P values and false discovery rates (FDR) with 10,000 phenotype permutations. Sparse partial least squares regression analysis (sPLS) in R  was used to estimate correlations between specific KEGG genes and methane yield, using the spls package . Genes selected as methane predictors were assessed using the optimum sparsity tuning parameter (eta) and number of hidden components (K) predicted in the mean squared prediction error plot (MSPE). KEGG gene predictors with confidence intervals (95%) shifting from positive to negative correlation or from negative to positive, were excluded from the predictor gene set. The selected KEGG gene sets of the sPLS regression analysis of each dataset were manually compared with the results of categorical statistical analysis using Wilcoxon Rank Sum (WRS) test (10,000 permutations)  and GSEA pathways enrichment scores. For all statistical analyses a significance threshold of Benjamini Hochberg (BH) corrected P < 0.05 was used. Subsequent analyses concentrated on KEGG gene sets or gene categories supported by two or more methods of analysis.
Read extraction and assembly
To gain an overview of the taxonomic origin of bacterial T3SS genes we reassembled these genes based on raw reads and contigs from existing assemblies, with hits to any of the KEGG genes encoding for the relevant subunits (K03219, K03221–K03230, K04056–K04059) from both metagenome and metatranscriptome data. Individual assemblies based on the 20 metagenome and 20 metatranscriptome datasets were constructed as previously described . Assemblies were combined, and duplicated contigs, or smaller contigs covered by larger ones, were removed using the clustering function in Vmatch (http://www.vmatch.de). The resulting contigs were submitted to IMG/Mer  for gene-calling and automatic annotation. All assembled genes with hits to the relevant KEGG genes were extracted and combined with the corresponding reads from metagenome and metatranscriptome data. Contigs were extended using the Distributed Nucleating Assembler function in Kmernator (https://github.com/JGI-Bioinformatics/Kmernator). Genes on the resulting contigs were predicted using MetaGeneMark .
Phylogenetic assignments of bacterial T3SS genes and prediction of effector proteins
Re-assembled contigs that contained partial and full length T3SS genes were imported into Geneious v.7.0.5 (Biomatters Ltd, Auckland, New Zealand) and de novo assembled into larger contigs where possible. The resulting 73 contigs were BLAST searched against the NCBI non-redundant database using BLASTX . The output was analysed using MEGAN5  for phylogenetic assignment using default settings. Matches were confirmed by homology searches of relevant genes in the IMG/M database . Protein coding sequences from the draft genomes of S. dextrinosolvens strains H5, 22B and ACV-10 were analysed using EffectiveDB , which uses various software tools to predict bacterial secreted proteins based on their amino acid sequences. These tools include EffectiveT3 for prediction of Type III secretion signals, EffectiveCCBD for detection of conserved binding domains of Type III chaperones, EffectiveELD for secretion system independent prediction of secreted proteins based on eukaryotic-like domains, T4SEpre for recognition of Type IV secreted proteins and Predotar to screen N-terminal targeting sequences to predict their subcellular localization in eukaryotic host cells.
Cultures of rumen methanogens were inoculated with Succinivibrio dextrinosolvens strain H5 to explore its effects on methane production in vitro. Two Methanobrevibacter species were used, representing hydrogenotrophic methanogens within the Mbb. ruminantium clade (Mbb. ruminantium M1) or Mbb. gottschalkii clade (Mbb. millerae SM9), along with two methylotrophs representing the Methanosphaera clade (Methanosphaera sp. ISO3-F5) and the order Methanomassiliicoccales, (methanogenic archaeon ISO4-H5). Each methanogen was grown with S. dextrinosolvens strain H5 in either RM02 medium  with hydrogen + carbon dioxide (for the hydrogenotrophic methanogens, 1 atm over-pressure of a 80:20, hydrogen:carbon dioxide) or BY medium  with methanol + hydrogen (1 atm over pressure of 80:20, hydrogen:carbon dioxide + 20 mM methanol). Both media contained acetate (20 mM) and coenzyme M (1 mM). The hydrogenotrophic and methylotrophic methanogen cultures were grown in their respective media until methane was first detected in the culture headspace, then were inoculated with 10% inoculum S. dextrinosolvens strain H5 with glucose (10 mM final concentration) or pectin (1% pectin final concentration) added as a substrate for the S. dextrinosolvens, respectively. Incubation was continued until maximal methane production was detected in the control tubes. Control tubes contained methanogen cultures which received only the growth substrate without the S. dextrinosolvens inoculum. The pHs of the cultures were measured after inoculation and after completion of growth using pH indicator strips (Merck KGaA, Darmstadt, Germany). Methane and hydrogen concentrations in the culture headspaces were determined by removing a sample of the gases and analysing by gas chromatography (Aerograph 660, Wilkins Instruments & Research Inc., Walnut Creek, CA, USA) against methane and hydrogen standards. Samples of cultures were taken before and after incubation for analysis by gas–liquid chromatography to determine production or use of volatile fatty acids. Samples were centrifuged (21,000×g at 4 °C for 10 min) and 0.9 mL of the supernatant was added to 0.1 mL of internal standard solution (19 mM ethyl butyrate in 20% (v/v) phosphoric acid). Samples were kept at −20 °C until analysis, when they were thawed, clarified by centrifugation (21,000×g at 4 °C for 10 min) and 0.8 mL of the supernatant was transferred into a 2-mL crimp cap vial for analysis by gas–liquid chromatography. Supernatant samples were analysed on a nitroterephthalic acid-modified polyethylene glycol column (DB-FFAP; 30 m × 0.53 mm × 1.0 µm film thickness; J & W Scientific, Folsom, CA, USA) attached to a Hewlett-Packard 6890 series gas chromatography system, using helium as the carrier gas (5 mL/min). The oven temperature started at 85 °C, ramped to 200 °C at 10 °C/min, was held at 200 °C for 10 min, and then decreased to 50 °C and held for 5 min before the next sample was injected. Peaks were detected with a flame ionization detector, identified by comparison with standards, and integrated with Hewlett-Packard ChemStation software (version 4.02).
Animal ethics approval
The collections of rumen contents from sheep were carried out under the approval of the AgResearch Ltd Grasslands Animal Ethics Committee (Approval 13606).
Consent to participate and publish
T3SS genes and transcripts are enriched in HMY animals
The metagenomic and metatranscriptomic data generated from rumen samples collected from LMY and HMY sheep were analysed by sPLS and the WRS test, and 60 genes and 36 transcripts were significantly (BH corrected P < 0.05) more abundant in the HMY animals in both analyses (Additional file 1). Those genes and transcripts that were twofold or greater more abundant in the HMY animals were enriched for KEGG genes related to bacterial T3SS subunits and are listed in Table 1. In the WRS test, 12 out of a total of 15 subunits known to be involved in bacterial T3SSs had significantly more reads in the HMY animals. Ten T3SS subunits were also identified as correlation predictors in the sPLS analysis (Table 2) with an adjusted multiple regression coefficient of R2 = 0.76, (P = 3.48 × 10−7, Fig. 1). At the transcriptome level, a similar but weaker trend was observed, with two subunits showing significantly higher transcript abundance in HMY animals and one subunit being selected as a predictor gene for methane yield in the sPLS regression analysis (Table 2). The GSEA identified Bacterial secretion systems (ko03070) as the top ranked pathway represented in the metagenome dataset associated with HMY animals (Table 3) along with Drug metabolism and other enzymes (ko00983) and Alanine, aspartate and glutamate metabolism (ko00250) which showed a similar trend. Hierarchical clustering analysis based on Z-scores for both metagenome and metatranscriptome data, provided further support for differential T3SS gene abundance and expression between HMY and LMY animals (Fig. 2).
Bacterial T3SS genes in the rumen are associated with unidentified Proteobacteria and the genus Succinivibrio
To identity the origin of the T3SS genes found in the rumen metagenome and metatranscriptome sequences, we screened for contigs containing these genes, and a total of 101 genes were found on the 73 contigs retrieved. Using the lowest common ancestor algorithm in MEGAN, 82% of contigs were assigned to Bacteria, while the remaining 18% could not be assigned to any phylogenetic grouping. Some Bacteria-associated contigs could be further assigned (Fig. 3), with the largest group falling into the Proteobacteria (40%, 24 contigs). The largest genus-level group within the Proteobacteria was Succinivibrio (9 contigs, 12%). The phylogenetic affiliations of the Succinivibrio contigs were explored using manual BLAST searches against the IMG database. T3SS genes were identified in draft genomes of Succinivibrio dextrinosolvens strains H5, 22B, ACV-10, and DSM 3072 (Additional file 2). Further searches for T3SS KOs in bacteria isolated from the rumen showed that T3SS genes were also present in Succinimonas amylolytica DSM 2873 (type III secretion genes yscC, J, L, N, R, S, T, U, and V) and Desulfotomaculum ruminis DSM 2154 (yscN, ATP synthase type III secretion protein N [EC:184.108.40.206]) but in no other rumen bacteria.
Analysis of the scaffold sequences of the draft genomes of S. dextrinosolvens strains H5, 22B, ACV-10, and DSM 3072, showed that the genes encoding the T3SS structural subunits are arranged in conserved structures (Fig. 4) with several other genes known to be associated with T3SSs. These include genes encoding low calcium response chaperones (LcrH/SycD), several Tir chaperone proteins (CesT), secretion system effector C-like family protein (SseC; involved in translocon formation) and TyeA proteins (involved in translocation of Yersinia outer proteins into eukaryotic cells). Analysis of the amino acid sequences of the predicted protein coding sequences of S. dextrinosolvens strains H5, 22B and ACV-10 using the EffectiveDB suite of software, identified proteins predicted to be secreted by T3SSs, proteins containing conserved chaperone binding domains in their N-terminal regions, eukaryotic-like N-terminal signal sequences, and also uncovered several T4SS effector proteins in each of the strains (summarised in Table 4, individual proteins are shown in Additional file 3).
Succinivibrio dextrinosolvens stimulates methane formation in the Methanomassiliicoccales affiliated methanogenic archaeon ISO4-H5
The preceding analyses indicated that T3SSs, some of which were similar to those found in ruminal strains of Succinivibrio spp., were more abundant within the microbiomes of HMY versus LMY animals, and indicate a possible association of Succinivibrio with the HMY phenotype. Therefore experiments were performed to explore the effects of S. dextrinosolvens strain H5 on methane formation by common hydrogenotrophic and methylotrophic methanogens from the rumen. Inoculation of S. dextrinosolvens H5 into methanogen cultures stimulated methane formation in the methanogenic archaeon sp. ISO4-H5 (Fig. 5). The stimulation of methane formation was observed 24 h after S. dextrinosolvens inoculation and became significant after 48 h of co-culture. S. dextrinosolvens inoculation had a slight inhibitory effect on Mbb. millerae SM9, reducing methane by 21% at 96 h relative to the control (P < 0.05). S. dextrinosolvens inoculation had no effect on methane formation in either Mbb. ruminantium M1 or Methanosphaera sp. ISO3-F5 (Fig. 5). Formate was the main VFA produced in the co-cultures (Table 5).
The reduction of methane emissions from ruminants is being addressed via several lines of research, including the selection of animals based on their methane yields . As the methane emission trait is heritable in sheep, there is scope to select for flocks with low methane emissions, and such breeding programmes need to understand the mechanisms underlying the trait to enhance the selection process and to avoid negative impacts on digestive processes and animal productivity. We have shown previously that the microbiome of LMY sheep has small differences in methanogen communities  but large differences in the expression of genes encoding the methanogenesis pathway . We have also recently shown that the bacterial community in LMY animals is enriched for Sharpea spp. and expresses metabolic pathways that lead to lower hydrogen production and therefore lower methane formation. Here, the microbiomes of HMY animals have been analysed and we report the unusual finding that gene and transcript abundances of bacterial T3SS genes are positively correlated with HMY in sheep.
The enrichment of T3SS genes in rumen metagenome sequences from HMY animals is surprising as T3SSs are uncommon in the rumen microorganisms that have been characterised to date. T3SS are most commonly found in pathogenic, Gram-negative bacteria where they are known as “injectisomes”, because they encode a protein complex that forms a needle-like appendage  which is used during infection to translocate effector proteins into the host to promote pathogen survival and resistance to the host immune system. The protein subunits which are typically found in T3SSs are the cytoplasmic subunits (SctQ, L, O, and N which make up the cytoplasmic ring protein, stator, stalk, accessory protein, and ATPase respectively); the export apparatus (SctV, U, R, S, and T which are the major export protein, switch protein, three minor export proteins, respectively) and the base, and needle proteins (SctJ, D, I, C, P, X and F which are the inner Membrane Supramembrane (MS) ring, the outer MS ring, the inner rod, the secretin, the needle length regulator, a secretion protein and the needle filament, respectively). Representatives of nearly all the genes encoding subunits required for T3SS assembly were identified in the metagenome of HMY animals, including the cytoplasmic components (SctQLON), the export proteins (SctVURST) and the base and needle proteins (SctJDICF) and their occurrence and level of statistical support are shown in Fig. 6. Interestingly, genes encoding the needle protein (SctF), the needle length regulator (SctP) and the putative animal-specific secretion protein (SctX) were present at very low levels and were not enriched in HMY animals. T3SS in human pathogenic bacteria are also known to allow invasion of other hosts, including free-living amoebae and protozoa, which forms important environmental reservoirs for these pathogens and may protect the internalized bacteria from detection and treatment with biocides . Although most common in pathogens, T3SS are also found in bacteria that engage in symbiotic relationships with their animal or plant hosts, and in this context the role of T3SS is thought to be in determining the specificity and maintenance of the symbiotic interaction .
The majority of T3SS genes found in the metagenome dataset could not be assigned a particular phylogenetic origin, and appear to belong to unidentified and/or uncultured members of the rumen microbiome. Of those sequences that could be assigned, most were affiliated with the phylum Proteobacteria, and at the genus level to Succinivibrio, which includes as its only species the succinate-producing S. dextrinosolvens . Overall, the occurrence of genes encoding T3SSs in rumen bacteria was low; using available genome sequences, T3SS genes were only found in four strains of Succinivibrio dextrinosolvens, in another succinate-forming bacterium Succinimonas amylolytica DSM 2873  and one gene (yscN) was found in the sulphate-reducing bacterium Desulfotomaculum ruminis DSM 2154 . None of the ruminal strains of Sharpea, or the human gut organism Akkermansia mucinophila (the only genome available from a member of the phylum Verrucomicrobia isolated from a gut environment), contained T3SS gene homologues, therefore it is likely that the T3SS assigned to these genera by MEGAN belong to unidentified ruminal species.
From examination of their draft genome sequences, strains of the rumen bacterial species S. dextrinosolvens appear to encode complete T3SS, and along with Succinimonas amylolytica are the only two characterised rumen bacteria that have this type of secretory system. S. dextrinosolvens is a Gram negative, curved rod-shaped bacterium that degrades starch, and produces succinate, acetate, formate and sometimes lactate as its main end products of fermentation . It is known to be part of the epimural (adherent to the rumen wall) community in the bovine [45, 46] and is typically enriched when the diet contains high levels of starch [47, 48]. Succinivibrio are considered to be a small, but consistent part of the normal rumen microbiome  and is not known to be pathogenic to ruminants, although it has been reported as the infective agent in 2 cases of human disease [50, 51]. It seems unlikely that Succinivibrio spp. use T3SS for mounting infections, and it is more plausible that it mediates non-lethal relationships, either with other rumen organisms, or possibly the animal host itself. In pathogenic bacteria, effector proteins secreted by the T3SS act to modulate host cell functions to help avoid immune detection and disable protective functions such as macrophages . An amino acid sequence-based analysis of the protein-coding genes of S. dextrinosolvens strains H5, 22B and ACV-10 identified a large number of proteins containing N-terminal motifs indicative of involvement in T3SSs (Table 4), including chaperone binding domains and eukaryotic-like signal sequences specific for the endoplasmic reticulum or mitochondria (Additional file 2). The annotated functions of these proteins were quite diverse but many of them have putative roles in transcriptional regulation or transport functions, indicating they may mediate changes in gene expression and transport of molecules in the target organisms.
The specific contribution of the T3SS genes detected in the microbiomes of the HMY animals to the HMY phenotype remains unclear. The HMY phenotype in sheep is thought to be related to rumen size, particle retention time and turnover rate, such that HMY animals have a larger rumen which retains feed particles for longer and therefore have a slower turnover rate. These conditions are thought to lead to low hydrogen partial pressure, but rapid hydrogen production, which results in elevated expression of genes encoding the hydrogenotrophic methanogenesis pathway . The assignment of several of the T3SS genes to Succinivibrio suggests that this bacterial genus may be important in the HMY microbiome. However, this does not appear to be the case as the 16S rRNA genes retrieved from the metagenome sequences or from amplicon sequencing did not show Succinivibrio spp. as being significantly enriched in either the LMY (0.065% relative abundance) or HMY animals (0.033%). Investigations of the microbiomes of cattle with differing methane yield, have shown that OTUs assigned at a higher family level to the Succinivibrionaceae were more abundant in the microbiomes of LMY animals , while another study of cattle has reported OTUs corresponding to Succinivibrionaceae were reduced in feed restricted animals compared to ad libitum fed animals . The majority of T3SS genes which make up the main component of the correlation with methane yield to appear to come from unidentified rumen bacteria, but without knowledge of their metabolisms or physiologies, it is not possible to predict how their T3SS may influence the methane yield phenotype.
The co-culture experiments showed that S. dextrinosolvens had either stimulatory, inhibitory or neutral effects, depending on the type of methanogens in the co-culture. Succinivibrio spp. have been observed on several occasions in enrichment cultures with methanogens of the order Methanomassiliicoccales, and have proven difficult to remove from these enrichments to allow purification of the methanogens . This suggests that Succinivibrio spp. form close relationships with Methanomassiliicoccales, and helps explain the stimulation of ISO4-H5 by S. dextrinosolvens H5 in co-cultures. In a global census of microbes from ruminant animals, Henderson et al.  reported a positive association between succinate-producing Succinivibrionaceae, and methanogens belonging to the family Methanomassiliicoccaceae, particularly the subgroups, Candidatus Methanomethylophilus alvus and Methanomassiliicoccales group 11 sp. The mechanism of stimulation is currently not known, but it is possible that Succinivibrio spp. provide substrates or growth factors that stimulate the growth and methane formation in the methanogenic archaeon ISO4-H5. Some Succinivibrio spp. are able to degrade pectin  and release methanol  which is a known substrate for methane formation by the Methanomassiliicoccaceae .
The inhibition of methane formation during S. dextrinosolvens co-culture with Mbb. millerae SM9 is interesting as the Succinivibrionaceae bacterium WG-1 isolate from the Wallaby gut has been implicated in lower methane emissions from starch-containing diets, presumably via close coupling of redox reactions which led to less methane being formed . The global census of rumen microbes also reported a negative association between Succinivibrionaceae and the Methanobrevibacter gottschalkii clade . In our co-cultures, hydrogen was provided in excess, so the weak inhibition of Methanobrevibacter sp. SM9 by S. dextrinosolvens may be via a direct effect, rather than via competition for hydrogen. This is possibly via the production of formate by S. dextrinosolvens in the co-culture, which cannot be used for methane formation by Mbb. millerae  and which can slow growth when added to Mbb. millerae SM9 cultures growing on H2 + CO2 (Peter Janssen, pers. comm.). Investigating the exact mechanisms by which Succinivibrio spp. influences methanogen growth and methane formation, and the potential involvement of their T3SSs, are the subject of on-going research.
Bacterial T3SS genes and transcripts, were found to be positively correlated with methane yield in sheep. Most of these genes could not be assigned phylogenetically, but several were affiliated with the genus Succinivibrio, and genes encoding complete T3SSs were found in the genome sequences of S. dextrinosolvens strains H5, 22B and ACV-10, and another rumen succinate-producing rumen bacterium, Succinimonas amylolytica. This is the first report of T3SS genes being associated with methane emissions in ruminants, and identifies these secretions systems as potential new targets for methane mitigation research. S. dextrinosolvens H5 was shown to have direct growth-enhancing effects on a member of the Methanomassiliicoccales, and an inhibitory effect on a member of the Mbb. gottschalkii clade in co-culture experiments, which point towards bacteria-methanogen interactions being important modulators of methane production in ruminant animals.
high methane yield
low methane yield
dry matter intake
intermediate methane yield
gene set enrichment analysis
sparse partial least squares
Wilcoxon Rank Sum
type III secretion system
Kyoto encyclopedia of genes and genomes
mean squared prediction error
normalised enrichment score
integrated microbial genomes and microbiomes (expert review)
false discovery rates
reads per million
- NOM P-val:
nominal P value
IPCC. In climate change 2007: mitigation of climate change: contribution of working group III to the fourth assessment report of the intergovernmental panel on climate change. In: Metz ORD B, Bosch PR, Dave R, Meyer LA, editors. Cambridge, United Kingdom and New York, NY, USA; 2007.
Yusuf RO, Noor ZZ, Abba AH, Hassan MAA, Din MFM. Renewable and sustainable energy reviews. Renew Sustain Energy Rev. 2012;16:5059–70.
Zinder SH. Physiological ecology of methanogens. In: Ferry JG, editor. Methanogenesis: ecology, physiology, biochemistry and genetics. London: Chapman and Hall; 1993. p. 128–206.
Hook SE, Wright A-DG, McBride BW. Methanogens: methane producers of the rumen and mitigation strategies. Archaea. 2010;2010:1–11.
Lassey KR, Ulyatt MJ, Martin RJ, Walker CF, Shelton ID. Methane emissions measured directly from grazing livestock in New Zealand. Atmos Environ. 1997;31:2905–14.
Pinares-Patiño CS, Ebrahimi SH, McEwan JC, Dodds KG, Clark H, Luo D. Is rumen retention time implicated in sheep differences in methane emission? Proc New Zeal Soc Anim Prod. 2011;71:219–22.
Pinares-Patino CS, McEwan JC, Dodds KG, Cardenas EA, Hegarty RS, Koolaard JP, et al. Repeatability of methane emissions from sheep. Anim Feed Sci Tech. 2011;166–67:210–8.
Pinares-Patiño CS, Hickey SM, Young EA, Dodds KG, MacLean S, Molano G, et al. Heritability estimates of methane emissions from sheep. Animal. 2013;7:316–21.
Bain WE, Bezuidenhout L, Jopson NB, Pinares-Patino C, McEwan JC. Rumen differences between sheep identified as being low or high methane emitters. Proc 10th World Cong Genet Appl Livest Prod. 2014. doi:10.13140/2.1.3702.7206.
Goopy JP, Donaldson A, Hegarty R, Vercoe PE, Haynes F, Barnett MVHO. Low-methane yield sheep have smaller rumens and shorter rumen retention time. Br J Nutr. 2014;111:578–85.
Herd RM, Arthur PF, Donoghue KA, Bird SH, Bird-Gardiner T, Hegarty RS. Measures of methane production and their phenotypic relationships with dry matter intake, growth, and body composition traits in beef cattle. J Anim Sci. 2014;92:5267–74.
Robinson DL, Goopy JP, Hegarty RS, Oddy VH, Thompson AN, Toovey AF, et al. Genetic and environmental variation in methane emissions of sheep at pasture. J Anim Sci. 2014;92:4349–63.
Robinson DL, Goopy JP, Donaldson AJ, Woodgate RT, Oddy VH, Hegarty RS. Sire and liveweight affect feed intake and methane emissions of sheep confined in respiration chambers. Animal. 2014;8:1935–44. doi:10.1017/S1751731114001773.
Pickering NK, Chagunda MG, Banos G, Mrode R, McEwan JC, Wall E. Genetic parameters for predicted methane production and laser methane detector measurements. J Anim Sci. 2015;93:11–20. doi:10.2527/jas.2014-8302.
Orskov ER, Ojwang I, Reid GW. A study on consistency of differences between cows in rumen outflow rate of fibrous particles and other substrates and consequences for digestibility and intake of roughages. Anim Prod. 1988;47:45–51.
Smuts M, Meissner HH, Cronje PB. Retention time of digesta in the rumen—its repeatability and relationship with wool production of merino rams. J Anim Sci. 1995;73:206–10.
Benchaar C, Pomar C, Chiquette J. Evaluation of dietary strategies to reduce methane production in ruminants: a modelling approach. Can J Anim Sci. 2001;81:563–74.
Pinares-Patino CS, Ulyatt MJ, Lassey KR, Barry TN, Holmes CW. Rumen function and digestion parameters associated with differences between sheep in methane emissions when fed chaffed lucerne hay. J Agric Sci. 2003;140:205–14.
Kittelmann S, Seedorf H, Walters WA, Clemente JC, Knight R, Gordon JI, et al. Simultaneous amplicon sequencing to explore co-occurrence patterns of bacterial, archaeal and eukaryotic microorganisms in rumen microbial communities. PLoS ONE. 2013;8:e47879.
Shi W, Moon CD, Leahy SC, Kang D, Froula J, Kittelmann S, et al. Methane yield phenotypes linked to differential gene expression in the sheep rumen microbiome. Gen Res. 2014;24:1517–25.
Kamke J, Kittelmann S, Soni P, Li Y, Tavendale M, Ganesh S, et al. Rumen metagenome and metatranscriptome analyses of low methane yield sheep reveals a Sharpea-enriched microbiome characterised by lactic acid formation and utilisation. Microbiome. 2016;4:56. doi:10.1186/s40168-016-0201-2.
Yu Z, Morrison M. Improved extraction of PCR-quality community DNA from digesta and fecal samples. BioTechn. 2004;36:808–12.
Rosewarne CP, Pope PB, Denman SE, McSweeney CS, O’Cuiv P, Morrison M. High-yield and phylogenetically robust methods of DNA recovery for analysis of microbial biofilms adherent to plant biomass in the herbivore gut. Microb Ecol. 2011;61:448–54.
Caporaso JC, Kuczynski J, Stombaugh J, Bittinger K, Bushman FD, Costello EK, et al. QIIME allows analysis of high-throughput community sequencing data. Nat Methods. 2010. doi:10.1038/nmeth.f.303.
McDonald D, Price MN, Goodrich J, Nawrocki EP, DeSantis TZ, Probst A, et al. An improved Greengenes taxonomy with explicit ranks for ecological and evolutionary analyses of bacteria and archaea. ISME J. 2012;6:610–8.
Magoc T, Salzberg SL. FLASH: fast length adjustment of short reads to improve genome assemblies. Bioinformatics. 2011;27:2957–63.
Kanehisa M, Goto S. KEGG: Kyoto encyclopedia of genes and genomes. Nucl Acids Res. 2000;28:27–30.
Edgar RC. Search and clustering orders of magnitude faster than BLAST. Bioinformatics. 2010;26:2460–1.
Subramanian A, Kuehn H, Gould J, Tamayo P, Mesirov JP. GSEA-P: a desktop application for gene set enrichment analysis. Bioinformatics. 2007;23:3251–3.
R Core Team. R: a language and environment for statistical computing. Vienna: R Foundation for Statistical Computing; 2013.
Chung D, Chun H, Keleş S. Spls: sparse partial least squares (SPLS) regression and classification. R package; 2012.
Markowitz VM, Chen IM, Chu K, Szeto E, Palaniappan K, Pillay M, et al. IMG/M 4 version of the integrated metagenome comparative analysis system. Nucl Acids Res. 2014;42:D568–73. doi:10.1093/nar/gkt919.
Zhu W, Lomsadze A, Borodovsky M. Ab initio gene identification in metagenomic sequences. Nucl Acids Res. 2010;38:e132.
Altschul S, Gish W, Miller W, Myers EW, Lipman DJ. Basic local alignment search tool. J Mol Biol. 1990;5(215):403–10.
Huson DH, Mitra S, Ruscheweyh HJ, Weber N, Schuster SC. Integrative analysis of environmental sequences using MEGAN4. Genome Res. 2011;21:1552–60. doi:10.1101/gr.120618.111.
Eichinger V, Nussbaumer T, Platzer A, Jehl M-A, Arnold R, Rattei T. EffectiveDB—updates and novel features for a better annotation of bacterial secreted proteins and Type III, IV, VI secretion systems. Nucl Acids Res. 2016;4:D669–74. doi:10.1093/nar/gkv1269.
Kenters N, Henderson G, Jeyanathan J, Kittelmann S, Janssen PH. Isolation of previously uncultured rumen bacteria by dilution to extinction using a new liquid culture medium. J Microbiol Meth. 2011;84:52–60.
Joblin KN, Naylor GE, Williams AG. Effect of Methanobrevibacter smithii on xylanolytic activity of anaerobic ruminal fungi. Appl Environ Microbiol. 1990;56:2287–95.
Galán JE, Lara-Tejero M, Marlovits TC, Wagner S. Bacterial type III secretion systems: specialized nanomachines for protein delivery into target cells. Annu Rev Microbiol. 2014;68:415–38. doi:10.1146/annurev-micro-092412-155725.
Barker J, Brown MRW. Trojan horses of the microbial world: protozoa and the survival of bacterial pathogens in the environment. Microbiology. 1994;140:1253–9.
Fauvart M, Michiels J. Rhizobial secreted proteins as determinants of host specificity in the rhizobium legume symbiosis. FEMS Microbiol Lett. 2008;285:1–9.
Bryant MP, Small N. Characteristics of two new genera of anaerobic curved rods isolated from the rumen of cattle. J Bacteriol. 1956;72:22–6.
Bryant MP, Small N, Bouma C, Chu H. Bacteroides ruminicola n. sp. and Succinimonas amylolytica; the new genus and species; species of succinic acid-producing anaerobic bacteria of the bovine rumen. J Bacteriol. 1958;76:15–23.
Campbell LL, Postgate JR. Classification of the spore-forming sulfate-reducing bacteria. Bacteriol Rev. 1965;29:359–63.
Petri RM, Schwaiger T, Penner GB, Beauchemin KA, Forster RJ, McKinnon JJ, et al. Changes in the rumen epimural bacterial diversity of beef cattle as affected by diet and induced ruminal acidosis. Appl Environ Microbiol. 2013;79:3744–55.
Mead LJ, Jones GA. Isolation and presumptive identification of adherent epithelial bacteria (“epimural” bacteria) from the ovine rumen wall. Appl Environ Microbiol. 1981;41:1020–8.
Bryant MP. Bacterial species of the rumen. Bacteriol Rev. 1959;23:125–53.
Bryant MP. Normal flora-rumen bacteria. Am J Clin Nutr. 1970;23:1440–50.
Henderson G, Cox F, Ganesh S, Jonker A, Young W, Collaborators Global Rumen Census, et al. Rumen microbial community composition varies with diet and host, but a core microbiome is found across a wide geographical range. Sci Rep. 2015;5:14567. doi:10.1038/srep14567.
Southern PM. Bacteremia due to Succinivibrio dextrinosolvens. Report of a case. Am J Clin Pathol. 1975;64:540–3.
Porschen RK, Chan P. Anaerobic vibrio-like organisms cultured from blood: Desulfovibrio desulfuricans and Succinivibrio species. J Clin Microbiol. 1977;5:444–7.
Dean P. Functional domains and motifs of bacterial type III effector proteins and their roles in infection. FEMS Microbiol Rev. 2011;35:1100–25.
Wallace RJ, Rooke JA, McKain N, Duthie C-A, Hyslop JJ, Ross DW, et al. The rumen microbial metagenome associated with high methane production in cattle. BMC Genom. 2015;16:1–14.
McCabe MS, Cormican P, Keogh K, O’Connor A, O’Hara E, Palladino RA, et al. Illumina MiSeq phylogenetic amplicon sequencing shows a large reduction of an uncharacterised Succinivibrionaceae and an increase of the Methanobrevibacter gottschalkii clade in feed restricted cattle. PLoS ONE. 2015;10:e0133234. doi:10.1371/journal.pone.0133234.
Jeyanathan J. Investigation of rumen methanogens in New Zealand livestock. PhD Thesis. New Zealand: Massey University; 2010.
Pol A, Demeyer DI. Fermentation of methanol in the sheep rumen. Appl Environ Microbiol. 1988;54:832–4.
Paul K, Nonoh JO, Mikulski L, Brune A. “Methanoplasmatales”, Thermoplasmatales-related archaea in termite guts and other environments, are the seventh order of methanogens. Appl Environ Microbiol. 2012;78:8245–53.
Pope PB, Denman SE, Jones M, Tringe SG, Barry K, Malfatti SA, et al. Adaptation to herbivory by the Tammar wallaby includes bacterial and glycoside hydrolase profiles different from other herbivores. Proc Natl Acad Sci USA. 2010;107:14793–8.
Miller TL, Lin C. Description of Methanobrevibacter gottschalkii sp. nov., Methanobrevibacter thaueri sp. nov., Methanobrevibacter woesei sp. nov. and Methanobrevibacter wolinii sp. nov. Int J Syst Evol Microbiol. 2002;52:819–22.
The following authors contributed to the study design: GA, ER. Data collection and experimental procedures were conducted by: JK, GA, BK, SL, PS, and YL. Data analysis and interpretation: JK, GA, BK, SL, SG, WS, JF. The manuscript was prepared by: JK and GA. All authors read and approved the final manuscript.
We thank Christina Moon for critiquing the manuscript.
The authors declare that they have no competing interests.
Animal ethics approval
The collections of rumen contents from sheep were carried out under the approval of the AgResearch Ltd Grasslands Animal Ethics Committee (Approval 13606).
Availability of data and materials
The metagenome and metatranscriptome datasets used in this study are accessible at the National Centre for Biotechnology Information Sequence Read Archive (SRA; http://www.ncbi.nlm.nih.gov/sra) accession number SRA075938, bioproject number PRJNA202380, and additional 16S rRNA gene amplicon sequence data was submitted under the same bioproject number under SRA experiment Accession Numbers: SRX1079958–SRX1079985.
Consent to participate
This work was supported by the New Zealand Fund for Global Partnerships in Livestock Emissions Research, supporting the objectives of the Livestock Research Group of the Global Research Alliance on Agricultural Greenhouse Gases. The US Department of Energy Joint Genome Institute contribution was supported by the Office of Science of the US Department of Energy under Contract No. DE-AC02-05CH11231.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.