Convergent Pathways in Idiopathic Autism Revealed by Time Course Transcriptomic Analysis of Patient-Derived Neurons

Potentially pathogenic alterations have been identified in individuals with autism spectrum disorders (ASDs) within a variety of key neurodevelopment genes. While this hints at a common ASD molecular etiology, gaps persist in our understanding of the neurodevelopmental mechanisms impacted by genetic variants enriched in ASD patients. Induced pluripotent stem cells (iPSCs) can model neurodevelopment in vitro, permitting the characterization of pathogenic mechanisms that manifest during corticogenesis. Taking this approach, we examined the transcriptional differences between iPSC-derived cortical neurons from patients with idiopathic ASD and unaffected controls over a 135-day course of neuronal differentiation. Our data show ASD-specific misregulation of genes involved in neuronal differentiation, axon guidance, cell migration, DNA and RNA metabolism, and neural region patterning. Furthermore, functional analysis revealed defects in neuronal migration and electrophysiological activity, providing compelling support for the transcriptome analysis data. This study reveals important and functionally validated insights into common processes altered in early neuronal development and corticogenesis and may contribute to ASD pathogenesis.

Identification of rare noncoding sequence variants in gamma-aminobutyric acid A receptor, alpha 4 subunit in autism spectrum disorder

Alterations of the gamma-aminobutyric acid (GABA) signaling system has been strongly linked to the pathophysiology of autism spectrum disorder (ASD). Genetic associations of common variants in GABA receptor subunits, in particular GABRA4 on chromosome 4p12, with ASD have been replicated by several studies. Moreover, molecular investigations have identified altered transcriptional and translational levels of this gene and protein in brains of ASD individuals. Since the genotyped common variants are likely not the functional variants contributing to the molecular consequences or underlying ASD phenotype, this study aims to examine rare sequence variants in GABRA4, including those outside the protein coding regions of the gene. We comprehensively re-sequenced the entire protein coding and noncoding portions of the gene and putative regulatory sequences in 82 ASD individuals and 55 developmentally typical pediatric controls, all homozygous for the most significant previously associated ASD risk allele (G/G at rs1912960). We identified only a single common, coding variant, and no association of any single marker or set of variants with ASD. Functional annotation of noncoding variants identified several rare variants in putative regulatory sites. Finally, a rare variant unique to ASD cases, in an evolutionary conserved site of the 3'UTR, shows a trend toward decreasing gene expression. Hence, GABRA4 rare variants in noncoding DNA may be variants of modest physiological effects in ASD etiology.

Early-Onset Alzheimer Disease and Candidate Risk Genes Involved in Endolysosomal Transport

Importance

Mutations in APP, PSEN1, and PSEN2 lead to early-onset Alzheimer disease (EOAD) but account for only approximately 11% of EOAD overall, leaving most of the genetic risk for the most severe form of Alzheimer disease unexplained. This extreme phenotype likely harbors highly penetrant risk variants, making it primed for discovery of novel risk genes and pathways for AD.

Objective

To search for rare variants contributing to the risk for EOAD.

Design, Setting, and Participants

In this case-control study, whole-exome sequencing (WES) was performed in 51 non-Hispanic white (NHW) patients with EOAD (age at onset <65 years) and 19 Caribbean Hispanic families previously screened as negative for established APP, PSEN1, and PSEN2 causal variants. Participants were recruited from John P. Hussman Institute for Human Genomics, Case Western Reserve University, and Columbia University. Rare, deleterious, nonsynonymous, or loss-of-function variants were filtered to identify variants in known and suspected AD genes, variants in multiple unrelated NHW patients, variants present in 19 Hispanic EOAD WES families, and genes with variants in multiple unrelated NHW patients. These variants/genes were tested for association in an independent cohort of 1524 patients with EOAD, 7046 patients with late-onset AD (LOAD), and 7001 cognitively intact controls (age at examination, >65 years) from the Alzheimer's Disease Genetics Consortium. The study was conducted from January 21, 2013, to October 13, 2016.

Main Outcomes and Measures

Alzheimer disease diagnosed according to standard National Institute of Neurological and Communicative Disorders and Stroke and the Alzheimer Disease and Related Disorders Association criteria. Association between Alzheimer disease and genetic variants and genes was measured using logistic regression and sequence kernel association test-optimal gene tests, respectively.

Results

Of the 1524 NHW patients with EOAD, 765 (50.2%) were women and mean (SD) age was 60.0 (4.9) years; of the 7046 NHW patients with LOAD, 4171 (59.2%) were women and mean (SD) age was 77.4 (8.6) years; and of the 7001 NHW controls, 4215 (60.2%) were women and mean (SD) age was 77.4 (8.6) years. The gene PSD2, for which multiple unrelated NHW cases had rare missense variants, was significantly associated with EOAD (P = 2.05 × 10-6; Bonferroni-corrected P value [BP] = 1.3 × 10-3) and LOAD (P = 6.22 × 10-6; BP = 4.1 × 10-3). A missense variant in TCIRG1, present in a NHW patient and segregating in 3 cases of a Hispanic family, was more frequent in EOAD cases (odds ratio [OR], 2.13; 95% CI, 0.99-4.55; P = .06; BP = 0.413), and significantly associated with LOAD (OR, 2.23; 95% CI, 1.37-3.62; P = 7.2 × 10-4; BP = 5.0 × 10-3). A missense variant in the LOAD risk gene RIN3 showed suggestive evidence of association with EOAD after Bonferroni correction (OR, 4.56; 95% CI, 1.26-16.48; P = .02, BP = 0.091). In addition, a missense variant in RUFY1 identified in 2 NHW EOAD cases showed suggestive evidence of an association with EOAD as well (OR, 18.63; 95% CI, 1.62-213.45; P = .003; BP = 0.129).

Conclusions and Relevance

The genes PSD2, TCIRG1, RIN3, and RUFY1 all may be involved in endolysosomal transport-a process known to be important to development of AD. Furthermore, this study identified shared risk genes between EOAD and LOAD similar to previously reported genes, such as SORL1, PSEN2, and TREM2.

