Skip to main content
Advertisement
Browse Subject Areas
?

Click through the PLOS taxonomy to find articles in your field.

For more information about PLOS Subject Areas, click here.

  • Loading metrics

Epigenetic Mechanisms Contribute to the Expression of Immune Related Genes in the Livers of Dairy Cows Fed a High Concentrate Diet

  • Guangjun Chang,

    Affiliation College of Veterinary Medicine, Nanjing Agricultural University, Nanjing, P. R. China

  • Kai Zhang,

    Affiliation College of Veterinary Medicine, Nanjing Agricultural University, Nanjing, P. R. China

  • Tianle Xu,

    Affiliation College of Veterinary Medicine, Nanjing Agricultural University, Nanjing, P. R. China

  • Di Jin,

    Affiliation College of Veterinary Medicine, Nanjing Agricultural University, Nanjing, P. R. China

  • Junfei Guo,

    Affiliation College of Veterinary Medicine, Nanjing Agricultural University, Nanjing, P. R. China

  • Su Zhuang,

    Affiliation College of Animal Science and Technology, Nanjing Agricultural University, Nanjing, P. R. China

  • Xiangzhen Shen

    xzshen@njau.edu.cn

    Affiliation College of Veterinary Medicine, Nanjing Agricultural University, Nanjing, P. R. China

Abstract

Purpose

Epigenetic modifications critically regulate the expression of immune-related genes in response to inflammatory stimuli. It has been extensively reported that a high concentrate (HC) diet can trigger systemic inflammation in dairy cows, yet it is unclear whether epigenetic regulation is involved in the expression of immune genes in the livers of dairy cows. This study aimed to investigate the impact of epigenetic modifications on the expression of immune-related genes.

Experimental Design

In eight mid-lactating cows, we installed a rumen cannula and catheters of the portal and hepatic veins. Cows were randomly assigned to either the treatment group fed a high concentrate (HC) diet (60% concentrate + 40% forage, n = 4) or a control group fed a low concentrate (LC) diet (40% concentrate + 60% forage, n = 4).

Results

After 10 weeks of feeding, the rumen pH was reduced, and levels of lipopolysaccharide (LPS) in the rumen, and portal and hepatic veins were notably increased in the HC group compared with the LC group. The expression levels of detected immune response-related genes, including Toll-like receptor 4 (TLR4), cytokines, chemokines, and acute phase proteins, were significantly up-regulated in the livers of cows fed a HC diet. Chromatin loosening at the promoter region of four candidate immune-related genes (TLR4, LPS-binding protein, haptoglobin, and serum amyloid A3) was elicited, and was strongly correlated with enhanced expression of these genes in the HC group. Demethylation at the promoter region of all four candidate immune-related genes was accompanied by chromatin decompaction.

Conclusion

After HC diet feeding, LPS derived from the digestive tract translocated to the liver via the portal vein, enhancing hepatic immune gene expression. The up-regulation of these immune genes was mediated by epigenetic mechanisms, which involve chromatin remodeling and DNA methylation. Our findings suggest that modulating epigenetic mechanisms could provide novel ways to treat systemic inflammatory responses elicited by the feeding of a HC diet.

Introduction

In the dairy industry, it is currently a common practice to feed a high concentrate (HC) diet to lactating cows to meet their energy requirements to support high milk production. Although this practice might increase short-term economic efficiency, it is not a beneficial practice for the health of cows. The effects of a HC diet on gastrointestinal tract or mammary gland health have received increasing attention [13]. However, to date, the effects of a HC diet on immune responses and liver metabolism have remained unclear, despite the extreme importance of these pathways in cow health and performance [4, 5]. The consumption of a HC diet causes a series of changes that involve a significant depression of pH and increased lipopolysacchrides (LPS) in the rumen and intestine, as well as the elevation of acute phase proteins (APPs), including LPS-binding protein (LBP), serum anyloid A (SAA), and haptoglobin (Hp), in the peripheral blood [6]. Changes in the blood concentrations of SAA and Hp, which have been reported in many studies, suggest that feeding a HC diet triggers an inflammatory response [7, 8]. Moreover, serum levels LBP are an indicator of systemic inflammation caused by circulating LPS that has translocated from the digestive tract [9]. Although many hazardous substances are released in the digestive tract and translocate into the bloodstream during feeding with a HC diet, LPS is thought to be the major contributor to systemic inflammation because most LPS is released by Gram-negative bacteria, which are the main bacterial community in the digestive tract of cows [10].

The liver is an important site where harmful substances (e.g., microbial products, LPS, and histamine) from the digestive tract encounter immune cells and stimulate the innate immune system, which is the first line of defense against invading pathogens. Innate immune cells, such as Kupffer cells, dendritic cells, and neutrophils, initiate and maintain hepatic inflammation by secreting cytokines [11]. A global gene expression analysis of liver tissues showed that the expression of cytokines and APPs was enhanced in experimental LPS intra-mammary gland infusion [12]. Moreover, consuming a HC corn straw diet caused a higher pro-inflammatory response in the mammary gland [3]. However, the effects of a HC diet on pro-inflammatory responses in the liver are largely unknown. Furthermore, consuming a HC diet leads to marked increases in LPS produced in the digestive tract and delivered directly to the liver via the portal vein [13]. These events trigger a liver inflammatory response via the TLR4-signaling pathway, resulting in the release of pro- and anti-inflammatory cytokines and APPs [14].

Excessive activation of the TLR4 signaling pathway by LPS induces endotoxin tolerance and cross-tolerance towards other pathogen-associated molecular patterns (PAMPs) [15]. Epigenetic mechanisms involved in chromatin remodeling are thought to mediate these phenomena [1618]. A recent study conducted in dairy cows showed that feeding a HC corn straw diet elicited epigenetic alterations in mammary glands [19]. Additionally, E. coli mastitis has been reported to induce chromatin compaction at the αS1-casein promoter to shut down casein synthesis, and to trigger chromatin decompaction at the β-defensin LAP promoter to facilitate NF-κB recruitment, which drives enhanced gene expression [20, 21]. Although many advances have been made in recent years, the role of epigenetic mechanisms in the regulation of immune-related genes in the liver in response to gastrointestinal tract-derived LPS is less well-known. Therefore, we hypothesized that feeding a high concentrated diet to lactating cows results in the entry of excessive LPS into the liver via the portal vein, thereby enhancing the expression of hepatic immune related genes. We predicted that epigenetic mechanisms contribute to the up-regulated expression of innate immune genes. We tested this hypothesis by analyzing chromatin remodeling and DNA methylation at the promoter regions of immune-related genes. TLR4 and LBP were selected as parameters to evaluate activation of the TLR4 signaling cascade. Hp and SAA3 were chosen as model biomarkers to assess the acute phase reaction (APR). Finally, αS1-casein was used as an unmodulated control gene in the liver based on published data [20].

