1. Introduction
Nonalcoholic fatty liver disease (NAFLD) is recognized as the most common chronic liver disease with an increasing prevalence, and is strongly associated with obesity, insulin resistance, cardiovascular disease, hypertriglyceridemia, and hyperuricemia [
1]. NAFLD encompasses a broad spectrum ranging from mild hepatic steatosis to more severe forms such as nonalcoholic hepatitis (NASH), hepatic cirrhosis, and hepatocarcinoma [
1]. Hepatic steatosis, defined as triglyceride (TG) accumulation in the cytoplasm of hepatocytes, is a driving force for NAFLD [
1]. Hepatic steatosis arises from the imbalance between TG acquisition involving de novo lipogenesis (DNL) and uptake of fatty acid (FA) from the adipose tissue, and TG removal involving β-oxidation of FA and export of very low-density lipoprotein from the liver [
1]. NASH develops from hepatic steatosis by a combination of TG accumulation, apoptosis, oxidative stress, inflammation, and fibrosis, leading to hepatic cirrhosis, hepatocarcinoma, and ultimately liver-related mortality [
1].
While various factors, including genetic disposition, contribute to the development of NAFLD, alteration of the epigenome by environmental factors can play a key role in this event. Several studies have reported the link between the regulation of hepatic steatosis and histone modifications [
2,
3,
4]. Histone deacetylase 3 (HDAC3) is an epigenetic regulator associated with hepatic steatosis [
2]. HDAC3 downregulates the expression of genes related to fatty acid uptake, TG synthesis, and lipolysis, thus preventing hepatic steatosis [
2]. Furthermore, the histone H3K4 methyltransferase MLL3/4 upregulates proliferator-activated receptor gamma isoform 2 (PPARγ2) and its target steatosis genes, thereby inducing hepatic steatosis [
3]. The histone demethylase plant homeodomain finger 2 (Phf2) was recently shown to induce hepatic steatosis through the upregulation of carbohydrate-responsive element-binding protein (ChREBP) target genes [
4]. Our previous study demonstrated that the histone H3K9 demethylase JMJD2B stimulated the expression of PPARγ2 and its target genes and contributed to hepatic steatosis [
5]. However, the epigenetic regulation of hepatic steatosis through histone methylation is still poorly understood.
The histone lysine demethylase KDM7A (also known as JHDM1D) is a plant homeodomain (PHD) finger protein (PHF) and Jumonji (Jmj)-containing histone demethylase [
6]. KDM7A is a dual histone demethylase that removes repressive di-methylation marks at lysine 9 of histone H3 (H3K9me2) and lysine 27 of histone H3 (H3K27me2), converting the marks to a mono-methylated state; thus, KDM7A functions as a transcription activator [
6]. Reportedly, KDM7A promotes neural differentiation through the upregulation of CCAAT-enhancer-binding protein β and fibroblast growth factor 4 expression [
7], and balanced adipogenic and osteogenic differentiation from progenitor cells by controlling Wnt signaling [
8]. Furthermore, KDM7A is involved in promoting prostate cancer [
9]. KDM7A was recently shown to control inflammatory responses through the regulation of nuclear factor kappa B (NF-κB)-dependent genes [
10]. However, the role of KDM7A in the development of hepatic steatosis has not been elucidated.
Two diacylglycerol acyltransferases (DGATs), namely, DGAT1 and DGAT2, are known to play an important role in hepatic TG synthesis and induce hepatic steatosis through the catalysis of the final step of TG synthesis [
11]. In the liver, these two enzymes differentially use exogenously or endogenously synthesized fatty acids. While DGAT1 preferentially utilizes exogenously fatty acids for TG synthesis, DGAT2 consumes fatty acids from DNL and increases in NAFLD. Using RNA-sequencing study, we found that the expression of DGAT2 was increased in KDM7A-overexpressing AML12 cells. Accordingly, we hypothesized that KDM7A epigenetically upregulates DGAT2 expression by erasing the transcriptionally repressive histone marks H3K9me2 and H3K27me2 on the promoter of DGAT2, eventually leading to hepatic steatosis through increasing TG synthesis.
In the present study, we investigated whether KDM7A plays a role in the development of hepatic steatosis through the upregulation of DGAT2 in vitro and in vivo. We measured intracellular TG accumulation and DGAT2 expression in KDM7A-overexpressing or knockdown AML12 cells. In addition, we examined the enrichment of repressive histone marks H3K9me2 and H3K27me2 on the promoter of DGAT2 in KDM7A-overexpressing AML12 cells and assessed the development of hepatic steatosis in vivo in KDM7A-overexpressing mice.
3. Discussion
Epigenetic modifiers can contribute to the development and progression of several metabolic diseases, including NAFLD [
12]. Understanding of the epigenetic mechanisms involved in NAFLD progression may reveal therapeutic targets; therefore, identifying epigenetic modifiers is imperative to make progress in the treatment of NAFLD. Aberrations in histone methylation are known to correlate with the development of hepatic steatosis [
2,
3,
4,
5]. However, histone methylation underlying NAFLD pathogenesis is poorly understood. In the current study, we identified a novel histone methylation-mediated mechanism involved in the development of hepatic steatosis.
Histone lysine demethylase KDM7A is a recently identified member of the PHD and JmjC domain-containing histone demethylase family [
6]. It demethylates both H3K9me2 and H3K27me2, converting transcriptionally repressive histone marks to active histone marks. KDM7A was shown to regulate brain development, cell cycle, cell proliferation, and inflammation response [
7,
8,
9]. However, its role in hepatic steatosis remains unclear. We used gain-of-function and loss-of-function studies in AML12 cells to determine whether KDM7A induces hepatic steatosis. Adenovirus-mediated KDM7A overexpression resulted in an increase in intracellular TG accumulation in AML12 cells, whereas shRNA-mediated KDM7A knockdown significantly decreased intracellular TG levels in AML12 cells. These findings suggest that KDM7A plays an important role in the development of hepatic steatosis. Then, we attempted to identify the target gene regulated by KDM7A by performing RNA-sequencing analysis. RNA-sequencing analysis revealed the upregulation in the expression of the genes associated with TG synthesis, including perilipin 5, monoacyl-transferase (MOAT), and DGAT1/2. Among these genes, DGATs catalyze the final step in the synthesis of TG through the esterification of diacylglycerol (DAG) with a fatty acid. Mammals have two isozymes of DGAT, DGAT1 and DGAT2 which have no sequence homology to each other. DGAT1 is expressed abundantly in the intestine and plays a role in the absorption of dietary lipids by reconstituting TG in a committed step. DGAT2 is mainly located in adipose tissue and liver and plays a central role in intracellular TG accumulation through esterification of endogenous fatty acid to DAG. Several reports have demonstrated that DGAT2 contributes to the development of hepatic steatosis [
13,
14,
15,
16]. DGAT2 antisense-treated mice and hepatic-specific DGAT2 knockout mice have lower TG accumulation and ameliorated hepatic steatosis [
12,
13,
14]. Therefore, we believe that DGAT2 may be a potential downstream target gene of KDM7A necessary for the development of hepatic steatosis. In the present study, we observed that KDM7A overexpression increased DGAT2 expression and intracellular TG accumulation in AML12 cells. Enzymatic inhibition of DGAT2 using PF-06427878 significantly blocked the KDM7A-mediated-increase in TG accumulation, suggesting that DGAT2 is a downstream target gene of KDM7A. In contrast, KDM7A knockdown resulted in decreased DGAT2 expression and TG accumulation in AML12 cells and significantly reduced FFAs-induced intracellular TG accumulation. Together, these results indicate that KDM7A promotes DGAT2 expression, leading to increased intracellular TG accumulation, and consequently results in the development of hepatic steatosis.
We ascertained the functional role of KDM7A in the pathogenesis of hepatic steatosis in vivo. We found that mice injected with recombinant Ad-KDM7A showed elevated hepatic TG accumulation as revealed from the quantification of TG levels and ORO and H&E staining. Thus, KDM7A contributes to the development of hepatic steatosis. Furthermore, we confirmed KDM7A-DGAT2 signaling in the liver tissue of Ad-KDM7A-injected mice. Consistent with the results in AML12 cells, the expression of DGAT2 was higher in the liver of Ad-KDM7A-injected mice than in the liver of Ad-GFP-injected mice. Therefore, the in vivo data show that KDM7A is an epigenetic regulator involved in the development of hepatic steatosis through the upregulation of DGAT2 expression.
NASH develops from hepatic steatosis by a combination of TG accumulation, apoptosis, oxidative stress, and inflammation [
1]. It is uncertain whether the KDM7A-induced hepatic steatosis shown in the present study favors the development of NASH and liver injuries. Hepatic steatosis is a driving force for NAFLD because increased intracellular TG accumulation induces inflammation, oxidative stress and insulin resistance by generating lipotoxic intermediates, which leads to NASH, fibrosis and cirrhosis [
2]. On the other hand, hepatic TG synthesis from toxic saturated fatty acids prevents NAFLD development [
4]. Recently, it was reported that overexpression of histone demethylase Phf2, which belongs to the KDM7 histone demethylase family, favored hepatic steatosis development, but protected the liver from inflammation and oxidative stress by desaturation of detrimental saturated fatty acid into DAG and TG, and fibrosis development [
4].
Furthermore, the contradictory role of DGAT2 during hepatic steatosis progression to NASH has been reported [
13,
14,
15]. Amin et al. showed that DGAT2 inhibition prevented the development of hepatic steatosis and steatohepatitis because it inhibited TG synthesis [
13]. However, other studies reported opposite results [
13,
14]. Hepatic deletion of DGAT2 reduced hepatic steatosis but showed no increase in inflammation and fibrosis in mice [
14]. Further, DGAT2 deficiency using antisense oligo improved hepatic steatosis but exacerbated liver damage and hepatic fibrosis [
15]. Therefore, to clarify whether KDM7A-DGAT2 signaling-induced hepatic steatosis exacerbates or prevents the development of NASH and fibrosis, further studies should be warranted.
We determined the mechanism underlying KDM7A-mediated upregulation in DGAT2 expression. Di-methylation of H3K9 (H3K9me2) and H3K27 (H3K27me2) is associated with transcriptional repression of genes. To characterize the mechanism by which KDM7A promotes DGAT2 expression, we examined the changes in H3K9me2 and H3K27me2 levels on the DGAT2 promoter in KDM7A-overexpressing AML12 cells using ChIP-PCR. The recruitment of KDM7A to the promoter of DGAT2 was increased, whereas the enrichment of H3K9me2 and H3K27me2 was reduced in KDM7A-overexpressing AML12 cells. These results suggest that KDM7A erases repressive histone marks H3K9me2 and H3K27me2 on the promoter of DGAT2, resulting in the stimulation of DGAT2 expression.
Based on these findings, we propose an epigenetic model of KDM7A associated with the development of hepatic steatosis. Under the normal state, di-methylation of H3K9 (H3K9me2) and H3K27 (H3K27me2) on the promoter of DGAT2 represses DGAT2 expression and reduces hepatic TG accumulation. However, under hepatic steatosis states such as obesity, KDM7A is overexpressed and erases repressive H3K9me2 and H3K27me2 on the promoter of DGAT2, which leads to increased TG accumulation and results in hepatic steatosis (
Figure 6). Accordingly, the present study provides a new insight into the potential role of epigenetic regulation of hepatic steatosis.
4. Materials and Methods
4.1. Reagents
Dulbecco’s modified Eagle’s medium/Ham’s Nutrient Mixture F-12 1:1(DMEM/F-12), penicillin–streptomycin, and fetal bovine serum were obtained from HyClone Laboratories Inc. (Logan, UT, USA). Oleic acid, palmitic acid and PF-06427878 were purchased from Sigma-Aldrich (St. Louis, MO, USA). The antibody against KDM7A was purchased from Aviva systems biology (San Diego, CA, USA). Antibodies against H3K9me2 and H3K27me2 were obtained from Millipore (Billerica, MA, USA). Antibodies against DGAT2 was obtained from Novus Biologicals (Littleton, CO, USA) and β-actin was supplied by Santa Cruz Biotechnology (Santa Cruz, CA, USA).
4.2. Cell Culture
HepG2 and AML12 cells were obtained from the American Type Culture Collection (Manassas, VA, USA) and cultured in DMEM/F12 1:1 supplemented with 10% heat-inactivated FBS, 1X Insulin-Transferrin-Selenium (GibcoTM, Paisley, UK), 40 ng/mL dexamethasone, 20 U/mL penicillin, and 20 μg/mL streptomycin. HepG2 cells and AML12 cells were treated with 1 mM a mixture of free fatty acids (FFAs) (oleic acid: palmitic acid, 2:1) for 24 h to induce hepatic steatosis cell model.
4.3. Infection of AML12 Cells with Recombinant Adenovirus
Adenovirus vectors encoding green fluorescent protein (GFP), KDM7A (Ad-GFP, Ad-KDM7A), and containing scramble short-hairpin RNA (shRNA) and KDM7A shRNA (Ad-GFP shRNA, Ad-KDM7A shRNA) were purchased from Vector Biolabs (Malvern, PA, USA). AML12 cells were infected with Ad-GFP, Ad-KDM7A, Ad-GFP shRNA or Ad-KDM7A shRNA at a multiplicity of infection (MOI) ranging from 1 to 100 PFU/cell. Thereafter, the cells were washed and incubated with fresh medium. At 48 h post-infection, cells were harvested and analyzed.
4.4. Total RNA Preparation and Real-Time Quantitative Polymerase Chain Reaction (qPCR)
Total RNA was extracted using TRIzol (Invitrogen, Carlsbad, CA, USA) according to the manufacturer’s instructions. cDNA was generated from 1 μg of total RNA using the GoScript Reverse Transcription System (Promega, Madison, WI, USA) according to the manufacturer’s protocol. Real-time qPCR was performed using a SYBR Green premixed Taq reaction mixture with gene-specific primers; KDM7A sense, 5′-GAAATCAATCAGAAGGCACAAAG-3′; KDM7A antisense, 5′-GATTCAGCGTGTACTTTCCATTC-3′, and DGAT2 sense, 5′-AACCTGCTGACCACCAGGAACTAT-3′; DGAT2 antisense, 5′-AGGGCCTTATGCCAGGAAACTTCT-3′.
4.5. Western Blot Analysis
Equal amounts of proteins (20 μg/lane) from cell lysates or liver tissues were resolved on 10% sodium dodecyl sulfate polyacrylamide gel and then transferred onto polyvinylidene difluoride membranes (Millipore, Billerica, MA, USA). The membranes were blocked in 5% nonfat skim milk and probed with primary antibodies against KDM7A, H3K9me2, H3K27me2, DGAT2, and β-actin. After washing with Tween 20/Tris-buffered saline, the membranes were incubated with a horseradish peroxidase-conjugated secondary antibody (1:1000) for 1 h at room temperature. The proteins were detected using an enhanced chemiluminescence (ECL) western blot detection kit (Amersham, Uppsala, Sweden).
4.6. Chromatin Immunoprecipitation (ChIP)-qPCR
Briefly, AML12 cells were fixed with 1% formaldehyde for 10 min at room temperature. The crosslinked chromatin was sonicated to shear into 400-bp fragments using a Bioruptor sonicator (Diagenode, Denville, NJ, USA). Samples were immunoprecipitated using 1–2 μg antibodies against KDM7A, H3K9me2, H3K27me2, or nonspecific IgG control in the presence of a secondary antibody conjugated to Dynabeads (Invitrogen, Carlsbad, CA, USA). Purified DNA was subjected to qPCR using the following primers: DGAT2 sense, 5-GTCCAGCTGTAACCTAACGA-3; DGAT2 antisense, 5-AAGCCTCTCTCTCAGACGTT-3. ChIP data were normalized as per control IgG or expressed as a percentage of input.
4.7. Animal Experiments
C57BL/6 mice (8 weeks of age) were purchased from Central Lab Animal Inc. (Seoul, Korea) and maintained on a regular chow diet (ND: 10 kcal% fat, SLACOM, #M01). The mice were injected with a total of 1 × 109 PFU recombinant adenoviruses (Ad-GFP or Ad-KDM7A) via tail vein injection. After injection, adenovirus-injected mice were fed a ND for 4 weeks. Four weeks after injection, the mice were fasted for 6 h, and blood was harvested from the tail vein. The animals were then sacrificed, and liver tissues were collected for further analysis. The animal protocol used in this study was reviewed and approved by the Pusan National University’s Institutional Animal Care and Use Committee in accordance with established ethical and scientific care procedures (PNU-2020-2718).
4.8. Histological Analysis
Liver tissues isolated from the mice were dissected and fixed in 10% buffered formalin. Fixed tissues were embedded in paraffin and sectioned at a 5 μm thickness using a frozen microtome (HM560H, Microm Laboratory, Walldorf, Germany). The liver sections were stained with hematoxylin and eosin (H&E) and Oil Red O (ORO) and subjected to microscopic observation.
4.9. Hepatic TG Measurement
Hepatic lipids were extracted from cell extracts and liver tissues according to the following procedure. Briefly, liver tissues were homogenized in a chloroform-methanol solution (2:1 v/v). The homogenates were incubated for 1 h at room temperature and centrifuged at 3000 rpm for 10 min. The bottom layer (organic phase) was removed and dried overnight. After dissolving the sample in ethanol, hepatic TG level was determined using an AM 157S-K TG kit (Asan Pharmaceuticals Co., Ltd., Seoul, Korea) and normalized to the protein concentration.
4.10. Blood Biochemistry Analysis
After fasting for 6 h, the mice were sacrificed, and blood was collected. The samples were centrifuged at 1000× g for 15 min at 4 °C to obtain the serum and stored at −80 °C until analysis. The levels of TG, free fatty acids and total cholesterol in serum were determined using a DOGEN Kit (DoGenBio. Co. Ltd., Seoul, Korea).
4.11. Statistical Analysis
Data are expressed as the mean ± standard deviation (SD). Statistically significant differences were determined using a two-tailed Student’s t-test. For all statistical analyses, values of p < 0.05 were considered significant.