Wei Zhou1†
Yong Xu2†
Qinyu Lv3
Yong-hui Sheng4Luan Chen1Mo Li1Lu Shen1
Cong Huai1
Zhenghui Yi3*
Donghong Cui3,5*
Shengying Qin1,6*- 1Key Laboratory for the Genetics of Developmental and Neuropsychiatric Disorders, Bio-X Institutes, Ministry of Education, Shanghai Jiao Tong University, Shanghai, China
- 2Department of Psychiatry, First Hospital of Shanxi Medical University, Taiyuan, China
- 3Shanghai Mental Health Center, Shanghai Jiao Tong University School of Medicine, Shanghai, China
- 4Shandong Mental Health Center, Jinan, China
- 5Shanghai Key Laboratory of Psychotic Disorders, Shanghai, China
- 6The Third Affiliated Hospital, Guangzhou Medical University, Guangzhou, China
Olanzapine, a second-generation antipsychotic medication, plays a critical role in current treatment of schizophrenia (SCZ). It has been observed that the olanzapine responses in schizophrenia treatment are different across individuals. However, prediction of this individual-specific olanzapine response requires in-depth knowledge of biomarkers of drug response. Here, we performed an integrative investigation on 238 Han Chinese SCZ patients to identify predictive biomarkers that were associated with the efficacy of olanzapine treatment. This study applied HaloPlex technology to sequence 143 genes from 79 Han Chinese SCZ patients. Our result suggested that there were 12 single nucleotide polymorphisms (SNPs) had significant association with olanzapine response in Han Chinese SCZ patients. Using MassARRAY platform, we tested that if these 12 SNPs were also statistically significant in 159 other SCZ patients (independent cohort) and the combined 238 SCZ patients (composed of two tested cohorts). The result of this analysis showed that 2 SNPs were significantly associated with the olanzapine response in both independent cohorts (rs324026, P = 0.023; rs12610827, P = 0.043) and combined SCZ patient population (rs324026, adjust P = 0.014; rs12610827, adjust P = 0.012). Our study provides systematic analyses of genetic variants associated with olanzapine responses of Han Chinese SCZ patients. The discovery of these novel biomarkers of olanzapine-response will facilitate to advance future olanzapine treatment specific for Han Chinese SCZ patients.
Introduction
Schizophrenia (SCZ) is a severe chronic neuropsychiatric illness. According to a survey carried out in 33 countries, about 15 out of 100,000 individuals were suffering from SCZ globally (McGrath et al., 2004). These SCZ patients were estimated to have higher risk of death (about 2.5 times more) compared to healthy individuals (McGrath et al., 2004; Saha et al., 2007). Patients with SCZ disorder require long-term treatments to prevent themselves from illness progression or symptom relapse (Howes et al., 2015; Chong et al., 2016). Although the development of SCZ has long been regarded as caused by a combination of genetic and environmental factors, detailed pathophysiological mechanism of SCZ still remains unclear.
Second-generation antipsychotic (SGA) medications are widely considered as the most advanced and effective treatment for SCZ patients nowadays. These SGA includes olanzapine (OLA), risperidone and quetiapine (Owen et al., 2016). However, SCZ patients who received SGA treatments often experienced severe adverse drug reactions (ADRs) (Zhang and Malhotra, 2011). In fact, current SGA therapies could be seen as a subjective “trial and error” process. For some of the SCZ patients, these SGA treatments did not exert any therapeutic effects on their symptoms. This individual-specific response to SGA may be caused by the association between individual-specific genetic variation and the efficacy of SGAs (Zhang and Malhotra, 2011). In order to address such specificity of SGA response across different patients, scientists in pharmacogenomics field is now exploring the possibility of predicting drug response using individual-specific genetic signatures (Arranz and de Leon, 2007).
Olanzapine is one of the most commonly used SGAs. It has been shown to have relatively superior efficacy in various clinical trials compared to other SGAs (Lieberman et al., 2005; Leucht et al., 2009). It has been reported that OLA treatment had relatively low extrapyramidal side effects and better efficacy to minimize negative symptoms of SCZ patients when it is used in clinically effective doses (Meltzer, 1999). OLA binds to serotonin type 2 (5-HT2) and dopamine (D2) receptors with high affinity in patients’ body. Diphosphate glucuronosyltransferases (UGT), a member of cytochrome P450 family and flavin-containing monooxygenase 1 (FMO1), catalyze the oxidative hepatic metabolism process of OLA (Ring et al., 1996; Kassahun et al., 1997; Linnet, 2002). On the other hand, due to the heterogeneity in different SCZ patients, not all patients respond to OLA treatment adequately well as we expected. Some patients who received OLA therapy even experienced severe adverse side effects that resulted in non-compliance with drug treatment (Zhang and Malhotra, 2011; Musil et al., 2015). If these SCZ patients with no other effective therapies specific for them they would have to face the coming disease progression, relapses and potential long-term hospitalizations (Robinson et al., 1999; King et al., 2014).
Although numerous studies have been performed on the factors that influence the therapeutic efficacy of OLA, there were very few of them focused on the individual-specific genetic biomarkers of OLA response (Söderberg and Dahl, 2013). In comparison, earlier attempts to search for biomarkers of OLA response focused mainly on the relationship between OLA response and its metabolic pathways, including glucuronidation, hydroxylation, N-demethylation and N-oxidation pathways (Laika et al., 2010; Haslemo et al., 2012; Mao et al., 2012; Söderberg et al., 2013; Brandl et al., 2015). A number of genetic variants, including UGT2B10 rs61750900 (UGT2B10∗2) (Erickson-Ridout et al., 2011), CYP1A2 rs762551 (CYP1A2∗1F) (Laika et al., 2010), DRD3 rs6280 (Adams et al., 2008), AHR rs4410790 (Söderberg et al., 2013), FMO3 K158–G308, FMO1 rs12720462 (FMO1∗6) and FMO1 rs7877 have been reported to play important roles to influence OLA metabolism (Soderberg et al., 2013). In addition, the drug response and the pharmacokinetics of OLA have also been found to associate with genetic elements that are not directly involved in the metabolic pathway of OLA (Lin et al., 2006; Meary et al., 2008; Cabaleiro et al., 2013; Yu et al., 2018). For example, P-glycoprotein, a membrane protein that pumps foreign substances out of cells and is regarded as element that is not directly related to OLA metabolic pathway, affects the penetration of OLA into the central nervous system (Lin et al., 2006). These discoveries suggested that a comprehensive study of OLA response requires clear understanding of the complicated biological network that is composed of enzymes involved in drug metabolism, drug transportation and drugs targeted receptors.
In this study, we investigated the associations between SNPs in 143 genes and the OLA response of 79 Han Chinese SCZ patients using target-sequencing technology. The newly found biomarkers was considered as genetic signature of drug responses to 8-week treatment with OLA and were validated in the other independent Han Chinese SCZ patient cohort.
Materials and Methods
Subjects
In this study, we collected 2 independent sets of OLA response data from Han Chinese SCZ patients in order to validate our discoveries. We named the first set as ‘discovery cohort’ and other one as ‘independent cohort.’ The demographics and clinical details of the both sets of patients are demonstrated in Table 1.
Table 1. Demographic and clinical details of patients suffering from schizophrenia in three cohorts.
The discovery cohort was composed of 79 recruited Han Chinese SCZ patients who had been treated with OLA from the Shanghai Mental Health Center of China. It comprised 37 males and 42 females. The mean age of them was 43.1 ± 18.3 years old (Table 1). SCZ of the patients was diagnosed according to the criteria of the Diagnostic and Statistical Manual of Mental Disorders, Fourth Edition (DSM-IV) and confirmed by at least two experienced psychiatrists. Patients who had physical complications or other substance abuse were excluded from this investigation. Our analysis only considered patients who had not been previously treated with atypical antipsychotics and had not received any medication for more than 4 weeks before their enrollment in this study.
The independent cohort contained 159 recruited Han Chinese SCZ patients who were undergoing OLA monotherapy from the Shanghai Mental Health Center of China and the First Hospital of Shanxi Medical University. They were composed of 75 males and 84 females. The mean age of independent cohort was 38.5 ± 16.4 years (Table 1). This cohort of patients was used to validate the novel biomarkers found by the analysis of discovery cohort. The selection criteria were similar to the discovery cohort. Finally, patient data from discovery cohort and independent cohort were combined in order to gain a greater power in statistical analysis. Therefore, a total of 238 patients were used to perform another validation of the newly found biomarkers.
Clinical Assessment
Clinical effects of patients were evaluated using the Positive and Negative Syndrome Scale (PANSS) by two fully qualified psychiatrists during the 8 weeks of OLA treatment. The inter-rater reliability between the two raters was found to be high since the intraclass correlation coefficients (ICCs) was larger than 0.8. Based to Obermeier’s method, patients were classified as good responders (reduction of PNASS score ≥ 50%) and poor responders (reduction of PNASS score < 50%) for analysis (Obermeier et al., 2010). The initial daily dose of OLA was 10 mg per day and then it gradually increased to 15 mg per day within the 1st week. After that, the dosage was adjusted based on individual tolerance to the treatment. During the medication period, nursing staff closely monitored any medication compliance occurred in patients. No other drugs were administered when OLA monotherapy was performed, except for sennoside for constipation, flunitrazepam or lorazepam for acute insomnia and biperiden for any extrapyramidal side effects.
Ethics Statement
The study was approved by the Ethical Committee of Human Genetic Resources in Shanghai, China. All subjects or their legal guardians understood the procedure and had given written informed consent to their participation in this study according to the Declaration of Helsinki (Human, 1999).
Targeted Genes and Capture Design
One hundred and forty-three genes were selected for targeted sequencing were based on their involvement in drug metabolism (including genes encode for drug-metabolizing enzymes, drug-transporting enzymes and the receptors mediating drug responses) from PharmGKB database1 and the relevant literatures (Arranz et al., 2011; Li and Bluth, 2011; Arranz et al., 2016) that reported potential genes related to SGAs efficacy. We aimed to study these genes in order to investigate novel biomarkers of OLA response. The details of these 143 genes are listed in Supplementary Table S1.
Sequencing probes for the 143 targeted genes were designed using Agilent’s SureDesign tool2. Targeted regions of these genes of interest included their coding regions ± 10 bp and untranslated regions (UTR) according to information from RefSeq, Ensembl, CCDS, and GENCODE databases (Harrow et al., 2012; Pruitt et al., 2014; Cunningham et al., 2015).
Library Preparation and Next Generation Sequencing
Genomic DNA was extracted from whole blood using a QIAamp DNA Blood Mini Kit (Qiagen GmbH, Hilden, Germany). The quantity and quality of the genomic DNA were measured by Nanodrop 2000 (Thermo Scientific, United States). Then we adjusted the genomic DNA to a final concentration of 100 ng/μl with high-purity water and stored at -20°C. Libraries were prepared with a HaloPlex Target Enrichment System Kit (Agilent Technologies, Santa Clara, CA, United States) following the manufacturer’s instructions. Libraries were then quantified using the Agilent 2100 Bioanalyzer (Agilent Technologies). Sequencing was performed with the HiSeq 2500 platform (Illumina, San Diego, CA, United States) using paired-end libraries (2 × 101-bp).
Raw data were processed following standard protocols used in earlier reports (Gaynor et al., 2016). In short, raw image files were first converted to the FASTQ format and the reads were aligned to the human reference genome (hg19, GRCh37). SNPs were identified according to GATK standard hard filtering parameters (DePristo et al., 2011). On average, 99% of reads covered >80× and 81% >200×, which suggested that the coverage was sufficiently high to detect variants with appropriate sensitivity. The program ANNOVAR was used to annotate SNVs that covered >20× according to the information from Ensembl Variation, dbSNP, and 1000genome database (Sherry et al., 2001; Abecasis et al., 2010; Wang et al., 2010; Flicek et al., 2012). Subsequently, individual and SNP-level quality controls were performed using PLINK (v1.07) software (Purcell et al., 2007). Data cleaning was performed according to the following criteria: genotypic call rate < 95%, Hardy–Weinberg equilibrium (HWE) < 0.001, and minor allele frequency (MAF) < 0.01. After we conducted quality controls of the sequences, there were 77 individuals (38 good responders and 39 poor responders) and 807 SNPs remained in our data for later analysis.
Validation Trial
We identified 12 SNPs that were significantly associated with OLA response using discovery cohort. Similar analysis was performed using both independent cohort and patient samples with two cohort combined in order to validate these newly found SNPs. These 12 SNPs were genotyped using the Sequenom MassARRAY platform (Agena Bioscience, San Diego, CA, United States) following manufacturer’s instructions. MassARRAY primers were designed using a semi-automated software Assay Design Suite v2.03. The primer sequences are listed in Supplementary Table S2. Data cleaning was performed according to the following criteria: genotypic call rate < 95%, Hardy–Weinberg equilibrium (HWE) < 0.001, and minor allele frequency (MAF) < 0.01. Similarly, to the earlier processing of discovery cohort, quality control was carried on and only the filtered data was used for our analysis.
Statistical Analyses
The demographic characteristics of the both ‘good responder’ and ‘poor responder’ groups were examined to confirm the homogeneity of the data used in our analysis. The data was found to have normal distribution, allowing student’s t-tests to be performed on the obtained data (age, PANSS score, etc.). Gender differences were analyzed using the Chi-square test. SPSS software (version 11.0, Chicago, IL, United States) was used for all the statistical analyses in this study. The association between genotype and OLA response was assessed using logistic regression model by PLINK vl.07 software (Purcell et al., 2007). P-values were corrected using Bonferroni method for multiple testing adjustments. Two-tailed P-values of 0.05 were considered to be statistically significant. Power analysis was performed by the software GPower 3.1.
Results
Patients Characteristics and Sequencing Profile
The demographic and clinical details of patient subjects included in this study are shown in Table 1. Among the 79 patients in the discovery cohort, 38 patients were defined to be good responders to OLA while 41 of them were poor responders. On the other hand, among 159 patients in independent cohort, 119 patients were good responders and 40 as poor responders to OLA. In the total cohort that comprised both sets, there were 157 good responders and 81 poor responders out of a total of 238 patient subjects. There was no statistically significant difference in the baseline characteristics between good responders and poor responders, except in the case of the PANSS total scores at the 8-week endpoint, meaning that the population was homogeneous (Table 1).
Effects of Individual Polymorphisms on the OLA Response in the Discovery Cohort
Twelve out of 807 tested SNPs were found to be significantly associated with OLA response of 79 Han Chinese SCZ patients in discovery cohort. Table 2 lists the results of the SNP association analysis of pharmacogenetic impact on OLA treatment response (P < 0.05). Two newly found variants were located on the exon’s region (rs6280, P = 0.026, OR = 3.0, 95% CI = 1.14–7.87; rs2011404, P = 0.04, OR = 5.4, 95% CI = 1.08–26.93). The other 10 variants, which were not located in exons regions, were also found to be significantly associated with OLA response. However, there were no variants remained statistically significant after multiple-testing corrections (data not shown).
Table 2. 12 SNPs that were found to be significantly associated with responses to olanzapine treatment, by targeted sequencing.
Verification of the Genetic Variants Associated With the Response to OLA in the Independent and Total Cohort
We used independent cohort and total sample population composed of both discovery cohort and independent cohort to validate the 12 SNP signatures found from the analysis of discovery cohort. The relevant clinical information of the data is shown in Table 1. A total of 12 SNPs was genotyped from patients in these two sets. In particular, SNP rs324026 and rs12610827 were found to be significantly associated with OLA treatment response in the independent cohort (P = 0.023 and P = 0.023). In the combined cohort, 4 SNPs displayed significant difference in OLA response between good responders and poor responders. We obtained strong evidence to conclude that these 2 variants (rs324026 and rs6280) in the dopamine receptor D3 (DRD3) gene were significantly associated with OLA response in Han Chinese SCZ patients (P = 0.001 and 0.0047). In addition, SNP rs12610827 (near to PLK5) and rs1543494 (located in SUPT16H) were also shown to be significantly associated with the OLA treatment response (P = 0.001 and 0.038). Detailed information of these significant SNPs is shown in Table 3.
Table 3. Validation of SNPs associated with the olanzapine response.
Power Analysis
Post hoc power analysis revealed that the statistical power of the discovery cohort size (n = 79) in detecting a significant association (P < 0.05) was 0.76 with a medium effect size (Odds ratio = 2.0). The power of independent cohort size (n = 159) was 0.96 with the same effect size. These results indicated that the sample size in our study was sufficient to achieve a considerably low risk of a type II error.
Discussion
To date, most pharmacogenomic studies on the OLA response focused on a few genes that are known to be relevant to OLA metabolism. Our study represents a more systematic survey of genetic biomarkers, including drug metabolic enzyme genes, receptor genes and other related genes. 143 genes of interest were sequenced using Next-generation sequencing technology for our association analysis. This study is one of the most comprehensive pharmacogenetic analyses of association between SNP variants and OLA response.
Our result suggested that SNP rs324026 in DRD3 gene had significant association with OLA response using independent cohort. This difference still remained significant in the total cohort comprised 2 cohorts even after Bonferroni correction. However, the other variant rs6280 was only found to be significantly associated with an 8-week treatment of OLA response in the combined cohort population and did not have evidence to have significant associations with OLA response in independent cohort. This inconsistent result may be caused by the small sample set we tested. Notably, these 2 SNPs both exhibited strong linkage disequilibrium (r2 > 0.9) in the HaploReg database (Ward and Kellis, 2012). Therefore, both rs6280 and rs324026 may serve as biomarkers of OLA treatment response.
It is known that rs6280 mutation leads to a glycine for serine substitution and is associated with altered dopamine binding affinity. This glycine variant had been suggested to be able to increase the densities of the dopamine receptor D3 (DRD-3) in some areas in human brain (Jeanneteau et al., 2006). Adams et al. (2008) reported that DRD-3 gly/gly genotype and other polymorphisms in linkage disequilibrium with ser-9-gly variant were significantly associated with an increase in PANSS total score. Therefore, we concluded that our result was consistent with the earlier discoveries of the association between ser-9-gly variant and clozapine, which is the most similar receptor binding profile to OLA. Cerrato et al. (2017) surveyed 65 papers and found that rs6280 was successfully replicated as prognostic biomarkers of clozapine efficacy. In contrast, rs324026 variant had never been reported to affect the therapeutic efficacy of OLA. Rs324026 is located next to exon 5. Our analysis results suggested that we could only find evidence of significant association between SNP rs324026 and OLA efficacy after Bonferroni correction in the combined sample with both discovery and independent cohorts. Additionally, individuals with C alleles of rs324026 generally experience significantly better efficacy of OLA treatment.
In this study, SNP rs12610827 was validated in patients from independent cohort and its association with OLA response remained significant for multiple testing in the total cohort after Bonferroni correction. Rs12610827 variant is located near the PLK5 gene. Polo-like kinases (Plks) family, consisted of 5 members (Plk1-Plk5), is traditionally regarded to play an important part only in cell cycle progression. However, mounting evidence showed that Plk2 and Plk5 are also closely involved in neuron biology (de Carcer et al., 2011a). It had been suggested that Plk2 modulates neurite formation in response to activities of brain-derived growth factor (BDGF) (Inglis et al., 2009). Additionally, Plk5 was highly expressed in the central nervous system and it serve as a Plk2-like role in the cerebellum according earlier report (de Carcer et al., 2011b). It had been suggested that Plk5 was regulated by CpG methylation of the promoter region on the transcriptional level (de Carcer et al., 2011a). The level of PLK5 gene expression may be influenced by the methylation status of this variant. In this study, Han Chinese SCZ patients who carried allele T in SNP rs12610827 showed more good response in OLA treatment. However, no previous reports have found such association between PLK5 variants and drug response. Therefore, we believe rs12610827 variant is worthy of further investigations in order to verify its influence on OLA response in the future.
This study has several limitations. First, our analysis did not consider some other genes that may have significant association with antipsychotic. Therefore, our results may have neglected some important biomarkers of OLA response due to this incomplete gene collection. Secondly, a number of identified SNP associations failed to stay statistically significant after Bonferroni corrections. This may be caused by over-correction because sample size was relatively small. Therefore, employing strict multiple corrections such as Bonferroni to the data may be too harsh for this specific study.
Conclusion
In sum, we performed a comprehensive study on 238 Han Chinese SCZ patients in order to identify potential biomarkers of Han Chinese-specific OLA responses. The result showed that 143 genes were significantly associated with OLA. In addition, 2 variants (rs324026 and rs12610827) were found to have significant association with the OLA response. Future investigations with larger sample sizes and high-throughput methods such as high-density SNP arrays and whole exome sequencing are warranted to find more biomarkers to predict the efficacy of OLA in the Han Chinese population.
Data Availability
The datasets used and/or analyzed during the current study are available from the corresponding author on reasonable request.
Author Contributions
WZ and YX performed the experiments and wrote the manuscript. CH and ML aided in processing the data. QL, Y-hS, and YX aided in the collection of the materials. LC and LS helped in revising the manuscript. ZY, DC, and SQ designed and revised the manuscript.
Funding
This work was supported by grants from the 863 Program (Grant Numbers 2012AA02A515 and 2012AA021802), the National Natural Science Foundation of China (Grant Numbers 81421061, 81273596, J1210047, 30900799, 81361120389, 30972823, 81671326, and 81671336), National Key Research and Development Program (Grant Numbers 2016YFC0905000, 2016YFC0905002, 2016YFC1200200, 2016YFC0906400, and 2017YFC0909200), The 4th Three-Year Action Plan for Public Health of Shanghai (the Project No. 15GWZK0101), Shanghai Key Laboratory of Psychotic Disorders (Grant Number 13dz2260500), Public Science and Technology Research Funds (Grant Number 201210056), The Fourth Round of Shanghai Three-Year Action Plan on Public Health Discipline and Talent Program: Women and Children’s Health (Grant Number 15GWZK0401), the Shanghai Jiao Tong University Interdisciplinary Research fund, and the Shanghai Leading Academic Discipline Project (Grant Number B205).
Conflict of Interest Statement
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Supplementary Material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fphar.2019.00177/full#supplementary-material
Footnotes
References
Abecasis, G. R., Altshuler, D., Auton, A., Brooks, L. D., Durbin, R. M., and Gibbs, R. A. (2010). A map of human genome variation from population-scale sequencing. Nature 7319, 1061–1073. doi: 10.1038/nature09534
Adams, D. H., Close, S., Farmen, M., Downing, A. M., Breier, A., and Houston, J. P. (2008). Dopamine receptor D3 genotype association with greater acute positive symptom remission with olanzapine therapy in predominately caucasian patients with chronic schizophrenia or schizoaffective disorder. Hum. Psychopharmacol. 4, 267–274. doi: 10.1002/hup.930
Arranz, M. J., and de Leon, J. (2007). Pharmacogenetics and pharmacogenomics of schizophrenia: a review of last decade of research. Mol. Psychiatry 8, 707–747. doi: 10.1038/sj.mp.4002009
Arranz, M. J., Gallego, C., Salazar, J., and Arias, B. (2016). Pharmacogenetic studies of drug response in schizophrenia. Expert Rev. Precis. Med. Drug Dev. 1, 79–91. doi: 10.1080/23808993.2016.1140554
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
Brandl, E. J., Chowdhury, N. I., Tiwari, A. K., Lett, T. A. P., Meltzer, H. Y., and Kennedy, J. L. (2015). Genetic variation in CYP3A43 is associated with response to antipsychotic medication. J. Neural Transm. 1, 29–34. doi: 10.1007/s00702-014-1298-8
Cabaleiro, T., Lopez-Rodriguez, R., Ochoa, D., Roman, M., Novalbos, J., and Abad-Santos, F. (2013). Polymorphisms influencing olanzapine metabolism and adverse effects in healthy subjects. Hum. Psychopharmacol. 3, 205–214. doi: 10.1002/hup.2308
Cerrato, F., Guizzaro, L., Scudellari, P., Atti, A. R., and De Ronchi, D. (2017). Predictive biomarkers in clozapine-treated patients: assessment of the evidences and suggestion for research methodology. Eur. Psychiatry 41:S378. doi: 10.1016/j.eurpsy.2017.02.403
Chong, H. Y., Teoh, S. L., Wu, D. B., Kotirum, S., Chiou, C. F., and Chaiyakunapruk, N. (2016). Global economic burden of schizophrenia: a systematic review. Neuropsychiatr. Dis. Treat. 12, 357–373. doi: 10.2147/NDT.S96649
Cunningham, F., Amode, M. R., Barrell, D., Beal, K., Billis, K., and Brent, S. (2015). Ensembl 2015. Nucleic Acids Res. 43, D662–D669. doi: 10.1093/nar/gku1010
de Carcer, G., Escobar, B., Higuero, A. M., Garcia, L., Anson, A., and Perez, G. (2011a). Plk5, a polo box domain-only protein with specific roles in neuron differentiation and glioblastoma suppression. Mol. Cell. Biol. 6, 1225–1239. doi: 10.1128/MCB.00607-10
de Carcer, G., Manning, G., and Malumbres, M. (2011b). From Plk1 to Plk5: functional evolution of polo-like kinases. Cell Cycle 14, 2255–2262.
DePristo, M. A., Banks, E., Poplin, R., Garimella, K. V., Maguire, J. R., and Hartl, C. (2011). A framework for variation discovery and genotyping using next-generation DNA sequencing data. Nat. Genet. 5, 491–498. doi: 10.1038/ng.806
Erickson-Ridout, K. K., Zhu, J., and Lazarus, P. (2011). Olanzapine metabolism and the significance of UGT1A448V and UGT2B1067Y variants. Pharmacogenet. Genomics 9, 539–551. doi: 10.1097/FPC.0b013e328348c76b
Flicek, P., Amode, M. R., Barrell, D., Beal, K., Brent, S., and Carvalho-Silva, D. (2012). Ensembl 2012. Nucleic Acids Res. 40, D84–D90. doi: 10.1093/nar/gkr991
Gaynor, S. C., Breen, M. E., Monson, E. T., de Klerk, K., Parsons, M., and DeLuca, A. P. (2016). A targeted sequencing study of glutamatergic candidate genes in suicide attempters with bipolar disorder. Am. J. Med. Genet. B Neuropsychiatr. Genet. 8, 1080–1087. doi: 10.1002/ajmg.b.32479
Harrow, J., Frankish, A., Gonzalez, J. M., Tapanari, E., Diekhans, M., and Kokocinski, F. (2012). GENCODE: the reference human genome annotation for The ENCODE Project. Genome Res. 9, 1760–1774. doi: 10.1101/gr.135350.111
Haslemo, T., Loryan, I., Ueda, N., Mannheimer, B., Bertilsson, L., and Ingelman-Sundberg, M. (2012). UGT1A4∗3 encodes significantly increased glucuronidation of olanzapine in patients on maintenance treatment and in recombinant systems. Clin. Pharmacol. Ther. 2, 221–227. doi: 10.1038/clpt.2012.46
Howes, O., McCutcheon, R., and Stone, J. (2015). Glutamate and dopamine in schizophrenia: an update for the 21st century. J. Psychopharmacol. 2, 97–115. doi: 10.1177/0269881114563634
Inglis, K. J., Chereau, D., Brigham, E. F., Chiou, S. S., Schobel, S., and Frigon, N. L. (2009). Polo-like kinase 2 (PLK2) phosphorylates alpha-synuclein at serine 129 in central nervous system. J. Biol. Chem. 5, 2598–2602. doi: 10.1074/jbc.C800206200
Jeanneteau, F., Funalot, B., Jankovic, J., Deng, H., Lagarde, J. P., and Lucotte, G. (2006). A functional variant of the dopamine D3 receptor is associated with risk and age-at-onset of essential tremor. Proc. Natl. Acad. Sci. U.S.A. 28, 10753–10758. doi: 10.1073/pnas.0508189103
Kassahun, K., Mattiuz, E., and Nyhart, E. (1997). Disposition and biotransformation of the antipsychotic agent olanzapine in humans. Drug Metab. Dispos. 25, 81–93.
King, D., Knapp, M., Patel, A., Amaddeo, F., Tansella, M., and Schene, A. (2014). The impact of non-adherence to medication in patients with schizophrenia on health, social care and societal costs. Analysis of the QUATRO study. Epidemiol. Psychiatr. Sci. 1, 61–70. doi: 10.1017/S2045796013000097
Laika, B., Leucht, S., and Heres, S. (2010). Pharmacogenetics and olanzapine treatment: CYP1A2∗1F and serotonergic polymorphisms influence therapeutic outcome. Pharmacogenomics J. 10, 20–29. doi: 10.1038/tpj.2009.32
Leucht, S., Corves, C., Arbter, D., Engel, R., Li, C., and Davis, J. (2009). Second-generation versus first-generation antipsychotic drugs for schizophrenia: a meta-analysis. Lancet 373, 31–41. doi: 10.1016/S0140-6736(08)61764-X
Li, J., and Bluth, M. H. (2011). Pharmacogenomics of drug metabolizing enzymes and transporters: implications for cancer therapy. Pharmgenomics Pers. Med. 4, 11–33. doi: 10.2147/PGPM.S18861
Lieberman, J., Stroup, T., and McEvoy, J. (2005). Effectiveness of antipsychotic drugs in patients with chronic schizophrenia. N. Engl. J. Med. 353, 1209–1223. doi: 10.1056/NEJMoa051688
Lin, Y., Ellingrod, V. L., Bishop, J. R., and Miller, D. D. (2006). The relationship between P-glycoprotein (PGP) polymorphisms and response to olanzapine treatment in schizophrenia. Ther. Drug Monit. 5, 668–672. doi: 10.1097/01.ftd.0000246761.82377.a6
Linnet, K. (2002). Glucuronidation of olanzapine by cDNA-expressed human UDP-glucuronosyltransferases and human liver microsomes. Hum. Psychopharmacol. 5, 233–238. doi: 10.1002/hup.403
Mao, M., Skogh, E., and Scordo, M. (2012). Interindividual variation in olanzapine concentration influenced by UGT1A4 L48V polymorphism in serum and upstream FMO polymorphisms in cerebrospinal fluid. J. Clin. Psychopharmacol. 32, 287–289. doi: 10.1097/JCP.0b013e31824997a8
McGrath, J., Saha, S., Welham, J., El, S. O., MacCauley, C., and Chant, D. (2004). A systematic review of the incidence of schizophrenia: the distribution of rates and the influence of sex, urbanicity, migrant status and methodology. BMC Med. 2:13. doi: 10.1186/1741-7015-2-13
Meary, A., Brousse, G., Jamain, S., Schmitt, A., Szoke, A., and Schurhoff, F. (2008). Pharmacogenetic study of atypical antipsychotic drug response: involvement of the norepinephrine transporter gene. Am. J. Med. Genet. B Neuropsychiatr. Genet. 4, 491–494. doi: 10.1002/ajmg.b.30635
Meltzer, H. (1999). The role of serotonin in antipsychotic drug action. Neuropsychopharmacology 21, 106S–115S. doi: 10.1016/S0893-133X(99)00046-9
Musil, R., Obermeier, M., Russ, P., and Hamerle, M. (2015). Weight gain and antipsychotics: a drug safety review. Expert Opin. Drug Saf. 1, 73–96. doi: 10.1517/14740338.2015.974549
Obermeier, M., Mayr, A., and Schennach-Wolff, R. (2010). Should the PANSS be rescaled? Schizophr. Bull. 36, 455–460. doi: 10.1093/schbul/sbp124
Owen, M. J., Sawa, A., and Mortensen, P. B. (2016). Schizophrenia. Lancet 388, 86–97. doi: 10.1016/S0140-6736(15)01121-6
Pruitt, K. D., Brown, G. R., Hiatt, S. M., Thibaud-Nissen, F., Astashyn, A., and Ermolaeva, O. (2014). RefSeq: an update on mammalian reference sequences. Nucleic Acids Res. 42, D756–D763. doi: 10.1093/nar/gkt1114
Purcell, S., Neale, B., Todd-Brown, K., Thomas, L., Ferreira, M. A. R., and Bender, D. (2007). PLINK: a tool set for whole-genome association and population-based linkage analyses. Am. J. Hum. Genet. 3, 559–575. doi: 10.1086/519795
Ring, B. J., Catlow, J., Lindsay, T. J., Gillespie, T., Roskos, L. K., and Cerimele, B. J. (1996). Identification of the human cytochromes P450 responsible for the in vitro formation of the major oxidative metabolites of the antipsychotic agent olanzapine. J. Pharmacol. Exp. Ther. 2, 658–666.
Robinson, D., Woerner, M., Alvir, J., Bilder, R., Goldman, R., and Geisler, S. (1999). Predictors of relapse following response from a first episode of schizophrenia or schizoaffective disorder. Arch. Gen. Psychiatry 56, 241–247. doi: 10.1001/archpsyc.56.3.241
Saha, S., Chant, D., and McGrath, J. (2007). A systematic review of mortality in schizophrenia - Is the differential mortality gap worsening over time? Arch. Gen. Psychiatry 10, 1123–1131. doi: 10.1001/archpsyc.64.10.1123
Sherry, S. T., Ward, M. H., Kholodov, M., Baker, J., Phan, L., and Smigielski, E. M. (2001). dbSNP: the NCBI database of genetic variation. Nucleic Acids Res. 1, 308–311. doi: 10.1093/nar/29.1.308
Söderberg, M. M., and Dahl, M. (2013). Pharmacogenetics of olanzapine metabolism. Pharmacogenomics 11, 1319–1336. doi: 10.2217/pgs.13.120
Söderberg, M. M., Haslemo, T., Molden, E., and Dahl, M. (2013). Influence of CYP1A1/CYP1A2 and AHR polymorphisms on systemic olanzapine exposure. Pharmacogenet. Genomics 5, 279–285. doi: 10.1097/FPC.0b013e3283602876
Soderberg, M. M., Haslemo, T., Molden, E., and Dahl, M. L. (2013). Influence of FMO1 and 3 polymorphisms on serum olanzapine and its N-oxide metabolite in psychiatric patients. Pharmacogenomics J. 6, 544–550. doi: 10.1038/tpj.2012.47
Wang, K., Li, M., and Hakonarson, H. (2010). ANNOVAR: functional annotation of genetic variants from high-throughput sequencing data. Nucleic Acids Res. 16:e164. doi: 10.1093/nar/gkq603
Ward, L. D., and Kellis, M. (2012). HaploReg: a resource for exploring chromatin states, conservation, and regulatory motif alterations within sets of genetically linked variants. Nucleic Acids Res. 40, D930–D934. doi: 10.1093/nar/gkr917
Yu, H., Yan, H., Wang, L., Li, J., Tan, L., and Deng, W. (2018). Five novel loci associated with antipsychotic treatment response in patients with schizophrenia: a genome-wide association study. Lancet Psychiatry 5, 327–338. doi: 10.1016/S2215-0366(18)30049-X
Keywords: olanzapine, polymorphism, schizophrenia, pharmacogenetics, biomarker, association study
Citation: Zhou W, Xu Y, Lv Q, Sheng Y-h, Chen L, Li M, Shen L, Huai C, Yi Z, Cui D and Qin S (2019) Genetic Association of Olanzapine Treatment Response in Han Chinese Schizophrenia Patients. Front. Pharmacol. 10:177. doi: 10.3389/fphar.2019.00177
Received: 03 September 2018; Accepted: 11 February 2019;
Published: 04 March 2019.
Edited by:
Amit V. Pandey, University of Bern, SwitzerlandReviewed by:
Ming Ta Michael Lee, Geisinger Health System, United StatesVindhya Udhane, Medical College of Wisconsin, United States
Julio Benitez, Universidad de Extremadura, Spain
Copyright © 2019 Zhou, Xu, Lv, Sheng, Chen, Li, Shen, Huai, Yi, Cui and Qin. 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: Zhenghui Yi, eWl6aGVuZ2h1aTE5NzFAMTYzLmNvbQ== Donghong Cui, bWFueXVjY0AxMjYuY29t Shengying Qin, Y2hpbnNpckBzanR1LmVkdS5jbg==
†These authors have contributed equally to this work