Cellular architecture of evolving neuroinflammatory lesions and multiple sclerosis pathology

Multiple sclerosis (MS) is a neurological disease characterized by multifocal lesions and smoldering pathology. Although single-cell analyses provided insights into cytopathology, evolving cellular processes underlying MS remain poorly understood. We investigated the cellular dynamics of MS by modeling temporal and regional rates of disease progression in mouse experimental autoimmune encephalomyelitis (EAE). By performing single-cell spatial expression profiling using in situ sequencing (ISS), we annotated disease neighborhoods and found centrifugal evolution of active lesions. We demonstrated that disease-associated (DA)-glia arise independently of lesions and are dynamically induced and resolved over the disease course. Single-cell spatial mapping of human archival MS spinal cords confirmed the differential distribution of homeostatic and DA-glia, enabled deconvolution of active and inactive lesions into sub-compartments, and identified new lesion areas. By establishing a spatial resource of mouse and human MS neuropathology at a single-cell resolution, our study unveils the intricate cellular dynamics underlying MS.

Spatial epigenome-transcriptome co-profiling of mammalian tissues

Emerging spatial technologies, including spatial transcriptomics and spatial epigenomics, are becoming powerful tools for profiling of cellular states in the tissue context1-5. However, current methods capture only one layer of omics information at a time, precluding the possibility of examining the mechanistic relationship across the central dogma of molecular biology. Here, we present two technologies for spatially resolved, genome-wide, joint profiling of the epigenome and transcriptome by cosequencing chromatin accessibility and gene expression, or histone modifications (H3K27me3, H3K27ac or H3K4me3) and gene expression on the same tissue section at near-single-cell resolution. These were applied to embryonic and juvenile mouse brain, as well as adult human brain, to map how epigenetic mechanisms control transcriptional phenotype and cell dynamics in tissue. Although highly concordant tissue features were identified by either spatial epigenome or spatial transcriptome we also observed distinct patterns, suggesting their differential roles in defining cell states. Linking epigenome to transcriptome pixel by pixel allows the uncovering of new insights in spatial epigenetic priming, differentiation and gene regulation within the tissue architecture. These technologies are of great interest in life science and biomedical research.

Epigenomic priming of immune genes implicates oligodendroglia in multiple sclerosis susceptibility

Multiple sclerosis (MS) is characterized by a targeted attack on oligodendroglia (OLG) and myelin by immune cells, which are thought to be the main drivers of MS susceptibility. We found that immune genes exhibit a primed chromatin state in single mouse and human OLG in a non-disease context, compatible with transitions to immune-competent states in MS. We identified BACH1 and STAT1 as transcription factors involved in immune gene regulation in oligodendrocyte precursor cells (OPCs). A subset of immune genes presents bivalency of H3K4me3/H3K27me3 in OPCs, with Polycomb inhibition leading to their increased activation upon interferon gamma (IFN-γ) treatment. Some MS susceptibility single-nucleotide polymorphisms (SNPs) overlap with these regulatory regions in mouse and human OLG. Treatment of mouse OPCs with IFN-γ leads to chromatin architecture remodeling at these loci and altered expression of interacting genes. Thus, the susceptibility for MS may involve OLG, which therefore constitutes novel targets for immunological-based therapies for MS.

Splicing-associated chromatin signatures: a combinatorial and position-dependent role for histone marks in splicing definition

Alternative splicing relies on the combinatorial recruitment of splicing regulators to specific RNA binding sites. Chromatin has been shown to impact this recruitment. However, a limited number of histone marks have been studied at a global level. In this work, a machine learning approach, applied to extensive epigenomics datasets in human H1 embryonic stem cells and IMR90 foetal fibroblasts, has identified eleven chromatin modifications that differentially mark alternatively spliced exons depending on the level of exon inclusion. These marks act in a combinatorial and position-dependent way, creating characteristic splicing-associated chromatin signatures (SACS). In support of a functional role for SACS in coordinating splicing regulation, changes in the alternative splicing of SACS-marked exons between ten different cell lines correlate with changes in SACS enrichment levels and recruitment of the splicing regulators predicted by RNA motif search analysis. We propose the dynamic nature of chromatin modifications as a mechanism to rapidly fine-tune alternative splicing when necessary.

