RET Alterations Differentiate Molecular Profile of Medullary Thyroid Cancer

PURPOSE

Medullary thyroid cancer (MTC) is a rare cancer originating from parafollicular C cells of the thyroid gland. Therapeutically relevant alterations in MTC are predominantly reported in RET oncogene, and lower-frequency alterations are reported in KRAS and BRAF. Nevertheless, there is an unmet need existing to analyze the MTC in the Indian cohort by using in-depth sequencing techniques that go beyond the identification of known therapeutic biomarkers.

MATERIALS AND METHODS

Here, we characterize MTC using integrative whole-exome and whole-transcriptome sequencing of 32 MTC tissue samples. We performed clinically relevant variant analysis, molecular pathway analysis, tumor immune-microenvironment analysis, and structural characterization of RET novel mutation.

RESULTS

Mutational landscape analysis shows expected RET mutations in 50% of the cases. Furthermore, we observed mutations in known cancer genes like KRAS, HRAS, SF3B1, and BRAF to be altered only in the RET-negative cohort. Pathway analysis showed differential enrichment of mutations in transcriptional deregulation genes in the RET-negative cohort. Furthermore, we observed novel RET kinase domain mutation Y900S showing affinity to RET inhibitors accessed via molecular docking and molecular dynamics simulation.

CONCLUSION

Altogether, this study provides a detailed genomic characterization of patients with MTC of Indian origin, highlighting the possible utility of targeted therapies in this disease.

Clinical characteristics, outcomes and prognostic factors in KRAS mutant lung cancers: experience from a tertiary care cancer center in India

Objectives

Kirsten rat sarcoma viral oncogene homologue (KRAS) mutations in lung cancers, long considered untargetable, have had a recent rise in interest due to promising data of agents targeting KRAS p.G12C. As Indian data are scarce, we sought to identify baseline clinical characteristics, prognostic factors and outcomes of lung cancer patients with KRAS mutations at our hospital.

Methods

Patients with KRAS mutant lung cancers treated at our institute from 2016 to 2022 were analysed.

Results

133 patients with KRAS mutant lung cancers were identified. Median age was 57 (interquartile range 28-78) years, and 58 (43.6%) were smokers. 17 (12.7%) had brain metastases. The commonest variant was p.G12C, seen in 53 (39.8%) patients. Six (4.5%) had programmed death ligand 1 (PDL-1) expression >50% by Ventana SP263 PDL-1 assay, and 13 (9.7%) had epidermal growth factor mutation. Of 92 patients with available treatment details, the majority received intravenous chemotherapy, nine (9.8%) received tyrosine kinase inhibitors and four (4.4%) received immunotherapy (pembrolizumab). Median progression-free survival (PFS) with first-line therapy was 6 (95% confidence interval (CI) 2.8-9.2) months and median overall survival (OS) was 12 (CI 9.2-14.8) months. The incidence of brain metastases was higher in patients with G12C mutations (p = 0.025). Brain metastases (HR: 3.57, p < 0.001), Eastern Cooperative Oncology Group performance status (PS) ≥ 2 (HR: 2.13, p = 0.002) and G12C mutation (HR: 1.84, p = 0.011) were associated with inferior PFS, while brain metastases (HR: 4.6, p < 0.001), PS ≥ 2 (HR: 2.33, p = 0.001) and G12C mutation (HR: 1.93, p = 0.01) were associated with inferior OS.

Conclusion

This is the largest dataset of KRAS mutant lung cancers from India. Brain metastases were higher in patients with G12C mutations and associated with poorer PFS and OS. G12C mutation and PS ≥ 2 were also associated with inferior PFS and OS. Experience with targeted therapy for KRAS mutations remains an area of future exploration due to the unavailability of these agents in India.

Natural History of Germline BRCA1 Mutated and BRCA Wild-type Triple-negative Breast Cancer

We report a deep next-generation sequencing analysis of 13 sequentially obtained tumor samples, eight sequentially obtained circulating tumor DNA (ctDNA) samples and three germline DNA samples over the life history of 3 patients with triple-negative breast cancer (TNBC), 2 of whom had germline pathogenic BRCA1 mutation, to unravel tumor evolution. Tumor tissue from all timepoints and germline DNA was subjected to whole-exome sequencing (WES), custom amplicon deep sequencing (30,000X) of a WES-derived somatic mutation panel, and SNP arrays for copy-number variation (CNV), while whole transcriptome sequencing (RNA-seq) was performed only on somatic tumor.There was enrichment of homologous recombination deficiency signature in all tumors and widespread CNV, which remained largely stable over time. Somatic tumor mutation numbers varied between patients and within each patient (range: 70-216, one outlier). There was minimal mutational overlap between patients with TP53 being the sole commonly mutated gene, but there was substantial overlap in sequential samples in each patient. Each patient's tumor contained a founding ("stem") clone at diagnosis, which persisted over time, from which all other clones ("subclone") were derived ("branching evolution"), which contained mutations in well-characterized cancer-related genes like PDGFRB, ARID2, TP53 (Patient_02), TP53, BRAF, BRIP1, CSF3R (Patient_04), and TP53, APC, EZH2 (Patient_07). Including stem and subclones, tumors from all patients were polyclonal at diagnosis and during disease progression. ctDNA recapitulated most tissue-derived stem clonal and subclonal mutations while detecting some additional subclonal mutations. RNA-seq revealed a stable basal-like pattern, with most highly expressed variants belonging to stem clone.

SIGNIFICANCE

