Long-Term Effects of Subacute Ruminal Acidosis (SARA) on Milk Quality and Hepatic Gene Expression in Lactating Goats Fed a High-Concentrate Diet

Purpose The mechanism underlying the decline in milk quality during periods of feeding high-concentrate diets to dairy ruminants is not well documented. The aim of this study was to investigate the metabolic changes in the liver that contribute to the input of substrate precursors to the mammary gland after feeding a high-concentrate diet to lactating goats for a long period. Experimental Design Eight mid-lactating goats with rumen cannulas were randomly assigned to two groups. For 9 weeks, the treatment group was fed a high-concentrate diet (60% concentrate of dry matter, HC) and the control group was fed a low-concentrate diet (40% concentrate of dry matter, LC). Ruminal fluid, plasma, and liver tissues were sampled, microarray techniques and real-time polymerase chain reaction were used to evaluate metabolic parameters and gene expression in liver. Results Feeding a 60%-concentrate diet for 9 weeks resulted in a significant decrease in rumen pH. Changes in fat and protein content also occurred, which negatively affected milk quality. Plasma levels of leptin (p = 0.058), non-esterified fatty acid (p = 0.071), and glucose (p = 0.014) increased markedly in HC group. Plasma cortisol concentration was significantly elevated in the treatment group (p<0.05). Expression of the glucocorticoid receptor protein gene was significantly down-regulated (p<0.05) in the liver. The expression of genes for interleukin 1β, serum amyloid A, C-reactive protein, and haptoglobin mRNA was significantly increased (p<0.05) in the HC group. GeneRelNet analysis showed that gene expression involved in inflammatory responses and the metabolism of lipids, protein, and carbohydrate were significantly altered by feeding a high-concentrate diet for 9 weeks. Conclusions Activation of the acute phase response and the inflammatory response may contribute to nutrient partitioning and re-distribution of energy in the liver, and ultimately lead to a decline in milk quality.


Introduction
Dairy ruminants are often fed a high-concentrate diet to meet the energy requirements of high milk production. However, metabolic disorders such as subacute ruminal acidosis (SARA) can result from feeding excessive amounts of high proportions of fermentable concentrate and forage with low physically effective fiber [1]. SARA is a common digestive disorder in most dairy herds. In intensive management systems, 19% of early lactation and 26% of mid-lactation cows can be affected by SARA [2]. The consequences of SARA include depression of food intake, diarrhea, laminitis, and inflammatory response, which ultimately result in the depression of milk quality and quantity [3] [4] [5] [6].
Production costs from SARA have been estimated to be as high as $1.12 per cow per day [7]. SARA is a major concern for the global dairy industry. Especially, in the regions shortage of high quality forage and lack of reasonable feeding management practices prevents dairy cows from achieving their genetic potential for milk production.
SARA is experimentally characterized by repeated bouts of depressed ruminal pH values between 5.2 and 5.6 for at least 3 hours per day or below5.8 [5] [6]. Overfeeding fast fermentable starch or feeding excessive concentrate to forage ratios for 2 to 3 weeks can induce SARA clinical signs [2]. Most studies have focused on the short-term (2-3 weeks) effects of SARA on inflammatory response, rumen fermentation, and milk production.
Fairfield et al. [7] and Khafipoor [8] reported that experimentally induced SARA decreased milk fat percentage but increased milk protein content. However, other reports have suggested that during the short term SARA did not affect milk fat or protein content [9] [10]. However, dairy production herds affected by SARA can experience negative outcomes as long as several weeks (.3 weeks) after an episode [11]. These diseases include laminitis, weight loss, and unexplained abscesses. The long-term metabolic effects of SARA have not been well studied.
Feeding high-concentrate diets to lactating cows results in the release of bacterial endotoxins, such as lipopolysaccharide (LPS), from the rumen or hind gut [12]. Free LPS can translocate into the bloodstream from the digestive tract under conditions of high permeability and after injury to the epithelial tissue [12]. Immune responses are then activated and pro-inflammatory cytokines are released from toll-like receptors 4 (TLR-4) that recognize the circulating LPS [13]. The process of LPS recognition via TLR-4 is facilitated by the accessory molecule LPS-binding protein (LBP) and cluster of differentiation antigen 14 (CD14) [14] [15]. LBP is an acute phase protein (APP). The acute phase response to infection, tissue injuries, or other disorders participates in the restoration of physiological homeostasis [16]. The liver is the major site for the synthesis of APPs, including LBP, C-reactive protein (CRP), serum amyloid A (SAA), and haptoglobin (HP) [17] [18] [19]. Monokines such as interleukin 6 (IL-6), IL-1, and tumor necrosis factor alpha (TNF-a) induce APP synthesis via their specific hepatic receptors [20] These monokines also stimulate the release of adrenocorticotropic hormone (ACTH) from the pituitary gland. ACTH signals an increased secretion of glucocorticoid hormone from the adrenal cortex [21]. During the restoration of physiological homeostasis, these mediators act on specific receptors on different target cells, which results in changes in metabolism and systemic responses.
As the vital organ that controls metabolism [22], the liver is stimulated by endogenous and exogenous factors. It is responsible for nutrient partitioning, including homeostatic and long-term homeorhetic adaptations [23]. Physical and psychological stress elevates plasma IL-6 and APP levels in humans and in experimental animals [24] [25]. Although the mechanism of APP induction in response to stress is yet to be elucidated, activation of the hypothalamic-pituitary-adrenal (HPA) axis by stressors may trigger systemic or local cytokine production. It's reported that, during stress response physiological processes are aimed on redistribution of energy utilization in specific organs, inhibiting or stimulating energy mobilization according to priority [26]. As a fact, defense against infectious diseases and stress response are resource expensive [27]. Therefore, the objective of this study was to investigate the changes of acute phase response, stress and inflammatory response, as well as the metabolic changes in liver that may contribute to the input of substrate precursors to the mammary gland after feeding a high-concentrate diet to lactating goats for a long period.

