Neuronal and glial 3D chromatin architecture informs the cellular etiology of brain disorders

SWI/SNF Complex Connects Signaling and Epigenetic State in Cells of Nervous System

SWI/SNF protein complexes are evolutionarily conserved epigenetic regulators described in all eukaryotes. In metameric animals, the complexes are involved in all processes occurring in the nervous system, from neurogenesis to higher brain functions. On the one hand, the range of roles is wide because the SWI/SNF complexes act universally by mobilizing the nucleosomes in a chromatin template at multiple loci throughout the genome. On the other hand, the complexes mediate the action of multiple signaling pathways that control most aspects of neural tissue development and function. The issues are discussed to provide insight into the molecular basis of the multifaceted role of SWI/SNFs in cell cycle regulation, DNA repair, activation of immediate-early genes, neurogenesis, and brain and connectome formation. An overview is additionally provided for the molecular basis of nervous system pathologies associated with the SWI/SNF complexes and their contribution to neuroinflammation and neurodegeneration. Finally, we discuss the idea that SWI/SNFs act as an integration platform to connect multiple signaling and genetic programs.

Extensive long-range polycomb interactions and weak compartmentalization are hallmarks of human neuronal 3D genome

Chromatin architecture regulates gene expression and shapes cellular identity, particularly in neuronal cells. Specifically, polycomb group (PcG) proteins enable establishment and maintenance of neuronal cell type by reorganizing chromatin into repressive domains that limit the expression of fate-determining genes and sustain distinct gene expression patterns in neurons. Here, we map the 3D genome architecture in neuronal and non-neuronal cells isolated from the Wernicke's area of four human brains and comprehensively analyze neuron-specific aspects of chromatin organization. We find that genome segregation into active and inactive compartments is greatly reduced in neurons compared to other brain cells. Furthermore, neuronal Hi-C maps reveal strong long-range interactions, forming a specific network of PcG-mediated contacts in neurons that is nearly absent in other brain cells. These interacting loci contain developmental transcription factors with repressed expression in neurons and other mature brain cells. But only in neurons, they are rich in bivalent promoters occupied by H3K4me3 histone modification together with H3K27me3, which points to a possible functional role of PcG contacts in neurons. Importantly, other layers of chromatin organization also exhibit a distinct structure in neurons, characterized by an increase in short-range interactions and a decrease in long-range ones.

Three-dimensional chromatin mapping of sensory neurons reveals that cohesin-dependent genomic domains are required for axonal regeneration

The in vivo three-dimensional genomic architecture of adult mature neurons at homeostasis and after medically relevant perturbations such as axonal injury remains elusive. Here we address this knowledge gap by mapping the three-dimensional chromatin architecture and gene expression programme at homeostasis and after sciatic nerve injury in wild-type and cohesin-deficient mouse sensory dorsal root ganglia neurons via combinatorial Hi-C and RNA-seq. We find that cohesin is required for the full induction of the regenerative transcriptional program, by organising 3D genomic domains required for the activation of regenerative genes. Importantly, loss of cohesin results in disruption of chromatin architecture at regenerative genes and severely impaired nerve regeneration. Together, these data provide an original three-dimensional chromatin map of adult sensory neurons in vivo and demonstrate a role for cohesin-dependent chromatin interactions in neuronal regeneration.

The Genesis of Schizophrenia: An Origin Story

Schizophrenia is routinely referred to as a neurodevelopmental disorder, but the role of brain development in a disorder typically diagnosed during early adult life is enigmatic. The authors revisit the neurodevelopmental model of schizophrenia with genomic insights from the most recent schizophrenia clinical genetic association studies, transcriptomic and epigenomic analyses from human postmortem brain studies, and analyses from cellular models that recapitulate neurodevelopment. Emerging insights into schizophrenia genetic risk continue to converge on brain development, particularly stages of early brain development, that may be perturbed to deviate from a typical, normative course, resulting in schizophrenia clinical symptomatology. As the authors explicate, schizophrenia genetic risk is likely dynamic and context dependent, with effects of genetic risk varying spatiotemporally, across the neurodevelopmental continuum. Optimizing therapeutic strategies for the heterogeneous collective of individuals with schizophrenia may likely be guided by leveraging markers of genetic risk and derivative functional insights, well before the emergence of psychosis. Ultimately, rather than a focus on therapeutic intervention during adolescence or adulthood, principles of prediction and prophylaxis in the pre- and perinatal and neonatal stages may best comport with the biology of schizophrenia to address the early-stage perturbations that alter the normative neurodevelopmental trajectory.

Single-cell mapping of cell-type specific chromatin architecture in the central nervous system

Determining how chromatin is structured in the nucleus is critical to studying its role in gene regulation. Recent advances in the analysis of single-cell chromatin architecture have considerably improved our understanding of cell-type-specific chromosome conformation and nuclear architecture. In this review, we discuss the methods used for analysis of 3D chromatin conformation, including sequencing-based methods, imaging-based techniques, and computational approaches. We further review the application of these methods in the study of the role of chromatin topology in neural development and disorders.

A multi-ancestry cerebral cortex transcriptome-wide association study identifies genes associated with smoking behaviors

Transcriptome-wide association studies (TWAS) have provided valuable insight in identifying genes that may impact cigarette smoking. Most of previous studies, however, mainly focused on European ancestry. Limited TWAS studies have been conducted across multiple ancestries to explore genes that may impact smoking behaviors. In this study, we used cis-eQTL data of cerebral cortex from multiple ancestries in MetaBrain, including European, East Asian, and African samples, as reference panels to perform multi-ancestry TWAS analyses on ancestry-matched GWASs of four smoking behaviors including smoking initiation, smoking cessation, age of smoking initiation, and number of cigarettes per day in GWAS & Sequencing Consortium of Alcohol and Nicotine use (GSCAN). Multiple-ancestry fine-mapping approach was conducted to identify credible gene sets associated with these four traits. Enrichment and module network analyses were further performed to explore the potential roles of these identified gene sets. A total of 719 unique genes were identified to be associated with at least one of the four smoking traits across ancestries. Among those, 249 genes were further prioritized as putative causal genes in multiple ancestry-based fine-mapping approach. Several well-known smoking-related genes, including PSMA4, IREB2, and CHRNA3, showed high confidence across ancestries. Some novel genes, e.g., TSPAN3 and ANK2, were also identified in the credible sets. The enrichment analysis identified a series of critical pathways related to smoking such as synaptic transmission and glutamate receptor activity. Leveraging the power of the latest multi-ancestry GWAS and eQTL data sources, this study revealed hundreds of genes and relevant biological processes related to smoking behaviors. These findings provide new insights for future functional studies on smoking behaviors.

