- 1Centre for Genetic Origins of Health and Disease, University of Western Australia, Perth, WA, Australia
- 2School of Chemistry and Biochemistry, University of Western Australia, Perth, WA, Australia
- 3Telethon Institute for Child Health Research, Perth, WA, Australia
Lay abstract: It has been proposed that autistic-like traits in the general population lie on a continuum, with clinical Autism Spectrum Disorder (ASD), representing the extreme end of this distribution. The current study undertook a genome-wide association (GWA) scan of 965 young Western Australian adults to identify novel risk variants associated with autistic-like traits. No associations reached genome-wide significance; however, a review of nominally associated single nucleotide polymorphisms (SNPs) indicated two positional candidate loci that have been previously implicated in autistic-like trait etiology.
Scientific abstract: Research has proposed that autistic-like traits in the general population lie on a continuum, with clinical ASD representing the extreme end of this distribution. Inherent in this proposal is that biological mechanisms associated with clinical ASD may also underpin variation in autistic-like traits within the general population. A GWA study using 2,462,046 SNPs was undertaken for ASD in 965 individuals from the Western Australian Pregnancy Cohort (Raine) Study. No SNP associations reached genome-wide significance (p < 5.0 × 10−8). However, investigations into nominal observed SNP associations (p < 1.0 × 10−5) add support to two positional candidate genes previously implicated in ASD etiology, PRKCB1, and CBLN1. The rs198198 SNP (p = 9.587 × 10−6), is located within an intron of the protein kinase C, beta 1 (PRKCB1) gene on chromosome 16p11. The PRKCB1 gene has been previously reported in linkage and association studies for ASD, and its mRNA expression has been shown to be significantly down regulated in ASD cases compared with controls. The rs16946931 SNP (p = 1.78 × 10−6) is located in a region flanking the Cerebellin 1 (CBLN1) gene on chromosome 16q12.1. The CBLN1 gene is involved with synaptogenesis and is part of a gene family previously implicated in ASD. This GWA study is only the second to examine SNPs associated with autistic-like traits in the general population, and provides evidence to support roles for the PRKCB1 and CBLN1 genes in risk of clinical ASD.
Introduction
Autism Spectrum Disorders (ASD) represent a group of neurodevelopmental conditions characterized by impairments in social interaction and communication, and repetitive interests and behaviors. A recent review estimated the median global prevalence of ASD at 62 cases per 10,000 children (Elsabbagh et al., 2012). The overall heritability of ASD is estimated at 90%, which is amongst the highest for any neuropsychiatric disorder (Hollander et al., 2003; Skuse et al., 2005). Evidence suggests that ASD have a complex inheritance where different, likely overlapping, groups of genetic variants may cause susceptibility to disease (Veenstra-VanderWeele et al., 2004). Increasing numbers of epidemiological and genetic studies have deepened the understanding of genetic contribution to ASD, and show that a variety of genetic mechanisms may be involved in the etiology (Li et al., 2012).
Autistic-like traits are sub-threshold deficits in socialization, communication, and restricted interests that do not meet formal criteria for ASD (Constantino and Todd, 2003). Several authors have suggested that ASD can be conceptualized as the conditions arising in individuals found at the extreme end of a normal distribution of autistic-like traits (Gillberg, 1992; Constantino and Todd, 2003; Ronald et al., 2005). Population-based studies have supported this by finding that, in addition to individuals with ASD, many others exhibit sub-threshold autistic or autistic-like traits (Constantino and Todd, 2003; Posserud et al., 2006; Lundström et al., 2012). Sub-threshold autistic-like traits have been examined in general population twin studies, with heritability estimates ranging from 36 to 87% in these studies (Ronald and Hoekstra, 2011).
With the discovery of common genetic variants associated with ASD, the question is raised whether common risk loci also contribute to variation in phenotypes such as autistic-like traits (St Pourcain et al., 2010). One methodological approach that is gaining influence in ASD research is the examination of the quantitative distribution of autistic-like traits in the broader community, rather than a dichotomy of ASD cases and controls. Many authors have postulated that common genetic variants that are present in a significant proportion of the general population may play a role in the etiology of ASD (Campbell et al., 2006; Alarcón et al., 2008; Chakrabarti et al., 2009; Wang et al., 2009; Anney et al., 2010; Ronald et al., 2010). Further to this, it is thought that understanding the etiology of individual differences in autistic-like traits in the general population will help in understanding the causes of ASD (Ronald and Hoekstra, 2011).
The aim of the current study was to use genome-wide association (GWA) methods to search for novel risk variants associated with autistic-like traits in a Western Australian population sample.
Materials and Methods
Ethics Statement
Participant recruitment of the study families were approved by the Human Ethics Committee at King Edward Memorial Hospital. The 20 year follow-up ethics approval was received from the Human Research Ethics Committee at the University of Western Australia. Participants provided written informed consent for data collection of autistic-like trait outcomes at approximately 20 years of age.
GWA Study Sample Population
Data were obtained from the Western Australian Pregnancy Cohort (Raine) Study based at the Telethon Institute for Child Health Research. The Raine Study is a longitudinal investigation of pregnant women and their offspring, who were recruited from King Edward Memorial Hospital, Perth, Western Australia or nearby private practices between 1989 and 1991 (Newnham et al., 1993). From the 2,900 pregnancies recruited into the Raine Study, 2,868 live-born children have been followed since the commencement of the study. The final sample consisted of 965 individuals from the Raine study with both genotype and outcome measures.
Measure of Autistic-Like Traits
At the 20 year follow-up, Raine Study participants who did not have a diagnosis of intellectual disability or ASD were asked to complete the Autism Spectrum Quotient (AQ) (Baron-Cohen et al., 2001). The AQ is a self-report questionnaire that provides a quantitative measure of autistic-like traits in the general population (Baron-Cohen et al., 2001). Individuals are provided with 50 statements and asked to indicate on a 4-point scale how well that statement applies to them (strongly agree, agree, disagree, strongly disagree). The items were scored on a scale ranging from 1 to 4, based on previous research that this scoring method retains more information about responses than the dichotomous scoring first proposed for this instrument (Baron-Cohen et al., 2001; Austin, 2005; Stewart and Austin, 2009; Russell-Smith et al., 2011). Scores on each item are summed to provide a Total AQ, with higher scores indicating greater autistic-like traits. The Total AQ is known to have good test-retest reliability (r = 0.7), and validation studies have found that scores in the general population follow a normal quantitative distribution (Baron-Cohen et al., 2001; Whitehouse et al., 2011). Factor analyses of the AQ in several countries have consistently identified three clear factors, relating to social ability, attention to detail/patterns, and the understanding of others. In the current study, we divided items into the subscales identified in a study of Western Australian adults (Russell-Smith et al., 2011), who were highly similar to the sample under investigation here: Social Skills, Details/Patterns, and Communication/Mindreading. There is minimal difference between the items pertaining to these subscales, and those reported in other factor analyses. For the current data set, internal reliability for the scales ranged from moderate (Communication/Mindreading: α = 0.63) to good (Details/Patterns: α = 0.78) and excellent (Social Skills: α = 0.85). In this study, Total AQ scores and the scores of the three subscales were used as four outcome measures.
GWA Study Genotyping
In the Raine Study, DNA was collected using standardized procedures from 74% of adolescents who attended the 14 year follow-up, and a further 5% at the 16 year follow-up. Of the 1,593 offspring for whom genome-wide data were available, 99 individuals were excluded for the following reasons: 16 had low genotyping success ( >3% missing); 3 had excessive heterozygosity, 68 were in high identical-by-descent with another Raine Study participant (related), 7 whose sex was ambiguous, and 5 individuals had mislabeled samples. There were 1,494 individuals whose DNA sample passed quality control criteria and were eligible for the genetic analyses (768 males, 726 females), see Figure 1.
Genome-wide data were generated using an Illumina Human 660W Quad array at the Center for Applied Genomics (Toronto, ON, Canada). The 660W Quad Array includes 657,366 genetic variants including ∼560,000 single nucleotide polymorphisms (SNPs) and ∼95,000 copy number variants. SNPs were excluded based on the Wellcome Trust Case Control Consortium thresholds:
HWE p-value must be >5.7 × 10−7 (919 markers excluded)
Call rate must be >0.95 (95%) (97,718 markers excluded)
MAF must be >0.01 (1%) (119,246 markers excluded – includes copy number variants)
A/T and G/C SNPs were removed due to possible strand ambiguity.
In total, 535,632 SNPs passed genotype quality control before genotype imputation. Imputation was performed using MACH (Li and Abecasis, 2006) software on the 22 autosomes. Once the data were cleaned, a two-step process was carried out using the CEU samples from HapMap phase2 build 36 release 22 as a reference panel. This imputation process imputed both individual genotypes for SNPs that were typed but not called for a particular individual, and SNPs that were not typed in all individuals. After imputation, genotype information was available for 2,543,887 SNPs.
Data Analysis
Power calculations
Power calculations were performed using the statistical software package Quanto (Gauderman, 2002) for a sample size of 965 and MAF ranging from 0 to 0.5. The means and standard deviations for Total AQ and each subscale were used for the power calculations. Power was calculated to detect a minimum change of one unit assuming an additive genetic model. An alpha level of 1.0 × 10−8 was used for the GWA studies.
Methods for genome-wide association
Allelic dosage scores generated from MACH were used in analyses. The GWA studies were run adjusting for both age and sex. To account for potential population substructure, first and second principal components were also included in the multivariate GWA model. According to the Eigenvectors the first principal component accounted for 6.5% of variation, and the second component accounted for 2.3% of variation. Analysis was run using the ProbABEL package (Aulchenko et al., 2010). The combined results and SNP details file were then filtered to remove SNPs with MAF <1% and an imputation quality value of <0.3 (suggesting poor imputation quality).
A quantile–quantile (Q–Q) plot depicting –log10 transformations of the observed p-values as a function of the expected p-values was generated using “qqunif” function of the R Genetic Analysis Package (GAP) (R Development Core Team, 2008). The genomic inflation factor, λ, was calculated from the observed p-values. A Manhattan plot displaying a –log10 transformation of the observed p-values was generated using “mhtplot” function of the R GAP (R Development Core Team, 2008).
Genome-wide significance of results was defined as a p-value ≤ 5.0 × 10−8 (Panagiotou and Ioannidis, 2011). Tables of SNPs were generated, based on a p-value < 1.0 × 10−5, which has been used by the National Human Genome Research Institute for the identification and archiving of suggestive associations (National Human Genome Research Institute, 2012). Information regarding genes and function for SNPs included in these tables were sourced from SNP Annotation and Proxy Search (SNAP, 2008). Genes and/or regions that appeared in the generated tables were searched using the Simon’s Foundation Autism Research Initiative (SFARI) database (Simons Foundation Autism Research Initiative, 2012) for any previous associations with ASD and biological plausibility. LocusZoom (Pruim et al., 2010) was used to generate a regional association plot for loci of interest (±400 kb) based on hg18 genome build and HapMap Phase II CEU133 as the linkage disequilibrium (LD) population.
Results
Power Calculations
Total AQ, Social Skills, and Details/Patterns did not reach sufficient power ( >80%) to detect a change by one unit. Communication/Mindreading reached sufficient power at 47% MAF.
Normal Distribution of Scores
Autism Spectrum Quotient scores have been described in previous literature to provide a quantitative measurement of autistic-like traits in the general population (Baron-Cohen et al., 2001). Histograms and Q–Q plots indicate that the scores of the total AQ and each subscale were normally distributed. The histograms and Q–Q plots can be found in Figures S1–S4 in Supplementary Material.
Descriptive Statistics
Characteristics of the current study sample can be found in Table 1. About 51.3% of the sample were female and the mean age at AQ completion was 19.68 (SD = 0.70). Mean total AQ score was 103.2 (SD = 12.60).
Genome-Wide Association Results
After filtering for effect allele frequencies ≤1 and ≥99% and imputation quality R squared values <0.3, 2,462,046 successfully genotyped or imputed SNPs remained for analysis. No SNP associations in any outcome reached genome-wide significance, defined as p-value < 5.0 × 10−8 (Panagiotou and Ioannidis, 2011). Putative associations, defined as p < 1.0 × 10−5, were examined with respect to candidate genes in the SFARI Gene database. None of the genes that were putatively associated with Total AQ, Details/Patterns and Communication/Mindreading were found to be represented in this gene database. One putative gene association found in the Social Skills GWA results was found to be reported in the SFARI Gene database, prompting further investigation.
Social skills genome-wide association results
In the set of 965 Raine Study participants for Social Skills scores, the observed distribution of p-values for the SNPs exhibited minimal deviation from the expected distribution (see Figure 2) indicating minimal test statistic bias or underlying population structure (λ = 0.994). The Manhattan plot –log10 transformation of observed p-values across the genome are displayed in Figure 3.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/54229/fnhum-07-00658-HTML/image_m/fnhum-07-00658-g002.jpg)
Figure 2. Quantile–Quantile Plot of association results for social skills. Q–Q plots compare the observed p-values [−log10(Observed)] to expected p-values. [−log10(Expected)] on the logarithmic scale under the null hypothesis of no significant association.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/54229/fnhum-07-00658-HTML/image_m/fnhum-07-00658-g003.jpg)
Figure 3. Manhattan plot for social skills. Statistical significance of each SNP on the −log10 scale as a function of the chromosome position.
Table 2 displays association results for Social Skills, selected on a p-value below 1.0 × 10−5. None of these associations reached genome-wide significance. However, the rs198198 SNP (p-value = 9.587 × 10−6, effect size = −1.411) is located within the protein kinase C, beta 1 (PRKCB1) gene that has been previously reported to be associated with ASD. Exploratory analysis of data revealed a nominal association of rs198198 with Total AQ score (p-value = 0.002), but no association with Details/Patterns (p-value = 0.434) or Communication/Mindreading (p-value = 0.202). In addition, the rs16946931 SNP (p = 1.78 × 10−6, effect size = 1.681) is located in a region flanking the Cerebellin 1 (CBLN1), a gene previously implicated in ASD. This SNP was nominally associated with the Total AQ score (p-value = 4.27 × 10−4), but no association with Details/Patterns (p-value = 0.472), or Communication/Mindreading (p-value = 0.138). When association analysis was done conditional on rs16946931 for Social Skills, p-values for the 20 additional CBLN1 SNPs became insignificant (p-value > 0.05), demonstrating that high LD between these SNPs was responsible for their inclusion in the model.
Figures 4 and 5 show regional association plots of SNPs rs198198 and rs16946931, respectively.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/54229/fnhum-07-00658-HTML/image_m/fnhum-07-00658-g004.jpg)
Figure 4. Regional association plot of SNP rs198198 (±400 kb). Statistical significance of each SNP on the –log10 scale as a function of the chromosome position. The top SNP is shown as the purple diamond; the correlations (r2) of each of the surrounding SNPs to the top SNP are shown in the indicated colors. Recombination rate is shown in pale blue.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/54229/fnhum-07-00658-HTML/image_m/fnhum-07-00658-g005.jpg)
Figure 5. Regional association plot of SNP rs16946931 (±400 kb). Statistical significance of each SNP on the –log10 scale as a function of the chromosome position. The top SNP is shown as the purple diamond; the correlations (r2) of each of the surrounding SNPs to the top SNP are shown in the indicated colors. Recombination rate is shown in pale blue.
Bioinformatic Evaluation
All SNPs with p-value below 1.0 × 10−5 were evaluated for possible functional potential. As all SNPs were in non-coding regions they were assessed for regulatory potential using available ENCODE data through Regulome DB (Boyle et al., 2012) and the UCSC Genome (2012) browser. No SNP was shown to be part of a transcription factor binding site or be located in a region of the genome likely to have regulatory potential except the nominally associated SNP rs198198 located within an intron of PRKCB1. The major allele (T) was found to be embedded within a near-consensus CCAAT/enhancer binding protein (C/EBP) gamma binding site whereas the minor allele (A) is predicted to ablate C/EBP gamma binding. Further evidence for function comes from ENCODE ChIP-seq data indicating that rs198198 is located in a region subject to histone H3K9 modification.
Discussion
In this paper we report only the second contemporary GWAS for autistic-like traits in the broader community and the first to use a quantitative distribution measured by the AQ. Although we did not find any SNP associations reaching genome-wide significance (defined as p-value < 5.0 × 10−8) this GWAS has provided some evidence to support the further investigation of two previously identified candidate genes, namely PRKCB1 and CBLN1.
PRKCB1 Gene
Protein kinase C enzymes play an important role in signal transduction, regulation of gene expression, and control of cell division and differentiation (Lintas et al., 2009; NCBI, 2012). The alternative splicing of PRKCB1 generates two mRNA isoforms named PRKCB1-1 and PRKCB1-2, coding for two isoenzymes βI and βII. These mRNA isoforms are expressed in a variety of tissues within the central nervous system, including the hippocampus, striatum, suprachiasmatic nucleus, and cerebellar granule cells (Lintas et al., 2009).
Several previous genetic studies have supported a role for PRKCB1 in ASD. A genome-wide linkage scan of 116 families from the Autism Genetic Resource Exchange cohort (Philippi et al., 2005) applied direct physical identity-by-descent mapping using a strict ASD phenotype. Their results demonstrated linkage to a region on chromosome 16p (p-value = 0.0027). Subsequent high-resolution SNP genotyping and analysis of this region showed that haplotypes in the PRKCB1 gene were strongly associated with ASD (Philippi et al., 2005), a finding that was independently confirmed in two other family studies (Philippi et al., 2005) (Lintas et al., 2009) but not in a study on 171 Irish ASD trios (Yang et al., 2007). These haplotypes were not significantly associated with ASD in the current study (results not shown).
These genetic association studies have been supported by a study of PRKCB1 gene expression in post mortem brain tissue from 11 autistic patients and controls (Lintas et al., 2009). The superior temporal gyrus (BA 41/42) region was chosen for study because of a known role in processing socially relevant information and well documented structural and functional abnormalities seen in ASD (Zilbovicius et al., 2006). A significant down regulation of PRKCB1 gene expression was found to be associated with ASD (Lintas et al., 2009) providing additional strong evidence in support of a role for the PRKCB1 gene in the etiology of ASD.
CBLN1 Gene
The CBLN1 gene encodes the cerebellum-specific precursor protein, precerebellin and is involved along with Cerebellin 2 (CBLN2) in encoding ligands for the neurexin-neuroligin trans-synaptic complex (NTSC) (Clarke et al., 2012; Cristino et al., 2013). CBLN1 belongs to the CBLN subfamily of the C1q/tumor necrosis factor superfamily, which plays a role in intercellular signaling, neuronal cell adhesion, brain development, and formation of synapses (Clarke et al., 2012). In a mouse study, Uemura et al. (2010) found that post-synaptic GluRδ2/GRID2 (a member of the δ-type glutamate receptor) interacts with pre-synaptic neurexins through CBLN1 to mediate parallel fiber-Purkinje cell synaptic formation in the cerebellum. A recent study has identified rare copy number variants in NTSC gene family members and proposed their potential biological function with ASD and Tourette syndrome (TS) (Clarke et al., 2012) and mutations in these and other synaptic genes are common in ASD (Arons et al., 2012; Clarke et al., 2012). Therefore, it is plausible that long-range genetic variants affecting CBLN1 expression may be associated with ASD traits.
While no previous genetic linkage or association studies have detected significant associations between CBLN1 and ASD, CBLN2 (belonging to the same CBLN subfamily) on Chromosome 18q has been deleted in TS and is located near two TS translocation breakpoints (Clarke et al., 2012). There is evidence to suggest an overlap between TS and ASD from epidemiological perspectives, with an increased prevalence of TS among autistic patients of 6.5% (Baron-Cohen et al., 1999), compared to between 2 and 3% in the general population (Mason et al., 1998). Common susceptibility genes (Sundaram et al., 2010; Fernandez et al., 2012) have also been reported between ASD and TS.
This study was conducted on a well characterized, population-based cohort of adolescents and ASD measurements were collected without ascertainment for autism risk. Despite this study’s small sample size (n = 965), we were independently able to identify two regions on chromosome 16 that have been previously implicated in the etiology of ASD. While no variants in these regions reached genome-wide significance, they do demonstrate suggestive evidence for association and provide good candidates for future follow-up studies. The difference in observed signal between the two regions is due to high LD around the CBLN1 SNP, rs16946931, which disappears with conditional association analysis. The lack of replication for the previously identified PRKCB1 haplotype is potentially due to different underlying LD structure between populations, as the previous study focused on families (Philippi et al., 2005) whereas our study is composed of unrelated individuals. These observed associations are either intergenic or intronic variants suggesting that they are tagging nearby rare variants that may be causal. In most GWAS, the SNPs showing evidence of association are not usually the actual causal variants at play but rather are in LD with such variants. The failure to replicate the previously reported haplotype associations with ASD and PRKCB1 and any specific SNP associations in both PRKCB1 and CBLN1 does not invalidate the findings in this study. Over the last 5 years, considerable attention has been paid to the potential reasons why true associations may not replicate across independent data sets and genetic heterogeneity, environmental interactions, age-dependent effects, epistasis, and inadequate statistical power have all been cited as possible explanations (Chanock et al., 2007; Shriner et al., 2007; Williams et al., 2007; Greene et al., 2009).
The identification and characterization of the causal variants at play in common complex human traits/diseases is presenting to be a great challenge in human genetics, in particular as a large proportion are not in coding regions directly altering protein structure/function but rather are in non-coding regions (i.e., intronic and intergenic regions) having a likely regulatory function, such as may be the case for the PRKCB1 intronic rs198198 SNP identified in this study. The further characterization of this and other genetic variation within the PRKCB1 and CBLN1 gene regions is beyond the scope of this study possibly requiring exhaustive DNA sequencing to identify all variants within the regions of association followed by a variety of detailed molecular biological analyses, like those we have recently described (Karimi et al., 2009; Kaskow et al., 2013).
Conflict of Interest Statement
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Acknowledgments
Grant sponsor: National Health and Medical Research Council; Grant number: 572613.
Supplementary Material
The Supplementary Material for this article can be found online at http://www.frontiersin.org/Human_Neuroscience/
10.3389/fnhum.2013.00658/abstract
Figure S1. Histogram and Quantile–Quantile Plot of Total AQ Scores. Histograms depict the frequency of the Total AQ observations. Q–Q plots compare the quantiles from the current study (sample quantiles) to the quantiles from a Normal distribution (theoretical quantiles).
Figure S2. Histogram and Quantile–Quantile Plot of Social Skills Scores. Histogram depicts the frequency of the Social Skills observations. Q–Q plots compare the quantiles from the current study (sample quantiles) to the quantiles from a Normal distribution (theoretical quantiles).
Figure S3. Histogram and Quantile–Quantile Plot of Details/Patterns Scores. Histogram depicts the frequency of the Details/Patterns observations. Q–Q plots compare the quantiles from the current study (sample quantiles) to the quantiles from a Normal distribution (theoretical quantiles).
Figure S4. Histogram and Quantile–Quantile Plot of Communication/Mindreading Scores. Histogram depicts the frequency of the Communication/Mindreading observations. Q–Q plots compare the quantiles from the current study (sample quantiles) to the quantiles from a Normal distribution (theoretical quantiles).
References
Alarcón, M., Abrahams, B. S., Stone, J. L., Duvall, J. A., Perederiy, J. V., Bomar, J. M., et al. (2008). Linkage, association, and gene-expression analyses identify CNTNAP2 as an autism-susceptibility gene. Am. J. Hum. Genet. 82, 150–159. doi:10.1016/j.ajhg.2007.09.005
Anney, R., Klei, L., Pinto, D., Regan, R., Conroy, J., Magalhaes, T. R., et al. (2010). A genome-wide scan for common alleles affecting risk for autism. Hum. Mol. Genet. 19, 4072–4082. doi:10.1093/hmg/ddq307
Arons, M. H., Thynne, C. J., Grabrucker, A. M., Li, D., Schoen, M., Cheyne, J. E., et al. (2012). Autism-associated mutations in ProSAP2/Shank3 impair synaptic transmission and neurexin-neuroligin-mediated transsynaptic signaling. J. Neurosci. 32, 14966–14978. doi:10.1523/JNEUROSCI.2215-12.2012
Aulchenko, Y. S., Struchalin, M. V., and van Duijn, C. M. (2010). ProbABEL package for genome-wide association analysis of imputed data. BMC Bioinformatics 11:134. doi:10.1186/1471-2105-11-134
Austin, E. J. (2005). Personality correlates of the broader autism phenotype as assessed by the Autism Spectrum Quotient (AQ). Pers. Individ. Dif. 38, 451–460. doi:10.1016/j.paid.2004.04.022
Baron-Cohen, S., Scahill, V. L., Izaguirre, J., Hornsey, H., and Robertson, M. M. (1999). The prevalence of Gilles de la Tourette syndrome in children and adolescents with autism: a large scale study. Psychol. Med. 29, 1151–1159. doi:10.1017/S003329179900896X
Baron-Cohen, S., Wheelwright, S., Skinner, R., Martin, J., and Clubley, E. (2001). The autism-spectrum quotient (AQ): evidence from Asperger syndrome/high-functioning autism, males and females, scientists and mathematicians. J. Autism Dev. Disord. 31, 5–17. doi:10.1023/A:1005653411471
Boyle, A. P., Hong, E. L., Hariharan, M., Cheng, Y., Schaub, M. A., Kasowski, M., et al. (2012). Annotation of functional variation in personal genomes using RegulomeDB. Genome Res. 22, 1790–1797. doi:10.1101/gr.137323.112
Campbell, D. B., Sutcliffe, J. S., Ebert, P. J., Militerni, R., Bravaccio, C., Trillo, S., et al. (2006). A genetic variant that disrupts MET transcription is associated with autism. Proc. Natl. Acad. Sci. U.S.A. 103, 16834–16839. doi:10.1073/pnas.0605296103
Chakrabarti, B., Dudbridge, F., Kent, L., Wheelwright, S., Hill-Cawthorne, G., Allison, C., et al. (2009). Genes related to sex steroids, neural growth, and social-emotional behavior are associated with autistic traits, empathy, and Asperger syndrome. Autism Res. 2, 157–177. doi:10.1002/aur.80
Chanock, S., Manolio, T., Boehnke, M., Boerwinkle, E., Hunter, D., Thomas, G., et al. (2007). Replicating genotype-phenotype associations. Nature 447, 655–660. doi:10.1038/447655a
Clarke, R. A., Lee, S., and Eapen, V. (2012). Pathogenetic model for Tourette syndrome delineates overlap with related neurodevelopmental disorders including Autism. Transl. Psychiatry 2, e158. doi:10.1038/tp.2012.75
Constantino, J., and Todd, R. (2003). Autistic traits in the general population: a twin study. Arch. Gen. Psychiatry 60, 524–530. doi:10.1001/archpsyc.60.5.524
Cristino, A. S., Williams, S. M., Hawi, Z., An, J. Y., Bellgrove, M. A., Schwartz, C. E., et al. (2013). Neurodevelopmental and neuropsychiatric disorders represent an interconnected molecular system. Mol. Psychiatry doi:10.1038/mp.2013.16
Elsabbagh, M., Divan, G., Koh, Y.-J., Kim, Y. S., Kauchali, S., Marcín, C., et al. (2012). Global prevalence of autism and other pervasive developmental disorders. Autism Res. 5, 160–179. doi:10.1002/aur.239
Fernandez, T. V., Sanders, S. J., Yurkiewicz, I. R., Ercan-Sencicek, A. G., Kim, Y. S., Fishman, D. O., et al. (2012). Rare copy number variants in Tourette syndrome disrupt genes in histaminergic pathways and overlap with autism. Biol. Psychiatry 71, 392–402. doi:10.1016/j.biopsych.2011.09.034
Gauderman, W. J. (2002). Sample size requirements for association studies of gene–gene interactions. Am. J. Epidemiol. 155, 478–484. doi:10.1093/aje/155.5.478
Gillberg, C. L. (1992). The Emanuel Miller Memorial Lecture 1991. Autism and autistic-like conditions: subclasses among disorders of empathy. J. Child Psychol. Psychiatry 33, 813–842. doi:10.1111/j.1469-7610.1992.tb01959.x
Greene, C. S., Penrod, N. M., Williams, S. M., and Moore, J. H. (2009). Failure to replicate a genetic association may provide important clues about genetic architecture. PLoS ONE 4:e5639. doi:10.1371/journal.pone.0005639
Hollander, E., King, A., Delaney, K., Smith, C. J., and Silverman, J. M. (2003). Obsessive-compulsive behaviors in parents of multiplex autism families. Psychiatry Res. 117, 11–16. doi:10.1016/S0165-1781(02)00304-9
Karimi, M., Goldie, L. C., Cruickshank, M. N., Moses, E. K., and Abraham, L. J. (2009). A critical assessment of the factors affecting reporter gene assays for promoter SNP function: a reassessment of −308 TNF polymorphism function using a novel integrated reporter system. Eur. J. Hum. Genet. 17, 1454–1462. doi:10.1038/ejhg.2009.80
Kaskow, B. J., Diepeveen, L. A., Proffitt, J. M., Rea, A. J., Ulgiati, D., Blangero, J., et al. (2013). Molecular prioritisation strategies to identify functional genetic variants in the cardiovascular disease-associated expression QTL Vanin-1. Eur. J. Hum. Genet. doi:10.1038/ejhg.2013.208
Li, X., Zou, H., and Brown, W. T. (2012). Genes associated with autism spectrum disorder. Brain Res. Bull. 88, 543–552. doi:10.1016/j.brainresbull.2012.05.017
Li, Y., and Abecasis, G. (2006). Mach 1.0: rapid haplotype reconstruction and missing genotype inference. Am. J. Hum. Genet. 79, 2290.
Lintas, C., Sacco, R., Garbett, K., Mirnics, K., Militerni, R., Bravaccio, C., et al. (2009). Involvement of the PRKCB1 gene in autistic disorder: significant genetic association and reduced neocortical gene expression. Mol. Psychiatry 14, 705–718. doi:10.1038/mp.2008.21
Lundström, S., Chang, Z., Råstam, M., Gillberg, C., Larsson, H., Anckarsäter, H., et al. (2012). Autism spectrum disorders and autistic like traits: similar etiology in the extreme end and the normal variation. Arch. Gen. Psychiatry 69, 46–52. doi:10.1001/archgenpsychiatry.2011.144
Mason, A., Banerjee, S., Eapen, V., Zeitlin, H., and Robertson, M. M. (1998). The prevalence of Tourette syndrome in a mainstream school population. Dev. Med. Child Neurol. 40, 292–296.
National Human Genome Research Institute. (2012). A Catalog of Genome-Wide Association Studies [Online]. Bethesda, MD. Available at: http://www.genome.gov/26525384 (accessed 2012).
NCBI. (2012). PRKCB Protein Kinase C, Beta [Online]. Bethesda: National Center for Biotechnology Information.
Newnham, J. P., Evans, S. F., Michael, C. A., Stanley, F. J., and Landau, L. I. (1993). Effects of frequent ultrasound during pregnancy: a randomised controlled trial. Lancet 342, 887–891. doi:10.1016/0140-6736(93)91944-H
Panagiotou, O. A., and Ioannidis, J. P. A. (2011). What should the genome-wide significance threshold be? Empirical replication of borderline genetic associations. Int. J. Epidemiol. 41, 273–286. doi:10.1093/ije/dyr178
Philippi, A., Roschmann, E., Tores, G., Lindenbaum, P., Benajou, A., Germain-Leclerc, L., et al. (2005). Haplotypes in the gene encoding protein kinase c-beta (PRKCB1) on chromosome 16 are associated with autism. Mol. Psychiatry 10, 950–960. doi:10.1038/sj.mp.4001704
Posserud, M.-B., Lundervold, A. J., and Gillberg, C. (2006). Autistic features in a total population of 7–9-year-old children assessed by the ASSQ (Autism Spectrum Screening Questionnaire). J. Child Psychol. Psychiatry 47, 167–175. doi:10.1111/j.1469-7610.2005.01462.x
Pruim, R. J., Welch, R. P., Sanna, S., Teslovich, T. M., Chines, P. S., Gliedt, T. P., et al. (2010). LocusZoom: regional visualization of genome-wide association scan results. Bioinformatics 26, 2336–2337. doi:10.1093/bioinformatics/btq419
R Development Core Team. (2008). R: A Language and Environment for Statistical Computing. Vienna: R Foundation for Statistical Computing.
Ronald, A., Butcher, L. M., Docherty, S., Davis, O. S. P., Schalkwyk, L. C., Craig, I. W., et al. (2010). A genome-wide association study of social and non-social autistic-like traits in the general population using pooled DNA, 500 K SNP microarrays and both community and diagnosed autism replication samples. Behav. Genet. 40, 31–45. doi:10.1007/s10519-009-9308-6
Ronald, A., Happé, F., and Plomin, R. (2005). The genetic relationship between individual differences in social and nonsocial behaviours characteristic of autism. Dev. Sci. 8, 444–458. doi:10.1111/j.1467-7687.2005.00433.x
Ronald, A., and Hoekstra, R. A. (2011). Autism spectrum disorders and autistic traits: a decade of new twin studies. Am. J. Med. Genet. B Neuropsychiatr. Genet. 156, 255–274. doi:10.1002/ajmg.b.31159
Russell-Smith, S. N., Maybery, M. T., and Bayliss, D. M. (2011). Relationships between autistic-like and schizotypy traits: an analysis using the Autism Spectrum Quotient and Oxford-Liverpool Inventory of Feelings and Experiences. Pers. Individ. Dif. 51, 128–132. doi:10.1016/j.paid.2011.03.027
Shriner, D., Vaughan, L. K., Padilla, M. A., and Tiwari, H. K. (2007). Problems with genome-wide association studies. Science 316, 1840–1842. doi:10.1126/science.316.5833.1840c
Simons Foundation Autism Research Initiative. (2012). SFARI Gene: A Modular Database for Autism Research [Online]. New York. Available at: https://gene.sfari.org/autdb/Welcome.do (accessed June 4, 2012).
Skuse, D. H., Mandy, W. P. L., and Scourfield, J. (2005). Measuring autistic traits: heritability, reliability and validity of the Social and Communication Disorders Checklist. Br. J. Psychiatry 187, 568–572. doi:10.1192/bjp.187.6.568
SNAP. (2008). SNP Annotation and Proxy Search [Online]. Cambridge. Available at: http://www.broadinstitute.org/mpg/snap/ldsearch.php (accessed 2012).
St Pourcain, B., Wang, K., Glessner, J., Golding, J., Steer, C., Ring, S., et al. (2010). Association between a high-risk autism locus on 5p14 and social communication spectrum phenotypes in the general population. Am. J. Psychiatry 167, 1364–1372. doi:10.1176/appi.ajp.2010.09121789
Stewart, M. E., and Austin, E. J. (2009). The structure of the Autism-Spectrum Quotient (AQ): evidence from a student sample in Scotland. Pers. Individ. Dif. 47, 224–228. doi:10.1016/j.paid.2009.03.004
Sundaram, S. K., Huq, A. M., Wilson, B. J., and Chugani, H. T. (2010). Tourette syndrome is associated with recurrent exonic copy number variants. Neurology 74, 1583–1590. doi:10.1212/WNL.0b013e3181e0f147
UCSC Genome. (2012). Bioinformatics, [Online]. Santa Cruz. Available at: http://genome.ucsc.edu/ [accessed 2012].
Uemura, T., Lee, S. J., Yasumura, M., Takeuchi, T., Yoshida, T., Ra, M., et al. (2010). Trans-synaptic interaction of GluRdelta2 and Neurexin through Cbln1 mediates synapse formation in the cerebellum. Cell 141, 1068–1079. doi:10.1016/j.cell.2010.04.035
Veenstra-VanderWeele, J., Christian, S. L., and Cook, E. H. (2004). Autism as a paradigmatic complex genetic disorder. Annu. Rev. Genomics Hum. Genet. 5, 379–405. doi:10.1146/annurev.genom.5.061903.180050
Wang, K., Zhang, H., Ma, D., Bucan, M., Glessner, J. T., Abrahams, B. S., et al. (2009). Common genetic variants on 5p14.1 associate with autism spectrum disorders. Nature 459, 528–533. doi:10.1038/nature07999
Whitehouse, A. J. O., Hickey, M., and Ronald, A. (2011). Are autistic traits in the general population stable across development? PLoS ONE 6:e23029. doi:10.1371/journal.pone.0023029
Williams, S. M., Canter, J. A., Crawford, D. C., Moore, J. H., Ritchie, M. D., and Haines, J. L. (2007). Problems with genome-wide association studies. Science 316, 1840–1842. doi:10.1126/science.316.5833.1840c
Yang, M. S., Cochrane, L., Conroy, J., Hawi, Z., Fitzgerald, M., Gallagher, L., et al. (2007). Protein kinase C-beta 1 gene variants are not associated with autism in the Irish population. Psychiatr. Genet. 17, 39–41. doi:10.1097/YPG.0b013e3280115428
Keywords: autistic-like traits, genome-wide association, PRKCB1, autism spectrum disorder, autism spectrum quotient, CBLN1
Citation: Jones RM, Cadby G, Melton PE, Abraham LJ, Whitehouse AJ and Moses EK (2013) Genome-wide association study of autistic-like traits in a general population study of young adults. Front. Hum. Neurosci. 7:658. doi: 10.3389/fnhum.2013.00658
Received: 25 April 2013; Accepted: 22 September 2013;
Published online: 11 October 2013.
Edited by:
Charles Claudianos, University of Queensland, AustraliaReviewed by:
Raymond Allan Clarke, Ingham Institute University of Western Sydney, AustraliaIrina Voineagu, The University of New South Wales, Australia
Copyright: © 2013 Jones, Cadby, Melton, Abraham, Whitehouse and Moses. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) or licensor are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Eric K. Moses, Centre for Genetic Origins of Health and Disease, University of Western Australia, 35 Stirling Highway, Crawley, Perth, WA 6009, Australia e-mail: eric.moses@uwa.edu.au
†Rachel Maree Jones and Gemma Cadby have contributed equally to this work.