Imaging of regional ventilation: Is CT ventilation imaging the answer? A systematic review of the validation data

Integrating local and distant radiation-induced lung injury: Development and validation of a predictive model for ventilation loss

BACKGROUND

Investigations on radiation-induced lung injury (RILI) have predominantly focused on local effects, primarily those associated with radiation damage to lung parenchyma. However, recent studies from our group and others have revealed that radiation-induced damage to branching serial structures such as airways and vessels may also have a substantial impact on post-radiotherapy (RT) lung function. Furthermore, recent results from multiple functional lung avoidance RT trials, although promising, have demonstrated only modest toxicity reduction, likely because they were primarily focused on dose avoidance to lung parenchyma. These observations emphasize the critical need for predictive dose-response models that effectively incorporate both local and distant RILI effects.

PURPOSE

We develop and validate a predictive model for ventilation loss after lung RT. This model, referred to as P+A, integrates local (parenchyma [P]) and distant (central and peripheral airways [A]) radiation-induced damage, modeling partial (narrowing) and complete (collapse) obstruction of airways.

METHODS

In an IRB-approved prospective study, pre-RT breath-hold CTs (BHCTs) and pre- and one-year post-RT 4DCTs were acquired from lung cancer patients treated with definitive RT. Up to 13 generations of airways were automatically segmented on the BHCTs using a research virtual bronchoscopy software. Ventilation maps derived from the 4DCT scans were utilized to quantify pre- and post-RT ventilation, serving, respectively, as input data and reference standard (RS) in model validation. To predict ventilation loss solely due to parenchymal damage (referred to as P model), we used a normal tissue complication probability (NTCP) model. Our model used this NTCP-based estimate and predicted additional loss due radiation-induced partial or complete occlusion of individual airways, applying fluid dynamics principles and a refined version of our previously developed airway radiosensitivity model. Predictions of post-RT ventilation were estimated in the sublobar volumes (SLVs) connected to the terminal airways. To validate the model, we conducted a k-fold cross-validation. Model parameters were optimized as the values that provided the lowest root mean square error (RMSE) between predicted post-RT ventilation and the RS for all SLVs in the training data. The performance of the P+A and the P models was evaluated by comparing their respective post-RT ventilation values with the RS predictions. Additional evaluation using various receiver operating characteristic (ROC) metrics was also performed.

RESULTS

We extracted a dataset of 560 SLVs from four enrolled patients. Our results demonstrated that the P+A model consistently outperformed the P model, exhibiting RMSEs that were nearly half as low across all patients (13 ± 3 percentile for the P+A model vs. 24 ± 3 percentile for the P model on average). Notably, the P+A model aligned closely with the RS in ventilation loss distributions per lobe, particularly in regions exposed to doses ≥13.5 Gy. The ROC analysis further supported the superior performance of the P+A model compared to the P model in sensitivity (0.98 vs. 0.07), accuracy (0.87 vs. 0.25), and balanced predictions.

CONCLUSIONS

These early findings indicate that airway damage is a crucial factor in RILI that should be included in dose-response modeling to enhance predictions of post-RT lung function.

Incorporation of Functional Lung Imaging Into Radiation Therapy Planning in Patients With Lung Cancer: A Systematic Review and Meta-Analysis

Our purpose was to provide an understanding of current functional lung imaging (FLI) techniques and their potential to improve dosimetry and outcomes for patients with lung cancer receiving radiation therapy (RT). Excerpta Medica dataBASE (EMBASE), PubMed, and Cochrane Library were searched from 1990 until April 2023. Articles were included if they reported on FLI in one of: techniques, incorporation into RT planning for lung cancer, or quantification of RT-related outcomes for patients with lung cancer. Studies involving all RT modalities, including stereotactic body RT and particle therapy, were included. Meta-analyses were conducted to investigate differences in dose-function parameters between anatomic and functional RT planning techniques, as well as to investigate correlations of dose-function parameters with grade 2+ radiation pneumonitis (RP). One hundred seventy-eight studies were included in the narrative synthesis. We report on FLI modalities, dose-response quantification, functional lung (FL) definitions, FL avoidance techniques, and correlations between FL irradiation and toxicity. Meta-analysis results show that FL avoidance planning gives statistically significant absolute reductions of 3.22% to the fraction of well-ventilated lung receiving 20 Gy or more, 3.52% to the fraction of well-perfused lung receiving 20 Gy or more, 1.3 Gy to the mean dose to the well-ventilated lung, and 2.41 Gy to the mean dose to the well-perfused lung. Increases in the threshold value for defining FL are associated with decreases in functional parameters. For intensity modulated RT and volumetric modulated arc therapy, avoidance planning results in a 13% rate of grade 2+ RP, which is reduced compared with results from conventional planning cohorts. A trend of increased predictive ability for grade 2+ RP was seen in models using FL information but was not statistically significant. FLI shows promise as a method to spare FL during thoracic RT, but interventional trials related to FL avoidance planning are sparse. Such trials are critical to understanding the effect of FL avoidance planning on toxicity reduction and patient outcomes.

Systematic Evaluation of the Impact of Lung Segmentation Methods on 4-Dimensional Computed Tomography Ventilation Imaging Using a Large Patient Database

