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.     

Bipolar I disorder and schizophrenia: a 440-single-nucleotide polymorphism screen of 64 candidate genes among Ashkenazi Jewish case-parent trios

Bipolar, schizophrenia, and schizoaffective disorders are common, highly heritable psychiatric disorders, for which familial coaggregation, as well as epidemiological and genetic evidence, suggests overlapping etiologies. No definitive susceptibility genes have yet been identified for any of these disorders. Genetic heterogeneity, combined with phenotypic imprecision and poor marker coverage, has contributed to the difficulty in defining risk variants. We focused on families of Ashkenazi Jewish descent, to reduce genetic heterogeneity, and, as a precursor to genomewide association studies, we undertook a single-nucleotide polymorphism (SNP) genotyping screen of 64 candidate genes (440 SNPs) chosen on the basis of previous linkage or of association and/or biological relevance. We genotyped an average of 6.9 SNPs per gene, with an average density of 1 SNP per 11.9 kb in 323 bipolar I disorder and 274 schizophrenia or schizoaffective Ashkenazi case-parent trios. Using single-SNP and haplotype-based transmission/disequilibrium tests, we ranked genes on the basis of strength of association (P<.01). Six genes (DAO, GRM3, GRM4, GRIN2B, IL2RB, and TUBA8) met this criterion for bipolar I disorder; only DAO has been previously associated with bipolar disorder. Six genes (RGS4, SCA1, GRM4, DPYSL2, NOS1, and GRID1) met this criterion for schizophrenia or schizoaffective disorder; five replicate previous associations, and one, GRID1, shows a novel association with schizophrenia. In addition, six genes (DPYSL2, DTNBP1, G30/G72, GRID1, GRM4, and NOS1) showed overlapping suggestive evidence of association in both disorders. These results may help to prioritize candidate genes for future study from among the many suspected/proposed for schizophrenia and bipolar disorders. They provide further support for shared genetic susceptibility between these two disorders that involve glutamate-signaling pathways.     

Functional polymorphisms of HSPA5: possible association with bipolar disorder

Altered endoplasmic reticulum stress (ER) response signaling is suggested in bipolar disorder. Previously, we preliminarily reported the genetic association of HSPA5 (GRP78/BiP) with bipolar disorder. Here, we extended our analysis by increasing the number of Japanese case-control samples and NIMH Genetics Initiative bipolar trio samples (NIMH trios), and also analyzed schizophrenia samples. In Japanese, nominally significant association of one haplotype was observed in extended samples of bipolar disorder but not in schizophrenia. In NIMH trios, no association was found in total samples. However, an exploratory analysis suggested that the other haplotype was significantly over-transmitted to probands only from the paternal side. The associated haplotype in Japanese or NIMH pedigrees shared three common polymorphisms in the promotor, which was found to alter promotor activity. These findings suggested promotor polymorphisms of HSPA5 may affect the interindividual variability of ER stress response and may confer a genetic risk factor for bipolar disorder.     

A Cytogenetic Abnormality and Rare Coding Variants Identify ABCA13 as a Candidate Gene in Schizophrenia, Bipolar Disorder, and Depression

Schizophrenia and bipolar disorder are leading causes of morbidity across all populations, with heritability estimates of ∼80% indicating a substantial genetic component. Population genetics and genome-wide association studies suggest an overlap of genetic risk factors between these illnesses but it is unclear how this genetic component is divided between common gene polymorphisms, rare genomic copy number variants, and rare gene sequence mutations. We report evidence that the lipid transporter gene ABCA13 is a susceptibility factor for both schizophrenia and bipolar disorder. After the initial discovery of its disruption by a chromosome abnormality in a person with schizophrenia, we resequenced ABCA13 exons in 100 cases with schizophrenia and 100 controls. Multiple rare coding variants were identified including one nonsense and nine missense mutations and compound heterozygosity/homozygosity in six cases. Variants were genotyped in additional schizophrenia, bipolar, depression (n > 1600), and control (n > 950) cohorts and the frequency of all rare variants combined was greater than controls in schizophrenia (OR = 1.93, p = 0.0057) and bipolar disorder (OR = 2.71, p = 0.00007). The population attributable risk of these mutations was 2.2% for schizophrenia and 4.0% for bipolar disorder. In a study of 21 families of mutation carriers, we genotyped affected and unaffected relatives and found significant linkage (LOD = 4.3) of rare variants with a phenotype including schizophrenia, bipolar disorder, and major depression. These data identify a candidate gene, highlight the genetic overlap between schizophrenia, bipolar disorder, and depression, and suggest that rare coding variants may contribute significantly to risk of these disorders.     

