High-resolution targeted 3C interrogation of cis-regulatory element organization at genome-wide scale

3D Enhancer-promoter networks provide predictive features for gene expression and coregulation in early embryonic lineages

Mammalian embryogenesis commences with two pivotal and binary cell fate decisions that give rise to three essential lineages: the trophectoderm, the epiblast and the primitive endoderm. Although key signaling pathways and transcription factors that control these early embryonic decisions have been identified, the non-coding regulatory elements through which transcriptional regulators enact these fates remain understudied. Here, we characterize, at a genome-wide scale, enhancer activity and 3D connectivity in embryo-derived stem cell lines that represent each of the early developmental fates. We observe extensive enhancer remodeling and fine-scale 3D chromatin rewiring among the three lineages, which strongly associate with transcriptional changes, although distinct groups of genes are irresponsive to topological changes. In each lineage, a high degree of connectivity, or 'hubness', positively correlates with levels of gene expression and enriches for cell-type specific and essential genes. Genes within 3D hubs also show a significantly stronger probability of coregulation across lineages compared to genes in linear proximity or within the same contact domains. By incorporating 3D chromatin features, we build a predictive model for transcriptional regulation (3D-HiChAT) that outperforms models using only 1D promoter or proximal variables to predict levels and cell-type specificity of gene expression. Using 3D-HiChAT, we identify, in silico, candidate functional enhancers and hubs in each cell lineage, and with CRISPRi experiments, we validate several enhancers that control gene expression in their respective lineages. Our study identifies 3D regulatory hubs associated with the earliest mammalian lineages and describes their relationship to gene expression and cell identity, providing a framework to comprehensively understand lineage-specific transcriptional behaviors.

Technological advances in probing 4D genome organization

The last two decades of work on chromosome conformation in eukaryotic nuclei have revealed a complex and highly regulated hierarchy of architectural features, from self-associating domains and compartmental interactions to locus-specific loops. Recent findings have shown that these structures are dynamic and heterogeneous, with emerging insights into the factors that shape them and implications for the control of transcription and other nuclear processes. Here, we review the latest advances in the DNA sequencing- and microscopy-based technologies for probing these features in space and time (4D) and discuss how they have been combined with complementary approaches such as genetic perturbations, protein and RNA measurements, and modeling to gain mechanistic insights about genome regulation across space and time.

A sensitive mass spectrometry-based method to identify common respiratory pathogens in children

Public health threats posed by emerging respiratory infections are a significant concern, particularly in children and infants. Traditional culture-based detection methods are time-consuming and typically require 1-3 days. Herein, we developed and evaluated a 23-plex common respiratory pathogen mass spectrometry assay that enables the simultaneous detection of 18 common respiratory pathogens in children. This assay combines matrix-assisted laser desorption/ionization time of flight mass spectrometry with multiplex reverse transcription-PCR and targets 11 bacterial and 7 viral pathogens (including 10 subtypes), and two internal controls. The detection limit of the common respiratory pathogen mass spectrometry assay was as low as 1 copy/µL, with no cross-reactivity with other organisms. We assessed the clinical performance of the common respiratory pathogen mass spectrometry assay using respiratory samples from 450 children. The total 450 clinical specimens underwent analysis via matrix-assisted laser desorption/ionization time of flight mass spectrometry, and the outcomes were juxtaposed with those derived from real-time reverse-transcriptase PCR conducted concurrently. The concordance between these methods was 96.0%, and the multiple infection identification rate was 7.1%. This innovative approach enables the simultaneous analysis of numerous outcomes from a solitary examination across 192 specimens within a timeframe of approximately 7 hours, with a dramatically reduced sample use and cost. In summary, the common respiratory pathogen mass spectrometry assay is a sensitive, accurate, and cost-effective method for detecting common respiratory pathogens in children and has the potential to revolutionize the diagnosis of respiratory tract infections. IMPORTANCE This study aimed to present and evaluate a novel co-detection method that enables the simultaneous identification of 11 bacterial and 7 viral pathogens in about 7 hours using matrix-assisted laser desorption/ionization time of flight mass spectrometry. Our approach utilizes a combination of multiplex reverse transcription-PCR and matrix-assisted laser desorption/ionization time of flight mass spectrometry, which overcomes the limitations of conventional assays, which include a long assessment time, technical difficulty, and high costs. As a screening method for common respiratory pathogens in children, common respiratory pathogen mass spectrometry assay has the potential to revolutionize the diagnosis of respiratory tract infections by providing an accurate etiological diagnosis. The common respiratory pathogen mass spectrometry assay is expected to be a critical tool for the diagnosis of respiratory infections in children, offering a more efficient, cost-effective, and accurate approach for the detection of common respiratory pathogens.

