The community structure and microbial linkage of rumen protozoa and methanogens in response to the addition of tea seed saponins in the diet of beef cattle

Cui Tan,Carlos A.Ramírez-Restrepo,Ali Mujtaba Shah,4,Rui Hu,5,Matt Bell,Zhisheng Wang,5*and Chris McSweeney

Abstract Background:This study investigated changes in rumen protozoal and methanogenic communities,along with the correlations among microbial taxa and methane(CH4)production of six Belmont Red Composite beef steers fed tea seed saponins(TSS).Animals were fed in three consecutive feeding periods,a high-grain basal diet for 14 d(BD period)then a period of progressive addition of TSS to the basal diet up to 30 g/d for 20 d(TSS period),followed by the basal diet for 13 d without TSS(BDP post-control period).Results:The study found that TSS supplementation decreased the amount of the protozoal genus Entodinium and increased Polyplastron and Eudiplodinium genera.During BDP period,the protozoa community of steers did not return to the protozoal profiles observed in BD period,with higher proportions of Metadinium and Eudiplodinium and lower Isotricha.The addition of TSS was found to change the structure of methanogen community at the subgenus level by decreasing the abundance of methanogens in the SGMT clade and increasing the abundance of methanogens in the RO clade.The correlation analysis indicated that the abundance of SGMT clade methanogens were positively correlated with Isotricha,and Isotricha genus and SGMT clade methanogens were positively correlated with CH4 production.While RO clade were positively correlated with the proportion of Metadinium genus,which was negatively correlated with CH4 emission.Conclusions:These results suggest that different genera of rumen protozoa ciliates appear to be selectively inhibited by TSS,and the change in methanogen community at the subgenus level may be due to the mutualistic relationships between methanogens and rumen ciliates.

Keywords:Beef steers,Methane,Methanogen,Protozoa,Tea seed saponin

Background

Specific ciliate protozoa colonize the rumen ecosystem of ruminant animals and contribute 40%～50% of the rumen microbial biomass,which is almost equal to that of rumen bacteria[1].Rumen protozoa play an important role in feed digestion and homeostasis of the rumen ecosystem due to the predator-prey relationship and the symbiotic relationship between rumen ciliates and prokaryotic cells[2].Studies using electron microscopy and sequencing revealed the presence of methanogens inside and outside of rumen ciliates[3-5].And the higher methanogens/bacteria ratio observed within the protozoa fractions suggested that anaerobic ciliates contain symbiotic methanogens,which strengthened the association between rumen ciliates and methanogens[6].The strong functional association between rumen protozoa and methanogens via interspecies hydrogen transfer contributes to overall methane(CH4)production(9%～37%)in the rumen as demonstrated by defaunation or elimination of certain protozoa[7].Thus,there is increasing evidence about the relationship between protozoa and methanogenic archaea,and their important role in rumen(CH4)formation[7-9].

Many rumen protozoa contain an organelle(hydrogenosome)which produces hydrogen that attracts methanogens as endosymbionts[10,11].Hydrogenosomes have been identified in species ofEpidinium,IsotrichaandDasytricha,but they were not found inEntodinium caudatumandDiploplastron affine[12].Tymensen et al.[4]reported thatMethanobrevibacterandMethanomicrobiumgenera are the major groups of methanogen in the free-living(FL)and protozoa-associated methanogen(PAM)fractions,but their population structure differed between FL and PAM pools.In addition,the populations of PAM were more variable than FL methanogens[5].It has been found that rumen methanogens closely related toM.smithiiwere primarily associated with the vestibuliferid protozoa(IsotrichaandDasytricha)[13].While

Polyplastron multivesiculatumharbored a large number of intracellular bacteria,no methanogen was observed[3].Methanogens also associate with protozoa as ectosymbionts which may be a more stable population compared with the internal dwelling methanogens[11].Due to the apparent variability in association between methanogens and different protozoal taxa,it is hypothesized that some protozoa species may have a greater effect on rumen methanogenesis and CH4production than others.

