Skip to main content

ORIGINAL RESEARCH article

Front. Oncol., 26 February 2021
Sec. Neuro-Oncology and Neurosurgical Oncology
This article is part of the Research Topic nTMS, Connectivity and Neuromodulation in Brain Tumor Patients View all 25 articles

Optimizing Adjuvant Stereotactic Radiotherapy of Motor-Eloquent Brain Metastases: Sparing the nTMS-Defined Motor Cortex and the Hippocampus

Yvonne Dzierma*&#x;Yvonne Dzierma1*†Michaela Schuermann&#x;Michaela Schuermann1†Patrick MelchiorPatrick Melchior1Frank NueskenFrank Nuesken1Joachim OertelJoachim Oertel2Christian RübeChristian Rübe1Philipp HendrixPhilipp Hendrix2
  • 1Department of Radiotherapy and Radiation Oncology, Saarland University Medical Centre, Homburg, Germany
  • 2Department of Neurosurgery, Saarland University Medical Centre and Saarland University Faculty of Medicine, Homburg, Germany

Brain metastases can effectively be treated with surgical resection and adjuvant stereotactic radiotherapy (SRT). Navigated transcranial magnetic stimulation (nTMS) has been used to non-invasively map the motor cortex prior to surgery of motor eloquent brain lesions. To date, few studies have reported the integration of such motor maps into radiotherapy planning. The hippocampus has been identified as an additional critical structure of radiation-induced deficits. The aim of this study is to assess the feasibility of selective dose reduction to both the nTMS-based motor cortex and the hippocampi in SRT of motor-eloquent brain metastases. Patients with motor-eloquent brain metastases undergoing surgical resection and adjuvant SRT between 07/2014 and 12/2018 were retrospectively analyzed. The radiotherapy treatment plans were retrieved from the treatment planning system (“original” plan). For each case, two intensity-modulated treatment plans were created: the “motor” plan aimed to reduce the dose to the motor cortex, the “motor & hipp” plan additionally reduce the dose to the hippocampus. The optimized plans were compared with the “original” plan regarding plan quality, planning target volume (PTV) coverage, and sparing of organs at risk (OAR). 69 plans were analyzed, all of which were clinically acceptable with no significant differences for PTV coverage. All OAR were protected according to standard protocols. Sparing of the nTMS motor map was feasible: mean dose 9.66 ± 5.97 Gy (original) to 6.32 ± 3.60 Gy (motor) and 6.49 ± 3.78 Gy (motor & hipp), p<0.001. In the “motor & hipp” plan, dose to the ipsilateral hippocampi could be significantly reduced (max 1.78 ± 1.44 Gy vs 2.49 ± 1.87 Gy in “original”, p = 0.003; mean 1.01 ± 0.92 Gy vs. 1.32 ± 1.07 Gy in “original”, p = 0.007). The study confirms the results from previous studies that inclusion of nTMS motor information into radiotherapy treatment planning is possible with a relatively straightforward workflow and can achieve reduced doses to the nTMS-defined motor area without compromising PTV coverage. Furthermore, we demonstrate the feasibility of selective dose reduction to the hippocampus at the same time. The clinical significance of these optimized plans yet remains to be determined. However, with no apparent disadvantages these optimized plans call for further and broader exploration.

Introduction

Brain metastases can effectively be treated with surgical resection and/or stereotactic radiotherapy (SRT). Lesions located within or adjacent to critical motor areas pose a challenge to both the neuro- and radiosurgeon. Increasing survival from effective interdisciplinary treatment regimens shifts attention to ameliorated secondary outcome rates such as improved motor function, cognitive function and quality of life. Preoperative neurosurgical planning and surgical resection itself predominantly aim to identify and preserve critical motor areas. In the past decade, navigated transcranial magnetic stimulation (nTMS) has been used to non-invasively map the motor cortex prior to surgery of motor eloquent brain lesions. Here, these preoperative motor maps appear to facilitate better resection rates while maintaining neurological function (15).

Whereas planning of SRT focusses on sparing distinct structures at risks, motor-eloquent areas have not routinely been integrated. Motor deficits have been observed to occur after high-dose Gamma Knife SRS to sites close to the motor cortex (6). Beyond direct motor-deficits, Pfeiffer et al. (7) have postulated a relationship between higher dose to the precentral gyrus and impaired verbal and working memory, attention and executive functions. To date, a small number of studies have reported the integration of motor maps into radiotherapy planning for dose reduction to motor areas, primarily in CyberKnife and GammaKnife treatment (811). Two studies (12, 13) have been carried out for linear accelerator (linac)-based radiotherapy, one of which considered patients with brain metastases. All studies showed that the inclusion of nTMS information into the radiotherapeutic planning workflow was possible and allowed for improved dose sparing of the motor-eloquent areas. However, additional confirmation from different medical centers and using different planning systems is still warranted. Furthermore, the hippocampus has been identified as an additional critical structure of radiation-induced cognitive deficits, which has not been considered in previous studies on nTMS in radiotherapy. In addition to the observed negative impact of whole-brain radiotherapy on cognitive outcome (14, 15), several studies have focussed on the hippocampus itself as one of the most critical structures for radiation injury owing to the ongoing neurogenesis in the subgranular zone of its dentate gyrus (1625). Higher dose to the hippocampus has particularly been associated with impaired verbal memory and higher executive functions (20, 25). Thus, hippocampal sparing in brain radiotherapy planning has also recently gained significant attention.

The aim of this study is hence to confirm the implementation of nTMS motor information into treatment planning of linac-based stereotactic radiotherapy of motor-eloquent brain metastases, and to assess the feasibility of selective dose reduction to both the nTMS-based motor cortex and the hippocampi at the same time.

Materials and Methods

Patients with motor-eloquent brain metastases undergoing nTMS-based surgical resection and adjuvant SRT between 07/2014 and 12/2018 were retrospectively analyzed. Metastases were regarded as motor-eloquent when infiltration of the precentral gyrus and/or pyramidal tract was presumed, or if the precentral gyrus and adjacent sulci could not be distinguished due to neuroanatomical distortion.

nTMS Mapping and Import Into the Radiotherapy Treatment Planning System

