Identification of EP300 as a Key Gene Involved in Antipsychotic-Induced Metabolic Dysregulation Based on Integrative Bioinformatics Analysis of Multi-Tissue Gene Expression Data

Antipsychotics (APs) are associated with weight gain and other metabolic abnormalities such as hyperglycemia, dyslipidemia and metabolic syndrome. This translational study aimed to uncover the underlying molecular mechanisms and identify the key genes involved in AP-induced metabolic effects. An integrative gene expression analysis was performed in four different mouse tissues (striatum, liver, pancreas and adipose) after risperidone or olanzapine treatment. The analytical approach combined the identification of the gene co-expression modules related to AP treatment, gene set enrichment analysis and protein-protein interaction network construction. We found several co-expression modules of genes involved in glucose and lipid homeostasis, hormone regulation and other processes related to metabolic impairment. Among these genes, EP300, which encodes an acetyltransferase involved in transcriptional regulation, was identified as the most important hub gene overlapping the networks of both APs. Then, we explored the genetically predicted EP300 expression levels in a cohort of 226 patients with first-episode psychosis who were being treated with APs to further assess the association of this gene with metabolic alterations. The EP300 expression levels were significantly associated with increases in body weight, body mass index, total cholesterol levels, low-density lipoprotein cholesterol levels and triglyceride concentrations after 6 months of AP treatment. Taken together, our analysis identified EP300 as a key gene in AP-induced metabolic abnormalities, indicating that the dysregulation of EP300 function could be important in the development of these side effects. However, more studies are needed to disentangle the role of this gene in the mechanism of action of APs.

Introduction

Antipsychotics (APs) are the primary medication used for the treatment of schizophrenia and other psychiatric disorders. These medications, especially second-generation APs, are associated with a wide range of side effects including weight gain and other metabolic alterations that can lead to the development of metabolic syndrome, type 2 diabetes (T2D) and serious cardiovascular diseases (Fernandez-Egea et al., 2011; Foley and Morley, 2011; De Hert et al., 2012; Ward et al., 2013; Tek et al., 2016). The risk of metabolic adverse events varies among the APs, with olanzapine and clozapine having the highest risk of causing such effects and others having an intermediate (e.g., risperidone) or low risk (e.g., aripiprazole) (Lett et al., 2012; Pillinger et al., 2020). Even in the absence of AP treatment, individuals with psychiatric diseases are more likely to suffer obesity and metabolic disorders than the general population, with drug therapy contributing to the increase in these rates and aggravating the metabolic profile in these patients (Holt and Peveler, 2009; Mitchell et al., 2013; Correll et al., 2014). The associated co-morbidities reduce quality of life and increase patient mortality (Kritharides et al., 2017). Furthermore, these side effects are a major cause of treatment discontinuation, thus, worsening prognosis and increasing the risk of relapse (Doane et al., 2020). Therefore, the identification of robust predictors of AP-induced side effects could contribute to the identification of those subjects at a higher risk of developing metabolic disturbances, facilitating better treatment selection and improving clinical outcomes.

Although the underlying factors involved in AP-induced metabolic dysregulation are not well understood, genetic factors are assumed to be important based on the findings from twin and sibling studies (Gebhardt et al., 2010). Thus, several pharmacogenetic studies have been conducted, most of them using a candidate gene approach, which is limited by our current understanding of the mechanism of action of APs (Mas et al., 2012, 2015). Evidence suggests that the antagonistic effects of APs on neurotransmitter receptors involved in the control of food intake, such as the histamine H1 receptors (H1R), the serotonin 2C receptors (5-HT_2C) and the dopamine D2 receptors (D2DR), are closely linked to cardiovascular risk and metabolic dysfunction (De Hert et al., 2012; Reynolds and McGowan, 2017). However, the exact downstream mechanisms that may promote the development of metabolic abnormalities in patients receiving APs remain elusive. To date, by using this candidate gene approach, several susceptibility genes of weight gain and other metabolic adverse effects induced by APs have been identified. These genes are mostly related to neurotransmitter systems, food intake regulation and energy metabolism, and include ADR2A, BDNF, CNR1, DRD2, DRD3, FTO, GNB3, HTR2C, INSIG2, LEP, LEPR, MC4R, MTHFR, NPY, SNAP25, and SREBF1, among others (Arranz et al., 2011; Shams and Müller, 2014; Devlin and Panagiotopoulos, 2015; MacNeil and Müller, 2016; Zhang et al., 2016; Zhang and Malhotra, 2018; Gassó et al., 2020; Yoshida and Müller, 2020). A meta-analysis of 72 articles with a combined sample of 6700 patients from 46 non-overlapping cohorts found that ADRA2A, DRD2, HTR2C, and MC4R polymorphisms were some of the variants significantly associated with AP-induced weight gain, and were those with the largest effect sizes (Zhang et al., 2016). Genome-wide association studies (GWAS) have also addressed this issue, suggesting novel variants associated with AP-induced metabolic effects, such as those located in CIDEA, DGKB, MC4R, MEIS2, PRKAR2B and PTPRD, among others (Adkins et al., 2011; Malhotra et al., 2012; Brandl et al., 2016; Yu et al., 2016; Maciukiewicz et al., 2019; Corfitsen et al., 2020; ter Hark et al., 2020). Despite this associations, only few genes have been confirmed in independent studies. Hence, further research is required to identify robust markers and validate their values as predictors of AP-induced metabolic disturbances before their implementation in daily clinical practice (Arranz et al., 2021; Islam et al., 2021).

Here, we aimed to investigate the molecular pathways underlying AP-induced metabolic effects and identify potential candidate genes for future pharmacogenetic studies using a translational gene expression study from mice to humans. First, whole-genome expression changes were analyzed in different mouse tissues after treatment with olanzapine or risperidone, two APs with high and intermediate metabolic risks associated, respectively (Lett et al., 2012; Pillinger et al., 2020), that are considered the most commonly prescribed APs worldwide (Hálfdánarson et al., 2017). We focused on analyzing the gene expression and their related biological processes that are modified by both APs in order to investigate the overlapping genes and mechanisms of action involved in AP-induced metabolic imbalance. Then, we explored the genetically predicted gene expression levels of the key genes identified in a cohort of patients with first-episode psychosis (FEP) who were being treated with APs to further assess the association of these genes with metabolic alterations.

Materials and Methods

A detailed description of the methods used in the present study is available in Supplementary Material.

Animals and Drug Treatment

Mice were chronically treated for 28 days with a daily subcutaneous dose of vehicle (saline containing 5% dimethylsulfoxide (DMSO) and 5% Tween 20), risperidone (1 mg/kg) or olanzapine (3.5 mg/kg). Each group included 8 C57BL/6JOlaHsd mice (4 males and 4 females). The animals were maintained with ad libitum access to food and water, using a standard chow diet.

All animal-related procedures were performed in accordance with the European Union guidelines for the care and use of laboratory animals and were approved by the Animal Care Committee of the University of Barcelona and by the Generalitat de Catalunya.

Metabolic Assessment in Mice

Several metabolic parameters (including body weight, blood glucose levels, total cholesterol levels and triglyceride levels) were assessed in all mice at baseline (1 week before the start of AP treatment) and after 28 days of treatment. Fasting glucose levels were measured using a glucometer AlphaTRAK^® 2 (Zoetis, Parsippany, NJ, United States). Total cholesterol levels and triglyceride levels were measured using the Cholesterol Liquid assay kit and the Triglycerides Liquid assay kit, respectively (Química Clínica Aplicada S.A., Amposta, Spain). The percentage of change for each parameter between both measures was calculated. A detailed description of the methods used for the evaluation of metabolic parameters is available in Supplementary Material.

