Mapping of long-range chromatin interactions by proximity ligation-assisted ChIP-seq

Identifying genetic variants associated with chromatin looping and genome function

Here we present a comprehensive HiChIP dataset on naïve CD4 T cells (nCD4) from 30 donors and identify QTLs that associate with genotype-dependent and/or allele-specific variation of HiChIP contacts defining loops between active regulatory regions (iQTLs). We observe a substantial overlap between iQTLs and previously defined eQTLs and histone QTLs, and an enrichment for fine-mapped QTLs and GWAS variants. Furthermore, we describe a distinct subset of nCD4 iQTLs, for which the significant variation of chromatin contacts in nCD4 are translated into significant eQTL trends in CD4 T cell memory subsets. Finally, we define connectivity-QTLs as iQTLs that are significantly associated with concordant genotype-dependent changes in chromatin contacts over a broad genomic region (e.g., GWAS SNP in the RNASET2 locus). Our results demonstrate the importance of chromatin contacts as a complementary modality for QTL mapping and their power in identifying previously uncharacterized QTLs linked to cell-specific gene expression and connectivity.

Improved simultaneous mapping of epigenetic features and 3D chromatin structure via ViCAR

Methods to measure chromatin contacts at genomic regions bound by histone modifications or proteins are important tools to investigate chromatin organization. However, such methods do not capture the possible involvement of other epigenomic features such as G-quadruplex DNA secondary structures (G4s). To bridge this gap, we introduce ViCAR (viewpoint HiCAR), for the direct antibody-based capture of chromatin interactions at folded G4s. Through ViCAR, we showcase the first G4-3D interaction landscape. Using histone marks, we also demonstrate how ViCAR improves on earlier approaches yielding increased signal-to-noise. ViCAR is a practical and powerful tool to explore epigenetic marks and 3D genome interactomes.

Variant-to-function mapping of late-onset Alzheimer's disease GWAS signals in human microglial cell models implicates RTFDC1 at the CASS4 locus

Late-onset Alzheimer's disease (LOAD) research has principally focused on neurons over the years due to their known role in the production of amyloid beta plaques and neurofibrillary tangles. In contrast, recent genomic studies of LOAD have implicated microglia as culprits of the prolonged inflammation exacerbating the neurodegeneration observed in patient brains. Indeed, recent LOAD genome-wide association studies (GWAS) have reported multiple loci near genes related to microglial function, including TREM2, ABI3, and CR1. However, GWAS alone cannot pinpoint underlying causal variants or effector genes at such loci, as most signals reside in non-coding regions of the genome and could presumably confer their influence frequently via long-range regulatory interactions. We elected to carry out a combination of ATAC-seq and high-resolution promoter-focused Capture-C in two human microglial cell models (iPSC-derived microglia and HMC3) in order to physically map interactions between LOAD GWAS-implicated candidate causal variants and their corresponding putative effector genes. Notably, we observed consistent evidence that rs6024870 at the GWAS CASS4 locus contacted the promoter of nearby gene, RTFDC1. We subsequently observed a directionallly consistent decrease in RTFDC1 expression with the the protective minor A allele of rs6024870 via both luciferase assays in HMC3 cells and expression studies in primary human microglia. Through CRISPR-Cas9-mediated deletion of the putative regulatory region harboring rs6024870 in HMC3 cells, we observed increased pro-inflammatory cytokine secretion and decreased DNA double strand break repair related, at least in part, to RTFDC1 expression levels. Our variant-to-function approach therefore reveals that the rs6024870-harboring regulatory element at the LOAD 'CASS4' GWAS locus influences both microglial inflammatory capacity and DNA damage resolution, along with cumulative evidence implicating RTFDC1 as a novel candidate effector gene.

Beyond the Chromosome: Recent Developments in Decoding the Significance of Extrachromosomal Circular DNA (eccDNA) in Human Malignancies

Extrachromosomal circular DNA (eccDNA) is a form of a circular double-stranded DNA that exists independently of conventional chromosomes. eccDNA exhibits a broad and random distribution across eukaryotic cells and has been associated with tumor-related properties due to its ability to harbor the complete gene information of oncogenes. The complex and multifaceted mechanisms underlying eccDNA formation include pathways such as DNA damage repair, breakage-fusion-bridge (BFB) mechanisms, chromothripsis, and cell apoptosis. Of note, eccDNA plays a pivotal role in tumor development, genetic heterogeneity, and therapeutic resistance. The high copy number and transcriptional activity of oncogenes carried by eccDNA contribute to the accelerated growth of tumors. Notably, the amplification of oncogenes on eccDNA is implicated in the malignant progression of cancer cells. The improvement of high-throughput sequencing techniques has greatly enhanced our knowledge of eccDNA by allowing for a detailed examination of its genetic structures and functions. However, we still lack a comprehensive and efficient annotation for eccDNA, while challenges persist in the study and understanding of the functional role of eccDNA, emphasizing the need for the development of robust methodologies. The potential clinical applications of eccDNA, such as its role as a measurable biomarker or therapeutic target in diseases, particularly within the spectrum of human malignancies, is a promising field for future research. In conclusion, eccDNA represents a quite dynamic and multifunctional genetic entity with far-reaching implications in cancer pathogenesis and beyond. Further research is essential to unravel the molecular pathways of eccDNA formation, elucidate its functional roles, and explore its clinical applications. Addressing these aspects is crucial for advancing our understanding of genomic instability and developing novel strategies for tailored therapeutics, especially in cancer.

3D genomic analysis reveals novel enhancer-hijacking caused by complex structural alterations that drive oncogene overexpression

Cancer genomes are composed of many complex structural alterations on chromosomes and extrachromosomal DNA (ecDNA), making it difficult to identify non-coding enhancer regions that are hijacked to activate oncogene expression. Here, we describe a 3D genomics-based analysis called HAPI (Highly Active Promoter Interactions) to characterize enhancer hijacking. HAPI analysis of HiChIP data from 34 cancer cell lines identified enhancer hijacking events that activate both known and potentially novel oncogenes such as MYC, CCND1, ETV1, CRKL, and ID4. Furthermore, we found enhancer hijacking among multiple oncogenes from different chromosomes, often including MYC, on the same complex amplicons such as ecDNA. We characterized a MYC-ERBB2 chimeric ecDNA, in which ERBB2 heavily hijacks MYC's enhancers. Notably, CRISPRi of the MYC promoter led to increased interaction of ERBB2 with MYC enhancers and elevated ERBB2 expression. Our HAPI analysis tool provides a robust strategy to detect enhancer hijacking and reveals novel insights into oncogene activation.

