#PAGE_PARAMS# #ADS_HEAD_SCRIPTS# #MICRODATA#

Studying individual risk factors for self-harm in the UK Biobank: A polygenic scoring and Mendelian randomisation study


Authors: Kai Xiang Lim aff001;  Frühling Rijsdijk aff001;  Saskia P. Hagenaars aff001;  Adam Socrates aff001;  Shing Wan Choi aff001;  Jonathan R. I. Coleman aff001;  Kylie P. Glanville aff001;  Cathryn M. Lewis aff001;  Jean-Baptiste Pingault aff001
Authors place of work: Social, Genetic and Developmental Psychiatry Centre, Institute of Psychiatry, Psychology & Neuroscience, King’s College London, London, United Kingdom aff001;  Icahn School of Medicine at Mount Sinai, New York, New York, United States of America aff002;  Division of Psychology and Language Sciences, University College London, London, United Kingdom aff003
Published in the journal: Studying individual risk factors for self-harm in the UK Biobank: A polygenic scoring and Mendelian randomisation study. PLoS Med 17(6): e32767. doi:10.1371/journal.pmed.1003137
Category: Research Article
doi: https://doi.org/10.1371/journal.pmed.1003137

Summary

Background

Identifying causal risk factors for self-harm is essential to inform preventive interventions. Epidemiological studies have identified risk factors associated with self-harm, but these associations can be subject to confounding. By implementing genetically informed methods to better account for confounding, this study aimed to better identify plausible causal risk factors for self-harm.

Methods and findings

Using summary statistics from 24 genome-wide association studies (GWASs) comprising 16,067 to 322,154 individuals, polygenic scores (PSs) were generated to index 24 possible individual risk factors for self-harm (i.e., mental health vulnerabilities, substance use, cognitive traits, personality traits, and physical traits) among a subset of UK Biobank participants (N = 125,925, 56.2% female) who completed an online mental health questionnaire in the period from 13 July 2016 to 27 July 2017. In total, 5,520 (4.4%) of these participants reported having self-harmed in their lifetime. In binomial regression models, PSs indexing 6 risk factors (major depressive disorder [MDD], attention deficit/hyperactivity disorder [ADHD], bipolar disorder, schizophrenia, alcohol dependence disorder, and lifetime cannabis use) predicted self-harm, with effect sizes ranging from odds ratio (OR) = 1.05 (95% CI 1.02 to 1.07, q = 0.008) for lifetime cannabis use to OR = 1.20 (95% CI 1.16 to 1.23, q = 1.33 × 10−35) for MDD. No systematic differences emerged between suicidal and non-suicidal self-harm. To further probe causal relationships, two-sample Mendelian randomisation (MR) analyses were conducted, with MDD, ADHD, and schizophrenia emerging as the most plausible causal risk factors for self-harm. The genetic liabilities for MDD and schizophrenia were associated with self-harm independently of diagnosis and medication. Main limitations include the lack of representativeness of the UK Biobank sample, that self-harm was self-reported, and the limited power of some of the included GWASs, potentially leading to possible type II error.

Conclusions

In addition to confirming the role of MDD, we demonstrate that ADHD and schizophrenia likely play a role in the aetiology of self-harm using multivariate genetic designs for causal inference. Among the many individual risk factors we simultaneously considered, our findings suggest that systematic detection and treatment of core psychiatric symptoms, including psychotic and impulsivity symptoms, may be beneficial among people at risk for self-harm.

Keywords:

cannabis – depression – Medical risk factors – Schizophrenia – Bipolar disorder – suicide – Self harm – ADHD

Introduction

Self-harm refers to any act of self-injury or self-poisoning carried out by an individual, regardless of intention or motivation [1]. It can be further categorised into suicidal self-harm (SSH) and non-suicidal self-harm (NSSH). To illustrate, NSSH is an intentional act of harming oneself without suicidal intention, such as self-cutting, self-burning, and self-poisoning [2]. SSH, on the other hand, is associated with the intention to die, such as overdosing, severe cutting, and hanging/strangulation, and is often more lethal than NSSH [3]. According to a meta-analysis, the cross-national prevalence rate for NSSH peaks during adolescence (17.3%) and decreases among adults (5.5%) [4]. For SSH, the cross-national prevalence rate is also highest among adolescents (9.7%) [5] and drops among adults (2.7%) [6]. Recently, the fifth edition of the Diagnostic and Statistical Manual of Mental Disorders (DSM-5) divided self-harm along the dimension of suicidal intention as 2 separate conditions for further study [7]. The distinction between SSH and NSSH may facilitate investigations of the aetiology and heterogeneity of self-harm.

Multiple individual vulnerabilities and traits can potentially lead to self-harm, such as psychiatric illnesses [8], substance use [911], cognitive abilities [12], personality traits [13], and physical traits [14]. The stress–diathesis model of suicidal behaviour proposes that the risk for suicidal behaviours is determined not just by stressors such as major psychiatric disorders, but also by a diathesis, rooted in susceptibilities such as cognitive abilities and impulsivity [1518]. Although associations between these risk factors and self-harm have been shown in numerous observational studies, causality is difficult to infer reliably. Genetically informed designs can help strengthen causal inference [19]. A polygenic score (PS) is a single individual-level score computed for a given trait, weighted using summary statistics from an independent genome-wide association study (GWAS) for that particular trait. A PS for an individual risk factor (e.g., schizophrenia) can be regarded as a genetic proxy for this risk factor [20]. To illustrate, if schizophrenia is causally related to self-harm, a PS for schizophrenia should also be associated with self-harm. A significant association between the PS for schizophrenia and self-harm can be regarded as an initial indication of a possible causal relationship between the 2 traits. The PS approach can be construed as a first step in a series of genetically informed methods to investigate the aetiology of complex phenotypes, with follow-up steps including Mendelian randomisation (MR), discussed below [1921].

In previous studies, a PS for major depressive disorder (MDD) predicted SSH in 2 clinical samples [22,23] and 1 non-clinical sample [24]. However, this finding was not replicated in a family-based sample after applying Bonferroni correction (although, at the nominal p-value threshold, the PS for MDD generated from the first Psychiatric Genomics Consortium GWAS for MDD predicted SSH among MDD cases in this sample) [25]. A PS for depressive symptoms predicted SSH but not NSSH in a twin sample [26]. A PS for schizophrenia predicted SSH among offspring of suicide attempters [25] and in a population sample [27], but not in another clinical sample [28]. A PS for bipolar disorder predicted SSH in 1 clinical sample [29] but did not predict SSH nor NSSH among offspring of suicide attempters [25] or relatives of bipolar disorder patients [30]. A recently published study showed that PSs for cannabis use and MDD predicted self-harm in the UK Biobank sample [31].