Materials and Methods

Ethics Statement

The experimental protocol was approved by the Animal Care Committee of Nanjing Agricultural University in accordance with the Guidelines for Experimental Animals of the Ministry of Science and Technology (2006, Beijing, China).

Animal, Diets and Experimental Design

A total of eight multiparous mid-lactating Holstein dairy cows (average body weight, 512 ± 36 kg; average milk yield, 29.18 ± 3.23 kg, 1 to 2 months post partum) were used in this experiment and were randomly assigned to one of two groups. Cows were fed a HC (Forage: Concentrate = 4:6, n = 4) diet in the treatment group, and were fed low concentrated (LC, Forage: Concentrate = 6:4, n = 4) diet in the control group. Quarter milk samples of these cows were measured weekly two times during entire experimental period to be free from mastitis (somatic cells <100000 /mL). Before the beginning of formal experiment, all of cows were received LC diet for 4 weeks (as an adaption period) to ensure the similarity of metabolic status. Dietary ingredients and nutrient composition are shown in Table 1. The formal experiment lasted for 10 weeks. All experimental cows were housed in individual tie-stalls at the Dairy Farm of Nanjing Agricultural University (Nanjing, China). Animals were fed and milked three times daily at 04:00, 12:00, 20:00 h, and were allowed free access to fresh water. In the first week of the adaption period, cows were installed rumen cannula, and indwelt with hepatic and portal vein catheters. After surgery, cows were observed for 2 weeks during recovery from the surgery. Sterilized heparin saline (500 IU/mL) was used to prevent catheter blockage at daily 8-hour intervals until the end of the experiment.

thumbnail
Table 1. Ingredients and nutritional composition of the diets.

https://doi.org/10.1371/journal.pone.0123942.t001

Sample Collection

Milk yields were recorded daily. Rumen pH values were measured on the last 3 days at 2th, 4th, 6th, 8th, 10th week. Rumen fluid and blood samples of the portal and hepatic veins were collected on the last 3 consecutive days of the 10th week at 0 (15 min before feeding), 4, and 8 h after feeding. Rumen fluid samples were centrifuged at 10,000 ×g for 45 min, and then the supernatant was aspirated and passaged through a disposable 0.22-μm filter. The filtrate was collected in a sterile, depyrogenated glass tube (previously heated at 250°C for 2 h) and heated at 100°C for 30 min. Samples were cooled at room temperature for 15 min and stored at—20°C for subsequent LPS measurements. Plasma was isolated from blood samples by centrifugation at 1900 ×g at 4°C for 15 min, then was injected into pyrogen-free glass tubes and stored at—20°C for LPS measurements. The day after the collection of rumen fluid and blood samples was completed, liver tissue samples were obtained by punch biopsy under local anesthesia, immediately frozen in liquid nitrogen, and subsequently stored at—70°C, as previously described [12].

LPS Measurements

LPS is the component of cell wall of Gram-negative bacteria that are the predominant bacterial group in the rumen. After the feeding of HC diet, rapid growth of Gram-negative bacteria can result in the shedding of LPS in the rumen in the early stage. Meanwhile, the rapid growth of those bacteria leads to the production of massive organic acid (for example VFA, lactic acid) repressing the rumen pH, which cause the death and cell lysis to release a large amount of free LPS [10]. In addition, when starch enters into colon or caecum, the same response will happen, and these changes (the reduction of pH and increase of LPS concentration in digestive tract) can cause the disruption of gastrointestinal barrier facilitating the translocation of LPS from digestive tract into circulation [1, 2].

The concentration of LPS was determined using a chromogenic endpoint assay (cat. CE64406, Chinese Horseshoe Crab Reagent Manufactory Co., Xiamen, China) with a minimum detection limit of 0.05 EU/mL (rumen liquid) or 0.01 EU/mL (plasma). Procedures were performed according to the manufacturer’s instructions, as described by Dong et al. [22].

RNA Extraction and Real-Time Quantitative PCR (RT-qPCR)

RT-qPCR was performed using an ABI 7300 system (Applied Biosystems, Foster City, CA, USA) to determine the relative copy numbers of different mRNA transcripts. Liver samples were ground into a powder with a mortar under liquid nitrogen, and then total RNA was extracted in Trizol (Takara, Dalian, China) in accordance with the manufacturer’s protocol. For cDNA synthesis, 1.5 μg total RNA was prepared in a reverse transcription reaction (cat. RR036A, Takara) with one specific primer for each of the different mRNAs and oligo (dT) for all mRNAs. Then, cDNAs were purified using a purification kit (cat. D0033, Beyotime, Shanghai). RT-qPCR was run using gene-specific primer pairs to amplify cDNA target segments using SYBR Premix EX Taq kit (cat. DRR420A, Takara). The relative copy numbers of individual mRNA were calculated from a dilution series of 106 to 102 copies of the respective sequenced plasmids. All samples were assayed twice from two independent cDNA preparations. Reverse transcription and amplification primers are listed in S1 Table.

Chromatin Precipitation and CHART-PCR (Restriction Enzyme Digestion, Purification of DNA, and q-PCR)

