An omnidirectional visualization model of personalized gene regulatory networks

High-order interaction modeling of tumor-microenvironment crosstalk for tumor growth

Signaling interactions between cancer cells and nonmalignant cells in the tumor microenvironment (TME) are believed to influence tumor progression and drug resistance. However, the genomic machineries mediating such an influence remain elusive, making it difficult to determine therapeutic targets on the tumor and its microenvironment. Here, we argue that a computational model, derived from the integration of evolutionary game theory and ecosystem theory through allometric scaling law, can chart the genomic atlas of high-order interaction networks involving tumor cells, TME, and tumor mass. We assess the application of this model to identify the causal influence of gene-induced tumor-TME crosstalk on tumor growth. The findings demonstrate that cooperation and competition between tumor cells and their infiltrating microenvironment promote or inhibit tumor growth in diverse ways. We identify specific genes that govern this promotion or inhibition, which can be used as genetic targets to alter tumor growth. This model opens up a new avenue to precisely infer the genomic underpinnings of tumor-TME interactions and their impact on tumor progression from any omics data.

idopNetwork Analysis of Salt-Responsive Transcriptomes Reveals Hub Regulatory Modules and Genes in Populus euphratica

Euphrates poplar (Populus euphratica) is known as a system model to study the genomic mechanisms underlying the salt resistance of woody species. To characterize how dynamic gene regulatory networks (GRNs) drive the defense response of this species to salt stress, we performed mRNA sequencing of P. euphratica roots under short-term (ST) and long-term (LT) salt stress treatments across multiple time points. Comparisons of these transcriptomes revealed the diverged gene expression patterns between the ST and LT treated samples. Based on the informative, dynamic, omnidirectional, and personalized networks model (idopNetwork), inter- and intra-module networks were constructed across different time points for both the ST and LT groups. Through the analysis of the inter-module network, we identified module 4 as the hub, containing the largest number of genes. Further analysis of the gene network within module 4 revealed that gene XM_011048240.1 had the most prominent interactions with other genes. Under short-term salt stress, gene interactions within the network were predominantly promoted, whereas under long-term stress, these interactions shifted towards inhibition. As for the gene ontology (GO) annotation of differentially expressed genes, the results suggest that P. euphratica may employ distinct response mechanisms during the early and late stages of salt stress. Taking together, these results offer valuable insights into the regulatory mechanism involved in P. euphratica's stress response, advancing our understanding of complex biological processes.

Identification of biomarkers and target drugs for melanoma: a topological and deep learning approach

Introduction

Melanoma, a highly aggressive malignancy characterized by rapid metastasis and elevated mortality rates, predominantly originates in cutaneous tissues. While surgical interventions, immunotherapy, and targeted therapies have advanced, the prognosis for advanced-stage melanoma remains dismal. Globally, melanoma incidence continues to rise, with the United States alone reporting over 100,000 new cases and 7,000 deaths annually. Despite the exponential growth of tumor data facilitated by next-generation sequencing (NGS), current analytical approaches predominantly emphasize single-gene analyses, neglecting critical insights into complex gene interaction networks. This study aims to address this gap by systematically exploring immune gene regulatory dynamics in melanoma progression.

Methods

We developed a bidirectional, weighted, signed, and directed topological immune gene regulatory network to compare transcriptional landscapes between benign melanocytic nevi and cutaneous melanoma. Advanced network analysis tools were employed to identify structural disparities and functional module shifts. Key driver genes were validated through topological centrality metrics. Additionally, deep learning models were implemented to predict drug-target interactions, leveraging molecular features derived from network analyses.

Results

Significant topological divergences emerged between nevi and melanoma networks, with dominant functional modules transitioning from cell cycle regulation in benign lesions to DNA repair and cell migration pathways in malignant tumors. A group of genes, including AURKA, CCNE1, APEX2, and EXOC8, were identified as potential orchestrators of immune microenvironment remodeling during malignant transformation. The deep learning framework successfully predicted 23 clinically actionable drug candidates targeting these molecular drivers.

Discussion

The observed module shift from cell cycle to invasion-related pathways provides mechanistic insights into melanoma progression, suggesting early therapeutic targeting of DNA repair machinery might mitigate metastatic potential. The identified hub genes, particularly AURKA and DDX19B, represent novel candidates for immunomodulatory interventions. Our computational drug prediction strategy bridges molecular network analysis with clinical translation, offering a paradigm for precision oncology in melanoma. Future studies should validate these targets in preclinical models and explore network-based biomarkers for early detection.

Cooperation of TRADD- and RIPK1-dependent cell death pathways in maintaining intestinal homeostasis

