Pathway analysis of complex diseases for GWAS, extending to consider rare variants, multi-omics and interactions

The association between single-nucleotide polymorphisms within type 1 interferon pathway genes and human immunodeficiency virus type 1 viral load in antiretroviral-naïve participants

BACKGROUND

Human genetic contribution to HIV progression remains inadequately explained. The type 1 interferon (IFN) pathway is important for host control of HIV and variation in type 1 IFN genes may contribute to disease progression. This study assessed the impact of variations at the gene and pathway level of type 1 IFN on HIV-1 viral load (VL).

METHODS

Two cohorts of antiretroviral (ART) naïve participants living with HIV (PLWH) with either early (START) or advanced infection (FIRST) were analysed separately. Type 1 IFN genes (n = 17) and receptor subunits (IFNAR1, IFNAR2) were examined for both cumulated type 1 IFN pathway analysis and individual gene analysis. SKAT-O was applied to detect associations between the genotype and HIV-1 study entry viral load (log10 transformed) as a proxy for set point VL; P-values were corrected using Bonferroni (P < 0.0025).

RESULTS

The analyses among those with early infection included 2429 individuals from five continents. The median study entry HIV VL was 14,623 (IQR 3460-45100) copies/mL. Across 673 SNPs within 19 type 1 IFN genes, no significant association with study entry VL was detected. Conversely, examining individual genes in START showed a borderline significant association between IFNW1, and study entry VL (P = 0.0025). This significance remained after separate adjustments for age, CD4+ T-cell count, CD4+/CD8+ T-cell ratio and recent infection. When controlling for population structure using linear mixed effects models (LME), in addition to principal components used in the main model, this was no longer significant (p = 0.0244). In subgroup analyses stratified by geographical region, the association between IFNW1 and study entry VL was only observed among African participants, although, the association was not significant when controlling for population structure using LME. Of the 17 SNPs within the IFNW1 region, only rs79876898 (A > G) was associated with study entry VL (p = 0.0020, beta = 0.32; G associated with higher study entry VL than A) in single SNP association analyses. The findings were not reproduced in FIRST participants.

CONCLUSION

Across 19 type 1 IFN genes, only IFNW1 was associated with HIV-1 study entry VL in a cohort of ART-naïve individuals in early stages of their infection, however, this was no longer significant in sensitivity analyses that controlled for population structures using LME.

Enhancing the anticancer potential of metformin: fabrication of efficient nanospanlastics, in vitro cytotoxic studies on HEP-2 cells and reactome enhanced pathway analysis

Metformin (MET), an oral antidiabetic drug, was reported to possess promising anticancer effects. We hypothesized that MET encapsulation in unique nanospanlastics would enhance its anticancer potential against HEP-2 cells. Our results showed the successful fabrication of Nano-MET spanlastics (d = 232.10 ± 0.20 nm; PDI = 0.25 ± 0.11; zeta potential = (-) 44.50 ± 0.96; drug content = 99.90 ± 0.11 and entrapment efficiency = 88.01 ± 2.50%). MTT assay revealed the enhanced Nano-MET cytotoxicity over MET with a calculated IC50 of 50 μg/mL and > 500 μg/mL, respectively. Annexin V/PI apoptosis assay showed that Nano-MET significantly decreased the percentage of live cells from 95.49 to 93.70 compared to MET and increased the percentage of cells arrested in the G0/G1 phase by 8.38%. Moreover, Nano-MET downregulated BCL-2 and upregulated BAX protein levels by 1.57 and 1.88 folds, respectively. RT-qPCR revealed that Nano-MET caused a significant 13.75, 4.15, and 2.23-fold increase in caspase-3, -8, and - 9 levels as well as a 100 and 43.47-fold decrease in cyclin D1 and mTOR levels, respectively. The proliferation marker Ki67 immunofluorescent staining revealed a 3-fold decrease in positive cells in Nano-MET compared to the control. Utilizing the combined Pathway-Enrichment Analysis (PEA) and Reactome analysis indicated high enrichment of certain pathways including nucleotides metabolism, Nudix-type hydrolase enzymes, carbon dioxide hydration, hemostasis, and the innate immune system. In summary, our results confirm MET cytotoxicity enhancement by its encapsulation in nanospanlastics. We also highlight, using PEA, that MET can modulate multiple pathways implicated in carcinogenesis.

Cauchy combination methods for the detection of gene-environment interactions for rare variants related to quantitative phenotypes

The characterization of gene-environment interactions (GEIs) can provide detailed insights into the biological mechanisms underlying complex diseases. Despite recent interest in GEIs for rare variants, published GEI tests are underpowered for an extremely small proportion of causal rare variants in a gene or a region. By extending the aggregated Cauchy association test (ACAT), we propose three GEI tests to address this issue: a Cauchy combination GEI test with fixed main effects (CCGEI-F), a Cauchy combination GEI test with random main effects (CCGEI-R), and an omnibus Cauchy combination GEI test (CCGEI-O). ACAT was applied to combine p values of single-variant GEI analyses to obtain CCGEI-F and CCGEI-R and p values of multiple GEI tests were combined in CCGEI-O. Through numerical simulations, for small numbers of causal variants, CCGEI-F, CCGEI-R and CCGEI-O provided approximately 5% higher power than the existing GEI tests INT-FIX and INT-RAN; however, they had slightly higher power than the existing GEI test TOW-GE. For large numbers of causal variants, although CCGEI-F and CCGEI-R exhibited comparable or slightly lower power values than the competing tests, the results were still satisfactory. Among all simulation conditions evaluated, CCGEI-O provided significantly higher power than that of competing GEI tests. We further applied our GEI tests in genome-wide analyses of systolic blood pressure or diastolic blood pressure to detect gene-body mass index (BMI) interactions, using whole-exome sequencing data from UK Biobank. At a suggestive significance level of 1.0 × 10-4, KCNC4, GAR1, FAM120AOS and NT5C3B showed interactions with BMI by our GEI tests.