Patients received pre-operative magnetic resonance imaging (MRI) on a 1.5 T or 3 T scanner (Magnetom Symphony-TIM 1.5 T, Magnetom Skyra 3.0 T, Siemens, Erlangen, Germany) with contrast-enhanced T1-weighted MPRage in the axial direction (repetition time TR = 0.9 ms, echo time TE = 3.52 ms, flip-angle 15°, slice thickness 1 mm), on which dataset navigated transcranial magnetic stimulation was performed. The nTMS motor mapping was carried out using the Nexstim NBS system 4.3 (Nexstim Oy, Helsinki, Finlany) as previously described (5, 26, 27). In brief, the patients were sitting reclining in a chair with open eyes with surface electromyography electrodes attached to the muscles used for mapping (m. first dorsal interosseus, abductor pollicis brevis, abductor digiti minimi). The presumed location of the hand knob was used as a starting point, then varying the coil location and orientation to determine the resting motor threshold (RMT), defined as the lowest stimulus intensity which will elicit a 50 μV peak-to-peak amplitude motor evoked potential in five out of ten stimulations. The hand area was then mapped using 110% of the RMT and 0.25 Hz, holding the coil perpendicular to the precentral gyrus. Where possible, the lower extremity was mapped as delineated by the anterior tibial and plantar muscles, using 130% RMT intensity and a coil orientation perpendicular to the midline/falx. However, since lower extremity mapping was not available for all patients in this study collective, we only included patients who suffered from a motor-eloquent lesion where imaging predominantly appeared to demonstrate jeopardy of the upper extremity/hand area. Consequently, retrospective dose planning was carried out only for sparing of the upper extremity motor area.

In all cases the nTMS motor maps for the upper limb (i.e. hand area) were exported as an additional secondary dataset into the original radiotherapy treatment plan in the Philips Pinnacle treatment planning system (TPS) V16.2 for each patient. The secondary image was rigidly co-registered to the primary data set (planning computed tomography (CT)) and planning MRI (acquired post-operatively in both T1 MPRage and T2 flair weighting) based on a mutual information algorithm and then manually shifted until optimal correspondence was achieved. Correspondence was verified independently by two radiation physicists and/or radiation oncologists. The original clinical target volume (CTV), planning target volume (PTV) and organs at risk (OAR: lenses, bulbi, optic nerves, chiasm, cochleae, brainstem) for brain irradiation as defined on the planning-CT and planning-MRI were re-checked and the additional organs at risk (OAR’s) were contoured (i.e. nTMS-based motor cortex, basal ganglia, thalamus, hippocampus) based on the T1 weighted planning MRI sequences according to (28, 29). Each hippocampus was expanded by 5 mm into all directions to create the “hippocampus avoidance zone” for plan optimization.

Treatment Planning

For the original plans, different radiotherapy techniques were used: static beams with 3D-conformal radiotherapy (3D-CRT), intensity-modulated radiotherapy (IMRT) or non-coplanar arcs (30); 6 MV or flattening-filter-free 7 MV photons were used (31, 32).

For each metastasis, two re-optimized treatment plans were created in addition to the original clinically treated plan (“original”). The “motor” plan spared the hand motor areas delineated by nTMS. The “motor & hipp” plan aimed to reduce the dose to the nTMS-based motor cortex and the hippocampi. If the plan was originally planned by IMRT or by 3D-CRT, the same beam arrangements were used for the re-optimization; if the original plan used conformal non-coplanar arcs, a new plan was established by IMRT planning with between 8 and 13 beams. The maximum number of segments allowed was 35. The optimization objective for the motor cortex was to lower the maximum as well as the mean dose as much as possible without reducing the coverage of the PTV. Regarding the hippocampi, the optimization objectives were iteratively reduced to lower the maximum and the mean dose as much as possible without reducing the coverage of the PTV or burden the motor cortex again.

Optimization was performed using direct machine parameter optimization (DMPO) on a 0.2 cm dose grid; the final dose distributions were calculated using a collapsed cone (CC) convolution algorithm. All plans were revised by an experienced radiation oncologist and were in agreement with the general guidelines of the DEGRO Working Group on SRS for clinical stereotactic treatment (33) and the dose limits for sensitive brain structures based on the criteria of the Quantitative Analysis of Normal Tissue Effects in the Clinic (QUANTEC, 34).

Plan Evaluation

To evaluate plan differences, several measures of quality are considered.

The Paddick conformity index (CI) (3436)

CI=OR·UR=TVPIV2PIV·TV

is the product of the Paddick overdose ratio (OR) and underdose ratio (UR).

The overdose ratio

OR=TVPIVPIV

estimates the ratio of the PTV volume inside the prescribed 80 % isodose (TVPIV) to the total volume encompassing the 80% isodose (PIV = V80 %). This relates the covered PTV volume to the total volume irradiated with the prescribed encompassing dose.

The underdose ratio

UR=TVPIVTV

estimates the ratio of the PTV volume inside the 80 % isodose (TVPIV) to the PTV volume (TV). This relates the covered PTV volume to the total PTV volume.

The homogeneity index (HI) (3436)

HI=PTV1%PTV99%PTVmean

measures the PTV homogeneity by considering PTV1% as a measure of the maximum and PTV99% as a measure of the minimum dose in the PTV.

The gradient index (GI)

GI=V40%PIV

indicates the steepness of dose fall-off by comparing the volume of the prescription encompassing isodose to the volume (80%) of half this dose (V40%).

Two important values are V12 Gy and V10 Gy as well as their relative value to the total brain volume, as they correlate with the risk for necrosis in the case of stereotactic radiosurgery (37).

Besides these plan quality parameters, the doses in the critical sensitive brain structures based on the QUANTEC recommendations (38) as well as the PTV are determined. For the PTV, D01 % is given as a measure of the relative maximum; it is considered relative to the prescription dose in target point to estimate the amount of overdosage. D99% of the PTV is considered as a measure of the relative mimimum PTV dose, it is considered relative to the prescribed encompassing dose of 80% to estimate the amount of underdosage.

For the motor cortex, the intersections with the PTV, 90 %, 80 %, 70 %, 50 % and 20 % isodose are determined. Furthermore, D01% as a measure for the maximum dose and the mean dose are evaluated for the motor cortex and hippocampus as well as the other OARs.

Statistical Analyses

Dose-volume histogram (DVH) values were exported by an in-house Pinnacle script. Each OAR as well as the motor cortex and the PTV were saved in a CSV (comma separated variables) table. The reorganization into one table for each OAR and all calculations were performed with MATLAB R2019b. A normal distribution could not be presumed, so Wilcoxon’s signed-rank test for paired data was used. A 5 % level of significance was applied. For multiple comparisons (three scenarios), a Bonferroni correction was applied, in which p values below 0.0167 were considered statistically significant.