Dysfunctional NF-κB signaling is critically involved in inflammatory bowel disease (IBD). We investigated the mechanism by which RIPK1 and TRADD, two key mediators of NF-κB signaling, in mediating intestinal pathology using TAK1 IEC deficient model. We show that phosphorylation of TRADD by TAK1 modulates RIPK1-dependent apoptosis. TRADD and RIPK1 act cooperatively to mediate cell death regulated by TNF and TLR signaling. We demonstrate the pathological evolution from RIPK1-dependent ileitis to RIPK1- and TRADD-co-dependent colitis in TAK1 IEC deficient condition. Combined RIPK1 inhibition and TRADD knockout completely protect against intestinal pathology and lethality in TAK1 IEC KO mice. Furthermore, we identify distinctive microbiota dysbiosis biomarkers for RIPK1-dependent ileitis and TRADD-dependent colitis. These findings reveal the cooperation between RIPK1 and TRADD in mediating cell death and inflammation in IBD with NF-κB deficiency and suggest the possibility of combined inhibition of RIPK1 kinase and TRADD as a new therapeutic strategy for IBD.

IdopNetwork as a genomic predictor of drug response

Despite being challenging, elucidating the systematic control mechanisms of multifactorial drug responses is crucial for pharmacogenomic research. We describe a new form of statistical mechanics to reconstruct informative, dynamic, omnidirectional, and personalized networks (idopNetworks), which cover all pharmacogenomic factors and their interconnections, interdependence, and mechanistic roles. IdopNetworks can characterize how cell-cell crosstalk is mediated by genes and proteins to shape body-drug interactions and identify key roadmaps of information flow and propagation for determining drug efficacy and toxicity. We argue that idopNetworks could potentially provide insight into the genomic machinery of drug responses and provide scientific guidance for the design of drugs whose potency is maximized at lower doses.

Network Modeling and Control of Dynamic Disease Pathways, Review and Perspectives

Dynamic disease pathways are a combination of complex dynamical processes among bio-molecules in a cell that leads to diseases. Network modeling of disease pathways considers disease-related bio-molecules (e.g. DNA, RNA, transcription factors, enzymes, proteins, and metabolites) and their interaction (e.g. DNA methylation, histone modification, alternative splicing, and protein modification) to study disease progression and predict therapeutic responses. These bio-molecules and their interactions are the basic elements in the study of the misregulation in the disease-related gene expression that lead to abnormal cellular responses. Gene regulatory networks, cell signaling networks, and metabolic networks are the three major types of intracellular networks for the study of the cellular responses elicited from extracellular signals. The disease-related cellular responses can be prevented or regulated by designing control strategies to manipulate these extracellular or other intracellular signals. The paper reviews the regulatory mechanisms, the dynamic models, and the control strategies for each intracellular network. The applications, limitations and the prospective for modeling and control are also discussed.

Topological change of soil microbiota networks for forest resilience under global warming

Forest management by thinning can mitigate the detrimental impact of increasing drought caused by global warming. Growing evidence shows that the soil microbiota can coordinate the dynamic relationship between forest functions and drought intensity, but how they function as a cohesive whole remains elusive. We outline a statistical topology model to chart the roadmap of how each microbe acts and interacts with every other microbe to shape the dynamic changes of microbial communities under forest management. To demonstrate its utility, we analyze a soil microbiota data collected from a two-way longitudinal factorial experiment involving three stand densities and three levels of rainfall over a growing season in artificial plantations of a forest tree - larix (Larix kaempferi). We reconstruct the most sophisticated soil microbiota networks that code maximally informative microbial interactions and trace their dynamic trajectories across time, space, and environmental signals. By integrating GLMY homology theory, we dissect the topological architecture of these so-called omnidirectional networks and identify key microbial interaction pathways that play a pivotal role in mediating the structure and function of soil microbial communities. The statistical topological model described provides a systems tool for studying how microbial community assembly alters its structure, function and evolution under climate change.

5-methylcytosine RNA modification regulators-based patterns and features of immune microenvironment in acute myeloid leukemia