The aforementioned PS studies with mixed results were limited in several ways. First, these studies mostly focused on PSs for psychiatric disorders or symptoms, and did not include potential individual risk factors from other domains, such as substance use [911], cognitive abilities [12], personality traits [13], and physical traits [14]. Second, with 2 exceptions [26,30], none of the studies investigated SSH and NSSH simultaneously. Third, these studies have a mixture of clinical and non-clinical samples with varying sample sizes ranging from 224 individuals [28] to 63,054 individuals [31], making any comparison difficult. Furthermore, none of these studies have implemented multivariate analyses with multiple PSs to better estimate their unique effect.

A caveat of the PS method is its proneness to unmediated (or horizontal) pleiotropy, arising from the inclusion of thousands of genetic variants [20]. Unmediated pleiotropy exists when a genetic variant associated with an exposure causes the outcome through an alternative pathway, instead of via the exposure. Unmediated pleiotropy can generate associations between PSs and outcomes in the absence of a causal relationship between the risk factors indexed by the PS, and the outcome. MR can more stringently address unmediated pleiotropy and further strengthen causal inference. In MR, individual genetic variants associated with an exposure of interest are used as instrumental variables to infer causality between exposure and outcome. A number of complementary analyses, further detailed in the Methods section, can be implemented to account for pleiotropy [19]. To date, to our knowledge, there is no published MR study that focuses on any risk factor for self-harm.

The current study addressed the aforementioned limitations by systematically using 24 PSs as proxies for risk factors from different domains to predict both NSSH and SSH, using a population-based sample of 125,925 individuals from the UK Biobank. Selected PSs reflect in part the 2 dimensions of the stress–diathesis model of suicidal behaviour [1518]. We also conducted follow-up MR analyses to strengthen causal inference. This study builds on previous genetic studies using UK Biobank samples that have identified 3 genome-wide significant loci for suicidality [32] and found significant genetic correlations between SSH and several psychiatric disorders [33].

Methods

This study did not have a prespecified analysis plan.

Participants

