High temporal resolution RNA-seq time course data reveals widespread synchronous activation between mammalian lncRNAs and neighboring protein-coding genes.

The advent of massively parallel sequencing revealed extensive transcription beyond protein-coding genes, identifying tens of thousands of long noncoding RNAs (lncRNAs). Selected functional examples raised the possibility that lncRNAs, as a class, may maintain broad regulatory roles. Expression of lncRNAs is strongly linked with adjacent protein-coding gene expression, suggesting potential cis-regulatory functions. A more detailed understanding of these regulatory roles may be obtained through careful examination of the precise timing of lncRNA expression relative to adjacent protein-coding genes. Despite the diversity of reported lncRNA regulatory mechanisms, where causal cis-regulatory relationships exist, lncRNA transcription is expected to precede changes in target gene expression. Using a high temporal resolution RNA-seq time course, we profiled the expression dynamics of several thousand lncRNAs and protein-coding genes in synchronized, transitioning human cells. Our findings reveal that lncRNAs are expressed synchronously with adjacent protein-coding genes. Analysis of lipopolysaccharide-activated mouse dendritic cells revealed the same temporal relationship observed in transitioning human cells. Our findings suggest broad-scale cis-regulatory roles for lncRNAs are not common. The strong association between lncRNAs and adjacent genes may instead indicate an origin as transcriptional by-products from active protein-coding gene promoters and enhancers.

Transcriptomic Profiling of Human Pluripotent Stem Cell-derived Retinal Pigment Epithelium over Time.

Human pluripotent stem cell (hPSC)-derived progenies are immature versions of cells, presenting a potential limitation to the accurate modelling of diseases associated with maturity or age. Hence, it is important to characterise how closely cells used in culture resemble their native counterparts. In order to select appropriate time points of retinal pigment epithelium (RPE) cultures that reflect native counterparts, we characterised the transcriptomic profiles of the hPSC-derived RPE cells from 1- and 12-month cultures. We differentiated the human embryonic stem cell line H9 into RPE cells, performed single-cell RNA-sequencing of a total of 16,576 cells to assess the molecular changes of the RPE cells across these two culture time points. Our results indicate the stability of the RPE transcriptomic signature, with no evidence of an epithelial-mesenchymal transition, and with the maturing populations of the RPE observed with time in culture. Assessment of Gene Ontology pathways revealed that as the cultures age, RPE cells upregulate expression of genes involved in metal binding and antioxidant functions. This might reflect an increased ability to handle oxidative stress as cells mature. Comparison with native human RPE data confirms a maturing transcriptional profile of RPE cells in culture. These results suggest that long-term in vitro culture of RPE cells allows the modelling of specific phenotypes observed in native mature tissues. Our work highlights the transcriptional landscape of hPSC-derived RPE cells as they age in culture, which provides a reference for native and patient samples to be benchmarked against.

Repeat RNAs associate with replication forks and post-replicative DNA.

Noncoding RNA has a proven ability to direct and regulate chromatin modifications by acting as scaffolds between DNA and histone-modifying complexes. However, it is unknown if ncRNA plays any role in DNA replication and epigenome maintenance, including histone eviction and reinstallment of histone modifications after genome duplication. Isolation of nascent chromatin has identified a large number of RNA-binding proteins in addition to unknown components of the replication and epigenetic maintenance machinery. Here, we isolated and characterized long and short RNAs associated with nascent chromatin at active replication forks and track RNA composition during chromatin maturation across the cell cycle. Shortly after fork passage, GA-rich-, alpha- and TElomeric Repeat-containing RNAs (TERRA) are associated with replicated DNA. These repeat containing RNAs arise from loci undergoing replication, suggesting an interaction in cis. Post-replication during chromatin maturation, and even after mitosis in G1, the repeats remain enriched on DNA. This suggests that specific types of repeat RNAs are transcribed shortly after DNA replication and stably associate with their loci of origin throughout the cell cycle. The presented method and data enable studies of RNA interactions with replication forks and post-replicative chromatin and provide insights into how repeat RNAs and their engagement with chromatin are regulated with respect to DNA replication and across the cell cycle.

Adar3 Is Involved in Learning and Memory in Mice.

