A gene-environment-induced epigenetic program initiates tumorigenesis.

Tissue damage increases the risk of cancer through poorly understood mechanisms1. In mouse models of pancreatic cancer, pancreatitis associated with tissue injury collaborates with activating mutations in the Kras oncogene to markedly accelerate the formation of early neoplastic lesions and, ultimately, adenocarcinoma2,3. Here, by integrating genomics, single-cell chromatin assays and spatiotemporally controlled functional perturbations in autochthonous mouse models, we show that the combination of Kras mutation and tissue damage promotes a unique chromatin state in the pancreatic epithelium that distinguishes neoplastic transformation from normal regeneration and is selected for throughout malignant evolution. This cancer-associated epigenetic state emerges within 48 hours of pancreatic injury, and involves an 'acinar-to-neoplasia' chromatin switch that contributes to the early dysregulation of genes that define human pancreatic cancer. Among the factors that are most rapidly activated after tissue damage in the pre-malignant pancreatic epithelium is the alarmin cytokine interleukin 33, which recapitulates the effects of injury in cooperating with mutant Kras to unleash the epigenetic remodelling program of early neoplasia and neoplastic transformation. Collectively, our study demonstrates how gene-environment interactions can rapidly produce gene-regulatory programs that dictate early neoplastic commitment, and provides a molecular framework for understanding the interplay between genetic and environmental cues in the initiation of cancer.

CCR6 Defines Memory B Cell Precursors in Mouse and Human Germinal Centers, Revealing Light-Zone Location and Predominant Low Antigen Affinity.

Memory B cells (MBCs) and plasma cells (PCs) constitute the two cellular outputs of germinal center (GC) responses that together facilitate long-term humoral immunity. Although expression of the transcription factor BLIMP-1 identifies cells undergoing PC differentiation, no such marker exists for cells committed to the MBC lineage. Here, we report that the chemokine receptor CCR6 uniquely marks MBC precursors in both mouse and human GCs. CCR6+ GC B cells were highly enriched within the GC light zone (LZ), were the most quiescent of all GC B cells, exhibited a cell-surface phenotype and gene expression signature indicative of an MBC transition, and possessed the augmented response characteristics of MBCs. MBC precursors within the GC LZ predominantly possessed a low affinity for antigen but also included cells from within the high-affinity pool. These data indicate a fundamental dichotomy between the processes that drive MBC and PC differentiation during GC responses.

Widespread promoter methylation of synaptic plasticity genes in long-term potentiation in the adult brain in vivo.

BACKGROUND: DNA methylation is a key modulator of gene expression in mammalian development and cellular differentiation, including neurons. To date, the role of DNA modifications in long-term potentiation (LTP) has not been explored. RESULTS: To investigate the occurrence of DNA methylation changes in LTP, we undertook the first detailed study to describe the methylation status of all known LTP-associated genes during LTP induction in the dentate gyrus of live rats. Using a methylated DNA immunoprecipitation (MeDIP)-array, together with previously published matched RNA-seq and public histone modification data, we discover widespread changes in methylation status of LTP-genes. We further show that the expression of many LTP-genes is correlated with their methylation status. We show that these correlated genes are enriched for RNA-processing, active histone marks, and specific transcription factors. These data reveal that the synaptic activity-evoked methylation changes correlates with pre-existing activation of the chromatin landscape. Finally, we show that methylation of Brain-derived neurotrophic factor (Bdnf) CpG-islands correlates with isoform switching from transcripts containing exon IV to exon I. CONCLUSIONS: Together, these data provide the first evidence of widespread regulation of methylation status in LTP-associated genes.

lncRNAdb v2.0: expanding the reference database for functional long noncoding RNAs.

Despite the prevalence of long noncoding RNA (lncRNA) genes in eukaryotic genomes, only a small proportion have been examined for biological function. lncRNAdb, available at http://lncrnadb.org, provides users with a comprehensive, manually curated reference database of 287 eukaryotic lncRNAs that have been described independently in the scientific literature. In addition to capturing a great proportion of the recent literature describing functions for individual lncRNAs, lncRNAdb now offers an improved user interface enabling greater accessibility to sequence information, expression data and the literature. The new features in lncRNAdb include the integration of Illumina Body Atlas expression profiles, nucleotide sequence information, a BLAST search tool and easy export of content via direct download or a REST API. lncRNAdb is now endorsed by RNAcentral and is in compliance with the International Nucleotide Sequence Database Collaboration.

