Search

Motivation A major challenge in analyzing cancer patient transcriptomes is that the tumors are inherently heterogeneous and evolving. We analyzed 214 bulk RNA samples of a longitudinal, prospective ovarian cancer cohort and found that the sample composition changes systematically due to chemotherapy and between the anatomical sites, preventing direct comparison of treatment-naive and treated samples. Results To overcome this, we developed PRISM, a latent statistical framework to simultaneously extract the sample composition and cell-type-specific whole-transcriptome profiles adapted to each individual sample. Our results indicate that the PRISM-derived composition-free transcriptomic profiles and signatures derived from them predict the patient response better than the composite raw bulk data. We validated our findings in independent ovarian cancer and melanoma cohorts, and verified that PRISM accurately estimates the composition and cell-type-specific expression through whole-genome sequencing and RNA in situ hybridization experiments. Availabilityand implementation https://bitbucket.org/anthakki/prism. Supplementary information Supplementary data are available at Bioinformatics online.

While host immune response is likely to be important for the prognosis of gastric cancer patients, detailed information on the T lymphocyte infiltration in different gastric cancer subtypes is lacking. Here, we studied the presence of CD3, CD8, and FOXP3 (Forkhead box p3) expressing T lymphocytes in a retrospective cohort of 190 intestinal gastric and gastroesophageal adenocarcinomas. The cancers represented four distinct molecular subtypes: Epstein-Barr virus–positive (EBV+), mismatch-repair-deficient (MMR-D), aberrant TP53, and the “other” subtype. The absolute numbers of CD3+, CD8+, and FOXP3+ T lymphocytes were analyzed in relation with these molecular subtypes and selected clinicopathological parameters. Overall, there was a large variation in the amount of infiltrating T lymphocyte in all molecular subtypes. Among the subtypes, EBV+ cancers differed from the other subtypes in increased lymphocyte infiltration and high CD8+/FOXP3+ ratio. While the TP53 aberrant subtype did not differ in the absolute amount of T lymphocyte, the ratio of CD8+/FOXP3+ and CD3+/FOXP3+ cells was highest in this subtype, possibly reflecting immunosuppression associated with genomic instability. Increased CD3+ and CD8+ T lymphocyte infiltrates were associated with better survival, and remained as independent prognostic factors in a multivariate analysis. This study is the first to investigate lymphocytic infiltration within four molecular subtypes of intestinal-type gastric cancer in a European cohort. The results provide an important addition to the current knowledge of T lymphocyte–dependent immune response in gastric cancer and its prognostic significance.

Chemotherapy resistance is a critical contributor to cancer mortality and thus an urgent unmet challenge in oncology. To characterize chemotherapy resistance processes in high-grade serous ovarian cancer, we prospectively collected tissue samples before and after chemotherapy and analyzed their transcriptomic profiles at a single-cell resolution. After removing patient-specific signals by a novel analysis approach, PRIMUS, we found a consistent increase in stress-associated cell state during chemotherapy, which was validated by RNA in situ hybridization and bulk RNA sequencing. The stress-associated state exists before chemotherapy, is subclonally enriched during the treatment, and associates with poor progression-free survival. Co-occurrence with an inflammatory cancer–associated fibroblast subtype in tumors implies that chemotherapy is associated with stress response in both cancer cells and stroma, driving a paracrine feed-forward loop. In summary, we have found a resistant state that integrates stromal signaling and subclonal evolution and offers targets to overcome chemotherapy resistance.

