Radiomics and artificial intelligence for precision medicine in lung cancer treatment.

Lung cancer is the leading cause of cancer-related deaths worldwide. It exhibits, at the mesoscopic scale, phenotypic characteristics that are generally indiscernible to the human eye but can be captured non-invasively on medical imaging as radiomic features, which can form a high dimensional data space amenable to machine learning. Radiomic features can be harnessed and used in an artificial intelligence paradigm to risk stratify patients, and predict for histological and molecular findings, and clinical outcome measures, thereby facilitating precision medicine for improving patient care. Compared to tissue sampling-driven approaches, radiomics-based methods are superior for being non-invasive, reproducible, cheaper, and less susceptible to intra-tumoral heterogeneity. This review focuses on the application of radiomics, combined with artificial intelligence, for delivering precision medicine in lung cancer treatment, with discussion centered on pioneering and groundbreaking works, and future research directions in the area.

Radiogenomics of C9orf72 Expansion Carriers Reveals Global Transposable Element Derepression and Enables Prediction of Thalamic Atrophy and Clinical Impairment.

Hexanucleotide repeat expansion (HRE) within C9orf72 is the most common genetic cause of frontotemporal dementia (FTD). Thalamic atrophy occurs in both sporadic and familial FTD but is thought to distinctly affect HRE carriers. Separately, emerging evidence suggests widespread derepression of transposable elements (TEs) in the brain in several neurodegenerative diseases, including C9orf72 HRE-mediated FTD (C9-FTD). Whether TE activation can be measured in peripheral blood and how the reduction in peripheral C9orf72 expression observed in HRE carriers relates to atrophy and clinical impairment remain unknown. We used FreeSurfer software to assess the effects of C9orf72 HRE and clinical diagnosis (n = 78 individuals, male and female) on atrophy of thalamic nuclei. We also generated a novel, human, whole-blood RNA-sequencing dataset to determine the relationships among peripheral C9orf72 expression, TE activation, thalamic atrophy, and clinical severity (n = 114 individuals, male and female). We confirmed global thalamic atrophy and reduced C9orf72 expression in HRE carriers. Moreover, we identified disproportionate atrophy of the right mediodorsal lateral nucleus in HRE carriers and showed that C9orf72 expression associated with clinical severity, independent of thalamic atrophy. Strikingly, we found global peripheral activation of TEs, including the human endogenous LINE-1 element L1HS L1HS levels were associated with atrophy of multiple pulvinar nuclei, a thalamic region implicated in C9-FTD. Integration of peripheral transcriptomic and neuroimaging data from human HRE carriers revealed atrophy of specific thalamic nuclei, demonstrated that C9orf72 levels relate to clinical severity, and identified marked derepression of TEs, including L1HS, which predicted atrophy of FTD-relevant thalamic nuclei.SIGNIFICANCE STATEMENT Pathogenic repeat expansion in C9orf72 is the most frequent genetic cause of FTD and amyotrophic lateral sclerosis (ALS; C9-FTD/ALS). The clinical, neuroimaging, and pathologic features of C9-FTD/ALS are well characterized, whereas the intersections of transcriptomic dysregulation and brain structure remain largely unexplored. Herein, we used a novel radiogenomic approach to examine the relationship between peripheral blood transcriptomics and thalamic atrophy, a neuroimaging feature disproportionately impacted in C9-FTD/ALS. We confirmed reduction of C9orf72 in blood and found broad dysregulation of transposable elements-genetic elements typically repressed in the human genome-in symptomatic C9orf72 expansion carriers, which associated with atrophy of thalamic nuclei relevant to FTD. C9orf72 expression was also associated with clinical severity, suggesting that peripheral C9orf72 levels capture disease-relevant information.

Radiogenomics and Texture Analysis to Detect von Hippel-Lindau (VHL) Mutation in Clear Cell Renal Cell Carcinoma.

Radiogenomics, a burgeoning field in biomedical research, explores the correlation between imaging features and genomic data, aiming to link macroscopic manifestations with molecular characteristics. In this review, we examine existing radiogenomics literature in clear cell renal cell carcinoma (ccRCC), the predominant renal cancer, and von Hippel-Lindau (VHL) gene mutation, the most frequent genetic mutation in ccRCC. A thorough examination of the literature was conducted through searches on the PubMed, Medline, Cochrane Library, Google Scholar, and Web of Science databases. Inclusion criteria encompassed articles published in English between 2014 and 2022, resulting in 10 articles meeting the criteria out of 39 initially retrieved articles. Most of these studies applied computed tomography (CT) images obtained from open source and institutional databases. This literature review investigates the role of radiogenomics, with and without texture analysis, in predicting VHL gene mutation in ccRCC patients. Radiogenomics leverages imaging modalities such as CT and magnetic resonance imaging (MRI), to analyze macroscopic features and establish connections with molecular elements, providing insights into tumor heterogeneity and biological behavior. The investigations explored diverse mutations, with a specific focus on VHL mutation, and applied CT imaging features for radiogenomic analysis. Moreover, radiomics and machine learning techniques were employed to predict VHL gene mutations based on CT features, demonstrating promising results. Additional studies delved into the relationship between VHL mutation and body composition, revealing significant associations with adipose tissue distribution. The review concludes by highlighting the potential role of radiogenomics in guiding targeted and selective therapies.

