Nucleotide variation,haplotype structure,and association with end-stage renal disease of the human interleukin-1 gene cluster

A dense gene-based SNP map was constructed across a 360-kb region containing the interleukin-1 gene cluster (IL1A, IL1B, and IL1RN), focusing on IL1RN. In total, 95 polymorphisms were confirmed or identified primarily by direct sequencing. Polymorphisms were precisely mapped to completed BAC and genomic sequences spanning this region. The polymorphisms were typed in 443 case-control subjects from Caucasian and African American groups. Consecutive pair-wise marker linkage disequilibrium was not strictly correlated with distance and ranged from D'=0.0079 to 1.000 and D'=0.0521 to 1.0000 in Caucasians and African Americans, respectively. Single markers and haplotypes in IL1 cluster genes were evaluated for association with end-stage renal disease (ESRD). Eleven SNPs show some evidence of association with ESRD, with the strongest associations in two IL1A variants, one SNP, rs1516792-3, in intron 5 (p=0.0015) and a 4-bp insertion/deletion within the 3'UTR, rs16347-2 (p=0.0024), among African Americans with non-T2DM-associated ESRD.     

Identification of a common risk haplotype for canine idiopathic epilepsy in the ADAM23 gene

Background

Idiopathic epilepsy is a common neurological disease in human and domestic dogs but relatively few risk genes have been identified to date. The seizure characteristics, including focal and generalised seizures, are similar between the two species, with gene discovery facilitated by the reduced genetic heterogeneity of purebred dogs. We have recently identified a risk locus for idiopathic epilepsy in the Belgian Shepherd breed on a 4.4 megabase region on CFA37.

Results

We have expanded a previous study replicating the association with a combined analysis of 157 cases and 179 controls in three additional breeds: Schipperke, Finnish Spitz and Beagle (p_c = 2.9e–07, p_GWAS = 1.74E-02). A targeted resequencing of the 4.4 megabase region in twelve Belgian Shepherd cases and twelve controls with opposite haplotypes identified 37 case-specific variants within the ADAM23 gene. Twenty-seven variants were validated in 285 cases and 355 controls from four breeds, resulting in a strong replication of the ADAM23 locus (p_raw = 2.76e–15) and the identification of a common 28 kb-risk haplotype in all four breeds. Risk haplotype was present in frequencies of 0.49–0.7 in the breeds, suggesting that ADAM23 is a low penetrance risk gene for canine epilepsy.

Conclusions

These results implicate ADAM23 in common canine idiopathic epilepsy, although the causative variant remains yet to be identified. ADAM23 plays a role in synaptic transmission and interacts with known epilepsy genes, LGI1 and LGI2, and should be considered as a candidate gene for human epilepsies.

Electronic supplementary material

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

Prioritization of SNPs for genome-wide association studies using an interaction model of genetic variation,gene expression,and trait variation

The identification of true causal loci to unravel the statistical evidence of genotype-phenotype correlations and the biological relevance of selected single-nucleotide polymorphisms (SNPs) is a challenging issue in genome-wide association studies (GWAS). Here, we introduced a novel method for the prioritization of SNPs based on p-values from GWAS. The method uses functional evidence from populations, including phenotype-associated gene expressions. Based on the concept of genetic interactions, such as perturbation of gene expression by genetic variation, phenotype and gene expression related SNPs were prioritized by adjusting the p-values of SNPs. We applied our method to GWAS data related to drug-induced cytotoxicity. Then, we prioritized loci that potentially play a role in druginduced cytotoxicity. By generating an interaction model, our approach allowed us not only to identify causal loci, but also to find intermediate nodes that regulate the flow of information among causal loci, perturbed gene expression, and resulting phenotypic variation.     

Effects of population structure on genetic association studies

