- Open Access
Autism risk assessment in siblings of affected children using sex-specific genetic scores
Molecular Autismvolume 2, Article number: 17 (2011)
The inheritance pattern in most cases of autism is complex. The risk of autism is increased in siblings of children with autism and previous studies have indicated that the level of risk can be further identified by the accumulation of multiple susceptibility single nucleotide polymorphisms (SNPs) allowing for the identification of a higher-risk subgroup among siblings. As a result of the sex difference in the prevalence of autism, we explored the potential for identifying sex-specific autism susceptibility SNPs in siblings of children with autism and the ability to develop a sex-specific risk assessment genetic scoring system.
SNPs were chosen from genes known to be associated with autism. These markers were evaluated using an exploratory sample of 480 families from the Autism Genetic Resource Exchange (AGRE) repository. A reproducibility index (RI) was proposed and calculated in all children with autism and in males and females separately. Differing genetic scoring models were then constructed to develop a sex-specific genetic score model designed to identify individuals with a higher risk of autism. The ability of the genetic scores to identify high-risk children was then evaluated and replicated in an independent sample of 351 affected and 90 unaffected siblings from families with at least 1 child with autism.
We identified three risk SNPs that had a high RI in males, two SNPs with a high RI in females, and three SNPs with a high RI in both sexes. Using these results, genetic scoring models for males and females were developed which demonstrated a significant association with autism (P = 2.2 × 10-6 and 1.9 × 10-5, respectively).
Our results demonstrate that individual susceptibility associated SNPs for autism may have important differential sex effects. We also show that a sex-specific risk score based on the presence of multiple susceptibility associated SNPs allow for the identification of subgroups of siblings of children with autism who have a significantly higher risk of autism.
Autistic disorder is the most severe form of a group of autism spectrum disorders (ASDs) characterized by impairments in social interaction, deficits in verbal and non-verbal communication, restricted interests, and repetitive behaviors . With a prevalence of 1 in 110 children, ASDs are among the most common forms of severe developmental disability . The average recurrence risk of autism in siblings of affected children is approximately 10% . This rate is much higher than the prevalence rate for ASDs in the general population, but lower than would be expected for a highly penetrant mutation in a mendelian disorder .
The inheritance pattern of autism in most families is complex and not compatible with simple Mendelian inheritance [5, 6]. 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 [7, 8]. Several common variants localized in biological and positional (that is, under known linkage peaks) candidate genes have been associated with autism and some have been replicated in independent studies . Further support for these associations comes from genes for which, in addition to autism-associated common variants, rare mutations and/or copy number variations (CNVs) have been shown to contribute to the disease, and/or for which gene-disrupted mice exhibited autism-like traits. These genes include CNTNAP2 [10–13], RELN [14–19] and GABRB3 [20–23].
When taken individually, the risk of autism associated with variants remains modest, but Carayol et al.  recently showed that the accumulation of multiple risk alleles markedly increases the risk of autism in siblings of children who have been diagnosed with autism. They proposed a genetic score (GS) that, compared with studying polymorphisms individually, improves the identification of subgroups of individuals at greater risk of autism . In the case of autism, tools for genetic risk assessment are highly desirable to complement available behavioral assessments.
Another important characteristic of autism is the sex difference with a 4.5:1 male to female ratio . Second, intellectual disability, a key clinical dimension associated with outcome, is more frequent in females than males . Third, the risk of epilepsy is 18 times higher in females than males . This sex difference may partly be explained by sex-specific risk alleles or genes with different expression or activity based on sex [27, 28].
In the present study we propose to improve the genetic risk score model developed by Carayol et al.  by adding additional SNPs filtered for their relative importance using internal validation process and by also developing separate sex-specific genetic risk scores for males and females using a first sample of families with children with autism (exploratory sample). Their ability to better identify siblings of children with autism who are at high risk of autism was then evaluated and replicated in an independent second sample of autism families (replication sample).
The study design involved two independent family samples. The first sample (the 'exploratory' sample) consisted of 480 families from the Autism Genetic Resource Exchange (AGRE; http://www.agre.org) repository with at least 1 sibling diagnosed with a 'strict' definition of autism according to the Autism Diagnostic Interview Revisited (ADI-R) and no unaffected siblings. A total of 844 affected siblings including 664 males and 179 females met the diagnostic criteria for 'strict' autism. Minimizing phenotypic heterogeneity can lead to an improvement of the study power . Shao et al.  demonstrated that the use of homogeneous phenotype increases the power of linkage studies in autism. Linkage signals have been observed in studies in which the samples were stratified according to specific phenotypes such as the sex [28, 31, 32], delayed onset of phrase speech [30, 33, 34], and severe obsessive-compulsive behaviors . Two genome-wide association studies using overlapping samples of children with autism identified two different common variants in CNTNAP2, a gene localized in the 7q34-7q36 region linked to language disability in autism ; one SNP has been associated with autism through the use of the quantitative trait 'age at first word'  and the other using a qualitative strict autism diagnosis . Similarly, a recent genome-wide association study (GWAS)  reported the largest association with autism in MACROD2 using the strict autism diagnosis. Therefore, as in Shao et al. , we studied individuals with a strict autism rather than the heterogeneous broad autism spectrum disorder phenotype. The second sample (the 'replication' sample) included 187 families consisting of the 2 parents, at least 1 child with autism and 1 unaffected sibling from a sample collection at the University of Pennsylvania. This replication sample led to 351 children with autism (291 males and 60 females) with the same strict definition of the disease and 90 unaffected children (39 males and 51 females). 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 ethnicity. Caucasians represented the major ethnicity, with more than two-thirds of families in each sample.
Ten autism susceptibility genes were selected for this study. Four of them (PITX1, EN2, SLC25A12 and ATP2B2) have been previously demonstrated to have a predictive ability and were used in a genetic score-based model . Genes shown to be statistically associated with autism in at least one study using AGRE collection, even at the nominal level, and for which additional data support their implication in autism, were also included. Six genes fulfilled the statistical association condition, four of which were replicated in one or more independent study: HOXA1 [38, 39], GRIK2 [40–42], ITGB3 [43–46] and CNTNAP2 [10, 11]; one gene, MARK1, was found to be significantly overexpressed in brain from individuals with autism compared to unaffected individuals  and the last gene, JARID2 was chosen since one SNP, rs7766973, displays the strongest association with autism (P = 6.8 × 10-7 ) among the three GWAS performed on AGRE family data [37, 42, 48]. Table 1 lists the genes selected for the study and the associated SNPs with their deleterious alleles and corresponding frequencies.
All parents and children from the exploratory sample were genotyped for these ten markers. Only SNPs that were selected for further investigation were genotyped in the replication sample. 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 min, followed by 50 cycles at 92°C for 15 s and 60°C for 90 s (9700 Gene Amp PCR System; Applied Biosystems). Plates were then subjected to endpoint reading (7900 Real-Time PCR System; Applied Biosystems). The alleles were called automatically using the SDS software (Applied Biosystems), and a visual inspection of genotype clusters was performed. Genotyping quality was assessed by signal intensity plots and missing genotype frequencies; any sample with poor clustering and missing fractions ≥5% per SNP were retyped. Parental genotypes were used to investigate Hardy-Weinberg equilibrium (HWE) and to check for Mendelian inconsistencies. Families with remaining inconsistencies were excluded.
The development of the genetic score model and the definition of the increased risk GS thresholds (that define the high-risk groups) were based on the exploratory sample with all affected children whereas, for the replication study using the second sample, the index cases were excluded.
A model that is efficient only in the sample in which it was developed does not have validity. To be valid, the results need to be reproduced in a separate independent population. A genetic score model, such as the one proposed in this paper, is generally built on the simple sum of deleterious alleles observed at each of the chosen genes. Thus, the reproducibility of the genetic score is conditioned by the reproducibility of the deleterious allele for each SNPs included in the model. Markers that are more reproducible carry stronger and more stable information. The reproducibility of the SNPs was analyzed using the bootstrap resampling process and a reproducibility index (RI) was estimated similarly to Ma  as follows: (1) generation of a 'pseudosample' consisting of 480 families by randomly sampling the 480 families of the exploratory population with replacement; (2) estimation of the genetic relative risk associated with the deleterious allele of each SNP as defined in Table 1; (3) repetition 1,000 times of steps 1 and 2; (4) estimation for each SNP of the RIs indicating the proportion of 'pseudosamples' in which the deleterious allele maintains a risk greater than 1.00 in males, in females or in both males and females.
A high RI indicates that the effect of a deleterious allele of a given SNP is maintained across the bootstrap pseudosamples and that this SNP is a good candidate for the reproducibility of the genetic score. A stringent RI = 0.80 in children with autism was set to select best SNPs. Then, the RI in males and females with autism was checked separately to discard SNPs that lack of stability in a particular sex. Since all variants have been associated with autism using AGRE family data, this internal validation process prevents from an optimistic evaluation of their association, that is, an overestimation of the effect of risk alleles, and a potential deterioration of this effect in an independent sample. The sex genetic scores (GS) was then constructed as follows:
where sex = (male, female); RSall and RSsex are the risk scores built as the sum of deleterious alleles from genes with a high RI in males only (RSmale), in females only (RSfemale) or in both sexes (RSall); and Wall, Wmale, and Wfemale are the integer values of the corresponding genetic relative risks (GRR) associated with the corresponding risk scores (RSall, RSmale and RSfemale, respectively). These weights were calculated following Lin et al.  who showed that a weighted genetic score provided more predictive value than an unweighted genetic score.
Because the exploratory sample did not include unaffected children, all genetic relative risks were estimated as described in Carayol et al.  using the case-pseudocontrol approach proposed by Cordell and Clayton  and implemented in the DGCgenetics R package (http://www-gene.cimr.cam.ac.uk/clayton/software/). Sensitivity and specificity values of the GSs were estimated in the exploratory and the replication samples as in Carayol et al. . Areas under the receiver operating curves (AUCs) were estimated in the exploratory sample and tested against the AUC = 0.5 null hypothesis to validate the discriminative power of the GSs. However, AUCs do not provide an informative tool of the clinical utility of the genetic score (here, the high-risk classification of siblings of children with autism). Cutoff values were chosen to define a high-risk group in the exploratory sample and the odds ratios were estimated. These high-risk thresholds (one for male and one for female) were selected considering a false positive rate lower than 20% (that is, specificity higher than 80%). External validation of the clinical utility of the high-risk GS group was then conducted in the replication sample. Positive predictive values in siblings of children with autism were estimated from the sensitivity, specificity and the sibling recurrence risk estimates in males and females. Since no data were available in the literature, we estimated the sibling recurrence risk to 0.16 in males and 0.04 in females assuming an overall 0.10 sibling recurrence risk  and a 4:1 male to female sex ratio .
None of the SNPs exhibited a departure from HWE and allele frequencies were similar between samples (Table 1). Table 2 lists the RI of each SNP based on the bootstrap analysis using the exploratory sample. Eight markers reached the stringent 80% RI threshold. SNPs rs2292813 (SLC25A12) and rs2235076 (GRIK2) were excluded because of their low reproducibility (RI = 52% and 36%, respectively). Among the eight remaining SNPs, two displayed low RI in males but RI of 100% in females, rs12410279 (MARK1, RImale = 47%) and rs5918 (ITGB3, RImale = 65%). Inversely, three SNPs displayed a low RI in females and RI greater than 95% in males, rs227855 (ATP2B2, RIfemale = 59%), rs6872664 (PITX1, RIfemale = 30%) and rs10951154 (HOXA1, RIfemale = 20%).
The three separate risk scores were then constructed based on the sum of deleterious alleles in their corresponding SNPs. These included rs7794745, rs1861972 and rs7766973 for RSall, rs12410279 and rs5918 for RSfemale, and rs2278556, rs6872664 and rs10951154 for RSmale. The GRRs associated to one point increase in the RS were estimated to be 1.23 for RSall (P = 2.3 × 10-5; 95% confidence interval (CI) 1.12 to 1.36), 1.25 for RSmale (P = 5.8 × 10-4; 95% CI 1.10 to 1.41) and 2.29 for RSfemale (P = 1.7 × 10-6; 95% CI 1.57 to 3.34). The overall P value of the three tested scores were 3.1 × 10-9 with corresponding weights of 1.00, 1.00 and 2.00 for RSall, RSmale and RSfemale, respectively. The two genetic scores (GSs) were then constructed. GSmale ranged between 3 and 12 with a GRR associated to 1 point increase in the score of 1.23 (P = 2.2 × 10-6; 95% CI 1.13 to 1.34) and GSfemale ranged between 4 and 14 with a GRR of 1.41 (P = 1.9 × 10-5; 95% CI 1.21 to 1.65) for a highly significant global test with P = 8.4 × 10-10. Table 3 displays the sensitivity and specificity values for the GS in males and females. To define the high-risk group, GS values were selected in males and females with the aim to minimize the number of false positive below 20% and to maximize the sensitivity as high as possible. A genetic score threshold of nine points for males was associated with a moderate 0.24 sensitivity (95% CI 0.19 to 0.28) and a 0.86 specificity (95% CI 0.82 to 0.90) that minimizes the number of false positive test to 0.14 and lead to a 0.23 positive predictive value (PPV). For females, a genetic score threshold of 12 was associated with a similar specificity of 0.86 (95% CI 0.80 to 0.92) but a higher sensitivity of 0.37 (95% CI 0.29 to 0.44) and a PPV of 0.09. These two GS values were chosen as thresholds to define the group of children with a high risk of autism. AUCs were estimated to be 0.59 and 0.66 in males and females, respectively. They are both significantly different from the 0.5 null hypothesis (P = 2 × 10-8 and 1.5 × 10-7) indicating a predictive ability of the GSs.
In the replication sample (Table 4), sensitivity and specificity associated with the high-risk group GS threshold (GSmale = 9) were slightly higher in males (but not significantly different as it can be seen from the overlapping 95% CIs) with a 0.26 (95% CI 0.18 to 0.35) sensitivity and 0.87 (95% CI 0.76 to 0.98) specificity. The PPV reached 0.28 for a 0.16 sibling recurrence risk. Differences were observed in females for the sensitivity with an estimated 0.28 (95% CI 0.12 to 0.44) instead of 0.37 and the specificity with a 0.76 specificity (95% CI 0.64 to 0.89) instead of 0.86 but the differences were not significant (overlapping confidence intervals). In females, variances for sensitivity and specificity values were larger in the replication sample than in the exploratory sample because of the smaller number of females in the replication sample. As a consequence, the PPV (estimated to 5%) was very small and close to the 4% sibling recurrence risk.
Extending the analysis to a broader definition of autism and including or excluding the index cases as was performed with the replication study did not change the characteristics of the genetic score or the associated significance levels.
Our results demonstrate that the sex difference in autism may have an important influence on the genetic score characteristics, and therefore, on the risk assessment. Taking sex and reproducibility of the SNPs into account led to two GSs with different characteristics that allowed the identification of a subgroup of siblings of children with autism with a high risk of autism in males. The genetic score model with four genes  was also tested on this large sample of families and its association was clearly lower (P = 7 × 10-4 in males and females as a whole) compared to those of the sex-specific GSs (P = 2.2 × 10-6 and 1.9 × 10-5 for males and females, respectively). The risk for males with a high GS to develop autism was 28%, almost three times higher than the reported 10% sibling recurrence risk. In females, the 10% recurrence risk seems overestimated and we estimate this value to 4% considering a 4.5:1 male to female sex ratio.
The GS model has been developed through the use of affected children and the pseudocontrol approach [52, 53]. This was confirmed by analyzing unaffected siblings of children with autism. The pseudocontrols approach has been validated for the estimation of diagnostic accuracy using only affected children compared to full population-based data . We cannot exclude an over-representation of deleterious alleles in unaffected siblings compared to pseudocontrols, which are genetically the opposite of affected children, nor the effect of population controls that may lower the risk ratio between affected and unaffected siblings and consequently affect the discriminative ability of the GS models. This does not seem to occur for males since the high-risk class replicates its predictive accuracy but would need further investigation for females.
Reproducibility of effects is of major interest to enter in a predictive model since it conditions the reproducibility of the predictive model outside the study sample, which is of primary importance to validate such a model. According to the replication of the performance of the risk assessment model in males in an independent sample and the ability to find support for female specific variants despite the relatively small number of samples, the proposed approach can be used for developing stable and reproducible models. SLC25A12 associated and replicated in different studies [55–58] did not reach the reproducibility thresholds, whereas JARID2 that reached a suggestive significant threshold in a unique GWAS  seems of more interest. Some markers were reproducible (high RI) in a specific sex only but did not show any statistically significant interaction with sex nor were reported as being sex specific in the literature. The SNP rs7794745 located within CNTNAP2 has a high RI in both sexes whereas a previous association with autism has been reported preferentially in males [10, 11]. Due to the low number of females analyzed, these studies lack power to observe any association in females . Another SNP, rs5918 located within ITGB3, has been shown to be associated with autism in both sexes but with different risk effect , which could explain the difference of reproducibility observed in males and females. The stability is not necessarily linked to the sex specificity of the SNP or to the strength of previous association results. This may be explained in part by a study of Jakobsdotir et al.  which showed that a highly significant association of genes with a disease does not guarantee an effective discrimination between cases and controls.
Several limits of the study may be identified. The moderate number of females with autism in the replication sample as a consequence of the significant sex ratio in autism led to a lack of power for the replication of the high-risk group characteristics. Sibling recurrence risk of males and females were not estimated or reported from real data but calculated assuming a sibling recurrence risk of 10%  and the widely observed 4.5:1 male to female sex ratio. Reported PPVs are intuitive estimates that quantify the increase in the risk for an individual (a sibling of a child with autism) who has a genetic score that falls in the high-risk class. Accurate PPVs could be estimated by using observed and reported data. The selection of the genes and the SNPs included in the genetic scores could be discussed. The methodology used to select the common variants and the internal validation approach performed in this study strongly support the implication of these SNPs in autism as well as their discriminative ability. The addition of other SNPs from the same genetic region would have led to a much more complicated model because of the linkage disequilibrium (LD) between these SNPs as well as the haplotypes resulting from the different combination of alleles. Finally, other approaches may be used to select genes to enter in a genetic score. Genes may be selected using statistically significant results from GWAS [60, 61] or a complementary approach as in convergent functional genomics (CFG) autism [62, 63], when none or few association results reach significance as it is frequently the case in complex disease and particularly in autism.
The recent paper of Lu and Cantor  together with the present results highlights the importance of the sex in genetic study of autism. They showed that using sex as a risk factor in GWAS of multiplex autism families increased the power of the study and identified one new gene implicated in calcium channel defect. Stone et al.  also suggested that sex is an important factor in the genetics of autism and could be used to decrease heterogeneity in genetic study.
The results of this study confirm previous results  that predictive models are of major interest in autism and may help to identify siblings of children with autism at high risk of disease. The choice of genes to enter in the model must be made with caution since association and replication of a particular SNP in different studies are not sufficient justification to enter a SNP in a genetic score and sex is an important factor that needs to be included in autism risk evaluation.
Johnson CP, Myers SM: Identification and evaluation of children with autism spectrum disorders. Pediatrics. 2007, 120: 1183-1215. 10.1542/peds.2007-2361.
Rice CE: Prevalence of autism spectrum disorders - Autism and Developmental Disabilities Monitoring Network, United States, 2006. MMWR Surveill Summ. 2009, 58: 1-20.
Constantino JN, Zhang Y, Frazier T, Abbacchi AM, Law P: Sibling recurrence and the genetic epidemiology of autism. Am J Psychiatry. 2010, 167: 1349-1356. 10.1176/appi.ajp.2010.09101470.
Muhle R, Trentacoste SV, Rapin I: The genetics of autism. Pediatrics. 2004, 113: e472-e486. 10.1542/peds.113.5.e472.
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.
Autism Genome Project Consortium, Szatmari P, Paterson AD, Zwaigenbaum L, Roberts W, Brian J, Liu XQ, Vincent JB, Skaug JL, Thompson AP, Senman L, Feuk L, Qian C, Bryson SE, Jones MB, Marshall CR, Scherer SW, Vieland VJ, Bartlett C, Mangin LV, Goedken R, Segre A, Pericak-Vance MA, Cuccaro ML, Gilbert JR, Wright HH, Abramson RK, Betancur C, Bourgeron T, Gillberg C: Mapping autism risk loci using genetic linkage and chromosomal rearrangements. Nat Genet. 2007, 39: 319-328. 10.1038/ng1985.
Fenske E, Zalenski S, Krants P, McClannahand L: Age at intervention and treatment outcome for autistic children in a comprehensive intervention program. Special issue: early intervention. Anal Interven Devel. 1985, 5: 49-58.
Rogers SJ: Empirically supported comprehensive treatments for young children with autism. J Clin Child Psychol. 1998, 27: 168-179. 10.1207/s15374424jccp2702_4.
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.
Alarcón M, Abrahams BS, Stone JL, Duvall JA, Perederiy JV, Bomar JM, Sebat J, Wigler M, Martin CL, Ledbetter DH: Linkage, association, and gene-expression analyses identify CNTNAP2 as an autism-susceptibility gene. Am J Hum Genet. 2008, 82: 150-159. 10.1016/j.ajhg.2007.09.005.
Arking DE, Cutler DJ, Brune CW, Teslovich TM, West K, Ikeda M, Rea A, Guy M, Lin S, Cook EH, Chakravarti A: A common genetic variant in the neurexin superfamily member CNTNAP2 increases familial risk of autism. Am J Hum Genet. 2008, 82: 160-164. 10.1016/j.ajhg.2007.09.015.
Bakkaloglu B, O'Roak BJ, Louvi A, Gupta AR, Abelson JF, Morgan TM, Chawarska K, Klin A, Ercan-Sencicek AG, Stillman AA, Tanriover G, Abrahams BS, Duvall JA, Robbins EM, Geschwind DH, Biederer T, Gunel M, Lifton RP, State MW: Molecular cytogenetic analysis and resequencing of contactin associated protein-like 2 in autism spectrum disorders. Am J Hum Genet. 2008, 82: 165-173. 10.1016/j.ajhg.2007.09.017.
Poot M, Beyer V, Schwaab I, Damatova N, Slot R, Prothero J, Holder SE, Haaf T: Disruption of CNTNAP2 and additional structural genome changes in a boy with speech delay and autism spectrum disorder. Neurogenetics. 2009, 11: 81-89.
Bonora E, Beyer KS, Lamb JA, Parr JR, Klauck SM, Benner A, Paolucci M, Abbott A, Ragoussis I, Poustka A, Bailey AJ, Monaco AP, International Molecular Genetic Study of Autism (IMGSAC): Analysis of reelin as a candidate gene for autism. Mol Psychiatry. 2003, 8: 885-892. 10.1038/sj.mp.4001310.
Li H, Li Y, Shao J, Li R, Qin Y, Xie C, Zhao Z: The association analysis of RELN and GRM8 genes with autistic spectrum disorder in Chinese Han population. Am J Med Genet B Neuropsychiatr Genet. 2008, 147B: 194-200. 10.1002/ajmg.b.30584.
Persico AM, D'Agruma L, Maiorano N, Totaro A, Militerni R, Bravaccio C, Wassink TH, Schneider C, Melmed R, Trillo S, Montecchi F, Palermo M, Pascucci T, Puglisi-Allegra S, Reichelt KL, Conciatori M, Marino R, Quattrocchi CC, Baldi A, Zelante L, Gasparini P, Keller F, Collaborative Linkage Study of Autism: Reelin gene alleles and haplotypes as a factor predisposing to autistic disorder. Mol Psychiatry. 2001, 6: 150-159. 10.1038/sj.mp.4000850.
Salinger WL, Ladrow P, Wheeler C: Behavioral phenotype of the reeler mutant mouse: effects of RELN gene dosage and social isolation. Behav Neurosci. 2003, 117: 1257-1275.
Serajee FJ, Zhong H, Mahbubul Huq AH: Association of Reelin gene polymorphisms with autism. Genomics. 2006, 87: 75-83. 10.1016/j.ygeno.2005.09.008.
Skaar DA, Shao Y, Haines JL, Stenger JE, Jaworski J, Martin ER, DeLong GR, Moore JH, McCauley JL, Sutcliffe JS, Ashley-Koch AE, Cuccaro ML, Folstein SE, Gilbert JR, Pericak-Vance MA: Analysis of the RELN gene as a genetic risk factor for autism. Mol Psychiatry. 2005, 10: 563-571. 10.1038/sj.mp.4001614.
Buxbaum JD, Silverman JM, Smith CJ, Greenberg DA, Kilifarski M, Reichert J, Cook EH, Fang Y, Song CY, Vitale R: Association between a GABRB3 polymorphism and autism. Mol Psychiatry. 2002, 7: 311-316. 10.1038/sj.mp.4001011.
Cook EH, Courchesne RY, Cox NJ, Lord C, Gonen D, Guter SJ, Lincoln A, Nix K, Haas R, Leventhal BL, Courchesne E: Linkage-disequilibrium mapping of autistic disorder, with 15q11-13 markers. Am J Hum Genet. 1998, 62: 1077-1083. 10.1086/301832.
DeLorey TM, Sahbaie P, Hashemi E, Homanics GE, Clark JD: Gabrb3 gene deficient mice exhibit impaired social and exploratory behaviors, deficits in non-selective attention and hypoplasia of cerebellar vermal lobules: a potential model of autism spectrum disorder. Behav Brain Res. 2008, 187: 207-220. 10.1016/j.bbr.2007.09.009.
Samaco RC, Hogart A, LaSalle JM: Epigenetic overlap in autism-spectrum neurodevelopmental disorders: MECP2 deficiency causes reduced expression of UBE3A and GABRB3. Hum Mol Genet. 2005, 14: 483-492.
Carayol J, Schellenberg GD, Tores F, Hager J, Ziegler A, Dawson G: Assessing the impact of a combined analysis of four common low-risk genetic variants on autism risk. Molecular Autism. 2010, 1: 4-10.1186/2040-2392-1-4.
Gillberg C, Steffenburg S, Schaumann H: Is autism more common now than ten years ago?. Br J Psychiatry. 1991, 158: 403-409. 10.1192/bjp.158.3.403.
Amiet C, Gourfinkel-An I, Bouzamondo A, Tordjman S, Baulac M, Lechat P, Mottron L, Cohen D: Epilepsy in autism is associated with intellectual disability and gender: evidence from a meta-analysis. Biol Psychiatry. 2008, 64: 577-582. 10.1016/j.biopsych.2008.04.030.
Harrison PJ, Tunbridge EM: Catechol-O-methyltransferase (COMT): a gene contributing to sex differences in brain function, and to sexual dimorphism in the predisposition to psychiatric disorders. Neuropsychopharmacology. 2008, 33: 3037-3045. 10.1038/sj.npp.1301543.
Stone JL, Merriman B, Cantor RM, Yonan AL, Gilliam TC, Geschwind DH, Nelson SF: Evidence for sex-specific risk alleles in autism spectrum disorder. Am J Hum Genet. 2004, 75: 1117-1123. 10.1086/426034.
McCarthy MI, Abecasis GR, Cardon LR, Goldstein DB, Little J, Ioannidis JP, Hirschhorn JN: Genome-wide association studies for complex traits: consensus, uncertainty and challenges. Nat Rev Genet. 2008, 9: 356-369. 10.1038/nrg2344.
Shao Y, Raiford KL, Wolpert CM, Cope HA, Ravan SA, Ashley-Koch AA, Abramson RK, Wright HH, DeLong RG, Gilbert JR, Cuccaro ML, Pericak-Vance MA: Phenotypic homogeneity provides increased support for linkage on chromosome 2 in autistic disorder. Am J Hum Genet. 2002, 70: 1058-1061. 10.1086/339765.
Cantor RM, Kono N, Duvall JA, Alvarez-Retuerto A, Stone JL, Alarcon M, Nelson SF, Geschwind DH: Replication of autism linkage: fine-mapping peak at 17q21. Am J Hum Genet. 2005, 76: 1050-1056. 10.1086/430278.
Lamb JA, Barnby G, Bonora E, Sykes N, Bacchelli E, Blasi F, Maestrini E, Broxholme J, Tzenova J, Weeks D, Bailey AJ, Monaco AP, International Molecular Genetic Study of Autism Consortium (IMGSAC): Analysis of IMGSAC autism susceptibility loci: evidence for sex limited and parent of origin specific effects. J Med Genet. 2005, 42: 132-137. 10.1136/jmg.2004.025668.
Buxbaum JD, Silverman JM, Smith CJ, Kilifarski M, Reichert J, Hollander E, Lawlor BA, Fitzgerald M, Greenberg DA, Davis KL: Evidence for a susceptibility gene for autism on chromosome 2 and for genetic heterogeneity. Am J Hum Genet. 2001, 68: 1514-1520. 10.1086/320588.
Shao Y, Wolpert CM, Raiford KL, Menold MM, Donnelly SL, Ravan SA, Bass MP, McClain C, von Wendt L, 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.
Buxbaum JD, Silverman J, Keddache M, Smith CJ, Hollander E, Ramoz N, Reichert JG: Linkage analysis for autism in a subset families with obsessive-compulsive behaviors: evidence for an autism susceptibility gene on chromosome 1 and further support for susceptibility genes on chromosome 6 and 19. Mol Psychiatry. 2004, 9: 144-150. 10.1038/sj.mp.4001465.
Alarcon M, Cantor RM, Liu J, Gilliam TC, Geschwind DH: Evidence for a language quantitative trait locus on chromosome 7q in multiplex autism families. Am J Hum Genet. 2002, 70: 60-71. 10.1086/338241.
Anney R, Klei L, Pinto D, Regan R, Conroy J, Magalhaes TR, Correia C, Abrahams BS, Sykes N, Pagnamenta AT, Almeida J, Bacchelli E, Bailey AJ, Baird G, Battaglia A, Berney T, Bolshakova N, Bölte S, Bolton PF, Bourgeron T, Brennan S, Brian J, Carson AR, Casallo G, Casey J, Chu SH, Cochrane L, Corsello C, Crawford EL, Crossett A: A genome-wide scan for common alleles affecting risk for autism. Hum Mol Genet. 2010, 19: 4072-4082. 10.1093/hmg/ddq307.
Conciatori M, Stodgell CJ, Hyman SL, O'Bara M, Militerni R, Bravaccio C, Trillo S, Montecchi F, Schneider C, Melmed R, Elia M, Crawford L, Spence SJ, Muscarella L, Guarnieri V, D'Agruma L, Quattrone A, Zelante L, Rabinowitz D, Pascucci T, Puglisi-Allegra S, Reichelt KL, Rodier PM, Persico AM: Association between the HOXA1 A218G polymorphism and increased head circumference in patients with autism. Biol Psychiatry. 2004, 55: 413-419. 10.1016/j.biopsych.2003.10.005.
Ingram JL, Stodgell CJ, Hyman SL, Figlewicz DA, Weitkamp LR, Rodier PM: Discovery of allelic variants of HOXA1 and HOXB1: genetic susceptibility to autism spectrum disorders. Teratology. 2000, 62: 393-405. 10.1002/1096-9926(200012)62:6<393::AID-TERA6>3.0.CO;2-V.
Dutta S, Das S, Guhathakurta S, Sen B, Sinha S, Chatterjee A, Ghosh S, Ahmed S, Usha R: Glutamate receptor 6 gene (GluR6 or GRIK2) polymorphisms in the Indian population: a genetic association study on autism spectrum disorder. Cell Mol Neurobiol. 2007, 27: 1035-1047. 10.1007/s10571-007-9193-6.
Kim SA, Kim JH, Park M, Cho IH, Yoo HJ: Family-based association study between GRIK2 polymorphisms and autism spectrum disorders in the Korean trios. Neurosci Res. 2007, 58: 332-335. 10.1016/j.neures.2007.03.002.
Wang K, Zhang H, Ma D, Bucan M, Glessner JT, Abrahams BS, Salyakina D, Imielinski M, Bradfield JP, Sleiman PM, Kim CE, Hou C, Frackelton E, Chiavacci R, Takahashi N, Sakurai T, Rappaport E, Lajonchere CM, Munson J, Estes A, Korvatska O, Piven J, Sonnenblick LI, Alvarez Retuerto AI, Herman EI, Dong H, Hutman T, Sigman M, Ozonoff S, Klin A: Common genetic variants on 5p14.1 associate with autism spectrum disorders. Nature. 2009, 459: 528-533. 10.1038/nature07999.
Coutinho AM, Sousa I, Martins M, Correia C, Morgadinho T, Bento C, Marques C, Ataíde A, Miguel TS, Moore JH, Oliveira G, Vicente AM: Evidence for epistasis between SLC6A4 and ITGB3 in autism etiology and in the determination of platelet serotonin levels. Hum Genet. 2007, 121: 243-256. 10.1007/s00439-006-0301-3.
Ma DQ, Rabionet R, Konidari I, Jaworski J, Cukier HN, Wright HH, Abramson RK, Gilbert JR, Cuccaro ML, Pericak-Vance MA, Martin ER: Association and gene-gene interaction of SLC6A4 and ITGB3 in autism. Am J Med Genet B Neuropsychiatr Genet. 2010, 153B: 477-483.
Napolioni V, Lombardi F, Sacco R, Curatolo P, Manzi B, Alessandrelli R, Militerni R, Bravaccio C, Lenti C, Saccani M, Schneider C, Melmed R, Pascucci T, Puglisi-Allegra S, Reichelt KL, Rousseau F, Lewin P, Persico AM: Family-based association study of ITGB3 in autism spectrum disorder and its endophenotypes. Eur J Hum Genet. 2011, 19: 353-359. 10.1038/ejhg.2010.180.
Weiss LA, Kosova G, Delahanty RJ, Jiang L, Cook EH, Ober C, Sutcliffe JS: Variation in ITGB3 is associated with whole-blood serotonin level and autism susceptibility. Eur J Hum Genet. 2006, 14: 923-931. 10.1038/sj.ejhg.5201644.
Maussion G, Carayol J, Lepagnol-Bestel AM, Tores F, Loe-Mie Y, Milbreta U, Rousseau F, Fontaine K, Renaud J, Moalic JM, Philippi A, Chedotal A, Gorwood P, Ramoz N, Hager J, Simonneau M: Convergent evidence identifying MAP/microtubule affinity-regulating kinase 1 (MARK1) as a susceptibility gene for autism. Hum Mol Genet. 2008, 17: 2541-2551. 10.1093/hmg/ddn154.
Weiss LA, Arking DE, Daly MJ, Chakravarti A: A genome-wide linkage and association scan reveals novel loci for autism. Nature. 2009, 461: 802-808. 10.1038/nature08490.
Ma S: Empirical study of supervised gene screening. BMC Bioinformatics. 2006, 7: 537-10.1186/1471-2105-7-537.
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.
Cordell HJ, Clayton DG: A unified stepwise regression procedure for evaluating the relative effects of polymorphisms within a gene using case/control or family data: application to HLA in type 1 diabetes. Am J Hum Genet. 2002, 70: 124-141. 10.1086/338007.
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.
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.
Carayol J, Tores F, Konig IR, Hager J, Ziegler A: Evaluating diagnostic accuracy of genetic profiles in affected offspring families. Stat Med. 2010, 29: 2359-2368. 10.1002/sim.4006.
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.
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.
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.
Silverman JM, Buxbaum JD, Ramoz N, Schmeidler J, Reichenberg A, 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.
Jakobsdottir J, Gorin MB, Conley YP, Ferrell RE, Weeks DE: Interpretation of genetic association studies: markers with replicated highly significant odds ratios may be poor classifiers. PLoS Genet. 2009, 5: e1000337-10.1371/journal.pgen.1000337.
De Jager PL, Chibnik LB, Cui J, Reischl J, Lehr S, Simon KC, Aubin C, Bauer D, Heubach JF, Sandbrink R, Tyblova M, Lelkova P, Steering committee of the BENEFIT study; Steering committee of the BEYOND study; Steering committee of the LTF study, Steering committee of the CCR1 study, Havrdova E, Pohl C, Horakova D, Ascherio A, Hafler DA, Karlson EW: Integration of genetic risk factors into a clinical algorithm for multiple sclerosis susceptibility: a weighted genetic risk score. Lancet Neurol. 2009, 8: 1111-1119. 10.1016/S1474-4422(09)70275-3.
Manolio TA: Genomewide association studies and assessment of the risk of disease. N Engl J Med. 2010, 363: 166-176. 10.1056/NEJMra0905980.
Le-Niculescu H, Patel SD, Bhat M, Kuczenski R, Faraone SV, Tsuang MT, McMahon FJ, Schork NJ, Nurnberger JI, Niculescu AB: Convergent functional genomics of genome-wide association data for bipolar disorder: comprehensive identification of candidate genes, pathways and mechanisms. Am J Med Genet B Neuropsychiatr Genet. 2009, 150B: 155-181. 10.1002/ajmg.b.30887.
Patel SD, Le-Niculescu H, Koller DL, Green SD, Lahiri DK, McMahon FJ, Nurnberger JI, Niculescu AB: Coming to grips with complex disorders: genetic risk prediction in bipolar disorder using panels of genes identified through convergent functional genomics. Am J Med Genet B Neuropsychiatr Genet. 2010, 153B: 850-877.
Lu AT, Cantor RM: Allowing for sex differences increases power in a GWAS of multiplex Autism families. Mol Psychiatry. 2010
Wigginton JE, Cutler DJ, Abecasis GR: A note on exact tests of Hardy-Weinberg equilibrium. Am J Hum Genet. 2005, 76: 887-893. 10.1086/429864.
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 Pennsylvania sample collection was funded by UW Autism Center for Excellence grant # 5-P50-HD055782 and a grant from Autism Speaks. We thank Dr Thomas Rio Frio and Dr Brett S Abrahams for their helpful critical review of the manuscript. 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, 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 FR are currently salaried employees of IntegraGen SA and have stock options and patent applications with IntegraGen. GDS, BD and GD declare that they have no competing interests. EG is a consultant for IntegraGen SA.
JC, FR and GDS conceived and designed the experiments. FR, BD and GDS performed the experiments. JC analyzed the data and draft the manuscript. EG validated the statistical method. JC, FR and GDS contributed reagents, materials and/or analysis tools. GD and GS contributed to the collection of the University of Pennsylvanian sample. All coauthors assisted with writing of the manuscript. All authors read and approved the final manuscript.