The first enhancer in an enhancer chain safeguards subsequent enhancer-promoter contacts from a distance

Supervised learning of enhancer-promoter specificity based on genome-wide perturbation studies highlights areas for improvement in learning

MOTIVATION

Understanding the rules that govern enhancer-driven transcription remains a central unsolved problem in genomics. Now with multiple massively parallel enhancer perturbation assays published, there are enough data that we can utilize to learn to predict enhancer-promoter (EP) relationships in a data-driven manner.

RESULTS

We applied machine learning to one of the largest enhancer perturbation studies integrated with transcription factor (TF) and histone modification ChIP-seq. The results uncovered a discrepancy in the prediction of genome-wide data compared to data from targeted experiments. Relative strength of contact was important for prediction, confirming the basic principle of EP regulation. Novel features such as the density of the enhancers/promoters in the genomic region was found to be important, highlighting our lack of understanding on how other elements in the region contribute to the regulation. Several TF peaks were identified that improved the prediction by identifying the negatives and reducing False Positives. In summary, integrating genomic assays with enhancer perturbation studies increased the accuracy of the model, and provided novel insights into the understanding of enhancer-driven transcription.

AVAILABILITY AND IMPLEMENTATION

The trained models, data, and the source code are available at http://doi.org/10.5281/zenodo.11290386 and https://github.com/HanLabUNLV/sleps.

Building Haplotype-Resolved 3D Genome Maps of Chicken Skeletal Muscle

Haplotype-resolved 3D chromatin architecture related to allelic differences in avian skeletal muscle development has not been addressed so far, although chicken husbandry for meat consumption has been prevalent feature of cultures on every continent for more than thousands of years. Here, high-resolution Hi-C diploid maps (1.2-kb maximum resolution) are generated for skeletal muscle tissues in chicken across three developmental stages (embryonic day 15 to day 30 post-hatching). The sequence features governing spatial arrangement of chromosomes and characterize homolog pairing in the nucleus, are identified. Multi-scale characterization of chromatin reorganization between stages from myogenesis in the fetus to myofiber hypertrophy after hatching show concordant changes in transcriptional regulation by relevant signaling pathways. Further interrogation of parent-of-origin-specific chromatin conformation supported that genomic imprinting is absent in birds. This study also reveals promoter-enhancer interaction (PEI) differences between broiler and layer haplotypes in skeletal muscle development-related genes are related to genetic variation between breeds, however, only a minority of breed-specific variations likely contribute to phenotypic divergence in skeletal muscle potentially via allelic PEI rewiring. Beyond defining the haplotype-specific 3D chromatin architecture in chicken, this study provides a rich resource for investigating allelic regulatory divergence among chicken breeds.

Cohesin-Dependent Loop Extrusion: Molecular Mechanics and Role in Cell Physiology

The most prominent representatives of multisubunit SMC complexes, cohesin and condensin, are best known as structural components of mitotic chromosomes. It turned out that these complexes, as well as their bacterial homologues, are molecular motors, the ATP-dependent movement of these complexes along DNA threads leads to the formation of DNA loops. In recent years, we have witnessed an avalanche-like accumulation of data on the process of SMC dependent DNA looping, also known as loop extrusion. This review briefly summarizes the current understanding of the place and role of cohesin-dependent extrusion in cell physiology and presents a number of models describing the potential molecular mechanism of extrusion in a most compelling way. We conclude the review with a discussion of how the capacity of cohesin to extrude DNA loops may be mechanistically linked to its involvement in sister chromatid cohesion.

Haplotype-resolved 3D chromatin architecture of the hybrid pig

In diploid mammals, allele-specific three-dimensional (3D) genome architecture may lead to imbalanced gene expression. Through ultradeep in situ Hi-C sequencing of three representative somatic tissues (liver, skeletal muscle, and brain) from hybrid pigs generated by reciprocal crosses of phenotypically and physiologically divergent Berkshire and Tibetan pigs, we uncover extensive chromatin reorganization between homologous chromosomes across multiple scales. Haplotype-based interrogation of multi-omic data revealed the tissue dependence of 3D chromatin conformation, suggesting that parent-of-origin-specific conformation may drive gene imprinting. We quantify the effects of genetic variations and histone modifications on allelic differences of long-range promoter-enhancer contacts, which likely contribute to the phenotypic differences between the parental pig breeds. We also observe the fine structure of somatically paired homologous chromosomes in the pig genome, which has a functional implication genome-wide. This work illustrates how allele-specific chromatin architecture facilitates concomitant shifts in allele-biased gene expression, as well as the possible consequential phenotypic changes in mammals.

MotifHub: Detection of trans-acting DNA motif group with probabilistic modeling algorithm

BACKGROUND

