Refine
Has Fulltext
- yes (15)
Is part of the Bibliography
- yes (15)
Year of publication
Document Type
- Journal article (14)
- Doctoral Thesis (1)
Keywords
- mouse (4)
- magnetic resonance imaging (3)
- 4D flow (2)
- WSS (2)
- aortic arch (2)
- atherosclerosis (2)
- in vivo imaging (2)
- macrophages (2)
- mice (2)
- quantification (2)
Institute
- Physikalisches Institut (14)
- Institut für diagnostische und interventionelle Neuroradiologie (ehem. Abteilung für Neuroradiologie) (5)
- Medizinische Klinik und Poliklinik I (5)
- 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)
- Deutsches Zentrum für Herzinsuffizienz (DZHI) (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.
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.
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.
Magnetic particle imaging is an emerging tomographic method used for evaluation of the spatial distribution of iron‐oxide nanoparticles. In this work, the effect of the polymer coating on the response of particles was studied. Particles with covalently crosslinked coating showed improved signal and image resolution.
CNS imaging characteristics in fibromyalgia patients with and without peripheral nerve involvement
(2022)
We tested the hypothesis that reduced skin innervation in fibromyalgia syndrome is associated with specific CNS changes. This prospective case–control study included 43 women diagnosed with fibromyalgia syndrome and 40 healthy controls. We further compared the fibromyalgia subgroups with reduced (n = 21) and normal (n = 22) skin innervation. Brains were analysed for cortical volume, for white matter integrity, and for functional connectivity. Compared to controls, cortical thickness was decreased in regions of the frontal, temporal and parietal cortex in the fibromyalgia group as a whole, and decreased in the bilateral pericalcarine cortices in the fibromyalgia subgroup with reduced skin innervation. Diffusion tensor imaging revealed a significant increase in fractional anisotropy in the corona radiata, the corpus callosum, cingulum and fornix in patients with fibromyalgia compared to healthy controls and decreased FA in parts of the internal capsule and thalamic radiation in the subgroup with reduced skin innervation. Using resting-state fMRI, the fibromyalgia group as a whole showed functional hypoconnectivity between the right midfrontal gyrus and the posterior cerebellum and the right crus cerebellum, respectively. The subgroup with reduced skin innervation showed hyperconnectivity between the inferior frontal gyrus, the angular gyrus and the posterior parietal gyrus. Our results suggest that the subgroup of fibromyalgia patients with pronounced pathology in the peripheral nervous system shows alterations in morphology, structural and functional connectivity also at the level of the encephalon. We propose considering these subgroups when conducting clinical trials.
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.
Simultaneous measurements of 3D wall shear stress and pulse wave velocity in the murine aortic arch
(2021)
Purpose
Wall shear stress (WSS) and pulse wave velocity (PWV) are important parameters to characterize blood flow in the vessel wall. Their quantification with flow-sensitive phase-contrast (PC) cardiovascular magnetic resonance (CMR), however, is time-consuming. Furthermore, the measurement of WSS requires high spatial resolution, whereas high temporal resolution is necessary for PWV measurements. For these reasons, PWV and WSS are challenging to measure in one CMR session, making it difficult to directly compare these parameters. By using a retrospective approach with a flexible reconstruction framework, we here aimed to simultaneously assess both PWV and WSS in the murine aortic arch from the same 4D flow measurement.
Methods
Flow was measured in the aortic arch of 18-week-old wildtype (n = 5) and ApoE\(^{−/−}\) mice (n = 5) with a self-navigated radial 4D-PC-CMR sequence. Retrospective data analysis was used to reconstruct the same dataset either at low spatial and high temporal resolution (PWV analysis) or high spatial and low temporal resolution (WSS analysis). To assess WSS, the aortic lumen was labeled by semi-automatically segmenting the reconstruction with high spatial resolution. WSS was determined from the spatial velocity gradients at the lumen surface. For calculation of the PWV, segmentation data was interpolated along the temporal dimension. Subsequently, PWV was quantified from the through-plane flow data using the multiple-points transit-time method. Reconstructions with varying frame rates and spatial resolutions were performed to investigate the influence of spatiotemporal resolution on the PWV and WSS quantification.
Results
4D flow measurements were conducted in an acquisition time of only 35 min. Increased peak flow and peak WSS values and lower errors in PWV estimation were observed in the reconstructions with high temporal resolution. Aortic PWV was significantly increased in ApoE\(^{−/−}\) mice compared to the control group (1.7 ± 0.2 versus 2.6 ± 0.2 m/s, p < 0.001). Mean WSS magnitude values averaged over the aortic arch were (1.17 ± 0.07) N/m\(^2\) in wildtype mice and (1.27 ± 0.10) N/m\(^2\) in ApoE\(^{−/−}\) mice.
Conclusion
The post processing algorithm using the flexible reconstruction framework developed in this study permitted quantification of global PWV and 3D-WSS in a single acquisition. The possibility to assess both parameters in only 35 min will markedly improve the analyses and information content of in vivo measurements.
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.
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.