Integrative approaches for finding modular structure in biological networks

Harnessing network pharmacology in drug discovery: an integrated approach

Traditional drug discovery approach is based on one drug-one target, that is associated with very lengthy timelines, high costs and very low success rates. Network pharmacology (NP) is a novel method of drug designing, that is based on a multiple-target approach. NP integrates systems such as biology, pharmacology and computational techniques to address the limitations of traditional methods of drug discovery. With help of mapping biological networks, it provides deep insights into biological molecules' interactions and enhances our understanding to the mechanism of drugs, polypharmacology and disease etiology. This review explores the theoretical framework of network pharmacology, discussing the principles and methodologies that enable the construction of drug-target and disease-gene networks. It highlights how data mining, bioinformatics tools and computational models are utilised to predict drug behaviour, repurpose existing drugs and identify novel therapeutic targets. Applications of network pharmacology in the treatment of complex diseases-such as cancer, neurodegenerative disorders, cardiovascular diseases and infectious diseases-are extensively covered, demonstrating its potential to identify multi-target drugs for multifaceted disease mechanisms. Despite the promising results, NP faces challenges due to incomplete and quality of biological data, computational complexities and biological system redundancy. It also faces regulatory challenges in drug approval, demanding revision in regulatory guidelines towards multi-target therapies. Advancements in AI and machine learning, dynamic network modelling and global collaboration can further enhance the efficacy of network pharmacology. This integrative approach has the potential to revolutionise drug discovery, offering new solutions for personalised medicine, drug repurposing and tackling the complexities of modern diseases.

From Data to Cure: A Comprehensive Exploration of Multi-omics Data Analysis for Targeted Therapies

In the dynamic landscape of targeted therapeutics, drug discovery has pivoted towards understanding underlying disease mechanisms, placing a strong emphasis on molecular perturbations and target identification. This paradigm shift, crucial for drug discovery, is underpinned by big data, a transformative force in the current era. Omics data, characterized by its heterogeneity and enormity, has ushered biological and biomedical research into the big data domain. Acknowledging the significance of integrating diverse omics data strata, known as multi-omics studies, researchers delve into the intricate interrelationships among various omics layers. This review navigates the expansive omics landscape, showcasing tailored assays for each molecular layer through genomes to metabolomes. The sheer volume of data generated necessitates sophisticated informatics techniques, with machine-learning (ML) algorithms emerging as robust tools. These datasets not only refine disease classification but also enhance diagnostics and foster the development of targeted therapeutic strategies. Through the integration of high-throughput data, the review focuses on targeting and modeling multiple disease-regulated networks, validating interactions with multiple targets, and enhancing therapeutic potential using network pharmacology approaches. Ultimately, this exploration aims to illuminate the transformative impact of multi-omics in the big data era, shaping the future of biological research.

A collaborative network analysis for the interpretation of transcriptomics data in Huntington's disease

Rare diseases may affect the quality of life of patients and be life-threatening. Therapeutic opportunities are often limited, in part because of the lack of understanding of the molecular mechanisms underlying these diseases. This can be ascribed to the low prevalence of rare diseases and therefore the lower sample sizes available for research. A way to overcome this is to integrate experimental rare disease data with prior knowledge using network-based methods. Taking this one step further, we hypothesized that combining and analyzing the results from multiple network-based methods could provide data-driven hypotheses of pathogenic mechanisms from multiple perspectives.We analyzed a Huntington's disease transcriptomics dataset using six network-based methods in a collaborative way. These methods either inherently reported enriched annotation terms or their results were fed into enrichment analyses. The resulting significantly enriched Reactome pathways were then summarized using the ontological hierarchy which allowed the integration and interpretation of outputs from multiple methods. Among the resulting enriched pathways, there are pathways that have been shown previously to be involved in Huntington's disease and pathways whose direct contribution to disease pathogenesis remains unclear and requires further investigation.In summary, our study shows that collaborative network analysis approaches are well-suited to study rare diseases, as they provide hypotheses for pathogenic mechanisms from multiple perspectives. Applying different methods to the same case study can uncover different disease mechanisms that would not be apparent with the application of a single method.

Distinct network patterns emerge from Cartesian and XOR epistasis models: a comparative network science analysis

BACKGROUND

Epistasis, the phenomenon where the effect of one gene (or variant) is masked or modified by one or more other genes, significantly contributes to the phenotypic variance of complex traits. Traditionally, epistasis has been modeled using the Cartesian epistatic model, a multiplicative approach based on standard statistical regression. However, a recent study investigating epistasis in obesity-related traits has identified potential limitations of the Cartesian epistatic model, revealing that it likely only detects a fraction of the genetic interactions occurring in natural systems. In contrast, the exclusive-or (XOR) epistatic model has shown promise in detecting a broader range of epistatic interactions and revealing more biologically relevant functions associated with interacting variants. To investigate whether the XOR epistatic model also forms distinct network structures compared to the Cartesian model, we applied network science to examine genetic interactions underlying body mass index (BMI) in rats (Rattus norvegicus).

RESULTS

Our comparative analysis of XOR and Cartesian epistatic models in rats reveals distinct topological characteristics. The XOR model exhibits enhanced sensitivity to epistatic interactions between the network communities found in the Cartesian epistatic network, facilitating the identification of novel trait-related biological functions via community-based enrichment analysis. Additionally, the XOR network features triangle network motifs, indicative of higher-order epistatic interactions. This research also evaluates the impact of linkage disequilibrium (LD)-based edge pruning on network-based epistasis analysis, finding that LD-based edge pruning may lead to increased network fragmentation, which may hinder the effectiveness of network analysis for the investigation of epistasis. We confirmed through network permutation analysis that most XOR and Cartesian epistatic networks derived from the data display distinct structural properties compared to randomly shuffled networks.

CONCLUSIONS

Collectively, these findings highlight the XOR model's ability to uncover meaningful biological associations and higher-order epistasis derived from lower-order network topologies. The introduction of community-based enrichment analysis and motif-based epistatic discovery emphasize network science as a critical approach for advancing epistasis research and understanding complex genetic architectures.

GNN4DM: a graph neural network-based method to identify overlapping functional disease modules

MOTIVATION

Identifying disease modules within molecular interaction networks is an essential exploratory step in computational biology, offering insights into disease mechanisms and potential therapeutic targets. Traditional methods often struggle with the inherent complexity and overlapping nature of biological networks, and they are limited in effectively leveraging the vast amount of available genomic data and biological knowledge. This limitation underscores the need for more effective, automated approaches to integrate these rich data sources.

RESULTS

In this work, we propose GNN4DM, a novel graph neural network-based structured model that automates the discovery of overlapping functional disease modules. GNN4DM effectively integrates network topology with genomic data to learn the representations of the genes corresponding to functional modules and align these with known biological pathways for enhanced interpretability. Following the DREAM benchmark evaluation setting and extending with three independent data sources (GWAS Atlas, FinnGen, and DisGeNET), we show that GNN4DM performs better than several state-of-the-art methods in detecting biologically meaningful modules. Moreover, we demonstrate the method's applicability by discovering two novel multimorbidity modules significantly enriched across a diverse range of seemingly unrelated diseases.

AVAILABILITY AND IMPLEMENTATION

Source code, all training data, and all identified disease modules are freely available for download at https://github.com/gezsi/gnn4dm. GNN4DM was implemented in Python.

Networks as tools for defining emergent properties of microbiomes and their stability

The potential promise of the microbiome to ameliorate a wide range of societal and ecological challenges, from disease prevention and treatment to the restoration of entire ecosystems, hinges not only on microbiome engineering but also on the stability of beneficial microbiomes. Yet the properties of microbiome stability remain elusive and challenging to discern due to the complexity of interactions and often intractable diversity within these communities of bacteria, archaea, fungi, and other microeukaryotes. Networks are powerful tools for the study of complex microbiomes, with the potential to elucidate structural patterns of stable communities and generate testable hypotheses for experimental validation. However, the implementation of these analyses introduces a cascade of dichotomies and decision trees due to the lack of consensus on best practices. Here, we provide a road map for network-based microbiome studies with an emphasis on discerning properties of stability. We identify important considerations for data preparation, network construction, and interpretation of network properties. We also highlight remaining limitations and outstanding needs for this field. This review also serves to clarify the varying schools of thought on the application of network theory for microbiome studies and to identify practices that enhance the reproducibility and validity of future work. Video Abstract.

