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.     

LRRK2 G2019S in families with Parkinson disease who originated from Europe and the Middle East: evidence of two distinct founding events beginning two millennia ago

The leucine-rich repeat kinase 2 (LRRK2) G2019S mutation is the most common genetic determinant of Parkinson disease (PD) identified to date. It accounts for 1%-7% of PD in patients of European origin and 20%-40% in Ashkenazi Jews and North African Arabs with PD. Previous studies concluded that patients from these populations all shared a common Middle Eastern founder who lived in the 13th century. We tested this hypothesis by genotyping 25 microsatellite and single-nucleotide-polymorphism markers in 22 families with G2019S and observed two distinct haplotypes. Haplotype 1 was present in 19 families of Ashkenazi Jewish and European ancestry, whereas haplotype 2 occurred in three European American families. Using a maximum-likelihood method, we estimated that the families with haplotype 1 shared a common ancestor 2,250 (95% confidence interval 1,650-3,120) years ago, whereas those with haplotype 2 appeared to share a more recent founder. Our data suggest two separate founding events for G2019S in these populations, beginning at a time that coincides with the Jewish Diasporas.     

Common SNP-based haplotype analysis of the 4p16.3 Huntington disease gene region

Age at the onset of motor symptoms in Huntington disease (HD) is determined largely by the length of a CAG repeat expansion in HTT but is also influenced by other genetic factors. We tested whether common genetic variation near the mutation site is associated with differences in the distribution of expanded CAG alleles or age at the onset of motor symptoms. To define disease-associated single-nucleotide polymorphisms (SNPs), we compared 4p16.3 SNPs in HD subjects with population controls in a case:control strategy, which revealed that the strongest signals occurred at a great distance from the HD mutation as a result of "synthetic association" with SNP alleles that are of low frequency in population controls. Detailed analysis delineated a prominent ancestral haplotype that accounted for ～50% of HD chromosomes and extended to at least 938 kb on about half of these. Together, the seven most abundant haplotypes accounted for ～83% of HD chromosomes. Neither the extended shared haplotype nor the individual local HTT haplotypes were associated with altered CAG-repeat length distribution or residual age at the onset of motor symptoms, arguing against modification of these disease features by common cis-regulatory elements. Similarly, the 11 most frequent control haplotypes showed no trans-modifier effect on age at the onset of motor symptoms. Our results argue against common local regulatory variation as a factor influencing HD pathogenesis, suggesting that genetic modifiers be sought elsewhere in the genome. They also indicate that genome-wide association analysis with a small number of cases can be effective for regional localization of genetic defects, even when a founder effect accounts for only a fraction of the disorder.     

Geographical distribution of haplotypes in Swedish families with Huntington's disease

This study was planned to determine the number of origins of the mutation underlying Huntington's disease (HD) in Sweden. Haplotypes were constructed for 23 different HD families, using six different polymorphisms [(CCG) _n , GT70, 674, BS1, E2 and 4.2], including two within the gene. In addition, extensive genealogical investigations were performed, and the geographical origin of the haplotypes was studied. Ten different haplotypes were observed suggesting multiple origins for the HD mutation in Sweden. Analysis of the two polymorphic markers within the HD gene (the CCG repeat and GT70) indicates that there are at least three origins for the HD mutation in Sweden. One of these haplotypes (7/A) accounts for 89% of the families, suggesting that the majority of the Swedish HD families are related through a single HD mutation of ancient origin. Furthermore, three of the families that were previously considered to be unrelated could be traced to a common ancestor in the 15th century, a finding that is consistent with this hypothesis.     

Dominant-Negative Effects of Adult-Onset Huntingtin Mutations Alter the Division of Human Embryonic Stem Cells-Derived Neural Cells

