0772

Pathological validation of MP-MRI intensity-based signatures in brain cancer patients using autopsy tissue samples
Samuel Bobholz1, Allison Lowman2, Michael Brehler2, Savannah Duenweg1, Fitzgerald Kyereme2, Elizabeth Cochran3, Jennifer Connelly4, Wade Mueller5, Mohit Agarwal2, Darren O'Neill2, Anjishnu Banerjee6, and Peter LaViolette2,7
1Biophysics, Medical College of Wisconsin, Milwaukee, WI, United States, 2Radiology, Medical College of Wisconsin, Milwaukee, WI, United States, 3Pathology, Medical College of Wisconsin, Milwaukee, WI, United States, 4Neurology, Medical College of Wisconsin, Milwaukee, WI, United States, 5Neurosurgery, Medical College of Wisconsin, Milwaukee, WI, United States, 6Biostatistics, Medical College of Wisconsin, Milwaukee, WI, United States, 7Biomedical Engineering, Medical College of Wisconsin, Milwaukee, WI, United States

Synopsis

This study investigated the relationship between MRI-intensity values and pathological features of brain cancer using autopsy tissues as ground truth. Mixed effect models were used to examine the association between T1, T1C, FLAIR, and ADC intensity and pathological features (cellularity, cytoplasm density, and extracellular fluid density), as well as to compare the strength of these associations between GBM and non-GBM patients. These analyses confirmed many of the associations seen in prior literature, but with decreased strength than expected. Additionally, this study found that radio-pathomic associations were weaker in GBM patients than non-GBM patients.

Introduction

MP-MRI is central to the diagnosis and non-invasive monitoring of progression for glioblastomas (GBM) and other brain cancers. T1-weighted images acquired pre- (T1) and post Gd contrast (T1C) are currently used to define the primary tumor area and are often used to define boundaries for surgical resection 1,2. Hyperintensities on FLAIR images are considered signatures of a combination of infiltrative glioma and vasogenic edema resulting from tumor growth3,4, and ADC maps computed from diffusion-weighted images are often used as an inverse correlate of tumor cellularity 5–7. Imaging signatures of brain cancer are often validated using autopsy tissue samples, which allow for pointwise associations between tissue and MRI signals. However, these samples are often taken prior to treatment, and are rarely taken from beyond the traditional contrast-enhancing region. Therefore, further investigation is warranted to assess imaging signatures in the presence of treatment and beyond the contrast-enhancing region, especially in the case of GBM patients, whose tumors present with much greater pathological heterogeneity.
This study used autopsy tissue samples in order to examine the association of T1, T1C, FLAIR, and ADC images with histopathological signatures of brain cancer in both GBM and non-GBM patients. Specifically, we tested the hypothesis that traditional imaging signatures are reflective of cellularity, cytoplasm density (CYT), and extracellular fluid density (ECF), with stronger radio-pathomic relationships observed in non-GBM patients.

Methods

This study included 45 patients (13, non-GBM, 32 GBM) with a pathologically confirmed brain cancer, whose demographic and clinical characteristics are presented in Table 1. A schematic of the data collection process is presented in Figure 2. Across these patients, 93 tissue samples were taken at autopsy, processed, and stained for hematoxylin and eosin (HE). Following staining, tissue samples were digitized using a sliding stage microscope. A color deconvolution-based algorithm was used to extract cellularity, CYT, and ECF across each slide. T1, T1C, FLAIR, and ADC images were collected from each patients’ most recent clinical scan session prior to death. Scans were coregistered to the FLAIR image, and qualitative images (T1, T1C, FLAIR) were intensity normalized by dividing the image by the whole-brain standard deviation. Grey and white matter probability maps were computed using SPM12 segmentations derived from the T1 image. Custom in-house software was used to align tissue samples to the FLAIR image using a non-linear warping computed from manually defined control points 7–9.
Linear mixed effect models were used to assess the relationship between voxel-wise image intensity values and histopathological features. Each model included a covariate for time between clinical imaging and death and grey matter probability, as well as a random effect of subject. Analogous models including an interaction between intensity value and diagnosis were also computed in order to assess the differences in radio-pathomic relationship strength between the two diagnostic groups. Due to the inflated sample size due to voxelwise assessment (n=591,502), p-values were not considered a reliable indicator of significant effects, thus results are presented in terms of measures of effect size.

Results

