劉亞，鄭罡，張龍江，盧光明

腦氧代謝率（cerebral metabolic rate of oxygen，CMRO2）表征機體組織的有氧代謝活動，是衡量腦組織的氧利用狀況的一個綜合性標志。測量此參數(shù)是研究腦組織在生理和病理狀態(tài)的重要手段，尤其在研究缺血性腦血管疾病時，該方法優(yōu)勢明顯[1]。對研究急性腦梗死疾病的發(fā)展過程、腦出血的預測以及各種內外科治療效果的判定也都有著非常重要的意義。

當前，定量CMRO2的技術主要是基于同位素自放射成像原理的正電子發(fā)射斷層顯像技術[2]（positron emission tomography，PET）。其基本過程如下：首先，吸入放射性同位素15O標記的氧氣作為示蹤劑；其次，氣態(tài)形式的15O被吸入人體后先與血液中的血紅蛋白相結合；再次，經(jīng)過人體的有氧代謝，轉換成含有15O的水；最后，放射性元素15O在人體內衰變放出正電子，正負電子湮滅后轉化為一對光子，最后光子被探測器檢測成像。然而，由于15O半衰期較短（123s），在成像過程中為了得到穩(wěn)定成像，一般需要持續(xù)給人體吸入高劑量的放射性氣體，會對人體造成較大的電離輻射傷害。

近年來，隨著磁共振代謝成像技術的發(fā)展，出現(xiàn)了一批無需外源示蹤劑的全腦及局域性代謝成像技術。不同于基于自放射成像的PET，這些技術多基于測量動靜脈含氧量差的菲克原理。這些技術目前大致可以分為兩個方案：一種方案是基于T2的靜脈血氧飽和度（venous oxygen saturation，Yv）測量，結合全腦水平的腦血流量（cerebral blood flow，CBF）來測量全腦的CMRO2；另一種方案是基于靜脈血的磁敏感值測量Yv，結合局部的CBF水平來測量局部的CMRO2。本文擬對這兩個方案及其子方案原理進行簡述。

基于菲克原理的CMRO2測量

目前，磁共振測量CMRO2是基于菲克原理，其核心是測量動脈和靜脈的氧飽和度，以及結合血流動力學，并通過離體血紅蛋白攜氧量，近似量化CMRO2（圖1），其原理如下公式所示：

其中，Ya為動脈血氧飽合度（arterial oxygen saturation，Ya），Ca是單位血液體積中的氧分子數(shù)，其常用值為833.7μmol O2／100mL[3]。Ya可以用脈搏血氧儀測量，一般接近為100%。

根據(jù)菲克原理，量化CMRO2的核心實質上是測量Yv和CBF。根據(jù)靜脈與腦組織的生理特性，測量Yv的方法可以分為根據(jù)靜脈與腦組織的 T2的差異[4－5]和磁 化 率 的差異[6－7]。根據(jù)這兩種特性，分別發(fā)展出了T2弛豫自旋標記成像（T2relaxation under spin tagging，TRUST）技術和基于相位的磁共振技術。結合相位對比技術（phase contrast，PC）或動脈自旋標記技術（arterial spin labeling，ASL）測量的 CBF，代入公式（1）可以求得CMRO2。

基于T2的腦氧代謝率測量

圖1 測量腦氧代謝的原理圖。

