Single-cell transcriptomic reveals network topology changes of cancer at the individual level

Network biology facilitates a better understanding of complex diseases. Single-sample networks retain individual information and have the potential to distinguish disease status. Previous studies mainly used bulk RNA sequencing data to construct single-sample networks, but different cell types in the tissue microenvironment perform significantly different functions. In this study, we investigated whether network topology features of cell-type-specific networks varied in different pathological states at the individual level. Protein-protein interaction network (PPI) and co-expression network of cancer and ulcerative colitis were established using four publicly single-cell RNA sequencing (scRNA-seq) datasets. We analyzed cell-cell interactions of epithelial cells and immune cells using CellChat R package. Network topology changes between normal tissues and pathological tissues were analyzed using Cytoscape software and QUACN R package. Results showed cell-cell interactions of epithelial cells were enhanced in carcinoma and adenoma. The average number of neighbors and graphindex of co-expression network increased in epithelial cells of adenoma, carcinoma and paracancer compared with normal tissues. The co-expression network density of T cells in tumors was significantly higher than that in normal tissues. The co-expression network complexity of epithelial cells in the benign tissues was associated with the grade group of paired tumors. This study suggests topological properties of cell-type-specific individual network vary in different pathological states, providing an insight into understanding complex diseases.

Uncovering critical transitions and molecule mechanisms in disease progressions using Gaussian graphical optimal transport

Understanding disease progression is crucial for detecting critical transitions and finding trigger molecules, facilitating early diagnosis interventions. However, the high dimensionality of data and the lack of aligned samples across disease stages have posed challenges in addressing these tasks. We present a computational framework, Gaussian Graphical Optimal Transport (GGOT), for analyzing disease progressions. The proposed GGOT uses Gaussian graphical models, incorporating protein interaction networks, to characterize the data distributions at different disease stages. Then we use population-level optimal transport to calculate the Wasserstein distances and transport between stages, enabling us to detect critical transitions. By analyzing the per-molecule transport distance, we quantify the importance of each molecule and identify trigger molecules. Moreover, GGOT predicts the occurrence of critical transitions in unseen samples and visualizes the disease progression process. We apply GGOT to the simulation dataset and six disease datasets with varying disease progression rates to substantiate its effectiveness. Compared to existing methods, our proposed GGOT exhibits superior performance in detecting critical transitions.

Identifying critical States of complex diseases by local network Wasserstein distance

Complex diseases often undergo abrupt transitions from pre-disease to disease states, with the pre-disease state is typically unstable but potentially reversible through timely intervention. Detecting these critical transitions is crucial. We propose a model-free method, Local Network Wasserstein Distance (LNWD), for identifying critical transitions/pre-disease states in complex diseases using single sample analysis. LNWD measures statistical perturbations in normal samples caused by diseased samples using the Wasserstein distance, and identifies critical states by observing LNWD score changes. Applied to KIRP, KIRC, LUAD, ESCA (TCGA datasets) and GSE2565, GSE13268 (GEO datasets), the method successfully identified critical states in six disease datasets. This single-sample, local network-based approach provides early warning signals for medical diagnosis and holds great potential for personalized disease diagnosis.

Disease prediction by network information gain on a single sample basis

There are critical transition phenomena during the progression of many diseases. Such critical transitions are usually accompanied by catastrophic disease deterioration, and their prediction is of significant importance for disease prevention and treatment. However, predicting disease deterioration solely based on a single sample is a difficult problem. In this study, we presented the network information gain (NIG) method, for predicting the critical transitions or disease state based on network flow entropy from omics data of each individual. NIG can not only efficiently predict disease deteriorations but also detect their dynamic network biomarkers on an individual basis and further identify potential therapeutic targets. The numerical simulation demonstrates the effectiveness of NIG. Moreover, our method was validated by successfully predicting disease deteriorations and identifying their potential therapeutic targets from four real omics datasets, i.e., an influenza dataset and three cancer datasets.

Gene expression ranking change based single sample pre-disease state detection

Introduction

To prevent disease, it is of great importance to detect the critical point (pre-disease state) when the biological system abruptly transforms from normal to disease state. However, rapid and accurate pre-disease state detection is still a challenge when there is only a single sample available. The state transition of the biological system is driven by the variation in regulations between genes.

Methods

In this study, we propose a rapid single-sample pre-disease state-identifying method based on the change in gene expression ranking, which can reflect the coordinated shifts between genes, that is, S-PCR. The R codes of S-PCR can be accessed at https://github.com/ZhenshenBao/S-PCR.

Results

This model-free method is validated by the successful identification of pre-disease state for both simulated and five real datasets. The functional analyses of the pre-disease state-related genes identified by S-PCR also demonstrate the effectiveness of this computational approach. Furthermore, the time efficiency of S-PCR is much better than that of its peers.

Discussion

Hence, the proposed S-PCR approach holds immense potential for clinical applications in personalized disease diagnosis.

Multi-omics with dynamic network biomarker algorithm prefigures organ-specific metastasis of lung adenocarcinoma

Efficacious strategies for early detection of lung cancer metastasis are of significance for improving the survival of lung cancer patients. Here we show the marker genes and serum secretome foreshadowing the lung cancer site-specific metastasis through dynamic network biomarker (DNB) algorithm, utilizing two clinical cohorts of four major types of lung cancer distant metastases, with single-cell RNA sequencing (scRNA-seq) of primary lesions and liquid chromatography-mass spectrometry data of sera. Also, we locate the intermediate status of cancer cells, along with its gene signatures, in each metastatic state trajectory that cancer cells at this stage still have no specific organotropism. Furthermore, an integrated neural network model based on the filtered scRNA-seq data is successfully constructed and validated to predict the metastatic state trajectory of cancer cells. Overall, our study provides an insight to locate the pre-metastasis status of lung cancer and primarily examines its clinical application value, contributing to the early detection of lung cancer metastasis in a more feasible and efficacious way.

