Identification of Drought Stress-Responsive Genes in Rice by Random Walk with Multi-Restart Probability on MultiPlex Biological Networks.

Exploring drought stress-responsive genes in rice is essential for breeding drought-resistant varieties. Rice drought resistance is controlled by multiple genes, and mining drought stress-responsive genes solely based on single omics data lacks stability and accuracy. Multi-omics correlation analysis and biological molecular network analysis provide robust solutions. This study proposed a random walk with a multi-restart probability (RWMRP) algorithm, based on the Restarted Random Walk (RWR) algorithm, to operate on rice MultiPlex biological networks. It explores the interactions between biological molecules across various levels and ranks potential genes. RWMRP uses eigenvector centrality to evaluate node importance in the network and adjusts the restart probabilities accordingly, diverging from the uniform restart probability employed in RWR. In the random walk process, it can be better to consider the global relationships in the network. Firstly, we constructed a MultiPlex biological network by integrating the rice protein-protein interaction, gene pathway, and gene co-expression network. Then, we employed RWMRP to predict the potential genes associated with rice tolerance to drought stress. Enrichment and correlation analyses resulted in the identification of 12 drought-related genes. We further conducted quantitative real-time polymerase chain reaction (qRT-PCR) analysis on these 12 genes, ultimately identifying 10 genes responsive to drought stress.

CTGF, FN1, IL-6, THBS1, and WISP1 genes and PI3K-Akt signaling pathway as prognostic and therapeutic targets in gastric cancer identified by gene network modeling.

OBJECTIVE: Gastric cancer (GC) is one of the most common malignancies worldwide and it is considered the fourth most common cause of cancer death. This study aimed to find critical genes/pathways in GC pathogenesis to be used as biomarkers or therapeutic targets. METHODS: Differentially expressed genes were explored between human gastric cancerous and noncancerous tissues, and Gene Ontology and Kyoto Encyclopedia of Genes and Genomes signaling pathway enrichment analyses were done. Hub genes were identified based on the protein-protein interaction network constructed in the STRING database with Cytoscape software. The hub genes were selected for further investigation using GEPIA2 and DrugBank databases. RESULTS: Ten overexpressed hub genes in GC were identified in the current study, including FN1, TP53, IL-6, CXCL5, ELN, ADAMTS2, WISP1, MMP2, CTGF, and THBS1. The study demonstrated the PI3K-Akt pathway's central involvement in GC, with pronounced alterations in essential components. Survival analysis revealed significant correlations between CTGF, FN1, IL-6, THBS1, and WISP1 overexpression and reduced overall survival times in GC patients. CONCLUSION: A mutual interplay emerged, where PI3K-Akt signaling could upregulate certain genes, forming feedback loops and intensifying cancer phenotypes. The interconnected overexpression of genes and the PI3K-Akt pathway fosters gastric tumorigenesis, suggesting therapeutic potential. DrugBank analysis identified limited FDA-approved drugs, advocating for further exploration while targeting these hub genes could reshape GC treatment. The identified genes could be novel diagnostic/prognostic biomarkers or potential therapeutic targets for GC, but further clinical validation is required.

A python library for the fast and scalable computation of biologically meaningful individual specific networks.

Individual Specific Networks (ISNs) are a tool used in computational biology to infer Individual Specific relationships between biological entities from omics data. ISNs provide insights into how the interactions among these entities affect their respective functions. To address the scarcity of solutions for efficiently computing ISNs on large biological datasets, we present ISN-tractor, a data-agnostic, highly optimized Python library to build and analyse ISNs. ISN-tractor demonstrates superior scalability and efficiency in generating Individual Specific Networks (ISNs) when compared to existing methods such as LionessR, both in terms of time and memory usage, allowing ISNs to be used on large datasets. We show how ISN-tractor can be applied to real-life datasets, including The Cancer Genome Atlas (TCGA) and HapMap, showcasing its versatility. ISN-tractor can be used to build ISNs from various -omics data types, including transcriptomics, proteomics, and genotype arrays, and can detect distinct patterns of gene interactions within and across cancer types. We also show how Filtration Curves provided valuable insights into ISN characteristics, revealing topological distinctions among individuals with different clinical outcomes. Additionally, ISN-tractor can effectively cluster populations based on genetic relationships, as demonstrated with Principal Component Analysis on HapMap data.