Conditional generative adversarial network driven radiomic prediction of mutation status based on magnetic resonance imaging of breast cancer.

BACKGROUND: Breast Cancer (BC) is a highly heterogeneous and complex disease. Personalized treatment options require the integration of multi-omic data and consideration of phenotypic variability. Radiogenomics aims to merge medical images with genomic measurements but encounter challenges due to unpaired data consisting of imaging, genomic, or clinical outcome data. In this study, we propose the utilization of a well-trained conditional generative adversarial network (cGAN) to address the unpaired data issue in radiogenomic analysis of BC. The generated images will then be used to predict the mutations status of key driver genes and BC subtypes. METHODS: We integrated the paired MRI and multi-omic (mRNA gene expression, DNA methylation, and copy number variation) profiles of 61 BC patients from The Cancer Imaging Archive (TCIA) and The Cancer Genome Atlas (TCGA). To facilitate this integration, we employed a Bayesian Tensor Factorization approach to factorize the multi-omic data into 17 latent features. Subsequently, a cGAN model was trained based on the matched side-view patient MRIs and their corresponding latent features to predict MRIs for BC patients who lack MRIs. Model performance was evaluated by calculating the distance between real and generated images using the Fréchet Inception Distance (FID) metric. BC subtype and mutation status of driver genes were obtained from the cBioPortal platform, where 3 genes were selected based on the number of mutated patients. A convolutional neural network (CNN) was constructed and trained using the generated MRIs for mutation status prediction. Receiver operating characteristic area under curve (ROC-AUC) and precision-recall area under curve (PR-AUC) were used to evaluate the performance of the CNN models for mutation status prediction. Precision, recall and F1 score were used to evaluate the performance of the CNN model in subtype classification. RESULTS: The FID of the images from the well-trained cGAN model based on the test set is 1.31. The CNN for TP53, PIK3CA, and CDH1 mutation prediction yielded ROC-AUC values 0.9508, 0.7515, and 0.8136 and PR-AUC are 0.9009, 0.7184, and 0.5007, respectively for the three genes. Multi-class subtype prediction achieved precision, recall and F1 scores of 0.8444, 0.8435 and 0.8336 respectively. The source code and related data implemented the algorithms can be found in the project GitHub at https://github.com/mattthuang/BC_RadiogenomicGAN . CONCLUSION: Our study establishes cGAN as a viable tool for generating synthetic BC MRIs for mutation status prediction and subtype classification to better characterize the heterogeneity of BC in patients. The synthetic images also have the potential to significantly augment existing MRI data and circumvent issues surrounding data sharing and patient privacy for future BC machine learning studies.

Novel model integrating computed tomography-based image markers with genetic markers for discriminating radiation pneumonitis in patients with unresectable stage III non-small cell lung cancer receiving radiotherapy: a retrospective multi-center radiogenomics study.

BACKGROUND: Chemoradiotherapy is a critical treatment for patients with locally advanced and unresectable non-small cell lung cancer (NSCLC), and it is essential to identify high-risk patients as early as possible owing to the high incidence of radiation pneumonitis (RP). Increasing attention is being paid to the effects of endogenous factors for RP. This study aimed to investigate the value of computed tomography (CT)-based radiomics combined with genomics in analyzing the risk of grade ≥ 2 RP in unresectable stage III NSCLC. METHODS: In this retrospective multi-center observational study, 100 patients with unresectable stage III NSCLC who were treated with chemoradiotherapy were analyzed. Radiomics features of the entire lung were extracted from pre-radiotherapy CT images. The least absolute shrinkage and selection operator algorithm was used for optimal feature selection to calculate the Rad-score for predicting grade ≥ 2 RP. Genomic DNA was extracted from formalin-fixed paraffin-embedded pretreatment biopsy tissues. Univariate and multivariate logistic regression analyses were performed to identify predictors of RP for model development. The area under the receiver operating characteristic curve was used to evaluate the predictive capacity of the model. Statistical comparisons of the area under the curve values between different models were performed using the DeLong test. Calibration and decision curves were used to demonstrate discriminatory and clinical benefit ratios, respectively. RESULTS: The Rad-score was constructed from nine radiomic features to predict grade ≥ 2 RP. Multivariate analysis demonstrated that histology, Rad-score, and XRCC1 (rs25487) allele mutation were independent high-risk factors correlated with RP. The area under the curve of the integrated model combining clinical factors, radiomics, and genomics was significantly higher than that of any single model (0.827 versus 0.594, 0.738, or 0.641). Calibration and decision curve analyses confirmed the satisfactory clinical feasibility and utility of the nomogram. CONCLUSION: Histology, Rad-score, and XRCC1 (rs25487) allele mutation could predict grade ≥ 2 RP in patients with locally advanced unresectable NSCLC after chemoradiotherapy, and the integrated model combining clinical factors, radiomics, and genomics demonstrated the best predictive efficacy.