PURPOSE

A novel form of lung functional imaging applied for functional avoidance radiation therapy has been developed that uses 4-dimensional computed tomography (4DCT) data and image processing techniques to calculate lung ventilation (4DCT-ventilation). Lung segmentation is a common step to define a region of interest for 4DCT-ventilation generation. The purpose of this study was to quantitatively evaluate the sensitivity of 4DCT-ventilation imaging using different lung segmentation methods.

METHODS AND MATERIALS

The 4DCT data of 350 patients from 2 institutions were used. Lung contours were generated using 3 methods: (1) reference segmentations that removed airways and pulmonary vasculature manually (Lung-Manual), (2) standard lung contours used for planning (Lung-RadOnc), and (3) artificial intelligence (AI)-based contours that removed the airways and pulmonary vasculature (Lung-AI). The AI model was based on a residual 3-dimensional U-Net and was trained using the Lung-Manual contours of 279 patients. We compared the Lung-RadOnc or Lung-AI with Lung-Manual contours for the entire 4DCT-ventilation functional avoidance process including lung segmentation (surface Dice similarity coefficient [Surface DSC]), 4DCT-ventilation generation (correlation), and subanalysis of 10 patients on a dosimetric endpoint (percentage of high functional volume of lung receiving ≥20 Gy [fV20{%}]).

RESULTS

Surface DSC comparing Lung-Manual/Lung-RadOnc and Lung-Manual/Lung-AI contours was 0.40 ± 0.06 and 0.86 ± 0.04, respectively. The correlation between 4DCT-ventilation images generated with Lung-Manual/Lung-RadOnc and Lung-Manual/Lung-AI were 0.48 ± 0.14 and 0.85 ± 0.14, respectively. The difference in fV20[%] between 4DCT-ventilation generated with Lung-Manual/Lung-RadOnc and Lung-Manual/Lung-AI was 2.5% ± 4.1% and 0.3% ± 0.5%, respectively.

CONCLUSIONS

Our work showed that using standard planning lung contours can result in significantly variable 4DCT-ventilation images. The study demonstrated that AI-based segmentations generate lung contours and 4DCT-ventilation images that are similar to those generated using manual methods. The significance of the study is that it characterizes the lung segmentation sensitivity of the 4DCT-ventilation process and develops methods that can facilitate the integration of this novel imaging in busy clinics.

Respiratory Invariant Textures From Static Computed Tomography Scans for Explainable Lung Function Characterization

PURPOSE

The inherent characteristics of lung tissue independent of breathing maneuvers may provide fundamental information for function assessment. This paper attempted to correlate textural signatures from computed tomography (CT) with pulmonary function measurements.

MATERIALS AND METHODS

Twenty-one lung cancer patients with thoracic 4-dimensional CT, DTPA-single-photon emission CT ventilation ( VNM ) scans, and available spirometry measurements (forced expiratory volume in 1 s, FEV 1 ; forced vital capacity, FVC; and FEV 1 /FVC) were collected. In subregional feature discovery, function-correlated candidates were identified from 79 radiomic features based on the statistical strength to differentiate defected/nondefected lung regions. Feature maps (FMs) of selected candidates were generated on 4-dimensional CT phases for a voxel-wise feature distribution study. Quantitative metrics were applied for validations, including the Spearman correlation coefficient (SCC) and the Dice similarity coefficient for FM- VNM spatial agreement assessments, intraclass correlation coefficient for FM interphase robustness evaluations, and FM-spirometry comparisons.

RESULTS

At the subregion level, 8 function-correlated features were identified (effect size>0.330). The FMs of candidates yielded moderate-to-strong voxel-wise correlations with the reference VNM . The FMs of gray level dependence matrix dependence nonuniformity showed the highest robust (intraclass correlation coefficient=0.96 and P <0.0001) spatial correlation, with median SCCs ranging from 0.54 to 0.59 throughout the 10 breathing phases. Its phase-averaged FM achieved a median SCC of 0.60, a median Dice similarity coefficient of 0.60 (0.65) for high (low) functional lung volumes, and a correlation of 0.565 (0.646) between the spatially averaged feature values and FEV 1 (FEV 1 /FVC).

CONCLUSIONS

The results provide further insight into the underlying association of specific pulmonary textures with both local ( VNM ) and global (FEV 1 /FVC, FEV 1 ) functions. Further validations of the FM generalizability and the standardization of implementation protocols are warranted before clinically relevant investigations.

Mid-treatment adaptive planning during thoracic radiation using 68 Ventilation-Perfusion Positron emission tomography

Four-Dimensional Gallium 68 Ventilation-Perfusion Positron Emission Tomography (⁶⁸Ga-4D-V/Q PET/CT) allows for dynamic imaging of lung function. To date there has been no assessment of the feasibility of adapting radiation therapy plans to changes in lung function imaged at mid-treatment function using ⁶⁸Ga-4D-V/Q PET/CT. This study assessed the potential reductions of dose to the functional lung when radiation therapy plans were adapted to avoid functional lung at the mid-treatment timepoint using volumetric arc radiotherapy (VMAT).

Methods