CPMI: comprehensive neighborhood-based perturbed mutual information for identifying critical states of complex biological processes

BACKGROUND

There exists a critical transition or tipping point during the complex biological process. Such critical transition is usually accompanied by the catastrophic consequences. Therefore, hunting for the tipping point or critical state is of significant importance to prevent or delay the occurrence of catastrophic consequences. However, predicting critical state based on the high-dimensional small sample data is a difficult problem, especially for single-cell expression data.

RESULTS

In this study, we propose the comprehensive neighbourhood-based perturbed mutual information (CPMI) method to detect the critical states of complex biological processes. The CPMI method takes into account the relationship between genes and neighbours, so as to reduce the noise and enhance the robustness. This method is applied to a simulated dataset and six real datasets, including an influenza dataset, two single-cell expression datasets and three bulk datasets. The method can not only successfully detect the tipping points, but also identify their dynamic network biomarkers (DNBs). In addition, the discovery of transcription factors (TFs) which can regulate DNB genes and nondifferential 'dark genes' validates the effectiveness of our method. The numerical simulation verifies that the CPMI method is robust under different noise strengths and is superior to the existing methods on identifying the critical states.

CONCLUSIONS

In conclusion, we propose a robust computational method, i.e., CPMI, which is applicable in both the bulk and single cell datasets. The CPMI method holds great potential in providing the early warning signals for complex biological processes and enabling early disease diagnosis.

Detecting tipping points of complex diseases by network information entropy

The progression of complex diseases often involves abrupt and non-linear changes characterized by sudden shifts that trigger critical transformations. Identifying these critical states or tipping points is crucial for understanding disease progression and developing effective interventions. To address this challenge, we have developed a model-free method named Network Information Entropy of Edges (NIEE). Leveraging dynamic network biomarkers, sample-specific networks, and information entropy theories, NIEE can detect critical states or tipping points in diverse data types, including bulk, single-sample expression data. By applying NIEE to real disease datasets, we successfully identified critical predisease stages and tipping points before disease onset. Our findings underscore NIEE's potential to enhance comprehension of complex disease development.

Edge-based relative entropy as a sensitive indicator of critical transitions in biological systems

BACKGROUND

Disease progression in biosystems is not always a steady process but is occasionally abrupt. It is important but challenging to signal critical transitions in complex biosystems.

METHODS

In this study, based on the theoretical framework of dynamic network biomarkers (DNBs), we propose a model-free method, edge-based relative entropy (ERE), to identify temporal key biomolecular associations/networks that may serve as DNBs and detect early-warning signals of the drastic state transition during disease progression in complex biological systems. Specifically, by combining gene‒gene interaction (edge) information with the relative entropy, the ERE method converts gene expression values into network entropy values, quantifying the dynamic change in a biomolecular network and indicating the qualitative shift in the system state.

RESULTS

The proposed method was validated using simulated data and real biological datasets of complex diseases. The applications show that for certain diseases, the ERE method helps to reveal so-called "dark genes" that are non-differentially expressed but with high ERE values and of essential importance in both gene regulation and prognosis.

CONCLUSIONS

The proposed method effectively identified the critical transition states of complex diseases at the network level. Our study not only identified the critical transition states of various cancers but also provided two types of new prognostic biomarkers, positive and negative edge biomarkers, for further practical application. The method in this study therefore has great potential in personalized disease diagnosis.

Data-driven energy landscape reveals critical genes in cancer progression

The evolution of cancer is a complex process characterized by stable states and transitions among them. Studying the dynamic evolution of cancer and revealing the mechanisms of cancer progression based on experimental data is an important topic. In this study, we aim to employ a data-driven energy landscape approach to analyze the dynamic evolution of cancer. We take Kidney renal clear cell carcinoma (KIRC) as an example. From the energy landscape, we introduce two quantitative indicators (transition probability and barrier height) to study critical shifts in KIRC cancer evolution, including cancer onset and progression, and identify critical genes involved in these transitions. Our results successfully identify crucial genes that either promote or inhibit these transition processes in KIRC. We also conduct a comprehensive biological function analysis on these genes, validating the accuracy and reliability of our predictions. This work has implications for discovering new biomarkers, drug targets, and cancer treatment strategies in KIRC.

MIWE: detecting the critical states of complex biological systems by the mutual information weighted entropy

Complex biological systems often undergo sudden qualitative changes during their dynamic evolution. These critical transitions are typically characterized by a catastrophic progression of the system. Identifying the critical point is critical to uncovering the underlying mechanisms of complex biological systems. However, the system may exhibit minimal changes in its state until the critical point is reached, and in the face of high throughput and strong noise data, traditional biomarkers may not be effective in distinguishing the critical state. In this study, we propose a novel approach, mutual information weighted entropy (MIWE), which uses mutual information between genes to build networks and identifies critical states by quantifying molecular dynamic differences at each stage through weighted differential entropy. The method is applied to one numerical simulation dataset and four real datasets, including bulk and single-cell expression datasets. The critical states of the system can be recognized and the robustness of MIWE method is verified by numerical simulation under the influence of different noises. Moreover, we identify two key transcription factors (TFs), CREB1 and CREB3, that regulate downstream signaling genes to coordinate cell fate commitment. The dark genes in the single-cell expression datasets are mined to reveal the potential pathway regulation mechanism.

mNFE: microbiome network flow entropy for detecting pre-disease states of type 1 diabetes

