- Research article
- Open Access
Transcriptome analysis of ruminal epithelia revealed potential regulatory mechanisms involved in host adaptation to gradual high fermentable dietary transition in beef cattle
BMC Genomics volume 18, Article number: 976 (2017)
The transition from a high forage to a highly fermentable diet can induce digestive disorders in the rumen. To date, the host mechanisms that regulate the adaption to such dietary transition are largely unknown. To understand the molecular mechanisms involved in such phenomena, RNA-sequencing was performed to identify the changes in the transcriptome of ruminal epithelia during gradual transition from a diet containing 0% to 89% grain.
In total, the expression of 11,044, 11,322 and 11,282 genes were detected in ruminal epithelia of beef heifers (n = 15) fed 0%, 72% and 89% barley grain diet, respectively. The transcriptome profiles of rumen epithelia differed between low grain diet (LGD) (0% grain) and high grain diet (HGD) (72% and 89%), and HGD tended to reduce the expression of genes involved in epithelial catalytic and binding activities. When diet was changed from 72% to 89% grain, the mean ruminal pH change was significantly different among individual heifers with five of them decreased (down group (DG); from 6.30±0.09 to 5.87±0.15, P < 0.01) and five of them increased (up group (UG); from 5.84±0.42 to 6.35±0.37, P < 0.05). The functional analysis of differentially expressed (DE) genes revealed inhibited “Immune response of leukocytes”, “Attraction of phagocytes”, and “Cell movement of leukocytes” (P < 0.05) functions (Z-score = −2.2, −2.2 and −2.0, respectively) in DG, and inhibited “Concentration of lipid” and “Proliferation of epithelial cells” functions in UG (Z-score = −2.0, and −1.8, respectively). In addition, the expression of genes involved in ketogenesis (HMGCL) and lipid synthesis (SREBF2, FABP4) was increased in DG, while the expression of ketogenesis (ACAT2, HMGCS) and cholesterol synthesis related genes (HMGC and FDPS) were deceased in UG. Furthermore, the upstream regulators were found to be involved in the regulation of immune response and cell cycle progress, and SNP (g.46834311A > G) in FABP4 was identified between two groups of animals (P < 0.1).
The identified genes, upstream regulators, and SNP could be potential genetic markers that may account for the varied individual ruminal pH responses to the dietary transition stress.
To meet the global demand for meat consumption, it has become a common practice in the beef industry to use intensive feeding strategies such as high energy and high concentrate diets to finish cattle . However, feeding high concentrate diets has been reported to be associated with digestive disorders that could lead to ruminal acidosis , laminitis , liver abscesses , and hindgut acidosis  in cattle.
When cattle are fed high fermentable diets, the increase and accumulation in short chain fatty acid (SCFA) production together with subsequent decrease in ruminal pH are usually observed due to increased microbial fermentation . In the meantime, ruminal epithelium plays an important role in response to the ruminal pH change and regulation via SCFA absorption [7,8,9]. Therefore, highly fermentable diets could partly elicit an adaptive response by rumen epithelium to maintain the normal function of rumen. Indeed, when a gradual transition strategy is applied, ruminants have strong ability to adapt to highly fermentable diet [2, 10]. However, the individual variation in the adaptation to high concentrate diet has been widely observed in both beef cattle [2, 11] and dairy cows [12, 13]. It is suggested that the differences in the rate and ways of SCFA absorption by the animal may explain the individual variation in the severity of subacute ruminal acidosis in sheep . Recently, Schlau et al.  reported that the expression of sodium hydrogen exchanger isoform 3 (NHE3) gene was different between acidosis-resistant and acidosis-susceptible steers during rapid high grain diet transition, indicating the differences in intracellular proton removal could be attributed to variation in the host response.
To date, the underlying regulatory mechanisms for the host response to ruminal pH change has not been well defined. In this study, we aimed to identify the molecular mechanisms for variation in the response to a gradual transition to a high grain diet within the same animal by characterizing the global gene expression pattern of ruminal epithelia using RNA-seq based transcriptome profiling. The identified mechanisms may help to explain the observed animal variation in maintaining a balanced ruminal pH, which provide new insight into decreasing risk of ruminal acidosis for the beef industry.
Rumen epithelial transcriptomes fed 0%, 72% and 89% grain diets
A total of 1130 million (25.11 ± 2.89 million reads per sample) high-quality 100-bp paired-end reads were obtained from 45 ruminal papillae samples collected from 15 heifers fed diets containing 0%, 72%, and 89% grain, respectively. Of these reads, ~86.1% of them were mapped to the bovine genome (UMD 3.1) and the expression of 11,044, 11,322 and 11,282 genes were detected (with reads per million (RPM) > 1 in 15 heifers fed each diet) in the ruminal epithelial tissue under each dietary condition, respectively. Among them, the expression of 10,880 genes was commonly detected from three dietary conditions (Fig. 1a). The most relevant gene ontology (GO) terms of these commonly expressed genes were “catalytic activity” and “bindings”, followed by “nucleic acid binding transcription factor activity”, “structural molecular activity”, “enzyme regulator activity”, and “receptor activity” (Fig. 1b). When the transcriptome profiles were further compared among three diets, the dietary-dependent expression of genes was detected, with 43, 119 and 102 genes expressed only in 0%, 72% and 89% grain, respectively (Fig. 1a). Further analysis on these diet-dependent genes using DAVID (Database for Annotation, Visualization and Integrated Discovery) showed different functional annotation clusters (Additional file 1 worksheet 2). For example, the enriched gene ontologies for 43 genes under 0% grain diet were “immune response” and “inflammatory response”, and for the 119 genes under 72% grain diet 102 genes under 89% grain, their enriched functions were “extracellular matrix organization”, “protein catabolic process”, respectively. Although the most of these dietary dependent genes were expressed at low level (Fig. 1c) with RPM < 5, the expression of 13 genes showed relatively high abundance (RPM > 5) when heifers were fed diets containing 0% (4 genes), 72% (3 genes), and 89% grain (6 genes), respectively (Table 1).
The differentially expressed genes in the ruminal epithelium between low and high grain diets
The PCA plot revealed that the rumen epithelial transcriptome profiles of heifers fed the low grain (0%) diet were different than those fed high grain diets (72% and 89%) except one outlier from 89% diet (Fig. 2a). When the expression of genes was further compared, 562 genes were differentially expressed (DE). Among them, 432 genes showed higher expression (the fold change (FC) in expression >1.5 with the false discovery rate (FDR) < 0.05) in the ruminal epithelia of heifers fed the high grain diet compared to those fed the low grain diet (Fig. 2b, Additional file 1 worksheet 1) with 130 from 72% grain diet, 191 from 89% grain diet and 111 from both high grain diets (Fig. 2b). These genes were defined as high grain diet enriched (HGD) DE genes. On the other hand, 104 genes had higher expression (the FC in expression >1.5 with the FDR < 0.05) in heifers fed the 0% grain diet comparing to those fed high grain diets with 11 compared with the those fed 72% grain diet, 57 compared with the those fed 89% grain diet, and 36 from both high grain diets (Fig. 2c; Additional file 1). These genes were defined as low grain diet enriched (LGD) DE genes. Although, the functional classification showed no difference in terms of the major functions between the HGD and LGD DE genes (Fig. 2d), the highly expressed HGD DE genes were mainly involved in “Paxillin signaling” and “Integrin signaling” pathways, while highly expressed LGD DE genes were enriched to the canonical pathway of “Complement system” and “Interferon signaling”, based on DAVID and Ingenuity Pathway Analysis (IPA) functional analysis (Additional file1 worksheets 2 and 3).
Variation in responses of heifers during dietary transition from 72% to 89% grain diets
When the ruminal pH change for each individual animal was further compared using mean ruminal pH from a previous study , mean ruminal pH became lower in all heifers after the first dietary transition, from 0% to 72% grain. However, after second dietary transition, from 72% to 89% grain, the mean ruminal pH of heifers had three patterns: lower, from 6.30±0.09 to 5.87±0.15, P < 0.01 (down group (DG), 5 heifers); similar, from 6.17±0.25 to 6.19±0.25, P = 0.27 (balanced group (BG), 5 heifers); and higher, from 5.84±0.42 to 6.35±0.37, P < 0.05 (up group (UG), 5 heifers) (Fig. 3a). When the acidosis index values (pH·min/kg) were further compared, they were increased in DG animals and decreased in UG animals when the diet transitioned from 72% to 89% grain (Fig. 3b). The DG and UG heifers were then selected for further transcriptome comparison to identify the genes that may be associated with the potential molecular mechanisms behind such varied responses.
The differentially expressed genes associated with varied responses to ruminal pH change during the second dietary transition
When the transcriptomes were further compared within DG and UG heifers during the second dietary transition period (89% vs. 72%), the gene expression profiles displayed different changes. In total, 67 DE genes (22 up-regulated and 45 down-regulated; FDR < 0.05, FC > 1.5 or < −1.5) were found in DG (Additional file 2 worksheet 1), and 285 DE genes (122 up-regulated and 163 down-regulated; FDR < 0.05, FC > 1.5 or < −1.5) were found in UG (Additional file 2 worksheet 2). Among them, 23 DE genes were commonly detected from both DG and UG, while 44 and 262 DE genes were only found in DG and UG, respectively (Fig. 3c). Most of the common DE genes (n = 18) showed the same change trend in both DG and UG (14 genes were down regulated and 4 gene were up regulated), while expression of 5 DE genes exhibited opposite change trend between the two groups. In particular, the expression of CYP1A1 and CYP1B1 were up regulated in UG (FC = 1.8 and 1.6) but were down regulated in DG (FC = −1.7 and −2.0), whereas the expression of CHAC1, CDSN and KLK10 were up regulated in DG (FC = 2.3, 2.1 and 1.5) but were down regulated in UG (FC = −2.0, −1.4 and −3.3) (Fig. 3d).
Functional analysis of DE genes detected for DG and UG during the second dietary transition
The functional prediction of 67 DE genes in the DG using IPA revealed that their most relevant functions were “Immune response of leukocytes”, “Attraction of phagocytes”, and “Cell movement of leukocytes” (P < 0.05), and were inhibited when fed 89% grain diet (Z-score = −2.2, −2.2 and −2.0, respectively) (Table 2). In addition, the function of “Quantity of Ca2+”, “Fatty acid metabolism”, and “Mobilization of Ca2+” were also inhibited (Z-score < −1.5) for DG after second dietary transition (Table 2). The predicted function of the down-regulated DE genes (using DAVID n = 45) in the DG was related to “Innate immune response” (P < 0.0001) (Table 3).
For 285 DE genes in the UG, functions of “Concentration of lipid” and “Proliferation of epithelial cells” were inhibited (Z-score = −2.0, and −1.8, respectively), while the functions of “Cell cycle progression” and “Transmigration of phagocytes” were activated (Z-score = 2.2 and 2.0, respectively) after second dietary transition (Table 2). The function prediction of down-regulated genes (n = 163) in UG using DAVID revealed functions related to activity in cell proliferation, such as “Cell division”, “Mitotic metaphase plate congression”, “Chromosome segregation”, and “mitotic nuclear division” (P < 0.001) (Table 3). Moreover, the enriched function of up-regulated genes (n = 122) in UG were “Positive regulation of JUN kinase activity”, “Toxin metabolic process”, “Cholesterol homeostasis”, and “Drug metabolic process” (P < 0.05) (Table 3).
Further functional pathways analysis using IPA showed eight DE genes in DG were found to be involved in two different pathways including “Role of pattern recognition receptors in recognition of bacteria and viruses” (Z-score = −2.0) and “TREM1 signaling” (Z-score = −2.0) (Table 4). For the UG, the “PTEN signaling” pathway was inhibited (Z-score = −2.0), while the pathways of “NF-κB signaling”, “PDGF signaling”, and “Cell cycle: G2/M DNA damage checkpoint regulation” were activated (Z-score = 2.4, 2.2 and 2.2, respectively). The KEGG pathways of those DE genes enriched by DAVID showed similar results (Table 5). Briefly, eight DE genes in DG were found to be involved in three KEGG pathways including “Staphylococcus aureus infection” and “complement and coagulation cascades” (5 down regulated genes), and “T cell receptor signaling pathway” (3 upregulated genes) (Table 5). In the UG, 25 DE genes were enriched in ten different KEGG pathways. Among them, the downregulated genes were involved in four pathways including “Cell cycle”, “Glycolysis/Gluconeogenesis”, “FoxO signaling pathway”, “p53 signaling pathway”, and “Mineral absorption”; the upregulated genes were mainly related to metabolic pathways including “tryptophan metabolism”, “Retinol metabolism”, “Drug metabolism – cytochrome P450”, and “metabolism of xenobiotics by cytochrome P450” (Table 5).
Upstream regulator and network analysis of DE genes detected for DG and UG during the second dietary transition
Four upstream regulators were identified in regulating the DE genes in DG, including Interferon gamma (IFNγ), transmembrane protein 173 (TMEM173), toll like receptor 3 (TLR3), and tumor necrosis factor (TNF) (Fig. 4a). For the DE genes in UG, nine upstream regulators were identified with colony stimulating factor 2 (CSF2) and prostaglandin E receptor 2 (PTGER2) being the hub nodes (Fig. 4b). Two and three networks were enriched (score > 20) in DG and UG, and the related functions were immune response and cell cycle, respectively (Additional file 2 worksheets 3 and 4).
Differential expression of genes involved in lipid transport, metabolism and intracellular homeostasis regulation in DG and UG
Among above genes identified in DG and UG animals, the expression of 14 lipid transport and 13 fatty acid metabolism related genes (RPM > 1 in 5 heifers of at least one diet) were further investigated in the rumen epithelia (Fig. 5). In DG heifers, the expression of SMCT1 (P < 0.1) was lower while the expression of FABP4 (P < 0.05), HMGCL (P < 0.1), and SREBF2 (P < 0.1) was higher when fed the 89% grain diet than the 72% grain diet. In the UG heifers, the expression of FABP5 (P < 0.05), ABCA2 (P < 0.05), and ABCA7 (P < 0.1), ACAT2 (P < 0.05), HMGCS1 (P < 0.1), HMGCR (P < 0.1), and FDPS (P < 0.1) was lower, whereas the expression of ABCA1 (P < 0.05), and ABCA5 (P < 0.1) was higher when fed the 89% grain diet than the 72% grain diet (Fig. 5a & b). In addition, the expression pattern of genes associated with ion transportation in rumen epithelium was further investigated. The expression of NHE8 (P < 0.1), and MCT4 (P < 0.05) was decreased in DG, while the expression of NHE3 (P < 0.05) was increased in UG after the second dietary transition (Additional file 3: Figure S1). Further validation of the differential expression using reverse transcription quantitative real-time PCR (RT-qPCR) showed the expression of FABP4, ABCA1, NHE3 and ACAT2 was consistent as detected with the RNA-seq data (Additional file 4: Figure S2).
In addition, 70 SNPs were detected in 9 genes that relate to lipid transport, metabolism and intracellular homeostasis regulation (Additional file 2 worksheet 5), and association analysis (Fisher’s exact test) of detected single nucleotide polymorphisms (SNPs) indicated that one intronic SNP (g.46834311A > G) in FABP4 was associated with the varied expression in the ruminal epithelial tissue between UG and DG animals (P < 0.1, Additional file 5: Figure S3).
The comparative analysis of the whole transcriptome profiles revealed the high grain diet can alter the gene expression in the rumen epithelia. Although gene expression of rumen epithelial tissues in response to high concentrate diets have been studied, most of the previous studies have focused on selected genes using qPCR/microarray [15,16,17] and/or compared the difference using different groups of animals [11, 18] as well as to study the dairy cows with small numbers of animals . This is the first study that applied RNA-seq based genome-wide transcriptome analysis to study the global gene expression changes in ruminal epithelia under three dietary conditions for the same individual using 15 beef heifers. Furthermore, our effort is the first to identify the varied expression pattern changes in the rumen epithelia of animals having different pH change patterns with the gradual high grain transition from 72% to 89% grain.
It has been known that feeding highly fermentable diets could be associated with the prevalence of ruminal acidosis [6, 19]. When ruminal pH duration time, a widely used diagnostic parameter for the acidosis [20, 21], was analyzed, 8 out of 15 heifers had pH below 5.8 for longer than 5.4 h/d when fed 89% grain, suggesting they may have developed subacute acidosis. In contrast to these heifers, the others had lower ruminal pH duration time when fed the same diet, emphasizing that response variation existed among individuals as previously reported [11, 12], which was reflected by our observation of DG and UG heifers. Together with the mean of continuous ruminal pH, it suggests that each heifer may have developed a different adaptive mechanism, with DG responding less favorably than UG. To understand the mechanism behind such individual variation, the following discussion will be mainly focused on the gene expression pattern variation between in DG and in UG animals.
The most enriched networks of DE genes in DG and UG were related to immune response and cell cycle, respectively, suggesting that the innate immune function might be differently regulated between DG and UG heifers during second dietary transition. Although scarce information is available about the organization of the mucosal immune system in the rumen epithelia, the Langerhans cells were identified in epithelium of bovine forestomach  and bovine forestomachs could receive, elaborate and produce signals and mediators of the innate immune response . Here, we found the innate immune response related functions such as “Immune response of leukocytes”, “Attraction of phagocytes”, and “Cell movement of leukocytes” were inhibited in DG but not in UG, suggesting that innate immune response could be one of the mechanisms involved in epithelial adaptation to high grain diet. In addition, the most abundant GO categories enriched by downregulated DE genes in DG were related to immune response. Decreased expression of immune-related genes (including IL-6, IL-10, and interferon) had been reported in ruminal epithelia of lambs when fed concentrate starter . Similar results were found in DG, in particular, the expression of genes encoding complement component 2/3 (C2 / C3), C-X-C chemokine ligand 5 (CXCL5), and NOD-like receptor C5 (NLRC5) were downregulated during the diet transition from 72% to 89% grain. The C2/C3 plays central role in the complement system and contributes to innate immunity . It’s well known that epithelial cells recognized microbial components through means of Toll-like and NOD-like receptors (TLRs and NLRs) [26, 27]. The decreased expression of NLRC5 might limit the recognition of commensal bacteria, which is essential for the development and function of the immune system in the mucosal and peripheral districts . Indeed, the epimural bacterial population differed between acidosis tolerant and susceptible beef cattle and the expression of TLR2 and TLR4 were lower in susceptible cattle . The previous study also showed varied microbial changes associated with the ruminal epithelia for the same animals . The observed potential inhibited innate immune responses further support our previous speculation that host-microbial interactions could play a role in affecting the host adaptation to the high grain diet transition. The lowered innate immune function may negatively affect the rumen function which result in the persistent decrease of ruminal pH during the second transition in DG.
In addition, the enriched pathway of “G2/M DNA damage checkpoint regulation” was activated through downregulation of genes related to cell cycle progression in UG [30,31,32], including cell division cycle 25 homolog B (CDC25B), cyclin-dependent kinases regulatory subunit 2 (CKS2), cyclin B1/2 (CCNB1/2), polo like kinase 1 (PLK1), and aurora kinase A (AURKA), suggesting the enhanced cell cycle arrest in UG. Furthermore, the enriched KEGG pathway of “p53 signaling” in UG also indicates the higher cell arrest since the activated p53 could activate cell repair and apoptosis procedure to induce cell cycle arrest [33, 34]. Moreover, the pathway of “Metabolism of xenobiotics by cytochrome P450” was enriched in UG heifers. The higher expression of three genes involved in this pathway (alcohol dehydrogenase 6 (ADH6), cytochrome p450 family 1 subfamily A member 1(CYP1A1), and subfamily B member 1 (CYP1B1)) in UG suggest the higher capacity to remove the ruminal toxins in these animals since these genes have been reported to be involved in the metabolism of the some ruminal toxins, such as ethanol and serials of xenobiotics [35, 36]. Taking these together, we speculate that the UG animals may have more activated innate immune responses, cell repair function and toxin removal activities to maintain the ruminal health during the second dietary transition. The upstream regulator analysis further supports our speculation since the identified IFNG, TNF, TLR3, and TMEM173 regulate genes involved in the immune function, while the CSF2 and PTGER2 regulate genes involved in cell proliferation.
Considering the fact that VFAs absorption accounts for up to 50% of the ruminal buffering capacity , the increase in intra-epithelial uptake and metabolism of SCFAs could promote the uptake of SCFAs and a stabilization of ruminal pH . We hypothesized that the expression patterns of genes involved in SCFAs absorption and metabolism may also play an important role in individualized animal rumen pH changes when diet was switched from 72% grain to 89% grain. The expression of sodium dependent monocarboxylate transporters 1 (SMCT1) tended to decrease in the DG heifers but no difference was observed in the UG heifers during the second transition period. The SMCT1 and SMCT2 play important roles in SCFA transport in addition to the SCFA−/HCO3− exchanger system [37, 38]. These results suggest that the absorption of SCFA via the SMCT1 in the DG heifers might be decreased, which may partly account for the lower mean pH in DG heifers. In addition to SCFA absorption, SCFA metabolism could also impact on the rumen pH environment. In this study, we focused on the genes involved in butyrate catabolism in the rumen epithelia, which occurs via ketogenesis and beta-oxidation to produce ketone bodies [39, 40]. The expression of ACATs and HMGCS, the genes encoding the rate-limiting enzymes for ketogenesis , decreased in UG after the second transition period. However, the expression of HMGCL increased in DG for the same dietary transition period, suggesting that UG heifers have decreased ketogenesis in the epithelia comparing to DG.
In addition, the expression of genes involved in fatty acids and lipids transport and metabolism were also found to be different between DG and UG heifers. For example, the increased expression of FABP4 was observed in DG heifers but not in UG heifers when diet transitioned from 72% to 89% grain. Fatty acid binding proteins (FABPs) usually bind to fatty acids and/or lipids to transport and deliver them to different sites for utilization . It has been reported that increased expression of FABP4 was responsible for the enhanced macrophage lipid accumulation in diabetic patients of human, and this effect was further paralleled with elevated intracellular total cholesterol and triacylglycerol levels . Therefore, the increased expression of FABP4 in the ruminal epithelia of DG heifers suggest the intracellular lipid accumulation, which could be partially supported by the tended increased expression of sterol regulatory element binding protein 2 (SREBF2) in the same group of animals. The SREBF family plays important role in regulating cholesterogenic and lipogenic gene expression in the bovine liver and mammary glands [44, 45], while the isoform SREBF2 preferentially activates cholesterol biosynthesis . Although cholesterol synthesis pathway has not been well characterized in the ruminal epithelia, it’s reported that butyrate can be used as a precursor for cholesterol biosynthesis [46, 47]. Compared with DG, the expression of cholesterol synthesis related genes (HMGCS1, HMGCR and FDPS) was decreased in UG heifers when diet transitioned from 72% to 89% grain. Steele et al.  also reported the coordinated downregulation of these genes in the rumen epithelia of dairy cows when fed a 65% grain diet from the first to third week. These suggest that the intracellular cholesterol synthesis could also be one of the mechanisms for the individualized responses when adapted to the high grain diet. It is known that cholesterol is the key component of the plasma membrane and plays crucial role in cellular signal transduction, cell growth, cell polarization, migration, and survival . The cellular content of cholesterol is determined not only by the synthesis but also by the efficiency of influx and efflux  transportation of them and their derivatives through ABC transporters [49, 50]. Among them, ABCA1 has been most intensively studied, which could regulate and maintain cellular cholesterol homeostasis through transferring cholesterol to high-density lipoprotein in various types of cells [51, 52]. In this study, the increased expression of ABCA1 was observed in the UG, but not in DG heifers. In addition, the function of “Concentration of lipid” was inhibited in the UG heifers, suggesting the UG might be more efficient in exporting cholesterol out of rumen epithelia and have less stored cholesterol. Taken together, our results suggest that the deposition of cholesterol in addition to cholesterol synthesis in the rumen epithelia is also one of the mechanisms to attribute to the variation observed between DG and UG animals (Additional file 6: Figure S4). Moreover, previous studies showed that increased intracellular cholesterol could inhibit various ion channels in the biological membrane of human [53, 54], and high cholesterol is closely related to membrane permeability and inflammation in human . Therefore, the inhibited function of “Quantity of Ca2+” and innate immunity in the DG animals could be resulted from higher concentration intracellular content of cholesterol during the dietary adaptation from 72% to 89% grain. Further measuring the cholesterol concentration in rumen epithelia is needed to verify the findings obtained from our transcriptome analysis.
Furthermore, this study also identified the genetic variation of genes related to lipid transport and metabolism. Association analysis indicated that 1 SNP (g.46834311A > G) in FABP4 was associated with the varied ruminal pH response between DG and UG heifers, which may be attributed to the changes in the expression of this gene between DG and UG animal after the secondary diet transition. This suggests that the genetic background between “Up” and “Down” heifers could be one of the factors that lead to the observed varied gene expression. To determine the role of genetic variation in regulating the genes involved in different responses to high fermentable dietary transition, future researches on genotyping of large population with varied phenotypic measures (ruminal pH, duration and area of ruminal pH below 5.8) using genome wide association study (GWAS) are needed.
This is the first study that has identified varied molecular mechanisms may contribute to the individual variations in response to the high grain adaptation through analyzing the genome wide transcriptome in the ruminal epithelia of beef heifers. Understanding the whole transcriptome is essential to reveal the molecular events within cells, and to elucidate the mechanisms regulating the adaptive function and physiology of rumen in response to dietary changes. Overall, we suggest that the different genes that controlled immune function, cellular repair function, and intracellular homeostasis (cholesterol) might be the molecular mechanism accounting for individual variation in the response to gradual transition to a high grain diet. As summarized in Fig. 6, feeding high fermentable diet usually results in increased concentration of ruminal VFA and toxins, which may induce the cell damage for both DG and UG. At transcriptomic level, the innate immunity, cell cycle, toxin metabolism, and cholesterol homeostasis was differentially regulated in DG and UG. In the DG, the inhibited innate immune response and accumulated cholesterol in the cells may aggravate the cell damage and finally decrease the cellular stability and homeostasis. However, in the UG, cell cycle arrest and xenobiotic metabolism was activated for repairing and protecting the cell damage, which contributes to increase cellular stability and homeostasis. It is noticeable that, gene expression can be regulated at both transcriptional and post-transcriptional level, and even post-translational modifications can also affect the activity of protein [56, 57]. Future studies on the protein, enzyme and metabolites are needed to comprehensively determine the molecular mechanisms of individualized host response to high grain feeding. Regardless, the transcriptomic information from this study provides the informative clue for further study such as the identified gene networks, upstream regulators (IFNG, TNF, TLR3, CSF2, and PTGER2), and FABP4-A/G SNP could be potential gene and genetic markers for selecting cattle with maintained ruminal pH during high fermentable diet transition.
Experimental design and sample collection
Twenty-four Angus-Hereford cross-bred yearling heifers (about 8 months old, weighing 244 kg to 369 kg) from the University of Alberta Kinsella Ranch (Kinsella, AB) were used for this study. Detailed information on the animal study, diets, feeding of diets, and intake has been reported previously . Briefly, heifers were blocked into 1 of 6 blocks based on body weight and, within block, were randomly assigned to 1 of 2 treatments. Each block consisted of 1 heifer assigned to the control (CON, n = 6)) treatment and 3 heifers assigned to the rapid grain adaptation (RA) treatment (n = 18). Unequal weighting of the treatment assignment was used as more variation was expected for the RA treatment than for the CON (Bevans et al., 2005). All the RA heifers we randomly assigned to a gradual grain transition treatment (n = 18) throughout a 29-d experiment period. Briefly, heifers were initially fed a diet containing 0% grain (days 1 to 4) and transitioned to a final diet containing 89% grain using the following intermediate diets: 40% grain (days 5 to 8), 60% grain (days 9 to 12), 72% grain (days 13 to 16), 85% grain (days 17 to 20), and 89% grain (days 21 to 29). Rumen papillae biopsies were collected from the ventral sac of the rumen when the heifers were fed 0% grain on day 3 (the 3rd day of a 4-day feeding period), 72% grain on day 15 (the 3rd day of a 4-day feeding period), and 89% grain on day 26 (the 6th day of a 9-day feeding period). The biopsies were immediately washed with sterile 0.01 M phosphate-buffered saline buffer (pH 6.8) and transferred into RNA-later solution (Invitrogen, Carlsbad, CA). After being kept at 4 °C overnight, all the samples were stored at −20 °C for further analysis. To identify the molecular mechanisms responsible for the varied response to a gradual high grain diet transition, the ruminal papillae samples from RA heifers were used in this study and due to three of animals did not have biopsy samples for all three treatment periods, only fifteen animals (under three dietary conditions) were subjected to transcriptome analysis.
Ruminal pH measurement
Ruminal pH was measured using a stand-alone submersible continuous ruminal pH measurement system following the procedures as reported previously . Briefly, the pH meter was inserted into the rumen through a cannula on day 1 and was removed every four days for standardization. Ruminal pH data were collected every 30 s for each diet treatment (0%, 72% and 89% grain diet), and the daily minimum, mean, and maximum pH values were determined. In addition, the number of episodes that ruminal pH was below the threshold of 5.8, as well as the duration and area of these episodes were calculated. The area was calculated as pH unit × min/day when the ruminal pH was lower than 5.8 and the duration was calculated the time (h/day) that ruminal pH was below 5.8 in the heifers.
The acidosis index, an indicator of the severity of ruminal acidosis, was calculated by the following formula: acidosis index (pH·min/kg) = area that ruminal pH below 5.8 (min)/dry matter intake (kg) .
Ruminal papillae samples were ground while immersed in liquid nitrogen using a frozen mortar and pestle prior to RNA extraction. Total RNA was extracted from 50 mg of the resultant papillae tissue powder using Trizol reagent (Invitrogen, CA, USA) following the manufacturer’s instructions. The quantity and quality of the RNA were determined using Qubit 2.0 Fluorometer (Invitrogen, Carlsbad, CA) and Agilent 2100 Bioanalyzer (Agilent Technologies, Santa Clara, CA), respectively. The samples with RNA integrity number (RIN) greater than 8.0 were used for library construction.
RNA-seq library construction and sequencing
The RNA-seq libraries were constructed using the TruSeq mRNA Sample Preparation Kit (Illumina, San Diego, CA) according to the manufacturer’s instructions. Total RNA (1.0 μg each) from 45 samples (each with a unique index) were mixed and pooled for transcriptome sequencing (paired ends sequencing, 100 bp) at Génome Québec (Montréal, Canada) using the Illumina HiSeq 2000 system (Illumina) with all 45 samples loaded on a single chip with 9 samples per lane: each lane contains samples from 0%, 72%, and 89% grain diet. All reads were demultiplexed according to their index sequences with CASAVA version 1.8 (Illumina) and reads that did not pass the Illumina chastity filter were discarded.
RNA-seq reads mapping and annotation
RNA-seq reads were aligned to the bovine genome (UMD 3.1) using Tophat 2.0.10 with default parameters . The number of reads mapped to each gene was counted by htseq-count (http://htseq.readthedocs.io/) based on the annotation from the ENSEMBL (http://uswest.ensembl.org/) bovine gene annotation software (v75.30). The expression levels of mRNA in each library was calculated by normalizing reads to reads per million (RPM) by the following formula: RPM = (gene reads number/total mapped reads number per library) × 1,000,000.
Identification of differentially expressed genes
Among the detected genes, those that were detected in at least one sampling point with more than 1 RPM in 15 heifers were considered as expressed genes. Then, differentially expressed (DE) genes were investigated by using bioinformatics tool edgeR . The DE genes affected by diet were identified by comparing any two diets (72% vs. 0%; 89% vs. 0%) (n = 15). After that, the individual DE genes were identified by comparing the two high grain diets (89% vs. 72%) of each group (n = 10). The DE mRNA were identified by false discovery rate (FDR) < 0.05 based on Benjamini and Hochberg multiple testing correction  as well as a fold change >1.5 or < −1.5. The expression of four DE genes that related to lipid synthesis (FABP4), lipid transport (ABCA1), Na+/H+ exchange (NHE3), and ketogenesis (ACAT2) was detected using RT-qPCR.
Quantitative real-time PCR (qRT-PCR) analysis
A total of 4 genes were selected to validate the DE genes and the primers used were listed in Additional file 7: Table S1. The RT-qPCR reactions were performed with SYBR Green (Fast SYBR® Green Master Mix; Applied Biosystems) using StepOnePlus™ Real-Time PCR System (Applied Biosystems, Foster City, CA, USA) with the fast cycle and the following program: 20 s pre-denaturalization at 95°C, followed by 40 cycles of 3 s denaturation at 95°C and 30s annealing and extension at 60°C. Gene expression values were normalized to reference gene of β-actin in the same sample. The relative changes in each gene expression were calculated using the 2-ΔΔCT (cycle threshold, CT) method.
The gene list analysis tool in PANTHER classification system was used for gene ontology (GO) terms analysis of the commonly and highly expressed genes among 0%, 72%, and 89% die, and the ratio was calculated according to the number of ‘hits’ to the terms over the total number of ‘class hits’ . The functional analysis of dietary enriched/DE genes was performed by Database for Annotation, Visualization and Integrated Discovery (DAVID, http://david.abcc.ncifcrf.gov) , and Ingenuity pathway analysis (IPA, Ingenuity Systems, www.ingenuity.com). In details, we used all detected genes in rumen epithelia (11,516 genes) as background of DAVID for functional analysis (Additional file 1), and selected the database of digestive tract tissues (including forestomach, stomach, small intestine, and large intestine) and epithelial/immune cells for IPA analysis in this study. A threshold of P < 0.05 for DAVID and Benjamini and Hochberg multiple testing correction P < 0.05 for IPA, and molecules number > 2 was applied to enrich significant biological functions and pathways. The absolute value of Z-score > 1.5 was used as cutoff for activation or inhibition of biological functions and pathways, and the cutoff for networks was score > 20.
Systematic analysis of lipid transport and fatty acid metabolism related genes in DG and UG during diet transition
The expression of sodium-linked monocarboxylate transporter (SMCT1 and SMCT2), proton-linked monocarboxylate transporter (MCT1, MCT2, MCT3, and MCT4), long chain fatty acid transporter (FABP1, FABP2, FABP3, FABP4, FABP5, FABP6, FABP7, FABP9, and FABP12), cholesterol efflux protein (ABCA1, ABCA2, ABCA3, ABCA4, ABCA5, ABCA6, ABCA7, ABCA9, ABCA10, ABCA12, and ABCA13), bicarbonate transporter (PAT1), anion exchanger (AE2, AE3, AE4, DRA) were analyzed by T-test in rumen epithelia of DG and UG heifers when the diet switched from 72% to 89% grain.
The SCFA metabolism related genes including acyl-CoA synthetases (ACSS1, ACSS2 and ACSS3) that activate the volatile fatty acids, acetyl-CoA acetyl transferases (ACAT1 and ACAT2) and 3-hydroxy, 3-methylglutaryl CoA synthase (HMGCS1) that convert acetyl-CoA to 3-hydroxy, 3-methylglutaryl CoA (HMG-CoA) ,HMG-CoA lyase (HMGCL) which catalyze the synthesis of ketone bodies acetoacetate and β-hydroxybutyrate [47, 64], cholesterolgenic genes: HMG-CoA reductase (HMGCR), farnesyl-diphosphate farnesyl-transferase 1 (FDFT1), farnesyl diphosphate synthase (FDPS), and lanosterol synthase (LSS), and sterol regulatory element binding transcription factors (SREBF1 and SREBF2) were analyzed by T-test in rumen epithelia of DG and UG heifers during diet transition (from 72% to 89%).
The 5 UG and 5 DG heifers were chose for SNP analysis. RNA-Seq reads of ruminal papillae fed 0%, 72% and 89% grain diets were first combined for each heifer to increase read coverage. Then the SNP calling was performed using VarScan2 . The minimum based quality of reads was 15, minimum reads depth at a position to call a SNP was 8, and the minimum variant allele frequency threshold was 0.1. The association between alleles and varied ruminal pH response was determined by Fisher exact test, and P < 0.1 was considered as significance difference. The SNP was discarded when less than 2 reads mapped to that location in at least one heifer in this study.
- ABCA :
Cholesterol efflux protein
- ACAT :
Acetyl-CoA acetyl transferases;
- ACSS :
- ADH6 :
Alcohol dehydrogenase 6
- AE2 :
- AURKA :
Aurora kinase A
- C2/C3 :
Complement component 2/3
- CDC25B :
Cell division cycle 25 homolog B
- CDSN :
- CHAC1 :
ChaC glutathione specific gamma-glutamylcyclotransferase 1
- CKS2 :
Cyclin-dependent kinases regulatory subunit 2
- CSF2 :
Colony stimulating factor 2
C-X-C chemokine ligand 5
- CYP1A1 :
Cytochrome p450 family 1 subfamily A member 1
- CYP1B1 :
Subfamily B member 1
Database for Annotation, Visualization and Integrated Discovery
- DRA :
MHC class II DR alpha chain
- FABP :
Fatty acid binding protein
- FDFT1 :
Farnesyl-diphosphate farnesyl-transferase 1
- FDPS :
Farnesyl diphosphate synthase
False discovery rate
Genome wide association study
High grain diet
- HMGCL :
3-methylglutaryl CoA lyase
- HMGCR :
3-methylglutaryl CoA reductase
- HMGCS :
3-hydroxy, 3-methylglutaryl CoA synthase
- IFNg :
Ingenuity Pathway Analysis
Kyoto Encyclopedia of Genes and Genomes
- KLK10 :
kallikrein related peptidase 10
Low grain diet
- LSS :
- MCT :
Proton-linked monocarboxylate transporter
- NHE :
Sodium hydrogen exchanger
- NLRC5 :
NLR FAMILY card domain containing 5
- NLRC5 :
NOD-like receptor C5
- PAT1 :
Principal component analysis
Platelet-derived growth factor
- PLK1 :
Polo like kinase 1
- PTGER2 :
Prostaglandin E receptor 2
Rapid grain adaptation
RNA integrity number
Reads per million
Reverse transcription quantitative real-time PCR
Short chain fatty acid
- SMCT1 :
Sodium dependent monocarboxylate transporters
Single nucleotide polymorphisms
- SREBF2 :
Sterol regulatory element binding protein 2
- TLR3 :
Toll like receptor 3
- TMEM173 :
Tumor necrosis factor
Triggering receptor expressed on myeloid cells 1
Wierenga KT, McAllister TA, Gibb DJ, Chaves AV, Okine EK, Beauchemin KA, Oba M. Evaluation of triticale dried distillers grains with solubles as a substitute for barley grain and barley silage in feedlot finishing diets. J Anim Sci. 2010;88(9):3018–29.
Bevans DW, Beauchemin KA, Schwartzkopf-Genswein KS, McKinnon JJ, McAllister TA. Effect of rapid or gradual grain adaptation on subacute acidosis and feed intake by feedlot cattle. J Anim Sci. 2005;83(5):1116–32.
Nocek JE. Bovine acidosis: implications on laminitis. J Dairy Sci. 1997;80(5):1005–28.
Ametaj BN, Zebli Q, Iqbal S. Nutrition, microbiota, and endotoxin-related diseases in dairy cows. R Bras Zootec. 2010;39:433–44.
Gressley TF, Hall MB, Armentano LE. Ruminant nutrition symposium: productivity, digestion, and health responses to hindgut acidosis in ruminants. J Anim Sci. 2011;89(4):1120–30.
Owens FN, Secrist DS, Hill WJ, Gill DR. Acidosis in cattle: a review. J Anim Sci. 1998;76(1):275–86.
Gäbel G, Bestmann M, Martens H. Influences of diet, short-chain fatty-acids, lactate and chloride on bicarbonate movement across the reticulorumen wall of sheep. J Vet Med A. 1991;38(7):523–9.
Penner GB, Aschenbach JR, Gäbel G, Rackwitz R, Oba M. Epithelial capacity for apical uptake of short chain fatty acids is a key determinant for intraruminal pH and the susceptibility to subacute ruminal acidosis in sheep. J Nutr. 2009;139(9):1714–20.
Aschenbach JR, Penner GB, Stumpff F, Gäbel G. Ruminant nutrition symposium: role of fermentation acid absorption in the regulation of ruminal pH. J Anim Sci. 2011;89(4):1092–107.
Castillo-Lopez E, Wiese BI, Hendrick S, McKinnon JJ, McAllister TA, Beauchemin KA, Penner GB. Incidence, prevalence, severity, and risk factors for ruminal acidosis in feedlot steers during backgrounding, diet transition, and finishing. J Anim Sci. 2014;92(7):3053–63.
Schlau N, Guan LL, Oba M. The relationship between rumen acidosis resistance and expression of genes involved in regulation of intracellular pH and butyrate metabolism of ruminal epithelial cells in steers. J Dairy Sci. 2012;95(10):5866–75.
Mohammed R, Stevenson DM, Weimer PJ, Penner GB, Beauchemin KA. Individual animal variability in ruminal bacterial communities and ruminal acidosis in primiparous Holstein cows during the periparturient period. J Dairy Sci. 2012;95(11):6716–30.
Humer E, Ghareeb K, Harder H, Mickdam E, Khol-Parisini A, Zebeli Q. Peripartal changes in reticuloruminal pH and temperature in dairy cows differing in the susceptibility to subacute rumen acidosis. J Dairy Sci. 2015;98(12):8788–99.
Chen YH, Penner GB, Li MJ, Oba M, Guan LL. Changes in bacterial diversity associated with epithelial tissue in the beef cow rumen during the transition to a high-grain diet. Appl Environ Microbiol. 2011;77(16):5770–81.
Steele MA, Vandervoort G, AlZahal O, Hook SE, Matthews JC, McBride BW. Rumen epithelial adaptation to high-grain diets involves the coordinated regulation of genes involved in cholesterol homeostasis. Physiol Genomics. 2011;43(6):308–16.
Steele MA, Dionissopoulos L, AlZahal O, Doelman J, McBride BW. Rumen epithelial adaptation to ruminal acidosis in lactating cattle involves the coordinated expression of insulin-like growth factor-binding proteins and a cholesterolgenic enzyme. J Dairy Sci. 2012;95(1):318–27.
Oba M, Mewis JL, Zhining Z. Effects of ruminal doses of sucrose, lactose, and corn starch on ruminal fermentation and expression of genes in ruminal epithelial cells. J Dairy Sci. 2015;98(1):586–94.
Penner GB, Taniguchi M, Guan LL, Beauchemin KA, Oba M. Effect of dietary forage to concentrate ratio on volatile fatty acid absorption and the expression of genes related to volatile fatty acid absorption and metabolism in ruminal tissue. J Dairy Sci. 2009;92(6):2767–81.
Kleen JL, Hooijer GA, Rehage J, Noordhuizen JPTM. Subacute ruminal acidosis (SARA): a review. J Vet Med A Physiol Pathol Clin Med. 2003;50(8):406–14.
Gozho GN, Plaizier JC, Krause DO, Kennedy AD, Wittenberg KM. Subacute ruminal acidosis induces ruminal lipopolysaccharide endotoxin release and triggers an inflammatory response. J Dairy Sci. 2005;88(4):1399–403.
Penner GB, Oba M, Gäbel G, Aschenbach JR. A single mild episode of subacute ruminal acidosis does not affect ruminal barrier function in the short term. J Dairy Sci. 2010;93(10):4838–45.
Gerneke WH. Langerhans cells in epithelium of bovine forestomach - their role in primary immune response. J S Afr Vet Assoc. 1977;48(3):187–92.
Trevisi E, Amadori M, Riva F, Bertoni G, Bani P. Evaluation of innate immune responses in bovine forestomachs. Res Vet Sci. 2014;96(1):69–78.
Liu J, Bian G, Sun D, Zhu W, Mao S. Starter feeding altered ruminal epithelial bacterial communities and some key immune-related genes’ expression before weaning in lambs. J Anim Sci. 2017;95(2):910–21.
Lubbers R, van Essen MF, van Kooten C, Trouw LA. Production of complement components by cells of the immune system. Clin Exp Immunol. 2017;188(2):183–94.
Rakoff-Nahoum S, Paglino J, Eslami-Varzaneh F, Edberg S, Medzhitov R. Recognition of commensal microflora by toll-like receptors is required for intestinal homeostasis. Cell. 2004;118(2):229–41.
Uehara A, Fujimoto Y, Fukase K, Takada H. Various human epithelial cells express functional toll-like receptors, NOD1 and NOD2 to produce anti-microbial peptides, but not proinflammatory cytokines. Mol Immunol. 2007;44(12):3100–11.
Macpherson AJ, Harris NL. Interactions between commensal intestinal bacteria and the immune system. Nat Rev Immunol. 2004;4:478–85.
Chen YH, Oba M, Guan LL. Variation of bacterial communities and expression of toll-like receptor genes in the rumen of steers differing in susceptibility to subacute ruminal acidosis. Vet Microbiol. 2012;159(3–4):451–9.
Sur S, Agrawal DK. Phosphatases and kinases regulating CDC25 activity in the cell cycle: clinical implications of CDC25 overexpression and potential treatment strategies. Mol Cell Biochem. 2016;416(1):33–46.
Lobjois V, Jullien D, Bouche JP, Ducommun B. The polo-like kinase 1 regulates CDC25B-dependent mitosis entry. Biochim Biophys Acta. 2009;1793(3):462–8.
Durlacher CT, Li ZL, Chen XW, He ZX, Zhou SF. An update on the pharmacokinetics and pharmacodynamics of alisertib, a selective aurora kinase a inhibitor. Clin Exp Pharmacol Physiol. 2016;43(6):585–601.
Levine AJ. p53, the cellular gatekeeper for growth and division. Cell. 1997;88(3):323–31.
Monte M, Benetti R, Buscemi G, Sandy P, Del Sal G, Schneider C. The cell cycle-regulated protein human GTSE-1 controls DNA damage-induced apoptosis by affecting p53 function. J Biol Chem. 2003;278(32):30356–64.
Nelson DR, Koymans L, Kamataki T, Stegeman JJ, Feyereisen R, Waxman DJ, Waterman MR, Gotoh O, Coon MJ, Estabrook RW, Gunsalus IC, Nebert DW. P450 superfamily: update on new sequences, gene mapping, accession numbers and nomenclature. Pharmacogenetics. 1996;6(1):1–42.
Höög JO, Strömberg P, Hedberg JJ, Griffiths WJ. The mammalian alcohol dehydrogenases interact in several metabolic pathways. Chem Biol Interact. 2003;143:175–81.
Miyauchi S, Gopal E, Fei YJ, Ganapathy V. Functional identification of SLC5A8, a tumor suppressor down-regulated in colon cancer, as a Na+-coupled transporter for short-chain fatty acids. J Biol Chem. 2004;279(14):13293–6.
Ganapathy V, Thangaraju M, Gopal E, Martin PM, Itagaki S, Miyauchi S, Prasad PD. Sodium-coupled monocarboxylate transporters in normal tissues and in cancer. AAPS J. 2008;10(1):193–9.
Britton R, Krehbiel C. Nutrient metabolism by gut tissues. J Dairy Sci. 1993;76(7):2125–31.
Kristensen NB, Harmon DL. Splanchnic metabolism of volatile fatty acids absorbed from the washed reticulorumen of steers. J Anim Sci. 2004;82(7):2033–42.
Lane MA, Baldwin RVI, Jesse BW. Developmental changes in ketogenic enzyme gene expression during sheep rumen development. J Anim Sci. 2002;80(6):1538–44.
Hertzel AV, Bernlohr DA. The mammalian fatty acid-binding protein multigene family: molecular and genetic insights into function. Trends Endocrinol Metab. 2000;11(5):175–80.
Wang XQ, Yang K, He YS, Lu L, Shen WF. Receptor mediated elevation in FABP4 levels by advanced glycation end products induces cholesterol and triacylglycerol accumulation in THP-1 macrophages. Lipids. 2011;46(6):479–86.
Harvatine KJ, Bauman DE. SREBP1 and thyroid hormone responsive spot 14 (S14) are involved in the regulation of bovine mammary lipid synthesis during diet-induced milk fat depression and treatment with CLA. J Nutr. 2006;136(10):2468–74.
Viturro E, Koenning M, Kroemer A, Schlamberger G, Wiedemann S, Kaske M, Meyer HHD. Cholesterol synthesis in the lactating cow: induced expression of candidate genes. J Steroid Biochem Mol Biol. 2009;115(1–2):62–7.
Dempsey ME. Regulation of steroid biosynthesis. Annu Rev Biochem. 1974;43:967–90.
Baldwin RL. Use of isolated ruminal epithelial cells in the study of rumen metabolism. J Nutr. 1998;128(2):293S–6S.
Simons K, Toomre D. Lipid rafts and signal transduction. Nat Rev Mol Cell Biol. 2000;1(1):31–9.
Glavinas H, Krajcsi P, Cserepes J, Sarkadi B. The role of ABC transporters in drug resistance, metabolism and toxicity. Curr Drug Deliv. 2004;1(1):27–42.
Sahoo S, Aurich MK, Jonsson JJ, Thiele I. Membrane transporters in a human genonne-scale metabolic knowledgebase and their implications for disease. Front Physiol. 2014;5:91.
Denis M, Bissonnette R, Haidar B, Krimbou L, Bouvier M, Genest J. Expression, regulation, and activity of ABCA1 in human cell lines. Mol Genet Metab. 2003;78(4):265–74.
Jessup W, Gelissen IC, Gaus K, Kritharides L. Roles of ATP binding cassette transporters A1 and G1, scavenger receptor BI and membrane lipid domains in cholesterol export from macrophages. Curr Opin Lipidol. 2006;17(3):247–57.
Chun YS, Shin S, Kim Y, Cho H, Park MK, Kim TW, Voronov SV, Di Paolo G, Suh BC, Chung S. Cholesterol modulates ion channels via down-regulation of phosphatidylinositol 4,5-bisphosphate. J Neurochem. 2010;112(5):1286–94.
Dopico AM, Bukiya AN, Singh AK. Large conductance, calcium- and voltage-gated potassium (BK) channels: regulation by cholesterol. Pharmacol Ther. 2012;135(2):133–50.
Liao JK, Laufs U. Pleiotropic effects of statins. Annu Rev Pharmacol. 2005;45:89–118.
Alberts B, Johnson A, Lewis J, Raff M, Roberts K, Walter P. Molecular biology of the cell. 5th ed. New York: Garland Science; 2007.
Phillips T. Regulation of transcription and gene expression in eukaryotes. Nat Educ. 2008;1(1):199.
Penner GB, Beauchemin KA, Mutsvangwa T. An evaluation of the accuracy and precision of a stand-alone submersible continuous ruminal pH measurement system. J Dairy Sci. 2006;89(6):2132–40.
Kim D, Pertea G, Trapnell C, Pimentel H, Kelley R, Salzberg SL. TopHat2: accurate alignment of transcriptomes in the presence of insertions, deletions and gene fusions. Genome Biol. 2013;14(4):R36.
Robinson MD, McCarthy DJ, Smyth GK. edgeR: a bioconductor package for differential expression analysis of digital gene expression data. Bioinformatics. 2010;26(1):139–40.
Benjamini Y, Hochberg Y. Controlling the false discovery rate-a practical and powerful approach to multiple testing. J R Statist Soc B. 1995;57(1):289–300.
Mi H, Muruganujan A, Casagrande JT, Thomas PD. Large-scale gene function analysis with the PANTHER classification system. Nat Protoc. 2013;8(8):1551–66.
Huang DW, Sherman BT, Lempicki RA. Systematic and integrative analysis of large gene lists using DAVID bioinformatics resources. Nat Protoc. 2009;4(1):44–57.
Rémond D, Ortigues I, Jouany JP. Energy substrates for the rumen epithelium. Proc Nutr Soc. 1995;54(1):95–105.
Koboldt DC, Zhang QY, Larson DE, Shen D, McLellan MD, Lin L, Miller CA, Mardis ER, Ding L, Wilson RK. VarScan 2: somatic mutation and copy number alteration discovery in cancer by exome sequencing. Genome Res. 2012;22(3):568–76.
The authors acknowledge the technical assistance by X. Sun for the RNA-seq library construction and the members of Drs. Guan and Oba labs for the sampling.
This study is supported by Alberta Livestock Meat Agency (ALMA 2011F129R and ALMA2015P008R), NSERC discovery grant for L. L. Guan and Chinese International Postdoctoral Exchange Fellowship Program (20140020) for K. Zhao.
Availability of data and materials
All the sequencing data sets supporting the results in the current study have been deposited in the publicly available NCBI’s Gene Expression Omnibus Database (http://www.ncbi.nlm.nih.gov/geo/) (GSE89874).
The animal study was approved by the Livestock Animal Care Committee of the University of Alberta (protocol number OBA077) and all procedures were conducted following the guidelines of Canadian Council on Animal Care.
Consent for publication
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Additional file 1:
Differentially expressed genes with higher expression in cattle fed low grain diets (72% vs. 0% and 89% vs. 0%); Functional annotation of HGD and LGD enriched genes using DAVID (72% vs. 0% and 89% vs. 0%) (worksheet 2); Functional annotation of HGD and LGD enriched genes using IPA (72% vs. 0% and 89% vs. 0%) (worksheet 3). (XLSX 465 kb)
Additional file 2:
Differentially expressed genes in Down group (89% vs. 72% grain) (worksheet 1); Differentially expressed genes in Up group (89% vs. 72% grain) (worksheet 2); Networks enriched using DE genes in Down group (89% vs. 72% grain) (worksheet 3); Networks enriched using DE genes in Up group (89% vs. 72% grain) (worksheet 4); The detected SNPs (worksheet 5). (XLSX 67 kb)
Additional file 3: Figure S1.
The patterns of ions transportation related genes expression in DG and UP. The legend represents the RPM value scaled by rows, the blue means highly expressed in the 72% grain and yellow means highly expressed in the 89% grain. The data were analyzed by T-test, * indicated P < 0.1, ** indicated P < 0.05, and *** indicated P < 0.01. ↓ means decreased and ↑ means increased when the diet transitioned from 72% grain to 89% grain within each group. (PDF 59 kb)
Additional file 4: Figure S2.
RT-qPCR validation of selected target genes identified by RNA-seq. The gene expressions detected by RT-qPCR are shown by line graphs on the top and values are shown on the right Y-axis as relative abundance. The gene expressions detected by RNA-seq are shown by bar graphs on the bottom and values are shown on the left Y-axis as log2RPM. The data were analyzed by T-test, * indicated P < 0.1, ** indicated P < 0.05. (PDF 105 kb)
Additional file 5: Figure S3.
Single nucleotide polymorphisms (SNPs) associated with the varied ruminal pH response. (A) Fisher exact test of the association between SNPs and varied ruminal pH response. (B) The sequences of the SNP (g46,834,311 A > G). (PDF 81 kb)
Additional file 6: Figure S4.
Illustrative comparison of the gene networks related to fatty acid transport and metabolism between Down group and Up group. ↓ means decreased and ↑ means increased when the diet transitioned from 72% grain to 89% grain within. (PDF 131 kb)
Additional file 7: Table S1.
Primers sequences used for RT-qPCR. (DOCX 61 kb)
Rights and permissions
Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.
About this article
Cite this article
Zhao, K., Chen, Y.H., Penner, G.B. et al. Transcriptome analysis of ruminal epithelia revealed potential regulatory mechanisms involved in host adaptation to gradual high fermentable dietary transition in beef cattle. BMC Genomics 18, 976 (2017). https://doi.org/10.1186/s12864-017-4317-y
- Ruminal epithelia
- Dietary transition
- Beef cattle