Genomic views of distant-acting enhancers

A functional schizophrenia-associated genetic variant near the TSNARE1 and ADGRB1 genes

Recent collaborative genome-wide association studies (GWAS) have identified >200 independent loci contributing to risk for schizophrenia (SCZ). The genes closest to these loci have diverse functions, supporting the potential involvement of multiple relevant biological processes, yet there is no direct evidence that individual variants are functional or directly linked to specific genes. Nevertheless, overlap with certain epigenetic marks suggest that most GWAS-implicated variants are regulatory. Based on the strength of association with SCZ and the presence of regulatory epigenetic marks, we chose one such variant near TSNARE1 and ADGRB1, rs4129585, to test for functional potential and assay differences that may drive the pathogenicity of the risk allele. We observed that the variant-containing sequence drives reporter expression in relevant neuronal populations in zebrafish. Next, we introduced each allele into human induced pluripotent cells and differentiated four isogenic clones homozygous for the risk allele and five clones homozygous for the non-risk allele into neural progenitor cells. Employing RNA sequencing, we found that the two alleles yield significant transcriptional differences in the expression of 109 genes at a false discovery rate (FDR) of <0.05 and 259 genes at a FDR of <0.1. We demonstrate that these genes are highly interconnected in pathways enriched for synaptic proteins, axon guidance, and regulation of synapse assembly. Exploration of genes near rs4129585 suggests that this variant does not regulate TSNARE1 transcripts, as previously thought, but may regulate the neighboring ADGRB1, a regulator of synaptogenesis. Our results suggest that rs4129585 is a functional common variant that functions in specific pathways likely involved in SCZ risk.

Unravelling the Cerebellar Involvement in Autism Spectrum Disorders: Insights into Genetic Mechanisms and Developmental Pathways

Autism spectrum disorders (ASDs) are complex neurodevelopmental conditions characterized by deficits in social interaction and communication, as well as repetitive behaviors. Although the etiology of ASD is multifactorial, with both genetic and environmental factors contributing to its development, a strong genetic basis is widely recognized. Recent research has identified numerous genetic mutations and genomic rearrangements associated with ASD-characterizing genes involved in brain development. Alterations in developmental programs are particularly harmful during critical periods of brain development. Notably, studies have indicated that genetic disruptions occurring during the second trimester of pregnancy affect cortical development, while disturbances in the perinatal and early postnatal period affect cerebellar development. The developmental defects must be viewed in the context of the role of the cerebellum in cognitive processes, which is now well established. The present review emphasizes the genetic complexity and neuropathological mechanisms underlying ASD and aims to provide insights into the cerebellar involvement in the disorder, focusing on recent advances in the molecular landscape governing its development in humans. Furthermore, we highlight when and in which cerebellar neurons the ASD-associated genes may play a role in the development of cortico-cerebellar circuits. Finally, we discuss improvements in protocols for generating cerebellar organoids to recapitulate the long period of development and maturation of this organ. These models, if generated from patient-induced pluripotent stem cells (iPSC), could provide a valuable approach to elucidate the contribution of defective genes to ASD pathology and inform diagnostic and therapeutic strategies.

The broken Alzheimer's disease genome

The complex pathobiology of late-onset Alzheimer's disease (AD) poses significant challenges to therapeutic and preventative interventions. Despite these difficulties, genomics and related disciplines are allowing fundamental mechanistic insights to emerge with clarity, particularly with the introduction of high-resolution sequencing technologies. After all, the disrupted processes at the interface between DNA and gene expression, which we call the broken AD genome, offer detailed quantitative evidence unrestrained by preconceived notions about the disease. In addition to highlighting biological pathways beyond the classical pathology hallmarks, these advances have revitalized drug discovery efforts and are driving improvements in clinical tools. We review genetic, epigenomic, and gene expression findings related to AD pathogenesis and explore how their integration enables a better understanding of the multicellular imbalances contributing to this heterogeneous condition. The frontiers opening on the back of these research milestones promise a future of AD care that is both more personalized and predictive.

Distinct Histone H3 Lysine 27 Modifications Dictate Different Outcomes of Gene Transcription

Site-specific histone modifications have long been recognized to play an important role in directing gene transcription in chromatin in biology of health and disease. However, concrete illustration of how different histone modifications in a site-specific manner dictate gene transcription outcomes, as postulated in the influential "Histone code hypothesis", introduced by Allis and colleagues in 2000, has been lacking. In this review, we summarize our latest understanding of the dynamic regulation of gene transcriptional activation, silence, and repression in chromatin that is directed distinctively by histone H3 lysine 27 acetylation, methylation, and crotonylation, respectively. This represents a special example of a long-anticipated verification of the "Histone code hypothesis."

Cis-regulatory control of mammalian Trps1 gene expression

TRPS1 serves as the causative gene for tricho-rhino phalangeal syndrome, known for its craniofacial and skeletal abnormalities. The Trps1 gene encodes a protein that represses Wnt signaling through strong interactions with Wnt signaling inhibitors. The identification of genomic cis-acting regulatory sequences governing Trps1 expression is crucial for understanding its role in embryogenesis. Nevertheless, to date, no investigations have been conducted concerning these aspects of Trps1. To identify deeply conserved noncoding elements (CNEs) within the Trps1 locus, we employed a comparative genomics approach, utilizing slowly evolving fish such as coelacanth and spotted gar. These analyses resulted in the identification of eight CNEs in the intronic region of the Trps1 gene. Functional characterization of these CNEs in zebrafish revealed their regulatory potential in various tissues, including pectoral fins, heart, and pharyngeal arches. RNA in-situ hybridization experiments revealed concordance between the reporter expression pattern induced by the identified set of CNEs and the spatial expression pattern of the trps1 gene in zebrafish. Comparative in vivo data from zebrafish and mice for CNE7/hs919 revealed conserved functions of these enhancers. Each of these eight CNEs was further investigated in cell line-based reporter assays, revealing their repressive potential. Taken together, in vivo and in vitro assays suggest a context-dependent dual functionality for the identified set of Trps1-associated CNE enhancers. This functionally characterized set of CNE-enhancers will contribute to a more comprehensive understanding of the developmental roles of Trps1 and can aid in the identification of noncoding DNA variants associated with human diseases.

Affinity-optimizing enhancer variants disrupt development

Enhancers control the location and timing of gene expression and contain the majority of variants associated with disease1-3. The ZRS is arguably the most well-studied vertebrate enhancer and mediates the expression of Shh in the developing limb4. Thirty-one human single-nucleotide variants (SNVs) within the ZRS are associated with polydactyly4-6. However, how this enhancer encodes tissue-specific activity, and the mechanisms by which SNVs alter the number of digits, are poorly understood. Here we show that the ETS sites within the ZRS are low affinity, and identify a functional ETS site, ETS-A, with extremely low affinity. Two human SNVs and a synthetic variant optimize the binding affinity of ETS-A subtly from 15% to around 25% relative to the strongest ETS binding sequence, and cause polydactyly with the same penetrance and severity. A greater increase in affinity results in phenotypes that are more penetrant and more severe. Affinity-optimizing SNVs in other ETS sites in the ZRS, as well as in ETS, interferon regulatory factor (IRF), HOX and activator protein 1 (AP-1) sites within a wide variety of enhancers, cause gain-of-function gene expression. The prevalence of binding sites with suboptimal affinity in enhancers creates a vulnerability in genomes whereby SNVs that optimize affinity, even slightly, can be pathogenic. Searching for affinity-optimizing SNVs in genomes could provide a mechanistic approach to identify causal variants that underlie enhanceropathies.

