磁気共鳴画像法と運動補償によるヒト胎児血流定量化
Please note that all translations are automatically generated.
Click here for the English version.
概要
ここでは、MRIで胎児の血流を迅速に測定し、動き補正と心同期を後方視的に行うためのプロトコルを紹介します。
Abstract
磁気共鳴画像法(MRI)は、心血管形態と心機能の臨床評価のための重要なツールです。また、位相差MRIに基づく血流定量化の標準治療として認められています。このような血流の測定は何十年にもわたって成人で可能でしたが、この能力を胎児の血流に拡張する方法は最近開発されたばかりです。
主要血管における胎児血流の定量化は、先天性心疾患(CHD)や胎児発育制限(FGR)などの胎児の病状を監視するために重要です。CHDは、胎児の血液の流れを変える心臓構造と血管系の変化を引き起こします。FGRでは、血流の経路はシャントの拡張によって変化し、脳への酸素化された血液供給が増加します。血流の定量化により、胎児の病状の重症度を評価することができ、子宮内患者の管理や産後ケアの計画に適しています。
位相差MRIをヒト胎児に適用する際の主な課題には、血管サイズが小さく、胎児の心拍数が高い、母体呼吸によるMRIデータの破損の可能性、予測不可能な胎児の動き、データ収集を同期させるための従来の心臓同期法の欠如などがあります。ここでは、位相差MRIを用いた胎児血流の定量化を可能にした当研究室の最近の技術開発について、加速イメージング、運動補償、心同期の進歩などについて述べる。
Introduction
胎児発育制限(FGR)や先天性心疾患(CHD)などの胎児の病状を監視するには、胎児循環の包括的な評価が必要です1,2,3。子宮内では、出生後ケアの患者管理と計画は、胎児の病状の重症度に依存します4,5,6,7。MRIによる胎児血流定量の実現可能性と胎児の病状の評価におけるその応用が最近実証されました3,8,9。しかし、イメージング法は、高い時空間分解能を達成するためのイメージング時間の増加、心臓同期法の欠如、予測不可能な胎児の動きなどの課題に直面しています10。
胎児血管系は小さな構造(下行大動脈、動脈管、上行大動脈、主肺動脈、上大静脈を構成する主要血管の直径~5mm)で構成されています11,12,13)。これらの構造を分解し、流れを定量化するためには、高い空間分解能でのイメージングが必要です。さらに、胎児の心拍数は成人の約2倍です。したがって、胎児の心周期全体の動的な心臓運動と血流を解決するには、高い時間分解能も必要です。この高い時空間分解能での従来のイメージングには、比較的長い取得時間が必要です。この問題に対処するために、加速胎児MRI14,15,16が導入されました。簡単に言うと、これらの加速技術には、データ収集中の周波数領域でのアンダーサンプリングと、反復技術を使用した遡及的な高忠実度再構成が含まれます。そのようなアプローチの1つは、圧縮センシング(CS)再構成であり、再構成された画像が既知のドメインでまばらであり、アンダーサンプリングアーティファクトがインコヒーレントである場合に、高度にアンダーサンプリングされたデータから画像を再構成することを可能にする17。
胎児イメージングにおける動きは大きな課題を提示します。運動の腐敗は、母親の呼吸運動、母親のバルク運動、または胎児の肉眼的運動から発生する可能性があります。母親の呼吸は胎児の周期的な翻訳につながりますが、胎児の動きはより複雑です。胎児の動きは、局所的または肉眼的10,18として分類できます。局所的な動きには、体のセグメントのみの動きが含まれます。それらは通常約10〜14秒間続き、その頻度は妊娠とともに増加します(正期産期では1時間あたり~90)10。これらの動きは一般に小さな破損を引き起こし、関心のあるイメージング領域には影響しません。ただし、胎児の全体的な動きは、平面運動コンポーネントを介した深刻な画像の破損につながる可能性があります。これらの動きは、脊椎によって媒介される全身の動きであり、60〜90秒間続きます。
胎児の動きによるアーティファクトを避けるために、最初に母親の動きを最小限に抑えるための措置が取られます。妊娠中の女性は、スキャナーベッドの支持枕を使用してよりリラックスし、快適なガウンを着て、閉所恐怖症を軽減するためにスキャナーの横にパートナーがいる場合があります19,20。母親の呼吸運動の影響を軽減するために、研究では母親の息止め下で胎児MR検査を実施しました21,22,23。ただし、そのような取得は、妊娠中の被験者の息止め耐性の低下を考えると、短く(~15秒)でなければなりません。.最近、胎児MRI14,15,16に対して遡及的動き補正法が導入されている。これらの方法は、登録ツールキットを使用して胎児の動きを追跡し、動きを補正するか、取得したデータの修正不可能な部分を破棄します。
最後に、出生後の心臓MR画像は、従来、心電図(ECG)ゲーティングを使用して取得され、データ収集を心周期に同期させます。ゲーティングなしでは、心周期全体からの心臓の動きと拍動の流れが組み合わされ、アーティファクトが生成されます。残念なことに、胎児心電図信号は、母体心電図信号24からの干渉および磁場25からの歪みに苦しんでいる。したがって、セルフゲーティング、メトリック最適化ゲーティング(MOG)およびドップラー超音波ゲーティング21,26,27,28を含む、胎児心臓ゲーティングに対する代替の非侵襲的アプローチが提案されている。
次のセクションで説明するように、胎児の血流を定量化する私たちのMRIアプローチは、私たちの研究室で開発され、加速されたMRI取得の動き補正と反復再構成を組み合わせた新しいゲーティング法MOGを活用しています。このアプローチは、以前に発表された研究14のパイプラインに基づいており、次の5つの段階で構成されています:(1)胎児血流取得、(2)リアルタイム再構成、(3)運動補正、(4)心臓ゲーティング、および(5)ゲート再構成。
Protocol
すべてのMRIスキャンは、私たちの機関研究倫理委員会によって承認された研究の一環として、ボランティアからのインフォームドコンセントで実施されました。 注意: 以下に説明する方法は、3T MRIシステムで使用されています。取得は、放射状位相差MRIシーケンスを用いて行われる。このシーケンスは、メーカーのデカルト位相差MRIの読み出し軌道を変更することによって(星状パターンを達成するために)準備されました。シーケンスおよびサンプルプロトコルは、当社のC2P交換プラットフォームを通じてリクエストに応じて入手できます。この作業のすべての再構築は、32 GBメモリ、8コアの3.40 GHzプロセッサ、および1024コンピューティング統合デバイスアーキテクチャ(CUDA)コアを備えた2GBグラフィックカードを備えた標準のデスクトップコンピューターで実行されました。画像再構成はMATLABで実施した。不一様高速フーリエ変換(NUFFT)29は、グラフィックスプロセッシングユニット(GPU)で実行されました。動き補正パラメータは、エラスティックス30を用いて計算した。図1は、プロトコルを時系列で示しており、取得した速度エンコード(図1では色分け)が再構成の各段階で代表的な画像でどのように処理されるかを追跡しています。再構築コードは https://github.com/datta-g/Fetal_PC_MRI で入手できます。ここではプロトコルの手順を提供しますが、これらのアルゴリズムの手順のほとんどはパイプラインで自動化されています。 1. 科目のポジショニングとローカライザー試験 MRI検査のために、母親がMRI台上で好みの快適な位置、通常は仰臥位または側褥瘡位に身を置くのを支援します。 心臓コイルを母親の腹部の上に置きます。 MRIテーブルをマグネットボアにロードし、スキャンが開始されようとしていることを母親に通知します。 ローカライザー検査を実行して、胎児の体を特定します(解像度:0.9 x 0.9 x 10 mm3、TE / TR:5.0 / 15.0ミリ秒、FOV:450 x 450 mm2、スライス:6)。 洗練されたローカライザー検査を実行して、胎児の心臓を中心としたスライスグループで胎児血管系を見つけます(解像度1.1 x 1.1 x 6.0 mm3、TE / TR:2.69 / 1335.4ミリ秒、FOV:350 x 350 mm2、スライス:10、向き:胎児に対して軸方向)。 胎児の血管をより明確に見るために、矢状および冠状の向きで洗練されたローカライザーを繰り返します。 全体的な胎児の動きの場合には洗練されたローカライザーを繰り返します。 2. 胎児血流データの取得 ローカライザー検査を使用して胎児血管を見つけます。例えば、下行大動脈は矢状面の脊椎近くの長いまっすぐな血管です。上行大動脈と主肺動脈は、それぞれ左心室と右心室を離れる血管として識別できます。動脈管は、下行大動脈の近位にある主肺動脈の下流セグメントとして追跡できます。上大静脈は、胎児の心臓の基部近くの軸面から、上行大動脈に隣接する血管として識別できます。 対象の胎児血管の軸に垂直なスライスを処方します。MRIコンソールコンピュータ上のスライスガイドラインを回転および移動して、ターゲット血管と垂直に交差させます。 スキャンパラメータを設定します(取得タイプ:ラジアル位相差MRI、分解能:1.3 x 1.3 x 5.0 mm 3、エコー時間(TE)/繰り返し時間(TR):3.25 / 5.75 ms、視野(FOV):240 x 240 mm2、スライス:1、速度エンコーディング:対象容器に応じて100-150 cm / s、速度エンコーディング方向:スループレーン、ラジアルビュー:エンコードあたり1500、フリーブリージング)。 スキャンを実行し、実行され、MRIコンソールコンピューターに表示された最初の時間平均再構成に基づいて処方を確認します。ローカライザースキャンと位相差スキャンを繰り返します ターゲット容器が存在しないか、最初の再構成から識別できない場合。集録された生データは、図1Aの回路図に、速度補償された平面取得とスループレーンの取得がそれぞれ赤と青で色分けされて表されます。 対象血管ごとに胎児血流データ取得を繰り返す。注:取得した生データ(形式:DATファイル)は、オフラインで再構築するために転送する必要があります。たとえば、シーメンスのスキャナーでは、「twix」を実行することでこれを実行できます。取得した生データをリストから右クリックし、「合計RAIDファイルをコピー」を選択します。 3.胎児測定値の動き補正 図1Bの回路図で表される空間的全変動(STV、重み:0.008)と時間的全変動(TTV、重み:0.08)の正則化を利用した共役勾配降下最適化を15回反復してCSを使用して、取得したデータからリアルタイム系列(時間分解能:370ミリ秒、放射状ビュー:64)を再構築します。 MATLABで開発されたグラフィックユーザーインターフェイスを使用して、この最初のリアルタイム再構成から関心のある容器を含む関心領域(ROI)を選択します。このステップでは、ユーザは、標的大血管または胎児の心臓などの胎児の解剖学的構造を囲む輪郭を描画しなければならない。 elastix30 を使用して剛体モーショントラッキングを実行します (経験的に最適化されたパラメーターを持つ正規化された相互情報に基づく: 4 つのピラミッド レベル、300 回の反復、並進変換)。 他のすべてのフレームと低い相互情報量(MI)を共有する追跡されたリアルタイムフレームを拒否します(MIは平均MIの四分位範囲の1.5倍未満です)。これらのフレームは、平面運動または総胎児運動によって表されると見なされます。 残りのフレームからの連続リアルタイムフレーム(ギャップなし)の最も長いシリーズに対応するMRIデータを、さらなる再構成に使用される静止期間として使用します。 並進モーション補正パラメータを、リアルタイム系列の時間分解能(370ms)から静止時の集録のTR(5.75ms)まで補間します。 補間されたパラメータをMRIデータの定義された静止期間に適用するには、次のように位相を変調します。ここで、s’ はモーション補正データ、k x と k y は k-空間の座標、s は取得した未補正データ、Δ x と Δ y は空間内の追跡変位、j は を表します。注:この作業の正則化係数のすべての数値は、以前の実験で最適化されました。これは、ブルートフォースグリッド検索を使用して、高度にサンプリングされた胎児参照データセットの再構成と同じデータセットからの遡及的にアンダーサンプリングされた症例との間の誤差を最小限に抑える正則化係数を見つけることで達成されました。 4.胎児の心拍数を解く CSを使用して取得したデータを使用して、より高い時間分解能(時間分解能:46ミリ秒、放射状ビュー:8)で2番目のリアルタイム画像シリーズを再構築し、図1Cの回路図で表されるように、STV(重み:0.008)およびTTV(重み:0.08)正則化を使用した共役勾配降下最適化を15回繰り返します。 目的の胎児血管を含むROIを再選択します。 リアルタイム系列でマルチパラメータMOGを実行して、時間依存の胎児心拍数を導き出します。 ビンモーションは、導出された心拍数波形を使用してMRIデータを15の心相に補正しました。このステップでは、前のステップの心拍数を使用して心相の時間的境界が計算されます。たとえば、k番目のハートビートのi番目のフェーズの境界は、次のように与えられます。ここで、HR(K) は k番目のハートビートが発生する時刻です。n番目のラジアル取得のタイムスタンプは、(n x TR)で与えられます。心臓相の境界内にあるタイムスタンプを持つデータは、その段階に割り当てられます。注:MOGは、関心領域にわたって画像メトリックを最適化するCINE画像を作成するために、マルチパラメータ胎児心拍数モデルに基づいて取得されたデータの反復ビニングを含むゲーティング技術26である。 5. 胎児CINEの再建 STV(重み:0.025)およびTTV(重量:0.01)正則化による共役勾配降下最適化の10回の反復で、ビンモーション補正MRIデータとCSを使用して胎児流CINEを再構築します。このステップでは、図1Dの回路図に示すように、フロー補償集録CFC用とフローエンコードデータCFE用の2つのCINEが生成されます。 C FEの元素単位積とCFCの複素共役の位相によって与えられる速度像を計算します。 バックグラウンド位相補正31を適用して、渦電流効果を補正します。簡単に言えば、この自動ステップでは、平面が静的な胎児および母体組織の相にフィットする。補正は、4.2で計算された速度に敏感な位相から平面を引くことによって実行されます。 再構築されたデータを DICOM ファイルに書き込みます。 DICOMをセグメントv2.232などのフロー解析ソフトウェアにロードします。 解剖学的および速度に敏感な画像を使用して、対象の血管の内腔を含むROIを描画します。 ROIをすべての心臓相に伝播し、血管の直径の変化を補正します。 流量測定値を記録します。
Representative Results
一般に、フローのフェーズMRI検査は、下行大動脈、上行大動脈、主肺動脈、動脈管、上大静脈、および臍帯静脈の6つの主要な胎児血管を対象としています。これらの血管は、CHDおよびFGRに関与していることが多く、胎児全体の血液分布に影響を与えるため、臨床医にとって興味深いものです9。放射状位相差MRIの一般的なスキャン時間は血管あたり17秒であるため、スキャン…
Discussion
この方法は、ヒト胎児大血管の血流の非侵襲的測定を可能にし、反復再構成技術を利用することにより、遡及的運動補正および心ゲーティングを可能にする。胎児の血流の定量化は、過去にMRIで行われています1,3,8,9。これらの研究では、スキャナーでの最初の再構成から…
開示
The authors have nothing to disclose.
Acknowledgements
何一つ。
Materials
| elastix | Image Sciences Institute, University Medical Center Utrecht | Image registration software | |
| Geforce GTX 960 | Nvidia | 04G-P4-3967-KR | |
| gpuNUFFT | CAI²R | Non-uniform fast Fourier transform | |
| MAGNETOM Prisma | Siemens | 10849583 | |
| MATLAB | MathWorks | ||
| Radial Phase Contrast MRI sequence | Trajectory modification of manufacturer's Cartesian Phase Contrast sequence | ||
| Segment | Medvisio | Data analysis | |
| VENGEANCE | Corsair | LPX DDR4-2666 |
参考文献
- Zhu, M. Y., et al. The hemodynamics of late-onset intrauterine growth restriction by MRI. American Journal of Obstetrics and Gynecology. 214 (3), 1-17 (2016).
- Zhu, M. Y., Jaeggi, E., Roy, C. W., Macgowan, C. K., Seed, M. Reduced combined ventricular output and increased oxygen extraction fraction in a fetus with complete heart block demonstrated by MRI. HeartRhythm Case Reports. 2 (2), 164-168 (2016).
- Sun, L., et al. Reduced Fetal Cerebral Oxygen Consumption is Associated With Smaller Brain Size in Fetuses With Congenital Heart Disease. Circulation. 131 (15), 1313-1323 (2015).
- Freud, L. R., et al. Fetal aortic valvuloplasty for evolving hypoplastic left heart syndrome: postnatal outcomes of the first 100 patients. Circulation. 130 (8), 638-645 (2014).
- Peleg, D., Kennedy, C. M., Hunter, S. K. Intrauterine growth restriction: identification and management. American Family Physician. 58 (2), 453-467 (1998).
- Krishna, U., Bhalerao, S. Placental Insufficiency and Fetal Growth Restriction. Journal of Obstetrics and Gynaecology of India. 61 (5), 505-511 (2011).
- Seravalli, V., Miller, J. L., Block-Abraham, D., Baschat, A. A. Ductus venosus Doppler in the assessment of fetal cardiovascular health: an updated practical approach. Acta Obstetricia et Gynecologica Scandinavica. 95 (6), 635-644 (2016).
- Seed, M., et al. Feasibility of quantification of the distribution of blood flow in the normal human fetal circulation using CMR: a cross-sectional study. Journal of Cardiovascular Magnetic Resonance. 14 (1), 79 (2012).
- Prsa, M., et al. Reference ranges of blood flow in the major vessels of the normal human fetal circulation at term by phase-contrast magnetic resonance imaging. Circulation. Cardiovascular Imaging. 7 (4), 663-670 (2014).
- Piontelli, A. . Development of Normal Fetal Movements: The Last 15 Weeks of Gestation. , (2015).
- Cartier, M., et al. The normal diameter of the fetal aorta and pulmonary artery: echocardiographic evaluation in utero. American Journal of Roentgenology. 149 (5), 1003-1007 (1987).
- Ruano, R., de Fátima Yukie Maeda, M., Niigaki, J. I., Zugaib, M. Pulmonary artery diameters in healthy fetuses from 19 to 40 weeks’ gestation. Journal of Ultrasound in Medicine. 26 (3), 309-316 (2007).
- Nowak, D., Kozłowska, H., Żurada, A., Gielecki, J. Diameter of the ductus arteriosus as a predictor of patent ductus arteriosus (PDA). Central European Journal of Medicine. 6 (4), 418-424 (2011).
- Goolaub, D. S., et al. Multidimensional fetal flow imaging with cardiovascular magnetic resonance: a feasibility study. Journal of Cardiovascular Magnetic Resonance. 20 (1), 77 (2018).
- Roy, C. W., Seed, M., Kingdom, J. C., Macgowan, C. K. Motion compensated cine CMR of the fetal heart using radial undersampling and compressed sensing. Journal of Cardiovascular Magnetic Resonance. 19 (1), 29 (2017).
- van Amerom, J. F. P., et al. Fetal cardiac cine imaging using highly accelerated dynamic MRI with retrospective motion correction and outlier rejection. Magnetic Resonance in Medicine. 79 (1), 327-338 (2018).
- Lustig, M., Donoho, D., Pauly, J. M. Sparse MRI: The application of compressed sensing for rapid MR imaging. Magnetic Resonance in Medicine. 58 (6), 1182-1195 (2007).
- Edwards, D. D., Edwards, J. S. Fetal movement: development and time course. Science. 169 (3940), 95-97 (1970).
- Malamateniou, C., et al. Motion-Compensation Techniques in Neonatal and Fetal MR Imaging. American Journal of Neuroradiology. 34 (6), 1124-1136 (2013).
- Rutherford, M., et al. MR imaging methods for assessing fetal brain development. Developmental Neurobiology. 68 (6), 700-711 (2008).
- Haris, K., et al. Self-gated fetal cardiac MRI with tiny golden angle iGRASP: A feasibility study: Self-Gated Fetal Cardiac MRI with iGRASP. Journal of Magnetic Resonance Imaging. 46 (1), 207-217 (2017).
- Glenn, O. A. MR imaging of the fetal brain. Pediatric Radiology. 40 (1), 68-81 (2010).
- Rodríguez-Soto, A. E., et al. MRI Quantification of Human Fetal O2 Delivery Rate in the Second and Third Trimesters of Pregnancy. Magnetic Resonance in Medicine. 80 (3), 1148-1157 (2018).
- Sameni, R., Clifford, G. D. A Review of Fetal ECG Signal Processing; Issues and Promising Directions. The Open Pacing, Electrophysiology & Therapy Journal. 3, 4-20 (2010).
- Millis, R. Advances in Electrocardiograms: Methods and Analysis. BoD – Books on Demand. , (2012).
- Jansz, M. S., et al. Metric optimized gating for fetal cardiac MRI. Magnetic Resonance in Medicine. 64 (5), 1304-1314 (2010).
- Yamamura, J., et al. Cardiac MRI of the fetal heart using a novel triggering method: initial results in an animal model. Journal of Magnetic Resonance Imaging: JMRI. 35 (5), 1071-1076 (2012).
- Larson, A. C., et al. Self-gated cardiac cine MRI. Magnetic Resonance in Medicine. 51 (1), 93-102 (2004).
- Knoll, F., Schwarzl, A., Diwoky, C., Sodickson, D. K. gpuNUFFT-An open source GPU library for 3D regridding with direct Matlab interface. Proceedings of the 22nd Annual Meeting of ISMRM. , (2014).
- Klein, S., Staring, M., Murphy, K., Viergever, M. A., Pluim, J. P. W. elastix: a toolbox for intensity-based medical image registration. IEEE Transactions on Medical Imaging. 29 (1), 196-205 (2010).
- Walker, P. G., et al. Semiautomated method for noise reduction and background phase error correction in MR phase velocity data. Journal of Magnetic Resonance Imaging. 3 (3), 521-530 (1993).
- Heiberg, E., et al. Design and validation of Segment – freely available software for cardiovascular image analysis. BMC Medical Imaging. 10 (1), 1 (2010).
- Inder, T. E., Volpe, J. J. Chapter 17 – Intrauterine, Intrapartum Assessments in the Term Infant. Volpe’s Neurology of the Newborn (Sixth Edition). , 458-483 (2018).
- Pelc, N. J., Herfkens, R. J., Shimakawa, A., Enzmann, D. R. Phase contrast cine magnetic resonance imaging. Magnetic Resonance Quarterly. 7 (4), 229-254 (1991).
- Steeden, J. A., Atkinson, D., Hansen, M. S., Taylor, A. M., Muthurangu, V. Rapid flow assessment of congenital heart disease with high-spatiotemporal-resolution gated spiral phase-contrast MR imaging. Radiology. 260 (1), 79-87 (2011).
- Kowalik, G. T., Knight, D., Steeden, J. A., Muthurangu, V. Perturbed spiral real-time phase-contrast MR with compressive sensing reconstruction for assessment of flow in children. Magnetic Resonance in Medicine. 83 (6), 2077-2091 (2020).
