Mesenchymal-epithelial transition in lymph node metastases of oral squamous cell carcinoma is accompanied by ZEB1 expression.

BACKGROUND: Oral squamous cell carcinoma (OSCC), an HPV-negative head and neck cancer, frequently metastasizes to the regional lymph nodes but only occasionally beyond. Initial phases of metastasis are associated with an epithelial-mesenchymal transition (EMT), while the consolidation phase is associated with mesenchymal-epithelial transition (MET). This dynamic is referred to as epithelial-mesenchymal plasticity (EMP). While it is known that EMP is essential for cancer cell invasion and metastatic spread, less is known about the heterogeneity of EMP states and even less about the heterogeneity between primary and metastatic lesions. METHODS: To assess both the heterogeneity of EMP states in OSCC cells and their effects on stromal cells, we performed single-cell RNA sequencing (scRNAseq) of 5 primary tumors, 9 matching metastatic and 5 tumor-free lymph nodes and re-analyzed publicly available scRNAseq data of 9 additional primary tumors. For examining the cell type composition, we performed bulk transcriptome sequencing. Protein expression of selected genes were confirmed by immunohistochemistry. RESULTS: From the 23 OSCC lesions, the single cell transcriptomes of a total of 7263 carcinoma cells were available for in-depth analyses. We initially focused on one lesion to avoid confounding inter-patient heterogeneity and identified OSCC cells expressing genes characteristic of different epithelial and partial EMT stages. RNA velocity and the increase in inferred copy number variations indicated a progressive trajectory towards epithelial differentiation in this metastatic lesion, i.e., cells likely underwent MET. Extension to all samples revealed a less stringent but essentially similar pattern. Interestingly, MET cells show increased activity of the EMT-activator ZEB1. Immunohistochemistry confirmed that ZEB1 was co-expressed with the epithelial marker cornifin B in individual tumor cells. The lack of E-cadherin mRNA expression suggests this is a partial MET. Within the tumor microenvironment we found immunomodulating fibroblasts that were maintained in primary and metastatic OSCC. CONCLUSIONS: This study reveals that EMP enables different partial EMT and epithelial phenotypes of OSCC cells, which are endowed with capabilities essential for the different stages of the metastatic process, including maintenance of cellular integrity. During MET, ZEB1 appears to be functionally active, indicating a more complex role of ZEB1 than mere induction of EMT.

The HDAC inhibitor domatinostat induces type I interferon α in Merkel cell carcinoma by HES1 repression.

BACKGROUND: Class I selective histone deacetylase inhibitors (HDACi) have been previously demonstrated to not only increase major histocompatibility complex class I surface expression in Merkel cell carcinoma (MCC) cells by restoring the antigen processing and presentation machinery, but also exert anti-tumoral effect by inducing apoptosis. Both phenomena could be due to induction of type I interferons (IFN), as has been described for HDACi. However, the mechanism of IFN induction under HDACi is not fully understood because the expression of IFNs is regulated by both activating and inhibitory signaling pathways. Our own preliminary observations suggest that this may be caused by suppression of HES1. METHODS: The effect of the class I selective HDACi domatinostat and IFNα on cell viability and the apoptosis of MCPyV-positive (WaGa, MKL-1) and -negative (UM-MCC 34) MCC cell lines, as well as, primary fibroblasts were assessed by colorimetric methods or measuring mitochondrial membrane potential and intracellular caspase-3/7, respectively. Next, the impact of domatinostat on IFNA and HES1 mRNA expression was measured by RT-qPCR; intracellular IFNα production was detected by flow cytometry. To confirm that the expression of IFNα induced by HDACi was due to the suppression of HES1, it was silenced by RNA interference and then mRNA expression of IFNA and IFN-stimulated genes was assessed. RESULTS: Our studies show that the previously reported reduction in viability of MCC cell lines after inhibition of HDAC by domatinostat is accompanied by an increase in IFNα expression, both of mRNA and at the protein level. We confirmed that treatment of MCC cells with external IFNα inhibited their proliferation and induced apoptosis. Re-analysis of existing single-cell RNA sequencing data indicated that induction of IFNα by domatinostat occurs through repression of HES1, a transcriptional inhibitor of IFNA; this was confirmed by RT-qPCR. Finally, siRNA-mediated silencing of HES1 in the MCC cell line WaGa not only increased mRNA expression of IFNA and IFN-stimulated genes but also decreased cell viability. CONCLUSION: Our results demonstrate that the direct anti-tumor effect of HDACi domatinostat on MCC cells is at least in part mediated via decreased HES1 expression allowing the induction of IFNα, which in turn causes apoptosis.

