Genomewide association study for susceptibility genes contributing to familial Parkinson disease

Five genes have been identified that contribute to Mendelian forms of Parkinson disease (PD); however, mutations have been found in fewer than 5% of patients, suggesting that additional genes contribute to disease risk. Unlike previous studies that focused primarily on sporadic PD, we have performed the first genomewide association study (GWAS) in familial PD. Genotyping was performed with the Illumina HumanCNV370Duo array in 857 familial PD cases and 867 controls. A logistic model was employed to test for association under additive and recessive modes of inheritance after adjusting for gender and age. No result met genomewide significance based on a conservative Bonferroni correction. The strongest association result was with SNPs in the GAK/DGKQ region on chromosome 4 (additive model: p = 3.4 × 10⁻⁶; OR = 1.69). Consistent evidence of association was also observed to the chromosomal regions containing SNCA (additive model: p = 5.5 × 10⁻⁵; OR = 1.35) and MAPT (recessive model: p = 2.0 × 10⁻⁵; OR = 0.56). Both of these genes have been implicated previously in PD susceptibility; however, neither was identified in previous GWAS studies of PD. Meta-analysis was performed using data from a previous case–control GWAS, and yielded improved p values for several regions, including GAK/DGKQ (additive model: p = 2.5 × 10⁻⁷) and the MAPT region (recessive model: p = 9.8 × 10⁻⁶; additive model: p = 4.8 × 10⁻⁵). These data suggest the identification of new susceptibility alleles for PD in the GAK/DGKQ region, and also provide further support for the role of SNCA and MAPT in PD susceptibility. Electronic supplementary material  The online version of this article (doi:) contains supplementary material, which is available to authorized users. N. Pankratz and J. B. Wilk are joint first authors.     

A whole-genome association study for pig reproductive traits

A whole-genome association study was performed for reproductive traits in commercial sows using the PorcineSNP60 BeadChip and Bayesian statistical methods. The traits included total number born (TNB), number born alive (NBA), number of stillborn (SB), number of mummified foetuses at birth (MUM) and gestation length (GL) in each of the first three parities. We report the associations of informative QTL and the genes within the QTL for each reproductive trait in different parities. These results provide evidence of gene effects having temporal impacts on reproductive traits in different parities. Many QTL identified in this study are new for pig reproductive traits. Around 48% of total genes located in the identified QTL regions were predicted to be involved in placental functions. The genomic regions containing genes important for foetal developmental (e.g. MEF2C) and uterine functions (e.g. PLSCR4) were associated with TNB and NBA in the first two parities. Similarly, QTL in other foetal developmental (e.g. HNRNPD and AHR) and placental (e.g. RELL1 and CD96) genes were associated with SB and MUM in different parities. The QTL with genes related to utero-placental blood flow (e.g. VEGFA) and hematopoiesis (e.g. MAFB) were associated with GL differences among sows in this population. Pathway analyses using genes within QTL identified some modest underlying biological pathways, which are interesting candidates (e.g. the nucleotide metabolism pathway for SB) for pig reproductive traits in different parities. Further validation studies on large populations are warranted to improve our understanding of the complex genetic architecture for pig reproductive traits.     

Replication of association between ELAVL4 and Parkinson disease: the GenePD study

Genetic variants in embryonic lethal, abnormal vision, Drosophila-like 4 (ELAVL4) have been reported to be associated with onset age of Parkinson disease (PD) or risk for PD affection in Caucasian populations. In the current study we genotyped three single nucleotide polymorphisms in ELAVL4 in a Caucasian study sample consisting of 712 PD patients and 312 unrelated controls from the GenePD study. The minor allele of rs967582 was associated with increased risk of PD (odds ratio = 1.46, nominal P value = 0.011) in the GenePD population. The minor allele of rs967582 was also the risk allele for PD affection or earlier onset age in the previously studied populations. This replication of association with rs967582 in a third cohort further implicates ELAVL4 as a PD susceptibility gene.     

A validated whole-genome association study of efficient food conversion in cattle

