Run Length Non-uniformity Research Articles

Machine learning approaches to detect hepatocyte chromatin alterations from iron oxide nanoparticle exposure

This study focuses on developing machine learning models to detect subtle alterations in hepatocyte chromatin organization due to Iron (II, III) oxide nanoparticle exposure, hypothesizing that exposure will significantly alter chromatin texture. A total of 2000 hepatocyte nuclear regions of interest (ROIs) from mouse liver tissue were analyzed, and for each ROI, 5 different parameters were calculated: Long Run Emphasis, Short Run Emphasis, Run Length Nonuniformity, and 2 wavelet coefficient energies obtained after the discrete wavelet transform. These parameters served as input for supervised machine learning models, specifically random forest and gradient boosting classifiers. The models demonstrated relatively robust performance in distinguishing hepatocyte chromatin structures belonging to the group exposed to IONPs from the controls. The study's findings suggest that iron oxide nanoparticles induce substantial changes in hepatocyte chromatin distribution and underscore the potential of AI techniques in advancing hepatocyte evaluation in physiological and pathological conditions.

Deep learning-derived splenic radiomics, genomics, and coronary artery disease.

Despite advances in managing traditional risk factors, coronary artery disease (CAD) remains the leading cause of mortality. Circulating hematopoietic cells influence risk for CAD, but the role of a key regulating organ, spleen, is unknown. The understudied spleen is a 3-dimensional structure of the hematopoietic system optimally suited for unbiased radiologic investigations toward novel mechanistic insights. Deep learning-based image segmentation and radiomics techniques were utilized to extract splenic radiomic features from abdominal MRIs of 42,059 UK Biobank participants. Regression analysis was used to identify splenic radiomics features associated with CAD. Genome-wide association analyses were applied to identify loci associated with these radiomics features. Overlap between loci associated with CAD and the splenic radiomics features was explored to understand the underlying genetic mechanisms of the role of the spleen in CAD. We extracted 107 splenic radiomics features from abdominal MRIs, and of these, 10 features were associated with CAD. Genome-wide association analysis of CAD-associated features identified 219 loci, including 35 previously reported CAD loci, 7 of which were not associated with conventional CAD risk factors. Notably, variants at 9p21 were associated with splenic features such as run length non-uniformity. Our study, combining deep learning with genomics, presents a new framework to uncover the splenic axis of CAD. Notably, our study provides evidence for the underlying genetic connection between the spleen as a candidate causal tissue-type and CAD with insight into the mechanisms of 9p21, whose mechanism is still elusive despite its initial discovery in 2007. More broadly, our study provides a unique application of deep learning radiomics to non-invasively find associations between imaging, genetics, and clinical outcomes.

Texture analysis of ultrasonography to differentiate metastatic from nonmetastatic cervical lymph nodes in mandibular gingival squamous cell carcinoma

AbstractAimThis study aimed to distinguish between metastatic and nonmetastatic cervical lymph nodes (CLNs) in patients with mandibular gingival squamous cell carcinoma using texture analysis of lymph nodes on ultrasonography images.MethodsA total of 29 metastatic and 34 nonmetastatic CLNs were selected for this study. Fifty‐six texture characteristics were retrieved from the ultrasonography images using the LIFEx software.The Mann–Whitney U test was used to evaluate measurable differences in texture characteristics between metastatic and nonmetastatic CLNs. The ability of the surface features to distinguish between metastatic and nonmetastatic CLNs was illustrated using receiver operating characteristic curves. Youden's J statistic was used to determine the receiver operating characteristic curve cutoff values with maximal sensitivity and specificity.ResultsGray‐level nonuniformity, run length nonuniformity, coarseness, strength, and zone size nonuniformity showed the most notable contrast between the metastatic and nonmetastatic CLNs (p < 0.001). Zone size nonuniformity showed an area under the curve value of 0.827, sensitivity of 0.645, and specificity of 0.972 at the cutoff value of 1010.538.ConclusionQualitative characteristics may be the optimal surface features to distinguish between metastatic and nonmetastatic CLNs and predict CLN metastasis in patients with mandibular gingival squamous cell carcinoma.

Comparative Analysis of LSTM and BiLSTM in Image Detection Processing

