The CD8+ T cell tolerance checkpoint triggers a distinct differentiation state defined by protein translation defects.

Peripheral CD8+ T cell tolerance is a checkpoint in both autoimmune disease and anti-cancer immunity. Despite its importance, the relationship between tolerance-induced states and other CD8+ T cell differentiation states remains unclear. Using flow cytometric phenotyping, single-cell RNA sequencing (scRNA-seq), and chromatin accessibility profiling, we demonstrated that in vivo peripheral tolerance to a self-antigen triggered a fundamentally distinct differentiation state separate from exhaustion, memory, and functional effector cells but analogous to cells defectively primed against tumors. Tolerant cells diverged early and progressively from effector cells, adopting a transcriptionally and epigenetically distinct state within 60 h of antigen encounter. Breaching tolerance required the synergistic actions of strong T cell receptor (TCR) signaling and inflammation, which cooperatively induced gene modules that enhanced protein translation. Weak TCR signaling during bystander infection failed to breach tolerance due to the uncoupling of effector gene expression from protein translation. Thus, tolerance engages a distinct differentiation trajectory enforced by protein translation defects.

INX-315, a Selective CDK2 Inhibitor, Induces Cell Cycle Arrest and Senescence in Solid Tumors.

Cyclin-dependent kinase 2 (CDK2) is thought to play an important role in driving proliferation of certain cancers, including those harboring CCNE1 amplification and breast cancers that have acquired resistance to CDK4/6 inhibitors (CDK4/6i). The precise impact of pharmacologic inhibition of CDK2 is not known due to the lack of selective CDK2 inhibitors. Here we describe INX-315, a novel and potent CDK2 inhibitor with high selectivity over other CDK family members. Using cell-based assays, patient-derived xenografts (PDX), and transgenic mouse models, we show that INX-315 (i) promotes retinoblastoma protein hypophosphorylation and therapy-induced senescence (TIS) in CCNE1-amplified tumors, leading to durable control of tumor growth; (ii) overcomes breast cancer resistance to CDK4/6i, restoring cell cycle control while reinstating the chromatin architecture of CDK4/6i-induced TIS; and (iii) delays the onset of CDK4/6i resistance in breast cancer by driving deeper suppression of E2F targets. Our results support the clinical development of selective CDK2 inhibitors. SIGNIFICANCE: INX-315 is a novel, selective inhibitor of CDK2. Our preclinical studies demonstrate activity for INX-315 in both CCNE1-amplified cancers and CDK4/6i-resistant breast cancer. In each case, CDK2 inhibition induces cell cycle arrest and a phenotype resembling cellular senescence. Our data support the development of selective CDK2 inhibitors in clinical trials. See related commentary by Watts and Spencer, p. 386. This article is featured in Selected Articles from This Issue, p. 384.

An epigenetic signature of advanced colorectal cancer metastasis.

Colorectal cancer (CRC) is a leading cause of morbidity and mortality worldwide. The majority of CRC deaths are caused by tumor metastasis, even following treatment. There is strong evidence for epigenetic changes, such as DNA methylation, accompanying CRC metastasis and poorer patient survival. Earlier detection and a better understanding of molecular drivers for CRC metastasis are of critical clinical importance. Here, we identify a signature of advanced CRC metastasis by performing whole genome-scale DNA methylation and full transcriptome analyses of paired primary cancers and liver metastases from CRC patients. We observed striking methylation differences between primary and metastatic pairs. A subset of loci showed coordinated methylation-expression changes, suggesting these are potentially epigenetic drivers that control the expression of critical genes in the metastatic cascade. The identification of CRC epigenomic markers of metastasis has the potential to enable better outcome prediction and lead to the discovery of new therapeutic targets.

Repurposing Melanoma Chemotherapy to Activate Inflammasomes in the Treatment of BRAF/MAPK Inhibitor Resistant Melanoma.

