Rare Variants in Pharmacogenes Influence Clozapine Metabolism in Individuals with Schizophrenia

Introduction

Schizophrenia is a chronic psychiatric disorder affecting 7 out of 1000 people worldwide ¹. Between 20-30% of individuals with schizophrenia develop treatment-resistant schizophrenia (TRS), a particularly severe presentation of the disorder in which symptoms do not respond to treatment after at least two courses of standard antipsychotics^{2, 3}. Clozapine is the first-line treatment for TRS, but its use is associated with adverse effects which require careful clinical management⁴. In particular, due to the risk of a haematological adverse drug reaction (ADR) – agranulocytosis – individuals treated with clozapine undergo regular mandatory haematological monitoring. This monitoring can also be used to ensure treatment adherence by measuring concentrations of the drug and its metabolites in plasma. Therapeutic Drug Monitoring (TDM) frameworks have also been proposed to ensure clozapine dosing leads to plasma concentrations within the standard “therapeutic range” (350-600 ng/mL) that aims to balance the likelihood of treatment response with the risk of dose-dependent ADRs^{5, 6}. However, even in high-income countries with extensive public health systems, TDM is not always used by clinicians for several reasons including the low availability of testing, and a lack of familiarity with how to apply it in clinical decision-making⁷.

Drug pharmacokinetics can vary substantially between people, with both environmental and genetic factors partially explaining these differences. While 40% of the inter-individual variability in clozapine metabolism has been related to demographic, lifestyle and clinical predictors (sex, age, body weight and smoking)⁸, the few genetic studies performed to date have already found a small number of common genetic variants explaining up to an extra 3%-8%⁹. As the effects of these markers are additive to those of environmental predictors¹⁰, their discovery is a first step towards eventual pre-emptive pharmacogenomic (PGx) interventions to personalise and improve clozapine treatment protocols.

Most PGx studies of clozapine and other psychiatric medications have been based on the “candidate study” experimental design, investigating well-characterized common polymorphisms in genes encoding cytochrome P450 enzymes (CYPs) that in many cases have been mapped to haplotypes commonly known as “star alleles”¹¹. These PGx markers are routinely assessed by genetic testing companies and, given their strong effects on enzyme activity¹², they have been proposed as useful to guide drug selection and dosing. The hepatic first-pass metabolism of clozapine is driven by CYP1A2 (with a minor role played by CYP3A4 and CYP2D6), and therefore lifestyle factors and concomitant medications influencing the activity of this enzyme are commonly discussed in the literature⁶. However, a clear relationship between CYP1A2 genetic variation and clozapine metabolism has not yet been established, and no genetic variants are considered relevant for clinical management on recent expert reviews^{13, 14} or in the PharmGKB database¹⁵. The star allele CYP1A2*1F has been the most consistently studied polymorphism as it is indexed by a single common SNP (rs762551; global minor allele frequency [MAF] 38%¹⁶) and leads to ultra-rapid metabolism in conjunction with enzyme inducers¹⁷. Despite multiple candidate studies likely sparked by its proposed relevance for other metabolites and conditions¹⁸, there is little evidence for CYP1A2*1F association with clozapine pharmacokinetics or response^19-21. Genome-wide association studies (GWAS) of clozapine plasma concentrations have also not implicated variants within the CYP1A2 gene, although a 13kb upstream variant potentially involved in its regulation, rs2472297, has a reported and replicated association^{10, 22}. It is important to consider that both GWAS and candidate gene studies are best suited to identify the influences of relatively common alleles. Sequencing studies are required to detect rare genetic variation, which is responsible for most of the known variability of drug-related enzymes and transporters (“pharmacogenes”)^{23, 24}. Rare variants with a functional effect on pharmacogenes can thus become consequential for a given carrier when a relevant drug is prescribed, and thereby contribute to individual-level variability in drug metabolism and ultimately become clinically actionable through downstream effects on medication response, tolerability, and other outcomes^{17, 25}.

