Single-cell chromatin accessibility reveals principles of regulatory variation.

Cell-to-cell variation is a universal feature of life that affects a wide range of biological phenomena, from developmental plasticity to tumour heterogeneity. Although recent advances have improved our ability to document cellular phenotypic variation, the fundamental mechanisms that generate variability from identical DNA sequences remain elusive. Here we reveal the landscape and principles of mammalian DNA regulatory variation by developing a robust method for mapping the accessible genome of individual cells by assay for transposase-accessible chromatin using sequencing (ATAC-seq) integrated into a programmable microfluidics platform. Single-cell ATAC-seq (scATAC-seq) maps from hundreds of single cells in aggregate closely resemble accessibility profiles from tens of millions of cells and provide insights into cell-to-cell variation. Accessibility variance is systematically associated with specific trans-factors and cis-elements, and we discover combinations of trans-factors associated with either induction or suppression of cell-to-cell variability. We further identify sets of trans-factors associated with cell-type-specific accessibility variance across eight cell types. Targeted perturbations of cell cycle or transcription factor signalling evoke stimulus-specific changes in this observed variability. The pattern of accessibility variation in cis across the genome recapitulates chromosome compartments de novo, linking single-cell accessibility variation to three-dimensional genome organization. Single-cell analysis of DNA accessibility provides new insight into cellular variation of the 'regulome'.

Mice with an isoform-ablating Mecp2 exon 1 mutation recapitulate the neurologic deficits of Rett syndrome.

Mutations in MECP2 cause the neurodevelopmental disorder Rett syndrome (RTT OMIM 312750). Alternative inclusion of MECP2/Mecp2 exon 1 with exons 3 and 4 encodes MeCP2-e1 or MeCP2-e2 protein isoforms with unique amino termini. While most MECP2 mutations are located in exons 3 and 4 thus affecting both isoforms, MECP2 exon 1 mutations but not exon 2 mutations have been identified in RTT patients, suggesting that MeCP2-e1 deficiency is sufficient to cause RTT. As expected, genetic deletion of Mecp2 exons 3 and/or 4 recapitulates RTT-like neurologic defects in mice. However, Mecp2 exon 2 knockout mice have normal neurologic function. Here, a naturally occurring MECP2 exon 1 mutation is recapitulated in a mouse model by genetic engineering. A point mutation in the translational start codon of Mecp2 exon 1, transmitted through the germline, ablates MeCP2-e1 translation while preserving MeCP2-e2 production in mouse brain. The resulting MeCP2-e1 deficient mice developed forelimb stereotypy, hindlimb clasping, excessive grooming and hypo-activity prior to death between 7 and 31 weeks. MeCP2-e1 deficient mice also exhibited abnormal anxiety, sociability and ambulation. Despite MeCP2-e1 and MeCP2-e2 sharing, 96% amino acid identity, differences were identified. A fraction of phosphorylated MeCP2-e1 differed from the bulk of MeCP2 in subnuclear localization and co-factor interaction. Furthermore, MeCP2-e1 exhibited enhanced stability compared with MeCP2-e2 in neurons. Therefore, MeCP2-e1 deficient mice implicate MeCP2-e1 as the sole contributor to RTT with non-redundant functions.

R-loop formation at Snord116 mediates topotecan inhibition of Ube3a-antisense and allele-specific chromatin decondensation.

Prader-Willi syndrome (PWS) and Angelman syndrome (AS) are oppositely imprinted autism-spectrum disorders with known genetic bases, but complex epigenetic mechanisms underlie their pathogenesis. The PWS/AS locus on 15q11-q13 is regulated by an imprinting control region that is maternally methylated and silenced. The PWS imprinting control region is the promoter for a one megabase paternal transcript encoding the ubiquitous protein-coding Snrpn gene and multiple neuron-specific noncoding RNAs, including the PWS-related Snord116 repetitive locus of small nucleolar RNAs and host genes, and the antisense transcript to AS-causing ubiquitin ligase encoding Ube3a (Ube3a-ATS). Neuron-specific transcriptional progression through Ube3a-ATS correlates with paternal Ube3a silencing and chromatin decondensation. Interestingly, topoisomerase inhibitors, including topotecan, were recently identified in an unbiased drug screen for compounds that could reverse the silent paternal allele of Ube3a in neurons, but the mechanism of topotecan action on the PWS/AS locus is unknown. Here, we demonstrate that topotecan treatment stabilizes the formation of RNA:DNA hybrids (R loops) at G-skewed repeat elements within paternal Snord116, corresponding to increased chromatin decondensation and inhibition of Ube3a-ATS expression. Neural precursor cells from paternal Snord116 deletion mice exhibit increased Ube3a-ATS levels in differentiated neurons and show a reduced effect of topotecan compared with wild-type neurons. These results demonstrate that the AS candidate drug topotecan acts predominantly through stabilizing R loops and chromatin decondensation at the paternally expressed PWS Snord116 locus. Our study holds promise for targeted therapies to the Snord116 locus for both AS and PWS.

