Search

In this work we discuss a spatial evolutionary model for a heterogeneous cancer cell population. We consider the gain-of-function mutations that not only change the fitness potential of the mutant phenotypes against normal background cells but may also increase the relative motility of the mutant cells. The spatial modeling is implemented as a stochastic evolutionary system on a structured grid (a lattice, with random neighborhoods, which is not necessarily bi-directional) or on a two-dimensional unstructured mesh, i.e. a bi-directional graph with random numbers of neighbors. We present a computational approach to investigate the fixation probability of mutants in these spatial models. Additionally, we examine the effect of the migration potential on the spatial dynamics of mutants on unstructured meshes. Our results suggest that the probability of fixation is negatively correlated with the width of the distribution of the neighborhood size. Also, the fixation probability increases given a migration potential for mutants. We find that the fixation probability (of advantaged, disadvantaged and neutral mutants) on unstructured meshes is relatively smaller than the corresponding results on regular grids. More importantly, in the case of neutral mutants the introduction of a migration potential has a critical effect on the fixation probability and increases this by orders of magnitude. Further, we examine the effect of boundaries and as intuitively expected, the fixation probability is smaller on the boundary of regular grids when compared to its value in the bulk. Based on these computational results, we speculate on possible better therapeutic strategies that may delay tumor progression to some extent.

Objectives: Radiomics can predict patient outcomes by automatically extracting a large number of features from medical images. This study is aimed to investigate the sensitivity of radiomics features extracted from 2 different pipelines, namely, Pyradiomics and RaCat, as well as the impact of gray-level discretization on the discovery of immune checkpoint inhibitors (ICIs) biomarkers., Methods: A retrospective cohort of 164 non-small cell lung cancer patients administered with ICIs was used in this study. Radiomic features were extracted from the pre-treatment CT scans. Univariate models were used to assess the association of common radiomics features between 2 libraries with progression-free survival (PFS), programmed death ligand 1 (PD-L1), and tumour infiltrating lymphocytes (CD8 counts). In addition, we also examined the impact of gray-level discretization incorporated in Pyradiomics on the robustness of features across various clinical endpoints., Results: We extracted 1224, 441 radiomic features using Pyradiomics and RaCat, respectively. Among these, 75 features were found to be common between the 2 libraries. Our analysis revealed that the directionality of association between radiomic features and clinical endpoints is highly dependent on the library. Notably, a larger number of Pyradiomics features were statistically associated with PFS, whereas RaCat features showed a stronger association with PD-L1 expression. Furthermore, intensity-based features were found to have a consistent association with clinical endpoints regardless of the gray-level discretization parameters in Pyradiomics-extracted features., Conclusions: This study highlights the heterogeneity of radiomics libraries and the gray-level discretization parameters that will impact the feature selection and predictive model development for biomarkers. Importantly, our work highlights the significance of standardizing radiomic features to facilitate translational studies that use imaging as an endpoint., Advances in Knowledge: Our study emphasizes the need to select stable CT-derived handcrafted features to build immunotherapy biomarkers, which is a necessary precursor for multi-institutional validation of imaging biomarkers., (© The Author(s) 2024. Published by Oxford University Press on behalf of the British Institute of Radiology. All rights reserved. For permissions, please email: journals.permissions@oup.com.)

Objective: Lung cancer remains the leading cause of cancer-related mortality worldwide, with most cases diagnosed at advanced stages. Hence, there is a need to develop effective predictive models for early detection. This study aims to investigate the impact of imaging parameters and delta radiomic features from temporal scans on lung cancer risk prediction., Methods: Using the National Lung Screening Trial (NLST) within a nested case-control study involving 462 positive screenings, radiomic features were extracted from temporal computed tomography (CT) scans and harmonized with ComBat method to adjust variations in slice thickness category (TC) and reconstruction kernel type (KT). Both harmonized and non-harmonized features from baseline (T0), delta features between T0 and a year later (T1), and combined T0 and delta features were utilized for the analysis. Feature reduction was done using LASSO, followed by five feature selection (FS) methods and nine machine learning (ML) models, evaluated with 5-fold cross-validation repeated 10 times. Synthetic Minority Oversampling Technique (SMOTE) was applied to address class imbalances for lung cancer risk prediction., Results: Models using delta features outperformed baseline features, with SMOTE consistently boosting performance when using combination of baseline and delta features. TC-based harmonized features improved performance with SMOTE, but overall, harmonization did not significantly enhance the model performance. The highest test score of 0.76 was achieved in three scenarios: delta features with a Gradient Boosting (GB) model (TC-based harmonization and MultiSurf FS); and T0 + delta features, with both a Support Vector Classifier (SVC) model (KT-based harmonization and F-test FS), and an XGBoost (XGB) model (TC-based harmonization and Mutual Information (MI) FS), all using SMOTE., Conclusions: This study underscores the significance of delta radiomic features and balanced datasets to improve lung cancer prediction. While our findings are based on a subsample of NLST data, they provide a valuable foundation for further exploration. Further research is needed to assess the impact of harmonization on imaging-derived models. Future investigations should explore advanced harmonization techniques and additional imaging parameters to develop robust radiomics-based biomarkers of lung cancer risk., (© 2024 The Author(s). Cancer Medicine published by John Wiley & Sons Ltd.)

