Integrated quantitative susceptibility mapping and radiomics for multiple system atrophy subtypes classification.

Mao Xiaohan, Zhang Di, Bai Jiayu, Liu Yumei, Yu Juncheng, et al.

View Original ↗
AI 導讀 technology general 重要性 3/5

CIT168 圖譜以 168 位在體 MRI、700 μm 解析度定義 16 個皮質下核團,並開放 OSF 眾包貢獻

  • 168 位 HCP 受試者、700 μm 等向性解析度,ANTs SyN 配準建出 CIT168 皮質下群組模板
  • 殼核 Dice 係數 0.94 vs 腹側被蓋區 0.55,核團體積越小邊界越模糊標記難度越高
  • 圖譜存放 OSF 並附 Python atlaskit 工具庫,支援眾包新核團標記持續演進

168 位健康成人(22–35 歲,84 男 84 女)的 T1w/T2w MRI 建出、空間解析度達 700 微米等向性的 CIT168 皮質下核團圖譜,把在個別腦部影像中幾乎看不見的 16 個核團邊界,全數清晰定義在群組模板上。這份 Caltech 團隊發表於 Scientific Data 的資料描述,以概率標記取代二元遮罩,並開放眾包貢獻,目標是成為在體(in vivo)人腦皮質下解剖學的持續演進資源。

建立 CIT168 的三個研究缺口

現有皮質下核團圖譜存在三個核心問題,推動了 CIT168 的誕生。其一,部分皮質下結構在單一研究的群組平均影像中根本無法清晰顯示;其二,各研究間的解剖標記命名不一致,導致跨研究比較困難;其三,因配準到標準空間可能存在空間誤差,維持概率性標記遠比使用二元遮罩更為嚴謹。

解決這些問題的策略,是從人腦連接組計畫(Human Connectome Project,HCP)S500 資料集中擷取大量個體影像,透過高精度微分同胚配準(diffeomorphic registration)平均建立群組模板,讓組織邊界比任何單一受試者的影像都更清晰。現有的 FreeSurfer 部分結構標記來自 10–15 個老年捐贈者的離體腦切片;組織學圖譜雖具細胞解析度,但死後組織的收縮與形變使其難以精確對應在體 MRI。CIT168 直接以在體資料為基礎,提供了不同性質的互補參考。

168 名受試者、700 μm 解析度:ANTs SyN 配準流程

篩選條件要求每位受試者的 T1w 與 T2w 影像各有兩次獨立掃描,以最大化訊雜比;套用此條件後,可用樣本從原始 208 人縮減,再經年齡與性別平衡後得到最終 168 位受試者。建構流程使用 ANTs(Advanced Normalization Toolbox,進階標準化工具)的 antsMultivariateTemplateConstruction.sh 腳本,執行雙變量對稱標準化(bivariate SyN)配準,T1w 與 T2w 影像被賦予相同權重,共同最佳化單一微分同胚映射,最終群組模板體素大小為 700 μm 等向性。

為建立概率標記,研究者從 168 個影像對中隨機抽取 84 對,平均建構驗證用模板,並以剩餘 84 對建立互補驗證模板,共生成 8 個驗證模板。三名觀察者在每個驗證模板的左半球分別手動標記,再利用反射翹曲(reflection warp)從左半球概率標記映射至右半球,以最小化左右腦觀察者偏差。

16 個核團的體積分布與左右腦對稱性

標記結果涵蓋 16 個皮質下灰質核團,分四組:(1)受多巴胺調控的紋狀體及輸出核,包含尾狀核(Ca)、殼核(Pu)、蒼白球(GPe/GPi)、丘腦底核(STH)及黑質網狀部(SNr);(2)驅動多巴胺中腦活動的結構,含下視丘(HTH)、韁核(HN)、腹側蒼白球(VeP)及伏隔核(NAC);(3)多巴胺中腦本身,含黑質緻密部(SNc)、臂旁色素核(PBP)及腹側被蓋區(VTA);(4)地標結構,含延伸杏仁核(EXA)、紅核(RN)及乳頭核(MN)。