A PtdIns4,5P2-regulated nuclear poly(A) polymerase controls expression of select mRNAs.

Phosphoinositides are a family of lipid signalling molecules that regulate many cellular functions in eukaryotes. Phosphatidylinositol-4,5-bisphosphate (PtdIns4,5P2), the central component in the phosphoinositide signalling circuitry, is generated primarily by type I phosphatidylinositol 4-phosphate 5-kinases (PIPKIalpha, PIPKIbeta and PIPKIgamma). In addition to functions in the cytosol, phosphoinositides are present in the nucleus, where they modulate several functions; however, the mechanism by which they directly regulate nuclear functions remains unknown. PIPKIs regulate cellular functions through interactions with protein partners, often PtdIns4,5P2 effectors, that target PIPKIs to discrete subcellular compartments, resulting in the spatial and temporal generation of PtdIns4,5P2 required for the regulation of specific signalling pathways. Therefore, to determine roles for nuclear PtdIns4,5P2 we set out to identify proteins that interacted with the nuclear PIPK, PIPKIalpha. Here we show that PIPKIalpha co-localizes at nuclear speckles and interacts with a newly identified non-canonical poly(A) polymerase, which we have termed Star-PAP (nuclear speckle targeted PIPKIalpha regulated-poly(A) polymerase) and that the activity of Star-PAP can be specifically regulated by PtdIns4,5P2. Star-PAP and PIPKIalpha function together in a complex to control the expression of select mRNAs, including the transcript encoding the key cytoprotective enzyme haem oxygenase-1 (refs 8, 9) and other oxidative stress response genes by regulating the 3'-end formation of their mRNAs. Taken together, the data demonstrate a model by which phosphoinositide signalling works in tandem with complement pathways to regulate the activity of Star-PAP and the subsequent biosynthesis of its target mRNA. The results reveal a mechanism for the integration of nuclear phosphoinositide signals and a method for regulating gene expression.

MeCP2 is required for global heterochromatic and nucleolar changes during activity-dependent neuronal maturation.

Mutations in MECP2, encoding methyl CpG binding protein 2, cause the neurodevelopmental disorder Rett syndrome. MeCP2 is an abundant nuclear protein that binds to chromatin and modulates transcription in response to neuronal activity. Prior studies of MeCP2 function have focused on specific gene targets of MeCP2, but a more global role for MeCP2 in neuronal nuclear maturation has remained unexplored. MeCP2 levels increase during postnatal brain development, coinciding with dynamic changes in neuronal chromatin architecture, particularly detectable as changes in size, number, and location of nucleoli and perinucleolar heterochromatic chromocenters. To determine a potential role for MeCP2 in neuronal chromatin maturational changes, we measured nucleoli and chromocenters in developing wild-type and Mecp2-deficient mouse cortical sections, as well as mouse primary cortical neurons and a human neuronal cell line following induced maturation. Mecp2-deficient mouse neurons exhibited significant differences in nucleolar and chromocenter number and size, as more abundant, smaller nucleoli in brain and primary neurons compared to wild-type, consistent with delayed neuronal nuclear maturation in the absence of MeCP2. Primary neurons increased chromocenter size following depolarization in wild-type, but not Mecp2-deficient cultures. Wild-type MECP2e1 over-expression in human SH-SY5Y cells was sufficient to induce significantly larger nucleoli, but not a T158M mutation of the methyl-binding domain. These results suggest that, in addition to the established role of MeCP2 in transcriptional regulation of specific target genes, the global chromatin-binding function of MeCP2 is essential for activity-dependent global chromatin dynamics during postnatal neuronal maturation.

The role of MeCP2 in brain development and neurodevelopmental disorders.

Methyl CpG binding protein-2 (MeCP2) is an essential epigenetic regulator in human brain development. Rett syndrome, the primary disorder caused by mutations in the X-linked MECP2 gene, is characterized by a period of cognitive decline and development of hand stereotypies and seizures following an apparently normal early infancy. In addition, MECP2 mutations and duplications are observed in a spectrum of neurodevelopmental disorders, including severe neonatal encephalopathy, X-linked mental retardation, and autism, implicating MeCP2 as an essential regulator of postnatal brain development. In this review, we compare the mutation types and inheritance patterns of the human disorders associated with MECP2. In addition, we summarize the current understanding of MeCP2 as a central epigenetic regulator of activity-dependent synaptic maturation. As MeCP2 occupies a central role in the pathogenesis of multiple neurodevelopmental disorders, continued investigation into MeCP2 function and regulatory pathways may show promise for developing broad-spectrum therapies.