In this study, we sought to complement the growing literature on candidate studies and GWAS of clozapine metabolism by assessing whole exome sequencing (WES) data from over 2000 individuals taking this drug. Our aim was to investigate for the first time, whether rare variants (MAF ≤1%) in genes involved in the clozapine metabolism pathway affect clozapine pharmacokinetics.

Subjects and Methods

Results

After phenotypic and genetic quality control, 2062 individuals linked to 6585 clozapine plasma concentration assays (1-42 observations per individual) were retained: 72.8% were male (1501 male; 561 female), and the mean age was 43.61 years (SD= 11.38) at the last available clozapine monitoring measure. Over 6000 rare damaging nonsynonymous alleles (MAF ≤ 1%) affecting 313 of 321 genes included in the four gene sets were identified (624 PTVs; 5435 missense SNVs).

Gene set Association Analyses

We observed significant associations between the burden of rare damaging mutations in the PharmaADME core set and both clozapine and norclozapine plasma concentrations (Figure 1). The combined PTVs + missense variant set included 837 rare damaging alleles affecting 29 genes in the core PharmaADME set. Each rare allele counted in this gene set was associated with lower clozapine (β= -0.054, SE= 0.019, P-value= 0.005) and norclozapine levels (β= -0.043, SE= 0.018, P-value= 0.015) when compared to individuals carrying zero damaging alleles, though only the clozapine association remained significant after FDR correction. Through comparison with other predictors in the statistical model (Table 2), we estimate that the effects in clozapine plasma concentrations of a single damaging allele in this gene set are akin to reducing the clozapine dose by about 35 mg/day (Figure 2). Supplementary Table 2 details the main results for all gene sets for both clozapine and norclozapine levels, including the percentage of pharmacokinetic variance explained by rare variants in these sets.

• Download figure
• Open in new tab

Figure 1:

Estimated effect sizes and confidence intervals for the different variant classes in each of the four gene sets analysed. Supplementary Table 2 gives the full results for all gene sets for both clozapine and norclozapine levels, including the percentage of variability explained by rare variants in each of the sets.

• Download figure
• Open in new tab

Figure 2:

Marginal effects for individuals carrying different numbers of rare damaging variants in the PharmaADME core gene set in the relationship between clozapine doses and plasma concentrations. Dots and whiskers indicate the estimated value of the regression effect size and its corresponding 95% confidence interval at specific doses.

Leave-one-gene-out set association analyses

Rare alleles in CYP1A2 variants had an opposite direction of effect to rare alleles in the full PharmaADME core set (Figure 3 – Supplementary Table 5), that is, they were associated with higher clozapine and norclozapine levels. Thus, we performed a sensitivity analysis by excluding CYP1A2 from the PharmaADME set, recalculating the burden scores, and conducting the GLMM regression again. This yielded the same result of rare variant burden in the PharmaADME core set being associated with clozapine and norclozapine metabolism, and as expected the effect sizes were slightly stronger after the removal of CYP1A2 (Supplementary Figure 4). We also tested for additive effects of separately considering rare damaging variants in CYP1A2 and the remainder of the PharmaADME gene set, showing that the addition of both burdens as individual predictors resulted in an 80.83% increase in the variance of clozapine levels explained by rare genetic variation (Table 3). We did not find any evidence of multiplicative (i.e., non-linear) relationships between these two predictors in our dataset (interaction test P= 0.470).

• Download figure
• Open in new tab

Figure 3:

Marginal effects representing the relationship between clozapine doses and plasma concentrations in individuals carrying or not rare damaging variants in the CYP1A2 gene. Shaded areas highlight a 95% confidence interval for the estimates.

Discussion

To our knowledge, this is the largest sequencing study seeking to evaluate the effects of rare pharmacogenomic variants in clozapine metabolism. Using longitudinal clozapine monitoring data, we evaluated exome data obtained by WES to discover genes and gene sets harbouring rare variants associated pharmacokinetic variation. We provide evidence that deleterious variants in pharmacogenes can significantly affect clozapine metabolism, leading to different metabolite plasma concentrations for individuals prescribed the same drug dose.

