Detection of pathological contrast enhancement with synthetic brain imaging from quantitative multiparametric MRI

BACKGROUND AND PURPOSE

We aimed to test whether synthetic T1-weighted imaging derived from a post-contrast Quantitative Transient-state Imaging (QTI) acquisition enabled revealing pathological contrast enhancement in intracranial lesions.

METHODS

The analysis included 141 patients who underwent a 3 Tesla-MRI brain exam with intravenous contrast media administration, with the post-contrast acquisition protocol comprising a three-dimensional fast spoiled gradient echo (FSPGR) sequence and a QTI acquisition. Synthetic T1-weighted images were generated from QTI-derived quantitative maps of relaxation times and proton density. Two neuroradiologists assessed synthetic and conventional post-contrast T1-weighted images for the presence and pattern of pathological contrast enhancement in intracranial lesions. Enhancement volumes were quantitatively compared.

RESULTS

Using conventional imaging as a reference, synthetic T1-weighted imaging was 93% sensitive in revealing the presence of contrast enhancing lesions. The agreement for the presence/absence of contrast enhancement was almost perfect both between readers (k = 1 for both conventional and synthetic imaging) and between sequences (k = 0.98 for both readers). In 91% of lesions, synthetic T1-weighted imaging showed the same pattern of contrast enhancement visible in conventional imaging. Differences in enhancement pattern in the remaining lesions can be due to the lower spatial resolution and the longer acquisition delay from contrast media administration of QTI compared to FSPGR. Overall, enhancement volumes appeared larger in synthetic imaging.

CONCLUSIONS

QTI-derived post-contrast synthetic T1-weighted imaging captures pathological contrast enhancement in most intracranial enhancing lesions. Further comparative studies employing quantitative imaging with higher spatial resolution is needed to support our data and explore possible future applications in clinical trials.

Magnetic resonance fingerprinting-based myelin water fraction mapping for the assessment of white matter maturation and integrity in typical development and leukodystrophies

A quantitative biomarker for myelination, such as myelin water fraction (MWF), would boost the understanding of normative and pathological neurodevelopment, improving patients' diagnosis and follow-up. We quantified the fraction of a rapidly relaxing pool identified as MW using multicomponent three-dimensional (3D) magnetic resonance fingerprinting (MRF) to evaluate white matter (WM) maturation in typically developing (TD) children and alterations in leukodystrophies (LDs). We acquired DTI and 3D MRF-based R1, R2 and MWF data of 15 TD children and 17 LD patients (9 months-12.5 years old) at 1.5 T. We computed normative maturation curves in corpus callosum and corona radiata and performed WM tract profile analysis, comparing MWF with R1, R2 and fractional anisotropy (FA). Normative maturation curves demonstrated a steep increase for all tissue parameters in the first 3 years of age, followed by slower growth for MWF while R1, R2R2 and FA reached a plateau. Unlike FA, MWF values were similar for regions of interest (ROIs) with different degrees of axonal packing, suggesting independence from fiber bundle macro-organization and higher myelin specificity. Tract profile analysis indicated a specific spatial pattern of myelination in the major fiber bundles, consistent across subjects. LD were better distinguished from TD by MWF rather than FA, showing reduced MWF with respect to age-matched controls in both ROI-based and tract analysis. In conclusion, MRF-based MWF provides myelin-specific WM maturation curves and is sensitive to alteration due to LDs, suggesting its potential as a biomarker for WM disorders. As MRF allows fast simultaneous acquisition of relaxometry and MWF, it can represent a valuable diagnostic tool to study and follow up developmental WM disorders in children.

Fast high-resolution electric properties tomography using three-dimensional quantitative transient-state imaging-based water fraction estimation

In this study, we aimed to develop a fast and robust high-resolution technique for clinically feasible electrical properties tomography based on water content maps (wEPT) using Quantitative Transient-state Imaging (QTI), a multiparametric transient state-based method that is similar to MR fingerprinting. Compared with the original wEPT implementation based on standard spin-echo acquisition, QTI provides robust electrical properties quantification towards B1 + inhomogeneities and full quantitative relaxometry data. To validate the proposed approach, 3D QTI data of 12 healthy volunteers were acquired on a 1.5 T scanner. QTI-provided T1 maps were used to compute water content maps of the tissues using an empirical relationship based on literature ex-vivo measurements. Assuming that electrical properties are modulated mainly by tissue water content, the water content maps were used to derive electrical conductivity and relative permittivity maps. The proposed technique was compared with a conventional phase-only Helmholtz EPT (HH-EPT) acquisition both within whole white matter, gray matter, and cerebrospinal fluid masks, and within different white and gray matter subregions. In addition, QTI-based wEPT was retrospectively applied to four multiple sclerosis adolescent and adult patients, compared with conventional contrast-weighted imaging in terms of lesion delineation, and quantitatively assessed by measuring the variation of electrical properties in lesions. Results obtained with the proposed approach agreed well with theoretical predictions and previous in vivo findings in both white and gray matter. The reconstructed maps showed greater anatomical detail and lower variability compared with standard phase-only HH-EPT. The technique can potentially improve delineation of pathology when compared with conventional contrast-weighted imaging and was able to detect significant variations in lesions with respect to normal-appearing tissues. In conclusion, QTI can reliably measure conductivity and relative permittivity of brain tissues within a short scan time, opening the way to the study of electric properties in clinical settings.

Quantitative T1 mapping detects blood-brain barrier breakdown in apparently non-enhancing multiple sclerosis lesions

OBJECTIVES

The disruption of the blood-brain barrier (BBB) is a key and early feature in the pathogenesis of demyelinating multiple sclerosis (MS) lesions and has been neuropathologically demonstrated in both active and chronic plaques. The local overt BBB disruption in acute demyelinating lesions is captured as signal hyperintensity in post-contrast T1-weighted images because of the contrast-related shortening of the T1 relaxation time. On the contrary, the subtle BBB disruption in chronic lesions is not visible at conventional radiological evaluation but it might be of clinical relevance. Indeed, persistent, subtle BBB leakage might be linked to low-grade inflammation and plaque evolution. Here we hypothesised that 3D Quantitative Transient-state Imaging (QTI) was able to reveal and measure T1 shortening (ΔT1) reflecting small amounts of contrast media leakage in apparently non-enhancing lesions (ANELs).

