Skip to main content
  • Research article
  • Open access
  • Published:

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

Abstract

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

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.

Background

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 [1], 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 [8], and the animal characteristics controlling feed particle retention time [15,16,17,18] and rumen volume [10] 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 [19] and methanogen gene expression [20] 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 [21]. 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.

Methods

Study aim and design

This report is part of a larger study [20], 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 [8] 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 [20]. 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 [20]. Briefly, DNA was extracted using the “Repeated Bead Beating and Column (RBB + C) purification” method [22]. For large paired-end insert libraries, high-molecular-weight DNA was extracted using the method of Rosewarne et al. [23]. RNA extraction was via a hot lysis-acid phenol extraction method [20]. 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 [24] and phylogenetically assigned using the Greengenes database (version gg_13_5) [25]. Metagenomic and metatranscriptome sequencing datasets were quality and artifact filtered and paired end reads merged using FLASH [26] and used for read count based comparative analysis by screening against the KEGG database [27], (version 58.1) using USEARCH 6.0 [28] 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.

Data analyses

Statistical analyses were conducted using reads per million (RPM)-normalised read count matrices as previously described [21]. Gene set enrichment analysis (GSEA, [29]) 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 [29] 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 [30] was used to estimate correlations between specific KEGG genes and methane yield, using the spls package [31]. 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) [20] 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 [20]. 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 [32] 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 [33].

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 [34]. The output was analysed using MEGAN5 [35] for phylogenetic assignment using default settings. Matches were confirmed by homology searches of relevant genes in the IMG/M database [32]. Protein coding sequences from the draft genomes of S. dextrinosolvens strains H5, 22B and ACV-10 were analysed using EffectiveDB [36], 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.

Co-culture experiments

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 [37] with hydrogen + carbon dioxide (for the hydrogenotrophic methanogens, 1 atm over-pressure of a 80:20, hydrogen:carbon dioxide) or BY medium [38] 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

Not applicable.

Results

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).

Table 1 T3SS genes significantly enriched in metagenome and metatranscriptome datasets from HMY animals
Table 2 T3SS subunits identified as correlation predictors in sPLS analysis
Fig. 1
figure 1

Multiple regression plot of 10 bacterial T3SS subunit genes. The multiple regression plot is based on partial least squares prediction of correlation coefficient over all 20 samples. Low methane yield animals are shown in green, intermediate methane yield animals in blue and high methane yield animals in red. Adjusted multiple regression coefficient R2 = 0.76, P = 3.48 × 10−7

Table 3 Gene set enrichment analysis of the metagenome dataset
Fig. 2
figure 2

Z-score based hierarchical clustering of bacterial T3SS genes. Hierarchical clustering based on metagenome (a) and metatranscriptome (b) data

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:3.6.3.14]) but in no other rumen bacteria.

Fig. 3
figure 3

Phylogenetic affiliations of T3SS genes. Sunburst plot showing phylogenetic affiliations of contigs containing T3SS genes within the domain 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).

Fig. 4
figure 4

T3SS gene arrangement in genomes of Succinivibrio dextrinosolvens strains. Sct genes are shown in yellow with their associated gene subunit letter. Scaffold boundaries are shown by square brackets. The remaining genes are identified by COG categories

Table 4 Effective T3 predictions of chaperone binding, secreted, N-terminal signal and T4SS proteins with protein coding sequences of Succinivibrio dextrinosolvens

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).

Fig. 5
figure 5

Co-cultures of rumen methanogens with Succinivibrio dextrinosolvens H5. Methane production and hydrogen consumption during co-culture experiments of Succinivibrio dextrinosolvens H5 with the rumen methanogens Methanobrevibacter ruminantium M1 (a), Methanobrevibacter sp. SM9 (b), Methanosphaera sp. ISO3-F5(c) and methanogenic archaeon ISO4-H5 (d). Red arrows indicate the time point of inoculation with S. dextrinosolvens H5

Table 5 Concentrations (mM) of volatile fatty acids in co-cultures

Discussion

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 [8]. 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 [19] but large differences in the expression of genes encoding the methanogenesis pathway [20]. 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 [39] 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 [40]. 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 [41].

Fig. 6
figure 6