A technique for preserving network structure in randomized Hi-C data.

Chromatin interaction data are frequently analyzed as a network to study several aspects of chromatin structure. Hi-C experiments are costly and there is a need to create simulated networks for quality assessment or result validation purposes. Existing tools do not maintain network properties during randomization. We propose an algorithm to modify an existing chromatin interaction graph while preserving the graphs most basic topological features - node degrees and interaction length distribution. The algorithm is implemented in Python and its open-source code as well as the data to reproduce the results are available on Github.

BPP: a platform for automatic biochemical pathway prediction.

A biochemical pathway consists of a series of interconnected biochemical reactions to accomplish specific life activities. The participating reactants and resultant products of a pathway, including gene fragments, proteins, and small molecules, coalesce to form a complex reaction network. Biochemical pathways play a critical role in the biochemical domain as they can reveal the flow of biochemical reactions in living organisms, making them essential for understanding life processes. Existing studies of biochemical pathway networks are mainly based on experimentation and pathway database analysis methods, which are plagued by substantial cost constraints. Inspired by the success of representation learning approaches in biomedicine, we develop the biochemical pathway prediction (BPP) platform, which is an automatic BPP platform to identify potential links or attributes within biochemical pathway networks. Our BPP platform incorporates a variety of representation learning models, including the latest hypergraph neural networks technology to model biochemical reactions in pathways. In particular, BPP contains the latest biochemical pathway-based datasets and enables the prediction of potential participants or products of biochemical reactions in biochemical pathways. Additionally, BPP is equipped with an SHAP explainer to explain the predicted results and to calculate the contributions of each participating element. We conduct extensive experiments on our collected biochemical pathway dataset to benchmark the effectiveness of all models available on BPP. Furthermore, our detailed case studies based on the chronological pattern of our dataset demonstrate the effectiveness of our platform. Our BPP web portal, source code and datasets are freely accessible at https://github.com/Glasgow-AI4BioMed/BPP.

Method to link medicines to diseases using multiplex networks.

The reuse of well-established medicines using computational modeling has gained a lot of attention due to its tremendous benefits. Based on this perspective, a new method for linking known medicines to diseases is proposed. The creation of a new treatment or medicine can be financially and temporally costly and the reuse of medicines is one possibility to accelerate this process efficiently. The main purpose of the reuse of medicines is to reduce some stages of the development of new medicines, motivating the proposition of several methods nowadays. In this work, a new method is developed aiming to connect known medicines to diseases based on available networks of protein interactions and available lists of medicines that affect protein action. The concepts of multiplex networks are used to connect subgraphs of vertices that represent medicines and proteins. The core of the procedure is determined by a weighting strategy constructed to define precisely the more relevant connections. The method was compared to other network link methods in the literature and a case study was presented and evaluated by the proposed method.

Genome­wide association study and polygenic risk scores predict psoriasis and its shared phenotypes in Taiwan.

Psoriasis is a chronic inflammatory dermatological disease, and there is a lack of understanding of the genetic factors involved in psoriasis in Taiwan. To establish associations between genetic variations and psoriasis, a genome­wide association study was performed in a cohort of 2,248 individuals with psoriasis and 67,440 individuals without psoriasis. Using the ingenuity pathway analysis software, biological networks were constructed. Human leukocyte antigen (HLA) diplotypes and haplotypes were analyzed using Attribute Bagging (HIBAG)­R software and chi­square analysis. The present study aimed to assess the potential risks associated with psoriasis using a polygenic risk score (PRS) analysis. The genetic association between single nucleotide polymorphisms (SNPs) in psoriasis and various human diseases was assessed by phenome­wide association study. METAL software was used to analyze datasets from China Medical University Hospital (CMUH) and BioBank Japan (BBJ). The results of the present study revealed 8,585 SNPs with a significance threshold of P<5x10­8, located within 153 genes strongly associated with the psoriasis phenotype, particularly on chromosomes 5 and 6. This specific genomic region has been identified by analyzing the biological networks associated with numerous pathways, including immune responses and inflammatory signaling. HLA genotype analysis indicated a strong association between HLA­A*02:07 and HLA­C*06:02 in a Taiwanese population. Based on our PRS analysis, the risk of psoriasis associated with the SNPs identified in the present study was quantified. These SNPs are associated with various dermatological, circulatory, endocrine, metabolic, musculoskeletal, hematopoietic and infectious diseases. The meta­analysis results indicated successful replication of a study conducted on psoriasis in the BBJ. Several genetic loci are significantly associated with susceptibility to psoriasis in Taiwanese individuals. The present study contributes to our understanding of the genetic determinants that play a role in susceptibility to psoriasis. Furthermore, it provides valuable insights into the underlying etiology of psoriasis in the Taiwanese community.