In germline BRCA1 mutated and BRCA wild-type patients, TNBC shows a branching evolutionary pattern of mutations with a single founding clone, are polyclonal throughout their disease course, and have widespread copy-number aberrations. This evolutionary pattern may be associated with treatment resistance or sensitivity and could be therapeutically exploited.

Whole exome sequencing uncovers HRAS mutations as potential mediators of resistance to metronomic chemotherapy

OBJECTIVES

The aim of this pilot study is to identify the genetic factors that contribute to the response of metronomic chemotherapy in head and neck squamous cell carcinoma (HNSCC) patients using whole-exome sequencing (WES). This study would facilitate the identification of predictive biomarkers, which would enable personalized treatment strategies and improve treatment outcomes for patients with HNSCC.

MATERIALS AND METHODS

We have selected patients with recurrent head and neck cancer who underwent metronomic chemotherapy. Sequential tumor biopsies were collected from the patients at different stages of treatment to capture the genomic alterations and tumor evolution during metronomic chemotherapy and sequenced using WES.

RESULTS

We identified several known HNSCC hallmark genes reported in COSMIC, including KMT2B, NOTCH1, FAT1, TP53, HRAS, CASP8, and CDKN2A. Copy number alteration analysis revealed amplifications and deletions in several oncogenic and tumor suppressor genes. COSMIC Mutational Signature 15 associated with defective DNA mismatch repair was enriched in 73% of HNSCC samples. Further, the comparison of genomic alterations between responders and non-responders identified HRAS gene uniquely mutated in non-responders that could potentially contribute to resistance against metronomic chemotherapy.

DISCUSSION

Our findings corroborate the molecular heterogeneity of recurrent HNSCC tumors and establish an association between HRAS mutations and resistance to metronomic chemotherapy, suggesting HRAS as a potential therapeutic target. Combining HRAS inhibitors with metronomic regimens could improve treatment sensitivity in HRAS-mutated HNSCC patients. Further studies are needed to fully elucidate the genomic mechanisms underlying the response to metronomic chemotherapy.

Association of Cutibacterium acnes with human thyroid cancer

Introduction

The diverse subtypes of thyroid carcinoma have distinct clinical outcomes despite a comparable spectrum of underlying genetic alterations. Beyond genetic alterations, sparse efforts have been made to characterize the microbes associated with thyroid cancer. In this study, we examine the microbial profile of thyroid cancer.

Methods

We sequenced the whole transcriptome of 70 thyroid cancers (40 papillary and 30 anaplastic). Using Infectious Pathogen Detector IPD 2.0, we analysed the relative abundance of 1060 microbes across 70 tumours from patients with thyroid cancer against 118 tumour samples from patients with breast, cervical, colorectal, and tongue cancer.

Results

Our analysis reveals a significant prevalence of Cutibacterium acnes in 58.6% thyroid cancer samples compared to other cancer types (p=0.00038). Immune cell fraction analysis between thyroid cancer samples with high and low Cutibacterium loads identify enrichment of immunosuppressive cells, including Tregs (p=0.015), and other anti-inflammatory cytokines in the tumour microenvironment, suggesting an immune evasion/immunosuppression milieu is associated with the infection. A higher burden of Cutibacterium acnes was also found to be associated with poor survival defining a distinct sub-group of thyroid cancer.

Conclusion

Cutibacterium acnes is associated with immune suppression and poor prognosis in a subpopulation of thyroid cancer. This study may help design novel therapeutic measures involving appropriate antibiotics to manage the disease better.

Real-World Evidence of EGFR Targeted Therapy in NSCLC- A Brief Report of Decade Long Single Center Experience

The significance of EGFR targeted therapy in the lung adenocarcinoma is paramount. Several controlled clinical trials have reported considerable survival of EGFR mutation positive patients on receiving the EGFR tyrosine kinase inhibitor (TKI). However, the real-world evidence of benefits of EGFR TKI would be further useful to understand how the designated therapeutic regimen benefits the patients. In this study, we report a decade long real-world evidence of EGFR molecular testing in lung cancer at Tata Memorial Hospital (Mumbai, India). Laboratory and hospital records containing basic demographic details, clinical characteristics, treatment regimen, survival outcome were collected retrospectively. Statistical association and survival analysis were performed using the R programming. The cohort includes 9,053 lung cancer patients tested for EGFR mutations during 2011 to 2019. Baseline T790M and compound mutations were the only mutations observed co-occurring while all other EGFR mutations were mutually exclusive. Furthermore, the baseline T790M were also observed to be associated with TTF1 positivity, smoking and local metastasis. Overall survival of the patients harboring co-occurring compound mutations was significantly lesser than the other EGFR positive patients. Overall, our study suggests that EGFR TKI may provide real-world benefit to the lung cancer patients harboring mutually exclusive EGFR mutations. On the other hand, further systematic study is essential to develop better therapeutic regimen for co-occurring baseline EGFR T790M and other compound mutations.

Understanding the Impact of Population and Cancer Type on Tumor Mutation Burden Scores: A Comprehensive Whole-Exome Study in Cancer Patients From India

PURPOSE

The purpose of this study was to understand the impact of population diversity and geographic variation on tumor mutation burden (TMB) scores across cancers and its implication on stratification of patients for immune checkpoint inhibitor (ICI) therapy.

MATERIALS AND METHODS

This retrospective study used whole-exome sequencing (WES) to profile 1,233 Indian patients with cancer across 30 different cancer types and to estimate their TMB scores. A WES-based pipeline was adopted, along with an indigenously developed strategy for arriving at true somatic mutations. A robust unsupervised machine learning approach was used to understand the distribution of TMB scores across different populations and within the population.