Trans-acting factors are of special importance in transcription regulation, which is a group of proteins that can directly or indirectly recognize or bind to the 8-12 bp core sequence of cis-acting elements and regulate the transcription efficiency of target genes. The progressive development in high-throughput chromatin capture technology (e.g., Hi-C) enables the identification of chromatin-interacting sequence groups where trans-acting DNA motif groups can be discovered. The problem difficulty lies in the combinatorial nature of DNA sequence pattern matching and its underlying sequence pattern search space.

METHOD

Here, we propose to develop MotifHub for trans-acting DNA motif group discovery on grouped sequences. Specifically, the main approach is to develop probabilistic modeling for accommodating the stochastic nature of DNA motif patterns.

RESULTS

Based on the modeling, we develop global sampling techniques based on EM and Gibbs sampling to address the global optimization challenge for model fitting with latent variables. The results reflect that our proposed approaches demonstrate promising performance with linear time complexities.

CONCLUSION

MotifHub is a novel algorithm considering the identification of both DNA co-binding motif groups and trans-acting TFs. Our study paves the way for identifying hub TFs of stem cell development (OCT4 and SOX2) and determining potential therapeutic targets of prostate cancer (FOXA1 and MYC). To ensure scientific reproducibility and long-term impact, its matrix-algebra-optimized source code is released at http://bioinfo.cs.cityu.edu.hk/MotifHub.

STAT3-Mediated Promoter-Enhancer Interaction Up-Regulates Inhibitor of DNA Binding 1 (ID1) to Promote Colon Cancer Progression

BACKGROUND

High expression of inhibitor of DNA binding 1 (ID1) correlates with poor prognosis in colorectal cancer (CRC). Aberrant enhancer activation in regulating ID1 transcription is limited.

METHODS

Immunohistochemistry (IHC), quantitative RT-PCR (RT-qPCR) and Western blotting (WB) were used to determine the expression of ID1. CRISPR-Cas9 was used to generate ID1 or enhancer E1 knockout cell lines. Dual-luciferase reporter assay, chromosome conformation capture assay and ChIP-qPCR were used to determine the active enhancers of ID1. Cell Counting Kit 8, colony-forming, transwell assays and tumorigenicity in nude mice were used to investigate the biological functions of ID1 and enhancer E1.

RESULTS

Human CRC tissues and cell lines expressed a higher level of ID1 than normal controls. ID1 promoted CRC cell proliferation and colony formation. Enhancer E1 actively regulated ID1 promoter activity. Signal transducer and activator of transcription 3 (STAT3) bound to ID1 promoter and enhancer E1 to regulate their activity. The inhibitor of STAT3 Stattic attenuated ID1 promoter and enhancer E1 activity and the expression of ID1. Enhancer E1 knockout down-regulated ID1 expression level and cell proliferation in vitro and in vivo.

CONCLUSIONS

Enhancer E1 is positively regulated by STAT3 and contributes to the regulation of ID1 to promote CRC cell progression and might be a potential target for anti-CRC drug studies.

New insight into copy number variations of goat SMAD2 gene and their associations with litter size and semen quality

Copy number variations (CNV) contribute significantly to genetic variations. Numerous studies have shown that CNV affects phenotypic traits in livestock. The SMAD family member 2 (SMAD2) is a leading candidate gene in reproduction and has a crucial effect on litter size. Additionally, SMAD2 is also required for male reproduction and influences male germ cell development. However, there are no reports on investigating the effect of CNVs in the SMAD2 gene on reproductive traits in goat. Therefore, the goal of this study was to explore associations between CNV of the SMAD2 gene and litter size and semen quality in Shaanbei white cashmere (SBWC) goats. In this study, two CNVs within the SMAD2 were identified in 352 SBWC goats (50 males and 302 females). The association analysis revealed that only CNV2 was significantly associated with female goat first-born litter size (P = 3.59 × 10^-4), male semen concentration (P < 0.01), ejaculation volume, live sperm count, and sperm deformity rate (P < 0.05). In terms of phenotypic performance, the individuals with loss genotypes outperformed those with other genotypes. CNV1 and CNV2 genotype combinations containing their dominant genotypes were also associated with goat litter size (P = 1.7 × 10^-5), but no differences in semen quality were found. In summary, CNV2 of the SMAD2 gene is useful for molecular marker-assisted selection breeding, as it is associated with essential goat reproductive traits.

Heterogeneity of enhancers embodies shared and representative functional groups underlying developmental and cell type-specific gene regulation