Exome Sequencing of Extended Families with Alzheimer's Disease Identifies Novel Genes Implicated in Cell Immunity and Neuronal Function

Objective

Alzheimer’s disease (AD) is a neurodegenerative disorder for which more than 20 genetic loci have been implicated to date. However, studies demonstrate not all genetic factors have been identified. Therefore, in this study we seek to identify additional rare variants and novel genes potentially contributing to AD.

Methods

Whole exome sequencing was performed on 23 multi-generational families with an average of eight affected subjects. Exome sequencing was filtered for rare, nonsynonymous and loss-of-function variants. Alterations predicted to have a functional consequence and located within either a previously reported AD gene, a linkage peak (LOD>2), or clustering in the same gene across multiple families, were prioritized.

Results

Rare variants were found in known AD risk genes including AKAP9, CD33, CR1, EPHA1, INPP5D, NME8, PSEN1, SORL1, TREM2 and UNC5C. Three families had five variants of interest in linkage regions with LOD>2. Genes with segregating alterations in these peaks include CD163L1 and CLECL1, two genes that have both been implicated in immunity, CTNNA1, which encodes a catenin in the cerebral cortex and MIEF1, a gene that may induce mitochondrial dysfunction and has the potential to damage neurons. Four genes were identified with alterations in more than one family include PLEKHG5, a gene that causes Charcot-Marie-Tooth disease and THBS2, which promotes synaptogenesis.

Conclusion

Utilizing large families with a heavy burden of disease allowed for the identification of rare variants co-segregating with disease. Variants were identified in both known AD risk genes and in novel genes.

SORL1 mutations in early- and late-onset Alzheimer disease

Objective:

To characterize the clinical and molecular effect of mutations in the sortilin-related receptor (SORL1) gene.

Methods:

We performed whole-exome sequencing in early-onset Alzheimer disease (EOAD) and late-onset Alzheimer disease (LOAD) families followed by functional studies of select variants. The phenotypic consequences associated with SORL1 mutations were characterized based on clinical reviews of medical records. Functional studies were completed to evaluate β-amyloid (Aβ) production and amyloid precursor protein (APP) trafficking associated with SORL1 mutations.

Results:

SORL1 alterations were present in 2 EOAD families. In one, a SORL1 T588I change was identified in 4 individuals with AD, 2 of whom had parkinsonian features. In the second, an SORL1 T2134 alteration was found in 3 of 4 AD cases, one of whom had postmortem Lewy bodies. Among LOAD cases, 4 individuals with either SORL1 A528T or T947M alterations had parkinsonian features. Functionally, the variants weaken the interaction of the SORL1 protein with full-length APP, altering levels of Aβ and interfering with APP trafficking.

Conclusions:

The findings from this study support an important role for SORL1 mutations in AD pathogenesis by way of altering Aβ levels and interfering with APP trafficking. In addition, the presence of parkinsonian features among select individuals with AD and SORL1 mutations merits further investigation.

ABCA7 frameshift deletion associated with Alzheimer disease in African Americans

Objective:

To identify a causative variant(s) that may contribute to Alzheimer disease (AD) in African Americans (AA) in the ATP-binding cassette, subfamily A (ABC1), member 7 (ABCA7) gene, a known risk factor for late-onset AD.

Methods:

Custom capture sequencing was performed on ∼150 kb encompassing ABCA7 in 40 AA cases and 37 AA controls carrying the AA risk allele (rs115550680). Association testing was performed for an ABCA7 deletion identified in large AA data sets (discovery n = 1,068; replication n = 1,749) and whole exome sequencing of Caribbean Hispanic (CH) AD families.

Results:

A 44-base pair deletion (rs142076058) was identified in all 77 risk genotype carriers, which shows that the deletion is in high linkage disequilibrium with the risk allele. The deletion was assessed in a large data set (531 cases and 527 controls) and, after adjustments for age, sex, and APOE status, was significantly associated with disease (p = 0.0002, odds ratio [OR] = 2.13 [95% confidence interval (CI): 1.42–3.20]). An independent data set replicated the association (447 cases and 880 controls, p = 0.0117, OR = 1.65 [95% CI: 1.12–2.44]), and joint analysis increased the significance (p = 1.414 × 10⁻⁵, OR = 1.81 [95% CI: 1.38–2.37]). The deletion is common in AA cases (15.2%) and AA controls (9.74%), but in only 0.12% of our non-Hispanic white cohort. Whole exome sequencing of multiplex, CH families identified the deletion cosegregating with disease in a large sibship. The deleted allele produces a stable, detectable RNA strand and is predicted to result in a frameshift mutation (p.Arg578Alafs) that could interfere with protein function.

Conclusions:

This common ABCA7 deletion could represent an ethnic-specific pathogenic alteration in AD.

Dopaminergic variants in siblings at high risk for autism: Associations with initiating joint attention

Younger siblings of children with autism spectrum disorder (ASD; high-risk siblings) exhibit lower levels of initiating joint attention (IJA; sharing an object or experience with a social partner through gaze and/or gesture) than low-risk siblings of children without ASD. However, high-risk siblings also exhibit substantial variability in this domain. The neurotransmitter dopamine is linked to brain areas associated with reward, motivation, and attention, and common dopaminergic variants have been associated with attention difficulties. We examined whether these common dopaminergic variants, DRD4 and DRD2, explain variability in IJA in high-risk (n = 55) and low-risk (n = 38) siblings. IJA was assessed in the first year during a semi-structured interaction with an examiner. DRD4 and DRD2 genotypes were coded according to associated dopaminergic functioning to create a gene score, with higher scores indicating more genotypes associated with less efficient dopaminergic functioning. Higher dopamine gene scores (indicative of less efficient dopaminergic functioning) were associated with lower levels of IJA in the first year for high-risk siblings, while the opposite pattern emerged in low-risk siblings. Findings suggest differential susceptibility-IJA was differentially associated with dopaminergic functioning depending on familial ASD risk. Understanding genes linked to ASD-relevant behaviors in high-risk siblings will aid in early identification of children at greatest risk for difficulties in these behavioral domains, facilitating targeted prevention and intervention. Autism Res 2016, 9: 1142-1150. © 2016 International Society for Autism Research, Wiley Periodicals, Inc.