In the development of Type 1 diabetes (T1D), there are critical states just before drastic changes, and identifying these pre-disease states may predict T1D or provide crucial early-warning signals. Unlike gene expression data, gut microbiome data can be collected noninvasively from stool samples. Gut microbiome sequencing data contain different levels of phylogenetic information that can be utilized to detect the tipping point or critical state in a reliable manner, thereby providing accurate and effective early-warning signals. However, it is still difficult to detect the critical state of T1D based on gut microbiome data due to generally non-significant differences between healthy and critical states. To address this problem, we proposed a new method - microbiome network flow entropy (mNFE) based on a single sample from each individual - for detecting the critical state before seroconversion and abrupt transitions of T1D at various taxonomic levels. The numerical simulation validated the robustness of mNFE under different noise levels. Furthermore, based on real datasets, mNFE successfully identified the critical states and their dynamic network biomarkers (DNBs) at different taxonomic levels. In addition, we found some high-frequency species, which are closely related to the unique clinical characteristics of autoantibodies at the four levels, and identified some non-differential 'dark species' play important roles during the T1D progression. mNFE can robustly and effectively detect the pre-disease states at various taxonomic levels and identify the corresponding DNBs with only a single sample for each individual. Therefore, our mNFE method provides a new approach not only for T1D pre-disease diagnosis or preventative treatment but also for preventative medicine of other diseases by gut microbiome.

SGAE: single-cell gene association entropy for revealing critical states of cell transitions during embryonic development

The critical point or pivotal threshold of cell transition occurs in early embryonic development when cell differentiation culminates in its transition to specific cell fates, at which the cell population undergoes an abrupt and qualitative shift. Revealing such critical points of cell transitions can track cellular heterogeneity and shed light on the molecular mechanisms of cell differentiation. However, precise detection of critical state transitions proves challenging when relying on single-cell RNA sequencing data due to their inherent sparsity, noise, and heterogeneity. In this study, diverging from conventional methods like differential gene analysis or static techniques that emphasize classification of cell types, an innovative computational approach, single-cell gene association entropy (SGAE), is designed for the analysis of single-cell RNA-seq data and utilizes gene association information to reveal critical states of cell transitions. More specifically, through the translation of gene expression data into local SGAE scores, the proposed SGAE can serve as an index to quantitatively assess the resilience and critical properties of genetic regulatory networks, consequently detecting the signal of cell transitions. Analyses of five single-cell datasets for embryonic development demonstrate that the SGAE method achieves better performance in facilitating the characterization of a critical phase transition compared with other existing methods. Moreover, the SGAE value can effectively discriminate cellular heterogeneity over time and performs well in the temporal clustering of cells. Besides, biological functional analysis also indicates the effectiveness of the proposed approach.

Identifying the critical state of cancers by single-sample Markov flow entropy

Background

The progression of complex diseases sometimes undergoes a drastic critical transition, at which the biological system abruptly shifts from a relatively healthy state (before-transition stage) to a disease state (after-transition stage). Searching for such a critical transition or critical state is crucial to provide timely and effective scientific treatment to patients. However, in most conditions where only a small sample size of clinical data is available, resulting in failure when detecting the critical states of complex diseases, particularly only single-sample data.

Methods

In this study, different from traditional methods that require multiple samples at each time, a model-free computational method, single-sample Markov flow entropy (sMFE), provides a solution to the identification problem of critical states/pre-disease states of complex diseases, solely based on a single-sample. Our proposed method was employed to characterize the dynamic changes of complex diseases from the perspective of network entropy.

Results

The proposed approach was verified by unmistakably identifying the critical state just before the occurrence of disease deterioration for four tumor datasets from The Cancer Genome Atlas (TCGA) database. In addition, two new prognostic biomarkers, optimistic sMFE (O-sMFE) and pessimistic sMFE (P-sMFE) biomarkers, were identified by our method and enable the prognosis evaluation of tumors.

Conclusions

The proposed method has shown its capability to accurately detect pre-disease states of four cancers and provide two novel prognostic biomarkers, O-sMFE and P-sMFE biomarkers, to facilitate the personalized prognosis of patients. This is a remarkable achievement that could have a major impact on the diagnosis and treatment of complex diseases.

SPNE: sample-perturbed network entropy for revealing critical states of complex biological systems

Complex biological systems do not always develop smoothly but occasionally undergo a sharp transition; i.e. there exists a critical transition or tipping point at which a drastic qualitative shift occurs. Hunting for such a critical transition is important to prevent or delay the occurrence of catastrophic consequences, such as disease deterioration. However, the identification of the critical state for complex biological systems is still a challenging problem when using high-dimensional small sample data, especially where only a certain sample is available, which often leads to the failure of most traditional statistical approaches. In this study, a novel quantitative method, sample-perturbed network entropy (SPNE), is developed based on the sample-perturbed directed network to reveal the critical state of complex biological systems at the single-sample level. Specifically, the SPNE approach effectively quantifies the perturbation effect caused by a specific sample on the directed network in terms of network entropy and thus captures the criticality of biological systems. This model-free method was applied to both bulk and single-cell expression data. Our approach was validated by successfully detecting the early warning signals of the critical states for six real datasets, including four tumor datasets from The Cancer Genome Atlas (TCGA) and two single-cell datasets of cell differentiation. In addition, the functional analyses of signaling biomarkers demonstrated the effectiveness of the analytical and computational results.

MAPKAPK2, a potential dynamic network biomarker of α-synuclein prior to its aggregation in PD patients

One of the important pathological features of Parkinson's disease (PD) is the pathological aggregation of α-synuclein (α-Syn) in the substantia nigra. Preventing the aggregation of α-Syn has become a potential strategy for treating PD. However, the molecular mechanism of α-Syn aggregation is unclear. In this study, using the dynamic network biomarker (DNB) method, we first identified the critical time point when α-Syn undergoes pathological aggregation based on a SH-SY5Y cell model and found that DNB genes encode transcription factors that regulated target genes that were differentially expressed. Interestingly, we found that these DNB genes and their neighbouring genes were significantly enriched in the cellular senescence pathway and thus proposed that the DNB genes HSF1 and MAPKAPK2 regulate the expression of the neighbouring gene SERPINE1. Notably, in Gene Expression Omnibus (GEO) data obtained from substantia nigra, prefrontal cortex and peripheral blood samples, the expression level of MAPKAPK2 was significantly higher in PD patients than in healthy people, suggesting that MAPKAPK2 has potential as an early diagnostic biomarker of diseases related to pathological aggregation of α-Syn, such as PD. These findings provide new insights into the mechanisms underlying the pathological aggregation of α-Syn.

