Synergistic combination of a topologically invariant imaging signature and a biomarker for the accurate prediction of symptomatic radiation pneumonitis before stereotactic ablative radiotherapy for lung cancer: A retrospective analysis

Objectives

We aimed to explore the synergistic combination of a topologically invariant Betti number (BN)-based signature and a biomarker for the accurate prediction of symptomatic (grade ≥2) radiation-induced pneumonitis (RP+) before stereotactic ablative radiotherapy (SABR) for lung cancer.

Methods

A total of 272 SABR cases with early-stage non-small cell lung cancer were chosen for this study. The occurrence of RP+ was predicted using a support vector machine (SVM) model trained with the combined features of the BN-based signature extracted from planning computed tomography (pCT) images and a pretreatment biomarker, serum Krebs von den Lungen-6 (BN+KL-6 model). In all, 242 (20 RP+ and 222 RP–(grade 1)) and 30 cases (8 RP+ and 22 RP–) were used for training and testing the model, respectively. The BN-based features were extracted from BN maps that characterize topologically invariant heterogeneous traits of potential RP+ lung regions on pCT images by applying histogram- and texture-based feature calculations to the maps. The SVM models were built to predict RP+ patients with a BN signature that was constructed based on the least absolute shrinkage and selection operator logistic regression model. The evaluation of the prediction models was performed based on the area under the receiver operating characteristic curves (AUCs) and accuracy in the test. The performance of the BN+KL-6 model was compared to the performance based on the BN, conventional original pCT, and wavelet decomposition (WD) models.

Results

The test AUCs obtained for the BN+KL-6, BN, pCT, and WD models were 0.825, 0.807, 0.642, and 0.545, respectively. The accuracies of the BN+KL-6, BN, pCT, and WD models were found to be 0.724, 0.708, 0.591, and 0.534, respectively.

Conclusion

This study demonstrated the comprehensive performance of the BN+KL-6 model for the prediction of potential RP+ patients before SABR for lung cancer.

Deexcitation Dynamics of Muonic Atoms Revealed by High-Precision Spectroscopy of Electronic K X Rays

We observed electronic K x rays emitted from muonic iron atoms using superconducting transition-edge sensor microcalorimeters. The energy resolution of 5.2 eV in FWHM allowed us to observe the asymmetric broad profile of the electronic characteristic Kα and Kβ x rays together with the hypersatellite K^{h}α x rays around 6 keV. This signature reflects the time-dependent screening of the nuclear charge by the negative muon and the L-shell electrons, accompanied by electron side feeding. Assisted by a simulation, these data clearly reveal the electronic K- and L-shell hole production and their temporal evolution on the 10-20 fs scale during the muon cascade process.

Robust radiogenomics approach to the identification of EGFR mutations among patients with NSCLC from three different countries using topologically invariant Betti numbers

Objectives

To propose a novel robust radiogenomics approach to the identification of epidermal growth factor receptor (EGFR) mutations among patients with non-small cell lung cancer (NSCLC) using Betti numbers (BNs).

Materials and methods

Contrast enhanced computed tomography (CT) images of 194 multi-racial NSCLC patients (79 EGFR mutants and 115 wildtypes) were collected from three different countries using 5 manufacturers’ scanners with a variety of scanning parameters. Ninety-nine cases obtained from the University of Malaya Medical Centre (UMMC) in Malaysia were used for training and validation procedures. Forty-one cases collected from the Kyushu University Hospital (KUH) in Japan and fifty-four cases obtained from The Cancer Imaging Archive (TCIA) in America were used for a test procedure. Radiomic features were obtained from BN maps, which represent topologically invariant heterogeneous characteristics of lung cancer on CT images, by applying histogram- and texture-based feature computations. A BN-based signature was determined using support vector machine (SVM) models with the best combination of features that maximized a robustness index (RI) which defined a higher total area under receiver operating characteristics curves (AUCs) and lower difference of AUCs between the training and the validation. The SVM model was built using the signature and optimized in a five-fold cross validation. The BN-based model was compared to conventional original image (OI)- and wavelet-decomposition (WD)-based models with respect to the RI between the validation and the test.

Results

The BN-based model showed a higher RI of 1.51 compared with the models based on the OI (RI: 1.33) and the WD (RI: 1.29).

Conclusion

The proposed model showed higher robustness than the conventional models in the identification of EGFR mutations among NSCLC patients. The results suggested the robustness of the BN-based approach against variations in image scanner/scanning parameters.

Radiomic features based on Hessian index for prediction of prognosis in head-and-neck cancer patients