Multi-omics research strategies in ischemic stroke: A multidimensional perspective

Ischemic stroke (IS) is a multifactorial and heterogeneous neurological disorder with high rate of death and long-term impairment. Despite years of studies, there are still no stroke biomarkers for clinical practice, and the molecular mechanisms of stroke remain largely unclear. The high-throughput omics approach provides new avenues for discovering biomarkers of IS and explaining its pathological mechanisms. However, single-omics approaches only provide a limited understanding of the biological pathways of diseases. The integration of multiple omics data means the simultaneous analysis of thousands of genes, RNAs, proteins and metabolites, revealing networks of interactions between multiple molecular levels. Integrated analysis of multi-omics approaches will provide helpful insights into stroke pathogenesis, therapeutic target identification and biomarker discovery. Here, we consider advances in genomics, transcriptomics, proteomics and metabolomics and outline their use in discovering the biomarkers and pathological mechanisms of IS. We then delineate strategies for achieving integration at the multi-omics level and discuss how integrative omics and systems biology can contribute to our understanding and management of IS.

Identifying genetic variants and pathways influencing daughter averages for twinning in North American Holstein cattle and evaluating the potential for genomic selection

Multiple birth in dairy cattle is a detrimental trait both economically for producers and for animal health. Genetics of twinning is complex and has led to several quantitative trait loci regions being associated with increased twinning. To identify variants associated with this trait, calving records from 2 time periods were used to estimate daughter averages for twinning for Holstein bulls. Multiple analyses were conducted and compared including GWAS, genomic prediction, and gene set enrichment analysis for pathway detection. Although pathway analysis did not yield many congruent pathways of interest between data sets, it did indicate two of interest. Both pathways have ties to the strong candidate region on BTA11 from the genome-wide association analysis across data sets. This region does not overlap with previously identified quantitative trait loci regions for twinning or ovulation rate in cattle. The strongest associated SNPs were upstream from 2 candidate genes LHCGR and FSHR, which are involved in folliculogenesis. Genomic prediction showed a moderate correlation accuracy (0.43) when predicting genomic breeding values for bulls with estimates from calving records from 2010 to 2016. Future analysis of the region on BTA11 and the relation of the candidate genes could improve this accuracy.

Integrative Pathway Analysis of SNP and Metabolite Data Using a Hierarchical Structural Component Model

Integrative multi-omics analysis has become a useful tool to understand molecular mechanisms and drug discovery for treatment. Especially, the couplings of genetics to metabolomics have been performed to identify the associations between SNP and metabolite. However, while the importance of integrative pathway analysis is increasing, there are few approaches to utilize pathway information to analyze phenotypes using SNP and metabolite. We propose an integrative pathway analysis of SNP and metabolite data using a hierarchical structural component model considering the structural relationships of SNPs, metabolites, pathways, and phenotypes. The proposed method utilizes genome-wide association studies on metabolites and constructs the genetic risk scores for metabolites referred to as genetic metabolomic scores. It is based on the hierarchical model using the genetic metabolomic scores and pathways. Furthermore, this method adopts a ridge penalty to consider the correlations between genetic metabolomic scores and between pathways. We apply our method to the SNP and metabolite data from the Korean population to identify pathways associated with type 2 diabetes (T2D). Through this application, we identified well-known pathways associated with T2D, demonstrating that this method adds biological insights into disease-related pathways using genetic predispositions of metabolites.

PlantPathMarks (PPMdb): an interactive hub for pathways-based markers in plant genomes

Over the past decade, the problem of finding an efficient gene-targeting marker set or signature for plant trait characterization has remained challenging. Many databases focusing on pathway mining have been released with one major deficiency, as they lack to develop marker sets that target only genes controlling a specific pathway or certain biological process. Herein, we present the PlantPathMarks database (PPMdb) as a comprehensive, web-based, user-friendly, and interactive hub for pathway-based markers in plant genomes. Based on our newly developed pathway gene set mining approach, two novel pathway-based marker systems called pathway gene-targeted markers (PGTMs) and pathway microsatellite-targeted markers (PMTMs) were developed as a novel class of annotation-based markers. In the PPMdb database, 2,690,742 pathway-based markers reflecting 9,894 marker panels were developed across 82 plant genomes. The markers include 691,555 PGTMs and 1,999,187 PMTMs. Across these genomes, 165,378 enzyme-coding genes were mapped against 126 KEGG reference pathway maps. PPMdb is furnished with three interactive visualization tools (Map Browse, JBrowse and Species Comparison) to visualize, map, and compare the developed markers over their KEGG reference pathway maps. All the stored marker panels can be freely downloaded. PPMdb promises to create a radical shift in the paradigm of the area of molecular marker research. The use of PPMdb as a mega-tool represents an impediment for non-bioinformatician plant scientists and breeders. PPMdb is freely available at http://ppmdb.easyomics.org.

A review on omics-based biomarkers discovery for Alzheimer's disease from the bioinformatics perspectives: Statistical approach vs machine learning approach

