Early diagnosis of complex diseases by molecular biomarkers, network biomarkers, and dynamical network biomarkers

Sensors and Devices Guided by Artificial Intelligence for Personalized Pain Medicine

Personalized pain medicine aims to tailor pain treatment strategies for the specific needs and characteristics of an individual patient, holding the potential for improving treatment outcomes, reducing side effects, and enhancing patient satisfaction. Despite existing pain markers and treatments, challenges remain in understanding, detecting, and treating complex pain conditions. Here, we review recent engineering efforts in developing various sensors and devices for addressing challenges in the personalized treatment of pain. We summarize the basics of pain pathology and introduce various sensors and devices for pain monitoring, assessment, and relief. We also discuss advancements taking advantage of rapidly developing medical artificial intelligence (AI), such as AI-based analgesia devices, wearable sensors, and healthcare systems. We believe that these innovative technologies may lead to more precise and responsive personalized medicine, greatly improved patient quality of life, increased efficiency of medical systems, and reducing the incidence of addiction and substance use disorders.

Inference of single-cell network using mutual information for scRNA-seq data analysis

BACKGROUND

With the advance in single-cell RNA sequencing (scRNA-seq) technology, deriving inherent biological system information from expression profiles at a single-cell resolution has become possible. It has been known that network modeling by estimating the associations between genes could better reveal dynamic changes in biological systems. However, accurately constructing a single-cell network (SCN) to capture the network architecture of each cell and further explore cell-to-cell heterogeneity remains challenging.

RESULTS

We introduce SINUM, a method for constructing the SIngle-cell Network Using Mutual information, which estimates mutual information between any two genes from scRNA-seq data to determine whether they are dependent or independent in a specific cell. Experiments on various scRNA-seq datasets with different cell numbers based on eight performance indexes (e.g., adjusted rand index and F-measure index) validated the accuracy and robustness of SINUM in cell type identification, superior to the state-of-the-art SCN inference method. Additionally, the SINUM SCNs exhibit high overlap with the human interactome and possess the scale-free property.

CONCLUSIONS

SINUM presents a view of biological systems at the network level to detect cell-type marker genes/gene pairs and investigate time-dependent changes in gene associations during embryo development. Codes for SINUM are freely available at https://github.com/SysMednet/SINUM .

Evaluation of serum mid-infrared spectroscopy as new prognostic marker for first-line bevacizumab-based chemotherapy in metastatic colorectal cancer

BACKGROUND AND AIMS

Bevacizumab-based chemotherapy is a recommended first-line treatment for metastatic colorectal cancer (mCRC). Robust biomarkers with clinical practice applicability have not been identified for patients with this treatment. We aimed to evaluate the prognostic yield of serum mid-infrared spectroscopy (MIRS) on patients receiving first-line bevacizumab-based chemotherapy for mCRC.

METHODS

We conducted an ancillary analysis from a multicentre prospective study (NCT00489697). All baseline serums were screened by attenuated total reflection method. Principal component analysis and unsupervised k-mean partitioning methods were performed blinded to all patients' data. Endpoints were progression-free survival (PFS) and overall survival (OS).

RESULTS

From the 108 included patients, MIRS discriminated two prognostic groups. First group patients had significantly lower body mass index (p = 0.026) and albumin levels (p < 0.001), and higher levels of angiogenic markers, lactate dehydrogenase and carcinoembryonic antigen (p < 0.001). In univariate analysis, their OS and PFS were shorter with respective medians: 17.6 vs 27.9 months (p = 0.02) and 8.7 vs 11.3 months (p = 0.03). In multivariate analysis, PFS was significantly shorter (HR = 1.74, p = 0.025) with a similar trend for OS (HR = 1.69, p = 0.061).

CONCLUSION

By metabolomic fingerprinting, MIRS proves to be a promising prognostic tool for patients receiving first-line bevacizumab-based chemotherapy for mCRC.

Characterization of Inflammation/Immune-, Acute Phase-, Extracellular Matrix-, Adhesion-, and Serine Protease-Related Proteins in the Amniotic Fluid of Women With Early Preterm Prelabor Rupture of Membranes

PROBLEM

To determine whether altered concentrations of various inflammation/immune-, acute phase-, extracellular matrix-, adhesion-, and serine protease-related proteins in the amniotic fluid (AF) are independently associated with microbial invasion of the amniotic cavity and/or intra-amniotic inflammation (MIAC/IAI), imminent spontaneous preterm delivery (SPTD; ≤7 days), and major neonatal morbidity/mortality (NMM) in women with early preterm prelabor rupture of membranes (PPROM).

METHOD OF STUDY

This was a retrospective cohort study involving 111 singleton pregnant women with PPROM (24-31 weeks) undergoing amniocentesis to diagnose MIAC/IAI. The following proteins were measured in stored AF samples by enzyme-linked immunosorbent assay (ELISA): APRIL, DKK-3, Gal-3BP, IGFBP-2, IL-8, VDBP, lumican, MMP-2, MMP-8, SPARC, TGFBI, TGF-β1, E-selectin, ICAM-5, P-selectin, haptoglobin, hepcidin, SAA1, kallistatin, and uPA.

RESULTS

Multivariate logistic regression analyses revealed that (i) elevated APRIL, IL-8, MMP-8, and TGFBI levels in the AF, reduced lumican and SPARC levels in the AF, and high percentages of samples above the lower limit of quantification for AF TGF-β1 and uPA were significantly associated with MIAC/IAI; (ii) elevated AF levels of IL-8 and MMP-8 were significantly associated with SPTD within 7 days; and (iii) elevated AF IL-6 levels were significantly associated with increased risk for major NMM, when adjusted for baseline covariates.

