Firing Rate Homeostasis Can Occur in the Absence of Neuronal Activity-Regulated Transcription.

Despite dynamic inputs, neuronal circuits maintain relatively stable firing rates over long periods. This maintenance of firing rate, or firing rate homeostasis, is likely mediated by homeostatic mechanisms such as synaptic scaling and regulation of intrinsic excitability. Because some of these homeostatic mechanisms depend on transcription of activity-regulated genes, including Arc and Homer1a, we hypothesized that activity-regulated transcription would be required for firing rate homeostasis. Surprisingly, however, we found that cultured mouse cortical neurons from both sexes grown on multi-electrode arrays homeostatically adapt their firing rates to persistent pharmacological stimulation even when activity-regulated transcription is disrupted. Specifically, we observed firing rate homeostasis in Arc knock-out neurons, as well as knock-out neurons lacking the activity-regulated transcription factors AP1 and SRF. Firing rate homeostasis also occurred normally during acute pharmacological blockade of transcription. Thus, firing rate homeostasis in response to increased neuronal activity can occur in the absence of neuronal-activity-regulated transcription.SIGNIFICANCE STATEMENT Neuronal circuits maintain relatively stable firing rates even in the face of dynamic circuit inputs. Understanding the molecular mechanisms that enable this firing rate homeostasis could potentially provide insight into neuronal diseases that present with an imbalance of excitation and inhibition. It has long been proposed that activity-regulated transcription could underlie firing rate homeostasis because activity-regulated genes turn on when neurons are above their target firing rates and include many genes that could regulate firing rate. Surprisingly, despite this prediction, we found that cortical neurons can undergo firing rate homeostasis in the absence of activity-regulated transcription, indicating that firing rate homeostasis can be controlled by non-transcriptional mechanisms.

Blue Light Increases Neuronal Activity-Regulated Gene Expression in the Absence of Optogenetic Proteins.

Optogenetics is widely used to control diverse cellular functions with light, requiring experimenters to expose cells to bright light. Because extended exposure to visible light can be toxic to cells, it is important to characterize the effects of light stimulation on cellular function in the absence of optogenetic proteins. Here we exposed mouse cortical cultures with no exogenous optogenetic proteins to several hours of flashing blue, red, or green light. We found that exposing these cultures to as short as 1 h of blue light, but not red or green light, results in an increase in the expression of neuronal activity-regulated genes. Our findings suggest that blue light stimulation is ill suited to long-term optogenetic experiments, especially those that measure transcription, and they emphasize the importance of performing light-only control experiments in samples without optogenetic proteins.

The neuronal stimulation-transcription coupling map.

Neurons transcribe different genes in response to different extracellular stimuli, and these genes regulate neuronal plasticity. Thus, understanding how different stimuli regulate different stimulus-dependent gene modules would deepen our understanding of plasticity. To systematically dissect the coupling between stimulation and transcription, we propose creating a 'stimulation-transcription coupling map' that describes the transcription response to each possible extracellular stimulus. While we are currently far from having a complete map, recent genomic experiments have begun to facilitate its creation. Here, we describe the current state of the stimulation-transcription coupling map as well as the transcriptional regulation that enables this coupling.

Cell-type-specific programs for activity-regulated gene expression.

Experience leaves a lasting mark on neural circuit function in part through activity-regulated gene (ARG) expression. New genome wide approaches have revealed that ARG programs are highly cell-type-specific, raising the possibility that they mediate different forms of experience-dependent plasticity in different cell types. The cell-type specificity of these gene programs is achieved by a combination of cell-intrinsic mechanisms that determine the transcriptional response of each neuronal subtype to a given stimulus and by cell-extrinsic mechanisms that influence the nature of the stimulus a cell receives. A better understanding of these mechanisms could usher in an era of molecular systems neuroscience in which genetic perturbations of cell-type-specific plasticities are assessed using electrophysiology and in vivo imaging to reveal the neural basis of adaptive behaviors.

Different Neuronal Activity Patterns Induce Different Gene Expression Programs.