A prospective clinical trial (U1111-1138-4421) was performed in patients undergoing conventionally fractionated radiation therapy for non-small cell lung cancer (NSCLC). A ⁶⁸Ga-4D-V/Q PET/CT was acquired at baseline and in the 4th week of treatment. Functional lung target volumes using the ventilated and perfused lung were created. Baseline functional volumes were compared to the week 4 V/Q functional volumes to describe the change in function over time. For each patient, 3 VMAT plans were created and optimised to spare ventilated, perfused or anatomical lung. All key dosimetry metrics were then compared including dose to target volumes, dose to organs at risk and dose to the anatomical and functional sub-units of lung.

Results

25 patients had both baseline and 4 week mid treatment ⁶⁸Ga-4D-V/Q PET/CT imaging. This resulted in a total of 75 adapted VMAT plans. The HPLung volume decreased in 16/25 patients with a mean of the change in volume (cc) -28 ± 515 cc [±SD, range -996 cc to 1496 cc]. The HVLung volume increased in 13/25 patients with mean of the change in volume (cc) + 112 ± 590 cc. [±SD, range -1424 cc to 950 cc]. The functional lung sparing technique was found to be feasible with no significant differences in dose to anatomically defined organs at risk. Most patients did derive a benefit with a reduction in functional volume receiving 20 Gy (fV20) and/or functional mean lung dose (fMLD) in either perfusion and/or ventilation. Patients with the most reduction in fV20 and fMLD were those with stage III NSCLC.

Conclusion

Functional lung volumes change during treatment. Some patients benefit from using ⁶⁸Ga-4D-V/Q PET/CT in the 4th week of radiation therapy to adapt radiation plans. In these patients, the role of mid-treatment adaptation requires further prospective investigation.

Four-Dimensional Computed Tomography Ventilation Image-Guided Lung Functional Avoidance Radiation Therapy: A Single-Arm Prospective Pilot Clinical Trial

PURPOSE

The primary objective of this prospective pilot trial was to assess the safety and feasibility of lung functional avoidance radiation therapy (RT) with 4-dimensional (4D) computed tomography (CT) ventilation imaging.

METHODS AND MATERIALS

Patients with primary lung cancer or metastatic disease to the lungs to receive conventionally fractionated RT (CFRT) or stereotactic body RT (SBRT) were eligible. Standard-of-care 4D-CT scans were used to generate ventilation images through image processing/analysis. Each patient required a standard intensity modulated RT plan and ventilation image guided functional avoidance plan. The primary endpoint was the safety of functional avoidance RT, defined as the rate of grade ≥3 adverse events (AEs) that occurred ≤12 months after treatment. Protocol treatment was considered safe if the rates of grade ≥3 pneumonitis and esophagitis were <13% and <21%, respectively for CFRT, and if the rate of any grade ≥3 AEs was <28% for SBRT. Feasibility of functional avoidance RT was assessed by comparison of dose metrics between the 2 plans using the Wilcoxon signed-rank test.

RESULTS

Between May 2015 and November 2019, 34 patients with non-small cell lung cancer were enrolled, and 33 patients were evaluable (n = 24 for CFRT; n = 9 for SBRT). Median follow-up was 14.7 months. For CFRT, the rates of grade ≥3 pneumonitis and esophagitis were 4.2% (95% confidence interval, 0.1%-21.1%) and 12.5% (2.7%-32.4%). For SBRT, no patients developed grade ≥3 AEs. Compared with the standard plans, the functional avoidance plans significantly (P < .01) reduced the lung dose-function metrics without compromising target coverage or adherence to standard organs at risk constraints.

CONCLUSIONS

This study, representing one of the first prospective investigations on lung functional avoidance RT, demonstrated that the 4D-CT ventilation image guided functional avoidance RT that significantly reduced dose to ventilated lung regions could be safely administered, adding to the growing body of evidence for its clinical utility.

Assessment of Regional Lung Ventilation with Positron Emission Tomography Using the Radiofluorinated Gas [18F]SF6: Application to an Animal Model of Impaired Ventilation

PURPOSE

Clinical ventilation studies are primarily performed with computerized tomography (CT) and more often with single-photon emission Computerized tomography (SPECT) using radiolabelled aerosols, both presenting certain limitations. Here, we investigate the use of the radiofluorinated gas [¹⁸F]SF₆ as a positron emission tomography (PET) ventilation marker in an animal model of impaired lung ventilation.

PROCEDURES

Sprague-Dawley rats (n = 15) were randomly assigned to spontaneous ventilation (sham group), endotracheal administration of phosphate-buffered saline (PBS group), or endotracheal administration of lipopolysaccharide (LPS group). PET-[¹⁸F]SF₆ images (10-min acquisition) were acquired at t = 48 h after LPS or PBS administration under mechanical ventilation. CT images were acquired after each PET session. Volumes of interest were manually delineated in the lungs on CT images, and voxel-by-voxel analysis was carried out on PET images to obtain the corresponding histograms. After the imaging sessions, lungs were harvested to conduct histological analysis.

RESULTS