Results

A total of 52 patients were identified. Of these, 24 patients received stereotactic radiation therapy at this department. The remaining patients either received a different radiotherapy regimen or received radiotherapy elsewhere. For all but one patient, the 80% isodose of the prescribed maximum dose in the isocentre was required to encompass the PTV. One patient received a different fractionation with the 95% isodose level surrounding the target volume and was therefore excluded from the analysis. Another patient received irradiation for three neighboring target volumes which were jointly optimized – this patient was also excluded, since re-optimization would have involved all three target volumes with different prescriptions and isocenters. Treatment details for the remaining 22 patients included in the analysis are given in Table 1.

TABLE 1
www.frontiersin.org

Table 1 Patient characteristics.

For all 22 patients, the “motor” and “motor & hipp” treatment plans were considered acceptable for treatment. An example of the resulting isodose distributions is shown in Figure 1. Metrics for plan quality and dose to organs at risk are given in Table 2.

FIGURE 1
www.frontiersin.org

Figure 1 Example dose distributions for the original (left), motor (middle) and motor&hipp (right) plans for three different patients and planning scenarios. (A) patient treated originally by static beams, (B) patient treated originally with intensity-modulated radiotherapy (IMRT), (C) patient treated originally with non-coplanar arcs. The planning target volume (PTV) is delineated by the red filled contour, the nTMS-based motor cortex in skin color.

TABLE 2
www.frontiersin.org

Table 2 Plan evaluation using metrics for plan quality, planning target volume (PTV) coverage, and dose to organs at risk.

PTV and Organs at Risk in the Re-Optimized Plans

No significant differences between the “original”, “motor” and “motor & hipp” plans were observed for the coverage of the planning target volume as assessed by the conformity index, PTV minimum and maximum, and overdose ratio (Figure 2). There was a small, but statistically significant improvement in the underdose ratios of the “motor” and “motor & hipp” plans when compared with the original plans (but not with each other). The gradient index assessing dose fall-off outside the PTV was slightly worse in the newly-optimized plans, however, this did not affect the clinical acceptability of the plans.

FIGURE 2
www.frontiersin.org

Figure 2 Quality measures of planning target volume (PTV) coverage (upper panel) and dose to motor cortex and ipsilateral hippocampus. Statistically significant differences are denoted by asterisk.

All organs at risk could be well protected in both the “motor” and “motor & hipp” plans. The volume of the brain receiving a dose of 10 Gy or 12 Gy was slightly increased in the “motor & hipp” plans relative to the original and “motor” plans by ca. 3–4 cm³. However, this parameter is only relevant for stereotactic radiosurgery, i.e. very high single dose fractionation regimes. If we consider only patients receiving stereotactic radiosurgery (five cases), the three planning scenarios do not exhibit a significant difference in the volume of the 12 Gy isodose inside healthy brain, although there appears to be a trend toward somewhat increased volume (9.2 ± 2.8 cm³ and 9.5 ± 2.4 cm³ in the “motor” and “motor & hipp” plans, respectively, compared with 5.9 ± 1.9 cm³ in the original plan, p = 0.0625 and p = 0.3125, respectively). Simultaneously, sparing of the hippocampus resulted in reduced dose to the brainstem, thalami and basal ganglia. Some variations in OAR sparing were observed between the planning scenarios, however, all these organs received very little dose as compared with the planning objectives remained far below the clinically acceptable limits.

Sparing of the Motor Cortex

Sparing of the motor cortex could significantly be improved by both scenarios with nTMS information included, with mean dose to the motor cortex reduced from 9.66 ± 5.97 Gy (original) to 6.32 ± 3.60 Gy (motor) and 6.49 ± 3.78 Gy (motor & hipp), respectively (p<0.001 for both re-optimized plans vs. “original”). Regarding the spatial relationship of the motor cortex with the isodose levels (Figure 3), a reduction of overlap with all isodoses from 20 to 90% of the prescribed dose was observed relative to the original plans; however, this difference only reached statistical significance for the 70% isodose, becoming more pronounced for the lower isodose levels. The volume of the nTMS-derived motor map covered inside the 70% isodose could hence be reduced from an average of 7.4% (max 30.5%) to 4.8% (max 21.9%) (p = 0.015). A larger reduction was observed for the volume contained within the 50% isodose (22.3% vs. 10.1%, p = 0.003) and 20% isodose (61.8% vs. 34.7%, p < 0.001).

FIGURE 3
www.frontiersin.org

Figure 3 Motor cortex volume included inside planning target volume (PTV) and different relative isodose lines. Statistically significant differences are denoted by asterisk.

Since the patients in the collective received the original radiation treatment using three different technical approaches (3D-CRT, IMRT, non-coplanar arcs), a separate analysis of the dose to the motor cortex is performed for planning techniques relying on static beams with the original geometry (“beams” plans including 3D-CRT and IMRT) and non-coplanar “arcs” plans. Dose to the motor cortex in the “motor” and “motor & hipp” plans is compared with the original plans separately for the “beams” and “arcs” techniques in Figure 4. A considerably greater improvement could be attained by re-optimizing the “beams” plans than the “arcs” plans. This can also been visually confirmed in the isodose distribution (Figure 1). In all cases, the re-optimized plans show a more asymmetric behavior than the original plans since the isodoses are deformed so as to avoid the nTMS volume. While the prescription isodoses still encompass the PTV at least as well as in the original plan, the intermediate isodoses show a relatively strong asymmetry, where the gradient toward the motor cortex is much steeper than in all other directions. In particular, the original “arcs” plans already had a very steep gradient from the combination of many non-coplanar arcs, which cannot be achieved by a relatively simple co-planar IMRT beam configuration. Therefore, this planning scenario suffers from a relatively strong change in gradient (also reflected in GI), while achieving relatively little additional sparing of the motor cortex.

FIGURE 4
www.frontiersin.org

Figure 4 Average dose to the motor cortex for plans using non-coplanar conformal arcs vs. plans with static beams or intensity-modulated radiotherapy (IMRT). Statistically significant differences are denoted by asterisk.

