James M. Jaworski

A Genome-wide Association Study of Autism Reveals a Common Novel Risk Locus at 5p14.1

Although autism is one of the most heritable neuropsychiatric disorders, its underlying genetic architecture has largely eluded description. To comprehensively examine the hypothesis that common variation is important in autism, we performed a genome‐wide association study (GWAS) using a discovery dataset of 438 autistic Caucasian families and the Illumina Human 1M beadchip. 96 single nucleotide polymorphisms (SNPs) demonstrated strong association with autism risk (p‐value < 0.0001). The validation of the top 96 SNPs was performed using an independent dataset of 487 Caucasian autism families genotyped on the 550K Illumina BeadChip. A novel region on chromosome 5p14.1 showed significance in both the discovery and validation datasets. Joint analysis of all SNPs in this region identified 8 SNPs having improved p‐values (3.24E‐04 to 3.40E‐06) than in either dataset alone. Our findings demonstrate that in addition to multiple rare variations, part of the complex genetic architecture of autism involves common variation.

Analysis of the RELN gene as a genetic risk factor for autism

Several genome-wide screens have indicated the presence of an autism susceptibility locus within the distal long arm of chromosome 7 (7q). Mapping at 7q22 within this region is the candidate gene reelin (RELN). RELN encodes a signaling protein that plays a pivotal role in the migration of several neuronal cell types and in the development of neural connections. Given these neurodevelopmental functions, recent reports that RELN influences genetic risk for autism are of significant interest. The total data set consists of 218 Caucasian families collected by our group, 85 Caucasian families collected by AGRE, and 68 Caucasian families collected at Tufts University were tested for genetic association of RELN variants to autism. Markers included five single-nucleotide polymorphisms (SNPs) and a repeat in the 5′-untranslated region (5′-UTR). Tests for association in Duke and AGRE families were also performed on four additional SNPs in the genes PSMC2 and ORC5L, which flank RELN. Family-based association analyses (PDT, Geno-PDT, and FBAT) were used to test for association of single-locus markers and multilocus haplotypes with autism. The most significant association identified from this combined data set was for the 5′-UTR repeat (PDT P-value=0.002). These analyses show the potential of RELN as an important contributor to genetic risk in autism.

A noise-reduction GWAS analysis implicates altered regulation of neurite outgrowth and guidance in autism

BackgroundGenome-wide Association Studies (GWAS) have proved invaluable for the identification of disease susceptibility genes. However, the prioritization of candidate genes and regions for follow-up studies often proves difficult due to false-positive associations caused by statistical noise and multiple-testing. In order to address this issue, we propose the novel GWAS noise reduction (GWAS-NR) method as a way to increase the power to detect true associations in GWAS, particularly in complex diseases such as autism.MethodsGWAS-NR utilizes a linear filter to identify genomic regions demonstrating correlation among association signals in multiple datasets. We used computer simulations to assess the ability of GWAS-NR to detect association against the commonly used joint analysis and Fishers methods. Furthermore, we applied GWAS-NR to a family-based autism GWAS of 597 families and a second existing autism GWAS of 696 families from the Autism Genetic Resource Exchange (AGRE) to arrive at a compendium of autism candidate genes. These genes were manually annotated and classified by a literature review and functional grouping in order to reveal biological pathways which might contribute to autism aetiology.ResultsComputer simulations indicate that GWAS-NR achieves a significantly higher classification rate for true positive association signals than either the joint analysis or Fishers methods and that it can also achieve this when there is imperfect marker overlap across datasets or when the closest disease-related polymorphism is not directly typed. In two autism datasets, GWAS-NR analysis resulted in 1535 significant linkage disequilibrium (LD) blocks overlapping 431 unique reference sequencing (RefSeq) genes. Moreover, we identified the nearest RefSeq gene to the non-gene overlapping LD blocks, producing a final candidate set of 860 genes. Functional categorization of these implicated genes indicates that a significant proportion of them cooperate in a coherent pathway that regulates the directional protrusion of axons and dendrites to their appropriate synaptic targets.ConclusionsAs statistical noise is likely to particularly affect studies of complex disorders, where genetic heterogeneity or interaction between genes may confound the ability to detect association, GWAS-NR offers a powerful method for prioritizing regions for follow-up studies. Applying this method to autism datasets, GWAS-NR analysis indicates that a large subset of genes involved in the outgrowth and guidance of axons and dendrites is implicated in the aetiology of autism.

Evaluation of copy number variations reveals novel candidate genes in autism spectrum disorder-associated pathways