Global and local ancestry in African-Americans: Implications for Alzheimer's disease risk

INTRODUCTION

African-American (AA) individuals have a higher risk for late-onset Alzheimer's disease (LOAD) than Americans of primarily European ancestry (EA). Recently, the largest genome-wide association study in AAs to date confirmed that six of the Alzheimer's disease (AD)-related genetic variants originally discovered in EA cohorts are also risk variants in AA; however, the risk attributable to many of the loci (e.g., APOE, ABCA7) differed substantially from previous studies in EA. There likely are risk variants of higher frequency in AAs that have not been discovered.

METHODS

We performed a comprehensive analysis of genetically determined local and global ancestry in AAs with regard to LOAD status.

RESULTS

Compared to controls, LOAD cases showed higher levels of African ancestry, both globally and at several LOAD relevant loci, which explained risk for AD beyond global differences.

DISCUSSION

Exploratory post hoc analyses highlight regions with greatest differences in ancestry as potential candidate regions for future genetic analyses.

Segregation of a rare TTC3 variant in an extended family with late-onset Alzheimer disease

OBJECTIVE

The genetic risk architecture of Alzheimer disease (AD) is complex with single pathogenic mutations leading to early-onset AD, while both rare and common genetic susceptibility variants contribute to the more widespread late-onset AD (LOAD); we sought to discover novel genes contributing to LOAD risk.

METHODS

Whole-exome sequencing and genome-wide genotyping were performed on 11 affected individuals in an extended family with an apparent autosomal dominant pattern of LOAD. Variants of interest were then evaluated in a large cohort of LOAD cases and aged controls.

RESULTS

We detected a single rare, nonsynonymous variant shared in all 11 LOAD individuals, a missense change in the tetratricopeptide repeat domain 3 (TTC3) gene. The missense variant, rs377155188 (p.S1038C), is predicted to be damaging. Affecteds-only multipoint linkage analysis demonstrated that this region of TTC3 has a LOD score of 2.66 in this family.

CONCLUSION

The TTC3 p.S1038C substitution may represent a segregating, rare LOAD risk variant. Previous studies have shown that TTC3 expression is consistently reduced in LOAD patients and negatively correlated with AD neuropathology and that TTC3 is a regulator of Akt signaling, a key pathway disrupted in LOAD. This study demonstrates how utilizing whole-exome sequencing in a large, multigenerational family with a high incidence of LOAD could reveal a novel candidate gene.

Targeted massively parallel sequencing of autism spectrum disorder-associated genes in a case control cohort reveals rare loss-of-function risk variants

Background

Autism spectrum disorder (ASD) is highly heritable, yet genome-wide association studies (GWAS), copy number variation screens, and candidate gene association studies have found no single factor accounting for a large percentage of genetic risk. ASD trio exome sequencing studies have revealed genes with recurrent de novo loss-of-function variants as strong risk factors, but there are relatively few recurrently affected genes while as many as 1000 genes are predicted to play a role. As such, it is critical to identify the remaining rare and low-frequency variants contributing to ASD.

Methods

We have utilized an approach of prioritization of genes by GWAS and follow-up with massively parallel sequencing in a case-control cohort. Using a previously reported ASD noise reduction GWAS analyses, we prioritized 837 RefSeq genes for custom targeting and sequencing. We sequenced the coding regions of those genes in 2071 ASD cases and 904 controls of European white ancestry. We applied comprehensive annotation to identify single variants which could confer ASD risk and also gene-based association analysis to identify sets of rare variants associated with ASD.

Results

We identified a significant over-representation of rare loss-of-function variants in genes previously associated with ASD, including a de novo premature stop variant in the well-established ASD candidate gene RBFOX1. Furthermore, ASD cases were more likely to have two damaging missense variants in candidate genes than controls. Finally, gene-based rare variant association implicates genes functioning in excitatory neurotransmission and neurite outgrowth and guidance pathways including CACNAD2, KCNH7, and NRXN1.

Conclusions

We find suggestive evidence that rare variants in synaptic genes are associated with ASD and that loss-of-function mutations in ASD candidate genes are a major risk factor, and we implicate damaging mutations in glutamate signaling receptors and neuronal adhesion and guidance molecules. Furthermore, the role of de novo mutations in ASD remains to be fully investigated as we identified the first reported protein-truncating variant in RBFOX1 in ASD. Overall, this work, combined with others in the field, suggests a convergence of genes and molecular pathways underlying ASD etiology.

Electronic supplementary material

The online version of this article (doi:10.1186/s13229-015-0034-z) contains supplementary material, which is available to authorized users.

Two knockdown models of the autism genes SYNGAP1 and SHANK3 in zebrafish produce similar behavioral phenotypes associated with embryonic disruptions of brain morphogenesis

Despite significant progress in the genetics of autism spectrum disorder (ASD), how genetic mutations translate to the behavioral changes characteristic of ASD remains largely unknown. ASD affects 1-2% of children and adults, and is characterized by deficits in verbal and non-verbal communication, and social interactions, as well as the presence of repetitive behaviors and/or stereotyped interests. ASD is clinically and etiologically heterogeneous, with a strong genetic component. Here, we present functional data from syngap1 and shank3 zebrafish loss-of-function models of ASD. SYNGAP1, a synaptic Ras GTPase activating protein, and SHANK3, a synaptic scaffolding protein, were chosen because of mounting evidence that haploinsufficiency in these genes is highly penetrant for ASD and intellectual disability (ID). Orthologs of both SYNGAP1 and SHANK3 are duplicated in the zebrafish genome and we find that all four transcripts (syngap1a, syngap1b, shank3a and shank3b) are expressed at the earliest stages of nervous system development with pronounced expression in the larval brain. Consistent with early expression of these genes, knockdown of syngap1b or shank3a cause common embryonic phenotypes including delayed mid- and hindbrain development, disruptions in motor behaviors that manifest as unproductive swim attempts, and spontaneous, seizure-like behaviors. Our findings indicate that both syngap1b and shank3a play novel roles in morphogenesis resulting in common brain and behavioral phenotypes.