B-cell lymphoma-extra large is a promising drug target in Merkel cell carcinoma.

BACKGROUND: Merkel cell carcinoma (MCC) is an aggressive skin tumour with neuroendocrine differentiation. Immunotherapies are effective in the treatment of patients with advanced-stage MCC, but for patients whose tumours cannot be controlled by the immune system, alternative approaches are urgently needed. OBJECTIVES: To identify overexpressed oncogenes as potential drug targets for MCC. METHODS: NanoString platform, digital droplet polymerase chain reaction (ddPCR) and fluorescence in situ hybridization (FISH) assays were used to determine copy number variations (CNVs); BCL2L1 and PARP1 mRNA expression levels were determined by quantitative real-time polymerase chain reaction (qRT-PCR), B-cell lymphoma extra-large (Bcl-xL) and poly (ADP-ribose) polymerase 1 (PARP1) protein by immuno-blot. Specific Bcl-xL inhibitors and a PARP1 inhibitor were used alone or in combination to test their antitumour effect. RESULTS: Screening for CNVs in 13 classic Merkel cell polyomavirus (MCPyV)-positive and MCPyV-negative MCC cell lines revealed BCL2L1 gains and amplifications, confirmed by ddPCR in 10 cell lines. By ddPCR and FISH, we demonstrated that BCL2L1 gains are present in tumour tissue. BCL2L1 copy number gains were associated with increased Bcl-xL mRNA and protein expression. However, high Bcl-xL expression was not restricted to MCC cells harbouring a BCL2L1 gain/amplification, suggesting additional epigenetic means of regulation. The functional relevance of Bcl-xL in MCC cells was demonstrated by the fact that specific Bcl-xL inhibitors (A1331852 and WEHI-539) led to the induction of apoptosis. Owing to the strong expression and activation of PARP1 in MCC cell lines, we next tested the combination of Bcl-xL inhibitors with the PARP1 inhibitor olaparib, which showed synergistic antitumour effects. CONCLUSIONS: Bcl-xL, which is highly expressed in MCC, appears to be an attractive therapeutic target for the treatment of this tumour, especially as the effect of specific Bcl-xL inhibitors is synergistically enhanced by simultaneous PARP inhibition.

Phenotypic plasticity of malignant T cells in blood and skin of a Sézary syndrome patient revealed by single cell transcriptomics.

Background: Sézary Syndrome (SS) is an aggressive leukemic variant of cutaneous T-cell lymphomas (CTCL). In SS patients, malignant T cells are circulating through the blood and cause erythroderma. Objective: To compare the transcriptome of single cells in blood and skin samples from a patient with advanced SS. Methods: We utilized combined single cell RNA and T-cell receptor (TCR) sequencing (scRNA-seq). Results: We scrutinized the malignant T cells in blood and skin in an unbiased manner without pre-sorting of cells. We observed different phenotypes of the same monoclonal malignant T-cell population, confirmed by TCR sequencing and inferred copy number variation analysis. Malignant T cells present in the circulating blood expressed genes resembling central memory T cells such as CCR7, IL7R and CD27. In the skin, we detected two major malignant T-cell populations: One subpopulation was closely related to the malignant T cells from the blood, while the other subpopulation expressed genes reminiscent of skin resident effector memory T cells including GZMB and NKG7. Pseudotime analysis indicated crucial transcriptomic changes in the transition of malignant T cells between blood and skin. These changes included the differential regulation of TXNIP, a putative tumor suppressor in CTCL, and the adaptation to the hypoxic conditions in the skin. Tumor cell proliferation in the skin was supported by stimulating interactions between myeloid cells and malignant T cells. Conclusions: Using scRNA-seq we detected a high degree of functional heterogeneity within the malignant T-cell population in SS and highlighted crucial differences between SS cells in blood and skin.

