Dynamic network biomarker indicates pulmonary metastasis at the tipping point of hepatocellular carcinoma

Time-series-analysis-based detection of critical transitions in real-world non-autonomous systems

Real-world non-autonomous systems are open, out-of-equilibrium systems that evolve in and are driven by temporally varying environments. Such systems can show multiple timescale and transient dynamics together with transitions to very different and, at times, even disastrous dynamical regimes. Since such critical transitions disrupt the systems' intended or desired functionality, it is crucial to understand the underlying mechanisms, to identify precursors of such transitions, and to reliably detect them in time series of suitable system observables to enable forecasts. This review critically assesses the various steps of investigation involved in time-series-analysis-based detection of critical transitions in real-world non-autonomous systems: from the data recording to evaluating the reliability of offline and online detections. It will highlight pros and cons to stimulate further developments, which would be necessary to advance understanding and forecasting nonlinear behavior such as critical transitions in complex systems.

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.

Interleukin-25 as a Potential Biomarker in Lung Metastasis of Hepatocellular Carcinoma with HBV History in Chinese Patients: A Single Center, Case-control Study

Background: Interleukin-25 (IL-25) has been proved to play a role in the pathogenesis and metastasis of Hepatocellular carcinoma (HCC), but the relationship between the level of IL-25 and the metastasis and prognosis of HCC is still not clear. This study aimed to investigate the expression of IL-25 and other potential biochemical indicators among healthy people, HBV-associated HCC patients without lung metastasis and HBV-associated HCC patients with lung metastasis. Methods: From September 2019 to November 2021, 33 HCC patients without lung metastasis, 37 HCC patients with lung metastasis and 29 healthy controls were included in the study. IL-25 and other commonly used biochemical markers were measured to establish predictors of overall survival (OS) and progression-free survival (PFS) after treatment. Results: The serum level of IL-25 was increased in HCC patients than healthy controls (p < 0.001) and HCC patients with lung metastasis had higher IL-25 level than HCC patients without metastasis (p = 0.035). Lung metastasis also indicated higher death rate (p < 0.001) by chi-square test, higher GGT level (p = 0.024) and higher AFP level (p = 0.049) by non-parametric test. Kaplan-Meier analysis demonstrated that IL-25 was negatively associated with PFS (p = 0.024). Multivariate Cox-regression analysis indicated IL-25 (p = 0.030) and GGT (p = 0.020) to be independent predictors of poorer PFS, while IL-25 showed no significant association with OS. Conclusion: The level of IL-25 was significantly associated with disease progression and lung metastasis of HBV-associated HCC. The high expression of IL-25 predicted high recurrence rate and death probability of HCC patients after treatment. Therefore, IL-25 may be an effective predictor of prognosis in HCC.

Exploring MiR-484 Regulation by Polyalthia longifolia: A Promising Biomarker and Therapeutic Target in Cervical Cancer through Integrated Bioinformatics and an In Vitro Analysis

BACKGROUND

MiR-484, implicated in various carcinomas, holds promise as a prognostic marker, yet its relevance to cervical cancer (CC) remains unclear. Our prior study demonstrated the Polyalthia longifolia downregulation of miR-484, inhibiting HeLa cells. This study investigates miR-484's potential as a biomarker and therapeutic target in CC through integrated bioinformatics and an in vitro analysis.

METHODS

MiR-484 levels were analyzed across cancers, including CC, from The Cancer Genome Atlas. The limma R package identified differentially expressed genes (DEGs) between high- and low-miR-484 CC cohorts. We assessed biological functions, tumor microenvironment (TME), immunotherapy, stemness, hypoxia, RNA methylation, and chemosensitivity differences. Prognostic genes relevant to miR-484 were identified through Cox regression and Kaplan-Meier analyses, and a prognostic model was captured via multivariate Cox regression. Single-cell RNA sequencing determined cell populations related to prognostic genes. qRT-PCR validated key genes, and the miR-484 effect on CC proliferation was assessed via an MTT assay.

RESULTS

MiR-484 was upregulated in most tumors, including CC, with DEGs enriched in skin development, PI3K signaling, and immune processes. High miR-484 expression correlated with specific immune cell infiltration, hypoxia, and drug sensitivity. Prognostic genes identified were predominantly epidermal and stratified patients with CC into risk groups, with the low-risk group showing enhanced survival and immunotherapeutic responses. qRT-PCR confirmed FGFR3 upregulation in CC cells, and an miR-484 mimic reversed the P. longifolia inhibitory effect on HeLa proliferation.

CONCLUSION

MiR-484 plays a crucial role in the CC progression and prognosis, suggesting its potential as a biomarker for targeted therapy.

Kinesin Family Member-18A (KIF18A) Promotes Cell Proliferation and Metastasis in Hepatocellular Carcinoma

BACKGROUND & AIMS

Kinesin family member 18A (KIF18A) is notable for its aberrant expression across various cancer types and its pivotal role is driving cancer progression. In this study, we aim to investigate the intricate molecular mechanisms underlying the impact of KIF18A on the progression of HCC.

METHODS

Western blotting assays, a quantitative real-time PCR and immunohistochemical analyses were performed to quantitatively assess KIF18A expression in HCC tissues. We then performed genetic manipulations within HCC cells by silencing endogenous KIF18A using short hairpin RNA (shRNA) and introducing exogenous plasmids to overexpress KIF18A. We monitored cell progression, analyzed cell cycle and cell apoptosis and assessed cell migration and invasion both in vitro and in vivo. Moreover, we conducted RNA-sequencing to explore KIF18A-related signaling pathways utilizing Reactome and KEGG enrichment methods and validated these critical mediators in these pathways.

RESULTS

Analysis of the TCGA-LIHC database revealed pronounced overexpression of KIF18A in HCC tissues, the finding was subsequently confirmed through the analysis of clinical samples obtained from HCC patients. Notably, silencing KIF18A in cells led to an obvious inhibition of cell proliferation, migration and invasion in vitro. Furthermore, in subcutaneous and orthotopic xenograft models, suppression of KIF18A sgnificantly redudce tumor weight and the number of lung metastatic nodules. Mechanistically, KIF18A appears to facilitate cell proliferation by upregulating MAD2 and CDK1/CyclinB1 expression levels, with the activation of SMAD2/3 signaling contributing to KIF18A-driven metastasis.

CONCLUSION

Our study elucidates the molecular mechanism by which KIF18A mediates proliferation and metastasis in HCC cells, offering new insights into potential therapeutic targets.

Dysfunction of CCT3-associated network signals for the critical state during progression of hepatocellular carcinoma