Exome sequencing of extended families with autism reveals genes shared across neurodevelopmental and neuropsychiatric disorders

BACKGROUND

Autism spectrum disorders (ASDs) comprise a range of neurodevelopmental conditions of varying severity, characterized by marked qualitative difficulties in social relatedness, communication, and behavior. Despite overwhelming evidence of high heritability, results from genetic studies to date show that ASD etiology is extremely heterogeneous and only a fraction of autism genes have been discovered.

METHODS

To help unravel this genetic complexity, we performed whole exome sequencing on 100 ASD individuals from 40 families with multiple distantly related affected individuals. All families contained a minimum of one pair of ASD cousins. Each individual was captured with the Agilent SureSelect Human All Exon kit, sequenced on the Illumina Hiseq 2000, and the resulting data processed and annotated with Burrows-Wheeler Aligner (BWA), Genome Analysis Toolkit (GATK), and SeattleSeq. Genotyping information on each family was utilized in order to determine genomic regions that were identical by descent (IBD). Variants identified by exome sequencing which occurred in IBD regions and present in all affected individuals within each family were then evaluated to determine which may potentially be disease related. Nucleotide alterations that were novel and rare (minor allele frequency, MAF, less than 0.05) and predicted to be detrimental, either by altering amino acids or splicing patterns, were prioritized.

RESULTS

We identified numerous potentially damaging, ASD associated risk variants in genes previously unrelated to autism. A subset of these genes has been implicated in other neurobehavioral disorders including depression (SLIT3), epilepsy (CLCN2, PRICKLE1), intellectual disability (AP4M1), schizophrenia (WDR60), and Tourette syndrome (OFCC1). Additional alterations were found in previously reported autism candidate genes, including three genes with alterations in multiple families (CEP290, CSMD1, FAT1, and STXBP5). Compiling a list of ASD candidate genes from the literature, we determined that variants occurred in ASD candidate genes 1.65 times more frequently than in random genes captured by exome sequencing (P = 8.55 × 10-5).

CONCLUSIONS

By studying these unique pedigrees, we have identified novel DNA variations related to ASD, demonstrated that exome sequencing in extended families is a powerful tool for ASD candidate gene discovery, and provided further evidence of an underlying genetic component to a wide range of neurodevelopmental and neuropsychiatric diseases.

Variation in oxytocin receptor gene (OXTR) polymorphisms is associated with emotional and behavioral reactions to betrayal

Variations in the gene that encodes the oxytocin receptor (OXTR) have been associated with many aspects of social cognition as well as several prosocial behaviors. However, potential associations of OXTR variants with reactions to betrayals of trust while cooperating for mutual benefit have not yet been explored. We examined how variations in 10 single-nucleotide polymorphisms on OXTR were associated with behavior and emotional reactions after a betrayal of trust in an iterated Prisoner's Dilemma Game. After correction for multiple testing, one haplotype (C-rs9840864, T-rs2268494) was significantly associated with faster retaliation post-betrayal-an association that appeared to be due to this haplotype's intermediate effect of exacerbating people's anger after they had been betrayed. Furthermore, a second haplotype (A-rs237887, C-rs2268490) was associated with higher levels of post-betrayal satisfaction, and a third haplotype (G-rs237887, C-rs2268490) was associated with lower levels of post-betrayal satisfaction.

Genome-wide association and linkage study in the Amish detects a novel candidate late-onset Alzheimer disease gene

To identify novel late-onset Alzheimer disease (LOAD) risk genes, we have analysed Amish populations of Ohio and Indiana. We performed genome-wide SNP linkage and association studies on 798 individuals (109 with LOAD). We tested association using the Modified Quasi-Likelihood Score test and also performed two-point and multipoint linkage analyses. We found that LOAD was significantly associated with APOE (P= 9.0 × 10-6) in all our ascertainment regions except for the Adams County, Indiana, community (P= 0.55). Genome-wide, the most strongly associated SNP was rs12361953 (P= 7.92 × 10-7). A very strong, genome-wide significant multipoint peak [recessive heterogeneity multipoint LOD (HLOD) = 6.14, dominant HLOD = 6.05] was detected on 2p12. Three additional loci with multipoint HLOD scores >3 were detected on 3q26, 9q31 and 18p11. Converging linkage and association results, the most significantly associated SNP under the 2p12 peak was at rs2974151 (P= 1.29 × 10-4). This SNP is located in CTNNA2, which encodes catenin alpha 2, a neuronal-specific catenin known to have function in the developing brain. These results identify CTNNA2 as a novel candidate LOAD gene, and implicate three other regions of the genome as novel LOAD loci. These results underscore the utility of using family-based linkage and association analyses in isolated populations to identify novel loci for traits with complex genetic architecture.

Evaluating mitochondrial DNA variation in autism spectrum disorders

Despite the increasing speculation that oxidative stress and abnormal energy metabolism may play a role in Autism Spectrum Disorders (ASD), and the observation that patients with mitochondrial defects have symptoms consistent with ASD, there are no comprehensive published studies examining the role of mitochondrial variation in autism. Therefore, we have sought to comprehensively examine the role of mitochondrial DNA (mtDNA) variation with regard to ASD risk, employing a multi-phase approach. In phase 1 of our experiment, we examined 132 mtDNA single-nucleotide polymorphisms (SNPs) genotyped as part of our genome-wide association studies of ASD. In phase 2 we genotyped the major European mitochondrial haplogroup-defining variants within an expanded set of autism probands and controls. Finally in phase 3, we resequenced the entire mtDNA in a subset of our Caucasian samples (∼400 proband-father pairs). In each phase we tested whether mitochondrial variation showed evidence of association to ASD. Despite a thorough interrogation of mtDNA variation, we found no evidence to suggest a major role for mtDNA variation in ASD susceptibility. Accordingly, while there may be attractive biological hints suggesting the role of mitochondria in ASD our data indicate that mtDNA variation is not a major contributing factor to the development of ASD.