For the “beams” plans, a strong positive correlation (Pearson correlation coefficient r = 0.903, p < 0.001) is observed between the distance PTV to nTMS-derived motor map and the relative change in nTMS mean dose, i.e., a greater improvement in nTMS mean dose is obtained the closer the motor cortex is located to the PTV.

Both the absolute dose reduction to the nTMS and its clinical relevance depend on the fractionation scheme. Figure 5 therefore shows the absolute difference in nTMS mean dose achieved for different fractionation regimes (one, three, or five fractions as detailed in Table 1). The amount of sparing achievable was higher in absolute dose for the more fractionated schedules. The single-fraction regimes included in our collective comprised three “arcs” plans. We have already seen that the “arcs” plans provided least improvement in nTMS dose by reoptimization, which will also influence the relatively low sparing achieved by the single-fraction regimes. Still, an average dose reduction (nTMS mean dose) of 0.88 Gy (both for “motor” and “motor & hipp”) plans could be achieved for this sub-collective, and for two patients around 2 Gy decrease were reached.

FIGURE 5
www.frontiersin.org

Figure 5 Absolute dose sparing in the navigated transcranial magnetic stimulation (nTMS) area (mean dose) achieved by the reoptimized plans as a function of fractionation regime.

Sparing of the Hippocampus

In the “motor” plans without hippocampus sparing, the dose to the ipsilateral hippocampus was increased over the original plan (in which also no attempt at hippocampus sparing was made): maximum dose 3.43 ± 2.47 Gy vs. 2.49 ± 1.87 Gy, p = 0.003; mean dose 1.82 ± 1.72 Gy vs. 1.32 ± 1.07 Gy, p = 0.002. After reoptimization in the “motor & hipp” plans, hippocampus ipsilateral mean and maximum dose (Figure 2) were significantly lower, even when compared with the original plan (maximum dose 1.78 ± 1.44 Gy in “motor & hipp” vs. 2.49 ± 1.87 Gy in “original”, p = 0.003; mean dose 1.01 ± 0.92 Gy in “motor & hipp” vs. 1.32 ± 1.07 Gy in “original”, p = 0.006). Also, hippocampus contralateral mean and maximum dose could be reduced by re-optimization, though not reaching statistical significance (max dose 0.56 ± 0.80 Gy in “motor & hipp” vs. 0.72 ± 1.11 Gy in “original”; mean dose 0.20 ± 0.21 Gy in “motor & hipp” vs. 0.32 ± 0.40 Gy in “original”).

Discussion

The study demonstrates the feasibility of selective dose reduction to both the nTMS-identified motor area and hippocampus for a comparison of treatment 69 plans, including three different planning approaches (static 3D-CRT, IMRT, non-coplanar conformal arcs). Even though this introduces some heterogeneity in the study collective, this wide selection of planning methods is representative of routine radiotherapy practice and differences could actually be identified in the amount of sparing achieved depending on the original planning technique.

The presented study cohort adds to hitherto only two reports investigating the implementation of nTMS into linac-based radiotherapy planning (12, 13), whereas the remaining studies investigated nTMS maps for CyberKnife (8, 9) or GammaKnife (10, 11) stereotactic radiosurgery. Linac-based radiotherapy, CyberKnife and GamaKnife differ in technical implementation, planning and dosimetry. As CyberKnife and GammaKnife availability is restricted to a small selection of specialized centers, whereas linear accelerators with stereotactic capability are much more wide-spread, an evaluation of motor cortex and hippocampus sparing achievable with linac stereotactic radiotherapy will be useful for a wider range of patients and treatment centers.

Diehl et al. (12) presented a planning study for 30 patients with high-grade gliomas treated by volumetric modulated arc therapy (VMAT) in a simultaneously integrated boost (SiB) concept. By including dose constraints to the nTMS motor maps (only areas outside the planning target volume), they could achieve a dose reduction of 12.8% (4.6 Gy), without compromising PTV coverage. In a collective more similar to the herein presented cohort, Schwendner et al. (13) presented a re-optimization of VMAT plans for a collective of 30 patients with supratentorial brain metastases. Again, the nTMS motor map not included in the PTV was spared in the re-optimized plans, resulting in a dose reduction of ca. 4.1 Gy (18.1%). Our study presents a similar collective of patients with brain metastases, but with a wider variation in original planning techniques, and with a larger range of distances of the motor cortex from the PTV. In line with those previous reports motor cortex sparing by about 3 Gy in mean dose (ca. 30%), with best results obtained for the static 3D-CRT and IMRT treatment plans was observed. Additionally, a correlation of dose reduction with proximity to the motor cortex was recognized.

A considerably greater improvement could be attained by re-optimizing the static 3D-CRT and IMRT plans than the plans using non-conformal arcs. However, this does not imply that non-conformal arcs should generally be the preferred planning scenario, since this is only applicable to relatively spherical and small target volumes. The dose to the motor cortex in the “beams” plans might be higher than in the “arcs” plans not primarily because of a presumably inferior planning technique, but possibly because a more complex PTV shape or volume, which might have been associated with a higher motor cortex dose and precluded the use of arcs. Yet, in the “arcs” plans, a small improvement could still be achieved by re-planning.

For the first time, the hippocampus is included as an additional organ at risk in the attempt to spare the motor map. Reducing radiation dose to the hippocampus and nTMS based motor map was feasible simultaneously without compromising the PTV or other organs at risk. Despite the inherent low dose to the ipsilateral hippocampi in the original plans (D01% 2.49 ± 1.87 Gy, mean dose 1.32 ± 1.07 Gy), a further dose reduction by 23%–28% could be achieved (D01% 1.78 ± 1.44 Gy, mean dose 1.01 ± 0.92 Gy). In particular, in the plans re-optimized for motorcortex sparing without hippocampus inclusion, the dose to the ipsilateral hippocampus was increased, so that the optimization of hippocampus dose reduced the dose from the “motor” plans by over 40%. Singular optimization of the motor cortex resulted in increased doses to the ipsilateral hippocampus avoidance zone and thalamus. A similar trend (without reaching statistical significance) could be found for the brainstem and ipsilateral basal ganglia. Although these changes would not have resulted in clinical rejection of the plans – partly because these structures are not routinely contoured and dosimetrically evaluated in clinical plans and partly because no clinically evaluated dose limits are yet available for these structures – the inclusion of hippocampus protection at the optimization stage could totally reverse this effect without loss of motor cortex sparing.