Inferring causal genes at type 2 diabetes GWAS loci through chromosome interactions in islet cells

Background: Resolving causal genes for type 2 diabetes at loci implicated by genome-wide association studies (GWAS) requires integrating functional genomic data from relevant cell types. Chromatin features in endocrine cells of the pancreatic islet are particularly informative and recent studies leveraging chromosome conformation capture (3C) with Hi-C based methods have elucidated regulatory mechanisms in human islets. However, these genome-wide approaches are less sensitive and afford lower resolution than methods that target specific loci. Methods: To gauge the extent to which targeted 3C further resolves chromatin-mediated regulatory mechanisms at GWAS loci, we generated interaction profiles at 23 loci using next-generation (NG) capture-C in a human beta cell model (EndoC-βH1) and contrasted these maps with Hi-C maps in EndoC-βH1 cells and human islets and a promoter capture Hi-C map in human islets. Results: We found improvements in assay sensitivity of up to 33-fold and resolved ~3.6X more chromatin interactions. At a subset of 18 loci with 25 co-localised GWAS and eQTL signals, NG Capture-C interactions implicated effector transcripts at five additional genetic signals relative to promoter capture Hi-C through physical contact with gene promoters. Conclusions: High resolution chromatin interaction profiles at selectively targeted loci can complement genome- and promoter-wide maps.

Systematic mapping and modeling of 3D enhancer-promoter interactions in early mouse embryonic lineages reveal regulatory principles that determine the levels and cell-type specificity of gene expression

Mammalian embryogenesis commences with two pivotal and binary cell fate decisions that give rise to three essential lineages, the trophectoderm (TE), the epiblast (EPI) and the primitive endoderm (PrE). Although key signaling pathways and transcription factors that control these early embryonic decisions have been identified, the non-coding regulatory elements via which transcriptional regulators enact these fates remain understudied. To address this gap, we have characterized, at a genome-wide scale, enhancer activity and 3D connectivity in embryo-derived stem cell lines that represent each of the early developmental fates. We observed extensive enhancer remodeling and fine-scale 3D chromatin rewiring among the three lineages, which strongly associate with transcriptional changes, although there are distinct groups of genes that are irresponsive to topological changes. In each lineage, a high degree of connectivity or "hubness" positively correlates with levels of gene expression and enriches for cell-type specific and essential genes. Genes within 3D hubs also show a significantly stronger probability of coregulation across lineages, compared to genes in linear proximity or within the same contact domains. By incorporating 3D chromatin features, we build a novel predictive model for transcriptional regulation (3D-HiChAT), which outperformed models that use only 1D promoter or proximal variables in predicting levels and cell-type specificity of gene expression. Using 3D-HiChAT, we performed genome-wide in silico perturbations to nominate candidate functional enhancers and hubs in each cell lineage, and with CRISPRi experiments we validated several novel enhancers that control expression of one or more genes in their respective lineages. Our study comprehensively identifies 3D regulatory hubs associated with the earliest mammalian lineages and describes their relationship to gene expression and cell identity, providing a framework to understand lineage-specific transcriptional behaviors.

The Mediator complex regulates enhancer-promoter interactions

