We investigate whether, in Swedish national registers, social and psychiatric outcomes for six major psychiatric and substance disorders – drug use disorder (DUD), alcohol use disorder (AUD), major depression (MD), bipolar disorder (BD), anxiety disorder (AD), and schizophrenia (SZ) – reflect the primary genetic risk for each disorder and the level of genetic heterogeneity.

We utilize Genetic Risk Ratios – defined as the ratio of the genetic risk for secondary disorders to the genetic risk for the primary disorder – derived from Family Genetic Risk Scores. Poor social outcome was defined by a common factor of four variables: receipt of social welfare, sick leave, early retirement pension, and residence in a socially deprived area. Psychiatric outcome was defined as days of inpatient psychiatric hospitalization.

With poorer social outcomes, the primary genetic risks rose robustly for all disorders except SZ, as did the secondary genetic risks for DUD, AUD, and attention-deficit hyperactivity disorder. With poorer psychiatric outcomes, available only for BD and SZ, the primary genetic risks increased sharply. Overall, MD, AD, and BD became substantially more genetically heterogenous as their social outcomes became poorer, while for AUD, DUD, and SZ, the increase in heterogeneity was more modest. By contrast, with poorer psychiatric outcome, genetic risks for SZ became substantially more genetically homogeneous, with a similar but less robust trend seen for BD.

Despite important differences between our primary disorders, social and psychiatric outcomes are often robust indices of genetic risk and can reflect the levels of genetic heterogeneity.

Functional disorders (FDs) are associated with internalizing disorders (IDs). Studies investigating the nature of these associations over time are limited. We tested the direction of causation between measures of IDs (major depressive disorder [MDD], generalized anxiety disorder [GAD]) and FDs (fibromyalgia [FM] and myalgic encephalomyelitis/chronic fatigue syndrome [ME/CFS]) measured across two waves of longitudinal data (N = 108,034 and N = 73,590).

The Lifelines Cohort Study is a large prospective population-based cohort study in the northeast of the Netherlands. We tested competing causal models for the longitudinal association between IDs and FDs and, to follow-up results from the model with all IDs and FDs, tested the direction of causation between MDD and FM.

FDs were more stable over time than IDs. Initial model comparisons support a bidirectional relationship between most IDs and FDs. Follow-up analyses support a unidirectional model where FM predicts MDD over time (β = 0.14, 95% confidence interval = [0.11, 0.18]), but not vice versa.

The cross-time associations between ME/CFS, MDD, and GAD appear bidirectional (causal in both directions). Our results are consistent with, but not demonstrative of, a causal relationship from FM to MDD. The consequences of specific FDs vary, underscoring the value of studying these conditions as distinct constructs.

Functional disorders (FDs) are characterized by persistent somatic symptoms and are highly comorbid with internalizing disorders (IDs). To provide much-needed insight into FD etiology, we evaluated FD and ID familial coaggregation and shared familiality.

Lifelines is a three-generation cohort study, which assessed three FDs (myalgic encephalomyelitis/chronic fatigue syndrome [ME/CFS], irritable bowel syndrome [IBS], and fibromyalgia [FM]) and six IDs (major depressive disorder [MDD], dysthymia [DYS], generalized anxiety disorder [GAD], agoraphobia [AGPH], social phobia [SPH], and panic disorder [PD]) according to diagnostic criteria. Based on 153,803 individuals, including 90,397 with a first-degree relative in Lifelines, we calculated recurrence risk ratios (λRs) and tetrachoric correlations to evaluate familial aggregation and coaggregation of these disorders in first-degree relatives. We then estimated their familiality and familial correlations.

Familial aggregation was observed across disorders, with λR ranging from 1.45 to 2.23 within disorders and from 1.17 to 1.94 across disorders. Familiality estimates ranged from 22% (95% confidence interval [CI]: 16–29) for IBS to 42% (95% CI: 33–50) for ME/CFS. Familial correlations ranged from +0.37 (95% CI: 0.24–0.51) between FM and AGPH to +0.97 (95% CI: 0.80–1) between ME/CFS and FM. The highest familial correlation between an ID and FD was +0.83 (95% CI: 0.66–0.99) for MDD and ME/CFS.

There is a clear familial component to FDs, which is partially shared with IDs. This suggests that IDs and FDs share both genetic and family-environmental risk factors. Of the FDs, ME/CFS is most closely related to IDs.

Functional Somatic Disorders (FSD) and Internalizing Psychiatric Disorders (IPD) are frequently comorbid and likely share familial/genetic risk factors.