RNA Isolation and Microarray Hybridization

After 28 days of treatment, animals were sacrificed by decapitation 1 h after the last injection. Then, their tissues were rapidly removed and placed on ice. Four tissues that are potentially relevant for the development of AP-induced metabolic effects were used: the liver, the pancreas, visceral adipose tissue and the striatum. Total RNA was isolated using Trizol reagent (Life Technologies, Foster City, CA, United States) and was further purified using the miTotal RNA Extraction Miniprep System (Viogene Biotek Corp, New Taipei City, Taiwan). The RNA samples were then submitted to the Kompetenzzentrum für Fluoreszente Bioanalytik Microarray Technology (KFB, BioPark Regensburg GmbH, Regensburg, Germany) for labeling and hybridization to the Affymetrix mouse Clariom S Arrays (Affymetrix, Inc., Santa Clara, CA, United States).

Microarray Data Analysis

A weighted gene co-expression network analysis (WGCNA) was used to identify modules of highly co-expressed genes significantly associated (p-value < 0.05) with AP treatment (risperidone or olanzapine) (Langfelder and Horvath, 2008). Then, a gene set enrichment analysis of these co-expression modules was performed using FatiGO (Al-Shahrour et al., 2007) and Gene Ontology (GO) terms (Ashburner et al., 2000). From the significantly enriched (false discovery rate (FDR)-corrected p-value < 0.05) biological processes, the GO terms related to metabolic processes and metabolic dysregulation were selected for further analysis. Genes included in the selected GO terms were used to create a single protein-protein interaction (PPI) network for each AP using the SNOW program (Minguez et al., 2009). A single PPI network integrating both APs was also constructed by merging the networks obtained independently for risperidone and olanzapine using the intersection option in Cytoscape 3.7.2 (Shannon et al., 2003). For each node, the degree of connectivity was calculated to identify the candidate hub genes from the constructed PPI networks.