Background/objectives: Pulmonary neuroendocrine neoplasms (NENs) account for 20% of malignant lung tumors. Their management is challenging due to their diverse clinical features and aggressive nature. Currently, metabolomics offers a range of potential cancer biomarkers for diagnosis, monitoring tumor progression, and assessing therapeutic response. However, a specific metabolomic profile for early diagnosis of lung NENs has yet to be identified. This study aims to identify specific metabolomic profiles that can serve as biomarkers for early diagnosis of lung NENs., Methods: We measured 153 metabolites using liquid chromatography combined with mass spectrometry (LC-MS) in the plasma of 120 NEN patients and compared them with those of 71 healthy individuals. Additionally, we compared these profiles with those of 466 patients with non-small-cell lung cancers (NSCLCs) to ensure clinical relevance., Results: We identified 21 metabolites with consistently altered plasma concentrations in NENs. Compared to healthy controls, 18 metabolites were specific to carcinoid tumors, 5 to small-cell lung carcinomas (SCLCs), and 10 to large-cell neuroendocrine carcinomas (LCNECs). These findings revealed alterations in various metabolic pathways, such as fatty acid biosynthesis and beta-oxidation, the Warburg effect, and the citric acid cycle., Conclusions: Our study identified biomarker metabolites in the plasma of patients with each subtype of lung NENs and demonstrated significant alterations in several metabolic pathways. These metabolomic profiles could potentially serve as biomarkers for early diagnosis and better management of lung NENs.

Background: Immune checkpoint inhibitors (ICIs) have revolutionized non-small cell lung cancers (NSCLCs) treatment, but only 20-30% of patients benefit from these treatments. Currently, PD-L1 expression in tumor cells is the only clinically approved predictor of ICI response in lung cancer, but concerns arise due to its low negative and positive predictive value. Recent studies suggest that CXCL13+ T cells in the tumor microenvironment (TME) may be a good predictor of response. We aimed to assess if CXCL13+ cell localization within the TME can predict ICI response in advanced NSCLC patients. Methods: This retrospective study included 65 advanced NSCLC patients treated with Nivolumab/Pembrolizumab at IUCPQ or CHUM and for whom a pretreatment surgical specimen was available. Good responders were defined as having a complete radiologic response at 1 year, and bad responders were defined as showing cancer progression at 1 year. IHC staining for CXCL13 was carried out on a representative slide from a resection specimen, and CXCL13+ cell density was evaluated in tumor (T), invasive margin (IM), non-tumor (NT), and tertiary lymphoid structure (TLS) compartments. Cox models were used to analyze progression-free survival (PFS) and overall survival (OS) probability, while the Mann-Whitney test was used to compare CXCL13+ cell density between responders and non-responders. Results: We showed that CXCL13+ cell density localization within the TME is associated with ICI efficacy. An increased density of CXCL13+ cells across all compartments was associated with a poorer prognostic (OS; HR = 1.22; 95%CI = 1.04-1.42; p = 0.01, PFS; HR = 1.16; p = 0.02), or a better prognostic when colocalized within TLSs (PFS; HR = 0.84, p = 0.03). Conclusion: Our results support the role of CXCL13+ cells in advanced NSCLC patients, with favorable prognosis when localized within TLSs and unfavorable prognosis when present elsewhere. The concomitant proximity of CXCL13+ and CD20+ cells within TLSs may favor antigen presentation to T cells, thus enhancing the effect of PD-1/PD-L1 axis inhibition. Further validation is warranted to confirm the potential relevance of this biomarker in a clinical setting.

Background: Recent advances in cancer biomarker development have led to a surge of distinct data modalities, such as medical imaging and histopathology. To develop predictive immunotherapy biomarkers, these modalities are leveraged independently, despite their orthogonality. This study aims to explore the cross-scale association between radiological scans and digitalized pathology images for immunotherapy-treated non-small cell lung cancer (NSCLC) patients., Methods: This study involves 36 NSCLC patients who were treated with immunotherapy and for whom both radiology and pathology images were available. A total of 851 and 260 features were extracted from CT scans and cell density maps of histology images at different resolutions. We investigated the radiopathomics relationship and their association with clinical and biological endpoints. We used the Kolmogorov-Smirnov (KS) method to test the differences between the distributions of correlation coefficients with the two imaging modality features. Unsupervised clustering was done to identify which imaging modality captures poor and good survival patients., Results: Our results demonstrated a significant correlation between cell density pathomics and radiomics features. Furthermore, we also found a varying distribution of correlation values between imaging-derived features and clinical endpoints. The KS test revealed that the two imaging feature distributions were different for PFS and CD8 counts, while similar for OS. In addition, clustering analysis resulted in significant differences in the two clusters generated from the radiomics and pathomics features with respect to patient survival and CD8 counts., Conclusion: The results of this study suggest a cross-scale association between CT scans and pathology H&E slides among ICI-treated patients. These relationships can be further explored to develop multimodal immunotherapy biomarkers to advance personalized lung cancer care.