Background Gastric cancer (GC) is the third most common cause of cancer death. Intestinal type GC is a molecularly diverse disease. Formins control cytoskeletal processes and have been implicated in the progression of many cancers. Their clinical significance in GC remains unclear. Here, we characterize the expression of formin proteins FHOD1 and FMNL1 in intestinal GC tissue samples and investigate their association with clinical parameters, GC molecular subtypes and intratumoral T lymphocytes. Methods The prognostic significance of FHOD1 and FMNL1 mRNA expression was studied with Kaplan–Meier analyses in an online database. The expression of FHOD1 and FMNL1 proteins was characterized in GC cells, and in non-neoplastic and malignant tissues utilizing tumor microarrays of intestinal GC representing different molecular subtypes. FHOD1 and FMNL1 expression was correlated with clinical parameters, molecular features and T lymphocyte infiltration. Immunohistochemical expression of neither formin correlated with survival. Results Kaplan–Meier analysis associated high FHOD1 and FMNL1 mRNA expression with reduced overall survival (OS). Characterization of FHOD1 and FMNL1 in GC cells showed cytoplasmic expression along the actin filaments. Similar pattern was recapitulated in GC tissue samples. Elevated FMNL1 was associated with larger tumor size and higher disease stage. Downregulation of FHOD1 associated with TP53-mutated GC tumors. Tumor cell FHOD1 expression strongly correlated with high numbers of tumor-infiltrating CD8 + lymphocytes. Conclusions FHOD1 and FMNL1 proteins are expressed in the tumor cells of intestinal GC and significantly associate with clinical parameters without direct prognostic significance. FHOD1 correlates with high intratumoral CD8 + T lymphocyte infiltration in this cohort.

Aim: Colorectal cancer (CRC) is one of the leading cancers in the world. Even though its mortality and pathophysiology are well documented in the US and the European countries, it is seldom studied in North African population. Recent studies have shown link of HER2 overexpression in oesophageal and gastric cancers. The aim of this study is to assess the HER2 protein and mRNA expression and its correlation with tumor pathogenesis in Libyan CRC patients. Methodology: A total of 17 FFPE tissue blocks were collected from patients with primary CRC. The HER2 protein expression was assessed by immunohistochemistry and the mRNA expression was assessed using qRT-PCR. Survival analysis of the role of HER2 overexpression on rectal adenocarcinoma was carried out on additional 165 patients. Results: From the CRC cohort, adenocarcinoma was found to be more frequent accounting for 88.2%, and 11.8% for mucinous adenocarcinomas. Almost 47% of the cases were positive for HER2 (score ≥ 2+) and about 50% adenocarcinoma cases with tumor grade II were positive for HER2. Moreover, 57.4% adenocarcinoma patients with grade-II tumor had undergone right hemicolectomy. Furthermore, significant correlation (p = 0.03) between the HER2 mRNA expression with the tumor grade was observed. In addition, poor overall all survival was observed with high HER2 expression in rectum adenocarcinoma. Conclusion: To our knowledge, this is the first study that HER2 overexpression correlates with more aggressive colorectal cancer in North African population. Our study shows that HER2 overexpression associates with right colon surgeries. Also, the correlation of mRNA and protein expression could warrant the implementation of a nationwide screening program for HER2 positivity in CRC patients. Taken together, stratifying patients according to HER2 expression can help in the diagnosis and prognosis of CRC patients from North African origin.

Chemotherapy resistance is the greatest contributor to cancer mortality and the most urgent unmet challenge in oncology. In order to reveal transcriptomics changes due to platinum-based chemotherapy we analyzed single-cell RNA-seq data from fresh tissue samples taken at the time of diagnosis and after neoadjuvant chemotherapy from 11 high-grade serous ovarian cancer (HGSOC) patients. With a novel clustering method accounting for patient-specific variability and technical confounders, we identified 12 clusters. Of these, a stress-related transcriptional profile was enriched during chemotherapy and associated significantly to poor progression-free survival (PFS) and disease-free interval (DFI) in deconvoluted bulk RNA-seq data analysis of treatment-naive samples in TCGA cohort. Pan-cancer cell line analysis suggests that patients with high stress-related transcriptional profile may benefit from MEK1/2 inhibitors instead of platinum. Further, high proportion of stromal components and high interaction score between tumor and stromal suggest the tumor cells with high-stress profile actively interact with and potentially recruit stromal cells to their microenvironment already prior to chemotherapy, potentially facilitating protection from chemotherapeutic treatments. In summary, we have identified a stress-related transcriptional profile, which is present at the time of diagnosis, enriched during platinum treatment, independent predictor for poor PFS and DFI, and, based on in vitro data, targetable with MEK1/2 inhibitors.Translational relevanceWe discovered a stress-related transcriptional profile that is significantly enriched in fresh tissue samples after chemotherapy and is significantly associated with poor progression-free survival in an independent patient cohort. The survival association is independent of age, tumor purity or BRCAness. Therefore, this chemotherapy resistance associated profile is intrinsic and could thus be targeted already in treatment-naive patients. The translation potential of the stress-related transcriptomics profile was further supported by pan-cancer cell line analysis that showed that cell lines with high stress-related transcriptional profile are not affected by platinum, corroborating our results, whereas they were more sensitive to MEK1/2 inhibitors. Taken together, the stress-related transcriptional profile, quantifiable with a set of 35 marker genes, provides a basis for improved prediction of platinum response as well as novel avenues to treat this patient group more effectively.

