Epigenetic plasticity potentiates a rapid cyclical shift to and from an aggressive cancer phenotype.

Highly tumorigenic, drug-resistant cancer stem-like cells drive cancer progression. These aggressive cells can arise repeatedly from bulk tumor cells independently of mutational events, suggesting an epigenetic mechanism. To test this possibility, we studied bladder cancer cells as they cyclically shifted to and from a cancer stem-like phenotype, and we discovered that these two states exhibit distinct DNA methylation and chromatin accessibility. Most differential chromatin accessibility was independent of methylation and affected the expression of driver genes such as E2F3, a cell cycle regulator associated with aggressive bladder cancer. Cancer stem-like cells exhibited increased E2F3 promoter accessibility and increased E2F3 expression that drove cell migration, invasiveness and drug resistance. Epigenetic interference using a DNA methylation inhibitor blocked the transition to a cancer stem-like state and reduced E2F3 expression. Our findings indicate that epigenetic plasticity plays a key role in the transition to and from an aggressive, drug-resistant phenotype.

SETD2 Deficiency Confers Sensitivity to Dual Inhibition of DNA Methylation and PARP in Kidney Cancer.

SETD2 deficiency alters the epigenetic landscape by causing depletion of H3K36me3 and plays an important role in diverse forms of cancer, most notably in aggressive and metastatic clear-cell renal cell carcinomas (ccRCC). Development of an effective treatment scheme targeting SETD2-compromised cancer is urgently needed. Considering that SETD2 is involved in DNA methylation and DNA repair, a combination treatment approach using DNA hypomethylating agents (HMA) and PARP inhibitors (PARPi) could have strong antitumor activity in SETD2-deficient kidney cancer. We tested the effects of the DNA HMA 5-aza-2'-dexoxydytidine (DAC), the PARPi talazoparib (BMN-673), and both in combination in human ccRCC models with or without SETD2 deficiency. The combination treatment of DAC and BMN-673 synergistically increased cytotoxicity in vitro in SETD2-deficient ccRCC cell lines but not in SETD2-proficient cell lines. DAC and BMN-673 led to apoptotic induction, increased DNA damage, insufficient DNA damage repair, and increased genomic instability. Furthermore, the combination treatment elevated immune responses, upregulated STING, and enhanced viral mimicry by activating transposable elements. Finally, the combination effectively suppressed the growth of SETD2-deficient ccRCC in in vivo mouse models. Together, these findings indicate that combining HMA and PARPi is a promising potential therapeutic strategy for treating SETD2-compromised ccRCC. SIGNIFICANCE: SETD2 deficiency creates a vulnerable epigenetic status that is targetable using a DNA hypomethylating agent and PARP inhibitor combination to suppress renal cell carcinoma, identifying a precision medicine-based approach for SETD2-compromised cancers.

RNA mis-splicing drives viral mimicry response after DNMTi therapy in SETD2-mutant kidney cancer.

Tumors with mutations in chromatin regulators present attractive targets for DNA hypomethylating agent 5-aza-2'-deoxycytidine (DAC) therapy, which further disrupts cancer cells' epigenomic fidelity and reactivates transposable element (TE) expression to drive viral mimicry responses. SETD2 encodes a histone methyltransferase (H3K36me3) and is prevalently mutated in advanced kidney cancers. Here, we show that SETD2-mutant kidney cancer cells are especially sensitive in vitro and in vivo to DAC treatment. We find that the viral mimicry response are direct consequences of mis-splicing events, such as exon inclusions or extensions, triggered by DAC treatment in an SETD2-loss context. Comprehensive epigenomic analysis reveals H3K9me3 deposition, rather than DNA methylation dynamics, across intronic TEs might contribute to elevated mis-splicing rates. Through epigenomic and transcriptomic analyses, we show that SETD2-deficient kidney cancers are prone to mis-splicing, which can be therapeutically exacerbated with DAC treatment to increase viral mimicry activation and provide synergy with combinatorial immunotherapy approaches.

