Refine
Has Fulltext
- yes (17)
Is part of the Bibliography
- yes (17)
Year of publication
Document Type
- Journal article (16)
- Doctoral Thesis (1)
Keywords
- mouse (4)
- MOLLI (3)
- magnetic resonance imaging (3)
- 4D flow (2)
- WSS (2)
- aortic arch (2)
- atherosclerosis (2)
- in vivo imaging (2)
- macrophages (2)
- mice (2)
Institute
- Physikalisches Institut (16)
- Medizinische Klinik und Poliklinik I (6)
- Institut für diagnostische und interventionelle Neuroradiologie (ehem. Abteilung für Neuroradiologie) (5)
- Deutsches Zentrum für Herzinsuffizienz (DZHI) (3)
- Institut für Experimentelle Biomedizin (2)
- Neurologische Klinik und Poliklinik (2)
- Abteilung für Funktionswerkstoffe der Medizin und der Zahnheilkunde (1)
- Abteilung für Molekulare Innere Medizin (in der Medizinischen Klinik und Poliklinik II) (1)
- Institut für Molekulare Infektionsbiologie (1)
- Institut für Virologie und Immunbiologie (1)
Sonstige beteiligte Institutionen
In Vivo Imaging of Stepwise Vessel Occlusion in Cerebral Photothrombosis of Mice by \(^{19}\)F MRI
(2011)
Background
\(^{19}\)F magnetic resonance imaging (MRI) was recently introduced as a promising technique for in vivo cell tracking. In the present study we compared \(^{19}\)F MRI with iron-enhanced MRI in mice with photothrombosis (PT) at 7 Tesla. PT represents a model of focal cerebral ischemia exhibiting acute vessel occlusion and delayed neuroinflammation.
Methods/Principal Findings
Perfluorocarbons (PFC) or superparamagnetic iron oxide particles (SPIO) were injected intravenously at different time points after photothrombotic infarction. While administration of PFC directly after PT induction led to a strong \(^{19}\)F signal throughout the entire lesion, two hours delayed application resulted in a rim-like \(^{19}\)F signal at the outer edge of the lesion. These findings closely resembled the distribution of signal loss on T2-weighted MRI seen after SPIO injection reflecting intravascular accumulation of iron particles trapped in vessel thrombi as confirmed histologically. By sequential administration of two chemically shifted PFC compounds 0 and 2 hours after illumination the different spatial distribution of the \(^{19}\)F markers (infarct core/rim) could be visualized in the same animal. When PFC were applied at day 6 the fluorine marker was only detected after long acquisition times ex vivo. SPIO-enhanced MRI showed slight signal loss in vivo which was much more prominent ex vivo indicative for neuroinflammation at this late lesion stage.
Conclusion
Our study shows that vessel occlusion can be followed in vivo by \(^{19}\)F and SPIO-enhanced high-field MRI while in vivo imaging of neuroinflammation remains challenging. The timing of contrast agent application was the major determinant of the underlying processes depicted by both imaging techniques. Importantly, sequential application of different PFC compounds allowed depiction of ongoing vessel occlusion from the core to the margin of the ischemic lesions in a single MRI measurement.
Background
Oncolytic virotherapy of tumors is an up-coming, promising therapeutic modality of cancer therapy. Unfortunately, non-invasive techniques to evaluate the inflammatory host response to treatment are rare. Here, we evaluate \(^{19}\)F magnetic resonance imaging (MRI) which enables the non-invasive visualization of inflammatory processes in pathological conditions by the use of perfluorocarbon nanoemulsions (PFC) for monitoring of oncolytic virotherapy.
Methodology/Principal Findings
The Vaccinia virus strain GLV-1h68 was used as an oncolytic agent for the treatment of different tumor models. Systemic application of PFC emulsions followed by \(^1H\)/\(^{19}\)F MRI of mock-infected and GLV-1h68-infected tumor-bearing mice revealed a significant accumulation of the \(^{19}\)F signal in the tumor rim of virus-treated mice. Histological examination of tumors confirmed a similar spatial distribution of the \(^{19}\)F signal hot spots and \(CD68^+\)-macrophages. Thereby, the \(CD68^+\)-macrophages encapsulate the GFP-positive viral infection foci. In multiple tumor models, we specifically visualized early inflammatory cell recruitment in Vaccinia virus colonized tumors. Furthermore, we documented that the \(^{19}\)F signal correlated with the extent of viral spreading within tumors.
Conclusions/Significance
These results suggest \(^{19}\)F MRI as a non-invasive methodology to document the tumor-associated host immune response as well as the extent of intratumoral viral replication. Thus, \(^{19}\)F MRI represents a new platform to non-invasively investigate the role of the host immune response for therapeutic outcome of oncolytic virotherapy and individual patient response.
Background
The aortic pulse-wave velocity (PWV) is an important indicator of cardiovascular risk. In recent studies MRI methods have been developed to measure this parameter noninvasively in mice. Present techniques require additional hardware for cardiac and respiratory gating. In this work a robust self-gated measurement of the local PWV in mice without the need of triggering probes is proposed.
Methods
The local PWV of 6-months-old wild-type C57BL/6J mice (n=6) was measured in the abdominal aorta with a retrospectively triggered radial Phase Contrast (PC) MR sequence using the flow-area (QA) method. A navigator signal was extracted from the CMR data of highly asymmetric radial projections with short repetition time (TR=3 ms) and post-processed with high-pass and low-pass filters for retrospective cardiac and respiratory gating. The self-gating signal was used for a reconstruction of high-resolution Cine frames of the aortic motion. To assess the local PWV the volume flow Q and the cross-sectional area A of the aorta were determined. The results were compared with the values measured with a triggered Cartesian and an undersampled triggered radial PC-Cine sequence.
Results
In all examined animals a self-gating signal could be extracted and used for retrospective breath-gating and PC-Cine reconstruction. With the non-triggered measurement PWV values of 2.3±0.2 m/s were determined. These values are in agreement with those measured with the triggered Cartesian (2.4±0.2 m/s) and the triggered radial (2.3±0.2 m/s) measurement. Due to the strong robustness of the radial trajectory against undersampling an acceleration of more than two relative to the prospectively triggered Cartesian sampling could be achieved with the retrospective method.
Conclusion
With the radial flow-encoding sequence the extraction of a self-gating signal is feasible. The retrospective method enables a robust and fast measurement of the local PWV without the need of additional trigger hardware.
Microstructural Analysis of Peripheral Lung Tissue through CPMG Inter-Echo Time R2 Dispersion
(2015)
Since changes in lung microstructure are important indicators for (early stage) lung pathology, there is a need for quantifiable information of diagnostically challenging cases in a clinical setting, e.g. to evaluate early emphysematous changes in peripheral lung tissue. Considering alveoli as spherical air-spaces surrounded by a thin film of lung tissue allows deriving an expression for Carr-Purcell-Meiboom-Gill transverse relaxation rates R-2 with a dependence on inter-echo time, local air-tissue volume fraction, diffusion coefficient and alveolar diameter, within a weak field approximation. The model relaxation rate exhibits the same hyperbolic tangent dependency as seen in the Luz-Meiboom model and limiting cases agree with Brooks et al. and Jensen et al. In addition, the model is tested against experimental data for passively deflated rat lungs: the resulting mean alveolar radius of RA = 31.46 \(\pm\) 13.15 \(\mu\)m is very close to the literature value (similar to 34 \(\mu\)m). Also, modeled radii obtained from relaxometer measurements of ageing hydrogel foam (that mimics peripheral lung tissue) are in good agreement with those obtained from mu CT images of the same foam (mean relative error: 0.06 \(\pm\) 0.01). The model's ability to determine the alveolar radius and/or air volume fraction will be useful in quantifying peripheral lung microstructure.
Orthogonality, Lommel integrals and cross product zeros of linear combinations of Bessel functions
(2015)
The cylindrical Bessel differential equation and the spherical Bessel differential equation in the interval R\(\leq\)r\(\leq\)\(\gamma\)R with Neumann boundary conditions are considered. The eigenfunctions are linear combinations of the Bessel function \(\Phi\)\(_{n,ν}\)(r) = Y'\(_{ν}\) (\(\lambda\)\(_{n,ν}\))J\(_{ν}\)(\(\lambda\)\(_{n,ν}\) r/R) - J'\(_{ν}\)(\(\lambda\)\(_{n,ν}\))Y\(_{ν}\)(\(\lambda\)\(_{n,ν}\)r/R) or linear combinations of the spherical Bessel functions \(\psi\)\(_{m,ν}\)(r) = y'\(_{ν}\)(\(\lambda\)\(_{m,ν}\))j\(_{ν}\)(\(\lambda\)\(_{m,ν}\)r/R) - j'\(_{ν}\)(\(\lambda\)\(_{m,ν}\))y\(_{ν}\)(\(\lambda\)\(_{m,ν}\)r/R). The orthogonality relations with analytical expressions for the normalization constant are given. Explicit expressions for the Lommel integrals in terms of Lommel functions are derived. The cross product zeros Y'\(_{ν}\)\(\lambda\)\(_{n,ν}\))J'\(_{ν}\)(\(\gamma\)\(\lambda\)\(_{n,ν}\))- J'\(_{ν}\)(\(\lambda\)\(_{n,ν}\))Y'\(_{ν}\)(\(\gamma\)\(\lambda\)\(_{n,ν}\)) = 0 and y'\(_{ν}\)(\(\lambda\)\(_{m,ν}\))j'\(_{ν}\)(\(\gamma\)\(\lambda\)\(_{m,ν}\)) - j'\(_{ν}\)(\(\lambda\)\(_{m,ν}\))y'\(_{ν}\)(\(\gamma\)\(\lambda\)\(_{m,ν}\)) = 0 are considered in the complex plane for real as well as complex values of the index ν and approximations for the exceptional zero \(\lambda\)\(_{1,ν}\) are obtained. A numerical scheme based on the discretization of the twodimensional and three-dimensional Laplace operator with Neumann boundary conditions is presented. Explicit representations of the radial part of the Laplace operator in form of a tridiagonal matrix allow the simple computation of the cross product zeros.
In biological tissue, an accumulation of similarly shaped objects with a susceptibility difference to the surrounding tissue generates a local distortion of the external magnetic field in magnetic resonance imaging. It induces stochastic field fluctuations that characteristically influence proton spin dephasing in the vicinity of these magnetic perturbers. The magnetic field correlation that is associated with such local magnetic field inhomogeneities can be expressed in the form of a dynamic frequency autocorrelation function that is related to the time evolution of the measured magnetization. Here, an eigenfunction expansion for two simple magnetic perturber shapes, that of spheres and cylinders, is considered for restricted spin diffusion in a simple model geometry. Then, the concept of generalized moment analysis, an approximation technique that is applied in the study of (non-)reactive processes that involve Brownian motion, allows deriving analytical expressions of the correlation function for different exponential decay forms. Results for the biexponential decay for both spherical and cylindrical magnetized objects are derived and compared with the frequently used (less accurate) monoexponential decay forms. They are in asymptotic agreement with the numerically exact value of the correlation function for long and short times.
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.
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.
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.
Purpose
4D flow cardiovascular magnetic resonance (CMR) and the assessment of wall shear stress (WSS) are non-invasive tools to study cardiovascular risks in vivo. Major limitations of conventional triggered methods are the long measurement times needed for high-resolution data sets and the necessity of stable electrocardiographic (ECG) triggering. In this work an ECG-free retrospectively synchronized method is presented that enables accelerated high-resolution measurements of 4D flow and WSS in the aortic arch of mice.
Methods
4D flow and WSS were measured in the aortic arch of 12-week-old wildtype C57BL/6 J mice (n = 7) with a radial 4D-phase-contrast (PC)-CMR sequence, which was validated in a flow phantom. Cardiac and respiratory motion signals were extracted from the radial CMR signal and were used for the reconstruction of 4D-flow data. Rigid motion correction and a first order B0 correction was used to improve the robustness of magnitude and velocity data.
The aortic lumen was segmented semi-automatically. Temporally averaged and time-resolved WSS and oscillatory shear index (OSI) were calculated from the spatial velocity gradients at the lumen surface at 14 locations along the aortic arch. Reproducibility was tested in 3 animals and the influence of subsampling was investigated.
Results
Volume flow, cross-sectional areas, WSS and the OSI were determined in a measurement time of only 32 min. Longitudinal and circumferential WSS and radial stress were assessed at 14 analysis planes along the aortic arch. The average longitudinal, circumferential and radial stress values were 1.52 ± 0.29 N/m2, 0.28 ± 0.24 N/m2 and − 0.21 ± 0.19 N/m2, respectively. Good reproducibility of WSS values was observed.
Conclusion
This work presents a robust measurement of 4D flow and WSS in mice without the need of ECG trigger signals. The retrospective approach provides fast flow quantification within 35 min and a flexible reconstruction framework.