Individual common variants exert weak effects on the risk for autism spectrum disorders

While it is apparent that rare variation can play an important role in the genetic architecture of autism spectrum disorders (ASDs), the contribution of common variation to the risk of developing ASD is less clear. To produce a more comprehensive picture, we report Stage 2 of the Autism Genome Project genome-wide association study, adding 1301 ASD families and bringing the total to 2705 families analysed (Stages 1 and 2). In addition to evaluating the association of individual single nucleotide polymorphisms (SNPs), we also sought evidence that common variants, en masse, might affect the risk. Despite genotyping over a million SNPs covering the genome, no single SNP shows significant association with ASD or selected phenotypes at a genome-wide level. The SNP that achieves the smallest P-value from secondary analyses is rs1718101. It falls in CNTNAP2, a gene previously implicated in susceptibility for ASD. This SNP also shows modest association with age of word/phrase acquisition in ASD subjects, of interest because features of language development are also associated with other variation in CNTNAP2. In contrast, allele scores derived from the transmission of common alleles to Stage 1 cases significantly predict case status in the independent Stage 2 sample. Despite being significant, the variance explained by these allele scores was small (Vm< 1%). Based on results from individual SNPs and their en masse effect on risk, as inferred from the allele score results, it is reasonable to conclude that common variants affect the risk for ASD but their individual effects are modest.

The expanding role of MBD genes in autism: identification of a MECP2 duplication and novel alterations in MBD5, MBD6, and SETDB1

The methyl-CpG-binding domain (MBD) gene family was first linked to autism over a decade ago when Rett syndrome, which falls under the umbrella of autism spectrum disorders (ASDs), was revealed to be predominantly caused by MECP2 mutations. Since that time, MECP2 alterations have been recognized in idiopathic ASD patients by us and others. Individuals with deletions across the MBD5 gene also present with ASDs, impaired speech, intellectual difficulties, repetitive behaviors, and epilepsy. These findings suggest that further investigations of the MBD gene family may reveal additional associations related to autism. We now describe the first study evaluating individuals with ASD for rare variants in four autosomal MBD family members, MBD5, MBD6, SETDB1, and SETDB2, and expand our initial screening in the MECP2 gene. Each gene was sequenced over all coding exons and evaluated for copy number variations in 287 patients with ASD and an equal number of ethnically matched control individuals. We identified 186 alterations through sequencing, approximately half of which were novel (96 variants, 51.6%). We identified 17 ASD specific, nonsynonymous variants, four of which were concordant in multiplex families: MBD5 Tyr1269Cys, MBD6 Arg883Trp, MECP2 Thr240Ser, and SETDB1 Pro1067del. Furthermore, a complex duplication spanning of the MECP2 gene was identified in two brothers who presented with developmental delay and intellectual disability. From our studies, we provide the first examples of autistic patients carrying potentially detrimental alterations in MBD6 and SETDB1, thereby demonstrating that the MBD gene family potentially plays a significant role in rare and private genetic causes 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.

Evidence of novel fine-scale structural variation at autism spectrum disorder candidate loci

Background

Autism spectrum disorders (ASD) represent a group of neurodevelopmental disorders characterized by a core set of social-communicative and behavioral impairments. Gamma-aminobutyric acid (GABA) is the major inhibitory neurotransmitter in the brain, acting primarily via the GABA receptors (GABR). Multiple lines of evidence, including altered GABA and GABA receptor expression in autistic patients, indicate that the GABAergic system may be involved in the etiology of autism.

Methods

As copy number variations (CNVs), particularly rare and de novo CNVs, have now been implicated in ASD risk, we examined the GABA receptors and genes in related pathways for structural variation that may be associated with autism. We further extended our candidate gene set to include 19 genes and regions that had either been directly implicated in the autism literature or were directly related (via function or ancestry) to these primary candidates. For the high resolution CNV screen we employed custom-designed 244 k comparative genomic hybridization (CGH) arrays. Collectively, our probes spanned a total of 11 Mb of GABA-related and additional candidate regions with a density of approximately one probe every 200 nucleotides, allowing a theoretical resolution for detection of CNVs of approximately 1 kb or greater on average. One hundred and sixty-eight autism cases and 149 control individuals were screened for structural variants. Prioritized CNV events were confirmed using quantitative PCR, and confirmed loci were evaluated on an additional set of 170 cases and 170 control individuals that were not included in the original discovery set. Loci that remained interesting were subsequently screened via quantitative PCR on an additional set of 755 cases and 1,809 unaffected family members.

Results

Results include rare deletions in autistic individuals at JAKMIP1, NRXN1, Neuroligin4Y, OXTR, and ABAT. Common insertion/deletion polymorphisms were detected at several loci, including GABBR2 and NRXN3. Overall, statistically significant enrichment in affected vs. unaffected individuals was observed for NRXN1 deletions.

Conclusions

These results provide additional support for the role of rare structural variation in ASD.

A novel approach of homozygous haplotype sharing identifies candidate genes in autism spectrum disorder

