Comparisons of seven algorithms for pathway analysis using the WTCCC Crohn's Disease dataset

Data-adaptive and pathway-based tests for association studies between somatic mutations and germline variations in human cancers

Cancer is a disease driven by a combination of inherited genetic variants and somatic mutations. Recently available large-scale sequencing data of cancer genomes have provided an unprecedented opportunity to study the interactions between them. However, previous studies on this topic have been limited by simple, low statistical power tests such as Fisher's exact test. In this paper, we design data-adaptive and pathway-based tests based on the score statistic for association studies between somatic mutations and germline variations. Previous research has shown that two single-nucleotide polymorphism (SNP)-set-based association tests, adaptive sum of powered score (aSPU) and data-adaptive pathway-based (aSPUpath) tests, increase the power in genome-wide association studies (GWASs) with a single disease trait in a case-control study. We extend aSPU and aSPUpath to multi-traits, that is, somatic mutations of multiple genes in a cohort study, allowing extensive information aggregation at both SNP and gene levels.p $p$ -values from different parameters assuming varying genetic architecture are combined to yield data-adaptive tests for somatic mutations and germline variations. Extensive simulations show that, in comparison with some commonly used methods, our data-adaptive somatic mutations/germline variations tests can be applied to multiple germline SNPs/genes/pathways, and generally have much higher statistical powers while maintaining the appropriate type I error. The proposed tests are applied to a large-scale real-world International Cancer Genome Consortium whole genome sequencing data set of 2583 subjects, detecting more significant and biologically relevant associations compared with the other existing methods on both gene and pathway levels. Our study has systematically identified the associations between various germline variations and somatic mutations across different cancer types, which potentially provides valuable utility for cancer risk prediction, prognosis, and therapeutics.

Hierarchical association of COPD to principal genetic components of biological systems

Many disease-causing genetic variants converge on common biological functions and pathways. Precisely how to incorporate pathway knowledge in genetic association studies is not yet clear, however. Previous approaches employ a two-step approach, in which a regular association test is first performed to identify variants associated with the disease phenotype, followed by a test for functional enrichment within the genes implicated by those variants. Here we introduce a concise one-step approach, Hierarchical Genetic Analysis (Higana), which directly computes phenotype associations against each function in the large hierarchy of biological functions documented by the Gene Ontology. Using this approach, we identify risk genes and functions for Chronic Obstructive Pulmonary Disease (COPD), highlighting microtubule transport, muscle adaptation, and nicotine receptor signaling pathways. Microtubule transport has not been previously linked to COPD, as it integrates genetic variants spread over numerous genes. All associations validate strongly in a second COPD cohort.

AI applications in functional genomics

We review the current applications of artificial intelligence (AI) in functional genomics. The recent explosion of AI follows the remarkable achievements made possible by "deep learning", along with a burst of "big data" that can meet its hunger. Biology is about to overthrow astronomy as the paradigmatic representative of big data producer. This has been made possible by huge advancements in the field of high throughput technologies, applied to determine how the individual components of a biological system work together to accomplish different processes. The disciplines contributing to this bulk of data are collectively known as functional genomics. They consist in studies of: i) the information contained in the DNA (genomics); ii) the modifications that DNA can reversibly undergo (epigenomics); iii) the RNA transcripts originated by a genome (transcriptomics); iv) the ensemble of chemical modifications decorating different types of RNA transcripts (epitranscriptomics); v) the products of protein-coding transcripts (proteomics); and vi) the small molecules produced from cell metabolism (metabolomics) present in an organism or system at a given time, in physiological or pathological conditions. After reviewing main applications of AI in functional genomics, we discuss important accompanying issues, including ethical, legal and economic issues and the importance of explainability.

Gene- and pathway-level analyses of iCOGS variants highlight novel signaling pathways underlying familial breast cancer susceptibility

Single‐nucleotide polymorphisms (SNPs) in over 180 loci have been associated with breast cancer (BC) through genome‐wide association studies involving mostly unselected population‐based case‐control series. Some of them modify BC risk of women carrying a BRCA1 or BRCA2 (BRCA1/2) mutation and may also explain BC risk variability in BC‐prone families with no BRCA1/2 mutation. Here, we assessed the contribution of SNPs of the iCOGS array in GENESIS consisting of BC cases with no BRCA1/2 mutation and a sister with BC, and population controls. Genotyping data were available for 1281 index cases, 731 sisters with BC, 457 unaffected sisters and 1272 controls. In addition to the standard SNP‐level analysis using index cases and controls, we performed pedigree‐based association tests to capture transmission information in the sibships. We also performed gene‐ and pathway‐level analyses to maximize the power to detect associations with lower‐frequency SNPs or those with modest effect sizes. While SNP‐level analyses identified 18 loci, gene‐level analyses identified 112 genes. Furthermore, 31 Kyoto Encyclopedia of Genes and Genomes and 7 Atlas of Cancer Signaling Network pathways were highlighted (false discovery rate of 5%). Using results from the “index case‐control” analysis, we built pathway‐derived polygenic risk scores (PRS) and assessed their performance in the population‐based CECILE study and in a data set composed of GENESIS‐affected sisters and CECILE controls. Although these PRS had poor predictive value in the general population, they performed better than a PRS built using our SNP‐level findings, and we found that the joint effect of family history and PRS needs to be considered in risk prediction models.

What's new?

Genetic studies have identified more than 180 single‐nucleotide polymorphisms (SNPs) associated with breast cancer susceptibility, but these studies are reaching their limits. Here, the authors evaluated SNPs in the iCOGS genotyping array using a multilevel approach, including single variant, gene, and pathway analyses. They measured the contribution of the SNPs to breast cancer in patients who have a sister with breast cancer but do not carry a BRCA1/2 mutation. They showed that a pathway‐derived polygenic risk score performed poorly in the general population, and that the best predictive model must include family history.

Association Between Polymorphisms in DNA Damage Repair Genes and Radiation Therapy-Induced Early Adverse Skin Reactions in a Breast Cancer Population: A Polygenic Risk Score Approach

PURPOSE

Genetic variations in DNA damage repair (DDR) genes may influence radiation therapy (RT)-induced acute normal tissue toxicity in patients with breast cancer. Identifying an individual or multiple single-nucleotide polymorphisms (SNPs) associated with RT-induced early adverse skin reactions (EASR) is critical for precision medicine in radiation oncology.

METHODS AND MATERIALS

At the completion of RT, EASR was assessed using the Oncology Nursing Society scale (0-6) in 416 patients with breast cancer, and Oncology Nursing Society score ≥4 was considered RT-induced EASR. PLINK set-based tests and subsequent individual SNP association analyses were conducted to identify genes and SNPs associated with EASR among the 53 DDR genes and 1968 SNPs. A weighted polygenic risk score (PRS) model was constructed to ascertain the association between the joint effect of risk alleles and EASR.

RESULTS

The study population consisted of 264 Hispanic whites, 86 blacks or African Americans, 55 non-Hispanic whites, and 11 others. A total of 115 patients (27.6%) developed EASR. Five genes (ATM, CHEK1, ERCC2, RAD51C, and TGFB1) were significantly associated with RT-induced EASR. Nine SNPs within these 5 genes were further identified: ATM rs61915066, CHEK1 rs11220184, RAD51C rs302877, rs405684, TBFB1 rs4803455, rs2241714, and ERCC2 rs60152947, rs10404465, rs1799786. In a multivariable-adjusted PRS model, patients in a higher quartile of PRS were more likely to develop EASR compared with patients in the lowest quartile (OR_{q2 vs.q1} = 1.94, 95% CI, 0.86-4.39; OR_{q3 vs.q1} = 3.46, 95% CI, 1.57-7.63; OR_{q4 vs.q1} = 8.64, 95% CI, 3.92-19.02; and P_trend < .0001).

CONCLUSIONS