Schematic overview of T3SS subunits based on KEGG genes. The T3SS subunits are labelled using the Yersinia gene annotation scheme. The colour shading of each subunit shows the level of support for the related genes (DNA) or transcripts (RNA) being significantly (P ≤ 0.05) more abundant in HMY animals based on the WRS test and sPLS

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 [42]. 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 [43] and one gene (yscN) was found in the sulphate-reducing bacterium Desulfotomaculum ruminis DSM 2154 [44]. 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 [42]. 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 [49] 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 [52]. 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 [20]. 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 [53], while another study of cattle has reported OTUs corresponding to Succinivibrionaceae were reduced in feed restricted animals compared to ad libitum fed animals [54]. 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 [55]. 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. [49] 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 [42] and release methanol [56] which is a known substrate for methane formation by the Methanomassiliicoccaceae [57].

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 [58]. The global census of rumen microbes also reported a negative association between Succinivibrionaceae and the Methanobrevibacter gottschalkii clade [49]. 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 [59] 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.

Conclusions

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.

Abbreviations

HMY:

high methane yield

LMY:

low methane yield

DMI:

dry matter intake

IMY:

intermediate methane yield

GSE:

gene set enrichment analysis

sPLS:

sparse partial least squares

WRS:

Wilcoxon Rank Sum

T3SS:

type III secretion system

KEGG:

Kyoto encyclopedia of genes and genomes

MSPE:

mean squared prediction error

NES:

normalised enrichment score

IMG/Mer:

integrated microbial genomes and microbiomes (expert review)

FDR:

false discovery rates

BH:

Benjamini–Hochberg procedure

MEGAN:

metagenome analyser

MS:

membrane supramembrane

RPM:

reads per million

NOM P-val:

nominal P value

FWER:

familywise-error rate