The genetic factors that contribute to efficient food conversion are largely unknown. Several physiological systems are likely to be important, including basal metabolic rate, the generation of ATP, the regulation of growth and development, and the homeostatic control of body mass. Using whole-genome association, we found that DNA variants in or near proteins contributing to the background use of energy of the cell were 10 times as common as those affecting appetite and body-mass homeostasis. In addition, there was a genic contribution from the extracellular matrix and tissue structure, suggesting a trade-off between efficiency and tissue construction. Nevertheless, the largest group consisted of those involved in gene regulation or control of the phenotype. We found that the distribution of micro-RNA motifs was significantly different for the genetic variants associated with residual feed intake than for the genetic variants in total, although the distribution of promoter sequence motifs was not different. This suggests that certain subsets of micro-RNA are more important for the regulation of this trait. Successful validation depended on the sign of the allelic association in different populations rather than on the strength of the initial association or its size of effect.     

Polymorphism in environment responsive genes and association with Parkinson disease

Attempts were made in the present case-control study to investigate the association of polymorphism in the genes encoding proteins involved in toxication–detoxication and dopaminergic pathways and susceptibility to Parkinson’s disease (PD). Seventy patients suffering from PD and one hundred healthy controls belonging to the same geographical location and same ethnicity were included in the study. PCR-RFLP and allele-specific PCR-based methodology were used to identify the genotypes. Multivariate logistic regression analysis revealed that heterozygous genotypes of cytochrome P4502D6*4(CYP2D6*4), CYP2E1*5B (RsaI) polymorphism and homozygous mutant genotypes of CYP2E1*6 (Dra1) were found to be overrepresented in PD cases when compared to the controls. Risk was also found to be increased in patients carrying glutathione S-transferase T1 (GSTT1) null or homozygous variant genotypes of GSTP1. Significant association was observed for monoamine oxidase-B(MAO-B) variant allele G and PD, whereas no difference in genotype and allele frequencies was observed for manganese-superoxide dismutase (MnSOD), dopamine receptor-D2(DRD2), and dopamine transporter (DAT) genes between controls and PD cases. Genotype combinations characterized by the presence of two variant genotypes on their corresponding loci revealed that four combinations of GSTT1 null and MnSOD(-9Val) or GST null and MAOB-G or CYP2E1*5B and MAO-B-AG or CYP2E1*5B and DRD2 (Taq1A-het) genotypes in the patients exhibited severalfold higher and significant association with risk to PD. Our data suggest that polymorphism in the genes involved in detoxification and dopamine regulation may modulate the susceptibility to PD and could be important risk factors in the pathogenesis of PD.     

Linkage disequilibrium and association of MAPT H1 in Parkinson disease

The MAPT H1 haplotype has been associated with four-repeat (4R) tauopathies, including progressive supranuclear palsy, corticobasal degeneration, and argyrophilic grain disease. More controversial is that the same haplotype has been associated with Parkinson disease (PD). Using H1-specific single-nucleotide polymorphisms, we demonstrate that MAPT H1 is a misnomer and consists of a family of recombining H1 alleles. Population genetics, linkage disequilibrium, and association analyses have shown that specific MAPT H1 subhaplotypes are preferentially associated with Parkinson disease. Using a sliding scale of MAPT H1-specific haplotypes--in age/sex-matched PD cases and controls from central Norway--we have refined the disease association to within an approximately 90-kb interval of the 5' end of the MAPT locus.     

ParaHaplo: A program package for haplotype-based whole-genome association study using parallel computing

Background  

Fungi from environmental samples are typically identified to species level through DNA sequencing of the nuclear ribosomal internal transcribed spacer (ITS) region for use in BLAST-based similarity searches in the International Nucleotide Sequence Databases. These searches are time-consuming and regularly require a significant amount of manual intervention and complementary analyses. We here present software – in the form of an identification pipeline for large sets of fungal ITS sequences – developed to automate the BLAST process and several additional analysis steps. The performance of the pipeline was evaluated on a dataset of 350 ITS sequences from fungi growing as epiphytes on building material.     

Type 2 diabetes whole-genome association study in four populations: the DiaGen consortium

Type 2 diabetes (T2D) is a common, polygenic chronic disease with high heritability. The purpose of this whole-genome association study was to discover novel T2D-associated genes. We genotyped 500 familial cases and 497 controls with >300,000 HapMap-derived tagging single-nucleotide-polymorphism (SNP) markers. When a stringent statistical correction for multiple testing was used, the only significant SNP was at TCF7L2, which has already been discovered and confirmed as a T2D-susceptibility gene. For a replication study, we selected 10 SNPs in six chromosomal regions with the strongest association (singly or as part of a haplotype) for retesting in an independent case-control set including 2,573 T2D cases and 2,776 controls. The most significant replicated result was found at the AHI1-LOC441171 gene region.     

A whole-genome association approach for large-scale interspecies traits

Science China Life Sciences -     