We newly identified the associations between 9 SNPs in ATM, CHEK1, RAD51C, TGFB1, and ERCC2 and RT-induced EASR. PRS modeling showed its potential in identifying populations at risk. Multiple SNPs in DDR genes may jointly contribute to interindividual variation in RT-induced EASR. Validation in an independent external cohort is required to determine the clinical significance of these predictive biomarkers.

A Bayesian hierarchical variable selection prior for pathway-based GWAS using summary statistics

While genome-wide association studies (GWASs) have been widely used to uncover associations between diseases and genetic variants, standard SNP-level GWASs often lack the power to identify SNPs that individually have a moderate effect size but jointly contribute to the disease. To overcome this problem, pathway-based GWASs methods have been developed as an alternative strategy that complements SNP-level approaches. We propose a Bayesian method that uses the generalized fused hierarchical structured variable selection prior to identify pathways associated with the disease using SNP-level summary statistics. Our prior has the flexibility to take in pathway structural information so that it can model the gene-level correlation based on prior biological knowledge, an important feature that makes it appealing compared to existing pathway-based methods. Using simulations, we show that our method outperforms competing methods in various scenarios, particularly when we have pathway structural information that involves complex gene-gene interactions. We apply our method to the Wellcome Trust Case Control Consortium Crohn's disease GWAS data, demonstrating its practical application to real data.

Genome-wide association analysis of hippocampal volume identifies enrichment of neurogenesis-related pathways

Adult neurogenesis occurs in the dentate gyrus of the hippocampus during adulthood and contributes to sustaining the hippocampal formation. To investigate whether neurogenesis-related pathways are associated with hippocampal volume, we performed gene-set enrichment analysis using summary statistics from a large-scale genome-wide association study (N = 13,163) of hippocampal volume from the Enhancing Neuro Imaging Genetics through Meta-Analysis (ENIGMA) Consortium and two year hippocampal volume changes from baseline in cognitively normal individuals from Alzheimer's Disease Neuroimaging Initiative Cohort (ADNI). Gene-set enrichment analysis of hippocampal volume identified 44 significantly enriched biological pathways (FDR corrected p-value < 0.05), of which 38 pathways were related to neurogenesis-related processes including neurogenesis, generation of new neurons, neuronal development, and neuronal migration and differentiation. For genes highly represented in the significantly enriched neurogenesis-related pathways, gene-based association analysis identified TESC, ACVR1, MSRB3, and DPP4 as significantly associated with hippocampal volume. Furthermore, co-expression network-based functional analysis of gene expression data in the hippocampal subfields, CA1 and CA3, from 32 normal controls showed that distinct co-expression modules were mostly enriched in neurogenesis related pathways. Our results suggest that neurogenesis-related pathways may be enriched for hippocampal volume and that hippocampal volume may serve as a potential phenotype for the investigation of human adult neurogenesis.

A non-invasive diagnostic model of immunoglobulin A nephropathy and serological markers for evaluating disease severity

BACKGROUND

Immunoglobulin A nephropathy (IgAN) is the most common pathological type of glomerular disease. Kidney biopsy, the gold standard for IgAN diagnosis, has not been routinely applied in hospitals worldwide due to its invasion nature. Thus, we aim to establish a non-invasive diagnostic model and determine markers to evaluate disease severity by analyzing the serological parameters and pathological stages of patients with IgAN.

METHODS

A total of 272 biopsy-diagnosed IgAN inpatients and 518 non-IgA nephropathy inpatients from the Department of Nephrology of Chinese People's Liberation Army General Hospital were recruited for this study. Routine blood examination, blood coagulation testing, immunoglobulin-complement testing, and clinical biochemistry testing were conducted and pathological stages were analyzed according to Lee grading system. The serological parameters and pathological stages were analyzed. The receiver operating characteristic (ROC) analysis was performed to estimate the diagnostic value of the clinical factors. Logistic regression was used to establish the diagnostic model.

RESULTS

There were 15 significantly different serological parameters between the IgAN and non-IgAN groups (all P < 0.05). The ROC analysis was performed to measure the diagnostic value for IgAN of these parameters and the results showed that the area under the ROC curve (AUC) of total protein (TP), total cholesterol (TC), fibrinogen (FIB), D-dimer (D2), immunoglobulin A (IgA), and immunoglobulin G (IgG) were more than 0.70. The AUC of the "TC + FIB + D2 + IgA + age" combination was 0.86, with a sensitivity of 85.98% and a specificity of 73.85%. Pathological grades of I, II, III, IV, and V accounted for 2.21%, 17.65%, 62.50%, 11.76%, and 5.88%, respectively, with grade III being the most prevalent. The levels of urea nitrogen (UN) (13.57 ± 5.95 vs. 6.06 ± 3.63, 5.92 ± 2.97, 5.41 ± 1.73, and 8.41 ± 3.72 mmol/L, respectively) and creatinine (Cr) (292.19 ± 162.21 vs. 80.42 ± 24.75, 103.79 ± 72.72, 96.41 ± 33.79, and 163.04 ± 47.51 μmol/L, respectively) were significantly higher in grade V than in the other grades, and the levels of TP (64.45 ± 7.56, 67.16 ± 6.94, 63.22 ± 8.56, and 61.41 ± 10.86 vs. 37.47 ± 5.6 mg/d, respectively), direct bilirubin (DB) (2.34 ± 1.23, 2.58 ± 1.40, 1.91 ± 0.97, and 1.81 ± 1.44 vs. 0.74 ± 0.57 μmol/L, respectively), and IgA (310.35 ± 103.78, 318.48 ± 107.54, 292.58 ± 81.85, and 323.29 ± 181.67 vs. 227.17 ± 68.12 g/L, respectively) were significantly increased in grades II-V compared with grade I (all P < 0.05).

CONCLUSIONS

The established diagnostic model that combined multiple factors (TC, FIB, D2, IgA, and age) might be used for IgAN non-invasive diagnosis. TP, DB, IgA, Cr, and UN have the potential to be used to evaluate IgAN disease severity.

Powerful gene set analysis in GWAS with the Generalized Berk-Jones statistic

A common complementary strategy in Genome-Wide Association Studies (GWAS) is to perform Gene Set Analysis (GSA), which tests for the association between one phenotype of interest and an entire set of Single Nucleotide Polymorphisms (SNPs) residing in selected genes. While there exist many tools for performing GSA, popular methods often include a number of ad-hoc steps that are difficult to justify statistically, provide complicated interpretations based on permutation inference, and demonstrate poor operating characteristics. Additionally, the lack of gold standard gene set lists can produce misleading results and create difficulties in comparing analyses even across the same phenotype. We introduce the Generalized Berk-Jones (GBJ) statistic for GSA, a permutation-free parametric framework that offers asymptotic power guarantees in certain set-based testing settings. To adjust for confounding introduced by different gene set lists, we further develop a GBJ step-down inference technique that can discriminate between gene sets driven to significance by single genes and those demonstrating group-level effects. We compare GBJ to popular alternatives through simulation and re-analysis of summary statistics from a large breast cancer GWAS, and we show how GBJ can increase power by incorporating information from multiple signals in the same gene. In addition, we illustrate how breast cancer pathway analysis can be confounded by the frequency of FGFR2 in pathway lists. Our approach is further validated on two other datasets of summary statistics generated from GWAS of height and schizophrenia.

Genetic influence on scar height and pliability after burn injury in individuals of European ancestry: A prospective cohort study