References

  1. 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.

  2. Yusuf RO, Noor ZZ, Abba AH, Hassan MAA, Din MFM. Renewable and sustainable energy reviews. Renew Sustain Energy Rev. 2012;16:5059–70.

    Article  CAS  Google Scholar 

  3. Zinder SH. Physiological ecology of methanogens. In: Ferry JG, editor. Methanogenesis: ecology, physiology, biochemistry and genetics. London: Chapman and Hall; 1993. p. 128–206.

    Chapter  Google Scholar 

  4. Hook SE, Wright A-DG, McBride BW. Methanogens: methane producers of the rumen and mitigation strategies. Archaea. 2010;2010:1–11.

    Article  Google Scholar 

  5. 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.

    Article  CAS  Google Scholar 

  6. 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.

    Google Scholar 

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

    Article  Google Scholar 

  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.

    Article  PubMed  PubMed Central  Google Scholar 

  9. 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.

    Google Scholar 

  10. 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.

    Article  CAS  PubMed  Google Scholar 

  11. 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.

    Article  CAS  PubMed  Google Scholar 

  12. 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.

    Article  CAS  PubMed  Google Scholar 

  13. 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.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  14. 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.

    Article  CAS  PubMed  Google Scholar 

  15. 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.

    Article  Google Scholar 

  16. 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.

    Article  CAS  PubMed  Google Scholar 

  17. 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.

    Article  Google Scholar 

  18. 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.

    Article  CAS  Google Scholar 

  19. 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.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  20. 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.

    Article  CAS  Google Scholar 

  21. 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.

    Article  PubMed  PubMed Central  Google Scholar 

  22. Yu Z, Morrison M. Improved extraction of PCR-quality community DNA from digesta and fecal samples. BioTechn. 2004;36:808–12.

    CAS  Google Scholar 

  23. 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.

    Article  CAS  PubMed  Google Scholar 

  24. 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.

    PubMed  PubMed Central  Google Scholar 

  25. 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.

    Article  CAS  PubMed  Google Scholar 

  26. Magoc T, Salzberg SL. FLASH: fast length adjustment of short reads to improve genome assemblies. Bioinformatics. 2011;27:2957–63.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  27. Kanehisa M, Goto S. KEGG: Kyoto encyclopedia of genes and genomes. Nucl Acids Res. 2000;28:27–30.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  28. Edgar RC. Search and clustering orders of magnitude faster than BLAST. Bioinformatics. 2010;26:2460–1.

    Article  CAS  PubMed  Google Scholar 

  29. 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.

    Article  CAS  PubMed  Google Scholar 

  30. R Core Team. R: a language and environment for statistical computing. Vienna: R Foundation for Statistical Computing; 2013.

    Google Scholar 

  31. Chung D, Chun H, Keleş S. Spls: sparse partial least squares (SPLS) regression and classification. R package; 2012.

  32. 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.

    Article  CAS  PubMed  Google Scholar 

  33. Zhu W, Lomsadze A, Borodovsky M. Ab initio gene identification in metagenomic sequences. Nucl Acids Res. 2010;38:e132.

    Article  PubMed  PubMed Central  Google Scholar 

  34. Altschul S, Gish W, Miller W, Myers EW, Lipman DJ. Basic local alignment search tool. J Mol Biol. 1990;5(215):403–10.

    Article  Google Scholar 

  35. 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.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  36. 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.

    Article  Google Scholar 

  37. 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.

    Article  Google Scholar 

  38. Joblin KN, Naylor GE, Williams AG. Effect of Methanobrevibacter smithii on xylanolytic activity of anaerobic ruminal fungi. Appl Environ Microbiol. 1990;56:2287–95.

    CAS  PubMed  PubMed Central  Google Scholar 

  39. 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.

    Article  PubMed  PubMed Central  Google Scholar 

  40. 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.

    Article  CAS  PubMed  Google Scholar 

  41. Fauvart M, Michiels J. Rhizobial secreted proteins as determinants of host specificity in the rhizobium legume symbiosis. FEMS Microbiol Lett. 2008;285:1–9.

    Article  CAS  PubMed  Google Scholar 

  42. 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.

    CAS  PubMed  PubMed Central  Google Scholar 

  43. 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.

    CAS  PubMed  PubMed Central  Google Scholar 

  44. Campbell LL, Postgate JR. Classification of the spore-forming sulfate-reducing bacteria. Bacteriol Rev. 1965;29:359–63.

    CAS  PubMed  PubMed Central  Google Scholar 

  45. 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.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  46. 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.

    CAS  PubMed  PubMed Central  Google Scholar 

  47. Bryant MP. Bacterial species of the rumen. Bacteriol Rev. 1959;23:125–53.

    CAS  PubMed  PubMed Central  Google Scholar 

  48. Bryant MP. Normal flora-rumen bacteria. Am J Clin Nutr. 1970;23:1440–50.

    CAS  PubMed  Google Scholar 

  49. 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.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  50. Southern PM. Bacteremia due to Succinivibrio dextrinosolvens. Report of a case. Am J Clin Pathol. 1975;64:540–3.

    Article  PubMed  Google Scholar 

  51. Porschen RK, Chan P. Anaerobic vibrio-like organisms cultured from blood: Desulfovibrio desulfuricans and Succinivibrio species. J Clin Microbiol. 1977;5:444–7.

    CAS  PubMed  PubMed Central  Google Scholar 

  52. Dean P. Functional domains and motifs of bacterial type III effector proteins and their roles in infection. FEMS Microbiol Rev. 2011;35:1100–25.

    Article  CAS  PubMed  Google Scholar 

  53. 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.

    Article  Google Scholar 

  54. 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.

    Article  PubMed  PubMed Central  Google Scholar 

  55. Jeyanathan J. Investigation of rumen methanogens in New Zealand livestock. PhD Thesis. New Zealand: Massey University; 2010.

  56. Pol A, Demeyer DI. Fermentation of methanol in the sheep rumen. Appl Environ Microbiol. 1988;54:832–4.

    CAS  PubMed  PubMed Central  Google Scholar 

  57. 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.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  58. 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.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  59. 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.

    CAS  PubMed  Google Scholar 

Download references

Authors’ contributions

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.

Acknowledgements

We thank Christina Moon for critiquing the manuscript.

Competing interests

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

Not applicable.

Funding

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.

Publisher’s Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Author information

Authors and Affiliations

Authors

Corresponding author

Correspondence to Graeme T. Attwood.

Additional files

Rights and permissions

Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Kamke, J., Soni, P., Li, Y. et al. Gene and transcript abundances of bacterial type III secretion systems from the rumen microbiome are correlated with methane yield in sheep. BMC Res Notes 10, 367 (2017). https://doi.org/10.1186/s13104-017-2671-0

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s13104-017-2671-0

Keywords