Tuberculosis is an infectious disease and requires serious treatment. Extrapulmonary Tuberculosis is detected using a microscope. Currently it will take a long time because the fluid preparations are viewed in a microscope one by one carefully and in the fluid preparations there are 150 fields of vision. Examination for Extra Pulmonary Tuberculosis by culture takes between 1-2 weeks or even more. Examination by biopsy will take a long time because the fluid preparations are looked at carefully under the microscope one by one. The image of Tuberculosis is expressed if in the image there is a bacillus object in red, and it turns out that apart from the bacillus object there are other objects also in red. So that examinations for tuberculosis can be more efficient, examinations using computer technology are needed. This research aims to compare the Long Short-Term Memory (LSTM) and Bidirectional Long Short-Term Memory (BiLSTM) classification methods in the detection of extra-pulmonary tuberculosis disease to obtain better accuracy results. This research carried out HSI color space transformation, segmentation using global thresholding, feature extraction using 13 features based on shape and texture using the Correlation Based Feature Selection (CFS) feature selection method. The results show that BiLSTM has the best accuracy with a value of 88.40% at the number of features = 3, namely Short Run High Gray-Level Emphasis, Run Length Nonuniformity, Minor axis length, while LSTM produces an accuracy of 63.19% at the number of features = 5. BiLSTM is capable of detecting opposite features, meaning that BiLSTM can detect opposite features in data sequences and BiLSTM's ability to understand multiple contexts, so it tends to provide more accurate results in some data classification tasks.

Three-dimensional texture analyses of multi-quantitative relaxation time maps for evaluating cartilage repair with the treatment of allogeneic human adipose-derived mesenchymal progenitor cells

BackgroundTo explore the ability of three-dimensional texture analyses based on gray-level run-length matrix (GLRLM) for examining the spatial distribution of pixel values on magnetic resonance imaging (MRI) relaxation time maps and detecting the compositional variation of cartilage repair following treatment with allogeneic human adipose-derived mesenchymal progenitor cells (haMPCs). MethodsParticipants with knee osteoarthritis were randomly divided into three groups with intra-articular haMPCs injections: low-, medium-, and high-dose groups. We analyzed five GLRLM parameters in the T1rho, T2 and T2star maps, including run length non-uniformity (RLNonUni), gray-level non-uniformity (GLevNonU), long run emphasis (LngREmph), short run emphasis (ShrtREmp), and fraction of images in runs. We used the relative D values (the ratio of difference values to baseline) as the objective to avoid errors caused by individual differences. We calculated the two-tailed Pearson's linear correlation coefficient (r) to investigate the correlations of the texture parameters with the Western Ontario and McMaster Universities Osteoarthritis Index (WOMAC) scores. ResultsCompared with the base time, significant reduction of WOMAC score was observed in both high and medium doses groups at terminal time, indicating relief of pain symptoms in high and medium groups with the treatment of allogeneic haMPCs. Significant differences were observed in the GLRLM parameters of cartilage MR relaxation time maps in different doses groups. In both T1rho and T2 relaxation time maps, the high-dose group showed significant increases in relative D values of RLNonUni, GLevNonU, LngREmph and ShrtREmp, which indicated significant changes in the uniformity of relaxation time maps. For T2star map, GLRLM parameters such as GLevNonU and ShrtREmp, especially LngREmph, showed significant increases in relative D values in high-dose group. Among all GLRLM features, LngREmph of three relaxation time maps had performed excellent linear correlations with WOMAC scores. ConclusionsTexture analysis of the cartilage may allow the detection of compositional variation in cartilage repair with the treatment of allogeneic haMPCs. This technique displays potential applications in understanding the mechanism of stem cell repair of the cartilage and assessing the treatment response.

Prognostic 18F-FDG Radiomic Features in Advanced High-Grade Serous Ovarian Cancer.