Genomic Convergence in Hibernating Mammals Elucidates the Genetics of Metabolic Regulation in the Hypothalamus

Elucidating the genetic basis of mammalian metabolism could help define mechanisms central to health and disease. Here, we define conserved cis-regulatory elements (CREs) and programs for mammalian metabolic control. We delineate gene expression and chromatin responses in the mouse hypothalamus for 7 steps of the Fed-to-Fasted-to-Refed (FFR) response process. Comparative genomics of hibernating versus non-hibernating lineages then illuminates cis-elements showing convergent changes in hibernators. Hibernators accumulated loss-of-function effects for specific CREs regulating hypothalamic FFR responses. Multi-omics approaches pinpoint key CREs, genes, regulatory programs, and cell types in the divergence of hibernating and homeothermic lineages. The refeeding period after extended fasting is revealed as one critical period of chromatin remodeling with convergent genomic changes. This genetic framework is a step toward harnessing hibernator adaptations in medicine.

3D genomic analysis reveals novel enhancer-hijacking caused by complex structural alterations that drive oncogene overexpression

Cancer genomes are composed of many complex structural alterations on chromosomes and extrachromosomal DNA (ecDNA), making it difficult to identify non-coding enhancer regions that are hijacked to activate oncogene expression. Here, we describe a 3D genomics-based analysis called HAPI (Highly Active Promoter Interactions) to characterize enhancer hijacking. HAPI analysis of HiChIP data from 34 cancer cell lines identified enhancer hijacking events that activate both known and potentially novel oncogenes such as MYC, CCND1 , ETV1 , CRKL , and ID4 . Furthermore, we found enhancer hijacking among multiple oncogenes from different chromosomes, often including MYC , on the same complex amplicons such as ecDNA. We characterized a MYC - ERBB2 chimeric ecDNA, in which ERBB2 heavily hijacks MYC 's enhancers. Notably, CRISPRi of the MYC promoter led to increased interaction of ERBB2 with MYC enhancers and elevated ERBB2 expression. Our HAPI analysis tool provides a robust strategy to detect enhancer hijacking and reveals novel insights into oncogene activation.

Identification of drug responsive enhancers by predicting chromatin accessibility change from perturbed gene expression profiles

Individual may response to drug treatment differently due to their genetic variants located in enhancers. These variants can alter transcription factor's (TF) binding strength, affect enhancer's chromatin activity or interaction, and eventually change expression level of downstream gene. Here, we propose a computational framework, PERD, to Predict the Enhancers Responsive to Drug. A machine learning model was trained to predict the genome-wide chromatin accessibility from transcriptome data using the paired expression and chromatin accessibility data collected from ENCODE and ROADMAP. Then the model was applied to the perturbed gene expression data from Connectivity Map (CMAP) and Cancer Drug-induced gene expression Signature DataBase (CDS-DB) and identify drug responsive enhancers with significantly altered chromatin accessibility. Furthermore, the drug responsive enhancers were related to the pharmacogenomics genome-wide association studies (PGx GWAS). Stepping on the traditional drug-associated gene signatures, PERD holds the promise to enhance the causality of drug perturbation by providing candidate regulatory element of those drug associated genes.

Regulation of Myc transcription by an enhancer cluster dedicated to pluripotency and early embryonic expression

MYC plays various roles in pluripotent stem cells, including the promotion of somatic cell reprogramming to pluripotency, the regulation of cell competition and the control of embryonic diapause. However, how Myc expression is regulated in this context remains unknown. The Myc gene lies within a ~ 3-megabase gene desert with multiple cis-regulatory elements. Here we use genomic rearrangements, transgenesis and targeted mutation to analyse Myc regulation in early mouse embryos and pluripotent stem cells. We identify a topologically-associated region that homes enhancers dedicated to Myc transcriptional regulation in stem cells of the pre-implantation and early post-implantation embryo. Within this region, we identify elements exclusively dedicated to Myc regulation in pluripotent cells, with distinct enhancers that sequentially activate during naive and formative pluripotency. Deletion of pluripotency-specific enhancers dampens embryonic stem cell competitive ability. These results identify a topologically defined enhancer cluster dedicated to early embryonic expression and uncover a modular mechanism for the regulation of Myc expression in different states of pluripotency.

H3K4me1 facilitates promoter-enhancer interactions and gene activation during embryonic stem cell differentiation

Histone H3 lysine 4 mono-methylation (H3K4me1) marks poised or active enhancers. KMT2C (MLL3) and KMT2D (MLL4) catalyze H3K4me1, but their histone methyltransferase activities are largely dispensable for transcription during early embryogenesis in mammals. To better understand the role of H3K4me1 in enhancer function, we analyze dynamic enhancer-promoter (E-P) interactions and gene expression during neural differentiation of the mouse embryonic stem cells. We found that KMT2C/D catalytic activities were only required for H3K4me1 and E-P contacts at a subset of candidate enhancers, induced upon neural differentiation. By contrast, a majority of enhancers retained H3K4me1 in KMT2C/D catalytic mutant cells. Surprisingly, H3K4me1 signals at these KMT2C/D-independent sites were reduced after acute depletion of KMT2B, resulting in aggravated transcriptional defects. Our observations therefore implicate KMT2B in the catalysis of H3K4me1 at enhancers and provide additional support for an active role of H3K4me1 in enhancer-promoter interactions and transcription in mammalian cells.

Loop Catalog: a comprehensive HiChIP database of human and mouse samples

