Polygenic risk for schizophrenia converges on alternative polyadenylation as molecular mechanism underlying synaptic impairment.

Schizophrenia (SCZ) is a genetically heterogenous psychiatric disorder of highly polygenic nature. Correlative evidence from genetic studies indicate that the aggregated effects of distinct genetic risk factor combinations found in each patient converge onto common molecular mechanisms. To prove this on a functional level, we employed a reductionistic cellular model system for polygenic risk by differentiating induced pluripotent stem cells (iPSCs) from 104 individuals with high polygenic risk load and controls into cortical glutamatergic neurons (iNs). Multi-omics profiling identified widespread differences in alternative polyadenylation (APA) in the 3' untranslated region of many synaptic transcripts between iNs from SCZ patients and healthy donors. On the cellular level, 3'APA was associated with a reduction in synaptic density of iNs. Importantly, differential APA was largely conserved between postmortem human prefrontal cortex from SCZ patients and healthy donors, and strongly enriched for transcripts related to synapse biology. 3'APA was highly correlated with SCZ polygenic risk and affected genes were significantly enriched for SCZ associated common genetic variation. Integrative functional genomic analysis identified the RNA binding protein and SCZ GWAS risk gene PTBP2 as a critical trans-acting factor mediating 3'APA of synaptic genes in SCZ subjects. Functional characterization of PTBP2 in iNs confirmed its key role in 3'APA of synaptic transcripts and regulation of synapse density. Jointly, our findings show that the aggregated effects of polygenic risk converge on 3'APA as one common molecular mechanism that underlies synaptic impairments in SCZ.

Massively parallel functional dissection of schizophrenia-associated noncoding genetic variants.

Schizophrenia (SCZ) is a highly heritable mental disorder with thousands of associated genetic variants located mostly in the noncoding space of the genome. Translating these associations into insights regarding the underlying pathomechanisms has been challenging because the causal variants, their mechanisms of action, and their target genes remain largely unknown. We implemented a massively parallel variant annotation pipeline (MVAP) to perform SCZ variant-to-function mapping at scale in disease-relevant neural cell types. This approach identified 620 functional variants (1.7%) that operate in a highly developmental context and neuronal-activity-dependent manner. Multimodal integration of epigenomic and CRISPRi screening data enabled us to link these functional variants to target genes, biological processes, and ultimately alterations of neuronal physiology. These results provide a multistage prioritization strategy to map functional single-nucleotide polymorphism (SNP)-to-gene-to-endophenotype relations and offer biological insights into the context-dependent molecular processes modulated by SCZ-associated genetic variation.

Generation and molecular characterization of human pluripotent stem cell-derived pharyngeal foregut endoderm.

Approaches to study human pharyngeal foregut endoderm-a developmental intermediate that is linked to various human syndromes involving pharynx development and organogenesis of tissues such as thymus, parathyroid, and thyroid-have been hampered by scarcity of tissue access and cellular models. We present an efficient stepwise differentiation method to generate human pharyngeal foregut endoderm from pluripotent stem cells. We determine dose and temporal requirements of signaling pathway engagement for optimized differentiation and characterize the differentiation products on cellular and integrated molecular level. We present a computational classification tool, "CellMatch," and transcriptomic classification of differentiation products on an integrated mouse scRNA-seq developmental roadmap confirms cellular maturation. Integrated transcriptomic and chromatin analyses infer differentiation stage-specific gene regulatory networks. Our work provides the method and integrated multiomic resource for the investigation of disease-relevant loci and gene regulatory networks and their role in developmental defects affecting the pharyngeal endoderm and its derivatives.

Cell type and condition specific functional annotation of schizophrenia associated non-coding genetic variants.