Dynamic Network Construction for Identifying Early Warning Signals Based On a Data-Driven Approach: Early Diagnosis Biomarker Discovery for Gastric Cancer

During the development of complex diseases, there is a critical transition from one status to another at a tipping point, which can be an early indicator of disease deterioration. To effectively enhance the performance of early risk identification, a novel dynamic network construction algorithm for identifying early warning signals based on a data-driven approach (EWS-DDA) was proposed. In EWS-DDA, the shrunken centroid was introduced to measure dynamic expression changes in assumed pathway reactions during the progression of complex disease for network construction and to define early warning signals by means of a data-driven approach. We applied EWS-DDA to perform a comprehensive analysis of gene expression profiles of gastric cancer (GC) from The Cancer Genome Atlas database and the Gene Expression Omnibus database. Six crucial genes were selected as potential biomarkers for the early diagnosis of GC. The experimental results of statistical analysis and biological analysis suggested that the six genes play important roles in GC occurrence and development. Then, EWS-DDA was compared with other state-of-the-art network methods to validate its performance. The theoretical analysis and comparison results suggested that EWS-DDA has great potential for a more complete presentation of disease deterioration and effective extraction of early warning information.

Detecting early-warning signals for influenza by dysregulated dynamic network biomarkers

As a dynamical system, complex disease always has a sudden state transition at the tipping point, which is the result of the long-term accumulation of abnormal regulations. This paper proposes a novel approach to detect the early-warning signals of influenza A (H3N2 and H1N1) outbreaks by dysregulated dynamic network biomarkers (dysregulated DNBs) for individuals. The results of cross-validation show that our approach can detect early-warning signals before the symptom appears successfully. Unlike the traditional DNBs, our dysregulated DNBs are anchored and very few, which is essential for disease early diagnosis in clinical practice. Moreover, the genes of dysregulated DNBs are significantly enriched in the influenza-related pathways. The source code of this paper can be freely downloaded from https://github.com/YanhaoHuo/dysregulated-DNBs.git.

Identifying the critical state of complex biological systems by the directed-network rank score method

MOTIVATION

Catastrophic transitions are ubiquitous in the dynamic progression of complex biological systems; that is, a critical transition at which complex systems suddenly shift from one stable state to another occurs. Identifying such a critical point or tipping point is essential for revealing the underlying mechanism of complex biological systems. However, it is difficult to identify the tipping point since few significant differences in the critical state are detected in terms of traditional static measurements.

RESULTS

In this study, by exploring the dynamic changes in gene cooperative effects between the before-transition and critical states, we presented a model-free approach, the directed-network rank score (DNRS), to detect the early-warning signal of critical transition in complex biological systems. The proposed method is applicable to both bulk and single-cell RNA-sequencing (scRNA-seq) data. This computational method was validated by the successful identification of the critical or pre-transition state for both simulated and six real datasets, including three scRNA-seq datasets of embryonic development and three tumor datasets. In addition, the functional and pathway enrichment analyses suggested that the corresponding DNRS signaling biomarkers were involved in key biological processes.

AVAILABILITY AND IMPLEMENTATION

The source code is freely available at https://github.com/zhongjiayuan/DNRS.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

TPD: a web tool for tipping-point detection based on dynamic network biomarker

Tipping points or critical transitions widely exist during the progression of many biological processes. It is of great importance to detect the tipping point with the measured omics data, which may be a key to achieving predictive or preventive medicine. We present the tipping point detector (TPD), a web tool for the detection of the tipping point during the dynamic process of biological systems, and further its leading molecules or network, based on the input high-dimensional time series or stage course data. With the solid theoretical background of dynamic network biomarker (DNB) and a series of computational methods for DNB detection, TPD detects the potential tipping point/critical state from the input omics data and outputs multifarious visualized results, including a suggested tipping point with a statistically significant P value, the identified key genes and their functional biological information, the dynamic change in the DNB/leading network that may drive the critical transition and the survival analysis based on DNB scores that may help to identify ‘dark’ genes (nondifferential in terms of expression but differential in terms of DNB scores). TPD fits all current browsers, such as Chrome, Firefox, Edge, Opera, Safari and Internet Explorer. TPD is freely accessible at http://www.rpcomputationalbiology.cn/TPD.

The single-sample network module biomarkers (sNMB) method reveals the pre-deterioration stage of disease progression

The progression of complex diseases generally involves a pre-deterioration stage that occurs during the transition from a healthy state to disease deterioration, at which a drastic and qualitative shift occurs. The development of an effective approach is urgently needed to identify such a pre-deterioration stage or critical state just before disease deterioration, which allows the timely implementation of appropriate measures to prevent a catastrophic transition. However, identifying the pre-deterioration stage is a challenging task in clinical medicine, especially when only a single sample is available for most patients, which is responsible for the failure of most statistical methods. In this study, a novel computational method, called single-sample network module biomarkers (sNMB), is presented to predict the pre-deterioration stage or critical point using only a single sample. Specifically, the proposed single-sample index effectively quantifies the disturbance caused by a single sample against a group of given reference samples. Our method successfully detected the early warning signal of the critical transitions when applied to both a numerical simulation and four real datasets, including acute lung injury, stomach adenocarcinoma, esophageal carcinoma, and rectum adenocarcinoma. In addition, it provides signaling biomarkers for further practical application, which helps to discover prognostic indicators and reveal the underlying molecular mechanisms of disease progression.

