Vessel architecture imaging using multiband gradient-echo/spin-echo EPI

Whole-Brain Vascular Architecture Mapping Identifies Region-Specific Microvascular Profiles in Vivo

BACKGROUND AND PURPOSE

The novel MR imaging technique of vascular architecture mapping allows in vivo characterization of local changes in cerebral microvasculature, but reference ranges for vascular architecture mapping parameters in healthy brain tissue are lacking, limiting its potential applicability as an MR imaging biomarker in clinical practice. We conducted whole-brain vascular architecture mapping in a large cohort to establish vascular architecture mapping parameter references ranges and identify region-specific cortical and subcortical microvascular profiles.

MATERIALS AND METHODS

This was a single-center examination of adult patients with unifocal, stable low-grade gliomas with multiband spin- and gradient-echo EPI sequence at 3T using parallel imaging. Voxelwise plotting of resulting values for gradient-echo (R2*) versus spin-echo (R2) relaxation rates during contrast agent bolus administration generates vessel vortex curves that allow the extraction of vascular architecture mapping parameters representative of, eg, vessel type, vessel radius, or CBV in the underlying voxel. Averaged whole-brain parametric maps were calculated for 9 parameters, and VOI analysis was conducted on the basis of a standardized brain atlas and individual cortical GM and WM segmentation.

RESULTS

Prevalence of vascular risk factors among subjects (n = 106; mean age, 39.2 [SD, 12.5] years; 56 women) was similar to those in the German population. Compared with WM, we found cortical GM to have larger mean vascular calibers (5.80 [SD, 0.59] versus 4.25 [SD, 0.62] P < .001), increased blood volume fraction (20.40 [SD, 4.49] s-1 versus 11.05 [SD, 2.44] s-1; P < .001), and a dominance of venous vessels. Distinct microvascular profiles emerged for cortical GM, where vascular architecture mapping vessel type indicator differed, eg, between the thalamus and cortical GM (mean, -2.47 [SD, 4.02] s-2 versus -5.41 [SD, 2.84] s-2; P < .001). Intraclass correlation coefficient values indicated overall high test-retest reliability for vascular architecture mapping parameter mean values when comparing multiple scans per subject.

CONCLUSIONS

Whole-brain vascular architecture mapping in the adult brain reveals region-specific microvascular profiles. The obtained parameter reference ranges for distinct anatomic and functional brain areas may be used for future vascular architecture mapping studies on cerebrovascular pathologies and might facilitate early discovery of microvascular changes, in, eg, neurodegeneration and neuro-oncology.

Non-contrast free-breathing liver perfusion imaging using velocity selective ASL combined with prospective motion compensation

PURPOSE

To apply velocity selective arterial spin labeling (VSASL) combined with a navigator-based (NAV) prospective motion compensation method for a free-breathing liver perfusion measurement without contrast agent.

METHODS

Sinc-modulated Velocity Selective Inversion (sinc-VSI) pulses were applied as labeling and control pulses. In order to account for respiratory motion, a navigator was employed in the form of a single gradient-echo projection readout, located at the diaphragm along the inferior-superior direction. Prior to each transverse imaging slice of the spin-echo EPI based readouts, navigator and fat suppression were incorporated. Motion data was obtained from the navigator and transmitted back to the sequence, allowing real-time adjustments to slice positioning. The sinc-VSI without velocity-selective gradients during the control condition but with velocity-selective gradients along all three directions during labeling was chosen for the VSASL. The VSASL was compared with pseudo-continuous ASL (pCASL) methods, which selectively tagged the moving spins using a tagging plane placed at the portal vein and hepatic artery.

RESULTS

The motion caused by respiratory activity was effectively computed using the navigator signal. The coefficients of variation (CoV) of average liver voxel in NAV were significantly decreased when compared to breath-hold (BH), with an average reduction of 29.4 ± 18.44% for control images, and 29.89 ± 20.83% for label images (p < 0.001). The resulting maps of normalized ASL signal (normalized to M0) showed significantly higher perfusion weightings in the NAV-compensated VSASL, when compared to the NAV-compensated pCASL techniques.

CONCLUSIONS

This study demonstrates the feasibility of using a navigator-based prospective motion compensation technique in conjunction with VSASL for the measurement of liver perfusion without the use of contrast agents while allowing for free-breathing.

