Gene regulatory networks and their applications: understanding biological and medical problems in terms of networks

Augusta: From RNA-Seq to gene regulatory networks and Boolean models

Computational models of gene regulations help to understand regulatory mechanisms and are extensively used in a wide range of areas, e.g., biotechnology or medicine, with significant benefits. Unfortunately, there are only a few computational gene regulatory models of whole genomes allowing static and dynamic analysis due to the lack of sophisticated tools for their reconstruction. Here, we describe Augusta, an open-source Python package for Gene Regulatory Network (GRN) and Boolean Network (BN) inference from the high-throughput gene expression data. Augusta can reconstruct genome-wide models suitable for static and dynamic analyses. Augusta uses a unique approach where the first estimation of a GRN inferred from expression data is further refined by predicting transcription factor binding motifs in promoters of regulated genes and by incorporating verified interactions obtained from databases. Moreover, a refined GRN is transformed into a draft BN by searching in the curated model database and setting logical rules to incoming edges of target genes, which can be further manually edited as the model is provided in the SBML file format. The approach is applicable even if information about the organism under study is not available in the databases, which is typically the case for non-model organisms including most microbes. Augusta can be operated from the command line and, thus, is easy to use for automated prediction of models for various genomes. The Augusta package is freely available at github.com/JanaMus/Augusta. Documentation and tutorials are available at augusta.readthedocs.io.

Investigating the Impact of Flavonoids on Aspergillus flavus: Insights into Cell Wall Damage and Biofilms

Aspergillus flavus, a fungus known for producing aflatoxins, poses significant threats to agriculture and global health. Flavonoids, plant-derived compounds, inhibit A. flavus proliferation and mitigate aflatoxin production, although the precise molecular and physical mechanisms underlying these effects remain poorly understood. In this study, we investigated three flavonoids-apigenin, luteolin, and quercetin-applied to A. flavus NRRL 3357. We determined the following: (1) glycosylated luteolin led to a 10% reduction in maximum fungal growth capacity; (2) quercetin affected cell wall integrity by triggering extreme mycelial collapse, while apigenin and luteolin caused peeling of the outer layer of cell wall; (3) luteolin exhibited the highest antioxidant capacity in the environment compared to apigenin and quercetin; (4) osmotic stress assays did not reveal morphological defects; (5) flavonoids promoted cell adherence, a precursor for biofilm formation; and (6) RNA sequencing analysis revealed that flavonoids impact expression of putative cell wall and plasma membrane biosynthesis genes. Our findings suggest that the differential effects of quercetin, luteolin, and apigenin on membrane integrity and biofilm formation may be driven by their interactions with fungal cell walls. These insights may inform the development of novel antifungal additives or plant breeding strategies focusing on plant-derived compounds in crop protection.

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.

Biological computations: Limitations of attractor-based formalisms and the need for transients

Living systems, from single cells to higher vertebrates, receive a continuous stream of non-stationary inputs that they sense, for e.g. via cell surface receptors or sensory organs. By integrating these time-varying, multi-sensory, and often noisy information with memory using complex molecular or neuronal networks, they generate a variety of responses beyond simple stimulus-response association, including avoidance behavior, life-long-learning or social interactions. In a broad sense, these processes can be understood as a type of biological computation. Taking as a basis generic features of biological computations, such as real-time responsiveness or robustness and flexibility of the computation, we highlight the limitations of the current attractor-based framework for understanding computations in biological systems. We argue that frameworks based on transient dynamics away from attractors are better suited for the description of computations performed by neuronal and signaling networks. In particular, we discuss how quasi-stable transient dynamics from ghost states that emerge at criticality have a promising potential for developing an integrated framework of computations, that can help us understand how living system actively process information and learn from their continuously changing environment.

Holimap: an accurate and efficient method for solving stochastic gene network dynamics

Gene-gene interactions are crucial to the control of sub-cellular processes but our understanding of their stochastic dynamics is hindered by the lack of simulation methods that can accurately and efficiently predict how the distributions of gene product numbers vary across parameter space. To overcome these difficulties, here we present Holimap (high-order linear-mapping approximation), an approach that approximates the protein or mRNA number distributions of a complex gene regulatory network by the distributions of a much simpler reaction system. We demonstrate Holimap's computational advantages over conventional methods by applying it to predict the stochastic time-dependent dynamics of various gene networks, including transcriptional networks ranging from simple autoregulatory loops to complex randomly connected networks, post-transcriptional networks, and post-translational networks. Holimap is ideally suited to study how the intricate network of gene-gene interactions results in precise coordination and control of gene expression.

Detecting responsible nodes in differential Bayesian networks

To study the roles that different nodes play in differentiating Bayesian networks under two states, such as control versus disease, we formulate two node-specific scores to facilitate such assessment. The first score is motivated by the prediction invariance property of a causal model. The second score results from modifying an existing score constructed for differential analysis of undirected networks. We develop strategies based on these scores to identify nodes responsible for topological differences between two Bayesian networks. Synthetic data and real-life data from designed experiments are used to demonstrate the efficacy of the proposed methods in detecting responsible nodes.

DeepIMAGER: Deeply Analyzing Gene Regulatory Networks from scRNA-seq Data

Understanding the dynamics of gene regulatory networks (GRNs) across diverse cell types poses a challenge yet holds immense value in unraveling the molecular mechanisms governing cellular processes. Current computational methods, which rely solely on expression changes from bulk RNA-seq and/or scRNA-seq data, often result in high rates of false positives and low precision. Here, we introduce an advanced computational tool, DeepIMAGER, for inferring cell-specific GRNs through deep learning and data integration. DeepIMAGER employs a supervised approach that transforms the co-expression patterns of gene pairs into image-like representations and leverages transcription factor (TF) binding information for model training. It is trained using comprehensive datasets that encompass scRNA-seq profiles and ChIP-seq data, capturing TF-gene pair information across various cell types. Comprehensive validations on six cell lines show DeepIMAGER exhibits superior performance in ten popular GRN inference tools and has remarkable robustness against dropout-zero events. DeepIMAGER was applied to scRNA-seq datasets of multiple myeloma (MM) and detected potential GRNs for TFs of RORC, MITF, and FOXD2 in MM dendritic cells. This technical innovation, combined with its capability to accurately decode GRNs from scRNA-seq, establishes DeepIMAGER as a valuable tool for unraveling complex regulatory networks in various cell types.

Developing a novel causal inference algorithm for personalized biomedical causal graph learning using meta machine learning

BACKGROUND

Modeling causality through graphs, referred to as causal graph learning, offers an appropriate description of the dynamics of causality. The majority of current machine learning models in clinical decision support systems only predict associations between variables, whereas causal graph learning models causality dynamics through graphs. However, building personalized causal graphs for each individual is challenging due to the limited amount of data available for each patient.

METHOD

