Association Between Polygenic Risk Scores and Outcome of ECT

The study population was derived from the Predictors for ECT (PREFECT) study, conducted in Sweden (12). Recruitment occurred between 2013 and 2017. Study participants had received at least one acute ECT series that was registered in the Swedish National Quality Register for ECT (Q-ECT) (http://ect.registercentrum.se). Q-ECT was launched nationally in Sweden in 2011, although some hospitals started registration in 2008. As of 2014, Q-ECT covered 89% of all ECT series in Sweden (20). Local hospital staff routinely enter data on each treatment series into the register via a web-based platform.

The study sample has been described previously (12, 21). Briefly, a letter of invitation was sent to patients ≥18 years old registered in Q-ECT. Those who volunteered to participate completed a telephone interview and donated blood that was sent via overnight mail to Karolinska Institutet Biobank. Additional participants were prospectively recruited at eight Swedish hospitals prior to receiving ECT for a major depressive episode. Blood samples were stored at −20°C pending shipment to Karolinska Institutet Biobank for DNA extraction and long-term storage.

A total of 2,880 participants provided DNA and had received at least one acute ECT treatment series for a major depressive episode. After exclusion of participants whose genotyping failed quality control (N=40), genetic ancestral outliers (N=138), and those with no registered primary outcome in Q-ECT (N=382), we categorized the remaining 2,320 eligible participants into two groups, as previously described (12): The broadly defined group included participants with a major depressive episode occurring in the context of any one of the following: 1) a unipolar depressive disorder (ICD-10 codes F32–F33, F41.2, F53.0), 2) bipolar disorder (code F31), 3) another severe mood disorder with a pretreatment score ≥20 on the self-rated Montgomery-Åsberg Depression Rating Scale (MADRS-S) (mixed anxiety and depressive disorder, code F41.2; schizoaffective disorder, code F25.9; or mood disorder not otherwise specified, code F34.9), or 4) major depressive episode as indicated by free text, or by a pretreatment MADRS-S score ≥20 if the specific indication for treatment was missing. The narrowly defined group consisted of the subset of participants who received ECT for a major depressive episode in the context of a unipolar depressive episode only. Psychotic features were considered present if the indication for ECT was coded with any of ICD-10 codes F323, F333, F315, and F259, or if the indication in free text implied delusions and/or hallucinations.

All participants provided written informed consent. The study was approved by the regional ethical review board in Stockholm (approval nos. 2012/1969-31/1 and 2020-10151).

Participants could have several registered treatment series in Q-ECT (Table 1). Here, outcome data were captured from the first or only ECT treatment series. Our primary outcome measure was the Clinical Global Impressions improvement (CGI-I) score (22), rated immediately after the ECT series. We chose the CGI-I because it had the lowest proportion of missing information. CGI-I scores were available for the first registered treatment series of 92.8% of participants (N=2,154). CGI-I score, which ranges from 1 (very much improved) to 7 (very much worse), was treated as an ordinal variable. Secondary outcome measures were response and remission according to the MADRS-S (23). Both pre- and posttreatment MADRS-S scores were available for 52.0% of the sample (N=1,207), of which 79.2% (N=956) provided MADRS-S ratings at their first registered treatment series. The MADRS-S includes nine items rated from 0 to 6 (maximum score=54), which correspond to the items in the observer-rated MADRS except “apparent sadness.” Response was defined as a reduction ≥50% in MADRS-S score from pre- to posttreatment assessment. Remission was defined as a posttreatment MADRS-S score ≤10.

A subset of 1,052 participants had both MADRS-S and CGI-I data available from the same treatment series. There was a moderate correlation between CGI-I and MADRS-S response (Spearman’s ρ=0.48, p<0.001) and MADRS-S remission (Spearman’s ρ=0.42, p<0.001).

ECT was administered using bidirectional, constant-current, brief-pulse devices from Mecta (Mecta Corp., Lake Oswego, Ore.) or Thymatron (Somatics, Inc., Lake Buff, Ill.). Propofol or thiopental was used for anesthesia. Suxamethonium was used for muscle relaxation. Seizure time was registered with electroencephalography. From Q-ECT, we retrieved data on electrode placement (bilateral or right unilateral at first or last ECT), and charge (mC) and pulse width (categorized into 0.25–0.49 ms, 0.5 ms, and 0.51–1.20 ms) used at the first session of each treatment series. All participating clinics administered ECT three times per week (Monday, Wednesday, and Friday).