A major challenge in analyzing cancer patient transcriptomes is that the tumors are inherently heterogeneous and evolving. We analyzed 214 bulk RNA samples of a longitudinal, prospective ovarian cancer cohort and found that the sample composition changes systematically due to chemotherapy and between the anatomical sites, preventing direct comparison of treatment-naive and treated samples. To overcome this, we developed PRISM, a latent statistical framework to simultaneously extract the sample composition and cell type specific whole-transcriptome profiles adapted to each individual sample. Our results indicate that the PRISM-derived composition-free transcriptomic profiles and signatures derived from them predict the patient response better than the composite raw bulk data. We validated our findings in independent ovarian cancer and melanoma cohorts, and verified that PRISM accurately estimates the composition and cell type specific expression through whole-genome sequencing and RNA in situ hybridization experiments. PRISM is freely available with full source code and documentation.

Tumor heterogeneity in high-grade serous ovarian cancer (HGSOC) stems from highly variable copy number changes that lack shared actionable drug targets and contribute to inevitable recurrence and death. To unveil the vast heterogeneity in HGSOC, and understand the molecular changes that occur during chemotherapy, we analyzed single-cell transcriptomes of 49 dissociated HGSOC tumor or ascites specimens, along with six normal postmenopausal fimbria/tuba specimens as controls. This set contains 13 pairs of samples from patients who underwent neoadjuvant chemotherapy, taken before treatment (primary), and after three cycles of standard carboplatin-paclitaxel chemotherapy (interval). Clustering single-cell transcriptomes based on acknowledged cell type markers revealed epithelial cancer cells, as well as diverse stromal and immune cell types. The cell type proportions displayed high variability across HGSOC tumors, intriguingly with lowest tumor purity in patients with short platinum-free interval. To identify the impact of chemotherapy on transcriptional programs, we extracted shared transcriptomic differences in the cancer cells of seven primary-interval pairs. Regardless of the clinical response, three cycles of chemotherapy induced an oxidative stress response program that inhibits apoptosis, while inducing genes involved in MHC II antigen presentation to T helper cells. Next, to study how cell states change during chemotherapy, we ordered the cells in pseudotime based on their transcriptomic similarity, revealing a mixture of primary and interval cells in all but two primary-interval pairs. Furthermore, based on inferred copy number analysis, we estimate that genomic differences contribute variably, between 20%, to the overall transcriptomic differences between the primary and interval pairs. These results suggest that pre- and post-chemotherapy states are not completely separate and that genomic evolution has a highly variable effect on the state differences, providing support for an intrinsic chemoresistance model with genomic and nongenomic components. To further the understanding of intrinsic resistance, we next focused on patients’ clinical response mechanisms in HGSOC cancer cells. Genes encoding well-known serum markers of HGSOC, WFDC2 and MUC16, were expressed at higher levels in primary specimens of patients with incomplete response, whereas patients with complete response displayed higher expression of genes related to chromatin repression/heterochromatin. Finally, we focused on tumor resectability to understand the transcriptomic patterns underlying HGSOC metastasis. Already prior to chemotherapy, patients with residual disease showed induced FOXO3-driven regulon linked to quiescence, stress resistance, and stemness, whereas patients with no residual disease express increased level of E2F1/DP1-driven regulon linked to BRCAness, high proliferation, and high metabolic rate. Citation Format: Erdogan Pekcan Erkan, Jun Dai, Kaiyang Zhang, Katja Kaipio, Tarja Lamminen, Naziha Mansuri, Lasse Suominen, Amjad Alkodsi, Kaisa Huhtinen, Sakari Hietanen, Johanna Hynninen, Seija Grénman, Olli Carpén, Sampsa Hautaniemi, Anna Vähärautio. Singling out tumor heterogeneity and chemoresistance in high-grade serous ovarian cancer [abstract]. In: Proceedings of the AACR Special Conference on Advances in Ovarian Cancer Research; 2019 Sep 13-16, 2019; Atlanta, GA. Philadelphia (PA): AACR; Clin Cancer Res 2020;26(13_Suppl):Abstract nr A40.

