Analysis Of Genome-wide Association Data Research Articles

A simple yet efficient method of local false discovery rate estimation designed for genome-wide association data analysis

A simple yet efficient method of local false discovery rate estimation designed for genome-wide association data analysis

Fifteen Years of Gene Set Analysis for High-Throughput Genomic Data: A Review of Statistical Approaches and Future Challenges.

Over the last decade, gene set analysis has become the first choice for gaining insights into underlying complex biology of diseases through gene expression and gene association studies. It also reduces the complexity of statistical analysis and enhances the explanatory power of the obtained results. Although gene set analysis approaches are extensively used in gene expression and genome wide association data analysis, the statistical structure and steps common to these approaches have not yet been comprehensively discussed, which limits their utility. In this article, we provide a comprehensive overview, statistical structure and steps of gene set analysis approaches used for microarrays, RNA-sequencing and genome wide association data analysis. Further, we also classify the gene set analysis approaches and tools by the type of genomic study, null hypothesis, sampling model and nature of the test statistic, etc. Rather than reviewing the gene set analysis approaches individually, we provide the generation-wise evolution of such approaches for microarrays, RNA-sequencing and genome wide association studies and discuss their relative merits and limitations. Here, we identify the key biological and statistical challenges in current gene set analysis, which will be addressed by statisticians and biologists collectively in order to develop the next generation of gene set analysis approaches. Further, this study will serve as a catalog and provide guidelines to genome researchers and experimental biologists for choosing the proper gene set analysis approach based on several factors.

Multi-resolution localization of causal variants across the genome

In the statistical analysis of genome-wide association data, it is challenging to precisely localize the variants that affect complex traits, due to linkage disequilibrium, and to maximize power while limiting spurious findings. Here we report on KnockoffZoom: a flexible method that localizes causal variants at multiple resolutions by testing the conditional associations of genetic segments of decreasing width, while provably controlling the false discovery rate. Our method utilizes artificial genotypes as negative controls and is equally valid for quantitative and binary phenotypes, without requiring any assumptions about their genetic architectures. Instead, we rely on well-established genetic models of linkage disequilibrium. We demonstrate that our method can detect more associations than mixed effects models and achieve fine-mapping precision, at comparable computational cost. Lastly, we apply KnockoffZoom to data from 350k subjects in the UK Biobank and report many new findings.

Agnostic Pathway/Gene Set Analysis of Genome-Wide Association Data Identifies Associations for Pancreatic Cancer.

Genome-wide association studies (GWAS) identify associations of individual single-nucleotide polymorphisms (SNPs) with cancer risk but usually only explain a fraction of the inherited variability. Pathway analysis of genetic variants is a powerful tool to identify networks of susceptibility genes. We conducted a large agnostic pathway-based meta-analysis of GWAS data using the summary-based adaptive rank truncated product method to identify gene sets and pathways associated with pancreatic ductal adenocarcinoma (PDAC) in 9040 cases and 12 496 controls. We performed expression quantitative trait loci (eQTL) analysis and functional annotation of the top SNPs in genes contributing to the top associated pathways and gene sets. All statistical tests were two-sided. We identified 14 pathways and gene sets associated with PDAC at a false discovery rate of less than 0.05. After Bonferroni correction (P ≤ 1.3 × 10-5), the strongest associations were detected in five pathways and gene sets, including maturity-onset diabetes of the young, regulation of beta-cell development, role of epidermal growth factor (EGF) receptor transactivation by G protein-coupled receptors in cardiac hypertrophy pathways, and the Nikolsky breast cancer chr17q11-q21 amplicon and Pujana ATM Pearson correlation coefficient (PCC) network gene sets. We identified and validated rs876493 and three correlating SNPs (PGAP3) and rs3124737 (CASP7) from the Pujana ATM PCC gene set as eQTLs in two normal derived pancreas tissue datasets. Our agnostic pathway and gene set analysis integrated with functional annotation and eQTL analysis provides insight into genes and pathways that may be biologically relevant for risk of PDAC, including those not previously identified.

Pathway-Based Analysis of Genome-Wide Association Data Identified SNPs in HMMR as Biomarker for Chemotherapy- Induced Neutropenia in Breast Cancer Patients.