Ventilation studies in sham animals showed uniform distribution of [¹⁸F]SF₆ and fast elimination of the radioactivity after discontinuation of the administration. For PBS- and LPS-treated rats, ventilation defects were observed on PET images in some animals, identified as regions with low presence of the radiolabelled gas. Hypoventilated areas co-localized with regions with higher x-ray attenuation than healthy lungs on the CT images, suggesting the presence of oedema and, in some cases, atelectasis. Histograms obtained from PET images showed quasi-Gaussian distributions for control animals, while PBS- and LPS-treated animals demonstrated the presence of hypoventilated voxels. Deviation of the histograms from Gaussian distribution correlated with histological score was obtained by ex vivo histological analysis.

CONCLUSIONS

[¹⁸F]SF₆ is an appropriate marker of regional lung ventilation and may find application in the early diagnose of acute lung disease.

CT ventilation image-guided helical Tomotherapy at sparing functional lungs for locally advanced lung cancer: analysis of dose-function metrics and the impact on pulmonary toxicity

PURPOSE

CT ventilation image (CTVI)-guided radiotherapy that selectively avoids irradiating highly-functional lung regions has potential to reduce pulmonary toxicity. Considering Helical TomoTherapy (HT) has higher modulation capabilities, we investigated the capability and characteristic of HT at sparing functional lungs for locally advanced lung cancer.

METHODS AND MATERIALS

Pretreatment 4DCT scans were carried out for 17 patients. Local lung volume expansion (or contraction) during inspiration is related to the volume change at a given lung voxel and is used as a surrogate for ventilation. The ventilation maps were generated from two sets of CT images (peak-exhale and peak-inhale) by deformable registration and a Jacobian-based algorithm. Each ventilation map was normalized to percentile images. Six plans were designed for each patient: one anatomical plan without ventilation map and five functional plans incorporating ventilation map which designed to spare varying degrees of high-functional lungs that were defined as the top 10%, 20%, 30%, 40%, and 50% of the percentile ventilation ranges, respectively. The dosimetric and evaluation factors were recorded regarding planning target volume (PTV) and other organs at risk (OARs), with particular attention to the dose delivered to total lung and functional lungs. An established dose-function-based normal tissue complication probability (NTCP) model was used to estimate risk of radiation pneumonitis (RP) for each scenario.

RESULTS

Patients were divided into a benefit group (8 patients) and a non-benefit group (9 patients) based on whether the RP-risk of functional plan was lower than that of anatomical plan. The distance between high-ventilated region and PTV, as well as tumor volume had significant differences between the two groups (P < 0.05). For patients in the benefit group, the mean value of fV5, fV10, fV20, and fMLD (functional V5, V10, V20, and mean lung dose, respectively) were significantly lower starting from top 30% functional plan than in anatomical plan (P < 0.05). With expand of avoidance region in functional plans, the dose coverage of PTV is not sacrificed (P > 0.05) but at the cost of increased dose received by OARs.

CONCLUSION

Ventilation image-guided HT plans can reduce the dose received by highly-functional lung regions with a range up to top 50% ventilated area. The spatial distribution of ventilation and tumor size were critical factors to better select patients who could benefit from the functional plan.

Radiation pneumonitis prediction model with integrating multiple dose-function features on 4DCT ventilation images

PURPOSE

Radiation pneumonitis (RP) is dose-limiting toxicity for non-small-cell cancer (NSCLC). This study developed an RP prediction model by integrating dose-function features from computed four-dimensional computed tomography (4DCT) ventilation using the least absolute shrinkage and selection operator (LASSO).

METHODS

Between 2013 and 2020, 126 NSCLC patients were included in this study who underwent a 4DCT scan to calculate ventilation images. We computed two sets of candidate dose-function features from (1) the percentage volume receiving > 20 Gy or the mean dose on the functioning zones determined with the lower cutoff percentile ventilation value, (2) the functioning zones determined with lower and upper cutoff percentile ventilation value using 4DCT ventilation images. An RP prediction model was developed by LASSO while simultaneously determining the regression coefficient and feature selection through fivefold cross-validation.

RESULTS

We found 39.3 % of our patients had a ≥ grade 2 RP. The mean area under the curve (AUC) values for the developed models using clinical, dose-volume, and dose-function features with a lower cutoff were 0.791, and the mean AUC values with lower and upper cutoffs were 0.814. The relative regression coefficient (RRC) on dose-function features with upper and lower cutoffs revealed a relative impact of dose to each functioning zone to RP. RRCs were 0.52 for the mean dose on the functioning zone, with top 20 % of all functioning zone was two times greater than that of 0.19 for these with 60 %-80 % and 0.17 with 40 %-60 % (P < 0.01).

CONCLUSIONS

The introduction of dose-function features computed from functioning zones with lower and upper cutoffs in a machine learning framework can improve RP prediction. The RRC given by LASSO using dose-function features allows for the quantification of the RP impact of dose on each functioning zones and having the potential to support treatment planning on functional image-guided radiotherapy.

Quantification of lung function on CT images based on pulmonary radiomic filtering

PURPOSE

To develop a radiomics filtering technique for characterizing spatial-encoded regional pulmonary ventilation information on lung computed tomography (CT).

METHODS