Introduction: To identify genes and pathways that drive platinum-taxane resistance in high-grade serous ovarian cancer (HGSOC), we analyzed 108 samples subjected to RNA-seq before and after chemotherapy. Experimental Procedure: We collected prospectively 108 tissue samples from various anatomic locations (omentum, ovary, mesothelium, peritoneum) for 53 HGSOC patients with high enough tumor purity (>25%) and follow-up information to calculate platinum free interval (PFI). These samples were subjected to RNA-seq to obtain gene expression data for known protein coding genes. Additionally, data from 16,826 samples (cells) with single-cell RNA-seq were used in the decomposition analysis to overcome cell heterogeneity in the samples. For validation we used 306 The Cancer Genome Atlas (TCGA) HGSOC samples that were analyzed with the decompose algorithm. Differential expressions for selected resistance associated genes were validated with immunohistochemistry (IHC). Results: We examined firstly whether gene expressions values are systematically influenced by anatomic location of the samples and chemotherapy. Our results show that bias due to anatomic locations can be controlled by relatively simple computational adjusting. However, chemotherapy-treated samples have significantly lower tumor purity than treatment-naive samples and samples have varying cell compositions. We developed a computational decomposition method to correct these biases. Secondly, we compared gene expression values of cancer cells between 1) well and poor responding patients, and 2) before and after chemotherapy followed by pathway analysis. We identified several genes that were significantly upregulated in the poor-responder group, such as FOXK1 and FBXO32, which also had significant survival association in the TCGA cohort. Our results highlight genes whose expression values are significantly enriched after chemotherapy, such as BTG2, CITED2, CTGF, FOS, DUSP1, and EGR1, which have been shown to promote resistance to platinum or paclitaxel in HGSOC or other cancers. In addition to these known chemoresistant-associated genes, we identified several hitherto unknown genes and validated them with IHC. Pathway analysis revealed actionable targets in the B-cell receptor and MAPK signaling pathways, which were the most significant enriched signaling pathways after chemotherapy. Conclusions: We have shown that successful comparison of gene expression values before and after chemotherapy requires heavy computational correction of the strong and systematic bias caused by chemotherapy and tumor heterogeneity. To overcome these biases, we have developed a novel decomposition algorithm. This method allows to obtain reliable data for the subsequent analyses. The gene expression analyses herein identified several genes and pathways that play central roles in chemoresistance and are rational targets for overcoming platinum-taxane resistance in HGSOC. Citation Format: Sanaz Jamalzadeh, Kaiyang Zhang, Antti Häkkinen, Kaisa Huhtinen, Johanna Hynninen, Noora Andersson, Naziha Mansuri, Olli Carpén, Erkan Erdogan, Jun Dai, Anna Vähärautio, Sakari Hietanen, Jaana Oikkonen, Rainer Lehtonen, Sampsa Hautaniemi. RNA-seq analysis of high-grade serous ovarian cancer patients before and after chemotherapy reveals chemoresistance-associated genes and pathways [abstract]. In: Proceedings of the AACR Special Conference on Advances in Ovarian Cancer Research; 2019 Sep 13-16, 2019; Atlanta, GA. Philadelphia (PA): AACR; Clin Cancer Res 2020;26(13_Suppl):Abstract nr A59.