Schizophrenia (SCZ) is a highly polygenic disease and genome wide association studies have identified thousands of genetic variants that are statistically associated with this psychiatric disorder. However, our ability to translate these associations into insights on the disease mechanisms has been challenging since the causal genetic variants, their molecular function and their target genes remain largely unknown. In order to address these questions, we established a functional genomics pipeline in combination with induced pluripotent stem cell technology to functionally characterize ~35,000 non-coding genetic variants associated with schizophrenia along with their target genes. This analysis identified a set of 620 (1.7%) single nucleotide polymorphisms as functional on a molecular level in a highly cell type and condition specific fashion. These results provide a high-resolution map of functional variant-gene combinations and offer comprehensive biological insights into the developmental context and stimulation dependent molecular processes modulated by SCZ associated genetic variation.

Distinct genetic liability profiles define clinically relevant patient strata across common diseases.

Genome-wide association studies have unearthed a wealth of genetic associations across many complex diseases. However, translating these associations into biological mechanisms contributing to disease etiology and heterogeneity has been challenging. Here, we hypothesize that the effects of disease-associated genetic variants converge onto distinct cell type specific molecular pathways within distinct subgroups of patients. In order to test this hypothesis, we develop the CASTom-iGEx pipeline to operationalize individual level genotype data to interpret personal polygenic risk and identify the genetic basis of clinical heterogeneity. The paradigmatic application of this approach to coronary artery disease and schizophrenia reveals a convergence of disease associated variant effects onto known and novel genes, pathways, and biological processes. The biological process specific genetic liabilities are not equally distributed across patients. Instead, they defined genetically distinct groups of patients, characterized by different profiles across pathways, endophenotypes, and disease severity. These results provide further evidence for a genetic contribution to clinical heterogeneity and point to the existence of partially distinct pathomechanisms across patient subgroups. Thus, the universally applicable approach presented here has the potential to constitute an important component of future personalized medicine concepts.

Associations of psychiatric disease and ageing with FKBP5 expression converge on superficial layer neurons of the neocortex.

Identification and characterisation of novel targets for treatment is a priority in the field of psychiatry. FKBP5 is a gene with decades of evidence suggesting its pathogenic role in a subset of psychiatric patients, with potential to be leveraged as a therapeutic target for these individuals. While it is widely reported that FKBP5/FKBP51 mRNA/protein (FKBP5/1) expression is impacted by psychiatric disease state, risk genotype and age, it is not known in which cell types and sub-anatomical areas of the human brain this occurs. This knowledge is critical to propel FKBP5/1-targeted treatment development. Here, we performed an extensive, large-scale postmortem study (n = 1024) of FKBP5/1, examining neocortical areas (BA9, BA11 and ventral BA24/BA24a) derived from subjects that lived with schizophrenia, major depression or bipolar disorder. With an extensive battery of RNA (bulk RNA sequencing, single-nucleus RNA sequencing, microarray, qPCR, RNAscope) and protein (immunoblot, immunohistochemistry) analysis approaches, we thoroughly investigated the effects of disease state, ageing and genotype on cortical FKBP5/1 expression including in a cell type-specific manner. We identified consistently heightened FKBP5/1 levels in psychopathology and with age, but not genotype, with these effects strongest in schizophrenia. Using single-nucleus RNA sequencing (snRNAseq; BA9 and BA11) and targeted histology (BA9, BA24a), we established that these disease and ageing effects on FKBP5/1 expression were most pronounced in excitatory superficial layer neurons of the neocortex, and this effect appeared to be consistent in both the granular and agranular areas examined. We then found that this increase in FKBP5 levels may impact on synaptic plasticity, as FKBP5 gex levels strongly and inversely correlated with dendritic mushroom spine density and brain-derived neurotrophic factor (BDNF) levels in superficial layer neurons in BA11. These findings pinpoint a novel cellular and molecular mechanism that has potential to open a new avenue of FKBP51 drug development to treat cognitive symptoms in psychiatric disorders.

Primary cilia and SHH signaling impairments in human and mouse models of Parkinson's disease.

