Activation of the imprinted Prader-Willi Syndrome locus by CRISPR-based epigenome editing.

Epigenome editing with DNA-targeting technologies such as CRISPR-dCas9 can be used to dissect gene regulatory mechanisms and potentially treat associated disorders. For example, Prader-Willi Syndrome (PWS) is caused by loss of paternally expressed imprinted genes on chromosome 15q11.2-q13.3, although the maternal allele is intact but epigenetically silenced. Using CRISPR repression and activation screens in human induced pluripotent stem cells (iPSCs), we identified genomic elements that control expression of the PWS gene SNRPN from the paternal and maternal chromosomes. We showed that either targeted transcriptional activation or DNA demethylation can activate the silenced maternal SNRPN and downstream PWS transcripts. However, these two approaches function at unique regions, preferentially activating different transcript variants and involving distinct epigenetic reprogramming mechanisms. Remarkably, transient expression of the targeted demethylase leads to stable, long-term maternal SNRPN expression in PWS iPSCs. This work uncovers targeted epigenetic manipulations to reprogram a disease-associated imprinted locus and suggests possible therapeutic interventions.

Parallelized engineering of mutational models using piggyBac transposon delivery of CRISPR libraries.

New technologies and large-cohort studies have enabled novel variant discovery and association at unprecedented scale, yet functional characterization of these variants remains paramount to deciphering disease mechanisms. Approaches that facilitate parallelized genome editing of cells of interest or induced pluripotent stem cells (iPSCs) have become critical tools toward this goal. Here, we developed an approach that incorporates libraries of CRISPR-Cas9 guide RNAs (gRNAs) together with inducible Cas9 into a piggyBac (PB) transposon system to engineer dozens to hundreds of genomic variants in parallel against isogenic cellular backgrounds. This method empowers loss-of-function (LoF) studies through the introduction of insertions or deletions (indels) and copy-number variants (CNVs), though generating specific nucleotide changes is possible with prime editing. The ability to rapidly establish high-quality mutational models at scale will facilitate the development of isogenic cellular collections and catalyze comparative functional genomic studies investigating the roles of hundreds of genes and mutations in development and disease.

Excitatory Dysfunction Drives Network and Calcium Handling Deficits in 16p11.2 Duplication Schizophrenia Induced Pluripotent Stem Cell-Derived Neurons.

BACKGROUND: Schizophrenia (SCZ) is a debilitating psychiatric disorder with a large genetic contribution; however, its neurodevelopmental substrates remain largely unknown. Modeling pathogenic processes in SCZ using human induced pluripotent stem cell-derived neurons (iNs) has emerged as a promising strategy. Copy number variants confer high genetic risk for SCZ, with duplication of the 16p11.2 locus increasing the risk 14.5-fold. METHODS: To dissect the contribution of induced excitatory neurons (iENs) versus GABAergic (gamma-aminobutyric acidergic) neurons (iGNs) to SCZ pathophysiology, we induced iNs from CRISPR (clustered regularly interspaced short palindromic repeats)-Cas9 isogenic and SCZ patient-derived induced pluripotent stem cells and analyzed SCZ-related phenotypes in iEN monocultures and iEN/iGN cocultures. RESULTS: In iEN/iGN cocultures, neuronal firing and synchrony were reduced at later, but not earlier, stages of in vitro development. These were fully recapitulated in iEN monocultures, indicating a primary role for iENs. Moreover, isogenic iENs showed reduced dendrite length and deficits in calcium handling. iENs from 16p11.2 duplication-carrying patients with SCZ displayed overlapping deficits in network synchrony, dendrite outgrowth, and calcium handling. Transcriptomic analysis of both iEN cohorts revealed molecular markers of disease related to the glutamatergic synapse, neuroarchitecture, and calcium regulation. CONCLUSIONS: Our results indicate the presence of 16p11.2 duplication-dependent alterations in SCZ patient-derived iENs. Transcriptomics and cellular phenotyping reveal overlap between isogenic and patient-derived iENs, suggesting a central role of glutamatergic, morphological, and calcium dysregulation in 16p11.2 duplication-mediated pathogenesis. Moreover, excitatory dysfunction during early neurodevelopment is implicated as the basis of SCZ pathogenesis in 16p11.2 duplication carriers. Our results support network synchrony and calcium handling as outcomes directly linked to this genetic risk variant.