Genes regulated by DNA methylation are involved in distinct phenotypes during melanoma progression and are prognostic factors for patients.

In addition to mutations, epigenetic alterations are important contributors to malignant transformation and tumor progression. The aim of this work was to identify epigenetic events in which promoter or gene body DNA methylation induces gene expression changes that drive melanocyte malignant transformation and metastasis. We previously developed a linear mouse model of melanoma progression consisting of spontaneously immortalized melanocytes, premalignant melanocytes, a nonmetastatic tumorigenic, and a metastatic cell line. Here, through the integrative analysis of methylome and transcriptome data, we identified the relationship between promoter and/or gene body DNA methylation alterations and gene expression in early, intermediate, and late stages of melanoma progression. We identified adenylate cyclase type 3 (Adcy3) and inositol polyphosphate 4-phosphatase type II (Inpp4b), which affect tumor growth and metastatic potential, respectively. Importantly, the gene expression and DNA methylation profiles found in this murine model of melanoma progression were correlated with available clinical data from large population-based primary melanoma cohorts, revealing potential prognostic markers.

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

BACKGROUND: Merkel cell carcinoma (MCC) is a rare but aggressive skin cancer with limited treatment possibilities. Merkel cell tumors display with neuroendocrine features and Merkel cell polyomavirus (MCPyV) infection in the majority (80%) of patients. Although loss of histone H3 lysine 27 trimethylation (H3K27me3) has been shown during MCC tumorigenesis, epigenetic dysregulation has largely been overlooked. METHODS: We conducted global DNA methylation profiling of clinically annotated MCC primary tumors, metastatic skin tumors, metastatic lymph node tumors, paired normal tissues, and two human MCC cell lines using the Illumina Infinium EPIC DNA methylation BeadArray platform. RESULTS: Significant differential DNA methylation patterns across the genome are revealed between the four tissue types, as well as based on MCPyV status. Furthermore, 964 genes directly regulated by promoter or gene body DNA methylation were identified with high enrichment in neuro-related pathways. Finally, our findings suggest that loss of H3K27me3 occupancy in MCC is attributed to KDM6B and EZHIP overexpression as a consequence of promoter DNA hypomethylation. CONCLUSIONS: We have demonstrated specific DNA methylation patterns for primary MCC tumors, metastatic MCCs, and adjacent-normal tissues. We have also identified DNA methylation markers that not only show potential diagnostic or prognostic utility in MCC management, but also correlate with MCC tumorigenesis, MCPyV expression, neuroendocrine features, and H3K27me3 status. The identification of DNA methylation alterations in MCC supports the need for further studies to understand the clinical implications of epigenetic dysregulation and potential therapeutic targets in MCC.

A Novel DNA Methylation Signature as an Independent Prognostic Factor in Muscle-Invasive Bladder Cancer.

BACKGROUND: Muscle-invasive bladder cancer (MIBC) accounts for approximately 20% of all urothelial bladder carcinomas (UBC) at time of diagnosis, and up to 30% of patients with non-muscle invasive UBC will progress to MIBC over time. An increasing body of evidence has revealed a strong correlation between aberrant DNA methylation and tumorigenesis in MIBC. RESULTS: Using The Cancer Genome Atlas (TCGA) molecular data for 413 patients, we described a DNA methylation-based signature as a prognostic factor for overall survival (OS) in MIBC patients. By using a least absolute shrinkage and selection operator (LASSO) model, differentially methylated regions were first identified using multiple criteria followed by survival and LASSO analyses to identify DNA methylation probes related to OS and build a classifier to stratify patients with MIBC. The prognostic value of the classifier, referred to as risk score (RS), was validated in a held-out testing set from the TCGA MIBC cohort. Finally, receiver operating characteristic (ROC) analysis was used to compare the prognostic accuracy of the models built with RS alone, RS plus clinicopathologic features, and clinicopathologic features alone. We found that our seven-probe classifier-based RS stratifies patients into high- and low-risk groups for overall survival (OS) in the testing set (n = 137) (AUC at 3 years, 0.65; AUC at 5 years, 0.65). In addition, RS significantly improved the prognostic model when it was combined with clinical information including age, smoking status, Tumor (T) stage, and Lymph node metastasis (N) stage. CONCLUSIONS: The DNA methylation-based RS can be a useful tool to predict the accuracy of preoperative and/or post-cystectomy models of OS in MIBC patients.

