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.

Spatial profiling of chromatin accessibility in mouse and human tissues.

Cellular function in tissue is dependent on the local environment, requiring new methods for spatial mapping of biomolecules and cells in the tissue context1. The emergence of spatial transcriptomics has enabled genome-scale gene expression mapping2-5, but the ability to capture spatial epigenetic information of tissue at the cellular level and genome scale is lacking. Here we describe a method for spatially resolved chromatin accessibility profiling of tissue sections using next-generation sequencing (spatial-ATAC-seq) by combining in situ Tn5 transposition chemistry6 and microfluidic deterministic barcoding5. Profiling mouse embryos using spatial-ATAC-seq delineated tissue-region-specific epigenetic landscapes and identified gene regulators involved in the development of the central nervous system. Mapping the accessible genome in the mouse and human brain revealed the intricate arealization of brain regions. Applying spatial-ATAC-seq to tonsil tissue resolved the spatially distinct organization of immune cell types and states in lymphoid follicles and extrafollicular zones. This technology progresses spatial biology by enabling spatially resolved chromatin accessibility profiling to improve our understanding of cell identity, cell state and cell fate decision in relation to epigenetic underpinnings in development and disease.

Cell-type specialization is encoded by specific chromatin topologies.

The three-dimensional (3D) structure of chromatin is intrinsically associated with gene regulation and cell function1-3. Methods based on chromatin conformation capture have mapped chromatin structures in neuronal systems such as in vitro differentiated neurons, neurons isolated through fluorescence-activated cell sorting from cortical tissues pooled from different animals and from dissociated whole hippocampi4-6. However, changes in chromatin organization captured by imaging, such as the relocation of Bdnf away from the nuclear periphery after activation7, are invisible with such approaches8. Here we developed immunoGAM, an extension of genome architecture mapping (GAM)2,9, to map 3D chromatin topology genome-wide in specific brain cell types, without tissue disruption, from single animals. GAM is a ligation-free technology that maps genome topology by sequencing the DNA content from thin (about 220 nm) nuclear cryosections. Chromatin interactions are identified from the increased probability of co-segregation of contacting loci across a collection of nuclear slices. ImmunoGAM expands the scope of GAM to enable the selection of specific cell types using low cell numbers (approximately 1,000 cells) within a complex tissue and avoids tissue dissociation2,10. We report cell-type specialized 3D chromatin structures at multiple genomic scales that relate to patterns of gene expression. We discover extensive 'melting' of long genes when they are highly expressed and/or have high chromatin accessibility. The contacts most specific of neuron subtypes contain genes associated with specialized processes, such as addiction and synaptic plasticity, which harbour putative binding sites for neuronal transcription factors within accessible chromatin regions. Moreover, sensory receptor genes are preferentially found in heterochromatic compartments in brain cells, which establish strong contacts across tens of megabases. Our results demonstrate that highly specific chromatin conformations in brain cells are tightly related to gene regulation mechanisms and specialized functions.

Altered human oligodendrocyte heterogeneity in multiple sclerosis.

Oligodendrocyte pathology is increasingly implicated in neurodegenerative diseases as oligodendrocytes both myelinate and provide metabolic support to axons. In multiple sclerosis (MS), demyelination in the central nervous system thus leads to neurodegeneration, but the severity of MS between patients is very variable. Disability does not correlate well with the extent of demyelination1, which suggests that other factors contribute to this variability. One such factor may be oligodendrocyte heterogeneity. Not all oligodendrocytes are the same-those from the mouse spinal cord inherently produce longer myelin sheaths than those from the cortex2, and single-cell analysis of the mouse central nervous system identified further differences3,4. However, the extent of human oligodendrocyte heterogeneity and its possible contribution to MS pathology remain unknown. Here we performed single-nucleus RNA sequencing from white matter areas of post-mortem human brain from patients with MS and from unaffected controls. We identified subclusters of oligodendroglia in control human white matter, some with similarities to mouse, and defined new markers for these cell states. Notably, some subclusters were underrepresented in MS tissue, whereas others were more prevalent. These differences in mature oligodendrocyte subclusters may indicate different functional states of oligodendrocytes in MS lesions. We found similar changes in normal-appearing white matter, showing that MS is a more diffuse disease than its focal demyelination suggests. Our findings of an altered oligodendroglial heterogeneity in MS may be important for understanding disease progression and developing therapeutic approaches.