After similar extent of injury there is considerable variability in scarring between individuals, in part due to genetic factors. This study aimed to identify genetic variants associated with scar height and pliability after burn injury. An exome-wide array association study and gene pathway analysis were performed on a prospective cohort of 665 patients treated for burn injury. Outcomes were scar height (SH) and scar pliability (SP) sub-scores of the modified Vancouver Scar Scale (mVSS). DNA was genotyped using the Infinium^® HumanCoreExome-24 BeadChip. Associations between genetic variants (single nucleotide polymorphisms) and SH and SP were estimated using an additive genetic model adjusting for age, sex, number of surgical procedures and % total body surface area of burn in subjects of European ancestry. No individual genetic variants achieved the cut-off threshold of significance. Gene regions were analysed for spatially correlated single nucleotide polymorphisms and significant regions identified using comb-p software. This gene list was subject to gene pathway analysis to find which biological process terms were over-represented. Using this approach biological processes related to the nervous system and cell adhesion were the predominant gene pathways associated with both SH and SP. This study suggests genes associated with innervation may be important in scar fibrosis. Further studies using similar and larger datasets will be essential to validate these findings.

Pharmacogenomics study of thiazide diuretics and QT interval in multi-ethnic populations: the cohorts for heart and aging research in genomic epidemiology

Thiazide diuretics, commonly used antihypertensives, may cause QT interval (QT) prolongation, a risk factor for highly fatal and difficult to predict ventricular arrhythmias. We examined whether common single-nucleotide polymorphisms (SNPs) modified the association between thiazide use and QT or its component parts (QRS interval, JT interval) by performing ancestry-specific, trans-ethnic and cross-phenotype genome-wide analyses of European (66%), African American (15%) and Hispanic (19%) populations (N=78 199), leveraging longitudinal data, incorporating corrected standard errors to account for underestimation of interaction estimate variances and evaluating evidence for pathway enrichment. Although no loci achieved genome-wide significance (P<5 × 10-8), we found suggestive evidence (P<5 × 10-6) for SNPs modifying the thiazide-QT association at 22 loci, including ion transport loci (for example, NELL1, KCNQ3). The biologic plausibility of our suggestive results and simulations demonstrating modest power to detect interaction effects at genome-wide significant levels indicate that larger studies and innovative statistical methods are warranted in future efforts evaluating thiazide-SNP interactions.

A Review on Methods for Detecting SNP Interactions in High-Dimensional Genomic Data

In this era of genome-wide association studies (GWAS), the quest for understanding the genetic architecture of complex diseases is rapidly increasing more than ever before. The development of high throughput genotyping and next generation sequencing technologies enables genetic epidemiological analysis of large scale data. These advances have led to the identification of a number of single nucleotide polymorphisms (SNPs) responsible for disease susceptibility. The interactions between SNPs associated with complex diseases are increasingly being explored in the current literature. These interaction studies are mathematically challenging and computationally complex. These challenges have been addressed by a number of data mining and machine learning approaches. This paper reviews the current methods and the related software packages to detect the SNP interactions that contribute to diseases. The issues that need to be considered when developing these models are addressed in this review. The paper also reviews the achievements in data simulation to evaluate the performance of these models. Further, it discusses the future of SNP interaction analysis.

A Review of Pathway-Based Analysis Tools That Visualize Genetic Variants

Pathway analysis is a powerful method for data analysis in genomics, most often applied to gene expression analysis. It is also promising for single-nucleotide polymorphism (SNP) data analysis, such as genome-wide association study data, because it allows the interpretation of variants with respect to the biological processes in which the affected genes and proteins are involved. Such analyses support an interactive evaluation of the possible effects of variations on function, regulation or interaction of gene products. Current pathway analysis software often does not support data visualization of variants in pathways as an alternate method to interpret genetic association results, and specific statistical methods for pathway analysis of SNP data are not combined with these visualization features. In this review, we first describe the visualization options of the tools that were identified by a literature review, in order to provide insight for improvements in this developing field. Tool evaluation was performed using a computational epistatic dataset of gene–gene interactions for obesity risk. Next, we report the necessity to include in these tools statistical methods designed for the pathway-based analysis with SNP data, expressly aiming to define features for more comprehensive pathway-based analysis tools. We conclude by recognizing that pathway analysis of genetic variations data requires a sophisticated combination of the most useful and informative visual aspects of the various tools evaluated.

Genetic Analysis of Mitochondrial Ribosomal Proteins and Cognitive Aging in Postmenopausal Women

Genes encoding mitochondrial ribosomal proteins (MRPs) have been linked to aging and longevity in model organisms (i.e., mice, Caenorhabditis elegans). Here we evaluated if the MRPs have conserved effects on aging traits in humans. We utilized data from 4,504 participants of the Women's Health Initiative Memory Study (WHIMS) who had both longitudinal cognitive data and genetic data. Two aging phenotypes were considered: (1) gross lifespan (time to all-cause mortality), and (2) cognitive aging (longitudinal rate of change in modified mini-mental state scores). We tested genetic association with variants in 78 members of the MRP gene family. Genetic association tests were done at the single nucleotide polymorphism (SNP) level, and at gene-set level using two distinct procedures (GATES and MAGMA). We included SNPs in APOE and adjusted the tests for the APOE-ε4 allele, a known risk factor for dementia. The strongest association signal is for the known cognitive aging SNP, rs429358, in APOE (p-value = 5 × 10-28 for cognitive aging; p-value = 0.03 for survival). We found no significant association between the MRPs and survival time. For cognitive aging, we detected SNP level association for rs189661478 in MRPL23 (p-value < 9 × 10-6). Furthermore, the gene-set analysis showed modest but significant association between the MRP family and cognitive aging. In conclusion, our results indicate a potential pathway-level association between the MRPs and cognitive aging that is independent of the APOE locus. We however did not detect association between the MRPs and lifespan.

Sharing of Genes and Pathways Across Complex Phenotypes: A Multilevel Genome-Wide Analysis

