Memory recovery is related to default mode network impairment and neurite density during brain tumours treatment
================================================================================================================

* R Romero-Garcia
* J Suckling
* M Owen
* M Assem
* R Sinha
* P Coelho
* E Woodberry
* SJ Price
* GAA Burke
* T Santarius
* Y Erez
* MG Hart

## ABSTRACT

Surgical resection is an effective treatment to delay brain tumour progression and improve survival. Nevertheless, a large proportion of patients have treatment-induced memory deficits that dramatically reduce their life quality. Our aim is to test brain tumour interactions with brain networks thereby identifying protective features and risk factors for memory recovery. Seventeen patients with diffuse glioma (aged 22-56 years) were MRI-scanned and cognitively assessed before and after surgery, and during a 12-months recovery period (47 MRI in total after exclusion). Using structural MRI and Neurite Orientation Dispersion and Density Imaging (NODDI) derived from diffusion-weighted images, we respectively estimated lesion overlap and Neurite Density with brain networks derived from normative data in healthy participants. Non-local effects of tumours were evident on neurite density, which was reduced not only within the tumour, but also beyond the tumour boundary. In contrast, high preoperative neurite density outside the tumour, but within the DMN, was associated with better memory recovery. Furthermore, follow-up neurite density within the DMN and Frontoparietal network were also associated with memory recovery. Preoperative tumour, and post-operative lesion, overlap with the DMN showed a significant negative association with memory recovery. Our findings suggest that brain tumours and their corresponding treatment affecting brain networks that are fundamental for memory functioning such as the DMN can have a major impact on patient’s memory recovery.

KEYWORDS
*   Brain Tumours
*   memory
*   Default Mode Network
*   diffusion MRI
*   Neurosurgery
*   Brain Networks

## INTRODUCTION

Every year more than 300,000 people worldwide face the diagnosis of a brain tumour. Patients who undergo surgical resection rather than biopsy have a better overall survival, and extending the resection beyond the abnormality seen on MRI may further improve prognosis. However, the extent of resection is only a worthwhile prognostic factor in the management of the tumour if subsequent cognitive functioning can be maintained. Cognitive function, despite being recognised as a fundamental outcome measure of treatment success, remains one of the most unpredictable aspects of patients’ prognosis. As a result of their tumours, a significant proportion of patients develop cognitive deficits ranging from 29% in patients with non-irradiated low-grade glioma, to about 90% in patients following chemo-radiotherapy for metastases[Meyers et al., 2004]. Among the wide range of acquired cognitive deficits, memory disabilities represent a major challenge for preserving patient’s quality of life, including job preservation[Rusbridge et al., 2013]. Working[Teixidor et al., 2007], verbal[De Groot et al., 2013] and non-verbal[Owen et al., 2020] memory deficits have been extensively reported in brain tumour patients. Pharmacological interventions on these patients improve memory recovery[Brown et al., 2013] which has been validated in animal models[Schültke et al., 2008]. Unfortunately, the impact of surgery on cognitive outcome and memory has traditionally been underappreciated, potentially due to the complex distributed networks involved in normal memory function, and the limited understanding on how focal lesions interact with these networks[Sagberg et al., 2017].

Recently developed MRI sequences that estimate tissue microstructure integrity have revolutionized *in-vivo* investigation of the human brain, but their applications to neuro-oncology have been limited. The capability of modern MRI protocols has been demonstrated for estimates of tumour infiltration, grading, heterogeneity, progression and patient survival[Villanueva-Meyer et al., 2017]. The incorporation of recent MRI protocols into the clinical routine may be useful in reducing treatment-induced neurocognitive dysfunction in brain tumours. Notwithstanding, most previous brain tumour studies have used tensor models to explore microstructural changes under the assumption of Gaussian diffusion processes. In contrast, the Neurite Orientation Dispersion and Density Imaging (NODDI) technique uses varying gradients strengths to provide more specific measurements of tissue microstructure than traditional Diffusion Tensor Imaging (DTI)[Zhang et al., 2012]. Intra-cellular diffusion estimation derived from NODDI as a marker of Neurite Density has been validated using histochemical analysis in mouse models[Wang et al., 2019] and has shown sensitivity to neurotypical ageing[Nazeri et al., 2015] and for detecting axonal degeneration in neurological (such as Alzheimer’s[Colgan et al., 2016], multiple sclerosis[Grussu et al., 2017]) and psychiatric conditions[Rae et al., 2017]. Caverzasi et al. (2016) showed that even for tumour lesions that appear to be homogeneous on corresponding fluid-attenuated inversion-recovery (FLAIR) images, NODDI has the potential to differentiate infiltrative tumour components[Caverzasi et al., 2016]. Additionally, NODDI has demonstrated to be an accurate predictor of glioma grade[Zhao et al., 2018a]. Despite these promising findings, no study has yet used NODDI to evaluate the impact of tumours and their treatment on memory.

A major challenge when predicting the consequences of surgery and chemo-radiotherapy is that higher-order cognitive functions such as memory are sub-served by communication across spatially extended neural circuits that will require whole-brain analysis approaches. Large-scale brain networks offer a holistic framework to analysing the brain as a circuit of interacting components that are critically related to cognition[Bressler and Menon, 2010]. In neuro-oncology, markers derived from brain network approaches have been demonstrated to be sensitive to the presence of low-grade gliomas, plasticity differences between low- and high-grade gliomas, and surgically induced alterations[Huang et al., 2014]. Consequently, the use of brain network data to explore the potential for memory disruption induced by brain tumours and their treatment could be of major clinical relevance.