RNA velocity of single cells.

RNA abundance is a powerful indicator of the state of individual cells. Single-cell RNA sequencing can reveal RNA abundance with high quantitative accuracy, sensitivity and throughput1. However, this approach captures only a static snapshot at a point in time, posing a challenge for the analysis of time-resolved phenomena such as embryogenesis or tissue regeneration. Here we show that RNA velocity-the time derivative of the gene expression state-can be directly estimated by distinguishing between unspliced and spliced mRNAs in common single-cell RNA sequencing protocols. RNA velocity is a high-dimensional vector that predicts the future state of individual cells on a timescale of hours. We validate its accuracy in the neural crest lineage, demonstrate its use on multiple published datasets and technical platforms, reveal the branching lineage tree of the developing mouse hippocampus, and examine the kinetics of transcription in human embryonic brain. We expect RNA velocity to greatly aid the analysis of developmental lineages and cellular dynamics, particularly in humans.

Crossing boundaries: Interplay between the immune system and oligodendrocyte lineage cells.

Oligodendrocytes and their progenitors are glial cells in the central nervous system, which have been mainly implicated with the homeostatic roles of axonal myelin ensheathment but serve as targets of the peripheral immune system attack in the context of diseases like multiple sclerosis. This view of oligodendroglia as passive bystanders with no immunological properties was first challenged in the 1980s when it was reported that the cytokine interferon γ could induce the gene expression of the major histocompatibility complexes (MHC) class I and II. While the physiological role of this induction was controversial for decades to follow, recent studies suggest that oligodendroglia survey their environment, respond to a larger array of cues and can indeed exert immunomodulatory functions, which are particularly relevant in the context of neurodegeneration and demyelinating diseases. The alternative functionality of oligodendroglia not only regulates immune cell responses, but also hinders remyelination, and might thereby be key to understanding MS disease pathology and promoting regeneration after immune-mediated demyelination.

Spatial and temporal heterogeneity in the lineage progression of fine oligodendrocyte subtypes.

BACKGROUND: Oligodendrocytes are glial cells that support and insulate axons in the central nervous system through the production of myelin. Oligodendrocytes arise throughout embryonic and early postnatal development from oligodendrocyte precursor cells (OPCs), and recent work demonstrated that they are a transcriptional heterogeneous cell population, but the regional and functional implications of this heterogeneity are less clear. Here, we apply in situ sequencing (ISS) to simultaneously probe the expression of 124 marker genes of distinct oligodendrocyte populations, providing comprehensive maps of the corpus callosum, cingulate, motor, and somatosensory cortex in the brain, as well as gray matter (GM) and white matter (WM) regions in the spinal cord, at postnatal (P10), juvenile (P20), and young adult (P60) stages. We systematically compare the abundances of these populations and investigate the neighboring preference of distinct oligodendrocyte populations. RESULTS: We observed that oligodendrocyte lineage progression is more advanced in the juvenile spinal cord compared to the brain, corroborating with previous studies. We found myelination still ongoing in the adult corpus callosum while it was more advanced in the cortex. Interestingly, we also observed a lateral-to-medial gradient of oligodendrocyte lineage progression in the juvenile cortex, which could be linked to arealization, as well as a deep-to-superficial gradient with mature oligodendrocytes preferentially accumulating in the deeper layers of the cortex. The ISS experiments also exposed differences in abundances and population dynamics over time between GM and WM regions in the brain and spinal cord, indicating regional differences within GM and WM, and we found that neighboring preferences of some oligodendroglia populations are altered from the juvenile to the adult CNS. CONCLUSIONS: Overall, our ISS experiments reveal spatial heterogeneity of oligodendrocyte lineage progression in the brain and spinal cord and uncover differences in the timing of oligodendrocyte differentiation and myelination, which could be relevant to further investigate functional heterogeneity of oligodendroglia, especially in the context of injury or disease.

Birth, coming of age and death: The intriguing life of long noncoding RNAs.