Long noncoding RNAs in cancer: mechanisms of action and technological advancements.

The previous decade has seen long non-coding RNAs (lncRNAs) rise from obscurity to being defined as a category of genetic elements, leaving its mark on the field of cancer biology. With the current number of curated lncRNAs increasing by 10,000 in the last five years, the field is moving from annotation of lncRNA expression in various tumours to understanding their importance in the key cancer signalling networks and characteristic behaviours. Here, we summarize the previously identified as well as recently discovered mechanisms of lncRNA function and their roles in the hallmarks of cancer. Furthermore, we identify novel technologies for investigation of lncRNA properties and their function in carcinogenesis, which will be important for their translation to the clinic as novel biomarkers and therapeutic targets.

CD151 Gene and Protein Expression Provides Independent Prognostic Information for Patients with Adenocarcinoma of the Esophagus and Gastroesophageal Junction Treated by Esophagectomy.

BACKGROUND: Esophageal and gastroesophageal junctional (GEJ) adenocarcinoma is one of the most fatal cancers and has the fastest rising incidence rate of all cancers. Identification of biomarkers is needed to tailor treatments to each patient's tumor biology and prognosis. METHODS: Gene expression profiling was performed in a test cohort of 80 chemoradiotherapy (CRTx)-naïve patients with external validation in a separate cohort of 62 CRTx-naïve patients and 169 patients with advanced-stage disease treated with CRTx. RESULTS: As a novel prognostic biomarker after external validation, CD151 showed promise. Patients exhibiting high levels of CD151 (≥median) had a longer median overall survival than patients with low CD151 tumor levels (median not reached vs. 30.9 months; p = 0.01). This effect persisted in a multivariable Cox-regression model with adjustment for tumor stage [adjusted hazard ratio (aHR), 0.33; 95 % confidence interval (CI), 0.14-0.78; p = 0.01] and was further corroborated through immunohistochemical analysis (aHR, 0.22; 95 % CI, 0.08-0.59; p = 0.003). This effect was not found in the separate cohort of CRTx-exposed patients. CONCLUSION: Tumoral expression levels of CD151 may provide independent prognostic information not gained by conventional staging of patients with esophageal and GEJ adenocarcinoma treated by esophagectomy alone.

Cannabinoid effects on ß amyloid fibril and aggregate formation, neuronal and microglial-activated neurotoxicity in vitro.

Cannabinoid (CB) ligands have demonstrated neuroprotective properties. In this study we compared the effects of a diverse set of CB ligands against ß amyloid-mediated neuronal toxicity and activated microglial-conditioned media-based neurotoxicity in vitro, and compared this with a capacity to directly alter ß amyloid (Aß) fibril or aggregate formation. Neuroblastoma (SH-SY5Y) cells were exposed to Aß1-42 directly or microglial (BV-2 cells) conditioned media activated with lipopolysaccharide (LPS) in the presence of the CB1 receptor-selective agonist ACEA, CB2 receptor-selective agonist JWH-015, phytocannabinoids Δ(9)-THC and cannabidiol (CBD), the endocannabinoids 2-arachidonoyl glycerol (2-AG) and anandamide or putative GPR18/GPR55 ligands O-1602 and abnormal-cannabidiol (Abn-CBD). TNF-α and nitrite production was measured in BV-2 cells to compare activation via LPS or albumin with Aß1-42. Aß1-42 evoked a concentration-dependent loss of cell viability in SH-SY5Y cells but negligible TNF-α and nitrite production in BV-2 cells compared to albumin or LPS. Both albumin and LPS-activated BV-2 conditioned media significantly reduced neuronal cell viability but were directly innocuous to SH-SY5Y cells. Of those CB ligands tested, only 2-AG and CBD were directly protective against Aß-evoked SH-SY5Y cell viability, whereas JWH-015, THC, CBD, Abn-CBD and O-1602 all protected SH-SY5Y cells from BV-2 conditioned media activated via LPS. While CB ligands variably altered the morphology of Aß fibrils and aggregates, there was no clear correlation between effects on Aß morphology and neuroprotective actions. These findings indicate a neuroprotective action of CB ligands via actions at microglial and neuronal cells.

PRC2-Inactivating Mutations in Cancer Enhance Cytotoxic Response to DNMT1-Targeted Therapy via Enhanced Viral Mimicry.