Abnormal vascular structure and function within brain metastases is linked to pembrolizumab resistance

BACKGROUND

We recently conducted a phase 2 trial (NCT028865685) evaluating intracranial efficacy of pembrolizumab for brain metastases (BM) of diverse histologies. Our study met its primary efficacy endpoint and illustrates that pembrolizumab exerts promising activity in a select group of patients with BM. Given the importance of aberrant vasculature in mediating immunosuppression, we explored the relationship between immune checkpoint inhibitor (ICI) efficacy and vascular architecture in the hopes of identifying potential mechanisms of intracranial ICI response or resistance for BM.

METHODS

Using Vessel Architectural Imaging, a histologically validated quantitative metric for in vivo tumor vascular physiology, we analyzed dual-echo DSC/DCE MRI for 44 patients on trial. Tumor and peri-tumor cerebral blood volume/flow, vessel size, arterial and venous dominance, and vascular permeability were measured before and after treatment with pembrolizumab.

RESULTS

BM that progressed on ICI were characterized by a highly aberrant vasculature dominated by large-caliber vessels. In contrast, ICI-responsive BM possessed a more structurally balanced vasculature consisting of both small and large vessels, and there was a trend toward a decrease in under-perfused tissue, suggesting a reversal of the negative effects of hypoxia. In the peri-tumor region, the development of smaller blood vessels, consistent with neo-angiogenesis, was associated with tumor growth before radiographic evidence of contrast enhancement on anatomical MRI.

CONCLUSIONS

This study, one of the largest functional imaging studies for BM, suggests that vascular architecture is linked with ICI efficacy. Studies identifying modulators of vascular architecture, and effects on immune activity, are warranted and may inform future combination treatments.

Simultaneous T1, T2 and T2⁎ mapping of the liver with multi-shot MI-SAGE

PURPOSE

To apply multi-shot high-resolution multi inversion spin and gradient echo (MI-SAGE) acquisition for simultaneous liver T1, T2 and T2* mapping.

METHODS

Inversion prepared spin- and gradient-echo EPI was developed with ascending slice order across measurements for efficient acquisition with T1, T2, and T2⁎ weighting. Multi-shot EPI was also implemented to minimize distortion and blurring while enabling high in-plane resolution. A dictionary-matching approach was used to fit the images to quantitative parameter maps, which were compared to T1 measured by modified Look-Locker (MOLLI), T1 measured by variable flip angle (VFA), T2 measured by multiple echo time-based Half Fourier Single-shot Turbo spin-Echo (HASTE), T2 measured by radial turbo-spin-echo (rTSE) and T2⁎ measured by multiple gradient echo (MGRE) sequences.

RESULTS

The multi-shot variant of the sequence achieved higher in-plane resolution of 1.7 × 1.7 mm2 with good image quality in 28 s. Derived quantitative maps showed comparable values to conventional mapping methods. As measured in phantom and in vivo, MOLLI, MESE and MGRE give closest values to MISAGE. VFA, HASTE and rTSE show obvious overestimation.

CONCLUSIONS

The proposed multi-shot inversion prepared spin- and gradient-echo EPI sequence allows for high-resolution quantitative T1, T2 and T2 liver tissue characterization in a single breath-hold scan.

Structural and functional vascular dysfunction within brain metastases is linked to pembrolizumab inefficacy

Structurally and functionally aberrant vasculature is a hallmark of tumor angiogenesis and treatment resistance. Given the synergistic link between aberrant tumor vasculature and immunosuppression, we analyzed perfusion MRI for 44 patients with brain metastases (BM) undergoing treatment with pembrolizumab. To date, vascular-immune communication, or the relationship between immune checkpoint inhibitor (ICI) efficacy and vascular architecture, has not been well-characterized in human imaging studies. We found that ICI-responsive BM possessed a structurally balanced vascular makeup, which was linked to improved vascular efficiency and an immune-stimulatory microenvironment. In contrast, ICI-resistant BM were characterized by a lack of immune cell infiltration and a highly aberrant vasculature dominated by large-caliber vessels. Peri-tumor region analysis revealed early functional changes predictive of ICI resistance before radiographic evidence on conventional MRI. This study was one of the largest functional imaging studies for BM and establishes a foundation for functional studies that illuminate the mechanisms linking patterns of vascular architecture with immunosuppression, as targeting these aspects of cancer biology may serve as the basis for future combination treatments.

