Plasma Oxylipin Profile Discriminates Ethnicities in Subjects with Non-Alcoholic Steatohepatitis: An Exploratory Analysis
Department of Nutrition, University of California Davis, Davis, CA 95616, USA; tmazi@ucdavis.edu (T.A.M.); pjhavel@ucdavis.edu (P.J.H.); john.newman2@usda.gov (J.W.N.)
Division of Clinical Nutrition, Department of Community Health Sciences, College of Applied Medical Sciences, King Saud University, Riyadh 11433, Saudi Arabia
West Coast Metabolomic Center, Genome Center, University of California Davis, Davis, CA 95616, USA; kborkowski@ucdavis.edu (K.B.); ofiehn@ucdavis.edu (O.F.)
Division of Gastroenterology and Hepatology, Department of Internal Medicine, University of California Davis, Sacramento, CA 95817, USA; clbowlus@ucdavis.edu (C.L.B.); ssarkar@ucdavis.edu (S.S.); dorothy.a.kieffer@gmail.com (D.A.K.)
Department of Pathology and Laboratory Medicine, University of California Davis, Sacramento, CA 95817, USA; kmatsukuma@ucdavis.edu (K.M.); yjywan@ucdavis.edu (Y.-J.Y.W.)
Department of Surgery, University of California Davis, Sacramento, CA 95817, USA; mrali@ucdavis.edu
Department of Molecular Biosciences, School of Veterinary Medicine, University of California Davis, Davis, CA 95616, USA; klstanhope@ucdavis.edu
United States Department of Agriculture-Agriculture Research Service-Western Human Nutrition Research Center, Davis, CA 95616, USA
Abstract
Non-alcoholic fatty liver disease (NAFLD) is a common liver pathology that includes steatosis, or non-alcoholic fatty liver (NAFL), and non-alcoholic steatohepatitis (NASH). Without a clear pathophysiological mechanism, it affects Hispanics disproportionately compared to other ethnicities. Polyunsaturated fatty acids (PUFAs) and inflammatory lipid mediators including oxylipin (OXL) and endocannabinoid (eCB) are altered in NAFLD and thought to contribute to its pathogenesis. However, the existence of ethnicity-related differences is not clear. We employed targeted lipidomic profiling for plasma PUFAs, non-esterified OXLs and eCBs in White Hispanics (HIS, n = 10) and Caucasians (CAU, n = 8) with biopsy-confirmed NAFL, compared with healthy control subjects (HC; n = 14 HIS; n = 8 CAU). NAFLD was associated with diminished long chain PUFA in HIS, independent of histological severity. Differences in plasma OXLs and eCBs characterized ethnicities in NASH, with lower arachidonic acid derived OXLs observed in HIS. The secondary analysis comparing ethnicities within NASH (n = 12 HIS; n = 17 CAU), confirms these ethnicity-related differences and suggests lower lipoxygenase(s) and higher soluble epoxide hydrolase(s) activities in HIS compared to CAU. While causes are not clear, these lipidomic differences might be with implications for NAFLD severity and are worth further investigation. We provide preliminary data indicating ethnicity-specific lipidomic signature characterizes NASH which requires further validation.
Untitled section
Keywords: arachidonic acid, Caucasian, cyclooxygenases, endocannabinoids, ethnicity, Hispanic, lipoxygenases, non-alcoholic steatohepatitis, oxylipins, polyunsaturated fatty acids, soluble epoxide hydrolases
Article notes
Untitled section
Received 2022 Feb 5; Accepted 2022 Feb 17; Collection date 2022 Feb.
1. Introduction
Non-alcoholic fatty liver disease (NAFLD) is a chronic liver condition affecting one in four adults worldwide and this rate increases with coexisting components of metabolic syndrome [1]. Its histological presentation includes hepatocellular steatosis, or non-alcoholic fatty liver (NAFL) with a range of necroinflammation with or without fibrosis. When hepatocellular damage and ballooning are present, this is clinically defined as non-alcoholic steatohepatitis (NASH) [2]. The pathogenesis of NAFLD is not fully elucidated. Whereas its onset involves an interplay between genetics and environmental factors with coexisting comorbidities, the progression to NASH appears to be provoked by multiple or parallel hits including oxidative stress and inflammation [3,4]. Oxidative stress modulates insulin signaling, lipid metabolism, inflammation, and fibrogenesis, and many oxidative stress biomarkers have been associated with NAFLD severity [5,6]. In the U.S., the risk and severity of NAFLD vary among ethnic/racial groups, with Hispanics (HIS) being affected disproportionately and presenting more frequently with advanced inflammation and fibrosis compared to other ethnicities [7,8,9]. The metabolic drivers underlying this disparity are not clear.
Polyunsaturated fatty acid (PUFAs) are bioactive lipids and precursors to inflammatory lipid mediators including oxylipins (OXLs) and endocannabinoids (eCBs). OXLs are produced from PUFAs by mono- and dioxygenases, including lipoxygenases (e.g., 5-LOX, 12-LOXs, and 15-LOXs); cyclooxygenases (i.e., COX-1 and -2), and a variety of cytochrome P450s (CYPs) [10]. PUFAs can also undergo non-enzymatic oxygenation mediated by free radicals and the rate of this production is increased under oxidative stress [11]. In general, OXLs from n-3 PUFAs have anti-inflammatory or less pro-inflammatory effect compared to those derived from n-6 PUFAs [12,13]. The fatty acid ethanolamides (i.e., N-acylethanolamides), one class of eCB, are synthesized by complex interactions of lipases and fatty acid amide hydrolase from PUFAs and membrane associated precursors [14]. Collectively, these lipids work through receptor-mediated mechanisms and likely contribute to NAFLD by modulating processes including lipogenesis, inflammation, and mitochondrial β-oxidation [10,15].
Previous lipidomic analyses showed that NAFLD is associated with dysregulated PUFAs metabolism [16,17,18,19,20]. Alterations in circulating OXL and eCB profiles are reported in NAFLD and other liver pathologies. In fact, numerous lipid mediators have been shown to predict NAFL or NASH [18,21,22,23,24,25,26]. However, metabolomic profiling in NAFLD with regards to ethnicity is limited. Our prior semi-quantitative lipidomic profiling study indicated ethnicity-specific differences in plasma PUFA profiles in subjects with NAFL, with higher abundance of linoleic acid (LA) and α-linolenic acid (ALA) seen in Caucasians (CAU) compared to ethnicity-matched lean subjects [27]. In the same study, the progression to NASH was characterized by ethnicity-specific differences in hepatic lipidomic profiles with higher levels of saturated and unsaturated fatty acids seen in NASH-HIS. Ethnicity was not previously addressed in OXL and eCB profiling efforts. Examining such lipidomic differences among ethnicities may shed light on potential mechanisms modulating the disparity in NAFLD prevalence and severity.
The objective of this “proof-of-concept” study is to examine ethnicity-related changes in PUFAs and their downstream inflammatory mediators in a group of subjects with obesity and biopsy-confirmed NAFL and NASH. We employed targeted lipidomic analysis of plasma PUFAs, OXLs, and the N-acylethanolamides class of eCBs to compare HIS and CAU subjects with medically complicated obesity to ethnicity-matched lean healthy controls (HC). Profiles in subjects diagnosed with NASH were also compared to ethnicity and BMI-matched participants without NASH (0-NASH). In addition, we conducted a secondary analysis including prospectively collected subjects to compare OXL profile between ethnicities in NASH.
2. Results
2.1. Subject Characteristics
The clinical and histological features of NAFL subjects from the primary cohort are presented in Table 1. The mean age in NAFL and HC, respectively, was 47 ± 15 and 43 ± 14 in HIS; 50 ± 18 and 44 ± 12 in CAU (n.s). The mean BMI in HIS was 46 ± 6 in NAFL and 26 ± 2 in HC (p-value < 0.05); in CAU, the mean BMI was 42 ± 8 in NAFL and 25 ± 3 in HC (p-value < 0.05). Within NAFL group, the mean NAS score was 3 ± 3 and 3 ± 1 for HIS and CAU, respectively (n.s). No difference in clinical and histological parameters was found between ethnicities. In subjects with NASH compared to 0-NASH, the mean NAS score was 5 ± 2 and 4 ± 1 for HIS and CAU, respectively (n.s) (data not shown).
| NAFL-HIS | NAFL-CAU | p-Value * | |
|---|---|---|---|
| Plasma, n (F/M) | 10 (7/3) | 8 (4/4) | - |
| Age (years) | 47 ± 15 | 50 ± 18 | 0.6 |
| DM, yes (%) | 4 (40) | 4 (50) | 1 |
| FBG mmol/L | 101 ± 14 | 94 ± 13 | 0.3 |
| Cholesterol (mg/dL) | 172 ± 28 | 166 ± 34 | 0.7 |
| TG (mg/dL) | 135 ± 79 | 129 ± 70 | 0.8 |
| HDL (mg/dL) | 44 ± 6 | 43 ± 8 | 0.7 |
| LDL (mg/dL) | 107 ± 24 | 99 ± 35 | 0.4 |
| HbA1c (%) | 6 ± 1 | 6 ± 1 | 0.7 |
| AST (U/L) | 31 ± 13 | 26 ± 9 | 0.6 |
| ALT (U/L) | 40 ± 27 | 31 ± 11 | 0.9 |
| Platelet | 280 ± 84 | 302 ± 110 | 0.8 |
| NAS | 3 ± 3 | 3 ± 1 | 0.7 |
| Steatosis (%) | |||
| <5% | 4 (40) | 1 (12.5) | 0.9 |
| 5 to ≤33% | 3 (30) | 5 (62.5) | |
| 34 to ≤66 % | 2 (20) | 2 (25) | |
| >66% | 1 (10.0) | 0 (0) | |
| Inflammation (%) | |||
| None | 3 (30) | 1 (12) | 0.6 |
| Mild | 1 (10) | 4 (50) | |
| Moderate | 4 (40) | 3 (38) | |
| Severe | 2 (20) | 0 (0) | |
| Ballooning (%) | |||
| None | 4 (40) | 5 (62) | 0.3 |
| Few | 5 (50) | 3 (38) | |
| Many | 1 (10) | 0 (0) | |
| Fibrosis (%) | |||
| None | 7 (70) | 7 (88) | 0.3 |
| 1A | 0 (0) | 1 (12) | |
| 2 | 1 (10) | 0 (0) | |
| 4 | 2 (20) | 0 (0) | |
The secondary cohort included NASH subject with various degrees of necroinflammation and fibrosis (Table 2). When comparing NASH-HIS and NASH-CAU, no difference was found with BMI and other clinical and histological parameters. The mean NAS score was 5 ± 2 and 5 ± 1 for NASH-HIS and NASH-CAU, respectively (n.s).
| NASH-HIS | NASH-CAU | p-Value * | |
|---|---|---|---|
| Plasma, n (F/M) | 12 (9/3) | 17 (11/6) | - |
| Age (years) | 49 ± 8 | 50 ± 13 | 0.8 |
| BMI (Kg/m2) | 41 ± 8 | 37 ± 7 | 0.6 |
| NAS | 5 ± 2 | 5 ± 1 | 0.9 |
| Steatosis (%) | |||
| <5% | 2 (17) | 0 (0) | 0.4 |
| 5 to ≤33% | 4 (33) | 7 (41) | |
| 34 to ≤66 % | 3 (25) | 4 (24) | |
| >66% | 3 (25) | 6 (35) | |
| Inflammation (%) | |||
| None | 0 (0) | 0 (0) | 0.1 |
| Mild | 2 (17) | 3 (18) | |
| Moderate | 6 (50) | 13 (76) | |
| Severe | 4 (33) | 1 (6) | |
| Ballooning (%) | |||
| None | 0 (0) | 0 (0) | 0.6 |
| Few | 8 (67) | 13 (76) | |
| Many | 4 (33) | 4 (24) | |
| Fibrosis (%) | |||
| None | 5 (41) | 4 (23) | 0.8 |
| 1a, b, c | 2 (17) | 6 (35) | |
| 2 | 2 (17) | 3 (18) | |
| 3 | 2 (17) | 2 (12) | |
| 4 | 1 (8) | 2 (12) | |
2.2. Ethnicity-Related Alterations in Plasma PUFAs and Lipid Mediator Profiles Characterize NAFL
We examined differences in plasma fatty acids and lipid mediators between NAFL and HC (Figure 1 and Table S1). Compared to corresponding HC, 25 (38% of total detected) and 7 (11%) lipid levels were different in NAFL-HIS and NAFL-CAU, respectively (FDR-adjusted p < 0.2). Ethnicity-specific changes observed in NAFL, with interaction (ethnicity × NAFL), include 8 lipids (15%) and 2 enzymatic ratios (raw p-Interaction < 0.05) but did not survive FDR-correction (q = 0.2). To rule out any effect of histological severity on the differences observed between ethnicities, the analysis was repeated on a subset of histology-matched subjects (Table S2). As a result, 12 lipids (19%) and one enzymatic ratio were found altered (raw p-Interaction < 0.05), with 3 lipids, i.e., ALA, LA, and 9-hydroperoxyoctadecadienoic acid (-HpODE) surviving FDR correction.
When compared to HC, there were overlapping differences seen in NAFL for both ethnicities as well as ethnicity-specific differences (Figure S1). In both ethnicities, NAFL showed higher levels of several eCBs derived from PUFA and other fatty acids. The 18 carbon (C18) PUFAs, ALA and LA showed similar higher trend with higher levels of downstream fatty acid alcohols, hydroperoxide, ketones, epoxides, and vicinal diols. Specific to NAFL-HIS, there were differentially higher levels of LA-triols (i.e., trihydroxyoctadecaenoic acids (TriHOMEs) (raw p-Interaction < 0.05). On histology-matched analysis, higher TriHOMEs and LA-epoxide, 12(13)-epoxyoctadecenoic acid (-EpOME) levels were found significant (raw p-Interaction < 0.05). Of note, the n-6 to n-3 ratio was higher, however, with no ethnicity × NAFL interaction. In NAFL-CAU, there was differentially higher ALA and LA and its hydroperoxide, 9-HpODE (raw p-Interaction < 0.05). On histology-matched analysis, these lipids retained significance with LA, and 9-HpODE passing the FDR-threshold.
The 20 carbon (C20) and longer chain PUFAs (LC-PUFA) showed opposite trends with higher levels in NAFL-CAU and lower in NAFL-HIS. The ratio of docosahexaenoic acid (DHA)/eicosapentanoic acid (EPA) + ALA was found lower in both ethnicities (FDR-adjusted p-value). In HIS, alcohols, ketones, thromboxane derived from arachidonic acid (AA) were higher, however with no interaction (ethnicity × NAFL). DHA and its vicinal diol, 19,20-dihydroxydocosapentaenoic acid (-DiHDoPA), levels were differentially lower (raw p-Interaction < 0.05) with tendency shown for lower EPA. On histology-matched analysis, lower levels of these lipids were found significant (raw p-Interaction < 0.05). Specific to CAU, there was differentially higher AA and its vicinal diol, 14,15-dihydroxyeicosatrienoic acid (-DiHETrE) that remained significant after histology adjustment (raw p < 0.05). While these findings show common alterations seen in NAFL for both ethnicities, they also highlight ethnicity-specific changes. This includes a divergence in LC-PUFA profile, mainly with lower EPA and DHA seen in HIS. Although these differences did not pass FDR adjustment, histology-matched analysis yielded consistent and stronger differences, suggesting ethnicity-specific differences characterized NAFL, independently of liver histology severity. It also suggests that fibrosis may weaken the differences between ethnicities.
2.4. The Progression to NASH Is Characterized by Ethnicity-Related Alterations in Plasma PUFAs and Lipid Mediator Profiles
We examined differences in plasma fatty acids and lipid mediators between 0-NASH and NASH (Figure 2 and Table S3). Compared to corresponding 0-NASH, 7 (11% of total detected) and 6 (9%) lipids were found different in NASH-HIS and NASH-CAU, respectively (raw p-value < 0.05). None passed the FDR-correction threshold. There were differentially altered lipids by ethnicity group (raw p-Interaction < 0.05), including 11 lipids (17%) and two enzymatic ratios. Three of these lipids and the two enzymatic ratios passed the FDR-correction threshold.
With the progression from 0-NASH to NASH, less marked differences in plasma PUFA profile were observed. Compared to 0-NASH, NASH-HIS showed a trend for lower plasma PUFAs, only affecting AA (raw p-value < 0.05). There was a trend for higher C18-PUFA derived alcohols, triol, epoxides and vicinal diols, with TriHOMEs being differentially higher (raw p-Interaction < 0.05). We also observed a trend for lower LC-PUFA derived lipid mediators, mainly affecting AA-alcohols, 5-hydroxyeicosatetraenoic acid (-HETE); thromboxane (TXB2); and prostaglandin (PGE2), with 5-HETE and TXB2 found differentially changed between ethnicities (raw p-Interaction < 0.05). Also, the oleic acid (OA)-derived N-oleoyl glycine levels were lower with NASH (raw p-Interaction < 0.05).
In NASH-CAU, there was a trend of higher C18, LC-PUFAs and downstream lipid mediators. Interaction (ethnicity × NASH) was shown with higher 9-HpODE, TXB2, and in EPA-epoxide, 17(18)-epoxyeicosatetraenoic acid (-EpETE) (raw p-Interaction < 0.05). There was an opposite trend for C18-PUFA derived vicinal diols that were higher in HIS and lower in CAU, compared to corresponding 0-NASH, with interaction (ethnicity × NASH) shown for 12,13- and 15,16-dihydroxyoctadecadienoic acid (-DiHODE) (raw p-Interaction < 0.05). Multiple sEH enzymatic indices were higher in HIS and lower in CAU, including 9_10-dihydroxyoctadecenoic acid (-DiHOME)/9(10)- epoxyoctadecenoic acid (-EpOME) and 9_10-DiHODE)/9(10)-epoxyoctadecadienoic acid (-EpODE) (raw p-Interaction < 0.05). Also, with NASH, there were higher levels of many eCBs, including dihomo-γ-linolenoylethanolamide, palmitoleoylethanolamide, palmitoleoylethanolamide, oleoyl-ethanolamide and N-oleoylglycine (raw p-Interaction < 0.05). Although PUFA changes are less marked in NASH, trends are consistent with changes seen in NAFL and support divergence in LC-PUFA profiles. It also highlights ethnicity-related differences in OXLs and eCBs associated with NASH progression. Given that NASH groups in both ethnicities had comparable NAS scores, this suggests the ethnicity-related differences observed with NASH are not likely driven by histological severity.
2.5. Plasma OXLs Profile Discriminates between Ethnicities with NASH
A supervised PLS-DA was performed including all profiled lipids to examine if plasma OXLs profile can discriminate between ethnicities with NASH. The model demonstrated separation between HIS and CAU with 22 (49%) contributing lipids having VIP > 1. This indicates differences in OXLs profile characterizes HIS and CAU with NASH (Figure S3 and Table S4). The Q2 and R2 for the model were 0.62 and 0.72, respectively, indicating a fair reliability of the model. Of note, an overlap between ethnic groups was observed and subjects within this area shared advanced fibrosis (grade 3 and 4), indicating that HIS and CAU subjects with NASH and advanced fibrosis share similar plasma OXLs profile. Also, it suggests that advanced fibrosis may be attenuating the multivariate model. Therefore, we repeated the analysis after excluding subjects with advance fibrosis (Figure 3 and Table S4). As a result, the model exhibited complete separation between ethnicities in NASH with 20 (44%) lipids contributing to this difference (VIP > 1.0). The Q2, R2 were 0.99 and 0.98, respectively, indicating optimal prediction and reliability of the multivariate model.
In subjects with NASH and mild to moderate fibrosis, OXLs profiles showed opposite direction of change between ethnicities in some AA-derived mediators, with TXB2, 15-keto PGE2, 5-HETE and 5,6 DiHETrE being lower in HIS compared to CAU. Many OXLs derived via LOX pathway and/or auto-oxidative routes were lower in NASH-HIS, including TriHOMEs, 9- and 13-keto-octadecadienoic acid (-KODE), 9- and 13- hydroxyoctadecatrienoic acid (-HOTE), and 5-HETE. Some CYP-derived OXLs were lower in NASH-HIS, as 9(10)-EpOME, 12(13) EpOME and 12(13)-EpODE. Multiple sEH enzymatic indices were found higher in NASH-HIS, compared to NASH-CAU, including 12,13-DiHOME/12(13)-EpOME; 9,10-DiHOME/9(10)-EpOME; 17_18-DiHETE/17(18)-EpETE along with higher levels of the vicinal diol, 14,15-DiHETrE. These findings further confirm lower AA-derived OXLs in HIS with NASH, which are also characterized by lower LOX and higher sEH-derived lipid mediators compared to CAU. They also indicate that plasma OXL profiles can discriminate between ethnicities in NASH.
3. Discussion
This study is the first to examine targeted plasma PUFA, OXL, and eCB profiles with regards to ethnicity in a group of HIS and CAU subjects with obesity and biopsy-diagnosed NAFL and NASH. Our findings indicate that: (1) NAFL and NASH are characterized by ethnicity-related differences in plasma PUFA profiles, independent of histological severity; (2) Ethnicity-related differences in plasma OXLs profiles characterize NASH, independent of histological scores; (3) Plasma PUFA profile is altered in apparently healthy HIS, independent of obesity.
The hepatic and serum/plasma PUFA profiles are dysregulated in NAFLD [16,17,18,19,20]. Our results expand on these findings and show ethnicity-related differences in plasma PUFA profile in NAFL and NASH. With NAFL, both ethnicities showed higher ALA and LA levels, which were pronounced in CAU but not in HIS. This can be attributed to differences in the levels of these PUFAs in lean and healthy subjects, as HIS showed higher levels compared to their CAU counterpart. There was also a divergence in LC-PUFA profiles between ethnicities. CAU showed higher levels mainly affecting AA, while HIS displayed lower levels mainly affecting DHA and EPA. Consistent trends were shown with the progression from 0-NASH to NASH with arachidonic almost reaching significance (p-Interaction = 0.07). This is in line with our previous untargeted profiling in NAFL done on the same subjects showing higher C18-PUFAs in CAU, and a trend for lower LC-PUFA with tendencies shown for DHA and EPA in HIS [27]. The lack of significant difference in LC-PUFA in our previous analysis may be due to the semi-quantitative nature and the clustering statistical approach. Also in agreement with our current finding is the lower serum/plasma DHA and EPA levels reported in obese HIS, compared to non-HIS [28,29]. Together, this implies diminished plasma LC-PUFA characterizes obese HIS with NAFL and NASH. As diet affects circulating and tissue PUFA levels [30], the ethnicity-dependent differences in dietary intake of n-3 PUFAs which is reportedly lower in HIS could be responsible for these observed changes [29,31]. While we did not account for diet, our findings suggest a possible etiological role for it as we observed higher linoleic and ALA levels independent of obesity in lean healthy HIS compared to CAU. Beside diet, genetic variants in cluster region of fatty acid desaturases (FADS) can predict LC-PUFA serum/blood levels [32,33]. Single nucleotide polymorphisms (SNP) in FADS1 and FADS2, which encode fatty acid desaturases, were robustly associated with NAFLD [34,35]. Lower Δ−5 desaturase levels are reported in both NAFL and NASH [16,17,20]. Notably, SNPs in FADS that are associated with insufficient LC-PUFA biosynthesis present with high frequency in Amerindians, a subgroup of HIS [36]. However, genotype was not examined in current study, and both ethnicities had lower estimated Δ−5 desaturase activity with NAFL. Therefore, diet and/or genetic factors may contribute to the observed ethnicity-related PUFA alterations but need further assessment.
Other key findings include the ethnicity-related differences in OXLs and eCBs profiles. The COX pathway exerts pro-inflammatory effects as it catalyzes the conversion of AA to prostaglandin PGE2, thromboxane TXB2, and other fatty acid alcohols [10]. In animal models of NASH, the expression and activity of COX-2 were upregulated, and its inhibition ameliorated NAFL and NASH [37,38]. Previously, high TXB2 and PGE2 levels were reported in subjects with NAFL and NASH [22]. Findings in NASH-CAU from our secondary analysis are consistent with this literature, as TXB2 and 15-Keto PGE2 discriminated between ethnicities with NASH with higher levels in CAU and lower levels in HIS. When comparing within ethnicities, the progression from 0-NASH to NASH in HIS was marked with a trend for lower AA, almost reaching statistical significance (p-Interaction = 0.07) and downstream OXLs with TXB2 being differentially lower (p-Interaction < 0.05). These findings suggest ethnicity-related alterations in AA metabolism and downstream COX-derived OXLs in NASH.
Animal studies indicate a role for LOX pathways in NAFL and inflammation [39,40]. LOX pathways lead to the synthesis of fatty acid alcohols, ketones, hydroperoxides, and the specialized pro-resolving mediators (SPMs). With possible exceptions, n-6 PUFA derived alcohols are pro-inflammatory [10]. Under oxidative stress, PUFAs can also undergo auto-oxidation to form alcohols, ketones, hydroperoxides [11]. Previous studies reported higher LOX and auto-oxidation metabolites in NAFL and increased AA metabolites via LOX with the progression to NASH [18,22,23,41]. In our results, compared to control groups, NAFL and NASH in both ethnicities presented higher alcohols and ketones derived from C18-PUFAs, indicating an upregulated LOX pathway(s). In NAFL, we observed a positive correlation between some fatty acid alcohols and the oxidative stress markers, F2-isoprostanes and 9-HETE (Table S5), implying a contribution of non-enzymatic auto-oxidation. Interestingly, our secondary analysis showed many LOX derived OXLs being higher in NASH-CAU compared to NASH-HIS, with a similar trend found for the oxidative stress marker, 9-HETE (VIP = 0.98). Together, while LOX and oxidative pathways are upregulated with NAFL in both ethnicities, the magnitude of these alterations is lesser in HIS with NASH, compared to CAU. Based on this finding, we reasoned that LOX, and possibly oxidative stress, may be pivotal for NASH severity in CAU, and to a lesser extent in HIS.
CYP enzymes catalyze the synthesis of fatty acids epoxides and alcohols. In general, fatty acid alcohols are pro-inflammatory, and epoxides are anti-inflammatory and transient, and are hydrolyzed by the action of sEH to form inactive or less active vicinal diols [10,42]. A role for sEH in NAFLD progression is indicated by animal studies, showing that sEH inhibition improves NAFL, NASH, and fibrosis [42]. In subjects with NASH, compared to NAFL, AA derived vicinal diols are higher [22]. Our results show, with the progression to NASH, an ethnicity-dependent opposite trend for vicinal diols derived from C18-PUFA, which were higher in HIS and lower in CAU. Some of these vicinal diols and sEH enzymatic indices showed interaction (ethnicity × NASH) and were found higher in NAHS-HIS and lower in NASH-CAU, compared corresponding 0-NASH. This may suggest higher activity of sEH in NASH-HIS. Our secondary analysis also shows higher ratios of multiple sEH enzymatic indices in HIS compared to CAU, and lower C18-PUFA epoxides possibly due higher hydrolysis rate. NASH-CAU showed higher levels of many PUFA epoxides, compared to NASH-HIS, indicating upregulated CYP pathway(s) and/or less hydrolysis. Together, our finding highlights ethnicity-related differences in sEH activity that was higher in HIS with NASH.
Extensive evidence from animal studies indicates a role for eCB system in NAFL, mitochondrial dysfunction and inflammation and fibrosis [43,44]. In NAFL, both ethnicities had higher levels of several eCBs. However, with the progression to NASH, many eCBs were higher in CAU and lower in HIS as compared to corresponding 0-NASH (raw p-Interaction < 0.05). We also observed levels of the OA-derived mediators N-oleoyl glycine, and oleoylethanolamide. These observations could not be examined in our secondary analysis as we detected limited numbers of eCBs and did not profile for fatty acids. Nevertheless, this may indicate ethnicity-related variations in eCBs profiles and OA metabolism with NAFLD in HIS that need to be further examined.
Our findings corroborate the epidemiological evidence indicating ethnicity as one variable affecting the association between PUFAs and cardiometabolic risks [29,45,46]. In fact, the observed ethnicity-related alterations may be relevant to NAFLD severity. EPA and DHA modulate hepatic fatty acid oxidation, de novo lipogenesis, redox balance and inflammation via direct interaction with nuclear receptors and transcription factors [47]. These LC-PUFAs are also precursors to potent SPMs which drive inflammatory resolution [10]. Also, the pro-inflammatory cascade of AA via COX is necessary for the biosynthesis of SPMs and initiating inflammatory resolution [13,48]. Therefore, a diminished level of these PUFAs may abolish anti-steatogenic and anti-inflammatory mechanisms. Likewise, a higher sEH activity may result in deactivation of anti-inflammatory PUFA epoxides [10,42]. Interestingly, our findings suggest that upregulated LOX pathway(s) may be imperative to NASH severity in CAU with a lesser extent in HIS. Collectively, we postulate that the observed ethnicity-related changes translate to the more advanced NASH histological presentation seen in HIS. Of note, these changes are independent of fibrosis or NAS scores, in fact, histology adjustment resulted in stronger differences in both analyses, implying that subjects with advanced fibrosis may share similar lipidomic profile.
Our findings have clinical/diagnostic implications. Given liver biopsy risks and limitations [49], there is an ongoing search for noninvasive biomarker for NAFLD, with multiple biomarkers have been recently proposed including betatrophin and fetuin-A [50,51]. Also, several AA- derived OXLs were shown to predict NASH including higher levels of 5- and 15-HETE, PGE2, and some vicinal diols [18,22]. While our findings in NASH-CAU show trends consistent with current literature, findings in HIS indicate otherwise. Ethnicity-related differences in plasma metabolomic profile have been reported before in diabetes, Alzheimer’s disease, and bladder cancer [52,53,54]. We propose that ethnicity-specific plasma signature may characterize NASH. In fact, utilizing ethnicity-related variations in plasma lipidomic profile may be instrumental for the enhanced precision of such diagnostic tools. If further verified, it will serve as a much-needed non-invasive tool aiding in clinical practice for early detection of NASH in both HIS and CAU populations. It can also pave the road for examination of ethnicity-specific lipidomic signatures in other ethnicities as the Asian and African American populations. On another note, a role of EPA and DHA supplementations in improving NAFLD and its risk factors is supported by clinical evidence [55,56]. Also, growing data indicate the utility of sEH inhibitors in NASH treatment [42,57,58]. Therefore, evaluating these interventions for NASH treatment seems warranted, particularly in the HIS population.
This “proof-of-concept” analysis is based on a small, single-center study. The limited sample size may have compromised the correction for multiple testing in the primary analysis. However, findings from the secondary analysis were consistent and the multivariate model is validated for overfitting and predictability. Other strengths include biopsy-characterized NAFL and NASH and analysis adjusted for BMI and histology. While NAFLD prevalence is reported to be higher in males compared to females [59,60,61], we could not examine sex differences due to small sample size. However, we did adjust for sex as a covariate.
In conclusion, we performed targeted lipidomic profiling for PUFAs and related lipid mediators with regards to ethnicity. Results show ethnicity-related divergence in LC-PUFA and downstream OXLs profiles with NAFL and NASH progression, independent of histological scores. Our secondary analysis indicates that in NASH and compared to CAU, HIS are characterized by lower levels of AA derived OXLs, lower LOX with an upregulated sEH pathway(s). These lipidomic differences may be relevant to the ethnicity-related disparity reported in NAFLD rate and severity and are worth further investigations. Our findings suggest ethnicity-specific lipidomic signature may characterize NASH. Although preliminary, these novel observations support the need for larger validation studies.
4. Subjects and Methods
4.1. Subjects and Samples
In this retrospective/ prospective cohort study (Figure S4), all subjects self-reported ethnicity as either HIS or CAU. HC subjects (n = 22) were recruited via public posts. Plasma and liver samples form bariatric surgery patients with medically complicated obesity were retrieved from the biobank repository of the Division of Gastroenterology and Hepatology, UC Davis Medical Center. The primary cohort (n = 18) consisted of subjects with NAFL and various degrees of necroinflammation. Only subjects with NASH were included in the secondary analysis (n = 9) and this cohort was expanded with prospectively collected subjects diagnosed with NASH (n = 20). Subject inclusion and exclusion criteria and details on data collection are described elsewhere [25]. Briefly, plasma samples were collected preoperatively after an overnight fast, and liver tissue samples were collected by biopsy performed during bariatric surgery. Liver histopathological evaluations were performed in a blinded fashion in the UC Davis Medical Center Department of Pathology, and samples were scored according to the NASH Clinical Research Network (NASH-CRN) histology system, The NAFLD Activity Score (NAS) and fibrosis scores were calculated [26]. NASH diagnosis was determined using a diagnostic algorithm based on steatosis, inflammation, and fibrosis scores [27]. All subjects were consented and the Institutional Review Board at the University of California, Davis approved the study protocol (# 856052).
4.2. Plasma Targeted Lipidomic Analysis
Quantitative lipidomic profiling of PUFAs, OXLs and NAEs was performed by ultra-high-performance liquid chromatography-electrospray ionization-tandem mass spectrometry (UPLC-ESI-MS/MS) (S2), as previously described [62]. Briefly, plasma samples were enriched with deuterated surrogates, isolated by liquid/liquid extraction, and separated and quantified by UPLC-ESI-MS/MS. In the primary analysis, ESI-polarity switching facilitated the simultaneous detection of eCBs (positive mode) and oxylipins and PUFAs (negative mode) on an API 6500 QTRAP (AB Sciex, Framingham, MA, USA). Metabolites were quantified against authentic analytical standards with 6-to-10-point calibration curves and calculated concentrations were corrected for analytical surrogate recovery. This method detected 5 PUFAs and 66 lipid mediators, including 10 eCBs and 46 OXLs.
The secondary analysis was performed on an API 4000 QTrap (AB Sciex) and restricted to the negative mode electrospray ionization to increase the power of the OXLs discovery. This approach detected 46 OXLs and two nitrolipids. Details on the analysis protocols and reported data are available on the Metabolomics Workbench (http://www.metabolomicsworkbench.org, accessed on 13 June 2021). ID numbers (ST000977 and ST001845). Analyses were carried out at the UC Davis West Coast Metabolomics Center. In this manuscript, abbreviations used for OXLs and eCBs follow standard consensus and are detailed with lipid identifiers in Table S6.
4.3. Statistical Analysis
Statistical analyses were performed using JMP Pro 14.1 (SAS Institute Inc., Cary, NC; http://www.jmp.com, accessed on 10 March 2021). Outliers were identified and excluded using “robust Huber M test”. Lipids with >30% missing data were excluded. Missing data were imputed by “multivariate normal imputation”. Data normality was achieved by Johnson’s transformation. After data processing, and to determine if subjects with stage 4 fibrosis (n = 2) are biological outliers, we employed principal component analysis (PCA). As a result, no outliers were detected (Figure S5).
Non-normalized data were used to calculate metabolite geometric means. Fold change (FC) was calculated for each ethnicity separately as (A − B)/B where A is the mean of (NAFL or NASH) and B is the mean of (HC or 0-NASH). A FC > 0 indicates an increase and <0 indicates decrease and ± 20% FC was set as a threshold. Student’s t-test of Johnson normalized data was used to examine differences between (NAFL vs. HC) and (0-NASH vs. NASH) in each ethnicity. Full factorial analysis of covariates (ANCOVA) was employed to evaluate the interaction of ethnicity × health status. This model included ethnicity (HIS or CAU), health status (NAFL or HC; 0-NASH or NASH), ethnicity × health status interaction as fixed effects, with age and sex as covariates. To check any effect of fibrosis or advance NAS score on the differences observed, we repeated the analysis on a subset of histology-matched subjects (n = 5 HIS and n = 5 CAU). Pathways/network visualization with fold change and p-values were plotted using Cytoscape 3.8.2 (https://cytoscape.org, accessed on 22 May 2021) [63]. Mean differences were considered likely at p < 0.05. To adjust for false discovery rate (FDR), Benjamini-Hochberg FDR correction was performed [64]. A q = 0.2 was set as a threshold, given the pilot nature and small sample size of the study.
For the secondary analysis, raw data were auto-scaled to correct for batch effect (Figure S6) [65]. Lipids affected by batch were excluded (three lipids) and data were normalized by Johnson’s transformation. Partial least square-discriminant analysis (PLS-DA) was performed to discriminate ethnicities in NASH subjects with leave- one-out cross validation (LOOCV) [66]. An R2 and Q2 > 0.5 are acceptable values to indicate reliability of the model in explaining differences between groups [67]. A variable importance in projection (VIP) score of >1.0 was set as a threshold for variable selection. To check any effect of advanced fibrosis, we repeated the analysis after excluding subjects with fibrosis grade 3 and 4 (n = 3 HIS and n = 4 CAU).
Supplementary Materials
The following supporting information can be downloaded at: https://www.mdpi.com/article/10.3390/metabo12020192/s1, Table S1: Targeted quantification data table for comparisons of HIS and CAU obese subjects with NAFL vs. lean HC, Table S2: Demographic, clinical, and histological characteristics of histology-matched subjects, Table S3: Targeted quantification data table for comparisons of HIS and CAU with NASH vs. 0-NASH, Table S4: Targeted quantification data table for comparisons of HIS and CAU with NASH in the secondary analysis, Table S5: Pearson’s correlation analysis between LOX metabolites and markers of oxidative stress, Table S6: Details on detected lipid class and identifiers, Figure S1: Venn-Diagram illustrating overlapping and unshared differences in polyunsaturated fatty acids (PUFAs) and lipid mediators between non-alcoholic fatty liver (NAFL) compared to healthy control (HC) in both ethnicities, Figure S2: Differences in plasma PUFA and lipid mediators between Hispanic (HIS) and Caucasian (CAU) in lean healthy subjects, Figure S3: Supervised multivariate clustering model demonstrating ethnicity specific oxylipins profile discriminates Hispanic (HIS) and Caucasian (CAU) with non-alcoholic steatohepatitis (NASH) with advanced fibrosis, Figure S4: Flow chart illustrating subject recruitment details for primary and secondary analysis, Figure S5: Principal component analysis (PCA) illustrating outliers before (left) and after (right) data normalization, Figure S6: Principal component analysis (PCA) illustrating the unsupervised clustering of samples from primary (red) and secondary cohort (green) shown before (left) and after (right) batch correction.
Funding
The National Institutes of Health (NIH)-West Coast Metabolomic Center (WCMC) grant number U24 DK097154, R01 DK104770 (to V.M.), 1R01 HL09133 (to P.H.) and 1R01 HL107256 (to P.H.). The U.S. Department of Agriculture (USDA) project 2032-51530-025-00D (to J.W.N). The USDA is an equal opportunity provider and employer. The content is exclusively the responsibility of the authors and does not represent the official views of NIH, WCMC or USDA.
Institutional Review Board Statement
The study was conducted according to the guidelines of the Declaration of Helsinki, and the Institutional Review Board at the University of California, Davis approved the study protocol under No. (# 856052).
Informed Consent Statement
Informed consent was obtained from all subjects involved in the study and a written informed consent has been obtained from the patient(s) to publish this paper.
Data Availability Statement
Analysis protocol and data are available at the Metabolomics Workbench (http://www.metabolomicsworkbench.org, accessed on 13 June 2021), study ID numbers (ST000977 and ST001845).
Conflicts of Interest
V.M. serves on the Advisory Board of Alexion Pharmaceuticals. The funders had no role in the design of the study; in the collection, analyses, or interpretation of data; in the writing of the manuscript, or in the decision to publish the results.
Footnotes
Footnote Group
References
Untitled section
References
- 1.Younossi Z.M., Koenig A.B., Abdelatif D., Fazel Y., Henry L., Wymer M. Global epidemiology of nonalcoholic fatty liver disease-Meta-analytic assessment of prevalence, incidence, and outcomes. Hepatology. 2016;64:73–84. doi: 10.1002/hep.28431.
- 2.Chalasani N., Younossi Z., LaVine J.E., Charlton M., Cusi K., Rinella M., Harrison S.A., Brunt E.M., Sanyal A.J. The diagnosis and management of nonalcoholic fatty liver disease: Practice guidance from the American Association for the Study of Liver Diseases. Hepatology. 2018;67:328–357. doi: 10.1002/hep.29367.
- 3.Yu J., Marsh S., Hu J., Feng W., Wu C. The Pathogenesis of Nonalcoholic Fatty Liver Disease: Interplay between Diet, Gut Microbiota, and Genetic Background. Gastroenterol. Res. Pract. 2016;2016:2862173. doi: 10.1155/2016/2862173.
- 4.Tilg H., Moschen A.R. Evolution of inflammation in nonalcoholic fatty liver disease: The multiple parallel hits hypothesis. Hepatology. 2010;52:1836–1846. doi: 10.1002/hep.24001.
- 5.Chen Z., Tian R., She Z., Cai J., Li H. Role of oxidative stress in the pathogenesis of nonalcoholic fatty liver disease. Free Radic. Biol. Med. 2020;152:116–141. doi: 10.1016/j.freeradbiomed.2020.02.025.
- 6.Rolo A.P., Teodoro J.S., Palmeira C.M. Role of oxidative stress in the pathogenesis of nonalcoholic steatohepatitis. Free Radic. Biol. Med. 2012;52:59–69. doi: 10.1016/j.freeradbiomed.2011.10.003.
- 7.Lazo M., Hernaez R., Eberhardt M.S., Bonekamp S., Kamel I., Guallar E., Koteish A., Brancati F.L., Clark J.M. Prevalence of Nonalcoholic Fatty Liver Disease in the United States: The Third National Health and Nutrition Examination Survey, 1988–1994. Am. J. Epidemiology. 2013;178:38–45. doi: 10.1093/aje/kws448.
- 8.Rich N.E., Oji S., Mufti A.R., Browning J.D., Parikh N.D., Odewole M., Mayo H., Singal A.G. Racial and Ethnic Disparities in Nonalcoholic Fatty Liver Disease Prevalence, Severity, and Outcomes in the United States: A Systematic Review and Meta-analysis. Clin. Gastroenterol. Hepatol. 2017;16:198–210.e2. doi: 10.1016/j.cgh.2017.09.041.
- 9.Mohanty S.R., Troy T.N., Huo D., O’Brien B.L., Jensen D.M., Hart J. Influence of ethnicity on histological differences in non-alcoholic fatty liver disease. J. Hepatol. 2009;50:797–804. doi: 10.1016/j.jhep.2008.11.017.
- 10.Gabbs M., Leng S., Devassy J.G., Monirujjaman M., Aukema H.M. Advances in Our Understanding of Oxylipins Derived from Dietary PUFAs. Adv. Nutr. Int. Rev. J. 2015;6:513–540. doi: 10.3945/an.114.007732.
- 11.Roy J., Le Guennec J.-Y., Galano J.-M., Thireau J., Bultel-Poncé V., Demion M., Oger C., Lee J.C.-Y., Durand T. Non-enzymatic cyclic oxygenated metabolites of omega-3 polyunsaturated fatty acid: Bioactive drugs? Biochimie. 2016;120:56–61. doi: 10.1016/j.biochi.2015.06.010.
- 12.Dennis E.A., Norris P.C. Eicosanoid storm in infection and inflammation. Nat. Rev. Immunol. 2015;15:511–523. doi: 10.1038/nri3859.
- 13.Serhan C.N. Treating inflammation and infection in the 21st century: New hints from decoding resolution mediators and mechanisms. FASEB J. 2016;31:1273–1288. doi: 10.1096/fj.201601222R.
- 14.Mallat A., Teixeira-Clerc F., Deveaux V., Manin S., Lotersztajn S. The endocannabinoid system as a key mediator during liver diseases: New insights and therapeutic openings. J. Cereb. Blood Flow Metab. 2011;163:1432–1440. doi: 10.1111/j.1476-5381.2011.01397.x.
- 15.Tam J., Liu J., Mukhopadhyay B., Cinar R., Godlewski G., Kunos G. Endocannabinoids in liver disease. Hepatology. 2010;53:346–355. doi: 10.1002/hep.24077.
- 16.Puri P., Baillie R.A., Wiest M.M., Mirshahi F., Choudhury J., Cheung O., Sargeant C., Contos M.J., Sanyal A.J. A lipidomic analysis of nonalcoholic fatty liver disease. Hepatology. 2007;46:1081–1090. doi: 10.1002/hep.21763.
- 17.Chiappini F., Coilly A., Kadar H., Gual P., Tran A., Desterke C., Samuel D., Duclos-Vallée J.-C., Touboul D., Bertrand-Michel J., et al. Metabolism dysregulation induces a specific lipid signature of nonalcoholic steatohepatitis in patients. Sci. Rep. 2017;7:46658. doi: 10.1038/srep46658.
- 18.Puri P., Wiest M.M., Cheung O., Mirshahi F., Sargeant C., Min H.-K., Contos M.J., Sterling R.K., Fuchs M., Zhou H., et al. The plasma lipidomic signature of nonalcoholic steatohepatitis. Hepatology. 2009;50:1827–1838. doi: 10.1002/hep.23229.
- 19.Kalhan S.C., Guo L., Edmison J., Dasarathy S., McCullough A.J., Hanson R.W., Milburn M. Plasma metabolomic profile in nonalcoholic fatty liver disease. Metab. Clin. Exp. 2011;60:404–413. doi: 10.1016/j.metabol.2010.03.006.
- 20.Walle P., Takkunen M., Männistö V., Vaittinen M., Lankinen M., Kärjä V., Käkelä P., Ågren J., Tiainen M., Schwab U., et al. Fatty acid metabolism is altered in non-alcoholic steatohepatitis independent of obesity. Metabolism. 2016;65:655–666. doi: 10.1016/j.metabol.2016.01.011.
- 21.Maciejewska D., Ossowski P., Drozd A., Ryterska K., Jamioł-Milc D., Banaszczak M., Kaczorowska M., Sabinicz A., Raszeja-Wyszomirska J., Stachowska E. Metabolites of arachidonic acid and linoleic acid in early stages of non-alcoholic fatty liver disease—A pilot study. Prostaglandins Other Lipid Mediat. 2015;121:184–189. doi: 10.1016/j.prostaglandins.2015.09.003.
- 22.Loomba R., Quehenberger O., Armando A., Dennis E.A. Polyunsaturated fatty acid metabolites as novel lipidomic biomarkers for noninvasive diagnosis of nonalcoholic steatohepatitis. J. Lipid Res. 2015;56:185–192. doi: 10.1194/jlr.P055640.
- 23.Feldstein A.E., Lopez R., Tamimi T.A., Yerian L., Chung Y.M., Berk M., Zhang R., McIntyre T.M., Hazen S.L. Mass spectrometric profiling of oxidized lipid products in human nonalcoholic fatty liver disease and nonalcoholic steatohepatitis. J. Lipid Res. 2010;51:3046–3054. doi: 10.1194/jlr.M007096.
- 24.Gorden D., Myers D.S., Ivanova P.T., Fahy E., Maurya M.R., Gupta S., Min J., Spann N.J., McDonald J.G., Kelly S.L., et al. Biomarkers of NAFLD progression: A lipidomics approach to an epidemic. J. Lipid Res. 2015;56:722–736. doi: 10.1194/jlr.P056002.
- 25.Kimberly W.T., O’Sullivan J., Nath A.K., Keyes M., Shi X., Larson M.G., Yang Q., Long M.T., Vasan R., Peterson R.T., et al. Metabolite profiling identifies anandamide as a biomarker of nonalcoholic steatohepatitis. JCI Insight. 2017;2:9. doi: 10.1172/jci.insight.92989.
- 26.Zelber-Sagi S., Azar S., Nemirovski A., Webb M., Halpern Z., Shibolet O., Tam J. Serum levels of endocannabinoids are independently associated with nonalcoholic fatty liver disease. Obesity. 2016;25:94–101. doi: 10.1002/oby.21687.
- 27.Mazi T.A., Borkowski K., Newman J.W., Fiehn O., Bowlus C.L., Sarkar S., Matsukuma K., Ali M.R., Kieffer D.A., Wan Y.-J.Y., et al. Ethnicity-specific alterations of plasma and hepatic lipidomic profiles are related to high NAFLD rate and severity in Hispanic Americans, a pilot study. Free Radic. Biol. Med. 2021;172:490–502. doi: 10.1016/j.freeradbiomed.2021.06.024.
- 28.Murphy R.A., Devarshi P.P., Ekimura S., Marshall K., Mitmesser S.H. Long-chain omega-3 fatty acid serum concentrations across life stages in the USA: An analysis of NHANES 2011–2012. BMJ Open. 2021;11:e043301. doi: 10.1136/bmjopen-2020-043301.
- 29.Steffen B.T., Steffen L.M., Tracy R., Siscovick D., Jacobs D., Liu K., He K., Hanson N.Q., Nettleton J.A., Tsai M.Y. Ethnicity, plasma phospholipid fatty acid composition and inflammatory/endothelial activation biomarkers in the Multi-Ethnic Study of Atherosclerosis (MESA) Eur. J. Clin. Nutr. 2012;66:600–605. doi: 10.1038/ejcn.2011.215.
- 30.Wood K., Mantzioris E., Gibson R., Ramsden C., Muhlhausler B. The effect of modifying dietary LA and ALA intakes on omega-3 long chain polyunsaturated fatty acid (n-3 LCPUFA) status in human adults: A systematic review and commentary. Prostaglandins Leukot. Essent. Fat. Acids. 2015;95:47–55. doi: 10.1016/j.plefa.2015.01.001.
- 31.Ramírez-Silva I., Villalpando S., Moreno-Saracho J.E., Bernal-Medina D. Fatty acids intake in the Mexican population. Results of the National Nutrition Survey 2006. Nutr. Metab. 2011;8:33. doi: 10.1186/1743-7075-8-33.
- 32.Schaeffer L., Gohlke H., Müller-Nurasyid M., Heid I.M., Palmer L., Kompauer I., Demmelmair H., Illig T., Koletzko B., Heinrich J. Common genetic variants of the FADS1 FADS2 gene cluster and their reconstructed haplotypes are associated with the fatty acid composition in phospholipids. Hum. Mol. Genet. 2006;15:1745–1756. doi: 10.1093/hmg/ddl117.
- 33.Koletzko B., Lattka E., Zeilinger S., Illig T., Steer C. Genetic variants of the fatty acid desaturase gene cluster predict amounts of red blood cell docosahexaenoic and other polyunsaturated fatty acids in pregnant women: Findings from the Avon Longitudinal Study of Parents and Children. Am. J. Clin. Nutr. 2010;93:211–219. doi: 10.3945/ajcn.110.006189.
- 34.Wang L., Athinarayanan S., Jiang G., Chalasani N., Zhang M., Liu W. Fatty acid desaturase 1 gene polymorphisms control human hepatic lipid composition. Hepatology. 2014;61:119–128. doi: 10.1002/hep.27373.
- 35.Vernekar M., Amarapurkar D., Joshi K., Singhal R. Gene polymorphisms of desaturase enzymes of polyunsaturated fatty acid metabolism and adiponutrin and the increased risk of nonalcoholic fatty liver disease. Meta Gene. 2017;11:152–156. doi: 10.1016/j.mgene.2016.08.009.
- 36.Yang C., Hallmark B., Chai J.C., O’Connor T.D., Reynolds L.M., Wood A.C., Seeds M., Chen Y.-D.I., Steffen L.M., Tsai M.Y., et al. Impact of Amerind ancestry and FADS genetic variation on omega-3 deficiency and cardiometabolic traits in Hispanic populations. Commun. Biol. 2021;4:1–12. doi: 10.1038/s42003-021-02431-4.
- 37.Cheng Q., Li N., Chen M., Zheng J., Qian Z., Wang X., Huang C., Xu S., Shi G. Cyclooxygenase-2 Promotes Hepatocellular Apoptosis by Interacting with TNF-α and IL-6 in the Pathogenesis of Nonalcoholic Steatohepatitis in Rats. Am. J. Dig. Dis. 2013;58:2895–2902. doi: 10.1007/s10620-013-2823-6.
- 38.Du Y., Chen J., Liu D., Bai Q., Song J., Guan J., Gao J., Liu B., Ma X. Celecoxib attenuates liver steatosis and inflammation in non-alcoholic steatohepatitis induced by high-fat diet in rats. Mol. Med. Rep. 2011;4:811–816. doi: 10.3892/mmr.2011.501.
- 39.Lazić M., Inzaugarat M.E., Povero D., Zhao I.C., Chen M., Nalbandian M., Miller Y.I., Cherñavsky A.C., Feldstein A.E., Sears D.D. Reduced Dietary Omega-6 to Omega-3 Fatty Acid Ratio and 12/15-Lipoxygenase Deficiency Are Protective against Chronic High Fat Diet-Induced Steatohepatitis. PLoS ONE. 2014;9:e107658. doi: 10.1371/journal.pone.0107658.
- 40.López-Parra M., Titos E., Horrillo R., Ferré N., González-Périz A., Martínez-Clemente M., Planagumà A., Masferrer J., Arroyo V., Claria J. Regulatory effects of arachidonate 5-lipoxygenase on hepatic microsomal TG transfer protein activity and VLDL-triglyceride and apoB secretion in obese mice. J. Lipid Res. 2008;49:2513–2523. doi: 10.1194/jlr.M800101-JLR200.
- 41.Li Q., Rempel J.D., Ball T.B., Aukema H., Minuk G.Y. Plasma Oxylipins Levels in Nonalcoholic Fatty Liver Disease. Am. J. Dig. Dis. 2020;65:3605–3613. doi: 10.1007/s10620-020-06095-8.
- 42.Warner J., Hardesty J., Zirnheld K., McClain C., Warner D., Kirpich I. Soluble Epoxide Hydrolase Inhibition in Liver Diseases: A Review of Current Research and Knowledge Gaps. Biology. 2020;9:124. doi: 10.3390/biology9060124.
- 43.Jourdan T., Djaouti L., Demizieux L., Gresti J., Vergès B., Degrace P. CB1 Antagonism Exerts Specific Molecular Effects on Visceral and Subcutaneous Fat and Reverses Liver Steatosis in Diet-Induced Obese Mice. Diabetes. 2010;59:926–934. doi: 10.2337/db09-1482.
- 44.Osei-Hyiaman U., Liu J., Zhou L., Godlewski G., Harvey-White J., Jeong W.-I., Bátkai S., Marsicano G., Lutz B., Buettner C., et al. Hepatic CB1 receptor is required for development of diet-induced steatosis, dyslipidemia, and insulin and leptin resistance in mice. J. Clin. Investig. 2008;118:3160–3169. doi: 10.1172/JCI34827.
- 45.Weir N.L., Nomura S.O., Steffen B.T., Guan W., Karger A., Klein R., Klein B.E., Cotch M.F., Tsai M.Y. Associations between omega-6 polyunsaturated fatty acids, hyperinsulinemia and incident diabetes by race/ethnicity: The Multi-Ethnic Study of Atherosclerosis. Clin. Nutr. 2020;39:3031–3041. doi: 10.1016/j.clnu.2020.01.003.
- 46.Garg P.K., Guan W., Nomura S., Weir N., Karger A.B., Duprez D., Heckbert S.R., Tsai M.Y. Plasma ω-3 and ω-6 PUFA Concentrations and Risk of Atrial Fibrillation: The Multi-Ethnic Study of Atherosclerosis. J. Nutr. 2021;151:1479–1486. doi: 10.1093/jn/nxab016.
- 47.De Castro G.S., Calder P.C. Non-alcoholic fatty liver disease and its treatment with n-3 polyunsaturated fatty acids. Clin. Nutr. 2018;37:37–55. doi: 10.1016/j.clnu.2017.01.006.
- 48.Musso G., Gambino R., Cassader M., Paschetta E., Sircana A. Specialized Proresolving Mediators: Enhancing Nonalcoholic Steatohepatitis and Fibrosis Resolution. Trends Pharmacol. Sci. 2018;39:387–401. doi: 10.1016/j.tips.2018.01.003.
- 49.Golabi P., Sayiner M., Fazel Y., Koenig A., Henry L., Younossi Z.M. Current complications and challenges in nonalcoholic steatohepatitis screening and diagnosis. Expert Rev. Gastroenterol. Hepatol. 2016;10:63–71. doi: 10.1586/17474124.2016.1099433.
- 50.Sonmez A., Dogru T., Ercin C., Genc H., Celebi G., Gurel H., Tapan S., Cicek A., Barcin C., Haymana C., et al. Betatrophin Levels Are Related to the Early Histological Findings in Nonalcoholic Fatty Liver Disease. Metabolites. 2021;11:425. doi: 10.3390/metabo11070425.
- 51.Dogru T., Kirik A., Gurel H., Rizvi A., Rizzo M., Sonmez A. The Evolving Role of Fetuin-A in Nonalcoholic Fatty Liver Disease: An Overview from Liver to the Heart. Int. J. Mol. Sci. 2021;22:6627. doi: 10.3390/ijms22126627.
- 52.Vardarajan B., Kalia V., Manly J., Brickman A., Reyes-Dumeyer D., Lantigua R., Ionita-Laza I., Jones D.P., Miller G.W., Mayeux R. Differences in plasma metabolites related to Alzheimer’s disease, APOE ε4 status, and ethnicity. Alzheimers Dementia Transl. Res. Clin. Interv. 2020;6:e12025. doi: 10.1002/trc2.12025.
- 53.Vantaku V., Donepudi S.R., Piyarathna D.W.B., Amara C.S., Ambati C.S.R., Tang W., Putluri V., Chandrashekar D.S., Varambally S., Terris M.K., et al. Large-scale profiling of serum metabolites in African American and European American patients with bladder cancer reveals metabolic pathways associated with patient survival. Cancer. 2019;125:921–932. doi: 10.1002/cncr.31890.
- 54.Van Valkengoed I.G., Argmann C., Ghauharali-van der Vlugt K., Aerts J.M., Brewster L.M., Peters R., Vaz F.M., Houtkooper R.H. Ethnic differences in metabolite signatures and type 2 diabetes: A nested case–control analysis among people of South Asian, African and European origin. Nutr. Diabetes. 2017;7:1–8. doi: 10.1038/s41387-017-0003-z.
- 55.Musa-Veloso K., Venditti C., Lee H.Y., Darch M., Floyd S., West S., Simon R. Systematic review and meta-analysis of controlled intervention studies on the effectiveness of long-chain omega-3 fatty acids in patients with nonalcoholic fatty liver disease. Nutr. Rev. 2018;76:581–602. doi: 10.1093/nutrit/nuy022.
- 56.Lu W., Li S., Li J., Wang J., Zhang R., Zhou Y., Yin Q., Zheng Y., Wang F., Xia Y., et al. Effects of Omega-3 Fatty Acid in Nonalcoholic Fatty Liver Disease: A Meta-Analysis. Gastroenterol. Res. Pract. 2016;2016:1–11. doi: 10.1155/2016/1459790.
- 57.Hammock B.D., McReynolds C.B., Wagner K., Buckpitt A., Cortes-Puch I., Croston G., Lee K.S.S., Yang J., Schmidt W.K., Hwang S.H. Movement to the Clinic of Soluble Epoxide Hydrolase Inhibitor EC5026 as an Analgesic for Neuropathic Pain and for Use as a Nonaddictive Opioid Alternative. J. Med. Chem. 2021;64:1856–1872. doi: 10.1021/acs.jmedchem.0c01886.
- 58.Lazaar A.L., Yang L., Boardley R.L., Goyal N.S., Robertson J., Baldwin S.J., Newby D.E., Wilkinson I.B., Tal-Singer R., Mayer R.J., et al. Pharmacokinetics, pharmacodynamics and adverse event profile of GSK2256294, a novel soluble epoxide hydrolase inhibitor. Br. J. Clin. Pharmacol. 2016;81:971–979. doi: 10.1111/bcp.12855.
- 59.Pan J.-J., Fallon M.B. Gender and racial differences in nonalcoholic fatty liver disease. World J. Hepatol. 2014;6:274. doi: 10.4254/wjh.v6.i5.274.
- 60.Browning J.D., Szczepaniak L.S., Dobbins R., Nuremberg P., Horton J.D., Cohen J.C., Grundy S.M., Hobbs H.H. Prevalence of hepatic steatosis in an urban population in the United States: Impact of ethnicity. Hepatology. 2004;40:1387–1395. doi: 10.1002/hep.20466.
- 61.Pramfalk C., Pavlides M., Banerjee R., McNeil C.A., Neubauer S., Karpe F., Hodson L. Sex-Specific Differences in Hepatic Fat Oxidation and Synthesis May Explain the Higher Propensity for NAFLD in Men. J. Clin. Endocrinol. Metab. 2015;100:4425–4433. doi: 10.1210/jc.2015-2649.
- 62.Pedersen T.L., Gray I.J., Newman J.W. Plasma and serum oxylipin, endocannabinoid, bile acid, steroid, fatty acid and nonsteroidal anti-inflammatory drug quantification in a 96-well plate format. Anal. Chim. Acta. 2020;1143:189–200. doi: 10.1016/j.aca.2020.11.019.
- 63.Shannon P., Markiel A., Ozier O., Baliga N.S., Wang J.T., Ramage D., Amin N., Schwikowski B., Ideker T. Cytoscape: A software environment for integrated models of Biomolecular Interaction Networks. Genome Res. 2003;13:2498–2504. doi: 10.1101/gr.1239303.
- 64.Benjamini Y., Yekutieli D. Quantitative Trait Loci Analysis Using the False Discovery Rate. Genetics. 2005;171:783–789. doi: 10.1534/genetics.104.036699.
- 65.Van den Berg R.A., Hoefsloot H.C., Westerhuis J.A., Smilde A.K., van der Werf M.J. Centering, scaling, and transformations: Improving the biological information content of metabolomics data. BMC Genom. 2006;7:142. doi: 10.1186/1471-2164-7-142.
- 66.Westerhuis J.A., Hoefsloot H.C., Smit S., Vis D.J., Smilde A.K., Van Velzen E.J.J., Van Duijnhoven J.P.M., Van Dorsten F.A. Assessment of PLSDA cross validation. Metabolomics. 2008;4:81–89. doi: 10.1007/s11306-007-0099-6.
- 67.Eriksson L., Byrne T., Johansson E., Trygg J., Vikström C. Multi-and Megavariate Data Analysis Basic Principles and Applications. Volume 1 Umetrics Academy; Umeå, Sweden: 2013.
Associated Data
Supplementary Materials
Data Availability Statement
Analysis protocol and data are available at the Metabolomics Workbench (http://www.metabolomicsworkbench.org, accessed on 13 June 2021), study ID numbers (ST000977 and ST001845).