Genome-wide gene and pathway analysis

Identification of a Predictive Model for Schizophrenia Based on SNPs in a Chinese Population

Background

Schizophrenia is a devastating mental disease with high heritability. A growing number of susceptibility genes associated with schizophrenia, as well as their corresponding SNPs loci, have been revealed by genome-wide association studies. However, using SNPs as predictors of disease and diagnosis remains difficult. Here, we aimed to uncover susceptibility SNPs in a Chinese population and to construct a prediction model for schizophrenia.

Methods

A total of 210 participants, including 70 patients with schizophrenia, 70 patients with bipolar disorder, and 70 healthy controls, were enrolled in this study. We estimated 14 SNPs using published risk loci of schizophrenia, and used these SNPs to build a model for predicting schizophrenia via comparison of genotype frequencies and regression. We evaluated the efficacy of the diagnostic model in schizophrenia and control patients using ROC curves and then used the 70 patients with bipolar disorder to evaluate the model's differential diagnostic efficacy.

Results

5 SNPs were selected to construct the model: rs148415900, rs71428218, rs4666990, rs112222723 and rs1716180. Correlation analysis results suggested that, compared with the risk SNP of 0, the risk SNP of 3 was associated with an increased risk of schizophrenia (OR = 13.00, 95% CI: 2.35-71.84, p = 0.003). The ROC-AUC of this prediction model for schizophrenia was 0.719 (95% CI: 0.634-0.804), with the greatest sensitivity and specificity being 60% and 80%, respectively. The ROC-AUC of the model in distinguishing between schizophrenia and bipolar disorder was 0.591 (95% CI: 0.497-0.686), with the greatest sensitivity and specificity being 60% and 55.7%, respectively.

Conclusion

The SNP risk score prediction model had good performance in predicting schizophrenia. To the best of our knowledge, previous studies have not applied SNP-based models to differentiate between cases of schizophrenia and other mental illnesses. It could have several potential clinical applications, including shaping disease diagnosis, treatment, and outcomes.

Systematic review of gastric cancer-associated genetic variants, gene-based meta-analysis, and gene-level functional analysis to identify candidate genes for drug development

Objective: Despite being a powerful tool to identify novel variants, genome-wide association studies (GWAS) are not sufficient to explain the biological function of variants. In this study, we aimed to elucidate at the gene level the biological mechanisms involved in gastric cancer (GC) development and to identify candidate drug target genes. Materials and methods: We conducted a systematic review for GWAS on GC following the PRISMA guidelines. Single nucleotide polymorphism (SNP)-level meta-analysis and gene-based analysis (GBA) were performed to identify SNPs and genes significantly associated with GC. Expression quantitative trait loci (eQTL), disease network, pathway enrichment, gene ontology, gene-drug, and chemical interaction analyses were conducted to elucidate the function of the genes identified by GBA. Results: A review of GWAS on GC identified 226 SNPs located in 91 genes. In the comprehensive GBA, 44 genes associated with GC were identified, among which 12 genes (THBS3, GBAP1, KRTCAP2, TRIM46, HCN3, MUC1, DAP3, EFNA1, MTX1, PRKAA1, PSCA, and ABO) were eQTL. Using disease network and pathway analyses, we identified that PRKAA, THBS3, and EFNA1 were significantly associated with the PI3K-Alt-mTOR-signaling pathway, which is involved in various oncogenic processes, and that MUC1 acts as a regulator in both the PI3K-Alt-mTOR and P53 signaling pathways. Furthermore, RPKAA1 had the highest number of interactions with drugs and chemicals. Conclusion: Our study suggests that PRKAA1, a gene in the PI3K-Alt-mTOR-signaling pathway, could be a potential target gene for drug development associated with GC in the future. Systematic Review Registration: website, identifier registration number.

Omnibus testing approach for gene-based gene-gene interaction

Genetic interaction is considered as one of the main heritable component of complex traits. With the emergence of genome-wide association studies (GWAS), a collection of statistical methods dedicated to the identification of interaction at the SNP level have been proposed. More recently, gene-based gene-gene interaction testing has emerged as an attractive alternative as they confer advantage in both statistical power and biological interpretation. Most of the gene-based interaction methods rely on a multidimensional modeling of the interaction, thus facing a lack of robustness against the huge space of interaction patterns. In this paper, we study a global testing approaches to address the issue of gene-based gene-gene interaction. Based on a logistic regression modeling framework, all SNP-SNP interaction tests are combined to produce a gene-level test for interaction. We propose an omnibus test that takes advantage of (1) the heterogeneity between existing global tests and (2) the complementarity between allele-based and genotype-based coding of SNPs. Through an extensive simulation study, it is demonstrated that the proposed omnibus test has the ability to detect with high power the most common interaction genetic models with one causal pair as well as more complex genetic models where more than one causal pair is involved. On the other hand, the flexibility of the proposed approach is shown to be robust and improves power compared to single global tests in replication studies. Furthermore, the application of our procedure to real datasets confirms the adaptability of our approach to replicate various gene-gene interactions.

Metagenomics and chemotherapy-induced nausea: A roadmap for future research

Uncontrolled chemotherapy-induced nausea and vomiting can reduce patients' quality of life and may result in premature discontinuation of chemotherapy. Although nausea and vomiting are commonly grouped together, research has shown that antiemetics are clinically effective against chemotherapy-induced vomiting (CIV) but less so against chemotherapy-induced nausea (CIN). Nausea remains a problem for up to 68% of patients who are prescribed guideline-consistent antiemetics. Despite the high prevalence of CIN, relatively little is known regarding its etiology independent of CIV. This review summarizes a metagenomics approach to the study and treatment of CIN with the goal of encouraging future research. Metagenomics focuses on genetic risk factors and encompasses both human (ie, host) and gut microbial genetic variation. Little work to date has focused on metagenomics as a putative biological mechanism of CIN. Metagenomics has the potential to be a powerful tool in advancing scientific understanding of CIN by identifying new biological pathways and intervention targets. The investigation of metagenomics in the context of well-established demographic, clinical, and patient-reported risk factors may help to identify patients at risk and facilitate the prevention and management of CIN.

Using common genetic variants to find drugs for common epilepsies