Background: Although substantial efforts have been made to build molecular biomarkers to predict radiation sensitivity, the ability to accurately stratify the patients is still limited. In this study, we aim to leverage large-scale radiogenomics datasets to build genomic predictors of radiation response using the integral of the radiation dose-response curve., Methods: Two radiogenomics datasets consisting of 511 and 60 cancer cell lines were utilized to develop genomic predictors of radiation sensitivity. The intrinsic radiation sensitivity, defined as the integral of the dose-response curve (AUC) was used as the radioresponse variable. The biological determinants driving AUC and SF2 were compared using pathway analysis. To build the predictive model, the largest and smallest datasets consisting of 511 and 60 cancer cell lines were used as the discovery and validation cohorts, respectively, with AUC as the response variable., Results: Utilizing a compendium of three pathway databases, we illustrated that integral of the radiobiological model provides a more comprehensive characterization of molecular processes underpinning radioresponse compared to SF2. Furthermore, more pathways were found to be unique to AUC than SF2-30, 288 and 38 in KEGG, REACTOME and WIKIPATHWAYS, respectively. Also, the leading-edge genes driving the biological pathways using AUC were unique and different compared to SF2. With regards to radiation sensitivity gene signature, we obtained a concordance index of 0.65 and 0.61 on the discovery and validation cohorts, respectively., Conclusion: We developed an integrated framework that quantifies the impact of physical radiation dose and the biological effect of radiation therapy in interventional pre-clinical model systems. With the availability of more data in the future, the clinical potential of this signature can be assessed, which will eventually provide a framework to integrate genomics into biologically-driven precision radiation oncology., (© 2023. The Author(s).)

Background: Although the immune checkpoint inhibitors, nivolumab and pembrolizumab, were found to be promising in patients with advanced NSCLC, some of them either do not respond or have recurrence after an initial response. It is still unclear who will benefit from these therapies, and, hence, there is an unmet clinical need to build robust biomarkers., Methods: Patients with advanced NSCLC (N = 323) who were treated with pembrolizumab or nivolumab were retrospectively identified from two institutions. Radiomics features extracted from baseline pretreatment computed tomography scans along with the clinical variables were used to build the predictive models for overall survival (OS), progression-free survival (PFS), and programmed death-ligand 1 (PD-L1). To develop the imaging and integrative clinical-imaging predictive models, we used the XGBoost learning algorithm with ReliefF feature selection method and validated them in an independent cohort. The concordance index for OS, PFS, and area under the curve for PD-L1 was used to evaluate model performance., Results: We developed radiomics and the ensemble radiomics-clinical predictive models for OS, PFS, and PD-L1 expression. The concordance indices of the radiomics model were 0.60 and 0.61 for predicting OS and PFS and area under the curve was 0.61 for predicting PD-L1 in the validation cohort, respectively. The combined radiomics-clinical model resulted in higher performance with 0.65, 0.63, and 0.68 to predict OS, PFS, and PD-L1 in the validation cohort, respectively., Conclusions: We found that pretreatment computed tomography imaging along with clinical data can aid as predictive biomarkers for PD-L1 and survival end points. These imaging-driven approaches may prove useful to expand the therapeutic options for nonresponders and improve the selection of patients who would benefit from immune checkpoint inhibitors., (© 2023 Published by the Elsevier Inc. on behalf of the International Association for the Study of Lung Cancer.)

Background: Immune checkpoint inhibitors (ICIs) are a great breakthrough in cancer treatments and provide improved long-term survival in a subset of non-small cell lung cancer (NSCLC) patients. However, prognostic and predictive biomarkers of immunotherapy still remain an unmet clinical need. In this work, we aim to leverage imaging data and clinical variables to develop survival risk models among advanced NSCLC patients treated with immunotherapy., Methods: This retrospective study includes a total of 385 patients from two institutions who were treated with ICIs. Radiomics features extracted from pretreatment CT scans were used to build predictive models. The objectives were to predict overall survival (OS) along with building a classifier for short- and long-term survival groups. We employed the XGBoost learning method to build radiomics and integrated clinical-radiomics predictive models. Feature selection and model building were developed and validated on a multicenter cohort., Results: We developed parsimonious models that were associated with OS and a classifier for short- and long-term survivor groups. The concordance indices (C-index) of the radiomics model were 0.61 and 0.57 to predict OS in the discovery and validation cohorts, respectively. While the area under the curve (AUC) values of the radiomic models for short- and long-term groups were found to be 0.65 and 0.58 in the discovery and validation cohorts. The accuracy of the combined radiomics-clinical model resulted in 0.63 and 0.62 to predict OS and in 0.77 and 0.62 to classify the survival groups in the discovery and validation cohorts, respectively., Conclusions: We developed and validated novel radiomics and integrated radiomics-clinical survival models among NSCLC patients treated with ICIs. This model has important translational implications, which can be used to identify a subset of patients who are not likely to benefit from immunotherapy. The developed imaging biomarkers may allow early prediction of low-group survivors, though additional validation of these radiomics models is warranted.

