Long-term high-grain diet alters ruminal pH, fermentation, and epithelial transcriptomes, leading to restored mitochondrial oxidative phosphorylation in Japanese Black cattle

To increase intramuscular fat accumulation, Japanese Black beef cattle are commonly fed a high-grain diet from 10 to 30 months of age. Castrated and fistulated cattle (n = 9) were fed a high-concentrate diets during the early, middle, and late stages consecutively (10–14, 15–22, 23–30 months of age, respectively). Ruminal pH was measured continuously, and rumen epithelium and fluid samples were collected on each stage. The 24-h mean ruminal pH during the late stage was significantly lower than that during the early stage. Total volatile fatty acid (VFA) and lactic acid levels during the late stage were significantly lower and higher, respectively, than those during the early and middle stages. In silico analysis of differentially expressed genes showed that “Oxidative Phosphorylation” was the pathway inhibited most between the middle and early stages in tandem with an inhibited upstream regulator (PPARGC1A, also called PGC-1α) but the most activated pathway between the late and middle stages. These results suggest that mitochondrial dysfunction and thereby impaired cell viability due to acidic irritation under the higher VFA concentration restored stable mitochondrial oxidative phosphorylation and cell viability by higher lactic acid levels used as cellular oxidative fuel under a different underlying mechanism in subacute ruminal acidosis.


Results
Daily dietary intake, ruminal pH, total VFA, LPS, and LBP. In this study, we re-analyzed our previously published data regarding rumen fermentation parameters 14 to evaluate the relationships between ruminal pH and candidate gene expression in the RE.
No adverse health condition throughout the study period and effects of ruminal cannulation after surgery were observed for any of the cattle. The daily intake and rumen fermentations were described previously 14 . Briefly, the organic matter basis concentrate intake was significantly (P < 0.05) higher during the Middle stage than during the early and late stages (7.6 vs. 6.0 and 6.1 kg/day, respectively), and the forage diet intake was significantly (P < 0.05) higher during the early stage than the middle and late stages (2.1 vs. 1.1 and 1.0 kg/day, respectively). Nutrient adequacy rates of dry matter (DM) and total digestible nutrients (TDN) were calculated based on the total consumption amounts of concentrate diet and forage.
No change in pH was observed during calibration. The 24 h mean ruminal pH was significantly (P < 0.05) lower during the middle and late stages than during the early stage (Fig. 1). Accordingly, the duration of pH <5.6 was significantly (P < 0.05) longer during the late stage than during the early stage, and the durations of pH <5.8 were significantly (P < 0.05) longer during the middle and late stages than during the early stage. The 1 h mean ruminal pH was significantly (P < 0.05) lower during the late stage than during the early and middle stages, and was significantly (P < 0.05) lower during the middle stage than during the early stage (Fig. S1).
The total VFA concentration was significantly (P < 0.05) lower during the late stage than during the early and middle stages (Fig. 2). The lactic acid concentration was significantly (P < 0.05) higher during the late stage than during the early and middle stages, and was significantly (P < 0.05) lower during the middle stage than during the early stage. The rumen lipopolysaccharide (LPS) activities were significantly (P < 0.05) higher during the middle and late stages than during the early stage. However, the proportions of individual VFAs and peripheral blood LPS-binding protein (LBP) concentration did not differ (P > 0.10) among the fattening stages (Fig. 2). (FDR corrected P < 0.05) were identified when comparing the middle and early (first period), late and middle (second period), and late and early (total period) stages, respectively. There were 873 DEGs (315 upregulated) during the first period, 1,216 DEGs (656 upregulated) during the second period, and 115 DEGs (34 upregulated) during the total period, and they were mapped to molecular functions (FDR corrected P < 0.05, fold change ≥2; Table S3). Twenty-two candidate genes were common to all stages (FDR corrected P < 0.05, fold change ≥2), and all of the upregulated or downregulated DEGs during the first period were reversely downregulated or upregulated during the second and total periods (Table S4). Transporter gene expression of the rumen epithelium (RE) included 22 different soluble carrier (SLC) family members, and they were significantly changed during the first (2 upregulated, 11 downregulated), second (12 upregulated, 7 downregulated), and total (1 upregulated, 2 downregulated) periods (Table 1). The principal component analyses (PCA) plots of the fattening stages were similar within and distinct between stages, particularly in the middle versus late and early stages (Fig. 3).

