Dynamic edge-based biomarker non-invasively predicts hepatocellular carcinoma with hepatitis B virus infection for individual patients based on blood testing

Emerging biomolecules for practical theranostics of liver hepatocellular carcinoma

Most cases of hepatocellular carcinoma (HCC) are able to be diagnosed through regular surveillance in an identifiable patient population with chronic hepatitis B or cirrhosis. Nevertheless, 50% of global cases might present incidentally owing to symptomatic advanced-stage HCC after worsening of liver dysfunction. A systematic search based on PUBMED was performed to identify relevant outcomes, covering newer surveillance modalities including secretory proteins, DNA methylation, miRNAs, and genome sequencing analysis which proposed molecular expression signatures as ideal tools in the early-stage HCC detection. In the face of low accuracy without harmonization on the analytical approaches and data interpretation for liquid biopsy, a more accurate incidence of HCC will be unveiled by using deep machine learning system and multiplex immunohistochemistry analysis. A combination of molecular-secretory biomarkers, high-definition imaging and bedside clinical indexes in a surveillance setting offers a comprehensive range of HCC potential indicators. In addition, the sequential use of numerous lines of systemic anti-HCC therapies will simultaneously benefit more patients in survival. This review provides an overview on the most recent developments in HCC theranostic platform.

Mechanisms of Ganweikang Tablets against Chronic Hepatitis B: A Comprehensive Study of Network Analysis, Molecular Docking, and Chemical Profiling

The hepatitis B virus (HBV) is one of the major viral infection problems worldwide in public health. The exclusive proprietary Chinese medicine Ganweikang (GWK) tablet has been marketed for years in the treatment of chronic hepatitis B (CHB). However, the pharmacodynamic material basis and underlying mechanism of GWK are not completely clear. This study is aimed at investigating the pharmacological mechanism of the GWK tablet in the treatment of CHB. The chemical ingredient information was obtained from the Traditional Chinese Medicine Database and Analysis Platform (TCMSP), Traditional Chinese Medicines Integrated Database (TCMID), and Shanghai Institute of Organic Chemistry of CAS. Ingredients and disease-related targets were defined by a combination of differentially expressed genes from CHB transcriptome data and open-source databases. Target-pathway-target (TPT) network analysis, molecular docking, and chemical composition analysis were adopted to further verify the key targets and corresponding active ingredients of GWK. Eight herbs of GWK were correlated to 330 compounds with positive oral bioavailability, and 199 correlated targets were identified. The TPT network was constructed based on the 146 enriched targets by KEGG pathway analysis, significantly associated with 95 pathways. Twenty-five nonvolatile components and 25 volatile components in GWK were identified in UPLC-QTOF/MS and GC-MS chromatograms. The key active ingredients of GWK include ferulic acid, oleanolic acid, ursolic acid, tormentic acid, 11-deoxyglycyrrhetic acid, dibenzoyl methane, anisaldehyde, wogonin, protocatechuic acid, psoralen, caffeate, dimethylcaffeic acid, vanillin, β-amyrenyl acetate, formonentin, aristololactam IIIa, and 7-methoxy-2-methyl isoflavone, associated with targets CA2, NFKB1, RELA, AKT1, JUN, CA1, CA6, IKBKG, FOS, EP300, CREB1, STAT1, MMP9, CDK2, ABCB1, and ABCG2.

The Variation of Transcriptomic Perturbations is Associated with the Development and Progression of Various Diseases

Background

Although transcriptomic data have been widely applied to explore various diseases, few studies have investigated the association between transcriptomic perturbations and disease development in a wide variety of diseases.

Methods

Based on a previously developed algorithm for quantifying intratumor heterogeneity at the transcriptomic level, we defined the variation of transcriptomic perturbations (VTP) of a disease relative to the health status. Based on publicly available transcriptome datasets, we compared VTP values between the disease and health status and analyzed correlations between VTP values and disease progression or severity in various diseases, including neurological disorders, infectious diseases, cardiovascular diseases, respiratory diseases, liver diseases, kidney diseases, digestive diseases, and endocrine diseases. We also identified the genes and pathways whose expression perturbations correlated positively with VTP across diverse diseases.

Results