Background: Recent developments in artificial intelligence suggest that radiomics may represent a promising non-invasive biomarker to predict response to immune checkpoint inhibitors (ICIs). Nevertheless, validation of radiomics algorithms in independent cohorts remains a challenge due to variations in image acquisition and reconstruction. Using radiomics, we investigated the importance of scan normalization as part of a broader machine learning framework to enable model external generalizability to predict ICI response in non-small cell lung cancer (NSCLC) patients across different centers., Methods: Radiomics features were extracted and compared from 642 advanced NSCLC patients on pre-ICI scans using established open-source PyRadiomics and a proprietary DeepRadiomics deep learning technology. The population was separated into two groups: a discovery cohort of 512 NSCLC patients from three academic centers and a validation cohort that included 130 NSCLC patients from a fourth center. We harmonized images to account for variations in reconstruction kernel, slice thicknesses, and device manufacturers. Multivariable models, evaluated using cross-validation, were used to estimate the predictive value of clinical variables, PD-L1 expression, and PyRadiomics or DeepRadiomics for progression-free survival at 6 months (PFS-6)., Results: The best prognostic factor for PFS-6, excluding radiomics features, was obtained with the combination of Clinical + PD-L1 expression (AUC = 0.66 in the discovery and 0.62 in the validation cohort). Without image harmonization, combining Clinical + PyRadiomics or DeepRadiomics delivered an AUC = 0.69 and 0.69, respectively, in the discovery cohort, but dropped to 0.57 and 0.52, in the validation cohort. This lack of generalizability was consistent with observations in principal component analysis clustered by CT scan parameters. Subsequently, image harmonization eliminated these clusters. The combination of Clinical + DeepRadiomics reached an AUC = 0.67 and 0.63 in the discovery and validation cohort, respectively. Conversely, the combination of Clinical + PyRadiomics failed generalizability validations, with AUC = 0.66 and 0.59., Conclusion: We demonstrated that a risk prediction model combining Clinical + DeepRadiomics was generalizable following CT scan harmonization and machine learning generalization methods. These results had similar performances to routine oncology practice using Clinical + PD-L1. This study supports the strong potential of radiomics as a future non-invasive strategy to predict ICI response in advanced NSCLC., Competing Interests: Author BR reports research funding with the company Imagia. Authors KP, CL-K, FD, LJ, FR, KK, DH, and FIC were Imagia employees at the time of the study. The remaining authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest., (Copyright © 2023 Tonneau, Phan, Manem, Low-Kam, Dutil, Kazandjian, Vanderweyen, Panasci, Malo, Coulombe, Gagné, Elkrief, Belkaïd, Di Jorio, Orain, Bouchard, Muanza, Rybicki, Kafi, Huntsman, Joubert, Chandelier and Routy.)

With the increasing use of immune checkpoint inhibitors (ICIs), there is an urgent need to identify biomarkers to stratify responders and non-responders using programmed death-ligand (PD-L1) expression, and to predict patient-specific outcomes such as progression free survival (PFS). The current study is aimed to determine the feasibility of building imaging-based predictive biomarkers for PD-L1 and PFS through systematically evaluating a combination of several machine learning algorithms with different feature selection methods. A retrospective, multicenter study of 385 advanced NSCLC patients amenable to ICIs was undertaken in two academic centers. Radiomic features extracted from pretreatment CT scans were used to build predictive models for PD-L1 and PFS (short-term vs. long-term survivors). We first employed the LASSO methodology followed by five feature selection methods and seven machine learning approaches to build the predictors. From our analyses, we found several combinations of feature selection methods and machine learning algorithms to achieve a similar performance. Logistic regression with ReliefF feature selection (AUC = 0.64, 0.59 in discovery and validation cohorts) and SVM with Anova F-test feature selection (AUC = 0.64, 0.63 in discovery and validation datasets) were the best-performing models to predict PD-L1 and PFS. This study elucidates the application of suitable feature selection approaches and machine learning algorithms to predict clinical endpoints using radiomics features. Through this study, we identified a subset of algorithms that should be considered in future investigations for building robust and clinically relevant predictive models., (© 2023. The Author(s).)

