Fast bound pool fraction imaging of the in vivo rat brain: Association with myelin content and validation in the C6 glioma model
Research Highlights
►The bound pool fraction, f, strongly correlates to in vivo myelin density. ►Time-efficient bound pool fraction imaging provides whole-brain parametric f maps. ►Glioma invasion reduces f in gray matter, white matter (WM), and WM fiber tracts.
Introduction
Quantitative magnetization transfer (qMT) extracts the parametric indices that govern magnetization exchange in tissues. In the two-pool model of magnetization transfer (Edzes and Samulski, 1978, Grad et al., 1991, Morrison and Henkelman, 1995), magnetization exchange occurs between protons bound to water (i.e. the free pool) and protons bound to macromolecules (i.e. the bound pool). Mapping of the various tissue-specific parameters that govern magnetization transfer has been accomplished by a variety of techniques (Cercignani et al., 2005, Gloor et al., 2008, Gochberg and Gore, 2003, Helms and Piringer, 2005, Ramani et al., 2002, Ropele et al., 2003, Sled and Pike, 2001, Tozer et al., 2003, Yarnykh, 2002, Yarnykh and Yuan, 2004) that make use of various pulse sequences, sampling schemes and reconstruction algorithms.
A key parameter of qMT captured by all approaches is the molar fraction of protons bound to macromolecules (i.e. bound pool fraction), f. The acquisition of f has been shown to be obtainable with reasonable accuracy regardless of methodology (Portnoy and Stanisz, 2007). In the brain, f has been associated with white matter (WM) fiber tracts (Yarnykh and Yuan, 2004). However, in contrast to quantitative diffusion imaging techniques, f provides data that is not associated with directional coherence of fibers (Underhill et al., 2009). In multiple sclerosis (MS), in vivo alterations to the bound pool fraction have been identified in WM (Davies et al., 2004, Sled and Pike, 2001, Tozer et al., 2003, Yarnykh, 2002). The decrease in bound pool fraction from within MS lesions has been attributed to a localized reduction in the rich macromolecular content of myelin (Fralix et al., 1991, Koenig, 1991, Wolff and Balaban, 1989) associated with demyelination. These in vivo observations have been supported by ex vivo human studies and in vivo animal models. Schmierer et al. (2007) imaged post-mortem ex vivo brains afflicted with MS and found significant differences in the bound pool fraction when measured from histologically identified normal-appearing WM, remyelinated WM lesions, and demyelinated WM lesions. Rausch et al. has produced demyelinating lesions in the rat brain by inducing an experimental autoimmune encephalitis (Rausch et al., 2009). A decrease in the bound pool fraction from regions of interest (ROIs) within demyelinated lesions compared to ROIs from unaffected WM within the same animal and corresponding WM from healthy animals was reported (Rausch et al., 2009). Samsonov et al. (2006, 2010) demonstrated a dramatic and uniform decrease of the bound pool fraction in demyelinated WM using the mutant canine model. While the evidence indicates an association between myelin content and the bound pool fraction, the data (Rausch et al., 2009, Samsonov et al., 2006, Samsonov et al., 2010, Schmierer et al., 2007) were derived from the assessment of abnormal tissues. As such, the dominant underlying tissue property that governs bound pool fraction variations in normal brain tissues remains ambiguous since factors such as inflammation, edema, or genetically-determined structural myelin disorganization may confound observations based on post-mortem diseased brain or animal models.
While the major field of application of qMT methods has been historically associated with demyelinating disorders (Davies et al., 2004, Rausch et al., 2009, Schmierer et al., 2007, Sled and Pike, 2001, Tozer et al., 2003, Yarnykh, 2002), another area of considerable interest is related to brain tumors, especially those of glial origin. The devastating mortality associated with gliomas, particularly high-grade gliomas (5-year survival < 3%) (Laws et al., 2003), has been attributed to tumor cell invasion of normal parenchyma, particularly WM (Giese and Westphal, 2001). Understanding glioma cell invasion is a key aspect necessary to improve survival, which has remained unchanged for more than 30 years (Polin et al., 2005, Scanlon and Taylor, 1979). Accordingly, animal models have been developed to characterize and study in vivo glioma cell invasion (Barth, 1998). Development of a technology that allows non-invasive detection of glioma spread across normal brain tissues would not only enrich animal studies, but may also improve treatment planning and prognostic guidance for clinical decision making (Giese et al., 2003). qMT methods have a promising potential for this purpose due to a high sensitivity of cross-relaxation parameters to WM organization (Underhill et al., 2009, Yarnykh and Yuan, 2004). Current knowledge about cross-relaxation in brain tumors and associated effects on surrounding tissues is very limited. The distinctions in the two-pool model parameters between glioma and normal brain tissues were first demonstrated using a well-established C6 rat glioma model based on continuous wave Z-spectroscopic data acquisition (Quesson et al., 1997). Parametric maps of the bound pool fraction and cross-relaxation rate constant in the single-case human observation (Yarnykh, 2002) identified clear contrast between the glioma tissue, surrounding edema, and normal brain. Recently, Garica et al. (2010) and Tozer et al. (2007) reported preliminary data identifying alterations in qMT parametric maps associated with low-grade and high-grade gliomas, respectively. In the former study, parametric maps appeared heterogenous compared to the relatively homogenous appearance of corresponding contrast-enhanced T1-weighted and T2-weighted images suggesting qMT may provide additional information for diagnostic tumor characterization (Garica et al., 2010). Additionally, the bound pool fraction appeared to differentiate better in gliomas between tumor and perifocal edema compared to other qMT parameters (Garica et al., 2010). These preliminary findings warrant further investigation of CRI in gliomas, and particularly, a detailed study of histological correlations using a well-established animal model.
From a technical standpoint, qMT remains a rather challenging experimental approach due to time-consuming multi-point data acquisition and sophisticated reconstruction algorithms utilizing voxel-based fit of several parameters (from two to five, depending on the model parameterization) (Cercignani et al., 2005, Gloor et al., 2008, Gochberg and Gore, 2003, Helms and Piringer, 2005, Ramani et al., 2002, Ropele et al., 2003, Sled and Pike, 2001, Tozer et al., 2003, Yarnykh, 2002, Yarnykh and Yuan, 2004). This circumstance to a large extent restricts clinical qMT applications. Several studies have been focused on the improvement of time efficiency of qMT techniques (Gloor et al., 2008, Ropele et al., 2003, Yarnykh and Yuan, 2004). Collectively, these approaches are aimed to reduce the parameter space of the two-pool MT model by focusing on the parameters which are assumed to be most clinically informative, while introducing certain model simplifications to avoid estimation of other parameters. In the cross-relaxation imaging (CRI) method (Yarnykh and Yuan, 2004), a reduction of the number of adjustable parameters to two (bound pool fraction, f, and forward cross-relaxation rate constant, k, defined for magnetization transfer from free to bound pool) in conjunction with constraining the remaining parameters (transverse relaxation times of the free and bound pool, T2F and T2B) has been shown to allow the use of four data points with variable off-resonance saturation, and therefore, provide considerable total scan time shortening as compared to multi-point Z-spectroscopic techniques. Ropele et al. (2003) proposed an inversion-based pulse sequence and a simplified theoretical model for determination of the bound pool fraction from two images. Recently, Gloor et al. (2008) introduced the technique based on the on-resonance MT effect induced by the fast SSFP sequence with multi-point data acquisition for reconstruction of the two MT parameters (bound pool fraction and cross-relaxation rate constant). The above techniques still suffer from a relatively long scan time if used for whole-brain imaging (at least 25–30 min for human brain studies) due to the need of either several data points with 3D acquisition (Gloor et al., 2008, Yarnykh and Yuan, 2004) or sequential 2D acquisition with a long repetition time (Ropele et al., 2003). To further improve time efficiency of qMT, a modified CRI technique for rapid acquisition of f maps has been recently proposed in human brain imaging at 3.0 T (Underhill et al., 2010b, Yarnykh et al., 2010). This method applies average-brain constraints for the parameter T2B and the product of the longitudinal relaxation rate and transverse relaxation time of the free pool, R1FT2F, similar to the earlier approach (Underhill et al., 2009, Yarnykh and Yuan, 2004), and additionally utilizes a constrained value for the parameter combination k(1 − f) / f, which physically corresponds to the inverse cross-relaxation rate constant describing magnetization transfer from bound to free pool. This approach (Underhill et al., 2010b, Yarnykh et al., 2010) is particularly fast, since only two off-resonance saturation data points are required, and all images are acquired using a spoiled gradient-echo sequence with a short repetition time. However, the accuracy of this method has not been studied in detail, and its further validation is needed.
In accordance with the above, we hypothesized that the bound pool fraction would be the most robust parameter for relating qMT to histology in the brain. In this study, we first aimed to determine associations between the parameters of the two-pool model (f, k, T2F, and T2B) and key histological characteristics (myelin density, axonal content, and cellularity) in the normal rat brain in vivo at 3.0 T. Second, we sought to identify histopathological correlations of cross-relaxation parameters in the C6 rat glioma model. Then, we aimed to validate with histology a recent time-efficient methodology for single-parameter bound pool fraction mapping in both normal and pathologic brain tissues. Lastly, this single-parameter bound pool fraction mapping method was compared to the magnetization transfer ratio (MTR) mapping, a widely used simple semi-quantitative magnetization transfer technique that describes the relative change in image contrast after a single off-resonance RF saturation pulse (Dousset et al., 1992). MTR has been previously used to study diseases associated with alterations to WM, such as MS (Dousset et al., 1992, Loevner et al., 1995), metastatic lesions (Boorstein et al., 1994), and progressive multifocal leukoencephalopathy (Dousset et al., 1997).
Section snippets
Animal procedures
Nine adult male Wistar rats (Charles River Laboratories, Wilmington, MA) were used in this study—five healthy rats and four rats imaged 2 weeks after intracranial inoculation with C6 cells. Weight range for all animals at the time of imaging was 315–335 g. The C6 cell line was established from a glioma generated by intravenous exposure of random-bred Wistar–Furth rats to N,N′-nitrosomethylurea (Benda et al., 1971), and has been widely used as an animal model of glioma (Bernstein et al., 1990,
CRI versus histology
Representative histology sections from each stain are presented in Fig. 1. Quantitative histology data and two-pool model parameters for each anatomic structure and tumors are reported in Table 1. Fits of sample experimental Z-spectra from GM, WM, and glioma ROIs are presented in Fig. 2. Correlations between two-pool model parameters and quantitative histology are summarized in Table 2, and associations for the bound pool fraction are also exemplified in Fig. 3.
Across all non-tumor anatomic
Discussion
In this study, the bound pool fraction obtainable with CRI in the in vivo normal rat brain at 3.0 T was found to correspond to myelin density on histology across all combined tissue types, and separately within GM and WM. This relationship was further substantiated by whole-brain bound pool fraction maps from animals with gliomas. Alterations in the in vivo bound pool fraction in WM structures consequent of bulk tumor mass effect, tumor cell invasion, or both were readily identifiable. These
Conclusions
The bound pool fraction obtainable with in vivo, whole-brain CRI in the rat at 3.0 T corresponds to myelin density. A time-efficient acquisition directed at solely acquiring the bound pool fraction yields a similar correspondence with myelin density without substantial errors in the estimation of the bound pool fraction induced by the application of a constrained parametric model. The improved time-efficiency may enable improvements in resolution, SNR, or both. The in vivo production of
Acknowledgments
The authors would like to thank Marina S. Ferguson, MT, and Randy Small, HT, for their assistance with tissue preparation and staining.
Grant Support: NIH (R21EB009908 and R21AG029406).
References (63)
- et al.
Three-dimensional quantitative magnetisation transfer imaging of the human brain
NeuroImage
(2005) - et al.
The measurement of cross-relaxation effects in the proton NMR spin-lattice relaxation of water in biological systems: hydrated collagen and muscle
J. Magn. Reson.
(1978) - et al.
Quantitative magnetization transfer by trains of radio frequency pulses in human brain: extension of a free evolution model to continuous-wave-like conditions
Magn. Reson. Imaging
(2005) - et al.
Insights into brain microstructure from the T2 distribution
Magn. Reson. Imaging
(2006) - et al.
Precise estimate of fundamental in-vivo MT parameters in human brain in clinically feasible times
Magn. Reson. Imaging
(2002) - et al.
Direct quantitative comparison between cross-relaxation imaging and diffusion tensor imaging of the human brain at 3.0 T
NeuroImage
(2009) - et al.
Cross-relaxation imaging reveals detailed anatomy of white matter fiber tracts in the human brain
NeuroImage
(2004) Rat brain tumor models in experimental neuro-oncology: the 9 L, C6, T9, F98, RG2 (D74), RT-2 and CNS-1 gliomas
J. Neurooncol.
(1998)- et al.
Morphological and immunochemical studies of rat glial tumors and clonal strains propagated in culture
J. Neurosurg.
(1971) - et al.
C6 glioma cell invasion and migration of rat brain after neural homografting: ultrastructure
Neurosurgery
(1990)
Characterization of the NMR behavior of white matter in bovine brain
Magn. Reson. Med.
Quantitative MRI-pathology correlations of brain white matter lesions developing in a non-human primate model of multiple sclerosis
NMR Biomed.
Metastatic lesions of the brain: imaging with magnetization transfer
Radiology
Invading C6 glioma cells maintaining tumorigenicity
J. Neurosurg.
Estimation of the macromolecular proton fraction and bound pool T2 in multiple sclerosis
Mult. Scler.
Magnetization transfer study of HIV encephalitis and progressive multifocal leukoencephalopathy. Groupe d'Epidemiologie Clinique du SIDA en Aquitaine
AJNR Am. J. Neuroradiol.
Experimental allergic encephalomyelitis and multiple sclerosis: lesion characterization with magnetization transfer imaging
Radiology
Multiexponential T2, magnetization transfer, and quantitative histology in white matter tracts of rat spinal cord
Magn. Reson. Med.
Lipid bilayer and water proton magnetization transfer: effect of cholesterol
Magn. Reson. Med.
Analysis of brain tumors and metastases by quantitative MT imaging with bSSFP: initial experiences
Cost of migration: invasion of malignant gliomas and implications for treatment
J. Clin. Oncol.
Treatment of malignant glioma: a problem beyond the margins of resection
J. Cancer Res. Clin. Oncol.
Quantitative magnetization transfer imaging using balanced SSFP
Magn. Reson. Med.
Quantitative imaging of magnetization transfer using an inversion recovery sequence
Magn. Reson. Med.
Digital Image Processing
Applications of nuclear magnetic cross-relaxation spectroscopy to tissues
Magn. Reson. Med.
Magnetic resonance imaging: physical principles and sequence design
Magnetization transfer imaging in multiple sclerosis
J. Neuroimaging
Cholesterol of myelin is the determinant of gray–white contrast in MRI of brain
Magn. Reson. Med.
Survival following surgery and prognostic factors for recently diagnosed malignant glioma: data from the Glioma Outcomes Project
J. Neurosurg.
Low power method for estimating the magnetization transfer bound-pool macromolecular fraction
J. Magn. Reson. Imaging
Cited by (111)
Brain myelination at 7 months of age predicts later language development
2022, NeuroImageCitation Excerpt :According to multiple histological studies (see Kisel et al. 2022 for a detailed review), MPF provides a strong linear relationship with myelin content in normal and demyelinated brain tissues typically characterized by correlation coefficients greater than 0.8. Practical advantages of MPF include straightforward implementation of fast MPF mapping protocols based on clinical pulse sequences (Korostyshevskaya et al., 2019; Smirnova et al., 2021; Yarnykh et al., 2018b) and insensitivity of MPF to certain confounders, such as iron deposition (Trujillo et al., 2017; Yarnykh et al., 2018a), microstructural organization of white matter fibers (Stikov et al., 2011; Underhill et al., 2009), cellularity (Underhill et al., 2011), and proliferation of microglia (Khodanovich et al., 2018) and astroglia (Khodanovich et al., 2021). At the same time, MPF provides a relative measure of myelin, which potentially can be affected by changes in the content of water (Khodanovich et al., 2018) and non-myelin macromolecules (West et al., 2018).
Robotic Customer Service System ALKETON
2024, Studies in Computational Intelligence