Some studies have shown that variation in CH4emissions are not necessarily due to the changes in methanogen cell density,but rather the community structure of the methanogens[14,15].Phytochemicals such as saponins,tannins and essential oils,have been used to manipulate rumen microbial populations and reduce CH4emissions from ruminants[16].Tea seed(Camellia sinensisL.)saponins(TSS)were reported to reduce methane production by an inhibitory effect on ciliate protozoa[17-19].Tannins,essential oils,and linseed oil are also inhibitory to protozoa,but species are impacted differently[20-22].Our previous study in beef cattle supplemented with TSS showed reduction in CH4emissions were associated with changes in methanogen population structure,while total protozoal numbers remained unchanged[23].It remains to be determined if there was an alteration of rumen ciliate composition in response to supplemental TSS,which further modified the methanogens composition and CH4emissions.The recent development of sequencing primers that can discriminate members of the rumen protozoal populations may provide an opportunity to characterize changes in the relative abundance of different protozoal taxa that previously could only be determined by microscopy[24].

Therefore to gain greater insight into the community structure of ruminal ciliate protozoa and their relationship to methanogen diversity and CH4production,we employed next-generation sequencing techniques to investigate these populations in the rumen of beef steers fed a high grain-based diet supplemented with TSS.

Methods

Animals,diet and experimental design

Care of animals and experimental procedures were approved by the CSIRO Animal Ethic Committee protocol No A12/2012.This research is part of a larger study that has been previously reported and see Ramirez-Restrepo et al.[23]for more details.Briefly,six Belmont Red Composite[Africander(African Sanga)×Brahman(Bos indicus)×Hereford-Shorthorn(3/4B.taurus)]rumencannulated steers[327±17.1 kg,initial body weight]were used in a repeated measures design.

After free grazing,the steers were individually allowed to adjust to a conventional-finishing feedlot basal diet(BD)fed ad libitum in two equal parts at 09:00 and 16:30 for 14 d(control BD period).This was followed by supplementing the BD in the morning feeding with 6,10,15,20,25 and 30 g/d of TSS light-yellow color powder dissolved in～250 mL cold water during 5,2,4,3,2 and 4 d,respectively(TSS period).After that,using only the BD,steers were fed for other 13 days(post-control period;BDP).The basal diet consisted of Rhodes grass(Chloris gayana;0.15)and a sorghum high-grain mixture(0.85;Coleman Stock Feeds Pty Ltd.,Charters Towers,Australia;Additional file 1:Table S1).The TSS per kg of commercial product(Zhejiang Oriental Tea Technology Co.,Ltd.,Changshan,China)included 570 g triterpenoid saponin,330 g crude fiber,54 g crude protein,50 g ash,and 1 g water-insoluble matter[25].Water was available ad libitum and steers were allowed to exercise daily in the cattle yard before the morning feeding to reduce stress,record experimental data and facilitate pen cleaning.

Methane emissions from all steers were measured in open-circuit respiration chambers for each period,and the results of the six cattle in our previous paper[23]showed that the average CH4production in the BDP period(14.2 g/kg DMI)was significantly lower than the CH4production in the BD period(18.0 g/kg DMI)and TSS period(18.7 g/kg DMI)respectively(P＜0.05).The lower CH4emissions in BDP period could be due to the reestablishment of a microbial community that had lower CH4producing potential than the controls.

Analysing protozoal and methanogenic communities

Rumen samples from each rumen-cannulated steer at the end of each feeding period were collected 2 h after the morning feeding.Raft and liquid components of the rumen content obtained from the dorsal and ventral sacks and were squeezed through two layers of cheese cloth.The fluid was immediately transferred to sterile containers placed on dry ice,and stored at-80°C for later DNA extraction.

Genomic DNA was extracted from 2 mL rumen fluid using the standard cetyl-trimethyl ammonium bromide(CTAB)method with some modifications to obtain higher concentrations and better purity of gDNA extracts from the rumen fluid[26].The procedure used for DNA extraction is described in detail in our previous study[23]using bead beating in CTAB lysis buffer followed by phenolchloroform separation.DNA extracts were dissolved in a 200μL buffer EB,and the quality and quantity of the extracted DNA were determined by UV spectrophotometric analysis,using a NanoDrop ND-1000 Spectrophotometer(Nyxor Biotech,Paris,France).

