- Open Access
Assessing the impact of a combined analysis of four common low-risk genetic variants on autism risk
Molecular Autism volume 1, Article number: 4 (2010)
Autism is a complex disorder characterized by deficits involving communication, social interaction, and repetitive and restrictive patterns of behavior. Twin studies have shown that autism is strongly heritable, suggesting a strong genetic component. In other disease states with a complex etiology, such as type 2 diabetes, cancer and cardiovascular disease, combined analysis of multiple genetic variants in a genetic score has helped to identify individuals at high risk of disease. Genetic scores are designed to test for association of genetic markers with disease.
The accumulation of multiple risk alleles markedly increases the risk of being affected, and compared with studying polymorphisms individually, it improves the identification of subgroups of individuals at greater risk. In the present study, we show that this approach can be applied to autism by specifically looking at a high-risk population of children who have siblings with autism. A two-sample study design and the generation of a genetic score using multiple independent genes were used to assess the risk of autism in a high-risk population.
In both samples, odds ratios (ORs) increased significantly as a function of the number of risk alleles, with a genetic score of 8 being associated with an OR of 5.54 (95% confidence interval [CI] 2.45 to 12.49). The sensitivities and specificities for each genetic score were similar in both analyses, and the resultant area under the receiver operating characteristic curves were identical (0.59).
These results suggest that the accumulation of multiple risk alleles in a genetic score is a useful strategy for assessing the risk of autism in siblings of affected individuals, and may be better than studying single polymorphisms for identifying subgroups of individuals with significantly greater risk.
Autism is a heterogeneous disorder characterized by impairments in social interaction, deficits in verbal and nonverbal communication, restricted interests and repetitive behaviors . Autism comprises the severe end of a group of autism spectrum disorders (ASD) . The prevalence of autism is estimated at 0.2%, with males being more likely than females to have a diagnosis of autism (ratio of approximately 4:1) .
There is compelling evidence from twin and family studies indicating a strong genetic component in autism. The average risk of recurrence of autism in siblings is approximately 10%  in families with one affected sibling, which is much higher than the prevalence in the general population, but much lower than would be expected for a single-gene disease . Indeed, cases of autism have only rarely been associated with mutations in single genes or with chromosomal duplications or other abnormalities [5–9]. Furthermore, the inheritance pattern in most cases of autism is complex and not compatible with simple Mendelian inheritance [10, 11]. A series of common autism susceptibility genes has been identified in the literature, but taken individually, the risk associated with these genes remains modest. Although spontaneous mutations may help explain isolated cases of autism, the inheritance pattern of common variants may be at the root of autism in a large proportion of multiplex families.
Studies using simulated data have demonstrated that the joint analysis of such common low-risk variants has the potential to identify individuals at higher risk for disease, and could be useful for complex disease prediction [12–14]. In complex disease states such as type 2 diabetes [15–19], cancer [20, 21] and cardiovascular disease [22–25], the accumulation of multiple risk alleles markedly increases the risk of being affected, and compared with studying single nucleotide polymorphisms (SNPs) independently, it is better at identifying subgroups of individuals with significantly greater risk. Previous studies have shown associations between several genetic variants and autism, but many of these studies could not be replicated in different populations.
For our study, we selected single nucleotide polymorphisms (SNPs) within four genes that have been previously shown to be associated with autism. Original publications and replication studies have shown that SLC25A12 [26–28] and EN2 [29, 30] are associated with autism. PITX1 has been identified through a linkage study  using the physical identity by descent method . The same method was used to identify ATP2B2 on chromosome 3  in a region in which different autism linkage peaks and association with microsatellite markers have been reported [34–36]. We sought to demonstrate that joint analysis of common low-risk variants would have the potential to identify a subgroup of individuals at increased risk for developing autism. To assess the effect of combined analysis of multiple genetic variants in autism, we conducted a primary analysis on the first family sample using the four genes mentioned above. The aims of this analysis were to assess the strength of the association of the combination of these four genes with autism in a population of affected children, and to define a genetic scoring metric. We followed up with a second analysis to confirm our findings in an independent set of families, including a sibling diagnosed with autism.
The study design involved two independent family samples. The first sample came from the Autism Genetic Resource Exchange (AGRE; http://www.agre.org) repository (AGRE sample), and consisted of 222 families with at least two affected offspring per family. In total, 527 affected siblings met the diagnostic criteria for 'narrow autism' according to the Autism Diagnostic Interview Revisited (ADI-R). The male to female ratio was 3.5:1, which reflected the gender distribution of autism in the general population. In this sample, all parents and affected children were genotyped at the four SNPs, whereas unaffected siblings were not genotyped.
The second (Seattle) sample was part of the University of Washington Family Study of Autism collection funded by the National Institute of Child Health and Human Development (NICHD), and consisted of 241 families with at least two offspring, at least one of whom was affected by autism. Additional criteria were used to categorize the 461 affected siblings with 'narrow autism,' including the ADI-R , the Autism Diagnostic Observation Schedule (ADOS-G)  scores, and a clinical diagnosis by an experienced clinician . The male to female ratio was 3.7:1 and included 127 unaffected siblings as controls. The unaffected siblings were classified based on parental report in the Family History Interview  and the Broader Phenotype of Autism Symptom Scale .
In the Seattle sample, the IQ score as measured by the Wechsler scale in affected siblings was (mean ± SD) 76.3 ± 26.1, and age at IQ measurement was 9.4 ± 4.1. In unaffected siblings, the IQ score was 111 ± 16.3 with age at IQ measurement being 11.6 ± 6.4, which was comparable with the mean IQ score in parents (112.0 ± 14.0). All parents, affected siblings and unaffected siblings were genotyped at the four SNPs.
In the AGRE sample, single SNP analyses were conducted primarily with the goal of verifying that the preselected risk alleles corresponded with previous reports in the literature. The primary goal with the Seattle sample was to assess if the risk scores first generated in the AGRE sample could be replicated in the Seattle sample using 'true' sibling controls.
For both samples, exclusion criteria for children included a diagnosis of Rett syndrome and childhood disintegrative disorder, as defined by the Diagnostic and Statistical Manual of Mental Disorders, fourth edition (DSM-IV) criteria for other pervasive developmental disorders, presence of a known genetic condition, history of serious head injury or neurologic disease, or significant sensory or motor impairment . Ethnicity was self-reported by parents as Caucasian, Asian, Hispanic or Latino, Black or African American, Native Hawaiian or other Pacific Islander, or of mixed ethnicities. Caucasians represented the major ethnicity, with more than two-thirds of families in each sample.
The study was approved by the respective institutional review boards (IRBs), and written informed consent was obtained from all authorized representatives. The AGRE sample had been approved by an IRB to ensure the protection of research participants. For the Seattle sample, all work with autism subjects was approved by the University of Washington Human Subjects Division Review Committees.
SNP selection and genotyping
Four SNPs located in four separate genes from a pool of genetic polymorphisms that have been associated with autism were included in this multigene study: rs2292813 in SLC25A12 [26, 27], rs1861972 in EN2 [29, 30], rs35678 in ATP2B2  and rs6872664 in PITX1 . Genotyping was performed using TaqMan allele discrimination assays (Applied Biosystems, Foster City, CA, USA). Genotyping was performed in 384-well plates with 5 ng genomic DNA, 0.075 μL of 20 × SNP TaqMan Assay mix, 1.5 μL of TaqMan Universal PCR Master Mix and 1.425 μL of dH2O in each well. PCR was performed at 95°C for 10 minutes, followed by 50 cycles at 92°C for 15 seconds and 60°C for 90 seconds (9700 Gene Amp PCR System; Applied Biosystems). Plates were then subjected to end-point reading (7900 Real-Time PCR System; Applied Biosystems). Allele calls were assigned automatically, and a visual inspection of genotype clusters was performed. Genotyping quality was assessed by signal intensity plots and missing genotype frequencies; any samples with poor clustering and missing fractions ≥ 5% per SNP were re-typed. The genotyping success rate was found to be 97.4%. Parents' genotypes were used to investigate Hardy-Weinberg equilibrium (HWE) and to check for Mendelian inconsistencies. Families with remaining inconsistencies were excluded.
Allelic frequencies were estimated for each variant in parents in both population samples. All analyses excluded the index case, defined as the oldest affected child in the family. Odds ratios (ORs), confidence intervals (CIs) and P-values were first calculated for all four variants individually. The sibling recurrence risk ratio attributable to the four markers was estimated as described previously . Genetic scores for combinations of SLC25A12, ATP2B2, EN2 and PITX1 SNPs were then obtained using the allele count model [15, 16, 20, 24], which takes into account all possible combinations of risk alleles by adding the number of alleles carried by each individual under an additive assumption expected for ATP2B2 (see Additional file 1). Individuals with an ATP2B2 risk allele had two points added to their genetic score if they were a homozygous carrier of the risk allele and zero points in all other cases, because the association of ATP2B2 with autism follows a recessive transmission model . A subgroup analysis was conducted for the Caucasian families only because of the small sample sizes for all other ethnicities. Another subgroup analysis was performed separately for males and females using the risk score as an independent variable. For both subgroup analyses, the 'narrow' phenotype definition was used, and the analyses were conducted using the logistic regression iterative estimating equation (IEE).
Logistic regression was used to examine the association between the genetic score and the risk of autism, calculating the OR and CI associated with each additional risk allele in the genetic score. Following the definition of the reference group for single SNP OR estimates, the reference group for the risk score was chosen based on previous data [18, 19, 43] as the lower risk score. The discriminative ability of the model was evaluated, estimating sensitivity and specificity. The sensitivity describes the proportion of affected individuals having a genetic score above a given γ threshold. The specificity describes the proportion of unaffected individuals not having a genetic score above a given γ threshold. The receiver operating characteristic (ROC) curve was created by plotting sensitivity against 1 minus specificity, and the area under the curve (AUC) was calculated using the C-statistic.
The AGRE sample did not include genotyped unaffected siblings; specificity was not directly derived from data from unaffected individuals. We combined the case-pseudocontrol approach proposed by Cordell and Clayton [44, 45] with a formula adapted from DeLong et al. , giving the specificity (spec) for a specific genetic score threshold γ as a function of the OR (ORγ) and the sensitivity (sensγ) as given in the equation below:
The OR associated with the γ threshold was calculated by conditional logistic regression using cases and their matched pseudocontrols. For each child with autism, a pseudocontrol was constructed from parental untransmitted alleles to the child with autism. Case-pseudocontrol statistical tests (developed by D. Clayton, Cambridge Institute for Medical Research; available at: http://www-gene.cimr.cam.ac.uk/clayton/software) were implemented within R (statistical program language). A version that was robust to non-independent siblings was implemented using the DGCgenetics R library, allowing for multiple affected offspring from the same nuclear family to be used in this study. In the Seattle sample, the analyses were performed as a case-control study in unaffected and affected siblings of the proband, using a generalized estimating equation (GEE) with a 'logit' link and the independence working correlation matrix to adjust for familial correlation . The case-pseudocontrol design provides estimates of genetic relative risk [44, 45] rather than OR, as in the case-control design. However, under the rare disease assumption, both estimates are equivalent.
A supplementary analysis using a GEE to correct for sibling relationship was conducted in the Seattle sample to check for a possible association between the genetic score and the Wechsler IQ scoring. Analyses were conducted first on the whole case-control Seattle sample, adjusting for gender and autism status, and in affected and unaffected siblings adjusting for gender.
None of the SNPs exhibited a departure from HWE, and allele frequencies were similar between samples (Table 1). All risk alleles were very common, with frequencies ranging from 0.4 (ATP2B2) to 0.9 (PITX), similar to those estimated in the literature [26, 27, 29–31, 33]. The distribution of case subjects in both samples was similar for each genetic score category (range 2 to 8), as was the distribution of pseudocontrol and control subjects. Control subjects were more often in the lower genetic score categories (≤ 4), and the case subjects were more often in the higher genetic score categories (> 6) (Figure 1).
A primary analysis using the AGRE sample was performed to assess the association with autism of each individual SNP used in these analyses. The ORs for development of autism in siblings of affected children were calculated for each SNP using data from this sample. Nominally significant ORs were observed for the risk alleles for EN2 (P = 0.02), PITX (P = 0.01) and ATP2B2 (P = 0.01) (Table 2). Despite its reported association with autism [26–28], the OR for SLC25A12 did not reach significance (P = 0.38). An analysis was performed in which a genetic score was calculated (see Additional file 1), using the number of risk alleles carried by a subject, and testing for value in predicting the risk of autism in siblings of affected children (Table 3). A one-point increase in genetic score significantly increased the OR for autism in siblings of affected children (OR = 1.33; P < 0.0001). Individuals with a genetic score of 8 had an increased risk of being affected by autism (OR = 5.54; 95% CI 2.45 to 12.49) compared with individuals with two risk alleles (reference group). The same analysis was performed in the independent Seattle sample to confirm these findings. A similar distribution of genetic scores and ORs was observed in this analysis, with a one-point genetic score increase giving OR = 1.32 (95% CI 1.07 to 1.63) (Table 3). More than 50% of cases in both analyses had a genetic score of five or six risk alleles. In the AGRE sample, a genetic score of 5 or 6 gave ORs of 2.35 (95% CI 1.57 to 3.53) and 3.13 (95% CI 1.82 to 5.38), respectively, and in the Seattle sample, gave ORs of 2.30 (95% CI 1.23 to 4.31) and 3.04 (95% CI 1.32 to 7.01), respectively. Odds ratios for a one-point increase in the genetic scores in the Caucasian subgroup from the AGRE sample (OR = 1.39; 95% CI 1.17 to 1.66; see Additional file 2) and the Seattle sample (OR = 1.23; 95% CI 1.02 to 1.50; see Additional file 2) were similar to the ORs (1.33 and 1.32, respectively) estimated for the complete samples (Table 3). The ORs for gender subgroups were very similar in both samples (see Additional file 2).
The sensitivities and specificities for each genetic score were similar in the two samples, with overlapping CIs (Table 4). The resultant ROC curves for the genetic scores estimated in the two analyses (Figure 2) gave an AUC of 0.59 for both the AGRE (95% CI 0.55 to 0.64) and Seattle (95% CI 0.53 to 0.65) samples.
In the Seattle sample, for which true controls were available, 17% of individuals with autism carried seven risk alleles (genetic score ≥ 7) compared with only 8% of children without autism (Figure 1). Low sensitivity of 0.17 (95% CI 0.12 to 0.22) was associated with a high specificity of 0.93 (95% CI 0.87 to 0.98), which translated into a 2.33-fold increased risk for autism (Table 4). In the Seattle sample, 8% of cases had a maximum genetic score of 8, compared with only 2% of unaffected siblings of individuals with autism. The specificities for an extreme genetic score of 8 reached 0.98 (95% CI 0.96 to 1.00), but this was due to reduction in sensitivity (0.08; 95% CI 0.04 to 0.12), leading to a 4.44 increased risk of autism.
Based on these analyses and the distribution of the frequency of genetic score, we defined a low-score group of ≤ four risk alleles, which included approximately 23% of siblings, regardless of phenotype, in both samples. An intermediate-score group was defined as five or six risk alleles covering 62% and 65% of siblings and a high-score group of ≥ seven risk alleles covered 15% and 12% of siblings in the AGRE and Seattle samples, respectively. Using the low-score group as the reference group, the ORs for autism in the intermediate-risk group were 1.61 (95% CI 1.23 to 2.12) and 1.78 (95% CI 1.12 to 2.83) in the AGRE and the Seattle samples, respectively. Similarly, ORs for the high-risk group were 2.59 (95% CI 1.50 to 4.47) in the AGRE sample and 3.15 (95% CI 1.25 to 7.99) in the Seattle sample (Table 5).
We did not observe any association between the genetic score and IQ score in the whole case-control Seattle sample (P = 0.44), when adjusted for autism status and gender. This was also the case in affected (P = 0.63) and unaffected (P = 0.68) siblings adjusted for gender.
Complex pathologies such as autism may involve multiple common low-risk variants that confer relatively small predisposing effects. In the present study, we hypothesized that the combined analysis of such low-risk variants would more accurately identify a group of individuals at higher risk of developing the disease than an analysis of independent variants. To our knowledge, this is the first study report of its kind in autism.
Although many variants have been described in the literature as being associated with autism, few associations have been confirmed, possibly due to differences in methods and sample sets used, and to the weak effect of the individual variants. Previous studies have shown that variants in SLC25A12 [26–28], EN2 [29, 30], ATP2B2  and PITX1 , have an association with autism. A few studies have physiologically linked these genes to autism as well.
Expression of SLC25A12 may be involved in the pathophysiology of autism  because postmortem samples of brain tissue from patients with autism have been found to have stronger expression of SLC25A12 than normal brain samples. The EN2 gene is specifically involved in patterning the region that gives rise to the cerebellum . EN2 knockout mice exhibit neuroanatomical and behavioral abnormalities that partly resemble those in patients with autism . PITX1 is a key regulator of hormonal genes in the pituitary-hypothalamic axis. Its putative involvement in autism is supported by evidence documenting abnormal levels of downstream hormones, such as adrenocorticotropic hormone (ACTH), beta-endorphin and cortisol in individuals with autism [51–53]. Deregulation of pro-opiomelanocortin and high levels of beta-endorphin in the morning, for example, have been shown to be involved in certain maladaptive behaviors, such as self-injurious behaviors, which are often seen in individuals with autism . The ACTH-cortisol system, which also plays an important role in stress-related responses, is impaired in individuals with autism in whom lower cortisol levels and higher ACTH levels have been reported [44, 52]. Finally, the ATP2B2 gene codes for a major calcium pump expressed at particularly high levels in Purkinje neurons. Purkinje cells are involved in motor coordination, working memory and learning. Even though the exact correlation between Purkinje neurons and autism remains unclear, selective loss of Purkinje cells and cerebellar atrophies are the most consistently found neurological abnormalities in people diagnosed with autism .
Odds ratios ranging from 1.16 to 1.74 were observed in this study, and P values supported a nominal association for individual variants of EN2, ATP2B2 and PITX1, but not SLC25A12. However, because of its reported genetic association with autism [26, 27, 29] and the possible relationship with autism of its overexpression , it was included in the combined analysis. In the first stage of the combined assessment, a primary analysis was performed to evaluate the effect of the combination of variants as determined by the genetic score. In the second stage, a case-control analysis was performed to confirm the findings of the AGRE sample analysis. The distribution of genetic scores and ORs in this analysis was similar to that observed in the AGRE sample analysis, with an OR of 1.32. The ORs for a one-point increase in genetic score for the Caucasian subgroup in both the AGRE and Seattle sample analyses were also similar to those of the overall respective sample sets. We estimated that the OR for siblings who carry > seven risk alleles (approximately 12% of siblings), regardless of phenotype, was 3.15, compared with siblings who carry < four risk alleles (approximately 23% of siblings) in the Seattle sample, which used true sibling controls. A subgroup analysis did not yield any significant differences based on gender. Because of sample size, only the Caucasian subgroup could be analyzed. The ORs for a one-point increase in genetic score for the Caucasian subgroup in both the AGRE and Seattle sample analyses were also similar to those of the overall respective sample sets. We did not observe any association between IQ and genetic score, thus, we believe that the observed association is not an effect of the different IQ distributions between cases and controls.
The sensitivity and specificity for each genetic score was similar in both analyses. The resultant ROC curves for the genetic score estimated in these analyses produced an AUC of 0.59, significantly different (P < 0.01) from a random prediction of the score (AUC = 0.50). These results may be considered a first step toward the goal of using common genetic variants at multiple loci to develop a multigene approach to predicting the risk of autism. Indeed, these data are in line with other multigene risk assessment models in other complex diseases, such as type 2 diabetes, with reported AUCs of approximately 0.6 [15, 16, 19, 56] using 3 to 18 common genetic variants. Results from these studies provide evidence that accumulation of multiple risk alleles markedly increases the risk of diabetes. However, compared with known nongenetic risk factors, such as body mass index, age and family history of the disease, combination of risk alleles provides only a relatively small increase in risk prediction [19, 56, 57]. In the case of autism, no early predictors exist to compare with genetic score. Nevertheless, there is significant interest in the early identification of infants at higher risk for autism because studies have shown that early intervention leads to significantly improved long-term outcome for the whole family [58, 59].
The AUC is not a clinically relevant way of summarizing predictive performances. Furthermore, the AUC is a poor metric for evaluating markers for disease diagnoses, screening or prognosis . Using a threshold of six risk alleles allows identification of 18% of cases (19% and 17% in the AGRE and Seattle samples, respectively) while reducing the likelihood of false-positive results (around 8%). Importantly, these sensitivity and specificity values translate into a > 2-fold increased risk for autism. In such a disorder, defined entirely on the basis of child behavior, this can be considered an important first step toward improved risk prediction.
Several studies support a strong genetic contribution to the etiology of autism, including twin studies that estimated a heritability of > 90% [61, 62] and a family study that indicated a sibling recurrence risk ratio of 22 , which indicate that there is much potential for genetic prediction in autism. In the current study, the four SNPs used in the score are associated with a sibling recurrence risk ratio of 1.8. These results demonstrate that up to 10% of the variation in risk could be explained by these four genes.
Several strategies exist to identify relevant variants to be combined in a genetic score. These range from conservative approaches using only well-established factors [16, 56, 64, 65] to more liberal approaches, which include polymorphisms from candidate genes or from regions that have already been associated with the disease or that display functional relevance [20, 22, 66]. We selected a subset of four variants from four candidate genes already associated with autism. An additional strategy would have been to include genes identified in a genome-wide association study . The approach described here could be applied to other established variants and be part of future studies. Only one recent paper reported genome-wide scan results, which identified a common genetic variant reaching genome-wide significance located on 5p14.1 in the intergenic region between the CDH10 and CDH9 genes . The reported ORs for the associated variant were similar to those of the variants in this study. However, it would be interesting to add this variant into the risk score. In our study, we chose the risk allele count as the most commonly used method of constructing genetic scores from low-effect genetic variants [15, 16, 20, 24]. We used a recessive model to calculate the effect of ATP2B2, thereby strengthening our genetic score. Weighted scores have been proposed to increase discriminative power, but models that introduce information from variants with similar, small or moderate effect size only slightly increased the AUC [19, 55].
One limitation of our study was that we assumed an individual effect for each gene. Gene by gene and/or gene by environment interaction should be tested and could be introduced into the genetic score for a more accurate assessment of risk associated with a given genetic score. Another limitation was that the polymorphisms used are probably not causal variants but rather markers in linkage disequilibrium. Causal variants may have different allele distributions in siblings with and without autism, and could have a stronger effect on disease risk. This is not yet clearly established and needs to be clarified .
Surprisingly, for three of the four loci investigated in this study, the major allele was overtransmitted in cases versus controls. Although it is unclear whether this same effect may hold for other autism-related risk loci, the observation is interesting and is one that merits further consideration in the field.
In conclusion, the results of this study suggest that the joint analysis of common-risk variants will be more effective in identifying subgroups of individuals at risk than using single polymorphisms for risk assessment. The developed score could be improved by adding additional genes and environmental factors in the future.
Johnson CP, Myers SM: Identification and evaluation of children with autism spectrum disorders. Pediatrics. 2007, 120: 1183-1215. 10.1542/peds.2007-2361.
Muhle R, Trentacoste SV, Rapin I: The genetics of autism. Pediatrics. 2004, 113: e472-e486. 10.1542/peds.113.5.e472.
Fombonne E: Epidemiology of pervasive developmental disorders. Pediatr Res. 2009, 65: 591-598. 10.1203/PDR.0b013e31819e7203.
Constantino JN, Lajonchere C, Lutz M, Gray T, Abbacchi A, McKenna K, Singh D, Todd RD: Autistic social impairment in the siblings of children with pervasive developmental disorders. Am J Psychiatry. 2006, 163: 294-296. 10.1176/appi.ajp.163.2.294.
Durand CM, Betancur C, Boeckers TM, Bockmann J, Chaste P, Fauchereau F, Nygren G, Rastam M, Gillberg IC, Anckarsater H, Sponheim E, Goubran-Botros H, Delorme R, Chabane N, Mouren-Simeoni MC, de Mas P, Bieth E, Rogé B, Héron D, Burglen L, Gillberg C, Leboyer M, Bourgeron T: Mutations in the gene encoding the synaptic scaffolding protein SHANK3 are associated with autism spectrum disorders. Nat Genet. 2007, 39: 25-27. 10.1038/ng1933.
Jamain S, Quach H, Betancur C, Rastam M, Colineaux C, Gillberg IC, Soderstrom H, Giros B, Leboyer M, Gillberg C, Bourgeron T, Paris Autism Research International Sibpair Study: Mutations of the X-linked genes encoding neuroligins NLGN3 and NLGN4 are associated with autism. Nat Genet. 2003, 34: 27-29. 10.1038/ng1136.
Veenstra-VanderWeele J, Cook EH: Molecular genetics of autism spectrum disorder. Mol Psychiatry. 2004, 9: 819-832. 10.1038/sj.mp.4001505.
Sebat J, Lakshmi B, Malhotra D, Troge J, Lese-Martin C, Walsh T, Yamrom B, Yoon S, Krasnitz A, Kendall J, Leotta A, Pai D, Zhang R, Lee YH, Hicks J, Spence SJ, Lee AT, Puura K, Lehtimäki T, Ledbetter D, Gregersen PK, Bregman J, Sutcliffe JS, Jobanputra V, Chung W, Warburton D, King MC, Skuse D, Geschwind DH, Gilliam TC, Ye K, Wigler M: Strong association of de novo copy number mutations with autism. Science. 2007, 316: 445-449. 10.1126/science.1138659.
Weiss LA, Shen Y, Korn JM, Arking DE, Miller DT, Fossdal R, Saemundsen E, Stefansson H, Ferreira MA, Green T, Platt OS, Ruderfer DM, Walsh CA, Altshuler D, Chakravarti A, Tanzi RE, Stefansson K, Santangelo SL, Gusella JF, Sklar P, Wu BL, Daly MJ, Autism Consortium: Association between microdeletion and microduplication at 16p11.2 and autism. N Engl J Med. 2008, 358: 667-675. 10.1056/NEJMoa075974.
Szatmari P, Paterson AD, Zwaigenbaum L, Roberts W, Brian J, Liu XQ, Vincent JB, Skaug JL, Thompson AP, Senman L: Mapping autism risk loci using genetic linkage and chromosomal rearrangements. Nat Genet. 2007, 39: 319-328. 10.1038/ng1985.
Jorde LB, Hasstedt SJ, Ritvo ER, Mason-Brothers A, Freeman BJ, Pingree C, McMahon WM, Petersen B, Jenson WR, Mo A: Complex segregation analysis of autism. Am J Hum Genet. 1991, 49: 932-938.
Pharoah PD, Antoniou A, Bobrow M, Zimmern RL, Easton DF, Ponder BA: Polygenic susceptibility to breast cancer and implications for prevention. Nat Genet. 2002, 31: 33-36. 10.1038/ng853.
Yang Q, Khoury MJ, Botto L, Friedman JM, Flanders WD: Improving the prediction of complex diseases by testing for multiple disease-susceptibility genes. Am J Hum Genet. 2003, 72: 636-649. 10.1086/367923.
Janssens AC, Aulchenko YS, Elefante S, Borsboom GJ, Steyerberg EW, van Duijn CM: Predictive testing for complex diseases using multiple genes: fact or fiction?. Genet Med. 2006, 8: 395-400. 10.1097/01.gim.0000229689.18263.f4.
Weedon MN, McCarthy MI, Hitman G, Walker M, Groves CJ, Zeggini E, Rayner NW, Shields B, Owen KR, Hattersley AT, Frayling TM: Combining information from common type 2 diabetes risk polymorphisms improves disease prediction. PLoS Med. 2006, 3: e374-10.1371/journal.pmed.0030374.
Lango H, Palmer CN, Morris AD, Zeggini E, Hattersley AT, McCarthy MI, Frayling TM, Weedon MN: Assessing the combined impact of 18 common genetic variants of modest effect sizes on type 2 diabetes risk. Diabetes. 2008, 57: 3129-3135. 10.2337/db08-0504.
Lyssenko V, Almgren P, Anevski D, Orho-Melander M, Sjogren M, Saloranta C, Tuomi T, Groop L: Genetic prediction of future type 2 diabetes. PLoS Med. 2005, 2: e345-10.1371/journal.pmed.0020345.
Lu Q, Elston RC: Using the optimal receiver operating characteristic curve to design a predictive genetic test, exemplified with type 2 diabetes. Am J Hum Genet. 2008, 82: 641-651. 10.1016/j.ajhg.2007.12.025.
Lin X, Song K, Lim N, Yuan X, Johnson T, Abderrahmani A, Vollenweider P, Stirnadel H, Sundseth SS, Lai E, Burns DK, Middleton LT, Roses AD, Matthews PM, Waeber G, Cardon L, Waterworth DM, Mooser V: Risk prediction of prevalent diabetes in a Swiss population using a weighted genetic score--the CoLaus Study. Diabetologia. 2009, 52: 600-608. 10.1007/s00125-008-1254-y.
Zheng SL, Sun J, Wiklund F, Smith S, Stattin P, Li G, Adami HO, Hsu FC, Zhu Y, Balter K, Kader AK, Turner AR, Liu W, Bleecker ER, Meyers DA, Duggan D, Carpten JD, Chang BL, Isaacs WB, Xu J, Grönberg H: Cumulative association of five genetic variants with prostate cancer. N Engl J Med. 2008, 358: 910-919. 10.1056/NEJMoa075819.
Gail MH: Discriminatory accuracy from single-nucleotide polymorphisms in models to predict breast cancer risk. J Natl Cancer Inst. 2008, 100: 1037-1041. 10.1093/jnci/djn180.
Kathiresan S, Melander O, Anevski D, Guiducci C, Burtt NP, Roos C, Hirschhorn JN, Berglund G, Hedblad B, Groop L, Altshuler DM, Newton-Cheh C, Orho-Melander M: Polymorphisms associated with cholesterol and risk of cardiovascular events. N Engl J Med. 2008, 358: 1240-1249. 10.1056/NEJMoa0706728.
Martinelli N, Trabetti E, Pinotti M, Olivieri O, Sandri M, Friso S, Pizzolo F, Bozzini C, Caruso PP, Cavallari U, Cheng S, Pignatti PF, Bernardi F, Corrocher R, Girelli D: Combined effect of hemostatic gene polymorphisms and the risk of myocardial infarction in patients with advanced coronary atherosclerosis. PLoS ONE. 2008, 3: e1523-10.1371/journal.pone.0001523.
Morrison AC, Bare LA, Chambless LE, Ellis SG, Malloy M, Kane JP, Pankow JS, Devlin JJ, Willerson JT, Boerwinkle E: Prediction of coronary heart disease risk using a genetic risk score: the Atherosclerosis Risk in Communities Study. Am J Epidemiol. 2007, 166: 28-35. 10.1093/aje/kwm060.
Humphries SE, Ridker PM, Talmud PJ: Genetic testing for cardiovascular disease susceptibility: a useful clinical management tool or possible misinformation?. Arterioscler Thromb Vasc Biol. 2004, 24: 628-636. 10.1161/01.ATV.0000116216.56511.39.
Segurado R, Conroy J, Meally E, Fitzgerald M, Gill M, Gallagher L: Confirmation of association between autism and the mitochondrial aspartate/glutamate carrier SLC25A12 gene on chromosome 2q31. Am J Psychiatry. 2005, 162: 2182-2184. 10.1176/appi.ajp.162.11.2182.
Ramoz N, Reichert JG, Smith CJ, Silverman JM, Bespalova IN, Davis KL, Buxbaum JD: Linkage and association of the mitochondrial aspartate/glutamate carrier SLC25A12 gene with autism. Am J Psychiatry. 2004, 161: 662-669. 10.1176/appi.ajp.161.4.662.
Silverman JM, Buxbaum JD, Ramoz N, Schmeidler J, Reichert JG, Hollander E, Angelo G, Smith CJ, Kryzak LA: Autism-related routines and rituals associated with a mitochondrial aspartate/glutamate carrier SLC25A12 polymorphism. Am J Med Genet B Neuropsychiatr Genet. 2008, 147: 408-410.
Gharani N, Benayed R, Mancuso V, Brzustowicz LM, Millonig JH: Association of the homeobox transcription factor, ENGRAILED 2, 3, with autism spectrum disorder. Mol Psychiatry. 2004, 9: 474-484. 10.1038/sj.mp.4001498.
Benayed R, Gharani N, Rossman I, Mancuso V, Lazar G, Kamdar S, Bruse SE, Tischfield S, Smith BJ, Zimmerman RA, Dicicco-Bloom E, Brzustowicz LM, Millonig JH: Support for the homeobox transcription factor gene ENGRAILED 2 as an autism spectrum disorder susceptibility locus. Am J Hum Genet. 2005, 77: 851-868. 10.1086/497705.
Philippi A, Tores F, Carayol J, Rousseau F, Letexier M, Roschmann E, Lindenbaum P, Benajjou A, Fontaine K, Vazart C, Gesnouin P, Brooks P, Hager J: Association of autism with polymorphisms in the paired-like homeodomain transcription factor 1 (PITX1) on chromosome 5q31: a candidate gene analysis. BMC Med Genet. 2007, 8: 74-10.1186/1471-2350-8-74.
Brooks P, Marcaillou C, Vanpeene M, Saraiva JP, Stockholm D, Francke S, Favis R, Cohen N, Rousseau F, Tores F, Lindenbaum P, Hager J, Philippi A: Robust physical methods that enrich genomic regions identical by descent for linkage studies: confirmation of a locus for osteogenesis imperfecta. BMC Genet. 2009, 10: 16-10.1186/1471-2156-10-16.
Carayol J, Letexier M, Tores F, Rousseau F, Sacco R, Persico AM, Hager J: Association study of a linkage region on chromosome 3P25 in 816 families shows strong evidence for association of ATP2B2 with autism and language delay. [http://imfar.confex.com/imfar/2009/webprogram/Paper4258.html]
Lauritsen MB, Als TD, Dahl HA, Flint TJ, Wang AG, Vang M, Kruse TA, Ewald H, Mors O: A genome-wide search for alleles and haplotypes associated with autism and related pervasive developmental disorders on the Faroe Islands. Mol Psychiatry. 2006, 11: 37-46. 10.1038/sj.mp.4001754.
McCauley JL, Li C, Jiang L, Olson LM, Crockett G, Gainer K, Folstein SE, Haines JL, Sutcliffe JS: Genome-wide and ordered-subset linkage analyses provide support for autism loci on 17 q and 19 p with evidence of phenotypic and interlocus genetic correlates. BMC Med Genet. 2005, 6: 1-10.1186/1471-2350-6-1.
Shao Y, Wolpert CM, Raiford KL, Menold MM, Donnelly SL, Ravan SA, Bass MP, McClain C, von WL, Vance JM, Abramson RH, Wright HH, Ashley-Koch A, Gilbert JR, DeLong RG, Cuccaro ML, Pericak-Vance MA: Genomic screen and follow-up analysis for autistic disorder. Am J Med Genet. 2002, 114: 99-105. 10.1002/ajmg.10153.
Lord C, Rutter M, Couteur A: Autism Diagnostic Interview-Revised: a revised version of a diagnostic interview for caregivers of individuals with possible pervasive developmental disorders. J Autism Dev Disord. 1994, 24: 659-685. 10.1007/BF02172145.
Lord C, Risi S, Lambrecht L, Cook EH, Leventhal BL, DiLavore PC, Pickles A, Rutter M: The autism diagnostic observation schedule-generic: a standard measure of social and communication deficits associated with the spectrum of autism. J Autism Dev Disord. 2000, 30: 205-223. 10.1023/A:1005592401947.
Schellenberg GD, Dawson G, Sung YJ, Estes A, Munson J, Rosenthal E, Rothstein J, Flodman P, Smith M, Coon H, Leong L, Yu CE, Stodgell C, Rodier PM, Spence MA, Minshew N, McMahon WM, Wijsman EM: Evidence for multiple loci from a genome scan of autism kindreds. Mol Psychiatry. 2006, 11: 1049-60. 10.1038/sj.mp.4001874.
Bolton PF, Pickles A, Murphy M, Rutter M: Autism, affective and other psychiatric disorders: patterns of familial aggregation. Psychol Med. 1998, 28: 385-395. 10.1017/S0033291797006004.
Dawson G, Estes A, Munson J, Schellenberg G, Bernier R, Abbott R: Quantitative assessment of autism symptom-related traits in probands and parents: Broader Phenotype Autism Symptom Scale. J Autism Dev Disord. 2007, 37: 523-536. 10.1007/s10803-006-0182-2.
Clayton DG: Prediction and interaction in complex disease genetics: Experience in type 1 diabetes. PLoS Genet. 2009, 5: e1000540-10.1371/journal.pgen.1000540.
Cauchi S, Meyre D, Durand E, Proenca C, Marre M, Hadjadj S, Choquet H, De GF, Gaget S, Allegaert F, Delplanque J, Permutt MA, Wasson J, Blech I, Charpentier G, Balkau B, Vergnaud AC, Czernichow S, Patsch W, Chikri M, Glaser B, Sladek R, Froguel P: Post genome-wide association studies of novel genes associated with type 2 diabetes show gene-gene interaction and high predictive value. PLoS ONE. 2008, 3: e2031-10.1371/journal.pone.0002031.
Cordell HJ, Barratt BJ, Clayton DG: Case/pseudocontrol analysis in genetic association studies: A unified framework for detection of genotype and haplotype associations, gene-gene and gene-environment interactions, and parent-of-origin effects. Genet Epidemiol. 2004, 26: 167-185. 10.1002/gepi.10307.
Cordell HJ: Properties of case/pseudocontrol analysis for genetic association studies: Effects of recombination, ascertainment, and multiple affected offspring. Genet Epidemiol. 2004, 26: 186-205. 10.1002/gepi.10306.
DeLong ER, Vernon WB, Bollinger RR: Sensitivity and specificity of a monitoring test. Biometrics. 1985, 41: 947-958. 10.2307/2530966.
Ziegler A, Blettner M, Kastner C, Chang-Claude J: Identifying influential families using regression diagnostics for generalized estimating equations. Genet Epidemiol. 1998, 15: 341-353. 10.1002/(SICI)1098-2272(1998)15:4<341::AID-GEPI2>3.0.CO;2-5.
Lepagnol-Bestel AM, Maussion G, Boda B, Cardona A, Iwayama Y, Delezoide AL, Moalic JM, Muller D, Dean B, Yoshikawa T, Gorwood P, Buxbaum JD, Ramoz N, Simonneau M: SLC25A12 expression is associated with neurite outgrowth and is upregulated in the prefrontal cortex of autistic subjects. Mol Psychiatry. 2008, 13: 385-397. 10.1038/sj.mp.4002120.
Yang P, Lung F, Jong Y, Hsieh H, Liang C, Juo S: Association of the homeobox transcription factor gene ENGRAILED 2 with autistic disorder in Chinese children. Neuropsychobiology. 2008, 57: 3-8. 10.1159/000123115.
Tripathi PP, Sgadò P, Scali M, Viaggi C, Casarosa S, Simon HH, Vaglini F, Corsini GU, Bozzi Y: Increased susceptibility to kainic acid-induced seizures in Engrailed-2 knockout mice. Neuroscience. 2009, 159: 842-849. 10.1016/j.neuroscience.2009.01.007.
Marinovic-Curin J, Marinovic-Terzic I, Bujas-Petkovic Z, Zekan L, Skrabic V, Dogas Z, Terzic J: Slower cortisol response during ACTH stimulation test in autistic children. Eur Child Adolesc Psychiatry. 2008, 17: 39-43. 10.1007/s00787-007-0632-1.
Curin J, Terzic J, Petkovic Z, Zekan L, Terzic I, Susnjara I: Lower cortisol and higher ACTH levels in individuals with autism. J Autism Dev Disord. 2003, 33: 443-448. 10.1023/A:1025019030121.
Hiroi N, Kino T, Bassett M, Rainey W, Phung M, Abu-Asab M, Fojo T, Briata P, Chorousos G, Bornstein S: Pituitary homeobox factor 1, a novel transcription factor in the adrenal regulating steroid 11ß-hydroxylase. Horm Metab Res. 2003, 35: 273-278. 10.1055/s-2003-41301.
Sandman CA, Touchette P, Marion S, Lenjavi M, Chicz-Demet A: Disregulation of proopiomelanocortin and contagious maladaptive behavior. Regulatory Peptides. 2002, 108: 179-185. 10.1016/S0167-0115(02)00097-6.
Rout UK, Dhossche DM: A pathogenetic model of autism involving Purkinje cell loss through anti-GAD antibodies. Medical Hypotheses. 2008, 71: 218-221. 10.1016/j.mehy.2007.11.012.
Meigs JB: Genotype score in addition to common risk factors for prediction of type 2 diabetes [abstract]. N Engl J Med. 2008, 359: 2208-10.1056/NEJMoa0804742.
Lyssenko V, Jonsson A, Almgren P, Pulizzi N, Isomaa B, Tuomi T, Berglund G, Altshuler D, Nilsson P, Groop L: Clinical risk factors, DNA variants, and the development of type 2 diabetes. N Engl J Med. 2008, 359: 2220-2232. 10.1056/NEJMoa0801869.
Fenske EC, Zalenski S, Krantz PJ, McClannahan LE: Age at intervention and treatment outcome for autistic children in a comprehensive intervention program. Analysis and Intervention in Developmental Disabilities. 1985, 5: 49-58. 10.1016/S0270-4684(85)80005-7.
Rogers SJ: Empirically supported comprehensive treatments for young children with autism. J Clin Child Psychol. 1998, 27: 138-145. 10.1207/s15374424jccp2702_1.
Pepe MS, Jane HE: Gauging the performance of SNPs, biomarkers, clinical factors for predicting risk of breast cancer. JNCI. 2008, 14: 978-979.
Abrahams BS, Geschwind DH: Advances in autism genetics: on the threshold of a new neurobiology. Nat Rev Genet. 2008, 9: 341-355. 10.1038/nrg2346.
Abrahams BS, Geschwind DH: Genetics of autism. Vogel and Motulsky's Human Genetics: Problems and Approaches. Edited by: Speicher M. 2009, New York, NY: Springer-Verlag, Chapter 23: 4
Lauritsen MB, Pedersen CB, Mortensen PB: Effects of familial risk factors and place of birth on the risk of autism: a nationwide register-based study. J Child Psychol Psychiatry. 2005, 46: 963-971. 10.1111/j.1469-7610.2004.00391.x.
Wang J, Ban MR, Zou GY, Cao H, Lin T, Kennedy BA, Anand S, Yusuf S, Huff MW, Pollex RL, Hegele RA: Polygenic determinants of severe hypertriglyceridemia. Hum Mol Genet. 2008, 17: 2894-2899. 10.1093/hmg/ddn188.
van Hoek M, Dehghan A, Witteman JCM, van Duijn CM, Uitterlinden AG, Oostra BA, Hofman A, Sijbrands EJG, Janssens AC: Predicting type 2 diabetes based on polymorphisms from genome-wide association studies. Diabetes. 2008, 57: 3122-3128. 10.2337/db08-0425.
Humphries SE, Cooper JA, Talmud PJ, Miller GJ: Candidate gene genotypes, along with conventional risk factor assessment, improve estimation of coronary heart disease risk in healthy UK men. Clin Chem. 2007, 53: 8-16. 10.1373/clinchem.2006.074591.
Wang K, Zhang H, Ma D, Bucan M, Glessner JT, Abrahams BS, Salyakina D, Imielinski M, Bradfield JP, Sleiman PM: Common genetic variants on 5p14.1 associate with autism spectrum disorders. Nature. 2009, 459: 528-533. 10.1038/nature07999.
Janssens AC, van Duijn CM: Genome-based prediction of common diseases: Methodological considerations for future research. Genome Med. 2009, 1: 20-10.1186/gm20.
We gratefully acknowledge the resources provided by the AGRE Consortium and the participating AGRE families. AGRE is a program of Autism Speaks, and is supported in part by grant 1U24MH081810 from the National Institute of Mental Health to Clara M. Lajonchere (PI). The University of Washington sample collection was funded by a grant from the NICHD (U19HD35465), which is part of the National Institutes of Health Collaborative Program of Excellence for Autism. We thank Drs Sommer L. Miller and Brett S. Abrahams and Patricia Lewin for their medical writing and editorial contributions. We also acknowledge Dr Inke R. König of the Institute of Medical Biometry and Statistics, University of Luebeck, Germany, for help with analyses of the Seattle sample data.
IntegraGen sponsored the design and statistical analysis of the AGRE sample analysis, and funded writing assistance in the form of preparation of the manuscript, references, figures, tables, formatting to journal style, and administrative support. The corresponding author had full access to the data in the study and final responsibility for interpretation of the data, and made the decision to submit for publication.
JC and FT are currently salaried employees of IntegraGen SA and JH is a former salaried employee. JC, FT and JH have stock options and patent applications with IntegraGen. GDS declares that he has no competing interests.
JC, FT, JH, AZ and GD conceived and designed the experiments. GDS and GD performed the experiments. JC and FT analyzed the AGRE sample data and AZ analyzed the Seattle sample data. JC, FT and GDS contributed reagents, materials and/or analysis tools, and JC, FT, JH, AZ and GD assisted with writing of the manuscript.
Electronic supplementary material
Additional file 1: Supplementary Table 1. Allelic variation used in the calculation of genetic score under an additive model. (PDF 8 KB)
Additional file 2: Supplementary Table 2. Odds ratios (ORs) and 95% confidence intervals (CIs) associated with the genetic score in subpopulations. (PDF 8 KB)
Authors’ original submitted files for images
Below are the links to the authors’ original submitted files for images.
Rights and permissions
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.
About this article
Cite this article
Carayol, J., Schellenberg, G.D., Tores, F. et al. Assessing the impact of a combined analysis of four common low-risk genetic variants on autism risk. Molecular Autism 1, 4 (2010). https://doi.org/10.1186/2040-2392-1-4
- Autism Spectrum Disorder
- Autism Spectrum Disorder
- Risk Allele
- Rett Syndrome
- Purkinje Neuron