Population-based case-control association is a promising approach for unravelling the genetic basis of complex diseases. One potential problem of this approach is the presence of population structure in the samples. Using the Collaborative Study on the Genetics of Alcoholism (COGA) single-nucleotide polymorphism (SNP) datasets, we addressed three questions: How can the degree of population structure be quantified, and how does the population structure affect association studies? How accurate and efficient is the genomic control method in correcting for population structure? The amount of population structure in the COGA SNP data was found to inflate the p-value in association tests. Genomic control was found to be effective only when the appropriate number of markers was used in the control group in order to correctly calibrate the test. The approach presented in this paper could be used to select the appropriate number of markers for use in the genomic control method of correcting population structure.     

Assessment of the functionality of genome‐wide canine SNP arrays and implications for canine disease association studies

Domestic dogs share a wide range of important disease conditions with humans, including cancers, diabetes and epilepsy. Many of these conditions have similar or identical underlying pathologies to their human counterparts and thus dogs represent physiologically relevant natural models of human disorders. Comparative genomic approaches whereby disease genes can be identified in dog diseases and then mapped onto the human genome are now recognized as a valid method and are increasing in popularity. The majority of dog breeds have been created over the past few hundred years and, as a consequence, the dog genome is characterized by extensive linkage disequilibrium (LD), extending usually from hundreds of kilobases to several megabases within a breed, rather than tens of kilobases observed in the human genome. Genome‐wide canine SNP arrays have been developed, and increasing success of using these arrays to map disease loci in dogs is emerging. No equivalent of the human HapMap currently exists for different canine breeds, and the LD structure for such breeds is far less understood than for humans. This study is a dedicated large‐scale assessment of the functionalities (LD and SNP tagging performance) of canine genome‐wide SNP arrays in multiple domestic dog breeds. We have used genotype data from 18 breeds as well as wolves and coyotes genotyped by the Illumina 22K canine SNP array and Affymetrix 50K canine SNP array. As expected, high tagging performance was observed with most of the breeds using both Illumina and Affymetrix arrays when multi‐marker tagging was applied. In contrast, however, large differences in population structure, LD coverage and pairwise tagging performance were found between breeds, suggesting that study designs should be carefully assessed for individual breeds before undertaking genome‐wide association studies (GWAS).     

Identification and genetic effect of haplotype in the bovine BMP7 gene

Bone morphogenetic proteins (BMPs) are peptide growth factors belonging to the transforming growth factor-beta (TGF-β) superfamily, and some members of the BMP family support white adipocyte differentiation. In this study, we focused on the BMP7 which singularly promotes the differentiation of brown preadipocytes. Haplotypes involving 5 single nucleotide polymorphism (SNP) sites in the bovine BMP7 gene were identified and their effect on body weight was analyzed. 16 haplotypes and 18 combined haplotypes were revealed and the linkage disequilibrium was assessed in the cattle population with 602 individuals representing three main cattle breeds from China. The results showed that haplotypes 3, 10 and 14 were predominant and accounted for 75.64%, 69.85%, and 83.36% in Nanyang, Qinchuan and Jiaxian cattle breeds, respectively. The statistical analyses indicated that the SNP 1, 4, and 5 are associated with the body weight, body length, and heart girth at 12 and 24 months in Nanyang cattle population (P < 0.05), whereas there is no significant association between their 16 haplotypes and 18 combined haplotypes. Our results provide evidence that some SNPs and haplotypes in BMP7 are associated with growth traits, and may be utilized as a genetic marker in marker-assisted selection for beef cattle breeding programs.     

A comprehensive survey of sequence variation in the ABCA4 (ABCR) gene in Stargardt disease and age-related macular degeneration

