Hepatic Transcriptome Analysis Revealing the Molecular Pathogenesis of Type 2 Diabetes Mellitus in Zucker Diabetic Fatty Rats

Around 9% of the adult population in the world (463 million) suffer from diabetes mellitus. Most of them (~90%) belong to type 2 diabetes mellitus (T2DM), which is a common chronic metabolic disorder, and the number of cases has been reported to increase each year. Zucker diabetic fatty (ZDF) rat provides a successful animal model to study the pathogenesis of T2DM. Although previous hepatic transcriptome studies revealed some novel genes associated with the occurrence and development of T2DM, there still lacks the comprehensive transcriptomic analysis for the liver tissues of ZDF rats. We performed comparative transcriptome analyses between the liver tissues of ZDF rats and healthy ZCL rats and also evaluated several clinical indices. We could identify 214 and 104 differentially expressed genes (DEGs) and lncRNAs in ZDF rats, respectively. Pathway and biofunction analyses showed a synergistic effect between mRNAs and lncRNAs. By comprehensively analyzing transcriptomic data and clinical indices, we detected some typical features of T2DM in ZDF rats, such as upregulated metabolism (significant increased lipid absorption/transport/utilization, gluconeogenesis, and protein hydrolysis), increased inflammation, liver injury and increased endoplasmic reticulum (ER) stress. In addition, of the 214 DEGs, 114 were known and 100 were putative T2DM-related genes, most of which have been associated with substance metabolism (particularly degradation), inflammation, liver injury and ER stress biofunctions. Our study provides an important reference and improves understanding of molecular pathogenesis of obesity-associated T2DM. Our data can also be used to identify potential diagnostic markers and therapeutic targets, which should strengthen the prevention and treatment of T2DM.

Introduction

Type 2 diabetes mellitus (T2DM) is a very common chronic metabolic disorder and mainly arises from insulin resistance, which accounts for about 90% of diabetes mellitus (DM) patients (1–4). According to International Diabetes Federation, in 2019, there were around 463 million patients with DM across the globe, corresponding to 9.09% of the adult population in the world (2). China ranks first among all countries, about 116.4 million patients with DM, accounting for approximately one-tenth of the Chinese population. T2DM is often accompanied by various forms of liver injury (5, 6), which facilitates the deterioration of T2DM (7). T2DM and liver injury interact with each other to form a vicious circle.

To study T2DM pathogenesis, many animal models have been successfully established (7). Among them, male Zucker diabetic fatty (ZDF) rats are widely considered to be a successful animal model of diabetic obesity; this model have been extensively applied in the research of T2DM as it aptly simulates the occurrence and development of obesity-associated T2DM (8). Male ZDF rats carry a spontaneous homozygous mutation in the leptin receptor gene; thus, the symptoms of obesity and insulin resistance appear in them at a young age, gradually developing into hyperglycemia and T2DM with age.

Some ZDF rat model-based “omics” studies on T2DM have been reported in recent years (9, 10). Suh YH et al. reported the downregulation of an antioxidant defense gene cluster and cytochrome P450 gene in the adipose and liver tissues of ZDF rats using a cDNA chip (11); Nawano M et al. performed transcriptomic analyses and reported the downregulation of some glucose utilization-related genes and upregulation of some gluconeogenic genes in the skeletal muscle tissue of ZDF rats (12). Slieker LJ et al. also used the ZDF rat model and found that the expression of two important glucose transporter proteins (GLUT2 in islet beta cells and GLUT4 in skeletal muscle cells) were downregulated (13).

Although earlier transcriptomic studies have used the ZDF rat model to discover some novel genes associated with T2DM occurrence and development, comprehensive transcriptomic analyses of the liver tissues of ZDF rats have not yet been conducted. Herein we performed comparative transcriptome analyses of the liver tissues of ZDF and Zucker control lean (ZCL) rats, in addition to clinical indices detection. In total, 214 differentially expressed genes (DEGs) and 104 differentially expressed lncRNAs were identified. Pathway and biofunction analyses revealed a synergistic effect between mRNAs and lncRNAs in the liver tissues of ZDF rats. The 214 DEGs included 114 known and 100 putative T2DM-related genes. Our results could provide an important reference and improve the understanding of the molecular pathogenesis of obesity-associated T2DM.

Materials and Methods

Animals and Tissues

A group of male ZDF/Gmi-fa/fa rats (n = 5) and their gender- and age-matched ZCL (ZCL/Gmi-+/fa) rats (n = 5) were purchased at 6 weeks of age from Beijing Vital River Laboratory Animal Technology Co. Ltd. (Beijing, China). The animals were maintained from weeks 6 to 22 in a P3 room with a 12:12-h light/dark cycle at 20°C–25°C and 60% ± 5% atmospheric humidity. The study was approved by the institutional review board of Xiyuan hospital. All animal manipulations and care were performed 1.5–3.5 h after lights were switched on.

The rats were narcotized with pentobarbital sodium (3%, 0.2 ml/100 g) and sacrificed by abdominal aortic exsanguination at week 22. Serum samples were obtained from the blood by centrifugation (2,000 ×g, 4°C, 10 min) and stored at −80°C. Liver tissues were harvested, immediately frozen in liquid nitrogen, and stored at −80°C. Histopathological analyses were performed by routine staining of the frozen sections with hematoxylin and eosin.

Clinical Indices

FBG levels were measured using Accu-Chek Sensor Comfort test strips (Roche, Germany). Serum levels of triglyceride, cholesterol, alanine aminotransferase, and glutamic oxaloacetic transaminase were measured using enzymatic methods with an automatic biochemistry analyzer (Hitachi 7020, Tokyo, Japan).

RNA Extraction, Purification, Library Preparation, and Sequencing

Total RNA was extracted and purified from the frozen liver tissues of two ZDF and two ZCL rats using TRIzol (Life Technologies), according to manufacture instructions. The RNA integrity number was calculated to analyze RNA integrity using an Agilent Bioanalyzer 2100 (Agilent Technologies). Subsequently, poly(A)⁺ mRNA was isolated using the Dynabeads^® mRNA Purification Kit (Invitrogen), following manufacturer instructions. We then fragmented and purified RNA, synthesized double-stranded cDNA, and prepared indexed Illumina libraries, as described for RNase H libraries. Finally, the library was sequenced on an Illumina HiSeq 4000 system (Illumina, USA), as per the corresponding user guide.

Sequencing Data Processing

