DOI: https://doi.org/10.21203/rs.3.rs-41136/v1
Background: Dietary cation–anion difference (DCAD) has been receiving increased attention in recent years; however, information on the rumen fermentation, cellulolytic bacteria populations, and microbiota of goats fed a low-DCAD diet is less. This study aimed to evaluate the feasibility of feeding a low-DCAD diet for goats with emphasis on rumen fermentation parameters, cellulolytic bacteria populations and microbiota. Growth performance, urine pH, and plasma metabolites were also analyzed as well.
Materials and method: Eighteen goats were randomly allocated to 3 treatments with six replicates of each treatment and 1 goat per replicate. Animals were fed diets with varying DCAD levels at +338 (High DCAD; HD), +152 (Control; CON), and −181 (Low DCAD; LD). This study includes 15-d experimental period and 30-d adaption period.
Results: The DCAD level did not affect the rumen fermentation parameters including pH, buffering capability, acetic acid, propionic acid, butyric acid, total volatile fatty acids, and ratio of acetic acid/propionic acid (P > 0.05). The 4 main ruminal cellulolytic bacteria populations including Fibrobacter succinogenes, Ruminococcus flavefaciens, Butyrivibrio fibrisolvens and Ruminococcus albus did not differ from DCAD treatments (P > 0.05). The DCAD levels did not affect bacterial richness and diversity indicated by the indices Chao, Ace and Simpson and Shannon, respectively (P > 0.05). Both weighted UniFrac and unweighted UniFrac showed no difference in the composition of rumen microbiota for CON, HD and LD (P > 0.05). At the phylum level, Bacteroidetes was the predominant phylum followed by Firmicutes, Synergistetes, Proteobacteria, Spirochaetae, and Tenericutes, and they showed no difference (P > 0.05) in relative abundances except for Firmicutes, which was higher in HD and LD compared to CON (P < 0.05). At the genus level, relative abundance of 11 genera were not affected by DCAD treatments (P > 0.05). Level of DCAD had no effect (P > 0.05) on growth performance including dry matter intake, average net gain, average daily gain, and feed conversion ratio; and nutrients digestibility of crude protein, neutral detergent fiber, acid detergent fiber, and organic matter (P > 0.05). Urine pH in LD was lower than HD and CON (P < 0.05). LD resulted in higher plasma calcium than HD and CON (P < 0.05) but not for other plasma metabolites (P > 0.05).
Conclusion: We conclude that, with regard to the great importance of rumen fermentation, these results suggest that reducing DCAD is unharmful for rumen status and provide the feasibility of feeding a low-DCAD to goats.
Dietary cation–anion difference (DCAD; mmol/kg dry matter [DM]) refers to the difference between the number of millimoles of major cations (Na+ and K+) and major anions (Cl− and S2−) of per kg of DM in the diet [1]. Variations in DCAD enables maintenance of electrolyte balance and osmotic pressure, improving the health status of animals. The initial benefit of feeding a low-DCAD diet was reported by Block [2], who observed that treatment with − 172.3 (mmol/kg DM) DCAD could prevent hypocalcemia compared to a control DCAD of + 448.6 (mmol/kg DM). The DCAD level has been an important parameter in the process of formulating diets in recent years [3–7]. Studies have shown that diets with a lower DCAD could enhance the health and extend the economic life of transition mammary animals [8–10]. Recently, attention was given to the interaction of DCAD with calcium (Ca) concentration[11], cholecalciferol/calcidiol[12], 5-hydroxy-l-tryptophan [13], vitamin D [14–16], and duration [17, 18] of blood Ca level increases. Results of these studies indicated that lower DCAD in association with the above factors was effective in improving homeostasis of peripheral blood Ca. Most recently, Santos et al. [19] and Lean et al. [20] used meta-analyses to determine the effects of varying DCAD on the performance, production, and health of cows.
Ruminant diets usually include up to 70% roughage and so are fiber-rich. As the most important digestive organ, the rumen is the main site for digestion and utilization of nutrient, especially for degradation of cellulose by cellulolytic bacteria. Among rumen measurements, pH is the key parameter determining rumen fermentation status and is often lowered when ruminants are fed high levels of concentrate, resulting in ruminal acidosis [21–23]. Bacteria account for 95% of the total amount of rumen microorganisms [24], and thus are the main factor to modulate digestive and metabolic activity of the rumen[25]. Butyrivibrio fibrisolvens, Fibrobacter succinogenes, Ruminococcus flavefaciens, and Ruminococcus albus are reported the four most important cellulolytic bacteria for digestion and utilization of fiber in the rumen [26–29].
However, to our knowledge, there is less information on the rumen fermentation, cellulolytic bacteria populations and microbiota for goats fed a low-DCAD diet, therefore, the present study was conducted to evaluate the effect of a low-DCAD on the above-mentioned parameters of rumen pH, buffering capability, volatile fatty acids of acid, propionic acid, butyric acid, total volatile fatty acids, and acetic acid/propionic acid profiles; ruminal cellulolytic bacteria populations of Fibrobacter succinogenes, Ruminococcus flavefaciens, Butyrivibrio fibrisolvens and Ruminococcus albus, and microbiota. Growth performance, acid-base balance, plasma calcium level and metabolites were also measured. The results should provide a comprehensive evaluation of the feasibility of feeding a low-DCAD diet to goats.
The animal treatment procedures were approved by the local animal care and use committee. Using a randomized block design, 18 Qianbei miscellaneous goats (a native goat breed in the southwest of China; 30.07 kg initial weight and aged 13 months) were allotted to 3 treatments of six replicates of 1 goat per replicate. Animals were fed one of three diets with different DCAD levels (mmol/kg DM): +350 (HD), +100 (CON), and −150 (LD). Diets consisted of peanut straw (Arachis hypogaea), faba bean straw (Vicia faba) and concentrate which were mixed with NaHCO3 (HD) or NH4Cl (LD) and were pelleted (4 mm diameter), with a concentrate: roughage ratio of 30:70.
Goats were fed in separate metabolic cages. The experiment duration was 45 d including a 30-d adaption period and 15-d trial period. The adaption period was divided into three stages. In the first stage (1–12 d), goats were observed for health condition and treated for parasites and disinfected. During the second stage (13–18 d), goats were allowed to adjust to their respective diets. In the third stage (19–30 d), steady DM intake (DMI) was determined for individual goats. After that, in the trial period (31–45 d), goats were fed the treatment diets strictly according to the established DMI at 09:00 and 18:00. All goats had free access to water during the whole experiment. Ingredients and nutrient levels of diets for goats are shown in Table 1. The actual DCAD levels were measured as +338, +152, and −181, showing slight differences from the designed values of +350, +100, and −150 for HD, CON, and LD, respectively.
Table 1 Ingredients and nutrient compositions of diet for goat
Items |
DCAD |
||
HD |
CON |
LD |
|
Ingredients |
|
|
|
Concentrate (%)1 |
30 |
30 |
30 |
Peanut straw (%) |
50 |
50 |
50 |
Faba bean straw (%) |
20 |
20 |
20 |
NaHCO3 (g/d) |
12 |
-- |
-- |
NH4Cl (g/d) |
-- |
-- |
15 |
Nutrient levels2 |
|
|
|
DM (%) |
93.70 |
93.77 |
93.90 |
CP (% DM) |
13.40 |
13.64 |
14.21 |
NDF (% DM) |
39.40 |
39.64 |
39.82 |
ADF (% DM) |
30.35 |
30.62 |
30.64 |
OM (% DM) |
88.02 |
87.03 |
87.71 |
Ca (% DM) |
0.67 |
0.64 |
0.68 |
P (% DM) |
0.48 |
0.49 |
0.49 |
Na (% DM) |
0.52 |
0.19 |
0.14 |
K (% DM) |
1.08 |
0.98 |
0.87 |
Cl (% DM) |
0.23 |
0.27 |
1.25 |
S (% DM) |
0.16 |
0.17 |
0.18 |
DCAD (mmol/kg DM) |
+338 |
+152 |
-181 |
1 Composition and proportion (%): corn 51.07, soybean meal 24.82, wheat bran 11.19, rapeseed meal 8.07, CaHPO4 0.45, lysine 0.27, methionine 0.29, NaCl 0.77, and premix 3.06.
2 Actually measured values.
Diet samples were collected daily during 19–45 d and stored at −20°C. At the end of the feeding experiment, dietary samples were composited and dried at 65°C and were ground to pass a 1-mm screen for proximate chemical composition determination of DM, crude protein (CP), organic matter (OM), Ca and P [30]; and neutral detergent fiber (NDF) and acid detergent fiber (ADF) [31]. An atomic absorption spectrophotometer (iCE 3000 SERIES, Thermo Fisher Scientific, USA) was used to measure Na and K contents. Silver nitrate titration was used to determine Cl concentration. The S level was determined using the magnesium nitrate method as described by Wang and Beede [32]. The DCAD was calculated using the following equation:
DCAD = Na (%)/0.0023 + K (%)/0.0039 − Cl (%)/0.00355 − S (%)/0.0016
All goats were weighed on d 32 as the initial weight and on 46 d as the final weight. The DMI was recorded daily for each goat calculated by allowance of refusals. The average net gain (ANG) was determined by subtraction of initial weight from final weight. The average daily gain (ADG) was determined by dividing ANG with the trial period (15 d). The feed conversion ratio (FCR) was the ratio of DMI to ADG. Total feces was collected during the last five consecutive days for nutrients digestibility analysis. Feces samples from individual goats were dried in an oven at 65°C for 48 h and ground to pass a 1-mm screen for approximate nutrient measurements: CP and OM [30] and NDF and ADF [31].
Rumen status
At 9:00, 13:00, and 17:00 on d 44, rumen fluid was collected through the esophageal cannula via a vacuum pump (VP30, Labtech Instrument Co. Ltd, Beijing, China) and was measured pH using a pH meter (PHS-3C, Youke Instrument Co. Ltd, Shanghai, China). Sample collected at 17:00 was used to assess rumen buffering capability (BC) as described by Tucker et al. [33], and the operational procedure was titrating a 10-mL aliquot of the sample from the original pH dropping to pH 5 using 1 N HCl.
Rumen fluid (20 mL) was used to detect volatile fatty acids (VFA), and the remaining sample was used for microbial high-throughput sequencing analysis. To determine the VFA, acidified samples were centrifuged at 4°C for 10 min (Thermo Fisher-ST 16R), and the supernatant fraction filtered through a 0.45-µm filter. The 1280 μL of filtrate was mixed with 600 μL of 20% metaphosphoric acid and 120 μL of crotonic acid (internal standard). The VFA concentrations in filtered samples were determined by gas chromatography (GC-2010-plus, Shimadzu, Japan) equipped with a flame ionization detector and a capillary column (SH-Rtx-Wax, Shimadzu), and nitrogen was used as the carrier gas. Rumen microbial high-throughput sequencing analysis was performed by TinyGen Bio-Tech (Shanghai) Co. Ltd.
Ruminal fluid samples were stored at −80°C and the DNA extracted from 200-mg samples using a QIAamp DNA Stool Mini Kit (Qiagen, Hilden, Germany) following the manufacturer’s instructions. The DNA concentration and purity were checked by running samples on 1.0% agarose gels.
The PCR amplification of 16S rRNA genes was performed using general bacterial primers: 515F 5ʹ-GTGCCAGCMGCCGCGGTAA-3ʹ and 926R 5ʹ-CCGTCAATTCMTTTGAGTTT-3ʹ. The primers also contained the Illumina 5ʹ-overhang adapter sequences for two-step amplicon library building, following the manufacturer’s instructions for the overhang sequences. The initial PCR reactions were carried out in 25-µL reaction volumes with 1–2 µL of DNA template, 250 mM dNTPs, 0.25 mM of each primer, 1× reaction buffer, and 0.5 U of Phusion DNA Polymerase (New England Biolabs, USA). The PCR conditions consisted of initial denaturation at 94°C for 2 min, followed by 25 cycles of denaturation at 94°C for 30 s, annealing at 56°C for 30 s, and extension at 72°C for 30 s, with a final extension of 72°C for 5 min. The second step of PCR with dual eight-base barcodes was used for multiplexing. Eight-cycle PCR reactions were used to incorporate two unique barcodes to either end of the 16S amplicons. Cycling conditions consisted of one cycle of 94°C for 3 min, followed by eight cycles of 94°C for 30 s, 56°C for 30 s, and 72°C for 30 s, and a final extension of 72°C for 5 min. Prior to library pooling, the barcoded PCR products were purified using a DNA gel extraction kit (Axygen, China) and quantified using FTC-3000 TM real-time PCR. The libraries were sequenced by 2 × 300 bp paired-end sequencing on the MiSeq platform using MiSeq v3 Reagent Kit (Illumina) at Tiny Gene Bio-Tech (Shanghai) Co. Ltd.
The raw fastq files were demultiplexed based on the barcode. The PE reads for all samples were run through Trimmomatic (version 0.35) to remove low-quality base pairs using parameters SLIDINGWINDOW: 50:20 and MINLEN: 50. Trimmed reads were then further merged using the FLASH program (version 1.2.11) with default parameters. The low-quality contigs were removed based on screen.seqs command using the following filtering parameters, maxambig = 0, minlength = 200, maxlength = 580, and maxhomop = 8. The 16S sequences were analyzed using a combination of software mothur (version 1.33.3), UPARSE (usearch version v8.1.1756, http://drive5.com/uparse/), and R (version 3.2.3). The demultiplexed reads were clustered at 97% sequence identity into operational taxonomic units (OTUs) using the UPARSE pipeline (http://drive5.com/usearch/manual/uparse cmds.html). The OTU representative sequences were used for taxonomic assignment against the Silva 128 database with a confidence score ≥0.6 by the classify.seqs command in mothur. The OTU taxonomies (from phylum to species) were determined based on NCBI. The four members of the rumen cellulolytic bacteria community (B. fibrisolvens, F. succinogenes, R. flavefaciens, and R. albus, % of total bacterial 16S rDNA) were selected from “Species” for statistical analysis.
For precise observation of the effect of DCAD level on acid–base balance in vivo, urine pH was measured once every 3 d during the first and second adaption periods (1–12 and 13–18 d), once every 2 days during the third stage (19–30 d), and daily for the trial period (31–45 d). Urine was immediately dipped with special indicator paper (5.4–7.0, SSSreagent Co. Ltd, Shanghai, China; 6.4–9.0, Fuyang Special Paper Co. Ltd, Hangzhou, China) when goats urinated.
Blood samples (10-mL) of every goat were collected from the jugular vein into heparinized plastic syringes at 45 d. Samples were centrifuged at 805×g for 15 min to harvest plasma. Plasma was subsequently analyzed for levels of Ca, glucose (Glu), urea nitrogen (UN), alanine aminotransferase (ALT), aspartate transaminase (AST), alkaline phosphatase (AKP), total protein (TP), albumin (Alb), superoxide dismutase (SOD), glutathione peroxidase (GSH-Px), malondialdehyde (MDA), and catalase (CAT).
The MIXED module in SAS 9.4 (SAS Institute Inc, Cary, NC, USA) was applied for analysis of experimental data. The DCAD levels (+338, +152, and −181) were designated as fixed effects, and goats as the random effect, then Tukey’s method was adopted to determine differences among means of the three DCAD treatments. The experiment results were expressed as mean ± standard deviation. Statistical significance was defined as P < 0.05.
There were no significant differences in rumen pH for the same sampling intervals (0, 4, and 8 h; P > 0.05; Table 2) and all the collected mean among the HD, CON, and LD treatments (0, 4, and 8 h; P > 0.05). The variation of DCAD had no effect on ruminal BC and levels of acetic acid, propionic acid, butyric acid, total VFA (TVFA), and acetic acid/propionic acid (A/P) in the goats (P > 0.05).
Table 2 Effect of dietary cation–anion difference on the rumen pH of goats
Items |
DCAD |
SEM1 |
P-value |
|||
HD (+338) |
CON +152) |
LD (-181) |
||||
pH |
|
|
|
|
|
|
Sampling time interval (h) |
Sampling time |
|
|
|
|
|
0 |
09:00 |
7.22 ± 0.19 |
7.21 ± 0.41 |
7.17 ± 0.36 |
0.14 |
0.97 |
|
13:00 |
7.10 ± 0.22 |
7.17 ± 0.57 |
6.99 ± 0.49 |
0.23 |
0.86 |
|
17:00 |
7.07 ± 0.25 |
7.20 ± 0.66 |
6.98 ± 0.54 |
0.26 |
0.83 |
4 |
13:00 |
7.05 ± 0.31 |
6.96 ± 0.36 |
7.04 ± 0.34 |
0.18 |
0.93 |
|
17:00 |
7.08 ± 0.38 |
6.96 ± 0.52 |
7.14 ± 0.47 |
0.24 |
0.87 |
8 |
17:00 |
7.15 ± 0.23 |
6.99 ± 0.38 |
7.08 ± 0.46 |
0.19 |
0.86 |
Mean |
7.20 ± 0.22 |
7.18 ± 0.41 |
7.17 ± 0.34 |
0.14 |
0.98 |
|
BC2 (mL/L) |
44.75 ± 0.88 |
43.35 ± 3.65 |
43.10 ± 3.81 |
1.38 |
0.67 |
|
Acetic acid (mmol/L) |
42.22 ± 4.86 |
43.45 ± 5.71 |
44.81 ± 6.36 |
2.32 |
0.74 |
|
Propionic acid (mmol/L) |
13.00 ± 1.96 |
13.07 ± 1.75 |
12.71 ± 2.37 |
0.84 |
0.95 |
|
Butyric acid (mmol/L) |
8.00 ± 1.78 |
8.48 ± 2.28 |
8.11 ± 2.02 |
0.83 |
0.91 |
|
TVFA3 (mmol/L) |
63.21 ± 8.08 |
65.01 ± 8.89 |
65.63 ± 8.10 |
3.41 |
0.87 |
|
A/P4 |
3.27 ± 0.34 |
3.34 ± 0.36 |
3.67 ± 1.13 |
0.29 |
0.60 |
1SEM=standard error of mean. The same as below.
2BC=Buffering capability.
3TVFA = acetic acid + propionic acid + butyric acid.
4A/P = acetic acid/propionic acid.
Rumen cellulolytic bacteria
The relative contents of F. succinogenes, R. flavefaciens, B. fibrisolvens, and R. albus were not significantly affected among goats fed HD, CON, and LD diets (P > 0.05; Table 3). The proportions of F. succinogenes and R. flavefaciens were markedly increased with lower DCAD compared to B. fibrisolvens and R. albus.
Table 3 Effect of dietary cation–anion difference on the rumen cellulolytic bacteria communities of goats
Items |
DCAD |
SEM |
P-value |
||
HD (+338) |
CON (+152) |
LD (-181) |
|||
Fibrobacter succinogenes (%) |
0.443±0.003 (48.63 %) |
0.362±0.003 (41.18 %) |
0.324±0.001 (62.91%) |
0.098 |
0.70 |
Ruminococcus flavefaciens (%) |
0.343±0.21 (37.65 %) |
0.414±0.86 (47.10 %) |
0.119±0.09 (23.11 %) |
0.210 |
0.51 |
Butyrivibrio fibrisolvens (%) |
0.12±0.12 (13.17 %) |
0.085±0.07 (9.67 %) |
0.065±0.03 (12.62 %) |
0.033 |
0.59 |
Ruminococcus albus (%) |
0.005±0.01 (0.55 %) |
0.018±0.03 (2.05 %) |
0.007±0.01 (1.36 %) |
0.008 |
0.50 |
After Illumina Miseq high-throughput sequencing, a total of 698,626 valid reads were obtained with an average length of 410 bp. The VENN graph showed that there were 1075 of total 1261 OTU in 3 groups shared. There were 22, 25, 5 individual OTUs, accounting for 1.89%, 2.33%, and 0.43% for LD, HD, CON, respectively (Fig. 1). Rarefaction curves were established to quantify the OUT coverage of sampling and each rarefaction tended to be gentle with the increase of sequence number, and meanwhile, the OTU rank abundance in the 3 groups exhibited a gentler slope and wider distribution on the horizontal axis (Fig. 2).
Fig. 1 The VENN graph of rumen bacterial of goats fed diets with varying cation-anion difference levels, which displaying that the disposition of operational taxonomic units (OTUs) among the control group (CON, n=6), high dietary cation-anion difference (HD, n=6) and low dietary cation-anion difference (LD, n=6) in rumen fluid microbiota.
Fig. 2 Rarefaction curves (A) and the OTU rank abundance (B) of rumen bacterial of goats fed diets with varying cation-anion difference levels. The CON group: CON1, CON2, CON3, CON4, CON5, CON6. The HD group: HD1, HD2, HD3, HD4, HD5, HD6. The LD group: LD1, LD2, LD3, LD4, LD5, LD6.
Alpha diversity results showed that DCAD levels did not affect Chao, Ace, Simpson, and Shannon as listed in Table 4 (P > 0.05). The Chao, Ace, Shanno and Simpson chart of each group was also tended to be gentle corresponding to the increase of sequence number (Fig. 3).
Table 4 Effect of varying dietary cation-anion difference levels on rumen bacterial community richness and diversity of goats
Items |
Dietary treatments1 |
SEM |
P-value |
|||
HD (+338) |
CON (+152) |
LD (−181) |
||||
Bacterial Richness |
Chao |
906±89 |
843±129 |
811±116 |
4.28 |
0.333 |
Ace2 |
900±82 |
834±117 |
815±96 |
4.33 |
0.356 |
|
Bacterial Diversity |
Simpson |
0.027±0.02 |
0.033±0.01 |
0.028±0.01 |
0.04 |
0.760 |
Shannon |
4.87±0.31 |
4.45±0.21 |
4.69±0.29 |
0.21 |
0.048 |
1Dietary treatments: HD: High dietary cation-anion difference; CON:Control group; LD: Low dietary cation-anion difference.
2Ace=abundance-based coverage estimator.
Fig. 3 Rarefaction curves (A) and the OTU rank abundance (B) of rumen bacterial of goats fed diets with varying cation-anion difference levels. The CON group: CON1, CON2, CON3, CON4, CON5, CON6. The HD group: HD1, HD2, HD3, HD4, HD5, HD6. The LD group: LD1, LD2, LD3, LD4, LD5, LD6.
According to Fig. 4, both weighted UniFrac (axis 1 + axis 2 = 66.7%, Fig. 4a) and unweighted UniFrac (axis 1 + axis 2 = 38.37%, Fig. 4b) were observed no difference in the composition of rumen microbiota for CON, HD and LD (P > 0.05).
Fig. 4 Principle coordinate analysis (PCoA) based on weighted UniFrac and unweighted UniFrac on rumen bacterial of goats fed diets with varying cation-anion difference levels. (a) PCoA based on weighted UniFrac of rumen bacterial of goats fed diets with varying cation-anion difference levels. (b) PCoA based on unweighted UniFrac of rumen bacterial of goats fed diets with varying dietary cation-anion difference levels. The percentage of variation explained by PC1 (first principal component) and PC2 (second principal component) was indicated in the axis. CON = the control group (n = 6); HD= high dietary cation-anion difference (n = 6); LD = low dietary cation-anion difference (n = 6).
Taxonomic classification summary indicated that 16 phyla were tested in all samples (Fig. 5A). At the phylum level, Bacteroidetes (61.60%) was the predominant phylum followed by Firmicutes, Synergistetes, Proteobacteria, Spirochaetae, Tenericutes with average relative abundances of 25.32%, 5.84%, 1.82%, 2.08%, 1.2%, respectively, but there was no difference (P > 0.05) among the groups on the above phylum levels except for Firmicutes, which was significantly higher in HD and LD compared to CON (P=0.008, Table 5).
At the genus level, taxon displayed that the relative abundance of 11 genera were not affected by DCAD among all samples (P > 0.05; Fig. 5B). At the same time, Prevotella, Paraprevotella, Selenomonas, Ruminococcus, Ruminococcus, Butyrivibrio, Quinella, Fretibacterium and Treponema showed no grouping difference of the genera across treatments (P > 0.05). Among the genera with relative abundance exceeded 0.1%, prevotella was the dominant genus in each group with the highest proportion (Table 5).
Fig. 5 Distributions of microbiota at phyla and genus level. (a) Relative abundances of phyla levels are depicted as mean values for the CON, the LD (Low dietary cation-anion difference) group and HD (High dietary cation-anion difference) group. (b) Distributions of genera in rumen fluid of the CON, the LD (Low dietary cation-anion difference) group and HD (High dietary cation-anion difference) group.
Table 5 Effect of varying dietary cation-cation difference on relative abundance (%) of bacteria taxa > 0.1% of average abundance in the rumen fluid of goats
Phylum |
Genus |
Dietary treatments1 |
SEM |
P-value |
||
HD |
CON |
LD |
||||
Bacteroidetes |
|
60.40±6.24 |
65.08±4.41 |
59.33±9.13 |
2.81 |
0.328 |
|
Prevotella |
18.92±3.66 |
24.16±2.72 |
16.09±2.68 |
3.05 |
0.223 |
|
Paraprevotella |
3.58±1.23 |
2.76±0.89 |
3.51±1.52 |
1.24 |
0.877 |
Firmicutes |
|
28.43±4.90a |
18.71±5.30b |
28.81±6.08a |
2.23 |
0.008 |
|
Selenomonas |
1.59±0.97 |
0.56±0.29 |
1.32±0.89 |
0.78 |
0.629 |
|
Ruminococcus |
1.07±0.16 |
1.54±0.74 |
1.29±0.39 |
0.49 |
0.786 |
|
Succiniclasticum |
1.34±0.36 |
0.69±0.09 |
0.68±0.24 |
0.25 |
0.143 |
|
Butyrivibrio |
1.02±0.32 |
0.46±0.13 |
0.48±0.01 |
0.20 |
0.118 |
|
Quinella |
2.29±0.37 |
0.71±0.30 |
3.73±2.15 |
1.27 |
0.277 |
Synergistetes |
|
5.20±1.87 |
6.34±1.91 |
5.99±2.29 |
2.03 |
0.923 |
|
Fretibacterium |
5.18±1.88 |
6.31±1.92 |
5.97±2.30 |
2.04 |
0.923 |
Spirochaetae |
|
1.44±0.25 |
2.13±0.74 |
1.89±0.35 |
0.49 |
0.610 |
|
Treponema |
0.72±0.16 |
0.98±0.52 |
1.39±0.35 |
0.37 |
0.458 |
Proteobacteria |
|
1.10±0.42 |
4.49±2.38 |
0.67±0.07 |
1.40 |
0.140 |
Tenericutes |
|
1.07±0.17 |
1.20±0.16 |
1.33±0.44 |
0.29 |
0.817 |
Levels of DMI were unaffected by DCAD variations (P > 0.05; Table 6). Lower DCAD had no effect (P > 0.05) on gwowth performance of final weight, ANG, ADG, and FCR and digestibility of crude protein, NDF, ADF, and OM for goats.
Table 6 Effect of dietary cation–anion difference on growth performance of goats
Items |
DCAD |
SEM |
P-value |
||
HD (+338) |
CON (+152) |
LD (-181) |
|||
Performance |
|
|
|
|
|
Initial weight (kg) |
30.35 ± 3.4 |
30.43 ± 2.6 |
29.44 ± 3.7 |
1.36 |
0.86 |
DMI (g/d) |
899.0 ± 213.6 |
857.5 ± 120.9 |
864.0 ± 124.3 |
64.85 |
0.89 |
Final weight (kg) |
31.33 ± 3.4 |
31.50 ± 3.8 |
30.41 ± 3.3 |
1.48 |
0.82 |
ANG (kg) |
0.98 ± 0.3 |
1.07 ± 0.4 |
0.97 ± 0.3 |
0.21 |
0.89 |
ADG (g/d) |
65.3 ± 12.2 |
71.3 ± 22.7 |
64.7 ± 18.8 |
8.28 |
0.63 |
FCR |
13.76 ± 6.6 |
12.02 ± 3.5 |
13.36 ± 5.9 |
3.00 |
0.81 |
Digestibility (%) |
|
|
|
|
|
CP (%) |
56.37 ± 2.41 |
58.66 ± 5.91 |
57.68 ± 5.42 |
1.97 |
0.75 |
NDF (%) |
50.26 ± 6.01 |
45.84 ± 8.39 |
45.58 ± 7.65 |
3.12 |
0.50 |
ADF (%) |
51.19 ± 6.84 |
48.35 ± 7.74 |
47.59 ± 8.27 |
3.20 |
0.71 |
OM (%) |
62.34 ± 3.59 |
57.73 ± 5.66 |
58.59 ± 3.92 |
2.11 |
0.36 |
There was no difference (P > 0.05) in urine pH for HD, CON, and LD during the observation period of 1–12 d, with pH values of 8.48, 8.43, and 8.46, respectively (Fig. 6a). Urine pH decreased slightly (8.45, 8.50, and 8.13 for HD, CON, and LD, respectively) during the dietary replacement period (13–18 d). Urine pH decreased significantly (8.43, 8.36, and 7.40 for HD, CON, and LD, respectively) with LD obviously lower than both HD and CON. During the trial period (31–45 d), compared with HD and CON, LD reduced urine pH over HD and CON (8.43, 8.42, and 6.75 for HD, CON, and LD, respectively; P < 0.05). Urine pH values were unaffected by DCAD variation between HD and CON (P > 0.05). Furthermore, urine pH was closely related to DCAD levels within the trial period (31–45 d; R2 = 0.9066, P < 0.05; Fig. 6b).
Fig. 6a Goat urine pH for dietary cation–anion difference levels throughout the experiment
Fig. 6b Association between urine pH and dietary cation-anion difference of goats in the trial period
Feeding of the LD diet resulted in the highest plasma Ca level (Table 7), which was significantly higher than both HD and CON (P < 0.05).There were no significant differences in plasma Glu, UN, ALT, AST, AKP, TP, Alb, GSH-Px, CAT, SOD, and MDA among the DCAD treatments (P > 0.05; Table 7).
Table 7 Effect of dietary cation–anion difference on the plasma metabolites of goats
Items |
DCAD |
SEM1 |
P-value |
||
HD (+338) |
CON (+152) |
LD (−181) |
|||
Ca (mmol/L) |
2.31 ± 0.21 |
2.38 ± 0.17 |
2.91 ± 0.17 |
0.08 |
<0.01 |
Glu (mmol/L) |
4.57 ± 1.06 |
5.07 ± 0.91 |
4.60 ± 0.36 |
0.34 |
0.56 |
UN (mmol/L) |
6.22 ± 0.66 |
5.63 ± 1.09 |
6.24 ± 1.16 |
0.41 |
0.50 |
ALT (IU/L) |
8.61 ± 2.55 |
10.22 ± 3.21 |
10.76 ± 3.57 |
1.70 |
0.85 |
AST (IU/L) |
11.25 ± 3.63 |
13.98 ± 4.77 |
11.91 ± 3.41 |
2.01 |
0.80 |
AKP (King unit/100 mL) |
20.85 ± 6.90 |
21.73 ± 12.57 |
17.81 ± 6.92 |
3.75 |
0.75 |
TP (g/L) |
99.38 ± 14.10 |
112.89 ± 28.72 |
92.62 ± 21.62 |
11.47 |
0.37 |
Alb (g/L) |
36.38 ± 4.34 |
37.11 ± 6.36 |
35.19 ± 6.31 |
2.35 |
0.85 |
GSH-Px (U/mL) |
670.59 ± 127.49 |
755.29±296.72 |
709.41 ± 164.46 |
152.50 |
0.59 |
CAT (U/mL) |
2.35 ± 0.79 |
1.95 ± 0.60 |
2.47 ± 1.03 |
0.66 |
0.84 |
SOD (U/mL) |
68.53 ± 9.27 |
63.59 ± 8.59 |
67.98 ± 6.79 |
3.38 |
0.43 |
MDA (nmol/mL) |
37.14 ± 5.57 |
34.61 ± 7.06 |
36.79 ± 5.94 |
2.54 |
0.75 |
Rumen pH has a key role in measuring rumen status and can reflect the composition and abundance of rumen microflora. Rumen health is maintained by an appropriate pH which is in range of 6.26 ~ 6.79 [34]. Rumen pH decreases when ruminants are fed high levels of concentrate. This is quite different from the acidified rumen status induced by increases in anion (Cl− and S2−) concentrations in the diet and which elevates Cl− and S2− in the urine and so reduces urine pH. In this study, rumen pH was unaffected for all DCAD treatments and ruminal fluid sampling time points. This was also shown in the study of Apper-Bossard [35].
Briggs et al. [36] found that rumen BC was closely related to rumen pH, and Church [37] argued that the rumen buffer system was controlled by pH, pCO2, and VFA. The rumen BC can maintain a stable rumen status by keeping any sudden rise or fall in rumen pH within a certain range [38]. The rumen has a relatively stable buffer system, which is closely related to the feed, saliva, and secretion of the rumen wall and maintains the rumen internal environment. In our study, the lack of significant difference in rumen BC and pH showed that DCAD reduction did not influence the rumen internal environment. However, the rumen BC was numerically elevated with increased DCAD level on cows when adding NaHCO3 [33, 39, 40]. The explanation could be attributed to the use of a different animal species and different dietary ingredients. Generally, up to now there has been insufficient information on rumen BC in goats, and further study is needed due to the importance of BC for rumen status.
The rumen VFA profile is mainly impacted by the proportion of concentrate and forage in the diet. Increasing DCAD, by adding K and Na, had no effect on rumen VFA concentration [41]. Tucker et al. [42] reported that the rumen VFA profile was unaffected by DCAD levels of − 100, 0, 100, and 200. Apper-Bossard et al. [35] found no significant difference in VFA concentration with varying DCAD. These results are consistent with the present study in which rumen VFA concentration was not significantly affected by DCAD level. Correspondingly, the VFA profiles and A/P levels were unaffected by DCAD. This indicates that the rumen fermentation pattern was unaltered by the three DCAD treatments.
In the current study, DCAD variation had no influence on the populations of B. fibrisolvens, F. succinogenes, R. flavefaciens, and R. albus. This result indicates that reducing DCAD would not affect the growth and colonization of rumen cellulolytic bacteria. This is likely associated with stable rumen pH and VFA maintained by the constant ratio of concentrate to roughage (30:70) for the three DCAD diets used in this study. Grilli et al. [43] reported that a high proportion of maize disturbed the ruminal bacterial ecosystem of goats. Wang et al. [44] and Li et al. [45] noted that the relative quantity of ruminal B. fibrisolvens, F. succinogenes, R. flavefaciens, and R. albus was improved with decreasing rumen pH. In our most recent study (unpublished data) with dairy goats as experimental animals fed 4 DCAD levels at 349, 120, and − 167, respectively, rumen BC was not influenced. This result also supports the conclusions in the present study.
The microbiota composition of the gastrointestinal tract (GIT) influence the health of animals as well as productivity [46], the structure and composition of ruminant microbial community are often devoted to evaluating the health status of the host and ensuring the healthy development of the rumen. Furthermore, the diversity and composition of the GIT microbiota can be influenced by many factors including age, diet, feeding management and feed additives [47, 48]. Microbial richness in rumen will be altered with dietary composition [49, 50]. Dietary nutrients are fermented by rumen rumen microbial, such as bacteria, fungi and protozoa, and then degraded into VFA and MCP to provide energy for ruminants [51] and to guarantee the healthy and stable rumen environment [52, 53]. In such reactions, rumen pH plays a decisive role in the composition and abundance of rumen microflora, and is an important indicator reflecting whether the composition and abundance of rumen microflora are normal [54, 55]. Similarly, rumen microbial is also influenced by animal species, diet composition and different ages [56–58].
The study of Zhang et al. [59] showed that rumen bacterial species, Chao index and Ace index were affected by rumen pH which can alter the bacterial community structure. Accordingly, Guo [60] found that decreased rumen pH would up-regulate bacterial diversity, composition and abundance of bacteria. In this experiment, rumen pH was not impacted by DCAD level, which exactly explained why the Chao index, Ace index, Simpson index and Beta diversity of the three groups were homogeneous. The Shannon index of the HD group was the highest, possibly because NaHCO3 was added as a buffer to neutralize gastric acid and was essential for stomach health by creating a suitable internal environment for rumen microorganisms. Reducing DCAD had little effect on Shannon index. In addition, the dominant bacteria of the three groups in this experiment were Bacteroidetes followed by Firmicutes and Synergistetes, which was coincide with the results of previous studies on the basis of phylum [61–63]. In terms of genus level, the relative abundance of Prevotella, one of the primary protein-degrading microorganisms of Qianbei miscellaneous goats, was the highest, which was supported by the results of other ruminant studies [64]. The Firmicutes are a kind of intestinal bacteria related to obesity can degrade insoluble fiber and play an important role in the process of substance metabolism in rumen. The abundance of Firmicutes in HD and LD were 51.95% and 53.98% higher than CON, respectively in the study. The up-regulation in Firmicutes abundance might be due to our application of a high forage diet. In summary, neither the reduction of DCAD have a negative effect on the bacterial diversity nor changed the content of the main bacterial genus.
As described above, rumen pH, BC, VFA, cellulolytic bacteria population and rumen microflora diversity are crucial in the fermentation status of ruminants. The unaffected these parameters in the present study provide further reliable information and the feasibility on feeding a low-DCAD diet to goats.
Level of feed intake is the most important prerequisite for animal growth performance. Generally, pure anionic salt exerts some reduction on feed intake when it was simply mixed into diet, due to its bitter taste and poor palatability [65, 66]. Therefore, improving palatability of anionic salt is important for DMI and growth performance. Our previous study [67] showed that DMI of goats fed a low-DCAD did not decrease because the anionic salts were mixed with molasses and dried distillers grains with solubles. Takagi and Block [68, 69] also observed that reducing DCAD did not impact DMI containing anionic salts when feeding a total mixture ration. Diets were pelleted in our study, and therefore, they were unaffected for goats fed diet HD, CON, and LD. This indicates that DMI is unaffected by anionic salts inclusion as long as the bitter taste is concealed.
The levels of DCAD had no effect on goat final weights. This can be attributed to the similar DMI level and possibly because the goats in this experiment were a local breed, with their adult steady body weight averaging as much as 35 kg, and thus there was limited potential for body weight gain. Accordingly, ADG and FCR did not show difference among goats provided HD, CON, or LD diets.
Nutrients digestibility can be used to measure the digestion and absorption degree of a diet. Correlative analysis in the study demonstrated the digestibility of CP, NDF, ANF, and OM were unaffected for goats fed the three DCAD diets. Apper-Bossard et al. [35] reached similar conclusions. Positive DCAD diet is beneficial to the growth and reproduction of rumen cellulolytic bacteria, and can maintain the activity of digestive enzymes in the digestive tract, thus improving nutrient digestibility [70]. However, different levels of DCAD did not alter the digestibility of CP, NDF, ADF and OM with the rumen pH and rumen cellulolytic bacteria community remained unchanged in this research.
Urine pH is a useful indicator to monitor the effect of a reduced DCAD diet on acid–base balance in goats and sheep [71, 72], dairy cows [8, 73, 74], and buffalo [4]. This phenomenon can be explained using the strong ion difference theory of Stewart [75], who argued that with reductions of DCAD, the concentration of anions in blood would increase and cause the kidney to expel redundant H+ in urine, resulting in lower urine pH. The recommended urine pH is 6.5–6.8, because too low a level would exert a burden on the kidneys [76–78]. Our results showed that urine pH value in goats fed LD was lower than HD and CON. This is accordant with the recommended level, and there is a strong association between DCAD and urine pH in the trial period, suggesting that the LD level is appropriate for the diet of goats.
Muscle contraction, conduction of nervous impulses, and signal transduction are closely dependent on blood Ca homeostasis. Following the study of Block [2] and subsequent results of ruminant researchers [71, 79, 80], reducing the DCAD level has been the most commonly used strategy to increase blood Ca levels in transition mammary animals [81]. Our previous study showed that reducing DCAD could increase the plasma Ca concentration of female goats [67]. In the current experiment, the LD caused higher plasma Ca level than HD and CON by 25.97% and 22.27%, respectively, indicating more stable blood Ca homeostasis. Horst et al. [82] and Goff and Horst [83] claimed this may be because LD-induced acidic status enhanced Ca absorption in the gastro-intestine and also increased Ca resorption in the bone, facilitating Ca matrix flow into blood for easier transfer from lumen to blood.
In the review of Khanal and Nemere [84] who described three steps for Ca absorption. First, Ca2+ influx occurs at the apical membrane via epithelial Ca2+ channels, which include the transient receptor potential vanilloid receptor 6 (TRPV6). This step is considered to be rate-limiting for transcellular Ca2+ transport [85]. Second, intracellular diffusion is facilitated by the vitamin D-dependent calcium binding protein D9k (CaBP-D9k) [86]. Last, extrusion at the basolateral membrane is achieved by either the Na+/Ca2+ exchanger1 (NCX1) or the plasma membrane calcium ATPase 1b (PMCA1b) [87]. The proteins of TRPV6, CaBP-D9k, and PMCA1b are believed to be the three key factors involved in the Ca absorption process, therefore, we speculate that TRPV6, CaBP-D9k, and PMCA1b expression levels in the intestine might be upregulated for animals fed the LD diet. However, a lack of data on TRPV6, CaBP-D9k, and PMCA1b expression in the intestine limits our further discussion. Future work is needed in this area.
Blood measurements are useful to reflect the metabolic status of animals. Our results showed that all plasma levels of Glu, UN, ALT, AST, AKP, TP, Alb, SOD, GSH-Px, MDA, and CAT were unaffected by DCAD variation. This indicates that reducing DCAD had no effect on nutrient metabolite processes in the plasma. This is supported by previous study. Melendez and Poock [9] reported that lowering DCAD had little effect on blood Alb. Wu et al. [67] found that DCAD level (+ 300, + 150, 0, and − 150) had no significant effect on plasma GSH-Px and MDA content in female goats.
Reducing DCAD has little influence on rumen status and rumen microbiota, showing no harmful to rumen fermentation of goats. Blood Ca level is increased and urine pH is decreased by DCAD reduction. These results indicate the feasibility of feeding a low-DCAD diets to goats.
Dietary cation-anion difference; HD:High DCAD; CON:Control; LD:Low DCAD; DM:Dry matter; DMI:Dry matter intake; NDF:Neutral detergent fiber; ADF:Acid detergent fiber; AIA:Acid insoluble ash; VFA:Volatile fatty acid; TVFA:Total volatile fatty acid; A/P:Acetate/Propionate; NH3-N:Ammonia nitrogen; CP:Crude protein; OM:Organic matter; ANG:Average net gain; ADG:Average daily gain; F/G:Feed/Gain; FCR:Feed conversion ratio; BC:Buffer capability; Glu:Glucose; UN:Urea nitrogen; ALT:Alanine aminotransferase; AST:Aspartate transaminase; SOD:Superoxide dismutase; AKP:Alkaline phosphatase; ALB:Albumin; TP:Total protein; GSH-Px:Glutathione peroxidase; MDA:Malondialdehyde; CAT:Catalase; SEM:Standard error of means
All data generated or analyzed during this study are included in this published article.
Ethics approval
All procedures with animals received prior approval from the Animal Care and Use Committee of the Guizhou University and followed the regulations and guidelines for animal care and welfare established by the committee.
Consent for publication
Not applicable.
Competing interests
All authors declare that there are no present or potential conflicts of interest among the authors and other people or organizations that could inappropriately bias their work.
This study is financially funded by National Natural Science Foundation of China (31960674), Science and Technology Project of Guizhou Province ([2018]1403), and the Construction Program of Biology First-class Discipline in Guizhou (GNYL[2017]009).
K.Y., X.T., and W.W. conceived and designed the experiment. S.J., L.S., E.H., F.W., Z.G., and Z.M. carried out the farm and the collection of samples. K.Y. and S.J. gathered and analyzed the data, and both took charge for data interpretation and sequencing data (figures and tables), and K.Y. and X.T. wrote the manuscript. T.H. made feasible suggestions for the manuscript. WW acted as overall study director. All authors have read and approved the final manuscript.
The authors thank the Animal Nutrition Laboratory of Guizhou University for providing the feeding trial site and chemical analysis facilities for this research. The authors also thank International Science Editing for editing this manuscript (http://www.internationalscienceediting.com).