Stargardt disease (STGD) is a common autosomal recessive maculopathy of early and young-adult onset and is caused by alterations in the gene encoding the photoreceptor-specific ATP-binding cassette (ABC) transporter (ABCA4). We have studied 144 patients with STGD and 220 unaffected individuals ascertained from the German population, to complete a comprehensive, population-specific survey of the sequence variation in the ABCA4 gene. In addition, we have assessed the proposed role for ABCA4 in age-related macular degeneration (AMD), a common cause of late-onset blindness, by studying 200 affected individuals with late-stage disease. Using a screening strategy based primarily on denaturing gradient gel electrophoresis, we have identified in the three study groups a total of 127 unique alterations, of which 90 have not been previously reported, and have classified 72 as probable pathogenic mutations. Of the 288 STGD chromosomes studied, mutations were identified in 166, resulting in a detection rate of approximately 58%. Eight different alleles account for 61% of the identified disease alleles, and at least one of these, the L541P-A1038V complex allele, appears to be a founder mutation in the German population. When the group with AMD and the control group were analyzed with the same methodology, 18 patients with AMD and 12 controls were found to harbor possible disease-associated alterations. This represents no significant difference between the two groups; however, for detection of modest effects of rare alleles in complex diseases, the analysis of larger cohorts of patients may be required.     

Allelic variation and haplotype structure of the dopamine receptor gene DRD2 in nine Indian populations

The human dopaminergic system is a significant focal point of study in the fields of neuropsychiatry and pharmacology, plus it is also a promising nuclear DNA marker in studies of human genome diversity. In this study, we assayed six polymorphic markers in the dopamine D2 receptor gene (DRD2) in 482 unrelated individuals from nine ethnic populations of India. Our results demonstrate that the six markers are highly polymorphic in all populations and the constructed haplotypes show a high level of heterozygosity. Out of the eight possible three-site haplotypes, all populations commonly shared only three haplotypes. The haplotypes exhibited fairly high frequencies across multiple populations; Kurumba population showed all eight three-site haplotypes. The ancestral haplotype (B2-D2-Al) was observed at high frequency only in the Siddi population. Haplotypes based on all six markers revealed 16 haplotypes, out of which only 6 are most common with a frequency of greater than 5% in at least one of the nine populations. But only three haplotypes were shared by all nine populations with the cumulative frequency ranging from 80.8% (Kurumba) to 96.6% (Onge). Great variation in levels of linkage disequilibrium (LD) was detected, ranging from complete LD in the Badaga to virtually no LD in the Siddi. This range of LD likely reflects different population histories, such as African ancestry in the Siddi and recent founding events in the population isolates, Badaga and Kota.     

Missing data imputation and haplotype phase inference for genome-wide association studies

Imputation of missing data and the use of haplotype-based association tests can improve the power of genome-wide association studies (GWAS). In this article, I review methods for haplotype inference and missing data imputation, and discuss their application to GWAS. I discuss common features of the best algorithms for haplotype phase inference and missing data imputation in large-scale data sets, as well as some important differences between classes of methods, and highlight the methods that provide the highest accuracy and fastest computational performance.     

Analysis of polymorphisms and haplotype structure of the human thymidylate synthase genetic region: a tool for pharmacogenetic studies

5-Fluorouracil (5FU), a widely used chemotherapeutic drug, inhibits the DNA replicative enzyme, thymidylate synthase (Tyms). Prior studies implicated a VNTR (variable numbers of tandem repeats) polymorphism in the 5'-untranslated region (5'-UTR) of the TYMS gene as a determinant of Tyms expression in tumors and normal tissues and proposed that these VNTR genotypes could help decide fluoropyrimidine dosing. Clinical associations between 5FU-related toxicity and the TYMS VNTR were reported, however, results were inconsistent, suggesting that additional genetic variation in the TYMS gene might influence Tyms expression. We thus conducted a detailed genetic analysis of this region, defining new polymorphisms in this gene including mononucleotide (poly A:T) repeats and novel single nucleotide polymorphisms (SNPs) flanking the VNTR in the TYMS genetic region. Our haplotype analysis of this region used data from both established and novel genetic variants and found nine SNP haplotypes accounting for more than 90% of the studied population. We observed non-exclusive relationships between the VNTR and adjacent SNP haplotypes, such that each type of VNTR commonly occurred on several haplotype backgrounds. Our results confirmed the expectation that the VNTR alleles exhibit homoplasy and lack the common ancestry required for a reliable marker of a linked adjacent locus that might govern toxicity. We propose that it may be necessary in a clinical trial to assay multiple types of genetic polymorphisms in the TYMS region to meaningfully model linkage of genetic markers to 5FU-related toxicity. The presence of multiple long (up to 26 nt), polymorphic monothymidine repeats in the promoter region of the sole human thymidylate synthetic enzyme is intriguing.     