In this study, we present a new algorithmic framework using meta-learning for learning personalized causal graphs in biomedicine. Our framework extracts common patterns from multiple patient graphs and applies this information to develop individualized graphs. In multi-task causal graph learning, the proposed optimized initial guess of shared commonality enables the rapid adoption of knowledge to new tasks for efficient causal graph learning.

RESULTS

Experiments on one real-world biomedical causal graph learning benchmark data and four synthetic benchmarks show that our algorithm outperformed the baseline methods. Our algorithm can better understand the underlying patterns in the data, leading to more accurate predictions of the causal graph. Specifically, we reduce the structural hamming distance by 50-75%, indicating an improvement in graph prediction accuracy. Additionally, the false discovery rate is decreased by 20-30%, demonstrating that our algorithm made fewer incorrect predictions compared to the baseline algorithms.

CONCLUSION

To the best of our knowledge, this is the first study to demonstrate the effectiveness of meta-learning in personalized causal graph learning and cause inference modeling for biomedicine. In addition, the proposed algorithm can also be generalized to transnational research areas where integrated analysis is necessary for various distributions of datasets, including different clinical institutions.

GRouNdGAN: GRN-guided simulation of single-cell RNA-seq data using causal generative adversarial networks

We introduce GRouNdGAN, a gene regulatory network (GRN)-guided reference-based causal implicit generative model for simulating single-cell RNA-seq data, in silico perturbation experiments, and benchmarking GRN inference methods. Through the imposition of a user-defined GRN in its architecture, GRouNdGAN simulates steady-state and transient-state single-cell datasets where genes are causally expressed under the control of their regulating transcription factors (TFs). Training on six experimental reference datasets, we show that our model captures non-linear TF-gene dependencies and preserves gene identities, cell trajectories, pseudo-time ordering, and technical and biological noise, with no user manipulation and only implicit parameterization. GRouNdGAN can synthesize cells under new conditions to perform in silico TF knockout experiments. Benchmarking various GRN inference algorithms reveals that GRouNdGAN effectively bridges the existing gap between simulated and biological data benchmarks of GRN inference algorithms, providing gold standard ground truth GRNs and realistic cells corresponding to the biological system of interest.

Investigating the mechanisms underlying resistance to chemoterapy and to CRISPR-Cas9 in cancer cell lines

Cancer is one of the major causes of death worldwide and the development of multidrug resistance (MDR) in cancer cells is the principal cause of chemotherapy failure. To gain insights into the specific mechanisms of MDR in cancer cell lines, we developed a novel method for the combined analysis of recently published datasets on drug sensitivity and CRISPR loss-of-function screens for the same set of cancer cell lines. For our analysis, we first selected cell lines that consistently exhibit drug resistance across several classes of compounds. We then identified putative resistance genes for each class of compound and used inferred gene regulatory networks (GRNs) to study possible mechanisms underlying the development of MDR in the identified cancer cell lines. We show that the same method of analysis can also be used to identify cell lines that consistently exhibit resistance to the gene knockout effect of the CRISPR-Cas9 technique and to study the possible underlying mechanisms. In the GRN associated to the drug resistant cell lines, we identify genes previously associated with resistance (UHMK1, RALYL, MGST3, USP9X, and ESRG), genes for which an indirect association can be identified (SPINK13, LINC00664, MRPL38, and EMILIN3), and genes that are found to be overexpressed in non-resistant cancer cell lines (MRPL38, EMILIN3 and RALYL). In the GRNs associated to the CRISPR-Cas9 resistance mechanism, none of the identified genes has been previously reported in the admittedly sparse literature on the subject. However, some of these genes have a common role: APBB2, RUNX1T1, ZBTB7C, and ISX regulate transcription, while APBB2, BTG3, ZBTB7C, SZRD1 and LEF1 have a function in regulating proliferation, suggesting a role for these two pathways. While our results are specific for the lung cancer cell lines we selected for this work, our method of analysis can be applied to cell lines from other tissues and for which the required data is available.

From local to global gene co-expression estimation using single-cell RNA-seq data

In genomics studies, the investigation of gene relationships often brings important biological insights. Currently, the large heterogeneous datasets impose new challenges for statisticians because gene relationships are often local. They change from one sample point to another, may only exist in a subset of the sample, and can be nonlinear or even nonmonotone. Most previous dependence measures do not specifically target local dependence relationships, and the ones that do are computationally costly. In this paper, we explore a state-of-the-art network estimation technique that characterizes gene relationships at the single cell level, under the name of cell-specific gene networks. We first show that averaging the cell-specific gene relationship over a population gives a novel univariate dependence measure, the averaged Local Density Gap (aLDG), that accumulates local dependence and can detect any nonlinear, nonmonotone relationship. Together with a consistent nonparametric estimator, we establish its robustness on both the population and empirical levels. Then, we show that averaging the cell-specific gene relationship over mini-batches determined by some external structure information (eg, spatial or temporal factor) better highlights meaningful local structure change points. We explore the application of aLDG and its minibatch variant in many scenarios, including pairwise gene relationship estimation, bifurcating point detection in cell trajectory, and spatial transcriptomics structure visualization. Both simulations and real data analysis show that aLDG outperforms existing ones.

Functional selectivity of Receptor Tyrosine Kinases regulates distinct cellular outputs

Functional selectivity refers to the activation of differential signalling and cellular outputs downstream of the same membrane-bound receptor when activated by two or more different ligands. Functional selectivity has been described and extensively studied for G-protein Coupled Receptors (GPCRs), leading to specific therapeutic options for dysregulated GPCRs functions. However, studies regarding the functional selectivity of Receptor Tyrosine Kinases (RTKs) remain sparse. Here, we will summarize recent data about RTK functional selectivity focusing on how the nature and the amount of RTK ligands and the crosstalk of RTKs with other membrane proteins regulate the specificity of RTK signalling. In addition, we will discuss how structural changes in RTKs upon ligand binding affects selective signalling pathways. Much remains to be known about the integration of different signals affecting RTK signalling specificity to orchestrate long-term cellular outcomes. Recent advancements in omics, specifically quantitative phosphoproteomics, and in systems biology methods to study, model and integrate different types of large-scale omics data have increased our ability to compare several signals affecting RTK functional selectivity in a global, system-wide fashion. We will discuss how such methods facilitate the exploration of important signalling hubs and enable data-driven predictions aiming at improving the efficacy of therapeutics for diseases like cancer, where redundant RTK signalling pathways often compromise treatment efficacy.

Sex diversity in the 21st century: Concepts, frameworks, and approaches for the future of neuroendocrinology