Polycomb repressive complex 2 (PRC2) has oncogenic and tumor-suppressive roles in cancer. There is clinical success of targeting this complex in PRC2-dependent cancers, but an unmet therapeutic need exists in PRC2-loss cancer. PRC2-inactivating mutations are a hallmark feature of high-grade malignant peripheral nerve sheath tumor (MPNST), an aggressive sarcoma with poor prognosis and no effective targeted therapy. Through RNAi screening in MPNST, we found that PRC2 inactivation increases sensitivity to genetic or small-molecule inhibition of DNA methyltransferase 1 (DNMT1), which results in enhanced cytotoxicity and antitumor response. Mechanistically, PRC2 inactivation amplifies DNMT inhibitor-mediated expression of retrotransposons, subsequent viral mimicry response, and robust cell death in part through a protein kinase R (PKR)-dependent double-stranded RNA sensor. Collectively, our observations posit DNA methylation as a safeguard against antitumorigenic cell-fate decisions in PRC2-loss cancer to promote cancer pathogenesis, which can be therapeutically exploited by DNMT1-targeted therapy. SIGNIFICANCE: PRC2 inactivation drives oncogenesis in various cancers, but therapeutically targeting PRC2 loss has remained challenging. Here we show that PRC2-inactivating mutations set up a tumor context-specific liability for therapeutic intervention via DNMT1 inhibitors, which leads to innate immune signaling mediated by sensing of derepressed retrotransposons and accompanied by enhanced cytotoxicity. See related commentary by Guil and Esteller, p. 2020. This article is highlighted in the In This Issue feature, p. 2007.

SWI/SNF Complex Mutations Promote Thyroid Tumor Progression and Insensitivity to Redifferentiation Therapies.

Mutations of subunits of the SWI/SNF chromatin remodeling complexes occur commonly in cancers of different lineages, including advanced thyroid cancers. Here we show that thyroid-specific loss of Arid1a, Arid2, or Smarcb1 in mouse BRAFV600E-mutant tumors promotes disease progression and decreased survival, associated with lesion-specific effects on chromatin accessibility and differentiation. As compared with normal thyrocytes, BRAFV600E-mutant mouse papillary thyroid cancers have decreased lineage transcription factor expression and accessibility to their target DNA binding sites, leading to impairment of thyroid-differentiated gene expression and radioiodine incorporation, which is rescued by MAPK inhibition. Loss of individual SWI/SNF subunits in BRAF tumors leads to a repressive chromatin state that cannot be reversed by MAPK pathway blockade, rendering them insensitive to its redifferentiation effects. Our results show that SWI/SNF complexes are central to the maintenance of differentiated function in thyroid cancers, and their loss confers radioiodine refractoriness and resistance to MAPK inhibitor-based redifferentiation therapies. SIGNIFICANCE: Reprogramming cancer differentiation confers therapeutic benefit in various disease contexts. Oncogenic BRAF silences genes required for radioiodine responsiveness in thyroid cancer. Mutations in SWI/SNF genes result in loss of chromatin accessibility at thyroid lineage specification genes in BRAF-mutant thyroid tumors, rendering them insensitive to the redifferentiation effects of MAPK blockade.This article is highlighted in the In This Issue feature, p. 995.

FOXA1 Mutations Reveal Distinct Chromatin Profiles and Influence Therapeutic Response in Breast Cancer.

Mutations in the pioneer transcription factor FOXA1 are a hallmark of estrogen receptor-positive (ER+) breast cancers. Examining FOXA1 in ∼5,000 breast cancer patients identifies several hotspot mutations in the Wing2 region and a breast cancer-specific mutation SY242CS, located in the third ß strand. Using a clinico-genomically curated cohort, together with breast cancer models, we find that FOXA1 mutations associate with a lower response to aromatase inhibitors. Mechanistically, Wing2 mutations display increased chromatin binding at ER loci upon estrogen stimulation, and an enhanced ER-mediated transcription without changes in chromatin accessibility. In contrast, SY242CS shows neomorphic properties that include the ability to open distinct chromatin regions and activate an alternative cistrome and transcriptome. Structural modeling predicts that SY242CS confers a conformational change that mediates stable binding to a non-canonical DNA motif. Taken together, our results provide insights into how FOXA1 mutations perturb its function to dictate cancer progression and therapeutic response.

gganatogram: An R package for modular visualisation of anatograms and tissues based on ggplot2.