Investigation of a genome wide association signal for obesity: synthetic association and haplotype analyses at the melanocortin 4 receptor gene locus

Background

Independent genome-wide association studies (GWAS) showed an obesogenic effect of two single nucleotide polymorphisms (SNP; rs12970134 and rs17782313) more than 150 kb downstream of the melanocortin 4 receptor gene (MC4R). It is unclear if the SNPs directly influence MC4R function or expression, or if the SNPs are on a haplotype that predisposes to obesity or includes functionally relevant genetic variation (synthetic association). As both exist, functionally relevant mutations and polymorphisms in the MC4R coding region and a robust association downstream of the gene, MC4R is an ideal model to explore synthetic association.

Methodology/Principal Findings

We analyzed a genomic region (364.9 kb) encompassing the MC4R in GWAS data of 424 obesity trios (extremely obese child/adolescent and both parents). SNP rs12970134 showed the lowest p-value (p = 0.004; relative risk for the obesity effect allele: 1.37); conditional analyses on this SNP revealed that 7 of 78 analyzed SNPs provided independent signals (p≤0.05). These 8 SNPs were used to derive two-marker haplotypes. The three best (according to p-value) haplotype combinations were chosen for confirmation in 363 independent obesity trios. The confirmed obesity effect haplotype includes SNPs 3′ and 5′ of the MC4R. Including MC4R coding variants in a joint model had almost no impact on the effect size estimators expected under synthetic association.

Conclusions/Significance

A haplotype reaching from a region 5′ of the MC4R to a region at least 150 kb from the 3′ end of the gene showed a stronger association to obesity than single SNPs. Synthetic association analyses revealed that MC4R coding variants had almost no impact on the association signal. Carriers of the haplotype should be enriched for relevant mutations outside the MC4R coding region and could thus be used for re-sequencing approaches. Our data also underscore the problems underlying the identification of relevant mutations depicted by GWAS derived SNPs.     

Sequence variation and haplotype structure at the human HFE locus

The HFE locus encodes an HLA class-I-type protein important in iron regulation and segregates replacement mutations that give rise to the most common form of genetic hemochromatosis. The high frequency of one disease-associated mutation, C282Y, and the nature of this disease have led some to suggest a selective advantage for this mutation. To investigate the context in which this mutation arose and gain a better understanding of HFE genetic variation, we surveyed nucleotide variability in 11.2 kb encompassing the HFE locus and experimentally determined haplotypes. We fully resequenced 60 chromosomes of African, Asian, or European ancestry as well as one chimpanzee, revealing 41 variable sites and a nucleotide diversity of 0.08%. This indicates that linkage to the HLA region has not substantially increased the level of HFE variation. Although several haplotypes are shared between populations, one haplotype predominates in Asia but is nearly absent elsewhere, causing higher than average genetic differentiation among the three major populations. Our samples show evidence of intragenic recombination, so the scarcity of recombination events within the C282Y allele class is consistent with selection increasing the frequency of a young allele. Otherwise, the pattern of variability in this region does not clearly indicate the action of positive selection at this or linked loci.     

Rapid and accurate haplotype phasing and missing-data inference for whole-genome association studies by use of localized haplotype clustering