Sex is ubiquitous and variable throughout the animal kingdom. Historically, scientists have used reductionist methodologies that rely on a priori sex categorizations, in which two discrete sexes are inextricably linked with gamete type. However, this binarized operationalization does not adequately reflect the diversity of sex observed in nature. This is due, in part, to the fact that sex exists across many levels of biological analysis, including genetic, molecular, cellular, morphological, behavioral, and population levels. Furthermore, the biological mechanisms governing sex are embedded in complex networks that dynamically interact with other systems. To produce the most accurate and scientifically rigorous work examining sex in neuroendocrinology and to capture the full range of sex variability and diversity present in animal systems, we must critically assess the frameworks, experimental designs, and analytical methods used in our research. In this perspective piece, we first propose a new conceptual framework to guide the integrative study of sex. Then, we provide practical guidance on research approaches for studying sex-associated variables, including factors to consider in study design, selection of model organisms, experimental methodologies, and statistical analyses. We invite fellow scientists to conscientiously apply these modernized approaches to advance our biological understanding of sex and to encourage academically and socially responsible outcomes of our work. By expanding our conceptual frameworks and methodological approaches to the study of sex, we will gain insight into the unique ways that sex exists across levels of biological organization to produce the vast array of variability and diversity observed in nature.

Prospects and perspectives: inferring physiological and regulatory targets for CAM from molecular and modelling approaches

BACKGROUND AND SCOPE

This review summarizes recent advances in our understanding of Crassulacean Acid Metabolism (CAM) by integrating evolutionary, ecological, physiological, metabolic and molecular perspectives. A number of key control loops which moderate the expression of CAM phases, and their metabolic and molecular control, are explored. These include nocturnal stomatal opening, activation of phosphoenolpyruvate carboxylase by a specific protein kinase, interactions with circadian clock control, as well as daytime decarboxylation and activation of Rubisco. The vacuolar storage and release of malic acid and the interplay between the supply and demand for carbohydrate reserves are also key metabolic control points.

FUTURE OPPORTUNITIES

We identify open questions and opportunities, with experimentation informed by top-down molecular modelling approaches allied with bottom-up mechanistic modelling systems. For example, mining transcriptomic datasets using high-speed systems approaches will help to identify targets for future genetic manipulation experiments to define the regulation of CAM (whether circadian or metabolic control). We emphasize that inferences arising from computational approaches or advanced nuclear sequencing techniques can identify potential genes and transcription factors as regulatory targets. However, these outputs then require systematic evaluation, using genetic manipulation in key model organisms over a developmental progression, combining gene silencing and metabolic flux analysis and modelling to define functionality across the CAM day-night cycle. From an evolutionary perspective, the origins and function of CAM succulents and responses to water deficits are set against the mesophyll and hydraulic limitations imposed by cell and tissue succulence in contrasting morphological lineages. We highlight the interplay between traits across shoots (3D vein density, mesophyll conductance and cell shrinkage) and roots (xylem embolism and segmentation). Thus, molecular, biophysical and biochemical processes help to curtail water losses and exploit rapid rehydration during restorative rain events. In the face of a changing climate, we hope such approaches will stimulate opportunities for future research.

Development of the RF-GSEA Method for Identifying Disulfidptosis-Related Genes and Application in Hepatocellular Carcinoma

Disulfidptosis is a newly discovered cellular programmed cell death mode. Presently, a considerable number of genes related to disulfidptosis remain undiscovered, and its significance in hepatocellular carcinoma remains unrevealed. We have developed a powerful analytical method called RF-GSEA for identifying potential genes associated with disulfidptosis. This method draws inspiration from gene regulation networks and graph theory, and it is implemented through a combination of random forest regression model and Gene Set Enrichment Analysis. Subsequently, to validate the practical application value of this method, we applied it to hepatocellular carcinoma. Based on the RF-GSEA method, we developed a disulfidptosis-related signature. Lastly, we looked into how the disulfidptosis-related signature is connected to HCC prognosis, the tumor microenvironment, the effectiveness of immunotherapy, and the sensitivity of chemotherapy drugs. The RF-GSEA method identified a total of 220 disulfidptosis-related genes, from which 7 were selected to construct the disulfidptosis-related signature. The high-disulfidptosis-related score group had a worse prognosis compared to the low-disulfidptosis-related score group and showed lower infiltration levels of immune-promoting cells. The high-disulfidptosis-related score group had a higher likelihood of benefiting from immunotherapy compared to the low-disulfidptosis-related score group. The RF-GSEA method is a powerful tool for identifying disulfidptosis-related genes. The disulfidptosis-related signature effectively predicts HCC prognosis, immunotherapy response, and drug sensitivity.

SIGNET: transcriptome-wide causal inference for gene regulatory networks