A vast number of different neuronal activity patterns could each induce a different set of activity-regulated genes. Mapping this coupling between activity pattern and gene induction would allow inference of a neuron's activity-pattern history from its gene expression and improve our understanding of activity-pattern-dependent synaptic plasticity. In genome-scale experiments comparing brief and sustained activity patterns, we reveal that activity-duration history can be inferred from gene expression profiles. Brief activity selectively induces a small subset of the activity-regulated gene program that corresponds to the first of three temporal waves of genes induced by sustained activity. Induction of these first-wave genes is mechanistically distinct from that of the later waves because it requires MAPK/ERK signaling but does not require de novo translation. Thus, the same mechanisms that establish the multi-wave temporal structure of gene induction also enable different gene sets to be induced by different activity durations.

Total RNA-seq to identify pharmacological effects on specific stages of mRNA synthesis.

Pharmacological perturbation is a powerful tool for understanding mRNA synthesis, but identification of the specific steps of this multi-step process that are targeted by small molecules remains challenging. Here we applied strand-specific total RNA sequencing (RNA-seq) to identify and distinguish specific pharmacological effects on transcription and pre-mRNA processing in human cells. We found unexpectedly that the natural product isoginkgetin, previously described as a splicing inhibitor, inhibits transcription elongation. Compared to well-characterized elongation inhibitors that target CDK9, isoginkgetin caused RNA polymerase accumulation within a broader promoter-proximal band, indicating that elongation inhibition by isoginkgetin occurs after release from promoter-proximal pause. RNA-seq distinguished isoginkgetin and CDK9 inhibitors from topoisomerase I inhibition, which alters elongation across gene bodies. We were able to detect these and other specific defects in mRNA synthesis at low sequencing depth using simple metagene-based metrics. These metrics now enable total-RNA-seq-based screening for high-throughput identification of pharmacological effects on individual stages of mRNA synthesis.

Brain-wide maps of Fos expression during fear learning and recall.

Fos induction during learning labels neuronal ensembles in the hippocampus that encode a specific physical environment, revealing a memory trace. In the cortex and other regions, the extent to which Fos induction during learning reveals specific sensory representations is unknown. Here we generate high-quality brain-wide maps of Fos mRNA expression during auditory fear conditioning and recall in the setting of the home cage. These maps reveal a brain-wide pattern of Fos induction that is remarkably similar among fear conditioning, shock-only, tone-only, and fear recall conditions, casting doubt on the idea that Fos reveals auditory-specific sensory representations. Indeed, novel auditory tones lead to as much gene induction in visual as in auditory cortex, while familiar (nonconditioned) tones do not appreciably induce Fos anywhere in the brain. Fos expression levels do not correlate with physical activity, suggesting that they are not determined by behavioral activity-driven alterations in sensory experience. In the thalamus, Fos is induced more prominently in limbic than in sensory relay nuclei, suggesting that Fos may be most sensitive to emotional state. Thus, our data suggest that Fos expression during simple associative learning labels ensembles activated generally by arousal rather than specifically by a particular sensory cue.

MPRAnator: a web-based tool for the design of massively parallel reporter assay experiments.

MOTIVATION: With the rapid advances in DNA synthesis and sequencing technologies and the continuing decline in the associated costs, high-throughput experiments can be performed to investigate the regulatory role of thousands of oligonucleotide sequences simultaneously. Nevertheless, designing high-throughput reporter assay experiments such as massively parallel reporter assays (MPRAs) and similar methods remains challenging. RESULTS: We introduce MPRAnator, a set of tools that facilitate rapid design of MPRA experiments. With MPRA Motif design, a set of variables provides fine control of how motifs are placed into sequences, thereby allowing the investigation of the rules that govern transcription factor (TF) occupancy. MPRA single-nucleotide polymorphism design can be used to systematically examine the functional effects of single or combinations of single-nucleotide polymorphisms at regulatory sequences. Finally, the Transmutation tool allows for the design of negative controls by permitting scrambling, reversing, complementing or introducing multiple random mutations in the input sequences or motifs. AVAILABILITY AND IMPLEMENTATION: MPRAnator tool set is implemented in Python, Perl and Javascript and is freely available at www.genomegeek.com and www.sanger.ac.uk/science/tools/mpranator The source code is available on www.github.com/hemberg-lab/MPRAnator/ under the MIT license. The REST API allows programmatic access to MPRAnator using simple URLs. CONTACT: igs@sanger.ac.uk or mh26@sanger.ac.ukSupplementary information: Supplementary data are available at Bioinformatics online.

