Robust Hi-C Maps of Enhancer-Promoter Interactions Reveal the Function of Non-coding Genome in Neural Development and Diseases

Unfolding neural diversity: how dynamic three-dimensional genome architecture regulates brain function and disease

The advent of single cell multi-omic technologies has ushered in a revolution in how we study the impact of three-dimensional genome organization on brain cellular composition and function. Transcriptomic and epigenomic studies reveal enormous cellular diversity that is present in mammalian nervous systems, raising the question, "how does this diversity arise and for what is its use?" Advances in the field of three-dimensional nuclear architecture have illuminated our understanding of how genome folding gives rise to dynamic gene expression programs important in healthy brain function and in disease. In this review we highlight recent work defining how neuronal identity, maturation, and plasticity are shaped by genome architecture. We discuss how newly identified genetic variations influence genome architecture and contribute to the evolution of species-unique neuronal and behavioral functional traits. We include examples for both humans and model organisms in which maladaptive genomic architecture is a causal agent in disease. Finally, we make conclusions and address future perspectives of dynamic three-dimensional genome (4D nucelome) research.

Heat shock induces alternative polyadenylation through dynamic DNA methylation and chromatin looping

Alternative cleavage and polyadenylation (APA) is a gene regulatory mechanism used by cells under stress to upregulate proteostasis-promoting transcripts, but how cells achieve this remains poorly understood. Previously, we elucidated a DNA methylation-regulated APA mechanism, in which gene body DNA methylation enhances distal poly(A) isoform expression by blocking CTCF binding and chromatin loop formation at APA control regions. We hypothesized that DNA methylation-regulated APA is one mechanism cells employ to induce proteostasis-promoting poly(A) isoforms. At the DNAJB6 co-chaperone locus, acute heat shock resulted in binding of stress response transcription factors HSF1, ATF6, and YY1 at the APA control region and an increase in the expression of the proximal poly(A) isoform known to prevent protein aggregation. Furthermore, TET1 was recruited to rapidly demethylate DNA, facilitating CTCF binding and chromatin loop formation, thereby reinforcing preferential proximal poly(A) isoform expression. As cells recovered, the transcription factors vacated the APA control region, and DNMT1 was recruited to remethylate the region. This process resolved chromatin looping and reset the poly(A) isoform expression pattern. Our findings unveil an epigenetic mechanism enabling cells to dynamically modulate poly(A) isoforms in response to stress while shedding light on the interplay between DNA methylation, transcription factor binding, and chromatin looping.

Spatiotemporal 3D chromatin organization across multiple brain regions during human fetal development

Elucidating the regulatory mechanisms underlying the development of different brain regions in humans is essential for understanding advanced cognition and neuropsychiatric disorders. However, the spatiotemporal organization of three-dimensional (3D) chromatin structure and its regulatory functions across different brain regions remain poorly understood. Here, we generated an atlas of high-resolution 3D chromatin structure across six developing human brain regions, including the prefrontal cortex (PFC), primary visual cortex (V1), cerebellum (CB), subcortical corpus striatum (CS), thalamus (TL), and hippocampus (HP), spanning gestational weeks 11-26. We found that the spatial and temporal dynamics of 3D chromatin organization play a key role in regulating brain region development. We also identified H3K27ac-marked super-enhancers as key contributors to shaping brain region-specific 3D chromatin structures and gene expression patterns. Finally, we uncovered hundreds of neuropsychiatric GWAS SNP-linked genes, shedding light on critical molecules in various neuropsychiatric disorders. In summary, our findings provide important insights into the 3D chromatin regulatory mechanisms governing brain region-specific development and can serve as a valuable resource for advancing our understanding of neuropsychiatric disorders.

Significance in scale space for Hi-C data

MOTIVATION

Hi-C technology has been developed to profile genome-wide chromosome conformation. So far Hi-C data have been generated from a large compendium of different cell types and different tissue types. Among different chromatin conformation units, chromatin loops were found to play a key role in gene regulation across different cell types. While many different loop calling algorithms have been developed, most loop callers identified shared loops as opposed to cell-type-specific loops.

RESULTS

We propose SSSHiC, a new loop calling algorithm based on significance in scale space, which can be used to understand data at different levels of resolution. By applying SSSHiC to neuronal and glial Hi-C data, we detected more loops that are potentially engaged in cell-type-specific gene regulation. Compared with other loop callers, such as Mustache, these loops were more frequently anchored to gene promoters of cellular marker genes and had better APA scores. Therefore, our results suggest that SSSHiC can effectively capture loops that contain more gene regulatory information.

AVAILABILITY AND IMPLEMENTATION

The Hi-C data used in this study can be accessed through the PsychENCODE Knowledge Portal at https://www.synapse.org/#! Synapse: syn21760712. The code utilized for Curvature SSS cited in this study is available at https://github.com/jsmarron/MarronMatlabSoftware/blob/master/Matlab9/Matlab9Combined.zip. All custom code used in this research can be found in the GitHub repository: https://github.com/jerryliu01998/HiC. The code has also been submitted to Code Ocean with the doi: 10.24433/CO.1912913.v1.

Transcriptional enhancers in human neuronal differentiation provide clues to neuronal disorders

Genome-wide association studies (GWASs) have identified thousands of variants associated with complex phenotypes, including neuropsychiatric disorders. To better understand their pathogenesis, it is necessary to identify the functional roles of these variants, which are largely located in non-coding DNA regions. Here, we employ a human mesencephalic neuronal cell differentiation model, LUHMES, with sensitive and high-resolution methods to discover enhancers (NET-CAGE), perform DNA conformation analysis (Capture Hi-C) to link enhancers to their target genes, and finally validate selected interactions. We expand the number of known enhancers active in differentiating human LUHMES neurons to 47,350, and find overlap with GWAS variants for Parkinson's disease and schizophrenia. Our findings reveal a fine-tuned regulation of human neuronal differentiation, even between adjacent developmental stages; provide a valuable resource for further studies on neuronal development, regulation, and disorders; and emphasize the importance of exploring the vast regulatory potential of non-coding DNA and enhancers.

Cell type-specific 3D-genome organization and transcription regulation in the brain

3D organization of the genome plays a critical role in regulating gene expression. How 3D-genome organization differs among different cell types and relates to cell type-dependent transcriptional regulation remains unclear. Here, we used genome-scale DNA and RNA imaging to investigate 3D-genome organization in transcriptionally distinct cell types in the mouse cerebral cortex. We uncovered a wide spectrum of differences in the nuclear architecture and 3D-genome organization among different cell types, ranging from the size of the cell nucleus to higher-order chromosome structures and radial positioning of chromatin loci within the nucleus. These cell type-dependent variations in nuclear architecture and chromatin organization exhibit strong correlations with both the total transcriptional activity of the cell and transcriptional regulation of cell type-specific marker genes. Moreover, we found that the methylated DNA binding protein MeCP2 promotes active-inactive chromatin segregation and regulates transcription in a nuclear radial position-dependent manner that is highly correlated with its function in modulating active-inactive chromatin compartmentalization.