Enhancer-mediated gene activation generally requires physical proximity between enhancers and their target gene promoters. However, the molecular mechanisms by which interactions between enhancers and promoters are formed are not well understood. Here, we investigate the function of the Mediator complex in the regulation of enhancer-promoter interactions, by combining rapid protein depletion and high-resolution MNase-based chromosome conformation capture approaches. We show that depletion of Mediator leads to reduced enhancer-promoter interaction frequencies, which are associated with a strong decrease in gene expression. In addition, we find increased interactions between CTCF-binding sites upon Mediator depletion. These changes in chromatin architecture are associated with a redistribution of the Cohesin complex on chromatin and a reduction in Cohesin occupancy at enhancers. Together, our results indicate that the Mediator and Cohesin complexes contribute to enhancer-promoter interactions and provide insights into the molecular mechanisms by which communication between enhancers and promoters is regulated.

Multiplex-GAM: genome-wide identification of chromatin contacts yields insights overlooked by Hi-C

Technology for measuring 3D genome topology is increasingly important for studying gene regulation, for genome assembly and for mapping of genome rearrangements. Hi-C and other ligation-based methods have become routine but have specific biases. Here, we develop multiplex-GAM, a faster and more affordable version of genome architecture mapping (GAM), a ligation-free technique that maps chromatin contacts genome-wide. We perform a detailed comparison of multiplex-GAM and Hi-C using mouse embryonic stem cells. When examining the strongest contacts detected by either method, we find that only one-third of these are shared. The strongest contacts specifically found in GAM often involve 'active' regions, including many transcribed genes and super-enhancers, whereas in Hi-C they more often contain 'inactive' regions. Our work shows that active genomic regions are involved in extensive complex contacts that are currently underestimated in ligation-based approaches, and highlights the need for orthogonal advances in genome-wide contact mapping technologies.

Information about immune cell proportions and tumor stage improves the prediction of recurrence in patients with colorectal cancer

Predicting cancer recurrence is essential to improving the clinical outcomes of patients with colorectal cancer (CRC). Although tumor stage information has been used as a guideline to predict CRC recurrence, patients with the same stage show different clinical outcomes. Therefore, there is a need to develop a method to identify additional features for CRC recurrence prediction. Here, we developed a network-integrated multiomics (NIMO) approach to select appropriate transcriptome signatures for better CRC recurrence prediction by comparing the methylation signatures of immune cells. We validated the performance of the CRC recurrence prediction based on two independent retrospective cohorts of 114 and 110 patients. Moreover, to confirm that the prediction was improved, we used both NIMO-based immune cell proportions and TNM (tumor, node, metastasis) stage data. This work demonstrates the importance of (1) using both immune cell composition and TNM stage data and (2) identifying robust immune cell marker genes to improve CRC recurrence prediction.

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.

New genetic and epigenetic insights into the chemokine system: the latest discoveries aiding progression toward precision medicine

Over the past thirty years, the importance of chemokines and their seven-transmembrane G protein-coupled receptors (GPCRs) has been increasingly recognized. Chemokine interactions with receptors trigger signaling pathway activity to form a network fundamental to diverse immune processes, including host homeostasis and responses to disease. Genetic and nongenetic regulation of both the expression and structure of chemokines and receptors conveys chemokine functional heterogeneity. Imbalances and defects in the system contribute to the pathogenesis of a variety of diseases, including cancer, immune and inflammatory diseases, and metabolic and neurological disorders, which render the system a focus of studies aiming to discover therapies and important biomarkers. The integrated view of chemokine biology underpinning divergence and plasticity has provided insights into immune dysfunction in disease states, including, among others, coronavirus disease 2019 (COVID-19). In this review, by reporting the latest advances in chemokine biology and results from analyses of a plethora of sequencing-based datasets, we outline recent advances in the understanding of the genetic variations and nongenetic heterogeneity of chemokines and receptors and provide an updated view of their contribution to the pathophysiological network, focusing on chemokine-mediated inflammation and cancer. Clarification of the molecular basis of dynamic chemokine-receptor interactions will help advance the understanding of chemokine biology to achieve precision medicine application in the clinic.

