李佳偉，時(shí)兆婷，郭 翌，胡雨舟，喬夢(mèng)云，余錦華，汪源源，周世崇，常 才

1. 復(fù)旦大學(xué)附屬腫瘤醫(yī)院超聲科，復(fù)旦大學(xué)上海醫(yī)學(xué)院腫瘤學(xué)系，上海 200032；2. 復(fù)旦大學(xué)電子工程系，上海 200433

超聲影像組學(xué)對(duì)浸潤(rùn)性乳腺癌激素受體表達(dá)預(yù)測(cè)價(jià)值的探索性研究

李佳偉1，時(shí)兆婷1，郭 翌2，胡雨舟2，喬夢(mèng)云2，余錦華2，汪源源2，周世崇1，常 才1

1. 復(fù)旦大學(xué)附屬腫瘤醫(yī)院超聲科，復(fù)旦大學(xué)上海醫(yī)學(xué)院腫瘤學(xué)系，上海 200032；2. 復(fù)旦大學(xué)電子工程系，上海 200433

李佳偉，復(fù)旦大學(xué)附屬腫瘤醫(yī)院主治醫(yī)師。2006年7月畢業(yè)于西安交通大學(xué)醫(yī)學(xué)院，獲取臨床醫(yī)學(xué)專業(yè)本科學(xué)歷，2009年2月碩士畢業(yè)于香港理工大學(xué)醫(yī)療科技與資訊學(xué)系，2012年博士畢業(yè)于香港中文大學(xué)麻醉學(xué)系。2013年1月至2014年7月于香港中文大學(xué)麻醉學(xué)系從事高級(jí)研究員工作。2015年7月年完成放射影像學(xué)住院醫(yī)師規(guī)范化培訓(xùn)，2016年7月完成腫瘤超聲住院醫(yī)師培訓(xùn)，現(xiàn)為復(fù)旦大學(xué)附屬腫瘤醫(yī)院超聲科主治醫(yī)師。對(duì)乳腺及甲狀腺等淺表器官常見疾病的超聲診斷積累了一定經(jīng)驗(yàn)。其教育和工作經(jīng)歷體現(xiàn)了醫(yī)工結(jié)合，多學(xué)科交叉的特點(diǎn)。主要社會(huì)任職有海醫(yī)會(huì)超聲醫(yī)學(xué)專家委員會(huì)青年委員會(huì)青年委員，美國(guó)超聲醫(yī)師協(xié)會(huì)會(huì)員。發(fā)表國(guó)內(nèi)外專業(yè)論文20多篇，會(huì)議論文30多篇，參與書籍《乳腺原位癌的超聲診斷》《Musculoskeletal Ultrasound for Regional Anaesthesia and Pain Medicine》等的編寫。

目的：探討超聲影像組學(xué)定量特征對(duì)浸潤(rùn)性乳腺癌激素受體表達(dá)的預(yù)測(cè)價(jià)值。方法：回顧性分析204例浸潤(rùn)性乳腺癌患者的術(shù)前超聲圖像及術(shù)后病理結(jié)果。根據(jù)雌激素受體(estrogen receptor，ER)、孕激素受體(progesterone receptor，PR)及人表皮生長(zhǎng)因子受體2(human epidermal growth factor receptor 2，HER-2)表達(dá)，將患者分為兩組：激素受體陽性組(ER+、PR+、HER-2-)，激素受體陰性組(ER-、PR-、HER-2-)。兩名具有5年以上臨床經(jīng)驗(yàn)的超聲科醫(yī)師對(duì)乳腺癌腫塊超聲圖像進(jìn)行回顧性特征分析與評(píng)估，評(píng)估內(nèi)容包括腫塊的形態(tài)、邊緣、內(nèi)部回聲、后方回聲改變及鈣化。然后，對(duì)同一個(gè)腫塊利用基于相位信息的動(dòng)態(tài)輪廓模型進(jìn)行邊緣分割。通過t檢驗(yàn)，篩選出與激素受體相關(guān)性最強(qiáng)的特征參數(shù)，通過支持向量機(jī)分類器，運(yùn)用徑向基核函數(shù)進(jìn)行分析與研究。為減少偏倚，采用留一法對(duì)分類性能進(jìn)行驗(yàn)證。結(jié)果：激素受體陽性組與陰性組在形態(tài)、邊緣毛刺成角、內(nèi)部回聲及后方回聲改變等二維特征方面存在顯著統(tǒng)計(jì)學(xué)差異(P＜0.05)。定量分析選出54個(gè)定量特征，對(duì)激素受體表達(dá)具有較高準(zhǔn)確率(準(zhǔn)確率67.7%，曲線下面積73.2%)。此外，邊緣、內(nèi)部回聲、后方回聲及鈣化等定量特征在激素受體陽性與陰性組之間均存在顯著統(tǒng)計(jì)學(xué)差異(P＜0.05)。結(jié)論：超聲影像組學(xué)定量特征分析降低了傳統(tǒng)超聲影像的主觀性，在預(yù)測(cè)浸潤(rùn)性乳腺癌激素受體表達(dá)方面具有較大優(yōu)勢(shì)，其在提高超聲影像學(xué)特征對(duì)乳腺癌精準(zhǔn)診斷及生物學(xué)行為預(yù)測(cè)方面的價(jià)值仍需進(jìn)一步研究。

影像組學(xué)；定量超聲特征；激素受體；浸潤(rùn)性乳腺癌