Genome-wide searching of rare genetic variants in WTCCC data

Although they have demonstrated success in searching for common variants for complex diseases, genome-wide association (GWA) studies are less successful in detecting rare genetic variants because of the poor statistical power of most of current methods. We developed a two-stage method that can apply to GWA studies for detecting rare variants. Here we report the results of applying this two-stage method to the Wellcome Trust Case Control Consortium (WTCCC) dataset that include seven complex diseases: bipolar disorder, cardiovascular disease, hypertension (HT), rheumatoid arthritis, Crohn’s disease, type 1 diabetes and type 2 diabetes (T2D). We identified 24 genes or regions that reach genome wide significance. Eight of them are novel and were not reported in the WTCCC study. The cumulative risk (or protective) haplotype frequency for each of the 8 genes or regions is small, being at most 11%. For each of the novel genes, the risk (or protective) haplotype set cannot be tagged by the common SNPs available in chips (r ² < 0.32). The gene identified in HT was further replicated in the Framingham Heart Study, and is also significantly associated with T2D. Our analysis suggests that searching for rare genetic variants is feasible in current GWA studies and candidate gene studies, and the results can severe as guides to future resequencing studies to identify the underlying rare functional variants.     

The impact of using related individuals for haplotype reconstruction in population studies

Recent studies have highlighted the dangers of using haplotypes reconstructed directly from population data for a fine-scale mapping analysis. Family data may help resolve ambiguity, yet can be costly to obtain. This study is concerned with the following question: How much family data (if any) should be used to facilitate haplotype reconstruction in a population study? We conduct a simulation study to evaluate how changes in family information can affect the accuracy of haplotype frequency estimates and phase reconstruction. To reconstruct haplotypes, we introduce an EM-based algorithm that can efficiently accommodate unrelated individuals, parent-child trios, and arbitrarily large half-sib pedigrees. Simulations are conducted for a diverse set of haplotype frequency distributions, all of which have been previously published in empirical studies. A wide variety of important results regarding the effectiveness of using pedigree data in a population study are presented in a coherent, unified framework. Insight into the different properties of the haplotype frequency distribution that can influence experimental design is provided. We show that a preliminary estimate of the haplotype frequency distribution can be valuable in large population studies with fixed resources.     

Candidate-gene screening and association analysis at the autism-susceptibility locus on chromosome 16p: evidence of association at GRIN2A and ABAT