Discussion
During the fattening period, Japanese Black beef cattle are raised on a high-grain diet between 10 and 30 months of age to increase intramuscular fat accumulation. In this study, samples were collected throughout the entire fattening period to explore the long-term changes in ruminal pH and fermentation as described previously 14 , and their consequences for rumen epithelial transcriptomic dynamics. In cattle, a ruminal pH below the physiological range of 5.8-6.5 6 can induce inflammatory responses 15 , and structural or transcriptomic changes in the RE [8][9][10] . Consistent with a previous study 14 , the 24 h mean ruminal pH decreased gradually as the lengths of time at pH <5.6 and <5.8 increased during the middle and late stages of fattening as a result of a long-term high-grain diet. www.nature.com/scientificreports www.nature.com/scientificreports/ These observations are sufficient to diagnose SARA during the middle and late stages (ruminal pH <5.6 or <5.8 for more than 3 or 5 h per day, respectively 16,17 ;). However, lactic acid is 10 times less protonated than VFA in the rumen (pKa 4.9 vs. 3.9; 2), and its concentration peaked during the late stage. Therefore, it is reasonable to postulate that different mechanisms underlie the low ruminal pH in the late (lactic acid accumulation) versus early and middle stages (total VFA accumulation) as an adaptation of rumen bacterial community to and endurance of a long-term, harsh ruminal environment in Japanese Black beef cattle 14 . However, lower ruminal pH during the late stage did not disrupt the gastrointestinal barrier or induce higher LBP levels in peripheral blood, despite inducing significantly higher LPS levels in the rumen 14 .
Transporter genes (SLC family) are considered to be responsible for pH regulation by rumen epithelial absorption 10 . A low ruminal pH or high-grain diet enhances rumen epithelial absorption via enhanced transporter gene expression, such as monocarboxylate transporter (MCT1; 5), sodium hydrogen exchanger isoform 3 (NHE3 18 ;), and downregulated in adenoma (DRA 19 ;). In the present study, we observed the opposite, i.e., most of the 22 different SLC family genes identified were downregulated (11/13) during the first period and then upregulated (12/19) during the second period. In detail, the expression of SLC9A6 (NHE6) and SLC26A3 (DRA) were downregulated (fold change = −4.08 and −6.71, respectively) during the first period, while they were upregulated (fold change = 4.08 and 10.80, respectively) during the second period. These results suggest that expression of genes encoding rumen epithelial absorption was suppressed by acidotic insult during the first period, whereas it was restored during the second, enhancing VFA absorption or transporter gene expression in the rumen.
Interestingly, different types of SARA in the earlier and later stages had completely opposite consequences on the candidate gene expression (DEGs) and its analysis results (canonical pathway and upstream regulators). Furthermore, PCA plots also suggested that the middle stage was the stage most affected by long-term high-grain diet. LaMonte et al. 20 suggested that acidosis alters cellular metabolism, such as the pentose phosphate pathway and glutaminolysis, to mitigate increased stress due to reactive oxygen species (ROS). A high-starch diet also leads to increased oxidative stress and changes in oxidative phosphorylation to alleviate the toxic effects of ROS 21 . Moreover, VFAs, particularly butyric acid, are catabolized within rumen epithelial cells via oxidative pathways, which generate ROS in the RE 22 . However, mitochondrial superoxide dismutase 2 (SOD2; fold change = −2.29, P = 2.16 × 10 −4 ) was significantly downregulated during the first period, and then it was upregulated (fold change = 2.32, P = 5.24 × 10 −5 ) again simultaneously with SOD1 (fold change = 1.52, P = 1.77 × 10 −3 ) during the second period. These results were supported by the fact that the top upregulated candidate gene during the first period, NADH dehydrogenase subunit 6 (ND6), was associated with the significant canonical pathway "Mitochondrial Dysfunction" (P = 2.69 × 10 −7 ; data not shown). Although excess oxidative stress is inferred from the literatures [22][23][24] , the responses of rumen epithelial cells to oxidative stress might be depressed during the first period in contrast to a short-term (17 days) SARA challenge in dairy cows 21 . These results also imply www.nature.com/scientificreports www.nature.com/scientificreports/ mitochondrial metabolic dysfunction (i.e., handling oxidative stress), suggesting greater rumen epithelial cell vulnerability to cellular damage by acidotic insult 25 .
Mitochondrial proficiency and ROS detoxification are critical for cancer cell viability 26 , and proliferating cancer cells must preserve intracellular ATP and NADPH levels through mitochondrial oxidative phosphorylation and the prevention of excessive ROS accumulation 27 . In silico analyses of the DEGs also suggested that the results were in line with mitochondrial dysfunction during the first period. For example, "Oxidative Phosphorylation", in response to peroxisome proliferator-activated receptor γ coactivator 1α (PGC-1α; 29), was the most inhibited canonical pathway with coordinated inhibition of the upstream regulator PPARGC1A (also called PGC-1α) during the first period. In addition, reduced expression of PGC-1 26,28 and decreased oxidative phosphorylation by coordinated downregulation of PGC-1α-responsive genes 29 are responsible for greater insulin resistance, which might be consistent with the inhibited "Insulin Receptor Signaling" canonical pathway in the present study. Therefore, long-term high-grain diet induced excessive oxidative stress in the RE due to the higher acidity and VFA metabolism, resulting in mitochondrial oxidative phosphorylation dysfunction via the downregulation of ubiquinone oxidoreductase subunits, which are cell survival factors 30 . Activation of upstream regulators such as CD437 (induction of cellular apoptosis and cell cycle arrest in cancer cells 31 ;) and inhibition such as NFE2L2 (regulator for normal mitochondrial biogenesis 32 ;) may indicate regulation of cell proliferation, apoptosis, and mitochondrial biogenesis in the RE. Collectively, the in silico analyses suggest impaired cellular viability and vulnerability to cellular damage due to acidotic irritation during the first period.
During the second period, however, the rumen epithelial transcriptomes restored cellular functions with changes in the underlying mechanism of SARA. The significantly higher lactic acid concentration during the late stage seemed to affect the responses of rumen epithelial cells to acidic irritation differently compared to the early and middle stages. Riemann et al. 33 reported that extracellular acidosis results in a rapid, sustained decrease in intracellular pH (pH i ) and tissue hypoxia, and induces increased transcription factor phosphorylation and tumor cells transcriptional activity 33 . The involvement of lactic acid in cancer has attracted attention only recently, and lactic acid directly contributes to tumor growth and progression (i.e., the lactate anion and protons) as an oxidative fuel and support for angiogenesis and metastasis 34,35 . For example, extracellular lactic acid acidosis causes a lower pH i and higher intracellular lactate concentration, and rescues 4T1 (murine breast cancer) cells from glucose deprivation 36 . In addition, lactic acidosis induces a higher oxidative phosphorylation rate in ATP production in nine different cell lines compared to without lactic acidosis 37 , consistent with the activation of "Oxidative Phosphorylation" during the second period. Although our hypothesis may be preliminary and it is necessary to clarify the positive correlation between the ruminal and intracellular lactic acid concentrations, we postulate that the rumen epithelial transcriptomes were influenced, at least in part, by the significantly increased level of lactic acid in the rumen, resulting in greater cell viability by restoring the mitochondrial oxidative phosphorylation pathway. Our hypothesis may also be supported by the most activated (genistein, tyrosine kinase inhibition that enhances apoptosis and cell cycle arrest in cancer cells 38,39 ;) and inhibited (IgG, malignant factor for colorectal cancer cells 40 ;) upstream regulators during the second period.
Finally, activation of upstream regulators such as PD98059 (a mitogen-activated protein kinases kinase inhibitor 41 ;), sirolimus (immunosuppressant that inhibits cell cycle progression 42 ;), bexarotene (anti-cancer action by inducing T-cell lymphoma cell apoptosis 43 ;), and KDM5B (critical regulator of genome stability and DNA repair 44 ;) and inhibition such as MYC (transcription regulator suppressed by intracellular butyrate 45 ;) throughout the overall study period implies moderate restoration of cellular functions in response to long-term high-grain diet in Japanese Black cattle during the fattening period. However, the effects of gene expression on SARA occurrence or long-term high-grain diet feeding remain unclear, and further studies are required to provide host-side elucidation of host-microbiome interactions.