Lu等[8]提出了可以用TRUST技術測量全腦的Yv，該技術原理是通過定量上矢狀竇中純血的T2弛豫時間來定量Yv。在TRUST技術中，靜脈回流血液被飽和脈沖標記，而成像面在標記面之下，通過標記和未標記影像的差來獲得靜脈血液信號。其掃描序列由數(shù)據(jù)采集瓣與控制掃描瓣交錯使用，每次采集包含四個范圍從0ms到160ms的有效的回波時間（time of echo，TE）。每次掃描前，要先加入一個預飽和脈沖來抑制靜態(tài)組織的信號，從而增加圖像的信噪比。在圖像數(shù)據(jù)采集之前，預先打入180°的反轉脈沖，將所有的縱向磁化反轉，使得T1弛豫是從－1開始，增強了縱向磁化信號。預先對成像區(qū)進行一次未標記血成像，成為控制像。隨后打入控制脈沖來磁化即將流入的靜脈血，然后需要等待一個恰當?shù)臅r間間隔（約1.2s），使標記后的血液完全流到標記層面并采集圖像成為標記像。把標記像與控制像相減得到灌注像，灌注像就只含有靜脈竇中的靜脈血的信號。在數(shù)據(jù)采集之前，先加入一系列非選擇層面的T2預飽和脈沖，最小限度地減少血液流動帶來的影響。接著施加非選擇性的T2預制脈沖序列，來獲得T2加權信號，其持續(xù)時間由有效回波時間（eTE）表示。使用不同數(shù)量且對血液流動不敏感的T2制備脈沖來調制不同的T2加權信號。得到了靜脈血液的T2值之后，再使用血液樣本擬合Yv、紅細胞比容（hematocrit，Hct）以及T2的三維矯正曲線。通過缺氧控制血氧值，與血氧儀測得的黃金標準值對比。樣本用溫水浴控制并穩(wěn)定溫度，氧氣閥控制血液含氧量，當血液到達穩(wěn)定狀態(tài)時，進行TRUST測試，得到T2預飽和值，TE值等相關方程的MR信號，并引出血液T2值。接著使用雙室交換模型[9]，建立有關Yv，Hct，T2值的分析方程，其中Hct的值是已知的。關系方程如下：

經(jīng)過校正擬合實驗推算出各個參數(shù)。T2制備時間與血流信號的衰減時間擬合然后得出了靜脈血的T2值，最后根據(jù)T2值和靜脈血氧飽和度曲線擬合，最終獲得靜脈中的Yv值[10]。得到了全腦的Yv，還需要利用PC技術測量全腦的CBF[3，11]。在PC技術中，需要測量頸內動脈和椎動脈共四根血管的血流值，相加得到總的供給大腦的CBF。最終利用公式（1）計算出全腦的CMRO2。

Xu等[3]招募36名22～62歲的志愿者利用TRUST技術和PC技術測算全腦的CMRO2，得出的結果為（132.1±20.0）mol／100g／min，與以往文獻記載的用PET技術做的結果有很強的一致性。Xu等[12]用TRUST等技術研究了血液中二氧化碳水平對大腦活動的影響，發(fā)現(xiàn)血碳酸過多會減少大腦的代謝活動，導致大腦自發(fā)連通性降低，腦電圖中低頻率功率譜相對增加，最終得出二氧化碳會對腦活動起到抑制作用的結論。Xu等[13]用TRUST技術和ASL技術做了10名22～35歲的健康年輕人攝入200mg咖啡因后CBF、CMRO2和氧攝取分數(shù)（oxygen extraction fraction，OEF）隨著時間變化的實驗，發(fā)現(xiàn)CBF減少，OEF顯著增加，CMRO2沒有顯著變化。同時發(fā)現(xiàn)不同腦區(qū)CBF變化不一，他做了兩種解釋，第一，不同腦區(qū)對咖啡因的依賴性可能不同；第二，咖啡因對靜脈血管的影響可能是區(qū)域性的。

基于相位的腦氧代謝率測量

測量局部腦氧代謝率首先要測量出局部腦區(qū)的Yv?；谙辔坏拇殴舱窦夹g根據(jù)靜脈血與周圍組織的磁敏感值的不同來成像（圖2），最后可以計算出局部腦區(qū)的Yv。然后再利用ASL技術測算出局部腦區(qū)的CBF，利用菲克原理可以測算出局部腦區(qū)的CMRO2。不同物質在磁場內的表現(xiàn)，可分為順磁性逆磁性及鐵磁性物質。相位圖對順磁性物質十分敏感，如脫氧血紅蛋白、含鐵血黃素、鈣沉積等[14]，這些物質使得在相位圖中顯示出很明顯的低信號。