Single-cell technologies have revealed the complexity of the tumour immune microenvironment with unparalleled resolution 1-9 . Most clinical strategies rely on histopathological stratification of tumour subtypes, yet the spatial context of single-cell phenotypes within these stratified subgroups is poorly understood. Here we apply imaging mass cytometry to characterize the tumour and immunological landscape of samples from 416 patients with lung adenocarcinoma across five histological patterns. We resolve more than 1.6 million cells, enabling spatial analysis of immune lineages and activation states with distinct clinical correlates, including survival. Using deep learning, we can predict with high accuracy those patients who will progress after surgery using a single 1-mm 2 tumour core, which could be informative for clinical management following surgical resection. Our dataset represents a valuable resource for the non-small cell lung cancer research community and exemplifies the utility of spatial resolution within single-cell analyses. This study also highlights how artificial intelligence can improve our understanding of microenvironmental features that underlie cancer progression and may influence future clinical practice., (© 2023. The Author(s).)

Background: Immunotherapy has revolutionized clinical outcomes for patients suffering from lung cancer, yet relatively few patients sustain long-term durable responses. Recent studies have demonstrated that the tumor immune microenvironment fosters tumorous heterogeneity and mediates both disease progression and response to immune checkpoint inhibitors (ICI). As such, there is an unmet need to elucidate the spatially defined single-cell landscape of the lung cancer microenvironment to understand the mechanisms of disease progression and identify biomarkers of response to ICI., Methods: Here, in this study, we applied imaging mass cytometry to characterize the tumor and immunological landscape of immunotherapy response in non-small cell lung cancer by describing activated cell states, cellular interactions and neighborhoods associated with improved efficacy. We functionally validated our findings using preclinical mouse models of cancer treated with anti-programmed cell death protein-1 (PD-1) immune checkpoint blockade., Results: We resolved 114,524 single cells in 27 patients treated with ICI, enabling spatial resolution of immune lineages and activation states with distinct clinical outcomes. We demonstrated that CXCL13 expression is associated with ICI efficacy in patients, and that recombinant CXCL13 potentiates anti-PD-1 response in vivo in association with increased antigen experienced T cell subsets and reduced CCR2+ monocytes., Discussion: Our results provide a high-resolution molecular resource and illustrate the importance of major immune lineages as well as their functional substates in understanding the role of the tumor immune microenvironment in response to ICIs., Competing Interests: Competing interests: None declared., (© Author(s) (or their employer(s)) 2023. Re-use permitted under CC BY-NC. No commercial re-use. See rights and permissions. Published by BMJ.)

Background: Glioblastoma mortality is driven by tumour progression or recurrence despite administering a therapeutic arsenal consisting of surgical resection, radiation, and alkylating chemotherapy. The genetic changes underlying tumour progression and chemotherapy resistance are poorly understood., Methods: In this study, we sought to define the relationship between EGFR amplification status, EGFR mRNA expression, and EGFR pathway activity. We compared RNA-sequencing data from matched primary and recurrent tumour samples ( n = 40 patients, 20 with EGFR amplification)., Results: In the setting of glioblastoma recurrence, the EGFR pathway was overexpressed regardless of EGFR-amplification status, suggesting a common genomic endpoint in recurrent glioblastoma, although EGFR amplification did associate with higher EGFR mRNA expression. Three of forty patients in the study cohort had EGFR-amplified tumours and received targeted EGFR therapy. Their molecular subtypes and clinical outcomes did not significantly differ from patients who received conventional chemotherapy., Conclusion: Our findings suggest that while the EGFR amplification may confer a unique molecular profile in primary glioblastoma, pathway analysis reveals upregulation of the EGFR pathway in recurrence, regardless of amplification status. As such, the EGFR pathway may be a key mediator of glioblastoma progression.

We recently identified Zinc-finger protein 768 (ZNF768) as a novel transcription factor controlling cell fate decision downstream of Rat sarcoma virus (RAS). We showed that ZNF768 depletion impairs cell cycle progression and triggers cellular senescence, while its overexpression allows cells to bypass oncogene-induced senescence. Elevated ZNF768 levels is common in tumors, suggesting that ZNF768 may help to escape cellular senescence, sustain proliferation and promote malignant transformation. Here, we discuss these recent findings and highlight key questions emerging from our work., Competing Interests: No potential conflict of interest was reported by the author(s)., (© 2021 The Author(s). Published with license by Taylor & Francis Group, LLC.)

Background: Radiation therapy is among the most effective and commonly used therapeutic modalities of cancer treatments in current clinical practice. The fundamental paradigm that has guided radiotherapeutic regimens are 'one-size-fits-all', which are not in line with the dogma of precision medicine. While there were efforts to build radioresponse signatures using OMICS data, their ability to accurately predict in patients is still limited., Methods: We proposed to integrate two large-scale radiogenomics datasets consisting of 511 with 23 tissues and 60 cancer cell lines with 9 tissues to build and validate radiation response biomarkers. We used intrinsic radiation sensitivity, i.e., surviving fraction of cells (SF2) as the radiation response indicator. Gene set enrichment analysis was used to examine the biological determinants driving SF2. Using SF2 as a continuous variable, we used five different approaches, univariate, rank gene ensemble, rank gene multivariate, mRMR and elasticNet to build genomic predictors of radiation response through a cross-validation framework., Results: Through the pathway analysis, we found 159 pathways to be statistically significant, out of which 54 and 105 were positively and negatively enriched with SF2. More importantly, we found cell cycle and repair pathways to be enriched with SF2, which are inline with the fundamental aspects of radiation biology. With regards to the radiation response gene signature, we found that all multivariate models outperformed the univariate model with a ranking based approach performing well compared to other models, indicating complex biological processes underpinning radiation response., Conclusion: To summarize, we found biological processes underpinning SF2 and systematically compared different machine learning approaches to develop and validate predictors of radiation response. With more patient data available in the future, the clinical value of these biomarkers can be assessed that would allow for personalization of radiotherapy., (© 2021. The Author(s).)