The development of resistance to treatments of melanoma is commonly associated with an upregulation of the MAPK pathway and the development of an undifferentiated state. Previous studies have suggested that melanoma with these resistance characteristics may be susceptible to innate death mechanisms such as pyroptosis triggered by the activation of inflammasomes. In this study, we have taken cell lines from patients before and after the development of resistance to BRAF V600 inhibitors and exposed the resistant melanoma to temozolomide (a commonly used chemotherapy) with and without chloroquine to inhibit autophagy. It was found that melanoma with an inflammatory undifferentiated state appeared susceptible to this combination when tested in vitro and in vivo against xenografts in nonobese diabetic scid gamma mice. Translation of the latter results into patients would promise durable responses in patients treated by the combination. The inflammasome and death mechanism involved appeared to vary between melanoma and involved either AIM2 or NLRP3 inflammasomes and gasdermin D or E. These preliminary studies have raised questions as to the selectivity for different inflammasomes in different melanoma and their selective targeting by chemotherapy. They also question whether the inflammatory state of melanoma may be used as biomarkers to select patients for inflammasome-targeted therapy.

Transcriptional Reprogramming and Constitutive PD-L1 Expression in Melanoma Are Associated with Dedifferentiation and Activation of Interferon and Tumour Necrosis Factor Signalling Pathways.

Melanoma is the most aggressive type of skin cancer, with increasing incidence worldwide. Advances in targeted therapy and immunotherapy have improved the survival of melanoma patients experiencing recurrent disease, but unfortunately treatment resistance frequently reduces patient survival. Resistance to targeted therapy is associated with transcriptomic changes and has also been shown to be accompanied by increased expression of programmed death ligand 1 (PD-L1), a potent inhibitor of immune response. Intrinsic upregulation of PD-L1 is associated with genome-wide DNA hypomethylation and widespread alterations in gene expression in melanoma cell lines. However, an in-depth analysis of the transcriptomic landscape of melanoma cells with intrinsically upregulated PD-L1 expression is lacking. To determine the transcriptomic landscape of intrinsically upregulated PD-L1 expression in melanoma, we investigated transcriptomes in melanomas with constitutive versus inducible PD-L1 expression (referred to as PD-L1CON and PD-L1IND). RNA-Seq analysis was performed on seven PD-L1CON melanoma cell lines and ten melanoma cell lines with low inducible PD-L1IND expression. We observed that PD-L1CON melanoma cells had a reprogrammed transcriptome with a characteristic pattern of dedifferentiated gene expression, together with active interferon (IFN) and tumour necrosis factor (TNF) signalling pathways. Furthermore, we identified key transcription factors that were also differentially expressed in PD-L1CON versus PD-L1IND melanoma cell lines. Overall, our studies describe transcriptomic reprogramming of melanomas with PD-L1CON expression.

Marked Global DNA Hypomethylation Is Associated with Constitutive PD-L1 Expression in Melanoma.

Constitutive expression of the immune checkpoint, PD-L1, inhibits anti-tumor immune responses in cancer, although the factors involved in PD-L1 regulation are poorly understood. Here we show that loss of global DNA methylation, particularly in intergenic regions and repeat elements, is associated with constitutive (PD-L1CON), versus inducible (PD-L1IND), PD-L1 expression in melanoma cell lines. We further show this is accompanied by transcriptomic up-regulation. De novo epigenetic regulators (e.g., DNMT3A) are strongly correlated with PD-L1 expression and methylome status. Accordingly, decitabine-mediated inhibition of global methylation in melanoma cells leads to increased PD-L1 expression. Moreover, viral mimicry and immune response genes are highly expressed in lymphocyte-negative plus PD-L1-positive melanomas, versus PD-L1-negative melanomas in The Cancer Genome Atlas (TCGA). In summary, using integrated genomic analysis we identified that global DNA methylation influences PD-L1 expression in melanoma, and hence melanoma's ability to evade anti-tumor immune responses. These results have implications for combining epigenetic therapy with immunotherapy.

A Guide for Designing and Analyzing RNA-Seq Data.

