- 1Department of Molecular Pharmacology, Graduate School of Pharmaceutical Sciences, Kyoto University, Kyoto, Japan
- 2Division of Pharmacology, Department of Biomedical Sciences, Nihon University School of Medicine, Tokyo, Japan
Olanzapine is an atypical antipsychotic drug that is clinically applied in patients with schizophrenia. It increases the risk of dyslipidemia, a disturbance of lipid metabolic homeostasis, usually characterized by increased low-density lipoprotein (LDL) cholesterol and triglycerides, and accompanied by decreased high-density lipoprotein (HDL) in the serum. In this study, analyzing the FDA Adverse Event Reporting System, JMDC insurance claims, and electronic medical records from Nihon University School of Medicine revealed that a co-treated drug, vitamin D, can reduce the incidence of olanzapine-induced dyslipidemia. In the following experimental validations of this hypothesis, short-term oral olanzapine administration in mice caused a simultaneous increase and decrease in the levels of LDL and HDL cholesterol, respectively, while the triglyceride level remained unaffected. Cholecalciferol supplementation attenuated these deteriorations in blood lipid profiles. RNA-seq analysis was conducted on three cell types that are closely related to maintaining cholesterol metabolic balance (hepatocytes, adipocytes, and C2C12) to verify the direct effects of olanzapine and the functional metabolites of cholecalciferol (calcifediol and calcitriol). Consequently, the expression of cholesterol-biosynthesis-related genes was reduced in calcifediol- and calcitriol-treated C2C12 cells, which was likely to be mediated by activating the vitamin D receptor that subsequently inhibited the cholesterol biosynthesis process via insulin-induced gene 2 regulation. This clinical big-data-based drug repurposing approach is effective in finding a novel treatment with high clinical predictability and a well-defined molecular mechanism.
Introduction
Schizophrenia is a group of chronic psychiatric disorders characterized by hallucinations, delusions, reduced motivation, and blunt affect (Saha et al., 2005). Atypical antipsychotics are currently used as first-line drugs in patients with schizophrenia because they potently improve positive symptoms without causing extrapyramidal symptoms (Campbell et al., 1999). Among these, olanzapine is one of the most widely used drugs (Kochi et al., 2017). However, olanzapine causes dyslipidemia as a significant adverse effect, characterized by an increase in total blood cholesterol and triglyceride levels (Lieberman et al., 2005). A recent meta-analysis further revealed a relationship between olanzapine treatment and an increase in low-density lipoprotein (LDL) cholesterol in the blood (Li et al., 2020). Although olanzapine has excellent cardiac safety (Kennedy et al., 2001), the obesity induced by olanzapine treatment has been proven to be a cardiovascular risk factor (Wang et al., 2014; Correll et al., 2015), thereby making its treatment a significant risk factor that contributes to excess premature mortality in patients with schizophrenia (Brown et al., 2000). At present, preventing this adverse effect is impossible due to the incomplete understanding of the underlying mechanisms of olanzapine-induced dyslipidemia (Yan et al., 2013). Thus, there exists a clinical demand for novel target-screening approaches.
Drug repurposing from market-accessible drugs is a classic approach for refining the potential therapeutic applications. Repurposing studies have shown that topiramate, a newer anticonvulsant, was able to attenuate the obesity-inducing effect of olanzapine (Narula et al., 2010), and diabetes-treating drugs, such as metformin (Chen et al., 2008) and liraglutide (Larsen et al., 2017) have relieved the excessive body weight gain due to olanzapine treatment. However, a recently published meta-analysis review on the clinical trials focusing on blood lipid levels concluded that the available lipid-lowering agents are not effective in treating patients who are prescribed with atypical antipsychotics, including olanzapine (Kanagasundaram et al., 2021). Therefore, current practices need to be optimized, considering that existing drug repurposing approaches are inefficient in refining interventions with high clinical efficiency for treating olanzapine-induced dyslipidemia.
In this study, an updated approach, “reverse translational drug repurposing” (Kaneko and Nagashima, 2020), was adopted, making it viable to conduct a retrospective analysis with clinical big data and screen for hypothetical drug-drug pairs that potently ameliorate olanzapine-induced dyslipidemia in patients. We have previously utilized the self-reports of adverse events extracted from the FDA Adverse Events Reporting System (FAERS) to study the underlying mechanisms of hyperglycemia induced by an atypical antipsychotic, quetiapine, and identified the agent with treatment potential (Nagashima et al., 2016). Furthermore, to improve the drawbacks of FAERS (lacking time stamps and the population size of patients who were prescribed the specified drugs), another source of real-world clinical big data, insurance claims from JMDC Inc., was utilized in our recently published studies (Nagaoka et al., 2021; Siswanto et al., 2021). In the present study, in addition to the FAERS and JMDC databases, electronic medical records stored in the clinical data warehouse of Nihon University School of Medicine (NUSM’s CDW) containing detailed demographic, diagnostic, and laboratory data from patients at three hospitals affiliated with NUSM (Nagashima et al., 2022) were also used. Electronic medical records can provide continuous clinical information, leading to unattainably precise outcomes of the onset pattern of specific drug-induced adverse effects, treatment significance, and validated medication safety. To summarize, the corroborating results from the analysis of the three data sources mentioned above enable a more sensitive retrospective analysis for detecting causal relationships between the resultant adverse events and specific drug treatment, and the lowered incidence of adverse events and potential rescuing drugs causing this change. This hypothesis has been validated in animal experiments.
Materials and methods
Analysis of the FAERS database
Adversary event reports from 2004 to 2019 were downloaded from the FDA website (https://www.fda.gov/drugs/drug-approvals-anddatabases/fda-adverse-event-reporting-system-faers). Duplicate reports were deleted as previously reported (Banda et al., 2016), and 11,438,031 preserved reports were analyzed in the present study. Arbitrary drug names, including trade names and abbreviations, were manually mapped to unified generic names using the Medical Subject Headings (MeSH) descriptor ID. Reports of dyslipidemia were defined according to the preferred terms listed in Supplementary Table S1 of MedDRA (version 23.0). The FAERS data analysis was performed as previously described (Nagashima et al., 2016). Adversary event risk was evaluated by calculating the reporting odds ratio (ROR) along with the 95% confidence interval (CI) and Z-score (see Supplementary Tables S2, S3 for details). Z-scores were applied instead of p-values to save space graphically.
Analyzing the JMDC claims database
Insurance claims data collected by JMDC Inc. from January 2005 to March 2018 were purchased, which contained the medical and prescription claims of 5,550,241 individuals and their dependents on a monthly basis. Health checkup data, including blood test results, body mass index (BMI), and waist circumferences, were provided by 2,278,697 individuals. The patients were mainly aged ≤ 65 years, and no patients aged ≥ 75 years were included. The drug names were coded using the Anatomical Therapeutic Chemical Classification System.
On retrospectively analyzing the profile of the first event onsets, the users of olanzapine were identified as those who were prescribed olanzapine more than 2 months after being included in the JMDC claims database. Patients without specific health checkup data were also excluded. Thus, 1,853 patients with olanzapine were included in the present study. Among them, vitamin D users were defined as patients whose first prescription of vitamin D was ahead of olanzapine treatment. The pre-prescription period for olanzapine was defined as the period within 12 months, while the post-prescription period was defined as the period within 12 months after initiating olanzapine treatment.
The blood test results included in the present study were triglyceride, LDL cholesterol, and high-density lipoprotein (HDL) cholesterol levels apart from BMI and waist circumstances. Test results were collected for each individual on the nearest date before initiating olanzapine treatment in the pre-prescription period and on the closest date in the post-prescription period after initiating olanzapine treatment. To reduce bias in population backgrounds, the propensity score matching method (greedy 1:1 matching) was used by balancing covariates between settings. The propensity score for vitamin D was obtained by fitting a logistic regression model that included all covariates of interest (Table 2). The propensity score of each patient with or without vitamin D treatment was subsequently matched using the nearest neighbor method. Using the matched outcomes, an unpaired two-tailed t-test with Welch’s correction for continuous variables and Fisher’s exact test for categorical data were conducted to compare the differences in baseline characteristics between patients with or without vitamin D exposure.
Analysis of the NUSM electronic medical records
Electronic medical records were obtained from the Nihon University School of Medicine’s Clinical Data Warehouse (NUSM’s CDW). This database includes detailed diagnostic, demographic, and laboratory data for inpatients and outpatients. We received written informed consent and agreement for the secondary use after anonymization at three hospitals affiliated with the NUSM (Nagashima et al., 2022). Similar to the JMDC claims data analysis, the pre-prescription period of olanzapine was defined as the period within 12 months of the first prescription. Contrastingly, the post-prescription period was defined as the period within 12 months after initiating olanzapine treatment. Among them, preceding users of vitamin D were defined as patients whose first prescription of vitamin D was ahead of the first prescription of olanzapine for at least 1 day. Matching was not performed in this study because there were no significant differences in the baseline characteristics (Table 3). The blood laboratory values were extracted from the database, and the baseline value (month 0) was selected from the nearest test prior to the first prescription of olanzapine. The values at each checkpoint (months 3, 6, 9, and 12) were selected from the furthest test results from the first prescription of olanzapine. The missing values were imputed using the last observation carried forward method. An unpaired two-tailed t-test with Welch’s correction for continuous variables and Fisher’s exact test for categorical data were conducted to compare the differences in baseline characteristics between patients with or without vitamin D exposure.
Animals
All animal experiments were approved by the Kyoto University Animal Research Committee in accordance with the ethical guidelines of the committee. All experiments were designed to minimize the use of animals and the number of required experiments. Male and female C57BL/6J mice were purchased from Japan SLC (Shizuoka, Japan) and housed at a constant ambient temperature (24 ± 1°C) and humidity (55% ± 10%) on a 12/12 h light/dark cycle. Mice were fed an ad libitum diet consisting of water and chow.
Reagents and treatments (in vivo)
Olanzapine was purchased from the Tokyo Chemical Industry (Tokyo, Japan). The medium-fat diet (containing 1.37 IU of cholecalciferol/g) and cholecalciferol-supplemented medium-fat diet (containing 200 IU of cholecalciferol/g) were purchased from Oriental Yeast (Tokyo, Japan). Olanzapine was dissolved in water with 0.5% of carboxymethyl cellulose before use. Cholecalciferol is usually administered orally as a chow supplement, and the intake dose may vary slightly among mice. Mice were first randomized into two groups and fed a medium-fat diet or a cholecalciferol-supplemented medium-fat diet for 1 week. Mice from each group were further randomized to be treated with olanzapine (10 mg/kg, orally administered) or vehicle (0.5% of carboxymethyl cellulose solution) for another 5 days. One day after the last dose of olanzapine, the mice were anesthetized and dissected for blood collection by cardiac puncture. Mouse serum was isolated from blood samples by centrifugation. Serum total cholesterol levels were measured using a LabAssay Cholesterol kit (Wako, Osaka, Japan), while serum LDL cholesterol, HDL cholesterol, and triglyceride levels were measured by Nagahama Lifescience (Oriental Yeast Co., Ltd., Shiga, Japan).
Cell culture
Primary mouse hepatocytes were prepared as previously described (Charni-Natan and Goldstein, 2020; Jung et al., 2020) with modifications. Livers from 6 to 8 weeks old male mice were perfused with Hanks’ balanced salt solution (HBSS) without calcium, magnesium, and phenol red, and subsequently supplemented with 0.5 mM of ethylenediaminetetraacetic acid (EDTA) and 25 mM of hydroxyethyl piperazineethanesulfonic acid (HEPES). The livers were then perfused with a digestive enzyme mix [collagenase, type I (Worthington, NJ, United States) 0.15 mg/mL; collagenase, type II (Worthington, NJ, United States) 0.15 mg/mL; Dispase, type II (Gibco, MA, United States) 0.15 mg/mL] solution, dissolved in standard HBSS with phenol red, and supplemented with 25 mM of HEPES. Hepatocytes were released into high-glucose Dulbecco’s modified Eagle’s medium (DMEM; D5796; Sigma-Aldrich, MO, United States) from digested livers, filtered through 70 μm cell strainers (Corning, NY, United States), and further purified using a 45% Percoll® cushion by density gradient centrifugation. The purified cells were resuspended in William’s E medium (A1217601; Gibco) supplemented with 5% of fetal bovine serum (FBS; Sigma-Aldrich), 1 μM of dexamethasone (Nacalai Tesque, Kyoto, Japan), 1% of Penicillin–Streptomycin Mixed Solution (P/S; Nacalai Tesque), 5 μg/mL of human recombinant insulin (Sigma-Aldrich), 2 mM of GlutaMAX supplement (Sigma-Aldrich), and 15 mM of HEPES. Cells were plated onto dishes coated with 0.1 mg/mL of collagen type I (Nippi, Tokyo, Japan) and incubated at 37°C in a humidified chamber containing 5% CO2 for 3 h. After letting cells attach to the surface, the medium was refreshed with serum-free William’s E supplemented with 1 μM of dexamethasone, 0.5% of P/S, 1% of ITS + premix (Corning), 2 mM of GlutaMAX, and 15 mM of HEPES. Cells were used for treatment within 48 h.
Primary mouse adipocytes were prepared as previously described (Lequeux et al., 2009; Galmozzi et al., 2021) with modifications. Subcutaneous adipose tissues were collected from neonatal mouse pups and digested with 0.1 mg/mL collagenase type II and 0.1 mg/mL dispase type II dissolved in standard HBSS. Digested tissues were homogenized by pipetting, filtered through 70 μm cell strainers, and resuspended in high-glucose DMEM supplemented with 20% of FBS, 1% of P/S, and 10 mM of HEPES. The collected preadipocytes were plated onto dishes coated with 0.1% of gelatin (Nacalai Tesque) in distilled water. Cells were flushed twice with HBSS and refreshed with the same medium. After the cells reached 90% confluence, they were transferred into high-glucose DMEM supplemented with 10% of FBS, 1% of P/S, 10 mM of HEPES, 170 nM of human recombinant insulin, 1 μM of dexamethasone, 0.5 mM of 3-isobutyl-1-methylxanthine, 1 nM of triiodothyronine, and 10 nM of hydrocortisone for differentiation into adipocytes. After 48 h, the cells were transferred to high-glucose DMEM supplemented with 10% of FBS, 1% of P/S, 10 mM of HEPES, 170 nM of human recombinant insulin, and 1 nM of triiodothyronine. After 5 days of maintenance, the cells were ready for treatment.
Mouse C2C12 cells were prepared as previously described (Nagashima et al., 2016) with some modifications. The cells were obtained from Prof. H. Takeshima (Kyoto University Graduate School of Pharmaceutical Sciences, Kyoto, Japan). Cells were thawed in a water bath at 37°C and plated in high-glucose DMEM supplemented with 10% of FBS, 1% of P/S, and 10 mM of HEPES. After reaching 70% confluence, the cells were digested with 0.25% of trypsin solution (Nacalai Tesque) and passaged in the same medium. When the cells reached 90% confluence, the medium was changed to high-glucose DMEM supplemented with 2% of horse serum (Sigma-Aldrich), 1% of P/S, and 10 mM of HEPES for differentiation. The medium was refreshed 3 days after the changes, and the cells were ready on day five.
Reagents and treatments (in vitro)
Calcifediol was purchased from Selleckchem (TX, United States) and calcitriol was purchased from Cayman (MI, United States. ZK159222 was purchased from Cayman. All reagents were reconstructed and preserved in DMSO as a vehicle at −80°C, and thus DMSO was added as a blank reference in control groups. Primary mouse hepatocytes were treated with maintenance medium. Primary mouse adipocytes and C2C12 cells were treated in Advanced DMEM/F12 medium (12634010, Gibco) supplemented with 1% of P/S, 2 mM of GlutaMAX supplement, and 10 mM of HEPES.
RNA-seq and quantitative reverse transcriptase polymerase chain reaction (qRT-PCR)
RNA expression levels in cultured cells were evaluated using RNA-seq and qRT-PCR. Total RNA was isolated from cultured cells following the standard protocols provided by Isogen reagents (Nippon Gene, Tokyo, Japan).
For RNA-seq analysis, poly(A)+ RNA was selected from total RNA and sequenced using DNBseq (BGI, Shenzhen, China). Total reads were filtered by SOAPnuke (Cock et al., 2010) and clean reads were mapped to the mouse reference genome GRCm38.p6 using HISAT, version 2.0.4 (Kim et al., 2015) and Bowtie, version 2.2.5 (Langmead and Salzberg, 2012) in parallel. Gene expression was calculated using RSEM, version 1.2.8 (Li and Dewey, 2011). The gene set related to cholesterol biosynthesis was obtained from MGI (http://www.informatics.jax.org/vocab/gene_ontology/GO:0006695) (Hayamizu et al., 2005). All the genes from the three types of cultured cells (hepatocytes, adipocytes, and C2C12 cells) were mapped to this gene set to generate a list of cell-specific cholesterol biosynthetic process-related genes. Among the mapped genes, those with low expression levels (transcript per million (TPM) < 1) in vehicle-treated cells were removed from the analysis. Differentially expressed genes were detected between vehicle and olanzapine, olanzapine and olanzapine-calcifediol co-treated cells, and olanzapine and olanzapine-calcitriol co-treated cells based on PossionDis (Audic and Claverie, 1997). Genes with fold-change (FC) ≥ 2 and false discovery rate (FDR) ≤ 0.05 were recognized as differentially expressed. The analyzed results were visualized using volcano plots drawn using Prism 9.4.1 (GraphPad Software, CA, United States), with −log10 (FDR) in the y-direction and log2(FC) in the x-direction. The original sequence datasets were deposited in the NCBI sequence read archive with the accession number GSE221683.
For qRT-PCR, the extracted RNA was translated to cDNA using ReverTra Ace (Toyobo, Osaka, Japan) and subjected to StepOne real-time PCR (Life Technologies, Carlsbad, CA, United States) with Thunderbird SYBR qPCR Mix (Toyobo). The amplification process was set as follows: 10 min at 95°C, followed by 40 cycles of looping from 15 s at 95°C to 60 s at 60°C. Oligonucleotide primers were purchased from Sigma-Aldrich, and the sequences were as follows (gene name:5′-forward-3′, 5′-reverse-3′): Mus ribosomal protein, large, P0 (Rplp0): GCT TCG TGT TCA CCA AGG A, GTC CTA GAC CAG TGT TCT GAG G; Mus 3-hydroxy-3-methylglutaryl-CoA reductase (Hmgcr): GCT CGT CTA CAG AAA CTC CAC G, GCT TCA GCA GTG CTT TCT CCG T; Mus insulin-induced gene 2 (Insig2): GCC CAT CCA GAA CCT CTG AC, AGA CGG GGC AAA AGG ACT TC. The expression levels of each mRNA were normalized to those of Rplp0 mRNA.
Western blotting
Western blotting was performed as previously described (Ohashi et al., 2018) with some modifications. Cells were lysed in radioimmunoprecipitation assay (RIPA) buffer (Nacalai Tesque), diluted to 1 μg/mL with 1% sodium dodecyl sulfate (SDS) aqueous solution, and then buffered with NuPAGE® lithium dodecyl sulfate (LDS) sample buffer (4×; Life Technologies, CA, United States). The samples were loaded onto a 10% SDS-polyacrylamide gel and blotted onto a ClearTrans® PVDF Membrane (Wako). After blocking with Blocking One (Nacalai Tesque), the membranes were split into two pieces to simultaneously blot the target protein, vitamin D receptor (VDR), and the endo-reference protein, glyceraldehyde 3-phosphate dehydrogenase (GAPDH). The membranes were incubated overnight at 4°C with anti-VDR antibody (1:500 dilution, 12,550, Cell Signaling Technology, MA, United States) and anti-GAPDH (1:50,000 dilution, 32,233, Santa Cruz, CA, United States) in Tris-buffered saline supplemented with 0.1% of Tween 20 (TBS-T) and 10% of Blocking One. After washing with TBS-T, the membranes were incubated with peroxidase-conjugated donkey anti-rabbit IgG (1:5000 dilution, NA934V, GE Healthcare, IL, United States) and Peroxidase AffiniPure Goat Anti-Mouse IgG (1:5000 dilution, 115-035-003; Jackson ImmunoResearch, PA, United States) for 2 h at RT. Specific bands were detected using Immobilon Western Chemiluminescent HRP Substrate (Millipore) and visualized using EZCapture MG (ATTO, Tokyo, Japan). The expression level of VDR was normalized to that of GAPDH.
Statistics
Statistical analysis of the FAERS database was performed using SQL Software (IBM Research Laboratory, NY, United States) and that of the JMDC database was performed using R version 3.5.1 Software (The R Foundation for Statistical Computing, Vienna, Austria). In the JMDC database, paired t-tests were used to compare the mean values within the pre-prescription and post-prescription periods. Differences in continuous variables between the two groups were compared using an unpaired two-tailed t-test with Welch’s correction. In the NUSM database, mixed-effects models were used to analyze repeated measures data.
Statistical analysis of the animal experiments was performed using GraphPad Prism 9.3.1, and the data are expressed as mean ± standard error of the mean (SEM). The blood profile data were analyzed using two-way analysis of variance (ANOVA) with post hoc Tukey’s multiple comparison test. Revalidation using qRT-PCR of the C2C12 cell line on RNA-seq resultant data was analyzed by two-way ANOVA. The dose-dependent changes in specific genes responding to calcifediol, calcitriol, and ZK159222 in the C2C12 cell line were analyzed by one-way ANOVA. The western blotting results of VDR cell-specific enrichment were analyzed using Student’s t-test. All reported p-values of < 0.05 were considered statistically significant.
Results
Analysis of FAERS
First, the association between the treatment with a particular drug and the incidence of dyslipidemia in the FAERS database was investigated using disproportionality analysis by calculating the ROR and Z-score of each drug (Figure 1A). The known reporting bias and lack of incidence denominators accompanied by self-reports (Alatawi and Hansen, 2017) only allowed these values to be limited demonstrations of the real-world incidence rate. Nevertheless, the high ROR and Z-score of olanzapine indicated that it was one of the drugs that exhibited the strongest association between its use and the onset of dyslipidemia (the complete set of data is provided in Supplementary Table S2). Olanzapine was selected for the present study because of its higher ROR than that of other atypical antipsychotics (Table 1).
FIGURE 1. Increased incidence of dyslipidemia with the prescription of drugs and confounding effects of concomitant drugs on the olanzapine-induced dyslipidemia in FDA Adverse Event Reporting System (FAERS) data. Volcano plots for visualizing the reporting odds ratio (ROR, on a log scale) and its statistical significance (absolute Z-score) are shown. Each circle indicates an individual drug, and the size of the circle reflects the number of patients taking the drug. (A) Strong and significant increases in the ROR for dyslipidemia were seen in patients using olanzapine. (B) Within the population taking olanzapine, confounding effects of concomitantly used drugs on the incidence of olanzapine-induced dyslipidemia were calculated thoroughly and plotted. Overall values are presented in (Supplementary Tables S2, S3).
TABLE 1. Association between atypical antipsychotics and the incidence of dyslipidemia in the FAERS data.
The confounding effects of all drug combinations used in the olanzapine-treated patients were also calculated using the ROR and Z-score (Figure 1B). Vitamin D was found to be impressively effective in lowering the incidence of olanzapine-induced dyslipidemia. Although vitamin D was not among the drugs with the lowest ROR, which suggested higher effects on dyslipidemia, vitamin D had the highest Z-score due to its sufficient number of cases (the complete set of data is provided in Supplementary Table S3). However, vitamin D itself was associated with a slightly increased risk of dyslipidemia according to the FAERS data (ROR = 2.55, Z-score = 39.2; also see Supplementary Table S2), which may be due to the lack of time information to establish an accurate causal relationship between the reports and the effects of vitamin D. Thus, further analyses to reinforce the present findings were performed using the JMDC data.
Analysis of JMDC insurance claims
To further investigate whether the resultant dyslipidemia is associated with the clinical consequences of olanzapine treatment, the JMDC insurance claims data were introduced for analysis. Blood test data from the JMDC insurance claims of olanzapine-exposed 1,853 patients with or without vitamin D treatment were extracted (n = 22 for with vitamin D, n = 1,831 for without vitamin D). Since vitamin D is predominately prescribed to female elderly, a balancing approach, propensity score matching, was applied to eliminate possible systemic bias. After adjusting for populations, no significant differences in baseline characteristics between patients with or without vitamin D treatment (Table 2) were detected, and the batch size was balanced (n = 20). The effect of vitamin D on propensity score-matched data was investigated. In both groups of patients exposed only to olanzapine and to olanzapine and vitamin D simultaneously, the triglyceride level was not influenced (Figure 2A); however, the level of LDL cholesterol significantly increased by 11 mg/dL in the olanzapine-only group, while this value in patients with vitamin D co-treatment did not significantly increase (Figure 2B). Furthermore, the HDL cholesterol level in the blood significantly decreased by 5 mg/dL in the olanzapine-only group, while vitamin D supplementation significantly reversed the declining trend and even slightly increased HDL cholesterol levels by 4 mg/dL (Figure 2C). These results indicate that vitamin D counteracted the influence of olanzapine on blood lipid profiles.
TABLE 2. Population characteristics of the patients in the propensity score-matched olanzapine cohort extracted from the JMDC data; quantifiable data are given as means ± SEM.
FIGURE 2. Deteriorative effects of olanzapine on blood lipid profiles and body weight maintenance and the rescuing effects of vitamin D co-treatment from the JMDC insurance claims. The data of pre-prescription and post-prescription periods were collected from the nearest test data before or after the olanzapine prescription within 1 year. Background characteristics were matched between groups (n = 20) as shown in Table 2. (A–C) The annual blood test results from patients included in JMDC were extracted. (A) Serum concentrations of triglycerides, (B) low-density lipoprotein (LDL) cholesterol, and (C) high-density lipoprotein (HDL) cholesterol were expressed in units of mg/dL. (D,E) The annual body weights, heights, and waist circumferences included in the JMDC were extracted. (D) Body mass index (BMI) is calculated from body weights and heights and (E) waist circumferences were obtained by direct measurement. (F) the level of HbA1c expressed in units of percentage of blood. Data are shown as means ± SEM. *p < 0.05; a paired t-test of the pre-prescription period against the post-prescription period. #p < 0.05, ##p < 0.01; an unpaired two-tailed t-test with Welch’s correction between the group with and without vitamin D supplementation.
In addition to blood lipid profiles, the effects of vitamin D on body mass, fat distribution, and blood glucose levels in olanzapine-treated patients were investigated by analyzing related data from JMDC insurance claims. Within the same groups providing blood test data, BMI, an indicator used to roughly estimate whether one is overweight (as the higher BMI refers to a higher risk of being overweight), significantly increased by olanzapine treatment, while the co-treatment of vitamin D diminished this trend (Figure 2D). The mean value of the waist circumference, an approximate indicator of body fat mass, also significantly increased in olanzapine-only treated patients but not in olanzapine-vitamin D co-treated individuals (Figure 2E); however, the blood glucose indicator, HbA1c, was neither influenced in the olanzapine-only treated patients nor in the olanzapine-vitamin D co-treated individuals (Figure 2F). To summarize, these results indicate that using olanzapine aggravates the lipid profile, condition of being overweight, and excessive fat accumulation, while the combined use of vitamin D has beneficial effects on all three health aspects to some extent or at least provides no unfavorable effects. The blood and body test data in the JMDC insurance claims were annually collected from medical examinations provided by the Japan Health Insurance Association, which limited the precision of the data of the onset and development of olanzapine-induced dyslipidemia. Thus, to improve this flaw, electronic medical records were obtained from NUSM’s CDW, which is recorded in units of days, allowing tracing of blood profile changes on a smaller but more precise scale than that of the JMDC insurance claims data.
Analysis of electronic medical records from NUSM’s CDW
The blood test data of patients exposed to olanzapine with or without vitamin D treatment, including TG, LDL cholesterol, and HDL cholesterol values, were obtained from the electronic medical records from NUSM’s CDW. No significant bias was detected between patients with and without vitamin D supplementation (Table 3). By continuously tracking the blood lipid profile changes within 1 year after the first exposure to olanzapine, a significant increase in the blood triglyceride level was detected as early as 6 months after initiating the treatment and kept increasing by 54 mg/dL at the end of the year, while co-treatment with vitamin D significantly diminished this trend and remained stable (Figure 3A). A significant increase in blood LDL cholesterol in patients without vitamin D supplementation was detected not earlier than 12 months from the first exposure to olanzapine and increased by 15 mg/dL, but this trend was not observed in patients with vitamin D supplementation (Figure 3B). Contrastingly, HDL cholesterol in the blood of patients without vitamin D supplementation significantly reduced by 5 mg/dL after 6 months of follow-up, while vitamin D supplementation especially prevented this risk from developing (Figure 3C). Based on NUSM’s CDW electronic medical recordings, the beneficial effects of vitamin D treatments on olanzapine-induced dyslipidemia can also be detected within a shorter time scale in contrast to the insurance claims from the JMDC.
TABLE 3. Population characteristics of the patients in the olanzapine cohort of the NUSM’s CDW database; quantifiable data are given as means ± SEM.
FIGURE 3. Deteriorative effects of olanzapine on blood lipid profiles and the rescuing effects of vitamin D co-treatment from electronic health records of the clinical data warehouse of Nihon University School of Medicine. Blood test results were extracted every 3 months after the first prescription of olanzapine. (A) the triglycerides (n = 406 in olanzapine without vitamin D group; n = 23 in olanzapine with vitamin D group), (B) the LDL cholesterol levels (n = 154 in olanzapine without vitamin D group; n = 11 in olanzapine with vitamin D group), and (C) the HDL cholesterol levels (n = 258 in olanzapine without vitamin D group; n = 16 in olanzapine with vitamin D group) were expressed in units of mg/dL. Data are shown as means ± SEM. *p < 0.05, **p < 0.01, ***p < 0.001; a paired t-test of data on any months after the olanzapine prescription against the data on month 0 (baseline). #p < 0.05, ##p < 0.01; an unpaired two-tailed t-test with Welch’s correction between the groups with and without vitamin D supplementation.
Effects of vitamin D3 on olanzapine-induced dyslipidemia mouse model
Previous studies have reported non-alcoholic fatty liver disease and LDL cholesterol or HDL cholesterol imbalance in the blood phenotypes of olanzapine-treated wild-type mouse models (Chen et al., 2018; Jiang et al., 2019; Liu et al., 2019; Zhu et al., 2022) under chronic treatment (at least 4 weeks). In the present study, in order to study the primary stage of dyslipidemia development, the treatment was limited to 5 days (Figure 4A). To mimic the calorie intake of human patients, mice were provided with 30% of fructose in drinking water during the acclimatization period of 1 week and this was continued in the following runs of treatment (Faeh et al., 2005). The results showed that neither a scheduled administration of olanzapine (10 mg/kg, orally administered by gavage, daily for 5 days) nor vitamin D supplementation (administrated via a cholecalciferol (vitamin D3)-supplemented diet containing 200 IU of cholecalciferol/g) caused changes in blood triglyceride (Figure 4B) and total cholesterol levels (Figure 4C). However, a significant increase in the blood LDL cholesterol level and the attenuating effects of vitamin D supplementation were detected (Figure 4D). In addition, the reduction in HDL cholesterol level caused by olanzapine treatment was attenuated by cholecalciferol supplementation (Figure 4E). These findings suggest that short-term administration of olanzapine is sufficient to induce dyslipidemia by challenging the balance between blood LDL and HDL cholesterol levels without influencing blood triglyceride and total cholesterol levels. Simultaneously, cholecalciferol as dietary supplementation diminished these effects.
FIGURE 4. Olanzapine treatment caused dyslipidemia in a rodent model while vitamin D3 supplements attenuated the worsening of the condition. (A) C57BL/6J mice were fed for 1 week with the medium-fat diet (containing 1.37 IU of cholecalciferol/g of chow) or cholecalciferol-supplemented medium-fat diet (containing 200 IU of cholecalciferol/g of chow). Each group of mice were divided into two groups and orally given 10 mg/kg/day of olanzapine or vehicle for 5 days; resulting 4 groups of mice (each n = 10) received normal chow and vehicle as control (C), normal chow and olanzapine (O), vitamin D3-supplemented chow and vehicle (D), and vitamin D3-supplemented chow and olanzapine (DO). Blood samples were collected from these mice 1 day after the last dose of olanzapine without fasting. (B) the triglycerides, (C) total cholesterol levels, (D) LDL cholesterol levels, and (E) HDL cholesterol levels were measured. Data are shown as means ± SEM. The multiple comparisons of olanzapine and cholecalciferol influences were compared by two-way analysis of variance (ANOVA) with post hoc Tukey’s multiple comparison test. *p < 0.05, **p < 0.01.
Molecular mechanisms of vitamin D on olanzapine-induced dyslipidemia
To investigate the molecular mechanisms underlying cholesterol dyshomeostasis in blood accompanied by sub-chronic olanzapine and the rescuing effects of cholecalciferol supplementation, changes in expression of genes in cells related to cholesterol metabolism were analyzed using RNA-seq. Purified primary cultures of mouse hepatocytes, adipocytes, and fully differentiated C2C12 cells (myotubes) were used instead of liver, adipose, and skeletal muscle tissues, since these tissues are generally the mixtures of various cells with different origins and functions, which may blur the direct effects of olanzapine and vitamin D metabolites. All three types of cells were treated under the same condition, a 24 h pre-exposure to the circulating form of vitamin D3, 25-hydroxycholecarciferol, (calcifediol at 10 μM) or the active form of vitamin D3, 1,25-dihydroxycholecarciferol, (calcitriol at 0.1 μM), followed by another 12 h co-treatment with olanzapine (1 μM). From the MGI database, only cholesterol biosynthesis-related genes were selected for analysis (Hayamizu et al., 2005) because only cholesterol homeostasis was disrupted in the blood according to in vivo modeling. Comparing these gene expressions under different treatments in various types of cells with PossionDis (Audic and Claverie, 1997), the expression of a set of genes was significantly changed in olanzapine-calcifediol-treated C2C12 cells and slightly influenced in olanzapine-calcitriol-treated C2C12 cells. However, olanzapine displayed passive effects in all three cell types (Figure 5A). The regulatory relationship within this set of genes is illustrated in Figure 5B. In this scheme, Insig2, a suppressive regulator of cholesterol biosynthesis (Yabe et al., 2002), was upregulated by vitamin D metabolites, while genes participating in the mevalonate pathway (DeBose-Boyd, 2008), such as 3-hydroxy-3-methylglutaryl-CoA reductase (Hmgcr), 3-hydroxy-3-methylglutaryl-coenzyme A synthase 1 (Hmgcs1), farnesyl diphosphate synthetase (Fdps), farnesyl diphosphate farnesyl transferase 1 (Fdft1), lanosterol synthase (Lss), cytochrome P450, family 51 (Cyp51), NAD(P) dependent steroid dehydrogenase-like (Nsdhl), and 24-dehydrocholesterol reductase (Dhcr24) were significantly downregulated.
FIGURE 5. Volcano plots of gene expression in the cultured mice cells treated with vehicle, olanzapine, or olanzapine co-treated with calcifediol or calcitriol. The changing ranges were expressed as log2-transformed fold-change (FC), and the significance was expressed as −log10-transformed false discovery rate (FDR). (A) Demonstrations of the cholesterol biosynthesis-related gene changes under olanzapine treatment and calcifediol or calcitriol-olanzapine co-treatment. The genes significantly downregulated in the olanzapine-calcifediol co-treated C2C12 cells are tagged with blue dots, and the genes upregulated are tagged with red dots in all the data sets. (B) A schematic diagram illustrating the functions of downregulated genes (blue) and upregulated gene (red) in the process of cholesterol biosynthesis. Genes plotted with FC absolute value > 2 and FDR value > 0.05 were considered to be significantly changed by treatment.
To validate these findings from the RNA-seq analysis, an enlarged batch of evaluations was performed on C2C12 cells, which displayed the highest sensitivity to vitamin D metabolites, using qRT-PCR. Hmgcr was selected as a representative gene for cholesterol biosynthesis because it is the most representative gene involved in the process of cholesterol biosynthesis (Horton et al., 2003). The qRT-PCR results revealed that both the suppression of Hmgcr expression (Figure 6A) and the induction of the expression of Insig2 (Figure 6B) by calcifediol and calcitriol were significant. A previous report claimed that calcifediol, but not calcitriol, independently reduced cholesterol biosynthesis without VDR actions (Asano et al., 2017). However, further validations by changing the concentration of either calcifediol or calcitriol in the treatment medium showed that the expression of Hmgcr (Figure 6C) and Insig2 (Figure 6D) changed in a dose-dependent manner under calcifediol treatment. In parallel, the changes in the expression of Hmgcr (Figure 6E) and Insig2 (Figure 6F) also displayed a dose-dependent effect in calcitriol treatment, with approximately 5.8 times higher sensitivity (IC50 of calcifediol on Hmgcr: 0.68 µM; IC50 of calcitriol on Hmgcr: 0.12 µM). Since calcitriol is the bioactive form of vitamin D, which acts only through VDR, the role of VDR in cholesterol synthesis suppression in this model should not be underestimated and requires further investigation.
FIGURE 6. Validation by quantitative RT-PCR results showing the expression of cholesterol biosynthesis-associated genes in the cultured C2C12 cell line. (A) 3-hydroxy-3-methylglutaryl-CoA reductase (Hmgcr) and (B) Insulin-induced gene 2 (Insig2) were validated in the same treating condition as that in the RNA-seq experiments (n = 12). Results were normalized to ribosomal protein, large, P0 (Rplp0). (C) Hmgcr and (D) Insig2 were validated in calcifediol with gradient-changed concentrations increasing from 0 to 10 μM (n = 8–9). (E) Hmgcr and (F) Insig2 were validated in calcitriol with gradient-changed concentrations increasing from 0.01 μM to 1 μM. Data are shown as means ± SEM (n = 9). **p < 0.01, ****p < 0.0001. The comparisons of olanzapine and vitamin D influences were compared by two-way analysis of variance (ANOVA). The validation of the dose-dependent effects of calcifediol and calcitriol were compared by one-way ANOVA.
To investigate whether the actions of calcifediol and calcitriol were mediated by VDR activity, western blotting of VDR levels in specific cell types was performed using RNA-seq. VDR was highly expressed in C2C12 cells, which was much higher than that in adipocytes and hepatocytes (Figure 7A). The varied tissue-specific enrichment of VDR indicates that suppression of Hmgcr expression is mediated by VDR activity. A partial agonist of VDR (Teske et al., 2016), ZK159222, was used to validate the action of VDR. ZK159222, like calcifediol and calcitriol, suppressed Hmgcr expression (Figure 7B) and induced Insig2 expression (Figure 7C) in a dose-dependent manner, further confirming the essential role of VDR activity.
FIGURE 7. Validation by quantitative RT-PCR results showing calcitriol regulating the cholesterol biosynthesis-associated genes mediated by VDR in the C2C12 cell line. (A) The results of vitamin D receptor (VDR) protein levels in C2C12 cells, adipocytes, and hepatocytes, normalized to glyceraldehyde 3-phosphate dehydrogenase (GAPDH) (n = 3). (B) Hmgcr and (C) Insig2 were validated in ZK159222 with gradient-changed concentrations increasing from 0 to 3 μM (n = 6). Data are shown as means ± SEM. *p < 0.05, **p < 0.01. The validation of the dose-dependent effects of ZK150222 was compared by one-way analysis of variance (ANOVA).
Discussion
In this study, we showed for the first time that vitamin D supplementation is effective in improving the lipid profiles in patients treated with olanzapine. This finding was revealed from three independent sources of clinical data. This hypothesis was validated by using the blood profiles of mice models fed with high-calorie diets, and the mechanism of action was investigated at the molecular level using cells that were metabolically active with cholesterol.
Drug repurposing using FAERS has been proven to be a reliable tool for refining unexpected drug-drug interactions with beneficial effects in treating adverse effects (Nagaoka et al., 2021; Siswanto et al., 2021). In the present study, olanzapine was one of the most dyslipidemia-inducing drugs causing a large ROR in more than tens of thousands of cases, according to the FAERS database, providing a solid foundation to further explore rescue drugs. Additionally, olanzapine had a greater tendency to induce severe dyslipidemia than other atypical antipsychotics that share similar molecular structures, such as clozapine and quetiapine, in the FAERS database (Table 1). Notably, patients are more heavily dosed with clozapine (Cmax > 1.07 μM) (Jann et al., 1993) and quetiapine (Cmax ≈ 0.21 μM) (DeVane and Nemeroff, 2001) than olanzapine (Cmax < 0.064 μM) (Callaghan et al., 1999) in clinical practice. These two facts prioritize screening drug targets for treating olanzapine-induced dyslipidemia from other atypical antipsychotics. Contrastingly, vitamin D treatment is recognized as a dyslipidemia risk factor other than olanzapine, which is likely caused by the confounding in identifying causes and results. Since patients with dyslipidemia are treated with this drug, the onset of the disease as per the FAERS data, is also likely to be positively related to the usage of this drug. This confusion is due to the lack of information regarding the prescription period to help define the time point of drug exposure and disease onset. However, cross-validations using health databases with records on prescription time can help solve this issue.
The hypothesis developed from the FAERS data analysis was verified in a retrospective cohort analysis based on the JMDC insurance claims data. The onset rate of dyslipidemia and the drug exposure period can be analyzed in units of years using the raw blood test results of patients under 75 years of age that the JMDC claims can annually provide as a supplement to the FAERS data. In the present study, propensity score matching revealed olanzapine-induced dyslipidemia and weight gain in the matched groups. However, the patient population characteristics are strongly sex-biased (primarily women), because vitamin D is predominantly prescribed to women in cases where they are more susceptible to osteoporosis (Bohon and Goolsby, 2013). Nevertheless, olanzapine has been reported to increase LDL cholesterol and triglyceride levels in adolescents (Kryzhanovskaya et al., 2009) and young adults (Narula et al., 2010), regardless of gender variance, indicating that the analysis results from JMDC claims still remain highly consistent with clinical trials. The present study found that vitamin D supplementation during the exposure period to olanzapine suppressed the increase in the LDL cholesterol level, decrease in the HDL cholesterol level, increase in BMI, and increase in waist circumference. This is the first study to demonstrate the clinical significance of vitamin D in improving blood lipid profiles, since only one previous trial has shown that vitamin D helps in maintaining waist circumference in olanzapine-treated patients (Kaviyani et al., 2017). Despite already being efficient in discovering the olanzapine-vitamin D interaction at the laboratory test level, the JMDC is still unable to precisely locate the onset time pattern of olanzapine-induced dyslipidemia due to the long data collection intervals. Further improvement was achieved by utilizing the electronic medical records from NUSM’s CDW, since the electronic records were updated daily and allowed for a flexible time-course setting. The blood test results were collected every 3 months after initiating the olanzapine treatment. The triglyceride and LDL cholesterol levels were reported to increase monotonically as the exposure time increased, and the appearance of abnormality in triglyceride levels was prior to that in the LDL and HDL cholesterol levels. Differentiating from JMDC claims, patients treated with olanzapine only exhibited a strong trend in developing triglyceride abnormality, while triglyceride levels from patients in JMDC claims were passively influenced. This inconsistency may be due to the variance in population background characteristics, since patients included in the electronic records from NUSM’s CDW patients were more balanced in terms of sex (female percentage = 61.3%) than those from the JMDC claims (female rate = 90.0%), and higher estrogen levels in females are generally believed to be the reason for better triglyceride control than males (Carr, 2003). Vitamin D treatment counteracted these deteriorative changes by maintaining the blood lipid profile at stable levels throughout the observation window, consistent with the outcomes of the JMDC claims. By combining these three outcomes based on different data sources, a causal relationship between lipid profiles’ deterioration and olanzapine treatment was detected in the present study, and the suppressive effects of vitamin D supplementation on these influences were also confirmed.
Using drug-induced animal models, significant efforts have been devoted to identifying the mechanism of olanzapine-induced dyslipidemia. In several rodent models, strong hepatoxicity, represented by non-alcoholic fatty liver diseases, has generally been reported after chronic treatments ranging from 4 to 12 weeks (Chen et al., 2018; Jiang et al., 2019; Liu et al., 2019; Zhu et al., 2022). However, these models cannot accurately represent clinical observations from patients, since olanzapine rarely causes clinically apparent hepatotoxicity (Larrey and Ripault, 2013). Furthermore, it generally takes years for non-alcoholic fatty liver disease to develop under constant environmental stress (Friedman et al., 2018). Contrastingly, according to the electronic medical records from NUSM’s CDW in the present study, the onset of dyslipidemia was as early as 6 months after initiating olanzapine treatment, which suggests that the development of olanzapine-induced dyslipidemia is not necessarily related to the onset of hepatic diseases. This assumption is supported by the in vivo observation in mouse models in the present study that olanzapine exposure as short as 1 week is already enough to cause the LDL and HDL cholesterol levels to increase and drop, respectively. A recent clinical study revealed an alternative explanation to olanzapine-induced dyslipidemia that an increased appetite is essential for developing olanzapine-induced dyslipidemia (Huang et al., 2020), and moreover, enhanced food reward circuitry has been detected in olanzapine-treated patients (Mathews et al., 2012). This phenotype was previously observed in rodent models and attributed to the antagonistic effects on histamine H1 receptors in the central nervous system by olanzapine (Kim et al., 2007; He et al., 2014). Nevertheless, a slight increase in the expression of genes related to cholesterol biosynthesis was detected in the hepatocytes in the present study. However, the changing ranges were too limited to be considered investigation-worthy. In conclusion, hepatotoxicity is not causally related to olanzapine-induced dyslipidemia and the direct influence of olanzapine on lipid metabolic processes in related tissues lacks physiological significance.
However, the role of vitamin D in lipid metabolism remains controversial. Several studies have shown that vitamin D is necessary for lipid synthesis and accumulation in adipose tissue. In VDR-knockout mouse models, the absence of VDR showed substantial protective effects against high-fat diet-induced adipose tissue accumulation (Narvaez et al., 2009; Wong et al., 2009), whereas in another knock-in model, mice with overexpressed human VDR in adipose tissue developed obesity even under a standard diet (Wong et al., 2011). Further investigation into this mechanism revealed that vitamin D potentiates subcutaneous preadipocyte differentiation via the VDR pathway. Simultaneously, the expression level of VDR significantly decreased as differentiation progressed until the full maturity of adipocytes (Nimitphong et al., 2012). Contrastingly, vitamin D alleviated atherosclerosis progression in patients by preventing foam cell formation (Oh et al., 2009) and suppressing oxidative stress in blood vessel endothelial cells (Kassi et al., 2013). However, the role of vitamin D in lipid homeostasis is not necessarily associated with maintaining the blood lipid profile. Clinical observations have reported improved blood triglyceride and LDL cholesterol levels in populations supplied with extra vitamin D, with varied background characteristics, according to two meta-analyses (Mirhosseini et al., 2018; Dibaba, 2019). The big data mining-based results in the present study are consistent with these findings. However, according to a previous study on the effects of calcium supplements on the blood lipid profile, the primary medical function of vitamin D, which improves calcium absorption, is not likely to be related to its role in treating dyslipidemia (Reid et al., 2010). To date, studies revealing the mechanisms of these clinical results are lacking. Nevertheless, a recent study established an inhibitory connection between vitamin D and lipid synthesis by proving that calcifediol, but not calcitriol, inhibits cholesterol and fatty acid biosynthesis by preventing sterol regulatory element binding transcription factors from maturing in a VDR-independent pattern (Asano et al., 2017). However, this result was not reproduced in primary hepatocytes and adipocytes, which may be due to the vast metabolic variance between highly differentiated primary cells and immortalized cells in that report (CHO cell line). Another study suggested that VDR-mediated Insig2 expression by calcitriol in HepG2 cells mainly contributed to its suppressive effects on cholesterol synthesis (Li et al., 2016). This result is doubtful because VDR is absent in normal hepatocytes (Gascon-Barre et al., 2003). This understanding is consistent with the lack of VDR expression in primary cultured hepatocytes in the western blot results of the present study. Nevertheless, the VDR activation-induced Insig2 expression increase is reproduced in this study. Insig2 is an essential regulator of cholesterol biosynthesis that not only prevents sterol regulatory element binding transcription factors from entering the cell nucleus to suppress their transcription function (Yabe et al., 2002) but also facilitates HMGCR degradation (Jo et al., 2011). A functional vitamin D response element has been localized in the promoter of the mouse Insig2 (Lee et al., 2005), and further comparisons between rodent and human Insig2 promoters have confirmed that this vitamin D response element is also included in the human genome (Fernandez-Alvarez et al., 2010). Although in contrast to the liver, the role of the muscle tissues is minor in the essentiality of cholesterol production and balance, cholesterol production in peripheral tissues should not be overlooked. A report has shown myopathy development in skeletal muscle-specific Hmgcr knockout mice (Osaki et al., 2015). Additionally, humans with skeletal muscle Hmgcr deficiency easily develop myopathy (Lorenzo-Villalba et al., 2021). Furthermore, as VDR is widely expressed in tissues throughout the body, vitamin D may improve the blood cholesterol condition by limiting tissue self-biosynthesis of cholesterol and simultaneously inducing cholesterol absorption in peripheral tissues. A primary limitation of this research is the failure to clarify the mechanism of olanzapine-induced dyslipidemia through investigations on lipid metabolism. Future explorations will focus on the neurological role of olanzapine in increasing appetites. Additionally, the investigation on vitamin D target tissues and organs is also limited, since the enterocytes in the intestine are also essential for lipid absorption and highly enriched with VDR. The possible role of enterocytes is also scheduled for future studies.
Although the exploration of the potential of VDR-mediated cholesterol biosynthesis suppression in new druggable target development still has a long way to cover, from the clinical findings in this research, it is safe to suggest vitamin D supplements in patients under olanzapine treatment, especially in those diagnosed as vitamin D deficient. Moreover, this treatment is also recommended to be simultaneously applied with a calorie-limited diet for patients under olanzapine treatment for a better outcome.
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 Gene Expression Omnibus (GEO) [https://www.ncbi.nlm.nih.gov/geo/], GSE221683.
Ethics statement
The animal study protocol was reviewed and approved by the Kyoto University Animal Research Committee.
Author contributions
TN and SK designed the study. TN, CT, and TS performed clinical data analysis. ZZ and MK conducted the experiments and analyzed the data. ZZ and SK wrote the manuscript. KN, HS, SA, and SK provided materials and technical advice.
Funding
This work was supported by the Japan Agency for Medical Research and Development (Grant Nos. 20ak0101088h0003 and 21ak0101153h0001 to SK) and the Japanese Society for the Promotion of Science Grants-in-Aid for Scientific Research to SK (20H00491).
Conflict of interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as potential conflicts 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.
Supplementary material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fphar.2023.1135516/full#supplementary-material
References
Alatawi, Y. M., and Hansen, R. A. (2017). Empirical estimation of under-reporting in the U.S. Food and drug administration adverse event reporting system (FAERS). Expert Opin. Drug Saf. 16, 761–767. doi:10.1080/14740338.2017.1323867
Asano, L., Watanabe, M., Ryoden, Y., Usuda, K., Yamaguchi, T., Khambu, B., et al. (2017). Vitamin D metabolite, 25-hydroxyvitamin D, regulates lipid metabolism by inducing degradation of SREBP/SCAP. Cell Chem. Biol. 24, 207–217. doi:10.1016/j.chembiol.2016.12.017
Audic, S., and Claverie, J. M. (1997). The significance of digital gene expression profiles. Genome Res. 7, 986–995. doi:10.1101/gr.7.10.986
Banda, J. M., Evans, L., Vanguri, R. S., Tatonetti, N. P., Ryan, P. B., and Shah, N. H. (2016). A curated and standardized adverse drug event resource to accelerate drug safety research. Sci. Data 3, 160026. doi:10.1038/sdata.2016.26
Bohon, T. M., and Goolsby, M. A. (2013). The role of vitamin D supplements in women's health. Clin. Med. Insights Womens Health 6, 67–70. doi:10.4137/CMWH.S11067
Brown, S., Inskip, H., and Barraclough, B. (2000). Causes of the excess mortality of schizophrenia. Br. J. Psychiatry 177, 212–217. doi:10.1192/bjp.177.3.212
Callaghan, J. T., Bergstrom, R. F., Ptak, L. R., and Beasley, C. M. (1999). Olanzapine. pharmacokinetic and pharmacodynamic profile. Clin. Pharmacokinet. 37, 177–193. doi:10.2165/00003088-199937030-00001
Campbell, M., Young, P. I., Bateman, D. N., Smith, J. M., and Thomas, S. H. (1999). The use of atypical antipsychotics in the management of schizophrenia. Br. J. Clin. Pharmacol. 47, 13–22. doi:10.1046/j.1365-2125.1999.00849.x
Carr, M. C. (2003). The emergence of the metabolic syndrome with menopause. J. Clin. Endocrinol. Metab. 88, 2404–2411. doi:10.1210/jc.2003-030242
Charni-Natan, M., and Goldstein, I. (2020). Protocol for primary mouse hepatocyte isolation. Star. Protoc. 1, 100086. doi:10.1016/j.xpro.2020.100086
Chen, C. H., Chiu, C. C., Huang, M. C., Wu, T. H., Liu, H. C., and Lu, M. L. (2008). Metformin for metabolic dysregulation in schizophrenic patients treated with olanzapine. Prog. Neuropsychopharmacol. Biol. Psychiatry 32, 925–931. doi:10.1016/j.pnpbp.2007.11.013
Chen, C. H., Shyue, S. K., Hsu, C. P., and Lee, T. S. (2018). Atypical antipsychotic drug olanzapine deregulates hepatic lipid metabolism and aortic inflammation and aggravates atherosclerosis. Cell Physiol. Biochem. 50, 1216–1229. doi:10.1159/000494573
Cock, P. J., Fields, C. J., Goto, N., Heuer, M. L., and Rice, P. M. (2010). The Sanger FASTQ file format for sequences with quality scores, and the Solexa/Illumina FASTQ variants. Nucleic Acids Res. 38, 1767–1771. doi:10.1093/nar/gkp1137
Correll, C. U., Joffe, B. I., Rosen, L. M., Sullivan, T. B., and Joffe, R. T. (2015). Cardiovascular and cerebrovascular risk factors and events associated with second-generation antipsychotic compared to antidepressant use in a non-elderly adult sample: Results from a claims-based inception cohort study. World Psychiatry 14, 56–63. doi:10.1002/wps.20187
DeBose-Boyd, R. A. (2008). Feedback regulation of cholesterol synthesis: Sterol-accelerated ubiquitination and degradation of HMG CoA reductase. Cell Res. 18, 609–621. doi:10.1038/cr.2008.61
DeVane, C. L., and Nemeroff, C. B. (2001). Clinical pharmacokinetics of quetiapine: An atypical antipsychotic. Clin. Pharmacokinet. 40, 509–522. doi:10.2165/00003088-200140070-00003
Dibaba, D. T. (2019). Effect of vitamin D supplementation on serum lipid profiles: A systematic review and meta-analysis. Nutr. Rev. 77, 890–902. doi:10.1093/nutrit/nuz037
Faeh, D., Minehira, K., Schwarz, J. M., Periasamy, R., Park, S., and Tappy, L. (2005). Effect of fructose overfeeding and fish oil administration on hepatic de novo lipogenesis and insulin sensitivity in healthy men. Diabetes 54, 1907–1913. doi:10.2337/diabetes.54.7.1907
Fernandez-Alvarez, A., Soledad Alvarez, M., Cucarella, C., and Casado, M. (2010). Characterization of the human insulin-induced gene 2 (INSIG2) promoter: The role of Ets-binding motifs. J. Biol. Chem. 285, 11765–11774. doi:10.1074/jbc.M109.067447
Friedman, S. L., Neuschwander-Tetri, B. A., Rinella, M., and Sanyal, A. J. (2018). Mechanisms of NAFLD development and therapeutic strategies. Nat. Med. 24, 908–922. doi:10.1038/s41591-018-0104-9
Galmozzi, A., Kok, B. P., and Saez, E. (2021). Isolation and differentiation of primary white and Brown preadipocytes from newborn Mice. J. Vis. Exp. 167, e62005. doi:10.3791/62005
Gascon-Barre, M., Demers, C., Mirshahi, A., Neron, S., Zalzal, S., and Nanci, A. (2003). The normal liver harbors the vitamin D nuclear receptor in nonparenchymal and biliary epithelial cells. Hepatology 37, 1034–1042. doi:10.1053/jhep.2003.50176
Hayamizu, T. F., Mangan, M., Corradi, J. P., Kadin, J. A., and Ringwald, M. (2005). The adult mouse anatomical dictionary: A tool for annotating and integrating data. Genome Biol. 6, R29. doi:10.1186/gb-2005-6-3-r29
He, M., Zhang, Q., Deng, C., Wang, H., Lian, J., and Huang, X. F. (2014). Hypothalamic histamine H1 receptor-AMPK signaling time-dependently mediates olanzapine-induced hyperphagia and weight gain in female rats. Psychoneuroendocrinology 42, 153–164. doi:10.1016/j.psyneuen.2014.01.018
Horton, J. D., Shah, N. A., Warrington, J. A., Anderson, N. N., Park, S. W., Brown, M. S., et al. (2003). Combined analysis of oligonucleotide microarray data from transgenic and knockout mice identifies direct SREBP target genes. Proc. Natl. Acad. Sci. U. S. A. 100, 12027–12032. doi:10.1073/pnas.1534923100
Huang, J., Hei, G.-R., Yang, Y., Liu, C.-C., Xiao, J.-M., Long, Y.-J., et al. (2020). Increased appetite plays a key role in olanzapine-induced weight gain in first-episode schizophrenia patients. Front. Pharmacol. 11, 739. doi:10.3389/fphar.2020.00739
Jann, M. W., Grimsley, S. R., Gray, E. C., and Chang, W. H. (1993). Pharmacokinetics and pharmacodynamics of clozapine. Clin. Pharmacokinet. 24, 161–176. doi:10.2165/00003088-199324020-00005
Jiang, T., Zhang, Y., Bai, M., Li, P., Wang, W., Chen, M., et al. (2019). Up-regulation of hepatic fatty acid transporters and inhibition/down-regulation of hepatic OCTN2 contribute to olanzapine-induced liver steatosis. Toxicol. Lett. 316, 183–193. doi:10.1016/j.toxlet.2019.08.013
Jo, Y., Lee, P. C., Sguigna, P. V., and DeBose-Boyd, R. A. (2011). Sterol-induced degradation of HMG CoA reductase depends on interplay of two Insigs and two ubiquitin ligases, gp78 and Trc8. Proc. Natl. Acad. Sci. U. S. A. 108, 20503–20508. doi:10.1073/pnas.1112831108
Jung, Y., Zhao, M., and Svensson, K. J. (2020). Isolation, culture, and functional analysis of hepatocytes from mice with fatty liver disease. Star. Protoc. 1, 100222. doi:10.1016/j.xpro.2020.100222
Kanagasundaram, P., Lee, J., Prasad, F., Costa-Dookhan, K. A., Hamel, L., Gordon, M., et al. (2021). Pharmacological interventions to treat antipsychotic-induced dyslipidemia in schizophrenia patients: A systematic review and meta analysis. Front. Psychiatry 12, 642403. doi:10.3389/fpsyt.2021.642403
Kaneko, S., and Nagashima, T. (2020). Drug repositioning and target finding based on clinical evidence. Biol. Pharm. Bull. 43, 362–365. doi:10.1248/bpb.b19-00929
Kassi, E., Adamopoulos, C., Basdra, E. K., and Papavassiliou, A. G. (2013). Role of vitamin D in atherosclerosis. Circulation 128, 2517–2531. doi:10.1161/CIRCULATIONAHA.113.002654
Kaviyani, S., Bahadoram, M., Houshmand, G., and Bahadoram, S. (2017). Therapeutic impact of cholecalciferol in patients with psychiatric disorders receiving olanzapine. J. Parathyr. Dis. 6, 19–22. doi:10.15171/jpd.2018.07
Kennedy, J. S., Bymaster, F. P., Schuh, L., Calligaro, D. O., Nomikos, G., Felder, C. C., et al. (2001). A current review of olanzapine's safety in the geriatric patient: From pre-clinical pharmacology to clinical data. Int. J. Geriatr. Psychiatry 16, S33–S61. doi:10.1002/1099-1166(200112)16:1+<::aid-gps571>3.0.co;2-5
Kim, S. F., Huang, A. S., Snowman, A. M., Teuscher, C., and Snyder, S. H. (2007). From the cover: Antipsychotic drug-induced weight gain mediated by histamine H1 receptor-linked activation of hypothalamic AMP-kinase. Proc. Natl. Acad. Sci. U. S. A. 104, 3456–3459. doi:10.1073/pnas.0611417104
Kim, D., Langmead, B., and Salzberg, S. L. (2015). HISAT: A fast spliced aligner with low memory requirements. Nat. Methods 12, 357–360. doi:10.1038/nmeth.3317
Kochi, K., Sato, I., Nishiyama, C., Tanaka-Mizuno, S., Doi, Y., Arai, M., et al. (2017). Trends in antipsychotic prescriptions for Japanese outpatients during 2006-2012: A descriptive epidemiological study. Pharmacoepidemiol. Drug Saf. 26, 642–656. doi:10.1002/pds.4187
Kryzhanovskaya, L., Schulz, S. C., McDougle, C., Frazier, J., Dittmann, R., Robertson-Plouch, C., et al. (2009). Olanzapine versus placebo in adolescents with schizophrenia: A 6-week, randomized, double-blind, placebo-controlled trial. J. Am. Acad. Child. Adolesc. Psychiatry 48, 60–70. doi:10.1097/CHI.0b013e3181900404
Langmead, B., and Salzberg, S. L. (2012). Fast gapped-read alignment with Bowtie 2. Nat. Methods 9, 357–359. doi:10.1038/nmeth.1923
Larrey, D., and Ripault, M.-P. (2013). “Hepatotoxicity of psychotropic drugs and drugs of abuse,” in Drug-induced liver disease, 443–462.
Larsen, J. R., Vedtofte, L., Jakobsen, M. S. L., Jespersen, H. R., Jakobsen, M. I., Svensson, C. K., et al. (2017). Effect of liraglutide treatment on prediabetes and overweight or obesity in clozapine- or olanzapine-treated patients with schizophrenia spectrum disorder: A randomized clinical trial. JAMA Psychiatry 74, 719–728. doi:10.1001/jamapsychiatry.2017.1220
Lee, S., Lee, D. K., Choi, E., and Lee, J. W. (2005). Identification of a functional vitamin D response element in the murine Insig-2 promoter and its potential role in the differentiation of 3T3-L1 preadipocytes. Mol. Endocrinol. 19, 399–408. doi:10.1210/me.2004-0324
Lequeux, C., Auxenfans, C., Mojallal, A., Sergent, M., and Damour, O. (2009). Optimization of a culture medium for the differentiation of preadipocytes into adipocytes in a monolayer. Biomed. Mater Eng. 19, 283–291. doi:10.3233/BME-2009-0593
Li, B., and Dewey, C. N. (2011). RSEM: Accurate transcript quantification from RNA-seq data with or without a reference genome. BMC Bioinforma. 12, 323. doi:10.1186/1471-2105-12-323
Li, S., He, Y., Lin, S., Hao, L., Ye, Y., Lv, L., et al. (2016). Increase of circulating cholesterol in vitamin D deficiency is linked to reduced vitamin D receptor activity via the Insig-2/SREBP-2 pathway. Mol. Nutr. Food Res. 60, 798–809. doi:10.1002/mnfr.201500425
Li, R., Zhang, Y., Zhu, W., Ding, C., Dai, W., Su, X., et al. (2020). Effects of olanzapine treatment on lipid profiles in patients with schizophrenia: A systematic review and meta-analysis. Sci. Rep. 10, 17028. doi:10.1038/s41598-020-73983-4
Lieberman, J. A., Stroup, T. S., McEvoy, J. P., Swartz, M. S., Rosenheck, R. A., Perkins, D. O., et al. (2005). Effectiveness of antipsychotic drugs in patients with chronic schizophrenia. N. Engl. J. Med. 353, 1209–1223. doi:10.1056/NEJMoa051688
Liu, X. M., Zhao, X. M., Deng, C., Zeng, Y. P., and Hu, C. H. (2019). Simvastatin improves olanzapine-induced dyslipidemia in rats through inhibiting hepatic mTOR signaling pathway. Acta Pharmacol. Sin. 40, 1049–1057. doi:10.1038/s41401-019-0212-1
Lorenzo-Villalba, N., Andres, E., and Meyer, A. (2021). Chronic onset form of anti-HMG-CoA reductase myopathy. Eur. J. Case Rep. Intern Med. 8, 002672. doi:10.12890/2021_002672
Mathews, J., Newcomer, J. W., Mathews, J. R., Fales, C. L., Pierce, K. J., Akers, B. K., et al. (2012). Neural correlates of weight gain with olanzapine. Arch. Gen. Psychiatry 69, 1226–1237. doi:10.1001/archgenpsychiatry.2012.934
Mirhosseini, N., Rainsbury, J., and Kimball, S. M. (2018). Vitamin D supplementation, serum 25(OH)D concentrations and cardiovascular disease risk factors: A systematic review and meta-analysis. Front. Cardiovasc Med. 5, 87. doi:10.3389/fcvm.2018.00087
Nagaoka, K., Nagashima, T., Asaoka, N., Yamamoto, H., Toda, C., Kayanuma, G., et al. (2021). Striatal TRPV1 activation by acetaminophen ameliorates dopamine D2 receptor antagonist-induced orofacial dyskinesia. JCI Insight 6, e145632. doi:10.1172/jci.insight.145632
Nagashima, T., Shirakawa, H., Nakagawa, T., and Kaneko, S. (2016). Prevention of antipsychotic-induced hyperglycaemia by vitamin D: A data mining prediction followed by experimental exploration of the molecular mechanism. Sci. Rep. 6, 26375. doi:10.1038/srep26375
Nagashima, T., Hayakawa, T., Akimoto, H., Minagawa, K., Takahashi, Y., and Asai, S. (2022). Identifying antidepressants less likely to cause hyponatremia: Triangulation of retrospective cohort, disproportionality, and pharmacodynamic studies. Clin. Pharmacol. Ther. 111, 1258–1267. doi:10.1002/cpt.2573
Narula, P. K., Rehan, H. S., Unni, K. E., and Gupta, N. (2010). Topiramate for prevention of olanzapine associated weight gain and metabolic dysfunction in schizophrenia: A double-blind, placebo-controlled trial. Schizophr. Res. 118, 218–223. doi:10.1016/j.schres.2010.02.001
Narvaez, C. J., Matthews, D., Broun, E., Chan, M., and Welsh, J. (2009). Lean phenotype and resistance to diet-induced obesity in vitamin D receptor knockout mice correlates with induction of uncoupling protein-1 in white adipose tissue. Endocrinology 150, 651–661. doi:10.1210/en.2008-1118
Nimitphong, H., Holick, M. F., Fried, S. K., and Lee, M. J. (2012). 25-hydroxyvitamin D₃ and 1,25-dihydroxyvitamin D₃ promote the differentiation of human subcutaneous preadipocytes. PLoS One 7, e52171. doi:10.1371/journal.pone.0052171
Oh, J., Weng, S., Felton, S. K., Bhandare, S., Riek, A., Butler, B., et al. (2009). 1,25(OH)2 vitamin d inhibits foam cell formation and suppresses macrophage cholesterol uptake in patients with type 2 diabetes mellitus. Circulation 120, 687–698. doi:10.1161/CIRCULATIONAHA.109.856070
Ohashi, K., Deyashiki, A., Miyake, T., Nagayasu, K., Shibasaki, K., Shirakawa, H., et al. (2018). TRPV4 is functionally expressed in oligodendrocyte precursor cells and increases their proliferation. Pflugers Arch. 470, 705–716. doi:10.1007/s00424-018-2130-3
Osaki, Y., Nakagawa, Y., Miyahara, S., Iwasaki, H., Ishii, A., Matsuzaka, T., et al. (2015). Skeletal muscle-specific HMG-CoA reductase knockout mice exhibit rhabdomyolysis: A model for statin-induced myopathy. Biochem. Biophys. Res. Commun. 466, 536–540. doi:10.1016/j.bbrc.2015.09.065
Reid, I. R., Ames, R., Mason, B., Bolland, M. J., Bacon, C. J., Reid, H. E., et al. (2010). Effects of calcium supplementation on lipids, blood pressure, and body composition in healthy older men: A randomized controlled trial. Am. J. Clin. Nutr. 91, 131–139. doi:10.3945/ajcn.2009.28097
Saha, S., Chant, D., Welham, J., and McGrath, J. (2005). A systematic review of the prevalence of schizophrenia. PLoS Med. 2, e141. doi:10.1371/journal.pmed.0020141
Siswanto, S., Yamamoto, H., Furuta, H., Kobayashi, M., Nagashima, T., Kayanuma, G., et al. (2021). Drug repurposing prediction and validation from clinical big data for the effective treatment of interstitial lung disease. Front. Pharmacol. 12, 635293. doi:10.3389/fphar.2021.635293
Teske, K. A., Yu, O., and Arnold, L. A. (2016). Inhibitors for the Vitamin D receptor-coregulator interaction. Vitam. Horm. 100, 45–82. doi:10.1016/bs.vh.2015.10.002
Wang, J., Liu, Y. S., Zhu, W. X., Zhang, F. Q., and Zhou, Z. H. (2014). Olanzapine-induced weight gain plays a key role in the potential cardiovascular risk: Evidence from heart rate variability analysis. Sci. Rep. 4, 7394. doi:10.1038/srep07394
Wong, K. E., Szeto, F. L., Zhang, W., Ye, H., Kong, J., Zhang, Z., et al. (2009). Involvement of the vitamin D receptor in energy metabolism: Regulation of uncoupling proteins. Am. J. Physiol. Endocrinol. Metab. 296, E820–E828. doi:10.1152/ajpendo.90763.2008
Wong, K. E., Kong, J., Zhang, W., Szeto, F. L., Ye, H., Deb, D. K., et al. (2011). Targeted expression of human vitamin D receptor in adipocytes decreases energy expenditure and induces obesity in mice. J. Biol. Chem. 286, 33804–33810. doi:10.1074/jbc.M111.257568
Yabe, D., Brown, M. S., and Goldstein, J. L. (2002). Insig-2, a second endoplasmic reticulum protein that binds SCAP and blocks export of sterol regulatory element-binding proteins. Proc. Natl. Acad. Sci. U. S. A. 99, 12753–12758. doi:10.1073/pnas.162488899
Yan, H., Chen, J. D., and Zheng, X. Y. (2013). Potential mechanisms of atypical antipsychotic-induced hypertriglyceridemia. Psychopharmacol. Berl. 229, 1–7. doi:10.1007/s00213-013-3193-7
Keywords: olanzapine, vitamin D, C2C12 cells, dyslipidemia, clinical big data, low-density lipoprotein, high-density lipoprotein
Citation: Zhou Z, Nagashima T, Toda C, Kobayashi M, Suzuki T, Nagayasu K, Shirakawa H, Asai S and Kaneko S (2023) Vitamin D supplementation is effective for olanzapine-induced dyslipidemia. Front. Pharmacol. 14:1135516. doi: 10.3389/fphar.2023.1135516
Received: 01 January 2023; Accepted: 13 February 2023;
Published: 21 February 2023.
Edited by:
Ali H. Eid, Qatar University, QatarReviewed by:
Wei Zhong, University of North Carolina at Greensboro, United StatesLana Nezic, University of Banja Luka, Bosnia and Herzegovina
Janne Hukkanen, University of Oulu, Finland
Chiara Bianca Maria Platania, University of Catania, Italy
Copyright © 2023 Zhou, Nagashima, Toda, Kobayashi, Suzuki, Nagayasu, Shirakawa, Asai and Kaneko. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Shuji Kaneko, c2thbmVrb0BwaGFybS5reW90by11LmFjLmpw