Hostname: page-component-669899f699-7tmb6 Total loading time: 0 Render date: 2025-05-05T09:15:28.771Z Has data issue: false hasContentIssue false
  • English
  • Français

Familial coaggregation and shared familiality of functional and internalizing disorders in the Lifelines cohort

Published online by Cambridge University Press:  02 May 2025

Martje Bos Open the ORCID record for Martje Bos [Opens in a new window] ,
Rei Monden Open the ORCID record for Rei Monden [Opens in a new window] ,
Naomi R. Wray Open the ORCID record for Naomi R. Wray [Opens in a new window] ,
Yiling Zhou ,
Kenneth S. Kendler Open the ORCID record for Kenneth S. Kendler [Opens in a new window] ,
Judith G. M. Rosmalen Open the ORCID record for Judith G. M. Rosmalen [Opens in a new window] ,
Hanna M. van Loo Open the ORCID record for Hanna M. van Loo [Opens in a new window]  and
Harold Snieder Open the ORCID record for Harold Snieder [Opens in a new window]
Martje Bos*
Affiliation:
Department of Psychiatry, University of Groningen, University Medical Center Groningen, Groningen, The Netherlands
Rei Monden
Affiliation:
Department of Psychiatry, University of Groningen, University Medical Center Groningen, Groningen, The Netherlands Informatics and Data Science Program, Graduate School of Advanced Science and Engineering, Hiroshima University, Hiroshima, Japan
Naomi R. Wray
Affiliation:
Institute for Molecular Bioscience, The University of Queensland, Brisbane, QLD, Australia Department of Psychiatry and Big Data Institute, University of Oxford, Oxford, UK
Yiling Zhou
Affiliation:
Department of Epidemiology, University of Groningen, University Medical Center Groningen, Groningen, The Netherlands
Kenneth S. Kendler
Affiliation:
Virginia Institute for Psychiatric and Behavioral Genetics, Virginia Commonwealth University, Richmond, VA, USA Department of Psychiatry, Virginia Commonwealth University, Richmond, VA, USA
Judith G. M. Rosmalen
Affiliation:
Department of Psychiatry, University of Groningen, University Medical Center Groningen, Groningen, The Netherlands Department of Internal Medicine, University of Groningen, University Medical Center Groningen, Groningen, The Netherlands
Hanna M. van Loo
Affiliation:
Department of Psychiatry, University of Groningen, University Medical Center Groningen, Groningen, The Netherlands
Harold Snieder
Affiliation:
Department of Epidemiology, University of Groningen, University Medical Center Groningen, Groningen, The Netherlands
*
Corresponding author: M. Bos; Email: [email protected]
Rights & Permissions [Opens in a new window]

Abstract

Background

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.

Methods

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.

Results

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.

Conclusions

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.

Keywords