CONCLUSION

ECM (lumican, SPRAC, TGFBI, and TGF-β1)- and serine protease (uPA)-associated proteins in the AF are involved in the regulation of the host response to infection/inflammation in the amniotic cavity, whereas AF inflammation (IL-8, MMP-8, and IL-6)-associated mediators are implicated in the development of preterm parturition and major NMM in early PPROM.

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.

Personalized decision support system for tailoring IgA nephropathy treatment strategies

BACKGROUND

The ongoing debate surrounding the use of immunosuppressive treatments for IgA nephropathy (IgAN) underscores the demand for personalized and effective strategies.

METHODS

Analyzed data from 807 IgAN patients over 5+ years using three methods: Random Forest with molecular biomarkers, network biomarkers with graph engineering, and an auto-encoder model. All models were trained using identical demographic, clinical, and pathological data, employing an 80-20 split for training and testing purposes.

RESULTS

In the comprehensive assessment of IgAN prognosis, the Random Forest model, employing molecular biomarkers, demonstrated strong performance metrics (AUC = 0.83, sensitivity = 0.51, specificity = 0.96). However, traditional graph feature engineering on patient-specific networks outperformed these results with an AUC of 0.90, sensitivity of 0.64, and specificity of 0.94. The Auto-encoder model showed the best accuracy (AUC = 0.91, sensitivity = 0.46, specificity = 0.96). The findings highlighted the superior predictive capabilities of network biomarkers over molecular biomarkers for adverse renal outcome prediction in IgAN. Consequently, we integrated Auto-encoder-derived Network Biomarkers with Random Forest Models to enhance prognostic precision in diverse IgAN treatment scenarios. The prediction for the prognosis of patients receiving supportive care, glucocorticoid therapy, and immunosuppressant treatment yielded AUC values of 0.95, 0.96, and 1, respectively, indicating high specificity. Drawing from these insights, we pioneered the development of an innovative decision support model for IgAN treatment. This model demonstrated the ability to make medical decisions comparable to those by experienced nephrologists, enabling the customization of personalized disease management strategies.

CONCLUSION

Our system accurately predicted IgAN prognosis and evaluated various treatment efficacies, aiding physicians in devising optimal therapeutic strategies for patients.

Single-cell biclustering for cell-specific transcriptomic perturbation detection in AD progression

The pathogenesis of Alzheimer disease (AD) involves complex gene regulatory changes across different cell types. To help decipher this complexity, we introduce single-cell Bayesian biclustering (scBC), a framework for identifying cell-specific gene network biomarkers in scRNA and snRNA-seq data. Through biclustering, scBC enables the analysis of perturbations in functional gene modules at the single-cell level. Applying the scBC framework to AD snRNA-seq data reveals the perturbations within gene modules across distinct cell groups and sheds light on gene-cell correlations during AD progression. Notably, our method helps to overcome common challenges in single-cell data analysis, including batch effects and dropout events. Incorporating prior knowledge further enables the framework to yield more biologically interpretable results. Comparative analyses on simulated and real-world datasets demonstrate the precision and robustness of our approach compared to other state-of-the-art biclustering methods. scBC holds potential for unraveling the mechanisms underlying polygenic diseases characterized by intricate gene coexpression patterns.

Machine-learning-based plasma metabolomic profiles for predicting long-term complications of cirrhosis

BACKGROUND AND AIMS

The complications of liver cirrhosis occur after long asymptomatic stages of progressive fibrosis and are generally diagnosed late. We aimed to develop a plasma metabolomic-based score tool to predict these events.

APPROACH AND RESULTS

We enrolled 64,005 UK biobank participants with metabolomic profiles. Participants were randomly divided into the training (n=43,734) and validation cohorts (n=20,271). Liver cirrhosis complications were defined as hospitalization for liver cirrhosis or presentation with HCC. An interpretable machine-learning framework was applied to learn the metabolomic states extracted from 168 circulating metabolites in the training cohort. An integrated nomogram was developed and compared to conventional and genetic risk scores. We created 3 groups: low-risk, middle-risk, and high-risk through selected cutoffs of the nomogram. The predictive performance was validated through the area under a time-dependent receiver operating characteristic curve (time-dependent AUC), calibration curves, and decision curve analysis. The metabolomic state model could accurately predict the 10-year risk of liver cirrhosis complications in the training cohort (time-dependent AUC: 0.84 [95% CI: 0.82-0.86]), and outperform the fibrosis-4 index (time-dependent AUC difference: 0.06 [0.03-0.10]) and polygenic risk score (0.25 [0.21-0.29]). The nomogram, integrating metabolomic state, aspartate aminotransferase, platelet count, waist/hip ratio, and smoking status showed a time-dependent AUC of 0.930 at 3 years, 0.889 at 5 years, and 0.861 at 10 years in the validation cohort, respectively. The HR in the high-risk group was 43.58 (95% CI: 27.08-70.12) compared with the low-risk group.

CONCLUSIONS

We developed a metabolomic state-integrated nomogram, which enables risk stratification and personalized administration of liver-related events.

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.

An assessment system for clinical and biological interpretability in ulcerative colitis