Correlation of gene expression with magnetic resonance imaging features of retinoblastoma: a multi-center radiogenomics validation study.

OBJECTIVES: To validate associations between MRI features and gene expression profiles in retinoblastoma, thereby evaluating the repeatability of radiogenomics in retinoblastoma. METHODS: In this retrospective multicenter cohort study, retinoblastoma patients with gene expression data and MRI were included. MRI features (scored blinded for clinical data) and matched genome-wide gene expression data were used to perform radiogenomic analysis. Expression data from each center were first separately processed and analyzed. The end product normalized expression values from different sites were subsequently merged by their Z-score to permit cross-sites validation analysis. The MRI features were non-parametrically correlated with expression of photoreceptorness (radiogenomic analysis), a gene expression signature informing on disease progression. Outcomes were compared to outcomes in a previous described cohort. RESULTS: Thirty-six retinoblastoma patients were included, 15 were female (42%), and mean age was 24 (SD 18) months. Similar to the prior evaluation, this validation study showed that low photoreceptorness gene expression was associated with advanced stage imaging features. Validated imaging features associated with low photoreceptorness were multifocality, a tumor encompassing the entire retina or entire globe, and a diffuse growth pattern (all p < 0.05). There were a number of radiogenomic associations that were also not validated. CONCLUSIONS: A part of the radiogenomic associations could not be validated, underlining the importance of validation studies. Nevertheless, cross-center validation of imaging features associated with photoreceptorness gene expression highlighted the capability radiogenomics to non-invasively inform on molecular subtypes in retinoblastoma. CLINICAL RELEVANCE STATEMENT: Radiogenomics may serve as a surrogate for molecular subtyping based on histopathology material in an era of eye-sparing retinoblastoma treatment strategies. KEY POINTS: • Since retinoblastoma is increasingly treated using eye-sparing methods, MRI features informing on molecular subtypes that do not rely on histopathology material are important. • A part of the associations between retinoblastoma MRI features and gene expression profiles (radiogenomics) were validated. • Radiogenomics could be a non-invasive technique providing information on the molecular make-up of retinoblastoma.

Radiogenomics as an Integrated Approach to Glioblastoma Precision Medicine.

PURPOSE OF REVIEW: Isocitrate dehydrogenase wild-type glioblastoma is the most aggressive primary brain tumour in adults. Its infiltrative nature and heterogeneity confer a dismal prognosis, despite multimodal treatment. Precision medicine is increasingly advocated to improve survival rates in glioblastoma management; however, conventional neuroimaging techniques are insufficient in providing the detail required for accurate diagnosis of this complex condition. RECENT FINDINGS: Advanced magnetic resonance imaging allows more comprehensive understanding of the tumour microenvironment. Combining diffusion and perfusion magnetic resonance imaging to create a multiparametric scan enhances diagnostic power and can overcome the unreliability of tumour characterisation by standard imaging. Recent progress in deep learning algorithms establishes their remarkable ability in image-recognition tasks. Integrating these with multiparametric scans could transform the diagnosis and monitoring of patients by ensuring that the entire tumour is captured. As a corollary, radiomics has emerged as a powerful approach to offer insights into diagnosis, prognosis, treatment, and tumour response through extraction of information from radiological scans, and transformation of these tumour characteristics into quantitative data. Radiogenomics, which links imaging features with genomic profiles, has exhibited its ability in characterising glioblastoma, and determining therapeutic response, with the potential to revolutionise management of glioblastoma. The integration of deep learning algorithms into radiogenomic models has established an automated, highly reproducible means to predict glioblastoma molecular signatures, further aiding prognosis and targeted therapy. However, challenges including lack of large cohorts, absence of standardised guidelines and the 'black-box' nature of deep learning algorithms, must first be overcome before this workflow can be applied in clinical practice.