Better drugs are needed for common epilepsies. Drug repurposing offers the potential of significant savings in the time and cost of developing new treatments. In order to select the best candidate drug(s) to repurpose for a disease, it is desirable to predict the relative clinical efficacy that drugs will have against the disease. Common epilepsy can be divided into different types and syndromes. Different antiseizure medications are most effective for different types and syndromes of common epilepsy. For predictions of antiepileptic efficacy to be clinically translatable, it is essential that the predictions are specific to each form of common epilepsy, and reflect the patterns of drug efficacy observed in clinical studies and practice. These requirements are not fulfilled by previously published drug predictions for epilepsy. We developed a novel method for predicting the relative efficacy of drugs against any common epilepsy, by using its Genome-Wide Association Study summary statistics and drugs' activity data. The methodological advancement in our technique is that the drug predictions for a disease are based upon drugs' effects on the function and abundance of proteins, and the magnitude and direction of those effects, relative to the importance, degree and direction of the proteins' dysregulation in the disease. We used this method to predict the relative efficacy of all drugs, licensed for any condition, against each of the major types and syndromes of common epilepsy. Our predictions are concordant with findings from real-world experience and randomized clinical trials. Our method predicts the efficacy of existing antiseizure medications against common epilepsies; in this prediction, our method outperforms the best alternative existing method: area under receiver operating characteristic curve (mean ± standard deviation) 0.83 ± 0.03 and 0.63 ± 0.04, respectively. Importantly, our method predicts which antiseizure medications are amongst the more efficacious in clinical practice, and which antiseizure medications are amongst the less efficacious in clinical practice, for each of the main syndromes of common epilepsy, and it predicts the distinct order of efficacy of individual antiseizure medications in clinical trials of different common epilepsies. We identify promising candidate drugs for each of the major syndromes of common epilepsy. We screen five promising predicted drugs in an animal model: each exerts a significant dose-dependent effect upon seizures. Our predictions are a novel resource for selecting suitable candidate drugs that could potentially be repurposed for each of the major syndromes of common epilepsy. Our method is potentially generalizable to other complex diseases.

DNA polymorphisms of FGFBP1, leptin, κ-casein, and αs1-casein genes and their association with reproductive performance in dromedary she-camels

In livestock, selection programs on reproductive traits are time consuming due to long generation interval and low heritability; however understanding the genetic factors that influence reproductive performance facilitates an accurate selection and genetic improvement of these quantitative traits. Therefore, the aim of this study was to identify the effect of lactation order and polymorphisms of FGFBP1, leptin, κ-casein, and αs1-casein genes on reproductive performance in dromedary she-camels. Blood samples were collected from one hundred she-camels using vacutainer tube containing EDTA as an anticoagulant for DNA extraction. PCR-DNA sequencing approach was performed to investigate the polymorphisms in fragments for coding sequence (CDS) parts of FGFBP1 (539-bp) and leptin (359-bp), 5'-flanking region of κ-casein (488-bp) and intervening non-coding fragment spanning from exon 4 to exon 6 (924-bp) of αs1-casein genes. Results revealed that there were single nucleotide polymorphisms (SNPs) in the investigated genes between enrolled she-camels except for leptin gene that elicited a monomorphic pattern. The identified SNPs allocated camels into different groups according to which animal harbored the explored SNPs. Multiple linear regression analysis (MLR) reveled that identified SNPs and parity order had significant effects on age at first calving (AFC), days open (DO), calving interval (CI), number of services/conception (NS/C) and milk yield traits. Determination coefficients (R²) value indicated also that variability in the phenotypic measurements of studied traits might be correlated with the identified SNPs in the reproduction related genes. This study highlights that parity order has a remarkable effect on reproductive performance. Single nucleotide polymorphisms in the investigated genes could be also used as proxy predictors and genetic markers that permit a marker assisted selection for reproductive performance traits in dromedary she-camels.

Targeting Common Signaling Pathways for the Treatment of Stroke and Alzheimer's: a Comprehensive Review

Neurodegenerative diseases such as stroke and Alzheimer's disease (AD) are two inter-related disorders that affect the neurons in the brain and central nervous system. Alzheimer's is a disease by undefined origin and causes. Stroke and its most common type, ischemic stroke (IS), occurs due to the blockade of cerebral blood vessels. As an important feature, both of disorders are associated with irreversible damages to the brain and nervous system. In this regard, finding common signaling pathways and the same molecular origin between these two diseases may be a promising way for their solution. On the basis of literature appraisal, the most common signaling cascades implicated in the pathogenesis of AD and stroke including notch, autophagy, inflammatory, and insulin signaling pathways were reviewed. Furthermore, current therapeutic strategies including natural and synthetic pharmaceuticals aiming modulation of respective signaling factors were scrutinized to ameliorate neural deficits in AD and stroke. Taken together, digging deeper in the common connections and signal targeting can be greatly helpful in understanding and unified treating of these disorders.

Genome-wide pathway-based quantitative multiple phenotypes analysis

For complex diseases, genome-wide pathway association studies have become increasingly promising. Currently, however, pathway-based association analysis mainly focus on a single phenotype, which may insufficient to describe the complex diseases and physiological processes. This work proposes a combination model to evaluate the association between a pathway and multiple phenotypes and to reduce the run time based on asymptotic results. For a single phenotype, we propose a semi-supervised maximum kernel-based U-statistics (mSKU) method to assess the pathway-based association analysis. For multiple phenotypes, we propose the fisher combination function with dependent phenotypes (FC) to transform the p-values between the pathway and each marginal phenotype individually to achieve pathway-based multiple phenotypes analysis. With real data from the Alzheimer Disease Neuroimaging Initiative (ADNI) study and Human Liver Cohort (HLC) study, the FC-mSKU method allows us to specify which pathways are specific to a single phenotype or contribute to common genetic constructions of multiple phenotypes. If we only focus on single-phenotype tests, we may miss some findings for etiology studies. Through extensive simulation studies, the FC-mSKU method demonstrates its advantages compared with its counterparts.

A network analysis to identify mediators of germline-driven differences in breast cancer prognosis

Identifying the underlying genetic drivers of the heritability of breast cancer prognosis remains elusive. We adapt a network-based approach to handle underpowered complex datasets to provide new insights into the potential function of germline variants in breast cancer prognosis. This network-based analysis studies ~7.3 million variants in 84,457 breast cancer patients in relation to breast cancer survival and confirms the results on 12,381 independent patients. Aggregating the prognostic effects of genetic variants across multiple genes, we identify four gene modules associated with survival in estrogen receptor (ER)-negative and one in ER-positive disease. The modules show biological enrichment for cancer-related processes such as G-alpha signaling, circadian clock, angiogenesis, and Rho-GTPases in apoptosis.

Shared Biological Pathways Between Alzheimer's Disease and Ischemic Stroke

Alzheimer's disease (AD) and ischemic stroke (IS) are an immense socioeconomic burden worldwide. There is a possibility that shared genetic factors lead to their links at epidemiological and pathophysiological levels. Although recent genome-wide association studies (GWAS) have provided profound insights into the genetics of AD and IS, no shared genetic variants have been identified to date. This prompted us to initiate this study, which sought to identify shared pathways linking AD and IS. We took advantage of large-scale GWAS summary data of AD (17,008 AD cases and 37,154 controls) and IS (10,307 cases and 19,326 controls) to conduct pathway analyses using genetic pathways from multiple well-studied databases, including GO, KEGG, PANTHER, Reactome, and Wikipathways. Collectively, we discovered that AD and IS shared 179 GO categories (56 biological processes, 95 cellular components, and 28 molecular functions); and the following pathways: six KEGG pathways; two PANTHER pathways; four Reactome pathways; and one in Wikipathways pathway. The more fine-grained GO terms were mainly summarized into different functional categories: transcriptional and post-transcriptional regulation, synapse, endocytic membrane traffic through the endosomal system, signaling transduction, immune process, multi-organism process, protein catabolic metabolism, and cell adhesion. The shared pathways were roughly classified into three categories: immune system; cancer (NSCLC and glioma); and signal transduction pathways involving the cadherin signaling pathway, Wnt signaling pathway, G-protein signaling and downstream signaling mediated by phosphoinositides (PIPs). The majority of these common pathways linked to both AD and IS were supported by convincing evidence from the literature. In conclusion, our findings contribute to a better understanding of common biological mechanisms underlying AD and IS and serve as a guide to direct future research.

