Search

Endometrial carcinoma, the most common gynaecological cancer, develops from endometrial epithelium which is composed of secretory and ciliated cells. Pathologic classification is unreliable and there is a need for prognostic tools. We used single cell sequencing to study organoid model systems derived from normal endometrial endometrium to discover novel markers specific for endometrial ciliated or secretory cells. A marker of secretory cells (MPST) and several markers of ciliated cells (FAM92B, WDR16, and DYDC2) were validated by immunohistochemistry on organoids and tissue sections. We performed single cell sequencing on endometrial and ovarian tumours and found both secretory-like and ciliated-like tumour cells. We found that ciliated cell markers (DYDC2, CTH, FOXJ1, and p73) and the secretory cell marker MPST were expressed in endometrial tumours and positively correlated with disease-specific and overall survival of endometrial cancer patients. These findings suggest that expression of differentiation markers in tumours correlates with less aggressive disease, as would be expected for tumours that retain differentiation capacity, albeit cryptic in the case of ciliated cells. These markers could be used to improve the risk stratification of endometrial cancer patients, thereby improving their management. We further assessed whether consideration of MPST expression could refine the ProMiSE molecular classification system for endometrial tumours. We found that higher expression levels of MPST could be used to refine stratification of three of the four ProMiSE molecular subgroups, and that any level of MPST expression was able to significantly refine risk stratification of the copy number high subgroup which has the worst prognosis. Taken together, this shows that single cell sequencing of putative cells of origin has the potential to uncover novel biomarkers that could be used to guide management of cancers. © 2020 Pathological Society of Great Britain and Ireland. Published by John Wiley & Sons, Ltd.

High-grade serous ovarian cancer (HGSOC) is an archetypal cancer of genomic instability1–4 patterned by distinct mutational processes5,6, tumour heterogeneity7–9 and intraperitoneal spread7,8,10. Immunotherapies have had limited efficacy in HGSOC11–13, highlighting an unmet need to assess how mutational processes and the anatomical sites of tumour foci determine the immunological states of the tumour microenvironment. Here we carried out an integrative analysis of whole-genome sequencing, single-cell RNA sequencing, digital histopathology and multiplexed immunofluorescence of 160 tumour sites from 42 treatment-naive patients with HGSOC. Homologous recombination-deficient HRD-Dup (BRCA1 mutant-like) and HRD-Del (BRCA2 mutant-like) tumours harboured inflammatory signalling and ongoing immunoediting, reflected in loss of HLA diversity and tumour infiltration with highly differentiated dysfunctional CD8+ T cells. By contrast, foldback-inversion-bearing tumours exhibited elevated immunosuppressive TGFβ signalling and immune exclusion, with predominantly naive/stem-like and memory T cells. Phenotypic state associations were specific to anatomical sites, highlighting compositional, topological and functional differences between adnexal tumours and distal peritoneal foci. Our findings implicate anatomical sites and mutational processes as determinants of evolutionary phenotypic divergence and immune resistance mechanisms in HGSOC. Our study provides a multi-omic cellular phenotype data substrate from which to develop and interpret future personalized immunotherapeutic approaches and early detection research.

High-grade serous ovarian cancer (HGSOC) is an archetypal cancer of genomic instability patterned by distinct mutational processes, intratumoral heterogeneity and intraperitoneal spread. We investigated determinants of immune recognition and evasion in HGSOC to elucidate co- evolutionary processes underlying malignant progression and tumor immunity. Mutational processes and anatomic sites of tumor foci were key determinants of tumor microenvironment cellular phenotypes, inferred from whole genome sequencing, single-cell RNA sequencing, digital histopathology and multiplexed immunofluorescence of 160 tumor sites from 42 treatment-naive HGSOC patients. Homologous recombination-deficient (HRD)-Dup (BRCA1 mutant-like) and HRD- Del (BRCA2 mutant-like) tumors harbored increased neoantigen burden, inflammatory signaling and ongoing immunoediting, reflected in loss of HLA diversity and tumor infiltration with highly- differentiated dysfunctional CD8+ T cells. Foldback inversion (FBI, non-HRD) tumors exhibited elevated TGFβ signaling and immune exclusion, with predominantly naive/stem-like and memory T cells. Our findings implicate distinct immune resistance mechanisms across HGSOC subtypes which can inform future immunotherapeutic strategies.HIGHLIGHTSMulti-region, multi-modal profiling of malignant and immune cell phenotypes in ovarian cancerAnatomic site specificity is a determinant of cancer cell and intratumoral immune phenotypesTumor mutational processes impact mechanisms of immune control and immune evasionSpatial topology of HR-deficient tumors is defined by immune interactions absent from immune inert HR-proficient subtypes