Detecting the Critical States of Type 2 Diabetes Mellitus Based on Degree Matrix Network Entropy by Cross-Tissue Analysis

Type 2 diabetes mellitus (T2DM) is a metabolic disease caused by multiple etiologies, the development of which can be divided into three states: normal state, critical state/pre-disease state, and disease state. To avoid irreversible development, it is important to detect the early warning signals before the onset of T2DM. However, detecting critical states of complex diseases based on high-throughput and strongly noisy data remains a challenging task. In this study, we developed a new method, i.e., degree matrix network entropy (DMNE), to detect the critical states of T2DM based on a sample-specific network (SSN). By applying the method to the datasets of three different tissues for experiments involving T2DM in rats, the critical states were detected, and the dynamic network biomarkers (DNBs) were successfully identified. Specifically, for liver and muscle, the critical transitions occur at 4 and 16 weeks. For adipose, the critical transition is at 8 weeks. In addition, we found some "dark genes" that did not exhibit differential expression but displayed sensitivity in terms of their DMNE score, which is closely related to the progression of T2DM. The information uncovered in our study not only provides further evidence regarding the molecular mechanisms of T2DM but may also assist in the development of strategies to prevent this disease.

A novel method to detect the early warning signal of COVID-19 transmission

BACKGROUND

Infectious illness outbreaks, particularly the corona-virus disease 2019 (COVID-19) pandemics in recent years, have wreaked havoc on human society, and the growing number of infected patients has put a strain on medical facilities. It's necessary to forecast early warning signals of potential outbreaks of COVID-19, which would facilitate the health ministry to take some suitable control measures timely to prevent or slow the spread of COVID-19. However, since the intricacy of COVID-19 transmission, which connects biological and social systems, it is a difficult task to predict outbreaks of COVID-19 epidemics timely.

RESULTS

In this work, we developed a new model-free approach, called, the landscape network entropy based on Auto-Reservoir Neural Network (ARNN-LNE), for quantitative analysis of COVID-19 propagation, by mining dynamic information from regional networks and short-term high-dimensional time-series data. Through this approach, we successfully identified the early warning signals in six nations or areas based on historical data of COVID-19 infections.

CONCLUSION

Based on the newly published data on new COVID-19 disease, the ARNN-LNE method can give early warning signals for the outbreak of COVID-19. It's worth noting that ARNN-LNE only relies on small samples data. Thus, it has great application potential for monitoring outbreaks of infectious diseases.

Research on Potential Network Markers and Signaling Pathways in Type 2 Diabetes Based on Conditional Cell-Specific Network

Traditional methods concerning type 2 diabetes (T2D) are limited to grouped cells instead of each single cell, and thus the heterogeneity of single cells is erased. Therefore, it is still challenging to study T2D based on a single-cell and network perspective. In this study, we construct a conditional cell-specific network (CCSN) for each single cell for the GSE86469 dataset which is a single-cell transcriptional set from nondiabetic (ND) and T2D human islet samples, and obtain a conditional network degree matrix (CNDM). Since beta cells are the key cells leading to T2D, we search for hub genes in CCSN of beta cells and find that ATP6AP2 is essential for regulation and storage of insulin, and the renin-angiotensin system involving ATP6AP2 is related to most pathological processes leading to diabetic nephropathy. The communication between beta cells and other endocrine cells is performed and three gene pairs with obvious interaction are found. In addition, different expression genes (DEGs) are found based on CNDM and the gene expression matrix (GEM), respectively. Finally, 'dark' genes are identified, and enrichment analysis shows that NFATC2 is involved in the VEGF signaling pathway and indirectly affects the production of Prostacyclin (PGI2), which may be a potential biomarker for diabetic nephropathy.

Identifying the critical states and dynamic network biomarkers of cancers based on network entropy

Background

There are sudden deterioration phenomena during the progression of many complex diseases, including most cancers; that is, the biological system may go through a critical transition from one stable state (the normal state) to another (the disease state). It is of great importance to predict this critical transition or the so-called pre-disease state so that patients can receive appropriate and timely medical care. In practice, however, this critical transition is usually difficult to identify due to the high nonlinearity and complexity of biological systems.

Methods

In this study, we employed a model-free computational method, local network entropy (LNE), to identify the critical transition/pre-disease states of complex diseases. From a network perspective, this method effectively explores the key associations among biomolecules and captures their dynamic abnormalities.

Results

Based on LNE, the pre-disease states of ten cancers were successfully detected. Two types of new prognostic biomarkers, optimistic LNE (O-LNE) and pessimistic LNE (P-LNE) biomarkers, were identified, enabling identification of the pre-disease state and evaluation of prognosis. In addition, LNE helps to find “dark genes” with nondifferential gene expression but differential LNE values.

Conclusions

The proposed method effectively identified the critical transition states of complex diseases at the single-sample level. Our study not only identified the critical transition states of ten cancers but also provides two types of new prognostic biomarkers, O-LNE and P-LNE biomarkers, for further practical application. The method in this study therefore has great potential in personalized disease diagnosis.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12967-022-03445-0.

Identifying the critical states of complex diseases by the dynamic change of multivariate distribution

The dynamics of complex diseases are not always smooth; they are occasionally abrupt, i.e. there is a critical state transition or tipping point at which the disease undergoes a sudden qualitative shift. There are generally a few significant differences in the critical state in terms of gene expressions or other static measurements, which may lead to the failure of traditional differential expression-based biomarkers to identify such a tipping point. In this study, we propose a computational method, the direct interaction network-based divergence, to detect the critical state of complex diseases by exploiting the dynamic changes in multivariable distributions inferred from observable samples and local biomolecular direct interaction networks. Such a method is model-free and applicable to both bulk and single-cell expression data. Our approach was validated by successfully identifying the tipping point just before the occurrence of a critical transition for both a simulated data set and seven real data sets, including those from The Cancer Genome Atlas and two single-cell RNA-sequencing data sets of cell differentiation. Functional and pathway enrichment analyses also validated the computational results from the perspectives of both molecules and networks.