Neutropenia secondary to chemotherapy in breast cancer patients can be life-threatening and there are no biomarkers available to predict the risk of drug-induced neutropenia in those patients. We previously performed a genome-wide association study (GWAS) for neutropenia events in women with breast cancer who were treated with 5-fluorouracil, epirubicin and cyclophosphamide and recruited to the SUCCESS-A trial. A genome-wide significant single-nucleotide polymorphism (SNP) signal in the tumor necrosis factor superfamily member 13B (TNFSF13B) gene, encoding the cytokine B-cell activating factor (BAFF), was identified in that GWAS. Taking advantage of these existing GWAS data, in the present study we utilized a pathway-based analysis approach by leveraging knowledge of the pharmacokinetics and pharmacodynamics of drugs and breast cancer pathophysiology to identify additional SNPs/genes associated with the underlying etiology of chemotherapy-induced neutropenia. We identified three SNPs in the hyaluronan mediated motility receptor (HMMR) gene that were significantly associated with neutropenia (p < 1.0E-04). Those three SNPs were trans-expression quantitative trait loci for the expression of TNFSF13B (p < 1.0E-04). The minor allele of these HMMR SNPs was associated with a decreased TNFSF13B mRNA level. Additional functional studies performed with lymphoblastoid cell lines (LCLs) demonstrated that LCLs possessing the minor allele for the HMMR SNPs were more sensitive to drug treatment. Knock-down of TNFSF13B in LCLs and HL-60 promyelocytic cells and treatment of those cells with BAFF modulated the cell sensitivity to chemotherapy treatment. These results demonstrate that HMMR SNP-dependent cytotoxicity of these chemotherapeutic agents might be related to TNFSF13B expression level. In summary, utilizing a pathway-based approach for the analysis of GWAS data, we identified additional SNPs in the HMMR gene that were associated with neutropenia and also were correlated with TNFSF13B expression.

Analysis of genome-wide association data highlights candidates for drug repositioning in psychiatry.

Knowledge of psychiatric disease genetics has advanced rapidly during the past decade with the advent of genome-wide association studies (GWAS). However, less progress has been made in harnessing these data to reveal new therapies. Here we propose a framework for drug repositioning by comparing transcriptomes imputed from GWAS data with drug-induced gene expression profiles from the Connectivity Map database and apply this approach to seven psychiatric disorders. We found a number of repositioning candidates, many supported by preclinical or clinical evidence. Repositioning candidates for a number of disorders were also significantly enriched for known psychiatric medications or therapies considered in clinical trials. For example, candidates for schizophrenia were enriched for antipsychotics, while those for bipolar disorder were enriched for both antipsychotics and antidepressants. These findings provide support for the usefulness of GWAS data in guiding drug discovery.

Deletions at 22q11.2 in idiopathic Parkinson's disease: a combined analysis of genome-wide association data

SummaryBackgroundParkinson's disease has been reported in a small number of patients with chromosome 22q11.2 deletion syndrome. In this study, we screened a series of large, independent Parkinson's disease case-control studies for deletions at 22q11.2.MethodsWe used data on deletions spanning the 22q11.2 locus from four independent case-control Parkinson's disease studies (UK Wellcome Trust Case Control Consortium 2, Dutch Parkinson's Disease Genetics Consortium, US National Institute on Aging, and International Parkinson's Disease Genomics Consortium studies), which were independent of the original reports of chromosome 22q11.2 deletion syndrome. We did case-control association analysis to compare the proportion of 22q11.2 deletions found, using the Fisher's exact test for the independent case-control studies and the Mantel-Haenszel test for the meta-analyses. We retrieved clinical details of patients with Parkinson's disease who had 22q11.2 deletions from the medical records of these patients.FindingsWe included array-based copy number variation data from 9387 patients with Parkinson's disease and 13 863 controls. Eight patients with Parkinson's disease and none of the controls had 22q11.2 deletions (p=0·00082). In the 8451 patients for whom age at onset data were available, deletions at 22q11.2 were associated with Parkinson's disease age at onset (Mann-Whitney U test p=0·001). Age at onset of Parkinson's disease was lower in patients carrying a 22q11.2 deletion (median 37 years, 95% CI 32·0–55·5; mean 42·1 years [SD 11·9]) than in those who did not carry a deletion (median 61 years, 95% CI 60·5–61·0; mean 60·3 years [SD 12·8]). A 22q11.2 deletion was present in more patients with early-onset (p<0·0001) and late-onset Parkinson's disease (p=0·016) than in controls, and in more patients with early-onset than late-onset Parkinson's disease (p=0·005).InterpretationClinicians should be alert to the possibility of 22q11.2 deletions in patients with Parkinson's disease who have early presentation or features associated with the chromosome 22q11.2 deletion syndrome, or both.FundingUK Medical Research Council, UK Wellcome Trust, Parkinson's UK, Patrick Berthoud Trust, National Institutes of Health, “Investissements d'Avenir” ANR-10-IAIHU-06, Dutch Parkinson Foundation (Parkinson Vereniging), Neuroscience Campus Amsterdam, National Institute for Health Research, National Institute on Aging, National Institutes of Health.