Gene silencing dynamics are modulated by transiently active regulatory elements

Transcriptional enhancers have been extensively characterized, but cis-regulatory elements involved in acute gene repression have received less attention. Transcription factor GATA1 promotes erythroid differentiation by activating and repressing distinct gene sets. Here, we study the mechanism by which GATA1 silences the proliferative gene Kit during murine erythroid cell maturation and define stages from initial loss of activation to heterochromatinization. We find that GATA1 inactivates a potent upstream enhancer but concomitantly creates a discrete intronic regulatory region marked by H3K27ac, short noncoding RNAs, and de novo chromatin looping. This enhancer-like element forms transiently and serves to delay Kit silencing. The element is ultimately erased via the FOG1/NuRD deacetylase complex, as revealed by the study of a disease-associated GATA1 variant. Hence, regulatory sites can be self-limiting by dynamic co-factor usage. Genome-wide analyses across cell types and species uncover transiently active elements at numerous genes during repression, suggesting that modulation of silencing kinetics is widespread.

Tissue-specific Grb10/Ddc insulator drives allelic architecture for cardiac development

Allele-specific expression of imprinted gene clusters is governed by gametic DNA methylation at master regulators called imprinting control regions (ICRs). Non-gametic or secondary differentially methylated regions (DMRs) at promoters and exonic regions reinforce monoallelic expression but do not control an entire cluster. Here, we unveil an unconventional secondary DMR that is indispensable for tissue-specific imprinting of two previously unlinked genes, Grb10 and Ddc. Using polymorphic mice, we mapped an intronic secondary DMR at Grb10 with paternal-specific CTCF binding (CBR2.3) that forms contacts with Ddc. Deletion of paternal CBR2.3 removed a critical insulator, resulting in substantial shifting of chromatin looping and ectopic enhancer-promoter contacts. Destabilized gene architecture precipitated abnormal Grb10-Ddc expression with developmental consequences in the heart and muscle. Thus, we redefine the Grb10-Ddc imprinting domain by uncovering an unconventional intronic secondary DMR that functions as an insulator to instruct the tissue-specific, monoallelic expression of multiple genes-a feature previously ICR exclusive.

Distinct roles for CKM-Mediator in controlling Polycomb-dependent chromosomal interactions and priming genes for induction

Precise control of gene expression underpins normal development. This relies on mechanisms that enable communication between gene promoters and other regulatory elements. In embryonic stem cells (ESCs), the cyclin-dependent kinase module Mediator complex (CKM-Mediator) has been reported to physically link gene regulatory elements to enable gene expression and also prime genes for induction during differentiation. Here, we show that CKM-Mediator contributes little to three-dimensional genome organization in ESCs, but it has a specific and essential role in controlling interactions between inactive gene regulatory elements bound by Polycomb repressive complexes (PRCs). These interactions are established by the canonical PRC1 (cPRC1) complex but rely on CKM-Mediator, which facilitates binding of cPRC1 to its target sites. Importantly, through separation-of-function experiments, we reveal that this collaboration between CKM-Mediator and cPRC1 in creating long-range interactions does not function to prime genes for induction during differentiation. Instead, we discover that priming relies on an interaction-independent mechanism whereby the CKM supports core Mediator engagement with gene promoters during differentiation to enable gene activation.

Complex chromosomal rearrangements induced by transposons in maize