Parkinson's disease (PD) as a progressive neurodegenerative disorder arises from multiple genetic and environmental factors. However, underlying pathological mechanisms remain poorly understood. Using multiplexed single-cell transcriptomics, we analyze human neural precursor cells (hNPCs) from sporadic PD (sPD) patients. Alterations in gene expression appear in pathways related to primary cilia (PC). Accordingly, in these hiPSC-derived hNPCs and neurons, we observe a shortening of PC. Additionally, we detect a shortening of PC in PINK1-deficient human cellular and mouse models of familial PD. Furthermore, in sPD models, the shortening of PC is accompanied by increased Sonic Hedgehog (SHH) signal transduction. Inhibition of this pathway rescues the alterations in PC morphology and mitochondrial dysfunction. Thus, increased SHH activity due to ciliary dysfunction may be required for the development of pathoetiological phenotypes observed in sPD like mitochondrial dysfunction. Inhibiting overactive SHH signaling may be a potential neuroprotective therapy for sPD.

Temporally divergent regulatory mechanisms govern neuronal diversification and maturation in the mouse and marmoset neocortex.

Mammalian neocortical neurons span one of the most diverse cell type spectra of any tissue. Cortical neurons are born during embryonic development, and their maturation extends into postnatal life. The regulatory strategies underlying progressive neuronal development and maturation remain unclear. Here we present an integrated single-cell epigenomic and transcriptional analysis of individual mouse and marmoset cortical neuron classes, spanning both early postmitotic stages of identity acquisition and later stages of neuronal plasticity and circuit integration. We found that, in both species, the regulatory strategies controlling early and late stages of pan-neuronal development diverge. Early postmitotic neurons use more widely shared and evolutionarily conserved molecular regulatory programs. In contrast, programs active during later neuronal maturation are more brain- and neuron-specific and more evolutionarily divergent. Our work uncovers a temporal shift in regulatory choices during neuronal diversification and maturation in both mice and marmosets, which likely reflects unique evolutionary constraints on distinct events of neuronal development in the neocortex.

Expression of Lineage Transcription Factors Identifies Differences in Transition States of Induced Human Oligodendrocyte Differentiation.

Oligodendrocytes (OLs) are critical for myelination and are implicated in several brain disorders. Directed differentiation of human-induced OLs (iOLs) from pluripotent stem cells can be achieved by forced expression of different combinations of the transcription factors SOX10 (S), OLIG2 (O), and NKX6.2 (N). Here, we applied quantitative image analysis and single-cell transcriptomics to compare different transcription factor (TF) combinations for their efficacy towards robust OL lineage conversion. Compared with S alone, the combination of SON increases the number of iOLs and generates iOLs with a more complex morphology and higher expression levels of myelin-marker genes. RNA velocity analysis of individual cells reveals that S generates a population of oligodendrocyte-precursor cells (OPCs) that appear to be more immature than those generated by SON and to display distinct molecular properties. Our work highlights that TFs for generating iOPCs or iOLs should be chosen depending on the intended application or research question, and that SON might be beneficial to study more mature iOLs while S might be better suited to investigate iOPC biology.

Single-Cell RNA-Sequencing-Based CRISPRi Screening Resolves Molecular Drivers of Early Human Endoderm Development.

Studies in vertebrates have outlined conserved molecular control of definitive endoderm (END) development. However, recent work also shows that key molecular aspects of human END regulation differ even from rodents. Differentiation of human embryonic stem cells (ESCs) to END offers a tractable system to study the molecular basis of normal and defective human-specific END development. Here, we interrogated dynamics in chromatin accessibility during differentiation of ESCs to END, predicting DNA-binding proteins that may drive this cell fate transition. We then combined single-cell RNA-seq with parallel CRISPR perturbations to comprehensively define the loss-of-function phenotype of those factors in END development. Following a few candidates, we revealed distinct impairments in the differentiation trajectories for mediators of TGFß signaling and expose a role for the FOXA2 transcription factor in priming human END competence for human foregut and hepatic END specification. Together, this single-cell functional genomics study provides high-resolution insight on human END development.