VTP values were upregulated in various diseases relative to their normal controls. VTP values were significantly greater in define than in possible or probable Alzheimer's disease. VTP values were significantly larger in intensive care unit (ICU) COVID-19 patients than in non-ICU patients, and in COVID-19 patients requiring mechanical ventilatory support (MVS) than in those not requiring MVS. VTP correlated positively with viral loads in acquired immune deficiency syndrome (AIDS) patients. Moreover, the AIDS patients treated with abacavir or zidovudine had lower VTP values than those without such therapies. In pulmonary tuberculosis (TB) patients, VTP values followed the pattern: active TB > latent TB > normal controls. VTP values were greater in clinically apparent than in presymptomatic malaria. VTP correlated negatively with the cardiac index of left ventricular ejection fraction (LVEF). In chronic obstructive pulmonary disease (COPD), VTP showed a negative correlation with forced expiratory volume in the first second (FEV1). VTP values increased with H. pylori infection and were upregulated in atrophic gastritis caused by H. pylori infection. The genes and pathways whose expression perturbations correlated positively with VTP scores across diseases were mainly involved in the regulation of immune, metabolic, and cellular activities.

Conclusions

VTP is upregulated in the disease versus health status, and its upregulation is associated with disease progression and severity in various diseases. Thus, VTP has potential clinical implications for disease diagnosis and prognosis.

Dynamic network biomarker analysis and system pharmacology methods to explore the therapeutic effects and targets of Xiaoyaosan against liver cirrhosis

ETHNOPHARMACOLOGICAL RELEVANCE

Xiaoyaosan is a traditional Chinese herbal formula that has long been used to treat liver cirrhosis, liver failure, and hepatocarcinoma (HCC). However, little is known about its mechanism of action and targets in treating chronic liver disease.

AIM OF THE STUDY

This study aimed to detect the critical transition of HCC progression and to explore the regulatory mechanism and targets of Xiaoyaosan treating liver cirrhosis (cirrhosis) using integrative medicinal research involving system biology and pharmacology.

MATERIALS AND METHODS

We recruited chronic liver disease participants to obtain gene expression data and applied the dynamic network biomarker (DNB) method to identify molecular markers and the critical transition. We combined network pharmacology and DNB analysis to locate the potential DNBs (targets). Then we validated the DNBs in the liver cirrhosis rat models using Xiaoyaosan treatment. The expression of genes encoding the four DNBs, including Cebpa, Csf1, Egfr, and Il7r, were further validated in rat liver tissue using Western blot analysis.

RESULTS

We found EGFR, CEBPA, Csf1, Ccnb1, Rrmm2, C3, Il7r, Ccna2, and Peg10 overlap in the DNB list and Xiaoyaosan-Target-Disease (XTD) network constructed using network pharmacology databases. We investigated the diagnostic ability of each member in the DNB cluster and found EGFR, CEBPA, CSF1, and IL7R had high diagnostic abilities with AUC >0.7 and P-value < 0.05. We validated these findings in rats and found that liver function improved significantly and fibrotic changes were relieved in the Xiaoyaosan treatment group. The expression levels of CSF1 and IL7R in the Xiaoyaosan group were significantly lower than those in the cirrhosis model group. In contrast, CEBPA expression in the Xiaoyaosan group was significantly higher than that in the cirrhosis model group. The expression of EGFR in the Xiaoyaosan group was slightly decreased than in the model group but not significantly.

CONCLUSION

Using the DNB method and network pharmacology approach, this study revealed that CEBPA, IL7R, EGFR, and CSF1 expression was remarkably altered in chronic liver disease and thus, may play an important role in driving the progression of cirrhosis. Therefore, CEBPA, IL7R, EGFR, and CSF1 may be important targets of Xiaoyaosan in treating cirrhosis and can be considered for developing novel therapeutics.

Interpreting Functional Impact of Genetic Variations by Network QTL for Genotype–Phenotype Association Study