Urbanisation-associated shifts in the avian metabolome within the annual cycle

While organisms have evolved to cope with predictable changes in the environment, the rapid rate of current global change presents numerous novel and unpredictable stressors to which organisms have had less time to adapt. To persist in the urban environment, organisms must modify their physiology, morphology and behaviour accordingly. Metabolomics offers great potential for characterising organismal responses to natural and anthropogenic stressors at the systems level and can be applied to any species, even without genomic knowledge. Using metabolomic profiling of blood, we investigated how two closely related species of passerine bird respond to the urban environment. Great tits Parus major and blue tits Cyanistes caeruleus residing in urban and forest habitats were sampled during the breeding (spring) and non-breeding (winter) seasons across replicated sites in southern Sweden. During breeding, differences in the plasma metabolome between urban and forest birds were characterised by higher levels of amino acids in urban-dwelling tits and higher levels of fatty acyls in forest-dwelling tits. The suggested higher rates of fatty acid oxidation in forest tits could be driven by habitat-associated differences in diet and could explain the higher reproductive investment and success of forest tits. High levels of amino acids in breeding urban tits could reflect the lack of lipid-rich caterpillars in the urban environment and a dietary switch to protein-rich spiders, which could be of benefit for tackling inflammation and oxidative stress associated with pollution. In winter, metabolomic profiles indicated lower overall levels of amino acids and fatty acyls in urban tits, which could reflect relaxed energetic demands in the urban environment. Our metabolomic profiling of two urban-adapted species suggests that their metabolism is modified by urban living, though whether these changes represent adaptative or non-adaptive mechanisms to cope with anthropogenic challenges remains to be determined.

AutoFocus: a hierarchical framework to explore multi-omic disease associations spanning multiple scales of biomolecular interaction

Recent advances in high-throughput measurement technologies have enabled the analysis of molecular perturbations associated with disease phenotypes at the multi-omic level. Such perturbations can range in scale from fluctuations of individual molecules to entire biological pathways. Data-driven clustering algorithms have long been used to group interactions into interpretable functional modules; however, these modules are typically constrained to a fixed size or statistical cutoff. Furthermore, modules are often analyzed independently of their broader biological context. Consequently, such clustering approaches limit the ability to explore functional module associations with disease phenotypes across multiple scales. Here, we introduce AutoFocus, a data-driven method that hierarchically organizes biomolecules and tests for phenotype enrichment at every level within the hierarchy. As a result, the method allows disease-associated modules to emerge at any scale. We evaluated this approach using two datasets: First, we explored associations of biomolecules from the multi-omic QMDiab dataset (n = 388) with the well-characterized type 2 diabetes phenotype. Secondly, we utilized the ROS/MAP Alzheimer's disease dataset (n = 500), consisting of high-throughput measurements of brain tissue to explore modules associated with multiple Alzheimer's Disease-related phenotypes. Our method identifies modules that are multi-omic, span multiple pathways, and vary in size. We provide an interactive tool to explore this hierarchy at different levels and probe enriched modules, empowering users to examine the full hierarchy, delve into biomolecular drivers of disease phenotype within a module, and incorporate functional annotations.

Current and future directions in network biology

Summary

Network biology is an interdisciplinary field bridging computational and biological sciences that has proved pivotal in advancing the understanding of cellular functions and diseases across biological systems and scales. Although the field has been around for two decades, it remains nascent. It has witnessed rapid evolution, accompanied by emerging challenges. These stem from various factors, notably the growing complexity and volume of data together with the increased diversity of data types describing different tiers of biological organization. We discuss prevailing research directions in network biology, focusing on molecular/cellular networks but also on other biological network types such as biomedical knowledge graphs, patient similarity networks, brain networks, and social/contact networks relevant to disease spread. In more detail, we highlight areas of inference and comparison of biological networks, multimodal data integration and heterogeneous networks, higher-order network analysis, machine learning on networks, and network-based personalized medicine. Following the overview of recent breakthroughs across these five areas, we offer a perspective on future directions of network biology. Additionally, we discuss scientific communities, educational initiatives, and the importance of fostering diversity within the field. This article establishes a roadmap for an immediate and long-term vision for network biology.

Availability and implementation

Not applicable.

DDN3.0: determining significant rewiring of biological network structure with differential dependency networks

MOTIVATION

Complex diseases are often caused and characterized by misregulation of multiple biological pathways. Differential network analysis aims to detect significant rewiring of biological network structures under different conditions and has become an important tool for understanding the molecular etiology of disease progression and therapeutic response. With few exceptions, most existing differential network analysis tools perform differential tests on separately learned network structures that are computationally expensive and prone to collapse when grouped samples are limited or less consistent.

RESULTS

We previously developed an accurate differential network analysis method-differential dependency networks (DDN), that enables joint learning of common and rewired network structures under different conditions. We now introduce the DDN3.0 tool that improves this framework with three new and highly efficient algorithms, namely, unbiased model estimation with a weighted error measure applicable to imbalance sample groups, multiple acceleration strategies to improve learning efficiency, and data-driven determination of proper hyperparameters. The comparative experimental results obtained from both realistic simulations and case studies show that DDN3.0 can help biologists more accurately identify, in a study-specific and often unknown conserved regulatory circuitry, a network of significantly rewired molecular players potentially responsible for phenotypic transitions.

AVAILABILITY AND IMPLEMENTATION

The Python package of DDN3.0 is freely available at https://github.com/cbil-vt/DDN3. A user's guide and a vignette are provided at https://ddn-30.readthedocs.io/.

Malignant peritoneal mesothelioma interactome with 417 novel protein-protein interactions

BACKGROUND

Malignant peritoneal mesothelioma (MPeM) is an aggressive cancer affecting the abdominal peritoneal lining and intra-abdominal organs, with a median survival of ~2.5 years.

METHODS

We constructed the protein interactome of 59 MPeM-associated genes with previously known protein-protein interactions (PPIs) as well as novel PPIs predicted using our previously developed HiPPIP computational model and analysed it for transcriptomic and functional associations and for repurposable drugs.

RESULTS

The MPeM interactome had over 400 computationally predicted PPIs and 4700 known PPIs. Transcriptomic evidence validated 75.6% of the genes in the interactome and 65% of the novel interactors. Some genes had tissue-specific expression in extramedullary hematopoietic sites and the expression of some genes could be correlated with unfavourable prognoses in various cancers. 39 out of 152 drugs that target the proteins in the interactome were identified as potentially repurposable for MPeM, with 29 having evidence from prior clinical trials, animal models or cell lines for effectiveness against peritoneal and pleural mesothelioma and primary peritoneal cancer. Functional modules related to chromosomal segregation, transcriptional dysregulation, IL-6 production and hematopoiesis were identified from the interactome. The MPeM interactome overlapped significantly with the malignant pleural mesothelioma interactome, revealing shared molecular pathways.

CONCLUSIONS

Our findings demonstrate the utility of the interactome in uncovering biological associations and in generating clinically translatable results.

Distinct Network Patterns Emerge from Cartesian and XOR Epistasis Models: A Comparative Network Science Analysis

Background

Epistasis, the phenomenon where the effect of one gene (or variant) is masked or modified by one or more other genes, can significantly contribute to the observed phenotypic variance of complex traits. To date, it has been generally assumed that genetic interactions can be detected using a Cartesian, or multiplicative, interaction model commonly utilized in standard regression approaches. However, a recent study investigating epistasis in obesity-related traits in rats and mice has identified potential limitations of the Cartesian model, revealing that it only detects some of the genetic interactions occurring in these systems. By applying an alternative approach, the exclusive-or (XOR) model, the researchers detected a greater number of epistatic interactions and identified more biologically relevant ontological terms associated with the interacting loci. This suggests that the XOR model may provide a more comprehensive understanding of epistasis in these species and phenotypes. To further explore these findings and determine if different interaction models also make up distinct epistatic networks, we leverage network science to provide a more comprehensive view into the genetic interactions underlying BMI in this system.