含氧血紅蛋白在磁場中的順磁性較弱，然而當含氧血紅蛋白經(jīng)過組織后變?yōu)槊撗跹t蛋白時，脫氧血紅蛋白中二價鐵含量升高使得脫氧血紅蛋白在磁場中具有較高的順磁性。因此，可以采用順磁性的不同來反映靜脈中脫氧血紅蛋白的含量。人腦在高強度磁場條件下，相位圖可以很好地記錄下局部腦組織順磁性物質的分布情況，這對人們分析順磁性物質在人腦分布情況是非常有用的。由于脫氧血紅蛋白在動脈中含量極少而主要存在于靜脈血管和組織微靜脈中，所以采集到的圖像主要包含了靜脈與組織微靜脈的信息。在腦組織中，運用該技術可以獲取這些信息，也就可以獲得腦組織中的脫氧血紅蛋白的含量。

Fernandez－Seara等[15]和 Haacke等[16]用靜脈與組織之間的磁敏感的差異測量了頸靜脈的Yv。由于脫氧血紅蛋白的順磁性，靜脈的磁敏感值與周圍組織的磁敏感值會有差異，用Δχvein－tissue表示。在相同的外加磁場下，它們產生的局部磁場也會不同，差異用 ΔΒvein－tissue表 示，由拉莫方程可得：

其中，Δω是靜脈與組織之間質子的角頻率差值，γ為回轉磁率。TE時間下的靜脈與組織之間質子的相位差Δφvein－tissue為：

由此可以看出，主磁場內兩種物質間的磁化率差異引起局部磁場的差異，進而產生質子角頻率的偏移，最終導致相位的差異。

最后根據(jù)ΔΒvein－tissue和 Δχvein－tissue的近似關系計算得出

然后利用Δχvein－tissue與 Yv的線性關系[16]：

計算出更精確的Yv值，公式如下：

其中，Δχdo為動靜脈之間的磁敏感差異值；Β0為主磁場場強；γ＝2.67×108rad／s·T。

Fan等[7]提出了采用多TE回波GRE序列采集相位圖，最后定量Yv，這樣可以減少誤差。以上都是在假設靜脈血管方向與主磁場方向平行，且靜脈血管無限長的情況下得出的Yv的值。

獲得局部腦區(qū)的Yv之后，可以采用ASL技術測算相應腦區(qū)的CBF。ASL的基本原理是對成像平面的上游動脈血液中的水分子進行標記后使其自旋狀態(tài)發(fā)生改變，改變后的血液整體磁化強度就會隨之發(fā)生改變，待磁化強度改變的血液對組織灌注后進行標記成像[17]。ASL技術標記的血液是動脈血，成像面在標記面的上面，標記的血液流向成像平面后對組織進行灌注，灌注像是標記像與控制像的差值，就只有動脈血的信號，從中可以計算出各腦區(qū)的CBF。把測得的Yv與CBF代入公式（1）中，從而計算出局部腦區(qū)的CMRO2。

圖2 不同物質在MRI中的絕對磁敏感值[7]。

Shen等[18]用六只大鼠做了一個創(chuàng)傷性腦損傷對CBF的影響，他用基于相位的磁共振技術檢測大鼠的Yv的變化，該方法新穎、安全、無創(chuàng)，為實驗起了重要的作用。Jain等[19]用新穎的磁共振技術研究了輕度血碳酸過多癥對腦代謝的影響，他用基于相位的磁共振技術監(jiān)測了Yv的變化，用PC技術測量整個大腦的CBF，得出的結果與以往文獻記載的結果相一致，同時該測量方案提高了測量精度和可靠性。Fan等[20]通過讓人體吸入二氧化碳，探索該技術能否檢測到由血碳酸過多引起的OEF的變化。結果表明用該技術測得的OEF的變化與其之前用偽連續(xù)式動脈自旋轉標記技術（pseudo continuous arterial spin labeling，PCASL）技術測的結果有顯著的相關性，這一發(fā)現(xiàn)表明用該技術對測量OEF具有一定的前景。

討論