Fine mapping versus replication in whole-genome association studies

Association replication studies have a poor track record and, even when successful, often claim association with different markers, alleles, and phenotypes than those reported in the primary study. It is unknown whether these outcomes reflect genuine associations or false-positive results. A greater understanding of these observations is essential for genomewide association (GWA) studies, since they have the potential to identify multiple new associations that that will require external validation. Theoretically, a repeat association with precisely the same variant in an independent sample is the gold standard for replication, but testing additional variants is commonplace in replication studies. Finding different associated SNPs within the same gene or region as that originally identified is often reported as confirmatory evidence. Here, we compare the probability of replicating a gene or region under two commonly used marker-selection strategies: an "exact" approach that involves only the originally significant markers and a "local" approach that involves both the originally significant markers and others in the same region. When a region of high intermarker linkage disequilibrium is tested to replicate an initial finding that is only weak association with disease, the local approach is a good strategy. Otherwise, the most powerful and efficient strategy for replication involves testing only the initially identified variants. Association with a marker other than that originally identified can occur frequently, even in the presence of real effects in a low-powered replication study, and instances of such association increase as the number of included variants increases. Our results provide a basis for the design and interpretation of GWA replication studies and point to the importance of a clear distinction between fine mapping and replication after GWA.     

An R package for analysis of whole-genome association studies

OBJECTIVE: To provide data classes and methods to facilitate the analysis of whole genome association studies in the R language for statistical computing. METHODS: We have implemented data classes in which each genotype call is stored as a single byte. At this density, data for single chromosomes derived from large studies and new high-throughput gene chip platforms can be handled in memory. We use the object-oriented programming model introduced with version 4 of the S-plus package, usually termed 'S4 methods'. RESULTS: At the current state of development the package only supports population-based studies, although we would hope to provide support for family-based studies soon. Both quantitative and qualitative phenotypes may be analysed. Flexible association testing functions are provided which can carry out single SNP tests which control for potential confounding by quantitative and qualitative covariates. Tests involving several SNPs taken together as 'tags' are also supported. Efficient calculation of pair-wise linkage disequilibrium measures is implemented and data input functions include a function which can download data directly from the international HapMap project website.     

Calibrating the performance of SNP arrays for whole-genome association studies

To facilitate whole-genome association studies (WGAS), several high-density SNP genotyping arrays have been developed. Genetic coverage and statistical power are the primary benchmark metrics in evaluating the performance of SNP arrays. Ideally, such evaluations would be done on a SNP set and a cohort of individuals that are both independently sampled from the original SNPs and individuals used in developing the arrays. Without utilization of an independent test set, previous estimates of genetic coverage and statistical power may be subject to an overfitting bias. Additionally, the SNP arrays' statistical power in WGAS has not been systematically assessed on real traits. One robust setting for doing so is to evaluate statistical power on thousands of traits measured from a single set of individuals. In this study, 359 newly sampled Americans of European descent were genotyped using both Affymetrix 500K (Affx500K) and Illumina 650Y (Ilmn650K) SNP arrays. From these data, we were able to obtain estimates of genetic coverage, which are robust to overfitting, by constructing an independent test set from among these genotypes and individuals. Furthermore, we collected liver tissue RNA from the participants and profiled these samples on a comprehensive gene expression microarray. The RNA levels were used as a large-scale set of quantitative traits to calibrate the relative statistical power of the commercial arrays. Our genetic coverage estimates are lower than previous reports, providing evidence that previous estimates may be inflated due to overfitting. The Ilmn650K platform showed reasonable power (50% or greater) to detect SNPs associated with quantitative traits when the signal-to-noise ratio (SNR) is greater than or equal to 0.5 and the causal SNP's minor allele frequency (MAF) is greater than or equal to 20% (N=359). In testing each of the more than 40,000 gene expression traits for association to each of the SNPs on the Ilmn650K and Affx500K arrays, we found that the Ilmn650K yielded 15% times more discoveries than the Affx500K at the same false discovery rate (FDR) level.     

ParaHaplo 2.0: a program package for haplotype-estimation and haplotype-based whole-genome association study using parallel computing

Background  

The use of haplotype-based association tests can improve the power of genome-wide association studies. Since the observed genotypes are unordered pairs of alleles, haplotype phase must be inferred. However, estimating haplotype phase is time consuming. When millions of single-nucleotide polymorphisms (SNPs) are analyzed in genome-wide association study, faster methods for haplotype estimation are required.