Mutations of the huntingtin protein (HTT) gene underlie both adult-onset and juvenile forms of Huntington’s disease (HD). HTT modulates mitotic spindle orientation and cell fate in mouse cortical progenitors from the ventricular zone. Using human embryonic stem cells (hESC) characterized as carrying mutations associated with adult-onset disease during pre-implantation genetic diagnosis, we investigated the influence of human HTT and of an adult-onset HD mutation on mitotic spindle orientation in human neural stem cells (NSCs) derived from hESCs. The RNAi-mediated silencing of both HTT alleles in neural stem cells derived from hESCs disrupted spindle orientation and led to the mislocalization of dynein, the p150^Glued subunit of dynactin and the large nuclear mitotic apparatus (NuMA) protein. We also investigated the effect of the adult-onset HD mutation on the role of HTT during spindle orientation in NSCs derived from HD-hESCs. By combining SNP-targeting allele-specific silencing and gain-of-function approaches, we showed that a 46-glutamine expansion in human HTT was sufficient for a dominant-negative effect on spindle orientation and changes in the distribution within the spindle pole and the cell cortex of dynein, p150^Glued and NuMA in neural cells. Thus, neural derivatives of disease-specific human pluripotent stem cells constitute a relevant biological resource for exploring the impact of adult-onset HD mutations of the HTT gene on the division of neural progenitors, with potential applications in HD drug discovery targeting HTT-dynein-p150^Glued complex interactions.     

Genomic control for association studies: a semiparametric test to detect excess-haplotype sharing

Individuals who share a disease mutation from a common ancestor often share alleles at genetic markers adjacent to the mutation, even if the common ancestor is remote. The alleles at these adjacent markers, called the haplotype, can be visualized as a string of realizations of random variables, which may be dependent when individuals are related in some fashion. Ideally, for a sample of individuals all having the same (genetic) disease, this dependence-measured as haplotype-sharing-will be greater in the vicinity of disease genes than in other regions of the genome. In this paper we present a semiparametric test for haplotype-sharing. We begin by developing a model assuming that the ancestral haplotype is known and thus the extent of haplotype-sharing from a common ancestor can be determined unambiguously. The amount of overlap at markers far from the disease is treated as a random variable with an unknown distribution F, which we estimate non-parametrically. Overlap of markers surrounding disease genes are modeled as a mixture pF(x - theta) + (1 - p)F(x), in which p is the fraction of subjects with the disease mutation. Testing for a disease gene then amounts to testing whether p = 0. Next we drop the assumption that the ancestral haplotype is known. To detect excess clustering of haplotypes, we measure the pairwise overlap of a set of haplotypes. As in the simpler scenario, this distribution is modeled as a location-shift mixture. To test the hypothesis we construct a score test with a simple limiting distribution.     

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.

     

The reliability of haplotyping inference in nuclear families: misassignment rates for SNPs and microsatellites

Single nucleotide polymorphisms (SNPs) are widely used when investigators try to map complex disease genes. Although biallelic SNP markers are less informative than microsatellite markers, one can increase their information content by using haplotypes. However, assigning haplotypes (i.e., assigning phase) correctly can be problematic in the presence of SNP heterozygosity. For example, a doubly heterozygous individual, with genotype 12, 12, could have haplotypes 1-1/2-2 or 1-2/2-1 with equal probability; in the absence of additional information, there is no way to determine which haplotype is correct. Thus an algorithm that assigns haplotypes to such an individual will assign the wrong one 50% of the time. We have studied the frequency of haplotype misassignments, i.e., haplotypes that are misassigned solely because of inherent marker ambiguity (not because of errors in genotyping or calculation). We examined both SNPs and microsatellite markers. We used the computer programs GENEHUNTER and SIMWALK to assign the haplotypes. We simulated (a) families with 1-5 children, (b) haplotypes involving different numbers of marker loci (3, 5, 7 and 10 loci, all in linkage equilibrium), and (c) different allele frequencies. Misassignment rates are highest (a) in small families, (b) with many SNP loci, and (c) for loci with the greatest heterozygosity (i.e., where both alleles have frequency 0.5). For example, for triads (i.e., one-child families with both parents genotyped), misassignment rates for SNPs can reach almost 50%. Family sizes of 4-5 children are required in order to ensure a misassignment frequency of < or = 5% for ten-SNP haplotypes with allele frequencies of 0.25-0.5. For microsatellites, a family size of at least 2-3 children is necessary to keep haplotyping misassignments < or = 5%. Finally, we point out that it is misleading for a computer program to yield haplotype assignments without indicating that they may have been misassigned, and we discuss the implications of these misassignments for association and linkage analysis.     

Cu, Zn-superoxide dismutase 1 (SOD1) is a novel target of Puromycin-sensitive aminopeptidase (PSA/NPEPPS): PSA/NPEPPS is a possible modifier of amyotrophic lateral sclerosis

Background

