Refine
Has Fulltext
- yes (3)
Is part of the Bibliography
- yes (3)
Year of publication
- 2018 (3)
Document Type
- Journal article (2)
- Doctoral Thesis (1)
Keywords
- MOLLI (3) (remove)
Institute
- Physikalisches Institut (3) (remove)
Sonstige beteiligte Institutionen
Quantitative nuclear magnetic resonance imaging (MRI) shifts more and more into the focus of clinical research. Especially determination of relaxation times without/and with contrast agents becomes the foundation of tissue characterization, e.g. in cardiac MRI for myocardial fibrosis. Techniques which assess longitudinal relaxation times rely on repetitive application of readout modules, which are interrupted by free relaxation periods, e.g. the Modified Look-Locker Inversion Recovery = MOLLI sequence. These discontinuous sequences reveal an apparent relaxation time, and, by techniques extrapolated from continuous readout sequences, a putative real T1 is determined. What is missing is a rigorous analysis of the dependence of the apparent relaxation time on its real partner, readout sequence parameters and biological parameters as heart rate. This is provided in this paper for the discontinuous balanced steady state free precession (bSSFP) and spoiled gradient echo readouts. It turns out that the apparent longitudinal relaxation rate is the time average of the relaxation rates during the readout module, and free relaxation period. Knowing the heart rate our results vice versa allow to determine the real T1 from its measured apparent partner.
Contrast and non-contrast MRI based characterization of myocardium by T1-mapping will be of paramount importance to obtain biomarkers, e.g. fibrosis, which determines the risk of heart failure patients.
T1-mapping by the standard post-processing of the modified look-locker inversion recovery (MOLLI) lacks of accuracy when trying to reduce its duration, which on the other hand, is highly desirable in patients with heart failure. The recently suggested inversion group fitting (IGF) technique, which considers more parameters for fitting, has a superior accuracy for long T1 times despite a shorter duration. However, for short T1 values, the standard method has a superior precision. A conditional fitting routine is proposed which ideally takes advantage of both algorithms.
Materials and methods
All measurements were performed on a 1.5 T clinical scanner (ACHIEVA, Philips Healthcare, The Netherlands) using a MOLLI 5(n)3(n)3 prototype with n(heart beats) being a variable waiting time between inversion experiments. Phantom experiments covered a broad range of T1 times, waiting times and heart rates. A saturation recovery experiment served as a gold standard for T1 measurement. All data were analyzed with the standard MOLLI, the IGF fit and the conditional fitting routine and the obtained T1 values were compared with the gold standard. In vivo measurements were performed in a healthy volunteer and a total of 34 patients with normal findings, dilative cardiomyopathy and amyloidosis.
Results
Theoretical analysis and phantom experiments provided a threshold value for an apparent IGF
determining processing with IGF post processing for values above, or switching to the standard technique for values below. This was validated in phantoms and patients measurements. A reduction of the waiting time to 1 instead of 3 heart beats between the inversion experiments showed reliable results. The acquisition time was reduced from 17 to 13 heart beats. The in vivo measurements showed ECV values between 25% (18–33%; SD 0.03) in the healthy, 30% (22–40%; SD 0.04) in patients with DCM and 45% (30–60%; SD 0.9) in patients with amyloidosis.
Conclusion
The adopted post-processing algorithm determines long T1 values with high accuracy and short T1 values while maintaining a high precision. Based on reduction of waiting time, and independence of heart rate, it shortens breath hold duration and allows fast T1-mapping, which is frequently a prerequisite in patients with cardiac diseases.
Quantifizierung myokardialer Mikrostruktur und Perfusion mittels longitudinaler NMR Relaxation
(2018)
Ziel der Arbeit war es die Quantifizierung funktioneller bzw. mikrostruktureller Parameter des Herzmuskels mit Hilfe T1-basierter Methoden zu verbessern. Diese Methoden basieren darauf, die gewünschte Information durch eine geeignete Präparation der Magnetisierung bzw. durch die Gabe von Kontrastmittel in den Zeitverlauf der longitudinalen Relaxation zu kodieren. Aus der Änderung der Relaxationszeit läßt sich dann die gewünschte Information bestimmen. Dafür sollte sowohl der Einfluß der Anatomie als auch derjenige der Meßmethodik auf die Bestimmung der longitudinalen Relaxationszeit und damit auf die Quantifizierung der Funktion bzw. Mikrostrukturparameter untersucht werden.
Speziell der Einfluß der Bildgebungssequenz führt dazu, daß nur eine scheinbare Relaxationszeit gemessen wird. Während dies keinen Einfluß auf die T1-basierte Bestimmung der untersuchten Mikrostrukturparameter hatte, ergab sich für die Perfusionsquantifizierung eine deutliche Abhängigkeit von den Parametern der verwendeten IRLL-Sequenz. Um diesen Einfluß gerecht zu werden, wurden an die Meßmethodik angepaßte Gleichungen zur Bestimmung der Perfusion gefunden mit denen die systematischen Abweichungen korrigiert werden können. Zusätzlich reduzieren die angepaßten Gleichungen die Anforderungen bezüglich der Inversionsqualität im schichtselektiven Experiment. Dies wurde in einem weiteren Projekt bei der Bestimmung der Nierenperfusion im Mausmodell ausgenutzt.
Neben der Untersuchung der Auswirkungen der Meßmethode wurde auch der Einfluß der anatomischen Besonderheiten des Blutkreislaufs am Herzen auf die Parameterquantifizierung mittels T1-basierter Methoden untersucht. Es konnte gezeigt werden, daß auf Grund der Anatomie des Herzens bei typischen Orientierungen der Bildgebungsschicht, auch bei der schichtselektiven Inversionspräparation der Magnetisierung des Herzmuskels ein Anteil des Blutpools invertiert wird. Daraus folgt, daß die vereinfachende Annahme, nach welcher bei schichtselektiver Präparation in Folge von Perfusion nur Blut mit Gleichgewichtsmagnetisierung den Herzmuskel erreicht, nicht erfüllt ist. Es konnte gezeigt werden, daß dies bei Perfusion zu einer deutlichen Unterschätzung der berechneten Perfusionswertes führt. Um mit diesem Problem umgehen zu können, wurde aufbauend auf einem vereinfachten Modell der zeitlichen Entwicklung der Blutmagnetisierung eine Korrektur für die Bestimmung der Perfusionswerte gefunden welche den Einfluß der anatomischen Besonderheiten berücksichtigt.
Das für die Perfusionskorrektur eingeführte Model prognostiziert ebenso, daß auch bei schichtselektiver Inversion die T1-basierte Bestimmung der untersuchten Mikrostrukturparameter von der Perfusion abhängig wird und eine systematische Überschätzung der quantifizierten Werte verursacht. Da die Perfusion im Kleintier deutlich höher ist als im Menschen, ist dieser Einfluß besonders in der präklinischen Forschung zu beachten. So können dort allein durch verminderte Perfusion deutliche Änderungen in den bestimmten Werten der Mikrostrukturparameter erzeugt werden, welche zu einer fehlerhaften Interpretation der Ergebnisse führen und somit ein falsches Bild für die Vorgänge im Herzmuskel suggerieren. Dabei bestätigt der Vergleich mit experimentellen Ergebnissen aus der Literatur die Vorhersagen für das Rattenmodell. Beim Menschen ist der prognostizierte Effekt deutlich kleiner. Der prognostizierte Fehler bspw. im RBV-Wert liegt in diesem Fall bei etwa 10% und wird üblicherweise in der aktuellen Forschung vernachlässigt. Inwieweit dies in er klinischen Forschung gerechtfertigt ist, muß in weiteren Untersuchungen geklärt werden.
Den untersuchten Methoden zur Bestimmung von funktionellen und mikrostrukturellen Parametern ist gemein, daß sie eine exakte Quantifizierung der longitudinalen Relaxationszeit T1 benötigen. Dabei ist im Kleintierbereich die klassische IRLL-Methode als zuverlässige Sequenz zur T1-Quantifizierung etabliert. In der klinischen Bildgebung werden auf Grund der unterschiedlichen Zeitskalen und anderer technischer Voraussetzungen andere Anforderungen an die Datenakquisition gestellt. Dabei hat in den letzten Jahren die MOLLI-Sequenz große Verbreitung gefunden. Sie ist eine Abwandlung der IRLL-Sequenz, bei der mit einer bSSFP-Bildgebungssequenz getriggert ganze Bilder während eines Herzschlages aufgenommen werden. Die MOLLI-Sequenz reagiert dabei empfindlich auf die Wartezeiten zwischen den einzelnen Transienten. Um mit diese Problematik in den Griff zu bekommen und gleichzeitig die Meßzeit verkürzen zu können wurde eine neue Methode zum Fitten der Daten entwickelt, welche die Abhängigkeit der scheinbaren Relaxationszeit von der Wartezeit zwischen den einzelnen Transienten, sowie der mittleren Herzrate fast vollständig eliminiert. Diese Methode liefert für das ganze klinisch Spektrum an erwarteten T1-Zeiten, vor und nach Kontrastmittelgabe, stabile Ergebnisse und erlaubte ein deutliche Verkürzung der Meßzeit, ohne die Anzahl der aufgenommenen Meßzeitpunkte zu reduzieren. Dies wurde in einer initialen klinischen Studie genutzt, um ECV-Werte in Patienten zu bestimmen.
Ein Nachteil der Verwendung der MOLLI-Sequenz ist, daß nur die scheinbare Relaxationszeit aus den Fit der Meßdaten bestimmt wird. Die standardmäßig genutzte Korrektur benutzt aber dem gefitteten Wert der Gleichgewichtsmagnetisierung um den wahren T1-Wert zu bestimmen. Somit ist es für die Bestimmung des T1-Wertes notwendig, die Qualität der Inversionspräparation zu kennen. Auf Basis der neuen Fitmethode wurde eine Anpassung der MOLLI-Sequenz demonstriert, welche die Bestimmung der Gleichgewichtsmagnetisierung unabhängig von der Qualität der Inversionspräparation erlaubt. Dafür verlängert sich die Meßdauer lediglich um einen Herzschlag um in geeigneter Weise ein zusätzliches Bild aufnehmen zu können.
Abschließend wurde in dieser Arbeit der Signal-Zeit-Verlauf der MOLLI-Sequenz eingehend theoretische untersucht um ein besseres Verständnis der getriggerten IRLL-Sequenzen zu entwickeln. In diesem Zusammenhang konnte eine einfache Interpretation der scheinbaren Relaxationszeit gefunden werden. Ebenso konnte erklärt werden, warum die für ungetriggerte IRLL-Sequenzen abgeleitete Korrekturgleichung auch im getriggerten Fall erstaunlich gute Ergebnisse liefert. Weiterhin konnten Fehlerquellen für die verbleibenden Abweichungen identifiziert werden, welche als Ausgangspunkt für die Ableitung verbesserter Korrekturgleichungen genutzt werden können.