Cite

Introduction

High grade glioma (HGG) still represent a challenge for the medical community worldwide.1, 2, 3, 4 Currently, average survival rates reported for anaplastic astrocytoma (AA) and glioblastoma (GBM) are 3–4 years and 10–16 months, respectively.5 One of the greatest concerns remains the lack of knowledge about the etiopathogenesis and pathophysiology of gliomas. The new WHO classification tried to explain the biological behavior of these tumors, analyzing their metabolic and genetic pathways.5 However, such classification has determined further fragmentation of our knowledge, in a time when we need a minimum common denominator to solve this terrible puzzle.6,7 In recent years, research focused on glioma neuroradiological features trying to relate these with their biological behavior.8, 9, 10, 11, 12, 13, 14, 15 As a consequence, many preoperative information about low grade gliomas (LGG) came from analysis of velocity of diametric expansion (VDE). Pallud and Mandonnet16, 17, 18, 19 reported a strict correlation between VDE and glioma aggressiveness: growth values under 4 mm between two consecutive MRI studies suggest a low invasiveness rate and thus lower malignancy.20,21 However, the same evidences are not available for HGG because of the lack of several preoperative MRI and no growth patterns have been standardized for HGG. We report a series of four patients in whom HGG growth and evolution were radiologically analyzed thanks to a long preoperative follow-up since tumor diagnosis until treatment.

Patients and methods

We describe a series of 4 patients followed by the first author from 2014 to January 2019. Two patients are male, two are female, two had a pathological diagnosis of GBM, one of AA and one received a neuroradiological diagnosis of GBM. Every patient was radiologically followed through almost three pre-operative MRI scans, focusing on T1 with gadolinium and T2/fluid attenuated inversion recovery (FLAIR) sequences. Tumor volume was calculated by neurosurgeons with experience in neuro-oncological field, using Horos software for MacOs for manual segmentation of T1 and FLAIR 3D volumetric images in each MRI scan.22 For each patient’s MRI scans mean tumor diameter (MTD)20 was calculated through the following formula, starting from FLAIR-altered region’s and enhancing nodule’s (EN) volume (V), and expressed in millimeters (mm):

MTD=2V3$$ \begin{equation}M T D=\sqrt[3]{2 * V} \end{equation}$$

It represents the reverse formula of ellipsoid volume formula calculation: D1D2D32 $\frac{D 1 * D 2 * D 3}{2}$assuming the three diameters as equivalent.

Variation of MTD (ΔMTD) $(\Delta M T D)$between two consecutive MRI scans of each patient, was calculated as follows:

ΔMTD=MTD(n+1)MTD(n)n0$$ \begin{equation}\Delta M T D=M T D_{(n+1)}-M T D_{(n)\cup}\, n \neq 0 \end{equation}$$

Where n is the ordinal number of MRI scan for each patient.

Velocity of diametric expansion (VDE)18,20, or MTD difference over time (days) between two consecutive MRI scans, was calculated as:

VDE=ΔMTDΔT$$ \begin{equation}\mathrm{VDE}=\frac{\Delta M T D}{\Delta T} \end{equation}$$

where ΔT is the number of days between two consecutive MRI scans. VDE is reported as millimeters on days (mm/dd). Assuming a linear VDE growth, VDE for each interval of time was then calculated respect to 1-solar-year, as follow:

 Annual VDE=VDE×365$$ \begin{equation}\text { Annual } \mathrm{VDE}=\mathrm{VDE} \times 365 \end{equation}$$

For each patient, mean acceleration (Acc) or VDE variation over time, was calculated as follows:

ACC=VDE(n+1)VDE(n)ΔT$$ \begin{equation}\mathrm{ACC}=\frac{V D E(n+1)-V D E(n)}{\Delta T} \end{equation}$$

Where n is the ordinal number of MRI scan for each patient; ΔT is the number of days between two consecutive MRI scans. VDE is expressed in mm/ dd, thus Acc results as mm/dd2.

In order to evaluate variation of VDE and Acc along observation time, further analysis were conducted for any patient. In particular, we observed peak of VDE and Acc values in all patients.