Evidence from genome-wide association studies (GWAS) suggest that pleiotropic effects on human complex phenotypes are very common. Recently, an atlas of genetic correlations among complex phenotypes has broadened our understanding of human diseases and traits. Here, we examine genetic overlap, from a gene-centric perspective, among the same 24 phenotypes previously investigated for genetic correlations. After adopting the multilevel pipeline (freely available at http://grass.cgs.hku.hk/limx/kgg/), which includes intragenic single nucleotide polymorphisms (SNPs), genes, and gene-sets, to estimate genetic similarities across phenotypes, a large amount of sharing of several biologically related phenotypes was confirmed. In addition, significant genetic overlaps were also found among phenotype pairs that were previously unidentified by SNP-level approaches. All these pairs with new genetic links are supported by earlier epidemiological evidence, although only a few of them have pleiotropic genes in the GWAS Catalog. Hence, our gene and gene-set analyses are able to provide new insights into cross-phenotype connections. The investigation on genetic sharing at three different levels presents a complementary picture of how common DNA sequence variations contribute to disease comorbidities and trait manifestations.

Genome-wide pathway analysis identifies VEGF pathway association with oral ulceration in systemic lupus erythematosus

Background

Systemic lupus erythematosus (SLE) is a genetically complex rheumatic disease characterized by heterogeneous clinical manifestations of unknown etiology. Recent studies have suggested the existence of a genetic basis for SLE heterogeneity. The objective of the present study was to identify new genetic variation associated with the clinically relevant phenotypes in SLE.

Methods

A two-stage pathway-based approach was used to identify the genetic variation associated with the main clinical phenotypes in SLE. In the discovery stage, 482 SLE patients were genotyped using Illumina Human Quad610 microarrays. Association between 798 reference genetic pathways from the Molecular Signatures Database and 11 SLE phenotypes was tested using the set-based method implemented in PLINK software. Pathways significantly associated after multiple test correction were subsequently tested for replication in an independent cohort of 425 SLE patients. Using an in silico approach, we analyzed the functional effects of common SLE therapies on the replicated genetic pathways. The association of known SLE risk variants with the development of the clinical phenotypes was also analyzed.

Results

In the discovery stage, we found a significant association between the vascular endothelial growth factor (VEGF) pathway and oral ulceration (P value for false discovery rate (P_FDR) < 0.05), and between the negative regulation signaling pathway of retinoic acid inducible gene-I/melanoma differentiation associated gene 5 and the production of antinuclear antibodies (P_FDR < 0.05). In the replication stage, we validated the association between the VEGF pathway and oral ulceration. Therapies commonly used to treat mucocutaneous phenotypes in SLE were found to strongly influence VEGF pathway gene expression (P = 4.60e-4 to 5.38e-14). Analysis of known SLE risk loci identified a strong association between PTPN22 and the risk of hematologic disorder and with the development of antinuclear antibodies.

Conclusions

The present study has identified VEGF genetic pathway association with the risk of oral ulceration in SLE. New therapies targeting the VEGF pathway could be more effective in reducing the severity of this phenotype. These findings represent a first step towards the understanding of the genetic basis of phenotype heterogeneity in SLE.

Electronic supplementary material

The online version of this article (doi:10.1186/s13075-017-1345-6) contains supplementary material, which is available to authorized users.

Improving the detection of pathways in genome-wide association studies by combined effects of SNPs from Linkage Disequilibrium blocks

Genome-wide association studies (GWAS) have successfully identified single variants associated with diseases. To increase the power of GWAS, gene-based and pathway-based tests are commonly employed to detect more risk factors. However, the gene- and pathway-based association tests may be biased towards genes or pathways containing a large number of single-nucleotide polymorphisms (SNPs) with small P-values caused by high linkage disequilibrium (LD) correlations. To address such bias, numerous pathway-based methods have been developed. Here we propose a novel method, DGAT-path, to divide all SNPs assigned to genes in each pathway into LD blocks, and to sum the chi-square statistics of LD blocks for assessing the significance of the pathway by permutation tests. The method was proven robust with the type I error rate >1.6 times lower than other methods. Meanwhile, the method displays a higher power and is not biased by the pathway size. The applications to the GWAS summary statistics for schizophrenia and breast cancer indicate that the detected top pathways contain more genes close to associated SNPs than other methods. As a result, the method identified 17 and 12 significant pathways containing 20 and 21 novel associated genes, respectively for two diseases. The method is available online by http://sparks-lab.org/server/DGAT-path.

The statistical properties of gene-set analysis

The rapid increase in loci discovered in genome-wide association studies has created a need to understand the biological implications of these results. Gene-set analysis provides a means of gaining such understanding, but the statistical properties of gene-set analysis are not well understood, which compromises our ability to interpret its results. In this Analysis article, we provide an extensive statistical evaluation of the core structure that is inherent to all gene- set analyses and we examine current implementations in available tools. We show which factors affect valid and successful detection of gene sets and which provide a solid foundation for performing and interpreting gene-set analysis.

Discovery and fine-mapping of loci associated with MUFAs through trans-ethnic meta-analysis in Chinese and European populations

MUFAs are unsaturated FAs with one double bond and are derived from endogenous synthesis and dietary intake. Accumulating evidence has suggested that plasma and erythrocyte MUFA levels are associated with cardiometabolic disorders, including CVD, T2D, and metabolic syndrome (MS). Previous genome-wide association studies (GWASs) have identified seven loci for plasma and erythrocyte palmitoleic and oleic acid levels in populations of European origin. To identify additional MUFA-associated loci and the potential functional variant at each locus, we performed ethnic-specific GWAS meta-analyses and trans-ethnic meta-analyses in more than 15,000 participants of Chinese and European ancestry. We identified novel genome-wide significant associations for vaccenic acid at FADS1/2 and PKD2L1 [log10(Bayes factor) ≥ 8.07] and for gondoic acid at FADS1/2 and GCKR [log10(Bayes factor) ≥ 6.22], and also observed improved fine-mapping resolutions at FADS1/2 and GCKR loci. The greatest improvement was observed at GCKR, where the number of variants in the 99% credible set was reduced from 16 (covering 94.8 kb) to 5 (covering 19.6 kb, including a missense variant rs1260326) after trans-ethnic meta-analysis. We also confirmed the previously reported associations of PKD2L1, FADS1/2, GCKR, and HIF1AN with palmitoleic acid and of FADS1/2 and LPCAT3 with oleic acid in the Chinese-specific GWAS and the trans-ethnic meta-analyses. Pathway-based analyses suggested that the identified loci were in unsaturated FA metabolism and signaling pathways. Our findings provide novel insight into the genetic basis relevant to MUFA metabolism and biology.

Integrative Analysis of Genetic, Genomic, and Phenotypic Data for Ethanol Behaviors: A Network-Based Pipeline for Identifying Mechanisms and Potential Drug Targets

Complex behavioral traits, such as alcohol abuse, are caused by an interplay of genetic and environmental factors, producing deleterious functional adaptations in the central nervous system. The long-term behavioral consequences of such changes are of substantial cost to both the individual and society. Substantial progress has been made in the last two decades in understanding elements of brain mechanisms underlying responses to ethanol in animal models and risk factors for alcohol use disorder (AUD) in humans. However, treatments for AUD remain largely ineffective and few medications for this disease state have been licensed. Genome-wide genetic polymorphism analysis (GWAS) in humans, behavioral genetic studies in animal models and brain gene expression studies produced by microarrays or RNA-seq have the potential to produce nonbiased and novel insight into the underlying neurobiology of AUD. However, the complexity of such information, both statistical and informational, has slowed progress toward identifying new targets for intervention in AUD. This chapter describes one approach for integrating behavioral, genetic, and genomic information across animal model and human studies. The goal of this approach is to identify networks of genes functioning in the brain that are most relevant to the underlying mechanisms of a complex disease such as AUD. We illustrate an example of how genomic studies in animal models can be used to produce robust gene networks that have functional implications, and to integrate such animal model genomic data with human genetic studies such as GWAS for AUD. We describe several useful analysis tools for such studies: ComBAT, WGCNA, and EW_dmGWAS. The end result of this analysis is a ranking of gene networks and identification of their cognate hub genes, which might provide eventual targets for future therapeutic development. Furthermore, this combined approach may also improve our understanding of basic mechanisms underlying gene x environmental interactions affecting brain functioning in health and disease.

Gene- and pathway-based association tests for multiple traits with GWAS summary statistics

To identify novel genetic variants associated with complex traits and to shed new insights on underlying biology, in addition to the most popular single SNP-single trait association analysis, it would be useful to explore multiple correlated (intermediate) traits at the gene- or pathway-level by mining existing single GWAS or meta-analyzed GWAS data. For this purpose, we present an adaptive gene-based test and a pathway-based test for association analysis of multiple traits with GWAS summary statistics. The proposed tests are adaptive at both the SNP- and trait-levels; that is, they account for possibly varying association patterns (e.g. signal sparsity levels) across SNPs and traits, thus maintaining high power across a wide range of situations. Furthermore, the proposed methods are general: they can be applied to mixed types of traits, and to Z-statistics or P-values as summary statistics obtained from either a single GWAS or a meta-analysis of multiple GWAS. Our numerical studies with simulated and real data demonstrated the promising performance of the proposed methods.

AVAILABILITY AND IMPLEMENTATION

The methods are implemented in R package aSPU, freely and publicly available at: https://cran.r-project.org/web/packages/aSPU/ CONTACT: weip@biostat.umn.eduSupplementary information: Supplementary data are available at Bioinformatics online.

Pathway analysis in attention deficit hyperactivity disorder: An ensemble approach

Despite a wealth of evidence for the role of genetics in attention deficit hyperactivity disorder (ADHD), specific and definitive genetic mechanisms have not been identified. Pathway analyses, a subset of gene-set analyses, extend the knowledge gained from genome-wide association studies (GWAS) by providing functional context for genetic associations. However, there are numerous methods for association testing of gene sets and no real consensus regarding the best approach. The present study applied six pathway analysis methods to identify pathways associated with ADHD in two GWAS datasets from the Psychiatric Genomics Consortium. Methods that utilize genotypes to model pathway-level effects identified more replicable pathway associations than methods using summary statistics. In addition, pathways implicated by more than one method were significantly more likely to replicate. A number of brain-relevant pathways, such as RhoA signaling, glycosaminoglycan biosynthesis, fibroblast growth factor receptor activity, and pathways containing potassium channel genes, were nominally significant by multiple methods in both datasets. These results support previous hypotheses about the role of regulation of neurotransmitter release, neurite outgrowth and axon guidance in contributing to the ADHD phenotype and suggest the value of cross-method convergence in evaluating pathway analysis results. © 2016 Wiley Periodicals, Inc.

Adaptive gene- and pathway-trait association testing with GWAS summary statistics

BACKGROUND

Gene- and pathway-based analyses offer a useful alternative and complement to the usual single SNP-based analysis for GWAS. On the other hand, most existing gene- and pathway-based tests are not highly adaptive, and/or require the availability of individual-level genotype and phenotype data. It would be desirable to have highly adaptive tests applicable to summary statistics for single SNPs. This has become increasingly important given the popularity of large-scale meta-analyses of multiple GWASs and the practical availability of either single GWAS or meta-analyzed GWAS summary statistics for single SNPs.

RESULTS

We extend two adaptive tests for gene- and pathway-level association with a univariate trait to the case with GWAS summary statistics without individual-level genotype and phenotype data. We use the WTCCC GWAS data to evaluate and compare the proposed methods and several existing methods. We further illustrate their applications to a meta-analyzed dataset to identify genes and pathways associated with blood pressure, demonstrating the potential usefulness of the proposed methods. The methods are implemented in R package aSPU, freely and publicly available.

AVAILABILITY AND IMPLEMENTATION

https://cran.r-project.org/web/packages/aSPU/ CONTACT: weip@biostat.umn.edu

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Genome-Wide Pathway Analysis Identifies Genetic Pathways Associated with Psoriasis

Psoriasis is a chronic inflammatory disease with a complex genetic architecture. To date, the psoriasis heritability is only partially explained. However, there is increasing evidence that the missing heritability in psoriasis could be explained by multiple genetic variants of low effect size from common genetic pathways. The objective of this study was to identify new genetic variation associated with psoriasis risk at the pathway level. We genotyped 598,258 single nucleotide polymorphisms in a discovery cohort of 2,281 case-control individuals from Spain. We performed a genome-wide pathway analysis using 1,053 reference biological pathways. A total of 14 genetic pathways (PFDR ≤ 2.55 × 10(-2)) were found to be significantly associated with psoriasis risk. Using an independent validation cohort of 7,353 individuals from the UK, a total of 6 genetic pathways were significantly replicated (PFDR ≤ 3.46 × 10(-2)). We found genetic pathways that had not been previously associated with psoriasis risk such as retinol metabolism (Pcombined = 1.84 × 10(-4)), the transport of inorganic ions and amino acids (Pcombined = 1.57 × 10(-7)), and post-translational protein modification (Pcombined = 1.57 × 10(-7)). In the latter pathway, MGAT5 showed a strong network centrality, and its association with psoriasis risk was further validated in an additional case-control cohort of 3,429 individuals (P < 0.05). These findings provide insights into the biological mechanisms associated with psoriasis susceptibility.

Pathway-based variant enrichment analysis on the example of dilated cardiomyopathy

Genome-wide association (GWA) studies have significantly contributed to the understanding of human genetic variation and its impact on clinical traits. Frequently only a limited number of highly significant associations were considered as biologically relevant. Increasingly, network analysis of affected genes is used to explore the potential role of the genetic background on disease mechanisms. Instead of first determining affected genes or calculating scores for genes and performing pathway analysis on the gene level, we integrated both steps and directly calculated enrichment on the genetic variant level. The respective approach has been tested on dilated cardiomyopathy (DCM) GWA data as showcase. To compute significance values, 5000 permutation tests were carried out and p values were adjusted for multiple testing. For 282 KEGG pathways, we computed variant enrichment scores and significance values. Of these, 65 were significant. Surprisingly, we discovered the "nucleotide excision repair" and "tuberculosis" pathways to be most significantly associated with DCM (p = 10(-9)). The latter pathway is driven by genes of the HLA-D antigen group, a finding that closely resembles previous discoveries made by expression quantitative trait locus analysis in the context of DCM-GWA. Next, we implemented a sub-network-based analysis, which searches for affected parts of KEGG, however, independent on the pre-defined pathways. Here, proteins of the contractile apparatus of cardiac cells as well as the FAS sub-network were found to be affected by common polymorphisms in DCM. In this work, we performed enrichment analysis directly on variants, leveraging the potential to discover biological information in thousands of published GWA studies. The applied approach is cutoff free and considers a ranked list of genetic variants as input.

Gene set analysis: A step-by-step guide

To maximize the potential of genome-wide association studies, many researchers are performing secondary analyses to identify sets of genes jointly associated with the trait of interest. Although methods for gene-set analyses (GSA), also called pathway analyses, have been around for more than a decade, the field is still evolving. There are numerous algorithms available for testing the cumulative effect of multiple SNPs, yet no real consensus in the field about the best way to perform a GSA. This paper provides an overview of the factors that can affect the results of a GSA, the lessons learned from past studies, and suggestions for how to make analysis choices that are most appropriate for different types of data. © 2015 Wiley Periodicals, Inc.

A Powerful Pathway-Based Adaptive Test for Genetic Association with Common or Rare Variants

In spite of the success of genome-wide association studies (GWASs), only a small proportion of heritability for each complex trait has been explained by identified genetic variants, mainly SNPs. Likely reasons include genetic heterogeneity (i.e., multiple causal genetic variants) and small effect sizes of causal variants, for which pathway analysis has been proposed as a promising alternative to the standard single-SNP-based analysis. A pathway contains a set of functionally related genes, each of which includes multiple SNPs. Here we propose a pathway-based test that is adaptive at both the gene and SNP levels, thus maintaining high power across a wide range of situations with varying numbers of the genes and SNPs associated with a trait. The proposed method is applicable to both common variants and rare variants and can incorporate biological knowledge on SNPs and genes to boost statistical power. We use extensively simulated data and a WTCCC GWAS dataset to compare our proposal with several existing pathway-based and SNP-set-based tests, demonstrating its promising performance and its potential use in practice.

JAG: A Computational Tool to Evaluate the Role of Gene-Sets in Complex Traits

Gene-set analysis has been proposed as a powerful tool to deal with the highly polygenic architecture of complex traits, as well as with the small effect sizes typically found in GWAS studies for complex traits. We developed a tool, Joint Association of Genetic variants (JAG), which can be applied to Genome Wide Association (GWA) data and tests for the joint effect of all single nucleotide polymorphisms (SNPs) located in a user-specified set of genes or biological pathway. JAG assigns SNPs to genes and incorporates self-contained and/or competitive tests for gene-set analysis. JAG uses permutation to evaluate gene-set significance, which implicitly controls for linkage disequilibrium, sample size, gene size, the number of SNPs per gene and the number of genes in the gene-set. We conducted a power analysis using the Wellcome Trust Case Control Consortium (WTCCC) Crohn’s disease data set and show that JAG correctly identifies validated gene-sets for Crohn’s disease and has more power than currently available tools for gene-set analysis. JAG is a powerful, novel tool for gene-set analysis, and can be freely downloaded from the CTG Lab website.

Gene-based and pathway-based genome-wide association study of alcohol dependence

BACKGROUND

The organization of risk genes within signaling pathways may provide clues about the converging neurobiological effects of risk genes for alcohol dependence.

AIM

Identify risk genes and risk gene pathways for alcohol dependence.

METHODS

We conducted a pathway-based genome-wide association study (GWAS) of alcohol dependence using a gene-set-rich analytic approach. Approximately one million genetic markers were tested in the discovery sample which included 1409 European-American (EA) alcohol dependent individuals and 1518 EA healthy comparison subjects. An additional 681 African-American (AA) cases and 508 AA healthy subjects served as the replication sample.

RESULTS

We identified several genome-wide replicable risk genes and risk pathways that were significantly associated with alcohol dependence. After applying the Bonferroni correction for multiple testing, the 'cellextracellular matrix interactions' pathway (p<2.0E-4 in EAs) and the PXN gene (which encodes paxillin) (p=3.9E-7 in EAs) within this pathway were the most promising risk factors for alcohol dependence. There were also two nominally replicable pathways enriched in alcohol dependence-related genes in both EAs (0.015≤p≤0.035) and AAs (0.025≤p≤0.050): the 'Na+/Cl- dependent neurotransmitter transporters' pathway and the 'other glycan degradation' pathway.

CONCLUSION

These findings provide new evidence highlighting several genes and biological signaling processes that may be related to the risk for alcohol dependence.

Effect of occupational exposures on lung cancer susceptibility: a study of gene-environment interaction analysis

BACKGROUND

Occupational exposures are known risk factors for lung cancer. Role of genetically determined host factors in occupational exposure-related lung cancer is unclear.

METHODS

We used genome-wide association (GWA) data from a case-control study conducted in 6 European countries from 1998 to 2002 to identify gene-occupation interactions and related pathways for lung cancer risk. GWA analysis was performed for each exposure using logistic regression and interaction term for genotypes, and exposure was included in this model. Both SNP-based and gene-based interaction P values were calculated. Pathway analysis was performed using three complementary methods, and analyses were adjusted for multiple comparisons. We analyzed 312,605 SNPs and occupational exposure to 70 agents from 1,802 lung cancer cases and 1,725 cancer-free controls.

RESULTS

Mean age of study participants was 60.1 ± 9.1 years and 75% were male. Largest number of significant associations (P ≤ 1 × 10(-5)) at SNP level was demonstrated for nickel, brick dust, concrete dust, and cement dust, and for brick dust and cement dust at the gene-level (P ≤ 1 × 10(-4)). Approximately 14 occupational exposures showed significant gene-occupation interactions with pathways related to response to environmental information processing via signal transduction (P < 0.001 and FDR < 0.05). Other pathways that showed significant enrichment were related to immune processes and xenobiotic metabolism.

CONCLUSION

Our findings suggest that pathways related to signal transduction, immune process, and xenobiotic metabolism may be involved in occupational exposure-related lung carcinogenesis.

IMPACT

Our study exemplifies an integrative approach using pathway-based analysis to demonstrate the role of genetic variants in occupational exposure-related lung cancer susceptibility. Cancer Epidemiol Biomarkers Prev; 24(3); 570-9. ©2015 AACR.

Pathway analysis of genome-wide association study and transcriptome data highlights new biological pathways in colorectal cancer

Colorectal cancer (CRC) is a common malignancy that meets the definition of a complex disease. Genome-wide association study (GWAS) has identified several loci of weak predictive value in CRC, however, these do not fully explain the occurrence risk. Recently, gene set analysis has allowed enhanced interpretation of GWAS data in CRC, identifying a number of metabolic pathways as important for disease pathogenesis. Whether there are other important pathways involved in CRC, however, remains unclear. We present a systems analysis of KEGG pathways in CRC using (1) a human CRC GWAS dataset and (2) a human whole transcriptome CRC case-control expression dataset. Analysis of the GWAS dataset revealed significantly enriched KEGG pathways related to metabolism, immune system and diseases, cellular processes, environmental information processing, genetic information processing, and neurodegenerative diseases. Altered gene expression was confirmed in these pathways using the transcriptome dataset. Taken together, these findings not only confirm previous work in this area, but also highlight new biological pathways whose deregulation is critical for CRC. These results contribute to our understanding of disease-causing mechanisms and will prove useful for future genetic and functional studies in CRC.

Systematic permutation testing in GWAS pathway analyses: identification of genetic networks in dilated cardiomyopathy and ulcerative colitis

Background

Genome wide association studies (GWAS) are applied to identify genetic loci, which are associated with complex traits and human diseases. Analogous to the evolution of gene expression analyses, pathway analyses have emerged as important tools to uncover functional networks of genome-wide association data. Usually, pathway analyses combine statistical methods with a priori available biological knowledge. To determine significance thresholds for associated pathways, correction for multiple testing and over-representation permutation testing is applied.

Results

We systematically investigated the impact of three different permutation test approaches for over-representation analysis to detect false positive pathway candidates and evaluate them on genome-wide association data of Dilated Cardiomyopathy (DCM) and Ulcerative Colitis (UC). Our results provide evidence that the gold standard - permuting the case–control status – effectively improves specificity of GWAS pathway analysis. Although permutation of SNPs does not maintain linkage disequilibrium (LD), these permutations represent an alternative for GWAS data when case–control permutations are not possible. Gene permutations, however, did not add significantly to the specificity. Finally, we provide estimates on the required number of permutations for the investigated approaches.

Conclusions

To discover potential false positive functional pathway candidates and to support the results from standard statistical tests such as the Hypergeometric test, permutation tests of case control data should be carried out. The most reasonable alternative was case–control permutation, if this is not possible, SNP permutations may be carried out. Our study also demonstrates that significance values converge rapidly with an increasing number of permutations. By applying the described statistical framework we were able to discover axon guidance, focal adhesion and calcium signaling as important DCM-related pathways and Intestinal immune network for IgA production as most significant UC pathway.

Set-based joint test of interaction between SNPs in the VEGF pathway and exogenous estrogen finds association with age-related macular degeneration

Purpose:Age-Related Macular Degeneration (AMD) is the leading cause of irreversible visual loss in developed countries. Its etiology includes genetic and environmental factors. Although VEGFA variants are associated with AMD, the joint action of variants within the VEGF pathway and their interaction with non-genetic factors has not been investigated. Methods:Affymetrix 6.0 chipsets were used to genotype 668,238 SNPs in 1,207 AMD cases and 686 controls. Environmental exposures were collected by questionnaire. A set-based test was conducted using the chi-square statistic at each SNP derived from Kraft's 2df joint test. Pathway and gene-based test statistics were calculated as the mean of all independent SNP statistics. Phenotype labels were permuted 10,000 times to generate an empirical p-value. Results: While a main effect of the VEGF pathway was not identified, the pathway was associated with neovascular AMD in women when accounting for birth control pill (BCP) use (P= 0.017). Analysis of VEGF's subpathways found that SNPs in the Proliferation subpathway were associated with neovascular AMD (P=0.029) when accounting for BCP use. Nominally significant genes within this subpathway were also observed. Stratification by BCP use revealed novel significant genetic effects in women who had taken BCPs. Conclusions: These results illustrate that some AMD genetic risk factors may only be revealed when considering complex relationships among risk factors. This shows the utility of exploring pathways of previously associated genes to find novel effects. It also demonstrates the importance of incorporating environmental exposures in tests of genetic association at the SNP, gene, or pathway level.

Integrative pathway-based approach for genome-wide association studies: identification of new pathways for rheumatoid arthritis and type 1 diabetes

Genome-wide association studies (GWAS) led to the identification of numerous novel loci for a number of complex diseases. Pathway-based approaches using genotypic data provide tangible leads which cannot be identified by single marker approaches as implemented in GWAS. The available pathway analysis approaches mainly differ in the employed databases and in the applied statistics for determining the significance of the associated disease markers.

So far, pathway-based approaches using GWAS data failed to consider the overlapping of genes among different pathways or the influence of protein–interactions. We performed a multistage integrative pathway (MIP) analysis on three common diseases - Crohn's disease (CD), rheumatoid arthritis (RA) and type 1 diabetes (T1D) - incorporating genotypic, pathway, protein- and domain-interaction data to identify novel associations between these diseases and pathways. Additionally, we assessed the sensitivity of our method by studying the influence of the most significant SNPs on the pathway analysis by removing those and comparing the corresponding pathway analysis results. Apart from confirming many previously published associations between pathways and RA, CD and T1D, our MIP approach was able to identify three new associations between disease phenotypes and pathways. This includes a relation between the influenza-A pathway and RA, as well as a relation between T1D and the phagosome and toxoplasmosis pathways. These results provide new leads to understand the molecular underpinnings of these diseases.

The developed software herein used is available at http://www.cogsys.cs.uni-tuebingen.de/software/GWASPathwayIdentifier/index.htm.

Measles contributes to rheumatoid arthritis: evidence from pathway and network analyses of genome-wide association studies

Growing evidence from epidemiological studies indicates the association between rheumatoid arthritis (RA) and measles. However, the exact mechanism for this association is still unclear now. We consider that the strong association between both diseases may be caused by shared genetic pathways. We performed a pathway analysis of large-scale RA genome-wide association studies (GWAS) dataset with 5,539 cases and 20,169 controls of European descent. Meanwhile, we evaluated our findings using previously identified RA loci, protein-protein interaction network and previous results from pathway analysis of RA and other autoimmune diseases GWAS. We confirmed four pathways including Cytokine-cytokine receptor interaction, Jak-STAT signaling, T cell receptor signaling and Cell adhesion molecules. Meanwhile, we highlighted for the first time the involvement of Measles and Intestinal immune network for IgA production pathways in RA. Our results may explain the strong association between RA and measles, which may be caused by the shared genetic pathway. We believe that our results will be helpful for future genetic studies in RA pathogenesis and may significantly assist in the development of therapeutic strategies.

Pathway-PDT: a flexible pathway analysis tool for nuclear families

Background

Pathway analysis based on Genome-Wide Association Studies (GWAS) data has become popular as a secondary analysis strategy. Although many pathway analysis tools have been developed for case–control studies, there is no tool that can use all information from raw genotypes in general nuclear families. We developed Pathway-PDT, which uses the framework of Pedigree Disequilibrium Test (PDT) for general family data, to perform pathway analysis based on raw genotypes in family-based GWAS.

Results

Simulation results showed that Pathway-PDT is more powerful than the p-value based method, ALIGATOR. Pathway-PDT also can be more powerful than the PLINK set-based test when analyzing general nuclear families with multiple siblings or missing parents. Additionally, Pathway-PDT has a flexible and convenient user interface, which allows users to modify their analysis parameters as well as to apply various types of gene and pathway definitions.

Conclusions

The Pathway-PDT method is implemented in C++ with POSIX threads and is computationally feasible for pathway analysis with large scale family GWAS datasets. The Windows binary along with Makefile and source codes for the Linux are available at https://sourceforge.net/projects/pathway-pdt/.

Pathway analysis using genome-wide association study data for coronary restenosis--a potential role for the PARVB gene

Background

Coronary restenosis after percutaneous coronary intervention (PCI) still remains a significant limitation of the procedure. The causative mechanisms of restenosis have not yet been fully identified. The goal of the current study was to perform gene-set analysis of biological pathways related to inflammation, proliferation, vascular function and transcriptional regulation on coronary restenosis to identify novel genes and pathways related to this condition.

Methods

The GENetic DEterminants of Restenosis (GENDER) databank contains genotypic data of 556,099SNPs of 295 cases with restenosis and 571 matched controls. Fifty-four pathways, related to known restenosis-related processes, were selected. Gene-set analysis was performed using PLINK, GRASS and ALIGATOR software. Pathways with a p<0.01 were fine-mapped and significantly associated SNPs were analyzed in an independent replication cohort.

Results

Six pathways (cell-extracellular matrix (ECM) interactions pathway, IL2 signaling pathway, IL6 signaling pathway, platelet derived growth factor pathway, vitamin D receptor pathway and the mitochondria pathway) were significantly associated in one or two of the software packages. Two SNPs in the cell-ECM interactions pathway were replicated in an independent restenosis cohort. No replication was obtained for the other pathways.

Conclusion

With these results we demonstrate a potential role of the cell-ECM interactions pathway in the development of coronary restenosis. These findings contribute to the increasing knowledge of the genetic etiology of restenosis formation and could serve as a hypothesis-generating effort for further functional studies.

Genome-wide pathway analysis of memory impairment in the Alzheimer's Disease Neuroimaging Initiative (ADNI) cohort implicates gene candidates, canonical pathways, and networks

Memory deficits are prominent features of mild cognitive impairment (MCI) and Alzheimer's disease (AD). The genetic architecture underlying these memory deficits likely involves the combined effects of multiple genetic variants operative within numerous biological pathways. In order to identify functional pathways associated with memory impairment, we performed a pathway enrichment analysis on genome-wide association data from 742 Alzheimer's Disease Neuroimaging Initiative (ADNI) participants. A composite measure of memory was generated as the phenotype for this analysis by applying modern psychometric theory to item-level data from the ADNI neuropsychological test battery. Using the GSA-SNP software tool, we identified 27 canonical, expertly-curated pathways with enrichment (FDR-corrected p-value < 0.05) against this composite memory score. Processes classically understood to be involved in memory consolidation, such as neurotransmitter receptor-mediated calcium signaling and long-term potentiation, were highly represented among the enriched pathways. In addition, pathways related to cell adhesion, neuronal differentiation and guided outgrowth, and glucose- and inflammation-related signaling were also enriched. Among genes that were highly-represented in these enriched pathways, we found indications of coordinated relationships, including one large gene set that is subject to regulation by the SP1 transcription factor, and another set that displays co-localized expression in normal brain tissue along with known AD risk genes. These results 1) demonstrate that psychometrically-derived composite memory scores are an effective phenotype for genetic investigations of memory impairment and 2) highlight the promise of pathway analysis in elucidating key mechanistic targets for future studies and for therapeutic interventions.