Displaying data onto anatomical structures is a convenient technique to quickly observe tissue related information. However, drawing tissues is a complex task that requires both expertise in anatomy and the arts. While web based applications exist for displaying gene expression on anatograms, other non-genetic disciplines lack similar tools. Moreover, web based tools often lack the modularity associated with packages in programming languages, such as R. Here I present gganatogram, an R package used to plot modular species anatograms based on a combination of the graphical grammar of ggplot2 and the publicly available anatograms from the Expression Atlas. This combination allows for quick and easy, modular, and reproducible generation of anatograms. Using only one command and a data frame with tissue name, group, colour, and  value, this tool enables the user to visualise specific human and mouse tissues with desired colours, grouped by a variable, or displaying a desired value, such as gene-expression, pharmacokinetics, or bacterial load across selected tissues. gganatogram consists of 5 highly annotated organisms, male/female human/mouse, and a cell anatogram. It further consists of 24 other less annotated organisms from the animal and plant kingdom. I hope that this tool will be useful by the wider community in biological sciences. Community members are welcome to submit additional anatograms, which can be incorporated into the package. A stable version gganatogram has been deposited to neuroconductor, and a development version can be found on  github/jespermaag/gganatogram. An interactive shiny app of gganatogram can be found on  https://jespermaag.shinyapps.io/gganatogram/, which allows for non-R users to create anatograms.

Benchmarking of RNA-sequencing analysis workflows using whole-transcriptome RT-qPCR expression data.

RNA-sequencing has become the gold standard for whole-transcriptome gene expression quantification. Multiple algorithms have been developed to derive gene counts from sequencing reads. While a number of benchmarking studies have been conducted, the question remains how individual methods perform at accurately quantifying gene expression levels from RNA-sequencing reads. We performed an independent benchmarking study using RNA-sequencing data from the well established MAQCA and MAQCB reference samples. RNA-sequencing reads were processed using five workflows (Tophat-HTSeq, Tophat-Cufflinks, STAR-HTSeq, Kallisto and Salmon) and resulting gene expression measurements were compared to expression data generated by wet-lab validated qPCR assays for all protein coding genes. All methods showed high gene expression correlations with qPCR data. When comparing gene expression fold changes between MAQCA and MAQCB samples, about 85% of the genes showed consistent results between RNA-sequencing and qPCR data. Of note, each method revealed a small but specific gene set with inconsistent expression measurements. A significant proportion of these method-specific inconsistent genes were reproducibly identified in independent datasets. These genes were typically smaller, had fewer exons, and were lower expressed compared to genes with consistent expression measurements. We propose that careful validation is warranted when evaluating RNA-seq based expression profiles for this specific gene set.

Novel Aberrations Uncovered in Barrett's Esophagus and Esophageal Adenocarcinoma Using Whole Transcriptome Sequencing.

Esophageal adenocarcinoma (EAC) has one of the fastest increases in incidence of any cancer, along with poor five-year survival rates. Barrett's esophagus (BE) is the main risk factor for EAC; however, the mechanisms driving EAC development remain poorly understood. Here, transcriptomic profiling was performed using RNA-sequencing (RNA-seq) on premalignant and malignant Barrett's tissues to better understand this disease. Machine-learning and network analysis methods were applied to discover novel driver genes for EAC development. Identified gene expression signatures for the distinction of EAC from BE were validated in separate datasets. An extensive analysis of the noncoding RNA (ncRNA) landscape was performed to determine the involvement of novel transcriptomic elements in Barrett's disease and EAC. Finally, transcriptomic mutational investigation of genes that are recurrently mutated in EAC was performed. Through these approaches, novel driver genes were discovered for EAC, which involved key cell cycle and DNA repair genes, such as BRCA1 and PRKDC. A novel 4-gene signature (CTSL, COL17A1, KLF4, and E2F3) was identified, externally validated, and shown to provide excellent distinction of EAC from BE. Furthermore, expression changes were observed in 685 long noncoding RNAs (lncRNA) and a systematic dysregulation of repeat elements across different stages of Barrett's disease, with wide-ranging downregulation of Alu elements in EAC. Mutational investigation revealed distinct pathways activated between EAC tissues with or without TP53 mutations compared with Barrett's disease. In summary, transcriptome sequencing revealed altered expression of numerous novel elements, processes, and networks in EAC and premalignant BE.Implications: This study identified opportunities to improve early detection and treatment of patients with BE and esophageal adenocarcinoma. Mol Cancer Res; 15(11); 1558-69. ©2017 AACR.

The long non-coding RNA NEAT1 is responsive to neuronal activity and is associated with hyperexcitability states.

