慢性肝病中医证候演变的分子特征:动态网络生物标志物分析

Molecular features of traditional Chinese medicine syndrome evolution in chronic liver diseases: a dynamic network biomarker analysis

  • 摘要:
    目的 从“异病同证”视角阐释慢性乙型肝炎(CHB)、肝硬化(LC)和肝细胞癌(HCC)患者中医证候的生物学基础,为慢性肝病(CLD)的诊断与治疗提供补充策略。
    方法 为探究中医证候在CLD中的动态演变,本研究对肝胆湿热证(LGDHS)、肝郁脾虚证(LDSDS)或肝肾阴虚证(LKYDS)的CHB、LC或HCC患者的外周血单个核细胞进行转录组学分析,受试者于2018年8月1日至2021年12月31日期间招募自上海中医药大学附属曙光医院。采用随机方差模型(RVM)F检验筛选差异表达基因(DEGs),并经错误发现率(FDR)校正。运用主成分分析(PCA)和无监督层次聚类对样本分组进行可视化。采用动态网络生物标志物(DNB)分析识别证候演变中的关键转折阶段,随后对DNB成员进行基因本体(GO)和京都基因与基因组百科全书(KEGG)通路富集分析,以阐明其功能角色及参与的通路。利用随机森林(RF)分析和受试者工作特征(ROC)曲线下面积(AUC)对候选基因的重要性进行排序。采用独立CLD队列的数据(GSE89377)以及癌症基因组图谱肝细胞癌(TCGA-LIHC)数据库中的RNA-seq数据进行外部验证。此外,在独立的LC患者队列中采用逆转录定量聚合酶链式反应(RT-qPCR)验证候选基因的表达水平。
    结果 本研究共纳入132例受试者。DNB分析表明,LDSDS阶段是CHB、LC和HCC中医证候演变的关键转折点。PI3K-AKT信号通路在所有3种CLD的DNB分析中均富集,提示其可能参与中医证候的关键转变。从PI3K-AKT通路的24个核心DNB成员中,通过RF(基尼系数 > 1)和ROC分析鉴定出4个基因:整合素亚基β1(ITGB1)、IV型胶原α1链(COL4A1)、IV型胶原α2链(COL4A2)和DNA损伤诱导转录因子3(DDIT3)。ROC分析显示,上述4个基因在区分CHB患者中的LGDHS与LKYDS时具有较高的判别能力,ITGB1的AUC为0.7891COL4A10.7070COL4A20.7148DDIT30.8945。在独立CLD队列(GSE89377)中,4个基因从正常至CHB、LC、HCC均呈显著阶梯性上调(P < 0.05)。在TCGA-LIHC数据集中,其表达随肿瘤分期逐渐升高。在独立LC队列(30例LGDHS vs. 30例LKYDS)中进行RT-qPCR实验验证,与LGDHS相比,LKYDS中ITGB1COL4A2DDIT3的表达显著上调(分别为P = 0.01520.01860.0247),而COL4A1呈上升趋势,但不显著(P = 0.1201)。
    结论 本研究为理解CLD中医证候演变的分子特征提供了一种新方法。PI3K-AKT通路及其4个基因(ITGB1COL4A1COL4A2DDIT3)在从实证(LGDHS)经关键LDSDS阶段向虚证(LKYDS)转变的过程中发挥重要作用。这些发现为中医证候分型提供了潜在的定量生物标志物和治疗靶点,并可能有助于阻断CLD的证候进展。

     

    Abstract:
    Objective To elucidate the biological basis of traditional Chinese medicine (TCM) syndromes from the perspective of “same syndrome, different diseases” in patients with chronic hepatitis B (CHB), liver cirrhosis (LC), and hepatocellular carcinoma (HCC), thereby providing a complementary approach for the diagnosis and treatment of chronic liver diseases (CLD).
    Methods To investigate the dynamic characteristics of TCM syndromes in CLD, transcriptomic profiling of peripheral blood mononuclear cells (PBMCs) was performed from patients with CHB, LC, or HCC presenting with three TCM syndromes: liver gallbladder dampness heat syndrome (LGDHS), liver depression spleen deficiency syndrome (LDSDS), and liver kidney Yin deficiency syndrome (LKYDS). These participants were recruited at Shuguang Hospital Affiliated to Shanghai University of Traditional Chinese Medicine between August 1, 2018 and December 31, 2021. Differentially expressed genes (DEGs) were identified using the random variance model (RVM) F test with false discovery rate (FDR) correction. Principal component analysis (PCA) and unsupervised hierarchical clustering were applied to visualize sample grouping. Dynamic network biomarkers (DNB) analysis was employed to detect critical transition stages during syndrome evolution, followed by Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment analyses to characterize the functional roles and pathway involvement of the DNB members. Random forest (RF) analysis and the area under the receiver operating characteristic (ROC) curves (AUC) were used to rank the importance of candidate genes. External validation was performed using microarray data from an independent CLD cohort (GSE89377) and RNA-seq data from The Cancer Genome Atlas Liver Hepatocellular Carcinoma (TCGA-LIHC) dataset. Additionally, reverse transcription quantitative polymerase chain reaction (RT-qPCR) was performed on an independent cohort of LC patients to validate the expression levels of the candidate genes.
    Results The study included a total of 132 participants. DNB analysis identified LDSDS stage as a critical tipping point during TCM syndrome evolution across CHB, LC, and HCC. The phosphoinositide 3-kinase/protein kinase B (PI3K-AKT) signaling pathway was consistently enriched in the DNB analysis across all three types of CLD, suggesting its potential involvement in the critical transition of TCM syndromes. Among the 24 core DNB members of the PI3K-AKT pathway, four genes—integrin subunit beta 1 (ITGB1), collagen type IV alpha 1 chain (COL4A1), collagen type IV alpha 2 chain (COL4A2), and DNA damage inducible transcript 3 (DDIT3)—were identified by RF analysis (Gini score > 1) and ROC analysis. ROC analysis demonstrated high discriminative ability for distinguishing LGDHS from LKYDS in CHB patients, with AUC of 0.7891 for ITGB1, 0.7070 for COL4A1, 0.7148 for COL4A2, and 0.8945 for DDIT3. In the independent CLD cohort (GSE89377), all four genes showed significant stepwise upregulation from normal to CHB, LC, and HCC (all P < 0.05). In the TCGA-LIHC dataset, their expression progressively increased with tumor stage. RT-qPCR validation in an independent LC cohort (30 LGDHS vs. 30 LKYDS) confirmed that ITGB1, COL4A2, and DDIT3 were significantly upregulated in LKYDS compared with LGDHS (P = 0.0152, 0.0186, and 0.0247, respectively), whereas COL4A1 showed a non-significant upward trend (P = 0.1201).
    Conclusion This study introduces a novel approach to understanding the molecular features underlying TCM syndrome evolution in CLD. The PI3K-AKT pathway and four identified genes (ITGB1, COL4A1, COL4A2, and DDIT3) play crucial roles in the transition from excess (LGDHS) to deficiency (LKYDS) via the critical LDSDS stage. These findings offer potential quantitative biomarkers and therapeutic targets for TCM syndrome differentiation and may help arrest syndrome progression in CLD.

     

/

返回文章
返回