Hepatocellular carcinoma (HCC) is one of the most common malignant tumors and is a serious threat to human health; thus, early diagnosis and adequate treatment are essential. However, there are still great challenges in identifying the tipping point and detecting early warning signals of early HCC. In this study, we aimed to identify the tipping point (critical state) of and key molecules involved in hepatocarcinogenesis based on time series transcriptome expression data of HCC patients. The phase from veHCC (very early HCC) to eHCC (early HCC) was identified as the critical state in HCC progression, with 143 genes identified as key candidate molecules by combining the DDRTree (dimensionality reduction via graph structure learning) and DNB (dynamic network biomarker) methods. Then, we ranked the candidate genes to verify their mRNA levels using the diethylnitrosamine (DEN)-induced HCC mouse model and identified five early warning signals, namely, CCT3, DSTYK, EIF3E, IARS2 and TXNRD1; these signals can be regarded as the potential early warning signals for the critical state of HCC. We identified CCT3 as an independent prognostic factor for HCC, and functions of CCT3 involving in the "MYCtargets_V1" and "E2F-Targets" are closely related to the progression of HCC. The predictive method combining the DDRTree and DNB methods can not only identify the key critical state before cancer but also determine candidate molecules of critical state, thus providing new insight into the early diagnosis and preemptive treatment of HCC.

Establishing Monoclonal Gammopathy of Undetermined Significance as an Independent Pre-Disease State of Multiple Myeloma Using Raman Spectroscopy, Dynamical Network Biomarker Theory, and Energy Landscape Analysis

Multiple myeloma (MM) is a cancer of plasma cells. Normal (NL) cells are considered to pass through a precancerous state, such as monoclonal gammopathy of undetermined significance (MGUS), before transitioning to MM. In the present study, we acquired Raman spectra at three stages-834 NL, 711 MGUS, and 970 MM spectra-and applied the dynamical network biomarker (DNB) theory to these spectra. The DNB analysis identified MGUS as the unstable pre-disease state of MM and extracted Raman shifts at 1149 and 1527-1530 cm-1 as DNB variables. The distribution of DNB scores for each patient showed a significant difference between the mean values for MGUS and MM patients. Furthermore, an energy landscape (EL) analysis showed that the NL and MM stages were likely to become stable states. Raman spectroscopy, the DNB theory, and, complementarily, the EL analysis will be applicable to the identification of the pre-disease state in clinical samples.

Rate-induced tipping in complex high-dimensional ecological networks

In an ecosystem, environmental changes as a result of natural and human processes can cause some key parameters of the system to change with time. Depending on how fast such a parameter changes, a tipping point can occur. Existing works on rate-induced tipping, or R-tipping, offered a theoretical way to study this phenomenon but from a local dynamical point of view, revealing, e.g., the existence of a critical rate for some specific initial condition above which a tipping point will occur. As ecosystems are subject to constant disturbances and can drift away from their equilibrium point, it is necessary to study R-tipping from a global perspective in terms of the initial conditions in the entire relevant phase space region. In particular, we introduce the notion of the probability of R-tipping defined for initial conditions taken from the whole relevant phase space. Using a number of real-world, complex mutualistic networks as a paradigm, we find a scaling law between this probability and the rate of parameter change and provide a geometric theory to explain the law. The real-world implication is that even a slow parameter change can lead to a system collapse with catastrophic consequences. In fact, to mitigate the environmental changes by merely slowing down the parameter drift may not always be effective: Only when the rate of parameter change is reduced to practically zero would the tipping be avoided. Our global dynamics approach offers a more complete and physically meaningful way to understand the important phenomenon of R-tipping.

Multiomics characterization of fatty acid metabolism for the clinical management of hepatocellular carcinoma

Hepatocellular carcinoma (HCC) is a prevalent malignancy and there is a lack of effective biomarkers for HCC diagnosis. Living organisms are complex, and different omics molecules interact with each other to implement various biological functions. Genomics and metabolomics, which are the top and bottom of systems biology, play an important role in HCC clinical management. Fatty acid metabolism is associated with malignancy, prognosis, and immune phenotype in cancer, which is a potential hallmark in malignant tumors. In this study, the genes and metabolites related to fatty acid metabolism were thoroughly investigated by a dynamic network construction algorithm named EWS-DDA for the early diagnosis and prognosis of HCC. Three gene ratios and eight metabolite ratios were identified by EWS-DDA as potential biomarkers for HCC clinical management. Further analysis using biological analysis, statistical analysis and document validation in the discovery and validation sets suggested that the selected potential biomarkers had great clinical prognostic value and helped to achieve effective early diagnosis of HCC. Experimental results suggested that in-depth evaluation of fatty acid metabolism from different omics viewpoints can facilitate the further understanding of pathological alterations associated with HCC characteristics, improving the performance of early diagnosis and clinical prognosis.

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.

Repurposing fluphenazine as an autophagy modulator for treating liver cancer

Hepatocellular carcinoma (HCC) is a common malignant tumor of the digestive system with a low early diagnosis rate. Owing to the side effects, tolerance, and patient contraindications of existing therapies, effective drug treatments for HCC remain a major clinical challenge. However, using approved or investigational drugs not initially intended for cancer therapy is a promising strategy for resolving this problem because their safety have been tested in clinic. Therefore, this study evaluated differentially expressed genes between liver cancer and normal tissues in a cohort of patients with HCC from The Cancer Genome Atlas and applied them to query a connectivity map to identify candidate anti-HCC drugs. As a result, fluphenazine was identified as a candidate for anti-HCC therapy in vitro and in vivo. Fluphenazine suppressed HCC cell proliferation and migration and induced cell cycle arrest and apoptosis, possibly owing to disrupted lysosomal function, blocking autophagy flux. Additionally, in vivo studies demonstrated that fluphenazine suppresses HCC subcutaneous xenografts growth without causing severe side effects. Strikingly, fluphenazine could be used as an analgesic to alleviate oxaliplatin-induced pain as well as pain related anxiety-like behavior. Therefore, fluphenazine could be a novel liver cancer treatment candidate.

Identifying the potential miRNA biomarkers based on multi-view networks and reinforcement learning for diseases

MicroRNAs (miRNAs) play important roles in the occurrence and development of diseases. However, it is still challenging to identify the effective miRNA biomarkers for improving the disease diagnosis and prognosis. In this study, we proposed the miRNA data analysis method based on multi-view miRNA networks and reinforcement learning, miRMarker, to define the potential miRNA disease biomarkers. miRMarker constructs the cooperative regulation network and functional similarity network based on the expression data and known miRNA-disease relations, respectively. The cooperative regulation of miRNAs was evaluated by measuring the changes of relative expression. Natural language processing was introduced for calculating the miRNA functional similarity. Then, miRMarker integrates the multi-view miRNA networks and defines the informative miRNA modules through a reinforcement learning strategy. We compared miRMarker with eight efficient data analysis methods on nine transcriptomics datasets to show its superiority in disease sample discrimination. The comparison results suggested that miRMarker outperformed other data analysis methods in receiver operating characteristic analysis. Furthermore, the defined miRNA modules of miRMarker on colorectal cancer data not only show the excellent performance of cancer sample discrimination but also play significant roles in the cancer-related pathway disturbances. The experimental results indicate that miRMarker can build the robust miRNA interaction network by integrating the multi-view networks. Besides, exploring the miRNA interaction network using reinforcement learning favors defining the important miRNA modules. In summary, miRMarker can be a hopeful tool in biomarker identification for human diseases.

Trehalose metabolism coordinates transcriptional regulatory control and metabolic requirements to trigger the onset of cassava storage root initiation