HiChIP enables cost-effective and high-resolution profiling of regulatory and structural loops. To leverage the increasing number of publicly available HiChIP datasets from diverse cell lines and primary cells, we developed the Loop Catalog (https://loopcatalog.lji.org), a web-based database featuring HiChIP loop calls for 1319 samples across 133 studies and 44 high-resolution Hi-C loop calls. We demonstrate its utility in interpreting fine-mapped GWAS variants (SNP-to-gene linking), in identifying enriched sequence motifs and motif pairs at loop anchors, and in network-level analysis of loops connecting regulatory elements (community detection). Our comprehensive catalog, spanning over 4M unique 5kb loops, along with the accompanying analysis modalities constitutes an important resource for studies in gene regulation and genome organization.

Predicting chromatin conformation contact maps

Over the past 15 years, a variety of next-generation sequencing assays have been developed for measuring the 3D conformation of DNA in the nucleus. Each of these assays gives, for a particular cell or tissue type, a distinct picture of 3D chromatin architecture. Accordingly, making sense of the relationship between genome structure and function requires teasing apart two closely related questions: how does chromatin 3D structure change from one cell type to the next, and how do different measurements of that structure differ from one another, even when the two assays are carried out in the same cell type? In this work, we assemble a collection of chromatin 3D datasets-each represented as a 2D contact map- spanning multiple assay types and cell types. We then build a machine learning model that predicts missing contact maps in this collection. We use the model to systematically explore how genome 3D architecture changes, at the level of compartments, domains, and loops, between cell type and between assay types.

Advances in the multimodal analysis of the 3D chromatin structure and gene regulation

Recent studies have demonstrated that the three-dimensional conformation of the chromatin plays a crucial role in gene regulation, with aberrations potentially leading to various diseases. Advanced methodologies have revealed a link between the chromatin conformation and biological function. This review divides these methodologies into sequencing-based and imaging-based methodologies, tracing their development over time. We particularly highlight innovative techniques that facilitate the simultaneous mapping of RNAs, histone modifications, and proteins within the context of the 3D architecture of chromatin. This multimodal integration substantially improves our ability to establish a robust connection between the spatial arrangement of molecular components in the nucleus and their functional roles. Achieving a comprehensive understanding of gene regulation requires capturing diverse data modalities within individual cells, enabling the direct inference of functional relationships between these components. In this context, imaging-based technologies have emerged as an especially promising approach for gathering spatial information across multiple components in the same cell.

3C methods in cancer research: recent advances and future prospects

In recent years, Hi-C technology has revolutionized cancer research by elucidating the mystery of three-dimensional chromatin organization and its role in gene regulation. This paper explored the impact of Hi-C advancements on cancer research by delving into high-resolution techniques, such as chromatin loops, structural variants, haplotype phasing, and extrachromosomal DNA (ecDNA). Distant regulatory elements interact with their target genes through chromatin loops. Structural variants contribute to the development and progression of cancer. Haplotype phasing is crucial for understanding allele-specific genomic rearrangements and somatic clonal evolution in cancer. The role of ecDNA in driving oncogene amplification and drug resistance in cancer cells has also been revealed. These innovations offer a deeper understanding of cancer biology and the potential for personalized therapies. Despite these advancements, challenges, such as the accurate mapping of repetitive sequences and precise identification of structural variants, persist. Integrating Hi-C with multiomics data is key to overcoming these challenges and comprehensively understanding complex cancer genomes. Thus, Hi-C is a powerful tool for guiding precision medicine in cancer research and treatment.

Multifaceted roles of cohesin in regulating transcriptional loops

Cohesin is required for chromatin loop formation. However, its precise role in regulating gene transcription remains largely unknown. We investigated the relationship between cohesin and RNA Polymerase II (RNAPII) using single-molecule mapping and live-cell imaging methods in human cells. Cohesin-mediated transcriptional loops were highly correlated with those of RNAPII and followed the direction of gene transcription. Depleting RAD21, a subunit of cohesin, resulted in the loss of long-range (>100 kb) loops between distal (super-)enhancers and promoters of cell-type-specific genes. By contrast, the short-range (<50 kb) loops were insensitive to RAD21 depletion and connected genes that are mostly housekeeping. This result explains why only a small fraction of genes are affected by the loss of long-range chromatin interactions due to cohesin depletion. Remarkably, RAD21 depletion appeared to up-regulate genes located in early initiation zones (EIZ) of DNA replication, and the EIZ signals were amplified drastically without RAD21. Our results revealed new mechanistic insights of cohesin's multifaceted roles in establishing transcriptional loops, preserving long-range chromatin interactions for cell-specific genes, and maintaining timely order of DNA replication.

MMCT-Loop: a mix model-based pipeline for calling targeted 3D chromatin loops

Protein-specific Chromatin Conformation Capture (3C)-based technologies have become essential for identifying distal genomic interactions with critical roles in gene regulation. The standard techniques include Chromatin Interaction Analysis by Paired-End Tag (ChIA-PET), in situ Hi-C followed by chromatin immunoprecipitation (HiChIP) also known as PLAC-seq. To identify chromatin interactions from these data, a variety of computational methods have emerged. Although these state-of-art methods address many issues with loop calling, only few methods can fit different data types simultaneously, and the accuracy as well as the efficiency these approaches remains limited. Here we have generated a pipeline, MMCT-Loop, which ensures the accurate identification of strong loops as well as dynamic or weak loops through a mixed model. MMCT-Loop outperforms existing methods in accuracy, and the detected loops show higher activation functionality. To highlight the utility of MMCT-Loop, we applied it to conformational data derived from neural stem cell (NSCs) and uncovered several previously unidentified regulatory regions for key master regulators of stem cell identity. MMCT-Loop is an accurate and efficient loop caller for targeted conformation capture data, which supports raw data or pre-processed valid pairs as input, the output interactions are formatted and easily uploaded to a genome browser for visualization.

Machine and deep learning methods for predicting 3D genome organization

Three-Dimensional (3D) chromatin interactions, such as enhancer-promoter interactions (EPIs), loops, Topologically Associating Domains (TADs), and A/B compartments play critical roles in a wide range of cellular processes by regulating gene expression. Recent development of chromatin conformation capture technologies has enabled genome-wide profiling of various 3D structures, even with single cells. However, current catalogs of 3D structures remain incomplete and unreliable due to differences in technology, tools, and low data resolution. Machine learning methods have emerged as an alternative to obtain missing 3D interactions and/or improve resolution. Such methods frequently use genome annotation data (ChIP-seq, DNAse-seq, etc.), DNA sequencing information (k-mers, Transcription Factor Binding Site (TFBS) motifs), and other genomic properties to learn the associations between genomic features and chromatin interactions. In this review, we discuss computational tools for predicting three types of 3D interactions (EPIs, chromatin interactions, TAD boundaries) and analyze their pros and cons. We also point out obstacles of computational prediction of 3D interactions and suggest future research directions.

An in silico procedure for generating protein-mediated chromatin interaction data and comparison of significant interaction calling methods

The ability to simulate high-throughput data with high fidelity to real experimental data is fundamental for benchmarking methods used to detect true long-range chromatin interactions mediated by a specific protein. Yet, such tools are not currently available. To fill this gap, we develop an in silico experimental procedure, ChIA-Sim, which imitates the experimental procedures that produce real ChIA-PET, Hi-ChIP, or PLAC-seq data. We show the fidelity of ChIA-Sim to real data by using guiding characteristics of several real datasets to generate data using the simulation procedure. We also used ChIA-Sim data to demonstrate the use of our in silico procedure in benchmarking methods for significant interactions analysis by evaluating four methods for significant interaction calling (SIC). In particular, we assessed each method's performance in terms of correct identification of long-range interactions. We further analyzed four experimental datasets from publicly available databases and shew that the trend of the results are consistent with those seen in data generated from ChIA-Sim. This serves as additional evidence that ChIA-Sim closely resembles data produced from the experimental protocols it models after.

Mapping nucleosome-resolution chromatin organization and enhancer-promoter loops in plants using Micro-C-XL

Although chromatin organizations in plants have been dissected at the scales of compartments and topologically associating domain (TAD)-like domains, there remains a gap in resolving fine-scale structures. Here, we use Micro-C-XL, a high-throughput chromosome conformation capture (Hi-C)-based technology that involves micrococcal nuclease (instead of restriction enzymes) and long cross-linkers, to dissect single nucleosome-resolution chromatin organization in Arabidopsis. Insulation analysis reveals more than 14,000 boundaries, which mostly include chromatin accessibility, epigenetic modifications, and transcription factors. Micro-C-XL reveals associations between RNA Pols and local chromatin organizations, suggesting that gene transcription substantially contributes to the establishment of local chromatin domains. By perturbing Pol II both genetically and chemically at the gene level, we confirm its function in regulating chromatin organization. Visible loops and stripes are assigned to super-enhancers and their targeted genes, thus providing direct insights for the identification and mechanistic analysis of distal CREs and their working modes in plants. We further investigate possible factors regulating these chromatin loops. Subsequently, we expand Micro-C-XL to soybean and rice. In summary, we use Micro-C-XL for analyses of plants, which reveal fine-scale chromatin organization and enhancer-promoter loops and provide insights regarding three-dimensional genomes in plants.

Nuclear RNA catabolism controls endogenous retroviruses, gene expression asymmetry, and dedifferentiation

Endogenous retroviruses (ERVs) are remnants of ancient parasitic infections and comprise sizable portions of most genomes. Although epigenetic mechanisms silence most ERVs by generating a repressive environment that prevents their expression (heterochromatin), little is known about mechanisms silencing ERVs residing in open regions of the genome (euchromatin). This is particularly important during embryonic development, where induction and repression of distinct classes of ERVs occur in short temporal windows. Here, we demonstrate that transcription-associated RNA degradation by the nuclear RNA exosome and Integrator is a regulatory mechanism that controls the productive transcription of most genes and many ERVs involved in preimplantation development. Disrupting nuclear RNA catabolism promotes dedifferentiation to a totipotent-like state characterized by defects in RNAPII elongation and decreased expression of long genes (gene-length asymmetry). Our results indicate that RNA catabolism is a core regulatory module of gene networks that safeguards RNAPII activity, ERV expression, cell identity, and developmental potency.

Spatial promoter-enhancer hubs in cancer: organization, regulation, and function

Transcriptional dysregulation is a hallmark of cancer and can be driven by altered enhancer landscapes. Recent studies in genome organization have revealed that multiple enhancers and promoters can spatially coalesce to form dynamic topological assemblies, known as promoter-enhancer hubs, which strongly correlate with elevated gene expression. In this review, we discuss the structure and complexity of promoter-enhancer hubs recently identified in multiple cancer types. We further discuss underlying mechanisms driving dysregulation of promoter-enhancer hubs and speculate on their functional role in pathogenesis. Understanding the role of promoter-enhancer hubs in transcriptional dysregulation can provide insight into new therapeutic approaches to target these complex features of genome organization.

MSL2 ensures biallelic gene expression in mammals

In diploid organisms, biallelic gene expression enables the production of adequate levels of mRNA1,2. This is essential for haploinsufficient genes, which require biallelic expression for optimal function to prevent the onset of developmental disorders1,3. Whether and how a biallelic or monoallelic state is determined in a cell-type-specific manner at individual loci remains unclear. MSL2 is known for dosage compensation of the male X chromosome in flies. Here we identify a role of MSL2 in regulating allelic expression in mammals. Allele-specific bulk and single-cell analyses in mouse neural progenitor cells revealed that, in addition to the targets showing biallelic downregulation, a class of genes transitions from biallelic to monoallelic expression after MSL2 loss. Many of these genes are haploinsufficient. In the absence of MSL2, one allele remains active, retaining active histone modifications and transcription factor binding, whereas the other allele is silenced, exhibiting loss of promoter-enhancer contacts and the acquisition of DNA methylation. Msl2-knockout mice show perinatal lethality and heterogeneous phenotypes during embryonic development, supporting a role for MSL2 in regulating gene dosage. The role of MSL2 in preserving biallelic expression of specific dosage-sensitive genes sets the stage for further investigation of other factors that are involved in allelic dosage compensation in mammalian cells, with considerable implications for human disease.

Ultrafine mapping of chromosome conformation at hundred basepair resolution reveals regulatory genome architecture

The resolution limit of chromatin conformation capture methodologies (3Cs) has restrained their application in detection of fine-level chromatin structure mediated by cis-regulatory elements (CREs). Here, we report two 3C-derived methods, Tri-4C and Tri-HiC, which utilize multirestriction enzyme digestions for ultrafine mapping of targeted and genome-wide chromatin interaction, respectively, at up to one hundred basepair resolution. Tri-4C identified CRE loop interaction networks and quantitatively revealed their alterations underlying dynamic gene control. Tri-HiC uncovered global fine-gauge regulatory interaction networks, identifying >20-fold more enhancer:promoter (E:P) loops than in situ Hi-C. In addition to vastly improved identification of subkilobase-sized E:P loops, Tri-HiC also uncovered interaction stripes and contact domain insulation from promoters and enhancers, revealing their loop extrusion behaviors resembling the topologically associating domain boundaries. Tri-4C and Tri-HiC provide robust approaches to achieve the high-resolution interactome maps required for characterizing fine-gauge regulatory chromatin interactions in analysis of development, homeostasis, and disease.

Foxp3 orchestrates reorganization of chromatin architecture to establish regulatory T cell identity

Chromatin conformation reorganization is emerging as an important layer of regulation for gene expression and lineage specification. Yet, how lineage-specific transcription factors contribute to the establishment of cell type-specific 3D chromatin architecture in the immune cells remains unclear, especially for the late stages of T cell subset differentiation and maturation. Regulatory T cells (Treg) are mainly generated in the thymus as a subpopulation of T cells specializing in suppressing excessive immune responses. Here, by comprehensively mapping 3D chromatin organization during Treg cell differentiation, we show that Treg-specific chromatin structures were progressively established during its lineage specification, and highly associated with Treg signature gene expression. Additionally, the binding sites of Foxp3, a Treg lineage specifying transcription factor, were highly enriched at Treg-specific chromatin loop anchors. Further comparison of the chromatin interactions between wide-type Tregs versus Treg cells from Foxp3 knock-in/knockout or newly-generated Foxp3 domain-swap mutant mouse revealed that Foxp3 was essential for the establishment of Treg-specific 3D chromatin architecture, although it was not dependent on the formation of the Foxp3 domain-swapped dimer. These results highlighted an underappreciated role of Foxp3 in modulating Treg-specific 3D chromatin structure formation.

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.

3D organization of regulatory elements for transcriptional regulation in Arabidopsis

BACKGROUND

Although spatial organization of compartments and topologically associating domains at large scale is relatively well studied, the spatial organization of regulatory elements at fine scale is poorly understood in plants.

RESULTS

Here we perform high-resolution chromatin interaction analysis using paired-end tag sequencing approach. We map chromatin interactions tethered with RNA polymerase II and associated with heterochromatic, transcriptionally active, and Polycomb-repressive histone modifications in Arabidopsis. Analysis of the regulatory repertoire shows that distal active cis-regulatory elements are linked to their target genes through long-range chromatin interactions with increased expression of the target genes, while poised cis-regulatory elements are linked to their target genes through long-range chromatin interactions with depressed expression of the target genes. Furthermore, we demonstrate that transcription factor MYC2 is critical for chromatin spatial organization, and propose that MYC2 occupancy and MYC2-mediated chromatin interactions coordinately facilitate transcription within the framework of 3D chromatin architecture. Analysis of functionally related gene-defined chromatin connectivity networks reveals that genes implicated in flowering-time control are functionally compartmentalized into separate subdomains via their spatial activity in the leaf or shoot apical meristem, linking active mark- or Polycomb-repressive mark-associated chromatin conformation to coordinated gene expression.

CONCLUSION

The results reveal that the regulation of gene transcription in Arabidopsis is not only by linear juxtaposition, but also by long-range chromatin interactions. Our study uncovers the fine scale genome organization of Arabidopsis and the potential roles of such organization in orchestrating transcription and development.

Spatial and temporal organization of the genome: Current state and future aims of the 4D nucleome project

The four-dimensional nucleome (4DN) consortium studies the architecture of the genome and the nucleus in space and time. We summarize progress by the consortium and highlight the development of technologies for (1) mapping genome folding and identifying roles of nuclear components and bodies, proteins, and RNA, (2) characterizing nuclear organization with time or single-cell resolution, and (3) imaging of nuclear organization. With these tools, the consortium has provided over 2,000 public datasets. Integrative computational models based on these data are starting to reveal connections between genome structure and function. We then present a forward-looking perspective and outline current aims to (1) delineate dynamics of nuclear architecture at different timescales, from minutes to weeks as cells differentiate, in populations and in single cells, (2) characterize cis-determinants and trans-modulators of genome organization, (3) test functional consequences of changes in cis- and trans-regulators, and (4) develop predictive models of genome structure and function.

Evaluating the role of the nuclear microenvironment in gene function by population-based modeling

The nuclear folding of chromosomes relative to nuclear bodies is an integral part of gene function. Here, we demonstrate that population-based modeling-from ensemble Hi-C data-provides a detailed description of the nuclear microenvironment of genes and its role in gene function. We define the microenvironment by the subnuclear positions of genomic regions with respect to nuclear bodies, local chromatin compaction, and preferences in chromatin compartmentalization. These structural descriptors are determined in single-cell models, thereby revealing the structural variability between cells. We demonstrate that the microenvironment of a genomic region is linked to its functional potential in gene transcription, replication, and chromatin compartmentalization. Some chromatin regions feature a strong preference for a single microenvironment, due to association with specific nuclear bodies in most cells. Other chromatin shows high structural variability, which is a strong indicator of functional heterogeneity. Moreover, we identify specialized nuclear microenvironments, which distinguish chromatin in different functional states and reveal a key role of nuclear speckles in chromosome organization. We demonstrate that our method produces highly predictive three-dimensional genome structures, which accurately reproduce data from a variety of orthogonal experiments, thus considerably expanding the range of Hi-C data analysis.

Quantitative control of Ets1 dosage by a multi-enhancer hub promotes Th1 cell differentiation and protects from allergic inflammation

Multi-enhancer hubs are spatial clusters of enhancers present across numerous developmental programs. Here, we studied the functional relevance of these three-dimensional structures in T cell biology. Mathematical modeling identified a highly connected multi-enhancer hub at the Ets1 locus, comprising a noncoding regulatory element that was a hotspot for sequence variation associated with allergic disease in humans. Deletion of this regulatory element in mice revealed that the multi-enhancer connectivity was dispensable for T cell development but required for CD4+ T helper 1 (Th1) differentiation. These mice were protected from Th1-mediated colitis but exhibited overt allergic responses. Mechanistically, the multi-enhancer hub controlled the dosage of Ets1 that was required for CTCF recruitment and assembly of Th1-specific genome topology. Our findings establish a paradigm wherein multi-enhancer hubs control cellular competence to respond to an inductive cue through quantitative control of gene dosage and provide insight into how sequence variation within noncoding elements at the Ets1 locus predisposes individuals to allergic responses.

SALL1 enforces microglia-specific DNA binding and function of SMADs to establish microglia identity

Spalt-like transcription factor 1 (SALL1) is a critical regulator of organogenesis and microglia identity. Here we demonstrate that disruption of a conserved microglia-specific super-enhancer interacting with the Sall1 promoter results in complete and specific loss of Sall1 expression in microglia. By determining the genomic binding sites of SALL1 and leveraging Sall1 enhancer knockout mice, we provide evidence for functional interactions between SALL1 and SMAD4 required for microglia-specific gene expression. SMAD4 binds directly to the Sall1 super-enhancer and is required for Sall1 expression, consistent with an evolutionarily conserved requirement of the TGFβ and SMAD homologs Dpp and Mad for cell-specific expression of Spalt in the Drosophila wing. Unexpectedly, SALL1 in turn promotes binding and function of SMAD4 at microglia-specific enhancers while simultaneously suppressing binding of SMAD4 to enhancers of genes that become inappropriately activated in enhancer knockout microglia, thereby enforcing microglia-specific functions of the TGFβ-SMAD signaling axis.

Recent advances in chromosome capture techniques unraveling 3D genome architecture in germ cells, health, and disease

In eukaryotes, the genome does not emerge in a specific shape but rather as a hierarchial bundle within the nucleus. This multifaceted genome organization consists of multiresolution cellular structures, such as chromosome territories, compartments, and topologically associating domains, which are frequently defined by architecture, design proteins including CTCF and cohesin, and chromatin loops. This review briefly discusses the advances in understanding the basic rules of control, chromatin folding, and functional areas in early embryogenesis. With the use of chromosome capture techniques, the latest advancements in technologies for visualizing chromatin interactions come close to revealing 3D genome formation frameworks with incredible detail throughout all genomic levels, including at single-cell resolution. The possibility of detecting variations in chromatin architecture might open up new opportunities for disease diagnosis and prevention, infertility treatments, therapeutic approaches, desired exploration, and many other application scenarios.

Combinatorial transcription factor binding encodes cis-regulatory wiring of forebrain GABAergic neurogenesis

Transcription factors (TFs) bind combinatorially to genomic cis-regulatory elements (cREs), orchestrating transcription programs. While studies of chromatin state and chromosomal interactions have revealed dynamic neurodevelopmental cRE landscapes, parallel understanding of the underlying TF binding lags. To elucidate the combinatorial TF-cRE interactions driving mouse basal ganglia development, we integrated ChIP-seq for twelve TFs, H3K4me3-associated enhancer-promoter interactions, chromatin and transcriptional state, and transgenic enhancer assays. We identified TF-cREs modules with distinct chromatin features and enhancer activity that have complementary roles driving GABAergic neurogenesis and suppressing other developmental fates. While the majority of distal cREs were bound by one or two TFs, a small proportion were extensively bound, and these enhancers also exhibited exceptional evolutionary conservation, motif density, and complex chromosomal interactions. Our results provide new insights into how modules of combinatorial TF-cRE interactions activate and repress developmental expression programs and demonstrate the value of TF binding data in modeling gene regulatory wiring.

Advancements in Focal Amplification Detection in Tumor/Liquid Biopsies and Emerging Clinical Applications

Focal amplifications (FAs) are crucial in cancer research due to their significant diagnostic, prognostic, and therapeutic implications. FAs manifest in various forms, such as episomes, double minute chromosomes, and homogeneously staining regions, arising through different mechanisms and mainly contributing to cancer cell heterogeneity, the leading cause of drug resistance in therapy. Numerous wet-lab, mainly FISH, PCR-based assays, next-generation sequencing, and bioinformatics approaches have been set up to detect FAs, unravel the internal structure of amplicons, assess their chromatin compaction status, and investigate the transcriptional landscape associated with their occurrence in cancer cells. Most of them are tailored for tumor samples, even at the single-cell level. Conversely, very limited approaches have been set up to detect FAs in liquid biopsies. This evidence suggests the need to improve these non-invasive investigations for early tumor detection, monitoring disease progression, and evaluating treatment response. Despite the potential therapeutic implications of FAs, such as, for example, the use of HER2-specific compounds for patients with ERBB2 amplification, challenges remain, including developing selective and effective FA-targeting agents and understanding the molecular mechanisms underlying FA maintenance and replication. This review details a state-of-the-art of FA investigation, with a particular focus on liquid biopsies and single-cell approaches in tumor samples, emphasizing their potential to revolutionize the future diagnosis, prognosis, and treatment of cancer patients.

Genetic insights into immune mechanisms of Alzheimer's and Parkinson's disease

Microglia, the macrophages of the brain, are vital for brain homeostasis and have been implicated in a broad range of brain disorders. Neuroinflammation has gained traction as a possible therapeutic target for neurodegeneration, however, the precise function of microglia in specific neurodegenerative disorders is an ongoing area of research. Genetic studies offer valuable insights into understanding causality, rather than merely observing a correlation. Genome-wide association studies (GWAS) have identified many genetic loci that are linked to susceptibility to neurodegenerative disorders. (Post)-GWAS studies have determined that microglia likely play an important role in the development of Alzheimer's disease (AD) and Parkinson's disease (PD). The process of understanding how individual GWAS risk loci affect microglia function and mediate susceptibility is complex. A rapidly growing number of publications with genomic datasets and computational tools have formulated new hypotheses that guide the biological interpretation of AD and PD genetic risk. In this review, we discuss the key concepts and challenges in the post-GWAS interpretation of AD and PD GWAS risk alleles. Post-GWAS challenges include the identification of target cell (sub)type(s), causal variants, and target genes. Crucially, the prediction of GWAS-identified disease-risk cell types, variants and genes require validation and functional testing to understand the biological consequences within the pathology of the disorders. Many AD and PD risk genes are highly pleiotropic and perform multiple important functions that might not be equally relevant for the mechanisms by which GWAS risk alleles exert their effect(s). Ultimately, many GWAS risk alleles exert their effect by changing microglia function, thereby altering the pathophysiology of these disorders, and hence, we believe that modelling this context is crucial for a deepened understanding of these disorders.

Determining chromatin architecture with Micro Capture-C

Micro Capture-C (MCC) is a chromatin conformation capture (3C) method for visualizing reproducible three-dimensional contacts of specified regions of the genome at base pair resolution. These methods are an established family of techniques that use proximity ligation to assay the topology of chromatin. MCC can generate data at substantially higher resolution than previous techniques through multiple refinements of the 3C method. Using a sequence agnostic nuclease, the maintenance of cellular integrity and full sequencing of the ligation junctions, MCC achieves subnucleosomal levels of resolution, which can be used to reveal transcription factor binding sites analogous to DNAse I footprinting. Gene dense regions, close-range enhancer-promoter contacts, individual enhancers within super-enhancers and multiple other types of loci or regulatory regions that were previously challenging to assay with conventional 3C techniques, are readily observed using MCC. MCC requires training in common molecular biology techniques and bioinformatics to perform the experiment and analyze the data. The protocol can be expected to be completed in a 3 week timeframe for experienced molecular biologists.

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.

SCN1A: bioinformatically informed revised boundaries for promoter and enhancer regions

Pathogenic variations in the sodium voltage-gated channel alpha subunit 1 (SCN1A) gene are responsible for multiple epilepsy phenotypes, including Dravet syndrome, febrile seizures (FS) and genetic epilepsy with FS plus. Phenotypic heterogeneity is a hallmark of SCN1A-related epilepsies, the causes of which are yet to be clarified. Genetic variation in the non-coding regulatory regions of SCN1A could be one potential causal factor. However, a comprehensive understanding of the SCN1A regulatory landscape is currently lacking. Here, we summarized the current state of knowledge of SCN1A regulation, providing details on its promoter and enhancer regions. We then integrated currently available data on SCN1A promoters by extracting information related to the SCN1A locus from genome-wide repositories and clearly defined the promoter and enhancer regions of SCN1A. Further, we explored the cellular specificity of differential SCN1A promoter usage. We also reviewed and integrated the available human brain-derived enhancer databases and mouse-derived data to provide a comprehensive computationally developed summary of SCN1A brain-active enhancers. By querying genome-wide data repositories, extracting SCN1A-specific data and integrating the different types of independent evidence, we created a comprehensive catalogue that better defines the regulatory landscape of SCN1A, which could be used to explore the role of SCN1A regulatory regions in disease.

Characterization of altered molecular mechanisms in Parkinson's disease through cell type-resolved multiomics analyses

Parkinson's disease (PD) is a progressive neurodegenerative disorder. However, cell type-dependent transcriptional regulatory programs responsible for PD pathogenesis remain elusive. Here, we establish transcriptomic and epigenomic landscapes of the substantia nigra by profiling 113,207 nuclei obtained from healthy controls and patients with PD. Our multiomics data integration provides cell type annotation of 128,724 cis-regulatory elements (cREs) and uncovers cell type-specific dysregulations in cREs with a strong transcriptional influence on genes implicated in PD. The establishment of high-resolution three-dimensional chromatin contact maps identifies 656 target genes of dysregulated cREs and genetic risk loci, uncovering both potential and known PD risk genes. Notably, these candidate genes exhibit modular gene expression patterns with unique molecular signatures in distinct cell types, highlighting altered molecular mechanisms in dopaminergic neurons and glial cells including oligodendrocytes and microglia. Together, our single-cell transcriptome and epigenome reveal cell type-specific disruption in transcriptional regulations related to PD.

Highly connected 3D chromatin networks established by an oncogenic fusion protein shape tumor cell identity

Cell fate transitions observed in embryonic development involve changes in three-dimensional genomic organization that provide proper lineage specification. Whether similar events occur within tumor cells and contribute to cancer evolution remains largely unexplored. We modeled this process in the pediatric cancer Ewing sarcoma and investigated high-resolution looping and large-scale nuclear conformation changes associated with the oncogenic fusion protein EWS-FLI1. We show that chromatin interactions in tumor cells are dominated by highly connected looping hubs centered on EWS-FLI1 binding sites, which directly control the activity of linked enhancers and promoters to establish oncogenic expression programs. Conversely, EWS-FLI1 depletion led to the disassembly of these looping networks and a widespread nuclear reorganization through the establishment of new looping patterns and large-scale compartment configuration matching those observed in mesenchymal stem cells, a candidate Ewing sarcoma progenitor. Our data demonstrate that major architectural features of nuclear organization in cancer cells can depend on single oncogenes and are readily reversed to reestablish latent differentiation programs.

Foxp3 Orchestrates Reorganization of Chromatin Architecture to Establish Regulatory T Cell Identity

Chromatin conformation reorganization is emerging as an important layer of regulation for gene expression and lineage specification. Yet, how lineage-specific transcription factors contribute to the establishment of cell type-specific 3D chromatin architecture in the immune cells remains unclear, especially for the late stages of T cell subset differentiation and maturation. Regulatory T cells (Treg) are mainly generated in the thymus as a subpopulation of T cells specializing in suppressing excessive immune responses. Here, by comprehensively mapping 3D chromatin organization during Treg cell differentiation, we show that Treg-specific chromatin structures were progressively established during its lineage specification, and highly associated with Treg signature gene expression. Additionally, the binding sites of Foxp3, a Treg lineage specifying transcription factor, were highly enriched at Treg-specific chromatin loop anchors. Further comparison of the chromatin interactions between wide-type Tregs versus Treg cells from Foxp3 knock-in/knockout or newly-generated Foxp3 domain-swap mutant mouse revealed that Foxp3 was essential for the establishment of Treg-specific 3D chromatin architecture, although it was not dependent on the formation of the Foxp3 domain-swapped dimer. These results highlighted an underappreciated role of Foxp3 in modulating Treg-specific 3D chromatin structure formation.

High-throughput Pore-C reveals the single-allele topology and cell type-specificity of 3D genome folding

Canonical three-dimensional (3D) genome structures represent the ensemble average of pairwise chromatin interactions but not the single-allele topologies in populations of cells. Recently developed Pore-C can capture multiway chromatin contacts that reflect regional topologies of single chromosomes. By carrying out high-throughput Pore-C, we reveal extensive but regionally restricted clusters of single-allele topologies that aggregate into canonical 3D genome structures in two human cell types. We show that fragments in multi-contact reads generally coexist in the same TAD. In contrast, a concurrent significant proportion of multi-contact reads span multiple compartments of the same chromatin type over megabase distances. Synergistic chromatin looping between multiple sites in multi-contact reads is rare compared to pairwise interactions. Interestingly, the single-allele topology clusters are cell type-specific even inside highly conserved TADs in different types of cells. In summary, HiPore-C enables global characterization of single-allele topologies at an unprecedented depth to reveal elusive genome folding principles.

3D genomics and its applications in precision medicine

Three-dimensional (3D) genomics is an emerging discipline that studies the three-dimensional structure of chromatin and the three-dimensional and functions of genomes. It mainly studies the three-dimensional conformation and functional regulation of intranuclear genomes, such as DNA replication, DNA recombination, genome folding, gene expression regulation, transcription factor regulation mechanism, and the maintenance of three-dimensional conformation of genomes. Self-chromosomal conformation capture (3C) technology has been developed, and 3D genomics and related fields have developed rapidly. In addition, chromatin interaction analysis techniques developed by 3C technologies, such as paired-end tag sequencing (ChIA-PET) and whole-genome chromosome conformation capture (Hi-C), enable scientists to further study the relationship between chromatin conformation and gene regulation in different species. Thus, the spatial conformation of plant, animal, and microbial genomes, transcriptional regulation mechanisms, interaction patterns of chromosomes, and the formation mechanism of spatiotemporal specificity of genomes are revealed. With the help of new experimental technologies, the identification of key genes and signal pathways related to life activities and diseases is sustaining the rapid development of life science, agriculture, and medicine. In this paper, the concept and development of 3D genomics and its application in agricultural science, life science, and medicine are introduced, which provides a theoretical basis for the study of biological life processes.

Elucidating the structure and function of the nucleus-The NIH Common Fund 4D Nucleome program

Genomic architecture appears to play crucial roles in health and a variety of diseases. How nuclear structures reorganize over different timescales is elusive, partly because the tools needed to probe and perturb them are not as advanced as needed by the field. To fill this gap, the National Institutes of Health Common Fund started a program in 2015, called the 4D Nucleome (4DN), with the goal of developing and ultimately applying technologies to interrogate the structure and function of nuclear organization in space and time.

ChIATAC is an efficient strategy for multi-omics mapping of 3D epigenomes from low-cell inputs

Connecting genes to their cis-regulatory elements has been enabled by genome-wide mapping of chromatin interactions using proximity ligation in ChIA-PET, Hi-C, and their derivatives. However, these methods require millions of input cells for high-quality data and thus are unsuitable for many studies when only limited cells are available. Conversely, epigenomic profiling via transposase digestion in ATAC-seq requires only hundreds to thousands of cells to robustly map open chromatin associated with transcription activity, but it cannot directly connect active genes to their distal enhancers. Here, we combine proximity ligation in ChIA-PET and transposase accessibility in ATAC-seq into ChIATAC to efficiently map interactions between open chromatin loci in low numbers of input cells. We validate ChIATAC in Drosophila cells and optimize it for mapping 3D epigenomes in human cells robustly. Applying ChIATAC to primary human T cells, we reveal mechanisms that topologically regulate transcriptional programs during T cell activation.

HPTAD: A computational method to identify topologically associating domains from HiChIP and PLAC-seq datasets

High-throughput chromatin conformation capture technologies, such as Hi-C and Micro-C, have enabled genome-wide view of chromatin spatial organization. Most recently, Hi-C-derived enrichment-based technologies, including HiChIP and PLAC-seq, offer attractive alternatives due to their high signal-to-noise ratio and low cost. While a series of computational tools have been developed for Hi-C data, methods tailored for HiChIP and PLAC-seq data are still under development. Here we present HPTAD, a computational method to identify topologically associating domains (TADs) from HiChIP and PLAC-seq data. We performed comprehensive benchmark analysis to demonstrate its superior performance over existing TAD callers designed for Hi-C data. HPTAD is freely available at https://github.com/yunliUNC/HPTAD.

Hooked Up from a Distance: Charting Genome-Wide Long-Range Interaction Maps in Neural Cells Chromatin to Identify Novel Candidate Genes for Neurodevelopmental Disorders

DNA sequence variants (single nucleotide polymorphisms or variants, SNPs/SNVs; copy number variants, CNVs) associated to neurodevelopmental disorders (NDD) and traits often map on putative transcriptional regulatory elements, including, in particular, enhancers. However, the genes controlled by these enhancers remain poorly defined. Traditionally, the activity of a given enhancer, and the effect of its possible alteration associated to the sequence variants, has been thought to influence the nearest gene promoter. However, the obtainment of genome-wide long-range interaction maps in neural cells chromatin challenged this view, showing that a given enhancer is very frequently not connected to the nearest promoter, but to a more distant one, skipping genes in between. In this Perspective, we review some recent papers, who generated long-range interaction maps (by HiC, RNApolII ChIA-PET, Capture-HiC, or PLACseq), and overlapped the identified long-range interacting DNA segments with DNA sequence variants associated to NDD (such as schizophrenia, bipolar disorder and autism) and traits (intelligence). This strategy allowed to attribute the function of enhancers, hosting the NDD-related sequence variants, to a connected gene promoter lying far away on the linear chromosome map. Some of these enhancer-connected genes had indeed been already identified as contributive to the diseases, by the identification of mutations within the gene's protein-coding regions (exons), validating the approach. Significantly, however, the connected genes also include many genes that were not previously found mutated in their exons, pointing to novel candidate contributors to NDD and traits. Thus, long-range interaction maps, in combination with DNA variants detected in association with NDD, can be used as "pointers" to identify novel candidate disease-relevant genes. Functional manipulation of the long-range interaction network involving enhancers and promoters by CRISPR-Cas9-based approaches is beginning to probe for the functional significance of the identified interactions, and the enhancers and the genes involved, improving our understanding of neural development and its pathology.

Аpplication of massive parallel reporter analysis in biotechnology and medicine

The development and functioning of an organism relies on tissue-specific gene programs. Genome regulatory elements play a key role in the regulation of such programs, and disruptions in their function can lead to the development of various pathologies, including cancers, malformations and autoimmune diseases. The emergence of high-throughput genomic studies has led to massively parallel reporter analysis (MPRA) methods, which allow the functional verification and identification of regulatory elements on a genome-wide scale. Initially MPRA was used as a tool to investigate fundamental aspects of epigenetics, but the approach also has great potential for clinical and practical biotechnology. Currently, MPRA is used for validation of clinically significant mutations, identification of tissue-specific regulatory elements, search for the most promising loci for transgene integration, and is an indispensable tool for creating highly efficient expression systems, the range of application of which extends from approaches for protein development and design of next-generation therapeutic antibody superproducers to gene therapy. In this review, the main principles and areas of practical application of high-throughput reporter assays will be discussed.

Cooperation of chromatin remodeling SWI/SNF complex and pioneer factor AP-1 shapes 3D enhancer landscapes

The mechanism controlling the dynamic targeting of SWI/SNF has long been postulated to be coordinated by transcription factors (TFs), yet demonstrating a specific TF influence has proven difficult. Here we take a multi-omics approach to interrogate transient SWI/SNF interactors, chromatin targeting and the resulting three-dimensional epigenetic landscape. We utilize the labeling technique TurboID to map the SWI/SNF interactome and identify the activator protein-1 (AP-1) family members as critical interacting partners for SWI/SNF complexes. CUT&RUN profiling demonstrates SWI/SNF targeting enrichment at AP-1 bound loci, as well as SWI/SNF-AP-1 cooperation in chromatin targeting. HiChIP reveals AP-1-SWI/SNF-dependent restructuring of the three-dimensional promoter-enhancer architecture and generation of enhancer hubs. Through interrogation of the SWI/SNF-AP-1 interaction, we demonstrate an SWI/SNF dependency on AP-1-mediated chromatin localization. We propose that pioneer factors, such as AP-1, bind and target SWI/SNF to inactive chromatin, where it restructures the genomic landscape into an active state through epigenetic rewiring spanning multiple dimensions.

The spatial organization of transcriptional control

In animals, the sequences for controlling gene expression do not concentrate just at the transcription start site of genes, but are frequently thousands to millions of base pairs distal to it. The interaction of these sequences with one another and their transcription start sites is regulated by factors that shape the three-dimensional (3D) organization of the genome within the nucleus. Over the past decade, indirect tools exploiting high-throughput DNA sequencing have helped to map this 3D organization, have identified multiple key regulators of its structure and, in the process, have substantially reshaped our view of how 3D genome architecture regulates transcription. Now, new tools for high-throughput super-resolution imaging of chromatin have directly visualized the 3D chromatin organization, settling some debates left unresolved by earlier indirect methods, challenging some earlier models of regulatory specificity and creating hypotheses about the role of chromatin structure in transcriptional regulation.

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.