Locus and gene-based GWAS meta-analysis identifies new diabetic nephropathy genes

Objective Assimilation of SNPs Interacting in Synchrony (OASIS) is a locus-based clustering algorithm recently described that can potentially address false positives and negatives in genome-wide association studies (GWAS) of complex disorders. Diabetic nephropathy (DN) is incompletely understood due to a paucity of genes identified despite several GWAS. OASIS was applied to three DN dbGAP GWAS datasets (4725 subjects; 1.06 million SNPs). OASIS identified 19 DN genes which were verified using single variant replication in a standard association study and gene-based analysis using GATES. CARS and FRMD3 were confirmed as DN genes, and five known diabetes-associated genes, viz. NLRP3, INPPL1, PIK3C2G, NRXN3, and TBC1D4, not previously identified using these datasets were discovered. Furthermore, three additional novel DN genes were found which replicated in two sets of analysis, viz. NTN1, EBF2, and DNAH11. Hence, composite analysis with OASIS, gene-based, and single variant association testing can be universally applied to existing GWAS datasets for the identification of new genes.

Locus and gene-based GWAS meta-analysis identifies new diabetic nephropathy genes

Objective Assimilation of SNPs Interacting in Synchrony (OASIS) is a locus-based clustering algorithm recently described that can potentially address false positives and negatives in genome-wide association studies (GWAS) of complex disorders. Diabetic nephropathy (DN) is incompletely understood due to a paucity of genes identified despite several GWAS. OASIS was applied to three DN dbGAP GWAS datasets (4725 subjects; 1.06 million SNPs). OASIS identified 19 DN genes which were verified using single variant replication in a standard association study and gene-based analysis using GATES. CARS and FRMD3 were confirmed as DN genes, and five known diabetes-associated genes, viz. NLRP3, INPPL1, PIK3C2G, NRXN3, and TBC1D4, not previously identified using these datasets were discovered. Furthermore, three additional novel DN genes were found which replicated in two sets of analysis, viz. NTN1, EBF2, and DNAH11. Hence, composite analysis with OASIS, gene-based, and single variant association testing can be universally applied to existing GWAS datasets for the identification of new genes.

OPATs: Omnibus P-value association tests

Combining statistical significances (P-values) from a set of single-locus association tests in genome-wide association studies is a proof-of-principle method for identifying disease-associated genomic segments, functional genes and biological pathways. We review P-value combinations for genome-wide association studies and introduce an integrated analysis tool, Omnibus P-value Association Tests (OPATs), which provides popular analysis methods of P-value combinations. The software OPATs programmed in R and R graphical user interface features a user-friendly interface. In addition to analysis modules for data quality control and single-locus association tests, OPATs provides three types of set-based association test: window-, gene- and biopathway-based association tests. P-value combinations with or without threshold and rank truncation are provided. The significance of a set-based association test is evaluated by using resampling procedures. Performance of the set-based association tests in OPATs has been evaluated by simulation studies and real data analyses. These set-based association tests help boost the statistical power, alleviate the multiple-testing problem, reduce the impact of genetic heterogeneity, increase the replication efficiency of association tests and facilitate the interpretation of association signals by streamlining the testing procedures and integrating the genetic effects of multiple variants in genomic regions of biological relevance. In summary, P-value combinations facilitate the identification of marker sets associated with disease susceptibility and uncover missing heritability in association studies, thereby establishing a foundation for the genetic dissection of complex diseases and traits. OPATs provides an easy-to-use and statistically powerful analysis tool for P-value combinations. OPATs, examples, and user guide can be downloaded from http://www.stat.sinica.edu.tw/hsinchou/genetics/association/OPATs.htm.

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.

MIDAS: Mining differentially activated subpaths of KEGG pathways from multi-class RNA-seq data

Pathway based analysis of high throughput transcriptome data is a widely used approach to investigate biological mechanisms. Since a pathway consists of multiple functions, the recent approach is to determine condition specific sub-pathways or subpaths. However, there are several challenges. First, few existing methods utilize explicit gene expression information from RNA-seq. More importantly, subpath activity is usually an average of statistical scores, e.g., correlations, of edges in a candidate subpath, which fails to reflect gene expression quantity information. In addition, none of existing methods can handle multiple phenotypes. To address these technical problems, we designed and implemented an algorithm, MIDAS, that determines condition specific subpaths, each of which has different activities across multiple phenotypes. MIDAS utilizes gene expression quantity information fully and the network centrality information to determine condition specific subpaths. To test performance of our tool, we used TCGA breast cancer RNA-seq gene expression profiles with five molecular subtypes. 36 differentially activate subpaths were determined. The utility of our method, MIDAS, was demonstrated in four ways. All 36 subpaths are well supported by the literature information. Subsequently, we showed that these subpaths had a good discriminant power for five cancer subtype classification and also had a prognostic power in terms of survival analysis. Finally, in a performance comparison of MIDAS to a recent subpath prediction method, PATHOME, our method identified more subpaths and much more genes that are well supported by the literature information.

AVAILABILITY

http://biohealth.snu.ac.kr/software/MIDAS/.

Gene-set meta-analysis of lung cancer identifies pathway related to systemic lupus erythematosus

INTRODUCTION

Gene-set analysis (GSA) is an approach using the results of single-marker genome-wide association studies when investigating pathways as a whole with respect to the genetic basis of a disease.

METHODS

We performed a meta-analysis of seven GSAs for lung cancer, applying the method META-GSA. Overall, the information taken from 11,365 cases and 22,505 controls from within the TRICL/ILCCO consortia was used to investigate a total of 234 pathways from the Kyoto Encyclopedia of Genes and Genomes (KEGG) database.

RESULTS

META-GSA reveals the systemic lupus erythematosus KEGG pathway hsa05322, driven by the gene region 6p21-22, as also implicated in lung cancer (p = 0.0306). This gene region is known to be associated with squamous cell lung carcinoma. The most important genes driving the significance of this pathway belong to the genomic areas HIST1-H4L, -1BN, -2BN, -H2AK, -H4K and C2/C4A/C4B. Within these areas, the markers most significantly associated with LC are rs13194781 (located within HIST12BN) and rs1270942 (located between C2 and C4A).

CONCLUSIONS

We have discovered a pathway currently marked as specific to systemic lupus erythematosus as being significantly implicated in lung cancer. The gene region 6p21-22 in this pathway appears to be more extensively associated with lung cancer than previously assumed. Given wide-stretched linkage disequilibrium to the area APOM/BAG6/MSH5, there is currently simply not enough information or evidence to conclude whether the potential pleiotropy of lung cancer and systemic lupus erythematosus is spurious, biological, or mediated. Further research into this pathway and gene region will be necessary.