Cassava storage roots (SR) are an important source of food energy and raw material for a wide range of applications. Understanding SR initiation and the associated regulation is critical to boosting tuber yield in cassava. Decades of transcriptome studies have identified key regulators relevant to SR formation, transcriptional regulation and sugar metabolism. However, there remain uncertainties over the roles of the regulators in modulating the onset of SR development owing to the limitation of the widely applied differential gene expression analysis. Here, we aimed to investigate the regulation underlying the transition from fibrous (FR) to SR based on Dynamic Network Biomarker (DNB) analysis. Gene expression analysis during cassava root initiation showed the transition period to SR happened in FR during 8 weeks after planting (FR8). Ninety-nine DNB genes associated with SR initiation and development were identified. Interestingly, the role of trehalose metabolism, especially trehalase1 (TRE1), in modulating metabolites abundance and coordinating regulatory signaling and carbon substrate availability via the connection of transcriptional regulation and sugar metabolism was highlighted. The results agree with the associated DNB characters of TRE1 reported in other transcriptome studies of cassava SR initiation and Attre1 loss of function in literature. The findings help fill the knowledge gap regarding the regulation underlying cassava SR initiation.

Review of T cell proliferation regulatory factors in treatment and prognostic prediction for solid tumors

T cell proliferation regulators (Tcprs), which are positive regulators that promote T cell function, have made great contributions to the development of therapies to improve T cell function. CAR (chimeric antigen receptor) -T cell therapy, a type of adoptive cell transfer therapy that targets tumor cells and enhances immune lethality, has led to significant progress in the treatment of hematologic tumors. However, the applications of CAR-T in solid tumor treatment remain limited. Therefore, in this review, we focus on the development of Tcprs for solid tumor therapy and prognostic prediction. We summarize potential strategies for targeting different Tcprs to enhance T cell proliferation and activation and inhibition of cancer progression, thereby improving the antitumor activity and persistence of CAR-T. In summary, we propose means of enhancing CAR-T cells by expressing different Tcprs, which may lead to the development of a new generation of cell therapies.

Dynamical hallmarks of cancer: Phenotypic switching in melanoma and epithelial-mesenchymal plasticity

Phenotypic plasticity was recently incorporated as a hallmark of cancer. This plasticity can manifest along many interconnected axes, such as stemness and differentiation, drug-sensitive and drug-resistant states, and between epithelial and mesenchymal cell-states. Despite growing acceptance for phenotypic plasticity as a hallmark of cancer, the dynamics of this process remains poorly understood. In particular, the knowledge necessary for a predictive understanding of how individual cancer cells and populations of cells dynamically switch their phenotypes in response to the intensity and/or duration of their current and past environmental stimuli remains far from complete. Here, we present recent investigations of phenotypic plasticity from a systems-level perspective using two exemplars: epithelial-mesenchymal plasticity in carcinomas and phenotypic switching in melanoma. We highlight how an integrated computational-experimental approach has helped unravel insights into specific dynamical hallmarks of phenotypic plasticity in different cancers to address the following questions: a) how many distinct cell-states or phenotypes exist?; b) how reversible are transitions among these cell-states, and what factors control the extent of reversibility?; and c) how might cell-cell communication be able to alter rates of cell-state switching and enable diverse patterns of phenotypic heterogeneity? Understanding these dynamic features of phenotypic plasticity may be a key component in shifting the paradigm of cancer treatment from reactionary to a more predictive, proactive approach.

Dynamic gene regulatory networks improving spike fertility through regulation of floret primordia fate in wheat

The developmental process of spike is critical for spike fertility through affecting floret primordia fate in wheat; however, the genetic regulation of this dynamic and complex developmental process remains unclear. Here, we conducted a high temporal-resolution analysis of spike transcriptomes and monitored the number and morphology of floret primordia within spike. The development of all floret primordia in a spike was clearly separated into three distinct phases: differentiation, pre-dimorphism and dimorphism. Notably, we identified that floret primordia with meiosis ability at the pre-dimorphism phase usually develop into fertile floret primordia in the next dimorphism phase. Compared to control, increasing plant space treatment achieved the maximum increasement range (i.e., 50%) in number of fertile florets by accelerating spike development. The process of spike fertility improvement was directed by a continuous and dynamic regulatory network involved in transcription factor and genes interaction. This was based on the coordination of genes related to heat shock protein and jasmonic acid biosynthesis during differentiation phase, and genes related to lignin, anthocyanin and chlorophyll biosynthesis during dimorphism phase. The multi-dimensional association with high temporal-resolution approach reported here allows rapid identification of genetic resource for future breeding studies to realise the maximum spike fertility potential in more cereal crops.

A clinically practical model for the preoperative prediction of lymph node metastasis in bladder cancer: a multicohort study

BACKGROUND

The aim of this study was to construct a clinically practical model to precisely predict lymph node (LN) metastasis in bladder cancer patients.

METHODS

Four independent cohorts were included. The least absolute shrinkage and selection operator regression with multivariate logistic regression were applied. The diagnostic efficacy of LN score and CT/MRI was compared by accuracy, sensitivity, specificity, and area under curve (AUC).

RESULTS

A total of 606 patients were included to develop a basic prediction model. After multistep gene selection, the LN metastasis prediction model was constructed with 5 genes. The model can accurately predict LN metastasis with an AUC of 0.781. For clinically practical use, we transformed the model into a Fast LN Scoring System using the SYSMH cohort (n = 105). High LN score patients exhibited a 72.2% LN metastasis rate, while low LN score patients showed a 3.4% LN metastasis rate. The LN score achieved a superior accuracy than CT/MRI (0.882 vs. 0.727). Application of LN score can correct the diagnosis of 88% (22/25) patients who were misdiagnosed by CT/MRI.

DISCUSSION

The clinically practical LN score can precisely, rapidly, and conveniently predict LN status, which will assist preoperative diagnosis for LN metastasis and guide precise therapy.

New Possibilities for Evaluating the Development of Age-Related Pathologies Using the Dynamical Network Biomarkers Theory

Aging is the slowest process in a living organism. During this process, mortality rate increases exponentially due to the accumulation of damage at the cellular level. Cellular senescence is a well-established hallmark of aging, as well as a promising target for preventing aging and age-related diseases. However, mapping the senescent cells in tissues is extremely challenging, as their low abundance, lack of specific markers, and variability arise from heterogeneity. Hence, methodologies for identifying or predicting the development of senescent cells are necessary for achieving healthy aging. A new wave of bioinformatic methodologies based on mathematics/physics theories have been proposed to be applied to aging biology, which is altering the way we approach our understand of aging. Here, we discuss the dynamical network biomarkers (DNB) theory, which allows for the prediction of state transition in complex systems such as living organisms, as well as usage of Raman spectroscopy that offers a non-invasive and label-free imaging, and provide a perspective on potential applications for the study of aging.

Hypoxia-induced exosomes facilitate lung pre-metastatic niche formation in hepatocellular carcinoma through the miR-4508-RFX1-IL17A-p38 MAPK-NF-κB pathway