About VDE, we quantified these peaks highest value of VDE /mean VDE value ratio (listed in Table 3) in every patient along observation period, expressed as percentage [eq. 6]:

 VDE peak =vmax×nΣinvi×100$$ \begin{equation}\text { VDE peak }=\left(\frac{v_{\max } \times n}{\Sigma_{i}^{n} v_{i}}\right) \times 100 \end{equation}$$

Where v is the VDE value; vmax is the highest value of VDE recorded for any patient along all observation period; n total number of MRI scans for each patient; i the ordinal number of MRI scan.

In Acc peak analysis, we evaluated ratio between highest Acc value and mean Acc value for any patient (listed in table 4), expressed as percentage, as follows:

 Acc peak =αmax×ninαi×100$$ \begin{equation}\text { Acc peak }=\left(\frac{\alpha_{\max } \times n}{\sum_{i}^{n} \alpha_{i}}\right) \times 100 \end{equation}$$

Where α is the Acc value; αmax is the highest value of Acc recorded for any patient along all observation period; n total number of MRI scans for each patient, i the ordinal number of MRI scan.

Further analysis evaluated EN/FLAIR volumes ratio (EFVR) (eq. 8), EN/FLAIR MTD ratio (EFMR) and their evolution over observation time (eq. 10 and 11), as follows:

EFVR=ENvolume(n)FLAIRvolume(n)$$ \begin{equation}\mathrm{EFVR}=\frac{E N \, { volume }_{(n)}}{ {FLAIR\, volume}_{(n)}} \end{equation}$$EFMR=ENMTD(n)FLAIRMTD(n)$$ \begin{equation}\mathrm{EFMR}=\frac{E N M T D_{(n)}}{F L A I R M T D_{(n)}} \end{equation}$$ EFVR on single interval of time between consequent MRI scan =EFVR/T(n+1)T(n)$$ \begin{equation}\begin{aligned} &\text { EFVR on single interval of time between consequent MRI scan }=E F V R /\left[T_{(n+1)}-T_{(n)}\right] \end{aligned} \end{equation}$$EFVR on total interval of time between first and last MRI scan of each patient=EFVR/TnfTn0$$ \text {EFVR on total interval of time between first and last MRI scan of each patient} =E F V R /\left[T_{\left(n_{f}\right)}-T_{\left(n_{0}\right)}\right]$$ EFMR on single interval of time between consequent MRI scan =EFMR/T(n+1)T(n)$$ \text { EFMR on single interval of time between consequent MRI scan }=E F M R /\left[T_{(n+1)}-T_{(n)}\right]$$ EFMR on total interval of time between first and last MRI scan of each patient =EFMR/TnfTn0$$ \begin{equation}\begin{aligned} &\text { EFMR on total interval of time between first and last MRI scan of each patient }=E F M R /\left[T_{\left(n_{f}\right)}-T_{\left(n_{0}\right)}\right] \end{aligned} \end{equation}$$

Where n is the ordinal number of MRI scan for each patient; n0 and nf are the first and last MRI scan for each patient, respectively.

Every analysis was conducted for whole-tumor FLAIR volume and for single EN. The analysis and graphic representation were carried out with Microsoft Office Excel® software.

In addition, a literature review was performed using PubMed MEDLINE database and searching for “radiological glioma growth”, “glioma radiologic follow-up”, “glioma volume growth”.

Ethical standards

All procedures performed in studies involving human participants were in accordance with the ethical standards of the institutional and/or national research committee and with the 1964 Helsinki declaration and its later amendments or comparable ethical standards. Informed consent was obtained from all individual participants included in the study.

Results
Case 1