Mammalian genomes are pervasively transcribed, with long noncoding RNAs being the most abundant fraction. Recent studies have highlighted the central role played by these transcripts in several physiological and pathological processes. Despite several metabolic features shared between coding and noncoding transcripts, these two classes of RNAs exhibit multiple differences regarding their biogenesis and processing. Here we review such distinctions, focusing on the unique features of specific long noncoding RNAs.

Epigenetic regulation of oligodendrocyte differentiation: From development to demyelinating disorders.

The maintenance of progenitor states or the differentiation of progenitors into specific lineages requires epigenetic remodeling of the gene expression program. In the central nervous system, oligodendrocyte progenitors (OPCs) give rise to oligodendrocytes (OLs), whose main function has been thought to be to produce myelin, a lipid-rich structure insulating the axons. However, recent findings suggest diverse OL transcriptional states, which might imply additional functions. The differentiation of OPCs into postmitotic OLs is a highly regulated and sensitive process and requires temporal waves of gene expression through epigenetic remodeling of the genome. In this review, we will discuss recent advances in understanding the events shaping the chromatin landscape through histone modifications and long noncoding RNAs during OPC differentiation, in physiological and pathological conditions. We suggest that epigenetic regulation plays a fundamental role in governing the accessibility of transcriptional machinery to DNA sequences, which ultimately determines functional outcomes in OLs.

Interaction of Sox2 with RNA binding proteins in mouse embryonic stem cells.

Sox2 is a master transcriptional regulator of embryonic development. In this study, we determined the protein interactome of Sox2 in the chromatin and nucleoplasm of mouse embryonic stem (mES) cells. Apart from canonical interactions with pluripotency-regulating transcription factors, we identified interactions with several chromatin modulators, including members of the heterochromatin protein 1 (HP1) family, suggesting a role for Sox2 in chromatin-mediated transcriptional repression. Sox2 was also found to interact with RNA binding proteins (RBPs), including proteins involved in RNA processing. RNA immunoprecipitation followed by sequencing revealed that Sox2 associates with different messenger RNAs, as well as small nucleolar RNA Snord34 and the non-coding RNA 7SK. 7SK has been shown to regulate transcription at gene regulatory regions, which could suggest a functional interaction with Sox2 for chromatin recruitment. Nevertheless, we found no evidence of Sox2 modulating recruitment of 7SK to chromatin when examining 7SK chromatin occupancy by Chromatin Isolation by RNA Purification (ChIRP) in Sox2 depleted mES cells. In addition, knockdown of 7SK in mES cells did not lead to any change in Sox2 occupancy at 7SK-regulated genes. Thus, our results show that Sox2 extensively interacts with RBPs, and suggest that Sox2 and 7SK co-exist in a ribonucleoprotein complex whose function is not to regulate chromatin recruitment, but could rather regulate other processes in the nucleoplasm.

Citrullination regulates pluripotency and histone H1 binding to chromatin.

Citrullination is the post-translational conversion of an arginine residue within a protein to the non-coded amino acid citrulline. This modification leads to the loss of a positive charge and reduction in hydrogen-bonding ability. It is carried out by a small family of tissue-specific vertebrate enzymes called peptidylarginine deiminases (PADIs) and is associated with the development of diverse pathological states such as autoimmunity, cancer, neurodegenerative disorders, prion diseases and thrombosis. Nevertheless, the physiological functions of citrullination remain ill-defined, although citrullination of core histones has been linked to transcriptional regulation and the DNA damage response. PADI4 (also called PAD4 or PADV), the only PADI with a nuclear localization signal, was previously shown to act in myeloid cells where it mediates profound chromatin decondensation during the innate immune response to infection. Here we show that the expression and enzymatic activity of Padi4 are also induced under conditions of ground-state pluripotency and during reprogramming in mouse. Padi4 is part of the pluripotency transcriptional network, binding to regulatory elements of key stem-cell genes and activating their expression. Its inhibition lowers the percentage of pluripotent cells in the early mouse embryo and significantly reduces reprogramming efficiency. Using an unbiased proteomic approach we identify linker histone H1 variants, which are involved in the generation of compact chromatin, as novel PADI4 substrates. Citrullination of a single arginine residue within the DNA-binding site of H1 results in its displacement from chromatin and global chromatin decondensation. Together, these results uncover a role for citrullination in the regulation of pluripotency and provide new mechanistic insights into how citrullination regulates chromatin compaction.

