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.

Targeted sequencing of ABCA7 identifies splicing, stop-gain and intronic risk variants for Alzheimer disease.

Several variants in the gene ABCA7 have been identified as potential causal variants for late-onset Alzheimer's disease (LOAD). In order to replicate these findings, and search for novel causal variants, we performed targeted sequencing of this gene in cohorts of non-Hispanic White (NHW) and African-American (AA) LOAD cases and controls. We sequenced the gene ABCA7 in 291 NHW LOAD cases and 103 controls. Variants were prioritized for rare, damaging variants and previously reported variants associated with LOAD, and were follow-up genotyped in 4076 NHW and 1157 AA cases and controls. We confirm three previously associated ABCA7 risk variants and extend two of these associations to other populations, an intronic variant in NHW (P=3.0×10-3) (originally reported in a Belgian population), and a splice variant originally associated in the Icelandic population, which was significantly associated in the NHW cohort (P=1.2×10-6) and nominally associated in the AA cohort (P=0.017). We also identify a 3'-UTR splice variant that segregates in four siblings of one family and is nominally associated with LOAD (P=0.040). Multiple variants in ABCA7 contribute to LOAD risk.

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.

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.

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.

Examination of association of genes in the serotonin system to autism.

Autism is characterized as one of the pervasive developmental disorders, a spectrum of often severe behavioral and cognitive disturbances of early development. The high heritability of autism has driven multiple efforts to identify genetic variation that increases autism susceptibility. Numerous studies have suggested that variation in peripheral and central metabolism of serotonin (5-hydroxytryptamine) may play a role in the pathophysiology of autism. We screened 403 autism families for 45 single nucleotide polymorphisms in ten serotonin pathway candidate genes. Although genome-wide linkage scans in autism have provided support for linkage to various loci located within the serotonin pathway, our study does not provide strong evidence for linkage to any specific gene within the pathway. The most significant association (p = 0.0002; p = 0.02 after correcting for multiple comparisons) was found at rs1150220 (HTR3A) located on chromosome 11 ( approximately 113 Mb). To test specifically for multilocus effects, multifactor dimensionality reduction was employed, and a significant two-way interaction (p value = 0.01) was found between rs10830962, near MTNR1B (chromosome11; 92,338,075 bp), and rs1007631, near SLC7A5 (chromosome16; 86,413,596 bp). These data suggest that variation within genes on the serotonin pathway, particularly HTR3A, may have modest effects on autism risk.

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.

Multifactor dimensionality reduction-phenomics: a novel method to capture genetic heterogeneity with use of phenotypic variables.

Complex human diseases do not have a clear inheritance pattern, and it is expected that risk involves multiple genes with modest effects acting independently or interacting. Major challenges for the identification of genetic effects are genetic heterogeneity and difficulty in analyzing high-order interactions. To address these challenges, we present MDR-Phenomics, a novel approach based on the multifactor dimensionality reduction (MDR) method, to detect genetic effects in pedigree data by integration of phenotypic covariates (PCs) that may reflect genetic heterogeneity. The P value of the test is calculated using a permutation test adjusted for multiple tests. To validate MDR-Phenomics, we compared it with two MDR-based methods: (1) traditional MDR pedigree disequilibrium test (PDT) without consideration of PCs (MDR-PDT) and (2) stratified phenotype (SP) analysis based on PCs, with use of MDR-PDT with a Bonferroni adjustment (SP-MDR). Using computer simulations, we examined the statistical power and type I error of the different approaches under several genetic models and sampling scenarios. We conclude that MDR-Phenomics is more powerful than MDR-PDT and SP-MDR when there is genetic heterogeneity, and the statistical power is affected by sample size and the number of PC levels. We further compared MDR-Phenomics with conditional logistic regression (CLR) for testing interactions across single or multiple loci with consideration of PC. The results show that CLR with PC has only slightly smaller power than does MDR-Phenomics for single-locus analysis but has considerably smaller power for multiple loci. Finally, by applying MDR-Phenomics to autism, a complex disease in which multiple genes are believed to confer risk, we attempted to identify multiple gene effects in two candidate genes of interest--the serotonin transporter gene (SLC6A4) and the integrin beta 3 gene (ITGB3) on chromosome 17. Analyzing four markers in SLC6A4 and four markers in ITGB3 in 117 white family triads with autism and using sex of the proband as a PC, we found significant interaction between two markers--rs1042173 in SLC6A4 and rs3809865 in ITGB3.

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.

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.

