The British Journal of Psychiatry

Abstract

Background Symptomatology in psychoses can be summarised as quantitative symptom dimensions, but their genetic basis is unknown.

Aims To investigate whether genes make an important contribution to symptom dimensions.

Method A total of 224 probandwise twin pairs (106 monozygotic, 118 same-gender dizygotic) where probands had psychosis were ascertained from the Maudsley Twin Register in London. Factor analysis was performed on lifetime symptoms rated on the Operational Checklist for Psychotic Disorders (OPCRIT). Correlations of dimension scores within monozygotic and dizygotic pairs concordant for Research Diagnostic Criteria psychoses were performed. Relationships between dimension scores and genetic loading for psychoses were assessed using logistic regression.

Results Patterns of familial aggregation consistent with a genetic effect were found for the disorganised dimension and for some measures of the negative, manic and general psychotic dimensions. Disorganised dimension scores were related significantly to genetic loading for psychoses.

Conclusions The disorganised dimension, and possibly other symptom dimensions, may be useful phenotypes for molecular genetic studies of psychoses.

It is well established that symptomatology in psychotic disorders can be ascribed by factor analysis to three main psychotic symptom dimensions (positive, negative and disorganised; Liddle, 1987; Andreasen et al, 1995) and two main mood dimensions (manic and depressive; Lindenmayer et al, 1994; van Os et al, 1997). In addition, studies that include a wide range of positive symptoms often find subdivisions of the positive dimension, e.g. first-rank delusions and paranoid dimensions (Cardno et al, 1997). It is not known whether genetic factors make an important contribution to variation in these symptom dimensions. Studies of affected sib-pairs show modest familial aggregation (Burke et al, 1996; Kendler et al, 1997; Loftus et al, 1998; Cardno et al, 1999a) and familial risk of psychoses has been associated with proband scores on the negative, disorganised, manic and depressive dimensions (Cardno et al, 1997; van Os et al, 1997). In this study, we extended such investigations to a sample of twins with a range of psychoses: we investigated the degree of familial aggregation of dimensions in monozygotic (MZ) and dizygotic (DZ) twin pairs concordant for psychoses; and we investigated whether dimensions are markers of the degree of genetic loading (or liability) for psychoses, by determining whether dimension scores in ill probands predict the presence of psychoses in their co-twins, and whether this effect is greater in MZ than DZ pairs. If the genes that influence a symptom dimension do not affect the risk of developing psychoses but only influence the clinical presentation (i.e. are modifying genes), there should be greater familial aggregation for the dimension in MZ than in DZ concordant pairs, but dimension scores in ill twin probands should not predict the presence of psychoses in their co-twins. However, if some genes that affect the risk of developing psychoses also influence what symptoms occur, relevant symptom dimension scores in probands should predict the presence of psychoses in co-twins.

METHOD

Subjects

The sample has been described in detail elsewhere (Cardno et al, 1999b). Briefly, probands were ascertained from the Maudsley Twin Register in London as patients of multiple birth who had attended the Maudsley and Bethlem Royal Hospitals between 1948 and 1993, had a same-gender co-twin surviving to 15 years of age and had suffered psychotic symptomatology (Sartorius et al, 1986) or an episode of Research Diagnostic Criteria (RDC) mania or hypomania (Spitzer et al, 1978) at some time in their lives. Thus, a broad concept of psychoses was adopted. Probands were excluded if they had dementia or if psychotic symptoms occurred only during acute organic states.

Zygosity determination was carried out blind to research diagnosis and was based on all available information, including analysis of genetic markers in 42% of pairs. Agreement between zygosity determination by genetic markers and by resemblance information was 95%.

The sample comprised 224 probandwise pairs (120 male and 104 female) from which 28 probands (12.5%) were doubly ascertained; 197 (87.9%) pairs were White. The mean age of the co-twins at last followup was 46.5 years (s.d.=15.4). There was no substantial ascertainment bias for zygosity (106 MZ and 118 same-gender DZ pairs), but the sample may have contained more severe cases than if we had ascertained via a population-based sample (Cardno et al, 1999b).

Clinical assessment