Radiogenomics uncovers an interplay between angiogenesis and clinical outcomes in bladder cancer.

BACKGROUND: Precision medicine has become a promising clinical treatment strategy for various cancers, including bladder cancer, where angiogenesis plays a critical role in cancer progression. However, the relationship between angiogenesis, immune cell infiltration, clinical outcomes, chemotherapy, and targeted therapy remains unclear. METHODS: We conducted a comprehensive evaluation of angiogenesis-related genes (ARGs) to identify their association with immune cell infiltration, transcription patterns, and clinical outcomes in bladder cancer. An ARG score was constructed to identify angiogenic subgroups in each sample and we evaluated their predictive performance for overall survival rate and treatment response. In addition, we optimized existing clinical detection protocols by performing image data processing. RESULTS: Our study revealed the genomic-level mutant landscape and expression patterns of ARGs in bladder cancer specimens. Using analysis, we identified three molecular subgroups where ARG mutations correlated with patients' pathological features, clinical outcomes, and immune cell infiltration. To facilitate clinical applicability, we constructed a precise nomogram based on the ARG score, which significantly correlated with stem cell index and drug sensitivity. Finally, we proposed the radiogenomics model, which combines the precision of genomics with the convenience of radiomics. CONCLUSION: Our study sheds light on the prognostic characteristics of ARGs in bladder cancer and provides insights into the tumor environment's characteristics to explore more effective immunotherapy strategies. The findings have significant implications for the development of personalized treatment approaches in bladder cancer and pave the way for future studies in this field.

Radiogenomics Map-Based Molecular and Imaging Phenotypical Characterization in Localised Prostate Cancer Using Pre-Biopsy Biparametric MR Imaging.

To create a radiogenomics map and evaluate the correlation between molecular and imaging phenotypes in localized prostate cancer (PCa), using radical prostatectomy histopathology as a reference standard. Radiomic features were extracted from T2-weighted (T2WI) and Apparent Diffusion Coefficient (ADC) images of clinically localized PCa patients (n = 15) across different Gleason score-based risk categories. DNA extraction was performed on formalin-fixed, paraffin-embedded (FFPE) samples. Gene expression analysis of androgen receptor expression, apoptosis, and hypoxia was conducted using the Chromosome Analysis Suite (ChAS) application and OSCHIP files. The relationship between gene expression alterations and textural features was assessed using Pearson's correlation analysis. Receiver operating characteristic (ROC) analysis was utilized to evaluate the predictive accuracy of the model. A significant correlation was observed between radiomic texture features and copy number variation (CNV) of genes associated with apoptosis, hypoxia, and androgen receptor (p-value ≤ 0.05). The identified radiomic features, including Sum Entropy ADC, Inverse Difference ADC, Sum Variance T2WI, Entropy T2WI, Difference Variance T2WI, and Angular Secondary Moment T2WI, exhibited potential for predicting cancer grade and biological processes such as apoptosis and hypoxia. Incorporating radiomics and genomics into a prediction model significantly improved the prediction of prostate cancer grade (clinically significant prostate cancer), yielding an AUC of 0.95. Radiomic texture features significantly correlate with genotypes for apoptosis, hypoxia, and androgen receptor expression in localised prostate cancer. Integration of these into the prediction model improved prediction accuracy of clinically significant prostate cancer.

Graph Neural Network Model for Prediction of Non-Small Cell Lung Cancer Lymph Node Metastasis Using Protein-Protein Interaction Network and 18F-FDG PET/CT Radiomics.