Background: Hypoxia plays an important role in the lung metastasis of hepatocellular carcinoma (HCC). However, the process by which hypoxia promotes the formation of a pre-metastatic niche (PMN) and its underlying mechanism remain unclear. Methods: Exosomes derived from normoxic and hypoxic HCC cells were collected to induce fibroblast activation in vitro and PMN formation in vivo. The micro RNA (miR) profiles of the exosomes were sequenced to identify differentially expressed miRNAs. Gain- and loss-of-function analyses were performed to investigate miR-4508 function. Dual-luciferase, western blotting, and real-time reverse transcription-PCR analyses were used to identify the direct targets of miR-4508 and its downstream signaling pathways. To demonstrate the roles of hypoxic tumor-derived exosomes (H-TDEs) and miR-4508 in the lung metastasis of liver cancer, H22 tumor cells were injected through the tail vein of mice. Blood plasma-derived exosomes from patients with HCC who underwent transarterial chemoembolization (TACE) were applied to determine clinical correlations. Results: We demonstrated that H-TDEs activated lung fibroblasts and facilitated PMN formation, thereby promoting lung metastasis in mice. Screening for upregulated exosomal miRNAs revealed that miR-4508 and its target, regulatory factor X1 (RFX1), were involved in H-TDE-induced lung PMN formation. Moreover, miR-4508 was significantly upregulated in plasma exosomes derived from patients with HCC after TACE. We confirmed that the p38 MAPK-NF-κB signaling pathway is involved in RFX1 knockdown-induced fibroblast activation and PMN formation. In addition, IL17A, a downstream target of RFX1, was identified as a link between RFX1 knockdown and p38 MAPK activation in fibroblasts. Conclusion: Hypoxia enhances the release of TDEs enriched with miR-4508, thereby promoting lung PMN formation by targeting the RFX1-IL17A-p38 MAPK-NF-κB pathway. These findings highlight a novel mechanism underlying hypoxia-induced pulmonary metastasis of HCC.

DA-SRN: Omics data analysis based on the sample network optimization for complex diseases

Effective biomarker identification and accurate sample label prediction are still challenging for complex diseases. Patient similarity network (PSN) analysis is a powerful tool in disease omics data analysis. The topology of PSN can reflect the discriminative ability of the corresponding feature space on which the sample network is built. In this study, a novel omics data analysis method based on the sample reference network (DA-SRN) is proposed to identify the potential biomarkers and predict the sample categories. DA-SRN defines the informative features and the sample reference network in optimizing the network structure by genetic algorithm. It labels the samples based on the graph neural network, the reference network and the selected informative features. DA-SRN was compared with nine efficient omics data analysis methods on the genomics, metabolomics and transcriptomics datasets to show its validation. The comparison results showed that it outperformed the other methods in area under receiver operating characteristic curve (AUROC), sensitivity, specificity and area under precision-recall curve (AUPRC) in most cases. Besides, the important metabolites identified by DA-SRN for the type 2 diabetes (T2D) metabolomics data were further examined. The pathway analysis revealed the close relationships between the identified metabolites and the critical metabolic pathways related to the occurrence and development of T2D. The experimental results illustrate that DA-SRN can extract the valuable information from the complex omics data by analyzing the sample relationship, and is promising in biomarker identification and sample discrimination for complex diseases.

Identification of CCDC115 as an adverse prognostic biomarker in liver cancer based on bioinformatics and experimental analyses

Liver hepatocellular carcinoma (LIHC) is a major malignant tumor of the digestive system with a high incidence rate and poor early diagnosis. Coiled-coil domain-containing protein 115 (CCDC115), an accessory component of vacuolar-ATPase with dramatically abnormal expression, is associated with survival outcomes of cancer patients. However, the role of CCDC115 in LIHC remains unclear. In this study, we aimed to determine the functional role of CCDC115 in LIHC by examining CCDC115 expression, and its influence on LIHC prognosis. Through extensive statistical analyses, using LIHC patient databases, we observed that CCDC115 expression significantly increased in tumor tissues of LIHC patients. In addition, CCDC115 expression correlated with the poor prognosis. Additionally, CCDC115 was found to be involved in several cancer-related pathways, specifically the PI3K-Akt pathway. The expression of CCDC115 was positively correlated with human leukocyte antigen molecules as well as with immune checkpoint molecules in LIHC patients. We performed in vitro experiments and confirmed that the expression of CCDC115 significantly affects the proliferation potential, metastasis and sorafenib resistance of liver cancer cells, as well as some key protein expression in PI3K-Akt pathway. These results indicate that CCDC115 could serve as a diagnostic and prognostic biomarker of LIHC, and targeting CCDC115 may provide a potential strategy to enhance the efficacy of liver cancer therapy.

Deciphering gastric inflammation-induced tumorigenesis through multi-omics data and AI methods

Gastric cancer (GC), the fifth most common cancer globally, remains the leading cause of cancer deaths worldwide. Inflammation-induced tumorigenesis is the predominant process in GC development; therefore, systematic research in this area should improve understanding of the biological mechanisms that initiate GC development and promote cancer hallmarks. Here, we summarize biological knowledge regarding gastric inflammation-induced tumorigenesis, and characterize the multi-omics data and systems biology methods for investigating GC development. Of note, we highlight pioneering studies in multi-omics data and state-of-the-art network-based algorithms used for dissecting the features of gastric inflammation-induced tumorigenesis, and we propose translational applications in early GC warning biomarkers and precise treatment strategies. This review offers integrative insights for GC research, with the goal of paving the way to novel paradigms for GC precision oncology and prevention.

Glucose-Responsive Antioxidant Hydrogel Accelerates Diabetic Wound Healing

Diabetic complications can be ameliorated by inhibiting excessive oxidative stress with antioxidants. To enhance therapeutic intervention, it is crucial to develop intelligent scaffolds for efficient delivery of antioxidants to diabetic wounds. This study introduces reversible boronic bonds to create an intelligent antioxidant hydrogel scaffold. This study modifies gelatin methacryloyl (GelMA) with 4-carboxyphenyboronic acid (CPBA) to synthesize a derivative of GelMA (GelMA-CPBA), and then photo cross-links GelMA-CPBA with (-)-epigallocatechin-3-gallate (EGCG) to form GelMA-CPBA/EGCG (GMPE) hydrogel. The GMPE hydrogel responds to changes in glucose levels, and more EGCG is released as glucose level increases due to the dissociation of boronic ester bonds. The GMPE hydrogel shows good biocompatibility and biodegradability, and its mechanical property is similar to that of the skin tissue. Both in vitro and in vivo results demonstrate that the GMPE hydrogel scaffolds effectively eliminate reactive oxygen species (ROS), reduce the inflammation, and promote angiogenesis, thereby improve collagen deposition and tissue remodeling during diabetic wound healing. This strategy offers new insight into glucose-responsive scaffolds, and this responsive antioxidan hydrogel scaffold holds great potential for the treatment of chronic diabetic wounds.

Early warning and diagnosis of liver cancer based on dynamic network biomarker and deep learning

Background

Early detection of complex diseases like hepatocellular carcinoma remains challenging due to their network-driven pathology. Dynamic network biomarkers (DNB) based on monitoring changes in molecular correlations may enable earlier predictions. However, DNB analysis often overlooks disease heterogeneity.