Our primary analysis show that a greater burden of rare damaging variants in the PharmaADME core gene set, which contained 30 genes associated with general drug metabolism (Table 1), was associated with significantly lower clozapine and norclozapine plasma concentrations (Figure 1). Furthermore, the burden of rare variants in this set explains 0.40% of the variability of this phenotype (Table 3), a similar figure to the 0.61% variance explained by a PRS generated from genome-wide significant loci in a recent clozapine metabolism GWAS⁹. In the CLOZUK2 sample, the effects of each rare damaging allele in this set were approximately equivalent to a 35 mg reduction in clozapine daily dose. As our models assume a linear additive effect for the mutation burden, carriers of multiple rare alleles (the maximum seen in our study was 3) would require further dose adjustments to reach the expected levels of the wider population of non-carriers (Figure 2). Furthermore, through ordinal mixed-model regression analysis (Supplementary Results), we also observed that carriers of at least one damaging allele in the PharmaADME core set were 30% less likely than non-carriers in CLOZUK2 to reach the therapeutic range (350–600 ng/mL) of plasma concentrations of clozapine (OR= 0.694; SE= 0.145; P= 0.012; Supplementary Figure 5). Metabolic networks have properties that make them resistant to deactivation by deleterious genetic variation³⁹, but this also makes them more likely to tolerate protein-damaging variants. Our findings suggest that individuals accumulating several rare damaging variants (as assessed in the gene set burden scores) might then have impaired biochemical functions relevant to drug response or even downstream disease risk⁴⁰. However, we note this study is not powered to assess the strength of specific rare variants, which could also have different directions of effect for a given phenotype depending on their protein consequence.

In further analyses, we analysed the burden of rare variants within specific genes included in the PharmaADME core gene set. Although no single gene was significantly associated with clozapine or norclozapine plasma levels after Bonferroni correction, the most significant gene in this analysis was CYP1A2 (P= 0.009). However, surprisingly, individuals carrying rare variants in CYP1A2 showed higher plasma clozapine levels when taking the same clozapine dose as non-carriers (Figure 3). This was contrary to the PharmaADME core gene set results and prompted us to test for individual contributions to clozapine metabolism of rare variants within and outside CYP1A2 (Table 3). CYP1A2 encodes the main catalyst enzyme of the clozapine to norclozapine biotransformation, and while common variants in its proximity have been previously associated with clozapine metabolite levels^{27, 41}, no such association yet exists for coding SNPs. Our results suggest that CYP1A2 rare coding variants, not included in the common haplotypes and star alleles frequently assessed by candidate studies, might also influence clozapine metabolism. Further investigation will be needed to identify specific variants and map them to new PGx star alleles, if possible, as well as to understand if these results could potentially support initiatives for pharmacogenetic testing as a decision-support tool for clozapine clinics.

Although genetic factors that account for a significant proportion of the interindividual variability in antipsychotic drug metabolism have been suggested as useful to improve therapeutic response and avoid ADRs ⁴², an effective model for the clinical implementation of this information in guiding treatment choices is not yet available. This study supports the view that such a precision psychiatry framework might benefit from taking account of genetic variation beyond the traditional PGx star allele definitions assessed by commercial genetic testing companies^{23, 27}. Our analyses included rare damaging variants affecting pharmacogenes based on in-silico functional approaches, and thus most of them have not yet been functionally assessed (Supplementary Table 6). Nevertheless, and while specific variants should be externally validated before being evaluated further, the effects we detect for the mutational burdens are of a sufficient magnitude to be relevant for particularly vulnerable subsets of the population of those on clozapine, such as people initiating treatment (when ADRs are most common⁴³) or who only tolerate low drug doses. This knowledge, added to that coming from common variant analyses, might one day be incorporated into clozapine prescribing guidelines that at present rely on non-genetic characteristics to personalise dosing¹⁴. In the meantime, proactive therapeutic drug monitoring strategies remain a tried-and-tested approach for navigating the large inter-individual variability of clozapine metabolism and achieving beneficial clinical outcomes^{6, 44}.