Transcriptional and functional consequences of alterations to MEF2C and its topological organization in neuronal models.

Point mutations and structural variants that directly disrupt the coding sequence of MEF2C have been associated with a spectrum of neurodevelopmental disorders (NDDs). However, the impact of MEF2C haploinsufficiency on neurodevelopmental pathways and synaptic processes is not well understood, nor are the complex mechanisms that govern its regulation. To explore the functional changes associated with structural variants that alter MEF2C expression and/or regulation, we generated an allelic series of 204 isogenic human induced pluripotent stem cell (hiPSC)-derived neural stem cells and glutamatergic induced neurons. These neuronal models harbored CRISPR-engineered mutations that involved direct deletion of MEF2C or deletion of the boundary points for topologically associating domains (TADs) and chromatin loops encompassing MEF2C. Systematic profiling of mutation-specific alterations, contrasted to unedited controls that were exposed to the same guide RNAs for each edit, revealed that deletion of MEF2C caused differential expression of genes associated with neurodevelopmental pathways and synaptic function. We also discovered significant reduction in synaptic activity measured by multielectrode arrays (MEAs) in neuronal cells. By contrast, we observed robust buffering against MEF2C regulatory disruption following deletion of a distal 5q14.3 TAD and loop boundary, whereas homozygous loss of a proximal loop boundary resulted in down-regulation of MEF2C expression and reduced electrophysiological activity on MEA that was comparable to direct gene disruption. Collectively, these studies highlight the considerable functional impact of MEF2C deletion in neuronal cells and systematically characterize the complex interactions that challenge a priori predictions of regulatory consequences from structural variants that disrupt three-dimensional genome organization.

Statistical and functional convergence of common and rare genetic influences on autism at chromosome 16p.

The canonical paradigm for converting genetic association to mechanism involves iteratively mapping individual associations to the proximal genes through which they act. In contrast, in the present study we demonstrate the feasibility of extracting biological insights from a very large region of the genome and leverage this strategy to study the genetic influences on autism. Using a new statistical approach, we identified the 33-Mb p-arm of chromosome 16 (16p) as harboring the greatest excess of autism's common polygenic influences. The region also includes the mechanistically cryptic and autism-associated 16p11.2 copy number variant. Analysis of RNA-sequencing data revealed that both the common polygenic influences within 16p and the 16p11.2 deletion were associated with decreased average gene expression across 16p. The transcriptional effects of the rare deletion and diffuse common variation were correlated at the level of individual genes and analysis of Hi-C data revealed patterns of chromatin contact that may explain this transcriptional convergence. These results reflect a new approach for extracting biological insight from genetic association data and suggest convergence of common and rare genetic influences on autism at 16p.

Tissue- and cell-type-specific molecular and functional signatures of 16p11.2 reciprocal genomic disorder across mouse brain and human neuronal models.