Using a comprehensive atlas and predictive models to reveal the complexity and evolution of brain-active regulatory elements

Most genetic variants associated with psychiatric disorders are located in noncoding regions of the genome. To investigate their functional implications, we integrate epigenetic data from the PsychENCODE Consortium and other published sources to construct a comprehensive atlas of candidate brain cis-regulatory elements. Using deep learning, we model these elements' sequence syntax and predict how binding sites for lineage-specific transcription factors contribute to cell type-specific gene regulation in various types of glia and neurons. The elements' evolutionary history suggests that new regulatory information in the brain emerges primarily via smaller sequence mutations within conserved mammalian elements rather than entirely new human- or primate-specific sequences. However, primate-specific candidate elements, particularly those active during fetal brain development and in excitatory neurons and astrocytes, are implicated in the heritability of brain-related human traits. Additionally, we introduce PsychSCREEN, a web-based platform offering interactive visualization of PsychENCODE-generated genetic and epigenetic data from diverse brain cell types in individuals with psychiatric disorders and healthy controls.

Single-cell genomics and regulatory networks for 388 human brains

Single-cell genomics is a powerful tool for studying heterogeneous tissues such as the brain. Yet little is understood about how genetic variants influence cell-level gene expression. Addressing this, we uniformly processed single-nuclei, multiomics datasets into a resource comprising >2.8 million nuclei from the prefrontal cortex across 388 individuals. For 28 cell types, we assessed population-level variation in expression and chromatin across gene families and drug targets. We identified >550,000 cell type-specific regulatory elements and >1.4 million single-cell expression quantitative trait loci, which we used to build cell-type regulatory and cell-to-cell communication networks. These networks manifest cellular changes in aging and neuropsychiatric disorders. We further constructed an integrative model accurately imputing single-cell expression and simulating perturbations; the model prioritized ~250 disease-risk genes and drug targets with associated cell types.

Genetic regulation of cell type-specific chromatin accessibility shapes brain disease etiology

Nucleotide variants in cell type-specific gene regulatory elements in the human brain are risk factors for human disease. We measured chromatin accessibility in 1932 aliquots of sorted neurons and non-neurons from 616 human postmortem brains and identified 34,539 open chromatin regions with chromatin accessibility quantitative trait loci (caQTLs). Only 10.4% of caQTLs are shared between neurons and non-neurons, which supports cell type-specific genetic regulation of the brain regulome. Incorporating allele-specific chromatin accessibility improves statistical fine-mapping and refines molecular mechanisms that underlie disease risk. Using massively parallel reporter assays in induced excitatory neurons, we screened 19,893 brain QTLs and identified the functional impact of 476 regulatory variants. Combined, this comprehensive resource captures variation in the human brain regulome and provides insights into disease etiology.

From compartments to loops: understanding the unique chromatin organization in neuronal cells

The three-dimensional organization of the genome plays a central role in the regulation of cellular functions, particularly in the human brain. This review explores the intricacies of chromatin organization, highlighting the distinct structural patterns observed between neuronal and non-neuronal brain cells. We integrate findings from recent studies to elucidate the characteristics of various levels of chromatin organization, from differential compartmentalization and topologically associating domains (TADs) to chromatin loop formation. By defining the unique chromatin landscapes of neuronal and non-neuronal brain cells, these distinct structures contribute to the regulation of gene expression specific to each cell type. In particular, we discuss potential functional implications of unique neuronal chromatin organization characteristics, such as weaker compartmentalization, neuron-specific TAD boundaries enriched with active histone marks, and an increased number of chromatin loops. Additionally, we explore the role of Polycomb group (PcG) proteins in shaping cell-type-specific chromatin patterns. This review further emphasizes the impact of variations in chromatin architecture between neuronal and non-neuronal cells on brain development and the onset of neurological disorders. It highlights the need for further research to elucidate the details of chromatin organization in the human brain in order to unravel the complexities of brain function and the genetic mechanisms underlying neurological disorders. This research will help bridge a significant gap in our comprehension of the interplay between chromatin structure and cell functions.

Putative Risk Biomarkers of Bipolar Disorder in At-risk Youth

Bipolar disorder is a highly heritable and functionally impairing disease. The recognition and intervention of BD especially that characterized by early onset remains challenging. Risk biomarkers for predicting BD transition among at-risk youth may improve disease prognosis. We reviewed the more recent clinical studies to find possible pre-diagnostic biomarkers in youth at familial or (and) clinical risk of BD. Here we found that putative biomarkers for predicting conversion to BD include findings from multiple sample sources based on different hypotheses. Putative risk biomarkers shown by perspective studies are higher bipolar polygenetic risk scores, epigenetic alterations, elevated immune parameters, front-limbic system deficits, and brain circuit dysfunction associated with emotion and reward processing. Future studies need to enhance machine learning integration, make clinical detection methods more objective, and improve the quality of cohort studies.

Antipsychotic-induced epigenomic reorganization in frontal cortex of individuals with schizophrenia

Genome-wide association studies have revealed >270 loci associated with schizophrenia risk, yet these genetic factors do not seem to be sufficient to fully explain the molecular determinants behind this psychiatric condition. Epigenetic marks such as post-translational histone modifications remain largely plastic during development and adulthood, allowing a dynamic impact of environmental factors, including antipsychotic medications, on access to genes and regulatory elements. However, few studies so far have profiled cell-specific genome-wide histone modifications in postmortem brain samples from schizophrenia subjects, or the effect of antipsychotic treatment on such epigenetic marks. Here, we conducted ChIP-seq analyses focusing on histone marks indicative of active enhancers (H3K27ac) and active promoters (H3K4me3), alongside RNA-seq, using frontal cortex samples from antipsychotic-free (AF) and antipsychotic-treated (AT) individuals with schizophrenia, as well as individually matched controls (n=58). Schizophrenia subjects exhibited thousands of neuronal and non-neuronal epigenetic differences at regions that included several susceptibility genetic loci, such as NRG1, DISC1, and DRD3. By analyzing the AF and AT cohorts separately, we identified schizophrenia-associated alterations in specific transcription factors, their regulatees, and epigenomic and transcriptomic features that were reversed by antipsychotic treatment; as well as those that represented a consequence of antipsychotic medication rather than a hallmark of schizophrenia in postmortem human brain samples. Notably, we also found that the effect of age on epigenomic landscapes was more pronounced in frontal cortex of AT-schizophrenics, as compared to AF-schizophrenics and controls. Together, these data provide important evidence of epigenetic alterations in the frontal cortex of individuals with schizophrenia, and remark for the first time on the impact of age and antipsychotic treatment on chromatin organization.