One limitation of this study is that rare variant analyses require a large sample size to allow the identification of single rare variants that are associated with complex traits. Although we had access to WES data, we chose to focus our efforts on a relatively small number of genes with robust a priori evidence of their involvement in drug metabolism, attempting to avoid limiting the power of our statistical tests through multiple testing. This strategy followed the rationale and results from other genomic studies of drug metabolism, showing that these traits are likely influenced by a small number of variants in relevant genes, each encompassing a relatively large effect^{27, 45}. However, unlike candidate studies, which often have included variants with experimentally derived functional activity information, we selected variants based on existing protein impact prediction algorithms. This is a PGx-informed approach to variant prioritisation, selecting rare variants putatively affecting enzyme and transporter function^{34, 46}, but can still include markers with little or no impact on gene product while omitting other potentially important variants. Thus, assessing and collecting more information on the functional genetic makeup of pharmacogenes, through either bioinformatics or direct experimental approaches^{47, 48}, will be valuable to update these analyses and for further research into developing PGx biomarkers for metabolite pharmacokinetics and drug response.

Another limitation is that, as with most statistical testing of burden scores, the regression models fitted for gene- and set-based analyses in the GLMM framework assume that all rare variants under testing follow a similar direction of effects⁴⁹. This is known not to be the case in disease studies (i.e., both risk-increasing and protective variants can exist) and is likely a simplification for metabolism as well. While, again, our sample size limits the testing we can reliably perform for individual genes, results from gene-based analyses, particularly for CYP1A2, suggest that future studies might benefit from statistical methods allowing for testing variants with multiple directions of effect.

Finally, due to the limited data collection performed for CLOZUK2, important environmental factors, such as cigarette smoking or the use of other medications, were not available to be included in our analyses. These are known to be responsible for a fraction of the inter-individual differences in clozapine metabolism and thus might mediate some of our results. Nevertheless, some of the covariates already included in our regression models might capture some of the variance due to these exposures. For example, previous studies have found that sex can be a partial marker of oral contraceptive use and smoking, both interactors of CYP1A2, capturing some of their effects in drug metabolism^{50, 51}. A related issue is that while the CLOZUK2 cohort is derived from the diverse current population of the UK, the majority of individuals for which both WES and clozapine levels monitoring data were available were of European ancestry. Given that individuals of other biogeographical ancestries are known to have different clozapine pharmacokinetic profiles⁹, our study cannot assess a potential rare variant contribution to those known ancestral differences, and its results might not be straightforwardly transferrable to non-European populations. We note, however, that consistency of variant effect sizes across populations is a reasonable assumption for pharmacogenes⁵² and that the occurrence of deleterious variation in these gene sets is frequent (though individually rare) across continental ancestries⁵³. Thus, the association between rare damaging burden and clozapine metabolism identified here could be of relevance beyond the European sample in which it was discovered.

In conclusion, we found that the burden of rare damaging variants within genes in drug metabolism pathways was associated with differences in clozapine metabolite plasma concentrations, explaining part of the genetic variability in clozapine metabolism phenotypes. Our results suggest that pharmacogenomic efforts trying to predict clozapine metabolism and personalise drug therapy could benefit from the inclusion of rare damaging variants in pharmacogenes beyond those already identified and catalogued as PGx star alleles. Further investigation is warranted to confirm and replicate these findings before they can be evaluated in prospective trials or considered in clinical settings.

Conflict of Interest

AK is a full-time employee of Magna Laboratories Ltd. MH and JJ are full-time employees of Leyden Delta B.V. MJO, MCOD, and JTRW are supported by a collaborative research grant from Takeda Pharmaceuticals Ltd. for a project unrelated to work presented here. ER, AFP, MJO, MCOD, and JTRW also reported receiving grants from Akrivia Health for a project unrelated to this submission. Takeda and Akrivia played no part in the conception, design, implementation, or interpretation of this study.