Chromosome 16p11.2 reciprocal genomic disorder, resulting from recurrent copy-number variants (CNVs), involves intellectual disability, autism spectrum disorder (ASD), and schizophrenia, but the responsible mechanisms are not known. To systemically dissect molecular effects, we performed transcriptome profiling of 350 libraries from six tissues (cortex, cerebellum, striatum, liver, brown fat, and white fat) in mouse models harboring CNVs of the syntenic 7qF3 region, as well as cellular, transcriptional, and single-cell analyses in 54 isogenic neural stem cell, induced neuron, and cerebral organoid models of CRISPR-engineered 16p11.2 CNVs. Transcriptome-wide differentially expressed genes were largely tissue-, cell-type-, and dosage-specific, although more effects were shared between deletion and duplication and across tissue than expected by chance. The broadest effects were observed in the cerebellum (2,163 differentially expressed genes), and the greatest enrichments were associated with synaptic pathways in mouse cerebellum and human induced neurons. Pathway and co-expression analyses identified energy and RNA metabolism as shared processes and enrichment for ASD-associated, loss-of-function constraint, and fragile X messenger ribonucleoprotein target gene sets. Intriguingly, reciprocal 16p11.2 dosage changes resulted in consistent decrements in neurite and electrophysiological features, and single-cell profiling of organoids showed reciprocal alterations to the proportions of excitatory and inhibitory GABAergic neurons. Changes both in neuronal ratios and in gene expression in our organoid analyses point most directly to calretinin GABAergic inhibitory neurons and the excitatory/inhibitory balance as targets of disruption that might contribute to changes in neurodevelopmental and cognitive function in 16p11.2 carriers. Collectively, our data indicate the genomic disorder involves disruption of multiple contributing biological processes and that this disruption has relative impacts that are context specific.

Orgo-Seq integrates single-cell and bulk transcriptomic data to identify cell type specific-driver genes associated with autism spectrum disorder.

Cerebral organoids can be used to gain insights into cell type specific processes perturbed by genetic variants associated with neuropsychiatric disorders. However, robust and scalable phenotyping of organoids remains challenging. Here, we perform RNA sequencing on 71 samples comprising 1,420 cerebral organoids from 25 donors, and describe a framework (Orgo-Seq) to integrate bulk RNA and single-cell RNA sequence data. We apply Orgo-Seq to 16p11.2 deletions and 15q11-13 duplications, two loci associated with autism spectrum disorder, to identify immature neurons and intermediate progenitor cells as critical cell types for 16p11.2 deletions. We further applied Orgo-Seq to identify cell type-specific driver genes. Our work presents a quantitative phenotyping framework to integrate multi-transcriptomic datasets for the identification of cell types and cell type-specific co-expressed driver genes associated with neuropsychiatric disorders.

16p11.2 deletion is associated with hyperactivation of human iPSC-derived dopaminergic neuron networks and is rescued by RHOA inhibition in vitro.

Reciprocal copy number variations (CNVs) of 16p11.2 are associated with a wide spectrum of neuropsychiatric and neurodevelopmental disorders. Here, we use human induced pluripotent stem cells (iPSCs)-derived dopaminergic (DA) neurons carrying CNVs of 16p11.2 duplication (16pdup) and 16p11.2 deletion (16pdel), engineered using CRISPR-Cas9. We show that 16pdel iPSC-derived DA neurons have increased soma size and synaptic marker expression compared to isogenic control lines, while 16pdup iPSC-derived DA neurons show deficits in neuronal differentiation and reduced synaptic marker expression. The 16pdel iPSC-derived DA neurons have impaired neurophysiological properties. The 16pdel iPSC-derived DA neuronal networks are hyperactive and have increased bursting in culture compared to controls. We also show that the expression of RHOA is increased in the 16pdel iPSC-derived DA neurons and that treatment with a specific RHOA-inhibitor, Rhosin, rescues the network activity of the 16pdel iPSC-derived DA neurons. Our data suggest that 16p11.2 deletion-associated iPSC-derived DA neuron hyperactivation can be rescued by RHOA inhibition.

Transcriptional consequences of MBD5 disruption in mouse brain and CRISPR-derived neurons.