Acute myeloid leukemia (AML) is a highly heterogeneous malignant disease of the blood cell. The current therapies for AML are unsatisfactory and the molecular mechanisms underlying AML are unclear. 5-methylcytosine (m5C) is an important posttranscriptional modification of mRNA, and is involved in the regulation of mRNA stability, translation, and other aspects of RNA metabolism. However, based on our knowledge of published literature, the role of the m5C regulators has not been explored in AML till date. In this study, we clarified the expression and gene variants of m5C regulators in AML and found that most m5C regulators were differentially expressed and correlated with disease prognosis. We also found that the methylation status of certain m5C regulators (e.g., DNMT3A, DNMT3B) affects the survival of AML patients. Two m5C modification subtypes, and high- and low-risk subgroups identified based on the expression of m5C regulators showed significant differences in the prognosis as well as immune cell infiltration. In addition, most of the m5C regulators were found to be correlated with miRNA expression in AML, as well as IC50 values of many drugs. The miRNA and GSVA analysis were used to identify the different miRNAs and KEGG or hallmark pathways between high- and low-risk subgroups. We also built a prognostic model based on m5C regulators, which was validated by two GSE databases. To verify the reliability of our analysis and conclusions, qPCR was used to identify the expressions of m5C regulators between normal and AML. In summary, we comprehensively explored the molecular characteristics of m5C regulators and built a prognostic model in AML. We proposed new mechanistic insights into the role of m5C in multiple databases and clinical data, which may pave novel ways for the development of therapeutic strategies.

A pleiotropic-epistatic entangelement model of drug response

Because drug response is multifactorial, graph models are uniquely powerful for comprehending its genetic architecture. We deconstruct drug response into many different and interdependent sub-traits, with each sub-trait controlled by multiple genes that act and interact in a complicated manner. The outcome of drug response is the consequence of multileveled intertwined interactions between pleiotropic effects and epistatic effects. Here, we propose a general statistical physics framework to chart the 3D geometric network that codes how epistasis pleiotropically influences a complete set of sub-traits to shape body-drug interactions. This model can dissect the topological architecture of epistatically induced pleiotropic networks (EiPN) and pleiotropically influenced epistatic networks (PiEN). We analyze and interpret the practical implications of the pleiotropic-epistatic entanglement model for pharmacogenomic studies.

The metabolomic physics of complex diseases

Human diseases involve metabolic alterations. Metabolomic profiles have served as a vital biomarker for the early identification of high-risk individuals and disease prevention. However, current approaches can only characterize individual key metabolites, without taking into account the reality that complex diseases are multifactorial, dynamic, heterogeneous, and interdependent. Here, we leverage a statistical physics model to combine all metabolites into bidirectional, signed, and weighted interaction networks and trace how the flow of information from one metabolite to the next causes changes in health state. Viewing a disease outcome as the consequence of complex interactions among its interconnected components (metabolites), we integrate concepts from ecosystem theory and evolutionary game theory to model how the health state-dependent alteration of a metabolite is shaped by its intrinsic properties and through extrinsic influences from its conspecifics. We code intrinsic contributions as nodes and extrinsic contributions as edges into quantitative networks and implement GLMY homology theory to analyze and interpret the topological change of health state from symbiosis to dysbiosis and vice versa. The application of this model to real data allows us to identify several hub metabolites and their interaction webs, which play a part in the formation of inflammatory bowel diseases. The findings by our model could provide important information on drug design to treat these diseases and beyond.

A Cartographic Tool to Predict Disease Risk-associated Pseudo-Dynamic Networks from Tissue-specific Gene Expression

Understanding how genes are differentially expressed across tissues is key to reveal the etiology of human diseases. Genes are never expressed in isolation, but rather co-expressed in a community; thus, they co-act through intricate but well-orchestrated networks. However, existing approaches cannot coalesce the full properties of gene-gene communication and interactions into networks. In particular, the unavailability of dynamic gene expression data might impair the application of existing network models to unleash the complexity of human diseases. To address this limitation, we developed a statistical pipeline named DRDNetPro to visualize and trace how genes dynamically interact with each other across diverse tissues, to ascertain health risk from static expression data. This protocol contains detailed tutorials designed to learn a series of networks, with the illustration example from the Genotype-Tissue Expression (GTEx) project. The proposed toolbox relies on the method developed in our published paper ( Chen et al., 2022 ), coding all genes into bidirectional, signed, weighted, and feedback looped networks, which will provide profound genomic information enabling medical doctors to design precise medicine. Graphical abstract Flowchart illustrating the use of DRDNetPro. The left panel contains the summarized pipeline of DRDNetPro and the right panel contains one pseudo-illustrative example. See the Equipment and Procedure sections for detailed explanations.

Vaginal microbiota networks as a mechanistic predictor of aerobic vaginitis