Autism spectrum disorder (ASD) is a highly heritable disorder of complex and heterogeneous aetiology. It is primarily characterized by altered cognitive ability including impaired language and communication skills and fundamental deficits in social reciprocity. Despite some notable successes in neuropsychiatric genetics, overall, the high heritability of ASD (~90%) remains poorly explained by common genetic risk variants. However, recent studies suggest that rare genomic variation, in particular copy number variation, may account for a significant proportion of the genetic basis of ASD. We present a large scale analysis to identify candidate genes which may contain low-frequency recessive variation contributing to ASD while taking into account the potential contribution of population differences to the genetic heterogeneity of ASD. Our strategy, homozygous haplotype (HH) mapping, aims to detect homozygous segments of identical haplotype structure that are shared at a higher frequency amongst ASD patients compared to parental controls. The analysis was performed on 1,402 Autism Genome Project trios genotyped for 1 million single nucleotide polymorphisms (SNPs). We identified 25 known and 1,218 novel ASD candidate genes in the discovery analysis including CADM2, ABHD14A, CHRFAM7A, GRIK2, GRM3, EPHA3, FGF10, KCND2, PDZK1, IMMP2L and FOXP2. Furthermore, 10 of the previously reported ASD genes and 300 of the novel candidates identified in the discovery analysis were replicated in an independent sample of 1,182 trios. Our results demonstrate that regions of HH are significantly enriched for previously reported ASD candidate genes and the observed association is independent of gene size (odds ratio 2.10). Our findings highlight the applicability of HH mapping in complex disorders such as ASD and offer an alternative approach to the analysis of genome-wide association data.

Derivation of autism spectrum disorder-specific induced pluripotent stem cells from peripheral blood mononuclear cells

Induced pluripotent stem cells (iPSCs) hold tremendous potential both as a biological tool to uncover the pathophysiology of disease by creating relevant cell models and as a source of stem cells for cell-based therapeutic applications. Typically, iPSCs have been derived by the transgenic overexpression of transcription factors associated with progenitor cell or stem cell function in fibroblasts derived from skin biopsies. However, the need for skin punch biopsies to derive fibroblasts for reprogramming can present a barrier to study participation among certain populations of individuals, including children with autism spectrum disorders (ASDs). In addition, the acquisition of skin punch biopsies in non-clinic settings presents a challenge. One potential mechanism to avoid these limitations would be the use of peripheral blood mononuclear cells (PBMCs) as the source of the cells for reprogramming. In this article we describe, for the first time, the derivation of iPSC lines from PBMCs isolated from the whole blood of autistic children, and their subsequent differentiation in GABAergic neurons.

An X chromosome-wide association study in autism families identifies TBL1X as a novel autism spectrum disorder candidate gene in males

BACKGROUND

Autism spectrum disorder (ASD) is a complex neurodevelopmental disorder with a strong genetic component. The skewed prevalence toward males and evidence suggestive of linkage to the X chromosome in some studies suggest the presence of X-linked susceptibility genes in people with ASD.

METHODS

We analyzed genome-wide association study (GWAS) data on the X chromosome in three independent autism GWAS data sets: two family data sets and one case-control data set. We performed meta- and joint analyses on the combined family and case-control data sets. In addition to the meta- and joint analyses, we performed replication analysis by using the two family data sets as a discovery data set and the case-control data set as a validation data set.

RESULTS

One SNP, rs17321050, in the transducin β-like 1X-linked (TBL1X) gene [OMIM:300196] showed chromosome-wide significance in the meta-analysis (P value = 4.86 × 10-6) and joint analysis (P value = 4.53 × 10-6) in males. The SNP was also close to the replication threshold of 0.0025 in the discovery data set (P = 5.89 × 10-3) and passed the replication threshold in the validation data set (P = 2.56 × 10-4). Two other SNPs in the same gene in linkage disequilibrium with rs17321050 also showed significance close to the chromosome-wide threshold in the meta-analysis.

CONCLUSIONS

TBL1X is in the Wnt signaling pathway, which has previously been implicated as having a role in autism. Deletions in the Xp22.2 to Xp22.3 region containing TBL1X and surrounding genes are associated with several genetic syndromes that include intellectual disability and autistic features. Our results, based on meta-analysis, joint analysis and replication analysis, suggest that TBL1X may play a role in ASD risk.

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.

Childhood adversity interacts separately with 5-HTTLPR and BDNF to predict lifetime depression diagnosis

The serotonin transporter polymorphism (5-HTTLPR) and the brain-derived neurotrophic factor (BDNF) val66met polymorphism have both been linked to depression symptoms and to depression diagnosis (MDD) in interaction with adversity; there have also been failures to find the effects. We reexamined both interactions for lifetime MDD in a college sample. Lifetime MDD was diagnosed by Structured Clinical Interview for DSM-IV in 133 undergraduates; genotypes for 5-HTTLPR and BDNF were assayed from blood, and self-reports were collected concerning childhood adversity (Risk). 5-HTTLPR interacted with Risk such that Risk predicted less likelihood of MDD among ll carriers and tended to predict greater likelihood of MDD among s carriers. BDNF interacted with Risk such that Risk predicted greater likelihood of MDD among met carriers and did not influence val/val carriers. These two interactions were additive: both were significant in a combined model.

Microduplications in an autism multiplex family narrow the region of susceptibility for developmental disorders on 15q24 and implicate 7p21

Copy number variations (CNVs) play a crucial role in the intricate genetics of autism spectrum disorders. A region on chromosome 15q24 vulnerable to both deletions and duplications has been previously implicated in a range of phenotypes including autism, Asperger's syndrome, delayed development, and mild to severe mental retardation. Prior studies have delineated a minimal critical region of approximately 1.33 Mb. In this study, a multiplex autism family was evaluated for CNVs using genotyping data from the Illumina 1 M BeadChip and analyzed with the PennCNV algorithm. Variants were then identified that co-segregate with autism features in this family. Here, we report autistic first cousins who carry two microduplications concordant with disease. Both duplications were inherited maternally and found to be identical by descent. The first is an approximately 10,000 base pair microduplication within the minimal region on 15q24 that falls across a single gene, ubiquitin-like 7. This is the smallest duplication in the region to result in a neuropsychiatric disorder, potentially narrowing the critical region for susceptibility to developmental and autism spectrum disorders. The second is a novel, 352 kb tandem duplication on 7p21 that replicates part of the neurexophilin 1 and islet cell autoantigen 1 genes. The breakpoint junction falls within the intronic regions of these genes and demonstrates a microhomology of four base pairs. Each of these microduplications may contribute to the complex etiology of autism spectrum disorders.