Non-homologous end-joining pathway associated with occurrence of myocardial infarction: gene set analysis of genome-wide association study data

PURPOSE

DNA repair deficiencies have been postulated to play a role in the development and progression of cardiovascular disease (CVD). The hypothesis is that DNA damage accumulating with age may induce cell death, which promotes formation of unstable plaques. Defects in DNA repair mechanisms may therefore increase the risk of CVD events. We examined whether the joints effect of common genetic variants in 5 DNA repair pathways may influence the risk of CVD events.

METHODS

The PLINK set-based test was used to examine the association to myocardial infarction (MI) of the DNA repair pathway in GWAS data of 866 subjects of the GENetic DEterminants of Restenosis (GENDER) study and 5,244 subjects of the PROspective Study of Pravastatin in the Elderly at Risk (PROSPER) study. We included the main DNA repair pathways (base excision repair, nucleotide excision repair, mismatch repair, homologous recombination and non-homologous end-joining (NHEJ)) in the analysis.

RESULTS

The NHEJ pathway was associated with the occurrence of MI in both GENDER (P = 0.0083) and PROSPER (P = 0.014). This association was mainly driven by genetic variation in the MRE11A gene (PGENDER = 0.0001 and PPROSPER = 0.002). The homologous recombination pathway was associated with MI in GENDER only (P = 0.011), for the other pathways no associations were observed.