Genome-Wide Association mapping for intelligence in military working dogs: unsupervised classification algorithm for genome-wide association data analysis

Genome-Wide Association mapping for intelligence in military working dogs: unsupervised classification algorithm for genome-wide association data analysis

MOBAS: identification of disease-associated protein subnetworks using modularity-based scoring.

Network-based analyses are commonly used as powerful tools to interpret the findings of genome-wide association studies (GWAS) in a functional context. In particular, identification of disease-associated functional modules, i.e., highly connected protein-protein interaction (PPI) subnetworks with high aggregate disease association, are shown to be promising in uncovering the functional relationships among genes and proteins associated with diseases. An important issue in this regard is the scoring of subnetworks by integrating two quantities: disease association of individual gene products and network connectivity among proteins. Current scoring schemes either disregard the level of connectivity and focus on the aggregate disease association of connected proteins or use a linear combination of these two quantities. However, such scoring schemes may produce arbitrarily large subnetworks which are often not statistically significant or require tuning of parameters that are used to weigh the contributions of network connectivity and disease association.Here, we propose a parameter-free scoring scheme that aims to score subnetworks by assessing the disease association of interactions between pairs of gene products. We also incorporate the statistical significance of network connectivity and disease association into the scoring function. We test the proposed scoring scheme on a GWAS dataset for two complex diseases type II diabetes (T2D) and psoriasis (PS). Our results suggest that subnetworks identified by commonly used methods may fail tests of statistical significance after correction for multiple hypothesis testing. In contrast, the proposed scoring scheme yields highly significant subnetworks, which contain biologically relevant proteins that cannot be identified by analysis of genome-wide association data alone. We also show that the proposed scoring scheme identifies subnetworks that are reproducible across different cohorts, and it can robustly recover relevant subnetworks at lower sampling rates.

Genetic analysis for a shared biological basis between migraine and coronary artery disease.

Objective:To apply genetic analysis of genome-wide association data to study the extent and nature of a shared biological basis between migraine and coronary artery disease (CAD).Methods:Four separate methods for cross-phenotype genetic analysis were applied on data from 2 large-scale genome-wide association studies of migraine (19,981 cases, 56,667 controls) and CAD (21,076 cases, 63,014 controls). The first 2 methods quantified the extent of overlapping risk variants and assessed the load of CAD risk loci in migraineurs. Genomic regions of shared risk were then identified by analysis of covariance patterns between the 2 phenotypes and by querying known genome-wide significant loci.Results:We found a significant overlap of genetic risk loci for migraine and CAD. When stratified by migraine subtype, this was limited to migraine without aura, and the overlap was protective in that patients with migraine had a lower load of CAD risk alleles than controls. Genes indicated by 16 shared risk loci point to mechanisms with potential roles in migraine pathogenesis and CAD, including endothelial dysfunction (PHACTR1) and insulin homeostasis (GIP).Conclusions:The results suggest that shared biological processes contribute to risk of migraine and CAD, but surprisingly this commonality is restricted to migraine without aura and the impact is in opposite directions. Understanding the mechanisms underlying these processes and their opposite relationship to migraine and CAD may improve our understanding of both disorders.

Genome-wide association data classification and SNPs selection using two-stage quality-based Random Forests.