Variants in several genomic regions associated with asperger disorder

Asperger disorder (ASP) is one of the autism spectrum disorders (ASD) and is differentiated from autism largely on the absence of clinically significant cognitive and language delays. Analysis of a homogenous subset of families with ASP may help to address the corresponding effect of genetic heterogeneity on identifying ASD genetic risk factors. To examine the hypothesis that common variation is important in ASD, we performed a genome-wide association study (GWAS) in 124 ASP families in a discovery data set and 110 ASP families in a validation data set. We prioritized the top 100 association results from both cohorts by employing a ranking strategy. Novel regions on 5q21.1 (P = 9.7 × 10(-7) ) and 15q22.1-q22.2 (P = 7.3 × 10(-6) ) were our most significant findings in the combined data set. Three chromosomal regions showing association, 3p14.2 (P = 3.6 × 10(-6) ), 3q25-26 (P = 6.0 × 10(-5) ) and 3p23 (P = 3.3 × 10(-4) ) overlapped linkage regions reported in Finnish ASP families, and eight association regions overlapped ASD linkage areas. Our findings suggest that ASP shares both ASD-related genetic risk factors, as well as has genetic risk factors unique to the ASP phenotype.

A de novo 1.5 Mb microdeletion on chromosome 14q23.2-23.3 in a patient with autism and spherocytosis

Autism is a neuro-developmental disorder characterized by deficits in social interaction and communication as well as restricted interests or repetitive behaviors. Cytogenetic studies have implicated large chromosomal aberrations in the etiology of approximately 5-7% of autism patients, and the recent advent of array-based techniques allows the exploration of submicroscopic copy number variations (CNVs). We genotyped a 14-year-old boy with autism, spherocytosis and other physical dysmorphia, his parents, and two non-autistic siblings with the Illumina Human 1M Beadchip as part of a study of the molecular genetics of autism and determined copy number variants using the PennCNV algorithm. We identified and validated a de novo 1.5 Mb microdeletion of 14q23.2-23.3 in our autistic patient. This region contains 15 genes, including spectrin beta (SPTB), encoding a cytoskeletal protein previously associated with spherocytosis, methylenetetrahydrofolate dehydrogenase 1 (MTHFD1), a folate metabolizing enzyme previously associated with bipoloar disorder and schizophrenia, pleckstrin homology domain-containing family G member 3 (PLEKHG3), a guanide nucleotide exchange enriched in the brain, and churchill domain containing protein 1 (CHURC1), homologs of which regulate neuronal development in model organisms. While a similar deletion has previously been reported in a family with spherocytosis, severe learning disabilities, and mild mental retardation, this is the first implication of chr14q23.2-23.3 in the etiology of autism and points to MTHFD1, PLEKHG3, and CHURC1 as potential candidate genes contributing to autism risk.

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

BACKGROUND

Genome-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.

METHODS

GWAS-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 Fisher's 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.

RESULTS

Computer simulations indicate that GWAS-NR achieves a significantly higher classification rate for true positive association signals than either the joint analysis or Fisher's 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.

CONCLUSIONS

As 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.

A genome-wide scan for common alleles affecting risk for autism

Although autism spectrum disorders (ASDs) have a substantial genetic basis, most of the known genetic risk has been traced to rare variants, principally copy number variants (CNVs). To identify common risk variation, the Autism Genome Project (AGP) Consortium genotyped 1558 rigorously defined ASD families for 1 million single-nucleotide polymorphisms (SNPs) and analyzed these SNP genotypes for association with ASD. In one of four primary association analyses, the association signal for marker rs4141463, located within MACROD2, crossed the genome-wide association significance threshold of P < 5 × 10⁻⁸. When a smaller replication sample was analyzed, the risk allele at rs4141463 was again over-transmitted; yet, consistent with the winner's curse, its effect size in the replication sample was much smaller; and, for the combined samples, the association signal barely fell below the P < 5 × 10⁻⁸ threshold. Exploratory analyses of phenotypic subtypes yielded no significant associations after correction for multiple testing. They did, however, yield strong signals within several genes, KIAA0564, PLD5, POU6F2, ST8SIA2 and TAF1C.

Novel variants identified in methyl-CpG-binding domain genes in autistic individuals

Misregulation of the methyl-CpG-binding protein 2 (MECP2) gene has been found to cause a myriad of neurological disorders including autism, mental retardation, seizures, learning disabilities, and Rett syndrome. We hypothesized that mutations in other members of the methyl-CpG-binding domain (MBD) family may also cause autistic features in individuals. We evaluated 226 autistic individuals for alterations in the four genes most homologous to MECP2: MBD1, MBD2, MBD3, and MBD4. A total of 46 alterations were identified in the four genes, including ten missense changes and two deletions that alter coding sequence. Several are either unique to our autistic population or cosegregate with affected individuals within a family, suggesting a possible relation of these variations to disease etiology. Variants include a R23M alteration in two affected half brothers which falls within the MBD domain of the MBD3 protein, as well as a frameshift in MBD4 that is predicted to truncate almost half of the protein. These results suggest that rare cases of autism may be influenced by mutations in members of the dynamic MBD protein family.

Functional impact of global rare copy number variation in autism spectrum disorders

The autism spectrum disorders (ASDs) are a group of conditions characterized by impairments in reciprocal social interaction and communication, and the presence of restricted and repetitive behaviours. Individuals with an ASD vary greatly in cognitive development, which can range from above average to intellectual disability. Although ASDs are known to be highly heritable ( approximately 90%), the underlying genetic determinants are still largely unknown. Here we analysed the genome-wide characteristics of rare (<1% frequency) copy number variation in ASD using dense genotyping arrays. When comparing 996 ASD individuals of European ancestry to 1,287 matched controls, cases were found to carry a higher global burden of rare, genic copy number variants (CNVs) (1.19 fold, P = 0.012), especially so for loci previously implicated in either ASD and/or intellectual disability (1.69 fold, P = 3.4 x 10(-4)). Among the CNVs there were numerous de novo and inherited events, sometimes in combination in a given family, implicating many novel ASD genes such as SHANK2, SYNGAP1, DLGAP2 and the X-linked DDX53-PTCHD1 locus. We also discovered an enrichment of CNVs disrupting functional gene sets involved in cellular proliferation, projection and motility, and GTPase/Ras signalling. Our results reveal many new genetic and functional targets in ASD that may lead to final connected pathways.

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.