Identification of Epigenetically Regulated Genes Distinguishing Intracranial from Extracranial Melanoma Metastases.

Despite remarkable advances in treating patients with metastatic melanoma, the management of melanoma brain metastases remains challenging. Recent evidence suggests that epigenetic reprogramming is an important mechanism for the adaptation of melanoma cells to the brain environment. In this study, the methylomes and transcriptomes of a cohort of matched melanoma metastases were evaluated by integrated omics data analysis. The identified 38 candidate genes displayed distinct promoter methylation and corresponding gene expression changes in intracranial compared with extracranial metastases. The 11 most promising genes were validated on protein level in both tumor and surrounding normal tissue using immunohistochemistry. In accordance with the underlying promoter methylation and gene expression changes, a significantly different protein expression was confirmed for STK10, PDXK, WDR24, CSSP1, NMB, RASL11B, phosphorylated PRKCZ, PRKCZ, and phosphorylated GRB10 in the intracranial metastases. The observed changes imply a distinct intracranial phenotype with increased protein kinase B phosphorylation and a higher frequency of proliferating cells. Knockdown of PRKCZ or GRB10 altered the expression of phosphorylated protein kinase B and decreased the viability of a brain-specific melanoma cell line. In summary, epigenetically regulated cancer-relevant alterations were identified that provide insights into the molecular mechanisms that discriminate brain metastases from other organ metastases, which could be exploited by targeting the affected signaling pathways.

Intratumoral heterogeneity of MYC drives medulloblastoma metastasis and angiogenesis.

BACKGROUND: Intratumoral heterogeneity is crucially involved in metastasis, resistance to therapy, and cancer relapse. Amplifications of the proto-oncogene MYC display notable heterogeneity at the single-cell level and are associated with a particularly dismal prognosis in high-risk medulloblastomas (MBs). The aim of this study was to establish the relevance of interclonal cross-talk between MYC-driven and non-MYC-driven MB cells. METHODS: We used fluorescence in situ hybridization, single-cell transcriptomics, and immunohistochemistry, in vitro isogenic cell models, non-targeted proteomics, mass spectrometry-based metabolite quantification, HUVECs tube formation assay, and orthotopic in vivo experiments to investigate interclonal cross-talk in MB. RESULTS: We found that the release of lactate dehydrogenase A (LDHA) from MYC-driven cells facilitates metastatic seeding and outgrowth, while secretion of dickkopf WNT signaling pathway inhibitor 3 from non-MYC-driven cells promotes tumor angiogenesis. This tumor-supporting interaction between both subclones was abrogated by targeting the secretome through pharmacological and genetic inhibition of LDHA, which significantly suppressed tumor cell migration. CONCLUSION: Our study reveals the functional relevance of clonal diversity and highlights the therapeutic potential of targeting the secretome to interrupt interclonal communication and progression in high-risk MB.

Knockout of Factor-Inhibiting HIF (Hif1an) in Colon Epithelium Attenuates Chronic Colitis but Does Not Reduce Colorectal Cancer in Mice.

Inflammatory bowel disease such as chronic colitis promotes colorectal cancer, which is a common cause of cancer mortality worldwide. Hypoxia is a characteristic of inflammation as well as of solid tumors and enforces a gene expression response controlled by hypoxia-inducible factors (HIFs). Once established, solid tumors are immunosuppressive to escape their abatement through immune cells. Although HIF activity is known to 1) promote cancer development and 2) drive tumor immune suppression through the secretion of adenosine, both prolyl hydroxylases and an asparaginyl hydroxylase termed factor-inhibiting HIF (FIH) negatively regulate HIF. Thus, FIH may act as a tumor suppressor in colorectal cancer development. In this study, we examined the role of colon epithelial FIH in a mouse model of colitis-induced colorectal cancer. We recapitulated colitis-associated colorectal cancer development in mice using the azoxymethane/dextran sodium sulfate model in Vil1-Cre/FIH+f/+f and wild-type siblings. Colon samples were analyzed regarding RNA and protein expression and histology. Vil1-Cre/FIH+f/+f mice showed a less severe colitis progress compared with FIH+f/+f animals and a lower number of infiltrating macrophages in the inflamed tissue. RNA sequencing analyses of colon tissue revealed a lower expression of genes associated with the immune response in Vil1-Cre/FIH+f/+f mice. However, tumor occurrence did not significantly differ between Vil1-Cre/FIH+f/+f and wild-type mice. Thus, FIH knockout in colon epithelial cells did not modulate colorectal cancer development but reduced the inflammatory response in chronic colitis.