Single-cell genomics and regulatory networks for 388 human brains

Single-cell genomics is a powerful tool for studying heterogeneous tissues such as the brain. Yet, little is understood about how genetic variants influence cell-level gene expression. Addressing this, we uniformly processed single-nuclei, multi-omics datasets into a resource comprising >2.8M nuclei from the prefrontal cortex across 388 individuals. For 28 cell types, we assessed population-level variation in expression and chromatin across gene families and drug targets. We identified >550K cell-type-specific regulatory elements and >1.4M single-cell expression-quantitative-trait loci, which we used to build cell-type regulatory and cell-to-cell communication networks. These networks manifest cellular changes in aging and neuropsychiatric disorders. We further constructed an integrative model accurately imputing single-cell expression and simulating perturbations; the model prioritized ~250 disease-risk genes and drug targets with associated cell types.

SATB2 organizes the 3D genome architecture of cognition in cortical neurons

The DNA-binding protein SATB2 is genetically linked to human intelligence. We studied its influence on the three-dimensional (3D) epigenome by mapping chromatin interactions and accessibility in control versus SATB2-deficient cortical neurons. We find that SATB2 affects the chromatin looping between enhancers and promoters of neuronal-activity-regulated genes, thus influencing their expression. It also alters A/B compartments, topologically associating domains, and frequently interacting regions. Genes linked to SATB2-dependent 3D genome changes are implicated in highly specialized neuronal functions and contribute to cognitive ability and risk for neuropsychiatric and neurodevelopmental disorders. Non-coding DNA regions with a SATB2-dependent structure are enriched for common variants associated with educational attainment, intelligence, and schizophrenia. Our data establish SATB2 as a cell-type-specific 3D genome modulator, which operates both independently and in cooperation with CCCTC-binding factor (CTCF) to set up the chromatin landscape of pyramidal neurons for cognitive processes.

Deciphering the Role of Rapidly Evolving Conserved Elements in Primate Brain Development and Exploring Their Potential Involvement in Alzheimer's Disease

Although previous studies have identified human-specific accelerated regions as playing a key role in the recent evolution of the human brain, the characteristics and cellular functions of rapidly evolving conserved elements (RECEs) in ancestral primate lineages remain largely unexplored. Here, based on large-scale primate genome assemblies, we identify 888 RECEs that have been highly conserved in primates that exhibit significantly accelerated substitution rates in the ancestor of the Simiiformes. This primate lineage exhibits remarkable morphological innovations, including an expanded brain mass. Integrative multiomic analyses reveal that RECEs harbor sequences with potential cis-regulatory functions that are activated in the adult human brain. Importantly, genes linked to RECEs exhibit pronounced expression trajectories in the adult brain relative to the fetal stage. Furthermore, we observed an increase in the chromatin accessibility of RECEs in oligodendrocytes from individuals with Alzheimer's disease (AD) compared to that of a control group, indicating that these RECEs may contribute to brain aging and AD. Our findings serve to expand our knowledge of the genetic underpinnings of brain function during primate evolution.

Epigenetic mechanisms underlying sex differences in the brain and behavior

Sex differences are found across brain regions, behaviors, and brain diseases. Sexual differentiation of the brain is initiated prenatally but it continues throughout life, as a result of the interaction of three major factors: gonadal hormones, sex chromosomes, and the environment. These factors are thought to act, in part, via epigenetic mechanisms which control chromatin and transcriptional states in brain cells. In this review, we discuss evidence that epigenetic mechanisms underlie sex-specific neurobehavioral changes during critical organizational periods, across the estrous cycle, and in response to diverse environments throughout life. We further identify future directions for the field that will provide novel mechanistic insights into brain sex differences, inform brain disease treatments and women's brain health in particular, and apply to people across genders.

Lineage specific 3D genome structure in the adult human brain and neurodevelopmental changes in the chromatin interactome

The human brain is a complex organ comprised of distinct cell types, and the contribution of the 3D genome to lineage specific gene expression remains poorly understood. To decipher cell type specific genome architecture, and characterize fine scale changes in the chromatin interactome across neural development, we compared the 3D genome of the human fetal cortical plate to that of neurons and glia isolated from the adult prefrontal cortex. We found that neurons have weaker genome compartmentalization compared to glia, but stronger TADs, which emerge during fetal development. Furthermore, relative to glia, the neuronal genome shifts more strongly towards repressive compartments. Neurons have differential TAD boundaries that are proximal to active promoters involved in neurodevelopmental processes. CRISPRi on CNTNAP2 in hIPSC-derived neurons reveals that transcriptional inactivation correlates with loss of insulation at the differential boundary. Finally, re-wiring of chromatin loops during neural development is associated with transcriptional and functional changes. Importantly, differential loops in the fetal cortex are associated with autism GWAS loci, suggesting a neuropsychiatric disease mechanism affecting the chromatin interactome. Furthermore, neural development involves gaining enhancer-promoter loops that upregulate genes that control synaptic activity. Altogether, our study provides multi-scale insights on the 3D genome in the human brain.

Deficiency of the Heterogeneous Nuclear Ribonucleoprotein U locus leads to delayed hindbrain neurogenesis

Genetic variants affecting Heterogeneous Nuclear Ribonucleoprotein U (HNRNPU) have been identified in several neurodevelopmental disorders (NDDs). HNRNPU is widely expressed in the human brain and shows the highest postnatal expression in the cerebellum. Recent studies have investigated the role of HNRNPU in cerebral cortical development, but the effects of HNRNPU deficiency on cerebellar development remain unknown. Here, we describe the molecular and cellular outcomes of HNRNPU locus deficiency during in vitro neural differentiation of patient-derived and isogenic neuroepithelial stem cells with a hindbrain profile. We demonstrate that HNRNPU deficiency leads to chromatin remodeling of A/B compartments, and transcriptional rewiring, partly by impacting exon inclusion during mRNA processing. Genomic regions affected by the chromatin restructuring and host genes of exon usage differences show a strong enrichment for genes implicated in epilepsies, intellectual disability, and autism. Lastly, we show that at the cellular level HNRNPU downregulation leads to an increased fraction of neural progenitors in the maturing neuronal population. We conclude that the HNRNPU locus is involved in delayed commitment of neural progenitors to differentiate in cell types with hindbrain profile.

Single-cell DNA methylation and 3D genome architecture in the human brain

