- Methodology
- Open access
- Published:
Method of computing direction-dependent margins for the development of consensus contouring guidelines
Radiation Oncology volume 16, Article number: 71 (2021)
Abstract
Background
Clinical target volume (CTV) contouring guidelines are frequently developed through studies in which experts contour the CTV for a representative set of cases for a given treatment site and the consensus CTVs are analyzed to generate margin recommendations. Measures of interobserver variability are used to quantify agreement between experts. In cases where an isotropic margin is not appropriate, however, there is no standard method to compute margins in specified directions that represent possible routes of tumor spread. Moreover, interobserver variability metrics are often measures of volume overlap that do not account for the dependence of disagreement on direction. To aid in the development of consensus contouring guidelines, this study demonstrates a novel method of quantifying CTV margins and interobserver variability in clinician-specified directions.
Methods
The proposed algorithm was applied to 11 cases of non-spine bone metastases to compute the consensus CTV margin in each direction of intraosseous and extraosseous disease. The median over all cases for each route of spread yielded the recommended margins. The disagreement between experts on the CTV margin was quantified by computing the median of the coefficients of variation for intraosseous and extraosseous margins.
Results
The recommended intraosseous and extraosseous margins were 7.0 mm and 8.0 mm, respectively. The median coefficient of variation quantifying the margin disagreement between experts was 0.59 and 0.48 for intraosseous and extraosseous disease.
Conclusions
The proposed algorithm permits the generation of margin recommendations in relation to adjacent anatomy and quantifies interobserver variability in specified directions. This method can be applied to future consensus CTV contouring studies.
Background
Conformal radiotherapy techniques, including stereotactic radiosurgery (SRS) and stereotactic body radiotherapy (SBRT), are increasingly being delivered for most malignancies. Consistent and accurate delineation of the clinical target volume (CTV) is essential for ensuring treatment efficacy and meaningful interpretation of dosimetric and clinical outcomes in clinical trials. The CTV is defined as an additional margin to the gross tumor volume (GTV), the demonstrable extent of the tumor, to encompass microscopic disease. Ideally, CTV margin recommendations are informed by pathological evaluation of microscopic tumor extension [1,2,3] and examination of patterns of disease recurrence. In the absence of such data, expert consensus can provide guidance. Over the last decade, a number of consensus contouring studies to formulate guidelines have emerged for various disease sites [4,5,6,7,8,9,10,11,12].
In these studies, experts are typically provided with patient images and asked to contour the CTV given a reference GTV. Measures of inter-observer variability are used to determine the agreement between experts [13]. A consensus CTV is often created from a collection of expert contours by applying the simultaneous truth and performance level estimation (STAPLE) algorithm [14]. Consensus contouring recommendations typically include CTV margins based, in part, on the STAPLE consensus CTV. In clinical contouring practice, the GTV would be expanded isotropically or anisotropically along the Cartesian axes using these margins.
However, the extent of tumor spread can vary in arbitrary directions because of anatomic barriers, anisotropic tissue structure, or preferential migration within the tissue of origin [15]. For certain tumor sites, CTV contouring guidelines that account for this complexity may be more appropriate than isotropic margin recommendations. In spine SRS, for example, the CTV is defined by vertebral sectors depending on GTV involvement instead of by a simple isotropic margin [5]. Automatic CTV delineation using vertebral sectors has already been implemented in treatment planning systems [16]. Future treatment planning software should also facilitate expansion in arbitrary directions using anisotropic margins for better conformation to the disease volume.
For making contouring recommendations by direction, a method of analyzing expert contours to generate margins in specified directions that represent potential routes of tumor spread is needed. Such analysis would also benefit from a metric for quantifying interobserver variability in specified directions to assess disagreement on the appropriate margin. Common measures of interobserver variability like the Dice coefficient and the kappa statistic [17] are based on volume overlap and therefore do not capture the dependence of disagreement on direction. Existing algorithms that could be applied for computation of anisotropic margins and contour variability rely on approximately spherical volumes [18, 19], which may not be a valid assumption, or do not make use of specified directions [20,21,22].
This study proposes a novel method to compute consensus CTV margins and margin variability in directions chosen based on clinical assessment of potential routes of spread and anatomical barriers. For generating clinical contouring recommendations, the consensus CTV margin is computed for specified routes of spread in multiple cases, then a margin recommendation is made based on a summary statistic across all cases. For demonstration, the method is applied to 11 cases of non-spine bone metastases, where margins are expected to be anisotropic, to determine separate intraosseous and extraosseous margins.
Methods and materials
Metrics
The upcoming metric definitions rely on a quantity called the directional margin \(M({\varvec{d}})\), where \({\varvec{d}}\) is any clinically relevant direction in 3D. The directional margin is computed from a GTV and a CTV in two steps. First, a set of \(N\) expansion vectors \(\{{{\varvec{v}}}_{i}{\}}_{i=1}^{N}\) are identified: these vectors are oriented parallel to \({\varvec{d}}\) and extend from the GTV surface to the CTV surface without intersecting the interior of the GTV. Second, the directional margin is defined as the median of the lengths of the expansion vectors: \(M({\varvec{d}})={\text{MEDIAN}}\left({\left\{\left|{{\varvec{v}}}_{i}\right|\right\}}_{i=1}^{N}\right)\), where \(|\cdot |\) is the Euclidean norm. This procedure is illustrated in Fig. 1. The diagram shows 2D contours for demonstration, but the algorithm can be applied to any contiguous 3D volume. Pseudocode is shown in an additional file (see Additional file 1).
The consensus CTV margin is the directional margin of the GTV to the consensus CTV. The margin deviation is the coefficient of variation of the directional margins from the GTV to each expert CTV. For \(K\) expert CTVs and a set of directional margins \(\{M({\varvec{d}}{)}_{j}{\}}_{j=1}^{K}\) in direction \({\varvec{d}}\), \({\text{CV}}(\mathbf{d})={\text{STD}}\left(\{M({\varvec{d}}{)}_{j}{\}}_{j=1}^{K}\right)/{\text{MEAN}}\left(\{M({\varvec{d}}{)}_{j}{\}}_{j=1}^{K}\right)\) is the margin deviation.
Application to clinical cases
The method was applied to 11 clinical cases, chosen specifically to represent a wide range of sites, and taken from a study on international practice patterns for the use of SBRT to treat non-spine bone metastases [23]. Each of nine radiation oncologists delineated the CTV contour, based on the provided GTV contour and the simulation CT and MRI scans. The contours were defined in the CT space (in-plane pixel size: 1.17 × 1.17 mm; slice thickness: 1–3 mm). Using all CTV contours, a consensus contour was computed using the STAPLE algorithm and approved by all radiation oncologists for use in making contouring recommendations [23].
For each case, directions of intraosseous and extraosseous extension were identified by two co-authors ([TKN] and [CLT]), both radiation oncologists. In each direction, the consensus margin was computed using the GTV and the consensus CTV, while the margin deviation was computed using the GTV and the individual participant CTVs. Computations were done using MATLAB (version R2016b, The Mathworks, Inc.). The median consensus margin and margin deviation for intraosseous and extraosseous extension were computed across all cases.
Results
Figures 2 and 3 illustrate the results. Figure 2 shows four directions of potential intraosseous extension for a metastasis in the pubic symphysis (the tumor showed no extraosseous component). The directions are termed superior-anterior (SA), inferior-posterior (IP), superior-right (SR), and posterior-right (PR). The consensus margins for the four directions were 6.4 (SA), 8.2 (IP), 6.4 (SR), and 7.6 (PR) mm. Figure 3 shows a metastasis in the iliac crest demonstrating cortical disruption. Three directions of potential intraosseous extension were identified: superior-posterior (SP), inferior-anterior (IA), and posterior-right (PR). Two directions of extraosseous extension were also selected: superior-right (SR) and anterior-right (AR). The intraosseous margins were 7.0 (SP), 7.0 (IA), and 8.0 (PR) mm, while the extraosseous margins were 5.5 (SR) and 6.5 (AR) mm. Note that computing a directional margin allowed separate estimates of intraosseous and extraosseous extension from the same case. The margin deviations for the pubic symphysis case (Fig. 2) were 0.62 (SA), 0.61 (IP), 0.55 (SR), and 0.59 (PR). For the iliac crest case (Fig. 3), the margin deviations for the intraosseous directions were 0.44 (SP), 0.40 (IA), and 0.47 (PR); for the extraosseous directions, the deviations were 0.43 (SR) and 0.47 (AR). The spread of margin deviation values over the different directions for both cases indicates that contouring disagreement can vary by direction.
Figure 4 shows the consensus margins and margin deviations over all 11 cases. The median intraosseous consensus margin was 7.0 mm (range 5–14 mm). The median extraosseous consensus margin was 8.0 mm (range 5.5–9.4 mm). The median margin deviations were 0.59 (range 0.40–0.73) for intraosseous disease and 0.48 (range 0.43–0.92) for extraosseous disease. These consensus margins can inform contouring recommendations for stereotactic body radiotherapy of non-spine bone metastases.
Discussion
The purpose of the directional margin algorithm is to quantify CTV expansions by direction for the development of contouring guidelines. The advantage of applying this algorithm over computing an isotropic margin is that the extent of tumor spread may depend on direction because of anatomical barriers or the nature of surrounding tissue. In an application to multiple cases of non-spine bone metastases the proposed method yielded separate estimates of intraosseous and extraosseous extension, which would not have been possible by computing an isotropic margin, as illustrated by the case in Fig. 3. This data can inform the development of subsequent consensus contouring recommendations. Since routes of spread will typically vary in different directions, anisotropic CTV expansions in arbitrary directions should be a feature of future treatment planning systems.
The agreement between experts on the appropriate margin for a given route of spread can be quantified using the margin deviation. The margin deviation was found to vary by direction (Figs. 2 and 3), illustrating the advantage of a direction-dependent measure of interobserver variability. The range of the margin deviation will need to be characterized in more subjects so that outliers can be identified in future studies on contouring variability that employ this technique.
One of the strengths of the proposed method is the ability to account for complex geometries, since the lengths of all expansion vectors are incorporated (Fig. 5a). The primary limitation of this method is that it does not account for the fact that only a subset of the GTV surface is relevant for margin computation when barriers to tumor spread are present in the specified direction (Fig. 5b). The method can be generalized by segmenting anatomical barriers and excluding inappropriate vectors. This strategy could be facilitated by recent advances in automated segmentation for CTV delineation (15, 24).
Future work will include the application of the presented methods to the same 11 cases to inform contouring recommendations for SBRT of non-spine bone metastases. Applying the method to other tumor sites is another direction for future studies.
Conclusions
This report described and demonstrated a novel method to compute CTV margins and margin disagreement between expert contours in any number of specified clinically relevant directions. The proposed approach allows for the generation of margin recommendations in relation to adjacent anatomy. The method was applied to 11 cases of non-spine bone metastases to compute recommended intraosseous and extraosseous margins. The margin disagreement was also quantified. The method can be applied in future consensus contouring studies.
Availability of data and materials
The datasets used and/or analyzed during the current study are available from the authors on reasonable request.
References
Teh BS, Bastasch MD, Wheeler TM, Mai W-Y, Frolov A, Uhl BM, et al. IMRT for prostate cancer: defining target volume based on correlated pathologic volume of disease. Int J Radiat Oncol Biol Phys. 2003;56(1):184–91.
Wang W, Feng X, Zhang T, Jin J, Wang S, Liu Y, et al. Prospective evaluation of microscopic extension using whole-mount preparation in patients with hepatocellular carcinoma: Definition of clinical target volume for radiotherapy. Radiat Oncol. 2010;5(73):1–7.
Jenkins P, Anjarwalla S, Gilbert H, Kinder R. Defining the clinical target volume for bladder cancer radiotherapy treatment planning. Int J Radiat Oncol Biol Phys. 2009;75(5):1379–84.
Soliman H, Ruschin M, Angelov L, Brown PD, Chiang VLS, Kirkpatrick JP, et al. Consensus contouring guidelines for postoperative completely resected cavity stereotactic radiosurgery for brain metastases. Int J Radiat Oncol Biol Phys. 2018;100(2):436–42.
Cox BW, Spratt DE, Lovelock M, Bilsky MH, Lis E, Ryu S, et al. International spine radiosurgery consortium consensus guidelines for target volume definition in spinal stereotactic radiosurgery. Int J Radiat Oncol Biol Phys. 2012;83(5):e597-605.
Redmond KJ, Robertson S, Lo SS, Soltys SG, Ryu S, McNutt T, et al. Consensus contouring guidelines for postoperative stereotactic body radiation therapy for metastatic solid tumor malignancies to the spine. Int J Radiat Oncol Biol Phys. 2017;97(1):64–74.
Dunne EM, Sahgal A, Lo SS, Bergman A, Kosztyla R, Dea N, et al. International consensus recommendations for target volume delineation specific to sacral metastases and spinal stereotactic body radiation therapy (SBRT). Radiother Oncol. 2020;145:21–9.
Lim K, Small W, Portelance L, Creutzberg C, Jürgenliemk-Schulz IM, Mundt A, et al. Consensus guidelines for delineation of clinical target volume for intensity-modulated pelvic radiotherapy for the definitive treatment of cervix cancer. Int J Radiat Oncol Biol Phys. 2011;79(2):348–55.
Myerson RJ, Garofalo MC, El Naqa I, Abrams RA, Apte A, Bosch WR, et al. Elective clinical target volumes for conformal therapy in anorectal cancer: a Radiation Therapy Oncology Group consensus panel contouring atlas. Int J Radiat Oncol Biol Phys. 2009;74(3):824–30.
Wu AJ, Bosch WR, Chang DT, Hong TS, Jabbour SK, Kleinberg LR, et al. Expert consensus contouring guidelines for intensity modulated radiation therapy in esophageal and gastroesophageal junction cancer. Int J Radiat Oncol Biol Phys. 2015;92(4):911–20.
Leung E, D’Souza D, Bachand F, Han K, Alfieri J, Huang F, et al. MRI-based interstitial brachytherapy for vaginal tumors: A multi-institutional study on practice patterns, contouring, and consensus definitions of target volumes. Brachytherapy. 2019;18(5):598–605.
Tseng C-L, Stewart J, Whitfield G, Verhoeff JJC, Bovi J, Soliman H, et al. Glioma consensus contouring recommendations from a MR-Linac International Consortium Research Group and evaluation of a CT-MRI and MRI-only workflow. J Neurooncol. 2020. https://0-doi-org.brum.beds.ac.uk/10.1007/s11060-020-03605-6.
Vinod SK, Jameson MG, Min M, Holloway LC. Uncertainties in volume delineation in radiation oncology: A systematic review and recommendations for future studies. Radiother Oncol. 2016;121:169–79.
Warfield SK, Zou KH, Wells WM. Simultaneous Truth and Performance Level Estimation (STAPLE): An algorithm for the validation of image segmentation. IEEE Trans Med Imaging. 2004;23(7):903–21.
Unkelbach J, Bortfeld T, Cardenas CE, Gregoire V, Hager W, Heijmen B, et al. The role of computational methods for automating and improving clinical target volume definition. Radiother Oncol. 2020;153:15–25.
Giaj-Levra N, Niyazi M, Figlia V, Napoli G, Mazzola R, Nicosia L, et al. Feasibility and preliminary clinical results of linac-based stereotactic body radiotherapy for spinal metastases using a dedicated contouring and planning system. Radiat Oncol. 2019;14:184.
Allozi R, Li XA, White J, Apte A, Tai A, Michalski JM, et al. Tools for consensus analysis of experts’ contours for radiotherapy structure definitions. Radiother Oncol. 2010;97(3):572–8.
Rasch C, Barillot I, Remeijer P, Touw A, van Herk M, Lebesque JV. Definition of the prostate in CT and MRI: a multi-observer study. Int J Radiat Oncol. 1999;43(1):57–66.
Song WY, Chiu B, Bauman GS, Lock M, Rodrigues G, Ash R, et al. Prostate contouring uncertainty in megavoltage computed tomography images acquired with a helical tomotherapy unit during image-guided radiation therapy. Int J Radiat Oncol Biol Phys. 2006;65(2):595–607.
Deurloo KEI, Steenbakkers RJHM, Zijp LJ, de Bois JA, Nowak PJCM, Rasch CRN, et al. Quantification of shape variation of prostate and seminal vesicles during external beam radiotherapy. Int J Radiat Oncol Biol Phys. 2005;61(1):228–38.
Jansen EPM, Nijkamp J, Gubanski M, Lind PARM, Verheij M. Interobserver variation of clinical target volume delineation in gastric cancer. Int J Radiat Oncol Biol Phys. 2010;77(4):1166–70.
Weiss E, Wu J, Sleeman W, Bryant J, Mitra P, Myers M, et al. Clinical evaluation of soft tissue organ boundary visualization on cone-beam computed tomographic imaging. Int J Radiat Oncol Biol Phys. 2010;78(3):929–36.
Nguyen TK, Sahgal A, Dagan R, Eppinga W, Guckenberger M, Kim JH, et al. Stereotactic body radiation therapy for nonspine bone metastases: international practice patterns to guide treatment planning. Pract Radiat Oncol. 2020;10(6):e452–60.
Shusharina N, Söderberg J, Edmunds D, Löfman F, Shih H, Bortfeld T. Automated delineation of the clinical target volume using anatomically constrained 3D expansion of the gross tumor volume. Radiother Oncol. 2020;146:37–43.
Acknowledgements
Not applicable.
Funding
This research received support from the Natural Sciences and Engineering Research Council (NSERC) (RGPIN-2017-06596, CRD 507521-16). LSPL was funded by an NSERC CGS D scholarship and a Walter C. Sumner Memorial Fellowship.
Author information
Authors and Affiliations
Contributions
LSPL contributed to the formal analysis, methodology, software, and visualization. LC, RWC, and AZL contributed to the methodology and supervision. AZL also contributed to funding acquisition. TKN and CLT contributed to the formal analysis and methodology. AS contributed to conceptualization, funding acquisition, and resources. All authors contributed to writing, reviewing, and editing the manuscript.
Corresponding author
Ethics declarations
Ethics approval and consent to participate
This study was approved by the Sunnybrook institutional research ethics board.
Consent for publication
Not applicable.
Competing interests
CLT has received travel accommodations/expenses and honoraria for past educational seminars from Elekta and belongs to the Elekta MR-Linac Research Consortium. AS has been an advisor/consultant with Abbvie, Merck, Roche, Varian (Medical Advisory Group), Elekta (Gamma Knife Icon), BrainLAB, and VieCure (Medical Advisor Board); ex officio Board Member to International Stereotactic Radiosurgery Society (ISRS); received honorarium for past educational seminars with Elekta AB, Accuray Inc, Varian (CNS Teaching Faculty), BrainLAB, and Medtronic Kyphon; research grant with Elekta AB; and travel accommodations/expenses by Elekta, Varian, and BrainLAB. AS also belongs to the Elekta MR-Linac Research Consortium, Elekta Spine, Oligometastases and Linac Based SRS Consortia.
Additional information
Publisher's Note
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Supplementary Information
Additional file 1.
Title: Directional margin algorithm pseudocode. Description: Pseudocode for the directional margin algorithm, which forms the basis of the proposed method.
Rights and permissions
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated in a credit line to the data.
About this article
Cite this article
Lawrence, L.S.P., Chin, L.C.L., Chan, R.W. et al. Method of computing direction-dependent margins for the development of consensus contouring guidelines. Radiat Oncol 16, 71 (2021). https://0-doi-org.brum.beds.ac.uk/10.1186/s13014-021-01799-1
Received:
Accepted:
Published:
DOI: https://0-doi-org.brum.beds.ac.uk/10.1186/s13014-021-01799-1