Advanced MR Techniques for Preoperative Glioma Characterization: Part 1

Preoperative clinical magnetic resonance imaging (MRI) protocols for gliomas, brain tumors with dismal outcomes due to their infiltrative properties, still rely on conventional structural MRI, which does not deliver information on tumor genotype and is limited in the delineation of diffuse gliomas. The GliMR COST action wants to raise awareness about the state of the art of advanced MRI techniques in gliomas and their possible clinical translation or lack thereof. This review describes current methods, limits, and applications of advanced MRI for the preoperative assessment of glioma, summarizing the level of clinical validation of different techniques. In this first part, we discuss dynamic susceptibility contrast and dynamic contrast-enhanced MRI, arterial spin labeling, diffusion-weighted MRI, vessel imaging, and magnetic resonance fingerprinting. The second part of this review addresses magnetic resonance spectroscopy, chemical exchange saturation transfer, susceptibility-weighted imaging, MRI-PET, MR elastography, and MR-based radiomics applications. Evidence Level: 3 Technical Efficacy: Stage 2.

Navigator-based slice tracking for prospective motion correction in kidney vessel architecture imaging

OBJECTIVES

To apply a navigator-based slice tracking method to prospectively compensate the respiratory motion for kidney vessel architecture imaging (VAI).

MATERIALS AND METHODS

A dual gradient echo spin echo 2D EPI sequence was developed for kidney VAI. A single gradient-echo slice selection and projection readout at the location of the diaphragm along the inferior-superior direction was applied as a navigator. Navigator acquisition and fat suppression were inserted before each transverse imaging slice. Motion information was calculated after exclusion of the signal saturation in the navigator signal caused by imaging slices. The motion information was then directly sent back to the sequence and slice positioning was adjusted in real-time. The whole sequence was applied during a contrast agent pass-through.

RESULTS

VAI parametric maps show the structural heterogeneity of the renal vasculature. The respiratory motion from the navigator signal was precisely calculated and slice positioning was changed in real-time based on the motion information. The vibration amplitude of the signal intensity of the liver tissue at the liver-lung interface in the case of prospective motion correction (PMC) on is about 28% of the PMC off case. Compared to the case of PMC off, the coefficient of variation was reduced 30% of the case of PMC on.

CONCLUSIONS

This study demonstrates the feasibility of the motion-compensating technique in kidney VAI. The sequence may improve the evaluation of microvasculature in kidney diseases.

Navigator-based slice tracking for kidney pCASL using spin-echo EPI acquisition

PURPOSE

To apply a navigator-based slice-tracking method to prospectively compensate respiratory motion for kidney pseudo-continuous arterial spin labeling (pCASL), using spin-echo (SE) EPI acquisition.

METHODS

A single gradient-echo slice selection and projection readout at the location of the diaphragm along the inferior-superior direction was applied as a navigator. Navigator acquisition and fat suppression were inserted before each transverse imaging slice of the readouts of a 2D-SE-EPI-based pCASL sequence. Motion information was calculated after exclusion of the signal saturation in the navigator signal caused by EPI excitations. The motion information was then used to directly adjust the slice positioning in real time.

RESULTS

The respiratory motion from the navigator signal was calculated, and slice positioning was changed in real time based on the motion information. We could show that motion compensation reduces kidney movement, and that the coefficients of variation across renal perfusion values were significantly reduced when motion correction was applied. The average reduction of coefficients of variation was approximately 20%, resulting in a more accurate and detailed structure of the respective perfusion maps.

CONCLUSIONS

This study demonstrates the feasibility of a navigator-based slice-tracking technique in kidney imaging with a SE-EPI readout pCASL sequence to reduce kidney motion.

Imaging in translational cancer research

This review is aimed at presenting some of the recent developments in translational cancer imaging research, with a focus on novel, recently established, or soon to be established cross-sectional imaging techniques for computed tomography (CT), magnetic resonance imaging (MRI), and positron-emission tomography (PET) imaging, including computational investigations based on machine-learning techniques.