Gene regulation plays an important role in understanding the mechanisms of human biology and diseases. However, inferring causal relationships between all genes is challenging due to the large number of genes in the transcriptome. Here, we present SIGNET (Statistical Inference on Gene Regulatory Networks), a flexible software package that reveals networks of causal regulation between genes built upon large-scale transcriptomic and genotypic data at the population level. Like Mendelian randomization, SIGNET uses genotypic variants as natural instrumental variables to establish such causal relationships but constructs a transcriptome-wide gene regulatory network with high confidence. SIGNET makes such a computationally heavy task feasible by deploying a well-designed statistical algorithm over a parallel computing environment. It also provides a user-friendly interface allowing for parameter tuning, efficient parallel computing scheduling, interactive network visualization, and confirmatory results retrieval. The Open source SIGNET software is freely available ( https://www.zstats.org/signet/ ).

A Differential Transcriptional Regulome Approach to Unpack Cancer Biology: Insights on Renal Cell Carcinoma Subtypes

Cancer research calls for new approaches that account for the regulatory complexities of biology. We present, in this study, the differential transcriptional regulome (DIFFREG) approach for the identification and prioritization of key transcriptional regulators and apply it to the case of renal cell carcinoma (RCC) biology. Of note, RCC has a poor prognosis and the biomarker and drug discovery studies to date have tended to focus on gene expression independent from mutations and/or post-translational modifications. DIFFREG focuses on the differential regulation between transcription factors (TFs) and their target genes rather than differential gene expression and integrates transcriptome profiling with the human transcriptional regulatory network to analyze differential gene regulation between healthy and RCC cases. In this study, RNA-seq tissue samples (n = 1020) from the Cancer Genome Atlas (TCGA), including healthy and tumor subjects, were integrated with a comprehensive human TF-gene interactome dataset (1122603 interactions between 1289 TFs and 25177 genes). Comparative analysis of DIFFREG profiles, consisting of perturbed TF-gene interactions, from three common subtypes (clear cell RCC, papillary RCC and chromophobe RCC) revealed subtype-specific alterations, supporting the hypothesis that these signatures in the transcriptional regulome profiles may be considered potential biomarkers that may play an important role in elucidating the molecular mechanisms of RCC development and translating knowledge about the genetic basis of RCC into the clinic. In addition, these indicators may help oncologists make the best decisions for diagnosis and prognosis management.

Inference of differential key regulatory networks and mechanistic drug repurposing candidates from scRNA-seq data with SCANet

MOTIVATION

The reconstruction of small key regulatory networks that explain the differences in the development of cell (sub)types from single-cell RNA sequencing is a yet unresolved computational problem.

RESULTS

To this end, we have developed SCANet, an all-in-one package for single-cell profiling that covers the whole differential mechanotyping workflow, from inference of trait/cell-type-specific gene co-expression modules, driver gene detection, and transcriptional gene regulatory network reconstruction to mechanistic drug repurposing candidate prediction. To illustrate the power of SCANet, we examined data from two studies. First, we identify the drivers of the mechanotype of a cytokine storm associated with increased mortality in patients with acute respiratory illness. Secondly, we find 20 drugs for eight potential pharmacological targets in cellular driver mechanisms in the intestinal stem cells of obese mice.

AVAILABILITY AND IMPLEMENTATION

SCANet is a free, open-source, and user-friendly Python package that can be seamlessly integrated into single-cell-based systems medicine research and mechanistic drug discovery.

iLSGRN: inference of large-scale gene regulatory networks based on multi-model fusion

MOTIVATION

Gene regulatory networks (GRNs) are a way of describing the interaction between genes, which contribute to revealing the different biological mechanisms in the cell. Reconstructing GRNs based on gene expression data has been a central computational problem in systems biology. However, due to the high dimensionality and non-linearity of large-scale GRNs, accurately and efficiently inferring GRNs is still a challenging task.

RESULTS

In this article, we propose a new approach, iLSGRN, to reconstruct large-scale GRNs from steady-state and time-series gene expression data based on non-linear ordinary differential equations. Firstly, the regulatory gene recognition algorithm calculates the Maximal Information Coefficient between genes and excludes redundant regulatory relationships to achieve dimensionality reduction. Then, the feature fusion algorithm constructs a model leveraging the feature importance derived from XGBoost (eXtreme Gradient Boosting) and RF (Random Forest) models, which can effectively train the non-linear ordinary differential equations model of GRNs and improve the accuracy and stability of the inference algorithm. The extensive experiments on different scale datasets show that our method makes sensible improvement compared with the state-of-the-art methods. Furthermore, we perform cross-validation experiments on the real gene datasets to validate the robustness and effectiveness of the proposed method.

AVAILABILITY AND IMPLEMENTATION

The proposed method is written in the Python language, and is available at: https://github.com/lab319/iLSGRN.

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.

Unraveling patient heterogeneity in complex diseases through individualized co-expression networks: a perspective

This perspective highlights the potential of individualized networks as a novel strategy for studying complex diseases through patient stratification, enabling advancements in precision medicine. We emphasize the impact of interpatient heterogeneity resulting from genetic and environmental factors and discuss how individualized networks improve our ability to develop treatments and enhance diagnostics. Integrating system biology, combining multimodal information such as genomic and clinical data has reached a tipping point, allowing the inference of biological networks at a single-individual resolution. This approach generates a specific biological network per sample, representing the individual from which the sample originated. The availability of individualized networks enables applications in personalized medicine, such as identifying malfunctions and selecting tailored treatments. In essence, reliable, individualized networks can expedite research progress in understanding drug response variability by modeling heterogeneity among individuals and enabling the personalized selection of pharmacological targets for treatment. Therefore, developing diverse and cost-effective approaches for generating these networks is crucial for widespread application in clinical services.

Transcriptome profiling of the Caenorhabditis elegans intestine reveals that ELT-2 negatively and positively regulates intestinal gene expression within the context of a gene regulatory network

ELT-2 is the major transcription factor (TF) required for Caenorhabditis elegans intestinal development. ELT-2 expression initiates in embryos to promote development and then persists after hatching through the larval and adult stages. Though the sites of ELT-2 binding are characterized and the transcriptional changes that result from ELT-2 depletion are known, an intestine-specific transcriptome profile spanning developmental time has been missing. We generated this dataset by performing Fluorescence Activated Cell Sorting on intestine cells at distinct developmental stages. We analyzed this dataset in conjunction with previously conducted ELT-2 studies to evaluate the role of ELT-2 in directing the intestinal gene regulatory network through development. We found that only 33% of intestine-enriched genes in the embryo were direct targets of ELT-2 but that number increased to 75% by the L3 stage. This suggests additional TFs promote intestinal transcription especially in the embryo. Furthermore, only half of ELT-2's direct target genes were dependent on ELT-2 for their proper expression levels, and an equal proportion of those responded to elt-2 depletion with over-expression as with under-expression. That is, ELT-2 can either activate or repress direct target genes. Additionally, we observed that ELT-2 repressed its own promoter, implicating new models for its autoregulation. Together, our results illustrate that ELT-2 impacts roughly 20-50% of intestine-specific genes, that ELT-2 both positively and negatively controls its direct targets, and that the current model of the intestinal regulatory network is incomplete as the factors responsible for directing the expression of many intestinal genes remain unknown.

From time-series transcriptomics to gene regulatory networks: A review on inference methods

Inference of gene regulatory networks has been an active area of research for around 20 years, leading to the development of sophisticated inference algorithms based on a variety of assumptions and approaches. With the ever increasing demand for more accurate and powerful models, the inference problem remains of broad scientific interest. The abstract representation of biological systems through gene regulatory networks represents a powerful method to study such systems, encoding different amounts and types of information. In this review, we summarize the different types of inference algorithms specifically based on time-series transcriptomics, giving an overview of the main applications of gene regulatory networks in computational biology. This review is intended to give an updated reference of regulatory networks inference tools to biologists and researchers new to the topic and guide them in selecting the appropriate inference method that best fits their questions, aims, and experimental data.

SIGNET: Transcriptome-wide Causal Inference for Gene Regulatory Networks

Gene regulation plays an important role in understanding the mechanisms of human biology and diseases. However, inferring causal relationships between all genes is challenging due to the large number of genes in the transcriptome. Here, we present SIGNET (Statistical Inference on Gene Regulatory Networks), a flexible software package that reveals networks of causal regulation between genes built upon large-scale transcriptomic and genotypic data at the population level. Like Mendelian randomization, SIGNET uses genotypic variants as natural instrumental variables to establish such causal relationships but constructs a transcriptome-wide gene regulatory network with high confidence. SIGNET makes such a computationally heavy task feasible by deploying a well-designed statistical algorithm over a parallel computing environment. It also provides a user-friendly interface allowing for parameter tuning, efficient parallel computing scheduling, interactive network visualization, and confirmatory results retrieval. The Open source SIGNET software is freely available (https://www.zstats.org/signet/).

SCING: Inference of robust, interpretable gene regulatory networks from single cell and spatial transcriptomics

Gene regulatory network (GRN) inference is an integral part of understanding physiology and disease. Single cell/nuclei RNA-seq (scRNA-seq/snRNA-seq) data has been used to elucidate cell-type GRNs; however, the accuracy and speed of current scRNAseq-based GRN approaches are suboptimal. Here, we present Single Cell INtegrative Gene regulatory network inference (SCING), a gradient boosting and mutual information-based approach for identifying robust GRNs from scRNA-seq, snRNA-seq, and spatial transcriptomics data. Performance evaluation using Perturb-seq datasets, held-out data, and the mouse cell atlas combined with the DisGeNET database demonstrates the improved accuracy and biological interpretability of SCING compared to existing methods. We applied SCING to the entire mouse single cell atlas, human Alzheimer's disease (AD), and mouse AD spatial transcriptomics. SCING GRNs reveal unique disease subnetwork modeling capabilities, have intrinsic capacity to correct for batch effects, retrieve disease relevant genes and pathways, and are informative on spatial specificity of disease pathogenesis.

Gene regulatory network inference based on a nonhomogeneous dynamic Bayesian network model with an improved Markov Monte Carlo sampling

A nonhomogeneous dynamic Bayesian network model, which combines the dynamic Bayesian network and the multi-change point process, solves the limitations of the dynamic Bayesian network in modeling non-stationary gene expression data to a certain extent. However, certain problems persist, such as the low network reconstruction accuracy and poor model convergence. Therefore, we propose an MD-birth move based on the Manhattan distance of the data points to increase the rationality of the multi-change point process. The underlying concept of the MD-birth move is that the direction of movement of the change point is assumed to have a larger Manhattan distance between the variance and the mean of its left and right data points. Considering the data instability characteristics, we propose a Markov chain Monte Carlo sampling method based on node-dependent particle filtering in addition to the multi-change point process. The candidate parent nodes to be sampled, which are close to the real state, are pushed to the high probability area through the particle filter, and the candidate parent node set to be sampled that is far from the real state is pushed to the low probability area and then sampled. In terms of reconstructing the gene regulatory network, the model proposed in this paper (FC-DBN) has better network reconstruction accuracy and model convergence speed than other corresponding models on the Saccharomyces cerevisiae data and RAF data.

Cellular signaling pathways as plastic, proto-cognitive systems: Implications for biomedicine

Many aspects of health and disease are modeled using the abstraction of a "pathway"-a set of protein or other subcellular activities with specified functional linkages between them. This metaphor is a paradigmatic case of a deterministic, mechanistic framework that focuses biomedical intervention strategies on altering the members of this network or the up-/down-regulation links between them-rewiring the molecular hardware. However, protein pathways and transcriptional networks exhibit interesting and unexpected capabilities such as trainability (memory) and information processing in a context-sensitive manner. Specifically, they may be amenable to manipulation via their history of stimuli (equivalent to experiences in behavioral science). If true, this would enable a new class of biomedical interventions that target aspects of the dynamic physiological "software" implemented by pathways and gene-regulatory networks. Here, we briefly review clinical and laboratory data that show how high-level cognitive inputs and mechanistic pathway modulation interact to determine outcomes in vivo. Further, we propose an expanded view of pathways from the perspective of basal cognition and argue that a broader understanding of pathways and how they process contextual information across scales will catalyze progress in many areas of physiology and neurobiology. We argue that this fuller understanding of the functionality and tractability of pathways must go beyond a focus on the mechanistic details of protein and drug structure to encompass their physiological history as well as their embedding within higher levels of organization in the organism, with numerous implications for data science addressing health and disease. Exploiting tools and concepts from behavioral and cognitive sciences to explore a proto-cognitive metaphor for the pathways underlying health and disease is more than a philosophical stance on biochemical processes; at stake is a new roadmap for overcoming the limitations of today's pharmacological strategies and for inferring future therapeutic interventions for a wide range of disease states.

A gene regulatory network inference model based on pseudo-siamese network

MOTIVATION

Gene regulatory networks (GRNs) arise from the intricate interactions between transcription factors (TFs) and their target genes during the growth and development of organisms. The inference of GRNs can unveil the underlying gene interactions in living systems and facilitate the investigation of the relationship between gene expression patterns and phenotypic traits. Although several machine-learning models have been proposed for inferring GRNs from single-cell RNA sequencing (scRNA-seq) data, some of these models, such as Boolean and tree-based networks, suffer from sensitivity to noise and may encounter difficulties in handling the high noise and dimensionality of actual scRNA-seq data, as well as the sparse nature of gene regulation relationships. Thus, inferring large-scale information from GRNs remains a formidable challenge.

RESULTS

This study proposes a multilevel, multi-structure framework called a pseudo-Siamese GRN (PSGRN) for inferring large-scale GRNs from time-series expression datasets. Based on the pseudo-Siamese network, we applied a gated recurrent unit to capture the time features of each TF and target matrix and learn the spatial features of the matrices after merging by applying the DenseNet framework. Finally, we applied a sigmoid function to evaluate interactions. We constructed two maize sub-datasets, including gene expression levels and GRNs, using existing open-source maize multi-omics data and compared them to other GRN inference methods, including GENIE3, GRNBoost2, nonlinear ordinary differential equations, CNNC, and DGRNS. Our results show that PSGRN outperforms state-of-the-art methods. This study proposed a new framework: a PSGRN that allows GRNs to be inferred from scRNA-seq data, elucidating the temporal and spatial features of TFs and their target genes. The results show the model's robustness and generalization, laying a theoretical foundation for maize genotype-phenotype associations with implications for breeding work.

Interactions between genes altered during cardiotoxicity and neurotoxicity in zebrafish revealed using induced network modules analysis

As the manufacturing and development of new synthetic compounds increase to keep pace with the expanding global demand, adverse health effects due to these compounds are emerging as critical public health concerns. Zebrafish have become a prominent model organism to study toxicology due to their genomic similarity to humans, optical clarity, well-defined developmental stages, short generation time, and cost-effective maintenance. It also provides a shorter time frame for in vivo toxicology evaluation compared to the mammalian experimental systems. Here, we used meta-analysis to examine the alteration in genes during cardiotoxicity and neurotoxicity in zebrafish, caused by chemical exposure of any kind. First, we searched the literature comprehensively for genes that are altered during neurotoxicity and cardiotoxicity followed by meta-analysis using ConsensusPathDB. Since constant communication between the heart and the brain is an important physiological phenomenon, we also analyzed interactions among genes altered simultaneously during cardiotoxicity and neurotoxicity using induced network modules analysis in ConsensusPathDB. We observed inflammation and regeneration as the major pathways involved in cardiotoxicity and neurotoxicity. A large number of intermediate genes and input genes anchored in these pathways are molecular regulators of cell cycle progression and cell death and are implicated in tumor manifestation. We propose potential predictive biomarkers for neurotoxicity and cardiotoxicity and the major pathways potentially implicated in the manifestation of a particular toxicity phenotype.

Network medicine: an approach to complex kidney disease phenotypes

Scientific reductionism has been the basis of disease classification and understanding for more than a century. However, the reductionist approach of characterizing diseases from a limited set of clinical observations and laboratory evaluations has proven insufficient in the face of an exponential growth in data generated from transcriptomics, proteomics, metabolomics and deep phenotyping. A new systematic method is necessary to organize these datasets and build new definitions of what constitutes a disease that incorporates both biological and environmental factors to more precisely describe the ever-growing complexity of phenotypes and their underlying molecular determinants. Network medicine provides such a conceptual framework to bridge these vast quantities of data while providing an individualized understanding of disease. The modern application of network medicine principles is yielding new insights into the pathobiology of chronic kidney diseases and renovascular disorders by expanding the understanding of pathogenic mediators, novel biomarkers and new options for renal therapeutics. These efforts affirm network medicine as a robust paradigm for elucidating new advances in the diagnosis and treatment of kidney disorders.

STGRNS: an interpretable transformer-based method for inferring gene regulatory networks from single-cell transcriptomic data

MOTIVATION

Single-cell RNA-sequencing (scRNA-seq) technologies provide an opportunity to infer cell-specific gene regulatory networks (GRNs), which is an important challenge in systems biology. Although numerous methods have been developed for inferring GRNs from scRNA-seq data, it is still a challenge to deal with cellular heterogeneity.

RESULTS

To address this challenge, we developed an interpretable transformer-based method namely STGRNS for inferring GRNs from scRNA-seq data. In this algorithm, gene expression motif technique was proposed to convert gene pairs into contiguous sub-vectors, which can be used as input for the transformer encoder. By avoiding missing phase-specific regulations in a network, gene expression motif can improve the accuracy of GRN inference for different types of scRNA-seq data. To assess the performance of STGRNS, we implemented the comparative experiments with some popular methods on extensive benchmark datasets including 21 static and 27 time-series scRNA-seq dataset. All the results show that STGRNS is superior to other comparative methods. In addition, STGRNS was also proved to be more interpretable than "black box" deep learning methods, which are well-known for the difficulty to explain the predictions clearly.

AVAILABILITY AND IMPLEMENTATION

The source code and data are available at https://github.com/zhanglab-wbgcas/STGRNS.

De-novo reconstruction and identification of transcriptional gene regulatory network modules differentiating single-cell clusters

Single-cell RNA sequencing (scRNA-seq) technology provides an unprecedented opportunity to understand gene functions and interactions at single-cell resolution. While computational tools for scRNA-seq data analysis to decipher differential gene expression profiles and differential pathway expression exist, we still lack methods to learn differential regulatory disease mechanisms directly from the single-cell data. Here, we provide a new methodology, named DiNiro, to unravel such mechanisms de novo and report them as small, easily interpretable transcriptional regulatory network modules. We demonstrate that DiNiro is able to uncover novel, relevant, and deep mechanistic models that not just predict but explain differential cellular gene expression programs. DiNiro is available at https://exbio.wzw.tum.de/diniro/.

On the gene expression landscape of cancer

Kauffman picture of normal and tumor states as attractors in an abstract state space is used in order to interpret gene expression data for 15 cancer localizations obtained from The Cancer Genome Atlas. A principal component analysis of this data unveils the following qualitative aspects about tumors: 1) The state of a tissue in gene expression space can be described by a few variables. In particular, there is a single variable describing the progression from a normal tissue to a tumor. 2) Each cancer localization is characterized by a gene expression profile, in which genes have specific weights in the definition of the cancer state. There are no less than 2500 differentially-expressed genes, which lead to power-like tails in the expression distribution functions. 3) Tumors in different localizations share hundreds or even thousands of differentially expressed genes. There are 6 genes common to the 15 studied tumor localizations. 4) The tumor region is a kind of attractor. Tumors in advanced stages converge to this region independently of patient age or genetic characteristics. 5) There is a landscape of cancer in gene expression space with an approximate border separating normal tissues from tumors.