A novel hypomorphic MECP2 point mutation is associated with a neuropsychiatric phenotype.

The MECP2 gene on Xq28 encodes a transcriptional repressor, which binds to and modulates expression of active genes. Mutations in MECP2 cause classic or preserved speech variant Rett syndrome and intellectual disability in females and early demise or marked neurodevelopmental handicap in males. The consequences of a hypomorphic Mecp2 allele were recently investigated in a mouse model, which developed obesity, motor, social, learning, and behavioral deficits, predicting a human neurobehavioral syndrome. Here, we describe mutation analysis of a nondysmorphic female proband and her father who presented with primarily neuropsychiatric manifestations and obesity with relative sparing of intelligence, language, growth, and gross motor skills. We identified and characterized a novel missense mutation (c.454C>G; p.P152A) in the critical methyl-binding domain of MeCP2 that disrupts MeCP2 functional activity. We show that a gradient of impairment is present when the p.P152A mutation is compared with an allelic p.P152R mutation, which causes classic Rett syndrome and another Rett syndrome-causing mutation, such that protein-heterochromatin binding observed by immunofluorescence and immunoblotting is wild-type > P152A > P152R > T158 M, consistent with the severity of the observed phenotype. Our findings provide evidence for very mild phenotypes in humans associated with partial reduction of MeCP2 function arising from subtle variation in MECP2.

Multimodal Single-Cell Analysis Reveals Physiological Maturation in the Developing Human Neocortex.

In the developing human neocortex, progenitor cells generate diverse cell types prenatally. Progenitor cells and newborn neurons respond to signaling cues, including neurotransmitters. While single-cell RNA sequencing has revealed cellular diversity, physiological heterogeneity has yet to be mapped onto these developing and diverse cell types. By combining measurements of intracellular Ca2+ elevations in response to neurotransmitter receptor agonists and RNA sequencing of the same single cells, we show that Ca2+ responses are cell-type-specific and change dynamically with lineage progression. Physiological response properties predict molecular cell identity and additionally reveal diversity not captured by single-cell transcriptomics. We find that the serotonin receptor HTR2A selectively activates radial glia cells in the developing human, but not mouse, neocortex, and inhibiting HTR2A receptors in human radial glia disrupts the radial glial scaffold. We show highly specific neurotransmitter signaling during neurogenesis in the developing human neocortex and highlight evolutionarily divergent mechanisms of physiological signaling.

Stress-ING out: phosphoinositides mediate the cellular stress response.

Phosphoinositides regulate numerous cellular processes required for growth, proliferation, and motility. Whereas phosphoinositide signal transduction pathways within the cytosol have been well characterized, nuclear signaling pathways remain poorly understood. Accumulating experimental data have now started to uncover critical functions for nuclear phosphoinositides. In particular, phosphoinositides modulate the activity of the tumor suppressor protein ING2 in response to extracellular stress. These findings highlight a previously uncharacterized function for phosphoinositides and implicate their metabolism in signaling pathways critical for cell survival.

Corrigendum: Fluidic Logic Used in a Systems Approach to Enable Integrated Single-Cell Functional Analysis.

[This corrects the article on p. 70 in vol. 4, PMID: 27709111.].

Fluidic Logic Used in a Systems Approach to Enable Integrated Single-Cell Functional Analysis.

The study of single cells has evolved over the past several years to include expression and genomic analysis of an increasing number of single cells. Several studies have demonstrated wide spread variation and heterogeneity within cell populations of similar phenotype. While the characterization of these populations will likely set the foundation for our understanding of genomic- and expression-based diversity, it will not be able to link the functional differences of a single cell to its underlying genomic structure and activity. Currently, it is difficult to perturb single cells in a controlled environment, monitor and measure the response due to perturbation, and link these response measurements to downstream genomic and transcriptomic analysis. In order to address this challenge, we developed a platform to integrate and miniaturize many of the experimental steps required to study single-cell function. The heart of this platform is an elastomer-based integrated fluidic circuit that uses fluidic logic to select and sequester specific single cells based on a phenotypic trait for downstream experimentation. Experiments with sequestered cells that have been performed include on-chip culture, exposure to various stimulants, and post-exposure image-based response analysis, followed by preparation of the mRNA transcriptome for massively parallel sequencing analysis. The flexible system embodies experimental design and execution that enable routine functional studies of single cells.

Single-Cell Genetic Analysis Using Automated Microfluidics to Resolve Somatic Mosaicism.