RNA sequencing from neural ensembles activated during fear conditioning in the mouse temporal association cortex.

The stable formation of remote fear memories is thought to require neuronal gene induction in cortical ensembles that are activated during learning. However, the set of genes expressed specifically in these activated ensembles is not known; knowledge of such transcriptional profiles may offer insights into the molecular program underlying stable memory formation. Here we use RNA-Seq to identify genes whose expression is enriched in activated cortical ensembles labeled during associative fear learning. We first establish that mouse temporal association cortex (TeA) is required for remote recall of auditory fear memories. We then perform RNA-Seq in TeA neurons that are labeled by the activity reporter Arc-dVenus during learning. We identify 944 genes with enriched expression in Arc-dVenus+ neurons. These genes include markers of L2/3, L5b, and L6 excitatory neurons but not glial or inhibitory markers, confirming Arc-dVenus to be an excitatory neuron-specific but non-layer-specific activity reporter. Cross comparisons to other transcriptional profiles show that 125 of the enriched genes are also activity-regulated in vitro or induced by visual stimulus in the visual cortex, suggesting that they may be induced generally in the cortex in an experience-dependent fashion. Prominent among the enriched genes are those encoding potassium channels that down-regulate neuronal activity, suggesting the possibility that part of the molecular program induced by fear conditioning may initiate homeostatic plasticity.

High-throughput functional comparison of promoter and enhancer activities.

Promoters initiate RNA synthesis, and enhancers stimulate promoter activity. Whether promoter and enhancer activities are encoded distinctly in DNA sequences is unknown. We measured the enhancer and promoter activities of thousands of DNA fragments transduced into mouse neurons. We focused on genomic loci bound by the neuronal activity-regulated coactivator CREBBP, and we measured enhancer and promoter activities both before and after neuronal activation. We find that the same sequences typically encode both enhancer and promoter activities. However, gene promoters generate more promoter activity than distal enhancers, despite generating similar enhancer activity. Surprisingly, the greater promoter activity of gene promoters is not due to conventional core promoter elements or splicing signals. Instead, we find that particular transcription factor binding motifs are intrinsically biased toward the generation of promoter activity, whereas others are not. Although the specific biases we observe may be dependent on experimental or cellular context, our results suggest that gene promoters are distinguished from distal enhancers by specific complements of transcriptional activators.

Genomic Views of Transcriptional Enhancers: Essential Determinants of Cellular Identity and Activity-Dependent Responses in the CNS.

Sprinkled throughout the genome are a million regulatory sequences called transcriptional enhancers that activate gene promoters in the right cells, at the right time. Enhancers endow the brain with its incredible diversity of cell types and also translate neural activity into gene induction. Thanks to rapid advances in genomic technologies, it is now possible to identify thousands of enhancers rapidly, test their transcriptional function en masse, and address their neurobiological functions via genome editing. Enhancers also promise to be a great technological opportunity for neuroscience, offering the potential for cell-type-specific genetic labeling and manipulation without the need for transgenesis. The objective of this review and the accompanying 2015 SfN mini-symposium is to highlight the use of new and emerging genomic technologies to probe enhancer function in the nervous system. SIGNIFICANCE STATEMENT: Transcriptional enhancers turn on genes in the right cells, at the right time. Enhancers are also the genomic sequences that encode the incredible diversity of cell types in the brain and enable the brain to turn genes on in response to new experiences. New technology enables enhancers to be found and manipulated. The study of enhancers promises to inform our understanding of brain development and function. The application of enhancer technology holds promise in accelerating basic neuroscience research and enabling gene therapies to be targeted to specific cell types in the brain.