The Roles of Human DNA Methyltransferases and Their Isoforms in Shaping the Epigenome.

A DNA sequence is the hard copy of the human genome and it is a driving force in determining the physiological processes in an organism. Concurrently, the chemical modification of the genome and its related histone proteins is dynamically involved in regulating physiological processes and diseases, which overall constitutes the epigenome network. Among the various forms of epigenetic modifications, DNA methylation at the C-5 position of cytosine in the cytosine⁻guanine (CpG) dinucleotide is one of the most well studied epigenetic modifications. DNA methyltransferases (DNMTs) are a family of enzymes involved in generating and maintaining CpG methylation across the genome. In mammalian systems, DNA methylation is performed by DNMT1 and DNMT3s (DNMT3A and 3B). DNMT1 is predominantly involved in the maintenance of DNA methylation during cell division, while DNMT3s are involved in establishing de novo cytosine methylation and maintenance in both embryonic and somatic cells. In general, all DNMTs require accessory proteins, such as ubiquitin-like containing plant homeodomain (PHD) and really interesting new gene (RING) finger domain 1 (UHRF1) or DNMT3-like (DNMT3L), for their biological function. This review mainly focuses on the role of DNMT3B and its isoforms in de novo methylation and maintenance of DNA methylation, especially with respect to their role as an accessory protein.

Profiling DNA methylation differences between inbred mouse strains on the Illumina Human Infinium MethylationEPIC microarray.

The Illumina Infinium MethylationEPIC provides an efficient platform for profiling DNA methylation in humans at over 850,000 CpGs. Model organisms such as mice do not currently benefit from an equivalent array. Here we used this array to measure DNA methylation in mice. We defined probes targeting conserved regions and performed differential methylation analysis and compared between the array-based assay and affinity-based DNA sequencing of methyl-CpGs (MBD-seq) and reduced representation bisulfite sequencing. Mouse samples consisted of 11 liver DNA from two strains, C57BL/6J (B6) and DBA/2J (D2), that varied widely in age. Linear regression was applied to detect differential methylation. In total, 13,665 probes (1.6% of total probes) aligned to conserved CpGs. Beta-values (ß-value) for these probes showed a distribution similar to that in humans. Overall, there was high concordance in methylation signal between the EPIC array and MBD-seq (Pearson correlation r = 0.70, p-value < 0.0001). However, the EPIC probes had higher quantitative sensitivity at CpGs that are hypo- (ß-value < 0.3) or hypermethylated (ß-value > 0.7). In terms of differential methylation, no EPIC probe detected a significant difference between age groups at a Benjamini-Hochberg threshold of 10%, and the MBD-seq performed better at detecting age-dependent change in methylation. However, the top most significant probe for age (cg13269407; uncorrected p-value = 1.8 x 10-5) is part of the clock CpGs used to estimate the human epigenetic age. For strain, 219 EPIC probes detected significant differential methylation (FDR cutoff 10%) with ~80% CpGs associated with higher methylation in D2. This higher methylation profile in D2 compared to B6 was also replicated by the MBD-seq data. To summarize, we found only a small subset of EPIC probes that target conserved sites. However, for this small subset the array provides a reliable assay of DNA methylation and can be effectively used to measure differential methylation in mice.

Juvenile hormone regulation of female reproduction in the common bed bug, Cimex lectularius.

