Interaction testing and polygenic risk scoring to estimate the association of common genetic variants with treatment resistance in schizophrenia

Interaction Testing and Polygenic Risk Scoring to Estimate the Association of Common Genetic Variants With Treatment Resistance in Schizophrenia Antonio F. Pardiñas, PhD; Sophie E. Smart, PhD; Isabella R. Willcocks, MSc; Peter A. Holmans, PhD; Charlotte A. Dennison, BSc; Amy J. Lynham, PhD; Sophie E. Legge, PhD; Bernhard T. Baune, MD, PhD; Tim B. Bigdeli, PhD; Murray J. Cairns, PhD; Aiden Corvin, MD, PhD; Ayman H. Fanous, MD; Josef Frank, PhD; Brian Kelly, MD, PhD; Andrew McQuillin, PhD; Ingrid Melle, MD, PhD; Preben B. Mortensen, DrMedSc; Bryan J. Mowry, MD; Carlos N. Pato, MD, PhD; Sathish Periyasamy, PhD; Marcella Rietschel, MD; Dan Rujescu, MD, PhD; Carmen Simonsen, PhD; David St Clair, PhD; Paul Tooney, PhD; Jing Qin Wu, PhD; Ole A. Andreassen, MD, PhD; Kaarina Kowalec, PhD; Patrick F. Sullivan, MD, PhD; Robin M. Murray, DSc; Michael J. Owen, PhD; James H. MacCabe, PhD; Michael C. O’Donovan, PhD; James T. R. Walters, PhD; and the Genetics Workstream of the Schizophrenia Treatment Resistance and Therapeutic Advances (STRATA) Consortium and the Schizophrenia Working Group of the Psychiatric Genomics Consortium (PGC) IMPORTANCE About 20% to 30% of people with schizophrenia have psychotic symptoms that do not respond adequately to first-line antipsychotic treatment. This clinical presentation, chronic and highly disabling, is known as treatment-resistant schizophrenia (TRS). The causes of treatment resistance and their relationships with causes underlying schizophrenia are largely unknown. Adequately powered genetic studies of TRS are scarce because of the difficulty in collecting data from well-characterized TRS cohorts. OBJECTIVE To examine the genetic architecture of TRS through the reassessment of genetic data from schizophrenia studies and its validation in carefully ascertained clinical samples. DESIGN, SETTING, AND PARTICIPANTS Two case-control genome-wide association studies (GWASs) of schizophrenia were performed in which the case samples were defined as individuals with TRS (n = 10 501) and individuals with non-TRS (n = 20 325). The differences in effect sizes for allelic associations were then determined between both studies, the reasoning being such differences reflect treatment resistance instead of schizophrenia. Genotype data were retrieved from the CLOZUK and Psychiatric Genomics Consortium (PGC) schizophrenia studies. The output was validated using polygenic risk score (PRS) profiling of 2 independent schizophrenia cohorts with TRS and non-TRS: a prevalence sample with 817 individuals (Cardiff Cognition in Schizophrenia [CardiffCOGS]) and an incidence sample with 563 individuals (Genetics Workstream of the Schizophrenia Treatment Resistance and Therapeutic Advances [STRATA-G]). MAIN OUTCOMES AND MEASURES GWAS of treatment resistance in schizophrenia. The results of the GWAS were compared with complex polygenic traits through a genetic correlation approach and were used for PRS analysis on the independent validation cohorts using the same TRS definition. RESULTS The study included a total of 85 490 participants (48635 [56.9%] male) in its GWAS stage and 1380 participants (859 [62.2%] male) in its PRS validation stage. Treatment resistance in schizophrenia emerged as a polygenic trait with detectable heritability (1% to 4%), and several traits related to intelligence and cognition were found to be genetically correlated with it (genetic correlation, 0.41-0.69). PRS analysis in the CardiffCOGS prevalence sample showed a positive association between TRS and a history of taking clozapine (r 2 =2.03%; P= .001), which was replicated in the STRATA-G incidence sample (r 2 = 1.09%; P= .04). CONCLUSIONS AND RELEVANCE In this GWAS, common genetic variants were differentially associated with TRS, and these associations may have been obscured through the amalgamation of large GWAS samples in previous studies of broadly defined schizophrenia. Findings of this study suggest the validity of meta-analytic approaches for studies on patient outcomes, including treatment resistance. JAMA Psychiatry. 2022;79(3):260-269. doi:10.1001/jamapsychiatry.2021.3799 Published online January 12, 2022. Supplemental content Author Affiliations: Author affiliations are listed at the end of this article. Group Information: Authors and collaborators in the Genetics Workstream of the Schizophrenia Treatment Resistance and Therapeutic Advances (STRATA) Consortium and collaborators in the Schizophrenia Working Group of the Psychiatric Genomics Consortium (PGC) are listed at the end of the article. Corresponding Authors: Antonio F. Pardiñas, PhD (pardinasa@cardiff. ac.uk), and James T. R. Walters, PhD ([email protected]), MRC Centre for Neuropsychiatric Genetics and Genomics, Division of Psychological Medicine and Clinical Neurosciences, Cardiff University, Maindy Road, Hadyn Ellis Building, Cardiff CF24 4HQ, United Kingdom. Research JAMA Psychiatry | Original Investigation 260 (Reprinted) jamapsychiatry.com Downloaded From: https://jamanetwork.com/ by a Universidad de Sevilla User on 11/22/2022 Precision psychiatry provides a potential pathway to improve psychiatric classification and develop treatments that are better tailored to specific patients. Major advances in determining the role of genetic variation in the risk of developing psychiatric disorders are helping realize this potential, 1 but the relevance of these findings to patient outcomes and response to treatment is unclear. Evidence from nonpsychiatric disorders points to distinct genetic bases for disorder susceptibility and prognosis, 2,3 although few sufficiently powered psychiatric genetic studies have been conducted investigating disorder progression or outcome. 4 Such study designs could offer valuable insights as to the feasibility of precision psychiatry approaches and their potential impact for patients who have poor outcomes because of the lack of effectiveness of current treatments. 5 Some of the most disadvantaged individuals in this respect are those with a diagnosis of schizophrenia whose symptoms do not respond adequately to conventional antipsychotic medication. 6-8 This clinical picture is known as treatment-resistant schizophrenia (TRS), which affects approximately 20% to 30% of people with this disorder. The only licensed treatment for TRS is clozapine, which is effective in approximately 60% of cases 9 and improves most indicators of morbidity and mortality. 10,11 While the mechanism of action of clozapine is still not fully understood, 12 it has been proposed that its efficacy might be related to the biological underpinnings of treatment resistance, suggesting a distinct neurobiological etiology to that of non-TRS. 13 Additionally, several studies have shown that delay in clozapine prescription is associated with resistance even to clozapine. 14,15 This makes early identification of TRS critical and the ascertainment of correlates of TRS a priority for the field of schizophrenia research. To date, there is considerable heterogeneity in the findings of genetic studies related to TRS. A family history of schizophrenia is likely associated with developing TRS, 16 but a recentsystematicreview didnot findanyindividualgenesrobustly and specifically associated with this condition. 17 Researchers have also investigated the long-standing hypothesis that TRS isa more severeformof schizophreniaresultingfroma largeburden of schizophrenia risk alleles. 18 While 2 studies reported that aggregating risk variants into a polygenic risk score (PRS) revealedsmalldifferences betweenTRSandnon-TRSsamples, 19,20 other analyses did not replicate this result. 16,21-23 However, because each study used slightly different recruitment criteria and definitions of TRS, heterogeneity in the results is expected, 24 making their collective interpretation difficult. This limitation likely reflects the previously discussed problem of individuals with treatment-resistant psychiatric symptoms being disproportionatelyunderrepresentedin research studiesowing topoor health, limited capacity to consent, and other causes of attrition, including therapeutic pessimism, 25 which have precluded the execution of well-powered genome-wide association studies (GWASs). 26 Inthis study, weaim tocharacterizethe contributionofcommon genetic variants to treatment resistance in schizophrenia by exploiting data generated by large consortium-based GWASs of schizophrenia, 27,28 in which individuals with TRS can be formally defined. Our main hypotheses are that schizophrenia risk alleles will show different genetic associations in the analysis of individuals with and without TRS and that these differences reflect the underlying genetics of treatment resistance. When assessed at a genome-wide level, such differences could reveal commonalities with other complex traits or be validated in clinical cohorts using the PRS approach, 29 helping us to better understand the biological and epidemiological characteristics of treatment resistance in schizophrenia. 17 Methods For our main analysis, we used samples from large genomic studies of schizophrenia. 27,28 Some of the largest cohorts in these research initiatives were recruited based on clozapine prescription (a proxy of TRS status), and forming a case-case data set from them would require avoiding confounding factors, such as GWAS batch effects 30 or population stratification, 31 which are difficult to control in a multiple-cohort design. 32 As a safeguard against these, we have used a meta-analytic procedure to assess the differences between GWAS in which individuals with TRS and non-TRS have been compared with matched sets of healthy controls, before comparing the allelic association effect sizes of these 2 GWASs on a genomewide basis to create a GWAS specific to treatment resistance. Genetic Samples and Analysis We used the CLOZUK1 and CLOZUK2 cohorts as our primary source of individuals with treatment-resistant symptoms, with a total sample size of 10 501 individuals with TRS and 24 542 controls. These cohorts have been described in previous studies. 28,33 All the individuals with TRS in these samples were prescribed clozapine in the UK after failure of at least 2 trials of antipsychotics, following National Institute for Health and Care Excellence guidelines for TRS. 34 The use of a history of taking clozapine as equivalent to a research diagnosis of TRS has been validated in these samples, 28 as well as in independent studies. 35,36 Control individuals were collected from public databases or through collaboration with population sequencing projects in the UK. To identify individuals with non-TRS, we used 34 studies included in the meta-analysis by the Schizophrenia Key Points Question Can common genetic variants be used to differentiate between treatment-resistant schizophrenia (TRS) and other forms of this disorder? Findings Data from this genome-wide association study including 85 490 participants were used to estimate genome-wide single-nucleotide variation effect size differences between individuals with and without TRS, which were compatible with a polygenic model of treatment resistance. Results were used to generate a polygenic risk score, which was significantly associated with TRS status in independent incidence and prevalence samples. Meaning Findings of this study based on common genetic variants indicate that TRS is heritable with a modest but significant single-nucleotide variation–based heritability. Contribution of Common Genetic Variants to Treatment Resistance in Schizophrenia Original Investigation Research jamapsychiatry.com (Reprinted) JAMA Psychiatry March 2022 Volume 79, Number 3 261 Downloaded From: https://jamanetwork.com/ by a Universidad de Sevilla User on 11/22/2022 Working Group of the Psychiatric Genomics Consortium (PGC), 27 with a total sample size of 20 325 individuals with schizophrenia and 30 122 controls. In 14 of these studies, available clinical records allowed us to identify and remove all individuals with TRS (eMethods in Supplement 1). The remaining 20 studies were not screened as comparable data were not available, and thus, we conservatively included these samples in the analysis as non–treatment-resistant cases (eTable 1 in Supplement 1). Healthy control individuals in this analysis were a mixture of publicly available samples and clinically ascertained (nonpsychiatric) individuals, and extensive analyses to discard population outliers and assess stratification were also carried out. 27 GWASs of individuals with TRS vs controls and individuals with non-TRS vs controls were carried out separately using the logistic regression model implemented in PLINK version 2. 37 Further details on the imputation, quality control, and association testing procedures are described in the eMethods in Supplement 1. All studies used in the GWAS were reviewed and approved by their local ethical committee. Written informed consent (or legal guardian consent and participant assent) was obtained for all study participants except for those in the CLOZUK cohort, as this study used anonymized blood samples as approved by the UK Multicenter Research Ethics Committee. Statistical Analysis Comparing the TRS and Non-TRS Association Studies To generate association statistics that reflect TRS vs non-TRS differences, we used the test for interaction proposed by Altman and Bland, 38 which is analogous to a fixed-effect test for moderators in the meta-analytic setting. 39 This test allows us to calculate an estimate of the difference between 2 odds ratios (ORs), which in our case were those of individuals with TRS vs controls (OR 1 ) and those with non-TRS vs controls (OR 2 ). The difference, in the scale of the regression β coefficient, equates to difference = log(OR 1 ) − log(OR 2 ), with standard error SE(d) = 兹(SE[OR 1 ] 2 + SE[OR 2 ] 2 ). We transformed this effect size into a zscore and calculated its associated Pvalue at each overlapping single-nucleotide variation (SNV; formerly single-nucleotide polymorphism) between the TRS and non-TRS GWAS, excluding 146 SNVs with a minor allele frequency difference greater than 20% between both data sets to mirror preimputation quality control. For consistency, we also used the same definition of OR 1 and OR 2 in all our tests to preserve directions of effect, and thus, positive zscores in the interaction analysis reflect stronger SNV associations in individuals with TRS compared with non-TRS. Estimating Heritability SNV-based heritability was estimated from the TRS interaction summary statistics using LD SCore version 1.01, 40 after restricting to markers present in the HapMap3 study 41 and a precomputed linkage disequilibrium reference panel based on the 1000 Genomes phase 3 samples. As point estimates of SNV-based heritability can differ based on model assumptions, we also estimated this quantity via the SumHer framework implemented in LDAK5, 42 using a default parameter to define the relationship between minor allele frequency and effect size (α = −.25) and multiplicative inflation correction (genomic control). The LD reference for LDAK was estimated from the European samples of the public Haplotype Reference Consortium panel. 43 Polygenic Validation We sought to investigate whether a TRS GWAS PRS was associated with treatment resistance status in 2 independent samples: one from a schizophrenia prevalence cohort (Cardiff Cognition in Schizophrenia [CardiffCOGS]) and the other from a first-episode incidence cohort (Genetics Workstream of the Schizophrenia Treatment Resistance and Therapeutic Advances [STRATA-G]). The CardiffCOGS prevalence cohort is a sample of individuals with schizophrenia that included both individuals with TRS and treatment-responsive schizophrenia, defined here again from a history of taking or not taking clozapine (n = 817; 315 with TRS and 502 with non-TRS). This is a cross-sectional study with detailed clinical ratings based on research diagnostic interviews and contemporaneous records, as previously described. 44 We also used a new multiancestry incidence cohort of people with first-episode psychosis, the STRATA-G consortium, in which the participants have been followed up for at least 1 year after initial presentation to ascertain diagnosis and treatment response (eMethods in Supplement 1). As both the first-episode and broad psychosis nature of this sample likely make it diagnostically heterogeneous, we first restricted our analyses to individuals with a diagnosis of schizophrenia at the last time of follow-up. TRS was rated as the presence or absence of a history of taking clozapine (n = 562; 71 with TRS and 492 with non-TRS) to best mirror our TRS GWAS and the definition used in the prevalence cohort. Further details about the recruitment and phenotyping of the STRATA-G sample as well as details of the genotyping and imputation procedure for both data sets are given in the eMethods and eTable 2 in Supplement 1. The CardiffCOGS study had the relevant UK National Health Serviceethical approval, andstudies includedin STRATA-G were reviewed and approved by their local ethical committees. Written informed consent was obtained for all participants. PRS Analysis Our primary PRS analysis used the TRS interaction GWAS summary statistics as training set, and the prevalence and incidence cohorts as testing sets where their treatmentresistance phenotype was tested for association. We also investigated the association of PRS and treatment resistance status based on (1) the CLOZUK TRS vs healthy control GWAS and (2) the PGC non-TRS GWAS. Polygenic scores were estimated using PRSice-2 version 2.35 45 and PRS-CS version June 4, 2021, 46 as detailed in the eMethods in Supplement 1. Within each pairing of training and testing set, PRS association Pvalues were corrected for multiple testing using the Benjamini and Hochberg false discovery rate (FDR) method 47 as implemented in the R framework stats package. Estimating Genetic Correlation and Colocalization Correlations between treatment resistance and other complex traits were computed using LD-Hub version 1.93, 48 again Research Original Investigation Contribution of Common Genetic Variants to Treatment Resistance in Schizophrenia 262 JAMA Psychiatry March 2022 Volume 79, Number 3 (Reprinted) jamapsychiatry.com Downloaded From: https://jamanetwork.com/ by a Universidad de Sevilla User on 11/22/2022 based on the LD-Score regression 40 framework. We only performed the analysis for GWAS categories that had significant phenome-wide associations with psychiatric disorders, based on previous research. 49 The categories chosen were education, cognitive, personality, psychiatric, and smoking, for a total of 28 summary statistics. Multiple testing correction of these results was also carried out using FDR, since most of the summary statistics tested have overlapping samples and thus yield partially dependent Pvalues. For all associations, statistical significance was set at P< .05. All Pvalues were 2-tailed. Results TRS/Non-TRS GWAS The study included a total of 85490 participants (48635 [56.9%] male) in its GWAS stage and 1380 participants (859 [62.2%] male) in its PRS validation stage. Discovery GWAS of individuals with TRS and non-TRS against independent controls showed highly consistent genome-wide significant signals (eFigure 1 in Supplement 1), all reported in previous schizophrenia studies. This consistency extended to most of the polygenic architecture, as shown by a high genetic correlation estimate between the data sets (r = 0.966; SE, 0.037; P=1.43×10 −147 ) in which the 95% CI includes 1. No genomewide SNVs were detected by the TRS interaction analysis (Figure 1), although a slight departure from normality was detected in the genome-wide Q-Q plot (λ = 1.062; λ 1000 = 1.002; LD-Score intercept, 1.032; SE, 0.007; LDAK scaling estimate, 1.017; SE, 0.010; eFigure 2 in Supplement 1). This is compatible with a polygenic signal, 50 as the 2 discovery GWASs did not show evidence of inflation caused by other confounding factors. 27,28 Indeed, the 2 methods we used to estimate the common variant contribution to the observed-scale heritability of TRS returned small values of similar magnitude (SNV- based heritability: LD-Score, 0.013; SE, 0.006; SNV-based heritability: LDAK, 0.040; SE, 0.014). Polygenic Score–Based Prediction of TRS in Independent Samples A PRS generated from the TRS interaction GWAS was positively associatedwithtreatmentresistance inour validationprevalence cohort (CardiffCOGS), explaining up to 2.03% of the variance on the liability scale of treatment resistance (P= .001; FDR P= .01; Figure 2;eTable3inSupplement 1). Analyzing the 2 GWASs used in the interaction analysis, we found that a PRS based on CLOZUK GWAS results (TRS vs healthy controls) also explained a significant proportion of the variance in the TRS phenotype (maximumr 2 = 1.63%; P= .004;FDRP= .03), very similarto that obtained by testing the larger PGC non-TRS sample (maximum r 2 =1.21%;P= .01;FDR P= .10). However, the effect sizes of these PRSs were generally in opposite directions (Figure 2; eTable 3 in Supplement 1); ie, those taking clozapine in CardiffCOGS had higher PRS derived from the TRS GWAS (OR, 1.22; 95% CI, 1.05- 1.41) but a lower PRS from the non-TRS GWAS (OR, 0.83; 95% CI, 0.72-0.96). Focusing on the polygenic prediction via Bayesian regression and continuous shrinkage priors estimates, which do not require LD clumping or Pvalue thresholding, 51 we also observed a significant and positive difference between the effect sizes of the TRS interaction and CLOZUK PRS, while the PGC PRS had a null effect (Figure 2). This pattern of results in an independent case-case comparison is consistent with our interaction analysis detecting a true polygenic signal for treatment resistance, with a greater burden of risk alleles from both this and the CLOZUK TRS GWAS associating with a history of clozapine prescription. Alternatively, the PRS derived from non-TRS samples is not as clearly associated with treatment resistance, and some resultsof theseanalyses areevencompatiblewithitsenrichment in individuals with non-TRS. Generation of PRS and their analysis within the firstepisode incidence sample (STRATA-G) showed that a greater Figure 1. Manhattan Plot of the Treatment-Resistant Schizophrenia (TRS) vs Non-TRS Interaction Genome-Wide Association Study 7 5 6 4 3 2 1 0 –log10 (P value) Chromosome 121 2 3 4 5 6 7 8 9 10 11 13 14 15 16 17 18 19 20 21 X22 Contribution of Common Genetic Variants to Treatment Resistance in Schizophrenia Original Investigation Research jamapsychiatry.com (Reprinted) JAMA Psychiatry March 2022 Volume 79, Number 3 263 Downloaded From: https://jamanetwork.com/ by a Universidad de Sevilla User on 11/22/2022 interaction PRS burden was also associated with having treatment resistance in this setting, explaining up to 1.09% of the variance (P= .04; FDR P= .21). In contrast, neither the CLOZUK nor PGC PRS were significantly associated with treatment resistance across any of the tested Pvalue thresholds or using polygenic prediction via Bayesian regression and continuous shrinkage priors (Figure 2; eTable 3 in Supplement 1). However, a meta-analysis of both CardiffCOGS and STRATA-G resulted in narrower confidence intervals around all the PRS effect size estimates previously estimated from CardiffCOGS, supporting the overall consistency in directions of association between these samples (eFigure 3 and eTable 3 in Supplement 1). To interpret the results from the prevalence and incidence cohorts alongside each other, we merged the genotyped SNVs from both samples (to a total of 199 074 SNVs) and estimated polygenic scores using all LD-independent SNVs (P< 1). As a reference point for this analysis, we retained a small subset of 242 screened unaffected control individuals that had been genotyped as part of the Genetics and Psychosis project, 26 one of the STRATA-G cohorts. Splitting the combined sample by TRS status showed that the polygenic profile of treatment resistance is largely consistent between CardiffCOGS and STRATA-G, with the TRS interaction PRS burden showing no difference between incident and prevalent TRS individuals (Figure 3). However, we noted a distinction in terms of general schizophrenia risk alleles, as indexed by the CLOZUK and PGC PRS, which were enriched within the individuals with firstepisode TRS compared with the individuals with TRS from our cross-sectional sample (CLOZUK: OR, 1.52; SE, 0.21; P= .047; PGC: OR, 3.48; SE, 0.39; P=7.69×10 −4 ). None of the scores we tested were significantly different between the 2 groups of individuals with non-TRS, and all individuals with schizophrenia showed larger mean scores than the unaffected controls. Genetic Correlation Between TRS and Other Traits We used the results of the TRS interaction GWAS to examine the genetic associations between the treatment-resistant phenotype and other disorders and traits. In these analyses, 8 of 28 publicly available GWAS summary statistics displayed nominally significant (P< .05) genetic correlations with TRS, of which 5 survived multiple testing correction (FDER p < 0.05; Figure 4;eTable4inSupplement 1). All 5 genetically correlated phenotypes were associated with cognitive measures Figure 2. Polygenic Score Analysis of the Cardiff Cognition in Schizophrenia (CardiffCOGS) Cohort and Genetics Workstream of the Schizophrenia Treatment Resistance and Therapeutic Advances (STRATA-G) Cohort Using 3 Different Schizophrenia-Related Polygenic Risk Scores (PRS) 0.20 0.15 0.10 0.05 0 Liability scale R2 CardiffCOGS A Interaction (TRS vs non-TRS) CLOZUK (TRS vs controls) PGC (non-TRS vs controls) 1.70 1.25 1.00 0.80 0.60 Standardized PRS odds ratio 1.70 1.25 1.00 0.80 0.60 Standardized PRS odds ratio 0.20 0.15 0.10 0.05 0 Liability scale R2 STRATA-G B Interaction (TRS vs non-TRS) CLOZUK (TRS vs controls) PGC (non-TRS vs controls) P value threshold 1 × 10–5 1 × 10–4 .001 .01 .05 .10 .50 >.99 PRS-CS a a a a a a a a a a a Associations between PRS and treatment resistance in schizophrenia (TRS) were defined as a history of taking clozapine in people with a schizophrenia diagnosis. There were 315 individuals with TRS and 502 individuals with non-TRS in the CardiffCOGS cohort and 71 individuals with TRS and 492 control individuals in the STRATA-G cohort. PGC indicates Psychiatric Genomics Consortium; PRS-CS, polygenic prediction via Bayesian regression and continuous shrinkage priors. Whiskers indicate 95% CIs. a P< .05. Research Original Investigation Contribution of Common Genetic Variants to Treatment Resistance in Schizophrenia 264 JAMA Psychiatry March 2022 Volume 79, Number 3 (Reprinted) jamapsychiatry.com Downloaded From: https://jamanetwork.com/ by a Universidad de Sevilla User on 11/22/2022 and educational attainment (genetic r, 0.41-0.69), showing our GWAS of TRS to be genetically correlated with lower cognitive ability and lower educational attainment. Discussion We report a large study of treatment resistance in schizophrenia, showing that common genetic variants were associated with this presentation, that these variants were also associated with general risk of schizophrenia, and that the difference in their associations with the 2 phenotypes may have influencedthe polygenic scoreanalyses. Furthermore,our results, particularly the large genetic correlation between the TRS and non-TRS GWASs, are consistent with the risk conferred by common SNVs being similar regardless of whether individuals respond to first-line treatments. This suggests that treatment resistance, at least to the extent that we have been able to define it and explore it within this study, is not likely to index a cluster of individuals with a range of etiologies and pathophysiologies fundamentally different from those that contribute to schizophrenia more widely. However, by also assessing differences at the level of individual allelic associations, our results also support the existence of a polygenic contribution associated with treatment resistance that seems largely distinct from liability to schizophrenia. While this approach has not identified specific SNVs or genes that could be followed up as potential drivers of treatment outcomes, we find genomewide correlations that recapitulate previous epidemiological findings, such as the association of treatment resistance with poorer cognitive performance. 52,53 Such polygenic overlaps validate our indirect approach for carrying out the GWAS of TRS in the absence of a single large-scale harmonized case-case sample.However, they should notbe interpreted asproof of causality in either direction and do not account for confounders, such as sociodemographic indicators, medication adherence, or antipsychotic effects in treatment response. Implications In revealing the existence of a common, heritable genetic signal for TRS, our study adds a new layer of evidence to the ongoing discussion of whether TRS is categorically distinct from treatment-responsive schizophrenia. 17 Despite a large genetic correlation between both conditions, we show that PRS derived from them perform differently when attempting out-of-sample prediction of TRS, which suggests that the heterogeneous results previously obtained in PRS analyses might have been influenced by the presence of individuals with TRS in the generic schizophrenia training sample. 19 While follow-up studies would likely require rich clinical data on large numbers of individuals with TRS and non-TRS to fully explore the implications of these results, they pose an important consideration for research seeking to understand whether Figure 4. Linkage Disequilibrium Score Regression Genetic Correlation Results of LD-HUB Phenotypes With the Treatment-Resistant Schizophrenia Interaction Genome-Wide Association Study Genetic correlation –1.0 1.00.50–0.5 –log10 (P value) 0 54321 Educational attainment Cognitive performance Educational attainment University degree Intelligence Former vs current smoker Cigarettes per day Ever vs never smoker Educational attainment Cognitive performance Educational attainment University degree Intelligence Former vs current smoker Cigarettes per day Ever vs never smoker Top, Genetic correlation coefficients shown with standard errors. Bottom, Genetic correlation Pvalues. The dashed line indicates a false discovery rate–corrected significant level of P= .05. Figure 3. Polygenic Profile Derived From Combined Genotype Data From the Cardiff Cognition in Schizophrenia (CardiffCOGS) and Genetics Workstream of the Schizophrenia Treatment Resistance and Therapeutic Advances (STRATA-G) Cohorts Using Scores Calculated With Linkage Disequilibrium–Independent (P< 1) Single-Nucleotide Variants Standardized PRS –0.50 –0.25 0 0.25 0.50 Interaction A TRS (prevalence) TRS (incidence) Non-TRS (prevalence) Non-TRS (incidence) Controls Standardized PRS –0.50 –0.25 0 0.25 0.50 CLOZUK B TRS (prevalence) TRS (incidence) Non-TRS (prevalence) Non-TRS (incidence) Controls Standardized PRS –0.50 –0.25 0 0.25 0.50 PGC C TRS (prevalence) TRS (incidence) Non-TRS (prevalence) Non-TRS (incidence) Controls Dots indicate the mean values of the standardized polygenic risk score (PRS; corrected for sex, relevant principal components, and ancestry-informative marker–based ancestry estimates). Whiskers delimit the 95% CIs of the PRS means. PGC indicates Psychiatric Genomics Consortium; TRS, treatment-resistant schizophrenia. Contribution of Common Genetic Variants to Treatment Resistance in Schizophrenia Original Investigation Research jamapsychiatry.com (Reprinted) JAMA Psychiatry March 2022 Volume 79, Number 3 265 Downloaded From: https://jamanetwork.com/ by a Universidad de Sevilla User on 11/22/2022 the genetics of disorder susceptibility are the same as disease course in schizophrenia 2 and whether treatment response can be influenced by the accumulation of genetic and clinical factors. 54,55 Finally, while our results show that TRS is associated with a polygenic signal, the variance in TRS explained remains modest, and associated area under the curve values are small (eTable 3 in Supplement 1). Thus, polygenic scores tapping into this signal are thus unlikely to be of clinical utility in predicting treatment resistance, although their contribution to multifactorial predictive models (integrating rare variants, neuroimaging biomarkers, environmental exposures, demographic factors, and clinical measures) is an interesting avenue for future research. Limitations This study has limitations. TRS is an underreported diagnosis, and while our research definition for this phenotype aligns with international criteria, 56 we acknowledge that some individuals with treatment-resistant symptoms might still be present in the non-TRS data set, particularly in those samples where no ascertainment of TRS could be carried out. The effect of such misclassification would be akin to a reduction of effective sample size, reducing SNV discovery power and adding noise to GWAS effect size estimates. 57 This further limits our power to detect differences in allelic associations through the interaction analysis, which is already only sensitive to relatively large differences in effect size between studies. 58 Additionally, the interaction test does not model or addresses potential between-study heterogeneity. This precludes the use of more sophisticated statistical methods for downstream analyses, such as tissue-specific enrichment analyses or transcriptome-wide association studies. 59 Imperfect phenotyping and potential misclassification might also have attenuated our results, and given the small total heritability detected for our TRS GWAS, we note that the magnitude of genetic correlations with other traits should be interpreted cautiously. Also, similar to what has been argued in the study of etiologically heterogeneous phenotypes through GWAS, 50,60 the existence of significant genetic correlations and consistent polygenic association in independent samples reassures us that the common genetic variants we detect have (en masse) an association with TRS liability, even if our results cannot quantify the degree of contribution from nongenetic causes. In this regard, given that most of our analyses are based on European-based and UK-based samples, our conclusions might not be generalizable to non- European countries, where the diagnosis and treatment pathways of TRS might be influenced by racial and ethnic or cultural backgrounds. 61 Conclusions In this GWAS, TRS had a small but detectable heritability associated with common risk alleles. Validation work, via the PRS method, showed that the contribution of these alleles were similar across incidence and prevalence samples. This is despite differences in cohort characteristics, and the use of drug prescription data as a proxy for TRS instead of quantitative metrics. Altogether, these results highlight the usefulness of well-controlled clinical phenotype data in psychiatric genetics to explore beyond diagnostic classifications and into treatment outcome and response to aid precision psychiatry. ARTICLE INFORMATION Accepted for Publication: November 2, 2021. Published Online: January 12, 2022. doi:10.1001/jamapsychiatry.2021.3799 Open Access: This is an open access article distributed under the terms of the CC-BY License. © 2022 Pardiñas AF et al. JAMA Psychiatry. Author Affiliations: MRC Centre for Neuropsychiatric Genetics and Genomics, Division of Psychological Medicine and Clinical Neurosciences, School of Medicine, Cardiff University, Cardiff, United Kingdom (Pardiñas, Smart, Willcocks, Holmans, Dennison, Lynham, Legge, Owen, O’Donovan, Walters); Department of Psychosis Studies, Institute of Psychiatry, Psychology & Neuroscience, King’s College London, London, United Kingdom (Smart, Murray, MacCabe); Department of Psychiatry, University of Münster, Münster, Germany (Baune); Department of Psychiatry, Melbourne Medical School, The University of Melbourne, Melbourne, Australia (Baune); The Florey Institute of Neuroscience and Mental Health, The University of Melbourne, Melbourne, Australia (Baune); Department of Psychiatry and the Behavioral Sciences, State University of New York Downstate Medical Center, Brooklyn (Bigdeli, Fanous); Institute for Genomic Health, State University of New York Downstate Medical Center, Brooklyn (Bigdeli, Fanous); Department of Psychiatry, Veterans Affairs New York Harbor Healthcare System, Brooklyn (Bigdeli); School of Biomedical Sciences and Pharmacy, University of Newcastle, Newcastle, Australia (Cairns, Tooney); Centre for Brain and Mental Health Research, University of Newcastle, Newcastle, Australia (Cairns); Hunter Medical Research Institute, Newcastle, Australia (Cairns, Tooney); Neuropsychiatric Genetics Research Group, Department of Psychiatry, Trinity College Dublin, Dublin, Ireland (Corvin); Department of Genetic Epidemiology in Psychiatry, Central Institute of Mental Health, Medical Faculty Mannheim, University of Heidelberg, Heidelberg, Mannheim, Germany (Frank, Rietschel); School of Medicine & Public Health, The University of Newcastle, Newcastle, Australia (Kelly); Molecular Psychiatry Laboratory, Division of Psychiatry, University College London, London, United Kingdom (McQuillin); Norwegian Centre for Mental Disorders Research, Institute of Clinical Medicine, University of Oslo, Oslo, Norway (Melle, Simonsen, Andreassen); Division of Mental Health and Addiction, Institute of Clinical Medicine, Oslo University Hospital, Oslo, Norway (Melle, Andreassen); National Centre for Register-based Research, Aarhus University, Aarhus, Denmark (Mortensen); The Lundbeck Foundation Initiative for Integrative Psychiatric Research, Aarhus, Denmark (Mortensen); Queensland Brain Institute, The University of Queensland, Brisbane, Australia (Mowry, Periyasamy); Queensland Centre for Mental Health Research, The University of Queensland, Brisbane, Australia (Mowry, Periyasamy); Department of Psychiatry and Behavioral Sciences, State University of New York Downstate Medical Center, Brooklyn (Pato); Department of Psychiatry and Zilkha Neurogenetics Institute, Keck School of Medicine, University of Southern California, Los Angeles (Pato); Institute for Genomic Health, State University of New York Downstate Medical Center, Brooklyn (Pato); University Clinic and Outpatient Clinic for Psychiatry, Psychotherapy and Psychosomatics, Martin Luther University of Halle-Wittenberg, Halle, Germany (Rujescu); Division of General Psychiatry, Department of Psychiatry and Psychotherapy, Medical University of Vienna, Vienna, Austria (Rujescu); Early Intervention in Psychosis Advisory Unit for South-East Norway, Division of Mental Health and Addiction, Oslo University Hospital, Oslo, Norway (Simonsen); Institute of Medical Sciences, University of Aberdeen, Aberdeen, United Kingdom (St Clair); Baker Heart and Diabetes Institute, Melbourne, Australia (Wu); College of Pharmacy, University of Manitoba, Winnipeg, Manitoba, Canada (Kowalec); Department of Medical Epidemiology and Biostatistics, Karolinska Institutet, Stockholm, Sweden (Kowalec, Sullivan); Department of Psychiatry, Icahn School of Medicine, Mount Sinai Hospital, New York, New York (Sullivan); Research Original Investigation Contribution of Common Genetic Variants to Treatment Resistance in Schizophrenia 266 JAMA Psychiatry March 2022 Volume 79, Number 3 (Reprinted) jamapsychiatry.com Downloaded From: https://jamanetwork.com/ by a Universidad de Sevilla User on 11/22/2022 Department of Genetics, University of North Carolina, Chapel Hill (Sullivan). Authors of the Genetics Workstream of the Schizophrenia Treatment Resistance and Therapeutic Advances (STRATA) Consortium: Olesya Ajnakina, PhD; Luis Alameda, MD; Thomas R. E. Barnes, DSc; Domenico Berardi, MD; Elena Bonora, PhD; Sara Camporesi, MD; Martine Cleusix, PhD; Philippe Conus, MD; Benedicto Crespo-Facorro, MD, PhD; Giuseppe D’Andrea, MD; Arsime Demjaha, MBBS, PhD; Kim Q. Do, PhD; Gillian A. Doody, MD; Chin B. Eap, PhD; Aziz Ferchiou, MD; Marta Di Forti, MD, PhD; Lorenzo Guidi, MSc; Lina Homman, PhD; Raoul Jenni, PhD; Eileen M. Joyce, MD, PhD; Laura Kassoumeri, BSc; Inès Khadimallah, PhD; Ornella Lastrina, PhD; Roberto Muratori, MD, PhD; Handan Noyan, PhD; Francis A. O’Neill, PhD; Baptiste Pignon, MD; Romeo Restellini, MD; Jean-Romain Richard, MSc; Franck Schürhoff, MD; Filip Španiel, MD, PhD; Andrei Szöke, MD, PhD; Ilaria Tarricone, MD, PhD; Andrea Tortelli, MD, PhD; Alp Üçok, MD; Javier Vázquez-Bourgon, MD, PhD. Affiliations of Authors of the Genetics Workstream of the Schizophrenia Treatment Resistance and Therapeutic Advances (STRATA) Consortium: Department of Biostatistics & Health Informatics, Institute of Psychiatry, Psychology & Neuroscience, King’s College London, University of London, London, United Kingdom (Ajnakina); Department of Behavioural Science and Health, Institute of Epidemiology and Health Care, University College London, London, United Kingdom (Ajnakina); Department of Psychosis Studies, Institute of Psychiatry, Psychology & Neuroscience, King’s College London, London, United Kingdom (Alameda, Demjaha, Kassoumeri); Centro de Investigacion Biomedica en Red de Salud Mental, Spanish Network for Research in Mental Health, Sevilla, Spain (Alameda, Crespo-Facorro); Instituto de Biomedicina de Sevilla, Hospital Universitario Virgen del Rocio, Departamento de Psiquiatria, Universidad de Sevilla, Sevilla, Spain (Alameda, Crespo-Facorro); Treatment and Early Intervention in Psychosis Program, Service of General Psychiatry, Department of Psychiatry, Lausanne University Hospital, Lausanne, Switzerland (Alameda, Camporesi, Conus, Restellini); Division of Psychiatry, Imperial College London, London, United Kingdom (Barnes); Department of Biomedical and Neuro-motor Sciences, Psychiatry Unit, Alma Mater Studiorum Università di Bologna, Bologna, Italy (Berardi, D’Andrea, Lastrina); Department of Medical and Surgical Sciences, Bologna Transcultural Psychosomatic Team, Alma Mater Studiorum, University of Bologna, Bologna, Italy (Bonora, Guidi, Muratori, Tarricone); Unit for Research in Schizophrenia, Center for Psychiatric Neuroscience, Department of Psychiatry, Lausanne University Hospital, Lausanne, Switzerland (Camporesi, Cleusix, Do, Jenni, Khadimallah, Restellini); Department of Medical Education, University of Nottingham Faculty of Medicine and Health Sciences, Nottingham, United Kingdom (Doody); Unit of Pharmacogenetics and Clinical Psychopharmacology, Centre for Psychiatric Neuroscience, Department of Psychiatry, Lausanne University Hospital, University of Lausanne, Prilly, Switzerland (Eap); School of Pharmaceutical Sciences, University of Geneva, Geneva, Switzerland (Eap); Center for Research and Innovation in Clinical Pharmaceutical Sciences, Lausanne University Hospital and University of Lausanne, Lausanne, Switzerland (Eap); Institute of Pharmaceutical Sciences of Western Switzerland, University of Geneva, Geneva, Switzerland (Eap); University Paris-Est Créteil, Institut national de la santé et de la recherche médicale, Mondor Institute for Biomedical Research, Translational Neuropsychiatry, Fondation FondaMental, Créteil, France (Ferchiou, Pignon, Richard, Schürhoff, Szöke, Tortelli); Social Genetics and Developmental Psychiatry, Institute of Psychiatry, Psychology & Neuroscience, King’s College London, London, United Kingdom (Di Forti); South London and Maudsley National Health Service Mental Health Foundation Trust, London, United Kingdom (Di Forti); Department of Social and Welfare Studies, Department of Behavioural Sciences and Learning, Linköping University, Linköping, Sweden (Homman); Centre For Public Health, Institute Of Clinical Sciences, Queens University Belfast, Belfast, United Kingdom (Homman, O’Neill); UCL Queen Square Institute of Neurology, University College London, London, United Kingdom (Joyce); Faculty of Social Sciences, Department of Psychology, Beykoz University, Istanbul, Turkey (Noyan); Assistance Publique–Hôpitaux de Paris, Hôpitaux Universitaires HMondor, Département Médico-Universitaire de Psychiatrie et d’Addictologie, Fédération Hospitalo-Universitaire de Médecine de Précision, Créteil, France (Pignon, Schürhoff, Szöke); Department of Applied Neuroscience and Neuroimaging, National Institute of Mental Health, Klecany, Czechia (Španiel); Department of Psychiatry and Medical Psychology, Third Faculty of Medicine, Charles University, Prague, Czechia (Španiel); Groupe Hospitalier Universitaire Psychiatrie Neurosciences Paris, Pôle Psychiatrie Précarité, Paris, France (Tortelli); Department of Psychiatry, Istanbul University, Istanbul, Turkey (Üçok); Department of Psychiatry, University Hospital Marques de Valdecilla-Instituto de Investigación Marques de Valdecilla, Santander, Spain (Vázquez-Bourgon); Department of Medicine and Psychiatry, School of Medicine, University of Cantabria, Santander, Spain (Vázquez-Bourgon); Centro de Investigacion Biomedica en Red de Salud Mental, Spanish Network for Research in Mental Health, Santander, Spain (Vázquez-Bourgon). Author Contributions: Dr Pardiñas had full access to all of the data in the study and takes responsibility for the integrity of the data and the accuracy of the data analysis. Concept and design: Pardiñas, Conus, Corvin, Kelly, McQuillin, O’Neill, Periyasamy, St Clair, Murray, MacCabe, O’Donovan, Walters. Acquisition, analysis, or interpretation of data: Pardiñas, Smart, Willcocks, Holmans, Dennison, Lynham, Legge, Kassoumeri, Ajnakina, Alameda, Barnes, Baune, Berardi, Bigdeli, Bonora, Cairns, Camporesi, Cleusix, Conus, Corvin, Crespo-Facorro, D’Andrea, Demjaha, Do, Doody, Eap, Fanous, Ferchiou, Di Forti, Frank, Guidi, Homman, Jenni, Joyce, Khadimallah, Lastrina, McQuillin, Melle, Mortensen, Mowry, Muratori, Noyan, Pato, Pignon, Restellini, Richard, Rietschel, Rujescu, Schurhoff, Simonsen, Spaniel, St Clair, Szöke, Tarricone, Tooney, Tortelli, Ucok, Vazquez-Bourgon, Wu, Andreassen, Kowalec, Sullivan, Murray, Owen, MacCabe, O’Donovan, Walters. Drafting of the manuscript: Pardiñas, Kassoumeri, Joyce, Khadimallah, Lastrina, O’Neill, Schurhoff, O’Donovan, Walters. Critical revision of the manuscript for important intellectual content: Pardiñas, Smart, Willcocks, Holmans, Dennison, Lynham, Legge, Ajnakina, Alameda, Barnes, Baune, Berardi, Bigdeli, Bonora, Cairns, Camporesi, Cleusix, Conus, Corvin, Crespo-Facorro, D’Andrea, Demjaha, Do, Doody, Eap, Fanous, Ferchiou, Di Forti, Frank, Guidi, Homman, Jenni, Kelly, McQuillin, Melle, Mortensen, Mowry, Muratori, Noyan, Pato, Pignon, Periyasamy, Restellini, Richard, Rietschel, Rujescu, Simonsen, Spaniel, St Clair, Szöke, Tarricone, Tooney, Tortelli, Ucok, Vazquez-Bourgon, Wu, Andreassen, Kowalec, Sullivan, Murray, Owen, MacCabe, O’Donovan, Walters. Statistical analysis: Pardiñas, Smart, Willcocks, Holmans, Lynham, Bigdeli, Frank. Obtained funding: Conus, Corvin, Do, Di Forti, Melle, Mortensen, Mowry, O’Neill, Pato, Rietschel, Spaniel, St Clair, Sullivan, Murray, Owen, MacCabe, O’Donovan, Walters. Administrative, technical, or material support: Pardiñas, Dennison, Legge, Kassoumeri, Ajnakina, Alameda, Berardi, Bonora, Cairns, Doody, Eap, Fanous, Homman, Joyce, Khadimallah, McQuillin, Melle, Noyan, Pato, Periyasamy, Restellini, Richard, Rietschel, Simonsen, Spaniel, Ucok, Wu, Andreassen, Kowalec, Murray, Walters. Supervision: Conus, Crespo-Facorro, Do, Doody, McQuillin, Pignon, Rujescu, Murray, Owen, MacCabe, O’Donovan, Walters. Conflict of Interest Disclosures: Dr Baune has received grants from EraPerMEd as well as personal fees from AstraZeneca, Lundbeck, Pfizer, Takeda Pharmaceuticals, Servier, Bristol Myers Squibb, Otsuka, LivaNova, and Janssen-Cilag. Dr Andreassen has received grants from KG Jebsen Stiftelsen, South East Norway Health Authority, Research Council of Norway, and the European Union Horizon 2020 as well as personal fees from HealthLytix, Lundbeck, and Sunovion. Dr Sullivan has received research funding from Lundbeck as well as personal fees from Pfizer, Neumora, Element Genomics, and Roche; is a shareholder of Neumora; and serves on the advisory committee of Lundbeck. Dr Murray has received grants from Medical Research Council as well as personal fees from Janssen, Merck, Otsuka, Sunovian, and Lundbeck. Drs Owen, O’Donovan, and Walters are investigators on a grant from Takeda Pharmaceuticals paid to Cardiff University; Dr Smart is employed by funds from this grant. Dr MacCabe has received research funding and nonfinancial support from Lundbeck. Dr O’Donovan reported grants from the National Institutes of Health, UK Medical Research Council, Commission of the European Union, and Takeda Pharmaceuticals. Dr Walters has received grants from Medical Research Council, National Institutes of Health, and European Union 7th Framework Programme for Research, and Takeda Pharmaceuticals. Dr Barnes is a member of an advisory board for Gedeon Richter. Dr Camporesi has received grants from the Swiss National Centre Competence in Research. Dr Crespo-Facorro has received personal fees from Adamed, Mylan, Angelini, Janssen, Johnson & Johnson, Lundbeck, and Otsuka Pharmaceuticals. Dr Do has received grants from the Swiss National Science Foundation, and Boehringer Ingelheim and has a patent for EP19218841.5 pending. Dr Eap has received grants from the Swiss National Science Foundation and personal fees from Janssen-Cilag, Lundbeck, Otsuka, Sandoz, Servier, Sunovion, Vifor Pharma, and Zeller. Dr Di Forti has received fees for educational seminars from Lundbeck and Janssen. Dr Szöke has received grants from the European Contribution of Common Genetic Variants to Treatment Resistance in Schizophrenia Original Investigation Research jamapsychiatry.com (Reprinted) JAMA Psychiatry March 2022 Volume 79, Number 3 267 Downloaded From: https://jamanetwork.com/ by a Universidad de Sevilla User on 11/22/2022 Community’s Seventh Framework Programme. No other disclosures were reported. Funding/Support: This work was supported by Medical Research Council Centre grant MR/ L010305/1, Medical Research Council Program grant MR/P005748/1, and Medical Research Council Project grants MR/L011794/1 and MC_PC_17212 to Cardiff University and by the National Centre for Mental Health, funded by the Welsh Government through Health and Care Research Wales. This work acknowledges the support of the Supercomputing Wales project, which is partially funded by the European Regional Development Fund via the Welsh Government. Dr Pardiñas was supported by an Academy of Medical Sciences Springboard Award (SBF005\1083). Dr Andreassen was supported by the Research Council of Norway (grants 283798, 262656, 248980, 273291, 248828, 248778, and 223273); KG Jebsen Stiftelsen, South-East Norway Health Authority, and the European Union’s Horizon 2020 Research and Innovation Programme (grant 847776). Dr Ajnakina was supported by an National Institute for Health Research postdoctoral fellowship (PDF-2018-11-ST2-020). Dr Joyce was supported by the University College London Hospitals/UCL University College London Biomedical Research Centre. Dr Kowalec received funding from the European Union’s Horizon 2020 Research and Innovation Programme under the Marie Skłodowska-Curie grant agreement (793530) from the government of Canada Banting postdoctoral fellowship programme and the University of Manitoba. Dr Sullivan was supported by the Swedish Research Council (Vetenskapsrådet, D0886501), the European Union’s Horizon 2020 programme (COSYN, 610307) and the US National Institute of Mental Health (U01 MH109528 and R01 MH077139). The Psychiatric Genomics Consortium was partly supported by the National Institute Of Mental Health (grants R01MH124873). The Sweden Schizophrenia Study was supported by the National Institute Of Mental Health (grant R01MH077139). The STRATA consortium was supported by a Stratified Medicine Programme grant to Dr MacCabe from the Medical Research Council (grant MR/L011794/1), which funded the research and supported Drs Pardiñas, Smart, Kassoumeri, Murray, Walters, and MacCabe. Dr Smart was supported by a Collaboration for Leadership in Applied Health Research and Care South London at King’s College Hospital National Health Service Foundation Trust. The AESOP (US) cohort was funded by the UK Medical Research Council (grant G0500817). The Belfast (UK) cohort was funded by the Research and Development Office of Northern Ireland. The Bologna (Italy) cohort was funded by the European Community’s Seventh Framework program (HEALTH-F2-2010–241909, project EU-GEI). The Genetics and Psychosis project (London, UK) cohort was funded by the UK National Institute of Health Research Specialist Biomedical Research Centre for Mental Health, South London and the Maudsley National Health Service Mental Health Foundation Trust (SLAM) and the Institute of Psychiatry, Psychology, and Neuroscience at King’s College London; Psychiatry Research Trust; Maudsley Charity Research Fund; and the European Community’s Seventh Framework program (HEALTH-F2-2009-241909, project EU-GEI). The Lausanne (Switzerland) cohort was funded by the Swiss National Science Foundation (grants 320030_135736/1, 320030-120686, 324730-144064, 320030-173211, and 171804); the National Center of Competence in Research Synaptic Bases of Mental Diseases from the Swiss National Science Foundation (grant 51AU40_125759); and Fondation Alamaya. The Oslo (Norway) cohort was funded by the Research Council of Norway (grant 223273/F50, under the Centers of Excellence funding scheme, 300309, 283798) and the South-Eastern Norway Regional Health Authority (grants 2006233, 2006258, 2011085, 2014102, 2015088, and 2017-112). The Paris (France) cohort was funded by European Community’s Seventh Framework program (HEALTH-F2-2010–241909, project EU-GEI). The Prague (Czech Republic) cohort was funded by the Ministry of Health of the Czech Republic (grant NU20-04-00393). The Santander (Spain) cohort was funded by the following grants to Dr Crespo-Facorro: Instituto de Salud Carlos III (grants FIS00/3095, PI020499, PI050427, and PI060507), Plan Nacional de Drogas Research (grant 2005-Orden sco/3246/2004), SENY Fundatio Research (grant 2005-0308007), Fundacion Marques de Valdecilla (grant A/02/07, API07/011) and Ministry of Economy and Competitiveness and the European Fund for Regional Development (grants SAF2016-76046-R and SAF2013-46292-R). The West London (UK) cohort was funded by The Wellcome Trust (grants 042025, 052247, and 064607). Role of the Funder/Sponsor: The funders had no role in the design and conduct of the study; collection, management, analysis, and interpretation of the data; preparation, review, or approval of the manuscript; and decision to submit the manuscript for publication. Group Information: The Genetics Workstream of the Schizophrenia Treatment Resistance and Therapeutic Advances (STRATA) Consortium and the Schizophrenia Working Group of the Psychiatric Genomics Consortium (PGC) collaborators appear in Supplement 2. Disclaimer: The research content reported in this publication is solely the responsibility of the authors and does not necessarily represent the official views of the National Institutes of Health. The views expressed are those of the authors and not necessarily those of the Medical Research Council, National Health Service, the National Institute for Health Research, or the Department of Health. Additional Information: Summary statistics from the treatment-resistant schizophrenia interaction genome-wide association study can be downloaded from https://walters.psycm.cf.ac.uk/. REFERENCES 1. Sullivan PF, Agrawal A, Bulik CM, et al; Psychiatric Genomics Consortium. Psychiatric genomics: an update and an agenda. Am J Psychiatry. 2018;175 (1):15-27. doi:10.1176/appi.ajp.2017.17030283 2. Paternoster L, Tilling K, Davey Smith G. Genetic epidemiology and mendelian randomization for informing disease therapeutics: conceptual and methodological challenges. PLoS Genet.2017;13 (10):e1006944. doi:10.1371/journal.pgen.1006944 3. Lee JC, Biasci D, Roberts R, et al; UK IBD Genetics Consortium. Genome-wide association study identifies distinct genetic contributions to prognosis and susceptibility in Crohn’s disease. Nat Genet. 2017;49(2):262-268. doi:10.1038/ng.3755 4. Rees E, Owen MJ. Translating insights from neuropsychiatric genetics and genomics for precision psychiatry. Genome Med. 2020;12(1):43. doi:10.1186/s13073-020-00734-5 5. Birnbaum R, Weinberger DR. Special article: translational science update. pharmacological implications of emerging schizophrenia genetics: can the bridge from ‘genomics’ to ‘therapeutics’ be defined and traversed? J Clin Psychopharmacol. 2020;40(4):323-329. doi:10.1097/JCP. 0000000000001215 6. Lally J, Gaughran F, Timms P, Curran SR. Treatment-resistant schizophrenia: current insights on the pharmacogenomics of antipsychotics. Pharmgenomics Pers Med. 2016;9:117-129. doi:10.2147/PGPM.S115741 7. Mørup MF, Kymes SM, Oudin Åström D. A modelling approach to estimate the prevalence of treatment-resistant schizophrenia in the United States. PLoS One. 2020;15(6):e0234121. doi:10.1371/journal.pone.0234121 8. Siskind D, Orr S, Sinha S, et al. Rates of treatment-resistant schizophrenia from first-episode cohorts: systematic review and meta-analysis. Br J Psychiatry. 2021:1-6. doi:10.1192/ bjp.2021.61 9. Kumra S, Kranzler H, Gerbino-Rosen G, et al. Clozapine and “high-dose” olanzapine in refractory early-onset schizophrenia: a 12-week randomized and double-blind comparison. Biol Psychiatry. 2008;63(5):524-529. doi:10.1016/j.biopsych.2007. 04.043 10. Conley RR, Kelly DL. Management of treatment resistance in schizophrenia. Biol Psychiatry. 2001; 50(11):898-911. doi:10.1016/S0006-3223(01)01271-9 11. Tiihonen J, Mittendorfer-Rutz E, Majak M, et al. Real-world effectiveness of antipsychotic treatments in a nationwide cohort of 29 823 patients with schizophrenia. JAMA Psychiatry.2017; 74(7):686-693. doi:10.1001/jamapsychiatry.2017.1322 12. Pardiñas AF, Owen MJ, Walters JTR. Pharmacogenomics: a road ahead for precision medicine in psychiatry. Neuron. 2021;S0896-6273 (21)00683-8. 13. Krivoy A, Hochman E, Sendt K-V, et al. Association between serum levels of glutamate and neurotrophic factors and response to clozapine treatment. Schizophr Res. 2018;192:226-231. doi:10.1016/j.schres.2017.05.040 14. Shah P, Iwata Y, Brown EE, et al. Clozapine response trajectories and predictors of non-response in treatment-resistant schizophrenia: a chart review study. Eur Arch Psychiatry Clin Neurosci. 2020;270(1):11-22. doi:10.1007/s00406- 019-01053-6 15. Shah P, Iwata Y, Plitman E, et al. The impact of delay in clozapine initiation on treatment outcomes in patients with treatment-resistant schizophrenia: a systematic review. Psychiatry Res. 2018;268:114- 122. doi:10.1016/j.psychres.2018.06.070 16. Kowalec K, Lu Y, Sariaslan A, et al. Increased schizophrenia family history burden and reduced premorbid IQ in treatment-resistant schizophrenia: a Swedish National Register and Genomic Study. Mol Psychiatry. 2021;26(8):4487-4495. doi:10.1038/s41380-019-0575-1 17. Gillespie AL, Samanaite R, Mill J, Egerton A, MacCabe JH. Is treatment-resistant schizophrenia categorically distinct from treatment-responsive schizophrenia? a systematic review. BMC Psychiatry. 2017;17(1):12. doi:10.1186/s12888-016-1177-y 18. Sobell JL, Mikesell MJ, McMurray CT. Genetics and etiopathophysiology of schizophrenia. Mayo Clin Proc. 2002;77(10):1068-1082. doi:10.4065/77. 10.1068 19. Frank J, Lang M, Witt SH, et al. Identification of increased genetic risk scores for schizophrenia in Research Original Investigation Contribution of Common Genetic Variants to Treatment Resistance in Schizophrenia 268 JAMA Psychiatry March 2022 Volume 79, Number 3 (Reprinted) jamapsychiatry.com Downloaded From: https://jamanetwork.com/ by a Universidad de Sevilla User on 11/22/2022