Novel Autism Subtype-Dependent Genetic Variants Are Revealed by Quantitative Trait and Subphenotype Association Analyses of Published GWAS Data

The heterogeneity of symptoms associated with autism spectrum disorders (ASDs) has presented a significant challenge to genetic analyses. Even when associations with genetic variants have been identified, it has been difficult to associate them with a specific trait or characteristic of autism. Here, we report that quantitative trait analyses of ASD symptoms combined with case-control association analyses using distinct ASD subphenotypes identified on the basis of symptomatic profiles result in the identification of highly significant associations with 18 novel single nucleotide polymorphisms (SNPs). The symptom categories included deficits in language usage, non-verbal communication, social development, and play skills, as well as insistence on sameness or ritualistic behaviors. Ten of the trait-associated SNPs, or quantitative trait loci (QTL), were associated with more than one subtype, providing partial replication of the identified QTL. Notably, none of the novel SNPs is located within an exonic region, suggesting that these hereditary components of ASDs are more likely related to gene regulatory processes (or gene expression) than to structural or functional changes in gene products. Seven of the QTL reside within intergenic chromosomal regions associated with rare copy number variants that have been previously reported in autistic samples. Pathway analyses of the genes associated with the QTL identified in this study implicate neurological functions and disorders associated with autism pathophysiology. This study underscores the advantage of incorporating both quantitative traits as well as subphenotypes into large-scale genome-wide analyses of complex disorders.