Detecting the critical states during disease development based on temporal network flow entropy

Complex diseases progression can be generally divided into three states, which are normal state, predisease/critical state and disease state. The sudden deterioration of diseases can be viewed as a bifurcation or a critical transition. Therefore, hunting for the tipping point or critical state is of great importance to prevent the disease deterioration. However, it is still a challenging task to detect the critical states of complex diseases with high-dimensional data, especially based on an individual. In this study, we develop a new method based on network fluctuation of molecules, temporal network flow entropy (TNFE) or temporal differential network flow entropy, to detect the critical states of complex diseases on the basis of each individual. By applying this method to a simulated dataset and six real diseases, including respiratory viral infections and tumors with four time-course and two stage-course high-dimensional omics datasets, the critical states before deterioration were detected and their dynamic network biomarkers were identified successfully. The results on the simulated dataset indicate that the TNFE method is robust under different noise strengths, and is also superior to the existing methods on detecting the critical states. Moreover, the analysis on the real datasets demonstrated the effectiveness of TNFE for providing early-warning signals on various diseases. In addition, we also predicted disease deterioration risk and identified drug targets for cancers based on stage-wise data.

Development of a dynamic network biomarkers method and its application for detecting the tipping point of prior disease development

The dynamic network biomarker (DNB) method has advanced since it was first proposed. This review discusses advances in the DNB method that can identify the dynamic change in the expression signature related to the critical time point of disease progression by utilizing different kinds of transcriptome data. The DNB method is good at identifying potential biomarkers for cancer and other disease development processes that are represented by a limited molecular profile change between the normal and critical stages. We highlight that the cancer tipping point or premalignant state has been widely discovered for different types of cancer by using the DNB method that utilizes bulk or single-cell RNA sequencing data. This method could also be applied to other dynamic research studies and help identify early warning signals, such as the prediction of a pre-outbreak of COVID-19. We also discuss how the identification of reliable biomarkers of cancer and the development of new methods can be utilized for early detection and intervention and provide insights into emerging paths of the widespread biomarker candidate pool for further validation and disease/health management.

Transcriptome analysis method based on differential distribution evaluation

Identifying differential genes over conditions provides insights into the mechanisms of biological processes and disease progression. Here we present an approach, the Kullback-Leibler divergence-based differential distribution (klDD), which provides a flexible framework for quantifying changes in higher-order statistical information of genes including mean and variance/covariation. The method can well detect subtle differences in gene expression distributions in contrast to mean or variance shifts of the existing methods. In addition to effectively identifying informational genes in terms of differential distribution, klDD can be directly applied to cancer subtyping, single-cell clustering and disease early-warning detection, which were all validated by various benchmark datasets.

Detect the early-warning signals of diseases based on signaling pathway perturbations on a single sample

BACKGROUND

During the pathogenesisof complex diseases, a sudden health deterioration will occur as results of the cumulative effect of various internal or external factors. The prediction of an early warning signal (pre-disease state) before such deterioration is very important in clinical practice, especially for a single sample. The single-sample landscape entropy (SLE) was proposed to tackle this issue. However, the PPI used in SLE was lack of definite biological meanings. Besides, the calculation of multiple correlations based on limited reference samples in SLE is time-consuming and suspect.

RESULTS

Abnormal signals generally exert their effect through the static definite biological functions in signaling pathways across the development of diseases. Thus, it is a natural way to study the propagation of the early-warning signals based on the signaling pathways in the KEGG database. In this paper, we propose a signaling perturbation method named SSP, to study the early-warning signal in signaling pathways for single dynamic time-series data. Results in three real datasets including the influenza virus infection, lung adenocarcinoma, and acute lung injury show that the proposed SSP outperformed the SLE. Moreover, the early-warning signal can be detected by one important signaling pathway PI3K-Akt.

CONCLUSIONS

These results all indicate that the static model in pathways could simplify the detection of the early-warning signals.

Dynamical network biomarkers: Theory and applications

This paper reviews theory of DNB (Dynamical Network Biomarkers) and its applications including both modern medicine and traditional medicine. We show that omics data such as gene/protein expression profiles can be effectively used to detect pre-disease states before critical transitions from healthy states to disease states by using the DNB theory. The DNB theory with big biological data is expected to lead to ultra-early precision and preventive medicine.

Single-cell transcriptomics reveal DHX9 in mature B cell as a dynamic network biomarker before lymph node metastasis in CRC

Increasing evidence indicates that mature B cells in the adjacent tumor tissue, both as an intermediate state, are vital in advanced colorectal cancer (CRC), which is associated with a low survival rate. Developing predictive biomarkers that detect the tipping point of mature B cells before lymph node metastasis in CRC is critical to prevent irreversible deterioration. We analyzed B cells in the adjacent tissues of CRC samples from different stages using the dynamic network biomarker (DNB) method. Single-cell profiling of 725 CRC-derived B cells revealed the emergence of a mature B cell subtype. Using the DNB method, we identified stage II as a critical period before lymph node metastasis and that reversed difference genes triggered by DNBs were enriched in the Janus kinase (JAK)-signal transducer and activator of transcription (STAT) pathway involving B cell immune capability. DHX9 (DEAH-box helicase 9) was a specific para-cancerous tissue DNB key gene. The dynamic expression levels of DHX9 and its proximate network genes involved in B cell-related pathways were reversed at the network level from stage I to III. In summary, DHX9 in mature B cells of CRC-adjacent tissues may serve as a predictable biomarker and a potential immune target in CRC progression.