Results

Our comparative analysis of networks derived from Cartesian and XOR interaction models in rats (Rattus norvegicus) uncovers distinct topological characteristics for each model-derived network. Notably, we discover that networks based on the XOR model exhibit an enhanced sensitivity to epistatic interactions. This sensitivity enables the identification of network communities, revealing novel trait-related biological functions through enrichment analysis. Furthermore, we identify triangle network motifs in the XOR epistatic network, suggestive of higher-order epistasis, based on the topology of lower-order epistasis.

Conclusions

These findings highlight the XOR model's ability to uncover meaningful biological associations as well as higher-order epistasis from lower-order epistatic networks. Additionally, our results demonstrate that network approaches not only enhance epistasis detection capabilities but also provide more nuanced understandings of genetic architectures underlying complex traits. The identification of community structures and motifs within these distinct networks, especially in XOR, points to the potential for network science to aid in the discovery of novel genetic pathways and regulatory networks. Such insights are important for advancing our understanding of phenotype-genotype relationships.

A protein network refinement method based on module discovery and biological information

BACKGROUND

The identification of essential proteins can help in understanding the minimum requirements for cell survival and development to discover drug targets and prevent disease. Nowadays, node ranking methods are a common way to identify essential proteins, but the poor data quality of the underlying PIN has somewhat hindered the identification accuracy of essential proteins for these methods in the PIN. Therefore, researchers constructed refinement networks by considering certain biological properties of interacting protein pairs to improve the performance of node ranking methods in the PIN. Studies show that proteins in a complex are more likely to be essential than proteins not present in the complex. However, the modularity is usually ignored for the refinement methods of the PINs.

METHODS

Based on this, we proposed a network refinement method based on module discovery and biological information. The idea is, first, to extract the maximal connected subgraph in the PIN, and to divide it into different modules by using Fast-unfolding algorithm; then, to detect critical modules according to the orthologous information, subcellular localization information and topology information within each module; finally, to construct a more refined network (CM-PIN) by using the identified critical modules.

RESULTS

To evaluate the effectiveness of the proposed method, we used 12 typical node ranking methods (LAC, DC, DMNC, NC, TP, LID, CC, BC, PR, LR, PeC, WDC) to compare the overall performance of the CM-PIN with those on the S-PIN, D-PIN and RD-PIN. The experimental results showed that the CM-PIN was optimal in terms of the identification number of essential proteins, precision-recall curve, Jackknifing method and other criteria, and can help to identify essential proteins more accurately.

Community cohesion looseness in gene networks reveals individualized drug targets and resistance

Community cohesion plays a critical role in the determination of an individual's health in social science. Intriguingly, a community structure of gene networks indicates that the concept of community cohesion could be applied between the genes as well to overcome the limitations of single gene-based biomarkers for precision oncology. Here, we develop community cohesion scores which precisely quantify the community ability to retain the interactions between the genes and their cellular functions in each individualized gene network. Using breast cancer as a proof-of-concept study, we measure the community cohesion score profiles of 950 case samples and predict the individualized therapeutic targets in 2-fold. First, we prioritize them by finding druggable genes present in the community with the most and relatively decreased scores in each individual. Then, we pinpoint more individualized therapeutic targets by discovering the genes which greatly contribute to the community cohesion looseness in each individualized gene network. Compared with the previous approaches, the community cohesion scores show at least four times higher performance in predicting effective individualized chemotherapy targets based on drug sensitivity data. Furthermore, the community cohesion scores successfully discover the known breast cancer subtypes and we suggest new targeted therapy targets for triple negative breast cancer (e.g. KIT and GABRP). Lastly, we demonstrate that the community cohesion scores can predict tamoxifen responses in ER+ breast cancer and suggest potential combination therapies (e.g. NAMPT and RXRA inhibitors) to reduce endocrine therapy resistance based on individualized characteristics. Our method opens new perspectives for the biomarker development in precision oncology.

Pungency related gene network in Allium sativum L., response to sulfur treatments

Pungency of garlic (Allium sativum L.) is generated from breakdown of the alk(en)yl cysteine sulphoxide (CSO), alliin and its subsequent breakdown to allicin under the activity of alliinase (All). Based on recent evidence, two other important genes including Sulfite reductase (SiR) and Superoxide dismutase (SOD) are thought to be related to sulfur metabolism. These three gene functions are in sulfate assimilation pathway. However, whether it is involved in stress response in crops is largely unknown. In this research, the order and priority of simultaneous expression of three genes including All, SiR and SOD were measured on some garlic ecotypes of Iran, collected from Zanjan, Hamedan and Gilan, provinces under sulfur concentrations (0, 6, 12, 24 and 60 g/ per experimental unit: pot) using real-time quantitative PCR (RT-qPCR) analysis. For understanding the network interactions between studied genes and other related genes, in silico gene network analysis was constructed to investigate various mechanisms underlying stimulation of A. sativum L. to cope with imposed sulfur. Complicated network including TF-TF, miRNA-TF, and miRNA-TF-gene, was split into sub-networks to have a deeper insight. Analysis of q-RT-PCR data revealed the highest expression in All and SiR genes respectively. To distinguish and select significant pathways in sulfur metabolism, RESNET Plant database of Pathway Studio software v.10 (Elsevier), and other relative data such as chemical reactions, TFs, miRNAs, enzymes, and small molecules were extracted. Complex sub-network exhibited plenty of routes between stress response and sulfate assimilation pathway. Even though Alliinase did not display any connectivity with other stress response genes, it showed binding relation with lectin functional class, as a result of which connected to leucine zipper, exocellulase, peroxidase and ARF functional class indirectly. Integration network of these genes revealed their involvement in various biological processes such as, RNA splicing, stress response, gene silencing by miRNAs, and epigenetic. The findings of this research can be used to extend further research on the garlic metabolic engineering, garlic stress related genes, and also reducing or enhancing the activity of the responsible genes for garlic pungency for health benefits and industry demands.

HetFCM: functional co-module discovery by heterogeneous network co-clustering

Functional molecular module (i.e., gene-miRNA co-modules and gene-miRNA-lncRNA triple-layer modules) analysis can dissect complex regulations underlying etiology or phenotypes. However, current module detection methods lack an appropriate usage and effective model of multi-omics data and cross-layer regulations of heterogeneous molecules, causing the loss of critical genetic information and corrupting the detection performance. In this study, we propose a heterogeneous network co-clustering framework (HetFCM) to detect functional co-modules. HetFCM introduces an attributed heterogeneous network to jointly model interplays and multi-type attributes of different molecules, and applies multiple variational graph autoencoders on the network to generate cross-layer association matrices, then it performs adaptive weighted co-clustering on association matrices and attribute data to identify co-modules of heterogeneous molecules. Empirical study on Human and Maize datasets reveals that HetFCM can find out co-modules characterized with denser topology and more significant functions, which are associated with human breast cancer (subtypes) and maize phenotypes (i.e., lipid storage, drought tolerance and oil content). HetFCM is a useful tool to detect co-modules and can be applied to multi-layer functional modules, yielding novel insights for analyzing molecular mechanisms. We also developed a user-friendly module detection and analysis tool and shared it at http://www.sdu-idea.cn/FMDTool.

Identification of common factors among fibrosarcoma, rhabdomyosarcoma, and osteosarcoma by network analysis