CONCLUSION

This is the first study analyzing the joint effect of common genetic variation in DNA repair pathways and the risk of CVD events, demonstrating an association between the NHEJ pathway and MI in 2 different cohorts.

Pathway-based analysis of primary biliary cirrhosis genome-wide association studies

Genome-wide association studies (GWAS) have successfully identified several loci associated with primary biliary cirrhosis (PBC) risk. Pathway analysis complements conventional GWAS analysis. We applied the recently developed linear combination test for pathways to datasets drawn from independent PBC GWAS in Italian and Canadian subjects. Of the Kyoto Encyclopedia of Genes and Genomes and BioCarta pathways tested, 25 pathways in the Italian dataset (449 cases, 940 controls) and 26 pathways in the Canadian dataset (530 cases, 398 controls) were associated with PBC susceptibility (P<0.05). After correcting for multiple comparisons, only the eight most significant pathways in the Italian dataset had FDR <0.25 with tumor necrosis factor/stress-related signaling emerging as the top pathway (P=7.38 × 10⁻⁴, FDR=0.18). Two pathways, phosphatidylinositol signaling and hedgehog signaling, were replicated in both datasets (P<0.05), and subjected to two additional complementary pathway tests. Both pathway signals remained significant in the Italian dataset on modified gene set enrichment analysis (P<0.05). In both GWAS, variants nominally associated with PBC were significantly overrepresented in the phosphatidylinositol pathway (Fisher exact P<0.05). These results point to established and novel pathway-level associations with inherited predisposition to PBC that, on further independent replication and functional validation, may provide fresh insights into PBC etiology.