Ulcerative colitis (UC) is a serious inflammatory bowel disease (IBD) with high morbidity and mortality worldwide. As the traditional diagnostic techniques have various limitations in the practice and diagnosis of early ulcerative colitis, it is necessary to develop new diagnostic models from molecular biology to supplement the existing methods. In this study, we developed a machine learning-based synthesis to construct an artificial intelligence diagnostic model for ulcerative colitis, and the correctness of the model is verified using an external independent dataset. According to the significantly expressed genes related to the occurrence of UC in the model, an unsupervised quantitative ulcerative colitis related score (UCRScore) based on principal coordinate analysis was established. The UCRScore is not only highly generalizable across UC bulk cohorts at different stages, but also highly generalizable across single-cell datasets, with the same effect in terms of cell numbers, activation pathways and mechanisms. As an important role of screening genes in disease occurrence, based on connectivity map analysis, 5 potential targeting molecular compounds were identified, which can be used as an additional supplement to the therapeutic of UC. Overall, this study provides a potential tool for differential diagnosis and assessment of bio-pathological changes in UC at the macroscopic level, providing an opportunity to optimize the diagnosis and treatment of UC.

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.

Identification and Characterization of CD8+ CD27+ CXCR3- T Cell Dysregulation and Progression-Associated Biomarkers in Systemic Lupus Erythematosus

Systemic Lupus Erythematosus (SLE) etiopathogenesis highlights the contributions of overproduction of CD4+ T cells and loss of immune tolerance. However, the involvement of CD8+ T cells in SLE pathology and disease progression remains unclear. Here, the comprehensive immune cell dysregulation in total 263 clinical peripheral blood samples composed of active SLE (aSLE), remission SLE (rSLE) and healthy controls (HCs) is investigated via mass cytometry, flow cytometry and single-cell RNA sequencing. This is observed that CD8+ CD27+ CXCR3- T cells are increased in rSLE compare to aSLE. Meanwhile, the effector function of CD8+ CD27+ CXCR3- T cells are overactive in aSLE compare to HCs and rSLE, and are positively associated with clinical SLE activity. In addition, the response of peripheral blood mononuclear cells (PBMCs) is monitored to interleukin-2 stimulation in aSLE and rSLE to construct dynamic network biomarker (DNB) model. It is demonstrated that DNB score-related parameters can faithfully predict the remission of aSLE and the flares of rSLE. The abundance and functional dysregulation of CD8+ CD27+ CXCR3- T cells can be potential biomarkers for SLE prognosis and concomitant diagnosis. The DNB score with accurate prediction to SLE disease progression can provide clinical treatment suggestions especially for drug dosage determination.

Specific network information gain for detecting the critical state of colorectal cancer based on gut microbiome

There generally exists a critical state or tipping point from a stable state to another in the development of colorectal cancer (CRC) beyond which a significant qualitative transition occurs. Gut microbiome sequencing data can be collected non-invasively from fecal samples, making it more convenient to obtain. Furthermore, intestinal microbiome sequencing data contain phylogenetic information at various levels, which can be used to reliably identify critical states, thereby providing early warning signals more accurately and effectively. Yet, pinpointing the critical states using gut microbiome data presents a formidable challenge due to the high dimension and strong noise of gut microbiome data. To address this challenge, we introduce a novel approach termed the specific network information gain (SNIG) method to detect CRC's critical states at various taxonomic levels via gut microbiome data. The numerical simulation indicates that the SNIG method is robust under different noise levels and that it is also superior to the existing methods on detecting the critical states. Moreover, utilizing SNIG on two real CRC datasets enabled us to discern the critical states preceding deterioration and to successfully identify their associated dynamic network biomarkers at different taxonomic levels. Notably, we discovered certain 'dark species' and pathways intimately linked to CRC progression. In addition, we accurately detected the tipping points on an individual dataset of type I diabetes.

Use of 3D Printing Techniques to Fabricate Implantable Microelectrodes for Electrochemical Detection of Biomarkers in the Early Diagnosis of Cardiovascular and Neurodegenerative Diseases

This Review provides a comprehensive overview of 3D printing techniques to fabricate implantable microelectrodes for the electrochemical detection of biomarkers in the early diagnosis of cardiovascular and neurodegenerative diseases. Early diagnosis of these diseases is crucial to improving patient outcomes and reducing healthcare systems' burden. Biomarkers serve as measurable indicators of these diseases, and implantable microelectrodes offer a promising tool for their electrochemical detection. Here, we discuss various 3D printing techniques, including stereolithography (SLA), digital light processing (DLP), fused deposition modeling (FDM), selective laser sintering (SLS), and two-photon polymerization (2PP), highlighting their advantages and limitations in microelectrode fabrication. We also explore the materials used in constructing implantable microelectrodes, emphasizing their biocompatibility and biodegradation properties. The principles of electrochemical detection and the types of sensors utilized are examined, with a focus on their applications in detecting biomarkers for cardiovascular and neurodegenerative diseases. Finally, we address the current challenges and future perspectives in the field of 3D-printed implantable microelectrodes, emphasizing their potential for improving early diagnosis and personalized treatment strategies.

Analyzing omics data by feature combinations based on kernel functions

Defining meaningful feature (molecule) combinations can enhance the study of disease diagnosis and prognosis. However, feature combinations are complex and various in biosystems, and the existing methods examine the feature cooperation in a single, fixed pattern for all feature pairs, such as linear combination. To identify the appropriate combination between two features and evaluate feature combination more comprehensively, this paper adopts kernel functions to study feature relationships and proposes a new omics data analysis method KF-[Formula: see text]-TSP. Besides linear combination, KF-[Formula: see text]-TSP also explores the nonlinear combination of features, and allows hybridizing multiple kernel functions to evaluate feature interaction from multiple views. KF-[Formula: see text]-TSP selects [Formula: see text] > 0 top-scoring pairs to build an ensemble classifier. Experimental results show that KF-[Formula: see text]-TSP with multiple kernel functions which evaluates feature combinations from multiple views is better than that with only one kernel function. Meanwhile, KF-[Formula: see text]-TSP performs better than TSP family algorithms and the previous methods based on conversion strategy in most cases. It performs similarly to the popular machine learning methods in omics data analysis, but involves fewer feature pairs. In the procedure of physiological and pathological changes, molecular interactions can be both linear and nonlinear. Hence, KF-[Formula: see text]-TSP, which can measure molecular combination from multiple perspectives, can help to mine information closely related to physiological and pathological changes and study disease mechanism.