Huntington disease (HD) is caused by a polyglutamine expansion of more than 35 units in the huntingtin protein. This expanded repeat length inversely correlates with the age-at-onset (AAO), however, additional genetic factors apart from the expanded CAG repeat size are thought to influence the course and the AAO in HD. Until now, among others, the gene encoding PCG-1α (PPARGC1A) was shown to modify the AAO in two independent, however small, populations. PGC-1α is involved in the induction of various mechanisms regulating mitochondrial biogenesis and oxidative stress defence. Furthermore, several studies have linked impairment of its function and/or its expression to HD pathogenesis. As the identification of distinct modifiers in association studies is largely dependent on the size of the observed population, we investigated nine different single nucleotide polymorphisms (SNPs) in PPARGC1A in order to replicate the disease modifying effect in more than 800 European HD patients and to identify an association with AAO in HD.

Results

Two SNPs, one in the promoter and one in the transcribed region of the gene, showed a significant effect on the AAO. While the minor allele of SNP rs7665116 (g.38570C), located in the transcribed gene region, was associated with a delay in disease onset, especially in HD patients with Italian ancestry, the minor allele of SNP rs2970870 (g.-1437C) in the promoter region leads to an earlier onset of HD in its homozygous state. Additionally, global testing of haplotype block 2, which covers the main part of the transcribed region of the gene, revealed an association between block 2 haplotypes and the disease onset.

Conclusion

Therefore, our results indicate opposing modifying influences of two SNPs within one gene on AAO and support the idea that PGC-1α dysfunction is involved in HD pathology.     

Nucleotide Diversity and Haplotype Structure of the Human Angiotensinogen Gene in Two Populations

Variation in the angiotensinogen gene, AGT, has been associated with variation in plasma angiotensinogen levels. In addition, the T235M polymorphism in the AGT product is associated with an increased risk of essential hypertension in multiple populations, making AGT a good example of a quantitative-trait locus underlying susceptibility to a common disease. To better understand genetic variation in AGT, we sequenced a 14.4-kb genomic region spanning the entire AGT and identified 44 single-nucleotide polymorphisms (SNPs). Forty-two SNPs were observed both in 88 white and in 77 Japanese unselected subjects. Six major haplotypes accounted for most of the variation in this region, indicating less allelic complexity than in many other genomic regions. Although the two populations were found to share all of the major AGT haplotypes, there were substantial differences in haplotype frequencies. Pairwise linkage disequilibrium (LD), measured by the D', r(2), and d(2) statistics, demonstrated a general pattern of decline with increasing distance, but, as expected in a small genomic region, individual LD values were highly variable. LD between T235M and each of the other 39 SNPs was assessed in order to model the usefulness of LD to detect a disease-associated mutation. Among the Japanese subjects, 13 (33%) of the SNPs had r(2) values >0.1, whereas this statistic was substantially higher for the white subjects (occurring in 35/39 [90%]). LD between a hypertension-associated promoter mutation, A-6G, and 39 SNPs was also measured. Similar results were obtained, with 33% of the SNPs showing r(2)>0.1 in the Japanese subjects and 92% of the SNPs showing r(2)>0.1 in the white subjects. This difference, which occurs despite an overall similarity in LD patterns in the two populations, reflects a much higher frequency of the M235-associated haplotype in the white sample. These results have important implications for the usefulness of LD approaches in the mapping of genes underlying susceptibility to complex diseases.     

PTPN22 genetic variation: evidence for multiple variants associated with rheumatoid arthritis

The minor allele of the R620W missense single-nucleotide polymorphism (SNP) (rs2476601) in the hematopoietic-specific protein tyrosine phosphatase gene, PTPN22, has been associated with multiple autoimmune diseases, including rheumatoid arthritis (RA). These genetic data, combined with biochemical evidence that this SNP affects PTPN22 function, suggest that this phosphatase is a key regulator of autoimmunity. To determine whether other genetic variants in PTPN22 contribute to the development of RA, we sequenced the coding regions of this gene in 48 white North American patients with RA and identified 15 previously unreported SNPs, including 2 coding SNPs in the catalytic domain. We then genotyped 37 SNPs in or near PTPN22 in 475 patients with RA and 475 individually matched controls (sample set 1) and selected a subset of markers for replication in an additional 661 patients with RA and 1,322 individually matched controls (sample set 2). Analyses of these results predict 10 common (frequency >1%) PTPN22 haplotypes in white North Americans. The sole haplotype found to carry the previously identified W620 risk allele was strongly associated with disease in both sample sets, whereas another haplotype, identical at all other SNPs but carrying the R620 allele, showed no association. R620W, however, does not fully explain the association between PTPN22 and RA, since significant differences between cases and controls persisted in both sample sets after the haplotype data were stratified by R620W. Additional analyses identified two SNPs on a single common haplotype that are associated with RA independent of R620W, suggesting that R620W and at least one additional variant in the PTPN22 gene region influence RA susceptibility.     