Chromatin is known to regulate splicing by modulating recruitment of splicing factors. Using machine learning approaches, the authors have underlined a chromatin code for alternative splicing regulation that is conserved amongst cell lines.

Oligodendrocyte Intrinsic miR-27a Controls Myelination and Remyelination

Remyelination requires the generation of new oligodendrocytes (OLs), which are derived from oligodendrocyte progenitor cells (OPCs). Maturation of OPCs into OLs is a multi-step process. Here, we describe a microRNA expressed by OLs, miR-27a, as a regulator of OL development and survival. Increased levels of miR-27a were found in OPCs associated with multiple sclerosis (MS) lesions and in animal models of demyelination. Increased levels of miR-27a led to inhibition of OPC proliferation by cell-cycle arrest, as well as impaired differentiation of human OPCs (hOPCs) and myelination by dysregulating the Wnt-β-catenin signaling pathway. In vivo administration of miR-27a led to suppression of myelinogenic signals, leading to loss of endogenous myelination and remyelination. Our findings provide evidence supporting a critical role for a steady-state level of OL-specific miR-27a in supporting multiple steps in the complex process of OPC maturation and remyelination.

Nuclear role for human Argonaute-1 as an estrogen-dependent transcription coactivator

In mammals, argonaute (AGO) proteins have been characterized for their roles in small RNA-mediated posttranscriptional and also in transcriptional gene silencing. Here, we report a different role for AGO1 in estradiol-triggered transcriptional activation in human cells. We show that in MCF-7 mammary gland cells, AGO1 associates with transcriptional enhancers of estrogen receptor α (ERα) and that this association is up-regulated by treating the cells with estrogen (E2), displaying a positive correlation with the activation of these enhancers. Moreover, we show that AGO1 interacts with ERα and that this interaction is also increased by E2 treatment, but occurs in the absence of RNA. We show that AGO1 acts positively as a coactivator in estradiol-triggered transcription regulation by promoting ERα binding to its enhancers. Consistently, AGO1 depletion decreases long-range contacts between ERα enhancers and their target promoters. Our results point to a role of AGO1 in transcriptional regulation in human cells that is independent from small RNA binding.

An Atlas of Vagal Sensory Neurons and Their Molecular Specialization

Sensory functions of the vagus nerve are critical for conscious perceptions and for monitoring visceral functions in the cardio-pulmonary and gastrointestinal systems. Here, we present a comprehensive identification, classification, and validation of the neuron types in the neural crest (jugular) and placode (nodose) derived vagal ganglia by single-cell RNA sequencing (scRNA-seq) transcriptomic analysis. Our results reveal major differences between neurons derived from different embryonic origins. Jugular neurons exhibit fundamental similarities to the somatosensory spinal neurons, including major types, such as C-low threshold mechanoreceptors (C-LTMRs), A-LTMRs, Aδ-nociceptors, and cold-, and mechano-heat C-nociceptors. In contrast, the nodose ganglion contains 18 distinct types dedicated to surveying the physiological state of the internal body. Our results reveal a vast diversity of vagal neuron types, including many previously unanticipated types, as well as proposed types that are consistent with chemoreceptors, nutrient detectors, baroreceptors, and stretch and volume mechanoreceptors of the respiratory, gastrointestinal, and cardiovascular systems.

•

A comprehensive molecular identification of neuronal types in vagal ganglion complex

•

Prdm12⁺ jugular ganglion neurons share features with spinal somatosensory neurons

•

Phox2b⁺ viscerosensory nodose neurons are molecularly versatile and highly specialized

•

Nodose neuron types are consistent with chemo-, baro-, stretch-, tension-, and volume-sensors

PAD2-Mediated Citrullination Contributes to Efficient Oligodendrocyte Differentiation and Myelination

