ABSTRACT
Background Major depression (MD) is a heterogeneous disorder; however, the extent to which genetic factors distinguish MD patient subgroups (genetic heterogeneity) remains uncertain. This study sought evidence for genetic heterogeneity in MD.
Methods Using UK Biobank cohort, the authors defined 16 MD subtypes within eight comparison groups (vegetative symptoms, symptom severity, comorbid anxiety disorder, age at onset, recurrence, suicidality, impairment and postpartum depression; N∼3,000-47,000). To compare genetic architecture of these subtypes, subtype-specific genome-wide association studies were performed to estimate SNP-heritability, and genetic correlations within subtype comparison and with other related disorders or traits.
Results MD subtypes were divergent in their SNP-heritability, and genetic correlations both within subtype comparisons and with other related disorders/traits. Three subtype comparisons (age at onset, suicidality, and impairment) showed significant differences in SNP-heritability; while genetic correlations within subtypes comparisons ranged from 0.55 to 0.86, suggesting genetic profiles are only partially shared among MD subtypes. Furthermore, subtypes that are more clinically challenging, e.g., early-onset, recurrent, suicidal, more severely impaired, had stronger genetic correlations with other psychiatric disorders. MD with atypical features showed a positive genetic correlation (+0.40) with BMI while a negative correlation (−0.09) was found in those with non-atypical symptoms. Novel genomic loci with subtype-specific effects were identified.
Conclusions These results provide the most comprehensive evidence to date for genetic heterogeneity within MD, and suggest that the phenotypic complexity of MD can be effectively reduced by studying the subtypes which share partially distinct etiologies.
INTRODUCTION
Major depression (MD) is a common psychiatric disorder that affects 15% of the population during lifetime.(1) Individuals with MD vary considerably in symptoms, severity, course, treatment response, and neurobiology.(2) MD heterogeneity is a major research and clinical challenge.(3) Despite major efforts in epidemiological, clinical, and biological psychiatry, this decade-long challenge remains largely unresolved.(4-6) MD subtypes have been proposed within five major categories that focused on: symptoms (typical versus atypical, with or without concomitant anxiety, etc.), etiology (with or without trauma or postpartum exposure), time of onset/time course (early-versus late-onset, recurrent), sex, and treatment outcome (treatment responsive versus resistant).(6) Many of these subtypes, however, exhibit unclear distinctions in underlying biology, psychosocial factors, and treatment efficacy.(6) One of the key biological component is genetics—the extent to which genetic factors distinguish these MD subtypes (i.e. genetic heterogeneity) is largely unknown.
Given its relatively low heritability (30-40%)(7, 8), identifying MD subtypes that are more heritable is of particular importance. Among the proposed subtypes, the sex difference in heritability is the most intensively studied, and current findings support that MD is more heritable in women than in men.(9) Early-onset, recurrent MD, and postpartum depression have been suggested to confer higher genetic liability from family-based studies, which was subsequently confirmed using polygenic risk scores (PRS) in recent MD genome-wide association studies (GWAS).(9-13) Comparisons of MD subtypes between early-versus late-onset, atypical versus non-atypical, with or without adversity have yielded interesting findings (e.g., the genetic overlap with metabolic traits was only found in MD with atypical features subtype, but not among those with non-atypical symptoms).(14) The studies to-date that have used genetic approaches to index the heterogeneity of MD subtypes are encouraging (summarized in Table 1) but overall impeded by a paucity of large cohorts with similar ascertainment, phenotyping, and genotyping.(5) As a result, a systematic comparison across the MD subtypes is lacking and overall evidence for genetic heterogeneity within MD is inconclusive.
The goal of this study was to test genetic heterogeneity in clinically-informed MD subtypes. To accomplish this, we systematically evaluated 16 subtypes in the unique UK Biobank (UKB) cohort with large-scale genomic data and a wide array of phenotypic measures uniformly assessed. In particular, we compared genetic architectures among subtypes by quantifying differences in heritability (i.e., measuring the relative importance of genetic effects on phenotypic variance) and estimating genetic correlations (i.e., to determine if underlying genetic risk factors are identical) within subtype comparisons and with other traits.
METHODS and MATERIALS
To identify MD subtypes and compare their genetic architectures, we carefully selected phenotypes and large-scale genotype data from the UKB. The full protocol and scripts are available via Github.
Participant and phenotype definitions
UKB is a population-based cohort of over 500,000 adults (age 37-73) from across the United Kingdom.(15) UKB has phenotypic data from questionnaires, health records, biological sampling, and physical measurements. Information about mental health including MD was collected using various sources, including touchscreen questionnaires, nurse interviews, hospital admission records, and web-based mental health questionnaires (MHQ) follow-up. The UKB data profile were available elsewhere (15) and briefly described in Supplementary S1.1.
MD case definition
Cases were identified using five MD definitions, including (i) lifetime MD based on the Composite International Diagnostic Interview (CIDI) Short Form; (ii) ICD-coded MD based on linked hospital admission records; (iii) Probable MD based on Smith et al.(16); (iv) Self-reported MD as part of past and current medical conditions; and (v) MD cardinal symptoms of anhedonia and dysphoria (Supplementary table S2.1). These MD definitions have been used in previous studies.(17-19) Because some definitions were available only for parts of the UKB samples, to maximize sample size for MD subtypes, we included individuals who met criteria for at least one of the five MD definitions as cases. MD subtypes were all nested in the broad MD group but coming from different MD definitions (Supplementary table S2.2).
MD subtypes
According to major clinical features in MD, we defined 16 MD subtypes within eight comparison dimensions including (i) MD with versus without atypical features based on vegetative symptoms of hypersomnia and weight gain; (ii) severe versus mild/moderate MD based on symptom severity defined in Smith et al.(16) or ICD codes; (iii) MD with or without comorbid anxiety disorder either self-reported or based on ICD codes; (iv) early-versus late-onset MD based on age at which first experienced a ≥2-week episode of cardinal symptoms; (v) recurrent MD vs single-episode MD based on the number of episodes self-reported or ICD codes; (vi) MD with or without suicidal thoughts or self-harm either experienced recently or during the worst episode; (vii) MD with mild, moderate, severe impairment on normal roles; and (viii) postpartum depression (PPD), either self-reported or based on ICD codes (Table 2). The majority of these subtypes are included in the five major categories proposed in the previous meta-review; while the subtypes on suicidality and on impairment—related to general outcomes of MD—are extensions of the category focused on treatment outcomes (Supplementary S1.1, table S2.3).(6)
Control group
We used a common control group without lifetime history of MD to compare with all but the subtypes of comorbid anxiety disorder and PPD. From the entire UKB population, we excluded those with any indications of MD using five MD case criteria described above, and two additional exclusion criteria, help-seeking MD and antidepressant use (medication list in Supplementary table S2.4). We further excluded those with ICD-diagnoses of anxiety disorders from the controls for the MD subtype with or without comorbid anxiety disorder. For PPD, we restricted controls to women who reported giving at least one live birth. (Supplementary table S2.1)
Exclusionary criteria for cases and controls
We excluded any case or control who met lifetime criteria for schizophrenia, schizoaffective disorder, and bipolar disorder I (including unipolar mania) (Supplementary table S2.1). Thus, anyone who had ICD-diagnosis of schizophrenia/psychosis, bipolar disorder, mania or reported any use of antipsychotics or lithium for psychiatric symptoms (Supplementary table S2.4) were excluded from analyses. Application of these criteria removed 2,385 MD cases and 231 controls (Supplementary figure S3.1).
Genotyping, quality control, imputation
Genotype data were available for 488,363 UKB participants, after a stringent quality control procedure and imputation using combined reference panels of Haplotype Reference Consortium (HRC) and UK10K merged with 1000 Genomes phase 3.(15) 459,590 individuals remained after the exclusion of subjects with low-quality genotype data, unmatched ID with phenotype data, consent withdrawal, and non-European ancestry outliers (Supplementary figure S3.1).
Statistical analysis
Genome-wide association studies (GWAS)
We generated GWAS summary statistics for MD subtypes to estimate SNP-heritability and genetic correlations for computational efficiency. In the UKB, about 30% of the participants were found to be related to at least one other person in the cohort up to the 3rd degree.(15) Cryptic relatedness within sample could bias results in GWAS, while restricting to the unrelated individuals would cause a major loss of statistical power. We therefore performed the mixed linear model-based GWAS analysis (fastGWA) to retain related individuals in the UKB.(20) We first constructed a sparse genetic relationship matrix (GRM) for all included individuals of European ancestries, and then conducted case-control GWAS for each subtype using fastGWA module in GCTA, adjusting for sex, age, and the first 10 PCs(20) (Supplementary S2.2).
For subtype-specific GWAS with genome-wide significant SNPs (p≤5×10−8), we identified independent genomic loci using SNP2GENE module in FUMA(21); then compared our loci with the latest published MD GWAS results which consisted of samples from the Psychiatric Genomics Consortium (PGC), UKB, and 23andMe.(19)
SNP-Heritability
We estimated SNP-heritability (h2SNP) on the observed scale for each MD subtype using linkage disequilibrium score regression (LDSC).(22) LDSC estimates h2SNP by regressing GWAS summary statistics on LD scores estimated from a reference population (1000 genomes European samples). To allow comparison between subtypes, we converted the observed h2SNP to the liability scale, and as previously suggested(23), we corrected for oversampling and extreme phenotyping using sample prevalence, and proportions of the population as cases and controls. Except for PPD, the subtype-specific population prevalence was calculated as MD lifetime prevalence (15%) scaled by the literature-based proportion of subtype in MD, and we used 85% as the proportion of the population as non-MD controls (details in Supplementary table S2.5). We provided a figure showing the impact of population case prevalence estimates on h2SNP. When comparing heritability estimates within subtype comparisons, we cannot directly test for the statistical significance of the difference in estimates due sample overlap; therefore, we considered that estimates are significantly different when non-overlapping confidence intervals are presented.
Genetic correlation
Genetic correlations (rg) were estimated using High-Definition Likelihood (HDL) method which yields more precise estimates of genetic correlations than LDSC (Supplementary S2.2).(24) We estimated rg within subtype comparisons using the LD reference computed from 335,265 Genomic British individuals in the UKB.
To examine whether the subtypes differ in their genetic overlap with other psychiatric disorders and traits, we also estimated genetic correlations between these MD subtypes and 11 traits (six psychiatric disorders, neuroticism, self-reported well-being, body mass index, and two cognitive traits) and compared results within subtype comparisons. These disorders and traits were chosen given the strong evidence for their genetic correlations with MD, or in some cases, for their causal effects on MD.(13, 19) We have used the summary statistics from the latest published GWAS for the calculations of rg using HDL.(19, 25-34)
Sensitivity analyses
To examine whether our broad MD definition that included less strictly defined cases may bias results, we further restricted the analyses to the CIDI-based definition—previously suggested as the closest to the gold standard for diagnosing MD in the UKB(17, 35, 36)—and performed similar analyses for all subtypes except impairment (Supplementary S2.2).
RESULTS
Of 459,590 individuals included in this study (54% females, mean age at recruitment 57 (SD 8.00)), 126,506 (27.53%) met at least one of the five definitions for MD (i.e., broad MD phenotype). After applying exclusion criteria, we retained 124,121 cases and 250,229 controls. Compared with controls, MD cases had more females (64% vs 47%), higher Townsend deprivation index (mean - 1.33, SD 3.02 vs -1.63, SD 2.90), more lifetime smokers (57% vs 52%), but did not differ in mean BMI (mean 27.3, SD 4.6 vs 27.3, SD 5.0).
The estimates of SNP-heritability varied across the five MD case definitions, and for the broad MD phenotype it was 7.38% (95% CI= 6.75-8.01%) (Supplementary figure S3.2).
Differences in genetic architecture reflect subtype heterogeneity
Overall, SNP-heritability estimates tended to be higher in MD subtypes with more severe manifestation (e.g., MD with atypical features, comorbid anxiety disorder, PPD, severe impairment and severe symptoms subtypes) (Figure 1a). All of the subtype comparisons had higher h2SNP estimates for the more severe manifestation, and three (age at onset, suicidality, and impairment) showed significant differences in h2SNP estimates (Figure 1a-b). All examined genetic correlations within comparisons were significantly less than one and the estimates ranged between 0.55-0.86 (Figure 1c).
The h2SNP estimate for MD with atypical features was the highest among all subtypes, and it was almost twice higher than the estimate for non-atypical MD even though the 95% confidence intervals overlapped (13.35%, CI=7.57-19.13% and 7.48%, CI=6.60-8.36%). The genetic correlation between MD subtype with atypical features and subtype without atypical features was the lowest among all comparisons (rg=0.55, CI=0.44-0.66) (Figure 1c). The two subtypes did not significantly differ in their genetic correlations with PGC MD (Figure 2); instead major differences were found in their correlations with anorexia nervosa and ADHD. Consistent with previous findings (14, 37, 38), MD with atypical features showed a strong positive rg with BMI (0.40, CI=0.34-0.46) while non-atypical MD showed a small negative rg instead (rg=-0.09, CI=-0.13 to -0.06). Furthermore, positive correlations with cognitive traits were observed in non-atypical MD (rg=0.36 and 0.35 with educational attainment and intelligence) which were not found in MD with atypical features (corresponding rg= 0.04 and 0.07).
The MD subtype with severe symptoms had slightly higher h2SNP estimate than the one with mild/moderate symptoms, although the two estimates were not significantly different. The rg within comparison was significantly lower than 1 (0.80, CI=0.68-0.92). However, the two subtypes did not differ in their correlations with other traits except for a stronger rg with schizophrenia found in the subtype with severe symptoms (Figure 1-2).
Assuming the proportions of MD cases with and without comorbid anxiety disorder at 55% and 45%, respectively(39), the former subtype was more heritable than the latter (h2SNP=11.35%, CI=10.12-12.58%, for MD with comorbid anxiety disorder, compared with 9.43%, CI=7.92-10.95%, for MD without anxiety disorder). The rg within this comparison was 0.80 (CI=0.73-0.88) (Figure 1). Furthermore, the subtype with comorbid anxiety disorder showed higher genetic correlations with MD, schizophrenia and neuroticism, as well as lower correlations with cognitive traits, when compared with the subtype without anxiety disorder (Figure 2).
The SNP-heritability of early-onset MD was three times higher than that of the late-onset subtype (9.97%, CI=8.89-11.04% compared with 3.25%, CI=2.46-4.04%). The rg within comparison was 0.76 (CI=0.68-0.85). (Figure 1). Significant differences in their rg with other traits were observed, including higher genetic correlations in early-onset MD with PGC MD, schizophrenia, anorexia nervosa, and autism spectrum disorder, than in late-onset MD (Figure 2).
Recurrent and single-episode MD had similar h2SNP estimates (7.94% and 7.47%). However, their rg was significantly lower than one (0.83, CI=0.74-0.92) (Figure 1). Compared with single-episode cases, recurrent MD had stronger positive correlations with schizophrenia, bipolar disorder, anorexia nervosa, while lower genetic correlation with BMI (Figure 2).
The MD subtype with suicidal thoughts was more heritable than the subtype without (8.14%, CI=7.20-9.07% and 6.25%, CI=5.46-7.03%). The rg within this comparison was 0.79 (CI=0.73-0.84). The two subtypes in this comparison significantly differed in their genetic correlations with the majority of the other traits considered. Compared with the subtype without suicidal thoughts, the suicidal subtype showed substantially higher positive rg with PGC MD, schizophrenia, neuroticism, and negative rg with well-being; while its rg with cognitive traits was much weaker (Figure 1-2).
For subtypes based on impairment, the h2SNP estimates increased with the degree of impairment, roughly in a dose-response relationship, i.e., mild impairment had the lowest h2SNP (3.72%, CI=3.09-4.34%), followed by moderate (5.62%, CI=4.81-6.44%) and severe impairment (10.42%, CI=9.08-11.76%). This dose-response relationship was also reflected in the pair-wise genetic correlation estimates, with the rg comparing mild and severe impairment (0.65, CI=0.58-0.72) markedly lower than the other two correlations (Figure 1). We observed a clear trend, that is, the more severe impairment in the subtype, the stronger genetic correlation it had with other psychiatric disorders and neuroticism (positive rg), and with self-reported well-being (negative rg), while less severe impairment was more strongly associated with cognitive traits (positive rg) and with BMI (negative rg) (Figure 2).
The h2SNP of PPD was estimated at 11.31% (CI=6.61-16.0%) which was higher compared with h2SNP of broad MD phenotype. PPD showed significant positive rg with other psychiatric disorders, with the strongest rg observed in PGC MD as expected (0.61, CI=0.53-0.69), and with neuroticism (rg=0.34) and cognitive traits (rg=0.35 and 0.41 with educational attainment and intelligence), and a negative rg with well-being (rg=-0.39) (Figure 2).
The broad MD definition was used above to allow sufficient statistical power in analyzing each subtype. We further assessed the impact of MD definition by performing a sensitivity analysis based on more strictly defined MD cases. The SNP heritability of the CIDI-based definition was in line with previously reported (h2SNP=15.7%, CI=13.4-18.1%, Supplementary figure S3.2)(17, 35). Restricting the analyses to the CIDI-based cases, the results were highly similar, except for the comparisons of symptom severity and recurrence, where the CIs of the rg estimates now included one due to markedly reduced sample sizes in these subtypes (Supplementary table S2.7).
Stratified GWAS reveal novel subtype-specific loci
Over all 16 subtype-specific GWAS, we identified 47 genome-wide significant loci (45 non-overlapping) associated with nine subtypes. Less than half (22 loci) were significant in our largest GWAS of broad MD. Comparing with the latest published MD GWAS (19), we found 14 loci that have not been reported on MD, with 3 for early-onset, 3 for recurrent, 3 for suicidal MD, 2 for non-suicidal, 1 for non-atypical symptoms, 1 for moderate impairment and 1 for PPD (Table 3; full results on the 45 loci in Supplementary table S2.6). The majority (64%) of these 14 loci showed no statistically significant association with the other subtype in comparison (P>0.05; Supplementary S2.6), suggesting subtype-specific effects. The chromosome 2 locus for recurrent MD, with the leading SNP rs6431690, was significant even after the stringent Bonferroni correction (P<3.125*10−9).
DISCUSSION
In this comprehensive report using the large-scale UKB data, we compared the genetic architectures of 16 MD subtypes and demonstrated that these subtypes were divergent in their SNP-heritability and genetic correlations both within subtype comparisons and with other related disorders/traits. Our results provide convincing evidence for genetic heterogeneity within MD, as indexed by its clinical subtypes. These findings suggest that the complexity in the phenotype of MD can be effectively reduced by studying the subtypes which share partially distinct etiologies. In particular, we note the following key findings:
First, clinically-informed subtypes are, in general, genetically more homogeneous than considering all types of MD together. Accurately identifying more homogenous forms is the first step to reduce heterogeneity in MD. The majority of the subtypes showed higher estimates of SNP-heritability compared with MD of all forms. Our results corroborated previous findings from family-based studies that early-onset, recurrent MD and PPD represent more heritable MD subtypes.(10, 12) We further extended the list to include MD with atypical features, MD with or without comorbid anxiety disorder, and with severe impairment. In contrast, subtypes with lower heritability than all-form MD are those with mild/moderate clinical manifestation or with late onset.
Second, we demonstrated subtype heterogeneity in both SNP heritability and genetic correlations. All subtype comparisons showed non-identical genetic sharing (i.e., rg between subtypes significantly differ from unity) and some had heritability differences (i.e., h2SNP significantly differ between subtypes). Interestingly, the subtype comparisons on vegetative symptoms, time onset, and impairment showed the strongest evidence for genetic heterogeneity, meaning these clinical features characterize major etiological differences within MD.
However, the observed genetic correlations across subtype comparisons were moderate to high, 0.55-0.86, revealing substantial genetic overlaps between subtypes. The level of genetic correlation can be translated into the proportion of genetic variance in one trait attributable to that of another (rg2).(17) Thus, it would suggest about 30-70% of genetic variances are shared within subtype comparisons. In line with previous estimates of genetic correlations between male versus female MD(9) and across MD symptoms(40), our findings confirm that the genetic profiles of MD subtypes are only partially distinct.
The estimates of genetic correlations between subtypes need to be benchmarked against genetic correlations between different psychiatric disorders (e.g., schizophrenia and bipolar disorder, two clinically distinct psychiatric disorders, had a rg of ∼0.70 (28)), between different datasets but with same phenotype (e.g., mean rg∼0.76 across the seven cohorts at PGC MD(13)), and between different populations (e.g., rg∼1 between schizophrenia samples of East Asian and European ancestries(41)). Genetic correlations can be found lower than one due to differences in phenotype definitions, populations, or technical factors(42). In this study, we minimized these potential differences by using the single large sample from UKB. We also restricted the estimation of genetic correlations to within subtype comparisons instead of pair-wise comparisons across all subtypes, to limit the impact of phenotypic differences between subtypes (e.g., we found mean rg across all subtypes was indeed lower than that within comparison groups). Our genetic correlation estimates are thus reliable for quantifying overall genetic sharing between MD subtypes.
Third, the MD subtypes preserve the overall pattern of genetic sharing found between MD (of all forms) and other psychiatric disorders, but differ in the relationships with other traits. MD was shown to be positively correlated with many psychiatric disorders (e.g., rg∼0.3 with schizophrenia and bipolar disorders) and with BMI (rg=0.09), and negatively correlated with educational attainment (rg=-0.13).(13, 19) A similar level of genetic correlations was found between MD subtypes and other psychiatric disorders; notably, we found stronger correlations in the MD subtypes that are more clinically challenging, especially early-onset, recurrent, suicidal, more severely impaired. Regarding their relationships with other traits, MD subtypes showed some differences compared with all MD. The positive correlation found between MD and BMI was only detected in MD with atypical features, but with a markedly higher estimate (rg∼0.5). This result concurred with previous findings mainly using PRS or other samples.(14, 37, 38, 43) In contrast with the negative value found in all MD, we found positive correlations with educational attainment in many MD subtypes. However, this finding might be specific to the UKB cohort as previous research have shown that participation in mental health survey and other optional components is genetically correlated with higher education and intelligence.(44)
Taken together, our findings provide an improved understanding on heritable MD subtypes and overall genetic sharing between subtypes. These results have strong implications in the gene mapping strategies for MD. Current efforts predominantly aim to maximize samples size. The alternative strategy—to reduce phenotypic heterogeneity through more homogeneous phenotype— has not been fully evaluated, potentially due to theoretical and methodological challenges.(45) This strategy relies on the premise that “clinical heterogeneity in MD emerges from an aggregation of different underlying liabilities expressed through partially distinct biological pathways” (45) which, to the best knowledge, was not proven. Limited by a lack of large-scale dataset with deep phenotyping, prior studies were only able to focus on a few key subtypes.(5, 45) Our comprehensive report, by contrast, convincingly demonstrated genetic heterogeneity in MD, and thus forms a strong theoretical basis for this strategy. We further illustrated the potential of such strategy by performing stratified GWAS on each subtype. This yielded the identification of 47 independent genomic loci, a third of which were undetected in the latest MD GWAS with about 5- to 10-fold more cases than in our subtype-specific analyses. These results warrant further replications in large biobanks with consistent genotyping and phenotyping.
Here we used the UKB data which provide the unique opportunity to evaluate multiple subtypes with sufficient statistical power. We, however, note the following limitations in the context of interpreting the results. First, we were unable to study all MD subtypes, especially the treatment-related subtypes, as more refined clinical and treatment data would be required. We also acknowledge that the quality of phenotypic definitions varied across the subtypes studied, with those relying on self-reported and retrospective recalls of symptoms more compromised than the others. Finally, “healthy volunteer bias” was known for UKB(46) and likely to contribute to part of our results.
Etiological heterogeneity hinders treatment efficacy. Our finding of ubiquitous subtype heterogeneity within MD underscores the potential of drug development and treatment optimization for patient subgroups to achieve precision psychiatry.
URLs
Full protocol and scripts available via Github: https://github.com/Thuy-Dung-Nguyen/MD-subtypes;
UK Biobank Showcase User Guide (2017): http://biobank.ctsu.ox.ac.uk/crystal/crystal/exinfo/ShowcaseUserGuide.pdf;
UK Biobank-Mental health web-based questionnaire (2017) http://biobank.ctsu.ox.ac.uk/crystal/crystal/docs/mental_health_online.pdf;
GCTA-fastGWA: https://cnsgenomics.com/software/gcta/#fastGWA;
FUMA: https://fuma.ctglab.nl;
LDSC: https://github.com/bulik/ldsc;
HDL: https://github.com/zhenin/HDL;
Howard et al. 2019 MD GWAS summary results: https://datashare.is.ed.ac.uk/handle/10283/3203.
Data Availability
Data used for this manuscript is not available for public access. Please contact the corresponding author for questions regarding data. Code that was used to generate the results is available on GitHub link below.
Disclosures
PFS reports the following potentially competing financial interests. Current: Lundbeck (advisory committee, grant recipient), RBNC Therapeutics (advisory committee, stock ownership). CMB reports: Shire (grant recipient, Scientific Advisory Board member); Idorsia (consultant); Pearson (author, royalty recipient)
Acknowledgement
This research has been conducted using the UK Biobank Resource under Application Number 22224. This study was funded by the US NIMH grant (R01 MH123724) and the European Union’s Horizon 2020 research and innovation program under grant agreement No 847776. PFS was supported by the Swedish Research Council (Vetenskapsrådet, award D0886501), the Horizon 2020 Program of the European Union (COSYN, RIA grant agreement n° 610307), and US NIMH (U01 MH109528 and R01 MH077139). YL is in part supported by a 2018 NARSAD Young Investigator Grant from the Brain & Behaviour Research Foundation and US NIMH (R01 MH123724).
The computations were enabled by resources provided by the Swedish National Infrastructure for Computing (SNIC) at UPPMAX server partially funded by the Swedish Research Council through grant agreement no. 2018-05973.