Delineating the gene-regulatory programs underlying complex cell types is fundamental for understanding brain function in health and disease. Here, we comprehensively examined human brain cell epigenomes by probing DNA methylation and chromatin conformation at single-cell resolution in 517 thousand cells (399 thousand neurons and 118 thousand non-neurons) from 46 regions of three adult male brains. We identified 188 cell types and characterized their molecular signatures. Integrative analyses revealed concordant changes in DNA methylation, chromatin accessibility, chromatin organization, and gene expression across cell types, cortical areas, and basal ganglia structures. We further developed single-cell methylation barcodes that reliably predict brain cell types using the methylation status of select genomic sites. This multimodal epigenomic brain cell atlas provides new insights into the complexity of cell-type-specific gene regulation in adult human brains.

Systematic investigation of allelic regulatory activity of schizophrenia-associated common variants

Genome-wide association studies (GWASs) have successfully identified 145 genomic regions that contribute to schizophrenia risk, but linkage disequilibrium makes it challenging to discern causal variants. We performed a massively parallel reporter assay (MPRA) on 5,173 fine-mapped schizophrenia GWAS variants in primary human neural progenitors and identified 439 variants with allelic regulatory effects (MPRA-positive variants). Transcription factor binding had modest predictive power, while fine-map posterior probability, enhancer overlap, and evolutionary conservation failed to predict MPRA-positive variants. Furthermore, 64% of MPRA-positive variants did not exhibit expressive quantitative trait loci signature, suggesting that MPRA could identify yet unexplored variants with regulatory potentials. To predict the combinatorial effect of MPRA-positive variants on gene regulation, we propose an accessibility-by-contact model that combines MPRA-measured allelic activity with neuronal chromatin architecture.

Cellular specificity is key to deciphering epigenetic changes underlying Alzheimer's disease

Different cell types in the brain play distinct roles in Alzheimer's disease (AD) progression. Late onset AD (LOAD) is a complex disease, with a large genetic component, but many risk loci fall in non-coding genome regions. Epigenetics implicates the non-coding genome with control of gene expression. The epigenome is highly cell-type specific and dynamically responds to the environment. Therefore, epigenetic mechanisms are well placed to explain genetic and environmental factors that are associated with AD. However, given this cellular specificity, purified cell populations or single cells need to be profiled to avoid effect masking. Here we review the current state of cell-type specific genome-wide profiling in LOAD, covering DNA methylation (CpG, CpH, and hydroxymethylation), histone modifications, and chromatin changes. To date, these data reveal that distinct cell types contribute and react differently to AD progression through epigenetic alterations. This review addresses the current gap in prior bulk-tissue derived work by spotlighting cell-specific changes that govern the complex interplay of cells throughout disease progression and are critical in understanding and developing effective treatments for AD.

Neuronal DNA double-strand breaks lead to genome structural variations and 3D genome disruption in neurodegeneration

Persistent DNA double-strand breaks (DSBs) in neurons are an early pathological hallmark of neurodegenerative diseases including Alzheimer's disease (AD), with the potential to disrupt genome integrity. We used single-nucleus RNA-seq in human postmortem prefrontal cortex samples and found that excitatory neurons in AD were enriched for somatic mosaic gene fusions. Gene fusions were particularly enriched in excitatory neurons with DNA damage repair and senescence gene signatures. In addition, somatic genome structural variations and gene fusions were enriched in neurons burdened with DSBs in the CK-p25 mouse model of neurodegeneration. Neurons enriched for DSBs also had elevated levels of cohesin along with progressive multiscale disruption of the 3D genome organization aligned with transcriptional changes in synaptic, neuronal development, and histone genes. Overall, this study demonstrates the disruption of genome stability and the 3D genome organization by DSBs in neurons as pathological steps in the progression of neurodegenerative diseases.

Sex Significantly Impacts the Function of Major Depression-Linked Variants In Vivo

BACKGROUND

Genome-wide association studies have discovered blocks of common variants-likely transcriptional-regulatory-associated with major depressive disorder (MDD), though the functional subset and their biological impacts remain unknown. Likewise, why depression occurs in females more frequently than males is unclear. We therefore tested the hypothesis that risk-associated functional variants interact with sex and produce greater impact in female brains.

METHODS

We developed techniques to directly measure regulatory variant activity and sex interactions using massively parallel reporter assays in the mouse brain in vivo, in a cell type-specific manner, and applied these approaches to measure activity of >1000 variants from >30 MDD loci.

RESULTS

We identified extensive sex-by-allele effects in mature hippocampal neurons, suggesting that sex-differentiated impacts of genetic risk may underlie sex bias in disease. Unbiased informatics approaches indicated that functional MDD variants recurrently disrupt a number of transcription factor binding motifs, including those of sex hormone receptors. We confirmed a role for the latter by performing massively parallel reporter assays in neonatal mice on the day of birth (during a sex-differentiating hormone surge) and hormonally quiescent juveniles.

CONCLUSIONS

Our study provides novel insights into the influence of age, biological sex, and cell type on regulatory variant function and provides a framework for in vivo parallel assays to functionally define interactions between organismal variables such as sex and regulatory variation. Moreover, we experimentally demonstrate that a portion of the sex differences seen in MDD occurrence may be a product of sex-differentiated effects at associated regulatory variants.

Genetic insights into human cortical organization and development through genome-wide analyses of 2,347 neuroimaging phenotypes

Our understanding of the genetics of the human cerebral cortex is limited both in terms of the diversity and the anatomical granularity of brain structural phenotypes. Here we conducted a genome-wide association meta-analysis of 13 structural and diffusion magnetic resonance imaging-derived cortical phenotypes, measured globally and at 180 bilaterally averaged regions in 36,663 individuals and identified 4,349 experiment-wide significant loci. These phenotypes include cortical thickness, surface area, gray matter volume, measures of folding, neurite density and water diffusion. We identified four genetic latent structures and causal relationships between surface area and some measures of cortical folding. These latent structures partly relate to different underlying gene expression trajectories during development and are enriched for different cell types. We also identified differential enrichment for neurodevelopmental and constrained genes and demonstrate that common genetic variants associated with cortical expansion are associated with cephalic disorders. Finally, we identified complex interphenotype and inter-regional genetic relationships among the 13 phenotypes, reflecting the developmental differences among them. Together, these analyses identify distinct genetic organizational principles of the cortex and their correlates with neurodevelopment.

Chromatin loop dynamics during cellular differentiation are associated with changes to both anchor and internal regulatory features