The lung volume was segmented on 46 CT images, and a 3D sliding window kernel was implemented across the lung volume to capture the spatial-encoded image information. Fifty-three radiomic features were extracted within the kernel, resulting in a fourth-order tensor object. As such, each voxel coordinate of the original lung was represented as a 53-dimensional feature vector, such that radiomic features could be viewed as feature maps within the lungs. To test the technique as a potential pulmonary ventilation biomarker, the radiomic feature maps were compared to paired functional images (Galligas PET or DTPA-SPECT) based on the Spearman correlation (ρ) analysis.

RESULTS

The radiomic feature maps GLRLM-based Run-Length Non-Uniformity and GLCOM-based Sum Average are found to be highly correlated with the functional imaging. The achieved ρ (median [range]) for the two features are 0.46 [0.05, 0.67] and 0.45 [0.21, 0.65] across 46 patients and 2 functional imaging modalities, respectively.

CONCLUSIONS

The results provide evidence that local regions of sparsely encoded heterogeneous lung parenchyma on CT are associated with diminished radiotracer uptake and measured lung ventilation defects on PET/SPECT imaging. These findings demonstrate the potential of radiomics to serve as a complementary tool to the current lung quantification techniques and provide hypothesis-generating data for future studies.

Computational lung modelling in respiratory medicine

Computational modelling of the lungs is an active field of study that integrates computational advances with lung biophysics, biomechanics, physiology and medical imaging to promote individualized diagnosis, prognosis and therapy evaluation in lung diseases. The complex and hierarchical architecture of the lung offers a rich, but also challenging, research area demanding a cross-scale understanding of lung mechanics and advanced computational tools to effectively model lung biomechanics in both health and disease. Various approaches have been proposed to study different aspects of respiration, ranging from compartmental to discrete micromechanical and continuum representations of the lungs. This article reviews several developments in computational lung modelling and how they are integrated with preclinical and clinical data. We begin with a description of lung anatomy and how different tissue components across multiple length scales affect lung mechanics at the organ level. We then review common physiological and imaging data acquisition methods used to inform modelling efforts. Building on these reviews, we next present a selection of model-based paradigms that integrate data acquisitions with modelling to understand, simulate and predict lung dynamics in health and disease. Finally, we highlight possible future directions where computational modelling can improve our understanding of the structure–function relationship in the lung.

Investigating the use of machine learning to generate ventilation images from CT scans

Background

Radiotherapy treatment planning incorporating ventilation imaging can reduce the incidence of radiation‐induced lung injury. The gold‐standard of ventilation imaging, using nuclear medicine, has limitations with respect to availability and cost.

Purpose

An alternative type of ventilation imaging to nuclear medicine uses 4DCT (or breath‐hold CT [BHCT] pair) with deformable image registration (DIR) and a ventilation metric to produce a CT ventilation image (CTVI). The purpose of this study is to investigate the application of machine learning as an alternative to DIR‐based methods when producing CTVIs.

Methods

A patient dataset of 15 inhale and exhale BHCTs and Galligas PET ventilation images were used to train and test a 2D U‐Net style convolutional neural network. The neural network established relationships between axial input BHCT image pairs and axial labeled Galligas PET images and was evaluated using eightfold cross‐validation. Once trained, the neural network could produce a CTVI from an input BHCT image pair. The CTVIs produced by the neural network were qualitatively assessed visually and quantitatively compared to a Galligas PET ventilation image using a Spearman correlation and Dice similarity coefficient (DSC). The DSC measured the spatial overlap between three segmented equal lung volumes by ventilation (high, medium, and low functioning lung [LFL]).

Results

The mean Spearman correlation between the CTVIs and the Galligas PET ventilation images was 0.58 ± 0.14. The mean DSC over high, medium, and LFL between the CTVIs and Galligas PET ventilation images was 0.55 ± 0.06. Visually, a systematic overprediction of ventilation within the lung was observed in the CTVIs with respect to the Galligas PET ventilation images, with jagged regions of ventilation in the sagittal and coronal planes.

Conclusions

A convolutional neural network was developed that could produce a CTVI from a BHCT image pair, which was then compared with a Galligas PET ventilation image. The performance of this machine learning method was comparable to previous benchmark studies investigating a DIR‐based CTVI, warranting future development, and investigation of applying machine learning to a CTVI.

Artificial intelligence in functional imaging of the lung

Artificial intelligence (AI) is transforming the way we perform advanced imaging. From high-resolution image reconstruction to predicting functional response from clinically acquired data, AI is promising to revolutionize clinical evaluation of lung performance, pushing the boundary in pulmonary functional imaging for patients suffering from respiratory conditions. In this review, we overview the current developments and expound on some of the encouraging new frontiers. We focus on the recent advances in machine learning and deep learning that enable reconstructing images, quantitating, and predicting functional responses of the lung. Finally, we shed light on the potential opportunities and challenges ahead in adopting AI for functional lung imaging in clinical settings.

Rapid Progress in Intelligent Radiotherapy and Future Implementation

Radiotherapy is one of the major approaches to cancer treatment. Artificial intelligence in radiotherapy (shortly, Intelligent radiotherapy) mainly involves big data, deep learning, extended reality, digital twin, radiomics, Internet plus and Internet of Things (IoT), which establish an automatic and intelligent network platform consisting of radiotherapy preparation, target volume delineation, treatment planning, radiation delivery, quality assurance (QA) and quality control (QC), prognosis judgment and post-treatment follow-up. Intelligent radiotherapy is an interdisciplinary frontier discipline in infancy. The review aims to summary the important implements of intelligent radiotherapy in various areas and put forward the future of unmanned radiotherapy center.