Deciphering sperm functions using biological networks

The global human population is exponentially increasing, which requires the production of quality food through efficient reproduction as well as sustainable production of livestock. Lack of knowledge and technology for assessing semen quality and predicting bull fertility is hindering advances in animal science and food animal production and causing millions of dollars of economic losses annually. The intent of this systemic review is to summarize methods from computational biology for analysis of gene, metabolite, and protein networks to identify potential markers that can be applied to improve livestock reproduction, with a focus on bull fertility. We provide examples of available gene, metabolic, and protein networks and computational biology methods to show how the interactions between genes, proteins, and metabolites together drive the complex process of spermatogenesis and regulate fertility in animals. We demonstrate the use of the National Center for Biotechnology Information (NCBI) and Ensembl for finding gene sequences, and then use them to create and understand gene, protein and metabolite networks for sperm associated factors to elucidate global cellular processes in sperm. This study highlights the value of mapping complex biological pathways among livestock and potential for conducting studies on promoting livestock improvement for global food security.

SPASCER: spatial transcriptomics annotation at single-cell resolution

In recent years, the explosive growth of spatial technologies has enabled the characterization of spatial heterogeneity of tissue architectures. Compared to traditional sequencing, spatial transcriptomics reserves the spatial information of each captured location and provides novel insights into diverse spatially related biological contexts. Even though two spatial transcriptomics databases exist, they provide limited analytical information. Information such as spatial heterogeneity of genes and cells, cell-cell communication activities in space, and the cell type compositions in the microenvironment are critical clues to unveil the mechanism of tumorigenesis and embryo differentiation. Therefore, we constructed a new spatial transcriptomics database, named SPASCER (https://ccsm.uth.edu/SPASCER), designed to help understand the heterogeneity of tissue organizations, region-specific microenvironment, and intercellular interactions across tissue architectures at multiple levels. SPASCER contains datasets from 43 studies, including 1082 sub-datasets from 16 organ types across four species. scRNA-seq was integrated to deconvolve/map spatial transcriptomics, and processed with spatial cell-cell interaction, gene pattern and pathway enrichment analysis. Cell-cell interactions and gene regulation network of scRNA-seq from matched spatial transcriptomics were performed as well. The application of SPASCER will provide new insights into tissue architecture and a solid foundation for the mechanistic understanding of many biological processes in healthy and diseased tissues.

A Consensus Gene Regulatory Network for Neurodegenerative Diseases Using Single-Cell RNA-Seq Data

Gene regulatory network (GRN) inference from gene expression data is a highly complex and challenging task in systems biology. Despite the challenges, GRNs have emerged, and for complex diseases such as neurodegenerative diseases, they have the potential to provide vital information and identify key regulators. However, every GRN method produced predicts results based on its assumptions, providing limited biological insights. For that reason, the current work focused on the development of an ensemble method from individual GRN methods to address this issue. Four state-of-the-art GRN algorithms were selected to form a consensus GRN from their common gene interactions. Each algorithm uses a different construction method, and for a more robust behavior, both static and dynamic methods were selected as well. The algorithms were applied to a scRNA-seq dataset from the CK-p25 mus musculus model during neurodegeneration. The top subnetworks were constructed from the consensus network, and potential key regulators were identified. The results also demonstrated the overlap between the algorithms for the current dataset and the necessity for an ensemble approach. This work aims to demonstrate the creation of an ensemble network and provide insights into whether a combination of different GRN methods can produce valuable results.

Network-based approaches for modeling disease regulation and progression

Molecular interaction networks lay the foundation for studying how biological functions are controlled by the complex interplay of genes and proteins. Investigating perturbed processes using biological networks has been instrumental in uncovering mechanisms that underlie complex disease phenotypes. Rapid advances in omics technologies have prompted the generation of high-throughput datasets, enabling large-scale, network-based analyses. Consequently, various modeling techniques, including network enrichment, differential network extraction, and network inference, have proven to be useful for gaining new mechanistic insights. We provide an overview of recent network-based methods and their core ideas to facilitate the discovery of disease modules or candidate mechanisms. Knowledge generated from these computational efforts will benefit biomedical research, especially drug development and precision medicine. We further discuss current challenges and provide perspectives in the field, highlighting the need for more integrative and dynamic network approaches to model disease development and progression.

Global Distribution and Diversity of Prevalent Sewage Water Plasmidomes

Sewage water from around the world contains an abundance of short plasmids, several of which harbor antimicrobial resistance genes (ARGs). The global dynamics of plasmid-derived antimicrobial resistance and functions are only starting to be unveiled. Here, we utilized a previously created data set of 159,332 assumed small plasmids from 24 different global sewage samples. The detailed phylogeny, as well as the interplay between their protein domains, ARGs, and predicted bacterial host genera, were investigated to understand sewage plasmidome dynamics globally. A total of 58,429 circular elements carried genes encoding plasmid-related features, and MASH distance analyses showed a high degree of diversity. A single (yet diverse) cluster of 520 predicted Acinetobacter plasmids was predominant among the European sewage water. Our results suggested a prevalence of plasmid-backbone gene combinations over others. This could be related to selected bacterial genera that act as bacterial hosts. These combinations also mirrored the geographical locations of the sewage samples. Our functional domain network analysis identified three groups of plasmids. However, these backbone domains were not exclusive to any given group, and Acinetobacter was the dominant host genus among the theta-replicating plasmids, which contained a reservoir of the macrolide resistance gene pair msr(E) and mph(E). Macrolide resistance genes were the most common in the sewage plasmidomes and were found in the largest number of unique plasmids. While msr(E) and mph(E) were limited to Acinetobacter, erm(B) was disseminated among a range of Firmicutes plasmids, including Staphylococcus and Streptococcus, highlighting a potential reservoir of antibiotic resistance for these pathogens from around the globe. IMPORTANCE Antimicrobial resistance is a global threat to human health, as it inhibits our ability to treat infectious diseases. This study utilizes sewage water plasmidomes to identify plasmid-derived features and highlights antimicrobial resistance genes, particularly macrolide resistance genes, as abundant in sewage water plasmidomes in Firmicutes and Acinetobacter hosts. The emergence of macrolide resistance in these bacteria suggests that macrolide selective pressure exists in sewage water and that the resident bacteria can readily acquire macrolide resistance via small plasmids.

HOXA9 has the hallmarks of a biological switch with implications in blood cancers

Blood malignancies arise from the dysregulation of haematopoiesis. The type of blood cell and the specific order of oncogenic events initiating abnormal growth ultimately determine the cancer subtype and subsequent clinical outcome. HOXA9 plays an important role in acute myeloid leukaemia (AML) prognosis by promoting blood cell expansion and altering differentiation; however, the function of HOXA9 in other blood malignancies is still unclear. Here, we highlight the biological switch and prognosis marker properties of HOXA9 in AML and chronic myeloproliferative neoplasms (MPN). First, we establish the ability of HOXA9 to stratify AML patients with distinct cellular and clinical outcomes. Then, through the use of a computational network model of MPN, we show that the self-activation of HOXA9 and its relationship to JAK2 and TET2 can explain the branching progression of JAK2/TET2 mutant MPN patients towards divergent clinical characteristics. Finally, we predict a connection between the RUNX1 and MYB genes and a suppressive role for the NOTCH pathway in MPN diseases.

Enhancer RNA (eRNA) in Human Diseases

Enhancer RNAs (eRNAs), a class of non-coding RNAs (ncRNAs) transcribed from enhancer regions, serve as a type of critical regulatory element in gene expression. There is increasing evidence demonstrating that the aberrant expression of eRNAs can be broadly detected in various human diseases. Some studies also revealed the potential clinical utility of eRNAs in these diseases. In this review, we summarized the recent studies regarding the pathological mechanisms of eRNAs as well as their potential utility across human diseases, including cancers, neurodegenerative disorders, cardiovascular diseases and metabolic diseases. It could help us to understand how eRNAs are engaged in the processes of diseases and to obtain better insight of eRNAs in diagnosis, prognosis or therapy. The studies we reviewed here indicate the enormous therapeutic potency of eRNAs across human diseases.

A Joint estimation approach to sparse additive ordinary differential equations

Ordinary differential equations (ODEs) are widely used to characterize the dynamics of complex systems in real applications. In this article, we propose a novel joint estimation approach for generalized sparse additive ODEs where observations are allowed to be non-Gaussian. The new method is unified with existing collocation methods by considering the likelihood, ODE fidelity and sparse regularization simultaneously. We design a block coordinate descent algorithm for optimizing the non-convex and non-differentiable objective function. The global convergence of the algorithm is established. The simulation study and two applications demonstrate the superior performance of the proposed method in estimation and improved performance of identifying the sparse structure.

Consensus clustering for Bayesian mixture models

BACKGROUND

Cluster analysis is an integral part of precision medicine and systems biology, used to define groups of patients or biomolecules. Consensus clustering is an ensemble approach that is widely used in these areas, which combines the output from multiple runs of a non-deterministic clustering algorithm. Here we consider the application of consensus clustering to a broad class of heuristic clustering algorithms that can be derived from Bayesian mixture models (and extensions thereof) by adopting an early stopping criterion when performing sampling-based inference for these models. While the resulting approach is non-Bayesian, it inherits the usual benefits of consensus clustering, particularly in terms of computational scalability and providing assessments of clustering stability/robustness.

RESULTS

In simulation studies, we show that our approach can successfully uncover the target clustering structure, while also exploring different plausible clusterings of the data. We show that, when a parallel computation environment is available, our approach offers significant reductions in runtime compared to performing sampling-based Bayesian inference for the underlying model, while retaining many of the practical benefits of the Bayesian approach, such as exploring different numbers of clusters. We propose a heuristic to decide upon ensemble size and the early stopping criterion, and then apply consensus clustering to a clustering algorithm derived from a Bayesian integrative clustering method. We use the resulting approach to perform an integrative analysis of three 'omics datasets for budding yeast and find clusters of co-expressed genes with shared regulatory proteins. We validate these clusters using data external to the analysis.

CONCLUSTIONS

Our approach can be used as a wrapper for essentially any existing sampling-based Bayesian clustering implementation, and enables meaningful clustering analyses to be performed using such implementations, even when computational Bayesian inference is not feasible, e.g. due to poor exploration of the target density (often as a result of increasing numbers of features) or a limited computational budget that does not along sufficient samples to drawn from a single chain. This enables researchers to straightforwardly extend the applicability of existing software to much larger datasets, including implementations of sophisticated models such as those that jointly model multiple datasets.

Discovering gene regulatory networks of multiple phenotypic groups using dynamic Bayesian networks

Dynamic Bayesian networks (DBNs) can be used for the discovery of gene regulatory networks (GRNs) from time series gene expression data. Here, we suggest a strategy for learning DBNs from gene expression data by employing a Bayesian approach that is scalable to large networks and is targeted at learning models with high predictive accuracy. Our framework can be used to learn DBNs for multiple groups of samples and highlight differences and similarities in their GRNs. We learn these DBN models based on different structural and parametric assumptions and select the optimal model based on the cross-validated predictive accuracy. We show in simulation studies that our approach is better equipped to prevent overfitting than techniques used in previous studies. We applied the proposed DBN-based approach to two time series transcriptomic datasets from the Gene Expression Omnibus database, each comprising data from distinct phenotypic groups of the same tissue type. In the first case, we used DBNs to characterize responders and non-responders to anti-cancer therapy. In the second case, we compared normal to tumor cells of colorectal tissue. The classification accuracy reached by the DBN-based classifier for both datasets was higher than reported previously. For the colorectal cancer dataset, our analysis suggested that GRNs for cancer and normal tissues have a lot of differences, which are most pronounced in the neighborhoods of oncogenes and known cancer tissue markers. The identified differences in gene networks of cancer and normal cells may be used for the discovery of targeted therapies.

Approaches in Gene Coexpression Analysis in Eukaryotes

Simple Summary

Genes whose expression levels rise and fall similarly in a large set of samples, may be considered coexpressed. Gene coexpression analysis refers to the en masse discovery of coexpressed genes from a large variety of transcriptomic experiments. The type of biological networks that studies gene coexpression, known as Gene Coexpression Networks, consist of an undirected graph depicting genes and their coexpression relationships. Coexpressed genes are clustered in smaller subnetworks, the predominant biological roles of which can be determined through enrichment analysis. By studying well-annotated gene partners, the attribution of new roles to genes of unknown function or assumption for participation in common metabolic pathways can be achieved, through a guilt-by-association approach. In this review, we present key issues in gene coexpression analysis, as well as the most popular tools that perform it.

Abstract

Gene coexpression analysis constitutes a widely used practice for gene partner identification and gene function prediction, consisting of many intricate procedures. The analysis begins with the collection of primary transcriptomic data and their preprocessing, continues with the calculation of the similarity between genes based on their expression values in the selected sample dataset and results in the construction and visualisation of a gene coexpression network (GCN) and its evaluation using biological term enrichment analysis. As gene coexpression analysis has been studied extensively, we present most parts of the methodology in a clear manner and the reasoning behind the selection of some of the techniques. In this review, we offer a comprehensive and comprehensible account of the steps required for performing a complete gene coexpression analysis in eukaryotic organisms. We comment on the use of RNA-Seq vs. microarrays, as well as the best practices for GCN construction. Furthermore, we recount the most popular webtools and standalone applications performing gene coexpression analysis, with details on their methods, features and outputs.

GRNbenchmark - a web server for benchmarking directed gene regulatory network inference methods

Accurate inference of gene regulatory networks (GRN) is an essential component of systems biology, and there is a constant development of new inference methods. The most common approach to assess accuracy for publications is to benchmark the new method against a selection of existing algorithms. This often leads to a very limited comparison, potentially biasing the results, which may stem from tuning the benchmark's properties or incorrect application of other methods. These issues can be avoided by a web server with a broad range of data properties and inference algorithms, that makes it easy to perform comprehensive benchmarking of new methods, and provides a more objective assessment. Here we present https://GRNbenchmark.org/ - a new web server for benchmarking GRN inference methods, which provides the user with a set of benchmarks with several datasets, each spanning a range of properties including multiple noise levels. As soon as the web server has performed the benchmarking, the accuracy results are made privately available to the user via interactive summary plots and underlying curves. The user can then download these results for any purpose, and decide whether or not to make them public to share with the community.

Network Approaches for Charting the Transcriptomic and Epigenetic Landscape of the Developmental Origins of Health and Disease

The early developmental phase is of critical importance for human health and disease later in life. To decipher the molecular mechanisms at play, current biomedical research is increasingly relying on large quantities of diverse omics data. The integration and interpretation of the different datasets pose a critical challenge towards the holistic understanding of the complex biological processes that are involved in early development. In this review, we outline the major transcriptomic and epigenetic processes and the respective datasets that are most relevant for studying the periconceptional period. We cover both basic data processing and analysis steps, as well as more advanced data integration methods. A particular focus is given to network-based methods. Finally, we review the medical applications of such integrative analyses.

Application of modular response analysis to medium- to large-size biological systems

The development of high-throughput genomic technologies associated with recent genetic perturbation techniques such as short hairpin RNA (shRNA), gene trapping, or gene editing (CRISPR/Cas9) has made it possible to obtain large perturbation data sets. These data sets are invaluable sources of information regarding the function of genes, and they offer unique opportunities to reverse engineer gene regulatory networks in specific cell types. Modular response analysis (MRA) is a well-accepted mathematical modeling method that is precisely aimed at such network inference tasks, but its use has been limited to rather small biological systems so far. In this study, we show that MRA can be employed on large systems with almost 1,000 network components. In particular, we show that MRA performance surpasses general-purpose mutual information-based algorithms. Part of these competitive results was obtained by the application of a novel heuristic that pruned MRA-inferred interactions a posteriori. We also exploited a block structure in MRA linear algebra to parallelize large system resolutions.

A one-dimensional parameter-free model for carcinogenesis in gene expression space

A small portion of a tissue defines a microstate in gene expression space. Mutations, epigenetic events or external factors cause microstate displacements which are modeled by combining small independent gene expression variations and large Levy jumps, resulting from the collective variations of a set of genes. The risk of cancer in a tissue is estimated as the microstate probability to transit from the normal to the tumor region in gene expression space. The formula coming from the contribution of large Levy jumps seems to provide a qualitatively correct description of the lifetime risk of cancer in 8 tissues, and reveals an interesting connection between the risk and the way the tissue is protected against infections.