To begin studies on reproduction in common bed bug, Cimex lectularius, we identified three genes coding for vitellogenin (Vg, a protein required for the reproductive success of insects) and studied their hormonal regulation. RNA interference studied showed that expression of Vg3 gene in the adult females is a prerequisite for successful completion of embryogenesis in the eggs laid by them. Juvenile hormone (JH) receptor, Methoprene-tolerant (Met), steroid receptor coactivator (SRC) and GATAa but not ecdysone receptor (EcR) or its partner, ultraspiracle (USP) are required for expression of Vg genes. Feeding and mating working through Vg, Met, SRC, EcR, and GATAa regulate oocyte development. Knockdown of the expression of Met, SRC, EcR, USP, BR-C (Broad-Complex), TOR (target of rapamycin), and GATAa in female adults resulted in a reduction in the number eggs laid by them. Interestingly, Kruppel homolog 1 (Kr-h1) knockdown in the adult females did not reduce their fecundity but affected the development of embryos in the eggs laid by females injected with Kr-h1 double-stranded RNA. These data suggest that JH functioning through Met and SRC regulate both vitellogenesis and oogenesis in C. lectularius. However, JH does not work through Kr-h1 but may work through transcription factors not yet identified.

Krüppel homolog 1 and E93 mediate Juvenile hormone regulation of metamorphosis in the common bed bug, Cimex lectularius.

The common bed bug is an obligate hematophagous parasite of humans. We studied the regulation of molting and metamorphosis in bed bugs with a goal to identify key players involved. qRT-PCR studies on the expression of genes known to be involved in molting and metamorphosis showed high levels of Krüppel homolog 1 [Kr-h1, a transcription factor that plays key roles in juvenile hormone (JH) action] mRNA in the penultimate nymphal stage (N4). However, low levels of Kr-h1 mRNA were detected in the fifth and last nymphal stage (N5). Knockdown of Kr-h1 in N4 resulted in a precocious development of adult structures. Kr-h1 maintains the immature stage by suppressing E93 (early ecdysone response gene) in N4. E93 expression increases during the N5 in the absence of Kr-h1 and promotes the development of adult structures. Knockdown of E93 in N5 results in the formation of supernumerary nymphs. The role of JH in the suppression of adult structures through interaction with Kr-h1 and E93 was also studied by the topical application of JH analog, methoprene, to N5. Methoprene induced Kr-h1 and suppressed E93 and induced formation of the supernumerary nymph. These data show interactions between Kr-h1, E93 and JH in the regulation of metamorphosis in the bed bugs.

Unique features of a global human ectoparasite identified through sequencing of the bed bug genome.

The bed bug, Cimex lectularius, has re-established itself as a ubiquitous human ectoparasite throughout much of the world during the past two decades. This global resurgence is likely linked to increased international travel and commerce in addition to widespread insecticide resistance. Analyses of the C. lectularius sequenced genome (650 Mb) and 14,220 predicted protein-coding genes provide a comprehensive representation of genes that are linked to traumatic insemination, a reduced chemosensory repertoire of genes related to obligate hematophagy, host-symbiont interactions, and several mechanisms of insecticide resistance. In addition, we document the presence of multiple putative lateral gene transfer events. Genome sequencing and annotation establish a solid foundation for future research on mechanisms of insecticide resistance, human-bed bug and symbiont-bed bug associations, and unique features of bed bug biology that contribute to the unprecedented success of C. lectularius as a human ectoparasite.

Bed bugs evolved unique adaptive strategy to resist pyrethroid insecticides.

Recent advances in genomic and post-genomic technologies have facilitated a genome-wide analysis of the insecticide resistance-associated genes in insects. Through bed bug, Cimex lectularius transcriptome analysis, we identified 14 molecular markers associated with pyrethroid resistance. Our studies revealed that most of the resistance-associated genes functioning in diverse mechanisms are expressed in the epidermal layer of the integument, which could prevent or slow down the toxin from reaching the target sites on nerve cells, where an additional layer of resistance (kdr) is possible. This strategy evolved in bed bugs is based on their unique morphological, physiological and behavioral characteristics and has not been reported in any other insect species. RNA interference-aided knockdown of resistance associated genes showed the relative contribution of each mechanism towards overall resistance development. Understanding the complexity of adaptive strategies employed by bed bugs will help in designing the most effective and sustainable bed bug control methods.