While enhancers in a particular tissue coordinately fulfill regulatory functions, these functions are heterogeneous in nature and comprise of multiple enhancer subclasses and the associated regulatory mechanisms. In this work, we used multiple cell lines to identify enhancer subclasses linked to development, differentiation, and cellular identity. We found that enhancer functional heterogeneity during development encompasses subclasses of ubiquitous functions (11%), development specific regulatory activity (62%), and chromatin interactions (12%). In differentiated cell lines, ubiquitous enhancers (10%) stay active across multiple cell lines.They are accompanied by a large enhancer subclass (ranging from 33% to 63%) with functions specific to the corresponding lineage. The remaining enhancers (27-40%) establish regulatory chromatin structure and facilitate interactions of cell type-specific enhancers with their target promoters. In addition to specialized functions of cell type-specific enhancers, we show that proper accounting of enhancer heterogeneity leads to a 10% increase in accuracy of enhancer classification, which significantly improves the modeling of enhancers and identification of underlying regulatory mechanisms. In summary, our observations suggest that although cell type-specific enhancers are heterogeneous and coordinate different regulatory programs, enhancers from different cell lines maintain common categories of functional groups across developmental and differentiation stages, indicating a higher order rule followed by enhancer-gene regulation.

Application of the 3C Method to Study the Developmental Genes in Drosophila Larvae

A transition from one developmental stage to another is accompanied by activation of developmental programs and corresponding gene ensembles. Changes in the spatial conformation of the corresponding loci are associated with this activation and can be investigated with the help of the Chromosome Conformation Capture (3C) methodology. Application of 3C to specific developmental stages is a sophisticated task. Here, we describe the use of the 3C method to study the spatial organization of developmental loci in Drosophila larvae. We critically analyzed the existing protocols and offered our own solutions and the optimized protocol to overcome limitations. To demonstrate the efficiency of our procedure, we studied the spatial organization of the developmental locus Dad in 3rd instar Drosophila larvae. Differences in locus conformation were found between embryonic cells and living wild-type larvae. We also observed the establishment of novel regulatory interactions in the presence of an adjacent transgene upon activation of its expression in larvae. Our work fills the gap in the application of the 3C method to Drosophila larvae and provides a useful guide for establishing 3C on an animal model.

Dynamic 3D genome reorganization during development and metabolic stress of the porcine liver

Liver development is a complex process that is regulated by a series of signaling pathways. Three-dimensional (3D) chromatin architecture plays an important role in transcriptional regulation; nonetheless, its dynamics and role in the rapid transition of core liver functions during development and obesity-induced metabolic stress remain largely unexplored. To investigate the dynamic chromatin architecture during liver development and under metabolic stress, we generated high-resolution maps of chromatin architecture for porcine livers across six major developmental stages (from embryonic day 38 to the adult stage) and under a high-fat diet-induced obesity. The characteristically loose chromatin architecture supports a highly plastic genome organization during early liver development, which fundamentally contributes to the rapid functional transitions in the liver after birth. We reveal the multi-scale reorganization of chromatin architecture and its influence on transcriptional regulation of critical signaling processes during liver development, and show its close association with transition in hepatic functions (i.e., from hematopoiesis in the fetus to metabolism and immunity after birth). The limited changes in chromatin structure help explain the observed metabolic adaptation to excessive energy intake in pigs. These results provide a global overview of chromatin architecture dynamics associated with the transition of physiological liver functions between prenatal development and postnatal maturation, and a foundational resource that allows for future in-depth functional characterization.

Duplications disrupt chromatin architecture and rewire GPR101-enhancer communication in X-linked acrogigantism

X-linked acrogigantism (X-LAG) is the most severe form of pituitary gigantism and is characterized by aggressive growth hormone (GH)-secreting pituitary tumors that occur in early childhood. X-LAG is associated with chromosome Xq26.3 duplications (the X-LAG locus typically includes VGLL1, CD40LG, ARHGEF6, RBMX, and GPR101) that lead to massive pituitary tumoral expression of GPR101, a novel regulator of GH secretion. The mechanism by which the duplications lead to marked pituitary misexpression of GPR101 alone was previously unclear. Using Hi-C and 4C-seq, we characterized the normal chromatin structure at the X-LAG locus. We showed that GPR101 is located within a topologically associating domain (TAD) delineated by a tissue-invariant border that separates it from centromeric genes and regulatory sequences. Next, using 4C-seq with GPR101, RBMX, and VGLL1 viewpoints, we showed that the duplications in multiple X-LAG-affected individuals led to ectopic interactions that crossed the invariant TAD border, indicating the existence of a similar and consistent mechanism of neo-TAD formation in X-LAG. We then identified several pituitary active cis-regulatory elements (CREs) within the neo-TAD and demonstrated in vitro that one of them significantly enhanced reporter gene expression. At the same time, we showed that the GPR101 promoter permits the incorporation of new regulatory information. Our results indicate that X-LAG is a TADopathy of the endocrine system in which Xq26.3 duplications disrupt the local chromatin architecture forming a neo-TAD. Rewiring GPR101-enhancer interaction within the new regulatory unit is likely to cause the high levels of aberrant expression of GPR101 in pituitary tumors caused by X-LAG.

CT-FOCS: a novel method for inferring cell type-specific enhancer–promoter maps