Investigating Recovery After Subarachnoid Hemorrhage With the Imaging, Cognition and Outcome of Neuropsychological Functioning After Subarachnoid Hemorrhage (ICONS) Study: Protocol for a Longitudinal, Prospective Cohort Study

Background

A subarachnoid hemorrhage is a hemorrhage in the subarachnoid space that is often caused by the rupture of an aneurysm. Patients who survive a subarachnoid hemorrhage have a high risk of complications and a negative long-term outcome.

Objective

The aim of the Imaging, Cognition and Outcome of Neuropsychological functioning after Subarachnoid hemorrhage (ICONS) study is to investigate whether and to what extent deficits exist in multiple domains after subarachnoid hemorrhage, including cognition, emotion and behavior, and to investigate whether brain damage can be detected in patients with subarachnoid hemorrhage. We aim to determine which early measures of cognition, emotion and behavior, and brain damage in the subacute stage play a role in long-term recovery after subarachnoid hemorrhage. Recovery is defined as functioning at a societal participation level, with a focus on resuming and maintaining work, leisure activities, and social relationships over the long term.

Methods

The ICONS study is an observational, prospective, single-center cohort study. The study includes patients with subarachnoid hemorrhage admitted to the Neurosurgery Unit of the University Medical Centre Groningen in the Netherlands. The inclusion criteria include diagnosis of an aneurysmal subarachnoid hemorrhage or an angiographically negative subarachnoid hemorrhage, sufficient ability in the Dutch language, and age older than 18 years. Patients will undergo neuropsychological assessment and magnetic resonance imaging 6 months after the subarachnoid hemorrhage. Furthermore, patients will be asked to fill in questionnaires on multiple psychosocial measures and undergo a structured interview at 6 months, 1 year, and 2 years after the subarachnoid hemorrhage. The primary outcome measure of the ICONS study is societal participation 1 year after the subarachnoid hemorrhage, measured with the Dutch version of the Impact on Participation and Autonomy questionnaire.

Results

The study was launched in December 2019 and recruitment is expected to continue until June 2023. At the time of the acceptance of this paper, 76 patients and 69 healthy controls have been included. The first results are expected in early 2023.

Conclusions

The ICONS study is the first to collect and combine data after subarachnoid hemorrhage in a variety of domains, including cognition, emotion and behavior, and brain damage. The results will contribute to a more comprehensive understanding of the consequences of both aneurysmal subarachnoid hemorrhage and angiographically negative subarachnoid hemorrhage, which may ultimately optimize timely treatment for this patient group by setting realistic and attainable goals to improve daily functioning.

Trial Registration

Netherlands Trial Register NL7803; https://trialsearch.who.int/Trial2.aspx?TrialID=NL7803

International Registered Report Identifier (IRRID)

DERR1-10.2196/38190

Brain tumor classification of virtual NMR voxels based on realistic blood vessel-induced spin dephasing using support vector machines

Remodeling of tissue microvasculature commonly promotes neoplastic growth; however, there is no imaging modality in oncology yet that noninvasively quantifies microvascular changes in clinical routine. Although blood capillaries cannot be resolved in typical magnetic resonance imaging (MRI) measurements, their geometry and distribution influence the integral nuclear magnetic resonance (NMR) signal from each macroscopic MRI voxel. We have numerically simulated the expected transverse relaxation in NMR voxels with different dimensions based on the realistic microvasculature in healthy and tumor-bearing mouse brains (U87 and GL261 glioblastoma). The 3D capillary structure in entire, undissected brains was acquired using light sheet fluorescence microscopy to produce large datasets of the highly resolved cerebrovasculature. Using this data, we trained support vector machines to classify virtual NMR voxels with different dimensions based on the simulated spin dephasing accountable to field inhomogeneities caused by the underlying vasculature. In prediction tests with previously blinded virtual voxels from healthy brain tissue and GL261 tumors, stable classification accuracies above 95% were reached. Our results indicate that high classification accuracies can be stably attained with achievable training set sizes and that larger MRI voxels facilitated increasingly successful classifications, even with small training datasets. We were able to prove that, theoretically, the transverse relaxation process can be harnessed to learn endogenous contrasts for single voxel tissue type classifications on tailored MRI acquisitions. If translatable to experimental MRI, this may augment diagnostic imaging in oncology with automated voxel-by-voxel signal interpretation to detect vascular pathologies.