Fitted model results for cellularity are presented in Figure 3, including standardized model coefficients (i.e. Increase in cells/mm, conditional R2 (cR2, variance explained by the model including random effects), and marginal R2 (mR2, variance explained by only fixed effects). CYT fitted model results are included in Figure 4, and results for ECF are included in Figure 5. Across each plot, associations generally followed expected trends seen in prior literature (i.e. negative ADC-cellularity association, positive FLAIR-ECF association), though models at most accounted for approximately 30% of variance, with fixed effects only accounting for at most 10% of variance. Across the majority of associations, the radio-pathomic association seen within non-GBM patients was stronger than that of the GBM patients.

Discussion

This study provides an autopsy-based pathological validation between image intensity values and histopathological features of brain cancer. These results generally confirmed relationships observed in prior literature, such as ADC as an inverse marker of cellularity and FLAIR as a direct marker of edema (heightened ECF). However, these relationships were much weaker than expected, with the majority of pathological variance left unexplained with regards to each pathological signature. These results generally indicate that traditional imaging signatures may be less robust beyond the non-contrast-enhancing region, and further study may be warranted in order to improve non-invasive signatures of tumor progression. These signatures also may be less strong due to the effect of radiation and other treatments, which influence the relationship between imaging values and tissue pathology 10,11. Additionally, radio-pathomic relationships within GBM patients were weaker than those within non-GBM patients. This indicates that the pathological heterogeneity associated with GBM may confound traditional imaging signatures. While this study represents the largest of its kind, a still relatively small subject-level sample size in relation to the degree of clinical heterogeneity leaves additional research questions. The time between MRI scan time and death may influence results beyond statistical control, and future studies are warranted to investigate the implications of this time period on radio-pathomic relationships.

Conclusion

MP-MRI signatures of cellularity, CYT, and ECF capture a modest degree of variance using autopsy samples as ground truth, with stronger radio-pathomic relationships observed in non-GBM patients than GBM patients.

Acknowledgements

No acknowledgement found.

References

1. Smits M, van den Bent MJ. Imaging Correlates of Adult Glioma Genotypes. Radiology. 2017;284(2):316-331. doi:10.1148/radiol.2017151930

2. Shukla G, Alexander GS, Bakas S, et al. Advanced magnetic resonance imaging in glioblastoma: a review. Chinese Clin Oncol. 2017;6(4):40. doi:10.21037/cco.2017.06.28

3. Husstedt HW, Sickert M, Köstler H, Haubitz B, Becker H. Diagnostic value of the fast-FLAIR sequence in MR imaging of intracranial tumors. Eur Radiol. 2000;10(5):745-752. doi:10.1007/s003300050997

4. Hudak AM, Peng L, Marquez de la Plata C, et al. Cytotoxic and vasogenic cerebral oedema in traumatic brain injury: assessment with FLAIR and DWI imaging. Brain Inj. 2014;28(12):1602-1609. doi:10.3109/02699052.2014.936039

5. Surov A, Meyer HJ, Wienke A. Correlation between apparent diffusion coefficient (ADC) and cellularity is different in several tumors: a meta-analysis. Oncotarget. 2017;8(35):59492-59499. doi:10.18632/oncotarget.17752

6. Chen L, Liu M, Bao J, et al. The correlation between apparent diffusion coefficient and tumor cellularity in patients: a meta-analysis. PLoS One. 2013;8(11):e79008. doi:10.1371/journal.pone.0079008

7. Ellingson BM, Malkin MG, Rand SD, et al. Validation of functional diffusion maps (fDMs) as a biomarker for human glioma cellularity. J Magn Reson Imaging. 2010;31(3):538-548. doi:10.1002/jmri.22068 8. Bobholz SA, Lowman AK, Barrington A, et al. Radiomic Features of Multiparametric MRI Present Stable Associations With Analogous Histological Features in Patients With Brain Cancer. Tomogr (Ann Arbor, Mich). 2020;6(2):160-169. doi:10.18383/j.tom.2019.00029

9. LaViolette PS, Mickevicius NJ, Cochran EJ, et al. Precise ex vivo histological validation of heightened cellularity and diffusion-restricted necrosis in regions of dark apparent diffusion coefficient in 7 cases of high-grade glioma. Neuro Oncol. 2014;16(12):1599-1606. doi:10.1093/neuonc/nou142

10. Hawkins-Daarud A, Rockne RC, Anderson ARA, Swanson KR. Modeling Tumor-Associated Edema in Gliomas during Anti-Angiogenic Therapy and Its Impact on Imageable Tumor. Front Oncol. 2013;3:66. doi:10.3389/fonc.2013.00066

11. Kessler AT, Bhatt AA. Brain tumour post-treatment imaging and treatment-related complications. Insights Imaging. 2018;9(6):1057-1075. doi:10.1007/s13244-018-0661-y

Figures

Demographic and clinical summary of included subjects.

Schematic representation of the study structure including A) data collection for MRI and tissue modalities and associated preprocessing, B) MR-histology coregistration, and C) mixed effect analysis structure

Mixed effect model results for cell density. The left-hand plots represent the first order models, whereas the right hand results represent the interaction between MR intensity and diagnosis, where blue = GBM and red = non-GBM.

Mixed effect model results for cytoplasm density. The left-hand plots represent the first order models, whereas the right hand results represent the interaction between MR intensity and diagnosis, where blue = GBM and red = non-GBM.

Mixed effect model results for ECF density. The left-hand plots represent the first order models, whereas the right hand results represent the interaction between MR intensity and diagnosis, where blue = GBM and red = non-GBM.

Proc. Intl. Soc. Mag. Reson. Med. 29 (2021)
0772