Increased regional ventilation as early imaging marker for future disease progression of interstitial lung disease: a feasibility study

OBJECTIVES

Idiopathic pulmonary fibrosis (IPF) is a disease with a poor prognosis and a highly variable course. Pathologically increased ventilation-accessible by functional CT-is discussed as a potential predecessor of lung fibrosis. The purpose of this feasibility study was to investigate whether increased regional ventilation at baseline CT and morphological changes in the follow-up CT suggestive for fibrosis indeed occur in spatial correspondence.

METHODS

In this retrospective study, CT scans were performed at two time points between September 2016 and November 2020. Baseline ventilation was divided into four categories ranging from low, normal to moderately, and severely increased (C1-C4). Correlation between baseline ventilation and volume and density change at follow-up was investigated in corresponding voxels. The significance of the difference of density and volume change per ventilation category was assessed using paired t-tests with a significance level of p ≤ 0.05. The analysis was performed separately for normal (NAA) and high attenuation areas (HAA).

RESULTS

The study group consisted of 41 patients (73 ± 10 years, 36 men). In both NAA and HAA, significant increases of density and loss of volume were seen in areas of severely increased ventilation (C4) at baseline compared to areas of normal ventilation (C2, p < 0.001). In HAA, morphological changes were more heterogeneous compared to NAA.

CONCLUSION

Functional CT assessing the extent and distribution of lung parenchyma with pathologically increased ventilation may serve as an imaging marker to prospectively identify lung parenchyma at risk for developing fibrosis.

KEY POINTS

• Voxelwise correlation of serial CT scans suggests spatial correspondence between increased ventilation at baseline and structural changes at follow-up. • Regional assessment of pathologically increased ventilation at baseline has the potential to prospectively identify tissue at risk for developing fibrosis. • Presence and extent of pathologically increased ventilation may serve as an early imaging marker of disease activity.

Synthetic pulmonary perfusion images from 4DCT for functional avoidance using deep learning

Purpose.To develop and evaluate the performance of a deep learning model to generate synthetic pulmonary perfusion images from clinical 4DCT images for patients undergoing radiotherapy for lung cancer.Methods. A clinical data set of 58 pre- and post-radiotherapy^99mTc-labeled MAA-SPECT perfusion studies (32 patients) each with contemporaneous 4DCT studies was collected. Using the inhale and exhale phases of the 4DCT, a 3D-residual network was trained to create synthetic perfusion images utilizing the MAA-SPECT as ground truth. The training process was repeated for a 50-imaging study, five-fold validation with twenty model instances trained per fold. The highest performing model instance from each fold was selected for inference upon the eight-study test set. A manual lung segmentation was used to compute correlation metrics constrained to the voxels within the lungs. From the pre-treatment test cases (N = 5), 50th percentile contours of well-perfused lung were generated from both the clinical and synthetic perfusion images and the agreement was quantified.Results. Across the hold-out test set, our deep learning model predicted perfusion with a Spearman correlation coefficient of 0.70 (IQR: 0.61-0.76) and a Pearson correlation coefficient of 0.66 (IQR: 0.49-0.73). The agreement of the functional avoidance contour pairs was Dice of 0.803 (IQR: 0.750-0.810) and average surface distance of 5.92 mm (IQR: 5.68-7.55).Conclusion. We demonstrate that from 4DCT alone, a deep learning model can generate synthetic perfusion images with potential application in functional avoidance treatment planning.

Radiation-induced Hounsfield unit change correlates with dynamic CT perfusion better than 4DCT-based ventilation measures in a novel-swine model

To analyze radiation induced changes in Hounsfield units and determine their correlation with changes in perfusion and ventilation. Additionally, to compare the post-RT changes in human subjects to those measured in a swine model used to quantify perfusion changes and validate their use as a preclinical model. A cohort of 5 Wisconsin Miniature Swine (WMS) were studied. Additionally, 19 human subjects were recruited as part of an IRB approved clinical trial studying functional avoidance radiation therapy for lung cancer and were treated with SBRT. Imaging (a contrast enhanced dynamic perfusion CT in the swine and 4DCT in the humans) was performed prior to and post-RT. Jacobian elasticity maps were calculated on all 4DCT images. Contours were created from the isodose lines to discretize analysis into 10 Gy dose bins. B-spline deformable image registration allowed for voxel-by-voxel comparative analysis in these contours between timepoints. The WMS underwent a research course of 60 Gy in 5 fractions delivered locally to a target in the lung using an MRI-LINAC system. In the WMS subjects, the dose-bin contours were copied onto the contralateral lung, which received < 5 Gy for comparison. Changes in HU and changes in Jacobian were analyzed in these contours. Statistically significant (p < 0.05) changes in the mean HU value post-RT compared to pre-RT were observed in both the human and WMS groups at all timepoints analyzed. The HU increased linearly with dose for both groups. Strong linear correlation was observed between the changes seen in the swine and humans (Pearson coefficient > 0.97, p < 0.05) at all timepoints. Changes seen in the swine closely modeled the changes seen in the humans at 12 months post RT (slope = 0.95). Jacobian analysis showed between 30 and 60% of voxels were damaged post-RT. Perfusion analysis in the swine showed a statistically significant (p < 0.05) reduction in contrast inside the vasculature 3 months post-RT compared to pre-RT. The increases in contrast outside the vasculature was strongly correlated (Pearson Correlation 0.88) with the reduction in HU inside the vasculature but were not correlated with the changes in Jacobians. Radiation induces changes in pulmonary anatomy at 3 months post-RT, with a strong linear correlation with dose. The change in HU seen in the non-vessel lung parenchyma suggests this metric is a potential biomarker for change in perfusion. Finally, this work suggests that the WMS swine model is a promising pre-clinical model for analyzing radiation-induced changes in humans and poses several benefits over conventional swine models.