Sarcoma cancers are uncommon malignant tumors, and there are many subgroups, including fibrosarcoma (FS), which mainly affects middle-aged and older adults in deep soft tissues. Rhabdomyosarcoma (RMS), on the other hand, is the most common soft-tissue sarcoma in children and is located in the head and neck area. Osteosarcomas (OS) is the predominant form of primary bone cancer among young adults, primarily resulting from sporadically random mutations. This frequently results in the dissemination of cancer cells to the lungs, commonly known as metastasis. Mesodermal cells are the origin of sarcoma cancers. In this study, a rather radical approach has been applied. Instead of comparing homogenous cancer types, we focus on three main subtypes of sarcoma: fibrosarcoma, rhabdomyosarcoma, and osteosarcoma, and compare their gene expression with normal cell groups to identify the differentially expressed genes (DEGs). Next, by applying protein-protein interaction (PPI) network analysis, we determine the hub genes and crucial factors, such as transcription factors (TFs), affected by these types of cancer. Our findings indicate a modification in a range of pathways associated with cell cycle, extracellular matrix, and DNA repair in these three malignancies. Results showed that fibrosarcoma (FS), rhabdomyosarcoma (RMS), and osteosarcoma (OS) had 653, 1270, and 2823 differentially expressed genes (DEGs), respectively. Interestingly, there were 24 DEGs common to all three types. Network analysis showed that the fibrosarcoma network had two sub-networks identified in FS that contributed to the catabolic process of collagen via the G-protein coupled receptor signaling pathway. The rhabdomyosarcoma network included nine sub-networks associated with cell division, extracellular matrix organization, mRNA splicing via spliceosome, and others. The osteosarcoma network has 13 sub-networks, including mRNA splicing, sister chromatid cohesion, DNA repair, etc. In conclusion, the common DEGs identified in this study have been shown to play significant and multiple roles in various other cancers based on the literature review, indicating their significance.

Integration of circulating microRNAs and transcriptome signatures identifies early-pregnancy biomarkers of preeclampsia

BACKGROUND

MicroRNAs (miRNAs) have been implicated in the pathobiology of preeclampsia, a common hypertensive disorder of pregnancy. In a nested matched case-control cohort within the Vitamin D Antenatal Asthma Reduction Trial (VDAART), we previously identified peripheral blood mRNA signatures related to preeclampsia and vitamin D status (≤30 ng/mL) during gestation from 10 to 18 weeks, using differential expression analysis.

METHODS

Using quantitative PCR arrays, we conducted profiling of circulating miRNAs at 10-18 weeks of gestation in the same VDAART cohort to identify differentially expressed (DE) miRNAs associated with preeclampsia and vitamin D status. For the validation of the expression of circulating miRNA signatures in the placenta, the HTR-8/SVneo trophoblast cell line was used. Targets of circulating miRNA signatures in the preeclampsia mRNA signatures were identified by consensus ranking of miRNA-target prediction scores from four sources. The connected component of target signatures was identified by mapping to the protein-protein interaction (PPI) network and hub targets were determined. As experimental validation, we examined the gene and protein expression of IGF1R, one of the key hub genes, as a target of the DE miRNA, miR-182-5p, in response to a miR-182-5p mimic in HTR-8/SVneo cells.

RESULTS

Pregnant women with preeclampsia had 16 circulating DE miRNAs relative to normal pregnancy controls that were also DE under vitamin D insufficiency (9/16 = 56% upregulated, FDR < .05). Thirteen miRNAs (13/16 = 81.3%) were detected in HTR-8/SVneo cells. Overall, 16 DE miRNAs had 122 targets, of which 87 were unique. Network analysis demonstrated that the 32 targets of DE miRNA signatures created a connected subnetwork in the preeclampsia module with CXCL8, CXCL10, CD274, MMP9 and IGF1R having the highest connectivity and centrality degree. In an in vitro validation experiment, the introduction of an hsa-miR-182-5p mimic resulted in significant reduction of its target IGF1R gene and protein expression within HTR-8/SVneo cells.

CONCLUSIONS

The integration of the circulating DE miRNA and mRNA signatures associated preeclampsia added additional insights into the subclinical molecular signature of preeclampsia. Our systems and network biology approach revealed several biological pathways, including IGF-1, that may play a role in the early pathophysiology of preeclampsia. These pathways and signatures also denote potential biomarkers for the early stages of preeclampsia and suggest possible preventive measures.

scapGNN: A graph neural network-based framework for active pathway and gene module inference from single-cell multi-omics data

Although advances in single-cell technologies have enabled the characterization of multiple omics profiles in individual cells, extracting functional and mechanistic insights from such information remains a major challenge. Here, we present scapGNN, a graph neural network (GNN)-based framework that creatively transforms sparse single-cell profile data into the stable gene-cell association network for inferring single-cell pathway activity scores and identifying cell phenotype-associated gene modules from single-cell multi-omics data. Systematic benchmarking demonstrated that scapGNN was more accurate, robust, and scalable than state-of-the-art methods in various downstream single-cell analyses such as cell denoising, batch effect removal, cell clustering, cell trajectory inference, and pathway or gene module identification. scapGNN was developed as a systematic R package that can be flexibly extended and enhanced for existing analysis processes. It provides a new analytical platform for studying single cells at the pathway and network levels.

A novel approach for predicting upstream regulators (PURE) that affect gene expression

External factors such as exposure to a chemical, drug, or toxicant (CDT), or conversely, the lack of certain chemicals can cause many diseases. The ability to identify such causal CDTs based on changes in the gene expression profile is extremely important in many studies. Furthermore, the ability to correctly infer CDTs that can revert the gene expression changes induced by a given disease phenotype is a crucial step in drug repurposing. We present an approach for Predicting Upstream REgulators (PURE) designed to tackle this challenge. PURE can correctly infer a CDT from the measured expression changes in a given phenotype, as well as correctly identify drugs that could revert disease-induced gene expression changes. We compared the proposed approach with four classical approaches as well as with the causal analysis used in Ingenuity Pathway Analysis (IPA) on 16 data sets (1 rat, 5 mouse, and 10 human data sets), involving 8 chemicals or drugs. We assessed the results based on the ability to correctly identify the CDT as indicated by its rank. We also considered the number of false positives, i.e. CDTs other than the correct CDT that were reported to be significant by each method. The proposed approach performed best in 11 out of the 16 experiments, reporting the correct CDT at the very top 7 times. IPA was the second best, reporting the correct CDT at the top 5 times, but was unable to identify the correct CDT at all in 5 out of the 16 experiments. The validation results showed that our approach, PURE, outperformed some of the most popular methods in the field. PURE could effectively infer the true CDTs responsible for the observed gene expression changes and could also be useful in drug repurposing applications.

Gene expression profile analysis unravelled the systems level association of renal cell carcinoma with diabetic nephropathy and Matrix-metalloproteinase-9 as a potential therapeutic target

Type 2 diabetes (T2D) and cancer share many common risk factors. However, the potential biological link that connects the two at the molecular level is still unclear. The experimental evidence suggests that several genes and their pathways may be involved in developing cancerous conditions associated with diabetes. In this study, we identified the protein-protein interaction (PPI) networks and the hub protein(s) that interlink T2D and cancer using genome-scale differential gene expression profiles. Further, the PPI network of AMP-activated protein kinase (AMPK) in cancer was analyzed to explore novel insights into the molecular association between the two conditions. The densely connected regions were analyzed by constructing the backbone and subnetworks with key nodes and shortest pathways, respectively. The PPI network studies identified Matrix-metalloproteinase-9 (MMP-9) as a hub protein playing a vital role in glomerulonephritis tubular diseases and some genetic kidney diseases. MMP-9 was also associated with different growth factors, like tumor necrosis factor (TNF-α), transforming growth factor 1 (TGF-1), and pathways like chemokine signaling, NOD-like receptor signaling, etc. Further, the molecular docking and molecular dynamic simulation studies supported the druggability of MMP-9, suggesting it as a potential therapeutic target in treating renal cell carcinoma linked with diabetic kidney disease.Communicated by Ramaswamy H. Sarma.

Identification of shared biological features in four different lung cell lines infected with SARS-CoV-2 virus through RNA-seq analysis