Citrullination, the deimination of peptidylarginine residues into peptidylcitrulline, has been implicated in the etiology of several diseases. In multiple sclerosis, citrullination is thought to be a major driver of pathology through hypercitrullination and destabilization of myelin. As such, inhibition of citrullination has been suggested as a therapeutic strategy for MS. Here, in contrast, we show that citrullination by peptidylarginine deiminase 2 (PAD2) contributes to normal oligodendrocyte differentiation, myelination, and motor function. We identify several targets for PAD2, including myelin and chromatin-related proteins, implicating PAD2 in epigenomic regulation. Accordingly, we observe that PAD2 inhibition and its knockdown affect chromatin accessibility and prevent the upregulation of oligodendrocyte differentiation genes. Moreover, mice lacking PAD2 display motor dysfunction and a decreased number of myelinated axons in the corpus callosum. We conclude that citrullination contributes to proper oligodendrocyte lineage progression and myelination.

•

PAD2 is increased upon OL differentiation

•

OL differentiation is facilitated by PAD2-mediated chromatin remodeling in myelin genes

•

PAD2 contributes to efficient myelination and motor and cognitive functions

•

Nuclear and myelin proteins interact and are citrullinated by PAD2

RADICL-seq identifies general and cell type-specific principles of genome-wide RNA-chromatin interactions

Mammalian genomes encode tens of thousands of noncoding RNAs. Most noncoding transcripts exhibit nuclear localization and several have been shown to play a role in the regulation of gene expression and chromatin remodeling. To investigate the function of such RNAs, methods to massively map the genomic interacting sites of multiple transcripts have been developed; however, these methods have some limitations. Here, we introduce RNA And DNA Interacting Complexes Ligated and sequenced (RADICL-seq), a technology that maps genome-wide RNA–chromatin interactions in intact nuclei. RADICL-seq is a proximity ligation-based methodology that reduces the bias for nascent transcription, while increasing genomic coverage and unique mapping rate efficiency compared with existing methods. RADICL-seq identifies distinct patterns of genome occupancy for different classes of transcripts as well as cell type–specific RNA-chromatin interactions, and highlights the role of transcription in the establishment of chromatin structure.

Mammalian genomes encode tens of thousands of ncRNAs that have important roles in regulation of gene expression and chromatin organization. Here, the authors present RADICLseq to map RNA-chromatin interactions in intact nuclei to shed light on these fine-tuned processes.

SETDB1-dependent heterochromatin stimulates alternative lengthening of telomeres

Alternative lengthening of telomeres, or ALT, is a recombination-based process that maintains telomeres to render some cancer cells immortal. The prevailing view is that ALT is inhibited by heterochromatin because heterochromatin prevents recombination. To test this model, we used telomere-specific quantitative proteomics on cells with heterochromatin deficiencies. In contrast to expectations, we found that ALT does not result from a lack of heterochromatin; rather, ALT is a consequence of heterochromatin formation at telomeres, which is seeded by the histone methyltransferase SETDB1. Heterochromatin stimulates transcriptional elongation at telomeres together with the recruitment of recombination factors, while disrupting heterochromatin had the opposite effect. Consistently, loss of SETDB1, disrupts telomeric heterochromatin and abrogates ALT. Thus, inhibiting telomeric heterochromatin formation in ALT cells might offer a new therapeutic approach to cancer treatment.

Disease-specific oligodendrocyte lineage cells arise in multiple sclerosis

Multiple sclerosis (MS) is characterized by an immune system attack targeting myelin, which is produced by oligodendrocytes (OLs). We performed single-cell transcriptomic analysis of OL lineage cells from the spinal cord of mice induced with experimental autoimmune encephalomyelitis (EAE), which mimics several aspects of MS. We found unique OLs and OL precursor cells (OPCs) in EAE and uncovered several genes specifically alternatively spliced in these cells. Surprisingly, EAE-specific OL lineage populations expressed genes involved in antigen processing and presentation via major histocompatibility complex class I and II (MHC-I and -II), and in immunoprotection, suggesting alternative functions of these cells in a disease context. Importantly, we found that disease-specific oligodendroglia are also present in human MS brains and that a substantial number of genes known to be susceptibility genes for MS, so far mainly associated with immune cells, are expressed in the OL lineage cells. Finally, we demonstrate that OPCs can phagocytose and that MHC-II-expressing OPCs can activate memory and effector CD4-positive T cells. Our results suggest that OLs and OPCs are not passive targets but instead active immunomodulators in MS. The disease-specific OL lineage cells, for which we identify several biomarkers, may represent novel direct targets for immunomodulatory therapeutic approaches in MS.