A comparison of two methods for segmentation of functional volumes in radiotherapy planning of lung cancer patients

BACKGROUND

In radiotherapy (RT) of lung cancer, dose to functional lung (FL) volumes segmented with two different methods (perfusion SPECT (Q-SPECT) and 4D-CT (4D) ventilation (V)) have been shown to correlate with the incidence of radiation pneumonitis (RP). This study aims to compare the FL volumes identified by both methods.

MATERIAL AND METHODS

Thirty lung cancer patients had a 4D and Q-SPECT prior to treatment. Seventeen of these patients also had a ventilation SPECT (V-SPECT). FL sub-volumes were segmented automatically, using cut-off values. The volumes were compared in terms of overlap fraction (OF) relative to the minimal volume, and intersection fraction (IF) of the FL volume relative to the total lung volume (V_Lung).

RESULTS

Cut-off values suggested in literature for Q-SPECT and 4D-V resulted in volumes differing in size by a median 18% [6%;31%], and a median OF and IF of 0.48 [0.23;0.70] and 0.09 [0.02;0.25], respectively. Segmenting volumes of comparable size of about 1/3 of V_Lung (FL-m(1/3), m = method) resulted in a median OF and IF of 0.43 [0.23;0.58] and 0.12 [0.06;0.19], respectively. Twenty-five patients (83%) had a reasonable overlap between FL-Q(1/3) and FL-4D-V(1/3) volumes, with OF values above 0.33. IF increased significantly (p = .036) compared to using fixed cut-off values. Similarly, volumes of comparable size of about 1/3 V_Lung were produced for V-SPECT, and FL-Q(1/3), FL-V(1/3), and FL-4D-V(1/3) were compared. The overlaps and intersections of FL-V(1/3) with FL-Q(1/3) volumes were significantly (p<.001) larger than the corresponding overlaps and intersections of FL-Q(1/3) with FL-4D(1/3) and FL-V(1/3) with FL-4D(1/3).

CONCLUSION

The Q-SPECT and 4D-V methods do not segment entirely the same FL volumes. A reasonable overlap of the volumes along with the findings of other studies that both correlate to RP incidence, suggests that a combination of both volumes, e.g. using the IF, may be useful in RT treatment planning.

Pulmonary Ventilation Maps Generated with Free-breathing Proton MRI and a Deep Convolutional Neural Network

Background Hyperpolarized noble gas MRI helps measure lung ventilation, but clinical translation remains limited. Free-breathing proton MRI may help quantify lung function using existing MRI systems without contrast material and may assist in providing information about ventilation not visible to the eye or easily extracted with segmentation methods. Purpose To explore the use of deep convolutional neural networks (DCNNs) to generate synthetic MRI ventilation scans from free-breathing MRI (deep learning [DL] ventilation MRI)-derived specific ventilation maps as a surrogate of noble gas MRI and to validate this approach across a wide range of lung diseases. Materials and Methods In this secondary analysis of prospective trials, 114 paired noble gas MRI and two-dimensional free-breathing MRI scans were obtained in healthy volunteers with no history of chronic or acute respiratory disease and in study participants with a range of different obstructive lung diseases, including asthma, bronchiectasis, chronic obstructive pulmonary disease, and non-small-cell lung cancer between September 2013 and April 2018 (ClinicalTrials.gov identifiers: NCT03169673, NCT02351141, NCT02263794, NCT02282202, NCT02279329, and NCT02002052). A U-Net-based DCNN model was trained to map free-breathing proton MRI to hyperpolarized helium 3 (³He) MRI ventilation and validated using a sixfold validation. During training, the DCNN ventilation maps were compared with noble gas MRI scans using the Pearson correlation coefficient (r) and mean absolute error. DCNN ventilation images were segmented for ventilation and ventilation defects and were compared with noble gas MRI scans using the Dice similarity coefficient (DSC). Relationships were evaluated with the Spearman correlation coefficient (r_S). Results One hundred fourteen study participants (mean age, 56 years ± 15 [standard deviation]; 66 women) were evaluated. As compared with ³He MRI, DCNN model ventilation maps had a mean r value of 0.87 ± 0.08. The mean DSC for DL ventilation MRI and ³He MRI ventilation was 0.91 ± 0.07. The ventilation defect percentage for DL ventilation MRI was highly correlated with ³He MRI ventilation defect percentage (r_S = 0.83, P < .001, mean bias = -2.0% ± 5). Both DL ventilation MRI (r_S = -0.51, P < .001) and ³He MRI (r_S = -0.61, P < .001) ventilation defect percentage were correlated with the forced expiratory volume in 1 second. The DCNN model required approximately 2 hours for training and approximately 1 second to generate a ventilation map. Conclusion In participants with diverse pulmonary pathologic findings, deep convolutional neural networks generated ventilation maps from free-breathing proton MRI trained with a hyperpolarized noble-gas MRI ventilation map data set. The maps showed correlation with noble gas MRI ventilation and pulmonary function measurements. © RSNA, 2020 See also the editorial by Vogel-Claussen in this issue.