The COVID-19 pandemic caused by SARS-CoV-2 has resulted in millions of confirmed cases and deaths worldwide. Understanding the biological mechanisms of SARS-CoV-2 infection is crucial for the development of effective therapies. This study conducts differential expression (DE) analysis, pathway analysis, and differential network (DN) analysis on RNA-seq data of four lung cell lines, NHBE, A549, A549.ACE2, and Calu3, to identify their common and unique biological features in response to SARS-CoV-2 infection. DE analysis shows that cell line A549.ACE2 has the highest number of DE genes, while cell line NHBE has the lowest. Among the DE genes identified for the four cell lines, 12 genes are overlapped, associated with various health conditions. The most significant signaling pathways varied among the four cell lines. Only one pathway, "cytokine-cytokine receptor interaction", is found to be significant among all four cell lines and is related to inflammation and immune response. The DN analysis reveals considerable variation in the differential connectivity of the most significant pathway shared among the four lung cell lines. These findings help to elucidate the mechanisms of SARS-CoV-2 infection and potential therapeutic targets.

Systems biology and metabolic modeling for cultivated meat: A promising approach for cell culture media optimization and cost reduction

The cultivated meat industry, also known as cell-based meat, cultured meat, lab-grown meat, or meat alternatives, is a growing field that aims to generate animal tissues ex-vivo in a cost-effective manner that achieves price parity with traditional agricultural products. However, cell culture media costs account for 55%-90% of production costs. To address this issue, efforts are aimed at optimizing media composition. Systems biology-driven approaches have been successfully used to improve the biomass and productivity of multiple bioproduction platforms, like Chinese hamster ovary cells, by accelerating the development of cell line-specific media and reducing research and development and production costs related to cell media and its optimization. In this review, we summarize systems biology modeling approaches, methods for cell culture media and bioprocess optimization, and metabolic studies done in animals of interest to the cultivated meat industry. More importantly, we identify current gaps in knowledge that prevent the identification of metabolic bottlenecks. These include the lack of genome-scale metabolic models for some species (pigs and ducks), a lack of accurate biomass composition studies for different growth conditions, and 13 C-metabolic flux analysis (MFA) studies for many of the species of interest for the cultivated meat industry (only shrimp and duck cells have been subjected to 13 C-MFA). We also highlight the importance of characterizing the metabolic requirements of cells at the organism, breed, and cell line-specific levels, and we outline future steps that this nascent field needs to take to achieve price parity and production efficiency similar to those of other bioproduction platforms. Practical Application: Our article summarizes systems biology techniques for cell culture media design and bioprocess optimization, which may be used to significantly reduce cell-based meat production costs. We also present the results of experimental studies done on some of the species of interest to the cultivated meat industry and highlight why modeling approaches are required for multiple species, cell-types, and cell lines.

Gemini: memory-efficient integration of hundreds of gene networks with high-order pooling

MOTIVATION

The exponential growth of genomic sequencing data has created ever-expanding repositories of gene networks. Unsupervised network integration methods are critical to learn informative representations for each gene, which are later used as features for downstream applications. However, these network integration methods must be scalable to account for the increasing number of networks and robust to an uneven distribution of network types within hundreds of gene networks.

RESULTS

To address these needs, we present Gemini, a novel network integration method that uses memory-efficient high-order pooling to represent and weight each network according to its uniqueness. Gemini then mitigates the uneven network distribution through mixing up existing networks to create many new networks. We find that Gemini leads to more than a 10% improvement in F1 score, 15% improvement in micro-AUPRC, and 63% improvement in macro-AUPRC for human protein function prediction by integrating hundreds of networks from BioGRID, and that Gemini's performance significantly improves when more networks are added to the input network collection, while Mashup and BIONIC embeddings' performance deteriorates. Gemini thereby enables memory-efficient and informative network integration for large gene networks and can be used to massively integrate and analyze networks in other domains.

AVAILABILITY AND IMPLEMENTATION

Gemini can be accessed at: https://github.com/MinxZ/Gemini.

Genetic epidemiology of resistance to M. tuberculosis Infection: importance of study design and recent findings

Resistance to M. tuberculosis, often referred to as "RSTR" in the literature, is being increasingly studied because of its potential relevance as a clinical outcome in vaccine studies. This review starts by addressing the importance of epidemiological characterization of this phenotype, and ongoing challenges in that characterization. Then, this review summarizes the extant genetic and genomic studies of this phenotype, including heritability studies, candidate gene studies, and genome-wide association studies, as well as whole transcriptome studies. Findings from recent studies that used longitudinal characterization of the RSTR phenotype are compared to those using a cross-sectional definition, and the challenges of using tuberculin skin test and interferon-gamma release assay are discussed. Finally, future directions are proposed. Since this is a rapidly evolving area of public health significance, this review will help frame future research questions and study designs.

SEMtree: tree-based structure learning methods with structural equation models

MOTIVATION

With the exponential growth of expression and protein-protein interaction (PPI) data, the identification of functional modules in PPI networks that show striking changes in molecular activity or phenotypic signatures becomes of particular interest to reveal process-specific information that is correlated with cellular or disease states. This requires both the identification of network nodes with reliability scores and the availability of an efficient technique to locate the network regions with the highest scores. In the literature, a number of heuristic methods have been suggested. We propose SEMtree(), a set of tree-based structure discovery algorithms, combining graph and statistically interpretable parameters together with a user-friendly R package based on structural equation models framework.

RESULTS

Condition-specific changes from differential expression and gene-gene co-expression are recovered with statistical testing of node, directed edge, and directed path difference between groups. In the end, from a list of seed (i.e. disease) genes or gene P-values, the perturbed modules with undirected edges are generated with five state-of-the-art active subnetwork detection methods. The latter are supplied to causal additive trees based on Chu-Liu-Edmonds' algorithm (Chow and Liu, Approximating discrete probability distributions with dependence trees. IEEE Trans Inform Theory 1968;14:462-7) in SEMtree() to be converted in directed trees. This conversion allows to compare the methods in terms of directed active subnetworks. We applied SEMtree() to both Coronavirus disease (COVID-19) RNA-seq dataset (GEO accession: GSE172114) and simulated datasets with various differential expression patterns. Compared to existing methods, SEMtree() is able to capture biologically relevant subnetworks with simple visualization of directed paths, good perturbation extraction, and classifier performance.

AVAILABILITY AND IMPLEMENTATION

SEMtree() function is implemented in the R package SEMgraph, easily available at https://CRAN.R-project.org/package=SEMgraph.

Rewiring of miRNA-mRNA bipartite co-expression network as a novel way to understand the prostate cancer related players

The differential expression and direct targeting of mRNA by miRNA are two main logics of the traditional approach to constructing the miRNA-mRNA network. This approach, could be led to the loss of considerable information and some challenges of direct targeting. To avoid these problems, we analyzed the rewiring network and constructed two miRNA-mRNA expression bipartite networks for both normal and primary prostate cancer tissue obtained from PRAD-TCGA. We then calculated beta-coefficient of the regression-model when miR was dependent and mRNA independent for each miR and mRNA and separately in both networks. We defined the rewired edges as a significant change in the regression coefficient between normal and cancer states. The rewired nodes through multinomial distribution were defined and network from rewired edges and nodes was analyzed and enriched. Of the 306 rewired edges, 112(37%) were new, 123(40%) were lost, 44(14%) were strengthened, and 27(9%) weakened connections were discovered. The highest centrality of 106 rewired mRNAs belonged to PGM5, BOD1L1, C1S, SEPG, TMEFF2, and CSNK2A1. The highest centrality of 68 rewired miRs belonged to miR-181d, miR-4677, miR-4662a, miR-9.3, and miR-1301. SMAD and beta-catenin binding were enriched as molecular functions. The regulation was a frequently repeated concept in the biological process. Our rewiring analysis highlighted the impact of β-catenin and SMAD signaling as also some transcript factors like TGFB1I1 in prostate cancer progression. Altogether, we developed a miRNA-mRNA co-expression bipartite network to identify the hidden aspects of the prostate cancer mechanism, which traditional analysis -like differential expression- was not detect it.

Mena regulates nesprin-2 to control actin-nuclear lamina associations, trans-nuclear membrane signalling and gene expression