The identity of a cell or an organism is at least in part defined by its gene expression and therefore analyzing gene expression remains one of the most frequently performed experimental techniques in molecular biology. The development of the RNA-Sequencing (RNA-Seq) method allows an unprecedented opportunity to analyze expression of protein-coding, noncoding RNA and also de novo transcript assembly of a new species or organism. However, the planning and design of RNA-Seq experiments has important implications for addressing the desired biological question and maximizing the value of the data obtained. In addition, RNA-Seq generates a huge volume of data and accurate analysis of this data involves several different steps and choices of tools. This can be challenging and overwhelming, especially for bench scientists. In this chapter, we describe an entire workflow for performing RNA-Seq experiments. We describe critical aspects of wet lab experiments such as RNA isolation, library preparation and the initial design of an experiment. Further, we provide a step-by-step description of the bioinformatics workflow for different steps involved in RNA-Seq data analysis. This includes power calculations, setting up a computational environment, acquisition and processing of publicly available data if desired, quality control measures, preprocessing steps for the raw data, differential expression analysis, and data visualization. We particularly mention important considerations for each step to provide a guide for designing and analyzing RNA-Seq data.

The Slow Cycling Phenotype: A Growing Problem for Treatment Resistance in Melanoma.

Treatment resistance in metastatic melanoma is a longstanding issue. Current targeted therapy regimes in melanoma largely target the proliferating cancer population, leaving slow-cycling cancer cells undamaged. Consequently, slow-cycling cells are enriched upon drug therapy and can remain in the body for years until acquiring proliferative potential that triggers cancer relapse. Here we overview the molecular mechanisms of slow-cycling cells that underlie treatment resistance in melanoma. Three main areas of molecular reprogramming are discussed that mediate slow cycling and treatment resistance. First, a low microphthalmia-associated transcription factor (MITF) dedifferentiated state activates various signaling pathways. This includes WNT5A, EGFR, as well as other signaling activators, such as AXL and NF-κB. Second, the chromatin-remodeling factor Jumonji/ARID domain-containing protein 1B (JARID1B, KDM5B) orchestrates and maintains slow cycling and treatment resistance in a small subpopulation of melanoma cells. Finally, a shift in metabolic state toward oxidative phosphorylation has been demonstrated to regulate treatment resistance in slow-cycling cells. Elucidation of the underlying processes of slow cycling and its utilization by melanoma cells may reveal new vulnerable characteristics as therapeutic targets. Moreover, combining current therapies with targeting slow-cycling subpopulations of melanoma cells may allow for more durable and greater treatment responses. Mol Cancer Ther; 16(6); 1002-9. ©2017 AACR.

Comparative assessment of DNA methylation patterns between reduced representation bisulfite sequencing and Sequenom EpiTyper methylation analysis.

AIM: Validation of sequencing-based DNA methylation data is an important step for meaningful translation of findings. However, there has been limited assessment of different platforms to validate methylation data from next generation sequencing. METHODS: We performed a comparative methylation analysis between the genome-wide platform of reduced representation bisulfite sequencing with a targeted, Sequenom EpiTyper platform (four genes were analyzed in 15 cell lines covering 52 CpG sites). RESULTS: We show that the accuracy of validation substantially improves if results from multiple and adjacent CpG sites are combined rather than at single CpG sites. We demonstrate increased read number improves accuracy of reduced representation bisulfite sequencing results. Further, by using series of replicates, we document variation in samples analyzed by Sequenom EpiTyper, indicating the importance of including replicates to increase precision. CONCLUSION: The results reveal potential sources of bias and provide a guideline for refining study design for DNA methylation analysis.

Genome-wide methylation sequencing of paired primary and metastatic cell lines identifies common DNA methylation changes and a role for EBF3 as a candidate epigenetic driver of melanoma metastasis.