A 39 y.o. male patient was under radiological follow-up since 5 years after meningioma resection. MRI scan was obtained once a year. After 5 years of negative post-operative MRI, the exam showed hyperintensity in the right hippocampal and para-hippocampal areas on T2/FLAIR sequences, without enhancement after gadolinium infusion. MTD was 35.5 mm. For the sake of a closer follow-up, a new MRI scan with spectroscopy was scheduled 3 months later. In this second MRI scan, the FLAIR hyperintense lesion MTD was 39.7 mm, with no enhancement after gadolinium injection. Spectroscopy showed a peak in N-acetyl-aspartate (NAA), choline (Cho) and Creatinine (Cr) values. Both data suggested a likely diagnosis of ex novo HGG. Despite strong recommendation to undergo surgery, the patient delayed surgical treatment by 3 months. The new MRI scan, obtained the day before surgery for surgical planning, showed a MTD of 45 mm and contrast enhancement on T1-weighted sequences. The post-operative MRI showed a complete resection of the enhancing tumor. After resection, histopathologic diagnosis revealed AA, wild-type (WT) for IDH1. All values of ΔMTD, VDE and Acc calculated according to equations [2], [4,] [5], [6], [7] are listed in Table 1. After tumor gross total resection, the patient underwent the Stupp regimen for adjuvant treatment and died 16 months later.

Case 2

A 37 y.o. male patient, after an episode of dysarthria, underwent a neuroradiological workout. MRI images showed a left fronto-parieto-insular enhancing lesion suspicious for ex novo GBM. FLAIR and EN MTD were respectively 56.59 mm and 42.02 mm. After a week, he repeated the MRI with spectroscopy showing elevation of Cho and reduction of NAA values, elevation of Cho/NAA ratio, mild elevation of lactate signal. At the same time, central necrosis was visible. A diagnosis of GBM was made. FLAIR and EN MTD were 59.83 cc and 48 cc respectively. VDE was 169.23 mm/year for FLAIR tumor volume, whilst for EN was 311.81 mm/year. Due to poor neurological condition and tumor location, surgery was avoided. Following further neurological deterioration, he underwent a new MRI scan 2 weeks later. FLAIR tumor MTD was 72.60 mm, EN MTD was 52.16 mm. At this point, chemotherapy with temozolmide (TMZ) started. About two months later, a repeat MRI study showed good response. MTD was 67.37 mm for FLAIR tumor volume and 43.80 mm for EN, in which TMZ effect was higher (Figure 1). In fact, ΔMTD was - 5.22 mm and - 9.36 on FLAIR and EN nodule volume, respectively. Values of ΔMTD, VDE and Acc calculated according to equations [2], [4,] [5], [6], [7] for every MRI scan are listed in Table 1, divided for FLAIR and EN volumes. 6 months after initial diagnosis the patient remains alive.

Figure 1

MRI scans in Patient 2 at 0 (A), 18 (B) and 54 (C) days. (T1 with gadolinium sequences).

Data about patient #1, #2, #3, #4

Acc = acceleration; MRI = magnetic resonance imaging; MTD = mean tumor diameter; ΔMTD = variation of MTD respect to previous value; VDE = velocity of diameter expansion; V EN = volume of enhancing nodule post gadolinium; V FLAIR = volume of fluid attenuated inversion recovery altered regions

Case 3

A 72 y.o. female patient, after a seizure, underwent MRI workout without pathological findings. She repeated a MRI scan after 15 months and a FLAIR alteration in the left motor areas was seen. MTD was 21.25 mm. She refused any treatment. MRI scan was repeated further 9 months later showing a bilateral alteration of FLAIR signal in the frontal lobes. MTD was 26.30 mm. However, the patient chose to undergo neuroradiological follow-up rather than surgery. After two months, a new MRI scan showed an EN in the right frontal lobe. MTD of EN was 8.43 mm, MTD of FLAIR altered regions was 32.01 mm. Neurosurgeons suggested excision of the EN, but the patient refused again. Further MRI scans at 2, 6 and 10 weeks were obtained. MTD of EN were 11.19, 17.54 and 21.40 mm, respectively. MTD of FLAIR-altered regions were 35.73, 41.68 and 41.68 mm, respectively (Figure 2). At the end complete resection of EN was obtained; histopathological diagnosis was compatible with GBM, IDH1 WT. She underwent adjuvant therapy according to the Stupp regimen and remains alive 31 months after the index surgery. All values of ΔMTD, VDE and Acc calculated according to equations [2], [4,] [5], [6], [7] for every MRI scan are listed in Table 1, divided for FLAIR and EN volumes.