Facial Skin Microbiota-Mediated Host Response to Pollution Stress Revealed by Microbiome Networks of Individual

Urban living has been reported to cause various skin disorders. As an integral part of the skin barrier, the skin microbiome is among the key factors associated with urbanization-related skin alterations. The role of skin microbiome in mediating the effect of urban stressors (e.g., air pollutants) on skin physiology is not well understood. We generated 16S sequencing data and constructed a microbiome network of individual (MNI) to analyze the effect of pollution stressors on the microbiome network and its downstream mediation effect on skin physiology in a personalized manner. In particular, we found that the connectivity and fragility of MNIs significantly mediated the adverse effects of air pollution on skin health, and a smoking lifestyle deepened the negative effects of pollution stress on facial skin microbiota. This is the first study that describes the mediation effect of the microbiome network on the skin’s physiological response toward environmental factors as revealed by our newly developed MNI approach and conditional process analysis.

IMPORTANCE The association between the skin microbiome and skin health has been widely reported. However, the role of the skin microbiome in mediating skin physiology remains a challenging and yet priority subject in the field. Through developing a novel MNI method followed by mediation analysis, we characterized the network signature of the skin microbiome at an individual level and revealed the role of the skin microbiome in mediating the skin’s responses toward environmental stressors. Our findings may shed new light on microbiome functions in skin health and lay the foundation for the design of a microbiome-based intervention strategy in the future.

scGET: Predicting Cell Fate Transition During Early Embryonic Development by Single-cell Graph Entropy

During early embryonic development, cell fate commitment represents a critical transition or "tipping point" of embryonic differentiation, at which there is a drastic and qualitative shift of the cell populations. In this study, we presented a computational approach, scGET, to explore the gene-gene associations based on single-cell RNA sequencing (scRNA-seq) data for critical transition prediction. Specifically, by transforming the gene expression data to the local network entropy, the single-cell graph entropy (SGE) value quantitatively characterizes the stability and criticality of gene regulatory networks among cell populations and thus can be employed to detect the critical signal of cell fate or lineage commitment at the single-cell level. Being applied to five scRNA-seq datasets of embryonic differentiation, scGET accurately predicts all the impending cell fate transitions. After identifying the "dark genes" that are non-differentially expressed genes but sensitive to the SGE value, the underlying signaling mechanisms were revealed, suggesting that the synergy of dark genes and their downstream targets may play a key role in various cell development processes.The application in all five datasets demonstrates the effectiveness of scGET in analyzing scRNA-seq data from a network perspective and its potential to track the dynamics of cell differentiation. The source code of scGET is accessible at https://github.com/zhongjiayuna/scGET_Project.

Resilience function uncovers the critical transitions in cancer initiation

Considerable evidence suggests that during the progression of cancer initiation, the state transition from wellness to disease is not necessarily smooth but manifests switch-like nonlinear behaviors, preventing the cancer prediction and early interventional therapy for patients. Understanding the mechanism of such wellness-to-disease transitions is a fundamental and challenging task. Despite the advances in flux theory of nonequilibrium dynamics and 'critical slowing down'-based system resilience theory, a system-level approach still lacks to fully describe this state transition. Here, we present a novel framework (called bioRFR) to quantify such wellness-to-disease transition during cancer initiation through uncovering the biological system's resilience function from gene expression data. We used bioRFR to reconstruct the biologically and dynamically significant resilience functions for cancer initiation processes (e.g. BRCA, LUSC and LUAD). The resilience functions display the similar resilience pattern with hysteresis feature but different numbers of tipping points, which implies that once the cell become cancerous, it is very difficult or even impossible to reverse to the normal state. More importantly, bioRFR can measure the severe degree of cancer patients and identify the personalized key genes that are associated with the individual system's state transition from normal to tumor in resilience perspective, indicating that bioRFR can contribute to personalized medicine and targeted cancer therapy.

c-CSN: Single-cell RNA Sequencing Data Analysis by Conditional Cell-specific Network

The rapid advancement of single-cell technologies has shed new light on the complex mechanisms of cellular heterogeneity. However, compared to bulk RNA sequencing (RNA-seq), single-cell RNA-seq (scRNA-seq) suffers from higher noise and lower coverage, which brings new computational difficulties. Based on statistical independence, cell-specific network (CSN) is able to quantify the overall associations between genes for each cell, yet suffering from a problem of overestimation related to indirect effects. To overcome this problem, we propose the c-CSN method, which can construct the conditional cell-specific network (CCSN) for each cell. c-CSN method can measure the direct associations between genes by eliminating the indirect associations. c-CSN can be used for cell clustering and dimension reduction on a network basis of single cells. Intuitively, each CCSN can be viewed as the transformation from less “reliable” gene expression to more “reliable” gene–gene associations in a cell. Based on CCSN, we further design network flow entropy (NFE) to estimate the differentiation potency of a single cell. A number of scRNA-seq datasets were used to demonstrate the advantages of our approach. 1) One direct association network is generated for one cell. 2) Most existing scRNA-seq methods designed for gene expression matrices are also applicable to c-CSN-transformed degree matrices. 3) CCSN-based NFE helps resolving the direction of differentiation trajectories by quantifying the potency of each cell. c-CSN is publicly available at https://github.com/LinLi-0909/c-CSN.

Real-Time Forecast of Influenza Outbreak Using Dynamic Network Marker Based on Minimum Spanning Tree