In this longitudinal study, we prospectively investigated a cohort of patients with diffuse glioma and whether MRI could be used to predict cognitive recovery postoperatively and during follow-up. Specifically, we wished to test the effect of tumours both locally and globally on the brain, in terms of their structural and functional network effects. We hypothesized that treatment-induced memory deficits are caused by disruption of brain networks that have been previously identified as fundamental for memory, and that this is mediated by global effects on normal brain microstructure.

## METHODS

### Sample

This study is a single centre prospective cohort design approved by the Cambridge Central Research Ethics Committee (protocol number 16/EE/0151). Patients deemed to have typical appearances of a diffuse glioma were identified at adult neuro-oncology multidisciplinary team (MDT) meetings at Addenbrooke’s Hospital (Cambridge, UK), and a consultant neurosurgeon directly involved in the study identified potential patients based on the outcome of the MDT discussion. All patients gave written informed consent. Inclusion criteria included: (*i*) participant is willing and able to give informed consent for participation in the study, (*ii*) imaging is evaluated by the MDT and judged to have typical appearances of a diffuse glioma, (*iii*) Stealth MRI is obtained (routine neuronavigation MRI scan performed prior to surgery), (*iv*) World Health Organisation (WHO) performance status 0 or 1, (*v)* age between 18 to 80 years, (*vi*) tumour located in or near eloquent areas of the brain thought to be important for speech and executive functions, and (*vii*) patient undergoing awake surgical resection of a diffuse glioma. This last inclusion criterion was adopted to collect additional intraoperative electrocorticography data that will be reported separately. Participants were excluded if any of the following applied: *i*) concomitant anti-cancer therapy (except for dexamethasone treatment), (*ii*) history of previous malignancy (except for adequately treated basal and squamous cell carcinoma or carcinoma in-situ of the skin) within 5 years, and (*iii*) previous severe head injury.

We recruited 17 patients aged 22-56 years (8 females). Final histological diagnoses revealed different grades of glioma: WHO-I n=2, WHO-II n=7, WHO-III n=5, WHO-IV n=3. Resection was complete (no residual FLAIR signal) in 9 patients, whereas 8 patients had partial resection. Adjuvant chemoradiotherapy was performed in 12 patients. Each patient was scanned up to four times: before surgery (preop), within 72 hours after surgery (postop), and at 3 and 12 months after surgery (month-3 and month-12). See **Table 1** for demographic details of participants.