Despite their abundance, the molecular functions of long non-coding RNAs in mammalian nervous systems remain poorly understood. Here we show that the long non-coding RNA, NEAT1, directly modulates neuronal excitability and is associated with pathological seizure states. Specifically, NEAT1 is dynamically regulated by neuronal activity in vitro and in vivo, binds epilepsy-associated potassium channel-interacting proteins including KCNAB2 and KCNIP1, and induces a neuronal hyper-potentiation phenotype in iPSC-derived human cortical neurons following antisense oligonucleotide knockdown. Next generation sequencing reveals a strong association of NEAT1 with increased ion channel gene expression upon activation of iPSC-derived neurons following NEAT1 knockdown. Furthermore, we show that while NEAT1 is acutely down-regulated in response to neuronal activity, repeated stimulation results in NEAT1 becoming chronically unresponsive in independent in vivo rat model systems relevant to temporal lobe epilepsy. We extended previous studies showing increased NEAT1 expression in resected cortical tissue from high spiking regions of patients suffering from intractable seizures. Our results indicate a role for NEAT1 in modulating human neuronal activity and suggest a novel mechanistic link between an activity-dependent long non-coding RNA and epilepsy.

Corrigendum: Dynamic expression of long noncoding RNAs and repeat elements in synaptic plasticity.

[This corrects the article on p. 351 in vol. 9, PMID: 26483626.].

Dynamic expression of long noncoding RNAs and repeat elements in synaptic plasticity.

Long-term potentiation (LTP) of synaptic transmission is recognized as a cellular mechanism for learning and memory storage. Although de novo gene transcription is known to be required in the formation of stable LTP, the molecular mechanisms underlying synaptic plasticity remain elusive. Noncoding RNAs have emerged as major regulatory molecules that are abundantly and specifically expressed in the mammalian brain. By combining RNA-seq analysis with LTP induction in the dentate gyrus of live rats, we provide the first global transcriptomic analysis of synaptic plasticity in the adult brain. Expression profiles of mRNAs and long noncoding RNAs (lncRNAs) were obtained at 30 min, 2 and 5 h after high-frequency stimulation of the perforant pathway. The temporal analysis revealed dynamic expression profiles of lncRNAs with many positively, and highly, correlated to protein-coding genes with known roles in synaptic plasticity, suggesting their possible involvement in LTP. In light of observations suggesting a role for retrotransposons in brain function, we examined the expression of various classes of repeat elements. Our analysis identifies dynamic regulation of LINE1 and SINE retrotransposons, and extensive regulation of tRNA. These experiments reveal a hitherto unknown complexity of gene expression in long-term synaptic plasticity involving the dynamic regulation of lncRNAs and repeat elements. These findings provide a broader foundation for elucidating the transcriptional and epigenetic regulation of synaptic plasticity in both the healthy brain and in neurodegenerative and neuropsychiatric disorders.

Contrasting protective effects of cannabinoids against oxidative stress and amyloid-ß evoked neurotoxicity in vitro.

Cannabinoids have been widely reported to have neuroprotective properties in vitro and in vivo. In this study we compared the effects of CB1 and CB2 receptor-selective ligands, the endocannabinoid anandamide and the phytocannabinoid cannabidiol, against oxidative stress and the toxic hallmark Alzheimer's protein, ß-amyloid (Aß) in neuronal cell lines. PC12 or SH-SY5Y cells were selectively exposed to either hydrogen peroxide, tert-butyl hydroperoxide or Aß, alone or in the presence of the CB1 specific agonist arachidonyl-2'-chloroethylamide (ACEA), CB2 specific agonist JWH-015, anandamide or cannabidiol. Cannabidiol improved cell viability in response to tert-butyl hydroperoxide in PC12 and SH-SY5Y cells, while hydrogen peroxide-mediated toxicity was unaffected by cannabidiol pretreatment. Aß exposure evoked a loss of cell viability in PC12 cells. Of the cannabinoids tested, only anandamide was able to inhibit Aß-evoked neurotoxicity. ACEA had no effect on Aß-evoked neurotoxicity, suggesting a CB1 receptor-independent effect of anandamide. JWH-015 pretreatment was also without protective influence on PC12 cells from either pro-oxidant or Aß exposure. None of the cannabinoids directly inhibited or disrupted preformed Aß fibrils and aggregates. In conclusion, the endocannabinoid anandamide protects neuronal cells from Aß exposure via a pathway unrelated to CB1 or CB2 receptor activation. The protective effect of cannabidiol against oxidative stress does not confer protection against Aß exposure, suggesting divergent pathways for neuroprotection of these two cannabinoids.