MATERIALS AND METHODS

Thirty-four patients with relapsing remitting MS were included in the study. All patients underwent a 3 T MRI exam of the brain including conventional sequences and QTI acquisitions (1.1 mm isotropic voxel) performed both before and after contrast media administration. For each patient, a ΔT1 map was obtained via voxel-wise subtraction of pre- and post- contrast QTI-derived T1 maps. ΔT1 values measured in ANELs were compared with those recorded in enhancing lesions and in the normal appearing white matter. A reference distribution of ΔT1 in the white matter was obtained from datasets acquired in 10 non-MS patients with unrevealing MR imaging.

RESULTS

Mean ΔT1 in ANELs (57.45 ± 48.27 ms) was significantly lower than in enhancing lesions (297.71 ± 177.52 ms; p < 0. 0001) and higher than in the normal appearing white matter (36.57 ± 10.53 ms; p < 0.005). Fifty-two percent of ANELs exhibited ΔT1 higher than those observed in the white matter of non-MS patients.

CONCLUSIONS

QTI-derived quantitative ΔT1 mapping enabled to measure contrast-related T1 shortening in ANELs. ANELs exhibiting ΔT1 values that deviate from the reference distribution in non-MS patients may indicate persistent, subtle, BBB disruption. Access to this information may be proved useful to better characterise pathology and objectively monitor disease activity and response to therapy.

Generating Synthetic Radiological Images with PySynthMRI: An Open-Source Cross-Platform Tool

Synthetic MR Imaging allows for the reconstruction of different image contrasts from a single acquisition, reducing scan times. Commercial products that implement synthetic MRI are used in research. They rely on vendor-specific acquisitions and do not include the possibility of using custom multiparametric imaging techniques. We introduce PySynthMRI, an open-source tool with a user-friendly interface that uses a set of input images to generate synthetic images with diverse radiological contrasts by varying representative parameters of the desired target sequence, including the echo time, repetition time and inversion time(s). PySynthMRI is written in Python 3.6, and it can be executed under Linux, Windows, or MacOS as a python script or an executable. The tool is free and open source and is developed while taking into consideration the possibility of software customization by the end user. PySynthMRI generates synthetic images by calculating the pixelwise signal intensity as a function of a set of input images (e.g., T1 and T2 maps) and simulated scanner parameters chosen by the user via a graphical interface. The distribution provides a set of default synthetic contrasts, including T1w gradient echo, T2w spin echo, FLAIR and Double Inversion Recovery. The synthetic images can be exported in DICOM or NiFTI format. PySynthMRI allows for the fast synthetization of differently weighted MR images based on quantitative maps. Specialists can use the provided signal models to retrospectively generate contrasts and add custom ones. The modular architecture of the tool can be exploited to add new features without impacting the codebase.

Simultaneous relaxometry and morphometry of human brain structures with 3D magnetic resonance fingerprinting: a multicenter, multiplatform, multifield-strength study

Relaxation times and morphological information are fundamental magnetic resonance imaging-derived metrics of the human brain that reflect the status of the underlying tissue. Magnetic resonance fingerprinting (MRF) enables simultaneous acquisition of T1 and T2 maps inherently aligned to the anatomy, allowing whole-brain relaxometry and morphometry in a single scan. In this study, we revealed the feasibility of 3D MRF for simultaneous brain structure-wise morphometry and relaxometry. Comprehensive test-retest scan analyses using five 1.5-T and three 3.0-T systems from a single vendor including different scanner types across 3 institutions demonstrated that 3D MRF-derived morphological information and relaxation times are highly repeatable at both 1.5 T and 3.0 T. Regional cortical thickness and subcortical volume values showed high agreement and low bias across different field strengths. The ability to acquire a set of regional T1, T2, thickness, and volume measurements of neuroanatomical structures with high repeatability and reproducibility facilitates the ability of longitudinal multicenter imaging studies to quantitatively monitor changes associated with underlying pathologies, disease progression, and treatments.

Quantitative imaging metrics derived from magnetic resonance fingerprinting using ISMRM/NIST MRI system phantom: An international multicenter repeatability and reproducibility study

PURPOSE

To compare the bias and inherent reliability of the quantitative (T₁ and T₂ ) imaging metrics generated from the magnetic resonance fingerprinting (MRF) technique using the ISMRM/NIST system phantom in an international multicenter setting.

METHOD

ISMRM/NIST MRI system phantom provides standard reference T₁ and T₂ relaxation values (vendor-provided) for each of the 14 vials in T₁ and T₂ arrays. MRF-SSFP scans repeated over 30 days on GE 1.5 and 3.0 T scanners at three collaborative centers. MRF estimated T_1, and T₂ values averaged over 30 days were compared with the phantom vendor-provided and spin-echo (SE) based convention gold standard (GS) method. Repeatability and reproducibility were characterized by the within-case coefficient of variation (wCV) of the MRF data acquired over 30 days, along with the biases.

RESULT