BackgroundSingle-nucleotide polymorphisms (SNPs) selection and identification are the most important tasks in Genome-wide association data analysis. The problem is difficult because genome-wide association data is very high dimensional and a large portion of SNPs in the data is irrelevant to the disease. Advanced machine learning methods have been successfully used in Genome-wide association studies (GWAS) for identification of genetic variants that have relatively big effects in some common, complex diseases. Among them, the most successful one is Random Forests (RF). Despite of performing well in terms of prediction accuracy in some data sets with moderate size, RF still suffers from working in GWAS for selecting informative SNPs and building accurate prediction models. In this paper, we propose to use a new two-stage quality-based sampling method in random forests, named ts-RF, for SNP subspace selection for GWAS. The method first applies p-value assessment to find a cut-off point that separates informative and irrelevant SNPs in two groups. The informative SNPs group is further divided into two sub-groups: highly informative and weak informative SNPs. When sampling the SNP subspace for building trees for the forest, only those SNPs from the two sub-groups are taken into account. The feature subspaces always contain highly informative SNPs when used to split a node at a tree.ResultsThis approach enables one to generate more accurate trees with a lower prediction error, meanwhile possibly avoiding overfitting. It allows one to detect interactions of multiple SNPs with the diseases, and to reduce the dimensionality and the amount of Genome-wide association data needed for learning the RF model. Extensive experiments on two genome-wide SNP data sets (Parkinson case-control data comprised of 408,803 SNPs and Alzheimer case-control data comprised of 380,157 SNPs) and 10 gene data sets have demonstrated that the proposed model significantly reduced prediction errors and outperformed most existing the-state-of-the-art random forests. The top 25 SNPs in Parkinson data set were identified by the proposed model including four interesting genes associated with neurological disorders.ConclusionThe presented approach has shown to be effective in selecting informative sub-groups of SNPs potentially associated with diseases that traditional statistical approaches might fail. The new RF works well for the data where the number of case-control objects is much smaller than the number of SNPs, which is a typical problem in gene data and GWAS. Experiment results demonstrated the effectiveness of the proposed RF model that outperformed the state-of-the-art RFs, including Breiman's RF, GRRF and wsRF methods.

Two novel pathway analysis methods based on a hierarchical model

Motivation: Over the past few years several pathway analysis methods have been proposed for exploring and enhancing the analysis of genome-wide association data. Hierarchical models have been advocated as a way to integrate SNP and pathway effects in the same model, but their computational complexity has prevented them being applied on a genome-wide scale to date.Methods: We present two novel methods for identifying associated pathways. In the proposed hierarchical model, the SNP effects are analytically integrated out of the analysis, allowing computationally tractable model fitting to genome-wide data. The first method uses Bayes factors for calculating the effect of the pathways, whereas the second method uses a machine learning algorithm and adaptive lasso for finding a sparse solution of associated pathways.Results: The performance of the proposed methods was explored on both simulated and real data. The results of the simulation study showed that the methods outperformed some well-established association methods: the commonly used Fisher’s method for combining P-values and also the recently published BGSA. The methods were applied to two genome-wide association study datasets that aimed to find the genetic structure of platelet function and body mass index, respectively. The results of the analyses replicated the results of previously published pathway analysis of these phenotypes but also identified novel pathways that are potentially involved.Availability: An R package is under preparation. In the meantime, the scripts of the methods are available on request from the authors.Contact: marina.evangelou@cimr.cam.ac.ukSupplementary Information: Supplementary data are available at Bioinformatics online.

A “Candidate-Interactome” Aggregate Analysis of Genome-Wide Association Data in Multiple Sclerosis

Though difficult, the study of gene-environment interactions in multifactorial diseases is crucial for interpreting the relevance of non-heritable factors and prevents from overlooking genetic associations with small but measurable effects. We propose a “candidate interactome” (i.e. a group of genes whose products are known to physically interact with environmental factors that may be relevant for disease pathogenesis) analysis of genome-wide association data in multiple sclerosis. We looked for statistical enrichment of associations among interactomes that, at the current state of knowledge, may be representative of gene-environment interactions of potential, uncertain or unlikely relevance for multiple sclerosis pathogenesis: Epstein-Barr virus, human immunodeficiency virus, hepatitis B virus, hepatitis C virus, cytomegalovirus, HHV8-Kaposi sarcoma, H1N1-influenza, JC virus, human innate immunity interactome for type I interferon, autoimmune regulator, vitamin D receptor, aryl hydrocarbon receptor and a panel of proteins targeted by 70 innate immune-modulating viral open reading frames from 30 viral species. Interactomes were either obtained from the literature or were manually curated. The P values of all single nucleotide polymorphism mapping to a given interactome were obtained from the last genome-wide association study of the International Multiple Sclerosis Genetics Consortium & the Wellcome Trust Case Control Consortium, 2. The interaction between genotype and Epstein Barr virus emerges as relevant for multiple sclerosis etiology. However, in line with recent data on the coexistence of common and unique strategies used by viruses to perturb the human molecular system, also other viruses have a similar potential, though probably less relevant in epidemiological terms.

A PLSPM-Based Test Statistic for Detecting Gene-Gene Co-Association in Genome-Wide Association Study with Case-Control Design