Pulmonary neuroendocrine tumors (NETs) are a heterogeneous family of malignancies whose classification relies on morphology and mitotic rate, unlike extrapulmonary neuroendocrine tumors that require both mitotic rate and Ki-67. As mitotic count is proportional to Ki-67, it is crucial to understand if Ki-67 can complement the existing diagnostic guidelines, as well as discover the benefit of these two markers to unravel the biological heterogeneity. In this study, we investigated the association of mitotic rate and Ki-67 at gene- and pathway-level using transcriptomic data in lung NET malignancies. Lung resection tumor specimens obtained from 28 patients diagnosed with NETs were selected. Mitotic rate, Ki-67 and transcriptomic data were obtained for all samples. The concordance between mitotic rate and Ki-67 was evaluated at gene-level and pathway-level using gene expression data. Our analysis revealed a strong association between mitotic rate and Ki-67 across all samples and cell cycle genes were found to be differentially ranked between them. Pathway analysis indicated that a greater number of pathways overlapped between these markers. Analyses based on lung NET subtypes revealed that mitotic rate in carcinoids and Ki-67 in large cell neuroendocrine carcinomas provided comprehensive characterization of pathways among these malignancies. Among the two subtypes, we found distinct leading-edge gene sets that drive the enrichment signal of commonly enriched pathways between mitotic index and Ki-67. Overall, our findings delineated the degree of benefit of the two proliferation markers, and offers new layer to predict the biological behavior and identify high-risk patients using a more comprehensive diagnostic workup., Competing Interests: CONFLICTS OF INTEREST Authors have no conflicts of interest to declare., (Copyright: © 2021 Manem et al.)