An enormous challenge in the post-genome era is to annotate and resolve the consequences of genetic variation on diverse phenotypes. The genome-wide association study (GWAS) is a well-known method to identify potential genetic loci for complex traits from huge genetic variations, following which it is crucial to identify expression quantitative trait loci (eQTL). However, the conventional eQTL methods usually disregard the systematical role of single-nucleotide polymorphisms (SNPs) or genes, thereby overlooking many network-associated phenotypic determinates. Such a problem motivates us to recognize the network-based quantitative trait loci (QTL), i.e., network QTL (nQTL), which is to detect the cascade association as genotype → network → phenotype rather than conventional genotype → expression → phenotype in eQTL. Specifically, we develop the nQTL framework on the theory and approach of single-sample networks, which can identify not only network traits (e.g., the gene subnetwork associated with genotype) for analyzing complex biological processes but also network signatures (e.g., the interactive gene biomarker candidates screened from network traits) for characterizing targeted phenotype and corresponding subtypes. Our results show that the nQTL framework can efficiently capture associations between SNPs and network traits (i.e., edge traits) in various simulated data scenarios, compared with traditional eQTL methods. Furthermore, we have carried out nQTL analysis on diverse biological and biomedical datasets. Our analysis is effective in detecting network traits for various biological problems and can discover many network signatures for discriminating phenotypes, which can help interpret the influence of nQTL on disease subtyping, disease prognosis, drug response, and pathogen factor association. Particularly, in contrast to the conventional approaches, the nQTL framework could also identify many network traits from human bulk expression data, validated by matched single-cell RNA-seq data in an independent or unsupervised manner. All these results strongly support that nQTL and its detection framework can simultaneously explore the global genotype–network–phenotype associations and the underlying network traits or network signatures with functional impact and importance.

OUP accepted manuscript

Recent advances of single-cell RNA sequencing (scRNA-seq) technologies have led to extensive study of cellular heterogeneity and cell-to-cell variation. However, the high frequency of dropout events and noise in scRNA-seq data confounds the accuracy of the downstream analysis, i.e. clustering analysis, whose accuracy depends heavily on the selected feature genes. Here, by deriving an entropy decomposition formula, we propose a feature selection method, i.e. an intrinsic entropy (IE) model, to identify the informative genes for accurately clustering analysis. Specifically, by eliminating the ‘noisy’ fluctuation or extrinsic entropy (EE), we extract the IE of each gene from the total entropy (TE), i.e. TE = IE + EE. We show that the IE of each gene actually reflects the regulatory fluctuation of this gene in a cellular process, and thus high-IE genes provide rich information on cell type or state analysis. To validate the performance of the high-IE genes, we conduct computational analysis on both simulated datasets and real single-cell datasets by comparing with other representative methods. The results show that our IE model is not only broadly applicable and robust for different clustering and classification methods, but also sensitive for novel cell types. Our results also demonstrate that the intrinsic entropy/fluctuation of a gene serves as information rather than noise in contrast to its total entropy/fluctuation.

The prominence of potential biomarkers in the diagnosis and management of hepatocellular carcinoma: Current scenario and future anticipation

Hepatocellular carcinoma (HCC) is one of the most aggressive and truculent types of cancer. Early detection of HCC is a massive concern that can boost the overall survival rates of HCC patients. As a result, there is a continual quest for advancements in screening, diagnosis, and treatment strategies to enhance the prognosis at its early stages. However, the confluence of inflammation and cirrhosis hampers the early detection of HCC. The analysis of different types of biomarkers such as tissue biomarkers, serum biomarkers, protein biomarkers, autoantibody markers, and improved imaging techniques has played a vital role in ameliorating HCC monitoring responses. Therefore biomarkers that can identify HCC early with a high degree of sensitivity and specificity might be prodigiously serviceable in the diagnosis and treatment of this notorious disorder. This study offers an overview of the contemporary understanding of several types of biomarkers implicated in hepatocarcinogenesis and their applications in monitoring, diagnosis, and prognosis presage. In additament, we address the role of image techniques associated with HCC diagnosis.

RDAClone: Deciphering Tumor Heterozygosity through Single-Cell Genomics Data Analysis with Robust Deep Autoencoder

Rapid advances in single-cell genomics sequencing (SCGS) have allowed researchers to characterize tumor heterozygosity with unprecedented resolution and reveal the phylogenetic relationships between tumor cells or clones. However, high sequencing error rates of current SCGS data, i.e., false positives, false negatives, and missing bases, severely limit its application. Here, we present a deep learning framework, RDAClone, to recover genotype matrices from noisy data with an extended robust deep autoencoder, cluster cells into subclones by the Louvain-Jaccard method, and further infer evolutionary relationships between subclones by the minimum spanning tree. Studies on both simulated and real datasets demonstrate its robustness and superiority in data denoising, cell clustering, and evolutionary tree reconstruction, particularly for large datasets.

Identification of ARGLU1 as a potential therapeutic target for gastric cancer based on genome-wide functional screening data