Prioritizing biological pathways by recognizing context in time-series gene expression data

Background

The primary goal of pathway analysis using transcriptome data is to find significantly perturbed pathways. However, pathway analysis is not always successful in identifying pathways that are truly relevant to the context under study. A major reason for this difficulty is that a single gene is involved in multiple pathways. In the KEGG pathway database, there are 146 genes, each of which is involved in more than 20 pathways. Thus activation of even a single gene will result in activation of many pathways. This complex relationship often makes the pathway analysis very difficult. While we need much more powerful pathway analysis methods, a readily available alternative way is to incorporate the literature information.

Results

In this study, we propose a novel approach for prioritizing pathways by combining results from both pathway analysis tools and literature information. The basic idea is as follows. Whenever there are enough articles that provide evidence on which pathways are relevant to the context, we can be assured that the pathways are indeed related to the context, which is termed as relevance in this paper. However, if there are few or no articles reported, then we should rely on the results from the pathway analysis tools, which is termed as significance in this paper. We realized this concept as an algorithm by introducing Context Score and Impact Score and then combining the two into a single score. Our method ranked truly relevant pathways significantly higher than existing pathway analysis tools in experiments with two data sets.

Conclusions

Our novel framework was implemented as ContextTRAP by utilizing two existing tools, TRAP and BEST. ContextTRAP will be a useful tool for the pathway based analysis of gene expression data since the user can specify the context of the biological experiment in a set of keywords. The web version of ContextTRAP is available at http://biohealth.snu.ac.kr/software/contextTRAP.

Electronic supplementary material

The online version of this article (doi:10.1186/s12859-016-1335-8) contains supplementary material, which is available to authorized users.

Multiple linear combination (MLC) regression tests for common variants adapted to linkage disequilibrium structure

By jointly analyzing multiple variants within a gene, instead of one at a time, gene‐based multiple regression can improve power, robustness, and interpretation in genetic association analysis. We investigate multiple linear combination (MLC) test statistics for analysis of common variants under realistic trait models with linkage disequilibrium (LD) based on HapMap Asian haplotypes. MLC is a directional test that exploits LD structure in a gene to construct clusters of closely correlated variants recoded such that the majority of pairwise correlations are positive. It combines variant effects within the same cluster linearly, and aggregates cluster‐specific effects in a quadratic sum of squares and cross‐products, producing a test statistic with reduced degrees of freedom (df) equal to the number of clusters. By simulation studies of 1000 genes from across the genome, we demonstrate that MLC is a well‐powered and robust choice among existing methods across a broad range of gene structures. Compared to minimum P‐value, variance‐component, and principal‐component methods, the mean power of MLC is never much lower than that of other methods, and can be higher, particularly with multiple causal variants. Moreover, the variation in gene‐specific MLC test size and power across 1000 genes is less than that of other methods, suggesting it is a complementary approach for discovery in genome‐wide analysis. The cluster construction of the MLC test statistics helps reveal within‐gene LD structure, allowing interpretation of clustered variants as haplotypic effects, while multiple regression helps to distinguish direct and indirect associations.

Bayesian variable selection for post-analytic interrogation of susceptibility loci

Understanding the complex interplay among protein coding genes and regulatory elements requires rigorous interrogation with analytic tools designed for discerning the relative contributions of overlapping genomic regions. To this aim, we offer a novel application of Bayesian variable selection (BVS) for classifying genomic class level associations using existing large meta-analysis summary level resources. This approach is applied using the expectation maximization variable selection (EMVS) algorithm to typed and imputed SNPs across 502 protein coding genes (PCGs) and 220 long intergenic non-coding RNAs (lncRNAs) that overlap 45 known loci for coronary artery disease (CAD) using publicly available Global Lipids Gentics Consortium (GLGC) (Teslovich et al., 2010; Willer et al., 2013) meta-analysis summary statistics for low-density lipoprotein cholesterol (LDL-C). The analysis reveals 33 PCGs and three lncRNAs across 11 loci with >50% posterior probabilities for inclusion in an additive model of association. The findings are consistent with previous reports, while providing some new insight into the architecture of LDL-cholesterol to be investigated further. As genomic taxonomies continue to evolve, additional classes such as enhancer elements and splicing regions, can easily be layered into the proposed analysis framework. Moreover, application of this approach to alternative publicly available meta-analysis resources, or more generally as a post-analytic strategy to further interrogate regions that are identified through single point analysis, is straightforward. All coding examples are implemented in R version 3.2.1 and provided as supplemental material.

Genome-wide pathway-based association analysis identifies risk pathways associated with Parkinson's disease

Parkinson's disease (PD) is the second most common neurodegenerative disease. It is generally believed that it is influenced by both genetic and environmental factors, but the precise pathogenesis of PD is unknown to date. In this study, we performed a pathway analysis based on genome-wide association study (GWAS) to detect risk pathways of PD in three GWAS datasets. We first mapped all SNP markers to autosomal genes in each GWAS dataset. Then, we evaluated gene risk values using the minimum P-value of the tagSNPs. We took a pathway as a unit to identify the risk pathways based on the cumulative risks of the genes in the pathway. Finally, we combine the analysis results of the three datasets to detect the high risk pathways associated with PD. We found there were five same pathways in the three datasets. Besides, we also found there were five pathways which were shared in two datasets. Most of these pathways are associated with nervoussystem. Five pathways had been reported to be PD-related pathways in the previous literature. Our findings also implied that there was a close association between immune response and PD. Continued investigation of these pathways will further help us explain the pathogenesis of PD.

Embracing Integrative Multiomics Approaches

As “-omics” data technology advances and becomes more readily accessible to address complex biological questions, increasing amount of cross “-omics” dataset is inspiring the use and development of integrative bioinformatics analysis. In the current review, we discuss multiple options for integrating data across “-omes” for a range of study designs. We discuss established methods for such analysis and point the reader to in-depth discussions for the various topics. Additionally, we discuss challenges and new directions in the area.

How far from the SNP may the causative genes be?

While GWAS identify many disease-associated SNPs, using them to decipher disease mechanisms is hindered by the difficulty in mapping SNPs to genes. Most SNPs are in non-coding regions and it is often hard to identify the genes they implicate. To explore how far the SNP may be from the affected genes we used a pathway-based approach. We found that affected genes are often up to 2 Mbps away from the associated SNP, and are not necessarily the closest genes to the SNP. Existing approaches for mapping SNPs to genes leave many SNPs unmapped to genes and reveal only 86 significant phenotype-pathway associations for all known GWAS hits combined. Using the pathway-based approach we propose here allows mapping of virtually all SNPs to genes and reveals 435 statistically significant phenotype-pathway associations. In search for mechanisms that may explain the relationships between SNPs and distant genes, we found that SNPs that are mapped to distant genes have significantly more large insertions/deletions around them than other SNPs, suggesting that these SNPs may sometimes be markers for large insertions/deletions that may affect large genomic regions.

Pathway analysis of bladder cancer genome-wide association study identifies novel pathways involved in bladder cancer development