Differentiation of Cerebral Neoplasms with Vessel Size Imaging (VSI)

PURPOSE

Cerebral neoplasms of various histological origins may show comparable appearances on conventional Magnetic Resonance Imaging (MRI). Vessel size imaging (VSI) is an MRI technique that enables noninvasive assessment of microvasculature by providing quantitative estimates of microvessel size and density. In this study, we evaluated the potential of VSI to differentiate between brain tumor types based on their microvascular morphology.

METHODS

Using a clinical 3T MRI scanner, VSI was performed on 25 patients with cerebral neoplasms, 10 with glioblastoma multiforme (GBM), 8 with primary CNS lymphoma (PCNSL) and 7 with cerebral lung cancer metastasis (MLC). Following the postprocessing of VSI maps, mean vessel diameter (vessel size index, vsi) and microvessel density (Q) were compared across tumors, peritumoral areas, and healthy tissues.

RESULTS

The MLC tumors have larger and less dense microvasculature compared to PCNSLs in terms of vsi and Q (p = 0.0004 and p < 0.0001, respectively). GBM tumors have higher yet non-significantly different vsi values than PCNSLs (p = 0.065) and non-significant differences in Q. No statistically significant differences in vsi or Q were present between GBMs and MLCs. GBM tumor volume was positively correlated with vsi (r = 0.502, p = 0.0017) and negatively correlated with Q (r = -0.531, p = 0.0007).

CONCLUSION

Conventional MRI parameters are helpful in differentiating between PCNSLs, GBMs, and MLCs. Additionally incorporating VSI parameters into the diagnostic protocol could help in further differentiating between PCNSLs and metastases and potentially between PCNSLs and GBMs. Future studies in larger patient cohorts are required to establish diagnostic cut-off values for VSI.

Advanced Imaging Techniques for Newly Diagnosed and Recurrent Gliomas

Management of gliomas following initial diagnosis requires thoughtful presurgical planning followed by regular imaging to monitor treatment response and survey for new tumor growth. Traditional MR imaging modalities such as T1 post-contrast and T2-weighted sequences have long been a staple of tumor diagnosis, surgical planning, and post-treatment surveillance. While these sequences remain integral in the management of gliomas, advances in imaging techniques have allowed for a more detailed characterization of tumor characteristics. Advanced MR sequences such as perfusion, diffusion, and susceptibility weighted imaging, as well as PET scans have emerged as valuable tools to inform clinical decision making and provide a non-invasive way to help distinguish between tumor recurrence and pseudoprogression. Furthermore, these advances in imaging have extended to the operating room and assist in making surgical resections safer. Nevertheless, surgery, chemotherapy, and radiation treatment continue to make the interpretation of MR changes difficult for glioma patients. As analytics and machine learning techniques improve, radiomics offers the potential to be more quantitative and personalized in the interpretation of imaging data for gliomas. In this review, we describe the role of these newer imaging modalities during the different stages of management for patients with gliomas, focusing on the pre-operative, post-operative, and surveillance periods. Finally, we discuss radiomics as a means of promoting personalized patient care in the future.

Radiomics and radiogenomics in gliomas: a contemporary update

The natural history and treatment landscape of primary brain tumours are complicated by the varied tumour behaviour of primary or secondary gliomas (high-grade transformation of low-grade lesions), as well as the dilemmas with identification of radiation necrosis, tumour progression, and pseudoprogression on MRI. Radiomics and radiogenomics promise to offer precise diagnosis, predict prognosis, and assess tumour response to modern chemotherapy/immunotherapy and radiation therapy. This is achieved by a triumvirate of morphological, textural, and functional signatures, derived from a high-throughput extraction of quantitative voxel-level MR image metrics. However, the lack of standardisation of acquisition parameters and inconsistent methodology between working groups have made validations unreliable, hence multi-centre studies involving heterogenous study populations are warranted. We elucidate novel radiomic and radiogenomic workflow concepts and state-of-the-art descriptors in sub-visual MR image processing, with relevant literature on applications of such machine learning techniques in glioma management.

Vessel Type Determined by Vessel Architectural Imaging Improves Differentiation between Early Tumor Progression and Pseudoprogression in Glioblastoma