Aerobic vaginitis (AV) is a complex vaginal dysbiosis that is thought to be caused by the micro-ecological change of the vaginal microbiota. While most studies have focused on how changes in the abundance of individual microbes are associated with the emergence of AV, we still do not have a complete mechanistic atlas of the microbe-AV link. Network modeling is central to understanding the structure and function of any microbial community assembly. By encapsulating the abundance of microbes as nodes and ecological interactions among microbes as edges, microbial networks can reveal how each microbe functions and how one microbe cooperate or compete with other microbes to mediate the dynamics of microbial communities. However, existing approaches can only estimate either the strength of microbe-microbe link or the direction of this link, failing to capture full topological characteristics of a network, especially from high-dimensional microbial data. We combine allometry scaling law and evolutionary game theory to derive a functional graph theory that can characterize bidirectional, signed, and weighted interaction networks from any data domain. We apply our theory to characterize the causal interdependence between microbial interactions and AV. From functional networks arising from different functional modules, we find that, as the only favorable genus from Firmicutes among all identified genera, the role of Lactobacillus in maintaining vaginal microbial symbiosis is enabled by upregulation from other microbes, rather than through any intrinsic capacity. Among Lactobacillus species, the proportion of L. crispatus to L. iners is positively associated with more healthy acid vaginal ecosystems. In a less healthy alkaline ecosystem, L. crispatus establishes a contradictory relationship with other microbes, leading to population decrease relative to L. iners. We identify topological changes of vaginal microbiota networks when the menstrual cycle of women changes from the follicular to luteal phases. Our network tool provides a mechanistic approach to disentangle the internal workings of the microbiota assembly and predict its causal relationships with human diseases including AV.

A statistical framework for recovering pseudo-dynamic networks from static data

MOTIVATION

The collection of temporal or perturbed data is often a prerequisite for reconstructing dynamic networks in most cases. However, these types of data are seldom available for genomic studies in medicine, thus significantly limiting the use of dynamic networks to characterize the biological principles underlying human health and diseases.

RESULTS

We proposed a statistical framework to recover disease risk-associated pseudo-dynamic networks (DRDNet) from steady-state data. We incorporated a varying coefficient model with multiple ordinary differential equations to learn a series of networks. We analyzed the publicly available Genotype-Tissue Expression data to construct networks associated with hypertension risk, and biological findings showed that key genes constituting these networks had pivotal and biologically relevant roles associated with the vascular system. We also provided the selection consistency of the proposed learning procedure and evaluated its utility through extensive simulations.

AVAILABILITY AND IMPLEMENTATION

DRDNet is implemented in the R language, and the source codes are available at https://github.com/chencxxy28/DRDnet/.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

A graph model of combination therapies

The simultaneous use of multiple medications causes drug-drug interactions (DDI) that impact therapeutic efficacy. Here, we argue that graph theory, in conjunction with game theory and ecosystem theory, can address this issue. We treat the coexistence of multiple drugs as a system in which DDI is modeled by game theory. We develop an ordinary differential equation model to characterize how the concentration of a drug changes as a result of its independent capacity and the dependent influence of other drugs through the metabolic response of the host. We coalesce all drugs into personalized and context-specific networks, which can reveal key DDI determinants of therapeutical efficacy. Our model can quantify drug synergy and antagonism and test the translational success of combination therapies to the clinic.

A Single-Cell Omics Network Model of Cell Crosstalk during the Formation of Primordial Follicles

The fate of fetal germ cells (FGCs) in primordial follicles is largely determined by how they interact with the surrounding granulosa cells. However, the molecular mechanisms underlying this interactive process remain poorly understood. Here, we develop a computational model to characterize how individual genes program and rewire cellular crosstalk across FGCs and somas, how gene regulatory networks mediate signaling pathways that functionally link these two cell types, and how different FGCs diversify and evolve through cooperation and competition during embryo development. We analyze single-cell RNA-seq data of human female embryos using the new model, identifying previously uncharacterized mechanisms behind follicle development. The majority of genes (70%) promote FGC–soma synergism, only with a small portion (4%) that incur antagonism; hub genes function reciprocally between the FGC network and soma network; and germ cells tend to cooperate between different stages of development but compete in the same stage within a developmental embryo. Our network model could serve as a powerful tool to unravel the genomic signatures that mediate folliculogenesis from single-cell omics data.

Modeling spatial interaction networks of the gut microbiota

How the gut microbiota is organized across space is postulated to influence microbial succession and its mutualistic relationships with the host. The lack of dynamic or perturbed abundance data poses considerable challenges for characterizing the spatial pattern of microbial interactions. We integrate allometric scaling theory, evolutionary game theory, and prey-predator theory into a unified framework under which quasi-dynamic microbial networks can be inferred from static abundance data. We illustrate that such networks can capture the full properties of microbial interactions, including causality, the sign of the causality, strength, and feedback loop, and are dynamically adaptive along spatial gradients, and context-specific, characterizing variability between individuals and within the same individual across time and space. We design and conduct a gut microbiota study to validate the model, characterizing key spatial determinants of the microbial differences between ulcerative colitis and healthy controls. Our model provides a sophisticated means of unraveling a complete atlas of how microbial interactions vary across space and quantifying causal relationships between such spatial variability and change in health state.