The external validity of the study is limited since nTMS mappings were not uniformly performed for the entire primary cortex. In a subset of included patients, lesion-specific mapping was performed to outline peri-lesional and critical motor areas to facilitate surgery. To increase internal validity, we included only nTMS-motor areas of the upper extremity into the optimization assessment. The omission of these structures means that for some patients at least the anatomic leg or face areas of the motor cortex may have received higher doses, since functional leg or face areas were not optimized in the plans. For future studies, comprehensive nTMS mapping of bilateral primary motor cortices including nTMS-based tractography is required to allow evaluation to which extent specific motor areas can be spared and whether sparing causes increased doses to other critical areas. Few authors have hitherto presented treatment plans with dose optimization for white-matter tracts (8, 9), and the clinical significance of this improvement is yet unclear.

Clinical Significance of the Dosimetric Improvements

The clinical significance of these dose-optimized plans for motor function, cognitive function and quality of life yet remains to be elucidated. The risk of brain radionecrosis in single-session radiosurgery is correlated with the volume receiving doses of 10 Gy and 12 Gy, and motor deficits are among the complications observed (39). However, a dose threshold has not yet been definitely established for functional impairment of the motor cortex. For SRS of the corticospinal tract, Maruyama et al. (40) proposed a 5% risk of motor complication when the volume receiving 20 and 25 Gy exceeded 58 mm³ and 21 mm², respectively. Based on a neuroanatomical target theory for a patient collective treated by whole-brain radiotherapy with boost in conventional fractionation (1.8–2Gy/fraction), Pfeiffer et al. (41) proposed that cognitive outcomes were affected by the volume of the left hippocampus receiving 10 Gy and by the volume of the left precentral gyrus receiving 40 Gy, among a range of other regions.

For the hippocampus, recent radiotherapy optimization studies have aimed at reducing maximum dose to 16–17 Gy and minimum dose to 9 Gy in hippocampus-sparing whole-brain radiotherapy (14, 25). However, while a significant correlation was observed between mean dose to the left hippocampus and cognitive deterioration (25) as well as hippocampal volume reduction (41), no dose cut-off was defined in these studies. Hippocampus volume loss was observed to be significant one year after high-dose radiotherapy (> 40 Gy), but not after low-dose radiotherapy (< 10 Gy) by Seibert et al. (42). The normal-tissue complication probability (NTCP) model proposed for the hippocampus by Gondi et al. (17) suggests that an EQD2 dose (dose equivalent to 2 Gy fractions) of greater than 7.3 Gy to 40% of the bilateral hippocampi may be associated with memory impairment as assessed by list-learning delayed verbal recall; however, this model could not be confirmed for low-grade glioma patients in the EORTC 22033 clinical trial (21). In the Hopkins Verbal Learning Test-Revised Delayed Recall, Ma et al. (19) reported 20% probability of decline for D100% hippocampus doses above 10.9 Gy and 50% probability for 59.3 Gy. Taken together, a relatively steep gradient appears to exist in this dose range, making any attempts of hippocampus sparing potentially relevant for cognitive performance.

In the present patient collective, doses to the hippocampus were already low in the original plans due to the location of the brain metastases near motor-eloquent cortical areas. However, even low doses to the hippocampus may influence neurogenesis and differentiation of the dendritic arbor, as has been shown in animal studies for single-fractions of 1 Gy (16) and fractionated low dose irradiation [5 to 20 fractions of 0.1 Gy (43)]. Although a reduction in hippocampus dose as attained in our study has not yet been proven to result in an observable change in clinical outcome, with no apparent disadvantages these optimized plans call for further and broader exploration. Prospective studies are required to assess whether dose sparing to the motor cortex and hippocampus can contribute to improved motor and/or cognitive outcomes and higher health-related quality of life.

In contrast to radiotherapy, surgical treatment of motor-eloquent lesions does not take brain areas into account that are deemed uncritical for the procedure or even contralateral. For radiotherapy, however, this information is essential due to the larger expansion of intermediate- and low-dose areas. From the neurosurgical perspective, uncritical areas may not be mapped or tractography not be performed unless explicitly requested for. The quality of the localization would be improved if post-operative mapping was performed at the same time with the planning MRI, so that a more precise fusion of the images unbiased by changes in tumor or edema distribution would become possible. On this basis, a clinical evaluation of the motor and cognitive performance could be correlated with the dose to the motor cortex and corticospinal tract, which will be a prerequisite for exploring dose-effect relationships. Furthermore, the distance of these structures from the high-dose volume may serve for risk stratification of radiation-induced effects. Rosenstock et al. (44). and Sollmann et al. (45). proposed nTMS-based risk-stratification models for neurosurgery. They correlated the proximity of motor-eloquent brain lesions to the motor cortex and corticospinal tract with the risk of motor deterioration due to treatment. Furthermore, the RMT was also reported as a marker of increased hazard for post-operative deficit. Possibly, a similar observation can be established for radiotherapy.

While exploring these possibilities, we would like to emphasize the importance of including in the considerations a critical structure such as the hippocampus, which might even be less well protected when a new optimization objective such as motor cortex sparing is added to the planning process without adequate hippocampus constraints. On-going neurogenesis in this structure renders it particularly sensitive to radiation-induced damage, so care should be taken to reduce hippocampus dose if this is possible. As we could show, sparing of the motor cortex and hippocampus is not mutually exclusive. Rather, adequate coverage of the PTV can be reconciled with combined dose reduction to the nTMS-defined hand area and hippocampus at no detriment in plan quality.

Conclusions

Selective dose reduction to the motor cortex and hippocampus is feasible without compromising PTV coverage or other organs at risk. The inclusion of the nTMS-based information on the motor cortex into plan optimization allowed for about 30% dose reduction (approximately 3 Gy in mean dose). However, singular optimization of the motor cortex causes an increase in dose to the hippocampus, thalamus and brain stem, which can be prevented by including hippocampus dose as an additional planning objective. In these plans, a further dose reduction to the hippocampus along with the motor cortex can be achieved, resulting in increased overall protection of these functional cortical areas.

Data Availability Statement

The raw data supporting the conclusions of this article will be made available by the authors, without undue reservation.

Author Contributions

