Skip to main content

Association of DRD2, DRD4 and COMT genes variants and their gene-gene interactions with antipsychotic treatment response in patients with schizophrenia


Antipsychotic drugs are the first line of treatment in schizophrenia; although antipsychotic responses indicate a wide interindividual variety in patients with schizophrenia. This study aimed to investigate the association between four polymorphisms in DRD2, DRD4 and COMT genes and their gene-gene interactions with antipsychotic treatment response in patients with schizophrenia. A total of 101 patients with schizophrenia were recruited and stratified in treatment responder and treatment resistant groups based on the published criteria of resistant to treatment using PANSS. Clinical and demographic factors were analyzed. Genomic DNA was extracted from whole blood and genotyping for the four polymorphisms were done by ARMS-PCR, PCR-RFLP and gap-PCR. Gene-gene interactions were analyzed by logistic regression. In case of DRD2 A-241G, G allele was significantly associated with resistant to treatment. Regarding DRD4 120-bp duplication, 240/240 genotype was significantly associated with resistant to treatment comparing to other genotypes in a dominant model. The genotype combination of DRD4 240/240 and COMT Val/Val was significantly associated with treatment resistant. Among DRD2 AA genotype, COMT met allele carriers which also had a 120 bp allele of DRD4 had a significantly better response to antipsychotics. Moreover, analysis of clinical and demographic factors demonstrated a significantly longer duration of hospitalization and higher chlorpromazine-equivalent daily dose in resistant to treatment patients. Discovering the polymorphisms which effect treatment response to antipsychotics will provide the possibility of genetic screening before starting an antipsychotic treatment which enhances the chance of responding to antipsychotics and decreases drugs side effects and costs.

Peer Review reports


Schizophrenia (SCZ) is a chronic, hereditary and disabling neuropsychiatric disorder with a worldwide prevalence of approximately 1% [1, 2]. In the etiology of SCZ, genetic factors are thought to play an important role, the heritability currently ranges from 64 to 81%; although, genetic mechanisms remain unclear [2,3,4]. The mainstay of schizophrenia treatment has been antipsychotic drugs over the past 60 years; however, clinical response differs significantly between patients, with an overall response rate of 50–70% [1, 5,6,7]. Many of patients with schizophrenia discontinue or switch drug regimens due to lack of treatment efficacy and/or drugs adverse side effects. In this view, there is a great need for the identification of predictive clinical and biological markers of treatment consequence [8, 9]. Pharmacogenetic biomarkers focus to predict which patients could improve with specified drugs according to genetic variants. Thus, genotype-based customized drug treatments may allow optimizing the antipsychotic treatment, while helping to minimize drugs side effects [10, 11].

To date, pharmacogenomics studies of response to treatment in schizophrenia, have typically focused on genes encoding for drug targets, called pharmacodynamics related genes. Many of the research investigating the association of pharmacodynamic genes with antipsychotic treatment response have concentrated on dopaminergic pathways, one of the primary mechanisms of function of antipsychotics, especially the gene coding for the dopamine D2 receptor(DRD2 gene), which is a binding object for all available antipsychotic drugs. Dopaminergic gene SNPs are strongly related to drug sensitivity of antipsychotics; for example, several studies have indicated positive associations between DRD2 gene and antipsychotic response [10,11,12,13,14].

Lingyue Ma and et al. by a systematic review and meta-analysis indicated that for Asian patients, at rs1799978(A241G) in DRD2 gene AA genotype had a significantly greater improvement after risperidone therapy [15]. Three studies analyzed the association between rs1801028 in DRD2 gene and antipsychotic or risperidone treatment in China and Slovenia; findings were inconsistent between the two studies [15,16,17].

Dopamine receptor D4(DRD4) is one of the main targets of antipsychotic drugs; furthermore, DRD4 polymorphisms were associated with antipsychotic treatment response. A 120-bp duplication in the promoter region of DRD4 is hypothesized to influence the clinical response to antipsychotic drugs [18,19,20].

Catechol-O-methyltransferase (COMT) is a major enzyme that inactivates dopamine in the prefrontal cortex by enzymatic degradation. The efficacy of antipsychotic drugs in patients with schizophrenia has been indicated to be related to rs4680 SNP in COMT gene; associated with a decreased therapeutic effect, leading to pharmacotherapy resistance. In addition, the COMT interacts with several genes to produce various phenotypes [20,21,22,23,24].

Based on our search in Pubmed and Google Scholar, this is the first study examining the frequencies of DRD2 rs1799978, DRD2 rs1801028 and DRD4 120-bp duplication in Iranian population. Although, the frequency of rs4680 in COMT gene is assessed in several studies with Iranian ethnicity. According to the mentioned studies, the overall frequency of G allele of rs4680 in a total sample size of n = 2200 is 51.57% in Iranian population [25,26,27,28,29,30].

Consequently, there is insufficient evidence of associations between polymorphisms in dopaminergic genes and antipsychotic treatment, particularly in Iran; thus, we carried out an association study between clinical factors and four polymorphisms in DRD2, DRD4 and COMT genes, and antipsychotic treatment response in Iranian population. To best of our knowledge this is the first study in Iran that is investigating the effects of genetic factors and antipsychotic treatment response.

Materials and methods