For the wide ranges of MRF estimated T₁ values, vials #1-8 (T₁ relaxation time between 2033 and 184 ms) exhibited a wCV less than 3% and 4%, respectively, on 3.0 and 1.5 T scanners. T₂ values in vials #1-8 (T₂ relaxation, 1044-45 ms) have shown wCV to be <7% on both 3.0 and 1.5 T scanners. A stronger linear correlation overall for T₁ (R²  = 0.9960 and 0.9963 at center-1 and center-2 on 3.0 T scanner, and R²  = 0.9951 and 0.9988 at center-1 and center-3 on 1.5 T scanner) compared to T₂ (R²  = 0.9971 and 0.9972 at center-1 and center-2 on 3.0 T scanner, and R²  = 0.9815 and 0.9754 at center-1 and center-3 on 1.5 T scanner). Bland-Altman (BA) analysis showed MRF based T₁ and T₂ values were within the limit of agreement (LOA) except for one data point. The mean difference or bias and 95% lower bound (LB) and upper bound (UB) LOA are reported in the format; mean bias: 95% LB LOA: 95% UB LOA. The biases for T₁ values were 21.34: -50.00: 92.69, 21.32: -47.29: 89.94 ms, and for T₂ values were -19.88: -42.37: 2.61, -19.06: -43.58: 5.45 ms on 3.0 T scanner at center-1 and center-2, respectively. Similarly, on 1.5 T scanner biases for T₁ values were 26.54: -53.41: 106.50, 9.997: -51.94: 71.94 ms, and for T₂ values were -23.84: -135.40: 87.76, -37.30: 134.30: 59.73 ms at center-1 and center-3, respectively. Additionally, the correlation between the SE based GS and MRF estimated T₁ and T₂ values (R²  = 0.9969 and 0.9977) showed a similar trend as we observed between vendor-provided and MRF estimated T₁ and T₂ values (R²  = 0.9963 and 0.9972). In addition to correlation, BA analysis showed that all the vials are within the LOA between the GS and vendor-provided for the T₁ values and except one vial for T₂ . All the vials are within the LOA between GS and MRF except one vial in T₁ and T₂ array. The wCV for reproducibility was <3% for both T₁ and T₂ values in vials #1-8, for all the 14 vials, wCV calculated for reproducibility was <4% for T₁ values and <5% for T₂ .

CONCLUSION