Methods

We integrated DNB analysis with graph convolutional neural networks (GCN) to identify critical transitions during hepatocellular carcinoma development in a mouse model. A DNB-GCN model was constructed using transcriptomic data and gene expression levels as node features.

Results

DNB analysis identified a critical transition point at 7 weeks of age despite histological examinations being unable to detect cancerous changes at that time point. The DNB-GCN model achieved 100% accuracy in classifying healthy and cancerous mice, and was able to accurately predict the health status of newly introduced mice.

Conclusion

The integration of DNB analysis and GCN demonstrates potential for the early detection of complex diseases by capturing network structures and molecular features that conventional biomarker discovery methods overlook. The approach warrants further development and validation.

Tumor type classification and candidate cancer-specific biomarkers discovery via semi-supervised learning

Identifying cancer-related differentially expressed genes provides significant information for diagnosing tumors, predicting prognoses, and effective treatments. Recently, deep learning methods have been used to perform gene differential expression analysis using microarray-based high-throughput gene profiling and have achieved good results. In this study, we proposed a new robust multiple-datasets-based semi-supervised learning model, MSSL, to perform tumor type classification and candidate cancer-specific biomarkers discovery across multiple tumor types and multiple datasets, which addressed the following long-lasting obstacles: (1) the data volume of the existing single dataset is not enough to fully exert the advantages of deep learning; (2) a large number of datasets from different research institutions cannot be effectively used due to inconsistent internal variances and low quality; (3) relatively uncommon cancers have limited effects on deep learning methods. In our article, we applied MSSL to The Cancer Genome Atlas (TCGA) and the Gene Expression Comprehensive Database (GEO) pan-cancer normalized-level3 RNA-seq data and got 97.6% final classification accuracy, which had a significant performance leap compared with previous approaches. Finally, we got the ranking of the importance of the corresponding genes for each cancer type based on classification results and validated that the top genes selected in this way were biologically meaningful for corresponding tumors and some of them had been used as biomarkers, which showed the efficacy of our method.

Personalised virtual brain models in epilepsy

Individuals with drug-resistant focal epilepsy are candidates for surgical treatment as a curative option. Before surgery can take place, the patient must have a presurgical evaluation to establish whether and how surgical treatment might stop their seizures without causing neurological deficits. Virtual brains are a new digital modelling technology that map the brain network of a person with epilepsy, using data derived from MRI. This technique produces a computer simulation of seizures and brain imaging signals, such as those that would be recorded with intracranial EEG. When combined with machine learning, virtual brains can be used to estimate the extent and organisation of the epileptogenic zone (ie, the brain regions related to seizure generation and the spatiotemporal dynamics during seizure onset). Virtual brains could, in the future, be used for clinical decision making, to improve precision in localisation of seizure activity, and for surgical planning, but at the moment these models have some limitations, such as low spatial resolution. As evidence accumulates in support of the predictive power of personalised virtual brain models, and as methods are tested in clinical trials, virtual brains might inform clinical practice in the near future.

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 biomarker to determine the critical point of breast cancer stage progression

BACKGROUND

The discovery of early warning signs and biomarkers in patients with early breast cancer is crucial for the prevention and treatment of breast cancer. Dynamic Network Biomarker (DNB) is an approach based on nonlinear dynamics theory, which we exploited to identify a set of DNB members and their key genes as early warning signals during breast cancer staging progression.

METHODS

First, based on the gene expression profile of breast cancer in the TCGA database, the DNB algorithm was used to calculate the composite index (CI) of each gene cluster in the process of breast cancer anatomical staging. Then we calculated gene modules associated with the clinical phenotype stage based on weighted gene co-expression network analysis (WGCNA), combined with DNB membership to identify key genes in the network.

RESULTS

We identified a set of gene clusters with the highest CI in Stage II as DNBs, whose roles in related pathways indicate the emergence of a tipping point and impact on breast cancer development. In addition, analysis of the key gene GPRIN1 showed that high expression of GPRIN1 predicts poor prognosis, and related immune analysis showed that GPRIN1 is involved in the development of breast cancer through immune aspects.

CONCLUSION

The discovery of DNBs and the key gene GPRIN1 can provide potential biomarkers and therapeutic targets for breast cancer.

Delineating epileptogenic networks using brain imaging data and personalized modeling in drug-resistant epilepsy

Precise estimates of epileptogenic zone networks (EZNs) are crucial for planning intervention strategies to treat drug-resistant focal epilepsy. Here, we present the virtual epileptic patient (VEP), a workflow that uses personalized brain models and machine learning methods to estimate EZNs and to aid surgical strategies. The structural scaffold of the patient-specific whole-brain network model is constructed from anatomical T1 and diffusion-weighted magnetic resonance imaging. Each network node is equipped with a mathematical dynamical model to simulate seizure activity. Bayesian inference methods sample and optimize key parameters of the personalized model using functional stereoelectroencephalography recordings of patients' seizures. These key parameters together with their personalized model determine a given patient's EZN. Personalized models were further used to predict the outcome of surgical intervention using virtual surgeries. We evaluated the VEP workflow retrospectively using 53 patients with drug-resistant focal epilepsy. VEPs reproduced the clinically defined EZNs with a precision of 0.6, where the physical distance between epileptogenic regions identified by VEP and the clinically defined EZNs was small. Compared with the resected brain regions of 25 patients who underwent surgery, VEP showed lower false discovery rates in seizure-free patients (mean, 0.028) than in non-seizure-free patients (mean, 0.407). VEP is now being evaluated in an ongoing clinical trial (EPINOV) with an expected 356 prospective patients with epilepsy.

Low expression of the dynamic network markers FOS/JUN in pre-deteriorated epithelial cells is associated with the progression of colorectal adenoma to carcinoma

BACKGROUND

Deterioration of normal intestinal epithelial cells is crucial for colorectal tumorigenesis. However, the process of epithelial cell deterioration and molecular networks that contribute to this process remain unclear.

METHODS

Single-cell data and clinical information were downloaded from the Gene Expression Omnibus (GEO) database. We used the recently proposed dynamic network biomarker (DNB) method to identify the critical stage of epithelial cell deterioration. Data analysis and visualization were performed using R and Cytoscape software. In addition, Single-Cell rEgulatory Network Inference and Clustering (SCENIC) analysis was used to identify potential transcription factors, and CellChat analysis was conducted to evaluate possible interactions among cell populations. Gene Ontology (GO), Kyoto Encyclopedia of Genes and Genomes (KEGG), and gene set variation analysis (GSVA) analyses were also performed.

RESULTS

The trajectory of epithelial cell deterioration in adenoma to carcinoma progression was delineated, and the subpopulation of pre-deteriorated epithelial cells during colorectal cancer (CRC) initialization was identified at the single-cell level. Additionally, FOS/JUN were identified as biomarkers for pre-deteriorated epithelial cell subpopulations in CRC. Notably, FOS/JUN triggered low expression of P53-regulated downstream pro-apoptotic genes and high expression of anti-apoptotic genes through suppression of P53 expression, which in turn inhibited P53-induced apoptosis. Furthermore, malignant epithelial cells contributed to the progression of pre-deteriorated epithelial cells through the GDF signaling pathway.