Data on fluid attenuated inversion recovery (FLAIR) tumor and enhancing nodule (EN) volume velocity of diameter expansion (VDE). Reported as mm/year

Figure 2

MRI scans in Patient 3 at 0 (A), 465 (B), 742 (C), 790 (D), 811 (E), 825 (F) and 854 (G) days (fluid attenuated inversion recovery [FLAIR] and T1 with gadolinium sequences).

Case 4

A 56 y.o. female patient underwent MRI for persisting headache. An incidental slight alteration of FLAIR signal in the left gyrus rectus was discovered. MTD was 7.37 mm. She repeated an MRI 5 months later to evaluate tumor evolution and an EN appeared in the previously FLAIR-altered area. MTD of the EN was 7.37 mm, whereas MTD of FLAIR-altered area was 12.6 mm and there were no further changes in the next MRIs. An MRI with spectroscopy was obtained one month later and this suggested a GBM. The new MTD of the EN was 8.43 mm. Surgical treatment was then proposed. She underwent preoperative MRI for neuronavigation one month later. MTD of the EN was 8.43 mm (Figure 3). Complete resection of EN was carried out and histopathological diagnosis revealed GBM IDH1 WT. She underwent adjuvant therapy with Stupp regimen and is still alive 12 months after surgery. All values of ΔMTD, VDE and Acc calculated according to equations [2], [4,] [5], [6], [7] for every MRI scan are listed in Table 1, divided for FLAIR and EN volumes.

Figure 3

MRI scans in Patient 4 at 0 (A), 160 (B), 196 (C) and 223 (D) days (fluid attenuated inversion recovery [FLAIR] and T1 with gadolinium sequences).

General evaluations

Each patient underwent sequential MRI studies, with an average number of 4 scans along a mean period of 413 days. ENs were evident only in patient 2 on the first MRI study, whereas appeared at fourth, fourth and second MRI scan in patients n. 1, 3 and 4, respectively. No common time correlation between evidence of first FLAIR alteration and EN appearance was found. Median appearance time after first MRI scan was 432 days.

Final mean tumor volume, evaluated on each patient’s last T2/FLAIR MRI sequences, was 59 cc, with a median volume variation along sequential MRIs of 6.4 cc.

Mean MTD at last MRI was 41,66 mm and 24.54 mm for FLAIR and EN volume, respectively, whilst median ΔMTD along the observation period was 5.4 mm and 2.11 mm, for FLAIR and EN volume, respectively (Table 1).

Mean VDE was 39.9 mm/year and 45.2 mm/year for FLAIR volume and EN, respectively (Table 1, 2).

VDE analysis revealed peak in variation along time. Mean percentage maximal VDE/mean VDE ratio, calculated according to eq. 6, was 221.49% and 259.47% for FLAIR and EN volumes, respectively (Table 3).

Raw data about velocity of diameter expansion (VDE) of fluid attenuated inversion recovery (FLAIR) and enhancing nodule (EN) volume. Mean percentage maximal VDE/mean VDE ratio could help in early high grade gliomas (HGG) diagnosis

Interval of consecutive MRIVDE Pt #1 FLAIRVDE Pt #2 FLAIRVDE Pt #2 ENVDE Pt #3 FLAIRVDE Pt #3 ENVDE Pt #4 FLAIRVDE Pt #4 ENMean Values
135.50169.23311.8116.680.0011.9316.81
217.03258.82104.546.660.000.0010.81
332.24-52.97-94.9143.3964.140.000.00
464.5947.84
5155.18165.74
60.0048.53
Mean Values28.26125.03107.1547.7554.373.989.21
FLAIR VDE Peak/ VDE mean value ratio (percentage)125.63135.36324.99300.00221.49
EN VDE Peak/ VDE mean value ratio (percentage)291.02304.81182.57259.47
Mean Values240.48

