Deciphering the multi-scale, quantitative cis-regulatory code

Pre-hypertrophic chondrogenic enhancer landscape of limb and axial skeleton development

Chondrocyte differentiation controls skeleton development and stature. Here we provide a comprehensive map of chondrocyte-specific enhancers and show that they provide a mechanistic framework through which non-coding genetic variants can influence skeletal development and human stature. Working with fetal chondrocytes isolated from mice bearing a Col2a1 fluorescent regulatory sensor, we identify 780 genes and 2'704 putative enhancers specifically active in chondrocytes using a combination of RNA-seq, ATAC-seq and H3K27ac ChIP-seq. Most of these enhancers (74%) show pan-chondrogenic activity, with smaller populations being restricted to limb (18%) or trunk (8%) chondrocytes only. Notably, genetic variations overlapping these enhancers better explain height differences than those overlapping non-chondrogenic enhancers. Finally, targeted deletions of identified enhancers at the Fgfr3, Col2a1, Hhip and, Nkx3-2 loci confirm their role in regulating cognate genes. This enhancer map provides a framework for understanding how genes and non-coding variations influence bone development and diseases.

Comparative cofactor screens show the influence of transactivation domains and core promoters on the mechanisms of transcription

Eukaryotic transcription factors (TFs) activate gene expression by recruiting cofactors to promoters. However, the relationships between TFs, promoters and their associated cofactors remain poorly understood. Here we combine GAL4-transactivation assays with comparative CRISPR-Cas9 screens to identify the cofactors used by nine different TFs and core promoters in human cells. Using this dataset, we associate TFs with cofactors, classify cofactors as ubiquitous or specific and discover transcriptional co-dependencies. Through a reductionistic, comparative approach, we demonstrate that TFs do not display discrete mechanisms of activation. Instead, each TF depends on a unique combination of cofactors, which influences distinct steps in transcription. By contrast, the influence of core promoters appears relatively discrete. Different promoter classes are constrained by either initiation or pause-release, which influences their dynamic range and compatibility with cofactors. Overall, our comparative cofactor screens characterize the interplay between TFs, cofactors and core promoters, identifying general principles by which they influence transcription.

A computational approach for deciphering the interactions between proximal and distal gene regulators in GC B-cell response

Delineating the intricate interplay between promoter-proximal and -distal regulators is crucial for understanding the function of transcriptional mediator complexes implicated in the regulation of gene expression. The present study aimed to develop a computational method for accurately modeling the spatial proximal and distal regulatory interactions. Our method combined regression-based models to identify key regulators through gene expression prediction and a graph-embedding approach to detect coregulated genes. This approach enabled a detailed investigation of the gene regulatory mechanisms for germinal center B cells, accompanied by dramatic rearrangements of the genome structure. We found that while the promoter-proximal regulatory elements were the principal regulators of gene expression, the distal regulators fine-tuned transcription. Moreover, our approach unveiled the presence of modular regulators, such as cofactors and proximal/distal transcription factors, which were co-expressed with their target genes. Some of these modules exhibited abnormal expression patterns in lymphoma. These findings suggest that the dysregulation of interactions between transcriptional and architectural factors is associated with chromatin reorganization failure, which may increase the risk of malignancy. Therefore, our computational approach helps decipher the transcriptional cis-regulatory code spatially interacting.

A cis-regulatory module underlies retinal ganglion cell genesis and axonogenesis

Atoh7 is transiently expressed in retinal progenitor cells (RPCs) and is required for retinal ganglion cell (RGC) differentiation. In humans, a deletion in a distal non-coding regulatory region upstream of ATOH7 is associated with optic nerve atrophy and blindness. Here, we functionally interrogate the significance of the Atoh7 regulatory landscape to retinogenesis in mice. Deletion of the Atoh7 enhancer structure leads to RGC deficiency, optic nerve hypoplasia, and retinal blood vascular abnormalities, phenocopying inactivation of Atoh7. Further, loss of the Atoh7 remote enhancer impacts ipsilaterally projecting RGCs and disrupts proper axonal projections to the visual thalamus. Deletion of the Atoh7 remote enhancer is also associated with the dysregulation of axonogenesis genes, including the derepression of the axon repulsive cue Robo3. Our data provide insights into how Atoh7 enhancer elements function to promote RGC development and optic nerve formation and highlight a key role of Atoh7 in the transcriptional control of axon guidance molecules.

Transfer learning reveals sequence determinants of the quantitative response to transcription factor dosage

Deep learning approaches have made significant advances in predicting cell type-specific chromatin patterns from the identity and arrangement of transcription factor (TF) binding motifs. However, most models have been applied in unperturbed contexts, precluding a predictive understanding of how chromatin state responds to TF perturbation. Here, we used transfer learning to train and interpret deep learning models that use DNA sequence to predict, with accuracy approaching experimental reproducibility, how the concentration of two dosage-sensitive TFs (TWIST1, SOX9) affects regulatory element (RE) chromatin accessibility in facial progenitor cells. High-affinity motifs that allow for heterotypic TF co-binding and are concentrated at the center of REs buffer against quantitative changes in TF dosage and strongly predict unperturbed accessibility. In contrast, motifs with low-affinity or homotypic binding distributed throughout REs lead to sensitive responses with minimal contributions to unperturbed accessibility. Both buffering and sensitizing features show signatures of purifying selection. We validated these predictive sequence features using reporter assays and showed that a biophysical model of TF-nucleosome competition can explain the sensitizing effect of low-affinity motifs. Our approach of combining transfer learning and quantitative measurements of the chromatin response to TF dosage therefore represents a powerful method to reveal additional layers of the cis-regulatory code.

YY1 binding is a gene-intrinsic barrier to Xist-mediated gene silencing

X chromosome inactivation (XCI) in mammals is mediated by Xist RNA which functions in cis to silence genes on a single X chromosome in XX female cells, thereby equalising levels of X-linked gene expression relative to XY males. XCI progresses over a period of several days, with some X-linked genes silencing faster than others. The chromosomal location of a gene is an important determinant of silencing rate, but uncharacterised gene-intrinsic features also mediate resistance or susceptibility to silencing. In this study, we examine mouse embryonic stem cell lines with an inducible Xist allele (iXist-ChrX mESCs) and integrate allele-specific data of gene silencing and decreasing inactive X (Xi) chromatin accessibility over time courses of Xist induction with cellular differentiation. Our analysis reveals that motifs bound by the transcription factor YY1 are associated with persistently accessible regulatory elements, including many promoters and enhancers of slow-silencing genes. We further show that YY1 is evicted relatively slowly from target sites on Xi, and that silencing of X-linked genes is increased upon YY1 degradation. Together our results suggest that YY1 acts as a barrier to Xist-mediated silencing until the late stages of the XCI process.