Interactions between cells and the extracellular matrix, mediated by integrin adhesion complexes, play key roles in fundamental cellular processes, including the sensing and transduction of mechanical cues. Here, we investigate systems-level changes in the integrin adhesome in patient-derived cutaneous squamous cell carcinoma cells and identify the actin regulatory protein Mena as a key node in the adhesion complex network. Mena is connected within a subnetwork of actin-binding proteins to the LINC complex component nesprin-2, with which it interacts and co-localises at the nuclear envelope. Moreover, Mena potentiates the interactions of nesprin-2 with the actin cytoskeleton and the nuclear lamina. CRISPR-mediated Mena depletion causes altered nuclear morphology, reduces tyrosine phosphorylation of the nuclear membrane protein emerin and downregulates expression of the immunomodulatory gene PTX3 via the recruitment of its enhancer to the nuclear periphery. We uncover an unexpected role for Mena at the nuclear membrane, where it controls nuclear architecture, chromatin repositioning and gene expression. Our findings identify an adhesion protein that regulates gene transcription via direct signalling across the nuclear envelope.

Gene Functional Networks from Time Expression Profiles: A Constructive Approach Demonstrated in Chili Pepper (Capsicum annuum L.)

Gene co-expression networks are powerful tools to understand functional interactions between genes. However, large co-expression networks are difficult to interpret and do not guarantee that the relations found will be true for different genotypes. Statistically verified time expression profiles give information about significant changes in expressions through time, and genes with highly correlated time expression profiles, which are annotated in the same biological process, are likely to be functionally connected. A method to obtain robust networks of functionally related genes will be useful to understand the complexity of the transcriptome, leading to biologically relevant insights. We present an algorithm to construct gene functional networks for genes annotated in a given biological process or other aspects of interest. We assume that there are genome-wide time expression profiles for a set of representative genotypes of the species of interest. The method is based on the correlation of time expression profiles, bound by a set of thresholds that assure both, a given false discovery rate, and the discard of correlation outliers. The novelty of the method consists in that a gene expression relation must be repeatedly found in a given set of independent genotypes to be considered valid. This automatically discards relations particular to specific genotypes, assuring a network robustness, which can be set a priori. Additionally, we present an algorithm to find transcription factors candidates for regulating hub genes within a network. The algorithms are demonstrated with data from a large experiment studying gene expression during the development of the fruit in a diverse set of chili pepper genotypes. The algorithm is implemented and demonstrated in a new version of the publicly available R package "Salsa" (version 1.0).

Genome-wide characterization of the common bean kinome: Catalog and insights into expression patterns and genetic organization

The protein kinase (PK) superfamily is one of the largest superfamilies in plants and is the core regulator of cellular signaling. Even considering this substantial importance, the kinome of common bean (Phaseolus vulgaris) has not been profiled yet. Here, we identified and characterised the complete set of kinases of common bean, performing an in-depth investigation with phylogenetic analyses and measurements of gene distribution, structural organization, protein properties, and expression patterns over a large set of RNA-Sequencing data. Being composed of 1,203 PKs distributed across all P. vulgaris chromosomes, this set represents 3.25% of all predicted proteins for the species. These PKs could be classified into 20 groups and 119 subfamilies, with a more pronounced abundance of subfamilies belonging to the receptor-like kinase (RLK)-Pelle group. In addition to provide a vast and rich reservoir of data, our study supplied insights into the compositional similarities between PK subfamilies, their evolutionary divergences, highly variable functional profile, structural diversity, and expression patterns, modeled with coexpression networks for investigating putative interactions associated with stress response.

Evolution of gene regulatory network of C4 photosynthesis in the genus Flaveria reveals the evolutionary status of C3-C4 intermediate species

C₄ photosynthesis evolved from ancestral C₃ photosynthesis by recruiting pre-existing genes to fulfill new functions. The enzymes and transporters required for the C₄ metabolic pathway have been intensively studied and well documented; however, the transcription factors (TFs) that regulate these C₄ metabolic genes are not yet well understood. In particular, how the TF regulatory network of C₄ metabolic genes was rewired during the evolutionary process is unclear. Here, we constructed gene regulatory networks (GRNs) for four closely evolutionarily related species from the genus Flaveria, which represent four different evolutionary stages of C₄ photosynthesis: C₃ (F. robusta), type I C₃-C₄ (F. sonorensis), type II C₃-C₄ (F. ramosissima), and C₄ (F. trinervia). Our results show that more than half of the co-regulatory relationships between TFs and core C₄ metabolic genes are species specific. The counterparts of the C₄ genes in C₃ species were already co-regulated with photosynthesis-related genes, whereas the required TFs for C₄ photosynthesis were recruited later. The TFs involved in C₄ photosynthesis were widely recruited in the type I C₃-C₄ species; nevertheless, type II C₃-C₄ species showed a divergent GRN from C₄ species. In line with these findings, a ¹³CO₂ pulse-labeling experiment showed that the CO₂ initially fixed into C₄ acid was not directly released to the Calvin-Benson-Bassham cycle in the type II C₃-C₄ species. Therefore, our study uncovered dynamic changes in C₄ genes and TF co-regulation during the evolutionary process; furthermore, we showed that the metabolic pathway of the type II C₃-C₄ species F. ramosissima represents an alternative evolutionary solution to the ammonia imbalance in C₃-C₄ intermediate species.

Systems Pharmacology and Pharmacokinetics Strategy to Decode Bioactive Ingredients and Molecular Mechanisms from Zingiber officinale as Phyto-therapeutics against Neurological Diseases

BACKGROUND

The bioactive constituents from Zingiber officinale (Z. officinale) have shown a positive effect on neurodegenerative diseases like Alzheimer's disease (AD), which manifests as progressive memory loss and cognitive impairment.

OBJECTIVE

This study investigates the binding ability and the pharmaco-therapeutic potential of Z. officinale with AD disease targets by molecular docking and molecular dynamic (MD) simulation approaches.

METHODS

By coupling enormous available phytochemical data and advanced computational technologies, the possible molecular mechanism of action of these bioactive compounds was deciphered by evaluating phytochemicals, target fishing, and network biological analysis.

RESULTS

As a result, 175 bioactive compounds and 264 human target proteins were identified. The gene ontology and Kyoto Encyclopaedia of Genes and Genomes pathway enrichment analysis and molecular docking were used to predict the basis of vital bioactive compounds and biomolecular mechanisms involved in the treatment of AD. Amongst selected bioactive compounds, 10- Gingerdione and 1-dehydro-[8]-gingerdione exhibited significant anti-neurological properties against AD targeting amyloid precursor protein with docking energy of -6.0 and -5.6, respectively.

CONCLUSION

This study suggests that 10-Gingerdione and 1-dehydro-[8]-gingerdione strongly modulates the anti-neurological activity and are associated with pathological features like amyloid-β plaques and hyperphosphorylated tau protein are found to be critically regulated by these two target proteins. This comprehensive analysis provides a clue for further investigation of these natural compounds' inhibitory activity in drug discovery for AD treatment.

Interdisciplinary Overview of Lipopeptide and Protein-Containing Biosurfactants

Biosurfactants are amphipathic molecules capable of lowering interfacial and superficial tensions. Produced by living organisms, these compounds act the same as chemical surfactants but with a series of improvements, the most notable being biodegradability. Biosurfactants have a wide diversity of categories. Within these, lipopeptides are some of the more abundant and widely known. Protein-containing biosurfactants are much less studied and could be an interesting and valuable alternative. The harsh temperature, pH, and salinity conditions that target organisms can sustain need to be understood for better implementation. Here, we will explore biotechnological applications via lipopeptide and protein-containing biosurfactants. Also, we discuss their natural role and the organisms that produce them, taking a glimpse into the possibilities of research via meta-omics and machine learning.

An Approach for Systems-Level Understanding of Prostate Cancer from High-Throughput Data Integration to Pathway Modeling and Simulation