YD and PH designed the study concept. PH and JO recruited the patient. JO was responsible for the surgical treatment of the patients and CR for the radiotherapy. nTMS mapping was performed by PH. Patient datasets were retrieved from the data base and retrospectively revised by YD. Import into the planning system and rigid registration with the planning CT was performed by YD and MS. PM and YD contoured the organs at risk. MS created the re-optimized radiotherapy treatment plans and performed the plan evaluation and statistical analysis. PM and FN evaluated the treatment plans. YD wrote the manuscript. MS prepared the figures and tables. All authors were involved in the interpretation of the results. All authors contributed to the article and approved the submitted version.

Conflict of Interest

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

References

1. Frey D, Schilt S, Strack V, Zdunczyk A, Rösler J, Niraula B, et al. Navigated transcranial magnetic stimulation improves the treatment outcome in patients with brain tumors in motor eloquent locations. Neuro Oncol (2014) 16(10):1365–72. doi: 10.1093/neuonc/nou110

PubMed Abstract | CrossRef Full Text | Google Scholar

2. Krieg SM, Sabih J, Bulubasova L, Obermueller T, Negwer C, Janssen I, et al. Preoperative motor mapping by navigated transcranial magnetic brain stimulation improves outcome for motor eloquent lesions. Neuro-Oncol (2014) 16(9):1274–82. doi: 10.1093/neuonc/nou007

PubMed Abstract | CrossRef Full Text | Google Scholar

3. Krieg SM, Picht T, Sollmann N, Bährend I, Ringel F, Nagarajan SS, et al. Resection of Motor Eloquent Metastases Aided by Preoperative nTMS-Based Motor Maps – Comparison of Two Observational Cohorts. Front Oncol (2016) 6:261. doi: 10.3389/fonc.2016.00261

PubMed Abstract | CrossRef Full Text | Google Scholar

4. Raffa G, Scibilia A, Conti A, Ricciardo G, Rizzo V, Morelli A, et al. The role of navigated transcranial magnetic stimulation for surgery of motor-eloquent brain tumors: a systematic review and meta-analysis. Clin Neurol Neurosurg (2019) 180:7–17. doi: 10.1016/j.clineuro.2019.03.003

PubMed Abstract | CrossRef Full Text | Google Scholar

5. Hendrix P, Dzierma Y, Burkhardt BW, Simgen A, Wagenpfeil G, Griessenauer CJ, et al. Preoperative navigated transcranial magnetic stimulation improves gross total resection rates in patients with motor eloquent high-grade gliomas: A matched cohort study. Neurosurgery (2020) 0:1–10. doi: 10.1093/neuros/nyaa486

CrossRef Full Text | Google Scholar

6. Park C-Z, Choi H-Y, Lee S-R, Roh TH, Seo M-R, Kim S-H. Neurological Change after Gamma Knife Radiosurgery for Brain Metastases Involving the Motor Cortex. Brain Tum Res Treat (2016) 4(2):111–5. doi: 10.14791/btrt.2016.4.2.111

CrossRef Full Text | Google Scholar

7. Pfeiffer AM, Leyrer CM, Greene-Schloesser DM, Shing E, Kearns WT, Hinson WH, et al. Neuroanatomical target theory as a predictive model for radiation-induced cognitive decline. Neurology (2013) 80:747–53. doi: 10.1212/WNL.0b013e318283bb0a

PubMed Abstract | CrossRef Full Text | Google Scholar

8. Picht T, Schilt S, Frey D, Vajkoczy P, Kufeld M. Integration of navigated brain stimulation data into radiosurgical planning: potential benefits and dangers. Acta Neurochir (2014) 156:1125–33. doi: 10.1007/s00701-014-2079-8

PubMed Abstract | CrossRef Full Text | Google Scholar

9. Conti A, Pontoriero A, Ricciardi GK, Granata F, Vinci S, Angileri FF, et al. Integration of functional neuroimaging in CyberKnife radiosurgery: feasibility and dosimetric results. Neurosurg Focus (2013) 34(4):E5. doi: 10.3171/2013.2.FOCUS12414

PubMed Abstract | CrossRef Full Text | Google Scholar

10. Islam M, Cooray G, Benmakhlouf H, Hatiboglu M, Sinclair G. Integrating navigated transcranial magnetic stimulation motor mapping in hypofractionated and single-dose gamma knife radiosurgery: A two-patient case series and a review of literature. Surg Neurol Int (2020) 11(29):1–11. doi: 10.25259/SNI_406_2019

PubMed Abstract | CrossRef Full Text | Google Scholar

11. Tokarev AS, Rak VA, Sinkin MV, Evdokimova EL, Stepanov VN, Koynash GV, et al. Appliance of Navigated Transcranial Magnetic Stimulation in Radiosurgery of Brain Metastases. J Clin Neurophysiol (2020) 37(1):50–5. doi: 10.1097/WNP.0000000000000621

PubMed Abstract | CrossRef Full Text | Google Scholar

12. Diehl CD, Schwendner MJ, Sollmann N, Oechsner M, Meyer B, Combs SE, et al. Application of presurgical navigated transcranial magnetic stimulation motor mapping for adjuvant radiotherapy planning in patients with high-grade gliomas. Radiother Oncol (2019) 138:30–7. doi: 10.1016/j.radonc.2019.04.029

PubMed Abstract | CrossRef Full Text | Google Scholar

13. Schwendner MJ, Sollmann N, Diehl CD, Oechsner M, Meyer B, Krieg SM, et al. The Role of Navigated Transcranial Magnetic Stimulation Motor Mapping in Adjuvant Radiotherapy Planning in Patients With Supratentorial Brain Metastases. Front Oncol (2018) 8:424. doi: 10.3389/fonc.2018.00424

PubMed Abstract | CrossRef Full Text | Google Scholar

14. Brown PD, Jaeckle K, Ballmann KV, Farace E, Cerhan JH, Anderson K, et al. Effect of Radiosurgery Alone vs Radiosurgery with Whole Brain Radiation Therapy on Cognitive Function in Patients with 1 to 3 Brain Metastases – A Randomized Clinical Trial. JAMA (2016) 316(4):401–9. doi: 10.1001/jama.2016.9839

PubMed Abstract | CrossRef Full Text | Google Scholar

15. Chang EL, Wegel JS, Hess KR, Allen PK, Lang FF, Kornguth DG, et al. Neurocognition in patients with brain metastases treated with radiosurgery or radiosurgery plus whole-brain irradiation: a randomies controlled trial. Lancet (2009) 10:1037–44. doi: 10.1016/S1470-2045(09)70263-3

CrossRef Full Text | Google Scholar