Positional differences of axon growth rates between sensory neurons encoded by Runx3.

The formation of functional connectivity in the nervous system is governed by axon guidance that instructs nerve growth and branching during development, implying a similarity between neuronal subtypes in terms of nerve extension. We demonstrate the molecular mechanism of another layer of complexity in vertebrates by defining a transcriptional program underlying growth differences between positionally different neurons. The rate of axon extension of the early subset of embryonic dorsal root ganglion sensory neurons is encoded in neurons at different axial levels. This code is determined by a segmental pattern of axial levels of Runx family transcription factor Runx3. Runx3 in turn determines transcription levels of genes encoding cytoskeletal proteins involved in axon extension, including Rock1 and Rock2 which have ongoing activities determining axon growth in early sensory neurons and blocking Rock activity reverses axon extension deficits of Runx3(-/-) neurons. Thus, Runx3 acts to regulate positional differences in axon extension properties apparently without affecting nerve guidance and branching, a principle that could be relevant to other parts of the nervous system.

Acute treatment with valproic acid and l-thyroxine ameliorates clinical signs of experimental autoimmune encephalomyelitis and prevents brain pathology in DA rats.

Multiple sclerosis (MS) is the most common chronic inflammatory demyelinating disease of the central nervous system (CNS) in young adults. Chronic treatments with histone deacetylase inhibitors (HDACis) have been reported to ameliorate experimental autoimmune encephalomyelitis (EAE), a rodent model of MS, by targeting immune responses. We have recently shown that the HDAC inhibition/knockdown in the presence of thyroid hormone (T3) can also promote oligodendrocyte (OL) differentiation and expression of myelin genes in neural stem cells (NSCs) and oligodendrocyte precursors (OPCs). In this study, we found that treatment with an HDACi, valproic acid (VPA), and T3, alone or in combination, directly affects encephalitogenic CD4+ T cells. VPA, but not T3, compromised their proliferation, while both molecules reduced the frequency of IL-17-producing cells. Transfer of T3, VPA and VPA/T3 treated encephalitogenic CD4+ T cells into naïve rats induced less severe EAE, indicating that the effects of these molecules are persistent and do not require their maintenance after the initial stimuli. Thus, we investigated the effect of acute treatment with VPA and l-thyroxine (T4), a precursor of T3, on myelin oligodendrocyte glycoprotein-induced EAE in Dark Agouti rats, a close mimic of MS. We found that a brief treatment after disease onset led to sustained amelioration of EAE and prevention of inflammatory demyelination in the CNS accompanied with a higher expression of myelin-related genes in the brain. Furthermore, the treatment modulated immune responses, reduced the number of CD4+ T cells and affected the Th1 differentiation program in the brain. Our data indicate that an acute treatment with VPA and T4 after the onset of EAE can produce persistent clinically relevant therapeutic effects by limiting the pathogenic immune reactions while promoting myelin gene expression.

Histone H2AX-dependent GABA(A) receptor regulation of stem cell proliferation.

Stem cell self-renewal implies proliferation under continued maintenance of multipotency. Small changes in numbers of stem cells may lead to large differences in differentiated cell numbers, resulting in significant physiological consequences. Proliferation is typically regulated in the G1 phase, which is associated with differentiation and cell cycle arrest. However, embryonic stem (ES) cells may lack a G1 checkpoint. Regulation of proliferation in the 'DNA damage' S/G2 cell cycle checkpoint pathway is known for its role in the maintenance of chromatin structural integrity. Here we show that autocrine/paracrine gamma-aminobutyric acid (GABA) signalling by means of GABA(A) receptors negatively controls ES cell and peripheral neural crest stem (NCS) cell proliferation, preimplantation embryonic growth and proliferation in the boundary-cap stem cell niche, resulting in an attenuation of neuronal progenies from this stem cell niche. Activation of GABA(A) receptors leads to hyperpolarization, increased cell volume and accumulation of stem cells in S phase, thereby causing a rapid decrease in cell proliferation. GABA(A) receptors signal through S-phase checkpoint kinases of the phosphatidylinositol-3-OH kinase-related kinase family and the histone variant H2AX. This signalling pathway critically regulates proliferation independently of differentiation, apoptosis and overt damage to DNA. These results indicate the presence of a fundamentally different mechanism of proliferation control in these stem cells, in comparison with most somatic cells, involving proteins in the DNA damage checkpoint pathway.

