Skip to main content Accessibility help


  • Access
  • Cited by 3


      • Send article to Kindle

        To send this article to your Kindle, first ensure is added to your Approved Personal Document E-mail List under your Personal Document Settings on the Manage Your Content and Devices page of your Amazon account. Then enter the ‘name’ part of your Kindle email address below. Find out more about sending to your Kindle. Find out more about sending to your Kindle.

        Note you can select to send to either the or variations. ‘’ emails are free but can only be sent to your device when it is connected to wi-fi. ‘’ emails can be delivered even when you are not connected to wi-fi, but note that service fees apply.

        Find out more about the Kindle Personal Document Service.

        Aetiological overlap between obsessive–compulsive related and anxiety disorder symptoms: Multivariate twin study
        Available formats

        Send article to Dropbox

        To send this article to your Dropbox account, please select one or more formats and confirm that you agree to abide by our usage policies. If this is the first time you use this feature, you will be asked to authorise Cambridge Core to connect with your <service> account. Find out more about sending content to Dropbox.

        Aetiological overlap between obsessive–compulsive related and anxiety disorder symptoms: Multivariate twin study
        Available formats

        Send article to Google Drive

        To send this article to your Google Drive account, please select one or more formats and confirm that you agree to abide by our usage policies. If this is the first time you use this feature, you will be asked to authorise Cambridge Core to connect with your <service> account. Find out more about sending content to Google Drive.

        Aetiological overlap between obsessive–compulsive related and anxiety disorder symptoms: Multivariate twin study
        Available formats
Export citation



The aetiological boundary between obsessive–compulsive related disorders (OCRDs) including obsessive–compulsive disorder (OCD) and anxiety disorders is unclear and continues to generate debate.


To determine the genetic overlap and the pattern of causal relationships among OCRDs and anxiety disorders.


Multivariate twin modelling methods and a new regression analysis to infer causation were used, involving 2495 male and female twins.


The amount of common genetic liability observed for OCD symptoms was higher when considering anxiety disorders and OCRDs in the model v. modelling OCRD symptoms alone. OCD symptoms emerged as risk factors for the presence of generalised anxiety, panic and hoarding symptoms, whereas social phobia appeared as a risk factor for OCD symptoms.


OCD represents a complex phenotype that includes important shared features with anxiety disorders and OCRDs. The novel patterns of risk identified between OCD and anxiety disorder may help to explain their frequent co-occurrence.


This study was funded by an Early Career Researcher Grant from The University of Melbourne to B.J.H. B.J.H. is supported by a National Health and Medical Research Council of Australia (NHMRC) Clinical Career Development Fellowship (I.D. 628509). C.L.-S. is supported by the Spanish Ministry of Education, Culture and Sport (FPU12/01636). C.P. was supported by a NHMRC Senior Principal Research Fellowship (ID: 628386) and NARSAD Distinguished Investigator Award. M.Y. was supported by a NHMRC Senior Research Fellowship (ID: 1021973). This research was facilitated through access to the Australian Twin Registry, a national resource supported by an Enabling Grant (ID 628911) from the National Health and Medical Research Council. None of these funding bodies had any role in the design and conduct of the study; collection, management, analysis and interpretation of the data; preparation, review, or approval of the manuscript; or decision to submit the manuscript for publication.

Declaration of interest


There has been much debate 15 regarding the optimal diagnostic classification of obsessive–compulsive disorder (OCD) in relation to the recently published DSM-5. 6 This revision has seen OCD removed from the broad category of anxiety disorders and placed at the centre of a new separate chapter – obsessive–compulsive and related disorders (OCRDs) – including body dysmorphic disorder (BDD), trichotillomania (hair-pulling disorder), as well as excoriation (skin picking) and hoarding disorder as new diagnoses. Based upon available evidence, it has been argued that OCD shares a stronger commonality with these disorders in terms of its core phenomenological, neurobiological and treatment characteristics. 2,710 However, by implication, the notion that OCD has less in common with anxiety disorders or that anxiety is less relevant to OCD remains contested. 3 Multivariate twin studies are particularly well-suited for addressing whether co-occurring mental disorders, such as OCRDs and anxiety disorders, demonstrate overlap in their genetic and environmental risk factors. 11 Of the few existing studies to have directly examined OCD with other anxiety disorders in adults, OCD was reported to show the highest percentage of specific genetic risk factors (45%), although OCD and anxiety disorders demonstrate a substantial common genetic liability (55%). 12 In a recent study that evaluated dimensional representations of OCD and OCRD symptoms (i.e. symptoms in a normative twin population), two distinct genetic liability factors were identified: one factor primarily representing OCD symptoms, hoarding disorder and BDD symptoms; and a second factor representing trichotillomania and skin picking symptoms. 13 In this same cohort it has been reported that 64% of the total covariance between OCD and BDD was explained by shared genetic factors. 14 In summary, prior twin studies suggest that OCD most likely shares genetic factors with both anxiety disorders and certain OCRDs. However, we highlight that no study to date has directly compared OCD with OCRDs and anxiety disorder symptoms in the same twin population assessed at the same period of time. Such comparisons will be important for addressing whether OCD is more or less aetiologically aligned with OCRDs v. anxiety disorders.