A Single-Molecule Bioelectronic Portable Array for Early Diagnosis of Pancreatic Cancer Precursors

A cohort of 47 patients is screened for pancreatic cancer precursors with a portable 96-well bioelectronic sensing-array for single-molecule assay in cysts fluid and blood plasma, deployable at point-of-care (POC). Pancreatic cancer precursors are mucinous cysts diagnosed with a sensitivity of at most 80% by state-of-the-art cytopathological molecular analyses (e.g., KRASmut DNA). Adding the simultaneous assay of proteins related to malignant transformation (e.g., MUC1 and CD55) is deemed essential to enhance diagnostic accuracy. The bioelectronic array proposed here, based on single-molecule-with-a-large-transistor (SiMoT) technology, can assay both nucleic acids and proteins at the single-molecule limit-of-identification (LOI) (1% of false-positives and false-negatives). It comprises an enzyme-linked immunosorbent assay (ELISA)-like 8 × 12-array organic-electronics disposable cartridge with an electrolyte-gated organic transistor sensor array, and a reusable reader, integrating a custom Si-IC chip, operating via software installed on a USB-connected smart device. The cartridge is complemented by a 3D-printed sensing gate cover plate. KRASmut , MUC1, and CD55 biomarkers either in plasma or cysts-fluid from 5 to 6 patients at a time, are multiplexed at single-molecule LOI in 1.5 h. The pancreatic cancer precursors are classified via a machine-learning analysis resulting in at least 96% diagnostic-sensitivity and 100% diagnostic-specificity. This preliminary study opens the way to POC liquid-biopsy-based early diagnosis of pancreatic-cancer precursors in plasma.

Practices, Potential, and Perspectives for Detecting Predisease Using Raman Spectroscopy

Raman spectroscopy shows great potential for practical clinical applications. By analyzing the structure and composition of molecules through real-time, non-destructive measurements of the scattered light from living cells and tissues, it offers valuable insights. The Raman spectral data directly link to the molecular composition of the cells and tissues and provides a "molecular fingerprint" for various disease states. This review focuses on the practical and clinical applications of Raman spectroscopy, especially in the early detection of human diseases. Identifying predisease, which marks the transition from a healthy to a disease state, is crucial for effective interventions to prevent disease onset. Raman spectroscopy can reveal biological processes occurring during the transition states and may eventually detect the molecular dynamics in predisease conditions.

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.

Extracellular matrix-related and serine protease proteins in the amniotic fluid of women with early preterm labor: Association with spontaneous preterm birth, intra-amniotic inflammation, and microbial invasion of the amniotic cavity

PROBLEM

We aimed to determine whether altered levels of various extracellular matrix (ECM)-related and serine protease proteins in the amniotic fluid (AF) are associated with imminent spontaneous preterm birth (SPTB; ≤7 days) and intra-amniotic inflammation and/or microbial invasion of the amniotic cavity (IAI/MIAC) in women with early preterm labor (PTL).

METHOD OF STUDY

This retrospective cohort study included 252 women with singleton pregnancies undergoing transabdominal amniocentesis who demonstrated PTL (24-31 weeks). The AF was cultured for microorganism detection to characterize MIAC. IL-6 concentrations were determined in the AF samples to identify IAI (≥2.6 ng/mL). The following mediators were measured in the AF samples using ELISA: kallistatin, lumican, MMP-2, SPARC, TGFBI, and uPA.

RESULTS

Kallistatin, MMP-2, TGFBI, and uPA levels were significantly higher and SPARC and lumican levels were significantly lower in the AF of women who spontaneously delivered within 7 days than in the AF of those who delivered after 7 days; the levels of the first five mediators were independent of baseline clinical variables. In the multivariate analysis, elevated levels of kallistatin, MMP-2, TGFBI, and uPA and low levels of lumican and SPARC in the AF were significantly associated with IAI/MIAC and MIAC, even after adjusting for the gestational age at sampling. The areas under the curves of the aforementioned biomarkers ranged from 0.58 to 0.87 for the diagnoses of each of the corresponding endpoints.

CONCLUSION

ECM-related (SPARC, TGFBI, lumican, and MMP-2) and serine protease (kallistatin and uPA) proteins in the AF are involved in preterm parturition and regulation of intra-amniotic inflammatory/infectious responses in PTL.

A patient-specific functional module and path identification technique from RNA-seq data

With the advancement of new technologies, a huge amount of high dimensional data is being generated which is opening new opportunities and challenges to the study of cancer and diseases. In particular, distinguishing the patient-specific key components and modules which drive tumorigenesis is necessary to analyze. A complex disease generally does not initiate from the dysregulation of a single component but it is the result of the dysfunction of a group of components and networks which differs from patient to patient. However, a patient-specific network is required to understand the disease and its molecular mechanism. We address this requirement by constructing a patient-specific network by sample-specific network theory with integrating cancer-specific differentially expressed genes and elite genes. By elucidating patient-specific networks, it can identify the regulatory modules, driver genes as well as personalized disease networks which can lead to personalized drug design. This method can provide insight into how genes are associating with each other and characterized the patient-specific disease subtypes. The results show that this method can be beneficial for the detection of patient-specific differential modules and interaction between genes. Extensive analysis using existing literature, gene enrichment and survival analysis for three cancer types STAD, PAAD and LUAD shows the effectiveness of this method over other existing methods. In addition, this method can be useful for personalized therapeutics and drug design. This methodology is implemented in the R language and is available at https://github.com/riasatazim/PatientSpecificRNANetwork.

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.

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.