Genotype quality control and imputation have been described in detail elsewhere (12). Briefly, DNA was extracted from peripheral blood, and samples were genotyped on the Illumina GSA-MD SNP array (version GSAMD-24v1-0_20011747_A1) at Life & Brain GmbH (Bonn, Germany). Standard quality control was applied using the PGC RICOPILI pipeline (24). All potential samples were included in quality control, but the final association analysis was performed on the phenotyped subsets outlined above that passed quality control. Samples were excluded (N=40) for genotype missingness >0.02 (after first filtering SNPs with call rate <0.95), genotypic sex ambiguous or not matching phenotypic data, or autosomal heterozygosity |F|>0.2. SNPs were excluded for call rate <0.99, difference in missingness between cases and controls >0.005, minor allele frequency <0.01, or deviating from Hardy-Weinberg equilibrium in cases or controls (p<10^–6). Ancestry outliers were identified by projecting study samples on principal components with respect to 1000 Genomes Project data (phase 3 v5) (25). Individuals further than three standard deviations from the European reference population mean for principal components 1 or 2 were excluded (N=138), and then the principal components were regenerated for use in the study analyses as covariates to capture residual confounding by genetic ancestry. Potentially related individuals were identified, and one from each pair was flagged for exclusion (estimated identity-by-descent sharing >0.2). Samples that passed quality control were imputed using the Haplotype Reference Consortium (HRC) r1.1 reference panel on the Sanger Imputation Service using Eagle2 and positional Burrows-Wheeler transform (PBWF) for phasing and imputation (26). The genome build was hg19.