conclusions
In Japanese Black beef cattle, long-term feeding of a high-grain diet (at 10-30 months of age) induced SARA due to the total VFA production during the early and middle stages and due to lactic acid production during the late stage, suggestive of different mechanisms underlying SARA. Total VFA induced a low ruminal pH, which inhibited the oxidative phosphorylation pathway, and was associated with mitochondrial dysfunction and thereby impaired cell viability of the RE during the earlier fattening stage. By contrast, higher lactic acid levels used as cellular oxidative fuel restored the oxidative phosphorylation pathway and cellular viability during the latter fattening stages. Therefore, the specialized fattening technique applied to Japanese Black beef cattle results in unique changes in rumen fermentation characteristics, influencing epithelial transcriptomes through the interaction with altered rumen fermentation as a response to a long-term high-grain diet.

Methods
The experimental protocol was approved by the Iwate University Laboratory Animal Care and Use Committee (A201720; Morioka, Japan). All animal experiments adhered to the animal experiment policy of Hyogo Prefectural Technology Center for Agriculture, Forestry and Fisheries (Hyogo Prefecture, Japan).
Animals and experimental design. The experimental animals and designs were described previously 14 .
Briefly, nine castrated (at age 5-6 months) and subsequently fistulated (at age 12 months) Japanese Black beef bull cattle were housed with free access to food and water throughout the study (10-30 months of age). The fattening period was subdivided into the early, middle, and late stages (10)(11)(12)(13)(14)(15)(16)(17)(18)(19)(20)(21)(22), and 23-30 months of age, respectively) according to general agreement in Japan 12,14 . Cattle were fed concentrate and forage (rice straw) diets during all three stages, and the amount of concentrate diet was increased gradually throughout the experimental period. The forage-to-concentrate ratio was 26:74, 13:87, and 14:86 during the early, middle, and late stages, respectively. The concentrate diet was composed of barely, steam-flaked corn, wheat bran, and soybean meal and contains 71.2% total digestible nutrient (TDN) and 15.7% crude protein (CP), 72.2% TDN and 13.9% CP, and 72.8% TDN and 12.0% CP during the early, middle, and late stage, respectively. The mean ± SE body weight of the cattle was 335 ± 4.4, 439 ± 7.6, 562 ± 11.6, and 712 ± 18.5 kg on prior to the experiment (10 months of age), and early (14 months of age), middle (21 months of age), and late (29 months of age) fattening stage sampling days, respectively. The forage diet was supplied daily at 0930 and 1530 h in two equal portions, and the concentrate was supplied 1 h after the rice straw to maximize forage diet intake and to prevent excessive consumption of concentrate diet during the early stage. Abnormalities of body condition (body temperature, appetite, hydration, and defecation) were observed daily throughout the study period. The intakes of concentrate and rice straw diets were recorded daily during the final 7 days of the early, middle, and late fattening stages. Table S1 shows the body weight, organic matter intake amount, and chemical composition of the fattening stage diets. Chemical composition of the diets was analyzed according to the official method analysis of the Association of Official Analytical Chemists (AOAC) that registered in the Official Method Feed Analysis of Japan 46 . The adequacy rate of diet was calculated based on the nutrient requirement of Japanese Feeding Standard for Beef cattle 47 . Sampling and measurements. The sample collections and measurements were described previously 14 .
Briefly, ruminal pH was measured continuously every 10 min during the final 7 days of the early, middle, and late fattening stages using a radio transmission system (YCOW-S; DKK-TOA, Yamagata, Japan), as described  Table 2. Canonical pathways activated (z-score ≥2) and inhibited (z-score ≤ −2) generated by Ingenuity Pathway Analysis (IPA) in comparisons of the early and middle stages and late and middle stages. a Values indicate a statistically significant match between up-and down-regulation patterns. (2020) 10:6381 | https://doi.org/10.1038/s41598-020-63471-0 www.nature.com/scientificreports www.nature.com/scientificreports/ previously 48 . A pH sensor was placed in the ventral sac of the rumen through the rumen fistula and calibrated at standard pH values of 4-7, before and after obtaining data in each fattening stage. Rumen fluid from the ventral sac of the rumen (adjacent to the pH sensor) and blood samples from the jugular vein were collected on day 4 of the pH measurements during the three stages, for analyses of total VFA, individual VFAs, lactic acid concentration, LPS activity, and LBP concentration. The fluid samples were immediately filtered through two layers of cheesecloth and stored at −80 °C until use. Blood samples were immediately centrifuged (1,500 × g, 15 min, 4 °C) to separate the plasma and then preserved at −80 °C until analyses.
For the VFA analyses, 1 mL 25% HO 3 P in 3 N H 2 SO 4 was added to 5 mL rumen fluid. Total VFA and individual VFAs (acetic, propionic, and butyric acids) were separated and quantified by gas chromatography (GC-2014; Shimadzu, Kyoto, Japan) using a packed glass column (Thermon-3000; 3%) with a Shimalite TPA 60-80 mesh support (Shinwa Chemical Industries, Kyoto, Japan). For lactic acid analyses, fluid samples were centrifuged at 2,000×g for 15 min at 4 °C, and the concentration in the supernatant was determined using a commercial F-kit (D-lactate/L-lactate) (J.K. International, Tokyo, Japan). To measure rumen LPS activity, rumen fluid samples were centrifuged at 11,000×g for 15 min at 4 °C and assayed using a kinetic Limulus amebocyte lysate assay (Pyrochrome with Glucashield; Seikagaku, Tokyo, Japan) as described previously 49 . The plasma LBP concentration was measured using a commercial kit (HK503; HyCult Biotech, Uden, The Netherlands) as described elsewhere 50 .
Transcriptome analyses of the rumen epithelium. The RE was biopsied from the ventral sac of the rumen at a site adjacent to the pH sensor simultaneously with rumen fluid sampling. A sterile, disposable biopsy punch 8 mm in diameter (Kai Industries, Tokyo, Japan) was used to scratch off the RE. Then the RE samples were washed three times in ice-cold PBS and immediately stored at −80 °C until use. Total RNA was extracted from the RE using TRIzol reagent (Invitrogen) as described previously 10 . The purity of the extracted RNA was increased using an RNeasy RNA Clean-up Kit (QIAGEN, Valencia, CA). Total RNA was quantified using a NanoDrop ND-1000 spectrophotometer (NanoDrop Technologies, Thermo Fisher Scientific, Waltham, MA), and its quality