Autism spectrum disorders (ASDs) are highly heritable, yet relatively few associated genetic loci have been replicated. Copy number variations (CNVs) have been implicated in autism; however, the majority of loci contribute to <1% of the disease population. Therefore, independent studies are important to refine associated CNV regions and discover novel susceptibility genes. In this study, a genome-wide SNP array was utilized for CNV detection by two distinct algorithms in a European ancestry case-control data set. We identify a significantly higher burden in the number and size of deletions, and disrupting more genes in ASD cases. Moreover, 18 deletions larger than 1 Mb were detected exclusively in cases, implicating novel regions at 2q22.1, 3p26.3, 4q12 and 14q23. Case-specific CNVs provided further evidence for pathways previously implicated in ASDs, revealing new candidate genes within the GABAergic signaling and neural development pathways. These include DBI, an allosteric binder of GABA receptors, GABARAPL1, the GABA receptor-associated protein, and SLC6A11, a postsynaptic GABA transporter. We also identified CNVs in COBL, deletions of which cause defects in neuronal cytoskeleton morphogenesis in model vertebrates, and DNER, a neuron-specific Notch ligand required for cerebellar development. Moreover, we found evidence of genetic overlap between ASDs and other neurodevelopmental and neuropsychiatric diseases. These genes include glutamate receptors (GRID1, GRIK2 and GRIK4), synaptic regulators (NRXN3, SLC6A8 and SYN3), transcription factor (ZNF804A) and RNA-binding protein FMR1. Taken together, these CNVs may be a few of the missing pieces of ASD heritability and lead to discovering novel etiological mechanisms.

An analysis paradigm for investigating multi-locus effects in complex disease: examination of three GABA receptor subunit genes on 15q11-q13 as risk factors for autistic disorder.

Gene‐gene interactions are likely involved in many complex genetic disorders and new statistical approaches for detecting such interactions are needed. We propose a multi‐analytic paradigm, relying on convergence of evidence across multiple analysis tools. Our paradigm tests for main and interactive effects, through allele, genotype and haplotype association. We applied our paradigm to genotype data from three GABAA receptor subunit genes (GABRB3, GABRA5, and GABRG3) on chromosome 15 in 470 Caucasian autism families. Previously implicated in autism, we hypothesized these genes interact to contribute to risk. We detected no evidence of main effects by allelic (PDT, FBAT) or genotypic (genotype‐PDT) association at individual markers. However, three two‐marker haplotypes in GABRG3 were significant (HBAT). We detected no significant multi‐locus associations using genotype‐PDT analysis or the EMDR data reduction program. However, consistent with the haplotype findings, the best single locus EMDR model selected a GABRG3 marker. Further, the best pairwise genotype‐PDT result involved GABRB3 and GABRG3, and all multi‐locus EMDR models also selected GABRB3 and GABRG3 markers. GABA receptor subunit genes do not significantly interact to contribute to autism risk in our overall data set. However, the consistency of results across analyses suggests that we have defined a useful framework for evaluating gene‐gene interactions.

Copy number variants in extended autism spectrum disorder families reveal candidates potentially involved in autism risk

Copy number variations (CNVs) are a major cause of genetic disruption in the human genome with far more nucleotides being altered by duplications and deletions than by single nucleotide polymorphisms (SNPs). In the multifaceted etiology of autism spectrum disorders (ASDs), CNVs appear to contribute significantly to our understanding of the pathogenesis of this complex disease. A unique resource of 42 extended ASD families was genotyped for over 1 million SNPs to detect CNVs that may contribute to ASD susceptibility. Each family has at least one avuncular or cousin pair with ASD. Families were then evaluated for co-segregation of CNVs in ASD patients. We identified a total of five deletions and seven duplications in eleven families that co-segregated with ASD. Two of the CNVs overlap with regions on 7p21.3 and 15q24.1 that have been previously reported in ASD individuals and two additional CNVs on 3p26.3 and 12q24.32 occur near regions associated with schizophrenia. These findings provide further evidence for the involvement of ICA1 and NXPH1 on 7p21.3 in ASD susceptibility and highlight novel ASD candidates, including CHL1, FGFBP3 and POUF41. These studies highlight the power of using extended families for gene discovery in traits with a complex etiology.

Association and gene-gene interaction of SLC6A4 and ITGB3 in autism.

Autism is a heritable neurodevelopmental disorder with substantial genetic heterogeneity. Studies point to possible links between autism and two serotonin related genes: SLC6A4 and ITGB3 with a sex‐specific genetic effect and interaction between the genes. Despite positive findings, inconsistent results have complicated interpretation. This study seeks to validate and clarify previous findings in an independent dataset taking into account sex, family‐history (FH) and gene–gene effects. Family‐based association analysis was performed within each gene. Gene–gene interactions were tested using extended multifactor dimensionality reduction (EMDR) and MDR‐phenomics (MDR‐P) using sex of affecteds and FH as covariates. No significant associations with individual SNPs were found in the datasets stratified by sex, but associations did emerge when we stratified by family history. While not significant in the overall dataset, nominally significant association was identified at RS2066713 (P = 0.006) within SLC6A4 in family‐history negative (FH−) families, at RS2066713 (P = 0.038) in family‐history positive (FH+) families but with the opposite risk allele as in the FH− families. For ITGB3, nominally significant association was identified at RS3809865 overall (P = 0.040) and within FH+ families (P = 0.031). However, none of the associations survived the multiple testing correction. MDR‐P confirmed gene–gene effects using sex of affecteds (P = 0.023) and family history (P = 0.014, survived the multiple testing corrections) as covariates. Our results indicate the extensive heterogeneity within these two genes among families. The potential interaction between SLC6A4 and ITGB3 may be clarified using family history as an indicator of genetic architecture, illustrating the importance of covariates as markers of heterogeneity in genetic analyses.