A 100-mg tissue sample was ground into a powder in liquid nitrogen and suspended in 3-mL pre-chilled resuspension buffer (RSB; 10 mM Tris [pH 8.0], 3-mM MgCl2, and 10-mM NaCl) that contained 0.5% Nonidet NP 40. A freshly prepared dilution (1/200) of a proteinase inhibitor cocktail (Roche) and phenyl-methyl-sulfonyl-fluoride (PMSF; 1 mM) were added immediately prior to use. After incubation for 5 min on ice, the tissue powder was homogenized using a dounce homogenizer (Sigma). The liquid was filtered through sterilized glass wool into a pre-cooled 15-mL Eppendorf tube, then the filtrate was centrifuged at 1000 ×g for 10 min at 4°C to pellet the nuclei. Pellets were washed once in RSB buffer containing 1 mM β-mercaptoethanol and then were transferred into a pre-cooled 1.5-mL Eppendorf vial. Nuclei were pelleted again (1000 ×g, 5 min, 4°C) and re-suspended in 100-μL RSB buffer containing 50% glycerol and stored at—20°C.

Chromatin compaction was measured using the Chromatin Accessibility by Real-time PCR assay (CHART) [23], as described Vanselow et al [20] with minor modifications. Briefly, 4-μL nuclear suspension was added to 46-μL restriction digest buffer supplemented with 20-U restriction enzyme and proteinase inhibitors (see above), and mixtures were incubated at 37°C for 30 min. A control sample was prepared by omitting the enzyme and treating it similarly. Subsequently, proteins were digested (56°C, 2 h) with proteinase K (0.5 mg/mL in 50 mM Tris [pH 8.0], 0.15 M EDTA, 0.1 M NaCl, 6 mM DTT, and 1% SDS). DNA was purified using a purification kit (cat. D0033, Beyotime) and the DNA concentration was measured using a Nanodrop spectrophotometer. The restriction enzymes used to determine the amount of chromatin compaction at the different promoters were MaeI (for TLR4, LBP, and SAA3) and DdeI (for αS1-casein and Hp). It has been reported that αS1-casein encoding gene would be expressed to a considerable amount in the liver of cows suffering E.coli mastitis, and the degree of chromatin compaction at its promoter region had no difference in the liver between infectied cows and control cows [20, 24]. Hence, this gene was treated as unregulated control gene to highlight the sensitivity of immune relevant genes (TLR4, LBP, Hp and SAA3) to chromatin remodeling in the present study.

The quantity of undigested target DNA was measured by real-time PCR (ABI 7300), similarly to the protocol described above for RT-qPCR. Primer sequences are listed in S2 Table. The degree of compaction was represented as the fraction of the copy numbers determined from the digested vs undigested control samples.

Methylation Assay of Target Gene Promoters

Genomic DNA was isolated from liver tissues and analyzed by Methyl-Profiler DNA Methylation qPCR Assays with some modifications [25, 26]. This method exploits the theory that HpaII digestion can be blocked by the methylation of C nucleotides in CpG dinucleotides, whereas its isoschizomer MspI is not blocked. Purified genomic DNA was predigested with EcoRI to facilitate accurate aliquoting. Subsequently, 3-μg purified DNA was distributed and digested with HpaII or MspI (10 U, 37°C, 2 h), or was similarly treated without enzyme addition as a control. Subsequently, DNA was re-purified and re-quantified again and then 200 ng (equivalent of ~6000 gene copies) was used to determine the remaining DNA copies for the qPCR assays. Primer sequences are listed in S3 Table. The degree of methylation was calculated based on the difference between the ratio of copy numbers obtained after HpaII or MspI digestion, and each relative value is presented compared with the value of the undigested control (set as 100% of input DNA).

Western Blotting

Total proteins were extracted from smashed liver tissue with RIPA Lysis Buffer (cat. SN338, Sunshine Biotechnology Co., Nanjing, China). Protein concentrations were determined using the BCA assay (Pierce, Rockford, IL, USA). A total of 50-μg protein extract from each sample was subjected to electrophoresis on a 7.5% SDS–PAGE gel, and separated proteins were transferred onto nitrocellulose membranes (Bio Trace, Pall Co., Port Washington, NY, USA). Western blotting analysis of TLR4 (cat. sc-293072, Santa Cruz Biotechnology, Santa Cruz, CA, USA, 1:200) was carried out with the primary antibody and a corresponding HRP-conjugated secondary antibody. Additionally, β-actin (KC-5A08, Kang Chen Bio-tech, Shanghai, China, 1:5000) was used as a reference protein for normalization purposes in western blotting analysis. Then, blots were washed and detected by enhanced chemiluminescence (ECL) using the LumiGlo substrate (Super Signal West Pico Trial Kit, Pierce). ECL signals were recorded using an imaging system (Bio-Rad, Hercules, CA, USA) and analyzed with Quantity One software (Bio-Rad). The intensity of TLR4 protein bands are presented as the relative fold-change compared to the average value of the control group.

Statistical Analysis

The mixed procedure of SAS (SAS version 9.2, SAS Institute Inc.) was used to analyze LPS concentrations with a repeated measures design. The effects of treatment and time were classified as fixed factors. The effects of cows were considered randomly. The amount of time within treatments and cows were considered as repeated measurements, and compound symmetry (CS) was used as the type of covariance. Analysis of rumen pH values, milk yields, selected mRNA expression levels, the degree of chromatin compaction, the percentage of promoter methylation, and TLR4 protein expression was performed using the ANOVA package of SAS. The coefficient correlations between the mRNA expression levels and the degree of chromatin compaction, as well as between the degree of chromatin compaction and the percentages of promoter methylation were analyzed using Pearson’s model in SAS. Significant differences that were detected were defined at two levels as follows: p<0.05 (significant) or p<0.01 (highly significant).

Results

Rumen pH Values and Milk Yields

Feeding a HC diet to cows linearly reduced the rumen pH in the treatment group compared with that in control group throughout the experimental period (p<0.01), and pH values were recorded from the 4th week to the end of this experiment (Table 2). Milk yields was increased in first 4 weeks, and then gradually decreased in subsequent 6 weeks in the treatment group compared with that in the control group and obtained significantly difference at 10th week (p<0.01, Table 2) between treatment and control groups.

thumbnail
Table 2. Rumen pH values and milk yields of cows from treatment and control groups.

https://doi.org/10.1371/journal.pone.0123942.t002

LPS Concentrations in the Rumen, Portal, and Hepatic Veins