RESULTS

The results of the study showed a biphasic distribution of TMB scores in most cancers, with different threshold scores across cancer types. Patients with cancer in India had higher TMB scores compared with the Caucasian patients. We also observed that the TMB score value at 90th percentile (predicting high efficacy to ICI) was high in four different cancer types (sarcoma, ovary, head and neck, and breast) in the Indian cohort as compared with The Cancer Genome Atlas or public cohort. However, in lung and colorectal cancers, the TMB score distribution was similar between the two population cohorts.

CONCLUSION

The findings of this study indicate that it is crucial to benchmark both cancer-specific and population-specific TMB distributions to establish a TMB threshold for each cancer in various populations. Additional prospective studies on much larger population across different cancers are warranted to validate this observation to become the standard of care.

In silico structural analysis of secretory clusterin to assess pathogenicity of mutations identified in the evolutionarily conserved regions

Clusterin (CLU) is a secreted glycoprotein, heterodimeric in nature, and is expressed in a wide variety of tissues and body fluids such as serum and plasma. CLU has also been known to be a promising biomarker for cell death, malignancy, cancer progression, and resistance development. However, the lack of a CLU crystal structure obstructs understanding the possible role of reported mutations on the structure, and the subsequent effects on downstream signaling pathways and cancer progression. Considering the importance of crystal structure, a model structure of the pre-secretory isoform of CLU was built to predict the effect of mutations at the molecular level. Ab initio model was built using RaptorX, and loop refinement and energy minimization were carried out with ModLoop, ModRefiner, and GalaxyWeb servers. The cancer associated mutational spectra of CLU was retrieved from the cBioPortal server and 117 unique missense mutations were identified. Evolutionarily conserved regions and pathogenicity of mutations identified in CLU were analyzed using ConSurf and Rhapsody, respectively. Furthermore, sequence and structure-based mutational analysis were carried out with iSTABLE, DynaMut and PremPS servers. Molecular dynamics simulations were carried out with GROMACS for 50 ns to determine the stability of the wild type and mutant protein structures. A dynamically stable model structure of pre-secretory CLU (psCLU) which has high concurrence with the sequence based secondary structure predictions has been explored. Changes in the intra-atomic interactions and folding pattern between wild type and mutant structures were observed. To our conclusion, eleven mutations with the highest structural and functional significance have been predicted to have pathogenic and deleterious effects.Communicated by Ramaswamy H. Sarma.

Progesterone modulates the DSCAM-AS1/miR-130a/ESR1 axis to suppress cell invasion and migration in breast cancer

BACKGROUND

A preoperative-progesterone intervention increases disease-free survival in patients with breast cancer, with an unknown underlying mechanism. We elucidated the role of non-coding RNAs in response to progesterone in human breast cancer.

METHODS

Whole transcriptome sequencing dataset of 30 breast primary tumors (10 tumors exposed to hydroxyprogesterone and 20 tumors as control) were re-analyzed to identify differentially expressed non-coding RNAs followed by real-time PCR analyses to validate the expression of candidates. Functional analyses were performed by genetic knockdown, biochemical, and cell-based assays.

RESULTS

We identified a significant downregulation in the expression of a long non-coding RNA, Down syndrome cell adhesion molecule antisense DSCAM-AS1, in response to progesterone treatment in breast cancer. The progesterone-induced expression of DSCAM-AS1 could be effectively blocked by the knockdown of progesterone receptor (PR) or treatment of cells with mifepristone (PR-antagonist). We further show that knockdown of DSCAM-AS1 mimics the effect of progesterone in impeding cell migration and invasion in PR-positive breast cancer cells, while its overexpression shows an opposite effect. Additionally, DSCAM-AS1 sponges the activity of miR-130a that regulates the expression of ESR1 by binding to its 3'-UTR to mediate the effect of progesterone in breast cancer cells. Consistent with our findings, TCGA analysis suggests that high levels of miR-130a correlate with a tendency toward better overall survival in patients with breast cancer.

CONCLUSION

This study presents a mechanism involving the DSCAM-AS1/miR-130a/ESR1 genomic axis through which progesterone impedes breast cancer cell invasion and migration. The findings highlight the utility of progesterone treatment in impeding metastasis and improving survival outcomes in patients with breast cancer.

Integrated proteomics and phosphoproteomics revealed druggable kinases in neoadjuvant chemotherapy resistant tongue cancer

Tongue squamous cell carcinoma is an aggressive oral cancer with a high incidence of metastasis and poor prognosis. Most of the oral cavity cancer patients present in clinics with locally advanced unresectable tumors. Neoadjuvant treatment is beneficial for these individuals as it reduces the tumor size aiding complete resection. However, patients develop therapy resistance to the drug regimen. In this study, we explored the differential expression of proteins and altered phosphorylation in the neoadjuvant chemotherapy resistant tongue cancer patients. We integrated the proteomic and phosphoproteomic profiles of resistant (n = 4) and sensitive cohorts (n = 4) and demonstrated the differential expression and phosphorylation of proteins in the primary tissue of the respective subject groups. We observed differential and extensive phosphorylation of keratins such as KRT10 and KRT1 between the two cohorts. Furthermore, our study revealed a kinase signature associated with neoadjuvant chemotherapy resistance. Kinases such as MAPK1, AKT1, and MAPK3 are predicted to regulate the resistance in non-responders. Pathway analysis showed enrichment of Rho GTPase signaling and hyperphosphosphorylation of proteins involved in cell motility, invasion, and drug resistance. Targeting the kinases could help with the clinical management of neoadjuvant chemotherapy-resistant tongue cancer.