We performed a Common Factor Multivariate Analysis of 2 FSDs, Fibromyalgia (FM) and Irritable Bowel Syndrome (IBS), and two IPDs, Major Depression (MD) and Anxiety Disorders (AD), in five kinds of Swedish female–female relative pairs: monozygotic (n = 8,052) dizygotic (n = 7216), full siblings (n = 712,762), half-siblings reared together (n = 23,623), and half-siblings reared apart (n = 53,873). Model fitting was by full information maximum likelihood using OpenMx.

The best-fit model included genetic, shared environmental, and unique environmental factors. The common factor, ~50% heritable with a small shared environmental effect, loaded more strongly on the two IPDs (~0.80) than the 2 FSDs (0.40). Disorder-specific genetic effects were larger for the 2 FSDs (~0.30) than the 2 IPDs (~0.03). Estimated genetic correlations were high for MD and AD (+0.91), moderate between IBS and IPDs (+0.62), and intermediate between FM and MD (+0.54), FM and AD (+0.28), and FM and IBS (+0.38). Shared environmental influences on all disorders were present but small.

In women, FSDs and IPDs shared a moderate proportion of their genetic risk factors, greater for IBS than for FM. However, the genetic sharing between IBS and FM was less than between MD and AD, suggesting that FSDs do not form a highly genetically coherent group of disorders. The shared environment made a modest contribution to the familial aggregation of FSDs and IPDs.

We utilize a novel contrastive genetic-epidemiological method, the Maternal Half-Sibling Families with Discordant Fathers (MHSFDF) design, to examine cross-generational genetic transmission of posttraumatic stress disorder (PTSD) and related internalizing major depression (MD), and externalizing disorders: alcohol use disorder (AUD) and drug use disorder (DUD).

Using Swedish national registries, we identified 72,467 maternal half-sibling pairs reared together whose biological fathers were discordant for the diagnoses of PTSD, MD, AUD, and DUD. Offspring selected had to have less than 1 year of contact with their affected fathers. We examined the differences in outcome for within- and cross-disorder risk of diagnosis in the half-siblings with an affected versus unaffected father.

Paternal PTSD increased the risk of PTSD (HR: 1.43, 95% CI: 1.05–1.96) and MD (HR: 1.55, CI: 1.28–1.88) in offspring. It did not, however, elevate the risk of externalizing disorders (AUD or DUD). Offspring of fathers with AUD, DUD, or MD had increased risk of PTSD, suggesting sharing of vertically transmitted genetic risk between these disorders. No sex effects were found for any studied diagnosis.

This study is the first to show cross-generation genetic transmission for PTSD using the MHSFDF design. The pattern of cross-disorder genetic risk broadly supported an internalizing versus externalizing disorder split.

Stressful life events (SLEs) increase the risk for subsequent major depression (MD) episodes. Most prior studies of SLEs utilized questionnaires or interview-based assessments. We sought to develop and evaluate an environmental risk score (ERS) for MD from multiple classes of SLEs obtained from national Swedish registries.

We assessed, in the entire adult population of Sweden (n = 7,105,712), the occurrence of 52 categories of SLEs derived from registry information for the 6 months prior to January 9, 2010 and the risk for MD registration over the subsequent 6 months. Weights for our ERS were obtained from a random half of our sample and ERS and its relationship to MD risk was evaluated in the second half.

The ERS was robustly related to risk for subsequent MD episodes. Women were more sensitive to the depressogenic effect of the ERS than men. Those with prior episodes of MD had larger absolute increases in MD risk from our ERS than those without previous episodes. Genetic risk for MD was associated with a greater sensitivity to the depressogenic effects of the ERS. A co-sibling control analysis suggested that most, but not all, of the association between ERS and subsequent risk for MD was causal.

Valid measures of SLEs that predispose to risk for MD can be assessed from high-quality registry data. While not all event categories (e.g. interpersonal or romantic difficulties) can be assessed, this method avoids problems with accurate dating and recall bias and can be performed in very large samples.

We seek to clarify how changes in the prevalence of drug use disorder (DUD) in Sweden in the 1950–1990 birth cohort impact the aggregation and co-aggregation in siblings of DUD and alcohol use disorder (AUD).

We examined risk for DUD and AUD in siblings of 102,624 DUD cases and matched control probands and 123,837 AUD case and matched control probands identified using Swedish registries. Flexible parametric survival models assessed the difference in disorder risk in siblings of case versus control probands.