This study demonstrated the usefulness of radiomic features based on the Hessian index of differential topology for the prediction of prognosis prior to treatment in head-and-neck (HN) cancer patients. The Hessian index, which can indicate tumor heterogeneity with convex, concave, and other points (saddle points), was calculated as the number of negative eigenvalues of the Hessian matrix at each voxel on computed tomography (CT) images. Three types of signatures were constructed in a training cohort (n = 126), one type each from CT conventional features, Hessian index features, and combined features from the conventional and index feature sets. The prognostic value of the signatures were evaluated using statistically significant difference (p value, log-rank test) to compare the survival curves of low- and high-risk groups. In a test cohort (n = 68), the p values of the models built with conventional, index, combined features, and clinical variables were 2.95 [Formula: see text] 10-2, 1.85 [Formula: see text] 10-2, 3.17 [Formula: see text] 10-2, and 1.87 [Formula: see text] 10-3, respectively. When the features were integrated with clinical variables, the p values of conventional, index, and combined features were 3.53 [Formula: see text] 10-3, 1.28 [Formula: see text] 10-3, and 1.45 [Formula: see text] 10-3, respectively. This result indicates that index features could provide more prognostic information than conventional features and further increase the prognostic value of clinical variables in HN cancer patients.

Semi-automated prediction approach of target shifts using machine learning with anatomical features between planning and pretreatment CT images in prostate radiotherapy

The goal of this study was to develop a semi-automated prediction approach of target shifts using machine learning architecture (MLA) with anatomical features for prostate radiotherapy. Our hypothesis was that anatomical features between planning computed tomography (pCT) and pretreatment cone-beam computed tomography (CBCT) images could be used to predict the target, i.e. clinical target volume (CTV) shifts, with small errors. The pCT and daily CBCT images of 20 patients with prostate cancer were selected. The first 10 patients were employed for the development, and the second 10 patients for a validation test. The CTV position errors between the pCT and CBCT images were determined as reference CTV shifts (teacher data) after an automated bone-based registration. The anatomical features associated with rectum, bladder and prostate were calculated from the pCT and CBCT images. The features were fed as the input with the teacher data into five MLAs, i.e. three types of artificial neural networks, support vector regression (SVR) and random forests. Since the CTV shifts along the left-right direction were negligible, the MLAs were developed along the superior-inferior and anterior-posterior directions. The proposed framework was evaluated from the residual errors between the reference and predicted CTV shifts. In the validation test, the mean residual error with its standard deviation was 1.01 ± 1.09 mm in SVR using only one feature (one click), which was associated with positional difference of the upper rectal wall. The results suggested that MLAs with anatomical features could be useful in prediction of CTV shifts for prostate radiotherapy.

Automated approach for estimation of grade groups for prostate cancer based on histological image feature analysis

BACKGROUND

There is a low reproducibility of the Gleason scores that determine the grade group of prostate cancer given the intra- and interobserver variability among pathologists. This study aimed to develop an automated approach for estimating prostate cancer grade groups based on features obtained from histological image analysis.

METHODS

Fifty-nine patients who underwent radical prostatectomy were selected under the approval of the institutional review board of our university hospital. For estimation, we followed the grade group criteria provided by the International Society of Urological Pathology in 2014. One hundred eight specimen slides obtained from the patients were digitized to extract 110 regions of interest (ROI) from hematoxylin and eosin-stained histological images using a digital whole slide scanner at ×20 magnification with a pixel size of 0.4 μm. Each color pixel value in the ROI was decomposed into six intensities corresponding to the RGB (red, green, and blue) and HSV (hue, saturation, and value) color models. Image features were extracted by histological image analysis, obtaining 54 features from the ROI based on histogram and texture analyses in the six types of decomposed histological images. Then, 40 representative features were selected from the 324 histological image features based on statistically significant differences (P < .05) between the mean image feature values for high (≥3, Gleason score ≥4 + 3) and low (≤2, Gleason score ≤3 + 4) grade groups. The relationship between grade groups and the most representative image feature (ie, complexity) was approximated using regression to estimate real-number grade groups defined by continuous numerical grading. Finally, the grade groups were expressed as the conventional grade groups (ie, integers from 1 to 5) using a piecewise step function.

RESULTS

The grade groups were correctly estimated by the proposed approach without errors on training (70 ROIs) and validation (40 ROIs) data.

CONCLUSIONS

Our results suggest that the proposed approach may support pathologists during the evaluation of grade groups for prostate cancer, thus mitigating intra- and interobserver variability.

Homological radiomics analysis for prognostic prediction in lung cancer patients

PURPOSE