16. Parihar VK, Limoli CL. Cranial irradiation compromises neuronal architecture in the hippocampus. PNAS (2013) 110(31):12822–7. doi: 10.1073/pnas.1307301110

PubMed Abstract | CrossRef Full Text | Google Scholar

17. Gondi V, Hermann BP, Mehta MP, Tomé WA. Hippocampal Dosimetry Predicts Neurocognitive Function Impairment After Fractionated Stereotactic Radiotherapy for Benign or Low-Grade Adult Brain Tumors. Int J Radiat Oncol Biol Phys (2013) 85(2):348–54. doi: 10.1016/j.ijrobp.2012.11.031

PubMed Abstract | CrossRef Full Text | Google Scholar

18. Gondi V, Pugh SL, Tome WA, Caine C, Corn B, Kanner A, et al. Preservation of Memory With Conformal Avoidance of the Hippocampal neural Stem-Cell Compartment During Whole-Brain Radiotherapy for Brain Metastases (RTOG 0933): A Phase II Multi-Institutional Trial. J Clin Onc (2014) 32(34):3810–6. doi: 10.1200/JCO.2014.57.2909

CrossRef Full Text | Google Scholar

19. Ma TM, Grimm J, McIntyre R, Anderson-Keightly H, Kleinberg LR, Hales RK, et al. A prospective evaluation of hippocampal radiation dose volume effects and memory deficits following cranial irradiation. Radiother Oncol (2017) 125:234–40. doi: 10.1016/j.radonc.2017.09.035

PubMed Abstract | CrossRef Full Text | Google Scholar

20. Tsai P-F, Yang C-C, Chuang C-C, Huang T-Y, Wu Y-M, Pai P-C, et al. Hippocampal dosimetry correlates with the change in neurocognitive function after hippocampal sparing during whole brain radiotherapy: a prospective study. Radiat Oncol (2015) 10:253. doi: 10.1186/s13014-015-0562-x

PubMed Abstract | CrossRef Full Text | Google Scholar

21. Jaspers J, Mèndez Romero A, Hoogeman MS, van den Bent M, Wiggenraad RGJ, Taphoorn MJB, et al. Evaluation of the Hippocampal Normal Tissue Complication Model in a Prospective Cohort of Low Grade Glioma Patients – An Analysis Within the EORTC 22033 Clinical Trial. Font Oncol (2019) 9:991. doi: 10.3389/fonc.2019.00991

CrossRef Full Text | Google Scholar

22. Birer SR, Olson AC, Adamson J, Hood R, Susen M, Kim G, et al. Hippocampal dose from stereotactic radiosurgery for 4 to 10 brain metastases: Risk factors, feasibility of dose reduction via re-optimization, and patient outcomes. Med Dosim (2017) 42(4):310–6. doi: 10.1016/j.meddos.2017.06.007

PubMed Abstract | CrossRef Full Text | Google Scholar

23. Okoukoni C, McTyre ER, Ayala Peacock DN, Pfeiffer AM, Strowd R, Cramer C, et al. Hippocampal dose volume histogram predicts Hopkins Verbal Learning Test scores after brain irradiation. Adv Radiat Oncol (2017) 2:624–9. doi: 10.1016/j.adro.2017.08.013

PubMed Abstract | CrossRef Full Text | Google Scholar

24. Brown PD, Gondy V, Pugh S, Tome WA, Wefel JS, Armstrong TS, et al. Hippocampal Avoidance During Whole-Brain Radiotherapy Plus Memantine for Patients With Brain Metastases: Phase III Trial NRG Oncology CC001. J Clin Oncol (2020) 38(10):1019. doi: 10.1200/JCO.19.02767

PubMed Abstract | CrossRef Full Text | Google Scholar

25. Kim KS, Wee CW, Seok J-Y, Hong JW, Chung J-B, Eom K-Y, et al. Hippocampus-sparing radiotherapy using volumetric modulated arc therapy (VMAT) to the primary brain tumor: the result of dosimetric study and neurocognitive function assessment. Radiat Oncol (2018) 13:29. doi: 10.1186/s13014-018-0975-4

PubMed Abstract | CrossRef Full Text | Google Scholar

26. Hendrix P, Senger S, Griessenauer CJ, Simgen A, Schwerdtfeger K, Oertel J. Preoperative Navigated Transcranial Magnetic Stimulation in Patients with Motor Eloquent Lesions with Emphasis on Metastasis. Clin Anat (2016) 29:925–31. doi: 10.1002/ca.22765

PubMed Abstract | CrossRef Full Text | Google Scholar

27. Hendrix P, Senger S, Griessenauer CJ, Simgen A, Linsler S, Oertel J. Preoperative navigated transcranial magnetic stimulation and tractography in transparietal approach to the trigone of the lateral ventricle. J Clin Neurosci (2017) 41:154–61. doi: 10.1016/j.jocn.2017.02.029

PubMed Abstract | CrossRef Full Text | Google Scholar

28. Eekers DBP, Ven L, Roelofs E, Postma A, Alapetite C, NG B, et al. The EPTN consensus-based atlas for CT- and MR-based contouring in neuro-oncology, 2018. Radiother Oncol 128:37–43. doi: 10.1016/j.radonc.2017.12.013

PubMed Abstract | CrossRef Full Text | Google Scholar

29. Scoccianti S, Detti B, Gadda D, Greto D, Furfaro I, Fiammetta M, et al. Organs at risk in the brain and their dose-constraints in adults and in children: A radiation oncologist’s guide for delineation in everyday practice. Radiother Oncol (2015) 2015:114:230–238. doi: 10.1016/j.radonc.2015.01.016

CrossRef Full Text | Google Scholar

30. Dzierma Y, Nuesken FG, Palm J, Licht NP, Ruebe C. Planning study and dose measurements of intracranial stereotactic radiation surgery with a flattening filter-free linac. Pract Radiat Oncol (2014) 4(2):e109–116. doi: 10.1016/j.prro.2013.04.004

PubMed Abstract | CrossRef Full Text | Google Scholar

31. Dzierma Y, Licht N, Nuesken F, Ruebe C. Beam properties and stability of a flattening-filter free 7 MV beam – an overview. Med Phys (2012) 39(5):2595–602. doi: 10.1118/1.3703835

PubMed Abstract | CrossRef Full Text | Google Scholar

