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.

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.

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.

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.

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.

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 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.

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.

Genome-wide analysis of ETS-family DNA-binding in vitro and in vivo.

Members of the large ETS family of transcription factors (TFs) have highly similar DNA-binding domains (DBDs)-yet they have diverse functions and activities in physiology and oncogenesis. Some differences in DNA-binding preferences within this family have been described, but they have not been analysed systematically, and their contributions to targeting remain largely uncharacterized. We report here the DNA-binding profiles for all human and mouse ETS factors, which we generated using two different methods: a high-throughput microwell-based TF DNA-binding specificity assay, and protein-binding microarrays (PBMs). Both approaches reveal that the ETS-binding profiles cluster into four distinct classes, and that all ETS factors linked to cancer, ERG, ETV1, ETV4 and FLI1, fall into just one of these classes. We identify amino-acid residues that are critical for the differences in specificity between all the classes, and confirm the specificities in vivo using chromatin immunoprecipitation followed by sequencing (ChIP-seq) for a member of each class. The results indicate that even relatively small differences in in vitro binding specificity of a TF contribute to site selectivity in vivo.

Counting absolute numbers of molecules using unique molecular identifiers.

Counting individual RNA or DNA molecules is difficult because they are hard to copy quantitatively for detection. To overcome this limitation, we applied unique molecular identifiers (UMIs), which make each molecule in a population distinct, to genome-scale human karyotyping and mRNA sequencing in Drosophila melanogaster. Use of this method can improve accuracy of almost any next-generation sequencing method, including chromatin immunoprecipitation-sequencing, genome assembly, diagnostics and manufacturing-process control and monitoring.

Protein kinase Cα (PKCα) regulates p53 localization and melanoma cell survival downstream of integrin αv in three-dimensional collagen and in vivo.

Protein kinase C α (PKCα) is overexpressed in numerous types of cancer. Importantly, PKCα has been linked to metastasis of malignant melanoma in patients. However, it has been unclear how PKCα may be regulated and how it exerts its role in melanoma. Here, we identified a role for PKCα in melanoma cell survival in a three-dimensional collagen model mimicking the in vivo pathophysiology of the dermis. A pathway was identified that involved integrin αv-mediated up-regulation of PKCα and PKCα-dependent regulation of p53 localization, which was connected to melanoma cell survival. Melanoma survival and growth in three-dimensional microenvironments requires the expression of integrin αv, which acts to suppress p53 activity. Interestingly, microarray analysis revealed that PKCα was up-regulated by integrin αv in a three-dimensional microenvironment-dependent manner. Integrin αv was observed to promote a relocalization of endogenous p53 from the nucleus to the cytoplasm upon growth in three-dimensional collagen as well as in vivo, whereas stable knockdown of PKCα inhibited the integrin αv-mediated relocalization of p53. Importantly, knockdown of PKCα also promoted apoptosis in three-dimensional collagen and in vivo, resulting in reduced tumor growth. This indicates that PKCα constitutes a crucial component of the integrin αv-mediated pathway(s) that promote p53 relocalization and melanoma survival.

Multiplexed massively parallel SELEX for characterization of human transcription factor binding specificities.

The genetic code-the binding specificity of all transfer-RNAs--defines how protein primary structure is determined by DNA sequence. DNA also dictates when and where proteins are expressed, and this information is encoded in a pattern of specific sequence motifs that are recognized by transcription factors. However, the DNA-binding specificity is only known for a small fraction of the approximately 1400 human transcription factors (TFs). We describe here a high-throughput method for analyzing transcription factor binding specificity that is based on systematic evolution of ligands by exponential enrichment (SELEX) and massively parallel sequencing. The method is optimized for analysis of large numbers of TFs in parallel through the use of affinity-tagged proteins, barcoded selection oligonucleotides, and multiplexed sequencing. Data are analyzed by a new bioinformatic platform that uses the hundreds of thousands of sequencing reads obtained to control the quality of the experiments and to generate binding motifs for the TFs. The described technology allows higher throughput and identification of much longer binding profiles than current microarray-based methods. In addition, as our method is based on proteins expressed in mammalian cells, it can also be used to characterize DNA-binding preferences of full-length proteins or proteins requiring post-translational modifications. We validate the method by determining binding specificities of 14 different classes of TFs and by confirming the specificities for NFATC1 and RFX3 using ChIP-seq. Our results reveal unexpected dimeric modes of binding for several factors that were thought to preferentially bind DNA as monomers.

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.

Inhibition of glycolytic enzymes mediated by pharmacologically activated p53: targeting Warburg effect to fight cancer.

Unique sensitivity of tumor cells to the inhibition of glycolysis is a good target for anticancer therapy. Here, we demonstrate that the pharmacologically activated tumor suppressor p53 mediates the inhibition of glycolytic enzymes in cancer cells in vitro and in vivo. We showed that p53 binds to the promoters of metabolic genes and represses their expression, including glucose transporters SLC2A12 (GLUT12) and SLC2A1 (GLUT1). Furthermore, p53-mediated repression of transcription factors c-Myc and HIF1α, key drivers of ATP-generating pathways in tumors, contributed to ATP production block. Inhibition of c-Myc by p53 mediated the ablation of several glycolytic genes in normoxia, whereas in hypoxia down-regulation of HIF1α contributed to this effect. We identified Sp1 as a transcription cofactor cooperating with p53 in the ablation of metabolic genes. Using different approaches, we demonstrated that glycolysis block contributes to the robust induction of apoptosis by p53 in cancer cells. Taken together, our data suggest that tumor-specific reinstatement of p53 function targets the "Achilles heel" of cancer cells (i.e. their dependence on glycolysis), which could contribute to the tumor-selective killing of cancer cells by pharmacologically activated p53.

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.

Small molecule RITA binds to p53, blocks p53-HDM-2 interaction and activates p53 function in tumors.

In tumors that retain wild-type p53, its tumor-suppressor function is often impaired as a result of the deregulation of HDM-2, which binds to p53 and targets it for proteasomal degradation. We have screened a chemical library and identified a small molecule named RITA (reactivation of p53 and induction of tumor cell apoptosis), which bound to p53 and induced its accumulation in tumor cells. RITA prevented p53-HDM-2 interaction in vitro and in vivo and affected p53 interaction with several negative regulators. RITA induced expression of p53 target genes and massive apoptosis in various tumor cells lines expressing wild-type p53. RITA suppressed the growth of human fibroblasts and lymphoblasts only upon oncogene expression and showed substantial p53-dependent antitumor effect in vivo. RITA may serve as a lead compound for the development of an anticancer drug that targets tumors with wild-type p53.