Next Article in Journal
Statistical Optimisation of Phenol Degradation and Pathway Identification through Whole Genome Sequencing of the Cold-Adapted Antarctic Bacterium, Rhodococcus sp. Strain AQ5-07
Next Article in Special Issue
Mechanism of Action of Atypical Antipsychotic Drugs in Mood Disorders
Previous Article in Journal
Light Pollution and Cancer
Previous Article in Special Issue
Enhancement in Phospholipase D Activity as a New Proposed Molecular Mechanism of Haloperidol-Induced Neurotoxicity
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Understanding Mechanisms Underlying Non-Alcoholic Fatty Liver Disease (NAFLD) in Mental Illness: Risperidone and Olanzapine Alter the Hepatic Proteomic Signature in Mice

1
Department of Biomedical Sciences, College of Osteopathic Medicine, University of New England, Biddeford, ME 04005, USA
2
Department of Psychology, University of Southern Maine, Portland, ME 04103, USA
3
Center for Molecular Medicine, Maine Medical Center Research Institute, Scarborough, ME 04074, USA
*
Author to whom correspondence should be addressed.
Int. J. Mol. Sci. 2020, 21(24), 9362; https://doi.org/10.3390/ijms21249362
Submission received: 2 November 2020 / Revised: 1 December 2020 / Accepted: 3 December 2020 / Published: 8 December 2020

Abstract

:
Patients with severe mental illness have increased mortality, often linked to cardio-metabolic disease. Non-alcoholic fatty liver disease (NAFLD) incidence is higher in patients with schizophrenia and is exacerbated with antipsychotic treatment. NAFLD is associated with obesity and insulin resistance, both of which are induced by several antipsychotic medications. NAFLD is considered an independent risk factor for cardiovascular disease, the leading cause of death for patients with severe mental illness. Although the clinical literature clearly defines increased risk of NAFLD with antipsychotic therapy, the underlying mechanisms are not understood. Given the complexity of the disorder as well as the complex pharmacology associated with atypical antipsychotic (AA) medications, we chose to use a proteomic approach in healthy mice treated with a low dose of risperidone (RIS) or olanzapine (OLAN) for 28 days to determine effects on development of NAFLD and to identify pathways impacted by AA medications, while removing confounding intrinsic effects of mental illness. Both AA drugs caused development of steatosis in comparison with vehicle controls (p < 0.01) and affected multiple pathways relating to energy metabolism, NAFLD, and immune function. AA-associated alteration in autonomic function appears to be a unifying theme in the regulation of hepatic pathology.

1. Introduction

Patients with severe mental illness, including schizophrenia, experience a significantly shortened life expectancy with high rates of all-cause mortality [1,2,3], most commonly associated with metabolic dysregulation and cardiovascular disease. Non-alcoholic fatty liver disease (NAFLD) is a growing epidemic worldwide and is associated with obesity, insulin resistance, and inflammation [4,5,6,7]. Incidence of NAFLD is elevated in patients with severe mental illness [3] and with atypical antipsychotic (AA) treatment [8]. NAFLD is progressive in nature, characterized by hepatic steatosis and fibrosis, and can culminate in cirrhosis and hepatocarcinoma if not managed [9,10,11]. NAFLD is the most prevalent chronic liver disease in the world, and given the exponential rise in incidence of NAFLD [12], this disease is now one of the leading indications for hepatic transplant [13]. Understanding drug-associated effects on NAFLD incidence and progression is therefore of growing importance.
Among patients with severe mental illness, metabolic dysregulation can be due to intrinsic factors as well as extrinsic factors, including antipsychotic medications [14]. Although the clinical literature describes relative risks and contributions of intrinsic and extrinsic effects on metabolic disorders in schizophrenia [15,16,17], underlying mechanisms in both schizophrenic and non-schizophrenic patients remain elusive. Patients prescribed antipsychotic medications experience significant endocrine and metabolic side effects including weight gain, insulin resistance, hyperglycemia, dyslipidemia, hyperprolactinemia, and bone loss, regardless of diagnosis [6,18]. Broadly stated, AA side effects appear to be mediated both by global (i.e., the central nervous system) and tissue-specific mechanisms, including those in the liver.
Multiple groups have explored underlying mechanisms using preclinical models in order to avoid confounding intrinsic risks present in patient populations. Given the tremendous complexity of the known AA pharmacology as potent antagonists at multiple G-protein coupled receptors (GCPR), as well as the complex endocrine, metabolic, and immune side effects observed clinically and pre-clinically across the lifespan, we chose to take a proteomic approach to evaluate the effects of AA medications on protein expression profiles in liver. The goal of this study was to elucidate AA-regulated pathways associated with the development of NAFLD. To address this aim, we chose two widely prescribed AA drugs with divergent pharmacology [19] and differences in metabolic risk [20]: risperidone (RIS) and olanzapine (OLAN). OLAN is associated with high clinical risk for obesity, dyslipidemia, and other aspects of insulin resistance/metabolic syndrome, while RIS is known to have intermediate metabolic liability.
We provide the first report (to our knowledge) of hepatic proteomic signatures following sub-chronic treatment with clinically relevant doses of RIS and OLAN, and report their diversity from the hepatic proteome of drug vehicle (VEH)-treated controls using healthy adult mice. These data identify canonical pathways impacted with AA treatment that are linked to metabolic and immunologic dysregulation and are consistent with drug-associated increase in the incidence of NAFLD and enhanced sympathetic tone.

2. Results

2.1. Animal Health

All animals used in this study were healthy and displayed no adverse effects of daily gavage for the 28-day study duration. Details of animal growth rates, locomotor activity, adiposity, and energy expenditure for this cohort were recently reported in a manuscript describing AA effects on the heart [21]. Overall, the animals in this study showed normal feeding behavior, weight gain, and locomotor activity. Although there was no treatment effect on body weight gain, animals treated with RIS or OLAN showed a significant increase in adiposity compared to VEH treatment as well as a reduction in plasma-free fatty acid concentrations [21]. Importantly, animals were maintained on standard mouse chow to avoid potentially confounding effects of high-fat feeding, an experimental paradigm associated with diet-induced obesity and NAFLD [22,23]. Consistent with this dietary regimen, there was no treatment effect on concentrations of total plasma triglycerides (299 µM ± 55.1, 291 µM ± 40.0, 247 µM ± 21.6; p > 0.05) or total cholesterol (200 µM ± 13.4, 205 µM ±7.93, 213 µM ± 10.7; p > 0.05) in VEH, RIS, or OLAN mice, respectively.

2.2. Atypical Antipsychotic Treatment Induced Histopathologic Changes Consistent with the Development of NAFLD

Moderate duration AA treatment induced microsteatosis, macrosteatosis, ballooning hepatocytes, inflammatory foci, and Mallory–Denk body formation in healthy mice, as shown across two representative field views in Figure 1A. Significant differences in histopathologic lesion scores between RIS or OLAN and VEH controls were found for microsteatosis, Mallory–Denk body formation, and formation of inflammatory foci (Figure 1B). Significant differences in hepatocyte ballooning and macrosteatosis were found for RIS and OLAN, respectively. There was early evidence of periportal fibrosis occurring sporadically in OLAN-treated mice, although this did not rise to statistical significance. Lesions generated by RIS and OLAN were largely consistent with early NAFLD.

2.3. Atypical Antipsychotic Treatment Altered the Hepatic Proteome

We evaluated the effect of sub-chronic (28 day) AA treatment on hepatic protein expression using an unbiased sequential window acquisition of all theoretical spectra (SWATH) proteomic approach to characterize the unique protein signature resulting from RIS or OLAN treatment, compared to VEH (see Supplementary File S1). The proteomic signature data were then curated by analyzing the Kyoto Encyclopedia of Genes and Genomes (KEGG) database pathway maps for the role of each protein and the resulting changes in predicted physiological functions, which were cataloged and sorted for up- or downregulation relative to vehicle controls (see Supplementary Files S2–S4). Importantly, the proteomic analysis was conducted using liver samples from the same animals that we evaluated histologically (Figure 1), and thus proteomic results are coincident with AA-associated NAFLD.
AA treatment was associated with a large number of proteins that were differentially expressed (DE) in the livers of AA-treated mice relative to VEH control (Figure 2A). Livers from RIS mice showed the greatest overall number of proteins changed, with 69 proteins upregulated and 76 proteins downregulated vs. VEH (p < 0.05). For OLAN mice, expression of 45 proteins were upregulated, while 46 proteins were downregulated (p < 0.05). Further analysis revealed that 20 were similarly changed in RIS and OLAN mice (up- or downregulated; Figure 2A), while no proteins were differentially changed (increased by one drug but decreased by the other).
The DE proteins common between the two treatments were predicted to affect mostly metabolic, mitogenic (growth signaling), and immune functions, on the basis of KEGG pathway analysis (Figure 2B). The highest cumulative number of predicted functional changes caused by both RIS- and OLAN-induced protein dysregulation were metabolic functions, especially disruptions in the ATP production pathways of glycolysis and mitochondrial oxidative phosphorylation. Similarly, many of the same disrupted ATP-generating pathways also predicted decreased amino and nucleic acid synthesis as a result of the dysregulated proteins. Interestingly, some predicted pathway functions related to mitogenic signaling, such as survival (resistance to apoptosis), were increased, and with them the predicted functions that indicated increased probabilities of various cancers, according to the KEGG pathway maps. Additionally common between the two AA treatments were predicted increases in generation of reactive oxygen species and oxidative stress, resulting from the protein changes (Table 1 and Table 2).
Setting a cutoff of 10 combined changes (“up” + “down”) for each KEGG pathway-predicted function, we found that OLAN caused 24 and RIS caused 17 functional pathway changes (Table 1 and Table 2). The complete list can be found in Supplementary File S4. Predicted functional pathway changes with OLAN treatment also fell into a larger number of categories than RIS, including immune effects, gastrointestinal (GI) functions, intracellular functions, and effects on regulation of blood pressure.
Many of the predicted functional changes resulting from the AA treatment-induced protein dysregulation were differentially altered (increased by one drug but decreased by the other), likely due to differences in pharmacological properties of the two medications [19,24]. For example, lipolysis pathway function was more frequently downregulated by OLAN, but was upregulated by RIS. Similarly, pathway changes that predicted increased cell migration by RIS were decreased by OLAN (Table 1 and Table 2).
As with the 20 DE proteins common to RIS and OLAN (Figure 2 panel B), the predicted functions of the other 216 dysregulated proteins also predominantly affected metabolic functions, with predicted reductions in ATP-generating pathways (glycolysis, mitochondrial oxidative phosphorylation, and generation of pyruvate). Nucleic acid synthesis and amino acid synthesis were also predicted to be reduced more frequently, as a result of the up- or downregulation of the proteins. Both RIS and OLAN showed an increase in predicted pathways that led to lipogenesis, which is in line with the observed increases in liver steatosis (Figure 1).
We found that there were some differentially expressed proteins (p < 0.05) that were represented in a large and diverse number of KEGG pathway maps (Table 3). Using cutoff criterion of a minimum of five representative KEGG pathways, we found that RIS had 25 proteins (16 upregulated, 9 downregulated) and OLAN had 15 proteins (6 upregulated, 9 downregulated) represented in multiple pathways. These impactful DE proteins featured prominently in signal transduction pathways, mitochondrial functions, and glycolysis.