The participants are a subset of the UK Biobank (http://www.ukbiobank.ac.uk). A total of 157,358 participants completed an online mental health questionnaire in a period from 13 July 2016 to 27 July 2017, including questions regarding their lifetime symptoms of mental disorders [34]. After genotyping and quality control processing (see S1 Method), the final sample size was 125,925 individuals (56.2% female), with ages ranging from 48 to 82 years (mean = 65.88, SD = 7.69).

The UK Biobank received ethical approval from the North West–Haydock Research Ethics Committee (REC reference 11/NW/0382). The current study was conducted under UK Biobank application 18177. Appropriate informed consent was obtained from the participants.

Defining self-harm phenotypes

To know whether the participants have ever self-harmed, participants were asked, ‘Have you deliberately harmed yourself, whether or not you meant to end your life?’ To ascertain whether their self-harm episodes were NSSH or SSH, they were asked, ‘Have you harmed yourself with the intention to end your life?’ In both questions, responses of ‘Prefer not to answer’ (0.43%) were recoded as missing values. A flowchart depicting exclusion of participants and the number of participants who answered each question is shown in Fig 1. We used the final sample with complete data (i.e., without missing genotype or self-harm data) for subsequent analyses.

Participant flowchart.
Fig. 1. Participant flowchart.
MHQ, mental health questionnaire; PS, polygenic score.

Statistical analyses

All statistical analyses were conducted in Linux environment using R version 3.5.0 [35].

PS analyses

PSs of UK Biobank participants were generated using PRSice-2 [36] based on their genotype data and 24 publicly available summary data from GWASs (see Table 1) selected based on the following criteria. First, we selected GWASs indexing individual vulnerabilities and traits that can potentially increase the risk of self-harm, including mental health vulnerabilities (e.g., MDD) [37], cognitive abilities (e.g., education attainment) [38], personality traits (e.g., neuroticism) [39], substance use phenotypes (e.g., cannabis use) [40], and physical traits (e.g., BMI) [14]. Second, we selected GWASs that only included participants of European ancestry and did not include UK Biobank participants (to avoid overlapping between the discovery and target samples). Finally, we excluded GWASs with effective sample sizes less than N = 15,000 to limit the use of underpowered PSs.

Tab. 1. Single and multiple PS prediction of self-harm.
Single and multiple PS prediction of self-harm.

Each participant had 24 PSs, which were each calculated as the sum of alleles associated with their respective phenotypes, weighted by their effect sizes, with p-values less than a threshold pT < 0.3. Similar to previous PS studies [27,31], we used 1 single liberal threshold, which allowed us to limit multiple testing while retaining the predictive ability of included PSs [58]. Clumping was used to remove SNPs in linkage equilibrium (r2 < 0.1 within a 250-kb window). All PSs in the final analytical sample were standardised.

Single PS binomial logistic regression

For each PS, a binomial logistic regression was conducted to test whether it predicted self-harm (i.e., ‘self-harmed’ versus ‘never self-harmed’).

Multiple PS binomial logistic regression

All PSs significantly associated with self-harm in single PS binomial logistic regressions after correction for multiple testing were then jointly modelled in a multivariate binomial logistic regression model to assess their unique effects.

Single PS multinomial logistic regressions

To investigate whether each PS differentially predicted NSSH versus SSH, we fitted a series of multinomial logistic regression models. We first compared each of the NSSH and SSH groups to the never-self-harmed group (i.e., ‘never self-harmed’ as the reference group). We then directly compared the NSSH and SSH groups by testing a model with NSSH as the reference group.

Covariates and multiple testing

All regression models were controlled for sex, age, and population stratification (by including assessment centre, genotyping batch, and the first 6 principal components as covariates in the models). To control for multiple testing in single PS binomial and multinomial regressions, we employed the false discovery rate (FDR) method [59], which controls the expected proportion of false positives among the rejected hypotheses. We used q < 0.05 as the significance threshold.

MR analyses

All MR analyses were conducted using R package TwoSampleMR [60]. Risk factors for which the PS significantly predicted self-harm were selected as exposures for follow-up MR analyses. For self-harm in the UK Biobank sample as the outcome for MR analyses, we obtained GWAS summary statistics from Neale Lab (http://www.nealelab.is/uk-biobank). Exposure SNPs that passed the p-value threshold of p < 5 × 10−5 were selected as instrumental variables. A liberal threshold was used to ensure that enough variants were available for all risk factors, including those with few genome-wide significant SNPs (e.g., ADHD). The strategy entails potential weak instrument bias. In two-sample MR, the resulting bias is towards the null, making estimates more conservative (see below for how this was dealt with) [61]. For more details on clumping and harmonisation, see Text A in S2 Method.

We selected 4 MR methods that have different strengths and limitations. We conducted univariable MR using (i) inverse variance weighting (IVW)—the most powerful method (but cannot account for directional pleiotropy) [62]; (ii) robust adjusted profile score (RAPS)—to account for the selection of weak instruments [63]; (iii) weighted median—more robust to directional pleiotropy than IVW, and more robust to individual genetic variants with outlying causal estimates than IVW and MR–Egger regression [64]; and (iv) MR–Egger regression—the significance of its intercept term informs on the presence of directional pleiotropy [65].

We calculated the I2GX statistics for the MR instruments to assess the suitability of carrying out the MR–Egger regression method [66]. We also implemented MR Steiger filtering [67] to address the possibility of reverse causation (see Text B in S2 Method). Finally, similar to PS analyses, exposures that were significant in univariable MR were assessed for their independent effect in a multivariable MR model using the IVW method. GWAS summary statistics for SSH are also available. We thus repeated MR analyses for SSH only.

For PS analyses, we conducted further complementary analyses excluding cases with MDD and schizophrenia diagnoses to investigate the effect of genetic liability on self-harm with the influence of these diagnoses excluded. We also calculated risk ratios of self-harm for medicated and non-medicated cases compared to those with median PS in the general population (see S1 Appendix).

This study is reported as per the Strengthening the Reporting of Observational Studies in Epidemiology (STROBE) guideline (see S1 Checklist).

Results

Descriptive statistics

Fig 1 shows the number of participants who never self-harmed, self-harmed, engaged in SSH, and engaged in NSSH. S1 Table shows the descriptive statistics of each group.

PS analyses

Single PS binomial logistic regression

Table 1 and Fig 2 show results from 24 single PS binomial logistic regression tests, using each PS as predictor. Out of the 24 PSs, 10 PSs were significant predictors of self-harm at the nominal level (p < 0.05). After FDR correction, 6 PSs had q-value < 0.05. In order of decreasing effect sizes, they are PSs for MDD, schizophrenia, ADHD, bipolar disorder, alcohol dependence disorder (ALC), and lifetime cannabis use, with effect sizes ranging from OR = 1.197 (95% CI 1.164 to 1.230) for MDD to OR = 1.045 (95% CI 1.016 to 1.074) for lifetime cannabis use. S1 Fig shows the pseudo R2 plots of these 6 PSs in accounting for the variance in self-harm.

Single and multiple PS analysis of self-harm.
Fig. 2. Single and multiple PS analysis of self-harm.
Estimates from single PS regressions and multiple PS regression for PSs of 24 risk factors in predicting self-harm, ordered according to decreasing effect sizes. ‘(1)’ indicates not significant in multiple PS regression. ‘(2)’ indicates not significant after FDR correction. ADHD, attention deficit/hyperactivity disorder; FDR, false discovery rate; IRT, item response theory; MDD, major depressive disorder; OR, odds ratio; PS, polygenic score.

Multiple PS binomial logistic regression

In the multiple PS model, all PSs except the PS for ALC had an independent effect on self-harm, as shown in Table 1 and Fig 2. Because this model controls for the effects of other PSs, the effect sizes of these PSs are slightly lower than those in the single PS binomial logistic regressions, ranging from OR = 1.032 (95% CI 1.002 to 1.061) for bipolar disorder to OR = 1.155 (95% CI 1.123 to 1.189) for MDD. These PSs were only weakly correlated, suggesting that multicollinearity was not an issue (see S2 Table).

Fig 3 shows the independent effects of the 5 significantly predictive PSs (MDD, schizophrenia, ADHD, bipolar disorder, and lifetime cannabis use) corresponding to the mean PS in each quintile, and their combined effect on the predicted risk of self-harm. To illustrate, for the independent effect of MDD, the predicted risk of self-harm for those whose PS for MDD was in the lowest 20% of the sample was 3.10% (95% CI 2.95% to 3.26%), as compared to 4.58% (95% CI 4.37% to 4.79%) for those whose PS for MDD was in the highest 20%. For combined effect prediction, among participants whose 5 PSs were all in the lowest 20% of the sample, their combined predicted risk of self-harm was 2.18% (95% CI 1.99% to 2.36%). For those whose 5 PSs were all in the top 20% of the sample, their combined predicted risk of self-harm was 6.47% (95% CI 6.00% to 6.95%; also see S3 Method, S3 and S4 Tables).

Predicted risk of self-harm for each of the 5 significantly predictive PSs and their combined predicted risk.
Fig. 3. Predicted risk of self-harm for each of the 5 significantly predictive PSs and their combined predicted risk.
Plot of the independently predicted risk of self-harm for each of the 5 significantly predictive PSs (for ADHD, bipolar disorder, lifetime cannabis use, MDD, and schizophrenia) and their combined predicted risk obtained from multivariate logistic regression. For independent prediction using a single PS, the mean of that particular PS in each quintile was used as predictor while the values of other PSs were set to 0 (i.e., the centred mean). For combined prediction, for each quintile, 5 mean values of the PS (mean PS for ADHD, bipolar disorder, lifetime cannabis use, MDD, and schizophrenia) were entered simultaneously as predictors. ADHD, attention deficit/hyperactivity disorder; MDD, major depressive disorder; PS, polygenic score.

Single PS multinomial logistic regressions

S5 Table shows results from 24 multinomial logistic regression tests, using PSs as predictor for 3 possible outcomes: never self-harmed, NSSH, and SSH. When ‘never self-harmed’ was used as the reference group, PSs for bipolar disorder, lifetime cannabis use, and extreme BMI predicted SSH but not NSSH, with q < 0.05. However, when NSSH was set as the reference group in order to directly compare NSSH versus SSH, none of the PSs significantly distinguished between NSSH and SSH.

MR analyses

Table 2 shows the results from MR analyses: ADHD, ALC, bipolar disorder, lifetime cannabis use, MDD, and schizophrenia were exposures in 6 separate univariable MR analyses, with self-harm as the outcome. Out of these 6 exposures, MDD, ADHD, and schizophrenia had MR estimates with p-values < 0.05. For the other exposures, none of their MR estimates had p < 0.05.

Tab. 2. Univariable and multivariable MR analyses.
Univariable and multivariable MR analyses.

For MDD, despite having the strongest IVW (β = 0.008, 95% CI 0.005 to 0.011, p = 2.84 × 10−8), MR RAPS (β = 0.008, 95% CI 0.005 to 0.011, p = 1.24 × 10−7), and weighted median (β = 0.006, 95% CI 0.001 to 0.011, p = 0.013) estimates among the 3 exposures with significant MR estimates, the MR–Egger regression estimate was not significant. For ADHD and schizophrenia, all MR estimates were significant. In all cases, the I2GX values were >90% (see S6 Table), suggesting that regression dilution has not substantially impacted the MR–Egger analyses [66].

Significance of intercept terms in MR–Egger analyses indicates the presence of pleiotropy. None of the MR–Egger intercept terms were significant, except for MDD (p = 0.023). MR Steiger filtering could only be conducted for MDD and schizophrenia (see Text B in S2 Method). Selected SNPs of these 2 exposures are more predictive of the respective exposures than self-harm, suggesting that reverse causation is unlikely to explain our findings. Hence, no SNPs were excluded following the MR Steiger filtering.

When ADHD, MDD, and schizophrenia were included as exposures in multivariable IVW MR analysis, only MDD (β = 0.011, 95% CI 0.007 to 0.015, p = 1.04 × 10−12) and schizophrenia (β = 0.002, 95% CI 4.00 × 10−5 to 0.004, p = 0.002) remained as independent predictors of self-harm. Due to the potential presence of pleiotropy between MDD and self-harm, another multivariable IVW MR analysis was conducted with only ADHD and schizophrenia as exposures. Both ADHD (β = 0.003, 95% CI 0.001 to 0.005, p = 2.21 × 10−4) and schizophrenia (β = 0.003, 95% CI 0.002 to 0.004, p = 7.60 × 10−7) were significant predictors in this model.

We repeated the multivariable IVW MR analysis using a genome-wide threshold in selecting the instruments for schizophrenia. Similar results were observed, whereby only MDD (β = 0.010, 95% CI 0.006 to 0.014, p = 5.48 × 10−8) and schizophrenia (β = 0.003, 95% CI 0.001 to 0.005, p = 0.002) remained as significant predictors for self-harm (see S7 Table). Results for MR analyses with SSH as the outcome are presented in S8 Table. Only MR estimates for schizophrenia were significant (possibly due to the lower prevalence of SSH, leading to a corresponding loss of power in the SSH GWAS).

In complementary analyses that excluded cases with MDD and schizophrenia diagnoses, PSs for MDD and schizophrenia predicted self-harm in a healthy screened cohort, indicating that genetic liabilities can predict self-harm when the influence of diagnoses is excluded (See Table A in S1 Appendix). As shown in Fig 4, individuals diagnosed with schizophrenia and MDD appear to be at much larger risk for self-harm than the rest of the population. For MDD, medicated cases were at higher risk of self-harm than non-medicated cases, which was not the case for schizophrenia.

Risk ratios of self-harm in the general population and in non-medicated and medicated MDD and schizophrenia cases.
Fig. 4. Risk ratios of self-harm in the general population and in non-medicated and medicated MDD and schizophrenia cases.
Risk ratios of self-harm in the general population, non-medicated cases, and medicated cases were compared to those with median polygenic score (11th quantile) for schizophrenia (A) and MDD (B). Out of 177 schizophrenia cases in the final analytical sample, 89 (50.3%) of them were medicated. Out of 34,680 MDD cases in the final analytical sample, 7,852 (22.6%) of them were medicated. MDD, major depressive disorder; RR, risk ratio.

Discussion

We used multiple PSs as genetic proxies to systematically investigate multiple individual vulnerabilities and traits as risk factors for self-harm in a population sample. In PS analyses, we identified 6 risk factors (MDD, schizophrenia, ADHD, bipolar disorder, ALC, and lifetime cannabis use) that predicted self-harm. Five among them (all except for ALC) remained significant in a multiple PS regression. We found little evidence of differential prediction for SSH versus NSSH. In follow-up MR analyses, MDD, schizophrenia, and ADHD emerged as plausible causal risk factors for self-harm, despite evidence of unmediated pleiotropy for MDD.

Insights into the aetiology of self-harm

Results from our PS methods corroborated previous observational findings where MDD [8], schizophrenia [68], ADHD [69], bipolar disorder [8], and ALC [10] were phenotypically associated with self-harm. Our results are also consistent with positive associations found in PS studies for MDD [2224], schizophrenia [25], and bipolar disorder [29]. Previous mixed findings for these PSs may have stemmed from lack of power, as sample sizes for those studies varied widely. The current study adds the PSs for ADHD and ALC as novel PSs genetically associated with self-harm. However, when controlling for other PSs, the PS for ALC did not significantly predict self-harm. This finding may suggest that the genetic liability for ALC does not independently predict self-harm when the effects of genetic liability for MDD, bipolar disorder, schizophrenia, ADHD, and lifetime cannabis use are accounted for. For example, ALC may be a marker for a true predictor such as impulsivity, which is more efficiently captured in the PS for ADHD [70]. Alternatively, null findings for ALC may also be due to a lack of power compared to other PSs. Hence, we cannot completely rule out that the PS for ALC has an independent effect on self-harm and there may be a corresponding causal effect of ALC on self-harm.

Most of the PSs that predicted self-harm in this study relate to psychiatric conditions, which confirms the prominence of psychiatric conditions in the aetiology of self-harm [71]. Cognitive traits, physical traits, and personality traits were not found to be associated with self-harm using the PS approach, although previous observational findings found significant phenotypic associations for these domains [1214]. The absence of significant findings in our study is unlikely to be solely due to lack of power, given that GWASs for some of these traits are more powerful than GWASs for psychiatric conditions (e.g., BMI and education attainment). These findings suggest that these traits and vulnerabilities are unlikely to have (strong) causal effects on self-harm.

Our MR analyses further supported the role of MDD, ADHD, and schizophrenia in the aetiology of self-harm. Psychiatric conditions such as schizophrenia and MDD are among possible stressors in the stress–diathesis model of suicidal behaviour, whereas ADHD, which includes symptoms of impulsivity, may be considered as part of the diathesis in the model [1518]. Our results showed that traits pertaining to both the stress and diathesis dimensions are contributing to the aetiology of self-harm. An intriguing finding is the presence of significant pleiotropy in the case of MDD. Rather than signifying that MDD does not have a causal effect on self-harm, this may reflect a measurement issue. Indeed, one of the diagnostic criteria for MDD is related to having suicidal thoughts and attempts, which could artificially introduce a pleiotropic effect [7]. To deal with this issue, future studies could rely on a GWAS for MDD excluding the diagnostic criterion related to suicidal thoughts and attempts. This might also explain why, in the multivariable MR analysis, the effect of ADHD was no longer significant—as we partially controlled for self-harm—whereas it was significant when considering only ADHD and schizophrenia.

The current study found mixed results for whether there are distinct aetiologies for SSH and NSSH. Most PSs that predicted self-harm also predicted both SSH and NSSH, except the PSs for bipolar disorder, lifetime cannabis use, and extreme BMI, which only predicted SSH but not NSSH. However, in a formal test comparing NSSH and SSH, the estimates of these 3 risk factors were not significantly different between NSSH and SSH. Hence, our findings do not provide evidence for marked differences in aetiology between SSH and NSSH.

Clinical implications

This study suggests that individual vulnerabilities and traits underlying self-harm most likely relate to psychiatric conditions (i.e., MDD and schizophrenia), rather than other domains such as personality traits. Hence, treatments focusing on the core symptoms of these psychiatric conditions are important in preventing or addressing the risk of self-harm. Findings from PS analyses suggest that genetic liabilities for these conditions increase the likelihood of self-harm even in those not clinically diagnosed, suggesting that subthreshold symptoms of these psychiatric conditions may increase the risk of self-harm. Clinicians may want to systematically test for such symptoms in self-harming patients. Future investigations could test whether drugs for such core conditions can be repurposed for treating self-harming patients with either full-blown or subthreshold conditions. For example, prescription of methylphenidate for ADHD treatment was found to be associated with reduction of suicide risk [72]. As a note of caution, medicated schizophrenia cases were not at less risk of self-harm than non-medicated patients, whereas medicated MDD patients were at substantially higher risk for self-harm in the present sample. This could be due to medicated patients having more severe symptoms than non-medicated patients, or to adverse effects of medication, particularly for MDD, where suicidality might be an adverse effect of antidepressant treatment [73].

Limitations

To avoid the overlapping of discovery and target samples, we excluded GWASs that contain the UK Biobank sample, resulting in our selecting older GWASs for generating PSs in some cases. This might have led to non-significant findings due to lack of power. The results should be generalised with caution because the UK Biobank is not representative of the UK population, as participants are more educated, older, wealthier, and healthier [74]. The questions asked in the mental health questionnaire were retrospective and their formulation led to an exclusive dichotomy between NSSH or SSH, although some might have engaged in both NSSH and SSH at different times.

PSs in this study explained only a small amount of the variance in self-harm. The maximum variance explained by any single PS was 0.4%, for the MDD PS. Importantly, such prediction of self-harm is likely indirect, i.e., we expect the PS for a given risk factor to influence this risk factor first, and then self-harm only to the extent that this risk factor is a causal risk for self-harm. As a result, given that, for example, the PS for MDD only explains 1.9% of the variance in the liability for MDD [37], it is unsurprising that its effect on self-harm is much lower. In addition, if the PS for MDD were to capture the full heritability of MDD (which is 37%, according to a meta-analysis of twin studies [75]), the explanatory power of the MDD PS could be expected to proportionally increase up to 7.8%. In turn, this would contribute to explaining the gap between the currently estimated heritability of self-harm in family-based studies (up to about 50% [76]) and the low percentage of variance explained in our analyses. Any remaining gap would indicate that other genetically influenced risk factors for self-harm remain to be identified beyond those uncovered in the present study.

The lack of evidence showing significant aetiological differences between NSSH and SSH may also be due to lack of power for subset analyses. However, the effect sizes of the PSs for both NSSH and SSH appeared similar overall, suggesting that the aetiology of NSSH and SSH may be similar. For MR analyses, although a range of sensitivity analyses were conducted to account for pleiotropy and weak instrument bias, they cannot completely rule out horizontal pleiotropy. The possibility thus remains that significant associations are non-causal. On the other hand, risk factors that were non-significant in MR analyses should not be entirely ruled out as potential causal risk factors for self-harm. Future studies with more powerful genetic instruments may uncover new risk factors that we were unable to identify.

Conclusion

Among 24 PSs used as genetic proxies for vulnerabilities and traits possibly associated with self-harm, we found that PSs for MDD, schizophrenia, ADHD, bipolar disorder, ALC, and lifetime cannabis use were statistically significant. After a series of complementary analyses to further strengthen causal inference, schizophrenia emerged as the most plausible causal risk factor, followed by MDD and ADHD. Detection and treatment of core symptoms of these conditions, such as psychotic or impulsivity symptoms, may benefit self-harming patients.

Supporting information

S1 Appendix [docx]
Complementary analyses.

S1 Checklist [docx]
STROBE checklist.

S1 Fig [docx]
Pseudo plots of 6 PSs in predicting self-harm.

S1 Method [docx]
Genotyping, imputation, and quality control.

S2 Method [docx]
Mendelian randomisation.

S3 Method [docx]
Predicted risk of self-harm.

S1 Table [docx]
Descriptive statistics for each group of self-harm-related phenotypes.

S2 Table [docx]
Correlations between PSs that were significant in single PS regression in predicting self-harm.

S3 Table [docx]
Single PS prediction of self-harm risk.

S4 Table [docx]
Multiple PS prediction of self-harm risk.

S5 Table [docx]
Multinomial regression models in predicting SSH and NSSH.

S6 Table [docx]
2GX statistics for the MR instruments in univariable MR analyses.

S7 Table [docx]
Multivariable MR with self-harm as outcome, using a genome-wide -value threshold ( < 5 × 10) to select the instruments.

S8 Table [docx]
Univariable MR results with SSH as the outcome.


Zdroje

1. Hawton K, Harriss L, Hall S, Simkin S, Bale E, Bond A. Deliberate self-harm in Oxford, 1990–2000: a time of change in patient characteristics. Psychol Med. 2003;33:987–95. doi: 10.1017/s0033291703007943 12946083

2. McManus S, Gunnell D, Cooper C, Bebbington PE, Howard LM, Brugha T, et al. Prevalence of non-suicidal self-harm and service contact in England, 2000–14: repeated cross-sectional surveys of the general population. Lancet Psychiatry. 2019;6:573–81. doi: 10.1016/S2215-0366(19)30188-9 31175059

3. Andover MS, Gibb BE. Non-suicidal self-injury, attempted suicide, and suicidal intent among psychiatric inpatients. Psychiatry Res. 2010;178:101–5. doi: 10.1016/j.psychres.2010.03.019 20444506

4. Swannell SV., Martin GE, Page A, Hasking P, St John NJ. Prevalence of nonsuicidal self-injury in nonclinical samples: systematic review, meta-analysis and meta-regression. Suicide Life Threat Behav. 2014;44:273–303. doi: 10.1111/sltb.12070 24422986

5. Evans E, Hawton K, Rodham K, Psychol C, Deeks J. The prevalence of suicidal phenomena in adolescents: a systematic review of population-based studies. Suicide Life Threat Behav. 2005;35:239–50. doi: 10.1521/suli.2005.35.3.239 16156486

6. Nock MK, Borges G, Bromet EJ, Alonso J, Angermeyer M, Beautrais A, et al. Cross-national prevalence and risk factors for suicidal ideation, plans and attempts. Br J Psychiatry. 2008;192:98–105. doi: 10.1192/bjp.bp.107.040113 18245022

7. American Psychiatric Association. Diagnostic and statistical manual of mental disorders (DSM-5). Washington (DC): American Psychiatric Association; 2013.

8. Vaughn MG, Salas-Wright CP, DeLisi M, Larson M. Deliberate self-harm and the nexus of violence, victimization, and mental health problems in the United States. Psychiatry Res. 2015;225:588–95. doi: 10.1016/j.psychres.2014.11.041 25500323

9. Borges G, Bagge CL, Orozco R. A literature review and meta-analyses of cannabis use and suicidality. J Affect Disord. 2016; 195:63–74. doi: 10.1016/j.jad.2016.02.007 26872332

10. Darvishi N, Farhadi M, Haghtalab T, Poorolajal J. Alcohol-related risk of suicidal ideation, suicide attempt, and completed suicide: a meta-analysis. PLoS ONE. 2015;10:e0126870. doi: 10.1371/journal.pone.0126870 25993344

11. Evins AE, Korhonen T, Kinnunen TH, Kaprio J. Prospective association between tobacco smoking and death by suicide: a competing risks hazard analysis in a large twin cohort with 35-year follow-up. Psychol Med. 2017;47:2143–54. doi: 10.1017/S0033291717000587 28399944

12. Rehkopf DH, Buka SL. The association between suicide and the socio-economic characteristics of geographical areas: a systematic review. Psychol Med. 2006;36:145–57. doi: 10.1017/S003329170500588X 16420711

13. Brezo J, Paris J, Turecki G. Personality traits as correlates of suicidal ideation, suicide attempts, and suicide completions: a systematic review. Acta Psychiatr Scand. 2006;113:180–206. doi: 10.1111/j.1600-0447.2005.00702.x 16466403

14. Perera S, Eisen RB, Dennis BB, Bawor M, Bhatt M, Bhatnagar N, et al. Body mass index is an important predictor for suicide: results from a systematic review and meta-analysis. Suicide Life Threat Behav. 2016;46:697–736. doi: 10.1111/sltb.12244 27094229

15. Mann JJ, Waternaux C, Haas GL, Malone KM. Toward a clinical model of suicidal behavior in psychiatric patients. Am J Psychiatry. 1999;156:181–9. 9989552

16. Van Heeringen K, Mann JJ. The neurobiology of suicide. Lancet Psychiatry. 2014;1:63–72. doi: 10.1016/S2215-0366(14)70220-2 26360403

17. Van Heeringen K. Stress-diathesis model of suicidal behavior. In: Dwivedi Y, editor. The neurobiological basis of suicide. Boca Raton (FL): CRC Press; 2012. pp. 113–123. doi: 10.1201/b12215

18. McGirr A, Turecki G. The relationship of impulsive aggressiveness to suicidality and other depression-linked behaviors. Curr Psychiatry Rep. 2007;9:460–6. doi: 10.1007/s11920-007-0062-2 18221625

19. Pingault J-B, O’Reilly PF, Schoeler T, Ploubidis GB, Rijsdijk F, Dudbridge F. Using genetic data to strengthen causal inference in observational research. Nat Rev Genet. 2018;19:566–80. doi: 10.1038/s41576-018-0020-3 29872216

20. Gage SH, Davey Smith G, Ware JJ, Flint J, Munafò MR. G = E: what GWAS can tell us about the environment. PLoS Genet. 2016;12:e1005765. doi: 10.1371/journal.pgen.1005765 26866486

21. Pingault J-B, Cecil C, Murray J, Munafo M, Viding E. Causal inference in psychopathology: a systematic review of Mendelian randomisation studies aiming to identify environmental risk factors for psychopathology. Psychopathol Rev. 2017;a4:4–25. doi: 10.5127/pr.038115

22. Mullins N, Perroud N, Uher R, Butler AW, Cohen-Woods S, Rivera M, et al. Genetic relationships between suicide attempts, suicidal ideation and major psychiatric disorders: a genome-wide association and polygenic scoring study. Am J Med Genet B Neuropsychiatr Genet. 2014;165:428–37. doi: 10.1002/ajmg.b.32247 24964207

23. Mullins N, Bigdeli TB, Børglum AD, Coleman JRIR, Demontis D, Mehta D, et al. GWAS of suicide attempt in psychiatric disorders and association with major depression polygenic risk scores. Am J Psychiatry. 2019;176:651–60. doi: 10.1176/appi.ajp.2019.18080957 31164008

24. Levey DF, Polimanti R, Cheng Z, Zhou H, Nuñez YZ, Jain S, et al. Genetic associations with suicide attempt severity and genetic overlap with major depression. Transl Psychiatry. 2019;9:22. doi: 10.1038/s41398-018-0340-2 30655502

25. Sokolowski M, Wasserman J, Wasserman D. Polygenic associations of neurodevelopmental genes in suicide attempt. Mol Psychiatry. 2016;21:1381–90. doi: 10.1038/mp.2015.187 26666204

26. Maciejewski DF, Renteria ME, Abdellaoui A, Medland SE, Few LR, Gordon SD, et al. The association of genetic predisposition to depressive symptoms with non-suicidal and suicidal self-injuries. Behav Genet. 2017;47:3–10. doi: 10.1007/s10519-016-9809-z 27590903

27. Laursen TM, Trabjerg BB, Mors O, Børglum AD, Hougaard DM, Mattheisen M, et al. Association of the polygenic risk score for schizophrenia with mortality and suicidal behavior—a Danish population-based study. Schizophr Res. 2017;184:122–7. doi: 10.1016/j.schres.2016.12.001 27939829

28. Bani-Fatemi A, Tasmim S, Wang KZ, Warsh J, Sibille E, De Luca V. No interaction between polygenic scores and childhood trauma in predicting suicide attempt in schizophrenia. Prog Neuropsychopharmacol Biol Psychiatry. 2019;89:169–73. doi: 10.1016/j.pnpbp.2018.08.027 30149093

29. Wiste A, Robinson EB, Milaneschi Y, Meier S, Ripke S, Clements CC, et al. Bipolar polygenic loading and bipolar spectrum features in major depressive disorder. Bipolar Disord. 2014;16:608–16. doi: 10.1111/bdi.12201 24725193

30. Wilcox HC, Fullerton JM, Glowinski AL, Benke K, Kamali M, Hulvershorn LA, et al. Traumatic stress interacts with bipolar disorder genetic risk to increase risk for suicide attempts. J Am Acad Child Adolesc Psychiatry. 2017;56:1073–80. doi: 10.1016/j.jaac.2017.09.428 29173741

31. Hodgson K, Coleman JRI, Hagenaars SP, Purves KL, Glanville K, Choi SW, et al. Cannabis use, depression and self-harm: phenotypic and genetic relationships. Addiction. 2019;115:482–92. doi: 10.1111/add.14845 31833150

32. Strawbridge RJ, Ward J, Ferguson A, Graham N, Shaw RJ, Cullen B, et al. Identification of novel genome-wide associations for suicidality in UK Biobank, genetic correlation with psychiatric disorders and polygenic association with completed suicide. EBioMedicine. 2019;41:517–25. doi: 10.1016/j.ebiom.2019.02.005 30745170

33. Ruderfer DM, Walsh CG, Aguirre MW, Tanigawa Y, Ribeiro JD, Franklin JC, et al. Significant shared heritability underlies suicide attempt and clinically predicted probability of attempting suicide. Mol Psychiatry. 2019 Jan 4. doi: 10.1038/s41380-018-0326-8 30610202

34. Davis KAS, Coleman JRI, Adams M, Allen N, Breen G, Cullen B, et al. Mental health in UK Biobank—development, implementation and results from an online questionnaire completed by 157 366 participants: a reanalysis. BJPsych Open. 2020;6:e18. doi: 10.1192/bjo.2019.100 32026800

35. R Core Team. R: a language and environment for statistical computing. Version 3.5.0. Vienna: R Foundation for Statistical Computing; 2017 [cited 2020 May 8]. https://www.r-project.org/.

36. Euesden J, Lewis CM, O’Reilly PF. PRSice: polygenic risk score software. Bioinformatics. 2015;31:1466–8. doi: 10.1093/bioinformatics/btu848 25550326

37. Wray NR, Ripke S, Mattheisen M, Trzaskowski M, Byrne EM, Abdellaoui A, et al. Genome-wide association analyses identify 44 risk variants and refine the genetic architecture of major depression. Nat Genet. 2018;50:668–81. doi: 10.1038/s41588-018-0090-3 29700475

38. Rietveld CA, Esko T, Davies G, Pers TH, Turley P, Benyamin B, et al. Common genetic variants associated with cognitive performance identified using the proxy-phenotype method. Proc Natl Acad Sci U S A. 2014;111:13790–4. doi: 10.1073/pnas.1404623111 25201988

39. van den Berg SM, de Moor MHM, McGue M, Pettersson E, Terracciano A, Verweij KJH, et al. Harmonization of neuroticism and extraversion phenotypes across inventories and cohorts in the Genetics of Personality Consortium: an application of Item Response Theory. Behav Genet. 2014;44:295–313. doi: 10.1007/s10519-014-9654-x 24828478

40. Stringer S, Minică CC, Verweij KJH, Mbarek H, Bernard M, Derringer J, et al. Genome-wide association study of lifetime cannabis use based on a large meta-analytic sample of 32 330 subjects from the International Cannabis Consortium. Transl Psychiatry. 2016;6:e769. doi: 10.1038/tp.2016.36 27023175

41. Middeldorp CM, Hammerschlag AR, Ouwens KG, Groen-Blokhuis MM, St. Pourcain B, Greven CU, et al. A genome-wide association meta-analysis of attention-deficit/hyperactivity disorder symptoms in population-based pediatric cohorts. J Am Acad Child Adolesc Psychiatry. 2016;55:896–905.e6. doi: 10.1016/j.jaac.2016.05.025 27663945

42. Demontis D, Walters RK, Martin J, Mattheisen M, Als TD, Agerbo E, et al. Discovery of the first genome-wide significant risk loci for attention deficit/hyperactivity disorder. Nat Genet. 2019;51:63–75. doi: 10.1038/s41588-018-0269-7 30478444

43. Walters RK, Polimanti R, Johnson ECEO, McClintick JN, Adams MJ, Adkins AE, et al. Transancestral GWAS of alcohol dependence reveals common genetic underpinnings with psychiatric disorders. Nat Neurosci. 2018;21:1656–69. doi: 10.1038/s41593-018-0275-1 30482948

44. Otowa T, Hek K, Lee M, Byrne EM, Mirza SS, Nivard MG, et al. Meta-analysis of genome-wide association studies of anxiety disorders. Mol Psychiatry. 2016;21:1391–9. doi: 10.1038/mp.2015.197 26754954

45. Psychiatric GWAS Consortium Bipolar Disorder Working Group. Large-scale genome-wide association analysis of bipolar disorder identifies a new susceptibility locus near ODZ4. Nat Genet. 2011;43:977–83. doi: 10.1038/ng.943 21926972

46. Schizophrenia Working Group of the Psychiatric Genomics Consortium. Biological insights from 108 schizophrenia-associated genetic loci. Nature. 2014;511:421–7. doi: 10.1038/nature13595 25056061

47. Tobacco and Genetics Consortium. Genome-wide meta-analyses identify multiple loci associated with smoking behavior. Nat Genet. 2010;42:441–7. doi: 10.1038/ng.571 20418890

48. Schumann G, Liu C, O’Reilly P, Gao H, Song P, Xu B, et al. KLB is associated with alcohol drinking, and its gene product β-Klotho is necessary for FGF21 regulation of alcohol preference. Proc Natl Acad Sci U S A. 2016;113:14372–7. doi: 10.1073/pnas.1611243113 27911795

49. de Moor MHM, Costa PT, Terracciano A, Krueger RF, de Geus EJC, Toshiko T, et al. Meta-analysis of genome-wide association studies for personality. Mol Psychiatry. 2012;17:337–49. doi: 10.1038/mp.2010.128 21173776

50. van den Berg SM, de Moor MHM, Verweij KJH, Krueger RF, Luciano M, Arias Vasquez A, et al. Meta-analysis of genome-wide association studies for extraversion: findings from the Genetics of Personality Consortium. Behav Genet. 2016;46:170–82. doi: 10.1007/s10519-015-9735-5 26362575

51. Pappa I, St Pourcain B, Benke K, Cavadino A, Hakulinen C, Nivard MG, et al. A genome-wide approach to children’s aggressive behavior: the EAGLE consortium. Am J Med Genet B Neuropsychiatr Genet. 2016;171:562–72. doi: 10.1002/ajmg.b.32333 26087016

52. Tielbeek JJ, Johansson A, Polderman TJC, Rautiainen MR, Jansen P, Taylor M, et al. Genome-wide association studies of a broad spectrum of antisocial behavior. JAMA Psychiatry. 2017;74:1242–50. doi: 10.1001/jamapsychiatry.2017.3069 28979981

53. van der Valk RJP, Kreiner-Møller E, Kooijman MN, Guxens M, Stergiakouli E, Sääf A, et al. A novel common variant in DCST2 is associated with length in early life and height in adulthood. Hum Mol Genet. 2015;24:1155–68. doi: 10.1093/hmg/ddu510 25281659

54. Horikoshi M, Yaghootkar H, Mook-Kanamori DO, Sovio U, Taal HR, Hennig BJ, et al. New loci associated with birth weight identify genetic links between intrauterine growth and adult height and metabolism. Nat Genet. 2013;45:76–82. doi: 10.1038/ng.2477 23202124

55. Wood AR, Esko T, Yang J, Vedantam S, Pers TH, Gustafsson S, et al. Defining the role of common variation in the genomic and biological architecture of adult human height. Nat Genet. 2014;46:1173–86. doi: 10.1038/ng.3097 25282103

56. Berndt SI, Gustafsson S, Mägi R, Ganna A, Wheeler E, Feitosa MF, et al. Genome-wide meta-analysis identifies 11 new loci for anthropometric traits and provides insights into genetic architecture. Nat Genet. 2013;45:501–12. doi: 10.1038/ng.2606 23563607

57. Locke AE, Kahali B, Berndt SI, Justice AE, Pers TH, Day FR, et al. Genetic studies of body mass index yield new insights for obesity biology. Nature. 2015;518:197–206. doi: 10.1038/nature14177 25673413

58. Chatterjee N, Shi J, García-Closas M. Developing and evaluating polygenic risk prediction models for stratified disease prevention. Nat Rev Genet. 2016;17:392–406. doi: 10.1038/nrg.2016.27 27140283

59. Benjamini Y, Hochberg Y. Controlling the false discovery rate: a practical and powerful approach to multiple testing. J R Stat Soc. 1995;57:289–300.

60. Hemani G, Zheng J, Elsworth B, Wade KH, Haberland V, Baird D, et al. The MR-Base platform supports systematic causal inference across the human phenome. Elife. 2018;7:e34408. doi: 10.7554/eLife.34408 29846171

61. Hemani G, Bowden J, Davey Smith G. Evaluating the potential role of pleiotropy in Mendelian randomization studies. Hum Mol Genet. 2018;27:R195–208. doi: 10.1093/hmg/ddy163 29771313

62. Burgess S, Butterworth A, Thompson SG. Mendelian randomization analysis with multiple genetic variants using summarized data. Genet Epidemiol. 2013;37:658–65. doi: 10.1002/gepi.21758 24114802

63. Zhao Q, Wang J, Hemani G, Bowden J, Small DS. Statistical inference in two-sample summary-data Mendelian randomization using robust adjusted profile score. arXiv: 1801.09652. 2018 Jan 29.

64. Bowden J, Davey Smith G, Haycock PC, Burgess S. Consistent estimation in Mendelian randomization with some invalid instruments using a weighted median estimator. Genet Epidemiol. 2016;40:304–14. doi: 10.1002/gepi.21965 27061298

65. Bowden J, Smith GD, Burgess S. Mendelian randomization with invalid instruments: effect estimation and bias detection through Egger regression. Int J Epidemiol. 2015;44:512–25. doi: 10.1093/ije/dyv080 26050253

66. Bowden J, Del Greco MF, Minelli C, Davey Smith G, Sheehan NA, Thompson JR. Assessing the suitability of summary data for two-sample Mendelian randomization analyses using MR-Egger regression: the role of the I2 statistic. Int J Epidemiol. 2016;45:1961–74. doi: 10.1093/ije/dyw220 27616674

67. Hemani G, Tilling K, Davey Smith G. Orienting the causal relationship between imprecisely measured traits using GWAS summary data. PLoS Genet. 2017;13:e1007081. doi: 10.1371/journal.pgen.1007081 29149188

68. Angelakis I, Gooding P, Tarrier N, Panagioti M. Suicidality in obsessive compulsive disorder (OCD): a systematic review and meta-analysis. Clin Psychol Rev. 2015;39:1–15. doi: 10.1016/j.cpr.2015.03.002 25875222

69. James A, Lai FH, Dahl C. Attention deficit hyperactivity disorder and suicide: a review of possible associations. Acta Psychiatr Scand. 2004;110:408–15. doi: 10.1111/j.1600-0447.2004.00384.x 15521824

70. Khemiri L, Kuja-Halkola R, Larsson H, Jayaram-Lindström N. Genetic overlap between impulsivity and alcohol dependence: a large-scale national twin study. Psychol Med. 2016;46:1091–102. doi: 10.1017/S0033291715002652 26671289

71. Franklin JC, Ribeiro JD, Fox KR, Bentley KH, Kleiman EM, Huang X, et al. Risk factors for suicidal thoughts and behaviors: a meta-analysis of 50 years of research. Psychol Bull. 2017;143:187–232. doi: 10.1037/bul0000084 27841450

72. Liang SH-Y, Yang Y-H, Kuo T-Y, Liao Y-T, Lin T-C, Lee Y, et al. Suicide risk reduction in youths with attention-deficit/hyperactivity disorder prescribed methylphenidate: a Taiwan nationwide population-based cohort study. Res Dev Disabil. 2018;72:96–105. doi: 10.1016/j.ridd.2017.10.023 29121517

73. Braun C, Bschor T, Franklin J, Baethge C. Suicides and suicide attempts during long-term treatment with antidepressants: a meta-analysis of 29 placebo-controlled studies including 6,934 patients with major depressive disorder. Psychother Psychosom. 2016;85:171–9. doi: 10.1159/000442293 27043848

74. Fry A, Littlejohns TJ, Sudlow C, Doherty N, Adamska L, Sprosen T, et al. Comparison of sociodemographic and health-related characteristics of UK Biobank participants with those of the general population. Am J Epidemiol. 2017;186:1026–34. doi: 10.1093/aje/kwx246 28641372

75. Sullivan PF, Neale MC, Kendler KS. Genetic epidemiology of major depression: review and meta-analysis. Am J Psychiatry. 2000;157:1552–62. doi: 10.1176/appi.ajp.157.10.1552 11007705

76. Voracek M, Loibl LM. Genetics of suicide: a systematic review of twin studies. Wien Klin Wochenschr. 2007;119:463–75. doi: 10.1007/s00508-007-0823-2 17721766


Článek vyšel v časopise

PLOS Medicine


2020 Číslo 6
Nejčtenější tento týden
Nejčtenější v tomto čísle
Kurzy

Zvyšte si kvalifikaci online z pohodlí domova

KOST
Koncepce osteologické péče pro gynekology a praktické lékaře
nový kurz
Autoři: MUDr. František Šenk

Sekvenční léčba schizofrenie
Autoři: MUDr. Jana Hořínková

Hypertenze a hypercholesterolémie – synergický efekt léčby
Autoři: prof. MUDr. Hana Rosolová, DrSc.

Svět praktické medicíny 5/2023 (znalostní test z časopisu)

Imunopatologie? … a co my s tím???
Autoři: doc. MUDr. Helena Lahoda Brodská, Ph.D.

Všechny kurzy
Kurzy Podcasty Doporučená témata Časopisy
Přihlášení
Zapomenuté heslo

Zadejte e-mailovou adresu, se kterou jste vytvářel(a) účet, budou Vám na ni zaslány informace k nastavení nového hesla.

Přihlášení

Nemáte účet?  Registrujte se

#ADS_BOTTOM_SCRIPTS#