Whole-genome association studies present many new statistical and computational challenges due to the large quantity of data obtained. One of these challenges is haplotype inference; methods for haplotype inference designed for small data sets from candidate-gene studies do not scale well to the large number of individuals genotyped in whole-genome association studies. We present a new method and software for inference of haplotype phase and missing data that can accurately phase data from whole-genome association studies, and we present the first comparison of haplotype-inference methods for real and simulated data sets with thousands of genotyped individuals. We find that our method outperforms existing methods in terms of both speed and accuracy for large data sets with thousands of individuals and densely spaced genetic markers, and we use our method to phase a real data set of 3,002 individuals genotyped for 490,032 markers in 3.1 days of computing time, with 99% of masked alleles imputed correctly. Our method is implemented in the Beagle software package, which is freely available.     

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.     

Identification of the genetic basis for complex disorders by use of pooling-based genomewide single-nucleotide-polymorphism association studies

We report the development and validation of experimental methods, study designs, and analysis software for pooling-based genomewide association (GWA) studies that use high-throughput single-nucleotide-polymorphism (SNP) genotyping microarrays. We first describe a theoretical framework for establishing the effectiveness of pooling genomic DNA as a low-cost alternative to individually genotyping thousands of samples on high-density SNP microarrays. Next, we describe software called "GenePool," which directly analyzes SNP microarray probe intensity data and ranks SNPs by increased likelihood of being genetically associated with a trait or disorder. Finally, we apply these methods to experimental case-control data and demonstrate successful identification of published genetic susceptibility loci for a rare monogenic disease (sudden infant death with dysgenesis of the testes syndrome), a rare complex disease (progressive supranuclear palsy), and a common complex disease (Alzheimer disease) across multiple SNP genotyping platforms. On the basis of these theoretical calculations and their experimental validation, our results suggest that pooling-based GWA studies are a logical first step for determining whether major genetic associations exist in diseases with high heritability.     

Towards compendia of negative genetic association studies: an example for Alzheimer disease

Most genetic sequence variants that contribute to variability in complex human traits will have small effects that are not readily detectable with population samples typically used in genetic association studies. A potentially valuable tool in the gene discovery process is meta-analysis of the accumulated published data, but in order to be valid these require a sample of studies representative of the true genetic effect and thus hypothetically should include some positive and an abundance of negative reports. A survey of the literature on association studies for Alzheimer disease (AD) from January 2004–April 2005, identified 138 studies, 86 of which reported positive findings other than for apolipoprotein E (APOE), strongly indicative of publication bias. We report here an analysis of 62 genetic markers, tested for association with AD risk as well as for possible effects upon quantitative indices of AD severity (mini-mental state examination scores, age-at-onset, and cerebrospinal fluid (CSF) β-amyloid (Aβ) and CSF tau proteins). Within this set, only modest signals were present that, with the exception of APOE are easily lost when corrections for multiple hypotheses are applied. In isolation, results are thus broadly negative. Genes studied encompass both novel candidates as well as several recently claimed to be associated with AD (e.g. urokinase plasminogen activator (PLAU) and acetyl-coenzyme A acetyltransferase 1 (ACAT1)). By reporting these data we hope to encourage the publication of gene compendia to guide further studies and aid future meta-analyses aimed at resolving the involvement of genes in complex human traits.     

CARD15 genetic variation in a Quebec population: prevalence,genotype-phenotype relationship,and haplotype structure