Structural genome alterations are determinants of cancer ontogeny and therapeutic response. While bulk genome sequencing has enabled delineation of structural variation (SV) mutational processes which generate patterns of DNA damage, we have little understanding of how these processes lead to cell-to-cell variations which underlie selection and rates of accrual of different genomic lesions. We analysed 309 high grade serous ovarian and triple negative breast cancer genomes to determine their mutational processes, selecting 22 from which we sequenced >22,000 single cell whole genomes across a spectrum of mutational processes. We show that distinct patterns of cell-to-cell variation in aneuploidy, copy number alteration (CNA) and segment length occur in homologous recombination deficiency (HRD) and fold-back inversion (FBI) phenotypes. Widespread aneuploidy through induction of HRD throughBRCA1andBRCA2inactivation was mirrored by continuous whole genome duplication in HRD tumours, contrasted with early ploidy fixation in FBI. FBI tumours exhibited copy number distributions skewed towards gains, widespread clone-specific variation in amplitude of high-level amplifications, often impacting oncogenes, and break-point variability consistent with progressive genomic diversification, which we termed serriform structural variation (SSV). SSVs were consistent with a CNA-based molecular clock reflecting a continual and distributed process across clones within tumours. These observations reveal previously obscured genome plasticity and evolutionary properties with implications for cancer evolution, therapeutic targeting and response.

How cell-to-cell copy number alterations that underpin genomic instability

BackgroundSingle-cell RNA sequencing (scRNAseq) is a powerful tool for studying complex biological systems, such as tumour heterogeneity and tissue microenvironments. However, the sources of technical and biological variation in primary solid tumour tissues and patient-derived mouse xenografts for scRNAseq, are not well understood. Here, we used low temperature (6°C) protease and collagenase (37°C) to identify the transcriptional signatures associated with tissue dissociation across a diverse scRNAseq dataset comprising 128,481 cells from patient cancer tissues, patient-derived breast cancer xenografts and cancer cell lines.ResultsWe observe substantial variation in standard quality control (QC) metrics of cell viability across conditions and tissues. From FACS sorted populations gated for cell viability, we identify a sub-population of dead cells that would pass standard data filtering practices, and quantify the extent to which their transcriptomes differ from live cells. We identify a further subpopulation of transcriptomically “dying” cells that exhibit up-regulation of MHC class I transcripts, in contrast with live and fully dead cells. From the contrast between tissue protease dissociation at 37°C or 6°C, we observe that collagenase digestion results in a stress response. We derive a core gene set of 512 heat shock and stress response genes, includingFOSandJUN, induced by collagenase (37°C), which are minimized by dissociation with a cold active protease (6°C). While induction of these genes was highly conserved across all cell types, cell type-specific responses to collagenase digestion were observed in patient tissues. We observe that the yield of cancer and non-cancer cell types varies between tissues and dissociation methods.ConclusionsThe method and conditions of tumour dissociation influence cell yield and transcriptome state and are both tissue and cell type dependent. Interpretation of stress pathway expression differences in cancer single cell studies, including components of surface immune recognition such as MHC class I, may be especially confounded. We define a core set of 512 genes that can assist with identification of such effects in dissociated scRNA-seq experiments.

Single-cell RNA sequencing (scRNA-seq) has transformed biomedical research, enabling decomposition of complex tissues into disaggregated, functionally distinct cell types. For many applications, investigators wish to identify cell types with known marker genes. Typically, such cell type assignments are performed through unsupervised clustering followed by manual annotation based on these marker genes, or via "mapping" procedures to existing data. However, the manual interpretation required in the former case scales poorly to large datasets, which are also often prone to batch effects, while existing data for purified cell types must be available for the latter. Furthermore, unsupervised clustering can be error-prone, leading to under- and over- clustering of the cell types of interest. To overcome these issues we present CellAssign, a probabilistic model that leverages prior knowledge of cell type marker genes to annotate scRNA-seq data into pre-defined and de novo cell types. CellAssign automates the process of assigning cells in a highly scalable manner across large datasets while simultaneously controlling for batch and patient effects. We demonstrate the analytical advantages of CellAssign through extensive simulations and exemplify real-world utility to profile the spatial dynamics of high-grade serous ovarian cancer and the temporal dynamics of follicular lymphoma. Our analysis reveals subclonal malignant phenotypes and points towards an evolutionary interplay between immune and cancer cell populations with cancer cells escaping immune recognition.