2.4. Antipsychotic Treatment Altered Protein Expression in Hepatic Energy Metabolism Pathways

AA medications alter whole-body energy metabolism and insulin sensitivity both clinically and pre-clinically, including liver-specific effects such as AA-associated increases in hepatic glucose production [25,26,27]. In this study, dysregulation of energy metabolism was a prominent predicted outcome of the AA-induced DE in the liver proteome. Utilization of glucose, lipids, and amino acids in various pathways were disturbed, and this contributed to insulin resistance.
Many of the DE proteins were predicted to impact mitochondrial energy metabolism, including oxidative phosphorylation, the tricarboxylic acid cycle, and the electron transport chain. RIS affected multiple proteins within mitochondria, including COX6B1, ACAA2, BCAT2, ALDH7A1, CS, ECI2, MECR, NDUFA8, NDUFB8, and NDUFS1, and OLAN affected COX6B1, CYCS, and AK2. Some of the proteins affected mitochondrial function from outside of the organelle, such as ACLY and MTHFD1 by RIS, and DECR2 and ENO3 by OLAN.
RIS and OLAN treatment caused different patterns of DE protein expression relating to glycolytic pathways, with RIS effects promoting increases in glycolytic flux pathways mainly by shifts away from mitochondrial function and increased pro-growth pathways, but OLAN had reduced glycolytic pathway modulation, mainly due to decreased glycolytic enzymes (Table 3).
The most frequent contributing factor to the predicted overall decrease in energy generation for both RIS and OLAN were the DE proteins predicting reduced generation of pyruvate via metabolism of amino acids and other metabolites. With RIS, this included downregulation of AGXT, ALDH7A1, ENO3, CS, ABAT, TALDO1, and AOX3, and with OLAN the downregulation of ENO3, PNP, AMY2B, VAPA, AK2, and ATP2B4.
Insulin resistance can lead to hyperglycemia, and several of the DE proteins in both AA treatment groups were associated with pathways linked to disrupted circulating glucose homeostasis and insulin resistance. RIS induction of INPP5D was predicted to decrease insulin receptor signaling, and its downregulation of ALR1 was predicted to affect multiple upstream glucose processing pathways. RIS also increased ADIPOQ, whose adipokine functions are predicted to improve insulin sensitivity. OLAN’s downregulation of the two important glycolytic enzymes ENO3 and AK2 were predicted to increase hyperglycemia.
Ectopic lipid accumulation (lipotoxicity) is associated with insulin resistance [28,29,30], and predicted dysregulation of lipid metabolism was evident in both the RIS and OLAN treatment groups. Lipotoxicity was predicted by RIS treatment via changes in processing of fatty acids and coenzyme-A by the DE proteins NUDT7, AOX3, ACAA2, and ECI2. OLAN’s effects on lipotoxicity were predicted by changes to fatty acid oxidation and cholesterol transport by the DE proteins CYCS, COX6B1, DECR2, NPC2, and GNAI1. RIS also had multiple predicted interactions with the peroxisome proliferator-activated receptors (PPARs), which are transcription factors regulating expression of genes involved in various aspects of metabolism including lipid metabolism [31,32]. RIS increased ADIPOQ, which is associated with regulation of adipocyte differentiation, lipid uptake, and insulin sensitivity [33,34,35], as well as SMARCC1, a co-transcription factor of PPARγ. RIS also increased SLC27A1, whose fatty acid transport function is predicted to activate all PPARs (alpha, beta/delta, and gamma). At the RNA level, OLAN treatment significantly altered expression of 26 genes (p < 0.05) that are regulated by PPARs (Supplementary File S5). Notably, FABP3, lipoprotein lipase and UCP1 were upregulated > 3-fold vs. VEH.

2.5. Antipsychotic Treatment Altered Expression of Proteins Associated with NAFLD

Some of the DE proteins among the RIS and OLAN lists were those with predicted functional effects associated with liver steatosis, fibrosis, and NAFLD (Figure 3). Metabolic syndrome and insulin resistance contribute greatly to the progression of AA-mediated liver disease, including steatosis and fibrosis [8,36]. Among the KEGG pathway-predicted functional changes, metabolism featured prominently for both RIS and OLAN.
Several DE proteins are linked to pathways that regulate fibrosis. In the RIS treatment group, increased PRKCA signaling was predicted to increase activation of the Transforming Growth Factor-beta (TGF-β) pathway, leading to increased matrix deposition, and decreased AOX3 would lead to decreased retinoic acid, both leading to increased fibrosis. In the OLAN treatment group, increased PLOD3 activity was predicted to increasingly hydroxylate and glycosylate collagen motif-containing proteins, leading to cross-linking and fiber formation [37].
The predicted effects of the DE of some proteins also indicated potential insulin resistance, an important cause of NAFLD. With RIS treatment, the decreased expression of AOX3 and increased expression of INPP5D are predicted to contribute to insulin resistance. With OLAN treatment, in addition to its induction of INPP5D, the suppression of AMY2B and GNAI1 are predicted to increase hyperglycemia and decrease insulin sensitivity. Some of the DE proteins predicted to affect mitochondrial function also cause dyslipidemia, including MECR, CS, ACLY, AOX3, and COX6B1 by RIS, and CYCS, COX6B1, and ENO3 by OLAN.
RIS and OLAN were also predicted to impact regeneration and repair pathways, but their functional profiles were different. One important predicted regenerative function was mitochondrial biogenesis, which was essentially absent in the OLAN proteome profile, but promoted by RIS-mediated increase of ADIPOQ, GNAS, PRKCA, PDGFRB, SMARCC1, and SLC27A1.
The suppression of apoptotic pathways led both AAs to have a predicted increase in cancer incidence or severity. RIS had more increases in proteins with pro-growth and pro-mitogenic signaling effects, such that both proliferation and survival were increased. These included increased receptor tyrosine kinases such as EPHA2 and PDGFRB, as well as GNAS (Gs G-protein) and PRKCA, which were prominently represented in many KEGG pathways. There were also protein DE that predicted suppression of apoptosis with RIS, such as GSTA3 (increased glutathione activity), and IL1RAP. OLAN, however, saw increases in predicted survival pathways mostly by suppression of apoptosis via decreased expression of CYCS and GNAI1 (Gi G-protein), as well as increased G6PD, which contributes to increased glutathione production.
The pathways that indicated increases in cancer incidence and severity (see Table 2 and Table 3) followed the same distinct trends for the two drugs, with RIS promoting growth and survival, and OLAN decreasing apoptotic pathways. However, both RIS and OLAN treatment caused increase in proteins that predicted DNA damage (and impairment in DNA repair). For RIS, these were from the downregulation of MTHFD1 and ATIC, and for OLAN, the downregulation of CTC1 and TOP2B. Some of the predicted functions that indirectly promoted cancer severity included angiogenesis, metastasis, and drug resistance.

2.6. Antipsychotic Treatment Disrupted Hepatic Immune/Inflammatory Signaling Pathways

Insulin resistance is an inflammatory condition [38,39,40], and disruption of the liver’s metabolic system activates immune and inflammatory pathways that can further exacerbate metabolic perturbations that feed forward to NAFLD (Figure 4). These perturbations include tissue damage, steatohepatitis, and fibrosis [41,42]. Hepatocyte; hepatic stellate cell; sinusoidal endothelial cell; and Kupffer cell stress responses to lipotoxicity, metabolic dysregulation, and oxidative stress can all increase inflammatory cytokine release, innate intracellular immune reactivity, and immune cell recruitment and activation [43,44,45,46].
Both RIS and OLAN increased predicted activation of inflammatory pathways (Figure 4). Various intracellular mechanisms predicted to be activated by the dysregulated proteins contributed to this activation, including increased Nuclear Factor kappa-B (NF-κB) expression, increased expression of receptors for advanced glycation end products (RAGE) and pattern-recognition receptors, increased intracellular calcium, increased stress fiber formation, increased senescence, and increased prostaglandin formation. RIS-mediated upregulation of PRKCA, CTSS, GNAS, ROCK2, KLKB1, ATIC, IL1RAP, DGKZ, HSPB1, and RNF31, and OLAN-mediated upregulation of GLB1, SGPL1, CD68, and FMOD predicted increases in inflammatory cytokine production. Some of the dysregulated proteins also had predicted anti-inflammatory functions, which may represent local attempts to restore tissue homeostasis. These DE proteins were fewer in number, and included upregulated GSTA3, and downregulated CARD9, LAP3, and PSME1 by RIS, and upregulated G6PD and INPP5D, and downregulated PIP4P2 and ELAVL1 by OLAN.
Pathways that impact immune cell activation and migration were also predicted to be dysregulated by RIS and OLAN, including functions such as antigen processing, phagocytosis, Ras and MAP kinase signaling, and Nod signaling. RIS-induced upregulation of AFDN, LAMP1, PRKCA, and IL1RAP, and OLAN-induced upregulation of CFL1 and downregulation of CNTN1 and GNAI1 were predicted to increase immune cell activation and migration. Some of the dysregulated proteins also countered these predicted immune cell activations, such as RIS-induced downregulation of MSN and SNRPB, and OLAN-induced upregulation of INPP5D.
Some predicted intracellular functions serving as innate immune responses were also changed. These included predicted responses to numerous infectious diseases, necroptosis, antigen processing, immunogenic phospholipids, and protein misfolding. RIS-mediated downregulation of HYOU1, RPN2, GANAB, DDOST, ATIC, and AOX3, and upregulated RNF31, and OLAN-mediated upregulation of G6PD and NUP98, and downregulation of SEPT2, ACE, CYCS, and NCL were predicted to affect these functions.
A few of the dysregulated proteins were also predicted to affect changes in platelet and complement activation, potentially affecting coagulation and thrombotic conditions [47] as well as innate immune responses. RIS-mediated upregulation of KLKB1 and MBL1 (also increased by OLAN), and OLAN-mediated downregulation of GNAI1 were predicted to impact the coagulation cascade, platelet activation, and activation of the complement cascade.

2.7. Antipsychotic Treatment Affected Pathways Regulated by the Autonomic Nervous System: A Unifying Mechanism