DNA-methylation patterns imply a common cellular origin of virus- and UV-associated Merkel cell carcinoma.

Merkel cell carcinoma (MCC) is a neuroendocrine tumor either induced by integration of the Merkel cell polyomavirus into the cell genome or by accumulation of UV-light-associated mutations (VP-MCC and UV-MCC). Whether VP- and UV-MCC have the same or different cellular origins is unclear; with mesenchymal or epidermal origins discussed. DNA-methylation patterns have a proven utility in determining cellular origins of cancers. Therefore, we used this approach to uncover evidence regarding the cell of origin of classical VP- and UV-MCC cell lines, i.e., cell lines with a neuroendocrine growth pattern (n = 9 and n = 4, respectively). Surprisingly, we observed high global similarities in the DNA-methylation of UV- and VP-MCC cell lines. CpGs of lower methylation in VP-MCC cell lines were associated with neuroendocrine marker genes such as SOX2 and INSM1, or linked to binding sites of EZH2 and SUZ12 of the polycomb repressive complex 2, i.e., genes with an impact on carcinogenesis and differentiation of neuroendocrine cancers. Thus, the observed differences appear to be rooted in viral compared to mutation-driven carcinogenesis rather than distinct cells of origin. To test this hypothesis, we used principal component analysis, to compare DNA-methylation data from different epithelial and non-epithelial neuroendocrine cancers and established a scoring model for epithelial and neuroendocrine characteristics. Subsequently, we applied this scoring model to the DNA-methylation data of the VP- and UV-MCC cell lines, revealing that both clearly scored as epithelial cancers. In summary, our comprehensive analysis of DNA-methylation suggests a common epithelial origin of UV- and VP-MCC cell lines.

Classical and Variant Merkel Cell Carcinoma Cell Lines Display Different Degrees of Neuroendocrine Differentiation and Epithelial-Mesenchymal Transition.

Merkel cell carcinoma (MCC) is an aggressive neuroendocrine skin cancer characterized by high invasiveness, early metastases, and high mortality. Because of the lack of suitable animal models, most functional studies are performed using cell lines, some of which lack classical neuroendocrine growth characteristics. Here, we scrutinized the molecular characteristics of classical MCC and variant MCC cell lines by differential gene expression and the respective epigenetic regulation by microRNAs and DNA methylation. Cutaneous squamous cell carcinoma cell lines were used for comparison. The most striking observation was a lower expression of epithelial-mesenchymal transition-related genes in classical MCCs, which was accompanied by higher expression of the epithelial-mesenchymal transition-regulating microRNA clusters miR-200c-141 and miR-183-96-182 and hypomethylation of the respective microRNA loci. Experimental expression of the MCC lineage factor ATOH1 in variant MCCs resulted in an increased expression of miR-200c-141 paralleled by a reduction of genes associated with epithelial-mesenchymal transition, thus demonstrating a connection between neuroendocrine characteristics and the lack of epithelial-mesenchymal transition. Together, our observations not only reinforce concerns about the use of variant MCCs as proper MCC representatives, but also suggest variant MCCs as cells locked in an intermediate state between neuroendocrine and epithelial differentiation.

Merkel cell carcinoma-derived exosome-shuttle miR-375 induces fibroblast polarization by inhibition of RBPJ and p53.

Merkel cell carcinoma (MCC) is a highly invasive and metastatic skin cancer. While high expression of miR-375 is a characteristic of MCC, it seems not to contribute to the malignant phenotype of MCC cells. miR-375 enrichment in MCC-derived extracellular vesicles suggests its intercellular signaling function. Here, we demonstrate that horizontally transferred miR-375 causes fibroblast polarization toward cancer-associated fibroblasts (CAFs). The polarization is evidenced by phenotypic changes and induction of α-SMA, CXCL2, and IL-1ß. Fibroblast polarization is inhibited by specific antagomirs and mimicked by experimental miR-375 expression. Mechanistically, miR-375 downregulates RBPJ and p53, two key players regulating fibroblast polarization. In clinical MCC samples, in situ hybridization located miR-375 in CAFs, which correlated with high α-SMA protein and low RBPJ and TP53 expression; single-cell RNAseq revealed a disparate fibroblast polarization negatively correlating with p53 pathway-related gene expression. Thus, the functional role of miR-375 in MCC is to generate a pro-tumorigenic microenvironment by inducing fibroblast polarization.