Dynamic network biomarker identifies cdkn1a-mediated bone mineralization in the triggering phase of osteoporosis

The identification of predictive markers to determine the triggering phase prior to the onset of osteoporosis is essential to mitigate further irrevocable deterioration. To determine the early warning signs before osteoporosis, we used the dynamic network biomarker (DNB) approach to analyze time-series gene expression data in a zebrafish osteoporosis model, which revealed that cyclin-dependent kinase inhibitor 1 A (cdkn1a) is a core DNB. We found that cdkn1a negatively regulates osteogenesis, as evidenced by loss-of-function and gain-of-function studies. Specifically, CRISPR/Cas9-mediated cdkn1a knockout in zebrafish significantly altered skeletal development and increased bone mineralization, whereas inducible cdkn1a expression significantly contributed to osteoclast differentiation. We also found several mechanistic clues that cdkn1a participates in osteoclast differentiation by regulating its upstream signaling cascades. To summarize, in this study, we provided new insights into the dynamic nature of osteoporosis and identified cdkn1a as an early-warning signal of osteoporosis onset.

Gene filtering strategies for machine learning guided biomarker discovery using neonatal sepsis RNA-seq data

Machine learning (ML) algorithms are powerful tools that are increasingly being used for sepsis biomarker discovery in RNA-Seq data. RNA-Seq datasets contain multiple sources and types of noise (operator, technical and non-systematic) that may bias ML classification. Normalisation and independent gene filtering approaches described in RNA-Seq workflows account for some of this variability and are typically only targeted at differential expression analysis rather than ML applications. Pre-processing normalisation steps significantly reduce the number of variables in the data and thereby increase the power of statistical testing, but can potentially discard valuable and insightful classification features. A systematic assessment of applying transcript level filtering on the robustness and stability of ML based RNA-seq classification remains to be fully explored. In this report we examine the impact of filtering out low count transcripts and those with influential outliers read counts on downstream ML analysis for sepsis biomarker discovery using elastic net regularised logistic regression, L1-reguarlised support vector machines and random forests. We demonstrate that applying a systematic objective strategy for removal of uninformative and potentially biasing biomarkers representing up to 60% of transcripts in different sample size datasets, including two illustrative neonatal sepsis cohorts, leads to substantial improvements in classification performance, higher stability of the resulting gene signatures, and better agreement with previously reported sepsis biomarkers. We also demonstrate that the performance uplift from gene filtering depends on the ML classifier chosen, with L1-regularlised support vector machines showing the greatest performance improvements with our experimental data.

Label-Free Multimetric Measurement of Molecular Binding Kinetics by Electrical Modulation of a Flexible Nanobiolayer

Most label-free techniques rely on measuring refractive index or mass change on the sensor surface. Thus, it is challenging for them to measure small molecules or enzymatic processes that only induce a minor mass change on the analyte molecules. Here, we have developed a technique by combining Surface Plasmon Resonance sensing with an Oscillating Biomolecule Layer approach (SPR-OBL) to enhance the sensitivity of traditional SPR. In addition to the inherent mass sensitivity, SPR-OBL is also sensitive to the charge and conformational change of the analyte; hence it overcomes the mass limit and is able to detect small molecules. We show that the multimetric SPR-OBL measurement allows for sensing any changes regarding mass, charge, and conformation, which expands the detection capability of SPR.

Fluorescence biosensing of the leukemia gene by combining Target-Programmed controllable signal inspiring engineering

Clinical diagnosis urgently requires ultrasensitive, accurate and rapid monitoring of low-abundance biomarkers. A biosensing strategy capable of detecting target genes at the femtomolar scale was designed in this work. In the biosensing strategy, the target can induce the specially designed hairpin probe H1 to self-fold and form a 3' blunt-ended structure. When there are the hybrid double-stranded P1-T1, ligase, polymerase and nickase, the target gene was recycled, and at the same time the system produces a lot of T1 and T2. T1 and T2 can simultaneously trigger HCR, causing the modified fluorophore FAM on the DNA strand to move away from the quencher group BHQ. The amplified fluorescent signal can be captured by a fluorescence instrument. It is exciting for us that three signal amplifications are involved to achieve femtomolar detection of target genes, namely target recycling, dual-triggered HCR of T1 and T2, and HCR. In addition, it still has good detection ability in actual samples simulated by serum. We expect that the sensing strategy proposed in this paper offers great potential for biomarker detection of leukemia for early clinical diagnosis.

Exploring drought-responsive crucial genes in Sorghum

Drought severely affects global food production. Sorghum is a typical drought-resistant model crop. Based on RNA-seq data for Sorghum with multiple time points and the gray correlation coefficient, this paper firstly selects candidate genes via mean variance test and constructs weighted gene differential co-expression networks (WGDCNs); then, based on guilt-by-rewiring principle, the WGDCNs and the hidden Markov random field model, drought-responsive crucial genes are identified for five developmental stages respectively. Enrichment and sequence alignment analysis reveal that the screened genes may play critical functional roles in drought responsiveness. A multilayer differential co-expression network for the screened genes reveals that Sorghum is very sensitive to pre-flowering drought. Furthermore, a crucial gene regulatory module is established, which regulates drought responsiveness via plant hormone signal transduction, MAPK cascades, and transcriptional regulations. The proposed method can well excavate crucial genes through RNA-seq data, which have implications in breeding of new varieties with improved drought tolerance.

•

We design a method that unites gene rewiring network and Markov random field model

•

Drought-responsive genes for five developmental stages of Sorghum are explored