在中國(guó)乃至世界范圍內(nèi)，乳腺癌發(fā)病率位居女性癌癥的首位。根據(jù)其分子標(biāo)記的表達(dá)及生物學(xué)行為，可分為幾種臨床亞型[1-2]，這些基于分子生物學(xué)的分型是臨床上選擇最適當(dāng)治療方法的重要參考，并在個(gè)體化治療方面頗有價(jià)值[2-5]。在所有分子標(biāo)記中，雌激素受體(estrogen receptor，ER)、孕激素受體(progesterone receptor，PR)及人表皮生長(zhǎng)因子受體2 (human epidermal growth factor receptor 2，HER-2)的表達(dá)狀態(tài)在個(gè)體化治療決策中最重要，已有針對(duì)這些分子標(biāo)記陽性表達(dá)的內(nèi)分泌治療及靶向治療。然而，對(duì)于某些類型乳腺癌，如三陰性乳腺癌(triple negative breast cancer，TNBC)，以上3個(gè)受體均缺乏表達(dá)，尚無針對(duì)性治療，僅限于一般的細(xì)胞毒性化療藥物。TNBC往往腫塊較大，分化級(jí)別較高，淋巴結(jié)轉(zhuǎn)移概率較大，預(yù)后較差[4]。因此，能及早發(fā)現(xiàn)并鑒別診斷此類惡性程度高的乳腺癌對(duì)患者預(yù)后及臨床診治均有較大指導(dǎo)價(jià)值。目前，明確這些分子標(biāo)記的表達(dá)主要依靠穿刺標(biāo)本或手術(shù)切除腫塊的免疫組織化學(xué)分析。

