Network.assisted analysis to prioritize GWAS results: principles, methods and perspectives

Genomic Signatures of Domestication in European Seabass (Dicentrarchus labrax L.) Reveal a Potential Role for Epigenetic Regulation in Adaptation to Captivity

Genome scans provide a comprehensive method to explore genome-wide variation associated with traits under study. However, linking individual genes to broader functional groupings and pathways is often challenging, yet crucial for understanding the evolutionary mechanisms underlying these traits. This task is particularly relevant for multi-trait processes such as domestication, which are influenced by complex interactions between numerous genetic and non-genetic factors, including epigenetic regulation. As various traits within the broader spectrum of domestication are selected in concert over time, this process offers an opportunity to identify broader functional overlaps and understand the integrated genetic architecture underlying these traits. In this study, we analyzed approximately 600,000 SNPs from a Pool-Seq experiment comparing eight natural-origin and 12 farmed populations of European seabass in the Mediterranean Sea region. We implemented two genome scan approaches and focused on genomic regions supported by both methods, resulting in the identification of 96 candidate genes, including nine CpG islands, which highligt potential epigenetic influences. Many of these genes and CpG islands are in linkage groups previously associated with domestication-related traits. The most significantly overrepresented molecular function was "oxidoreductase activity". Furthermore, a dense network of interactions was identified, connecting 22 of the candidate genes. Within this network, the most significantly enriched pathways and central genes were involved in "chromatin organization", highlighting another potential epigenetic mechanism. Altogether, our findings underscore the utility of interactome-assisted pathway analysis in elucidating the genomic architecture of polygenic traits and suggest that epigenetic regulation may play a crucial role in the domestication of European seabass.

GWAS breakthroughs: mapping the journey from one locus to 393 significant coronary artery disease associations

Coronary artery disease (CAD) poses a substantial threat to global health, leading to significant morbidity and mortality worldwide. It has a significant genetic component that has been studied through genome-wide association studies (GWAS) over the past 17 years. These studies have made progress with larger sample sizes, diverse ancestral backgrounds, and the discovery of multiple genomic regions related to CAD risk. In this review, we provide a comprehensive overview of CAD GWAS, including information about the genetic makeup of the disease and the importance of ethnic diversity in these studies. We also discuss challenges of identifying causal genes and variants within GWAS loci with a focus on non-coding regions. Additionally, we highlight tissues and cell types relevant to CAD, and discuss clinical implications of GWAS findings including polygenic risk scores, sex-specific differences in CAD genetics, ethnical aspects of personalized interventions, and GWAS guided drug development.

Pharmacogenomics of clinical response to Natalizumab in multiple sclerosis: a genome-wide multi-centric association study

BACKGROUND

Inter-individual differences in treatment response are marked in multiple sclerosis (MS). This is true for Natalizumab (NTZ), to which a subset of patients displays sub-optimal treatment response. We conducted a multi-centric genome-wide association study (GWAS), with additional pathway and network analysis to identify genetic predictors of response to NTZ.

METHODS

MS patients from three different centers were included. Response to NTZ was dichotomized, nominating responders (R) relapse-free patients and non-responders (NR) all the others, over a follow-up of 4 years. Association analysis on ~ 4.7 M imputed autosomal common single-nucleotide polymorphisms (SNPs) was performed fitting logistic regression models, adjusted for baseline covariates, followed by meta-analysis at SNP and gene level. Finally, these signals were projected onto STRING interactome, to elicit modules and hub genes linked to response.

RESULTS

Overall, 1834 patients were included: 119 from Italy (R = 94, NR = 25), 81 from Germany (R = 61, NR = 20), and 1634 from Sweden (R = 1349, NR = 285). The top-associated variant was rs11132400T (p = 1.33 × 10-6, OR = 0.58), affecting expression of several genes in the locus, like KLKB1. The interactome analysis implicated a module of 135 genes, with over-representation of terms like canonical WNT signaling pathway (padjust = 7.08 × 10-6). Response-associated genes like GRB2 and LRP6, already implicated in MS pathogenesis, were topologically prioritized within the module.

CONCLUSION

This GWAS, the largest pharmacogenomic study of response to NTZ, suggested MS-implicated genes and Wnt/β-catenin signaling pathway, an essential component for blood-brain barrier formation and maintenance, to be related to treatment response.

Are trait-associated genes clustered together in a gene network?

Genome-wide association studies (GWAS) have provided an abundance of information about the genetic variants and their loci that are associated to complex traits and diseases. However, due to linkage disequilibrium (LD) and noncoding regions of loci, it remains a challenge to pinpoint the causal genes. Gene network-based approaches, paired with network diffusion methods, have been proposed to prioritize causal genes and to boost statistical power in GWAS based on the assumption that trait-associated genes are clustered in a gene network. Due to the difficulty in mapping trait-associated variants to genes in GWAS, this assumption has never been directly or rigorously tested empirically. On the other hand, whole exome sequencing (WES) data focuses on the protein-coding regions, directly identifying trait-associated genes. In this study, we tested the assumption by leveraging the recently available exome-based association statistics from the UK Biobank WES data along with two types of networks. We found that almost all trait-associated genes were significantly more proximal to each other than randomly selected genes within both networks. These results support the assumption that trait-associated genes are clustered in gene networks, which can be further leveraged to boost the power of GWAS such as by introducing less stringent p value thresholds.

Rare and Common Variants Associated with Alcohol Consumption Identify a Conserved Molecular Network

Genome-wide association studies (GWAS) have identified hundreds of common variants associated with alcohol consumption. In contrast, rare variants have only begun to be studied for their role in alcohol consumption. No studies have examined whether common and rare variants implicate the same genes and molecular networks. To address this knowledge gap, we used publicly available alcohol consumption GWAS summary statistics (GSCAN, N=666,978) and whole exome sequencing data (Genebass, N=393,099) to identify a set of common and rare variants for alcohol consumption. Gene-based analysis of each dataset have implicated 294 (common variants) and 35 (rare variants) genes, including ethanol metabolizing genes ADH1B and ADH1C, which were identified by both analyses, and ANKRD12, GIGYF1, KIF21B, and STK31, which were identified only by rare variant analysis, but have been associated with related psychiatric traits. We then used a network colocalization procedure to propagate the common and rare gene sets onto a shared molecular network, revealing significant overlap. The shared network identified gene families that function in alcohol metabolism, including ADH, ALDH, CYP, and UGT. 74 of the genes in the network were previously implicated in comorbid psychiatric or substance use disorders, but had not previously been identified for alcohol-related behaviors, including EXOC2, EPM2A, CACNB3, and CACNG4. Differential gene expression analysis showed enrichment in the liver and several brain regions supporting the role of network genes in alcohol consumption. Thus, genes implicated by common and rare variants identify shared functions relevant to alcohol consumption, which also underlie psychiatric traits and substance use disorders that are comorbid with alcohol use.

Protein-protein interaction network-based integration of GWAS and functional data for blood pressure regulation analysis

BACKGROUND