The autonomic nervous system (ANS) exerts control over multiple physiological pathways in the liver, including those relating to insulin resistance and NAFLD [48]. Furthermore, AA medications have been implicated in regulating sympathetic nervous system (SNS) tone, and AA effects on heart and bone appear to be due at least in part to enhanced SNS activity [21,49,50,51]. In addition to pathways relating to insulin resistance, energy metabolism, fibrosis, and regeneration/repair described above, AA treatment altered the expression of proteins in other pathways known to be regulated by the ANS, including those relating to biliary physiology [48,52].
RIS and OLAN each had multiple DE proteins predicted to interfere with bile synthesis, secretion, and ductal motility. RIS-mediated DE of GNAS, BLVRB, AFDN, VAPA, LRP1, and HDS3B7 were all predicted to decrease proper processing and function of bile. With OLAN treatment, bile dysfunction was predicted due to the dysregulation of GNAI1, NPC2, VAPA, GLB1, AQP1, ATP2B4, TNNC2, and P2RX4.
Given the overarching regulatory role of the ANS in whole-body and hepatic functions linked to antipsychotic side effect profiles, AA-associated enhanced SNS tone is a powerful unifying regulatory system for drug associated metabolic side effects including NAFLD (Figure 5).

3. Discussion

Patients prescribed AA medications carry significant metabolic and cardiovascular risk burden, across the lifespan and regardless of diagnosis. Although the cardiovascular/metabolic risk varies in severity across the drug class, given the large patient population (on- and off-label prescribing) and added cardiovascular risk associated with severe mental illness, the overall public health risk and impact is high. Incidence of NAFLD is comorbid with obesity, insulin resistance/metabolic syndrome, as well as with severe mental illness [3,10]. NAFLD is a distinct risk factor for cardiovascular disease, and AA medications have been reported to increase NAFLD, clinically. Despite the extensive clinical literature characterizing the adverse endocrine and metabolic side effects of AA medications, the underlying molecular/pharmacologic mechanisms remain elusive. In order to specifically address the role of AA medications in the pathogenesis of NAFLD, we chose a proteomic approach, given the complexity of AA pharmacology as well as the complexity of cardiovascular disease and its underlying risk factors. This approach has been employed by others to evaluate the multifactorial etiology of NAFLD progression spanning multi-organ, multi-system effects as well as comorbidities and sex differences associated with NAFLD [53].
To our knowledge, this is the first report that examines the impact of moderate duration, clinically relevant dosing of AA medications on the liver proteome. AA medications display complex pharmacology, potently antagonizing diverse G protein-coupled receptors and eliciting complex phenotypes clinically; thus, a pathways analysis approach is needed to help tease out canonical and tissue-specific effects associated with these medications. Use of a preclinical model to evaluate potential mechanisms correlated with clinical phenotypes is a powerful approach as it (a) removes confounding effects of the disease on hepatic function/metabolic disease; (b) allows control of extrinsic factors such as nutrition, which is a known contributor to NAFLD (high fat diet); and (c) allows evaluation of complex, interconnected pathways involved in the chronic disease phenotypes induced by antipsychotic treatment. Our intentional use of C57BL/6 J mice fed a low-fat (chow) diet was to minimize the possibility that any AA-induced changes could be attributed to confounding effects of diet on dyslipidemia and metabolic disease.
Taking a proteomic approach allows for a broad, wide-ranging analysis in order to begin proposing mechanisms of AA-induced NAFLD. It should be noted that the DE proteins measured are almost certainly an underestimation of those affected, because this approach does not account for posttranslational modifications. Many of the pathways highlighted are regulated and mediated by phosphorylation events, and thus changes in fine-scale rheostasis of a pathway would not be detected by this approach. However, overt changes in protein level in AA-treated animals relative to VEH-treated controls are almost certainly biologically meaningful, particularly when the proteins themselves initiate downstream cascades. One such protein that is not activated by phosphorylation, mannose-binding lectin (MBL1), was significantly overexpressed during both RIS and OLAN treatment. The downstream impact of MBL1 interactions include activation of numerous phosphorylation-dependent signaling cascades (e.g., NF-κB signaling), some of which were not seen as changes in total proteins but whose impact could be seen phenotypically in the current report (Figure 1), in our prior reports demonstrating global immune dysregulation during AA treatment [54,55], and in several small-scale studies of patients taking AA medications [56]. Similarly, significantly depressed levels of cytochrome oxidase subunit VI B (COX6B) during both RIS and OLAN treatment would lead to profound changes in oxidative phosphorylation and metabolic rate that may not be reflected as level changes of downstream proteins because differentially phosphorylated states are not detected. However, metabolic changes in AA-treated patients and our previous report phenotypically demonstrate mitochondrial dysfunction that is highly consistent with reduced oxidative phosphorylation [21]. Finally, DE by both RIS and OLAN of proteins whose function depends on their phosphorylation state are highly likely to be relevant in the absence of posttranslational modification detection when the predicted phenotype is consistent with known AA-induced side effects. One such protein is smoothelin (SMTN1), which mediates contractility of smooth muscle such as that found in the vasculature. Significant changes in SMTN1 would be predictive of altered vascular contractility, leading to changes in blood pressure, and acute-onset changes in blood pressure are known side effects of AA, including both OLAN and RIS. All of these proposed mechanisms required further analysis to determine the impact of posttranslational modifications; however, the goal of this report was to describe an overarching framework with which to generate and evaluate specific mechanistic hypotheses addressing the development of AA-induced NAFLD.
It is well documented in the clinical literature that patients treated with AA medications often have increased adiposity, insulin resistance, dyslipidemia, and other related endocrine and metabolic side effects (see reviews [15,57,58]). Medication effects vary in severity across the AA drug class, with OLAN and clozapine noted for significant metabolic liability while RIS has intermediate metabolic liability [20]. These effects are observed across the lifespan in patients with severe mental illness, including schizophrenia, as well as in patients who are treated for other reasons [3,57,59]. This results in significant diminution in quality of life, shortening of lifespan, and increased burden on the healthcare system. Here, we focused on AA effects on the liver as insulin resistance is a pivotal factor underlying NAFLD progression [3,10]. The experimental design we employed culminated in significant hepatic steatosis, hepatocyte ballooning, formation of Mallory–Denk bodies, and development of inflammatory foci consistent with NAFLD (Figure 1). Coincident with this liver histopathology, we quantified significant changes in the liver proteomes of these mice, consistent with NAFLD. Low-dose, moderate duration of AA treatment resulted in a large number of DE proteins in RIS and OLAN mice. In this report, we chose to highlight DE profiles of proteins involved in pathways that are reported to be involved in NAFLD pathophysiology, including those regulating insulin sensitivity, energy metabolism, inflammatory/immune function, biliary function, fibrosis, and carcinogenesis.
It is not surprising, given the well-documented metabolic side effects observed clinically, that canonical pathways relating to energy homeostasis (including glucose, lipid, and amino acid metabolism and regulation of insulin sensitivity) were significantly altered in the livers of mice treated with RIS or OLAN. What is somewhat surprising, however, is the sheer number of DE proteins and canonical metabolic pathways perturbed by AA treatment. The magnitude of AA-associated effects we report here, observed in the absence of other extrinsic effects (such as high-fat feeding) or intrinsic effects of severe mental illness, provides mechanistic context to the well-documented, rapid, and profound metabolic dysregulation observed both clinically and preclinically.
There is growing evidence that schizophrenia has an inflammatory component [59], and previous clinical reports have suggested altered cytokine profiles in patients taking AAs [60,61,62]. This has led to the hypothesis that dampening of inflammatory responses may contribute to the efficacy of AA medications for psychosis. A unified mechanistic understanding of immune dysregulation has been difficult to decipher from these clinical reports because they are largely confounded by examining patients with psychiatric diagnoses and thus atypical inflammatory backgrounds. We previously reported global immune dysregulation with and without challenge, as measured by circulating cytokines and antibody levels using our preclinical murine model [54,55]. This current study bolsters these findings by identifying tissue-specific inflammatory dysfunction that could lead to NAFLD in the context of histopathologic findings consistent with its development.
Drug effects on hepatic function are likely to be both indirect (CNS mediated) and direct (tissue/cell mediated), as AA target receptors are expressed in the periphery as well as in the CNS. Direct and indirect regulation of AA-associated effects have been reported using preclinical models evaluating AA-associated whole-body and hepatic insulin resistance [25,26,27] and AA-associated effects on bone biology [63,64]. In the case of the liver, AA receptor targets are expressed on hepatocytes, and our proteomic analysis shows alterations in expression of proteins associated with GCPR signaling, consistent with direct effects of AA on the liver. This was most noticeable in the direct upregulation of the Gs G-protein (GNAS) by RIS, and the direct downregulation of the Gi G-protein (GNAI1) by OLAN, both of which were featured in numerous KEGG pathway maps. Globally impaired metabolic pathways, especially deficient production of amino acids, nucleic acids, and ATP, also limit the mitigation of NAFLD by inhibiting cell replication and repair mechanisms [65].
The liver is innervated by both afferent and efferent autonomic nerves, and myriad hepatic functions are regulated by the autonomic nervous system [48,66,67]. ANS-mediated hepatic functions include regulation of glucose and lipid metabolism, nutrient sensing, ion sensing and osmotic pressure, drug metabolism (CYP enzymes), blood flow, the biliary system, circadian rhythm, and hepatic regeneration and repair. Furthermore, the ANS has been implicated in the pathogenesis and progression of NAFLD [7]. Autonomic dysfunction is a hallmark of psychiatric as well as metabolic disease [68,69,70,71], and some side effects associated with AA therapy are consistent with medication-associated increase in sympathetic tone [50,58]. Some AA drugs have been reported to act via SNS mechanisms to cause adverse effects on glycemia [72], cardiac function [21,49,73], and bone biology [51,63]. In this report, both RIS and OLAN cause robust histological and proteomic changes in the liver that are consistent with NAFLD and dysregulation of myriad pathways associated with NAFLD and known to be regulated by the ANS. Thus, AA-associated disruption in ANS function is likely to be a pivotal overarching cause of AA-associated NAFLD (Figure 5).
In summary, low dose, moderate duration AA treatment leads to hepatic steatosis and fibrosis, preclinically, and is associated with hepatic proteomic signatures consistent with NAFLD. Patients with NAFLD have enhanced sympathetic nervous system activity (as do patients with schizophrenia and patients with diabetes). AA medications also are associated with increased sympathetic tone that has been implicated in AA-associated side effects of tissues including bone and heart, and effects on the liver proteome are consistent with increased sympathetic tone. Given the high prevalence of metabolic disease, including NAFLD, in psychiatric patients and the high rate of off-label prescribing of AA medications to non-psychotic patients including children and older adults, the likely impact of AA prescribing on population incidence of NAFLD is significant. Furthermore, as there is currently no US Food and Drug Administration (FDA)-approved therapeutic for NAFLD, increasing our understanding of the organ-specific as well as the multi-system mechanisms underlying disease pathology as well as elucidating pharmacological mechanisms associated with AA-induced NAFLD can aid in identifying novel therapeutic targets for this rapidly growing, global patient population.

4. Materials and Methods

4.1. Animals and Experimental Design