Structural basis of differential gene expression at eQTLs loci from high-resolution ensemble models of 3D single-cell chromatin conformations

MOTIVATION

Techniques such as high-throughput chromosome conformation capture (Hi-C) have provided a wealth of information on nucleus organization and genome important for understanding gene expression regulation. Genome-Wide Association Studies have identified numerous loci associated with complex traits. Expression quantitative trait loci (eQTL) studies have further linked the genetic variants to alteration in expression levels of associated target genes across individuals. However, the functional roles of many eQTLs in noncoding regions remain unclear. Current joint analyses of Hi-C and eQTLs data lack advanced computational tools, limiting what can be learned from these data.

RESULTS

We developed a computational method for simultaneous analysis of Hi-C and eQTL data, capable of identifying a small set of nonrandom interactions from all Hi-C interactions. Using these nonrandom interactions, we reconstructed large ensembles (×105) of high-resolution single-cell 3D chromatin conformations with thorough sampling, accurately replicating Hi-C measurements. Our results revealed many-body interactions in chromatin conformation at the single-cell level within eQTL loci, providing a detailed view of how 3D chromatin structures form the physical foundation for gene regulation, including how genetic variants of eQTLs affect the expression of associated eGenes. Furthermore, our method can deconvolve chromatin heterogeneity and investigate the spatial associations of eQTLs and eGenes at subpopulation level, revealing their regulatory impacts on gene expression. Together, ensemble modeling of thoroughly sampled single-cell chromatin conformations combined with eQTL data, helps decipher how 3D chromatin structures provide the physical basis for gene regulation, expression control, and aid in understanding the overall structure-function relationships of genome organization.

AVAILABILITY AND IMPLEMENTATION

It is available at https://github.com/uic-liang-lab/3DChromFolding-eQTL-Loci.

Regulation of epigenetics and chromosome structure by human ORC2

The six subunit Origin Recognition Complex (ORC) is a DNA replication initiator that also promotes heterochromatinization in some species. A multi-omics study in a human cell line with mutations in three subunits of ORC, reveals that the subunits bind to DNA independent of each other rather than as part of a common six-subunit ORC. While DNA-bound ORC2 was seen to compact chromatin and attract repressive histone marks, the activation of chromatin and protection from repressive marks was seen at a large number of sites. The epigenetic changes regulate hundreds of genes, including some epigenetic regulators, adding an indirect mechanism by which ORC2 regulates epigenetics without local binding. DNA-bound ORC2 also prevents the acquisition of CTCF at focal sites in the genome to regulate chromatin loops. Thus, individual ORC subunits are major regulators, in both directions, of epigenetics, gene expression and chromosome structure, independent of the role of ORC in replication.

Convergent pairs of highly transcribed genes restrict chromatin looping in Dictyostelium discoideum

Dictyostelium discoideum is a unicellular slime mold, developing into a multicellular fruiting body upon starvation. Development is accompanied by large-scale shifts in gene expression program, but underlying features of chromatin spatial organization remain unknown. Here, we report that the Dictyostelium 3D genome is organized into positionally conserved, largely consecutive, non-hierarchical and weakly insulated loops at the onset of multicellular development. The transcription level within the loop interior tends to be higher than in adjacent regions. Loop interiors frequently contain functionally linked genes and genes which coherently change expression level during development. Loop anchors are predominantly positioned by the genes in convergent orientation. Results of polymer simulations and Hi-C-based observations suggest that the loop profile may arise from the interplay between transcription and extrusion-driven chromatin folding. In this scenario, a convergent gene pair serves as a bidirectional extrusion barrier or a 'diode' that controls passage of the cohesin extruder by relative transcription level of paired genes.

The interplay between histone modifications and nuclear lamina in genome regulation

Gene expression is regulated by chromatin architecture and epigenetic remodeling in cell homeostasis and pathologies. Histone modifications act as the key factors to modulate the chromatin accessibility. Different histone modifications are strongly associated with the localization of chromatin. Heterochromatin primarily localizes at the nuclear periphery, where it interacts with lamina proteins to suppress gene expression. In this review, we summarize the potential bridges that have regulatory functions of histone modifications in chromatin organization and transcriptional regulation at the nuclear periphery. We use lamina-associated domains (LADs) as examples to elucidate the biological roles of the interactions between histone modifications and nuclear lamina in cell differentiation and development. In the end, we highlight the technologies that are currently used to identify and visualize histone modifications and LADs, which could provide spatiotemporal information for understanding their regulatory functions in gene expression and discovering new targets for diseases.

Technological advancements in functional interpretation of genome-wide association studies (GWAS) findings: bridging the gap to clinical translation

Genome-wide association studies (GWAS) significantly advanced our understanding of the genetic underpinnings of diseases. However, challenges persist, particularly in interpreting non-coding variants in linkage disequilibrium that affect genes in disease-relevant cells. Addressing key obstacles-identifying causal variants, uncovering target genes, and understanding their network impact-is crucial. This graphical review navigates advanced techniques to fully leverage GWAS for future therapeutic breakthroughs.

The 3D Genome in Brain Development: An Exploration of Molecular Mechanisms and Experimental Methods

The human brain contains multiple cell types that are spatially organized into functionally distinct regions. The proper development of the brain requires complex gene regulation mechanisms in both neurons and the non-neuronal cell types that support neuronal function. Studies across the last decade have discovered that the 3D nuclear organization of the genome is instrumental in the regulation of gene expression in the diverse cell types of the brain. In this review, we describe the fundamental biochemical mechanisms that regulate the 3D genome, and comprehensively describe in vitro and ex vivo studies on mouse and human brain development that have characterized the roles of the 3D genome in gene regulation. We highlight the significance of the 3D genome in linking distal enhancers to their target promoters, which provides insights on the etiology of psychiatric and neurological disorders, as the genetic variants associated with these disorders are primarily located in noncoding regulatory regions. We also describe the molecular mechanisms that regulate chromatin folding and gene expression in neurons. Furthermore, we describe studies with an evolutionary perspective, which have investigated features that are conserved from mice to human, as well as human gained 3D chromatin features. Although most of the insights on disease and molecular mechanisms have been obtained from bulk 3C based experiments, we also highlight other approaches that have been developed recently, such as single cell 3C approaches, as well as non-3C based approaches. In our future perspectives, we highlight the gaps in our current knowledge and emphasize the need for 3D genome engineering and live cell imaging approaches to elucidate mechanisms and temporal dynamics of chromatin interactions, respectively.

Multi-omic analysis of longitudinal acute myeloid leukemia patient samples reveals potential prognostic markers linked to disease progression