LPS concentrations in rumen fluid were higher in the treatment group than that in the control group (P = 0.04), and the highest concentration of rumen LPS was observed at 4 h after feeding in the treatment group. The kinetic effects and temporal interactions with treatment on rumen LPS were not significant (Table 3). The concentrations of LPS in the portal (p<0.01) and hepatic (p<0.01) veins were significantly increased in the treatment group compared with those in the control group. The LPS concentrations in the portal and hepatic veins were both highest 4 h after feeding in the treatment group. The effects of treatment × time on LPS concentrations in the portal vein (p = 0.02) were significant, but were not in the hepatic vein (p = 0.15). The effects of time on LPS concentrations in the portal (p = 0.07) and hepatic (p = 0.24) veins were not significantly different between the treatment and control groups (Table 3).

thumbnail
Table 3. LPS concentrations in the rumen, and portal and hepatic veins of cows from treatment and control groups.

https://doi.org/10.1371/journal.pone.0123942.t003

The Expression of Immune-Related Genes in the Liver

The expression levels of immune genes, encoding cytokines, chemokines, APPs, and TLR4, were up-regulated in the liver of the treated cows compared with control cows (Table 4). We found that the expression of the pro-inflammatory cytokines IL-1α, IL-1β, and TNF-α were similarly increased by three-fold, whereas the expression levels of IL-6 were elevated by 47-fold in the treatment group compared with those in the control group. The expression levels of acute phase proteins (LBP, SAA3, and Hp) showed highly significant differences (p<0.01) between the treatment and control groups, but the fold changes in expression varied. The expression levels of LBP, SAA3, and Hp were increased by 10.32-, 8.64-, and 3.97-fold, respectively, in the treatment group compared with those in the control group. TLR4, the specific receptor for LPS, was significantly elevated by 4.48-fold (p<0.05). Additionally, the expression of αS1-casein was not significantly different between the treatment and control groups, which we used as a control gene that was unmodulated for subsequent epigenetic analyses.

thumbnail
Table 4. The levels of mRNA expression (relative copy number) in the livers of cows from treatment and control groups.

https://doi.org/10.1371/journal.pone.0123942.t004

A Map of the Binding Sites for Relevant Transcription Factors in Promoters of Candidate Genes

The target promoter sequences of all candidate genes, including TLR4, LBP, Hp, SAA3, and αS1-casein, were analyzed using the MatInspector program (Genomatix, online tools) to search for potential binding sites for relevant transcription factors. Subsequently, a map of all target gene promoters was prepared to show the positions of relevant binding sites of transcription factors and restriction enzymes that are involved in CHART and methylation assays, as well as the positions of primer pairs associated with those assays (Fig 1).

thumbnail
Fig 1. A map of the target promoter region of candidate genes and the distribution of relevant transcription factor binding sites in the corresponding promoter regions of those genes.

Numbers refer to the transcriptional start position (TSP, +1), indicated by bold black arrows, and the 5’-position of exon 1 (black box). The gray boxes in the TLR4, Hp, and αS1-casein promoter regions indicate repeat elements. The position of the various transcription factors and restriction enzyme recognition sites are indicated by the respective symbols. The positions of primers used for CHART-PCR and the methylation assays are denoted by light black arrows. The promoter structure of αS1-casein has been described by Vanselow et al. The TSP of the SAA3 promoter was reported by Marilynn et al. The target promoter regions of candidate genes were identified by BLAST analyses as DNA-sequences that are 5’-upstream of the mRNA sequence deposited in the following NCBI files: NM_174198 (TLR4), NM_001038674 (LBP), NM_001040470 (Hp), and NW_003104637 (SAA3).

https://doi.org/10.1371/journal.pone.0123942.g001

Chromatin Remodeling of Immune Gene Promoters

Loosening of the chromatin at the promoter region of the four immune candidate genes was observed in the treatment group (Table 5). The average degree of chromatin compaction in the liver samples from the treated cows (0.26 ± 0.04) was much lower than that in samples from the control cows (0.91 ± 0.02), and this difference was highly significant (p<0.01). By contrast, the average degree of chromatin compaction at the aS1-casein promoter was very similar in liver samples from both groups (0.40 ± 0.06 vs 0.39 ± 0.08; treatment vs control, respectively). Moreover, plotting the individual degree of chromatin compaction against the relative mRNA copy numbers showed a distinct interdependence of both parameters for the four immune genes (Fig 2). The coefficients of correlation (R2) for TLR4 (R2 = -0.82, p = 0.013), LBP (R2 = -0.89, p<0.01), and Hp (R2 = -0.88, p<0.01) were all robust and significant. While the R2 of SAA3 was lower than the others, it still showed a significant correlation (R2 = -0.67, p = 0.05). Additionally, there was no significant correlation (R2 = 0.07, p = 0.88) between the degree of chromatin compaction at the aS1-casein promoter and the relative mRNA copy numbers of this gene.

thumbnail
Fig 2. The correlation between chromatin compaction and mRNA expression.

The graph shows the individual degrees of chromatin compaction (axis of abscissa) compared with the relevant mRNA copy numbers (axis of ordinate) for each respective gene, which are presented in panels A–E. Black rhombs, treated cows fed a HC diet; white circles, control cows fed a LC diet; R2, coefficient of correlation; p, significance of correlation.

https://doi.org/10.1371/journal.pone.0123942.g002

thumbnail
Table 5. A comparison of chromatin compaction in the liver of cows between treatment and control groups.

https://doi.org/10.1371/journal.pone.0123942.t005

Promoter Demethylation for Candidate Immune Genes

A single HpaII/MspI restriction site in the promoter region of these candidate genes (Fig 1) was used to probe for an eventual association of CpG demethylation with chromatin loosening. This site in all immune candidate gene promoters was hypomethylated in liver samples obtained from cows in the treatment group (Fig 3). The average degree of promoter methylation in liver samples from treatment group cows (0.14 ± 0.02) was significantly lower than that in liver samplers from control group cows (0.52 ± 0.01; p<0.01; Table 6). Compared with these immune genes, there was no significant difference in αS1-casein promoter methylation between the treatment and control groups, and the average degree of αS1-casein promoter methylation in the samples from both groups was much more similar (0.27 ± 0.06 vs 0.22 ± 0.04; treatment vs controls, respectively).

thumbnail
Fig 3. The correlation between the degree of chromatin compaction and the percentage of promoter methylation.