This study utilized 8-week-old male C57BL/6 J mice (The Jackson Laboratory, Bar Harbor, ME, USA) fed a standard chow diet for the duration of the study (28 days). We deliberately chose to maintain these animals on a standard, low-fat mouse chow diet in order to test effects of AA treatment per se, as high fat diets are known to induce insulin resistance and metabolic disease in this model. Animals were randomly assigned to 1 of 3 treatments: (a) vehicle (VEH) control, (b) risperidone (RIS), or (c) olanzapine (OLAN). At the culmination of the study, animals were euthanized by intraperitoneal (IP) injection of 2.5% avertin (in PBS) followed by cardiac perfusion. Livers were extracted following perfusion and were either snap frozen in liquid nitrogen for RNA and protein analysis or fixed in 4% paraformaldehyde and paraffin embedded. All animal experiments were carried out in compliance with the federal animal welfare laws and policies using protocol # 052918-003 (approved for operation from June 2018 until June 2021) approved by the Institutional Animal Care and Use Committee (IACUC) at the University of New England, Biddeford, ME, USA, following the National Institutes of Health guide for the care and use of laboratory animals [74].

4.2. Drug Formulation and Dosing Strategy

The dosing strategy employed was designed to achieve peak plasma drug concentrations that fall within the range of plasma drug exposures observed clinically in patients [24,75]. Briefly, animals were administered drug or vehicle via once daily oral (PO) gavage for 4 weeks (28 days). Doses included VEH (0.1 % acetic acid), RIS (Sigma-Aldrich, St. Louis, MO, USA; 1.0 mg/kg), or OLAN (Sigma-Aldrich, St. Louis, MO, USA; 5 mg/kg). Our dosing strategy was based on pharmacokinetic studies previously conducted in our laboratory using male and female mice [54,61].

4.3. Histopathology

Sections of each liver lobe were collected during necropsy and post-fixed in 4% paraformaldehyde overnight. Tissues were embedded in paraffin and 5 mm sections were taken for staining using hematoxylin and eosin (H&E) in order to qualitatively assess changes in pathology. Sections from each lobe were examined via brightfield microscopy at 20× magnification using a Keyence BZ-X710 inverted wide field digital microscope. Images were collected using BZ-X analyzer software. Histopathologic lesions were evaluated using a modification of the scoring criteria of Jensen et al. [76]. Specific modifications included (1) evaluation of tissues for Mallory–Denk bodies (present = 1; absent = 0), (2) scoring of “steatosis” parameters as microsteatosis, and (3) evaluation of macrosteatosis (i.e., intercellular steatosis) (lipid < 5% fieldview = 0; lipid 5%–33% fieldview = 1; lipid 33%–66% of fieldview = 2; lipid > 66% of fieldview = 3). Lesions scores were averaged across lobes for each animal and reported as means per treatment group.

4.4. Proteomic Methodology

Livers from VEH-, RIS-, and OLAN-treated mice were homogenized in HEPES-Triton-Sodium-Glycerol (HTNG) lysis buffer (20 mM hydroxyethyl-piperazineethanesulfonic acid (HEPES), 150 mM NaCl, 1.5 mM MgCl2, 10% glycerol, 1% Triton-X 100, 1 mM ethylenediamine tretraacetic acid (EDTA), protease-inhibitor cocktail (Calbiochem)). Protein concentration of the supernatant was determined by Bicinchoninic acid assay (BCA, Pierce). A total of 40 μg of protein was used from each sample. Tryptic digests of protein samples were performed using the ProteoExtract digestion kit (Calbiochem). Tryptic peptides were then separated on an Ultimate RSLC system 3000 (ThermoFisher/Dionex) nanoscale liquid chromatograph and infused onto a 5600 TripleTOF mass spectrometer (Sciex). Sequential window acquisition of all theoretical spectra (SWATH) was used to profile all proteins in each sample using a data-independent acquisition method. A mouse-specific ion library comprising 4091 proteins was constructed using ProteinPilot software (Sciex). For identification of peptides, we retrieved multiple fragment ion chromatograms from the spectral library for each peptide of interest. These spectra were compared with the extracted fragment ion traces for the corresponding isolation window to identify the transitions that best identify the target peptide. SWATH analysis was performed using PeakView software, and MarkerView software was utilized for principal component analysis and t-test comparisons [77]. Detailed proteomic analysis methods are available at PeptideAtlas (identifier: PASS01349). The proteomic data were analyzed in the following pairwise comparisons: liver RIS to liver VEH and liver OLAN to liver VEH. An adjusted p-value was calculated for each protein in each pairwise dataset using the False Data Recovery (FDR) method. Two significant protein lists were then made for each pairwise dataset using the raw and adjusted p-values with a threshold of 0.05. Proteins associated or phenotypes showing altered expression during AA treatment relative to VEH-treated controls were binned and tabulated for this study.

4.5. Proteomic Analysis