In vivo tumor immune microenvironment phenotypes correlate with inflammation and vasculature to predict immunotherapy response

Response to immunotherapies can be variable and unpredictable. Pathology-based phenotyping of tumors into ‘hot’ and ‘cold’ is static, relying solely on T-cell infiltration in single-time single-site biopsies, resulting in suboptimal treatment response prediction. Dynamic vascular events (tumor angiogenesis, leukocyte trafficking) within tumor immune microenvironment (TiME) also influence anti-tumor immunity and treatment response. Here, we report dynamic cellular-level TiME phenotyping in vivo that combines inflammation profiles with vascular features through non-invasive reflectance confocal microscopic imaging. In skin cancer patients, we demonstrate three main TiME phenotypes that correlate with gene and protein expression, and response to toll-like receptor agonist immune-therapy. Notably, phenotypes with high inflammation associate with immunostimulatory signatures and those with high vasculature with angiogenic and endothelial anergy signatures. Moreover, phenotypes with high inflammation and low vasculature demonstrate the best treatment response. This non-invasive in vivo phenotyping approach integrating dynamic vasculature with inflammation serves as a reliable predictor of response to topical immune-therapy in patients.

Standard assessment of immune infiltration of biopsies is not sufficient to accurately predict response to immunotherapy. Here, the authors show that reflectance confocal microscopy can be used to quantify dynamic vasculature and inflammatory features to better predict treatment response in skin cancers.

Impact of Molecular Tumor Board on the Clinical Management of Patients With Cancer

PURPOSE

Multidisciplinary molecular tumor boards (MTBs) help in interpreting complex genomic data generated by molecular tumor profiling and improve patients' access to targeted therapies. The purpose of this study was to assess the impact of our institution's MTB on the clinical management of patients with cancer.

METHODS

This study was conducted at a tertiary cancer center in India. Cases to be discussed in the MTB were identified by molecular pathologists, scientists, or oncologists. On the basis of the clinical data and molecular test reports, a course of clinical management was recommended and made available to the treating oncologist. We determined the proportion of patients who were recommended a change in the clinical management. We also assessed compliance of the treating oncologists with MTB recommendations.

RESULTS

There were 339 discussions for 328 unique patients. The median age of the cohort was 54 years (range 17-87), and the majority of the patients were men (65.1%). Of 339 cases, 133 (39.2%) were recommended continuation of ongoing therapy while the remaining 206 (60.7%) were recommended a change in clinical management. Compliance with MTB recommendations for a change in clinical management was 58.5% (79 of 138 evaluable cases). Compliance and implementation for MTB's recommendation to start a new therapy in 104 evaluable cases were 60.5% and 44.2%, respectively. A total of 248 biopsies had at least one actionable mutation. A total of 646 mutations were identified in the cohort, with EGFR being the most frequently altered gene.

CONCLUSION

MTBs help in interpreting results of molecular tests, understanding the significance of molecular abnormalities, and assessing the benefits of available targeted therapies and clinical trials in the management of patients with targetable genetic alterations.

Treatment pattern and outcomes in de novo T790M-mutated non-small cell lung cancer

Introduction

Methods

Results

Conclusion

Fusobacterium nucleatum is associated with inflammation and poor survival in early-stage HPV-negative tongue cancer

Persistent pathogen infection is a known cause of malignancy, although with sparse systematic evaluation across tumor types. We present a comprehensive landscape of 1060 infectious pathogens across 239 whole exomes and 1168 transcriptomes of breast, lung, gallbladder, cervical, colorectal, and head and neck tumors. We identify known cancer-associated pathogens consistent with the literature. In addition, we identify a significant prevalence of Fusobacterium in head and neck tumors, comparable to colorectal tumors. The Fusobacterium-high subgroup of head and neck tumors occurs mutually exclusive to human papillomavirus, and is characterized by overexpression of miRNAs associated with inflammation, elevated innate immune cell fraction and nodal metastases. We validate the association of Fusobacterium with the inflammatory markers IL1B, IL6 and IL8, miRNAs hsa-mir-451a, hsa-mir-675 and hsa-mir-486-1, and MMP10 in the tongue tumor samples. A higher burden of Fusobacterium is also associated with poor survival, nodal metastases and extracapsular spread in tongue tumors defining a distinct subgroup of head and neck cancer.

Abstract 2814: Dynamic imaging of tumor-immune microenvironment (TiME) and microvasculature identifies ‘hot' and ‘cold' tumor phenotypes in vivo in patients

Investigating the dynamic crosstalk between the tumor-immune microenvironment (TiME) and microvasculature in vivo in patients can lead to important insights into the underlying biology, help identify tumor phenotypes and reveal attractive druggable targets.