活檢或手術(shù)標(biāo)本檢測(cè)只采取一部分腫瘤樣本組織，存在一定的局限性，因?yàn)槠浜雎粤四[瘤尤其是較大腫塊可能存在的異質(zhì)性[6]。影像學(xué)成像可提供腫瘤組織整體性解剖和功能特性，如大小、形狀、邊緣形態(tài)、內(nèi)部回聲或增強(qiáng)程度[7]。文獻(xiàn)報(bào)道，ER陽性表達(dá)的乳腺癌往往表現(xiàn)為邊緣模糊不光整及后方回聲衰減，TNBC通常表現(xiàn)為邊緣較光整和后方回聲增強(qiáng)[8-14]。但文獻(xiàn)報(bào)道結(jié)果不盡相同，且這些二維特征均為主觀性判斷，缺乏影像組學(xué)的定量分析。

目前，影像學(xué)已逐漸向自動(dòng)化數(shù)據(jù)特征分析及高通量特征提取方向發(fā)展，即影像組學(xué)[15-17]，其理論依據(jù)是不同患者腫瘤之間、同一患者不同腫瘤組織之間或同一腫瘤內(nèi)部均存在基因異質(zhì)性，這些基因異質(zhì)性表現(xiàn)為不同的影像學(xué)表型。關(guān)于乳腺癌影像組學(xué)研究主要集中于MRI紋理分析對(duì)新輔助化療的效果預(yù)測(cè)等方面[18-19]，尚未有關(guān)于乳腺癌超聲影像組學(xué)與免疫組織化學(xué)分子標(biāo)記相關(guān)性研究。本研究旨在探索定量超聲影像學(xué)特征對(duì)浸潤(rùn)性乳腺癌激素受體表達(dá)的預(yù)測(cè)價(jià)值。

1 資料和方法

1.1 研究對(duì)象

收集2014年4月—2016年6月于復(fù)旦大學(xué)附屬腫瘤醫(yī)院接受乳腺癌手術(shù)的411例患者并對(duì)其臨床資料進(jìn)行回顧性研究。收集內(nèi)容包括術(shù)前超聲圖像及描述、手術(shù)方式和術(shù)后病理診斷資料。入選標(biāo)準(zhǔn)：女性，初發(fā)及單發(fā)乳腺癌結(jié)節(jié)，未接受術(shù)前新輔助化療，術(shù)后石蠟病理證實(shí)為浸潤(rùn)性乳腺癌。根據(jù)腫塊的免疫組織化學(xué)表達(dá)，選取HER-2陰性、ER和PR均陽性或均陰性的患者，最終204例入組。所有臨床報(bào)告均由同一研究者分析校對(duì)。本研究獲得復(fù)旦大學(xué)附屬腫瘤醫(yī)院倫理委員會(huì)的批準(zhǔn)及對(duì)患者知情同意書的豁免。

1.2 儀器與二維超聲圖像特征分析

選用SUPERSONIC公司Aixplorer及GE公司Logiq E9超聲儀器，探頭頻率6～13 MHz，超聲檢查按乳腺檢查規(guī)范進(jìn)行。根據(jù)每個(gè)患者的條件調(diào)節(jié)儀器頻率、增益及時(shí)間增益補(bǔ)償(time gain compensation，TGC)曲線以達(dá)到最佳灰階超聲成像質(zhì)量，深度范圍及焦點(diǎn)根據(jù)結(jié)節(jié)的深度調(diào)節(jié)。每個(gè)乳腺腫塊的超聲圖像選取一幀以DICOM形式存儲(chǔ)。兩名具有5年以上臨床經(jīng)驗(yàn)的超聲科醫(yī)師對(duì)超聲圖像進(jìn)行回顧性特征分析與評(píng)估，其對(duì)患者其他臨床資料均不知情。對(duì)腫塊的評(píng)價(jià)指標(biāo)，按2013年乳腺影像報(bào)告和數(shù)據(jù)系統(tǒng)(Breast Imaging Reporting and Data System，BI-RADS)超聲影像詞典：形態(tài)、邊緣、內(nèi)部回聲、后方回聲改變及鈣化[20]。對(duì)腫塊定性分析出現(xiàn)不一致的患者，由兩名醫(yī)師共同討論達(dá)成一致后記錄數(shù)據(jù)資料。

1.3 超聲影像組學(xué)特征分析

乳腺癌超聲圖像由計(jì)算機(jī)輔助定量分析。每個(gè)腫塊先利用基于相位信息的動(dòng)態(tài)輪廓模型(phase-based active contour，PBAC)進(jìn)行分割[21]。該模型利用具有相位信息的能量方程，使其最小化可自動(dòng)提取腫塊的邊界信息。對(duì)某些極其不規(guī)則及對(duì)比度差的腫塊，如果自動(dòng)分割腫塊邊界的效果較差，則改用手動(dòng)分割方法來提取腫塊邊界。每一幅超聲圖像提取的特征共463個(gè)，分為四大類：形態(tài)、強(qiáng)度、質(zhì)地、小波特征[22-24]。通過t檢驗(yàn)，篩選出54個(gè)與激素受體相關(guān)性最強(qiáng)的特征參數(shù)，涵蓋乳腺癌大小、形態(tài)、邊緣、內(nèi)部回聲、后方回聲及鈣化等二維特征分類。具體描述如下：① 大?。耗[塊的面積。② 形態(tài)：腫塊的堅(jiān)固度，為腫瘤面積與其外接凸包面積的比，比值越大提示腫塊形態(tài)規(guī)則的可能性越大。③ 邊緣：5個(gè)特征用來描述腫塊的邊緣，其中邊緣粗糙度最具特征性，為歸一化的縱向長(zhǎng)度直方圖的熵[25]。粗糙度越大，邊緣不光整程度越高。④ 內(nèi)部回聲：40個(gè)基于灰度共生矩陣(gray level co-occurrence matrix，GLCM)、灰度區(qū)域矩陣(gray level size zone matrix，GLSZM)、灰度游程矩陣(gray level run length matrix，GLRLM)、局部灰度差分矩陣(neighborhood gray-tone dif f erence matrix，NGTDM)及其小波分解的圖像特征用來描述腫塊的內(nèi)部回聲，其中NGTDM粗糙度與激素受體表達(dá)最相關(guān)，值越大提示腫塊的不均質(zhì)程度越高。⑤ 后方回聲：腫塊后方回聲與周圍區(qū)域回聲的差值。⑥ 鈣化：6個(gè)基于直方圖的特征，其中定量分析中小波分解垂直分量的直方圖能量相關(guān)性最高。

1.4 乳腺癌組織學(xué)特征

所有乳腺癌標(biāo)本記錄大小后經(jīng)石蠟包埋，進(jìn)行 4 μm 厚切片，對(duì)其進(jìn)行常規(guī)蘇木精-伊紅(hematoxylin-eosin，HE)染色及免疫組織化學(xué)檢查。HE染色標(biāo)本記錄內(nèi)容包括病理類型、組織學(xué)分級(jí)、脈管內(nèi)癌栓、乳頭內(nèi)癌栓。組織學(xué)分級(jí)分3類：Ⅰ級(jí)(高分化)、Ⅱ級(jí)(中分化)、Ⅲ級(jí)(低分化)[26-27]。ER、PR、HER-2、Ki-67表達(dá)的檢測(cè)采用SP法，實(shí)驗(yàn)步驟嚴(yán)格按說明書進(jìn)行。ER、PR及Ki-67 表達(dá)情況依細(xì)胞核著色情況進(jìn)行判定，細(xì)胞核中出現(xiàn)棕黃色顆粒者為著色陽性細(xì)胞，視野中腫瘤細(xì)胞核著色細(xì)胞數(shù)＜1%為陰性，≥1%為陽性[2]。HER-2 表達(dá)情況根據(jù)CerbB-2 結(jié)果進(jìn)行判定，結(jié)果為－或+者為 HER-2 陰性，+++者為 HER-2 陽性。++為HER-2不確定患者，進(jìn)一步用熒光原位雜交(fluorescence in situ hybridization，F(xiàn)ISH)進(jìn)行HER-2基因擴(kuò)增狀態(tài)檢測(cè)。FISH技術(shù)通過熒光標(biāo)記的DNA探針與細(xì)胞核內(nèi)的DNA靶序列雜交，在熒光顯微鏡下觀察并分析細(xì)胞核內(nèi)雜交于DNA靶序列的探針信號(hào)，以獲得細(xì)胞核內(nèi)染色體(或染色體片段)中基因狀態(tài)的信息。所有病理信息及免疫組織化學(xué)信息均從患者術(shù)后病理報(bào)告中獲取并記錄。Ki-67陽性表達(dá)以細(xì)胞核內(nèi)出現(xiàn)棕黃色顆粒為標(biāo)準(zhǔn)，＜20%為低表達(dá)，≥20%為高表達(dá)[2]。

根據(jù)ER、PR、HER-2 表達(dá)，將患者分為兩組：激素受體陽性組(ER+、PR+、HER-2-，包括所有Luminal A型及部分Luminal B型)，激素受體陰性組(ER-、PR-、HER-2-，即三陰性)。

1.5 統(tǒng)計(jì)學(xué)處理

統(tǒng)計(jì)學(xué)分析采用SPSS 18.0及Matlab R2015b。對(duì)連續(xù)性變量，首先用Kolmogorov-Smirnov檢驗(yàn)數(shù)據(jù)的正態(tài)性，然后根據(jù)數(shù)據(jù)分布特點(diǎn)選用均數(shù)(標(biāo)準(zhǔn)差)或中位數(shù)(四分位數(shù))。分類變量用頻數(shù)(%，百分比)表達(dá)。t檢驗(yàn)用于定量數(shù)據(jù)的比較，卡方檢驗(yàn)用于分類數(shù)據(jù)的比較，P＜0.05為差異有統(tǒng)計(jì)學(xué)意義。計(jì)算機(jī)定量分析中，t檢驗(yàn)用于篩選與激素受體表達(dá)最相關(guān)的特征，然后將篩選出來的特征輸入支持向量機(jī)(support vector machine，SVM)分類器，采用徑向基核函數(shù)進(jìn)行分析與研究。由于所使用的數(shù)據(jù)庫(kù)相對(duì)較小，為減少偏倚，采用留一法對(duì)分類性能進(jìn)行驗(yàn)證。預(yù)測(cè)效能用準(zhǔn)確率、靈敏度、特異度及受試者工作特征(receiver operating characteristic，ROC)曲線的曲線下面積(area under curve，AUC)來評(píng)估[28]。

2 結(jié) 果

2.1 患者基本信息

表1為患者基本信息、手術(shù)方式及術(shù)后病理總結(jié)。204例患者平均年齡為52.0歲(標(biāo)準(zhǔn)差11.3歲)。激素受體陽性組153例，占75%，平均年齡52.4歲(標(biāo)準(zhǔn)差11.5歲)，其中Luminal A型77例、Luminal B型76例；激素受體陰性組51例，占25%，平均年齡50.8歲(標(biāo)準(zhǔn)差10.9歲)。激素受體陽性組與陰性組患者年齡差異無統(tǒng)計(jì)學(xué)意義(P=0.387)。兩組患者的手術(shù)方式、病理類型、淋巴結(jié)轉(zhuǎn)移情況、脈管內(nèi)癌栓及乳頭內(nèi)癌栓均無統(tǒng)計(jì)學(xué)差異(P＞0.05)。

表1 患者基本信息、手術(shù)方式及術(shù)后病理

兩組患者的腫瘤大小有統(tǒng)計(jì)學(xué)差異，激素受體陰性組腫塊(均值2.4 cm)大于陽性組(均值2.1 cm，P=0.04)。激素受體陽性組腫瘤76.1%為組織學(xué)Ⅰ～Ⅱ級(jí)，23.9%為組織學(xué)Ⅲ級(jí)；激素受體陰性組腫瘤18.0%為組織學(xué)Ⅰ～Ⅱ級(jí)，82.0%為組織學(xué)Ⅲ級(jí)，差異有統(tǒng)計(jì)學(xué)意義(P＜0.05)。此外，激素受體陰性組中Ki-67高表達(dá)患者比例(96.1%)明顯高于激素受體陽性組(49.7%，P＜0.05)。

2.2 乳腺癌超聲特征分析及其與激素受體表達(dá)的關(guān)系

表2總結(jié)了乳腺癌腫塊形態(tài)、邊緣、內(nèi)部回聲、后方回聲、鈣化等定性特征在激素受體陽性組和陰性組患者中的分布。兩組在形態(tài)、邊緣毛刺成角、內(nèi)部回聲及后方回聲改變等二維特征方面存在顯著統(tǒng)計(jì)學(xué)差異 (P＜0.05)。運(yùn)用定量超聲圖像分析為每個(gè)二維特征選取一個(gè)最具代表性的定量特征并與激素受體表達(dá)相關(guān)聯(lián)(表3)。除外形態(tài)，所有特征性指標(biāo)在激素受體陽性與陰性組之間均存在顯著統(tǒng)計(jì)學(xué)差異(P＜0.05)。表4總結(jié)了所有定量超聲特征對(duì)激素受體表達(dá)的綜合評(píng)估(準(zhǔn)確率57.4%～66.7%，AUC 56.0%～72.2%)，內(nèi)部回聲的40個(gè)代表性特征對(duì)激素受體表達(dá)預(yù)測(cè)具有最高的準(zhǔn)確率(66.7%)及AUC (72.2%)，另外54個(gè)定量特征的聯(lián)合分析更提高了對(duì)激素受體表達(dá)預(yù)測(cè)的準(zhǔn)確率(67.7%)及AUC (73.2%)。

二維特征結(jié)果如下：① 大?。菏荏w陽性組腫塊的病理標(biāo)本小于陰性組，有顯著性差異[(陽性組(2.1±0.9) cm，陰性組(2.4±1.0) cm，P=0.04)]。定量分析也發(fā)現(xiàn)陽性組腫塊面積小于陰性組(P=0.004)。② 形態(tài)：不規(guī)則形腫塊在激素受體陽性組和陰性組中分別占91.5%、74.5%，有統(tǒng)計(jì)學(xué)差異(P=0.007)，提示激素受體陰性組腫塊更易表現(xiàn)為形態(tài)規(guī)則。定量分析中，兩組的堅(jiān)固度無統(tǒng)計(jì)學(xué)差異(P=0.181)。③邊緣：腫塊邊緣不光整在激素受體陽性組和陰性組中分別占94.1%、90.2%，差異無統(tǒng)計(jì)學(xué)意義(P=0.337)。但邊緣出現(xiàn)毛刺成角或微小分葉的腫塊在激素受體陽性組所占比例(26.8%)高于陰性組(3.9%)，差異有統(tǒng)計(jì)學(xué)意義(P=0.001)。定量分析中，激素受體陽性組邊緣粗糙度高于受體陰性組(P=0.015)。④ 內(nèi)部回聲：兩組中大部分腫塊表現(xiàn)為低回聲，實(shí)質(zhì)混合性及囊實(shí)性腫塊在激素受體陰性組出現(xiàn)的比例(17.6%)高于陽性組(8.5%，P=0.042)。定量分析中，激素受體陽性組腫塊粗糙度高于受體陰性組(P＜0.001)。⑤ 后方回聲改變：兩組大部分腫塊表現(xiàn)為后方回聲無改變，但激素受體陽性組后方回聲衰減的比例(16.3%)高于后方回聲增強(qiáng)(10.5%)，激素受體陰性組后方回聲增強(qiáng)的比例(33.3%)高于后方回聲衰減(3.9%)，差異均有統(tǒng)計(jì)學(xué)意義(P＜0.000 5)。定量分析中，激素受體陽性組腫塊后方回聲低于周圍區(qū)域，而陰性組腫塊后方回聲高于周圍區(qū)域(P＜0.000 1)。⑥ 鈣化：鈣化在激素受體陽性組和陰性組出現(xiàn)的比例分別為25.5%、25.5%，差異無統(tǒng)計(jì)學(xué)意義(P=1.0)。但定量分析中，激素受體陽性組中小波分解垂直分量的直方圖能量低于陰性組(P=0.006)。

表2 二維超聲圖像特征與激素受體表達(dá)的關(guān)系

表3 定量超聲圖像特征與激素受體表達(dá)的關(guān)系

表4 運(yùn)用留一法對(duì)乳腺癌腫塊各個(gè)形態(tài)特征的分類性能評(píng)價(jià)

3 討 論

隨著化療、內(nèi)分泌治療及靶向治療的發(fā)展，乳腺癌的死亡率明顯下降。除治療方法的發(fā)展，影像學(xué)檢查對(duì)乳腺癌的早期發(fā)現(xiàn)也有不可忽視的作用。影像學(xué)檢查為乳腺癌診斷首選方法，但臨床應(yīng)用中其僅提供主觀性、半定量信息。目前，影像學(xué)的發(fā)展已趨向于對(duì)腫瘤圖像進(jìn)行自動(dòng)化數(shù)據(jù)特征分析及高通量特征提取，即影像組學(xué)。通過大量自動(dòng)化數(shù)據(jù)特征化算法，將感興趣區(qū)的影像數(shù)據(jù)轉(zhuǎn)化為具有高分辨率的可發(fā)掘的特征空間數(shù)據(jù)。影像組學(xué)研究主要集中于對(duì)CT、MRI及PET/CT定量圖像的定量特征分析[29-33]。研究表明，其可在很多方面輔助臨床決策，如輔助前列腺癌明確診斷[30]，預(yù)測(cè)肝癌[31]、腦膠質(zhì)細(xì)胞瘤[34]及肺癌的預(yù)后[15]，輔助肝癌個(gè)體化治療選擇[33]，評(píng)估乳腺癌新輔助化療效果[18-19]等。本研究初步探索了超聲影像組學(xué)特征分析在浸潤(rùn)性乳腺癌中的應(yīng)用，通過分析超聲影像組學(xué)特征與分子標(biāo)記的關(guān)系，發(fā)現(xiàn)超聲影像組學(xué)定量特征對(duì)激素受體表達(dá)有一定預(yù)測(cè)價(jià)值。

本研究發(fā)現(xiàn)激素受體陽性組腫塊較陰性組小，組織學(xué)分級(jí)低，Ki-67表達(dá)水平低，與文獻(xiàn)報(bào)道一致[8,13]。臨床研究也表明，激素受體陽性及HER-2陰性乳腺癌亞型的預(yù)后較好，激素受體陰性及HER-2陰性乳腺癌亞型的預(yù)后較差[4]。分子分型是乳腺癌治療及預(yù)后預(yù)測(cè)的重要參考，在個(gè)體化治療方面頗有價(jià)值[2-3]。分子標(biāo)記表達(dá)與否及表達(dá)水平可導(dǎo)致腫瘤圖像特征差異[14,35-36]，這些均構(gòu)成超聲影像組學(xué)的發(fā)展基礎(chǔ)。

本研究利用超聲影像組學(xué)特征分析，對(duì)浸潤(rùn)性乳腺癌的形態(tài)、邊緣、內(nèi)部回聲、后方回聲、鈣化特征進(jìn)行定量分析。二維超聲圖像中，腫塊的特征分類均為主觀性，例如形態(tài)只分規(guī)則及不規(guī)則，邊緣只分光整及不光整，內(nèi)部回聲及后方回聲判讀帶有很大的主觀性。而影像組學(xué)定量分析可克服二維超聲主觀性強(qiáng)的缺點(diǎn)，將腫瘤的圖像特征定量性地提取出來，尤其經(jīng)小波分析可提取隱藏在超聲二維圖像信號(hào)中的信息。本研究顯示，大部分定量特征在激素受體表達(dá)陽性組與陰性組之間有顯著性差異，所有特征的結(jié)合增加了對(duì)激素受體表達(dá)的預(yù)測(cè)價(jià)值，充分顯示了影像組學(xué)定量特征分析相較于傳統(tǒng)超聲的優(yōu)勢(shì)所在。激素受體陽性組更多表現(xiàn)為形態(tài)不規(guī)則、邊緣不光整、內(nèi)部回聲偏高、后方回聲衰減的腫塊，而陰性組更多表現(xiàn)為形態(tài)規(guī)則、邊緣光整、內(nèi)部低回聲、后方回聲增強(qiáng)的腫塊，與文獻(xiàn)報(bào)道基本一致[8,11-12,14,37-38]。

激素受體陰性組腫塊更多表現(xiàn)為邊緣光整[11-12]，激素受體陽性組腫塊更多表現(xiàn)為邊緣毛刺成角或微小分葉[8,14]?？赡茉?yàn)榧に厥荏w陽性腫塊往往表現(xiàn)為低級(jí)別，細(xì)胞增殖慢，腫瘤與宿主相互作用，周圍纖維結(jié)締組織增生形成，導(dǎo)致腫塊邊緣不光整，甚至毛刺成角；相反，激素受體陰性腫塊往往為高級(jí)別，細(xì)胞增殖快，腫瘤細(xì)胞與宿主缺乏相互作用[13]。腫塊內(nèi)部回聲在二維超聲中很難準(zhǔn)確評(píng)估與分類，尤其是一些內(nèi)部回聲不均勻的腫塊，導(dǎo)致文獻(xiàn)報(bào)道不一致[8,11,13]。計(jì)算機(jī)輔助技術(shù)在評(píng)估腫塊紋理特征方面獨(dú)具優(yōu)勢(shì)，可提取出40個(gè)反映腫塊內(nèi)部回聲的影像學(xué)特征。腫塊的粗糙程度在受體陽性組高于陰性組，結(jié)果的可重復(fù)性有待于大樣本研究來驗(yàn)證。后方回聲增強(qiáng)與ER陽性及低級(jí)別腫瘤密切相關(guān)[8,11-13,35]，可能原因是低級(jí)別腫瘤生長(zhǎng)慢，能引起更多間質(zhì)反應(yīng)及結(jié)蹄組織增生，產(chǎn)生多個(gè)不同的聲阻抗差別較大的反射界面，聲反射及衰減腫塊后方回聲失落[14]。

目前，乳腺癌的個(gè)體化治療主要著眼于應(yīng)用基因組學(xué)和蛋白組學(xué)等組學(xué)技術(shù)檢測(cè)其分子學(xué)特性，但基因組學(xué)與蛋白組學(xué)常規(guī)用于臨床實(shí)踐仍面臨著很大挑戰(zhàn)。臨床對(duì)腫瘤治療效果的評(píng)價(jià)主要依賴腫瘤大小的變化，忽略了其細(xì)胞代謝及增殖的變化。此外，由于腫瘤發(fā)展過程中空間和時(shí)間的異質(zhì)性，這類方法的局限性愈加凸顯，因?yàn)樾杌顧z以獲取標(biāo)本并分析小部分腫瘤組織，而這部分組織不能代表全部腫瘤組織特點(diǎn)。影像組學(xué)從某種程度上可彌補(bǔ)蛋白組學(xué)和基因組學(xué)的局限性，不但能為整個(gè)腫瘤提供全面的觀察視角，還能基于進(jìn)展期的腫瘤及其對(duì)治療的反應(yīng)進(jìn)行持續(xù)監(jiān)測(cè)，因此至關(guān)重要。2015年，St. Gallen國(guó)際乳腺癌研討會(huì)提出采用ER、PR、HER-2、Ki-67作為乳腺癌的臨床分型方法，但不夠完善，主要爭(zhēng)議在于對(duì)Luminal型的區(qū)分。目前，主要采用Ki-67來反映腫瘤增殖狀態(tài)，但界定其表達(dá)的高低非常困難。對(duì)此類乳腺癌，臨床決策遇到了瓶頸。對(duì)于Luminal分型，乳腺癌專家們認(rèn)為多基因譜的方法更為可信。如果能將乳腺癌影像組學(xué)與基因表達(dá)相關(guān)聯(lián)，將為臨床醫(yī)師的診療決策提供更多的信息。因此，進(jìn)一步的研究將對(duì)乳腺癌基因進(jìn)行測(cè)序，以期望提高影像學(xué)特征對(duì)腫瘤生物學(xué)行為預(yù)測(cè)的特異度。