The individual degrees of chromatin compaction (axis of abscissa) are plotted against the percentages of promoter methylation for the respective genes, as shown in panels A–E. Black squares, treated cows fed a HC diet; white circles, control cows fed a LC diet; R2, coefficient of correlation; p, significance of correlation.

https://doi.org/10.1371/journal.pone.0123942.g003

thumbnail
Table 6. A comparison of promoter methylation in the liver of cows from treatment and control groups.

https://doi.org/10.1371/journal.pone.0123942.t006

Plotting the individual degree of chromatin compaction against the degree of promoter methylation showed a strong relationship between these two parameters for the immune candidate genes. Our data analysis of promoter methylation showed a significant correlation for the four candidate immune genes as follows: TLR4 (R2 = 0.82, p = 0.01), LBP (R2 = 0.89, p<0.01), Hp (R2 = 0.92, p<0.01), and SAA3 (R2 = 0.85, p<0.01; Fig 3). By contrast, no significant correlation (R2 = -0.49, p = 0.22) was detected between the degree of chromatin compaction and the degree of promoter methylation for the αS1-casein encoding gene (Fig 3).

TLR4 Protein Expression

Western blotting analysis showed that the expression of TLR4 protein, which initiates the TLR4-dependent innate immune signaling pathway, was significantly elevated (p<0.01) in liver samples from treatment group cows compared with control group cows (Fig 4).

thumbnail
Fig 4. The expression of TLR4 protein in the livers of cows fed two different diets.

The average relative TLR4 protein expression levels compared with reference β-actin protein levels is presented. Black filled bar, treatment group fed a HC diet (n = 4, mean ± SEM); white filled bar, control group fed a LC diet (n = 4, mean ± SEM). The significance of the changes in TLR4 protein expression levels is indicated (*p<0.05, **p<0.01).

https://doi.org/10.1371/journal.pone.0123942.g004

Discussion

Our study shows that during the feeding of HC diet, the up-regulated expression of immune related genes could be triggered by the increased translocation of LPS into the liver. Epigenetic modifications were involved in the regulation of the expression of immune related genes. These findings yield new insights into the contributions of epigenetic mechanism to hepatic inflammatory responses in dairy cows.

In recent decades, the practice of feeding HC diets to cows has been extensively applied to increase milk yields, thereby improving cost-efficiency. However, the consumption of a HC diet is harmful to the health of dairy cows [27]. Based on practical experience, two different diets, HC and LC diets, were designed and used in our present study. The main difference between these two diets was the different contents of NDF, ADF, and NFC, as indicated. The percentage of NDF and ADF in the nutritional composition was lower, but NFC was higher in the treatment diet than that in the control diet. Many studies have indicated that feeding this easily fermentable diet to cows is associated with a low rumen pH [8, 9]. In accord with previous studies, in our present study the rumen pH of cows fed a HC diet decreased linearly compared with that of cows fed a LC diet. A reduction in rumen pH during HC diet feeding caused the collapse of Gram-negative bacteria populations and resulted in increased free rumen LPS concentrations. In the study of Khafipour et al., replacing 21% of the dry matter of the total mixed ratio with pellets containing 50% wheat and 50% barley results in reduced milk yields and the percentage of milk fat, accompanied by reductions in rumen pH. The results of our study showed that milk yields increased in the first 4 weeks, and then declined in the subsequent 6 weeks in cows fed a HC diet compared with cows fed a LC diet, although the levels of crude protein and net energy were similar between these two diets. One reason for this finding could be that the buffering capacity of the rumen might counteract the decline of rumen pH in cows fed a HC diet in the short-term, whereas the extent of rumen pH depression could have exceeded the buffering capacity of the rumen, resulting in the decline of milk yield. Furthermore, in this and previous studies [6, 28], feeding of a HC diet induced the massive release of LPS in the rumen, which could trigger a local or systemic inflammatory response after the translocation of LPS into the bloodstream. In this setting, more nutrients would be used to sustain the inflammatory response, thereby reducing nutrient availability for milk synthesis. In a recent study, feeding a HC corn straw diet to cows elicited a pro-inflammatory response in the mammary gland and reduced milk yields [3]. Circulating LPS has been reported to depress the activity of key enzymes in milk component synthesis, such as acetyl-CoA carboxylase, lipoprotein lipase, and fatty acid synthase, and led to lesions in the mammary epithelia [10], thereby resulting in reduced milk yield.

Our data also showed that concentrations of LPS in the portal and hepatic veins were higher in cows fed a HC diet than that in cows fed a LC diet. Taking the concentrations of LPS in the rumen and portal veins into consideration, our findings indicated that a HC diet induced the translocation of LPS from the digestive tract into the bloodstream, and the amount of LPS delivered directly into the liver via the portal vein was higher in cows fed a HC diet compared with that in cows fed a LC diet. Recent experiments showed that during HC diet feeding, LPS produced in the digestive tract can contribute to local inflammatory responses in the rumen and intestine, thereby resulting in injury of the gastrointestinal barrier and facilitating the translocation of LPS and other hazardous substances into the circulating bloodstream [1, 29]. Although severe local inflammation in peripheral tissues might result in pro- and anti- inflammatory responses in the liver [12, 30], LPS in circulating blood, especially in the portal vein, is thought to be a major inducer of liver inflammation. Our measurements of the relative mRNA copy numbers by RT-qPCR confirmed that all of the liver samples from cows fed a HC diet showed signs of inflammation. Additionally, circulating LPS bound to LBP can be transferred to TLR4 located on the surface of innate immune cells, especially Kupffer cells (liver macrophages), to activate the intracellular TLR4-dependent signaling pathway, which results in the secretion of cytokines and chemokines [31]. We found that the feeding of a HC diet significantly increased the levels of mRNA encoding TLR4 and LBP, as well as the levels of TLR4 protein in the liver. These findings suggested that the increased translocation of LPS into the liver strongly activated the TLR4-dependent signaling pathway, which subsequently induced increased expression levels of pro- and anti-inflammatory cytokines (IL-1α, IL-1β, IL-6, and IL-10) and chemokines (IL-8, CCL5, and CCL20); such pro-inflammatory cytokines further induced the secretion of APPs by hepatocytes [32, 33]. Overall, our findings indicate that the liver of cows fed a HC diet are inflamed because of the excessive translocation of LPS into the liver via the portal vein.