Dynamic non-invasive label-free imaging of TiME and microvasculature in real-time directly in patients using reflectance confocal microscopy (RCM) was investigated on 60 skin cancer patients (basal cell carcinoma, BCC; squamous cell carcinoma, SCC), followed by automated and machine-learning based quantification of TiME and microvasculature features such as vascular density, leukocyte trafficking and immune cell density. Manual (two readers) and histopathological evaluation (dermatopathologist) of these features was also performed. Molecular correlation of imaging features and phenotypes was performed using anti-CD3/anti-CD20 IHC staining for tertiary lymphoid structures (TLS) and total lymphocyte density (n=33), flow cytometry for immune cells (n=3), and differential RNA expression (n=14). Correlation of RCM features and phenotypes at baseline (before treatment) with treatment response was also evaluated on 9 cancer lesions undergoing topical immunotherapy imiquimod. High agreement for feature presence on RCM and Histology, and manual and automated RCM features was observed. Unsupervised clustering on total TiME and microvasculature features on RCM using principal component analysis (PCA) indicates four distinct tumor phenotypes (PCA 1). The phenotype with high inflammation, high trafficking and higher density of vessels or the denoted ‘hot' phenotype correlated with higher activated CD8+ Granzyme B+ cells (67% of total CD8+cells). The clustering pattern on RCM was compared to TLS and lymphocyte density (PCA 2) and gene expression following CIBERSORT analysis (PCA 3). The clustering in RCM correlated better with gene expression (PCA 1 and 3, 100% agreement) than TLS and lymphocyte density (PCA 1 and 2, 86% agreement). The ‘hot' phenotype in RCM correlated with higher immune gene signatures and higher TLS/lymphocyte density. Increased plasma, CD8, activated CD4 memory and activated NK cells, M1 macrophages and monocytes, along with up-regulation of JAK-STAT, chemokine and cell adhesion signaling cascade were found in the ‘hot' RCM phenotype. Statistical modeling for correlating phenotypes with treatment outcomes was performed using principal component-linear discriminant analysis (PC-LDA). Two responders with tumor regression were predicted as ‘hot' phenotype while the non-responding patients (remaining 7) were classified as cold phenotype; suggesting that RCM 'hot' phenotype correlates with better treatment response. Thus, we demonstrate the potential utility of noninvasive RCM imaging in identifying ‘hot' and ‘cold' tumor phenotypes directly in patients.

Citation Format: Aditi Sahu, Melissa Gill, Miguel Cordova, Anthony Santella, Kivanc Kose, Teguru Tembo, Anabel Alfonso, Pratik Chandrani, Christi Fox, Salvador Gonzalez, Nicholas Kurtansky, Melissa Pulitzer, William Phillips, Madison Li, Kimeil King, Stephen Dusza, Shuaitong Liu, Ning Yang, Haaris Jilani, Paras Mehta, Ashfaq Marghoob, Allan Halpern, Anthony Rossi, Liang Deng, Chih-Shan Jason Chen, Milind Rajadhyaksha. Dynamic imaging of tumor-immune microenvironment (TiME) and microvasculature identifies ‘hot' and ‘cold' tumor phenotypes in vivo in patients [abstract]. In: Proceedings of the American Association for Cancer Research Annual Meeting 2021; 2021 Apr 10-15 and May 17-21. Philadelphia (PA): AACR; Cancer Res 2021;81(13_Suppl):Abstract nr 2814.

Landscape of toll-like receptors expression in tumor microenvironment of triple negative breast cancer (TNBC): Distinct roles of TLR4 and TLR8

TNBC is the most aggressive and hormone receptor-negative subtype of breast cancer with molecular heterogeneity in bulk tumors hindering effective treatment. Toll-like receptors (TLRs) have the potential to ignite diverse immune responses in the tumor microenvironment (TME). This encouraged us to screen their transcript expression in the publically available TCGA datasets. Reported molecular subtypes of TNBC may represent different TMEs and we observed differentially expressed TLRs (DETs) i.e. TLR3/4/6/8/9 have unique expression pattern in the TNBC subtypes, particularly in Immunomodulatory (IM) TNBC subtype. We then dissected expression of the DETs in immune and other components of the TME. TLR4 and TLR8 showed significant (p-value ≤ 0.05) negative partial correlation with tumor purity compared to other DETs. Interestingly, TLR4 and TLR8 expression showed a significant (adjusted p-value ≤ 0.05) correlation with different subsets of immune infiltrating cells having the highest correlation with monocytes/macrophage/dendritic cell populations mediating both innate and adaptive response in TNBC. The co-expression network identified genes correlated with these immune cells. Further, GSEA analysis of co-expressed genes showed a significant association of TLR8 partners with 'Peptide ligand binding', 'Gά-signaling', and 'Cytokine-cytokine interaction' while TLR4 associated genes correlated with 'Adaptive immune system' and 'Systemic lupus erythematosus' interactome. Finally, the expression of TLR4 protein was validated in a panel of TNBC cell lines. TLR4 expression in chemoresponsive TNBC was also validated in TNBC cell lines upon Paclitaxel (PTX) treatment. Collectively, the present study identified specific DETs in TNBC and discovered a prospective role of TLR4 and TLR8 in the maintenance of tumor-immune-microenvironment.

Molecular characterization of lung squamous cell carcinoma tumors reveals therapeutically relevant alterations

INTRODUCTION

Unlike lung adenocarcinoma patients, there is no FDA-approved targeted-therapy likely to benefit lung squamous cell carcinoma patients.

MATERIALS AND METHODS

We performed survival analyses of lung squamous cell carcinoma patients harboring therapeutically relevant alterations identified by whole exome sequencing and mass spectrometry-based validation across 430 lung squamous tumors.

RESULTS

We report a mean of 11.6 mutations/Mb with a characteristic smoking signature along with mutations in TP53 (65%), CDKN2A (20%), NFE2L2 (20%), FAT1 (15%), KMT2C (15%), LRP1B (15%), FGFR1 (14%), PTEN (10%) and PREX2 (5%) among lung squamous cell carcinoma patients of Indian descent. In addition, therapeutically relevant EGFR mutations occur in 5.8% patients, significantly higher than as reported among Caucasians. In overall, our data suggests 13.5% lung squamous patients harboring druggable mutations have lower median overall survival, and 19% patients with a mutation in at least one gene, known to be associated with cancer, result in significantly shorter median overall survival compared to those without mutations.