CONCLUSIONS

We demonstrated the trajectory of epithelial cell deterioration and used DNB to characterize pre-deteriorated epithelial cells at the single-cell level. The expression of DNB-neighboring genes and cellular communication were triggered by DNB genes, which may be involved in epithelial cell deterioration. The DNB genes FOS/JUN provide new insights into early intervention in CRC.

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.

Dynamic network biomarker factors orchestrate cell-fate determination at tipping points during hESC differentiation

The generation of ectoderm, mesoderm, and endoderm layers is the most critical biological process during the gastrulation of embryo development. Such a differentiation process in human embryonic stem cells (hESCs) is an inherently nonlinear multi-stage dynamical process which contain multiple tipping points playing crucial roles in the cell-fate decision. However, the tipping points of the process are largely unknown, letting alone the understanding of the molecular regulation on these critical events. Here by designing a module-based dynamic network biomarker (M-DNB) model, we quantitatively pinpointed two tipping points of the differentiation of hESCs toward definitive endoderm, which leads to the identification of M-DNB factors (FOS, HSF1, MYCN, TP53, and MYC) of this process. We demonstrate that before the tipping points, M-DNB factors are able to maintain the cell states and orchestrate cell-fate determination during hESC (ES)-to-ME and ME-to-DE differentiation processes, which not only leads to better understanding of endodermal specification of hESCs but also reveals the power of the M-DNB model to identify critical transition points with their key factors in diverse biological processes, including cell differentiation and transdifferentiation dynamics.

Construction of multifunctional hydrogel with metal-polyphenol capsules for infected full-thickness skin wound healing

Damaged skin cannot prevent harmful bacteria from invading tissues, causing infected wounds or even severe tissue damage. In this study, we developed a controlled-release antibacterial composite hydrogel system that can promote wound angiogenesis and inhibit inflammation by sustained releasing Cu-Epigallocatechin-3-gallate (Cu-EGCG) nano-capsules. The prepared SilMA/HAMA/Cu-EGCG hydrogel showed an obvious inhibitory effect on Escherichia coli (E. coli) and Staphylococcus aureus (S. aureus). It could also promote the proliferation and migration of L929 fibroblasts. In vivo full-thickness infected wound healing experiments confirmed the angiogenesis and inflammation regulating effect. Accelerate collagen deposition and wound healing speed were also observed in the SilMA/HAMA/Cu-EGCG hydrogel treated group. The findings of this study show the great potential of this controlled-release antibacterial composite hydrogel in the application of chronic wound healing.

ARHGEF37 overexpression promotes extravasation and metastasis of hepatocellular carcinoma via directly activating Cdc42

Background

The extravasation capability of hepatocellular carcinoma (HCC) cells plays a vital role in distant metastasis. However, the underlying mechanism of extravasation in HCC lung metastasis remains largely unclear.

Methods

The expression of ARHGEF37 in human HCC specimens and HCC cell lines was examined by quantitative RT-PCR, western blot, and immunohistochemistry (IHC) analyses. The biological roles and mechanisms of ARHGEF37/Cdc42 in promoting lung metastasis were investigated in vitro and in vivo using cell lines, patient samples, xenograft models.

Results

In the current study, we found that Rho guanine nucleotide exchange factor 37 (ARHGEF37) was upregulated in human HCC samples and was associated with tumor invasiveness, pulmonary metastasis and poor prognosis. Overexpressing ARHGEF37 significantly enhanced the extravasation and metastatic capability of HCC cells via facilitating tumor cell adhesion to endothelial cells and trans-endothelial migration. Mechanistically, ARHGEF37 directly interacted with and activated Cdc42 to promote the invadopodia formation in HCC cells, which consequently disrupted the interaction between endothelial cells and pericytes. Importantly, treatment with ZCL278, a specific inhibitor of Cdc42, dramatically inhibited the attachment of ARHGEF37-overexpressing HCC cells to endothelial cells, and the adherence and extravasation in the lung alveoli, resulting in suppression of lung metastasis in mice.

Conclusion

Our findings provide a new insight into the underlying mechanisms on the ARHGEF37 overexpression-mediated extravasation and pulmonary metastasis of HCC cells, and provided a potential therapeutic target for the prevention and treatment of HCC pulmonary metastasis.

Supplementary Information

The online version contains supplementary material available at 10.1186/s13046-022-02441-y.

Modular networks and genomic variation during progression from stable angina pectoris through ischemic cardiomyopathy to chronic heart failure

Background

Analyzing disease–disease relationships plays an important role for understanding etiology, disease classification, and drug repositioning. However, as cardiovascular diseases with causative links, the molecular relationship among stable angina pectoris (SAP), ischemic cardiomyopathy (ICM) and chronic heart failure (CHF) is not clear.

Methods

In this study, by integrating the multi-database data, we constructed paired disease progression modules (PDPMs) to identified relationship among SAP, ICM and CHF based on module reconstruction pairs (MRPs) of K-value calculation (a Euclidean distance optimization by integrating module topology parameters and their weights) methods. Finally, enrichment analysis, literature validation and structural variation (SV) were performed to verify the relationship between the three diseases in PDPMs.

Results

Total 16 PDPMs were found with K > 0.3777 among SAP, ICM and CHF, in which 6 pairs in SAP–ICM, 5 pairs for both ICM–CHF and SAP–CHF. SAP–ICM was the most closely related by having the smallest average K-value (K = 0.3899) while the maximum is SAP–CHF (K = 0.4006). According to the function of the validation gene, inflammatory response were through each stage of SAP–ICM–CHF, while SAP–ICM was uniquely involved in fibrosis, and genes were related in affecting the upstream of PI3K–Akt signaling pathway. 4 of the 11 genes (FLT1, KDR, ANGPT2 and PGF) in SAP–ICM–CHF related to angiogenesis in HIF-1 signaling pathway. Furthermore, we identified 62.96% SVs were protein deletion in SAP–ICM–CHF, and 53.85% SVs were defined as protein replication in SAP–ICM, while ICM–CHF genes were mainly affected by protein deletion.

Conclusion

The PDPMs analysis approach combined with genomic structural variation provides a new avenue for determining target associations contributing to disease progression and reveals that inflammation and angiogenesis may be important links among SAP, ICM and CHF progression.

Application of the Dynamical Network Biomarker Theory to Raman Spectra

The dynamical network biomarker (DNB) theory detects the early warning signals of state transitions utilizing fluctuations in and correlations between variables in complex systems. Although the DNB theory has been applied to gene expression in several diseases, destructive testing by microarrays is a critical issue. Therefore, other biological information obtained by non-destructive testing is desirable; one such piece of information is Raman spectra measured by Raman spectroscopy. Raman spectroscopy is a powerful tool in life sciences and many other fields that enable the label-free non-invasive imaging of live cells and tissues along with detailed molecular fingerprints. Naïve and activated T cells have recently been successfully distinguished from each other using Raman spectroscopy without labeling. In the present study, we applied the DNB theory to Raman spectra of T cell activation as a model case. The dataset consisted of Raman spectra of the T cell activation process observed at 0 (naïve T cells), 2, 6, 12, 24 and 48 h (fully activated T cells). In the DNB analysis, the F-test and hierarchical clustering were used to detect the transition state and identify DNB Raman shifts. We successfully detected the transition state at 6 h and related DNB Raman shifts during the T cell activation process. The present results suggest novel applications of the DNB theory to Raman spectra ranging from fundamental research on cellular mechanisms to clinical examinations.