The image texture features obtained from 18F-fluorodeoxyglucose positron emission tomography/computed tomography (18F-FDG PET/CT) images of non-small cell lung cancer (NSCLC) have revealed tumor heterogeneity. A combination of genomic data and radiomics may improve the prediction of tumor prognosis. This study aimed to predict NSCLC metastasis using a graph neural network (GNN) obtained by combining a protein-protein interaction (PPI) network based on gene expression data and image texture features. 18F-FDG PET/CT images and RNA sequencing data of 93 patients with NSCLC were acquired from The Cancer Imaging Archive. Image texture features were extracted from 18F-FDG PET/CT images and area under the curve receiver operating characteristic curve (AUC) of each image feature was calculated. Weighted gene co-expression network analysis (WGCNA) was used to construct gene modules, followed by functional enrichment analysis and identification of differentially expressed genes. The PPI of each gene module and genes belonging to metastasis-related processes were converted via a graph attention network. Images and genomic features were concatenated. The GNN model using PPI modules from WGCNA and metastasis-related functions combined with image texture features was evaluated quantitatively. Fifty-five image texture features were extracted from 18F-FDG PET/CT, and radiomic features were selected based on AUC (n = 10). Eighty-six gene modules were clustered by WGCNA. Genes (n = 19) enriched in the metastasis-related pathways were filtered using DEG analysis. The accuracy of the PPI network, derived from WGCNA modules and metastasis-related genes, improved from 0.4795 to 0.5830 (p < 2.75 × 10-12). Integrating PPI of four metastasis-related genes with 18F-FDG PET/CT image features in a GNN model elevated its accuracy over a without image feature model to 0.8545 (95% CI = 0.8401-0.8689, p-value < 0.02). This model demonstrated significant enhancement compared to the model using PPI and 18F-FDG PET/CT derived from WGCNA (p-value < 0.02), underscoring the critical role of metastasis-related genes in prediction model. The enhanced predictive capability of the lymph node metastasis prediction GNN model for NSCLC, achieved through the integration of comprehensive image features with genomic data, demonstrates promise for clinical implementation.

Integrated imaging and molecular analysis to decipher tumor microenvironment in the era of immunotherapy.

Radiological imaging is an integral component of cancer care, including diagnosis, staging, and treatment response monitoring. It contains rich information about tumor phenotypes that are governed not only by cancer cellintrinsic biological processes but also by the tumor microenvironment, such as the composition and function of tumor-infiltrating immune cells. By analyzing the radiological scans using a quantitative radiomics approach, robust relations between specific imaging and molecular phenotypes can be established. Indeed, a number of studies have demonstrated the feasibility of radiogenomics for predicting intrinsic molecular subtypes and gene expression signatures in breast cancer based on MRI. In parallel, promising results have been shown for inferring the amount of tumor-infiltrating lymphocytes, a key factor for the efficacy of cancer immunotherapy, from standard-of-care radiological images. Compared with the biopsy-based approach, radiogenomics offers a unique avenue to profile the molecular makeup of the tumor and immune microenvironment as well as its evolution in a noninvasive and holistic manner through longitudinal imaging scans. Here, we provide a systematic review of the state of the art radiogenomics studies in the era of immunotherapy and discuss emerging paradigms and opportunities in AI and deep learning approaches. These technical advances are expected to transform the radiogenomics field, leading to the discovery of reliable imaging biomarkers. This will pave the way for their clinical translation to guide precision cancer therapy.

Deep learning, radiomics and radiogenomics applications in the digital breast tomosynthesis: a systematic review.

BACKGROUND: Recent advancements in computing power and state-of-the-art algorithms have helped in more accessible and accurate diagnosis of numerous diseases. In addition, the development of de novo areas in imaging science, such as radiomics and radiogenomics, have been adding more to personalize healthcare to stratify patients better. These techniques associate imaging phenotypes with the related disease genes. Various imaging modalities have been used for years to diagnose breast cancer. Nonetheless, digital breast tomosynthesis (DBT), a state-of-the-art technique, has produced promising results comparatively. DBT, a 3D mammography, is replacing conventional 2D mammography rapidly. This technological advancement is key to AI algorithms for accurately interpreting medical images. OBJECTIVE AND METHODS: This paper presents a comprehensive review of deep learning (DL), radiomics and radiogenomics in breast image analysis. This review focuses on DBT, its extracted synthetic mammography (SM), and full-field digital mammography (FFDM). Furthermore, this survey provides systematic knowledge about DL, radiomics, and radiogenomics for beginners and advanced-level researchers. RESULTS: A total of 500 articles were identified, with 30 studies included as the set criteria. Parallel benchmarking of radiomics, radiogenomics, and DL models applied to the DBT images could allow clinicians and researchers alike to have greater awareness as they consider clinical deployment or development of new models. This review provides a comprehensive guide to understanding the current state of early breast cancer detection using DBT images. CONCLUSION: Using this survey, investigators with various backgrounds can easily seek interdisciplinary science and new DL, radiomics, and radiogenomics directions towards DBT.

MRI-based breast cancer radiogenomics using RNA profiling: association with subtypes in a single-center prospective study.