Ruminal fluid pH
After 1 week of feeding a high-concentrate diet, the mean ruminal pH value in the HC group was pH ,5.8 from 1 to 5 hours after feeding, which was consistent with the experimental definition of SARA (figure 1A). During week 4 to week 8, enhanced ruminal buffering capacity and diet adaptation resulted in an increase in pH of HC group (mean pH .6.0), but it was still lower than that of LC group (figure 1B, C and D).

Plasma biochemical parameters and milk composition
In HC goats compared with LC counterparts, the levels of plasma non-esterified fatty acid (NEFA; p = 0.071) and leptin (p = 0.058) had a tendency to increase and blood glucose (p = 0.014) was significantly increased after feeding high concentrate diet for 9 weeks ( Table 1).
The percentage of milk fat and protein decreased significantly from 5-9 weeks in HC group compared with LC group (p,0.05). However, milk lactose concentration increased significantly in HC group from the 4th week of treatment compared to LC group (p,0.05; Table 2).

LPS content in ruminal fluid and plasma
LPS levels in ruminal fluids of HC group increased from 38,854 to 48,064 EU/mL before feeding in the morning and increased from 63,359 to 70,418 EU/mL after 8 h feeding compared to LC group. Plasma LPS concentration of HC goats was increased from 0.24 to 0.32 EU/mL before feeding in the morning and from 0.09 to 0.14 EU/mL after 4 h feeding compared with LC group. However, there was no significant difference in plasma LPS concentration and ruminal LPS content between HC and LC group (p.0.05; Table 3).

Cytokines, acute phase proteins
Plasma concentrations of IL-6 and TNF-a increased in HC group but were not significantly different compared to the LC group (p.0.05; Table 4). Compared with LC group, plasma SAA concentration in HC group was increased significantly (p,0.05), and plasma HP concentration also showed a tendency to increase (p = 0.08). Results for plasma LBP were not significantly different for the HC group compared to the LC group (p = 0.163).

Cortisol and glucocorticoid receptor
The level of plasma cortisol was significantly elevated in HC group (p,0.05). Glucocorticoid receptor (GR) protein expression in the liver was greatly down-regulated in HC group compared with LC group (p,0.05; figure 2).

Liver mRNA Expression Profiling by Microarray
Approximately 512 genes were observed to be differentially expressed in the HC group with approximately two-thirds of the genes exhibiting increased expression in liver. After removal of redundant and unannotated sequences (p,0.01), 43 genes were significantly down-regulated and 112 genes were significantly upregulated in treated goats compared to LC group (figure 3). Functional annotation of enriched Gene Ontology terms indicates that a number of biological processes within the liver are affected by consumption of high-concentrate diets, including those related to control of cell surface receptor linked signal-transduction, Gprotein coupled receptor protein signaling, regulation of transcription, immune response, and proteolysis.
The raw data and processed files for the two libraries have been deposited in NCBI's Gene Expression Omnibus database (http:// www.ncbi.nlm.nih.gov/geo/) under accession number GPL17623 and GPL17624.