The V3-V4 region of the protozoan 18S rRNA gene was PCR amplified using PSSU316F-GIC758R[24]barcoded primers set,and the V3-V4 region of the methanogen 16S rRNA gene was amplified using 915F-1386R[27]barcoded primers set.The two-step PCR amplification procedure with the barcoded primers(Additional file 2:Table S2)were used to obtain the PCR amplicons for paired-end sequencing on an Illumina MiSeq platform.The first-step PCR amplification setup consisted of 10μL 5×Phusion HF buffer(including MgCl2),1μL dNTP(10 mmol/L),1μL inner primer(10μmol/L),20-50 ng of DNA template 1 U of Phusion High Fidelity DNA Polymerase(New England Biolabs,Beijing,China)and ddH2O up to 50μL.Amplification was initiated with denaturation for 2 min at 94°C,followed by 30 cycles of 94°C for 30 s,annealing at 55°C for protozoa or 56°C for methanogen amplifying for 1 min and 30 s and extension at 72°C for 30 s,with a final elongation for 2 min.After purification of PCR product,the cleaned PCR fragments were used in the second step PCR amplification.The 2nd-step PCR setup consisted of 8μL HF buffer,1μL dNTP,1μL outer primer,5μL cleaned PCR product,0.8 U of DNA Polymerase and ddH2O up to 40μL.PCR cycles were as follows:1 cycle at 94°C for 2 min;8 cycles at 94°C for 30 s,56°C for 1.5 min and 72°C for 30 s;and a final extension step at 72°C for 2 min.