•

A multilayer network reveals that Sorghum is very sensitive to pre-flowering drought

•

A drought-responsive crucial gene regulatory module is established for Sorghum

Molecular biomarkers, network biomarkers, and dynamic network biomarkers for diagnosis and prediction of rare diseases

The difficulty of converting scientific research findings into novel pharmacological treatments for rare and life-threatening diseases is enormous. Biomarkers related to multiple biological processes involved in cell growth, proliferation, and disease occurrence have been identified in recent years with the development of immunology, molecular biology, and genomics technologies. Biomarkers are capable of reflecting normal physiological processes, pathological processes, and the response to therapeutic intervention; as such, they play vital roles in disease diagnosis, prevention, drug response, and other aspects of biomedicine. The discovery of valuable biomarkers has become a focal point of current research. Numerous studies have identified molecular biomarkers based on the differential expression/concentration of molecules (e.g., genes/proteins) for disease state diagnosis, characterization, and treatment. Although technological breakthroughs in molecular analysis platforms have enabled the identification of a large number of candidate biomarkers for rare diseases, only a small number of these candidates have been properly validated for use in patient treatment. The traditional molecular biomarkers may lose vital information by ignoring molecular associations/interactions, and thus the concept of network biomarkers based on differential associations/correlations of molecule pairs has been established. This approach promises to be more stable and reliable in diagnosing disease states. Furthermore, the newly-emerged dynamic network biomarkers (DNBs) based on differential fluctuations/correlations of molecular groups are able to recognize pre-disease states or critical states instead of disease states, thereby achieving rare disease prediction or predictive/preventative medicine and providing deep insight into the dynamic characteristics of disease initiation and progression.

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.

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.

Interpretable machine learning-derived nomogram model for early detection of diabetic retinopathy in type 2 diabetes mellitus: a widely targeted metabolomics study

Objective

Early identification of diabetic retinopathy (DR) is key to prioritizing therapy and preventing permanent blindness. This study aims to propose a machine learning model for DR early diagnosis using metabolomics and clinical indicators.

Methods

From 2017 to 2018, 950 participants were enrolled from two affiliated hospitals of Wenzhou Medical University and Anhui Medical University. A total of 69 matched blocks including healthy volunteers, type 2 diabetes, and DR patients were obtained from a propensity score matching-based metabolomics study. UPLC-ESI-MS/MS system was utilized for serum metabolic fingerprint data. CART decision trees (DT) were used to identify the potential biomarkers. Finally, the nomogram model was developed using the multivariable conditional logistic regression models. The calibration curve, Hosmer–Lemeshow test, receiver operating characteristic curve, and decision curve analysis were applied to evaluate the performance of this predictive model.

Results

The mean age of enrolled subjects was 56.7 years with a standard deviation of 9.2, and 61.4% were males. Based on the DT model, 2-pyrrolidone completely separated healthy controls from diabetic patients, and thiamine triphosphate (ThTP) might be a principal metabolite for DR detection. The developed nomogram model (including diabetes duration, systolic blood pressure and ThTP) shows an excellent quality of classification, with AUCs (95% CI) of 0.99 (0.97–1.00) and 0.99 (0.95–1.00) in training and testing sets, respectively. Furthermore, the predictive model also has a reasonable degree of calibration.

Conclusions

The nomogram presents an accurate and favorable prediction for DR detection. Further research with larger study populations is needed to confirm our findings.

Transcriptional progressive patterns from mild to severe renal ischemia/reperfusion-induced kidney injury in mice

The renal ischemia/reperfusion (I/R)-induced acute kidney injury incidence after nephron-sparing surgery for localized renal tumors is 20%, but the biological determinant process of postoperative acute kidney injury remains unclear. Using Gene Expression Omnibus database (GSE192883) and several bioinformatics analyses (discrete time points analysis, gene set enrichment analysis, dynamic network biomarker analysis, etc), combined with the establishment of the I/R model for verification, we identified three progressive patterns involving five core pathways confirmed using gene set enrichment analysis and six key genes (S100a10, Pcna, Abat, Kmo, Acadm, and Adhfe1) verified using quantitative polymerase chain reaction The dynamic network biomarker (DNB) subnetwork composite index value is the highest in the 22-min ischemia group, suggesting the transcriptome expression level fluctuated sharply in this group, which means 22-min ischemia is an critical warning point. This study illustrates the core molecular progressive patterns from mild to severe I/R kidney injury, laying the foundation for precautionary biomarkers and molecular intervention targets for exploration. In addition, the safe renal artery blocking time of nephron-sparing surgery that we currently accept may not be safe anymore.

A Dual Ranking Algorithm Based on the Multiplex Network for Heterogeneous Complex Disease Analysis

Identifying biomarkers of heterogeneous complex diseases has always been one of the focuses in medical research. In previous studies, the powerful network propagation methods have been applied to finding marker genes related to specific diseases, but existing methods are mostly based on a single network, which may be greatly affected by the incompleteness of the network and the ignorance of a large amount of information about physical and functional interactions between biological components. Other methods that directly integrate multiple types of interactions into an aggregate network have the risks that different types of data may conflict with each other and the characteristics and topologies of each individual network are lost. Meanwhile, biomarkers used in clinical trials should have the characteristics of small quantity and strong discriminate ability. In this study, we developed a multiplex network-based dual ranking framework (DualRank) for heterogeneous complex disease analysis. We applied the proposed method to heterogeneous complex diseases for diagnosis, prognosis, and classification. The results showed that DualRank outperformed competing methods and could identify biomarkers with the small quantity, great prediction performance (average AUC = 0.818) and biological interpretability.