Nucleosome-bound NR5A2 structure reveals pioneer factor mechanism by DNA minor groove anchor competition

Gene expression during natural and induced reprogramming is controlled by pioneer transcription factors that initiate transcription from closed chromatin. Nr5a2 is a key pioneer factor that regulates zygotic genome activation in totipotent embryos, pluripotency in embryonic stem cells and metabolism in adult tissues, but the mechanism of its pioneer activity remains poorly understood. Here, we present a cryo-electron microscopy structure of human NR5A2 bound to a nucleosome. The structure shows that the conserved carboxy-terminal extension (CTE) loop of the NR5A2 DNA-binding domain competes with a DNA minor groove anchor of the nucleosome and releases entry-exit site DNA. Mutational analysis showed that NR5A2 D159 of the CTE is dispensable for DNA binding but required for stable nucleosome association and persistent DNA 'unwrapping'. These findings suggest that NR5A2 belongs to an emerging class of pioneer factors that can use DNA minor groove anchor competition to destabilize nucleosomes and facilitate gene expression during reprogramming.

Emergence of enhancers at late DNA replicating regions

Enhancers are fast-evolving genomic sequences that control spatiotemporal gene expression patterns. By examining enhancer turnover across mammalian species and in multiple tissue types, we uncover a relationship between the emergence of enhancers and genome organization as a function of germline DNA replication time. While enhancers are most abundant in euchromatic regions, enhancers emerge almost twice as often in late compared to early germline replicating regions, independent of transposable elements. Using a deep learning sequence model, we demonstrate that new enhancers are enriched for mutations that alter transcription factor (TF) binding. Recently evolved enhancers appear to be mostly neutrally evolving and enriched in eQTLs. They also show more tissue specificity than conserved enhancers, and the TFs that bind to these elements, as inferred by binding sequences, also show increased tissue-specific gene expression. We find a similar relationship with DNA replication time in cancer, suggesting that these observations may be time-invariant principles of genome evolution. Our work underscores that genome organization has a profound impact in shaping mammalian gene regulation.

Circadian regulation of stereotypic chromatin conformations at enhancers

Cooperation between the circadian transcription factor (TF) CLOCK:BMAL1 and other TFs at cis-regulatory elements (CREs) is critical to daily rhythms of transcription. Yet, the modalities of this cooperation are unclear. Here, we analyzed the co-binding of multiple TFs on single DNA molecules in mouse liver using single molecule footprinting (SMF). We found that SMF reads clustered in stereotypic chromatin states that reflect distinguishable organization of TFs and nucleosomes, and that were remarkably conserved between all samples. DNA protection at CLOCK:BMAL1 binding motif (E-box) varied between CREs, from E-boxes being solely bound by CLOCK:BMAL1 to situations where other TFs competed with CLOCK:BMAL1 for E-box binding. SMF also uncovered CLOCK:BMAL1 cooperative binding at E-boxes separated by 250 bp, which structurally altered the CLOCK:BMAL1-DNA interface. Importantly, we discovered multiple nucleosomes with E-boxes at entry/exit sites that were removed upon CLOCK:BMAL1 DNA binding, thereby promoting the formation of open chromatin states that facilitate DNA binding of other TFs and that were associated with rhythmic transcription. These results demonstrate the utility of SMF for studying how CLOCK:BMAL1 and other TFs regulate stereotypical chromatin states at CREs to promote transcription.

Yeast heterochromatin stably silences only weak regulatory elements by altering burst duration

Transcriptional silencing in Saccharomyces cerevisiae involves the generation of a chromatin state that stably represses transcription. Using multiple reporter assays, a diverse set of upstream activating sequence enhancers and core promoters were investigated for their susceptibility to silencing. We show that heterochromatin stably silences only weak and stress-induced regulatory elements but is unable to stably repress housekeeping gene regulatory elements, and the partial repression of these elements did not result in bistable expression states. Permutation analysis of enhancers and promoters indicates that both elements are targets of repression. Chromatin remodelers help specific regulatory elements to resist repression, most probably by altering nucleosome mobility and changing transcription burst duration. The strong enhancers/promoters can be repressed if silencer-bound Sir1 is increased. Together, our data suggest that the heterochromatic locus has been optimized to stably silence the weak mating-type gene regulatory elements but not strong housekeeping gene regulatory sequences.

Physical and functional interaction among Irf8 enhancers during dendritic cell differentiation

The production of type 1 conventional dendritic cells (cDC1s) requires high expression of the transcription factor IRF8. Three enhancers at the Irf8 3' region function in a differentiation stage-specific manner. However, whether and how these enhancers interact physically and functionally remains unclear. Here, we show that the Irf8 3' enhancers directly interact with each other and contact the Irf8 gene body during cDC1 differentiation. The +56 kb enhancer, which functions from multipotent progenitor stages, activates the other 3' enhancers through an IRF8-dependent transcription factor program, that is, in trans. Then, the +32 kb enhancer, which operates in cDC1-committed cells, reversely acts in cis on the other 3' enhancers to maintain the high expression of Irf8. Indeed, mice with compound heterozygous deletion of the +56 and +32 kb enhancers are unable to generate cDC1s. These results illustrate how multiple enhancers cooperate to induce a lineage-determining transcription factor gene during cell differentiation.

Epigenomic insights into common human disease pathology

The epigenome-the chemical modifications and chromatin-related packaging of the genome-enables the same genetic template to be activated or repressed in different cellular settings. This multi-layered mechanism facilitates cell-type specific function by setting the local sequence and 3D interactive activity level. Gene transcription is further modulated through the interplay with transcription factors and co-regulators. The human body requires this epigenomic apparatus to be precisely installed throughout development and then adequately maintained during the lifespan. The causal role of the epigenome in human pathology, beyond imprinting disorders and specific tumour suppressor genes, was further brought into the spotlight by large-scale sequencing projects identifying that mutations in epigenomic machinery genes could be critical drivers in both cancer and developmental disorders. Abrogation of this cellular mechanism is providing new molecular insights into pathogenesis. However, deciphering the full breadth and implications of these epigenomic changes remains challenging. Knowledge is accruing regarding disease mechanisms and clinical biomarkers, through pathogenically relevant and surrogate tissue analyses, respectively. Advances include consortia generated cell-type specific reference epigenomes, high-throughput DNA methylome association studies, as well as insights into ageing-related diseases from biological 'clocks' constructed by machine learning algorithms. Also, 3rd-generation sequencing is beginning to disentangle the complexity of genetic and DNA modification haplotypes. Cell-free DNA methylation as a cancer biomarker has clear clinical utility and further potential to assess organ damage across many disorders. Finally, molecular understanding of disease aetiology brings with it the opportunity for exact therapeutic alteration of the epigenome through CRISPR-activation or inhibition.