High-grade serous ovarian cancer (HGSOC) is an aggressive disease with different clinical outcomes and poor prognosis. This could be due to tumor heterogeneity. The 18F-FDG PET radiomic parameters permit addressing tumor heterogeneity. Nevertheless, this has been not well studied in ovarian cancer. The aim of our work was to assess the prognostic value of pretreatment 18F-FDG PET radiomic features in patients with HGSOC. A review of 36 patients diagnosed with advanced HGSOC between 2016 and 2020 in our center was performed. Radiomic features were obtained from pretreatment 18F-FDGPET. Disease-free survival (DFS) and overall survival (OS) were calculated. Optimal cutoff values with receiver operating characteristic curve/median values were used. A correlation between radiomic features and DFS/OS was made. The mean DFS was 19.6 months and OS was 37.1 months. Total Lesion Glycolysis (TLG), GLSZM_ Zone Size Non-Uniformity (GLSZM_ZSNU), and GLRLM_Run Length Non-Uniformity (GLRLM_RLNU) were significantly associated with DFS. The survival-curves analysis showed a significant difference of DSF in patients with GLRLM_RLNU > 7388.3 versus patients with lower values (19.7 months vs. 31.7 months, p = 0.035), maintaining signification in the multivariate analysis (p = 0.048). Moreover, Intensity-based Kurtosis was associated with OS (p = 0.027). Pretreatment 18F-FDG PET radiomic features GLRLM_RLNU, GLSZM_ZSNU, and Kurtosis may have prognostic value in patients with advanced HGSOC.

Long-term reproducibility of opportunistically assessed vertebral bone mineral density and texture features in routine clinical multi-detector computed tomography using an automated segmentation framework.

To investigate reproducibility of texture features and volumetric bone mineral density (vBMD) extracted from trabecular bone in the thoracolumbar spine in routine clinical multi-detector computed tomography (MDCT) data in a single scanner environment. Patients who underwent two routine clinical thoraco-abdominal MDCT exams at a single scanner with a time interval of 6 to 26 months (n=203, 131 males; time interval mean, 13 months; median, 12 months) were included in this observational study. Exclusion criteria were metabolic and hematological disorders, bone metastases, use of bone-active medications, and history of osteoporotic vertebral fractures (VFs) or prior diagnosis of osteoporosis. A convolutional neural network (CNN)-based framework was used for automated spine labeling and segmentation (T5-L5), asynchronous Hounsfield unit (HU)-to-BMD calibration, and correction for the intravenous contrast medium phase. Vertebral vBMD and six texture features [varianceglobal, entropy, short-run emphasis (SRE), long-run emphasis (LRE), run-length non-uniformity (RLN), and run percentage (RP)] were extracted for mid- (T5-T8) and lower thoracic (T9-T12), and lumbar vertebrae (L1-L5), respectively. Relative annual changes were calculated in texture features and vBMD for each vertebral level and sorted by sex, and changes were checked for statistical significance (P<0.05) using paired t-tests. Root mean square coefficient of variation (RMSCV) and root mean square error (RMSE) were calculated as measures of variability. SRE, LRE, RLN, and RP exhibited substantial reproducibility with RMSCV-values below 2%, for both sexes and at all spine levels, while vBMD was less reproducible (RMSCV =11.9-16.2%). Entropy showed highest variability (RMSCV =4.34-7.69%) due to statistically significant increases [range, mean ± standard deviation: (4.40±5.78)% to (8.36±8.66)%, P<0.001]. RMSCV of varianceglobal ranged from 1.60% to 3.03%. Opportunistic assessment of texture features in a single scanner environment using the presented CNN-based framework yields substantial reproducibility, outperforming vBMD reproducibility. Lowest scan-rescan variability was found for higher-order texture features. Further studies are warranted to determine, whether microarchitectural changes to the trabecular bone may be assessed through texture features.

Differentiation of benign and malignant vertebral fractures using a convolutional neural network to extract CT-based texture features