A stably self-renewing adult blood-derived induced neural stem cell exhibiting patternability and epigenetic rejuvenation.

Recent reports suggest that induced neurons (iNs), but not induced pluripotent stem cell (iPSC)-derived neurons, largely preserve age-associated traits. Here, we report on the extent of preserved epigenetic and transcriptional aging signatures in directly converted induced neural stem cells (iNSCs). Employing restricted and integration-free expression of SOX2 and c-MYC, we generated a fully functional, bona fide NSC population from adult blood cells that remains highly responsive to regional patterning cues. Upon conversion, low passage iNSCs display a profound loss of age-related DNA methylation signatures, which further erode across extended passaging, thereby approximating the DNA methylation age of isogenic iPSC-derived neural precursors. This epigenetic rejuvenation is accompanied by a lack of age-associated transcriptional signatures and absence of cellular aging hallmarks. We find iNSCs to be competent for modeling pathological protein aggregation and for neurotransplantation, depicting blood-to-NSC conversion as a rapid alternative route for both disease modeling and neuroregeneration.

Global delay in nascent strand DNA methylation.

Cytosine methylation is widespread among organisms and essential for mammalian development. In line with early postulations of an epigenetic role in gene regulation, symmetric CpG methylation can be mitotically propagated over many generations with extraordinarily high fidelity. Here, we combine BrdU labeling and immunoprecipitation with genome-wide bisulfite sequencing to explore the inheritance of cytosine methylation onto newly replicated DNA in human cells. Globally, we observe a pronounced lag between the copying of genetic and epigenetic information in embryonic stem cells that is reconsolidated within hours to accomplish faithful mitotic transmission. Populations of arrested cells show a global reduction of lag-induced intermediate CpG methylation when compared to proliferating cells, whereas sites of transcription factor engagement appear cell-cycle invariant. Alternatively, the cancer cell line HCT116 preserves global epigenetic heterogeneity independently of cell-cycle arrest. Taken together, our data suggest that heterogeneous methylation largely reflects asynchronous proliferation, but is intrinsic to actively engaged cis-regulatory elements and cancer.

Author Correction: Global delay in nascent strand DNA methylation.

Following online publication of this article, the Gene Expression Omnibus records corresponding to accession codes GSM2406773, MN-d6, and GSM2406772, MN-d14, listed in the data availability statement were deleted. The data are now available under accession codes GSM3039355, WGBS_hESC_WT_D6_R4 (MN day 6), and GSM3039351, WGBS_hESC_WT_D14_R4 (MN day 14), and the data availability statement has been updated with the new accession codes in the HTML and PDF versions of the article.

Dissecting the Functional Consequences of De Novo DNA Methylation Dynamics in Human Motor Neuron Differentiation and Physiology.

The somatic DNA methylation (DNAme) landscape is established early in development but remains highly dynamic within focal regions that overlap with gene regulatory elements. The significance of these dynamic changes, particularly in the central nervous system, remains unresolved. Here, we utilize a powerful human embryonic stem cell differentiation model for the generation of motor neurons (MNs) in combination with genetic mutations in the de novo DNAme machinery. We quantitatively dissect the role of DNAme in directing somatic cell fate with high-resolution genome-wide bisulfite-, bulk-, and single-cell-RNA sequencing. We find defects in neuralization and MN differentiation in DNMT3A knockouts (KO) that can be rescued by the targeting of DNAme to key developmental loci using catalytically inactive dCas9. We also find decreased dendritic arborization and altered electrophysiological properties in DNMT3A KO MNs. Our work provides a list of DNMT3A-regulated targets and a mechanistic link between de novo DNAme, cellular differentiation, and human MN function.