Our primary aim was to therefore more thoroughly assess the structure of genetic and environmental risk factors for dimensional representations of OCD, other OCRD and anxiety disorder symptoms in an adult twin population using classical multivariate twin modelling. Our second aim was to complement this approach with a recently introduced regression-based twin analysis, which allows inferences to be made about potential causal relationships between predictor and outcome variables. 15 With both approaches we sought to address the debated question of whether OCD symptoms are more or less aetiologically aligned with the symptoms of other OCRDs or anxiety disorders. If OCD is more aligned with certain OCRDs v. anxiety disorders, then evidence for a strong common genetic liability should be minimal in a multivariate analysis that combines OCRD and anxiety disorder domains. Further, if this ‘weak’ pattern of liability is evident, there should be minimal evidence for potential causal influences between OCD and anxiety disorder symptoms.


Participants and measures

A total of 2495 voluntary twin members (18 to 45 years old; mean 34.5 years (s.d. = 7.8) and 33.9 (s.d. = 8) monozygotic (MZ) and dizygotic (DZ) respectively) of the Australian Twin Registry (ATR) were recruited for this email-based online survey (1281 MZ and 1214 DZ twins). Online Fig. DS1 provides a schematic overview of the whole sample. All participants gave informed consent after receiving complete information on the study and before starting to fill out the survey. Full recruitment details are provided in López-Solà et al. 16 To address this study's aims, we focused specifically on twins' responses to six validated dimensional self-report measures of OCD, OCRD and anxiety disorder symptoms. To address our particular aims, the inclusion of three symptom domains per diagnostic category was considered optimal with respect to the planned multivariate model fitting analyses. OCD symptoms were estimated with the Obsessive Compulsive Inventory-Revised (OCI-R). 17 For OCRDs, we assessed hoarding disorder symptoms with the Hoarding Rating Scale-Self Report (HRS-SR) 18 and BDD symptoms with the Dysmorphic Concern Questionnaire (DCQ). 19 For anxiety disorders, we assessed social phobia symptoms with the Social Phobia Inventory (SPIN); 20 panic disorder symptoms with the Anxiety Sensitivity Index (ASI); 21 and generalised anxiety disorder (GAD) symptoms with the ‘Stress’ subscale of the Depression, Anxiety and Stress Scale-21 (DASS-21). 22

Statistical analysis

To ensure data normality, all questionnaire responses underwent Box–Cox transformations (yλ t = (yλ−1)/λ). 16,23 To examine the phenotypic structure of OCRD and anxiety disorder symptoms prior to classical twin modelling analysis, a varimax-rotated principal component analysis (PCA) was performed in SPSS (version 20), adjusting the transformed total scores of each scale for age and gender. In order not to confound the comparison of the OCD symptoms and hoarding disorder symptoms measured by the HRS-SR, estimated total scores on the OCI-R excluded twins' responses to the hoarding subscale of this measure (items 1, 7 and 13).

Classical multivariate twin modelling

Structural equation models were conducted on transformed continuous variables fitted by maximum likelihood. Because univariate twin modelling of this data indicated the presence of genetic gender differences, 16 all multivariate twin models were performed using standardised residual values for each symptom domain adjusted for age and gender. 24 Although controlling for the influence of gender in this manner has been a useful approach in multivariate twin studies, this is not the same as testing for multivariate genetic gender differences in which the variance–covariance structure of the model is allowed to be different for males and females. We chose to control rather than test for multivariate genetic gender differences (a) because of an absence of specific hypothesis regarding gender-related multivariate heterogeneity; and (b) because of the acknowledged additional complexity in fitting such models. 25