Decoding dynamic miRNA:ceRNA interactions unveils therapeutic insights and targets across predominant cancer landscapes.

Competing endogenous RNAs play key roles in cellular molecular mechanisms through cross-talk in post-transcriptional interactions. Studies on ceRNA cross-talk, which is particularly dependent on the abundance of free transcripts, generally involve large- and small-scale studies involving the integration of transcriptomic data from tissues and correlation analyses. This abundance-dependent nature of ceRNA interactions suggests that tissue- and condition-specific ceRNA dynamics may fluctuate. However, there are no comprehensive studies investigating the ceRNA interactions in normal tissue, ceRNAs that are lost and/or appear in cancerous tissues or their interactions. In this study, we comprehensively analyzed the tumor-specific ceRNA fluctuations observed in the three highest-incidence cancers, LUAD, PRAD, and BRCA, compared to healthy lung, prostate, and breast tissues, respectively. Our observations pertaining to tumor-specific competing endogenous RNA (ceRNA) interactions revealed that, in the cases of lung adenocarcinoma (LUAD), prostate adenocarcinoma (PRAD), and breast invasive carcinoma (BRCA), 3,204, 1,233, and 406 ceRNAs, respectively, engage in post-transcriptional intercommunication within tumor tissues, in contrast to their absence in corresponding healthy samples. We also found that 90 ceRNAs are shared by the three cancer types and that these ceRNAs participate in ceRNA interactions in tumor tissues compared to those in normal tissues. Among the 90 ceRNAs that directly interact with miRNAs, we uncovered a core network of 165 miRNAs and 63 ceRNAs that should be considered in RNA-targeted and RNA-mediated approaches in future studies and could be used in these three aggressive cancer types. More specifically, in this core interaction network, ceRNAs such as GALNT7, KLF9, and DAB2 and miRNAs like miR-106a/b-5p, miR-20a-5p, and miR-519d-3p may have potential as common targets in the three critical cancers. In contrast to conventional methods that construct ceRNA networks using differentially expressed genes compared to normal tissues, our proposed approach identifies ceRNA players by considering their context within the ceRNA:miRNA interactions. Our results have the potential to reveal distinct and common ceRNA interactions in cancer types and to pinpoint critical RNAs, thereby paving the way for RNA-based strategies in the battle against cancer.

Possible New Candidates Involved to Thalidomide-Related Limbs and Cardiac Defects: A Systems Biology Approach.

Thalidomide is a known teratogen that causes malformations especially in heart and limbs. Its mechanism of teratogenicity is still not fully elucidated. Recently, a new target of thalidomide was described, TBX5, and was observed a new interaction between HAND2 and TBX5 that is disrupted in the presence of thalidomide. Therefore, our study aimed to raise potential candidates for thalidomide teratogenesis, through systems biology, evaluating HAND2 and TBX5 interaction and heart and limbs malformations of thalidomide. Genes and proteins related to TBX5 and HAND2 were selected through TF2DNA, REACTOME, Human Phenotype Ontology, and InterPro databases. Networks were assembled using STRING © database. Network analysis were performed in Cytoscape © and R v3.6.2. Differential gene expression (DGE) analysis was performed through gene expression omnibus. We constructed a network for HAND2 and TBX5 interaction; a network for heart and limbs malformations of TE; and the two joined networks. We observed that EP300 protein seemed to be important in all networks. We also looked for proteins containing C2H2 domain in the assembled networks. ZIC3, GLI1, GLI3, ZNF148, and PRDM16 were the ones present in both heart and limbs malformations of TE networks. Furthermore, in the DGE analysis after treatment with thalidomide, we observed that FANCB, ESCO2, and XRCC2 were downregulated and present both in heart and limbs networks. Through systems biology, we were able to point to different new proteins and genes, and selected specially EP300, which was important in all the analyzed networks, to be further evaluated in the TE teratogenicity.

Mechanistic interactions as the origin of modularity in biological networks.