Motivation: High-throughput molecular profiles of human cells have been used in predictive computational approaches for stratification of healthy and malignant phenotypes and identification of their biological states. In this regard, pathway activities have been used as biological features in unsupervised and supervised learning schemes., Results: We developed SIGN (Similarity Identification in Gene expressioN), a flexible open-source R package facilitating the use of pathway activities and their expression patterns to identify similarities between biological samples. We defined a new measure, the transcriptional similarity coefficient, which captures similarity of gene expression patterns, instead of quantifying overall activity, in biological pathways between the samples. To demonstrate the utility of SIGN in biomedical research, we establish that SIGN discriminates subtypes of breast tumors and patients with good or poor overall survival. SIGN outperforms the best models in DREAM challenge in predicting survival of breast cancer patients using the data from the Molecular Taxonomy of Breast Cancer International Consortium. In summary, SIGN can be used as a new tool for interrogating pathway activity and gene expression patterns in unsupervised and supervised learning schemes to improve prognostic risk estimation for cancer patients by the biomedical research community., Availability and Implementation: An open-source R package is available (https://cran.r-project.org/web/packages/SIGN/)., (© The Author(s) 2019. Published by Oxford University Press. All rights reserved. For permissions, please e-mail: journals.permissions@oup.com.)

A wealth of transcriptomic and clinical data on solid tumours are under-utilized due to unharmonized data storage and format. We have developed the MetaGxData package compendium, which includes manually-curated and standardized clinical, pathological, survival, and treatment metadata across breast, ovarian, and pancreatic cancer data. MetaGxData is the largest compendium of curated transcriptomic data for these cancer types to date, spanning 86 datasets and encompassing 15,249 samples. Open access to standardized metadata across cancer types promotes use of their transcriptomic and clinical data in a variety of cross-tumour analyses, including identification of common biomarkers, and assessing the validity of prognostic signatures. Here, we demonstrate that MetaGxData is a flexible framework that facilitates meta-analyses by using it to identify common prognostic genes in ovarian and breast cancer. Furthermore, we use the data compendium to create the first gene signature that is prognostic in a meta-analysis across 3 cancer types. These findings demonstrate the potential of MetaGxData to serve as an important resource in oncology research, and provide a foundation for future development of cancer-specific compendia.

Understanding the tumor immune microenvironment (TIME) promises to be key for optimal cancer therapy, especially in triple-negative breast cancer (TNBC). Integrating spatial resolution of immune cells with laser capture microdissection gene expression profiles, we defined distinct TIME stratification in TNBC, with implications for current therapies including immune checkpoint blockade. TNBCs with an immunoreactive microenvironment exhibited tumoral infiltration of granzyme B+CD8+ T cells (GzmB+CD8+ T cells), a type 1 IFN signature, and elevated expression of multiple immune inhibitory molecules including indoleamine 2,3-dioxygenase (IDO) and programmed cell death ligand 1 (PD-L1), and resulted in good outcomes. An "immune-cold" microenvironment with an absence of tumoral CD8+ T cells was defined by elevated expression of the immunosuppressive marker B7-H4, signatures of fibrotic stroma, and poor outcomes. A distinct poor-outcome immunomodulatory microenvironment, hitherto poorly characterized, exhibited stromal restriction of CD8+ T cells, stromal expression of PD-L1, and enrichment for signatures of cholesterol biosynthesis. Metasignatures defining these TIME subtypes allowed us to stratify TNBCs, predict outcomes, and identify potential therapeutic targets for TNBC.

Background: Immune checkpoint inhibitors (ICIs) demonstrate unprecedented efficacy in multiple malignancies; however, the mechanisms of sensitivity and resistance are poorly understood and predictive biomarkers are scarce. INSPIRE is a phase 2 basket study to evaluate the genomic and immune landscapes of peripheral blood and tumors following pembrolizumab treatment., Methods: Patients with incurable, locally advanced or metastatic solid tumors that have progressed on standard therapy, or for whom no standard therapy exists or standard therapy was not deemed appropriate, received 200 mg pembrolizumab intravenously every three weeks. Blood and tissue samples were collected at baseline, during treatment, and at progression. One core biopsy was used for immunohistochemistry and the remaining cores were pooled and divided for genomic and immune analyses. Univariable analysis of clinical, genomic, and immunophenotyping parameters was conducted to evaluate associations with treatment response in this exploratory analysis., Results: Eighty patients were enrolled from March 21, 2016 to June 1, 2017, and 129 tumor and 382 blood samples were collected. Immune biomarkers were significantly different between the blood and tissue. T cell PD-1 was blocked (≥98%) in the blood of all patients by the third week of treatment. In the tumor, 5/11 (45%) and 11/14 (79%) patients had T cell surface PD-1 occupance at weeks six and nine, respectively. The proportion of genome copy number alterations and abundance of intratumoral 4-1BB+ PD-1+ CD8 T cells at baseline (P < 0.05), and fold-expansion of intratumoral CD8 T cells from baseline to cycle 2-3 (P < 0.05) were associated with treatment response., Conclusion: This study provides technical feasibility data for correlative studies. Tissue biopsies provide distinct data from the blood and may predict response to pembrolizumab.

There has been a paradigm shift in translational oncology with the advent of novel molecular diagnostic tools in the clinic. However, several challenges are associated with the integration of these sophisticated tools into clinical oncology and daily practice. High-throughput profiling at the DNA, RNA and protein levels (omics) generate a massive amount of data. The analysis and interpretation of these is non-trivial but will allow a more thorough understanding of cancer. Linear modelling of the data as it is often used today is likely to limit our understanding of cancer as a complex disease, and at times under-performs to capture a phenotype of interest. Network science and systems biology-based approaches, using machine learning and network science principles, that integrate multiple data sources, can uncover complex changes in a biological system. This approach will integrate a large number of potential biomarkers in preclinical studies to better inform therapeutic decisions and ultimately make substantial progress towards precision medicine. It will however require development of a new generation of clinical trials. Beyond discussing the challenges of high-throughput technologies, this review will develop a framework on how to implement a network science approach in new clinical trial designs in order to advance cancer care., (Copyright © 2017 Elsevier Ltd. All rights reserved.)

In the past few years, proton therapy has taken the centre stage in treating various tumour types. The primary contribution of this study is to investigate the tumour control probability (TCP), relapse time and the corresponding secondary cancer risks induced by proton beam radiation therapy. We incorporate tumour relapse kinetics into the TCP framework and calculate the associated second cancer risks. To calculate proton therapy-induced secondary cancer induction, we used the well-known biologically motivated mathematical model, initiation-inactivation-proliferation formalism. We used the available in vitro data for the linear energy transfer (LET) dependence of cell killing and mutation induction parameters. We evaluated the TCP and radiation-induced second cancer risks for protons in the clinical range of LETs, i.e. approximately 8 $\mathrm{keV/\mu m}$ for the tumour volume and 1-3 $\mathrm{keV/\mu m}$ for the organs at risk. This study may serve as a framework for further work in this field and elucidates proton-induced TCP and the associated secondary cancer risks, not previously reported in the literature. Although studies with a greater number of cell lines would reduce uncertainties within the model parameters, we argue that the theoretical framework presented within is a sufficient rationale to assess proton radiation TCP, relapse and carcinogenic effects in various treatment plans. We show that compared with photon therapy, proton therapy markedly reduces the risk of secondary malignancies and for equivalent dosing regimens achieves better tumour control as well as a reduced primary recurrence outcome, especially within a hypo-fractionated regimen.

Drug combinations have been proposed as a promising therapeutic strategy to overcome drug resistance and improve efficacy of monotherapy regimens in cancer. This strategy aims at targeting multiple components of this complex disease. Despite the increasing number of drug combinations in use, many of them were empirically found in the clinic, and the molecular mechanisms underlying these drug combinations are often unclear. These challenges call for rational, systematic approaches for drug combination discovery. Although high-throughput screening of single-agent therapeutics has been successfully implemented, it is not feasible to test all possible drug combinations, even for a reduced subset of anticancer drugs. Hence, in vitro and in vivo screening of a large number of drug combinations are not practical. Therefore, devising computational methods to efficiently explore the space of drug combinations and to discover efficacious combinations has attracted a lot of attention from the scientific community in the past few years. Nevertheless, in the absence of consensus regarding the computational approaches used to predict efficacious drug combinations, a plethora of methods, techniques and hypotheses have been developed to date, while the research field lacks an elaborate categorization of the existing computational methods and the available data sources. In this manuscript, we review and categorize the state-of-the-art computational approaches for drug combination prediction, and elaborate on the limitations of these methods and the existing challenges. We also discuss about the recent pan-cancer drug combination data sets and their importance in revising the available methods or developing more performant approaches.

Objectives: We sought to investigate the tissue specificity of drug sensitivities in large-scale pharmacological studies and compare these associations to those found in drug clinical indications., Materials and Methods: We leveraged the curated cell line response data from PharmacoGx and applied an enrichment algorithm on drug sensitivity values' area under the drug dose-response curves (AUCs) with and without adjustment for general level of drug sensitivity., Results: We observed tissue specificity in 63% of tested drugs, with 8% of total interactions deemed significant (false discovery rate <0.05). By restricting the drug-tissue interactions to those with AUC > 0.2, we found that in 52% of interactions, the tissue was predictive of drug sensitivity (concordance index > 0.65). When compared with clinical indications, the observed overlap was weak (Matthew correlation coefficient, MCC = 0.0003, P > .10)., Discussion: While drugs exhibit significant tissue specificity in vitro, there is little overlap with clinical indications. This can be attributed to factors such as underlying biological differences between in vitro models and patient tumors, or the inability of tissue-specific drugs to bring additional benefits beyond gold standard treatments during clinical trials., Conclusion: Our meta-analysis of pan-cancer drug screening datasets indicates that most tested drugs exhibit tissue-specific sensitivities in a large panel of cancer cell lines. However, the observed preclinical results do not translate to the clinical setting. Our results suggest that additional research into showing parallels between preclinical and clinical data is required to increase the translational potential of in vitro drug screening., (© The Author 2017. Published by Oxford University Press on behalf of the American Medical Informatics Association. All rights reserved. For Permissions, please email: journals.permissions@oup.com)

Several studies have shown that pediatric patients have an increased risk of developing a secondary malignancy several decades after treatment with radiotherapy and chemotherapy. In this work, we use a biologically motivated mathematical formalism to estimate the relative risks of breast, lung and thyroid cancers in childhood cancer survivors due to concurrent therapy regimen. This model specifically includes possible organ-specific interaction between radiotherapy and chemotherapy. The model predicts relative risks for developing secondary cancers after chemotherapy in breast, lung and thyroid tissues, and compared with the epidemiological data. For a concurrent therapy protocol, our model predicted relative risks of 3.2, 9.3, 4.5 as compared to the clinical data, i.e., 1.4, 8.0, 2.3 for secondary breast, lung and thyroid cancer risks, respectively. The extracted chemotherapy mutation induction rates for breast, lung and thyroid are 10 &#8722;9 , 0.5 &#215; 10 &#8722;6 , 0.9 &#215; 10 &#8722;7 respectively. We found that there exists no synergistic interaction between radiation and chemotherapy for neither mutation induction nor cell kill in lung tissue, but there is an interaction in cell kill for the breast and thyroid organs. These findings help understand the risks of current clinical protocols and might provide rational guidance to develop future multi-modality treatment protocols to minimize secondary cancer risks., Competing Interests: The authors declare no conflict of interest.

Although the survival rate of cancer patients has significantly increased due to advances in anti-cancer therapeutics, one of the major side effects of these therapies, particularly radiotherapy, is the potential manifestation of radiation-induced secondary malignancies. In this work, a novel evolutionary stochastic model is introduced that couples short-term formalism (during radiotherapy) and long-term formalism (post-treatment). This framework is used to estimate the risks of second cancer as a function of spontaneous background and radiation-induced mutation rates of normal and pre-malignant cells. By fitting the model to available clinical data for spontaneous background risk together with data of Hodgkin's lymphoma survivors (for various organs), the second cancer mutation rate is estimated. The model predicts a significant increase in mutation rate for some cancer types, which may be a sign of genomic instability. Finally, it is shown that the model results are in agreement with the measured results for excess relative risk (ERR) as a function of exposure age and that the model predicts a negative correlation of ERR with increase in attained age. This novel approach can be used to analyze several radiotherapy protocols in current clinical practice and to forecast the second cancer risks over time for individual patients.