Eukaryotic genomes are large and complex, and gene expression can be affected by multiple regulatory elements and their positions within the dynamic chromatin architecture. Transposable elements are known to play important roles in genome evolution, yet questions remain as to how transposable elements alter genome structure and affect gene expression. Previous studies have shown that genome rearrangements can be induced by Reversed Ends Transposition involving termini of Activator and related transposable elements in maize and other plants. Here, we show that complex alleles can be formed by the rapid and progressive accumulation of Activator-induced duplications and rearrangements. The p1 gene enhancer in maize can induce ectopic expression of the nearby p2 gene in pericarp tissue when placed near it via different structural rearrangements. By screening for p2 expression, we identified and studied 5 cases in which multiple sequential transposition events occurred and increased the p1 enhancer copy number. We see active p2 expression due to multiple copies of the p1 enhancer present near p2 in all 5 cases. The p1 enhancer effects are confirmed by the observation that loss of p2 expression is correlated with transposition-induced excision of the p1 enhancers. We also performed a targeted Chromosome Conformation Capture experiment to test the physical interaction between the p1 enhancer and p2 promoter region. Together, our results show that transposon-induced rearrangements can accumulate rapidly and progressively increase genetic variation important for genomic evolution.

Natural and Experimental Rewiring of Gene Regulatory Regions

The successful development and ongoing functioning of complex organisms depend on the faithful execution of the genetic code. A critical step in this process is the correct spatial and temporal expression of genes. The highly orchestrated transcription of genes is controlled primarily by cis-regulatory elements: promoters, enhancers, and insulators. The medical importance of this key biological process can be seen by the frequency with which mutations and inherited variants that alter cis-regulatory elements lead to monogenic and complex diseases and cancer. Here, we provide an overview of the methods available to characterize and perturb gene regulatory circuits. We then highlight mechanisms through which regulatory rewiring contributes to disease, and conclude with a perspective on how our understanding of gene regulation can be used to improve human health.

Extensive co-binding and rapid redistribution of NANOG and GATA6 during emergence of divergent lineages

Fate-determining transcription factors (TFs) can promote lineage-restricted transcriptional programs from common progenitor states. The inner cell mass (ICM) of mouse blastocysts co-expresses the TFs NANOG and GATA6, which drive the bifurcation of the ICM into either the epiblast (Epi) or the primitive endoderm (PrE), respectively. Here, we induce GATA6 in embryonic stem cells-that also express NANOG-to characterize how a state of co-expression of opposing TFs resolves into divergent lineages. Surprisingly, we find that GATA6 and NANOG co-bind at the vast majority of Epi and PrE enhancers, a phenomenon we also observe in blastocysts. The co-bound state is followed by eviction and repression of Epi TFs, and quick remodeling of chromatin and enhancer-promoter contacts thus establishing the PrE lineage while repressing the Epi fate. We propose that co-binding of GATA6 and NANOG at shared enhancers maintains ICM plasticity and promotes the rapid establishment of Epi- and PrE-specific transcriptional programs.

Making connections: enhancers in cellular differentiation

Deciphering the process by which hundreds of distinct cell types emerge from a single zygote to form a complex multicellular organism remains one of the greatest challenges in biological research. Enhancers are known to be central to cell type-specific gene expression, yet many questions regarding how these genomic elements interact both temporally and spatially with other cis- and trans-acting factors to control transcriptional activity during differentiation and development remain unanswered. Here, we review our current understanding of the role of enhancers and their interactions in this context and highlight recent progress achieved with experimental methods of unprecedented resolution.

Capture-C: a modular and flexible approach for high-resolution chromosome conformation capture

Chromosome conformation capture (3C) methods measure the spatial proximity between DNA elements in the cell nucleus. Many methods have been developed to sample 3C material, including the Capture-C family of protocols. Capture-C methods use oligonucleotides to enrich for interactions of interest from sequencing-ready 3C libraries. This approach is modular and has been adapted and optimized to work for sampling of disperse DNA elements (NuTi Capture-C), including from low cell inputs (LI Capture-C), as well as to generate Hi-C like maps for specific regions of interest (Tiled-C) and to interrogate multiway interactions (Tri-C). We present the design, experimental protocol and analysis pipeline for NuTi Capture-C in addition to the variations for generation of LI Capture-C, Tiled-C and Tri-C data. The entire procedure can be performed in 3 weeks and requires standard molecular biology skills and equipment, access to a next-generation sequencing platform, and basic bioinformatic skills. Implemented with other sequencing technologies, these methods can be used to identify regulatory interactions and to compare the structural organization of the genome in different cell types and genetic models.