This was based on all available clinical information, including a research interview (Gottesman & Shields' (1972) cued-questions interview) and/or the Schedule for Affective Disorders and Schizophrenia — Lifetime version (SADS—L; Spitzer & Endicott, 1978), in 72.8% of probands and 59.4% of co-twins. Interviewed twins gave informed consent to be studied. The RDC lifetime-ever diagnoses were made by separate raters for each member of the concordant twin pairs. A non-hierarchical diagnostic approach was employed, that is, twins could qualify for more than one psychotic diagnosis on a lifetime basis.

Lifetime psychotic and affective symptoms were rated as present/absent using the Operational Criteria Checklist for Psychotic Disorders (OPCRIT; McGuffin et al, 1991). Interrater reliability was assessed on 30 cases and expressed as the mean kappa coefficient (κ) between raters for symptoms and the mean Spearman correlation (ρ) for symptom dimension scores. Variables with a reliability coefficient of less than 0.41, which is the lower limit of moderate agreement suggested by Landis & Koch (1977), were excluded from further analysis.

Previous studies have defined dimensions as either symptom scores (Loftus et al, 1998) or factor scores (Burke et al, 1996; Kendler et al, 1997). Because it was not clear which approach would be most useful, both were employed.

Dimensions based on symptom scores were derived from the summation of OPCRIT symptoms that were characteristic of dimensions in previous studies. The following dimensions were included (with their score range, symptoms included and reliability correlations): positive (0-2; delusions, hallucinations; ρ=0.69); negative (0-2; restricted affect, negative formal thought disorder; ρ=0.65); disorganised (0-2; positive formal thought disorder, inappropriate affect; ρ=0.66); manic (0-5; mood change, enhanced abilities, decreased need for sleep, speech/thought change, disinhibition; ρ =0.75); depressive (0-5; mood change, appetite change, sleep change, negative cognition/suicidal, agitation/slowness/energy loss; ρ=0.66); general psychotic (0-6; sum of positive, negative and disorganised dimensions).

Statistical analysis

Dimensions based on factor scores were derived from two factor analyses of OPCRIT symptoms. The first analysis followed the procedure adopted in our previous factor analyses of OPCRIT symptoms (Cardno et al, 1999a). The 26 individual psychotic symptom items were considered initially in 256 twins rated as having at least one of these symptoms, after excluding twins with only partial clinical information. The symptoms of incoherent speech, blunted affect, delusional perception, thought withdrawal, nihilistic delusions and thought echo were excluded because they were found in fewer than 10% of twins, and bizarre delusions and other primary delusions were excluded because they had insufficient reliability. This left 18 symptoms for further analysis. A principal components analysis was performed using SPSS Windows Version 6.0 software. The number of substantive factors was determined by the scree test (Cattell, 1966), and these underwent orthogonal (varimax) rotation. Dimension scores then were derived as regression factor scores.

The second factor analysis followed a similar approach but was based on a combination of psychotic and affective symptoms without including an excessively large number of symptoms for analysis. The six psychotic, five manic and five depressive items used to compute the symptom scores were chosen, and all 16 symptoms were included in the analysis, which involved 291 twins.

Spearman correlations were used to assess familial aggregation in twin pairs concordant for psychoses. Correlations were performed for MZ and DZ pairs concordant for any psychosis and for MZ pairs concordant for schizophrenia. Too few pairs (< 10) were concordant for other diagnoses to allow a meaningful analysis.

Logistic regression analysis was used to analyse associations between symptom dimensions and either familial or genetic loading for psychoses. Pairs where the proband had a particular psychotic diagnosis were selected. Presence/absence of the same diagnosis in co-twins was treated as the dependent variable and the particular dimension as the independent variable, with zygosity and age of co-twins at follow-up as covariates. An effect of a dimension at P < 0.05 was taken as evidence of an association with familial loading for the psychosis. In such cases a zygosity-by-dimension interaction term was added, and an effect of this term at P < 0.05 was taken as evidence of an association between the dimension and genetic loading for the psychosis (i.e. a greater effect in MZ than DZ pairs). Effects were investigated for lifetime-ever RDC schizophrenia, schizoaffective disorder, mania and all RDC psychoses combined.

RESULTS

Factor analyses

In the first factor analysis (of psychotic symptoms), the scree test suggested six substantive factors: disorganised, negative, first-rank delusions, paranoid, other hallucinations and first-rank hallucinations (Table 1). In the second factor analysis (of both psychotic and affective symptoms) there were three substantive factors: manic, general psychotic and depressive (Table 2).

View this table:
Table 1

First factor analysis: OPCRIT psychotic symptoms (loadings of > 0.50 in bold)

View this table:
Table 2

Second factor analysis: OPCRIT psychotic and affective symptoms (loadings of >0.50 in bold)

Within-pair correlations of symptom dimensions

These are shown in Table 3. In pairs concordant for any psychosis, there were trends towards higher MZ than DZ correlations for negative, disorganised, manic and general psychotic symptom scores, and relatively substantial MZ correlations for these dimensions based on factor scores. An atypical feature of the negative factor was the high loading for inappropriate affect; therefore, the first factor analysis was repeated excluding this symptom. Six similar factors emerged, with the highest loadings on the negative factor for negative formal thought disorder (0.75), catatonia (0.61) and restricted affect (0.45). Within-pair correlations of negative factor scores remained similar (for MZ pairs with any psychosis, ρ=0.52).

View this table:
Table 3

Within-pair Spearman correlations of symptom dimensions in concordant twin pairs

Relationships with genetic loading for psychoses

The odds ratios (OR) from analysis of relationships between the symptom dimensions and familial loading for psychoses are shown in Table 4. For all psychoses, higher familial loading was associated significantly with higher manic symptom and factor scores and higher disorganised factor scores. In addition, higher paranoid factor scores were associated with lower familial loading for mania. When each of the dimensions showing a significant association with familial loading for all psychoses was included in an interaction term with zygosity (i.e. investigating an association with genetic loading for psychoses), the only significant result was for disorganised factor scores (OR=1.47; P=0.03). The odds ratio for the interaction term remained similar after adjusting for gender (OR=1.46; P=0.04) and age of onset (OR=1.42; P=0.053). Higher disorganised scores are associated with earlier age of onset (ρ=-0.17; P=0.01), so the two variables are not entirely independent, but disorganised scores are better markers of genetic loading for psychoses because age of onset is not associated significantly with genetic loading for psychoses in this sample (interaction OR=0.97; P=0.44).

View this table:
Table 4

Logistic regression of RDC lifetime-ever psychotic diagnoses in co-twins on symptom dimensions in probands, conditional on zygosity and age of co-twins at last follow-up

DISCUSSION

Familial aggregation of symptom dimensions

Positive dimension

The results of this study are not consistent with an important genetic contribution to variation in scores on the positive dimension or its subdivisions within psychoses. Where zygosities could be compared directly, MZ correlations were no higher than DZ correlations. Furthermore, the MZ correlations in general were similar to correlations from studies of affected sib-pairs, including those for factor (up to 0.21) and symptom scores (0.21) derived from OPCRIT in the same way as for the current study (Cardno et al, 1999a).

Comparison with the results of the sib-pair studies is of interest but should be viewed with caution because of methodological differences from the current study. The closest point of comparison is for schizophrenia because the sib-pair studies have been based on schizophrenia plus related disorders in some cases. Although various definitions of schizophrenia have been employed in the sib-pair studies, this is unlikely to affect substantially the comparison with the current study because all of the concordant MZ twins fulfilled criteria for RDC, DSM—III—R (American Psychiatric Association, 1987) and ICD—10 (World Health Organization, 1993) schizophrenia.

Negative dimension

The correlations for symptom scores in MZ pairs with schizophrenia and DZ pairs with any psychosis were close to zero. This was unexpectedly lower than the correlations found in studies of sib-pairs (range 0.09-0.35) and may be an artefact of the small numbers analysed. The correlation for symptom scores in MZ pairs with any psychosis was modest. However, the correlations for factor scores in MZ pairs with schizophrenia and any psychosis were more substantial, and greater than twice the correlations for sib-pairs using OPCRIT (0.16; Cardno et al, 1999a). The discrepancy between the correlations for symptom and factor scores may be due to the larger number of symptoms and symptom weightings contributing to the factor scores. It is not accounted for by the atypically high loading for inappropriate affect on the negative factor. These results are consistent with a genetic contribution to variation in the negative dimension, at least as defined by factor scores.

Disorganised dimension

Correlations of both symptom and factor scores for MZ pairs with schizophrenia and any psychosis were substantial and at least twice as great as the correlations for sib-pairs using OPCRIT (symptom scores=0.25; factor scores=0.32; Cardno et al, 1999a). The symptom score correlation for DZ pairs with any psychosis was less than half that of the MZ correlation and similar to the OPCRIT-derived sib-pair correlation for schizophrenia/schizoaffective disorder (Cardno et al, 1999a). This is consistent with a genetic contribution to variation in the disorganised dimension.

Manic dimension

There were markedly negative correlations for MZ pairs with schizophrenia. These were unexpectedly lower than sib-pair correlations (factor score=0.27 and symptom score=0.43; Kendler et al, 1997) and may be an artefact of the small numbers analysed. The same applies to the low correlations for DZ pairs with any psychosis. In contrast, the correlations for MZ pairs with any psychosis were more substantial and approximately 1.5-fold higher than the sib-pair correlation for all non-affective psychoses (0.27 and 0.26; Kendler et al, 1997). These results are consistent with a genetic contribution to variation in manic dimension scores within psychoses.

Depressive dimension

Correlations were uniformly low. In some cases this could have been an artefact of small samples, but even for MZ pairs with any psychosis the correlation was marginally lower than for sib-pairs with any non-affective psychosis (0.28; Kendler et al, 1997). These results do not suggest that genetic effects make an important contribution to variation in depressive dimension scores within psychoses. Although major depressive disorder has a heritability of 40% or higher (Kendler & Prescott, 1999), there may be differences between this disorder in general and variation of depressive symptomatology in psychoses.

General psychotic dimension

Correlations were low in MZ twins with schizophrenia and similar to the correlation in our sib-pair sample (0.12; details available from P.M. upon request). This could be because there is little genetic contribution to variation in this dimension in schizophrenia, but equally it could be an artefact of the small numbers in this part of the sample. In contrast, for all psychoses the MZ correlations were consistently more substantial (0.47 and 0.53) and considerably greater than the corresponding DZ correlations (0.27 and 0.17). Some caution again is required in interpreting the DZ correlations because of the small numbers, but the pattern is consistent with a genetic contribution to variation in the general psychotic dimension in individuals with any psychosis.

Relationships with genetic loading for psychoses

Only the disorganised dimension, defined by factor scores, showed a significant association with genetic loading for all psychoses combined. This is partly consistent with previous family history studies. Disorganised scores in probands with RDC schizophrenia have been associated with familial risk of psychosis (Cardno et al, 1997; van Os et al, 1997), but in one study (van Os et al, 1997) negative dimension scores in probands with any psychosis best predicted the familial risk of psychosis. It is also notable that there are similarities between the disorganised dimension and the hebephrenic schizophrenia subtype, and hebephrenia has been associated with relatively high familial loading for schizophrenia in previous studies of Maudsley and other twin samples (Farmer et al, 1984; McGuffin et al, 1987; Onstad et al, 1991).

The disorganised dimension, defined by symptom scores, showed a trend towards an association with familial loading for all psychoses that failed to reach statistical significance (OR=1.23; P=0.32). This may be because the 0-2 symptom scale resulted in less statistical power than the continuous factor scores. The pattern of symptom weightings that contributed to disorganised factor scores also potentially could account for the difference. In addition to positive formal thought disorder, the disorganised factor included a high loading for grandiose delusions. This symptom has not been part of the disorganised dimension in most previous factor analytical studies, where only global delusion scores have been considered (Liddle, 1987; Andreasen et al, 1995), but has had moderate or high loadings in previous studies using OPCRIT in subjects with schizophrenia or schizo-affective disorder (Cardno et al, 1999a) and using the Scales for the Assessment of Positive and Negative Symptoms (SAPS/SANS; Andreasen, 1984a,b) for individual symptoms in subjects with a range of psychoses (Toomey et al, 1997).

The presence of grandiose delusions in the disorganised dimension raises the question of the relationship between disorganised and manic symptoms. Disorganised and manic factor scores were correlated significantly even when grandiose delusions were omitted from the disorganised dimension (ρ=0.29; P<0.0005), which is consistent with previous findings (Johnstone & Frith, 1996). The apparent relationship between disorganised and manic dimensions could be due to mislabelling of symptoms, because it can be difficult to distinguish reliably the knight's-move thinking from a flight of ideas (Andreasen, 1979). It is possible also that the two dimensions comprise distinct symptoms that are associated with common pathological processes; for example, disorganisation has been associated with cognitive disinhibition (Liddle & Morris, 1991) and disinhibition is an important clinical feature of mania.

The lack of relationship between the positive dimension and its subdivisions with familial loading for psychoses is consistent with previous studies (Kendler et al, 1994; Cardno et al, 1997; van Os et al, 1997). The lack of effect for the general psychotic dimension provides some support for the approach of dividing clinical variation within psychoses into a number of symptom dimensions. The lack of effect for the depressive dimension for all psychoses is consistent with the results of van Os et al (1997), who found an association only in schizoaffective disorder. The lack of effect for the negative dimension is consistent with some studies (Kendler et al, 1994; Cardno et al, 1997) but not others (Dworkin & Lenzenweger, 1984; van Os et al, 1997). One of the positive studies found the effect for clinical schizophrenia in MZ twins (Dworkin & Lenzenweger, 1984). In the current study there was a trend towards the largest (non-significant) odds ratio with schizophrenia (1.21), so a small effect cannot be excluded. The other family history study (van Os et al, 1997) found the effect for all psychoses using OPCRIT, but the dimensions differed from those in the current study in that they included illness history variables as well as symptoms. The current negative factor included inappropriate affect, which is more commonly associated with disorganisation (and indeed had a moderate loading on this dimension). However, we performed a further analysis of negative factor scores excluding the inappropriate affect, and the results for the relationship between the negative dimension and familial liability to all psychoses (OR=1.07; P=0.68) and schizophrenia (OR=1.24; P=0.38) remained similar. Therefore, the relationship between the negative dimension and familial loading for psychoses remains inconclusive.

Limitations

As noted above, the twin sample is probably selected for severity, and the small numbers in some categories limited the number of analyses that could be performed reasonably and limited direct MZ—DZ comparison with all psychoses combined. There is also a need for caution when comparing the results of the current study with the results of sib-pair studies: there are general methodological differences (e.g. concerning ascertainment); the twin sample includes same-gender pairs only; and there are potential twin versus sibling differences regarding familial aggregation of dimensions. In addition, comparison of factor scores across studies should be done cautiously because of variation in the patterns of factor loadings in different samples. Different results also could have occurred if dimensions were based on other rating scales. However, factor scores derived from OPCRIT and SAPS/SANS in our sib-pair study were correlated strongly (Cardno et al, 1999a).

Because the data were either ordinal or not easily transformed to a normal distribution, we employed non-parametric Spearman correlations. Confidence intervals therefore could not be calculated. It would be possible technically to apply biometric model fitting to make heritability estimates for the dimensions, but the sample size was not sufficient for this to be done with confidence.

The MZ—DZ differences in familial aggregation of dimensions will tend to be underestimated if genetic risk factors for psychotic disorders contribute to the variation in symptom dimensions. This is because pairs of DZ twins who both have psychosis will tend to share susceptibility alleles more frequently than the expected 50%, and so will be more like MZ twins at loci relevant to the aetiology of psychoses. However, the intrapair correlations in affected sib-pairs are modest, suggesting that genetic effects are due predominantly to genes with modifying effects on illness expression, with susceptibility genes themselves probably having a more minor role.

In the logistic regression analysis, numerous statistical tests were performed and so the result for disorganisation could be a type I error. However, the support from previous studies (Cardno et al, 1997; van Os et al, 1997) argues against this. It seems more likely that some of the non-significant findings in the current study were type II errors. For example, the association between manic factor scores and familial loading for mania was non-significant for an odds ratio of 3.54.

Future genetic studies that include cognitive, physiological and brain imaging data may shed further light on the genetic basis of symptom dimensions and their relationship with both pathological processes in psychotic disorders and extremes of variation in the general population.

Clinical Implications and Limitations

CLINICAL IMPLICATIONS

  • There is consistent evidence for a genetic contribution to variation in the disorganised symptom dimension, and to a lesser extent to the negative, manic and general psychotic dimensions.

  • Individuals with high disorganised dimension scores have relatively high genetic loading for psychotic disorders.

  • The disorganised dimension therefore may be a useful phenotype for molecular genetic studies of psychotic disorders.

LIMITATIONS

  • Although this is the largest twin study of operationally defined psychotic disorders to date, numbers of concordant pairs in some categories were small.

  • Dimensions based on other rating scales could produce different results.

  • The sample may have contained more severe cases than if it had been ascertained via a population-based sample.

Acknowledgments

This work was carried out while one of the authors (A.G.C.) was a Medical Research Council (UK) Clinical Training Fellow. The National Institute of Mental Health (MH44359) and the Stanley Foundation provided grant support that enabled collection of the clinical data. The authors thank Jane Marshall, Bina Coid, Alison Macdonald, Tracy Ribchester, Nadia Davies, Piero Venturi, Lisa Jones, Shôn Lewis, Irving Gottesman, Anne Farmer and Adrianne Reveley for their work on this study.

  • Received September 21, 2000.
  • Revision received February 13, 2001.
  • Accepted February 15, 2001.

References

View Abstract