Genome-wide association studies (GWAS) are designed to identify individual regions associated with cancer risk, but only explain a small fraction of the inherited variability. Alternative approach analyzing genetic variants within biological pathways has been proposed to discover networks of susceptibility genes with additional effects. The gene set enrichment analysis (GSEA) may complement and expand traditional GWAS analysis to identify novel genes and pathways associated with bladder cancer risk. We selected three GSEA methods: Gen-Gen, Aligator, and the SNP Ratio Test to evaluate cellular signaling pathways involved in bladder cancer susceptibility in a Texas GWAS population. The candidate genetic polymorphisms from the significant pathway selected by GSEA were validated in an independent NCI GWAS. We identified 18 novel pathways (P < 0.05) significantly associated with bladder cancer risk. Five of the most promising pathways (P ≤ 0.001 in any of the three GSEA methods) among the 18 pathways included two cell cycle pathways and neural cell adhesion molecule (NCAM), platelet-derived growth factor (PDGF), and unfolded protein response pathways. We validated the candidate polymorphisms in the NCI GWAS and found variants of RAPGEF1, SKP1, HERPUD1, CACNB2, CACNA1C, CACNA1S, COL4A2, SRC, and CACNA1C were associated with bladder cancer risk. Two CCNE1 variants, rs8102137 and rs997669, from cell cycle pathways showed the strongest associations; the CCNE1 signal at 19q12 has already been reported in previous GWAS. These findings offer additional etiologic insights highlighting the specific genes and pathways associated with bladder cancer development. GSEA may be a complementary tool to GWAS to identify additional loci of cancer susceptibility.

Genome-wide interrogation reveals hundreds of long intergenic noncoding RNAs that associate with cardiometabolic traits

Long intergenic noncoding RNAs (lincRNAs) play important roles in disease, but the vast majority of these transcripts remain uncharacterized. We defined a set of 54 944 human lincRNAs by drawing on four publicly available lincRNA datasets, and annotated ∼2.5 million single nucleotide polymorphisms (SNPs) from each of 15 cardiometabolic genome-wide association study datasets into these lincRNAs. We identified hundreds of lincRNAs with at least one trait-associated SNP: 898 SNPs in 343 unique lincRNAs at 5% false discovery rate, and 469 SNPs in 146 unique lincRNAs meeting Bonferroni-corrected P < 0.05. An additional 64 trait-associated lincRNAs were identified using a class-level testing strategy at Bonferroni-corrected P < 0.05. To better understand the genomic context and prioritize trait-associated lincRNAs, we examined the pattern of linkage disequilibrium between SNPs in the lincRNAs and SNPs that met genome-wide-significance in the region (±500 kb of lincRNAs). A subset of the lincRNA-trait association findings was replicated in independent Genome-wide association studies data from the Pakistan Risk of Myocardial Infarction Study study. For trait-associated lincRNAs, we also investigated synteny and conservation relative to mouse, expression patterns in five cardiometabolic-relevant tissues, and allele-specific expression in RNA sequencing data for adipose tissue and leukocytes. Finally, we revealed a functional role in human adipocytes for linc-NFE2L3-1, which is expressed in adipose and is associated with waist-hip ratio adjusted for BMI. This comprehensive profile of trait-associated lincRNAs provides novel insights into disease mechanism and serves as a launching point for interrogation of the biology of specific lincRNAs in cardiometabolic disease.

A Simple Test of Class-Level Genetic Association Can Reveal Novel Cardiometabolic Trait Loci

BACKGROUND

Characterizing the genetic determinants of complex diseases can be further augmented by incorporating knowledge of underlying structure or classifications of the genome, such as newly developed mappings of protein-coding genes, epigenetic marks, enhancer elements and non-coding RNAs.

METHODS

We apply a simple class-level testing framework, termed Genetic Class Association Testing (GenCAT), to identify protein-coding gene association with 14 cardiometabolic (CMD) related traits across 6 publicly available genome wide association (GWA) meta-analysis data resources. GenCAT uses SNP-level meta-analysis test statistics across all SNPs within a class of elements, as well as the size of the class and its unique correlation structure, to determine if the class is statistically meaningful. The novelty of findings is evaluated through investigation of regional signals. A subset of findings are validated using recently updated, larger meta-analysis resources. A simulation study is presented to characterize overall performance with respect to power, control of family-wise error and computational efficiency. All analysis is performed using the GenCAT package, R version 3.2.1.

RESULTS

We demonstrate that class-level testing complements the common first stage minP approach that involves individual SNP-level testing followed by post-hoc ascribing of statistically significant SNPs to genes and loci. GenCAT suggests 54 protein-coding genes at 41 distinct loci for the 13 CMD traits investigated in the discovery analysis, that are beyond the discoveries of minP alone. An additional application to biological pathways demonstrates flexibility in defining genetic classes.

CONCLUSIONS

We conclude that it would be prudent to include class-level testing as standard practice in GWA analysis. GenCAT, for example, can be used as a simple, complementary and efficient strategy for class-level testing that leverages existing data resources, requires only summary level data in the form of test statistics, and adds significant value with respect to its potential for identifying multiple novel and clinically relevant trait associations.

Trade-offs in the effects of the apolipoprotein E polymorphism on risks of diseases of the heart, cancer, and neurodegenerative disorders: insights on mechanisms from the Long Life Family Study

The lack of evolutionary established mechanisms linking genes to age-related traits makes the problem of genetic susceptibility to health span inherently complex. One complicating factor is genetic trade-off. Here we focused on long-living participants of the Long Life Family Study (LLFS), their offspring, and spouses to: (1) Elucidate whether trade-offs in the effect of the apolipoprotein E e4 allele documented in the Framingham Heart Study (FHS) are a more general phenomenon, and (2) explore potential mechanisms generating age- and gender-specific trade-offs in the effect of the e4 allele on cancer, diseases of the heart, and neurodegenerative disorders assessed retrospectively in the LLFS populations. The e4 allele can diminish risks of cancer and diseases of the heart and confer risks of diseases of the heart in a sex-, age-, and LLFS-population-specific manner. A protective effect against cancer is seen in older long-living men and, potentially, their sons (>75 years, relative risk [RR]>75=0.48, p=0.086), which resembles our findings in the FHS. The protective effect against diseases of the heart is limited to long-living older men (RR>76=0.50, p=0.016), as well. A detrimental effect against diseases of the heart is characteristic for a normal LLFS population of male spouses and is specific for myocardial infarction (RR=3.07, p=2.1×10(-3)). These trade-offs are likely associated with two inherently different mechanisms, including disease-specific (detrimental; characteristic for a normal male population) and systemic, aging-related (protective; characteristic for older long-living men) mechanisms. The e4 allele confers risks of neurological disorders in men and women (RR=1.98, p=0.046). The results highlight the complex role of the e4 allele in genetic susceptibility to health span.

META-GSA: Combining Findings from Gene-Set Analyses across Several Genome-Wide Association Studies

INTRODUCTION