BACKGROUND

Due to the molecular mechanism complexity and heterogeneity of gastric cancer (GC), mechanistically interpretable biomarkers were required for predicting prognosis and discovering therapeutic targets for GC patients.

METHODS

Based on a total of 824 GC-specific fitness genes from the Project Score database, LASSOCox regression was performed in TCGA-STAD cohort to construct a GC Prognostic (GCP) model which was then evaluated on 7 independent GC datasets. Targets prioritization was performed in GC organoids. ARGLU1 was selected to further explore the biological function and molecular mechanism. We evaluated the potential of ARGLU1 serving as a promising therapeutic target for GC using patients derived xenograft (PDX) model.

FINDINGS

The 9-gene GCP model showed a statistically significant prognostic performance for GC patients in 7 validation cohorts. Perturbation of SSX4, DDX24, ARGLU1 and TTF2 inhibited GC organoids tumor growth. The results of tissue microarray indicated lower expression of ARGLU1 was correlated with advanced TNM stage and worse overall survival. Over-expression ARGLU1 significantly inhibited GC cells viability in vitro and in vivo. ARGLU1 could enhance the transcriptional level of mismatch repair genes including MLH3, MSH2, MSH3 and MSH6 by potentiating the recruitment of SP1 and YY1 on their promoters. Moreover, inducing ARGLU1 by LNP-formulated saRNA significantly inhibited tumor growth in PDX model.

INTERPRETATION

Based on genome-wide functional screening data, we constructed a 9-gene GCP model with satisfactory predictive accuracy and mechanistic interpretability. Out of nine prognostic genes, ARGLU1 was verified to be a potential therapeutic target for GC.

FUNDING

National Natural Science Foundation of China.

Genomics of Human Fibrotic Diseases: Disordered Wound Healing Response

Fibrotic disease, which is implicated in almost half of all deaths worldwide, is the result of an uncontrolled wound healing response to injury in which tissue is replaced by deposition of excess extracellular matrix, leading to fibrosis and loss of organ function. A plethora of genome-wide association studies, microarrays, exome sequencing studies, DNA methylation arrays, next-generation sequencing, and profiling of noncoding RNAs have been performed in patient-derived fibrotic tissue, with the shared goal of utilizing genomics to identify the transcriptional networks and biological pathways underlying the development of fibrotic diseases. In this review, we discuss fibrosing disorders of the skin, liver, kidney, lung, and heart, systematically (1) characterizing the initial acute injury that drives unresolved inflammation, (2) identifying genomic studies that have defined the pathologic gene changes leading to excess matrix deposition and fibrogenesis, and (3) summarizing therapies targeting pro-fibrotic genes and networks identified in the genomic studies. Ultimately, successful bench-to-bedside translation of observations from genomic studies will result in the development of novel anti-fibrotic therapeutics that improve functional quality of life for patients and decrease mortality from fibrotic diseases.

New Blood Biomarkers for the Diagnosis of AFP-Negative Hepatocellular Carcinoma

An early diagnosis of hepatocellular carcinoma (HCC) followed by effective treatment is currently critical for improving the prognosis and reducing the associated economic burden. Alpha-fetoprotein (AFP) is the most widely used biomarker for HCC diagnosis. Based on elevated serum AFP levels as well as typical imaging features, AFP-positive HCC (APHC) can be easily diagnosed, but AFP-negative HCC (ANHC) is not easily detected due to lack of ideal biomarkers and thus mainly reliance on imaging. Imaging for the diagnosis of ANHC is probably insufficient in sensitivity and/or specificity because most ANHC tumors are small and early-stage HCC, and it is involved in sophisticated techniques and high costs. Moreover, ANHC accounts for nearly half of HCC and exhibits a better prognosis compared with APHC. Therefore, the diagnosis of ANHC in clinical practice has been a critical issue for the early treatment and prognosis improvement of HCC. In recent years, tremendous efforts have been made to discover new biomarkers complementary to AFP for HCC diagnosis. In this review, we systematically review and discuss the recent advances of blood biomarkers for HCC diagnosis, including DNA biomarkers, RNA biomarkers, protein biomarkers, and conventional laboratory metrics, focusing on their diagnostic evaluation alone and in combination, in particular on their diagnostic performance for ANHC.

Gene dysregulation analysis builds a mechanistic signature for prognosis and therapeutic benefit in colorectal cancer