Direct observation of a condensate effect on super-enhancer controlled gene bursting

Enhancers are distal DNA elements believed to loop and contact promoters to control gene expression. Recently, we found diffraction-sized transcriptional condensates at genes controlled by clusters of enhancers (super-enhancers). However, a direct function of endogenous condensates in controlling gene expression remains elusive. Here, we develop live-cell super-resolution and multi-color 3D-imaging approaches to investigate putative roles of endogenous condensates in the regulation of super-enhancer controlled gene Sox2. In contrast to enhancer distance, we find instead that the condensate's positional dynamics are a better predictor of gene expression. A basal gene bursting occurs when the condensate is far (>1 μm), but burst size and frequency are enhanced when the condensate moves in proximity (<1 μm). Perturbations of cohesin and local DNA elements do not prevent basal bursting but affect the condensate and its burst enhancement. We propose a three-way kissing model whereby the condensate interacts transiently with gene locus and regulatory DNA elements to control gene bursting.

Technologies to Study Genetics and Molecular Pathways

Over the last few decades, the study of congenital heart disease (CHD) has benefited from various model systems and the development of molecular biological techniques enabling the analysis of single gene as well as global effects. In this chapter, we first describe different models including CHD patients and their families, animal models ranging from invertebrates to mammals, and various cell culture systems. Moreover, techniques to experimentally manipulate these models are discussed. Second, we introduce cardiac phenotyping technologies comprising the analysis of mouse and cell culture models, live imaging of cardiogenesis, and histological methods for fixed hearts. Finally, the most important and latest molecular biotechniques are described. These include genotyping technologies, different applications of next-generation sequencing, and the analysis of transcriptome, epigenome, proteome, and metabolome. In summary, the models and technologies presented in this chapter are essential to study the function and development of the heart and to understand the molecular pathways underlying CHD.

A Functional Schizophrenia-associated genetic variant near the TSNARE1 and ADGRB1 genes

Recent collaborative genome wide association studies (GWAS) have identified >200 independent loci contributing to risk for schizophrenia (SCZ). The genes closest to these loci have diverse functions, supporting the potential involvement of multiple relevant biological processes; yet there is no direct evidence that individual variants are functional or directly linked to specific genes. Nevertheless, overlap with certain epigenetic marks suggest that most GWAS-implicated variants are regulatory. Based on the strength of association with SCZ and the presence of regulatory epigenetic marks, we chose one such variant near TSNARE1 and ADGRB1, rs4129585, to test for functional potential and assay differences that may drive the pathogenicity of the risk allele. We observed that the variant-containing sequence drives reporter expression in relevant neuronal populations in zebrafish. Next, we introduced each allele into human induced pluripotent cells and differentiated 4 isogenic clones homozygous for the risk allele and 5 clones homozygous for the non-risk allele into neural precursor cells. Employing RNA-seq, we found that the two alleles yield significant transcriptional differences in the expression of 109 genes at FDR <0.05 and 259 genes at FDR <0.1. We demonstrate that these genes are highly interconnected in pathways enriched for synaptic proteins, axon guidance, and regulation of synapse assembly. Exploration of genes near rs4129585 suggests that this variant does not regulate TSNARE1 transcripts, as previously thought, but may regulate the neighboring ADGRB1, a regulator of synaptogenesis. Our results suggest that rs4129585 is a functional common variant that functions in specific pathways likely involved in SCZ risk.

Genome-wide chromatin interaction profiling reveals a vital role of super-enhancers and rearrangements in host enhancer contacts during BmNPV infection

As influential regulatory elements in the genome, enhancers control gene expression under specific cellular conditions, and such connections are dynamic under different conditions. However, because of the lack of a genome-wide enhancer-gene connection map, the roles and regulatory pattern of enhancers were poorly investigated in insects, and the dynamic changes of enhancer contacts and functions under different conditions remain elusive. Here, combining Hi-C, ATAC-seq, and H3K27ac ChIP-seq data, we generate the genome-wide enhancer-gene map of silkworm and identify super-enhancers with a role in regulating the expression of vital genes related to cell state maintenance through a sophisticated interaction network. Additionally, a radical attenuation of chromatin interactions is found after infection of Bombyx mori nucleopolyhedrovirus (BmNPV), the main pathogen of silkworm, which directly rearranges the enhancer contacts. Such a rearrangement disturbs the intrinsic enhancer-gene connections in several antiviral genes, resulting in reduced expression of these genes, which accelerates viral infection. Overall, our results reveal the regulatory role of super-enhancers and shed new light on the mechanisms and dynamic changes of the genome-wide enhancer regulatory network.

Single-nucleotide variants within heart enhancers increase binding affinity and disrupt heart development

Transcriptional enhancers direct precise gene expression patterns during development and harbor the majority of variants associated with phenotypic diversity, evolutionary adaptations, and disease. Pinpointing which enhancer variants contribute to changes in gene expression and phenotypes is a major challenge. Here, we find that suboptimal or low-affinity binding sites are necessary for precise gene expression during heart development. Single-nucleotide variants (SNVs) can optimize the affinity of ETS binding sites, causing gain-of-function (GOF) gene expression, cell migration defects, and phenotypes as severe as extra beating hearts in the marine chordate Ciona robusta. In human induced pluripotent stem cell (iPSC)-derived cardiomyocytes, a SNV within a human GATA4 enhancer increases ETS binding affinity and causes GOF enhancer activity. The prevalence of suboptimal-affinity sites within enhancers creates a vulnerability whereby affinity-optimizing SNVs can lead to GOF gene expression, changes in cellular identity, and organismal-level phenotypes that could contribute to the evolution of novel traits or diseases.

CRISPR-based epigenome editing: mechanisms and applications

Epigenomic anomalies contribute significantly to the development of numerous human disorders. The development of epigenetic research tools is essential for understanding how epigenetic marks contribute to gene expression. A gene-editing technique known as CRISPR (clustered regularly interspaced short palindromic repeats) typically targets a particular DNA sequence using a guide RNA (gRNA). CRISPR/Cas9 technology has been remodeled for epigenome editing by generating a 'dead' Cas9 protein (dCas9) that lacks nuclease activity and juxtaposing it with an epigenetic effector domain. Based on fusion partners of dCas9, a specific epigenetic state can be achieved. CRISPR-based epigenome editing has widespread application in drug screening, cancer treatment and regenerative medicine. This paper discusses the tools developed for CRISPR-based epigenome editing and their applications.

Computational analysis of five neurodegenerative diseases reveals shared and specific genetic loci