Gene-set analysis (GSA) methods are used as complementary approaches to genome-wide association studies (GWASs). The single marker association estimates of a predefined set of genes are either contrasted with those of all remaining genes or with a null non-associated background. To pool the p-values from several GSAs, it is important to take into account the concordance of the observed patterns resulting from single marker association point estimates across any given gene set. Here we propose an enhanced version of Fisher's inverse χ2-method META-GSA, however weighting each study to account for imperfect correlation between association patterns.

SIMULATION AND POWER

We investigated the performance of META-GSA by simulating GWASs with 500 cases and 500 controls at 100 diallelic markers in 20 different scenarios, simulating different relative risks between 1 and 1.5 in gene sets of 10 genes. Wilcoxon's rank sum test was applied as GSA for each study. We found that META-GSA has greater power to discover truly associated gene sets than simple pooling of the p-values, by e.g. 59% versus 37%, when the true relative risk for 5 of 10 genes was assume to be 1.5. Under the null hypothesis of no difference in the true association pattern between the gene set of interest and the set of remaining genes, the results of both approaches are almost uncorrelated. We recommend not relying on p-values alone when combining the results of independent GSAs.

APPLICATION

We applied META-GSA to pool the results of four case-control GWASs of lung cancer risk (Central European Study and Toronto/Lunenfeld-Tanenbaum Research Institute Study; German Lung Cancer Study and MD Anderson Cancer Center Study), which had already been analyzed separately with four different GSA methods (EASE; SLAT, mSUMSTAT and GenGen). This application revealed the pathway GO0015291 "transmembrane transporter activity" as significantly enriched with associated genes (GSA-method: EASE, p = 0.0315 corrected for multiple testing). Similar results were found for GO0015464 "acetylcholine receptor activity" but only when not corrected for multiple testing (all GSA-methods applied; p ≈ 0.02).

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.

Genome-wide association study of perioperative myocardial infarction after coronary artery bypass surgery

OBJECTIVES

Identification of patient subpopulations susceptible to develop myocardial infarction (MI) or, conversely, those displaying either intrinsic cardioprotective phenotypes or highly responsive to protective interventions remain high-priority knowledge gaps. We sought to identify novel common genetic variants associated with perioperative MI in patients undergoing coronary artery bypass grafting using genome-wide association methodology.

SETTING

107 secondary and tertiary cardiac surgery centres across the USA.

PARTICIPANTS

We conducted a stage I genome-wide association study (GWAS) in 1433 ethnically diverse patients of both genders (112 cases/1321 controls) from the Genetics of Myocardial Adverse Outcomes and Graft Failure (GeneMAGIC) study, and a stage II analysis in an expanded population of 2055 patients (225 cases/1830 controls) combined from the GeneMAGIC and Duke Perioperative Genetics and Safety Outcomes (PEGASUS) studies. Patients undergoing primary non-emergent coronary bypass grafting were included.

PRIMARY AND SECONDARY OUTCOME MEASURES

The primary outcome variable was perioperative MI, defined as creatine kinase MB isoenzyme (CK-MB) values ≥10× upper limit of normal during the first postoperative day, and not attributable to preoperative MI. Secondary outcomes included postoperative CK-MB as a quantitative trait, or a dichotomised phenotype based on extreme quartiles of the CK-MB distribution.

RESULTS

Following quality control and adjustment for clinical covariates, we identified 521 single nucleotide polymorphisms in the stage I GWAS analysis. Among these, 8 common variants in 3 genes or intergenic regions met p<10(-5) in stage II. A secondary analysis using CK-MB as a quantitative trait (minimum p=1.26×10(-3) for rs609418), or a dichotomised phenotype based on extreme CK-MB values (minimum p=7.72×10(-6) for rs4834703) supported these findings. Pathway analysis revealed that genes harbouring top-scoring variants cluster in pathways of biological relevance to extracellular matrix remodelling, endoplasmic reticulum-to-Golgi transport and inflammation.

CONCLUSIONS

Using a two-stage GWAS and pathway analysis, we identified and prioritised several potential susceptibility loci for perioperative MI.

The oxytocin receptor gene (OXTR) is associated with autism spectrum disorder: a meta-analysis

The oxytocin receptor gene (OXTR) has been studied as a risk factor for autism spectrum disorder (ASD) owing to converging evidence from multiple levels of analysis that oxytocin (OXT) has an important role in the regulation of affiliative behavior and social bonding in both nonhuman mammals and humans. Inconsistency in the effect sizes of the OXTR variants included in association studies render it unclear whether OXTR is truly associated with ASD, and, if so, which OXTR single-nucleotide polymorphisms (SNPs) are associated. Thus, a meta-analytic review of extant studies is needed to determine whether OXTR shows association with ASD, and to elucidate which specific SNPs have a significant effect on ASD. The current meta-analysis of 16 OXTR SNPs included 3941 individuals with ASD from 11 independent samples, although analyses of each individual SNP included a subset of this total. We found significant associations between ASD and the SNPs rs7632287, rs237887, rs2268491 and rs2254298. OXTR was also significantly associated with ASD in a gene-based test. The current meta-analysis is the largest and most comprehensive investigation of the association of OXTR with ASD and the findings suggest directions for future studies of the etiology of ASD.

Genotyping in rheumatoid arthritis: a game changer in clinical management?

Rheumatoid arthritis (RA) is a genetically complex disease of immune dysregulation characterized by painful inflammation of synovial joints. Despite advances in its management afforded by biologic drug development, efforts to improve outcomes for patients are confounded by the condition's heterogeneous pathobiology, and consequent variability in therapeutic responses. Great strides have been made in understanding the genetic epidemiology of rheumatoid arthritis since its association with the HLA locus was established in the 1980s, with over 100 additional disease-associated variants now confirmed through cumulative genome-wide association studies. Yet translation of this new knowledge for patient benefit - whether as a route to predicting disease risk, drug development or personalized medicine - has been slow. To address this, collaborating teams of interdisciplinary scientists will need to pool resources, including ever larger, well-characterized patient cohorts and sophisticated biostatistical approaches. Recent advances suggest that the fruits of these endeavors are beginning to come within reach.

Unravelling genes and pathways implicated in working memory of schizophrenia in Han Chinese

Working memory deficit is the core neurocognitive disorder in schizophrenia patients. To identify the factors underlying working memory deficit in schizophrenia patients and to explore the implication of possible genes in the working memory using genome-wide association study (GWAS) of schizophrenia, computerized delay-matching-to-sample (DMS) and whole genome genotyping data were obtained from 100 first-episode, treatment-naïve patients with schizophrenia and 140 healthy controls from the Mental Health Centre of the West China Hospital, Sichuan University. A composite score, delay-matching-to-sample total correct numbers (DMS-TC), was found to be significantly different between the patients and control. On associating quantitative DMS-TC with interactive variables of groups × genotype, one SNP (rs1411832), located downstream of YWHAZP5 in chromosome 10, was found to be associated with the working memory deficit in schizophrenia patients with lowest p-value (p = 2.02 × 10(-7)). ConsensusPathDB identified that genes with SNPs for which p values below the threshold of 5 × 10(-5) were significantly enriched in GO:0007155 (cell adhesion, p < 0.001). This study indicates that working memory, as an endophenotype of schizophrenia, could improve the efficacy of GWAS in schizophrenia. However, further study is required to replicate the results from our study.