View this table:
[Table 1.](http://medrxiv.org/content/early/2020/10/22/19008581/T1)

Table 1. Demographic and pathological information.
# Number of patient, SFG, Superior Frontal Gyrus; MFG, Middle Frontal Gyrus; IFG, Inferior Frontal Gyrus; ITG, Inferior Temporal Gyrus; MTG, Middle Temporal Gyrus; SMA, Supplementary Motor Area; RT, radiotherapy; Astro – astrocytoma; GBM – glioblastoma multiforme; Oligo – oligodendroglioma; Lesion Vol Pre/Post – Total volume occupied by the tumour (Pre) and total amount of damaged tissue (Post) according the lesion mask manually drawn on the MPRAGE image and refined with Unified Segmentation Lesion toolbox.

### MRI and NODDI data acquisition and pre-processing

MRI data were acquired using a Siemens Magnetom Prisma-fit 3 Tesla MRI scanner and 16-channel receive-only head coil (Siemens AG, Erlangen, Germany). A T1-weighted MRI MPRAGE sequence was acquired using the following parameters: repetition time (TR) = 2300 ms, echo time (TE) = 2.98 ms, flip angle (FA) = 9 deg, 1 mm3 resolution, Field of View (FOV) = 256×240 mm2, 192 contiguous slices and acquisition time of 9 minutes and 14 seconds. During the same scanning session, we used a recently developed MRI multi-shell diffusion technique, NODDI with 30 directions (b-value=800 mm/s), 60 directions (b-value=2000 mm/s) and 10 unweighted B0 images. Other acquisition parameters were: TR = 8200 ms, TE = 95 ms, 2.5 mm3 voxel resolution, 60 slices, FOV = 240 mm and acquisition time of 15 minutes and 19 seconds. NODDI Matlab Toolbox ([http://mig.cs.ucl.ac.uk/index.php?n=Tutorial.NODDImatlab](http://mig.cs.ucl.ac.uk/index.php?n=Tutorial.NODDImatlab)) was used to quantify the microstructural complexity of dendrites and axons *in vivo*[Zhang et al., 2012]. Compared with traditional DTI, this multi-compartment tissue model disentangles two key contributing factors of Fractional Anisotropy: the Gaussian contribution from water molecules located in the extracellular space (defined as the space around neurites), and the restricted non-Gaussian diffusion that takes place in the intra-cellular space that is bounded by axonal and dendritic membranes. The apparent intra-cellular volume fraction that represents the fraction of dendrites and axons was used here as a measurement of Neurite Density.

### Lesion masking and image transformation to standard space

Images were corrected for B0 field inhomogeneity, Gibbs artefacts and eddy-current distortions using MRtrix v3 ([https://www.mrtrix.org/](https://www.mrtrix.org/)) and FSL v5.0 ([http://fsl.fmrib.ox.ac.uk](http://fsl.fmrib.ox.ac.uk)). Scans were visually inspected by an experienced researcher (RRG). Masks of the pre-op tumour and follow-up lesion (reflecting for example resected tissue, residual tumour, post surgical oedema, or gliosis) masks were created using a semi-automated procedure. An experienced neurosurgeon (MH) initially did a manual delineation for each participant on the preoperative T1 image slices that included the tumour, and the resection site on the follow-up T1 images. Resulting masks were refined by the Unified Segmentation with Lesion toolbox ([https://github.com/CyclotronResearchCentre/USwithLesion](https://github.com/CyclotronResearchCentre/USwithLesion)) that uses tissue probability maps to create a posterior tumour/lesion probability map. Inter-regional distances to the tumour boundary as defined in the tumour mask were estimated as the geodesic distance of the shortest path constrained by the white matter.

For each scan, the first B0 image of the diffusion sensitive sequence was linearly coregistered to the T1 image using ANTs ([http://stnava.github.io/ANTs/](http://stnava.github.io/ANTs/)). The resulting inverse transformation was used to map the Neurite Density map into the T1 image space. Each T1 image was non-linearly coregistered to standard space using Symmetric Normalization (ANTs-SyN), but excluding the tumour/lesion mask from the non-linear step of the wrapping to avoid distorting the spatial distribution of the tumour/lesion. The resulting transformation was additionally used to map the lesion mask and the Neurite Density map from T1 space to standard template space. The ICBM 2009a symmetric brain, an unbiased non-linear average of the MNI152, was used here as a standard template for normalization of Neurite Density using the contralateral values hemisphere that contained the tumour.

### Networks atlas based on normative data

We additionally utilised the map of large-scale networks defined in Yeo et al. (2011). This atlas was created by clustering functionally coupled regions in 1000 young, healthy adults. Regions delimited on the 7-Network liberal version of the Yeo atlas were used as Regions-of-Interest (ROIs) for calculating tumour overlap and Neurite Density (**Figure 1**, top). The visual network was not included in the analyses because none of the patients presented a tumour in the occipital lobe. The limbic network was also excluded from further analyses due to its poor signal to noise ratio and lack of reproducibility[Gordon et al., 2016].

![Figure 1.](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2020/10/22/19008581/F1.medium.gif)

[Figure 1.](http://medrxiv.org/content/early/2020/10/22/19008581/F1)

Figure 1. Flowchart of pipeline analysis.
After transforming lesion masks to standard space, tumour/lesion spatial overlap and Neurite Density were calculated for each Yeo network (Yeo et al. 2011). For each assessment (preoperative, postoperative, 3 months and 12 months), lesion overlap and Neurite Density were compared with memory recovery. DMN – default mode network

### Tumour overlapping and Neurite Density estimation

As location of the tumour varied across patients, we calculated, for each network and participant, the “Tumour Overlap Index” as the amount of regional volume (in mm3) that spatially overlapped with the tumour (for preoperative assessments) or lesioned tissue (for postoperative and follow-up assessments) according to the tumour/lesion masks, after being transformed to standard space.

Median Neurite Densities in each Yeo network were calculated for each participant. Only voxels of the network that were not overlapping with the tumour/lesion masks were included in the analysis to reduce the impact of tumour volume on Neurite Density estimation. See **Figure 1** for an illustrative flowchart.

### Cognitive assessment

Immediately after each MRI scanning session, memory performance was evaluated using the Oxford Cognitive Screen (OCS)-bridge tablet-based a tool ([https://ocs-bridge.com/](https://ocs-bridge.com/)), which is specifically designed for clinical setting. OCS-bridge is an app-based extension of the Birmingham Cognitive Screen and Oxford Cognitive Screen, two paper and pencil assessment tests that are widely used in health services in UK and that have been translated into multiple languages[Kong et al., 2016].

We performed 6 memory screening tasks that included: free verbal memory, overall verbal memory, episodic memory, orientation, forward digit span and backwards digit span. For each individual task, a z-score value was computed by subtracting the mean and dividing by the standard deviation of task scores across participants. Tasks defined by an inverted scale where high values represent low memory performance were flipped by multiplying the z-score by −1. Total memory score was calculated as the average z-score of all tasks. By subtracting the preoperative z-score from the z-score of each subsequent task, we defined the longitudinal trajectory of each assessment (Δ) where negative scores represent worse performance than before surgery (*i*.*e*. memory deficit) and positive scores are associated with increased performance than before surgery (*i*.*e*. memory improvement). We have recently demonstrated that memory scores of these brain tumour patients are consistent between OCS-bridge and a traditional 2-3 hours neuropsychological interview[Owen et al., 2020].

### Statistical analysis

Associations between cognition and imaging markers were tested using Linear Mixed Effect Models (LMMs) to account for the dependency of repeated measurements: ![Formula][1]</img>  where Δ *Memory* represents the total memory score normalized to preoperative performance, *X* is the predictor variable based on imaging data (*i*.*e*. tumour overlap or Neurite Density for a given network) and (1 |*subject*) represents the random intercept of the model. The model included as covariates: age, gender, grade (low-vs. high-grade glioma), treatment (surgery alone vs. surgery plus chemo-radiotherapy), tumour volume and total Neurite Density (*ND*, only for predictions based on ND). All statistical tests were corrected for multiple comparisons using Benjamini–Hochberg False Discovery Rate (FDR < 0.05).

## RESULTS

### Memory recovery trajectories

Tumours were located on frontal (5 Left Hemisphere -LH- and 4 Right hemisphere -RH-), temporal (6 LH) and insular (2 RH) lobes (**Figure 2A, Table 1**). OCS-Bridge cognitive assessment was completed by 17 patients before surgery, 8 after surgery, 11 after 3 months and 11 after 12 months (47 in total). Seven patients could not tolerate post-operative screening and the remain missing assessments could not be collected due to logistic and technical reasons (14 out of 69). Surgical resection and treatment had an impact on memory in most participants. Memory recovery after surgery (ΔMemory) showed a variety of trajectories, including progressive impairment, impairment followed by recovery, no change and improvement after surgery and during recovery.

![Figure 2.](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2020/10/22/19008581/F2.medium.gif)

[Figure 2.](http://medrxiv.org/content/early/2020/10/22/19008581/F2)

Figure 2. Patients’ memory recovery trajectories as a function of tumour location and type of treatment.
**(a)** Tumour location in the cohort. Colour scale indicates the number of participants that had a tumour in a given cortical region. Cortical surface representations were plotted using BrainsForPublication v0.2.1 ([https://doi.org/10.5281/zenodo.1069156](https://doi.org/10.5281/zenodo.1069156)). **(b)** Memory recovery trajectories normalised to pre-operative values coloured according to tumour location (frontal, blue; temporal, red; insular, green). **(c)** Recovery trajectories of patients receiving no further treatment beyond surgery (solid lines) and patients that also underwent chemo-radiotherapy (dashed lines). **(d)** Association between total lesion volume and memory recovery (postoperative, 3 months and 12 months) across patients. Memory recovery was normalized to pre-surgical values. Each color represents an individual participant as presented in **Table 1**. Positive scores correspond with post-surgical memory recovery and negative scores to post-surgical deficits.

### No impact of tumour location or volume on memory recovery

Follow-up recovery trajectories were not associated with tumour location based on traditional lobe classification (LMM, *P*=0.11, **Figure 2B**) nor with the type of treatment received by patients (*i*.*e*. surgery alone or surgery with adjuvant chemoradiotherapy; LMM, P=0.98, **Figure 2C**). The influence of the tumour and its treatment on recovery was initially evaluated by comparing the total lesion volume with changes in memory. Memory recovery was not correlated with tumour volume before surgery (ρ=-0.27, *P*=0.40), nor with lesion volume during recovery (LMM, ρ=-0.51, *P*=0.097, **Figure 2D**).

### Peri-tumoural neurite density within brain networks is associated with memory recovery

The impact of tumours and their treatment on brain structure were explored using a Neurite Density marker derived from diffusion imaging based on NODDI. We found that average whole-brain Neurite Density outside the tumour was negatively associated with tumour volume (*R**2*= 0.28, *P*=0.029, **Figure 3A**), suggesting that tumours may mediate global effects on brain structural integrity. In support of this hypothesis, we found a distance-effect on Neurite Density as a function of distance to the tumour boundary. Peri-tumoural regions located between 0 and 20 mm outward from the tumour boundary had up to half of the Neurite Density compared to contralateral regions (**Figure 3B**), suggesting a extended tumour and brain structural integrity interactions. Preoperative Neurite Density within the DMN in both hemispheres (*i*.*e*. not normalized) was positively correlated with memory deficits during recovery, suggesting that Neurite Density has an initial protective effect on outcome (ρ=0.64, *P**fdr*<10−4, **Figure 3C**). Moreover, Neurite Density during recovery was also associated with memory performance. When postoperative and follow-up Neurite Density values were considered, Memory scores were correlated with Neurite Density within Frontoparietal (ρ=0.60, *P**fdr*=0.0015) and DMN (ρ=0.71, *P**fdr*=0.002) (**Figure 3D**). Neither preoperative nor follow-up Neurite Density with any of the remain networks were associated with memory recovery (all *P**fdr*>0.05). All associations were tested using LMM after regressing out age, gender, tumour grade, total lesion volume, type of treatment and total Neurite Density effects.

![Figure 3.](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2020/10/22/19008581/F3.medium.gif)

[Figure 3.](http://medrxiv.org/content/early/2020/10/22/19008581/F3)

Figure 3. Impact of tumour and tumour treatment in Neurite Density and memory recovery.
**(a)** Association between total tumour volume and Neurite Density outside of the tumour. **(b)** Neurite Density as a function of distance to the tumour boundary. Values are normalized to the contralateral hemisphere (*i*.*e* values less than one represent regions with reduced neurite density). Zero distance (x=0 mm) corresponds to average Neurite Density values within each tumour. Each color represents an individual participant. Due to high-dimensional data, error bars are too small to be displayed. **(c)** Memory recovery (postoperative, 3 months and 12 months) as a function of *preoperative* Neurite Density within the DMN. For a given participant, the same Neurite Density (preoperative) value was used in all assessments, resulting in vertical lines between them. **(d)** Memory recovery as a function of *postoperartive and follow-up* Neurite Density within the DMN Fronto-Parietal Network. Neurite Density corresponds with the median density within the White Matter of each network, excluding lesioned regions. None of the other networks showed any significant association that survived correction for multiple comparisons. Colors illustrate each individual participant with lines connecting them.

### Tumour overlap with the DMN is associated with memory recovery

Preoperative tumour overlap with the DMN was negatively correlated with memory recovery (ρ=-0.55, *P**fdr*=0.005; **Figure 4A**). When follow-up imaging data were included in the analyses (post-operative, 3 months and 12 months), we also found a significant association between lesion overlapping with the Default Mode Network (DMN) and memory performance (ρ=-0.56, *P**fdr*<10−4; **Figure 4B**). Neither (preoperative) tumour nor (follow-up) lesion overlapping with any of the remain networks were associated with memory recovery (all *P**fdr*>0.05). All associations were tested using LMM after regressing out age, gender, tumour grade, type of treatment and total lesion volume effects.

![Figure 4.](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2020/10/22/19008581/F4.medium.gif)

[Figure 4.](http://medrxiv.org/content/early/2020/10/22/19008581/F4)

Figure 4. Memory recovery as a function of lesion overlapping with DMN.
**(a)** Memory recovery as a function of tumour overlap with DMN. For a given participant, the (preoperative) tumour overlap value was used in all assessments, resulting in vertical lines between them. **(b)** Memory recovery a function of (follow-up) lesion overlap with DMN. None of the other Yeo networks showed any significant association that survived correction for multiple comparisons. Overlap values were defined as the number of voxels (equivalent to mm3) of the tumour mask that overlapped with each network. Lines link assessments that correspond to the same participant.

## DISCUSSION

In this study, we combined structural and NODDI diffusion MRI with normative brain network data from healthy participants and whole-brain analyses to determine whether memory recovery trajectories are affected by the impact that the tumour and its treatment have on specific circuits. Overall, our results suggest that the effect of a tumour on the brain and its consequences in memory depends on interactions with the DMN at both global levels. Acquiring high-quality longitudinal memory outcome data required nuanced study design. Using a tablet-based cognitive assessment app afforded specific advantages, such as avoiding reliance on having a trained neuropsychologist available (with associated time and financial costs), more easily facilitating follow-up screening, and collecting specific data that is difficult to acquire using traditional clinical interviews tests, such as accurate reaction times and interactive visuospatial paradigms. We observed a variety of memory recovery trajectories after surgery and during subsequent treatment that included: no change, postsurgical deterioration and then improvement, and postsurgical improvement and then deterioration. Although the mechanisms behind cognitive improvement after major surgery are not understood, these findings are not without precedent[Habets et al., 2014]. Nevertheless, given that pre- and postoperative assessments were performed in a relatively short period of time, we cannot discard the possibility of practice and learning effects on the tasks.

Tumour location is one of the most relevant features to be considered when estimating the cognitive risks of surgical resection. In our cohort, gliomas were mainly located in frontal, temporal and insular cortices on the left hemisphere. The higher prevalence of left-hemisphere tumours is because this study only included patients undergoing awake brain surgery. On the other hand, the higher prevalence of frontal, temporal and insular tumours is consistent with previous studies showing that low- and high-grade gliomas are relatively scarce in primary cortices and occipital lobes[De Witt Hamer et al., 2013]. Although several developmental, cytomyeloarchitectonic, neurochemical, metabolic, and functional reasons have been proposed, the mechanisms behind the preferential location of gliomas in frontal, temporal and insular cortices is still an ongoing debate[Mandal et al., 2020]. The presence of gliomas in secondary and association cortices that have been traditionally associated with cognitive processing[Goldman-Rakic, 1988] may be an important factor to understand cognitive deficits induced by the tumour and its treatment.

Beyond the impact that age, gender, tumour grade, tumour volume and type of treatment have on recovery, we hypothesized that memory functioning was associated with treatment-related interactions with brain networks that have been previously identified as fundamental for cognition and memory. In support of this hypothesis, we found significant associations between memory recovery and lesion overlap with the DMN. The DMN constitute a set of regions that are active during passive rest and it is considered the neuronal correlate of internally-generated self-referential cognition. Accordingly, a large body of literature has accumulated evidence about the role of the DMN in memory functioning, including semantic[Wirth et al., 2011], episodic[McCormick et al., 2013], autobiographical[Philippi et al., 2013] and working memory[Sambataro et al., 2010]. Accordingly, damage to the DMN has been associated with memory deficits in patients with brain lesions[Philippi et al., 2013] and traumatic brain injury[Sours et al., 2013]. To the best of our knowledge, there are no studies linking DMN disruption with memory recovery in brain tumour patients to date, but a decrease of functional connectivity in this network have been extensively reported in the literature. DMN functional connectivity is reduced in glioma patients when compared with controls in both the hemisphere ipsilateral and contralateral to the tumour, which is particularly prominent for tumours located on the left side of the brain[Ghumman et al., 2016].

Not surprisingly, we found a consistent Neurite Density decrease within the tumour. Given that Neurite Density derived from NODDI and Fractional Anisotropy are strongly correlated[Zhang et al., 2012], our results align with previous evidence found from DTI studies showing that glioblastomas have reduced Fractional Anisotropy compared with the rest of the brain[Sinha et al., 2002]. For its part, NODDI has shown higher sensitivity for glioma grade differentiation than other diffusion sequences[Vellmer et al., 2018]. Zhao et al. (2018) recently reported that NODDI in combination with patient age can predict glioma grade with a sensitivity and specificity of 92% and 89%, respectively[Zhao et al., 2018b]. Despite these promising findings, the potential of NODDI as a predictor of patient cognitive recovery has not yet been explored in the literature.

Despite peri-tumoural regions being identified on T1 images as unaffected by the experienced neurosurgeon that delineated the mask, and by the semi-automatic segmentation procedure, these regions had decreased neurite density compared with the contralateral hemisphere. Long-distance tumour effects have been observed with regard to functional connectivity and functional complexity[Hart et al., 2019]. Disrupted white matter integrity has been reported in DTI studies that show decreased FA in peri-tumoural regions[Holly et al., 2017] and white matter tracts[Miller et al., 2012]. Extratumoural pre- and post-operative Neurite Density within the Fronto-parietal and DMN had a protective factor in memory recovery that may be mediated by the prominent role of these network in this cognitive domain. From the connectomic perspective, cognitive performance of brain tumour patients has been associated with hub-related structural connectivity within the DMN and Frontoparietal Network in hemispheres contralateral to the tumour[Douw et al., 2019]. By using connectomic metrics derived from both DTI and fMRI, Liu et al. (2016) achieved a 75% accuracy when predicting survival of high-grade glioma patients[Liu et al., 2016]. However, further research is needed to better understand the mechanisms by which low- and high-grade glioma, and treatment disruption mediate memory and, more generally, cognitive recovery.

### Limitations

In the present work, brain networks were defined using normative data from healthy individuals[Yeo et al., 2011]. By using brain network templates in standard space, we neglected the potential spatial shift of brain functioning in our participants. Thus, our results should be interpreted as the potential of the tumour to disrupt the corresponding healthy network, not the actual tumour/lesion overlap with each participant’s network.

The heterogeneity of this condition limits the generalisability of the results. Treatment of the patients was decided solely on clinical criteria. As a consequence, 6 patients had only a surgical intervention while 11 additionally had different chemo-radiotherapy regimes. Chemo-radiotherapy has a dose-dependent effect on white matter structural integrity that has been associated with poor cognitive performance[Deprez et al., 2012]. Additionally, all our patients had the pre-operative imaging appearances of a diffuse glioma (for example, non-enhancing and without oedema or mass effect), subsequent pathological examination revealed a range of histological diagnoses, in keeping with previous studies demonstrating the diagnostic limitations of standard MRI sequences. Although the LMMs considered here regressed out the effect of age, gender, tumour grade, tumour volume and treatment, the individual contribution of these factors cannot be untangled here due to the limited sample size.

### Conclusion

Our findings highlight that memory recovery is not associated with the type of treatment, tumour volume nor the tumour location based on traditional lobes. Interestingly, Neurite Density is decreased beyond tumour boundary and high values within Frontoparietal and DMN are associated with better memory recovery. On the other hand, tumour and lesion overlapping with DMN has a negative impact in memory.. Taken together, these results reveal the potential of combining brain network analyses, normative connectomes, and advanced MRI sequences to better understand and predict the impact of brain tumours and their treatment on patients’ cognitive outcomes.

## Data Availability

Data is not currently available

## CONFLICT OF INTEREST

The authors report no competing interests.

## ACKNOWLEDGEMENTS

We thank all patients for generous involvement in the study. We also thank to Luca Villa, Jessica Ingham, Alexa Mcdonald for their contribution to the study. This research was supported by The Brain Tumour Charity and the Guarantors of Brain.

*   Received October 6, 2019.
*   Revision received October 22, 2020.
*   Accepted October 22, 2020.


*   © 2020, Posted by Cold Spring Harbor Laboratory

This pre-print is available under a Creative Commons License (Attribution-NonCommercial-NoDerivs 4.0 International), CC BY-NC-ND 4.0, as described at [http://creativecommons.org/licenses/by-nc-nd/4.0/](http://creativecommons.org/licenses/by-nc-nd/4.0/)

## REFERENCES

1.  Bressler SL, Menon V (2010): Large-scale brain networks in cognition: emerging methods and principles. Trends in Cognitive Sciences.
    
    

2.  Brown PD, Pugh S, Laack NN, Wefel JS, Khuntia D, Meyers C, Choucair A, Fox S, Suh JH, Roberge D, Kavadi V, Bentzen SM, Mehta MP, Watkins-Bruner D (2013): Memantine for the prevention of cognitive dysfunction in patients receiving whole-brain radiotherapy: A randomized, double-blind, placebo-controlled trial. Neuro Oncol.
    
    

3.  Caverzasi E, Papinutto N, Castellano A, Zhu AH, Scifo P, Riva M, Bello L, Falini A, Bharatha A, Henry RG (2016): Neurite Orientation Dispersion and Density Imaging Color Maps to Characterize Brain Diffusion in Neurologic Disorders. J Neuroimaging 26:494–498.
    
    

4.  Colgan N, Siow B, O’Callaghan JM, Harrison IF, Wells JA, Holmes HE, Ismail O, Richardson S, Alexander DC, Collins EC, Fisher EM, Johnson R, Schwarz AJ, Ahmed Z, O’Neill MJ, Murray TK, Zhang H, Lythgoe MF (2016): Application of neurite orientation dispersion and density imaging (NODDI) to a tau pathology model of Alzheimer’s disease. Neuroimage 125:739–744.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.neuroimage.2015.10.043&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=26505297&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F10%2F22%2F19008581.atom) 

5.  Deprez S, Amant F, Smeets A, Peeters R, Leemans A, Van Hecke W, Verhoeven JS, Christiaens MR, Vandenberghe J, Vandenbulcke M, Sunaert S (2012): Longitudinal assessment of chemotherapy-induced structural changes in cerebral white matter and its correlation with impaired cognitive functioning. J Clin Oncol 30:274–281.
    
    [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6MzoiamNvIjtzOjU6InJlc2lkIjtzOjg6IjMwLzMvMjc0IjtzOjQ6ImF0b20iO3M6Mzk6Ii9tZWRyeGl2L2Vhcmx5LzIwMjAvMTAvMjIvMTkwMDg1ODEuYXRvbSI7fXM6ODoiZnJhZ21lbnQiO3M6MDoiIjt9) 

6.  Douw L, Miller JJ, Steenwijk MD, Stufflebeam SM, Gerstner ER (2019): Altered structural hub connectivity and its clinical relevance in glioma. [http://dx.doi.org/10.1101/610618](http://dx.doi.org/10.1101/610618).
    
    

7.  Ghumman S, Fortin D, Noel-Lamy M, Cunnane SC, Whittingstall K (2016): Exploratory study of the effect of brain tumors on the default mode network. J Neurooncol 128:437– 444.
    
    

8.  Goldman-Rakic P (1988): Topography Of Cognition: Parallel Distributed Networks In Primate Association Cortex. Annu Rev Neurosci 11:137–156.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1146/annurev.ne.11.030188.001033&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=3284439&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F10%2F22%2F19008581.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=A1988M413100009&link_type=ISI) 

9.  Gordon EM, Laumann TO, Adeyemo B, Huckins JF, Kelley WM, Petersen SE (2016): Generation and Evaluation of a Cortical Area Parcellation from Resting-State Correlations. Cereb Cortex.
    
    

10. De Groot M, Douw L, Sizoo EM, Bosma I, Froklage FE, Heimans JJ, Postma TJ, Klein M, Reijneveld JC (2013): Levetiracetam improves verbal memory in high-grade glioma patients. Neuro Oncol.
    
    

11. Grussu F, Schneider T, Tur C, Yates RL, Tachrount M, Ianuş A, Yiannakas MC, Newcombe J, Zhang H, Alexander DC, DeLuca GC, Gandini Wheeler-Kingshott CAM (2017): Neurite dispersion: a new marker of multiple sclerosis spinal cord pathology?Ann Clin Transl Neurol.
    
    

12. Habets EJJ, Kloet A, Walchenbach R, Vecht CJ, Klein M, Taphoorn MJB (2014): Tumour and surgery effects on cognitive functioning in high-grade glioma patients. Acta Neurochir (Wien) 156:1451–1459.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1007/s00701-014-2115-8&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=24879620&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F10%2F22%2F19008581.atom) 

13. Hart MG, Romero-Garcia R, Price SJ, Suckling J (2019): Global Effects of Focal Brain Tumors on Functional Complexity and Network Robustness: A Prospective Cohort Study. Neurosurgery 84:1201–1213.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1093/neuros/nyy378&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=30137556&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F10%2F22%2F19008581.atom) 

14. Holly KS, Barker BJ, Murcia D, Bennett R, Kalakoti P, Ledbetter C, Gonzalez-Toledo E, Nanda A, Sun H (2017): High-grade Gliomas Exhibit Higher Peritumoral Fractional Anisotropy and Lower Mean Diffusivity than Intracranial Metastases. Front Surg 4:1– 10.
    
    

15. Huang Q, Zhang R, Hu X, Ding S, Qian J, Lei T, Cao X, Tao L, Qian Z, Liu H (2014): Disturbed small-world networks and neurocognitive function in frontal lobe low-grade glioma patients. PLoS One 9.
    
    

16. Klein M, Heimans JJ, Aaronson NK, Van Der Ploeg HM, Grit J, Muller M, Postma TJ, Mooij JJ, Boerman RH, Beute GN, Ossenkoppele GJ, Van Imhoff GW, Dekker AW, Jolles J, Slotman BJ, Struikmans H, Taphoorn MJB (2002): Effect of radiotherapy and other treatment-related factors on mid-term to long-term cognitive sequelae in low-grade gliomas: A comparative study. Lancet 360:1361–1368.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/S0140-6736(02)11398-5&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=12423981&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F10%2F22%2F19008581.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000179081600008&link_type=ISI) 

17. Kong APH, Lam PHP, Ho DWL, Lau JK, Humphreys GW, Riddoch J, Weekes B (2016): The Hong Kong version of the Oxford Cognitive Screen (HK-OCS): validation study for Cantonese-speaking chronic stroke survivors. Aging, Neuropsychol Cogn 23:530–548.
    
    

18. Liu L, Zhang H, Rekik I, Chen X, Wang Q, Shen D (2016): Outcome prediction for patient with high-grade gliomas from brain functional and structural networks. In: Lecture Notes in Computer Science (including subseries Lecture Notes in Artificial Intelligence and Lecture Notes in Bioinformatics) Vol. 9901 LNCS, pp 26–34.
    
    

19. Mandal AS, Romero-Garcia R, Hart M, Suckling J (2020): Genetic, Cellular, and Connectomic Characterization of the Adult Human Brain Regions Commonly Plagued by Glioma. medRxiv:2020.02.18.20018531. [http://medrxiv.org/content/early/2020/02/20/2020.02.18.20018531.1.abstract](http://medrxiv.org/content/early/2020/02/20/2020.02.18.20018531.1.abstract).
    
    

20. McCormick C, Quraan M, Cohn M, Valiante TA, McAndrews MP (2013): Default mode network connectivity indicates episodic memory capacity in mesial temporal lobe epilepsy. Epilepsia.
    
    

21. Meyers CA, Smith JA, Bezjak A, Mehta MP, Liebmann J, Illidge T, Kunkler I, Caudrelier JM, Eisenberg PD, Meerwaldt J, Siemers R, Carrie C, Gaspar LE, Curran W, Phan SC, Miller RA, Renschler MF (2004): Neurocognitive function and progression in patients with brain metastases treated with whole-brain radiation and motexafin gadolinium: Results of a randomized phase III trial. J Clin Oncol 22:157–165.
    
    [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6MzoiamNvIjtzOjU6InJlc2lkIjtzOjg6IjIyLzEvMTU3IjtzOjQ6ImF0b20iO3M6Mzk6Ii9tZWRyeGl2L2Vhcmx5LzIwMjAvMTAvMjIvMTkwMDg1ODEuYXRvbSI7fXM6ODoiZnJhZ21lbnQiO3M6MDoiIjt9) 

22. Miller P, Coope D, Thompson G, Jackson A, Herholz K (2012): Quantitative evaluation of white matter tract DTI parameter changes in gliomas using nonlinear registration.Neuroimage 60:2309–2315. [http://dx.doi.org/10.1016/j.neuroimage.2012.02.033](http://dx.doi.org/10.1016/j.neuroimage.2012.02.033).
    
    

23. Nazeri X, Chakravart M, Rotenberg DJ, Rajji TK, Rathi X, Michailovich O V., Voineskos AN (2015): Functional consequences of neurite orientation dispersion and density in humans across the adult lifespan. J Neurosci.
    
    

24. Owen M, Romero-Garcia R, McDonald A, Woodberry E, Coelho P, Morris RC, Price SJ, Santarius T, Suckling J, Erez Y, Hart MG (2020): Assessment of neuropsychological function during early treatment of diffuse glioma. medRxiv:2020.06.03.20119255. [http://medrxiv.org/content/early/2020/06/05/2020.06.03.20119255.abstract](http://medrxiv.org/content/early/2020/06/05/2020.06.03.20119255.abstract).
    
    

25. Philippi CL, Tranel D, Duff M, Rudrauf D (2013): Damage to the default mode network disrupts autobiographical memory retrieval. Soc Cogn Affect Neurosci.
    
    

26. Rae CL, Davies G, Garfinkel SN, Gabel MC, Dowell NG, Cercignani M, Seth AK, Greenwood KE, Medford N, Critchley HD (2017): Deficits in Neurite Density Underlie White Matter Structure Abnormalities in First-Episode Psychosis. Biol Psychiatry.
    
    

27. Rusbridge SL, Walmsley NC, Griffiths SB, Wilford PA, Rees JH (2013): Predicting outcomes of vocational rehabilitation in patients with brain tumours. Psychooncology.
    
    

28. Sagberg LM, Drewes C, Jakola AS, Solheim O (2017): Accuracy of operating neurosurgeons’ prediction of functional levels after intracranial tumor surgery. J Neurosurg 126:1173–1180.
    
    

29. Sambataro F, Murty VP, Callicott JH, Tan HY, Das S, Weinberger DR, Mattay VS (2010): Age-related alterations in default mode network: Impact on working memory performance. Neurobiol Aging.
    
    

30. Schültke E, Juurlink BHJ, Ataelmannan K, Laissue J, Blattmann H, Bräuer-Krisch E, Bravin A, Minczewska J, Crosbie J, Taherian H, Frangou E, Wysokinsky T, Chapman LD, Griebel R, Fourney D (2008): Memory and survival after microbeam radiation therapy. Eur J Radiol.
    
    

31. Sinha S, Bastin ME, Whittle IR, Wardlaw JM (2002): Diffusion tensor MR imaging of high-grade cerebral gliomas. Am J Neuroradiol 23:520–527.
    
    [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6NDoiYWpuciI7czo1OiJyZXNpZCI7czo4OiIyMy80LzUyMCI7czo0OiJhdG9tIjtzOjM5OiIvbWVkcnhpdi9lYXJseS8yMDIwLzEwLzIyLzE5MDA4NTgxLmF0b20iO31zOjg6ImZyYWdtZW50IjtzOjA6IiI7fQ==) 

32. Sours C, Zhuo J, Janowich J, Aarabi B, Shanmuganathan K, Gullapalli RP (2013): Default mode network interference in mild traumatic brain injury - A pilot resting state study. Brain Res.
    
    

33. Teixidor P, Gatignol P, Leroy M, Masuet-Aumatell C, Capelle L, Duffau H (2007): Assessment of verbal working memory before and after surgery for low-grade glioma. J Neurooncol.
    
    

34. Vellmer S, Tonoyan AS, Suter D, Pronin IN, Maximov II (2018): Validation of DWI pre-processing procedures for reliable differentiation between human brain gliomas. Z Med Phys 28:14–24.
    
    

35. Villanueva-Meyer JE, Mabray MC, Cha S (2017): Current clinical brain tumor imaging. Clin Neurosurg.
    
    

36. Wang N, Zhang J, Cofer G, Qi Y, Anderson RJ, White LE,  Allan Johnson G (2019): Neurite orientation dispersion and density imaging of mouse brain microstructure. Brain Struct Funct 224:1797–1813.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1007/s00429-019-01877-x&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=31006072&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F10%2F22%2F19008581.atom) 

37. Wirth M, Jann K, Dierks T, Federspiel A, Wiest R, Horn H (2011): Semantic memory involvement in the default mode network: A functional neuroimaging study using independent component analysis. Neuroimage.
    
    

38. De Witt Hamer PC, Hendriks EJ, Mandonnet E, Barkhof F, Zwinderman AH, Duffau H (2013): Resection Probability Maps for Quality Assessment of Glioma Surgery without Brain Location Bias. PLoS One 8.
    
    

39. Yeo BT, Krienen FM, Sepulcre J, Sabuncu MR, Lashkari D, Hollinshead M, Roffman JL, Smoller JW, Zöllei L, Polimeni JR, Fisch B, Liu H, Buckner RL (2011): The organization of the human cerebral cortex estimated by intrinsic functional connectivity. J Neurophysiol 106:1125–1165.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1152/jn.00338.2011&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=21653723&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F10%2F22%2F19008581.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000294775500007&link_type=ISI) 

40. Zhang H, Schneider T, Wheeler-Kingshott CA, Alexander DC (2012): NODDI: Practical in vivo neurite orientation dispersion and density imaging of the human brain. Neuroimage 61:1000–1016. [http://dx.doi.org/10.1016/j.neuroimage.2012.03.072](http://dx.doi.org/10.1016/j.neuroimage.2012.03.072).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.neuroimage.2012.03.072&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=22484410&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2020%2F10%2F22%2F19008581.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000305920600028&link_type=ISI) 

41. Zhao J,  Li J bin,  Wang J yan,  Wang Y liang,  Liu D wei,  Li X bei,  Song Y kun,  Tian Y su, Yan X,  Li Z hao,  He S fu,  Huang X long, Jiang L,  Yang Z yun,  Chu J ping (2018a): Quantitative analysis of neurite orientation dispersion and density imaging in grading gliomas and detecting IDH-1 gene mutation status. NeuroImage Clin 19:174–181. [https://doi.org/10.1016/j.nicl.2018.04.011](https://doi.org/10.1016/j.nicl.2018.04.011).
    
    

42. Zhao J,  Li J bin,  Wang J yan,  Wang Y liang,  Liu D wei,  Li X bei,  Song Y kun,  Tian Y su, Yan X,  Li Z hao,  He S fu,  Huang X long, Jiang L,  Yang Z yun,  Chu J ping (2018b): Quantitative analysis of neurite orientation dispersion and density imaging in grading gliomas and detecting IDH-1 gene mutation status. NeuroImage Clin.

 [1]: /embed/graphic-4.gif