BACKGROUND: There are few prospective studies on the correlations between MRI features and whole RNA-sequencing data in breast cancer according to molecular subtypes. The purpose of our study was to explore the association between genetic profiles and MRI phenotypes of breast cancer and to identify imaging markers that influences the prognosis and treatment according to subtypes. METHODS: From June 2017 to August 2018, MRIs of 95 women with invasive breast cancer were prospectively analyzed, using the breast imaging-reporting and data system and texture analysis. Whole RNA obtained from surgical specimens was analyzed using next-generation sequencing. The association between MRI features and gene expression profiles was analyzed in the entire tumor and subtypes. Gene networks, enriched functions, and canonical pathways were analyzed using Ingenuity Pathway Analysis. The P value for differential expression was obtained using a parametric F test comparing nested linear models and adjusted for multiple testing by reporting Q value. RESULTS: In 95 participants (mean age, 53 years ± 11 [standard deviation]), mass lesion type was associated with upregulation of CCL3L1 (sevenfold) and irregular mass shape was associated with downregulation of MIR421 (sixfold). In estrogen receptor-positive cancer with mass lesion type, CCL3L1 (21-fold), SNHG12 (11-fold), and MIR206 (sevenfold) were upregulated, and MIR597 (265-fold), MIR126 (12-fold), and SOX17 (fivefold) were downregulated. In triple-negative breast cancer with increased standard deviation of texture analysis on precontrast T1-weighted imaging, CLEC3A (23-fold), SRGN (13-fold), HSPG2 (sevenfold), KMT2D (fivefold), and VMP1 (fivefold) were upregulated, and IGLC2 (73-fold) and PRDX4 (sevenfold) were downregulated (all, P < 0.05 and Q < 0.1). Gene network and functional analysis showed that mass type estrogen receptor-positive cancers were associated with cell growth, anti-estrogen resistance, and poor survival. CONCLUSION: MRI characteristics are associated with the different expressions of genes related to metastasis, anti-drug resistance, and prognosis, depending on the molecular subtypes of breast cancer.

Radiogenomic analysis of cellular tumor-stroma heterogeneity as a prognostic predictor in breast cancer.

BACKGROUND: The tumor microenvironment and intercellular communication between solid tumors and the surrounding stroma play crucial roles in cancer initiation, progression, and prognosis. Radiomics provides clinically relevant information from radiological images; however, its biological implications in uncovering tumor pathophysiology driven by cellular heterogeneity between the tumor and stroma are largely unknown. We aimed to identify radiogenomic signatures of cellular tumor-stroma heterogeneity (TSH) to improve breast cancer management and prognosis analysis. METHODS: This retrospective multicohort study included five datasets. Cell subpopulations were estimated using bulk gene expression data, and the relative difference in cell subpopulations between the tumor and stroma was used as a biomarker to categorize patients into good- and poor-survival groups. A radiogenomic signature-based model utilizing dynamic contrast-enhanced magnetic resonance imaging (DCE-MRI) was developed to target TSH, and its clinical significance in relation to survival outcomes was independently validated. RESULTS: The final cohorts of 1330 women were included for cellular TSH biomarker identification (n = 112, mean age, 57.3 years ± 14.6) and validation (n = 886, mean age, 58.9 years ± 13.1), radiogenomic signature of TSH identification (n = 91, mean age, 55.5 years ± 11.4), and prognostic (n = 241) assessments. The cytotoxic lymphocyte biomarker differentiated patients into good- and poor-survival groups (p < 0.0001) and was independently validated (p = 0.014). The good survival group exhibited denser cell interconnections. The radiogenomic signature of TSH was identified and showed a positive association with overall survival (p = 0.038) and recurrence-free survival (p = 3 × 10-4). CONCLUSION: Radiogenomic signatures provide insights into prognostic factors that reflect the imbalanced tumor-stroma environment, thereby presenting breast cancer-specific biological implications and prognostic significance.

Computed tomography-detected extramural venous invasion-related gene signature: a potential negative biomarker of immune checkpoint inhibitor treatment in patients with gastric cancer.

BACKGROUND: To investigate the association between computed tomography (CT)-detected extramural venous invasion (EMVI)-related genes and immunotherapy resistance and immune escape in patients with gastric cancer (GC). METHODS: Thirteen patients with pathologically proven locally advanced GC who had undergone preoperative abdominal contrast-enhanced CT and radical resection surgery were included in this study. Transcriptome sequencing was multidetector performed on the cancerous tissue obtained during surgery, and EMVI-related genes (P value for association < 0.001) were selected. A single-sample gene set enrichment analysis algorithm was also used to divide all GC samples (n = 377) in The Cancer Genome Atlas (TCGA) database into high and low EMVI-immune related groups based on immune-related differential genes. Cluster analysis was used to classify EMVI-immune-related genotypes, and survival among patients was validated in TCGA and Gene Expression Omnibus (GEO) cohorts. The EMVI scores were calculated using principal component analysis (PCA), and GC samples were divided into high and low EMVI score groups. Microsatellite instability (MSI) status, tumor mutation burden (TMB), response rate to immune checkpoint inhibitors (ICIs), immune escape were compared between the high and low EMVI score groups. Hub gene of the model in pan-cancer analysis was also performed. RESULTS: There were 17 EMVI-immune-related genes used for cluster analysis. PCA identified 8 genes (PCH17, SEMA6B, GJA4, CD34, ACVRL1, SOX17, CXCL12, DYSF) that were used to calculate EMVI scores. High EMVI score groups had lower MSI, TMB and response rate of ICIs, status but higher immune escape status. Among the 8 genes used for EMVI scores, CXCL12 and SOX17 were at the core of the protein-protein interaction (PPI) network and had a higher priority in pan-cancer analysis. Immunohistochemical analysis showed that the expression of CXCL12 and SOX17 was significantly higher in CT-detected EMVI-positive samples than in EMVI-negative samples (P < 0.0001). CONCLUSION: A CT-detected EMVI gene signature could be a potential negative biomarker for ICIs treatment, as the signature is negatively correlated with TMB, and MSI, resulting in poorer prognosis.