The implementation of cancer precision medicine requires biomarkers or signatures for predicting prognosis and therapeutic benefits. Most of current efforts in this field are paying much more attention to predictive accuracy than to molecular mechanistic interpretability. Mechanism-driven strategy has recently emerged, aiming to build signatures with both predictive power and explanatory power. Driven by this strategy, we developed a robust gene dysregulation analysis framework with machine learning algorithms, which is capable of exploring gene dysregulations underlying carcinogenesis from high-dimensional data with cooperativity and synergy between regulators and several other transcriptional regulation rules taken into consideration. We then applied the framework to a colorectal cancer (CRC) cohort from The Cancer Genome Atlas. The identified CRC-related dysregulations significantly covered known carcinogenic processes and exhibited good prognostic effect. By choosing dysregulations with greedy strategy, we built a four-dysregulation (4-DysReg) signature, which has the capability of predicting prognosis and adjuvant chemotherapy benefit. 4-DysReg has the potential to explain carcinogenesis in terms of dysfunctional transcriptional regulation. These results demonstrate that our gene dysregulation analysis framework could be used to develop predictive signature with mechanistic interpretability for cancer precision medicine, and furthermore, elucidate the mechanisms of carcinogenesis.

Identifying critical states of hepatocellular carcinoma based on landscape dynamic network biomarkers

Hepatocellular carcinoma (HCC) is the major histological form of primary liver cancer. It has usually reached the disease state once the patient is diagnosed since there are no specific symptoms in the early stages of HCC. This fact increases the difficulty of curing HCC. Recently, quantities of evidence have shown that many mathematical methods (such as dynamic network biomarkers, DNB) can be used to detect critical states or tipping points of complex diseases. However, it is difficult to apply the DNB theory to the clinic since multiple samples are generally unavailable for individual patient. This paper constructs a novel method based on landscape dynamic network biomarkers (L-DNB), which aims to detect early warning signals from cirrhosis state to very advanced HCC state in individual patient. The selected dataset contains multiple samples for each HCC state. A score that indicates the disease characteristics is calculated for each sample by RNA-seq data, and several scores constitute a distribution in the same state. Quantifying the statistical characteristics of these distributions and determining that low-grade dysplastic and high-grade dysplastic are the critical states of HCC. These results can provide scientific advice for early warning indicators and optimal treatment time for HCC.

Chronic hepatitis B: dynamic change in Traditional Chinese Medicine syndrome by dynamic network biomarkers

Background

In traditional Chinese medicine (TCM) clinical practice, TCM syndromes help to understand human homeostasis and guide individualized treatment. However, the TCM syndrome changes with disease progression, of which the scientific basis and mechanism remain unclear.

Methods

To demonstrate the underlying mechanism of dynamic changes in the TCM syndrome, we applied a dynamic network biomarker (DNB) algorithm to obtain the DNBs of changes in the TCM syndrome, based on the transcriptomic data of patients with chronic hepatitis B and typical TCM syndromes, including healthy controls and patients with liver-gallbladder dampness-heat syndrome (LGDHS), liver-depression spleen-deficiency syndrome (LDSDS), and liver-kidney yin-deficiency syndrome (LKYDS). The DNB model exploits collective fluctuations and correlations of the observed genes, then diagnoses the critical state.

Results

Our results showed that the DNBs of TCM syndromes were comprised of 52 genes and the tipping point occurred at the LDSDS stage. Meanwhile, there were numerous differentially expressed genes between LGDHS and LKYDS, which highlighted the drastic changes before and after the tipping point, implying the 52 DNBs could serve as early-warning signals of the upcoming change in the TCM syndrome. Next, we validated DNBs by cytokine profiling and isobaric tags for relative and absolute quantitation (iTRAQ). The results showed that PLG (plasminogen) and coagulation factor XII (F12) were significantly expressed during the progression of TCM syndrome from LGDHS to LKYDS.

Conclusions

This study provides a scientific understanding of changes in the TCM syndrome. During this process, the cytokine system was involved all the time. The DNBs PLG and F12 were confirmed to significantly change during TCM-syndrome progression and indicated a potential value of DNBs in auxiliary diagnosis of TCM syndrome in CHB.

Trial registration Identifier: NCT03189992. Registered on June 4, 2017. Retrospectively registered (http://www.clinicaltrials.gov)