The HDAC Inhibitor Domatinostat Promotes Cell-Cycle Arrest, Induces Apoptosis, and Increases Immunogenicity of Merkel Cell Carcinoma Cells.

Merkel cell carcinoma (MCC) is a rare, highly aggressive skin cancer for which immune modulation by immune checkpoint inhibitors shows remarkable response rates. However, primary or secondary resistance to immunotherapy prevents benefits in a significant proportion of patients. For MCC, one immune escape mechanism is insufficient for recognition by T cells owing to the downregulation of major histocompatibility complex I surface expression. Histone deacetylase inhibitors have been demonstrated to epigenetically reverse the low major histocompatibility complex I expression caused by the downregulation of the antigen-processing machinery. Domatinostat, an orally available small-molecule inhibitor targeting histone deacetylase class I, is currently in clinical evaluation to overcome resistance to immunotherapy. In this study, we present preclinical data on domatinostat's efficacy and mode of action in MCC. Single-cell RNA sequencing revealed a distinct gene expression signature of antigen processing and presentation, cell-cycle arrest, and execution phase of apoptosis on treatment. Accordingly, functional assays showed that domatinostat induced G2M arrest and apoptosis. In the surviving cells, antigen-processing machinery component gene transcription and translation were upregulated, consequently resulting in increased major histocompatibility complex I surface expression. Altogether, domatinostat not only exerts direct antitumoral effects but also restores HLA class I surface expression on MCC cells, therefore, restoring surviving MCC cells' susceptibility to recognition and elimination by cognate cytotoxic T cells.

Patched 1 expression in Merkel cell carcinoma.

The relevance of Hedgehog signaling in Merkel cell carcinoma has only been addressed by a few studies with conflicting results. Thus, we aimed to establish the expression of Hedgehog signaling molecules in Merkel cell carcinoma to characterize causes of aberrant expression and to correlate these findings with the clinical course of the patients. Immunohistochemistry was performed for Sonic, Indian, Patched 1 (PTCH1) and Smoothened on patients' tumor tissue. Respective mRNA expression was analyzed in 10 Merkel cell carcinoma cell lines using quantitative real-time polymerase chain reaction. PTCH1 sequencing and DNA methylation microarray analyses were carried out on tumor tissues as well as cell lines. PTCH1 immunoreactivity in Merkel cell carcinoma was similar to that of basal cell carcinomas, which both significantly differed from PTCH1 immunoreactivity in healthy skin. Most PTCH1 mutations found were synonymous or without known functional impact. However, on average, the promoter regions of both PTCH1 were hypomethylated independently from PTCH1 gene expression or Merkel cell polyomavirus status. PTCH1 and GLI1/2/3 genes were differently expressed in different cell lines; notably, there was a significant correlation between GLI2 and PTCH1 mRNA expression. Similar to PTCH1 protein expression in patient tissues, PTCH1 gene expression in Merkel cell carcinoma cell lines is highly variable, but due to the similar methylation pattern across Merkel cell carcinoma cell lines, effects other than methylation seem to be the reason for the differential expression and PTCH1 appears to be upregulated by GLI as a classical Hedgehog target gene.

MCPyV Large T Antigen-Induced Atonal Homolog 1 Is a Lineage-Dependency Oncogene in Merkel Cell Carcinoma.