Radiogenomic analysis of prediction HER2 status in breast cancer by linking ultrasound radiomic feature module with biological functions.

BACKGROUND: Human epidermal growth factor receptor 2 (HER2) overexpressed associated with poor prognosis in breast cancer and HER2 has been defined as a therapeutic target for breast cancer treatment. We aimed to explore the molecular biological information in ultrasound radiomic features (URFs) of HER2-positive breast cancer using radiogenomic analysis. Moreover, a radiomics model was developed to predict the status of HER2 in breast cancer. METHODS: This retrospective study included 489 patients who were diagnosed with breast cancer. URFs were extracted from a radiomics analysis set using PyRadiomics. The correlations between differential URFs and HER2-related genes were calculated using Pearson correlation analysis. Functional enrichment of the identified URFs-correlated HER2 positive-specific genes was performed. Lastly, the radiomics model was developed based on the URF-module mined from auxiliary differential URFs to assess the HER2 status of breast cancer. RESULTS: Eight differential URFs (p < 0.05) were identified among the 86 URFs extracted by Pyradiomics. 25 genes that were found to be the most closely associated with URFs. Then, the relevant biological functions of each differential URF were obtained through functional enrichment analysis. Among them, Zone Entropy is related to immune cell activity, which regulate the generation of calcification in breast cancer. The radiomics model based on the Logistic classifier and URF-module showed good discriminative ability (AUC = 0.80, 95% CI). CONCLUSION: We searched for the URFs of HER2-positive breast cancer, and explored the underlying genes and biological functions of these URFs. Furthermore, the radiomics model based on the Logistic classifier and URF-module relatively accurately predicted the HER2 status in breast cancer.

Radiogenomic Associations Clear Cell Renal Cell Carcinoma: An Exploratory Study.

INTRODUCTION: This study investigates how quantitative texture analysis can be used to non-invasively identify novel radiogenomic correlations with clear cell renal cell carcinoma (ccRCC) biomarkers. METHODS: The Cancer Genome Atlas-Kidney Renal Clear Cell Carcinoma open-source database was used to identify 190 sets of patient genomic data that had corresponding multiphase contrast-enhanced CT images in The Cancer Imaging Archive. 2,824 radiomic features spanning fifteen texture families were extracted from CT images using a custom-built MATLAB software package. Robust radiomic features with strong inter-scanner reproducibility were selected. Random forest, AdaBoost, and elastic net machine learning (ML) algorithms evaluated the ability of the selected radiomic features to predict the presence of 12 clinically relevant molecular biomarkers identified from the literature. ML analysis was repeated with cases stratified by stage (I/II vs. III/IV) and grade (1/2 vs. 3/4). 10-fold cross validation was used to evaluate model performance. RESULTS: Before stratification by tumor grade and stage, radiomics predicted the presence of several biomarkers with weak discrimination (AUC 0.60-0.68). Once stratified, radiomics predicted KDM5C, SETD2, PBRM1, and mTOR mutation status with acceptable to excellent predictive discrimination (AUC ranges from 0.70 to 0.86). CONCLUSIONS: Radiomic texture analysis can potentially identify a variety of clinically relevant biomarkers in patients with ccRCC and may have a prognostic implication.

Construction of a radiogenomic association map of pancreatic ductal adenocarcinoma.