PRSs were calculated for major depression and bipolar disorder using discovery GWAS summary statistics from the Psychiatric Genomics Consortium (https://www.med.unc.edu/pgc/download-results/) based on large GWASs for each phenotype (major depressive disorder: “mdd2019edinborough” [27]; bipolar disorder: “bip2019” [28]; schizophrenia: “scz2022” [29]). Although the PREFECT sample has not been included in any previous GWAS, we nevertheless removed all Swedish samples from the GWAS summary statistics to reduce the possibility of spurious associations.

PRSs were generated for the PREFECT samples as the sum of the risk allele scores, weighted by their effect size in the discovery samples. We performed linkage disequilibrium clumping (R² <0.1 in 1-Mb windows) on any overlapping SNPs with the 1,000 Genomes Project European samples for the reference (phase 3 v5) (25). PRSs were calculated using PLINK, version 1.9 (30). PRSs were coded as risk increasing and standardized to a mean of 0 and a standard deviation of 1 for interpretability. To reduce the number of comparisons, we used the PRSs calculated at a p threshold of ≤0.05 as exposures in our main analyses, in alignment with previous research (17). Table S1 in the online supplement contains the details on the number of SNPs used in the calculation of each PRS.

We present the descriptive characteristics of the sample using frequency (and percent) or median (and interquartile range). Our primary analyses investigated the associations of PRSs for major depression, bipolar disorder, and schizophrenia with CGI-I score. We used a multivariable proportional odds ordinal logistic regression model adjusting for the first five genetic ancestry principal components. To facilitate interpretation, we reversed the CGI-I values such that odds ratios >1 represent improvement. We also examined the association between quintiles of each PRS and CGI-I, using the lowest quintile as the reference, in multivariable ordinal logistic regressions, given that higher PRS values may carry greater risks. We conducted approximate likelihood-ratio tests to examine violation of the proportional odds assumption. For the secondary analyses of the association between each PRS and MADRS-S remission and response (defined above), we used multivariable logistic regression analyses adjusted for MADRS-S prior to ECT, and the first five genetic ancestry principal components. We present all results as odds ratios with 95% confidence intervals. For all analyses, we calculated the proportion of variance explained by each PRS as the difference between the Nagelkerke pseudo R² of the full multivariable model and that of a model excluding the PRS.

We performed four sensitivity analyses of the primary outcome (CGI-I). First, we repeated the analyses for the subset of participants with a narrowly defined major depressive episode. Second, we repeated the analyses using PRS calculated on the basis of alternative p-value inclusion thresholds (≤5E−8, 1E−5, 1E−3, 0.01, 0.1, 0.5, 1.0). Finally, we repeated the analysis separately for participants with unilateral electrode placement only (vs. bilateral placement) at first or last ECT, and for participants with versus without psychotic features.

We applied a two-tailed Bonferroni-adjusted significance level (p=0.05/3=0.017) in our primary analyses involving CGI-I scores because of the analyses of three PRSs. In our secondary and sensitivity analyses, which were exploratory, we applied the uncorrected statistical significance level (p<0.05). We used SPSS, version 26 (IBM Corp., Armonk, N.Y.) for data management, and STATA, version 16 (StataCorp, College Station, Tex.) for statistical analyses. Figures were produced in R, version 4.0.5 (R Foundation for Statistical Computing, Vienna) using ggplot2, version 3.3.3 (31).

We included 2,320 individuals of European genetic ancestry in the primary outcome analysis (CGI-I) (Table 1). The median age was 51 years, 62.8% were women, and 77.1% had a major depressive episode in the context of major depression (i.e., belonged to the narrowly defined group). The median CGI-I rating after ECT was 2 (corresponding to “much improved”). The distribution of CGI-I ratings is presented in Figure S1 in the online supplement.

Inheritance of a greater burden of common, risk-increasing genetic variants associated with major depression (major depression PRS) was significantly and negatively associated with improvement after ECT (CGI-I, odds ratio for improvement, 0.89 per SD, 95% CI=0.82, 0.96, p=0.002, Nagelkerke R²=0.004) (Figure 1A; see also Table S2 in the online supplement). Participants in the highest quintile of major depression PRS had 31% lower odds of improvement compared to those with the lowest burden (quintile 5 vs. 1, odds ratio=0.69, 95% CI=0.54, 0.87, p=0.002) (Figure 2A; see also Table S3 in the online supplement).

Higher bipolar disorder PRS was significantly and positively associated with improvement after ECT (odds ratio per SD, 1.14, 95% CI=1.05, 1.23, p=0.003, Nagelkerke R²=0.005). Participants in the highest quintile for genetic burden of bipolar disorder had 44% higher odds of improvement after ECT than those with the lowest burden (quintile 5 vs. 1, odds ratio=1.44, 95% CI=1.13, 1.84, p=0.003) (Figure 2B; see also Table S3 in the online supplement). Schizophrenia PRS was not associated with improvement after ECT (odds ratio per SD, 1.04, 95% CI=0.97, 1.14, p=0.247). All models using PRS as a continuous variable met the proportional odds assumption, but not those of major depression PRS and bipolar disorder quintiles, which therefore should be interpreted with caution (see Table S3 in the online supplement). In a post hoc analysis, we also found that a PRS (at p≤0.05) of percentage improvement on antidepressants (32) was not associated with improvement following ECT (see Table S4 in the online supplement).

For the MADRS-S analysis, 1,207 participants were eligible. Participant characteristics were similar to those of the CGI-I sample: 64% were female and 67.5% belonged to the narrowly defined group (a major depressive episode in the context of major depression) (Table 1). After ECT, 60.1% (N=725) met the criterion for response (Table 1). Similar to the primary analysis, higher major depression PRS was associated with lower odds of response (odds ratio per SD, 0.85, 95% CI=0.76, 0.96, p=0.008, Nagelkerke R²=0.008) (Figure 1B; see also Table S2 in the online supplement), while bipolar disorder PRS was associated with higher odds of response (odds ratio per SD, 1.13, 95% CI=1.00, 1.27, p=0.044, Nagelkerke R²=0.004). Schizophrenia PRS was not associated with response (odds ratio per SD, 1.05, 95% CI=0.93, 1.19, p=0.401).

In the MADRS-S sample, 40.1% (N=484) met the criterion for remission after ECT. Higher major depression PRS was associated with lower odds of remission (odds ratio per SD, 0.83, 95% CI=0.73, 0.94, Nagelkerke R²=0.01, p=0.002) (Figure 1C; see also Table S2 in the online supplement), while bipolar disorder PRS was associated with higher odds of remission (odds ratio per SD, 1.15, 95% CI=1.02, 1.29, p=0.023, Nagelkerke R²=0.006), as was schizophrenia PRS (odds ratio per SD, 1.16, 95% CI=1.02, 1.31, p=0.020, Nagelkerke R²=0.006).

We repeated the primary analyses restricting the sample to participants with narrowly defined major depressive episode in the context of major depression and obtained similar results (Figures 1 and 2; see also Tables S2 and S3 in the online supplement). The results were not sensitive to the choice of p-value threshold (see Figure S1 in the online supplement). Results were similar to those of the primary analysis in participants with unilateral electrode placement, but in the limited subsample treated with bilateral electrode placement (N=306), only major depression PRS was associated with improvement (see Table S5 in the online supplement). Among participants without psychotic features, results were similar to those of the main analysis (see Table S6 in the online supplement). No significant effect of any PRS was observed among participants with psychotic features (N=348), although point estimates were similar to those of the primary analysis.