The amount of regulatory RNA encoded in the genome and the extent of RNA editing by the post-transcriptional deamination of adenosine to inosine (A-I) have increased with developmental complexity and may be an important factor in the cognitive evolution of animals. The newest member of the A-I editing family of ADAR proteins, the vertebrate-specific ADAR3, is highly expressed in the brain, but its functional significance is unknown. In vitro studies have suggested that ADAR3 acts as a negative regulator of A-I RNA editing but the scope and underlying mechanisms are also unknown. Meta-analysis of published data indicates that mouse Adar3 expression is highest in the hippocampus, thalamus, amygdala, and olfactory region. Consistent with this, we show that mice lacking exon 3 of Adar3 (which encodes two double stranded RNA binding domains) have increased levels of anxiety and deficits in hippocampus-dependent short- and long-term memory formation. RNA sequencing revealed a dysregulation of genes involved in synaptic function in the hippocampi of Adar3-deficient mice. We also show that ADAR3 transiently translocates from the cytoplasm to the nucleus upon KCl-mediated activation in SH-SY5Y cells. These results indicate that ADAR3 contributes to cognitive processes in mammals.

THC exposure of human iPSC neurons impacts genes associated with neuropsychiatric disorders.

There is a strong association between cannabis use and schizophrenia but the underlying cellular links are poorly understood. Neurons derived from human-induced pluripotent stem cells (hiPSCs) offer a platform for investigating both baseline and dynamic changes in human neural cells. Here, we exposed neurons derived from hiPSCs to Δ9-tetrahydrocannabinol (THC), and identified diagnosis-specific differences not detectable in vehicle-controls. RNA transcriptomic analyses revealed that THC administration, either by acute or chronic exposure, dampened the neuronal transcriptional response following potassium chloride (KCl)-induced neuronal depolarization. THC-treated neurons displayed significant synaptic, mitochondrial, and glutamate signaling alterations that may underlie their failure to activate appropriately; this blunted response resembles effects previously observed in schizophrenia hiPSC- derived neurons. Furthermore, we show a significant alteration in THC-related genes associated with autism and intellectual disability, suggesting shared molecular pathways perturbed in neuropsychiatric disorders that are exacerbated by THC.

Evidence that TLR4 Is Not a Receptor for Saturated Fatty Acids but Mediates Lipid-Induced Inflammation by Reprogramming Macrophage Metabolism.

Chronic inflammation is a hallmark of obesity and is linked to the development of numerous diseases. The activation of toll-like receptor 4 (TLR4) by long-chain saturated fatty acids (lcSFAs) is an important process in understanding how obesity initiates inflammation. While experimental evidence supports an important role for TLR4 in obesity-induced inflammation in vivo, via a mechanism thought to involve direct binding to and activation of TLR4 by lcSFAs, several lines of evidence argue against lcSFAs being direct TLR4 agonists. Using multiple orthogonal approaches, we herein provide evidence that while loss-of-function models confirm that TLR4 does, indeed, regulate lcSFA-induced inflammation, TLR4 is not a receptor for lcSFAs. Rather, we show that TLR4-dependent priming alters cellular metabolism, gene expression, lipid metabolic pathways, and membrane lipid composition, changes that are necessary for lcSFA-induced inflammation. These results reconcile previous discordant observations and challenge the prevailing view of TLR4's role in initiating obesity-induced inflammation.

High resolution temporal transcriptomics of mouse embryoid body development reveals complex expression dynamics of coding and noncoding loci.

Cellular responses to stimuli are rapid and continuous and yet the vast majority of investigations of transcriptional responses during developmental transitions typically use long interval time courses; limiting the available interpretive power. Moreover, such experiments typically focus on protein-coding transcripts, ignoring the important impact of long noncoding RNAs. We therefore evaluated coding and noncoding expression dynamics at unprecedented temporal resolution (6-hourly) in differentiating mouse embryonic stem cells and report new insight into molecular processes and genome organization. We present a highly resolved differentiation cascade that exhibits coding and noncoding transcriptional alterations, transcription factor network interactions and alternative splicing events, little of which can be resolved by long-interval developmental time-courses. We describe novel short lived and cycling patterns of gene expression and dissect temporally ordered gene expression changes in response to transcription factors. We elucidate patterns in gene co-expression across the genome, describe asynchronous transcription at bidirectional promoters and functionally annotate known and novel regulatory lncRNAs. These findings highlight the complex and dynamic molecular events underlying mammalian differentiation that can only be observed though a temporally resolved time course.

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.

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.

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.