The caspase recruitment domain gene (CARD15) was recently identified as the underlying gene associated with the IBD1 locus that confers susceptibility to Crohn disease (CD). CARD15 is related to the NOD1/Apaf-1 family of apoptosis regulators, and three sequence variants (Arg702Trp, Gly908Arg, and Leu1007fsinsC) in the gene were demonstrated to be associated with CD. We collected a cohort of 231 patients with CD and 71 healthy control individuals from the Canadian province of Quebec, to determine the prevalence of these sequence variants in an independent population. Clinical records of all patients were systematically reviewed, and detailed phenotypic information was obtained. All patient DNA samples were genotyped for the three variants, thus enabling an analysis of genotype-phenotype correlations. In this cohort, 45.0% of patients with CD carried at least one variant in the CARD15 gene, compared with 9.0% of control individuals (P<10-7). Allele frequencies of Arg702Trp, Gly908Arg, and Leu1007fsinsC were 12.9%, 5.2%, and 10.3% in patients with CD, compared with 4.2%, 0.7%, and 0.7% in control individuals, respectively. Importantly, CARD15 mutants were seen with equal frequency in patients with familial and sporadic CD. Analysis of the relationship between genotype and phenotype convincingly demonstrates that CARD15 variants are significantly associated with ileal disease involvement, as opposed to strictly colonic disease (P<.001). Moreover, we were able to determine the haplotype structure surrounding this disease gene by genotyping 45 single-nucleotide polymorphisms (SNPs) in a 177-kb region that contained the CARD15 gene. This structure helps clarify the history of these causal mutations. Finally, this analysis shows that CARD15 involvement with CD is detectable by use of publicly available SNPs alone.     

Meta-analysis on association between the ATP-binding cassette transporter A1 gene (ABCA1) and Alzheimer's disease

Purpose

In the past decade, a number of case–control studies have been carried out to investigate the relationship between ABCA1 polymorphisms and Alzheimer's disease (AD). However, these studies have yielded contradictory results. To investigate this inconsistency, a meta-analysis was performed.

Methods

Databases including PubMed, Web of Science, EMBASE and CNKI were searched to find relevant studies. Odds ratios (ORs) with 95% confidence intervals (CIs) were used to assess the strength of association.

Results

A total of 13 case–control studies, involving 6214 patients and 6034 controls for ABCA1 polymorphisms were included. In a combined analysis, the summary per-allele odds ratio for AD of the 219 K was 1.03 (95% CI: 0.93–1.14, p = 0.56). A meta-analysis of studies on the 883 M and 1587 K variant showed no significant overall association with AD, yielding a per-allele odds ratio of 1.10 (95% CI: 0.96–1.26, p = 0.16), and 1.09 (95% CI: 0.97–1.24, p = 0.16) respectively. Similar results were also found for heterozygous and homozygous. In the subgroup analysis by ethnicity, sample size, APOE status and onset type, no significant associations were found in almost all genetic models.

Conclusions

In summary, there was no significant association detected between ABCA1 R219K, I883M and R1587K polymorphisms and risk for AD.     

Cladistic structure within the human Lipoprotein lipase gene and its implications for phenotypic association studies

Haplotype variation in 9.7 kb of genomic DNA sequence from the human lipoprotein lipase (LPL) gene was scored in three populations: African-Americans from Jackson, Mississippi (24 individuals), Finns from North Karelia, Finland (24), and non-Hispanic whites from Rochester, Minnesota (23). Earlier analyses had indicated that recombination was common but concentrated into a hotspot and that recurrent mutations at multiple sites may have occurred. We show that much evolutionary structure exists in the haplotype variation on either side of the recombinational hotspot. By peeling off significant recombination events from a tree estimated under the null hypothesis of no recombination, we also reveal some cladistic structure not disrupted by recombination during the time to coalescence of this variation. Additional cladistic structure is estimated to have emerged after recombination. Many apparent multiple mutational events at sites still remain after removing the effects of the detected recombination/gene conversion events. These apparent multiple events are found primarily at sites identified as highly mutable by previous studies, strengthening the conclusion that they are true multiple events. This analysis portrays the complexity of the interplay among many recombinational and mutational events that would be needed to explain the patterns of haplotype diversity in this gene. The cladistic structure in this region is used to identify four to six single-nucleotide polymorphisms (SNPs) that would provide disequilibrium coverage over much of this region. These sites may be useful in identifying phenotypic associations with variable sites in this gene. Evolutionary considerations also imply that the SNPs in the 3' region should have general utility in most human populations, but the 5' SNPs may be more population specific. Choosing SNPs at random would generally not provide adequate disequilibrium coverage of the sequenced region.