familial transmissionfunctional somatic syndromesgenetic correlationsheritabilitymood and anxiety disorders
Type
Original Article
Information
Psychological Medicine , Volume 55 , 2025 , e126
Check for updates
Creative Commons
Creative Common License - CCCreative Common License - BY
This is an Open Access article, distributed under the terms of the Creative Commons Attribution licence (http://creativecommons.org/licenses/by/4.0), which permits unrestricted re-use, distribution and reproduction, provided the original article is properly cited.
Copyright
© The Author(s), 2025. Published by Cambridge University Press

Introduction

Functional disorders (FDs) are characterized by persistent somatic symptoms of unknown origin. In the absence of reproducibly observable pathophysiological processes, these disorders are diagnosed solely by symptoms. FDs are common (Haller, Cramer, Lauche, & Dobos, Reference Haller, Cramer, Lauche and Dobos2015; Nimnuan, Hotopf, & Wessely, Reference Nimnuan, Hotopf and Wessely2001), costly (Konnopka et al., Reference Konnopka, Schaefert, Heinrich, Kaufmann, Luppa, Herzog and König2012) and disabling (Joustra, Janssens, Bültmann, & Rosmalen, Reference Joustra, Janssens, Bültmann and Rosmalen2015). The three best-known FDs are myalgic encephalomyelitis/chronic fatigue syndrome (ME/CFS), irritable bowel syndrome (IBS), and fibromyalgia (FM), with estimated prevalences of 1.5%, 9.1%, and 1.8%, respectively (Heidari, Afshari, & Moosazadeh, Reference Heidari, Afshari and Moosazadeh2017; Lim et al., Reference Lim, Ahn, Jang, Lee, Lee and Son2020; Rometsch et al., Reference Rometsch, Mansueto, Maas Genannt Bermpohl, Martin and Cosci2024). Understanding of FD etiology is still limited, but both biological and psychosocial factors have been associated with these disorders (Kleinstäuber et al., Reference Kleinstäuber, Schröder, Daehler, Pallesen, Rask, Sanyer and Rosmalen2023).

There is evidence for a genetic component to FDs, with twin and genetic studies finding moderate levels of heritability for ME/CFS (Buchwald et al., Reference Buchwald, Herrell, Ashton, Belcourt, Schmaling, Sullivan and Goldberg2001), IBS (Svedberg, Johansson, Wallander, & Pedersen, Reference Svedberg, Johansson, Wallander and Pedersen2008), and FM (Dutta et al., Reference Dutta, Brummett, Moser, Fritsche, Tsodikov, Lee and Scott2020; Magnusson, Turkiewicz, Rydén, & Englund, Reference Magnusson, Turkiewicz, Rydén and Englund2024). Moreover, a family study observed heritable overlap between these FDs (Allen-Brady, Fyer, & Weissman, Reference Allen-Brady, Fyer and Weissman2023). To increase understanding of FD etiology, further clarification of the genetic liability to FDs is important. However, genome-wide association studies (GWASs) of FDs are still in the early stages and have limited power (Bonfiglio et al., Reference Bonfiglio, Henström, Nag, Hadizadeh, Zheng, Cenit and D’Amato2018; Hajdarevic et al., Reference Hajdarevic, Lande, Mehlsen, Rydland, Sosa, Strand and Viken2022; Moscati et al., Reference Moscati, Faucon, Arnaiz-Yépez, Lönn, Sundquist, Sundquist and Kendler2023). Studying genetic liability shared with related disorders can provide further insight into the genetics of FDs.

Due to their high comorbidity with FDs, internalizing disorders (IDs) are good candidates to study alongside FDs. For instance, individuals who meet diagnostic criteria for ME/CFS, IBS, or FM show higher rates of major depressive disorder (MDD) (odds ratios [ORs] = 3.87–12.62) and generalized anxiety disorder (GAD) (ORs = 3.19–9.81) than those who do not meet FD diagnostic criteria (Thomas et al., Reference Thomas, Gillespie, Kendler, Oldehinkel, Rosmalen and van Loo2024). Furthermore, a twin study suggests that FD-ID comorbidity is partly due to shared genetic factors (Kato, Sullivan, Evengård, & Pedersen, Reference Kato, Sullivan, Evengård and Pedersen2009). Similarly, a study using Swedish registry data found that individuals diagnosed with ME/CFS, IBS, or FM have an increased familial genetic risk for IDs (Kendler et al., Reference Kendler, Rosmalen, Ohlsson, Sundquist and Sundquist2023). Moreover, genetic loci associated with IBS have also been associated with depression and anxiety (Eijsbouts et al., Reference Eijsbouts, Zheng, Kennedy, Bonfiglio, Anderson, Moutsianas and Parkes2021; Tavares et al., Reference Tavares, Lopera-Maya, Bonfiglio, Zheng, Sinha, Marques and D’Amato2024; Tesfaye et al., Reference Tesfaye, Jaholkowski, Hindley, Shadrin, Rahman, Bahrami and Andreassen2023). Thus, studying both the genetics and the genetic relationships of IDs and FDs may advance the understanding of the etiology of both types of disorders.

Family studies provide a useful approach to study the vulnerability shared between disorders. Two such studies explored familial coaggregation of multiple FDs and IDs, both finding coaggregation between FM and MDD (Allen-Brady et al., Reference Allen-Brady, Fyer and Weissman2023; Hudson et al., Reference Hudson, Arnold, Keck, Auchenbach and Pope2004). One study also observed coaggregation across ME/CFS, IBS, FM, and panic disorder (PD) (Allen-Brady et al., Reference Allen-Brady, Fyer and Weissman2023). However, these studies either used medical records to define cases (Allen-Brady et al., Reference Allen-Brady, Fyer and Weissman2023) or recruited patients from medical centers (Hudson et al., Reference Hudson, Arnold, Keck, Auchenbach and Pope2004). Both methods risk oversampling severe cases and may be influenced by help-seeking behavior and diagnostic biases (Tattan et al., Reference Tattan, Ørnbøl, Wellnitz, Hanssen, Dantoft, Rosmalen and Petersen2024). The current study aims to assess familial aggregation, coaggregation, familiality, and familial correlations of IDs and FDs. Familial aggregation and coaggregation examine whether diseases cluster together in families. Familiality, a concept closely related to heritability, is the proportion of phenotypic variation attributable to familial effects, that is genetic and shared environmental effects combined. Familial correlation is a measure of how much familial effects on one trait overlap with familial effects on another trait. This study is performed in the large population-based Lifelines cohort (N = 153,803), which assessed FDs and IDs according to official diagnostic criteria.

Methods

Data

This study was conducted within the Lifelines cohort study. Lifelines is a multidisciplinary, prospective, population-based cohort study examining in a unique three-generation design the health and health-related behaviors of 167,729 persons living in the North of the Netherlands. It employs a broad range of investigative procedures to assess biomedical, sociodemographic, behavioral, physical, and psychological factors that contribute to health and disease in the general population, with a special focus on multimorbidity and complex genetics. Since 2006, three assessment waves have been completed (Sijtsma et al., Reference Sijtsma, Rienks, van der Harst, Navis, Rosmalen and Dotinga2022) in 2006–2013 (wave 1), 2014–2017 (wave 2), and 2019–2023 (wave 3). The current study used complete data from waves 1 and 2 as well as all data released from wave 3 up to 1 March 2024. The 153,803 adult participants with measurements on IDs or FDs during any of the three assessment waves were included in this study (Supplementary Figure 1).

The Lifelines cohort study followed the guidelines of the Declaration of Helsinki, and all procedures involving human subjects were approved by the Medical Ethical Committee of the University Medical Center Groningen. Furthermore, written informed consent was obtained from all participants.

Measurements

Internalizing disorders

In all three waves, current MDD, dysthymia (DYS), GAD, agoraphobia (AGPH), social phobia (SPH), and PD were evaluated using the Mini-International Neuropsychiatric Interview (MINI) (Sheehan et al., Reference Sheehan, Lecrubier, Sheehan, Amorim, Janavs, Weiller and Dunbar1998), which assesses these disorders per DSM-IV criteria (American Psychiatric Association, 1998). In wave 1, the MINI was administered as a face-to-face interview by a trained research nurse at a Lifelines research facility. In wave 2, participants completed a digital MINI questionnaire at the facility. In wave 3 used a digital MINI questionnaire emailed to participants for completion at home. The presence of MDD, DYS, and GAD required symptoms in the past two weeks, two years, and six months, respectively. The presence of AGPH, SPH, and PD required symptoms in the past month (van Loo et al., Reference van Loo, Beijers, Wieling, de Jong, Schoevers and Kendler2023), conforming to DSM-IV-TR duration criteria (American Psychiatric Association, 2000). Participants meeting the diagnostic criteria during any wave were labeled as cases; those who never met the criteria were labeled as controls. Thus, diagnostic status was determined based on assessments at up to three specific time points and does not include lifetime diagnoses.

Functional disorders

ME/CFS disease status was assessed according to the 1994 Centers for Disease Control and Prevention (CDC) diagnostic criteria (Fukuda et al., Reference Fukuda, Straus, Hickie, Sharpe, Dobbins and Komaroff1994). IBS status was assessed using ROME III criteria (Drossman, Reference Drossman2006), but the criteria related to symptom occurrence were adjusted to align with ROME IV criteria (Drossman, Reference Drossman2016). Thus, IBS was present if abdominal discomfort occurred more than once weekly, rather than three days a month, for over six months, in addition to two additional symptoms. FM was assessed according to the 2010 American College of Rheumatology (ACR) criteria (Wolfe et al., Reference Wolfe, Clauw, Fitzcharles, Goldenberg, Katz, Mease and Yunus2010). FDs were only assessed in Lifelines’ second and third waves. Participants meeting the diagnostic criteria during any wave were labeled as cases; those who did not were labeled as controls. Thus, diagnostic status was determined based on assessments at up to two specific time points and does not include lifetime diagnoses.

Statistical analyses

Familial aggregation and coaggregation

Recurrence risk ratios (λRs) were calculated to assess familial aggregation and coaggregation (Risch, Reference Risch1990). λR is the ratio of disease prevalence in relatives of affected participants to the prevalence in the general population, that is the Lifelines population. We estimated these prevalences using plug-in methods from the “marginaleffects” package in R (Arel-Bundock, Greifer, & Heiss, Reference Arel-Bundock, Greifer and Heiss2024). First, we fitted a logistic regression model in which disease status (case or control) was regressed on whether a person had an affected relative while adjusting for age, age2, sex, and number of relatives in the data. This model quantified how these factors are associated with disease risk. Second, we estimated disease prevalences by plugging-in the fitted model to our study population: once with their existing covariate values to estimate prevalence in the general population, and once with each participant’s relative affected status set to ‘true’ (while keeping existing values for age, sex, and family size) to estimate prevalence among those with affected relatives. The ratio of these two prevalence estimates gives the marginal λR, which represents how much having an affected relative increases disease risk, averaged across the total population’s distribution of age, sex, and family size. See supplementary methods for details. To ensure that we capture the full complexity of familial relationships without artificially separating comorbid cases, individuals with multiple diagnoses were included in the analysis for each of their conditions. To minimize reporting bias, proband reports of relatives’ disease status were not used (Milne et al., Reference Milne, Caspi, Crump, Poulton, Rutter, Sears and Moffitt2009; Saito et al., Reference Saito, Zimmerman, Harmsen, De Andrade, Locke, Petersen and Talley2008). To account for correlated observations due to familial clustering, 95% confidence intervals (CIs) were estimated using a robust clustered sandwich method (Zeileis, Köll, & Graham, Reference Zeileis, Köll and Graham2020). In addition to calculating the λR for first-degree relatives in general, λR was also calculated for siblings, parents, and offspring separately and for second-degree relatives. Furthermore, λR was calculated for cohabiting spouses as a measure of spousal resemblance, which can arise as the result of assortative mating and/or shared environmental influences.

Unlike λR, tetrachoric correlations are relatively insensitive to differences in the prevalence of the disorders involved in the calculations (Babchishin & Helmus, Reference Babchishin and Helmus2016; Cummings, Reference Cummings2009). Therefore, tetrachoric correlations were calculated for all parent-offspring pairs/trios and sibling pairs within families as an additional measures of familial aggregation and coaggregation. These correlations accounted for nonindependence within nuclear families within MPlus (Muthén & Muthén, Reference Muthén and Muthén1998).

Familiality and familial correlation

The terms heritability and genetic correlation imply that familial resemblance is solely due to genetics. Unlike twin studies, our family study did not disentangle genetic from shared environmental influences because an estimation of shared environment based on, for example, a comparison of first- and second-degree relatives would have limited power. To address this conflation and the resulting incomparability to heritability estimates from twin studies, we use the terms familiality and familial correlation instead of heritability and genetic correlation (Kendler & Neale, Reference Kendler and Neale2009).

To estimate the familiality of and familial correlations between FDs and IDs, we used the methods of Falconer (Reference Falconer1965) and Reich, James, and Morris (Reference Reich, James and Morris1972), as adapted by Wray and Gottesman (Reference Wray and Gottesman2012) (Baselmans, Yengo, van Rheenen, & Wray, Reference Baselmans, Yengo, van Rheenen and Wray2021). These methods are based on the liability threshold model, which assumes that a normally distributed liability underlies disease status. Individuals exceeding the critical liability threshold are affected. Although unobservable, this threshold can be determined using normal distribution theory, given the proportion of affected individuals in the population. The prevalence of disease in the general population and in thefirst- and second-degree relatives of affected individuals was used to estimate both the familiality of the disorders and the familial correlations between disorders. See supplementary methods for details.

Sensitivity analyses

We performed three sensitivity analyses. First, given the functional limitations associated with IDs and FDs (Buist-Bouwman et al., Reference Buist-Bouwman, De Graaf, Vollebergh, Alonso, Bruffaerts and Ormel2006; Joustra et al., Reference Joustra, Janssens, Bültmann and Rosmalen2015), we hypothesized a higher dropout rate between assessment waves for individuals with affected relatives. This selective dropout could lead to underestimating disease prevalence in relatives of affected individuals, consequently resulting in an underestimation of λR and familiality. We used logistic regression to compare participation in the second and third waves between participants with and without affected first-degree relatives in wave 1 or 2, adjusting for age, sex, and the number of participating relatives in wave 2. For IDs, the participant’s disease status in wave 1 was also considered, which was not possible for FDs as they were not assessed in wave 1. For wave 3, we additionally evaluated if individuals with affected relatives in any of the three waves were less likely to participate. Models were adjusted for age, sex, number of participating relatives in wave 3, and the individual’s disease status in previous waves.

Second, to ensure specificity to FDs, participants who met the diagnostic criteria for an FD and reported a medical condition with similar symptoms were excluded from the familial coaggregation analyses. For ME/CFS, participants with multiple sclerosis (MS), dementia, schizophrenia, or an eating disorder were excluded. For IBS, participants with ulcerative colitis, Crohn’s disease, or coeliac disease were excluded. For FM, participants with rheumatoid arthritis were excluded. Additionally, hepatitis, cancer, or heart failure were exclusion criteria for all FDs.

Third, we assessed the impact of differences in diagnostic criteria strictness for FDs. The diagnostic criteria for ME/CFS are stricter than those for IBS and FM; ME/CFS requires interference with daily tasks, while IBS and FM do not. Moreover, ME/CFS and IBS require a six-month duration, in contrast to three months for FM. Familial coaggregation is often stronger for more severe phenotypes (Steinhausen, Jakobsen, & Munk-Jørgensen, Reference Steinhausen, Jakobsen and Munk-Jørgensen2017; Wang, Snieder, & Hartman, Reference Wang, Snieder and Hartman2022). Therefore, the duration and interference criteria of FM and IBS were aligned with ME/CFS criteria, meaning that a symptom duration of six months and interference with daily life activities were required for all three FDs. This helped assess to which extent differences between FDs in λR and tetrachoric correlations were due to (arbitrary) diagnostic criteria or due to a difference in the type of symptoms.

Reporting and software

Results are considered significant if p < .01. Mplus version 8.2 (Muthén & Muthén, Reference Muthén and Muthén1998) was used to calculate tetrachoric correlations. R Version 4.2.1 (R Core Team, 2022) was used for all other analyses. R scripts are available at Open Science Framework (OSF), via doi: 10.17605/OSF.IO/7RCVT

Results

Descriptives

The diagnostic status for IDs or FDs could be determined for 153,803 Lifelines participants, based on up to three assessments between 2006 and 2023 for IDs and up to two assessments between 2014 and 2023 for FDs. The descriptive characteristics of these individuals are presented in Table 1. For 90,397 (58.8%) individuals, it was possible to determine disease status for at least one first-degree relative. The disease status of a second-degree relative could be assessed for 23,978 (15.6%) individuals. The data included 37,184 sibling pairs from 23,185 nuclear families, 16,455 parent-offspring pairs (a child with only one parent) from 11,916 nuclear families, and 17,046 parent-offspring trios (a child with both parents) from 11,357 nuclear families.

Table 1. Descriptive characteristics of dataset

Abbreviations: SD, standard deviation; MDD, major depressive disorder; DYS, dysthymia; GAD, generalized anxiety disorder; AGPH, agoraphobia; SPH, social phobia; PD, panic disorder; ME/CFS, myalgic encephalomyelitis/chronic fatigue syndrome; IBS, irritable bowel syndrome; FM, fibromyalgia.

a Age during the first adult assessment an individual participated in.

b Pairs were considered spouses if they lived together during the baseline assessment and indicated to be spouses or share a child together. This includes same-sex spouses.

c Cases aggregated across assessment waves 1, 2, and 3. Lifetime diagnoses were not established.

d Sample size of DYS is smaller than for other IDs as in some assessment waves DYS questions were skipped if MDD was present.

e Cases aggregated across assessment waves 2 and 3. Functional disorders were not assessed in first assessment wave. Lifetime diagnoses were not established.

f Sample size of FM is smaller than for other FDs as FM questions were spread across two questionnaires.

Familial aggregation and coaggregation

For all disorders except PD and DYS, having a first-degree relative affected by the disorder was significantly associated with a higher risk for all of the other disorders (Figure 1). Estimates of λR ranged from 1.17 (95% CI: 1.06–1.27) for IBS with a relative affected by AGPH to 2.23 for both PD and ME/CFS with a relative affected by the same disorder (95% CI: 1.26–3.20 for PD, 1.89–2.58 for ME/CFS). Familial coaggregation within IDs was similar to that within FDs. Among FDs, ME/CFS showed the strongest coaggregation with IDs, while IBS showed the weakest. Generally, familial coaggregation did not exhibit different patterns by type of first-degree relative (Supplementary Table 1).

Figure 1. Familial coaggregation of internalizing and functional disorders amongst first-degree relatives expressed in λR. Note: FDR, first-degree relative; MDD, major depressive disorder; DYS, dysthymia; GAD, generalized anxiety disorder; AGPH, agoraphobia; SPH, social phobia; PD, panic disorder; ME/CFS, myalgic encephalomyelitis/chronic fatigue syndrome; IBS, irritable bowel syndrome; FM, fibromyalgia; λR, recurrence risk ratios. λR adjusted for age, age2, sex, and number of relatives present in the data. Estimates in bold are significant at p < .01. For 95% confidence intervals, see Supplementary Table 1.

Unlike the λR estimates for all first-degree relatives combined, tetrachoric correlations calculated between sibling pairs and between parent-offspring pairs/trios did not consistently reach statistical significance (Figure 2). Patterns of familial coaggregation expressed in tetrachoric correlations did align with λR of sibling and parents-offspring pairs (Supplementary Tables 1–3). For instance, IBS was less strongly associated with IDs than ME/CFS and FM.

Figure 2. Tetrachoric correlations between internalizing and functional disorders, for (a) siblings and (b) parent-offspring pairs/trios. Note: MDD, major depressive disorder; DYS, dysthymia; GAD, generalized anxiety disorder; AGPH, agoraphobia; SPH, social phobia; PD, panic disorder; ME/CFS, myalgic encephalomyelitis/chronic fatigue syndrome; IBS, irritable bowel syndrome; FM, fibromyalgia; rtet, tetrachoric correlation. Estimates in bold are significant at p < .01. For 95% confidence intervals, see Supplementary Tables 2 & 3.

Within disorders, significant spousal resemblance was observed for MDD, DYS, GAD, AGPH, ME/CFS, and FM. Furthermore, all disorders showed spousal coaggregation with at least one other disorder (Figure 3). In some cases, λR for spouses was larger than λR for first-degree relatives. For instance, for MDD λR for first-degree relatives was 1.76 (95% CI: 1.54–1.98) while λR for spouses was 2.10 (95% CI: 1.63–2.56).

Figure 3. Familial coaggregation of internalizing and functional disorders amongst spouses, expressed in λR. Note: MDD, major depressive disorder; DYS, dysthymia; GAD, generalized anxiety disorder; AGPH, agoraphobia; SPH, social phobia; PD, panic disorder; ME/CFS, myalgic encephalomyelitis/chronic fatigue syndrome; IBS, irritable bowel syndrome; FM, fibromyalgia; λR, recurrence risk ratio. λR adjusted for age, age2, sex, and presence of spouse in the data. Estimates in bold are significant at p < .01. For 95% confidence intervals, see Supplementary Table 1.

Familiality and familial correlation

In general, familiality estimates were modest to moderate, ranging from 22% (95% CI: 16–29) for IBS to 42% (95% CI: 33–50) for ME/CFS (Figure 4a). Familiality estimates based on second-degree relatives were similar to those based on first-degree relatives, except for SPH and ME/CFS (Supplementary Table 4). This suggests limited influences from shared environment for most disorders as second-degree relatives typically share much less of their family and community environment than do first-degree relatives. Moderate to high familial correlations were observed across disorders, ranging from +0.37 (95% CI: 0.24–0.51) between AGPH and FM to +0.97 (95% CI 0.80–1) between ME/CFS and FM. The familial correlation between MDD and ME/CFS was +0.83 (95% CI: 0.66–0.99), which was the strongest familial correlation observed between an ID and FD (Figure 4b).

Figure 4 (a) Familiality and (b) familial correlation estimates of internalizing and functional disorders based on both first- and second-degree relatives. Note: MDD, major depressive disorder; DYS, dysthymia; GAD, generalized anxiety disorder; AGPH, agoraphobia; SPH, social phobia; PD, panic disorder; ME/CFS, myalgic encephalomyelitis/chronic fatigue syndrome; IBS, irritable bowel syndrome; FM, fibromyalgia; rf, familial correlation. See Supplementary Table 4 for familiality estimates for first- and second-degree relatives separately and Supplementary Table 5 for 95% confidence intervals of familial correlation estimates.

Sensitivity analysis

We assessed whether dropout between waves was higher among participants whose first-degree relative was affected by any of the disorders, as this could lead to an underestimation of disease prevalence in relatives of affected individuals. Individuals with a first-degree relative affected by MDD, GAD, AGPH, PD, or FM during the first or second assessment waves were less likely to participate in the MINI of wave 2. Those with a first-degree relative affected by PD were least likely to participate (OR: 0.86, 95% CI: 0.77–0.97). Individuals with a first-degree relative affected by GAD in wave 1 or 2 were also less likely to participate in the FD questionnaire of wave 2 (OR: 0.90, 95% CI: 0.86–0.94). For wave 3, individuals with a first-degree relative affected by MDD, GAD, AGPH, CFS, or FM in any of the assessment waves were less likely to participate in both the MINI and the FD questionnaire. Individuals whose first-degree relative was affected by FM were least likely to participate in the wave 3 MINI (OR: 0.86, 95% CI: 0.81–0.91) and the FD questionnaire (OR: 0.87, 95% CI: 0.81–0.92). See Supplementary Table 6 for all logistic regression results.

To ensure specificity to FDs, subjects with somatic disorders were excluded from familial aggregation and coaggregation analyses, reducing the prevalence of ME/CFS, IBS, and FM to 3.9%, 6.3%, and 6.6%, respectively. This exclusion had no to a minimal impact on familial coaggregation estimates (Supplementary Tables 7–9).

To account for differences in strictness of diagnostic criteria, IBS and FM duration and interference criteria were aligned with ME/CFS to include 6-6-month duration and interference with daily activities. This reduced IBS prevalence from 7.5% to 1.1% and FM prevalence from 8.2% to 4.5%. Severity-aligned IBS coaggregated more with IDs than the original criteria. Although ME/CFS still showed stronger coaggregation with IDs than severity-aligned FM, differences did attenuate, indicating that criteria strictness partly explains the familial aggregation differences (Supplementary Tables 10–12).

Discussion

This study assessed the familial aggregation and coaggregation, as well as the familiality of and familial correlations between FDs and IDs in a large general population sample. We observed significant aggregation within and between studied disorders. Furthermore, each disorder showed spousal coaggregation with at least one other disorder. Whether this results from assortative mating or shared environment is unknown. We also observed moderate familiality across all disorders, with familial correlations indicating that this familiality is moderately to highly shared between FDs and IDs.

Comparison with familial coaggregation studies

Our results show that individuals with a first-degree relative affected by an ID or an FD are at increased risk of developing the same or other studied disorders, consistent with previous studies on IDs (Hettema, Neale, & Kendler, Reference Hettema, Neale and Kendler2001; Sullivan, Neale, & Kendler, Reference Sullivan, Neale and Kendler2000), FDs (Albright et al., Reference Albright, Light, Light, Bateman and Cannon-Albright2011; Saito et al., Reference Saito, Petersen, Larson, Atkinson, Fridley, de Andrade and Talley2010) and across IDs and FDs (Allen-Brady et al., Reference Allen-Brady, Fyer and Weissman2023; Hudson et al., Reference Hudson, Arnold, Keck, Auchenbach and Pope2004). Our coaggregation patterns largely resemble those of a medical records-based study in a population registry database (Allen-Brady et al., Reference Allen-Brady, Fyer and Weissman2023). Minor discrepancies, like IBS exhibiting stronger familial aggregation than FM in the previous study (Allen-Brady et al., Reference Allen-Brady, Fyer and Weissman2023) might be due to the severity of IBS cases in medical records. Our sensitivity analysis supports this explanation, as severe IBS that interferes with daily activities showed stronger familial aggregation than noninterfering IBS. A study focusing on FM mirrored our findings with MDD but reported no significant coaggregation with GAD, SPH, PD, or IBS, possibly due to limited sample size and the small number of individuals affected by these disorders (Hudson et al., Reference Hudson, Arnold, Keck, Auchenbach and Pope2004). In contrast to λR in all first-degree relatives, tetrachoric correlations between sibling and parent-offspring pairs in the current study were not always significant, likely due to reduced statistical power and methodological differences. For instance, unlike λR, tetrachoric correlations assume an underlying normally distributed liability to disease. Although this is a common assumption in genetic epidemiology, it may not fully reflect the complex familiality of these disorders.

Comparison with genetic studies

We use the terms familiality and familial correlation rather than heritability and genetic correlation because our methods do not disentangle genetic from shared environmental effects. However, evidence suggests that for most disorders, familial aggregation results from genetic factors, with limited impact from shared environmental factors. For instance, in twin studies on IDs, IBS, and FM, models not including a shared environmental component are a better fit for the data than models that do include a shared environmental component (Hettema et al., Reference Hettema, Neale and Kendler2001; Magnusson et al., Reference Magnusson, Turkiewicz, Rydén and Englund2024; Markkula et al., Reference Markkula, Järvinen, Leino-Arjas, Koskenvuo, Kalso and Kaprio2009; Sullivan et al., Reference Sullivan, Neale and Kendler2000; Svedberg et al., Reference Svedberg, Johansson, Wallander and Pedersen2008). Additionally, except for SPH and ME/CFS, our familiality estimates from second-degree relatives were very close to those from first-degree relatives (Supplementary Table 4), despite second-degree relatives sharing less of their environment than first-degree relatives. This also suggests a limited role of the shared environment. Moreover, for well-studied disorders like MDD and GAD, our familiality estimates are lower than the heritability estimates from meta-analyses of twin studies (Hettema et al., Reference Hettema, Neale and Kendler2001; Sullivan et al., Reference Sullivan, Neale and Kendler2000), which would be unlikely if there were strong shared environmental effects. The familial correlation between MDD and GAD also aligns with the genetic correlation found in a major twin study (Kendler, Gardner, Gatz, & Pedersen, Reference Kendler, Gardner, Gatz and Pedersen2007). Therefore, we believe that the current study contributes to the growing body of evidence supporting a genetic component in FDs (Ablin & Buskila, Reference Ablin and Buskila2015; Dibble, McGrath, & Ponting, Reference Dibble, McGrath and Ponting2020; Saito, Reference Saito2011).

Our familiality estimates for ME/CFS and IBS are similar to heritability estimates from twin studies (Buchwald et al., Reference Buchwald, Herrell, Ashton, Belcourt, Schmaling, Sullivan and Goldberg2001; Svedberg et al., Reference Svedberg, Johansson, Wallander and Pedersen2008). Our familiality estimate for FM is higher than the heritability estimate (23%, 95% CI 14–32) of a twin study based on ICD-10 codes for myalgia (muscle pain) and FM (Magnusson et al., Reference Magnusson, Turkiewicz, Rydén and Englund2024). In contrast, our estimate is lower than the estimate (51%, 95% CI: 45–56) of a twin study that identified FM through questionnaire items related to symptoms like morning and evening stiffness, neck pain and stiffness, tender points, daytime fatigue, and numbness (Markkula et al., Reference Markkula, Järvinen, Leino-Arjas, Koskenvuo, Kalso and Kaprio2009).

Our findings also suggest that the liability to disease is shared between FDs and IDs, corroborating previous studies (Allen-Brady et al., Reference Allen-Brady, Fyer and Weissman2023; Kato et al., Reference Kato, Sullivan, Evengård and Pedersen2009; Kendler et al., Reference Kendler, Rosmalen, Ohlsson, Sundquist and Sundquist2023). A Swedish twin study identified two latent genetic variables that underlie the comorbidity of MDD, GAD, and FDs (Kato et al., Reference Kato, Sullivan, Evengård and Pedersen2009). MDD, GAD, and FDs loaded on one latent genetic variable, which is in line with the familial correlations we found between IDs and FDs. The second latent genetic variable was exclusively related to FDs (Kato et al., Reference Kato, Sullivan, Evengård and Pedersen2009), fitting the high familial correlations we found between FDs. Strong correlations between FDs could be because individuals meeting criteria for one FD often report symptoms included in the diagnostic criteria for other FDs (van der Meulen et al., Reference van der Meulen, Bos, Bakker, Gans and Rosmalen2024).

ME/CFS and FM exhibited a closer familial relationship than either had with IBS, while IDs were more linked to ME/CFS than to FM or IBS. These findings parallel the factor loadings of the Swedish twin study (Kato et al., Reference Kato, Sullivan, Evengård and Pedersen2009). However, they differ from a family-based Swedish registry study, where FM and ME/CFS were more related to IBS than to each other, and IDs were most associated with FM (Kendler et al., Reference Kendler, Rosmalen, Ohlsson, Sundquist and Sundquist2023). Similarly, a medical records-based study found that IBS, FM, and MDD showed more extensive overlap with other FDs or MDD in distant relatives than ME/CFS (Allen-Brady et al., Reference Allen-Brady, Fyer and Weissman2023). These differences could stem from diagnostic biases (Tattan et al., Reference Tattan, Ørnbøl, Wellnitz, Hanssen, Dantoft, Rosmalen and Petersen2024) and overrepresentation of severe cases in registry/records-based studies, unlike the self-reported symptom assessments in our general population study and the Swedish twin study (Kato et al., Reference Kato, Sullivan, Evengård and Pedersen2009).

Relevance and future research

The observed familial coaggregation and familial correlations between FDs and IDs suggest shared etiological mechanisms between these groups of disorders. The specific shared mechanisms remain unknown and most likely reflect multilevel, interlocking etiological pathways (Kendler, Reference Kendler2012). For instance, emotion regulation problems have been linked to IDs, ME/CFS, and FM (Bram, Gottschalk, & Leeds, Reference Bram, Gottschalk and Leeds2018; Picó-Pérez et al., Reference Picó-Pérez, Radua, Steward, Menchón and Soriano-Mas2017; Pinto et al., Reference Pinto, Geenen, Wager, Lumley, Häuser, Kosek and da Silva2023), while inflammation has been associated with FDs and MDD (Andrés-Rodríguez et al., Reference Andrés-Rodríguez, Borràs, Feliu-Soler, Pérez-Aranda, Angarita-Osorio, Moreno-Peral and Luciano2020; Burns et al., Reference Burns, Carroll, Mathe, Horvat, Foster, Walker and Keely2019; Harsanyi, Kupcova, Danisovic, & Klein, Reference Harsanyi, Kupcova, Danisovic and Klein2023; Strawbridge, Sartor, Scott, & Cleare, Reference Strawbridge, Sartor, Scott and Cleare2019). Genetic overlap previously found between FDs and somatic disorders further emphasizes the multifaceted nature of FD etiology (Kendler et al., Reference Kendler, Rosmalen, Ohlsson, Sundquist and Sundquist2023).

To uncover specific shared biological mechanisms, molecular genetic studies are needed. A GWAS revealed shared loci between IBS and anxiety disorders (Eijsbouts et al., Reference Eijsbouts, Zheng, Kennedy, Bonfiglio, Anderson, Moutsianas and Parkes2021). Genes mapped to shared loci between IBS and anxiety disorders regulate neural circuits and influence white matter microstructure (Eijsbouts et al., Reference Eijsbouts, Zheng, Kennedy, Bonfiglio, Anderson, Moutsianas and Parkes2021) and were enriched for pathways relevant to the nervous and immune system (Tesfaye et al., Reference Tesfaye, Jaholkowski, Hindley, Shadrin, Rahman, Bahrami and Andreassen2023). While our findings advance the understanding of FDs, extensive molecular genetic studies, including ME/CFS and FM, are needed to identify shared genomic loci, which could unveil the underlying mechanisms common to IDs and FDs.

Genetic correlations may also stem from a causal link on the symptom level. For instance, sleep dysfunction, which is a symptom of MDD, is a potential trigger for FM (Choy, Reference Choy2015). Mendelian randomization studies could provide additional understanding of the causal relationship between FDs and IDs (Davey Smith & Ebrahim, Reference Davey Smith and Ebrahim2004).

Strengths and limitations

This study’s strengths lie in its use of a large population-based cohort and a diagnostic algorithm for assessing self-reported symptoms against official diagnostic criteria. By avoiding registry data, medical center recruitment, or self-reported diagnoses, we minimized the impact of diagnostic biases and help-seeking behavior (Kendler, Reference Kendler1995; Talley & Spiller, Reference Talley and Spiller2002; Tattan et al., Reference Tattan, Ørnbøl, Wellnitz, Hanssen, Dantoft, Rosmalen and Petersen2024; Wolfe et al., Reference Wolfe, Schmukler, Jamal, Castrejon, Gibson, Srinivasan and Pincus2019), prevented oversampling of severe cases, and captured cases otherwise being overlooked (Warren & Clauw, Reference Warren and Clauw2012). Furthermore, the direct assessment of disease status in relatives minimizes potential biases associated with participants reporting on the health of their relatives, addressing concerns related to incomplete knowledge about family members’ health (Milne et al., Reference Milne, Caspi, Crump, Poulton, Rutter, Sears and Moffitt2009; Saito et al., Reference Saito, Zimmerman, Harmsen, De Andrade, Locke, Petersen and Talley2008).

The primary limitation of this study is the potential inaccuracy in estimating disease prevalence among individuals with an affected relative. This arises from several factors. Firstly, we lacked complete data on first-degree relatives for all Lifelines participants, with 41% of participants having no first-degree relatives in the data. These individuals could not be considered for the population of individuals with an affected relative. By including these individuals in the general population estimate, we implicitly assumed they were representative of the general population. However, our data suggest that this might not be the case as disease prevalence was lower among those with relatives in the dataset (Table 1), suggesting that individuals without relatives were less healthy. Recognizing that the number of family members in the data is associated with both our outcome (disease prevalence) and our exposure (having an affected relative), we included the number of relatives in the data as a covariate in our analyses, addressing some of the inaccuracies in prevalence estimation. However, a second limitation is participation bias, as individuals with relatives affected by MDD, GAD, AGPH, PD, IBS, or FM were less likely to participate in subsequent assessment waves. This likely resulted in an underestimation of disease prevalence among those with affected relatives. Moreover, in the absence of lifetime diagnoses, we relied on point-in-time diagnoses. These limitations have possibly led to an underestimation of λR and familiality estimates.

Conclusion

In conclusion, our study underscores familial and likely genetic overlap between FDs and IDs, suggesting potential shared etiological mechanisms. To clarify the specific nature of these shared mechanisms and explore potential causal relationships between FDs and IDs, additional studies incorporating genotype data are needed.

Supplementary material

The supplementary materials for this article can be found at https://doi.org/10.1017/S003329172500100X.

Data availability statement

Code for estimating λR, familiality, and familial correlations are available on OSF via doi: 10.17605/OSF.IO/7RCVT.

Acknowledgments

We thank the members of the Pathways of Functional and Internalizing Disorders (PATHFINDER) consortium for their contributions to this project. The Lifelines Biobank initiative has been made possible by a subsidy from the Dutch Ministry of Health, Welfare and Sport, the Dutch Ministry of Economic Affairs, the University Medical Center Groningen (UMCG the Netherlands), University of Groningen, and the Northern Provinces of the Netherlands. The authors wish to acknowledge the services of the Lifelines Cohort Study, the contributing researchers involved in data collection, and all the study participants.

Funding statement

This project was supported by NIMH grant R01MH125902. Hanna van Loo was supported by a Veni grant from the Talent Program of the Netherlands Organization of Scientific Research (NWOZonMW 09150161810021). Judith Rosmalen was supported by a Vici grant from the Talent Program of the Netherlands Organization of Scientific Research (NWO Vici Vi.C.191.021).

Competing interests

We declare no conflicts of interest.

Footnotes

Judith G. M. Rosmalen, Hanna M. van Loo and Harold Snieder supervised this work.

References

Ablin, J. N., & Buskila, D. (2015). Update on the genetics of the fibromyalgia syndrome. Best Practice & Research Clinical Rheumatology, 29(1), 20–28. https://doi.org/10.1016/j.berh.2015.04.018CrossRefGoogle ScholarPubMed
Albright, F., Light, K., Light, A., Bateman, L., & Cannon-Albright, L. A. (2011). Evidence for a heritable predisposition to chronic fatigue syndrome. BMC Neurology, 11, 62. https://doi.org/10.1186/1471-2377-11-62CrossRefGoogle ScholarPubMed
Allen-Brady, K., Fyer, A. J., & Weissman, M. (2023). The multi-generational familial aggregation of interstitial cystitis, other chronic nociplastic pain disorders, depression, and panic disorder. Psychological Medicine, 53(16), 7847–7856. https://doi.org/10.1017/S0033291723001885CrossRefGoogle ScholarPubMed
American Psychiatric Association (Ed.) (1998). Diagnostic and statistical manual of mental disorders: DSM-IV ; includes ICD-9-CM codes effective 1. (4th ed., 7th print). Washington, DC.: APAGoogle Scholar
American Psychiatric Association (Ed.) (2000). Diagnostic and statistical manual of mental disorders fourth edition text revision (DSM-IV-TR) (4th ed., 9. print). Washington, DC: APAGoogle Scholar
Andrés-Rodríguez, L., Borràs, X., Feliu-Soler, A., Pérez-Aranda, A., Angarita-Osorio, N., Moreno-Peral, P., … Luciano, J. V. (2020). Peripheral immune aberrations in fibromyalgia: A systematic review, meta-analysis and meta-regression. Brain, Behavior, and Immunity, 87, 881–889. https://doi.org/10.1016/j.bbi.2019.12.020CrossRefGoogle ScholarPubMed
Arel-Bundock, V., Greifer, N., & Heiss, A. (2024). How to Interpret Statistical Models Using marginaleffects for R and Python. Journal of Statistical Software, 111(9), 1–32. https://doi.org/10.18637/jss.v111.i09CrossRefGoogle Scholar
Babchishin, K. M., & Helmus, L.-M. (2016). The influence of base rates on correlations: An evaluation of proposed alternative effect sizes with real-world data. Behavior Research Methods, 48(3), 1021–1031. https://doi.org/10.3758/s13428-015-0627-7CrossRefGoogle ScholarPubMed
Baselmans, B. M. L., Yengo, L., van Rheenen, W., & Wray, N. R. (2021). Risk in relatives, heritability, SNP-based heritability, and genetic correlations in psychiatric disorders: A review. Biological Psychiatry, 89(1), 11–19. https://doi.org/10.1016/j.biopsych.2020.05.034CrossRefGoogle ScholarPubMed
Bonfiglio, F., Henström, M., Nag, A., Hadizadeh, F., Zheng, T., Cenit, M. C., … D’Amato, M. (2018). A GWAS meta-analysis from 5 population-based cohorts implicates ion channel genes in the pathogenesis of irritable bowel syndrome. Neurogastroenterology and Motility, 30(9), e13358. https://doi.org/10.1111/nmo.13358CrossRefGoogle ScholarPubMed
Bram, A. D., Gottschalk, K. A., & Leeds, W. M. (2018). Emotional regulation in women with chronic fatigue syndrome and depression: Internal representations and adaptive defenses. Journal of the American Psychoanalytic Association, 66(4), 701–741. https://doi.org/10.1177/0003065118798043CrossRefGoogle ScholarPubMed
Buchwald, D., Herrell, R., Ashton, S., Belcourt, M., Schmaling, K., Sullivan, P., … Goldberg, J. (2001). A twin study of chronic fatigue. Psychosomatic Medicine, 63(6), 936–943. https://doi.org/10.1097/00006842-200111000-00012CrossRefGoogle ScholarPubMed
Buist-Bouwman, M. A., De Graaf, R., Vollebergh, W. a. M, Alonso, J., Bruffaerts, R., Ormel, J., & ESEMeD/MHEDEA 2000 Investigators. (2006). Functional disability of mental disorders and comparison with physical disorders: A study among the general population of six European countries. Acta Psychiatrica Scandinavica, 113(6), 492–500. https://doi.org/10.1111/j.1600-0447.2005.00684.xCrossRefGoogle ScholarPubMed
Burns, G., Carroll, G., Mathe, A., Horvat, J., Foster, P., Walker, M. M., … Keely, S. (2019). Evidence for local and systemic immune activation in functional dyspepsia and the irritable bowel syndrome: A systematic review. Official Journal of the American College of Gastroenterology|ACG, 114(3), 429. https://doi.org/10.1038/s41395-018-0377-0CrossRefGoogle ScholarPubMed
Choy, E. H. S. (2015). The role of sleep in pain and fibromyalgia. Nature Reviews Rheumatology, 11(9), 513–520. https://doi.org/10.1038/nrrheum.2015.56CrossRefGoogle ScholarPubMed
Cummings, P. (2009). The relative merits of risk ratios and odds ratios. Archives of Pediatrics & Adolescent Medicine, 163(5), 438–445. https://doi.org/10.1001/archpediatrics.2009.31CrossRefGoogle ScholarPubMed
Davey Smith, G., & Ebrahim, S. (2004). Mendelian randomization: Prospects, potentials, and limitations. International Journal of Epidemiology, 33(1), 30–42. https://doi.org/10.1093/ije/dyh132CrossRefGoogle Scholar
Dibble, J. J., McGrath, S. J., & Ponting, C. P. (2020). Genetic risk factors of ME/CFS: A critical review. Human Molecular Genetics, 29(R1), R117–R124. https://doi.org/10.1093/hmg/ddaa169CrossRefGoogle ScholarPubMed
Drossman, D. A. (2006). The functional gastrointestinal disorders and the Rome III process. Gastroenterology, 130(5), 1377–1390. https://doi.org/10.1053/j.gastro.2006.03.008CrossRefGoogle ScholarPubMed
Drossman, D. A. (2016). Functional gastrointestinal disorders: History, pathophysiology, clinical features and Rome IV. Gastroenterology, S0016-5085(16) 00223-7. https://doi.org/10.1053/j.gastro.2016.02.032Google Scholar
Dutta, D., Brummett, C. M., Moser, S. E., Fritsche, L. G., Tsodikov, A., Lee, S., … Scott, L. J. (2020). Heritability of the fibromyalgia phenotype varies by Age. Arthritis & Rheumatology, 72(5), 815–823. https://doi.org/10.1002/art.41171CrossRefGoogle ScholarPubMed
Eijsbouts, C., Zheng, T., Kennedy, N. A., Bonfiglio, F., Anderson, C. A., Moutsianas, L., … Parkes, M. (2021). Genome-wide analysis of 53,400 people with irritable bowel syndrome highlights shared genetic pathways with mood and anxiety disorders. Nature Genetics, 53(11), 1543–1552. https://doi.org/10.1038/s41588-021-00950-8CrossRefGoogle ScholarPubMed
Falconer, D. S. (1965). The inheritance of liability to certain diseases, estimated from the incidence among relatives. Annals of Human Genetics, 29(1), 51–76. https://doi.org/10.1111/j.1469-1809.1965.tb00500.xCrossRefGoogle Scholar
Fukuda, K., Straus, S. E., Hickie, I., Sharpe, M. C., Dobbins, J. G., & Komaroff, A. (1994). The chronic fatigue syndrome: A comprehensive approach to its definition and study. International Chronic Fatigue Syndrome Study Group. Annals of Internal Medicine, 121(12), 953–959. https://doi.org/10.7326/0003-4819-121-12-199412150-00009Google ScholarPubMed
Hajdarevic, R., Lande, A., Mehlsen, J., Rydland, A., Sosa, D. D., Strand, E. B., … Viken, M. K. (2022). Genetic association study in myalgic encephalomyelitis/chronic fatigue syndrome (ME/CFS) identifies several potential risk loci. Brain, Behavior, and Immunity, 102, 362–369. https://doi.org/10.1016/j.bbi.2022.03.010CrossRefGoogle ScholarPubMed
Haller, H., Cramer, H., Lauche, R., & Dobos, G. (2015). Somatoform disorders and medically unexplained symptoms in primary care. Deutsches Arzteblatt International, 112(16), 279–287. https://doi.org/10.3238/arztebl.2015.0279Google ScholarPubMed
Harsanyi, S., Kupcova, I., Danisovic, L., & Klein, M. (2023). Selected biomarkers of depression: What are the effects of cytokines and inflammation? International Journal of Molecular Sciences, 24(1), 578. https://doi.org/10.3390/ijms24010578CrossRefGoogle Scholar
Heidari, F., Afshari, M., & Moosazadeh, M. (2017). Prevalence of fibromyalgia in general population and patients, a systematic review and meta-analysis. Rheumatology International, 37(9), 1527–1539. https://doi.org/10.1007/s00296-017-3725-2CrossRefGoogle ScholarPubMed
Hettema, J. M., Neale, M. C., & Kendler, K. S. (2001). A review and meta-analysis of the genetic epidemiology of anxiety disorders. American Journal of Psychiatry, 158(10), 1568–1578. https://doi.org/10.1176/appi.ajp.158.10.1568CrossRefGoogle ScholarPubMed
Hudson, J. I., Arnold, L. M., Keck, P. E., Auchenbach, M. B., & Pope, H. G. (2004). Family study of fibromyalgia and affective spectrum disorder. Biological Psychiatry, 56(11), 884–891. https://doi.org/10.1016/j.biopsych.2004.08.009CrossRefGoogle Scholar
Joustra, M. L., Janssens, K. A. M., Bültmann, U., & Rosmalen, J. G. M. (2015). Functional limitations in functional somatic syndromes and well-defined medical diseases. Results from the general population cohort LifeLines. Journal of Psychosomatic Research, 79(2), 94–99. https://doi.org/10.1016/j.jpsychores.2015.05.004CrossRefGoogle ScholarPubMed
Kato, K., Sullivan, P. F., Evengård, B., & Pedersen, N. L. (2009). A population-based twin study of functional somatic syndromes. Psychological Medicine, 39(3), 497–505. https://doi.org/10.1017/S0033291708003784CrossRefGoogle ScholarPubMed
Kendler, K. S. (1995). Is seeking treatment for depression predicted by a history of depression in relatives? Implications for family studies of affective disorder. Psychological Medicine, 25(4), 807–814. https://doi.org/10.1017/S0033291700035054CrossRefGoogle ScholarPubMed
Kendler, K. S. (2012). The dappled nature of causes of psychiatric illness: Replacing the organic–functional/hardware–software dichotomy with empirically based pluralism. Molecular Psychiatry, 17(4), 377–388. https://doi.org/10.1038/mp.2011.182CrossRefGoogle ScholarPubMed
Kendler, K. S., Gardner, C. O., Gatz, M., & Pedersen, N. L. (2007). The sources of co-morbidity between major depression and generalized anxiety disorder in a Swedish national twin sample. Psychological Medicine, 37(3), 453–462. https://doi.org/10.1017/S0033291706009135CrossRefGoogle Scholar
Kendler, K. S., & Neale, M. C. (2009). “Familiality” or heritability. Archives of General Psychiatry, 66(4), 452–453. https://doi.org/10.1001/archgenpsychiatry.2009.14CrossRefGoogle ScholarPubMed
Kendler, K. S., Rosmalen, J. G. M., Ohlsson, H., Sundquist, J., & Sundquist, K. (2023). A distinctive profile of family genetic risk scores in a Swedish national sample of cases of fibromyalgia, irritable bowel syndrome, and chronic fatigue syndrome compared to rheumatoid arthritis and major depression. Psychological Medicine, 53(9), 3879–3886. https://doi.org/10.1017/S0033291722000526CrossRefGoogle Scholar
Kleinstäuber, M., Schröder, A., Daehler, S., Pallesen, K. J., Rask, C. U., Sanyer, M., … Rosmalen, J. G. M. (2023). Aetiological understanding of fibromyalgia, irritable bowel syndrome, chronic fatigue syndrome and classificatory analogues: A systematic umbrella review. Clinical Psychology in Europe, 5(3), 1–25. https://doi.org/10.32872/cpe.11179CrossRefGoogle ScholarPubMed
Konnopka, A., Schaefert, R., Heinrich, S., Kaufmann, C., Luppa, M., Herzog, W., & König, H.-H. (2012). Economics of medically unexplained symptoms: A systematic review of the literature. Psychotherapy and Psychosomatics, 81(5), 265–275. https://doi.org/10.1159/000337349CrossRefGoogle ScholarPubMed
Lim, E.-J., Ahn, Y.-C., Jang, E.-S., Lee, S.-W., Lee, S.-H., & Son, C.-G. (2020). Systematic review and meta-analysis of the prevalence of chronic fatigue syndrome/myalgic encephalomyelitis (CFS/ME). Journal of Translational Medicine, 18(1), 100. https://doi.org/10.1186/s12967-020-02269-0CrossRefGoogle ScholarPubMed
Magnusson, K., Turkiewicz, A., Rydén, M., & Englund, M. (2024). Genetic influence on osteoarthritis versus other rheumatic diseases. Arthritis & Rheumatology, 76(2), 206–215. https://doi.org/10.1002/art.42696CrossRefGoogle ScholarPubMed
Markkula, R., Järvinen, P., Leino-Arjas, P., Koskenvuo, M., Kalso, E., & Kaprio, J. (2009). Clustering of symptoms associated with fibromyalgia in a Finnish twin cohort. European Journal of Pain (London, England), 13(7), 744–750. https://doi.org/10.1016/j.ejpain.2008.09.007CrossRefGoogle Scholar
Milne, B. J., Caspi, A., Crump, R., Poulton, R., Rutter, M., Sears, M. R., & Moffitt, T. E. (2009). The validity of the family history screen for assessing family history of mental disorders. American Journal of Medical Genetics. Part B, Neuropsychiatric Genetics, 150(1), 41–49. https://doi.org/10.1002/ajmg.b.30764CrossRefGoogle Scholar
Moscati, A., Faucon, A. B., Arnaiz-Yépez, C., Lönn, S. L., Sundquist, J., Sundquist, K., … Kendler, K. S. (2023). Life is pain: Fibromyalgia as a nexus of multiple liability distributions. American Journal of Medical Genetics. Part B, Neuropsychiatric Genetics: The Official Publication of the International Society of Psychiatric Genetics, 192(7–8), 171–182. https://doi.org/10.1002/ajmg.b.32949CrossRefGoogle ScholarPubMed
Muthén, L., & Muthén, B. (1998). Mplus user’s guide. (8th ed). Los Angeles, CA: Muthén & Muthén.Google Scholar
Nimnuan, C., Hotopf, M., & Wessely, S. (2001). Medically unexplained symptoms: An epidemiological study in seven specialities. Journal of Psychosomatic Research, 51(1), 361–367. https://doi.org/10.1016/s0022-3999(01)00223-9CrossRefGoogle ScholarPubMed
Picó-Pérez, M., Radua, J., Steward, T., Menchón, J. M., & Soriano-Mas, C. (2017). Emotion regulation in mood and anxiety disorders: A meta-analysis of fMRI cognitive reappraisal studies. Progress in Neuro-Psychopharmacology and Biological Psychiatry, 79, 96–104. https://doi.org/10.1016/j.pnpbp.2017.06.001CrossRefGoogle ScholarPubMed
Pinto, A. M., Geenen, R., Wager, T. D., Lumley, M. A., Häuser, W., Kosek, E., … da Silva, J. A. P. (2023). Emotion regulation and the salience network: A hypothetical integrative model of fibromyalgia. Nature Reviews Rheumatology, 19(1), 44–60. https://doi.org/10.1038/s41584-022-00873-6CrossRefGoogle ScholarPubMed
R Core Team. (2022). R: The R project for statistical computing. Vienna, Austria.: R Foundation for Statistical Computing. https://www.r-project.org/Google Scholar
Reich, T., James, J. W., & Morris, C. A. (1972). The use of multiple thresholds in determining the mode of transmission of semi-continuous traits. Annals of Human Genetics, 36(2), 163–184. https://doi.org/10.1111/j.1469-1809.1972.tb00767.xCrossRefGoogle ScholarPubMed
Risch, N. (1990). Linkage strategies for genetically complex traits. I. Multilocus models. American Journal of Human Genetics, 46(2), 222–228.Google ScholarPubMed
Rometsch, C., Mansueto, G., Maas Genannt Bermpohl, F., Martin, A., & Cosci, F. (2024). Prevalence of functional disorders across Europe: A systematic review and meta-analysis. European Journal of Epidemiology. https://doi.org/10.1007/s10654-024-01109-5CrossRefGoogle ScholarPubMed
Saito, Y. A. (2011). The role of genetics in IBS. Gastroenterology Clinics of North America, 40(1), 45–67. https://doi.org/10.1016/j.gtc.2010.12.011CrossRefGoogle ScholarPubMed
Saito, Y. A., Petersen, G. M., Larson, J. J., Atkinson, E. J., Fridley, B. L., de Andrade, M., … Talley, N. J. (2010). Familial aggregation of irritable bowel syndrome: A family case–control study. The American Journal of Gastroenterology, 105(4), 833–841. https://doi.org/10.1038/ajg.2010.116CrossRefGoogle ScholarPubMed
Saito, Y. A., Zimmerman, J. M., Harmsen, S., W., De Andrade, M., Locke, G. R. Iii, Petersen, G. M., & Talley, N. J. (2008). Irritable bowel syndrome aggregates strongly in families: A family-based case-control study. Neurogastroenterology & Motility, 20(7), 790–797. https://doi.org/10.1111/j.1365-2982.2007.01077.xCrossRefGoogle ScholarPubMed
Sheehan, D. V., Lecrubier, Y., Sheehan, K. H., Amorim, P., Janavs, J., Weiller, E., … Dunbar, G. C. (1998). The mini-international neuropsychiatric interview (M.I.N.I.): The development and validation of a structured diagnostic psychiatric interview for DSM-IV and ICD-10. The Journal of Clinical Psychiatry, 59(suppl 20), 11980.Google ScholarPubMed
Sijtsma, A., Rienks, J., van der Harst, P., Navis, G., Rosmalen, J. G. M., & Dotinga, A. (2022). Cohort profile update: Lifelines, a three-generation cohort study and biobank. International Journal of Epidemiology, 51(5), e295–e302. https://doi.org/10.1093/ije/dyab257CrossRefGoogle ScholarPubMed
Steinhausen, H.-C., Jakobsen, H., & Munk-Jørgensen, P. (2017). Family aggregation and risk factors in substance use disorders over three generations in a nation-wide study. PloS One, 12(5), e0177700. https://doi.org/10.1371/journal.pone.0177700CrossRefGoogle ScholarPubMed
Strawbridge, R., Sartor, M.-L., Scott, F., & Cleare, A. J. (2019). Inflammatory proteins are altered in chronic fatigue syndrome—A systematic review and meta-analysis. Neuroscience & Biobehavioral Reviews, 107, 69–83. https://doi.org/10.1016/j.neubiorev.2019.08.011CrossRefGoogle ScholarPubMed
Sullivan, P. F., Neale, M. C., & Kendler, K. S. (2000). Genetic epidemiology of major depression: Review and meta-analysis. American Journal of Psychiatry, 157(10), 1552–1562. https://doi.org/10.1176/appi.ajp.157.10.1552CrossRefGoogle ScholarPubMed
Svedberg, P., Johansson, S., Wallander, M.-A., & Pedersen, N. L. (2008). No evidence of sex differences in heritability of irritable bowel syndrome in Swedish twins. Twin Research and Human Genetics, 11(2), 197–203. https://doi.org/10.1375/twin.11.2.197CrossRefGoogle ScholarPubMed
Talley, N. J., & Spiller, R. (2002). Irritable bowel syndrome: A little understood organic bowel disease? Lancet (London, England), 360(9332), 555–564. https://doi.org/10.1016/S0140-6736(02)09712-XCrossRefGoogle ScholarPubMed
Tattan, M., Ørnbøl, E., Wellnitz, K. B., Hanssen, D. J. C., Dantoft, T. M., Rosmalen, J. G. M., … Petersen, M. W. (2024). Factors associated with having previously received a diagnosis of fibromyalgia, chronic fatigue syndrome and irritable bowel syndrome: A cross sectional DanFunD study. Journal of Psychosomatic Research, 181, 111693. https://doi.org/10.1016/j.jpsychores.2024.111693CrossRefGoogle ScholarPubMed
Tavares, L. C., Lopera-Maya, E. A., Bonfiglio, F., Zheng, T., Sinha, T., Marques, F. Z., … D’Amato, M. (2024). Rome III criteria capture higher irritable bowel syndrome SNP-heritability and highlight a novel genetic link with cardiovascular traits, Cellular and Molecular Gastroenterology and Hepatology, 18(2), 101345. https://doi.org/10.1016/j.jcmgh.2024.04.002CrossRefGoogle Scholar
Tesfaye, M., Jaholkowski, P., Hindley, G. F. L., Shadrin, A. A., Rahman, Z., Bahrami, S., … Andreassen, O. A. (2023). Shared genetic architecture between irritable bowel syndrome and psychiatric disorders reveals molecular pathways of the gut-brain axis. Genome Medicine, 15(1), 60. https://doi.org/10.1186/s13073-023-01212-4CrossRefGoogle ScholarPubMed
Thomas, N. S., Gillespie, N. A., Kendler, K. S., Oldehinkel, A. J., Rosmalen, J. G. M., & van Loo, H. M. (2024). Comorbidity and sex differences in functional disorders and internalizing disorders. General Hospital Psychiatry, 90, 91–98. https://doi.org/10.1016/j.genhosppsych.2024.07.013CrossRefGoogle ScholarPubMed
van der Meulen, M. L., Bos, M., Bakker, S. J. L., Gans, R. O. B., & Rosmalen, J. G. M. (2024). Validity and diagnostic overlap of functional somatic syndrome diagnoses. Journal of Psychosomatic Research, 181, 111673. https://doi.org/10.1016/j.jpsychores.2024.111673CrossRefGoogle ScholarPubMed
van Loo, H. M., Beijers, L., Wieling, M., de Jong, T. R., Schoevers, R. A., & Kendler, K. S. (2023). Prevalence of internalizing disorders, symptoms, and traits across age using advanced nonlinear models. Psychological Medicine, 53(1), 78–87. https://doi.org/10.1017/S0033291721001148CrossRefGoogle ScholarPubMed
Wang, R., Snieder, H., & Hartman, C. A. (2022). Familial co-aggregation and shared heritability between depression, anxiety, obesity and substance use. Translational Psychiatry, 12(1), 108. https://doi.org/10.1038/s41398-022-01868-3CrossRefGoogle ScholarPubMed
Warren, J. W., & Clauw, D. J. (2012). Functional somatic syndromes: Sensitivities and specificities of self-reports of physician diagnosis. Psychosomatic Medicine, 74(9), 891–895. https://doi.org/10.1097/PSY.0b013e31827264aaCrossRefGoogle ScholarPubMed
Wolfe, F., Clauw, D. J., Fitzcharles, M.-A., Goldenberg, D. L., Katz, R. S., Mease, P., … Yunus, M. B. (2010). The American College of Rheumatology preliminary diagnostic criteria for fibromyalgia and measurement of symptom severity. Arthritis Care & Research, 62(5), 600–610. https://doi.org/10.1002/acr.20140CrossRefGoogle ScholarPubMed
Wolfe, F., Schmukler, J., Jamal, S., Castrejon, I., Gibson, K. A., Srinivasan, S., … Pincus, T. (2019). Diagnosis of fibromyalgia: Disagreement between fibromyalgia criteria and clinician-based fibromyalgia diagnosis in a university clinic. Arthritis Care & Research, 71(3), 343–351. https://doi.org/10.1002/acr.23731CrossRefGoogle ScholarPubMed
Wray, N. R., & Gottesman, I. I. (2012). Using summary data from the Danish national registers to estimate heritabilities for schizophrenia, bipolar disorder, and major depressive disorder. Frontiers in Genetics, 3, 118. https://doi.org/10.3389/fgene.2012.00118CrossRefGoogle ScholarPubMed
Zeileis, A., Köll, S., & Graham, N. (2020). Various versatile variances: An object-oriented implementation of clustered covariances in R. Journal of Statistical Software, 95, 1–36. https://doi.org/10.18637/jss.v095.i01CrossRefGoogle Scholar
Figure 0

Table 1. Descriptive characteristics of dataset

Figure 1

Figure 1. Familial coaggregation of internalizing and functional disorders amongst first-degree relatives expressed in λR. Note: FDR, first-degree relative; MDD, major depressive disorder; DYS, dysthymia; GAD, generalized anxiety disorder; AGPH, agoraphobia; SPH, social phobia; PD, panic disorder; ME/CFS, myalgic encephalomyelitis/chronic fatigue syndrome; IBS, irritable bowel syndrome; FM, fibromyalgia; λR, recurrence risk ratios. λR adjusted for age, age2, sex, and number of relatives present in the data. Estimates in bold are significant at p < .01. For 95% confidence intervals, see Supplementary Table 1.

Figure 2

Figure 2. Tetrachoric correlations between internalizing and functional disorders, for (a) siblings and (b) parent-offspring pairs/trios. Note: MDD, major depressive disorder; DYS, dysthymia; GAD, generalized anxiety disorder; AGPH, agoraphobia; SPH, social phobia; PD, panic disorder; ME/CFS, myalgic encephalomyelitis/chronic fatigue syndrome; IBS, irritable bowel syndrome; FM, fibromyalgia; rtet, tetrachoric correlation. Estimates in bold are significant at p < .01. For 95% confidence intervals, see Supplementary Tables 2 & 3.

Figure 3

Figure 3. Familial coaggregation of internalizing and functional disorders amongst spouses, expressed in λR. Note: MDD, major depressive disorder; DYS, dysthymia; GAD, generalized anxiety disorder; AGPH, agoraphobia; SPH, social phobia; PD, panic disorder; ME/CFS, myalgic encephalomyelitis/chronic fatigue syndrome; IBS, irritable bowel syndrome; FM, fibromyalgia; λR, recurrence risk ratio. λR adjusted for age, age2, sex, and presence of spouse in the data. Estimates in bold are significant at p < .01. For 95% confidence intervals, see Supplementary Table 1.

Figure 4

Figure 4 (a) Familiality and (b) familial correlation estimates of internalizing and functional disorders based on both first- and second-degree relatives. Note: MDD, major depressive disorder; DYS, dysthymia; GAD, generalized anxiety disorder; AGPH, agoraphobia; SPH, social phobia; PD, panic disorder; ME/CFS, myalgic encephalomyelitis/chronic fatigue syndrome; IBS, irritable bowel syndrome; FM, fibromyalgia; rf, familial correlation. See Supplementary Table 4 for familiality estimates for first- and second-degree relatives separately and Supplementary Table 5 for 95% confidence intervals of familial correlation estimates.

Supplementary material: File

Bos et al. supplementary material 1

Bos et al. supplementary material
Download Bos et al. supplementary material 1(File)
File 495.7 KB
Supplementary material: File

Bos et al. supplementary material 2

Bos et al. supplementary material
Download Bos et al. supplementary material 2(File)
File 55.3 KB