Three-dimensional (3D) chromatin structure has been shown to play a role in regulating gene transcription during biological transitions. Although our understanding of loop formation and maintenance is rapidly improving, much less is known about the mechanisms driving changes in looping and the impact of differential looping on gene transcription. One limitation has been a lack of well-powered differential looping data sets. To address this, we conducted a deeply sequenced Hi-C time course of megakaryocyte development comprising four biological replicates and 6 billion reads per time point. Statistical analysis revealed 1503 differential loops. Gained loop anchors were enriched for AP-1 occupancy and were characterized by large increases in histone H3K27ac (over 11-fold) but relatively small increases in CTCF and RAD21 binding (1.26- and 1.23-fold, respectively). Linear modeling revealed that changes in histone H3K27ac, chromatin accessibility, and JUN binding were better correlated with changes in looping than RAD21 and almost as well correlated as CTCF. Changes to epigenetic features between-rather than at-boundaries were highly predictive of changes in looping. Together these data suggest that although CTCF and RAD21 may be the core machinery dictating where loops form, other features (both at the anchors and within the loop boundaries) may play a larger role than previously anticipated in determining the relative loop strength across cell types and conditions.

A chromEM-staining protocol optimized for cardiac tissue

Three-dimensional (3D) chromatin organization has a key role in defining the transcription program of cells during development. Its alteration is the cause of gene expression changes responsible for several diseases. Thus, we need new tools to study this aspect of gene expression regulation. To this end, ChromEM was recently developed: this is an electron-microscopy staining technique that selectively marks nuclear DNA without altering its structure and, thus, allows better visualization of 3D chromatin conformation. However, despite increasingly frequent application of this staining technique on cells, it has not yet been applied to visualize chromatin ultrastructure in tissues. Here, we provide a protocol to carry out ChromEM on myocardial tissue harvested from the left ventricles of C57BL/6J mice and use this in combination with transmission electron microscopy (TEM) to measure some morphological parameters of peripheral heterochromatin in cardiomyocytes. This protocol could also be used, in combination with electron tomography, to study 3D chromatin organization in cardiomyocytes in different aspects of heart pathobiology (e.g., heart development, cardiac aging, and heart failure) as well as help to set-up ChromEM in other tissues.

Evaluation of C4 Gene Copy Number in Pediatric Acute Neuropsychiatric Syndrome

Pediatric acute-onset neuropsychiatric syndrome (PANS) is an abrupt-onset neuropsychiatric disorder. PANS patients have an increased prevalence of comorbid autoimmune illness, most commonly arthritis. In addition, an estimated one-third of PANS patients present with low serum C4 protein, suggesting decreased production or increased consumption of C4 protein. To test the possibility that copy number (CN) variation contributes to risk of PANS illness, we compared mean total C4A and total C4B CN in ethnically matched subjects from PANS DNA samples and controls (192 cases and 182 controls). Longitudinal data from the Stanford PANS cohort (n = 121) were used to assess whether the time to juvenile idiopathic arthritis (JIA) or autoimmune disease (AI) onset was a function of total C4A or C4B CN. Lastly, we performed several hypothesis-generating analyses to explore the correlation between individual C4 gene variants, sex, specific genotypes, and age of PANS onset. Although the mean total C4A or C4B CN did not differ in PANS compared to controls, PANS patients with low C4B CN were at increased risk for subsequent JIA diagnosis (hazard ratio = 2.7, p value = 0.004). We also observed a possible increase in risk for AI in PANS patients and a possible correlation between lower C4B and PANS age of onset. An association between rheumatoid arthritis and low C4B CN has been reported previously. However, patients with PANS develop different types of JIA: enthesitis-related arthritis, spondyloarthritis, and psoriatic arthritis. This suggests that C4B plays a role that spans these arthritis types.

Network models of chromatin structure

Increasing numbers of datasets and experimental assays that capture the organization of chromatin inside the nucleus warrant an effort to develop tools to visualize and analyze these structures. Alongside polymer physics or constraint-based modeling, network theory approaches to describe 3D epigenome organization have gained in popularity. Representing genomic regions as nodes in a network enables visualization of 1D epigenomics datasets in the context of chromatin structure maps, while network theory metrics can be used to describe 3D epigenome organization and dynamics. In this review, we summarize the most salient applications of network theory to the study of chromatin contact maps, demonstrating its potential in revealing epigenomic patterns and relating them to cellular phenotypes.

Chromatin remodeler Activity-Dependent Neuroprotective Protein (ADNP) contributes to syndromic autism

BACKGROUND

Individuals affected with autism often suffer additional co-morbidities such as intellectual disability. The genes contributing to autism cluster on a relatively limited number of cellular pathways, including chromatin remodeling. However, limited information is available on how mutations in single genes can result in such pleiotropic clinical features in affected individuals. In this review, we summarize available information on one of the most frequently mutated genes in syndromic autism the Activity-Dependent Neuroprotective Protein (ADNP).

RESULTS

Heterozygous and predicted loss-of-function ADNP mutations in individuals inevitably result in the clinical presentation with the Helsmoortel-Van der Aa syndrome, a frequent form of syndromic autism. ADNP, a zinc finger DNA-binding protein has a role in chromatin remodeling: The protein is associated with the pericentromeric protein HP1, the SWI/SNF core complex protein BRG1, and other members of this chromatin remodeling complex and, in murine stem cells, with the chromodomain helicase CHD4 in a ChAHP complex. ADNP has recently been shown to possess R-loop processing activity. In addition, many additional functions, for instance, in association with cytoskeletal proteins have been linked to ADNP.

CONCLUSIONS

We here present an integrated evaluation of all current aspects of gene function and evaluate how abnormalities in chromatin remodeling might relate to the pleiotropic clinical presentation in individual"s" with Helsmoortel-Van der Aa syndrome.

Single nucleus multiomics identifies ZEB1 and MAFB as candidate regulators of Alzheimer's disease-specific cis-regulatory elements

Cell type-specific transcriptional differences between brain tissues from donors with Alzheimer's disease (AD) and unaffected controls have been well documented, but few studies have rigorously interrogated the regulatory mechanisms responsible for these alterations. We performed single nucleus multiomics (snRNA-seq plus snATAC-seq) on 105,332 nuclei isolated from cortical tissues from 7 AD and 8 unaffected donors to identify candidate cis-regulatory elements (CREs) involved in AD-associated transcriptional changes. We detected 319,861 significant correlations, or links, between gene expression and cell type-specific transposase accessible regions enriched for active CREs. Among these, 40,831 were unique to AD tissues. Validation experiments confirmed the activity of many regions, including several candidate regulators of APP expression. We identified ZEB1 and MAFB as candidate transcription factors playing important roles in AD-specific gene regulation in neurons and microglia, respectively. Microglia links were globally enriched for heritability of AD risk and previously identified active regulatory regions.