CONCLUSIONS

We present the first comprehensive landscape of genetic alterations underlying Indian lung squamous cell carcinoma patients and identify EGFR, PIK3CA, KRAS and FGFR1 as potentially important therapeutic and prognostic target.

Metabolic rewiring in drug resistant cells exhibit higher OXPHOS and fatty acids as preferred major source to cellular energetics

Alteration in metabolic repertoire is associated with resistance phenotype. Although a common phenotype, not much efforts have been undertaken to design effective strategies to target the metabolic drift in cancerous cells with drug resistant properties. Here, we identified that drug resistant AML cell line HL-60/MX2 did not follow classical Warburg effect, instead these cells exhibited drastically low levels of aerobic glycolysis. Biochemical analysis confirmed reduced glucose consumption and lactic acid production by resistant population with no differences in glutamine consumption. Raman spectroscopy revealed increased lipid and cytochrome content in resistant cells which were also visualized as lipid droplets by Raman mapping, electron microscopy and lipid specific staining. Gene set enrichment analysis data from sensitive and resistant cell lines revealed significant enrichment of lipid metabolic pathways in HL-60/MX2 cells. Further, HL-60/MX2 possessed higher mitochondrial activity and increased OXPHOS suggesting the role of fatty acid metabolism as energy source which was confirmed by increased rate of fatty acid oxidation. Accordingly, OXPHOS inhibitor increased sensitivity of resistant cells to chemotherapeutic drug and fatty acid oxidation inhibitor Etomoxir reduced colony formation ability of resistant cells demonstrating the requirement of fatty acid metabolism and dependency on OXPHOS by resistant leukemic cells for survival and tumorigenicity.

Merkel Cell Polyomavirus in Merkel Cell Carcinoma: Integration Sites and Involvement of the KMT2D Tumor Suppressor Gene

Merkel cell carcinoma (MCC) is an uncommon, lethal cancer of the skin caused by either Merkel cell polyomavirus (MCPyV) or UV-linked mutations. MCPyV is found integrated into MCC tumor genomes, accompanied by truncation mutations that render the MCPyV large T antigen replication incompetent. We used the open access HPV Detector/Cancer-virus Detector tool to determine MCPyV integration sites in whole-exome sequencing data from five MCC cases, thereby adding to the limited published MCPyV integration site junction data. We also systematically reviewed published data on integration for MCPyV in the human genome, presenting a collation of 123 MCC cases and their linked chromosomal sites. We confirmed that there were no highly recurrent specific sites of integration. We found that chromosome 5 was most frequently involved in MCPyV integration and that integration sites were significantly enriched for genes with binding sites for oncogenic transcription factors such as LEF1 and ZEB1, suggesting the possibility of increased open chromatin in these gene sets. Additionally, in one case we found, for the first time, integration involving the tumor suppressor gene KMT2D, adding to previous reports of rare MCPyV integration into host tumor suppressor genes in MCC.

Chromatin protein PC4 is downregulated in breast cancer to promote disease progression: Implications of miR-29a

The human transcriptional coactivator PC4 has numerous roles to play in the cell. Other than its transcriptional coactivation function, it facilitates chromatin organization, DNA damage repair, viral DNA replication, etc. Although it was found to be an essential protein in vivo, the importance of this multifunctional protein in the regulation of different cellular pathways has not been investigated in details, particularly in oncogenesis. In this study, PC4 downregulation was observed in a significant proportion of mammary tissues obtained from Breast cancer patient samples as well as in a subset of highly invasive and metastatic Breast cancer patient-derived cell lines. We have identified a miRNA, miR-29a which potentially reduce the expression of PC4 both in RNA and protein level. This miR-29a was found to be indeed overexpressed in a substantial number of Breast cancer patient samples and cell lines as well, suggesting one of the key mechanisms of PC4 downregulation. Stable Knockdown of PC4 in MCF7 cells induced its migratory as well as invasive properties. Furthermore, in an orthotopic breast cancer mice model system; we have shown that reduced expression of PC4 enhances the tumorigenic potential substantially. Absence of PC4 led to the upregulation of several genes involved in Epithelial to Mesenchymal Transition (EMT), indicating the possible mechanism of uniform tumour progression in the orthotropic mice. Collectively these data establish the role of PC4 in tumour suppression.

Merkel cell polyomavirus is implicated in a subset of Merkel cell carcinomas, in the Indian subcontinent

Merkel cell carcinoma is a rare, lethal cancer histopathologically composed of cells showing similarity with mechanoreceptor Merkel cells. Merkel cell tumors manifest in two distinct forms. While a virus called Merkel cell polyomavirus is involved in the pathogenesis of one form of Merkel tumors, the other is driven by ultraviolet (UV)-linked mutations. In this study we investigated 18 cases, from the Indian population, of Merkel cell carcinoma for immunohistochemical (IHC) expression of Merkel cell polyomavirus (MCV) T antigen, including 12 cases tested by PCR, to identify viral etiopathology. We tested the tumors with two sensitive antibodies (CM2B4 and Ab3), targeting the viral large T antigen protein and with PCR primers targeting the N terminus of T antigen. Overall, we observed 38.8% (7/18) tumors displaying positive IHC expression of Merkel cell polyomavirus T antigen and 25% (3/12) tumors showing positive results, by both, immunohistochemistry and PCR. This constitutes the first report from India showing implication of MCV in Merkel cell carcinomas. Moreover, this is one of the larger series of Merkel cell carcinomas, tested for MCV, by both immunohistochemistry and PCR, in this part of the world. These results further indicate that a slightly more number of such cases in India are likely to be caused by UV-linked damage, as opposed to Merkel cell polyomavirus mediated tumorigenesis, which is definitely implicated in a subset of cases.