The relationship between restrictive and repetitive behaviors in individuals with autism and obsessive compulsive symptoms in parents.

This study investigated the relationship between repetitive behaviors in individuals with autism and obsessive-compulsive behaviors in parents. We hypothesized that repetitive behaviors in probands with autism would be associated with increased obsessive-compulsive behaviors in parents in sporadic families (1 known case of autism per family and no known history of autism). Parents with clinically significant Y-BOCS scores were more likely to have a family history of obsessive-compulsive disorder. The empirically derived Autism Diagnostic Interview-R (ADI-R) factor, Insistence on Sameness, was positively correlated with obsessive-compulsive behaviors in parents. Further, when probands were grouped on the basis of parental Y-BOCS scores (clinically significant versus non-clinically significant), probands whose parents had clinically significant Y-BOCS scores had higher ADI-R Insistence on Sameness factor scores. The findings of the current study of sporadic families extend previous work that has shown an association between restrictive/repetitive behaviors in probands with autism and obsessive-compulsive features in parents.

Identification of significant association and gene-gene interaction of GABA receptor subunit genes in autism.

Autism is a common neurodevelopmental disorder with a significant genetic component. Existing research suggests that multiple genes contribute to autism and that epigenetic effects or gene-gene interactions are likely contributors to autism risk. However, these effects have not yet been identified. Gamma-aminobutyric acid (GABA), the primary inhibitory neurotransmitter in the adult brain, has been implicated in autism etiology. Fourteen known autosomal GABA receptor subunit genes were studied to look for the genes associated with autism and their possible interactions. Single-nucleotide polymorphisms (SNPs) were screened in the following genes: GABRG1, GABRA2, GABRA4, and GABRB1 on chromosome 4p12; GABRB2, GABRA6, GABRA1, GABRG2, and GABRP on 5q34-q35.1; GABRR1 and GABRR2 on 6q15; and GABRA5, GABRB3, and GABRG3 on 15q12. Intronic and/or silent mutation SNPs within each gene were analyzed in 470 white families with autism. Initially, SNPs were used in a family-based study for allelic association analysis--with the pedigree disequilibrium test and the family-based association test--and for genotypic and haplotypic association analysis--with the genotype-pedigree disequilibrium test (geno-PDT), the association in the presence of linkage (APL) test, and the haplotype family-based association test. Next, with the use of five refined independent marker sets, extended multifactor-dimensionality reduction (EMDR) analysis was employed to identify the models with locus joint effects, and interaction was further verified by conditional logistic regression. Significant allelic association was found for markers RS1912960 (in GABRA4; P = .01) and HCV9866022 (in GABRR2; P = .04). The geno-PDT found significant genotypic association for HCV8262334 (in GABRA2), RS1912960 and RS2280073 (in GABRA4), and RS2617503 and RS12187676 (in GABRB2). Consistent with the allelic and genotypic association results, EMDR confirmed the main effect at RS1912960 (in GABRA4). EMDR also identified a significant two-locus gene-gene effect model involving RS1912960 in GABRA4 and RS2351299 in GABRB1. Further support for this two-locus model came from both the multilocus geno-PDT and the APL test, which indicated a common genotype and haplotype combination positively associated with disease. Finally, these results were also consistent with the results from the conditional logistic regression, which confirmed the interaction between GABRA4 and GABRB1 (odds ratio = 2.9 for interaction term; P = .002). Through the convergence of all analyses, we conclude that GABRA4 is involved in the etiology of autism and potentially increases autism risk through interaction with GABRB1. These results support the hypothesis that GABA receptor subunit genes are involved in autism, most likely via complex gene-gene interactions.

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.

No association between the APOE gene and autism.

Autism is a neurodevelopmental disorder characterized by stereotypic and repetitive behavior and interests, together with social and communicative deficiencies. The results of several genomic screens suggest the presence of an autism susceptibility locus on chromosome 19p13.2-q13.4. The apolipoprotein E (APOE) gene on chromosome 19 encodes for a protein, apoE, whose different isoforms (E2, E3, E4) influence neuronal growth. APOE participates in lipid transport and metabolism, repair, growth, and maintenance of axons and myelin during neuronal development. The APOE protein competes with the Reelin protein for VLDL/APOER2 receptor binding. Several studies have reported evidence for an association between autism and the Reelin gene. Based on these data we tested for association between APOE and autism using family-based association methods in a data set of 322 autism families. Three promoter, one intronic, and one 3' UTR single nucleotide polymorphisms (SNPs) in the APOE gene (-491a/t, -427c/t, -219g/t, 113c/g, and 5361c/t) as well as the APOE functional polymorphism (E2, E3, E4) were examined and failed to reveal significant evidence that autism is associated with APOE.