For genome-wide association data analysis, two genes in any pathway, two SNPs in the two linked gene regions respectively or in the two linked exons respectively within one gene are often correlated with each other. We therefore proposed the concept of gene-gene co-association, which refers to the effects not only due to the traditional interaction under nearly independent condition but the correlation between two genes. Furthermore, we constructed a novel statistic for detecting gene-gene co-association based on Partial Least Squares Path Modeling (PLSPM). Through simulation, the relationship between traditional interaction and co-association was highlighted under three different types of co-association. Both simulation and real data analysis demonstrated that the proposed PLSPM-based statistic has better performance than single SNP-based logistic model, PCA-based logistic model, and other gene-based methods.

A bias-reducing pathway enrichment analysis of genome-wide association data confirmed association of the MHC region with schizophrenia

BackgroundAfter the recent successes of genome-wide association studies (GWAS), one key challenge is to identify genetic variants that might have a significant joint effect on complex diseases but have failed...

Melatonin pathway genes and breast cancer risk among Chinese women

Previous studies suggest that melatonin may act on cancer growth through a variety of mechanisms, most notably by direct anti-proliferative effects on breast cancer cells and via interactions with the estrogen pathway. Three genes are largely responsible for mediating the downstream effects of melatonin: melatonin receptors 1a and 1b (MTNR1a and MTNR1b), and arylalkylamine N-acetyltransferase (AANAT). It is hypothesized that genetic variation in these genes may lead to altered protein production or function. To address this question, we conducted a comprehensive evaluation of the association between common single nucleotide polymorphisms (SNPs) in the MTNR1a, MTNR1b, and AANAT genes and breast cancer risk among 2,073 cases and 2,083 controls, using a two-stage analysis of genome-wide association data among women of the Shanghai Breast Cancer Study. Results demonstrate two SNPs were consistently associated with breast cancer risk across both study stages. Compared with MTNR1b rs10765576 major allele carriers (GG or GA), a decreased risk of breast cancer was associated with the AA genotype (OR = 0.78, 95% CI = 0.62-0.97, P = 0.0281). Although no overall association was seen in the combined analysis, the effect of MTNR1a rs7665392 was found to vary by menopausal status (P-value for interaction = 0.001). Premenopausal women with the GG genotype were at increased risk for breast cancer compared with major allele carriers (TT or TG) (OR = 1.57, 95% CI = 1.07-2.31, P = 0.020), while postmenopausal women were at decreased risk (OR = 0.58, 95% 0.36-0.95, P = 0.030). No significant breast cancer associations were found for variants in the AANAT gene. These results suggest that common genetic variation in the MTNR1a and 1b genes may contribute to breast cancer susceptibility, and that associations may vary by menopausal status. Given that multiple variants in high linkage disequilibrium with MTNR1b rs76653292 have been associated with altered function or expression of insulin and glucose family members, further research may focus on clarifying this relationship.

MI-GWAS: a SAS platform for the analysis of inherited and maternal genetic effects in genome-wide association studies using log-linear models

BackgroundSeveral platforms for the analysis of genome-wide association data are available. However, these platforms focus on the evaluation of the genotype inherited by affected (i.e. case) individuals, whereas for some conditions (e.g. birth defects) the genotype of the mothers of affected individuals may also contribute to risk. For such conditions, it is critical to evaluate associations with both the maternal and the inherited (i.e. case) genotype. When genotype data are available for case-parent triads, a likelihood-based approach using log-linear modeling can be used to assess both the maternal and inherited genotypes. However, available software packages for log-linear analyses are not well suited to the analysis of typical genome-wide association data (e.g. including missing data).ResultsAn integrated platform, Maternal and Inherited Analyses for Genome-wide Association Studies (MI-GWAS) for log-linear analyses of maternal and inherited genetic effects in large, genome-wide datasets, is described. MI-GWAS uses SAS and LEM software in combination to appropriately format data, perform the log-linear analyses and summarize the results. This platform was evaluated using existing genome-wide data and was shown to perform accurately and relatively efficiently.ConclusionsThe MI-GWAS platform provides a valuable tool for the analysis of association of a phenotype or condition with maternal and inherited genotypes using genome-wide data from case-parent triads. The source code for this platform is freely available at http://www.sph.uth.tmc.edu/sbrr/mi-gwas.htm.

Causal graph-based analysis of genome-wide association data in rheumatoid arthritis