Epigenetic alterations are increasingly implicated in metastasis, whereas very few genetic mutations have been identified as authentic drivers of cancer metastasis. Yet, to date, few studies have identified metastasis-related epigenetic drivers, in part because a framework for identifying driver epigenetic changes in metastasis has not been established. Using reduced representation bisulfite sequencing (RRBS), we mapped genome-wide DNA methylation patterns in three cutaneous primary and metastatic melanoma cell line pairs to identify metastasis-related epigenetic drivers. Globally, metastatic melanoma cell lines were hypomethylated compared to the matched primary melanoma cell lines. Using whole genome RRBS we identified 75 shared (10 hyper- and 65 hypomethylated) differentially methylated fragments (DMFs), which were associated with 68 genes showing significant methylation differences. One gene, Early B Cell Factor 3 (EBF3), exhibited promoter hypermethylation in metastatic cell lines, and was validated with bisulfite sequencing and in two publicly available independent melanoma cohorts (n = 40 and 458 melanomas, respectively). We found that hypermethylation of the EBF3 promoter was associated with increased EBF3 mRNA levels in metastatic melanomas and subsequent inhibition of DNA methylation reduced EBF3 expression. RNAi-mediated knockdown of EBF3 mRNA levels decreased proliferation, migration and invasion in primary and metastatic melanoma cell lines. Overall, we have identified numerous epigenetic changes characterising metastatic melanoma cell lines, including EBF3-induced aggressive phenotypic behaviour with elevated EBF3 expression in metastatic melanoma, suggesting that EBF3 promoter hypermethylation may be a candidate epigenetic driver of metastasis.

NRAS and EPHB6 mutation rates differ in metastatic melanomas of patients in the North Island versus South Island of New Zealand.

Melanoma, the most aggressive skin cancer type, is responsible for 75% of skin cancer related deaths worldwide. Given that New Zealand (NZ) has the world's highest melanoma incidence, we sought to determine the frequency of mutations in NZ melanomas in recurrently mutated genes. NZ melanomas were from localities distributed between North (35°S-42°S) and South Islands (41°S-47°S). A total of 529 melanomas were analyzed for BRAF exon 15 mutations by Sanger sequencing, and also by Sequenom MelaCarta MassARRAY. While, a relatively low incidence of BRAFV600E mutations (23.4%) was observed overall in NZ melanomas, the incidence of NRAS mutations in South Island melanomas was high compared to North Island melanomas (38.3% vs. 21.9%, P=0.0005), and to The Cancer Genome Atlas database (TCGA) (38.3% vs. 22%, P=0.0004). In contrast, the incidence of EPHB6G404S mutations was 0% in South Island melanomas, and was 7.8% in North Island (P=0.0002). Overall, these data suggest that melanomas from geographically different regions in NZ have markedly different mutation frequencies, in particular in the NRAS and EPHB6 genes, when compared to TCGA or other populations. These data have implications for the causation and treatment of malignant melanoma in NZ.

The rs11515 Polymorphism Is More Frequent and Associated With Aggressive Breast Tumors with Increased ANRIL and Decreased p16 (INK4a) Expression.

Chromosome position 9p21 encodes three-tumor suppressors p16(INK4a), p14(ARF), and p15(INK4b) and the long non-coding RNA ANRIL (antisense non-coding RNA in the INK4 locus). The rs11515 single-nucleotide polymorphism in the p16 (INK4a) /p14 (ARF) 3'-untranslated region is associated with glioblastoma, melanoma, and other cancers. This study investigated the frequency and effect of rs11515 genotypes in breast cancer. Genomic DNA samples from 400 women (200 with and 200 without a diagnosis of breast cancer) were genotyped for the rs11515 major (C) and minor (G) alleles. The rs11515 polymorphism was also investigated in 108 heart tissues to test for tissue-specific effects. Four 9p21 transcripts, p16 (INK4a) , p14 (ARF) , p15 (INK4b) , and ANRIL were measured in breast tumors and myocardium using quantitative PCR. Heterozygotes (CG genotype) were more frequent in women with breast cancer compared to the control population (P = 0.0039). In those with breast cancer, the CG genotype was associated with an older age (P = 0.016) and increased lymph node involvement (P = 0.007) compared to homozygotes for the major allele (CC genotype). In breast tumors, the CG genotype had higher ANRIL (P = 0.031) and lower p16 (INK4a) (P = 0.006) expression compared to the CC genotype. The CG genotype was not associated with altered 9p21 transcripts in heart tissue. In breast cancer, the rs11515 CG genotype is more frequent and associated with a more aggressive tumor that could be due to increased ANRIL and reduced p16 (INK4a) expression. The absence of association between rs11515 genotypes and 9p21 transcripts in heart tissue suggests this polymorphism has tissue- or disease-specific functions.