Relapse remains a determinant of treatment failure and contributes significantly to mortality in acute myeloid leukemia (AML) patients. Despite efforts to understand AML progression and relapse mechanisms, findings on acquired gene mutations in relapse vary, suggesting inherent genetic heterogeneity and emphasizing the role of epigenetic modifications. We conducted a multi-omic analysis using Omni-C, ATAC-seq, and RNA-seq on longitudinal samples from two adult AML patients at diagnosis and relapse. Herein, we characterized genetic and epigenetic changes in AML progression to elucidate the underlying mechanisms of relapse. Differential interaction analysis showed significant 3D chromatin landscape reorganization between relapse and diagnosis samples. Comparing global open chromatin profiles revealed that relapse samples had significantly fewer accessible chromatin regions than diagnosis samples. In addition, we discovered that relapse-related upregulation was achieved either by forming new active enhancer contacts or by losing interactions with poised enhancers/potential silencers. Altogether, our study highlights the impact of genetic and epigenetic changes on AML progression, underlining the importance of multi-omic approaches in understanding disease relapse mechanisms and guiding potential therapeutic interventions.

Condensate interfacial forces reposition DNA loci and probe chromatin viscoelasticity

Biomolecular condensates assemble in living cells through phase separation and related phase transitions. An underappreciated feature of these dynamic molecular assemblies is that they form interfaces with other cellular structures, including membranes, cytoskeleton, DNA and RNA, and other membraneless compartments. These interfaces are expected to give rise to capillary forces, but there are few ways of quantifying and harnessing these forces in living cells. Here, we introduce viscoelastic chromatin tethering and organization (VECTOR), which uses light-inducible biomolecular condensates to generate capillary forces at targeted DNA loci. VECTOR can be utilized to programmably reposition genomic loci on a timescale of seconds to minutes, quantitatively revealing local heterogeneity in the viscoelastic material properties of chromatin. These synthetic condensates are built from components that naturally form liquid-like structures in living cells, highlighting the potential role for native condensates to generate forces and do work to reorganize the genome and impact chromatin architecture.

DKK1-SE recruits AP1 to activate the target gene DKK1 thereby promoting pancreatic cancer progression

Super-enhancers are a class of DNA cis-regulatory elements that can regulate cell identity, cell fate, stem cell pluripotency, and even tumorigenesis. Increasing evidence shows that epigenetic modifications play an important role in the pathogenesis of various types of cancer. However, the current research is far from enough to reveal the complex mechanism behind it. This study found a super-enhancer enriched with abnormally active histone modifications in pancreatic ductal adenocarcinoma (PDAC), called DKK1-super-enhancer (DKK1-SE). The major active component of DKK1-SE is component enhancer e1. Mechanistically, AP1 induces chromatin remodeling in component enhancer e1 and activates the transcriptional activity of DKK1. Moreover, DKK1 was closely related to the malignant clinical features of PDAC. Deletion or knockdown of DKK1-SE significantly inhibited the proliferation, colony formation, motility, migration, and invasion of PDAC cells in vitro, and these phenomena were partly mitigated upon rescuing DKK1 expression. In vivo, DKK1-SE deficiency not only inhibited tumor proliferation but also reduced the complexity of the tumor microenvironment. This study identifies that DKK1-SE drives DKK1 expression by recruiting AP1 transcription factors, exerting oncogenic effects in PDAC, and enhancing the complexity of the tumor microenvironment.

Extended genome-wide association study employing the African genome resources panel identifies novel susceptibility loci for Alzheimer's disease in individuals of African ancestry

INTRODUCTION

Despite a two-fold risk, individuals of African ancestry have been underrepresented in Alzheimer's disease (AD) genomics efforts.

METHODS

Genome-wide association studies (GWAS) of 2,903 AD cases and 6,265 controls of African ancestry. Within-dataset results were meta-analyzed, followed by functional genomics analyses.

RESULTS

A novel AD-risk locus was identified in MPDZ on chromosome (chr) 9p23 (rs141610415, MAF = 0.002, p = 3.68×10-9). Two additional novel common and nine rare loci were identified with suggestive associations (P < 9×10-7). Comparison of association and linkage disequilibrium (LD) patterns between datasets with higher and lower degrees of African ancestry showed differential association patterns at chr12q23.2 (ASCL1), suggesting that this association is modulated by regional origin of local African ancestry.

DISCUSSION

These analyses identified novel AD-associated loci in individuals of African ancestry and suggest that degree of African ancestry modulates some associations. Increased sample sets covering as much African genetic diversity as possible will be critical to identify additional loci and deconvolute local genetic ancestry effects.

HIGHLIGHTS

Genetic ancestry significantly impacts risk of Alzheimer's Disease (AD). Although individuals of African ancestry are twice as likely to develop AD, they are vastly underrepresented in AD genomics studies. The Alzheimer's Disease Genetics Consortium has previously identified 16 common and rare genetic loci associated with AD in African American individuals. The current analyses significantly expand this effort by increasing the sample size and extending ancestral diversity by including populations from continental Africa. Single variant meta-analysis identified a novel genome-wide significant AD-risk locus in individuals of African ancestry at the MPDZ gene, and 11 additional novel loci with suggestive genome-wide significance at p < 9×10-7. Comparison of African American datasets with samples of higher degree of African ancestry demonstrated differing patterns of association and linkage disequilibrium at one of these loci, suggesting that degree and/or geographic origin of African ancestry modulates the effect at this locus. These findings illustrate the importance of increasing number and ancestral diversity of African ancestry samples in AD genomics studies to fully disentangle the genetic architecture underlying AD, and yield more effective ancestry-informed genetic screening tools and therapeutic interventions.

Puzzle Hi-C: An accurate scaffolding software

High-quality, chromosome-scale genomes are essential for genomic analyses. Analyses, including 3D genomics, epigenetics, and comparative genomics rely on a high-quality genome assembly, which is often accomplished with the assistance of Hi-C data. Curation of genomes reveal that current Hi-C-assisted scaffolding algorithms either generate ordering and orientation errors or fail to assemble high-quality chromosome-level scaffolds. Here, we offer the software Puzzle Hi-C, which uses Hi-C reads to accurately assign contigs or scaffolds to chromosomes. Puzzle Hi-C uses the triangle region instead of the square region to count interactions in a Hi-C heatmap. This strategy dramatically diminishes scaffolding interference caused by long-range interactions. This software also introduces a dynamic, triangle window strategy during assembly. Initially small, the window expands with interactions to produce more effective clustering. Puzzle Hi-C outperforms available scaffolding tools.

Accelerating 3D genomics data analysis with Microcket

The three-dimensional (3D) organization of genome is fundamental to cell biology. To explore 3D genome, emerging high-throughput approaches have produced billions of sequencing reads, which is challenging and time-consuming to analyze. Here we present Microcket, a package for mapping and extracting interacting pairs from 3D genomics data, including Hi-C, Micro-C, and derivant protocols. Microcket utilizes a unique read-stitch strategy that takes advantage of the long read cycles in modern DNA sequencers; benchmark evaluations reveal that Microcket runs much faster than the current tools along with improved mapping efficiency, and thus shows high potential in accelerating and enhancing the biological investigations into 3D genome. Microcket is freely available at https://github.com/hellosunking/Microcket .

STAG2 mutations regulate 3D genome organization, chromatin loops, and Polycomb signaling in glioblastoma multiforme

