Stem cell-derived brainstem mouse astrocytes obtain a neurotoxic phenotype in vitro upon neuroinflammation.

BACKGROUND: Astrocytes respond to injury and disease through a process known as reactive astrogliosis, of which inflammatory signaling is one subset. This inflammatory response is heterogeneous with respect to the inductive stimuli and the afflicted central nervous system region. This is of plausible importance in e.g. traumatic axonal injury (TAI), where lesions in the brainstem carries a particularly poor prognosis. In fact, astrogliotic forebrain astrocytes were recently suggested to cause neuronal death following axotomy. We therefore sought to assess if ventral brainstem- or rostroventral spinal astrocytes exert similar effects on motor neurons in vitro. METHODS: We derived brainstem/rostroventral spinal astrocyte-like cells (ES-astrocytes) and motor neurons using directed differentiation of mouse embryonic stem cells (ES). We activated the ES-astrocytes using the neurotoxicity-eliciting cytokines interleukin- (IL-) 1α and tumor necrosis factor-(TNF-)α and clinically relevant inflammatory mediators. In co-cultures with reactive ES-astrocytes and motor neurons, we assessed neurotoxic ES-astrocyte activity, similarly to what has previously been shown for other central nervous system (CNS) regions. RESULTS: We confirmed the brainstem/rostroventral ES-astrocyte identity using RNA-sequencing, immunocytochemistry, and by comparison with primary subventricular zone-astrocytes. Following cytokine stimulation, the c-Jun N-terminal kinase pathway down-stream product phosphorylated c-Jun was increased, thus demonstrating ES-astrocyte reactivity. These reactive ES-astrocytes conferred a contact-dependent neurotoxic effect upon co-culture with motor neurons. When exposed to IL-1ß and IL-6, two neuroinflammatory cytokines found in the cerebrospinal fluid and serum proteome following human severe traumatic brain injury (TBI), ES-astrocytes exerted similar effects on motor neurons. Activation of ES-astrocytes by these cytokines was associated with pathways relating to endoplasmic reticulum stress and altered regulation of MYC. CONCLUSIONS: Ventral brainstem and rostroventral spinal cord astrocytes differentiated from mouse ES can exert neurotoxic effects in vitro. This highlights how neuroinflammation following CNS lesions can exert region- and cell-specific effects. Our in vitro model system, which uniquely portrays astrocytes and neurons from one niche, allows for a detailed and translationally relevant model system for future studies on how to improve neuronal survival in particularly vulnerable CNS regions following e.g. TAI.

Single-cell transcriptome and accessible chromatin dynamics during endocrine pancreas development.

Delineating gene regulatory networks that orchestrate cell-type specification is a continuing challenge for developmental biologists. Single-cell analyses offer opportunities to address these challenges and accelerate discovery of rare cell lineage relationships and mechanisms underlying hierarchical lineage decisions. Here, we describe the molecular analysis of mouse pancreatic endocrine cell differentiation using single-cell transcriptomics, chromatin accessibility assays coupled to genetic labeling, and cytometry-based cell purification. We uncover transcription factor networks that delineate ß-, α-, and δ-cell lineages. Through genomic footprint analysis, we identify transcription factor-regulatory DNA interactions governing pancreatic cell development at unprecedented resolution. Our analysis suggests that the transcription factor Neurog3 may act as a pioneer transcription factor to specify the pancreatic endocrine lineage. These findings could improve protocols to generate replacement endocrine cells from renewable sources, like stem cells, for diabetes therapy.

Elastic dosage compensation by X-chromosome upregulation.

X-chromosome inactivation and X-upregulation are the fundamental modes of chromosome-wide gene regulation that collectively achieve dosage compensation in mammals, but the regulatory link between the two remains elusive and the X-upregulation dynamics are unknown. Here, we use allele-resolved single-cell RNA-seq combined with chromatin accessibility profiling and finely dissect their separate effects on RNA levels during mouse development. Surprisingly, we uncover that X-upregulation elastically tunes expression dosage in a sex- and lineage-specific manner, and moreover along varying degrees of X-inactivation progression. Male blastomeres achieve X-upregulation upon zygotic genome activation while females experience two distinct waves of upregulation, upon imprinted and random X-inactivation; and ablation of Xist impedes female X-upregulation. Female cells carrying two active X chromosomes lack upregulation, yet their collective RNA output exceeds that of a single hyperactive allele. Importantly, this conflicts the conventional dosage compensation model in which naïve female cells are initially subject to biallelic X-upregulation followed by X-inactivation of one allele to correct the X dosage. Together, our study provides key insights to the chain of events of dosage compensation, explaining how transcript copy numbers can remain remarkably stable across developmental windows wherein severe dose imbalance would otherwise be experienced by the cell.