The influenza pandemic is a wide-ranging threat to people's health and property all over the world. Developing effective strategies for predicting the influenza outbreak which may prevent or at least get ready for a new influenza pandemic is now a top global public health priority. Owing to the complexity of influenza outbreaks that are usually involved with spatial and temporal characteristics of both biological and social systems, however, it is a challenging task to achieve the real-time monitoring of influenza outbreaks. In this study, by exploring the rich dynamical information of the city network during influenza outbreaks, we developed a computational method, the minimum-spanning-tree-based dynamical network marker (MST-DNM), to identify the tipping point or critical stage prior to the influenza outbreak. With historical records of influenza outpatients between 2009 and 2018, the MST-DNM strategy has been validated by accurate predictions of the influenza outbreaks in three Japanese cities/regions, respectively, i.e., Tokyo, Osaka, and Hokkaido. These successful applications show that the early-warning signal was detected 4 weeks on average ahead of each influenza outbreak. The results show that our method is of considerable potential in the practice of public health surveillance.

Autoreservoir computing for multistep ahead prediction based on the spatiotemporal information transformation

We develop an auto-reservoir computing framework, Auto-Reservoir Neural Network (ARNN), to efficiently and accurately make multi-step-ahead predictions based on a short-term high-dimensional time series. Different from traditional reservoir computing whose reservoir is an external dynamical system irrelevant to the target system, ARNN directly transforms the observed high-dimensional dynamics as its reservoir, which maps the high-dimensional/spatial data to the future temporal values of a target variable based on our spatiotemporal information (STI) transformation. Thus, the multi-step prediction of the target variable is achieved in an accurate and computationally efficient manner. ARNN is successfully applied to both representative models and real-world datasets, all of which show satisfactory performance in the multi-step-ahead prediction, even when the data are perturbed by noise and when the system is time-varying. Actually, such ARNN transformation equivalently expands the sample size and thus has great potential in practical applications in artificial intelligence and machine learning.

Detecting the outbreak of influenza based on the shortest path of dynamic city network

The influenza pandemic causes a large number of hospitalizations and even deaths. There is an urgent need for an efficient and effective method for detecting the outbreak of influenza so that timely, appropriate interventions can be made to prevent or at least prepare for catastrophic epidemics. In this study, we proposed a computational method, the shortest-path-based dynamical network marker (SP-DNM), to detect the pre-outbreak state of influenza epidemics by monitoring the dynamical change of the shortest path in a city network. Specifically, by mapping the real-time information to a properly constructed city network, our method detects the early-warning signal prior to the influenza outbreak in both Tokyo and Hokkaido for consecutive 9 years, which demonstrate the effectiveness and robustness of the proposed method.

Single-Sample Node Entropy for Molecular Transition in Pre-deterioration Stage of Cancer

A complex disease, especially cancer, always has pre-deterioration stage during its progression, which is difficult to identify but crucial to drug research and clinical intervention. However, using a few samples to find mechanisms that propel cancer crossing the pre-deterioration stage is still a complex problem. In this study, we successfully developed a novel single-sample model based on node entropy with a priori established protein interaction network. Using this model, critical stages were successfully detected in simulation data and four TCGA datasets, indicating its sensitivity and robustness. Besides, compared with the results of the differential analysis, our results showed that most of dynamic network biomarkers identified by node entropy, such as NKD2 or DAAM1, located in upstream in many important cancer-related signaling pathways regulated intergenic signaling within pathways. We also identified some novel prognostic biomarkers such as PER2, TNFSF4, MMP13 and ENO4 using node entropy rather than expression level. More importantly, we found the switch of non-specific pathways related to DNA damage repairing was the main driven force for cancer progression. In conclusion, we have successfully developed a dynamic node entropy model based on single case data to find out tipping point and possible mechanism for cancer progression. These findings may provide new target genes in therapeutic intervention tactics.

Disease characterization using a partial correlation-based sample-specific network

A single-sample network (SSN) is a biological molecular network constructed from single-sample data given a reference dataset and can provide insights into the mechanisms of individual diseases and aid in the development of personalized medicine. In this study, we proposed a computational method, a partial correlation-based single-sample network (P-SSN), which not only infers a network from each single-sample data given a reference dataset but also retains the direct interactions by excluding indirect interactions (https://github.com/hyhRise/P-SSN). By applying P-SSN to analyze tumor data from the Cancer Genome Atlas and single cell data, we validated the effectiveness of P-SSN in predicting driver mutation genes (DMGs), producing network distance, identifying subtypes and further classifying single cells. In particular, P-SSN is highly effective in predicting DMGs based on single-sample data. P-SSN is also efficient for subtyping complex diseases and for clustering single cells by introducing network distance between any two samples.

Identifying critical state of complex diseases by single-sample Kullback-Leibler divergence

Background

Developing effective strategies for signaling the pre-disease state of complex diseases, a state with high susceptibility before the disease onset or deterioration, is urgently needed because such state usually followed by a catastrophic transition into a worse stage of disease. However, it is a challenging task to identify such pre-disease state or tipping point in clinics, where only one single sample is available and thus results in the failure of most statistic approaches.

Methods

In this study, we presented a single-sample-based computational method to detect the early-warning signal of critical transition during the progression of complex diseases. Specifically, given a set of reference samples which were regarded as background, a novel index called single-sample Kullback–Leibler divergence (sKLD), was proposed to explore and quantify the disturbance on the background caused by a case sample. The pre-disease state is then signaled by the significant change of sKLD.

Results

The novel algorithm was developed and applied to both numerical simulation and real datasets, including lung squamous cell carcinoma, lung adenocarcinoma, stomach adenocarcinoma, thyroid carcinoma, colon adenocarcinoma, and acute lung injury. The successful identification of pre-disease states and the corresponding dynamical network biomarkers for all six datasets validated the effectiveness and accuracy of our method.

Conclusions

The proposed method effectively explores and quantifies the disturbance on the background caused by a case sample, and thus characterizes the criticality of a biological system. Our method not only identifies the critical state or tipping point at a single sample level, but also provides the sKLD-signaling markers for further practical application. It is therefore of great potential in personalized pre-disease diagnosis.