Inactivating mutations of genes encoding the cohesin complex are common in a wide range of human cancers. STAG2 is the most commonly mutated subunit. Here we report the impact of stable correction of endogenous, naturally occurring STAG2 mutations on gene expression, 3D genome organization, chromatin loops, and Polycomb signaling in glioblastoma multiforme (GBM). In two GBM cell lines, correction of their STAG2 mutations significantly altered the expression of ∼10% of all expressed genes. Virtually all the most highly regulated genes were negatively regulated by STAG2 (i.e., expressed higher in STAG2-mutant cells), and one of them-HEPH-was regulated by STAG2 in uncultured GBM tumors as well. While STAG2 correction had little effect on large-scale features of 3D genome organization (A/B compartments, TADs), STAG2 correction did alter thousands of individual chromatin loops, some of which controlled the expression of adjacent genes. Loops specific to STAG2-mutant cells, which were regulated by STAG1-containing cohesin complexes, were very large, supporting prior findings that STAG1-containing cohesin complexes have greater loop extrusion processivity than STAG2-containing cohesin complexes and suggesting that long loops may be a general feature of STAG2-mutant cancers. Finally, STAG2 mutation activated Polycomb activity leading to increased H3K27me3 marks, identifying Polycomb signaling as a potential target for therapeutic intervention in STAG2-mutant GBM tumors. Together, these findings illuminate the landscape of STAG2-regulated genes, A/B compartments, chromatin loops, and pathways in GBM, providing important clues into the largely still unknown mechanism of STAG2 tumor suppression.

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.

Enhancer selectivity across cell types delineates three functionally distinct enhancer-promoter regulation patterns

BACKGROUND

Multiple enhancers co-regulating the same gene is prevalent and plays a crucial role during development and disease. However, how multiple enhancers coordinate the same gene expression across various cell types remains largely unexplored at genome scale.

RESULTS

We develop a computational approach that enables the quantitative assessment of enhancer specificity and selectivity across diverse cell types, leveraging enhancer-promoter (E-P) interactions data. We observe two well-known gene regulation patterns controlled by enhancer clusters, which regulate the same gene either in a limited number of cell types (Specific pattern, Spe) or in the majority of cell types (Conserved pattern, Con), both of which are enriched for super-enhancers (SEs). We identify a previously overlooked pattern (Variable pattern, Var) that multiple enhancers link to the same gene, but rarely coexist in the same cell type. These three patterns control the genes associating with distinct biological function and exhibit unique epigenetic features. Specifically, we discover a subset of Var patterns contains Shared enhancers with stable enhancer-promoter interactions in the majority of cell types, which might contribute to maintaining gene expression by recruiting abundant CTCF.

CONCLUSIONS

Together, our findings reveal three distinct E-P regulation patterns across different cell types, providing insights into deciphering the complexity of gene transcriptional regulation.

Analysis of long-range chromatin contacts, compartments and looping between mouse embryonic stem cells, lens epithelium and lens fibers

BACKGROUND

Nuclear organization of interphase chromosomes involves individual chromosome territories, "open" and "closed" chromatin compartments, topologically associated domains (TADs) and chromatin loops. The DNA- and RNA-binding transcription factor CTCF together with the cohesin complex serve as major organizers of chromatin architecture. Cellular differentiation is driven by temporally and spatially coordinated gene expression that requires chromatin changes of individual loci of various complexities. Lens differentiation represents an advantageous system to probe transcriptional mechanisms underlying tissue-specific gene expression including high transcriptional outputs of individual crystallin genes until the mature lens fiber cells degrade their nuclei.

RESULTS

Chromatin organization between mouse embryonic stem (ES) cells, newborn (P0.5) lens epithelium and fiber cells were analyzed using Hi-C. Localization of CTCF in both lens chromatins was determined by ChIP-seq and compared with ES cells. Quantitative analyses show major differences between number and size of TADs and chromatin loop size between these three cell types. In depth analyses show similarities between lens samples exemplified by overlaps between compartments A and B. Lens epithelium-specific CTCF peaks are found in mostly methylated genomic regions while lens fiber-specific and shared peaks occur mostly within unmethylated DNA regions. Major differences in TADs and loops are illustrated at the ~ 500 kb Pax6 locus, encoding the critical lens regulatory transcription factor and within a larger ~ 15 Mb WAGR locus, containing Pax6 and other loci linked to human congenital diseases. Lens and ES cell Hi-C data (TADs and loops) together with ATAC-seq, CTCF, H3K27ac, H3K27me3 and ENCODE cis-regulatory sites are shown in detail for the Pax6, Sox1 and Hif1a loci, multiple crystallin genes and other important loci required for lens morphogenesis. The majority of crystallin loci are marked by unexpectedly high CTCF-binding across their transcribed regions.

CONCLUSIONS

Our study has generated the first data on 3-dimensional (3D) nuclear organization in lens epithelium and lens fibers and directly compared these data with ES cells. These findings generate novel insights into lens-specific transcriptional gene control, open new research avenues to study transcriptional condensates in lens fiber cells, and enable studies of non-coding genetic variants linked to cataract and other lens and ocular abnormalities.

Unravelling the genetic basis of Schizophrenia

Neuronal development is a highly regulated mechanism that is central to organismal function in animals. In humans, disruptions to this process can lead to a range of neurodevelopmental phenotypes, including Schizophrenia (SCZ). SCZ has a significant genetic component, whereby an individual with an SCZ affected family member is eight times more likely to develop the disease than someone with no family history of SCZ. By examining a combination of genomic, transcriptomic and epigenomic datasets, large-scale 'omics' studies aim to delineate the relationship between genetic variation and abnormal cellular activity in the SCZ brain. Herein, we provide a brief overview of some of the key omics methods currently being used in SCZ research, including RNA-seq, the assay for transposase-accessible chromatin with high-throughput sequencing (ATAC-seq) and high-throughput chromosome conformation capture (3C) approaches (e.g., Hi-C), as well as single-cell/nuclei iterations of these methods. We also discuss how these techniques are being employed to further our understanding of the genetic basis of SCZ, and to identify associated molecular pathways, biomarkers, and candidate drug targets.

Integrative high-throughput enhancer surveying and functional verification divulges a YY2-condensed regulatory axis conferring risk for osteoporosis

The precise roles of chromatin organization at osteoporosis risk loci remain largely elusive. Here, we combined chromatin interaction conformation (Hi-C) profiling and self-transcribing active regulatory region sequencing (STARR-seq) to qualify enhancer activities of prioritized osteoporosis-associated single-nucleotide polymorphisms (SNPs). We identified 319 SNPs with biased allelic enhancer activity effect (baaSNPs) that linked to hundreds of candidate target genes through chromatin interactions across 146 loci. Functional characterizations revealed active epigenetic enrichment for baaSNPs and prevailing osteoporosis-relevant regulatory roles for their chromatin interaction genes. Further motif enrichment and network mapping prioritized several putative, key transcription factors (TFs) controlling osteoporosis binding to baaSNPs. Specifically, we selected one top-ranked TF and deciphered that an intronic baaSNP (rs11202530) could allele-preferentially bind to YY2 to augment PAPSS2 expression through chromatin interactions and promote osteoblast differentiation. Our results underline the roles of TF-mediated enhancer-promoter contacts for osteoporosis, which may help to better understand the intricate molecular regulatory mechanisms underlying osteoporosis risk loci.

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.