An overlooked subset of Cx3cr1wt/wt microglia in the Cx3cr1CreER-Eyfp/wt mouse has a repopulation advantage over Cx3cr1CreER-Eyfp/wt microglia following microglial depletion.

BACKGROUND: Fluorescent reporter labeling and promoter-driven Cre-recombinant technologies have facilitated cellular investigations of physiological and pathological processes, including the widespread use of the Cx3cr1CreER-Eyfp/wt mouse strain for studies of microglia. METHODS: Immunohistochemistry, Flow Cytometry, RNA sequencing and whole-genome sequencing were used to identify the subpopulation of microglia in Cx3cr1CreER-Eyfp/wt mouse brains. Genetically mediated microglia depletion using Cx3cr1CreER-Eyfp/wtRosa26DTA/wt mice and CSF1 receptor inhibitor PLX3397 were used to deplete microglia. Primary microglia proliferation and migration assay were used for in vitro studies. RESULTS: We unexpectedly identified a subpopulation of microglia devoid of genetic modification, exhibiting higher Cx3cr1 and CX3CR1 expression than Cx3cr1CreER-Eyfp/wtCre+Eyfp+ microglia in Cx3cr1CreER-Eyfp/wt mouse brains, thus termed Cx3cr1highCre-Eyfp- microglia. This subpopulation constituted less than 1% of all microglia under homeostatic conditions, but after Cre-driven DTA-mediated microglial depletion, Cx3cr1highCre-Eyfp- microglia escaped depletion and proliferated extensively, eventually occupying one-third of the total microglial pool. We further demonstrated that the Cx3cr1highCre-Eyfp- microglia had lost their genetic heterozygosity and become homozygous for wild-type Cx3cr1. Therefore, Cx3cr1highCre-Eyfp- microglia are Cx3cr1wt/wtCre-Eyfp-. Finally, we demonstrated that CX3CL1-CX3CR1 signaling regulates microglial repopulation both in vivo and in vitro. CONCLUSIONS: Our results raise a cautionary note regarding the use of Cx3cr1CreER-Eyfp/wt mouse strains, particularly when interpreting the results of fate mapping, and microglial depletion and repopulation studies.

Spinal Cord Injury Induces Permanent Reprogramming of Microglia into a Disease-Associated State Which Contributes to Functional Recovery.

Microglia are resident myeloid cells of the CNS. Recently, single-cell RNA sequencing (scRNAseq) has enabled description of a disease-associated microglia (DAM) with a role in neurodegeneration and demyelination. In this study, we use scRNAseq to investigate the temporal dynamics of immune cells harvested from the epicenter of traumatic spinal cord injury (SCI) induced in female mice. We find that as a consequence of SCI, baseline microglia undergo permanent transcriptional reprogramming into a previously uncharacterized subtype of microglia with striking similarities to previously reported DAM as well as a distinct microglial state found during development. Using a microglia depletion model we showed that DAM in SCI are derived from baseline microglia and strongly enhance recovery of hindlimb locomotor function following injury.SIGNIFICANCE STATEMENT Although disease-associated microglia (DAM) have been the subject of strong research interest during recent years (Keren-Shaul, 2017; Jordão, 2019), their cellular origin and their role in "normal" acute injury processes is not well understood. Our work directly addresses the origin and the role of DAM in traumatic injury response. Further, we use a microglia depletion model to prove that DAM in spinal cord injury (SCI) are indeed derived from homeostatic microglia, and that they strongly enhance recovery. Thus, in this work we significantly expand the knowledge of immune response to traumatic injury, demonstrate the applicability to human injury via our unique access to injured human spinal cord tissue, and provide the community with a comprehensive dataset for further exploration.

An unsupervised method for physical cell interaction profiling of complex tissues.