本研究存在一定的不足，僅包括HER-2陰性表達(dá)患者，以后的研究將納入更多乳腺癌免疫組織化學(xué)分子分型，并與定量超聲特征相關(guān)聯(lián)；另外，本研究為探索性研究，樣本量較小，需大樣本來進(jìn)一步驗(yàn)證超聲影像組學(xué)在乳腺癌精準(zhǔn)診斷及個(gè)體化治療中的價(jià)值。

超聲影像組學(xué)定量特征分析降低了傳統(tǒng)超聲影像的主觀性，在預(yù)測(cè)浸潤(rùn)性乳腺癌的激素受體表達(dá)方面具有較大優(yōu)勢(shì)，預(yù)期能提高超聲影像學(xué)特征對(duì)乳腺癌精準(zhǔn)診斷及生物學(xué)行為預(yù)測(cè)的價(jià)值，但仍需進(jìn)一步研究。

[1] CURTIS C, SHAH S P, CHIN S F, et al. The genomic and transcriptomic architecture of 2 000 breast tumours reveals novel subgroups [J]. Nature, 2012, 486(7403): 346-352.

[2] GOLDHIRSCH A, WINER E P, COATES A S, et al. Personalizing the treatment of women with early breast cancer: highlights of the St Gallen International Expert Consensus on the Primary Therapy of Early Breast Cancer 2013 [J]. Ann Oncol, 2013, 24(9): 2206-2223.