Rapid profiling of transcription factor-cofactor interaction networks reveals principles of epigenetic regulation

Transcription factor (TF)-cofactor (COF) interactions define dynamic, cell-specific networks that govern gene expression; however, these networks are understudied due to a lack of methods for high-throughput profiling of DNA-bound TF-COF complexes. Here we describe the Cofactor Recruitment (CoRec) method for rapid profiling of cell-specific TF-COF complexes. We define a lysine acetyltransferase (KAT)-TF network in resting and stimulated T cells. We find promiscuous recruitment of KATs for many TFs and that 35% of KAT-TF interactions are condition specific. KAT-TF interactions identify NF-κB as a primary regulator of acutely induced H3K27ac. Finally, we find that heterotypic clustering of CBP/P300-recruiting TFs is a strong predictor of total promoter H3K27ac. Our data supports clustering of TF sites that broadly recruit KATs as a mechanism for widespread co-occurring histone acetylation marks. CoRec can be readily applied to different cell systems and provides a powerful approach to define TF-COF networks impacting chromatin state and gene regulation.

Visualization of the Association of Dimeric Protein Complexes on Specific Enhancers in the Salivary Gland Nuclei of Drosophila Larva

Transcription factors (TFs) regulate gene expression by recognizing specific target enhancers in the genome. The DNA-binding and regulatory activity of TFs depend on the presence of additional protein partners, leading to the formation of versatile and dynamic multimeric protein complexes. Visualizing these protein-protein interactions (PPIs) in the nucleus is key for decrypting the molecular cues underlying TF specificity in vivo. Over the last few years, Bimolecular Fluorescence Complementation (BiFC) has been developed in several model systems and applied in the analysis of different types of PPIs. In particular, BiFC has been applied when analyzing PPIs with hundreds of TFs in the nucleus of live Drosophila embryos. However, the visualization of PPIs at the level of specific target enhancers or genomic regions of interest awaits the advent of DNA-labelling methods that can be coupled with BiFC. Here, we present a novel experimental strategy that we have called BiFOR and that is based on the coupling of BiFC with the bacterial ANCHOR DNA-labelling system. We demonstrate that BiFOR enables the precise quantification of the enrichment of specific dimeric protein complexes on target enhancers in Drosophila salivary gland nuclei. Given its versatility and sensitivity, BiFOR could be applied more widely to other tissues during Drosophila development. Our work sets up the experimental basis for future applications of this strategy.

SWI/SNF-dependent genes are defined by their chromatin landscape

SWI/SNF complexes are evolutionarily conserved, ATP-dependent chromatin remodeling machines. Here, we characterize the features of SWI/SNF-dependent genes using BRM014, an inhibitor of the ATPase activity of the complexes. We find that SWI/SNF activity is required to maintain chromatin accessibility and nucleosome occupancy for most enhancers but not for most promoters. SWI/SNF activity is needed for expression of genes with low to medium levels of expression that have promoters with (1) low chromatin accessibility, (2) low levels of active histone marks, (3) high H3K4me1/H3K4me3 ratio, (4) low nucleosomal phasing, and (5) enrichment in TATA-box motifs. These promoters are mostly occupied by the canonical Brahma-related gene 1/Brahma-associated factor (BAF) complex. These genes are surrounded by SWI/SNF-dependent enhancers and mainly encode signal transduction, developmental, and cell identity genes (with almost no housekeeping genes). Machine-learning models trained with different chromatin characteristics of promoters and their surrounding regulatory regions indicate that the chromatin landscape is a determinant for establishing SWI/SNF dependency.

PACS allows comprehensive dissection of multiple factors governing chromatin accessibility from snATAC-seq data

Single nucleus ATAC-seq (snATAC-seq) experimental designs have become increasingly complex with multiple factors that might affect chromatin accessibility, including genotype, cell type, tissue of origin, sample location, batch, etc., whose compound effects are difficult to test by existing methods. In addition, current snATAC-seq data present statistical difficulties due to their sparsity and variations in individual sequence capture. To address these problems, we present a zero-adjusted statistical model, Probability model of Accessible Chromatin of Single cells (PACS), that can allow complex hypothesis testing of factors that affect accessibility while accounting for sparse and incomplete data. For differential accessibility analysis, PACS controls the false positive rate and achieves on average a 17% to 122% higher power than existing tools. We demonstrate the effectiveness of PACS through several analysis tasks including supervised cell type annotation, compound hypothesis testing, batch effect correction, and spatiotemporal modeling. We apply PACS to several datasets from a variety of tissues and show its ability to reveal previously undiscovered insights in snATAC-seq data.

Interpreting cis-regulatory mechanisms from genomic deep neural networks using surrogate models

Deep neural networks (DNNs) have greatly advanced the ability to predict genome function from sequence. Interpreting genomic DNNs in terms of biological mechanisms, however, remains difficult. Here we introduce SQUID, a genomic DNN interpretability framework based on surrogate modeling. SQUID approximates genomic DNNs in user-specified regions of sequence space using surrogate models, i.e., simpler models that are mechanistically interpretable. Importantly, SQUID removes the confounding effects that nonlinearities and heteroscedastic noise in functional genomics data can have on model interpretation. Benchmarking analysis on multiple genomic DNNs shows that SQUID, when compared to established interpretability methods, identifies motifs that are more consistent across genomic loci and yields improved single-nucleotide variant-effect predictions. SQUID also supports surrogate models that quantify epistatic interactions within and between cis-regulatory elements. SQUID thus advances the ability to mechanistically interpret genomic DNNs.

Cell-free genomics: transcription factor interactions in reconstituted naïve embryonic chromatin

Extracts from Drosophila preblastoderm embryos (DREX) form the basis of a powerful in vitro chromatin reconstitution system that assembles entire genomes into complex chromatin with physiological nucleosome spacing and polymer condensation. As the zygotic genome has not yet been activated in preblastoderm embryos, the reconstitution extract lacks endogenous transcription factors (TFs) and the RNA polymerase machinery. At the same time, it contains high levels of ATP-dependent nucleosome sliding enzymes that render the reconstituted chromatin dynamic. The naïve chromatin can be used to determine the intrinsic DNA binding properties of exogenous, usually recombinant TFs (or DNA binding proteins in general) in a complex chromatin context. Recent applications of the system include the description of cooperation and competition of Drosophila pioneer TFs for composite binding sites, and the characterization of nucleosome interactions of mammalian pioneer TFs in the heterologous system.