Cellular identity in complex multicellular organisms is determined in part by the physical organization of cells. However, large-scale investigation of the cellular interactome remains technically challenging. Here we develop cell interaction by multiplet sequencing (CIM-seq), an unsupervised and high-throughput method to analyze direct physical cell-cell interactions between cell types present in a tissue. CIM-seq is based on RNA sequencing of incompletely dissociated cells, followed by computational deconvolution into constituent cell types. CIM-seq estimates parameters such as number of cells and cell types in each multiplet directly from sequencing data, making it compatible with high-throughput droplet-based methods. When applied to gut epithelium or whole dissociated lung and spleen, CIM-seq correctly identifies known interactions, including those between different cell lineages and immune cells. In the colon, CIM-seq identifies a previously unrecognized goblet cell subtype expressing the wound-healing marker Plet1, which is directly adjacent to colonic stem cells. Our results demonstrate that CIM-seq is broadly applicable to unsupervised profiling of cell-type interactions in different tissue types.

Tuning MPL signaling to influence hematopoietic stem cell differentiation and inhibit essential thrombocythemia progenitors.

Thrombopoietin (TPO) and the TPO-receptor (TPO-R, or c-MPL) are essential for hematopoietic stem cell (HSC) maintenance and megakaryocyte differentiation. Agents that can modulate TPO-R signaling are highly desirable for both basic research and clinical utility. We developed a series of surrogate protein ligands for TPO-R, in the form of diabodies (DBs), that homodimerize TPO-R on the cell surface in geometries that are dictated by the DB receptor binding epitope, in effect "tuning" downstream signaling responses. These surrogate ligands exhibit diverse pharmacological properties, inducing graded signaling outputs, from full to partial TPO agonism, thus decoupling the dual functions of TPO/TPO-R. Using single-cell RNA sequencing and HSC self-renewal assays we find that partial agonistic diabodies preserved the stem-like properties of cultured HSCs, but also blocked oncogenic colony formation in essential thrombocythemia (ET) through inverse agonism. Our data suggest that dampening downstream TPO signaling is a powerful approach not only for HSC preservation in culture, but also for inhibiting oncogenic signaling through the TPO-R.

A Highly Scalable Method for Joint Whole-Genome Sequencing and Gene-Expression Profiling of Single Cells.

To address how genetic variation alters gene expression in complex cell mixtures, we developed direct nuclear tagmentation and RNA sequencing (DNTR-seq), which enables whole-genome and mRNA sequencing jointly in single cells. DNTR-seq readily identified minor subclones within leukemia patients. In a large-scale DNA damage screen, DNTR-seq was used to detect regions under purifying selection and identified genes where mRNA abundance was resistant to copy-number alteration, suggesting strong genetic compensation. mRNA sequencing (mRNA-seq) quality equals RNA-only methods, and the low positional bias of genomic libraries allowed detection of sub-megabase aberrations at ultra-low coverage. Each cell library is individually addressable and can be re-sequenced at increased depth, allowing multi-tiered study designs. Additionally, the direct tagmentation protocol enables coverage-independent estimation of ploidy, which can be used to identify cell singlets. Thus, DNTR-seq directly links each cell's state to its corresponding genome at scale, enabling routine analysis of heterogeneous tumors and other complex tissues.

Modeling Spatial Correlation of Transcripts with Application to Developing Pancreas.

Recently high-throughput image-based transcriptomic methods were developed and enabled researchers to spatially resolve gene expression variation at the molecular level for the first time. In this work, we develop a general analysis tool to quantitatively study the spatial correlations of gene expression in fixed tissue sections. As an illustration, we analyze the spatial distribution of single mRNA molecules measured by in situ sequencing on human fetal pancreas at three developmental time points-80, 87 and 117 days post-fertilization. We develop a density profile-based method to capture the spatial relationship between gene expression and other morphological features of the tissue sample such as position of nuclei and endocrine cells of the pancreas. In addition, we build a statistical model to characterize correlations in the spatial distribution of the expression level among different genes. This model enables us to infer the inhibitory and clustering effects throughout different time points. Our analysis framework is applicable to a wide variety of spatially-resolved transcriptomic data to derive biological insights.

Single-Cell Analysis of Human Pancreas Reveals Transcriptional Signatures of Aging and Somatic Mutation Patterns.