Single-cell transcriptomic analysis of oligodendrocyte lineage cells

Oligodendrocytes (OLs) are glial cells in the central nervous system (CNS), which produce myelin, a lipid-rich membrane that insulates neuronal axons. The main function ascribed to OLs is to regulate the speed of electric pulse transmission, and as such OLs have been widely considered as a single and discrete population. Nevertheless, OLs and their precursor cells (OPCs) throughout the CNS have different morphologies and regional functional differences have been observed. Moreover, OLs have recently been involved in other functional processes such as metabolic coupling with axons. In this review, we focus on recent advances in single-cell transcriptomics suggesting that OLs are more heterogeneous than previously thought, with defined subpopulations and cell states that are associated with different stages of lineage progression and might also represent distinct functional states.

A chromatin code for alternative splicing involving a putative association between CTCF and HP1α proteins

Background

Alternative splicing is primarily controlled by the activity of splicing factors and by the elongation of the RNA polymerase II (RNAPII). Recent experiments have suggested a new complex network of splicing regulation involving chromatin, transcription and multiple protein factors. In particular, the CCCTC-binding factor (CTCF), the Argonaute protein AGO1, and members of the heterochromatin protein 1 (HP1) family have been implicated in the regulation of splicing associated with chromatin and the elongation of RNAPII. These results raise the question of whether these proteins may associate at the chromatin level to modulate alternative splicing.

Results

Using chromatin immunoprecipitation sequencing (ChIP-Seq) data for CTCF, AGO1, HP1α, H3K27me3, H3K9me2, H3K36me3, RNAPII, total H3 and 5metC and alternative splicing arrays from two cell lines, we have analyzed the combinatorial code of their binding to chromatin in relation to the alternative splicing patterns between two cell lines, MCF7 and MCF10. Using Machine Learning techniques, we identified the changes in chromatin signals that are most significantly associated with splicing regulation between these two cell lines. Moreover, we have built a map of the chromatin signals on the pre-mRNA, that is, a chromatin-based RNA-map, which can explain 606 (68.55%) of the regulated events between MCF7 and MCF10. This chromatin code involves the presence of HP1α, CTCF, AGO1, RNAPII and histone marks around regulated exons and can differentiate patterns of skipping and inclusion. Additionally, we found a significant association of HP1α and CTCF activities around the regulated exons and a putative DNA binding site for HP1α.

Conclusions

Our results show that a considerable number of alternative splicing events could have a chromatin-dependent regulation involving the association of HP1α and CTCF near regulated exons. Additionally, we find further evidence for the involvement of HP1α and AGO1 in chromatin-related splicing regulation.

Electronic supplementary material

The online version of this article (doi:10.1186/s12915-015-0141-5) contains supplementary material, which is available to authorized users.

Methods to study splicing from high-throughput RNA sequencing data

The development of novel high-throughput sequencing (HTS) methods for RNA (RNA-Seq) has provided a very powerful mean to study splicing under multiple conditions at unprecedented depth. However, the complexity of the information to be analyzed has turned this into a challenging task. In the last few years, a plethora of tools have been developed, allowing researchers to process RNA-Seq data to study the expression of isoforms and splicing events, and their relative changes under different conditions. We provide an overview of the methods available to study splicing from short RNA-Seq data, which could serve as an entry point for users who need to decide on a suitable tool for a specific analysis. We also attempt to propose a classification of the tools according to the operations they do, to facilitate the comparison and choice of methods.

Databases and resources for human small non-coding RNAs