PurposeTo assess the diagnostic performance of three-dimensional (3D) CT-based texture features (TFs) using a convolutional neural network (CNN)-based framework to differentiate benign (osteoporotic) and malignant vertebral fractures (VFs).MethodsA total of 409 patients who underwent routine thoracolumbar spine CT at two institutions were included. VFs were categorized as benign or malignant using either biopsy or imaging follow-up of at least three months as standard of reference. Automated detection, labelling, and segmentation of the vertebrae were performed using a CNN-based framework (https://anduin.bonescreen.de). Eight TFs were extracted: Varianceglobal, Skewnessglobal, energy, entropy, short-run emphasis (SRE), long-run emphasis (LRE), run-length non-uniformity (RLN), and run percentage (RP). Multivariate regression models adjusted for age and sex were used to compare TFs between benign and malignant VFs.ResultsSkewnessglobal showed a significant difference between the two groups when analyzing fractured vertebrae from T1 to L6 (benign fracture group: 0.70 [0.64–0.76]; malignant fracture group: 0.59 [0.56–0.63]; and p = 0.017), suggesting a higher skewness in benign VFs compared to malignant VFs.ConclusionThree-dimensional CT-based global TF skewness assessed using a CNN-based framework showed significant difference between benign and malignant thoracolumbar VFs and may therefore contribute to the clinical diagnostic work-up of patients with VFs.

A machine learning-based sonomics for prediction of thyroid nodule malignancies.

This study aims to use ultrasound derived features as biomarkers to assess the malignancy of thyroid nodules in patients who were candidates for FNA according to the ACR TI-RADS guidelines. Two hundred and ten patients who met the selection criteria were enrolled in the study and subjected to ultrasound-guided FNA of thyroid nodules. Different radiomics features were extracted from sonographic images, including intensity, shape, and texture feature sets. Least Absolute Shrinkage and Selection Operator (LASSO), Minimum Redundancy Maximum Relevance (MRMR), and Random Forests/Extreme Gradient Boosting Machine (XGBoost) algorithms were used for feature selection and classification of the univariate and multivariate modeling, respectively. Evaluation of models performed using accuracy, sensitivity, specificity, and area under the receiver operating characteristic curve (AUC). In the univariate analysis, Gray Level Run Length Matrix - Run-Length Non-Uniformity (GLRLM-RLNU) and gray-level zone length matrix - Run-Length Non-Uniformity (GLZLM-GLNU) (both with an AUC of 0.67) were top-performing for predicting nodules malignancy. In the multivariate analysis of the training dataset, the AUC of all combinations of feature selection algorithms and classifiers was 0.99, and the highest sensitivity was for XGBoost classifier and MRMR feature selection algorithms (0.99). Finally, the test dataset was used to evaluate our model in which XGBoost classifier with MRMR and LASSO feature selection algorithms had the highest performance (AUC = 0.95). Ultrasound-extracted features can be used as non-invasive biomarkers for thyroid nodules' malignancy prediction.

Left Ventricular Vertical Run-Length Nonuniformity MRI Adds Prognostic Value to MACE in Patients with End-Stage Renal Disease.

Vertical run-length nonuniformity (VRLN) is a texture feature representing heterogeneity within native T1 images and reflects the extent of cardiac fibrosis. In uremic cardiomyopathy, interstitial fibrosis was the major histological alteration. The prognostic value of VRLN in patients with end-stage renal disease (ESRD) remains unclear. To evaluate the prognostic value of VRLN MRI in patients with ESRD. Prospective. A total of 127 ESRD patients (30 participants in the major adverse cardiac events, MACE group). 3.0 T/steady-state free precession sequence, modified Look-Locker imaging. MRI image qualities were assessed by three independent radiologists. VRLN values were measured in the myocardium on the mid-ventricular short-axis slice of T1 mapping. Left ventricular (LV) mass, LV end-diastolic and end-systolic volume, as well as LV global strain cardiac parameters were measured. The primary endpoint was the incident of MACE from enrollment time to January 2023. MACE is a composite endpoint consisting of all-cause mortality, acute myocardial infarction, stroke, heart failure hospitalization, and life-threatening arrhythmia. Cox proportional-hazards regression was performed to test whether VRLN independently correlated with MACE. The intraclass correlation coefficients of VRLN were calculated to evaluate intraobserver and interobserver reproducibility. The C-index was computed to examine the prognostic value of VRLN. P-value <0.05 were considered statistically significant. Participants were followed for a median of 26 months. VRLN, age, LV end-systolic volume index, and global longitudinal strain remained significantly associated with MACE in the multivariable model. Adding VRLN to a baseline model containing clinical and conventional cardiac MRI parameters significantly improved the accuracy of the predictive model (C-index of the baseline model: 0.781 vs. the model added VRLN: 0.814). VRLN is a novel marker for risk stratification toward MACE in patients with ESRD, superior to native T1 mapping and LV ejection fraction. 2 TECHNICAL EFFICACY STAGE: 2.

Efficacy of texture analysis of ultrasonographic images in the differentiation of metastatic and non-metastatic cervical lymph nodes in patients with squamous cell carcinoma of the tongue

We investigated the efficacy of using texture analysis of ultrasonographic images of the cervical lymph nodes of patients with squamous cell carcinoma of the tongue to differentiate between metastatic and non-metastatic lymph nodes. We analyzed 32 metastatic and 28 non-metastatic lymph nodes diagnosed by histopathologic examination on presurgical US images. Using the LIFEx texture analysis program, we extracted 36 texture features from the images and calculated the statistical significance of differences in texture features between metastatic and non-metastatic lymph nodes using the t test. To assess the diagnostic ability of the significantly different texture features to discriminate between metastatic and non-metastatic nodes, we performed receiver operating characteristic curve analysis and calculated the area under the curve. We set the cutoff points that maximized the sensitivity and specificity for each curve according to the Youden J statistic. We found that 20 texture features significantly differed between metastatic and non-metastatic lymph nodes. Among them, only the gray-level run length matrix feature of run length non-uniformity and the gray-level zone length matrix features of gray-level non-uniformity and zone length non-uniformity showed an excellent ability to discriminate between metastatic and non-metastatic lymph nodes as indicated by the area under the curve and the sum of sensitivity and specificity. Analysis of the texture features of run length non-uniformity, gray-level non-uniformity, and zone length non-uniformity values allows for differentiation between metastatic and non-metastatic lymph nodes, with the use of gray-level non-uniformity appearing to be the best means of predicting metastatic lymph nodes.

FDG PET texture indices as imaging biomarkers for epidermal growth factor receptor mutation status in lung adenocarcinoma

Identifying the epidermal growth factor receptor (EGFR) mutation status is important for the optimal treatment of patients with EGFR mutations. We investigated the relationship between 18F-fluorodeoxyglucose (FDG) positron emission tomography (PET) texture indices and EGFR mutation status in patients with newly diagnosed lung adenocarcinoma. We retrospectively analyzed data of patients with newly diagnosed lung adenocarcinoma who underwent pretreatment FDG PET/computed tomography and EGFR mutation testing between August 2014 and November 2020. Patients were divided into mutated EGFR and wild-type EGFR groups. The maximum standardized uptake value (SUVmax) and 31 texture indices for the primary tumor were calculated from PET images and compared between the two groups. Of the 66 patients included, 22 had mutated EGFR and 44 had wild-type EGFR. The SUVmax did not significantly differ between the two groups. Among the 31 evaluated texture indices, the following five showed a statistically significant difference between the groups: correlation (P = 0.003), gray-level nonuniformity for run (P = 0.042), run length nonuniformity (P = 0.02), coarseness (P = 0.006), and gray-level nonuniformity for zone (P = 0.04). Based on the preliminary results of this study in a small patient population, FDG PET texture indices may be potential imaging biomarkers for the EGFR mutation status in patients with newly diagnosed lung adenocarcinoma.

Efficacy of magnetic resonance imaging texture features of the lateral pterygoid muscle in distinguishing rheumatoid arthritis and osteoarthritis of the temporomandibular joint.

The aim of this study was to assess whether magnetic resonance imaging (MRI) texture features of the lateral pterygoid muscle can distinguish between rheumatoid arthritis (RA) and osteoarthritis (OA) of the temporomandibular joint (TMJ). The authors extracted 279 texture features from 36 patients with RA and OA from the region of interest set for the lateral pterygoid muscle on short tau inversion recovery (STIR) images using MaZda Ver.3.3. A total of 10 texture features were selected using Fisher's coefficients, as well as probability of error and average correlation coefficients. Data observed to have a non-normal distribution using the Kolmogorov-Smirnov test were compared using the Mann-Whitney U-test. Receiver operating characteristic (ROC) curves were used to assess the ability of the 10 texture features to distinguish RA and OA of the TMJ. A total of 10 features (5 Correlation, 3 Run Length Nonuniformity, 1 Sigma, and 1 Teta) were selected from 279 texture features. These texture features revealed significant differences between the RA and OA groups (p < 0.01). The sensitivity, specificity, accuracy, and area under the ROC curve of the texture features for distinguishing RA from OA were 0.78-0.94, 0.89-1.0, 0.86-0.92, and 0.89-0.95, respectively. MRI texture analysis of the lateral pterygoid muscle may be useful for distinguishing between RA and OA of the TMJ.

Systematic analysis of changes in radiomics features during dynamic breast-MRI: Evaluation of specific biomarkers

ObjectivesIn this retrospective, single-center study we investigate the changes of radiomics features during dynamic breast-MRI for healthy tissue compared to benign and malignant lesions. Methods60 patients underwent breast-MRI using a dynamic 3D gradient-echo sequence. Changes of 34 texture features (TF) in 30 benign and 30 malignant lesions were calculated for 5 dynamic datasets and corresponding 4 subtraction datasets. Statistical analysis was performed with ANOVA, and systematic changes in features were described by linear and polynomial regression models. ResultsANOVA revealed significant differences (p < 0.05) between normal tissue and lesions in 13 TF, compared to 9 TF between benign and malignant lesions. Most TF showed significant differences in early dynamic and subtraction datasets. TF associated with homogeneity were suitable to discriminate between healthy parenchyma and lesions, whereas run-length features were more suitable to discriminate between benign and malignant lesions. Run length nonuniformity (RLN) was the only feature able to distinguish between all three classes with an AUC of 88.3%. Characteristic changes were observed with a systematic increase or decrease for most TF with mostly polynomial behavior. Slopes showed earlier peaks in malignant lesions, compared to benign lesions. Mean values for the coefficient of determination were higher during subtraction sequences, compared to dynamic sequences (benign: 0.98 vs 0. 72; malignant: 0.94 vs 0.74). ConclusionsTF of breast lesions follow characteristic patterns during dynamic breast-MRI, distinguishing benign from malignant lesions. Early dynamic and subtraction datasets are particularly suitable for texture analysis in breast-MRI. Features associated with tissue homogeneity seem to be indicative of benign lesions.

Radiomics in photon-counting dedicated breast CT: potential of texture analysis for breast density classification

BackgroundWe investigated whether features derived from texture analysis (TA) can distinguish breast density (BD) in spiral photon-counting breast computed tomography (PC-BCT).MethodsIn this retrospective single-centre study, we analysed 10,000 images from 400 PC-BCT examinations of 200 patients. Images were categorised into four-level density scale (a–d) using Breast Imaging Reporting and Data System (BI-RADS)-like criteria. After manual definition of representative regions of interest, 19 texture features (TFs) were calculated to analyse the voxel grey-level distribution in the included image area. ANOVA, cluster analysis, and multinomial logistic regression statistics were used. A human readout then was performed on a subset of 60 images to evaluate the reliability of the proposed feature set.ResultsOf the 19 TFs, 4 first-order features and 7 second-order features showed significant correlation with BD and were selected for further analysis. Multinomial logistic regression revealed an overall accuracy of 80% for BD assessment. The majority of TFs systematically increased or decreased with BD. Skewness (rho -0.81), as a first-order feature, and grey-level nonuniformity (GLN, -0.59), as a second-order feature, showed the strongest correlation with BD, independently of other TFs. Mean skewness and GLN decreased linearly from density a to d. Run-length nonuniformity (RLN), as a second-order feature, showed moderate correlation with BD, but resulted in redundant being correlated with GLN. All other TFs showed only weak correlation with BD (range -0.49 to 0.49, p < 0.001) and were neglected.ConclusionTA of PC-BCT images might be a useful approach to assess BD and may serve as an observer-independent tool.

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

To develop a radiomics filtering technique for characterizing spatial-encoded regional pulmonary ventilation information on lung computed tomography (CT). 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. 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. 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.

Pretherapy 18F-fluorodeoxyglucose positron emission tomography/computed tomography robust radiomic features predict overall survival in non-small cell lung cancer.

To extract robust radiomic features from staging positron emission tomography/computed tomography (18F- fluroodeoxyglucose PET/CT) in patients with non-small cell lung cancer from different segmentation methods and to assess their association with 2-year overall survival. Eighty-one patients with stage I-IV non-small cell lung cancer were included. All patients underwent a pretherapy 18F-FDG PET/CT. Primary tumors were delineated using four different segmentation methods: method 1, manual; method 2: manual with peripheral 1 mm erosion; method 3: absolute threshold at standardized uptake value (SUV) 2.5; and method 4: relative threshold at 40% SUVmax. Radiomic features from each method were extracted using Image Biomarker Standardization Initiative-compliant process. The study cohort was divided into two groups (exploratory and testing) in a ratio of 1:2 (n = 25 and n = 56, respectively). Exploratory cohort was used to identify robust radiomic features, defined as having a minimum concordance correlation coefficient ≥0.75 among all the 4-segmentation methods. The resulting texture features were evaluated for association with 2-year overall survival in the testing cohort (n = 56). All patients in the testing cohort had a follow-up for 2 years from the date of staging 18F-FDG PET/CT scan or till death. Cox proportional hazard models were used to evaluate the independent prognostic factors. Exploratory and validation cohorts were equivalent regarding their basic characteristics (age, sex, and tumor stage). Ten radiomic features were deemed robust to the described four segmentation methods: SUV SD, SUVmax, SUVQ3, SUVpeak in 0.5 ml, total lesion glycolysis, histogram entropy log 2, histogram entropy log 10, histogram energy uniformity, gray level run length matrix-gray level non-uniformity, and gray level zone length matrix-gray level non-uniformity. At the end of 2-year follow-up, 41 patients were dead and 15 were still alive (overall survival = 26.8%; median survival = 14.7 months, 95% confidence interval: 10.2-19.2 months). Three texture features, regardless the segmentation method, were associated with 2-year overall survival: total lesion glycolysis, gray level run length matrix_gray level non-uniformity, and gray level zone length matrix_run-length non-uniformity. In the final Cox-regression model: total lesion glycolysis, and gray level zone length matrix_gray level non-uniformity were independent prognostic factors. The quartiles from the two features were combined with clinical staging in a prognostic model that allowed better risk stratification of patients for overall survival. Ten radiomic features were robust to segmentation methods and two of them (total lesion glycolysis and gray level zone length matrix_gray level non-uniformity) were independently associated with 2-year overall survival. Together with the clinical staging, these features could be utilized towards improved risk stratification of lung cancer patients.

Gender-, Age- and Region-Specific Characterization of Vertebral Bone Microstructure Through Automated Segmentation and 3D Texture Analysis of Routine Abdominal CT.

PurposeTo identify long-term reproducible texture features (TFs) of spinal computed tomography (CT), and characterize variations with regard to gender, age and vertebral level using our automated quantification framework.MethodsWe performed texture analysis (TA) on baseline and follow-up CT (follow-up duration: 30–90 days) of 21 subjects (8 females, 13 males, age at baseline 61.2 ± 9.2 years) to determine long-term reproducibility. TFs with a long-term reproducibility error Δrel<5% were further analyzed for an association with age and vertebral level in a cohort of 376 patients (129 females, 247 males, age 62.5 ± 9.2 years). Automated analysis comprised labeling and segmentation of vertebrae into subregions using a convolutional neural network, calculation of volumetric bone mineral density (vBMD) with asynchronous calibration and TF extraction. Varianceglobal measures the spread of the gray-level distribution in an image while Entropy reflects the uniformity of gray-levels. Short-run emphasis (SRE), Long-run emphasis (LRE), Run-length non-uniformity (RLN) and Run percentage (RP) contain information on consecutive voxels of a particular grey-level, or grey-level range, in a particular direction. Long runs (LRE) represent coarse texture while short runs (SRE) represent fine texture. RLN reflects similarities in the length of runs while RP reflects distribution and homogeneity of runs with a specific direction.ResultsSix of the 24 extracted TFs had Δrel<5% (Varianceglobal, Entropy, SRE, LRE, RLN, RP), and were analyzed further in 4716 thoracolumbar vertebrae. Five TFs (Varianceglobal,SRE,LRE, RLN,RP) showed a significant difference between genders (p<0.001), potentially being caused by a finer and more directional vertebral trabecular microstructure in females compared to males. Varianceglobal and Entropy showed a significant increase from the thoracic to the lumbar spine (p<0.001), indicating a higher degree and earlier initiation of trabecular microstructure deterioration at lower spinal levels. The four higher-order TFs showed significant variations between spine regions without a clear directional gradient (p ≤ 0.001-0.012). No TF showed a clear age dependence. vBMD differed significantly between genders, age groups and spine regions (p ≤ 0.001–0.002).ConclusionLong-term reproducible CT-based TFs of the thoracolumbar spine were established and characterized in a predominantly older adult study population. The gender-, age- and vertebral-level-specific values may serve as foundation for osteoporosis diagnostics and facilitate future studies investigating vertebral microstructure.

Predicting Progression-Free Survival Using Primary and Nodal Radiomic Features in Head and Neck Cancer

Imaging quantified by radiomics has yielded reproducible predictions of both oncologic outcomes and treatment-related toxicity. Much of the correlation with oncologic outcomes has focused on image features from the primary tumor. However, both primary tumor and nodal features can impact outcomes, and it is known that differences in gene expression in primary tumors and nodal metastases provide distinct information regarding survival in head and neck squamous cell carcinomas (HNSCC). Nevertheless, current radiomic biomarkers often make no distinction between primary and nodal radiomic features. Therefore, we hypothesize that combining CT image features of primary tumors and nodal metastases in HNSCC will improve our ability to predict progression-free survival (PFS) compared to primary features alone.Pre-treatment CT images were prospectively collected from 115 HNSCC patients undergoing RT at our institution from 2008-2018. PFS was determined at time of last follow-up prior to 2018. Clinical data were systemically captured from patient medical charts. 107 radiomic features were extracted for the primary tumor and the ipsilateral/contralateral nodal regions. Features were pre-selected by literature review for those possibly significant in predicting loco-regional control of HNSCC. Tumor volume, correlation of gray level co-occurrence matrix (Corr-GLCM), run length non-uniformity of gray level of gray level run length matrix (RLN-GLRLM), and coarseness of neighboring gray tone difference matrix (Coarse-NGTDM) for the primary contours; and tumor volume and Corr-GLCM of the nodal contours were included in the model. A regression analysis of pre-selected features was performed using ridge regularization in the training set (n = 93 primary; n = 48 primary & nodal), and a generalized linear model was built with repeated ten-fold cross validation of patients treated from 2008-2017. Model performance was evaluated in a test set (patients treated 2017-2018) using the area under the receiver operating characteristic curve (AUC). Models were built separately for patients with primary tumor only and primary + nodal radiomic features.The AUC/sensitivity/specificity on the test set for the model with primary features (n = 22) was 0.91/1.00/0.75, and for the model with primary and nodal features (n = 18) was 0.75/0.60/0.92. The features with the highest coefficients in the combined model included primary Coarse-NGTDM (β = -0.21), primary Corr-GLCM (β = 0.0082), and nodal Corr-GLCM (β = -0.00041).In a small cohort of HNSCC patients, combining nodal features with primary features did not improve the prediction performance compared to primary features alone. This work suggests that primary tumor radiomic features may influence PFS more than nodal features. Future work is required in a larger cohort to validate these findings.

Establishing ADC-Based Histogram and Texture Features for Early Treatment-Induced Changes in Head and Neck Squamous Cell Carcinoma.

The purpose of this study was to assess baseline variability in histogram and texture features derived from apparent diffusion coefficient (ADC) maps from diffusion-weighted MRI (DW-MRI) examinations and to identify early treatment-induced changes to these features in patients with head and neck squamous cell carcinoma (HNSCC) undergoing definitive chemoradiation. Patients with American Joint Committee on Cancer Stage III–IV (7th edition) HNSCC were prospectively enrolled on an IRB-approved study to undergo two pre-treatment baseline DW-MRI examinations, performed 1 week apart, and a third early intra-treatment DW-MRI examination during the second week of chemoradiation. Forty texture and six histogram features were derived from ADC maps. Repeatability of the features from the baseline ADC maps was assessed with the intra-class correlation coefficient (ICC). A Wilcoxon signed-rank test compared average baseline and early treatment feature changes. Data from nine patients were used for this study. Comparison of the two baseline ADC maps yielded 11 features with an ICC ≥ 0.80, indicating that these features had excellent repeatability: Run Gray-Level Non-Uniformity, Coarseness, Long Zone High Gray-Level, Variance (Histogram Feature), Cluster Shade, Long Zone, Variance (Texture Feature), Run Length Non-Uniformity, Correlation, Cluster Tendency, and ADC Median. The Wilcoxon signed-rank test resulted in four features with significantly different early treatment-induced changes compared to the baseline values: Run Gray-Level Non-Uniformity (p = 0.005), Run Length Non-Uniformity (p = 0.005), Coarseness (p = 0.006), and Variance (Histogram) (p = 0.006). The feasibility of histogram and texture analysis as a potential biomarker is dependent on the baseline variability of each metric, which disqualifies many features.