Nano-CUT&Tag for multimodal chromatin profiling at single-cell resolution.

The ability to comprehensively analyze the chromatin state with single-cell resolution is crucial for understanding gene regulatory principles in heterogenous tissues or during development. Recently, we developed a nanobody-based single-cell CUT&Tag (nano-CT) protocol to simultaneously profile three epigenetic modalities-two histone marks and open chromatin state-from the same single cell. Nano-CT implements a new set of secondary nanobody-Tn5 fusion proteins to direct barcoded tagmentation by Tn5 transposase to genomic targets labeled by primary antibodies raised in different species. Such nanobody-Tn5 fusion proteins are currently not commercially available, and their in-house production and purification can be completed in 3-4 d by following our detailed protocol. The single-cell indexing in nano-CT is performed on a commercially available platform, making it widely accessible to the community. In comparison to other multimodal methods, nano-CT stands out in data complexity, low sample requirements and the flexibility to choose two of the three modalities. In addition, nano-CT works efficiently with fresh brain samples, generating multimodal epigenomic profiles for thousands of brain cells at single-cell resolution. The nano-CT protocol can be completed in just 3 d by users with basic skills in standard molecular biology and bioinformatics, although previous experience with single-cell assay for transposase-accessible chromatin using sequencing (scATAC-seq) is beneficial for more in-depth data analysis. As a multimodal assay, nano-CT holds immense potential to reveal interactions of various chromatin modalities, to explore epigenetic heterogeneity and to increase our understanding of the role and interplay that chromatin dynamics has in cellular development.

Spatial dynamics of mammalian brain development and neuroinflammation by multimodal tri-omics mapping.

The ability to spatially map multiple layers of the omics information over different time points allows for exploring the mechanisms driving brain development, differentiation, arealization, and alterations in disease. Herein we developed and applied spatial tri-omic sequencing technologies, DBiT ARP-seq (spatial ATAC-RNA-Protein-seq) and DBiT CTRP-seq (spatial CUT&Tag- RNA-Protein-seq) together with multiplexed immunofluorescence imaging (CODEX) to map spatial dynamic remodeling in brain development and neuroinflammation. A spatiotemporal tri-omic atlas of the mouse brain was obtained at different stages from postnatal day P0 to P21, and compared to the regions of interest in the human developing brains. Specifically, in the cortical area, we discovered temporal persistence and spatial spreading of chromatin accessibility for the layer-defining transcription factors. In corpus callosum, we observed dynamic chromatin priming of myelin genes across the subregions. Together, it suggests a role for layer specific projection neurons to coordinate axonogenesis and myelination. We further mapped the brain of a lysolecithin (LPC) neuroinflammation mouse model and observed common molecular programs in development and neuroinflammation. Microglia, exhibiting both conserved and distinct programs for inflammation and resolution, are transiently activated not only at the core of the LPC lesion, but also at distal locations presumably through neuronal circuitry. Thus, this work unveiled common and differential mechanisms in brain development and neuroinflammation, resulting in a valuable data resource to investigate brain development, function and disease.

Developmental origin of oligodendrocytes determines their function in the adult brain.

In the mouse embryonic forebrain, developmentally distinct oligodendrocyte progenitor cell populations and their progeny, oligodendrocytes, emerge from three distinct regions in a spatiotemporal gradient from ventral to dorsal. However, the functional importance of this oligodendrocyte developmental heterogeneity is unknown. Using a genetic strategy to ablate dorsally derived oligodendrocyte lineage cells (OLCs), we show here that the areas in which dorsally derived OLCs normally reside in the adult central nervous system become populated and myelinated by OLCs of ventral origin. These ectopic oligodendrocytes (eOLs) have a distinctive gene expression profile as well as subtle myelination abnormalities. The failure of eOLs to fully assume the role of the original dorsally derived cells results in locomotor and cognitive deficits in the adult animal. This study reveals the importance of developmental heterogeneity within the oligodendrocyte lineage and its importance for homeostatic brain function.