Power-law behavior of transcriptional bursting regulated by enhancer-promoter communication

Revealing how transcriptional bursting kinetics are genomically encoded is challenging because genome structures are stochastic at the organization level and are suggestively linked to gene transcription. To address this challenge, we develop a generic theoretical framework that integrates chromatin dynamics, enhancer-promoter (E-P) communication, and gene-state switching to study transcriptional bursting. The theory predicts that power law can be a general rule to quantitatively describe bursting modulations by E-P spatial communication. Specifically, burst frequency and burst size are up-regulated by E-P communication strength, following power laws with positive exponents. Analysis of the scaling exponents further reveals that burst frequency is preferentially regulated. Bursting kinetics are down-regulated by E-P genomic distance with negative power-law exponents, and this negative modulation desensitizes at large distances. The mutual information between burst frequency (or burst size) and E-P spatial distance further reveals essential characteristics of the information transfer from E-P communication to transcriptional bursting kinetics. These findings, which are in agreement with experimental observations, not only reveal fundamental principles of E-P communication in transcriptional bursting but also are essential for understanding cellular decision-making.

A trade-off in evolution: the adaptive landscape of spiders without venom glands

BACKGROUND

Venom glands play a key role in the predation and defense strategies of almost all spider groups. However, the spider family Uloboridae lacks venom glands and has evolved an adaptive strategy: they excessively wrap their prey directly with spider silk instead of paralyzing it first with toxins. This shift in survival strategy is very fascinating, but the genetic underpinnings behind it are poorly understood.

RESULTS

Spanning multiple spider groups, we conducted multiomics analyses on Octonoba sinensis and described the adaptive evolution of the Uloboridae family at the genome level. We observed the coding genes of myosin and twitchin in muscles are under positive selection, energy metabolism functions are enhanced, and gene families related to tracheal development and tissue mechanical strength are expanded or emerged, all of which are related to the unique anatomical structure and predatory behavior of spiders in the family Uloboridae. In addition, we also scanned the elements that are absent or under relaxed purifying selection, as well as toxin gene homologs in the genomes of 2 species in this family. The results show that the absence of regions and regions under relaxed selection in these spiders' genomes are concentrated in areas related to development and neurosystem. The search for toxin homologs reveals possible gene function shift between toxins and nontoxins and confirms that there are no reliable toxin genes in the genome of this group.

CONCLUSIONS

This study demonstrates the trade-off between different predation strategies in spiders, using either chemical or physical strategy, and provides insights into the possible mechanism underlying this trade-off. Venomless spiders need to mobilize multiple developmental and metabolic pathways related to motor function and limb mechanical strength to cover the decline in adaptability caused by the absence of venom glands.

Mapping enhancer and chromatin accessibility landscapes charts the regulatory network of Alzheimer's disease

BACKGROUND

Enhancers are regulatory elements that target and modulate gene expression and play a role in human health and disease. However, the roles of enhancer regulatory circuit abnormalities driven by epigenetic alterations in Alzheimer's disease (AD) are unclear.

METHODS

In this study, a multiomic integrative analysis was performed to map enhancer and chromatin accessibility landscapes and identify regulatory network abnormalities in AD. We identified differentially methylated enhancers and constructed regulatory networks across brain regions using AD brain tissue samples. Through the integration of snATAC-seq and snRNA-seq datasets, we mapped enhancers with DNA methylation alterations (DMA) and chromatin accessibility landscapes. Core regulatory triplets that contributed to AD neuropathology in specific cell types were further prioritized.

RESULTS

We revealed widespread DNA methylation alterations (DMA) in the enhancers of AD patients across different brain regions. In addition, the genome-wide transcription factor (TF) binding profiles showed that enhancers with DMA are pervasively regulated by TFs. The TF-enhancer-gene regulatory network analysis identified core regulatory triplets that are associated with brain and immune cell proportions and play important roles in AD pathogenesis. Enhancer regulatory circuits with DMA exhibited distinct chromatin accessibility patterns, which were further characterized at single-cell resolutions.

CONCLUSIONS

Our study comprehensively investigated DNA methylation-mediated regulatory circuit abnormalities and provided novel insights into the potential pathogenesis of AD.

Genetic analyses in multiplex families confirms chromosome 5q35 as a risk locus for Alzheimer's Disease in individuals of African Ancestry

There is a paucity of genetic studies of Alzheimer Disease (AD) in individuals of African Ancestry, despite evidence suggesting increased risk of AD in the African American (AA) population. We performed whole-genome sequencing (WGS) and multipoint linkage analyses in 51 multi-generational AA AD families ascertained through the Research in African American Alzheimer Disease Initiative (REAAADI) and the National Institute on Aging Late Onset Alzheimer's disease (NIA-LOAD) Family Based Study. Variants were prioritized on minor allele frequency (<0.01), functional potential of coding and noncoding variants, co-segregation with AD and presence in multi-ancestry ADSP release 3 WGS data. We identified a significant linkage signal on chromosome 5q35 (HLOD=3.3) driven by nine families. Haplotype segregation analysis in the family with highest LOD score identified a 3'UTR variant in INSYN2B with the most functional evidence. Four other linked AA families harbor within-family shared variants located in INSYN2B's promoter or enhancer regions. This AA family-based finding shows the importance of diversifying population-level genetic data to better understand the genetic determinants of AD on a global scale.

Cell-type-specific 3D-genome organization and transcription regulation in the brain

3D organization of the genome plays a critical role in regulating gene expression. However, it remains unclear how chromatin organization differs among different cell types in the brain. Here we used genome-scale DNA and RNA imaging to investigate 3D-genome organization in transcriptionally distinct cell types in the primary motor cortex of the mouse brain. We uncovered a wide spectrum of differences in the nuclear architecture and 3D-genome organization among different cell types, ranging from the physical size of the cell nucleus to the active-inactive chromatin compartmentalization and radial positioning of chromatin loci within the nucleus. These cell-type-dependent variations in nuclear architecture and chromatin organization exhibited strong correlation with both total transcriptional activity of the cell and transcriptional regulation of cell-type-specific marker genes. Moreover, we found that the methylated-DNA-binding protein MeCP2 regulates transcription in a divergent manner, depending on the nuclear radial positions of chromatin loci, through modulating active-inactive chromatin compartmentalization.

Joint epigenome profiling reveals cell-type-specific gene regulatory programmes in human cortical organoids

Gene expression is regulated by multiple epigenetic mechanisms, which are coordinated in development and disease. However, current multiomics methods are frequently limited to one or two modalities at a time, making it challenging to obtain a comprehensive gene regulatory signature. Here, we describe a method-3D genome, RNA, accessibility and methylation sequencing (3DRAM-seq)-that simultaneously interrogates spatial genome organization, chromatin accessibility and DNA methylation genome-wide and at high resolution. We combine 3DRAM-seq with immunoFACS and RNA sequencing in cortical organoids to map the cell-type-specific regulatory landscape of human neural development across multiple epigenetic layers. Finally, we apply a massively parallel reporter assay to profile cell-type-specific enhancer activity in organoids and to functionally assess the role of key transcription factors for human enhancer activation and function. More broadly, 3DRAM-seq can be used to profile the multimodal epigenetic landscape in rare cell types and different tissues.