Biological networks are often modular. Explanations for this peculiarity either assume an adaptive advantage of a modular design such as higher robustness, or attribute it to neutral factors such as constraints underlying network assembly. Interestingly, most insights on the origin of modularity stem from models in which interactions are either determined by highly simplistic mechanisms, or have no mechanistic basis at all. Yet, empirical knowledge suggests that biological interactions are often mediated by complex structural or behavioural traits. Here, we investigate the origins of modularity using a model in which interactions are determined by potentially complex traits. Specifically, we model system elements-such as the species in an ecosystem-as finite-state machines (FSMs), and determine their interactions by means of communication between the corresponding FSMs. Using this model, we show that modularity probably emerges for free. We further find that the more modular an interaction network is, the less complex are the traits that mediate the interactions. Altogether, our results suggest that the conditions for modularity to evolve may be much broader than previously thought.

Replication efficiencies of human cytomegalovirus-infected epithelial cells are dependent on source of virus production.

Human cytomegalovirus (HCMV) is a prevalent betaherpesvirus, and infection can lead to a range of symptomatology from mononucleosis to sepsis in immunocompromised individuals. HCMV is also the leading viral cause of congenital birth defects. Lytic replication is supported by many cell types with different kinetics and efficiencies leading to a plethora of pathologies. The goal of these studies was to elucidate HCMV replication efficiencies for viruses produced on different cell types upon infection of epithelial cells by combining experimental approaches with data-driven computational modeling. HCMV was generated from a common genetic background of TB40-BAC4, propagated on fibroblasts (TB40Fb) or epithelial cells (TB40Epi), and used to infect epithelial cells. We quantified cell-associated viral genomes (vDNA), protein levels (pUL44, pp28), and cell-free titers over time for each virus at different multiplicities of infection. We combined experimental quantification with data-driven simulations and determined that parameters describing vDNA synthesis were similar between sources. We found that pUL44 accumulation was higher in TB40Fb than TB40Epi. In contrast, pp28 accumulation was higher in TB40Epi which coincided with a significant increase in titer for TB40Epi over TB40Fb. These differences were most evident during live-cell imaging, which revealed syncytia-like formation during infection by TB40Epi. Simulations of the late lytic replication cycle yielded a larger synthesis constant for pp28 in TB40Epi along with increase in virus output despite similar rates of genome synthesis. By combining experimental and computational modeling approaches, our studies demonstrate that the cellular source of propagated virus impacts viral replication efficiency in target cell types.

Cell reprogramming design by transfer learning of functional transcriptional networks.

Recent developments in synthetic biology, next-generation sequencing, and machine learning provide an unprecedented opportunity to rationally design new disease treatments based on measured responses to gene perturbations and drugs to reprogram cells. The main challenges to seizing this opportunity are the incomplete knowledge of the cellular network and the combinatorial explosion of possible interventions, both of which are insurmountable by experiments. To address these challenges, we develop a transfer learning approach to control cell behavior that is pre-trained on transcriptomic data associated with human cell fates, thereby generating a model of the network dynamics that can be transferred to specific reprogramming goals. The approach combines transcriptional responses to gene perturbations to minimize the difference between a given pair of initial and target transcriptional states. We demonstrate our approach's versatility by applying it to a microarray dataset comprising >9,000 microarrays across 54 cell types and 227 unique perturbations, and an RNASeq dataset consisting of >10,000 sequencing runs across 36 cell types and 138 perturbations. Our approach reproduces known reprogramming protocols with an AUROC of 0.91 while innovating over existing methods by pre-training an adaptable model that can be tailored to specific reprogramming transitions. We show that the number of gene perturbations required to steer from one fate to another increases with decreasing developmental relatedness and that fewer genes are needed to progress along developmental paths than to regress. These findings establish a proof-of-concept for our approach to computationally design control strategies and provide insights into how gene regulatory networks govern phenotype.

Redox organization of living systems.

Redox organization governs an underlying simplicity in living systems. Critically, redox reactions enable the essential characteristics of life: extraction of energy from the environment, use of energy to support metabolic and structural organization, use of dynamic redox responses to defend against environmental threats, and use of redox mechanisms to direct differentiation of cells and organ systems essential for reproduction. These processes are sustained through a redox context in which electron donor/acceptor couples are poised at substantially different steady-state redox potentials, some with relatively reducing steady states and others with relatively oxidizing steady states. Redox-sensitive thiols of the redox proteome, as well as low molecular weight redox-active molecules, are maintained individually by the kinetics of oxidation-reduction within this redox system. Recent research has revealed opposing network interactions of the metallome, redox proteome, metabolome and transcriptome, which appear to be an evolved redox response structure to maintain stability of an organism in the presence of variable oxidative environments. Considerable opportunity exists to improve human health through detailed understanding of these redox networks so that targeted interventions can be developed to support new avenues for redox medicine.