Gene-environment interactions and obesity: recent developments and future directions

Obesity, a major public health concern, is a multifactorial disease caused by both environmental and genetic factors. Although recent genome-wide association studies have identified many loci related to obesity or body mass index, the identified variants explain only a small proportion of the heritability of obesity. Better understanding of the interplay between genetic and environmental factors is the basis for developing effective personalized obesity prevention and management strategies. This article reviews recent advances in identifying gene-environment interactions related to obesity and describes epidemiological designs and newly developed statistical approaches to characterizing and discovering gene-environment interactions on obesity risk.

Pathway-based analysis tools for complex diseases: a review

Genetic studies are traditionally based on single-gene analysis. The use of these analyses can pose tremendous challenges for elucidating complicated genetic interplays involved in complex human diseases. Modern pathway-based analysis provides a technique, which allows a comprehensive understanding of the molecular mechanisms underlying complex diseases. Extensive studies utilizing the methods and applications for pathway-based analysis have significantly advanced our capacity to explore large-scale omics data, which has rapidly accumulated in biomedical fields. This article is a comprehensive review of the pathway-based analysis methods—the powerful methods with the potential to uncover the biological depths of the complex diseases. The general concepts and procedures for the pathway-based analysis methods are introduced and then, a comprehensive review of the major approaches for this analysis is presented. In addition, a list of available pathway-based analysis software and databases is provided. Finally, future directions and challenges for the methodological development and applications of pathway-based analysis techniques are discussed. This review will provide a useful guide to dissect complex diseases.

Oesophageal squamous cell carcinoma in high-risk Chinese populations: Possible role for vascular epithelial growth factor A

BACKGROUND

Mechanisms involved in wound healing play some role in carcinogenesis in multiple organs, likely by creating a chronic inflammatory milieu. This study sought to assess the role of genetic markers in selected inflammation-related genes involved in wound healing (interleukin (IL)-1a, IL-1b, IL-1 Receptor type I (IL-1Ra), IL-1 Receptor type II (IL-1Rb), tumour necrosis factor (TNF)-α, tumour necrosis factor receptor superfamily member (TNFRSF)1A, nuclear factor kappa beta (NF-kB)1, NF-kB2, inducible nitric oxide synthase (iNOS), cyclooxygenase (COX)-2, hypoxia induced factor (HIF)-1α, vascular endothelial growth factor (VEGF)A and P-53) in risk to oesophageal squamous cell carcinoma (OSCC).

METHODS

We genotyped 125 tag single nucleotide polymorphism (SNP)s in 410 cases and 377 age and sex matched disease-free individuals from Nutritional Intervention Trial (NIT) cohort, and 546 cases and 556 controls individually matched for age, sex and neighbourhood from Shanxi case-control study, both conducted in high-risk areas of north-central China (1985-2007). Cox proportional-hazard models and conditional logistic regression models were used for SNPs analyses for NIT and Shanxi, respectively. Fisher's inverse test statistics were used to obtain gene-level significance.

RESULTS

Multiple SNPs were significantly associated with OSCC in both studies, however, none retained their significance after a conservative Bonferroni adjustment. Empiric p-values for tag SNPs in VEGFA in NIT were highly concentrated in the lower tail of the distribution, suggesting this gene may be influencing risk. Permutation tests confirmed the significance of this pattern. At the gene level, VEGFA yielded an empiric significance (P=0.027) in NIT. We also observed some evidence for interaction between environmental factors and some VEGFA tag SNPs.

CONCLUSION

Our finding adds further evidence for a potential role for markers in the VEGFA gene in the development and progression of early precancerous lesions of oesophagus.

Functional and genomic context in pathway analysis of GWAS data

Gene set analysis (GSA) is a promising tool for uncovering the polygenic effects associated with complex diseases. However, the available techniques reflect a wide variety of hypotheses about how genetic effects interact to contribute to disease susceptibility. The lack of consensus about the best way to perform GSA has led to confusion in the field and has made it difficult to compare results across methods. A clear understanding of the various choices made during GSA - such as how gene sets are defined, how single-nucleotide polymorphisms (SNPs) are assigned to genes, and how individual SNP-level effects are aggregated to produce gene- or pathway-level effects - will improve the interpretability and comparability of results across methods and studies. In this review we provide an overview of the various data sources used to construct gene sets and the statistical methods used to test for gene set association, as well as provide guidelines for ensuring the comparability of results.

Pathway analysis with next-generation sequencing data

Although pathway analysis methods have been developed and successfully applied to association studies of common variants, the statistical methods for pathway-based association analysis of rare variants have not been well developed. Many investigators observed highly inflated false-positive rates and low power in pathway-based tests of association of rare variants. The inflated false-positive rates and low true-positive rates of the current methods are mainly due to their lack of ability to account for gametic phase disequilibrium. To overcome these serious limitations, we develop a novel statistic that is based on the smoothed functional principal component analysis (SFPCA) for pathway association tests with next-generation sequencing data. The developed statistic has the ability to capture position-level variant information and account for gametic phase disequilibrium. By intensive simulations, we demonstrate that the SFPCA-based statistic for testing pathway association with either rare or common or both rare and common variants has the correct type 1 error rates. Also the power of the SFPCA-based statistic and 22 additional existing statistics are evaluated. We found that the SFPCA-based statistic has a much higher power than other existing statistics in all the scenarios considered. To further evaluate its performance, the SFPCA-based statistic is applied to pathway analysis of exome sequencing data in the early-onset myocardial infarction (EOMI) project. We identify three pathways significantly associated with EOMI after the Bonferroni correction. In addition, our preliminary results show that the SFPCA-based statistic has much smaller P-values to identify pathway association than other existing methods.

Integrative genomic analysis for the discovery of biomarkers in prostate cancer

Genome-wide association studies (GWAS) have achieved great success in identifying single nucleotide polymorphisms (SNPs, herein called genetic variants) and genes associated with risk of developing prostate cancer. However, GWAS do not typically link the genetic variants to the disease state or inform the broader context in which the genetic variants operate. Here, we present a novel integrative genomics approach that combines GWAS information with gene expression data to infer the causal association between gene expression and the disease and to identify the network states and biological pathways enriched for genetic variants. We identified gene regulatory networks and biological pathways enriched for genetic variants, including the prostate cancer, IGF-1, JAK2, androgen, and prolactin signaling pathways. The integration of GWAS information with gene expression data provides insights about the broader context in which genetic variants associated with an increased risk of developing prostate cancer operate.

A goodness-of-fit association test for whole genome sequencing data