4D nucleome equation predicts gene expression controlled by long-range enhancer-promoter interaction

Recent experimental evidence strongly supports that three-dimensional (3D) long-range enhancer-promoter (E-P) interactions have important influences on gene-expression dynamics, but it is unclear how the interaction information is translated into gene expression over time (4D). To address this question, we developed a general theoretical framework (named as a 4D nucleome equation), which integrates E-P interactions on chromatin and biochemical reactions of gene transcription. With this equation, we first present the distribution of mRNA counts as a function of the E-P genomic distance and then reveal a power-law scaling of the expression level in this distance. Interestingly, we find that long-range E-P interactions can induce bimodal and trimodal mRNA distributions. The 4D nucleome equation also allows for model selection and parameter inference. When this equation is applied to the mouse embryonic stem cell smRNA-FISH data and the E-P genomic-distance data, the predicted E-P contact probability and mRNA distribution are in good agreement with experimental results. Further statistical inference indicates that the E-P interactions prefer to modulate the mRNA level by controlling promoter activation and transcription initiation rates. Our model and results provide quantitative insights into both spatiotemporal gene-expression determinants (i.e., long-range E-P interactions) and cellular fates during development.

Leveraging iPSC technology to assess neuro-immune interactions in neurological and psychiatric disorders

Communication between the immune and the nervous system is essential for human brain development and homeostasis. Disruption of this intricately regulated crosstalk can lead to neurodevelopmental, psychiatric, or neurodegenerative disorders. While animal models have been essential in characterizing the role of neuroimmunity in development and disease, they come with inherent limitations due to species specific differences, particularly with regard to microglia, the major subset of brain resident immune cells. The advent of induced pluripotent stem cell (iPSC) technology now allows the development of clinically relevant models of the central nervous system that adequately reflect human genetic architecture. This article will review recent publications that have leveraged iPSC technology to assess neuro-immune interactions. First, we will discuss the role of environmental stressors such as neurotropic viruses or pro-inflammatory cytokines on neuronal and glial function. Next, we will review how iPSC models can be used to study genetic risk factors in neurological and psychiatric disorders. Lastly, we will evaluate current challenges and future potential for iPSC models in the field of neuroimmunity.

Genome-Wide Association Study of Pericardial Fat Area in 28 161 UK Biobank Participants

BACKGROUND Pericardial adipose tissue (PAT) is the visceral adipose tissue compartment surrounding the heart. Experimental and observational research has suggested that greater PAT deposition might mediate cardiovascular disease, independent of general or subcutaneous adiposity. We characterize the genetic architecture of adiposity-adjusted PAT and identify causal associations between PAT and adverse cardiac magnetic resonance imaging measures of cardiac structure and function in 28 161 UK Biobank participants. METHODS AND RESULTS The PAT phenotype was extracted from cardiac magnetic resonance images using an automated image analysis tool previously developed and validated in this cohort. A genome-wide association study was performed with PAT area set as the phenotype, adjusting for age, sex, and other measures of obesity. Functional mapping and Bayesian colocalization were used to understand the biologic role of identified variants. Mendelian randomization analysis was used to examine potential causal links between genetically determined PAT and cardiac magnetic resonance-derived measures of left ventricular structure and function. We discovered 12 genome-wide significant variants, with 2 independent sentinel variants (rs6428792, P=4.20×10-9 and rs11992444, P=1.30×10-12) at 2 distinct genomic loci, that were mapped to 3 potentially causal genes: T-box transcription factor 15 (TBX15), tryptophanyl tRNA synthetase 2, mitochondrial (WARS2) and early B-cell factor-2 (EBF2) through functional annotation. Bayesian colocalization additionally suggested a role of RP4-712E4.1. Genetically predicted differences in adiposity-adjusted PAT were causally associated with adverse left ventricular remodeling. CONCLUSIONS This study provides insights into the genetic architecture determining differential PAT deposition, identifies causal links with left structural and functional parameters, and provides novel data about the pathophysiological importance of adiposity distribution.

Genome-wide analysis of the interplay between chromatin-associated RNA and 3D genome organization in human cells

The interphase genome is dynamically organized in the nucleus and decorated with chromatin-associated RNA (caRNA). It remains unclear whether the genome architecture modulates the spatial distribution of caRNA and vice versa. Here, we generate a resource of genome-wide RNA-DNA and DNA-DNA contact maps in human cells. These maps reveal the chromosomal domains demarcated by locally transcribed RNA, hereafter termed RNA-defined chromosomal domains. Further, the spreading of caRNA is constrained by the boundaries of topologically associating domains (TADs), demonstrating the role of the 3D genome structure in modulating the spatial distribution of RNA. Conversely, stopping transcription or acute depletion of RNA induces thousands of chromatin loops genome-wide. Activation or suppression of the transcription of specific genes suppresses or creates chromatin loops straddling these genes. Deletion of a specific caRNA-producing genomic sequence promotes chromatin loops that straddle the interchromosomal target sequences of this caRNA. These data suggest a feedback loop where the 3D genome modulates the spatial distribution of RNA, which in turn affects the dynamic 3D genome organization.

HiCAR: Analysis of Open Chromatin Associated Long-range Chromatin Interaction Using Low-Input Materials

Cis-regulatory elements (cREs) and their long-range interactions are crucial for spatial-temporal gene regulation. While cREs can be characterized as accessible chromatin sequences, comprehensively identifying their spatial interactions remains a challenge. We recently developed a method, HiCAR (Hi-C on Accessible Regulatory DNA), which combines Tn5 transposase and chromatin proximity ligation to analyze open chromatin-anchored interactions in low-input cells. Application of HiCAR in human embryonic stem cells and lymphoblastoid cells reveals high-resolution chromatin contacts with efficiency comparable to in situ Hi-C across various distance ranges. Moreover, HiCAR was successfully applied to 30,000 primary human muscle stem cells, showcasing its potential for analyzing chromatin accessibility and looping in low-input primary cells and clinical samples. Here, we provide a detailed step-by-step protocol to perform the updated HiCAR experiments. © 2023 Wiley Periodicals LLC. Basic Protocol 1: Tn5 Transposase Assembly Basic Protocol 2: HiCAR Library Preparation.

Single cell multiomic analysis reveals diabetes-associated β-cell heterogeneity driven by HNF1A

Broad heterogeneity in pancreatic β-cell function and morphology has been widely reported. However, determining which components of this cellular heterogeneity serve a diabetes-relevant function remains challenging. Here, we integrate single-cell transcriptome, single-nuclei chromatin accessibility, and cell-type specific 3D genome profiles from human islets and identify Type II Diabetes (T2D)-associated β-cell heterogeneity at both transcriptomic and epigenomic levels. We develop a computational method to explicitly dissect the intra-donor and inter-donor heterogeneity between single β-cells, which reflect distinct mechanisms of T2D pathogenesis. Integrative transcriptomic and epigenomic analysis identifies HNF1A as a principal driver of intra-donor heterogeneity between β-cells from the same donors; HNF1A expression is also reduced in β-cells from T2D donors. Interestingly, HNF1A activity in single β-cells is significantly associated with lower Na+ currents and we nominate a HNF1A target, FXYD2, as the primary mitigator. Our study demonstrates the value of investigating disease-associated single-cell heterogeneity and provides new insights into the pathogenesis of T2D.