Alzheimer's Disease (AD) is a neurodegenerative disease that affects cognition and is the most common cause of dementia in the elderly. As the number of elderly individuals increases globally, the incidence and prevalence of AD are expected to increase. At present, AD is diagnosed clinically, according to accepted criteria. The essential elements in the diagnosis of AD include a patients history, a physical examination and neuropsychological testing, in addition to appropriate investigations such as neuroimaging. The omics-based approach is an emerging field of study that may not only aid in the diagnosis of AD but also facilitate the exploration of factors that influence the development of the disease. Omics techniques, including genomics, transcriptomics, proteomics and metabolomics, may reveal the pathways that lead to neuronal death and identify biomolecular markers associated with AD. This will further facilitate an understanding of AD neuropathology. In this review, omics-based approaches that were implemented in studies on AD were assessed from a bioinformatics perspective. Current state-of-the-art statistical and machine learning approaches used in the single omics analysis of AD were compared based on correlations of variants, differential expression, functional analysis and network analysis. This was followed by a review of the approaches used in the integration and analysis of multi-omics of AD. The strengths and limitations of multi-omics analysis methods were explored and the issues and challenges associated with omics studies of AD were highlighted. Lastly, future studies in this area of research were justified.

Variance-component-based meta-analysis of gene-environment interactions for rare variants

Complex diseases are often caused by interplay between genetic and environmental factors. Existing gene-environment interaction (G × E) tests for rare variants largely focus on detecting gene-based G × E effects in a single study; thus, their statistical power is limited by the sample size of the study. Meta-analysis methods that synthesize summary statistics of G × E effects from multiple studies for rare variants are still limited. Based on variance component models, we propose four meta-analysis methods of testing G × E effects for rare variants: HOM-INT-FIX, HET-INT-FIX, HOM-INT-RAN, and HET-INT-RAN. Our methods consider homogeneous or heterogeneous G × E effects across studies and treat the main genetic effect as either fixed or random. Through simulations, we show that the empirical distributions of the four meta-statistics under the null hypothesis align with their expected theoretical distributions. When the interaction effect is homogeneous across studies, HOM-INT-FIX and HOM-INT-RAN have as much statistical power as a pooled analysis conducted on a single interaction test with individual-level data from all studies. When the interaction effect is heterogeneous across studies, HET-INT-FIX and HET-INT-RAN provide higher power than pooled analysis. Our methods are further validated via testing 12 candidate gene-age interactions in blood pressure traits using whole-exome sequencing data from UK Biobank.

Distinct genetic regions are associated with differential population susceptibility to chemical exposures