Recent advances in functional assays of transcriptional enhancers.

In this special edition of Genomics, we present reviews of the current state of the field in identifying and functionally understanding transcriptional enhancers in cells and developing tissues. Typically several enhancers coordinate the expression of an individual target gene, each controlling that gene's expression in specific cell types at specific times. Until recently, identifying each gene's enhancers had been challenging because enhancers do not occupy prescribed locations relative to their target genes. Recently there have been powerful advances in DNA sequencing and other technologies that make it possible to identify the majority of enhancers in virtually any cell type of interest. The reviews in this edition of Genomics highlight some of these new and powerful approaches.

A circuit mechanism for differentiating positive and negative associations.

The ability to differentiate stimuli predicting positive or negative outcomes is critical for survival, and perturbations of emotional processing underlie many psychiatric disease states. Synaptic plasticity in the basolateral amygdala complex (BLA) mediates the acquisition of associative memories, both positive and negative. Different populations of BLA neurons may encode fearful or rewarding associations, but the identifying features of these populations and the synaptic mechanisms of differentiating positive and negative emotional valence have remained unknown. Here we show that BLA neurons projecting to the nucleus accumbens (NAc projectors) or the centromedial amygdala (CeM projectors) undergo opposing synaptic changes following fear or reward conditioning. We find that photostimulation of NAc projectors supports positive reinforcement while photostimulation of CeM projectors mediates negative reinforcement. Photoinhibition of CeM projectors impairs fear conditioning and enhances reward conditioning. We characterize these functionally distinct neuronal populations by comparing their electrophysiological, morphological and genetic features. Overall, we provide a mechanistic explanation for the representation of positive and negative associations within the amygdala.

Enhancer RNAs: a class of long noncoding RNAs synthesized at enhancers.

Recent studies have revealed that active enhancers are transcribed, producing a class of noncoding RNAs called enhancer RNAs (eRNAs). eRNAs are distinct from long noncoding RNAs (lncRNAs), but these two species of noncoding RNAs may share a similar role in the activation of mRNA transcription. Emerging studies, showing that eRNAs function in controlling mRNA transcription, challenge the idea that enhancers are merely sites of transcription factor assembly. Instead, communication between promoters and enhancers can be bidirectional with promoters required to activate enhancer transcription. Reciprocally, eRNAs may then facilitate enhancer-promoter interaction or activate promoter-driven transcription.

SnapShot-Seq: a method for extracting genome-wide, in vivo mRNA dynamics from a single total RNA sample.

mRNA synthesis, processing, and destruction involve a complex series of molecular steps that are incompletely understood. Because the RNA intermediates in each of these steps have finite lifetimes, extensive mechanistic and dynamical information is encoded in total cellular RNA. Here we report the development of SnapShot-Seq, a set of computational methods that allow the determination of in vivo rates of pre-mRNA synthesis, splicing, intron degradation, and mRNA decay from a single RNA-Seq snapshot of total cellular RNA. SnapShot-Seq can detect in vivo changes in the rates of specific steps of splicing, and it provides genome-wide estimates of pre-mRNA synthesis rates comparable to those obtained via labeling of newly synthesized RNA. We used SnapShot-Seq to investigate the origins of the intrinsic bimodality of metazoan gene expression levels, and our results suggest that this bimodality is partly due to spillover of transcriptional activation from highly expressed genes to their poorly expressed neighbors. SnapShot-Seq dramatically expands the information obtainable from a standard RNA-Seq experiment.

Integrated genome analysis suggests that most conserved non-coding sequences are regulatory factor binding sites.