Identification of chromosome 7 inversion breakpoints in an autistic family narrows candidate region for autism susceptibility

Chromosomal breaks and rearrangements have been observed in conjunction with autism and autistic spectrum disorders. A chromosomal inversion has been previously reported in autistic siblings, spanning the region from approximately 7q22.1 to 7q31. This family is distinguished by having multiple individuals with autism and associated disabilities. The region containing the inversion has been strongly implicated in autism by multiple linkage studies, and has been particularly associated with language defects in autism as well as in other disorders with language components. Mapping of the inversion breakpoints by FISH has localized the inversion to the region spanning approximately 99-108.75 Mb of chromosome 7. The proximal breakpoint has the potential to disrupt either the coding sequence or regulatory regions of a number of cytochrome P450 genes while the distal region falls in a relative gene desert. Copy number variant analysis of the breakpoint regions detected no duplication or deletion that could clearly be associated with disease status. Association analysis in our autism data set using single nucleotide polymorphisms located near the breakpoints showed no significant association with proximal breakpoint markers, but has identified markers near the distal breakpoint ( approximately 108-110 Mb) with significant associations to autism. The chromosomal abnormality in this family strengthens the case for an autism susceptibility gene in the chromosome 7q22-31 region and targets a candidate region for further investigation.

Examination of tetrahydrobiopterin pathway genes in autism

Autism is a complex disorder with a high degree of heritability and significant phenotypic and genotypic heterogeneity. Although candidate gene studies and genome-wide screens have failed to identify major causal loci associated with autism, numerous studies have proposed association with several variations in genes in the dopaminergic and serotonergic pathways. Because tetrahydrobiopterin (BH4) is the essential cofactor in the synthesis of these two neurotransmitters, we genotyped 25 SNPs in nine genes of the BH4 pathway in a total of 403 families. Significant nominal association was detected in the gene for 6-pyruvoyl-tetrahydropterin synthase, PTS (chromosome 11), with P = 0.009; this result was not restricted to an affected male-only subset. Multilocus interaction was detected in the BH4 pathway alone, but not across the serotonin, dopamine and BH4 pathways.

Genomic and epigenetic evidence for oxytocin receptor deficiency in autism

BACKGROUND

Autism comprises a spectrum of behavioral and cognitive disturbances of childhood development and is known to be highly heritable. Although numerous approaches have been used to identify genes implicated in the development of autism, less than 10% of autism cases have been attributed to single gene disorders.

METHODS

We describe the use of high-resolution genome-wide tilepath microarrays and comparative genomic hybridization to identify copy number variants within 119 probands from multiplex autism families. We next carried out DNA methylation analysis by bisulfite sequencing in a proband and his family, expanding this analysis to methylation analysis of peripheral blood and temporal cortex DNA of autism cases and matched controls from independent datasets. We also assessed oxytocin receptor (OXTR) gene expression within the temporal cortex tissue by quantitative real-time polymerase chain reaction (PCR).

RESULTS

Our analysis revealed a genomic deletion containing the oxytocin receptor gene, OXTR (MIM accession no.: 167055), previously implicated in autism, was present in an autism proband and his mother who exhibits symptoms of obsessive-compulsive disorder. The proband's affected sibling did not harbor this deletion but instead may exhibit epigenetic misregulation of this gene through aberrant gene silencing by DNA methylation. Further DNA methylation analysis of the CpG island known to regulate OXTR expression identified several CpG dinucleotides that show independent statistically significant increases in the DNA methylation status in the peripheral blood cells and temporal cortex in independent datasets of individuals with autism as compared to control samples. Associated with the increase in methylation of these CpG dinucleotides is our finding that OXTR mRNA showed decreased expression in the temporal cortex tissue of autism cases matched for age and sex compared to controls.

CONCLUSION

Together, these data provide further evidence for the role of OXTR and the oxytocin signaling pathway in the etiology of autism and, for the first time, implicate the epigenetic regulation of OXTR in the development of the disorder.See the related commentary by Gurrieri and Neri: http://www.biomedcentral.com/1741-7015/7/63.

Examination of association to autism of common genetic variationin genes related to dopamine

Autism is a severe neurodevelopmental disorder characterized by a triad of complications. Autistic individuals display significant disturbances in language and reciprocal social interactions, combined with repetitive and stereotypic behaviors. Prevalence studies suggest that autism is more common than originally believed, with recent estimates citing a rate of one in 150. Although multiple genetic linkage and association studies have yielded multiple suggestive genes or chromosomal regions, a specific risk locus has yet to be identified and widely confirmed. Because many etiologies have been suggested for this complex syndrome, we hypothesize that one of the difficulties in identifying autism genes is that multiple genetic variants may be required to significantly increase the risk of developing autism. Thus, we took the alternative approach of examining 14 prominent dopamine pathway candidate genes for detailed study by genotyping 28 single nucleotide polymorphisms. Although we did observe a nominally significant association for rs2239535 (P=0.008) on chromosome 20, single-locus analysis did not reveal any results as significant after correction for multiple comparisons. No significant interaction was identified when Multifactor Dimensionality Reduction was employed to test specifically for multilocus effects. Although genome-wide linkage scans in autism have provided support for linkage to various loci along the dopamine pathway, our study does not provide strong evidence of linkage or association to any specific gene or combination of genes within the pathway. These results demonstrate that common genetic variation within the tested genes located within this pathway at most play a minor to moderate role in overall autism pathogenesis.

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.