Interactions between environmental factors and genetics underlie the majority of chronic human diseases. Chemical exposures are likely an underestimated contributor, yet gene-environment (GxE) interaction studies rarely assess their modifying effects. Here, we describe a novel method to profile the human genome and identify regions associated with differential population susceptibility to chemical exposures. Single nucleotide polymorphisms (SNPs) implicated in enriched chemical-disease intersections were identified and validated for three chemical classes with expected GxE interaction potential (neuroactive, hepatoactive, and cardioactive compounds). The same approach was then used to characterize consumer product classes with unknown risk for GxE interactions (washing products, cosmetics, and adhesives). Additionally, high-risk variant sets that may confer differential population susceptibility were identified for these consumer product groups through frequent itemset mining and pathway analysis. A dataset of 2454 consumer product chemical-disease linkages, with risk values, SNPs, and pathways for each association was developed, describing the interplay between environmental factors and genetics in human disease progression. We found that genetic hotspots implicated in GxE interactions differ across chemical classes (e.g., washing products had high-risk SNPs implicated in nervous system disease) and illustrate how this approach can discover new associations (e.g., washing product n-butoxyethanol implicated SNPs in the PI3K-Akt signaling pathway for Alzheimer's disease). Hence, our approach can predict high-risk genetic regions for differential population susceptibility to chemical exposures and characterize chemical modifying factors in specific diseases. These methods show promise for describing how chemical exposures can lead to varied health outcomes in a population and for incorporating inter-individual variability into chemical risk assessment.

The Genetics of Parkinson's Disease and Implications for Clinical Practice

The genetic landscape of Parkinson’s disease (PD) is characterised by rare high penetrance pathogenic variants causing familial disease, genetic risk factor variants driving PD risk in a significant minority in PD cases and high frequency, low penetrance variants, which contribute a small increase of the risk of developing sporadic PD. This knowledge has the potential to have a major impact in the clinical care of people with PD. We summarise these genetic influences and discuss the implications for therapeutics and clinical trial design.

Two novel oxetane containing lignans and a new megastigmane from Paronychia arabica and in silico analysis of them as prospective SARS-CoV-2 inhibitors

The chemical characterization of the extract of the aerial parts of Paronychia arabica afforded two oxetane containing lignans, paronychiarabicine A (1) and B (2), and one new megastigmane, paronychiarabicastigmane A (3), alongside a known lignan (4), eight known phenolic compounds (5–12), one known elemene sesquiterpene (13) and one steroid glycoside (14). The chemical structures of the isolated compounds were constructed based upon the HRMS, 1D, and 2D-NMR results. The absolute configurations were established via NOESY experiments as well as experimental and TDDFT-calculated electronic circular dichroism (ECD). Utilizing molecular docking, the binding scores and modes of compounds 1–3 towards the SARS-CoV-2 main protease (M^pro), papain-like protease (PL^pro), and RNA-dependent RNA polymerase (RdRp) were revealed. Compound 3 exhibited a promising docking score (−9.8 kcal mol⁻¹) against SARS-CoV-2 M^pro by forming seven hydrogen bonds inside the active site with the key amino acids. The reactome pathway enrichment analysis revealed a correlation between the inhibition of GSK3 and GSK3B genes (identified as the main targets of megastigmane treatment) and significant inhibition of SARS-CoV-1 viral replication in infected Vero E6 cells. Our results manifest a novel understanding of genes, proteins and corresponding pathways against SARS-CoV-2 infection and could facilitate the identification and characterization of novel therapeutic targets as treatments of SARS-CoV-2 infection.

The hydromethanolic extract of Paronychia arabica aerial parts afforded two oxetane containing lignans, paronychiarabicine A (1) and B (2), and one new megastigmane, paronychiarabicastigmane A (3), alongside a known secondary metabolites (4–14).

Multi-omics approaches to improve malaria therapy

Malaria contributes to the most widespread infectious diseases worldwide. Even though current drugs are commercially available, the ever-increasing drug resistance problem by malaria parasites poses new challenges in malaria therapy. Hence, searching for efficient therapeutic strategies is of high priority in malaria control. In recent years, multi-omics technologies have been extensively applied to provide a more holistic view of functional principles and dynamics of biological mechanisms. We briefly review multi-omics technologies and focus on recent malaria progress conducted with the help of various omics methods. Then, we present up-to-date advances for multi-omics approaches in malaria. Next, we describe resistance phenomena to established antimalarial drugs and underlying mechanisms. Finally, we provide insight into novel multi-omics approaches, new drugs and vaccine developments and analyze current gaps in multi-omics research. Although multi-omics approaches have been successfully used in malaria studies, they are still limited. Many gaps need to be filled to bridge the gap between basic research and treatment of malaria patients. Multi-omics approaches will foster a better understanding of the molecular mechanisms of Plasmodium that are essential for the development of novel drugs and vaccines to fight this disastrous disease.

Emerging roles of rare and low-frequency genetic variants in type 1 diabetes mellitus

Type 1 diabetes mellitus (T1DM) is defined as an autoimmune disorder and has enormous complexity and heterogeneity. Although its precise pathogenic mechanisms are obscure, this disease is widely acknowledged to be precipitated by environmental factors in individuals with genetic susceptibility. To date, the known susceptibility loci, which have mostly been identified by genome-wide association studies, can explain 80%–85% of the heritability of T1DM. Researchers believe that at least a part of its missing genetic component is caused by undetected rare and low-frequency variants. Most common variants have only small to modest effect sizes, which increases the difficulty of dissecting their functions and restricts their potential clinical application. Intriguingly, many studies have indicated that rare and low-frequency variants have larger effect sizes and play more significant roles in susceptibility to common diseases, including T1DM, than common variants do. Therefore, better recognition of rare and low-frequency variants is beneficial for revealing the genetic architecture of T1DM and for providing new and potent therapeutic targets for this disease. Here, we will discuss existing challenges as well as the great significance of this field and review current knowledge of the contributions of rare and low-frequency variants to T1DM.

Machine and Deep Learning in Molecular and Genetic Aspects of Sleep Research

Epidemiological sleep research strives to identify the interactions and causal mechanisms by which sleep affects human health, and to design intervention strategies for improving sleep throughout the lifespan. These goals can be advanced by further focusing on the environmental and genetic etiology of sleep disorders, and by development of risk stratification algorithms, to identify people who are at risk or are affected by, sleep disorders. These studies rely on comprehensive sleep-related data which often contains complex multi-dimensional physiological and molecular measurements across multiple timepoints. Thus, sleep research is well-suited for the application of computational approaches that can handle high-dimensional data. Here, we survey recent advances in machine and deep learning together with the availability of large human cohort studies with sleep data that can jointly drive the next breakthroughs in the sleep-research field. We describe sleep-related data types and datasets, and present some of the tasks in the field that can be targets for algorithmic approaches, as well as the challenges and opportunities in pursuing them.

Mechanistic insight into the progressive retinal atrophy disease in dogs via pathway-based genome-wide association analysis

The retinal degenerative disease, progressive retinal atrophy (PRA) is a major reason of vision impairment in canine population. Canine PRA signifies an inherently dissimilar category of retinal dystrophies which has solid resemblances to human retinis pigmentosa. Even though much is known about the biology of PRA, the knowledge about the intricate connection among genetic loci, genes and pathways associated to this disease in dogs are still remain unknown. Therefore, we have performed a genome wide association study (GWAS) to identify susceptibility single nucleotide polymorphisms (SNPs) of PRA. The GWAS was performed using a case–control based association analysis method on PRA dataset of 129 dogs and 135,553 markers. Further, the gene-set and pathway analysis were conducted in this study. A total of 1,114 markers associations with PRA trait at p < 0.01 were extracted and mapped to 640 unique genes, and then selected significant (p < 0.05) enriched 35 gene ontology (GO) terms and 5 Kyoto Encyclopedia of Genes and Genomes (KEGG) pathways contain these genes. In particular, apoptosis process, homophilic cell adhesion, calcium ion binding, and endoplasmic reticulum GO terms as well as pathways related to focal adhesion, cyclic guanosine monophosphate)-protein kinase G signaling, and axon guidance were more likely associated to the PRA disease in dogs. These data could provide new insight for further research on identification of potential genes and causative pathways for PRA in dogs.

Identification of Candidate Genes and Pathways Associated with Obesity-Related Traits in Canines via Gene-Set Enrichment and Pathway-Based GWAS Analysis