Acceleration analysis according to [eq. 5] revealed a mean value of 0.0077 mm/dd2 and 0.0119 mm/days2 for FLAIR and EN volume, respectively. Mean percentage ratio between peak values and mean value of acceleration was 282.7% for FLAIR volume and 257.52% for EN, respectively (table 4). In patient n. 1 late occurrence of the EN didn’t allow complete evaluation of MTD variation, VDE and Acc (Figure 4).

Figure 4

Velocity of diameter expansion (VDE) variation of mean tumor diameter (MTD) for fluid attenuated inversion recovery [FLAIR] and enhancing nodule (EN) volumes for each patient. Plus the case reported by Mandonnet et al. as specified in Discussion paragraph. Graphic shows variation of acceleration (Acc) based on VDE (mm/days2) for FLAIR tumor volume and EN volume along time. Except obvious similar variation for two types of volume considered, peaks are easily notable. VDE of patient #1 is covered by patient #4.

Raw data of mean acceleration (Acc) based on velocity of diameter expansion (VDE) of fluid attenuated inversion recovery (FLAIR) and enhancing nodule (EN) volume. Mean percentage maximal Acc/mean Acc ratio of each patient could help in early high grade gliomas (HGG) diagnosis

Interval of consecutive MRIAcc Pt #1 FLAIRAcc Pt #2 FLAIRAcc Pt #2 ENAcc Pt #3 FLAIRAcc Pt #3 ENAcc Pt #4 FLAIRAcc Pt #4 ENMean Values
10.00030.06620.12200.0001000.00020.0003
20.00050.03940.01590.0000700.00000.0008
30.0015-0.0040-0.00720.002480.00370.00000.0000
40.008430.0062
50.030370.0324
60.000000.0046
Mean Values0.00080.03390.04360.00690.00780.00010.0004
FLAIR Acc Peak/ Acc mean value ratio (percentage)195.6695195.5808439.7477300282.7495
EN Acc Peak/ Acc mean value ratio (percentage)280.0617414.751677.74935257.5209
Mean Values270.1352

Analysis of EFVR and EFMR, according to eq. 813 did not show significant results.

Discussion