One hundred and one Iranian patients with schizophrenia were recruited to the study from Salamat psychiatric hospital of Khuzestan and Derakhshesh Rehabilitation center of Ahvaz. Inclusion criteria were as follows: 1) clinically confirmed diagnosis of schizophrenia according to the Diagnostic and Statistical Manual 5th Edition DSM-V criteria.2) hospitalized patients with schizophrenia whose drug intake were under supervision.3) at least 6 weeks under controlled antipsychotic treatment with adequate dose. exclusion criteria were as follows: 1) psychiatric comorbidities)such as mood disorders)2)other mental disorders 3)severe side effect of previous antipsychotic treatment 4)poor compliance to treatment5)substance abuse 6) substance-related disorder. The ethics committee of Ahvaz Jundishapur University of Medical Sciences approved this study (IR.AJUMS.MEDICINE.REC.1399.048). Informed consent was received from all patients/legal representatives enrolled in this study.

Clinical and demographic information

Clinical and demographic information containing age, sex, marital status, education level, consumption of tobacco, start age of schizophrenia, duration of treatment and hospitalization, history of hospitalization and family history of psychiatric disorder were obtained. Furthermore, clinical symptoms before and after antipsychotic treatment, typical and atypical antipsychotics taken with daily doses were noted from patients medical files. The chlorpromazine-equivalent daily dose of antipsychotic drugs administered to patients was calculated based on published updated guidelines [31].

Response and resistance to treatment

Response to treatment was evaluated by Persian version of PANSS- Positive and Negative Syndrome Scale- and psychiatrist diagnosis of remission in schizophrenia symptoms. The PANSS assessment was carried out by two trained and accurate psychologists after at least 6 weeks of under controlled antipsychotic treatment. PANSS is an internationally accepted scale used to measure symptoms severity in patients with schizophrenia and quantify clinical signs of schizophrenia. This scale contained 30 items separated into 3 scales; 7 items for the positive scale, 7 items for the negative scale and 16 items for the general psychopathological scale. The final score that indicates the severity of the disorder; is calculated by summing the answers for each question. Resistance to treatment was defined when the PANSS score was \(\ge\)4 in at least two of categories: P2 (conceptual disorganization), P3 (hallucinatory behavior), P6 (suspiciousness/persecution), G9 (unusual thought content) [17, 32].

DNA sampling and genotyping

Genomic DNA was extracted from whole blood using the SinaClone DNA kit following the manufacturer instructions. Four polymorphisms in three candidate genes were analyzed in the present study: rs1799978 and rs1801028 in DRD2, rs4680 in COMT and 120-bp duplication polymorphism in DRD4.Genotyping of DRD2 polymorphisms were conducted by using the amplification-refractory mutation system polymerase chain reaction (ARMS-PCR). In the case of rs1799978, two forward primers and one common reverse primer were as follow: F(allele A): 5’-CAGCCTGCAATCACAGCTTA-3’, F(allele G): 5’-CAGCCTGCAATCACAGCTTG-3 and R:5’-TGAAGCTGGACAGCTCTGC-3’. The two forward primers and one common reverse primer for rs1801028 were as follow: F (allele G): 5’-TGACTCTCCCCGACCCGTC-3’, F(allele C): 5’- TGACTCTCCCCGACCCGTG-3’ and R: 5’-GTTTGCCCATCTGTAAAGTGAGCAC-3’(primers used from Zalina Zahari et al [33]). The ARMS-PCRs were performed in two 15ul reactions for each patient, which contained an initial denaturation at 95 °C for 2 min followed by 35–36 cycles of denaturing at 95 °C for 30 s, annealing at 65 °C for 30s and extending at 72 °C for 45s. After that, a final extension at 72 °C for 5 min was applied.

COMT rs4680 was analyzed by polymerase chain reaction-restriction fragment length polymorphism (PCR- RFLP). The PCR was conducted in a 15ul reaction system which included an initial denaturation at 95 °C for 2 min followed by 35 cycles of denaturing at 95 °C for 30 s, annealing at 61 °C for 30 s and extending at 72 °C for 30s. After 35 cycles, it experienced a final extension at 72 °C for 5 min. NIaIII restriction enzyme was used for digestion; as explained by Qianqian He and et al [34]. Genotyping of 120-bp duplication polymorphism in DRD4 was performed using the Gap-PCR (primers by [35]). The PCR was carried out in a15ul reaction system which contained a first denaturation at 95 °C for 2 min followed by 30 cycles of denaturing at 95 °C for 30 s, annealing at 66 °C for 30 s and extending at 72 °C for 55s. After 30 cycles, a final extension at 72 °C for 5 min was done.

Statistical analysis

All data analysis was conducted using Statistical Package for the Social Sciences SPSS, Version 26. Continuous variables were expressed in the form of mean ± standard deviation. Kolmogorov–Smirnov test was applied to check whether the data were normally distributed; which normally distributed and abnormally distributed data between two groups were calculated respectively by t-test and Mann-Whitney U test. In order to assess the association between categorical variables of clinical parameters and genotype associations chi-square test was performed. Associations between polymorphisms and antipsychotic treatment response were analyzed under five genetic models including homozygous model, heterozygous model, recessive model, dominant model and co-dominant model. Logistic regression analysis was carried out in order to analyze these five genetic models. Adjusted odds ratios, adjusted P-values and 95% confident interval (95%CI) were calculated. Furthermore, logistic regression analysis was used to examine the effect of interaction among polymorphisms on antipsychotic treatment response. P < 0.05 was considered as statistical significance in all tests.