Identification of Key Genes as Early Warning Signals of Acute Myocardial Infarction Based on Weighted Gene Correlation Network Analysis and Dynamic Network Biomarker Algorithm

Purpose

The specific mechanisms and biomarkersunderlying the progression of stable coronary artery disease (CAD) to acute myocardial infarction (AMI) remain unclear. The current study aims to explore novel gene biomarkers associated with CAD progression by analyzing the transcriptomic sequencing data of peripheral blood monocytes in different stages of CAD.

Material and Methods

A total of 24 age- and sex- matched patients at different CAD stages who received coronary angiography were enrolled, which included 8 patients with normal coronary angiography, 8 patients with angiographic intermediate lesion, and 8 patients with AMI. The RNA from peripheral blood monocytes was extracted and transcriptome sequenced to analyze the gene expression and the differentially expressed genes (DEG). A Gene Oncology (GO) enrichment analysis was performed to analyze the biological function of genes. Weighted gene correlation network analysis (WGCNA) was performed to classify genes into several gene modules with similar expression profiles, and correlation analysis was carried out to explore the association of each gene module with a clinical trait. The dynamic network biomarker (DNB) algorithm was used to calculate the key genes that promote disease progression. Finally, the overlapping genes between different analytic methods were explored.

Results

WGCNA analysis identified a total of nine gene modules, of which two modules have the highest positive association with CAD stages. GO enrichment analysis indicated that the biological function of genes in these two gene modules was closely related to inflammatory response, which included T-cell activation, cell response to inflammatory stimuli, lymphocyte activation, cytokine production, and the apoptotic signaling pathway. DNB analysis identified a total of 103 genes that may play key roles in the progression of atherosclerosis plaque. The overlapping genes between DEG/WGCAN and DNB analysis identified the following 13 genes that may play key roles in the progression of atherosclerosis disease: SGPP2, DAZAP2, INSIG1, CD82, OLR1, ARL6IP1, LIMS1, CCL5, CDK7, HBP1, PLAU, SELENOS, and DNAJB6.

Conclusions

The current study identified a total of 13 genes that may play key roles in the progression of atherosclerotic plaque and provides new insights for early warning biomarkers and underlying mechanisms underlying the progression of CAD.

Identifying network biomarkers of cancer by sample-specific differential network

Abundant datasets generated from various big science projects on diseases have presented great challenges and opportunities, which contributed to unfolding the complexity of diseases. The discovery of disease-associated molecular networks for each individual plays an important role in personalized therapy and precision treatment of cancer-based on the reference networks. However, there are no effective ways to distinguish the consistency of different reference networks. In this study, we developed a statistical method, i.e. a sample-specific differential network (SSDN), to construct and analyze such networks based on gene expression of a single sample against a reference dataset. We proved that the SSDN is structurally consistent even with different reference datasets if the reference dataset can follow certain conditions. The SSDN also can be used to identify patient-specific disease modules or network biomarkers as well as predict the potential driver genes of a tumor sample.

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.

Critical Parameters in Dynamic Network Modeling of Sepsis

In this work, we propose a dynamical systems perspective on the modeling of sepsis and its organ-damaging consequences. We develop a functional two-layer network model for sepsis based upon the interaction of parenchymal cells and immune cells via cytokines, and the coevolutionary dynamics of parenchymal, immune cells, and cytokines. By means of the simple paradigmatic model of phase oscillators in a two-layer system, we analyze the emergence of organ threatening interactions between the dysregulated immune system and the parenchyma. We demonstrate that the complex cellular cooperation between parenchyma and stroma (immune layer) either in the physiological or in the pathological case can be related to dynamical patterns of the network. In this way we explain sepsis by the dysregulation of the healthy homeostatic state (frequency synchronized) leading to a pathological state (desynchronized or multifrequency cluster) in the parenchyma. We provide insight into the complex stabilizing and destabilizing interplay of parenchyma and stroma by determining critical interaction parameters. The coupled dynamics of parenchymal cells (metabolism) and nonspecific immune cells (response of the innate immune system) is represented by nodes of a duplex layer. Cytokine interaction is modeled by adaptive coupling weights between nodes representing immune cells (with fast adaptation timescale) and parenchymal cells (slow adaptation timescale), and between pairs of parenchymal and immune cells in the duplex network (fixed bidirectional coupling). The proposed model allows for a functional description of organ dysfunction in sepsis and the recurrence risk in a plausible pathophysiological context.

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.

Clusterin as a serum biomarker candidate contributes to the lung fibroblasts activation in chronic obstructive pulmonary disease

BACKGROUND

Fibrosis in the peripheral airways contributes to airflow limitation in patients with chronic obstructive pulmonary disease (COPD). However, the key proteins involved in its development are still poorly understood. Thus, we aimed to identify the differentially expressed proteins (DEPs) between smoker patients with and without COPD and elucidate the molecular mechanisms involved by investigating the effects of the identified biomarker candidate on lung fibroblasts.

METHODS

The potential DEPs were identified by isobaric tags for relative and absolute quantitation (iTRAQ)-based proteomic analysis. The messenger RNA and protein levels of clusterin (CLU) in COPD patients and 12% cigarette smoke extract (CSE)-treated human bronchial epithelial cells were determined at the indicated time points. Furthermore, an in vitro COPD model was established via the administration of 8% CSE to normal human lung fibroblasts (NHLFs) at indicated time points. The effects of CSE treatment and CLU silencing on proliferation and activation of lung fibroblasts were analyzed.

RESULTS