BACKGROUND: MBD5, encoding the methyl-CpG-binding domain 5 protein, has been proposed as a necessary and sufficient driver of the 2q23.1 microdeletion syndrome. De novo missense and protein-truncating variants from exome sequencing studies have directly implicated MBD5 in the etiology of autism spectrum disorder (ASD) and related neurodevelopmental disorders (NDDs). However, little is known concerning the specific function(s) of MBD5. METHODS: To gain insight into the complex interactions associated with alteration of MBD5 in individuals with ASD and related NDDs, we explored the transcriptional landscape of MBD5 haploinsufficiency across multiple mouse brain regions of a heterozygous hypomorphic Mbd5+/GT mouse model, and compared these results to CRISPR-mediated mutations of MBD5 in human iPSC-derived neuronal models. RESULTS: Gene expression analyses across three brain regions from Mbd5+/GT mice showed subtle transcriptional changes, with cortex displaying the most widespread changes following Mbd5 reduction, indicating context-dependent effects. Comparison with MBD5 reduction in human neuronal cells reinforced the context-dependence of gene expression changes due to MBD5 deficiency. Gene co-expression network analyses revealed gene clusters that were associated with reduced MBD5 expression and enriched for terms related to ciliary function. LIMITATIONS: These analyses included a limited number of mouse brain regions and neuronal models, and the effects of the gene knockdown are subtle. As such, these results will not reflect the full extent of MBD5 disruption across human brain regions during early neurodevelopment in ASD, or capture the diverse spectrum of cell-type-specific changes associated with MBD5 alterations. CONCLUSIONS: Our study points to modest and context-dependent transcriptional consequences of Mbd5 disruption in the brain. It also suggests a possible link between MBD5 and perturbations in ciliary function, which is an established pathogenic mechanism in developmental disorders and syndromes.

Restoring Wnt6 signaling ameliorates behavioral deficits in MeCP2 T158A mouse model of Rett syndrome.

The methyl-CpG-binding protein 2 gene, MECP2, is an X chromosome-linked gene encoding the MeCP2 protein, and mutations of MECP2 cause Rett syndrome (RTT). Previous study has shown that re-expression of SUMO-modified MeCP2 in Mecp2-null neurons rescues synaptic and behavioral deficits in Mecp2 conditional knockout mice, whereas about 12-fold decrease in Wnt6 mRNA level was found in MeCP2K412R sumo-mutant mice. Here, we examined the role of Wnt6 in MeCP2 T158A mouse model of RTT. Results show that lentiviral delivery of Wnt6 to the amygdala ameliorates locomotor impairment and social behavioral deficits in these animals. MeCP2 T158A mice show decreased level of GSK-3ß phosphorylation and increased level of ß-catenin phosphorylation. They also show reduced level of MeCP2 SUMOylation. These alterations were also restored by lenti-Wnt6 transduction. Further, both BDNF and IGF-1 expressions are decreased in MeCP2 T158A mice. Overexpression of Wnt6 increases Bdnf and Igf-1 promoter activity in HEK293T cells in a dose-dependent manner. Lenti-Wnt6 transduction to the amygdala similarly increases the mRNA level and protein expression of BDNF and IGF-1 in MeCP2 T158A mice. Moreover, environmental enrichment (EE) similarly ameliorates the locomotor and social behavioral deficits in MeCP2 T158A mice. One of the mechanisms underlying EE is mediated through enhanced MeCP2 SUMOylation and increased Wnt6 expression in these animals by EE.

Biallelic mutation of FBXL7 suggests a novel form of Hennekam syndrome.

Hennekam lymphangiectasia-lymphedema syndrome is an autosomal recessive disorder characterized by congenital lymphedema, intestinal lymphangiectasia, facial dysmorphism, and variable intellectual disability. Known disease genes include CCBE1, FAT4, and ADAMTS3. In a patient with clinically diagnosed Hennekam syndrome but without mutations or copy-number changes in the three known disease genes, we identified a homozygous single-exon deletion affecting FBXL7. Specifically, exon 3, which encodes the F-box domain and several leucine-rich repeats of FBXL7, is eliminated. Our analyses of databases representing >100,000 control individuals failed to identify biallelic loss-of-function variants in FBXL7. Published studies in Drosophila indicate Fbxl7 interacts with Fat, of which human FAT4 is an ortholog, and mutation of either gene yields similar morphological consequences. These data suggest that FBXL7 may be the fourth gene for Hennekam syndrome, acting via a shared pathway with FAT4.