ERBB2 and KRAS alterations mediate response to EGFR inhibitors in early stage gallbladder cancer

The uncommonness of gallbladder cancer in the developed world has contributed to the generally poor understanding of the disease. Our integrated analysis of whole exome sequencing, copy number alterations, immunohistochemical, and phospho-proteome array profiling indicates ERBB2 alterations in 40% early-stage rare gallbladder tumors, among an ethnically distinct population not studied before, that occurs through overexpression in 24% (n = 25) and recurrent mutations in 14% tumors (n = 44); along with co-occurring KRAS mutation in 7% tumors (n = 44). We demonstrate that ERBB2 heterodimerizes with EGFR to constitutively activate the ErbB signaling pathway in gallbladder cells. Consistent with this, treatment with ERBB2-specific, EGFR-specific shRNA or with a covalent EGFR family inhibitor Afatinib inhibits tumor-associated characteristics of the gallbladder cancer cells. Furthermore, we observe an in vivo reduction in tumor size of gallbladder xenografts in response to Afatinib is paralleled by a reduction in the amounts of phospho-ERK, in tumors harboring KRAS (G13D) mutation but not in KRAS (G12V) mutation, supporting an essential role of the ErbB pathway. In overall, besides implicating ERBB2 as an important therapeutic target under neo-adjuvant or adjuvant settings, we present the first evidence that the presence of KRAS mutations may preclude gallbladder cancer patients to respond to anti-EGFR treatment, similar to a clinical algorithm commonly practiced to opt for anti-EGFR treatment in colorectal cancer.

Up-regulation of the kinase gene SGK1 by progesterone activates the AP-1-NDRG1 axis in both PR-positive and -negative breast cancer cells

Preoperative progesterone intervention has been shown to confer a survival benefit to breast cancer patients independently of their progesterone receptor (PR) status. This observation raises the question how progesterone affects the outcome of PR-negative cancer. Here, using microarray and RNA-Seq-based gene expression profiling and ChIP-Seq analyses of breast cancer cells, we observed that the serum- and glucocorticoid-regulated kinase gene (SGK1) and the tumor metastasis-suppressor gene N-Myc downstream regulated gene 1 (NDRG1) are up-regulated and that the microRNAs miR-29a and miR-101-1 targeting the 3'-UTR of SGK1 are down-regulated in response to progesterone. We further demonstrate a dual-phase transcriptional and post-transcriptional regulation of SGK1 in response to progesterone, leading to an up-regulation of NDRG1 that is mediated by a set of genes regulated by the transcription factor AP-1. We found that NDRG1, in turn, inactivates a set of kinases, impeding the invasion and migration of breast cancer cells. In summary, we propose a model for the mode of action of progesterone in breast cancer. This model helps decipher the molecular basis of observations in a randomized clinical trial of the effect of progesterone on breast cancer and has therefore the potential to improve the prognosis of breast cancer patients receiving preoperative progesterone treatment.

CytoPred: 7-gene pair metric for AML cytogenetic risk prediction

Cytogenetic-based subjective prognostication of acute myeloid leukemia (AML) patients is a cumbersome process. Top scoring pair (TSP)-based decision tree using a robust analytical algorithm with statistical rigor offers a promising alternative. We describe CytoPred as a 7-gene pair signature based on the analysis of 2547 AML patient sample gene expression data using a modified TSP algorithm to estimate cytogenetic risk. The essential modification in TSP that helped computational encumbrance includes the filtration of gene pairs above random weighted guessers as well as sampling the gene pairs from the original gene pair pool to reduce overfitting issue. The CytoPred classifies AML cohort into clinically relevant `good' and `Int_poor' prognosis groups with distinct survival differences. The 7-gene pair was derived using 1248 AML patient samples in training set and 675 samples used for internal testing of the algorithm. The finest classifier 7-gene pair was picked from an initial pool size of 6.1 × 107 gene pairs that generated 57 687 decision trees. Further, for unbiased evaluation of CytoPred performance, we did an independent validation in 624 AML patient cohort. The CytoPred well qualifies the cutoffs for diagnostic application with 98.27% sensitivity and 99.27% specificity to predictive value in Int_poor class while 97.09% sensitivity and 91.74% specificity to predictive value for good class. Furthermore, CytoPred predicts almost identical survival probabilities like cytogenetics and its performance is not much influenced by various recurrent mutations as well as individual French-American-British (FAB) subtypes. In summary, we present a robust 7-gene pair-based metric to clinically prognosticate AML patients.

miR-129-2 mediates down-regulation of progesterone receptor in response to progesterone in breast cancer cells

OBJECTIVE

Hormonal therapy is an important component of first line of treatment for breast cancer. Response to hormonal therapy is influenced by the progesterone receptor (PR)-status of breast cancer patients. However as an early effect, exposure to progesterone decreases expression of PR in breast cancer cells. An understanding of the mechanism underlying down-regulation of PR could help improve response to hormonal therapy.

METHODS

We performed small RNA sequencing of breast cancer cells for identification of microRNAs targeting PR in response to progesterone treatment. Biochemical approaches were used to validate the findings in breast cancer cells.

RESULTS