A total of 144 DEPs were identified between COPD patients and normal smokers. The iTRAQ-based proteomics and bioinformatics analyses identified CLU as a serum biomarker candidate. We also discovered that CLU levels were significantly increased ( P   <  0.0001) in Global Initiative for Obstructive Lung Disease II, III, and IV patients and correlated ( P   <  0.0001) with forced expiratory volume in 1 s ( R  = -0.7705), residual volume (RV) ( R  = 0.6281), RV/total lung capacity ( R  = 0.5454), and computerized tomography emphysema ( R  = 0.7878). Similarly, CLU levels were significantly increased in CSE-treated cells at indicated time points ( P  < 0.0001). The CSE treatment significantly inhibited the proliferation, promoted the inflammatory response, differentiation of NHLFs, and collagen matrix deposition, and induced the apoptosis of NHLFs; however, these effects were partially reversed by CLU silencing.

CONCLUSION

Our findings suggest that CLU may play significant roles during airway fibrosis in COPD by regulating lung fibroblast activation.

Molecular Subtyping of Cancer Based on Distinguishing Co-Expression Modules and Machine Learning

Molecular subtyping of cancer is recognized as a critical and challenging step towards individualized therapy. Most existing computational methods solve this problem via multi-classification of gene-expressions of cancer samples. Although these methods, especially deep learning, perform well in data classification, they usually require large amounts of data for model training and have limitations in interpretability. Besides, as cancer is a complex systemic disease, the phenotypic difference between cancer samples can hardly be fully understood by only analyzing single molecules, and differential expression-based molecular subtyping methods are reportedly not conserved. To address the above issues, we present here a new framework for molecular subtyping of cancer through identifying a robust specific co-expression module for each subtype of cancer, generating network features for each sample by perturbing correlation levels of specific edges, and then training a deep neural network for multi-class classification. When applied to breast cancer (BRCA) and stomach adenocarcinoma (STAD) molecular subtyping, it has superior classification performance over existing methods. In addition to improving classification performance, we consider the specific co-expressed modules selected for subtyping to be biologically meaningful, which potentially offers new insight for diagnostic biomarker design, mechanistic studies of cancer, and individualized treatment plan selection.

Determination of tipping point in course of PM2.5 organic extracts-induced malignant transformation by dynamic network biomarkers

The dynamic network biomarkers (DNBs) are designed to identify the tipping point and specific molecules in initiation of PM_2.5-induced lung cancers. To discover early-warning signals, we analyzed time-series gene expression datasets over a course of PM_2.5 organic extraction-induced human bronchial epithelial (HBE) cell transformation (0^th~16^th week). A composition index of DNB (CI_DNB) was calculated to determine correlations and fluctuations in molecule clusters at each timepoint. We identified a group of genes with the highest CI_DNB at the 10^th week, implicating a tipping point and corresponding DNBs. Functional experiments revealed that manipulating respective DNB genes at the tipping point led to remarkable changes in malignant phenotypes, including four promoters (GAB2, NCF1, MMP25, LAPTM5) and three suppressors (BATF2, DOK3, DAP3). Notably, co-altered expression of seven core DNB genes resulted in an enhanced activity of malignant transformation compared to effects of single-gene manipulation. Perturbation of pathways (EMT, HMGB1, STAT3, NF-κB, PTEN) appeared in HBE cells at the tipping point. The core DNB genes were involved in regulating lung cancer cell growth and associated with poor survival, indicating their synergistic effects in initiation and development of lung cancers. These findings provided novel insights into the mechanism of dynamic networks attributable to PM_2.5-induced cell transformation.

Tipping the Balance: A Criticality Perspective

Cell populations are often characterised by phenotypic heterogeneity in the form of two distinct subpopulations. We consider a model of tumour cells consisting of two subpopulations: non-cancer promoting (NCP) and cancer-promoting (CP). Under steady state conditions, the model has similarities with a well-known model of population genetics which exhibits a purely noise-induced transition from unimodality to bimodality at a critical value of the noise intensity σ2. The noise is associated with the parameter λ representing the system-environment coupling. In the case of the tumour model, λ has a natural interpretation in terms of the tissue microenvironment which has considerable influence on the phenotypic composition of the tumour. Oncogenic transformations give rise to considerable fluctuations in the parameter. We compute the λ−σ2 phase diagram in a stochastic setting, drawing analogies between bifurcations and phase transitions. In the region of bimodality, a transition from a state of balance to a state of dominance, in terms of the competing subpopulations, occurs at λ = 0. Away from this point, the NCP (CP) subpopulation becomes dominant as λ changes towards positive (negative) values. The variance of the steady state probability density function as well as two entropic measures provide characteristic signatures at the transition point.

A Novel Method for Constructing Classification Models by Combining Different Biomarker Patterns

Different biomarker patterns, such as those of molecular biomarkers and ratio biomarkers, have their own merits in clinical applications. In this study, a novel machine learning method used in biomedical data analysis for constructing classification models by combining different biomarker patterns (CDBP)is proposed. CDBP uses relative expression reversals to measure the discriminative ability of different biomarker patterns, and selects the pattern with the higher score for classifier construction. The decision boundary of CDBP can be characterized in simple and biologically meaningful manners. The CDBP method was compared with eight state-of-the-art methods on eight gene expression datasets to test its performance. CDBP, with fewer features or ratio features, had the highest classification performance. Subsequently, CDBP was employed to extract crucial diagnostic information from a rat hepatocarcinogenesis metabolomics dataset. The potential biomarkers selected by CDBP provided better classification of hepatocellular carcinoma (HCC)and non-HCC stages than previous works in the animal model. The statistical analyses of these potential biomarkers in an independent human dataset confirmed their discriminative abilities of different liver diseases. These experimental results highlight the potential of CDBP for biomarker identification from high-dimensional biomedical datasets and demonstrate that it can be a useful tool for disease classification.