Although many genetic factors have been successfully identified for human diseases in genome-wide association studies, genes discovered to date only account for a small proportion of overall genetic contributions to many complex traits. Association studies have difficulty in detecting the remaining true genetic variants that are either common variants with weak allelic effects, or rare variants that have strong allelic effects but are weakly associated at the population level. In this work, we applied a goodness-of-fit test for detecting sets of common and rare variants associated with quantitative or binary traits by using whole genome sequencing data. This test has been proved optimal for detecting weak and sparse signals in the literature, which fits the requirements for targeting the genetic components of missing heritability. Furthermore, this p value-combining method allows one to incorporate different data and/or research results for meta-analysis. The method was used to simultaneously analyse the whole genome sequencing and genome-wide association studies data of Genetic Analysis Workshop 18 for detecting true genetic variants. The results show that goodness-of-fit test is comparable or better than the influential sequence kernel association test in many cases.

Unraveling genetic origin of aging-related traits: evolving concepts

Discovering the genetic origin of aging-related traits could greatly advance strategies aiming to extend health span. The results of genome-wide association studies (GWAS) addressing this problem are controversial, and new genetic concepts have been fostered to advance the progress in the field. A limitation of GWAS and new genetic concepts is that they do not thoroughly address specifics of aging-related traits. Integration of theoretical concepts in genetics and aging research with empirical evidence from different disciplines highlights the conceptual problems in studies of genetic origin of aging-related traits. To address these problems, novel approaches of systemic nature are required. These approaches should adopt the non-deterministic nature of linkage of genes with aging-related traits and, consequently, reinforce research strategies for improving our understanding of mechanisms shaping genetic effects on these traits. Investigation of mechanisms will help determine conditions that activate specific genetic variants or profiles and explore to what extent these conditions that shape genetic effects are conserved across human lives and generations.

Supervised categorical principal component analysis for genome-wide association analyses

In order to have a better understanding of unexplained heritability for complex diseases in conventional Genome-Wide Association Studies (GWAS), aggregated association analyses based on predefined functional regions, such as genes and pathways, become popular recently as they enable evaluating joint effect of multiple Single-Nucleotide Polymorphisms (SNPs), which helps increase the detection power, especially when investigating genetic variants with weak individual effects. In this paper, we focus on aggregated analysis methods based on the idea of Principal Component Analysis (PCA). The past approaches using PCA mostly make some inherent genotype data and/or risk effect model assumptions, which may hinder the accurate detection of potential disease SNPs that influence disease phenotypes. In this paper, we derive a general Supervised Categorical Principal Component Analysis (SCPCA), which explicitly models categorical SNP data without imposing any risk effect model assumption. We have evaluated the efficacy of SCPCA with the comparison to a traditional Supervised PCA (SPCA) and a previously developed Supervised Logistic Principal Component Analysis (SLPCA) based on both the simulated genotype data by HAPGEN2 and the genotype data of Crohn's Disease (CD) from Wellcome Trust Case Control Consortium (WTCCC). Our preliminary results have demonstrated the superiority of SCPCA over both SPCA and SLPCA due to its modeling explicitly designed for categorical SNP data as well as its flexibility on the risk effect model assumption.

Whole-genome pathway analysis on 132,497 individuals identifies novel gene-sets associated with body mass index

Whole genome pathway analysis is a powerful tool for the exploration of the combined effects of gene-sets within biological pathways. This study applied Interval Based Enrichment Analysis (INRICH) to perform whole-genome pathway analysis of body-mass index (BMI). We used a discovery set composed of summary statistics from a meta-analysis of 123,865 subjects performed by the GIANT Consortium, and an independent sample of 8,632 subjects to assess replication of significant pathways. We examined SNPs within nominally significant pathways using linear mixed models to estimate their contribution to overall BMI heritability. Six pathways replicated as having significant enrichment for association after correcting for multiple testing, including the previously unknown relationships between BMI and the Reactome regulation of ornithine decarboxylase pathway, the KEGG lysosome pathway, and the Reactome stabilization of P53 pathway. Two non-overlapping sets of genes emerged from the six significant pathways. The clustering of shared genes based on previously identified protein-protein interactions listed in PubMed and OMIM supported the relatively independent biological effects of these two gene-sets. We estimate that the SNPs located in examined pathways explain ∼20% of the heritability for BMI that is tagged by common SNPs (3.35% of the 16.93% total).

Incorporating prior knowledge to increase the power of genome-wide association studies

Typical methods of analyzing genome-wide single nucleotide variant (SNV) data in cases and controls involve testing each variant's genotypes separately for phenotype association, and then using a substantial multiple-testing penalty to minimize the rate of false positives. This approach, however, can result in low power for modestly associated SNVs. Furthermore, simply looking at the most associated SNVs may not directly yield biological insights about disease etiology. SNVset methods attempt to address both limitations of the traditional approach by testing biologically meaningful sets of SNVs (e.g., genes or pathways). The number of tests run in a SNVset analysis is typically much lower (hundreds or thousands instead of millions) than in a traditional analysis, so the false-positive rate is lower. Additionally, by testing SNVsets that are biologically meaningful finding a significant set may more quickly yield insights into disease etiology.In this chapter we summarize the short history of SNVset testing and provide an overview of the many recently proposed methods. Furthermore, we provide detailed step-by-step instructions on how to perform a SNVset analysis, including a substantial number of practical tips and questions that researchers should consider before undertaking a SNVset analysis. Lastly, we describe a companion R package (snvset) that implements recently proposed SNVset methods. While SNVset testing is a new approach, with many new methods still being developed and many open questions, the promise of the approach is worth serious consideration when considering analytic methods for GWAS.

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.

Integrated enrichment analysis of variants and pathways in genome-wide association studies indicates central role for IL-2 signaling genes in type 1 diabetes, and cytokine signaling genes in Crohn's disease

Pathway analyses of genome-wide association studies aggregate information over sets of related genes, such as genes in common pathways, to identify gene sets that are enriched for variants associated with disease. We develop a model-based approach to pathway analysis, and apply this approach to data from the Wellcome Trust Case Control Consortium (WTCCC) studies. Our method offers several benefits over existing approaches. First, our method not only interrogates pathways for enrichment of disease associations, but also estimates the level of enrichment, which yields a coherent way to promote variants in enriched pathways, enhancing discovery of genes underlying disease. Second, our approach allows for multiple enriched pathways, a feature that leads to novel findings in two diseases where the major histocompatibility complex (MHC) is a major determinant of disease susceptibility. Third, by modeling disease as the combined effect of multiple markers, our method automatically accounts for linkage disequilibrium among variants. Interrogation of pathways from eight pathway databases yields strong support for enriched pathways, indicating links between Crohn's disease (CD) and cytokine-driven networks that modulate immune responses; between rheumatoid arthritis (RA) and "Measles" pathway genes involved in immune responses triggered by measles infection; and between type 1 diabetes (T1D) and IL2-mediated signaling genes. Prioritizing variants in these enriched pathways yields many additional putative disease associations compared to analyses without enrichment. For CD and RA, 7 of 8 additional non-MHC associations are corroborated by other studies, providing validation for our approach. For T1D, prioritization of IL-2 signaling genes yields strong evidence for 7 additional non-MHC candidate disease loci, as well as suggestive evidence for several more. Of the 7 strongest associations, 4 are validated by other studies, and 3 (near IL-2 signaling genes RAF1, MAPK14, and FYN) constitute novel putative T1D loci for further study.