The GeneRelNet (Co-expression network) analysis
Co-expression network analysis was used to find K-core regulatory factors (genes) (figure 4A-D). For each pair of genes in a cluster, we calculated the Pearson correlation coefficient and chose the significant correlation pairs to construct the network (absolute value of interaction .0.999). These genes were mainly involved in immune and inflammatory responses (n = 8), lipid metabolism (n = 5), protein metabolism (n = 26), carbohydrate metabolism (n = 7), out of a total of 155 differential expression genes (p,0.01).

Real-time polymerase chain reaction (PCR) analysis involved in inflammatory response, glucose metabolism and lipid metabolism
The mRNA expression of hepatic TLR-4, IL-1b, SAA, CRP, and HP were significantly increased in HC group compared to LC counterparts (p,0.05; figure 5A). The expressions of hepatic TLR2, CD14, MyD88, TNF-a, and LBP did not show significant difference between HC and LC group (p.0.05). Among these genes, only HP mRNA expression showed a more than 2-fold difference between the two groups, which was consistent with the data from the transcriptome microarray.
Expression of hepatic phosphoenolpyruvate carboxykinase (PC) mRNA was significantly down-regulated in HC group compared to LC group (p,0.05) (figure 5B). Hexokinase (HK) mRNA expression had a tendency to decrease in HC group (p = 0.07). However, other gluconeogenesis-and glucolysis-related genes transcription was not significantly altered in HC group compared to LC group (p.0.05).
The mRNA expression of hepatic acetyl-CoA acyltransferase 2 (ACAA2), acyl-CoA synthetase long-chain family member-1 (ACSL1), and hydroxyacyl-CoA dehydrogenase (HADH) genes were significantly increased (p,0.05) in HC group compared to LC group (figure 5C). There was a .2-fold difference in hepatic fatty acid desaturase 1 (FADS1) mRNA expression, but did not reach statistical significance due to the high variability between individual values. This result was consistent with the data from the microarray high-throughput data analysis.

Discussion
Dynamic changes in ruminal pH during SARA Dairy ruminants are often fed high grain diets to meet the energy demand for high milk production. As a result, ruminal acidosis (especially SARA) often occurs during the periods of early and mid-lactation in dairy production herds [13] [11]. Due to the rapid fermentation and the accumulation of volatile fatty acid (VFA) in the rumen, as well as the increase of lactic acid, the value of rumen pH markedly decreased and prolonged for long time [28] [29]. In experimental study, rumen pH value has been considered as the most effective indicator of SARA. The thresholds of ruminal pH for the diagnosis of SARA have been set as ruminal pH ,5.6 [10] [30] [31] and ,5.8 [28] lasted for 3 hours per day using ruminal fluids detected by pH meter system in dairy cows [10]. Consistently, ruminal pH value ,5.8 lasted more than 3 hours was considered as the critical value of SARA diagnosis [32].
In this study, goats fed a high-concentrate diet for 2 weeks exhibited a significantly lower ruminal pH, which decreased to ,5.8 and persisted for more than 3 hours per day after feeding. According to the definition of experimental SARA, HC goats were suffered SARA disease. However, although ruminal pH value of HC goats remained significantly lower than LC goats from week 4 (4W) to week 8 (8W), the ruminal pH in HC goats was markedly increased to more than 6.0. The increase in ruminal pH of HC goats indicated that the rumen was adapting to the highconcentrate diet, which was consistent with the study conducted in dairy cows by Steele et al [33].