Autism is a highly heritable neurodevelopmental disorder whose underlying genetic causes have yet to be identified. To date, there have been eight genome screens for autism, two of which identified a putative susceptibility locus on chromosome 16p. In the present study, 10 positional candidate genes that map to 16p11-13 were examined for coding variants: A2BP1, ABAT, BFAR, CREBBP, EMP2, GRIN2A, MRTF-B, SSTR5, TBX6, and UBN1. Screening of all coding and regulatory regions by denaturing high-performance liquid chromatography identified seven nonsynonymous changes. Five of these mutations were found to cosegregate with autism, but the mutations are not predicted to have deleterious effects on protein structure and are unlikely to represent significant etiological variants. Selected variants from candidate genes were genotyped in the entire International Molecular Genetics Study of Autism Consortium collection of 239 multiplex families and were tested for association with autism by use of the pedigree disequilibrium test. Additionally, genotype frequencies were compared between 239 unrelated affected individuals and 192 controls. Patterns of linkage disequilibrium were investigated, and the transmission of haplotypes across candidate genes was tested for association. Evidence of single-marker association was found for variants in ABAT, CREBBP, and GRIN2A. Within these genes, 12 single-nucleotide polymorphisms (SNPs) were subsequently genotyped in 91 autism trios (one affected individual and two unaffected parents), and the association was replicated within GRIN2A (Fisher's exact test, P<.0001). Logistic regression analysis of SNP data across GRIN2A and ABAT showed a trend toward haplotypic differences between cases and controls.     

HAPLORE: a program for haplotype reconstruction in general pedigrees without recombination

MOTIVATION: Haplotype reconstruction is an essential step in genetic linkage and association studies. Although many methods have been developed to estimate haplotype frequencies and reconstruct haplotypes for a sample of unrelated individuals, haplotype reconstruction in large pedigrees with a large number of genetic markers remains a challenging problem. METHODS: We have developed an efficient computer program, HAPLORE (HAPLOtype REconstruction), to identify all haplotype sets that are compatible with the observed genotypes in a pedigree for tightly linked genetic markers. HAPLORE consists of three steps that can serve different needs in applications. In the first step, a set of logic rules is used to reduce the number of compatible haplotypes of each individual in the pedigree as much as possible. After this step, the haplotypes of all individuals in the pedigree can be completely or partially determined. These logic rules are applicable to completely linked markers and they can be used to impute missing data and check genotyping errors. In the second step, a haplotype-elimination algorithm similar to the genotype-elimination algorithms used in linkage analysis is applied to delete incompatible haplotypes derived from the first step. All superfluous haplotypes of the pedigree members will be excluded after this step. In the third step, the expectation-maximization (EM) algorithm combined with the partition and ligation technique is used to estimate haplotype frequencies based on the inferred haplotype configurations through the first two steps. Only compatible haplotype configurations with haplotypes having frequencies greater than a threshold are retained. RESULTS: We test the effectiveness and the efficiency of HAPLORE using both simulated and real datasets. Our results show that, the rule-based algorithm is very efficient for completely genotyped pedigree. In this case, almost all of the families have one unique haplotype configuration. In the presence of missing data, the number of compatible haplotypes can be substantially reduced by HAPLORE, and the program will provide all possible haplotype configurations of a pedigree under different circumstances, if such multiple configurations exist. These inferred haplotype configurations, as well as the haplotype frequencies estimated by the EM algorithm, can be used in genetic linkage and association studies. AVAILABILITY: The program can be downloaded from http://bioinformatics.med.yale.edu.     

Evolutionary signatures of common human cis-regulatory haplotypes

Variation in gene expression may give rise to a significant fraction of inter-individual phenotypic variation. Studies searching for the underlying genetic controls for such variation have been conducted in model organisms and humans in recent years. In our previous effort of assessing conserved underlying haplotype patterns across ethnic populations, we constructed common haplotypes using SNPs having conserved linkage disequilibrium (LD) across ethnic populations. These common haplotypes cluster into a simple evolutionary structure based on their frequencies, defining only up to three conserved clusters termed 'haplotype frameworks'. One intriguing preliminary finding was that a significant portion of reported variants strongly associated with cis-regulation tags these globally conserved haplotype frameworks. Here we expand the investigation by collecting genes showing stringently determined cis-association between genotypes and expression phenotypes from major studies. We conducted phylogenetic analysis of current major haplotypes along with the corresponding haplotypes derived from chimpanzee reference sequences. Our analysis reveals that, for the vast majority of such cis-regulatory genes, the tagging SNPs showing the strongest association also tag the haplotype lineages directly separated from ancestry, inferred from either chimpanzee reference sequences or the allele frequency-derived haplotype frameworks, suggesting that the differentially expressed phenotypes were evolved relatively early in human history. Such evolutionary signatures provide keys for a more effective identification of globally-conserved candidate regulatory haplotypes across human genes in future epidemiologic and pharmacogenetic studies.     

Selection and evaluation of tagging SNPs in the neuronal-sodium-channel gene SCN1A: implications for linkage-disequilibrium gene mapping

Association studies are widely seen as the most promising approach for finding polymorphisms that influence genetically complex traits, such as common diseases and responses to their treatment. Considerable interest has therefore recently focused on the development of methods that efficiently screen genomic regions or whole genomes for gene variants associated with complex phenotypes. One key element in this search is the use of linkage disequilibrium to gain maximal information from typing a selected subset of highly informative single-nucleotide polymorphism (SNP) markers, now often called "tagging SNPs" (tSNPs). Probably the most common approach to linkage-disequilibrium gene mapping involves a three-step program: (1) characterization of the haplotype structure in candidate genes or genomic regions of interest, (2) identification of tSNPs sufficient to represent the most common haplotypes, and (3) typing of tSNPs in clinical material. Early definitions of tSNPs focused on the amount of haplotype diversity that they explained. To select tSNPs that would have maximal power in a genetic association study, however, we have developed optimization criteria based on the r2 measure of association and have compared these with other criteria based on the haplotype diversity. To evaluate the full program and to assess how well the selected tags are likely to perform, we have determined the haplotype structure and have assessed tSNPs in the SCN1A gene, an important candidate gene for sporadic epilepsy. We find that as few as four tSNPs are predicted to maintain a consistently high r2 value with all other common SNPs in the gene, indicating that the tags could be used in an association study with only a modest reduction in power relative to direct assays of all common SNPs. This implies that very large case-control studies can be screened for variation in hundreds of candidate genes with manageable experimental effort, once tSNPs are identified. However, our results also show that tSNPs identified in one population may not necessarily perform well in another, indicating that the preliminary study to identify tSNPs and the later case-control study should be performed in the same population. Our results also indicate that tSNPs will not easily identify discrepant SNPs, which lie on importantly discriminating but apparently short genealogical branches. This could significantly complicate tagging approaches for phenotypes influenced by variants that have experienced positive selection.     

Pathways-based analyses of whole-genome association study data in bipolar disorder reveal genes mediating ion channel activity and synaptic neurotransmission

Despite known heritability, the complex genetic architecture of bipolar disorder (likely including trait, locus and allelic heterogeneity, as well as genetic interactions) has confounded genetic discovery for many years. Even modern day whole genome association studies (WGAS) using over half a million common SNPs have implicated only a handful of genes at the genomewide level. Temporally coincident with this series of WGAS, a host of pathways-based analyses (PBAs) have emerged as novel computational approaches in the examination of large-scale datasets, but thus far rarely have been applied to WGAS data in psychiatric disorders. Here, we report a series of PBAs conducted using exploratory visual analysis, an analytic and visualization software tool for examining genomic data, to examine results from the National Institutes of Mental Health and Wellcome-Trust Case Control Consortium WGAS in bipolar disorder. Consistent with a host of prior linkage findings, some candidate gene association studies, and recent WGAS, our strongest findings suggest involvement of ion channel structural and regulatory genes, including voltage-gated ion channels and the broader ion channel group that comprises both voltage- and ligand-gated channels. Moreover, we found only modest overlap in the particular genes driving the significance of these gene sets across the analyses. This observation strongly suggests that variation in ion channel genes, as a class of genes, may contribute to the susceptibility of bipolar disorder and that heterogeneity may figure prominently in the genetic architecture of this susceptibility. Electronic supplementary material  The online version of this article (doi:) contains supplementary material, which is available to authorized users.     

Genomewide scan for affective disorder susceptibility Loci in families of a northern Swedish isolated population

We analyzed nine multigenerational families with ascertained affective spectrum disorders in northern Sweden's geographically isolated population of Vasterbotten. This northern Swedish population, which originated from a limited number of early settlers approximately 8,000 years ago, is genetically more homogeneous than outbred populations. In a genomewide linkage analysis, we identified three chromosomal loci with multipoint LOD scores (MPLOD) >/=2 at 9q31.1-q34.1 (MPLOD 3.24), 6q22.2-q24.2 (MPLOD 2.48), and 2q33-q36 (MPLOD 2.26) under a recessive affected-only model. Follow-up genotyping with application of a 2-cM density simple-tandem-repeat (STR) map confirmed linkage at 9q31.1-q34.1 (MPLOD 3.22), 6q23-q24 (MPLOD 3.25), and 2q33-q36 (MPLOD 2.2). In an initial analysis aimed at identification of the underlying susceptibility genes, we focused our attention on the 9q locus. We fine mapped this region at a 200-kb STR density, with the result of an MPLOD of 3.70. Genealogical studies showed that three families linked to chromosome 9q descended from common founder couples approximately 10 generations ago. In this approximately 10-generation pedigree, a common ancestral haplotype was inherited by the patients, which reduced the 9q candidate region to 1.6 Mb. Further, the shared haplotype was observed in 4.2% of patients with bipolar disorder with alternating episodes of depression and mania, but it was not observed in control individuals in a patient-control sample from the Vasterbotten isolate. These results suggest a susceptibility locus on 9q31-q33 for affective disorder in this common ancestral region.     

Family-based analysis using a dense single-nucleotide polymorphism-based map defines genetic variation at PSORS1, the major psoriasis-susceptibility locus

Psoriasis is a common skin disorder of multifactorial origin. Genomewide scans for disease susceptibility have repeatedly demonstrated the existence of a major locus, PSORS1 (psoriasis susceptibility 1), contained within the major histocompatibility complex (MHC), on chromosome 6p21. Subsequent refinement studies have highlighted linkage disequilibrium (LD) with psoriasis, along a 150-kb segment that includes at least three candidate genes (encoding human leukocyte antigen-C [HLA-C], alpha-helix-coiled-coil-rod homologue, and corneodesmosin), each of which has been shown to harbor disease-associated alleles. However, the boundaries of the minimal PSORS1 region remain poorly defined. Moreover, interpretations of allelic association with psoriasis are compounded by limited insight of LD conservation within MHC class I interval. To address these issues, we have pursued a high-resolution genetic characterization of the PSORS1 locus. We resequenced genomic segments along a 220-kb region at chromosome 6p21 and identified a total of 119 high-frequency SNPs. Using 59 SNPs (18 coding and 41 noncoding SNPs) whose position was representative of the overall marker distribution, we genotyped a data set of 171 independently ascertained parent-affected offspring trios. Family-based association analysis of this cohort highlighted two SNPs (n.7 and n.9) respectively lying 7 and 4 kb proximal to HLA-C. These markers generated highly significant evidence of disease association (P<10-9), several orders of magnitude greater than the observed significance displayed by any other SNP that has previously been associated with disease susceptibility. This observation was replicated in a Gujarati Indian case/control data set. Haplotype-based analysis detected overtransmission of a cluster of chromosomes, which probably originated by ancestral mutation of a common disease-bearing haplotype. The only markers exclusive to the overtransmitted chromosomes are SNPs n.7 and n.9, which define a 10-kb PSORS1 core risk haplotype. These data demonstrate the power of SNP haplotype-based association analyses and provide high-resolution dissection of genetic variation across the PSORS1 interval, the major susceptibility locus for psoriasis.     

SimPed: a simulation program to generate haplotype and genotype data for pedigree structures

With the widespread availability of SNP genotype data, there is great interest in analyzing pedigree haplotype data. Intermarker linkage disequilibrium for microsatellite markers is usually low due to their physical distance; however, for dense maps of SNP markers, there can be strong linkage disequilibrium between marker loci. Linkage analysis (parametric and nonparametric) and family-based association studies are currently being carried out using dense maps of SNP marker loci. Monte Carlo methods are often used for both linkage and association studies; however, to date there are no programs available which can generate haplotype and/or genotype data consisting of a large number of loci for pedigree structures. SimPed is a program that quickly generates haplotype and/or genotype data for pedigrees of virtually any size and complexity. Marker data either in linkage disequilibrium or equilibrium can be generated for greater than 20,000 diallelic or multiallelic marker loci. Haplotypes and/or genotypes are generated for pedigree structures using specified genetic map distances and haplotype and/or allele frequencies. The simulated data generated by SimPed is useful for a variety of purposes, including evaluating methods that estimate haplotype frequencies for pedigree data, evaluating type I error due to intermarker linkage disequilibrium and estimating empirical p values for linkage and family-based association studies.     

A 6.9-Mb high-resolution BAC/PAC contig of human 4p15.3-p16.1, a candidate region for bipolar affective disorder

Bipolar affective disorder (BPAD) is a complex disease with a significant genetic component and a population lifetime risk of 1%. Our previous work identified a region of human chromosome 4p that showed significant linkage to BPAD in a large pedigree. Here, we report the construction of an accurate, high-resolution physical map of 6.9 Mb of human chromosome 4p15.3-p16.1, which includes an 11-cM (5.8 Mb) critical region for BPAD. The map consists of 460 PAC and BAC clones ordered by a combination of STS content analysis and restriction fragment fingerprinting, with a single approximately 300-kb gap remaining. A total of 289 new and existing markers from a wide range of sources have been localized on the contig, giving an average marker resolution of 1 marker/23 kb. The STSs include 57 ESTs, 9 of which represent known genes. This contig is an essential preliminary to the identification of candidate genes that predispose to bipolar affective disorder, to the completion of the sequence of the region, and to the development of a high-density SNP map.     

Haplotypes at ATM identify coding-sequence variation and indicate a region of extensive linkage disequilibrium

Genetic variation in the human population may lead to functional variants of genes that contribute to risk for common chronic diseases such as cancer. In an effort to detect such possible predisposing variants, we constructed haplotypes for a candidate gene and tested their efficacy in association studies. We developed haplotypes consisting of 14 biallelic neutral-sequence variants that span 142 kb of the ATM locus. ATM is the gene responsible for the autosomal recessive disease ataxia-telangiectasia (AT). These ATM noncoding single-nucleotide polymorphisms (SNPs) were genotyped in nine CEPH families (89 individuals) and in 260 DNA samples from four different ethnic origins. Analysis of these data with an expectation-maximization algorithm revealed 22 haplotypes at this locus, with three major haplotypes having frequencies > or = .10. Tests for recombination and linkage disequilibrium (LD) show reduced recombination and extensive LD at the ATM locus, in all four ethnic groups studied. The most striking example was found in the study population of European ancestry, in which no evidence for recombination could be discerned. The potential of ATM haplotypes for detection of genetic variants through association studies was tested by analysis of 84 individuals carrying one of three ATM coding SNPs. Each coding SNP was detected by association with an ATM haplotype. We demonstrate that association studies with haplotypes for candidate genes have significant potential for the detection of genetic backgrounds that contribute to disease.     

Pleiotropic effects of a chromosome 3 locus on speech-sound disorder and reading

Speech-sound disorder (SSD) is a complex behavioral disorder characterized by speech-sound production errors associated with deficits in articulation, phonological processes, and cognitive linguistic processes. SSD is prevalent in childhood and is comorbid with disorders of language, spelling, and reading disability, or dyslexia. Previous research suggests that developmental problems in domains associated with speech and language acquisition place a child at risk for dyslexia. Recent genetic studies have identified several candidate regions for dyslexia, including one on chromosome 3 segregating in a large Finnish pedigree. To explore common genetic influences on SSD and reading, we examined linkage for several quantitative traits to markers in the pericentrometric region of chromosome 3 in 77 families ascertained through a child with SSD. The quantitative scores measured several processes underlying speech-sound production, including phonological memory, phonological representation, articulation, receptive and expressive vocabulary, and reading decoding and comprehension skills. Model-free linkage analysis was followed by identification of sib pairs with linkage and construction of core shared haplotypes. In our multipoint analyses, measures of phonological memory demonstrated the strongest linkage (marker D3S2465, P=5.6 x 10(-5), and marker D3S3716, P=6.8 x 10(-4)). Tests for single-word decoding also demonstrated linkage (real word reading: marker D3S2465, P=.004; nonsense word reading: marker D3S1595, P=.005). The minimum shared haplotype in sib pairs with similar trait values spans 4.9 cM and is bounded by markers D3S3049 and D3S3045. Our results suggest that domains common to SSD and dyslexia are pleiotropically influenced by a putative quantitative trait locus on chromosome 3.     

Haplotype-specific <Emphasis Type="Italic">MAPT</Emphasis> exon 3 expression regulated by common intronic polymorphisms associated with Parkinsonian disorders

Background

Genome wide association studies have identified microtubule associated protein tau (MAPT) H1 haplotype single nucleotide polymorphisms (SNPs) as leading common risk variants for Parkinson’s disease, progressive supranuclear palsy and corticobasal degeneration. The MAPT risk variants fall within a large 1.8 Mb region of high linkage disequilibrium, making it difficult to discern the functionally important risk variants. Here, we leverage the strong haplotype-specific expression of MAPT exon 3 to investigate the functionality of SNPs that fall within this H1 haplotype region of linkage disequilibrium.

Methods

In this study, we dissect the molecular mechanisms by which haplotype-specific SNPs confer allele-specific effects on the alternative splicing of MAPT exon 3. Firstly, we use haplotype-hybrid whole-locus genomic MAPT vectors studies to identify functional SNPs. Next, we characterise the RNA-protein interactions at two loci by mass spectrometry. Lastly, we knockdown candidate splice factors to determine their effect on MAPT exon 3 using a novel allele-specific qPCR assay.

Results

Using whole-locus genomic DNA expression vectors to express MAPT haplotype variants, we demonstrate that rs17651213 regulates exon 3 inclusion in a haplotype-specific manner. We further investigated the functionality of this region using RNA-electrophoretic mobility shift assays to show differential RNA-protein complex formation at the H1 and H2 sequence variants of SNP rs17651213 and rs1800547 and subsequently identified candidate trans-acting splicing factors interacting with these functional SNPs sequences by RNA-protein pull-down experiment and mass spectrometry. Finally, gene knockdown of candidate splice factors identified by mass spectrometry demonstrate a role for hnRNP F and hnRNP Q in the haplotype-specific regulation of exon 3 inclusion.

Conclusions

We identified common splice factors hnRNP F and hnRNP Q regulating the haplotype-specific splicing of MAPT exon 3 through intronic variants rs1800547 and rs17651213. This work demonstrates an integrated approach to characterise the functionality of risk variants in large regions of linkage disequilibrium.

     

Genomic Landscape of a Three-Generation Pedigree Segregating Affective Disorder

Bipolar disorder (BPD) is a common psychiatric illness with a complex mode of inheritance. Besides traditional linkage and association studies, which require large sample sizes, analysis of common and rare chromosomal copy number variants (CNVs) in extended families may provide novel insights into the genetic susceptibility of complex disorders. Using the Illumina HumanHap550 BeadChip with over 550,000 SNP markers, we genotyped 46 individuals in a three-generation Old Order Amish pedigree with 19 affected (16 BPD and three major depression) and 27 unaffected subjects. Using the PennCNV algorithm, we identified 50 CNV regions that ranged in size from 12 to 885 kb and encompassed at least 10 single nucleotide polymorphisms (SNPs). Of 19 well characterized CNV regions that were available for combined genotype-expression analysis 11 (58%) were associated with expression changes of genes within, partially within or near these CNV regions in fibroblasts or lymphoblastoid cell lines at a nominal P value <0.05. To further investigate the mode of inheritance of CNVs in the large pedigree, we analyzed a set of four CNVs, located at 6q27, 9q21.11, 12p13.31 and 15q11, all of which were enriched in subjects with affective disorders. We additionally show that these variants affect the expression of neuronal genes within or near the rearrangement. Our analysis suggests that family based studies of the combined effect of common and rare CNVs at many loci may represent a useful approach in the genetic analysis of disease susceptibility of mental disorders.