Biomechanical imaging biomarker during chemoradiotherapy predicts treatment response in head and neck squamous cell carcinoma.

Objective. For response-adapted adaptive radiotherapy (R-ART), promising biomarkers are needed to predict post-radiotherapy (post-RT) responses using routine clinical information obtained during RT. In this study, a patient-specific biomechanical model (BM) of the head and neck squamous cell carcinoma (HNSCC) was proposed using the pre-RT maximum standardized uptake value (SUVmax) of18F-fluorodeoxyglucose (FDG) and tumor structural changes during RT as evaluated using computed tomography (CT). In addition, we evaluated the predictive performance of BM-driven imaging biomarkers for the treatment response of patients with HNSCC who underwent concurrent chemoradiotherapy (CCRT).Approach. Patients with histologically confirmed HNSCC treated with definitive CCRT were enrolled in this study. All patients underwent CT two times as follows: before the start of RT (pre-RT) and 3 weeks after the start of RT (mid-RT). Among these patients, 67 patients who underwent positron emission tomography/CT during the pre-RT period were included in the final analysis. The locoregional control (LC), progression-free survival (PFS), and overall survival (OS) prediction performances of whole tumor stress change (TS) between pre- and mid-RT computed using BM were assessed using univariate, multivariate, and Kaplan-Meier survival curve analyses, respectively. Furthermore, performance was compared with the pre and post-RT SUVmax, tumor volume reduction rate (TVRR) during RT, and other clinical prognostic factors.Main results. For both univariate, multivariate, and survival curve analyses, the significant prognostic factors were as follows (p< 0.05): TS and TVRR for LC; TS and pre-RT FDG-SUVmaxfor PFS; and TS only for OS. In addition, for 2 year LC, PFS, and OS prediction, TS showed a comparable predictive performance to post-RT FDG-SUVmax.Significance. BM-driven TS is an effective prognostic factor for tumor treatment response after CCRT. The proposed method can be a feasible functional imaging biomarker that can be acquired during RT using only routine clinical data and may provide useful information for decision-making during R-ART.

Assessing four-dimensional CT stress maps derived from patient-specific biomechanical models of the lung with pulmonary function test data in lung cancer patients.

OBJECTIVE: This study aims to retrospectively compare the stress map of the lung with pulmonary function test (PFT) results in lung cancer patients and to evaluate the potential of the stress map as an imaging biomarker for chronic obstructive pulmonary disease (COPD). METHODS: 25 lung cancer patients with pre-treatment four-dimensional CT (4DCT) and PFT data were retrospectively analysed. PFT metrics were used to diagnose obstructive lung disease. For each patient, forced expiratory volume in 1 s (FEV1 % predicted) and the ratio of FEV1 and forced vital capacity (FEV1/FVC) were recorded. 4DCT and biomechanical model-deformable image registration (BM-DIR) were used to obtain the lung stress map. The relationship between the mean of the total lung stress and PFT data was evaluated, and the COPD classification grade was also evaluated. RESULTS: The mean values of the total lung stress and FEV1 % predicted showed a significant strong correlation [R = 0.833, (p < 0.001)]. The mean values and FEV1/FVC showed a significant strong correlation [R = 0.805, (p < 0.001)]. For the total lung stress, the area under the curve and the optimal cut-off value were 0.94 and 510.8 Pa for the classification of normal or abnormal lung function, respectively. CONCLUSION: This study has demonstrated the potential of lung stress maps based on BM-DIR to accurately assess lung function by comparing them with PFT data. ADVANCES IN KNOWLEDGE: The derivation of stress map directly from 4DCT is novel method. The BM-DIR-based lung stress map can provide an accurate assessment of lung function.

Investigation of the combination of intratumoral and peritumoral radiomic signatures for predicting epidermal growth factor receptor mutation in lung adenocarcinoma.

PURPOSE: We investigated optimal peritumoral size and constructed predictive models for epidermal growth factor receptor (EGFR) mutation. METHODS: A total of 164 patients with lung adenocarcinoma were retrospectively analyzed. Radiomic signatures for the intratumoral region and combinations of intratumoral and peritumoral regions (3, 5, and 7 mm) from computed tomography images were extracted using analysis of variance and least absolute shrinkage. The optimal peritumoral region was determined by radiomics score (rad-score). Intratumoral radiomic signatures with clinical features (IRS) were used to construct predictive models for EGFR mutation. Combinations of intratumoral and 3, 5, or 7 mm-peritumoral signatures with clinical features (IPRS3, IPRS5, and IPRS7, respectively) were also used to construct predictive models. Support vector machine (SVM), logistic regression (LR), and LightGBM models with five-fold cross-validation were constructed, and the receiver operating characteristics were evaluated. Area under the curve (AUC) of the training and test cohorts values were calculated. Brier scores (BS) and decision curve analysis (DCA) were used to evaluate the predictive models. RESULTS: The AUC values of the SVM, LR, and LightGBM models derived from IRS were 0.783 (95% confidence interval: 0.602-0.956), 0.789 (0.654-0.927), and 0.735 (0.613-0.958) for training, and 0.791 (0.641-0.920), 0.781 (0.538-0.930), and 0.734 (0.538-0.930) for test cohort, respectively. Rad-score confirmed that the 3 mm-peritumoral size was optimal (IPRS3), and AUCs values of SVM, LR, and lightGBM models derived from IPRS3 were 0.831 (0.666-0.984), 0.804 (0.622-0.908), and 0.769 (0.628-0.921) for training and 0.765 (0.644-0.921), 0.783 (0.583-0.921), and 0.796 (0.583-0.949) for test cohort, respectively. The BS and DCA of the LR and LightGBM models derived from IPRS3 were better than those from IRS. CONCLUSION: Accordingly, the combination of intratumoral and 3 mm-peritumoral radiomic signatures may be helpful for predicting EGFR mutations.

