Chemotherapy induces myeloid-driven spatial T-cell exhaustion in ovarian cancer.

To uncover the intricate, chemotherapy-induced spatiotemporal remodeling of the tumor microenvironment, we conducted integrative spatial and molecular characterization of 97 high-grade serous ovarian cancer (HGSC) samples collected before and after chemotherapy. Using single-cell and spatial analyses, we identify increasingly versatile immune cell states, which form spatiotemporally dynamic microcommunities at the tumor-stroma interface. We demonstrate that chemotherapy triggers spatial redistribution and exhaustion of CD8+ T cells due to prolonged antigen presentation by macrophages, both within interconnected myeloid networks termed "Myelonets" and at the tumor stroma interface. Single-cell and spatial transcriptomics identifies prominent TIGIT-NECTIN2 ligand-receptor interactions induced by chemotherapy. Using a functional patient-derived immuno-oncology platform, we show that CD8+T-cell activity can be boosted by combining immune checkpoint blockade with chemotherapy. Our discovery of chemotherapy-induced myeloid-driven spatial T-cell exhaustion paves the way for novel immunotherapeutic strategies to unleash CD8+ T-cell-mediated anti-tumor immunity in HGSC.

Genome-wide quantification of copy-number aberration impact on gene expression in ovarian high-grade serous carcinoma.

Copy-number alterations (CNAs) are a hallmark of cancer and can regulate cancer cell states via altered gene expression values. Herein, we have developed a copy-number impact (CNI) analysis method that quantifies the degree to which a gene expression value is impacted by CNAs and leveraged this analysis at the pathway level. Our results show that a high CNA is not necessarily reflected at the gene expression level, and our method is capable of detecting genes and pathways whose activity is strongly influenced by CNAs. Furthermore, the CNI analysis enables unbiased categorization of CNA categories, such as deletions and amplifications. We identified six CNI-driven pathways associated with poor treatment response in ovarian high-grade serous carcinoma (HGSC), which we found to be the most CNA-driven cancer across 14 cancer types. The key driver in most of these pathways was amplified wild-type KRAS, which we validated functionally using CRISPR modulation. Our results suggest that wild-type KRAS amplification is a driver of chemotherapy resistance in HGSC and may serve as a potential treatment target.

Tracing back primed resistance in cancer via sister cells.

Exploring non-genetic evolution of cell states during cancer treatments has become attainable by recent advances in lineage-tracing methods. However, transcriptional changes that drive cells into resistant fates may be subtle, necessitating high resolution analysis. Here, we present ReSisTrace that uses shared transcriptomic features of sister cells to predict the states priming treatment resistance. Applying ReSisTrace in ovarian cancer cells perturbed with olaparib, carboplatin or natural killer (NK) cells reveals pre-resistant phenotypes defined by proteostatic and mRNA surveillance features, reflecting traits enriched in the upcoming subclonal selection. Furthermore, we show that DNA repair deficiency renders cells susceptible to both DNA damaging agents and NK killing in a context-dependent manner. Finally, we leverage the obtained pre-resistance profiles to predict and validate small molecules driving cells to sensitive states prior to treatment. In summary, ReSisTrace resolves pre-existing transcriptional features of treatment vulnerability, facilitating both molecular patient stratification and discovery of synergistic pre-sensitizing therapies.

Locus-specific LINE-1 expression in clinical ovarian cancer specimens at the single-cell level.

Long interspersed nuclear elements (LINE-1s/L1s) are a group of retrotransposons that can copy themselves within a genome. In humans, it is the most successful transposon in nucleotide content. L1 expression is generally mild in normal human tissues, but the activity has been shown to increase significantly in many cancers. Few studies have examined L1 expression at single-cell resolution, thus it is undetermined whether L1 reactivation occurs solely in malignant cells within tumors. One of the cancer types with frequent L1 activity is high-grade serous ovarian carcinoma (HGSOC). Here, we identified locus-specific L1 expression with 3' single-cell RNA sequencing in pre- and post-chemotherapy HGSOC sample pairs from 11 patients, and in fallopian tube samples from five healthy women. Although L1 expression quantification with the chosen technique was challenging due to the repetitive nature of the element, we found evidence of L1 expression primarily in cancer cells, but also in other cell types, e.g. cancer-associated fibroblasts. The expression levels were similar in samples taken before and after neoadjuvant chemotherapy, indicating that L1 transcriptional activity was unaffected by clinical platinum-taxane treatment. Furthermore, L1 activity was negatively associated with the expression of MYC target genes, a finding that supports earlier literature of MYC being an L1 suppressor.

The Emerging Role of the Single-Cell and Spatial Tumor Microenvironment in High-Grade Serous Ovarian Cancer.

The development of single-cell and spatial technologies has enabled a more detailed understanding of the tumor microenvironment and its role in therapy response and clinical outcome of high-grade serous ovarian cancer (HGSC). Interestingly, emerging evidence suggests that HGSCs with different genetic drivers harbor distinct tumor-immune microenvironments. Further, spatial cell-cell interactions have been shown to shape the CD8+ T-cell phenotypes and responses to immune checkpoint blockade therapies. The heterogeneous stroma consisting of cancer-associated fibroblast (CAF) subtypes, endothelia, and site-specific stromal types such as mesothelium modulates treatment responses via increasing stiffness and by producing ligands that promote drug resistance, angiogenesis, or immune escape. Chemotherapy itself shifts CAFs toward an inflammatory phenotype that associates with poor survival and immune-suppressive signaling. New emerging immunotherapies include combinational approaches and agents targeting, for example, the tumor-intrinsic endoplasmic reticulum pathway. A more detailed understanding of the spatial interplay of tumor, immune, and stromal cells in the tumor microenvironment is needed to develop more efficient immunotherapeutic strategies for HGSC.

A platform for efficient establishment and drug-response profiling of high-grade serous ovarian cancer organoids.

The broad research use of organoids from high-grade serous ovarian cancer (HGSC) has been hampered by low culture success rates and limited availability of fresh tumor material. Here, we describe a method for generation and long-term expansion of HGSC organoids with efficacy markedly improved over previous reports (53% vs. 23%-38%). We established organoids from cryopreserved material, demonstrating the feasibility of using viably biobanked tissue for HGSC organoid derivation. Genomic, histologic, and single-cell transcriptomic analyses revealed that organoids recapitulated genetic and phenotypic features of original tumors. Organoid drug responses correlated with clinical treatment outcomes, although in a culture conditions-dependent manner and only in organoids maintained in human plasma-like medium (HPLM). Organoids from consenting patients are available to the research community through a public biobank and organoid genomic data are explorable through an interactive online tool. Taken together, this resource facilitates the application of HGSC organoids in basic and translational ovarian cancer research.

Transgelin 2 guards T cell lipid metabolic programming and anti-tumor function.

Mounting effective immunity against pathogens and tumors relies on the successful metabolic programming of T cells by extracellular fatty acids1-3. During this process, fatty-acid-binding protein 5 (FABP5) imports lipids that fuel mitochondrial respiration and sustain the bioenergetic requirements of protective CD8+ T cells4,5. Importantly, however, the mechanisms governing this crucial immunometabolic axis remain unexplored. Here we report that the cytoskeletal organizer Transgelin 2 (TAGLN2) is necessary for optimal CD8+ T cell fatty acid uptake, mitochondrial respiration, and anti-cancer function. We found that TAGLN2 interacts with FABP5, enabling the surface localization of this lipid importer on activated CD8+ T cells. Analysis of ovarian cancer specimens revealed that endoplasmic reticulum (ER) stress responses elicited by the tumor microenvironment repress TAGLN2 in infiltrating CD8+ T cells, enforcing their dysfunctional state. Restoring TAGLN2 expression in ER-stressed CD8+ T cells bolstered their lipid uptake, mitochondrial respiration, and cytotoxic capacity. Accordingly, chimeric antigen receptor T cells overexpressing TAGLN2 bypassed the detrimental effects of tumor-induced ER stress and demonstrated superior therapeutic efficacy in mice with metastatic ovarian cancer. Our study unveils the role of cytoskeletal TAGLN2 in T cell lipid metabolism and highlights the potential to enhance cellular immunotherapy in solid malignancies by preserving the TAGLN2-FABP5 axis.

Upregulation of ribosome biogenesis via canonical E-boxes is required for Myc-driven proliferation.

The transcription factor Myc drives cell growth across animal phyla and is activated in most forms of human cancer. However, it is unclear which Myc target genes need to be regulated to induce growth and whether multiple targets act additively or if induction of each target is individually necessary. Here, we identified Myc target genes whose regulation is conserved between humans and flies and deleted Myc-binding sites (E-boxes) in the promoters of fourteen of these genes in Drosophila. E-box mutants of essential genes were homozygous viable, indicating that the E-boxes are not required for basal expression. Eight E-box mutations led to Myc-like phenotypes; the strongest mutant, ppanEbox-/-, also made the flies resistant to Myc-induced cell growth without affecting Myc-induced apoptosis. The ppanEbox-/- flies are healthy and display only a minor developmental delay, suggesting that it may be possible to treat or prevent tumorigenesis by targeting individual downstream targets of Myc.

Longitudinal single-cell RNA-seq analysis reveals stress-promoted chemoresistance in metastatic ovarian cancer.

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.

Network-guided identification of cancer-selective combinatorial therapies in ovarian cancer.

Each patient's cancer consists of multiple cell subpopulations that are inherently heterogeneous and may develop differing phenotypes such as drug sensitivity or resistance. A personalized treatment regimen should therefore target multiple oncoproteins in the cancer cell populations that are driving the treatment resistance or disease progression in a given patient to provide maximal therapeutic effect, while avoiding severe co-inhibition of non-malignant cells that would lead to toxic side effects. To address the intra- and inter-tumoral heterogeneity when designing combinatorial treatment regimens for cancer patients, we have implemented a machine learning-based platform to guide identification of safe and effective combinatorial treatments that selectively inhibit cancer-related dysfunctions or resistance mechanisms in individual patients. In this case study, we show how the platform enables prediction of cancer-selective drug combinations for patients with high-grade serous ovarian cancer using single-cell imaging cytometry drug response assay, combined with genome-wide transcriptomic and genetic profiles. The platform makes use of drug-target interaction networks to prioritize those combinations that warrant further preclinical testing in scarce patient-derived primary cells. During the case study in ovarian cancer patients, we investigated (i) the relative performance of various ensemble learning algorithms for drug response prediction, (ii) the use of matched single-cell RNA-sequencing data to deconvolute cell population-specific transcriptome profiles from bulk RNA-seq data, (iii) and whether multi-patient or patient-specific predictive models lead to better predictive accuracy. The general platform and the comparison results are expected to become useful for future studies that use similar predictive approaches also in other cancer types.

miQC: An adaptive probabilistic framework for quality control of single-cell RNA-sequencing data.

Single-cell RNA-sequencing (scRNA-seq) has made it possible to profile gene expression in tissues at high resolution. An important preprocessing step prior to performing downstream analyses is to identify and remove cells with poor or degraded sample quality using quality control (QC) metrics. Two widely used QC metrics to identify a 'low-quality' cell are (i) if the cell includes a high proportion of reads that map to mitochondrial DNA (mtDNA) encoded genes and (ii) if a small number of genes are detected. Current best practices use these QC metrics independently with either arbitrary, uniform thresholds (e.g. 5%) or biological context-dependent (e.g. species) thresholds, and fail to jointly model these metrics in a data-driven manner. Current practices are often overly stringent and especially untenable on certain types of tissues, such as archived tumor tissues, or tissues associated with mitochondrial function, such as kidney tissue [1]. We propose a data-driven QC metric (miQC) that jointly models both the proportion of reads mapping to mtDNA genes and the number of detected genes with mixture models in a probabilistic framework to predict the low-quality cells in a given dataset. We demonstrate how our QC metric easily adapts to different types of single-cell datasets to remove low-quality cells while preserving high-quality cells that can be used for downstream analyses. Our software package is available at https://bioconductor.org/packages/miQC.

PRISM: recovering cell-type-specific expression profiles from individual composite RNA-seq samples.

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.

Contribution of allelic imbalance to colorectal cancer.

Point mutations in cancer have been extensively studied but chromosomal gains and losses have been more challenging to interpret due to their unspecific nature. Here we examine high-resolution allelic imbalance (AI) landscape in 1699 colorectal cancers, 256 of which have been whole-genome sequenced (WGSed). The imbalances pinpoint 38 genes as plausible AI targets based on previous knowledge. Unbiased CRISPR-Cas9 knockout and activation screens identified in total 79 genes within AI peaks regulating cell growth. Genetic and functional data implicate loss of TP53 as a sufficient driver of AI. The WGS highlights an influence of copy number aberrations on the rate of detected somatic point mutations. Importantly, the data reveal several associations between AI target genes, suggesting a role for a network of lineage-determining transcription factors in colorectal tumorigenesis. Overall, the results unravel the contribution of AI in colorectal cancer and provide a plausible explanation why so few genes are commonly affected by point mutations in cancers.

A Functional Homologous Recombination Assay Predicts Primary Chemotherapy Response and Long-Term Survival in Ovarian Cancer Patients.

Purpose: Homologous recombination deficiency (HRD) correlates with platinum sensitivity in patients with ovarian cancer, which clinically is the most useful predictor of sensitivity to PARPi. To date, there are no reliable diagnostic tools to anticipate response to platinum-based chemotherapy, thus we aimed to develop an ex vivo functional HRD detection test that could predict both platinum-sensitivity and patient eligibility to targeted drug treatments.Experimental Design: We obtained a functional HR score by quantifying homologous recombination (HR) repair after ionizing radiation-induced DNA damage in primary ovarian cancer samples (n = 32). Samples clustered in 3 categories: HR-deficient, HR-low, and HR-proficient. We analyzed the HR score association with platinum sensitivity and treatment response, platinum-free interval (PFI) and overall survival (OS), and compared it with other clinical parameters. In parallel, we performed DNA-sequencing of HR genes to assess if functional HRD can be predicted by currently offered genetic screening.Results: Low HR scores predicted primary platinum sensitivity with high statistical significance (P = 0.0103), associated with longer PFI (HR-deficient vs. HR-proficient: 531 vs. 53 days), and significantly correlated with improved OS (HR score <35 vs. ≥35, hazard ratio = 0.08, P = 0.0116). At the genomic level, we identified a few unclear mutations in HR genes and the mutational signature associated with HRD, but, overall, genetic screening failed to predict functional HRD.Conclusions: We developed an ex vivo assay that detects tumor functional HRD and an HR score able to predict platinum sensitivity, which holds the clinically relevant potential to become the routine companion diagnostic in the management of patients with ovarian cancer. Clin Cancer Res; 24(18); 4482-93. ©2018 AACR.

Genome-wide screen of cell-cycle regulators in normal and tumor cells identifies a differential response to nucleosome depletion.

To identify cell cycle regulators that enable cancer cells to replicate DNA and divide in an unrestricted manner, we performed a parallel genome-wide RNAi screen in normal and cancer cell lines. In addition to many shared regulators, we found that tumor and normal cells are differentially sensitive to loss of the histone genes transcriptional regulator CASP8AP2. In cancer cells, loss of CASP8AP2 leads to a failure to synthesize sufficient amount of histones in the S-phase of the cell cycle, resulting in slowing of individual replication forks. Despite this, DNA replication fails to arrest, and tumor cells progress in an elongated S-phase that lasts several days, finally resulting in death of most of the affected cells. In contrast, depletion of CASP8AP2 in normal cells triggers a response that arrests viable cells in S-phase. The arrest is dependent on p53, and preceded by accumulation of markers of DNA damage, indicating that nucleosome depletion is sensed in normal cells via a DNA-damage -like response that is defective in tumor cells.

Integrated data analysis reveals uterine leiomyoma subtypes with distinct driver pathways and biomarkers.

Uterine leiomyomas are common benign smooth muscle tumors that impose a major burden on women's health. Recent sequencing studies have revealed recurrent and mutually exclusive mutations in leiomyomas, suggesting the involvement of molecularly distinct pathways. In this study, we explored transcriptional differences among leiomyomas harboring different genetic drivers, including high mobility group AT-hook 2 (HMGA2) rearrangements, mediator complex subunit 12 (MED12) mutations, biallelic inactivation of fumarate hydratase (FH), and collagen, type IV, alpha 5 and collagen, type IV, alpha 6 (COL4A5-COL4A6) deletions. We also explored the transcriptional consequences of 7q22, 22q, and 1p deletions, aiming to identify possible target genes. We investigated 94 leiomyomas and 60 corresponding myometrial tissues using exon arrays, whole genome sequencing, and SNP arrays. This integrative approach revealed subtype-specific expression changes in key driver pathways, including Wnt/ß-catenin, Prolactin, and insulin-like growth factor (IGF)1 signaling. Leiomyomas with HMGA2 aberrations displayed highly significant up-regulation of the proto-oncogene pleomorphic adenoma gene 1 (PLAG1), suggesting that HMGA2 promotes tumorigenesis through PLAG1 activation. This was supported by the identification of genetic PLAG1 alterations resulting in expression signatures as seen in leiomyomas with HMGA2 aberrations. RAD51 paralog B (RAD51B), the preferential translocation partner of HMGA2, was up-regulated in MED12 mutant lesions, suggesting a role for this gene in the genesis of leiomyomas. FH-deficient leiomyomas were uniquely characterized by activation of nuclear factor erythroid 2-related factor 2 (NRF2) target genes, supporting the hypothesis that accumulation of fumarate leads to activation of the oncogenic transcription factor NRF2. This study emphasizes the need for molecular stratification in leiomyoma research and possibly in clinical practice as well. Further research is needed to determine whether the candidate biomarkers presented herein can provide guidance for managing the millions of patients affected by these lesions.

The Glanville fritillary genome retains an ancient karyotype and reveals selective chromosomal fusions in Lepidoptera.

Previous studies have reported that chromosome synteny in Lepidoptera has been well conserved, yet the number of haploid chromosomes varies widely from 5 to 223. Here we report the genome (393 Mb) of the Glanville fritillary butterfly (Melitaea cinxia; Nymphalidae), a widely recognized model species in metapopulation biology and eco-evolutionary research, which has the putative ancestral karyotype of n=31. Using a phylogenetic analyses of Nymphalidae and of other Lepidoptera, combined with orthologue-level comparisons of chromosomes, we conclude that the ancestral lepidopteran karyotype has been n=31 for at least 140 My. We show that fusion chromosomes have retained the ancestral chromosome segments and very few rearrangements have occurred across the fusion sites. The same, shortest ancestral chromosomes have independently participated in fusion events in species with smaller karyotypes. The short chromosomes have higher rearrangement rate than long ones. These characteristics highlight distinctive features of the evolutionary dynamics of butterflies and moths.

Uterine leiomyoma-linked MED12 mutations disrupt mediator-associated CDK activity.

Somatic mutations in exon 2 of the RNA polymerase II transcriptional Mediator subunit MED12 occur at very high frequency (∼70%) in uterine leiomyomas. However, the influence of these mutations on Mediator function and the molecular basis for their tumorigenic potential remain unknown. To clarify the impact of these mutations, we used affinity-purification mass spectrometry to establish the global protein-protein interaction profiles for both wild-type and mutant MED12. We found that uterine leiomyoma-linked mutations in MED12 led to a highly specific decrease in its association with Cyclin C-CDK8/CDK19 and loss of Mediator-associated CDK activity. Mechanistically, this occurs through disruption of a MED12-Cyclin C binding interface that we also show is required for MED12-mediated stimulation of Cyclin C-dependent CDK8 kinase activity. These findings indicate that uterine leiomyoma-linked mutations in MED12 uncouple Cyclin C-CDK8/19 from core Mediator and further identify the MED12/Cyclin C interface as a prospective therapeutic target in CDK8-driven cancers.

Transcriptional networks controlling the cell cycle.

In this work, we map the transcriptional targets of 107 previously identified Drosophila genes whose loss caused the strongest cell-cycle phenotypes in a genome-wide RNA interference screen and mine the resulting data computationally. Besides confirming existing knowledge, the analysis revealed several regulatory systems, among which were two highly-specific and interconnected feedback circuits, one between the ribosome and the proteasome that controls overall protein homeostasis, and the other between the ribosome and Myc/Max that regulates the protein synthesis capacity of cells. We also identified a set of genes that alter the timing of mitosis without affecting gene expression, indicating that the cyclic transcriptional program that produces the components required for cell division can be partially uncoupled from the cell division process itself. These genes all have a function in a pathway that regulates the phosphorylation state of Cdk1. We provide evidence showing that this pathway is involved in regulation of cell size, indicating that a Cdk1-regulated cell size checkpoint exists in metazoans.