MZ and DZ cross-twin–within/cross-symptom correlations were estimated for each symptom domain by fitting the data to a constrained saturated model. The structure of genetic and environmental influences on OCRD and anxiety disorder symptoms was then estimated using classical multivariate twin models. 11 Model 1 corresponded to a fully saturated Cholesky decomposition that estimated 1A (additive genetic), 1C (shared environment) and 1E (non-shared environment) factors for each phenotype, making no assumptions about the nature of their underlying covariance. Model 2 corresponded to an independent pathway model, which seeks to estimate a set of common Ac, Cc, and Ec factors hypothesised to directly influence all phenotypes v. specific As, Cs and Es factors that may explain remaining phenotypic variance. Model 3 corresponded to a common pathway model, which estimated whether the covariance among phenotypes was influenced via one latent factor taking into account the shared contribution of common A, C and E factors. The Akaike information criterion (AIC) value was used to provide a measure of the goodness of fit of these models. Reduced submodels were systematically tested to derive the most parsimonious model fitting results. Classical twin modelling was carried out in R ( using the OpenMx package. 26

Identifying potential causal influences

Inference on causation from examination of familial confounding (ICE FALCON) is a regression-based approach for analysing twin pair data on a continuously or dichotomously distributed outcome and a familial predictor measured for both the twin and his or her co-twin. 15 The underlying statistical model allows one to make inferences about (but of course does not ‘prove’) the existence of causal relationships between predictor and outcome variables via the elimination of familial confounding. 27 This can also be thought of as using the co-twin as a ‘negative control’. If the predictor is familial – that is, it is strongly correlated in twins – and there is at least in part a causal relationship between the predictor and the outcome, the association between the predictor of twin A with the outcome of twin B will decrease in absolute strength towards the null after including the predictor of twin B in the model. 27 In other words, this model considers the evidence for confounding because of genetic and/or environmental factors shared by both twins (Ac, C and Ec factors) v. non-shared factors (As and Es). Accordingly, when evidence consistent with a ‘casual’ association is identified, the role of participant-specific factors is emphasised, although the relative contribution of specific genetic v. environmental influences is not determined. If the associations between the outcome of twin A and the predictors of both the twin A and co-twin B are the same before and after adjusting for each other, then, under this model, there is no evidence consistent with a potential causal relationship. On the other hand, if there is a significant attenuation of the cross-trait cross-pair association after conditioning on twin A, there is evidence ‘consistent with’ some causation. In our analyses, we focused on the relationship between OCD symptoms and each of the OCRD and anxiety disorder symptom domains. A step-by-step explanation of this approach is provided in online Fig. DS2 and online supplement DS1.


Correlation and PCA

Moderate-strength phenotypic correlations were observed across all symptom domains. The strongest correlations with OCD were found with anxiety disorder symptoms; the weakest associations were observed between hoarding disorder and BDD symptoms, and hoarding disorder and anxiety disorder symptoms (0.32–0.37). The pattern of cross-twin cross-trait correlations in MZ compared with DZ twins supports a relevant genetic component to the liability of each domain and their co-occurrence (Table 1). PCA retained one phenotypic factor with an eigenvalue greater than 1 and explaining 55.6% of the total variance. When forcing it to retain 2 factors, the total explained variance increased to 67.8%, and hoarding disorder symptoms emerged as a distinct factor (see online Table DS1).

Table 1 Phenotypic, cross-twin within-trait (diagonal) and cross-twin cross-trait (off-diagonal) correlations for monozygotic (MD) and dizygotic (DZ) male and female twins with 95% confidence intervals a


disorder symptoms
Hoarding disorder

Body dysmorphic disorder

Panic disorder

Generalised anxiety

disorder symptoms
Social phobia

Phenotypic correlations
    Obsessive–compulsive disorder symptoms 1
    Hoarding disorder symptoms 0.42 (0.39 to 0.46) 1
    Body dysmorphic disorder symptoms 0.44 (0.41 to 0.47) 0.32 (0.28 to 0.35) 1
    Panic disorder symptoms 0.56 (0.54 to 0.59) 0.36 (0.32 to 0.39) 0.48 (0.45 to 0.51) 1
    Generalised anxiety disorder symptoms 0.55 (0.52 to 0.57) 0.36 (0.32 to 0.39) 0.42 (0.39 to 0.45) 0.56 (0.53 to 0.58) 1
    Social phobia symptoms 0.51 (0.48 to 0.54) 0.37 (0.33 to 0.40) 0.48 (0.45 to 0.51) 0.61 (0.59 to 0.64) 0.51 (0.47 to 0.53) 1
Twin correlations
    Obsessive–compulsive disorder symptoms 0.38 (0.31 to 0.45)

0.14 (0.04 to 0.24)
0.20 (0.15 to 0.26) 0.23 (0.17 to 0.29) 0.25 (0.18 to 0.30) 0.28 (0.22 to 0.33) 0.25 (0.19 to 0.31)
    Hoarding disorder symptoms 0.02 (−0.05 to 0.09) 0.33 (0.25 to 0.40)

0.13 (0.04 to 0.22)
0.16 (0.10 to 0.22) 0.16 (0.10 to 0.22) 0.19 (0.13 to 0.25) 0.18 (0.13 to 0.24)
    Body dysmorphic disorder symptoms 0.08 (0.001 to 0.15) 0.03 (−0.04 to 0.10) 0.40 (0.32 to 0.47)

0.08 (−0.02 to 0.17)
0.22 (0.16 to 0.28) 0.19 (0.13 to 0.25) 0.24 (0.18 to 0.30)
    Panic disorder symptoms 0.12 (0.04 to 0.19) 0.05 (−0.02 to 0.12) 0.07 (−0.002 to 0.15) 0.30 (0.22 to 0.37)

0.14 (0.05 to 0.23)
0.24 (0.17 to 0.29) 0.26 (0.20 to 0.32)
    Generalised anxiety disorder symptoms 0.13 (0.05 to 0.20) 0.05 (−0.02 to 0.12) 0.08 (0.01 to 0.15) 0.11 (0.04 to 0.18) 0.33 (0.25 to 0.40)

0.17 (0.07 to 0.25)
0.21 (0.15 to 0.27)
    Social phobia symptoms 0.09 (0.01 to 0.16) 0.06 (−0.01 to 0.13) 0.06 (−0.02 to 0.13) 0.16 (0.08 to 0.23) 0.12 (0.04 to 0.18) 0.42 (0.35 to 0.48)

0.18 (0.08 to 0.26)

a. Cross-twin cross-trait correlations for MZ are given in the area above the diagonal and for DZ twins in the area below the diagonal.

Multivariate twin modelling

The independent pathway model provided an improved fit, compared with the fully saturated Cholesky with reduced parameters (Table 2). However, the more restrictive common pathway model resulted in a significantly worse fit with an increased AIC value. A series of independent pathway nested submodels were fitted to test the importance of specific parameters compared with the fully saturated Cholesky and the independent pathway model. In these submodels, the genetic and environmental liabilities were either forced to be entirely common/shared (models 4–6) or entirely independent/specific (models 7–9). Model 4 was the best-fitting such that the covariation between phenotypes was explained by a set of common Ac, Cc and Ec factors, and the remaining variance by As and Es effects specific to each phenotype. Removing the Cs factor did not lead to worse fit, suggesting that these factors were less important in explaining individual differences. Removing As +/− Cs (models 5 and 6) and forcing all genetic risk to be common led to significantly worse fit, suggesting that some of the genetic liability is specific to each phenotype. Models 7 to 9 were forced to have an As, Cs and Es specific to each phenotype, and none of them demonstrated good fit, indicating the existence of common liability to all phenotypes.

Table 2 Model-fitting results

Estimated parameter Fit statistic Comparison of nested models
Model Common factors Specific factors −2LL d.f. AIC χ2 Δd.f. P Compared with Δχ2 (Δd.f.) P
1 Cholesky saturated ACE ACE 30357.3 14901 555.3
2 Independent pathway ACE ACE 30397.4 14928 541.4 40.1 27 0.05
3 Common pathway ACE ACE 30450.1 14938 574.1 92.9 37 <0.001 2 52.8 (10) <0.001
4 Independent pathway a ACE AE 30397.7 14934 529.7 40.5 33 0.17 2 0.35 (6) 1
5 Independent pathway ACE CE 30427.3 14934 559.3 70.1 33 <0.001 2 29.95 (6) <0.001
6 Independent pathway ACE E 30510.3 14940 630.3 153.0 39 <0.001 2 112.89 (12) <0.001
7 Independent pathway AE ACE 30435.3 14934 567.3 78.0 33 <0.001 2 37.91 (6) <0.001
8 Independent pathway CE ACE 30463.8 14934 595.8 106.6 33 <0.001 2 66.47 (6) <0.001
9 Independent pathway E ACE 30557.1 14940 674.1 196.8 39 <0.001 2 156.69 (12) <0.001

−2LL, minus twice the log-likelihood; Δd.f., change in degrees of freedom between the submodel and the full model; Δχ2, difference in goodness-of-fit statistic between the submodel and the full model; A, additive genetic factor; C, shared environmental factor; E, non-shared environmental factor.

a. The best-fitting model based on the lower Akaike information criterion (AIC).

Parameter estimates for the best-fitting model are presented in Table 3 and the independent pathway model in Fig. 1(a). Although the Cc factor could not be dropped from the model, it accounted for a minor fraction of the overall variance (0.2–14%). Figure 1(b) shows that OCD symptoms were the only phenotype to share almost all of its additive genetic influence with the remaining OCRD and anxiety disorder symptoms, suggesting that OCD symptoms are aetiologically related to all of these phenotypes. For hoarding disorder and BDD symptoms, 55% and 61% of the total genetic variance, respectively, was as a result of As, suggesting more specific genetic risk factors for these symptoms v. other anxiety disorder and/or OCRD symptoms. For these reasons, we conducted a separate multivariate analysis with two Ac latent factors: one loading on all symptom domains and another loading only on OCD, hoarding disorder and BDD symptoms. Similar to the former model, we observed that the majority of the genetic variance for OCD symptoms still loaded onto the first Ac factor sharing genetic effects with other anxiety disorders and OCRDs. Only 9% of the total genetic variance of OCD was as a result of genetic factors shared with other OCRDs. For BDD symptoms 50% of its total genetic variance loaded onto the specific genetic factor (As), whereas 1.4% loaded onto shared genetic factors with the other OCRDs. A total of 37.5% of the total genetic variance of hoarding symptoms loaded onto shared genetic factors with the other OCRDs, with 31.3% as a result of specific genetic influences. Results of the two-factor independent pathway model are given in online Table DS2.

Table 3 Standardised parameters for the best-fitting Model 4 (with 95% confidence intervals)

Additive genetic Shared environment Non-shared environment
Common Specific Total Common Specific Common Specific Total
Obsessive–compulsive disorder symptoms 0.37 (0.24–0.44) 0.004 (0.00–0.08) 0.37 (0.25–0.44) 0.002 (0.00–0.10) 0.25 (0.19–0.31) 0.37 (0.33–0.42) 0.62 (0.56–0.69)
Hoarding disorder symptoms 0.15 (0.07–0.22) 0.18 (0.12–0.24) 0.33 (0.24–0.40) 0.006 (0.00–0.06) 0.12 (0.07–0.17) 0.55(0.49–0.62) 0.668 (0.60–0.74)
Body dysmorphic disorder symptoms 0.12 (0.03–0.22) 0.18 (0.13–0.24) 0.30 (0.19–0.40) 0.05 (0.00–0.13) 0.22 (0.16–0.28) 0.44 (0.39–0.50) 0.657 (0.58–0.74)
Panic disorder symptoms 0.15 (0.03–0.27) 0.04 (0.00–0.08) 0.19 (0.04–0.31) 0.09 (0.003–0.22) 0.396 (0.32–0.47) 0.33 (0.28–0.38) 0.725 (0.65–0.80)
Generalised anxiety disorder symptoms 0.20 (0.09–0.28) 0.10 (0.06–0.15) 0.30 (0.19–0.38) 0.02 (0.00–0.11) 0.28 (0.23–0.35) 0.40 (0.35–0.45) 0.68 (0.61–0.75)
Social phobia symptoms 0.13 (0.01–0.30) 0.12 (0.05–0.18) 0.25 (0.10–0.40) 0.14 (0.02–0.26) 0.326 (0.26–0.40) 0.29 (0.25–0.34) 0.616 (0.55–0.69)

Fig. 1 (a) Independent pathway (best-fitting) model. (b) and (c) The percentage of the variance accounted for by common and specific genetic and non-shared environmental factors.

(a) Ac (symptom-common genetic influence), Cc (symptom-common shared environmental influence) and Ec (symptom-common non-shared environmental influence). The breakdown of the genetic and non-shared environmental variance into common and specific factors is shown in (b) A (Ac and As) and (c) E (Ec and Es), respectively. OCD, obsessive–compulsive disorder symptoms; HD, hoarding disorder symptoms; BDD, body dysmorphic disorder symptoms; PD, panic disorder symptoms; GAD, generalised anxiety disorder symptoms; SP, social phobia symptoms; OCRD, obsessive–compulsive and related disorders.

Potential causal influences

As depicted in Fig. 2 and Table 4, there was evidence for significant causal influences between OCD and anxiety disorder symptom domains. Specifically, OCD symptoms demonstrate a significant causal influence on GAD and panic disorder symptoms, respectively (P<0.0001). In other words, it can be inferred that there is a high probability of observing changes in the severity of these latter domains when a person's OCD symptom severity changes, but not vice versa. By comparison, we observed a significant causal influence of social phobia symptoms on OCD symptoms (P = 0.03), suggesting that OCD symptoms themselves show some dependency on social phobia symptoms. Social phobia was the only symptom domain to demonstrate a causal influence on OCD symptoms.

Fig. 2 Causal modelling with inference on causation from examination of familial confounding (ICE FALCON).

P-values refer to the significance of the regression coefficient change between Model II to Model III. The direction of the arrows indicates the direction of estimated ‘causality’. The light and dark colouring symbolise DSM-5 representations of obsessive–compulsive and related disorders and anxiety disorders respectively. ns, not significant.

Table 4 Inference on causation from examination of familial confounding (ICE FALCON) after testing the probability for both directions of causation between each pair of symptoms a

Model I Model II Model III Change
Estimate s.e. P Estimate s.e. P Estimate s.e. P Absolute (%) P (1-side)
OCD (X)–hoarding disorder symptoms (Y) b
    Self −0.438 0.03 <0.0001 −0.430 0.030 <0.0001
    Co-twin −0.132 0.036 0.0002 −0.050 0.030 0.10 −0.082 (61) 0.01
    Self 0.440 0.029 <0.0001 0.434 0.028 <0.0001
    Co-twin 0.122 0.036 0.0007 0.069 0.029 0.02 −0.052 (43) 0.07
OCD (X)–panic disorder symptoms (Y)
    Self 0.592 0.027 <0.0001 0.584 0.027 <0.0001
    Co-twin 0.253 0.039 <0.0001 0.037 0.028 0.19 −0.215 (85.2) <0.0001
    Self 0.557 0.025 <0.0001 0.539 0.027 <0.0001
    Co-twin 0.266 0.035 <0.0001 0.081 0.027 0.003 −0.186 (69.7) 0.0001
OCD (X)–social phobia symptoms (Y)
    Self 0.495 0.029 <0.0001 0.493 0.029 <0.0001
    Co-twin 0.092 0.039 0.02 0.069 0.028 0.01 −0.023 (25) 0.33
Hoarding disorder symptoms (X)–OCD (Y) b
    Self −0.413 0.030 <0.0001 −0.414 0.029 <0.0001
    Co-twin −0.074 0.035 0.03 −0.082 0.028 0.004 0.008 (10.4) 0.42
    Self 0.452 0.029 <0.0001 0.447 0.028 <0.0001
    Co-twin 0.117 0.037 0.001 0.077 0.028 0.006 −0.040 (34.5) 0.16
Panic disorder symptoms (X)–OCD (Y)
    Self 0.570 0.027 <0.0001 0.569 0.027 <0.0001
    Co-twin 0.109 0.040 0.007 0.094 0.027 0.0005 −0.015 (13.7) 0.42
    Self 0.544 0.027 <0.0001 0.534 0.027 <0.0001
    Co-twin 0.190 0.037 <0.0001 0.118 0.027 <0.0001 −0.072 (37.8) 0.11
Social phobia symptoms (X)–OCD (Y)
    Self 0.511 0.030 <0.0001 0.505 0.031 <0.0001
    Co-twin 0.139 0.039 0.0004 0.047 0.030 0.11 −0.092 (66.2) 0.03

OCD, obsessive–compulsive disorder symptoms; X, independent variable; Y, outcome variable; BDD, body dysmorphic disorder symptoms; GAD, generalised anxiety disorder symptoms.

a. Regression estimates, standard error and P-value from regression models on monozygotic twin pairs.

b. The inverse relationship between OCD and hoarding disorder symptoms was because of the fact that we transformed the raw data using reciprocal function for hoarding disorder symptoms (1/hoarding disorder symptoms0.55) and logarithm function for OCD.

Within the OCRD domain, evidence of a significant causal influence of OCD symptoms was observed on hoarding disorder symptoms (P = 0.01), suggesting that any change in OCD symptom severity would predict a corresponding change in hoarding disorder severity, but not vice versa. There was trend evidence suggesting a causal influence of OCD symptoms on BDD symptoms (P = 0.07), but not vice versa.


In summary, our results do not support the contention that OCD symptoms are less aetiologically aligned with the symptoms of anxiety disorders compared with some of the revised DSM-5 OCRDs. On the basis of classical twin modelling, evidence of a genetic commonality between OCD and anxiety disorder symptoms was observed, such that the genetic liability to OCD symptoms was better explained when modelling its shared liability with OCRD and anxiety disorder symptoms compared with when modelling an additional OCRD latent genetic factor. On the basis of causal inference analysis we observed evidence consistent with OCD symptoms being a potential causal risk factor for panic disorder, GAD and hoarding disorder symptoms, in the sense that having OCD symptoms appears to increase the probability of having panic disorder and GAD symptoms, but not vice versa. By comparison, we observed evidence consistent with social phobia symptoms being a potential causal risk factor for OCD symptoms. We could reject the alternate (null) hypothesis that there are no direct causal relations.

A non-specific genetic vulnerability to OCD symptoms

Our twin modelling results indicate that the genetic liability to OCD symptoms was almost entirely shared with the five other symptom domains. Of these domains, panic disorder, GAD and social phobia symptoms also demonstrated greater common vs. specific genetic liabilities. When taken together with previous twin studies, there is now good evidence to suggest that OCD is influenced by moderately heritable genetic factors that are mostly shared with other OCRDs and anxiety disorders. 1214,28,29 Thus, it stands to reason that this common liability may partly underlie the co-occurrence of these disorders in terms of their high rates of comorbidity 30 and familial aggregation. 5,31 Phenomenologically, in addition to the general characteristic of heightened threat estimation, 3 there are other underlying features that link OCD and anxiety disorder symptoms. For instance, high levels of self-blaming emotions, such as guilt and shame, appear to be shared between OCD and other anxiety disorders, including social phobia. 32 Additionally, although most robustly linked to panic disorder, heightened anxiety sensitivity is observed in patients with OCD, GAD and social phobia, which may reflect a common cognitive bias towards ‘over-importance of thoughts’ – a recognised dimension of anxiety sensitivity. 33,34 Common deficits of attentional control have also been emphasised in relation to OCD and other anxiety disorders, particularly GAD, as a feature that may explain the pervasive negative cognitions (i.e. obsessions, worry) that characterise these disorders. 35 Of course, although such features are not characteristic of all patients with OCDs, nor patients with other anxiety disorders, they nonetheless appear to represent important transdiagnostic ‘common threads’ that may in part arise from such estimated common liabilities.

OCD and anxiety disorder symptoms as putative causal risk factors

Extending the classical modelling approach, causal pathway modelling has identified novel and potentially important relationships between OCD and anxiety disorder symptoms. With this approach, the results were more consistent with OCD symptoms being a causal risk factor for panic disorder and GAD symptoms (having OCD symptoms increases the probability of having panic disorder and GAD symptoms, not vice versa), rather than the association between traits being as a result of unmeasured familial factors such as genes and shared environment. In other studies without twins, OCD checking symptoms – one of the most common OCD symptom dimensions – have been linked to the increased probability of comorbid panic disorder and GAD diagnoses. 36 It was proposed that an intolerance of uncertainty, which is present in OCD checking symptoms, but more strongly associated with general worry and anxiety, might explain the link between these domains. Thus, one possibility is that intolerance of uncertainty represents an underlying trait dimension through which OCD aetiologically enhances the risk of panic disorder and GAD symptoms.

By contrast, social phobia symptoms emerged as a potential causal risk factor for OCD symptoms, although this was a more moderate finding compared with the estimated strength of the above associations. Nevertheless, adopting similar logic, one possibility is that a relevant underlying trait factor may explain these findings. In the case of social phobia symptoms, one obvious candidate would be ‘behavioural inhibition’. Although behavioural inhibition has been most strongly characterised as a childhood predictor of social phobia, 37 it has also been linked to the development of OCD symptoms in adulthood. 38 Supporting an early aetiological link between these domains, mother-reported levels of inhibition/shyness in preschool-aged twins were reported to show substantial overlap with other anxiety-related behaviours, including OCD-like behaviours. 39 One hypothesis may therefore be that behavioural inhibition, as a core social anxiety trait, partly underlies the development of OCD and potentially represents an important endophenotype related to the co-occurring nature of OCD and anxiety disorders.

Relationships for OCRD symptoms

With respect to the initial multivariate analysis, specific genetic influences were more apparent with regard to hoarding disorder and BDD symptoms. Considering that no multivariate twin studies have examined them together with anxiety disorder symptoms, these results are novel. In a recent study of five OCRD domains, OCD and hoarding disorder symptoms were characterised as sharing more common liability, followed by BDD, then trichotillomania and skin-picking symptoms. 13 In the independent pathway model with two common genetic factors, only a small percentage of the total genetic effect of BDD and OCD symptoms loaded onto the OCRD latent factor, as compared with hoarding disorder symptoms. The greater loading for hoarding disorder symptoms may be explained by the following two points: (a) hoarding disorder demonstrated a ‘weak’ phenotypic correlation with BDD symptoms – the lowest among all domains assessed (see also Monzani et al 13 ); and (b) hoarding disorder symptoms emerged from the exploratory PCA as a second distinct factor, suggesting it contains some unique variance with regard to the other OCRD domains.

Despite these seemingly complex associations, OCD symptoms were identified as a potential causal risk factor for hoarding disorder symptoms. Previous twin studies have documented a close association between both symptom domains, with hoarding disorder symptoms being reported to share a substantial common genetic liability with other major OCD symptom dimensions. 40 Some authors have suggested that the characteristic feature of indecision in OCD, 41 which is also observed in hoarding disorder, 42,43 may be a significant risk factor for hoarding disorder that is genetically transmitted with OCD. 43 Our results potentially add weight to this hypothesis by demonstrating a putative causal link between these symptom domains. Nonetheless, the unique variance estimated for hoarding disorder symptoms at the multivariate level suggests there are likely to be distinct aetiological factors underlying aspects of hoarding disorder that are not present in OCD, such as the inability to discard.


There are certain limitations to this study. First, all symptoms were assessed by self-report measures including the DCQ, DASS-21 Stress subscale and ASI, which do not perfectly match the diagnostic criteria for BDD, GAD and panic disorder. For example, although the ASI has shown validity in distinguishing between individuals experiencing panic attacks v. those with panic disorder, 44 it is also predictive of other psychiatric disorders. 34 Second, a reliance on self-report limits the generalisation of these findings to dimensional representations of symptoms rather than disorders, and does not allow one to rule out whether symptoms may be as a result of unmeasured third-party factors, such as other mental or medical conditions. Third, it will be important in future multivariate twin studies to assess relationships between OCRDs, anxiety disorders and depression, tic and somatoform disorder symptoms, given ongoing interest in clarifying the aetiological links between these domains. Indeed, it is highly likely that the aetiological common threads suggested by the current results will extend beyond the specific OCRD and anxiety disorder domains studied here. Finally, extension of the current findings in a prospective longitudinal twin study 45 will be important for validating inferences of direct causal associations between these domains.


In conclusion, the current findings suggest that ongoing aetiological (for example molecular genetic) and treatment-focused studies of OCD are likely to benefit from the consideration of a more diverse phenotype that represents its important links with some OCRDs (including tic disorders) but also with certain anxiety disorders. Parallels can be drawn between this sentiment and recent efforts to identify common molecular genetic risk factors that cut across other major psychiatric diagnoses, including schizophrenia, bipolar disorder and major depression, among others. 46,47 Importantly, if confirmed by future studies, the observed causal influences identified here may encourage novel approaches to treatment intervention, with potential to reduce the overall burden of these disorders when co-occurring in individual patients.


1 Hollander, E, Braun, A, Simeon, D. Should OCD leave the anxiety disorders in DSM-V? The case for obsessive compulsive-related disorders. Depress Anxiety 2008; 25: 317–29.
2 Bartz, JA, Hollander, E. Is obsessive-compulsive disorder an anxiety disorder? Prog Neuropsychopharmacol Biol Psychiatry 2006; 30: 338–52.
3 Storch, EA, Abramowitz, J, Goodman, WK. Where does obsessive-compulsive disorder belong in DSM-V? Depress Anxiety 2008; 25: 336–47.
4 Stein, DJ, Fineberg, NA, Bienvenu, OJ, Denys, D, Lochner, C, Nestadt, G, et al. Should OCD be classified as an anxiety disorder in DSM-V? Depress Anxiety 2010; 27: 495506.
5 Bienvenu, OJ, Samuels, JF, Wuyek, LA, Liang, K-Y, Wang, Y, Grados, MA, et al. Is obsessive-compulsive disorder an anxiety disorder, and what, if any, are spectrum conditions? A family study perspective. Psychol Med 2012; 42: 113.
6 American Psychiatric Association. Diagnostic and Statistical Manual of Mental Disorders (5th edn) (DSM-5). APA, 2013.
7 Lochner, C, Stein, DJ. Obsessive-compulsive spectrum disorders in obsessive-compulsive disorder and other anxiety disorders. Psychopathology 2010; 43: 389–96.
8 Hollander, E, Zohar, J. Beyond refractory obsessions and anxiety states: toward remission. J Clin Psychiatry 2004; 65 (suppl 1): 35.
9 Stein, DJ, Lochner, C. Obsessive-compulsive spectrum disorders: a multidimensional approach. Psychiatr Clin North Am 2006; 29: 343–51.
10 Pallanti, S, Hollander, E. Obsessive-compulsive disorder spectrum as a scientific “metaphor”. CNS Spectr 2008; 13: 615.
11 Neale, MC, Maes, HM. Methodology for Genetic Studies of Twins and Families. Kluwer Academic Publishers, 2001.
12 Tambs, K, Czajkowsky, N, Røysamb, E, Neale, MC, Reichborn-Kjennerud, T, Aggen, SH, et al. Structure of genetic and environmental risk factors for dimensional representations of DSM-IV anxiety disorders. Br J Psychiatry 2009; 195: 301–7.
13 Monzani, B, Rijsdijk, F, Harris, J, Mataix-Cols, D. The structure of genetic and environmental risk factors for dimensional representations of DSM-5 obsessive-compulsive spectrum disorders. JAMA Psychiatry 2014; 71: 182–9.
14 Monzani, B, Rijsdijk, F, Iervolino, AC, Anson, M, Cherkas, L, Mataix-Cols, D. Evidence for a genetic overlap between body dysmorphic concerns and obsessive-compulsive symptoms in an adult female community twin sample. Am J Med Genet B Neuropsychiatr Genet 2012; 159B: 376–82.
15 Hopper, JL, Bui, QM, Erbas, B, Matheson, MC, Gurrin, LC, Burgess, JA, et al. Does eczema in infancy cause hay fever, asthma, or both in childhood? Insights from a novel regression model of sibling data. J Allergy Clin Immunol 2012; 130: 1117–22.
16 López-Solà, C, Fontenelle, LF, Alonso, P, Cuadras, D, Foley, D, Pantelis, C, et al. Prevalence and heritability of obsessive-compulsive spectrum and anxiety disorder symptoms: a survey of the Australian Twin Registry. Am J Med Genet B Neuropsychiatr Genet 2014; 165B: 314–25.
17 Foa, EB, Huppert, JD, Leiberg, S, Langner, R, Kichic, R, Hajcak, G, et al. The Obsessive-Complusive Inventory: development and validation of a short version. Psychol Assess 2002; 14: 485–95.
18 Tolin, DF, Frost, RO, Steketee, G. A brief interview for assessing compulsive hoarding: the Hoarding Rating Scale-Interview. Psychiatry Res 2010; 178: 147–52.
19 Oosthuizen, P, Lambert, T, Castle, DJ. Dysmorphic concern: prevalence and associations with clinical variables. Aust N Z J Psychiatry 1998; 32: 129–32.
20 Connor, KM, Davidson, JR, Churchill, LE, Sherwood, A, Foa, E, Weisler, RH. Psychometric properties of the Social Phobia Inventory (SPIN). New self-rating scale. Br J Psychiatry 2000; 176: 379–86.
21 Reiss, S, Peterson, RA, Gursky, DM, McNally, RJ. Anxiety sensitivity, anxiety frequency and the prediction of fearfulness. Behav Res Ther 1986; 24: 18.
22 Antony, MM, Bieling, PJ, Cox, BJ, Enns, MW, Swinson, RD. Psychometric properties of the 42-item and 21-item versions of the Depression Anxiety Stress Scales in clinical groups and a community sample. Psychol Assess 1998; 10: 176–81.
23 Box, G, Cox, D. An Analysis of Transformations. J Roy Stat Soc Ser B 1964; 26: 211–52.
24 Zeegers, M, Rijsdijk, F, Sham, P. Adjusting for covariates in variance components QTL linkage analysis. Behav Genet 2004; 34: 127–33.
25 Neale, MC, Røysamb, E, Jacobson, K. Multivariate genetic analysis of sex limitation and G x E interaction. Twin Res Hum Genet 2006; 9: 481–9.
26 Boker, S, Neale, M, Maes, H, Wilde, M, Spiegel, M, Brick, T, et al. OpenMx: an open source extended structural equation modeling framework. Psychometrika 2011; 76: 306–17.
27 Dite, GS, Gurrin, LC, Byrnes, GB, Stone, J, Gunasekara, A, McCredie, MRE, et al. Predictors of mammographic density: insights gained from a novel regression analysis of a twin study. Cancer Epidemiol Biomarkers Prev 2008; 17: 3474–81.
28 Kirk, KM, Birley, AJ, Statham, DJ, Haddon, B, Lake, RI, Andrews, JG, et al. Anxiety and depression in twin and sib pairs extremely discordant and concordant for neuroticism: prodromus to a linkage study. Twin Res 2000; 3: 299309.
29 Taylor, S. Etiology of obsessions and compulsions: a meta-analysis and narrative review of twin studies. Clin Psychol Rev 2011; 31: 1361–72.
30 Nestadt, G, Samuels, J, Riddle, MA, Liang, KY, Bienvenu, OJ, Hoehn-Saric, R, et al. The relationship between obsessive-compulsive disorder and anxiety and affective disorders: results from the Johns Hopkins OCD Family Study. Psychol Med 2001; 31: 481–7.
31 Brakoulias, V, Starcevic, V, Sammut, P, Berle, D, Milicevic, D, Moses, K, et al. Obsessive-compulsive spectrum disorders: a comorbidity and family history perspective. Australas Psychiatry 2011; 19: 151–5.
32 Pallanti, S, Quercioli, L. Shame and psychopathology. CNS Spectr 2000; 5: 2843.
33 Wheaton, MG, Mahaffey, B, Timpano, KR, Berman, NC, Abramowitz, JS. The relationship between anxiety sensitivity and obsessive-compulsive symptom dimensions. J Behav Ther Exp Psychiatry 2012; 43: 891–6.
34 Calamari, JE, Rector, NA, Woodard, JL, Cohen, RJ, Chik, HM. Anxiety sensitivity and obsessive–compulsive disorder. Assessment 2008; 15: 351–63.
35 Armstrong, T, Zald, DH, Olatunji, BO. Attentional control in OCD and GAD: specificity and associations with core cognitive symptoms. Behav Res Ther 2011; 49: 756–62.
36 Stasik, SM, Naragon-Gainey, K, Chmielewski, M, Watson, D. Core OCD symptoms: exploration of specificity and relations with psychopathology. J Anxiety Disord 2012; 26: 859–70.
37 Rosenbaum, JC, Biederman, J, Faraone, SV, Hirshfeld, DR, Kagan, J. Behavioral inhibition in childhood: a risk factor for anxiety disorders. Biol Psychiatry 1997; 42: 145S.
38 Coles, ME, Schofield, CA, Pietrefesa, AS. Behavioral inhibition and obsessive-compulsive disorder. J Anxiety Disord 2006; 20: 1118–32.
39 Eley, TC, Bolton, D, O'Connor, TG, Perrin, S, Smith, P, Plomin, R. A twin study of anxiety-related behaviours in pre-school children. J Child Psychol Psychiatry 2003; 44: 945–60.
40 Iervolino, AC, Rijsdijk, FV, Cherkas, L, Fullana, MA, Mataix-Cols, D. A multivariate twin study of obsessive-compulsive symptom dimensions. Arch Gen Psychiatry 2011; 68: 637–44.
41 Sachdev, PS, Malhi, GS. Obsessive-compulsive behaviour: a disorder of decision-making. Aust N Z J Psychiatry 2005; 39: 757–63.
42 Frost, RO, Gross, RC. The hoarding of possessions. Behav Res Ther 1993; 31: 367–81.
43 Samuels, JF, Bienvenu, OJ, Pinto, A, Fyer, AJ, McCracken, JT, Rauch, SL, et al. Hoarding in obsessive-compulsive disorder: results from the OCD Collaborative Genetics Study. Behav Res Ther 2007; 45: 673–86.
44 Rector, NA, Szacun-Shimizu, K, Leybman, M. Anxiety sensitivity within the anxiety disorders: disorder-specific sensitivities and depression comorbidity. Behav Res Ther 2007; 45: 1967–75.
45 Stone, J, Dite, GS, Giles, GG, Cawson, J, English, DR, Hopper, JL. Inference about causation from examination of familial confounding: application to longitudinal twin data on mammographic density measures that predict breast cancer risk. Cancer Epidemiol Biomarkers Prev 2012; 21: 1149–55.
46 Lee, SH, Ripke, S, Neale, BM, Faraone, SV, Purcell, SM, Perlis, RH, et al. Genetic relationship between five psychiatric disorders estimated from genome-wide SNPs. Nat Genet 2013; 45: 984–94.
47 Smoller, JW, Ripke, S, Lee, P-H, Neale, B, Nurnberger, JI, Santangelo, S, et al. Identification of risk loci with shared effects on five major psychiatric disorders: a genome-wide analysis. Lancet 2013; 381: 1371–9.