GRAND: a database of gene regulatory network models across human conditions

Knowledge-guided learning methods for integrative analysis of multi-omics data

Integrative analysis of multi-omics data has the potential to yield valuable and comprehensive insights into the molecular mechanisms underlying complex diseases such as cancer and Alzheimer's disease. However, a number of analytical challenges complicate multi-omics data integration. For instance, -omics data are usually high-dimensional, and sample sizes in multi-omics studies tend to be modest. Furthermore, when genes in an important pathway have relatively weak signal, it can be difficult to detect them individually. There is a growing body of literature on knowledge-guided learning methods that can address these challenges by incorporating biological knowledge such as functional genomics and functional proteomics into multi-omics data analysis. These methods have been shown to outperform their counterparts that do not utilize biological knowledge in tasks including prediction, feature selection, clustering, and dimension reduction. In this review, we survey recently developed methods and applications of knowledge-guided multi-omics data integration methods and discuss future research directions.

Network-based drug repurposing for schizophrenia

Despite recent progress, the challenges in drug discovery for schizophrenia persist. However, computational drug repurposing has gained popularity as it leverages the wealth of expanding biomedical databases. Network analyses provide a comprehensive understanding of transcription factor (TF) regulatory effects through gene regulatory networks, which capture the interactions between TFs and target genes by integrating various lines of evidence. Using the PANDA algorithm, we examined the topological variances in TF-gene regulatory networks between individuals with schizophrenia and healthy controls. This algorithm incorporates binding motifs, protein interactions, and gene co-expression data. To identify these differences, we subtracted the edge weights of the healthy control network from those of the schizophrenia network. The resulting differential network was then analysed using the CLUEreg tool in the GRAND database. This tool employs differential network signatures to identify drugs that potentially target the gene signature associated with the disease. Our analysis utilised a large RNA-seq dataset comprising 532 post-mortem brain samples from the CommonMind project. We constructed co-expression gene regulatory networks for both schizophrenia cases and healthy control subjects, incorporating 15,831 genes and 413 overlapping TFs. Through drug repurposing, we identified 18 promising candidates for repurposing as potential treatments for schizophrenia. The analysis of TF-gene regulatory networks revealed that the TFs in schizophrenia predominantly regulate pathways associated with energy metabolism, immune response, cell adhesion, and thyroid hormone signalling. These pathways represent significant targets for therapeutic intervention. The identified drug repurposing candidates likely act through TF-targeted pathways. These promising candidates, particularly those with preclinical evidence such as rimonabant and kaempferol, warrant further investigation into their potential mechanisms of action and efficacy in alleviating the symptoms of schizophrenia.

SAFER: sub-hypergraph attention-based neural network for predicting effective responses to dose combinations

Background

The potential benefits of drug combination synergy in cancer medicine are significant, yet the risks must be carefully managed due to the possibility of increased toxicity. Although artificial intelligence applications have demonstrated notable success in predicting drug combination synergy, several key challenges persist: (1) Existing models often predict average synergy values across a restricted range of testing dosages, neglecting crucial dose amounts and the mechanisms of action of the drugs involved. (2) Many graph-based models rely on static protein-protein interactions, failing to adapt to dynamic and context-dependent networks. This limitation constrains the applicability of current methods.

Results

We introduced SAFER, a Sub-hypergraph Attention-based graph model, addressing these issues by incorporating complex relationships among biological knowledge networks and considering dosing effects on subject-specific networks. SAFER outperformed previous models on the benchmark and the independent test set. The analysis of subgraph attention weight for the lung cancer cell line highlighted JAK-STAT signaling pathway, PRDM12, ZNF781, and CDC5L that have been implicated in lung fibrosis.

Conclusions

SAFER presents an interpretable framework designed to identify drug-responsive signals. Tailored for comprehending dose effects on subject-specific molecular contexts, our model uniquely captures dose-level drug combination responses. This capability unlocks previously inaccessible avenues of investigation compared to earlier models. Finally, the SAFER framework can be leveraged by future inquiries to investigate molecular networks that uniquely characterize individual patients.

Use of gene regulatory network analysis to repurpose drugs to treat bipolar disorder

BACKGROUND