As organisms age, cells accumulate genetic and epigenetic errors that eventually lead to impaired organ function or catastrophic transformation such as cancer. Because aging reflects a stochastic process of increasing disorder, cells in an organ will be individually affected in different ways, thus rendering bulk analyses of postmitotic adult cells difficult to interpret. Here, we directly measure the effects of aging in human tissue by performing single-cell transcriptome analysis of 2,544 human pancreas cells from eight donors spanning six decades of life. We find that islet endocrine cells from older donors display increased levels of transcriptional noise and potential fate drift. By determining the mutational history of individual cells, we uncover a novel mutational signature in healthy aging endocrine cells. Our results demonstrate the feasibility of using single-cell RNA sequencing (RNA-seq) data from primary cells to derive insights into genetic and transcriptional processes that operate on aging human tissue.

Converting Adult Pancreatic Islet α Cells into ß Cells by Targeting Both Dnmt1 and Arx.

Insulin-producing pancreatic ß cells in mice can slowly regenerate from glucagon-producing α cells in settings like ß cell loss, but the basis of this conversion is unknown. Moreover, it remains unclear if this intra-islet cell conversion is relevant to diseases like type 1 diabetes (T1D). We show that the α cell regulators Aristaless-related homeobox (Arx) and DNA methyltransferase 1 (Dnmt1) maintain α cell identity in mice. Within 3 months of Dnmt1 and Arx loss, lineage tracing and single-cell RNA sequencing revealed extensive α cell conversion into progeny resembling native ß cells. Physiological studies demonstrated that converted α cells acquire hallmark ß cell electrophysiology and show glucose-stimulated insulin secretion. In T1D patients, subsets of glucagon-expressing cells show loss of DNMT1 and ARX and produce insulin and other ß cell factors, suggesting that DNMT1 and ARX maintain α cell identity in humans. Our work reveals pathways regulated by Arx and Dnmt1 that are sufficient for achieving targeted generation of ß cells from adult pancreatic α cells.

DNA-dependent formation of transcription factor pairs alters their binding specificity.

Gene expression is regulated by transcription factors (TFs), proteins that recognize short DNA sequence motifs. Such sequences are very common in the human genome, and an important determinant of the specificity of gene expression is the cooperative binding of multiple TFs to closely located motifs. However, interactions between DNA-bound TFs have not been systematically characterized. To identify TF pairs that bind cooperatively to DNA, and to characterize their spacing and orientation preferences, we have performed consecutive affinity-purification systematic evolution of ligands by exponential enrichment (CAP-SELEX) analysis of 9,400 TF-TF-DNA interactions. This analysis revealed 315 TF-TF interactions recognizing 618 heterodimeric motifs, most of which have not been previously described. The observed cooperativity occurred promiscuously between TFs from diverse structural families. Structural analysis of the TF pairs, including a novel crystal structure of MEIS1 and DLX3 bound to their identified recognition site, revealed that the interactions between the TFs were predominantly mediated by DNA. Most TF pair sites identified involved a large overlap between individual TF recognition motifs, and resulted in recognition of composite sites that were markedly different from the individual TF's motifs. Together, our results indicate that the DNA molecule commonly plays an active role in cooperative interactions that define the gene regulatory lexicon.

A survey of human brain transcriptome diversity at the single cell level.

The human brain is a tissue of vast complexity in terms of the cell types it comprises. Conventional approaches to classifying cell types in the human brain at single cell resolution have been limited to exploring relatively few markers and therefore have provided a limited molecular characterization of any given cell type. We used single cell RNA sequencing on 466 cells to capture the cellular complexity of the adult and fetal human brain at a whole transcriptome level. Healthy adult temporal lobe tissue was obtained during surgical procedures where otherwise normal tissue was removed to gain access to deeper hippocampal pathology in patients with medical refractory seizures. We were able to classify individual cells into all of the major neuronal, glial, and vascular cell types in the brain. We were able to divide neurons into individual communities and show that these communities preserve the categorization of interneuron subtypes that is typically observed with the use of classic interneuron markers. We then used single cell RNA sequencing on fetal human cortical neurons to identify genes that are differentially expressed between fetal and adult neurons and those genes that display an expression gradient that reflects the transition between replicating and quiescent fetal neuronal populations. Finally, we observed the expression of major histocompatibility complex type I genes in a subset of adult neurons, but not fetal neurons. The work presented here demonstrates the applicability of single cell RNA sequencing on the study of the adult human brain and constitutes a first step toward a comprehensive cellular atlas of the human brain.

CTCF/cohesin-binding sites are frequently mutated in cancer.