Recent advances in high-throughput sequencing have facilitated the genome-wide studies of small non-coding RNAs (sRNAs). Numerous studies have highlighted the role of various classes of sRNAs at different levels of gene regulation and disease. The fast growth of sequence data and the diversity of sRNA species have prompted the need to organise them in annotation databases. There are currently several databases that collect sRNA data. Various tools are provided for access, with special emphasis on the well-characterised family of micro-RNAs. The striking heterogeneity of the new classes of sRNAs and the lack of sufficient functional annotation, however, make integration of these datasets a difficult task. This review describes the currently available databases for human sRNAs that are accessible via the internet, and some of the large datasets for human sRNAs from high-throughput sequencing experiments that are so far only available as supplementary data in publications. Some of the main issues related to the integration and annotation of sRNA datasets are also discussed.

Chromatin and alternative splicing

Alternative splicing affects more than 90% of human genes. Coupling between transcription and splicing has become crucial in the complex network underlying alternative splicing regulation. Because chromatin is the real template for nuclear transcription, changes in its structure, but also in the "reading" and "writing" of the histone code, could modulate splicing choices. Here, we discuss the evidence supporting these ideas, from the first proposal of chromatin affecting alternative splicing, performed 20 years ago, to the latest findings including genome-wide evidence that nucleosomes are preferentially positioned in exons. We focus on two recent reports from our laboratories that add new evidence to this field. The first report shows that a physiological stimulus such as neuron depolarization promotes intragenic histone acetylation (H3K9ac) and chromatin relaxation, causing the skipping of exon 18 of the neural cell adhesion molecule gene. In the second report, we show how specific histone modifications can be created at targeted gene regions as a way to affect alternative splicing: Using small interfering RNAs (siRNAs), we increased the levels of H3K9me2 and H3K27me3 in the proximity of alternative exon 33 of the human fibronectin gene, favoring its inclusion into mature messenger RNA (mRNA) through a mechanism that recalls RNA-mediated transcriptional gene silencing.

Control of alternative splicing through siRNA-mediated transcriptional gene silencing

When targeting promoter regions, small interfering RNAs (siRNAs) trigger a previously proposed pathway known as transcriptional gene silencing by promoting heterochromatin formation. Here we show that siRNAs targeting intronic or exonic sequences close to an alternative exon regulate the splicing of that exon. The effect occurred in hepatoma and HeLa cells with siRNA antisense strands designed to enter the silencing pathway, suggesting hybridization with nascent pre-mRNA. Unexpectedly, in HeLa cells the sense strands were also effective, suggesting that an endogenous antisense transcript, detectable in HeLa but not in hepatoma cells, acts as a target. The effect depends on Argonaute-1 and is counterbalanced by factors favoring chromatin opening or transcriptional elongation. The increase in heterochromatin marks (dimethylation at Lys9 and trimethylation at Lys27 of histone H3) at the target site, the need for the heterochromatin-associated protein HP1alpha and the reduction in RNA polymerase II processivity suggest a mechanism involving the kinetic coupling of transcription and alternative splicing.

Co-evolution of the branch site and SR proteins in eukaryotes

Serine-arginine-rich (SR) proteins are essential for splicing in metazoans but are absent in yeast. By contrast, many fungi have SR protein homologs with variable arginine-rich regions analogous to the arginine-serine-rich (RS) domain in metazoans. The density of RS repeats in these regions correlates with the conservation of the branch site signal, providing evidence for an ancestral origin of SR proteins and indicating that the SR proteins and the branch site co-evolved.

On the Use of Automatically Acquired Examples for All-Nouns Word Sense Disambiguation

This article focuses on Word Sense Disambiguation (WSD), which is a Natural Language Processing task that is thought to be important for many Language Technology applications, such as Information Retrieval, Information Extraction, or Machine Translation. One of the main issues preventing the deployment of WSD technology is the lack of training examples for Machine Learning systems, also known as the Knowledge Acquisition Bottleneck. A method which has been shown to work for small samples of words is the automatic acquisition of examples. We have previously shown that one of the most promising example acquisition methods scales up and produces a freely available database of 150 million examples from Web snippets for all polysemous nouns in WordNet. This paper focuses on the issues that arise when using those examples, all alone or in addition to manually tagged examples, to train a supervised WSD system for all nouns. The extensive evaluation on both lexical-sample and all-words Senseval benchmarks shows that we are able to improve over commonly used baselines and to achieve top-rank performance. The good use of the prior distributions from the senses proved to be a crucial factor.