Unraveling the mechanisms of NK cell dysfunction in aging and Alzheimer's disease: insights from GWAS and single-cell transcriptomics

Background

Aging is an important factor in the development of Alzheimer's disease (AD). The senescent cells can be recognized and removed by NK cells. However, NK cell function is gradually inactivated with age. Therefore, this study used senescence as an entry point to investigate how NK cells affect AD.

Methods

The study validated the correlation between cognition and aging through a prospective cohort of the National Health and Nutrition Examination Survey database. A cellular trajectory analysis of the aging population was performed using single-cell nuclear transcriptome sequencing data from patients with AD and different ages. The genome-wide association study (GWAS) cohort of AD patients was used as the outcome event, and the expression quantitative trait locus was used as an instrumental variable. Causal associations between genes and AD were analyzed by bidirectional Mendelian randomization (MR) and co-localization. Finally, clinical cohorts were constructed to validate the expression of key genes.

Results

A correlation between cognition and aging was demonstrated using 2,171 older adults over 60 years of age. Gene regulation analysis revealed that most of the highly active transcription factors were concentrated in the NK cell subpopulation of AD. NK cell trajectories were constructed for different age populations. MR and co-localization analyses revealed that CHD6 may be one of the factors influencing AD.

Conclusion

We explored different levels of AD and aging from population cohorts, single-cell data, and GWAS cohorts and found that there may be some correlations of NK cells between aging and AD. It also provides some basis for potential causation.

Exploring noncoding variants in genetic diseases: from detection to functional insights

Previous studies on genetic diseases predominantly focused on protein-coding variations, overlooking the vast noncoding regions in the human genome. The development of high-throughput sequencing technologies and functional genomics tools has enabled the systematic identification of functional noncoding variants. These variants can impact gene expression, regulation, and chromatin conformation, thereby contributing to disease pathogenesis. Understanding the mechanisms that underlie the impact of noncoding variants on genetic diseases is indispensable for the development of precisely targeted therapies and the implementation of personalized medicine strategies. The intricacies of noncoding regions introduce a multitude of challenges and research opportunities. In this review, we introduce a spectrum of noncoding variants involved in genetic diseases, along with research strategies and advanced technologies for their precise identification and in-depth understanding of the complexity of the noncoding genome. We will delve into the research challenges and propose potential solutions for unraveling the genetic basis of rare and complex diseases.

DNA-guided transcription factor cooperativity shapes face and limb mesenchyme

Transcription factors (TFs) can define distinct cellular identities despite nearly identical DNA-binding specificities. One mechanism for achieving regulatory specificity is DNA-guided TF cooperativity. Although in vitro studies suggest that it may be common, examples of such cooperativity remain scarce in cellular contexts. Here, we demonstrate how "Coordinator," a long DNA motif composed of common motifs bound by many basic helix-loop-helix (bHLH) and homeodomain (HD) TFs, uniquely defines the regulatory regions of embryonic face and limb mesenchyme. Coordinator guides cooperative and selective binding between the bHLH family mesenchymal regulator TWIST1 and a collective of HD factors associated with regional identities in the face and limb. TWIST1 is required for HD binding and open chromatin at Coordinator sites, whereas HD factors stabilize TWIST1 occupancy at Coordinator and titrate it away from HD-independent sites. This cooperativity results in the shared regulation of genes involved in cell-type and positional identities and ultimately shapes facial morphology and evolution.

TFCheckpoint database update, a cross-referencing system for transcription factors from human, mouse and rat