Network integration of thermal proteome profiling with multi-omics data decodes PARP inhibition.

Complex disease phenotypes often span multiple molecular processes. Functional characterization of these processes can shed light on disease mechanisms and drug effects. Thermal Proteome Profiling (TPP) is a mass-spectrometry (MS) based technique assessing changes in thermal protein stability that can serve as proxies of functional protein changes. These unique insights of TPP can complement those obtained by other omics technologies. Here, we show how TPP can be integrated with phosphoproteomics and transcriptomics in a network-based approach using COSMOS, a multi-omics integration framework, to provide an integrated view of transcription factors, kinases and proteins with altered thermal stability. This allowed us to recover consequences of Poly (ADP-ribose) polymerase (PARP) inhibition in ovarian cancer cells on cell cycle and DNA damage response as well as interferon and hippo signaling. We found that TPP offers a complementary perspective to other omics data modalities, and that its integration allowed us to obtain a more complete molecular overview of PARP inhibition. We anticipate that this strategy can be used to integrate functional proteomics with other omics to study molecular processes.

A computational model of circRNA-associated diseases based on a graph neural network: prediction and case studies for follow-up experimental validation.

BACKGROUND: Circular RNAs (circRNAs) have been confirmed to play a vital role in the occurrence and development of diseases. Exploring the relationship between circRNAs and diseases is of far-reaching significance for studying etiopathogenesis and treating diseases. To this end, based on the graph Markov neural network algorithm (GMNN) constructed in our previous work GMNN2CD, we further considered the multisource biological data that affects the association between circRNA and disease and developed an updated web server CircDA and based on the human hepatocellular carcinoma (HCC) tissue data to verify the prediction results of CircDA. RESULTS: CircDA is built on a Tumarkov-based deep learning framework. The algorithm regards biomolecules as nodes and the interactions between molecules as edges, reasonably abstracts multiomics data, and models them as a heterogeneous biomolecular association network, which can reflect the complex relationship between different biomolecules. Case studies using literature data from HCC, cervical, and gastric cancers demonstrate that the CircDA predictor can identify missing associations between known circRNAs and diseases, and using the quantitative real-time PCR (RT-qPCR) experiment of HCC in human tissue samples, it was found that five circRNAs were significantly differentially expressed, which proved that CircDA can predict diseases related to new circRNAs. CONCLUSIONS: This efficient computational prediction and case analysis with sufficient feedback allows us to identify circRNA-associated diseases and disease-associated circRNAs. Our work provides a method to predict circRNA-associated diseases and can provide guidance for the association of diseases with certain circRNAs. For ease of use, an online prediction server ( http://server.malab.cn/CircDA ) is provided, and the code is open-sourced ( https://github.com/nmt315320/CircDA.git ) for the convenience of algorithm improvement.

Computational aspects of two important biochemical networks with respect to some novel molecular descriptors.

Quantitative structure-activity relationship (QSAR) represents quantitative correlation of biological structural features (called as topological indices) and pharmacological activity as response endpoints. Topological index is a molecular descriptor extensively used to study QSAR of pharmaceutical to assess their molecular characteristics by numerical computation. Meanwhile, the topological indices are numerical functions which are used to predict the growth rate of microorganisms in biological networks. Theoretical assessment of microorganism, such as bacteria and viruses help to expedite the vaccine design and discovery process by rationalizing the lead identification, lead optimization and understanding their mechanism of actions. Hypertree, a network structure derived from graph theory, has a great importance in biological networks for growth of microorganisms, such as bacteria and viruses. In this article, some novel eccentric and degree based topological features of two important biological networks (hypertree and its corona product) are obtained on h-level and derived closed formulas for them. Based on the obtained topological features, the biological properties of these networks are investigated.Communicated by Ramaswamy H. Sarma.

Phenotypic plasticity for improved light harvesting, in tandem with methylome repatterning in reef-building corals.

Acclimatization through phenotypic plasticity represents a more rapid response to environmental change than adaptation and is vital to optimize organisms' performance in different conditions. Generally, animals are less phenotypically plastic than plants, but reef-building corals exhibit plant-like properties. They are light dependent with a sessile and modular construction that facilitates rapid morphological changes within their lifetime. We induced phenotypic changes by altering light exposure in a reciprocal transplant experiment and found that coral plasticity is a colony trait emerging from comprehensive morphological and physiological changes within the colony. Plasticity in skeletal features optimized coral light harvesting and utilization and paralleled significant methylome and transcriptome modifications. Network-associated responses resulted in the identification of hub genes and clusters associated to the change in phenotype: inter-partner recognition and phagocytosis, soft tissue growth and biomineralization. Furthermore, we identified hub genes putatively involved in animal photoreception-phototransduction. These findings fundamentally advance our understanding of how reef-building corals repattern the methylome and adjust a phenotype, revealing an important role of light sensing by the coral animal to optimize photosynthetic performance of the symbionts.

A guide to ERK dynamics, part 1: mechanisms and models.

Extracellular signal-regulated kinase (ERK) has long been studied as a key driver of both essential cellular processes and disease. A persistent question has been how this single pathway is able to direct multiple cell behaviors, including growth, proliferation, and death. Modern biosensor studies have revealed that the temporal pattern of ERK activity is highly variable and heterogeneous, and critically, that these dynamic differences modulate cell fate. This two-part review discusses the current understanding of dynamic activity in the ERK pathway, how it regulates cellular decisions, and how these cell fates lead to tissue regulation and pathology. In part 1, we cover the optogenetic and live-cell imaging technologies that first revealed the dynamic nature of ERK, as well as current challenges in biosensor data analysis. We also discuss advances in mathematical models for the mechanisms of ERK dynamics, including receptor-level regulation, negative feedback, cooperativity, and paracrine signaling. While hurdles still remain, it is clear that higher temporal and spatial resolution provide mechanistic insights into pathway circuitry. Exciting new algorithms and advanced computational tools enable quantitative measurements of single-cell ERK activation, which in turn inform better models of pathway behavior. However, the fact that current models still cannot fully recapitulate the diversity of ERK responses calls for a deeper understanding of network structure and signal transduction in general.

Biocaiv: an integrative webserver for motif-based clustering analysis and interactive visualization of biological networks.

BACKGROUND: As an important task in bioinformatics, clustering analysis plays a critical role in understanding the functional mechanisms of many complex biological systems, which can be modeled as biological networks. The purpose of clustering analysis in biological networks is to identify functional modules of interest, but there is a lack of online clustering tools that visualize biological networks and provide in-depth biological analysis for discovered clusters. RESULTS: Here we present BioCAIV, a novel webserver dedicated to maximize its accessibility and applicability on the clustering analysis of biological networks. This, together with its user-friendly interface, assists biological researchers to perform an accurate clustering analysis for biological networks and identify functionally significant modules for further assessment. CONCLUSIONS: BioCAIV is an efficient clustering analysis webserver designed for a variety of biological networks. BioCAIV is freely available without registration requirements at http://bioinformatics.tianshanzw.cn:8888/BioCAIV/ .

PanomiR: a systems biology framework for analysis of multi-pathway targeting by miRNAs.

Charting microRNA (miRNA) regulation across pathways is key to characterizing their function. Yet, no method currently exists that can quantify how miRNAs regulate multiple interconnected pathways or prioritize them for their ability to regulate coordinate transcriptional programs. Existing methods primarily infer one-to-one relationships between miRNAs and pathways using differentially expressed genes. We introduce PanomiR, an in silico framework for studying the interplay of miRNAs and disease functions. PanomiR integrates gene expression, mRNA-miRNA interactions and known biological pathways to reveal coordinated multi-pathway targeting by miRNAs. PanomiR utilizes pathway-activity profiling approaches, a pathway co-expression network and network clustering algorithms to prioritize miRNAs that target broad-scale transcriptional disease phenotypes. It directly resolves differential regulation of pathways, irrespective of their differential gene expression, and captures co-activity to establish functional pathway groupings and the miRNAs that may regulate them. PanomiR uses a systems biology approach to provide broad but precise insights into miRNA-regulated functional programs. It is available at https://bioconductor.org/packages/PanomiR.