Genetic association test for multiple traits at gene level

Genome-wide association studies (GWASs) at the gene level are commonly used to understand biological mechanisms underlying complex diseases. In general, one response or outcome is used to present a disease of interest in such studies. In this study, we consider a multiple traits association test from the gene level. We propose and examine a class of test statistics that summarizes the association information between single nucleotide polymorphisms (SNPs) and each of the traits. Our simulation studies demonstrate the advantage of gene-based multiple traits association tests when multiple traits share common genes. Using our proposed tests, we reanalyze the dataset from the Study of Addiction: Genetics and Environment (SAGE). Our result validates previous findings while presenting stronger evidence for consideration of multiple traits.

Integrative pathway analysis of genome-wide association studies and gene expression data in prostate cancer

BACKGROUND

Pathway analysis of large-scale omics data assists us with the examination of the cumulative effects of multiple functionally related genes, which are difficult to detect using the traditional single gene/marker analysis. So far, most of the genomic studies have been conducted in a single domain, e.g., by genome-wide association studies (GWAS) or microarray gene expression investigation. A combined analysis of disease susceptibility genes across multiple platforms at the pathway level is an urgent need because it can reveal more reliable and more biologically important information.

RESULTS

We performed an integrative pathway analysis of a GWAS dataset and a microarray gene expression dataset in prostate cancer. We obtained a comprehensive pathway annotation set from knowledge-based public resources, including KEGG pathways and the prostate cancer candidate gene set, and gene sets specifically defined based on cross-platform information. By leveraging on this pathway collection, we first searched for significant pathways in the GWAS dataset using four methods, which represent two broad groups of pathway analysis approaches. The significant pathways identified by each method varied greatly, but the results were more consistent within each method group than between groups. Next, we conducted a gene set enrichment analysis of the microarray gene expression data and found 13 pathways with cross-platform evidence, including "Fc gamma R-mediated phagocytosis" (P GWAS = 0.003, P expr < 0.001, and P combined = 6.18 × 10(-8)), "regulation of actin cytoskeleton" (P GWAS = 0.003, P expr = 0.009, and P combined = 3.34 × 10(-4)), and "Jak-STAT signaling pathway" (P GWAS = 0.001, P expr = 0.084, and P combined = 8.79 × 10(-4)).