Cohesin is present in almost all active enhancer regions, where it is associated with transcription factors. Cohesin frequently colocalizes with CTCF (CCCTC-binding factor), affecting genomic stability, expression and epigenetic homeostasis. Cohesin subunits are mutated in cancer, but CTCF/cohesin-binding sites (CBSs) in DNA have not been examined for mutations. Here we report frequent mutations at CBSs in cancers displaying a mutational signature where mutations in A•T base pairs predominate. Integration of whole-genome sequencing data from 213 colorectal cancer (CRC) samples and chromatin immunoprecipitation sequencing (ChIP-exo) data identified frequent point mutations at CBSs. In contrast, CRCs showing an ultramutator phenotype caused by defects in the exonuclease domain of DNA polymerase ɛ (POLE) displayed significantly fewer mutations at and adjacent to CBSs. Analysis of public data showed that multiple cancer types accumulate CBS mutations. CBSs are a major mutational hotspot in the noncoding cancer genome.

Capture Hi-C identifies the chromatin interactome of colorectal cancer risk loci.

Multiple regulatory elements distant from their targets on the linear genome can influence the expression of a single gene through chromatin looping. Chromosome conformation capture implemented in Hi-C allows for genome-wide agnostic characterization of chromatin contacts. However, detection of functional enhancer-promoter interactions is precluded by its effective resolution that is determined by both restriction fragmentation and sensitivity of the experiment. Here we develop a capture Hi-C (cHi-C) approach to allow an agnostic characterization of these physical interactions on a genome-wide scale. Single-nucleotide polymorphisms associated with complex diseases often reside within regulatory elements and exert effects through long-range regulation of gene expression. Applying this cHi-C approach to 14 colorectal cancer risk loci allows us to identify key long-range chromatin interactions in cis and trans involving these loci.

Transcription factor binding in human cells occurs in dense clusters formed around cohesin anchor sites.

During cell division, transcription factors (TFs) are removed from chromatin twice, during DNA synthesis and during condensation of chromosomes. How TFs can efficiently find their sites following these stages has been unclear. Here, we have analyzed the binding pattern of expressed TFs in human colorectal cancer cells. We find that binding of TFs is highly clustered and that the clusters are enriched in binding motifs for several major TF classes. Strikingly, almost all clusters are formed around cohesin, and loss of cohesin decreases both DNA accessibility and binding of TFs to clusters. We show that cohesin remains bound in S phase, holding the nascent sister chromatids together at the TF cluster sites. Furthermore, cohesin remains bound to the cluster sites when TFs are evicted in early M phase. These results suggest that cohesin-binding functions as a cellular memory that promotes re-establishment of TF clusters after DNA replication and chromatin condensation.

DNA-binding specificities of human transcription factors.

Although the proteins that read the gene regulatory code, transcription factors (TFs), have been largely identified, it is not well known which sequences TFs can recognize. We have analyzed the sequence-specific binding of human TFs using high-throughput SELEX and ChIP sequencing. A total of 830 binding profiles were obtained, describing 239 distinctly different binding specificities. The models represent the majority of human TFs, approximately doubling the coverage compared to existing systematic studies. Our results reveal additional specificity determinants for a large number of factors for which a partial specificity was known, including a commonly observed A- or T-rich stretch that flanks the core motifs. Global analysis of the data revealed that homodimer orientation and spacing preferences, and base-stacking interactions, have a larger role in TF-DNA binding than previously appreciated. We further describe a binding model incorporating these features that is required to understand binding of TFs to DNA.

Mice lacking a Myc enhancer that includes human SNP rs6983267 are resistant to intestinal tumors.

Multiple cancer-associated single-nucleotide polymorphisms (SNPs) have been mapped to conserved sequences within a 500-kilobase region upstream of the MYC oncogene on human chromosome 8q24. These SNPs may affect cancer development through altered regulation of MYC expression, but this hypothesis has been difficult to confirm. We generated mice deficient in Myc-335, a putative MYC regulatory element that contains rs6983267, a SNP accounting for more human cancer-related morbidity than any other genetic variant or mutation. In Myc-335 null mice, Myc transcripts were expressed in the intestinal crypts in a pattern similar to that in wild-type mice but at modestly reduced levels. The mutant mice displayed no overt phenotype but were markedly resistant to intestinal tumorigenesis induced by the APCmin mutation. These results establish that a cancer-associated SNP identified in human genome-wide association studies has a functional effect in vivo.