Analysis of small RNA sequencing of four breast cancer cell lines treated with progesterone revealed an up-regulation of miR-129-2 independent of the PR status of the cells. We show that miR-129-2 targets 3'UTR of PR to down-regulate its expression. Furthermore, inhibition of miR-129-2 expression rescues the down-regulation of PR in breast cancer cells. Also, the expression levels of miR-129-2 was observed to be elevated in patients with low expression of PR in the TCGA cohort (n = 359).

CONCLUSION

miR-129-2 mediates down-regulation of PR in breast cancer cells in response to progesterone, while anti-miR-129-2 could potentiate PR expression levels among patients with inadequate PR levels. Thus, modulation of activity of miR-129-2 could stabilize PR expression and potentially improve response to hormonal therapy under adjuvant or neo-adjuvant settings.

Genomic characterization of tobacco/nut chewing HPV-negative early stage tongue tumors identify MMP10 asa candidate to predict metastases

OBJECTIVES

Nodal metastases status among early stage tongue squamous cell cancer patients plays a decisive role in the choice of treatment, wherein about 70% patients can be spared from surgery with an accurate prediction of negative pathological lymph node status. This underscores an unmet need for prognostic biomarkers to stratify the patients who are likely to develop metastases.

MATERIALS AND METHODS

We performed high throughput sequencing of fifty four samples derived from HPV negative early stage tongue cancer patients habitual of chewing betel nuts, areca nuts, lime or tobacco using whole exome (n=47) and transcriptome (n=17) sequencing that were analyzed using in-house computational tools. Additionally, gene expression meta-analyses were carried out for 253 tongue cancer samples. The candidate genes were validated using qPCR and immuno-histochemical analysis in an extended set of 50 early primary tongue cancer samples.

RESULTS AND CONCLUSION

Somatic analysis revealed a classical tobacco mutational signature C:G>A:T transversion in 53% patients that were mutated in TP53, NOTCH1, CDKN2A, HRAS, USP6, PIK3CA, CASP8, FAT1, APC, and JAK1. Similarly, significant gains at genomic locus 11q13.3 (CCND1, FGF19, ORAOV1, FADD), 5p15.33 (SHANK2, MMP16, TERT), and 8q24.3 (BOP1); and, losses at 5q22.2 (APC), 6q25.3 (GTF2H2) and 5q13.2 (SMN1) were observed in these samples. Furthermore, an integrated gene-expression analysis of 253 tongue tumors suggested an upregulation of metastases-related pathways and over-expression of MMP10 in 48% tumors that may be crucial to predict nodal metastases in early tongue cancer patients. In overall, we present the first descriptive portrait of somatic alterations underlying the genome of tobacco/nut chewing HPV-negative early tongue cancer, and identify MMP10 asa potential prognostic biomarker to stratify those likely to develop metastases.

Abstract 4397: Deciphering the diversity of somatic alterations and Salmonella infection in gallbladder cancer by whole exome sequencing

Background: Gallbladder cancer is fifth most common cancers among the gastrointestinal cancers with the majority of patients presented at an advanced state of the disease. In India, gallbladder cancer is a major problem in the northern part of the country with its highest incidence of 22/1,00,000 women and risk factors ranging from gallstones, female gender, ethnicity, Salmonella infections and genetic alterations. Despite its high incidence in our country, there is only few candidate gene based studies and systematic genome-wide studies are far in dismal. Hence there is an unmet need to understand the genomic landscape of Indian gallbladder cancer genome.

Materials and methods: We interrogated the coding region of gallbladder cancer genome of 27 samples (10 paired and 7 unpaired tumors) using whole exome sequencing at an average coverage of 100X and above. Further, we validated our findings from whole exome sequencing in an extended cohort of 27 FFPE (Formalin fixed paraffin embedded) samples using other sequencing technologies. In addition we used a computational subtraction tool to identify Salmonella DNA sequences in the whole exome sequencing data.

Results: We identified 383 somatic alterations across 17 tumors, which includes an average 112 synonymous, 245 missense, 8 nonsense, 8 indels and 8 splice site changes. The average mutation rate considering the paired tumors is about 14 mutations/Mb. We found recurrent alterations in TP53, CTNNB1, SF3B1, ATM, AKAP11 and other genes by exome sequencing analysis. In addition, we examined our exome sequencing data for identifying Salmonella sequences as well as presence of 143 HPV types using computational subtraction based on HPVDetector. In our analysis we found association of typhoidal Salmonella strains in 11 of 26 gall bladder cancer samples and non-typhoidal Salmonella species in 12 of 26 sample, 6 samples were co-infected with both. Moreover, we observed co-occurrence of TP53 alterations in 4 of 16 Salmonella positive samples while we did not observe TP53 alterations in Salmonella negative samples.

Conclusion: Taken together, we present the landscape of somatic alterations in Indian gallbladder cancer genome and identification of non-typhoidal Salmonella species along with co- occurrence of TP53 alterations could aid in the treatment of gallbladder cancer.

Note: This abstract was not presented at the meeting.

Citation Format: Prajish Sundaram Iyer, Nilesh laxman Gardi, Malika Ranjan, Bikram Sahoo, Pratik Chandrani, Pawan Upadhyay, Mukta R. Ramadwar, Shailesh V. Shrikhande, Amit Dutt. Deciphering the diversity of somatic alterations and Salmonella infection in gallbladder cancer by whole exome sequencing [abstract]. In: Proceedings of the American Association for Cancer Research Annual Meeting 2017; 2017 Apr 1-5; Washington, DC. Philadelphia (PA): AACR; Cancer Res 2017;77(13 Suppl):Abstract nr 4397. doi:10.1158/1538-7445.AM2017-4397