Poor chemotherapy response remains a major treatment challenge for high-grade serous ovarian cancer (HGSC). Cancer stem cells are the major contributors to relapse and treatment failure as they can survive conventional therapy. Our objectives were to characterise stemness features in primary patient-derived cell lines, correlate stemness markers with clinical outcome and test the response of our cells to both conventional and exploratory drugs. Tissue and ascites samples, treatment-naive and/or after neoadjuvant chemotherapy, were prospectively collected. Primary cancer cells, cultured under conditions favouring either adherent or spheroid growth, were tested for stemness markers; the same markers were analysed in tissue and correlated with chemotherapy response and survival. Drug sensitivity and resistance testing was performed with 306 oncology compounds. Spheroid growth condition HGSC cells showed increased stemness marker expression (including aldehyde dehydrogenase isoform I; ALDH1A1) as compared with adherent growth condition cells, and increased resistance to platinum and taxane. A set of eight stemness markers separated treatment-naive tumours into two clusters and identified a distinct subgroup of HGSC with enriched stemness features. Expression of ALDH1A1, but not most other stemness markers, was increased after neoadjuvant chemotherapy and its expression in treatment-naive tumours correlated with chemoresistance and reduced survival. In drug sensitivity and resistance testing, five compounds, including two PI3K-mTOR inhibitors, demonstrated significant activity in both cell culture conditions. Thirteen compounds, including EGFR, PI3K-mTOR and aurora kinase inhibitors, were more toxic to spheroid cells than adherent cells. Our results identify stemness markers in HGSC that are associated with a decreased response to conventional chemotherapy and reduced survival if expressed by treatment-naive tumours. EGFR, mTOR-PI3K and aurora kinase inhibitors are candidates for targeting this cell population. (c) 2019 Pathological Society of Great Britain and Ireland. Published by John Wiley & Sons, Ltd.

Gastric cancer is traditionally divided into intestinal and diffuse histological subtypes, but recent molecular analyses have led to novel classification proposals based on genomic alterations. While the intestinal- and diffuse-type tumours are distinguishable from each other at the molecular level, intestinal-type tumours have more diverse molecular profile. The technology required for comprehensive molecular analysis is expensive and not applicable for routine clinical diagnostics. In this study, we have used immunohistochemistry and in situ hybridisation in molecular classification of gastric adenocarcinomas with an emphasis on the intestinal subtype. A tissue microarray consisting of 244 gastric adenocarcinomas was constructed, and the tumours were divided into four subgroups based on the presence of Epstein-Barr virus, TP53 aberrations and microsatellite instability. The intestinal- and diffuse-type tumours were separately examined. The distribution of EGFR and HER2 gene amplifications was studied in the intestinal-type tumours. Epstein-Barr virus positive intestinal-type tumours were more common in male patients (p = 0.035) and most often found in the gastric corpus (p = 0.011). The majority of the intestinal-type tumours with TP53 aberrations were proximally located (p = 0.010). All tumours with microsatellite instability showed intestinal-type histology (p = 0.017) and were associated with increased overall survival both in the univariate (p = 0.040) and multivariate analysis (p = 0.015). In conclusion, this study shows that gastric adenocarcinomas can be classified into biologically and clinically different subgroups by using a simple method also applicable for clinical diagnostics.

Background: Among the patients with bladder cancer, a group is still at risk of disease recurrence, progression, and death from their cancer after curative treatment. Angiogenesis is a crucial pathogenic mechanism for this type of urothelial carcinoma and is a potential therapeutic target. Objectives: To quantify tumor angiogenesis in bladder cancer and determine whether it correlates with tumor stage and grade. Patients and methods: A series of 42 archival samples from carcinomas of the urinary bladder were graded, staged, and analyzed for microvessel density (MVD) by a double immunohistochemical technique using Factor VIII (FVIII) and CD31 antibodies. The correlation between MVD and histopathological grade and tumor stage was evaluated. Results: FVIII and CD31 immunoreactivity was observed in 100% of cases and more intensely with CD31. Significantly higher MVD was determined in invasive tumors than in superficial tumors (pB0.05). MVD increased with tumor grade and stage (pB0.05); MVD was not affected by age or sex of the patients. Conclusion: These data demonstrate that MVD in bladder carcinoma correlates with the tumor grade and stage. Quantification of tumor angiogenesis may allow selection of the type of treatment for bladder cancer patients. Keywords: angiogenesis; bladder carcinoma; Factor VIII; CD31; microvessel density (Published: 31 March 2011) Citation: Libyan J Med 2011, 6 : 6016 - DOI: 10.3402/ljm.v6i0.6016