To understand complex diseases, high-throughput data are generated at large and multiple levels. However, extracting meaningful information from large datasets for comprehensive understanding of cell phenotypes and disease pathophysiology remains a major challenge. Despite tremendous advances in understanding molecular mechanisms of cancer and its progression, current knowledge appears discrete and fragmented. In order to render this wealth of data more integrated and thus informative, we have developed a GECIP toolbox to investigate the crosstalk and the responsible genes'/proteins' connectivity of enriched pathways from gene expression data. To implement this toolbox, we used mainly gene expression datasets of prostate cancer, and the three datasets were GSE17951, GSE8218, and GSE1431. The raw samples were processed for normalization, prediction of differentially expressed genes, and the prediction of enriched pathways for the differentially expressed genes. The enriched pathways have been processed for crosstalk degree calculations for which number connections per gene, the frequency of genes in the pathways, sharing frequency, and the connectivity have been used. For network prediction, protein-protein interaction network database FunCoup2.0 was used, and cytoscape software was used for the network visualization. In our results, we found that there were enriched pathways 27, 45, and 22 for GSE17951, GSE8218, and GSE1431, respectively, and 11 pathways in common between all of them. From the crosstalk results, we observe that focal adhesion and PI3K pathways, both experimentally proven central for cellular output upon perturbation of numerous individual/distinct signaling pathways, displayed highest crosstalk degree. Moreover, we also observe that there were more critical pathways which appear to be highly significant, and these pathways are HIF1a, hippo, AMPK, and Ras. In terms of the pathways' components, GSK3B, YWHAE, HIF1A, ATP1A3, and PRKCA are shared between the aforementioned pathways and have higher connectivity with the pathways and the other pathway components. Finally, we conclude that the focal adhesion and PI3K pathways are the most critical pathways, and since for many other pathways, high-rank enrichment did not translate to high crosstalk degree, the global impact of one pathway on others appears distinct from enrichment.

Exploring key molecular signatures of immune responses and pathways associated with tuberculosis in comorbid diabetes mellitus: a systems biology approach

Background

Comorbid type 2 diabetes mellitus (T2DM) increases the risk for tuberculosis (TB) and its associated complications, although the pathological connections between T2DM and TB are unknown. The current research aims to identify shared molecular gene signatures and pathways that affirm the epidemiological association of T2DM and TB and afford clues on mechanistic basis of their association through integrative systems biology and bioinformatics approaches. Earlier research has found specific molecular markers linked to T2DM and TB, but, despite their importance, only offered a limited understanding of the genesis of this comorbidity. Our investigation used a network medicine method to find possible T2DM-TB molecular mediators.

Results

Functional annotation clustering, interaction networks, network cluster analysis, and network topology were part of our systematic investigation of T2DM-TB linked with 1603 differentially expressed genes (DEGs). The functional enrichment and gene interaction network analysis emphasized the importance of cytokine/chemokine signalling, T cell receptor signalling route, NF-kappa B signalling pathway and Jak-STAT signalling system. Furthermore, network analysis revealed significant DEGs such as ITGAM and STAT1, which may be necessary for T2DM-TB immune responses. Furthermore, these two genes are modulators in clusters C4 and C5, abundant in cytokine/chemokine signalling and Jak-STAT signalling pathways.

Conclusions

Our analyses highlight the role of ITGAM and STAT1 in T2DM-TB-associated pathways and advances our knowledge of the genetic processes driving this comorbidity.

Selecting the right therapeutic target for kidney disease

Kidney disease is a complex disease with several different etiologies and underlying associated pathophysiology. This is reflected by the lack of effective treatment therapies in chronic kidney disease (CKD) that stop disease progression. However, novel strategies, recent scientific breakthroughs, and technological advances have revealed new possibilities for finding novel disease drivers in CKD. This review describes some of the latest advances in the field and brings them together in a more holistic framework as applied to identification and validation of disease drivers in CKD. It uses high-resolution 'patient-centric' omics data sets, advanced in silico tools (systems biology, connectivity mapping, and machine learning) and 'state-of-the-art' experimental systems (complex 3D systems in vitro, CRISPR gene editing, and various model biological systems in vivo). Application of such a framework is expected to increase the likelihood of successful identification of novel drug candidates based on strong human target validation and a better scientific understanding of underlying mechanisms.

Validating methods for testing natural molecules on molecular pathways of interest in silico and in vitro

Differentially expressed genes can serve as drug targets and are used to predict drug response and disease progression. In silico drug analysis based on the expression of these genetic biomarkers allows the detection of putative therapeutic agents, which could be used to reverse a pathological gene expression signature. Indeed, a set of bioinformatics tools can increase the accuracy of drug discovery, helping in biomarker identification. Once a drug target is identified, in vitro cell line models of disease are used to evaluate and validate the therapeutic potential of putative drugs and novel natural molecules. This study describes the development of efficacious PCR primers that can be used to identify gene expression of specific genetic pathways, which can lead to the identification of natural molecules as therapeutic agents in specific molecular pathways. For this study, genes involved in health conditions and processes were considered. In particular, the expression of genes involved in obesity, xenobiotics metabolism, endocannabinoid pathway, leukotriene B4 metabolism and signaling, inflammation, endocytosis, hypoxia, lifespan, and neurotrophins were evaluated. Exploiting the expression of specific genes in different cell lines can be useful in in vitro to evaluate the therapeutic effects of small natural molecules.

Applying differential network analysis to longitudinal gene expression in response to perturbations

Differential Network (DN) analysis is a method that has long been used to interpret changes in gene expression data and provide biological insights. The method identifies the rewiring of gene networks in response to external perturbations. Our study applies the DN method to the analysis of RNA-sequencing (RNA-seq) time series datasets. We focus on expression changes: (i) in saliva of a human subject after pneumococcal vaccination (PPSV23) and (ii) in primary B cells treated ex vivo with a monoclonal antibody drug (Rituximab). The DN method enabled us to identify the activation of biological pathways consistent with the mechanisms of action of the PPSV23 vaccine and target pathways of Rituximab. The community detection algorithm on the DN revealed clusters of genes characterized by collective temporal behavior. All saliva and some B cell DN communities showed characteristic time signatures, outlining a chronological order in pathway activation in response to the perturbation. Moreover, we identified early and delayed responses within network modules in the saliva dataset and three temporal patterns in the B cell data.

In search of a Drosophila core cellular network with single-cell transcriptome data

Along with specialized functions, cells of multicellular organisms also perform essential functions common to most if not all cells. Whether diverse cells do this by using the same set of genes, interacting in a fixed coordinated fashion to execute essential functions, or a subset of genes specific to certain cells, remains a central question in biology. Here, we focus on gene coexpression to search for a core cellular network across a whole organism. Single-cell RNA-sequencing measures gene expression of individual cells, enabling researchers to discover gene expression patterns that contribute to the diversity of cell functions. Current efforts to study cellular functions focus primarily on identifying differentially expressed genes across cells. However, patterns of coexpression between genes are probably more indicative of biological processes than are the expression of individual genes. We constructed cell-type-specific gene coexpression networks using single-cell transcriptome datasets covering diverse cell types from the fruit fly, Drosophila melanogaster. We detected a set of highly coordinated genes preserved across cell types and present this as the best estimate of a core cellular network. This core is very small compared with cell-type-specific gene coexpression networks and shows dense connectivity. Gene members of this core tend to be ancient genes and are enriched for those encoding ribosomal proteins. Overall, we find evidence for a core cellular network in diverse cell types of the fruit fly. The topological, structural, functional, and evolutionary properties of this core indicate that it accounts for only a minority of essential functions.

Mendelian randomization of circulating proteome identifies actionable targets in heart failure

Background

Heart failure (HF) is a prevalent cause of mortality and morbidity. The molecular drivers of HF are still largely unknown.

Results

We aimed to identify circulating proteins causally associated with HF by leveraging genome-wide genetic association data for HF including 47,309 cases and 930,014 controls. We performed two-sample Mendelian randomization (MR) with multiple cis instruments as well as network and enrichment analysis using data from blood protein quantitative trait loci (pQTL) (2,965 blood proteins) measured in 3,301 individuals. Nineteen blood proteins were causally associated with HF, were not subject to reverse causality and were enriched in ligand-receptor and glycosylation molecules. Network pathway analysis of the blood proteins showed enrichment in NF-kappa B, TGF beta, lipid in atherosclerosis and fluid shear stress. Cross-phenotype analysis of HF identified genetic overlap with cardiovascular drugs, myocardial infarction, parental longevity and low-density cholesterol. Multi-trait MR identified causal associations between HF-associated blood proteins and cardiovascular outcomes. Multivariable MR showed that association of BAG3, MIF and APOA5 with HF were mediated by the blood pressure and coronary artery disease. According to the directional effect and biological action, 7 blood proteins are targets of existing drugs or are tractable for the development of novel therapeutics. Among the pathways, sialyl Lewis x and the activin type II receptor are potential druggable candidates.

