Skip to main content
Advertisement

Main menu

  • Home
  • Content
    • Current Issue
    • Publication Preview--Ahead of Print
    • Past Issue Archive
    • Case of the Week Archive
    • Classic Case Archive
    • Case of the Month Archive
    • COVID-19 Content and Resources
  • For Authors
  • About Us
    • About AJNR
    • Editors
    • American Society of Neuroradiology
  • Submit a Manuscript
  • Podcasts
    • Subscribe on iTunes
    • Subscribe on Stitcher
  • More
    • Subscribers
    • Permissions
    • Advertisers
    • Alerts
    • Feedback
  • Other Publications
    • ajnr

User menu

  • Subscribe
  • Alerts
  • Log in

Search

  • Advanced search
American Journal of Neuroradiology
American Journal of Neuroradiology

American Journal of Neuroradiology

  • Subscribe
  • Alerts
  • Log in

Advanced Search

  • Home
  • Content
    • Current Issue
    • Publication Preview--Ahead of Print
    • Past Issue Archive
    • Case of the Week Archive
    • Classic Case Archive
    • Case of the Month Archive
    • COVID-19 Content and Resources
  • For Authors
  • About Us
    • About AJNR
    • Editors
    • American Society of Neuroradiology
  • Submit a Manuscript
  • Podcasts
    • Subscribe on iTunes
    • Subscribe on Stitcher
  • More
    • Subscribers
    • Permissions
    • Advertisers
    • Alerts
    • Feedback
  • Follow AJNR on Twitter
  • Visit AJNR on Facebook
  • Follow AJNR on Instagram
  • Join AJNR on LinkedIn
  • RSS Feeds
Research ArticleAdult Brain

Hypovascular Cellular Tumor in Primary Central Nervous System Lymphoma is Associated with Treatment Resistance: Tumor Habitat Analysis Using Physiologic MRI

S.Y. Jeong, J.E. Park, N. Kim and H.S. Kim
American Journal of Neuroradiology January 2022, 43 (1) 40-47; DOI: https://doi.org/10.3174/ajnr.A7351
S.Y. Jeong
aFrom the Department of Radiology and Research Institute of Radiology (S.Y.J., J.E.P., H.S.K.), University of Ulsan College of Medicine, Asan Medical Center, Seoul, Republic of Korea
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
  • ORCID record for S.Y. Jeong
J.E. Park
aFrom the Department of Radiology and Research Institute of Radiology (S.Y.J., J.E.P., H.S.K.), University of Ulsan College of Medicine, Asan Medical Center, Seoul, Republic of Korea
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
  • ORCID record for J.E. Park
N. Kim
bJAPEX LLC (N.K.), Seoul, Korea
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
  • ORCID record for N. Kim
H.S. Kim
aFrom the Department of Radiology and Research Institute of Radiology (S.Y.J., J.E.P., H.S.K.), University of Ulsan College of Medicine, Asan Medical Center, Seoul, Republic of Korea
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
  • ORCID record for H.S. Kim
  • Article
  • Figures & Data
  • Supplemental
  • Info & Metrics
  • References
  • PDF
Loading

Abstract

BACKGROUND AND PURPOSE: The microenvironment of lymphomas is known to be highly variable and closely associated with treatment resistance and survival. We tried to develop a physiologic MR imaging–based spatial habitat analysis to identify regions associated with treatment resistance to facilitate the prediction of tumor response after initial chemotherapy in patients with primary central nervous system lymphoma.

MATERIALS AND METHODS: Eighty-one patients with pathologically confirmed primary central nervous system lymphoma were enrolled. Pretreatment physiologic MR imaging was performed, and K-means clustering was used to separate voxels into 3 spatial habitats according to ADC and CBV values. Associations of spatial habitats and clinical and conventional imaging predictors with time to progression were analyzed using Cox proportional hazards modeling. The performance of statistically significant predictors for time to progression was assessed using the concordance probability index.

RESULTS: The 3 spatial habitats of hypervascular cellular tumor, hypovascular cellular tumor, and hypovascular hypocellular tumor were identified. A large hypovascular cellular habitat was most significantly associated with short time to progression (hazard ratio, 2.83; P = . 017). The presence of an atypical finding (hazard ratio, 4.41; P = . 016), high performance score (hazard ratio, 5.82; P = . 04), and high serum lactate dehydrogenase level (hazard ratio, 1.01; P = .013) was significantly associated with time to progression. A predictive model constructed using the habitat score and other imaging parameters showed a concordance probability index for prediction of time to progression of 0.70 (95% CI, 0.54–0.87).

CONCLUSIONS: A hypovascular cellular tumor habitat is associated with treatment resistance in primary central nervous system lymphoma, and its assessment may refine prechemotherapy imaging–based response prediction for patients with primary central nervous system lymphoma.

ABBREVIATIONS:

CE
contrast-enhanced
C-index
concordance probability index
ECOG
Eastern Cooperative Oncology Group
LDH
lactate dehydrogenase
MTX
methotrexate
nCBV
normalized cerebral blood volume
PCNSL
primary central nervous system lymphoma
TTP
time to progression

Primary central nervous system lymphoma (PCNSL) represents a histologically and immunohistochemically homogeneous type of lymphoma.1 However, the microenvironment of B-cell lymphomas was recently shown to be highly variable2 with regard to both the spatial arrangement and composition of cells, including immune and inflammatory cells, blood and lymphatic vascular networks, and the extracellular matrix. The tumor microenvironment is closely associated with treatment resistance and survival2⇓-4 and partly reflects the tumor cell content and genetic aberrations. The mainstay of treatment for PCNSL is high-dose methotrexate (MTX) in combination with other chemotherapeutic agents.5,6 Because drug delivery, distribution, and interaction with the tumor are determined by the tumor microenvironment,7,8 depiction of the spatial arrangement of tumor cellularity and vascularity may provide valuable information for predicting treatment response.