Comparison of the radiomics-based predictive models using machine learning and nomogram for epidermal growth factor receptor mutation status and subtypes in lung adenocarcinoma.

The purpose of this study is to develop the predictive models for epidermal growth factor receptor (EGFR) mutation status and subtypes [exon 21-point mutation (L858R) and exon 19 deletion mutation (19Del)] and evaluate their clinical usefulness. Total 172 patients with lung adenocarcinoma were retrospectively analyzed. The analysis of variance and the least absolute shrinkage were used for feature selection from plain computed tomography images. Then, radiomic score (rad-score) was calculated for the training and test cohorts. Two machine learning (ML) models with 5-fold were applied to construct the predictive models with rad-score, clinical features, and the combination of rad-score and clinical features. The nomogram was developed using rad-score and clinical features. The prediction performance was evaluated by the area under the receiver operating characteristic curve (AUC). Finally, decision curve analysis (DCA) was performed using the best ML and nomogram models. In the test cohorts, the AUC of the best ML and the nomogram model were 0.73 (95% confidence interval, 0.59-0.87) and 0.79 (0.65-0.92) in the EGFR mutation groups, 0.83 (0.67-0.99) and 0.85 (0.72-0.97) in the L858R mutation groups, as well as 0.77 (0.58-0.97) and 0.77 (0.60-0.95) in the 19Del groups. The DCA showed that the nomogram models have comparable results with ML models. We constructed two predictive models for EGFR mutation status and subtypes. The nomogram models had comparable results to the ML models. Because the superiority of the performance of ML and nomogram models varied depending on the prediction groups, appropriate model selection is necessary.

Assessing liver fibrosis distribution through liver elasticity estimates obtained using a biomechanical model of respiratory motion with magnetic resonance elastography.

Objective. This study aimed to produce a three-dimensional liver elasticity map using the finite element method (FEM) and respiration-induced motion captured by T1-weighted magnetic resonance images (FEM-E-map) and to evaluate whether FEM-E-maps can be an imaging biomarker comparable to magnetic resonance elastography (MRE) for assessing the distribution and severity of liver fibrosis.Approach. We enrolled 14 patients who underwent MRI and MRE. T1-weighted MR images were acquired during shallow inspiration and expiration breath-holding, and the displacement vector field (DVF) between two images was calculated using deformable image registration. FEM-E-maps were constructed using FEM and DVF. First, three Poisson's ratio settings (0.45, 0.49, and 0.499995) were validated and optimized to minimize the difference in liver elasticity between the FEM-E-map and MRE. Then, the whole and regional liver elasticity values estimated using FEM-E-maps were compared with those obtained from MRE using Pearson's correlation coefficients. Spearman rank correlations and chi-square histograms were used to compare the voxel-level elasticity distribution.Main results. The optimal Poisson's ratio was 0.49. Whole liver elasticity estimated using FEM-E-maps was strongly correlated with that measured using MRE (r = 0.96). For regional liver elasticity, the correlation was 0.84 for the right lobe and 0.82 for the left lobe. Spearman analysis revealed a moderate correlation for the voxel-level elasticity distribution between FEM-E-maps and MRE (0.61 ± 0.10). The small chi-square distances between the two histograms (0.11 ± 0.07) indicated good agreement.Significance. FEM-E-maps represent a potential imaging biomarker for visualizing the distribution of liver fibrosis using only T1-weighted images obtained with a common MR scanner, without any additional examination or special elastography equipment. However, additional studies including comparisons with biopsy findings are required to verify the reliability of this method for clinical application.

Predicting programmed death-ligand 1 expression level in non-small cell lung cancer using a combination of peritumoral and intratumoral radiomic features on computed tomography.