Genome-wide tracking of dCas9-methyltransferase footprints.

In normal mammalian development cytosine methylation is essential and is directed to specific regions of the genome. Despite notable advances through mapping its genome-wide distribution, studying the direct contribution of DNA methylation to gene and genome regulation has been limited by the lack of tools for its precise manipulation. Thus, combining the targeting capability of the CRISPR-Cas9 system with an epigenetic modifier has attracted interest in the scientific community. In contrast to profiling the genome-wide cleavage of a nuclease competent Cas9, tracing the global activity of a dead Cas9 (dCas9) methyltransferase fusion protein is challenging within a highly methylated genome. Here, we report the generation and use of an engineered, methylation depleted but maintenance competent mouse ES cell line and find surprisingly ubiquitous nuclear activity of dCas9-methyltransferases. Subsequent experiments in human somatic cells refine these observations and point to an important difference between genetic and epigenetic editing tools that require unique experimental considerations.

Genetic determinants and epigenetic effects of pioneer-factor occupancy.

Transcription factors (TFs) direct developmental transitions by binding to target DNA sequences, influencing gene expression and establishing complex gene-regultory networks. To systematically determine the molecular components that enable or constrain TF activity, we investigated the genomic occupancy of FOXA2, GATA4 and OCT4 in several cell types. Despite their classification as pioneer factors, all three TFs exhibit cell-type-specific binding, even when supraphysiologically and ectopically expressed. However, FOXA2 and GATA4 can be distinguished by low enrichment at loci that are highly occupied by these factors in alternative cell types. We find that expression of additional cofactors increases enrichment at a subset of these sites. Finally, FOXA2 occupancy and changes to DNA accessibility can occur in G1-arrested cells, but subsequent loss of DNA methylation requires DNA replication.

Targeted bisulfite sequencing of the dynamic DNA methylome.

BACKGROUND: The ability to measure DNA methylation precisely and efficiently continues to drive our understanding of this modification in development and disease. Whole genome bisulfite sequencing has the advantage of theoretically capturing all cytosines in the genome at single-nucleotide resolution, but it has a number of significant practical drawbacks that become amplified with increasing sample numbers. All other technologies capture only a fraction of the cytosines that show dynamic regulation across cell and tissue types. RESULTS: Here, we present a novel hybrid selection design focusing on loci with dynamic methylation that captures a large number of differentially methylated gene-regulatory elements. We benchmarked this assay against matched whole genome data and profiled 25 human tissue samples to explore its ability to detect differentially methylated regions. CONCLUSIONS: Our target capture design fills a major gap left by all other assays that exist to map DNA methylation. It maintains the ability to link cytosine methylation to genetic differences, the single-base resolution and the analysis of neighboring cytosines while notably reducing the cost per sample by focusing the sequencing effort on the most informative and relevant regions of the genome.

Epigenetic Memory Underlies Cell-Autonomous Heterogeneous Behavior of Hematopoietic Stem Cells.

Stem cells determine homeostasis and repair of many tissues and are increasingly recognized as functionally heterogeneous. To define the extent of-and molecular basis for-heterogeneity, we overlaid functional, transcriptional, and epigenetic attributes of hematopoietic stem cells (HSCs) at a clonal level using endogenous fluorescent tagging. Endogenous HSC had clone-specific functional attributes over time in vivo. The intra-clonal behaviors were highly stereotypic, conserved under the stress of transplantation, inflammation, and genotoxic injury, and associated with distinctive transcriptional, DNA methylation, and chromatin accessibility patterns. Further, HSC function corresponded to epigenetic configuration but not always to transcriptional state. Therefore, hematopoiesis under homeostatic and stress conditions represents the integrated action of highly heterogeneous clones of HSC with epigenetically scripted behaviors. This high degree of epigenetically driven cell autonomy among HSCs implies that refinement of the concepts of stem cell plasticity and of the stem cell niche is warranted.