The tumor microenvironment can be depicted using imaging parameters derived from physiologic MR imaging such as ADC mapping, CBV mapping, and permeability parameters.9,10 In PCNSL, high cellularity with a low ADC is associated with poor prognosis,11 whereas a low permeability pattern 7,8 is associated with shorter progression-free survival. However, these previous studies were limited in their depiction of the spatial arrangement, with a single parameter providing unidirectional information on high or low parametric values. Physiologic MR imaging–based spatial habitat analysis refers to the parcellation of voxels on the basis of imaging parameters obtained from physiologic MR imaging of ADC and CBV to explain the habitat-based tumor environment. The parcellation of voxels using physiologic parameters enables depiction of regions of relative hyper- or hypocellularity and relative hyper- or hypovascularity. The combined spatial information from cellularity- and vascularity-reflecting physiologic MR imaging parameters will enable more precise prediction of treatment response,12 especially identification of any tumor fraction with treatment resistance, which may exhibit high cellularity and low vascularity.

Clustering methods can be applied to measurements of multiple parameters to separate different tumor habitats through parcellation.12 We hypothesized that physiologic MR imaging–based tumor parcellation could potentially provide spatial information on pretreatment PCNSL and identify regions of treatment resistance in patients with PCNSL. The purpose of this study was to develop spatial habitat analysis based on multiparametric physiologic MR imaging and identify regions associated with treatment resistance to facilitate the prediction of tumor response after initial chemotherapy in patients with PCNSL.

MATERIALS AND METHODS

Study Population

This retrospective study was reviewed and approved by the institutional review board of Asan Medical Center (institutional review board No. 2021–0265), and the requirement for informed consent was waived. A review of the database of Asan Medical Center identified 103 patients who were diagnosed with PCNSL between July 2011 and September 2018. The inclusion criteria were as follows: 1) pretreatment acquisition of 3T MR imaging with a brain tumor protocol including anatomic MR imaging and physiologic MR imaging (DWI and DSC perfusion imaging); 2) histopathologic confirmation of PCNSL by stereotactic biopsy; and 3) adequate follow-up examinations (no longer than 2 weeks after completion of all planned therapy) to determine the treatment response according to clinicoradiologic consensus after MTX-based chemotherapy. Patients were excluded if they were lost to follow-up visits (n = 6), the tumor was surgically resected (n = 15), or any imaging was unreadable (because of artifacts, n = 1). Therefore, 81 patients were included in this study (age range, 29–80 years; mean age, 61.6 [SD, 11.8] years; 43 women). A flow chart of patient selection is shown in Fig 1.

FIG 1.
  • Download figure
  • Open in new tab
  • Download powerpoint
FIG 1.

Flow chart of patient recruitment.

Clinical Variables and Treatment

Clinical and pathologic data were retrieved from the medical records. Baseline characteristics including age, sex, Eastern Cooperative Oncology Group (ECOG) performance status at diagnosis,13 initial tumor volume, mean serum lactate dehydrogenase (LDH) level, mean CSF total protein, pretreatment imaging characteristics (location and presence of atypical findings), and initial posttreatment response based on a clinicoradiologic examination were collected.

All included patients had a pathologic diagnosis of diffuse large B-cell PCNSL and received MTX-based induction chemotherapy following 1 of 3 protocols: 1) MTX only (3–8 g/m2/day; n = 54); 2) a combination therapy with MTX and procarbazine, lomustine (CCNU), and vincristine regimen (PCV): procarbazine (0.1–0.2 g/m2/day)/vincristine (1.4 mg/m2/week; n = 15); or 3) a combination therapy with MTX, PCV, and rituximab (375 mg/m2/week; n = 12). After the patient was diagnosed with diffuse large B-cell PCNSL, dexamethasone was used in all patients for supportive care to minimize symptoms due to vasogenic edema or tumor mass effect.

Response Assessment and Outcome Definition

The posttreatment tumor response was assessed according to the International Primary CNS Lymphoma Collaborative Group criteria14 based on MR imaging, corticosteroid use, CSF cytology, and slit lamp examination. A complete response was defined as resolution of contrast-enhancing lesions on follow-up MR imaging and, if indicated, by CSF cytologic analysis (performed if CSF cytology was positive for malignant cells at the time of diagnostic staging). A partial response was defined as an interval decrease in contrast-enhancing lesion volume, and progressive disease was defined as an interval increase in contrast-enhancing lesion volume or development of new enhancing lesions on follow-up contrast-enhanced MR imaging or involvement of the eye or CSF. All follow-up imaging to assess the overall treatment response was performed within 2 weeks after completion of all planned chemotherapy of the first or second cycle. Clinicoradiologic diagnoses of posttreatment responses were made by consensus between a neurooncologist (Shin-kyo Yoon, with 10 years of experience in neuro-oncology practice) and a neuroradiologist (Ho Sung Kim, with 18 years of experience in neuro-oncologic imaging) after complete imaging and medical chart review.15 Both complete response and partial response were used as a composite end point to define “treatment response,” and both stable disease and progressive disease were used as a composite end point to define “treatment resistance.” Time to progression (TTP) was defined5 as the interval between the initiation of induction therapy until documentation of resistance to induction therapy (MR imaging after 6 cycles or stable disease after at least 2 cycles and progression or death from any cause), relapse, or death from any cause after having achieved a partial or complete response.