MeCP2 SUMOylation rescues Mecp2-mutant-induced behavioural deficits in a mouse model of Rett syndrome.

The methyl-CpG-binding protein 2 (MeCP2) gene, MECP2, is an X-linked gene encoding the MeCP2 protein, and mutations of MECP2 cause Rett syndrome (RTT). However, the molecular mechanism of MECP2-mutation-caused RTT is less known. Here we find that MeCP2 could be SUMO-modified by the E3 ligase PIAS1 at Lys-412. MeCP2 phosphorylation (at Ser-421 and Thr-308) facilitates MeCP2 SUMOylation, and MeCP2 SUMOylation is induced by NMDA, IGF-1 and CRF in the rat brain. MeCP2 SUMOylation releases CREB from the repressor complex and enhances Bdnf mRNA expression. Several MECP2 mutations identified in RTT patients show decreased MeCP2 SUMOylation. Re-expression of wild-type MeCP2 or SUMO-modified MeCP2 in Mecp2-null neurons rescues the deficits of social interaction, fear memory and LTP observed in Mecp2 conditional knockout (cKO) mice. These results together reveal an important role of MeCP2 SUMOylation in social interaction, memory and synaptic plasticity, and that abnormal MeCP2 SUMOylation is implicated in RTT.

Engineering microdeletions and microduplications by targeting segmental duplications with CRISPR.

Recurrent, reciprocal genomic disorders resulting from non-allelic homologous recombination (NAHR) between near-identical segmental duplications (SDs) are a major cause of human disease, often producing phenotypically distinct syndromes. The genomic architecture of flanking SDs presents a challenge for modeling these syndromes; however, the capability to efficiently generate reciprocal copy number variants (CNVs) that mimic NAHR would represent a valuable modeling tool. We describe here a CRISPR/Cas9 genome engineering method, single-guide CRISPR/Cas targeting of repetitive elements (SCORE), to model reciprocal genomic disorders and demonstrate its capabilities by generating reciprocal CNVs of 16p11.2 and 15q13.3, including alteration of one copy-equivalent of the SDs that mediate NAHR in vivo. The method is reproducible, and RNA sequencing reliably clusters transcriptional signatures from human subjects with in vivo CNVs and their corresponding in vitro models. This new approach will provide broad applicability for the study of genomic disorders and, with further development, may also permit efficient correction of these defects.

CREB SUMOylation by the E3 ligase PIAS1 enhances spatial memory.

cAMP-responsive element binding protein (CREB) phosphorylation and signaling plays an important role in long-term memory formation, but other posttranslational modifications of CREB are less known. Here, we found that CREB1Δ, the short isoform of CREB, could be sumoylated by the small ubiquitin-like modifier (SUMO) E3 ligase protein inhibitor of activated STAT1 (PIAS1) at Lys271 and Lys290 and PIAS1 SUMOylation of CREB1Δ increased the expression level of CREB1Δ. CREB1Δ could also be sumoylated by other PIAS family proteins, but not by the E3 ligases RanBP2 and Pc2 or by the E2 ligase Ubc9. Furthermore, water maze training increased the level of endogenous CREB SUMOylation in rat CA1 neurons determined by in vitro SUMOylation assay, but this effect was not observed in other brain areas. Moreover, transduction of Lenti-CREBWT to rat CA1 area facilitated, whereas transduction of Lenti-CREB double sumo-mutant (CREBK271RK290R) impaired, spatial learning and memory performance. Transduction of Lenti-CREBWT-SUMO1 fusion vector to rat CA1 area showed a more significant effect in enhancing spatial learning and memory and CREB SUMOylation. Lenti-CREBWT transduction increased, whereas Lenti-CREBK271RK290R transduction decreased, CREB DNA binding to the brain-derived neurotrophic factor (bdnf) promoter and decreased bdnf mRNA expression. Knock-down of PIAS1 expression in CA1 area by PIAS1 siRNA transfection impaired spatial learning and memory and decreased endogenous CREB SUMOylation. In addition, CREB SUMOylation was CREB phosphorylation dependent and lasted longer. Therefore, CREB phosphorylation may be responsible for signal transduction during the early phase of long-term memory formation, whereas CREB SUMOylation sustains long-term memory.