Over birthyears 1950–1990, rates of DUD increased substantially in the Swedish population. In siblings of DUD cases versus controls, the risk for DUD increased dramatically starting in birthyear 1965 while their risk for AUD fell moderately. A similar, but less pronounced pattern, was seen in the siblings of AUD versus control probands. These differences were much larger in male than in female siblings.

The factors that drove upward population rates of DUD in Sweden (e.g. increased availability, reduced stigma) produced much stronger effects in high-risk subjects (siblings of DUD and AUD probands) than in normal risk groups (siblings of controls), thereby increasing familial aggregation of DUD. However, parallel declines in AUD rates in high-risk versus normal-risk siblings were observed, likely due to ‘competitive effects’ reducing coaggregation of DUD and AUD. Results of genetic studies of substance use disorders can be substantially impacted by changes in availability and stigma of psychoactive substance use and indirectly by ‘competition’ as predicted by behavioral economic models, between abusable substances.

Drug use Disorder (DUD), the risk for which is substantially influenced by both genetic and social factors, is geographically concentrated in high-risk regions. An important step toward understanding this pattern is to examine geographical distributions of the genetic liability to DUD and a key demographic risk factor – social deprivation.

We calculated the mean family genetic risk score (FGRS) for DUD ((FGRSDUD) and social deprivation for each of the 5983 areas Demographic Statistical Areas (DeSO) for all of Sweden and used geospatial techniques to analyze and map these factors.

Using 2018 data, substantial spatial heterogeneity was seen in the distribution of the genetic risk for DUD in Sweden as a whole and in its three major urban centers which was confirmed by hot-spot analyses. Across DeSOs, FGRSDUD and s.d. levels were substantially but imperfectly correlated (r = + 0.63), with more scattering at higher FGRSDUD and s.d. scores. Joint mapping across DeSOs for FGRSDUD and s.d. revealed a diversity of patterns across Sweden. The stability of the distributions of FGRSDUD and s.d. in DeSOs within Sweden over the years 2012–2018 was quite high.

The geographical distribution of the genetic risk to DUD is quite variable in Sweden. DeSO levels of s.d. and FRGSDUD were substantially correlated but also disassociated in a number of regions. The observed patterns were largely consistent with known trends in the human geography of Sweden. This effort lays the groundwork for further studies of the sources of geographic variation in rates of DUD.

Robust clinical indices of etiologic heterogeneity for psychiatric disorders are rare. We investigate whether age at onset (AAO) reflects genetic heterogeneity, utilizing Genetic Risk Ratios (GRR) derived from Family Genetic Risk Scores (FGRS).

We examined, in individuals born in Sweden 1940–2003, whether AAO for five primary disorders -- drug use disorder (DUD), alcohol use disorder (AUD), major depression (MD), bipolar disorder (BD), and schizophrenia (SZ)-- was associated with varying levels of GRRs with a range of informative secondary disorders and traits.

Our disorders displayed a varying pattern of change of GRRs with increasing AAO. At one end was SZ, where all GRRs rose with increasing AAO meaning that SZ became increasing genetically heterogeneous with later AAO. The most balanced disorder was AUD where, with increasing AAO, GRRs rose for AD, BD, and MD and declined for DUD, CB, and ADHD. That is, at young AAO, AUD had high levels of genetic risk for other externalizing disorders while at older AAO, high genetic risk for internalizing disorders were more prominent. MD was at the continuum's other end where all GRRs, except for AD, decreased with higher AAO, meaning that MD became increasingly genetically homogeneous with later AAO.

Genetic heterogeneity was robustly associated with AAO across our five primary disorders with substantial inter-disorder differences in the observed patterns. In particular, young AAO was associated with maximal genetic homogeneity for SZ and DUD while older AAO had greater genetic homogeneity for MD with AUD falling in between.

Shared genetic risk between schizophrenia (SCZ) and bipolar disorder (BD) is well-established, yet the extent to which they share environmental risk factors remains unclear. We compare the associations between environmental exposures during childhood/prior to disorder onset with the risk of developing SCZ and BD.

We conducted a Swedish register-based nested case–control study using 4184 SCZ cases and 18 681 BD cases diagnosed 1988–2013. Cases were matched to five controls by birth year, birth region, and sex. Conditional logistic regression was used to estimate incidence rate ratios (IRR) for SCZ and BD for each exposure (severe childhood infections, adverse childhood experiences (ACEs), substance use disorders (SUDs), urban birth/longest residence).

All SUD types were associated with very high risk (IRR 4.9–25.5), and all forms of ACEs with higher risk (IRR 1.5–4.3) for both disorders. In the mutually adjusted models, ACEs demonstrated slightly higher risk for BD (SCZ IRR 1.30, 1.19-1.42; BD IRR 1.49, 1.44–1.55), while for SUD, risk was higher for SCZ (SCZ IRR 9.43, 8.15–10.92; BD IRR 5.50, 5.15–5.88). Infections were associated with increased risk of BD (IRR 1.21, 1.17–1.26) but not SCZ. Urban birth and urban longest residence were associated with higher risk of SCZ (IRR 1.19, 1.03–1.37), while only the combination of urban birth and rural longest residence showed higher risk for BD (IRR 1.24, 1.13–1.35).

There were both shared and unique environmental risk factors: SUDs and ACEs were risk factors for both disorders, while infections were more strongly associated with BD and urbanicity with SCZ.

The concept of schizotypal personality disorder (SPD) emerged from observations of personality characteristics common in relatives of schizophrenic patients. While often studied in family designs, few studies and none with genetic measures, have examined SPD in epidemiological samples.

We studied individuals born in Sweden 1940–2000 with an ICD-10 diagnosis of SPD with no prior schizophrenia (SZ) diagnosis (n = 2292). Demographic features, patterns of comorbidity, and Family Genetic Risk Scores (FGRS) were assessed from multiple Swedish registries. Prediction of progression to SZ was assessed by Cox models.

SPD was rare, with a prevalence of 0.044%, and had high levels of comorbidity with autism spectrum disorder (ASD), OCD, ADHD, and major depression (MD), and increased rates of being single, unemployed and in receipt of welfare. Affected individuals had elevated levels of FGRS for SZ (+0.42), ASD (+0.30), MD (+0.29), and ADHD (+0.20). Compared to cases of schizophrenia, they had significantly lower rates of FGRSSZ, but significantly elevated rates of genetic risk for ASD, MD, and ADHD. Over a mean follow-up of 8.7 years, 14.6% of SPD cases received a first diagnosis of SZ, the risk for which was significantly increased by levels of FGRSSZ, male sex, young age at SPD diagnosis and an in-patient SPD diagnosis and significantly decreased by comorbidity with MD, ASD, and ADHD.

Our results not only support the designation of SPD as a schizophrenia spectrum disorder but also suggest potentially important etiologic links between SPD and ASD and, to a lesser extent, ADHD, OCD, and MD.

One potential cause of comorbidity is the direct causal effect of one disorder – A – on risk for subsequent onset of disorder B. Could genetic risk scores be utilized to test for such an effect? If disorder A causally impacts on risk for disorder B, then genetic risk for disorder A should be lower in cases of disorder A with v. without a prior onset of B.

In all individuals (n = 905 736) born in Sweden from 1980 to 1990, from six psychiatric and drug use disorders (major depression, anxiety disorders, alcohol use disorder, drug use disorder, bipolar disorder, and schizophrenia), we formed 14 pairs of disorders A and B. In these pairs, we compared, using Cox proportional hazards models, the predictive effect of the familial-genetic risk score (FGRS) for disorder B in those who had v. had not had a prior onset of disorder A.

In all pairs, the impact of the FGRS for disorder B was significantly stronger in cases without v. with a prior history of disorder A. These effects were similar across sex, stable across levels of FGRS and not likely due to clinician bias. In many of our disorder pairs, previous clinical studies suggest a mechanism for a causal effect of disorder A on B.

Our findings provide indirect evidence that the occurrence of one psychiatric or substance use disorder often has a causal effect on risk for subsequent disorders. This mechanism may substantially contribute to the widespread comorbidity among psychiatric conditions.

To determine whether genetic risk factors for major depression (MD) and alcohol use disorder (AUD) interact with a potent stressor – death of spouse, parent, and sibling – in predicting episodes of, respectively, MD and AUD.

MD and AUD registrations were assessed from national Swedish registries. In individuals born in Sweden 1960–1970, we identified 7586, 388 459, and 34 370 with the loss of, respectively, a spouse, parent, and sibling. We started following subjects at age 18 or the year 2002 with end of follow-up in 2018. We examined time to event – a registration for MD within 6 months or AUD within a year – on an additive scale, using the Nelson–Aalen estimator. Genetic risk was assessed by the Family Genetic Risk Score (FGRS).

In separate models controlling for the main effects of death of spouse, parent, and sibling, FGRS, and sex, significant interactions were seen in all analyses between genetic risk for MD and death of relative in prediction of subsequent MD registration. A similar pattern of results, albeit with weaker interaction effects, was seen for genetic risk for AUD and risk for AUD registration. Genetic risk for bipolar disorder (BD) and anxiety disorders (AD) also interacted with event exposure in predicting MD.

Genetic risk for both MD and AUD act in part by increasing the sensitivity of individuals to the pathogenic effects of environmental stressors. For prediction of MD, similar effects are also seen for genetic risk for AD and BD.

Prior research has reported an association between divorce and suicide attempt. We aimed to clarify this complex relationship, considering sex differences, temporal factors, and underlying etiologic pathways.

We used Swedish longitudinal national registry data for a cohort born 1960–1990 that was registered as married between 1978 and 2018 (N = 1 601 075). We used Cox proportional hazards models to estimate the association between divorce and suicide attempt. To assess whether observed associations were attributable to familial confounders or potentially causal in nature, we conducted co-relative analyses.

In the overall sample and in sex-stratified analyses, divorce was associated with increased risk of suicide attempt (adjusted hazard ratios [HRs] 1.66–1.77). Risk was highest in the year immediately following divorce (HRs 2.20–2.91) and declined thereafter, but remained elevated 5 or more years later (HRs 1.41–1.51). Divorcees from shorter marriages were at higher risk for suicide attempt than those from longer marriages (HRs 3.33–3.40 and 1.20–1.36, respectively). In general, HRs were higher for divorced females than for divorced males. Co-relative analyses suggested that familial confounders and a causal pathway contribute to the observed associations.

The association between divorce and risk of suicide attempt is complex, varying as a function of sex and time-related variables. Given evidence that the observed association is due in part to a causal pathway from divorce to suicide attempt, intervention or prevention efforts, such as behavioral therapy, could be most effective early in the divorce process, and in particular among females and those whose marriages were of short duration.

It is clinically important to predict the conversion of major depression (MD) to bipolar disorder (BD). Therefore, we sought to identify related conversion rates and risk factors.

This cohort study included the Swedish population born from 1941 onward. Data were collected from Swedish population-based registers. Potential risk factors, including family genetic risk scores (FGRS), which were calculated based on the phenotypes of relatives in the extended family and not molecular data, and demographic/clinical characteristics from these registers were retrieved. Those with first MD registrations from 2006 were followed up until 2018. The conversion rate to BD and related risk factors were analyzed using Cox proportional hazards models. Additional analyses were performed for late converters and with stratification by sex.

The cumulative incidence of conversion was 5.84% [95% confidence interval (95% CI) 5.72–5.96] for 13 years. In the multivariable analysis, the strongest risk factors for conversion were high FGRS of BD [hazard ratio (HR) = 2.73, 95% CI 2.43–3.08], inpatient treatment settings (HR = 2.64, 95% CI 2.44–2.84), and psychotic depression (HR = 2.58, 95% CI 2.14–3.11). For late converters, the first registration of MD during the teenage years was a stronger risk factor when compared with the baseline model. When the interactions between risk factors and sex were significant, stratification by sex revealed that they were more predictive in females.

Family history of BD, inpatient treatment, and psychotic symptoms were the strongest predictors of conversion from MD to BD.

Are genetic risk factors for current depressive symptoms good proxies for genetic risk factors for syndromal major depression (MD)?

In over 9000 twins from the population-based Virginia Adult Twin Study of Psychiatric and Substance Use Disorders, the occurrence of all nine DSM symptomatic criteria for MD in the last year was assessed at personal interview and then grouped by their temporal co-occurrence. The DSM criteria which occurred outside (OUT) v. inside of (IN) MD episodes were then separated. We calculated tetrachoric correlations for OUT and IN depressive criteria in monozygotic (MZ) and dizygotic (DZ) pairs and fitted univariate and bivariate ACE twin models using OpenMx.

The mean twin correlations (±95% CIs) for IN depressive criteria were substantially higher than for OUT depressive criteria in both MZ [+0.35 (0.32–0.38) v. 0.20 (0.17–0.24)] and DZ pairs [0.20 (0.17–0.24) v. 0.10 (0.04–0.16]. The mean IN–OUT cross-correlation in MZ and DZ pairs was modest [+0.15 (0.07–0.24) and +0.07 (0.03–0.12)]. The mean heritability estimates for the nine In v. Out depressive criteria was 0.31 (0.22–0.41) and 0.15 (0.08–0.21), in MZ and DZ pairs, respectively. The mean genetic correlation between the nine IN and OUT depressive criteria was +0.07 (−0.07 to 0.21).

Depressive criteria occurring outside depressive episodes are less heritable than those occurring within. These two ways criteria can manifest are not closely genetically related. Current depressive symptoms – most of which are occurring outside of depressive episodes – are not, for genetic studies, good proxies for MD.