[3] BAGARIA S P, RAY P S, SIM M S, et al. Personalizing breast cancer staging by the inclusion of ER, PR, and HER2 [J]. JAMA Surg, 2014, 149(2): 125-129.

[4] VODUC K D, CHEANG M C, TYLDESLEY S, et al. Breast cancer subtypes and the risk of local and regional relapse [J]. J Clin Oncol, 2010, 28(10): 1684-1691.

[5] GABOS Z, THOMS J, GHOSH S, et al. The association between biological subtype and locoregional recurrence in newly diagnosed breast cancer [J]. Breast Cancer Res Treat, 2010, 124(1): 187-194.

[6] HARRIS J R, LIPPMAN M E, MORROW M, et al. Diseases of the breast [M]. 5th ed. Philadelphia: Wolters Kluwer/Lippincott Williams & Wilkins Healths, 2014.

[7] KITAGAWA K, SAKUMA H, ISHIDA N, et al. Contrast-enhanced high-resolution MRI of invasive breast cancer: correlation with histopathologic subtypes [J]. AJR Am J Roentgenol, 2004, 183(6): 1805-1809.

[8] ZHANG L, LIU Y J, JIANG S Q, et al. Ultrasound utility for predicting biological behavior of invasive ductal breast cancers [J]. Asian Pac J Cancer Prev, 2014, 15(19): 8057-8062.