Spatiotemporal gene expression patterns are governed to a large extent by the activity of enhancer elements, which engage in physical contacts with their target genes. Identification of enhancer–promoter (EP) links that are functional only in a specific subset of cell types is a key challenge in understanding gene regulation. We introduce CT-FOCS (cell type FOCS), a statistical inference method that uses linear mixed effect models to infer EP links that show marked activity only in a single or a small subset of cell types out of a large panel of probed cell types. Analyzing 808 samples from FANTOM5, covering 472 cell lines, primary cells and tissues, CT-FOCS inferred such EP links more accurately than recent state-of-the-art methods. Furthermore, we show that strictly cell type-specific EP links are very uncommon in the human genome.

Transcriptional enhancers and their communication with gene promoters

Transcriptional enhancers play a key role in the initiation and maintenance of gene expression programmes, particularly in metazoa. How these elements control their target genes in the right place and time is one of the most pertinent questions in functional genomics, with wide implications for most areas of biology. Here, we synthesise classic and recent evidence on the regulatory logic of enhancers, including the principles of enhancer organisation, factors that facilitate and delimit enhancer-promoter communication, and the joint effects of multiple enhancers. We show how modern approaches building on classic insights have begun to unravel the complexity of enhancer-promoter relationships, paving the way towards a quantitative understanding of gene control.

Mechanisms of enhancer action: the known and the unknown

Differential gene expression mechanisms ensure cellular differentiation and plasticity to shape ontogenetic and phylogenetic diversity of cell types. A key regulator of differential gene expression programs are the enhancers, the gene-distal cis-regulatory sequences that govern spatiotemporal and quantitative expression dynamics of target genes. Enhancers are widely believed to physically contact the target promoters to effect transcriptional activation. However, our understanding of the full complement of regulatory proteins and the definitive mechanics of enhancer action is incomplete. Here, we review recent findings to present some emerging concepts on enhancer action and also outline a set of outstanding questions.

Cohesin-Dependent and -Independent Mechanisms Mediate Chromosomal Contacts between Promoters and Enhancers

It is currently assumed that 3D chromosomal organization plays a central role in transcriptional control. However, depletion of cohesin and CTCF affects the steady-state levels of only a minority of transcripts. Here, we use high-resolution Capture Hi-C to interrogate the dynamics of chromosomal contacts of all annotated human gene promoters upon degradation of cohesin and CTCF. We show that a majority of promoter-anchored contacts are lost in these conditions, but many contacts with distinct properties are maintained, and some new ones are gained. The rewiring of contacts between promoters and active enhancers upon cohesin degradation associates with rapid changes in target gene transcription as detected by SLAM sequencing (SLAM-seq). These results provide a mechanistic explanation for the limited, but consistent, effects of cohesin and CTCF depletion on steady-state transcription and suggest the existence of both cohesin-dependent and -independent mechanisms of enhancer-promoter pairing.

Evolution of Regulated Transcription

The genomes of all organisms abound with various cis-regulatory elements, which control gene activity. Transcriptional enhancers are a key group of such elements in eukaryotes and are DNA regions that form physical contacts with gene promoters and precisely orchestrate gene expression programs. Here, we follow gradual evolution of this regulatory system and discuss its features in different organisms. In eubacteria, an enhancer-like element is often a single regulatory element, is usually proximal to the core promoter, and is occupied by one or a few activators. Activation of gene expression in archaea is accompanied by the recruitment of an activator to several enhancer-like sites in the upstream promoter region. In eukaryotes, activation of expression is accompanied by the recruitment of activators to multiple enhancers, which may be distant from the core promoter, and the activators act through coactivators. The role of the general DNA architecture in transcription control increases in evolution. As a whole, it can be seen that enhancers of multicellular eukaryotes evolved from the corresponding prototypic enhancer-like regulatory elements with the gradually increasing genome size of organisms.

Transcriptional activation during cell reprogramming correlates with the formation of 3D open chromatin hubs

Chromosome structure is a crucial regulatory factor for a wide range of nuclear processes. Chromosome conformation capture (3C)-based experiments combined with computational modelling are pivotal for unveiling 3D chromosome structure. Here, we introduce TADdyn, a tool that integrates time-course 3C data, restraint-based modelling, and molecular dynamics to simulate the structural rearrangements of genomic loci in a completely data-driven way. We apply TADdyn on in situ Hi-C time-course experiments studying the reprogramming of murine B cells to pluripotent cells, and characterize the structural rearrangements that take place upon changes in the transcriptional state of 21 genomic loci of diverse expression dynamics. By measuring various structural and dynamical properties, we find that during gene activation, the transcription starting site contacts with open and active regions in 3D chromatin domains. We propose that these 3D hubs of open and active chromatin may constitute a general feature to trigger and maintain gene transcription.