Gpower software version was applied in order to analyze the power of logistic regression. According to Chen et al., 2010 [36] effect size in logistic regression can be classified as follows:

  • Odds Ratio < 1.68 - Very small.

  • 1.68 ≤ Odds Ratio < 3.47 - Small.

  • 3.47 ≤ Odds Ratio < 6.71 - Medium.

  • Odds Ratio ≥ 6.71 – Large.

The power analysis demonstrated that a sample size of 101 is sufficient to reveal a medium to large effect size with a minimum power of 80% at a significance level of 5%.


Clinical and demographic characteristics of patients

One hundred and one patients with schizophrenia were included in this study (75 men and 26 women), which among them 51 patients were classified in the treatment-responder group and 50 patients met the criteria for treatment-resistance. There were significant differences between two groups in marital status, smoking, duration of hospitalization, chlorpromazine-equivalent daily dose and total PANSS score(Table 1). With regarding to the clinical characteristics, significantly longer duration of hospitalization, higher PANSS score and also higher chlorpromazine-equivalent daily dose were observed in treatment-resistant group(Table 1). Furthermore, married and smoker patients were significantly more in the treatment-responder group comparing to the treatment-resistant group.

Table 1 clinical and demographic characteristics of patients

The frequently prescribed antipsychotics were perphenazine(47%), olanzapine(37%), risperidone(32%), quetiapine(26%) and haloperidol (24%).

Genetic analysis

Genotypic and allelic associations

The population was in Hardy-Weinberg equilibrium for the four polymorphisms genotyped in this study (p > 0.05). The G allele of DRD2 A-241G was associated with increased risk of resistant to treatment when compared to A allele (OR(95%CI): 3.661,P = 0.02,Table 2). As for DRD4 120-bp duplication, there existed a considerable difference in allele distribution(P = 0.064,Table 2), indicating a higher frequency of 240-bp allele in treatment-resistant patients; although, the P value was not significant.

Table 2 Genotype and allele frequencies of polymorphisms within DRD2, COMT and DRD4

