- Open Access
CyberKnife® radiosurgery in the treatment of complex skull base tumors: analysis of treatment planning parameters
Radiation Oncologyvolume 1, Article number: 46 (2006)
Tumors of the skull base pose unique challenges to radiosurgical treatment because of their irregular shapes, proximity to critical structures and variable tumor volumes. In this study, we investigate whether acceptable treatment plans with excellent conformity and homogeneity can be generated for complex skull base tumors using the Cyberknife® radiosurgical system.
At Georgetown University Hospital from March 2002 through May 2005, the CyberKnife® was used to treat 80 patients with 82 base of skull lesions. Tumors were classified as simple or complex based on their proximity to adjacent critical structures. All planning and treatments were performed by the same radiosurgery team with the goal of minimizing dosage to adjacent critical structures and maximizing target coverage. Treatments were fractionated to allow for safer delivery of radiation to both large tumors and tumors in close proximity to critical structures.
The CyberKnife® treatment planning system was capable of generating highly conformal and homogeneous plans for complex skull base tumors. The treatment planning parameters did not significantly vary between spherical and non-spherical target volumes. The treatment parameters obtained from the plans of the complex base of skull group, including new conformity index, homogeneity index and percentage tumor coverage, were not significantly different from those of the simple group.
Our data indicate that CyberKnife® treatment plans with excellent homogeneity, conformity and percent target coverage can be obtained for complex skull base tumors. Longer follow-up will be required to determine the safety and efficacy of fractionated treatment of these lesions with this radiosurgical system.
Lesions of the base of skull are typically slow growing, but potentially morbid tumors . They rarely metastasize making local control the primary determinant of long-term survival . Although surgical resection may still be the treatment "gold-standard" [3, 4], radiosurgery is an appropriate treatment option for many patients . However, single-fraction radiosurgical treatment may be difficult because of the potentially large size and irregular shapes of these tumors. Their proximity to critical structures also leads to a risk of radiation-induced, long-term, neurological complication .
The CyberKnife® is a newly FDA approved radiosurgical devise for the treatment of brain lesions. Unlike the LINAC and Gamma Knife, the CyberKnife® is an image-guided, frameless radiosurgery system. Treatment is delivered by a linear accelerator mounted on a flexible robotic arm. Several-hundred treatment beams are chosen out of a repertoire of greater than one thousand possible beam directions using inverse treatment planning. These beams are delivered in a non-isocentric manner via circular collimators of varying size without intensity modulation. Non-isocentric treatment allows for simultaneous irradiation of multiple lesions. The lack of a requirement for the use of a head-frame allows for staged treatment. Since the planning system has access to a large number of potential non-isocentric beams, the CyberKnife® should theoretically be able to deliver a highly conformal, uniform dose with steep dose gradients . Therefore, treatment with the CyberKnife® radiosurgical system should minimize toxicity to surrounding structures. When compared to commonly used radiosurgical devices, such as the Gamma Knife, linear-accelerator based stereotactic radiosurgery systems with multiple arcs (LINAC), or intensity modulated radiation therapy, dosimetric studies of ellipsoid phantoms have shown that the CyberKnife® radiosurgical system has the best homogeneity within the target volume and comparable conformity .
A dose-volume histogram (DVH) is the tool most commonly used to compare radiosurgical plans. Unfortunately, the large volume of data in these histograms does not allow for simple differentiation between multiple plans and systems [9, 10]. Thus, an effort has been made to determine simple measurements for plan optimization. A conformity index is a single measure of how well the treatment dose distribution of a specific radiation treatment plan conforms to the size and shape of the target volume. In general, the conformity index of a given radiosurgical plan is dependent on target shape , target volume , collimator size , type of collimation (circular vs multileaf) and radiosurgical delivery system.
The new conformity index (NCI) and homogeneity index (HI) allow for the quick and simple comparison of different radiosurgical treatment plans, whether within the same radiosurgical system, or across diverse systems such as between the LINAC and Gamma Knife . Conformity indices have been reported in the literature, ranging from 1.0 to 3.0 for varying radiosurgical systems [14–18]. Typically, multiple iso-center plans generated with the Gamma Knife have homogeneity indices (HI) of 2.0 to 3.0 while the LINAC plans generate homogeneity indices (HI) of 1.0 to 1.2 . The significance of these differences between systems is controversial.
We determined the NCI and HI for the first 82 base of skull lesions treated at Georgetown University Hospital using the CyberKnife® radiosurgical system (Accuray, Sunnyvale, CA). We undertook this study to determine the effect of target shape, target volume and proximity to critical structures on radiosurgical treatment parameters. This is the first study that we are aware of that investigates these parameters in patients treated with the CyberKnife® radiosurgery system.
Patients and methods
We performed a retrospective review of 262 patients with intracranial tumors, who were treated with CyberKnife® stereotactic radiosurgery at Georgetown University Hospital between March 2002 and May 2005. Eighty-one patients were classified to have tumors of the skull base resulting in a total of 84 treated lesions. Thirty-three percent of these lesions had been previously irradiated. One patient was excluded from analysis because two tumor volumes were treated simultaneously making it impossible to calculate indices for each individual lesion.
Of the remaining lesions, 46 were categorized into the complex, skull base tumor group. A complex skull base tumor was defined as one that completely encircles, partially circumscribes, or directly contacts the brainstem, optic chiasm, hypophysis, or cranial nerves with meaningful remaining function. This complex tumor group consisted of 18 men and 26 women, with a median age of 53 (range: 29 – 88). These tumors were further categorized by histopathology as follows: 21 meningiomas, 6 metastatic tumors, 8 schwannomas, 7 pituitary adenomas, 1 chordoma, 2 sarcomas, and 1 glioma. The median tumor size was 7.27 cc (range: 0.62 – 98.3 cc) (Table 1 &2).
The data from the group with complex skull base tumors were compared with data from two control groups. The first group consisted of 36 patients with skull base tumors that were classified as simple. Although still located in the region of the skull base, tumors in this group had at least a 2 mm separation from the nearest critical structure. This group consisted of 16 men and 20 women, with a median age of 55 (range: 17 – 18). These tumors were also categorized by histopathology as follows: 5 meningiomas, 13 metastatic tumors, 10 schwannomas, 3 pituitary adenomas, 1 chordoma, 2 sarcomas, and 2 malignant gliomas. The median tumor size in this group was 8.83 cc (range: 0.19 – 206.3 cc) (Table 1 &2).
A second control group used for comparison consisted of 43 patients with metastatic tumors of the cerebral and cerebellar hemispheres. These lesions represented volumes that were spherical, with smooth borders, and relatively distant from critical neurovascular structures. This group consisted of 23 men and 20 women, with a median age of 58 (range: 21 – 85). These tumors were further categorized by histopatholgy as 33 metastatic carcinomas and 10 melanomas. The median tumor size in this group was 1.43 cc (range: 0.12 – 66 cc) (Table 1 &2).
Radiosurgical treatment planning
The basic technical aspects of CyberKnife® radiosurgery for cranial tumors have been described in detail (CyberKnife® Radiosurgery, A Practical Guide). Briefly, the patient was placed in a supine position on a vacuum bag and a malleable thermoplastic mask was molded to the head and attached to the head support. Thin-sliced (1.25 mm) high-resolution CT images were obtained through the region of interest with the patient in the treatment position. Target volumes and critical structures were delineated by the treating neurosurgeon. The treating neurosurgeon and radiation oncologist determined the minimal tumor margin dose of the target volume and the treatment isodose. This discussion was influenced by various factors, including previous radiation to the area, tumor volume, and extent of contact and compression of critical neurological structures. In most cases, the dose was prescribed to the isodose surface that encompassed the margin of the tumor. Twelve collimator sizes are available with the CyberKnife® radiosurgical system ranging from 5 mm to 60 mm. In general, a collimator size less than the maximum length of the prescribed target volume (PTV) was chosen for treatment planning . An inverse planning method with non-isocenteric technique was used for all cases. The treating physician and physicist input the specific treatment criteria, limiting the maximum dose to structures such as the optic chiasm and brainstem. The majority of the treatments were given in five fractions. In general, for non-previously treated cases, treatment plans were deemed acceptable if the maximum dose to critical structures was less than 2000 cGy in five fractions. Non-anatomical dose constraint structures were commonly incorporated to aid the optimization process in minimizing the dose to critical structures. The planning software calculated the optimal solution for treatment. The DVH of each plan was evaluated until an acceptable plan was generated.
Treatment planning parameters
Target volume was defined as the volume contoured on the planning CT scan by the treating neurosurgeon. No margin was added to the target volume. In this study, there was no limit set on the treatable target volumes.
The homogeneity index (HI) describes the uniformity of dose within a treated target volume, and is directly calculated from the prescription isodose line chosen to cover the margin of the tumor:
New Conformity Index
The new conformity index (NCI) as formulated by Paddick , and modified by Nakamura  describes the degree to which the prescribed isodose volume conforms to the shape and size of the target volume. It also takes into account avoidance of surrounding normal tissue.
Percent Target Coverage
PTC = The percentage of the target volume covered by the prescription isodose.
Radiosurgical treatment delivery
Image-guided radiosurgery was employed to eliminate the need for stereotactic frame fixation. Using computed tomography planning, target volume locations were related to radiographic landmarks of the cranium. With the assumption that the target position is fixed within the cranium, cranial tracking allows for anatomy based tracking relatively independent of patient's daily setup. Position verification was validated several times per minute during treatment using paired, orthogonal, x-ray images.
Chi-square test or two-sample t-test was used to test the distributions of the characteristics between the simple and complex groups. To assess the association between radiation treatment parameters and the tumor volume, simple linear regressions on tumor volume for each of the three indices were performed. The estimates of the slopes and their 95% confidence intervals were determined. Pearson's correlation coefficients and their 95% confidence intervals were calculated for the whole cohort.
Patient and tumor characteristics
The characteristics of the two treatment groups including their gender, age, tumor histology and locations are detailed below and summarized in Tables 1 and 2. The simple group was composed predominantly of malignant lesions and vestibular schwannomas, while the complex group consisted primarily of cavernous sinus meningiomas and pituitary adenomas.
Overall radiosurgical parameters: effect of tumor shape
Overall, compared to previously reported conformity indices for LINAC and GammaKnife systems, the CyberKnife® radiosurgical system compared favorably with a mean NCI of 1.6–1.8 and a mean HI of 1.2–1.3 (Table 3). The standard percentage target coverage of 95% was not compromised to obtain these values.
Base of skull lesions commonly have irregular, non-spherical shapes due to the presence of dural tails and the anatomy of the region. To determine the effect of tumor shape on radiosurgical parameters, a group of spherical cerebellar and cerebral hemisphere metastases were analyzed for comparison (Control Group II (metastases)). The calculated indices for this group were similar to the indices obtained for the base of skull lesions: mean NCI of 1.73 and a mean HI of 1.21 (Table 3). These data suggest that the CyberKnife® radiosurgical system generates conformal and homogeneous plans independent of tumor shape.
Comparison of radiosurgical parameters between complex and simple base of skull lesions
Complex base of skull lesions were defined as one that completely encircles, partially circumscribes, or directly contacts the brainstem, optic chiasm, hypophysis, or cranial nerves with meaningful remaining function (see Figure 1 for example). All other lesions were classified as simple base of skull lesions (see Figure 2 for example). Table 4 gives the distribution of tumor volume, homogeneity index, new conformity index, and percentage target coverage for the simple and complex groups, respectively.
Overall, there is no statistically significant difference in homogeneity index, new conformity index and percentage target coverage between the two groups at the 5% level. There was a trend towards lower percent target coverage in the complex group, however this was not statistically significant. These data suggest that the CyberKnife® radiosurgical system generates acceptable plans independent of the proximity of adjacent critical structures to the target volume.
Relationship between tumor volume and radiosurgical parameters
Previous radiosurgical series have shown that radiosurgical indices can be influenced by target volume . In our study, the mean tumor volumes differed significantly between the simple and complex groups (p = 0.0059) (Table 4). For the simple group, the mean tumor volume was 45.6 cm3. The mean tumor volume for the complex group was smaller at 12.5 cm3. Hence, we explored the relationship between target volume and radiosurgical indices using the CyberKnife® treatment planning system. To assess the association between the three radiosurgical treatment parameters (new conformity index, homogeneity index, and percentage of tumor coverage) and the target volume, scatterplots were constructed from the data obtained from all skull base tumors (Figure 3, 4, 5).
Simple linear regressions on the tumor volume for each of the three indices were then performed. The estimates of the slopes are given in Table 5. The estimated slopes for all indices are near zero. Pearson's correlation coefficients were also calculated as seen in Table 5. All Pearson correlation coefficients were less than ± 0.4 suggesting a poor correlation between the examined variables. Therefore, tumor volume does not appear to markedly effect radiosurgical parameters when using the CyberKnife® radiosurgical treatment planning system in our patient population.
The CyberKnife® radiosurgical system has several advantages over conventional radiosurgical systems. Cranial tracking, using skeletal anatomy to position the radiation beam, is as precise as frame-based approaches and eliminates the need for headframes . In phantom studies, the system's precision has been shown to compare favorably to frame-based systems . Its sub-millimeter clinical accuracy is due both to improvements in radiation delivery and target localization [21, 22]. In addition, most LINAC and Gamma Knife systems use forward planning with user-selected arcs and beams. The CyberKnife® radiosurgical system employs inverse planning algorithms based on specific constraints to critical structures. In theory, inverse planning should allow for easily obtainable, optimized plans. The appropriate measure(s) of plan optimization is still debated .
Assessment of success in radiosurgery requires time for data to mature. But treatment-planning parameters, including conformity and homgeneity, can be assessed much earlier. In this study, we demonstrate that the CyberKnife® radiosurgical system generates plans with excellent conformity and homogeneity. Theoretically, improvements in conformity should improve local control and decrease complications in the treatment of skull base lesions with adjacent critical structures. These general principles have found acceptance in the treatment of other sites with radiation therapy [23, 24].
When irradiating complex skull base tumors that abut or displace critical normal structures the dose constraints to those normal structures may cause areas of under-dosing within the target volume. Of particular concern is that the resulting low dose regions within the tumor volume will increase the rate of local failure. In two radiosurgical series, the majority of local failures were due to tumor progression just outside the prescribed isodose volume [25, 26]. At least one report in the literature has documented that increased conformity is paradoxically associated with poorer outcomes . It has been suggested that improved conformity may lead to underdosing microscopic disease, not visible with current imaging modalities. However, in the study cited above, the poorer outcomes were likely due to the fact that conformity improves with increasing size of the lesion and is not related to an intrinsic and pure relationship between conformity and outcome. As logic dictates, increased rate of local failure is predicted to be dependent on both the dose minimum and the volume of this dose. Currently, percent target coverage is used as a surrogate for quantifying these low dose areas. In this study, percent target coverage was maintained across all groups. Longer follow-up is required to judge the effectiveness of this system in terms of local tumor control.
Dose homogeneity is a second measure by which radiosurgical plans are compared. The homogeneity index (HI), the maximum dose within the target volume divided by the prescription isodose (MDPD), is a commonly used measure of dose homogeneity. The importance of dose homogeneity in radiosurgical outcomes is controversial. Inhomogeneous high central doses achieved with some radiosurgical treatment systems may provide improved local control ; however, this increased local control may come with an increased risk of neurologic complications . A homogeneity index of less than 2.0 is felt to balance the risk of local failure and neurologic injury (RTOG guidelines) . Homogeneity indices less than 2.0 are especially important in treating large tumors or tumors in close proximity to critical structures . Even though we did not place limitations on target volume or proximity of critical structures, we were able to obtain homogeneity indices less than 2.0 for every plan. Homogeneity of dose distributions for the CyberKnife® was favorable compared with devices using multiple isocenters which are typically 2.0. In the opinion of the authors, allowable target volumes and proximity to critical structures need to be determined in the context of the homogeneity index. Larger target volumes and smaller separation from critical structures may be acceptable for systems that consistently generate low homogeneity indices .
Federal Drug Administration
Dose Volume Histogram
New Conformity Index
Planning Treatment Volume
Percent Target Coverage
Magnetic Resonance Imaging
Van Havenbergh T, Carvalho G, Tatagiba M, Plets C, Samii M: Natural history of petroclival meningiomas. Neurosurgery 2003,52(1):55-62; discussion 62-4. 10.1097/00006123-200301000-00006
Condra KS, Buatti JM, Mendenhall WM, Friedman WA, Marcus RB Jr., Rhoton AL: Benign meningiomas: primary treatment selection affects survival. Int J Radiat Oncol Biol Phys 1997,39(2):427-436. 10.1016/S0360-3016(97)00317-9
Mathiesen T, Lindquist C, Kihlstrom L, Karlsson B: Recurrence of cranial base meningiomas. Neurosurgery 1996,39(1):2-7; discussion 8-9. 10.1097/00006123-199607000-00002
De Jesus O, Sekhar LN, Parikh HK, Wright DC, Wagner DP: Long-term follow-up of patients with meningiomas involving the cavernous sinus: recurrence, progression, and quality of life. Neurosurgery 1996,39(5):915-9; discussion 919-20. 10.1097/00006123-199611000-00005
Pollock BE, Stafford SL, Utter A, Giannini C, Schreiner SA: Stereotactic radiosurgery provides equivalent tumor control to Simpson Grade 1 resection for patients with small- to medium-size meningiomas. Int J Radiat Oncol Biol Phys 2003,55(4):1000-1005. 10.1016/S0360-3016(02)04356-0
Stafford SL, Pollock BE, Leavitt JA, Foote RL, Brown PD, Link MJ, Gorman DA, Schomberg PJ: A study on the radiation tolerance of the optic nerves and chiasm after stereotactic radiosurgery. Int J Radiat Oncol Biol Phys 2003,55(5):1177-1181. 10.1016/S0360-3016(02)04380-8
Webb S: Conformal intensity-modulated radiotherapy (IMRT) delivered by robotic linac--testing IMRT to the limit? Phys Med Biol 1999,44(7):1639-1654. 10.1088/0031-9155/44/7/305
Yu C, Jozsef G, Apuzzo ML, Petrovich Z: Dosimetric comparison of CyberKnife with other radiosurgical modalities for an ellipsoidal target. Neurosurgery 2003,53(5):1155-62; discussion 1162-3. 10.1227/01.NEU.0000088805.01793.5A
Lomax NJ, Scheib SG: Quantifying the degree of conformity in radiosurgery treatment planning. Int J Radiat Oncol Biol Phys 2003,55(5):1409-1419. 10.1016/S0360-3016(02)04599-6
Theodorou K, Platoni K, Lefkopoulos D, Kappas C, Schlienger M, Dahl O: Dose-volume analysis of different stereotactic radiotherapy mono-isocentric techniques. Acta Oncol 2000,39(2):157-163. 10.1080/028418600430716
Wagner TH, Yi T, Meeks SL, Bova FJ, Brechner BL, Chen Y, Buatti JM, Friedman WA, Foote KD, Bouchet LG: A geometrically based method for automated radiosurgery planning. Int J Radiat Oncol Biol Phys 2000,48(5):1599-1611. 10.1016/S0360-3016(00)00790-2
Hamamoto Y, Manabe T, Nishizaki O, Takahashi T, Isshiki N, Murayama S, Nishina K, Umeda M: Influence of collimator size on three-dimensional conformal radiotherapy of the cyberknife. Radiat Med 2004,22(6):442-448.
Paddick I: A simple scoring ratio to index the conformity of radiosurgical treatment plans. Technical note. J Neurosurg 2000, 93 Suppl 3: 219-222.
Nedzi LA, Kooy HM, Alexander E 3rd, Svensson GK, Loeffler JS: Dynamic field shaping for stereotactic radiosurgery: a modeling study. Int J Radiat Oncol Biol Phys 1993,25(5):859-869.
Kubo HD, Wilder RB, Pappas CT: Impact of collimator leaf width on stereotactic radiosurgery and 3D conformal radiotherapy treatment plans. Int J Radiat Oncol Biol Phys 1999,44(4):937-945. 10.1016/S0360-3016(99)00041-3
Nakamura JL, Verhey LJ, Smith V, Petti PL, Lamborn KR, Larson DA, Wara WM, McDermott MW, Sneed PK: Dose conformity of gamma knife radiosurgery and risk factors for complications. Int J Radiat Oncol Biol Phys 2001,51(5):1313-1319. 10.1016/S0360-3016(01)01757-6
Shaw E, Scott C, Souhami L, Dinapoli R, Kline R, Loeffler J, Farnan N: Single dose radiosurgical treatment of recurrent previously irradiated primary brain tumors and brain metastases: final report of RTOG protocol 90-05. Int J Radiat Oncol Biol Phys 2000,47(2):291-298. 10.1016/S0360-3016(99)00507-6
Andrews DW, Scott CB, Sperduto PW, Flanders AE, Gaspar LE, Schell MC, Werner-Wasik M, Demas W, Ryu J, Bahary JP, Souhami L, Rotman M, Mehta MP, Curran WJ Jr.: Whole brain radiation therapy with or without stereotactic radiosurgery boost for patients with one to three brain metastases: phase III results of the RTOG 9508 randomised trial. Lancet 2004,363(9422):1665-1672. 10.1016/S0140-6736(04)16250-8
Adler JR Jr., Murphy MJ, Chang SD, Hancock SL: Image-guided robotic radiosurgery. Neurosurgery 1999,44(6):1299-306; discussion 1306-7. 10.1097/00006123-199906000-00079
Yu C, Main W, Taylor D, Kuduvalli G, Apuzzo ML, Adler JR Jr.: An anthropomorphic phantom study of the accuracy of Cyberknife spinal radiosurgery. Neurosurgery 2004,55(5):1138-1149. 10.1227/01.NEU.0000141080.54647.11
Chang SD, Gibbs IC, Sakamoto GT, Lee E, Oyelese A, Adler JR Jr.: Staged stereotactic irradiation for acoustic neuroma. Neurosurgery 2005,56(6):1254-61; discussion 1261-3. 10.1227/01.NEU.0000159650.79833.2B
Chang SD, Main W, Martin DP, Gibbs IC, Heilbrun MP: An analysis of the accuracy of the CyberKnife: a robotic frameless stereotactic radiosurgical system. Neurosurgery 2003,52(1):140-6; discussion 146-7. 10.1097/00006123-200301000-00018
Leibel SA, Fuks Z, Zelefsky MJ, Hunt M, Burman CM, Mageras GS, Chui CS, Jackson A, Amols HI, Ling CC: Technological advances in external-beam radiation therapy for the treatment of localized prostate cancer. Semin Oncol 2003,30(5):596-615. 10.1016/S0093-7754(03)00354-3
Ozyigit G, Yang T, Chao KS: Intensity-modulated radiation therapy for head and neck cancer. Curr Treat Options Oncol 2004,5(1):3-9.
Hakim R, Alexander E 3rd, Loeffler JS, Shrieve DC, Wen P, Fallon MP, Stieg PE, Black PM: Results of linear accelerator-based radiosurgery for intracranial meningiomas. Neurosurgery 1998,42(3):446-53; discussion 453-4. 10.1097/00006123-199803000-00002
Stafford SL, Pollock BE, Foote RL, Link MJ, Gorman DA, Schomberg PJ, Leavitt JA: Meningioma radiosurgery: tumor control, outcomes, and complications among 190 consecutive patients. Neurosurgery 2001,49(5):1029-37; discussion 1037-8. 10.1097/00006123-200111000-00001
DiBiase SJ, Kwok Y, Yovino S, Arena C, Naqvi S, Temple R, Regine WF, Amin P, Guo C, Chin LS: Factors predicting local tumor control after gamma knife stereotactic radiosurgery for benign intracranial meningiomas. Int J Radiat Oncol Biol Phys 2004,60(5):1515-1519. 10.1016/j.ijrobp.2004.05.073
Shaw E, Kline R, Gillin M, Souhami L, Hirschfeld A, Dinapoli R, Martin L: Radiation Therapy Oncology Group: radiosurgery quality assurance guidelines. Int J Radiat Oncol Biol Phys 1993,27(5):1231-1239.
Shaw E, Scott C, Souhami L, Dinapoli R, Bahary JP, Kline R, Wharam M, Schultz C, Davey P, Loeffler J, Del Rowe J, Marks L, Fisher B, Shin K: Radiosurgery for the treatment of previously irradiated recurrent primary brain tumors and brain metastases: initial report of radiation therapy oncology group protocol (90-05). Int J Radiat Oncol Biol Phys 1996,34(3):647-654. 10.1016/0360-3016(95)02106-X
The author(s) declare that they have no competing interests.
SC: Drafted the manuscript and participated in data analysis, prepared the manuscript for submission, created tables and results section
NC: Drafted the manuscript and participated in data analysis, prepared the manuscript for submission, created tables and results section
YZ: Biostatistical analysis
BC: Participated in treatment planning and manuscript revision
DM: Extracted data from treatment planning systems; manuscript revision
WJ: Participated in treatment planning and manuscript revision; corresponding author
Sean P Collins, Nicholas D Coppa contributed equally to this work.
Authors’ original submitted files for images
Below are the links to the authors’ original submitted files for images.