Inflammatory response and HPA axis activation
It is widely accepted that during experimentally-induced SARA, the free LPS contents were increased in rumen fluids, yet LPS concentration in blood was vague after feeding HC diet that has been enriched with starch [15]. A decline in ruminal pH during SARA causes death and cell lysis of Gram-negative bacteria, the predominant bacterial group in the rumen, which results in an increase in free ruminal LPS concentration [12] [15,17]. Endotoxin produced in the digestive tract can be translocated into the bloodstream, which increases the concentration of blood LPS [5] [7]. It's generally accepted that LPS translocation into the bloodstream occurs as a consequence of grain-induced SARA, and that the intestines are likely the main site of LPS translocation [15] [16]. In grain-based SARA, more starch may flow into the hind gut especially into the large intestine where LPS may be produced in a significant amount as in the rumen via microorganism fermentation [9]. After LPS translocates into the peripheral circulatory system, the inflammatory and acute phase response will be activated [12]. Our results showed that, LPS concentrations were increased from 38,854 to 48,064 EU/mL before feeding in the morning and increased from 63,359 to 70,418 EU/mL after 8h feeding compared to LC group. Plasma LPS concentration of HC goats was increased from 0.24 to 0.32 EU/mL before feeding in the morning and from 0.09 to 0.14 EU/mL after 4 h feeding compared with LC group. However, due to the limited animal sizes and the individual variations the difference of LPS concentrations in ruminal fluids and plasma between HC and LC goats did not reach the statistical significance, it can't be discarded whether this moderate increase of LPS in body fluids involved in activation of inflammatory response. The activation of TLR4 and IL-1b mRNA expression in liver of HC goats also supported the above speculation. Compared to LC group, HC goats showed a significant up-regulation in mRNA expression of APPs, including hepatic SAA, CRP, and HP mRNA, and paralleled a marked increase of SAA and HP in the plasma. These results were similar to results reported in dairy cattle [7] [13] [15].
The previous studies showed that the stress axis is activated when SARA is induced [34]. Cytokines and APPs are strong inducers that signal the HPA axis to release glucocorticoids from the adrenal gland [11] [35] [36]. In the present study, we found that goats fed a high-concentrate diet for 9 weeks had a significantly higher concentration of peripheral blood cortisol, which indicated that the HPA axis was activated during the longterm SARA challenge. This is consistent with limited, preliminary evidence that high concentrate diet (70% dry matter) adjustment to one-year-old Holstein heifers for 30 days increased blood total glucocorticoids from 7.9 to 14.8 ng/mL [37]. It's suggested that the increased concentration of glucocorticoids may be reflecting a metabolic adaptation to high concentrate feed [37]. In this study, hepatic GR protein expression was significantly down-regulated in liver of HC goats indicating a negative feedback from the higher level of plasma cortisol. As a primary and vital organ that is involved in the immune response and in nutrient metabolism, the liver controls the input of precursors of milk components to the mammary gland [24] [38]. These changes of HPA axis may be aimed on redistribution of energy utilization in specific organs, inhibiting or stimulating energy mobilization according to priority [26].