Genomic profiling of HIV-1 integration in microglia cells links viral integration to the topologically associated domains

HIV-1 encounters the hierarchically organized host chromatin to stably integrate and persist in anatomically distinct latent reservoirs. The contribution of genome organization in HIV-1 infection has been largely understudied across different HIV-1 targets. Here, we determine HIV-1 integration sites (ISs), associate them with chromatin and expression signatures at different genomic scales in a microglia cell model, and profile them together with the primary T cell reservoir. HIV-1 insertions into introns of actively transcribed genes with IS hotspots in genic and super-enhancers, characteristic of blood cells, are maintained in the microglia cell model. Genome organization analysis reveals dynamic CCCTC-binding factor (CTCF) clusters in cells with active and repressed HIV-1 transcription, whereas CTCF removal impairs viral integration. We identify CTCF-enriched topologically associated domain (TAD) boundaries with signatures of transcriptionally active chromatin as HIV-1 integration determinants in microglia and CD4+ T cells, highlighting the importance of host genome organization in HIV-1 infection.

Annotating genetic variants to target genes using H-MAGMA

An outstanding goal in modern genomics is to systematically predict the functional outcome of noncoding variation associated with complex traits. To address this, we developed Hi-C-coupled multi-marker analysis of genomic annotation (H-MAGMA), which builds on traditional MAGMA-a gene-based analysis tool that assigns variants to their target genes-by incorporating 3D chromatin configuration in assigning variants to their putative target genes. Applying this approach, we identified key biological pathways implicated in a wide range of brain disorders and showed its utility in complementing other functional genomic resources such as expression quantitative trait loci-based variant annotation. Here, we provide a detailed protocol for generating the H-MAGMA variant-gene annotation file by using chromatin interaction data from the adult human brain. In addition, we provide an example of how H-MAGMA is run by using genome-wide association study summary statistics of Parkinson's disease. Lastly, we generated variant-gene annotation files for 28 tissues and cell types, with the hope of contributing a resource for researchers studying a broad set of complex genetic disorders. H-MAGMA can be performed in <2 h for any cell type in which Hi-C data are available.

Genomic regulatory sequences in the pathogenesis of bipolar disorder

The lifetime prevalence of bipolar disorder is estimated to be about 2%. Epigenetics defines regulatory mechanisms that determine relatively stable patterns of gene expression by controlling all key steps, from DNA to messenger RNA to protein. This Mini Review highlights recent discoveries of modified epigenetic control resulting from genetic variants associated with bipolar disorder in genome-wide association studies. The revealed epigenetic abnormalities implicate gene transcription and post-transcriptional regulation. In the light of these discoveries, the Mini Review focuses on the genes PACS1, MCHR1, DCLK3, HAPLN4, LMAN2L, TMEM258, GNL3, LRRC57, CACNA1C, CACNA1D, and NOVA2 and their potential biological role in the pathogenesis of bipolar disorder. Molecular mechanisms under control of these genes do not translate into a unified picture and substantially more research is needed to fill the gaps in knowledge and to solve current limitations in prognosis and treatment of bipolar disorder. In conclusion, the genetic and functional studies confirm the complex nature of bipolar disorder and indicate future research directions to explore possible targeted treatment options, eventually working toward a personalized approach.

The Role of Epigenetics in Neuroinflammatory-Driven Diseases

Neurodegenerative disorders are characterized by the progressive loss of central and/or peripheral nervous system neurons. Within this context, neuroinflammation comes up as one of the main factors linked to neurodegeneration progression. In fact, neuroinflammation has been recognized as an outstanding factor for Alzheimer's disease (AD), amyotrophic lateral sclerosis (ALS), Parkinson's disease (PD), and multiple sclerosis (MS). Interestingly, neuroinflammatory diseases are characterized by dramatic changes in the epigenetic profile, which might provide novel prognostic and therapeutic factors towards neuroinflammatory treatment. Deep changes in DNA and histone methylation, along with histone acetylation and altered non-coding RNA expression, have been reported at the onset of inflammatory diseases. The aim of this work is to review the current knowledge on this field.

3D Genome Plasticity in Normal and Diseased Neurodevelopment

Non-random spatial organization of the chromosomal material inside the nuclei of brain cells emerges as an important regulatory layer of genome organization and function in health and disease. Here, we discuss how integrative approaches assessing chromatin in context of the 3D genome is providing new insights into normal and diseased neurodevelopment. Studies in primate (incl. human) and rodent brain have confirmed that chromosomal organization in neurons and glia undergoes highly dynamic changes during pre- and early postnatal development, with potential for plasticity across a much wider age window. For example, neuronal 3D genomes from juvenile and adult cerebral cortex and hippocampus undergo chromosomal conformation changes at hundreds of loci in the context of learning and environmental enrichment, viral infection, and neuroinflammation. Furthermore, locus-specific structural DNA variations, such as micro-deletions, duplications, repeat expansions, and retroelement insertions carry the potential to disrupt the broader epigenomic and transcriptional landscape far beyond the boundaries of the site-specific variation, highlighting the critical importance of long-range intra- and inter-chromosomal contacts for neuronal and glial function.

High-throughput CRISPRi and CRISPRa technologies in 3D genome regulation for neuropsychiatric diseases

Advances in genomics have led to the identification of many risk loci with hundreds of genes and thousands of DNA variants associated with neuropsychiatric disorders. A significant barrier to understanding the genetic underpinnings of complex diseases is the lack of functional characterization of risk genes and variants in biological systems relevant to human health and connecting disease-associated variants to pathological phenotypes. Characterizing gene and DNA variant functions requires genetic perturbations followed by molecular and cellular assays of neurobiological phenotypes. However, generating null or mutant alleles is low throughput, making it impossible to characterize disease-associated variants in large quantities efficiently. CRISPR interference (CRISPRi) and CRISPR activation (CRISPRa) screens can be leveraged to dissect the biological consequences of the tested genes and variants in their native context. Nevertheless, testing non-coding variants associated with complex diseases remains non-trivial. In this review, we first discuss the current challenges of interpreting the function of the non-coding genome and approaches to prioritizing disease-associated variants in the context of the 3D epigenome. Second, we provide a brief overview of high-throughput CRISPRi and CRISPRa screening strategies applicable for characterizing non-coding sequences in appropriate biological systems. Lastly, we discuss the promising prospects of using CRISPR-based technologies to dissect DNA sequences associated with neuropsychiatric diseases.

The three-dimensional landscape of cortical chromatin accessibility in Alzheimer's disease