For each sample, the quality of raw sequence data was assessed using FastQC (https://www.bioinformatics.babraham.ac.uk/projects/fastqc/). Reads with mapping quality < 30 were discarded using SAMtools v1.0 (14). Only paired reads were used for subsequent analyses. Clean reads were then aligned to the UCSC rat reference genome (rn6) using TopHat v2.0.9, with default parameters (15, 16). Cufflinks v2.2.1 was used for assembling transcripts (16, 17). DEGs were identified using Cuffdiff v2.2.1 (16). The GTF annotation file for mRNAs was downloaded from the Ensemble website for transcript assembly and gene annotation and that for long noncoding RNAs was downloaded from the NONCODE v5 website.

Canonical Pathway and Biofunction Analyses

We analyzed differential biofunctions and canonical pathways by IPA (QIAGEN, USA) using differentially expressed mRNAs (FPKM ≥1, p < 0.05, |fold-change| ≥ 2), as well as differentially expressed lncRNA (FPKM ≥ 1, p < 0.05, |fold-change| ≥ 2) targeted mRNAs. Considering the observed relative abundance of each gene in the cluster and literature-derived regulation direction, the activation states of the biofunctions and canonical pathways were evaluated using the activation Z-score. Z-score > 0 indicated that the cluster was activated and Z-score < 0 indicated that the cluster was inhibited in ZDF rats.

Statistical Analysis

Statistical analyses were performed using R software v3.4.2. p < 0.05 indicated statistical significance.

Results

ZDF Rats Displayed Typical Features of T2DM

To verify the clinical indices of T2DM in male ZDF rats, we measured pertinent physiological and biochemical indices such as body weight, blood glucose levels, blood lipid levels, and serum hepatic marker enzymes (Figure 1). Healthy ZCL rats served as controls. Body weight and fasting blood glucose (FBG) levels were constantly monitored from weeks 6 to 22; blood lipid levels and serum hepatic marker enzymes were measured at weeks 7, 10, 12, 14, 18, and 22. In comparison with ZCL rats, ZDF rats showed significantly elevated FBG levels with time, which developed into persistent hyperglycemia from week 9 (Figure 1A). Moreover, ZDF rats gained weight from weeks 6 to 17, which reached a plateau at around week 18 (Figure 1B). At all measured timepoints, serum triglyceride and serum cholesterol levels in ZDF rats were significantly higher than those in ZCL rats (Figure 1C). The aforementioned physiological and biochemical indices demonstrated some typical features of T2DM in ZDF rats (6).

FIGURE 1

Figure 1 Clinical indices revealing the progression of T2DM in Zucker diabetic fatty (ZDF) rats. Physiological and biochemical indices, including body weight (A), fasting blood glucose (B), serum lipid metabolism indexes (C) and hepatic function indexes (D), were constantly monitored in ZDF and healthy ZCL rats from weeks 6 to 22. Five rats were present within each group, and values indicate mean ± SEM. TG, triglycerides; CHOL, total cholesterol; ALT, glutamic pyruvic transaminase; AST, glutamic oxaloacetic transaminase.

Figure 1D shows that in ZDF rats, the levels of two serum liver enzymes (ALT, glutamic pyruvic transaminase; AST, glutamic oxaloacetic transaminase) were increased, indicating T2DM-induced progressive liver injury. To further verify the presence of liver injury, we performed Hepatic histopathology of 22-week-old ZDF rats and ZCL rats (Figure 2). In comparison with ZCL rats, the liver tissue of ZDF rats showed extensive cell turgescence, cytoplasmic degeneration, and narrowed hepatic sinus (Figure 2), demonstrating diffuse hepatocyte edema (18). A few hepatic cells with microvesicular steatosis were also observed in the hepatic lobules (Figure 2), suggestive of disorders of lipid metabolism (19). These results confirmed the occurrence of liver injury in ZDF rats.

FIGURE 2

Figure 2 HE, Masson trichrome and Oil red staining of the liver samples from the 22-week-old Zucker diabetic fatty (ZDF) and Zucker control lean (ZCL) rats (Scale bar: 50 µm; 400x).

DEGs and Differentially Expressed lncRNAs in ZDF and ZCL Rats

Through transcriptomic analyses of liver tissues from 22-week-old ZDF and ZCL rats, we identified 214 DEGs (p < 0.05, |fold-change| ≥ 2, FPKM ≥ 1; Figure 3A). Among them, 28 DEGs were exclusively expressed in ZDF rats, some of which have been reported to be associated with diabetes, pancreatitis, and pancreatic calcification (Figure 3B) (20). On the other hand, three DEGs (Adm2, MGC108823, and Gimd1) were specifically expressed in ZCL rats (Figure 3B).

FIGURE 3

Figure 3 Profiles of differentially expressed mRNAs between Zucker control lean (ZCL) and Zucker diabetic fatty (ZDF) rats. (A) Heatmap showing the DEGs between ZCL and ZDF rats (p < 0.05, |fold-change| ≥ 2, FPKM ≥ 1). The color key indicates the relative expression level (log10 (FPKM)) between ZCL and ZDF rats, ranging from green (low expression) to red (high expression). (B) Heatmaps showing the representative DEGs between ZCL and ZDF rats. Specifically expressed DEGs and top 20 up- and down-regulated DEGs are separately exhibited (B). The colors indicate the expression level (log10 (FPKM)), ranging from yellow (low expression) to red (high expression). The value of are shown in the cells.

Except the above 31 specifically expressed genes, 183 DEGs (85.5%) were identified to express in both ZDF and ZCL samples. In comparison with ZCL rats, 129 and 54 DEGs were up- and downregulated in ZDF rats, respectively (Figure 3B). Among them, 89 upregulated and 24 downregulated DEGs (Supplementary Table S1) were found to be associated with diabetes progression (e.g., Ctrb1, Cpa2, and Insig1) and diabetes (e.g., Car3, Dio1, and Mup5), respectively (21, 22).

We then analyzed lncRNAs as they have been previously linked to diabetes (23). As shown in Figure 4, 104 differentially expressed lncRNAs (p < 0.05, |fold-change| ≥ 2, FPKM ≥ 1) were identified; of them, 23 and 17 lncRNAs were specifically expressed in ZDF and ZCL rats, respectively. The remaining 64 lncRNAs were common; 30 and 34 lncRNAs were up- and downregulated in ZDF rats, respectively (Figures 4A, B). In addition, according to the genomic position of lncRNAs, we classified them into five types: exonic, lincRNA, sense non-exonic, antisense and unclassified lncRNAs (Figure 4C) (24). The results indicated the differential distribution of lncRNA-types between the ZDF and ZCL groups. In comparison with ZCL rats, the proportion of four types of lncRNAs (exonic, linc, sense non-exonic, and antisense) was lower in ZDF rats.

FIGURE 4

Figure 4 Profiles of differentially expressed lncRNAs between Zucker control lean (ZCL) and Zucker diabetic fatty (ZDF) rats. (A) Heatmap showing the differentially expressed lncRNAs between ZCL and ZDF rats (p < 0.05, |fold-change| ≥ 2, FPKM ≥ 1). The color key indicates the relative expression level (log10 (FPKM)) between ZCL and ZDF rats, ranging from green (low expression) to red (high expression). (B) Heatmaps showing the representative differentially expressed lncRNAs between ZCL and ZDF rats. Specifically expressed differentially expressed lncRNAs and top 20 up- and down-regulated DEGs are separately exhibited (B). The colors indicate the expression level (log10 (FPKM)), ranging from yellow (low expression) to red (high expression). The value of are shown in the cells. (C) The proportions of lncRNA types between ZCL and ZDF rats. The lncRNAs were classified into five types based on the position in the genome, including exonic, lincRNA, sense non-exonic, antisense and unclassified lncRNAs. Each type of circRNA is represented by the indicated color box.

Pathway and Biofunction Analyses of DEGs and Differentially Expressed lncRNAs

To explore the molecular pathogenesis of T2DM, we used liver samples obtained from ZDF rats and employed the interpretative phenomenological analysis (IPA) approach to subject DEGs and differentially expressed lncRNAs to “canonical pathway” and “biofunction” analyses; all significantly enriched items were then screened (p < 0.05). We first conducted “canonical pathway” analyses of the 214 DEGs using the IPA approach, identifying 84 enriched canonical pathways (p < 0.05). Further analyses of the top 20 significantly enriched canonical pathways revealed the enrichment of hepatic function-associated substance metabolism (particularly degradation) pathways and inflammatory pathways (Figure 5A). Sixteen substance metabolism-associated pathways were identified, including 12 degradation-associated pathways (e.g., “serotonin degradation”, “dopamine degradation” and “ethanol degradation II”); four inflammatory pathways were detected (“LPS/IL-1-mediated inhibition of RXR function,” “PXR/RXR activation”, “aryl hydrocarbon receptor signaling” and “hepatic fibrosis/hepatic stellate cell activation”).

FIGURE 5

Figure 5 Most significantly affected canonical pathways of differentially expressed genes (DEGs) between Zucker control lean (ZCL) and Zucker diabetic fatty (ZDF) rats. (A) Top 20 most significantly affected IPA canonical pathways of differentially expressed mRNAs. (B) Top 20 most significantly affected IPA canonical pathways of mRNAs targeted by differentially expressed lncRNAs. (C) Top 20 most significantly affected IPA canonical pathways of differentially expressed mRNAs, as well as mRNAs targeted by differentially expressed lncRNAs. The y-axis indicates the significance (-log10(p-value)>1.3) of DEG-enriched pathways.

We then conducted “canonical pathway” analyses of the 104 differentially expressed lncRNA-targeted mRNAs using the IPA approach. Seventy-nine canonical pathways (p < 0.05) were enriched, including 41 hepatic function-associated substance metabolism (specifically degradation) pathways and 38 inflammatory pathways; these results concurred with the aforementioned findings. The top 20 significantly enriched pathway items (Figure 5B) comprised 12 substance metabolism-associated pathways (e.g., “nicotine degradation III”, “nicotine degradation II” and “serotonin degradation”) and eight inflammatory pathways (e.g., “PXR/RXR activation”, “LPS/IL-1-mediated inhibition of RXR function” and “LXR/RXR activation”). Among them, 11 pathways were shared with the top 20 pathway items enriched from DEGs (e.g., “LPS/IL-1-mediated inhibition of RXR function”, “PXR/RXR activation” and “serotonin degradation”; Figure 5C), indicating a synergistic effect between mRNAs and lncRNAs.

Further, we used the IPA approach for “biofunction” analyses of DEGs; we identified top 5 significantly enriched items belonging to three categories (“molecular and cellular functions”, “physiological system development and function” and “diseases and disorders”; Table 1). Moreover, we observed the enrichment of some hepatic function-associated metabolism and inflammatory biofunction items (p < 0.05), such as “carbohydrate metabolism,” “hepatic system development and function,” and “inflammatory response.” These findings were further supported by “biofunction” analyses of differentially expressed lncRNA-targeted mRNAs using the IPA approach (Supplementary Table S2). Importantly, seven biofunctions were common between DEGs and differentially expressed lncRNA-targeted mRNAs (e.g., “lipid metabolism”, “digestive system development and function” and “organismal injury and abnormalities”), indicating a synergistic effect between mRNAs and lncRNAs.

TABLE 1

Table 1 Top-5 IPA Bio Function items of each categories using DEGs.

Molecular Pathogenesis Analysis of T2DM in ZDF Rats

To explore the molecular pathogenesis of T2DM in the liver tissues of ZDF rats, we analyzed all significant differentially expressed canonical pathways and biofunction items using activation Z-score values as well as clinical indices. The comprehensive analysis indicated some typical features of T2DM in ZDF rats, such as abnormal metabolism, increased inflammation, liver injury and increased endoplasmic reticulum (ER) stress (Figure 6 and Supplementary Figure S1).

FIGURE 6

Figure 6 Assessment of the biofunction activity in Zucker diabetic fatty (ZDF) rats. Using the IPA software, we performed functional enrichment analyses of DEGs and evaluated the activation Z-score of the biofunction items, such as “lipid metabolism”- (A), “carbohydrate metabolism”- (B), “protein metabolism”- (C), “immune cell trafficking”- (D), “liver injury”- (E), and “ER stress”- (F) related functions. Center icons in networks, significantly activated/inhibited biofunctions; peripheral icons, DEGs between ZCL and ZDF rats and significantly activated/inhibited canonical pathways. Colors of DEGs indicate the estimate of expression levels, (log10 (FPKM)) ranging from green (low expression) to red (high expression). Colors of biofunctions and canonical pathways indicate the activation Z-score, ranging from blue (inhibited) to orange (activated). Orange lines indicate that DEGs support the activation of biofunctions/canonical pathways. Blue lines indicate that DEGs support the inhibition of biofunctions/canonical pathways. Gray lines indicate that DEGs are involved in biofunctions/canonical pathways, but do not affect their activity. DEGs, differentially expressed genes; CP, canonical pathways.

Abnormal Metabolism of Lipid, Glucose, and Protein

In 22-week-old ZDF rats, 94 DEGs were enriched for biofunction items related to lipid, glucose, and protein metabolisms. Most (29 out of 41) lipid metabolism-related biofunctions were activated (Figure 6A), such as “uptake of lipid” (Z-score = 1.94), “synthesis of lipid” (Z-score = 0.84), “conversion of fatty acid” (Z-score = 1.95), “hydrolysis of lipid” (Z-score = 1.13), and “export of lipid” (Z-score = 1.56). Correspondingly, “concentration of lipid” (Z-score = −0.82) and “concentration of cholesterol” (Z-score = −1.78) biofunctions were inhibited. These results supported the differences in body weight and hepatic histopathology result between ZDF and ZCL rats at week 22 (Figures 1B and 2).

As for glucose metabolism, “synthesis of D-glucose” (Z-score = 1.72), “gluconeogenesis” (Z-score = 1.62), and “hyperglycemia” (Z-score = 1.67) biofunctions were activated in ZDF rats, and these findings supported the elevation of FBG levels (Figures 1A and 6B).4 protein metabolism-related biofunctions were also activated in ZDF rats, such as “proteolysis” (Z-score = 1.63) and “metabolism of protein” (Z-score = 1.51) (Figure 6C). In addition, 78 and 17 lipid/glucose/protein metabolism-related genes were up- and downregulated in ZDF rats, respectively, including glucokinase (GCK), insulin-like growth factor-binding protein-1 (IGFBP1), phosphoenolpyruvate carboxykinase-1 (PCK1), pancreatic lipase (PNLIP), CPA1/2, CTRL, CTRB2, and CELA2A (Figures 6A–C, Supplementary Table S1).

Increased Inflammation

Canonical pathway and biofunction analyses, as well as clinical indices, indicated increased inflammation in ZDF rats. The physiological results revealed the mild infiltration of inflammatory cells in the hepatic lobules and portal area in 22-week-old ZDF rats, indicating increased inflammation (Figure 2).

Comparative transcriptome analyses also indicated increased inflammation in ZDF rats upon analyzing related DEGs, canonical pathways, and biofunction items. The expression levels of 92 inflammation response-related DEGs, including EPHX2, G6PC, GSN, SPP1, DMBT1, CDKN1A, ID2, SLC22A8, and SLC7A2, in ZDF rat liver tissues were significantly altered. In addition, 17 inflammation response-related biofunction items were activated (Z-score > 0) in ZDF rats, including the upregulation of some “immune cell trafficking” biofunction items (Figure 6D), such as “cell movement of phagocyte” (Z-score = 2.51), “leukocyte migration” (z-score = 2.83), “cellular infiltration by granulocytes” (Z-score = 2.13), and “chemotaxis of neutrophils” (Z-score = 2.41). Moreover, canonical pathway analyses revealed four significantly altered inflammation-related pathways, including “ILK signaling” (Z-score = 2.24), “p38 MAPK signaling” (Z-score = 2.00), “LPS1/IL1-mediated inhibition of RXR function” (Z-score = 0.38), and “PPARα/RXRα activation” (Z-score = −0.45) (Figure 6D).

Liver Injury

Figure 6E shows the upregulation of eight liver injury-related biofunction items, such as “damage of liver” (Z-score = 1.66), “cell death of liver cells” (Z-score = 0.84), and “apoptosis of hepatocytes” (Z-score = 0.82), in liver tissue samples obtained from 22-week-old ZDF rats. Ninety liver injury-related DEGs were also identified, such as ID2, PLK3, IRF1, TNFRSF1B, and XBP1 (Figure 1D). This is consistent with the abnormality in serum liver enzymes and pathological observations (Figures 1D and 2).

Increased ER Stress

ER stress reportedly can inhibit the insulin signaling pathway, promote insulin resistance, and further accelerate T2DM progression (25). We found that the biofunction “endoplasmic reticulum stress response of cells” (Z-score = 1.08) was upregulated in ZDF rats (Figure 6F). Correspondingly, seven (RNF186, MYH11, ALOX15, CTH, PNLIPRP2, ATP2A1, and EPHX2) and three (XBP1, SDF2L1, and CREB3L1) ER stress-related genes were up- and downregulated, respectively, in ZDF rats.

Discussion

In this study, we performed comparative transcriptome analyses using liver tissue samples, which led to the identification of 214 DEGs and 104 differentially expressed lncRNAs in ZDF and ZCL rats. Pathway and biofunction analyses of DEGs and differentially expressed lncRNAs using the IPA approach indicated a synergistic effect between mRNAs and lncRNAs. Through comprehensive analyses of bioinformatics data and indices, we found some typical features of T2DM in ZDF rats, such as upregulated metabolism (significantly increased lipid absorption/transport/utilization, gluconeogenesis, and protein hydrolysis), increased inflammation, liver injury and increased ER stress.

Among the 214 DEGs, 114 have been reported to be closely related to T2DM occurrence and development (26–30) (Supplementary Table S3). Of these 114 DEGs, 83 were found to be involved in T2DM-related substance metabolism; 35 genes were noted to be involved in glycometabolism, 51 in lipid metabolism, and 33 in proteometabolism (Supplementary Figure S2). The abnormal expression of GCK, IGFBP1, and PCK1, for example, can cause the occurrence of hyperglycemia (26–28); moreover, the abnormal expression of PNLIP in patients with T2DM has been reported to influence lipid uptake (31) and that of some protein hydrolase/metabolism-related genes (such as CPA1/2, CTRL, CTRB2, and CELA2A) has been reported to be associated with T2DM (21, 29). In addition, in this study, the “ILK signaling” pathway was significantly altered, which is evidently associated with T2DM-related substance metabolism (Supplementary Figure S1). As a scaffold, ILK protein can regulate blood glucose levels by recruiting integrin, insulin receptor, and various downstream kinases (32). Of the 114 known T2DM-related DEGs, 65 have been reported to be closely related to some inflammation response-related biofunction items (e.g., “chemotaxis of neutrophils,” “leukocyte migration,” and “cellular infiltration by granulocytes”), such as EPHX2, G6PC, GSN, and SPP1. In addition, four significantly changed pathways were found to be involved in inflammatory responses: upregulation of “ILK signaling” pathway can lead to insulin resistance in T2DM (33) and that of “p38 MAPK signaling” pathway can upregulate immunoinflammatory responses in diabetes by activating cytokines. Moreover, the downregulation of “PPARα/RXRα activation” pathway can enhance inflammatory responses (34–36), and the activation of “LPS1/IL1-mediated inhibition of RXR function” pathway can inhibit the “PPARα/RXRα activation” pathway (Supplementary Figure S1), further facilitating the induction of inflammatory responses.

Fifty-eight DEGs of the 114 known T2DM-related DEGs have been reported to be associated with liver injury. IRF1, NR0B2, and TNFRSF1B are related to a significant increase in liver tissue/cell injury (37–39), and the abnormal expression of XBP1 and ATP2A1 can activate some downstream genes (e.g., Bcl2, Bax, and Bim) and cause hepatocyte death (40, 41). Six DEGs of the 114 known T2DM-related DEGs have been reported to be related to ER stress, such as XBP1, MYH11, ATP2A1, EPHX2, and CTH (40–43). We also noted that one ER stress-related pathway was significantly upregulated (“p38 MAPK signaling” pathway): the upregulation of this pathway evidently promotes ER stress and causes insulin resistance (44).

Moving on, we also identified 100 putative T2DM-related DEGs (Supplementary Figure S2, Supplementary Table S4): 62 were found to be involved in substance metabolism, inflammatory response, liver injury, and ER stress. Moreover, 34 DEGs have been reported to be associated with substance metabolism, such as EEPD1, Abca2, Acnat2, Acot12, Adh4, and Adh6. EEPD1, for instance, could regulate cellular cholesterol levels by upregulating the expression of the cholesterol efflux transporters ABCA1 and ABCG1 (45). Supplementary Table S1 shows that EEPD1 was specifically expressed in the liver tissue of ZDF rats, suggesting its involvement in T2DM progression.

Of the 100 putative T2DM-associated DEGs, 27 were inflammation response-related DEGs, including DMBT1, CDKN1A (p21), ID2, SLC22A8, and SLC7A2. p21 is an important proinflammatory factor in lymphotoxin-driven pancreatic injury. It could regulate innate immune cell recruitment by influencing the secretion of inflammatory mediators (46).

In addition, we identified 32 liver injury-related genes among the 100 putative T2DM-associated DEGs, including PLK3, CRELD2, DMBT1, FCHSD2, and FNIP2. It is noteworthy that LRH-1 can reportedly induce PLK3 upregulation, leading to the activation of the transcription factor ATF2, consequently causing hepatocyte death (47).

We also realized that four putative T2DM-associated DEGs have been reported to be involved in ER stress (48–51): PNLIPRP2, SDF2L1, ALOX15, and RNF186. The RING finger E3 ligase RNF186 evidently regulates ER stress-mediated apoptosis through interaction with BNip1 in HeLa cells (36, 51–54). As one of the specifically expressed genes in ZDF rats, we deduce that RNF186 plays a key role in T2DM progression.

Interestingly, we found that two significantly enriched disease-associated pathways might be associated with T2DM (“neuroprotective role of THOP1 in Alzheimer’s disease” and “osteoarthritis pathway”) (Supplementary Figure S1). A previous proteomic analysis on diabetic retinas (55) revealed that diabetic retinal neurodegeneration shared some common pathogenic pathways with brain neurodegenerative conditions, such as Alzheimer’s and Parkinson’s diseases. Herein we also observed that the “neuroprotective role of THOP1 in Alzheimer’s disease” pathway was upregulated in the liver tissue of ZDF rats [Z-score = 2.45, −log(p) = 3.09], suggesting a neuroprotective role in T2DM progression.

Previous research indicated that the increase of inflammatory cytokines could cause insulin resistance by influencing the phosphorylation of IRS in insulin signal transduction (56, 57). This is in agreement with our findings, which showed the increased inflammation in canonical pathways, biofunction items and DEGs for the ZDF samples (Figure 6D). On the other hand, abnormal metabolism of glucose is one of the main phenotypes of insulin resistance (56, 57). Our study indicated the higher FBG level in the ZDF samples (Figure 1A). Biofunction analysis also showed the activation of “synthesis of D-glucose” (Z-score = 1.72), “gluconeogenesis” (Z-score = 1.62), and “hyperglycemia” (Z-score = 1.67) (Figure 6B). Here liver insulin resistance contributes to the excessive HGO, which is highly correlated with hyperglycemia and the increase of FBG level in type 2 diabetes (58).

In addition to insulin resistance, previous studies and our work indicated the phenotypes of obesity and fatty liver in ZDF rats (12) (Figures 1B and 2). Our findings also revealed 15 hepatic steatosis associated DEGs, including CD36, CIDEA, GPT, MFSD2A, MYC, NR0B2, NR1I3, PNLIP, RGS16, XBP1, SULTLE1, FABP5, GCK, PNPLA3 and SPP1 (59–68). Besides, biofunction analysis showed the abnormal metabolism of lipid (Figure 6A) and upregulation of liver injury (Figure 1D) in ZDF rats. Overall, both phenotype and transcriptome analyses indicated the fatty liver in ZDF rats.

There is a limitation in this study. Our research only study the liver transcriptomes of ZDF rats at one time-point (22-week-old), since they had the stable phenotype of T2DM at that time-point. Certainly, this warrants further study to use more tissues/organs and sampling time-points in the future.

In conclusion, we identified 214 genes that were differentially expressed between ZDF and ZCL rats: 114 known and 100 putative T2DM-related genes. The 114 known T2DM-related genes included 69 substance (glucose/lipid/protein) metabolism-, 65 inflammation response-, 58 liver injury-, 15 hepatic steatosis- and six ER stress-related genes. Further, the 100 putative T2DM-related genes included 25 substance (glucose/lipid/protein) metabolism-, 27 inflammation response-, 32 liver injury-, and four ER stress-related genes. In addition, we observed 84 significantly enriched pathways associated with T2DM progression, including “ILK signaling,” “p38 MAPK signaling,” and “PPARα/RXRα activation” pathways. Our results should serve as an important reference and improve the understanding of obesity-associated T2DM; moreover, the data reported herein should help in the identification of some potential diagnostic markers and therapeutic targets.

Data Availability Statement

The datasets presented in this study can be found in online repositories. The names of the repository/repositories and accession number(s) can be found below: https://www.ncbi.nlm.nih.gov/geo/, GSE117447.

Ethics Statement

The animal study was reviewed and approved by Xiyuan hospital, China Academy of Chinese Medical Sciences.

Author Contributions

FC and GM designed the study. CX, GM, JW, YS, and KN conducted the experiments. XZ, TC, CL, and LY performed the analyses and interpreted the results. FC and XZ wrote the manuscript. All authors contributed to the article and approved the submitted version.

Funding

This work was supported by the National Natural Science Foundation of China (81373581). The funding body had no role in the design of the study and collection, analysis, and interpretation of data and in writing the manuscript.

Conflict of Interest

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Supplementary Material

The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fendo.2020.565858/full#supplementary-material

Abbreviations

ALT, alanine aminotransferase; AST, glutamic-oxaloacetic transaminase; CHOL, total cholesterol; DEGs, differentially expressed genes; DM, diabetes mellitus; ER, endoplasmic reticulum; HGO, hepatic glucose output; ILK, integrin linked kinase; IRS, insulin receptor substrate; IPA, Ingenuity Pathway Analysis; lncRNA, long noncoding RNA; LPS, lipopolysaccharide; MAPK, mitogen activated protein kinases; TG, triglycerides; T2DM, type 2 diabetes mellitus; ZDF, Zucker diabetic fatty.

References

1. Shoelson SE, Herrero L, Naaz A. Obesity, inflammation, and insulin resistance. Gastroenterology (2007) 132:2169–80. doi: 10.1053/j.gastro.2007.03.059

PubMed Abstract | CrossRef Full Text | Google Scholar

2. Saeedi P, Petersohn I, Salpea P, Malanda B, Karuranga S, Unwin N, et al. Global and regional diabetes prevalence estimates for 2019 and projections for 2030 and 2045: Results from the International Diabetes Federation Diabetes Atlas, 9(th) edition. Diabetes Res Clin Pract (2019) 157:107843. doi: 10.1016/j.diabres.2019.107843

PubMed Abstract | CrossRef Full Text | Google Scholar

3. Zimmet P, Alberti KG, Shaw J. Global and societal implications of the diabetes epidemic. Nature (2001) 414:782–7. doi: 10.1038/414782a

PubMed Abstract | CrossRef Full Text | Google Scholar

4. Rathmann W, Giani G. Global prevalence of diabetes: estimates for the year 2000 and projections for 2030. Diabetes Care (2004) 27:2568–9; author reply 2569. doi: 10.2337/diacare.27.10.2568

PubMed Abstract | CrossRef Full Text | Google Scholar

5. Ashraf NU, Sheikh TA. Endoplasmic reticulum stress and Oxidative stress in the pathogenesis of Non-alcoholic fatty liver disease. Free Radical Res (2015) 49:1405–18. doi: 10.3109/10715762.2015.1078461

CrossRef Full Text | Google Scholar

6. Fujii J, Homma T, Kobayashi S, Seo HG. Mutual interaction between oxidative stress and endoplasmic reticulum stress in the pathogenesis of diseases specifically focusing on non-alcoholic fatty liver disease. World J Biol Chem (2018) 9:1–15. doi: 10.4331/wjbc.v9.i1.1

PubMed Abstract | CrossRef Full Text | Google Scholar

7. Targher G, Byrne CD. Clinical Review: Nonalcoholic fatty liver disease: a novel cardiometabolic risk factor for type 2 diabetes and its complications. J Clin Endocrinol Metab (2013) 98:483–95. doi: 10.1210/jc.2012-3093

PubMed Abstract | CrossRef Full Text | Google Scholar

8. Shiota M, Printz RL. Diabetes in Zucker diabetic fatty rat. Methods Mol Biol (2012) 933:103–23. doi: 10.1007/978-1-62703-068-7_8

PubMed Abstract | CrossRef Full Text | Google Scholar

9. Yang X, Pratley RE, Tokraks S, Bogardus C, Permana PA. Microarray profiling of skeletal muscle tissues from equally obese, non-diabetic insulin-sensitive and insulin-resistant Pima Indians. Diabetologia (2002) 45:1584–93. doi: 10.1007/s00125-002-0905-7

PubMed Abstract | CrossRef Full Text | Google Scholar

10. Sreekumar R, Halvatsiotis P, Schimke JC, Nair KS. Gene expression profile in skeletal muscle of type 2 diabetes and the effect of insulin treatment. Diabetes (2002) 51:1913–20. doi: 10.2337/diabetes.51.6.1913

PubMed Abstract | CrossRef Full Text | Google Scholar

11. Suh YH, Kim Y, Bang JH, Choi KS, Lee JW, Kim WH, et al. Analysis of gene expression profiles in insulin-sensitive tissues from pre-diabetic and diabetic Zucker diabetic fatty rats. J Mol Endocrinol (2005) 34:299–315. doi: 10.1677/jme.1.01679

PubMed Abstract | CrossRef Full Text | Google Scholar

12. Nawano M, Oku A, Ueta K, Umebayashi T, Ishirahara T, Arakawa A, et al. Hyperglycemia contributes insulin resistance in hepatic and adipose tissue but not skeletal muscle of ZDF rats. Am J Physiol Endocrinol Metab (2000) 278:E535–43. doi: 10.1152/ajpendo.2000.278.3.E535

PubMed Abstract | CrossRef Full Text | Google Scholar

13. Slieker LJ, Sundell KL, Heath WF, Osborne HE, Bue J, Manetta J, et al. Glucose transporter levels in tissues of spontaneously diabetic Zucker fa/fa rat (ZDF/drt) and viable yellow mouse (Avy/a). Diabetes (1992) 41:187–93. doi: 10.2337/diab.41.2.187

PubMed Abstract | CrossRef Full Text | Google Scholar

14. Li H, Handsaker B, Wysoker A, Fennell T, Ruan J, Homer N, et al. The Sequence Alignment/Map format and SAMtools. Bioinformatics (2009) 25:2078–9. doi: 10.1093/bioinformatics/btp352

PubMed Abstract | CrossRef Full Text | Google Scholar

15. Trapnell C, Pachter L, Salzberg SL. TopHat: discovering splice junctions with RNA-Seq. Bioinformatics (2009) 25:1105–11. doi: 10.1093/bioinformatics/btp120

PubMed Abstract | CrossRef Full Text | Google Scholar

16. Kim D, Pertea G, Trapnell C, Pimentel H, Kelley R, Salzberg SL, et al. TopHat2: accurate alignment of transcriptomes in the presence of insertions, deletions and gene fusions. Genome Biol (2013) 14:R36. doi: 10.1186/gb-2013-14-4-r36

PubMed Abstract | CrossRef Full Text | Google Scholar

17. Trapnell C, Williams BA, Pertea G, Mortazavi A, Kwan G, van Baren MJ, et al. Transcript assembly and quantification by RNA-Seq reveals unannotated transcripts and isoform switching during cell differentiation. Nat Biotechnol (2010) 28:511–5. doi: 10.1038/nbt.1621

PubMed Abstract | CrossRef Full Text | Google Scholar

18. Ellis JA, Allan G, Krakowka S. Effect of coinfection with genogroup 1 porcine torque teno virus on porcine circovirus type 2-associated postweaning multisystemic wasting syndrome in gnotobiotic pigs. Am J Vet Res (2008) 69:1608–14. doi: 10.2460/ajvr.69.12.1608

PubMed Abstract | CrossRef Full Text | Google Scholar

19. Tandra S, Matthew MY, Elizabeth MB, Vuppalanchi R, Oscar WC, Aynur ÜA, et al. Presence and significance of microvesicular steatosis in nonalcoholic fatty liver disease. J Hepatol (2011) 55:654–9. doi: 10.1016/j.jhep.2010.11.021

PubMed Abstract | CrossRef Full Text | Google Scholar

20. Madro A, Ciesielka M, Celinski K, Slomka M, Czechowska G, Kurzepa J, et al. The genetic predisposition and its impact on the diabetes mellitus development in patients with alcoholic chronic pancreatitis. Gastroenterol Res Pract (2015) 2015:309156. doi: 10.1155/2015/309156

PubMed Abstract | CrossRef Full Text | Google Scholar

21. M 't Hart L, Fritsche A, Nijpels G, van Leeuwen N, Donnelly LA, Dekker JM, et al. The CTRB1/2 locus affects diabetes susceptibility and treatment via the incretin pathway. Diabetes (2013) 62:3275–81. doi: 10.2337/db13-0227

PubMed Abstract | CrossRef Full Text | Google Scholar

22. Sun SH, Liu SQ, Cai CP, Cai R, Chen L, Zhang QB, et al. Down-regulation of alpha-2u globulin in renal mitochondria of STZ-induced diabetic rats observed by a proteomic method. Annales d’endocrinol (2012) 73:530–41. doi: 10.1016/j.ando.2012.09.007

CrossRef Full Text | Google Scholar

23. Feng SD, Yang JH, Yao CH, Yang SS, Zhu ZM, Wu D, et al. Potential regulatory mechanisms of lncRNA in diabetes and its complications. Biochem Cell Biol = Biochim Biol Cell (2017) 95:361–7. doi: 10.1139/bcb-2016-0110

CrossRef Full Text | Google Scholar

24. Kopp F, Mendell JT. Functional Classification and Experimental Dissection of Long Noncoding RNAs. Cell (2018) 172:393–407. doi: 10.1016/j.cell.2018.01.011

PubMed Abstract | CrossRef Full Text | Google Scholar

25. Kim OK, Jun W, Lee J. Mechanism of ER Stress and Inflammation for Hepatic Insulin Resistance in Obesity. Ann Nutr Metab (2015) 67:218–27. doi: 10.1159/000440905

PubMed Abstract | CrossRef Full Text | Google Scholar

26. Park JM, Kim TH, Jo SH, Kim MY, Ahn YH. Acetylation of glucokinase regulatory protein decreases glucose metabolism by suppressing glucokinase activity. Sci Rep (2015) 5:17395. doi: 10.1038/srep17395

PubMed Abstract | CrossRef Full Text | Google Scholar

27. van Dijk PR, Logtenberg SJJ, Groenier KH, Kleefstra N, Bilo HJG, Arnqvist HJ, et al. Effect of i.p. insulin administration on IGF1 and IGFBP1 in type 1 diabetes. Endocr Connect (2014) 3:17–23. doi: 10.1530/EC-13-0089

PubMed Abstract | CrossRef Full Text | Google Scholar

28. Beale EG, Hammer RE, Antoine B, Forest C. Disregulated glyceroneogenesis: PCK1 as a candidate diabetes and obesity gene. Trends Endocrinol Metabol: TEM (2004) 15:129–35. doi: 10.1016/j.tem.2004.02.006

CrossRef Full Text | Google Scholar

29. Han D, Moon S, Kim H, Choi S-E, Lee S-J, Park KS, et al. Detection of differential proteomes associated with the development of type 2 diabetes in the Zucker rat model using the iTRAQ technique. J Proteome Res (2011) 10:564–77. doi: 10.1021/pr100759a

PubMed Abstract | CrossRef Full Text | Google Scholar

30. Dziewulska A, Dobosz AM, Dobrzyn A. High-Throughput Approaches onto Uncover (Epi)Genomic Architecture of Type 2 Diabetes. Genes (2018) 9(8):374. doi: 10.3390/genes9080374

CrossRef Full Text | Google Scholar

31. Lindner I, Helwig U, Rubin D, Li Y, Fisher E, Boeing H, et al. Putative association between a new polymorphism in exon 3 (Arg109Cys) of the pancreatic colipase gene and type 2 diabetes mellitus in two independent Caucasian study populations. Mol Nutr Food Res (2005) 49:972–6. doi: 10.1002/mnfr.200500087

PubMed Abstract | CrossRef Full Text | Google Scholar

32. Dongiovanni P, Meroni M, Baselli GA, Bassani GA, Rametta R, Pietrelli A, et al. Insulin resistance promotes Lysyl Oxidase Like 2 induction and fibrosis accumulation in non-alcoholic fatty liver disease. Clin Sci (2017) 131:1301–15. doi: 10.1042/CS20170175

CrossRef Full Text | Google Scholar

33. Williams AS, Trefts E, Lantier L, Grueter CA, Bracy DP, James FD, et al. Integrin-Linked Kinase Is Necessary for the Development of Diet-Induced Hepatic Insulin Resistance. Diabetes (2017) 66:325–34. doi: 10.2337/db16-0484

PubMed Abstract | CrossRef Full Text | Google Scholar

34. Cave MC, López-Velázquez JA, Carrillo-Córdova LD, Chávez-Tapia NC, Uribe M, Méndez-Sánchez N, et al. Nuclear receptors and nonalcoholic fatty liver disease. Biochim Biophys Acta (2016) 1859:1083–99. doi: 10.1016/j.bbagrm.2016.03.002

PubMed Abstract | CrossRef Full Text | Google Scholar

35. Jenniskens M, Langouche L, Vanwijngaerden YM, Mesotten D, Van den Berghe G. Cholestatic liver (dys)function during sepsis and other critical illnesses. Intensive Care Med (2016) 42:16–27. doi: 10.1007/s00134-015-4054-0

PubMed Abstract | CrossRef Full Text | Google Scholar

36. Kallwitz ER, McLachlan A, Cotler SJ. Role of peroxisome proliferators-activated receptors in the pathogenesis and treatment of nonalcoholic fatty liver disease. World J Gastroenterol (2008) 14:22–8.

PubMed Abstract | Google Scholar

37. Friesen M, Camahort R, Lee Y-K, Xia F, Gerszten RE, Rhee EP, et al. Activation of IRF1 in Human Adipocytes Leads to Phenotypes Associated with Metabolic Disease. Stem Cell Rep (2017) 8:1164–73. doi: 10.1016/j.stemcr.2017.03.014

CrossRef Full Text | Google Scholar

38. Mitchell SMS, Weedon MN, Owen KR, Shields B, Wilkins-Wall B, Walker M, et al. Genetic variation in the small heterodimer partner gene and young-onset type 2 diabetes, obesity, and birth weight in U.K. subjects. Diabetes (2003) 52:1276–9. doi: 10.2337/diabetes.52.5.1276

PubMed Abstract | CrossRef Full Text | Google Scholar

39. Benjafield AV, Glenn CL, Wang XL, Colagiuri S, Morris BJ. TNFRSF1B in genetic predisposition to clinical neuropathy and effect on HDL cholesterol and glycosylated hemoglobin in type 2 diabetes. Diabetes Care (2001) 24:753–7. doi: 10.2337/diacare.24.4.753

PubMed Abstract | CrossRef Full Text | Google Scholar

40. Yang J, Wu X, Wu X, Zhou D, Lin T, Ding S, et al. The Multiple Roles of XBP1 in Regulation of Glucose and Lipid Metabolism. Curr Protein Pept Sci (2017) 18:630–5. doi: 10.2174/1389203717666160627085011

PubMed Abstract | CrossRef Full Text | Google Scholar

41. Takeuchi F, Yamamoto K, Katsuya T, Nabika T, Sugiyama T, Fujioka A, et al. Association of genetic variants for susceptibility to obesity with type 2 diabetes in Japanese individuals. Diabetologia (2011) 54:1350–9. doi: 10.1007/s00125-011-2086-8

PubMed Abstract | CrossRef Full Text | Google Scholar

42. Gong G, Yuan H, Liu Y, Qi L. Investigation of the Effects and Mechanisms of Mai Tong Formula on Lower Limb Macroangiopathy in a Spontaneous Diabetic Rat Model. J Diabetes Res (2016) 2016:8076796. doi: 10.1155/2016/8076796

PubMed Abstract | CrossRef Full Text | Google Scholar

43. Ma L, Jiang Y, Kong X, Yan M, Zhao T, Zhao H, et al. Synergistic Effect of the MTHFR C677T and EPHX2 G860A Polymorphism on the Increased Risk of Ischemic Stroke in Chinese Type 2 Diabetic Patients. J Diabetes Res (2017) 2017:6216205. doi: 10.1155/2017/6216205

PubMed Abstract | CrossRef Full Text | Google Scholar

44. Jeong K, Oh Y, Kim S-J, Kim H, Park K-C, Kim SS, et al. Apelin is transcriptionally regulated by ER stress-induced ATF4 expression viaa p38 MAPK-dependent pathway. Apoptosis (2014) 19:1399–410. doi: 10.1007/s10495-014-1013-0

PubMed Abstract | CrossRef Full Text | Google Scholar

45. Henriksen EJ, Diamond-Stanic MK, Marchionne EM. Oxidative stress and the etiology of insulin resistance and type 2 diabetes. Free Radical Biol Med (2011) 51:993–9. doi: 10.1016/j.freeradbiomed.2010.12.005

CrossRef Full Text | Google Scholar

46. Seleznik GM, Reding T, Peter L, Gupta A, Steiner SG, Sonda S, et al. Development of autoimmune pancreatitis is independent of CDKN1A/p21-mediated pancreatic inflammation. Gut (2018) 67:1663–73. doi: 10.1136/gutjnl-2016-313458

PubMed Abstract | CrossRef Full Text | Google Scholar

47. Wang L, Gonzalez S, Dai W, Deng S, Lu L. Effect of Hypoxia-regulated Polo-like Kinase 3 (Plk3) on Human Limbal Stem Cell Differentiation. J Biol Chem (2016) 291:16519–29. doi: 10.1074/jbc.M116.725747

PubMed Abstract | CrossRef Full Text | Google Scholar

48. Németh BC, Pesei ZG, Hegyi E, Szücs Á, Szentesi A, Hegyi P, et al. The common truncation variant in pancreatic lipase related protein 2 (PNLIPRP2) is expressed poorly and does not alter risk for chronic pancreatitis. PLoS One (2018) 13:e0206869. doi: 10.1371/journal.pone.0206869

PubMed Abstract | CrossRef Full Text | Google Scholar

49. Fukuda S, Sumii M, Masuda Y, Takahashi M, Koike N, Teishima J, et al. Murine and human SDF2L1 is an endoplasmic reticulum stress-inducible gene and encodes a new member of the Pmt/rt protein family. Biochem Biophys Res Commun (2001) 280:407–14. doi: 10.1006/bbrc.2000.4111

PubMed Abstract | CrossRef Full Text | Google Scholar

50. Zhang W, Zhong W, Sun Q, Sun X, Zhou Z. Hepatic overproduction of 13-HODE due to ALOX15 upregulation contributes to alcohol-induced liver injury in mice. Sci Rep (2017) 7:8976. doi: 10.1038/s41598-017-02759-0

PubMed Abstract | CrossRef Full Text | Google Scholar

51. Wang P, Wu Y, Li Y, Zheng J, Tang J. A novel RING finger E3 ligase RNF186 regulate ER stress-mediated apoptosis through interaction with BNip1. Cell Signal (2013) 25:2320–33. doi: 10.1016/j.cellsig.2013.07.016

PubMed Abstract | CrossRef Full Text | Google Scholar

52. Fujimoto K, Kinoshita M, Tanaka H, Okuzaki D, Shimada Y, Kayama H, et al. Regulation of intestinal homeostasis by the ulcerative colitis-associated gene RNF186. Mucosal Immunol (2017) 10:446–59. doi: 10.1038/mi.2016.58

PubMed Abstract | CrossRef Full Text | Google Scholar

53. Beaudoin M, Goyette P, Boucher G, Lo KS, Rivas MA, Stevens C, et al. Deep resequencing of GWAS loci identifies rare variants in CARD9, IL23R and RNF186 that are associated with ulcerative colitis. PLoS Genet (2013) 9:e1003723. doi: 10.1371/journal.pgen.1003723

PubMed Abstract | CrossRef Full Text | Google Scholar

54. Yang S-K, Hong M, Zhao W, Jung Y, Tayebi N, Ye BD, et al. Genome-wide association study of ulcerative colitis in Koreans suggests extensive overlapping of genetic susceptibility with Caucasians. Inflamm Bowel Dis (2013) 19:954–66. doi: 10.1097/MIB.0b013e3182802ab6

PubMed Abstract | CrossRef Full Text | Google Scholar

55. Sundstrom JM, Hernández C, Weber SR, Zhao Y, Dunklebarger M, Tiberti N, et al. Proteomic Analysis of Early Diabetic Retinopathy Reveals Mediators of Neurodegenerative Brain Diseases. Invest Ophthalmol Visual Sci (2018) 59:2264–74. doi: 10.1167/iovs.17-23678

CrossRef Full Text | Google Scholar

56. Hirosumi J, Tuncman G, Chang L, Görgün CZ, Uysal KT, Maeda K, et al. A central role for JNK in obesity and insulin resistance. Nature (2002) 420:333–6. doi: 10.1038/nature01137

PubMed Abstract | CrossRef Full Text | Google Scholar

57. Wellen KE, Hotamisligil GS. Inflammation, stress, and diabetes. J Clin Invest (2005) 115:1111–9. doi: 10.1172/JCI25102

PubMed Abstract | CrossRef Full Text | Google Scholar

58. Leonard BL, Watson RN, Loomes KM, Phillips AR, Cooper GJ. Insulin resistance in the Zucker diabetic fatty rat: a metabolic characterisation of obese and lean phenotypes. Acta Diabetol (2005) 42:162–70. doi: 10.1007/s00592-005-0197-8

PubMed Abstract | CrossRef Full Text | Google Scholar

59. Pashkov V, Huang J, Parameswara VK, Kedzierski W, Kurrasch DM, Tall GG, et al. Regulator of G protein signaling (RGS16) inhibits hepatic fatty acid oxidation in a carbohydrate response element-binding protein (ChREBP)-dependent manner. J Biol Chem (2011) 286:15116–25. doi: 10.1074/jbc.M110.216234

PubMed Abstract | CrossRef Full Text | Google Scholar

60. Abe T, Takahashi M, Kano M, Amaike Y, Ishii C, Maeda K, et al. Activation of nuclear receptor CAR by an environmental pollutant perfluorooctanoic acid. Arch Toxicol (2017) 91:2365–74. doi: 10.1007/s00204-016-1888-3

PubMed Abstract | CrossRef Full Text | Google Scholar

61. Berger JH, Charron MJ, Silver DL. Major facilitator superfamily domain-containing protein 2a (MFSD2A) has roles in body growth, motor function, and lipid metabolism. PLoS One (2012) 7:e50629. doi: 10.1371/journal.pone.0050629

PubMed Abstract | CrossRef Full Text | Google Scholar

62. Garbacz WG, Jiang M, Xu M, Yamauchi J, Dong HH, Xie W. Sex- and Tissue-Specific Role of Estrogen Sulfotransferase in Energy Homeostasis and Insulin Sensitivity. Endocrinology (2017) 158:4093–104. doi: 10.1210/en.2017-00571

PubMed Abstract | CrossRef Full Text | Google Scholar

63. Wang A-J, Yang Z, Grinchuk Y, Smith A, Qin B, Lu N, et al. IL-25 or IL-17E Protects against High-Fat Diet-Induced Hepatic Steatosis in Mice Dependent upon IL-13 Activation of STAT6. J Immunol (2015) 195:4771–80. doi: 10.4049/jimmunol.1500337

PubMed Abstract | CrossRef Full Text | Google Scholar

64. McGettigan BM, McMahan RH, Luo Y, Wang XX, Orlicky DJ, Porsche C, et al. Sevelamer Improves Steatohepatitis, Inhibits Liver and Intestinal Farnesoid X Receptor (FXR), and Reverses Innate Immune Dysregulation in a Mouse Model of Non-alcoholic Fatty Liver Disease. J Biol Chem (2016) 291:23058–67. doi: 10.1074/jbc.M116.731042

PubMed Abstract | CrossRef Full Text | Google Scholar

65. Jolly RA, Ciurlionis R, Morfitt D, Helgren M, Patterson R, Ulrich RG, et al. Microvesicular steatosis induced by a short chain fatty acid: effects on mitochondrial function and correlation with gene expression. Toxicol Pathol (2004) 32(Suppl 2):19–25. doi: 10.1080/01926230490451699

PubMed Abstract | CrossRef Full Text | Google Scholar

66. Garbacz WG, Lu P, Miller TM, Poloyac SM, Eyre NS, Mayrhofer G, et al. Hepatic Overexpression of CD36 Improves Glycogen Homeostasis and Attenuates High-Fat Diet-Induced Hepatic Steatosis and Insulin Resistance. Mol Cell Biol (2016) 36:2715–27. doi: 10.1128/MCB.00138-16

PubMed Abstract | CrossRef Full Text | Google Scholar

67. Caballería L, Pera G, Auladell MA, Torán P, Muñoz L, Miranda D, et al. Prevalence and factors associated with the presence of nonalcoholic fatty liver disease in an adult population in Spain. Eur J Gastroenterol Hepatol (2010) 22:24–32. doi: 10.1097/MEG.0b013e32832fcdf0

PubMed Abstract | CrossRef Full Text | Google Scholar

68. Lake AD, Novak P, Hardwick RN, Flores-Keown B, Zhao F, Klimecki WT, et al. The adaptive endoplasmic reticulum stress response to lipotoxicity in progressive human nonalcoholic fatty liver disease. Toxicol Sci (2014) 137:26–35. doi: 10.1093/toxsci/kft230

PubMed Abstract | CrossRef Full Text | Google Scholar