It is valuable to analyze the genome-wide association studies (GWAS) data for a complex disease phenotype in the context of the protein-protein interaction (PPI) network, as the related pathophysiology results from the function of interacting polyprotein pathways. The analysis may include the design and curation of a phenotype-specific GWAS meta-database incorporating genotypic and eQTL data linking to PPI and other biological datasets, and the development of systematic workflows for PPI network-based data integration toward protein and pathway prioritization. Here, we pursued this analysis for blood pressure (BP) regulation.

METHODS

The relational scheme of the implemented in Microsoft SQL Server BP-GWAS meta-database enabled the combined storage of: GWAS data and attributes mined from GWAS Catalog and the literature, Ensembl-defined SNP-transcript associations, and GTEx eQTL data. The BP-protein interactome was reconstructed from the PICKLE PPI meta-database, extending the GWAS-deduced network with the shortest paths connecting all GWAS-proteins into one component. The shortest-path intermediates were considered as BP-related. For protein prioritization, we combined a new integrated GWAS-based scoring scheme with two network-based criteria: one considering the protein role in the reconstructed by shortest-path (RbSP) interactome and one novel promoting the common neighbors of GWAS-prioritized proteins. Prioritized proteins were ranked by the number of satisfied criteria.

RESULTS

The meta-database includes 6687 variants linked with 1167 BP-associated protein-coding genes. The GWAS-deduced PPI network includes 1065 proteins, with 672 forming a connected component. The RbSP interactome contains 1443 additional, network-deduced proteins and indicated that essentially all BP-GWAS proteins are at most second neighbors. The prioritized BP-protein set was derived from the union of the most BP-significant by any of the GWAS-based or the network-based criteria. It included 335 proteins, with ~ 2/3 deduced from the BP PPI network extension and 126 prioritized by at least two criteria. ESR1 was the only protein satisfying all three criteria, followed in the top-10 by INSR, PTN11, CDK6, CSK, NOS3, SH2B3, ATP2B1, FES and FINC, satisfying two. Pathway analysis of the RbSP interactome revealed numerous bioprocesses, which are indeed functionally supported as BP-associated, extending our understanding about BP regulation.

CONCLUSIONS

The implemented workflow could be used for other multifactorial diseases.

scGraph2Vec: a deep generative model for gene embedding augmented by graph neural network and single-cell omics data

BACKGROUND

Exploring the cellular processes of genes from the aspects of biological networks is of great interest to understanding the properties of complex diseases and biological systems. Biological networks, such as protein-protein interaction networks and gene regulatory networks, provide insights into the molecular basis of cellular processes and often form functional clusters in different tissue and disease contexts.

RESULTS