[9] ZHANG L, LI J, XIAO Y, et al. Identifying ultrasound and clinical features of breast cancer molecular subtypes by ensemble decision [J]. Sci Rep, 2015, 5: 11085.

[10] YANG Q, LIU HY, LIU D, et al. Ultrasonographic features of triple-negative breast cancer: a comparison with other breast cancer subtypes [J]. Asian Pac J Cancer Prev, 2015, 16(8): 3229-3232.

[11] AHO M, IRSHAD A, ACKERMAN S J, et al. Correlation of sonographic features of invasive ductal mammary carcinoma with age, tumor grade, and hormone-receptor status [J]. J Clin Ultrasound, 2013, 41(1): 10-17.

[12] CELEBI F, PILANCI K N, ORDU C, et al. The role of ultrasonographic findings to predict molecular subtype, histologic grade, and hormone receptor status of breast cancer [J]. Diagn Interv Radiol, 2015, 21(6): 448-453.

[13] COSTANTINI M, BELLI P, BUFI E, et al. Association between sonographic appearances of breast cancers and their histopathologic features and biomarkers [J]. J Clin Ultrasound, 2016, 44(1): 26-33.

[14] IRSHAD A, LEDDY R, PISANO E, et al. Assessing the role of ultrasound in predicting the biological behavior of breast cancer [J]. AJR Am J Roentgenol, 2013, 200(2): 284-290.