Novel role and mechanism of protein inhibitor of activated STAT1 in spatial learning.

By using differential display PCR, we have previously identified 98 cDNA fragments from rat dorsal hippocampus, which are expressed differentially between the fast learners and slow learners from water-maze learning task. One cDNA fragment, which showed a higher expression level in fast learners, encodes the rat protein inhibitor of activated STAT1 (pias1) gene. Spatial training induced a significant increase in PIAS1 expression in rat hippocampus. Transient transfection of the wild-type (WT) PIAS1 plasmid to CA1 neurons facilitated, whereas transfection of PIAS1 siRNA impaired spatial learning in rats. Meanwhile, PIAS1WT increased STAT1 sumoylation, decreased STAT1 DNA binding and decreased STAT1 phosphorylation at Tyr-701 associated with spatial learning facilitation. But PIAS1 siRNA transfection produced an opposite effect on these measures associated with spatial learning impairment. Further, transfection of STAT1 sumoylation mutant impaired spatial acquisition, whereas transfection of STAT1 phosphorylation mutant blocked the impairing effect of PIAS1 siRNA on spatial learning. In this study, we first demonstrate the role of PIAS1 in spatial learning. Both posttranslational modifications (increased sumoylation and decreased phosphorylation) mediate the effect of PIAS1 on spatial learning facilitation.

SGK1 phosphorylation of IkappaB Kinase alpha and p300 Up-regulates NF-kappaB activity and increases N-Methyl-D-aspartate receptor NR2A and NR2B expression.

Serum- and glucocorticoid-inducible kinase 1 (SGK1) is a downstream target of phosphatidylinositol 3-kinase signaling, and it regulates various cellular and physiological functions, but the SGK1 substrate proteins and genes regulated by SGK1 are less known. Here we have identified IkappaB kinase alpha (IKKalpha) as a novel substrate of SGK1 by using biochemical and bioinformatic approaches. SGK1 directly phosphorylates IKKalpha at Thr-23 and indirectly activates IKKalpha at Ser-180. Furthermore, SGK1 enhanced nuclear factor kappaB (NF-kappaB) activity and up-regulated N-methyl-d-aspartate receptor NR2A and NR2B expression through activation of IKKalpha at Thr-23 and Ser-180, and these two residues play an equally important role in mediating these effects of SGK1. Although SGK1 does not phosphorylate IKKbeta, IKKbeta activity is still required for IKK complex activation and for SGK1 phosphorylation and activation of NF-kappaB. In addition, SGK1 increased the acetylation of NF-kappaB through phosphorylation of p300 at Ser-1834, and this also leads to NF-kappaB activation and NR2A and NR2B expression. Moreover, an endogenous stimulus of SGK1, insulin, increased IKKalpha and NF-kappaB phosphorylation as well as NF-kappaB acetylation and NF-kappaB activity, but SGK1 small interfering RNA transfection blocked these effects of insulin. In examination of the functional significance of the SGK1-IKKalpha-NF-kappaB signaling pathway, we found that transfection of the IKKalpha double mutant (IKKalphaT23A/S180A) to rat hippocampus antagonized SGK-1-mediated spatial memory facilitation. Our results together demonstrated novel substrate proteins of SGK1 and novel SGK1 signaling pathways. Activation of these signaling pathways enhances NR2A and NR2B expression that is implicated in neuronal plasticity.