Neurodegenerative diseases (ND) are heterogeneous disorders of the central nervous system that share a chronic and selective process of neuronal cell death. A computational approach to investigate shared genetic and specific loci was applied to 5 different ND: Amyotrophic lateral sclerosis (ALS), Alzheimer's disease (AD), Parkinson's disease (PD), Multiple sclerosis (MS), and Lewy body dementia (LBD). The datasets were analyzed separately, and then we compared the obtained results. For this purpose, we applied a genetic correlation analysis to genome-wide association datasets and revealed different genetic correlations with several human traits and diseases. In addition, a clumping analysis was carried out to identify SNPs genetically associated with each disease. We found 27 SNPs in AD, 6 SNPs in ALS, 10 SNPs in PD, 17 SNPs in MS, and 3 SNPs in LBD. Most of them are located in non-coding regions, with the exception of 5 SNPs on which a protein structure and stability prediction was performed to verify their impact on disease. Furthermore, an analysis of the differentially expressed miRNAs of the 5 examined pathologies was performed to reveal regulatory mechanisms that could involve genes associated with selected SNPs. In conclusion, the results obtained constitute an important step toward the discovery of diagnostic biomarkers and a better understanding of the diseases.

Context-specific synthetic T cell promoters from assembled transcriptional elements

Genetic engineering of human lymphocytes for therapeutic applications is constrained by a lack of transgene transcriptional control, resulting in a compromised therapeutic index. Incomplete understanding of transcriptional logic limits the rational design of contextually responsive genetic modules1. Here, we juxtaposed rationally curated transcriptional response element (TRE) oligonucleotides by random concatemerization to generate a library from which we selected context-specific inducible synthetic promoters (iSynPros). Through functional selection, we screened an iSynPro library for "IF-THEN" logic-gated transcriptional responses in human CD8+ T cells expressing a 4-1BB second generation chimeric antigen receptor (CAR). iSynPros exhibiting stringent off-states in quiescent T cells and CAR activation-dependent transcriptional responsiveness were cloned and subjected to TRE composition and pattern analysis, as well as performance in regulating candidate antitumor potency enhancement modules. These data reveal synthetic TRE grammar can mediate logic-gated transgene transcription in human T cells that, when applied to CAR T cell engineering, enhance potency and improve therapeutic indices.

Genomic Signatures of Local Adaptation under High Gene Flow in Lumpfish-Implications for Broodstock Provenance Sourcing and Larval Production

Aquaculture of the lumpfish (Cyclopterus lumpus L.) has become a large, lucrative industry owing to the escalating demand for "cleaner fish" to minimise sea lice infestations in Atlantic salmon mariculture farms. We used over 10K genome-wide single nucleotide polymorphisms (SNPs) to investigate the spatial patterns of genomic variation in the lumpfish along the coast of Norway and across the North Atlantic. Moreover, we applied three genome scans for outliers and two genotype-environment association tests to assess the signatures and patterns of local adaptation under extensive gene flow. With our 'global' sampling regime, we found two major genetic groups of lumpfish, i.e., the western and eastern Atlantic. Regionally in Norway, we found marginal evidence of population structure, where the population genomic analysis revealed a small portion of individuals with a different genetic ancestry. Nevertheless, we found strong support for local adaption under high gene flow in the Norwegian lumpfish and identified over 380 high-confidence environment-associated loci linked to gene sets with a key role in biological processes associated with environmental pressures and embryonic development. Our results bridge population genetic/genomics studies with seascape genomics studies and will facilitate genome-enabled monitoring of the genetic impacts of escapees and allow for genetic-informed broodstock selection and management in Norway.

Integrative PTEN Enhancer Discovery Reveals a New Model of Enhancer Organization

Enhancers possess both structural elements mediating promoter looping and functional elements mediating gene expression. Traditional models of enhancer-mediated gene regulation imply genomic overlap or immediate adjacency of these elements. We test this model by combining densely-tiled CRISPRa screening with nucleosome-resolution Region Capture Micro-C topology analysis. Using this integrated approach, we comprehensively define the cis-regulatory landscape for the tumor suppressor PTEN, identifying and validating 10 distinct enhancers and defining their 3D spatial organization. Unexpectedly, we identify several long-range functional enhancers whose promoter proximity is facilitated by chromatin loop anchors several kilobases away, and demonstrate that accounting for this spatial separation improves the computational prediction of validated enhancers. Thus, we propose a new model of enhancer organization incorporating spatial separation of essential functional and structural components.

PorcineAI-Enhancer: Prediction of Pig Enhancer Sequences Using Convolutional Neural Networks

Understanding the mechanisms of gene expression regulation is crucial in animal breeding. Cis-regulatory DNA sequences, such as enhancers, play a key role in regulating gene expression. Identifying enhancers is challenging, despite the use of experimental techniques and computational methods. Enhancer prediction in the pig genome is particularly significant due to the costliness of high-throughput experimental techniques. The study constructed a high-quality database of pig enhancers by integrating information from multiple sources. A deep learning prediction framework called PorcineAI-enhancer was developed for the prediction of pig enhancers. This framework employs convolutional neural networks for feature extraction and classification. PorcineAI-enhancer showed excellent performance in predicting pig enhancers, validated on an independent test dataset. The model demonstrated reliable prediction capability for unknown enhancer sequences and performed remarkably well on tissue-specific enhancer sequences.The study developed a deep learning prediction framework, PorcineAI-enhancer, for predicting pig enhancers. The model demonstrated significant predictive performance and potential for tissue-specific enhancers. This research provides valuable resources for future studies on gene expression regulation in pigs.

Molecular basis of FAAH-OUT-associated human pain insensitivity

Chronic pain affects millions of people worldwide and new treatments are needed urgently. One way to identify novel analgesic strategies is to understand the biological dysfunctions that lead to human inherited pain insensitivity disorders. Here we report how the recently discovered brain and dorsal root ganglia-expressed FAAH-OUT long non-coding RNA (lncRNA) gene, which was found from studying a pain-insensitive patient with reduced anxiety and fast wound healing, regulates the adjacent key endocannabinoid system gene FAAH, which encodes the anandamide-degrading fatty acid amide hydrolase enzyme. We demonstrate that the disruption in FAAH-OUT lncRNA transcription leads to DNMT1-dependent DNA methylation within the FAAH promoter. In addition, FAAH-OUT contains a conserved regulatory element, FAAH-AMP, that acts as an enhancer for FAAH expression. Furthermore, using transcriptomic analyses in patient-derived cells we have uncovered a network of genes that are dysregulated from disruption of the FAAH-FAAH-OUT axis, thus providing a coherent mechanistic basis to understand the human phenotype observed. Given that FAAH is a potential target for the treatment of pain, anxiety, depression and other neurological disorders, this new understanding of the regulatory role of the FAAH-OUT gene provides a platform for the development of future gene and small molecule therapies.

Cis-regulatory landscapes in the evolution and development of the mammalian skull

Extensive morphological variation found in mammals reflects the wide spectrum of their ecological adaptations. The highest morphological diversity is present in the craniofacial region, where geometry is mainly dictated by the bony skull. Mammalian craniofacial development represents complex multistep processes governed by numerous conserved genes that require precise spatio-temporal control. A central question in contemporary evolutionary biology is how a defined set of conserved genes can orchestrate formation of fundamentally different structures, and therefore how morphological variability arises. In principle, differential gene expression patterns during development are the source of morphological variation. With the emergence of multicellular organisms, precise regulation of gene expression in time and space is attributed to cis-regulatory elements. These elements contribute to higher-order chromatin structure and together with trans-acting factors control transcriptional landscapes that underlie intricate morphogenetic processes. Consequently, divergence in cis-regulation is believed to rewire existing gene regulatory networks and form the core of morphological evolution. This review outlines the fundamental principles of the genetic code and genomic regulation interplay during development. Recent work that deepened our comprehension of cis-regulatory element origin, divergence and function is presented here to illustrate the state-of-the-art research that uncovered the principles of morphological novelty. This article is part of the theme issue 'The mammalian skull: development, structure and function'.

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

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

Protein arginine methyltransferase 5 (Prmt5) localizes to chromatin loop anchors and modulates expression of genes at TAD boundaries during early adipogenesis

Protein arginine methyltransferase 5 (Prmt5) is an essential regulator of embryonic development and adult progenitor cell functions. Prmt5 expression is mis-regulated in many cancers, and the development of Prmt5 inhibitors as cancer therapeutics is an active area of research. Prmt5 functions via effects on gene expression, splicing, DNA repair, and other critical cellular processes. We examined whether Prmt5 functions broadly as a genome-wide regulator of gene transcription and higher-order chromatin interactions during the initial stages of adipogenesis using ChIP-Seq, RNA-seq, and Hi-C using 3T3-L1 cells, a frequently utilized model for adipogenesis. We observed robust genome-wide Prmt5 chromatin-binding at the onset of differentiation. Prmt5 localized to transcriptionally active genomic regions, acting as both a positive and a negative regulator. A subset of Prmt5 binding sites co-localized with mediators of chromatin organization at chromatin loop anchors. Prmt5 knockdown decreased insulation strength at the boundaries of topologically associating domains (TADs) adjacent to sites with Prmt5 and CTCF co-localization. Genes overlapping such weakened TAD boundaries showed transcriptional dysregulation. This study identifies Prmt5 as a broad regulator of gene expression, including regulation of early adipogenic factors, and reveals an unappreciated requirement for Prmt5 in maintaining strong insulation at TAD boundaries and overall chromatin organization.

The role of noncoding genetic variants in cardiomyopathy

Cardiomyopathies remain one of the leading causes of morbidity and mortality worldwide. Environmental risk factors and genetic predisposition account for most cardiomyopathy cases. As with all complex diseases, there are significant challenges in the interpretation of the molecular mechanisms underlying cardiomyopathy-associated genetic variants. Given the technical improvements and reduced costs of DNA sequence technologies, an increasing number of patients are now undergoing genetic testing, resulting in a continuously expanding list of novel mutations. However, many patients carry noncoding genetic variants, and although emerging evidence supports their contribution to cardiac disease, their role in cardiomyopathies remains largely understudied. In this review, we summarize published studies reporting on the association of different types of noncoding variants with various types of cardiomyopathies. We focus on variants within transcriptional enhancers, promoters, intronic sites, and untranslated regions that are likely associated with cardiac disease. Given the broad nature of this topic, we provide an overview of studies that are relatively recent and have sufficient evidence to support a significant degree of causality. We believe that more research with additional validation of noncoding genetic variants will provide further mechanistic insights on the development of cardiac disease, and noncoding variants will be increasingly incorporated in future genetic screening tests.

Atypical histone targets of PHD fingers

Plant homeodomain (PHD) fingers are structurally conserved zinc fingers that selectively bind unmodified or methylated at lysine 4 histone H3 tails. This binding stabilizes transcription factors and chromatin-modifying proteins at specific genomic sites, which is required for vital cellular processes, including gene expression and DNA repair. Several PHD fingers have recently been shown to recognize other regions of H3 or histone H4. In this review, we detail molecular mechanisms and structural features of the non-canonical histone recognition, discuss biological implications of the atypical interactions, highlight therapeutic potential of PHD fingers, and compare inhibition strategies.

Enhancers display constrained sequence flexibility and context-specific modulation of motif function

The information about when and where each gene is to be expressed is mainly encoded in the DNA sequence of enhancers, sequence elements that comprise binding sites (motifs) for different transcription factors (TFs). Most of the research on enhancer sequences has been focused on TF motif presence, whereas the enhancer syntax, that is, the flexibility of important motif positions and how the sequence context modulates the activity of TF motifs, remains poorly understood. Here, we explore the rules of enhancer syntax by a two-pronged approach in Drosophila melanogaster S2 cells: we (1) replace important TF motifs by all possible 65,536 eight-nucleotide-long sequences and (2) paste eight important TF motif types into 763 positions within 496 enhancers. These complementary strategies reveal that enhancers display constrained sequence flexibility and the context-specific modulation of motif function. Important motifs can be functionally replaced by hundreds of sequences constituting several distinct motif types, but these are only a fraction of all possible sequences and motif types. Moreover, TF motifs contribute with different intrinsic strengths that are strongly modulated by the enhancer sequence context (the flanking sequence, the presence and diversity of other motif types, and the distance between motifs), such that not all motif types can work in all positions. The context-specific modulation of motif function is also a hallmark of human enhancers, as we demonstrate experimentally. Overall, these two general principles of enhancer sequences are important to understand and predict enhancer function during development, evolution, and in disease.

Diverse logics and grammar encode notochord enhancers

The notochord is a defining feature of all chordates. The transcription factors Zic and ETS regulate enhancer activity within the notochord. We conduct high-throughput screens of genomic elements within developing Ciona embryos to understand how Zic and ETS sites encode notochord activity. Our screen discovers an enhancer located near Lama, a gene critical for notochord development. Reversing the orientation of an ETS site within this enhancer abolishes expression, indicating that enhancer grammar is critical for notochord activity. Similarly organized clusters of Zic and ETS sites occur within mouse and human Lama1 introns. Within a Brachyury (Bra) enhancer, FoxA and Bra, in combination with Zic and ETS binding sites, are necessary and sufficient for notochord expression. This binding site logic also occurs within other Ciona and vertebrate Bra enhancers. Collectively, this study uncovers the importance of grammar within notochord enhancers and discovers signatures of enhancer logic and grammar conserved across chordates.

Splashed E-box and AP-1 motifs cooperatively drive regeneration response and shape regeneration abilities

Injury triggers a genetic program that induces gene expression for regeneration. Recent studies have identified regeneration-response enhancers (RREs); however, it remains unclear whether a common mechanism operates in these RREs. We identified three RREs from the zebrafish fn1b promoter by searching for conserved sequences within the surrounding genomic regions of regeneration-induced genes and performed a transgenic assay for regeneration response. Two regions contained in the transposons displayed RRE activity when combined with the -0.7 kb fn1b promoter. Another non-transposon element functioned as a stand-alone enhancer in combination with a minimum promoter. By searching for transcription factor-binding motifs and validation by transgenic assays, we revealed that the cooperation of E-box and activator protein 1 motifs is necessary and sufficient for regenerative response. Such RREs respond to variety of tissue injuries, including those in the zebrafish heart and Xenopus limb buds. Our findings suggest that the fidelity of regeneration response is ensured by the two signals evoked by tissue injuries. It is speculated that a large pool of potential enhancers in the genome has helped shape the regenerative capacities during evolution.

Molecular characterization of Richter syndrome identifies de novo diffuse large B-cell lymphomas with poor prognosis

Richter syndrome (RS) is the transformation of chronic lymphocytic leukemia (CLL) into aggressive lymphoma, most commonly diffuse large B-cell lymphoma (DLBCL). We characterize 58 primary human RS samples by genome-wide DNA methylation and whole-transcriptome profiling. Our comprehensive approach determines RS DNA methylation profile and unravels a CLL epigenetic imprint, allowing CLL-RS clonal relationship assessment without the need of the initial CLL tumor DNA. DNA methylation- and transcriptomic-based classifiers were developed, and testing on landmark DLBCL datasets identifies a poor-prognosis, activated B-cell-like DLBCL subset in 111/1772 samples. The classification robustly identifies phenotypes very similar to RS with a specific genomic profile, accounting for 4.3-8.3% of de novo DLBCLs. In this work, RS multi-omics characterization determines oncogenic mechanisms, establishes a surrogate marker for CLL-RS clonal relationship, and provides a clinically relevant classifier for a subset of primary "RS-type DLBCL" with unfavorable prognosis.

ChromLoops: a comprehensive database for specific protein-mediated chromatin loops in diverse organisms

Chromatin loops (or chromatin interactions) are important elements of chromatin structures. Disruption of chromatin loops is associated with many diseases, such as cancer and polydactyly. A few methods, including ChIA-PET, HiChIP and PLAC-Seq, have been proposed to detect high-resolution, specific protein-mediated chromatin loops. With rapid progress in 3D genomic research, ChIA-PET, HiChIP and PLAC-Seq datasets continue to accumulate, and effective collection and processing for these datasets are urgently needed. Here, we developed a comprehensive, multispecies and specific protein-mediated chromatin loop database (ChromLoops, https://3dgenomics.hzau.edu.cn/chromloops), which integrated 1030 ChIA-PET, HiChIP and PLAC-Seq datasets from 13 species, and documented 1 491 416 813 high-quality chromatin loops. We annotated genes and regions overlapping with chromatin loop anchors with rich functional annotations, such as regulatory elements (enhancers, super-enhancers and silencers), variations (common SNPs, somatic SNPs and eQTLs), and transcription factor binding sites. Moreover, we identified genes with high-frequency chromatin interactions in the collected species. In particular, we identified genes with high-frequency interactions in cancer samples. We hope that ChromLoops will provide a new platform for studying chromatin interaction regulation in relation to biological processes and disease.

Current advances in primate genomics: novel approaches for understanding evolution and disease

Primate genomics holds the key to understanding fundamental aspects of human evolution and disease. However, genetic diversity and functional genomics data sets are currently available for only a few of the more than 500 extant primate species. Concerted efforts are under way to characterize primate genomes, genetic polymorphism and divergence, and functional landscapes across the primate phylogeny. The resulting data sets will enable the connection of genotypes to phenotypes and provide new insight into aspects of the genetics of primate traits, including human diseases. In this Review, we describe the existing genome assemblies as well as genetic variation and functional genomic data sets. We highlight some of the challenges with sample acquisition. Finally, we explore how technological advances in single-cell functional genomics and induced pluripotent stem cell-derived organoids will facilitate our understanding of the molecular foundations of primate biology.

Bioinformatics Pipelines for Identification of Super-Enhancers and 3D Chromatin Contacts

Precise regulation of gene expression is integral in development. Emerging studies have highlighted that super-enhancers (SEs), which are clusters of multiple enhancers, play critical roles in regulating cell type-specific gene expression via 3D chromatin, thereby defining the cellular identities of given cells. Here we provide optimized bioinformatics pipelines to identify SEs and 3D chromatin contacts. Our pipelines encompass the processing of chromatin immunoprecipitation sequencing (ChIP-seq) data to identify SEs and the processing of genome-wide chromosome conformation capture (Hi-C) data. We can then infer long-range chromatin contacts between SEs and other genomic regions. This integrative computational approach, which can be applied to CUT&RUN and CUT&Tag, alternative technologies to ChIP-seq, allows us to identify genomic locations of SEs and their 3D genome configuration, whereby multiple SEs act in concert. We show an analysis of mouse spermatogenesis as an example of this application.

The landscape of enhancer RNA identify prognosis-related molecular subtypes in gastric cancer

BACKGROUND

Enhancer RNAs (eRNAs), the transcriptional products of active enhancers, are of great significance in the initial progression of cancers. However, the biological function and bioinformatics profiles of eRNA in gastric cancer remains largely enigmatic.

METHODS

Firstly, STAD were clustered into three subtypes with the data of eRNA expression from TCeA. Then we explored the difference of the tumor immune microenvironment, transcription levels, and transcription regulation among the three clusters. Finally, samples collected from 12 patients diagnosed with STAD were used to conduct qRT-PCR, verifying the conclusion based on network database.

RESULTS

The three clusters were detected to have different tumor microenvironments: Cluster A has an immune "cold" microenvironment. While cluster B features as more infiltration of immune cells, accompanied with higher expression of immune checkpoints such as PDCD1, LAG3, and TIGIT. Besides, Cluster C shows a higher stromal feature with B lineage, neutrophils, and fibroblasts. Further analyses indicated that CpG island methylation level of Cluster B is different from the other two clusters. Meanwhile, Cluster A and B showed significant enrichment of TP53 and KRAS mutation respectively while Cluster C has higher tumor mutation burden (TMB) and microsatellite instability (MSI). With the elaboration of transcriptional regulation of epigenetic clustering, we detected that Cluster A enriched in epithelial phenotype pathways. Cluster B enriched in cell-cell adhesion. Cluster C enriched in fibroblast proliferation. The clinical cohort show that Cluster B patients have lower interstitial cell characteristics and CAF infiltration.

CONCLUSION

We identified three unique epigenetic clusters of STAD through the differential activation of super-enhancers, and identified Cluster B with a higher immune infiltrating and a better prognosis, which provides a novel understanding of eRNAs and potential clinical applicability of eRNA-based molecular subtypes in gastric cancer.

Integration of 3D genome topology and local chromatin features uncovers enhancers underlying craniofacial-specific cartilage defects

Aberrations in tissue-specific enhancers underlie many developmental defects. Disrupting a noncoding region distal from the human SOX9 gene causes the Pierre Robin sequence (PRS) characterized by the undersized lower jaw. Such a craniofacial-specific defect has been previously linked to enhancers transiently active in cranial neural crest cells (CNCCs). We demonstrate that the PRS region also strongly regulates Sox9 in CNCC-derived Meckel's cartilage (MC), but not in limb cartilages, even after decommissioning of CNCC enhancers. Such an MC-specific regulatory effect correlates with the MC-specific chromatin contacts between the PRS region and Sox9, highlighting the importance of lineage-dependent chromatin topology in instructing enhancer usage. By integrating the enhancer signatures and chromatin topology, we uncovered >10,000 enhancers that function differentially between MC and limb cartilages and demonstrated their association with human diseases. Our findings provide critical insights for understanding the choreography of gene regulation during development and interpreting the genetic basis of craniofacial pathologies.

XLID syndrome gene Med12 promotes Ig isotype switching through chromatin modification and enhancer RNA regulation

The transcriptional coactivator Med12 regulates gene expression through its kinase module. Here, we show a kinase module-independent function of Med12 in CSR. Med12 is essential for super-enhancer activation by collaborating with p300-Jmjd6/Carm1 coactivator complexes. Med12 loss decreases H3K27 acetylation and eRNA transcription with concomitant impairment of AID-induced DNA breaks, S-S synapse formation, and 3'RR-Eμ interaction. CRISPR-dCas9-mediated enhancer activation reestablishes the epigenomic and transcriptional hallmarks of the super-enhancer and fully restores the Med12 depletion defects. Moreover, 3'RR-derived eRNAs are critical for promoting S region epigenetic regulation, synapse formation, and recruitment of Med12 and AID to the IgH locus. We find that XLID syndrome-associated Med12 mutations are defective in both 3'RR eRNA transcription and CSR, suggesting that B and neuronal cells may have cell-specific super-enhancer dysfunctions. We conclude that Med12 is essential for IgH 3'RR activation/eRNA transcription and plays a central role in AID-induced antibody gene diversification and genomic instability in B cells.

Shared cis-regulatory modules control expression of the tandem paralogs midline and H15 in the follicular epithelium

The posterior end of the follicular epithelium is patterned by midline (MID) and its paralog H15, the Drosophila homologs of the mammalian Tbx20 transcription factor. We have previously identified two cis-regulatory modules (CRMs) that recapitulate the endogenous pattern of mid in the follicular epithelium. Here, using CRISPR/Cas9 genome editing, we demonstrate redundant activity of these mid CRMs. Although the deletion of either CRM alone generated marginal change in mid expression, the deletion of both CRMs reduced expression by 60%. Unexpectedly, the deletion of the 5' proximal CRM of mid eliminated H15 expression. Interestingly, expression of these paralogs in other tissues remained unaffected in the CRM deletion backgrounds. These results suggest that the paralogs are regulated by a shared CRM that coordinates gene expression during posterior fate determination. The consistent overlapping expression of mid and H15 in various tissues may indicate that the paralogs could also be under shared regulation by other CRMs in these tissues.

Construction of a transposase accessible chromatin landscape reveals chromatin state of repeat elements and potential causal variant for complex traits in pigs

Background

A comprehensive landscape of chromatin states for multiple mammalian tissues is essential for elucidating the molecular mechanism underlying regulatory variants on complex traits. However, the genome-wide chromatin accessibility has been only reported in limited tissue types in pigs.

Results

Here we report a genome-wide landscape of chromatin accessibility of 20 tissues in two female pigs at ages of 6 months using ATAC-seq, and identified 557,273 merged peaks, which greatly expanded the pig regulatory element repository. We revealed tissue-specific regulatory elements which were associated with tissue-relevant biological functions. We identified both positive and negative significant correlations between the regulatory elements and gene transcripts, which showed distinct distributions in terms of their strength and distances from corresponding genes. We investigated the presence of transposable elements (TEs) in open chromatin regions across all tissues, these included identifications of porcine endogenous retroviruses (PERVs) exhibiting high accessibility in liver and homology of porcine specific virus sequences to universally accessible transposable elements. Furthermore, we prioritized a potential causal variant for polyunsaturated fatty acid in the muscle.

Conclusions

Our data provides a novel multi-tissues accessible chromatin landscape that serve as an important resource for interpreting regulatory sequences in tissue-specific and conserved biological functions, as well as regulatory variants of loci associated with complex traits in pigs.

Supplementary Information

The online version contains supplementary material available at 10.1186/s40104-022-00767-3.

"Stripe" transcription factors provide accessibility to co-binding partners in mammalian genomes

Regulatory elements activate promoters by recruiting transcription factors (TFs) to specific motifs. Notably, TF-DNA interactions often depend on cooperativity with colocalized partners, suggesting an underlying cis-regulatory syntax. To explore TF cooperativity in mammals, we analyze ∼500 mouse and human primary cells by combining an atlas of TF motifs, footprints, ChIP-seq, transcriptomes, and accessibility. We uncover two TF groups that colocalize with most expressed factors, forming stripes in hierarchical clustering maps. The first group includes lineage-determining factors that occupy DNA elements broadly, consistent with their key role in tissue-specific transcription. The second one, dubbed universal stripe factors (USFs), comprises ∼30 SP, KLF, EGR, and ZBTB family members that recognize overlapping GC-rich sequences in all tissues analyzed. Knockouts and single-molecule tracking reveal that USFs impart accessibility to colocalized partners and increase their residence time. Mammalian cells have thus evolved a TF superfamily with overlapping DNA binding that facilitate chromatin accessibility.

A machine learning technique for identifying DNA enhancer regions utilizing CIS-regulatory element patterns

Enhancers regulate gene expression, by playing a crucial role in the synthesis of RNAs and proteins. They do not directly encode proteins or RNA molecules. In order to control gene expression, it is important to predict enhancers and their potency. Given their distance from the target gene, lack of common motifs, and tissue/cell specificity, enhancer regions are thought to be difficult to predict in DNA sequences. Recently, a number of bioinformatics tools were created to distinguish enhancers from other regulatory components and to pinpoint their advantages. However, because the quality of its prediction method needs to be improved, its practical application value must also be improved. Based on nucleotide composition and statistical moment-based features, the current study suggests a novel method for identifying enhancers and non-enhancers and evaluating their strength. The proposed study outperformed state-of-the-art techniques using fivefold and tenfold cross-validation in terms of accuracy. The accuracy from the current study results in 86.5% and 72.3% in enhancer site and its strength prediction respectively. The results of the suggested methodology point to the potential for more efficient and successful outcomes when statistical moment-based features are used. The current study's source code is available to the research community at https://github.com/csbioinfopk/enpred.

Transcriptomics, regulatory syntax, and enhancer identification in mesoderm-induced ESCs at single-cell resolution

Embryonic stem cells (ESCs) can adopt lineage-specific gene-expression programs by stepwise exposure to defined factors, resulting in the generation of functional cell types. Bulk and single-cell-based assays were employed to catalog gene expression, histone modifications, chromatin conformation, and accessibility transitions in ESC populations and individual cells acquiring a presomitic mesoderm fate and undergoing further specification toward myogenic and neurogenic lineages. These assays identified cis-regulatory regions and transcription factors presiding over gene-expression programs occurring at defined ESC transitions and revealed the presence of heterogeneous cell populations within discrete ESC developmental stages. The datasets were employed to identify previously unappreciated genomic elements directing the initial activation of Pax7 and myogenic and neurogenic gene-expression programs. This study provides a resource for the discovery of genomic and transcriptional features of pluripotent, mesoderm-induced ESCs and ESC-derived cell lineages.

Temporal analysis of enhancers during mouse cerebellar development reveals dynamic and novel regulatory functions

We have identified active enhancers in the mouse cerebellum at embryonic and postnatal stages which provides a view of novel enhancers active during cerebellar development. The majority of cerebellar enhancers have dynamic activity between embryonic and postnatal development. Cerebellar enhancers were enriched for neural transcription factor binding sites with temporally specific expression. Putative gene targets displayed spatially restricted expression patterns, indicating cell-type specific expression regulation. Functional analysis of target genes indicated that enhancers regulate processes spanning several developmental epochs such as specification, differentiation and maturation. We use these analyses to discover one novel regulator and one novel marker of cerebellar development: Bhlhe22 and Pax3, respectively. We identified an enrichment of de novo mutations and variants associated with autism spectrum disorder in cerebellar enhancers. Furthermore, by comparing our data with relevant brain development ENCODE histone profiles and cerebellar single-cell datasets we have been able to generalize and expand on the presented analyses, respectively. We have made the results of our analyses available online in the Developing Mouse Cerebellum Enhancer Atlas (https://goldowitzlab.shinyapps.io/developing_mouse_cerebellum_enhancer_atlas/), where our dataset can be efficiently queried, curated and exported by the scientific community to facilitate future research efforts. Our study provides a valuable resource for studying the dynamics of gene expression regulation by enhancers in the developing cerebellum and delivers a rich dataset of novel gene-enhancer associations providing a basis for future in-depth studies in the cerebellum.

Multi-layered transcriptional control of cranial neural crest development

The neural crest (NC) is an emblematic population of embryonic stem-like cells with remarkable migratory ability. These distinctive attributes have inspired the curiosity of developmental biologists for over 150 years, however only recently the regulatory mechanisms controlling the complex features of the NC have started to become elucidated at genomic scales. Regulatory control of NC development is achieved through combinatorial transcription factor binding and recruitment of associated transcriptional complexes to distal cis-regulatory elements. Together, they regulate when, where and to what extent transcriptional programmes are actively deployed, ultimately shaping ontogenetic processes. Here, we discuss how transcriptional networks control NC ontogeny, with a special emphasis on the molecular mechanisms underlying specification of the cephalic NC. We also cover emerging properties of transcriptional regulation revealed in diverse developmental systems, such as the role of three-dimensional conformation of chromatin, and how they are involved in the regulation of NC ontogeny. Finally, we highlight how advances in deciphering the NC transcriptional network have afforded new insights into the molecular basis of human diseases.

Bridging between Mouse and Human Enhancer-Promoter Long-Range Interactions in Neural Stem Cells, to Understand Enhancer Function in Neurodevelopmental Disease

Non-coding variation in complex human disease has been well established by genome-wide association studies, and is thought to involve regulatory elements, such as enhancers, whose variation affects the expression of the gene responsible for the disease. The regulatory elements often lie far from the gene they regulate, or within introns of genes differing from the regulated gene, making it difficult to identify the gene whose function is affected by a given enhancer variation. Enhancers are connected to their target gene promoters via long-range physical interactions (loops). In our study, we re-mapped, onto the human genome, more than 10,000 enhancers connected to promoters via long-range interactions, that we had previously identified in mouse brain-derived neural stem cells by RNApolII-ChIA-PET analysis, coupled to ChIP-seq mapping of DNA/chromatin regions carrying epigenetic enhancer marks. These interactions are thought to be functionally relevant. We discovered, in the human genome, thousands of DNA regions syntenic with the interacting mouse DNA regions (enhancers and connected promoters). We further annotated these human regions regarding their overlap with sequence variants (single nucleotide polymorphisms, SNPs; copy number variants, CNVs), that were previously associated with neurodevelopmental disease in humans. We document various cases in which the genetic variant, associated in humans to neurodevelopmental disease, affects an enhancer involved in long-range interactions: SNPs, previously identified by genome-wide association studies to be associated with schizophrenia, bipolar disorder, and intelligence, are located within our human syntenic enhancers, and alter transcription factor recognition sites. Similarly, CNVs associated to autism spectrum disease and other neurodevelopmental disorders overlap with our human syntenic enhancers. Some of these enhancers are connected (in mice) to homologs of genes already associated to the human disease, strengthening the hypothesis that the gene is indeed involved in the disease. Other enhancers are connected to genes not previously associated with the disease, pointing to their possible pathogenetic involvement. Our observations provide a resource for further exploration of neural disease, in parallel with the now widespread genome-wide identification of DNA variants in patients with neural disease.

Identification of enhancers responsible for the coordinated expression of myosin heavy chain isoforms in skeletal muscle

Background

Skeletal muscles consist of fibers of differing contractility and metabolic properties, which are primarily determined by the content of myosin heavy chain (MYH) isoforms (MYH7, MYH2, MYH1, and MYH4). The regulation of Myh genes transcription depends on three-dimensional chromatin conformation interaction, but the mechanistic details remain to be determined.

Results

In this study, we characterized the interaction profiles of Myh genes using 4C-seq (circular chromosome conformation capture coupled to high-throughput sequencing). The interaction profile of Myh genes changed between fast quadriceps and slow soleus muscles. Combining chromatin immunoprecipitation-sequencing (ChIP-seq) and transposase accessible chromatin with high-throughput sequencing (ATAC-seq), we found that a 38 kb intergenic region interacting simultaneously with fast Myh genes promoters controlled the coordinated expression of fast Myh genes. We also identified four active enhancers of Myh7, and revealed that binding of MYOG and MYOD increased the activity of Myh7 enhancers.

Conclusions

This study provides new insight into the chromatin interactions that regulate Myh genes expression.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-022-08737-9.

StackEPI: identification of cell line-specific enhancer-promoter interactions based on stacking ensemble learning

Background

Understanding the regulatory role of enhancer–promoter interactions (EPIs) on specific gene expression in cells contributes to the understanding of gene regulation, cell differentiation, etc., and its identification has been a challenging task. On the one hand, using traditional wet experimental methods to identify EPIs often means a lot of human labor and time costs. On the other hand, although the currently proposed computational methods have good recognition effects, they generally require a long training time.

Results

In this study, we studied the EPIs of six human cell lines and designed a cell line-specific EPIs prediction method based on a stacking ensemble learning strategy, which has better prediction performance and faster training speed, called StackEPI. Specifically, by combining different encoding schemes and machine learning methods, our prediction method can extract the cell line-specific effective information of enhancer and promoter gene sequences comprehensively and in many directions, and make accurate recognition of cell line-specific EPIs. Ultimately, the source code to implement StackEPI and experimental data involved in the experiment are available at https://github.com/20032303092/StackEPI.git.

Conclusions

The comparison results show that our model can deliver better performance on the problem of identifying cell line-specific EPIs and outperform other state-of-the-art models. In addition, our model also has a more efficient computation speed.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12859-022-04821-9.

Marsupials and Multi-Omics: Establishing New Comparative Models of Neural Crest Patterning and Craniofacial Development

Studies across vertebrates have revealed significant insights into the processes that drive craniofacial morphogenesis, yet we still know little about how distinct facial morphologies are patterned during development. Studies largely point to evolution in GRNs of cranial progenitor cell types such as neural crest cells, as the major driver underlying adaptive cranial shapes. However, this hypothesis requires further validation, particularly within suitable models amenable to manipulation. By utilizing comparative models between related species, we can begin to disentangle complex developmental systems and identify the origin of species-specific patterning. Mammals present excellent evolutionary examples to scrutinize how these differences arise, as sister clades of eutherians and marsupials possess suitable divergence times, conserved cranial anatomies, modular evolutionary patterns, and distinct developmental heterochrony in their NCC behaviours and craniofacial patterning. In this review, I lend perspectives into the current state of mammalian craniofacial biology and discuss the importance of establishing a new marsupial model, the fat-tailed dunnart, for comparative research. Through detailed comparisons with the mouse, we can begin to decipher mammalian conserved, and species-specific processes and their contribution to craniofacial patterning and shape disparity. Recent advances in single-cell multi-omics allow high-resolution investigations into the cellular and molecular basis of key developmental processes. As such, I discuss how comparative evolutionary application of these tools can provide detailed insights into complex cellular behaviours and expression dynamics underlying adaptive craniofacial evolution. Though in its infancy, the field of "comparative evo-devo-omics" presents unparalleled opportunities to precisely uncover how phenotypic differences arise during development.

A synthetic synthesis to explore animal evolution and development

Identifying the general principles by which genotypes are converted into phenotypes remains a challenge in the post-genomic era. We still lack a predictive understanding of how genes shape interactions among cells and tissues in response to signalling and environmental cues, and hence how regulatory networks generate the phenotypic variation required for adaptive evolution. Here, we discuss how techniques borrowed from synthetic biology may facilitate a systematic exploration of evolvability across biological scales. Synthetic approaches permit controlled manipulation of both endogenous and fully engineered systems, providing a flexible platform for investigating causal mechanisms in vivo. Combining synthetic approaches with multi-level phenotyping (phenomics) will supply a detailed, quantitative characterization of how internal and external stimuli shape the morphology and behaviour of living organisms. We advocate integrating high-throughput experimental data with mathematical and computational techniques from a variety of disciplines in order to pursue a comprehensive theory of evolution.

This article is part of the theme issue ‘Genetic basis of adaptation and speciation: from loci to causative mutations’.

A high-density genome-wide association with absolute blood monocyte count in domestic sheep identifies novel loci

Monocytes are a core component of the immune system that arise from bone marrow and differentiate into cells responsible for phagocytosis and antigen presentation. Their derivatives are often responsible for the initiation of the adaptive immune response. Monocytes and macrophages are central in both controlling and propagating infectious diseases such as infection by Coxiella burnetii and small ruminant lentivirus in sheep. Genotypes from 513 Rambouillet, Polypay, and Columbia sheep (Ovis aries) were generated using the Ovine SNP50 BeadChip. Of these sheep, 222 animals were subsequently genotyped with the Ovine Infinium^® HD SNP BeadChip to increase SNP coverage. Data from the 222 HD genotyped sheep were combined with the data from an additional 258 unique sheep to form a 480-sheep reference panel; this panel was used to impute the low-density genotypes to the HD genotyping density. Then, a genome-wide association analysis was conducted to identify loci associated with absolute monocyte counts from blood. The analysis used a single-locus mixed linear model implementing EMMAX with age and ten principal components as fixed effects. Two genome-wide significant peaks (p < 5x10^-7) were identified on chromosomes 9 and 1, and ten genome-wide suggestive peaks (p < 1x10^-5) were identified on chromosomes 1, 2, 3, 4, 9, 10, 15, and 16. The identified loci were within or near genes including KCNK9, involved into cytokine production, LY6D, a member of a superfamily of genes, some of which subset monocyte lineages, and HMGN1, which encodes a chromatin regulator associated with myeloid cell differentiation. Further investigation of these loci is being conducted to understand their contributions to monocyte counts. Investigating the genetic basis of monocyte lineages and numbers may in turn provide information about pathogens of veterinary importance and elucidate fundamental immunology.

ΔNp63 regulates a common landscape of enhancer associated genes in non-small cell lung cancer

Distinct lung stem cells give rise to lung adenocarcinoma (LUAD) and squamous cell carcinoma (LUSC). ΔNp63, the p53 family member and p63 isoform, guides the maturation of these stem cells through the regulation of their self-renewal and terminal differentiation; however, the underlying mechanistic role regulated by ∆Np63 in lung cancer development has remained elusive. By utilizing a ΔNp63-specific conditional knockout mouse model and xenograft models of LUAD and LUSC, we found that ∆Np63 promotes non-small cell lung cancer by maintaining the lung stem cells necessary for lung cancer cell initiation and progression in quiescence. ChIP-seq analysis of lung basal cells, alveolar type 2 (AT2) cells, and LUAD reveals robust ∆Np63 regulation of a common landscape of enhancers of cell identity genes. Importantly, one of these genes, BCL9L, is among the enhancer associated genes regulated by ∆Np63 in Kras-driven LUAD and mediates the oncogenic effects of ∆Np63 in both LUAD and LUSC. Accordingly, high BCL9L levels correlate with poor prognosis in LUAD patients. Taken together, our findings provide a unifying oncogenic role for ∆Np63 in both LUAD and LUSC through the regulation of a common landscape of enhancer associated genes.

The mechanistic role regulated by the oncogene ∆Np63 in lung cancer development is currently unclear. Here, the authors show that ΔNp63 is pro-tumorigenic in lung adenocarcinoma as well as squamous cell carcinoma, and maintains lung cancer progenitor cells via regulation of super-enhancer-associated genes, including BCL9L

Finding and Verifying Enhancers for Endothelial-Expressed Genes

Identification and analysis of enhancers for endothelial-expressed genes can provide crucial information regarding their upstream transcriptional regulators. However, enhancer identification can be challenging, particularly for people with limited access or experience of bioinformatics, and transgenic analysis of enhancer activity patterns can be prohibitively expensive. Here we describe how to use publicly available datasets displayed on the UCSC Genome Browser to identify putative endothelial enhancers for mammalian genes. Furthermore, we detail how to utilize mosaic Tol2-mediated transgenesis in zebrafish to verify whether a putative enhancer is capable of directing endothelial-specific patterns of gene expression.

Regulation of the Late Onset alzheimer's Disease Associated HLA-DQA1/DRB1 Expression

(Genome-wide Association Studies) GWAS have identified ∼42 late-onset Alzheimer's disease (LOAD)-associated loci, each of which contains multiple single nucleotide polymorphisms (SNPs) in linkage disequilibrium (LD) and most of these SNPs are in the non-coding region of human genome. However, how these SNPs regulate risk gene expression remains unknown. In this work, by using a set of novel techniques, we identified 6 functional SNPs (fSNPs) rs9271198, rs9271200, rs9281945, rs9271243, and rs9271247 on the LOAD-associated HLA-DRB1/DQA1 locus and 42 proteins specifically binding to five of these 6 fSNPs. As a proof of evidence, we verified the allele-specific binding of GATA2 and GATA3, ELAVL1 and HNRNPA0, ILF2 and ILF3, NFIB and NFIC, as well as CUX1 to these five fSNPs, respectively. Moreover, we demonstrate that all these nine proteins regulate the expression of both HLA-DQA1 and HLA-DRB1 in human microglial cells. The contribution of HLA class II to the susceptibility of LOAD is discussed.

Targeting SWI/SNF ATPases in enhancer-addicted prostate cancer

The switch/sucrose non-fermentable (SWI/SNF) complex has a crucial role in chromatin remodelling1 and is altered in over 20% of cancers2,3. Here we developed a proteolysis-targeting chimera (PROTAC) degrader of the SWI/SNF ATPase subunits, SMARCA2 and SMARCA4, called AU-15330. Androgen receptor (AR)+ forkhead box A1 (FOXA1)+ prostate cancer cells are exquisitely sensitive to dual SMARCA2 and SMARCA4 degradation relative to normal and other cancer cell lines. SWI/SNF ATPase degradation rapidly compacts cis-regulatory elements bound by transcription factors that drive prostate cancer cell proliferation, namely AR, FOXA1, ERG and MYC, which dislodges them from chromatin, disables their core enhancer circuitry, and abolishes the downstream oncogenic gene programs. SWI/SNF ATPase degradation also disrupts super-enhancer and promoter looping interactions that wire supra-physiologic expression of the AR, FOXA1 and MYC oncogenes themselves. AU-15330 induces potent inhibition of tumour growth in xenograft models of prostate cancer and synergizes with the AR antagonist enzalutamide, even inducing disease remission in castration-resistant prostate cancer (CRPC) models without toxicity. Thus, impeding SWI/SNF-mediated enhancer accessibility represents a promising therapeutic approach for enhancer-addicted cancers.