The logistic regression analysis revealed that regarding DRD4 120-bp duplication, patients with 120/120 and 240/120 had a lower risk of developing resistant to treatment as compared to patients with 240/240 genotype(AOR(95%CI): 0.196, P value: 0.033, Table 3). Moreover, regarding DRD4 120-bp duplication, homozygous and heterozygous genetic models indicated relations with antipsychotic treatment response, however it did not reach the significance level(P = 0.055 and P = 0.053 Respectively.

Table 3 Associations between DRD2, COMT and DRD4 polymorphisms and antipsychotic treatment response under five genetic models by using chi-square and logistic regression for the AORs

Gene-gene interaction analysis

Whether the presence of three polymorphism’s genotypes could influence the risk for treatment resistant to antipsychotic drugs was determined between DRD4 120-bp duplication, COMT rs4680 and DRD2 A-241G. We carried out all possible subgroup analyses; which the significant interactions are indicated in Table 4. In the COMT Val/Val subset, we found significant association of the DRD4 genotype with antipsychotics treatment response; where the combination of COMT Val/Val genotype and DRD4 240/240 genotype had a high risk for developing treatment- resistance(OR(95%CI) = 3.232(1.056–9.892), P = 0.04). Also, among patients with COMT Val/Met - Met/Met genotypes(Met allele carriers) those whose genotypes where AA for DRD2 A-241G were significantly more likely to respond to antipsychotic drugs as compared to other genotype combinations(OR(95%CI) = 2.540(1.138–5.668), P = 0.023).

Table 4 Gene-gene Interaction analysis for 2 and 3 locus models by using logistic regression

Furthermore, analyzing the interactions of DRD2 A-241G and DRD4 120-bp duplication polymorphisms, revealed a significant association between DRD2 AA genotype and DRD4 120 bp allele carriers(DRD4 120/240 − 120/120), patients with this genotype combination had a significantly better respond to antipsychotics(OR(95%CI) = 3.000(1.279–7.035), P = 0.012).

Additionally, logistic regression analysis indicated a significant interaction among DRD2 A-241G, DRD4 120-bp duplication and COMT rs4680 polymorphism’s; where patients with AA − 120/240 or 120/120 - Val/Met or Met/Met showed a significantly better respond to antipsychotics when comparing to patients with GA or GG-240/240-Val/Val genotype(OR(95%CI) = 2.363(1.057–5.281), P = 0.036).


The key findings of the present study were as follows. First, our genetic analysis for DRD2 A-241G(rs1799978) polymorphism detected a significantly higher frequency of G allele in resistant to treatment patients in comparison with responders. A possible explanation for this association could be that since DRD2 binds to dopamine and is a G-protein coupled receptor, A-241G polymorphism is considered to be related to DRD2 density and affinity [37]. Furthermore, regarding DRD2, it is recorded that this receptor lonely could adjust effects of atypical antipsychotics; suggesting that DRD2 plays a substantial role in patients response to atypical antipsychotics [15]. In addition, Mingzhe Zhao et al. carried out a genome-wide and whole exome sequencing joint study analyzing both the common and rare genetic variants associated with risperidone treatment response in patients with schizophrenia. Significantly, after combining the findings of GWAS and WES they identified signaling pathways which were all associated with DRD2 and HTR2A, targets of risperidone [38]. Interestingly, consistent with our result, Escamilla R et al. in Mexico reported a significantly higher frequency of G allele in resistant to antipsychotic treatment patients [10]. Also, a meta-analysis reported a significant association between AA genotype and greater improvement after risperidone therapy in Asian patients [15]. However, inconsistently to our result, Yan P et al. studying 267 Han Chinese patients, documented a significantly superior response to olanzapine in G allele carriers comparing to wild AA patients with schizophrenia [37]. There are reasonable arguments to explain the disagreement observed in this study. First, diverse ethnicities as genetic background plays an important role in antipsychotic treatment response. Second, response phenotype is a complicated trait, differences observed in different studies may be related to heterogeneity of definitions or criteria of response phenotype [15, 39]. Third, different sample sizes may result in various outcomes.

Second, in the case of the Ser311Cys (rs1801028) polymorphism in DRD2 gene, Cys allele displayed a very low frequency in Iranian population, therefore we did not found any type of association related to this polymorphism. Consistent to our result, Terzić T et al. did not observe any association between DRD2 Ser311Cys polymorphism and antipsychotic treatment response in 138 Slovenian patients with schizophrenia [17]. Although, a study which was carried out in 690 Han Chinese patients, revealed a significant association between DRD2 Ser311Cys polymorphism and PANSS total improvement rates [16]. Since the genetic background is completely different in Han Chinese ethnicity compared to Iranian ethnicity and the sample size of the last study was larger than our study, the differences observed in these two studies could be explained.

Our third result, regarding DRD4 120-bp duplication polymorphism, documented a significant association between 240/240 genotype with resistant to treatment comparing to 120/240 and 120/120 genotypes. To the best of our knowledge this is the first time that this association is reported. The 120-bp duplication region which is located in the upstream of the first codon, includes common sequences for various transcriptional factors and the longer allele decreases transcription of DRD4 gene [20, 40]. Limited studies exist around the association of DRD4 120-bp duplication polymorphism with antipsychotic response; and a review study mentioned this polymorphism as a newly discovered marker of antipsychotic treatment response [41]. Rudi Hwang et al. reported a significant association between 120 bp allele and non-responder group in African Americans, although they did not observe this association in Caucasians [42]. Again the variety in ethnicities results in different outcomes.

Fourth, in case of COMT Val158Met (rs4680) polymorphism, we did not observe genotypic and allelic association with antipsychotic treatment response. Consistent to our result, Hajj A et al. studying 100 Lebanese patients, did not detect association between COMT Val158Met and antipsychotic treatment response in their whole sample. Although, they noted a gender-related difference for COMT SNP. In 27 men involved in their study, they reported that Met allele carriers were more prone to be resistant to treatment in comparison with men with Val/Val genotype [32].

Fifth, our results revealed three two-way and one three-way gene-gene interaction in DRD2 A-241G, COMT rs4680 and DRD4 120-bp duplication polymorphisms with antipsychotic treatment response. In case of these four gene-gene interactions, COMT Val158Met polymorphism was involved in three of them. First, we observed that in the COMT Val/Val subset, patients with DRD4 240/240 genotype had a significantly higher risk for resistant to treatment. To our knowledge this is the first time that this interaction is reported. Second, in the DRD2 A-241G AA homozygotes, COMT Met allele carriers demonstrated a significantly better response to antipsychotics. Third, our three way gene-gene interaction analysis indicated that DRD4 120-bp allele carriers and COMT Met allele carriers among DRD2 AA homozygotes had a significantly superior response to antipsychotic treatment. Interestingly, Rajagopal VM et al. observed that patients with COMT Val/Met or Met/Met genotype who also had DRD4 120-bp allele revealed a significantly better response to clozapine [20]. Additionally, a meta-analysis reported a significantly improved response in Met-allele carrier Asian and Caucasian patients [39]. COMT is a catabolic enzyme with functional connection to dopaminergic transmission, the first target of most antipsychotics [11, 39]. Protein function and structure studies indicate that COMT Val-Met substitution causes a disruption in enzyme stability and decreases enzymatic activity [11, 39]. Importantly, tonic-phasic dopamine hypothesis supports our findings; according to this hypothesis, the lower activity form of COMT enzyme encoded by Met-allele may increase tonic dopamine and reduce phasic dopamine subcortically. Thus, this hypothesis suggest that the Met allele would predict a greater response to atypical antipsychotics, as they are considered to increase tonic dopamine transmission [7, 43]. Finally, the fourth gene-gene interaction found in this study was between DRD2 A-241G and DR4 120-bp duplication polymorphisms, where the combination of DRD2 AA homozygotes and 120-bp allele carriers was found to be associated with an improved response to antipsychotics. To the best of our knowledge this is the first study analyzing the three way gene-gene interaction between DRD2 A-241G, DRD4 120-bp duplication and COMT Val158Met polymorphisms together. Also, the first to analyze the interaction between DRD2 A-241G with DRD4 120-bp duplication, and DRD2 A-241G with COMT Val158Met.

Sixth, regarding our clinical and demographic results, we demonstrated that marriage and smoking were significantly associated with better response to antipsychotics. Furthermore, duration of hospitalization, total PANSS score and chlorpromazine equivalent daily dose were significantly higher in resistant to treatment patients. Interestingly, Escamilla R et al. in Mexico reported significantly more weeks of disease evolution in ultra-resistant to treatment patients. Also, they observed the most prolonged durations of untreated psychosis and an earlier age of schizophrenia onset in ultra-resistant to treatment group. Besides, they observed more married patients in the responder group, however it did not reach the significance level [10]. Furthermore, consistent to our result, Aline Hajj et al. in Lebanon which studied 100 patients with schizophrenia, revealed a significantly higher chlorpromazine equivalent daily dose in resistant to treatment patients; either when the resistant to treatment was according to BRPS scale or when the PANSS scale was applied as the criteria of treatment resistant [32].Different ethnicities and sample sizes may be a reasonable explanation for diverse results obtained in these studies.

However, considerable work has been carried out in order to recognize valid biomarkers to predict antipsychotic treatment response in patients with schizophrenia, data acquired is still inadequate and with the exclusion of singular samples, no general guidelines regarding possible utilization of pharmacogenetics in antipsychotic consumption in daily clinical practice, exist. The exceptions are related to CYP2D6 and CYP3A4 enzymes. Clinical Pharmacogenetics Implementation Consortium and the Dutch Pharmacogenetics Working Group published guidelines recommending CYP2D6 genotyping for patients under treatment with risperidone; DPWG recommends a decreased dose of risperidone in patients with the poor metabolizer phenotype of CYP2D6 and altering treatment to another drug or titration of the dose in patients who are classified as CYP2D6 ultrarapid metabolizers. Again DPWG together with the FDA, in treatment with aripiprazole, recommends a reduced dose for patients who are classified as CYP2D6 poor metabolizers. Furthermore, in cases where aripiprazole and CYP3A4 inhibitors are taken simultaneously, the FDA suggests a decrease in aripiprazole dose, with no determination of the patient’s phenotype. Nevertheless, these guidelines are advantageous, they are in no way adequate, therefore, more and more emphasizing the need for further research in this field [44].

Our sample size was relatively small especially the subgroups in gene-gene interactions were small and a few genotypes had limited carriers(GG genotype of DRD2 A-241G polymorphism and GC of DRD2 rs1801028); also CC genotype of DRD2 rs1801028 was not observed in our sample. Furthermore, no multiple testing correction was used for P-values in gene-gene interaction analyses and we should mention the type I error possibility as a limitation. Besides, the studied population was only from one single ethnicity. Consequently, further investigations with larger sample sizes and meta-analyses, from various ethnicities analyzing several polymorphisms involved in pathways related to antipsychotics actions are warranted; in order to move toward personalized medicine in schizophrenia.


In summary, our results suggest that COMT, DRD2 and DRD4 genes together and DRD2 and DRD4 genes separately, may effect and predict antipsychotic treatment response in Iranian population. This kind of study may provide the possibility of genetic screening before starting a new antipsychotic trial, resulting in a better chance to achieve the most effective treatment for each patient in a shorter period of time, decreasing costs and minimizing adverse side effects of drugs.

Data availability

All data generated or analyzed during this study are included in this published article.


  1. Yan P, Gao B, Wang S, Wang S, Li J, Song M. Association of 5-HTR2A T102C and A-1438G polymorphisms with clinical response to atypical antipsychotic treatment in schizophrenia: a meta-analysis. Neurosci Lett. 2022;770:136395. Epub 2021 Dec 14. PMID: 34919991.

    Article  CAS  PubMed  Google Scholar 

  2. Liu Y, Fu X, Tang Z, Li C, Xu Y, Zhang F, Zhou D, Zhu C. Altered expression of the CSMD1 gene in the peripheral blood of patients with schizophrenia. BMC Psychiatry. 2019;19(1):113. PMID: 30987620; PMCID: PMC6466712.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  3. Fu X, Zhang G, Liu Y, Zhang L, Zhang F, Zhou C. Altered expression of the DISC1 gene in peripheral blood of patients with schizophrenia. BMC Med Genet. 2020;21(1):194. PMID: 33008326; PMCID: PMC7532617.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  4. Jauhar S, Johnstone M, McKenna PJ, Schizophrenia. Lancet. 2022;399(10323):473–486. PMID: 35093231.

  5. Zhang Y, Ren H, Wang Q, Deng W, Yue W, Yan H, Tan L, Chen Q, Yang G, Lu T, Wang L, Zhang F, Yang J, Li K, Lv L, Tan Q, Zhang H, Ma X, Yang F, Li L, Wang C, Zhang D, Zhao L, Wang H, Li X, Guo W, Hu X, Tian Y, Ma X, Li T. Chinese Antipsychotics Pharmacogenomics Consortium. Testing the role of genetic variation of the MC4R gene in Chinese population in antipsychotic-induced metabolic disturbance. Sci China Life Sci. 2019;62(4):535–43. Epub 2019 Mar 26. PMID: 30929193.

    Article  CAS  PubMed  Google Scholar 

  6. Han M, Deng C. BDNF as a pharmacogenetic target for antipsychotic treatment of schizophrenia. Neurosci Lett. 2020;726:133870. Epub 2018 Oct 9. PMID: 30312750.

    Article  CAS  PubMed  Google Scholar 

  7. Huang E, Zai CC, Lisoway A, Maciukiewicz M, Felsky D, Tiwari AK, Bishop JR, Ikeda M, Molero P, Ortuno F, Porcelli S, Samochowiec J, Mierzejewski P, Gao S, Crespo-Facorro B, Pelayo-Terán JM, Kaur H, Kukreti R, Meltzer HY, Lieberman JA, Potkin SG, Müller DJ, Kennedy JL. Catechol-O-Methyltransferase Val158Met polymorphism and clinical response to antipsychotic treatment in Schizophrenia and Schizo-Affective disorder patients: a Meta-analysis. Int J Neuropsychopharmacol. 2016;19(5):pyv132. PMID: 26745992; PMCID: PMC4886669.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  8. Maffioletti E, Valsecchi P, Minelli A, Magri C, Bonvicini C, Barlati S, Sacchetti E, Vita A, Gennarelli M. Association study between HTR2A rs6313 polymorphism and early response to risperidone and olanzapine in patients with schizophrenia. Drug Dev Res. 2020;81(6):754–61. Epub 2020 May 27. PMID: 32462699.

    Article  CAS  PubMed  Google Scholar 

  9. Wang J, Su P, Yang J, Xu L, Yuan A, Li C, Zhang T, Dong F, Zhou J, Samsom J, Wong AHC, Liu F. The D2R-DISC1 protein complex and associated proteins are altered in schizophrenia and normalized with antipsychotic treatment. J Psychiatry Neurosci. 2022;47(2):E134–47. PMID: 35361701; PMCID: PMC8979657.

    Article  PubMed  PubMed Central  Google Scholar 

  10. Escamilla R, Camarena B, Saracco-Alvarez R, Fresán A, Hernández S, Aguilar-García A. Association study between COMT, DRD2, and DRD3 gene variants and antipsychotic treatment response in Mexican patients with schizophrenia. Neuropsychiatr Dis Treat. 2018;14:2981–7. PMID: 30464483; PMCID: PMC6223330.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  11. Ye J, Ji F, Jiang D, Lin X, Chen G, Zhang W, Shan P, Zhang L, Zhuo C. Polymorphisms in dopaminergic genes in Schizophrenia and their implications in Motor deficits and antipsychotic treatment. Front Neurosci. 2019;13:355. PMID: 31057354; PMCID: PMC6479209.

    Article  PubMed  PubMed Central  Google Scholar 

  12. Islam F, Men X, Yoshida K, Zai CC, Müller DJ. Pharmacogenetics-guided advances in antipsychotic treatment. Clin Pharmacol Ther. 2021;110(3):582–8. Epub 2021 Aug 18. PMID: 34129738.

    Article  PubMed  Google Scholar 

  13. McClay JL, Adkins DE, Aberg K, Bukszár J, Khachane AN, Keefe RS, Perkins DO, McEvoy JP, Stroup TS, Vann RE, Beardsley PM, Lieberman JA, Sullivan PF, van den Oord EJ. Genome-wide pharmacogenomic study of neurocognition as an indicator of antipsychotic treatment response in schizophrenia. Neuropsychopharmacology. 2011;36(3):616–26. Epub 2010 Nov 24. PMID: 21107309; PMCID: PMC3055694.

    Article  CAS  PubMed  Google Scholar 

  14. Zhang JP, Robinson DG, Gallego JA, John M, Yu J, Addington J, Tohen M, Kane JM, Malhotra AK, Lencz T. Association of a Schizophrenia risk variant at the DRD2 locus with antipsychotic treatment response in First-Episode Psychosis. Schizophr Bull. 2015;41(6):1248–55. Epub 2015 Aug 28. PMID: 26320194; PMCID: PMC4601717.

    Article  PubMed  PubMed Central  Google Scholar 

  15. Ma L, Zhang X, Xiang Q, Zhou S, Zhao N, Xie Q, Zhao X, Zhou Y, Cui Y. Association between dopamine receptor gene polymorphisms and effects of risperidone treatment: A systematic review and meta-analysis. Basic Clin Pharmacol Toxicol. 2019;124(1):94–104. Epub 2018 Sep 11. PMID: 30103286.

  16. Han J, Li Y, Wang X. Potential link between genetic polymorphisms of catechol-O-methyltransferase and dopamine receptors and treatment efficacy of risperidone on schizophrenia. Neuropsychiatr Dis Treat. 2017;13:2935–43. PMID: 29255361; PMCID: PMC5722007.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  17. Terzić T, Kastelic M, Dolžan V, Plesničar BK. Genetic polymorphisms in dopaminergic system and treatment-resistant schizophrenia. Psychiatr Danub. 2016;28(2):127–31. PMID: 27287786.

    PubMed  Google Scholar 

  18. Taraskina AE, Nasyrova RF, Zabotina AM, Sosin DN, Sosina КА, Ershov EE, Grunina MN, Krupitsky EM. Potential diagnostic markers of olanzapine efficiency for acute psychosis: a focus on peripheral biogenic amines. BMC Psychiatry. 2017;17(1):394. PMID: 29221470; PMCID: PMC5723030.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  19. Zakharyan R, Ghazaryan H, Kocourkova L, Chavushyan A, Mkrtchyan A, Zizkova V, Arakelyan A, Petrek M. Association of Genetic Variants of Dopamine and serotonin in Schizophrenia. Arch Med Res. 2020;51(1):13–20. Epub 2020 Feb 18. PMID: 32086104.

    Article  CAS  PubMed  Google Scholar 

  20. Rajagopal VM, Rajkumar AP, Jacob KS, Jacob M. Gene-gene interaction between DRD4 and COMT modulates clinical response to clozapine in treatment-resistant schizophrenia. Pharmacogenet Genomics. 2018;28(1):31–5. PMID: 29087970.

    Article  CAS  PubMed  Google Scholar 

  21. Tybura P, Samochowiec A, Beszlej A, Grzywacz A, Mak M, Frydecka D, Bieńkowski P, Mierzejewski P, Potemkowski A, Samochowiec J. Some dopaminergic genes polymorphisms are not associated with response to antipsychotic drugs in patients with schizophrenia. Pharmacol Rep. 2012;64(3):528 – 35. PMID: 22814006.

  22. Vehof J, Burger H, Wilffert B, Al Hadithy A, Alizadeh BZ, Snieder H, GROUP investigators. ;. Clinical response to antipsychotic drug treatment: association study of polymorphisms in six candidate genes. Eur Neuropsychopharmacol. 2012;22(9):625 – 31. Epub 2012 Mar 3. PMID: 22386772.

  23. Zhuo C, Cheng L, Li G, Xu Y, Jing R, Li S, Zhang L, Lin X, Zhou C. COMT-Val158Met polymorphism modulates antipsychotic effects on auditory verbal hallucinations and temporal lobe gray matter volumes in healthy individuals-symptom relief accompanied by worrisome volume reductions. Brain Imaging Behav. 2020;14(5):1373–81. PMID: 30712251; PMCID: PMC7572342.

    Article  PubMed  Google Scholar 

  24. Werner FM, Coveñas R. Risk Genes in Schizophrenia and Their Importance in Choosing the Appropriate Antipsychotic Treatment. Curr Pharm Des. 2021;27(30):3281–3292. PMID: 33588721.

  25. Saravani R, Galavi HR, Lotfian Sargazi M. Catechol-O-Methyltransferase (COMT) gene (Val158Met) and brain-derived neurotropic factor (BDNF) (Val66Met) genes polymorphism in Schizophrenia: a case-control study. Iran J Psychiatry. 2017;12(4):265–70. PMID: 29472953; PMCID: PMC5816916.

    PubMed  PubMed Central  Google Scholar 

  26. Ahmadi L, Kazemi Nezhad SR, Behbahani P, Khajeddin N,Pourmehdi-Boroujeni. M. Genetic variations of DAOA (rs947267 and rs3918342) and COMT genes (rs165599 and rs4680) in Schizophrenia and Bipolar I disorder. Basic and Clinical Neuroscience. 2018;9(6):429–38.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  27. Torkaman-Boutorabi A, Ali Shahidi G, Choopani S, Reza Zarrindast M. Association of monoamine oxidase B and catechol-O-methyltransferase polymorphisms with sporadic Parkinson’s disease in an Iranian population. Folia Neuropathol. 2012;50(4):382-9. PMID: 23319194.

  28. Heydari M, Ghorbian S, Manizheh Sayyah Melli. and. The catechol-methyltransferase rs4680 G > A polymorphism is associated with uterine leiomyoma susceptibility. Gene, Cell and Tissue 6.1 (2019).

  29. تأیید et al. بررسی فراوانی پلی‌مورفیسم rs4680 ژن COMT در عملکرد رقابتی ورزشکاران نخبه‌ی رزمی‌کار زن و مرد ایرانی.مجله دانشکده پزشکی اصفهان41.705(2023):9–17.

  30. Azadmarzabadi E. Detection of six novel de novo mutations in individuals with low resilience to psychological stress. PLoS ONE. 2021;16(9):e0256285.

    Article  CAS  PubMed  PubMed Central  Google Scholar 


  32. Hajj A, Obeid S, Sahyoun S, Haddad C, Azar J, Rabbaa Khabbaz L, Hallit S. Clinical and Genetic Factors Associated with Resistance to treatment in patients with Schizophrenia: a case-control study. Int J Mol Sci. 2019;20(19):4753. PMID: 31557839; PMCID: PMC6801865.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  33. Zahari Z, Salleh MR, Zahri Johari MK, Musa N, Ismail R. A nested allele-specific multiplex polymerase chain reaction method for the detection of DRD2 polymorphisms. Malays J Med Sci. 2011;18(4):44–57. PMID: 22589672; PMCID: PMC3328935.

    PubMed  PubMed Central  Google Scholar 

  34. He Q, Shen Z, Ren L, Wang X, Qian M, Zhu J, Shen X. The association of catechol-O-methyltransferase (COMT) rs4680 polymorphisms and generalized anxiety disorder in the Chinese Han population. Int J Clin Exp Pathol. 2020;13(7):1712–9. PMID: 32782694; PMCID: PMC7414458.

    PubMed  PubMed Central  Google Scholar 

  35. Vereczkei A, Demetrovics Z, Szekely A, Sarkozy P, Antal P, Szilagyi A, Sasvari-Szekely M, Barta C. Multivariate analysis of dopaminergic gene variants as risk factors of heroin dependence. PLoS ONE. 2013;8(6):e66592. PMID: 23840506; PMCID: PMC3696122.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  36. Chen H, Cohen P, Chen S. How big is a big odds ratio? Interpreting the magnitudes of odds ratios in Epidemiological studies. Commun Stat - Simul Comput. 2010;39(4):860–4.

    Article  Google Scholar 

  37. Yan P, Song M, Gao B, Wang S, Wang S, Li J, Fang H, Wang C, Shi J. Association of the genetic polymorphisms of metabolizing enzymes, transporters, target receptors and their interactions with treatment response to olanzapine in Chinese Han patients with schizophrenia. Psychiatry Res. 2020;293:113470. Epub 2020 Sep 20. PMID: 32992097.

    Article  CAS  PubMed  Google Scholar 

  38. Zhao M, Ma J, Li M, Zhu W, Zhou W, Shen L, Wu H, Zhang N, Wu S, Fu C, Li X, Yang K, Tang T, Shen R, He L, Huai C, Qin S. Different responses to risperidone treatment in Schizophrenia: a multicenter genome-wide association and whole exome sequencing joint study. Transl Psychiatry. 2022;12(1):173. PMID: 35484098; PMCID: PMC9050705.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  39. Ma J, Zhao M, Zhou W, Li M, Huai C, Shen L, Wang T, Wu H, Zhang N, Zhang Z, He L. 2021. Association between the COMT Val158Met polymorphism and antipsychotic efficacy in schizophrenia: An updated meta-analysis. Current neuropharmacology, 19(10), p.1780.

  40. Sonia JA, Kabir T, Islam MMT, Kabir Y. Catechol-O-methyltransferase and dopamine receptor D4 gene variants: possible association with substance abuse in Bangladeshi male. PLoS ONE. 2021;16(2):e0246462. PMID: 33544778; PMCID: PMC7864466.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  41. Sun Y, Zhou W, Chen L, Huai C, Huang H, He L, Qin S. Omics in schizophrenia: current progress and future directions of antipsychotic treatments. J Bio-X Res. 2019;2(04):145–52.

    Google Scholar 

  42. Hwang R, Tiwari AK, Zai CC, Felsky D, Remington E, Wallace T, Tong RP, Souza RP, Oh G, Potkin SG, Lieberman JA. Dopamine D4 and D5 receptor gene variant effects on clozapine response in schizophrenia: replication and exploration. Prog Neuropsychopharmacol Biol Psychiatry. 2012;37(1):62–75.

    Article  CAS  PubMed  Google Scholar 

  43. Bilder R, Volavka J, Lachman H, et al. The Catechol-O-Methyltransferase polymorphism: relations to the Tonic–phasic dopamine hypothesis and neuropsychiatric phenotypes. Neuropsychopharmacol. 2004;29:1943–61.

    Article  CAS  Google Scholar 

  44. Płaza O, Gałecki P, Orzechowska A, Gałecka M, Sobolewska-Nowak J, Szulc A. Pharmacogenetics and Schizophrenia-Can Genomics improve the treatment with second-generation. Antipsychotics? Biomedicines. 2022;10(12):3165. PMID: 36551925; PMCID: PMC9775397.

    Article  CAS  PubMed  Google Scholar 

Download references


We would like to thank patients, their families and Salamat hospital personnel, especially Mrs Farokhmanesh for their co-operation in our study. We would like to thank Dr Morteza Mohammadzadeh from Iran University of Medical Sciences, for his assistance in a part of the statistical analysis. The study was supported by the Cellular and Molecular Research Center, Medical Basic Sciences Research Institute, Ahvaz Jundishapur University of Medical Sciences (CMRC-9726). The article is extracted from Narges Taheri thesis to obtain a master’s degree in genetics from Ahvaz Jundishapur University of Medical Sciences.


This work was supported by Cellular and Molecular Research Center, Medical Basic Sciences Research Institute, Ahvaz Jundishapur University of Medical Sciences, Ahvaz, Iran. [grant number CMRC-0001].

Author information

Authors and Affiliations



N.T. conducted sample collection, all of the genetic experiments and statistical analysis. Furthermore, N.T. wrote the manuscript text. P.G. was the leader of the study and was involved in the guidance of all of the study’s steps. Also, P.G. edited the manuscript text. M.B., M.S. and R.P. were involved in finding appropriate patients with the study’s inclusion criteria and providing the possibility of sample collection. In addition, R.P. was the advisor of all the psychiatry sections of the study. All authors read and approved the final manuscript.

Corresponding author

Correspondence to Pegah Ghandil.

Ethics declarations

Conflict of interest

We have no conflicts of interest to disclose.

Ethics approval and consent to participate

We confirm that all experiments were performed in accordance with relevant guidelines and regulations such as the Declaration of Helsinki.

Informed consent was received from all patients/legal representatives enrolled in this study.

The ethics committee of Ahvaz Jundishapur University of Medical Sciences approved this study (IR.AJUMS.MEDICINE.REC.1399.048).

Consent for publication

Not applicable.

Competing interests

The authors declare that they have no competing interests.

Additional information

Publisher’s Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article’s Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article’s Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit The Creative Commons Public Domain Dedication waiver ( applies to the data made available in this article, unless otherwise stated in a credit line to the data.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Taheri, N., Pirboveiri, R., Sayyah, M. et al. Association of DRD2, DRD4 and COMT genes variants and their gene-gene interactions with antipsychotic treatment response in patients with schizophrenia. BMC Psychiatry 23, 781 (2023).

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: