Enhancer-gene specificity in development and disease

Genetic factors mediating long-range enhancer-promoter communication in mammalian development

Enhancers are remotely located noncoding DNA sequences that regulate gene expression in response to developmental, homeostatic, and environmental cues. Canonical short-range enhancers located <50 kb from their cognate promoters function by binding transcription factors, coactivators, and chromatin modifiers. In this review, we discuss recent evidence that medium-range (50-400 kb) and long-range (>400 kb) enhancers rely on additional mechanisms, including cohesin, CCCTC-binding factor, and high-affinity protein-protein interactions. These mechanisms are crucial for establishing the physical proximity and interaction between enhancers and their target promoters over extended genomic distances and ensuring robust gene activation during mammalian development. Future studies will be critical to unravel their prevalence and evolutionary significance across various genomic loci, cell types, and species.

Enhancer cooperativity can compensate for loss of activity over large genomic distances

Enhancers are short DNA sequences that activate their target promoter from a distance; however, increasing the genomic distance between the enhancer and the promoter decreases expression levels. Many genes are controlled by combinations of multiple enhancers, yet the interaction and cooperation of individual enhancer elements are not well understood. Here, we developed a synthetic platform in mouse embryonic stem cells that allows building complex regulatory landscapes from the bottom up. We tested the system by integrating individual enhancers at different distances and confirmed that the strength of an enhancer contributes to how strongly it is affected by increased genomic distance. Furthermore, synergy between two enhancer elements depends on the distance at which the two elements are integrated: introducing a weak enhancer between a strong enhancer and the promoter strongly increases reporter gene expression, allowing enhancers to activate from increased genomic distances.

Rapid and quantitative functional interrogation of human enhancer variant activity in live mice

Functional analysis of non-coding variants associated with congenital disorders remains challenging due to the lack of efficient in vivo models. Here we introduce dual-enSERT, a robust Cas9-based two-color fluorescent reporter system which enables rapid, quantitative comparison of enhancer allele activities in live mice in less than two weeks. We use this technology to examine and measure the gain- and loss-of-function effects of enhancer variants previously linked to limb polydactyly, autism spectrum disorder, and craniofacial malformation. By combining dual-enSERT with single-cell transcriptomics, we characterise gene expression in cells where the enhancer is normally and ectopically active, revealing candidate pathways that may lead to enhancer misregulation. Finally, we demonstrate the widespread utility of dual-enSERT by testing the effects of fifteen previously uncharacterised rare and common non-coding variants linked to neurodevelopmental disorders. In doing so we identify variants that reproducibly alter the in vivo activity of OTX2 and MIR9-2 brain enhancers, implicating them in autism. Dual-enSERT thus allows researchers to go from identifying candidate enhancer variants to analysis of comparative enhancer activity in live embryos in under two weeks.

WT1 and DNMT3A Mutations in Prognostic Significance of Acute Myeloid Leukemia: A Meta-Analysis

Background: Adult acute leukemia most commonly manifests as acute myeloid leukemia (AML), a highly heterogeneous malignant tumor of the blood system. The application of genetic diagnostic technology is currently prevalent in numerous clinical sectors. According to recent research, the presence of specific gene mutations or rearrangements in leukemia cells is the primary cause of the disease. As different types of leukemia are caused by atypical mutated genes, testing for these mutations or rearrangements can help diagnose leukemia and identify the disease's molecular targets for treatment. Methods: Using the search fields "WT1," "DNMT3A," "Acute myeloid leukemia," and "survival," the CBM, Cochrane Library, Scopus, EMBASE, and PUBMED databases were separately reviewed. The methodology for evaluating the risk of bias developed by the Cochrane Collaboration was used in conjunction with a methodical evaluation of pertinent literature. Excluded studies with the following characteristics: (1) incomplete and repetitive publications, (2) unable to retrieve or convert data, (3) non-English or Chinese articles. Results: This analysis included 13 studies covering a total of 3478 subjects. The frequency of Wilms' Tumor 1 (WT1) mutations is 6.7%-35.73%, and the frequency of DNMT3A mutations is 12.06%-51.1%. The remission rate of patients with WT1 mutations was less than that of patients without WT1 mutations (OR = 0.22; 95% confidence interval [CI]: 0.14, 0.36; p < 0.00001; I2 = 55%). The DNMT3A mutation has no statistical significance for the prognosis of AML (OR = 1.21; 95% CI: 0.93, 1.58; p = 0.16; I2 = 80%). After removing one study, the heterogeneity of the indicator (mitigation rate) among other studies of DNMT3A mutation was dramatically reduced (OR = 0.63; 95% CI: 0.43, 0.93; p = 0.02; I2 = 0%). Conclusions: Our meta-analysis shows that WT1 mutations hurt the remission rate of AML. Moreover, the impact of DNMT3A mutations on AML needs to be treated with caution. Gene diagnosis is critical for the prognosis and clinical management of AML.

Active enhancers: recent research advances and insights into disease

Precise regulation of gene expression is crucial to development. Enhancers, the core of gene regulation, determine the spatiotemporal pattern of gene transcription. Since many disease-associated mutations are characterized in enhancers, the research on enhancer will provide clues to precise medicine. Rapid advances in high-throughput sequencing technology facilitate the characterization of enhancers at genome wide, but understanding the functional mechanisms of enhancers remains challenging. Herein, we provide a panorama of enhancer characteristics, including epigenetic modifications, enhancer transcripts, and enhancer-promoter interaction patterns. Furthermore, we outline the applications of high-throughput sequencing technology and functional genomics methods in enhancer research. Finally, we discuss the role of enhancers in human disease and their potential as targets for disease prevention and treatment strategies.

A gene desert required for regulatory control of pleiotropic Shox2 expression and embryonic survival

Approximately a quarter of the human genome consists of gene deserts, large regions devoid of genes often located adjacent to developmental genes and thought to contribute to their regulation. However, defining the regulatory functions embedded within these deserts is challenging due to their large size. Here, we explore the cis-regulatory architecture of a gene desert flanking the Shox2 gene, which encodes a transcription factor indispensable for proximal limb, craniofacial, and cardiac pacemaker development. We identify the gene desert as a regulatory hub containing more than 15 distinct enhancers recapitulating anatomical subdomains of Shox2 expression. Ablation of the gene desert leads to embryonic lethality due to Shox2 depletion in the cardiac sinus venosus, caused in part by the loss of a specific distal enhancer. The gene desert is also required for stylopod morphogenesis, mediated via distributed proximal limb enhancers. In summary, our study establishes a multi-layered role of the Shox2 gene desert in orchestrating pleiotropic developmental expression through modular arrangement and coordinated dynamics of tissue-specific enhancers.

A Ctnnb1 enhancer transcriptionally regulates Wnt signaling dosage to balance homeostasis and tumorigenesis of intestinal epithelia

The β-catenin-dependent canonical Wnt signaling is pivotal in organ development, tissue homeostasis, and cancer. Here, we identified an upstream enhancer of Ctnnb1 - the coding gene for β-catenin, named ieCtnnb1 (intestinal enhancer of Ctnnb1), which is crucial for intestinal homeostasis. ieCtnnb1 is predominantly active in the base of small intestinal crypts and throughout the epithelia of large intestine. Knockout of ieCtnnb1 led to a reduction in Ctnnb1 transcription, compromising the canonical Wnt signaling in intestinal crypts. Single-cell sequencing revealed that ieCtnnb1 knockout altered epithelial compositions and potentially compromised functions of small intestinal crypts. While deletion of ieCtnnb1 hampered epithelial turnovers in physiologic conditions, it prevented occurrence and progression of Wnt/β-catenin-driven colorectal cancers. Human ieCTNNB1 drove reporter gene expression in a pattern highly similar to mouse ieCtnnb1. ieCTNNB1 contains a single-nucleotide polymorphism associated with CTNNB1 expression levels in human gastrointestinal epithelia. The enhancer activity of ieCTNNB1 in colorectal cancer tissues was stronger than that in adjacent normal tissues. HNF4α and phosphorylated CREB1 were identified as key trans-factors binding to ieCTNNB1 and regulating CTNNB1 transcription. Together, these findings unveil an enhancer-dependent mechanism controlling the dosage of Wnt signaling and homeostasis in intestinal epithelia.

Predictive Prioritization of Enhancers Associated with Pancreas Disease Risk

Genetic and epigenetic variations in regulatory enhancer elements increase susceptibility to a range of pathologies. Despite recent advances, linking enhancer elements to target genes and predicting transcriptional outcomes of enhancer dysfunction remain significant challenges. Using 3D chromatin conformation assays, we generated an extensive enhancer interaction dataset for the human pancreas, encompassing more than 20 donors and five major cell types, including both exocrine and endocrine compartments. We employed a network approach to parse chromatin interactions into enhancer-promoter tree models, facilitating a quantitative, genome-wide analysis of enhancer connectivity. With these tree models, we developed a machine learning algorithm to estimate the impact of enhancer perturbations on cell type-specific gene expression in the human pancreas. Orthogonal to our computational approach, we perturbed enhancer function in primary human pancreas cells using CRISPR interference and quantified the effects at the single-cell level through RNA FISH coupled with high-throughput imaging. Our enhancer tree models enabled the annotation of common germline risk variants associated with pancreas diseases, linking them to putative target genes in specific cell types. For pancreatic ductal adenocarcinoma, we found a stronger enrichment of disease susceptibility variants within acinar cell regulatory elements, despite ductal cells historically being assumed as the primary cell-of-origin. Our integrative approach-combining cell type-specific enhancer-promoter interaction mapping, computational models, and single-cell enhancer perturbation assays-produced a robust resource for studying the genetic basis of pancreas disorders.

Integration of chromosome locations and functional aspects of enhancers and topologically associating domains in knowledge graphs enables versatile queries about gene regulation

Knowledge about transcription factor binding and regulation, target genes, cis-regulatory modules and topologically associating domains is not only defined by functional associations like biological processes or diseases but also has a determinative genome location aspect. Here, we exploit these location and functional aspects together to develop new strategies to enable advanced data querying. Many databases have been developed to provide information about enhancers, but a schema that allows the standardized representation of data, securing interoperability between resources, has been lacking. In this work, we use knowledge graphs for the standardized representation of enhancers and topologically associating domains, together with data about their target genes, transcription factors, location on the human genome, and functional data about diseases and gene ontology annotations. We used this schema to integrate twenty-five enhancer datasets and two domain datasets, creating the most powerful integrative resource in this field to date. The knowledge graphs have been implemented using the Resource Description Framework and integrated within the open-access BioGateway knowledge network, generating a resource that contains an interoperable set of knowledge graphs (enhancers, TADs, genes, proteins, diseases, GO terms, and interactions between domains). We show how advanced queries, which combine functional and location restrictions, can be used to develop new hypotheses about functional aspects of gene expression regulation.

Cooperative insulation of regulatory domains by CTCF-dependent physical insulation and promoter competition

The specificity of gene expression during development requires the insulation of regulatory domains to avoid inappropriate enhancer-gene interactions. In vertebrates, this insulator function is mostly attributed to clusters of CTCF sites located at topologically associating domain (TAD) boundaries. However, TAD boundaries allow some physical crosstalk across regulatory domains, which is at odds with the specific and precise expression of developmental genes. Here we show that developmental genes and nearby clusters of CTCF sites cooperatively foster the robust insulation of regulatory domains. By genetically dissecting a couple of representative loci in mouse embryonic stem cells, we show that CTCF sites prevent undesirable enhancer-gene contacts (i.e. physical insulation), while developmental genes preferentially contribute to regulatory insulation through non-structural mechanisms involving promoter competition rather than enhancer blocking. Overall, our work provides important insights into the insulation of regulatory domains, which in turn might help interpreting the pathological consequences of certain structural variants.

High-throughput affinity measurements of direct interactions between activation domains and co-activators

Sequence-specific activation by transcription factors is essential for gene regulation1,2. Key to this are activation domains, which often fall within disordered regions of transcription factors3,4 and recruit co-activators to initiate transcription5. These interactions are difficult to characterize via most experimental techniques because they are typically weak and transient6,7. Consequently, we know very little about whether these interactions are promiscuous or specific, the mechanisms of binding, and how these interactions tune the strength of gene activation. To address these questions, we developed a microfluidic platform for expression and purification of hundreds of activation domains in parallel followed by direct measurement of co-activator binding affinities (STAMMPPING, for Simultaneous Trapping of Affinity Measurements via a Microfluidic Protein-Protein INteraction Generator). By applying STAMMPPING to quantify direct interactions between eight co-activators and 204 human activation domains (>1,500 K ds), we provide the first quantitative map of these interactions and reveal 334 novel binding pairs. We find that the metazoan-specific co-activator P300 directly binds >100 activation domains, potentially explaining its widespread recruitment across the genome to influence transcriptional activation. Despite sharing similar molecular properties (e.g. enrichment of negative and hydrophobic residues), activation domains utilize distinct biophysical properties to recruit certain co-activator domains. Co-activator domain affinity and occupancy are well-predicted by analytical models that account for multivalency, and in vitro affinities quantitatively predict activation in cells with an ultrasensitive response. Not only do our results demonstrate the ability to measure affinities between even weak protein-protein interactions in high throughput, but they also provide a necessary resource of over 1,500 activation domain/co-activator affinities which lays the foundation for understanding the molecular basis of transcriptional activation.

Spatial orchestration of the genome: topological reorganisation during X-chromosome inactivation

Genomes are organised through hierarchical structures, ranging from local kilobase-scale cis-regulatory contacts to large chromosome territories. Most notably, (sub)-compartments partition chromosomes according to transcriptional activity, while topologically associating domains (TADs) define cis-regulatory landscapes. The inactive X chromosome in mammals has provided unique insights into the regulation and function of the three-dimensional (3D) genome. Concurrent with silencing of the majority of genes and major alterations of its chromatin state, the X chromosome undergoes profound spatial rearrangements at multiple scales. These include the emergence of megadomains, alterations of the compartment structure and loss of the majority of TADs. Moreover, the Xist locus, which orchestrates X-chromosome inactivation, has provided key insights into regulation and function of regulatory domains. This review provides an overview of recent insights into the control of these structural rearrangements and contextualises them within a broader understanding of 3D genome organisation.

Dual-role transcription factors stabilize intermediate expression levels

Precise control of gene expression levels is essential for normal cell functions, yet how they are defined and tightly maintained, particularly at intermediate levels, remains elusive. Here, using a series of newly developed sequencing, imaging, and functional assays, we uncover a class of transcription factors with dual roles as activators and repressors, referred to as condensate-forming level-regulating dual-action transcription factors (TFs). They reduce high expression but increase low expression to achieve stable intermediate levels. Dual-action TFs directly exert activating and repressing functions via condensate-forming domains that compartmentalize core transcriptional unit selectively. Clinically relevant mutations in these domains, which are linked to a range of developmental disorders, impair condensate selectivity and dual-action TF activity. These results collectively address a fundamental question in expression regulation and demonstrate the potential of level-regulating dual-action TFs as powerful effectors for engineering controlled expression levels.

Comparative 3D genome analysis between neural retina and retinal pigment epithelium reveals differential cis-regulatory interactions at retinal disease loci

BACKGROUND

Vision depends on the interplay between photoreceptor cells of the neural retina and the underlying retinal pigment epithelium (RPE). Most genes involved in inherited retinal diseases display specific spatiotemporal expression within these interconnected retinal components through the local recruitment of cis-regulatory elements (CREs) in 3D nuclear space.

RESULTS

To understand the role of differential chromatin architecture in establishing tissue-specific expression at inherited retinal disease loci, we mapped genome-wide chromatin interactions using in situ Hi-C and H3K4me3 HiChIP on neural retina and RPE/choroid from human adult donor eyes. We observed chromatin looping between active promoters and 32,425 and 8060 candidate CREs in the neural retina and RPE/choroid, respectively. A comparative 3D genome analysis between these two retinal tissues revealed that 56% of 290 known inherited retinal disease genes were marked by differential chromatin interactions. One of these was ABCA4, which is implicated in the most common autosomal recessive inherited retinal disease. We zoomed in on retina- and RPE-specific cis-regulatory interactions at the ABCA4 locus using high-resolution UMI-4C. Integration with bulk and single-cell epigenomic datasets and in vivo enhancer assays in zebrafish revealed tissue-specific CREs interacting with ABCA4.

CONCLUSIONS

Through comparative 3D genome mapping, based on genome-wide, promoter-centric, and locus-specific assays of human neural retina and RPE, we have shown that gene regulation at key inherited retinal disease loci is likely mediated by tissue-specific chromatin interactions. These findings do not only provide insight into tissue-specific regulatory landscapes at retinal disease loci, but also delineate the search space for non-coding genomic variation underlying unsolved inherited retinal diseases.

Genomic context sensitizes regulatory elements to genetic disruption

Genomic context critically modulates regulatory function but is difficult to manipulate systematically. The murine insulin-like growth factor 2 (Igf2)/H19 locus is a paradigmatic model of enhancer selectivity, whereby CTCF occupancy at an imprinting control region directs downstream enhancers to activate either H19 or Igf2. We used synthetic regulatory genomics to repeatedly replace the native locus with 157-kb payloads, and we systematically dissected its architecture. Enhancer deletion and ectopic delivery revealed previously uncharacterized long-range regulatory dependencies at the native locus. Exchanging the H19 enhancer cluster with the Sox2 locus control region (LCR) showed that the H19 enhancers relied on their native surroundings while the Sox2 LCR functioned autonomously. Analysis of regulatory DNA actuation across cell types revealed that these enhancer clusters typify broader classes of context sensitivity genome wide. These results show that unexpected dependencies influence even well-studied loci, and our approach permits large-scale manipulation of complete loci to investigate the relationship between regulatory architecture and function.

Genomic context sensitizes regulatory elements to genetic disruption

Enhancer function is frequently investigated piecemeal using truncated reporter assays or single deletion analysis. Thus it remains unclear to what extent enhancer function at native loci relies on surrounding genomic context. Using the Big-IN technology for targeted integration of large DNAs, we analyzed the regulatory architecture of the murine Igf2/H19 locus, a paradigmatic model of enhancer selectivity. We assembled payloads containing a 157-kb functional Igf2/H19 locus and engineered mutations to genetically direct CTCF occupancy at the imprinting control region (ICR) that switches the target gene of the H19 enhancer cluster. Contrasting activity of payloads delivered at the endogenous Igf2/H19 locus or ectopically at Hprt revealed that the Igf2/H19 locus includes additional, previously unknown long-range regulatory elements. Exchanging components of the Igf2/H19 locus with the well-studied Sox2 locus showed that the H19 enhancer cluster functioned poorly out of context, and required its native surroundings to activate Sox2 expression. Conversely, the Sox2 locus control region (LCR) could activate both Igf2 and H19 outside its native context, but its activity was only partially modulated by CTCF occupancy at the ICR. Analysis of regulatory DNA actuation across different cell types revealed that, while the H19 enhancers are tightly coordinated within their native locus, the Sox2 LCR acts more independently. We show that these enhancer clusters typify broader classes of loci genome-wide. Our results show that unexpected dependencies may influence even the most studied functional elements, and our synthetic regulatory genomics approach permits large-scale manipulation of complete loci to investigate the relationship between locus architecture and function.

Enhancer-promoter interactions can form independently of genomic distance and be functional across TAD boundaries

Topologically Associating Domains (TADs) have been suggested to facilitate and constrain enhancer-promoter interactions. However, the role of TAD boundaries in effectively restricting these interactions remains unclear. Here, we show that a significant proportion of enhancer-promoter interactions are established across TAD boundaries in Drosophila embryos, but that developmental genes are strikingly enriched in intra- but not inter-TAD interactions. We pursued this observation using the twist locus, a master regulator of mesoderm development, and systematically relocated one of its enhancers to various genomic locations. While this developmental gene can establish inter-TAD interactions with its enhancer, the functionality of these interactions remains limited, highlighting the existence of topological constraints. Furthermore, contrary to intra-TAD interactions, the formation of inter-TAD enhancer-promoter interactions is not solely driven by genomic distance, with distal interactions sometimes favored over proximal ones. These observations suggest that other general mechanisms must exist to establish and maintain specific enhancer-promoter interactions across large distances.

Multi-omics analysis in human retina uncovers ultraconserved cis-regulatory elements at rare eye disease loci

Cross-species genome comparisons have revealed a substantial number of ultraconserved non-coding elements (UCNEs). Several of these elements have proved to be essential tissue- and cell type-specific cis-regulators of developmental gene expression. Here, we characterize a set of UCNEs as candidate CREs (cCREs) during retinal development and evaluate the contribution of their genomic variation to rare eye diseases, for which pathogenic non-coding variants are emerging. Integration of bulk and single-cell retinal multi-omics data reveals 594 genes under potential cis-regulatory control of UCNEs, of which 45 are implicated in rare eye disease. Mining of candidate cis-regulatory UCNEs in WGS data derived from the rare eye disease cohort of Genomics England reveals 178 ultrarare variants within 84 UCNEs associated with 29 disease genes. Overall, we provide a comprehensive annotation of ultraconserved non-coding regions acting as cCREs during retinal development which can be targets of non-coding variation underlying rare eye diseases.

De novo enhancer deletion of LMX1B produces a mild nail-patella clinical phenotype

Critical genes involved in embryonic development are often transcription factors, regulating many downstream genes. LMX1B is a homeobox gene that is involved in formation of the limbs, eyes and kidneys, heterozygous loss-of-function sequence variants and deletions cause Nail-Patella syndrome. Most of the reported variants are localised within the gene's coding sequence, however, approximately 5%-10% of affected individuals do not have a pathogenic variant identified within this region. In this study, we present a family with four affected individuals across two generations with a deletion spanning a conserved upstream LMX1B-binding sequence. This deletion is de novo in the mother of three affected children. Furthermore, in this family, the manifestations appear limited to the nails and limbs, and therefore may reflect an attenuated phenotype of the classic Nail-Patella phenotype that includes ophthalmological and renal manifestations.

Rapid and Quantitative Functional Interrogation of Human Enhancer Variant Activity in Live Mice

Functional analysis of non-coding variants associated with human congenital disorders remains challenging due to the lack of efficient in vivo models. Here we introduce dual-enSERT, a robust Cas9-based two-color fluorescent reporter system which enables rapid, quantitative comparison of enhancer allele activities in live mice of any genetic background. We use this new technology to examine and measure the gain- and loss-of-function effects of enhancer variants linked to limb polydactyly, autism, and craniofacial malformation. By combining dual-enSERT with single-cell transcriptomics, we characterize variant enhancer alleles at cellular resolution, thereby implicating candidate molecular pathways in pathogenic enhancer misregulation. We further show that independent, polydactyly-linked enhancer variants lead to ectopic expression in the same cell populations, indicating shared genetic mechanisms underlying non-coding variant pathogenesis. Finally, we streamline dual-enSERT for analysis in F0 animals by placing both reporters on the same transgene separated by a synthetic insulator. Dual-enSERT allows researchers to go from identifying candidate enhancer variants to analysis of comparative enhancer activity in live embryos in under two weeks.

Transcriptional condensates and phase separation: condensing information across scales and mechanisms

Transcription is the fundamental process of gene expression, which in eukaryotes occurs within the complex physicochemical environment of the nucleus. Decades of research have provided extreme detail in the molecular and functional mechanisms of transcription, but the spatial and genomic organization of transcription remains mysterious. Recent discoveries show that transcriptional components can undergo phase separation and create distinct compartments inside the nucleus, providing new models through which to view the transcription process in eukaryotes. In this review, we focus on transcriptional condensates and their phase separation-like behaviors. We suggest differentiation between physical descriptions of phase separation and the complex and dynamic biomolecular assemblies required for productive gene expression, and we discuss how transcriptional condensates are central to organizing the three-dimensional genome across spatial and temporal scales. Finally, we map approaches for therapeutic manipulation of transcriptional condensates and ask what technical advances are needed to understand transcriptional condensates more completely.

What is an enhancer?

Tight control of the transcription process is essential for the correct spatial and temporal gene expression pattern during development and in homeostasis. Enhancers are at the core of correct transcriptional activation. The original definition of an enhancer is straightforward: a DNA sequence that activates transcription independent of orientation and direction. Dissection of numerous enhancer loci has shown that many enhancer-like elements might not conform to the original definition, suggesting that enhancers and enhancer-like elements might use multiple different mechanisms to contribute to transcriptional activation. Here, we review methodologies to identify enhancers and enhancer-like elements and discuss pitfalls and consequences for our understanding of transcriptional regulation.

Shaping faces: genetic and epigenetic control of craniofacial morphogenesis

Major differences in facial morphology distinguish vertebrate species. Variation of facial traits underlies the uniqueness of human individuals, and abnormal craniofacial morphogenesis during development leads to birth defects that significantly affect quality of life. Studies during the past 40 years have advanced our understanding of the molecular mechanisms that establish facial form during development, highlighting the crucial roles in this process of a multipotent cell type known as the cranial neural crest cell. In this Review, we discuss recent advances in multi-omics and single-cell technologies that enable genes, transcriptional regulatory networks and epigenetic landscapes to be closely linked to the establishment of facial patterning and its variation, with an emphasis on normal and abnormal craniofacial morphogenesis. Advancing our knowledge of these processes will support important developments in tissue engineering, as well as the repair and reconstruction of the abnormal craniofacial complex.

Transcription factor expression is the main determinant of variability in gene co-activity

Many genes are co-expressed and form genomic domains of coordinated gene activity. However, the regulatory determinants of domain co-activity remain unclear. Here, we leverage human individual variation in gene expression to characterize the co-regulatory processes underlying domain co-activity and systematically quantify their effect sizes. We employ transcriptional decomposition to extract from RNA expression data an expression component related to co-activity revealed by genomic positioning. This strategy reveals close to 1,500 co-activity domains, covering most expressed genes, of which the large majority are invariable across individuals. Focusing specifically on domains with high variability in co-activity reveals that contained genes have a higher sharing of eQTLs, a higher variability in enhancer interactions, and an enrichment of binding by variably expressed transcription factors, compared to genes within non-variable domains. Through careful quantification of the relative contributions of regulatory processes underlying co-activity, we find transcription factor expression levels to be the main determinant of gene co-activity. Our results indicate that distal trans effects contribute more than local genetic variation to individual variation in co-activity domains.

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.

POSTRE: a tool to predict the pathological effects of human structural variants

Understanding the pathological impact of non-coding genetic variation is a major challenge in medical genetics. Accumulating evidences indicate that a significant fraction of genetic alterations, including structural variants (SVs), can cause human disease by altering the function of non-coding regulatory elements, such as enhancers. In the case of SVs, described pathomechanisms include changes in enhancer dosage and long-range enhancer-gene communication. However, there is still a clear gap between the need to predict and interpret the medical impact of non-coding variants, and the existence of tools to properly perform these tasks. To reduce this gap, we have developed POSTRE (Prediction Of STRuctural variant Effects), a computational tool to predict the pathogenicity of SVs implicated in a broad range of human congenital disorders. By considering disease-relevant cellular contexts, POSTRE identifies SVs with either coding or long-range pathological consequences with high specificity and sensitivity. Furthermore, POSTRE not only identifies pathogenic SVs, but also predicts the disease-causative genes and the underlying pathological mechanism (e.g, gene deletion, enhancer disconnection, enhancer adoption, etc.). POSTRE is available at https://github.com/vicsanga/Postre.