Functional and genomic context in pathway analysis of GWAS data

Polygenic regression uncovers trait-relevant cellular contexts through pathway activation transformation of single-cell RNA sequencing data

Advances in single-cell RNA sequencing (scRNA-seq) techniques have accelerated functional interpretation of disease-associated variants discovered from genome-wide association studies (GWASs). However, identification of trait-relevant cell populations is often impeded by inherent technical noise and high sparsity in scRNA-seq data. Here, we developed scPagwas, a computational approach that uncovers trait-relevant cellular context by integrating pathway activation transformation of scRNA-seq data and GWAS summary statistics. scPagwas effectively prioritizes trait-relevant genes, which facilitates identification of trait-relevant cell types/populations with high accuracy in extensive simulated and real datasets. Cellular-level association results identified a novel subpopulation of naive CD8+ T cells related to COVID-19 severity and oligodendrocyte progenitor cell and microglia subsets with critical pathways by which genetic variants influence Alzheimer's disease. Overall, our approach provides new insights for the discovery of trait-relevant cell types and improves the mechanistic understanding of disease variants from a pathway perspective.

Hierarchical association of COPD to principal genetic components of biological systems

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

Modern drug discovery for inflammatory bowel disease: The role of computational methods

Inflammatory bowel diseases (IBDs) comprising ulcerative colitis, Crohn’s disease and microscopic colitis are characterized by chronic inflammation of the gastrointestinal tract. IBD has spread around the world and is becoming more prevalent at an alarming rate in developing countries whose societies have become more westernized. Cell therapy, intestinal microecology, apheresis therapy, exosome therapy and small molecules are emerging therapeutic options for IBD. Currently, it is thought that low-molecular-mass substances with good oral bio-availability and the ability to permeate the cell membrane to regulate the action of elements of the inflammatory signaling pathway are effective therapeutic options for the treatment of IBD. Several small molecule inhibitors are being developed as a promising alternative for IBD therapy. The use of highly efficient and time-saving techniques, such as computational methods, is still a viable option for the development of these small molecule drugs. The computer-aided (in silico) discovery approach is one drug development technique that has mostly proven efficacy. Computational approaches when combined with traditional drug development methodology dramatically boost the likelihood of drug discovery in a sustainable and cost-effective manner. This review focuses on the modern drug discovery approaches for the design of novel IBD drugs with an emphasis on the role of computational methods. Some computational approaches to IBD genomic studies, target identification, and virtual screening for the discovery of new drugs and in the repurposing of existing drugs are discussed.

Boosting of tau protein aggregation by CD40 and CD48 gene expression in Alzheimer's disease

Neurodegenerative diseases result from the interplay of abnormal gene expression and various pathological factors. Therefore, a disease-specific integrative genetic approach is required to understand the complexities and causes of target diseases. Recent studies have identified the correlation between genes encoding several transmembrane proteins, such as the cluster of differentiation (CD) and Alzheimer's disease (AD) pathogenesis. In this study, CD48 and CD40 gene expression in AD, a neurodegenerative disease, was analyzed to infer this link. Total RNA sequencing was performed using an Alzheimer's disease mouse model brain and blood, and gene expression was determined using a genome-wide association study (GWAS). We observed a marked elevation of CD48 and CD40 genes in Alzheimer's disease. Indeed, the upregulation of both CD48 and CD40 genes was significantly increased in the severe Alzheimer's disease group. With the elevation of CD48 and CD40 genes in Alzheimer's disease, associations of protein levels were also markedly increased in tissues. In addition, overexpression of CD48 and CD40 genes triggered tau aggregation, and co-expression of these genes accelerated aggregation. The nuclear factor kappa B (NF-ĸB) signaling pathway was enriched by CD48 and CD40 gene expression: it was also associated with tau pathology. Our data suggested that the CD48 and CD40 genes are novel AD-related genes, and this approach may be useful as a diagnostic or therapeutic target for the disease.

Overview of methods for characterization and visualization of a protein–protein interaction network in a multi-omics integration context

At the heart of the cellular machinery through the regulation of cellular functions, protein–protein interactions (PPIs) have a significant role. PPIs can be analyzed with network approaches. Construction of a PPI network requires prediction of the interactions. All PPIs form a network. Different biases such as lack of data, recurrence of information, and false interactions make the network unstable. Integrated strategies allow solving these different challenges. These approaches have shown encouraging results for the understanding of molecular mechanisms, drug action mechanisms, and identification of target genes. In order to give more importance to an interaction, it is evaluated by different confidence scores. These scores allow the filtration of the network and thus facilitate the representation of the network, essential steps to the identification and understanding of molecular mechanisms. In this review, we will discuss the main computational methods for predicting PPI, including ones confirming an interaction as well as the integration of PPIs into a network, and we will discuss visualization of these complex data.

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.

Bioinspired NO release coating enhances endothelial cells and inhibits smooth muscle cells

Thrombus and restenosis after stent implantation are the major complications because traditional drugs such as rapamycin delay the process of endothelialization. Nitric oxide (NO) is mainly produced by endothelial nitric oxide synthase (eNOS) on the membrane of endothelial cells (ECs) in the cardiovascular system and plays an important role in vasomotor function. It strongly inhibits the proliferation of smooth muscle cells (SMCs) and ameliorates endothelial function when ECs get hurt. Inspired by this, introducing NO to traditional stent coating may alleviate endothelial insufficiency caused by rapamycin. Here, we introduced SNAP as the NO donor, mimicking how NO affects in vivo, into rapamycin coating to alleviate endothelial damage while inhibiting SMC proliferation. Through wicking effects, SNAP was absorbed into a hierarchical coating that had an upper porous layer and a dense polymer layer with rapamycin at the bottom. Cells were cultured on the coatings, and it was observed that the injured ECs were restored while the growth of SMCs further diminished. Genome analysis was conducted to further clarify possible signaling pathways: the effect of cell growth attenuated by NO may cause by affecting cell cycle and enhancing inflammation. These findings supported the idea that introducing NO to traditional drug-eluting stents alleviates incomplete endothelialization and further inhibits the stenosis caused by the proliferation of SMCs.

The misuse of distributional assumptions in functional class scoring gene-set and pathway analysis

