Characterizing DNA methylation signatures and their potential functional roles in Merkel cell carcinoma

Comprehensive DNA Methylation Profiling of Medullary Thyroid Carcinoma: Molecular Classification, Potential Therapeutic Target, and Classifier System

PURPOSE

Medullary thyroid carcinoma (MTC) presents a distinct biological context from other thyroid cancers due to its specific cellular origin. This heterogeneous and rare tumor has a high prevalence of advanced diseases, making it crucial to address the limited therapeutic options and enhance complex clinical management. Given the high clinical accessibility of methylation information, we construct the largest MTC methylation cohort to date.

EXPERIMENTAL DESIGN

Seventy-eight fresh-frozen MTC samples constituted our methylation cohort. The comprehensive study process incorporated machine learning, statistical analysis, and in vitro experiments.

RESULTS

Our study pioneered the identification of a three-class clustering system for risk stratification, exhibiting pronounced epigenomic heterogeneity. The elevated overall methylation status in MTC-B, combined with the "mutual exclusivity" of hypomethylated sites displayed by MTC-A and MTC-C, distinctively characterized the MTC-specific methylation pattern. Integrating with the transcriptome, we further depicted the features of these three clusters to scrutinize biological properties. Several MTC-specific aberrant DNA methylation events were emphasized in our study. NNAT expression was found to be notably reduced in poor-prognostic MTC-C, with its promoter region overlapping with an upregulated differentially methylated region. In vitro experiments further affirmed NNAT's therapeutic potential. Moreover, we built an elastic-net logistic regression model with a relatively high AUC encompassing 68 probes, intended for future validation and systematic clinical application.

CONCLUSIONS

Conducting research on diseases with low incidence poses significant challenges, and we provide a robust resource and comprehensive research framework to assist in ongoing MTC case inclusion and facilitate in-depth dissection of its molecular biological features.

Reversal of dual epigenetic repression of non-canonical Wnt-5a normalises diabetic corneal epithelial wound healing and stem cells

AIMS/HYPOTHESIS

Diabetes is associated with epigenetic modifications including DNA methylation and miRNA changes. Diabetic complications in the cornea can cause persistent epithelial defects and impaired wound healing due to limbal epithelial stem cell (LESC) dysfunction. In this study, we aimed to uncover epigenetic alterations in diabetic vs non-diabetic human limbal epithelial cells (LEC) enriched in LESC and identify new diabetic markers that can be targeted for therapy to normalise corneal epithelial wound healing and stem cell expression.

METHODS

Human LEC were isolated, or organ-cultured corneas were obtained, from autopsy eyes from non-diabetic (59.87±20.89 years) and diabetic (71.93±9.29 years) donors. The groups were not statistically different in age. DNA was extracted from LEC for methylation analysis using Illumina Infinium 850K MethylationEPIC BeadChip and protein was extracted for Wnt phospho array analysis. Wound healing was studied using a scratch assay in LEC or 1-heptanol wounds in organ-cultured corneas. Organ-cultured corneas and LEC were transfected with WNT5A siRNA, miR-203a mimic or miR-203a inhibitor or were treated with recombinant Wnt-5a (200 ng/ml), DNA methylation inhibitor zebularine (1-20 µmol/l) or biodegradable nanobioconjugates (NBCs) based on polymalic acid scaffold containing antisense oligonucleotide (AON) to miR-203a or a control scrambled AON (15-20 µmol/l).

RESULTS

There was significant differential DNA methylation between diabetic and non-diabetic LEC. WNT5A promoter was hypermethylated in diabetic LEC accompanied with markedly decreased Wnt-5a protein. Treatment of diabetic LEC and organ-cultured corneas with exogenous Wnt-5a accelerated wound healing by 1.4-fold (p<0.05) and 37% (p<0.05), respectively, and increased LESC and diabetic marker expression. Wnt-5a treatment in diabetic LEC increased the phosphorylation of members of the Ca2+-dependent non-canonical pathway (phospholipase Cγ1 and protein kinase Cβ; by 1.15-fold [p<0.05] and 1.36-fold [p<0.05], respectively). In diabetic LEC, zebularine treatment increased the levels of Wnt-5a by 1.37-fold (p<0.01)and stimulated wound healing in a dose-dependent manner with a 1.6-fold (p<0.01) increase by 24 h. Moreover, zebularine also improved wound healing by 30% (p<0.01) in diabetic organ-cultured corneas and increased LESC and diabetic marker expression. Transfection of these cells with WNT5A siRNA abrogated wound healing stimulation by zebularine, suggesting that its effect was primarily due to inhibition of WNT5A hypermethylation. Treatment of diabetic LEC and organ-cultured corneas with NBC enhanced wound healing by 1.4-fold (p<0.01) and 23.3% (p<0.05), respectively, with increased expression of LESC and diabetic markers.