Three-dimensional genome structure and function

Linear DNA undergoes a series of compression and folding events, forming various three-dimensional (3D) structural units in mammalian cells, including chromosomal territory, compartment, topologically associating domain, and chromatin loop. These structures play crucial roles in regulating gene expression, cell differentiation, and disease progression. Deciphering the principles underlying 3D genome folding and the molecular mechanisms governing cell fate determination remains a challenge. With advancements in high-throughput sequencing and imaging techniques, the hierarchical organization and functional roles of higher-order chromatin structures have been gradually illuminated. This review systematically discussed the structural hierarchy of the 3D genome, the effects and mechanisms of cis-regulatory elements interaction in the 3D genome for regulating spatiotemporally specific gene expression, the roles and mechanisms of dynamic changes in 3D chromatin conformation during embryonic development, and the pathological mechanisms of diseases such as congenital developmental abnormalities and cancer, which are attributed to alterations in 3D genome organization and aberrations in key structural proteins. Finally, prospects were made for the research about 3D genome structure, function, and genetic intervention, and the roles in disease development, prevention, and treatment, which may offer some clues for precise diagnosis and treatment of related diseases.

Breaking enhancers to gain insights into developmental defects

Despite ground-breaking genetic studies that have identified thousands of risk variants for developmental diseases, how these variants lead to molecular and cellular phenotypes remains a gap in knowledge. Many of these variants are non-coding and occur at enhancers, which orchestrate key regulatory programs during development. The prevailing paradigm is that non-coding variants alter the activity of enhancers, impacting gene expression programs, and ultimately contributing to disease risk. A key obstacle to progress is the systematic functional characterization of non-coding variants at scale, especially since enhancer activity is highly specific to cell type and developmental stage. Here, we review the foundational studies of enhancers in developmental disease and current genomic approaches to functionally characterize developmental enhancers and their variants at scale. In the coming decade, we anticipate systematic enhancer perturbation studies to link non-coding variants to molecular mechanisms, changes in cell state, and disease phenotypes.

Improving comparative analyses of Hi-C data via contrastive self-supervised learning

Hi-C is a widely applied chromosome conformation capture (3C)-based technique, which has produced a large number of genomic contact maps with high sequencing depths for a wide range of cell types, enabling comprehensive analyses of the relationships between biological functionalities (e.g. gene regulation and expression) and the three-dimensional genome structure. Comparative analyses play significant roles in Hi-C data studies, which are designed to make comparisons between Hi-C contact maps, thus evaluating the consistency of replicate Hi-C experiments (i.e. reproducibility measurement) and detecting statistically differential interacting regions with biological significance (i.e. differential chromatin interaction detection). However, due to the complex and hierarchical nature of Hi-C contact maps, it remains challenging to conduct systematic and reliable comparative analyses of Hi-C data. Here, we proposed sslHiC, a contrastive self-supervised representation learning framework, for precisely modeling the multi-level features of chromosome conformation and automatically producing informative feature embeddings for genomic loci and their interactions to facilitate comparative analyses of Hi-C contact maps. Comprehensive computational experiments on both simulated and real datasets demonstrated that our method consistently outperformed the state-of-the-art baseline methods in providing reliable measurements of reproducibility and detecting differential interactions with biological meanings.

Setting the clock of neural progenitor cells during mammalian corticogenesis

Radial glial cells (RGCs) as primary neural stem cells in the developing mammalian cortex give rise to diverse types of neurons and glial cells according to sophisticated developmental programs with remarkable spatiotemporal precision. Recent studies suggest that regulation of the temporal competence of RGCs is a key mechanism for the highly conserved and predictable development of the cerebral cortex. Various types of epigenetic regulations, such as DNA methylation, histone modifications, and 3D chromatin architecture, play a key role in shaping the gene expression pattern of RGCs. In addition, epitranscriptomic modifications regulate temporal pre-patterning of RGCs by affecting the turnover rate and function of cell-type-specific transcripts. In this review, we summarize epigenetic and epitranscriptomic regulatory mechanisms that control the temporal competence of RGCs during mammalian corticogenesis. Furthermore, we discuss various developmental elements that also dynamically regulate the temporal competence of RGCs, including biochemical reaction speed, local environmental changes, and subcellular organelle remodeling. Finally, we discuss the underlying mechanisms that regulate the interspecies developmental tempo contributing to human-specific features of brain development.

Characterization of Gene Regulatory Elements in Human Fetal Cortical Development: Enhancing Our Understanding of Neurodevelopmental Disorders and Evolution

The neocortex is the region that most distinguishes human brain from other mammals and primates [Annu Rev Genet. 2021 Nov;55(1):555-81]. Studying the development of human cortex is important in understanding the evolutionary changes occurring in humans relative to other primates, as well as in elucidating mechanisms underlying neurodevelopmental disorders. Cortical development is a highly regulated process, spatially and temporally coordinated by expression of essential transcriptional factors in response to signaling pathways [Neuron. 2019 Sep;103(6):980-1004]. Enhancers are the most well-understood cis-acting, non-protein-coding regulatory elements that regulate gene expression [Nat Rev Genet. 2014 Apr;15(4):272-86]. Importantly, given the conservation of both DNA sequence and molecular function of the majority of proteins across mammals [Genome Res. 2003 Dec;13(12):2507-18], enhancers [Science. 2015 Mar;347(6226):1155-9], which are far more divergent at the sequence level, likely account for the phenotypes that distinguish the human brain by changing the regulation of gene expression. In this review, we will revisit the conceptual framework of gene regulation during human brain development, as well as the evolution of technologies to study transcriptional regulation, with recent advances in genome biology that open a window allowing us to systematically characterize cis-regulatory elements in developing human brain [Hum Mol Genet. 2022 Oct;31(R1):R84-96]. We provide an update on work to characterize the suite of all enhancers in the developing human brain and the implications for understanding neuropsychiatric disorders. Finally, we discuss emerging therapeutic ideas that utilize our emerging knowledge of enhancer function.

Novel correlative analysis identifies multiple genomic variations impacting ASD with macrocephaly

Autism spectrum disorders (ASD) display both phenotypic and genetic heterogeneity, impeding the understanding of ASD and development of effective means of diagnosis and potential treatments. Genes affected by genomic variations for ASD converge in dozens of gene ontologies (GOs), but the relationship between the variations at the GO level have not been well elucidated. In the current study, multiple types of genomic variations were mapped to GOs and correlations among GOs were measured in ASD and control samples. Several ASD-unique GO correlations were found, suggesting the importance of co-occurrence of genomic variations in genes from different functional categories in ASD etiology. Combined with experimental data, several variations related to WNT signaling, neuron development, synapse morphology/function and organ morphogenesis were found to be important for ASD with macrocephaly, and novel co-occurrence patterns of them in ASD patients were found. Furthermore, we applied this gene ontology correlation analysis method to find genomic variations that contribute to ASD etiology in combination with changes in gene expression and transcription factor binding, providing novel insights into ASD with macrocephaly and a new methodology for the analysis of genomic variation.