Activity-Dependent Changes in Gene Expression in Schizophrenia Human-Induced Pluripotent Stem Cell Neurons.

IMPORTANCE: Schizophrenia candidate genes participate in common molecular pathways that are regulated by activity-dependent changes in neurons. One important next step is to further our understanding on the role of activity-dependent changes of gene expression in the etiopathogenesis of schizophrenia. OBJECTIVE: To examine whether neuronal activity-dependent changes of gene expression are dysregulated in schizophrenia. DESIGN, SETTING, AND PARTICIPANTS: Neurons differentiated from human-induced pluripotent stem cells derived from 4 individuals with schizophrenia and 4 unaffected control individuals were depolarized using potassium chloride. RNA was extracted followed by genome-wide profiling of the transcriptome. Neurons were planted on June 21, 2013, and harvested on August 2, 2013. MAIN OUTCOMES AND MEASURES: We performed differential expression analysis and gene coexpression analysis to identify activity-dependent or disease-specific changes of the transcriptome. Gene expression differences were assessed with linear models. Furthermore, we used gene set analyses to identify coexpressed modules that are enriched for schizophrenia risk genes. RESULTS: We identified 1669 genes that were significantly different in schizophrenia-associated vs control human-induced pluripotent stem cell-derived neurons and 1199 genes that are altered in these cells in response to depolarization (linear models at false discovery rate ≤0.05). The effect of activity-dependent changes of gene expression in schizophrenia-associated neurons (59 significant genes at false discovery rate ≤0.05) was attenuated compared with control samples (594 significant genes at false discovery rate ≤0.05). Using gene coexpression analysis, we identified 2 modules (turquoise and brown) that were associated with diagnosis status and 2 modules (yellow and green) that were associated with depolarization at a false discovery rate of ≤0.05. For 3 of the 4 modules, we found enrichment with schizophrenia-associated variants: brown (χ2 = 20.68; P = .002), turquoise (χ2 = 12.95; P = .04), and yellow (χ2 = 15.34; P = .02). CONCLUSIONS AND RELEVANCE: In this analysis, candidate genes clustered within gene networks that were associated with a blunted effect of activity-dependent changes of gene expression in schizophrenia-associated neurons. Overall, these findings link schizophrenia candidate genes with specific molecular functions in neurons, which could be used to examine underlying mechanisms and therapeutic interventions related to schizophrenia.

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.

Sphingosine kinase 1 isoform-specific interactions in breast cancer.

Sphingosine kinase 1 (SK1) is a signaling enzyme that catalyzes the formation of sphingosine-1-phosphate. Overexpression of SK1 is causally associated with breast cancer progression and resistance to therapy. SK1 inhibitors are currently being investigated as promising breast cancer therapies. Two major transcriptional isoforms, SK143 kDa and SK151 kDa, have been identified; however, the 51 kDa variant is predominant in breast cancer cells. No studies have investigated the protein-protein interactions of the 51 kDa isoform and whether the two SK1 isoforms differ significantly in their interactions. Seeking an understanding of the regulation and role of SK1, we used a triple-labeling stable isotope labeling by amino acids in cell culture-based approach to identify SK1-interacting proteins common and unique to both isoforms. Of approximately 850 quantified proteins in SK1 immunoprecipitates, a high-confidence list of 30 protein interactions with each SK1 isoform was generated via a meta-analysis of multiple experimental replicates. Many of the novel identified SK1 interaction partners such as supervillin, drebrin, and the myristoylated alanine-rich C-kinase substrate-related protein supported and highlighted previously implicated roles of SK1 in breast cancer cell migration, adhesion, and cytoskeletal remodeling. Of these interactions, several were found to be exclusive to the 43 kDa isoform of SK1, including the protein phosphatase 2A, a previously identified SK1-interacting protein. Other proteins such as allograft inflammatory factor 1-like protein, the latent-transforming growth factor ß-binding protein, and dipeptidyl peptidase 2 were found to associate exclusively with the 51 kDa isoform of SK1. In this report, we have identified common and isoform-specific SK1-interacting partners that provide insight into the molecular mechanisms that drive SK1-mediated oncogenicity.