Gene-set analysis (GSA) is a standard procedure for exploring potential biological functions of a group of genes. The development of its methodology has been an active research topic in recent decades. Many GSA methods, when newly proposed, rely on simulation studies to evaluate their performance with an implicit assumption that the multivariate expression values are normally distributed. This assumption is commonly adopted in GSAs, particularly those in the group of functional class scoring (FCS) methods. The validity of the normality assumption, however, has been disputed in several studies, yet no systematic analysis has been carried out to assess the effect of this distributional assumption. Our goal in this study is not to propose a new GSA method but to first examine if the multi-dimensional gene expression data in gene sets follow a multivariate normal (MVN) distribution. Six statistical methods in three categories of MVN tests were considered and applied to a total of 24 RNA data sets. These RNA values were collected from cancer patients as well as normal subjects, and the values were derived from microarray experiments, RNA sequencing, and single-cell RNA sequencing. Our first finding suggests that the MVN assumption is not always satisfied. This assumption does not hold true in many applications tested here. In the second part of this research, we evaluated the influence of non-normality on the statistical power of current FCS methods, both parametric and nonparametric ones. Specifically, the scenario of mixture distributions representing more than one population for the RNA values was considered. This second investigation demonstrates that the non-normality distribution of the RNA values causes a loss in the statistical power of these GSA tests, especially when subtypes exist. Among the FCS GSA tools examined here and among the scenarios studied in this research, the N-statistics outperform the others. Based on the results from these two investigations, we conclude that the assumption of MVN should be used with caution when evaluating new GSA tools, since this assumption cannot be guaranteed and violation may lead to spurious results, loss of power, and incorrect comparison between methods. If a newly proposed GSA tool is to be evaluated, we recommend the incorporation of a wide range of multivariate non-normal distributions or sampling from large databases if available.

Bioinformatics analysis of epigenetic and SNP-related molecular markers in systemic lupus erythematosus

We analyzed gene expression in peripheral blood mononuclear cells (PBMCs) from patients with systemic lupus erythematosus (SLE) using public databases. The goal was to identify lupus biomarkers by determining whether differentially expressed genes are mediated by methylation, miRNA, or SNP. Two cDNA microarrays were subjected to integration analysis, and we calculated the mutually differentially expressed genes (|log2fold change (FC)| > 1, P < 0.05). These genes were analyzed using gene otology (GO), Kyoto Encyclopedia of Genes and Genomes (KEGG) pathways, and protein-protein interaction (PPI) networks. The differences in methylation sites for two methylation chips were calculated and the differentially methylated sites were annotated. These genes were compared to the differentially expressed genes. We obtained 135 differentially expressed microRNAs from the microRNA-chip results using PBMCs from SLE and healthy individuals. Predictive microRNA target genes were identified using GO, KEGG pathways, and PPI networks. The target genes identified were compared to the differentially expressed genes. We downloaded Chinese SLE genome-wide association study data from SLE-related literature, analyzed the loci with a P value < 0.05, and used annotated SLE-associated SNPs. We selected the genes corresponding to an SNP located on an exon and determined the intersection with the differentially expressed genes. We found 18 differentially expressed genes in both cDNA microarrays. The methylation chips had 50 corresponding methylation sites. On the basis of these results, we identified two genes, IFI44 and IFI44L. We further identified 135 differentially expressed microRNAs predicted to affect 5766 target genes. Two identified genes were in common with the differentially expressed genes. Finally, SNP annotated genes and cDNA chip genes overlap with identified MX1. Therefore, we used existing data to analyze the causes of differential gene expression in SLE, introducing new methods for determining biomarkers and therapeutic targets.

ACSNI: An unsupervised machine-learning tool for prediction of tissue-specific pathway components using gene expression profiles