This study shows that MRF is highly repeatable (wCV <4% for T₁ and <7% for T₂ ) and reproducible (wCV < 3% for both T₁ and T₂ ) in certain vials (vials #1-8).

Three dimensional MRF obtains highly repeatable and reproducible multi-parametric estimations in the healthy human brain at 1.5T and 3T

Magnetic resonance fingerprinting (MRF) is highly promising as a quantitative MRI technique due to its accuracy, robustness, and efficiency. Previous studies have found high repeatability and reproducibility of 2D MRF acquisitions in the brain. Here, we have extended our investigations to 3D MRF acquisitions covering the whole brain using spiral projection k-space trajectories. Our travelling head study acquired test/retest data from the brains of 12 healthy volunteers and 8 MRI systems (3 systems at 3 T and 5 at 1.5 T, all from a single vendor), using a study design not requiring all subjects to be scanned at all sites. The pulse sequence and reconstruction algorithm were the same for all acquisitions. After registration of the MRF-derived PD T1 and T2 maps to an anatomical atlas, coefficients of variation (CVs) were computed to assess test/retest repeatability and inter-site reproducibility in each voxel, while a General Linear Model (GLM) was used to determine the voxel-wise variability between all confounders, which included test/retest, subject, field strength and site. Our analysis demonstrated a high repeatability (CVs 0.7-1.3% for T1, 2.0-7.8% for T2, 1.4-2.5% for normalized PD) and reproducibility (CVs of 2.0-5.8% for T1, 7.4-10.2% for T2, 5.2-9.2% for normalized PD) in gray and white matter. Both repeatability and reproducibility improved when compared to similar experiments using 2D acquisitions. Three-dimensional MRF obtains highly repeatable and reproducible estimations of T1 and T2, supporting the translation of MRF-based fast quantitative imaging into clinical applications.

The effect of gadolinium-based contrast agent administration on magnetic resonance fingerprinting-based T1 relaxometry in patients with prostate cancer

Magnetic resonance fingerprinting (MRF) is a rapidly developing fast quantitative mapping technique able to produce multiple property maps with reduced sensitivity to motion. MRF has shown promise in improving the diagnosis of clinically significant prostate cancer but requires further validation as part of a prostate multiparametric (mp) MRI protocol. mpMRI protocol mandates the inclusion of dynamic contrast enhanced (DCE) imaging, known for its significant T1 shortening effect. MRF could be used to measure both pre- and post-contrast T1 values, but its utility must be assessed. In this proof-of-concept study, we sought to evaluate the variation in MRF T1 measurements post gadolinium-based contrast agent (GBCA) injection and the utility of such T1 measurements to differentiate peripheral and transition zone tumours from normal prostatic tissue. We found that the T1 variation in all tissues increased considerably post-GBCA following the expected significant T1 shortening effect, compromising the ability of MRF T1 to identify transition zone lesions. We, therefore, recommend performing MRF T1 prior to DCE imaging to maintain its benefit for improving detection of both peripheral and transition zone lesions while reducing additional scanning time. Demonstrating the effect of GBCA on MRF T1 relaxometry in patients also paves the way for future clinical studies investigating the added value of post-GBCA MRF in PCa, including its dynamic analysis as in DCE-MRF.

Repeatability and reproducibility of human brain morphometry using three-dimensional magnetic resonance fingerprinting

Three-dimensional (3D) Magnetic resonance fingerprinting (MRF) permits whole-brain volumetric quantification of T1 and T2 relaxation values, potentially replacing conventional T1-weighted structural imaging for common brain imaging analysis. The aim of this study was to evaluate the repeatability and reproducibility of 3D MRF in evaluating brain cortical thickness and subcortical volumetric analysis in healthy volunteers using conventional 3D T1-weighted images as a reference standard. Scan-rescan tests of both 3D MRF and conventional 3D fast spoiled gradient recalled echo (FSPGR) were performed. For each sequence, the regional cortical thickness and volume of the subcortical structures were measured using standard automatic brain segmentation software. Repeatability and reproducibility were assessed using the within-subject coefficient of variation (wCV), intraclass correlation coefficient (ICC), and mean percent difference and ICC, respectively. The wCV and ICC of cortical thickness were similar across all regions with both 3D MRF and FSPGR. The percent relative difference in cortical thickness between 3D MRF and FSPGR across all regions was 8.0 ± 3.2%. The wCV and ICC of the volume of subcortical structures across all structures were similar between 3D MRF and FSPGR. The percent relative difference in the volume of subcortical structures between 3D MRF and FSPGR across all structures was 7.1 ± 3.6%. 3D MRF measurements of human brain cortical thickness and subcortical volumes are highly repeatable, and consistent with measurements taken on conventional 3D T1-weighted images. A slight, consistent bias was evident between the two, and thus careful attention is required when combining data from MRF and conventional acquisitions.

Magnetic resonance fingerprinting of the pancreas at 1.5 T and 3.0 T

Magnetic resonance imaging of the pancreas is increasingly used as an important diagnostic modality for characterisation of pancreatic lesions. Pancreatic MRI protocols are mostly qualitative due to time constraints and motion sensitivity. MR Fingerprinting is an innovative acquisition technique that provides qualitative data and quantitative parameter maps from a single free-breathing acquisition with the potential to reduce exam times. This work investigates the feasibility of MRF parameter mapping for pancreatic imaging in the presence of free-breathing exam. Sixteen healthy participants were prospectively imaged using MRF framework. Regions-of-interest were drawn in multiple solid organs including the pancreas and T1 and T2 values determined. MRF T1 and T2 mapping was performed successfully in all participants (acquisition time:2.4-3.6 min). Mean pancreatic T1 values were 37-43% lower than those of the muscle, spleen, and kidney at both 1.5 and 3.0 T. For these organs, the mean pancreatic T2 values were nearly 40% at 1.5 T and < 12% at 3.0 T. The feasibility of MRF at 1.5 T and 3 T was demonstrated in the pancreas. By enabling fast and free-breathing quantitation, MRF has the potential to add value during the clinical characterisation and grading of pathological conditions, such as pancreatitis or cancer.

7-T MRI of explanted liver and ex-vivo pancreatic specimens: prospective study protocol of radiological-pathological correlation feasibility (the EXLIPSE project)

The study focuses on radiological-pathological correlation between imaging of ex vivo samples obtained by a 7-T scanner and histological examination. The specimens will be derived from native explanted cirrhotic livers, liver grafts excluded from donation because of severe steatosis, and primary pancreatic tumours. Magnetic resonance imaging (MRI) examinations will be performed within 24 h from liver or pancreatic lesion surgical removal. The MRI protocol will include morphological sequences, quantitative T1, T2, and fat-, water-fraction maps with Cartesian k-space acquisition, and multiparametric methods based on a transient-state “MRI fingerprinting”. Finally, the specimen will be fixed by formalin. Qualitative imaging analysis will be performed by two independent blinded radiologists to assess image consistency score. Quantitative analysis will be performed by drawing regions of interest on different tissue zones to measure T1 and T2 relaxation times as well as fat- and water-fraction. The same tissue areas will be analysed by the pathologists. This study will provide the possibility to improve our knowledge about qualitative and quantitative abdominal imaging assessment at 7 T, by correlating imaging characteristics and the corresponding histological composition of ex vivo specimens, in order to identify imaging biomarkers. Trial registration: ClinicalTrials.gov: 13646. Registered 9 July 2019—retrospectively registered

Rapid three-dimensional multiparametric MRI with quantitative transient-state imaging

Novel methods for quantitative, transient-state multiparametric imaging are increasingly being demonstrated for assessment of disease and treatment efficacy. Here, we build on these by assessing the most common Non-Cartesian readout trajectories (2D/3D radials and spirals), demonstrating efficient anti-aliasing with a k-space view-sharing technique, and proposing novel methods for parameter inference with neural networks that incorporate the estimation of proton density. Our results show good agreement with gold standard and phantom references for all readout trajectories at 1.5 T and 3 T. Parameters inferred with the neural network were within 6.58% difference from the parameters inferred with a high-resolution dictionary. Concordance correlation coefficients were above 0.92 and the normalized root mean squared error ranged between 4.2 and 12.7% with respect to gold-standard phantom references for T1 and T2. In vivo acquisitions demonstrate sub-millimetric isotropic resolution in under five minutes with reconstruction and inference times < 7 min. Our 3D quantitative transient-state imaging approach could enable high-resolution multiparametric tissue quantification within clinically acceptable acquisition and reconstruction times.

Retrospective rigid motion correction of three-dimensional magnetic resonance fingerprinting of the human brain

PURPOSE

To obtain three-dimensional (3D), quantitative and motion-robust imaging with magnetic resonance fingerprinting (MRF).

METHODS

Our acquisition is based on a 3D spiral projection k-space scheme. We compared different orderings of trajectory interleaves in terms of rigid motion-correction robustness. In all tested orderings, we considered the whole dataset as a sum of 56 segments of 7-s duration, acquired sequentially with the same flip angle schedule. We performed a separate image reconstruction for each segment, producing whole-brain navigators that were aligned to the first segment using normalized correlation. The estimated rigid motion was used to correct the k-space data, and the aligned data were matched with the dictionary to obtain motion-corrected maps.

RESULTS

A significant improvement on the motion-affected maps after motion correction is evident with the suppression of motion artifacts. Correlation with the motionless baseline improved by 20% on average for both T₁ and T₂ estimations after motion correction. In addition, the average motion-induced quantification bias of 70 ms for T₁ and 18 ms for T₂ values was reduced to 12 ms and 6 ms, respectively, improving the reliability of quantitative estimations.

CONCLUSION

We established a method that allows correcting 3D rigid motion on a 7-s timescale during the reconstruction of MRF data using self-navigators, improving the image quality and the quantification robustness.

Flexible and efficient optimization of quantitative sequences using automatic differentiation of Bloch simulations

PURPOSE

To investigate a computationally efficient method for optimizing the Cramér-Rao Lower Bound (CRLB) of quantitative sequences without using approximations or an analytical expression of the signal.

METHODS

Automatic differentiation was applied to Bloch simulations and used to optimize several quantitative sequences without the need for approximations or an analytical expression. The results were validated with in vivo measurements and comparisons to prior art. Multi-echo spin echo and DESPO T 1 were used as benchmarks to verify the CRLB implementation. The CRLB of the Magnetic Resonance Fingerprinting (MRF) sequence, which has a complicated analytical formulation, was also optimized using automatic differentiation.

RESULTS

The sequence parameters obtained for multi-echo spin echo and DESPO T 1 matched results obtained using conventional methods. In vivo, MRF scans demonstrate that the CRLB optimization obtained with automatic differentiation can improve performance in presence of white noise. For MRF, the CRLB optimization converges in 1.1 CPU hours for N TR = 400 and has O ( N TR ) asymptotic runtime scaling for the calculation of the CRLB objective and gradient.

CONCLUSIONS

Automatic differentiation can be used to optimize the CRLB of quantitative sequences without using approximations or analytical expressions. For MRF, the runtime is computationally efficient and can be used to investigate confounding factors as well as MRF sequences with a greater number of repetitions.

Multi-site repeatability and reproducibility of MR fingerprinting of the healthy brain at 1.5 and 3.0 T

Fully-quantitative MR imaging methods are useful for longitudinal characterization of disease and assessment of treatment efficacy. However, current quantitative MRI protocols have not been widely adopted in the clinic, mostly due to lengthy scan times. Magnetic Resonance Fingerprinting (MRF) is a new technique that can reconstruct multiple parametric maps from a single fast acquisition in the transient state of the MR signal. Due to the relative novelty of this technique, the repeatability and reproducibility of quantitative measurements obtained using MRF has not been extensively studied. Our study acquired test/retest data from the brains of nine healthy volunteers, each scanned on five MRI systems (two at 3.0 T and three at 1.5 T, all from a single vendor) located at two different centers. The pulse sequence and reconstruction algorithm were the same for all acquisitions. After registration of the MRF-derived M0, T1 and T2 maps to an anatomical atlas, coefficients-of-variation (CVs) were computed to assess test/retest repeatability and inter-site reproducibility in each voxel, while a General Linear Model (GLM) was used to determine the voxel-wise variability between all confounders, which included test/retest, subject, field strength and site. Our analysis demonstrated an excellent repeatability (CVs of 2-3% for T1, 5-8% for T2, 3% for normalized-M0) and a good reproducibility (CVs of 3-8% for T1, 8-14% for T2, 5% for normalized-M0) in grey and white matter.

Feasibility of Quantitative Magnetic Resonance Fingerprinting in Ovarian Tumors for T1 and T2 Mapping in a PET/MR Setting

Multiparametric magnetic resonance imaging (MRI) can be used to characterize many cancer subtypes including ovarian cancer. Quantitative mapping of MRI relaxation values, such as T 1 and T 2 mapping, is promising for improving tumor assessment beyond conventional qualitative T 1- and T 2-weighted images. However, quantitative MRI relaxation mapping methods often involve long scan times due to sequentially measuring many parameters. Magnetic resonance fingerprinting (MRF) is a new method that enables fast quantitative MRI by exploiting the transient signals caused by the variation of pseudorandom sequence parameters. These transient signals are then matched to a simulated dictionary of T 1 and T 2 values to create quantitative maps. The ability of MRF to simultaneously measure multiple parameters, could represent a new approach to characterizing cancer and assessing treatment response. This feasibility study investigates MRF for simultaneous T 1, T 2, and relative proton density (rPD) mapping using ovarian cancer as a model system.

Designing contrasts for rapid, simultaneous parameter quantification and flow visualization with quantitative transient-state imaging

Magnetic resonance imaging (MRI) has evolved into an outstandingly versatile diagnostic modality, as it has the ability to non-invasively produce detailed information on a tissue's structure and function. Complementary data is normally obtained in separate measurements, either as contrast-weighted images, which are fast and simple to acquire, or as quantitative parametric maps, which offer an absolute quantification of underlying biophysical effects, such as relaxation times or flow. Here, we demonstrate how to acquire and reconstruct data in a transient-state with a dual purpose: 1 - to generate contrast-weighted images that can be adjusted to emphasise clinically relevant image biomarkers; exemplified with signal modulation according to flow to obtain angiography information, and 2 - to simultaneously infer multiple quantitative parameters with a single, highly accelerated acquisition. This is achieved by introducing three novel elements: a model that accounts for flowing blood, a method for sequence design using smooth flip angle excitation patterns that incorporates both parameter encoding and signal contrast, and the reconstruction of temporally resolved contrast-weighted images. From these images we simultaneously obtain angiography projections and multiple quantitative maps. By doing so, we increase the amount of clinically relevant data without adding measurement time, creating new dimensions for biomarker exploration and adding value to MR examinations for patients and clinicians alike.

A novel cyclic biased agonist of the apelin receptor, MM07, is disease modifying in the rat monocrotaline model of pulmonary arterial hypertension

BACKGROUND AND PURPOSE

Apelin is an endogenous vasodilatory and inotropic peptide that is down-regulated in human pulmonary arterial hypertension, although the density of the apelin receptor is not significantly attenuated. We hypothesised that a G protein-biased apelin analogue MM07, which is more stable than the endogenous apelin peptide, may be beneficial in this condition with the advantage of reduced β-arrestin-mediated receptor internalisation with chronic use.

EXPERIMENTAL APPROACH

Male Sprague-Dawley rats received either monocrotaline to induce pulmonary arterial hypertension or saline and then daily i.p. injections of either MM07 or saline for 21 days. The extent of disease was assessed by right ventricular catheterisation, cardiac MRI, and histological analysis of the pulmonary vasculature. The effect of MM07 on signalling, proliferation, and apoptosis of human pulmonary artery endothelial cells was investigated.

KEY RESULTS

MM07 significantly reduced the elevation of right ventricular systolic pressure and hypertrophy induced by monocrotaline. Monocrotaline-induced changes in cardiac structure and function, including right ventricular end-systolic and end-diastolic volumes, ejection fraction, and left ventricular end-diastolic volume, were attenuated by MM07. MM07 also significantly reduced monocrotaline-induced muscularisation of small pulmonary blood vessels. MM07 stimulated endothelial NOS phosphorylation and expression, promoted proliferation, and attenuated apoptosis of human pulmonary arterial endothelial cells in vitro.

CONCLUSION AND IMPLICATIONS

Our findings suggest that chronic treatment with MM07 is beneficial in this animal model of pulmonary arterial hypertension by addressing disease aetiology. These data support the development of G protein-biased apelin receptor agonists with improved pharmacokinetic profiles for use in human disease.

Magnetic resonance fingerprinting with dictionary-based fat and water separation (DBFW MRF): A multi-component approach

Purpose

To obtain a fast and robust fat‐water separation with simultaneous estimation of water T₁, fat T₁, and fat fraction maps.

Methods

We modified an MR fingerprinting (MRF) framework to use a single dictionary combination of a water and fat dictionary. A variable TE acquisition pattern with maximum TE = 4.8 ms was used to increase the fat–water separability. Radiofrequency (RF) spoiling was used to reduce the size of the dictionary by reducing T₂ sensitivity. The technique was compared both in vitro and in vivo to an MRF method that incorporated 3‐point Dixon (DIXON MRF), as well as Cartesian IDEAL with different acquisition parameters.

Results

The proposed dictionary‐based fat–water separation technique (DBFW MRF) successfully provided fat fraction, water, and fat T₁, B₀, and B₁₊ maps both in vitro and in vivo. The fat fraction and water T₁ values obtained with DBFW MRF show excellent agreement with DIXON MRF as well as with the reference values obtained using a Cartesian IDEAL with a long TR (concordance correlation coefficient: 0.97/0.99 for fat fraction–water T₁). Whereas fat fraction values with Cartesian IDEAL were degraded in the presence of T₁ saturation, MRF methods successfully estimated and accounted for T₁ in the fat fraction estimates.

Conclusion

The DBFW MRF technique can successfully provide T₁ and fat fraction quantification in under 20 s per slice, intrinsically correcting T₁ biases typical of fast Dixon techniques. These features could improve the diagnostic quality and use of images in presence of fat.

Spiral MR fingerprinting at 7T with simultaneous B1 estimation

Magnetic resonance fingerprinting is an efficient, new approach for quantitative imaging with MR. We aimed to extend this technique to cases with B1+ inhomogeneities within the imaging volume. Previous approaches have used abrupt changes in flip angles to estimate the B1+ field simultaneously with T1 and T2, using a Cartesian approach in a small-animal scanner at 4.7T. Here, we evaluated whether a similar approach would be suitable for imaging human brains using spiral readouts with a 7T scanner. We found that our modified scheme could significantly reduce the adverse effects of B1+ inhomogeneities even in extreme cases, reducing both the bias and the variance in T2 estimations by an order of magnitude when compared to literature methods. Acquisitions used less than 1.5W/kg SAR and could be performed in 12s per slice. In conclusion, our approach can be used to perform quantitative imaging of the brain at 7T in a short time, simultaneously estimating the B1+ profile.

Elabela/Toddler Is an Endogenous Agonist of the Apelin APJ Receptor in the Adult Cardiovascular System, and Exogenous Administration of the Peptide Compensates for the Downregulation of Its Expression in Pulmonary Arterial Hypertension

Supplemental Digital Content is available in the text.

Background:

Elabela/toddler (ELA) is a critical cardiac developmental peptide that acts through the G-protein–coupled apelin receptor, despite lack of sequence similarity to the established ligand apelin. Our aim was to investigate the receptor pharmacology, expression pattern, and in vivo function of ELA peptides in the adult cardiovascular system, to seek evidence for alteration in pulmonary arterial hypertension (PAH) in which apelin signaling is downregulated, and to demonstrate attenuation of PAH severity with exogenous administration of ELA in a rat model.

Methods:

In silico docking analysis, competition binding experiments, and downstream assays were used to characterize ELA receptor binding in human heart and signaling in cells expressing the apelin receptor. ELA expression in human cardiovascular tissues and plasma was determined using real-time quantitative polymerase chain reaction, dual-labeling immunofluorescent staining, and immunoassays. Acute cardiac effects of ELA-32 and [Pyr¹]apelin-13 were assessed by MRI and cardiac catheterization in anesthetized rats. Cardiopulmonary human and rat tissues from PAH patients and monocrotaline- and Sugen/hypoxia-exposed rats were used to show changes in ELA expression in PAH. The effect of ELA treatment on cardiopulmonary remodeling in PAH was investigated in the monocrotaline rat model.

Results:

ELA competed for binding of apelin in human heart with overlap for the 2 peptides indicated by in silico modeling. ELA activated G-protein– and β-arrestin–dependent pathways. We detected ELA expression in human vascular endothelium and plasma. Comparable to apelin, ELA increased cardiac contractility, ejection fraction, and cardiac output and elicited vasodilatation in rat in vivo. ELA expression was reduced in cardiopulmonary tissues from PAH patients and PAH rat models, respectively. ELA treatment significantly attenuated elevation of right ventricular systolic pressure and right ventricular hypertrophy and pulmonary vascular remodeling in monocrotaline-exposed rats.

Conclusions:

These results show that ELA is an endogenous agonist of the human apelin receptor, exhibits a cardiovascular profile comparable to apelin, and is downregulated in human disease and rodent PAH models, and exogenous peptide can reduce the severity of cardiopulmonary remodeling and function in PAH in rats. This study provides additional proof of principle that an apelin receptor agonist may be of therapeutic use in PAH in humans.

Direct evidence of progressive cardiac dysfunction in a transgenic mouse model of Huntington's disease

HD is a progressive genetic neurological disorder, characterized by motor as well as cognitive impairments. The gene carrying the mutation causing Huntington's disease (HD) is not brain specific, and there is increasing evidence for peripheral, as well as brain pathology in this disorder. Here, we used in vivo and ex vivo techniques to assess the cardiac function of mice transgenic for the HD mutation. Using magnetic resonance imaging (MRI) of the beating heart, we show that abnormalities previously reported in end-stage mice are present by mid-stages of the disease. We also found abnormalities that have not been hitherto reported, including changes in cardiac efficiency and a mechanical distortion of the beating heart. Using the Langendorff preparation, we show reduced coronary blood flow, impaired myocardial contractility and reduced left ventricular developed pressure in HD mouse hearts. Together, our findings suggest that there is significant pathology of the HD mouse heart, even by mid stages of disease. Previous clinical research has demonstrated that the risk of cognitive symptoms increases markedly in patients with heart failure. R6/2 mice show significant progressive cognitive abnormalities, so we hypothesize that cardiac pathology in the R6/2 mouse may contribute, not only to their progressive decline and death, but also to their cognitive dysfunction. We suggest that closer attention should be paid to cardiovascular symptoms in HD patients.

MR fingerprinting with simultaneous B1 estimation

PURPOSE

MR fingerprinting (MRF) can be used for quantitative estimation of physical parameters in MRI. Here, we extend the method to incorporate B1 estimation.

METHODS

The acquisition is based on steady state free precession MR fingerprinting with a Cartesian trajectory. To increase the sensitivity to the B1 profile, abrupt changes in flip angle were introduced in the sequence. Slice profile and B1 effects were included in the dictionary and the results from two- and three-dimensional (3D) acquisitions were compared. Acceleration was demonstrated using retrospective undersampling in the phase encode directions of 3D data exploiting redundancy between MRF frames at the edges of k-space.

RESULTS

Without B1 estimation, T2 and B1 were inaccurate by more than 20%. Abrupt changes in flip angle improved B1 maps. T1 and T2 values obtained with the new MRF methods agree with classical spin echo measurements and are independent of the B1 field profile. When using view sharing reconstruction, results remained accurate (error <10%) when sampling under 10% of k-space from the 3D data.

CONCLUSION

The methods demonstrated here can successfully measure T1, T2, and B1. Errors due to slice profile can be substantially reduced by including its effect in the dictionary or acquiring data in 3D. Magn Reson Med 76:1127-1135, 2016. © 2015 The Authors Magnetic Resonance in Medicine published by Wiley Periodicals, Inc. on behalf of International Society for Magnetic Resonance in Medicine. This is an open access article under the terms of the Creative Commons Attribution-NonCommercial-NoDerivs License, which permits use and distribution in any medium, provided the original work is properly cited, the use is non-commercial and no modifications or adaptations are made.

Combining MRI with PET for partial volume correction improves image-derived input functions in mice

Accurate kinetic modelling using dynamic PET requires knowledge of the tracer concentration in plasma, known as the arterial input function (AIF). AIFs are usually determined by invasive blood sampling, but this is prohibitive in murine studies due to low total blood volumes. As a result of the low spatial resolution of PET, image-derived input functions (IDIFs) must be extracted from left ventricular blood pool (LVBP) ROIs of the mouse heart. This is challenging because of partial volume and spillover effects between the LVBP and myocardium, contaminating IDIFs with tissue signal. We have applied the geometric transfer matrix (GTM) method of partial volume correction (PVC) to 12 mice injected with 18F-FDG affected by a Myocardial Infarction (MI), of which 6 were treated with a drug which reduced infarction size [1]. We utilised high resolution MRI to assist in segmenting mouse hearts into 5 classes: LVBP, infarcted myocardium, healthy myocardium, lungs/body and background. The signal contribution from these 5 classes was convolved with the point spread function (PSF) of the Cambridge split magnet PET scanner and a non-linear fit was performed on the 5 measured signal components. The corrected IDIF was taken as the fitted LVBP component. It was found that the GTM PVC method could recover an IDIF with less contamination from spillover than an IDIF extracted from PET data alone. More realistic values of Ki were achieved using GTM IDIFs, which were shown to be significantly different (p<0.05) between the treated and untreated groups.

Impaired limbic cortico-striatal structure and sustained visual attention in a rodent model of schizophrenia

BACKGROUND

N-methyl-d-aspartate receptor (NMDAR) dysfunction is thought to contribute to the pathophysiology of schizophrenia. Accordingly, NMDAR antagonists such as phencyclidine (PCP) are used widely in experimental animals to model cognitive impairment associated with this disorder. However, it is unclear whether PCP disrupts the structural integrity of brain areas relevant to the profile of cognitive impairment in schizophrenia.

METHODS

Here we used high-resolution magnetic resonance imaging and voxel-based morphometry to investigate structural alterations associated with sub-chronic PCP treatment in rats.

RESULTS

Sub-chronic exposure of rats to PCP (5mg/kg twice daily for 7 days) impaired sustained visual attention on a 5-choice serial reaction time task, notably when the attentional load was increased. In contrast, sub-chronic PCP had no significant effect on the attentional filtering of a pre-pulse auditory stimulus in an acoustic startle paradigm. Voxel-based morphometry revealed significantly reduced grey matter density bilaterally in the hippocampus, anterior cingulate cortex, ventral striatum, and amygdala. PCP-treated rats also exhibited reduced cortical thickness in the insular cortex.

CONCLUSIONS

These findings demonstrate that sub-chronic NMDA receptor antagonism is sufficient to produce highly-localized morphological abnormalities in brain areas implicated in the pathogenesis of schizophrenia. Furthermore, PCP exposure resulted in dissociable impairments in attentional function.

Right ventricular dysfunction in the R6/2 transgenic mouse model of Huntington's disease is unmasked by dobutamine

BACKGROUND

Increasingly, evidence from studies in both animal models and patients suggests that cardiovascular dysfunction is important in HD. Previous studies measuring function of the left ventricle (LV) in the R6/2 model have found a clear cardiac abnormality, albeit with preserved LV systolic function. It was hypothesized that an impairment of RV function might play a role in this condition via mechanisms of ventricular interdependence.

OBJECTIVE

To investigate RV function in the R6/2 mouse model of Huntington's disease (HD).

METHODS

Cardiac cine-magnetic resonance imaging (MRI) was used to determine functional parameters in R6/2 mice. In a first experiment, these parameters were derived longitudinally to determine deterioration of cardiac function with disease progression. A second experiment compared the response to a stress test (using dobutamine) of wildtype and early-symptomatic R6/2 mice.

RESULTS

There was progressive deterioration of RV systolic function with age in R6/2 mice. Furthermore, beta-adrenergic stimulation with dobutamine revealed RV dysfunction in R6/2 mice before any overt symptoms of the disease were apparent.

CONCLUSIONS

This work adds to accumulating evidence of cardiovascular dysfunction in R6/2 mice, describing for the first time the involvement of the right ventricle. Cardiovascular dysfunction should be considered, both when treatment strategies are being designed, and when searching for biomarkers for HD.

Functional assessment of the mouse heart by MRI with a 1-min acquisition

In vivo assessment of heart function in mice is important for basic and translational research in cardiology. MRI is an accurate tool for the investigation of the anatomy and function in the preclinical setting; however, the long scan duration limits its usage. We aimed to reduce the acquisition time of cine MRI to 1 min. We employed spatiotemporal compressed sensing and parallel imaging to accelerate retrospectively gated cine MRI. We compared the functional parameters derived from full and undersampled data in Cartesian and radial MRI by means of Bland-Altman plots. We found that the scan time for the whole heart could be reduced to 2 min with Cartesian sampling and to 1 min with radial sampling. Despite a reduction in the signal-to-noise ratio, the accuracy in the estimation of left and right ventricular volumes was preserved for all tested subjects. This method can be used to perform accurate functional MRI examinations in mice for high-throughput phenotyping or translational studies.

Mitochondria selective S-nitrosation by mitochondria-targeted S-nitrosothiol protects against post-infarct heart failure in mouse hearts

Aims

Recently it has been shown that the mitochondria-targeted S-nitrosothiol MitoSNO protects against acute ischaemia/reperfusion (IR) injury by inhibiting the reactivation of mitochondrial complex I in the first minutes of reperfusion of ischaemic tissue, thereby preventing free radical formation that underlies IR injury. However, it remains unclear how this transient inhibition of mitochondrial complex I-mediated free radicals at reperfusion affects the long-term recovery of the heart following IR injury. Here we determined whether the acute protection by MitoSNO at reperfusion prevented the subsequent development of post-myocardial infarction heart failure.

Methods and results

Mice were subjected to 30 min left coronary artery occlusion followed by reperfusion and recovery over 28 days. MitoSNO (100 ng/kg) was applied 5 min before the onset of reperfusion followed by 20 min infusion (1 ng/kg/min). Infarct size and cardiac function were measured by magnetic resonance imaging (MRI) 24 h after infarction. MitoSNO-treated mice exhibited reduced infarct size and preserved function. In addition, MitoSNO at reperfusion improved outcome measures 28 days post-IR, including preserved systolic function (63.7 ±1.8% LVEF vs. 53.7 ± 2.1% in controls, P = 0.01) and tissue fibrosis.

Conclusions

MitoSNO action acutely at reperfusion reduces infarct size and protects from post-myocardial infarction heart failure. Therefore, targeted inhibition of mitochondrial complex I in the first minutes of reperfusion by MitoSNO is a rational therapeutic strategy for preventing subsequent heart failure in patients undergoing IR injury.

Complex I deficiency due to selective loss of Ndufs4 in the mouse heart results in severe hypertrophic cardiomyopathy

Mitochondrial complex I, the primary entry point for electrons into the mitochondrial respiratory chain, is both critical for aerobic respiration and a major source of reactive oxygen species. In the heart, chronic dysfunction driving cardiomyopathy is frequently associated with decreased complex I activity, from both genetic and environmental causes. To examine the functional relationship between complex I disruption and cardiac dysfunction we used an established mouse model of mild and chronic complex I inhibition through heart-specific Ndufs4 gene ablation. Heart-specific Ndufs4-null mice had a decrease of ∼50% in complex I activity within the heart, and developed severe hypertrophic cardiomyopathy as assessed by magnetic resonance imaging. The decrease in complex I activity, and associated cardiac dysfunction, occurred absent an increase in mitochondrial hydrogen peroxide levels in vivo, accumulation of markers of oxidative damage, induction of apoptosis, or tissue fibrosis. Taken together, these results indicate that diminished complex I activity in the heart alone is sufficient to drive hypertrophic cardiomyopathy independently of alterations in levels of mitochondrial hydrogen peroxide or oxidative damage.

PET/MRI assessment of the infarcted mouse heart

Heart failure originating from myocardial infarction (MI) is a leading cause of death worldwide. Mouse models of ischaemia and reperfusion injury (I/R) are used to study the effects of novel treatment strategies targeting MI, however staging disease and treatment efficacy is a challenge. Damage and recovery can be assessed on the cellular, tissue or whole-organ scale but these are rarely measured in concert. Here, for the first time, we present data showing measures of injury in infarcted mice using complementary techniques for multi-modal characterisation of the heart. We use in vivo magnetic resonance imaging (MRI) to assess heart function with cine-MRI, hindered perfusion with late gadolinium enhancement imaging and muscular function with displacement encoded with stimulated echoes (DENSE) MRI. These measures are followed by positron emission tomography (PET) with 18-F-fluorodeoxyglucose to assess cellular metabolism. We demonstrate a protocol combining each of these measures for the same animal in the same imaging session and compare how the different markers can be used to quantify cardiac recovery on different scales following injury.

MRI and PET in mouse models of myocardial infarction

Myocardial infarction is one of the leading causes of death in the Western world. The similarity of the mouse heart to the human heart has made it an ideal model for testing novel therapeutic strategies. In vivo magnetic resonance imaging (MRI) gives excellent views of the heart noninvasively with clear anatomical detail, which can be used for accurate functional assessment. Contrast agents can provide basic measures of tissue viability but these are nonspecific. Positron emission tomography (PET) is a complementary technique that is highly specific for molecular imaging, but lacks the anatomical detail of MRI. Used together, these techniques offer a sensitive, specific and quantitative tool for the assessment of the heart in disease and recovery following treatment. In this paper we explain how these methods are carried out in mouse models of acute myocardial infarction. The procedures described here were designed for the assessment of putative protective drug treatments. We used MRI to measure systolic function and infarct size with late gadolinium enhancement, and PET with fluorodeoxyglucose (FDG) to assess metabolic function in the infarcted region. The paper focuses on practical aspects such as slice planning, accurate gating, drug delivery, segmentation of images, and multimodal coregistration. The methods presented here achieve good repeatability and accuracy maintaining a high throughput.

PET/MRI in the infarcted mouse heart with the Cambridge split magnet

Chronic heart failure, as a result of acute myocardial infarction, is a leading cause of death worldwide. Combining diagnostic imaging modalities may aid the direct assessment of experimental treatments targeting heart failure in vivo. Here we present preliminary data using the Cambridge combined 18FDG PET/MRI imaging system in a mouse model of acute myocardial infarction. The split-magnet design can deliver uncompromised MRI and PET performance, for better assessment of disease and treatment in a preclinical environment.