CSP I-plus modified rEndostatin inhibits hepatocellular carcinoma metastasis via down-regulation of VEGFA and integrinβ1

BACKGROUND

In our previous study, N end of the Circumsporozoite protein (CSP I-plus) modified recombinant human Endostatin (rEndostatin, endostar) (rES-CSP) was constructed, which had antiangiogenic capability and bound to hepatocellular carcinoma in vivo and in vitro. In this study, the inhibition of rES-CSP on hepatocellular carcinoma metastasis was verified in vivo and in vitro, and its possible mechanism was explored.

METHODS

Firstly, the impact of rES-CSP on the migration, adhesion of hepatoma cell HCCLM3 was identified by wound healing, transwell, and on metastasis of orthotopic xenograft model was identified in nude mouse. Then the expression of metastasis-associated molecules (MMP2, E-cadherin, integrinβ1) and angiogenesis-related factors (VEGFA) in vitro and in vivo were detected by real-time PCR, western blotting, immunohistochemistry.

RESULTS

Finally, we found that rES-CSP could inhibit the migration and invasion of HCCLM3, and decrease tumor metastasis and growth in nude mouse orthotopic xenograft models. The tumor inhibiting rates of rES-CSP and Endostar were 42.46 ± 5.39% and 11.1 ± 1.88%. The lung metastasis rates of the control, Endostar and rES-CSP were 71, 50, and 42.8%, respectively. Compared with Endostar, rES-CSP significantly down-regulated the expression of VEGFA and integrinβ1. Heparin, a competitive inhibitor of CSP I-plus, which can be bind to the highly-sulfated heparan sulfate proteoglycans (HSPGs) over-expressed in liver and hepatocellular carcinoma, alleviated the down-regulation of VEGFA and integrinβ1.

CONCLUSIONS

These indicate that rES-CSP may play a role in inhibiting tumor growth and metastasis by down-regulating the angiogenic factor VEGF and the metastasis-related molecules or by interfering with HSPGs-mediated tumor metastasis.

Medical care of rare and undiagnosed diseases: Prospects and challenges

Rare and undiagnosed diseases tend to be diverse, misdiagnosed, and difficult to diagnose. In some cases, the disease is progressive and life-threatening. Yet, to date, an estimated 95% of rare diseases have no approved therapy. Therefore, rare and undiagnosed diseases are considered the ultimate challenges for understanding human diseases. Here, we review the research progress, research frontiers, and important scientific issues related to rare and undiagnosed diseases. We mainly focus on five topics: (1) the identification and functional analysis of disease-causing genes; (2) the construction of cells, organoids, and animal models for mechanism validation; (3) subtyping and diagnosis; (4) treatment and drug screening based on causative genes and mutations; and (5) new technologies and methods for studying rare and undiagnosed diseases. In this review, we briefly update and discuss the pathogenic mechanisms and precision medicine for rare and undiagnosed diseases.

Impact of preanalytical freezing delay time on the stability of metabolites in oral squamous cell carcinoma tissue samples

INTRODUCTION

Metabolite stability is critical for tissue metabolomics. However, changes in metabolites in tissues over time from the operating room to the laboratory remain underexplored.

OBJECTIVES

In this study, we evaluated the effect of postoperative freezing delay time on the stability of metabolites in normal and oral squamous cell carcinoma (OSCC) tissues.

METHODS

Tumor and paired normal tissues from five OSCC patients were collected after surgical resection, and samples was sequentially quenched in liquid nitrogen at 30, 40, 50, 60, 70, 80, 90 and 120 min (80 samples). Untargeted metabolic analysis by liquid chromatography-mass spectrometry/mass spectrometry in positive and negative ion modes was used to identify metabolic changes associated with delayed freezing time. The trends of metabolite changes at 30-120 and 30-60 min of delayed freezing were analyzed.

RESULTS

190 metabolites in 36 chemical classes were detected. After delayed freezing for 120 min, approximately 20% of the metabolites changed significantly in normal and tumor tissues, and differences in the metabolites were found in normal and tumor tissues. After a delay of 60 min, 29 metabolites had changed significantly in normal tissues, and 84 metabolites had changed significantly in tumor tissues. In addition, we constructed three tissue freezing schemes based on the observed variation trends in the metabolites.

CONCLUSION

Delayed freezing of tissue samples has a certain impact on the stability of metabolites. For metabolites with significant changes, we suggest that the freezing time of tissues be reasonably selected according to the freezing schemes and the actual clinical situation.

APEX1 promotes the oncogenicity of hepatocellular carcinoma via regulation of MAP2K6

Objective: Apurinic/apyrimidinic endonuclease 1 (APEX1), a key enzyme responsible for DNA base excision repair, has been linked to development and progression of cancers. In this work, we aimed to explore the role of APEX1 in hepatocellular carcinoma (HCC) and elucidate its molecular mechanism.

Methods: The expression of APEX1 in HCC tissues and matched adjacent normal tissues (n = 80 cases) was evaluated by immunohistochemistry. Web-based tools UALCAN and the Kaplan-Meier plotter were used to analyze the Cancer Genome Atlas database to compare expression of APEX1 mRNA to 5-year overall survival. APEX1 was stably silenced in two HCC cell lines, Hep 3B and Bel-7402, with shRNA technology. An in vivo tumorigenesis model was established by subcutaneously injecting sh-APEX1-transfected Bel-7402 cells into mice, and tumor growth was determined. We performed high-throughput transcriptome sequencing in sh-APEX1-treated HCC cells to identify the key KEGG signaling pathways induced by silencing of APEX1.

Results: APEX1 was significantly upregulated and predicted poor clinical overall survival in HCC patients. Silencing APEX1 inhibited the proliferation of HCC cells in vivo and in vitro, and it repressed invasion and migration and increased apoptosis and the percentage of cells in G1. Differentially expressed genes upon APEX1 silencing included genes involved in TNF signaling. A positive correlation between the expression of APEX1 and MAP2K6 was noted, and overexpressing MAP2K6 overcame cancer-related phenotypes associated with APEX1 silencing.

Conclusion: APEX1 enhances the malignant properties of HCC via MAP2K6. APEX1 may represent a valuable prognostic biomarker and therapeutic target in HCC.

scMTD: a statistical multidimensional imputation method for single-cell RNA-seq data leveraging transcriptome dynamic information

Background

Single-cell RNA sequencing (scRNA-seq) provides a powerful tool to capture transcriptomes at single-cell resolution. However, dropout events distort the gene expression levels and underlying biological signals, misleading the downstream analysis of scRNA-seq data.

Results