Deregulated Splicing Is a Major Mechanism of RNA-Induced Toxicity in Huntington's Disease

Huntington's disease (HD) is caused by an expanded CAG repeat in the huntingtin (HTT) gene, translating into an elongated polyglutamine stretch. In addition to the neurotoxic mutant HTT protein, the mutant CAG repeat RNA can exert toxic functions by trapping RNA-binding proteins. While few examples of proteins that aberrantly bind to mutant HTT RNA and execute abnormal function in conjunction with the CAG repeat RNA have been described, an unbiased approach to identify the interactome of mutant HTT RNA is missing. Here, we describe the analysis of proteins that preferentially bind mutant HTT RNA using a mass spectrometry approach. We show that (I) the majority of proteins captured by mutant HTT RNA belong to the spliceosome pathway, (II) expression of mutant CAG repeat RNA induces mis-splicing in a HD cell model, (III) overexpression of one of the splice factors trapped by mutant HTT ameliorates the HD phenotype in a fly model and (VI) deregulated splicing occurs in human HD brain. Our data suggest that deregulated splicing is a prominent mechanism of RNA-induced toxicity in HD.     

Association Analysis of GABRB3 Promoter Variants with Heroin Dependence

GABRB3 encoding the β3 subunit of GABAA receptor has been implicated in multiple neuropsychiatric disorders, including substance abuse. Previous studies reported that SNPs at the 5′ regulatory region of GABRB3 could regulate GABRB3 gene expression and associated with childhood absence epilepsy (CAE). The study aimed to investigate whether SNPs at the 5′ regulatory region of GABRB3 were associated with heroin dependence in our population. We first re-sequenced 1.5 kb of the 5′regulatory region of GABRB3 gene to examine the SNP profile in the genomic DNA of 365 control subjects. Then, we conducted a case-control association analysis between 576 subjects with heroin dependence (549 males, 27 females) and 886 controls (472 males, 414 females) by genotyping the rs4906902 as a tag SNP. We also conducted a reporter gene assay to assess the promoter activity of two major haplotypes derived from SNPs at this region. We detected 3 common SNPs (rs4906902, rs8179184 and rs20317) at this region that had strong pair-wise linkage disequilibrium. The C allele of rs4906902 was found to be associated with increased risk of heroin dependence (odds ratio:1.27, p = 0.002). Two major haplotypes (C-A-G and T-G-C) derived from these 3 SNPs accounted for 99% of this sample, and reporter gene activity assay showed that haplotype C-A-G that contained the C allele of the tag SNP rs4906902 had higher activity than haplotype T-G-C. Our data suggest that GABRB3 might be associated with heroin dependence, and increased expression of GABRB3 might contribute to the pathogenesis of heroin dependence.     

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.     

Structure of tumor necrosis factor-alpha haploblocks in European populations

DNA variants in the tumor necrosis factor-α (TNF) and linked lymphotoxin-α genes, and specific alleles of the highly polymorphic human leukocyte antigen B (HLA-B) gene have been implicated in a plethora of immune and infectious diseases. However, the tight linkage disequilibrium characterizing the central region of the human major histocompatibility complex (MHC) containing these gene loci has made difficult the unequivocal interpretation of genetic association data. To alleviate these difficulties and facilitate the design of more focused follow-up studies, we investigated the structure and distribution of HLA-B-specific MHC haplotypes reconstructed in a European population from unphased genotypes at a set of 25 single nucleotide polymorphism sites spanning a 66-kilobase long region across TNF. Consistent with the published data, we found limited genetic diversity across the so-called TNF block, with the emergence of seven common MHC haplotypes, termed TNF block super-haplotypes. We also found that the ancestral haplotype 8.1 shares a TNF block haplotype with HLA-B*4402. HLA-B*5701, a known protective allele in HIV-1 pathogenesis, occurred in a unique TNF block haplotype.