2024 , Vol. 21 >Issue 03: 319 - 326
DOI: https://doi.org/10.3877/cma.j.issn.1672-6448.2024.03.011
基于超声纹理影像转录组学预测前列腺癌
Copy editor: 吴春凤
收稿日期: 2023-03-27
网络出版日期: 2024-06-05
基金资助
国家自然科学基金(81471682)
中国博士后基金特别资助项目(2021T140795)
陕西省自然科学基础研究计划(2023-JC-QN-0912)
西安市科技计划项目(21YXYJ0134)
版权
A radiotranscriptomics approach for prediction of prostate cancer based on ultrasound image texture features
Received date: 2023-03-27
Online published: 2024-06-05
Copyright
探讨灰阶超声图像和超声造影图像的影像转录组学分析方法在前列腺结节良恶性预测诊断中的应用价值。
分析2020年12月至2022年3月在解放军总医院第一医学中心接受超声引导下可疑前列腺癌进行穿刺活检的68例患者。患者均采用经直肠灰阶超声及超声造影2种显像模式,并对图像进行自动分割及纹理特征分析。活检标本进行RNA测序及前列腺癌相关的基因表达谱、功能富集和通路分析。绘制随机森林、贝叶斯和支持向量机3种数据集模型的受试者操作特征(ROC)曲线和校准曲线评价模型的预测效能。
影像组学分析得到2个关键纹理特征:ca2-GLSZM-LZHGE和GLSZM-ZSV。RNA测序发现120个与前列腺癌相关的差异基因,利用影像转录组学方法得到区分前列腺癌和前列腺增生的生物标志物:ITGB3、CAV1、miR141-3p、let-7a-5p、miR25-5p和miR200c-3p,功能富集和通路分析发现上述生物标志物与雄激素受体状态、耐药性、增殖和凋亡相关的转录组改变有关。随机森林、贝叶斯和支持向量机3种模型联合数据集ROC的曲线下面积(AUC)分别为0.99、0.98和0.99,3种模型影像组学数据集AUC分别为0.99、0.95和0.99,分别优于临床数据集AUC(0.79、0.85和0.92)及分子生物标志物数据集(转录组学)AUC(0.66、0.80和0.86)。联合数据集组模型ROC曲线和校准曲线均显示模型区分度和准确度良好。
超声图像纹理特征在评估前列腺癌的生物标志物方面具有潜在的应用价值,并且基于超声图像构建的影像转录组学联合模型比影像组学模型具有更好的预测效能。
杨倩 , 李秋洋 , 李楠 , 罗渝昆 , 唐杰 . 基于超声纹理影像转录组学预测前列腺癌[J]. 中华医学超声杂志(电子版), 2024 , 21(03) : 319 -326 . DOI: 10.3877/cma.j.issn.1672-6448.2024.03.011
To evaluate the value of radiotranscriptomics analysis of gray-scale ultrasound images and contrast-enhanced ultrasound (CEUS) images in the diagnosis of benign and malignant prostate nodules.
A total of 68 patients who underwent ultrasound-guided biopsy for suspected prostate cancer (PCa) at the First Medical Center of Chinese PLA General Hospital from December 2020 to March 2022 were analyzed. All patients underwent transrectal gray-scale ultrasound and CEUS, and the images were automatically segmented and texture features were analyzed. Biopsy specimens were subjected to RNA sequencing and prostate cancer-related gene expression profiling as well as functional enrichment and pathway analysis. Random forest, Bayesian, and support vector machine (SVM) methods were used to draw the receiver operating characteristic (ROC) curve and calibration curve to evaluate the prediction efficiency of the model.
Two key texture features, ca2-GLSZM-LZHGE and GLSZM-ZSV, were obtained by radiomics. RNA sequencing identified 120 differentially expressed genes related to PCa, and the biomarkers to distinguish PCa from benign prostatic hyperplasia (BPH) were obtained by correlation analysis: ITGB3, CAV1, miR141-3p, let-7a-5p, miR25-5p, and miR200c-3p. Functional enrichment and pathway analysis identified transcriptomic alterations associated with androgen receptor status, drug resistance, proliferation, and apoptosis. The area under the ROC curve (AUC) values of the three combined dataset models (random forest, naive Bayes, and SVM) and radiomics dataset models were 0.99, 0.98, and 0.99, and 0.99, 0.95, and 0.99, respectively, which were better than those of the clinical dataset models (0.79, 0.85, and 0.92) and molecular biomarker dataset (transcriptomics) models (0.66, 0.80, and 0.86). The ROC curve and calibration curve of the combined dataset group showed that the model had good discrimination and accuracy.
Ultrasound image texture features have potential application value in the evaluation of biomarkers of PCa, and the combined radiotranscriptomics model has better predictive efficiency than the radiomics model.
Key words: Radiotranscriptomics; Radiomics; Texture feature; Prostate cancer; Ultrasound
表1 前列腺癌及良性前列腺增生2组患者一般临床资料比较 |
资料 | 前列腺癌组(n=47) | 良性前列腺增生组(n=21) | 统计值 | P值 |
---|---|---|---|---|
年龄(岁,![]() | 69.91±7.52 | 66.67±7.39 | t=1.702 | 0.093 |
总PSA(ng/ml,![]() | 14.74±9.67 | 11.18±5.13 | t=2.199 | 0.003 |
游离PSA/总PSA(![]() | 0.11±0.08 | 0.16±0.05 | t=2.694 | 0.009 |
前列腺体积(ml,![]() | 43.79±20.30 | 68.09±31.32 | t=3.576 | 0.001 |
PSA密度(ng/ml2,![]() | 0.43±0.35 | 0.17±0.10 | t=2.651 | 0.010 |
穿刺针数(针,![]() | 12±1 | 12±0 | t=1.902 | 0.060 |
阳性穿刺针数(针,![]() | 7±3 | 0 | t=9.986 | 0.001 |
Gleason评分(例) | ||||
3+3=6 | 3 | 0 | ||
3+4=7 | 18 | 0 | ||
4+3=7 | 14 | 0 | ||
≥4+4 | 12 | 0 |
注:PSA为前列腺特异性抗原 |
图3 超声图像纹理特征相关差异基因功能与表达。图a:网络图显示灰阶超声图像和超声造影纹理特征调节mRNA及miRNA相关的主要富集的功能和KEGG通路(绿色代表mRNA,红色代表miRNA,蓝色代表功能,紫色代表通路);图b:热图及火山图显示根据关键纹理特征有2个显著差异表达的mRNA;CAV1和ITGB3(P<0.05),上端热图中橙色代表前列腺癌患者,蓝色代表良性前列腺增生患者,橙色代表差异基因上调,蓝色代表差异基因下调;下端火山图中粉色代表上调基因,绿色代表下调基因。图c:热图及火山图显示关键纹理特征有4个显著差异表达miRNA:miR141-3p、let-7a-5p、miR25-5p和miR200c-3p(P<0.05)。上端热图中橙色代表前列腺癌患者,蓝色代表良性前列腺增生患者,橙色代表差异miRAN上调,蓝色代表差异miRAN下调;下端火山图中粉色代表上调miRAN,绿色代表下调miRAN |
1 |
|
2 |
|
3 |
梁梓南, 杨薇. 影像及影像组学评价肝细胞癌微血管侵犯的应用现状[J/OL].中华医学超声杂志(电子版), 2022, 19(9): 1003-1007.
|
4 |
|
5 |
|
6 |
|
7 |
|
8 |
|
9 |
|
10 |
|
11 |
|
12 |
|
13 |
|
14 |
|
15 |
|
16 |
|
17 |
|
18 |
|
19 |
|
20 |
|
21 |
|
22 |
|
23 |
|
24 |
|
25 |
|
26 |
|
/
〈 |
|
〉 |