體積估算以概率空間積分計算,體積最大的是殼核(Pu,左側約 5,224 μl),最小的是韁核(HN,左側約 29 μl),兩者相差逾 180 倍。多數核團左右腦體積高度對稱,側化指數普遍接近零,韁核的側化程度最高(4.5%,右偏)。

CIT168 各皮質下核團體積(μl)與左右腦側化指數
核團縮寫左側 (μl)右側 (μl)側化指數 (%)
殼核Pu522451960.3
尾狀核Ca44934543-0.6
下視丘HTH604617-1.1
蒼白球外段GPe6966781.3
蒼白球內段GPi3833830.0
伏隔核NAC397399-0.3
黑質網狀部SNr261269-1.4
紅核RN3012980.5
丘腦底核STH1351282.9
黑質緻密部SNc132136-1.5
延伸杏仁核EXA134136-1.0
臂旁色素核PBP99980.5
腹側蒼白球VeP6874-4.3
乳頭核MN64622.0
腹側被蓋區VTA33330.1
韁核HN29274.5

Dice 係數從 0.55 到 0.94:結構大小決定標記難度

評估標記可靠性採用兩個互補指標:Dice 係數(標記重疊程度,值域 0–1)及 Hausdorff 距離(兩標記間的最大最近距離,mm)。觀察者間(inter-rater)相似性數據顯示,邊界明確的大體積結構表現最佳——殼核(Pu)Dice 達 0.94,尾狀核(Ca)為 0.93;邊界模糊的小體積結構則明顯較低,腹側被蓋區(VTA)僅 0.55,臂旁色素核(PBP)為 0.62。觀察者內(intra-rater)Dice 普遍略高,VTA 觀察者內提升至 0.71,說明邊界難以辨識是各觀察者共同面臨的問題。作者指出,Dice 係數對小體積結構較不靈敏,Hausdorff 距離是評估小核團標記品質的更適當指標。

各核團觀察者間標記可靠性(Dice 係數與 Hausdorff 距離)
核團Dice 均值Dice SDHausdorff 距離均值 (mm)Hausdorff SD (mm)
Pu 殼核0.940.042.672.12
Ca 尾狀核0.930.053.643.39
RN 紅核0.920.070.790.61
GPi 蒼白球內段0.890.081.110.91
GPe 蒼白球外段0.880.091.251.03
HN 韁核0.870.100.690.58
MN 乳頭核0.850.120.810.65
STH 丘腦底核0.840.131.090.95
NAC 伏隔核0.810.152.111.85
HTH 下視丘0.810.142.512.16
SNr 黑質網狀部0.790.161.721.40
SNc 黑質緻密部0.770.171.611.36
EXA 延伸杏仁核0.750.182.112.07
VeP 腹側蒼白球0.690.251.200.97
PBP 臂旁色素核0.620.282.572.53
VTA 腹側被蓋區0.550.352.572.31

開放眾包機制與推薦應用場景

CIT168 圖譜已存放於開放科學框架(OSF),包含 MNI152 2009c 非線性空間的 T1w/T2w 模板及概率標記,並附有 CIT168 至 MNI152 的仿射與 SyN 轉換資料。Python 工具庫 atlaskit 完全開源,支援概率圖譜建立與眾包貢獻流程,研究者可下載驗證模板後提交自己的標記版本。

建議應用包括:fMRI 研究中預設感興趣區域(ROI)、測量疾病與健康狀態的皮質下體積或形態差異、定義擴散張量影像(DTI)纖維束追蹤的種子點,以及電生理記錄的目標定位。使用時需注意,極端年齡族群因髓鞘化程度差異或老年萎縮,可能影響配準精確度,應謹慎解讀。

CIT168 以 168 位在體 MRI、700 μm 解析度標記 16 個皮質下核團,以概率圖譜形式開放眾包,是組織學圖譜與在體 MRI 之間的高解析度橋接資源。