The profile of mRNA expression in liver
In this study, microarray technology was used to investigate the hepatic gene expression profile and high-throughput data analysis was used to compare the differential expression of functional genes between HC and LC group. The Co-expression network of LC group shows in Figures S1, S2, S3, and S4. The activation of inflammatory and stress response in HC group was also confirmed from the GeneRelNet (Co-expression network) dataset. The results showed that the significant up-regulation (FC .2, p,0.05) of keycore gene regulated genes included CD1 (NM_001123002) and CD4 (NM_001129902) molecules, IL-2 (X55641.1) and corticotropin releasing hormone binding protein (CRHBP, NM_001009339) in liver of HC goats. These genes are important for immune response, cell adhesion, and antigen presentation. It's well known that CD4 molecule plays a key role for IL-2 produce [39] and CD1 protein is important for the innate and adaptive immune response [40]. Moreover, the significant up-regulation of 3-hydroxybutyrate dehydrogenase 1 (BDH1, NM_001034600.1) and Dipeptidase 1 (DPEP1, NM_001034472.1), and the downregulation of glutathione peroxidase 2 (GPX2, NM_001163139.1) and aldehyde dehydrogenase 18A1 (ALDH18A1, NM_001046181.1) indicated the decrease of anti-oxidative capacity in liver of goat feeding high concentrate diet [41] [42] [43]. In a good agreement, feeding high concentrate diets enrich starch to dairy cow induced SARA disease and led to systemic inflammatory response due to the translocation of LPS from gastrointestinal tract into the peripheral circulation [10] [41] [42].
In addition to being an energy source, protein plays a functional role in many organs including the liver. The changes in protein metabolism are critical to maintaining cell size and directly impacts cell survive. Gene ontology analysis also identified 26 differential expression genes related to protein metabolism. Majority of these genes are involved in the biological processes of proteolysis, collagen catabolism, and purine nucleotide metabolism. Among these genes, two ATPase (NM_001083653.1, XM_001788345.1), branched chain amino-acid transaminase 1 (BCAT1, NM_001083644.1) and fascin homolog 3 (FSCN3, NM_001075543.2) were significantly up-regulated in HC goats, while RNA polymerase III (NM_001083750.1), and three kinases (XM_866451.4, NM_001109962.1, NM_001078008.1) were markedly down-regulated. These changes likely indicated the priority of protein catabolism in liver of HC goats compared to LC. It's reported that the elevation of plasma cortisol levels inhibits protein synthesis while increasing protein breakdown [43]. In this study, the relationship between increased plasma cortisol level and protein breakdown still warrants further investigation.
The priority of catabolism was also observed in lipids metabolism. Our results showed that the transcription of genes encoding three key enzymes acting in fatty acids b-oxidation [44], ACAA2, ACSL1 and HADH were up-regulated in liver of HC goats, which was consistent with the increase of NEFA level in blood. As a metabolic hormone, leptin plays vital roles on food intake, energy expenditure, and body weight mediated by leptin receptors situated in the hypothalamus [45]. Leptin has also been shown to have effects on whole body glucose and lipid metabolism [46]. In this study, the significant increase of leptin in HC goats paralleled the increase of NEFA, which indicates the stimulatory effect of leptin on lipids catabolism. With respect to glucose metabolism, key-core gene glycogen synthase (GYS2, NM_001192905.1) mRNA transcription was down-regulated in HC goats. In addition, Q-PCR results showed that mRNA expression of two important gluconeogenic genes hexokinase (HK) and pyruvate carboxylase (PC) genes was significantly downregulated in liver of HC group. These data suggests the decrease of glucogen synthesis and gluconeogenesis in liver of HC group. However, this is contradictory to the increase of blood glucose concentrations, which was significantly elevated in HC group. This increase of blood glucose might be related to the production of lactose in these lactating animals [47]. Previous study showed that fatty acids have shown to produce an insulin-resistant state via cascade signaling pathway involving protein kinase C (PKC) and AKT phospholation [48]. Also, high levels of free fatty acids in the plasma have been shown to induce insulin resistance [49]. However, whether the increase of blood glucose in HC group is due to insulin-resistance is still unclear.   Taken as a whole, our data suggest that feeding HC diet to lactating goats for long period caused the activation of the inflammatory and stress response, and altered substrates metabolism in liver. Based on the data presented in this study, limited conclusions can be drawn from just measuring mRNA abundance because the post-transcription, protein translation, and the amount of active protein are dependent on many factors [50]. To our knowledge, it's the first time to report the changes of metabolism in liver after feeding HC diet to lactating ruminants, and the determination of wide-genome expression in liver will contribute to find the specific biochemical pathways.

Conclusion
In summary, feeding a high-concentrate diet to lactating goats for the long-term will result in the depression of milk fat and milk protein. Activation of the stress response, the acute phase response stimulation and the inflammatory response may contribute to nutrient partitioning and re-distribution of energy in the liver, and ultimately lead to a decline in milk quality.

Animals
Eight health multiparous mid-locating goats with average initial BW 4768 kg (mean 6 SD) were housed in individual stalls in a standard animal feeding house at Nanjing Agricultural University (Nanjing, China). Animals were cannulated in rumen and were randomly allocated to two groups, one receiving diets with low concentrate (40% of dry matter) as the LC group (n = 4), and another receiving high concentrate diet (60% of dry matter) as the HC group (n = 4) for 9 weeks. Ingredient and nutrient composition of experimental diets were presented in Table 5. Goats were free access to fresh water throughout the experimental time.

Milk composition analysis
Fifty milliliter of milk were sampled every two weeks into the vials with potassium dichromate at 8:00 in the morning, and stored at the 4uC until milk fat, protein and lactose analysis conducted by the commercial company of Nanjing Weigang dairy industry Co., Ltd. The data were analyzed by Univariate using the General Linear. Models of SPSS 11.0 for Windows (StatSoft Inc, Tulsa, OK, USA).