Introduction
Autism spectrum disorders (ASDs) represent a group of neurodevelopmental disorders that are characterized by impaired reciprocal social interactions, delayed or aberrant communication, and stereotyped, repetitive behaviors, often with restricted interests [1,2]. With a concordance rate as high as 90% based on twin studies [3], ASDs are among the most heritable of neuropsychiatric conditions. Although specific genes have been found to be causal for several syndromes that are sometimes associated with autistic symptoms (for example, Fragile X [4], Retts [5,6], and tuberous sclerosis [7,8], there are no unequivocal genetic markers for non-syndromic idiopathic autism. Thus, a considerable amount of effort has been devoted to identifying genetic mutations or variants that associate with these perplexing and often devastating, life-long disorders. Because of the relatively high prevalence of ASDs in the general population (,1:110), genomewide association (GWA) analyses have been used recently to search for common variants that may associate with increased susceptibility to this set of disorders [9][10][11][12]. However, despite case-control studies that have now exceeded many thousands of subjects and more than 500,000 SNPs, only a few significantly associated SNPs have been identified. In addition, replication of these SNPs in independent studies has not been successful. The inability to replicate findings from GWA analyses may be in part due to the genetic heterogeneity of the ASD population, thus giving rise to increased ''noise'' in the data. This genetic heterogeneity is likely responsible for the well-noted phenotypic and symptomatic heterogeneity among individuals with autism.
In a recent paper, we demonstrated that the ASD population can be divided into at least 4 phenotypic subgroups on the basis of cluster analyses of 123 severity scores taken from each individual's diagnostic assessment using the Autism Diagnostic Interview-Revised [13]. The resulting subgroups included one with severe language impairment, another with mild severity across all items, a third of intermediate severity, and a fourth of moderate severity with a higher frequency of savant skills. We further demonstrated by gene expression profiling of lymphoblastoid cell lines from 3 of these subgroups (excluding the intermediate) and nonautistic controls that cells from each of these subgroups exhibited differentially expressed genes relative to those of the controls, but also were distinguishable from each other in terms of unique, subtype-specific differentially expressed genes [14]. These studies thus support the concept that different subgroups of autistic individuals may exhibit subtype-dependent biological differences due to genetic variation. We therefore hypothesized that genetic association analyses of such ASD subtypes with SNPs that are identified and filtered according to their association with quantitative traits relevant to ASDs should reveal more significant SNPs with increased statistical power. Here, we report the identification of 18 novel and highly significant SNPs that are associated with at least one of 4 different subtypes of ASDs, based on the combination of quantitative trait association analyses and subtype-dependent genetic association analyses using trait-associated SNPs, with 10 of these SNPs replicated in a second ASD subtype representing an independent case cohort.

Results
The overarching goal of these studies was to identify single nucleotide polymorphisms (SNPs) that are associated with both autistic traits and clinical subtypes of autism that are manifested by separate case cohorts. To accomplish this, we combined quantitative trait analysis and subphenotype association analyses using the wealth of genome-wide association (GWA) data published by Wang et al. in 2009 [9].

Quantitative Trait Association Analyses
The flowchart in Fig. 1 describes the experimental design and analyses that were used to derive the combined set of 18 novel and statistically significant SNPs that associate with 4 separate subtypes of ASDs. Raw item scores from the ADI-R score sheets of 2939 ASD cases were summed for spoken language skills, non-verbal communication, play skills, social development, and insistence on sameness/rituals, as described previously [13]. The specific items used to obtain the total score per ''trait'' category for each individual are shown in Table S1 and the profiles of total scores for each category are shown for the 2939 individuals in Fig. S1.
We then conducted quantitative trait association analyses using the distribution of scores in each of the categories to identify sets of SNPs that associate with symptomatic severity of each of the behaviors listed in Table S1. A nominal p-value#10 25 was used in this ''discovery phase'' designed to screen candidate SNPs that are potentially of functional significance with respect to the ASD traits. These sets of symptom-associated SNPs (or quantitative trait loci, QTL) are shown in Table S2. ASD subtype-dependent genetic association analyses with trait-associated SNPs Next, we performed cluster analyses as described by Hu and Steinberg [13] to divide the autistic cases into 4 non-overlapping phenotypic subgroups according to symptomatic severity profiles derived from 123 items on the ADI-R assessments. We postulated that this subtyping procedure, which reduces the behavioral/ symptomatic heterogeneity among the cases within each subgroup, is likely to also restrict the genetic heterogeneity within each subgroup. The results of K-means cluster analyses (K = 4) of the ADIR data from the 2939 individuals (which included 1867 genotyped cases from the GWA study) is shown in Fig. S2A. An unsupervised principal components analysis (PCA) was used to demonstrate that the 4 subgroups of individuals with ASDs from the K-means cluster analyses results in effective separation of cases according to similarity of symptomatic profile as defined by the 123 item scores from the ADI-R (Fig. S2B). The resulting phenotypic subgroups were then used in genetic association analyses with each set of QTL derived from the five quantitative trait association analyses, where the cases were either divided into Figure 1. Diagram of study design illustrating sequential application of quantitative trait and subphenotype association analyses. The quantitative trait association analyses were performed with 1867 cases for which the ADI-R assessments were available, using the complete set of genetic polymorphisms to identify SNPs that may have functional relevance to traits associated with ASDs. This had the net effect of filtering the set of SNPs from 513,312 to 167 QTL. In the second phase, each set of trait-associated SNPs (QTL) were employed in case-control association analyses with 2438 controls and the combined cases (1867) as well as with 4 ASD subgroups, each representing an independent cohort of cases (363-639 cases per subtype). From these analyses, only 18 QTL with Bonferroni-adjusted p-values,0.05 in at least one subtype for each trait were combined for the final set of genetic association analyses using combined cases as well as ASD subtypes against the controls. doi:10.1371/journal.pone.0019067.g001 the 4 ASD subtypes or used as a combined case group against the 2438 nonautistic controls. These analyses produced 5 sets of QTL that were associated with specific subtypes of ASDs (Tables 1, 2 , 3, 4, 5). It is interesting to note that SNPs associated with certain genes (such as HTR4, BACH2, and TINAG among Language QTL, and OR5B3, OR5B17, and ERBB4 among Nonverbal QTL) are replicated in more than one ASD subtype. Finally, significant SNPs with Bonferroni-adjusted p-values#0.05 from each of the 4 separate subtype-dependent association analyses with QTL were combined into a single set containing 18 unique SNPs, and case-control association analyses were performed for which the cases were either combined into one group or divided according to subphenotype. Table 6 shows the SNPs associated with each subtype of ASDs that resulted from the final association analyses using the combined highly significant QTL and 4 non-overlapping subgroups of ASD cases. Eighteen of the SNPs have pvalues,0.05 in at least one subtype even after using the stringent Bonferroni correction for multiple comparisons. Note that 10 of the SNPs, including rs317985, rs7785107, rs11671930, rs7950390, rs12266938, rs3861787, rs7725785, rs1827924, rs1231339, and rs757099, are associated with more than one subtype. Two of the replicated SNPs (rs317985 and rs7785107) are significant in two subtypes after Bonferroni adjustment (p,0.05), while the remaining 8 exhibit lower levels of significance (nominal p-values from 0.0037-0.051 or FDR_BH adjusted p-values of 0.0087-0.088) in the second (or third) subtype. Association of these QTL with more than one subtype of ASD serves as a replication for these 10 SNPs. While the replication of SNPs in different subtypes may be interpreted as evidence for those SNPs being highly penetrant for ASDs, the subtype-dependent differences in minor allele frequency (MAF) and odds ratios (OR) associated with the shared SNPs suggest that the subtypes are genetically heterogeneous. Figure 2 summarizes the extent of SNP overlap among the 4 ASD subtypes (case cohorts) and clearly demonstrates that the odds ratios are distinctly different for different subtypes that share the same SNP. All of the QTL associated with specific genes are present in non-exonic (promoter, intronic, or intergenic) regions. Interestingly, all but one of the SNPs residing within intergenic regions can be associated by band position to rare copy number variants (CNV) that have been recently identified for ASDs [15], thus providing further support for the probable relevance of these SNPs to ASDs. These CNVs are noted in Table 6. Effect of subtyping of ASDs on association of previously identified SNPs within Chr5p14.1

Partial replication of SNPs between subtypes of ASDs
Because none of the SNPs identified in this re-analysis overlapped with those of the previously published genome-wide association study [9], we examined the association of the 6 SNPs that were reported in the published study with our subphenotypes. Table 7 shows that only the ''Moderate'' ASD subtype (363 cases) is associated with two of the SNPs, with Bonferroni-adjusted pvalues of 0.035 and 0.053. Interestingly, these 2 SNPs have the least significant combined p-values in the published study. The remaining 4 SNPs were suggestively significant with FDR_BH- adjusted p-values of 0.074 in this subtype. The combined cases (1867 individuals in all) as well as the other 3 ASD subtypes show no association with any of the 6 SNPs even though there are more cases in each of these groups than in the Moderate group. This finding further illustrates the value of analyzing subphenotypes of ASDs in genome-wide association analyses in order to reduce the genetic heterogeneity of the population studied.

Pathway analyses of SNP-containing genes
To obtain a better understanding of how the novel SNPs identified in this study potentially relate to the biology of autism, pathway analysis was conducted to develop a better sense of the relationships among the SNP-associated genes and their impact on higher level functions and diseases. Fig. 3 shows a gene network constructed using Pathway Studio 7 which includes seven of the 9 genes associated with SNPs found within gene promoters or introns. Of the 7 genes, HTR4 and GCH1 are ''hubs'' showing the highest ''connectivity'' with other components within the network. The specific relationships between these two genes and other network components are more clearly illustrated in Figures  S3 and S4. It is noted that many of the cellular and higher level processes in this network, such as neurogenesis, axonogenesis, steroid metabolism, cell proliferation, long-term synaptic potentiation, learning and memory are relevant to identified deficits in ASDs.

Discussion
We have previously shown that the autistic population can be divided into subgroups of affected individuals according to symptomatic profile through cluster analyses of severity scores from the ADI-R assessment for each individual with an ASD [13]. Three of the 4 resulting subgroups were shown to exhibit distinct, though partially overlapping, gene expression profiles, each relative to a group of nonautistic controls, implying that both unique and shared genes are associated with the respective phenotypes [14]. We therefore applied the same rationale and methods in subtyping individuals with ASDs for this re-analysis of previously published genome-wide association data [9]. We first employed quantitative trait association analyses to the .500,000 SNPs tested in order to prioritize SNPs that might correlate with a behavioral or symptomatic ''trait'' relevant to ASDs. These quantitative traits for each individual with an ASD were derived from the sums of severity scores from the ADI-R items that described severity of language impairment, deficits in nonverbal communication, impaired play skills, insistence on sameness/ rituals, and delayed social development. The specific items used to establish each quantitative trait (listed in Table S1) were shown by Hu and Steinberg to exhibit differential severity among the several subtypes of ASDs [13]. In this ''discovery'' stage of the experiment, quantitative trait association analyses across the 5 selected traits produced a filtered set of 167 unique SNPs from the original 513,312 SNPs (Table S2). Subsequent association analyses of these QTL with both combined and subtyped individuals with ASDs revealed 18 novel SNPs which were found to be highly significant (Bonferroni-adjusted p-value,0.05) in at least one subtype of autism (Table 6). Interestingly, many of the language QTL from Table S2 are strongly associated with the severely language-impaired ASD subtype. Of the 18 novel SNPs, 10 SNPs were replicated, albeit at a more modest significance level, in at least one of the other subtypes. Two are replicated with Bonferroni-significant p-values while the other 8 are suggestively significant (FDR-BH adjusted p-value,0.09). While it may be argued that the standard statistical methods used here to correct for multiple testing are not stringent enough, there is still no clear consensus on the significance level that is appropriate for a candidate gene (or SNP) study such as this, where the candidate SNPs are determined by prior quantitative trait association analyses. More significantly, different minor allele frequencies and odds ratios are observed for the SNPs that are associated with more than one subtype (Table 6 and Fig. 2), thus reflecting the genetic heterogeneity among the ASD phenotypes, which is teased apart by the subtyping procedures employed here. It is noteworthy that no significant SNPs were identified when all 1867 individuals with ASDs were analyzed against 2438 nonautistic controls, thus underscoring the importance of phenotypic subtyping to unearthing SNP associations with ASDs.
By comparison, the original genome-wide study which was based on the combined analysis of 2503 cases and more than 7000 controls across 2 independent datasets in the ''discovery'' phase, identified 1 SNP that reached genome-wide significance and 5 additional SNPs of nominal significance in one intergenic region on chr5p14.1 that was located between 2 cadherin genes [9]. In this study, we replicated the association of 2 of the SNPs but only in the ''Moderate'' subtype of ASDs which consisted of only 363 cases (Table 7). Neither the other 3 subtypes (with 387-639 cases per subgroup) nor the combined cases group (consisting of 1867 cases) showed any association with these 6 SNPs. Furthermore, while none of the 6 SNPs identified in the published study correlated with expression level of either cadherin (CDH9 or CDH10) in the cortical brain of 93 genotyped human subjects, we note that 2 of the SNP-associated genes identified in the current study, HTR4 and CCL25, were found to be differentially expressed (FDR,5%) in lymphoblastoid cell lines from individuals with ASDs in a previous study by our laboratory [14]. This overlap of differentially expressed genes with those associated with at least some of the novel SNPs lends support to the potential functional relevance of these SNPs. Aside from the SNPs located within the promoter or intron regions of genes, seven of the other eight significant SNPs (Table 6) are located in intergenic regions that are linked by band position to rare copy number variants (CNVs) that have been recently associated with autism [15]. The finding that all 18 of the novel SNPs are in non-exonic (promoter and noncoding) regions further suggests that gene regulatory processes likely leading to differential gene expression, rather than to alterations of gene products, may contribute to the hereditary component of ASDs. The presence of the identified QTL in ASDrelated CNVs provides additional support for the potential relevance of these novel SNPs to ASDs. The current study thus illustrates the advantage of combining both quantitative traits and defined ASD phenotypes in analyzing genome-wide genetic data from individuals with this complex disorder.

Quantitative trait and phenotype incorporation into genetic analyses
There have been relatively few large-scale genetic studies on ASDs that incorporate quantitative traits or phenotypes into genome-wide analyses. The value of using phenotypic subsets of ASDs in genetic analyses was demonstrated in an early linkage study by Nurmi et al. who found that subsetting the subjects (and families) according to presence or absence of savant skills increased the HLOD score of the D15S511 marker within the GABRB3 gene on chr15q11-q13 from 0.8 in the entire study population to 2.6 for the savant-skills-positive subset of individuals [16]. Later, Alarcon et al. used quantitative trait analysis of 2 language scores on the ADI-R combined with ordered subsets analysis for subject stratification to reveal a strong language QTL on 7q35, and observed additional linkage peaks on several other chromosomes without stratification [17]. Chen et al. used a combination of 7 ADI-R scores describing nonverbal communication and the same approach to uncover QTL for nonverbal communication [18]. Similarly, Duvall et al. used quantitative trait analysis based upon scores from the Social Responsiveness Scale to discover QTL for social endophenotypes with strong linkage to chromosomes 11 and 17 across all samples and additional linkages to chromosomes 4, 8, and 10 in ''male-only'' families [19]. More recently, Nijmeijer et al. identified loci overlapping ASDs and attention-deficit/ hyperactivity disorder (ADHD) in a genome-wide study of children with ADHD and autistic traits (but not classical or atypical autism) by using total and subscale scores from the Social Communication Questionnaire (an autism diagnostic instrument) as well as by including scores related to severity of ADHD as covariates [20]. Along a similar vein, Ronald et al. examined both social and nonsocial autistic-like traits in the general population, quantified on the basis of DSM-IV social and non-social scales, to tease out 22 and 20 SNPs that associated with these traits, respectively [21]. However, only one SNP survived testing in an independent sample and its p-value was still only nominally significant in an all male sample. Recently, quantitative trait analysis was applied in a genome-wide linkage study of two repetitive ASD behavior phenotypes, ''insistence on sameness'' (IS) and ''repetitive and stereotyped behavior'' (RSMA) [22]. In this study, Cannon et al. demonstrated IS and RSMA were linked to 2 separate loci, each exhibiting respective HLOD scores of 3.42 and 3.93, suggesting that more narrowly-defined phenotypes improve linkage signals. When social behavior, as reflected by severity scores on the Social Responsiveness Scale was used as a quantitative trait, Coon et al. not only replicated a susceptibility locus identified previously in the Duvall et al. study [19], but also identified additional linkage peaks with HLOD scores ranging from 2.23-3.64 [23]. Interestingly, St. Pourcain et al. recently reported that one of the SNPs identified at Chr5p14 by GWA analyses [9] also associates by quantitative trait analyses with the social communication phenotype [24]. They suggested that this SNP, rs4307059, may also associate with multiple subthreshold impairments across a variety of ASD-related functions.
In the present study, we observe that rs4307059 is suggestively associated only with the Moderate ASD subtype, with an unadjusted p-value of 0.074 (Table 7). Thus, the above-mentioned studies, in combination with ours, reinforce the concept that the incorporation of quantitative traits into genetic analyses greatly improves the ability to elicit significant genetic information on complex disorders such as ASDs. Furthermore, the genes associated with the QTL in Tables 1, 2, 3,4, 5 and Table S2 may provide clues to the biology underlying the various autistic traits. Our study, which uniquely combines QTL with ASD subtypedependent genetic association analyses, provides an additional refinement to genome-wide studies that results not only in the identification of a greater number of significant SNPs, but also concomitant and internal validation of SNPs that are replicated between and among the ASD subtypes which represent independent cohorts of affected individuals. Since the number of individuals in each subgroup ranges from 363-639, our study also demonstrates the increased statistical power inherent in this approach. Clearly, there are many ways that individuals with ASDs can be divided into subgroups (such as by regression following normal development, co-morbid symptoms, IQ, etc.), and some of these methods may prove to be superior to the method described here and in our previous study [13]. However, the importance of using phenotypically distinct cohorts in large-scale genetic association analyses is illustrated here using a method that is capable of distinguishing meaningful phenotypes across the whole spectrum of ASDs (severe to mild), based on severity profiles across a broad range of behaviors and symptoms characteristic of ASDs. ChiSquare value; f Odds ratio with respect to major allele; g p-value adjusted using Benjamini-Hochberg False Discovery Rate method; h Bonferroni-adjusted p-value*CNV located in Chromosomal band as reported by Pinto et al. [15]. doi:10.1371/journal.pone.0019067.t007

Biological relevance of SNP-containing genes
To examine the biological processes and pathways that might be impacted by the SNP-associated genes, we performed pathway analyses on the genes using Pathway Studio 7 software. Fig. 3 shows that 7 of the 9 SNP-containing genes are included in a gene network with many genes leading to ''axonogenesis,'' a function which has been reported to be enriched among genes overlapping copy number variants in autistic patients [25]. HTR4 and GCH1 are network ''hubs'' connecting with many other genes, cellular processes and disorders (see Figs. S3 and S4 for specific connections). As illustrated in Fig. S3, HTR4 [5-hydroxytryptamine (serotonin) receptor 4] regulates neurogenesis, long-term synaptic potentiation and, in turn, learning and memory, as well as the release of neurotransmitters (dopamine, acetylcholine), peptide hormones (AVP, OXT, PRL, VIP) and steroid compounds (cortisol, corticosterone). Thus, any alteration in the expression or function of this gene can be expected to have wide-ranging consequences on processes known to be affected by ASDs. It is notable that one of the SNPs associated with HTR4, rs7725785, is associated with three ASD subtypes (Table 6). However, the odds ratio is 1.44 for the severe language-impaired subtype while it is 0.68 and 0.74 for the moderate and mild subtypes, respectively. Interestingly, a reduction in HTR4 expression was observed only for the language-impaired subtype of ASDs and not for the moderate (savant) or mild subtypes in our earlier gene expression study [14], suggesting that this SNP, which is located at an intron boundary within HTR4, may play a role in level of expression. Genetic variants in HTR4 have also been associated with schizophrenia [26], bipolar disorder [27] and attention deficit/ hyperactivity disorder [28], further suggesting some functional overlap among these neurological disorders. More recently, a de novo translocation on chromosome 5 close to HTR4 has been identified in an autistic boy [29], lending additional support for the relevance of this gene in ASDs.
The other hub gene, GCH1 [GTP cyclohydrolase I], is the ratelimiting enzyme in the de novo biosynthesis of tetrahydrobiopterin which is, in turn, required for the biosynthesis of folate, serotonin, dopamine, and catecholamines (Fig. S4), all of which are important for neural development and functions. It is interesting to note that elevated expression of GCH1 has been implicated in mood disorders [30], while genetic polymorphisms or mutations in GCH1 have been associated with pain sensitivity [31][32][33], and dystonia [34], which are often associated with ASDs. Although these genes are not likely to be causal for ASDs, genetic polymorphisms in them may be associated with some of the comorbid symptoms or pathobiology of ASDs.  Table 6. Genes from Table 6  Together, these analyses provide support for the possible biological relevance of these QTL to ASDs and identify additional candidate genes, such as HTR4 and GCH1, for functional testing. Importantly, this study also reveals genetic biomarkers which not only may be used for diagnostic screening of ASDs, but also offer the additional advantage of being associated with ASD subtypes that may be linked to specific and targeted therapies through pharmacogenomics studies. Finally, the association of different SNPs with the 4 subtypes of ASDs described here reinforces the idea that there are multiple genetic etiologies giving rise to the autistic spectrum, while the shared SNPs between different subtypes may reveal common genetic mechanisms responsible for core symptoms.

Summary
This study is the first to demonstrate the value of using a combination of quantitative trait analysis and subphenotyping of individuals with ASDs to identify genetic variants (SNPs) that associate with specific behavioral phenotypes of ASDs. It is noted that no Bonferroni-significant SNPs are detected when all 1867 autistic cases are combined into a single group and compared against 2438 non-autistic controls. Thus, even though the number of cases is lower in each of the subgroups, there is more power to detect statistically significant SNPs associated with the more homogeneous subgroups of individuals with ASDs than with the combined group of cases. Subtyping also creates independent case cohorts which are shown to replicate 10 of the novel SNPs identified in this study, thus providing a form of internal validation. Differences in minor allele frequencies of these 10 SNPs in the different cohorts further demonstrate the genetic heterogeneity among the subtypes. Together, these findings not only reveal novel subtype-dependent candidate genes for ASDs, but also identify potential genetic markers for diagnostic screening.

GWA data used for this study
Genome-wide association data from the study by Wang et al. [9] were downloaded from the Autism Speaks website at ftp://ftp. autismspeaks.org/Data/CHOP_PLINK/AGRERELEASE.ped. For this study, we used the file named CHOP.clean100121 where the data were '

ADI-R data
ADI-R assessments for 2939 individuals were obtained from Autism Speaks through Dr. Vlad Kustanovich of the Autism Genetics Resource Exchange. Of these, 1867 individuals were among the cases genotyped by Wang et al. [9]. Scores of 123 items on the ADI-R score sheets of each individual were analyzed as described by Hu and Steinberg [13] to identify ASD subtypes (that is, phenotypic subgroups) which are represented in Fig. S2.

Determination of quantitative traits
Raw item scores from the Autism Diagnostic Interview-Revised (ADI-R) score sheets of 2939 ASD cases were summed for spoken language skills, non-verbal communication, play skills, social development, and insistence on sameness/ritualistic behaviors. The specific items used to obtain the total score per category for each individual, shown in Table S1, were noted in an earlier study [13] to exhibit average differences in severity among several subtypes of ASDs, described below. The sums of item scores within each of the 5 categories were used as quantitative traits for genetic association analyses using the genotype data reported by Wang et al. [9]. Profiles of the traits across the 2939 individuals are shown in Fig. S1.

Subtyping of individuals with ASDs
Two thousand nine hundred and thirty nine (2939) individuals with ASDs were divided into phenotypic subgroups using clustering tools within the MeV software package [35], as previously described by Hu and Steinberg [13]. Briefly, subtyping of individuals with ASDs involved K-means cluster analysis (with K = 4) of scores from 123 items from the ADI-R score sheets of each individual which were adjusted as described to fall into a severity range of 0 (normal) to 3 (highest severity). A figure of merit analysis (not shown) indicated that the individuals with ASDs were optimally represented by 4 subgroups. A principal components analysis (PCA) was then used to verify that the 4 subgroups of individuals identified by the K-means cluster analyses were distinguishable by this unsupervised test. Fig. S2 shows the symptomatic profiles of the 4 ASD subtypes as well as their separation into discernible clusters by PCA. The subtypes are named ''Language-impaired,'' ''Intermediate,'' ''Moderate,'' and ''Mild,'' and contain 639, 478, 363, and 387 cases, respectively.

Quantitative trait association analyses
Using the score totals for the 5 categories of autistic symptoms exhibited by each of 1867 cases as quantitative traits, we utilized PLINK [36], a program to analyze whole genome association data, to perform quantitative trait association analyses with the genotype data reported by Wang et al. [9]. From the results of each of the 5 analyses (Table S2) in this discovery phase, we selected SNPs with unadjusted p-values,10 25 , which prioritized SNPs filtered by association with quantitative traits relevant to ASDs. These filtered sets of SNPs were then used in case-control association analyses, as described below.

Case-control association tests
Using the 5 sets of QTL generated for each symptom category, association analyses using PLINK were performed with each set where the cases, in addition to being treated as a single group of 1867 cases vs. 2438 controls, were also divided into the 4 ASD subtypes that were determined by the ADI-R cluster analyses described above. It should be noted that each ASD subgroup represents an entirely separate cohort of cases. From each of the 5 sets of genetic association analyses with subtypes and QTL (Tables 1, 2, 3, 4, 5), we selected SNPs associated with each ASD subtype with a Bonferroni-adjusted p-value#0.05 and combined them (a total of 18 unique SNPs) for a second case-control association analysis using the combined and subtyped ASD cases against controls. The results of this final analysis are shown in Table 6.

Pathway analysis
Pathway Studio 7 software (Ariadne Genomics, Inc.) was used to generate relational gene networks using the SNP-containing genes listed in Table 6. Figure S1 Quantitative trait profiles generated by summing the severity scores for ADI-R items for each trait listed in Table S1. (TIF) Figure S2 Identification of ASD subtypes by cluster analyses. A) Symptomatic profiles of the 4 ASD subtypes that resulted from K-means cluster analyses of 123 ADI-R severity scores per individual. In this figure, each row represents an individual and each column represents an item on the ADI-R. Black represents a score of 0 which is considered ''normal,'' while the intensity of red indicates severity scores ranging from 1-3. Gray represents unavailable data. The wide band of intensely red items in the language-impaired subgroup corresponds to spoken language. The 12 columns at the extreme right in each block represent items corresponding to ''Savant skills,'' which appear to be present at a slightly higher frequency in the group labeled ''Moderate.'' This group had been labeled ''Savant'' in our previous study [13]. Note that each cluster contains an independent cohort of cases. B) Principal components analysis (PCA) of the individuals based on the 123 ADI-R severity scores. Each subgroup of individuals identified in (A) is assigned a color, which identifies individuals from that subgroup in the PCA. Red: Language-impaired; Green: Intermediate; Yellow: Moderate; Blue: Mild. Each point on the PCA represents an individual with an ASD whose position is defined by his/her scores for the 123 ADI-R items. (TIF) Figure S3 Network connections centered on HTR4 from

Author Contributions
Conceived and designed the experiments: VWH. Performed the experiments: VWH. Analyzed the data: VWH. Wrote the paper: VWH. Made suggestions and provided advice on the study, reviewed and edited the manuscript: AA. Involved in preliminary quantitative trait association analyses related to this study: AH.