Ex novo HGG diagnosis usually becomes evident when the tumor mass is wide and resection must be carried out as soon as possible. Thus, for each patient no many preoperative MRI scans are available for mathematical analysis of tumor growth. In the past decades, imaging studies on LGG evolution23,24 let neurosurgeons to change their mind on treatment: from a “wait-and-see” approach to an early resection one, due to the understanding of malignant transformation risk.17,18 Indeed, Mandonnet et al. published in 200819 a pioneering paper on the role of VDE in LGG and Pallud, few years later, reported that a VDE over 8 mm/year is an independent negative prognostic factor for survival in LGG patients.20 In the past, few studies reported mathematical models of growth based on simple exponential growth of Volume Doubling Time (VDT)25,26 or linear growth of tumor radius.27 In 2015 Stensjoen et al. analyzed tumor’s growth curve in 106 patients, considering two consecutive MRI scans, and reported the Gompertzian growth to be the most reliable mathematical model to predict tumor mass evolution; and this had already been reported by previous studies.9,28,29 As seen in LGG, valuable evaluations of malignancy and evolutions could come from MTD analysis and its variation among sequential MRI studies. Mandonnet et al. reported an incidental insular tumor, followed for 9 months: VDE was 12 and 14 mm/year at 5 and 9 months, respectively. Histopathologic analysis then revealed an AA.21 Cochereau et al. reported a secondary GBM, in which VDE after evolution was 13 mm/year.30 In such cases, authors suggest that an annual VDE value over 8 mm/year should represent an alarm for HGG suspect. Our experience with ex novo HGG, shows even higher mean values for annual VDE: 39.91 and 45.02 mm/year for FLAIR and EN volume, respectively, enforcing previous findings and suggesting a more aggressive behaviour of an ex novo HGG respect the secondary ones. Wide differences in VDE values between our cases and cases reported by other authors, might come from different pathogenisis (primary versus secondary HGG), different observation times and biological features of the brain tumors. However, we want to underline how the VDE variate along the observation time. FLAIR VDE shows peaks, particularly in patients n.3 and n.4, with a ratio between peak and mean VDE value of 324.99% (n = 3.25 mm/dd) and 300% (n = 3 mm/dd), respectively. As shown in Figure 1, VDE variations follow a similar pattern in FLAIR and EN volume: progressive growth reach the highest value and then a rapid decrease is registered. Even if such peaks happen in different moments for each patients, such difference may be caused by different starting observation time. Speculating the possibility to follow growth curves of different HGG, all with the same onset time, peaks and following decrease could all happen in the same interval of time. Similar pattern, as in our cohort, could be found even in Mandonnet et al. case report.21 Moreover the reported mean VDE of the EN is higher than FLAIR volume VDE suggesting the higher biological aggressiveness of the enhancing regions. In addition, not yet reported to our knowledge, mean Acc values show different peaks along tumor life. In our experience mean percentage ratio between peak and mean Acc values was 270.13%. Last data could improve early diagnosis of suspected primary or secondary HGG, even if a consensus about MRI timing doesn’t exist. Follow-up MRI scans could be proposed every 4 months31 to patients with early-detected LGG. Pallud et al. proposed a timing between 6 weeks and 3 months according to Chang score.16 However, the wait and see policy in suspected HGG can let tumors grow and induce deficits. Thus, 3 months interval between two consecutive MRIs could be proposed when tumors could be misdiagnosed. Conversely, if strong suspect for HGG exists, observation could be carried out monthly. Indeed, Acc peak over 270% or VDE peak over 240% in comparison to mean values, for both FLAIR and EN volume, could help in diagnosis. Moreover, no long interval time is needed between MRIs: in our experience Acc variation could be seen even after 15–20 days. Obviously, at least three consecutive MRI scans are needed for VDE and Acc analysis. In summary, even if large series should confirm these data, we strongly propose surgical treatment in case of Acc peak and rapid VDE evolution, due to the high probability that MRI alterations refer to HGG. Such analysis could be helpful as prognostic factor and evolution parameter in case of surgical untreatable patients, i.e. due to poor general conditions or patient refusing surgery. However, sudden changes in VDE and Acc values do not have known causes. A clue could come from evolution of secondary gliomas, in which a previous LGG has been reported. By now, accumulation of mutation in LGG genes, particularly on chromatin modelling system, let a more aggressive phenotype acquisition. Maybe, as no long observation time has been obtained for HGG, our data could report a crucial moment of HGG evolution. Such peaks could refer to gain of mutation in high percentage of proliferating cells. Thus, MTD could increase in EN and FLAIR due to peripheral disposition of viable glioma cells. Such an expansion would not be compensated by core necrosis. However, too many different factors, i.e. mechanical, growth-factors-related, neo-angiogenesis, supervene in HGG infiltration and evolution. This could also explain our different results in time of peak presentation since first diagnosis, associated with the unknown time of tumor onset. It would be interesting to analyze such variation in animal models, even if inoculation of glioma cells is an important bias in infiltration pattern study. A sudden increase in MTD as reported by VDE and Acc analysis, a biopsy and genetic analysis could be investigated and related if possible.

Conclusions

To our knowledge, this is the first study describing parameters of acceleration and VDE peak in a retrospective series of ex novo HGG. Through sequential MRIs and evaluation of these parameters, an early diagnosis of HGG could be done. In addition, time between consecutive MRI scans could be reduced to 15–20 days. This work represent a pilot study on new findings of HGG growth curve, never investigated, due to short observation time in case of HGG suspect. Despite limitations in number of patients and observation time, if our data are confirmed in larger series, such values could be helpful as prognostic factors in non-surgical patient. In an attempt to obtain further confirmations, a multicenter register could be create, where all data about long-preoperative-follow-up of HGG could be reported. Moreover, studies on in-vivo and invitro models could associate genome mutations with such phenotypical variation, looking for new biological tumor pathways of interest.

eISSN:
1581-3207
Language:
English
Publication timeframe:
4 times per year
Journal Subjects:
Medicine, Clinical Medicine, Internal Medicine, Haematology, Oncology, Radiology