Single-cell RNA sequencing (scRNA-seq) has enabled decomposition of complex tissues into functionally distinct cell types. Often, investigators wish to assign cells to cell types, performed through unsupervised clustering followed by manual annotation, or via “mapping” procedures to existing data. However, manual interpretation scales poorly to large datasets, mapping approaches require purified or pre-annotated data, and both are prone to batch effects. To overcome these issues we present CellAssign (www.github.com/irrationone/cellassign), a probabilistic model that leverages prior knowledge of cell type marker genes to annotate scRNA-seq data into pre-defined or de novo cell types. CellAssign automates the process of assigning cells in a highly scalable manner across large datasets while controlling for batch and sample effects. We demonstrate the advantages of CellAssign through extensive simulations and analysis of tumor microenvironment composition in high grade serous ovarian cancer and follicular lymphoma.

Endometrial epithelium gives rise to both endometrial and ovarian cancers (of clear-cell and endometrioid subtypes), the latter arising from ectopic endometrium (endometriosis). Endometrial epithelium comprises mainly secretory cells, with a minor ciliated cell population. Due to their scarcity, little is known about the biology or function of endometrial ciliated cells. To understand the biology of endometrial epithelium, and by extension the cancers that arise from it, organoids derived from normal endometrial tissue were cultured. Notch signaling inhibitors were used to induce ciliated cell differentiation. Through single-cell RNA sequencing, distinct secretory and ciliated cell populations were observed, with the ciliated cell population increasing with Notch signaling inhibition. Many novel markers of ciliated cells were observed, but no highly specific markers of secretory cell differentiation. A marker of secretory cells (MST) and several markers of ciliated cells (FAM92B, WDR16 and DYDC2) were validated by immunohistochemistry on organoids and tissue sections. In endometrial tumors, both MST and FAM92B exhibited diffuse staining and were markers of better prognosis. This suggests that tumors expressing differentiation markers have better prognosis, whether it is a marker of secretory or ciliated cells. Interestingly, a small number of endometrial tumors stained positive for DYDC2; however, these tumors exhibited a variable staining pattern with 25-50% tumor cells staining intensely, and the remaining tumor cells not staining at all. A similar variable staining pattern had been observed previously with CTH, another ciliated cell marker. Endometrial and ovarian tumor tissue microarrays were stained with DYDC2, CTH and two ciliated cell markers, FOXJ1 and p73. For all these markers, a subset of tumors displayed a variable staining pattern and for endometrial cancers, the variable staining was a good prognostic indicator. Single-cell sequencing of endometrial tumors has been able to capture these two populations of tumor cells. In ovarian tumors, only variable CTH staining was a significant prognostic indicator. Normal endometrial secretory cells are able to differentiate into ciliated cells, and the variable staining pattern suggests that a subset of tumors retains this ability, and these are clinically less aggressive. Using single-cell sequencing technology on normal tissues to guide development of prognostic markers and provide insight into the biology of the tumors arising from these tissues may be useful for many other tumor types. Citation Format: Dawn R. Cochrane, Kieran R. Campbell, Kendall Greening, Germain C. Ho, James Hopkins, Minh Bui, Vassilena Sharlandjieva, Daniel Lai, Maya DeGrood, Evan W. Gibbard, Samuel Leung, Angela S. Cheng, Jamie L.P. Lim, Samantha Neilson, David Farnell, Friedrich Kommoss, Jessica N. McAlpine, Sohrab P. Shah, David G. Huntsman. Single-cell RNA sequencing of normal endometrial organoids uncovers novel cell-type markers for prognostication of primary tumor samples [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 B09.

The telomerase reverse transcriptase (TERT) gene is highly expressed in stem cells and silenced upon differentiation. Cancer cells can attain immortality by activating TERT to maintain telomere length and telomerase activity, which is a crucial step of tumorigenesis. Two somatic mutations in the TERT promoter (C228T; C250T) have been identified as gain-of-function mutations that promote transcriptional activation of TERT in multiple cancers, such as melanoma and glioblastoma. A recent study investigating TERT promoter mutations in ovarian carcinomas found C228T and C250T mutations in 15.9% of clear cell carcinomas. However, it is unknown whether these mutations are frequent in other ovarian cancer subtypes, in particular, sex cordstromal tumors including adult granulosa cell tumors (AGCTs). We performed whole genome sequencing on ten AGCTs with matched normal blood and identified a TERT C228T promoter mutation in 50% of tumors. We found that AGCT with mutated TERT promoter have increased expression of TERT mRNA compared to those with wild-type TERT promoter. AGCT with TERT C228T mutation exhibited significantly longer telomeres compared to AGCT with TERT wild-type promoter. Extension cohort analysis using allelic discrimination revealed the TERT C228T mutation in 51 of 229 primary AGCTs (22%), 24 of 58 recurrent AGCTs (41%), and 1 of 22 other sex cord-stromal tumors (5%). There was a significant difference in overall survival between patients with TERT C228T promoter mutation in the primary tumors and those without it (p = 0.00253, log rank test). In seven AGCTs, we found the TERT C228T mutation present in recurrent tumors and absent in the corresponding primary tumor. Our data suggests that TERT C228T mutations may have an important role in progression of AGCT. Telomeres are conserved, repetitive (TTAGGG) DNA-protein complexes that are added to the ends of chromosomes by the enzyme telomerase to prevent DNA damage and maintain replicative potential. Telomere attrition during DNA replication induces genomic instability that can result in tumorigenesis. Telomerase consists of a catalytic protein subunit known as telomerase reverse transcriptase (TERT) and a functional RNA called telomerase RNA component (TERC). TERT is highly expressed in stem cells and is silenced upon differentiation in somatic cells. Most cancer cells attain proliferative immortality by upregulating the TERT gene to maintain telomere length and telomerase activity. The known mechanisms of telomerase activation include mutations in the TERT promoter, TERT gene amplification, CpG methylation at the TERT promoter, changes in alternative splicing of TERT pre-mRNA and upregulation of transcriptional activators. Approximately 90% of cancers express TERT, while the remaining 10-15% of cancers maintain their telomere length through a telomerase-independent method called alternative lengthening of telomeres. TERT promoter mutations were first reported in familial melanoma and subsequently in sporadic melanoma. There are two hot-spot TERT promoter mutations, C228T and C250T, each generates an identical 11 base pair sequence containing a consensus binding motif for ETS transcription factors, and functions as either a transcriptional activator or repressor to regulate telomerase expression. These two mutations are implicated in the activation of telomerase in other malignances such as central nervous system tumors, hepatocellular carcinomas, bladder cancers and thyroid cancers. A recent study on TERT promoter mutations in gynecological malignancies, including ovarian and uterine carcinomas, reported TERT hot-spot mutations in 15.9% of ovarian clear cell carcinomas. However, it is unknown whether TERT promoter mutations are frequent in sex cord-stromal tumors, including adult granulosa cell tumors (AGCTs). In this study, we evaluated the biological and clinical significance of TERT promoter mutations, specifically C228T, in total of 251 primary ovarian sex cord-stromal tumors.

SummaryHigh-grade serous ovarian cancer exhibits extensive intratumoral heterogeneity coupled with widespread intraperitoneal disease. Despite this, metastatic spread of tumor clones is non-random, implying the existence of local microenvironmental factors that shape tumor progression. We interrogated the molecular interface between tumor-infiltrating lymphocytes (TIL) and cancer cells in 143 samples from 21 patients using whole-genome sequencing, immunohistochemistry, histologic image analysis, gene expression profiling, and T- and B-cell receptor sequencing. We identify 3 immunologic response categories, which frequently co-exist within individual patients. Furthermore, epithelial CD8+ TIL were inversely associated with malignant cell diversity, evidenced by subclonal neoepitope elimination and spatial tracking between tumor and T-cell clones. Intersecting mutational signatures and immune analysis showed that foldback inversion genomic aberrations lead to worse outcomes even in the presence of cytotoxic TIL (n=433). Thus, regional variation in immune contexture mirrors the pattern of intraperitoneal malignant spread, provoking new perspectives for treatment of this challenging disease.

High-grade serous ovarian cancer (HGSC) exhibits extensive malignant clonal diversity with widespread but non-random patterns of disease dissemination. We investigated whether local immune microenvironment factors shape tumor progression properties at the interface of tumor-infiltrating lymphocytes (TILs) and cancer cells. Through multi-region study of 212 samples from 38 patients with whole-genome sequencing, immunohistochemistry, histologic image analysis, gene expression profiling, and T and B cell receptor sequencing, we identified three immunologic subtypes across samples and extensive within-patient diversity. Epithelial CD8+ TILs negatively associated with malignant diversity, reflecting immunological pruning of tumor clones inferred by neoantigen depletion, HLA I loss of heterozygosity, and spatial tracking between T cell and tumor clones. In addition, combinatorial prognostic effects of mutational processes and immune properties were observed, illuminating how specific genomic aberration types associate with immune response and impact survival. We conclude that within-patient spatial immune microenvironment variation shapes intraperitoneal malignant spread, provoking new evolutionary perspectives on HGSC clonal dispersion.

The telomerase reverse transcriptase (TERT) gene is highly expressed in stem cells and silenced upon differentiation. Cancer cells can attain immortality by activating TERT to maintain telomere length and telomerase activity, a crucial step of tumorigenesis. Two somatic mutations in the TERT promoter (C228T; C250T) have been identified in multiple cancers, such as melanoma and glioblastoma, as gain-of-function mutations that promote transcriptional activation of TERT. A recent study investigating TERT promoter mutations in ovarian carcinomas found mutations in 15% of clear cell carcinomas. However, it is unknown whether these mutations are prevalent in adult-type granulosa cell tumors (AGCTs) of the ovary, Sertoli-Leydig cell tumors (SLCTs), and other malignant sex cord-stromal tumors. We performed whole-genome sequencing on ten AGCT cases with matched normal and identified the TERT C228T promoter mutation in 50% of cases. We found that AGCT cases with mutated TERT promoter have increased expression of TERT mRNA compared to those with wild-type TERT promoter. All five TERT promoter mutated cases had high levels of TERT mRNA expression, whereas three of the five wild-type TERT cases had no measurable TERT mRNA expression. There was a tendency towards longer telomere lengths in AGCT cases with the TERT promoter mutation relative to those without, although it was not statistically significant. These results suggest that telomerase may be activated by a different method in the cases with no TERT promoter mutations but have TERT mRNA expression. The remaining cases with neither TERT promoter mutations nor TERT mRNA expression likely maintain their telomeres using a telomerase-independent method, such as the alternative lengthening of telomeres pathway. TERT C228T allelic discrimination analysis of 331 AGCTs, 5 SLCTs, and 18 other malignant sex cord-stromal tumors detected the mutation in 56/247 (23%) of primary AGCTs, 22/84 (26%) of recurrent AGCTs, 1/5 (20%) of SLCTs and (0/18) 0% of other malignant sex cord-stromal tumors. The single SLCT case with the TERT promoter mutation was poorly differentiated and harbored the pathognomonic FOXL2 mutation of AGCT, suggesting this SLCT case may actually be an AGCT. In 204 AGCT cases with available survival data, there was a trend towards worse disease-specific survival in patients with the TERT promoter mutation compared to those without; however, statistical significance was not reached (p = 0.128, log-ranked test). In 5 AGCT cases with primary and recurrent tissues, we found that the TERT promoter mutation was absent in the primary tumors but present in the recurrent tumors, suggesting that TERT C228T mutation may play an active role in progression of AGCTs. Overall, we found that TERT C228T promoter mutation was most common in AGCTs among the different malignant sex cord-stromal tumors. Our data confirm the activation of telomerase in AGCTs via TERT C228T promoter mutation, although alternative telomerase activation methods in AGCTs may exist. Our results suggest that TERT activation may play a role in AGCT recurrence. As such, telomere biology may be important for the progression of AGCTs. This abstract is also being presented as Poster B54. Citation Format: Jessica A. Pilsworth, Dawn R. Cochrane, Zhouchunyang Xia, Geraldine Aubert, Anniina E. M. Färkkilä, Hugo M. Horlings, Satoshi Yanagida, Winnie Yang, Jamie L. P. Lim, Yikan Wang, Ali Bashashati, Jacqueline Keul, Adele Wong, Esther Oliva, Sohrab P. Shah, Stefan Kommoss, Friedrich Kommoss, Peter M. Lansdorp, Duncan M. Baird, David G. Huntsman. TERT is frequently mutated in adult-type granulosa cell tumors of the ovary compared to other malignant sex cord-stromal tumors. [abstract]. In: Proceedings of the AACR Conference: Addressing Critical Questions in Ovarian Cancer Research and Treatment; Oct 1-4, 2017; Pittsburgh, PA. Philadelphia (PA): AACR; Clin Cancer Res 2018;24(15_Suppl):Abstract nr PR13.