Determining the tissue- and disease-specific circuit of biological pathways remains a fundamental goal of molecular biology. Many components of these biological pathways still remain unknown, hindering the full and accurate characterization of biological processes of interest. Here we describe ACSNI, an algorithm that combines prior knowledge of biological processes with a deep neural network to effectively decompose gene expression profiles (GEPs) into multi-variable pathway activities and identify unknown pathway components. Experiments on public GEP data show that ACSNI predicts cogent components of mTOR, ATF2, and HOTAIRM1 signaling that recapitulate regulatory information from genetic perturbation and transcription factor binding datasets. Our framework provides a fast and easy-to-use method to identify components of signaling pathways as a tool for molecular mechanism discovery and to prioritize genes for designing future targeted experiments (https://github.com/caanene1/ACSNI).

Mantis: flexible and consensus-driven genome annotation

Background

The rapid development of the (meta-)omics fields has produced an unprecedented amount of high-resolution and high-fidelity data. Through the use of these datasets we can infer the role of previously functionally unannotated proteins from single organisms and consortia. In this context, protein function annotation can be described as the identification of regions of interest (i.e., domains) in protein sequences and the assignment of biological functions. Despite the existence of numerous tools, challenges remain in terms of speed, flexibility, and reproducibility. In the big data era, it is also increasingly important to cease limiting our findings to a single reference, coalescing knowledge from different data sources, and thus overcoming some limitations in overly relying on computationally generated data from single sources.

Results

We implemented a protein annotation tool, Mantis, which uses database identifiers intersection and text mining to integrate knowledge from multiple reference data sources into a single consensus-driven output. Mantis is flexible, allowing for the customization of reference data and execution parameters, and is reproducible across different research goals and user environments. We implemented a depth-first search algorithm for domain-specific annotation, which significantly improved annotation performance compared to sequence-wide annotation. The parallelized implementation of Mantis results in short runtimes while also outputting high coverage and high-quality protein function annotations.

Conclusions

Mantis is a protein function annotation tool that produces high-quality consensus-driven protein annotations. It is easy to set up, customize, and use, scaling from single genomes to large metagenomes. Mantis is available under the MIT license at https://github.com/PedroMTQ/mantis.

Recommendations for Statistical Reporting in Cardiovascular Medicine: A Special Report From the American Heart Association

Statistical analyses are a crucial component of the biomedical research process and are necessary to draw inferences from biomedical research data. The application of sound statistical methodology is a prerequisite for publication in the American Heart Association (AHA) journal portfolio. The objective of this document is to summarize key aspects of statistical reporting that might be most relevant to the authors, reviewers, and readership of AHA journals. The AHA Scientific Publication Committee convened a task force to inventory existing statistical standards for publication in biomedical journals and to identify approaches suitable for the AHA journal portfolio. The experts on the task force were selected by the AHA Scientific Publication Committee, who identified 12 key topics that serve as the section headers for this document. For each topic, the members of the writing group identified relevant references and evaluated them as a resource to make the standards summarized herein. Each section was independently reviewed by an expert reviewer who was not part of the task force. Expert reviewers were also permitted to comment on other sections if they chose. Differences of opinion were adjudicated by consensus. The standards presented in this report are intended to serve as a guide for high-quality reporting of statistical analyses methods and results.

Pathway-based analysis of anthocyanin diversity in diploid potato

Anthocyanin biosynthesis is one of the most studied pathways in plants due to the important ecological role played by these compounds and the potential health benefits of anthocyanin consumption. Given the interest in identifying new genetic factors underlying anthocyanin content we studied a diverse collection of diploid potatoes by combining a genome-wide association study and pathway-based analyses. By using an expanded SNP dataset, we identified candidate genes that had not been associated with anthocyanin variation in potatoes, namely a Myb transcription factor, a Leucoanthocyanidin dioxygenase gene and a vacuolar membrane protein. Importantly, a genomic region in chromosome 10 harbored the SNPs with strongest associations with anthocyanin content in GWAS. Some of these SNPs were associated with multiple anthocyanin compounds and therefore could underline the existence of pleiotropic genes or anthocyanin biosynthetic clusters. We identified multiple anthocyanin homologs in this genomic region, including four transcription factors and five enzymes that could be governing anthocyanin variation. For instance, a SNP linked to the phenylalanine ammonia-lyase gene, encoding the first enzyme in the phenylpropanoid biosynthetic pathway, was associated with all of the five anthocyanins measured. Finally, we combined a pathway analysis and GWAS of other agronomic traits to identify pathways related to anthocyanin biosynthesis in potatoes. We found that methionine metabolism and the production of sugars and hydroxycinnamic acids are genetically correlated to anthocyanin biosynthesis. The results contribute to the understanding of anthocyanins regulation in potatoes and can be used in future breeding programs focused on nutraceutical food.

A powerful subset-based method identifies gene set associations and improves interpretation in UK Biobank

Tests of association between a phenotype and a set of genes in a biological pathway can provide insights into the genetic architecture of complex phenotypes beyond those obtained from single-variant or single-gene association analysis. However, most existing gene set tests have limited power to detect gene set-phenotype association when a small fraction of the genes are associated with the phenotype and cannot identify the potentially "active" genes that might drive a gene set-based association. To address these issues, we have developed Gene set analysis Association Using Sparse Signals (GAUSS), a method for gene set association analysis that requires only GWAS summary statistics. For each significantly associated gene set, GAUSS identifies the subset of genes that have the maximal evidence of association and can best account for the gene set association. Using pre-computed correlation structure among test statistics from a reference panel, our p value calculation is substantially faster than other permutation- or simulation-based approaches. In simulations with varying proportions of causal genes, we find that GAUSS effectively controls type 1 error rate and has greater power than several existing methods, particularly when a small proportion of genes account for the gene set signal. Using GAUSS, we analyzed UK Biobank GWAS summary statistics for 10,679 gene sets and 1,403 binary phenotypes. We found that GAUSS is scalable and identified 13,466 phenotype and gene set association pairs. Within these gene sets, we identify an average of 17.2 (max = 405) genes that underlie these gene set associations.

CTD: An information-theoretic algorithm to interpret sets of metabolomic and transcriptomic perturbations in the context of graphical models

We consider the following general family of algorithmic problems that arises in transcriptomics, metabolomics and other fields: given a weighted graph G and a subset of its nodes S, find subsets of S that show significant connectedness within G. A specific solution to this problem may be defined by devising a scoring function, the Maximum Clique problem being a classic example, where S includes all nodes in G and where the score is defined by the size of the largest subset of S fully connected within G. Major practical obstacles for the plethora of algorithms addressing this type of problem include computational efficiency and, particularly for more complex scores which take edge weights into account, the computational cost of permutation testing, a statistical procedure required to obtain a bound on the p-value for a connectedness score. To address these problems, we developed CTD, "Connect the Dots", a fast algorithm based on data compression that detects highly connected subsets within S. CTD provides information-theoretic upper bounds on p-values when S contains a small fraction of nodes in G without requiring computationally costly permutation testing. We apply the CTD algorithm to interpret multi-metabolite perturbations due to inborn errors of metabolism and multi-transcript perturbations associated with breast cancer in the context of disease-specific Gaussian Markov Random Field networks learned directly from respective molecular profiling data.

An integrative multi-omics network-based approach identifies key regulators for breast cancer

Although genome-wide association studies (GWASs) have successfully identified thousands of risk variants for human complex diseases, understanding the biological function and molecular mechanisms of the associated SNPs involved in complex diseases is challenging. Here we developed a framework named integrative multi-omics network-based approach (IMNA), aiming to identify potential key genes in regulatory networks by integrating molecular interactions across multiple biological scales, including GWAS signals, gene expression-based signatures, chromatin interactions and protein interactions from the network topology. We applied this approach to breast cancer, and prioritized key genes involved in regulatory networks. We also developed an abnormal gene expression score (AGES) signature based on the gene expression deviation of the top 20 rank-ordered genes in breast cancer. The AGES values are associated with genetic variants, tumor properties and patient survival outcomes. Among the top 20 genes, RNASEH2A was identified as a new candidate gene for breast cancer. Thus, our integrative network-based approach provides a genetic-driven framework to unveil tissue-specific interactions from multiple biological scales and reveal potential key regulatory genes for breast cancer. This approach can also be applied in other complex diseases such as ovarian cancer to unravel underlying mechanisms and help for developing therapeutic targets.

Graph-based exploitation of gene ontology using GOxploreR for scrutinizing biological significance

Gene ontology (GO) is an eminent knowledge base frequently used for providing biological interpretations for the analysis of genes or gene sets from biological, medical and clinical problems. Unfortunately, the interpretation of such results is challenging due to the large number of GO terms, their hierarchical and connected organization as directed acyclic graphs (DAGs) and the lack of tools allowing to exploit this structural information explicitly. For this reason, we developed the R package GOxploreR. The main features of GOxploreR are (I) easy and direct access to structural features of GO, (II) structure-based ranking of GO-terms, (III) mapping to reduced GO-DAGs including visualization capabilities and (IV) prioritizing of GO-terms. The underlying idea of GOxploreR is to exploit a graph-theoretical perspective of GO as manifested by its DAG-structure and the containing hierarchy levels for cumulating semantic information. That means all these features enhance the utilization of structural information of GO and complement existing analysis tools. Overall, GOxploreR provides exploratory as well as confirmatory tools for complementing any kind of analysis resulting in a list of GO-terms, e.g., from differentially expressed genes or gene sets, GWAS or biomarkers. Our R package GOxploreR is freely available from CRAN.

Genomic Evidence for Sensorial Adaptations to a Nocturnal Predatory Lifestyle in Owls

Owls (Strigiformes) evolved specific adaptations to their nocturnal predatory lifestyle, such as asymmetrical ears, a facial disk, and a feather structure allowing silent flight. Owls also share some traits with diurnal raptors and other nocturnal birds, such as cryptic plumage patterns, reversed sexual size dimorphism, and acute vision and hearing. The genetic basis of some of these adaptations to a nocturnal predatory lifestyle has been studied by candidate gene approaches but rarely with genome-wide scans. Here, we used a genome-wide comparative analysis to test for selection in the early history of the owls. We estimated the substitution rates in the coding regions of 20 bird genomes, including 11 owls of which five were newly sequenced. Then, we tested for functional overrepresentation across the genes that showed signals of selection. In the ancestral branch of the owls, we found traces of positive selection in the evolution of genes functionally related to visual perception, especially to phototransduction, and to chromosome packaging. Several genes that have been previously linked to acoustic perception, circadian rhythm, and feather structure also showed signals of an accelerated evolution in the origin of the owls. We discuss the functions of the genes under positive selection and their putative association with the adaptation to the nocturnal predatory lifestyle of the owls.

Integrative Analysis of MAPK14 as a Potential Biomarker for Cardioembolic Stroke

The aim of this study was to obtain the candidate genes and biomarkers that are significantly related to cardioembolic stroke (CS) by applying bioinformatics analysis. In accordance with the results of the weighted gene coexpression network analysis (WGCNA) in the GSE58294 dataset, 11 CS-related coexpression network modules were identified in this study. Correlation analysis showed that the black and pink modules are significantly associated with CS. A total of 18 core genes in the black module and one core gene in the pink module were determined. We then identified differentially expressed genes (DEGs) of CS at 3 h, 5 h, and 24 h postonset. After performing intersection, it was found that 311 genes were coexpressed at these three time points. These genes were majorly enriched in positive regulation of transferase activity and regulation of peptidase activity. The abovementioned coexpressed DEGs were subjected to protein-protein interaction analysis and subnetwork module analysis. Subsequently, we used cytoHubba to obtain 11 key genes from DEGs. The intersection of the core genes screened from WGCNA and the key genes selected from DEGs yielded the MAPK14 gene. The expression level of MAPK14 on the receiver operating characteristic (ROC) curves of CS at 3 h, 5 h, and 24 h showed that the area under the ROC curve (AUC) was 0.923, 0.934, and 0.941, respectively. In a nutshell, MAPK14 screened out by using WGCNA showed differential expression in CS. We conclude that MAPK14 can be used as a potential biological marker of CS and exhibits potential to predict the physiopathological condition of CS patients.

Hyperketonemia GWAS and parity-dependent SNP associations in Holstein dairy cows intensively sampled for blood β-hydroxybutyrate concentration

Hyperketonemia (HYK) is a metabolic disorder that affects early postpartum dairy cows; however, there has been limited success in identifying genomic variants contributing to HYK susceptibility. We conducted a genome-wide association study (GWAS) using HYK phenotypes based on an intensive screening protocol, interrogated genotype interactions with parity group (GWIS), and evaluated the enrichment of annotated metabolic pathways. Holstein cows were enrolled into the experiment after parturition, and blood samples were collected at four timepoints between 5 and 18 days postpartum. Concentration of blood β-hydroxybutyrate (BHB) was quantified cow-side via a handheld BHB meter. Cows were labeled as a HYK case when at least one blood sample had BHB ≥ 1.2 mmol/L, and all other cows were considered non-HYK controls. After quality control procedures, 1,710 cows and 58,699 genotypes were available for further analysis. The GWAS and GWIS were performed using the forward feature select linear mixed model method. There was evidence for an association between ARS-BFGL-NGS-91238 and HYK susceptibility, as well as parity-dependent associations to HYK for BovineHD0600024247 and BovineHD1400023753. Candidate genes annotated to these single nuclear polymorphism associations have been previously associated with obesity, diabetes, insulin resistance, and fatty liver in humans and rodent models. Enrichment analysis revealed focal adhesion and axon guidance as metabolic pathways contributing to HYK etiology, while genetic variation in pathways related to insulin secretion and sensitivity may affect HYK susceptibility in a parity-dependent matter. In conclusion, the present work proposes several novel marker associations and metabolic pathways contributing to genetic risk for HYK susceptibility.

The Complex Association of FcγRIIb With Autoimmune Susceptibility

FcγRIIb is the only inhibitory Fc receptor and controls many aspects of immune and inflammatory responses. The observation 19 years ago that Fc γ RIIb ^-/- mice generated by gene targeting in 129 derived ES cells developed severe lupus like disease when backcrossed more than 7 generations into C57BL/6 background initiated extensive research on the functional understanding of this strong autoimmune phenotype. The genomic region in the distal part of Chr1 both in human and mice in which the Fc γ R gene cluster is located shows a high level of complexity in relation to the susceptibility to SLE. Specific haplotypes of closely linked genes including the Fc γ RIIb and Slamf genes are associated with increased susceptibility to SLE both in mice and human. Using forward and reverse genetic approaches including in human GWAS and in mice congenic strains, KO mice (germline and cell type specific, on different genetic background), knockin mice, overexpressing transgenic mice combined with immunological models such as adoptive transfer of B cells from Ig transgenic mice the involved genes and the causal mutations and their associated functional alterations were analyzed. In this review the results of this 19 years extensive research are discussed with a focus on (genetically modified) mouse models.

Gene set enrichment analysis to create polygenic scores: a developmental examination of aggression

Previous approaches for creating polygenic risk scores (PRSs) do not explicitly consider the biological or developmental relevance of the genetic variants selected for inclusion. We applied gene set enrichment analysis to meta-GWAS data to create developmentally targeted, functionally informed PRSs. Using two developmentally matched meta-GWAS discovery samples, separate PRSs were formed, then examined in time-varying effect models of aggression in a second, longitudinal sample of children (n = 515, 49% female) in early childhood (2-5 years old), and middle childhood (7.5-10.5 years old). Functional PRSs were associated with aggression in both the early and middle childhood models.

Genomic Regions Associated with IgE Levels against Culicoides spp. Antigens in Three Horse Breeds

Insect bite hypersensitivity (IBH), which is a cutaneous allergic reaction to antigens from Culicoides spp., is the most prevalent skin disorder in horses. Misdiagnosis is possible, as IBH is usually diagnosed based on clinical signs. Our study is the first to employ IgE levels against several recombinant Culicoides spp. allergens as an objective, independent, and quantitative phenotype to improve the power to detect genetic variants that underlie IBH. Genotypes of 200 Shetland ponies, 127 Icelandic horses, and 223 Belgian Warmblood horses were analyzed while using a mixed model approach. No single-nucleotide polymorphism (SNP) passed the Bonferroni corrected significance threshold, but several regions were identified within and across breeds, which confirmed previously identified regions of interest and, in addition, identifying new regions of interest. Allergen-specific IgE levels are a continuous and objective phenotype that allow for more powerful analyses when compared to a case-control set-up, as more significant associations were obtained. However, the use of a higher density array seems necessary to fully employ the use of IgE levels as a phenotype. While these results still require validation in a large independent dataset, the use of allergen-specific IgE levels showed value as an objective and continuous phenotype that can deepen our understanding of the biology underlying IBH.

Simultaneous inference of phenotype-associated genes and relevant tissues from GWAS data via Bayesian integration of multiple tissue-specific gene networks

Although genome-wide association studies (GWAS) have successfully identified thousands of genomic loci associated with hundreds of complex traits in the past decade, the debate about such problems as missing heritability and weak interpretability has been appealing for effective computational methods to facilitate the advanced analysis of the vast volume of existing and anticipated genetic data. Towards this goal, gene-level integrative GWAS analysis with the assumption that genes associated with a phenotype tend to be enriched in biological gene sets or gene networks has recently attracted much attention, due to such advantages as straightforward interpretation, less multiple testing burdens, and robustness across studies. However, existing methods in this category usually exploit non-tissue-specific gene networks and thus lack the ability to utilize informative tissue-specific characteristics. To overcome this limitation, we proposed a Bayesian approach called SIGNET (Simultaneously Inference of GeNEs and Tissues) to integrate GWAS data and multiple tissue-specific gene networks for the simultaneous inference of phenotype-associated genes and relevant tissues. Through extensive simulation studies, we showed the effectiveness of our method in finding both associated genes and relevant tissues for a phenotype. In applications to real GWAS data of 14 complex phenotypes, we demonstrated the power of our method in both deciphering genetic basis and discovering biological insights of a phenotype. With this understanding, we expect to see SIGNET as a valuable tool for integrative GWAS analysis, thereby boosting the prevention, diagnosis, and treatment of human inherited diseases and eventually facilitating precision medicine.

Disease classification: from phenotypic similarity to integrative genomics and beyond

A fundamental challenge of modern biomedical research is understanding how diseases that are similar on the phenotypic level are similar on the molecular level. Integration of various genomic data sets with the traditionally used phenotypic disease similarity revealed novel genetic and molecular mechanisms and blurred the distinction between monogenic (Mendelian) and complex diseases. Network-based medicine has emerged as a complementary approach for identifying disease-causing genes, genetic mediators, disruptions in the underlying cellular functions and for drug repositioning. The recent development of machine and deep learning methods allow for leveraging real-life information about diseases to refine genetic and phenotypic disease relationships. This review describes the historical development and recent methodological advancements for studying disease classification (nosology).

Muscle Gene Sets: a versatile methodological aid to functional genomics in the neuromuscular field

BACKGROUND

The approach of building large collections of gene sets and then systematically testing hypotheses across these collections is a powerful tool in functional genomics, both in the pathway analysis of omics data and to uncover the polygenic effects associated with complex diseases in genome-wide association study. The Molecular Signatures Database includes collections of oncogenic and immunologic signatures enabling researchers to compare transcriptional datasets across hundreds of previous studies and leading to important insights in these fields, but such a resource does not currently exist for neuromuscular research. In previous work, we have shown the utility of gene set approaches to understand muscle cell physiology and pathology.

METHODS

Following a systematic survey of public muscle data, we passed gene expression profiles from 4305 samples through a robust pre-processing and standardized data analysis pipeline. Two hundred eighty-two samples were discarded based on a battery of rigorous global quality controls. From among the remaining studies, 578 comparisons of interest were identified by a combination of text mining and manual curation of the study meta-data. For each comparison, significantly dysregulated genes (FDR adjusted p < 0.05) were identified.

RESULTS

Lists of dysregulated genes were divided between upregulated and downregulated to give 1156 Muscle Gene Sets (MGS). This resource is available for download ( www.sys-myo.com/muscle_gene_sets ) and is accessible through three commonly used functional genomics platforms (GSEA, EnrichR, and WebGestalt). Basic guidance and recommendations are provided for the use of MGS through these platforms. In addition, consensus muscle gene sets were created to capture the overlap between the results of similar studies, and analysis of these highlighted the potential for novel disease-relevant findings.

CONCLUSIONS

The MGS resource can be used to investigate the behaviour of any list of genes across previous comparisons of muscle conditions, to compare previous studies to one another, and to explore the functional relationship of muscle dysregulation to the Gene Ontology. Its major intended use is in enrichment testing for functional genomics analysis.

Pathway analysis of GWAS loci identifies novel drug targets and repurposing opportunities

Genome-wide association studies (GWAS) have made considerable progress and there is emerging evidence that genetics-based targets can lead to 28% more launched drugs. We analyzed 1589 GWAS across 1456 pathways to translate these often imprecise genetic loci into therapeutic hypotheses for 182 diseases. These pathway-based genetic targets were validated by testing whether current drug targets were enriched in the pathway space for the same indication. Remarkably, 30% of diseases had significantly more targets in these pathways than expected by chance; the comparable number for GWAS alone (without pathway analysis) was zero. This study shows that a systematic global pathway analysis can translate genetic findings into therapeutic hypotheses for both new drug discovery and repositioning opportunities for current drugs.

Genetic interactions between Ror2 and Wnt9a, Ror1 and Wnt9a and Ror2 and Ror1: Phenotypic analysis of the limb skeleton and palate in compound mutants

Mutations in the human receptor tyrosine kinase ROR2 are associated with Robinow syndrome (RRS) and brachydactyly type B1. Amongst others, the shortened limb phenotype associated with RRS is recapitulated in Ror2^-/- mutant mice. In contrast, Ror1^-/- mutant mice are viable and show no limb phenotype. Ror1^-/- ;Ror2^-/- double mutants are embryonic lethal, whereas double mutants containing a hypomorphic Ror1 allele (Ror1^hyp ) survive up to birth and display a more severe shortened limb phenotype. Both orphan receptors have been shown to act as possible Wnt coreceptors and to mediate the Wnt5a signal. Here, we analyzed genetic interactions between the Wnt ligand, Wnt9a, and Ror2 or Ror1, as Wnt9a has also been implicated in skeletal development. Wnt9a^-/- single mutants display a mild shortening of the long bones, whereas these are severely shortened in Ror2^-/- mutants. Ror2^-/- ;Wnt9a^-/- double mutants displayed even more severely shortened long bones, and intermediate phenotypes were observed in compound Ror2;Wnt9a mutants. Long bones were also shorter in Ror1^hyp/hyp ;Wnt9a^-/- double mutants. In addition, Ror1^hyp/hyp ;Wnt9a^-/- double mutants displayed a secondary palate cleft phenotype, which was not present in the respective single mutants. Interestingly, 50% of compound mutant pups heterozygous for Ror2 and homozygous mutant for Ror1 also developed a secondary palate cleft phenotype.

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

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

The Genomic Basis of Tumor Regression in Tasmanian Devils (Sarcophilus harrisii)

Understanding the genetic basis of disease-related phenotypes, such as cancer susceptibility, is crucial for the advancement of personalized medicine. Although most cancers are somatic in origin, a small number of transmissible cancers have been documented. Two such cancers have emerged in the Tasmanian devil (Sarcophilus harrisii) and now threaten the species with extinction. Recently, cases of natural tumor regression in Tasmanian devils infected with the clonally contagious cancer have been detected. We used whole-genome sequencing and FST-based approaches to identify the genetic basis of tumor regression by comparing the genomes of seven individuals that underwent tumor regression with those of three infected individuals that did not. We found three highly differentiated candidate genomic regions containing several genes related to immune response and/or cancer risk, indicating that the genomic basis of tumor regression was polygenic. Within these genomic regions, we identified putative regulatory variation in candidate genes but no nonsynonymous variation, suggesting that natural tumor regression may be driven, at least in part, by differential host expression of key loci. Comparative oncology can provide insight into the genetic basis of cancer risk, tumor development, and the pathogenicity of cancer, particularly due to our limited ability to monitor natural, untreated tumor progression in human patients. Our results support the hypothesis that host immune response is necessary for triggering tumor regression, providing candidate genes that may translate to novel treatments in human and nonhuman cancers.

Genome-Wide Association Study in African Americans with Acute Respiratory Distress Syndrome Identifies the Selectin P Ligand Gene as a Risk Factor

RATIONALE

Genetic factors are involved in acute respiratory distress syndrome (ARDS) susceptibility. Identification of novel candidate genes associated with increased risk and severity will improve our understanding of ARDS pathophysiology and enhance efforts to develop novel preventive and therapeutic approaches.

OBJECTIVES

To identify genetic susceptibility targets for ARDS.

METHODS

A genome-wide association study was performed on 232 African American patients with ARDS and 162 at-risk control subjects. The Identify Candidate Causal SNPs and Pathways platform was used to infer the association of known gene sets with the top prioritized intragenic SNPs. Preclinical validation of SELPLG (selectin P ligand gene) was performed using mouse models of LPS- and ventilator-induced lung injury. Exonic variation within SELPLG distinguishing patients with ARDS from sepsis control subjects was confirmed in an independent cohort.

MEASUREMENTS AND MAIN RESULTS

Pathway prioritization analysis identified a nonsynonymous coding SNP (rs2228315) within SELPLG, encoding P-selectin glycoprotein ligand 1, to be associated with increased susceptibility. In an independent cohort, two exonic SELPLG SNPs were significantly associated with ARDS susceptibility. Additional support for SELPLG as an ARDS candidate gene was derived from preclinical ARDS models where SELPLG gene expression in lung tissues was significantly increased in both ventilator-induced (twofold increase) and LPS-induced (5.7-fold increase) murine lung injury models compared with controls. Furthermore, Selplg-/- mice exhibited significantly reduced LPS-induced inflammatory lung injury compared with wild-type C57/B6 mice. Finally, an antibody that neutralizes P-selectin glycoprotein ligand 1 significantly attenuated LPS-induced lung inflammation.

CONCLUSIONS

These findings identify SELPLG as a novel ARDS susceptibility gene among individuals of European and African descent.

The Complex Interaction of Mitochondrial Genetics and Mitochondrial Pathways in Psychiatric Disease

While accounting for only 2% of the body's weight, the brain utilizes up to 20% of the body's total energy. Not surprisingly, metabolic dysfunction and energy supply-and-demand mismatch have been implicated in a variety of neurological and psychiatric disorders. Mitochondria are responsible for providing the brain with most of its energetic demands, and the brain uses glucose as its exclusive energy source. Exploring the role of mitochondrial dysfunction in the etiology of psychiatric disease is a promising avenue to investigate further. Genetic analysis of mitochondrial activity is a cornerstone in understanding disease pathogenesis related to metabolic dysfunction. In concert with neuroimaging and pathological study, genetics provides an important bridge between biochemical findings and clinical correlates in psychiatric disease. Mitochondrial genetics has several unique aspects to its analysis, and corresponding special considerations. Here, we review the components of mitochondrial genetic analysis - nuclear DNA, mitochon-drial DNA, mitochondrial pathways, pseudogenes, nuclear-mitochondrial mismatch, and microRNAs - that could contribute to an observable clinical phenotype. Throughout, we highlight psychiatric diseases that can arise due to dysfunction in these processes, with a focus on schizophrenia and bipolar disorder.

Functional Validation of Candidate Genes Detected by Genomic Feature Models

Understanding the genetic underpinnings of complex traits requires knowledge of the genetic variants that contribute to phenotypic variability. Reliable statistical approaches are needed to obtain such knowledge. In genome-wide association studies, variants are tested for association with trait variability to pinpoint loci that contribute to the quantitative trait. Because stringent genome-wide significance thresholds are applied to control the false positive rate, many true causal variants can remain undetected. To ameliorate this problem, many alternative approaches have been developed, such as genomic feature models (GFM). The GFM approach tests for association of set of genomic markers, and predicts genomic values from genomic data utilizing prior biological knowledge. We investigated to what degree the findings from GFM have biological relevance. We used the Drosophila Genetic Reference Panel to investigate locomotor activity, and applied genomic feature prediction models to identify gene ontology (GO) categories predictive of this phenotype. Next, we applied the covariance association test to partition the genomic variance of the predictive GO terms to the genes within these terms. We then functionally assessed whether the identified candidate genes affected locomotor activity by reducing gene expression using RNA interference. In five of the seven candidate genes tested, reduced gene expression altered the phenotype. The ranking of genes within the predictive GO term was highly correlated with the magnitude of the phenotypic consequence of gene knockdown. This study provides evidence for five new candidate genes for locomotor activity, and provides support for the reliability of the GFM approach.

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.

Calcium channel blockers as drug repurposing candidates for gestational diabetes: Mining large scale genomic and electronic health records data to repurpose medications

New therapeutic approaches are needed for gestational diabetes mellitus (GDM), but must show safety and efficacy in a historically understudied population. We studied associations between electronic medical record (EMR) phenotypes and genetic variants to uncover drugs currently considered safe in pregnancy that could treat or prevent GDM. We identified 129 systemically active drugs considered safe in pregnancy targeting the proteins produced from 196 genes. We tested for associations between GDM and/or type 2 diabetes (DM2) and 306 SNPs in 130 genes represented on the Illumina Infinium Human Exome Bead Chip (DM2 was included due to shared pathophysiological features with GDM). In parallel, we tested the association between drugs and glucose tolerance during pregnancy as measured by the glucose recorded during a routine 50-g glucose tolerance test (GTT). We found an association between GDM/DM2 and the genes targeted by 11 drug classes. In the EMR analysis, 6 drug classes were associated with changes in GTT. Two classes were identified in both analyses. L-type calcium channel blocking antihypertensives (CCBs), were associated with a 3.18 mg/dL (95% CI -6.18 to -0.18) decrease in glucose during GTT, and serotonin receptor type 3 (5HT-3) antagonist antinausea medications were associated with a 3.54 mg/dL (95% CI 1.86-5.23) increase in glucose during GTT. CCBs were identified as a class of drugs considered safe in pregnancy could have efficacy in treating or preventing GDM. 5HT-3 antagonists may be associated with worse glucose tolerance.

Methods and Tools in Genome-wide Association Studies

Many traits, such as height, the response to a given drug, or the susceptibility to certain diseases are presumably co-determined by genetics. Especially in the field of medicine, it is of major interest to identify genetic aberrations that alter an individual's risk to develop a certain phenotypic trait. Addressing this question requires the availability of comprehensive, high-quality genetic datasets. The technological advancements and the decreasing cost of genotyping in the last decade led to an increase in such datasets. Parallel to and in line with this technological progress, an analysis framework under the name of genome-wide association studies was developed to properly collect and analyze these data. Genome-wide association studies aim at finding statistical dependencies-or associations-between a trait of interest and point-mutations in the DNA. The statistical models used to detect such associations are diverse, spanning the whole range from the frequentist to the Bayesian setting.Since genetic datasets are inherently high-dimensional, the search for associations poses not only a statistical but also a computational challenge. As a result, a variety of toolboxes and software packages have been developed, each implementing different statistical methods while using various optimizations and mathematical techniques to enhance the computations.This chapter is devoted to the discussion of widely used methods and tools in genome-wide association studies. We present the different statistical models and the assumptions on which they are based, explain peculiarities of the data that have to be accounted for and, most importantly, introduce commonly used tools and software packages for the different tasks in a genome-wide association study, complemented with examples for their application.

Nonmonotonic Pathway Gene Expression Analysis Reveals Oncogenic Role of p27/Kip1 at Intermediate Dose

The mechanistic basis by which the level of p27^Kip1 expression influences tumor aggressiveness and patient mortality remains unclear. To elucidate the competing tumor-suppressing and oncogenic effects of p27^Kip1 on gene expression in tumors, we analyzed the transcriptomes of squamous cell papilloma derived from Cdkn1b nullizygous, heterozygous, and wild-type mice. We developed a novel functional pathway analysis method capable of testing directional and nonmonotonic dose response. This analysis can reveal potential causal relationships that might have been missed by other nondirectional pathway analysis methods. Applying this method to capture dose-response curves in papilloma gene expression data, we show that several known cancer pathways are dominated by low-high-low gene expression responses to increasing p27 gene doses. The oncogene cyclin D1, whose expression is elevated at an intermediate p27 dose, is the most responsive gene shared by these cancer pathways. Therefore, intermediate levels of p27 may promote cellular processes favoring tumorigenesis-strikingly consistent with the dominance of heterozygous mutations in CDKN1B seen in human cancers. Our findings shed new light on regulatory mechanisms for both pro- and anti-tumorigenic roles of p27^Kip1. Functional pathway dose-response analysis provides a unique opportunity to uncover nonmonotonic patterns in biological systems.

HPA Axis Genes, and Their Interaction with Childhood Maltreatment, are Related to Cortisol Levels and Stress-Related Phenotypes

Stress responses are controlled by the hypothalamus pituitary adrenal (HPA)-axis and maladaptive stress responses are associated with the onset and maintenance of stress-related disorders such as major depressive disorder (MDD). Genes that play a role in the HPA-axis regulation may likely contribute to the relation between relevant neurobiological substrates and stress-related disorders. Therefore, we performed gene-wide analyses for 30 a priori literature-based genes involved in HPA-axis regulation in 2014 subjects (34% male; mean age: 42.5) to study the relations with lifetime MDD diagnosis, cortisol awakening response, and dexamethasone suppression test (DST) levels (subsample N=1472) and hippocampal and amygdala volume (3T MR images; subsample N=225). Additionally, gene by childhood maltreatment (CM) interactions were investigated. Gene-wide significant results were found for dexamethasone suppression (CYP11A1, CYP17A1, POU1F1, AKR1D1), hippocampal volume (CYP17A1, CYP11A1, HSD3B2, PROP1, AVPRA1, SRD5A1), amygdala volume (POMC, CRH, HSD3B2), and lifetime MDD diagnosis (FKBP5 and CRH), all permutation p-values<0.05. Interactions with CM were found for several genes; the strongest interactions were found for NR3C2, where the minor allele of SNP rs17581262 was related to smaller hippocampal volume, smaller amygdala volume, higher DST levels, and higher odds of MDD diagnosis only in participants with CM. As hypothesized, several HPA-axis genes are associated with stress-related endophenotypes including cortisol response and reduced brain volumes. Furthermore, we found a pleiotropic interaction between CM and the mineralocorticoid receptor gene, suggesting that this gene plays an important moderating role in stress and stress-related disorders.

Genome-wide imaging association study implicates functional activity and glial homeostasis of the caudate in smoking addiction

Background

Nearly 6 million deaths and over a half trillion dollars in healthcare costs worldwide are attributed to tobacco smoking each year. Extensive research efforts have been pursued to elucidate the molecular underpinnings of smoking addiction and facilitate cessation. In this study, we genotyped and obtained both resting state and task-based functional magnetic resonance imaging from 64 non-smokers and 42 smokers. Smokers were imaged after having smoked normally (“sated”) and after having not smoked for at least 12 h (“abstinent”).

Results

While abstinent smokers did not differ from non-smokers with respect to pairwise resting state functional connectivities (RSFCs) between 12 brain regions of interest, RSFCs involving the caudate and putamen of sated smokers significantly differed from those of non-smokers (P < 0.01). Further analyses of caudate and putamen activity during elicited experiences of reward and disappointment show that caudate activity during reward (CR) correlated with smoking status (P = 0.015). Moreover, abstinent smokers with lower CR experienced greater withdrawal symptoms (P = 0.024), which suggests CR may be related to smoking urges. Associations between genetic variants and CR, adjusted for smoking status, were identified by genome-wide association study (GWAS). Genes containing or exhibiting caudate-specific expression regulation by these variants were enriched within Gene Ontology terms that describe cytoskeleton functions, synaptic organization, and injury response (P < 0.001, FDR < 0.05).

Conclusions

By integrating genomic and imaging data, novel insights into potential mechanisms of caudate activation and homeostasis are revealed that may guide new directions of research toward improving our understanding of addiction pathology.

Electronic supplementary material

The online version of this article (10.1186/s12864-017-4124-5) contains supplementary material, which is available to authorized users.

Shared genetic regulatory networks for cardiovascular disease and type 2 diabetes in multiple populations of diverse ethnicities in the United States

Cardiovascular diseases (CVD) and type 2 diabetes (T2D) are closely interrelated complex diseases likely sharing overlapping pathogenesis driven by aberrant activities in gene networks. However, the molecular circuitries underlying the pathogenic commonalities remain poorly understood. We sought to identify the shared gene networks and their key intervening drivers for both CVD and T2D by conducting a comprehensive integrative analysis driven by five multi-ethnic genome-wide association studies (GWAS) for CVD and T2D, expression quantitative trait loci (eQTLs), ENCODE, and tissue-specific gene network models (both co-expression and graphical models) from CVD and T2D relevant tissues. We identified pathways regulating the metabolism of lipids, glucose, and branched-chain amino acids, along with those governing oxidation, extracellular matrix, immune response, and neuronal system as shared pathogenic processes for both diseases. Further, we uncovered 15 key drivers including HMGCR, CAV1, IGF1 and PCOLCE, whose network neighbors collectively account for approximately 35% of known GWAS hits for CVD and 22% for T2D. Finally, we cross-validated the regulatory role of the top key drivers using in vitro siRNA knockdown, in vivo gene knockout, and two Hybrid Mouse Diversity Panels each comprised of >100 strains. Findings from this in-depth assessment of genetic and functional data from multiple human cohorts provide strong support that common sets of tissue-specific molecular networks drive the pathogenesis of both CVD and T2D across ethnicities and help prioritize new therapeutic avenues for both CVD and T2D.

Genomic Analysis of Genotype-by-Social Environment Interaction for Drosophila melanogaster Aggressive Behavior

Human psychiatric disorders such as schizophrenia, bipolar disorder, and attention-deficit/hyperactivity disorder often include adverse behaviors including increased aggressiveness. Individuals with psychiatric disorders often exhibit social withdrawal, which can further increase the probability of conducting a violent act. Here, we used the inbred, sequenced lines of the Drosophila Genetic Reference Panel (DGRP) to investigate the genetic basis of variation in male aggressive behavior for flies reared in a socialized and socially isolated environment. We identified genetic variation for aggressive behavior, as well as significant genotype-by-social environmental interaction (GSEI); i.e., variation among DGRP genotypes in the degree to which social isolation affected aggression. We performed genome-wide association (GWA) analyses to identify genetic variants associated with aggression within each environment. We used genomic prediction to partition genetic variants into gene ontology (GO) terms and constituent genes, and identified GO terms and genes with high prediction accuracies in both social environments and for GSEI. The top predictive GO terms significantly increased the proportion of variance explained, compared to prediction models based on all segregating variants. We performed genomic prediction across environments, and identified genes in common between the social environments that turned out to be enriched for genome-wide associated variants. A large proportion of the associated genes have previously been associated with aggressive behavior in Drosophila and mice. Further, many of these genes have human orthologs that have been associated with neurological disorders, indicating partially shared genetic mechanisms underlying aggression in animal models and human psychiatric disorders.

A Pathway-Centered Analysis of Pig Domestication and Breeding in Eurasia

Ascertaining the molecular and physiological basis of domestication and breeding is an active area of research. Due to the current wide distribution of its wild ancestor, the wild boar, the pig (Sus scrofa) is an excellent model to study these processes, which occurred independently in East Asia and Europe ca. 9000 yr ago. Analyzing genome variability patterns in terms of metabolic pathways is attractive since it considers the impact of interrelated functions of genes, in contrast to genome-wide scans that treat genes or genome windows in isolation. To that end, we studied 40 wild boars and 123 domestic pig genomes from Asia and Europe when metabolic pathway was the unit of analysis. We computed statistical significance for differentiation (Fst) and linkage disequilibrium (nSL) statistics at the pathway level. In terms of Fst, we found 21 and 12 pathways significantly differentiated at a q-value < 0.05 in Asia and Europe, respectively; five were shared across continents. In Asia, we found six significant pathways related to behavior, which involved essential neurotransmitters like dopamine and serotonin. Several significant pathways were interrelated and shared a variable percentage of genes. There were 12 genes present in >10 significant pathways (in terms of Fst), comprising genes involved in the transduction of a large number of signals, like phospholipase PCLB1, which is expressed in the brain, or ITPR3, which has an important role in taste transduction. In terms of nSL, significant pathways were mainly related to reproductive performance (ovarian steroidogenesis), a similarly important target trait during domestication and modern animal breeding. Different levels of recombination cannot explain these results, since we found no correlation between Fst and recombination rate. However, we did find an increased ratio of deleterious mutations in domestic vs. wild populations, suggesting a relaxed functional constraint associated with the domestication and breeding processes. Purifying selection was, nevertheless, stronger in significantly differentiated pathways than in random pathways, mainly in Europe. We conclude that pathway analysis facilitates the biological interpretation of genome-wide studies. Notably, in the case of pig, behavior played an important role, among other physiological and developmental processes.

Susceptibility to Childhood Pneumonia: A Genome-Wide Analysis

Previous studies have indicated that in adult smokers, a history of childhood pneumonia is associated with reduced lung function and chronic obstructive pulmonary disease. There have been few previous investigations using genome-wide association studies to investigate genetic predisposition to pneumonia. This study aims to identify the genetic variants associated with the development of pneumonia during childhood and over the course of the lifetime. Study subjects included current and former smokers with and without chronic obstructive pulmonary disease participating in the COPDGene Study. Pneumonia was defined by subject self-report, with childhood pneumonia categorized as having the first episode at <16 years. Genome-wide association studies for childhood pneumonia (843 cases, 9,091 control subjects) and lifetime pneumonia (3,766 cases, 5,659 control subjects) were performed separately in non-Hispanic whites and African Americans. Non-Hispanic white and African American populations were combined in the meta-analysis. Top genetic variants from childhood pneumonia were assessed in network analysis. No single-nucleotide polymorphisms reached genome-wide significance, although we identified potential regions of interest. In the childhood pneumonia analysis, this included variants in NGR1 (P = 6.3 × 10^-8), PAK6 (P = 3.3 × 10^-7), and near MATN1 (P = 2.8 × 10^-7). In the lifetime pneumonia analysis, this included variants in LOC339862 (P = 8.7 × 10^-7), RAPGEF2 (P = 8.4 × 10^-7), PHACTR1 (P = 6.1 × 10^-7), near PRR27 (P = 4.3 × 10^-7), and near MCPH1 (P = 2.7 × 10^-7). Network analysis of the genes associated with childhood pneumonia included top networks related to development, blood vessel morphogenesis, muscle contraction, WNT signaling, DNA damage, apoptosis, inflammation, and immune response (P ≤ 0.05). We have identified genes potentially associated with the risk of pneumonia. Further research will be required to confirm these associations and to determine biological mechanisms.

CLINICAL TRIAL REGISTRATION

NCT00608764.

Downstream targets of GWAS-detected genes for breast, lung, and prostate and colon cancer converge to G1/S transition pathway

Genome-wide association studies (GWASs) identified over 500 single nucleotide polymorphisms (SNPs) influencing cancer risk. It is logical to expect the cancer-associated genes to cluster in pathways directly involved in carcinogenesis, e.g. cell cycle. Nevertheless, analyses of the GWAS-detected cancer risk genes usually show no or weak enrichment by known cancer genes.We hypothesized that GWAS-detected cancer risk-associated genes function as upstream regulators of the genes directly involved in carcinogenesis. We have analyzed four common cancers: breast, colon, lung, and prostate. To identify downstream targets of GWAS-detected cancer risk genes we used MedScan, which is a text mining tool offered by PathwayStudio. We also used data on protein/protein interactions reported by BioGRID database. Among all identified targets we have selected common downstream targets. A gene was considered a common downstream target if it was a downstream target for at least three GWAS-detected genes for a given cancer type. Common downstream targets were identified separately for each cancer type. We found that common downstream targets for all four cancer types were enriched by cell cycle genes, more specifically, the genes involved in G1/S transition. Common downstream targets for bipolar disorder, Crohn's disease, and type 2 diabetes did not show G1/S transition enrichment.The results of this analysis suggest that many cancer risk genes function as upstream regulators of the genes directly involved in G1/S transition and exert their risk effects by reducing threshold for G1/S transition, elevating the background level of cell proliferation and cancer risk.

The statistical properties of gene-set analysis

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

Systems Genetics as a Tool to Identify Master Genetic Regulators in Complex Disease

Systems genetics stems from systems biology and similarly employs integrative modeling approaches to describe the perturbations and phenotypic effects observed in a complex system. However, in the case of systems genetics the main source of perturbation is naturally occurring genetic variation, which can be analyzed at the systems-level to explain the observed variation in phenotypic traits. In contrast with conventional single-variant association approaches, the success of systems genetics has been in the identification of gene networks and molecular pathways that underlie complex disease. In addition, systems genetics has proven useful in the discovery of master trans-acting genetic regulators of functional networks and pathways, which in many cases revealed unexpected gene targets for disease. Here we detail the central components of a fully integrated systems genetics approach to complex disease, starting from assessment of genetic and gene expression variation, linking DNA sequence variation to mRNA (expression QTL mapping), gene regulatory network analysis and mapping the genetic control of regulatory networks. By summarizing a few illustrative (and successful) examples, we highlight how different data-modeling strategies can be effectively integrated in a systems genetics study.