BACKGROUND AND PURPOSE

Currently available perfusion parameters are limited in differentiating early tumor progression and pseudoprogression with no insight about vessel size and type. We aimed to investigate differences in vessel size and type between early tumor progression and pseudoprogression in posttreatment glioblastoma and to demonstrate diagnostic performance using vessel architectural imaging.

MATERIALS AND METHODS

Fifty-eight patients with enlarging contrast-enhancing masses in posttreatment glioblastomas underwent simultaneous gradient recalled-echo and spin-echo dynamic susceptibility contrast imaging. Relative CBV and vessel architectural imaging parameters, including the relative vessel size index, peak shift between gradient recalled echo and spin-echo bolus signal peaks, and arterial dominance scores using spatial dominance of arterial/venous vessel type, were calculated and compared between the 2 conditions. The area under the curve and cross-validation were performed to compare the diagnostic performance of the relative CBV, vessel architectural imaging parameters, and their combinations.

RESULTS

There were 41 patients with early tumor progression and 17 patients with pseudoprogression. Relative to pseudoprogression, early tumor progression showed a lower peak shift (-0.02 versus 0.33, P = .02) and a lower arterial dominance score (1.46 versus 2.11, P = .001), indicating venous dominance. Patients with early tumor progression had higher relative CBV (1.88 versus 1.38, P = .02) and a tendency toward a larger relative vessel size index (99.67 versus 83.17, P = .15) than those with pseudoprogression. Combining arterial dominance scores and relative CBV showed significantly higher diagnostic performance (area under the curve = 0.82; 95% CI, 0.70-0.94; P = .02) than relative CBV alone (area under the curve = 0.64; 95% CI, 0.49-0.79) in distinguishing early tumor progression from pseudoprogression.

CONCLUSIONS

Vessel architectural imaging significantly improved the diagnostic performance of relative CBV by demonstrating venous dominance and a tendency toward larger vessel size in early tumor progression.

Dependence of the frequency distribution around a sphere on the voxel orientation

Microscopically small magnetic field inhomogeneities within an external static magnetic field cause a free induction decay in magnetic resonance imaging that generally exhibits two transverse components that are usually summarized to a complex entity. The Fourier transform of the complex-valued free induction decay is the purely real and positive-valued frequency distribution which allows an easy interpretation of the underlying dephasing mechanism. Typically, the frequency distribution inside a cubic voxel as caused by a spherical magnetic field inhomogeneity is determined by a histogram technique in terms of subdivision of the whole voxel into smaller subvoxels. A faster and more accurate computation is achieved by analytical expressions for the frequency distribution that are derived in this work. In contrast to the usually assumed simplified case of a spherical voxel, we also consider the tilt angles of the cubic voxel to the external magnetic field. The typical asymmetric form of the frequency distribution is reproduced and analyzed for the more realistic case of a cubic voxel. We observe a splitting of frequency distribution peaks for increasing tilt of the cubic voxel against the direction of the external magnetic field in analogy to the case for dephasing around cylindrical, vessel-like objects inside cubic voxels. These results are of value, e.g., for the analysis of susceptibility-weighted images or in quantitative susceptibility imaging since the reconstruction of these images is performed in cubic-shaped voxels.

Gibbs point field model quantifies disorder in microvasculature of U87-glioblastoma

Microvascular proliferation in glioblastoma multiforme is a biological key mechanism to facilitate tumor growth and infiltration and a main target for treatment interventions. The vascular architecture can be obtained by Single Plane Illumination Microscopy (SPIM) to evaluate vascular heterogeneity in tumorous tissue. We make use of the Gibbs point field model to quantify the order of regularity in capillary distributions found in the U87 glioblastoma model in a murine model and to compare tumorous and healthy brain tissue. A single model parameter Γ was assigned that is linked to tissue-specific vascular topology through Monte-Carlo simulations. Distributions of the model parameter Γ differ significantly between glioblastoma tissue with mean 〈Γ_G〉=2.1±0.4, as compared to healthy brain tissue with mean 〈Γ_H〉=4.9±0.4, suggesting that the average Γ-value allows for tissue differentiation. These results may be used for diagnostic magnetic resonance imaging, where it has been shown recently that Γ is linked to tissue-inherent relaxation parameters.