Despite the fact that the transcription factor ATOH1 is a master regulator of Merkel cell development, its role in Merkel cell carcinoma (MCC) carcinogenesis remains controversial. Here, we provide several lines of evidence that ATOH1 is a lineage-dependent oncogene in MCC. Luciferase assays revealed binding of ATOH1 and subsequent activation to the promoter of miR-375, which is one of the most abundant microRNAs in MCCs. Overexpression of ATOH1 in variant MCC cell lines and fibroblasts induced miR-375 expression, whereas ATOH1 knockdown in classical MCC cell lines reduced miR-375 expression. Moreover, ATOH1 overexpression in these cells changed their growth characteristics from adherent to suspension and/orspheroidal growth, that is, resembling the neuroendocrine growth pattern of classical MCC cell lines. Notably, ectopic expression of different Merkel cell polyomavirus (MCPyV)-derived truncated large T antigens induced ATOH1 expression in fibroblasts, which was paralleled by miR-375 expression and similar morphologic changes. In summary, MCPyV-associated carcinogenesis is likely to induce the characteristic neuroendocrine features of MCC via induction of ATOH1; thus, ATOH1 can be regarded as a lineage-dependent oncogene in MCC.

Conversion of Sox2-dependent Merkel cell carcinoma to a differentiated neuron-like phenotype by T antigen inhibition.

Viral cancers show oncogene addiction to viral oncoproteins, which are required for survival and proliferation of the dedifferentiated cancer cell. Human Merkel cell carcinomas (MCCs) that harbor a clonally integrated Merkel cell polyomavirus (MCV) genome have low mutation burden and require viral T antigen expression for tumor growth. Here, we showed that MCV+ MCC cells cocultured with keratinocytes undergo neuron-like differentiation with neurite outgrowth, secretory vesicle accumulation, and the generation of sodium-dependent action potentials, hallmarks of a neuronal cell lineage. Cocultured keratinocytes are essential for induction of the neuronal phenotype. Keratinocyte-conditioned medium was insufficient to induce this phenotype. Single-cell RNA sequencing revealed that T antigen knockdown inhibited cell cycle gene expression and reduced expression of key Merkel cell lineage/MCC marker genes, including HES6, SOX2, ATOH1, and KRT20 Of these, T antigen knockdown directly inhibited Sox2 and Atoh1 expression. MCV large T up-regulated Sox2 through its retinoblastoma protein-inhibition domain, which in turn activated Atoh1 expression. The knockdown of Sox2 in MCV+ MCCs mimicked T antigen knockdown by inducing MCC cell growth arrest and neuron-like differentiation. These results show Sox2-dependent conversion of an undifferentiated, aggressive cancer cell to a differentiated neuron-like phenotype and suggest that the ontology of MCC arises from a neuronal cell precursor.

Microsyntenic Clusters Reveal Conservation of lncRNAs in Chordates Despite Absence of Sequence Conservation.

Homologous long non-coding RNAs (lncRNAs) are elusive to identify by sequence similarity due to their fast-evolutionary rate. Here we develop LincOFinder, a pipeline that finds conserved intergenic lncRNAs (lincRNAs) between distant related species by means of microsynteny analyses. Using this tool, we have identified 16 bona fide homologous lincRNAs between the amphioxus and human genomes. We characterized and compared in amphioxus and Xenopus the expression domain of one of them, Hotairm1, located in the anterior part of the Hox cluster. In addition, we analyzed the function of this lincRNA in Xenopus, showing that its disruption produces a severe headless phenotype, most probably by interfering with the regulation of the Hox cluster. Our results strongly suggest that this lincRNA has probably been regulating the Hox cluster since the early origin of chordates. Our work pioneers the use of syntenic searches to identify non-coding genes over long evolutionary distances and helps to further understand lncRNA evolution.

Ultra Large Gene Families: A Matter of Adaptation or Genomic Parasites?

Gene duplication is an important mechanism of molecular evolution. It offers a fast track to modification, diversification, redundancy or rescue of gene function. However, duplication may also be neutral or (slightly) deleterious, and often ends in pseudo-geneisation. Here, we investigate the phylogenetic distribution of ultra large gene families on long and short evolutionary time scales. In particular, we focus on a family of NACHT-domain and leucine-rich-repeat-containing (NLR)-genes, which we previously found in large numbers to occupy one chromosome arm of the zebrafish genome. We were interested to see whether such a tight clustering is characteristic for ultra large gene families. Our data reconfirm that most gene family inflations are lineage-specific, but we can only identify very few gene clusters. Based on our observations we hypothesise that, beyond a certain size threshold, ultra large gene families continue to proliferate in a mechanism we term "run-away evolution". This process might ultimately lead to the failure of genomic integrity and drive species to extinction.