Bipolar disorder (BD) presents significant challenges in drug discovery, necessitating alternative approaches. Drug repurposing, leveraging computational techniques and expanding biomedical data, holds promise for identifying novel treatment strategies.

METHODS

This study utilized gene regulatory networks (GRNs) to identify significant regulatory changes in BD, using network-based signatures for drug repurposing. Employing the PANDA algorithm, we investigated the variations in transcription factor-GRNs between individuals with BD and unaffected individuals, incorporating binding motifs, protein interactions, and gene co-expression data. The differences in edge weights between BD and controls were then used as differential network signatures to identify drugs potentially targeting the disease-associated gene signature, employing the CLUEreg tool in the GRAND database.

RESULTS

Using a large RNA-seq dataset of 216 post-mortem brain samples from the CommonMind consortium, we constructed GRNs based on co-expression for individuals with BD and unaffected controls, involving 15,271 genes and 405 TFs. Our analysis highlighted significant influences of these TFs on immune response, energy metabolism, cell signalling, and cell adhesion pathways in the disorder. By employing drug repurposing, we identified 10 promising candidates potentially repurposed as BD treatments.

LIMITATIONS

Non-drug-naïve transcriptomics data, bulk analysis of BD samples, potential bias of GRNs towards well-studied genes.

CONCLUSIONS

Further investigation into repurposing candidates, especially those with preclinical evidence supporting their efficacy, like kaempferol and pramocaine, is warranted to understand their mechanisms of action and effectiveness in treating BD. Additionally, novel targets such as PARP1 and A2b offer opportunities for future research on their relevance to the disorder.

Gene Network Analyses Identify Co-regulated Transcription Factors and BACH1 as a Key Driver in Rheumatoid Arthritis Fibroblast-like Synoviocytes

RNA-sequencing and differential gene expression studies have significantly advanced our understanding of pathogenic pathways underlying Rheumatoid Arthritis (RA). Yet, little is known about cell-specific regulatory networks and their contributions to disease. In this study, we focused on fibroblast-like synoviocytes (FLS), a cell type central to disease pathogenesis and joint damage in RA. We used a strategy that computed sample-specific gene regulatory networks (GRNs) to compare network properties between RA and osteoarthritis FLS. We identified 28 transcription factors (TFs) as key regulators central to the signatures of RA FLS. Six of these TFs are new and have not been previously implicated in RA, and included BACH1, HLX, and TGIF1. Several of these TFs were found to be co-regulated, and BACH1 emerged as the most significant TF and regulator. The main BACH1 targets included those implicated in fatty acid metabolism and ferroptosis. The discovery of BACH1 was validated in experiments with RA FLS. Knockdown of BACH1 in RA FLS significantly affected the gene expression signatures, reduced cell adhesion and mobility, interfered with the formation of thick actin fibers, and prevented the polarized formation of lamellipodia, all required for the RA destructive behavior of FLS. This is the first time that BACH1 is shown to have a central role in the regulation of FLS phenotypes, and gene expression signatures, as well as in ferroptosis and fatty acid metabolism. These new discoveries have the potential to become new targets for treatments aimed at selectively targeting the RA FLS.

scGRN: a comprehensive single-cell gene regulatory network platform of human and mouse