Dissecting the locus heterogeneity of autism: significant linkage to chromosome 12q14

Autism is a common neurodevelopmental disorder with a significant genetic component and locus heterogeneity. To date, 12 microsatellite genome screens have been performed using various data sets of sib-pair families (parents and affected children) resulting in numerous regions of potential linkage across the genome. However, no universal region or consistent candidate gene from these regions has emerged. The use of large, extended pedigrees is a recognized powerful approach to identify significant linkage results, as these families potentially contain more potential linkage information than sib-pair families. A genome-wide linkage analysis was performed on 26 extended autism families (65 affected, 184 total individuals). Each family had two to four affected individuals comprised of either avuncular or cousin pairs. For analysis, we used a high-density single-nucleotide polymorphism genotyping assay, the Affymetrix GeneChip Human Mapping 10K array. Two-point analysis gave peak heterogeneity limit of detection (HLOD) of 2.82 at rs2877739 on chromosome 14q. Suggestive linkage evidence (HLOD>2) from a two-point analysis was also found on chromosomes 1q, 2q, 5q, 6p,11q and 12q. Chromosome 12q was the only region showing significant linkage evidence by multipoint analysis with a peak HLOD=3.02 at rs1445442. In addition, this linkage evidence was enhanced significantly in the families with only male affected (multipoint HLOD=4.51), suggesting a significant gender-specific effect in the etiology of autism. Chromosome-wide haplotype analyses on chromosome 12 localized the potential autism gene to a 4 cM region shared among the affected individuals across linked families. This novel linkage peak on chromosome 12q further supports the hypothesis of substantial locus heterogeneity in autism.

Investigation of potential gene-gene interactions between APOE and RELN contributing to autism risk

Background Several candidate gene studies support RELN as susceptibility gene for autism. Given the complex inheritance pattern of autism, it is expected that gene–gene interactions will exist. A logical starting point for examining potential gene–gene interactions is to evaluate the joint effects of genes involved in a common biological pathway. RELN shares a common biological pathway with APOE, and Persico et al. have observed transmission distortion of the APOE-2 allele in autism families. Objective We evaluated RELN and APOE for joint effects in autism susceptibility. Methods A total of 470 Caucasian autism families were analyzed (265 multiplex; 168 trios with no family history; 37 positive family history but only one sampled affected). These families were genotyped for 11 RELN polymorphisms, including the 5′ untranslated region repeat previously associated with autism, as well as for the APOE functional allele. We evaluated single locus allelic and genotypic association with the pedigree disequilibrium test and geno-PDT, respectively. Multilocus effects were evaluated using the extended version of the multifactorial dimensionality reduction method. Results For the single locus analyses, there was no evidence for an effect of APOE in our data set. Evidence for association with RELN (rs2073559; trio subset P=0.07 PDT; P=0.001 geno-PDT; overall geno-PDT P=0.05), however, was found. For multilocus geno-PDT analysis, the joint genotype of APOE and RELN rs2073559 was highly significant (trio subset, global P=0.0001), probably driven by the RELN single locus effect. Using the extended version of the multifactorial dimensionality reduction method to detect multilocus effects, there were no statistically significant associations for any of the n-locus combinations involving RELN or APOE in the overall or multiplex subset. In the trio subset, 1-locus and 2-locus models selected only markers in RELN as best models for predicting autism case status. Conclusion Thus, we conclude that there is no main effect of APOE in our autism data set, nor is there any evidence for a joint effect of APOE with RELN. RELN, however, remains a good candidate for autism susceptibility.

Analysis of the autism chromosome 2 linkage region: GAD1 and other candidate genes

Autism has a strong and complex genetic component, involving several genes. Genomic screens, including our own, have shown suggestive evidence for linkage over a 20-30 cM region on chromosome 2q31-q33. Two subsequent reports showed that the linkage evidence increased in the subset of families with phrase speech delay (PSD), defined as onset of phrase speech later than 3 years of age. To further investigate the linkage in the presumptive candidate region, microsatellite markers in a 2 cM grid covering the interval from 164 to 203 cM were analyzed in 110 multiplex (2 or more sampled autism patients) families. A maximum heterogeneity LOD (HLOD) score of 1.54 was detected at D2S1776 (173 cM) in the overall dataset (dominant model), increasing to 1.71 in the PSD subset. While not conclusive, these data continue to provide suggestive evidence for linkage, particularly considering replication by multiple independent groups. Positive LOD scores extended over the entire region, continuing to define a broad candidate interval. Association studies were performed on several functional candidates mapping within the region. These included GAD1, encoding GAD67, whose levels are reduced in autopsy brain material from autistic subjects, and STK17B, ABI2, CTLA4, CD28, NEUROD1, PDE1A, HOXD1 and DLX2. We found no evidence for significant allelic association between autism and any of these candidates, suggesting that they do not play a major role in the genetics of autism or that substantial allelic heterogeneity at any one of these loci dilutes potential disease-allele association.