The transcriptional start position (TSP) of each gene is a site that identifies the promoter region. The TSP of the gene encoding TLR4 was identified based on the results of the 5’-RACE experiment that we performed in a previous study, which showed that the predicted exon 1 in the database also acts as exon 1 in the liver. So, the 5’-adjacent region was considered to be its promoter. In the study of Marilynn et al., the TSP of SAA3 was also described [34]. For the genes encoding LBP and Hp, we identified that the TSP of both genes by blasting genomic DNA sequences with the respective mRNA sequences. The results of a MatInspector analysis of the promoters of LBP and Hp showed that both promoters indeed featured highly-conserved binding sites for the NF-κB (LBP) and C/EBP (LBP and Hp) transcription factors, which had previously been reported to be associated with expression of these genes [35, 36]. Moreover, an analysis of the promoter of the gene encoding aS1-casein was previously reported [20]. Finally, we confirmed the core promoter region of these genes, according to published data and our MatInspector analysis, and used them for further analysis by CHART and methylation.

Chromatin accessibility PCR was used to analyze the status of chromatin [23]. The local chromatin structure at promoters is a main factor that is highly associated with the amount of gene transcription. Our analysis of four immune related candidate genes firmly indicated that the loosening of chromatin at the promoter target region was significantly correlated with enhanced gene expression of immune genes. A previous study clearly demonstrated that the up-regulation of LAP expression could be attributed to chromatin decompaction at the promoter of bovine LAP, which facilitates NF-κB recruitment [21]. Additionally, a scatter analysis of the individual values of the expression of each candidate immune gene clearly demonstrated that chromatin loosening might be only one factor among several determinants that is important for inducing increased gene expression. Apart from chromatin remodeling, histone and DNA modifications are also necessary to create a proper chromatin environment to enhance or depress gene expression [37]. Furthermore, our data showed that the degree of chromatin compaction at the αS1-casein promoter was not different between the two groups fed different diets. Moreover, no significant correlation was detected between the degree of chromatin compaction at the αS1-casein promoter and its mRNA expression levels. These results showed that chromatin remodeling is highly related and specific to the expression level of immune-related genes. It is well-known that during sepsis, a severe illness caused by excessive endotoxin levels, chromatin remodeling is a key mechanism that drives the induction of endotoxin tolerance, which is induced to control excessive LPS-induced inflammation [18]. Herein, our study indicated that milder inflammation, which could be induced by circulating LPS derived from the digestive tract, could also trigger similar epigenetic mechanisms.

DNA methylation has been extensively studied and is thought to be one of the epigenetic mechanisms that can trigger chromatin compaction via methyl-CpG binding protein 2 (MeCP2) [38, 39]. Previously, the regulation of TLR2 and TLR4 expression in inflammation was shown to be associated with DNA methylation [40, 41]. In this present study, single HpaII restriction sites at the promoters of five candidate genes were analyzed to characterize the interdependence of CpG methylation and chromatin remodeling. Our data demonstrated that this site at the promoter of four immune relevant candidate genes was indeed hypomethylated in liver samples from cows fed a HC diet compared with those in liver samples from cows fed a LC diet. This finding indicated that after feeding a HC diet to cows, demethylation at the promoter of immune-related genes resulted in the loosening of chromatin at the corresponding promoter regions of those genes, thereby enhancing the expression of these genes. An earlier study showed that promoter hypomethylation can elicit histone hyperacetylation and release of the transcriptional repressor MeCP2, resulting in the activation of transcription [42]. Several studies have been preformed to examine the roles of epigenetic markers, such as DNA methylation and histone acetylation, on the regulation of LPS-induced innate immune gene expression [43, 44]. Additionally, the increased levels of LPS in the bloodstream of dairy cows fed a high proportion of concentrated diet have been well documented [3, 6]. In this present study, the content of LPS in the portal vein was significantly increased during feeding of a HC diet. Hence, the loosening of chromatin and hypomethylation at the promoter of immune-related genes in the liver of cows fed a HC diet could be attributed to increased translocation of LPS into the liver via the portal vein. Furthermore, the difference in the percentage of methylation at the αS1-casein promoter between the cows fed the two different diets was not significant. Thus, these data indicated that the expression of αS1-casein in the liver was not affected by the increased translocation of LPS into the liver, which highlights the sensitivity of immune-related genes to LPS.

Herein, we focused on the contribution of chromatin decompaction and DNA demethylation to the expression of immune genes. Further studies will be needed to elucidate the epigenetic regulation of histone modification and microRNAs. Considering all of the diverse modulators of chromatin remodeling, such as histone deacetylases [45, 46] or vitamin D3 [47, 48], these data might yield novel ways to regulate systemic inflammation, which can be induced by LPS that translocates from the digestive tract into the blood, and to improve high milk performance during productive practice of feeding a HC diet to cows.

In summary, our data indicate that feeding a HC diet to cows depresses the rumen pH value, results in the increase of LPS concentration in rumen. It has been previously reported that the feeding of HC diet caused the disruption of gastrointestinal barrier [1], this contributes to the translocation of LPS from digestive tract into portal vein. Hence, the increased translocation of LPS into the liver via the portal vein enhances the hepatic expression of immune-related genes, which are mediated by epigenetic mechanisms involving in chromatin decompacion and DNA demethylation. Our data confirmedly showed that chromatin decompaction and DNA demethylation in relevant areas of the promoter of candidate immune genes are strongly correlated with the enhanced expression of those immune genes. These findings demonstrate that epigenetic mechanisms contribute to the expression of immune-related genes in the livers of dairy cows fed a HC diet.

Supporting Information

S1 Fig. Gel images of promoter methylation.

https://doi.org/10.1371/journal.pone.0123942.s001

(TIF)

S3 Table. Primers for methylation analysis.

https://doi.org/10.1371/journal.pone.0123942.s004

(DOCX)

Acknowledgments

We thank Prof. Hans-Martin Seyfert (Dept. of Genome Biology, Leibniz Institute for Farm Animal Biology, Germany) for his valuable advice on the epigenetic analysis.

Author Contributions