MITF and PAX3 Play Distinct Roles in Melanoma Cell Migration; Outline of a "Genetic Switch" Theory Involving MITF and PAX3 in Proliferative and Invasive Phenotypes of Melanoma.

Melanoma is a very aggressive neoplasm with a propensity to undergo progression and invasion early in its evolution. The molecular pathways underpinning invasion in melanoma are now just beginning to be elucidated, but a clear understanding of the transition from non-invasive to invasive melanoma cells remains elusive. Microphthalmia-associated transcription factor (MITF), is thought to be a central player in melanoma biology, and it controls many aspects of the phenotypic expression of the melanocytic lineage. However, recently the paired box transcription factor PAX3 was shown to transcriptionally activate POU3F2/BRN2, leading to direct repression of MITF expression. Here we present a theory to explain melanoma phenotype switching and discuss the predictions that this theory makes. One prediction is that independent and opposing roles for MITF and PAX3 in melanoma would be expected, and we present empirical evidence supporting this: in melanoma tissues PAX3 expression occurs independently of MITF, and PAX3 does not play a key role in melanoma cell proliferation. Furthermore, we show that knockdown of PAX3 inhibits cell migration in a group of "lower MITF" melanoma cell lines, while knockdown of MITF promotes cell migration in a complementary "higher MITF" group of melanoma cell lines. Moreover, the morphological effects of knocking down PAX3 versus MITF in melanoma cells were found to differ. While these data support the notion of independent roles for MITF and PAX3, additional experiments are required to provide robust examination of the proposed genetic switch theory. Only upon clear delineation of the mechanisms associated with progression and invasion of melanoma cells will successful treatments for invasive melanoma be developed.

The CDKN2A G500 allele is more frequent in GBM patients with no defined telomere maintenance mechanism tumors and is associated with poorer survival.

Prognostic markers for glioblastoma multiforme (GBM) are important for patient management. Recent advances have identified prognostic markers for GBMs that use telomerase or the alternative lengthening of telomeres (ALT) mechanism for telomere maintenance. Approximately 40% of GBMs have no defined telomere maintenance mechanism (NDTMM), with a mixed survival for affected individuals. This study examined genetic variants in the cyclin-dependent kinase inhibitor 2A (CDKN2A) gene that encodes the p16(INK4a) and p14(ARF) tumor suppressors, and the isocitrate dehydrogenase 1 (IDH1) gene as potential markers of survival for 40 individuals with NDTMM GBMs (telomerase negative and ALT negative by standard assays), 50 individuals with telomerase, and 17 individuals with ALT positive tumors. The analysis of CDKN2A showed NDTMM GBMs had an increased minor allele frequency for the C500G (rs11515) polymorphism compared to those with telomerase and ALT positive GBMs (p = 0.002). Patients with the G500 allele had reduced survival that was independent of age, extent of surgery, and treatment. In the NDTMM group G500 allele carriers had increased loss of CDKN2A gene dosage compared to C500 homozygotes. An analysis of IDH1 mutations showed the R132H mutation was associated with ALT positive tumors, and was largely absent in NDTMM and telomerase positive tumors. In the ALT positive tumors cohort, IDH1 mutations were associated with a younger age for the affected individual. In conclusion, the G500 CDKN2A allele was associated with NDTMM GBMs from older individuals with poorer survival. Mutations in IDH1 were not associated with NDTMM GBMs, and instead were a marker for ALT positive tumors in younger individuals.