Considerations and caveats for analyzing chromatin compartments

Genomes are organized into nuclear compartments, separating active from inactive chromatin. Chromatin compartments are readily visible in a large number of species by experiments that map chromatin conformation genome-wide. When analyzing these maps, a common step is the identification of genomic intervals that interact within A (active) and B (inactive) compartments. It has also become increasingly common to identify and analyze subcompartments. We review different strategies to identify A/B and subcompartment intervals, including a discussion of various machine-learning approaches to predict these features. We then discuss the strengths and limitations of current strategies and examine how these aspects of analysis may have impacted our understanding of chromatin compartments.

Development and Application of Brain Region-Specific Organoids for Investigating Psychiatric Disorders

Human society has been burdened by psychiatric disorders throughout the course of its history. The emergence and rapid advances of human brain organoid technology provide unprecedented opportunities for investigation of potential disease mechanisms and development of targeted or even personalized treatments for various psychiatric disorders. In this review, we summarize recent advances for generating organoids from human pluripotent stem cells to model distinct brain regions and diverse cell types. We also highlight recent progress, discuss limitations, and propose potential improvements in using patient-derived or genetically engineered brain region-specific organoids for investigating various psychiatric disorders.

MYT1L is required for suppressing earlier neuronal development programs in the adult mouse brain

In vitro studies indicate the neurodevelopmental disorder gene myelin transcription factor 1-like (MYT1L) suppresses non-neuronal lineage genes during fibroblast-to-neuron direct differentiation. However, MYT1L's molecular and cellular functions in the adult mammalian brain have not been fully characterized. Here, we found that MYT1L loss leads to up-regulated deep layer (DL) gene expression, corresponding to an increased ratio of DL/UL neurons in the adult mouse cortex. To define potential mechanisms, we conducted Cleavage Under Targets & Release Using Nuclease (CUT&RUN) to map MYT1L binding targets and epigenetic changes following MYT1L loss in mouse developing cortex and adult prefrontal cortex (PFC). We found MYT1L mainly binds to open chromatin, but with different transcription factor co-occupancies between promoters and enhancers. Likewise, multiomic data set integration revealed that, at promoters, MYT1L loss does not change chromatin accessibility but increases H3K4me3 and H3K27ac, activating both a subset of earlier neuronal development genes as well as Bcl11b, a key regulator for DL neuron development. Meanwhile, we discovered that MYT1L normally represses the activity of neurogenic enhancers associated with neuronal migration and neuronal projection development by closing chromatin structures and promoting removal of active histone marks. Further, we showed that MYT1L interacts with HDAC2 and transcriptional repressor SIN3B in vivo, providing potential mechanisms underlying repressive effects on histone acetylation and gene expression. Overall, our findings provide a comprehensive map of MYT1L binding in vivo and mechanistic insights into how MYT1L loss leads to aberrant activation of earlier neuronal development programs in the adult mouse brain.

HiConfidence: a novel approach uncovering the biological signal in Hi-C data affected by technical biases

The chromatin interaction assays, particularly Hi-C, enable detailed studies of genome architecture in multiple organisms and model systems, resulting in a deeper understanding of gene expression regulation mechanisms mediated by epigenetics. However, the analysis and interpretation of Hi-C data remain challenging due to technical biases, limiting direct comparisons of datasets obtained in different experiments and laboratories. As a result, removing biases from Hi-C-generated chromatin contact matrices is a critical data analysis step. Our novel approach, HiConfidence, eliminates biases from the Hi-C data by weighing chromatin contacts according to their consistency between replicates so that low-quality replicates do not substantially influence the result. The algorithm is effective for the analysis of global changes in chromatin structures such as compartments and topologically associating domains. We apply the HiConfidence approach to several Hi-C datasets with significant technical biases, that could not be analyzed effectively using existing methods, and obtain meaningful biological conclusions. In particular, HiConfidence aids in the study of how changes in histone acetylation pattern affect chromatin organization in Drosophila melanogaster S2 cells. The method is freely available at GitHub: https://github.com/victorykobets/HiConfidence.

Chromatin modules and their implication in genomic organization and gene regulation

Regulation of gene expression is a complex but highly guided process. While genomic technologies and computational approaches have allowed high-throughput mapping of cis-regulatory elements (CREs) and their interactions in 3D, their precise role in regulating gene expression remains obscure. Recent complementary observations revealed that interactions between CREs frequently result in the formation of small-scale functional modules within topologically associating domains. Such chromatin modules likely emerge from a complex interplay between regulatory machineries assembled at CREs, including site-specific binding of transcription factors. Here, we review the methods that allow identifying chromatin modules, summarize possible mechanisms that steer CRE interactions within these modules, and discuss outstanding challenges to uncover how chromatin modules fit in our current understanding of the functional 3D genome.

RNAcontacts: A Pipeline for Predicting Contacts from RNA Proximity Ligation Assays

High-throughput RNA proximity ligation assays are molecular methods that are used to simultaneously analyze the spatial proximity of many RNAs in living cells. Their principle is based on cross-linking, fragmentation, and subsequent religation of RNAs, followed by high-throughput sequencing. The generated fragments have two different types of splits, one resulting from pre-mRNA splicing and the other formed by the ligation of spatially close RNA strands. Here, we present RNAcontacts, a universal pipeline for detecting RNA-RNA contacts in high-throughput RNA proximity ligation assays. RNAcontacts circumvents the inherent problem of mapping sequences with two distinct types of splits using a two-pass alignment, in which splice junctions are inferred from a control RNA-seq experiment on the first pass and then provided to the aligner as bona fide introns on the second pass. Compared to previously developed methods, our approach allows for a more sensitive detection of RNA contacts and has a higher specificity with respect to splice junctions that are present in the biological sample. RNAcontacts automatically extracts contacts, clusters their ligation points, computes the read support, and generates tracks for visualizing through the UCSC Genome Browser. The pipeline is implemented in Snakemake, a reproducible and scalable workflow management system for rapid and uniform processing of multiple datasets. RNAcontacts is a generic pipeline for the detection of RNA contacts that can be used with any proximity ligation method as long as one of the interacting partners is RNA. RNAcontacts is available via the GitHub repository https://github.com/smargasyuk/ RNAcontacts/.

Easy Hi-C: A Low-Input Method for Capturing Genome Organization

Chromosome conformation capture technology and its derivatives have been widely used to study genome organization. Among them, Hi-C (chromosome conformation capture coupling with high-throughput sequencing) is popular in dissecting chromatin architecture on the genome-wide level. However, the intrinsic limitations prevent its application when it comes to rare samples. Here, we present easy Hi-C, a biotin-free technology that dramatically reduces DNA loss and is suitable for low-input samples.