The Genomic Physics of COVID-19 Pathogenesis and Spread

Coronavirus disease (COVID-19) spreads mainly through close contact of infected persons, but the molecular mechanisms underlying its pathogenesis and transmission remain unknown. Here, we propose a statistical physics model to coalesce all molecular entities into a cohesive network in which the roadmap of how each entity mediates the disease can be characterized. We argue that the process of how a transmitter transforms the virus into a recipient constitutes a triad unit that propagates COVID-19 along reticulate paths. Intrinsically, person-to-person transmissibility may be mediated by how genes interact transversely across transmitter, recipient, and viral genomes. We integrate quantitative genetic theory into hypergraph theory to code the main effects of the three genomes as nodes, pairwise cross-genome epistasis as edges, and high-order cross-genome epistasis as hyperedges in a series of mobile hypergraphs. Charting a genome-wide atlas of horizontally epistatic hypergraphs can facilitate the systematic characterization of the community genetic mechanisms underlying COVID-19 spread. This atlas can typically help design effective containment and mitigation strategies and screen and triage those more susceptible persons and those asymptomatic carriers who are incubation virus transmitters.

Genetic dissection of growth trajectories in forest trees: From FunMap to FunGraph

Growth is the developmental process involving important genetic components. Functional mapping (FunMap) has been used as an approach to map quantitative trait loci (QTLs) governing growth trajectories by incorporating growth equations. FunMap is based on reductionism thinking, with a power to identify a small set of significant QTLs from the whole pool of genome-wide markers. Yet, increasing evidence shows that a complex trait is controlled by all genes the organism may possibly carry. Here, we describe and demonstrate a different mapping approach that encapsulates all markers into genetic interaction networks. This approach, symbolized as FunGraph, combines functional mapping, evolutionary game theory, and prey-predator theory into mathematical graphs, allowing the observed genetic effect of a locus to be decomposed into its independent component (resulting from this locus' intrinsic capacity) and dependent component (due to extrinsic regulation by other loci). Using FunGraph, we can visualize and trace the roadmap of how each locus interact with every other locus to impact growth. In a population-based association study of Euphrates poplar, we use FunGraph to identify the previously neglected genetic interaction effects that contribute to the genetic architecture of juvenile stem growth. FunGraph could open up a novel gateway to comprehend the global genetic control mechanisms of complex traits.

Modeling genome-wide by environment interactions through omnigenic interactome networks

How genes interact with the environment to shape phenotypic variation and evolution is a fundamental question intriguing to biologists from various fields. Existing linear models built on single genes are inadequate to reveal the complexity of genotype-environment (G-E) interactions. Here, we develop a conceptual model for mechanistically dissecting G-E interplay by integrating previously disconnected theories and methods. Under this integration, evolutionary game theory, developmental modularity theory, and a variable selection method allow us to reconstruct environment-induced, maximally informative, sparse, and casual multilayer genetic networks. We design and conduct two mapping experiments by using a desert-adapted tree species to validate the biological application of the model proposed. The model identifies previously uncharacterized molecular mechanisms that mediate trees' response to saline stress. Our model provides a tool to comprehend the genetic architecture of trait variation and evolution and trace the information flow of each gene toward phenotypes within omnigenic networks.

Computational Identification of Gene Networks as a Biomarker of Neuroblastoma Risk

Neuroblastoma is a common cancer in children, affected by a number of genes that interact with each other through intricate but coordinated networks. Traditional approaches can only reconstruct a single regulatory network that is topologically not informative enough to explain the complexity of neuroblastoma risk. We implemented and modified an advanced model for recovering informative, omnidirectional, dynamic, and personalized networks (idopNetworks) from static gene expression data for neuroblastoma risk. We analyzed 3439 immune genes of neuroblastoma for 217 high-risk patients and 30 low-risk patients by which to reconstruct large patient-specific idopNetworks. By converting these networks into risk-specific representations, we found that the shift in patients from a low to high risk or from a high to low risk might be due to the reciprocal change of hub regulators. By altering the directions of regulation exerted by these hubs, it may be possible to reduce a high risk to a low risk. Results from a holistic, systems-oriented paradigm through idopNetworks can potentially enable oncologists to experimentally identify the biomarkers of neuroblastoma and other cancers.