MR Imaging Acquisitions

The brain tumor imaging protocol was acquired on a 3T scanner (Ingenia 3.0 CX; Philips Healthcare) and included conventional and advanced sequences, including T2WI, T2-weighted FLAIR, T1WI, DWI, DSC perfusion imaging, and contrast-enhanced (CE)-T1WI. The detailed imaging protocols are provided in the Online Supplemental Data, and the preprocessing procedures for the DWI and DSC imaging are summarized in the Online Supplemental Data.

Conventional Image Analysis

All MR images were read by an experienced neuroradiologist (S.Y. Jeong with 6 years of experience) who was blinded to the clinical outcomes. Conventional MR imaging parameters included deep or superficial location and the presence of atypical features of necrosis or a hemorrhagic component.16,17 A deep location meant involvement of the basal ganglia, corpus callosum, periventricular areas, brainstem, and cerebellum. A superficial location meant involvement of the cerebral hemisphere with the exception of deep brain structures. The Online Supplemental Data show cases demonstrating atypical findings of PCNSL on conventional imaging.

Deep Learning Segmentation

The overall process of the deep learning segmentation and tumor habitat analysis is shown in Fig 2. The 3D CE-T1WI of each patient was registered and resampled into 1 × 1 × 1 mm isometric voxel sizes. T2WI, ADC, and normalized CBV (nCBV) maps were then coregistered and resampled to the isovoxel 3D enhanced T1WI using rigid transformations with 6 df in the SPM package (Version 12; http://www.fil.ion.ucl.ac.uk/spm/).

FIG 2.
  • Download figure
  • Open in new tab
  • Download powerpoint
FIG 2.

The overall process of the deep learning segmentation and tumor habitat analysis. A, Image acquisition, registration, and deep learning segmentation for a contrast-enhancing lesion (CEL). B, Extraction of ADC and nCBV values from the CEL and voxel classifications based on ADC and nCBV values. The individual voxels in each cluster are grouped according to their similarities and differences using a K-means clustering algorithm. C, The cluster number is set to 3 to depict 3 different habitats according to the combinations of ADC and nCBV parameters: hypervascular cellular, hypovascular cellular, and hypovascular hypocellular. D, Voxels are shown as spatial habitats in the original image space. Associations of pretreatment tumor habitats with TTP were analyzed.

Deep learning segmentation was performed using an implementation of a 3D UNet-based method (https://github.com/MIC-DKFZ/nnUNet)18 applied to the FLAIR and 3D CE-T1WI using the PyTorch package, Version 1.1 in Python 3.7 (http://www.python.org). This deep learning segmentation package is optimized for processing heterogeneous MR imaging results with various pathologic characteristics. The segmentation process included identification of contrast-enhancing lesions, necrosis, and peritumoral edema, with the contrast-enhancing lesions being chosen for subsequent investigation.

Tumor Habitat Analysis

The final voxel classifications based on ADC and nCBV values were implemented using a K-means clustering module in the scikit-learn Python package (https://scikit-learn.org/stable/index.html). The individual voxels in each cluster were grouped according to their similarities and differences using a K-means clustering algorithm with squared Euclidean distances between voxel intensities as the similarity metric. The cluster number was set to 3 to depict 3 different habitats according to the combinations of ADC and nCBV parameters, with the intention of avoiding overparameterized models19 and producing a model that was easily explainable. For example, low or high ADC and CBV values reflected “hyper-/hypovascularity” and “cellularity/hypocellularity,” respectively, and 3 clusters was the lowest number to demonstrate the biologic differences.

Population-Level Clustering.

Three clusters were set using ADC and nCBV values: Cluster 1 represented hypervascular cellular tumor with high CBV values and low ADC values; cluster 2 represented hypovascular cellular tumor with low CBV values and low ADC values; and cluster 3 represented hypovascular hypocellular tumor with low CBV values and high ADC values. All voxels were allocated to 1 of the 3 clusters and were displayed as spatial habitats in the original image space. The parameter value ranges for the boundary of the spatial habitats were 4.37–4.9×10−6 mm2/s for nCBV and 103–121×10−6 mm2/s for ADC.

Statistical Analysis

All results are reported as median with range or 95% CI for continuous variables and as frequency or percentage for categoric variables. All statistical tests were 2-sided, and a P value < .05 was considered to indicate statistical significance. Statistical analyses were performed by using R software, Version 3.4.3 (https://www.r-project.org).

The associations between TTP and clinical and conventional imaging predictors were determined by univariable Cox analysis. Univariate analyses were performed to analyze the associations of tumor habitats with TTP using Cox regression or the Kaplan-Meier method (log-rank test). Hazard ratios indicate the relative change in hazard incurred by a unit increase in each parameter, and 20,000 voxels were considered as 1 U in this study.

The Pearson correlation was used to further analyze the correlation of each tumor habitat with the overall tumor burden of the contrast-enhancing lesion.

To compare clinical predictors, conventional imaging findings, and tumor habitats, random survival forest analysis was performed. The variable importance of predictors for TTP was calculated with 10-fold cross-validation using the randomForestSRC (https://cran.r-project.org/web/packages/randomForestSRC/index.html) module in R.

Habitat Risk Score.

For the significant predictors identified in the univariate Cox regression, an optimal cutoff for stratifying high- and low-risk groups was estimated using the maxstat algorithm in R with 10-fold cross-validation, which ensured unbiased prediction within the sample.20 Using this cutoff from the exploratory analysis, we developed a habitat risk score to stratify patients according to risk, with each predictor being assigned a discrete score of 1 if it was above its cutoff and zero if it was below it.

Prediction Performance for TTP.

A model combining the habitat risk score and clinical and conventional imaging predictors was created, and its concordance probability index (C-index) was calculated. The performance of the model was quantified with respect to discrimination and calibration.21 Discrimination was measured with the Harrell C-index, and calibration was tested using the version of the Hosmer-Lemeshow test by D’Agostino and Nam.22

RESULTS

Patient and Imaging Characteristics of the Study Population

The baseline demographics, clinical characteristics, and conventional imaging characteristics of the study patients and tumors are summarized in Table 1. All the study patients underwent initial treatment with MTX-based chemotherapy. The mean duration between pretreatment MR imaging and the initial chemotherapy was 4.8 (SD, 5.1 ) days. The mean follow-up duration between pretreatment MR imaging and the first posttreatment MR imaging was 50.8 (SD, 27.1) days. No patient died from treatment-related toxicity.

View this table:
  • View inline
  • View popup
Table 1:

Patients and imaging characteristics of the PCNSLs included in this study

Clinical and Conventional Imaging Predictors for TTP

The results of the univariate analysis to evaluate clinical and conventional imaging predictors are summarized in Table 2. Among the potential clinical and imaging predictors, the presence of an atypical finding (hazard ratio, 4.41; 95% CI, 1.32–14.71; P = .016), a high ECOG score of 3 or 4 (hazard ratio, 5.82; 95% CI, 1.07–31.66; P = .04), and a high serum LDH level (hazard ratio, 1.01; 95% CI, 1.00–1.02; P = .013) were significantly associated with TTP in PCNSL.

View this table:
  • View inline
  • View popup
Table 2:

Clinical and imaging predictors associated with TTP in PCNSL

Tumor Habitats Associated with TTP

The results of the univariate analysis to evaluate the associations of the 3 tumor habitats with TTP are summarized in Table 3. Among the tumor habitats, a large hypovascular cellular habitat (20,000 voxels) was significantly associated with a short TTP (hazard ratio, 2.83; 95% CI, 1.20–6.65; P = .017; Fig 3). No other tumor habitats or percentages of the tumor habitats were associated with TTP.

FIG 3.
  • Download figure
  • Open in new tab
  • Download powerpoint
FIG 3.

A, Demonstration of the 3 spatial habitats defined by clustering of voxels using normalized ADC and nCBV maps in a 53-year-old male patient. The hypervascular cellular habitat (red) shows high nCBV and low ADC, the hypovascular cellular habitat (green) shows low nCBV and low ADC, and the hypovascular hypocellular habitat (blue) shows low nCBV and high ADC. The tumor exhibits a large hypovascular cellular habitat (green), and a persistent enhancing mass was associated with a short TTP after initial chemotherapy. B, Demonstration of the 3 spatial habitats defined by clustering of voxels using normalized ADC and nCBV maps in a 57-year-old male patient. The hypervascular cellular habitat (red) shows high nCBV and low ADC, the hypovascular cellular habitat (green) shows low nCBV and low ADC, and the hypovascular hypocellular habitat (blue) shows low nCBV and high ADC. The tumor has a small hypovascular cellular habitat (green) and showed a complete response at 53 days after initial chemotherapy.

View this table:
  • View inline
  • View popup
Table 3:

Exploratory analysis of spatial habitats for predicting TTP in patients with PCNSL

The single physiologic parameters of ADC and CBV showed no significant association with TTP (P > .05).

The correlations between the overall tumor burden and each habitat were moderate (hypervascular cellular habitat, r = 0.51, P < .001; and hypovascular hypocellular habitat, r = 0.55, P < .001) to strong (hypovascular cellular habitat, r = 0.81; P < .001).

The variable importance values of the clinical predictors, conventional imaging findings, and tumor habitats are shown in the Online Supplemental Data. The variable importance of tumor habitats of the hypovascular cellular habitat (variable importance, 0.0049) was higher than that of age (variable importance, 0.0037) and lower than that of atypical findings (variable importance, 0.0185).

Risk Stratification Using the Tumor Habitat Score

The optimal tumor habitat cutoff for stratifying patients according to treatment response and treatment resistance was >9000 voxels in the hypovascular cellular habitat. This cutoff value stratified patients according to treatment response and treatment resistance with a significant difference in the log-rank test (P = .001). Figure 4 shows the Kaplan-Meier survival curve and the risk table generated on the basis of the hypovascular cellular habitat. A habitat risk score based on discrete increments in the hypovascular cellular habitat was assigned for the patients.

FIG 4.
  • Download figure
  • Open in new tab
  • Download powerpoint
FIG 4.

Kaplan-Meier analysis of TTP in patients with PCNSL stratified by hypovascular cellular habitat (log-rank, P = .01).

Performance of the Prediction Model for TTP

The results of the prediction of TTP in PCNSL using tumor habitat and single physiologic MR imaging parameters are summarized in Table 4. The C-index for prediction of TTP using tumor habitats was 0.65 (95% CI, 0.52–0.78), which was significantly higher than that of the single physiologic MR imaging parameters.

Combination of Clinical, Conventional Imaging Predictors, and Tumor HabitatsCombined Clinical and Conventional PredictorsTumor HabitatsClinical Predictors OnlyConventional Imaging Predictors Only
C-index0.730.710.650.680.63
95% CI0.67–0.800.54–0.780.52–0.780.54–0.820.50–0.76
P valueReference.01.012.81.62
  • ↵a Combined clinical predictors were age, ECOG score, and mean serum LDH level; the conventional imaging predictor was the presence of atypical image findings. P value refers to the significance in the difference of the C-indices between the combined model and the single model assessed using “CompareC” (https://cran.r-project.org/web/packages/compareC/index.html) in the R statistical and computing software.

Table 4:

Prediction of TTP in PCNSLs according to tumor habitat, ADC, and CBVa

A predictive model was constructed using the habitat score, the presence of an atypical finding, the ECOG score, and the serum LDH level. The C-index of this model combining selected clinical predictors, conventional imaging predictors, and tumor habitat for predicting TTP was 0.73 (95% CI, 0.67–0.80), which was higher than that of models based on clinical predictors only (C-index, 0.68; 95% CI, 0.54–0.82) or conventional imaging predictors only (C-index, 0.63; 95% CI, 0.50–0.76). The Hosmer-Lemeshow test indicated that the model had good calibration (P > .5). The C-index of the combined model of clinical predictors and conventional imaging predictors for predicting TTP was 0.71 (95% CI, 0.54–0.78), which was slightly lower than the combined habitat score model, but it did not reach statistical significance (P = .10).

DISCUSSION

In this study, we demonstrated that in PCNSL, the TTP after initial chemotherapy could be predicted by physiologic MR imaging–based tumor habitat analysis. A hypovascular cellular habitat with both low CBV and ADC values was significantly associated with treatment resistance, with a large extent of hypovascular cellular habitat being associated with a short TTP. Furthermore, a model combining selected clinical parameters, conventional imaging predictors, and tumor habitat analysis had great predictive power for identifying treatment resistance to initial chemotherapy.

The tumor microenvironment is composed of both cellular and molecular components that continuously interact. The tumor microenvironments in lymphomas are highly variable with regard to both the spatial arrangement and composition of cells. The tumor microenvironment plays critical roles in the regulation of tumor cell survival, proliferation, and immune escape, which leads to treatment resistance.2,4,8 Imaging techniques such as texture analysis, fractal analysis, and habitat analysis can quantify these spatial heterogeneities of microenvironments within tumors,12 and in this study, we used tumor habitat analysis to identify tumor subregions and then used them to predict treatment resistance. Tumor habitat analysis, which parcellates similar voxels using a clustering method, can reflect the heterogeneity of tumor cellularity and vascularity and identify subregions that can be used to predict the response to chemotherapy in PCNSL.

A significant inverse correlation exists between tumor cellularity and ADC values,23 suggesting that ADC values can be used as a surrogate marker of tumor proliferation.24⇓-26 ADC values have shown potential as a prognostic imaging biomarker for the response to treatment of non-Hodgkin lymphoma.27 In our study, low ADC values in PCNSL with a cellular habitat indicating a high tumor burden correlated with adverse clinical outcomes.

CBV measurements obtained using DSC MR imaging are one of the indirect modalities used to evaluate the functional vascularity of the tumor and its oxygenation status,28,29 which are closely related to the response to chemotherapy or radiation therapy.30 The hypovascular habitat identified in our study may reflect a hypoxic microenvironment that is resistant to treatment, with changes in gene and molecular expression and evolution toward increased malignancy and an increasingly aggressive phenotype.30 Furthermore, low CBV values may represent a decrease in the patent vessels that deliver the chemotherapeutic agent to the tumor bed, which, therefore, induces resistance to treatment. This is supported in a previous study31 that found that high CBV values were associated with a good prognosis with longer overall survival and progression-free survival in PCNSL. The hypovascular cellular habitat showed the strongest correlation with overall tumor burden, and it may indicate the region of tumor core to be treated rather than the region responding well to therapy.

Five clinical parameters are well-known to correlate with a worse prognosis in PCNSL: elevated LDH, age older than 60 years, poor performance status, elevated CSF protein, and location of tumor within the deep regions of the brain such as periventricular regions, basal ganglia, brainstem, and/or cerebellum.16 The presence of atypical imaging findings that reflect necrosis or a hemorrhagic component is associated with short TTP.16,17 These findings also indicate poor tumor perfusion and decreased perfusion-related drug delivery because of necrosis and hemorrhage within the tumor.32,33

This study has several limitations. The first is the small sample size resulting from the rarity of the disease and the lack of external validation. Second, the chemotherapy regimens used in this study were not unified. Because there is no consensus regarding the other agents to be used in combination with MTX, several combinations of MTX-based chemotherapy were tried. We did not perform subgroup analysis per regimen because there were only a small number of patients in each combination regimen. Also, although high-dose radiation therapy plus consolidation chemotherapy or radiation therapy or both are used, there is no consensus on their combination with MTX.34 Our results should be confirmed by other studies with strict control of the chemotherapy regimen. Finally, the prediction performance was not high, with a C-index of only 0.70. Nonetheless, our study may be meaningful because it demonstrates the heterogeneity of PCNSL using physiologic MR imaging and the possibility of a future imaging biomarker–directed therapeutic plan. A well-designed prospective study is warranted.

CONCLUSIONS

The analysis of spatial habitats derived from multiparametric physiologic MR imaging may provide added value for predicting treatment response after initial chemotherapy. A hypovascular cellular tumor habitat was associated with treatment resistance in PCNSL and may be useful for refining imaging-based prediction for patients with PCNSL before they undergo chemotherapy.

ACKNOWLEDGMENT

We would like to thank Shin-Kyo Yoon, MD, for helping us in assessment of treatment response in this study.

Footnotes

  • This work was supported by the Basic Science Research Program through the National Research Foundation of Korea (NRF) funded by the Ministry of Education, Science and Technology (grant number: NRF-2020R1A2B5B01001707 and NRF-2020R1A2C4001748) and by a grant (2020IP0097) from the Asan Institute for Life Sciences, Asan Medical Center, Seoul, Korea.

  • Disclosure forms provided by the authors are available with the full text and PDF of this article at www.ajnr.org.

References

  1. 1.↵
    1. Cai Q,
    2. Fang Y,
    3. Young KH
    . Primary central nervous system lymphoma: molecular pathogenesis and advances in treatment. Transl Oncol 2019;12:523–38 doi:10.1016/j.tranon.2018.11.011 pmid:30616219
    CrossRefPubMed
  2. 2.↵
    1. Fowler NH,
    2. Cheah CY,
    3. Gascoyne RD, et al
    . Role of the tumor microenvironment in mature B-cell lymphoid malignancies. Haematologica 2016;101:531–40 doi:10.3324/haematol.2015.139493 pmid:27132279
    Abstract/FREE Full Text
  3. 3.↵
    1. Alame M,
    2. Pirel M,
    3. Costes-Martineau V, et al
    . Characterisation of tumour microenvironment and immune checkpoints in primary central nervous system diffuse large B cell lymphomas. Virchows Arch 2020;476:891–902 doi:10.1007/s00428-019-02695-6 pmid:31811434
    CrossRefPubMed
  4. 4.↵
    1. Scott DW,
    2. Gascoyne RD
    . The tumour microenvironment in B cell lymphomas. Nat Rev Cancer 2014;14:517–34 doi:10.1038/nrc3774 pmid:25008267
    CrossRefPubMed
  5. 5.↵
    1. Lenz G,
    2. Dreyling M,
    3. Hoster E, et al
    . Immunochemotherapy with rituximab and cyclophosphamide, doxorubicin, vincristine, and prednisone significantly improves response and time to treatment failure, but not long-term outcome in patients with previously untreated mantle cell lymphoma: results of a prospective randomized trial of the German Low Grade Lymphoma Study Group (GLSG). J Clin Oncol 2005;23:1984–92 doi:10.1200/JCO.2005.08.133 pmid:15668467
    Abstract/FREE Full Text
  6. 6.↵
    1. Nayak L,
    2. Batchelor TT
    . Recent advances in treatment of primary central nervous system lymphoma. Curr Treat Options Oncol 2013;14:539–52 doi:10.1007/s11864-013-0252-6 pmid:23934511
    CrossRefPubMed
  7. 7.↵
    1. Chung SR,
    2. Choi YJ,
    3. Kim HS, et al
    . Tumor vascular permeability pattern is associated with complete response in immunocompetent patients with newly diagnosed primary central nervous system lymphoma: retrospective cohort study. Medicine 2016;95:e2624 doi:10.1097/MD.0000000000002624 pmid:26871782
    CrossRefPubMed
  8. 8.↵
    1. Hatzoglou V,
    2. Oh JH,
    3. Buck O, et al
    . Pretreatment dynamic contrast-enhanced MRI biomarkers correlate with progression-free survival in primary central nervous system lymphoma. J Neurooncol 2018;140:351–58 doi:10.1007/s11060-018-2960-z pmid:30073640
    CrossRefPubMed
  9. 9.↵
    1. Kim M,
    2. Kim HS
    . Emerging techniques in brain tumor imaging: what radiologists need to know. Korean J Radiol 2016;17:598–619 doi:10.3348/kjr.2016.17.5.598 pmid:27587949
    CrossRefPubMed
  10. 10.↵
    1. Provenzale JM,
    2. Mukundan S,
    3. Barboriak DP
    . Diffusion-weighted and perfusion MR imaging for brain tumor characterization and assessment of treatment response. Radiology 2006;239:632–49 doi:10.1148/radiol.2393042031 pmid:16714455
    CrossRefPubMed
  11. 11.↵
    1. Barajas RF Jr.,
    2. Rubenstein JL,
    3. Chang JS, et al
    . Diffusion-weighted MR imaging derived apparent diffusion coefficient is predictive of clinical outcome in primary central nervous system lymphoma. AJNR Am J Neuroradiol 2010;31:60–66 doi:10.3174/ajnr.A1750 pmid:19729544
    Abstract/FREE Full Text
  12. 12.↵
    1. O'Connor JP,
    2. Rose CJ,
    3. Waterton JC, et al
    . Imaging intratumor heterogeneity: role in therapy response, resistance, and clinical outcome. Clin Cancer Res 2015;21:249–57 doi:10.1158/1078-0432.CCR-14-0990 pmid:25421725
    Abstract/FREE Full Text
  13. 13.↵
    1. Oken MM,
    2. Creech RH,
    3. Tormey DC, et al
    . Toxicity and response criteria of the Eastern Cooperative Oncology Group. Am J Clin Oncol 1982;5:649–55 doi:10.1097/00000421-198212000-00014 pmid:7165009
    CrossRefPubMed
  14. 14.↵
    1. Abrey LE,
    2. Batchelor TT,
    3. Ferreri AJ, et al
    ; International Primary CNS Lymphoma Collaborative Group. Report of an international workshop to standardize baseline evaluation and response criteria for primary CNS lymphoma. J Clin Oncol 2005;23:5034–43 doi:10.1200/JCO.2005.13.524 pmid:15955902
    Abstract/FREE Full Text
  15. 15.↵
    1. Wen PY,
    2. Macdonald DR,
    3. Reardon DA, et al
    . Updated response assessment criteria for high-grade gliomas: Response Assessment in Neuro-Oncology Working Group. J Clin Oncol 2010;28:1963–72 doi:10.1200/JCO.2009.26.3541 pmid:20231676
    Abstract/FREE Full Text
  16. 16.↵
    1. Ferreri AJ,
    2. Blay JY,
    3. Reni M, et al
    . Prognostic scoring system for primary CNS lymphomas: the International Extranodal Lymphoma Study Group experience. J Clin Oncol 2003;21:266–72 doi:10.1200/JCO.2003.09.139 pmid:12525518
    Abstract/FREE Full Text
  17. 17.↵
    1. Thurnher MM,
    2. Rieger A,
    3. Kleibl-Popov C, et al
    . Primary central nervous system lymphoma in AIDS: a wider spectrum of CT and MRI findings. Neuroradiology 2001;43:29–35 doi:10.1007/s002340000480 pmid:11214644
    CrossRefPubMed
  18. 18.↵
    1. Isensee F,
    2. Jaeger PF,
    3. Kohl SA, et al
    . nnU-Net: a self-configuring method for deep learning-based biomedical image segmentation. Nat Methods 2021;18:203–11 doi:10.1038/s41592-020-01008-z pmid:33288961
    CrossRefPubMed
  19. 19.↵
    1. Erickson GJ,
    2. Smith CR
    1. Gull SF
    . Bayesian inductive inference and maximum entropy. In: Erickson GJ, Smith CR, eds. Maximum-Entropy and Bayesian Methods in Science and Engineering: Fundamental Theories of Physics (An International Book Series on The Fundamental Theories of Physics: Their Clarification, Development and Application), vol 31-32. Springer-Verlag Netherlands; 1988:53–74
  20. 20.↵
    1. Ingrisch M,
    2. Schneider MJ,
    3. Norenberg D, et al
    . Radiomic analysis reveals prognostic information in T1-weighted baseline magnetic resonance imaging in patients with glioblastoma. Invest Radiol 2017;52:360–66 doi:10.1097/RLI.0000000000000349 pmid:28079702
    CrossRefPubMed
  21. 21.↵
    1. Han K,
    2. Song K,
    3. Choi BW
    . How to develop, validate, and compare clinical prediction models involving radiological parameters: study design and statistical methods. Korean J Radiol 2016;17:339–50 doi:10.3348/kjr.2016.17.3.339 pmid:27134523
    CrossRefPubMed
  22. 22.↵
    1. D'agostino RB,
    2. Nam BH
    . Evaluation of the performance of survival analysis models: discrimination and calibration measures. Advances in Survival Analysis 2003;23:1–25
  23. 23.↵
    1. Africa E,
    2. Pauciulo A,
    3. Vadalà R, et al
    . Primary central nervous system lymphoma: role of DWI in the differential diagnosis. Rays 2005;30:221–26 pmid:16512067
    PubMed
  24. 24.↵
    1. Guo AC,
    2. Cummings TJ,
    3. Dash RC, et al
    . Lymphomas and high-grade astrocytomas: comparison of water diffusibility and histologic characteristics. Radiology 2002;224:177–83 doi:10.1148/radiol.2241010637 pmid:12091680
    CrossRefPubMed
  25. 25.↵
    1. Herneth AM,
    2. Guccione S,
    3. Bednarski M
    . Apparent diffusion coefficient: a quantitative parameter for in vivo tumor characterization. Eur J Radiol 2003;45:208–13 doi:10.1016/s0720-048x(02)00310-8 pmid:12595105
    CrossRefPubMed
  26. 26.↵
    1. Khayal IS,
    2. Vandenberg SR,
    3. Smith KJ, et al
    . MRI apparent diffusion coefficient reflects histopathologic subtype, axonal disruption, and tumor fraction in diffuse-type grade II gliomas. Neuro Oncol 2011;13:1192–201 doi:10.1093/neuonc/nor122 pmid:21865401
    CrossRefPubMed
  27. 27.↵
    1. Wu X,
    2. Kellokumpu-Lehtinen PL,
    3. Pertovaara H, et al
    . Diffusion-weighted MRI in early chemotherapy response evaluation of patients with diffuse large B-cell lymphoma–a pilot study: comparison with 2-deoxy-2-fluoro- D-glucose-positron emission tomography/computed tomography. NMR Biomed 2011;24:1181–90 doi:10.1002/nbm.1689 pmid:21387451
    CrossRefPubMed
  28. 28.↵
    1. Yetkin FZ,
    2. Mendelsohn D
    . Hypoxia imaging in brain tumors. Neuroimaging Clin N Am 2002;12:537–52 doi:10.1016/S1052-5149(02)00029-1 pmid:12687910
    CrossRefPubMed
  29. 29.↵
    1. Le Duc G,
    2. Péoc’h M,
    3. Rémy C, et al
    . Use of T(2)-weighted susceptibility contrast MRI for mapping the blood volume in the glioma-bearing rat brain. Magn Reson Med 1999;42:754–61 doi:10.1002/(sici)1522-2594(199910)42:4<754::aid-mrm18>3.0.co;2-q pmid:10502765
    CrossRefPubMed
  30. 30.↵
    1. Knisely JP,
    2. Rockwell S
    . Importance of hypoxia in the biology and treatment of brain tumors. Neuroimaging Clin N Am 2002;12:525–36 doi:10.1016/S1052-5149(02)00032-1 pmid:12687909
    CrossRefPubMed
  31. 31.↵
    1. Valles FE,
    2. Perez-Valles CL,
    3. Regalado S, et al
    . Combined diffusion and perfusion MR imaging as biomarkers of prognosis in immunocompetent patients with primary central nervous system lymphoma. AJNR Am J Neuroradiol 2013;34:35–40 doi:10.3174/ajnr.A3165 pmid:22936096
    Abstract/FREE Full Text
  32. 32.↵
    1. Buttini M,
    2. Appel K,
    3. Sauter A, et al
    . Expression of tumor necrosis factor alpha after focal cerebral ischaemia in the rat. Neuroscience 1996;71:1–16 doi:10.1016/0306-4522(95)00414-9 pmid:8834388
    CrossRefPubMed
  33. 33.↵
    1. Zazulia AR,
    2. Diringer MN,
    3. Videen TO, et al
    . Hypoperfusion without ischemia surrounding acute intracerebral hemorrhage. J Cereb Blood Flow Metab 2001;21:804–10 doi:10.1097/00004647-200107000-00005 pmid:11435792
    CrossRefPubMed
  34. 34.↵
    1. Kim SK,
    2. Park JE,
    3. Kim KH
    , KSNO Guideline Working Group, et al. A national consensus survey for current practice in brain tumor management, III: brain metastasis and primary central nervous system lymphoma. Brain Tumor Res Treat 2020;8:20–28 doi:10.14791/btrt.2020.8.e7 pmid:32390350
    CrossRefPubMed
  • Received May 10, 2021.
  • Accepted after revision September 20, 2021.
  • © 2022 by American Journal of Neuroradiology
PreviousNext
Back to top

In this issue

American Journal of Neuroradiology: 43 (1)
American Journal of Neuroradiology
Vol. 43, Issue 1
1 Jan 2022
  • Table of Contents
  • Index by author
  • Complete Issue (PDF)
Advertisement
Print
Download PDF
Email Article

Thank you for your interest in spreading the word on American Journal of Neuroradiology.

NOTE: We only request your email address so that the person you are recommending the page to knows that you wanted them to see it, and that it is not junk mail. We do not capture any email address.

Enter multiple addresses on separate lines or separate them with commas.
Hypovascular Cellular Tumor in Primary Central Nervous System Lymphoma is Associated with Treatment Resistance: Tumor Habitat Analysis Using Physiologic MRI
(Your Name) has sent you a message from American Journal of Neuroradiology
(Your Name) thought you would like to see the American Journal of Neuroradiology web site.
CAPTCHA
This question is for testing whether or not you are a human visitor and to prevent automated spam submissions.
Citation Tools
Hypovascular Cellular Tumor in Primary Central Nervous System Lymphoma is Associated with Treatment Resistance: Tumor Habitat Analysis Using Physiologic MRI
S.Y. Jeong, J.E. Park, N. Kim, H.S. Kim
American Journal of Neuroradiology Jan 2022, 43 (1) 40-47; DOI: 10.3174/ajnr.A7351

Citation Manager Formats

  • BibTeX
  • Bookends
  • EasyBib
  • EndNote (tagged)
  • EndNote 8 (xml)
  • Medlars
  • Mendeley
  • Papers
  • RefWorks Tagged
  • Ref Manager
  • RIS
  • Zotero
Share
Hypovascular Cellular Tumor in Primary Central Nervous System Lymphoma is Associated with Treatment Resistance: Tumor Habitat Analysis Using Physiologic MRI
S.Y. Jeong, J.E. Park, N. Kim, H.S. Kim
American Journal of Neuroradiology Jan 2022, 43 (1) 40-47; DOI: 10.3174/ajnr.A7351
del.icio.us logo Digg logo Reddit logo Twitter logo CiteULike logo Facebook logo Google logo Mendeley logo
  • Tweet Widget
  • Facebook Like
  • Google Plus One
Purchase

Jump to section

  • Article
    • Abstract
    • ABBREVIATIONS:
    • MATERIALS AND METHODS
    • RESULTS
    • DISCUSSION
    • CONCLUSIONS
    • ACKNOWLEDGMENT
    • Footnotes
    • References
  • Figures & Data
  • Supplemental
  • Info & Metrics
  • References
  • PDF

Related Articles

  • PubMed
  • Google Scholar

Cited By...

  • No citing articles found.
  • Crossref
  • Google Scholar

This article has not yet been cited by articles in journals that are participating in Crossref Cited-by Linking.

More in this TOC Section

Adult Brain

  • Incidental Findings from 16,400 Brain MRI Examinations of Research Volunteers
  • Leptomeningeal Interfoliar Enhancement on Vessel Wall MR Imaging as a Unique Radiologic Finding of Susac Syndrome
  • Edge-Enhancing Gradient-Echo MP2RAGE for Clinical Epilepsy Imaging at 7T
Show more Adult Brain

Functional

  • Phenotyping Superagers Using Resting-State fMRI
  • Validity Assessment of an Automated Brain Morphometry Tool for Patients with De Novo Memory Symptoms
  • Identification of the Language Network from Resting-State fMRI in Patients with Brain Tumors: How Accurate Are Experts?
Show more Functional

Similar Articles

Advertisement

News and Updates

  • Lucien Levy Best Research Article Award
  • Thanks to our 2022 Distinguished Reviewers
  • Press Releases

Resources

  • Evidence-Based Medicine Level Guide
  • How to Participate in a Tweet Chat
  • AJNR Podcast Archive
  • Ideas for Publicizing Your Research
  • Librarian Resources
  • Terms and Conditions

Opportunities

  • Share Your Art in Perspectives
  • Get Peer Review Credit from Publons
  • Moderate a Tweet Chat

American Society of Neuroradiology

  • Neurographics
  • ASNR Annual Meeting
  • Fellowship Portal
  • Position Statements

© 2023 by the American Society of Neuroradiology | Print ISSN: 0195-6108 Online ISSN: 1936-959X

Powered by HighWire