Full details of the extraction, labeling and hybridization protocols as well as the raw array data (.cel files) and the pre-processed data matrix are available at the Gene Expression Omnibus database (http://www.ncbi.nlm.nih.gov/geo/; accession number GSE180473).

Exploration of the Hub Genes in a Naturalistic Cohort of Patients With FEP

Participants

The sample comprised 226 patients with FEP (age 23.6 ± 6.0 years; 66.8% were males) recruited in the PEPs study. The participants were prescribed at least one second-generation AP during the 6-month follow-up period. During the follow-up period, about 40% of the patients were treated with olanzapine or clozapine, which have the highest risk of inducing weight gain or worsening metabolic parameters, while more than 50% were administered APs with an intermediate or low risk of causing such side effects, with risperidone and aripiprazole being the most frequently prescribed in this group. Further details of the study population are described in the Supplementary Material. The detailed protocol of the PEPs study has been published elsewhere (Bernardo et al., 2013; Bioque et al., 2016). The study was approved by the ethics committees of all the participating clinical centers. Informed consent was obtained from all the participants. In the case of children under 16 years of age, parents or legal guardians gave written informed consent before study participation, while the patients themselves also agreed to participate.

Metabolic Assessment

Several anthropometric and metabolic traits were assessed at baseline and at 6 months, including body weight, body mass index (BMI), blood glucose levels, total cholesterol levels, low-density lipoprotein (LDL) cholesterol levels, high-density lipoprotein (HDL) cholesterol levels and triglyceride (TG) levels. The percentage of change for each parameter between both visits was calculated (Table 1).

TABLE 1

TABLE 1. Metabolic data of the study participants at baseline and at the 6-months follow-up.

Genotyping and Gene Expression Prediction

Samples from all the individuals in the study were genotyped at the Centro Nacional de Genotipado (CeGen, Santiago de Compostela, Spain), using the Affymetrix Axiom Spain Biobank Array. The genotyping data were called using the Axiom Analysis Suite (Mas et al., 2020) and were then submitted to the Michigan Imputation Server (Das et al., 2016), following the standard pipeline for the Minimac4 software. Genotyping data were used to predict the genetically regulated gene expression levels for each individual through the gene expression imputation method, PrediXcan (Gamazon et al., 2015). Gene expression of the previously identified hub genes was predicted in five different tissues that are potentially relevant for the development of AP-induced metabolic alterations, including the small intestine, the pancreas, the liver, visceral adipose tissue and subcutaneous adipose tissue, following the standard procedure found at https://github.com/hakyimlab/PrediXcan.

Statistics

Data were analyzed using IBM SPSS Statistics version 20.0 (IBM Corp, Chicago, IL, United States). Normality was assessed using Kolmogorov-Smirnov and Shapiro-Wilk tests. Means ± standard deviations (SD) were computed for continuous variables. The differences between groups in the metabolic parameters assessed in mice were analyzed using analysis of variance. As some differences between males and females were observed in some of the parameters evaluated, the analysis was adjusted for gender.

For the analysis of the hub genes identified in the cohort of patients with FEP, the predicted gene expression levels of the EP300 gene were categorized as low, medium or high according to the three values of gene expression obtained in each tissue (−0.13, −0.07, and 0.00, respectively). Differences in metabolic variables between these categories were tested using analysis of variance. The analysis was adjusted for sociodemographic and clinical variables that might affect the parameters assessed, including gender, age, and the type of AP according to its potency (high, low, or no risk of increasing weight or worsening metabolic parameters). A further pairwise comparison analysis was performed using Bonferroni’s post-hoc test.

Results

Figure 1 shows an overview of the analysis procedure followed in the present study.

FIGURE 1

FIGURE 1. Overview of the analysis procedure followed in the study. AP: antipsychotic; BP: biological processes; GO: Gene Ontology; GSEA: gene set enrichment analysis; PPI: protein-protein interaction; WGCNA: weighted gene co-expression network analysis.

An increase in body weight (15.56% ± 5.25), fasting glucose levels (28.10% ± 43.95) and triglyceride levels (96.52% ± 89.50) was observed in mice after 28 days of treatment while cholesterol levels remained mainly unaltered (−1.56 ± 13.15). No significant differences were identified between the different treatment groups.

Regarding the gene expression analysis, 423 modules of co-expressed genes were obtained in the risperidone-treated mice: 114 in the striatum, 230 in the liver, 73 in the pancreas and six in the adipose tissue (Supplementary Figure S1). Twenty of these modules, containing between 37 and 179 genes, were significantly associated with risperidone treatment: three of those in the striatum, 13 of those in the liver and 4 of those in the pancreas (Supplementary Table S1). Regarding the olanzapine-treated mice, 483 modules of co-expressed genes were identified in the first analysis: 65 in the striatum, 223 in the liver, 140 in the pancreas and 55 in the adipose tissue (Supplementary Figure S1). Forty-six of these modules, containing between 33 and 1841 genes, were significantly associated with olanzapine treatment: five of those in the striatum, 21 of those in the liver, 15 of those in the pancreas and five of those in the adipose tissue (Supplementary Table S1). Two of the modules (red and yellow) significantly associated with olanzapine treatment in the striatum contained a large number of co-expressed genes (1,105 and 1,841 genes, respectively); hence, from these modules, only the genes with a module membership (MM) above 0.5 and a gene significance (GS) above 0.6 were selected for further analysis. Thus, 245 and 346 genes were considered for the red and yellow modules, respectively.

To gain further insight into the biological functions of the genes included in each co-expression module that was significantly associated with AP treatment, a functional enrichment analysis was performed. We found significantly enriched GO terms in 14 of the 20 co-expression modules associated with risperidone treatment (2 in the striatum, 11 in the liver and one in the pancreas) and 26 of the 46 modules associated with olanzapine treatment (4 in the striatum, 13 in the liver, five in the pancreas and 4 in the adipose tissue). For both APs, the enriched biological processes were related to different functions including tissue development, cell proliferation and differentiation, intracellular signaling, immune system processes and metabolic and biosynthetic processes, among others (Supplementary Tables S2, S3). The significantly enriched GO terms related to metabolic processes and metabolic dysregulation were selected for further analysis. These biological processes were related to glucose and lipid homeostasis, substance metabolism, digestion and the regulation of nutrient levels, hormone and secretion regulation, intracellular signaling through second messengers, the regulation of vascular tone and responses to different stimuli (Figure 2). Out of all the significant GO terms identified, 61 processes (21.7%) were selected for risperidone and 341 (25.5%) for olanzapine.

FIGURE 2

FIGURE 2. Selected biological processes related to metabolic dysregulation in each tissue of risperidone- or olanzapine-treated mice. The Y-axis represents the different categories of biological functions considered. The X-axis represents the percentage of processes (% GO terms) related to each biological function out of all the selected GO terms in each tissue for each antipsychotic. The color gradient represents the percentage of genes (% genes) related to each biological function out of all the genes included in the selected GO terms in each tissue for each antipsychotic. No significant gene modules were identified in the adipose tissue of risperidone-treated mice. None of the significant GO terms identified in the pancreas of the risperidone-treated mice were related to metabolic dysregulation.

The genes included in the selected GO terms (60 genes for risperidone and 263 for olanzapine) were used to create a PPI network for each AP (Supplementary Figure S2). The network constructed with the gene set for the risperidone-treated mice comprised 46 genes and included 29 (48.3%) of the 60 genes from the original list whereas the network constructed with the gene set for the olanzapine-treated mice comprised 295 genes, including 140 (53.2%) of the 263 genes from the original list. The nodes of both PPI networks showed more connections (degree of connectivity p-value ≤ 1 × 10⁻³), greater connectivity (clustering coefficient p-value < 0.01) and more hub nodes (betweenness centrality p-value < 0.01) than would be expected by chance. The degree of connectivity was calculated for each node and the top 10% of the genes with the highest degree were considered the hub genes of each network. Thus, seven hub genes with a degree of connectivity ranging between 5 and 10 were identified in the risperidone network, while 30 hub genes with a node degree of 11–25 were obtained in the olanzapine network.

To find the overlapping genes involved in the downstream processes regulated by both APs, a single PPI network was constructed by merging the networks obtained independently for risperidone and olanzapine (Figure 3). The resulting network showed 12 common genes, 10 of which were interconnected. The most connected gene was EP300, with a degree of connectivity of 8, thereby becoming the central gene of the network. EP300 was also identified as one of the most important hub genes in the networks of both APs, showing a node degree of 10 in the risperidone network (thus, becoming the most connected gene of the network) and a degree of 21 in the olanzapine network (becoming the second most connected gene of the network).

FIGURE 3

FIGURE 3. Protein-protein interaction network constructed by merging the networks derived independently for risperidone and olanzapine. Only the genes overlapping the networks of both antipsychotics are shown. Node size and color represent the number of connections (degree) of each gene.

To obtain more evidence of the role of EP300 as a key gene in AP-induced metabolic effects, we used the genetically predicted gene expression levels in a naturalistic cohort of patients with FEP receiving APs. Two predictive models were used to impute the gene expression levels of EP300 in different tissues based on the genotyping data: one for the prediction in the small intestine, the liver and subcutaneous adipose tissue and another one for the prediction in the pancreas and visceral adipose tissue. Then, we tested the association between the predicted gene expression values and the metabolic status of the patients at baseline and after the 6-months follow-up. No significant association was observed between the EP300 expression levels and the anthropometric and metabolic parameters measured at baseline (data not shown). By contrast, the gene expression levels of EP300, as predicted with the prediction model for the small intestine, the liver and the adipose tissue, were significantly associated with the changes in body weight (F = 4.14, p = 0.016), BMI (F = 4.48, p = 0.012), total cholesterol levels (F = 6.22, p = 0.002), LDL cholesterol levels (F = 12.53, p < 0.001) and TG levels (F = 4.01, p = 0.02) after 6 months of AP treatment. Post-hoc analysis revealed that patients with the lowest levels of EP300 expression showed higher increases in the abovementioned anthropometric and metabolic traits (p < 0.05) (Figure 4). The increase in body weight and the BMI was almost two times higher in these patients (weight increase of 14.56%; BMI increase of 15.05%) than in individuals with medium (weight increase of 7.07%, p = 0.006; BMI increase of 7.09%, p = 0.004) or high EP300 expression (weight increase of 8.76%, p = 0.054; BMI increase of 8.64%, p = 0.029). Regarding the lipid metabolic parameters, subjects with low EP300 expression levels showed an increase in cholesterol concentrations (22.81%) that was four times higher than those observed in patients with medium (5.08%, p = 0.002) or high gene expression levels (5.43%, p = 0.002). Furthermore, those with low EP300 expression levels presented an increase in LDL cholesterol levels (36.81%) that was almost 12 times higher than that of patients with medium (3.16%, p < 0.001) or high (3.06%, p < 0.001) EP300 expression. Moreover, the rise in TG concentrations was more than 2.5 times higher in individuals with the lowest EP300 expression levels (69.98%) compared to the other participants (medium EP300 expression: 25.75%, p = 0.019; high EP300 expression: 28.27%, p = 0.032).

FIGURE 4

FIGURE 4. Exploration of EP300 as a key gene in antipsychotic-induced metabolic dysregulation in patients with FEP receiving second-generation antipsychotics. Scatter box plots of the increases in the anthropometric and metabolic traits assessed according to the predicted EP300 expression levels using the prediction model of the liver, the adipose tissue and the small intestine. *p < 0.05; **p < 0.01; ***p < 0.001.

Discussion

While accumulating evidence has linked AP therapy to metabolic disturbances, the mechanistic bases and susceptibility factors of these side effects remain unclear. In the present study, we performed an integrated analysis of gene expression profiles in four different mouse tissues after treatment with risperidone or olanzapine to identify important genes involved in the development of metabolic alterations. The key genes identified were further assessed in a naturalistic cohort of patients with FEP who were being treated with APs. Our findings provide new insights into the molecular mechanism involved in AP-induced metabolic adverse events, suggesting new candidate genes for future pharmacogenetic studies.

We first performed a WGCNA of gene expression in mice that allowed us to identify several modules of co-expressed genes that were significantly associated with risperidone or olanzapine treatment in each tissue. This approach considers the overall gene–gene correlation structure instead of focusing on individual genes, providing modules of highly interconnected genes that are enriched with genes involved in known biological pathways (Zhao et al., 2010). In our study, we found that the genes in the modules were related to biological functions that were similar for the two APs, including a high number of processes related to energy metabolism, hormone regulation and the control of vascular tone, among others. The modulation of these processes by APs at a transcriptomic level may lead to the development of AP-induced metabolic disturbances.

Most proteins carry out their functions by interacting with other components that are found either in the same cell or among different cells and even among different tissues (Barabási et al., 2011). We found that the genes of these metabolic-related pathways were highly interconnected, as demonstrated by the PPI networks constructed for each AP. The risperidone and olanzapine networks were merged to identify the overlapping genes in order to obtain a general picture of the molecular effects of both APs. The central gene of this merged network was EP300, which showed connections with almost all the other genes in the network. In addition, EP300 was found to be one of the most important hub genes in the networks constructed independently for each AP. In a PPI network, nodes with a higher degree of connectivity are considered central proteins of functional importance that are essential for the maintenance of the entire network (Barabási et al., 2011). Thus, our findings suggested that EP300 could be a crucial gene in the downstream mechanism of APs that leads to metabolic dysregulation. EP300 encodes the protein p300, which functions as an acetyltransferase for histones and other proteins, thereby regulating transcriptional activity (Yao et al., 2018). This protein is involved in the downstream signaling of glucagon and insulin and participates in the regulation of energy homeostasis in major metabolic organs. Therefore, it has a crucial role in maintaining blood glucose levels in both fasting and postprandial states through the regulation of glycolysis, gluconeogenesis and glycogen synthesis (He et al., 2012; Wondisford et al., 2014; Huang et al., 2018; Yano et al., 2018; Yao et al., 2018; Namwanje et al., 2019). In addition, p300 is important for the maintenance of β-cell function, promotes glucose-induced insulin secretion and is associated with the disruption of insulin signaling in obesity (Mosley et al., 2004; Bompada et al., 2016; Cao et al., 2017; Wong et al., 2018). Therefore, p300 has been proposed to be a prime factor leading to the development of hyperglycemia, impaired insulin levels and insulin resistance (Cao et al., 2017; Yao et al., 2018), which frequently occur in individuals treated with second-generation APs (Newcomer et al., 2002; Henderson et al., 2005; Teff et al., 2013). Interestingly, a previous study revealed that some T2D-associated genetic variants identified by GWAS were located in sequences encoding protein binding sites, including those for p300 (Cheng et al., 2017). In addition, as observed in our study, earlier analyses have revealed EP300 and its paralog CBP to be the most connected genes of the networks associated with T2D (Morris et al., 2012; Randhawa et al., 2013).

By regulating gene expression, p300 is also involved in the regulation of lipid homeostasis and adipogenesis. Thus, p300 activity may enhance triglyceride synthesis and alter lipid export, thereby increasing body weight and fat mass (Yao et al., 2018; Zhou et al., 2020). APs are associated with a higher risk of these metabolic disturbances, as patients on AP treatment frequently experience dyslipidemia, mainly due to an elevation in TG levels, as well as body weight increases and central obesity (Bioque et al., 2018; Nicol et al., 2018; Wedervang-Resell et al., 2020).

We further explored the reliability of EP300 as a candidate gene for AP-induced metabolic dysregulation in a naturalistic cohort of patients with FEP receiving APs. To this end, we used the genetically predicted gene expression levels to assess the influence of EP300 on different anthropometric and metabolic traits. Gene expression is an intermediate phenotype between genetic variability and disease manifestation. Most disease-associated genetic variants are located in non-coding regions, suggesting that these variants may exert their effects by modulating gene expression (Nicolae et al., 2010). In this scenario, a new promising approach has been developed that enables the prediction of gene expression based on genetic data. The imputed gene expression can then be used in place of genotypes to assess the effect of genes on a given phenotype (Gamazon et al., 2015). This strategy has been previously used to identify candidate genes associated with schizophrenia, attention deficit hyperactivity disorder and other complex traits (Gusev et al., 2018; Liao et al., 2019; Lindström et al., 2019; Wu et al., 2021). Here, we found that the predicted expression levels of EP300 in different metabolic-related tissues of patients with FEP receiving APs were significantly associated with increases in body weight, BMI, total cholesterol levels, LDL cholesterol levels and TG concentrations at the 6-months follow-up. Although metabolic alterations have been described in schizophrenia patients even in the absence of AP treatment, we found no association between the genetically predicted EP300 expression levels and the metabolic status of FEP patients at baseline, suggesting that this gene might be related to the AP effects instead of being a risk factor for the metabolic dysregulation associated with the disease. Thus, the genetic background would affect EP300 expression levels which in turn could make patients more susceptible to the effects exerted by AP on this gene. Hence, these results give more support to the findings obtained in mice and strongly indicate the influence of EP300 on AP-induced metabolic impairment.

Previous studies have found that common and rare variants of EP300 are associated with schizophrenia risk (Ripke et al., 2014; Girard et al., 2015). However, to the best of our knowledge, this is the first time this gene has been associated with AP treatment. To date, the main mechanism of action used to explain the AP-induced metabolic adverse effects is the blockade of DRD2, H1R and 5-HT_2C receptors, leading to increased food intake and other metabolic impairments (De Hert et al., 2012; Reynolds and McGowan, 2017). However, a more complicated mechanism is expected to explain these effects. Recently, Chen et al. (2020) proposed a molecular mechanism of olanzapine-induced hyperphagia and obesity that could be in agreement with our results. They reported that the antagonistic effect of olanzapine on H1R activates GHSR1a (ghrelin receptor) downstream signaling in hypothalamic neurons that includes AMPK-FOXO1-pCREB signaling, which involves p300. This action increases the expression of neuropeptide Y, leading to excessive food intake and weight gain (Chen et al., 2020). Our results suggested that the regulation of EP300 by APs in other tissues could also contribute to the appearance of metabolic disturbances. Further studies are needed to elucidate the pleiotropic effect of EP300 on schizophrenia and AP effects and to disentangle the specific role of this gene in the mechanism of action of APs and, consequently, in AP-induced metabolic dysregulation.

The main strength of our study was the translational design, which included the integration of transcriptomic data from various mouse tissues after treatment with different APs and the subsequent exploration of the EP300 effects in a human cohort. Nevertheless, some limitations should be considered when interpreting our findings. First, the strategy used to select the biological processes that could be related to AP-induced metabolic disturbances might have been affected by bias and, therefore, we cannot rule out the possible contribution of other genes that were not included in our analysis. Secondly, although drug dosage regimens were selected based on previous studies, here we do not control for drug plasma levels as we were unable to obtain enough blood volume of each animal to measure plasmatic drug levels with the methodology available in our facilities. Besides, the lack of a mouse group treated with an AP with a low risk of metabolic dysregulation meant that the specific effect of high and intermediate risk APs could not be explored. In addition, as no significant differences were observed in the metabolic parameters assessed between treatment groups in mice, we were not able to correlate these variables with AP-induced gene expression changes. This lack of significant differences may be due to the experimental design, in which mice were fed using a standard chow diet. Although a high-fat diet could mimic better the clinical effects of APs on metabolic parameters, which are in part due to changes in the food intake profile, we maintain the animals on a standard diet to test specifically the AP-induced gene expression changes, without the influence of external factors. Despite these limitations, we could ensure the influence of the key gene identified on metabolic impairment by assessing the effect of the genetically predicted gene expression levels in a human cohort. This exploration was conducted in a cohort of patients in a naturalistic study, and hence the sociodemographic and clinical characteristics of these individuals, such as the type of antipsychotic prescribed, could have affected the parameters assessed. For this reason, as the sample size did not allow the study of subgroups stratified by AP treatment, the statistical analysis was adjusted for several variables, including gender, age and AP potency. However, other confounders such as diagnosis, concomitant treatment or diet that we did not control for could have had an effect on the metabolic traits assessed. Finally, the prediction of the gene expression levels in this population was limited by the prediction models currently available, which in the future could include more genetic variants based on scientific advances. Hence, our results should be interpreted as exploratory and require replication.

In summary, our findings provide novel insights into the molecular mechanisms involved in the metabolic side effects of APs. Based on our results, we propose that the dysregulation of EP300 functioning could be important in the development of AP-induced metabolic disturbances. Hence, both EP300 and the genes involved in the EP300 network could be considered potential candidate genes for future pharmacogenetic studies on these side effects of APs.

Data Availability Statement

The datasets presented in this study can be found in online repositories. The names of the repository/repositories and accession number(s) can be found below: NCBI GEO, accession no: GSE180473.

Ethics Statement

The studies involving human participants were reviewed and approved by the Comité Ético de Investigación Clínica del Hospital Clínic de Barcelona. Written informed consent to participate in this study was provided by the participants or their legal guardian/next of kin. The animal study was reviewed and approved by the Animal Care Committee of the University of Barcelona and Generalitat de Catalunya.

Author Contributions

AM-P participated designing and carrying out the experimental procedures, contributed to data acquisition and analysis and helped in the drafting of the manuscript. PG participated in data analysis and the interpretation of results and writing the first version of the manuscript. LP and AS participated in data analysis and helped in the drafting of the manuscript. MP, JS-R, MC, and MB provided clinical samples, collected clinical information and helped in the drafting of the manuscript. MB and AL acquired funding for the research, participated in designing the study and helped in the drafting of the manuscript. SM and NR participated in designing the study, supervised the data analysis and contributed to the interpretation of results and the writing of the first version of the manuscript. All authors revised the manuscript and approved the final version for publication.

Funding

This study was supported by the Ministerio de Economía y Competitividad, Instituto de Salud Carlos III, Fondo de Investigación Sanitaria (FIS) (PI16/01,222 and PI08/0208); and Fondo Europeo de Desarrollo Regional (FEDER)-Unión Europea. Support was also given by the “Agència de Gestió d'Ajuts Universitaris i Recerca” (AGAUR) of the “Generalitat de Catalunya” (2017SGR1562, 2017SGR1355 and 2017SGR881) and with the grant FI-SDUR 2019-2020 (AGAUR) (2020 FI-SDUR 00327) to AS.

Conflict of Interest

MP has been a consultant for and received educational honoraria and travel grants from Exeltis, Janssen-Cilag, Otsuka and Servier. JS-R has been a speaker for and on the advisory boards of Lilly, GlaxoSmithKline, Lundbeck, Janssen, Servier, and Pfizer and has received grant/honoraria from Lilly and Astra-Zeneca. MB has been a consultant for, received grant/research support and honoraria from, and been on the speakers/advisory board of ABBiotics, Adamed, Almirall, AMGEN, Boehringer, Eli Lilly, Ferrer, Forum Pharmaceuticals, Gedeon, Hersill, Janssen-Cilag, Lundbeck, Otsuka, Pfizer, Roche, and Servier.

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

Publisher’s Note

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.

Acknowledgments

The authors thank the Language Advisory Service at the University of Barcelona for manuscript revision. The authors also thank all subjects and their families for the time and effort spent on this study as well as Ana Meseguer for sample collection assistance.

Supplementary Material

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

References

Adkins, D. E., Åberg, K., McClay, J. L., Bukszár, J., Zhao, Z., Jia, P., et al. (2011). Genomewide Pharmacogenomic Study of Metabolic Side Effects to Antipsychotic Drugs. Mol. Psychiatry 16, 321–332. doi:10.1038/mp.2010.14

PubMed Abstract | CrossRef Full Text | Google Scholar

Al-Shahrour, F., Minguez, P., Tárraga, J., Medina, I., Alloza, E., Montaner, D., et al. (2007). FatiGO +: A Functional Profiling Tool for Genomic Data. Integration of Functional Annotation, Regulatory Motifs and Interaction Data with Microarray Experiments. Nucleic Acids Res. 35, W91–W96. doi:10.1093/nar/gkm260

PubMed Abstract | CrossRef Full Text | Google Scholar

Arranz, M. J., Rivera, M., and Munro, J. C. (2011). Pharmacogenetics of Response to Antipsychotics in Patients with Schizophrenia. CNS Drugs 25, 933–969. doi:10.2165/11595380-000000000-00000

PubMed Abstract | CrossRef Full Text | Google Scholar

Arranz, M. J., Salazar, J., and Hernández, M. H. (2021). Pharmacogenetics of Antipsychotics: Clinical Utility and Implementation. Behav. Brain Res. 401, 113058. doi:10.1016/j.bbr.2020.113058

PubMed Abstract | CrossRef Full Text | Google Scholar

Ashburner, M., Ball, C. A., Blake, J. A., Botstein, D., Butler, H., Cherry, J. M., et al. (2000). Gene Ontology: Tool for the Unification of Biology. The Gene Ontology Consortium. Nat. Genet. 25, 25–29. doi:10.1038/75556

PubMed Abstract | CrossRef Full Text | Google Scholar

Barabási, A. L., Gulbahce, N., and Loscalzo, J. (2011). Network Medicine: A Network-Based Approach to Human Disease. Nat. Rev. Genet. 12, 56–68. doi:10.1038/nrg2918

PubMed Abstract | CrossRef Full Text | Google Scholar

Bernardo, M., Bioque, M., Parellada, M., Saiz Ruiz, J., Cuesta, M. J., Llerena, A., et al. (2013). Assessing Clinical and Functional Outcomes in a Gene-Environment Interaction Study in First Episode of Psychosis (PEPs). Rev. Psiquiatr Salud Ment. 6, 4–16. doi:10.1016/j.rpsm.2012.11.001

PubMed Abstract | CrossRef Full Text | Google Scholar

Bioque, M., García-Portilla, M. A. P., García-Rizo, C., Cabrera, B., Lobo, A., González-Pinto, A., et al. (2018). Evolution of Metabolic Risk Factors over a Two-Year Period in a Cohort of First Episodes of Psychosis. Schizophr. Res. 193, 188–196. doi:10.1016/j.schres.2017.06.032

PubMed Abstract | CrossRef Full Text | Google Scholar

Bioque, M., Llerena, A., Cabrera, B., Mezquida, G., Lobo, A., González-Pinto, A., et al. (2016). A Pharmacovigilance Study in First Episode of Psychosis: Psychopharmacological Interventions and Safety Profiles in the PEPs Project. Int. J. Neuropsychopharmacol. 19, 1–10. doi:10.1093/ijnp/pyv121

PubMed Abstract | CrossRef Full Text | Google Scholar

Bompada, P., Atac, D., Luan, C., Andersson, R., Omella, J. D., Laakso, E. O., et al. (2016). Histone Acetylation of Glucose-Induced Thioredoxin-Interacting Protein Gene Expression in Pancreatic Islets. Int. J. Biochem. Cel Biol. 81, 82–91. doi:10.1016/j.biocel.2016.10.022

CrossRef Full Text | Google Scholar

Brandl, E. J., Tiwari, A. K., Zai, C. C., Nurmi, E. L., Chowdhury, N. I., Arenovich, T., et al. (2016). Genome-wide Association Study on Antipsychotic-Induced Weight Gain in the CATIE Sample. Pharmacogenomics J. 16, 352–356. doi:10.1038/tpj.2015.59

PubMed Abstract | CrossRef Full Text | Google Scholar

Cao, J., Peng, J., An, H., He, Q., Boronina, T., Guo, S., et al. (2017). Endotoxemia-mediated Activation of Acetyltransferase P300 Impairs Insulin Signaling in Obesity. Nat. Commun. 8, 131. doi:10.1038/s41467-017-00163-w

PubMed Abstract | CrossRef Full Text | Google Scholar

Chen, X., Yu, Y., Zheng, P., Jin, T., He, M., Zheng, M., et al. (2020). Olanzapine Increases AMPK-NPY Orexigenic Signaling by Disrupting H1R-GHSR1a Interaction in the Hypothalamic Neurons of Mice. Psychoneuroendocrinology 114, 104594. doi:10.1016/j.psyneuen.2020.104594

PubMed Abstract | CrossRef Full Text | Google Scholar

Cheng, M., Liu, X., Yang, M., Han, L., Xu, A., and Huang, Q. (2017). Computational Analyses of Type 2 Diabetes-Associated Loci Identified by Genome-wide Association Studies. J. Diabetes 9, 362–377. doi:10.1111/1753-0407.12421

PubMed Abstract | CrossRef Full Text | Google Scholar

Corfitsen, H. T., Krantz, B., Larsen, A., and Drago, A. (2020). Molecular Pathway Analysis Associates Alterations in Obesity-Related Genes and Antipsychotic-Induced Weight Gain. Acta Neuropsychiatr. 32, 72–83. doi:10.1017/neu.2019.41

PubMed Abstract | CrossRef Full Text | Google Scholar

Correll, C. U., Robinson, D. G., Schooler, N. R., Brunette, M. F., Mueser, K. T., Rosenheck, R. A., et al. (2014). Cardiometabolic Risk in Patients with First-Episode Schizophrenia Spectrum Disorders: Baseline Results from the RAISE-ETP Study. JAMA Psychiatry 71, 1350–1363. doi:10.1001/jamapsychiatry.2014.1314

PubMed Abstract | CrossRef Full Text | Google Scholar

Das, S., Forer, L., Schönherr, S., Sidore, C., Locke, A. E., Kwong, A., et al. (2016). Next-generation Genotype Imputation Service and Methods. Nat. Genet. 48, 1284–1287. doi:10.1038/ng.3656

PubMed Abstract | CrossRef Full Text | Google Scholar

De Hert, M., Detraux, J., Van Winkel, R., Yu, W., and Correll, C. U. (2012). Metabolic and Cardiovascular Adverse Effects Associated with Antipsychotic Drugs. Nat. Rev. Endocrinol. 8, 114–126. doi:10.1038/nrendo.2011.156

CrossRef Full Text | Google Scholar

Devlin, A. M., and Panagiotopoulos, C. (2015). Metabolic Side Effects and Pharmacogenetics of Second-Generation Antipsychotics in Children. Pharmacogenomics 16, 981–996. doi:10.2217/pgs.15.55

PubMed Abstract | CrossRef Full Text | Google Scholar

Doane, M. J., Sajatovic, M., Weiden, P. J., O'Sullivan, A. K., Maher, S., Bjorner, J. B., et al. (2020). Antipsychotic Treatment Experiences of People with Schizophrenia: Patient Perspectives from an Online Survey. Patient Prefer. Adherence 14, 2043–2054. doi:10.2147/PPA.S270020

PubMed Abstract | CrossRef Full Text | Google Scholar

Fernandez-Egea, E., Miller, B., Garcia-Rizo, C., Bernardo, M., and Kirkpatrick, B. (2011). Metabolic Effects of Olanzapine in Patients with Newly Diagnosed Psychosis. J. Clin. Psychopharmacol. 31, 154–159. doi:10.1097/JCP.0b013e31820fcea3

PubMed Abstract | CrossRef Full Text | Google Scholar

Foley, D. L., and Morley, K. I. (2011). Systematic Review of Early Cardiometabolic Outcomes of the First Treated Episode of Psychosis. Arch. Gen. Psychiatry 68, 609–616. doi:10.1001/archgenpsychiatry.2011.2

PubMed Abstract | CrossRef Full Text | Google Scholar

Gamazon, E. R., Wheeler, H. E., Shah, K. P., Mozaffari, S. V., Aquino-Michaels, K., Carroll, R. J., et al. (2015). A Gene-Based Association Method for Mapping Traits Using Reference Transcriptome Data. Nat. Genet. 47, 1091–1098. doi:10.1038/ng.3367

PubMed Abstract | CrossRef Full Text | Google Scholar

Gassó, P., Arnaiz, J. A., Mas, S., Lafuente, A., Bioque, M., Cuesta, M. J., et al. (2020). Association Study of Candidate Genes with Obesity and Metabolic Traits in Antipsychotic-Treated Patients with First-Episode Psychosis over a 2-year Period. J. Psychopharmacol. 34, 514–523. doi:10.1177/0269881120903462

PubMed Abstract | CrossRef Full Text | Google Scholar

Gebhardt, S., Theisen, F. M., Haberhausen, M., Heinzel-Gutenbrunner, M., Wehmeier, P. M., Krieg, J. C., et al. (2010). Body Weight Gain Induced by Atypical Antipsychotics: an Extension of the Monozygotic Twin and Sib Pair Study. J. Clin. Pharm. Ther. 35, 207–211. doi:10.1111/j.1365-2710.2009.01084.x

PubMed Abstract | CrossRef Full Text | Google Scholar

Girard, S. L., Dion, P. A., Bourassa, C. V., Geoffroy, S., Lachance-Touchette, P., Barhdadi, A., et al. (2015). Mutation burden of Rare Variants in Schizophrenia Candidate Genes. PLoS One 10, e0128988. doi:10.1371/journal.pone.0128988

PubMed Abstract | CrossRef Full Text | Google Scholar

Gusev, A., Mancuso, N., Won, H., Kousi, M., Finucane, H. K., Reshef, Y., et al. (2018). Transcriptome-wide Association Study of Schizophrenia and Chromatin Activity Yields Mechanistic Disease Insights. Nat. Genet. 50, 538–548. doi:10.1038/s41588-018-0092-1

PubMed Abstract | CrossRef Full Text | Google Scholar

Hálfdánarson, Ó., Zoëga, H., Aagaard, L., Bernardo, M., Brandt, L., Fusté, A. C., et al. (2017). International Trends in Antipsychotic Use: A Study in 16 Countries, 2005-2014. Eur. Neuropsychopharmacol. 27, 1064–1076. doi:10.1016/j.euroneuro.2017.07.001

PubMed Abstract | CrossRef Full Text | Google Scholar

He, L., Naik, K., Meng, S., Cao, J., Sidhaye, A. R., Ma, A., et al. (2012). Transcriptional Co-activator P300 Maintains Basal Hepatic Gluconeogenesis. J. Biol. Chem. 287, 32069–32077. doi:10.1074/jbc.M112.385864

PubMed Abstract | CrossRef Full Text | Google Scholar

Henderson, D. C., Cagliero, E., Copeland, P. M., Borba, C. P., Evins, A. E., Hayden, D., et al. (2005). Glucose Metabolism in Patients with Schizophrenia Treated with Atypical Antipsychotic Agents: a Frequently Sampled Intravenous Glucose Tolerance Test and Minimal Model Analysis. Arch. Gen. Psychiatry 62, 19–28. doi:10.1001/archpsyc.62.1.19

PubMed Abstract | CrossRef Full Text | Google Scholar

Holt, R. I., and Peveler, R. C. (2009). Obesity, Serious Mental Illness and Antipsychotic Drugs. Diabetes Obes. Metab. 11, 665–679. doi:10.1111/j.1463-1326.2009.01038.x

PubMed Abstract | CrossRef Full Text | Google Scholar

Huang, H., Tang, S., Ji, M., Tang, Z., Shimada, M., Liu, X., et al. (2018). p300-mediated Lysine 2-hydroxyisobutyrylation Regulates Glycolysis. Mol. Cel. 70, 984–678. doi:10.1016/j.molcel.2018.05.035

PubMed Abstract | CrossRef Full Text | Google Scholar

Islam, F., Men, X., Yoshida, K., Zai, C. C., and Müller, D. J. (2021). Pharmacogenetics‐guided Advances in Antipsychotic Treatment. Clin. Pharmacol. Ther. 2021, 1. doi:10.1002/cpt.2339

CrossRef Full Text | Google Scholar

Kritharides, L., Chow, V., and Lambert, T. J. (2017). Cardiovascular Disease in Patients with Schizophrenia. Med. J. Aust. 207, 179–195. doi:10.5694/mja17.00258

PubMed Abstract | CrossRef Full Text | Google Scholar

Langfelder, P., and Horvath, S. (2008). WGCNA: An R Package for Weighted Correlation Network Analysis. BMC Bioinformatics 9, 559. doi:10.1186/1471-2105-9-559

PubMed Abstract | CrossRef Full Text | Google Scholar

Lett, T. A., Wallace, T. J., Chowdhury, N. I., Tiwari, A. K., Kennedy, J. L., and Müller, D. J. (2012). Pharmacogenetics of Antipsychotic-Induced Weight Gain: Review and Clinical Implications. Mol. Psychiatry 17, 242–266. doi:10.1038/mp.2011.109

PubMed Abstract | CrossRef Full Text | Google Scholar

Liao, C., Laporte, A. D., Spiegelman, D., Akçimen, F., Joober, R., Dion, P. A., et al. (2019). Transcriptome-wide Association Study of Attention Deficit Hyperactivity Disorder Identifies Associated Genes and Phenotypes. Nat. Commun. 10, 4450–4457. doi:10.1038/s41467-019-12450-9

PubMed Abstract | CrossRef Full Text | Google Scholar

Lindström, S., Wang, L., Smith, E. N., Gordon, W., Van Hylckama Vlieg, A., De Andrade, M., et al. (2019). Genomic and Transcriptomic Association Studies Identify 16 Novel Susceptibility Loci for Venous Thromboembolism. Blood 134, 1645–1657. doi:10.1182/blood.2019000435

PubMed Abstract | CrossRef Full Text | Google Scholar

Maciukiewicz, M., Tiwari, A. K., Zai, C. C., Gorbovskaya, I., Laughlin, C. P., Nurmi, E. L., et al. (2019). Genome-wide Association Study on Antipsychotic-Induced Weight Gain in Europeans and African-Americans. Schizophr. Res. 212, 204–212. doi:10.1016/j.schres.2019.07.022

PubMed Abstract | CrossRef Full Text | Google Scholar

MacNeil, R. R., and Müller, D. J. (2016). Genetics of Common Antipsychotic-Induced Adverse Effects. Mol. Neuropsychiatry 2, 61–78. doi:10.1159/000445802

PubMed Abstract | CrossRef Full Text | Google Scholar

Malhotra, A. K., Correll, C. U., Chowdhury, N. I., Müller, D. J., Gregersen, P. K., Lee, A. T., et al. (2012). Association between Common Variants Near the Melanocortin 4 Receptor Gene and Severe Antipsychotic Drug-Induced Weight Gain. Arch. Gen. Psychiatry 69, 904–912. doi:10.1001/archgenpsychiatry.2012.191

PubMed Abstract | CrossRef Full Text | Google Scholar

Mas, S., Boloc, D., Rodríguez, N., Mezquida, G., Amoretti, S., Cuesta, M. J., et al. (2020). Examining Gene-Environment Interactions Using Aggregate Scores in a First-Episode Psychosis Cohort. Schizophr. Bull. 46, 1019–1025. doi:10.1093/schbul/sbaa012

PubMed Abstract | CrossRef Full Text | Google Scholar

Mas, S., Gassó, P., and Lafuente, A. (2015). Applicability of Gene Expression and Systems Biology to Develop Pharmacogenetic Predictors; Antipsychotic-Induced Extrapyramidal Symptoms as an Example. Pharmacogenomics 16, 1975–1988. doi:10.2217/pgs.15.134

PubMed Abstract | CrossRef Full Text | Google Scholar

Mas, S., Llerena, A., Saíz, J., Bernardo, M., and Lafuente, A. (2012). Strengths and Weaknesses of Pharmacogenetic Studies of Antipsychotic Drugs: The Potential Value of the PEPs Study. Pharmacogenomics 13, 1773–1782. doi:10.2217/pgs.12.159

PubMed Abstract | CrossRef Full Text | Google Scholar

Minguez, P., Götz, S., Montaner, D., Al-Shahrour, F., and Dopazo, J. (2009). SNOW, a Web-Based Tool for the Statistical Analysis of Protein-Protein Interaction Networks. Nucleic Acids Res. 37, W109–W114. doi:10.1093/nar/gkp402

PubMed Abstract | CrossRef Full Text | Google Scholar

Mitchell, A. J., Vancampfort, D., De Herdt, A., Yu, W., and De Hert, M. (2013). Is the Prevalence of Metabolic Syndrome and Metabolic Abnormalities Increased in Early Schizophrenia? A Comparative Meta-Analysis of First Episode, Untreated and Treated Patients. Schizophr. Bull. 39, 295–305. doi:10.1093/schbul/sbs082

PubMed Abstract | CrossRef Full Text | Google Scholar

Morris, A. P., Voight, B. F., Teslovich, T. M., Ferreira, T., Segrè, A. V., Steinthorsdottir, V., et al. (2012). Large-scale Association Analysis Provides Insights into the Genetic Architecture and Pathophysiology of Type 2 Diabetes. Nat. Genet. 44, 981–990. doi:10.1038/ng.2383

PubMed Abstract | CrossRef Full Text | Google Scholar

Mosley, A. L., Corbett, J. A., and Özcan, S. (2004). Glucose Regulation of Insulin Gene Expression Requires the Recruitment of P300 by the Beta-cell-specific Transcription Factor Pdx-1. Mol. Endocrinol. 18, 2279–2290. doi:10.1210/me.2003-0463

PubMed Abstract | CrossRef Full Text | Google Scholar

Namwanje, M., Liu, L., Chan, M., Aaron, N., Kraakman, M. J., and Qiang, L. (2019). The Depot-specific and Essential Roles of CBP/p300 in Regulating Adipose Plasticity. J. Endocrinol. 240, 257–269. doi:10.1530/JOE-18-0361

PubMed Abstract | CrossRef Full Text | Google Scholar

Newcomer, J. W., Haupt, D. W., Fucetola, R., Melson, A. K., Schweiger, J. A., Cooper, B. P., et al. (2002). Abnormalities in Glucose Regulation during Antipsychotic Treatment of Schizophrenia. Arch. Gen. Psychiatry 59, 337–345. doi:10.1001/archpsyc.59.4.337

PubMed Abstract | CrossRef Full Text | Google Scholar

Nicol, G. E., Yingling, M. D., Flavin, K. S., Schweiger, J. A., Patterson, B. W., Schechtman, K. B., et al. (2018). Metabolic Effects of Antipsychotics on Adiposity and Insulin Sensitivity in Youths: A Randomized Clinical Trial. JAMA Psychiatry 75, 788–796. doi:10.1001/jamapsychiatry.2018.1088

PubMed Abstract | CrossRef Full Text | Google Scholar

Nicolae, D. L., Gamazon, E., Zhang, W., Duan, S., Dolan, M. E., and Cox, N. J. (2010). Trait-associated SNPs Are More Likely to Be eQTLs: Annotation to Enhance Discovery from GWAS. Plos Genet. 6, e1000888. doi:10.1371/journal.pgen.1000888

PubMed Abstract | CrossRef Full Text | Google Scholar

Pillinger, T., McCutcheon, R. A., Vano, L., Mizuno, Y., Arumuham, A., Hindley, G., et al. (2020). 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 7, 64–77. doi:10.1016/S2215-0366(19)30416-X

PubMed Abstract | CrossRef Full Text | Google Scholar

Randhawa, V., Sharma, P., Bhushan, S., and Bagler, G. (2013). Identification of Key Nodes of Type 2 Diabetes Mellitus Protein Interactome and Study of Their Interactions with Phloridzin. OMICS 17, 302–317. doi:10.1089/omi.2012.0115

PubMed Abstract | CrossRef Full Text | Google Scholar

Reynolds, G. P., and McGowan, O. O. (2017). Mechanisms Underlying Metabolic Disturbances Associated with Psychosis and Antipsychotic Drug Treatment. J. Psychopharmacol. 31, 1430–1436. doi:10.1177/0269881117722987

PubMed Abstract | CrossRef Full Text | Google Scholar

Ripke, S., Neale, B. M., Corvin, A., Walters, J. T., Farth, K.-H., Holmans, P. A., et al. (2014). Biological Insights from 108 Schizophrenia-Associated Genetic Loci. Nature 511, 421–427. doi:10.1038/nature13595

PubMed Abstract | CrossRef Full Text | Google Scholar

Shams, T. A., and Müller, D. J. (2014). Antipsychotic Induced Weight Gain: Genetics, Epigenetics, and Biomarkers Reviewed. Curr. Psychiatry Rep. 16, 473. doi:10.1007/s11920-014-0473-9

PubMed Abstract | CrossRef Full Text | Google Scholar

Shannon, P., Markiel, A., Ozier, O., Baliga, N. S., Wang, J. T., Ramage, D., et al. (2003). Cytoscape: a Software Environment for Integrated Models of Biomolecular Interaction Networks. Genome Res. 13, 2498–2504. doi:10.1101/gr.1239303

PubMed Abstract | CrossRef Full Text | Google Scholar

Teff, K. L., Rickels, M. R., Grudziak, J., Fuller, C., Nguyen, H. L., and Rickels, K. (2013). Antipsychotic-induced Insulin Resistance and Postprandial Hormonal Dysregulation Independent of Weight Gain or Psychiatric Disease. Diabetes 62, 3232–3240. doi:10.2337/db13-0430

PubMed Abstract | CrossRef Full Text | Google Scholar

Tek, C., Kucukgoncu, S., Guloksuz, S., Woods, S. W., Srihari, V. H., and Annamalai, A. (2016). Antipsychotic-induced Weight Gain in First-Episode Psychosis Patients: A Meta-Analysis of Differential Effects of Antipsychotic Medications. Early Interv. Psychiatry 10, 193–202. doi:10.1111/eip.12251

PubMed Abstract | CrossRef Full Text | Google Scholar

ter Hark, S. E., Jamain, S., Schijven, D., Lin, B. D., Bakker, M. K., Boland-Auge, A., et al. (2020). A New Genetic Locus for Antipsychotic-Induced Weight Gain: A Genome-wide Study of First-Episode Psychosis Patients Using Amisulpride (From the OPTiMiSE Cohort). J. Psychopharmacol. 34, 524–531. doi:10.1177/0269881120907972

PubMed Abstract | CrossRef Full Text | Google Scholar

Ward, A., Quon, P., Abouzaid, S., Haber, N., Ahmed, S., and Kim, E. (2013). Cardiometabolic Consequences of Therapy for Chronic Schizophrenia Using Second-Generation Antipsychotic Agents in a Medicaid Population: Clinical and Economic Evaluation. P T 38, 109–115.

PubMed Abstract | Google Scholar

Wedervang-Resell, K., Friis, S., Lonning, V., Smelror, R. E., Johannessen, C., Agartz, I., et al. (2020). Lipid Alterations in Adolescents with Early-Onset Psychosis May Be Independent of Antipsychotic Medication. Schizophr. Res. 216, 295–301. doi:10.1016/j.schres.2019.11.039

PubMed Abstract | CrossRef Full Text | Google Scholar

Wondisford, A. R., Xiong, L., Chang, E., Meng, S., Meyers, D. J., Li, M., et al. (2014). Control of Foxo1 Gene Expression by Co-activator P300. J. Biol. Chem. 289, 4326–4333. doi:10.1074/jbc.M113.540500

PubMed Abstract | CrossRef Full Text | Google Scholar

Wong, C. K., Wade-Vallance, A. K., Luciani, D. S., Brindle, P. K., Lynn, F. C., and Gibson, W. T. (2018). The P300 and CBP Transcriptional Coactivators Are Required for β-Cell and α-Cell Proliferation. Diabetes 67, 412–422. doi:10.2337/db17-0237

PubMed Abstract | CrossRef Full Text | Google Scholar

Wu, C., Tan, S., Liu, L., Cheng, S., Li, P., Li, W., et al. (2021). Transcriptome-wide Association Study Identifies Susceptibility Genes for Rheumatoid Arthritis. Arthritis Res. Ther. 23, 38. doi:10.1186/s13075-021-02419-9

PubMed Abstract | CrossRef Full Text | Google Scholar

Yano, H., Sakai, M., Matsukawa, T., Yagi, T., Naganuma, T., Mitsushima, M., et al. (2018). PHD3 Regulates Glucose Metabolism by Suppressing Stress-Induced Signalling and Optimising Gluconeogenesis and Insulin Signalling in Hepatocytes. Sci. Rep. 8, 14290. doi:10.1038/s41598-018-32575-z

PubMed Abstract | CrossRef Full Text | Google Scholar

Yao, W., Wang, T., and Huang, F. (2018). P300/CBP as a Key Nutritional Sensor for Hepatic Energy Homeostasis and Liver Fibrosis. Biomed. Res. Int. 2018, 8168791. doi:10.1155/2018/8168791

PubMed Abstract | CrossRef Full Text | Google Scholar

Yoshida, K., and Müller, D. J. (2020). Pharmacogenetics of Antipsychotic Drug Treatment: Update and Clinical Implications. Mol. Neuropsychiatry 5, 1–26. doi:10.1159/000492332

PubMed Abstract | CrossRef Full Text | Google Scholar

Yu, H., Wang, L., Lv, L., Ma, C., Du, B., Lu, T., et al. (2016). Genome-wide Association Study Suggested the PTPRD Polymorphisms Were Associated with Weight Gain Effects of Atypical Antipsychotic Medications. Schizophr. Bull. 42, 814–823. doi:10.1093/schbul/sbv179

PubMed Abstract | CrossRef Full Text | Google Scholar

Zhang, J. P., Lencz, T., Zhang, R. X., Nitta, M., Maayan, L., John, M., et al. (2016). Pharmacogenetic Associations of Antipsychotic Drug-Related Weight Gain: a Systematic Review and Meta-Analysis. Schizophr. Bull. 42, 1418–1437. doi:10.1093/schbul/sbw058

PubMed Abstract | CrossRef Full Text | Google Scholar

Zhang, J. P., and Malhotra, A. K. (2018). Recent Progress in Pharmacogenomics of Antipsychotic Drug Response. Curr. Psychiatry Rep. 20, 24. doi:10.1007/s11920-018-0886-y

PubMed Abstract | CrossRef Full Text | Google Scholar

Zhao, W., Langfelder, P., Fuller, T., Dong, J., Li, A., and Hovarth, S. (2010). Weighted Gene Coexpression Network Analysis: State of the Art. J. Biopharm. Stat. 20, 281–300. doi:10.1080/10543400903572753

PubMed Abstract | CrossRef Full Text | Google Scholar

Zhou, F., Liu, Q., Zhang, L., Zhu, Q., Wang, S., Zhu, K., et al. (2020). Selective Inhibition of CBP/p300 HAT by A-485 Results in Suppression of Lipogenesis and Hepatic Gluconeogenesis. Cell Death Dis. 11, 745. doi:10.1038/s41419-020-02960-6

PubMed Abstract | CrossRef Full Text | Google Scholar