To characterize the dysregulation of chromatin accessibility in Alzheimer's disease (AD), we generated 636 ATAC-seq libraries from neuronal and nonneuronal nuclei isolated from the superior temporal gyrus and entorhinal cortex of 153 AD cases and 56 controls. By analyzing a total of ~20 billion read pairs, we expanded the repertoire of known open chromatin regions (OCRs) in the human brain and identified cell-type-specific enhancer-promoter interactions. We show that interindividual variability in OCRs can be leveraged to identify cis-regulatory domains (CRDs) that capture the three-dimensional structure of the genome (3D genome). We identified AD-associated effects on chromatin accessibility, the 3D genome and transcription factor (TF) regulatory networks. For one of the most AD-perturbed TFs, USF2, we validated its regulatory effect on lysosomal genes. Overall, we applied a systematic approach to understanding the role of the 3D genome in AD. We provide all data as an online resource for widespread community-based analysis.

Implications of glial metabolic dysregulation in the pathophysiology of neurodegenerative diseases

Glial cells are the most abundant cells of the brain, outnumbering neurons. These multifunctional cells are crucial for maintaining brain homeostasis by providing trophic and nutritional support to neurons, sculpting synapses, and providing an immune defense. Glia are highly plastic and undergo both structural and functional alterations in response to changes in the brain microenvironment. Glial phenotypes are intimately regulated by underlying metabolic machinery, which dictates the effector functions of these cells. Altered brain energy metabolism and chronic neuroinflammation are common features of several neurodegenerative diseases. Microglia and astrocytes are the major glial cells fueling the ongoing neuroinflammatory process, exacerbating neurodegeneration. Distinct metabolic perturbations in microglia and astrocytes, including altered carbohydrate, lipid, and amino acid metabolism have been documented in neurodegenerative diseases. These disturbances aggravate the neurodegenerative process by potentiating the inflammatory activation of glial cells. This review covers the recent advances in the molecular aspects of glial metabolic changes in the pathophysiology of neurodegenerative diseases. Finally, we discuss studies exploiting glial metabolism as a potential therapeutic avenue in neurodegenerative diseases.

Focus on your locus with a massively parallel reporter assay

A growing number of variants associated with risk for neurodevelopmental disorders have been identified by genome-wide association and whole genome sequencing studies. As common risk variants often fall within large haplotype blocks covering long stretches of the noncoding genome, the causal variants within an associated locus are often unknown. Similarly, the effect of rare noncoding risk variants identified by whole genome sequencing on molecular traits is seldom known without functional assays. A massively parallel reporter assay (MPRA) is an assay that can functionally validate thousands of regulatory elements simultaneously using high-throughput sequencing and barcode technology. MPRA has been adapted to various experimental designs that measure gene regulatory effects of genetic variants within cis- and trans-regulatory elements as well as posttranscriptional processes. This review discusses different MPRA designs that have been or could be used in the future to experimentally validate genetic variants associated with neurodevelopmental disorders. Though MPRA has limitations such as it does not model genomic context, this assay can help narrow down the underlying genetic causes of neurodevelopmental disorders by screening thousands of sequences in one experiment. We conclude by describing future directions of this technique such as applications of MPRA for gene-by-environment interactions and pharmacogenetics.

A mini-review of the role of condensin in human nervous system diseases

Mitosis and meiosis are crucial life activities that transmit eukaryotic genetic information to progeny in a stable and orderly manner. The formation and appearance of chromosomes, which are derived from chromatin, are the preconditions and signs of mitosis. When entering mitosis, interphase loose chromatin is highly spiralized and folded to form compact chromosomes. In recent years, it has been found that in addition to the well-known DNA, histones, and topoisomerase, a large protein complex called condensin plays an important role in the process of chromosome formation. Numerous studies have shown that the abnormal function of condensin can lead to incomplete or excessive concentration of chromatin, as well as disorder of genome organization process, abnormal transmission of genetic information, and ultimately lead to various diseases of individual, especially in nervous system diseases. In this review, the biological function of condensin and the potential pathogenic mechanism of condensin in nervous system diseases are briefly summarized. Therefore, the investigation of these mechanisms makes a significant contribution to the understanding of those related diseases and provides new ideas for clinical treatments.

Convergence and Divergence of Rare Genetic Disorders on Brain Phenotypes: A Review

IMPORTANCE

Rare genetic disorders modulating gene expression-as exemplified by gene dosage disorders (GDDs)-represent a collectively common set of high-risk factors for neuropsychiatric illness. Research on GDDs is rapidly expanding because these variants have high effect sizes and a known genetic basis. Moreover, the prevalence of recurrent GDDs (encompassing aneuploidies and certain copy number variations) enables genetic-first phenotypic characterization of the same GDD across multiple individuals, thereby offering a unique window into genetic influences on the human brain and behavior. However, the rapid growth of GDD research has unveiled perplexing phenotypic convergences and divergences across genomic loci; while phenotypic profiles may be specifically associated with a genomic variant, individual behavioral and neuroimaging traits appear to be nonspecifically influenced by most GDDs.

OBSERVATIONS

This complexity is addressed by (1) providing an accessible survey of genotype-phenotype mappings across different GDDs, focusing on psychopathology, cognition, and brain anatomy, and (2) detailing both methodological and mechanistic sources for observed phenotypic convergences and divergences. This effort yields methodological recommendations for future comparative phenotypic research on GDDs as well as a set of new testable hypotheses regarding aspects of early brain patterning that might govern the complex mapping of genetic risk onto phenotypic variation in neuropsychiatric disorders.

CONCLUSIONS AND RELEVANCE

A roadmap is provided to boost accurate measurement and mechanistic interrogation of phenotypic convergence and divergence across multiple GDDs. Pursuing the questions posed by GDDs could substantially improve our taxonomical, neurobiological, and translational understanding of neuropsychiatric illness.

An atlas of endogenous DNA double-strand breaks arising during human neural cell fate determination

Endogenous DNA double-strand breaks (DSBs) occurring in neural cells have been implicated in the pathogenesis of neurodevelopmental disorders (NDDs). Currently, a genomic map of endogenous DSBs arising during human neurogenesis is missing. Here, we applied in-suspension Breaks Labeling In Situ and Sequencing (sBLISS), RNA-Seq, and Hi-C to chart the genomic landscape of DSBs and relate it to gene expression and genome architecture in 2D cultures of human neuroepithelial stem cells (NES), neural progenitor cells (NPC), and post-mitotic neural cells (NEU). Endogenous DSBs were enriched at the promoter and along the gene body of transcriptionally active genes, at the borders of topologically associating domains (TADs), and around chromatin loop anchors. NDD risk genes harbored significantly more DSBs in comparison to other protein-coding genes, especially in NEU cells. We provide sBLISS, RNA-Seq, and Hi-C datasets for each differentiation stage, and all the scripts needed to reproduce our analyses. Our datasets and tools represent a unique resource that can be harnessed to investigate the role of genome fragility in the pathogenesis of NDDs.