[15] AERTS H, VELAZQUEZ E R, LEIJENAAR R T H, et al. Decoding tumour phenotype by noninvasive imaging using a quantitative radiomics approach [J]. Nat Commun, 2014, 5: 4006.

[16] KUMAR V, GU Y, BASU S, et al. Radiomics: the process and the challenges [J]. Magn Reson Imaging, 2012, 30(9): 1234-1248.

[17] GILLIES R J, KINAHAN P E, HRICAK H. Radiomics: Images are more than pictures, they are data [J]. Radiology, 2016, 278(2): 563-577.

[18] TERUEL J R, HELDAHL M G, GOA P E, et al. Dynamic contrast-enhanced MRI texture analysis for pretreatment prediction of clinical and pathological response to neoadjuvant chemotherapy in patients with locally advanced breast cancer [J]. NMR Biomed, 2014, 27(8): 887-896.

[19] GOLDEN D I, LIPSON J A, TELLI M L, et al. Dynamiccontrast-enhanced MRI-based biomarkers of therapeutic response in triple-negative breast cancer [J]. J Am Med Inform Assoc, 2013, 20(6): 1059-1066.

[20] 周建橋, 詹維偉. 乳腺超聲影像報(bào)告與數(shù)據(jù)系統(tǒng)解讀[M]. 北京: 人民衛(wèi)生出版社, 2015.

[21] CAI L Y, WANG Y Y. A phase-based active contour model for segmentation of breast ultrasound images [C]. Hangzhou: [出版者不像], 2013.

[22] HUANG Y L, CHEN D R, JIANG Y R, et al. Computer-aided diagnosis using morphological features for classifying breast lesions on ultrasound [J]. Ultrasound Obstet Gynecol, 2008, 32(4): 565-572.

[23] CHEN C M, CHOU Y H, HAN K C, et al. Breast lesions on sonograms: computer-aided diagnosis with nearly setting-independent features and artificial neural networks [J]. Radiology, 2003, 226(2): 504-514.

[24] SHAN J, ALAM S K, GARRA B, et al. Computeraided diagnosis for breast ultrasound using computerized BI-RADS features and machine learning methods [J]. Ultrasound Med Biol, 2016, 42(4): 980-988.

[25] CHOU Y H, TIU C M, HUNG G S, et al. Stepwise logistic regression analysis of tumor contour features for breast ultrasound diagnosis [J]. Ultra Med Biol, 2001, 7: 1493-1498.