Conceived and designed the experiments: XS SZ. Performed the experiments: GC KZ TX DJ JG. Analyzed the data: GC. Contributed reagents/materials/analysis tools: XS. Wrote the paper: GC.

References

  1. 1. Liu JH, Xu TT, Zhu WY, Mao SY (2014) A high-grain diet alters the omasal epithelial structure and expression of tight junction proteins in a goat model. Vet J 201:95–100. pmid:24906500
  2. 2. Tao S, Duanmu Y, Dong H, Tian J, Ni Y, Zhao R (2014) A high-concentrate diet induced colonic epithelial barrier disruption is associated with the activating of cell apoptosis in lactating goats. BMC Vet Res 10:235. pmid:25256013
  3. 3. Zhou J, Dong GZ, Ao CJ, Zhang S, Qiu M, Wang X, et al. (2014) Feeding a high-concentrate corn straw diet increased the release of endotoxin in the rumen and pro-inflammatory cytokines in the mammary gland of dairy cows. BMC Vet Res 10:10. pmid:24410882
  4. 4. Murata H, Shimada N, Yoshioka M (2004) Current research on acute phase proteins in veterinary diagnosis: an overview. Vet J 168:28–40. pmid:15158206
  5. 5. Jourdan T, Djaouti L, Demizieux L, Gresti J, Verges B, Degrace P (2009) Liver Carbohydrate and Lipid Metabolism of Insulin-Deficient Mice Is Altered by trans-10, cis-12 Conjugated Linoleic Acid. J Nutr 139:1901–1907. pmid:19692528
  6. 6. Khafipour E, Krause DO, Plaizier JC (2009) A grain-based subacute ruminal acidosis challenge causes translocation of lipopolysaccharide and triggers inflammation. J Dairy Sci 92:1060–1070. pmid:19233799
  7. 7. Gozho GN, Plaizier JC, Krause DO, Kennedy AD, Wittenberg KM (2005) Subacute ruminal acidosis induces ruminal lipopolysaccharide endotoxin release and triggers an inflammatory response. J Dairy Sci 88:1399–1403. pmid:15778308
  8. 8. Gozho GN, Krause DO, Plaizier JC (2007) Ruminal lipopolysaccharide concentration and inflammatory response during grain-induced subacute ruminal acidosis in dairy cows. J Dairy Sci 90:856–866. pmid:17235162
  9. 9. Emmanuel DGV, Dunn SM, Ametaj BN (2008) Feeding high proportions of barley grain stimulates an inflammatory response in dairy cows. J Dairy Sci 91:606–14. pmid:18218747
  10. 10. Dong GZ, Liu SM, Wu YX, Lei CL, Zhou J, Zhang S (2011) Diet-induced bacterial immunogens in the gastrointestinal tract of dairy cows: Impacts on immunity and metabolism. Acta Vet Scand 53:48. pmid:21824438
  11. 11. Liaskou E, Wilson DV, Oo YH (2012) Innate immune cells in liver inflammation. Mediat Inflamm 21.
  12. 12. Vels L, Rontved CM, Bjerring M, Ingvartsen KL (2009) Cytokine and acute phase protein gene expression in repeated liver biopsies of dairy cows with a lipopolysaccharide-induced mastitis. J Dairy Sci92:922–934. pmid:19233785
  13. 13. Mimura Y, Sakisaka S, Harada M, Sata M, Tanikawa K (1995) Role of hepatocytes in direct clearance of lipopolysaccharide in rats. Gastroenterology 109:1969–1976. pmid:7498663
  14. 14. Bieghs V, Trautwein C (2013) The innate immune response during liver inflammation and metabolic disease. Trends Immunol 34:446–452. pmid:23668977
  15. 15. Saturnino SF, Prado RO, Cunha-Melo JR, Andrade MV (2010) Endotoxin tolerance and cross-tolerance in mast cells involves TLR4, TLR2 and Fc epsilon R1 interactions and SOCS expression: perspectives on immunomodulation in infectious and allergic diseases. BMC Infect Dis 10:8. pmid:20067628
  16. 16. Foster SL, Hargreaves DC, Medzhitov R (2007) Gene-specific control of inflammation by TLR-induced chromatin modifications. Nature 447:972–978. pmid:17538624
  17. 17. Park SH, Park-Min KH, Chen J, Hu XY, Ivashkiv LB (2011) Tumor necrosis factor induces GSK3 kinase-mediated cross-tolerance to endotoxin in macrophages. Nat Immunol 12:607–615. pmid:21602809
  18. 18. El Gazzar M, Liu TF, Yoza BK, McCall CE (2010) Dynamic and selective nucleosome repositioning during endotoxin tolerance. J Biol Chem 285:1259–1271. pmid:19901031
  19. 19. Dong GZ, Qiu M, Ao CJ, Zhou J, Khas E, Wang X (2014) Feeding a high-concentrate corn straw diet induced epigenetic alterations in the mammary tissue of dairy cows. Plos One 9: e107659. pmid:25222274
  20. 20. Vanselow J, Yang W, Herrmann J, Zerbe H, Schuberth HJ, Petzl W, et al. (2006) DNA-remethylation around a STAT5-binding enhancer in the alpha S1-casein promoter is associated with abrupt shutdown of alpha S1-casein synthesis during acute mastitis. J Mol Endocrinol 37:463–477. pmid:17170087
  21. 21. Liu SZ, Shi XM, Bauer I, Gunther J, Seyfert HM (2011) Lingual antimicrobial peptide and IL-8 expression are oppositely regulated by the antagonistic effects of NF-kappa B p65 and C/EBP beta in mammary epithelial cells. Mol Immunol 48:895–908. pmid:21255844
  22. 22. Dong HB, Wang SQ, Jia YY, Ni YD, Zhang YS, Zhuang S, et al. (2013) Long-term effects of subacute ruminal acidosis (SARA) on milk quality and hepatic gene expression in lactating goats fed a high-concentrate diet. Plos One 8: e82850. pmid:24376594
  23. 23. Rao SH, Procko E, Shannon MF (2001) Chromatin remodeling, measured by a novel real-time polymerase chain reaction assay, across the proximal promoter region of the IL-2 gene. J Immunol 167:4494–4503. pmid:11591776
  24. 24. Chang GJ, Petzl W, Vanselow J, Günther J, Shen XZ, Seyfert HM (2015) Epigenetic mechanisms contribute to enhanced expression of immune response genes in the liver of cows after experimentally induced escherichia coli mastitis. The Veterinary Journal In Press.
  25. 25. De Oliveira NFP, Andia DC, Planello AC, Pasetto S, Marques MR, Nociti FH, et al. (2011) TLR2 and TLR4 gene promoter methylation status during chronic periodontitis. J Clin Periodontol 38:975–83. pmid:21899586
  26. 26. de Faria Amormino SA, Arao TC, Saraiva AM, Gomez RS, Dutra WO, da Costa JE, et al. (2013) Hypermethylation and low transcription of TLR2 gene in chronic periodontitis. Hum Immunol 74:1231–1236. pmid:23747679
  27. 27. Plaizier JC, Khafipour E, Li S, Gozho GN, Krause DO (2012) Subacute ruminal acidosis (SARA), endotoxins and health consequences. Anim Feed Sci Tech 172:9–21.
  28. 28. Zebeli Q, Ametaj BN (2009) Relationships between rumen lipopolysaccharide and mediators of inflammatory response with milk fat production and efficiency in dairy cows. J Dairy Sci 92:3800–3809. pmid:19620662
  29. 29. Liu JH, Xu TT, Zhu WY, Mao SY (2014) High-grain feeding alters caecal bacterial microbiota composition and fermentation and results in caecal mucosal injury in goats. Brit J Nutr 112:416–427. pmid:24846282
  30. 30. Zhong J, Deaciuc JV, Burikhanov R, de Villiers WJS (2006) Lipopolysaccharide-induced liver apoptosis is increased in interleukin-10 knockout mice. Bba-mol Basis Dis 1762:468–477. pmid:16497487
  31. 31. Takeuchi O, Akira S (2010) Pattern Recognition Receptors and Inflammation. Cell 140:805–820. pmid:20303872
  32. 32. Yoshioka M, Ito T, Miyazaki S, Nakajima Y (1998) The release of tumor necrosis factor-alpha, interleukin-1, interleukin-6 and prostaglandin E-2 in bovine Kupffer cells stimulated with bacterial lipopolysaccharide. Vet Immunol Immunop 66:301–307. pmid:9880106
  33. 33. Jensen LE, Whitehead AS (1998) Regulation of serum amyloid A protein expression during the acute-phase response. Biochem J 334:489–503. pmid:9729453
  34. 34. Larson MA, Weber A, McDonald TL (2006) Bovine serum amyloid A3 gene structure and promoter analysis: Induced transcriptional expression by bacterial components and the hormone prolactin. Gene 380:104–110. pmid:16837143
  35. 35. Arambasic J, Poznanovic G, Ivanovic-Matic S, Bogojevic D, Mihailovic M, Uskokovic A, et al. (2010) Association of the glucocorticoid receptor with STAT3, C/EBP beta, and the hormone-responsive element within the rat haptoglobin gene promoter during the acute phase response. Iubmb Life 62:227–236. pmid:20155810
  36. 36. Kirschning CJ, Unbehaun A, Fiedler G, Hallatschek W, Lamping N, Pfeil D, et al. (1997) The transcriptional activation pattern of lipopolysaccharide binding protein (LBP) involving transcription factors AP-1 and C/EBP beta. Immunobiology 198:124–135. pmid:9442384
  37. 37. Bayarsaihan D (2011) Epigenetic mechanisms in inflammation. J Dent Res 90:9–17. pmid:21178119
  38. 38. Chahrour M, Jung SY, Shaw C, Zhou XB, Wong STC, Qin J, et al. (2008) MeCP2, a key contributor to neurological disease, activates and represses transcription. Science 320:1224–1229. pmid:18511691
  39. 39. Martinowich K, Hattori D, Wu H, Fouse S, He F, Hu Y, et al. (2003) DNA methylation-related chromatin remodeling in activity-dependent Bdnf gene regulation. Science 302:890–893. pmid:14593184
  40. 40. Shuto T, Furuta T, Oba M, Xu H, Li J, Cheung J, et al. (2006) Promoter hypomethylation of Toll-like receptor-2 gene is associated with increased proinflammatory response toward bacterial peptidoglycan in cystic fibrosis bronchial epithelial cells. Faseb J 20:782–800. pmid:16478769
  41. 41. Takahashi K, Sugi Y, Hosono A, Kaminogawa S (2009) Epigenetic regulation of TLR4 gene expression in intestinal epithelial cells for the maintenance of intestinal homeostasis. J Immunol 183:6522–6529. pmid:19846881
  42. 42. Kudo S (1998) Methyl-CpG-binding protein MeCP2 represses Sp1-activated transcription of the human leukosialin gene when the promoter is methylated. Mol Cell Biol 18:5492–5499. pmid:9710633
  43. 43. Doherty R, O'Farrelly C, Meade KG (2013) Epigenetic regulation of the innate immune response to LPS in bovine peripheral blood mononuclear cells (PBMC). Vet Immunol Immunop 154:102–110. pmid:23764468
  44. 44. Green BB, Kerr DE (2014) Epigenetic contribution to individual variation in response to lipopolysaccharide in bovine dermal fibroblasts. Vet Immunol Immunop 157:49–58. pmid:24268632
  45. 45. Villagra A, Sotomayor EM, Seto E (2010) Histone deacetylases and the immunological network: implications in cancer and inflammation. Oncogene 29:157–173. pmid:19855430
  46. 46. Delcuve GP, Khan DH, Davie JR (2012) Roles of histone deacetylases in epigenetic regulation: emerging paradigms from studies with inhibitors. Clin Epigenetics 4:13. pmid:22938091
  47. 47. Pereira F, Barbachano A, Silva J, Bonilla F, Campbell MJ, Munoz A, et al. (2011) KDM6B/JMJD3 histone demethylase is induced by vitamin D and modulates its effects in colon cancer cells. Hum Mol Genet 20:4655–4665. pmid:21890490
  48. 48. Pereira F, Barbachano A, Singh PK, Campbell MJ, Munoz A, Larriba MJ (2012) Vitamin D has wide regulatory effects on histone demethylase genes. Cell Cycle 11:1081–1089. pmid:22370479