Learn more: PMC Disclaimer | PMC Copyright Notice
Genetic Heritability and Shared Environmental Factors Among Twin Pairs With Autism
Abstract
Context
Autism is considered the most heritable of neurodevelopmental disorders, mainly because of the large difference in concordance rates between monozygotic and dizygotic twins.
Objective
To provide rigorous quantitative estimates of genetic heritability of autism and the effects of shared environment.
Design, Setting, and Participants
Twin pairs with at least 1 twin with an autism spectrum disorder (ASD) born between 1987 and 2004were identified through the California Department of Developmental Services.
Main Outcome Measures
Structured diagnostic assessments (Autism Diagnostic Interview–Revised and Autism Diagnostic Observation Schedule) were completed on 192 twin pairs. Concordance rates were calculated and parametric models were fitted for 2 definitions, 1 narrow (strict autism) and 1 broad (ASD).
Results
For strict autism, probandwise concordance for male twins was 0.58 for 40 monozygotic pairs (95% confidence interval [CI], 0.42–0.74) and 0.21 for 31 dizygotic pairs (95% CI, 0.09–0.43); for female twins, the concordance was 0.60 for 7 monozygotic pairs (95% CI, 0.28–0.90) and 0.27 for 10 dizygotic pairs (95% CI, 0.09–0.69). For ASD, the probandwise concordance for male twins was 0.77 for 45 monozygotic pairs (95% CI, 0.65–0.86) and 0.31 for 45 dizygotic pairs (95% CI, 0.16–0.46); for female twins, the concordance was 0.50 for 9 monozygotic pairs (95% CI, 0.16–0.84) and 0.36 for 13 dizygotic pairs (95% CI, 0.11–0.60). A large proportion of the variance in liability can be explained by shared environmental factors (55%; 95% CI, 9%–81% for autism and 58%; 95% CI, 30%–80% for ASD) in addition to moderate genetic heritability (37%; 95% CI, 8%–84% for autism and 38%; 95% CI, 14%–67% for ASD).
Conclusion
Susceptibility to ASD has moderate genetic heritability and a substantial shared twin environmental component.
Autism is a complex neurodevelopmental disorder that interferes with the normal course of social, communicative, and cognitive development. The diagnosis of autism is made in early childhood, with symptoms manifesting within the first 3 years of life. Over the last 3 decades, a substantial increase in the prevalence of autism has been reported, from 4 to 5 per 10 000 in the 1960s to around 40 per 10 000 children1 today. For autism spectrum disorders (ASDs), prevalence rates have been estimated to be about 1%. This 10-fold increase in prevalence has received considerable attention, with concomitant increases in research funding.2
There have been numerous studies designed to characterize the contribution of genetic factors to the development of this disorder. Concordance rates in siblings range from about 3% to 14%3–5 and linkage studies are consistent with a polygenic mode of transmission.6 Three studies7–9 of twins ascertained from clinical samples with a total of 36 monozygotic pairs (concordance rate of 72%) and 30 dizygotic pairs (concordance rate of 0%) have estimated the heritability of autism, or proportion of liability attributable to genetic factors, at about 90%.7 The dizygotic concordance (0%) is substantially lower than the expected rate based on estimates of sibling recurrence rates. More recent twin studies of autism based on either case vignettes10 or telephone interviews11 have also yielded high heritability estimates. In addition, a number of twin studies that used dimensional measures of some core features of autism, such as social responsiveness, have yielded a broad range of heritability estimates of autism.12 However, none of these more recent studies included structured clinical assessments by both parental interview and direct child observation, which is the contemporary standard for establishing the diagnosis of autism or ASD.13
This report describes the results of a twin study of autism: the California Autism Twins Study. The chief aims were to (1) collect a sample of twins with validated diagnoses of autism (ASD) from a population-based sample; (2) examine sex-specific concordance rates for narrow and broad definitions of autism (ASD); and (3) determine the extent to which genetic and shared environmental factors underlie susceptibility.
METHODS
PARTICIPANTS
Twin pairs were identified from the records of the Department of Developmental Services (DDS). The DDS operates a system of 21 regional centers throughout California that coordinate services for persons with autism, mental retardation, and other developmental disabilities. Referral sources for regional centers include primary care providers, educators, public health clinics, other service agencies, and parents. The electronic DDS client files were linked by California Center for Autism and Developmental Disabilities Research and Epidemiology (CADDRE) staff to California live birth records to identify twin births and obtain demographic data. The birth years of the twin pairs included in our study were 1987 to 2004. Twin pairs were considered eligible for the study if they met the following criteria:
- At least 1 child in the pair had a qualifying diagnosis (see later) in the DDS electronic file or in client records; each twin in a pair who met this criterion was considered a proband.
- Maternal residence in California at the time of delivery.
- Both co-twins were alive and residing in California at the time of enrollment.
- No history among proband(s) of neurogenetic conditions that might account for autism (eg, fragile X syndrome, Down syndrome, tuberous sclerosis, and neurofibromatosis).
- Proband(s) had a mental age more than 18 months.
- At least 1 parent with sufficient verbal and reading ability in English or Spanish was able to respond to interviews and checklists.
Using DDS electronic files, CADDRE staff ascertained twins receiving services for autism (DDS autism codes A1 or A2) and twins with a possible ASD based on other recorded information indicating a suspected ASD (DDS autism code A9), DSM-IV diagnosis of ASD (code A4, derived by CADDRE), mental retardation of unknown etiology (MRnoET, derived by CADDRE), or other developmental disability (other, derived by CADDRE). Trained CADDRE staff then abstracted client records to record presence or absence of an ASD diagnosis, criteria used for diagnosis, specific autistic behaviors, and other clinical characteristics. Expert clinical review was then conducted to verify that the children met eligibility based on qualifying diagnosis and other criteria (mentioned earlier). The diagnosis of autism can be most reliably established if a child is at least 4 years or older, so the birth years 1987 to 2004 were included.
A letter from the local regional center and a recruitment package were sent to parents explaining the study and asking their consent to send identifying information to study staff at Stanford University who then contacted them to further describe the study and obtain written consent. The twin pairs described in this article were assessed starting in August 2005 and continuing until December 2009. Assessments were carried out by the Autism Genetic Resource Exchange in Southern California and Stanford University in Northern California. The statistical analysis of the data started in March 2010 and was completed in June 2010.
The full study protocol, including procedures for identifying and contacting families, consent and assent forms, blood draw and use, clinical assessments, and procedures for maintaining confidentiality of data, was reviewed and approved by the California Health and Welfare Agency Committee for the Protection of Human Subjects and by the institutional review board at Stanford University. Annual reviews were conducted, as were any revisions to protocol, and approval was secured from both institutional review boards. The protocol was reviewed as a research project with the expectation that results would be published in peer-reviewed medical journals. Data on self-identified race, Hispanic ethnicity, and other demographic variables were obtained from California birth certificates as provided by the mother.
ZYGOSITY DETERMINATION
Zygosity of sex-concordant twin pairs was determined in batches of 10 to 20 pairs concomitantly. Nine short tandem repeat loci and the X/Y amelogenin locus were amplified. Twin pairs discordant for at least 1 marker were considered dizygotic. Twin pairs concordant on all markers were considered monozygotic.
DIAGNOSTIC ASSESSMENTS
All twin participants were assessed using the Autism Diagnostic Interview–Revised (ADI-R),14 a structured parent interview, and the Autism Diagnostic Observation Schedule (ADOS),15 a standardized play and interview scale administered to children and adolescents. The ADI-R and ADOS are the most widely used and well-validated research measures for diagnosing autism and ASD. To assess verbal and nonverbal cognitive abilities, we used the Stanford-Binet Intelligence Scales, 5th edition. For children who did not achieve a basal score on the Stanford-Binet, the Mullen Scales of Early Learning Skills were used.
DIAGNOSTIC CLASSES
We defined 2 different diagnostic classifications based on the ADI-R and ADOS, one narrow (strict autism) and one broader (ASD). The diagnostic criteria for autism are divided into 4 areas: social impairments; communication impairments; repetitive, ritualistic, and stereotyped behaviors; and age at onset. Each of the 3 core areas (ie, social, communication, and repetitive behaviors) is defined by a minimum of 5 criteria. Individuals in the narrow category had to meet the following standard criteria for autism on both the ADI-R14 and the ADOS15: on the ADI-R, a child must score above a prespecified cut point on the 3 symptom areas and also have an age at onset prior to 3 years of age, and on the ADOS, the child must exceed specified thresholds on the domain scores for communication (including verbal and nonverbal behaviors) and reciprocal social interaction. The ASD category included individuals in the narrow category (strict autism), along with individuals who met a broader definition of ASD based on published criteria for combining information from the ADI-R and ADOS.16 According to these criteria, individuals were considered affected if they met criteria for ASD according to either the standard15 or the revised17 ADOS algorithm and met criteria for 1 of the following 3 categories: (1) met the ADI-R autism cutoffs on the domains of social reciprocity and either communication or restricted, repetitive behaviors or (2) were within 1 point of meeting the ADI-R autism cutoffs on the social and communication domains or met the autism cutoff on 1 of these domains and scored within 2 points on the other or (3) met full criteria on the ADI-R for autism on all of the domains.
STATISTICAL ANALYSIS
Definition of Probands
Twin pairs were ascertained through the DDS electronic files and client records, which constituted our sampling frame. In some cases, both twins in the pair were identified as having a high likelihood of autism in these records. In other cases, only 1 twin was so identified. On direct examination of the twins, we determined final diagnoses of strict autism or ASD as described earlier. For the purpose of our concordance and model-fitting analyses, we defined a proband as any individual who was identified in the DDS files or client records with qualifying DDS data and whose diagnosis of autism or ASD was subsequently confirmed by us through expert clinical review of abstracted data. Some individuals were diagnosed by us with autism or ASD but were not ascertained in the DDS file or client records; these individuals were not considered probands. Any individual who was determined by us to be unaffected was not counted as a proband. Also, in the analysis of strict autism, some individuals with a diagnosis of ASD but not strict autism were considered unaffected and not probands, even if they were identified through DDS records. Pairs where neither child met our criteria of strict autism were not included in the concordance analysis of strict autism.
Calculation of Concordance Rates
Pairwise concordance is calculated as the proportion of all ascertained pairs with 2 affected individuals and can be impacted by method of ascertainment of the twins.18,19 Proband-wise concordance is calculated as the proportion affected among the co-twins of probands. Probandwise concordance accounts for ascertainment and represents the recurrence risk to a co-twin of an affected individual.19 In numerical terms, let N = total number of twin pairs, R = number of pairs with 2 affected twins, t = number of pairs among the R that have 2 probands, S = number of pairs among the R that have 1 proband, and U = number of discordant pairs (with 1 proband). Then S + t = R and R + U = N. The pairwise concordance is simply given by R/N. The total number of co-twins of probands (which is the same as the number of probands because each proband has 1 co-twin) is 2 × t + S + U = N + t. The number of these co-twins who are also affected is 2 × t + S = R + t. In both of these cases, we doubled t because both affected twins in the pair are co-twins of probands. By contrast, we did not double S because only 1 affected twin in the pair is the co-twin of a proband. Hence, the probandwise concordance is given by (R + t)/(N + t).
For opposite-sex dizygotic twin pairs, the pairwise concordance is straightforward to calculate, as described earlier. For the probandwise concordance of opposite-sex twin pairs, we calculated the proportion of female co-twins of male probands who are affected and the proportion of male co-twins of female probands who are affected.
For pairwise concordance, 95% confidence intervals (CIs) were calculated exactly assuming a binomial distribution. For probandwise concordance, 95% CIs were calculated using the modified variance formula given by Davie.19
Examination of Potential Demographic, Ascertainment, or Clinical Biases
To test for differences (ie, potential biases) across several demographic factors between twin pairs or individuals who were included in the study vs those from the original DDS files who met entry criteria but were not included, a univariate analysis of variance was performed for continuous variables and χ2 tests were used for dichotomous traits. For both analyses, a P value of .05 for a 2-sided test was considered significant. To determine whether demographic factors influenced the concordance rates of monozygotic and dizygotic twin pairs, we used logistic regression separately for monozygotic and dizygotic pairs, with concordance vs discordance of the twin pair as the dichotomous outcome (dependent) variable. Independent variables included twin type (male or female for monozygotic twins; male, female, or male-female for dizygotic twins), age of mother and father, years of paternal and maternal education, ethnicity of mother and father, difference in birth weight between twin individuals, and gestational age (in completed weeks). For significance, a 2-sided P value of .05 adjusted for number of tests performed was considered. A critical assumption of the logistic regression is independence of the outcomes (in this case, concordance vs discordance among twin pairs). Because the twin pairs in our study were not related, independence is a reasonable assumption. All study-specific statistical analyses were performed using SPSS Statistics 19 (IBM SPSS, Chicago, Illinois).
To explore possible ascertainment bias, we compared the proportion of twin pairs concordant for being affected according to the qualifying DDS diagnostic categories among the twins included in the final study vs the twins from the original list who were ultimately not included. This analysis was stratified on type of twin pair (male-male, female-female, male-female). We could not stratify on zygosity because zygosity was unknown for pairs not included. For this analysis, we calculated odds ratios and χ2 tests for each of the 3 strata (defined by sex) and then a summary odds ratio and χ2 test across the 3 strata by Mantel-Haenszel analysis. The Mantel-Haenszel χ2 was assumed to have 1 df and was considered significant if greater than 3.84 (P < .05).
To determine whether diagnostic shifts may have occurred during the recruitment process, we compared the distribution of DDS diagnostic categories between affected twins included in the study vs those from the original list who were not included using a χ2 test, with df equal to the number of categories minus 1. Results were considered significant if the χ2 exceeded the threshold corresponding to a P value of .05.
Model Fitting
To determine the relative importance of genetic (heritable) vs shared co-twin environmental factors to the development of autism and ASD, we performed a classic ACE twin analysis and fit parametric models to the twin concordance data using maximum likelihood, incorporating information on proband status of the twins.20 The model assumes that numerous genetic and environmental factors of small effect accumulate into a continuous, normally distributed unobserved trait called liability. Liability values above a given threshold value lead to disease. Hence, this model is usually termed a multifactorial or polygenic threshold model.21 The critical parameters for this model are the components of variance of liability and the threshold value. The threshold value reflects disease prevalence. We assume different threshold values for males and females, which are denoted by Tm and Tf, respectively. For the modeling, we assumed a prevalence of 1% for males and 0.3% for females for ASD and 0.5% for males and 0.15% for females for strict autism. Because autism may be more prevalent in twins, we also considered the effect of assuming higher as well as lower prevalence estimates on our results.
Aside from the thresholds, 3 components of liability variance are included, a genetic component, an environmental component shared by twins, and a random environmental component. The parameter A represents the proportion of liability variance due to the genetic component; the parameter C represents the proportion of liability variance due to common twin environment; and the parameter E represents the proportion of liability variance due to random environment (ie, unshared by twins). When considering the sexes separately, we define 6 total parameters, Am and Af, Cm and Cf, and Em and Ef, corresponding to the values of A, C, and E for males and females, respectively. Note that Am + Cm + Em = 1 and Af + Cf + Ef = 1. Male values were estimated from the monozygotic and dizygotic concordance data in males, while female values were estimated from the monozygotic and dizygotic concordance data in females. Concordance in the opposite-sex dizygotic pairs allowed us to estimate the correlation between ½Am + Cm and ½Af + Cf, that is, the correlation between males and females for their combined genetic and shared environmental components. This correlation is denoted rmf. Because all monozygotic twin pairs are identical for sex, it is not possible to estimate separate correlations for Am and Af or Cm and Cf. For male monozygotic pairs, their correlation in liability is Am + Cm; for female monozygotic pairs, it is Af + Cf. For male dizygotic pairs, the correlation in liability is ½Am + Cm, and for dizygotic female pairs, the correlation in liability is ½Af + Cf. For opposite-sex dizygotic pairs, the correlation in liability is rmf [(½Am + Cm)(½Af + Cf)]½. Liability correlation (and hence variance components) was estimated from the twin concordance data. The model also includes the ascertainment probability π, which represents the probability that an affected individual is a proband (as described earlier). Models were fit to data using maximum likelihood with a grid search algorithm as implemented in Excel (Microsoft, Redmond, Washington), and hypotheses about model parameters were based on likelihood ratio tests. Best-fitting models were determined by the Akaike Information Criterion. Models with a low Akaike Information Criterion are the most parsimonious. Details are provided in the eAppendix (http://www.archgenpsychiatry.com).
RESULTS
A total of 1156 pairs fulfilled our initial eligibility criteria. We were unable to establish contact with 384 pairs, primarily because the families had lost contact with the regional center. Of the 772 remaining pairs, 330 declined participation and 10 pairs were ineligible. Of the 432 pairs whose contact information was sent to the Stanford team, 90 pairs declined participation in the study. Of the remaining 333 pairs, we completed the assessments for 202 pairs.
Among the 404 twin individuals, 242 met criteria for ASD (171 of whom also met criteria for strict autism), with the remaining 162 failing to meet criteria for an ASD. Comparing the study-qualifying diagnosis based on the DDS data with the diagnosis based on the ADI-R and ADOS, we calculated the sensitivity for ASD as 229/242 = 94.6% and the specificity as 137/162 = 84.6%. For strict autism, the sensitivity was 166/171 = 97.1% and the specificity, 158/233 = 67.8%. The DDS-based diagnosis was therefore a good predictor of the study ASD diagnosis based on the ADI-R and ADOS. For 10 twin pairs, neither of the twins met study criteria for the broader definition of ASD; these twin pairs were excluded from the genetic analyses.
To consider potential bias or differences among the 192 twin pairs who were included in the genetic analysis vs those who were not, we compared the 2 groups on a number of demographic criteria (eTable 1). The 2 groups were comparable for most variables, except that the age of the twins included in the genetic analyses was on average slightly older, birth weight of males was somewhat higher, and their mothers and fathers were, on average, slightly more educated and more likely to be white and less likely to be African American.
We also examined clinical differences by examining the distribution of DDS categories between the 2 groups (eTable 2). The proportional representation of the 6 diagnostic categories in the included vs nonincluded twins was comparable (; P = .12). We also looked to see if there was any correspondence between the 6 DDS diagnostic categories and our research diagnoses of strict autism vs ASD (eTable 2). Most DDS categories showed a similar distribution of strict autism vs ASD vs unaffected with the possible exception of “other,” for which a somewhat higher proportion (50%) was unaffected, although the number of individuals in this group was quite small (n = 6).
Of the 192 twin pairs included in the final analysis, 54 were monozygotic (45 male and 9 female) and 138 were dizygotic (45 male, 13 female, and 80 sex-discordant). The monozygotic twins were slightly older and had slightly shorter gestation periods (Table 1). In addition, the mothers of the dizygotic twins were older than the mothers of the monozygotic twins (Table 1), consistent with the known increase in dizygotic twinning with maternal age,22 and more likely to be white and non-Hispanic.
Table 1
Demographics of Twin Pairs Included in Genetic Analysis, California Autism Twins Study, 1987–2004 Birthsa
| Characteristic | Monozygotic | Dizygotic | P Value | ||
|---|---|---|---|---|---|
| Sample Size (No. Missing) | Mean (95% CI) | Sample Size (No. Missing) | Mean (95% CI) | ||
| Age, y | 54 (0) | 13.44 (12.45 to 14.43) | 138 (0) | 12.01 (11.38 to 12.65) | .02 |
| Mother’s age, y | 54 (0) | 30.98 (29.80 to 32.17) | 138 (0) | 33.42 (32.53 to 34.31) | .003 |
| Father’s age, y | 52 (2) | 34.46 (32.26 to 36.67) | 136 (2) | 36.01 (34.47 to 37.18) | .19 |
| Mother’s education, y | 53 (1) | 14.40 (13.49 to 15.30) | 135 (3) | 14.67 (14.25 to 15.10) | .53 |
| Father’s education, y | 50 (4) | 14.02 (13.05 to 14.99) | 133 (5) | 14.58 (14.16 to 15.00) | .22 |
| Gestational age, wk | 53 (1) | 35.38 (34.22 to 36.54) | 130 (8) | 36.52 (35.98 to 37.06) | .04 |
| Differential birth weight, g | 54 (0) | 294.44 (235.07 to 353.82) | 138 (0) | 307.25 (263.55 to 350.94) | .75 |
| Sample Size | No. (%) [95% CI] | Sample Size | No. (%) [95% CI] | P Value | |
| Ethnicity | |||||
| Mother | 54 | 138 | .01 | ||
| White | 25 (46.3) [33.00 to 59.60] | 85 (61.59) [53.47 to 69.71] | |||
| Hispanic | 23 (42.59) [29.40 to 55.78] | 27 (19.57) [12.95 to 26.19] | |||
| African American | 0 | 9 (6.52) [2.40 to 10.64] | |||
| Asian | 5 (9.26) [1.53 to 16.99] | 13 (9.42) [4.55 to 14.29] | |||
| Other | 1 (1.85) [−1.74 to 5.44] | 4 (2.80) [0.10 to 5.70] | |||
| Father | 54 | 138 | .44 | ||
| White | 30 (55.56) [42.31 to 66.81] | 89 (64.49) [56.51 to 72.47] | |||
| Hispanic | 15 (27.78) [15.83 to 39.73] | 27 (19.57) [12.95 to 26.19] | |||
| African American | 1 (1.85) [−1.74 to 5.44] | 8 (5.80) [1.90 to 9.70] | |||
| Asian | 6 (11.11) [2.73 to 19.48] | 11 (7.97) [3.45 to 12.49] | |||
| Other | 2 (3.70) [−1.33 to 8.73] | 3 (2.17) [−0.26 to 4.60] | |||
Abbreviation: CI confidence interval.
For strict autism (143 pairs), probandwise concordance rates (Table 2) for monozygotic twins were similar for 40 male pairs (58%; 95% CI, 42%–74%) and 7 female pairs (60%; 95% CI, 28%–90%), comparable with those reported in prior studies. However, probandwise concordance rates for dizygotic pairs were higher than previously reported, with rates of 21% for 31 male pairs (95% CI, 9%–43%) and 27% for 10 female pairs (95% CI, 9%–69%) (Table 2). The probandwise concordance rate for 54 female dizygotic co-twins of male probands was 3.7% (95% CI, 0.5%–13%), which was lower than for sexconcordant pairs, while the rate for male dizygotic co-twins of female probands was 50% (95% CI, 1%–99%). This latter finding was based on small numbers (n = 2).
Table 2
Twin Concordance Rates, California Autism Twins Study, 1987–2004 Births
| Diagnosis | Sex of Proband | Sex of Co-twin | Concordance Rate (Proportion) [95% CI] | |||
|---|---|---|---|---|---|---|
| Monozygotic | Dizygotic | |||||
| Pairwise | Probandwisea | Pairwise | Probandwisea | |||
| Strict autism (narrow) | M | M | 0.425 (17/40) [0.270–0.591] | 0.582 (32/55) [0.415–0.737] | 0.129 (4/31) [0.036–0.298] | 0.206 (7/34) [0.086–0.433] |
| F | F | 0.429 (3/7) [0.100–0.816] | 0.600 (6/10) [0.283–0.895] | 0.200 (2/10) [0.025–0.556] | 0.273 (3/11) [0.093–0.685] | |
| M | F | … | … | 0.036 (2/55) [0.004–0.125] | 0.037 (2/54) [0.005–0.128] | |
| F | M | … | … | … | 0.500 (1/2) [0.013–0.987] | |
| ASD (broad) | M | M | 0.644 (29/45) [0.488–0.781] | 0.771 (54/70) [0.650–0.864] | 0.200 (9/45) [0.096–0.346] | 0.308 (16/52) [0.159–0.463] |
| F | F | 0.333 (3/9) [0.075–0.701] | 0.500 (6/12) [0.157–0.843] | 0.308 (4/13) [0.091–0.614] | 0.357 (5/14) [0.108–0.604] | |
| M | F | … | … | 0.063 (5/80) [0.021–0.140] | 0.053 (4/76) [0.015–0.129] | |
| F | M | … | … | … | 0.500 (3/6) [0.118–0.882] | |
Abbreviations: ASD, autism spectrum disorder; CI, confidence interval.
For ASD (192 pairs), concordance estimates for both monozygotic and dizygotic twin pairs were generally higher. Probandwise concordance for monozygotic twins was 77% (95% CI, 65%–86%) for 45 male pairs and 50% (95% CI, 16%–84%) for 9 female pairs. Concordance rates for dizygotic twins were 31% (95% CI, 16%–46%) and 36% (95% CI, 11%–60%) for 45 male and 13 female pairs, respectively. Probandwise concordance for female dizygotic co-twins of 76 male probands was 5.3% (95% CI, 1.5%–13.0%), while the rate for male dizygotic co-twins of 6 female probands was 50.0% (95% CI, 11.8%–88.2%). Again, these dizygotic concordance rates are higher than previously reported and have a significant impact on the heritability analysis. From logistic regression analysis of concordance, the only factor impacting concordance rates was sex of co-twins in the dizygotic pairs (P value = .045), which was formally not significant if adjusted for multiple testing.
We also considered possible ascertainment bias in terms of concordance. Using the DDS categories to define affected, we examined the proportion of twin pairs concordant among the included pairs vs the nonincluded pairs, stratified by sex (eTable 1). Overall, there was no significant difference in proportion concordant (Mantel-Haenszel odds ratio = 1.17; 95% CI, 0.80–1.70) between the included and nonincluded pairs.
As described in the eAppendix, twin concordance rates were used to obtain parameter estimates for the genetic models; results are presented in Table 3 and likelihood ratio tests of different models, in Table 4. In all models, both genetic and shared environmental components were significant. The best-fitting models for both strict autism and ASD suggest that a large proportion of the variance in liability is due to shared environmental factors in addition to genetic heritability. For strict autism, we could conclude that heritabilities in males and females were equal and that the shared environmental components in males and females were also equal. The best-fitting model had a genetic heritability of 37% (95% CI, 8%–84%) and a shared environmental variance component of 55% (95% CI, 9%–81%). For the broader ASD phenotype, again we could conclude that heritabilities and shared twin environmental components were equal in males and females. Heritability was estimated to be 38% (95% CI, 14%–67%) and the shared environmental component, to be 58% (95% CI, 30%–80%). The shared environment component was estimated to be larger than the genetic heritability component. The poor fit of a pure heritability model can be attributed primarily to the high dizygotic twin concordance relative to the monozygotic twin concordance and population prevalence.
Table 3
Results of Model Fitting, California Autism Twins Study, 1987–2004 Birthsa
| Model | Parameter Constraints | π | Am | Cm | Af | Cf | rmf | −2lnL + C | AIC |
|---|---|---|---|---|---|---|---|---|---|
| Strict Autism (Narrow) | |||||||||
| Equal heritability for M and F; no shared twin environment | Am = Af; Cm = Cf = 0 | 0.896 | 0.931 | (0) | 0.931 | (0) | 1.00 | 7.83 | 11.83 |
| No shared twin environment | Cm = Cf = 0 | 0.895 | 0.926 | (0) | 0.951 | (0) | 1.00 | 7.64 | 13.64 |
| No heritability, equal twin environment for M and F, perfectly correlated | Am = Af = 0; Cm = Cf; rmf = 1 | 0.904 | (0) | 0.819 | (0) | 0.819 | (1.0) | 19.35 | 23.35 |
| No heritability, equal twin environment for M and F | Am = Af = 0; Cm = Cf | 0.902 | (0) | 0.862 | (0) | 0.862 | 0.578 | 9.01 | 13.01 |
| No heritability; twin environment perfectly correlated for M and F | Am = Af = 0; rmf = 1 | 0.904 | (0) | 0.837 | (0) | 0.757 | (1.0) | 18.86 | 22.86 |
| No heritability | Am = Af = 0 | 0.903 | (0) | 0.858 | (0) | 0.885 | 0.573 | 8.86 | 14.86 |
| Equal heritability for M and F; equal twin environment for M and F, perfectly correlated | Am = Af; Cm = Cf; rmf = 1 | 0.902 | 0.558 | 0.366 | 0.558 | 0.366 | (1.0) | 2.98 | 6.98 |
| Equal heritability for M and F; equal twin environment for M and F | Am = Af; Cm = Cf | 0.902 | 0.369 | 0.554 | 0.369 | 0.554 | 0.676 | 0.81 | 6.81 |
| Perfect correlation in heritability and twin environment between M and F | rmf = 1 | 0.902 | 0.638 | 0.282 | 0.444 | 0.492 | (1.0) | 2.68 | 10.68 |
| Full model | None | 0.902 | 0.479 | 0.442 | 0.235 | 0.703 | 0.668 | 0 | 10.00 |
| ASD (Broad) | |||||||||
| Equal heritability for M and F; no shared twin environment | Am = Af; Cm = Cf = 0 | 0.878 | 0.969 | (0) | 0.969 | (0) | 1.00 | 19.93 | 23.93 |
| No shared twin environment | Cm = Cf = 0 | 0.879 | 0.978 | (0) | 0.930 | (0) | 1.00 | 18.61 | 24.61 |
| No heritability, equal twin environment for M and F, perfectly correlated | Am = Af = 0; Cm = Cf; rmf = 1 | 0.894 | (0) | 0.864 | (0) | 0.864 | (1.0) | 35.00 | 37.00 |
| No heritability, equal twin environment for M and F | Am = Af = 0; Cm = Cf | 0.889 | (0) | 0.901 | (0) | 0.901 | 0.622 | 18.73 | 22.73 |
| No heritability; twin environment perfectly correlated for M and F | Am = Af = 0; rmf = 1 | 0.892 | (0) | 0.895 | 0) | 0.723 | (1.0) | 31.32 | 35.32 |
| No heritability | Am = Af = 0; | 0.889 | (0) | 0.905 | (0) | 0.882 | 0.627 | 18.55 | 24.55 |
| Equal heritability for M and F; equal twin environment for M and F; perfectly correlated | Am = Af; Cm = Cf; rmf = 1 | 0.889 | 0.513 | 0.450 | 0.513 | 0.450 | (1.0) | 7.28 | 11.28 |
| Equal heritability for M and F; equal twin environment for M and F | Am = Af; Cm = Cf | 0.888 | 0.383 | 0.581 | 0.383 | 0.581 | 0.724 | 3.28 | 9.28 |
| Perfect correlation in heritability and twin environment between M and F | rmf = 1 | 0.889 | 0.633 | 0.342 | 0.208 | 0.682 | (1.0) | 4.62 | 12.62 |
| Full model | None | 0.887 | 0.572 | 0.404 | 0.029 | 0.861 | 0.720 | 0.00 | 10.00 |
Abbreviations: Af, additive genetic variance for females; AIC, Akaike information criterion; Am, additive genetic variance for males; ASD, autism spectrum disorder; Cf, twin environmental variance for females; Cm, twin environmental variance for males; F, females; M, males; rmf, male-female correlation in variance components; π, ascertainment probability; −2lnL + C, minus twice the natural logarithm of the likelihood plus a constant.
Table 4
Tests of Hypotheses, California Autism Twins Study, 1987–2004 Births
| Hypothesis Tested | Parameter Constraint | df | Strict Autism (Narrow) | ASD (Broad) | ||
|---|---|---|---|---|---|---|
| χ2 | P Value | χ2 | P Value | |||
| Perfect correlation in liability between M and F | rmf = 1 | 1 | 2.68 | .10 | 4.62 | .03 |
| Equal heritability and equal twin environment for M and F | Am = Af; Cm = Cf | 2 | 0.81 | .67 | 3.28 | .19 |
| Equal heritability and equal twin environment for M and F, perfectly correlated | Am = Af; Cm = Cf; rmf = 1 | 3 | 2.98 | .40 | 7.28 | .06 |
| No heritability | Am = Af = 0 | 2 | 8.86 | .01 | 18.55 | <.001 |
| No twin environment | Cm = Cf = 0 | 2 | 7.64 | .02 | 18.61 | <.001 |
Abbreviations: Af, additive genetic variance for females; Am, additive genetic variance for males; ASD, autism spectrum disorder; Cf, shared twin environmental variance for females; Cm, twin environmental variance for males; F, females; M, males; rmf, male-female correlation in variance components.
COMMENT
To our knowledge, this study is the largest population-based twin study of autism that used contemporary standards for the diagnosis of autism. All twin individuals underwent a thorough diagnostic and cognitive examination that included a structured interview and observation that allowed differentiation between autism and other delays of development.13,14 The results suggest that environmental factors common to twins explain about 55% of the liability to autism. Although genetic factors also play an important role, they are of substantially lower magnitude than estimates from prior twin studies of autism. Nearly identical estimates emerged for ASD, suggesting that ASD presents the same liability spectrum as strict autism.
The California population, as represented in our twin sample, is highly diverse regarding ethnicity, socioeconomic status, and other demographic factors. Hence, our results should be readily generalizable, especially as compared with previous twin studies, which were based exclusively on individuals from Northern Europe. Even though this study was population based, the participation rate was only about 17% of eligible twin pairs. However, comparison of demographic and clinical characteristics of twins who did and did not participate showed only modest differences, primarily in parental education (eTable 1). Furthermore, concordance rates were not influenced by a range of potential confounding factors.
One possible concern is that a higher proportion of concordant dizygotic pairs participated than concordant monozygotic pairs. Examination of proband status among the monozygotic vs dizygotic pairs suggests this is not the case. Among 28 concordant monozygotic pairs, 24 had 2 probands, while 4 had 1 proband. By contrast, for the same-sex dizygotic pairs, among 13 concordant pairs, 8 had 2 probands and 5 had 1 proband. Thus, if anything, there might have been some underascertainment of concordant dizygotic pairs.
The ACE model we used has several inherent assumptions. First, it assumes no gene × environment interaction. If such a gene × environment effect does exist, it would be confounded with the A parameter in our analysis, implying that as an estimate of pure genetic effect, A may actually be an overestimate. Similarly, a critical assumption in the model is that the shared twin environmental effect is the same for monozygotic and dizygotic twins. If, in fact, monozygotic twins share the relevant environment to a greater degree than the dizygotic twins, some of the effect included in the parameter A would actually be environmental rather than genetic; again, A may actually overestimate the true genetic heritability.
Another potential limitation is the validity of the assumptions regarding prevalence rates of autism and ASD. We therefore compared heritability and shared twin environment estimates obtained by varying assumptions about prevalence. As shown in the eAppendix, we examined the impact of both doubling and halving our prevalence assumptions on derived parameter estimates. The range of heritability (eg, 41%–56% in males; 13%–16% in females) and shared twin environment (eg, 57%–41% in males; 78%–72% in females) estimates were similar. Hence, our conclusions regarding the relative importance of genetic and shared twin environment are quite robust to prevalence assumptions.
Our study provides evidence that the rate of concordance in dizygotic twins may have been seriously underestimated in previous studies and the influence of genetic factors on the susceptibility to develop autism, overestimated. Because of the reported high heritability of autism, a major focus of research in autism has been on finding the underlying genetic causes, with less emphasis on potential environmental triggers or causes. The finding of significant influence of the shared environment, experiences that are common to both twin individuals, may be important for future research paradigms.23 Increasingly, evidence is accumulating that overt symptoms of autism emerge around the end of the first year of life. Because the prenatal environment and early postnatal environment are shared between twin individuals, we hypothesize that at least some of the environmental factors impacting susceptibility to autism exert their effect during this critical period of life. Nongenetic risk factors that may index environmental influences include parental age,24 low birth weight,25 multiple births,26 and maternal infections during pregnancy.27 Future studies that seek to elucidate such factors and their role in enhancing or suppressing genetic susceptibility are likely to enhance our understanding of autism.23
Acknowledgments
Funding/Support: The study was supported by grant R01MH067005 from the National Institute of Mental Health and a grant from Autism Speaks.
Role of the Sponsors: None of the foundations had a role in the design and conduct of the study; the collection, management, analysis, and interpretation of the data; or the preparation, review, or approval of the manuscript.
Footnotes
Author Contributions: Dr Hallmayer had full access to all the data in the study and takes responsibility for the integrity of the data and the accuracy of the data analysis.
Financial Disclosure: None reported.
Online-Only Material: The eAppendix and eTables are available at http://www.archgenpsychiatry.com.