Prior knowledge about DNA-binding transcription factors (dbTFs), transcription co-regulators (coTFs) and general transcriptional factors (GTFs) is crucial for the study and understanding of the regulation of transcription. This is reflected by the many publications and database resources describing knowledge about TFs. We previously launched the TFCheckpoint database, an integrated resource focused on human, mouse and rat dbTFs, providing users access to a comprehensive overview of these proteins. Here, we describe TFCheckpoint 2.0 (https://www.tfcheckpoint.org/index.php), comprising 13 collections of dbTFs, coTFs and GTFs. TFCheckpoint 2.0 provides an easy and versatile cross-referencing system for users to view and download collections that may otherwise be cumbersome to find, compare and retrieve.

Hold out the genome: a roadmap to solving the cis-regulatory code

Gene expression is regulated by transcription factors that work together to read cis-regulatory DNA sequences. The 'cis-regulatory code' - how cells interpret DNA sequences to determine when, where and how much genes should be expressed - has proven to be exceedingly complex. Recently, advances in the scale and resolution of functional genomics assays and machine learning have enabled substantial progress towards deciphering this code. However, the cis-regulatory code will probably never be solved if models are trained only on genomic sequences; regions of homology can easily lead to overestimation of predictive performance, and our genome is too short and has insufficient sequence diversity to learn all relevant parameters. Fortunately, randomly synthesized DNA sequences enable testing a far larger sequence space than exists in our genomes, and designed DNA sequences enable targeted queries to maximally improve the models. As the same biochemical principles are used to interpret DNA regardless of its source, models trained on these synthetic data can predict genomic activity, often better than genome-trained models. Here we provide an outlook on the field, and propose a roadmap towards solving the cis-regulatory code by a combination of machine learning and massively parallel assays using synthetic DNA.

Stage-specific dynamic reorganization of genome topology shapes transcriptional neighborhoods in developing human retinal organoids

We have generated a high-resolution Hi-C map of developing human retinal organoids to elucidate spatiotemporal dynamics of genomic architecture and its relationship with gene expression patterns. We demonstrate progressive stage-specific alterations in DNA topology and correlate these changes with transcription of cell-type-restricted gene markers during retinal differentiation. Temporal Hi-C reveals a shift toward A compartment for protein-coding genes and B compartment for non-coding RNAs, displaying high and low expression, respectively. Notably, retina-enriched genes are clustered near lost boundaries of topologically associated domains (TADs), and higher-order assemblages (i.e., TAD cliques) localize in active chromatin regions with binding sites for eye-field transcription factors. These genes gain chromatin contacts at their transcription start site as organoid differentiation proceeds. Our study provides a global view of chromatin architecture dynamics associated with diversification of cell types during retinal development and serves as a foundational resource for in-depth functional investigations of retinal developmental traits.

Sequential deregulation of histone marks, chromatin accessibility and gene expression in response to PROTAC-induced degradation of ASH2L

The trithorax protein ASH2L is essential for organismal and tissue development. As a subunit of COMPASS/KMT2 complexes, ASH2L is necessary for methylation of histone H3 lysine 4 (H3K4). Mono- and tri-methylation at this site mark active enhancers and promoters, respectively, although the functional relevance of H3K4 methylation is only partially understood. ASH2L has a long half-life, which results in a slow decrease upon knockout. This has made it difficult to define direct consequences. To overcome this limitation, we employed a PROTAC system to rapidly degrade ASH2L and address direct effects. ASH2L loss resulted in inhibition of proliferation of mouse embryo fibroblasts. Shortly after ASH2L degradation H3K4me3 decreased with its half-life varying between promoters. Subsequently, H3K4me1 increased at promoters and decreased at some enhancers. H3K27ac and H3K27me3, histone marks closely linked to H3K4 methylation, were affected with considerable delay. In parallel, chromatin compaction increased at promoters. Of note, nascent gene transcription was not affected early but overall RNA expression was deregulated late after ASH2L loss. Together, these findings suggest that downstream effects are ordered but relatively slow, despite the rapid loss of ASH2L and inactivation of KMT2 complexes. It appears that the systems that control gene transcription are well buffered and strong effects are only beginning to unfold after considerable delay.

'Enhancing' skeletal muscle and stem cells in three-dimensions: genome regulation of skeletal muscle in development and disease

The noncoding genome imparts important regulatory control over gene expression. In particular, gene enhancers represent a critical layer of control that integrates developmental and differentiation signals outside the cell into transcriptional outputs inside the cell. Recently, there has been an explosion in genomic techniques to probe enhancer control, function, and regulation. How enhancers are regulated and integrate signals in stem cell development and differentiation is largely an open question. In this review, we focus on the role gene enhancers play in muscle stem cell specification, differentiation, and progression. We pay specific attention toward the identification of muscle-specific enhancers, the binding of transcription factors to these enhancers, and how enhancers communicate to their target genes via three-dimensional looping.

Topology regulatory elements: From shaping genome architecture to gene regulation

The importance of 3D genome topology in the control of gene expression is becoming increasingly apparent, while regulatory mechanisms remain incompletely understood. Several recent studies have identified architectural elements that influence developmental gene expression by shaping locus topology. We refer to these elements as topological regulatory elements (TopoREs) to reflect their dual roles in genome organisation and gene expression. Importantly, these elements do not harbour autonomous transcriptional activation capacity, and instead appear to facilitate enhancer-promoter interactions, contributing to robust and precise timing of transcription. We discuss examples of TopoREs from two classes that are either dependent or independent of CTCF binding. Importantly, identification and interpretation of TopoRE function may shed light on multiple aspects of gene regulation, including the relationship between enhancer-promoter proximity and transcription, and enhancer-promoter specificity. Ultimately, understanding TopoRE diversity and function will aid in the interpretation of how human sequence variation can impact transcription and contribute to disease phenotypes.

Single-cell DNA methylome and 3D multi-omic atlas of the adult mouse brain

Cytosine DNA methylation is essential in brain development and is implicated in various neurological disorders. Understanding DNA methylation diversity across the entire brain in a spatial context is fundamental for a complete molecular atlas of brain cell types and their gene regulatory landscapes. Here we used single-nucleus methylome sequencing (snmC-seq3) and multi-omic sequencing (snm3C-seq)1 technologies to generate 301,626 methylomes and 176,003 chromatin conformation-methylome joint profiles from 117 dissected regions throughout the adult mouse brain. Using iterative clustering and integrating with companion whole-brain transcriptome and chromatin accessibility datasets, we constructed a methylation-based cell taxonomy with 4,673 cell groups and 274 cross-modality-annotated subclasses. We identified 2.6 million differentially methylated regions across the genome that represent potential gene regulation elements. Notably, we observed spatial cytosine methylation patterns on both genes and regulatory elements in cell types within and across brain regions. Brain-wide spatial transcriptomics data validated the association of spatial epigenetic diversity with transcription and improved the anatomical mapping of our epigenetic datasets. Furthermore, chromatin conformation diversities occurred in important neuronal genes and were highly associated with DNA methylation and transcription changes. Brain-wide cell-type comparisons enabled the construction of regulatory networks that incorporate transcription factors, regulatory elements and their potential downstream gene targets. Finally, intragenic DNA methylation and chromatin conformation patterns predicted alternative gene isoform expression observed in a whole-brain SMART-seq2 dataset. Our study establishes a brain-wide, single-cell DNA methylome and 3D multi-omic atlas and provides a valuable resource for comprehending the cellular-spatial and regulatory genome diversity of the mouse brain.

Expanding the coverage of regulons from high-confidence prior knowledge for accurate estimation of transcription factor activities

Gene regulation plays a critical role in the cellular processes that underlie human health and disease. The regulatory relationship between transcription factors (TFs), key regulators of gene expression, and their target genes, the so called TF regulons, can be coupled with computational algorithms to estimate the activity of TFs. However, to interpret these findings accurately, regulons of high reliability and coverage are needed. In this study, we present and evaluate a collection of regulons created using the CollecTRI meta-resource containing signed TF-gene interactions for 1186 TFs. In this context, we introduce a workflow to integrate information from multiple resources and assign the sign of regulation to TF-gene interactions that could be applied to other comprehensive knowledge bases. We find that the signed CollecTRI-derived regulons outperform other public collections of regulatory interactions in accurately inferring changes in TF activities in perturbation experiments. Furthermore, we showcase the value of the regulons by examining TF activity profiles in three different cancer types and exploring TF activities at the level of single-cells. Overall, the CollecTRI-derived TF regulons enable the accurate and comprehensive estimation of TF activities and thereby help to interpret transcriptomics data.

Epstein-Barr virus and host cell 3D genome organization

The human genome is organized in an extremely complexed yet ordered way within the nucleus. Genome organization plays a critical role in the regulation of gene expression. Viruses manipulate the host machinery to influence host genome organization to favor their survival and promote disease development. Epstein-Barr virus (EBV) is a common human virus, whose infection is associated with various diseases, including infectious mononucleosis, cancer, and autoimmune disorders. This review summarizes our current knowledge of how EBV uses different strategies to control the cellular 3D genome organization to affect cell gene expression to transform normal cells into lymphoblasts.

Gene regulatory network inference in the era of single-cell multi-omics

The interplay between chromatin, transcription factors and genes generates complex regulatory circuits that can be represented as gene regulatory networks (GRNs). The study of GRNs is useful to understand how cellular identity is established, maintained and disrupted in disease. GRNs can be inferred from experimental data - historically, bulk omics data - and/or from the literature. The advent of single-cell multi-omics technologies has led to the development of novel computational methods that leverage genomic, transcriptomic and chromatin accessibility information to infer GRNs at an unprecedented resolution. Here, we review the key principles of inferring GRNs that encompass transcription factor-gene interactions from transcriptomics and chromatin accessibility data. We focus on the comparison and classification of methods that use single-cell multimodal data. We highlight challenges in GRN inference, in particular with respect to benchmarking, and potential further developments using additional data modalities.

Idiosyncratic and dose-dependent epistasis drives variation in tomato fruit size

Epistasis between genes is traditionally studied with mutations that eliminate protein activity, but most natural genetic variation is in cis-regulatory DNA and influences gene expression and function quantitatively. In this study, we used natural and engineered cis-regulatory alleles in a plant stem-cell circuit to systematically evaluate epistatic relationships controlling tomato fruit size. Combining a promoter allelic series with two other loci, we collected over 30,000 phenotypic data points from 46 genotypes to quantify how allele strength transforms epistasis. We revealed a saturating dose-dependent relationship but also allele-specific idiosyncratic interactions, including between alleles driving a step change in fruit size during domestication. Our approach and findings expose an underexplored dimension of epistasis, in which cis-regulatory allelic diversity within gene regulatory networks elicits nonlinear, unpredictable interactions that shape phenotypes.

An enhancer RNA recruits MLL1 to regulate transcription of Myb

The Myb proto-oncogene encodes the transcription factor c-MYB, which is critical for hematopoiesis. Distant enhancers of Myb form a hub of interactions with the Myb promoter. We identified a long non-coding RNA (Myrlin) originating from the -81 kb murine Myb enhancer. Myrlin and Myb are coordinately regulated during erythroid differentiation. Myrlin TSS deletion using CRISPR/Cas9 reduced Myrlin and Myb expression and LDB1 complex occupancy at the Myb enhancers, compromising enhancer contacts and reducing RNA Pol II occupancy in the locus. In contrast, CRISPRi silencing of Myrlin left LDB1 and the Myb enhancer hub unperturbed, although Myrlin and Myb expression were downregulated, decoupling transcription and chromatin looping. Myrlin interacts with the MLL1 complex. Myrlin CRISPRi compromised MLL1 occupancy in the Myb locus, decreasing CDK9 and RNA Pol II binding and resulting in Pol II pausing in the Myb first exon/intron. Thus, Myrlin directly participates in activating Myb transcription by recruiting MLL1.

Yeast Heterochromatin Only Stably Silences Weak Regulatory Elements by Altering Burst Duration

The interplay between nucleosomes and transcription factors leads to programs of gene expression. Transcriptional silencing involves the generation of a chromatin state that represses transcription and is faithfully propagated through DNA replication and cell division. Using multiple reporter assays, including directly visualizing transcription in single cells, we investigated a diverse set of UAS enhancers and core promoters for their susceptibility to heterochromatic gene silencing. These results show that heterochromatin only stably silences weak and stress induced regulatory elements but is unable to stably repress housekeeping gene regulatory elements and the partial repression did not result in bistable expression states. Permutation analysis of different UAS enhancers and core promoters indicate that both elements function together to determine the susceptibility of regulatory sequences to repression. Specific histone modifiers and chromatin remodellers function in an enhancer specific manner to aid these elements to resist repression suggesting that Sir proteins likely function in part by reducing nucleosome mobility. We also show that the strong housekeeping regulatory elements can be repressed if silencer bound Sir1 is increased, suggesting that Sir1 is a limiting component in silencing. Together, our data suggest that the heterochromatic locus has been optimized to stably silence the weak mating type gene regulatory elements but not strong housekeeping gene regulatory sequences which could help explain why these genes are often found at the boundaries of silenced domains.

Runt-related Transcription Factors and Gene Regulatory Mechanisms in Skeletal Development and Diseases

PURPOSE OF REVIEW

Runt-related transcription factors (RUNX) play critical roles in skeletal development, metabolism, and diseases. In mammals, three RUNX members, namely RUNX1, RUNX2, and RUNX3, play distinct and redundant roles, although RUNX2 is a dominant factor in skeletal development and several skeletal diseases. This review is to provide an overview of the current understanding of RUNX-mediated transcriptional regulation in different skeletal cell types.

RECENT FINDINGS

Advances in chromatin immunoprecipitation and next-generation sequencing (ChIP-seq) have revealed genome-wide RUNX-mediated gene regulatory mechanisms, including their association with cis-regulatory elements and putative target genes. Further studies with genome-wide analysis and biochemical assays have shed light on RUNX-mediated pioneering action and involvements of RUNX2 in lipid-lipid phase separation. Emerging multi-layered mechanisms of RUNX-mediated gene regulations help us better understanding of skeletal development and diseases, which also provides clues to think how genome-wide studies can help develop therapeutic strategies for skeletal diseases.

May the force be with you: Nuclear condensates function beyond transcription control: Potential nongenetic functions of nuclear condensates in physiological and pathological conditions

Over the past decade, research has revealed biomolecular condensates' relevance in diverse cellular functions. Through a phase separation process, they concentrate macromolecules in subcompartments shaping the cellular organization and physiology. In the nucleus, biomolecular condensates assemble relevant biomolecules that orchestrate gene expression. We here hypothesize that chromatin condensates can also modulate the nongenetic functions of the genome, including the nuclear mechanical properties. The importance of chromatin condensates is supported by the genetic evidence indicating that mutations in their members are causative of a group of rare Mendelian diseases named chromatinopathies (CPs). Despite a broad spectrum of clinical features and the perturbations of the epigenetic machinery characterizing the CPs, recent findings highlighted negligible changes in gene expression. These data argue in favor of possible noncanonical functions of chromatin condensates in regulating the genome's spatial organization and, consequently, the nuclear mechanics. In this review, we discuss how condensates may impact nuclear mechanical properties, thus affecting the cellular response to mechanical cues and, eventually, cell fate and identity. Chromatin condensates organize macromolecules in the nucleus orchestrating the transcription regulation and mutations in their members are responsible for rare diseases named chromatinopathies. We argue that chromatin condensates, in concert with the nuclear lamina, may also govern the nuclear mechanical properties affecting the cellular response to external cues.

Frontiers and techniques in plant gene regulation

Understanding plant gene regulation has been a priority for generations of plant scientists. However, due to its complex nature, the regulatory code governing plant gene expression has yet to be deciphered comprehensively. Recently developed methods-often relying on next-generation sequencing technology and state-of-the-art computational approaches-have started to further our understanding of the gene regulatory logic used by plants. In this review, we discuss these methods and the insights into the regulatory code of plants that they can yield.

Metazoan nuclear pore complexes in gene regulation and genome stability

The nuclear pore complexes (NPCs), one of the hallmarks of eukaryotic nuclei, allow selective transport of macromolecules between the cytoplasm and the nucleus. Besides this canonical function, an increasing number of additional roles have been attributed to the NPCs and their constituents, the nucleoporins. Here we review recent insights into the mechanisms by which NPCs and nucleoporins affect transcription and DNA repair in metazoans. In the first part, we discuss how gene expression can be affected by the localization of genome-nucleoporin interactions at pores or "off-pores", by the role of nucleoporins in chromatin organization at different scales, or by the physical properties of nucleoporins. In the second part, we review the contribution of NPCs to genome stability, including transport-dependent and -independent functions and the role of positioning at NPCs in the repair of heterochromatic breaks and the regulation of replication stress.

Epromoters are new players in the regulatory landscape with potential pleiotropic roles

Precise spatiotemporal control of gene expression during normal development and cell differentiation is achieved by the combined action of proximal (promoters) and distal (enhancers) cis-regulatory elements. Recent studies have reported that a subset of promoters, termed Epromoters, works also as enhancers to regulate distal genes. This new paradigm opened novel questions regarding the complexity of our genome and raises the possibility that genetic variation within Epromoters has pleiotropic effects on various physiological and pathological traits by differentially impacting multiple proximal and distal genes. Here, we discuss the different observations pointing to an important role of Epromoters in the regulatory landscape and summarize the evidence supporting a pleiotropic impact of these elements in disease. We further hypothesize that Epromoter might represent a major contributor to phenotypic variation and disease.

Physical interaction between MSL2 and CLAMP assures direct cooperativity and prevents competition at composite binding sites

MSL2, the DNA-binding subunit of the Drosophila dosage compensation complex, cooperates with the ubiquitous protein CLAMP to bind MSL recognition elements (MREs) on the X chromosome. We explore the nature of the cooperative binding to these GA-rich, composite sequence elements in reconstituted naïve embryonic chromatin. We found that the cooperativity requires physical interaction between both proteins. Remarkably, disruption of this interaction does not lead to indirect, nucleosome-mediated cooperativity as expected, but to competition. The protein interaction apparently not only increases the affinity for composite binding sites, but also locks both proteins in a defined dimeric state that prevents competition. High Affinity Sites of MSL2 on the X chromosome contain variable numbers of MREs. We find that the cooperation between MSL2/CLAMP is not influenced by MRE clustering or arrangement, but happens largely at the level of individual MREs. The sites where MSL2/CLAMP bind strongly in vitro locate to all chromosomes and show little overlap to an expanded set of X-chromosomal MSL2 in vivo binding sites generated by CUT&RUN. Apparently, the intrinsic MSL2/CLAMP cooperativity is limited to a small selection of potential sites in vivo. This restriction must be due to components missing in our reconstitution, such as roX2 lncRNA.

Sedentary Behavior Impacts on the Epigenome and Transcriptome: Lessons from Muscle Inactivation in Drosophila Larvae

The biological mechanisms linking sedentary lifestyles and metabolic derangements are incompletely understood. In this study, temporal muscle inactivation in Drosophila larvae carrying a temperature-sensitive mutation in the shibire (shi1) gene was induced to mimic sedentary behavior during early life and study its transcriptional outcome. Our findings indicated a significant change in the epigenetic profile, as well as the genomic profile, of RNA Pol II binding in the inactive muscles relative to control, within a relatively short time period. Whole-genome analysis of RNA-Pol II binding to DNA by muscle-specific targeted DamID (TaDa) protocol revealed that muscle inactivity altered Pol II binding in 121 out of 2010 genes (6%), with a three-fold enrichment of genes coding for lncRNAs. The suppressed protein-coding genes included genes associated with longevity, DNA repair, muscle function, and ubiquitin-dependent proteostasis. Moreover, inducing muscle inactivation exerted a multi-level impact upon chromatin modifications, triggering an altered epigenetic balance of active versus inactive marks. The downregulated genes in the inactive muscles included genes essential for muscle structure and function, carbohydrate metabolism, longevity, and others. Given the multiple analogous genes in Drosophila for many human genes, extrapolating our findings to humans may hold promise for establishing a molecular link between sedentary behavior and metabolic diseases.

The architecture of binding cooperativity between densely bound transcription factors

The binding of transcription factors (TFs) along genomes is restricted to a subset of sites containing their preferred motifs. TF-binding specificity is often attributed to the co-binding of interacting TFs; however, apart from specific examples, this model remains untested. Here, we define dependencies among budding yeast TFs that localize to overlapping promoters by profiling the genome-wide consequences of co-depleting multiple TFs. We describe unidirectional interactions, revealing Msn2 as a central factor allowing TF binding at its target promoters. By contrast, no case of mutual cooperation was observed. Particularly, Msn2 retained binding at its preferred promoters upon co-depletion of fourteen similarly bound TFs. Overall, the consequences of TF co-depletions were moderate, limited to a subset of promoters, and failed to explain the role of regions outside the DNA-binding domain in directing TF-binding preferences. Our results call for re-evaluating the role of cooperative interactions in directing TF-binding preferences.

Active learning of enhancer and silencer regulatory grammar in photoreceptors

Cis-regulatory elements (CREs) direct gene expression in health and disease, and models that can accurately predict their activities from DNA sequences are crucial for biomedicine. Deep learning represents one emerging strategy to model the regulatory grammar that relates CRE sequence to function. However, these models require training data on a scale that exceeds the number of CREs in the genome. We address this problem using active machine learning to iteratively train models on multiple rounds of synthetic DNA sequences assayed in live mammalian retinas. During each round of training the model actively selects sequence perturbations to assay, thereby efficiently generating informative training data. We iteratively trained a model that predicts the activities of sequences containing binding motifs for the photoreceptor transcription factor Cone-rod homeobox (CRX) using an order of magnitude less training data than current approaches. The model's internal confidence estimates of its predictions are reliable guides for designing sequences with high activity. The model correctly identified critical sequence differences between active and inactive sequences with nearly identical transcription factor binding sites, and revealed order and spacing preferences for combinations of motifs. Our results establish active learning as an effective method to train accurate deep learning models of cis-regulatory function after exhausting naturally occurring training examples in the genome.

Dynamic 3D genome reorganization during senescence: defining cell states through chromatin

Cellular senescence, a cell state characterized by growth arrest and insensitivity to growth stimulatory hormones, is accompanied by a massive change in chromatin organization. Senescence can be induced by a range of physiological signals and pathological stresses and was originally thought to be an irreversible state, implicated in normal development, wound healing, tumor suppression and aging. Recently cellular senescence was shown to be reversible in some cases, with exit being triggered by the modulation of the cell's transcriptional program by the four Yamanaka factors, the suppression of p53 or H3K9me3, PDK1, and/or depletion of AP-1. Coincident with senescence reversal are changes in chromatin organization, most notably the loss of senescence-associated heterochromatin foci (SAHF) found in oncogene-induced senescence. In addition to fixed-cell imaging, chromatin conformation capture and multi-omics have been used to examine chromatin reorganization at different spatial resolutions during senescence. They identify determinants of SAHF formation and other key features that differentiate distinct types of senescence. Not surprisingly, multiple factors, including the time of induction, the type of stress experienced, and the type of cell involved, influence the global reorganization of chromatin in senescence. Here we discuss how changes in the three-dimensional organization of the genome contribute to the regulation of transcription at different stages of senescence. In particular, the distinct contributions of heterochromatin- and lamina-mediated interactions, changes in gene expression, and other cellular control mechanisms are discussed. We propose that high-resolution temporal and spatial analyses of the chromatin landscape during senescence will identify early markers of the different senescence states to help guide clinical diagnosis.

Shared Gene Targets of the ATF4 and p53 Transcriptional Networks

The master tumor suppressor p53 regulates multiple cell fate decisions, such as cell cycle arrest and apoptosis, via transcriptional control of a broad gene network. Dysfunction in the p53 network is common in cancer, often through mutations that inactivate p53 or other members of the pathway. Induction of tumor-specific cell death by restoration of p53 activity without off-target effects has gained significant interest in the field. In this study, we explore the gene regulatory mechanisms underlying a putative anticancer strategy involving stimulation of the p53-independent integrated stress response (ISR). Our data demonstrate the p53 and ISR pathways converge to independently regulate common metabolic and proapoptotic genes. We investigated the architecture of multiple gene regulatory elements bound by p53 and the ISR effector ATF4 controlling this shared regulation. We identified additional key transcription factors that control basal and stress-induced regulation of these shared p53 and ATF4 target genes. Thus, our results provide significant new molecular and genetic insight into gene regulatory networks and transcription factors that are the target of numerous antitumor therapies.

A variegated model of transcription factor function in the immune system

Specific combinations of transcription factors (TFs) control the gene expression programs that underlie specialized immune responses. Previous models of TF function in immunocytes had restricted each TF to a single functional categorization [e.g., lineage-defining (LDTFs) vs. signal-dependent TFs (SDTFs)] within one cell type. Synthesizing recent results, we instead propose a variegated model of immunological TF function, whereby many TFs have flexible and different roles across distinct cell states, contributing to cell phenotypic diversity. We discuss evidence in support of this variegated model, describe contextual inputs that enable TF diversification, and look to the future to imagine warranted experimental and computational tools to build quantitative and predictive models of immunocyte gene regulatory networks.

Shared gene targets of the ATF4 and p53 transcriptional networks

The master tumor suppressor p53 regulates multiple cell fate decisions, like cell cycle arrest and apoptosis, via transcriptional control of a broad gene network. Dysfunction in the p53 network is common in cancer, often through mutations that inactivate p53 or other members of the pathway. Induction of tumor-specific cell death by restoration of p53 activity without off-target effects has gained significant interest in the field. In this study, we explore the gene regulatory mechanisms underlying a putative anti-cancer strategy involving stimulation of the p53-independent Integrated Stress Response (ISR). Our data demonstrate the p53 and ISR pathways converge to independently regulate common metabolic and pro-apoptotic genes. We investigated the architecture of multiple gene regulatory elements bound by p53 and the ISR effector ATF4 controlling this shared regulation. We identified additional key transcription factors that control basal and stress-induced regulation of these shared p53 and ATF4 target genes. Thus, our results provide significant new molecular and genetic insight into gene regulatory networks and transcription factors that are the target of numerous antitumor therapies.

Integrative approaches to study enhancer-promoter communication

The spatiotemporal control of gene expression in complex multicellular organisms relies on noncoding regulatory sequences such as enhancers, which activate transcription of target genes often over large genomic distances. Despite the advances in the identification and characterization of enhancers, the principles and mechanisms by which enhancers select and control their target genes remain largely unknown. Here, we review recent interdisciplinary and quantitative approaches based on emerging techniques that aim to address open questions in the field, notably how regulatory information is encoded in the DNA sequence, how this information is transferred from enhancers to promoters, and how these processes are regulated in time.

DNA-guided transcription factor cooperativity shapes face and limb mesenchyme

Transcription factors (TFs) can define distinct cellular identities despite nearly identical DNA-binding specificities. One mechanism for achieving regulatory specificity is DNA-guided TF cooperativity. Although in vitro studies suggest it may be common, examples of such cooperativity remain scarce in cellular contexts. Here, we demonstrate how 'Coordinator', a long DNA motif comprised of common motifs bound by many basic helix-loop-helix (bHLH) and homeodomain (HD) TFs, uniquely defines regulatory regions of embryonic face and limb mesenchyme. Coordinator guides cooperative and selective binding between the bHLH family mesenchymal regulator TWIST1 and a collective of HD factors associated with regional identities in the face and limb. TWIST1 is required for HD binding and open chromatin at Coordinator sites, while HD factors stabilize TWIST1 occupancy at Coordinator and titrate it away from HD-independent sites. This cooperativity results in shared regulation of genes involved in cell-type and positional identities, and ultimately shapes facial morphology and evolution.

Multiomics Analysis Reveals Novel Genetic Determinants for Lens Differentiation, Structure, and Transparency

Recent advances in next-generation sequencing and data analysis have provided new gateways for identification of novel genome-wide genetic determinants governing tissue development and disease. These advances have revolutionized our understanding of cellular differentiation, homeostasis, and specialized function in multiple tissues. Bioinformatic and functional analysis of these genetic determinants and the pathways they regulate have provided a novel basis for the design of functional experiments to answer a wide range of long-sought biological questions. A well-characterized model for the application of these emerging technologies is the development and differentiation of the ocular lens and how individual pathways regulate lens morphogenesis, gene expression, transparency, and refraction. Recent applications of next-generation sequencing analysis on well-characterized chicken and mouse lens differentiation models using a variety of omics techniques including RNA-seq, ATAC-seq, whole-genome bisulfite sequencing (WGBS), chip-seq, and CUT&RUN have revealed a wide range of essential biological pathways and chromatin features governing lens structure and function. Multiomics integration of these data has established new gene functions and cellular processes essential for lens formation, homeostasis, and transparency including the identification of novel transcription control pathways, autophagy remodeling pathways, and signal transduction pathways, among others. This review summarizes recent omics technologies applied to the lens, methods for integrating multiomics data, and how these recent technologies have advanced our understanding ocular biology and function. The approach and analysis are relevant to identifying the features and functional requirements of more complex tissues and disease states.