The proteomic data were analyzed comparing RIS or OLAN to vehicle control after log transformation to achieve normality. A differentially expressed (DE) protein list with a rounded p-value threshold of 0.05 was established for pathway analysis. Pathway analysis was performed using the Kyoto Encyclopedia of Genes and Genomes (KEGG) database, releases 95.0, 95.1, and 95.2 (https://www.genome.jp/kegg) [78,79,80]. DE proteins were queried on the KEGG, and all pathway maps and corresponding object identifier numbers attributed to each protein were recorded. DE proteins without KEGG pathway maps were excluded from subsequent evaluations. The role of each DE protein on each KEGG pathway map was determined on the basis of the relationships with other proteins and molecules with which it interacted, according to the KEGG notations key (https://www.genome.jp/kegg/document/help_pathway.html). Changes in predicted physiological function noted in each KEGG pathway map, upstream and/or downstream in the signaling pathways relative to the DE protein, and the effects on those predicted functions caused by the upregulation or downregulation of the DE proteins were recorded. The specific predicted functional changes caused by the change in expression of each DE protein were recorded separately for each combination DE protein/KEGG pathway map. Subsequent to all DE proteins being analyzed, we tabulated the frequency of each predicted physiological function’s increase (“up”) or decrease (“down”) for each drug. For the evaluation of larger systemic changes, we sorted the predicted physiological functions into classifications (categories) on the basis of similar biological functions (e.g., metabolic functions).

4.6. RNA Analysis

Liver tissue was homogenized in liquid nitrogen, total RNA was isolated by TriReagent/chloroform extraction, and RNA quality was assessed by Nanodrop (Thermo Fisher Scientific, Waltham, MA, USA). Total RNA was then analyzed using RT2 Profiler PCR Array (Mouse PPAR Targets #330231; Qiagen).

4.7. Statistical Analysis

All data are expressed as the mean ± standard error of the mean (SEM), unless otherwise specified. All statistical analysis was performed using Student’s t-test or ANOVA with post hoc analyses, as indicated unless specified. All statistical analysis was performed using Prism 8 statistical software (GraphPad Software, Inc., La Jolla, CA, USA).

4.8. Data Statement

The proteomic datasets analyzed for this study can be found in the PeptideAtlas (identifier: PASS01349) http://www.peptideatlas.org/PASS/PASS01349. The raw data supporting the conclusions of this manuscript will be made available by the authors, without undue reservation, to any qualified researcher.

Supplementary Materials

The following are available online at https://www.mdpi.com/1422-0067/21/24/9362/s1. Five supplementary data files are provided in support of this manuscript.

Author Contributions

Conceptualization, K.L.H., M.M. and M.B.; methodology, B.R., M.B., C.B., C.P.H.V., E.B., M.M. and K.L.H.; formal analysis, M.B., C.B., B.R., M.M., K.L.H. and C.P.H.V.; resources, K.L.H. and C.P.H.V.; data curation, M.M., B.R., M.B. and K.L.H.; writing—original draft preparation, K.L.H., M.M. and B.R.; writing—review and editing, B.R., E.B., C.P.H.V., M.M. and K.L.H.; visualization, B.R., M.M., C.B. and K.L.H.; supervision, K.L.H. and M.M.; project administration, K.L.H.; funding acquisition, K.L.H. All authors have read and agreed to the published version of the manuscript.

Funding

This work was supported by the National Institutes of Health NIDDK award number DK095143 to KLH. The Proteomics and Lipidomics Core at Maine Medical Center Research Institute receives funding support from National Institutes of Health grants P30GM103392 (Robert Friesel/Donald St. Germain), P20GM121301 (Lucy Liaw), and U54GM115516 (Clifford Rosen).

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Chesney, E.; Goodwin, G.M.; Fazel, S. Risks of All-Cause and Suicide Mortality in Mental Disorders: A Meta-Review. World Psychiatry 2014, 13, 153–160. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  2. Vancampfort, D.; Stubbs, B.; Mitchell, A.J.; De Hert, M.; Wampers, M.; Ward, P.B.; Rosenbaum, S.; Correll, C.U. Risk of Metabolic Syndrome and Its Components in People with Schizophrenia and Related Psychotic Disorders, Bipolar Disorder and Major Depressive Disorder: A Systematic Review and Meta-Analysis. World Psychiatry 2015, 14, 339–347. [Google Scholar] [CrossRef] [PubMed]
  3. Soto-Angona, Ó.; Anmella, G.; Valdés-Florido, M.J.; De Uribe-Viloria, N.; Carvalho, A.F.; Penninx, B.W.J.H.; Berk, M. Non-Alcoholic Fatty Liver Disease (NAFLD) as a Neglected Metabolic Companion of Psychiatric Disorders: Common Pathways and Future Approaches. BMC Med. 2020, 18, 261. [Google Scholar] [CrossRef] [PubMed]
  4. Mikolasevic, I.; Milic, S.; Turk Wensveen, T.; Grgic, I.; Jakopcic, I.; Stimac, D.; Wensveen, F.; Orlic, L. Nonalcoholic Fatty Liver Disease—A Multisystem Disease? World J. Gastroenterol. 2016, 22, 9488–9505. [Google Scholar] [CrossRef] [PubMed]
  5. Engin, A. Non-Alcoholic Fatty Liver Disease. Adv. Exp. Med. Biol. 2017, 960, 443–467. [Google Scholar] [CrossRef] [PubMed]
  6. Mouzaki, M.; Yodoshi, T.; Arce-Clachar, A.C.; Bramlage, K.; Fei, L.; Ley, S.L.; Xanthakos, S.A. Psychotropic Medications Are Associated with Increased Liver Disease Severity in Pediatric Nonalcoholic Fatty Liver Disease. J. Pediatric Gastroenterol. Nutr. 2019, 69, 339–343. [Google Scholar] [CrossRef] [PubMed]
  7. Amir, M.; Yu, M.; He, P.; Srinivasan, S. Hepatic Autonomic Nervous System and Neurotrophic Factors Regulate the Pathogenesis and Progression of Non-Alcoholic Fatty Liver Disease. Front. Med. 2020, 7, 62. [Google Scholar] [CrossRef] [Green Version]
  8. Xu, H.; Zhuang, X. Atypical Antipsychotics-Induced Metabolic Syndrome and Nonalcoholic Fatty Liver Disease: A Critical Review. Neuropsychiatr. Dis. Treat 2019, 15, 2087–2099. [Google Scholar] [CrossRef] [Green Version]
  9. Zhu, Y.; Zhang, C.; Xu, F.; Zhao, M.; Bergquist, J.; Yang, C.; Liu, X.; Tan, Y.; Wang, X.; Li, S.; et al. System Biology Analysis Reveals the Role of Voltage-Dependent Anion Channel in Mitochondrial Dysfunction during Non-Alcoholic Fatty Liver Disease Progression into Hepatocellular Carcinoma. Cancer Sci. 2020, 111, 4288–4302. [Google Scholar] [CrossRef]
  10. Cobbina, E.; Akhlaghi, F. Non-Alcoholic Fatty Liver Disease (NAFLD)—Pathogenesis, Classification, and Effect on Drug Metabolizing Enzymes and Transporters. Drug Metab. Rev. 2017, 49, 197–211. [Google Scholar] [CrossRef]
  11. Anty, R.; Gual, P. Pathogenesis of Non-Alcoholic Fatty Liver Disease. Presse Med. 2019, 48, 1468–1483. [Google Scholar] [CrossRef] [PubMed]
  12. Estes, C.; Razavi, H.; Loomba, R.; Younossi, Z.; Sanyal, A.J. Modeling the Epidemic of Nonalcoholic Fatty Liver Disease Demonstrates an Exponential Increase in Burden of Disease. Hepatology 2018, 67, 123–133. [Google Scholar] [CrossRef] [PubMed]
  13. Benedict, M.; Zhang, X. Non-Alcoholic Fatty Liver Disease: An Expanded Review. World J. Hepatol. 2017, 9, 715–732. [Google Scholar] [CrossRef]
  14. Harris, L.W.; Guest, P.C.; Wayland, M.T.; Umrania, Y.; Krishnamurthy, D.; Rahmoune, H.; Bahn, S. Schizophrenia: Metabolic Aspects of Aetiology, Diagnosis and Future Treatment Strategies. Psychoneuroendocrinology 2013, 38, 752–766. [Google Scholar] [CrossRef] [PubMed]
  15. Mitchell, A.J.; Vancampfort, D.; Sweers, K.; van Winkel, R.; Yu, W.; De Hert, M. Prevalence of Metabolic Syndrome and Metabolic Abnormalities in Schizophrenia and Related Disorders—A Systematic Review and Meta-Analysis. Schizophr. Bull. 2013, 39, 306–318. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  16. Papanastasiou, E. The Prevalence and Mechanisms of Metabolic Syndrome in Schizophrenia: A Review. Ther. Adv. Psychopharmacol. 2013, 3, 33–51. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  17. Freyberg, Z.; Aslanoglou, D.; Shah, R.; Ballon, J.S. Intrinsic and Antipsychotic Drug-Induced Metabolic Dysfunction in Schizophrenia. Front. Neurosci. 2017, 11, 432. [Google Scholar] [CrossRef]
  18. Pillinger, T.; McCutcheon, R.A.; Vano, L.; Mizuno, Y.; Arumuham, A.; Hindley, G.; Beck, K.; Natesan, S.; Efthimiou, O.; Cipriani, A.; et al. Comparative Effects of 18 Antipsychotics on Metabolic Function in Patients with Schizophrenia, Predictors of Metabolic Dysregulation, and Association with Psychopathology: A Systematic Review and Network Meta-Analysis. Lancet Psychiatry 2020, 7, 64–77. [Google Scholar] [CrossRef]
  19. Siafis, S.; Tzachanis, D.; Samara, M.; Papazisis, G. Antipsychotic Drugs: From Receptor-Binding Profiles to Metabolic Side Effects. Curr. Neuropharmacol. 2018, 16, 1210–1223. [Google Scholar] [CrossRef]
  20. American Diabetes Association; American Psychiatric Association; American Association of Clinical Endocrinologists; North American Association for the Study of Obesity. Consensus Development Conference on Antipsychotic Drugs and Obesity and Diabetes. Diabetes Care 2004, 27, 596–601. [Google Scholar] [CrossRef] [Green Version]
  21. Beauchemin, M.; Geguchadze, R.; Guntur, A.R.; Nevola, K.; Le, P.T.; Barlow, D.; Rue, M.; Vary, C.P.H.; Lary, C.W.; Motyl, K.J.; et al. Exploring Mechanisms of Increased Cardiovascular Disease Risk with Antipsychotic Medications: Risperidone Alters the Cardiac Proteomic Signature in Mice. Pharmacol. Res. 2020, 152, 104589. [Google Scholar] [CrossRef] [PubMed]
  22. Velázquez, K.T.; Enos, R.T.; Bader, J.E.; Sougiannis, A.T.; Carson, M.S.; Chatzistamou, I.; Carson, J.A.; Nagarkatti, P.S.; Nagarkatti, M.; Murphy, E.A. Prolonged High-Fat-Diet Feeding Promotes Non-Alcoholic Fatty Liver Disease and Alters Gut Microbiota in Mice. World J. Hepatol. 2019, 11, 619–637. [Google Scholar] [CrossRef] [PubMed]
  23. Recena Aydos, L.; Aparecida do Amaral, L.; Serafim de Souza, R.; Jacobowski, A.C.; Freitas dos Santos, E.; Rodrigues Macedo, M.L. Nonalcoholic Fatty Liver Disease Induced by High-Fat Diet in C57bl/6 Models. Nutrients 2019, 11, 3067. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  24. Mauri, M.C.; Paletta, S.; Maffini, M.; Colasanti, A.; Dragogna, F.; Di Pace, C.; Altamura, A.C. Clinical Pharmacology of Atypical Antipsychotics: An Update. EXCLI J. 2014, 13, 1163–1191. [Google Scholar]
  25. Ader, M.; Kim, S.P.; Catalano, K.J.; Ionut, V.; Hucking, K.; Richey, J.M.; Kabir, M.; Bergman, R.N. Metabolic Dysregulation with Atypical Antipsychotics Occurs in the Absence of Underlying Disease A Placebo-Controlled Study of Olanzapine and Risperidone in Dogs. Diabetes 2005, 54, 862–871. [Google Scholar] [CrossRef] [Green Version]
  26. Houseknecht, K.L.; Robertson, A.S.; Zavadoski, W.; Gibbs, E.M.; Johnson, D.E.; Rollema, H. Acute Effects of Atypical Antipsychotics on Whole-Body Insulin Resistance in Rats: Implications for Adverse Metabolic Effects. Neuropsychopharmacology 2007, 32, 289–297. [Google Scholar] [CrossRef] [Green Version]
  27. Martins, P.J.F.; Haas, M.; Obici, S. Central Nervous System Delivery of the Antipsychotic Olanzapine Induces Hepatic Insulin Resistance. Diabetes 2010, 59, 2418–2425. [Google Scholar] [CrossRef] [Green Version]
  28. DeFronzo, R.A. Insulin Resistance, Lipotoxicity, Type 2 Diabetes and Atherosclerosis: The Missing Links. The Claude Bernard Lecture 2009. Diabetologia 2010, 53, 1270–1287. [Google Scholar] [CrossRef] [Green Version]
  29. Schaffer, J.E. Lipotoxicity: Many Roads to Cell Dysfunction and Cell Death: Introduction to a Thematic Review Series. J. Lipid Res. 2016, 57, 1327–1328. [Google Scholar] [CrossRef] [Green Version]
  30. Czech, M.P. Insulin Action and Resistance in Obesity and Type 2 Diabetes. Nat. Med. 2017, 23, 804–814. [Google Scholar] [CrossRef]
  31. Hong, F.; Pan, S.; Guo, Y.; Xu, P.; Zhai, Y. PPARs as Nuclear Receptors for Nutrient and Energy Metabolism. Molecules 2019, 24, 2545. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  32. Lamas Bervejillo, M.; Ferreira, A.M. Understanding Peroxisome Proliferator-Activated Receptors: From the Structure to the Regulatory Actions on Metabolism. Adv. Exp. Med. Biol. 2019, 1127, 39–57. [Google Scholar] [CrossRef] [PubMed]
  33. Astapova, O.; Leff, T. Adiponectin and PPARγ: Cooperative and Interdependent Actions of Two Key Regulators of Metabolism. Vitam. Horm. 2012, 90, 143–162. [Google Scholar] [CrossRef] [PubMed]
  34. Iwabu, M.; Okada-Iwabu, M.; Yamauchi, T.; Kadowaki, T. Adiponectin/AdipoR Research and Its Implications for Lifestyle-Related Diseases. Front. Cardiovasc. Med. 2019, 6, 116. [Google Scholar] [CrossRef] [Green Version]
  35. Straub, L.G.; Scherer, P.E. Metabolic Messengers: Adiponectin. Nat. Metab. 2019, 1, 334–339. [Google Scholar] [CrossRef]
  36. Albaugh, V.L.; Vary, T.C.; Ilkayeva, O.; Wenner, B.R.; Maresca, K.P.; Joyal, J.L.; Breazeale, S.; Elich, T.D.; Lang, C.H.; Lynch, C.J. Atypical Antipsychotics Rapidly and Inappropriately Switch Peripheral Fuel Utilization to Lipids, Impairing Metabolic Flexibility in Rodents. Schizophr. Bull. 2012, 38, 153–166. [Google Scholar] [CrossRef] [Green Version]
  37. Qi, Y.; Xu, R. Roles of PLODs in Collagen Synthesis and Cancer Progression. Front. Cell. Dev. Biol. 2018, 6, 66. [Google Scholar] [CrossRef]
  38. Shoelson, S.E.; Lee, J.; Goldfine, A.B. Inflammation and Insulin Resistance. J. Clin. Invest. 2006, 116, 1793–1801. [Google Scholar] [CrossRef]
  39. Chen, Z.; Yu, R.; Xiong, Y.; Du, F.; Zhu, S. A Vicious Circle between Insulin Resistance and Inflammation in Nonalcoholic Fatty Liver Disease. Lipids Health Dis. 2017, 16, 203. [Google Scholar] [CrossRef] [Green Version]
  40. Zatterale, F.; Longo, M.; Naderi, J.; Raciti, G.A.; Desiderio, A.; Miele, C.; Beguinot, F. Chronic Adipose Tissue Inflammation Linking Obesity to Insulin Resistance and Type 2 Diabetes. Front. Physiol. 2020, 10, 1607. [Google Scholar] [CrossRef]
  41. Katsarou, A.; Moustakas, I.I.; Pyrina, I.; Lembessis, P.; Koutsilieris, M.; Chatzigeorgiou, A. Metabolic Inflammation as an Instigator of Fibrosis during Non-Alcoholic Fatty Liver Disease. World J. Gastroenterol. 2020, 26, 1993–2011. [Google Scholar] [CrossRef] [PubMed]
  42. Tanwar, S.; Rhodes, F.; Srivastava, A.; Trembling, P.M.; Rosenberg, W.M. Inflammation and Fibrosis in Chronic Liver Diseases Including Non-Alcoholic Fatty Liver Disease and Hepatitis C. World J. Gastroenterol. 2020, 26, 109–133. [Google Scholar] [CrossRef] [PubMed]
  43. Arrese, M.; Cabrera, D.; Kalergis, A.M.; Feldstein, A.E. Innate Immunity and Inflammation in NAFLD/NASH. Dig Dis Sci 2016, 61, 1294–1303. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  44. Dawood, R.M.; El-Meguid, M.A.; Salum, G.M.; El Awady, M.K. Key Players of Hepatic Fibrosis. J Interferon Cytokine Res 2020, 40, 472–489. [Google Scholar] [CrossRef] [PubMed]
  45. Franceschetti, L.; Bonomini, F.; Rodella, L.F.; Rezzani, R. Critical Role of NFκB in the Pathogenesis of Non-Alcoholic Fatty Liver Disease: A Widespread Key Regulator. Curr. Mol. Med. 2020. [Google Scholar] [CrossRef] [PubMed]
  46. Sutti, S.; Albano, E. Adaptive Immunity: An Emerging Player in the Progression of NAFLD. Nat. Rev. Gastroenterol. Hepatol. 2020, 17, 81–92. [Google Scholar] [CrossRef] [PubMed]
  47. Stine, J.G.; Intagliata, N.; Northup, P.G.; Caldwell, S.H. Nonalcoholic Fatty Liver Disease, Portal Vein Thrombosis and Coagulation: More Questions Than Answers? Transplantation 2017, 101, e281–e282. [Google Scholar] [CrossRef]
  48. Mizuno, K.; Ueno, Y. Autonomic Nervous System and the Liver. Hepatol. Res. 2017, 47, 160–165. [Google Scholar] [CrossRef] [Green Version]
  49. Leung, J.Y.T.; Barr, A.M.; Procyshyn, R.M.; Honer, W.G.; Pang, C.C.Y. Cardiovascular Side-Effects of Antipsychotic Drugs: The Role of the Autonomic Nervous System. Pharmacol. Ther. 2012, 135, 113–122. [Google Scholar] [CrossRef]
  50. Scigliano, G.; Ronchetti, G. Antipsychotic-Induced Metabolic and Cardiovascular Side Effects in Schizophrenia: A Novel Mechanistic Hypothesis. CNS Drugs 2013, 27, 249–257. [Google Scholar] [CrossRef] [Green Version]
  51. Motyl, K.J.; DeMambro, V.E.; Barlow, D.; Olshan, D.; Nagano, K.; Baron, R.; Rosen, C.J.; Houseknecht, K.L. Propranolol Attenuates Risperidone-Induced Trabecular Bone Loss in Female Mice. Endocrinology 2015, 156, 2374–2383. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  52. Ehrlich, L.; Scrushy, M.; Meng, F.; Lairmore, T.C.; Alpini, G.; Glaser, S. Biliary Epithelium: A Neuroendocrine Compartment in Cholestatic Liver Disease. Clin. Res. Hepatol. Gastroenterol. 2018, 42, 296–305. [Google Scholar] [CrossRef] [PubMed]
  53. Blencowe, M.; Karunanayake, T.; Wier, J.; Hsu, N.; Yang, X. Network Modeling Approaches and Applications to Unravelling Non-Alcoholic Fatty Liver Disease. Genes 2019, 10, 966. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  54. May, M.; Beauchemin, M.; Vary, C.; Barlow, D.; Houseknecht, K.L. The Antipsychotic Medication, Risperidone, Causes Global Immunosuppression in Healthy Mice. PLoS ONE 2019, 14, e0218937. [Google Scholar] [CrossRef] [PubMed]
  55. May, M.; Slitzky, M.; Rostama, B.; Barlow, D.; Houseknecht, K.L. Antipsychotic-Induced Immune Dysfunction: A Consideration for COVID-19 risk. Brain Behav. Immun. Health 2020, 6, 100097. [Google Scholar] [CrossRef] [PubMed]
  56. Drzyzga, L.; Obuchowicz, E.; Marcinowska, A.; Herman, Z.S. Cytokines in Schizophrenia and the Effects of Antipsychotic Drugs. Brain Behav. Immun. 2006, 20, 532–545. [Google Scholar] [CrossRef]
  57. Correll, C.U. Cardiometabolic Risk of Second-Generation Antipsychotic Medications during First-Time Use in Children and Adolescents. JAMA 2009, 302, 1765–1773. [Google Scholar] [CrossRef] [Green Version]
  58. Kowalchuk, C.; Castellani, L.N.; Chintoh, A.; Remington, G.; Giacca, A.; Hahn, M.K. Antipsychotics and Glucose Metabolism: How Brain and Body Collide. Am. J. Physiol. Endocrinol. Metab. 2018, 316, E1–E15. [Google Scholar] [CrossRef] [Green Version]
  59. Kim, H.J.; Wilson, C.; Deusen, T.V.; Millard, H.; Qayyum, Z.; Parke, S. Metabolic Syndrome in Child and Adolescent Psychiatry. Psychiatr. Ann. 2020, 50, 326–333. [Google Scholar] [CrossRef]
  60. Müller, N.; Weidinger, E.; Leitner, B.; Schwarz, M.J. The Role of Inflammation in Schizophrenia. Front. Neurosci. 2015, 9, 372. [Google Scholar] [CrossRef] [Green Version]
  61. Gaughran, F. Immunity and Schizophrenia: Autoimmunity, Cytokines, and Immune Responses. Int. Rev. Neurobiol. 2002, 52, 275–302. [Google Scholar] [CrossRef] [PubMed]
  62. Schuld, A.; Hinze-Selch, D.; Pollmächer, T. Cytokine Network in Patients with Schizophrenia and Its Significance for the Pathophysiology of the Illness. Nervenarzt 2004, 75, 215–226. [Google Scholar] [CrossRef] [PubMed]
  63. Calarge, C.A.; Ivins, S.D.; Motyl, K.J.; Shibli-Rahhal, A.A.; Bliziotes, M.M.; Schlechte, J.A. Possible Mechanisms for the Skeletal Effects of Antipsychotics in Children and Adolescents. Ther. Adv. Psychopharmacol. 2013, 3, 278–293. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  64. Motyl, K.J.; Beauchemin, M.; Barlow, D.; Le, P.T.; Nagano, K.; Treyball, A.; Contractor, A.; Baron, R.; Rosen, C.J.; Houseknecht, K.L. A Novel Role for Dopamine Signaling in the Pathogenesis of Bone Loss from the Atypical Antipsychotic Drug Risperidone in Female Mice. Bone 2017, 103, 168–176. [Google Scholar] [CrossRef]
  65. Caldez, M.J.; Bjorklund, M.; Kaldis, P. Cell Cycle Regulation in NAFLD: When Imbalanced Metabolism Limits Cell Division. Hepatol. Int. 2020, 14, 463–474. [Google Scholar] [CrossRef]
  66. Wojcikowski, J.; Anna Daniel, W. The Role of the Nervous System in the Regulation of Liver Cytochrome P450. Curr. Drug Metab. 2011, 12, 124–138. [Google Scholar] [CrossRef]
  67. Jensen, K.J.; Alpini, G.; Glaser, S. Hepatic Nervous System and Neurobiology of the Liver. Compr. Physiol. 2013, 3, 655–665. [Google Scholar] [CrossRef] [Green Version]
  68. Bär, K.-J. Cardiac Autonomic Dysfunction in Patients with Schizophrenia and Their Healthy Relatives—A Small Review. Front. Neurol. 2015, 6, 139. [Google Scholar] [CrossRef] [Green Version]
  69. Esler, M.; Rumantir, M.; Wiesner, G.; Kaye, D.; Hastings, J.; Lambert, G. Sympathetic Nervous System and Insulin Resistance: From Obesity to Diabetes. Am. J. Hypertens. 2001, 14, 304S–309S. [Google Scholar] [CrossRef] [Green Version]
  70. Oliveira, C.; Silveira, E.A.; Rosa, L.; Santos, A.; Rodrigues, A.P.; Mendonça, C.; Silva, L.; Gentil, P.; Rebelo, A.C. Risk Factors Associated with Cardiac Autonomic Modulation in Obese Individuals. J. Obes. 2020, 2020. [Google Scholar] [CrossRef]
  71. Grassi, G.; Biffi, A.; Dell’Oro, R.; Quarti Trevano, F.; Seravalle, G.; Corrao, G.; Perseghin, G.; Mancia, G. Sympathetic Neural Abnormalities in Type 1 and Type 2 Diabetes: A Systematic Review and Meta-Analysis. J. Hypertens. 2020, 38, 1436–1442. [Google Scholar] [CrossRef] [PubMed]
  72. Savoy, Y.E.; Ashton, M.A.; Miller, M.W.; Nedza, F.M.; Spracklin, D.K.; Hawthorn, M.H.; Rollema, H.; Matos, F.F.; Hajos-Korcsok, E. Differential Effects of Various Typical and Atypical Antipsychotics on Plasma Glucose and Insulin Levels in the Mouse: Evidence for the Involvement of Sympathetic Regulation. Schizophr. Bull. 2010, 36, 410–418. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  73. Hattori, S.; Kishida, I.; Suda, A.; Miyauchi, M.; Shiraishi, Y.; Fujibayashi, M.; Tsujita, N.; Ishii, C.; Ishii, N.; Moritani, T.; et al. Effects of Four Atypical Antipsychotics on Autonomic Nervous System Activity in Schizophrenia. Schizophr. Res. 2018, 193, 134–138. [Google Scholar] [CrossRef] [PubMed]
  74. National Research Council of the National Academies; Institute for Laboratory Animals Research; Committee for the Update of the Guide for the Care and Use of Laboratory Animals. Guide for the Care and Use of Laboratory Animals: Eighth Edition; National Academies Press: Washington, DC, USA, 2010; ISBN 978-0-309-18663-6. [Google Scholar]
  75. Aman, M.G.; Vinks, A.A.; Remmerie, B.; Mannaert, E.; Ramadan, Y.; Masty, J.; Lindsay, R.L.; Malone, K. Plasma Pharmacokinetic Characteristics of Risperidone and Their Relationship to Saliva Concentrations in Children with Psychiatric or Neurodevelopmental Disorders. Clin. Ther. 2007, 29, 1476–1486. [Google Scholar] [CrossRef] [Green Version]
  76. Jensen, V.S.; Tveden-Nyborg, P.; Zacho-Rasmussen, C.; Quaade, M.L.; Ipsen, D.H.; Hvid, H.; Fledelius, C.; Wulff, E.M.; Lykkesfeldt, J. Variation in Diagnostic NAFLD/NASH Read-Outs in Paired Liver Samples from Rodent Models. J. Pharmacol. Toxicol. Methods 2020, 101, 106651. [Google Scholar] [CrossRef]
  77. Gillet, L.; Navarro, P.; Tate, S.; Rost, H.; Selevsek, N.; Reiter, L.; Bonner, R.; Aebersold, R. Targeted Data Extraction of the MS/MS Spectra Generated by Data-Independent Acquisition: A New Concept for Consistent and Accurate Proteome Analysis. Mol. Cell. Proteomics 2012, 11, O111.016717. [Google Scholar] [CrossRef] [Green Version]
  78. Kanehisa, M.; Goto, S. KEGG: Kyoto Encyclopedia of Genes and Genomes. Nucleic Acids Res. 2000, 28, 27–30. [Google Scholar] [CrossRef]
  79. Kanehisa, M.; Sato, Y.; Furumichi, M.; Morishima, K.; Tanabe, M. New Approach for Understanding Genome Variations in KEGG. Nucleic Acids Res. 2019, 47, D590–D595. [Google Scholar] [CrossRef] [Green Version]
  80. Kanehisa, M. Toward Understanding the Origin and Evolution of Cellular Organisms. Protein Sci. 2019, 28, 1947–1951. [Google Scholar] [CrossRef]
Figure 1. Atypical antipsychotic (AA) treatment led to histopathologic changes consistent with early non-alcoholic fatty liver disease (NAFLD). Hematoxylin and eosin staining showed several histopathologic changes at 20× magnification following AA treatment relative to vehicle (VEH)-treated controls (A). These included micro- and macrosteatosis, ballooned hepatocytes (green arrows), Mallory–Denk bodies (yellow arrows), and small areas of early fibrosis (boxed). Two representative field views are shown. Statistical significance (*, **) in lesion scores between AA-treated mice and VEH-treated mice was found for microsteatosis, Mallory–Denk body formation, and formation of inflammatory foci (B).
Figure 1. Atypical antipsychotic (AA) treatment led to histopathologic changes consistent with early non-alcoholic fatty liver disease (NAFLD). Hematoxylin and eosin staining showed several histopathologic changes at 20× magnification following AA treatment relative to vehicle (VEH)-treated controls (A). These included micro- and macrosteatosis, ballooned hepatocytes (green arrows), Mallory–Denk bodies (yellow arrows), and small areas of early fibrosis (boxed). Two representative field views are shown. Statistical significance (*, **) in lesion scores between AA-treated mice and VEH-treated mice was found for microsteatosis, Mallory–Denk body formation, and formation of inflammatory foci (B).
Ijms 21 09362 g001
Figure 2. Proteomic analysis of hepatic proteins differentially expressed (DE) by both risperidone and olanzapine. (A) Venn diagrams showing shared proteins that are significantly upregulated (top) or downregulated (bottom) in risperidone (RIS)- or olanzapine (OLAN)-treated mice (p < 0.05). (B) Of the 20 shared proteins, 9 were found to have KEGG (Kyoto Encyclopedia of Genes and Genomes) pathway maps from which predicted functional changes could be determined. Statistical significance (p < 0.5) of differentially expressed (DE) proteins relative to vehicle (VEH). RIS n = 5; OLAN n = 4; VEH n = 4. Abbreviations: NAFLD, non-alcoholic fatty liver disease; GI, gastrointestinal.
Figure 2. Proteomic analysis of hepatic proteins differentially expressed (DE) by both risperidone and olanzapine. (A) Venn diagrams showing shared proteins that are significantly upregulated (top) or downregulated (bottom) in risperidone (RIS)- or olanzapine (OLAN)-treated mice (p < 0.05). (B) Of the 20 shared proteins, 9 were found to have KEGG (Kyoto Encyclopedia of Genes and Genomes) pathway maps from which predicted functional changes could be determined. Statistical significance (p < 0.5) of differentially expressed (DE) proteins relative to vehicle (VEH). RIS n = 5; OLAN n = 4; VEH n = 4. Abbreviations: NAFLD, non-alcoholic fatty liver disease; GI, gastrointestinal.
Ijms 21 09362 g002
Figure 3. Dysregulation of pathways leading to NAFLD during risperidone and olanzapine treatment. DE proteins measured in this study are integral parts of several pathways in the KEGG database known to contribute to NAFLD. The number of proteins altered in RIS- or OLAN-treated mice (relative to VEH-treated controls) in each pathway group were tabulated (A). The intensity of red (λ = 100–255 nm) reflects the number of proteins per pathway significantly different from VEH-treated mice during RIS treatment (left panel) and OLAN treatment (right panel). The aggregate changes in functionality of pathways across proteins compared to VEH-treated mice for RIS and OLAN are depicted (B). RIS and OLAN dysregulated inflammation, fibrosis, insulin resistance, and lipogenesis at comparable levels, but vary markedly in their impact on pathways related to lipolysis, lipid transport, glycogen storage, cholesterol processing, and mitochondrial function.
Figure 3. Dysregulation of pathways leading to NAFLD during risperidone and olanzapine treatment. DE proteins measured in this study are integral parts of several pathways in the KEGG database known to contribute to NAFLD. The number of proteins altered in RIS- or OLAN-treated mice (relative to VEH-treated controls) in each pathway group were tabulated (A). The intensity of red (λ = 100–255 nm) reflects the number of proteins per pathway significantly different from VEH-treated mice during RIS treatment (left panel) and OLAN treatment (right panel). The aggregate changes in functionality of pathways across proteins compared to VEH-treated mice for RIS and OLAN are depicted (B). RIS and OLAN dysregulated inflammation, fibrosis, insulin resistance, and lipogenesis at comparable levels, but vary markedly in their impact on pathways related to lipolysis, lipid transport, glycogen storage, cholesterol processing, and mitochondrial function.
Ijms 21 09362 g003
Figure 4. Dysregulation of immune function pathways during risperidone and olanzapine treatment. DE proteins measured in this study are integral parts of several immune function pathways in the KEGG database. The number of proteins altered in RIS- or OLAN-treated mice (relative to VEH-treated controls) in each pathway group were tabulated. The intensity of red (λ = 100–255 nm) reflects the number of proteins per pathway significantly different from VEH-treated mice during RIS treatment (left panel) and OLAN treatment (right panel). Absence of any altered proteins in a particular pathway is depicted in black (λ = 0 nm).
Figure 4. Dysregulation of immune function pathways during risperidone and olanzapine treatment. DE proteins measured in this study are integral parts of several immune function pathways in the KEGG database. The number of proteins altered in RIS- or OLAN-treated mice (relative to VEH-treated controls) in each pathway group were tabulated. The intensity of red (λ = 100–255 nm) reflects the number of proteins per pathway significantly different from VEH-treated mice during RIS treatment (left panel) and OLAN treatment (right panel). Absence of any altered proteins in a particular pathway is depicted in black (λ = 0 nm).
Ijms 21 09362 g004
Figure 5. Model of antipsychotic regulation of NAFLD via sympathetic nervous system (SNS), metabolic, inflammatory, and endocrine factors. We propose that antipsychotic (AA) associated increase in sympathetic activity impacts the liver by causing alterations in the hepatic proteome and associated pathways leading to hepatic steatosis and NAFLD. A hematoxylin and eosin stained liver from OLAN-treated mouse at 20× magnification and pathway analysis heat maps (red-shading—RIS; green-shading—OLAN) are shown in support of this model.
Figure 5. Model of antipsychotic regulation of NAFLD via sympathetic nervous system (SNS), metabolic, inflammatory, and endocrine factors. We propose that antipsychotic (AA) associated increase in sympathetic activity impacts the liver by causing alterations in the hepatic proteome and associated pathways leading to hepatic steatosis and NAFLD. A hematoxylin and eosin stained liver from OLAN-treated mouse at 20× magnification and pathway analysis heat maps (red-shading—RIS; green-shading—OLAN) are shown in support of this model.
Ijms 21 09362 g005
Table 1. Selected results from the KEGG function categories for differentially expressed hepatic proteins from mice treated with risperidone.
Table 1. Selected results from the KEGG function categories for differentially expressed hepatic proteins from mice treated with risperidone.
Predicted Dysregulated Functions with RISUpDownAssociated Proteins
Predicted metabolic effectsEnergy production (all ATP-generating pathways combined)4162BCAT2, CS, AGXT, ABAT, NDUFA8, NDUFB8, NDUFS1, COX6B1, CYB5R3, ALDH7A1, TALDO1, ACAA2, ACLY, PRKCA, PDGFRB, GNAS, PGAM2, ENO3, LRP1, ALR1, AOX3, ADIPOQ, HIBADH, COMT, MTHFD1, ECI2, NUDT7, DGKZ, EPHA2, ITGA3, SLC27A1, BLVRB, SEC22B, SMARCC1, ASGR1, BTD
Amino acid production3237BCAT2, CS, ADIPOQ, LAP3, PGAM2, TALDO1, ENO3, ABAT, ACLY, PCBD2, ALR1, DGKZ, GNMT, AGXT, ALDH7A1, PRKCA, GNAS, FTCD, AMDHD1, ATIC, MTHFD1, EPHA2, PDGFRB, ITGA3, AOX3, COMT, ACAA2, HIBADH, TCN2
Glycolysis367ABAT, ACLY, ALDH7A1, AOX3, COX6B1, CS, EPHA2, GNAS, ITGA3, LAP3, NDUFS1, PDGFRB, PGAM2, PRKCA, TALDO1, ALR1, ENO3,
Lipolysis279ABAT, ACAA2, ACLY, ADIPOQ, ALDH7A1, BLVRB, COX6B1, CS, ECI2, GNAS, MECR, PDGFRB, PRKCA, SLC27A1, AOX3, NUDT7
Lipogenesis242ABAT, ACAA2, ADIPOQ, ALDH7A1, ALR1, AOX3, DGKZ, GNAS, MECR, PRKCA, SLC27A1, TCN2, SEC22B
TCA/OxPhos/ETC169ADIPOQ, BCAT2, NDUFA8, NDUFB8, PGAM2, ALDH7A1, ALR1, COX6B1, ENO3, NDUFS1
Nucleic acid synthesis713EPHA2, GNAS, ITGA3, PDGFRB, PRKCA, SMARCC1, TCN2, ALR1, AMDHD1, ATIC, FTCD, MTHFD1, NUDT7, TALDO1
Lysosomal function103ADIPOQ, CTSS, GSTA3, HGS, LAMP1, RAB11FIP5, SUMF1, RPN2, SEC22B, VAPA
Mitochondrial biogenesis110ADIPOQ, GNAS, PDGFRB, PRKCA, SLC27A1, SMARCC1
Predicted mitogenic effectsMigration/motility5521DGKZ, ENAH, EPHA2, GNAS, HSPB1, IL1RAP, ITGA3, MSN, NDUFA8, PDGFRB, PRKCA, ROCK2, ADIPOQ, AFDN, ASGR1, ELAVL1, INPP5D, MYH9, NDUFS1
Cell survival496CRYAB, DGKZ, EPHA2, GNAS, GSTA3, HSPB1, IL1RAP, ITGA3, NDUFA8, PDGFRB, PRKCA, RNF31, ADIPOQ, ASGR1, ELAVL1, INPP5D, NDUFS1
Cell proliferation466DGKZ, EPHA2, GNAS, HSPB1, IL1RAP, ITGA3, NDUFA8, NDUFB8, PDGFRB, PRKCA, ADIPOQ, ASGR1, ELAVL1, INPP5D, NDUFS1
Protein production1337EIF5, EPHA2, FAU, GNAS, HNRNPC, HSPB1, ITGA3, PDGFRB, PRKCA, TRA2B, AARS, ABAT, AGXT, ALDH7A1, ALR1, AMDHD1, ATIC, CRYAB, DDOST, EIF5B, FTCD, GANAB, HGS, HYOU1, NOP56, PABPC1, PUF60, RAB11FIP5, REX02, RPN2, RPS3A1, SEC22B, SNRPB, SNRPG, TALDO1, TARS
Cancer257DGKZ, GNAS, GSTA3, ITGA3, KLK3, PDGFRB, PRKCA, ROCK2, EPHA2, HPGD, SMARCC1
Predicted tissue-damaging effectsOxidation/ROS generation3410ABAT, AGXT, ALDH7A1, ALR1, AOX3, BLVRB, COMT, COX6B1, DDOST, EPHA2, GANAB, GNAS, HGS, YOU1, ITGA3, NDUFA8, NDUFB8, PDGFRB, PRKCA, RAB11FIP5, ROCK2, RPN2, ADIPOQ, COX6B1, GSTA3, LAP3, NDUFS1, SMARCC1,
Disrupted elimination of toxins/xenobiotics87GSTA3, LAP3, ALR1, AOX3, COX6B1, NUDT7, TALDO1
Protein misfolding91DDOST, GANAB, LAP3, RPN2, VCP
Table 2. Selected results from the KEGG function categories for differentially expressed hepatic proteins from mice treated with olanzapine.
Table 2. Selected results from the KEGG function categories for differentially expressed hepatic proteins from mice treated with olanzapine.
Predicted Dysregulated Functions of OLANUpDownAssociated Proteins
Predicted metabolic functionsEnergy production (all ATP-generating pathways combined)727G6PD, GLB1, NPC2, RPL34, AK2, AMY2B, ATP2B4, COX6B1, ENO3, GNAI, PNP, TNNC2, VAPA
Nucleic acid synthesis715G6PD, MTHFD1L, AK2, ATIC, ENO3, MTAP, PNP
Amino acid production710G6PD, GLB1, MTHFD1L, ATIC, ENO3, MTAP, PNP
Glycolysis610ATP2B4, G6PD, GLB1, P2RX4, AK2, AMY2B, ENO3
Lipogenesis125AGPAT3, G6PD, GNAI, HMGA2, HNRNPD, NPC2, ATP2B4, ENO3, VAPA
Lipolysis510AGPAT3, GNAI1, NPC2, AMY2B, COX6B1, CYCS, DECR2, ENO3, VAPA
TCA/OxPhos/ETC49G6PD, GLB1, RPL34, COX6B1, CYCS, ENO3, PNP
Glycogen production29ENO3, GNAI, ATP2B4, P2RX4
Predicted mitogenic effectsCell survival2114AGPAT3, ATP2B4, CYCS, FMOD, GNAI1, P2RX4, SGPL1, GLB1, GNAI1, INPP5D, LUM, PIP4P2, PNP, TOP2B
Cell proliferation815AGPAT3, FMOD, GLB1, GNAI1, SGPL1, ATP2B4, CTC1, ELAVL1, INPP5D, P2RX4, PIP4P2, SGPL1, TOP2B
Cancer135CYCS, G6PD, GNAI, HMGA2, PLOD3, LUM, TOP2B, TRIM32
Migration/motility69AGPAT3, CFL1, GLB1, SGPL1, ATP2B4, GNAI1, INPP5D, PIP4P2
Protein production76FAU, HNRNPC, NUP98, RPL34, RPL8, RPS15, RPS4Y1, ATP2B4, GNAI, RNPS1, TOP2B
Predicted immune effectsInflammation/cytokines239ACE, AMY2B, ATIC, CD68, CFL1, CNTN1, CYCS, G6PD, GLB1, GNAI1, MBL1, SGPL1, ELAVL1, INPP5D, MTHFD1L
Immune activation1013CD68, CFL1, DYNCH1H1, FMOD, G6PD, GLB1, GNAI, MBL1, AK2, CNTN1, CYCS, ELAVL1, INPP5D, PIP4P2, PNP
Cell microbe response147ACE, CFL1, CYCS, GNAI1, NCL, NUP98, SEPT2
Predicted GI/nutrient handlingGI function/nutrient absorption517GLB1, GNAI1, NPC2, AMY2B, AQP1, ATP2B4, TNNC2, VAPA
Bile processing510GLB1, GNAI1, NPC2, AQP1, ATP2B4, P2RX4, TNNC2, VAPA
Cholesterol processing55AMY2B, HNRNPD, NPC2, AQP1, VAPA
Predicted tissue damaging effectsOxidation/ROS generation173AMY2B, CNTN1, COX6B1, CYCS, GLB1, GNAi1, PLOD3, PNP, G6PD
Fibrosis91COX6B1, FMOD, GLB1, GNAI, PLOD3, NPC2
Predicted intracellular functionsCytoskeletal dynamics/intracellular transport135CFL1, CMT4F, DYNCH1H1, GNAI1, RAB31, SRP19, VPS13C, PIP4P2, SEPT2, VAPA
Export of intracellular calcium ions82ATP2B4, P2RX4, GNAI1
Predicted Blood physiologyBlood pressure77ACE, DYNCH1H1, GNAI1, AQP1, ATP2B4
Table 3. Significantly overrepresented proteins differentially expressed in livers from animals treated with risperidone or olanzapine.
Table 3. Significantly overrepresented proteins differentially expressed in livers from animals treated with risperidone or olanzapine.
Risperidone
UpregulatedName# KEGG PathwaysDownregulatedName# KEGG Pathways
PRKCAProtein kinase C alpha53ALDH7A1Aldehyde dehydrogenase 7, member a114
GNASGuanine nucleotide-binding protein, alpha-stimulating40COX6B1Cytochrome C oxidase subunit 6B111
PDGFRBPlatelet-derived growth factor receptor beta18AOX3Aldehyde oxidase 39
ROCK2Rho-associated coiled-coil-containing protein kinase 215ENO3Enolase 38
ITGA3Integrin alpha 311AGXTAlanine-glyoxylate and serine-pyruvate aminotransferase7
GSTA3Glutathione-S-transferase A310CSCitrate synthase oxidoreductase core subunit S17
NDUFA8/ NDUFB8Inner mitochondrial membrane complex 1 NADH:ubiquinone oxidoreductase subunit A8/B87/7NDUFS1NADH ubiquinone7
PGAM2Phosphoglycerate mutase 27ABAT4-Aminobutyrate aminotransferase7
BCAT2Branched-chain amino acid transaminase 27ALR1Aldo-keto reductase family 1 member B6
AFDNAfadin6
INPP5DInositol polyphosphate 5-phosphatase D6
ACAA2Acetyl CoA acyltransferase 26
ADIPOQAdiponectin, C1Q, and collagen domain-containing6
DGKZDiacylglycerol kinase-zeta6
EPHA2EPH-ephrin receptor A25
Olanzapine
UpregulatedName# KEGG PathwaysDownregulatedName# KEGG Pathways
GLB1Beta-galactosidase12GNAI1G-protein subunit alpha inhibitory-137
G6PDGlucose 6 phosphate dehydrogenase10CYCSSomatic cytochrome C26
ATP2B4ATPase plasma membrane Ca+2 transporting 19COX6B1Cytochrome C oxidase subunit 6B111
INPP5DInositol polyphosphate 5-phosphatase D6ENO3Enolase 38
AGPAT31-Acylglycerol-3-phosphate O-acyltransferase-35ACEAngiotensin I-converting enzyme7
CFL1Cofilin 15PNPPurine nucleoside phosphorylase6
AK2Adenylate kinase 25
ATIC5-Aminoimidazole-4-carboxamide ribonucleotide formyltransferase/IMP cyclohydrolase5
AMY2BAmylase alpha 2B5
KEGG (Kyoto Encyclopedia of Genes and Genomes) pathway analysis of risperidone- and olanzapine-associated changes in liver protein expression compared to VEH controls. #—indicates “number”. “Up” and “Down” respectively indicate increased or decreased instances of each function on KEGG pathways as a result of changes in the associated protein expression (see Supplementary File S3 for detailed correlations). RIS n = 5; OLAN n = 4; VEH n = 4. Abbreviations: RIS, risperidone; OLAN, olanzapine; OxPhos, oxidative phosphorylation; TCA, tricarboxylic acid cycle; ETC, electron transport chain; GI, gastrointestinal; ROS, reactive oxygen species.
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Rostama, B.; Beauchemin, M.; Bouchard, C.; Bernier, E.; Vary, C.P.H.; May, M.; Houseknecht, K.L. Understanding Mechanisms Underlying Non-Alcoholic Fatty Liver Disease (NAFLD) in Mental Illness: Risperidone and Olanzapine Alter the Hepatic Proteomic Signature in Mice. Int. J. Mol. Sci. 2020, 21, 9362. https://doi.org/10.3390/ijms21249362

AMA Style

Rostama B, Beauchemin M, Bouchard C, Bernier E, Vary CPH, May M, Houseknecht KL. Understanding Mechanisms Underlying Non-Alcoholic Fatty Liver Disease (NAFLD) in Mental Illness: Risperidone and Olanzapine Alter the Hepatic Proteomic Signature in Mice. International Journal of Molecular Sciences. 2020; 21(24):9362. https://doi.org/10.3390/ijms21249362

Chicago/Turabian Style

Rostama, Bahman, Megan Beauchemin, Celeste Bouchard, Elizabeth Bernier, Calvin P. H. Vary, Meghan May, and Karen L. Houseknecht. 2020. "Understanding Mechanisms Underlying Non-Alcoholic Fatty Liver Disease (NAFLD) in Mental Illness: Risperidone and Olanzapine Alter the Hepatic Proteomic Signature in Mice" International Journal of Molecular Sciences 21, no. 24: 9362. https://doi.org/10.3390/ijms21249362

APA Style

Rostama, B., Beauchemin, M., Bouchard, C., Bernier, E., Vary, C. P. H., May, M., & Houseknecht, K. L. (2020). Understanding Mechanisms Underlying Non-Alcoholic Fatty Liver Disease (NAFLD) in Mental Illness: Risperidone and Olanzapine Alter the Hepatic Proteomic Signature in Mice. International Journal of Molecular Sciences, 21(24), 9362. https://doi.org/10.3390/ijms21249362

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop