Identification of EMT signaling cross-talk and gene regulatory networks by single-cell RNA sequencing.

The epithelial-to-mesenchymal transition (EMT) plays a critical role during normal development and in cancer progression. EMT is induced by various signaling pathways, including TGF-ß, BMP, Wnt-ß-catenin, NOTCH, Shh, and receptor tyrosine kinases. In this study, we performed single-cell RNA sequencing on MCF10A cells undergoing EMT by TGF-ß1 stimulation. Our comprehensive analysis revealed that cells progress through EMT at different paces. Using pseudotime clustering reconstruction of gene-expression profiles during EMT, we found sequential and parallel activation of EMT signaling pathways. We also observed various transitional cellular states during EMT. We identified regulatory signaling nodes that drive EMT with the expression of important microRNAs and transcription factors. Using a random circuit perturbation methodology, we demonstrate that the NOTCH signaling pathway acts as a key driver of TGF-ß-induced EMT. Furthermore, we demonstrate that the gene signatures of pseudotime clusters corresponding to the intermediate hybrid EMT state are associated with poor patient outcome. Overall, this study provides insight into context-specific drivers of cancer progression and highlights the complexities of the EMT process.

Chromatin state dynamics confers specific therapeutic strategies in enhancer subtypes of colorectal cancer.

OBJECTIVE: Enhancer aberrations are beginning to emerge as a key epigenetic feature of colorectal cancers (CRC), however, a comprehensive knowledge of chromatin state patterns in tumour progression, heterogeneity of these patterns and imparted therapeutic opportunities remain poorly described. DESIGN: We performed comprehensive epigenomic characterisation by mapping 222 chromatin profiles from 69 samples (33 colorectal adenocarcinomas, 4 adenomas, 21 matched normal tissues and 11 colon cancer cell lines) for six histone modification marks: H3K4me3 for Pol II-bound and CpG-rich promoters, H3K4me1 for poised enhancers, H3K27ac for enhancers and transcriptionally active promoters, H3K79me2 for transcribed regions, H3K27me3 for polycomb repressed regions and H3K9me3 for heterochromatin. RESULTS: We demonstrate that H3K27ac-marked active enhancer state could distinguish between different stages of CRC progression. By epigenomic editing, we present evidence that gains of tumour-specific enhancers for crucial oncogenes, such as ASCL2 and FZD10, was required for excessive proliferation. Consistently, combination of MEK plus bromodomain inhibition was found to have synergistic effects in CRC patient-derived xenograft models. Probing intertumour heterogeneity, we identified four distinct enhancer subtypes (EPIgenome-based Classification, EpiC), three of which correlate well with previously defined transcriptomic subtypes (consensus molecular subtypes, CMSs). Importantly, CMS2 can be divided into two EpiC subgroups with significant survival differences. Leveraging such correlation, we devised a combinatorial therapeutic strategy of enhancer-blocking bromodomain inhibitors with pathway-specific inhibitors (PARPi, EGFRi, TGFßi, mTORi and SRCi) for EpiC groups. CONCLUSION: Our data suggest that the dynamics of active enhancer underlies CRC progression and the patient-specific enhancer patterns can be leveraged for precision combination therapy.

The Cancer Genome Atlas (TCGA): an immeasurable source of knowledge.

The Cancer Genome Atlas (TCGA) is a public funded project that aims to catalogue and discover major cancer-causing genomic alterations to create a comprehensive "atlas" of cancer genomic profiles. So far, TCGA researchers have analysed large cohorts of over 30 human tumours through large-scale genome sequencing and integrated multi-dimensional analyses. Studies of individual cancer types, as well as comprehensive pan-cancer analyses have extended current knowledge of tumorigenesis. A major goal of the project was to provide publicly available datasets to help improve diagnostic methods, treatment standards, and finally to prevent cancer. This review discusses the current status of TCGA Research Network structure, purpose, and achievements.

Computational characterisation of cancer molecular profiles derived using next generation sequencing.

Our current understanding of cancer genetics is grounded on the principle that cancer arises from a clone that has accumulated the requisite somatically acquired genetic aberrations, leading to the malignant transformation. It also results in aberrent of gene and protein expression. Next generation sequencing (NGS) or deep sequencing platforms are being used to create large catalogues of changes in copy numbers, mutations, structural variations, gene fusions, gene expression, and other types of information for cancer patients. However, inferring different types of biological changes from raw reads generated using the sequencing experiments is algorithmically and computationally challenging. In this article, we outline common steps for the quality control and processing of NGS data. We highlight the importance of accurate and application-specific alignment of these reads and the methodological steps and challenges in obtaining different types of information. We comment on the importance of integrating these data and building infrastructure to analyse it. We also provide exhaustive lists of available software to obtain information and point the readers to articles comparing software for deeper insight in specialised areas. We hope that the article will guide readers in choosing the right tools for analysing oncogenomic datasets.

Dysregulation and Epigenetic Reprogramming of NRF2 Signaling Axis Promote Acquisition of Cisplatin Resistance and Metastasis in Head and Neck Squamous Cell Carcinoma.

PURPOSE: Cisplatin (CDDP)-based chemotherapy is a first-line treatment for patients with advanced head and neck squamous cell carcinomas (HNSCC), despite a high rate of treatment failures, acquired resistance, and subsequent aggressive behavior. The purpose of this study was to study the mechanism of CDDP resistance and metastasis in HNSCC. We investigated the role of NRF2 pathway activation as a driven event for tumor progression and metastasis of HNSCC. EXPERIMENTAL DESIGN: Human HNSCC cell lines that are highly resistant to CDDP were generated. Clonogenic survival assays and a mouse model of oral cancer were used to examine the impact of NRF2 activation in vitro and in vivo on CDDP sensitivity and development of metastasis. Western blotting, immunostaining, whole-exome sequencing, single-cell transcriptomic and epigenomic profiling platforms were performed to dissect clonal evolution and molecular mechanisms. RESULTS: Implantation of CDDP-resistant HNSCC cells into the tongues of nude mice resulted in a very high rate of distant metastases. The CDDP-resistant cells had significantly higher expression of NRF2 pathway genes in the presence of newly acquired KEAP1 mutations, or via epigenomic activation of target genes. Knockdown of NRF2 or restoration of the wild-type KEAP1 genes resensitized resistant cells to CDDP and decreased distant metastasis (DM). Finally, treatment with inhibitor of glutaminase-1, a NRF2 target gene, alleviated CDDP resistance. CONCLUSIONS: CDDP resistance and development of DM are associated with dysregulated and epigenetically reprogrammed KEAP1-NRF2 signaling pathway. A strategy targeting KEAP1/NRF2 pathway or glutamine metabolism deserves further clinical investigation in patients with CDDP-resistant head and neck tumors.

Single-cell chromatin accessibility profiling of acute myeloid leukemia reveals heterogeneous lineage composition upon therapy-resistance.

Acute myeloid leukemia (AML) is a heterogeneous disease characterized by high rate of therapy resistance. Since the cell of origin can impact response to therapy, it is crucial to understand the lineage composition of AML cells at time of therapy resistance. Here we leverage single-cell chromatin accessibility profiling of 22 AML bone marrow aspirates from eight patients at time of therapy resistance and following subsequent therapy to characterize their lineage landscape. Our findings reveal a complex lineage architecture of therapy-resistant AML cells that are primed for stem and progenitor lineages and spanning quiescent, activated and late stem cell/progenitor states. Remarkably, therapy-resistant AML cells are also composed of cells primed for differentiated myeloid, erythroid and even lymphoid lineages. The heterogeneous lineage composition persists following subsequent therapy, with early progenitor-driven features marking unfavorable prognosis in The Cancer Genome Atlas AML cohort. Pseudotime analysis further confirms the vast degree of heterogeneity driven by the dynamic changes in chromatin accessibility. Our findings suggest that therapy-resistant AML cells are characterized not only by stem and progenitor states, but also by a continuum of differentiated cellular lineages. The heterogeneity in lineages likely contributes to their therapy resistance by harboring different degrees of lineage-specific susceptibilities to therapy.

Radiographic response to neoadjuvant therapy in pleural mesothelioma should serve as a guide for patient selection for cytoreductive operations.

Background: Malignant pleural mesothelioma (MPM) is associated with poor prognosis despite advances in multimodal therapeutic strategies. While patients with resectable disease may benefit from added survival with oncologic resection, patient selection for mesothelioma operations often relies on both objective and subjective evaluation metrics. We sought to evaluate factors associated with improved overall survival (OS) in patients with mesothelioma who underwent macroscopic complete resection (MCR). Methods: Patients with MPM who received neoadjuvant therapy and underwent MCR were identified in a prospectively maintained departmental database. Clinicopathologic, blood-based, and radiographic variables were collected and included in a Cox regression analysis (CRA). Response to neoadjuvant therapy was characterized by a change in tumor thickness from pretherapy to preoperative scans using the modified RECIST criteria. Results: In this study, 99 patients met the inclusion criteria. The median age of the included patients was 64.7 years, who were predominantly men, had smoking and asbestos exposure, and who received neoadjuvant therapy. The median change in tumor thickness following neoadjuvant therapy was -16.5% (interquartile range of -49.7% to +14.2%). CRA demonstrated reduced OS associated with non-epithelioid histology [hazard ratio (HR): 3.06, 95% confidence interval (CI): 1.62-5.78, p < 0.001] and a response to neoadjuvant therapy inferior to the median (HR: 2.70, CI: 1.55-4.72, p < 0.001). Patients who responded poorly (below median) to neoadjuvant therapy had lower median survival (15.8 months compared to 38.2 months, p < 0.001). Conclusion: Poor response to neoadjuvant therapy in patients with MPM is associated with poor outcomes even following maximum surgical cytoreduction and should warrant a patient-centered discussion regarding goals of care and may therefore help guide further therapeutic decisions.

Immune landscape of a genetically engineered murine model of glioma compared with human glioma.

Novel therapeutic strategies targeting glioblastoma (GBM) often fail in the clinic, partly because preclinical models in which hypotheses are being tested do not recapitulate human disease. To address this challenge, we took advantage of our previously developed spontaneous Qk/Trp53/Pten (QPP) triple-knockout model of human GBM, comparing the immune microenvironment of QPP mice with that of patient-derived tumors to determine whether this model provides opportunity for gaining insights into tumor physiopathology and preclinical evaluation of therapeutic agents. Immune profiling analyses and single-cell sequencing of implanted and spontaneous tumors from QPP mice and from patients with glioma revealed intratumoral immune components that were predominantly myeloid cells (e.g., monocytes, macrophages, and microglia), with minor populations of T, B, and NK cells. When comparing spontaneous and implanted mouse samples, we found more neutrophils and T and NK cells in the implanted model. Neutrophils and T and NK cells were increased in abundance in samples derived from human high-grade glioma compared with those derived from low-grade glioma. Overall, our data demonstrate that our implanted and spontaneous QPP models recapitulate the immunosuppressive myeloid-dominant nature of the tumor microenvironment of human gliomas. Our model provides a suitable tool for investigating the complex immune compartment of gliomas.

Genomic Correlates of Outcome in Tumor-Infiltrating Lymphocyte Therapy for Metastatic Melanoma.

PURPOSE: Adoptive cell therapy (ACT) of tumor-infiltrating lymphocytes (TIL) historically yields a 40%-50% response rate in metastatic melanoma. However, the determinants of outcome are largely unknown. EXPERIMENTAL DESIGN: We investigated tumor-based genomic correlates of overall survival (OS), progression-free survival (PFS), and response to therapy by interrogating tumor samples initially collected to generate TIL infusion products. RESULTS: Whole-exome sequencing (WES) data from 64 samples indicated a positive correlation between neoantigen load and OS, but not PFS or response to therapy. RNA sequencing analysis of 34 samples showed that expression of PDE1C, RTKN2, and NGFR was enriched in responders who had improved PFS and OS. In contrast, the expression of ELFN1 was enriched in patients with unfavorable response, poor PFS and OS, whereas enhanced methylation of ELFN1 was observed in patients with favorable outcomes. Expression of ELFN1, NGFR, and PDE1C was mainly found in cancer-associated fibroblasts and endothelial cells in tumor tissues across different cancer types in publicly available single-cell RNA sequencing datasets, suggesting a role for elements of the tumor microenvironment in defining the outcome of TIL therapy. CONCLUSIONS: Our findings suggest that transcriptional features of melanomas correlate with outcomes after TIL therapy and may provide candidates to guide patient selection.

Inhibition of mitochondrial complex I reverses NOTCH1-driven metabolic reprogramming in T-cell acute lymphoblastic leukemia.

T-cell acute lymphoblastic leukemia (T-ALL) is commonly driven by activating mutations in NOTCH1 that facilitate glutamine oxidation. Here we identify oxidative phosphorylation (OxPhos) as a critical pathway for leukemia cell survival and demonstrate a direct relationship between NOTCH1, elevated OxPhos gene expression, and acquired chemoresistance in pre-leukemic and leukemic models. Disrupting OxPhos with IACS-010759, an inhibitor of mitochondrial complex I, causes potent growth inhibition through induction of metabolic shut-down and redox imbalance in NOTCH1-mutated and less so in NOTCH1-wt T-ALL cells. Mechanistically, inhibition of OxPhos induces a metabolic reprogramming into glutaminolysis. We show that pharmacological blockade of OxPhos combined with inducible knock-down of glutaminase, the key glutamine enzyme, confers synthetic lethality in mice harboring NOTCH1-mutated T-ALL. We leverage on this synthetic lethal interaction to demonstrate that IACS-010759 in combination with chemotherapy containing L-asparaginase, an enzyme that uncovers the glutamine dependency of leukemic cells, causes reduced glutaminolysis and profound tumor reduction in pre-clinical models of human T-ALL. In summary, this metabolic dependency of T-ALL on OxPhos provides a rational therapeutic target.

Reprogramming of bivalent chromatin states in NRAS mutant melanoma suggests PRC2 inhibition as a therapeutic strategy.

The dynamic evolution of chromatin state patterns during metastasis, their relationship with bona fide genetic drivers, and their therapeutic vulnerabilities are not completely understood. Combinatorial chromatin state profiling of 46 melanoma samples reveals an association of NRAS mutants with bivalent histone H3 lysine 27 trimethylation (H3K27me3) and Polycomb repressive complex 2. Reprogramming of bivalent domains during metastasis occurs on master transcription factors of a mesenchymal phenotype, including ZEB1, TWIST1, and CDH1. Resolution of bivalency using pharmacological inhibition of EZH2 decreases invasive capacity of melanoma cells and markedly reduces tumor burden in vivo, specifically in NRAS mutants. Coincident with bivalent reprogramming, the increased expression of pro-metastatic and melanocyte-specific cell-identity genes is associated with exceptionally wide H3K4me3 domains, suggesting a role for this epigenetic element. Overall, we demonstrate that reprogramming of bivalent and broad domains represents key epigenetic alterations in metastatic melanoma and that EZH2 plus MEK inhibition may provide a promising therapeutic strategy for NRAS mutant melanoma patients.

Single cell T cell landscape and T cell receptor repertoire profiling of AML in context of PD-1 blockade therapy.

In contrast to the curative effect of allogenic stem cell transplantation in acute myeloid leukemia via T cell activity, only modest responses are achieved with checkpoint-blockade therapy, which might be explained by T cell phenotypes and T cell receptor (TCR) repertoires. Here, we show by paired single-cell RNA analysis and TCR repertoire profiling of bone marrow cells in relapsed/refractory acute myeloid leukemia patients pre/post azacytidine+nivolumab treatment that the disease-related T cell subsets are highly heterogeneous, and their abundance changes following PD-1 blockade-based treatment. TCR repertoires expand and primarily emerge from CD8+ cells in patients responding to treatment or having a stable disease, while TCR repertoires contract in therapy-resistant patients. Trajectory analysis reveals a continuum of CD8+ T cell phenotypes, characterized by differential expression of granzyme B and a bone marrow-residing memory CD8+ T cell subset, in which a population with stem-like properties expressing granzyme K is enriched in responders. Chromosome 7/7q loss, on the other hand, is a cancer-intrinsic genomic marker of PD-1 blockade resistance in AML. In summary, our study reveals that adaptive T cell plasticity and genomic alterations determine responses to PD-1 blockade in acute myeloid leukemia.

Next-Generation Immunotherapies to Improve Anticancer Immunity.

Checkpoint inhibitors are widely used immunotherapies for advanced cancer. Nonetheless, checkpoint inhibitors have a relatively low response rate, work in a limited range of cancers, and have some unignorable side effects. Checkpoint inhibitors aim to reinvigorate exhausted or suppressed T cells in the tumor microenvironment (TME). However, the TME contains various other immune cell subsets that interact to determine the fate of cytotoxic T cells. Activation of cytotoxic T cells is initiated by antigen cross-presentation of dendritic cells. Dendritic cells could also release chemokines and cytokines to recruit and foster T cells. B cells, another type of antigen-presenting cell, also foster T cells and can produce tumor-specific antibodies. Neutrophils, a granulocyte cell subset in the TME, impede the proliferation and activation of T cells. The TME also consists of cytotoxic innate natural killer cells, which kill tumor cells efficiently. Natural killer cells can eradicate major histocompatibility complex I-negative tumor cells, which escape cytotoxic T cell-mediated destruction. A thorough understanding of the immune mechanism of the TME, as reviewed here, will lead to further development of more powerful therapeutic strategies. We have also reviewed the clinical outcomes of patients treated with drugs targeting these immune cells to identify strategies for improvement and possible immunotherapy combinations.