Assessment of Multiway Interactions with Tri-C

Tri-C is a chromosome conformation capture (3C) approach that can efficiently identify multiway chromatin interactions with viewpoints of interest. As opposed to pair-wise interactions identified in methods such as Hi-C, 4C, and Capture-C, the detection of multiway interactions allows researchers to investigate how multiple cis-regulatory elements interact together in higher-order structures in single nuclei and address questions regarding structural cooperation between these elements. Here, we describe the procedure for designing and performing a Tri-C experiment.

Identification of LZTFL1 as a candidate effector gene at a COVID-19 risk locus

The severe acute respiratory syndrome coronavirus 2 (SARS‑CoV‑2) disease (COVID-19) pandemic has caused millions of deaths worldwide. Genome-wide association studies identified the 3p21.31 region as conferring a twofold increased risk of respiratory failure. Here, using a combined multiomics and machine learning approach, we identify the gain-of-function risk A allele of an SNP, rs17713054G>A, as a probable causative variant. We show with chromosome conformation capture and gene-expression analysis that the rs17713054-affected enhancer upregulates the interacting gene, leucine zipper transcription factor like 1 (LZTFL1). Selective spatial transcriptomic analysis of lung biopsies from patients with COVID-19 shows the presence of signals associated with epithelial-mesenchymal transition (EMT), a viral response pathway that is regulated by LZTFL1. We conclude that pulmonary epithelial cells undergoing EMT, rather than immune cells, are likely responsible for the 3p21.31-associated risk. Since the 3p21.31 effect is conferred by a gain-of-function, LZTFL1 may represent a therapeutic target.

Investigation of the Basic Steps in the Chromosome Conformation Capture Procedure

A constellation of chromosome conformation capture methods (С-methods) are an important tool for biochemical analysis of the spatial interactions between DNA regions that are separated in the primary sequence. All these methods are based on the long sequence of basic steps of treating cells, nuclei, chromatin, and finally DNA, thus representing a significant technical challenge. Here, we present an in-depth study of the basic steps in the chromatin conformation capture procedure (3С), which was performed using Drosophila Schneider 2 cells as a model. We investigated the steps of cell lysis, nuclei washing, nucleoplasm extraction, chromatin treatment with SDS/Triton X-100, restriction enzyme digestion, chromatin ligation, reversion of cross-links, DNA extraction, treatment of a 3C library with RNases, and purification of the 3C library. Several options were studied, and optimal conditions were found. Our work contributes to the understanding of the 3C basic steps and provides a useful guide to the 3C procedure.

Activation of γ-globin gene expression by GATA1 and NF-Y in hereditary persistence of fetal hemoglobin

Hereditary persistence of fetal hemoglobin (HPFH) ameliorates β-hemoglobinopathies by inhibiting the developmental switch from γ-globin (HBG1/HBG2) to β-globin (HBB) gene expression. Some forms of HPFH are associated with γ-globin promoter variants that either disrupt binding motifs for transcriptional repressors or create new motifs for transcriptional activators. How these variants sustain γ-globin gene expression postnatally remains undefined. We mapped γ-globin promoter sequences functionally in erythroid cells harboring different HPFH variants. Those that disrupt a BCL11A repressor binding element induce γ-globin expression by facilitating the recruitment of transcription factors NF-Y to a nearby proximal CCAAT box and GATA1 to an upstream motif. The proximal CCAAT element becomes dispensable for HPFH variants that generate new binding motifs for activators NF-Y or KLF1, but GATA1 recruitment remains essential. Our findings define distinct mechanisms through which transcription factors and their cis-regulatory elements activate γ-globin expression in different forms of HPFH, some of which are being recreated by therapeutic genome editing.