Somatic mosaicism occurs throughout normal development and contributes to numerous disease etiologies, including tumorigenesis and neurological disorders. Intratumor genetic heterogeneity is inherent to many cancers, creating challenges for effective treatments. Unfortunately, analysis of bulk DNA masks subclonal phylogenetic architectures created by the acquisition and distribution of somatic mutations amongst cells. As a result, single-cell genetic analysis is becoming recognized as vital for accurately characterizing cancers. Despite this, methods for single-cell genetics are lacking. Here we present an automated microfluidic workflow enabling efficient cell capture, lysis, and whole genome amplification (WGA). We find that ~90% of the genome is accessible in single cells with improved uniformity relative to current single-cell WGA methods. Allelic dropout (ADO) rates were limited to 13.75% and variant false discovery rates (SNV FDR) were 4.11x10(-6), on average. Application to ER-/PR-/HER2+ breast cancer cells and matched normal controls identified novel mutations that arose in a subpopulation of cells and effectively resolved the segregation of known cancer-related mutations with single-cell resolution. Finally, we demonstrate effective cell classification using mutation profiles with 10X average exome coverage depth per cell. Our data demonstrate an efficient automated microfluidic platform for single-cell WGA that enables the resolution of somatic mutation patterns in single cells.

Phosphorylation of distinct sites in MeCP2 modifies cofactor associations and the dynamics of transcriptional regulation.

Mutations in the gene encoding methyl-CpG-binding protein 2 (MeCP2) lead to disrupted neuronal function and can cause the neurodevelopmental disorder Rett syndrome. MeCP2 is a transcriptional regulator that binds to methylated DNA and is most abundant in neuronal nuclei. The mechanisms by which MeCP2 regulates gene expression remain ambiguous, as it has been reported to function as a transcriptional silencer or activator and to execute these activities through both gene-specific and genome-wide mechanisms. We hypothesized that posttranslational modifications of MeCP2 may be important for reconciling these apparently contradictory functions. Our results demonstrate that MeCP2 contains multiple posttranslational modifications, including phosphorylation, acetylation, and ubiquitylation. Phosphorylation of MeCP2 at S229 or S80 influenced selective in vivo interactions with the chromatin factors HP1 and SMC3 and the cofactors Sin3A and YB-1. pS229 MeCP2 was specifically enriched at the RET promoter, and phosphorylation of MeCP2 was necessary for differentiation-induced activation and repression of the MeCP2 target genes RET and EGR2. These results demonstrate that phosphorylation is one of several factors that are important for interpreting the complexities of MeCP2 transcriptional modulation.

CKIalpha is associated with and phosphorylates star-PAP and is also required for expression of select star-PAP target messenger RNAs.

We have recently identified Star-PAP, a nuclear poly(A) polymerase that associates with phosphatidylinositol-4-phosphate 5-kinase Ialpha (PIPKIalpha) and is required for the expression of a specific subset of mRNAs. Star-PAP activity is directly modulated by the PIPKIalpha product phosphatidylinositol 4,5-bisphosphate (PI-4,5-P(2)), linking nuclear phosphoinositide signaling to gene expression. Here, we show that PI-4,5-P(2)-dependent protein kinase activity is also a part of the Star-PAP protein complex. We identify the PI-4,5-P(2)-sensitive casein kinase Ialpha (CKIalpha) as a protein kinase responsible for this activity and further show that CKIalpha is capable of directly phosphorylating Star-PAP. Both CKIalpha and PIPKIalpha are required for the synthesis of some but not all Star-PAP target mRNA, and like Star-PAP, CKIalpha is associated with these messages in vivo. Taken together, these data indicate that CKIalpha, PIPKIalpha, and Star-PAP function together to modulate the production of specific Star-PAP messages. The Star-PAP complex therefore represents a location where multiple signaling pathways converge to regulate the expression of specific mRNAs.

Nuclear phosphoinositide kinases and inositol phospholipids.

The presence of inositol phospholipids in the nuclei of mammalian cells has by now been well established, as has the presence of the enzymes responsible for their metabolism. However, our understanding of the role of these nuclear phosphoinositides in regulating cellular events has lagged far behind that for its cytosolic counterpart. It is clear, though, that the nuclear phosphoinositide pool is independent of the cytosolic pool and is, therefore, likely to be regulating a unique set of cellular events. As with its cytosolic phosphoinositides, many nuclear phosphoinositides and their metabolic enzymes are located at distinct sub-cellular structures. This arrangement spatially limits the production and activity of inositol phospholipids and is believed to be a major mechanism for regulating their function. Here, we will introduce the components of nuclear inositol phospholipid signal transduction and discuss how their spatial arrangement may dictate which nuclear functions they are modulating.