This study explored a novel homological analysis method for prognostic prediction in lung cancer patients.

MATERIALS AND METHODS

The potential of homology-based radiomic features (HFs) was investigated by comparing HFs to conventional wavelet-based radiomic features (WFs) and combined radiomic features consisting of HFs and WFs (HWFs), using training (n = 135) and validation (n = 70) datasets, and Kaplan-Meier analysis. A total of 13,824 HFs were derived through homology-based texture analysis using Betti numbers, which represent the topologically invariant morphological characteristics of lung cancer. The prognostic potential of HFs was evaluated using statistically significant differences (p-values, log-rank test) to compare the survival curves of high- and low-risk patients. Those patients were stratified into high- and low-risk groups using the medians of the radiomic scores of signatures constructed with an elastic-net-regularized Cox proportional hazard model. Furthermore, deep learning (DL) based on AlexNet was utilized to compare HFs by stratifying patients into the two groups using a network that was pre-trained with over one million natural images from an ImageNet database.

RESULTS

For the training dataset, the p-values between the two survival curves were 6.7 × 10^-6 (HF), 5.9 × 10^-3 (WF), 7.4 × 10^-6 (HWF), and 1.1 × 10^-3 (DL). The p-values for the validation dataset were 3.4 × 10^-5 (HF), 6.7 × 10^-1 (WF), 1.7 × 10^-7 (HWF), and 1.2 × 10^-1 (DL).

CONCLUSION

This study demonstrates the excellent potential of HFs for prognostic prediction in lung cancer patients.

Radiomics with artificial intelligence for precision medicine in radiation therapy

Recently, the concept of radiomics has emerged from radiation oncology. It is a novel approach for solving the issues of precision medicine and how it can be performed, based on multimodality medical images that are non-invasive, fast and low in cost. Radiomics is the comprehensive analysis of massive numbers of medical images in order to extract a large number of phenotypic features (radiomic biomarkers) reflecting cancer traits, and it explores the associations between the features and patients' prognoses in order to improve decision-making in precision medicine. Individual patients can be stratified into subtypes based on radiomic biomarkers that contain information about cancer traits that determine the patient's prognosis. Machine-learning algorithms of AI are boosting the powers of radiomics for prediction of prognoses or factors associated with treatment strategies, such as survival time, recurrence, adverse events, and subtypes. Therefore, radiomic approaches, in combination with AI, may potentially enable practical use of precision medicine in radiation therapy by predicting outcomes and toxicity for individual patients.

Potentials of radiomics for cancer diagnosis and treatment in comparison with computer-aided diagnosis

Computer-aided diagnosis (CAD) is a field that is essentially based on pattern recognition that improves the accuracy of a diagnosis made by a physician who takes into account the computer's "opinion" derived from the quantitative analysis of radiological images. Radiomics is a field based on data science that massively and comprehensively analyzes a large number of medical images to extract a large number of phenotypic features reflecting disease traits, and explores the associations between the features and patients' prognoses for precision medicine. According to the definitions for both, you may think that radiomics is not a paraphrase of CAD, but you may also think that these definitions are "image manipulation". However, there are common and different features between the two fields. This review paper elaborates on these common and different features and introduces the potential of radiomics for cancer diagnosis and treatment by comparing it with CAD.

Diel patterns of larval drift of honmoroko Gnathopogon caerulescens in an inlet of Ibanaiko Lagoon, Lake Biwa, Japan

Diel drift patterns of larvae of the endangered cyprinid Gnathopogon caerelescens in an inlet of the Ibanaiko Lagoon, connected to Lake Biwa in Japan, were assessed in April 2012. Peak occurrence of yolk-sac larvae was within a few hours after dark. Drift of newly hatched larvae is considered to be an important biological mechanism that ensures larval dispersal and recruitment from the inlets (spawning grounds) to the lagoon which functions as a nursery ground.

Measurements of the neutron activation cross sections for Bi and Co at 386 MeV

Neutron activation cross sections for Bi and Co at 386 MeV were measured by activation method. A quasi-monoenergetic neutron beam was produced using the (7)Li(p,n) reaction. The energy spectrum of these neutrons has a high-energy peak (386 MeV) and a low-energy tail. Two neutron beams, 0° and 25° from the proton beam axis, were used for sample irradiation, enabling a correction for the contribution of the low-energy neutrons. The neutron-induced activation cross sections were estimated by subtracting the reaction rates of irradiated samples for 25° irradiation from those of 0° irradiation. The measured cross sections were compared with the findings of other studies, evaluated in relation to nuclear data files and the calculated data by Particle and Heavy Ion Transport code System code.