We present scGraph2Vec, a deep learning framework for generating informative gene embeddings. scGraph2Vec extends the variational graph autoencoder framework and integrates single-cell datasets and gene-gene interaction networks. We demonstrate that the gene embeddings are biologically interpretable and enable the identification of gene clusters representing functional or tissue-specific cellular processes. By comparing similar tools, we showed that scGraph2Vec clearly distinguished different gene clusters and aggregated more biologically functional genes. scGraph2Vec can be widely applied in diverse biological contexts. We illustrated that the embeddings generated by scGraph2Vec can infer disease-associated genes from genome-wide association study data (e.g., COVID-19 and Alzheimer's disease), identify additional driver genes in lung adenocarcinoma, and reveal regulatory genes responsible for maintaining or transitioning melanoma cell states.

CONCLUSIONS

scGraph2Vec not only reconstructs tissue-specific gene networks but also obtains a latent representation of genes implying their biological functions.

Deciphering the genetic landscape of obesity: a data-driven approach to identifying plausible causal genes and therapeutic targets

OBJECTIVES

Genome-wide association studies (GWAS) have successfully revealed numerous susceptibility loci for obesity. However, identifying the causal genes, pathways, and tissues/cell types responsible for these associations remains a challenge, and standardized analysis workflows are lacking. Additionally, due to limited treatment options for obesity, there is a need for the development of new pharmacological therapies. This study aimed to address these issues by performing step-wise utilization of knowledgebase for gene prioritization and assessing the potential relevance of key obesity genes as therapeutic targets.

METHODS AND RESULTS

First, we generated a list of 28,787 obesity-associated SNPs from the publicly available GWAS dataset (approximately 800,000 individuals in the GIANT meta-analysis). Then, we prioritized 1372 genes with significant in silico evidence against genomic and transcriptomic data, including transcriptionally regulated genes in the brain from transcriptome-wide association studies. In further narrowing down the gene list, we selected key genes, which we found to be useful for the discovery of potential drug seeds as demonstrated in lipid GWAS separately. We thus identified 74 key genes for obesity, which are highly interconnected and enriched in several biological processes that contribute to obesity, including energy expenditure and homeostasis. Of 74 key genes, 37 had not been reported for the pathophysiology of obesity. Finally, by drug-gene interaction analysis, we detected 23 (of 74) key genes that are potential targets for 78 approved and marketed drugs.

CONCLUSIONS

Our results provide valuable insights into new treatment options for obesity through a data-driven approach that integrates multiple up-to-date knowledgebases.

Allelic variation in transcription factor PtoWRKY68 contributes to drought tolerance in Populus

Drought stress limits woody species productivity and influences tree distribution. However, dissecting the molecular mechanisms that underpin drought responses in forest trees can be challenging due to trait complexity. Here, using a panel of 300 Chinese white poplar (Populus tomentosa) accessions collected from different geographical climatic regions in China, we performed a genome-wide association study (GWAS) on seven drought-related traits and identified PtoWRKY68 as a candidate gene involved in the response to drought stress. A 12-bp insertion and/or deletion and three nonsynonymous variants in the PtoWRKY68 coding sequence categorized natural populations of P. tomentosa into two haplotype groups, PtoWRKY68hap1 and PtoWRKY68hap2. The allelic variation in these two PtoWRKY68 haplotypes conferred differential transcriptional regulatory activities and binding to the promoters of downstream abscisic acid (ABA) efflux and signaling genes. Overexpression of PtoWRKY68hap1 and PtoWRKY68hap2 in Arabidopsis (Arabidopsis thaliana) ameliorated the drought tolerance of two transgenic lines and increased ABA content by 42.7% and 14.3% compared to wild-type plants, respectively. Notably, PtoWRKY68hap1 (associated with drought tolerance) is ubiquitous in accessions in water-deficient environments, whereas the drought-sensitive allele PtoWRKY68hap2 is widely distributed in well-watered regions, consistent with the trends in local precipitation, suggesting that these alleles correspond to geographical adaptation in Populus. Moreover, quantitative trait loci analysis and an electrophoretic mobility shift assay showed that SHORT VEGETATIVE PHASE (PtoSVP.3) positively regulates the expression of PtoWRKY68 under drought stress. We propose a drought tolerance regulatory module in which PtoWRKY68 modulates ABA signaling and accumulation, providing insight into the genetic basis of drought tolerance in trees. Our findings will facilitate molecular breeding to improve the drought tolerance of forest trees.

Genome-wide association studies of human and rat BMI converge on synapse, epigenome, and hormone signaling networks

A vexing observation in genome-wide association studies (GWASs) is that parallel analyses in different species may not identify orthologous genes. Here, we demonstrate that cross-species translation of GWASs can be greatly improved by an analysis of co-localization within molecular networks. Using body mass index (BMI) as an example, we show that the genes associated with BMI in humans lack significant agreement with those identified in rats. However, the networks interconnecting these genes show substantial overlap, highlighting common mechanisms including synaptic signaling, epigenetic modification, and hormonal regulation. Genetic perturbations within these networks cause abnormal BMI phenotypes in mice, too, supporting their broad conservation across mammals. Other mechanisms appear species specific, including carbohydrate biosynthesis (humans) and glycerolipid metabolism (rodents). Finally, network co-localization also identifies cross-species convergence for height/body length. This study advances a general paradigm for determining whether and how phenotypes measured in model species recapitulate human biology.

COVID-GWAB: A Web-Based Prediction of COVID-19 Host Genes via Network Boosting of Genome-Wide Association Data

Host genetics affect both the susceptibility and response to viral infection. Searching for host genes that contribute to COVID-19, the Host Genetics Initiative (HGI) was formed to investigate the genetic factors involved in COVID-19 via genome-wide association studies (GWAS). The GWAS suffer from limited statistical power and in general, only a few genes can pass the conventional significance thresholds. This statistical limitation may be overcome by boosting weak association signals through integrating independent functional information such as molecular interactions. Additionally, the boosted results can be evaluated by various independent data for further connections to COVID-19. We present COVID-GWAB, a web-based tool to boost original GWAS signals from COVID-19 patients by taking the signals of the interactome neighbors. COVID-GWAB takes summary statistics from the COVID-19 HGI or user input data and reprioritizes candidate host genes for COVID-19 using HumanNet, a co-functional human gene network. The current version of COVID-GWAB provides the pre-processed data of releases 5, 6, and 7 of the HGI. Additionally, COVID-GWAB provides web interfaces for a summary of augmented GWAS signals, prediction evaluations by appearance frequency in COVID-19 literature, single-cell transcriptome data, and associated pathways. The web server also enables browsing the candidate gene networks.

From random to predictive: a context-specific interaction framework improves selection of drug protein-protein interactions for unknown drug pathways

With high drug attrition, protein-protein interaction (PPI) network models are attractive as efficient methods for predicting drug outcomes by analyzing proteins downstream of drug targets. Unfortunately, these methods tend to overpredict associations and they have low precision and prediction performance; performance is often no better than random (AUROC ~0.5). Typically, PPI models identify ranked phenotypes associated with downstream proteins, yet methods differ in prioritization of downstream proteins. Most methods apply global approaches for assessing all phenotypes. We hypothesized that a per-phenotype analysis could improve prediction performance. We compared two global approaches-statistical and distance-based-and our novel per-phenotype approach, 'context-specific interaction' (CSI) analysis, on severe side effect prediction. We used a novel dataset of adverse events (or designated medical events, DMEs) and discovered that CSI had a 50% improvement over global approaches (AUROC 0.77 compared to 0.51), and a 76-95% improvement in average precision (0.499 compared to 0.284, 0.256). Our results provide a quantitative rationale for considering downstream proteins on a per-phenotype basis when using PPI network methods to predict drug phenotypes.

The neurobiological basis of divergent thinking: Insight from gene co-expression network-based analysis

Although many efforts have been made to explore the genetic basis of divergent thinking (DT), there is still a gap in the understanding of how these findings relate to the neurobiology of DT. In a combined sample of 1,682 Chinese participants, by integrating GWAS with previously identified brain-specific gene co-expression network modules, this study explored for the first time the functional brain-specific gene co-expression networks underlying DT. The results showed that gene co-expression network modules in anterior cingulate cortex, caudate, amygdala and substantia nigra were enriched with DT association signals. Further functional enrichment analysis showed that these DT-related gene co-expression network modules were enriched for key biological process and cellular component related to myelination, suggesting that cortical and sub-cortical grey matter myelination may serve as important neurobiological basis of DT. Although the underlying mechanisms need to be further refined, this exploratory study may provide new insight into the neurobiology of DT.

Distinct effect of prenatal and postnatal brain expression across 20 brain disorders and anthropometric social traits: a systematic study of spatiotemporal modularity

Different spatiotemporal abnormalities have been implicated in different neuropsychiatric disorders and anthropometric social traits, yet an investigation in the temporal network modularity with brain tissue transcriptomics has been lacking. We developed a supervised network approach to investigate the genome-wide association study (GWAS) results in the spatial and temporal contexts and demonstrated it in 20 brain disorders and anthropometric social traits. BrainSpan transcriptome profiles were used to discover significant modules enriched with trait susceptibility genes in a developmental stage-stratified manner. We investigated whether, and in which developmental stages, GWAS-implicated genes are coordinately expressed in brain transcriptome. We identified significant network modules for each disorder and trait at different developmental stages, providing a systematic view of network modularity at specific developmental stages for a myriad of brain disorders and traits. Specifically, we observed a strong pattern of the fetal origin for most psychiatric disorders and traits [such as schizophrenia (SCZ), bipolar disorder, obsessive-compulsive disorder and neuroticism], whereas increased co-expression activities of genes were more strongly associated with neurological diseases [such as Alzheimer's disease (AD) and amyotrophic lateral sclerosis] and anthropometric traits (such as college completion, education and subjective well-being) in postnatal brains. Further analyses revealed enriched cell types and functional features that were supported and corroborated prior knowledge in specific brain disorders, such as clathrin-mediated endocytosis in AD, myelin sheath in multiple sclerosis and regulation of synaptic plasticity in both college completion and education. Our study provides a landscape view of the spatiotemporal features in a myriad of brain-related disorders and traits.

A computational approach for detecting physiological homogeneity in the midst of genetic heterogeneity

The human genetic dissection of clinical phenotypes is complicated by genetic heterogeneity. Gene burden approaches that detect genetic signals in case-control studies are underpowered in genetically heterogeneous cohorts. We therefore developed a genome-wide computational method, network-based heterogeneity clustering (NHC), to detect physiological homogeneity in the midst of genetic heterogeneity. Simulation studies showed our method to be capable of systematically converging genes in biological proximity on the background biological interaction network, and capturing gene clusters harboring presumably deleterious variants, in an efficient and unbiased manner. We applied NHC to whole-exome sequencing data from a cohort of 122 individuals with herpes simplex encephalitis (HSE), including 13 individuals with previously published monogenic inborn errors of TLR3-dependent IFN-α/β immunity. The top gene cluster identified by our approach successfully detected and prioritized all causal variants of five TLR3 pathway genes in the 13 previously reported individuals. This approach also suggested candidate variants of three reported genes and four candidate genes from the same pathway in another ten previously unstudied individuals. TLR3 responsiveness was impaired in dermal fibroblasts from four of the five individuals tested, suggesting that the variants detected were causal for HSE. NHC is, therefore, an effective and unbiased approach for unraveling genetic heterogeneity by detecting physiological homogeneity.

EyeDiseases: an integrated resource for dedicating to genetic variants, gene expression and epigenetic factors of human eye diseases

Eye diseases are remarkably common and encompass a large and diverse range of morbidities that affect different components of the visual system and visual function. With advances in omics technology of eye disorders, genome-scale datasets have been rapidly accumulated in genetics and epigenetics field. However, the efficient collection and comprehensive analysis of different kinds of omics data are lacking. Herein, we developed EyeDiseases (https://eyediseases.bio-data.cn/), the first database for multi-omics data integration and interpretation of human eyes diseases. It contains 1344 disease-associated genes with genetic variation, 1774 transcription files of bulk cell expression and single-cell RNA-seq, 105 epigenomics data across 185 kinds of human eye diseases. Using EyeDiseases, we investigated SARS-CoV-2 potential tropism in eye infection and found that the SARS-CoV-2 entry factors, ACE2 and TMPRSS2 are highly correlated with cornea and keratoconus, suggest that ocular surface cells are susceptible to infection by SARS-CoV-2. Additionally, integrating analysis of Age-related macular degeneration (AMD) GWAS loci and co-expression data revealed 9 associated genes involved in HIF-1 signaling pathway and voltage-gate potassium channel complex. The EyeDiseases provides a valuable resource for accelerating the discovery and validation of candidate loci and genes contributed to the molecular diagnosis and therapeutic vulnerabilities with various eyes diseases.

Pathway analysis for genome-wide genetic variation data: Analytic principles, latest developments, and new opportunities

Pathway analysis, also known as gene-set enrichment analysis, is a multilocus analytic strategy that integrates a priori, biological knowledge into the statistical analysis of high-throughput genetics data. Originally developed for the studies of gene expression data, it has become a powerful analytic procedure for in-depth mining of genome-wide genetic variation data. Astonishing discoveries were made in the past years, uncovering genes and biological mechanisms underlying common and complex disorders. However, as massive amounts of diverse functional genomics data accrue, there is a pressing need for newer generations of pathway analysis methods that can utilize multiple layers of high-throughput genomics data. In this review, we provide an intellectual foundation of this powerful analytic strategy, as well as an update of the state-of-the-art in recent method developments. The goal of this review is threefold: (1) introduce the motivation and basic steps of pathway analysis for genome-wide genetic variation data; (2) review the merits and the shortcomings of classic and newly emerging integrative pathway analysis tools; and (3) discuss remaining challenges and future directions for further method developments.

A Review of Statistical Methods for Identifying Trait-Relevant Tissues and Cell Types

Genome-wide association studies (GWASs) have identified and replicated many genetic variants that are associated with diseases and disease-related complex traits. However, the biological mechanisms underlying these identified associations remain largely elusive. Exploring the biological mechanisms underlying these associations requires identifying trait-relevant tissues and cell types, as genetic variants likely influence complex traits in a tissue- and cell type-specific manner. Recently, several statistical methods have been developed to integrate genomic data with GWASs for identifying trait-relevant tissues and cell types. These methods often rely on different genomic information and use different statistical models for trait-tissue relevance inference. Here, we present a comprehensive technical review to summarize ten existing methods for trait-tissue relevance inference. These methods make use of different genomic information that include functional annotation information, expression quantitative trait loci information, genetically regulated gene expression information, as well as gene co-expression network information. These methods also use different statistical models that range from linear mixed models to covariance network models. We hope that this review can serve as a useful reference both for methodologists who develop methods and for applied analysts who apply these methods for identifying trait relevant tissues and cell types.

A framework for pathway knowledge driven prioritization in genome-wide association studies

Many variants with low frequencies or with low to modest effects likely remain unidentified in genome-wide association studies (GWAS) because of stringent genome-wide thresholds for detection. To improve the power of detection, variant prioritization based on their functional annotations and epigenetic landmarks has been used successfully. Here, we propose a novel method of prioritization of a GWAS by exploiting gene-level knowledge (e.g., annotations to pathways and ontologies) and show that it further improves power. Often, disease associated variants are found near genes that are coinvolved in specific biological pathways relevant to disease process. Utilization of this knowledge to conduct a prioritized scan increases the power to detect loci that map to genes clustered in a few specific pathways. We have developed a computationally scalable framework based on penalized logistic regression (termed GKnowMTest-Genomic Knowledge-guided Multiplte Testing) to enable a prioritized pathway-guided GWAS scan with a very large number of gene-level annotations. We demonstrate that the proposed strategy improves overall power and maintains the Type 1 error globally. Our method works on genome-wide summary level data and a user-specified list of pathways (e.g., those extracted from large pathway databases without reference to biology of a specific disease). It automatically reweights the input p values by incorporating the pathway enrichments as "adaptively learned" from the data using a cross-validation technique to avoid overfitting. We used whole-genome simulations and some publicly available GWAS data sets to illustrate the application of our method. The GKnowMTest framework has been implemented as a user-friendly open-source R package.

Genetic Basis of Maize Resistance to Multiple Insect Pests: Integrated Genome-Wide Comparative Mapping and Candidate Gene Prioritization

Several species of herbivores feed on maize in field and storage setups, making the development of multiple insect resistance a critical breeding target. In this study, an association mapping panel of 341 tropical maize lines was evaluated in three field environments for resistance to fall armyworm (FAW), whilst bulked grains were subjected to a maize weevil (MW) bioassay and genotyped with Diversity Array Technology's single nucleotide polymorphisms (SNPs) markers. A multi-locus genome-wide association study (GWAS) revealed 62 quantitative trait nucleotides (QTNs) associated with FAW and MW resistance traits on all 10 maize chromosomes, of which, 47 and 31 were discovered at stringent Bonferroni genome-wide significance levels of 0.05 and 0.01, respectively, and located within or close to multiple insect resistance genomic regions (MIRGRs) concerning FAW, SB, and MW. Sixteen QTNs influenced multiple traits, of which, six were associated with resistance to both FAW and MW, suggesting a pleiotropic genetic control. Functional prioritization of candidate genes (CGs) located within 10-30 kb of the QTNs revealed 64 putative GWAS-based CGs (GbCGs) showing evidence of involvement in plant defense mechanisms. Only one GbCG was associated with each of the five of the six combined resistance QTNs, thus reinforcing the pleiotropy hypothesis. In addition, through in silico co-functional network inferences, an additional 107 network-based CGs (NbCGs), biologically connected to the 64 GbCGs, and differentially expressed under biotic or abiotic stress, were revealed within MIRGRs. The provided multiple insect resistance physical map should contribute to the development of combined insect resistance in maize.

High-Resolution Regulatory Maps Connect Vascular Risk Variants to Disease-Related Pathways

Supplemental Digital Content is available in the text.

Genetic variant landscape of coronary artery disease is dominated by noncoding variants among which many occur within putative enhancers regulating the expression levels of relevant genes. It is crucial to assign the genetic variants to their correct genes both to gain insights into perturbed functions and better assess the risk of disease.

The Integrated Landscape of Biological Candidate Causal Genes in Coronary Artery Disease

Background

Genome-wide association studies (GWASs) have identified more than 150 genetic loci that demonstrate robust association with coronary artery disease (CAD). In contrast to the success of GWAS, the translation from statistical signals to biological mechanism and exploration of causal genes for drug development remain difficult, owing to the complexity of gene regulatory and linkage disequilibrium patterns. We aim to prioritize the plausible causal genes for CAD at a genome-wide level.

Methods

We integrated the latest GWAS summary statistics with other omics data from different layers and utilized eight different computational methods to predict CAD potential causal genes. The prioritized candidate genes were further characterized by pathway enrichment analysis, tissue-specific expression analysis, and pathway crosstalk analysis.

Results

Our analysis identified 55 high-confidence causal genes for CAD, among which 15 genes (LPL, COL4A2, PLG, CDKN2B, COL4A1, FES, FLT1, FN1, IL6R, LPA, PCSK9, PSRC1, SMAD3, SWAP70, and VAMP8) ranked the highest priority because of consistent evidence from different data-driven approaches. GO analysis showed that these plausible causal genes were enriched in lipid metabolic and extracellular regions. Tissue-specific enrichment analysis revealed that these genes were significantly overexpressed in adipose and liver tissues. Further, KEGG and crosstalk analysis also revealed several key pathways involved in the pathogenesis of CAD.

Conclusion

Our study delineated the landscape of CAD potential causal genes and highlighted several biological processes involved in CAD pathogenesis. Further studies and experimental validations of these genes may shed light on mechanistic insights into CAD development and provide potential drug targets for future therapeutics.

Dense module searching for gene networks associated with multiple sclerosis

BACKGROUND

Multiple sclerosis (MS) is a complex disease in which the immune system attacks the central nervous system. The molecular mechanisms contributing to the etiology of MS remain poorly understood. Genome-wide association studies (GWAS) of MS have identified a small number of genetic loci significant at the genome level, but they are mainly non-coding variants. Network-assisted analysis may help better interpret the functional roles of the variants with association signals and potential translational medicine application. The Dense Module Searching of GWAS tool (dmGWAS version 2.4) developed in our team is applied to 2 MS GWAS datasets (GeneMSA and IMSGC GWAS) using the human protein interactome as the reference network. A dual evaluation strategy is used to generate results with reproducibility.

RESULTS

Approximately 7500 significant network modules were identified for each independent GWAS dataset, and 20 significant modules were identified from the dual evaluation. The top modules included GRB2, HDAC1, JAK2, MAPK1, and STAT3 as central genes. Top module genes were enriched with functional terms such as "regulation of glial cell differentiation" (adjusted p-value = 2.58 × 10- 3), "T-cell costimulation" (adjusted p-value = 2.11 × 10- 6) and "virus receptor activity" (adjusted p-value = 1.67 × 10- 3). Interestingly, top gene networks included several MS FDA approved drug target genes HDAC1, IL2RA, KEAP1, and RELA, CONCLUSIONS: Our dmGWAS network analyses highlighted several genes (GRB2, HDAC1, IL2RA, JAK2, KEAP1, MAPK1, RELA and STAT3) in top modules that are promising to interpret GWAS signals and link to MS drug targets. The genes enriched with glial cell differentiation are important for understanding neurodegenerative processes in MS and for remyelination therapy investigation. Importantly, our identified genetic signals enriched in T cell costimulation and viral receptor activity supported the viral infection onset hypothesis for MS.

Leveraging gene co-expression patterns to infer trait-relevant tissues in genome-wide association studies

Genome-wide association studies (GWASs) have identified many SNPs associated with various common diseases. Understanding the biological functions of these identified SNP associations requires identifying disease/trait relevant tissues or cell types. Here, we develop a network method, CoCoNet, to facilitate the identification of trait-relevant tissues or cell types. Different from existing approaches, CoCoNet incorporates tissue-specific gene co-expression networks constructed from either bulk or single cell RNA sequencing (RNAseq) studies with GWAS data for trait-tissue inference. In particular, CoCoNet relies on a covariance regression network model to express gene-level effect measurements for the given GWAS trait as a function of the tissue-specific co-expression adjacency matrix. With a composite likelihood-based inference algorithm, CoCoNet is scalable to tens of thousands of genes. We validate the performance of CoCoNet through extensive simulations. We apply CoCoNet for an in-depth analysis of four neurological disorders and four autoimmune diseases, where we integrate the corresponding GWASs with bulk RNAseq data from 38 tissues and single cell RNAseq data from 10 cell types. In the real data applications, we show how CoCoNet can help identify specific glial cell types relevant for neurological disorders and identify disease-targeted colon tissues as relevant for autoimmune diseases.

Using networks to identify and interpret natural variation

Studies on natural variation and network biology inherently work to summarize vast amounts of information and data. The combination of these two areas of study while creating datasets of immense complexity is critical to their mutual progress. Networks are necessary as a way to work to reduce the dimensionality inherent in natural variation with 100 s to 1000 s of genotypes. Correspondingly natural variation is essential for testing how networks may or may not be shared across individuals or species. Advances in this area of cross-fertilization including using networks directly as phenotypes and the use of networks to help in prioritizing candidate gene validation efforts. Interesting new observations on frequent presence-absence variation in gene content and adaptation is beginning to highlight the potential for natural variation in network presence-absence. This review attempts to delve into these new insights.

Identification of disease-associated loci using machine learning for genotype and network data integration

MOTIVATION

Integration of different omics data could markedly help to identify biological signatures, understand the missing heritability of complex diseases and ultimately achieve personalized medicine. Standard regression models used in Genome-Wide Association Studies (GWAS) identify loci with a strong effect size, whereas GWAS meta-analyses are often needed to capture weak loci contributing to the missing heritability. Development of novel machine learning algorithms for merging genotype data with other omics data is highly needed as it could enhance the prioritization of weak loci.

RESULTS

We developed cNMTF (corrected non-negative matrix tri-factorization), an integrative algorithm based on clustering techniques of biological data. This method assesses the inter-relatedness between genotypes, phenotypes, the damaging effect of the variants and gene networks in order to identify loci-trait associations. cNMTF was used to prioritize genes associated with lipid traits in two population cohorts. We replicated 129 genes reported in GWAS world-wide and provided evidence that supports 85% of our findings (226 out of 265 genes), including recent associations in literature (NLGN1), regulators of lipid metabolism (DAB1) and pleiotropic genes for lipid traits (CARM1). Moreover, cNMTF performed efficiently against strong population structures by accounting for the individuals' ancestry. As the method is flexible in the incorporation of diverse omics data sources, it can be easily adapted to the user's research needs.

AVAILABILITY AND IMPLEMENTATION

An R package (cnmtf) is available at https://lgl15.github.io/cnmtf_web/index.html.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Benchmarking network propagation methods for disease gene identification

In-silico identification of potential target genes for disease is an essential aspect of drug target discovery. Recent studies suggest that successful targets can be found through by leveraging genetic, genomic and protein interaction information. Here, we systematically tested the ability of 12 varied algorithms, based on network propagation, to identify genes that have been targeted by any drug, on gene-disease data from 22 common non-cancerous diseases in OpenTargets. We considered two biological networks, six performance metrics and compared two types of input gene-disease association scores. The impact of the design factors in performance was quantified through additive explanatory models. Standard cross-validation led to over-optimistic performance estimates due to the presence of protein complexes. In order to obtain realistic estimates, we introduced two novel protein complex-aware cross-validation schemes. When seeding biological networks with known drug targets, machine learning and diffusion-based methods found around 2-4 true targets within the top 20 suggestions. Seeding the networks with genes associated to disease by genetics decreased performance below 1 true hit on average. The use of a larger network, although noisier, improved overall performance. We conclude that diffusion-based prioritisers and machine learning applied to diffusion-based features are suited for drug discovery in practice and improve over simpler neighbour-voting methods. We also demonstrate the large impact of choosing an adequate validation strategy and the definition of seed disease genes.

The use of biological network data has proven its effectiveness in many areas from computational biology. Networks consist of nodes, usually genes or proteins, and edges that connect pairs of nodes, representing information such as physical interactions, regulatory roles or co-occurrence. In order to find new candidate nodes for a given biological property, the so-called network propagation algorithms start from the set of known nodes with that property and leverage the connections from the biological network to make predictions. Here, we assess the performance of several network propagation algorithms to find sensible gene targets for 22 common non-cancerous diseases, i.e. those that have been found promising enough to start the clinical trials with any compound. We focus on obtaining performance metrics that reflect a practical scenario in drug development where only a small set of genes can be essayed. We found that the presence of protein complexes biased the performance estimates, leading to over-optimistic conclusions, and introduced two novel strategies to address it. Our results support that network propagation is still a viable approach to find drug targets, but that special care needs to be put on the validation strategy. Algorithms benefitted from the use of a larger -although noisier- network and of direct evidence data, rather than indirect genetic associations to disease.

Plant Networks as Traits and Hypotheses: Moving Beyond Description

Biology relies on the central thesis that the genes in an organism encode molecular mechanisms that combine with stimuli and raw materials from the environment to create a final phenotypic expression representative of the genomic programming. While conceptually simple, the genotype-to-phenotype linkage in a eukaryotic organism relies on the interactions of thousands of genes and an environment with a potentially unknowable level of complexity. Modern biology has moved to the use of networks in systems biology to try to simplify this complexity to decode how an organism's genome works. Previously, biological networks were basic ways to organize, simplify, and analyze data. However, recent advances are allowing networks to move beyond description and become phenotypes or hypotheses in their own right. This review discusses these efforts, like mapping responses across biological scales, including relationships among cellular entities, and the direct use of networks as traits or hypotheses.

Benchmarker: An Unbiased, Association-Data-Driven Strategy to Evaluate Gene Prioritization Algorithms

Genome-wide association studies (GWASs) are valuable for understanding human biology, but associated loci typically contain multiple associated variants and genes. Thus, algorithms that prioritize likely causal genes and variants for a given phenotype can provide biological interpretations of association data. However, a critical, currently missing capability is to objectively compare performance of such algorithms. Typical comparisons rely on "gold standard" genes harboring causal coding variants, but such gold standards may be biased and incomplete. To address this issue, we developed Benchmarker, an unbiased, data-driven benchmarking method that compares performance of similarity-based prioritization strategies to each other (and to random chance) by leave-one-chromosome-out cross-validation with stratified linkage disequilibrium (LD) score regression. We first applied Benchmarker to 20 well-powered GWASs and compared gene prioritization based on strategies employing three different data sources, including annotated gene sets and gene expression; genes prioritized based on gene sets had higher per-SNP heritability than those prioritized based on gene expression. Additionally, in a direct comparison of three methods, DEPICT and MAGMA outperformed NetWAS. We also evaluated combinations of methods; our results indicated that combining data sources and algorithms can help prioritize higher-quality genes for follow-up. Benchmarker provides an unbiased approach to evaluate any similarity-based method that provides genome-wide prioritization of genes, variants, or gene sets and can determine the best such method for any particular GWAS. Our method addresses an important unmet need for rigorous tool assessment and can assist in mapping genetic associations to causal function.

A Fast and Flexible Framework for Network-Assisted Genomic Association

We present an accessible, fast, and customizable network propagation system for pathway boosting and interpretation of genome-wide association studies. This system-NAGA (Network Assisted Genomic Association)-taps the NDEx biological network resource to gain access to thousands of protein networks and select those most relevant and performative for a specific association study. The method works efficiently, completing genome-wide analysis in under 5 minutes on a modern laptop computer. We show that NAGA recovers many known disease genes from analysis of schizophrenia genetic data, and it substantially boosts associations with previously unappreciated genes such as amyloid beta precursor. On this and seven other gene-disease association tasks, NAGA outperforms conventional approaches in recovery of known disease genes and replicability of results. Protein interactions associated with disease are visualized and annotated in Cytoscape, which, in addition to standard programmatic interfaces, allows for downstream analysis.

MicroRNA-655-3p and microRNA-497-5p inhibit cell proliferation in cultured human lip cells through the regulation of genes related to human cleft lip

Background

The etiology of cleft lip with or without cleft palate (CL/P), a common congenital birth defect, is complex and involves the contribution of genetic and environmental factors. Although many candidate genes have been identified, the regulation and interaction of these genes in CL/P remain unclear. In addition, the contribution of microRNAs (miRNAs), non-coding RNAs that regulate the expression of multiple genes, to the etiology of CL/P is largely unknown.

Methods

To identify the signatures of causative biological pathways for human CL/P, we conducted a systematic literature review for human CL/P candidate genes and subsequent bioinformatics analyses. Functional enrichment analyses of the candidate CL/P genes were conducted using the pathway databases GO and KEGG. The miRNA-mediated post-transcriptional regulation of the CL/P candidate genes was analyzed with miRanda, PITA, and TargetScan, and miRTarbase. Genotype-phenotype association analysis was conducted using GWAS. The functional significance of the candidate miRNAs was evaluated experimentally in cell proliferation and target gene regulation assays in human lip fibroblasts.

Results

Through an extensive search of the main biomedical databases, we mined 177 genes with mutations or association/linkage reported in individuals with CL/P, and considered them as candidate genes for human CL/P. The genotype-phenotype association study revealed that mutations in 12 genes (ABCA4, ADAM3A, FOXE1, IRF6, MSX2, MTHFR, NTN1, PAX7, TP63, TPM1, VAX1, and WNT9B) were significantly associated with CL/P. In addition, our bioinformatics analysis predicted 16 microRNAs (miRNAs) to be post-transcriptional regulators of CL/P genes. To validate the bioinformatics results, the top six candidate miRNAs (miR-124-3p, miR-369-3p, miR-374a-5p, miR-374b-5p, miR-497-5p, and miR-655-3p) were evaluated by cell proliferation/survival assays and miRNA-gene regulation assays in cultured human lip fibroblasts. We found that miR-497-5p and miR-655-3p significantly suppressed cell proliferation in these cells. Furthermore, the expression of the predicted miRNA-target genes was significantly downregulated by either miR-497-5p or miR-655-3p mimic.

Conclusion

Expression of miR-497-5p and miR-655-3p suppresses cell proliferation through the regulation of human CL/P-candidate genes. This study provides insights into the role of miRNAs in the etiology of CL/P and suggests possible strategies for the diagnosis of CL/P.

Electronic supplementary material

The online version of this article (10.1186/s12920-019-0535-2) contains supplementary material, which is available to authorized users.

Translational bioinformatics in mental health: open access data sources and computational biomarker discovery

Mental illness is increasingly recognized as both a significant cost to society and a significant area of opportunity for biological breakthrough. As -omics and imaging technologies enable researchers to probe molecular and physiological underpinnings of multiple diseases, opportunities arise to explore the biological basis for behavioral health and disease. From individual investigators to large international consortia, researchers have generated rich data sets in the area of mental health, including genomic, transcriptomic, metabolomic, proteomic, clinical and imaging resources. General data repositories such as the Gene Expression Omnibus (GEO) and Database of Genotypes and Phenotypes (dbGaP) and mental health (MH)-specific initiatives, such as the Psychiatric Genomics Consortium, MH Research Network and PsychENCODE represent a wealth of information yet to be gleaned. At the same time, novel approaches to integrate and analyze data sets are enabling important discoveries in the area of mental and behavioral health. This review will discuss and catalog into an organizing framework the increasingly diverse set of MH data resources available, using schizophrenia as a focus area, and will describe novel and integrative approaches to molecular biomarker discovery that make use of mental health data.

Network and Pathway-Based Analysis of Single-Nucleotide Polymorphism of miRNA in Temporal Lobe Epilepsy

Temporal lobe epilepsy (TLE) is a complex disease with its pathogenetic mechanism still unclear. Single-nucleotide polymorphisms (SNPs) of miRNA (miRSNPs) are SNPs located on miRNA genes or target sites of miRNAs, which have been proved to be associated with neuropsychic disease development by interfering with miRNA-mediated regulatory function. In this study, we integrated TLE-related risk genes and risk pathways multi-dimensionally based on public data resources. Furthermore, we systematically screened candidate functional miRSNPs for TLE and constructed a TLE-associated pathway-based miRSNP switching network, which included 92 miRNAs that target 12 TLE risk pathways. Moreover, we dissected thoroughly the correlation between 5 risk genes of 4 risk pathways and TLE development. Additionally, the biological function of several candidate miRSNPs were validated by luciferase reporter assay. In silico approach facilitates to select potential "miRSNP-miRNA-risk gene-pathway" axis for experimental validation, which provided new insights into the mechanism of miRSNPs as potential genetic risk factors of TLE.

Puberty and the Evolution of Developmental Science

In recent decades, theoretical and methodological advances have operated synergistically to advance understanding of puberty and prompt increasingly comprehensive models that engage with the temporal, psychosocial, and biological dimensions of this maturational milepost. This integrative overview discusses these theoretical and methodological advances and their implications for research and intervention to promote human development in the context of changing maturational schedules and massive ongoing social transformations.

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.

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.

Mouse Models and Online Resources for Functional Analysis of Osteoporosis Genome-Wide Association Studies

Osteoporosis is a complex genetic disease in which the number of loci associated with the bone mineral density, a clinical risk factor for fracture, has increased at an exponential rate in the last decade. The identification of the causative variants and candidate genes underlying these loci has not been able to keep pace with the rate of locus discovery. A large number of tools and data resources have been built around the use of the mouse as model of human genetic disease. Herein, we describe resources available for functional validation of human Genome Wide Association Study (GWAS) loci using mouse models. We specifically focus on large-scale phenotyping efforts focused on bone relevant phenotypes and repositories of genotype-phenotype data that exist for transgenic and mutant mice, which can be readily mined as a first step toward more targeted efforts designed to deeply characterize the role of a gene in bone biology.

Re-thinking Alzheimer's disease therapeutic targets using gene-based tests

Background

Alzheimer's disease (AD) is a devastating condition with no known effective drug treatments. Existing drugs only alleviate symptoms. Given repeated expensive drug failures, we assessed systematically whether approved and investigational AD drugs are targeting products of genes strongly associated with AD and whether these genes are targeted by existing drugs for other indications which could be re-purposed.

Methods

We identified genes strongly associated with late-onset AD from the loci of genetic variants associated with AD at genome-wide-significance and from a gene-based test applied to the most extensively genotyped late-onset AD case (n = 17,008)-control (n = 37,154) study, the International Genomics of Alzheimer's Project. We used three gene-to-drug cross-references, Kyoto Encyclopedia of Genes and Genomes, Drugbank and Drug Repurposing Hub, to identify genetically validated targets of AD drugs and any existing drugs or nutraceuticals targeting products of the genes strongly associated with late-onset AD.

Findings

A total of 67 autosomal genes (forming 9 gene clusters) were identified as strongly associated with late-onset AD, 28 from the loci of single genetic variants, 51 from the gene-based test and 12 by both methods. Existing approved or investigational AD drugs did not target products of any of these 67 genes. Drugs for other indications targeted 11 of these genes, including immunosuppressive disease-modifying anti-rheumatic drugs targeting PTK2B gene products.

Interpretation

Approved and investigational AD drugs are not targeting products of genes strongly associated with late-onset AD. However, other drugs targeting products of these genes exist and could perhaps be re-purposing to combat late-onset AD after further scrutiny.

Integrating RNA-seq and GWAS reveals novel genetic mutations for buffalo reproductive traits

Genome-wide association study (GWAS) has been applied in buffalo breeding programs and been used to identify a number of candidate genes associated with buffalo reproductive traits. The genetic code of specific genes underlying buffalo reproductive traits remains unclear. Association study that measures both genetic and transcriptional variation has been applied for the investigation of complex traits. To investigate genes involved in buffalo reproductive traits, integrated RNA-seq results were investigated of buffalo granulosa cells and candidate genes which were reported to be associated with buffalo reproductive traits in a previous GWAS. A large number of variants were detected by RNA-seq, and 214 variants were located within the buffalo reproductive candidate genes identified by GWAS. A further association study in 462 Italian Mediterranean buffalo indicated that 25 SNPs distributed in 13 genes were associated with reproductive traits. Of the 13 genes, 11 were expressed in granulosa cells of all antral follicle development stages, and significant difference was found in the expression of NDUFS2 between follicles of diameter <8 mm and > 8 mm. These findings extend the results of GWAS by expanding the knowledge about new and potentially functional single-nucleotide polymorphisms and provide useful information about regulatory genes affecting buffalo reproductive traits.

Network-based genetic investigation of virulence-associated phenotypes in methicillin-resistant Staphylococcus aureus

Staphylococcus aureus is a gram-positive bacterium that causes a wide range of infections. Recently, the spread of methicillin-resistant S. aureus (MRSA) strains has seriously reduced antibiotic treatment options. Anti-virulence strategies, the objective of which is to target the virulence instead of the viability of the pathogen, have become widely accepted as a means of avoiding the emergence of new antibiotic-resistant strains. To increase the number of anti-virulence therapeutic options, it is necessary to identify as many novel virulence-associated genes as possible in MRSA. Co-functional networks have proved useful for mapping gene-to-phenotype associations in various organisms. Herein, we present StaphNet (www.inetbio.org/staphnet), a genome-scale co-functional network for an MRSA strain, S. aureus subsp. USA300_FPR3757. StaphNet, which was constructed by the integration of seven distinct types of genomics data within a Bayesian statistics framework, covers approximately 94% of the coding genome with a high degree of accuracy. We implemented a companion web server for network-based gene prioritization of the phenotypes of 31 different S. aureus strains. We demonstrated that StaphNet can effectively identify genes for virulence-associated phenotypes in MRSA. These results suggest that StaphNet can facilitate target discovery for the development of anti-virulence drugs to treat MRSA infection.

Co-expression network of transcription factors reveal ethylene-responsive element-binding factor as key regulator of wood phenotype in Eucalyptus tereticornis

Suitability of wood biomass for pulp production is dependent on the cellular architecture and composition of secondary cell wall. Presently, systems genetics approach is being employed to understand the molecular basis of trait variation and co-expression network analysis has enabled holistic understanding of complex trait such as secondary development. Transcription factors (TFs) are reported as key regulators of meristematic growth and wood formation. The hierarchical TF network is a multi-layered system which interacts with downstream structural genes involved in biosynthesis of cellulose, hemicelluloses and lignin. Several TFs have been associated with wood formation in tree species such as Populus, Eucalyptus, Picea and Pinus. However, TF-specific co-expression networks to understand the interaction between these regulators are not reported. In the present study, co-expression network was developed for TFs expressed during wood formation in Eucalyptus tereticornis and ethylene-responsive element-binding factor, EtERF2, was identified as the major hub transcript which co-expressed with other secondary cell wall biogenesis-specific TFs such as EtSND2, EtVND1, EtVND4, EtVND6, EtMYB70, EtGRAS and EtSCL8. This study reveals a probable role of ethylene in determining natural variation in wood properties in Eucalyptus species. Understanding this transcriptional regulation underpinning the complex bio-processing trait of wood biomass will complement the Eucalyptus breeding program through selection of industrially suitable phenotypes by marker-assisted selection.

One for all and all for One: Improving replication of genetic studies through network diffusion

Improving accuracy in genetic studies would greatly accelerate understanding the genetic basis of complex diseases. One approach to achieve such an improvement for risk variants identified by the genome wide association study (GWAS) approach is to incorporate previously known biology when screening variants across the genome. We developed a simple approach for improving the prioritization of candidate disease genes that incorporates a network diffusion of scores from known disease genes using a protein network and a novel integration with GWAS risk scores, and tested this approach on a large Alzheimer disease (AD) GWAS dataset. Using a statistical bootstrap approach, we cross-validated the method and for the first time showed that a network approach improves the expected replication rates in GWAS studies. Several novel AD genes were predicted including CR2, SHARPIN, and PTPN2. Our re-prioritized results are enriched for established known AD-associated biological pathways including inflammation, immune response, and metabolism, whereas standard non-prioritized results were not. Our findings support a strategy of considering network information when investigating genetic risk factors.

Integrating multiple types of -omics data is a rapidly growing research area due in part to the increasing amount of diverse and publicly accessible data. In this study, we demonstrated that integration of genetic association and protein interaction data using a network diffusion approach measurably improves reproducibility of top candidate genes. Application of this approach to Alzheimer disease (AD) using a large dataset assembled by the Alzheimer’s Disease Genetics Consortium identified several novel candidate AD genes that are supported by pre-existing knowledge of AD pathobiology. Our findings support a strategy of considering network information when investigating genetic risk factors. Finally, we developed a transparent and easy-to-use R package that can facilitate the extension of our methodology to other phenotypes for which genetic data are available.

Convergence between biological, behavioural and genetic determinants of obesity

Multiple biological, behavioural and genetic determinants or correlates of obesity have been identified to date. Genome-wide association studies (GWAS) have contributed to the identification of more than 100 obesity-associated genetic variants, but their roles in causal processes leading to obesity remain largely unknown. Most variants are likely to have tissue-specific regulatory roles through joint contributions to biological pathways and networks, through changes in gene expression that influence quantitative traits, or through the regulation of the epigenome. The recent availability of large-scale functional genomics resources provides an opportunity to re-examine obesity GWAS data to begin elucidating the function of genetic variants. Interrogation of knockout mouse phenotype resources provides a further avenue to test for evidence of convergence between genetic variation and biological or behavioural determinants of obesity.

GWAB: a web server for the network-based boosting of human genome-wide association data

During the last decade, genome-wide association studies (GWAS) have represented a major approach to dissect complex human genetic diseases. Due in part to limited statistical power, most studies identify only small numbers of candidate genes that pass the conventional significance thresholds (e.g. P ≤ 5 × 10-8). This limitation can be partly overcome by increasing the sample size, but this comes at a higher cost. Alternatively, weak association signals can be boosted by incorporating independent data. Previously, we demonstrated the feasibility of boosting GWAS disease associations using gene networks. Here, we present a web server, GWAB (www.inetbio.org/gwab), for the network-based boosting of human GWAS data. Using GWAS summary statistics (P-values) for SNPs along with reference genes for a disease of interest, GWAB reprioritizes candidate disease genes by integrating the GWAS and network data. We found that GWAB could more effectively retrieve disease-associated reference genes than GWAS could alone. As an example, we describe GWAB-boosted candidate genes for coronary artery disease and supporting data in the literature. These results highlight the inherent value in sub-threshold GWAS associations, which are often not publicly released. GWAB offers a feasible general approach to boost such associations for human disease genetics.

Mouse Genome-Wide Association Study of Preclinical Group II Pulmonary Hypertension Identifies Epidermal Growth Factor Receptor

Pulmonary hypertension (PH) is associated with features of obesity and metabolic syndrome that translate to the induction of PH by chronic high-fat diet (HFD) in some inbred mouse strains. We conducted a genome-wide association study (GWAS) to identify candidate genes associated with susceptibility to HFD-induced PH. Mice from 36 inbred and wild-derived strains were fed with regular diet or HFD for 20 weeks beginning at 6-12 weeks of age, after which right ventricular (RV) and left ventricular (LV) end-systolic pressure (ESP) and maximum pressure (MaxP) were measured by cardiac catheterization. We tested for association of RV MaxP and RV ESP and identified genomic regions enriched with nominal associations to both of these phenotypes. We excluded genomic regions if they were also associated with LV MaxP, LV ESP, or body weight. Genes within significant regions were scored based on the shortest-path betweenness centrality, a measure of network connectivity, of their human orthologs in a gene interaction network of human PH-related genes. WSB/EiJ, NON/ShiLtJ, and AKR/J mice had the largest increases in RV MaxP after high-fat feeding. Network-based scoring of GWAS candidates identified epidermal growth factor receptor (Egfr) as having the highest shortest-path betweenness centrality of GWAS candidates. Expression studies of lung homogenate showed that EGFR expression is increased in the AKR/J strain, which developed a significant increase in RV MaxP after high-fat feeding as compared with C57BL/6J, which did not. Our combined GWAS and network-based approach adds evidence for a role for Egfr in murine PH.

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.