Identification of novel biomarkers for rheumatoid arthritis with liver-kidney deficiency pattern by a "GSEA-WGCNA-validation" integrated strategy
Special Theme: TCM Prevention and Treatment of Rheumatism|更新时间:2023-06-05
|
Identification of novel biomarkers for rheumatoid arthritis with liver-kidney deficiency pattern by a "GSEA-WGCNA-validation" integrated strategy
Journal of Beijing University of Traditional Chinese MedicineVol. 46, Issue 5, Pages: 599-606(2023)
作者机构:
1.中国中医科学院中药研究所 北京 100700
2.中国中医科学院广安门医院
3.天津中医药大学第一附属医院
作者简介:
Prof. ZHANG Yanqiong, Ph.D., Researcher, Doctoral Supervisor. Institute of Chinese Materia Medica, China Academy of Chinese Medical Sciences, No.16, Dongzhimennei Nanxiaojie Road, Dongcheng District, Beijing 100700. E-mail: yqzhang@icmm.ac.cn
基金信息:
National Key R&D Program of China(2018YFC1705201);National Natural Science Foundation of China(81630107)
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:
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.
Identification of novel biomarkers for rheumatoid arthritis with liver-kidney deficiency pattern by a "GSEA-WGCNA-validation" integrated strategy
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.
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/full https://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.0160225 https://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-484 https://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 .
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.11733 https://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 .