32. Sadrollahi A, Nuesken F, Licht N, Rübe C, Dzierma Y. Monte-Carlo simulation of the Siemens Artiste linear accelerator flat 6 MV and flattening-filter-free 7 MV beam line. PloS One (2019) 14(1):e0210069. doi: 10.1371/journal.pone.0210069

PubMed Abstract | CrossRef Full Text | Google Scholar

33. Kocher M, Wittig A, Piroth MD, Treuer H, Seegenschmiedt H, Ruge M, et al. Stereotactic radiosurgery for brain metastases. A report of the DEGRO Working Group on Stereotactic Radiotherapy. Strahlenther onkol (2014) 190(6):521–32. doi: 10.1007/s00066-014-0648-7

PubMed Abstract | CrossRef Full Text | Google Scholar

34. Shaw E, Kline R, Gillin M, Souhami L, Hirschfeld A, Dinapoli R, et al. Radiation Therapy Oncology Group: radiosurgery quality assurance guidelines. Int J Radiat Oncol Biol Phys (1993) 27:1231–9. doi: 10.1016/0360-3016(93)90548-A

PubMed Abstract | CrossRef Full Text | Google Scholar

35. Paddick I. A simple scoring ratio to index the conformity of radiosurgical treatment plans. Technical note. J Neurosurg (2000) 93(Suppl 3):219–22. doi: 10.3171/jns.2000.93.supplement_3.0219

PubMed Abstract | CrossRef Full Text | Google Scholar

36. Paddick I, Lippitz B. A simple dose gradient measurement tool to complement the conformity index. J Neurosurg (2006) 105(Suppl):194–201. doi: 10.3171/sup.2006.105.7.194

PubMed Abstract | CrossRef Full Text | Google Scholar

37. Lawrence YR, Allen Li X, el Naqa I, Hahn CA, Marks LB, Merchant TE, et al. Radiation dose-volume effects in the brain. Int J Radiat Oncol Biol Phys (2010) 76(3):20–7. doi: 10.1016/j.ijrobp.2009.02.091

CrossRef Full Text | Google Scholar

38. Marks LB, Yorke ED, Jackson A, Ten Haken RK, Constine LS, Eisbruch A, et al. Use of normal tissue complication probability models in the clinic. Int J Radiat Oncol Biol Phys (2010) 76(3):S10–9. doi: 10.1016/j.ijrobp.2009.07.1754

PubMed Abstract | CrossRef Full Text | Google Scholar

39. Minniti G, Clarke E, Lanzetta G, Osti MF, Trasimeni G, Bozzao A, et al. Stereotactic radiosurgery for brain metastases: analysis of outcome and risk of brain radionecrosis. Radiat Oncol (2011) 6:48. doi: 10.1186/1748-717X-6-48

PubMed Abstract | CrossRef Full Text | Google Scholar

40. Maruyama K, Kamada K, Ota T, Koga T, Itoh D, Ino K, et al. Tolerance of pyramidal tract to gamma knife radiosurgery based on diffusion-tensor tractography. Int J Radiat Oncol Biol Phys (2008) 70(5):1330–5. doi: 10.1016/j.ijrobp.2007.08.010

PubMed Abstract | CrossRef Full Text | Google Scholar

41. Takeshita Y, Watanabe K, Kakeda S, Hamamura T, Sugimoto K, Masaki H, et al. Early volume reduction of the hippocampus after whole-brain radiation therapy: an automated brain structure segmentation study. Jap J Radiol (2020) 38:118–25. doi: 10.1007/s11604-019-00895-3

CrossRef Full Text | Google Scholar

42. Seibert TM, Karunamuni R, Bartsch H, Kaifi S, Krishnan AP, Dalia Y, et al. Radiation Dose-Dependent Hippocampal Atrophy Detected With Longitudinal Volumetric Magnetic Resonance Imaging. Int J Radiat Oncol Biol Physis (2017) 97(2):263–9. doi: 10.1016/j.ijrobp.2016.10.035

CrossRef Full Text | Google Scholar

43. Schmal Z, Isermann A, Hladik D, von Toerne C, Tapio S, Rübe CE. DNA damage accumulation during fractionated low-dose radiation compromises hippocampal neurogenesis. Radiother Oncol (2019) 137:45–54. doi: 10.1016/j.radonc.2019.04.021

PubMed Abstract | CrossRef Full Text | Google Scholar

44. Rosenstock T, Grittner U, Acker G, Schwarzer V, Kulchytska N, Vajkoczy P, et al. Risk stratification in motor area-related glioma surgery based on navigated transcranial magnetic stimulation data. J Neurosurg (2017) 126:1227–37. doi: 10.3171/2016.4.JNS152896

PubMed Abstract | CrossRef Full Text | Google Scholar

45. Sollmann N, Wildschuetz N, Kelm A, Conway N, Moser T, Bulubas L, et al. Associations between clinical outcome and navigated transcranial magnetic stimulation characteristics in patients with motor-eloquent brain lesions: a combined navigated transcranial magnetic stimulation-diffusion tensor imaging fiber tracking approach. J Neurosurg (2018) 128:800–10. doi: 10.3171/2016.11.JNS162322

PubMed Abstract | CrossRef Full Text | Google Scholar

Keywords: nTMS mapping, motor cortex, hippocampus sparing, treatment planning, functional optimization, IMRT, stereotactic radiation, brain metastases

Citation: Dzierma Y, Schuermann M, Melchior P, Nuesken F, Oertel J, Rübe C and Hendrix P (2021) Optimizing Adjuvant Stereotactic Radiotherapy of Motor-Eloquent Brain Metastases: Sparing the nTMS-Defined Motor Cortex and the Hippocampus. Front. Oncol. 11:628007. doi: 10.3389/fonc.2021.628007

Received: 10 November 2020; Accepted: 06 January 2021;
Published: 26 February 2021.

Edited by:

Thomas Picht, Charité – Universitätsmedizin Berlin, Germany

Reviewed by:

Tizian Rosenstock, Charité – Universitätsmedizin Berlin, Germany
Sandro M. Krieg, Technical University of Munich, Germany
Ina Bährend, Charité-Klinik für Neurochirurgie, Germany

Copyright © 2021 Dzierma, Schuermann, Melchior, Nuesken, Oertel, Rübe and Hendrix. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.

*Correspondence: Yvonne Dzierma, WXZvbm5lLmR6aWVybWFAdWtzLmV1

†These authors share first authorship

Disclaimer: All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.