TRUST技術通過定量上矢狀竇中純血的T2弛豫時間來定量計算Yv，它與PC技術或者ASL技術相結合可以非侵入地測量全腦的CMRO2。該技術可以檢測到整個大腦生理狀態(tài)的改變，如老化、睡眠、高碳酸血癥和氧過多等癥，以及了解靜息狀態(tài)的大腦活動。由于TRUST采用回波平面成像（echo－planar imaging，EPI）序列來采集，空間分辨率低，無法分辨腦中小靜脈的T2信號，區(qū)域化測量CMRO2難以實現(xiàn)。

基于相位的磁共振方法不僅可以測量區(qū)域腦區(qū)的Yv，也可以根據(jù)上矢狀竇與周圍腦組織的的磁敏感差異定量全腦水平的Yv[21]，該技術與PC或者ASL技術相結合可以測量全腦水平的CMRO2?；谙辔坏拇殴舱穹椒y量Yv與其他磁共振方法相比，它提供了一個更高時間分辨力[22]；第二，靜脈血管周圍組織的磁化率接近于水[23]，因而無需校準。

TRUST技術與基于相位的磁共振技術相比，TRUST技術測量Yv的方法更為直接[24]。TRUST技術是通過T2加權測得標記的靜脈血與未標記的靜脈血之間的圖像信號之差直接測量Yv的信息，測得的結果穩(wěn)定性較高。該技術使用自旋標記原理從周圍組織中分離出血信號，因此，避免了手動選擇興趣區(qū)（regions of interest，ROI）區(qū)域，或者要使圖像的體素比血管管腔小等問題；而基于相位的磁共振技術是利用靜脈與周圍組織間的磁敏感差異形成圖像對比來解算Yv的信息。該技術測得的Yv的結果與靜脈血管走向及靜脈擴張程度有復雜的關系，雖然有學者[21]假設血管為無限長圓柱并與主磁場平行來簡化該技術對Yv的測量，但是這種情況會對實際血管的分析帶來影響，加之在計算過程中對靜脈和組織間磁感應強度和磁敏感之間的關系進行了近似，所以采用基于相位的磁共振技術測量Yv會有誤差。Fan等[7]通過測量大腦血管和周圍腦軟組織之間的磁化率差異，依據(jù)CMRO2和CBF之間的線性關系近似計算得出了區(qū)域CMRO2，但這種近似具有很強的假設性[21]，所以用基于相位的磁共振技術測得的CMRO2的值不如用TRUST技術測的精確。

總結

基于T2與相位的定量CMRO2的磁共振技術是近年來新提出來的，它具有非侵襲、安全可靠、成本低等優(yōu)點。以后隨著高場強設備的引入，用TRUST技術測量小血管的Yv將會成為可能，隨著TRUST技術理論進一步完善，以后會縮短它的掃描時間，更好的應用于臨床與研究；基于相位的磁共振技術分辨力會進一步提高，以后隨著理論進一步完善，比如把血管長度、傾斜度等限制問題解決，它的測量精度會進一步提高，以后會更好地應用于臨床診斷、鑒別診斷及科學研究之中。

[1]付占立，張春麗，王榮福.PET局部腦氧代謝率測定[J].同位素，2006，19（3）：184－189.

[2]Ito H，Ibaraki M，Kanno I，et al.Changes in cerebral blood flow and cerebral oxygen metabolism during neural activation measured by positron emission tomography：comparisonwith blood oxygenation level－dependent contrast measured by functional magnetic resonance imaging[J].J Cereb Blood Flow Metab，2005，25（3）：371－377.

[3]Xu F，Ge Y，Lu H.Noninvasive quantification of whole－brain cerebral metabolic rate of oxygen（CMRO2）by MRI[J].Magn Reson Med，2009，62（1）：141－148.

[4]Lisa C.Krishnamurthy，Peiying Liu，Yulin Ge，et al.Vessel－specific quantification of blood oxygenation with T2－relaxation－underphase－contrast MRI[J].Magn Reson Med，2013，71（3）：978－989.