BACKGROUND: Pancreatic adenocarcinoma (PDAC) persists as a malignancy with high morbidity and mortality that can benefit from new means to characterize and detect these tumors, such as radiogenomics. In order to address this gap in the literature, constructed a transcriptomic-CT radiogenomic (RG) map for PDAC. METHODS: In this Institutional Review Board approved study, a cohort of subjects (n = 50) with gene expression profile data paired with histopathologically confirmed resectable or borderline resectable PDAC were identified. Studies with pre-operative contrast-enhanced CT images were independently assessed for a set of 88 predefined imaging features. Microarray gene expression profiling was then carried out on the histopathologically confirmed pancreatic adenocarcinomas and gene networks were constructed using Weighted Gene Correlation Network Analysis (WCGNA) (n = 37). Data were analyzed with bioinformatics analyses, multivariate regression-based methods, and Kaplan-Meier survival analyses. RESULTS: Survival analyses identified multiple features of interest that were significantly associated with overall survival, including Tumor Height (P = 0.014), Tumor Contour (P = 0.033), Tumor-stroma Interface (P = 0.014), and the Tumor Enhancement Ratio (P = 0.047). Gene networks for these imaging features were then constructed using WCGNA and further annotated according to the Gene Ontology (GO) annotation framework for a biologically coherent interpretation of the imaging trait-associated gene networks, ultimately resulting in a PDAC RG CT-transcriptome map composed of 3 stage-independent imaging traits enriched in metabolic processes, telomerase activity, and podosome assembly (P < 0.05). CONCLUSIONS: A CT-transcriptomic RG map for PDAC composed of semantic and quantitative traits with associated biology processes predictive of overall survival, was constructed, that serves as a reference for further mechanistic studies for non-invasive phenotyping of pancreatic tumors.

Radiomics for the non-invasive prediction of PD-L1 expression in patients with brain metastases secondary to non-small cell lung cancer.

BACKGROUND: The expression level of the programmed cell death ligand 1 (PD-L1) appears to be a predictor for response to immunotherapy using checkpoint inhibitors in patients with non-small cell lung cancer (NSCLC). As differences in terms of PD-L1 expression levels in the extracranial primary tumor and the brain metastases may occur, a reliable method for the non-invasive assessment of the intracranial PD-L1 expression is, therefore of clinical value. Here, we evaluated the potential of radiomics for a non-invasive prediction of PD-L1 expression in patients with brain metastases secondary to NSCLC. PATIENTS AND METHODS: Fifty-three NSCLC patients with brain metastases from two academic neuro-oncological centers (group 1, n = 36 patients; group 2, n = 17 patients) underwent tumor resection with a subsequent immunohistochemical evaluation of the PD-L1 expression. Brain metastases were manually segmented on preoperative T1-weighted contrast-enhanced MRI. Group 1 was used for model training and validation, group 2 for model testing. After image pre-processing and radiomics feature extraction, a test-retest analysis was performed to identify robust features prior to feature selection. The radiomics model was trained and validated using random stratified cross-validation. Finally, the best-performing radiomics model was applied to the test data. Diagnostic performance was evaluated using receiver operating characteristic (ROC) analyses. RESULTS: An intracranial PD-L1 expression (i.e., staining of at least 1% or more of tumor cells) was present in 18 of 36 patients (50%) in group 1, and 7 of 17 patients (41%) in group 2. Univariate analysis identified the contrast-enhancing tumor volume as a significant predictor for PD-L1 expression (area under the ROC curve (AUC), 0.77). A random forest classifier using a four-parameter radiomics signature, including tumor volume, yielded an AUC of 0.83 ± 0.18 in the training data (group 1), and an AUC of 0.84 in the external test data (group 2). CONCLUSION: The developed radiomics classifiers allows for a non-invasive assessment of the intracranial PD-L1 expression in patients with brain metastases secondary to NSCLC with high accuracy.

Integrative Bayesian models using Post-selective inference: A case study in radiogenomics.

Integrative analyses based on statistically relevant associations between genomics and a wealth of intermediary phenotypes (such as imaging) provide vital insights into their clinical relevance in terms of the disease mechanisms. Estimates for uncertainty in the resulting integrative models are however unreliable unless inference accounts for the selection of these associations with accuracy. In this paper, we develop selection-aware Bayesian methods, which (1) counteract the impact of model selection bias through a "selection-aware posterior" in a flexible class of integrative Bayesian models post a selection of promising variables via ℓ1 -regularized algorithms; (2) strike an inevitable trade-off between the quality of model selection and inferential power when the same data set is used for both selection and uncertainty estimation. Central to our methodological development, a carefully constructed conditional likelihood function deployed with a reparameterization mapping provides tractable updates when gradient-based Markov chain Monte Carlo (MCMC) sampling is used for estimating uncertainties from the selection-aware posterior. Applying our methods to a radiogenomic analysis, we successfully recover several important gene pathways and estimate uncertainties for their associations with patient survival times.