The inheritance of migraine with aura estimated by means of structural equation modelling
- aDepartment of Neurology, Glostrup Hospital, University of Copenhagen, DK-2600 Glostrup, Copenhagen, Denmark, bThe Danish Twin Register, Department of Genetic Epidemiology, Institute of Community Health, Odense University, DK-5000 Odense, Denmark
- Dr Ulrich.
- Received 23 January 1998
- Revised 10 August 1998
Studies of migraine with aura (MA) have shown familial aggregation of the disorder, which cannot be explained by simple mendelian inheritance. The interest in a genetic basis for the disorder has increased after identification of three genetic loci for familial hemiplegic migraine (FHM), which is a rare subtype of MA with autosomal dominant inheritance. Both genetic and environmental factors seem to be important in the expression of MA. To elucidate the molecular pathogenesis of MA, knowledge of the relative role of genetic and environmental factors is essential. Twin studies are a classic way to analyse this. We applied structural equation modelling on MA with twin data obtained from a population based twin register in order to evaluate the effects of genes and environment. The correlation in liability of MA was 0.68 in monozygotic (MZ) and 0.22 in dizygotic (DZ) twin pairs, indicating a high degree of genetic determination in the total variance of liability. The best fitting model combined additive genetic effects and environmental effects that were not shared by the twins. The estimate of heritability was 0.65 and similar in males and females.
Migraine with aura (MA, formerly classical migraine) is a common neurological disorder with age at onset before 40 years in the majority of cases.1 2 Attacks are characterised by an initial aura, that is, reversible visual, sensory, aphasic, or motor symptoms usually lasting less than 60 minutes, and followed by a headache.3 The attacks may be accompanied by nausea, photophobia, and phonophobia.3 Familial hemiplegic migraine (FHM) is a rare subtype of MA with autosomal dominant inheritance. The discovery that FHM is a genetically heterogeneous disorder with at least three different loci4-7 has stimulated the search for a genetic basis of MA.
Studies of families with MA have suggested a familial susceptibility in the aetiology,8 9 which does not follow simple mendelian inheritance. An analysis of 31 high risk families with affected subjects in more than one generation found that autosomal dominant inheritance with reduced penetrance was unlikely even in these families.10 Autosomal recessive inheritance is unlikely because of the unequal sex distribution, females being more often affected than males.1 2 11 However, the female preponderance is too low to suggest sex influenced inheritance.10 Mitochondrial and X linked inheritance may occur in subgroups of affected subjects, but are excluded in many families because of paternal transmission.
MA is likely to be caused by interactions of genetic and environmental factors. A complex segregation analysis of 127 probands and their first degree relatives from the general population indicated multifactorial inheritance without generational differences.12 Studies of twins may help to clarify the relative importance of genetic and environmental factors. A recent twin study has emphasised the importance of genetic factors, since monozygotic (MZ) twin pairs had significantly higher concordance than dizygotic (DZ) twin pairs.13
We applied structural equation modelling on MA with twin data obtained from a population based twin register, thereby evaluating the influence of genetic and environmental factors and providing a quantitative estimate of the heritability.
Materials and methods
The study population originated from the population based Danish Twin Register which is representative of Danish twins.14Twins born between 1953 and 1960 were included and comprised 2026 MZ and 3334 same sex DZ twins. A total of 87% of the twins completed a screening questionnaire about migraine.13 The level of ascertainment was 85%.13 All twin pairs with at least one twin who had self-reported migraine or self-reported severe headache with accompanying symptoms were interviewed by one of two neurological residents experienced in headache diagnosis (VU, MG). The participation rate in the telephone interview was 90%.13 The criteria of the International Headache Society were used.3 A more detailed description of the study design and methods has been given elsewhere.13
The structural equation modelling approach is based on a classical biometric analysis of data. The underlying assumption is that the variance of a quantitative phenotype can be broken down into two components, the genotypic variance and the environmental variance.15 The genotypic variance can be subdivided into the variance owing to additive genetic effects (A) and the variance owing to non-additive genetic effects (D).15 The additive genetic effects are the effects of genes taken singly and added over multiple loci. The non-additive genetic effects are the effects of intralocus gene interactions. Correspondingly, the environmental variance is divided into effects shared by subjects (C) and effects not shared by subjects (E).15 The shared environmental effects contribute to the subjects’ phenotypic similarity and the non-shared environmental effects contribute to the subjects’ dissimilarity. The model fit is assessed by combining the various parameters A, C, D, and E. Possible models are ACE, ADE, AE, CE, and a model with E as the single parameter. The effects of additive genetic effects A and shared environmental effects C cannot be combined,16 leaving the ACE model a theoretical possibility. The criteria for best fitting model was Akaike’s information criterion (AIC) which combines the goodness of fit χ2 with degrees of freedom.16 The model with the lowest value of AIC is considered to have the best fit.16
Heritability in the broad sense is defined as the proportion of the total variance owing to the genetic variance and is a quantitative measure.15 To estimate the different components of the total variance for a non-quantitative trait like MA by means of structural equation modelling, it is assumed that there is an underlying normally distributed liability to MA. The concept of liability implies a graded continuum of unmeasured, continuously distributed latent traits.15 17 A high degree of correlation of liability between MZ twins expresses a high degree of genetic determination of the total variance in liability to a disease. The correlations of liability are calculated for MZ and DZ twins by means of the probandwise concordance rates and the lifetime prevalence of MA. The probandwise concordance rate is the proportion of affected twin partners of probands in relation to the total number of affected twins.16
The structural equation modelling was performed by means of the MX software computer program.
A total of 211 twin pairs (77 MZ, 134 DZ) were identified.13 The number of concordant pairs was 42 (26 MZ, 16 DZ) and the number of discordant pairs was 169 (51 MZ, 118 DZ).13 The probandwise concordance rates, the correlations in liability, and the lifetime prevalence of MA are shown in table 1. The results of the structural equation modelling are shown in table 2. The best fitting model was an AE model including additive genetic effects and environmental factors that were not shared by the twins. This was also the best fitting model for males analysed separately. In females, both an AE model with additive genetic effects and non-shared environmental effects, and an ADE model with dominant genetic effects and non-shared environmental effects, fitted the data. The heritability estimate was 0.65 under the AE model.
Our results can be considered representative since the ascertainment of twins was optimal.13 14 Only twins born between 1953 and 1960 were included since MA usually occurs before the age of 40 years.1 2 The probandwise concordance rate was significantly higher in MZ than in DZ twin pairs which emphasises that genetic factors are indeed important in the aetiology of MA. The correlation of liability was 0.68 in MZ twin pairs with no significant difference between males and females. This indicates a high degree of genetic determination of MA in the total variance in liability. The correlation of liability was higher in MZ twin pairs than in DZ twin pairs, as expected.
The best fitting model was the AE model which combines additive genetic effects and effects of environment that are not shared by the twins. The ADE model including both dominant genetic and additive genetic effects in combination with non-shared effects of environment also fitted the data in females, while in males the AE model had a better fit than the ADE model. In females the ADE model included dominant genetic effects only. However, a model including dominant genetic effects tends to have additive genetic effects.15Furthermore, the estimate of the non-additive genetic effects had very wide confidence intervals, which makes it difficult to separate the effects. The ADE model is thus unlikely in females and only suitable from a theoretical point of view. The models totally excluding genetic effects fitted the data poorly. In the ACE model including additive genetic effects and effects of environment, only the environmental effects that were not shared by the twins counted, since additive genetic effects and shared environmental effects are impossible to separate.16 The proportion of the non-shared effects in the ACE model was similar to the proportion of non-shared effects of environment in the ADE and AE models. Thus, the total effect of the environmental variance was constant and accounted for approximately 0.35 of the total variance, and included only non-shared effects. Correspondingly, the total effect of the genetic variance, that is, the heritability, is considered constant amounting to approximately 0.65 of the total variance, although the proportion of additive genetic and dominant genetic effects may vary. Our heritability estimate is in accordance with the estimate of 0.79 in a complex segregation analysis of 127 families with MA from the general population.12
Our results indicate that the liability to MA has a high degree of genetic determination and that the phenotypic variation is a combination of genetic and individual specific environmental factors.