Gene regulatory networks (GRNs) are interpretable graph models encompassing the regulatory interactions between transcription factors (TFs) and their downstream target genes. Making sense of the topology and dynamics of GRNs is fundamental to interpreting the mechanisms of disease etiology and translating corresponding findings into novel therapies. Recent advances in single-cell multi-omics techniques have prompted the computational inference of GRNs from single-cell transcriptomic and epigenomic data at an unprecedented resolution. Here, we present scGRN (https://bio.liclab.net/scGRN/), a comprehensive single-cell multi-omics gene regulatory network platform of human and mouse. The current version of scGRN catalogs 237 051 cell type-specific GRNs (62 999 692 TF-target gene pairs), covering 160 tissues/cell lines and 1324 single-cell samples. scGRN is the first resource documenting large-scale cell type-specific GRN information of diverse human and mouse conditions inferred from single-cell multi-omics data. We have implemented multiple online tools for effective GRN analysis, including differential TF-target network analysis, TF enrichment analysis, and pathway downstream analysis. We also provided details about TF binding to promoters, super-enhancers and typical enhancers of target genes in GRNs. Taken together, scGRN is an integrative and useful platform for searching, browsing, analyzing, visualizing and downloading GRNs of interest, enabling insight into the differences in regulatory mechanisms across diverse conditions.

Sex-biased gene expression and gene-regulatory networks of sex-biased adverse event drug targets and drug metabolism genes

BACKGROUND

Previous pharmacovigilance studies and a retroactive review of cancer clinical trial studies identified that women were more likely to experience drug adverse events (i.e., any unintended effects of medication), and men were more likely to experience adverse events that resulted in hospitalization or death. These sex-biased adverse events (SBAEs) are due to many factors not entirely understood, including differences in body mass, hormones, pharmacokinetics, and liver drug metabolism enzymes and transporters.

METHODS

We first identified drugs associated with SBAEs from the FDA Adverse Event Reporting System (FAERS) database. Next, we evaluated sex-specific gene expression of the known drug targets and metabolism enzymes for those SBAE-associated drugs. We also constructed sex-specific tissue gene-regulatory networks to determine if these known drug targets and metabolism enzymes from the SBAE-associated drugs had sex-specific gene-regulatory network properties and predicted regulatory relationships.

RESULTS

We identified liver-specific gene-regulatory differences for drug metabolism genes between males and females, which could explain observed sex differences in pharmacokinetics and pharmacodynamics. In addition, we found that ~ 85% of SBAE-associated drug targets had sex-biased gene expression or were core genes of sex- and tissue-specific network communities, significantly higher than randomly selected drug targets. Lastly, we provide the sex-biased drug-adverse event pairs, drug targets, and drug metabolism enzymes as a resource for the research community.

CONCLUSIONS

Overall, we provide evidence that many SBAEs are associated with drug targets and drug metabolism genes that are differentially expressed and regulated between males and females. These SBAE-associated drug metabolism enzymes and drug targets may be useful for future studies seeking to explain or predict SBAEs.

BIC-LP: A Hybrid Higher-Order Dynamic Bayesian Network Score Function for Gene Regulatory Network Reconstruction

Reconstructing gene regulatory networks(GRNs) is an increasingly hot topic in bioinformatics. Dynamic Bayesian network(DBN) is a stochastic graph model commonly used as a vital model for GRN reconstruction. But probabilistic characteristics of biological networks and the existence of data noise bring great challenges to GRN reconstruction and always lead to many false positive/negative edges. ScoreLasso is a hybrid DBN score function combining DBN and linear regression with good performance. Its performance is, however, limited by first-order assumption and ignorance of the initial network of DBN. In this article, an integrated model based on higher-order DBN model, higher-order Lasso linear regression model and Pearson correlation model is proposed. Based on this, a hybrid higher-order DBN score function for GRN reconstruction is proposed, namely BIC-LP. BIC-LP score function is constructed by adding terms based on Lasso linear regression coefficients and Pearson correlation coefficients on classical BIC score function. Therefore, it could capture more information from dataset and curb information loss, compared with both many existing Bayesian family score functions and many state-of-the-art methods for GRN reconstruction. Experimental results show that BIC-LP can reasonably eliminate some false positive edges while retaining most true positive edges, so as to achieve better GRN reconstruction performance.

Interactome Analysis of Visceral Adipose Tissue Elucidates Gene Regulatory Networks and Novel Gene Candidates in Obesity

Objective

Visceral adiposity is associated with increased proinflammatory activity, insulin resistance, diabetes risk and mortality rate. Numerous individual genes have been associated with obesity, but studies investigating gene-regulatory networks in human visceral obesity are lacking.

Methods

We analyzed gene-regulatory networks in human visceral adipose tissue (VAT) from 48 obese and 11 non-obese Chinese subjects using gene co-expression and network construction with RNA-sequencing data. We also conducted RNA interference-based tests on selected genes for adipocyte differentiation effects.

Results

A scale-free gene co-expression network was constructed from 360 differentially expressed genes between obese and non-obese VAT (absolute log fold-change >1, FDR<0.05) with edge probability >0.8. Gene regulatory network analysis identified candidate transcription factors associated with differentially expressed genes. Fifteen subnetworks (communities) displayed altered connectivity patterns between obese and non-obese networks. Genes in pro-inflammatory pathways showed increased network connectivities in obese VAT whereas the oxidative phosphorylation pathway displayed reduced connections (enrichment FDR<0.05). Functional screening via RNA interference identified SOX30 and OSBPL3 as potential network-derived gene candidates influencing adipocyte differentiation.

Conclusions

This interactome-based approach highlights the network architecture, identifies novel candidate genes, and leads to new hypotheses regarding network-assisted gene regulation in obese vs. non-obese VAT.What is already known about this subject?: Visceral adipose tissue (VAT) is associated with increased levels of proinflammatory activity, insulin resistance, diabetes risk and mortality rate.Gene expression studies have identified candidate genes associated with proinflammatory function in VAT.What are the new findings in your manuscript?: Using integrative network-science, we identified co-expression and gene regulatory networks that are differentially regulated in VAT samples from subjects with and without obesityWe used functional testing (adipocyte differentiation) to validate a subset of novel candidate genes with minimal prior reported associations to obesityHow might your results change the direction of research or the focus of clinical practice: Network biology-based investigation provides a new avenue to our understanding of gene function in visceral adiposityFunctional validation screen allows for the identification of novel gene candidates that may be targeted for the treatment of adipose tissue dysfunction in obesity.

A framework for considering prior information in network-based approaches to omics data analysis

For decades, molecular biologists have been uncovering the mechanics of biological systems. Efforts to bring their findings together have led to the development of multiple databases and information systems that capture and present pathway information in a computable network format. Concurrently, the advent of modern omics technologies has empowered researchers to systematically profile cellular processes across different modalities. Numerous algorithms, methodologies, and tools have been developed to use prior knowledge networks (PKNs) in the analysis of omics datasets. Interestingly, it has been repeatedly demonstrated that the source of prior knowledge can greatly impact the results of a given analysis. For these methods to be successful it is paramount that their selection of PKNs is amenable to the data type and the computational task they aim to accomplish. Here we present a five-level framework that broadly describes network models in terms of their scope, level of detail, and ability to inform causal predictions. To contextualize this framework, we review a handful of network-based omics analysis methods at each level, while also describing the computational tasks they aim to accomplish.

Architecture and topologies of gene regulatory networks associated with breast cancer, adjacent normal, and normal tissues

Most cancer studies employ adjacent normal tissues to tumors (ANTs) as controls, which are not completely normal and represent a pre-cancerous state. However, the regulatory landscape of ANTs compared to tumor and non-tumor-bearing normal tissues is largely unexplored. Among cancers, breast cancer is the most commonly diagnosed cancer and a leading cause of death in women worldwide, with a lack of sufficient treatment regimens for various reasons. Hence, we aimed to gain deeper insights into normal, pre-cancerous, and cancerous regulatory systems of breast tissues towards identifying ANT and subtype-specific candidate genes. For this, we constructed and analyzed eight gene regulatory networks (GRNs), including five subtypes (viz., Basal, Her2, Luminal A, Luminal B, and Normal-Like), one ANT, and two normal tissue networks. Whereas several topological properties of these GRNs enabled us to identify tumor-related features of ANT, escape velocity centrality (EVC+) identified 24 functionally significant common genes, including well-known genes such as E2F1, FOXA1, JUN, BRCA1, GATA3, ERBB2, and ERBB3 across all six tissues including subtypes and ANT. Similarly, the EVC+ also helped us to identify tissue-specific key genes (Basal: 18, Her2: 6, Luminal A: 5, Luminal B: 5, Normal-Like: 2, and ANT: 7). Additionally, differentially correlated switching gene pairs along with functional, pathway, and disease annotations highlighted the cancer-associated role of these genes. In a nutshell, the present study revealed ANT and subtype-specific regulatory features and key candidate genes, which can be explored further using in vitro and in vivo experiments for better and effective disease management at an early stage.

Gene regulatory Networks Reveal Sex Difference in Lung Adenocarcinoma

Lung adenocarcinoma (LUAD) has been observed to have significant sex differences in incidence, prognosis, and response to therapy. However, the molecular mechanisms responsible for these disparities have not been investigated extensively. Sample-specific gene regulatory network methods were used to analyze RNA sequencing data from non-cancerous human lung samples from The Genotype Tissue Expression Project (GTEx) and lung adenocarcinoma primary tumor samples from The Cancer Genome Atlas (TCGA); results were validated on independent data. We observe that genes associated with key biological pathways including cell proliferation, immune response and drug metabolism are differentially regulated between males and females in both healthy lung tissue, as well as in tumor, and that these regulatory differences are further perturbed by tobacco smoking. We also uncovered significant sex bias in transcription factor targeting patterns of clinically actionable oncogenes and tumor suppressor genes, including AKT2 and KRAS. Using differentially regulated genes between healthy and tumor samples in conjunction with a drug repurposing tool, we identified several small-molecule drugs that might have sex-biased efficacy as cancer therapeutics and further validated this observation using an independent cell line database. These findings underscore the importance of including sex as a biological variable and considering gene regulatory processes in developing strategies for disease prevention and management.

Heterogeneity in the gene regulatory landscape of leiomyosarcoma

Characterizing inter-tumor heterogeneity is crucial for selecting suitable cancer therapy, as the presence of diverse molecular subgroups of patients can be associated with disease outcome or response to treatment. While cancer subtypes are often characterized by differences in gene expression, the mechanisms driving these differences are generally unknown. We set out to model the regulatory mechanisms driving sarcoma heterogeneity based on patient-specific, genome-wide gene regulatory networks. We developed a new computational framework, PORCUPINE, which combines knowledge on biological pathways with permutation-based network analysis to identify pathways that exhibit significant regulatory heterogeneity across a patient population. We applied PORCUPINE to patient-specific leiomyosarcoma networks modeled on data from The Cancer Genome Atlas and validated our results in an independent dataset from the German Cancer Research Center. PORCUPINE identified 37 heterogeneously regulated pathways, including pathways representing potential targets for treatment of subgroups of leiomyosarcoma patients, such as FGFR and CTLA4 inhibitory signaling. We validated the detected regulatory heterogeneity through analysis of networks and chromatin states in leiomyosarcoma cell lines. We showed that the heterogeneity identified with PORCUPINE is not associated with methylation profiles or clinical features, thereby suggesting an independent mechanism of patient heterogeneity driven by the complex landscape of gene regulatory interactions.

The Network Zoo: a multilingual package for the inference and analysis of gene regulatory networks

Inference and analysis of gene regulatory networks (GRNs) require software that integrates multi-omic data from various sources. The Network Zoo (netZoo; netzoo.github.io) is a collection of open-source methods to infer GRNs, conduct differential network analyses, estimate community structure, and explore the transitions between biological states. The netZoo builds on our ongoing development of network methods, harmonizing the implementations in various computing languages and between methods to allow better integration of these tools into analytical pipelines. We demonstrate the utility using multi-omic data from the Cancer Cell Line Encyclopedia. We will continue to expand the netZoo to incorporate additional methods.

Mathematical modeling of regulatory networks of intracellular processes - Aims and selected methods

Regulatory networks structure and signaling pathways dynamics are uncovered in time- and resource consuming experimental work. However, it is increasingly supported by modeling, analytical and computational techniques as well as discrete mathematics and artificial intelligence applied to to extract knowledge from existing databases. This review is focused on mathematical modeling used to analyze dynamics and robustness of these networks. This paper presents a review of selected modeling methods that facilitate advances in molecular biology.

Functional impact of multi-omic interactions in breast cancer subtypes

Multi-omic approaches are expected to deliver a broader molecular view of cancer. However, the promised mechanistic explanations have not quite settled yet. Here, we propose a theoretical and computational analysis framework to semi-automatically produce network models of the regulatory constraints influencing a biological function. This way, we identified functions significantly enriched on the analyzed omics and described associated features, for each of the four breast cancer molecular subtypes. For instance, we identified functions sustaining over-representation of invasion-related processes in the basal subtype and DNA modification processes in the normal tissue. We found limited overlap on the omics-associated functions between subtypes; however, a startling feature intersection within subtype functions also emerged. The examples presented highlight new, potentially regulatory features, with sound biological reasons to expect a connection with the functions. Multi-omic regulatory networks thus constitute reliable models of the way omics are connected, demonstrating a capability for systematic generation of mechanistic hypothesis.

Drug screening and biomarker gene investigation in cancer therapy through the human transcriptional regulatory network

A complex and vast biological network regulates all biological functions in the human body in a sophisticated manner, and abnormalities in this network can lead to disease and even cancer. The construction of a high-quality human molecular interaction network is possible with the development of experimental techniques that facilitate the interpretation of the mechanisms of drug treatment for cancer. We collected 11 molecular interaction databases based on experimental sources and constructed a human protein-protein interaction (PPI) network and a human transcriptional regulatory network (HTRN). A random walk-based graph embedding method was used to calculate the diffusion profiles of drugs and cancers, and a pipeline was constructed by using five similarity comparison metrics combined with a rank aggregation algorithm, which can be implemented for drug screening and biomarker gene prediction. Taking NSCLC as an example, curcumin was identified as a potentially promising anticancer drug from 5450 natural small molecules, and combined with differentially expressed genes, survival analysis, and topological ranking, we obtained BIRC5 (survivin), which is both a biomarker for NSCLC and a key target for curcumin. Finally, the binding mode of curcumin and survivin was explored using molecular docking. This work has a guiding significance for antitumor drug screening and the identification of tumor markers.

Proceedings of the fifth international Molecular Pathological Epidemiology (MPE) meeting

Cancer heterogeneities hold the key to a deeper understanding of cancer etiology and progression and the discovery of more precise cancer therapy. Modern pathological and molecular technologies offer a powerful set of tools to profile tumor heterogeneities at multiple levels in large patient populations, from DNA to RNA, protein and epigenetics, and from tumor tissues to tumor microenvironment and liquid biopsy. When coupled with well-validated epidemiologic methodology and well-characterized epidemiologic resources, the rich tumor pathological and molecular tumor information provide new research opportunities at an unprecedented breadth and depth. This is the research space where Molecular Pathological Epidemiology (MPE) emerged over a decade ago and has been thriving since then. As a truly multidisciplinary field, MPE embraces collaborations from diverse fields including epidemiology, pathology, immunology, genetics, biostatistics, bioinformatics, and data science. Since first convened in 2013, the International MPE Meeting series has grown into a dynamic and dedicated platform for experts from these disciplines to communicate novel findings, discuss new research opportunities and challenges, build professional networks, and educate the next-generation scientists. Herein, we share the proceedings of the Fifth International MPE meeting, held virtually online, on May 24 and 25, 2021. The meeting consisted of 21 presentations organized into the three main themes, which were recent integrative MPE studies, novel cancer profiling technologies, and new statistical and data science approaches. Looking forward to the near future, the meeting attendees anticipated continuous expansion and fruition of MPE research in many research fronts, particularly immune-epidemiology, mutational signatures, liquid biopsy, and health disparities.

Repurposing Drugs via Network Analysis: Opportunities for Psychiatric Disorders

Despite advances in pharmacology and neuroscience, the path to new medications for psychiatric disorders largely remains stagnated. Drug repurposing offers a more efficient pathway compared with de novo drug discovery with lower cost and less risk. Various computational approaches have been applied to mine the vast amount of biomedical data generated over recent decades. Among these methods, network-based drug repurposing stands out as a potent tool for the comprehension of multiple domains of knowledge considering the interactions or associations of various factors. Aligned well with the poly-pharmacology paradigm shift in drug discovery, network-based approaches offer great opportunities to discover repurposing candidates for complex psychiatric disorders. In this review, we present the potential of network-based drug repurposing in psychiatry focusing on the incentives for using network-centric repurposing, major network-based repurposing strategies and data resources, applications in psychiatry and challenges of network-based drug repurposing. This review aims to provide readers with an update on network-based drug repurposing in psychiatry. We expect the repurposing approach to become a pivotal tool in the coming years to battle debilitating psychiatric disorders.

Predicting genotype-specific gene regulatory networks

Understanding how each person's unique genotype influences their individual patterns of gene regulation has the potential to improve our understanding of human health and development, and to refine genotype-specific disease risk assessments and treatments. However, the effects of genetic variants are not typically considered when constructing gene regulatory networks, despite the fact that many disease-associated genetic variants are thought to have regulatory effects, including the disruption of transcription factor (TF) binding. We developed EGRET (Estimating the Genetic Regulatory Effect on TFs), which infers a genotype-specific gene regulatory network for each individual in a study population. EGRET begins by constructing a genotype-informed TF-gene prior network derived using TF motif predictions, expression quantitative trait locus (eQTL) data, individual genotypes, and the predicted effects of genetic variants on TF binding. It then uses a technique known as message passing to integrate this prior network with gene expression and TF protein–protein interaction data to produce a refined, genotype-specific regulatory network. We used EGRET to infer gene regulatory networks for two blood-derived cell lines and identified genotype-associated, cell line–specific regulatory differences that we subsequently validated using allele-specific expression, chromatin accessibility QTLs, and differential ChIP-seq TF binding. We also inferred EGRET networks for three cell types from each of 119 individuals and identified cell type–specific regulatory differences associated with diseases related to those cell types. EGRET is, to our knowledge, the first method that infers networks reflective of individual genetic variation in a way that provides insight into the genetic regulatory associations driving complex phenotypes.

gpuZoo: Cost-effective estimation of gene regulatory networks using the Graphics Processing Unit

Gene regulatory network inference allows for the modeling of genome-scale regulatory processes that are altered during development, in disease, and in response to perturbations. Our group has developed a collection of tools to model various regulatory processes, including transcriptional (PANDA, SPIDER) and post-transcriptional (PUMA) gene regulation, as well as gene regulation in individual samples (LIONESS). These methods work by postulating a network structure and then optimizing that structure to be consistent with multiple lines of biological evidence through repeated operations on data matrices. Although our methods are widely used, the corresponding computational complexity, and the associated costs and run times, do limit some applications. To improve the cost/time performance of these algorithms, we developed gpuZoo which implements GPU-accelerated calculations, dramatically improving the performance of these algorithms. The runtime of the gpuZoo implementation in MATLAB and Python is up to 61 times faster and 28 times less expensive than multi-core CPU implementation of the same methods. gpuZoo is available in MATLAB through the netZooM package https://github.com/netZoo/netZooM and in Python through the netZooPy package https://github.com/netZoo/netZooPy.

TissueNexus: a database of human tissue functional gene networks built with a large compendium of curated RNA-seq data

Mapping gene interactions within tissues/cell types plays a crucial role in understanding the genetic basis of human physiology and disease. Tissue functional gene networks (FGNs) are essential models for mapping complex gene interactions. We present TissueNexus, a database of 49 human tissue/cell line FGNs constructed by integrating heterogeneous genomic data. We adopted an advanced machine learning approach for data integration because Bayesian classifiers, which is the main approach used for constructing existing tissue gene networks, cannot capture the interaction and nonlinearity of genomic features well. A total of 1,341 RNA-seq datasets containing 52,087 samples were integrated for all of these networks. Because the tissue label for RNA-seq data may be annotated with different names or be missing, we performed intensive hand-curation to improve quality. We further developed a user-friendly database for network search, visualization, and functional analysis. We illustrate the application of TissueNexus in prioritizing disease genes. The database is publicly available at https://www.diseaselinks.com/TissueNexus/.

Gene regulatory network inference as relaxed graph matching

Bipartite network inference is a ubiquitous problem across disciplines. One important example in the field molecular biology is gene regulatory network inference. Gene regulatory networks are an instrumental tool aiding in the discovery of the molecular mechanisms driving diverse diseases, including cancer. However, only noisy observations of the projections of these regulatory networks are typically assayed. In an effort to better estimate regulatory networks from their noisy projections, we formulate a non-convex but analytically tractable optimization problem called OTTER. This problem can be interpreted as relaxed graph matching between the two projections of the bipartite network. OTTER's solutions can be derived explicitly and inspire a spectral algorithm, for which we provide network recovery guarantees. We also provide an alternative approach based on gradient descent that is more robust to noise compared to the spectral algorithm. Interestingly, this gradient descent approach resembles the message passing equations of an established gene regulatory network inference method, PANDA. Using three cancer-related data sets, we show that OTTER outperforms state-of-the-art inference methods in predicting transcription factor binding to gene regulatory regions. To encourage new graph matching applications to this problem, we have made all networks and validation data publicly available.