Functionally weighted airway sparing (FWAS): a functional avoidance method for preserving post-treatment ventilation in lung radiotherapy

Recent changes to the guidelines for screening and early diagnosis of lung cancer have increased the interest in preserving post-radiotherapy lung function. Current investigational approaches are based on spatially mapping functional regions and generating regional avoidance plans that preferentially spare highly ventilated/perfused lung. A potentially critical, yet overlooked, aspect of functional avoidance is radiation injury to peripheral airways, which serve as gas conduits to and from functional lung regions. Dose redistribution based solely on regional function may cause irreparable damage to the 'supply chain'. To address this deficiency, we propose the functionally weighted airway sparing (FWAS) method. FWAS (i) maps the bronchial pathways to each functional sub-lobar lung volume; (ii) assigns a weighting factor to each airway based on the relative contribution of the sub-volume to overall lung function; and (iii) creates a treatment plan that aims to preserve these functional pathways. To evaluate it, we used four cases from a retrospective cohort of SAbR patients treated for lung cancer. Each patient's airways were auto-segmented from a diagnostic-quality breath-hold CT using a research virtual bronchoscopy software. A ventilation map was generated from the planning 4DCT to map regional lung function. For each terminal airway, as resolved by the segmentation software, the total ventilation within the sub-lobar volume supported by that airway was estimated and used as a function-based weighting factor. Upstream airways were weighted based on the cumulative volumetric ventilation supported by corresponding downstream airways. Using a previously developed model for airway radiosensitivity, dose constraints were determined for each airway corresponding to a <5% probability of airway collapse. Airway dose constraints, ventilation scores, and clinical dose constraints were input to a swarm optimization-based inverse planning engine to create a 3D conformal SAbR plan (CRT). The FWAS plans were compared to the patients' prescribed CRT clinical plans and the inverse-optimized clinical plans. Depending on the size and location of the tumour, the FWAS plan showed superior preservation of ventilation due to airflow preservation through open pathways (i.e. cumulative ventilation score from the sub-lobar volumes of open pathways). Improvements ranged between 3% and 23%, when comparing to the prescribed clinical plans, and between 3% and 35%, when comparing to the inverse-optimized clinical plans. The three plans satisfied clinical requirements for PTV coverage and OAR dose constraints. These initial results suggest that by sparing pathways to high-functioning lung subregions it is possible to reduce post-SAbR loss of respiratory function.

Variations Between Dose-Ventilation and Dose-Perfusion Metrics in Radiation Therapy Planning for Lung Cancer

Purpose

Currently, several active clinical trials of functional lung avoidance radiation therapy using different imaging modalities for ventilation or perfusion are underway. Patients with lung cancer often show ventilation-perfusion mismatch, whereas the significance of dose-function metric remains unclear. The aim of the present study was to compare dose-ventilation metrics with dose-perfusion metrics for radiation therapy plan evaluation.

Methods and Materials

Pretreatment 4-dimensional computed tomography and ^99mTc-macroaggregated albumin single-photon emission computed tomography perfusion images of 60 patients with lung cancer treated with radiation therapy were analyzed. Ventilation images were created using the deformable image registration of 4-dimensional computed tomography image sets and image analysis for regional volume changes as a surrogate for ventilation. Ventilation and perfusion images were converted into percentile distribution images. Analyses included Pearson’s correlation coefficient and comparison of agreements between the following dose-ventilation and dose-perfusion metrics: functional mean lung dose and functional percent lung function receiving 5, 10, 20, 30, and 40 Gy (fV₅, fV₁₀, fV₂₀, fV₃₀, and fV₄₀, respectively).

Results

Overall, the dose-ventilation metrics were greater than the dose-perfusion metrics (ie, fV₂₀, 26.3% ± 9.9% vs 23.9% ± 9.8%). Correlations between the dose-ventilation and dose-perfusion metrics were strong (range, r = 0.94-0.97), whereas the agreements widely varied among patients, with differences as large as 6.6 Gy for functional mean lung dose and 11.1% for fV₂₀. Paired t test indicated that the dose-ventilation and dose-perfusion metrics were significantly different.

Conclusions

Strong correlations were present between the dose-ventilation and dose-perfusion metrics. However, the agreement between the dose-ventilation and dose-perfusion metrics widely varied among patients, suggesting that ventilation-based radiation therapy plan evaluation may not be comparable to that based on perfusion. Future studies should elucidate the correlation of dose-function metrics with clinical pulmonary toxicity metrics.