BackgroundGWAS owe their popularity to the expectation that they will make a major impact on diagnosis, prognosis and management of disease by uncovering genetics underlying clinical phenotypes. The dominant paradigm in GWAS data analysis so far consists of extensive reliance on methods that emphasize contribution of individual SNPs to statistical association with phenotypes. Multivariate methods, however, can extract more information by considering associations of multiple SNPs simultaneously. Recent advances in other genomics domains pinpoint multivariate causal graph-based inference as a promising principled analysis framework for high-throughput data. Designed to discover biomarkers in the local causal pathway of the phenotype, these methods lead to accurate and highly parsimonious multivariate predictive models. In this paper, we investigate the applicability of causal graph-based method TIE* to analysis of GWAS data. To test the utility of TIE*, we focus on anti-CCP positive rheumatoid arthritis (RA) GWAS datasets, where there is a general consensus in the community about the major genetic determinants of the disease.ResultsApplication of TIE* to the North American Rheumatoid Arthritis Cohort (NARAC) GWAS data results in six SNPs, mostly from the MHC locus. Using these SNPs we develop two predictive models that can classify cases and disease-free controls with an accuracy of 0.81 area under the ROC curve, as verified in independent testing data from the same cohort. The predictive performance of these models generalizes reasonably well to Swedish subjects from the closely related but not identical Epidemiological Investigation of Rheumatoid Arthritis (EIRA) cohort with 0.71-0.78 area under the ROC curve. Moreover, the SNPs identified by the TIE* method render many other previously known SNP associations conditionally independent of the phenotype.ConclusionsOur experiments demonstrate that application of TIE* captures maximum amount of genetic information about RA in the data and recapitulates the major consensus findings about the genetic factors of this disease. In addition, TIE* yields reproducible markers and signatures of RA. This suggests that principled multivariate causal and predictive framework for GWAS analysis empowers the community with a new tool for high-quality and more efficient discovery.ReviewersThis article was reviewed by Prof. Anthony Almudevar, Dr. Eugene V. Koonin, and Prof. Marianthi Markatou.

Use of genome-wide expression data to mine the "Gray Zone" of GWA studies leads to novel candidate obesity genes.

To get beyond the “low-hanging fruits” so far identified by genome-wide association (GWA) studies, new methods must be developed in order to discover the numerous remaining genes that estimates of heritability indicate should be contributing to complex human phenotypes, such as obesity. Here we describe a novel integrative method for complex disease gene identification utilizing both genome-wide transcript profiling of adipose tissue samples and consequent analysis of genome-wide association data generated in large SNP scans. We infer causality of genes with obesity by employing a unique set of monozygotic twin pairs discordant for BMI (n = 13 pairs, age 24–28 years, 15.4 kg mean weight difference) and contrast the transcript profiles with those from a larger sample of non-related adult individuals (N = 77). Using this approach, we were able to identify 27 genes with possibly causal roles in determining the degree of human adiposity. Testing for association of SNP variants in these 27 genes in the population samples of the large ENGAGE consortium (N = 21,000) revealed a significant deviation of P-values from the expected (P = 4×10−4). A total of 13 genes contained SNPs nominally associated with BMI. The top finding was blood coagulation factor F13A1 identified as a novel obesity gene also replicated in a second GWA set of ∼2,000 individuals. This study presents a new approach to utilizing gene expression studies for informing choice of candidate genes for complex human phenotypes, such as obesity.

Genome-wide association analyses of North American Rheumatoid Arthritis Consortium and Framingham Heart Study data utilizing genome-wide linkage results

The power of genome-wide association studies can be improved by incorporating information from previous study findings, for example, results of genome-wide linkage analyses. Weighted false-discovery rate (FDR) control can incorporate genome-wide linkage scan results into the analysis of genome-wide association data by assigning single-nucleotide polymorphism (SNP) specific weights. Stratified FDR control can also be applied by stratifying the SNPs into high and low linkage strata. We applied these two FDR control methods to the data of North American Rheumatoid Arthritis Consortium (NARAC) study and the Framingham Heart Study (FHS), combining both association and linkage analysis results. For the NARAC study, we used linkage results from a previous genome scan of rheumatoid arthritis (RA) phenotype. For the FHS study, we obtained genome-wide linkage scores from the same 550 k SNP data used for the association analyses of three lipids phenotypes (HDL, LDL, TG). We confirmed some genes previously reported for association with RA and lipid phenotypes. Stratified and weighted FDR methods appear to give improved ranks to some of the replicated SNPs for the RA data, suggesting linkage scan results could provide useful information to improve genome-wide association studies.