1.中国中医科学院中药研究所 北京 100700
2.中国中医科学院广安门医院
3.天津中医药大学第一附属医院
陈文佳,女,博士
# 张彦琼,女,博士,教授,研究员,博士生导师,主要研究方向:中医药系统生物学,E-mail:yqzhang@icmm.ac.cn
纸质出版日期:2023-05-30,
网络出版日期:2023-03-15,
收稿日期:2023-01-30,
移动端阅览
陈文佳, 巩勋, 刘蔚翔, 等. 应用“GSEA-WGCNA-验证”整合策略分析类风湿关节炎肝肾亏虚证的生物标志物研究[J]. 北京中医药大学学报, 2023,46(5):599-606.
CHEN Wenjia, GONG Xun, LIU Weixiang, et al. Identification of novel biomarkers for rheumatoid arthritis with liver-kidney deficiency pattern by a "GSEA-WGCNA-validation" integrated strategy[J]. Journal of Beijing University of Traditional Chinese Medicine, 2023,46(5):599-606.
陈文佳, 巩勋, 刘蔚翔, 等. 应用“GSEA-WGCNA-验证”整合策略分析类风湿关节炎肝肾亏虚证的生物标志物研究[J]. 北京中医药大学学报, 2023,46(5):599-606. DOI: 10.3969/j.issn.1006-2157.2023.05.002.
CHEN Wenjia, GONG Xun, LIU Weixiang, et al. Identification of novel biomarkers for rheumatoid arthritis with liver-kidney deficiency pattern by a "GSEA-WGCNA-validation" integrated strategy[J]. Journal of Beijing University of Traditional Chinese Medicine, 2023,46(5):599-606. DOI: 10.3969/j.issn.1006-2157.2023.05.002.
目的
2
应用基因集富集分析(GSEA)和加权基因共表达网络分析(WGCNA)发现类风湿关节炎(RA)肝肾亏虚证的候选生物标志物,并进行临床验证,探讨RA肝肾亏虚证的生物学基础。
方法
2
临床收集RA肝肾亏虚证(3例)、7个非肝肾亏虚证(各3例)及健康志愿者(4名)的全血样本,开展转录组学测序。以健康志愿者为对照,筛选RA肝肾亏虚证的差异表达基因集进行富集分析和功能注释。以RA非肝肾亏虚证患者和健康志愿者为对照,对转录组表达谱开展GSEA和WGCNA联用挖掘,将获取的关键差异表达基因作为RA肝肾亏虚证的候选生物标志物。利用独立临床验证集样本(每组≥12例)对候选生物标志物的表达水平进行qPCR检测,采用受试者操作特征(ROC)曲线等评价其辨证效能。
结果
2
RA肝肾亏虚证的差异表达基因富集于“免疫-炎症”相关通路、细胞调控相关通路和代谢相关通路,同时,还参与肝肾发育和代谢等生物过程。转录组表达谱的GSEA富集结果表明,与非肝肾亏虚证组相比,肝肾亏虚证组的差异表达基因更明显地参与肝功能(脂质、血液)代谢调节,肾功能(水盐、激素)代谢调节和神经系统调节相关的作用通路。WGCNA分析使转录组表达谱的17 010个基因被分为了19个特征模块,其中3个特征模块与肝肾亏虚证呈明显正相关(
r
>
0.300,
P
<
0.05),其生物功能以“免疫-炎症”调节为主。整合GSEA和WGCNA分析结果后,选取变异系数、作用通路和生物模块代表性均排名前50%的3个关键基因[花生四烯酸5-脂氧合酶(ALOX5)、含Patatin样磷脂酶结构域蛋白8(PNPLA8)及抗沉默功能1(ASF1A)]作为RA肝肾亏虚证的候选生物标志物。验证结果显示:ALOX5、PNPLA8和ASF1A的辨证敏感度分别为88.89%、100.00%和100.00%,特异度分别为84.51%、76.47%和78.69%,准确度分别为85.00%、79.49%和80.00%,精确度分别为88.89%、100.00%和100.00%,ROC曲线下面积值分别为0.860、0.910和0.900。
结论
2
应用“GSEA-WGCNA-验证”整合策略,发现了RA肝肾亏虚证的新型生物标志物,有助于提高RA核心证候临床精准化诊断的水平和证候客观化研究的深度。
Objective
2
We aimed to (i) explore the biological basis of rheumatoid arthritis (RA) with liver-kidney deficiency (LKD) pattern using gene set enrichment analysis (GSEA) and weighted gene co-expression network analysis (WGCNA) method and (ii) identify and clinically validate candidate biomarkers.
Methods
2
Transcriptome sequencing was carried out on whole blood samples from RA patients with LKD pattern(3 cases)
RA patients with other seven TCM patterns(3 cases each)
and healthy volunteers(4 persons). Differentially expressed gene (DEG) sets of RA with LKD pattern were screened using healthy control samples and other TCM patterns as controls. Then
biological functions of DEGs were investigated by enrichment analysis and functional annotation. After that
the gene expression profiles were mined by GSEA and WGCNA to obtain key DEGs as candidate biomarkers for RA with LKD pattern. The expression levels of the candidate biomarkers were experimentally determined by qPCR using an independent clinical cohort (not less than 12 cases/group)
and their clinical efficacy was assessed using receiver operating characteristic (ROC)curve analysis.
Results
2
DEGs of RA with LKD pattern were most significantly enriched in the "inflammation-immune"
cell regulation
and metabolic pathways
and also involved in biological processes such as liver and kidney development and metabolism. The GSEA result of the gene expression profiles indicated that the DEGs of RA with LKD pattern were more significantly involved in hepatic function (lipid
blood) metabolic regulation-
renal function (water
salt
hormone) metabolic regulation-
and neurological regulation-related pathways than those of RA with other TCM patterns. The expression profiles of 17 010 genes were categorized into 19 functional modules through WGCNA
three of which were significantly positively correlated with LKD (
r
>
0.300
P
<
0.05)
and their biological functions mainly included "immune-inflammatory" regulation. After integrating the GSEA and WGCNA result
three key genes (ALOX5
PNPLA8
ASF1A) that ranked in the top 50% in terms of coefficient of variation and representativeness of pathway and biological modules were selected as candidate biomarkers for RA with LKD pattern. Further validation of the clinical independent samples and evaluation of the ROC model showed that the sensitivities of ALOX5
PNPLA8
and ASF1A were 88.89%
100.00%
and 100.00%
their specificities were 84.51%
76.47%
and 78.69%
their accuracies were 85.00%
79.49%
and 80.00%
their precision was 88.89%
100.00%
and 100.00%
and their values of area of ROC curve were 0.860
0.910
and 0.900
respectively.
Conclusion
2
This study applied the "GSEA-WGCNA-validation" integration strategy to identify novel biomarkers of RA with LKD patterns. The validation result of the independent sample set showed that they have good clinical efficacy
which may help improve the accuracy of clinical diagnosis of core RA patterns and the depth of objective research on the TCM patterns.
类风湿关节炎肝肾亏虚证转录组学基因集富集分析加权基因共表达网络分析花生四烯酸5-脂氧合酶含Patatin样磷脂酶结构域蛋白8抗沉默功能1
rheumatoid arthritisliver-kidney deficiency patterntranscriptomicsgene set enrichment analysisweighted gene co-expression network analysisALOX5PNPLA8ASF1A
陈灏珠. 实用内科学[M]. 12版. 北京:人民卫生出版社,2005:36-38.
黄高孝. 基于真实世界的类风湿关节炎中医证治规律研究[D]. 成都:成都中医药大学,2020.
许文亚,刘秋红,李晓茹. 类风湿关节炎从虚论治[J]. 中外医疗,2016, 35(34):195-198.
孟庆红,曲淑琴. 从肝肾论治类风湿关节炎[J]. 光明中医,2014, 29(4):713-714.
RAKSASUK S, UNGPRASERT P. Patients with rheumatoid arthritis have an increased risk of incident chronic kidney disease: a systematic review and meta-analysis of cohort studies[J]. Int Urol Nephrol, 2020, 52(1):147-154.
KOCHI M, KOHAGURA K, SHIOHIRA Y, et al. Chronic kidney disease, inflammation, and cardiovascular disease risk in rheumatoid arthritis[J]. J Cardiol, 2018, 71(3):277-283.
TOKOROYAMA T, ANDO M, SETOGUCHI K, et al. Prevalence, incidence and prognosis of chronic kidney disease classified according to current guidelines: a large retrospective cohort study of rheumatoid arthritis patients[J]. Nephrol Dial Transplant, 2017, 32(12):2035-2042.
BARBARROJA N, RUIZ-PONCE M, CUESTA-LÓPEZ L, et al. Nonalcoholic fatty liver disease in inflammatory arthritis: Relationship with cardiovascular risk[J/OL]. Front Immunol, 2022, 13: 997270[2022-11-07]. https://www.frontiersin.org/articles/10.3389/fimmu.2022.997270/fullhttps://www.frontiersin.org/articles/10.3389/fimmu.2022.997270/full.
KOCHI M, KOHAGURA K, SHIOHIRA Y, et al. Inflammation as a Risk of Developing Chronic Kidney Disease in Rheumatoid Arthritis [J/OL]. PLoS One, 2016, 11(8): 0160255[2022-11-07].https://journals.plos.org/plosone/article?id=10.1371/journal.pone.0160225https://journals.plos.org/plosone/article?id=10.1371/journal.pone.0160225.
ROBLES JA, QURESHI SE, STEPHEN SJ, et al. Efficient experimental design and analysis strategies for the detection of differential expression using RNA-sequencing [J/OL]. BMC Genomics, 2012, 13: 484[2022-11-07]. https://bmcgenomics.biomedcentral.com/articles/10.1186/1471-2164-13-484https://bmcgenomics.biomedcentral.com/articles/10.1186/1471-2164-13-484.
HANSEN KD, WU Z, IRIZARRY RA, et al. Sequencing technology does not eliminate biological variability[J]. Nat Biotechnol, 2011, 29(7):572-573.
ARNETT FC, EDWORTHY SM, BLOCH DA, et al. The American Rheumatism Association 1987 revised criteria for the classification of rheumatoid arthritis[J].Arthritis Rheum, 1988, 31 (3):315-324.
ALETAHA D, NEOGI T, SILMAN AJ, et al. 2010 Rheumatoid arthritis classification criteria: an American College of Rheumatology/European League Against Rheumatism collaborative initiative[J]. Arthritis Rheum, 2010, 62(9):2569-2581.
中华中医药学会风湿病分会. 类风湿关节炎病证结合诊疗指南[J].中医杂志,2018, 59(20):1794-1800.
TODD EV, BLACK MA, GEMMELL NJ. The power and promise of RNA-seq in ecology and evolution[J]. Molecul Ecol, 2016: 25(6), 1224-1241.
LIU YW, ZHOU J, WHITE KP. RNA-seq differential expression studies: more sequence or more replication[J]. Bioinformatics, 2014, 30(3):301-304.
陈彬.医学多因素分析设计样本例数估算:多因素分析设计样本例数综合估算法[J].伤害医学,2012,1(4):49-56.
陈文佳,巩勋,刘蔚翔,等. 从“病-证-症”关联网络探索类风湿关节炎中医证候的生物内涵[J].中国中药杂志,2022, 47(3):796-806.
PREVETE N, LIOTTI F, AMORESANO A, et al. New perspectives in cancer: Modulation of lipid metabolism and inflammation resolution[J]. Pharmacol Res, 2018, 128:80-87.
CHEN XT, XIE HW, LIU Y, et al. Interference of ALOX5 alleviates inflammation and fibrosis in high glucose induced renal mesangial cells[J/OL]. Exp Ther Med, 2022, 25(1): 34[2022-12-23].https://www.spandidos-publications.com/10.3892/etm.2022.11733https://www.spandidos-publications.com/10.3892/etm.2022.11733.
KIM HJ, KIM SW, LEE SH, et al. Inhibiting 5-lipoxygenase prevents skeletal muscle atrophy by targeting organogenesis signalling and insulin-like growth factor-1[J]. J Cachexia Sarcopenia Muscle, 2022, 13(6):3062-3077.
LIU GY, MOON SH, JENKINS CM, et al. Synthesis of oxidized phospholipids by sn-1 acyltransferase using 2-15-HETE lysophospholipids[J]. J Biol Chem, 2019, 294(26): 10146-10159.
XIE YC, SAHIN M, SINHA S, et al. SETD2 loss perturbs the kidney cancer epigenetic landscape to promote metastasis and engenders actionable dependencies on histone chaperone complexes[J]. Nat Cancer, 2022, 3(2):188-202.
0
浏览量
18
下载量
3
CSCD
关联资源
相关文章
相关作者
相关机构