Rumen fluid collecting and measurements
Fifteen minutes prior to feed delivery and 1, 2, 3 and 5 h (or 1, 3, 5 h) after feed delivery on 2 consecutive days during 2, 4, 6 and 8 week of experimental period, 15 mL rumen fluids was collected with nylon bag and divided into 5 mL for pH measurement immediately with pH-meter. Another 10 mL were transferred into tubes and kept on ice and then centrifuged at 1,9006g at 4uC for 30 min. After that, the supernatant was collected and stored at 220uC until for LPS analysis. The concentration of LPS in rumen fluid was determined by a chromogenic endpoint assay (CE64406, Chinese Horseshoe Crab Reagent Manufactory Co., Ltd., Xiamen, China) with a minimum detection limit of 0.1 EU/mL. The procedures were performed according to the manufacture's instruction. The data of pH and LPS concentrations were analyzed by Univariate using the General Linear Models of SPASS 11.0 for windows (StatSoft Inc, Tulsa, OK, USA).

Measurement of plasma biochemical parameters
At the end of experiment, plasma was sampled fifteen minutes prior to feed delivery using EDTA-containing vacuum tubes from jugular vein. Plasma leptin concentration was measured with radioimmunoassay (RIA) using commercial kits purchased from Shanghai Institute of Biological Products. The ranges of assay sensitivity were between 0.5 and 24 ng/mL. The inter-and intraassay coefficients of variation were 5% and 10%, respectively. The concentration of NEFA (A042) and glucose was detected by commercial kits purchased from Nanjing Jiancheng Bioengineering institution, and the procedures were performed according to the manufacture's instruction. The concentrations of leptin and glucose were carried out in Nanjing General Hospital of Nanjing Military. Inflammatory cytokines such as IL-6 (RGB-60023G, Beijing Rigorbio Science Development Co., Ltd.), TNFa (RGB-60080G) and acute phase protein such as LBP (RGB-60266G), HP (H086-09, abcam), SAA (ab100635, abcam) were detected by ELISA kits. The concentration of LPS in plasma was determined by a chromogenic endpoint assay with diazo coupling reagent (CE80545, Chinese Horseshoe Crab Reagent Manufactory Co., Ltd., Xiamen, China) with a minimum detection limit of 0.01 EU/mL. The procedures were performed according to the manufacture's instruction. The data were analyzed by Independent-Samples T test using the Compare Means of SPASS 11.0 for Windows (StaSoft Inc, Tulsa, OK, USA).

RNA extraction and real-time quantitative PCR
The liver tissues were sampled by a punch biopsy with a local anesthesia. Total RNA was extracted from each dairy goat using TRIZOL reagent (Takara, Dalian, China) according to the manufacturer's protocols. Total RNA concentration was then quantified by measuring the absorbance at 260 nm in a spectrophotometer (Eppendorf Biotechnology, Hamburg, Germany). Aliquots of RNA samples were subjected to electrophoresis with 1.4% agarose-formaldehyde gels stained with ethidium bromide to verify their integrity. Reverse transcription (RT) was performed using the total RNA (2 mg) in a final volume of 25 mL containing 16RT-buffer, 100 U reverse transcriptase, 8 U RNase inhibitor (Promega, USA), 5.3 mmol/L random hexamer primers and 0.8 mmol/L dNTP (TaKaRa, Dalian, China). After incubation at 37uC for 1 h, the reaction was terminated by heating at 95uC for 5 min and quickly cooling on ice. Real-time PCR was performed in Mx3000P (Stratagene, La Jolla, CA, USA). Mock RT and No Template Controls (NTC) were set to monitor the possible contamination of genomic DNA both at RT and PCR. Melting curves were performed to insure a single specific PCR product for each gene. Two microliter of 16-fold dilution of RT product was used for PCR in a final volume of 25 mL containing 12.5 mL SYBR Green Real-time PCR Master Mix (TOYOBO Ltd., Shanghai, China) and 0.6-0.8 mM of each forward and reverse primers for target genes (as shown in Table 6). Goat bactin mRNA was used as a reference gene for normalization purposes. The following PCR protocols were initial denaturation (1 min at 95uC), then a three-step amplification program (20 s at 95uC, 20-30 s at 60-64uC, 30 s at 72uC) was repeated 45 times. The method of 2 2DDCt was used to analyze qPCR data with Independent-Samples T test using Compare means of SPSS 11.0 for Windows to statistical analysis. All samples were included in the same run of RT-PCR and repeated for at least 3 times. The primers were shown in Table 6. The data was analyzed by Independent-Samples T test using the Compare Means of SPASS 11.0 for Windows (StaSoft Inc, Tulsa, OK, USA).