CONCLUSIONS/INTERPRETATION

We provide the first account of epigenetic changes in diabetic corneas including dual inhibition of WNT5A by DNA methylation and miRNA action. Overall, Wnt-5a is a new corneal epithelial wound healing stimulator that can be targeted to improve wound healing and stem cells in the diabetic cornea.

DATA AVAILABILITY

The DNA methylation dataset is available from the public GEO repository under accession no. GSE229328 ( https://www.ncbi.nlm.nih.gov/geo/query/acc.cgi?acc=GSE229328 ).

Distinct Regulation of EZH2 and its Repressive H3K27me3 Mark in Polyomavirus-Positive and -Negative Merkel Cell Carcinoma

Merkel cell carcinoma (MCC) is an aggressive skin cancer for which Merkel cell polyomavirus integration and expression of viral oncogenes small T and Large T have been identified as major oncogenic determinants. Recently, a component of the PRC2 complex, the histone methyltransferase enhancer of zeste homolog 2 (EZH2) that induces H3K27 trimethylation as a repressive mark has been proposed as a potential therapeutic target in MCC. Because divergent results have been reported for the levels of EZH2 and trimethylation of lysine 27 on histone 3, we analyzed these factors in a large MCC cohort to identify the molecular determinants of EZH2 activity in MCC and to establish MCC cell lines' sensitivity to EZH2 inhibitors. Immunohistochemical expression of EZH2 was observed in 92% of MCC tumors (156 of 170), with higher expression levels in virus-positive than virus-negative tumors (P = 0.026). For the latter, we showed overexpression of EZHIP, a negative regulator of the PRC2 complex. In vitro, ectopic expression of the large T antigen in fibroblasts led to the induction of EZH2 expression, whereas the knockdown of T antigens in MCC cell lines resulted in decreased EZH2 expression. EZH2 inhibition led to selective cytotoxicity on virus-positive MCC cell lines. This study highlights the distinct mechanisms of EZH2 induction between virus-negative and -positive MCC.

Chromosomal Aberrations Accumulate during Metastasis of Virus-Negative Merkel Cell Carcinoma

Merkel cell carcinoma is a rare, aggressive skin tumor initiated by polyomavirus integration or UV light DNA damage. In New Zealand, there is a propensity toward the UV-driven form (31 of 107, 29% virus positive). Using archival formalin-fixed, paraffin-embedded tissues, we report targeted DNA sequencing covering 246 cancer genes on 71 tumor tissues and 38 nonmalignant tissues from 37 individuals, with 33 of 37 being negative for the virus. Somatic variants of New Zealand virus-negative Merkel cell carcinomas partially overlapped with those reported overseas, including TP53 variants in all tumors and RB1, LRP1B, NOTCH1, and EPHA3/7 variants each found in over half of the cohort. Variants in genes not analyzed or reported in previous studies were also found. Cataloging variants in TP53 and RB1 from published datasets revealed a broad distribution across these genes. Chr 1p gain and Chr 3p loss were identified in around 50% of New Zealand virus-negative Merkel cell carcinomas, and RB1 loss of heterozygosity was found in 90% of cases. Copy number variants accumulate in most metastases. Virus-negative Merkel cell carcinomas have complex combinations of somatic DNA-sequence variants and copy number variants. They likely carry the small genomic changes permissive for metastasis from early tumor development; however, chromosomal alterations may contribute to driving metastatic progression.

Transcriptomic analyses reveal three distinct molecular subgroups of Merkel cell carcinoma with differing prognoses

Merkel cell carcinoma (MCC) is a cutaneous neuroendocrine malignancy with a poor prognosis and an unknown cell of origin. Proffered cells of origin include epithelial stem cells of the hair follicle or interfollicular epidermis, dermal stem cells and pro/pre- or pre-B cells. MCC has also been proposed to have more than one cell of origin and indeed to represent more than one type of carcinoma, currently grouped together due to phenotypic similarities. We explored the heterogeneous nature of MCC by studying the most variably expressed genes with the goal of identifying gene expression patterns that are either clinically relevant or have implications regarding the cell(s) of origin. We performed RNA sequencing on primary tumor samples from 102 patients and identified the top 200 most variably expressed genes. These genes and the tumor samples were hierarchically clustered based on their expression. The functions of three gene clusters exhibiting clearly divergent expression between samples were studied by cross-referencing the lists of genes with online databases. High expression of a gene cluster related to embryonic developmental processes and low expression of a gene cluster related to neuroendocrine processes distinguished Merkel cell polyomavirus (MCPyV)-negative tumors from MCPyV-positive tumors. Furthermore, two prognostically relevant subgroups of MCPyV-positive MCC were identified based on dichotomic expression of genes related to epidermal structures and processes. We identified three distinct molecular subgroups of MCC with prognostic relevance. We propose that the dichotomic expression of epidermis-related genes might reflect both an epidermal and a nonepidermal origin for MCPyV-positive MCC.

Neurocan expression associates with better survival and viral positivity in Merkel cell carcinoma

Merkel cell carcinoma (MCC) is a rare cutaneous neuroendocrine carcinoma that is frequently divided into Merkel cell polyomavirus negative and positive tumors due their distinct genomic and transcriptomic profiles, and disease outcomes. Although some prognostic factors in MCC are known, tumorigenic pathways, which that explain outcome differences in MCC are not fully understood. We investigated transcriptomes of 110 tissue samples of a formalin-fixed, paraffin-embedded MCC series by RNA sequencing to identify genes showing a bimodal expression pattern and predicting outcome in cancer and that potentially could play a role in tumorigenesis. We discovered 19 genes among which IGHM, IGKC, NCAN, OTOF, and USH2A were associated also with overall survival (all p-values < 0.05). From these genes, NCAN (neurocan) expression was detected in all 144 MCC samples by immunohistochemistry. Increased NCAN expression was associated with presence of Merkel cell polyomavirus DNA (p = 0.001) and viral large T antigen expression in tumor tissue (p = 0.004) and with improved MCC-specific survival (p = 0.027) and overall survival (p = 0.034). We conclude that NCAN expression is common in MCC, and further studies are warranted to investigate its role in MCC tumorigenesis.

Ross J, Scott Herron G. Genetic and Epigenetic Influences on Cutaneous Cellular Senescence

Skin is the largest human organ system, and its protective function is critical to survival. The epithelial, dermal, and subcutaneous compartments are heterogeneous mixtures of cell types, yet they all display age-related skin dysfunction through the accumulation of an altered phenotypic cellular state called senescence. Cellular senescence is triggered by complex and dynamic genetic and epigenetic processes. A senescence steady state is achieved in different cell types under various and overlapping conditions of chronological age, toxic injury, oxidative stress, replicative exhaustion, DNA damage, metabolic dysfunction, and chromosomal structural changes. These inputs lead to outputs of cell-cycle withdrawal and the appearance of a senescence-associated secretory phenotype, both of which accumulate as tissue pathology observed clinically in aged skin. This review details the influence of genetic and epigenetic factors that converge on normal cutaneous cellular processes to create the senescent state, thereby dictating the response of the skin to the forces of both intrinsic and extrinsic aging. From this work, it is clear that no single biomarker or process leads to senescence, but that it is a convergence of factors resulting in an overt aging phenotype.

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.

Epigenetic Dysregulations in Merkel Cell Polyomavirus-Driven Merkel Cell Carcinoma

Merkel cell polyomavirus (MCPyV) is a small DNA virus with oncogenic potential. MCPyV is the causative agent of Merkel Cell Carcinoma (MCC), a rare but aggressive tumor of the skin. The role of epigenetic mechanisms, such as histone posttranslational modifications (HPTMs), DNA methylation, and microRNA (miRNA) regulation on MCPyV-driven MCC has recently been highlighted. In this review, we aim to describe and discuss the latest insights into HPTMs, DNA methylation, and miRNA regulation, as well as their regulative factors in the context of MCPyV-driven MCC, to provide an overview of current findings on how MCPyV is involved in the dysregulation of these epigenetic processes. The current state of the art is also described as far as potentially using epigenetic dysregulations and related factors as diagnostic and prognostic tools is concerned, in addition to targets for MCPyV-driven MCC therapy. Growing evidence suggests that the dysregulation of HPTMs, DNA methylation, and miRNA pathways plays a role in MCPyV-driven MCC etiopathogenesis, which, therefore, may potentially be clinically significant for this deadly tumor. A deeper understanding of these mechanisms and related factors may improve diagnosis, prognosis, and therapy for MCPyV-driven MCC.