Aline F. O. Ramos1
Stephanie A. Terry2
Devin B. Holman3
Gerhard Breves4Luiz G. R. Pereira5
André G. M. Silva1
Alexandre V. Chaves2*- 1Animal Science Graduate Course, Veterinary Medicine Institute, Federal University of Pará, Belém, Brazil
- 2School of Life and Environmental Sciences, Faculty of Science, The University of Sydney, Sydney, NSW, Australia
- 3Lacombe Research and Development Centre, Agriculture and Agri-Food Canada, Lacombe, AB, Canada
- 4Institute of Physiology, University of Veterinary Medicine, Hanover, Germany
- 5Embrapa Dairy Cattle, Juiz de Fora, Brazil
Tucumã oil is sourced from the fruit pulp of the tucumã tree and contains high concentrations of unsaturated fatty acids and carotenoids. Due to these properties it may have the potential to decrease enteric methane (CH4) from ruminants when included in the diet. The objective of this study was to determine the effect of oil mechanically extracted from the fruit pulp of tucumã on fermentation characteristics, CH4 production and the microbial community using the rumen stimulation technique. Treatments consisted of a control diet (forage:concentrate; 70:30), and tucumã oil included at 0.5 or 1.0% (v/v). Addition of tucumã oil linearly decreased (P < 0.01) dry matter disappearance. Total gas (mL/d) and carbon dioxide (CO2) production (mL/d, mL/g DM) were unaffected (P ≥ 0.36) to increasing addition of tucumã oil where 0.5% (v/v) of Tucumã oil numerically increased both variables. Acetate and butyrate percentages of total VFA were linearly decreased (P ≤ 0.01) and propionate and valerate percentages of total VFA were linearly increased (P < 0.01) by increasing concentrations of tucumã oil added to the substrate. The ratio of acetate to propionate was linearly decreased (P < 0.01) with increasing concentration of tucumã oil. Methane production (mL/d) was linearly decreased (P = 0.04) with increasing addition of tucumã oil to the substrate. Tucumã oil reduced the bacterial richness and diversity when included at 1.0% (v/v) in both solid- and liquid- associated microbes. The abundance of the genera Fibrobacter and Rikenellaceae RC9 gut group were decreased and Pyramidobacter, Megasphaera, Anaerovibrio, and Selenomonas were enriched by the addition of 1.0% tucumã oil. In conclusion, tucumã oil resulted in the favorable shift in fermentation products away from acetate toward propionate, decreasing the production of CH4 when tucumã oil was included at 1.0% (v/v), however, substrate digestibility was also inhibited. The rumen microbiota was also altered by the addition of tucumã oil.
Introduction
The enteric production of methane (CH4) from livestock accounts for approximately 40% of total greenhouse gases (GHG) emitted from livestock production systems (Gerber et al., 2013). Brazil has the third largest cattle population in the world, making significant contributions to increasing GHGs emitted into the atmosphere. Methane emitted from cattle also constitutes as a 2–12% loss of gross energy intake (Johnson and Johnson, 1995) providing incentive for producers to decrease enteric CH4 production.
Enteric production of CH4 in ruminants is directly affected by microbial communities present in the rumen producing varing amounts of CH4 depending on the composition of the diet being fed (Henderson et al., 2015). Dietary supplementation with plant derived oils is considered one of the most effective methods for suppressing ruminal methanogenesis (Beauchemin et al., 2007). The mechanisms of action depend on the composition and type of fat but include depression of ciliate protozoa and methanogen populations, dilution through replacement of fermentable carbohydrates, biohydrogenation of free unsaturated fatty acids and reduction of ruminal organic matter fermentation (Eugene et al., 2008; Beauchemin et al., 2009; Knapp et al., 2014). Not only does oil supplementation provide energy to the diet, many sources of fat contain a range of secondary compounds which have the potential to further inhibit ruminal methanogenesis (Delgado et al., 2012). Fats may also be employed to alter the fatty acid composition of ruminant products (milk, meat) to improve its associated healthiness (Adeyemi et al., 2016).
Tucumã is one of the many oleaginous palms spread throughout the Amazonian region of South America. It belongs to the genus Astrocaryum of which 26 species are native to the southeast Amazon basin (Alexandre et al., 2015). It produces fruit in which the pulp can be pressed to extract oil for human consumption, with the dried oilseed press cake being used for animal feed (Ramos et al., 2011). Its extractivism plays an important economic role for people within these areas, especially with interest in its use as a biofuel (Bora et al., 2001). Tucumã oil is known to be rich in low molecular weight fatty acids as well as carotenoids which are naturally occurring antioxidant pigments (Bora et al., 2001). Oleic acid is the main constituent in tucumã oil extracted from the fruit pulp, followed by palmitic acid (Ferreira et al., 2008). Oleic acid is an unsaturated fatty acid, which can act as a hydrogen sink through biohydrogenation (Johnson and Johnson, 1995), however, this role in decreasing CH4 emissions is considered minor (Jenkins et al., 2008).
We hypothesized that the inclusion of tucumã oil would decrease CH4 emissions and alter the ruminal bacterial community. As such, the objective of this study was to investigate the effects of tucumã oil supplementation on in vitro rumen fermentation parameters, CH4 production and the rumen microbiota using the rumen stimulation technique (RUSITEC).
Materials and Methods
The donor cows used in this experiment were cared for in accordance with the guidelines of the German Animal Welfare Act approved by the Lower Saxony State Office for Consumer Protection and Food Safety (LAVES, approval number AZ 33.4-42505-04-13A373).
Experiment Design and Treatments
The experiment was conducted as a completely randomized design with three treatments duplicated in two runs with two replicates per run. The three treatments consisted of a control (no tucumã oil inclusion) and two different inclusions of tucumã oil in the substrate DM at 0.53% (0.5%) and 1.07% (1.0%) of liquid vessel volume (4 and 8 mL/d, respectively; 0.38 and 0.77 mL/g substrate DM). The experimental period was 15 days with day 1–7 used for adaptation and day 8–15 used for measurements. Tucumã oil was administered throughout the whole experimental period.
The substrate used was a hay:concentrate (70:30 DM basis; 92.3% DM, 13.5% CP, 62% NDF, and 5.2% ash – DM basis; Table 1) using hay obtained from natural grassland of Lower Saxony, Germany. Hay was prepared using an electrical clipper with a 76-mm blade (Duarte et al., 2017a). The commercial concentrate was pelleted (Deuka Schaffutter, Deutsche Tiernahrung Cremer, Düsseldorf, Germany) and contained 0.9% calcium, 0.55% phosphorus and 0.2% sodium. Both the hay and substrate were weighed into the same nylon bag (10 cm × 5 cm, pore size 50 ± 10 μm; in situ nylon bags Ankom Technology, Macedon NY, United States) for a total mass of 10.4 g of substrate (DM-basis).
TABLE 1. Chemical composition of the hay and concentrate (%DM) and fatty composition of Tucumã oil.
The fruit of the tucumã palm (Astrocaryum vulgare Mart) was obtained from a closed forest at Santo Antônio do Tauá in the state of Pará, Brazil during the fruiting period. It is typical palm of the Amazon, used as food source and income of the local communities. The crude or refined oi is commercialized internationally. Fruits were sorted for quality, washed, sanitized, and then dried at room temperature in polyethylene plastic packaging. Pulp was removed from the fruit and distributed in trays to dehydrate at 40°C for 72 h. The dried pulp was then cooled at room temperature and then crushed into a fine powder, using a processor. To obtain the oil, the powder was then pressed using a commercial hydraulic press (Marconi, model ME 098/A, Piracicaba, SP, Brazil) at room temperature with the initial and final pressure at 7500 and 30,000 kg/cm2, respectively (Ferreira et al., 2008). After pressing, press cake (pulp) is separated from the oil. Tucumã oil has a dark orange color, rich in beta-carotene, 0.982 g/L density, melting point 27–35°C which is in the liquid form, and the saponification index of 188.4 mg KOH/g. A profile of the major fatty acids in the oil were provided by the Amazon Oil Industry (Ananindeua, Pará, Brazil) and are shown in Table 1.
Inoculum Sampling and Incubation Procedure
Rumen inoculum was obtained from two ruminally cannulated Holstein cows, 2 h after morning feeding. Cattle were fed grass hay ad libitum and 600 g/d of a commercial concentrate. These were the same feeds used as substrate. Rumen contents were separated into rumen fluid and solid content by straining through gauze. Samples for DNA were taken from the solid and liquid fractions from each cow and frozen in liquid nitrogen (-196°C). Samples were then stored at -40°C until they were placed in a freeze dryer (over 48 h). Once dried, samples were kept at -20°C until DNA extraction.
Fluid samples from each cow were pooled together and the pH and redox potential were recorded. Samples (2 mL) were also taken and stored at -20°C for determination of volatile fatty acids (VFA; Duarte et al., 2017a) and ammonia (NH3-N; Duarte et al., 2017a).
Prewarmed 800 mL fermentation vessels were placed in the RUSITEC apparatus (Czerkawski and Breckenridge, 1977) and water was kept at 39°C. Each fermentation vessel had an inner vessel which contained one nylon bag with 70 g of solid digesta, and one bag containing the basal substrate in which the tucumã oil was added immediately before incubation. Tucumã oil was inserted into the bags using a pipette (4 and 8 mL for 0.5 and 1.0% (v/v) treatments, respectively). Each fermenter was filled with approximately 750 mL of rumen fluid and infused with McDougall’s buffer at a dilution rate of 30 mL/h. The inner vessels were continuously moved up and down by an electric motor to ensure adequate mixing between fluid and particles. After the first 24 h of incubation, the bag with the solid rumen digesta was replaced with a bag containing the substrate. Bags were replaced with a fresh bag containing feed after 48 h of incubation, replacing one bag per day. Bags from day 15 were not used for DM determination as they were only incubated for 24 h. Effluent was collected in 2 L glass flasks which were kept on ice to arrest bacterial growth and fermentation.
Sample Collection
Dry matter disappearance (DMD) at 48 h was determined on day 8 and day 10–13 when bags were not used for DNA extraction. After removal from the vessel, feed bags were washed in 50 mL of pre-warmed buffer in a small plastic bag, gently squeezed and the residual buffer was placed back into the fermenter to ensure transfer of solid-phase-associated microorganisms. The residual feed bag was rinsed under cold water until the water was clear and then dried at 55°C for 48 h for the determination of DMD (Duarte et al., 2017a).
Total daily gas production was collected in gas-tight bags (Plastigas, Linde AG, München, Germany). From day 8 to 15, before measurement of total gas, two 20 mL aliquots were taken from the septum of each gas bag and transferred into evacuated tubes for the analysis of CH4 and carbon dioxide (CO2). Total daily gas production was measured using a drum-type meter (Ritter Apparatebau, Bochum, Germany).
During bag exchange, fermenter pH, gas production and effluent volume for each fermenter was measured. The pH and redox potential of the vessel was measured daily during bag exchange using a Knick pH meter (digital pH meter 646, Knick, Berlin, Germany). Effluent from each fermenter was measured and two samples (2 mL) of effluent were taken and stored at -40°C until analyzed for VFA and NH3-N. Ammonia and VFA data are presented as concentrations, molar percentages of total VFA as well as total production. Daily ammonia and VFA production were calculated by multiplying NH3-H and VFA concentrations by the effluent volume (Duarte et al., 2017a,b).
DNA Extraction
On day 5, 10, and 15 nylon bags (solid associated microbes; SAM), as well as 30 mL of fermenter liquid (liquid associated microbes; LAM – d15 only) were removed from each vessel and immediately placed in liquid nitrogen for later extraction of DNA. Samples were taken across 3 days (e.g., day 5, 10, and 15) to evaluate differences in microbial communities across the experiment. Samples were stored at -40°C until they were placed in a freeze dryer (over 48 h). Samples were then finely ground using a coffee grinder and placed back into the freezer until DNA extraction. The liquid samples were freeze dried for 4 day and then ground using a mortar and pestle.
Total DNA was extracted from each sample using a QIAamp Fast DNA stool mini kit (Qiagen, Hilden, Germany), according to the manufacturer’s instructions. DNA yield and purity were measured using a NanoDrop spectrophotometer (ThermoFisher Scientific, Waltham, MA, United States). Extracted DNA was stored at -20°C until 16S rRNA gene library preparation and sequencing.
Sequencing of the 16S rRNA Gene
The modified 515-F and 806-R primers as found in Walters et al. (2016) were used to PCR amplify the V4 hypervariable of both the archaeal and bacterial 16S rRNA gene. The PCR conditions and 16S rRNA gene sequencing was as previously described (Duarte et al., 2017a). Briefly, a two-step PCR was used to generate the 16S rRNA gene amplicons and these amplicons were then sequenced on an Illumina MiSeq instrument (Illumina, Inc., San Diego, CA, United States) using the MiSeq Reagent Kit v2 (500 cycles) (Illumina, Inc.), according to manufacturer’s instructions.
The 16S rRNA gene sequences were processed using the R-package DADA2 (v. 1.4) (Callahan et al., 2016) and included primer removal and truncating both the forward and reverse reads at 225 bp. The pair-end reads were then merged, and chimera sequences removed. The SILVA SSU database v. 128 (Quast et al., 2012) and the RDP naïve Bayesian classifier (Wang et al., 2007) with a 50% confidence threshold were used to assign taxonomy to each inferred 16S rRNA gene sequence; defined here as an operational taxonomic unit (OTU) with 100% sequence similarity. QIIME v. 1.9.1 (Caporaso et al., 2010) was used to measure richness (number of OTUs) and the Shannon diversity index. Bray–Curtis dissimilarities were calculated using the R packages vegan (v. 2.4.4) (Oksanen et al., 2017) and phyloseq (v. 1.20.0) (McMurdie and Holmes, 2013).
All 16S rRNA gene sequences were deposited into the NCBI Sequence Read Archive under BioProject accessions PRJNA416148.
Chemical Composition
Feed was analyzed following the Association Official Analytical Chemistry [AOAC], 1998 method for DM (method 967.03). Neutral detergent fiber (NDF) content was analyzed according to Van Soest et al. (1991) with the use of sodium sulfite and heat-stable α-amylase. Methane and CO2 was measured by using gas chromatography (GC, 2014, Shimadzu Europa GmbH, Duisburg, Germany) and CH4 and CO2 production was calculated by multiplying the total gas volume by the percentage of CH4 with correction for temperature and pressure (0°C, 101.3 kPa; Riede et al., 2013). Determination of NH3-N concentrations was carried out as described previously by Riede et al. (2013). Briefly, 1 mL of sample was centrifuged at 4600 × g for 10 min. From the supernatant, 50 mL were mixed with 5 mL phenol solution (106 mM phenol, 0.17 mM sodium nitroprusside dehydrate) and 5 mL sodium hypochlorite solution (1% v/v sodium hypochlorite; 125 mM NaOH) and kept at room temperature for 10 min. After an incubation step at 60°C for 10 min, NH3-N concentration was determined photometrically at 546 nm in a spectrometer (DU 640, Beckman Coulter GmbH, Krefeld, Germany) using a NH4Cl standard solution (5 mM). Volatile fatty acids concentrations were analyzed by a gas chromatography system (model 5890 II, Hewlett Packard, Böblingen, Germany) equipped with a 1.8 m × 2 mm glass column packed with Chromosorb WAW (mesh 80/100) with 20% neopentyl glycol succinate and 2% orthophosphoric acid. Helium was used as a carrier gas with a flow rate of 25 mL/min. Injection port, detector and oven temperatures were 220°C, 250°C, and 130°C, respectively. The daily VFA production was estimated by multiplying VFA concentration by the volume of effluent (Duarte et al., 2017b).
Statistical Analysis
Fermentation Data
The univariate procedure of SAS (SAS, Inc. 2018; SAS Online Doc 9.1.4) was used to test for normal distribution of the data. Fermentation data were analyzed using the MIXED procedure of SAS (SAS, Inc., 2018; SAS Online Doc 9.1.4). The model included the fixed effects of treatment, day and treatment × day interactions with the day of sampling (e.g., day 8–15) from each fermenter treated as a repeated measure. Therefore, the individual fermenter (n = 4 per treatment) was used as the experimental unit for statistical analysis (Avila-Stagno et al., 2014). The method for computing denominator degrees of freedom was Kenwardroger. The DDFM = KENWARDROGER option performs the degrees of freedom calculations detailed by Kenward and Roger (1997, 2009). The minimum values of Akaike’s information criterion were used to select the covariance structure among Compound Symmetry, Heterogeneous Compound Symmetry, Autoregressive, Heterogeneous Autoregressive, Toeplitz, Heterogeneous Toeplitz, Ante-dependence, Unstructured and Banded for each parameter. When P ≤ 0.05 for Type III fixed effects for treatment, orthogonal polynomial contrasts were carried out to test for linear and quadratic responses to increasing concentrations of tucumã oil [0, 0.5%, and 1.0% (v/v)]. Since there was no quadratic effect (P ≥ 0.05) observed in any parameter tested, these P-values were omitted from the tables. Treatment means were compared using the least squares mean linear hypothesis test (LSMEANS/DIFF). Significance was declared at P ≤ 0.05 and tendencies 0.05 < P ≤ 0.10.
Microbiome Data
All samples were randomly subsampled to 44,500 sequences prior to analysis to account for differences in sequencing depth. The number of OTUs per sample and the Shannon index were compared by treatment group using a linear mixed model in R v. 3.4.2 (R Core Team, 2017) using lme4 v. 1.1.12 (Bates et al., 2014) with fermenter as the random effect and tucumã oil treatment, and sampling time (e.g., day 5, 10, and 15) as fixed effects, followed by Tukey’s Honest Significant Difference (Lenth, 2016). Both the archaeal and bacterial community structure was analyzed using permutational multivariate analysis of variance (PERMANOVA) and the Adonis function with 10,000 permutations in vegan. The betadisper function in vegan was used to assess the homogeneity of dispersion for each time point. Linear discriminant analysis (LDA) effect size (LEfSe) (Segata et al., 2011) was used to identify genera with a relative abundance of greater than 0.1% that were differentially abundant between the control and 1% tucumã oil treatments for both LAM and SAM samples at day 15. A minimum LDA score of 4.0 was used as the threshold for classifying differentially abundant genera.
Results
Effect of Tucumã Oil on in Vitro Fermentation
Dry matter disappearance was linearly decreased (P < 0.01) with increasing addition of tucumã oil to the substrate (Table 2). Acetate and butyrate percentages of total VFA (mmol/100 mmol) linearly decreased (P ≤ 0.01) and propionate and valerate percentages of the total VFA were linearly increased (P < 0.01) with increasing concentration of tucumã oil. Consequently, the ratio of acetate to propionate (P < 0.01) linearly decreased with increasing addition of tucumã oil to the substrate. In agreement with the molar percentages of VFA data (Table 2), propionate production (mmol/d) tended (P = 0.08) to increase, and butyrate and valerate production (mmol/d) increased (P ≤ 0.05) with increasing concentration of tucumã oil (Supplementary Table S1). There was a significant interaction between treatment and day for acetate production (mmol/d), where 1% tucumã supplementation decreased (P < 0.01) acetate production at days 10 and 12–14 compared to control (interaction means not shown). There was also a significant treatment by day interaction for ammonia production (mmol/d), where 1% tucumã supplementation decreased ammonia production compared to the control at day 10 only (interaction means not shown)
TABLE 2. Effect of Tucumã oil on dry matter digestibility, pH, redox, the molar percentages of individual volatile fatty acids (VFA), ammonia and total VFA concentrations in a Rusitec fed a mixed hay – concentrate diet.
Total gas production (mL/d) was unaffected (P = 0.45) by tucumã oil supplementation (Table 3). Similarly, tucumã oil addition did not affect (P ≥ 0.36) production of CO2 (mL/d, mg/g DM, mg/g DMD). Methane production (mL/d) was linearly decreased (P = 0.05) where production at 1% tucumã oil was halved compared to the control and 0.5% inclusion of tucumã oil (Table 3). When CH4 was expressed as mg CH4/g DM and mg CH4/g DM disappeared, treatment × day interaction was significant. Only on day 12 and 15, mg CH4/g DM was lower in treatment 1% tucumã compared to the control. For mg CH4/g DM disappeared, CH4 was only different on day 12 and 13 within the 0.5% (v/v) treatment, with no other differences observed between treatments.
TABLE 3. Effect of Tucumã oil on gas production in a Rusitec fed a mixed hay – concentrate diet.
Effect of Tucumã Oil on the Rumen Microbiota
The LAM and SAM samples had significantly different bacterial community structures (R2 = 0.08; P ≤ 0.001); however, most of the abundant genera were shared between the two sample types (Supplementary Figures S1, S2). This included Prevotella, Megasphaera, Anaerovibrio, Rikenellaceae RC9 gut group, Fibrobacter, Lactobacillus, and Selenomonas.
Among SAM samples, sampling time had a slightly larger effect (Figure 1; R2 = 0.099; P < 0.001) than tucumã oil supplementation on the structure of the rumen microbial community (R2 = 0.089; P < 0.001). The greatest effect of tucumã treatment was observed at day 15 (Supplementary Figures S3, S4; R2 = 0.240; P = 0.007) as the control samples were most dissimilar from both tucumã oil treatment samples at this point (Supplementary Figure S3).
FIGURE 1. Principal coordinates analysis plots of the Bray-Curtis dissimilarities for solid-associated microbes (SAM) by treatment and sampling time. Percentages of variation explained by the principal coordinates are indicated on the axes.
Tucumã oil also affected the richness (number of OTUs) and diversity (Shannon index) of the ruminal microbiota. The number of OTUs was significantly decreased in the 1.0% tucumã oil samples at day 10 and day 15 for SAM (Figure 2A) and day 15 for LAM (Figure 3A) compared with the control treatment. The Shannon diversity index was also lower in the 1.0% tucumã oil LAM samples compared with both the 0.5% tucumã oil and control treatment at day 15 (Figure 3B). Among the SAM samples, only at day 10 did the 1.0% tucumã treatment significantly affect the Shannon diversity index (Figure 2B). The addition of tucumã oil also had a comparable effect in the LAM samples (day 15) (Figure 4; R2 = 0.250; P = 0.006). In fact, the three treatment groups were more dissimilar to each other in terms of Bray–Curtis dissimilarities than at any sampling time for the LAM samples (Supplementary Figure S5).
FIGURE 2. Box plots of the (A) Number of OTUs and (B) Shannon diversity index by treatment and sampling time for solid-associated microbes (SAM). Different lowercase letters within each sampling time indicate significantly different means (P < 0.05) n = 4.
FIGURE 3. Box plots of the (A) number of OTUs and (B) Shannon diversity index by treatment for liquid-associated microbes (LAM; day 15). Different lowercase letters indicate significantly different means (P < 0.05).
FIGURE 4. Principal coordinates analysis plots of the Bray-Curtis dissimilarities for liquid-associated microbes (LAM) samples by treatment at day 15. Percentages of variation explained by the principal coordinates are indicated on the axes n = 4 (P < 0.05; LDA score ≥ 4.0).
Given that the greatest dissimilarity among samples was observed at day 15 between the 1.0% tucumã oil and control treatments, we identified genera that were most strongly associated with each treatment group at day 15 for both LAM and SAM samples (LDA score > 4.0; Table 4). Many of the differentially abundant genera between the two treatments were the same for both LAM and SAM samples including; Fibrobacter, Rikenellaceae RC9 gut group, Megasphaera, Selenomonas, Anaerovibrio, and Pyramidobacter. Tucumã oil appeared to have the greatest positive effect on Olsenella (33.6-fold increase) and Pyramidobacter (23.5-fold increase) in the LAM samples.
TABLE 4. Differentially abundant genera identified between the control and 1% Tucumã oil treatment for liquid- and solid- associated microbes (LAM and SAM, respectively) in a RUSITEC system with a mixed hay – concentrate diet.
Only three methanogenic genera were detected including Methanobrevibacter, Methanosphaera, and Methanomicrobium. There was no effect on the relative abundance on any of the detected methanogenic genera in either LAM and SAM samples at any sampling time (Supplementary Table S2).
Discussion
Tucumã oil was included at 0.50 and 1.0% (v/v) of total fermenter volume due to the dilution factor of the oil in the RUSITEC. Duarte et al. (2017b) suggested when they included Pequi oil into a Rusitec fermenter at 6% dietary fat content of the TMR, that when expressed as the total fermenter volume, oil was only included at 0.18% of total fermenter volume. It was suggested that this may have contributed to conflicting results between the Rusitec and a previous in vitro batch culture they had conducted. Whilst necessary for in vitro experiments, this expression of treatment is not extrapolable to the rumen diet as when calculated on a DM basis these concentrations constitute as 27.4 and 43.0% of total added substrate (i.e., TMR and Tucumã oil).
Dry matter disappearance was decreased by up to 25.6% in the 1.0% tucumã oil treatment compared to the control substrate which is not a desirable effect. This reduction in DM disappearance coincides with a decrease in the richness and diversity of the ruminal microbiota on day 10 and 15. Fibrobacter, a genus of fibrolytic bacteria were also lower in relative abundance in the 1.0% tucumã oil treatment which is consistent with the negative effect that oils have on fiber-digesting microbes within the rumen (McIntosh et al., 2003). The relative abundance of Rikenellaceae RC9 gut group was also negatively affected by tucumã oil. Similarly, Zened et al. (2013) found that the abundance of Rikenellaceae RC9 gut group was decreased by 79% in the rumen of cattle fed a diet supplemented with 5% DM sunflower oil. Rikenellaceae is a relatively new bacterial family and therefore its metabolic function in the rumen has not yet been defined; however, Pitta et al. (2010) suggested that members of this family are associated with either the primary or secondary degradation of carbohydrates as its abundance was decreased as increased proportions of bermuda grass were substituted for wheat in a cattle diet.
Tucumã oil has been evaluated for its anti-inflammatory properties due to its high carotenoid content and has been shown to effectively improve the immune system in mice (Baldissera et al., 2017). The Tucumã plant also possesses chemical antioxidant compounds which could be involved in the decrease in the relative abundance of Fibrobacter, Rikenellaceae RC9 gut group and Lachnospiraceae. Similarly, Jobim et al. (2014) indicated that tucumã extracts presented antimicrobial activities against gram positive bacteria including Enterococcus faecalis, Bacillus cereus, and Listeria monocytogenes and this was related to a combination of secondary compounds including polyphenols and carotenes and their ability to disrupt redox reactions (Jobim et al., 2014).
Oleic acid is the main fatty acid constituent in tucumã oil followed by palmitic, linolenic, stearic and linoleic acid. Oleic acid is a monounsaturated omega-9 fatty acid and is found in a variety of oils including olive, canola, peanut, and high oleic sunflower and safflower oils (He et al., 2012). The hydrogenation of unsaturated fatty acids from oil is suggested to contribute to the decrease in the ability of microorganisms to saturate fatty acids, which causes unsaturated fatty acids to accumulate and interfere with microbial digestion. Selenomonas spp., which were enriched in the 1.0% tucumã oil treatment, can hydrate oleic acid resulting in enhanced growth for members of this genus (Maczulak et al., 1981; Hudson et al., 1995). Oleic acid can decrease ruminal CH4 production in vitro (Wu et al., 2016); however, as with the present study, there have also been reductions in DM degradability observed when oils high in oleic acids have been fed to cattle (Hristov et al., 2011). This may be related to the inhibitory effect of a variety of oils on fibrolytic bacteria (Enjalbert et al., 2017). Wu et al. (2016) demonstrated that total gas and CH4 production were decreased when oleic acid was included at up to 60 mg/50 mL culture solution in a batch fermentation system. In the current study, CH4 (mL/d) was decreased by 49.7% in the 1.0% tucumã oil treatment compared to the control; however, this was in part a result of an undesirable decrease in DM disappearance. Despite the change in CH4 production (mL/d), there was no effect on the relative abundance of any the methanogenic genera (Supplementary Table S2).
The decrease in CH4 production was associated with shifts in rumen fermentation favoring the production of propionate over acetate. It is well established that an increase in the ratio of propionate to acetate favors a decrease in the production of CH4. This is because propionate acts as alternative hydrogen sinks which decreases the availability of hydrogen to CH4 producing Archaea, known as methanogens (McAllister and Newbold, 2008). The observed increase in propionate production is consistent with an increase in the abundance of Anaerovibrio in both LAM and SAM samples as Hobson et al. (1981) indicated that these bacteria play an important role in lipolytic activity in sheep, producing glycerol for propionate synthesis (Jenkins et al., 2008). The relative abundance of Anaerovibrio was also increased by pequi oil in a RUSITEC study, oil which also has a high oleic acid content (Duarte et al., 2017b). Despite the change in CH4 production (mL/d), there was no effect on the relative abundance of any the methanogenic genera identified in this study, that is, Methanobrevibacter, Methanosphaera, and Methanomicrobium (Supplementary Table S2). However, although the PCR primers used in the present study target both the archaeal and bacterial 16S rRNA gene, they are not specific for the methanogenic populations and therefore it is possible that their abundance is underestimated in the present study.
Interestingly, Pyramidobacter spp., which belong to the Synergistetes phylum were increased by 23.5 fold in LAM and 7.5 fold in SAM samples from the 1.0% tucumã oil treatment. These bacteria are obligatory anaerobic, produce major amounts of acetic acid, and are associated with increased fiber digestion (Downes et al., 2009; Pan et al., 2017). They have also been known to increase in abundance in high CH4 emitting cattle (Wallace et al., 2015). This contradicts the observed decrease in the proportion of acetic acid and DM disappearance when 1.0% of tucumã oil was included in the substrate, and it remains unclear as to why this group of bacteria was increased as a result of Tucumã oil supplementation. Prevotella which are involved in the metabolism of proteins and peptides in the rumen (Wallace et al., 1997), were also increased in the SAM samples. Prevotella spp. also produce propionate (Russell and Hespell, 1981; McCabe et al., 2015) and this is consistent with the increase in propionic acid production with increasing addition of tucumã oil.
The addition of tucumã oil also decreased the production of butyrate (Supplementary Table S1). Butyrate production in the bovine rumen is associated with certain members of the bacterial and protozoal community (Demeyer and Van Nevel, 1979; Carberry et al., 2012). The protozoa are positively linked with an increase in CH4 production as the principal fermentation end products of protozoal activity are acetate and butyrate, and these two VFA provide H2 for methanogenesis (McAllister and Newbold, 2008). The reduction in butyrate production was simultaneous with an increase in valerate production and both VFA have competitive metabolic pathways (Kristensen et al., 2000; Kristensen and Harmon, 2004). Many members of the Lachnospiraceae family are major producers of butyrate (Meehan and Beiko, 2014) as evidenced by the lower relative abundance of Lachnospiraceae XPB1014 group in the SAM samples from the 1.0% tucumã oil treatment. The concentration of NH3-N was not affected when tucumã oil was included.
Conclusion
The addition of tucumã oil resulted in the favorable shift in fermentation products away from acetate toward propionate, decreasing the production of CH4 when tucumã oil was included at 1.0%. However, dry matter digestibility was also negatively affected by tucumã oil addition into the substrate which is undesirable. The structure of the rumen microbiota was significantly altered and the bacterial richness and diversity decreased in the 1.0% tucumã oil treatment for both LAM and SAM samples. Tucumã oil lowered the abundance of the genera Fibrobacter and Rikenellaceae RC9 gut group and enriched Pyramidobacter, Megasphaera, Anaerovibrio, and Selenomonas. Whilst tucumã oil may be used in the ruminant diet to suppress CH4 production, it needs to be further evaluated using in vitro techniques at lower concentrations to investigate whether diet digestibility can be affected.
Author Contributions
AC, GB, and AR: study design. AR, ST, and AC: conducting RUSITEC study. AR, ST, DH, GB, and AC: lab analysis. AR and ST: DNA extraction. DH: bioinformatics. ST, AR, DH, GB, and AC: writing the manuscript. ST, AR, DH, GB, LP, AS, and AC: manuscript revision. AS: paying for publication.
Funding
We thank CAPES and Meat and Livestock Australia (MLA) for financial support of AR and ST, respectively, during their Ph.D. candidatures. We thank the Alexander von Humboldt Foundation for supporting AC during his fellowship. Special thanks to Marion Burmester for her lab assistance, as well as students and staff at the Department of Physiology, University of Veterinary Medicine (Ti-Ho), Hannover (Germany). We also thank the Pró-Reitoria de Pesquisa e Pós-Graduação da Universidade Federal do Pará (PROPESP UFPA) for the payment of this publication.
Conflict of Interest Statement
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Supplementary Material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fmicb.2018.01647/full#supplementary-material
FIGURE S1 | Relative abundance (%) of the ten most abundant genera in the solid-associated microbe (SAM) samples by treatment and sampling time.
FIGURE S2 | Relative abundance (%) of the ten most abundant genera in the liquid-associated microbe (LAM) samples by treatment at day 15.
FIGURE S3 | Principal coordinates analysis plot of the Bray–Curtis dissimilarities for SAM samples by treatment at day 15. Percentages of variation explained by the principal coordinates are indicated on the axes.
FIGURE S4 | Bray–Curtis dissimilarities between treatments by sampling time for SAM samples. Error bars represent standard error of the mean (n = 16 comparisons). Different lowercase letters within each sampling time indicate significantly different means (P < 0.05).
FIGURE S5 | Bray–Curtis dissimilarities between treatmentsfor LAM samples. Error bars represent standard error of the mean (n = 16 comparisons).
TABLE S1 | Effect of tucumã oil on volatile fatty acids and ammonia production in a RUSITEC fed a mixed hay – concentrate diet.
TABLE S2 | Percent relative abundance of methanogenic genera detected in LAM and SAM samples by sampling time and treatment (n = 4). ND = not detected.
References
Adeyemi, K. D., Sazili, A. Q., Ebrahimi, M., Samsudin, A. A., Alimon, A. R., Karim, R., et al. (2016). Effects of blend of canola oil and palm oil on nutrient intake and digestibility, growth performance, rumen fermentation and fatty acids in goats. Anim. Sci. J. 87, 1137–1147. doi: 10.1111/asj.12549
Alexandre, E. C. F., Silveira, E. V., Castro, C. F. S., Sales, J. F., De Oliveira, L. C. S., Viana, L. H., et al. (2015). Synthesis: characterization and study of the thermal behavior of methylic and ethylic biodiesel produced from tucumã (Astrocaryum huaimi Mart.) seed oil. Fuel 161, 233–238. doi: 10.1016/j.fuel.2015.08.062
Association Official Analytical Chemistry [AOAC] (1998). Official Methods of analysis, 16th Edn, Vol. 2. Arlington, TX: AOAC, 99.
Avila-Stagno, J., Chaves, A. V., Ribeiro, G. O. Jr., Ungerfeld, E. M., and McAllister, T. A. (2014). Inclusion of glycerol in forage diets increases methane production in a rumen simulation technique system. Br. J. Nutr. 111, 829–835. doi: 10.1017/S0007114513003206
Baldissera, M. D., Souza, C. F., Doleski, P. H., Grando, T. H., Sagrillo, M. R., da Silva, A. S., et al. (2017). Treatment with tucumã oil (Astrocaryum vulgare) for diabetic mice prevents changes in seric enzymes of the purinergic system: improvement of immune system. Biomed. Pharmacother. 94, 374–379. doi: 10.1016/j.biopha.2017.07.113
Bates, D., Maechler, M., Bolker, B., and Walker, S. (2014). Lme4: Linear Mixed-Effects Models Using Eigen and S4. R package version, Vol. 1. 1–23.
Beauchemin, K. A., McGinn, S. M., Benchaar, C., and Holtshausen, L. (2009). Crushed sunflower, flax, or canola seeds in lactating dairy cow diets: effects on methane production, rumen fermentation, and milk production. J. Dairy Sci. 92, 2118–2127. doi: 10.3168/jds.2008-1903
Beauchemin, K. A., McGinn, S. M., and Petit, H. V. (2007). Methane abatement strategies for cattle: lipid supplementation of diets. Can. J. Anim. Sci. 87, 431–440. doi: 10.1155/2010/945785
Bora, P. S., Narain, N., Rocha, R. V. M., De Oliveira Monteiro, A. C., and De Azevedo Moreira, R. (2001). Characterisation of the oil and protein fractions of tucuma (Astrocaryum Vulgare Mart.) fruit pulp and seed kernel. Cienc. Tecnol. Aliment. 3, 111–116. doi: 10.1080/11358120109487654
Callahan, B. J., McMurdie, P. J., Rosen, M. J., Han, A. W., Johnson, A. J. A., and Holmes, S. P. (2016). DADA2: high-resolution sample inference from Illumina amplicon data. Nat. Methods 13, 581–583. doi: 10.1038/nmeth.3869
Caporaso, J. G., Kuczynski, J., Stombaugh, J., Bittinger, K., Bushman, F. D., Costello, E. K., et al. (2010). QIIME allows analysis of high-throughput community sequencing data. Nat. Methods 7, 335–336. doi: 10.1038/nmeth.f.303
Carberry, C. A., Kenny, D. A., Han, S., McCabe, M. S., and Waters, S. M. (2012). Effect of phenotypic residual feed intake and dietary forage content on the rumen microbial community of beef cattle. Appl. Environ. Microbiol. 78, 4949–4958. doi: 10.1128/AEM.07759-11
Czerkawski, J. W., and Breckenridge, G. (1977). Design and development of a long-term rumen stimulation technique (Rusitec). Br. J. Nutr. 38, 371–384. doi: 10.1079/BJN19770102
Delgado, D. C., Galindo, J., GonzaÂlez, R., GonzaÂlez, N., Scull, I., Dihigo, L., et al. (2012). Feeding of tropical trees and shrub foliages as a strategy to reduce ruminal methanogenesis: studies conducted in Cuba. Trop. Anim. Health Prod. 44, 1097–1104. doi: 10.1007/s11250-011-0045-5
Demeyer, D. I., and Van Nevel, C. J. (1979). Effect of defaunation on the metabolism of rumen micro-organisms. Br. J. Nutr. 42, 515–524. doi: 10.1079/BJN19790143
Downes, J., Vartoukian, S. R., Dewhirst, F. E., Izard, J., Chen, T., Yu, W. H., et al. (2009). Pyramidobacter piscolens gen. nov., sp. nov., a member of the phylum “Synergistetes” isolated from the human oral cavity. Int. J. Syst. Evol. Microbiol. 59, 972–980. doi: 10.1099/ijs.0.000364-0
Duarte, A. C., Holman, D. B., Alexander, T. W., Durmic, Z., Vercoe, P. E., and Chaves, A. V. (2017a). The type of forage substrate preparation included as substrate in a RUSITEC system affects the ruminal microbiota and fermentation characteristics. Front. Microbiol. 8:704. doi: 10.3389/fmicb.2017.00704
Duarte, A. C., Holman, D. B., Alexander, T. W., Kiri, K., Breves, G., and Chaves, A. V. (2017b). Incubation temperature, but not pequi oil supplementation, affects methane production, and the ruminal microbiota in a rumen simulation technique (Rusitec) System. Front. Microbiol. 8:704. doi: 10.3389/fmicb.2017.00704
Enjalbert, F., Combes, S., Zened, A., and Meynadier, A. (2017). Rumen microbiota and dietary fat: a mutual shaping. J. Appl. Microbiol. 123, 782–797. doi: 10.1111/jam.13501
Eugene, M., Masse, D., Chiquette, J., and Benchaar, C. (2008). Meta-analysis on the effects of lipid supplementation on methane production in lactating dairy cows. Can. J. Anim. Sci. 88, 331–334. doi: 10.4141/CJAS07112
Ferreira, E. S., Lucien, V. G., Amaral, A. S., and Silveira, C. S. (2008). Physicochemical characterization of the fruit and oil extracted from tucumã (Astrocaryum vulgare Mart.). Alim. Nutr. Araraquara 19, 427–433.
Gerber, P. J., Steinfeld, H., Henderson, B., Mottet, A., Opio, C., Dijkman, J., et al. (2013). Tackling Climate Change Through Livestock – a Global Assessment of Emissions and Mitigation Opportunities. Rome: Food and Agriculture Organization of the United Nations (FAO).
He, M., Perfield, K. L., Green, H. B., and Armentano, L. E. (2012). Effect of dietary fat blend enriched in oleic or linoleic acid and monensin supplementation on dairy cattle performance, milk fatty acid profiles, and milk fat depression. J. Dairy Sci. 95, 1447–1461. doi: 10.3168/jds.2011-4635
Henderson, G., Cox, F., Ganesh, S., Jonker, A., Young, W., Global Rumen Census, C., et al. (2015). Rumen microbial community composition varies with diet and host, but a core microbiome is found across a wide geographical range. Sci. Rep. 5:14567. doi: 10.1038/srep14567
Hobson, P. N., Mann S. O., and Stewart C. S. (1981). Growth and rumen function of gnotobiotic lambs fed on starchy diets. J. Gen. Microbiol. 126, 219–230. doi: 10.1099/00221287-126-1-219
Hristov, A. N., Domitrovich, C., Wachter, A., Cassidy, T., Lee, C., Shingfield, K. J., et al. (2011). Effect of replacing solvent-extracted canola meal with high-oil traditional canola, high-oleic acid canola, or high-erucic acid rapeseed meals on rumen fermentation, digestibility, milk production, and milk fatty acid composition in lactating dairy cows. J. Dairy Sci. 94, 4057–4074. doi: 10.3168/jds.2011-4283
Hudson, J. A., MacKenzie, C. A., and Joblin, K. N. (1995). Conversion of oleic acid to 10-hydroxystearic acid by two species of ruminal bacteria. Appl. Microbiol. Biotechnol. 44, 1–6. doi: 10.1007/BF00164472
Jenkins, T. C., Wallace, R. J., Moate, P. J., and Mosley, E. E. (2008). Board-invited review: recent advances in biohydrogenation of unsaturated fatty acids within the rumen microbial ecosystem. J. Anim. Sci. 86, 397–412. doi: 10.2527/jas.2007-0588
Jobim, M. L., Santos, R. C. V., dos Santos Alves, C. F., Oliveira, R. M., Mostardeiro, C. P., Sagrillo, M. R., et al. (2014). Antimicrobial activity of Amazon Astrocaryum aculeatum extracts and its association to oxidative metabolism. Microbiol. Res. 169, 314–323. doi: 10.1016/j.micres.2013.06.006
Johnson, K. A., and Johnson, D. E. (1995). Methane emissions from cattle. J. Anim. Sci. 73, 2483–2492. doi: 10.2527/1995.7382483x
Kenward, M. G., and Roger, J. H. (1997). Small sample inference for fixed effects from restricted maximum likelihood. Biometrics 53, 983–997. doi: 10.2307/2533558
Kenward, M. G., and Roger, J. H. (2009). An improved approximation to the precision of fixed effects from restricted maximum likelihood. Comput. Stat. Data Anal. 53, 2583–2595. doi: 10.1016/j.csda.2008.12.013
Knapp, J. R., Laur, G. L., Vadas, P. A., Weiss, W. P., and Tricarico, J. M. (2014). Invited review: enteric methane in dairy cattle production: quantifying the opportunities and impact of reducing emissions. J. Dairy Sci. 97, 3231–3261. doi: 10.3168/jds.2013-7234
Kristensen, N. B., and Harmon, D. L. (2004). Splanchnic metabolism of VFA absorbed from the washed reticulorumen of steers. J. Anim. Sci. 82, 2033–2042. doi: 10.2527/2004.8272033x
Kristensen, N. B., Pierzynowski, S. G., and Danfær, A. (2000). Net portal appearance of volatile fatty acids in sheep intraruminally infused with mixtures of acetate, propionate, isobutyrate, butyrate, and valerate. J. Anim. Sci. 78, 1372–1379. doi: 10.2527/2000.7851372x
Lenth, R. V. (2016). Least-squares means: the R package lsmeans. J. Stat. Softw. 69, 1–33. doi: 10.18637/jss.v069.i01
Maczulak, A. E., Dehority, B. A., and Palmquist, D. L. (1981). Effects of long-chain Fatty acids on growth of rumen bacteria. Appl. Environ. Microbiol. 42, 856–862.
McAllister, T. A., and Newbold, C. J. (2008). Redirecting rumen fermentation to reduce methanogenesis. Aust. J. Exp. Agric. 48, 7–13. doi: 10.1071/EA07218
McCabe, M. S., Cormican, P., Keogh, K., O’Connor, A., O’Hara, E., Palladino, R. E., et al. (2015). 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 10:e0133234. doi: 10.1371/journal.pone.0133234
McIntosh, F. M., Williams, P., Losa, R., Wallace, R. J., Beever, D. A., and Newbold, C. J. (2003). Effects of essential oils on ruminal microorganisms and their protein metabolism. Appl. Environ. Microbiol. 69, 5011–5014. doi: 10.1128/AEM.69.8.5011-5014.2003
McMurdie, P. J., and Holmes, S. (2013). phyloseq: an R package for reproducible interactive analysis and graphics of microbiome census data. PLoS One 8:e61217. doi: 10.1371/journal.pone.0061217
Meehan, C. J., and Beiko, R. G. (2014). A phylogenomic view of ecological specialization in the Lachnospiraceae, a family of digestive tract-associated bacteria. Genome Biol. Evol. 6, 703–713. doi: 10.1093/gbe/evu050
Oksanen, J., Blanchet, F. G., Kindt, R., Legendre, P., Minchin, P. R., O’hara, R., et al. (2017). vegan: Community Ecology Package. R Package version 2.4-4.
Pan, X., Xue, F., Nan, X., Tang, Z., Wang, K., Beckers, Y., et al. (2017). Illumina sequencing approach to characterize thiamine metabolism related bacteria and the impacts of thiamine supplementation on ruminal microbiota in dairy cows fed high-grain diets. Front. Microbiol. 8:1818. doi: 10.3389/fmicb.2017.01818
Pitta, D. W., Pinchak, W. E., Dowd, S. E., Osterstock, J., Gontcharova, V., Youn, E., et al. (2010). Rumen bacterial diversity dynamics associated with changing from bermudagrass hay to grazed winter wheat diets. Microb. Ecol. 59, 511–522. doi: 10.1007/s00248-009-9609-6
R Core Team (2017). R: A Language and Environment for Statistical Computing, R Software version 3.4. 2.
Quast, C., Pruesse, E., Yilmaz, P., Gerken, J., Schweer, T., Yarza, P., et al. (2012). The SILVA ribosomal RNA gene database project: improved data processing and web-based tools. Nucleic Acids Res. 41, D590–D596. doi: 10.1093/nar/gks1219
Ramos, S. L. F., Lopes, M. T. G., Lopes, R., Cunha, R. N. V., da Macêdo, J. L. V., de Contim, L. A. S., et al. (2011). Determination of the mating system of Tucumã palm using microsatellite markers. Crop Breed. Appl. Biotechnol. 11, 181–185. doi: 10.1590/S1984-70332011000200011
Riede, S., Boguhn, J., and Breves, G. (2013). Studies on potential effects of fumaric acid on rumen microbial fermentation, methane production and microbial community. Arch. Anim. Nutr. 67, 368–380. doi: 10.1080/1745039X.2013.830518
Russell, J. B., and Hespell, R. B. (1981). Microbial rumen fermentation. J. Dairy Sci. 64, 1153–1169. doi: 10.3168/jds.S0022-0302(81)82694-X
Segata, N., Izard, J., Waldron, L., Gevers, D., Miropolsky, L., Garrett, W. S., et al. (2011). Metagenomic biomarker discovery and explanation. Genome Biol. 12:R60. doi: 10.1186/gb-2011-12-6-r60
Van Soest, P. J., Robertson, J. B., and Lewis, B. A. (1991). Methods for dietary fiber, neutral detergent fiber, and nonstarch polysaccharides in relation to animal nutrition. J. Dairy Sci. 74, 3583–3597. doi: 10.3168/jds.S0022-0302(91)78551-2
Wallace, R. J., Rooke, J. A., McKain, N., Duthie, C., Hyslop, J. J., Ross, D. W., et al. (2015). The rumen microbial metagenome associated with high methane production in cattle. BMC Genomics 16:839. doi: 10.1186/s12864-015-2032-0
Wallace, R. J., McKain, N., Broderick, G. A., Rode, L. M., Walker, N. D., Newbold, C. J., et al. (1997). Peptidases of the rumen bacterium, Prevotella ruminicola. Anaerobe 3, 35–42. doi: 10.1006/anae.1996.0065
Walters, W., Hyde, E. R., Berg-Lyons, D., Ackermann, G., Humphrey, G., Parada, A., et al. (2016). Improved bacterial 16S rRNA gene (V4 and V4-5) and fungal internal transcribed spacer marker gene primers for microbial community surveys. mSystems 1:e00009-15. doi: 10.1128/mSystems.00009-15
Wang, Q., Garrity, G. M., Tiedje, J. M., and Cole, J. R. (2007). Naive Bayesian classifier for rapid assignment of rRNA sequences into the new bacterial taxonomy. Appl. Environ. Microbiol. 73, 5261–5267. doi: 10.1128/AEM.00062-07
Wu, D., Xu, L., Tang, S., Guan, L., He, Z., Guan, Y., et al. (2016). Influence of oleic acid on rumen fermentation and fatty acid formation in vitro. PLoS One 11:e0156835. doi: 10.1371/journal.pone.0156835
Keywords: oil supplementation, cattle, Amazonia, rumen stimulation technique, rumen microbiome
Citation: Ramos AFO, Terry SA, Holman DB, Breves G, Pereira LGR, Silva AGM and Chaves AV (2018) Tucumã Oil Shifted Ruminal Fermentation, Reducing Methane Production and Altering the Microbiome but Decreased Substrate Digestibility Within a RUSITEC Fed a Mixed Hay – Concentrate Diet. Front. Microbiol. 9:1647. doi: 10.3389/fmicb.2018.01647
Received: 05 April 2018; Accepted: 02 July 2018;
Published: 26 July 2018.
Edited by:
Emilio M. Ungerfeld, Instituto de Investigaciones Agropecuarias (INIA), ChileCopyright © 2018 Ramos, Terry, Holman, Breves, Pereira, Silva and Chaves. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Alexandre V. Chaves, alex.chaves@sydney.edu.au