We develop a statistical model-based multidimensional imputation algorithm, scMTD, that identifies local cell neighbors and specific gene co-expression networks based on the pseudo-time of cells, leveraging information on cell-level, gene-level, and transcriptome dynamic to recover scRNA-seq data. Compared with the state-of-the-art imputation methods through several real-data-based analytical experiments, scMTD effectively recovers biological signals of transcriptomes and consistently outperforms the other algorithms in improving FISH validation, trajectory inference, differential expression analysis, clustering analysis, and identification of cell types.

Conclusions

scMTD maintains the gene expression characteristics, enhances the clustering of cell subpopulations, assists the study of gene expression dynamics, contributes to the discovery of rare cell types, and applies to both UMI-based and non-UMI-based data. Overall, scMTD’s reliability, applicability, and scalability make it a promising imputation approach for scRNA-seq data.

Supplementary Information

The online version contains supplementary material available at 10.1186/s13578-022-00886-4.

Buffering variability in cell regulation motifs close to criticality

Bistable biological regulatory systems need to cope with stochastic noise to fine tune their function close to bifurcation points. Here, we study stability properties of this regime in generic systems to demonstrate that cooperative interactions buffer system variability, hampering noise-induced regime shifts. Our analysis also shows that, in the considered cooperativity range, impending regime shifts can be generically detected by statistical early warning signals from distributional data. Our generic framework, based on minimal models, can be used to extract robustness and variability properties of more complex models and empirical data close to criticality.

The fundamental benefits of multiplexity in ecological networks

A tipping point presents perhaps the single most significant threat to an ecological system as it can lead to abrupt species extinction on a massive scale. Climate changes leading to the species decay parameter drifts can drive various ecological systems towards a tipping point. We investigate the tipping-point dynamics in multi-layer ecological networks supported by mutualism. We unveil a natural mechanism by which the occurrence of tipping points can be delayed by multiplexity that broadly describes the diversity of the species abundances, the complexity of the interspecific relationships, and the topology of linkages in ecological networks. For a double-layer system of pollinators and plants, coupling between the network layers occurs when there is dispersal of pollinator species. Multiplexity emerges as the dispersing species establish their presence in the destination layer and have a simultaneous presence in both. We demonstrate that the new mutualistic links induced by the dispersing species with the residence species have fundamental benefits to the well-being of the ecosystem in delaying the tipping point and facilitating species recovery. Articulating and implementing control mechanisms to induce multiplexity can thus help sustain certain types of ecosystems that are in danger of extinction as the result of environmental changes.

Single-cell entropy network detects the activity of immune cells based on ribosomal protein genes

We developed a new computational method, Single-Cell Entropy Network (SCEN) to analyze single-cell RNA-seq data, which used the information of gene-gene associations to discover new heterogeneity of immune cells as well as identify existing cell types. Based on SCEN, we defined association-entropy (AE) for each cell and each gene through single-cell gene co-expression networks to measure the strength of association between each gene and all other genes at a single-cell resolution. Analyses of public datasets indicated that the AE of ribosomal protein genes (RP genes) varied greatly even in the same cell type of immune cells and the average AE of RP genes of immune cells in each person was significantly associated with the healthy/disease state of this person. Based on existing research and theory, we inferred that the AE of RP genes represented the heterogeneity of ribosomes and reflected the activity of immune cells. We believe SCEN can provide more biological insights into the heterogeneity and diversity of immune cells, especially the change of immune cells in the diseases.

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.

Energy landscape decomposition for cell differentiation with proliferation effect

Complex interactions between genes determine the development and differentiation of cells. We establish a landscape theory for cell differentiation with proliferation effect, in which the developmental process is modeled as a stochastic dynamical system with a birth-death term. We find that two different energy landscapes, denoted U and V, collectively contribute to the establishment of non-equilibrium steady differentiation. The potential U is known as the energy landscape leading to the steady distribution, whose metastable states stand for cell types, while V indicates the differentiation direction from pluripotent to differentiated cells. This interpretation of cell differentiation is different from the previous landscape theory without the proliferation effect. We propose feasible numerical methods and a mean-field approximation for constructing landscapes U and V. Successful applications to typical biological models demonstrate the energy landscape decomposition's validity and reveal biological insights into the considered processes.

Detecting critical transition signals from single-cell transcriptomes to infer lineage-determining transcription factors

Analyzing single-cell transcriptomes promises to decipher the plasticity, heterogeneity, and rapid switches in developmental cellular state transitions. Such analyses require the identification of gene markers for semi-stable transition states. However, there are nontrivial challenges such as unexplainable stochasticity, variable population sizes, and alternative trajectory constructions. By advancing current tipping-point theory-based models with feature selection, network decomposition, accurate estimation of correlations, and optimization, we developed BioTIP to overcome these challenges. BioTIP identifies a small group of genes, called critical transition signal (CTS), to characterize regulated stochasticity during semi-stable transitions. Although methods rooted in different theories converged at the same transition events in two benchmark datasets, BioTIP is unique in inferring lineage-determining transcription factors governing critical transition. Applying BioTIP to mouse gastrulation data, we identify multiple CTSs from one dataset and validated their significance in another independent dataset. We detect the established regulator Etv2 whose expression change drives the haemato-endothelial bifurcation, and its targets together in CTS across three datasets. After comparing to three current methods using six datasets, we show that BioTIP is accurate, user-friendly, independent of pseudo-temporal trajectory, and captures significantly interconnected and reproducible CTSs. We expect BioTIP to provide great insight into dynamic regulations of lineage-determining factors.

Determination of WWOX Function in Modulating Cellular Pathways Activated by AP-2α and AP-2γ Transcription Factors in Bladder Cancer

Following the invention of high-throughput sequencing, cancer research focused on investigating disease-related alterations, often inadvertently omitting tumor heterogeneity. This research was intended to limit the impact of heterogeneity on conclusions related to WWOX/AP-2α/AP-2γ in bladder cancer which differently influenced carcinogenesis. The study examined the signaling pathways regulated by WWOX-dependent AP-2 targets in cell lines as biological replicates using high-throughput sequencing. RT-112, HT-1376 and CAL-29 cell lines were subjected to two stable lentiviral transductions. Following CAGE-seq and differential expression analysis, the most important genes were identified and functionally annotated. Western blot was performed to validate the selected observations. The role of genes in biological processes was assessed and networks were visualized. Ultimately, principal component analysis was performed. The studied genes were found to be implicated in MAPK, Wnt, Ras, PI3K-Akt or Rap1 signaling. Data from pathways were collected, explaining the differences/similarities between phenotypes. FGFR3, STAT6, EFNA1, GSK3B, PIK3CB and SOS1 were successfully validated at the protein level. Afterwards, a definitive network was built using 173 genes. Principal component analysis revealed that the various expression of these genes explains the phenotypes. In conclusion, the current study certified that the signaling pathways regulated by WWOX and AP-2α have more in common than that regulated by AP-2γ. This is because WWOX acts as an EMT inhibitor, AP-2γ as an EMT enhancer while AP-2α as a MET inducer. Therefore, the relevance of AP-2γ in targeted therapy is now more evident. Some of the differently regulated genes can find application in bladder cancer treatment.