[26] TAVASSOLI F A, DEVILEE P. World Health Organization classification of tumours. Pathology and genetics of tumours of the breast and female genital organs [M]. Lyon: IARC Press, 2003.

[27] ELSTON C W. Classification and grading of invasive breast carcinoma [J]. Verh Dtsch Ges Pathol, 2005, 89: 35-44.

[28] BALDI P, BRUNAK S, CGAYVIN Y, et al. Assessing the accuracy of prediction algorithms for classification: an overview [J]. Bioinformatics, 2000, 16: 412-424.

[29] COROLLER T P, GROSSMANN P, HOU Y, et al. CT-based radiomic signature predicts distant metastasis in lung adenocarcinoma [J]. Radiother Oncol, 2015, 114(3): 345-350.

[30] WIBMER A, HRICAK H, GONDO T, et al. Haralick texture analysis of prostate MRI: utility for differentiating non-cancerous prostate from prostate cancer and differentiating prostate cancers with different Gleason scores [J]. Eur Radiol, 2015, 25(10): 2840-2850.

[31] SEGAL E, SIRLIN C B, OOI C, et al. Decoding global gene expression programs in liver cancer by noninvasive imaging [J]. Nat Biotechnol, 2014, 25(6): 675-680.

[32] WANG J, KATO F, OYAMA-MANABE N, et al. Identifying triple-negative breast cancer using background parenchymal enhancement heterogeneity on dynamic contrast-enhanced MRI: A pilot radiomics study [J]. PLoS One, 2015, 10(11): e0143308.

[33] KUO M D, GOLLUB J, SIRLIN C B, et al. Radiogenomic analysis to identify imaging phenotypes associated with drug response gene expression programs in hepatocellular carcinoma [J]. J Vasc Interv Radiol, 2007, 18(7): 821-831.

[34] DIEHN M, NARDINI C, WANG D S, et al. Identification of noninvasive imaging surrogates for brain tumor gene-expression modules [J]. Proc Natl Acad Sci USA, 2008, 105(13): 5213-5218.

[35] KO E S, LEE B H, KIM H A, et al. Triple-negative breast cancer: correlation between imaging and pathological findings [J]. Eur Radiol, 2010, 20(5): 1111-1117.

[36] AU-YONG I T, EVANS A J, TANEJA S, et al. Sonographic correlations with the new molecular classification of invasive breast cancer [J]. Eur Radiol, 2009, 19(10): 2342-2348.

[37] 張曉曉, 周建橋, 朱櫻, 等. 乳腺癌超聲征象與分子亞型相關(guān)性的研究 [J]. 診斷學(xué)理論與實(shí)踐, 2011, 10(2): 153-157.

[38] 陳慧琪, 郭順華, 過新民. 乳腺癌超聲征象與分子分型相關(guān)性分析 [J]. 中國(guó)超聲醫(yī)學(xué)雜志, 2015, 31(6): 498-500.

Value of ultrasound radiomics analysis in predicting expression of hormone receptors in invasive breast carcinoma

LI Jiawei1, SHI Zhaoting1, GUO Yi2, HU Yuzhou2, QIAO Mengyun2, YU Jinhua2, WANG Yuanyuan2, ZHOU Shichong1, CHANG Cai1
(1. Department of Medical Ultrasound, Fudan University Shanghai Cancer Center; Department of Oncology, Shanghai Medical College, Fudan University, Shanghai 200032, China; 2. Department of Electronic Engineering, Fudan University, Shanghai 200433, China)

CHANG Cai E-mail: changcai@163.com

Objective:To investigate the automatic radiomics approach in predicting the association between quantitativeultrasound features and hormone receptor status in invasive breast carcinoma.Methods:A total of 204 patients who accepted breast cancer surgery were retrospectively reviewed for pre-operative ultrasound images and post-operative pathological reports. Based on the expressions of estrogen receptor (ER), progesterone receptor (PR) and human epidermal growth factor receptor 2 (HER-2), all patients were divided into hormone receptor positive (ER+, PR+, HER-2-) group and hormone receptor negative (ER-, PR-, HER-2-) group. Two dimensional features of shape, margin, echo pattern, posterior acoustic feature and calcif i cation were assessed by two experienced radiologists and correlated with hormone receptor status. The same ultrasound images were then segmented for using a phase-based active contour model. The high-throughput radiomics features were extracted based on the two dimensional sonographic features. Target features were selected using Students’ t-test. The support vector machine classif i er with radial basis function and leave-oneout-cross-validation were used to correlate quantitative sonographic features with the status of hormone receptor.Results:The two groups had signif i cant dif f erences in the objective sonographic characteristics of shape, angular/spiculated margin, echo pattern and posterior acoustic feature (P＜0.05). In the quantitative radiomics analysis, 54 features were selected with high accuracy in predicting the status of hormone receptor (accuracy 67.7%, area under the curve 73.2%). In addition, in the quantitative analysis, the two groups showed signif i cant dif f erence in margin, echo pattern, posterior acoustic pattern and calcif i cation (P＜0.05).Conclusion:The quantitative features of ultrasound radiomics were well correlated with hormone receptor status of invasive breast carcinoma. Further study is warranted to validate its value in the precise diagnosis and biological behavior prediction of breast carcinoma.

Radiomics; Quantitative ultrasonographic feature; Hormone receptor; Invasive breast carcinoma

R445.1

A

1008-617X(2017)02-0128-08

2017-04-01）

國(guó)家自然科學(xué)基金項(xiàng)目(No：61401102、81371575、81627804)

常才 E-mail：changcai@163.com