The purified barcoded amplicons were pooled in equimolar quantity to generate the Illumina library and quantified on Qubit 3.0 Fluorometer(Invitrogen,Waltham,USA).Sequencing was conducted on the Illumina multiplex 2×300-bp paired-end MiSeq-platform in TINY gene biotechnology Co.Ltd.(Shanghai,China).The sequencing procedure was followed as per Kozich et al.[28].The raw sequencing data were submitted to the NCBI Sequence Read Archive(SRA),which are available from NCBI under BioProject accession number PRJNA574788(https://www.ncbi.nlm.nih.gov/bioproject/PRJNA574788)for both protozoal(n=18)and methanogenic(n=18)community datasets.

The paired-end reads were assembled by using Mothur software package[29],and then trimmed and qualityfiltered in QIIME using default settings[30].Briefly,Paired-end reads with ambiguous bases(n＞0),homologous regions(n≥8),mismatches in the primer(n＞0),nontarget areas and errors in barcodes(n＞0)were culled,as well as sequences less than 200 bp,longer than 580 bp.Subsequently,chimeric sequences were removed by UCHIME[31].And the quality-filtered sequences were clustered into operational taxonomic units(OTUs)for species classification with 96% similarity for ciliate protozoa[24].An 18S rRNA gene sequence reference,providing the alignment for candidate sequences of ciliate protozoa,contained 131 sequences(≥1,500 bp in length)belonging to the gastrointestinal tract trichostome(subclass Trichostomatia)ciliates and downloaded from NCBI.The sequences of methanogen dataset were assigned to OTUs at a 3% species-level cutoff[32],with the RDP(Ribosomal Database Project)database as the reference.For OTUs grouping,only one or two reads were removed before taxonomic assignment.A good depth of coverage,community diversity indicators,the Shannon-Wiener index,Simpson’s index and community richness estimates,such as the Chao 1 index,were also calculated using Mothur and self-compiled employing program in BioLinker.Phylogenetic tree was constructed to visualize the phylogeny of dominant methanogen OTUs,representing sequences across all samples,using the neighbor-joining method in the MEGA6[33].

Statistical analyses

All statistical analysis was performed using the R Package(version 3.5.1).Normality test of all the data was analyzed first.The normally distributed data,included OTU richness and diversity indices,were analyzed by using one-way analysis of variance(ANOVA)followed by the application of Tukey test for multiple comparison.The non-normally distributed data,taxonomic richness,were analyzed by using Kruskal-Wallis test to assess overall significant differences among the groups,followed by the post hoc Wilcoxon rank sum test to determine which paired groups differed from each other.Unless otherwise noted,the mean and pooled standard error of the mean(SEM)are reported.Spearman correlation analysis was used to assess microbial community interactions between protozoa and methanogen,as well as the relationship between microbial taxa abundance and CH4.Correlation matrix and graphics were generated using the online resource(https://www.omicshare.com/).P-values of＜0.05 for the overall significant differences andP-values of＜0.05 after FDR correction using the Benjamini-Hochberg procedure for the multiple comparison were considered significant.

Results

Effect on rumen protozoa communities

A total of 457,733 high-quality sequences mapped to protozoan 18S rRNA references generated 1,473 representative OTUs at 96%similarity,for an average of 870 OTUs per sample(range 593-1208).TSS supplementation reduced the number of observed protozoa OTUs and the species richness metrics,Chao 1 and ACE(FDR-P＜0.05;Table 1).This indicates a reduction in protozoal species richness for beef steers fed on a diet containing TSS.However,the alpha diversity indices of rumen protozoa,Shannon-Wiener and Simpson’s indices,were not influenced by TSS supplementation(P＞0.05,Table 1),so the reduction of species richness and similarity of alpha diversity indices with TSS indicated a shift in the evenness of protozoa taxa due to TSS.Good’s Coverage index showed good overall sampling and sequencing,ranging between 96.9%and 98.5%for rumen protozoa.

Across all samples,taxonomic classification based on the observed OTUs at the family and genus level resulted in the identification of 3 families and 10 genera of protozoa.The families Ophryoscolecidae and Isotrichidaeaccounted for 85.9% and 8.7% of the sequences,respectively(Additional file 3:Figure S1A).Entodinium(44.1%),Metadinium(10.5%),Eudiplodinium(5.9%)andPolyplastron(4.7%)were the dominant genera identified in the family Ophryoscolecidae.Another two genera in the family Isotrichidae,Isotricha(8.0%)andDasytricha(0.7%)were also identified in this study(Additional file 3:Figure S1B).And the community composition of rumen protozoa at the order,family and genus level among three periods are shown in Additional file 3:Figure S2.At the order and family level,less relative abundance of Entodiniomorphida and Ophryoscolecidae were observed in TSS supplementation steers(FDR-P＜0.05)compared with BDP period,while the abundance of Vestibuliferida order and Isotrichidae family in BDP period were much lower than that of BD and TSS periods(FDR-P＜0.05,Table 2).At the genus level,TSS supplementation significantly reduced the abundance ofEntodinium(FDR-P＜0.05),while significantly increased the abundance of thePolyplastrongenus(FDR-P＜0.05,Table 2).Meanwhile,the abundance ofEudiplodiniumin BDP period was higher than the BD period(FDR-P＜0.05).As for theMetadiniumgenus,its abundance in BDP period was higher than that in other two periods(FDR-P＜0.05),but the abundance ofIsotrichain BDP period was lower(FDR-P＜0.05,Table 2).

Table 1 Effect of tea seed saponins(TSS)supplementation on OTU richness and alpha diversity of ciliate protozoa at 96%similarity

Effect on rumen methanogenic communities

A total of 370,065 high-quality sequences were mapped to methanogen 16S rRNA references and generated 1,441 representative OTUs at 97% similarity,for anaverage of 459 OTUs per sample.The number of observed methanogen OTUs in the rumen of beef steers fed on a basal diet was not affected by TSS supplementation,while higher ACE index for rumen methanogens was observed in TSS supplementation steers compared with BDP period(FDR-P=0.08,Table 3).However,the alpha diversity indices of rumen methanogen,Shannon-Wiener and Simpson’s indices,were not affected by TSS supplementation(Table 3).

Table 2 Effect of tea seed saponins(TSS)on the protozoal communities(%of total sequences)

Across all samples,the methanogen OTUs annotated to 2 families and 5 genera(Table 4).Methanobacteriaceae andMethanobrevibacterspp.were the predominant family and genera of methanogen in the rumen of steers during the three periods,accounting for more than 98%of the sequences.The abundance of methanogens at the family and genus levels were not significantly affected by TSS supplementation,while the abundance ofMethanosphaerain BDP steers tended to be higher than in the BD period(FDR-P=0.08;Table 4).

Effect on the dominant methanogen OTUs

There were 15 dominant methanogen OTUs,averaging greater than 1% of total methanogen sequences across all samples,were used to construct the phylogenetic tree(Fig.1).The results showed that 86.7%(13 OTUs)of the methanogen OTUs were located in the SGMT clade of the genusMethanobrevibacter(M.smithii,M.gottschalkii,M.milleraeandM.thauerii),only 13.3%(2 OTUs)were located in the RO cluster of the genusMethanobrevibacter(M.ruminentiumandM.olleyae).

The two predominant OTUs,OTU 1 and OTU 2,both representing approximately 30% of all sequences,were located in the SGMT clade of theMethanobrevibactergenus.The relative abundances of the 15 dominant methanogen OTUs were compared among the three feeding periods(Table 5).Compared with BD period,TSS supplementation significantly reduced the relative abundance of OTU 4 and OTU 7(FDR-P＜0.05),whileincreased the proportion of OTU 5 and OTU 6(FDRP＜0.05).However,the abundance of SGMT clade and OTU 1 in BDP period was much lower than that in the BD period(FDR-P＜0.05),while the proportion of RO clade of theMethanobrevibactergenus and OTU 6 in BDP period was higher than BD period(FDR-P＜0.05).

Table 3 Effect of tea seed saponins(TSS)supplementation on OTU richness and alpha diversity of methanogen at 97%similarity

Table 4 Effect of tea seed saponins(TSS)on the methanogen communities(%of total sequences)

Within-and between-domain associations among protozoa and methanogen

Pearson’s correlations indicated significant associations among protozoa genera,and between protozoa genera and methanogen taxa.Within ciliate protozoa,positive correlations were observed between the abundance ofMetadiniumandEudiplodinium(r=0.55;P=0.02),betweenIsotrichaandDasytricha(r=0.55;P=0.02),as well as betweenPolyplastrongenus withOstracodiniumandOphryoscolex(r=0.54 and 0.51,respectively;P＜0.05),and strong negative correlation was noticed betweenEntodiniumandEudiplodinium(r=-0.66;P=0.01;Fig.2).In relation to protozoa-methanogen interactions,negative correlations were observed between theMethanobrevibactergenus withPolyplastronandOphryoscolex(r=-0.48 and-0.49,respectively;P＜0.05).While positive correlations occurred between theMethanosphaeragenus withEudiplodiniumandOstracodinium(r=0.62 and 0.56,respectively;P＜0.05),as well as between theMethanomethylovoransgenus withPolyplastronandOphryoscolex(r=0.48 and 0.54,respectively;P＜0.05;Fig.2).More significant correlations between the abundance of protozoal genera and the two clades of the genusMethanobrevibacter.SGMT clade was positively correlated withIsotricha(r=0.51;P=0.03),and negatively withEudiplodiniumandEudiplodinium(r=-0.49 and 0.51,respectively;P=0.03),the RO clade was positively correlated withMetadiniumandEudiplodinium(r=0.57 and 0.56,respectively;P＜0.05),and negatively withIsotricha(r=-0.51;P=0.01;Fig.2).

Fig.1 Phylogenetic relationships between the 15 dominant methanogenic OTUs of all steers and Methanobrevibacter spp.The sequences of 15 dominant methanogenic OTUs obtained from 16S rRNA illumina library prepared from the community DNA of rumen of Belmont Red Composite steers,and the reference sequences of the Methanobrevibacter genus obtained from GenBank.The phylogenetic tree was constructed using the Neighbor-Joining method with the Kimura two-parameter model

Table 5 Effect of tea seed saponins(TSS)on the abundance of dominant methanogen clade/OTUs

Correlations between microbial taxa with rumen CH4 emission

The yield of CH4per unit of DMI(dry matter intake)was positively correlated with the relative abundance ofIsotrichagenus(r=0.59,respectively;P＜0.05),and negatively correlated with the genusMetadinium(r=-0.50;P＜0.05).As for methanogens,only SGMT clade of the genusMethanobrevibacterwas positively correlated with CH4production(r=0.46;P=0.05),while no significant correlation between CH4emissions and methanogen abundance at the genus and family level were observed(Fig.2).

Discussion

Fig.2 Correlation matrix between protozoal taxa with methanogens and methane.Positive correlations are shown in red and negative correlations in blue.Color intensity represents r-values of correlation,*indicates P＜0.05.The scale of the colors is denoted as follows:the more positive the correlation(closer to 1),the darker the shade of red;the more negative the correlation(closer to-1),the darker the shade of blue.*indicates P＜0.05

Rumen microbial community is shaped by complex relationships between bacteria,ciliate protozoa,fungi and archaea.Previous studies have shown that TTS as a feed additive can potentially mitigate CH4emissions from ruminants through their inhibitory effect on protozoa,which subsequently may alter the abundance/diversity of methanogens[18,34].In the present study,TSS supplementation reduced the protozoal diversity metrics,but diversity still remained high throughout the study.We speculate that the high value diversity metrics is probably linked with the high-concentrate diet(concentrate to forage,75:25)used,which can promote the growth of rumen ciliates[35,36].Furthermore,the high 18S rRNA gene copies number for large entodiniomorphids protozoa can partially explain the high diversity of rumen protozoa in this study[37].In this study,the rumen ciliates were taxonomically assigned to 3 protozoan families and 10 genera.The genusEntodiniumspp.was the predominant genus in the most dominant family,Ophryoscolecidae,comprising up to 85%of the total ciliate population,which is in agreement with other studies with cattle and sheep[27,38].In the present study,the TSS supplementation in a grain-based diet modified the community composition of rumen ciliate protozoa.The relative abundance of the family Ophryoscolecidae and theEntodiniumgenus in particular were reduced by TSS,however,concurrently there was an increase in the relative proportion ofPolyplastronandEudiplodiniumgenera,which demonstrated that there was a significant selective inhibitory effect of TSS on the ciliate protozoa community at the family and genus level.This selective inhibitory effect of TSS on rumen protozoa could be due to the physical or chemical differences in the cilia,the pellicle or metabolic activity of rumen ciliates.The concentration and composition of ciliate protozoa in the rumen can be influenced by diet type,pH,supplements,concentrate and roughage level[35,39].The most consistent effect of tea saponin supplementation is a decrease in the abundance of rumen protozoa in vitro[40],as well as in the rumen of goats and sheep[17,18].

Most Ophryoscolecidae ciliates seem to be able to harbor both intracellular and extracellular methanogens,whereas only endosymbiont methanogens have been observed with the Isotrichidae ciliates[41].In addition,previous studies have reported thatEntodiniumspp.andIsotrichaspp.contain intracellular methanogens,whereasPolyplastron multivesiculatumharbored mainly intracellular bacteria[42].In the present study,after the withdrawal of TSS supplementation,the ciliate community structure was also dissimilar with a higher proportion ofMetadinium,PolyplastronandEudiplodiniumand lowerIsotricha,which may also affect other rumen microbial communities due to the predator-prey relationship and the symbiotic relationship between rumen ciliates and prokaryotes(rumen bacteria and methanogens)[43].The CH4production in the rumen is significantly influenced by the protozoa population due to inter-species transfer of hydrogen from the protozoa to methanogens[44].It has also been suggested that the number and activity of endosymbiotic population of methanogens are greater than the ectosymbionts attached to the surface of protozoa[45].In most cases,the variation in CH4emissions,caused by the addition of plant secondary compounds or microbial modifiers,could not be attributed to differences in methanogen cell density,but rather the composition of the methanogen community[46].These results suggest that the shift of the ciliate community composition may result in dissimilarity in methanogen community that could vary in their potential for CH4production.

In this study,the Methanobacteriaceae andMethanobrevibacterspp.were the predominant family and genus of methanogen in the rumen of steers throughout the trial,which is consistent with the previous findings[47,48].We suggest that this may be related to the high-concentrate diet used in our study.As Hook et al.[35]reported,approximately 99.3%of methanogen clones were assigned to the genusMethanobrevibacterin the rumen of cattles fed 65% concentrate diet.Furthermore,Hook et al.[35]and Franzolin[36]found that a high-concentrate diet promoted the abundance of protozoa.And other studies have found that theMethanobrevibacterwas the predominant genus associated with protozoa[4,5].Therefore,the high density of rumen protozoa may contribute to the dominance ofMethanobrevibacterin the rumen of catle fed the high-concentrate diet.St-Pierre et al.[46]observed thatMethanobrevibacterrelated methanogen 16S rRNA gene sequences are mainly distributed between two large clades:SGMT clade,closely related toM.smithii,M.gottschalkii,M.milleraeandM.thaueri,and RO clade,closely related toM.ruminantiumandM.olleyae.The phylogenetic tree of the 15 dominant methanogen OTUs in this study showed that 86.7% methanogen OTUs,accounted for more than 70% of the sequences,and were located in the SGMT clade of the genusMethanobrevibacterwith only 13.3% OTUs located in the RO clade.A reexamining report compared the sequence distribution between the SGMT clade and the RO clade and found that methanogens in the SGMT clade were distinctively more highly represented than the RO clade in impalas,wallabies,Holstein dairy cows and water buffaloes[44],which is consistent with the results of the present study.In the current study,the lower abundance of SGMT clade and higher proportion of RO clade were observed after the withdrawal of TSS supplementation,which was also associated with the lower proportion ofIsotrichaand higherMetadinium.What’s more,a positive correlation was observed between SGMT clade methanogens andIsotricha,as well as between the RO clade forMetadinium.Therefore,changes in the structure and composition of protozoan community should alter the structure of the methanogen community at the subgenus level.

The correlations between microbial taxa with CH4emissions in the present study showed that CH4production was positively correlated withIsotricha,but negatively correlated withMetadinium.This finding is supported by a previous study which reported that CH4production and methanogen population in sheep inoculated with Isotrichidae ciliates were similar to the total protozoa of inoculated sheep[49].The correlation analysis also showed that CH4production was positively correlated with the abundance of SGMT clade methanogens.The results indicated that the shift in the methanogen community profile at the subgenus level were more effective in producing CH4,which could help explain the lower amount of CH4emissions in the steers after withdrawing TTS from the diet,with a reduction in predominant methanogen species(SGMT clade methanogens)and increase in minor species(RO clade methanogens).

Conclusions

Adding TSS in the basal diet of beef steers caused significant alterations in the rumen protozoa community structure at the family and genus level,as well as rumen methanogen community structure.After withdrawing TSS,the ciliate community of steers did not return to the same protozoa community prior to TSS supplementation,and the new structure of the rumen microbial community maybe shaped by complex relationship between ciliate protozoa and methanogenc archaea.A number of protozoa genera and methanogen species were found to be correlated with CH4emissions.Therefore,we speculated that decreasing the abundance of protozoa belonging to theIsotrichagenus,and increasing members of theMetadiniumgenus,are likely to reduce CH4emissions.

Supplementary information

Supplementary informationaccompanies this paper at https://doi.org/10.1186/s40104-020-00491-w.

Additional file 1 Table S1.Table S1 Composition and nutrient concentration of the basal diet.

Additional file 2 Table S2.The primers used in the construction of Illumina sequencing library.

Additional file 3 Figure S1.Relative abundance of ruminal protozoa genera of steers at the family level(A)and the genus level(B).

Additional file 3 Figure S2.The community composition of rumen protozoa at the order(A),family(B)and genus(C)level among three periods.

Abbreviations

TSS:Tea seed saponins;BD:Basal diet;BDP:Basal diet post;OTU:Operational taxonomic unit;SGMT:M.smithii,M.gottschalkii,M.millerae and M.thauerii;RO:M.ruminentium and M.olleyae

Acknowledgements

We would like to thank scholarships from the Chinese Scholarship Council(CSC)awarded to Cui Tan for 12 months of cooperative research at CSIRO.We also thank Jagadish.Padmanabha for assistance on sample analyzing.

Authors’contributions

CT,ZW and CM designed the experiments.CT and CR performed the experiments.CT,AMS and RH analyzed the experimental data.CT,CR and CM wrote the draft manuscript,ZW,MB,CR and CM reviewed and edited the manuscript.The author(s)read and approved the final manuscript.

Funding

This study was supported by grants from the Natural Science Foundation of China(31272472),and project grants from Meat &Livestock Australia,Fisheries and Forestry in Australia and Commonwealth Scientific and Industrial Research Organization(CSIRO).

Availability of data and materials

All sequencing data are available in the NCBI Sequence Read Archive(SRA),under the bioproject number PRJNA574788.

Ethics approval and consent to participate

Care of animals and experimental procedures were carried out following the approved CSIRO Animal Ethic Committee protocol No A12/2012.

Consent for publication

Not applicable.

Competing interests

All the authors declare they have no competing interest.

Author details

1Institute of Animal Nutrition,Sichuan Agricultural University,Chengdu 611130,Sichuan,China.2Commonwealth Scientific and Industrial Research Organisation,CSIRO Agriculture and Food,Australian Tropical Sciences and Innovation Precinct,James Cook University,Townsville QLD 4811,Australia.3Present address:CR Eco-efficient Agriculture Consultancy(CREAC),46 Bilbao Place,Bushland Beach QLD 4818,Australia.4Department of Livestock Production,Shaheed Benazir Bhutto University of Veterinary and Animal Sciences,Sakrand,Sindh 67210,Pakistan.5“Low Carbon Breeding Cattle and Safety Production”,University Key Laboratory of Sichuan Province,Ya’an 625014,Sichuan,China.6School of Biosciences,University of Nottingham,Sutton Bonington,Loughborough LE12 5RD,UK.7CSIRO Agriculture,Queensland BioScience Precinct,St Lucia,Brisbane QLD 4067,Australia.

Received:1 March 2020 Accepted:25 June 2020