The present study aimed to identify causative loci and genes enriched in pathways associated with canine obesity using a genome-wide association study (GWAS). The GWAS was first performed to identify candidate single-nucleotide polymorphisms (SNPs) associated with obesity and obesity-related traits including body weight and blood sugar in 18 different breeds of 153 dogs. A total of 10 and 2 SNPs were found to be significantly (p < 3.74 × 10^-7) associated with body weight and blood sugar, respectively. None of the SNPs were identified to be significantly associated with obesity trait. We subsequently followed up the GWAS analysis with gene-set enrichment and pathway analyses. A gene-set with 1057, 1409, and 1243 SNPs annotated to 449, 933 and 820 genes for obesity, body weight, and blood sugar, respectively was created by sub-setting the GWAS result at a threshold of p < 0.01 for the gene-set enrichment analysis. In total, 84 GO and 21 KEGG pathways for obesity, 114 GO and 44 KEGG pathways for blood sugar, 120 GO and 24 KEGG pathways for body weight were found to be enriched. Among the pathways and GO terms, we highlighted five enriched pathways (Wnt signaling pathway, adherens junction, pathways in cancer, axon guidance, and insulin secretion) and seven GO terms (fat cell differentiation, calcium ion binding, cytoplasm, nucleus, phospholipid transport, central nervous system development, and cell surface) that were found to be shared among all the traits. Our data provide insights into the genes and pathways associated with obesity and obesity-related traits.

The Promises and the Challenges of Integrating Multi-Omics and Systems Biology in Comparative Stress Biology

Comparative stress biology is inherently a systems biology approach with the goal of integrating the molecular, cellular, and physiological responses with fitness outcomes. In this way, the systems biology approach is expected to provide a holistic understanding of how different stressors result in different fitness outcomes, and how different individuals (or populations or species) respond to stressors differently. In this perceptive article, I focus on the use of multiple types of -omics data in stress biology. Targeting students and those researchers who are considering integrating -omics approaches in their comparative stress biology studies, I discuss the promise of the integration of these measures for furthering our holistic understanding of how organisms respond to different stressors. I also discuss the logistical and conceptual challenges encountered when working with -omics data and the current hurdles to fully utilize these data in studies of stress biology in non-model organisms.

A Bayesian Hierarchical Framework for Pathway Analysis in Genome-Wide Association Studies

BACKGROUND

Pathway analysis allows joint consideration of multiple SNPs belonging to multiple genes, which in turn belong to a biologically defined pathway. This type of analysis is usually more powerful than single-SNP analyses for detecting joint effects of variants in a pathway.

METHODS

We develop a Bayesian hierarchical model by fully modeling the 3-level hierarchy, namely, SNP-gene-pathway that is naturally inherent in the structure of the pathways, unlike the currently used ad hoc ways of combining such information. We model the effects at each level conditional on the effects of the levels preceding them within the generalized linear model framework. To deal with the high dimensionality, we regularize the regression coefficients through an appropriate choice of priors. The model is fit using a combination of iteratively weighted least squares and expectation-maximization algorithms to estimate the posterior modes and their standard errors. A normal approximation is used for inference.

RESULTS

We conduct simulations to study the proposed method and find that our method has higher power than some standard approaches in several settings for identifying pathways with multiple modest-sized variants. We illustrate the method by analyzing data from two genome-wide association studies on breast and renal cancers.

CONCLUSION

Our method can be helpful in detecting pathway association.

A gene prioritization method based on a swine multi-omics knowledgebase and a deep learning model