[5]Varsha Jain，Jeremy Magland，Michael Langham，et al.High temporal resolution in vivo blood oximetry via projection－based T2measurement[J].Magn Reson Med，2014，70（3）：785－790.

[6]Jain V，Langham MC，Wehrli FW.MRI estimation of global brain oxygen consumption rate[J].J Cereb Blood Flow Metab，2010，30（9）：1598－1607.

[7]Fan AP.Phase－based regional oxygen metabolism in magnetice imaging at high field[D].Palo Alto：Stanford Unversity，2010.

[8]Lu H，Ge Y.Quantitative evaluation of oxygenation in venous vessels using T2－relaxation－under－spin－tagging MRI[J].Magn Reson Med，2008，60（2）：357－363.

[9]Barth M，Moser E.Proton NMR relaxation times of human blood samples at 1.5Tand implications for functional MRI[J].Cell Mol Biol，1997，43（5）：783－791.

[10]Lu H，Xu F，Grgac K，et al.Calibration and validation of TRUST MRI for the estimation of cerebral blood oxygenation[J].Magn Reson Med，2012，67（1）：42－49.

[11]Haacke EM，Brown RW，Thompson MR，et al.Magnetic resonance imaging：physical principles and sequence design[M].New York：Wiley－Liss，1999.

[12]Xu F，Jinsoo Uh，Brier M.The influence of carbon dioxide on brain activity and metabolism in conscious humans[J].J Cereb Blood Flow Metab，2011，31（1）：58－67.

[13]Xu F，Liu P，Pekar JJ，et al.Does acute caffeine ingestion alter brain metabolism in young adults[J].Neuroimage，2015，110（1）：39－47.

[14]Sehgal V，Delproposto Z，Haacke EM，et al.Clinical applications of neuroimaging with susceptibility－weighted imaging[J].J Magn Reson Imaging，2005，22（4）：439－450.

[15]Fernandez－Seara MA，Techawiboonwong A，Detre JA，et al.MR susceptometry for measuring global brain oxygen extraction[J].Magn Reson Med，2006，55（5）：967－973.

[16]Haacke EM，Lai S，Reichenbach JR，et al.In vivo measurement of blood oxygen saturation using magnetic resonance imaging：a direct validation of the blood oxygen level dependent concept in functional brain imaging[J].Hum Brain Mapp，1997，5（5）：341－346.

[17]Golay X，Hendrikse J，Lim TC.Perfusion imaging using arterial spin labeling[J].Top Magn Reson Imaging，2004，15（1）：10－27.

[18]Shena Y，Kou Z，Kreipke CW，et al.In vivo measurement of tissue damage，oxygen saturation changes and blood flow changes after experimental traumatic brain injury in rats using susceptibility weighted imaging[J].Magn Reson Imaging，2007，25（2）：219－227.

[19]Jain V，Langham MC，F(xiàn)loyd TF，et al.Rapid magnetic resonance measurement of global cerebral metabolic rate of oxygen consumption in humans during rest and hypercapnia[J].J Cereb Blood Flow Metab，2011，31（7）：1504－1512.

[20]Fan AP，Evans KC，Stout JN，et al.Regional quantification of cerebral venous oxygenation from MRI susceptibility during hypercapnia[J].Neuroimage，2015，104（1）：146－155.

[21]Fan AP，Benner T，Bolar DS，et al.Phase－based regional oxygen metabolism （PROM）using MRI[J].Magn Reson Med，2012，67（3）：669－678.

[22]Langham MC，F(xiàn)loyd TF，Mohler ER III，et al.Evaluation of cuffinduced ischemia in the lower extremity by magnetic resonance oximetry[J].J Am Coll Cardiol，2010，55（6）：598－606.

[23]Schenck JF.The role of magnetic susceptibility in magnetic resonance imaging：MRI magnetic compatibility of the first and second kinds[J].Med Phys，1996，23（6）：815－850.

[24]Ge Y，Zhang Z，Lu H，et al.Characterizing brain oxygen metabolism in patients with multiple sclerosis with T2relaxation under spin tagging MRI[J].J Cereb Blood Flow Metab，2012，32（3）：403－412.