Conclusions

Integrative MR analyses of the blood proteins identified causally-associated proteins with HF and revealed pleiotropy of the blood proteome with cardiovascular risk factors. Some of the proteins or pathway related mechanisms could be targeted as novel treatment approach in HF.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-022-08811-2.

Adrenergic signaling controls early transcriptional programs during CD8+ T cell responses to viral infection

Norepinephrine is a key sympathetic neurotransmitter, which acts to suppress CD8 + T cell cytokine secretion and lytic activity by signaling through the β2-adrenergic receptor (ADRB2). Although ADRB2 signaling is considered generally immunosuppressive, its role in regulating the differentiation of effector T cells in response to infection has not been investigated. Using an adoptive transfer approach, we compared the expansion and differentiation of wild type (WT) to Adrb2-/- CD8 + T cells throughout the primary response to vesicular stomatitis virus (VSV) infection in vivo. We measured the dynamic changes in transcriptome profiles of antigen-specific CD8 + T cells as they responded to VSV. Within the first 7 days of infection, WT cells out-paced the expansion of Adrb2-/- cells, which correlated with reduced expression of IL-2 and the IL-2Rα in the absence of ADRB2. RNASeq analysis identified over 300 differentially expressed genes that were both temporally regulated following infection and selectively regulated in WT vs Adrb2-/- cells. These genes contributed to major transcriptional pathways including cytokine receptor activation, signaling in cancer, immune deficiency, and neurotransmitter pathways. By parsing genes within groups that were either induced or repressed over time in response to infection, we identified three main branches of genes that were differentially regulated by the ADRB2. These gene sets were predicted to be regulated by specific transcription factors involved in effector T cell development, such as Tbx21 and Eomes. Collectively, these data demonstrate a significant role for ADRB2 signaling in regulating key transcriptional pathways during CD8 + T cells responses to infection that may dramatically impact their functional capabilities and downstream memory cell development.

A Path-Based Analysis of Infected Cell Line and COVID-19 Patient Transcriptome Reveals Novel Potential Targets and Drugs Against SARS-CoV-2

Most transcriptomic studies of SARS-CoV-2 infection have focused on differentially expressed genes, which do not necessarily reveal the genes mediating the transcriptomic changes. In contrast, exploiting curated biological network, our PathExt tool identifies central genes from the differentially active paths mediating global transcriptomic response. Here we apply PathExt to multiple cell line infection models of SARS-CoV-2 and other viruses, as well as to COVID-19 patient-derived PBMCs. The central genes mediating SARS-CoV-2 response in cell lines were uniquely enriched for ATP metabolic process, G1/S transition, leukocyte activation and migration. In contrast, PBMC response reveals dysregulated cell-cycle processes. In PBMC, the most frequently central genes are associated with COVID-19 severity. Importantly, relative to differential genes, PathExt-identified genes show greater concordance with several benchmark anti-COVID-19 target gene sets. We propose six novel anti-SARS-CoV-2 targets ADCY2, ADSL, OCRL, TIAM1, PBK, and BUB1, and potential drugs targeting these genes, such as Bemcentinib, Phthalocyanine, and Conivaptan.

Systems seed biology to understand and manipulate rice grain quality and nutrition

Rice is one of the most essential crops since it meets the calorific needs of 3 billion people around the world. Rice seed development initiates upon fertilization, leading to the establishment of two distinct filial tissues, the endosperm and embryo, which accumulate distinct seed storage products, such as starch, storage proteins, and lipids. A range of systems biology tools deployed in dissecting the spatiotemporal dynamics of transcriptome data, methylation, and small RNA based regulation operative during seed development, influencing the accumulation of storage products was reviewed. Studies of other model systems are also considered due to the limited information on the rice transcriptome. This review highlights key genes identified through a holistic view of systems biology targeted to modify biochemical composition and influence rice grain quality and nutritional value with the target of improving rice as a functional food.

Network-Based Methods for Approaching Human Pathologies from a Phenotypic Point of View

Network and systemic approaches to studying human pathologies are helping us to gain insight into the molecular mechanisms of and potential therapeutic interventions for human diseases, especially for complex diseases where large numbers of genes are involved. The complex human pathological landscape is traditionally partitioned into discrete “diseases”; however, that partition is sometimes problematic, as diseases are highly heterogeneous and can differ greatly from one patient to another. Moreover, for many pathological states, the set of symptoms (phenotypes) manifested by the patient is not enough to diagnose a particular disease. On the contrary, phenotypes, by definition, are directly observable and can be closer to the molecular basis of the pathology. These clinical phenotypes are also important for personalised medicine, as they can help stratify patients and design personalised interventions. For these reasons, network and systemic approaches to pathologies are gradually incorporating phenotypic information. This review covers the current landscape of phenotype-centred network approaches to study different aspects of human diseases.

Benefits and Challenges of Pre-clustered Network-Based Pathway Analysis

Functional analysis of gene sets derived from experiments is typically done by pathway annotation. Although many algorithms exist for analyzing the association between a gene set and a pathway, an issue which is generally ignored is that gene sets often represent multiple pathways. In such cases an association to a pathway is weakened by the presence of genes associated with other pathways. A way to counteract this is to cluster the gene set into more homogenous parts before performing pathway analysis on each module. We explored whether network-based pre-clustering of a query gene set can improve pathway analysis. The methods MCL, Infomap, and MGclus were used to cluster the gene set projected onto the FunCoup network. We characterized how well these methods are able to detect individual pathways in multi-pathway gene sets, and applied each of the clustering methods in combination with four pathway analysis methods: Gene Enrichment Analysis, BinoX, NEAT, and ANUBIX. Using benchmarks constructed from the KEGG pathway database we found that clustering can be beneficial by increasing the sensitivity of pathway analysis methods and by providing deeper insights of biological mechanisms related to the phenotype under study. However, keeping a high specificity is a challenge. For ANUBIX, clustering caused a minor loss of specificity, while for BinoX and NEAT it caused an unacceptable loss of specificity. GEA had very low sensitivity both before and after clustering. The choice of clustering method only had a minor effect on the results. We show examples of this approach and conclude that clustering can improve overall pathway annotation performance, but should only be used if the used enrichment method has a low false positive rate.

De novo identification of maximally deregulated subnetworks based on multi-omics data with DeRegNet

Background

With a growing amount of (multi-)omics data being available, the extraction of knowledge from these datasets is still a difficult problem. Classical enrichment-style analyses require predefined pathways or gene sets that are tested for significant deregulation to assess whether the pathway is functionally involved in the biological process under study. De novo identification of these pathways can reduce the bias inherent in predefined pathways or gene sets. At the same time, the definition and efficient identification of these pathways de novo from large biological networks is a challenging problem.

Results

We present a novel algorithm, DeRegNet, for the identification of maximally deregulated subnetworks on directed graphs based on deregulation scores derived from (multi-)omics data. DeRegNet can be interpreted as maximum likelihood estimation given a certain probabilistic model for de-novo subgraph identification. We use fractional integer programming to solve the resulting combinatorial optimization problem. We can show that the approach outperforms related algorithms on simulated data with known ground truths. On a publicly available liver cancer dataset we can show that DeRegNet can identify biologically meaningful subgraphs suitable for patient stratification. DeRegNet can also be used to find explicitly multi-omics subgraphs which we demonstrate by presenting subgraphs with consistent methylation-transcription patterns. DeRegNet is freely available as open-source software.

Conclusion

The proposed algorithmic framework and its available implementation can serve as a valuable heuristic hypothesis generation tool contextualizing omics data within biomolecular networks.