In this study, we investigated the possibility of predicting expression levels of programmed death-ligand 1 (PD-L1) using radiomic features of intratumoral and peritumoral tumors on computed tomography (CT) images. We retrospectively analyzed 161 patients with non-small cell lung cancer. We extracted radiomic features for intratumoral and peritumoral regions on CT images. The null importance, least absolute shrinkage, and selection operator model were used to select the optimized feature subset to build the prediction models for the PD-L1 expression level. LightGBM with five-fold cross-validation was used to construct the prediction model and evaluate the receiver operating characteristics. The corresponding area under the curve (AUC) was calculated for the training and testing cohorts. The proportion of ambiguously clustered pairs was calculated based on consensus clustering to evaluate the validity of the selected features. In addition, Radscore was calculated for the training and test cohorts. For expression level of PD-L1 above 1%, prediction models that included radiomic features from the intratumoral region and a combination of radiomic features from intratumoral and peritumoral regions yielded an AUC of 0.83 and 0.87 and 0.64 and 0.74 in the training and test cohorts, respectively. In contrast, the models above 50% prediction yielded an AUC of 0.80, 0.97, and 0.74, 0.83, respectively. The selected features were divided into two subgroups based on PD-L1 expression levels≥50% or≥1%. Radscore was statistically higher for subgroup one than subgroup two when radiomic features for intratumoral and peritumoral regions were combined. We constructed a predictive model for PD-L1 expression level using CT images. The model using a combination of intratumoral and peritumoral radiomic features had a higher accuracy than the model with only intratumoral radiomic features.

Sex determination based on a thoracic vertebra and ribs evaluation using clinical chest radiography.

Our aim was to investigate whether sex can be determined from a combination of geometric features obtained from the 10th thoracic vertebra, 6th rib, and 7th rib. Six hundred chest radiographs (300 males and 300 females) were randomly selected to include patients of six age groups (20s, 30s, 40s, 50s, 60s, and 70s). Each group included 100 images (50 males and 50 females). A total of 14 features, including 7 lengths, 5 indices for the vertebra, and 2 types of widths for ribs, were utilized and analyzed for sex determination. Dominant features contributing to sex determination were selected by stepwise discriminant analysis after checking the variance inflation factors for multicollinearity. The accuracy of sex determination using a combination of the vertebra and ribs was evaluated from the selected features by the stepwise discriminant analysis. The accuracies in each age group were also evaluated in this study. The accuracy of sex determination based on a combination of features of the vertebra and ribs was 88.8% (533/600). This performance was superior to that of the vertebra or ribs only. Moreover, sex determination of subjects in their 20s demonstrated the highest accuracy (96.0%, 96/100). The features selected in the stepwise discriminant analysis included some features in both the vertebra and ribs. These results indicate the usefulness of combined information obtained from the vertebra and ribs for sex determination. We conclude that a combination of geometric characteristics obtained from the vertebra and ribs could be useful for determining sex.

Osteodystrophy in Cholestatic Liver Diseases Is Attenuated by Anti-γ-Glutamyl Transpeptidase Antibody.

BACKGROUND: Cholestatic liver diseases exhibit higher levels of serum γ-glutamyl transpeptidase (GGT) and incidence of secondary osteoporosis. GGT has been identified as a novel bone-resorbing factor that stimulates osteoclast formation. The aim of this study was to elucidate the interaction of elevated GGT levels and cholestatic liver disease-induced bone loss. METHODS: Wistar rats were divided into three groups: sham-operated control (SO) rats, bile duct ligation (BDL) rats, and anti-GGT antibody-treated BDL rats (AGT). Serum GGT level was measured. Bone mineral density (BMD) was analyzed by dual-energy X-ray absorptiometry. Bone morphometric parameters and microarchitectural properties were determined by micro-computed tomography and histomorphometry of the distal metaphysis of femurs. Alterations of bone metabolism-related factors were evaluated by cytokine array. Effects of GGT on osteoblasts or stromal cells were evaluated by RT-PCR, enzyme activity, and mineralization ability. RESULTS: Serum levels of GGT were significantly elevated in the BDL-group. In the BDL group, BMD, bone mass percentage, and osteoblast number were significantly decreased, whereas osteoclast number was significantly increased. These alterations were markedly attenuated in the AGT group. The mRNA levels of vascular endothelial growth factor-A, LPS-induced CXC chemokine, monocyte chemoattractant protein-1, tumor necrosis factor-α interleukin-1ß and receptor activator of nuclear factor-kappa B ligand were upregulated, and those of interferon-γ and osteoprotegerin were downregulated in the GGT-treated stromal cells. Furthermore, GGT inhibited mineral nodule formation and expression of alkaline phosphatase and bone sialo-protein in osteoblastic cells. CONCLUSION: Our results indicate that elevated GGT level is involved in hepatic osteodystrophy through secretion of bone resorbing factor from GGT-stimulated osteoblasts/bone marrow stromal cells. In addition, GGT also possesses suppressive effects on bone formation. Managing elevated GGT levels by anti-GGT antibody may become a novel therapeutic agent for hepatic osteodystrophy in chronic liver diseases.