Western blotting
Total protein was extracted from frozen liver tissue, and the protein concentration was determined by the BCA assay (Pierce, Rockford, IL, USA). Fifty micrograms of protein extract from each sample were subjected to electrophoresis on a 10% SDS-PAGE gel, and the separated proteins were transferred onto the nitrocellulose membranes (Bio Trace, Pall Co., USA). Western blot analysis for GR (sc-1004, Santa Cruz Biotechnology, 1:500) is performed with primary antibodies and corresponding HRPconjugated secondary antibodies. b-actin (KC-5A08, KangChen Bio-tech, China, 1:5000) was used as a reference protein for normalization purposes in the Western blot analysis. Finally, the blot was washed and detected by enhanced chemiluminescence (ECL) using the LumiGlo substrate (Super Signal West Pico Trial Kit, Pierce, USA). ECL signals were recorded by an imaging system (Bio-Rad, USA) and analyzed with Quantity One software (Bio-Rad, USA). Values of GR protein was presented as fold change relative to the average value of LC group. The data was analyzed by Independent-Samples T test using the Compare Means of SPASS 11.0 for Windows (StaSoft Inc, Tulsa, OK, USA).

mRNA Microarray Experiment
Transcriptome microarrays were performed for seven dairy goats (3 goats from the LC group and 4 goats from the HC group) for hepatic gene expression profiling. Because complete genome sequence data was not available, gene sequences from three animal species were employed to design the gene chip (goat, sheep, and Bos Taurus). Goat expressed sequence tag (EST) sequence (n = 13898), goat gene sequence (n = 720), sheep gene sequence (n = 2718), and Bos Taurus gene sequence (n.35000) transcripts were investigated using the microarray. There were 61012 probe IDs and 2-fold difference probe IDs (n = 2140). Reference sequences were cited from http://www.ncbi.nlm.nih.gov/gene/ and EST database was from http://www.ncbi.nlm.nih.gov/ nucest/. The FDR (False Discovery Rate) were calculated. Differentially expressed genes (DEGs) were selected with FDR ,5% and FDR ,10%. All data were MIAME compliant and have been deposited in GEO (accession number GPL17623 and GPL17624).

GO analysis
GO analysis was applied to analyze the main function of the differential expression genes according to the Gene Ontology which is the key functional classification of NCBI, which can organize genes into hierarchical categories and uncover the gene regulatory network on the basis of biological process and molecular function [51]. Specifically, two-side Fisher's exact test and x 2 test were used to classify the GO category, and the false discovery rate (FDR) [52] was calculated to correct the P-value,the smaller the FDR, the small the error in judging the p-value. The FDR was defined as FDR~1{ N k T , where N k refers to the number of Fisher's test P-values less than x 2 test P-values. We computed P-values for the GOs of all the differential genes. Enrichment provides a measure of the significance of the function: as the enrichment increases, the corresponding function is more specific, which helps us to find those GOs with more concrete function description in the experiment. Within the significant category, the enrichment Re was given by: Re~(n f =n)=(N f =N) where ''n f '' is the number of flagged genes within the particular category, ''n'' is the total number of genes within the same category, ''N f '' is the number of flagged genes in the entire microarray, and ''N'' is the total number of genes in the microarray [53].

GeneRelNet (Co-expression network)
Gene co-expression Networks was established to identify the interactions among genes. Gene co-expression Networks were built according to the normalized signal intensity of specific expression genes [54]. For each pair of genes, we calculate the Pearson correlation and choose the significant correlation pairs with which to construct the network [55]. In a network analysis, degree centrality is the simplest and most important measures of a gene centrality within a network that determining the relative importance. Degree centrality is defined as the link numbers one node has to the other [56]. Moreover, to study a variety of properties of networks, k-cores were introduced in graph theory as a method of simplifying graph topology analysis. A k-core of a network is a subnetwork in which all nodes are connected to at least k other genes in the subnetwork. The Purpose of Network Structure Analysis is to locate core regulatory factors (genes). In one network, core regulatory factors connect most adjacent genes and have the biggest degrees. While considering different networks, Core regulatory factors were determined by the degree differences between two class samples [57].