TRIM28 multi-domain protein regulates cancer stem cell population in breast tumor development.

The expression of Tripartite motif-containing protein 28 (TRIM28)/Krüppel-associated box (KRAB)-associated protein 1 (KAP1), is elevated in at least 14 tumor types, including solid and hematopoietic tumors. High level of TRIM28 is associated with triple-negative subtype of breast cancer (TNBC), which shows higher aggressiveness and lower survival rates. Interestingly, TRIM28 is essential for maintaining the pluripotent phenotype in embryonic stem cells. Following on that finding, we evaluated the role of TRIM28 protein in the regulation of breast cancer stem cells (CSC) populations and tumorigenesis in vitro and in vivo. Downregulation of TRIM28 expression in xenografts led to deceased expression of pluripotency and mesenchymal markers, as well as inhibition of signaling pathways involved in the complex mechanism of CSC maintenance. Moreover, TRIM28 depletion reduced the ability of cancer cells to induce tumor growth when subcutaneously injected in limiting dilutions. Our data demonstrate that the downregulation of TRIM28 gene expression reduced the ability of CSCs to self-renew that resulted in significant reduction of tumor growth. Loss of function of TRIM28 leads to dysregulation of cell cycle, cellular response to stress, cancer cell metabolism, and inhibition of oxidative phosphorylation. All these mechanisms directly regulate maintenance of CSC population. Our original results revealed the role of the TRIM28 in regulating the CSC population in breast cancer. These findings may pave the way to novel and more effective therapies targeting cancer stem cells in breast tumors.

Measles virus infection of human keratinocytes: Possible link between measles and atopic dermatitis.

BACKGROUND: Measles virus (MV) infection is marked with a skin rash in the acute phase of the disease, which pathogenesis remains poorly understood. Moreover, the association between measles and progression of skin diseases, such as atopic dermatitis (AD), is still elusive. OBJECTIVE: We have thus analysed the susceptibility of human keratinocytes to MV infection and explore the potential relationship between MV vaccination and the pathogenesis the AD. METHODS: We performed immunovirological characterisation of MV infection in human keratinocytes and then tested the effect of live attenuated measles vaccine on the progression of AD in adult patients, in a prospective, double-blind study. RESULTS: We showed that both human primary keratinocytes and the keratinocyte cell line HaCaT express MV receptors and could be infected by MV. The infection significantly modulated the expression of several keratinocyte-produced cytokines, known to be implicated in the pathogenesis of inflammatory allergic diseases, including AD. We then analysed the relationship between exposure to MV by vaccination and the progression of AD in 20 adults during six weeks. We found a significant decrease in CCL26 and thymic stromal lymphopoietin (TSLP) mRNA in biopsies from acute lesions of vaccinated patients, suggesting MV-induced modulation of skin cytokine expression. Clinical analysis revealed a transient improvement of SCORAD index in vaccinated compared to placebo-treated patients, two weeks after vaccination. CONCLUSIONS: Altogether, these results clearly demonstrate that keratinocytes are susceptible to MV infection, which could consequently modulate their cytokine production, resulting with a beneficial effect in the progression of AD. This study provides thus a proof of concept for the vaccination therapy in AD and may open new avenues for the development of novel strategies in the treatment of this allergic disease.

TRIM28 and Interacting KRAB-ZNFs Control Self-Renewal of Human Pluripotent Stem Cells through Epigenetic Repression of Pro-differentiation Genes.

Reprogramming to induced pluripotent stem cells (iPSCs) and differentiation of pluripotent stem cells (PSCs) are regulated by epigenetic machinery. Tripartite motif protein 28 (TRIM28), a universal mediator of Krüppel-associated box domain zinc fingers (KRAB-ZNFs), is known to regulate both processes; however, the exact mechanism and identity of participating KRAB-ZNF genes remain unknown. Here, using a reporter system, we show that TRIM28/KRAB-ZNFs alter DNA methylation patterns in addition to H3K9me3 to cause stable gene repression during reprogramming. Using several expression datasets, we identified KRAB-ZNFs (ZNF114, ZNF483, ZNF589) in the human genome that maintain pluripotency. Moreover, we identified target genes repressed by these KRAB-ZNFs. Mechanistically, we demonstrated that these KRAB-ZNFs directly alter gene expression of important developmental genes by modulating H3K9me3 and DNA methylation of their promoters. In summary, TRIM28 employs KRAB-ZNFs to evoke epigenetic silencing of its target differentiation genes via H3K9me3 and DNA methylation.