Microarray and Ingenuity pathway Analyses approaches.
A customized bovine oligonucleotide microarray comprising 15,268 genes (Agilent Technologies) was used for one-color microarray analyses to detect genes expressed in the RE 51 . Fluorescently labeled (cyanine 3) complementary cRNA probes were hybridized to the samples and washed using the Gene Expression Hybridization Kit and Gene Expression Wash Pack Kit (Agilent Technologies). The arrays were scanned using an Agilent Microarray Scanner (Agilent Technologies). Feature Extraction ver. 9.5 (Agilent Technologies) was used to process microarray images, align spots, and create raw numerical total spot intensity data. The microarray data from each sample were imported into GeneSpring 12.0 (Agilent Technologies) for use in the software's normalization algorithm and for detection of candidate genes. Normalization was performed by dividing each measurement of each array by the median of all measurements in that array (per-chip normalization). The entire microarray data set has been deposited at the Gene Expression Omnibus (GEO) database. The GEO accession numbers are as follows: Platform, GPL22091; samples, GSM 3901089 to GSM 3901115; series, GSE133152.
Pathway and network analyses of DEGs were performed using Ingenuity Pathway Analyses (IPA) software (Ingenuity Systems, Redwood City, CA 52 ;). Lists of DEGs identified by GeneSpring 12.0 (Agilent Technologies) that corresponded to raw estimated fold changes ≥2 were uploaded into the software application. The IPA knowledge base was used for the DEG enrichment analyses, and the canonical pathways, top upstream regulators, and statistical calculation were analyzed. Quantitative real-time PCR. Nine genes (both upregulated and downregulated) representing a range of fold changes in the microarray analyses were selected for validation by qPCR. Total RNA samples were treated with TURBO DNase (Applied Biosystems, Foster City, CA, USA) to remove contaminating DNA. Total RNA (800 ng) was converted into first-strand cDNA using the High-Capacity cDNA Reverse Transcription Kit (Applied Biosystems). The qPCR was performed using iQ SYBR Green Supermix (Bio-Rad, Hercules, CA) and the StepOne ™ Plus Real-Time PCR system (Applied Biosystems) as described elsewhere 10 . Primers (Table S2) were designed using Primer Express software 3.0 (Applied Biosystems). Each sample contained 2 μl cDNA, 2× SYBR green, and 0.6 μM each primer in a final volume of 20 μl. Amplification conditions were as follows: 95 °C for 10 min followed by 40 cycles of 15 s at 95 °C and 1 min at 60 °C, with collection of fluorescence signal at the end of each cycle. For melting curve data, the temperature was increased from 65 to 95 °C in 0.5 °C increments 10 . The results were recorded as relative changes in gene expression normalized to glyceraldehyde-3-phosphate dehydrogenase (GAPDH), ribosomal protein L27 (RPL27), and β-actin (ACTB) by the 2 −△△Ct method 53 . We examined GAPDH, RPL27, and ACTB to evaluate their use as reference genes for qPCR and there were no significant differences among them. All qPCR experiments were performed in accordance to the guidelines for minimum information for publication of quantitative real-time PCR experiments 5354 .
Statistical analyses. The pH data were summarized 1 day before the sample collection to minimize the influence of opening the cannula stopper on the ruminal pH. The normality of the distributions of variables was assessed using the Shapiro-Wilk test. Significant differences in ruminal pH, duration of time where pH was  Table 3. Genes (raw estimated fold change ≥2, FDR corrected P < 0.05) integrating the oxidative phosphorylation canonical pathway generated by Ingenuity Pathway Analysis (IPA) in comparisons of the early and middle stages and late and middle stages. Figure 5. Heatmap of upstream regulators (P < 0.05) generated by IPA software (QIAGEN Inc., https://www. qiagenbioinformatics.com/products/ingenuitypathway-analysis) in comparisons of the early, middle, and late stages in Japanese Black beef cattle. Cell colors are based on the activation z-score. Blue represents a negative z-score and red represents a positive z-score. N/A, not applicable.