CONCLUSIONS

Our results provide evidence at both the genetic variation and expression levels that several key pathways might have been involved in the pathological development of prostate cancer. Our framework that employs gene expression data to facilitate pathway analysis of GWAS data is not only feasible but also much needed in studying complex disease.

Pathway analysis of smoking quantity in multiple GWAS identifies cholinergic and sensory pathways

Cigarette smoking is a common addiction that increases the risk for many diseases, including lung cancer and chronic obstructive pulmonary disease. Genome-wide association studies (GWAS) have successfully identified and validated several susceptibility loci for nicotine consumption and dependence. However, the trait variance explained by these genes is only a small fraction of the estimated genetic risk. Pathway analysis complements single marker methods by including biological knowledge into the evaluation of GWAS, under the assumption that causal variants lie in functionally related genes, enabling the evaluation of a broad range of signals. Our approach to the identification of pathways enriched for multiple genes associated with smoking quantity includes the analysis of two studies and the replication of common findings in a third dataset. This study identified pathways for the cholinergic receptors, which included SNPs known to be genome-wide significant; as well as novel pathways, such as genes involved in the sensory perception of smell, that do not contain any single SNP that achieves that stringent threshold.

Systematic testing of literature reported genetic variation associated with coronary restenosis: results of the GENDER Study

Background

Coronary restenosis after percutaneous coronary intervention still remains a significant problem, despite all medical advances. Unraveling the mechanisms leading to restenosis development remains challenging. Many studies have identified genetic markers associated with restenosis, but consistent replication of the reported markers is scarce. The aim of the current study was to analyze the joined effect of previously in literature reported candidate genes for restenosis in the GENetic DEterminants of Restenosis (GENDER) databank.

Methodology/Principal Findings

Candidate genes were selected using a MEDLINE search including the terms ‘genetic polymorphism’ and ‘coronary restenosis’. The final set included 36 genes. Subsequently, all single nucleotide polymorphisms (SNPs) in the genomic region of these genes were analyzed in GENDER using set-based analysis in PLINK. The GENDER databank contains genotypic data of 2,571,586 SNPs of 295 cases with restenosis and 571 matched controls. The set, including all 36 literature reported genes, was, indeed, significantly associated with restenosis, p = 0.024 in the GENDER study. Subsequent analyses of the individual genes demonstrated that the observed association of the complete set was determined by 6 of the 36 genes.

Conclusion

Despite overt inconsistencies in literature, with regard to individual candidate gene studies, this is the first study demonstrating that the joint effect of all these genes together, indeed, is associated with restenosis.

A novel differential diagnostic model based on multiple biological parameters for immunoglobulin A nephropathy

BACKGROUND

Immunoglobulin A nephropathy (IgAN) is the most common form of glomerulonephritis in China. An accurate diagnosis of IgAN is dependent on renal biopsies, and there is lack of non-invasive and practical classification methods for discriminating IgAN from other primary kidney diseases. The objective of this study was to develop a classification model for the auxiliary diagnosis of IgAN using multiparameter analysis with various biological parameters.

METHODS

To establish an optimal classification model, 121 cases (58 IgAN vs. 63 non-IgAN) were recruited and statistically analyzed. The model was then validated in another 180 cases.

RESULTS

Of the 57 biological parameters, there were 16 parameters that were significantly different (P < 0.05) between IgAN and non-IgAN. The combination of fibrinogen, serum immunoglobulin A level, and manifestation was found to be significant in predicting IgAN. The validation accuracies of the logistic regression and discriminant analysis models were 77.5 and 77.0%, respectively at a predictive probability cut-off of 0.5, and 81.1 and 79.9%, respectively, at a predictive probability cut-off of 0.40. When the predicted probability of the equation containing the combination of fibrinogen, serum IgA level, and manifestation was more than 0.59, a patient had at least an 85.0% probability of having IgAN. When the predicted probability was lower than 0.26, a patient had at least an 88.5% probability of having non-IgAN. The results of the net reclassification improvement certificated serum Immunoglobulin A and fibrinogen had classification power for discriminating IgAN from non-IgAN.

CONCLUSIONS

These models possess potential clinical applications in distinguishing IgAN from other primary kidney diseases.

Pathway analysis of genomic data: concepts, methods, and prospects for future development

Genome-wide data sets are increasingly being used to identify biological pathways and networks underlying complex diseases. In particular, analyzing genomic data through sets defined by functional pathways offers the potential of greater power for discovery and natural connections to biological mechanisms. With the burgeoning availability of next-generation sequencing, this is an opportune moment to revisit strategies for pathway-based analysis of genomic data. Here, we synthesize relevant concepts and extant methodologies to guide investigators in study design and execution. We also highlight ongoing challenges and proposed solutions. As relevant analytical strategies mature, pathways and networks will be ideally placed to integrate data from diverse -omics sources to harness the extensive, rich information related to disease and treatment mechanisms.