Fine mapping of autistic disorder to chromosome 15q11-q13 by use of phenotypic subtypes.

Autistic disorder (AutD) is a complex genetic disease. Available evidence suggests that several genes contribute to the underlying genetic risk for the development of AutD. However, both etiologic heterogeneity and genetic heterogeneity confound the discovery of AutD-susceptibility genes. Chromosome 15q11-q13 has been identified as a strong candidate region on the basis of both the frequent occurrence of chromosomal abnormalities in that region and numerous suggestive linkage and association findings. Ordered-subset analysis (OSA) is a novel statistical method to identify a homogeneous subset of families that contribute to overall linkage at a given chromosomal location and thus to potentially help in the fine mapping and localization of the susceptibility gene within a chromosomal area. For the present analysis, a factor that represents insistence on sameness (IS)--derived from a principal-component factor analysis using data on 221 patients with AutD from the repetitive behaviors/stereotyped patterns domain in the Autism Diagnostic Interview-Revised--was used as a covariate in OSA. Analysis of families sharing high scores on the IS factor increased linkage evidence for the 15q11-q13 region, at the GABRB3 locus, from a LOD score of 1.45 to a LOD score of 4.71. These results narrow our region of interest on chromosome 15 to an area surrounding the gamma-aminobutyric acid-receptor subunit genes, in AutD, and support the hypothesis that the analysis of phenotypic homogeneous subtypes may be a powerful tool for the mapping of disease-susceptibility genes in complex traits.

De novo partial duplication of chromosome 7p in a male with autistic disorder.

We describe a de novo partial duplication of 7p in a 25-year-old male with autistic disorder (AD). High-resolution chromosome analysis revealed an extra segment added to the proximal short arm of chromosome 7. The G-band pattern was consistent with an inverted duplication of 7p11.2-p14.1. Fluorescent in situ hybridization (FISH), using a whole chromosome 7 DNA probe (Cytocell, Inc., UK), confirmed that the extra chromosome material is derived from chromosome 7, indicating that the patient is partially trisomic for a region of the short arm of chromosome 7. Partial duplication of the short arm of chromosome 7 is uncommon with little more than 30 cases in the literature. This is the first report of an individual with a 7p duplication who also has AD.

Association analysis of chromosome 15 gabaa receptor subunit genes in autistic disorder.

Gamma-aminobutyric acid (GABA) is the major inhibitory neurotransmitter in the brain, acting via the GABAA receptors. The GABAA receptors are comprised of several different homologous subunits, forming a group of receptors that are both structurally and functionally diverse. Three of the GABAA receptor subunit genes (GABRB3, GABRA5 and GABRG3) form a cluster on chromosome 15q11-q13, in a region that has been genetically associated with autistic disorder (AutD). Based on these data, we examined 16 single nucleotide polymorphisms (SNPs) located within GABRB3, GABRA5 and GABRG3 for linkage disequilibrium (LD) in 226 AutD families (AutD patients and parents). Genotyping was performed using either OLA (oligonucleotide ligation assay), or SSCP (single strand conformation polymorphism) followed by DNA sequencing. We tested for LD using the Pedigree Disequilibrium Test (PDT). PDT results gave significant evidence that AutD is associated with two SNPs located within the GABRG3 gene (exon5_539T/C, p=0.02 and intron5_687T/C, p=0.03), suggesting that the GABRG3 gene or a gene nearby contributes to genetic risk in AutD.

Genetic studies in autistic disorder and chromosome 15.

Autistic disorder (AD) is a developmental disorder affecting social interactions, communication, and behavior. AD is a disease of complex genetic architecture. It is postulated that several genes contribute to the underlying etiology of AD. Chromosome 15 is of particular interest due to numerous reports of AD in the presence of chromosomal abnormalities, located mainly in the 15q11-q13 region. There are also a number of plausible candidate genes in this area, including the gamma-aminobutyric acidA (GABA(A)) receptor gene complex. We have undertaken a study of this region of chromosome 15 in a data set of 63 multiplex families (with 2 or more AD affected individuals per family). We found evidence in support of linkage to the 15q11-q13 region, as well as evidence of increased recombination in this region. These findings provide further support for the involvement of chromosome 15q11-q13 in the genetic etiology of AD.