Measurement(s)	DNA Double Strand Break • Transcriptome • Whole Genome Sequencing • Chromosome conformation capture assay • Tumor Suppressor p53-Binding Protein 1
Technology Type(s)	sBLISS - Genome-wide detection of DNA double-strand breaks by in-suspension breaks labeling in situ and sequencing • RNA-seq of total RNA • copy number variation profiling assay • Hi-C • Immunocytochemistry
Factor Type(s)	Human Neural Cell Fate Specification
Sample Characteristic - Organism	Homo sapiens
Sample Characteristic - Environment	cell culture

Sex hormone fluctuation and increased female risk for depression and anxiety disorders: From clinical evidence to molecular mechanisms

Women are at twice the risk for anxiety and depression disorders as men are, although the underlying biological factors and mechanisms are largely unknown. In this review, we address this sex disparity at both the etiological and mechanistic level. We dissect the role of fluctuating sex hormones as a critical biological factor contributing to the increased depression and anxiety risk in women. We provide parallel evidence in humans and rodents that brain structure and function vary with naturally-cycling ovarian hormones. This female-unique brain plasticity and associated vulnerability are primarily driven by estrogen level changes. For the first time, we provide a sex hormone-driven molecular mechanism, namely chromatin organizational changes, that regulates neuronal gene expression and brain plasticity but may also prime the (epi)genome for psychopathology. Finally, we map out future directions including experimental and clinical studies that will facilitate novel sex- and gender-informed approaches to treat depression and anxiety disorders.

Chromatin architecture in addiction circuitry identifies risk genes and potential biological mechanisms underlying cigarette smoking and alcohol use traits

Cigarette smoking and alcohol use are among the most prevalent substances used worldwide and account for a substantial proportion of preventable morbidity and mortality, underscoring the public health significance of understanding their etiology. Genome-wide association studies (GWAS) have successfully identified genetic variants associated with cigarette smoking and alcohol use traits. However, the vast majority of risk variants reside in non-coding regions of the genome, and their target genes and neurobiological mechanisms are unknown. Chromosomal conformation mappings can address this knowledge gap by charting the interaction profiles of risk-associated regulatory variants with target genes. To investigate the functional impact of common variants associated with cigarette smoking and alcohol use traits, we applied Hi-C coupled MAGMA (H-MAGMA) built upon cortical and newly generated midbrain dopaminergic neuronal Hi-C datasets to GWAS summary statistics of nicotine dependence, cigarettes per day, problematic alcohol use, and drinks per week. The identified risk genes mapped to key pathways associated with cigarette smoking and alcohol use traits, including drug metabolic processes and neuronal apoptosis. Risk genes were highly expressed in cortical glutamatergic, midbrain dopaminergic, GABAergic, and serotonergic neurons, suggesting them as relevant cell types in understanding the mechanisms by which genetic risk factors influence cigarette smoking and alcohol use. Lastly, we identified pleiotropic genes between cigarette smoking and alcohol use traits under the assumption that they may reveal substance-agnostic, shared neurobiological mechanisms of addiction. The number of pleiotropic genes was ~26-fold higher in dopaminergic neurons than in cortical neurons, emphasizing the critical role of ascending dopaminergic pathways in mediating general addiction phenotypes. Collectively, brain region- and neuronal subtype-specific 3D genome architecture helps refine neurobiological hypotheses for smoking, alcohol, and general addiction phenotypes by linking genetic risk factors to their target genes.

Understanding Regulatory Mechanisms of Brain Function and Disease through 3D Genome Organization

The human genome has a complex and dynamic three-dimensional (3D) organization, which plays a critical role for gene regulation and genome function. The importance of 3D genome organization in brain development and function has been well characterized in a region- and cell-type-specific fashion. Recent technological advances in chromosome conformation capture (3C)-based techniques, imaging approaches, and ligation-free methods, along with computational methods to analyze the data generated, have revealed 3D genome features at different scales in the brain that contribute to our understanding of genetic mechanisms underlying neuropsychiatric diseases and other brain-related traits. In this review, we discuss how these advances aid in the genetic dissection of brain-related traits.

Advances in profiling chromatin architecture shed light on the regulatory dynamics underlying brain disorders

Understanding the exquisitely complex nature of the three-dimensional organization of the genome and how it affects gene regulation remains a central question in biology. Recent advances in sequencing- and imaging-based approaches in decoding the three-dimensional chromatin landscape have enabled a systematic characterization of gene regulatory architecture. In this review, we outline how chromatin architecture provides a reference atlas to predict the functional consequences of non-coding variants associated with human traits and disease. High-throughput perturbation assays such as massively parallel reporter assays (MPRA) and CRISPR-based genome engineering in combination with a reference atlas opened an avenue for going beyond observational studies to experimentally validating the regulatory principles of the genome. We conclude by providing a suggested path forward by calling attention to barriers that can be addressed for a more complete understanding of the regulatory landscape of the human brain.

Limited Association between Schizophrenia Genetic Risk Factors and Transcriptomic Features

Schizophrenia is a polygenic disorder with many genomic regions contributing to schizophrenia risk. The majority of genetic variants associated with schizophrenia lie in the non-coding genome and are thought to contribute to transcriptional regulation. Extensive transcriptomic dysregulation has been detected from postmortem brain samples of schizophrenia-affected individuals. However, the relationship between schizophrenia genetic risk factors and transcriptomic features has yet to be explored. Herein, we examined whether varying gene expression features, including differentially expressed genes (DEGs), co-expression networks, and central hubness of genes, contribute to the heritability of schizophrenia. We leveraged quantitative trait loci and chromatin interaction profiles to identify schizophrenia risk variants assigned to the genes that represent different transcriptomic features. We then performed stratified linkage disequilibrium score regression analysis on these variants to estimate schizophrenia heritability enrichment for different gene expression features. Notably, DEGs and co-expression networks showed nominal heritability enrichment. This nominal association can be partly explained by cellular heterogeneity, as DEGs were associated with the genetic risk of schizophrenia in a cell type-specific manner. Moreover, DEGs were enriched for target genes of schizophrenia-associated transcription factors, suggesting that the transcriptomic signatures of schizophrenia are the result of transcriptional regulatory cascades elicited by genetic risk factors.