More than 98% of a typical vertebrate genome does not code for proteins. Although non-coding regions are sprinkled with short (<200 bp) islands of evolutionarily conserved sequences, the function of most of these unannotated conserved islands remains unknown. One possibility is that unannotated conserved islands could encode non-coding RNAs (ncRNAs); alternatively, unannotated conserved islands could serve as promoter-distal regulatory factor binding sites (RFBSs) like enhancers. Here we assess these possibilities by comparing unannotated conserved islands in the human and mouse genomes to transcribed regions and to RFBSs, relying on a detailed case study of one human and one mouse cell type. We define transcribed regions by applying a novel transcript-calling algorithm to RNA-Seq data obtained from total cellular RNA, and we define RFBSs using ChIP-Seq and DNAse-hypersensitivity assays. We find that unannotated conserved islands are four times more likely to coincide with RFBSs than with unannotated ncRNAs. Thousands of conserved RFBSs can be categorized as insulators based on the presence of CTCF or as enhancers based on the presence of p300/CBP and H3K4me1. While many unannotated conserved RFBSs are transcriptionally active to some extent, the transcripts produced tend to be unspliced, non-polyadenylated and expressed at levels 10 to 100-fold lower than annotated coding or ncRNAs. Extending these findings across multiple cell types and tissues, we propose that most conserved non-coding genomic DNA in vertebrate genomes corresponds to promoter-distal regulatory elements.

Expression divergence measured by transcriptome sequencing of four yeast species.

BACKGROUND: The evolution of gene expression is a challenging problem in evolutionary biology, for which accurate, well-calibrated measurements and methods are crucial. RESULTS: We quantified gene expression with whole-transcriptome sequencing in four diploid, prototrophic strains of Saccharomyces species grown under the same condition to investigate the evolution of gene expression. We found that variation in expression is gene-dependent with large variations in each gene's expression between replicates of the same species. This confounds the identification of genes differentially expressed across species. To address this, we developed a statistical approach to establish significance bounds for inter-species differential expression in RNA-Seq data based on the variance measured across biological replicates. This metric estimates the combined effects of technical and environmental variance, as well as Poisson sampling noise by isolating each component. Despite a paucity of large expression changes, we found a strong correlation between the variance of gene expression change and species divergence (R² = 0.90). CONCLUSION: We provide an improved methodology for measuring gene expression changes in evolutionary diverged species using RNA Seq, where experimental artifacts can mimic evolutionary effects.GEO Accession Number: GSE32679.

Widespread transcription at neuronal activity-regulated enhancers.

We used genome-wide sequencing methods to study stimulus-dependent enhancer function in mouse cortical neurons. We identified approximately 12,000 neuronal activity-regulated enhancers that are bound by the general transcriptional co-activator CBP in an activity-dependent manner. A function of CBP at enhancers may be to recruit RNA polymerase II (RNAPII), as we also observed activity-regulated RNAPII binding to thousands of enhancers. Notably, RNAPII at enhancers transcribes bi-directionally a novel class of enhancer RNAs (eRNAs) within enhancer domains defined by the presence of histone H3 monomethylated at lysine 4. The level of eRNA expression at neuronal enhancers positively correlates with the level of messenger RNA synthesis at nearby genes, suggesting that eRNA synthesis occurs specifically at enhancers that are actively engaged in promoting mRNA synthesis. These findings reveal that a widespread mechanism of enhancer activation involves RNAPII binding and eRNA synthesis.

Neurons detect increases and decreases in oxygen levels using distinct guanylate cyclases.

Homeostatic sensory systems detect small deviations in temperature, water balance, pH, and energy needs to regulate adaptive behavior and physiology. In C. elegans, a homeostatic preference for intermediate oxygen (O2) levels requires cGMP signaling through soluble guanylate cyclases (sGCs), proteins that bind gases through an associated heme group. Here we use behavioral analysis, functional imaging, and genetics to show that reciprocal changes in O2 levels are encoded by sensory neurons that express alternative sets of sGCs. URX sensory neurons are activated by increases in O2 levels, and require the sGCs gcy-35 and gcy-36. BAG sensory neurons are activated by decreases in O2 levels, and require the sGCs gcy-31 and gcy-33. The sGCs are instructive O2 sensors, as forced expression of URX sGC genes causes BAG neurons to detect O2 increases. Both sGC expression and cell-intrinsic dynamics contribute to the differential roles of URX and BAG in O2-dependent behaviors.