The analyses of multi-omics data have revealed candidate genes for objective traits. However, they are integrated poorly, especially in non-model organisms, and they pose a great challenge for prioritizing candidate genes for follow-up experimental verification. Here, we present a general convolutional neural network model that integrates multi-omics information to prioritize the candidate genes of objective traits. By applying this model to Sus scrofa, which is a non-model organism, but one of the most important livestock animals, the model precision was 72.9%, recall 73.5%, and F1-Measure 73.4%, demonstrating a good prediction performance compared with previous studies in Arabidopsis thaliana and Oryza sativa. Additionally, to facilitate the use of the model, we present ISwine (http://iswine.iomics.pro/), which is an online comprehensive knowledgebase in which we incorporated almost all the published swine multi-omics data. Overall, the results suggest that the deep learning strategy will greatly facilitate analyses of multi-omics integration in the future.

Yuhua Fu et al. develop a CNN model that integrates multi-omics information to prioritize candidate genes of objective traits. Their model performs well when applied to important livestock non-model animals like Sus scrofa. Finally, the authors present ISwine, an online comprehensive knowledgebase which includes all published swine omics data to facilitate the integration of heterogeneous data.

A unified method for rare variant analysis of gene-environment interactions

Advanced technology in whole-genome sequencing has offered the opportunity to comprehensively investigate the genetic contribution, particularly rare variants, to complex traits. Several region-based tests have been developed to jointly model the marginal effect of rare variants, but methods to detect gene-environment (GE) interactions are underdeveloped. Identifying the modification effects of environmental factors on genetic risk poses a considerable challenge. To tackle this challenge, we develop a method to detect GE interactions for rare variants using generalized linear mixed effect model. The proposed method can accommodate either binary or continuous traits in related or unrelated samples. Under this model, genetic main effects, GE interactions, and sample relatedness are modeled as random effects. We adopt a kernel-based method to leverage the joint information across rare variants and implement variance component score tests to reduce the computational burden. Our simulation studies of continuous and binary traits show that the proposed method maintains correct type I error rates and appropriate power under various scenarios, such as genotype main effects and GE interaction effects in opposite directions and varying the proportion of causal variants in the model. We apply our method in the Framingham Heart Study to test GE interaction of smoking on body mass index or overweight status and replicate the Cholinergic Receptor Nicotinic Beta 4 gene association reported in previous large consortium meta-analysis of single nucleotide polymorphism-smoking interaction. Our proposed set-based GE test is computationally efficient and is applicable to both binary and continuous phenotypes, while appropriately accounting for familial or cryptic relatedness.

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.

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.

HisCoM-PAGE: Hierarchical Structural Component Models for Pathway Analysis of Gene Expression Data

Although there have been several analyses for identifying cancer-associated pathways, based on gene expression data, most of these are based on single pathway analyses, and thus do not consider correlations between pathways. In this paper, we propose a hierarchical structural component model for pathway analysis of gene expression data (HisCoM-PAGE), which accounts for the hierarchical structure of genes and pathways, as well as the correlations among pathways. Specifically, HisCoM-PAGE focuses on the survival phenotype and identifies its associated pathways. Moreover, its application to real biological data analysis of pancreatic cancer data demonstrated that HisCoM-PAGE could successfully identify pathways associated with pancreatic cancer prognosis. Simulation studies comparing the performance of HisCoM-PAGE with other competing methods such as Gene Set Enrichment Analysis (GSEA), Global Test, and Wald-type Test showed HisCoM-PAGE to have the highest power to detect causal pathways in most simulation scenarios.

Myopia genetics in genome-wide association and post-genome-wide association study era

Genome-wide association studies (GWAS) of myopia and refractive error have generated exciting results and identified novel risk-associated loci. However, the interpretation of the findings of GWAS of complex diseases is not straightforward and has remained challenging. This review provides a brief summary of the main focus on the advantages and limitations of GWAS of myopia, with potential strategies that may contribute to further insight into the genetics of myopia in the post-GWAS or omics era.

Genome-wide association studies of severe P. falciparum malaria susceptibility: progress, pitfalls and prospects

BACKGROUND

P. falciparum malaria has been recognized as one of the prominent evolutionary selective forces of human genome that led to the emergence of multiple host protective alleles. A comprehensive understanding of the genetic bases of severe malaria susceptibility and resistance can potentially pave ways to the development of new therapeutics and vaccines. Genome-wide association studies (GWASs) have recently been implemented in malaria endemic areas and identified a number of novel association genetic variants. However, there are several open questions around heritability, epistatic interactions, genetic correlations and associated molecular pathways among others. Here, we assess the progress and pitfalls of severe malaria susceptibility GWASs and discuss the biology of the novel variants.

RESULTS

We obtained all severe malaria susceptibility GWASs published thus far and accessed GWAS dataset of Gambian populations from European Phenome Genome Archive (EGA) through the MalariaGen consortium standard data access protocols. We noticed that, while some of the well-known variants including HbS and ABO blood group were replicated across endemic populations, only few novel variants were convincingly identified and their biological functions remain to be understood. We estimated SNP-heritability of severe malaria at 20.1% in Gambian populations and showed how advanced statistical genetic analytic methods can potentially be implemented in malaria susceptibility studies to provide useful functional insights.

CONCLUSIONS

The ultimate goal of malaria susceptibility study is to discover a novel causal biological pathway that provide protections against severe malaria; a fundamental step towards translational medicine such as development of vaccine and new therapeutics. Beyond singe locus analysis, the future direction of malaria susceptibility requires a paradigm shift from single -omics to multi-stage and multi-dimensional integrative functional studies that combines multiple data types from the human host, the parasite, the mosquitoes and the environment. The current biotechnological and statistical advances may eventually lead to the feasibility of systems biology studies and revolutionize malaria research.

Genome-wide association study of germline variants and breast cancer-specific mortality

BACKGROUND

We examined the associations between germline variants and breast cancer mortality using a large meta-analysis of women of European ancestry.

METHODS

Meta-analyses included summary estimates based on Cox models of twelve datasets using ~10.4 million variants for 96,661 women with breast cancer and 7697 events (breast cancer-specific deaths). Oestrogen receptor (ER)-specific analyses were based on 64,171 ER-positive (4116) and 16,172 ER-negative (2125) patients. We evaluated the probability of a signal to be a true positive using the Bayesian false discovery probability (BFDP).

RESULTS

We did not find any variant associated with breast cancer-specific mortality at P < 5 × 10-8. For ER-positive disease, the most significantly associated variant was chr7:rs4717568 (BFDP = 7%, P = 1.28 × 10-7, hazard ratio [HR] = 0.88, 95% confidence interval [CI] = 0.84-0.92); the closest gene is AUTS2. For ER-negative disease, the most significant variant was chr7:rs67918676 (BFDP = 11%, P = 1.38 × 10-7, HR = 1.27, 95% CI = 1.16-1.39); located within a long intergenic non-coding RNA gene (AC004009.3), close to the HOXA gene cluster.

CONCLUSIONS

We uncovered germline variants on chromosome 7 at BFDP < 15% close to genes for which there is biological evidence related to breast cancer outcome. However, the paucity of variants associated with mortality at genome-wide significance underpins the challenge in providing genetic-based individualised prognostic information for breast cancer patients.

Benefits and Challenges of Rare Genetic Variation in Alzheimer's Disease

Purpose of Review

It is well established that sporadic Alzheimer's disease (AD) is polygenic with common and rare genetic variation alongside environmental factors contributing to disease. Here, we review our current understanding of the genetic architecture of disease, paying specific attention to rare susceptibility variants, and explore some of the limitations in rare variant detection and analysis.

Recent Findings

Rare variation has been shown to robustly associate with disease. These include potentially damaging and loss of function mutations that are easily modelled in silico, in vitro and in vivo, and represent potentially druggable targets. A number of risk genes, including TREM2, SORL1 and ABCA7 show multiple independent associations suggesting that they may influence disease via multiple mechanisms. With transcriptional regulation, inflammatory response and modification of protein production suggested to be of primary importance.

Summary

We are at the beginning of our journey of rare variant detection in AD. Whole exome sequencing has been the predominant technology of choice. While fruitful, this has introduced a number of challenges with regard to data integration. Ultimately the future of disease-associated rare variant identification lies in whole genome sequencing projects that will allow the testing of the full range of genomic variation.

Using Openly Accessible Resources to Strengthen Causal Inference in Epigenetic Epidemiology of Neurodevelopment and Mental Health

The recent focus on the role of epigenetic mechanisms in mental health has led to several studies examining the association of epigenetic processes with psychiatric conditions and neurodevelopmental traits. Some studies suggest that epigenetic changes might be causal in the development of the psychiatric condition under investigation. However, other scenarios are possible, e.g., statistical confounding or reverse causation, making it particularly challenging to derive conclusions on causality. In the present review, we examine the evidence from human population studies for a possible role of epigenetic mechanisms in neurodevelopment and mental health and discuss methodological approaches on how to strengthen causal inference, including the need for replication, (quasi-)experimental approaches and Mendelian randomization. We signpost openly accessible resources (e.g., "MR-Base" "EWAS catalog" as well as tissue-specific methylation and gene expression databases) to aid the application of these approaches.

Efficient Signal Inclusion With Genomic Applications

This paper addresses the challenge of efficiently capturing a high proportion of true signals for subsequent data analyses when sample sizes are relatively limited with respect to data dimension. We propose the signal missing rate as a new measure for false negative control to account for the variability of false negative proportion. Novel data-adaptive procedures are developed to control signal missing rate without incurring many unnecessary false positives under dependence. We justify the efficiency and adaptivity of the proposed methods via theory and simulation. The proposed methods are applied to GWAS on human height to effectively remove irrelevant SNPs while retaining a high proportion of relevant SNPs for subsequent polygenic analysis.

Genetic risk factors for eating disorders: an update and insights into pathophysiology

Genome-wide-association studies (GWASs), epigenetic, gene-expression and gene-gene interaction projects, nutritional genomics and investigations of the gut microbiota have increased our knowledge of the pathophysiology of eating disorders (EDs). However, compared with anorexia nervosa, genetic studies in patients with bulimia nervosa and binge-eating disorder are relatively scarce, with the exception of a few formal genetic and small-sized candidate-gene-association studies. In this article, we review important findings derived from formal and molecular genetics in order to outline a genetics-based pathophysiological model of EDs. This model takes into account environmental and nutritional factors, genetic factors related to the microbiome, the metabolic and endocrine system, the immune system, and the brain, in addition to phenotypical traits of EDs. Shortcomings and advantages of genetic research in EDs are discussed against the historical background, but also in light of potential future treatment options for patients with EDs.

GWAS of Behavioral Traits

Over the past decade, genome-wide association studies (GWAS) have evolved into a powerful tool to investigate genetic risk factors for human diseases via a hypothesis-free scan of the genome. The success of GWAS for psychiatric disorders and behavioral traits have been somewhat mixed, partly owing to the complexity and heterogeneity of these traits. Significant progress has been made in the last few years in the development and implementation of complex statistical methods and algorithms incorporating GWAS. Such advanced statistical methods applied to GWAS hits in combination with incorporation of different layers of genomics data have catapulted the search for novel genes for behavioral traits and improved our understanding of the complex polygenic architecture of these traits.This chapter will give a brief overview on GWAS and statistical methods currently used in GWAS. The chapter will focus on reviewing the current literature and highlight some of the most important GWAS on psychiatric and other behavioral traits and will conclude with a discussion on future directions.

Sepsis: Personalized Medicine Utilizing 'Omic' Technologies-A Paradigm Shift?

Sepsis has over the years proven a considerable challenge to physicians and researchers. Numerous pharmacological and non-pharmacological interventions have been tested in trials, but have unfortunately failed to improve the general prognosis. This has led to the speculation that the sepsis population may be too heterogeneous to be targeted with the traditional one treatment suits all’ approach. Recent advances in genetic and biochemical analyses now allow genotyping and biochemical characterisation of large groups of patients via the ‘omics’ technologies. These new opportunities could lead to a paradigm shift in the approach to sepsis towards personalised treatments with interventions targeted towards specific pathophysiological mechanisms activated in the patient. In this article, we review the potentials and pitfalls of using new advanced technologies to deepen our understanding of the clinical syndrome of sepsis.

Rigor and reproducibility in genetic research on eating disorders

OBJECTIVE

We explored both within-method and between-method rigor and reproducibility in the field of eating disorders genetics.

METHOD

We present critical evaluation and commentary on component methods of genetic research (family studies, twin studies, molecular genetic studies) and discuss both successful and unsuccessful efforts in the field.

RESULTS

Eating disorders genetics has had a number of robust results that converge across component methodologies. Familial aggregation of eating disorders, twin-based heritability estimates of eating disorders, and genome-wide association studies (GWAS) all point toward a substantial role for genetics in eating disorders etiology and support the premise that genes do not act alone. Candidate gene and linkage studies have been less informative historically.

DISCUSSION

The eating disorders field has entered the GWAS era with studies of anorexia nervosa. Continued growth of sample sizes is essential for rigorous discovery of actionable variation. Molecular genetic studies of bulimia nervosa, binge-eating disorder, and other eating disorders are virtually nonexistent and lag seriously behind other major psychiatric disorders. Expanded efforts are necessary to reveal the fundamental biology of eating disorders, inform clinical practice, and deliver new therapeutic targets.

From SNPs to pathways: Biological interpretation of type 2 diabetes (T2DM) genome wide association study (GWAS) results

Genome-wide association studies (GWAS) have become a common method for discovery of gene-disease relationships, in particular for complex diseases like Type 2 Diabetes Mellitus (T2DM). The experience with GWAS analysis has revealed that the genetic risk for complex diseases involves cumulative, small effects of many genes and only some genes with a moderate effect. In order to explore the complexity of the relationships between T2DM genes and their potential function at the process level as effected by polymorphism effects, a secondary analysis of a GWAS meta-analysis is presented. Network analysis, pathway information and integration of different types of biological information such as eQTLs and gene-environment interactions are used to elucidate the biological context of the genetic variants and to perform an analysis based on data visualization. We selected a T2DM dataset from a GWAS meta-analysis, and extracted 1,971 SNPs associated with T2DM. We mapped 580 SNPs to 360 genes, and then selected 460 pathways containing these genes from the curated collection of WikiPathways. We then created and analyzed SNP-gene and SNP-gene-pathway network modules in Cytoscape. A focus on genes with robust connections to pathways permitted identification of many T2DM pertinent pathways. However, numerous genes lack literature evidence of association with T2DM. We also speculate on the genes in specific network structures obtained in the SNP-gene network, such as gene-SNP-gene modules. Finally, we selected genes relevant to T2DM from our SNP-gene-pathway network, using different sources that reveal gene-environment interactions and eQTLs. We confirmed functions relevant to T2DM for many genes and have identified some-LPL and APOB-that require further validation to clarify their involvement in T2DM.

Circadian pathway genetic variation and cancer risk: evidence from genome-wide association studies

BACKGROUND

Dysfunction of the circadian clock and single polymorphisms of some circadian genes have been linked to cancer susceptibility, although data are scarce and findings inconsistent. We aimed to investigate the association between circadian pathway genetic variation and risk of developing common cancers based on the findings of genome-wide association studies (GWASs).

METHODS

Single nucleotide polymorphisms (SNPs) of 17 circadian genes reported by three GWAS meta-analyses dedicated to breast (Discovery, Biology, and Risk of Inherited Variants in Breast Cancer (DRIVE) Consortium; cases, n = 15,748; controls, n = 18,084), prostate (Elucidating Loci Involved in Prostate Cancer Susceptibility (ELLIPSE) Consortium; cases, n = 14,160; controls, n = 12,724) and lung carcinoma (Transdisciplinary Research In Cancer of the Lung (TRICL) Consortium; cases, n = 12,160; controls, n = 16,838) in patients of European ancestry were utilized to perform pathway analysis by means of the adaptive rank truncated product (ARTP) method. Data were also available for the following subgroups: estrogen receptor negative breast cancer, aggressive prostate cancer, squamous lung carcinoma and lung adenocarcinoma.

RESULTS

We found a highly significant statistical association between circadian pathway genetic variation and the risk of breast (pathway P value = 1.9 × 10^-6; top gene RORA, gene P value = 0.0003), prostate (pathway P value = 4.1 × 10^-6; top gene ARNTL, gene P value = 0.0002) and lung cancer (pathway P value = 6.9 × 10^-7; top gene RORA, gene P value = 2.0 × 10^-6), as well as all their subgroups. Out of 17 genes investigated, 15 were found to be significantly associated with the risk of cancer: four genes were shared by all three malignancies (ARNTL, CLOCK, RORA and RORB), two by breast and lung cancer (CRY1 and CRY2) and three by prostate and lung cancer (NPAS2, NR1D1 and PER3), whereas four genes were specific for lung cancer (ARNTL2, CSNK1E, NR1D2 and PER2) and two for breast cancer (PER1, RORC).

CONCLUSIONS

Our findings, based on the largest series ever utilized for ARTP-based gene and pathway analysis, support the hypothesis that circadian pathway genetic variation is involved in cancer predisposition.

Association study of genetic variation in the autophagy lysosome pathway genes and risk of eight kinds of cancers

The autophagy lysosome pathway is essential to maintain cell viability and homeostasis in response to many stressful environments, which is reported to play a vital role in cancer development and therapy. However, the association of genetic alterations of this pathway with risk of cancer remains unclear. Based on genome-wide association study data of eight kinds of cancers, we used an adaptive rank truncated product approach to perform a pathway-level and gene-level analysis, and used a logistic model to calculate SNP-level associations to examine whether an altered autophagy lysosome pathway contributes to cancer susceptibility. Among eight kinds of cancers, four of them showed significant statistics in the pathway-level analysis, including breast cancer (p = 0.00705), gastric cancer (p = 0.00880), lung cancer (p = 0.000100) and renal cell carcinoma (p = 0.00190). We also found that some autophagy lysosome genes had signals of association with cancer risk. Our results demonstrated that inherited genetic variants in the overall autophagy lysosome pathway and certain associated genes might contribute to cancer susceptibility, which warrant further evaluation in other independent datasets.

The Fas Signaling Pathway Is a Common Genetic Risk Factor for Severe Cutaneous Drug Adverse Reactions Across Diverse Drugs

PURPOSE

Human leukocyte antigen (HLA) has been recognized as the most important genetic risk factor for severe cutaneous adverse drug reactions (SCARs) caused by certain drugs. However, cumulated observations suggest the presence of genetic risk factors for SCARs other than drug-specific HLA. We aimed to identify a common genetic risk factor of SCARs across multiple drugs.

METHODS

We performed 2 independent genome-wide association studies (GWASs). A total of 68 and 38 subjects with a diagnosis of SCAR were enrolled in each GWAS. Their allele frequencies were compared to those of healthy subjects in Korea.

RESULTS

No single nucleotide polymorphism (SNP) with genome-wide significance was found in either GWAS. We next selected and annotated the 200 top-ranked SNPs from each GWAS. These 2 sets of annotated genes were then entered into the web interface of ConsensusPathDB for a pathway-level analysis. The Fas signaling pathway was significantly over-represented in each gene set from the 2 GWASs.

CONCLUSIONS

Our observations suggest that the Fas signaling pathway may be a common genetic risk factor for SCARs across multiple drugs.

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.