EVOLUTIONARY CO-OPTION OF AN ANCESTRAL CLOACAL REGULATORY LANDSCAPE DURING THE EMERGENCE OF DIGITS AND GENITALS

The transition from fins to limbs has been a rich source of discussion for more than a century. One open and important issue is understanding how the mechanisms that pattern digits arose during vertebrate evolution. In this context, the analysis of Hox gene expression and functions to infer evolutionary scenarios has been a productive approach to explain the changes in organ formation, particularly in limbs. In tetrapods, the transcription of Hoxd genes in developing digits depends on a well-characterized set of enhancers forming a large regulatory landscape1,2. This control system has a syntenic counterpart in zebrafish, even though they lack bona fide digits, suggestive of deep homology3 between distal fin and limb developmental mechanisms. We tested the global function of this landscape to assess ancestry and source of limb and fin variation. In contrast to results in mice, we show here that the deletion of the homologous control region in zebrafish has a limited effect on the transcription of hoxd genes during fin development. However, it fully abrogates hoxd expression within the developing cloaca, an ancestral structure related to the mammalian urogenital sinus. We show that similar to the limb, Hoxd gene function in the urogenital sinus of the mouse also depends on enhancers located in this same genomic domain. Thus, we conclude that the current regulation underlying Hoxd gene expression in distal limbs was co-opted in tetrapods from a preexisting cloacal program. The orthologous chromatin domain in fishes may illustrate a rudimentary or partial step in this evolutionary co-option.

Shaping gene expression and its evolution by chromatin architecture and enhancer activity

Transcriptional regulation plays a pivotal role in orchestrating the intricate genetic programs governing embryonic development. The expression of developmental genes relies on the combined activity of several cis-regulatory elements (CREs), such as enhancers and silencers, which can be located at long linear distances from the genes that they regulate and that interact with them through establishment of chromatin loops. Mutations affecting their activity or interaction with their target genes can lead to developmental disorders and are thought to have importantly contributed to the evolution of the animal body plan. The income of next-generation-sequencing approaches has allowed identifying over a million of sequences with putative regulatory potential in the human genome. Characterizing their function and establishing gene-CREs maps is essential to decode the logic governing developmental gene expression and is one of the major challenges of the post-genomic era. Chromatin 3D organization plays an essential role in determining how CREs specifically contact their target genes while avoiding deleterious off-target interactions. Our understanding of these aspects has greatly advanced with the income of chromatin conformation capture techniques and fluorescence microscopy approaches to visualize the organization of DNA elements in the nucleus. Here we will summarize relevant aspects of how the interplay between CRE activity and chromatin 3D organization regulates developmental gene expression and how it relates to pathological conditions and the evolution of animal body plan.

The Drosophila Fab-7 boundary modulates Abd-B gene activity by guiding an inversion of collinear chromatin organization and alternate promoter use

Hox genes encode transcription factors that specify segmental identities along the anteroposterior body axis. These genes are organized in clusters, where their order corresponds to their activity along the body axis, a feature known as collinearity. In Drosophila, the BX-C cluster contains the three most posterior Hox genes, where their collinear activation incorporates progressive changes in histone modifications, chromatin architecture, and use of boundary elements and cis-regulatory regions. To dissect functional hierarchies, we compare chromatin organization in cell lines and larvae, with a focus on the Abd-B gene. Our work establishes the importance of the Fab-7 boundary for insulation between 3D domains carrying different histone modifications. Interestingly, we detect a non-canonical inversion of collinear chromatin dynamics at Abd-B, with the domain of active histone modifications progressively decreasing in size. This dynamic chromatin organization differentially activates the alternative promoters of the Abd-B gene, thereby expanding the possibilities for fine-tuning of transcriptional output.

What snakes and caecilians have in common? Molecular interaction units and the independent origins of similar morphotypes in Tetrapoda

Developmental pathways encompass transcription factors and cis-regulatory elements that interact as transcription factor-regulatory element (TF-RE) units. Independent origins of similar phenotypes likely involve changes in different parts of these units, a hypothesis promisingly tested addressing the evolution of the rib-associated lumbar (RAL) morphotype that characterizes emblematic animals such as snakes and elephants. Previous investigation in these lineages identified a polymorphism in the Homology region 1 [H1] enhancer of the Myogenic factor-5 [Myf5], which interacts with HOX10 proteins to modulate rib development. Here we address the evolution of TF-RE units focusing on independent origins of RAL morphotypes. We compiled an extensive database for H1-Myf5 and HOX10 sequences with two goals: (i) evaluate if the enhancer polymorphism is present in amphibians exhibiting the RAL morphotype and (ii) test a hypothesis of enhanced evolutionary flexibility mediated by TF-RE units, according to which independent origins of the RAL morphotype might involve changes in either component of the interaction unit. We identified the H1-Myf5 polymorphism in lineages that diverged around 340 Ma, including Lissamphibia. Independent origins of the RAL morphotype in Tetrapoda involved sequence variation in either component of the TF-RE unit, confirming that different changes may similarly affect the phenotypic outcome of a given developmental pathway.

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.

Synthetic regulatory reconstitution reveals principles of mammalian Hox cluster regulation

Precise Hox gene expression is crucial for embryonic patterning. Intra-Hox transcription factor binding and distal enhancer elements have emerged as the major regulatory modules controlling Hox gene expression. However, quantifying their relative contributions has remained elusive. Here, we introduce "synthetic regulatory reconstitution," a conceptual framework for studying gene regulation, and apply it to the HoxA cluster. We synthesized and delivered variant rat HoxA clusters (130 to 170 kilobases) to an ectopic location in the mouse genome. We found that a minimal HoxA cluster recapitulated correct patterns of chromatin remodeling and transcription in response to patterning signals, whereas the addition of distal enhancers was needed for full transcriptional output. Synthetic regulatory reconstitution could provide a generalizable strategy for deciphering the regulatory logic of gene expression in complex genomes.

Context-dependent enhancer function revealed by targeted inter-TAD relocation

The expression of some genes depends on large, adjacent regions of the genome that contain multiple enhancers. These regulatory landscapes frequently align with Topologically Associating Domains (TADs), where they integrate the function of multiple similar enhancers to produce a global, TAD-specific regulation. We asked if an individual enhancer could overcome the influence of one of these landscapes, to drive gene transcription. To test this, we transferred an enhancer from its native location, into a nearby TAD with a related yet different functional specificity. We used the biphasic regulation of Hoxd genes during limb development as a paradigm. These genes are first activated in proximal limb cells by enhancers located in one TAD, which is then silenced when the neighboring TAD activates its enhancers in distal limb cells. We transferred a distal limb enhancer into the proximal limb TAD and found that its new context suppresses its normal distal specificity, even though it is bound by HOX13 transcription factors, which are responsible for the distal activity. This activity can be rescued only when a large portion of the surrounding environment is removed. These results indicate that, at least in some cases, the functioning of enhancer elements is subordinated to the host chromatin context, which can exert a dominant control over its activity.

Here the authors show that a strong enhancer sequence can be controlled by the chromatin environment provided by a topologically associated domain (TAD) located nearby. An enhancer relocated by homologous recombination takes all the hallmarks of its new neighboring enhancers located in the recipient TAD.

Pleiotropic enhancers are ubiquitous regulatory elements in the human genome

Enhancers are regulatory elements of genomes that determine spatio-temporal patterns of gene expression. The human genome contains a vast number of enhancers, which largely outnumber protein-coding genes. Historically, enhancers have been regarded as highly tissue-specific. However, recent evidence has demonstrated that many enhancers are pleiotropic, with activity in multiple developmental contexts. Yet, the extent and impact of pleiotropy remain largely unexplored. In this study we analyzed active enhancers across human organs based on the analysis of both eRNA transcription (FANTOM5 consortium data sets) and chromatin architecture (ENCODE consortium data sets). We show that pleiotropic enhancers are pervasive in the human genome and that most enhancers active in a particular organ are also active in other organs. In addition, our analysis suggests that the proportion of context-specific enhancers of a given organ is explained, at least in part, by the proportion of context-specific genes in that same organ. The notion that such a high proportion of human enhancers can be pleiotropic suggests that small regions of regulatory DNA contain abundant regulatory information and that these regions evolve under important evolutionary constraints.

HiCAR is a robust and sensitive method to analyze open-chromatin-associated genome organization

The long-range interactions of cis-regulatory elements (cREs) play a central role in gene regulation. cREs can be characterized as accessible chromatin sequences. However, it remains technically challenging to comprehensively identify their spatial interactions. Here, we report a new method HiCAR (Hi-C on accessible regulatory DNA), which utilizes Tn5 transposase and chromatin proximity ligation, for the analysis of open-chromatin-anchored interactions with low-input cells. By applying HiCAR in human embryonic stem cells and lymphoblastoid cells, we demonstrate that HiCAR identifies high-resolution chromatin contacts with an efficiency comparable with that of in situ Hi-C over all distance ranges. Interestingly, we found that the "poised" gene promoters exhibit silencer-like function to repress the expression of distal genes via promoter-promoter interactions. Lastly, we applied HiCAR to 30,000 primary human muscle stem cells and demonstrated that HiCAR is capable of analyzing chromatin accessibility and looping using low-input primary cells and clinical samples.

Sex-specific evolution of a Drosophila sensory system via interacting cis- and trans-regulatory changes

The evolution of gene expression via cis-regulatory changes is well established as a major driver of phenotypic evolution. However, relatively little is known about the influence of enhancer architecture and intergenic interactions on regulatory evolution. We address this question by examining chemosensory system evolution in Drosophila. Drosophila prolongata males show a massively increased number of chemosensory bristles compared to females and males of sibling species. This increase is driven by sex-specific transformation of ancestrally mechanosensory organs. Consistent with this phenotype, the Pox neuro transcription factor (Poxn), which specifies chemosensory bristle identity, shows expanded expression in D. prolongata males. Poxn expression is controlled by nonadditive interactions among widely dispersed enhancers. Although some D. prolongata Poxn enhancers show increased activity, the additive component of this increase is slight, suggesting that most changes in Poxn expression are due to epistatic interactions between Poxn enhancers and trans-regulatory factors. Indeed, the expansion of D. prolongata Poxn enhancer activity is only observed in cells that express doublesex (dsx), the gene that controls sexual differentiation in Drosophila and also shows increased expression in D. prolongata males due to cis-regulatory changes. Although expanded dsx expression may contribute to increased activity of D. prolongata Poxn enhancers, this interaction is not sufficient to explain the full expansion of Poxn expression, suggesting that cis-trans interactions between Poxn, dsx, and additional unknown genes are necessary to produce the derived D. prolongata phenotype. Overall, our results demonstrate the importance of epistatic gene interactions for evolution, particularly when pivotal genes have complex regulatory architecture.

Contribution of 3D genome topological domains to genetic risk of cancers: a genome-wide computational study

Background

Genome-wide association studies have identified statistical associations between various diseases, including cancers, and a large number of single-nucleotide polymorphisms (SNPs). However, they provide no direct explanation of the mechanisms underlying the association. Based on the recent discovery that changes in three-dimensional genome organization may have functional consequences on gene regulation favoring diseases, we investigated systematically the genome-wide distribution of disease-associated SNPs with respect to a specific feature of 3D genome organization: topologically associating domains (TADs) and their borders.

Results

For each of 449 diseases, we tested whether the associated SNPs are present in TAD borders more often than observed by chance, where chance (i.e., the null model in statistical terms) corresponds to the same number of pointwise loci drawn at random either in the entire genome, or in the entire set of disease-associated SNPs listed in the GWAS catalog. Our analysis shows that a fraction of diseases displays such a preferential localization of their risk loci. Moreover, cancers are relatively more frequent among these diseases, and this predominance is generally enhanced when considering only intergenic SNPs. The structure of SNP-based diseasome networks confirms that localization of risk loci in TAD borders differs between cancers and non-cancer diseases. Furthermore, different TAD border enrichments are observed in embryonic stem cells and differentiated cells, consistent with changes in topological domains along embryogenesis and delineating their contribution to disease risk.

Conclusions

Our results suggest that, for certain diseases, part of the genetic risk lies in a local genetic variation affecting the genome partitioning in topologically insulated domains. Investigating this possible contribution to genetic risk is particularly relevant in cancers. This study thus opens a way of interpreting genome-wide association studies, by distinguishing two types of disease-associated SNPs: one with an effect on an individual gene, the other acting in interplay with 3D genome organization.

Supplementary Information

The online version contains supplementary material available at 10.1186/s40246-022-00375-2.

Transcriptional enhancers and their communication with gene promoters

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

Analyses of Promoter , Enhancer, and Nucleosome Organization in Mammalian Cells by MNase-Seq

MNase-Seq is a genome-wide procedure that allows mapping of DNA associated to nucleosomes following micrococcal nuclease digestion. It is a rapid and robust technology useful for the analysis of chromatin properties genome-wide at the resolution of mono-nucleosomes. Here, we describe how to produce high-resolution nucleosome maps of cells grown in suspension or adherent mammalian cells. After only three steps: nuclei or cell preparation, native MNase digestion and DNA purification, libraries for high-throughput sequencing can be prepared. Genome-wide nucleosome maps allow analyzing chromatin opening at promoters or enhancers, nucleosome displacement, or labile nucleosome occupancy depending on the digestion condition used. As presented, MNase-Seq is a versatile tool for investigating chromatin dynamics, regulation, and to define open chromatin regions of regulatory elements in mammalian genomes.

Chromatin Conformation in Development and Disease

Chromatin domains and loops are important elements of chromatin structure and dynamics, but much remains to be learned about their exact biological role and nature. Topological associated domains and functional loops are key to gene expression and hold the answer to many questions regarding developmental decisions and diseases. Here, we discuss new findings, which have linked chromatin conformation with development, differentiation and diseases and hypothesized on various models while integrating all recent findings on how chromatin architecture affects gene expression during development, evolution and disease.

Topologically Associating Domains and Regulatory Landscapes in Development, Evolution and Disease

Animal genomes are folded in topologically associating domains (TADs) that have been linked to the regulation of the genes they contain by constraining regulatory interactions between cis-regulatory elements and promoters. Therefore, TADs are proposed as structural scaffolds for the establishment of regulatory landscapes (RLs). In this review, we discuss recent advances in the connection between TADs and gene regulation, their relationship with gene RLs and their dynamics during development and differentiation. Moreover, we describe how restructuring TADs may lead to pathological conditions, which explains their high evolutionary conservation, but at the same time it provides a substrate for the emergence of evolutionary innovations that lay at the origin of vertebrates and other phylogenetic clades.

Dbx2 regulation in limbs suggests interTAD sharing of enhancers

BACKGROUND

During tetrapod limb development, the HOXA13 and HOXD13 transcription factors are critical for the emergence and organization of the autopod, the most distal aspect where digits will develop. Since previous work had suggested that the Dbx2 gene is a target of these factors, we set up to analyze in detail this potential regulatory interaction.

RESULTS

We show that HOX13 proteins bind to mammalian-specific sequences at the vicinity of the Dbx2 locus that have enhancer activity in developing digits. However, the functional inactivation of the DBX2 protein did not elicit any particular phenotype related to Hox genes inactivation in digits, suggesting either redundant or compensatory mechanisms. We report that the neighboring Nell2 and Ano6 genes are also expressed in distal limb buds and are in part controlled by the same Dbx2 enhancers despite being localized into two different topologically associating domains (TADs) flanking the Dbx2 locus.

CONCLUSIONS

We conclude that Hoxa13 and Hoxd genes cooperatively activate Dbx2 expression in developing digits through binding to mammalian specific regulatory sequences in the Dbx2 neighborhood. Furthermore, these enhancers can overcome TAD boundaries in either direction to co-regulate a set of genes located in distinct chromatin domains.

A developmental perspective of homology and evolutionary novelty

The development and evolution of multicellular body plans is complex. Many distinct organs and body parts must be reproduced at each generation, and those that are traceable over long time scales are considered homologous. Among the most pressing and least understood phenomena in evolutionary biology is the mode by which new homologs, or "novelties" are introduced to the body plan and whether the developmental changes associated with such evolution deserve special treatment. In this chapter, we address the concepts of homology and evolutionary novelty through the lens of development. We present a series of case studies, within insects and vertebrates, from which we propose a developmental model of multicellular organ identity. With this model in hand, we make predictions regarding the developmental evolution of body plans and highlight the need for more integrative analysis of developing systems.

Chromatin topology and the timing of enhancer function at the HoxD locus

The HoxD gene cluster is critical for proper limb formation in tetrapods. In the emerging limb buds, different subgroups of Hoxd genes respond first to a proximal regulatory signal, then to a distal signal that organizes digits. These two regulations are exclusive from one another and emanate from two distinct topologically associating domains (TADs) flanking HoxD, both containing a range of appropriate enhancer sequences. The telomeric TAD (T-DOM) contains several enhancers active in presumptive forearm cells and is divided into two sub-TADs separated by a CTCF-rich boundary, which defines two regulatory submodules. To understand the importance of this particular regulatory topology to control Hoxd gene transcription in time and space, we either deleted or inverted this sub-TAD boundary, eliminated the CTCF binding sites, or inverted the entire T-DOM to exchange the respective positions of the two sub-TADs. The effects of such perturbations on the transcriptional regulation of Hoxd genes illustrate the requirement of this regulatory topology for the precise timing of gene activation. However, the spatial distribution of transcripts was eventually resumed, showing that the presence of enhancer sequences, rather than either their exact topology or a particular chromatin architecture, is the key factor. We also show that the affinity of enhancers to find their natural target genes can overcome the presence of both a strong TAD border and an unfavorable orientation of CTCF sites.

Transcriptional Trajectories in Mouse Limb Buds Reveal the Transition from Anterior-Posterior to Proximal-Distal Patterning at Early Limb Bud Stage

Limb patterning relies in large part on the function of the Hox family of developmental genes. While the differential expression of Hox genes shifts from the anterior-posterior (A-P) to the proximal-distal (P-D) axis around embryonic day 11 (E11), whether this shift coincides with a more global change of A-P to P-D patterning program remains unclear. By performing and analyzing the transcriptome of the developing limb bud from E10.5 to E12.5, at single-cell resolution, we have uncovered transcriptional trajectories that revealed a general switch from A-P to P-D genetic program between E10.5 and E11.5. Interestingly, all the transcriptional trajectories at E10.5 end with cells expressing either proximal or distal markers suggesting a progressive acquisition of P-D identity. Moreover, we identified three categories of genes expressed in the distal limb mesenchyme characterized by distinct temporal expression dynamics. Among these are Hoxa13 and Hoxd13 (Hox13 hereafter), which start to be expressed around E10.5, and importantly the binding of the HOX13 factors was observed within or in the neighborhood of several of the distal limb genes. Our data are consistent with previous evidence suggesting that the transition from the early/proximal to the late/distal transcriptome of the limb mesenchyme largely relies on HOX13 function. Based on these results and the evidence that HOX13 factors restrict Hoxa11 expression to the proximal limb, in progenitor cells of the zeugopod, we propose that HOX13 act as a key determinant of P-D patterning.

Lysine demethylase 7a regulates murine anterior-posterior development by modulating the transcription of Hox gene cluster

Temporal and spatial colinear expression of the Hox genes determines the specification of positional identities during vertebrate development. Post-translational modifications of histones contribute to transcriptional regulation. Lysine demethylase 7A (Kdm7a) demethylates lysine 9 or 27 di-methylation of histone H3 (H3K9me2, H3K27me2) and participates in the transcriptional activation of developmental genes. However, the role of Kdm7a during mouse embryonic development remains to be elucidated. Herein, we show that Kdm7a^−/− mouse exhibits an anterior homeotic transformation of the axial skeleton, including an increased number of presacral elements. Importantly, posterior Hox genes (caudally from Hox9) are specifically downregulated in the Kdm7a^−/− embryo, which correlates with increased levels of H3K9me2, not H3K27me2. These observations suggest that Kdm7a controls the transcription of posterior Hox genes, likely via its demethylating activity, and thereby regulating the murine anterior-posterior development. Such epigenetic regulatory mechanisms may be harnessed for proper control of coordinate body patterning in vertebrates.

Higashijima et al show that mice lacking the Kdm7a demethylase exhibits anterior homeotic transformation of the axial skeleton and downregulation of posterior Hox gene transcription and these changes are associated with increased H3K9me2 at posterior Hox loci. These findings provide insights into the epigenetic control of Hox-mediated patterning in embryogenesis.

Three distinct mechanisms, Notch instructive, permissive, and independent, regulate the expression of two different pericardial genes to specify cardiac cell subtypes

The development of a complex organ involves the specification and differentiation of diverse cell types constituting that organ. Two major cell subtypes, contractile cardial cells (CCs) and nephrocytic pericardial cells (PCs), comprise the Drosophila heart. Binding sites for Suppressor of Hairless [Su(H)], an integral transcription factor in the Notch signaling pathway, are enriched in the enhancers of PC-specific genes. Here we show three distinct mechanisms regulating the expression of two different PC-specific genes, Holes in muscle (Him), and Zn finger homeodomain 1 (zfh1). Him transcription is activated in PCs in a permissive manner by Notch signaling: in the absence of Notch signaling, Su(H) forms a repressor complex with co-repressors and binds to the Him enhancer, repressing its transcription; upon alleviation of this repression by Notch signaling, Him transcription is activated. In contrast, zfh1 is transcribed by a Notch-instructive mechanism in most PCs, where mere alleviation of repression by preventing the binding of Su(H)-co-repressor complex is not sufficient to activate transcription. Our results suggest that upon activation of Notch signaling, the Notch intracellular domain associates with Su(H) to form an activator complex that binds to the zfh1 enhancer, and that this activator complex is necessary for bringing about zfh1 transcription in these PCs. Finally, a third, Notch-independent mechanism activates zfh1 transcription in the remaining, even skipped-expressing, PCs. Collectively, our data show how the same feature, enrichment of Su(H) binding sites in PC-specific gene enhancers, is utilized by two very distinct mechanisms, one permissive, the other instructive, to contribute to the same overall goal: the specification and differentiation of a cardiac cell subtype by activation of the pericardial gene program. Furthermore, our results demonstrate that the zfh1 enhancer drives expression in two different domains using distinct Notch-instructive and Notch-independent mechanisms.

Exonuclease combinations reduce noises in 3D genomics technologies

Chromosome conformation-capture technologies are widely used in 3D genomics; however, experimentally, such methods have high-noise limitations and, therefore, require significant bioinformatics efforts to extract reliable distal interactions. Miscellaneous undesired linear DNAs, present during proximity-ligation, represent a main noise source, which needs to be minimized or eliminated. In this study, different exonuclease combinations were tested to remove linear DNA fragments from a circularized DNA preparation. This method efficiently removed linear DNAs, raised the proportion of annulation and increased the valid-pairs ratio from ∼40% to ∼80% for enhanced interaction detection in standard Hi-C. This strategy is applicable for development of various 3D genomics technologies, or optimization of Hi-C sequencing efficiency.

Enhancer and super-enhancer dynamics in repair after ischemic acute kidney injury

The endogenous repair process can result in recovery after acute kidney injury (AKI) with adaptive proliferation of tubular epithelial cells, but repair can also lead to fibrosis and progressive kidney disease. There is currently limited knowledge about transcriptional regulators regulating these repair programs. Herein we establish the enhancer and super-enhancer landscape after AKI by ChIP-seq in uninjured and repairing kidneys on day two after ischemia reperfusion injury (IRI). We identify key transcription factors including HNF4A, GR, STAT3 and STAT5, which show specific binding at enhancer and super-enhancer sites, revealing enhancer dynamics and transcriptional changes during kidney repair. Loss of bromodomain-containing protein 4 function before IRI leads to impaired recovery after AKI and increased mortality. Our comprehensive analysis of epigenetic changes after kidney injury in vivo has the potential to identify targets for therapeutic intervention. Importantly, our data also call attention to potential caveats involved in use of BET inhibitors in patients at risk for AKI.

Developmental constraints on fin diversity

The fish fin is a breathtaking repository full of evolutionary diversity, novelty, and convergence. Over 500 million years, the adaptation to novel habitats has provided landscapes of fin diversity. Although comparative anatomy of evolutionarily divergent patterns over centuries has highlighted the fundamental architectures and evolutionary trends of fins, including convergent evolution, the developmental constraints on fin evolution, which bias the evolutionary trajectories of fin morphology, largely remain elusive. Here, we review the evolutionary history, developmental mechanisms, and evolutionary underpinnings of paired fins, illuminating possible developmental constraints on fin evolution. Our compilation of anatomical and genetic knowledge of fin development sheds light on the canalized and the unpredictable aspects of fin shape in evolution. Leveraged by an arsenal of genomic and genetic tools within the working arena of spectacular fin diversity, evolutionary developmental biology embarks on the establishment of conceptual framework for developmental constraints, previously enigmatic properties of evolution.

A complex regulatory landscape involved in the development of mammalian external genitals

Developmental genes are often controlled by large regulatory landscapes matching topologically associating domains (TADs). In various contexts, the associated chromatin backbone is modified by specific enhancer-enhancer and enhancer-promoter interactions. We used a TAD flanking the mouse HoxD cluster to study how these regulatory architectures are formed and deconstructed once their function achieved. We describe this TAD as a functional unit, with several regulatory sequences acting together to elicit a transcriptional response. With one exception, deletion of these sequences didn't modify the transcriptional outcome, a result at odds with a conventional view of enhancer function. The deletion and inversion of a CTCF site located near these regulatory sequences did not affect transcription of the target gene. Slight modifications were nevertheless observed, in agreement with the loop extrusion model. We discuss these unexpected results considering both conventional and alternative explanations relying on the accumulation of poorly specific factors within the TAD backbone.

Why a Constant Number of Vertebrae? Digital Control of Segmental Identity during Vertebrate Development: The Somite Cycle Controls a Digital, Chromatin-Based Counter That Defines Segmental Identity and Body Plans in Vertebrate Animals

It is not understood how the numbers and identities of vertebrae are controlled during mammalian development. The remarkable robustness and conservation of segmental numbers may suggest the digital nature of the underlying process. The study proposes a mechanism that allows cells to obtain and store the segmental information in digital form, and to produce a pattern of chromatin accessibility that in turn regulates Hox gene expression specific to the metameric segment. The model requires that a regulatory element be present such that the number of occurrences of the motif between two consecutive Hox genes equals the number of segments under the control of the anterior gene. This is true for the recently discovered hydroxyl radical cleavage 3bp-periodic (HRC3) motif, associated with histone modifications and developmental genes. The finding not only allows the correct prediction of the numbers of segments using only sequence information, but also resolves the 40-year-old enigma of the function of temporal and spatial collinearity of Hox genes. The logic of the mechanism is illustrated in the attached animated video. How different aspects of the proposed mechanism can be tested experimentally is also discussed.

Ultraconserved Elements Occupy Specific Arenas of Three-Dimensional Mammalian Genome Organization

This study explores the relationship between three-dimensional genome organization and ultraconserved elements (UCEs), an enigmatic set of DNA elements that are perfectly conserved between the reference genomes of distantly related species. Examining both human and mouse genomes, we interrogate the relationship of UCEs to three features of chromosome organization derived from Hi-C studies. We find that UCEs are enriched within contact domains and, further, that the subset of UCEs within domains shared across diverse cell types are linked to kidney-related and neuronal processes. In boundaries, UCEs are generally depleted, with those that do overlap boundaries being overrepresented in exonic UCEs. Regarding loop anchors, UCEs are neither overrepresented nor underrepresented, but those present in loop anchors are enriched for splice sites. Finally, as the relationships between UCEs and human Hi-C features are conserved in mouse, our findings suggest that UCEs contribute to interspecies conservation of genome organization and, thus, genome stability.

The Role of Insulation in Patterning Gene Expression

Development is orchestrated by regulatory elements that turn genes ON or OFF in precise spatial and temporal patterns. Many safety mechanisms prevent inappropriate action of a regulatory element on the wrong gene promoter. In flies and mammals, dedicated DNA elements (insulators) recruit protein factors (insulator binding proteins, or IBPs) to shield promoters from regulatory elements. In mammals, a single IBP called CCCTC-binding factor (CTCF) is known, whereas genetic and biochemical analyses in Drosophila have identified a larger repertoire of IBPs. How insulators function at the molecular level is not fully understood, but it is currently thought that they fold chromosomes into conformations that affect regulatory element-promoter communication. Here, we review the discovery of insulators and describe their properties. We discuss recent genetic studies in flies and mice to address the question: Is gene insulation important for animal development? Comparing and contrasting observations in these two species reveal that they have different requirements for insulation, but that insulation is a conserved and critical gene regulation strategy.

The living organism: evolutionary design or an accident

The presented work discusses some evolutionary phenomena underlining the complexity of organism creation and surprisingly the short evolutionary time of this process in particular. Uncommonness of this process ensued from the necessary simultaneous combining of highly complicated biological mechanisms, of which some were generated independently before the direct evolutionary demand. This in conclusion points to still not fully understood biological program ensuring superiority of the permanent evolutionary progress over effects of purely random mutational changes as the driving mechanism in evolution.

Preformed chromatin topology assists transcriptional robustness of Shh during limb development

Long-range gene regulation involves physical proximity between enhancers and promoters to generate precise patterns of gene expression in space and time. However, in some cases, proximity coincides with gene activation, whereas, in others, preformed topologies already exist before activation. In this study, we investigate the preformed configuration underlying the regulation of the Shh gene by its unique limb enhancer, the ZRS, in vivo during mouse development. Abrogating the constitutive transcription covering the ZRS region led to a shift within the Shh-ZRS contacts and a moderate reduction in Shh transcription. Deletion of the CTCF binding sites around the ZRS resulted in the loss of the Shh-ZRS preformed interaction and a 50% decrease in Shh expression but no phenotype, suggesting an additional, CTCF-independent mechanism of promoter-enhancer communication. This residual activity, however, was diminished by combining the loss of CTCF binding with a hypomorphic ZRS allele, resulting in severe Shh loss of function and digit agenesis. Our results indicate that the preformed chromatin structure of the Shh locus is sustained by multiple components and acts to reinforce enhancer-promoter communication for robust transcription.

Inference of the ancestral vertebrate phenotype through vestiges of the whole-genome duplications

Inferring the phenotype of the last common ancestor of living vertebrates is a challenging problem because of several unresolvable factors. They include the lack of reliable out-groups of living vertebrates, poor information about less fossilizable organs and specialized traits of phylogenetically important species, such as lampreys and hagfishes (e.g. secondary loss of vertebrae in adult hagfishes). These factors undermine the reliability of ancestral reconstruction by traditional character mapping approaches based on maximum parsimony. In this article, we formulate an approach to hypothesizing ancestral vertebrate phenotypes using information from the phylogenetic and functional properties of genes duplicated by genome expansions in early vertebrate evolution. We named the conjecture as ‘chronological reconstruction of ohnolog functions (CHROF)’. This CHROF conjecture raises the possibility that the last common ancestor of living vertebrates may have had more complex traits than currently thought.

'Building a perfect body': control of vertebrate organogenesis by PBX-dependent regulatory networks

Pbx genes encode transcription factors that belong to the TALE (three-amino-acid loop extension) superclass of homeodomain proteins. We have witnessed a surge in information about the roles of this gene family as leading actors in the transcriptional control of development. PBX proteins represent a clear example of how transcription factors can regulate developmental processes by combinatorial properties, acting within multimeric complexes to implement activation or repression of transcription depending on their interaction partners. Here, we revisit long-emphasized functions of PBX transcription factors as cofactors for HOX proteins, major architects of the body plan. We further discuss new knowledge on roles of PBX proteins in different developmental contexts as upstream regulators of Hox genes-as factors that interact with non-HOX proteins and can work independently of HOX-as well as potential pioneer factors. Committed to building a perfect body, PBX proteins govern regulatory networks that direct essential morphogenetic processes and organogenesis in vertebrate development. Perturbations of PBX-dependent networks can cause human congenital disease and cancer.

A spontaneous mouse deletion in Mctp1 uncovers a long-range cis-regulatory region crucial for NR2F1 function during inner ear development

Hundreds of thousands of cis-regulatory DNA sequences are predicted in vertebrate genomes, but unlike genes themselves, few have been characterized at the functional level or even unambiguously paired with a target gene. Here we serendipitously identified and started investigating the first reported long-range regulatory region for the Nr2f1 (Coup-TFI) transcription factor gene. NR2F1 is temporally and spatially regulated during development and required for patterning and regionalization in the nervous system, including sensory hair cell organization in the auditory epithelium of the cochlea. Analyzing the deaf wanderer (dwnd) spontaneous mouse mutation, we traced back the cause of its associated circling behavior to a 53 kb deletion removing five exons and adjacent intronic regions of the poorly characterized Mctp1 gene. Interestingly, loss of Mctp1 function cannot account for the hearing loss, inner ear dysmorphology and sensory hair cell disorganization observed in dwnd mutants. Instead, we found that the Mctp1dwnd deletion affects the Nr2f1 gene located 1.4 Mb away, downregulating transcription and protein expression in the embryonic cochlea. Remarkably, the Mctp1dwnd allele failed to complement a targeted inactivation allele of Nr2f1, and transheterozygotes or Mctp1dwnd homozygotes exhibit the same morphological defects observed in inner ears of Nr2f1 mutants without sharing their early life lethality. Defects include improper separation of the utricle and saccule in the vestibule not described previously, which can explain the circling behavior that first brought the spontaneous mutation to attention. By contrast, mice homozygous for a targeted inactivation of Mctp1 have normal hearing and inner ear structures. We conclude that the 53 kb Mctp1dwnd deletion encompasses a long-range cis-regulatory region essential for proper Nr2f1 expression in the embryonic inner ear, providing a first opportunity to investigate Nr2f1 function in postnatal inner ears. This work adds to the short list of long-range regulatory regions characterized as essential to drive expression of key developmental control genes.

Correlated Evolution of Two Copulatory Organs via a Single cis-Regulatory Nucleotide Change

Diverse traits often covary between species [1-3]. The possibility that a single mutation could contribute to the evolution of several characters between species [3] is rarely investigated as relatively few cases are dissected at the nucleotide level. Drosophila santomea has evolved additional sex comb sensory teeth on its legs and has lost two sensory bristles on its genitalia. We present evidence that a single nucleotide substitution in an enhancer of the scute gene contributes to both changes. The mutation alters a binding site for the Hox protein Abdominal-B in the developing genitalia, leading to bristle loss, and for another factor in the developing leg, leading to bristle gain. Our study suggests that morphological evolution between species can occur through a single nucleotide change affecting several sexually dimorphic traits. VIDEO ABSTRACT.

The role of 3D genome organization in disease: From compartments to single nucleotides

Since the advent of the chromosome conformation capture technology, our understanding of the human genome 3D organization has grown rapidly and we now know that human interphase chromosomes are folded into multiple layers of hierarchical structures and each layer can play a critical role in transcriptional regulation. Alterations in any one of these finely-tuned layers can lead to unwanted cascade of molecular events and ultimately drive the manifestation of diseases and phenotypes. Here we discuss, starting from chromosome level organization going down to single nucleotide changes, recent studies linking diseases or phenotypes to changes in the 3D genome architecture.

An evolutionary perspective of regulatory landscape dynamics in development and disease

The organization of animal genomes into topologically associating domains (TADs) provides a structural scaffold in which cis-regulatory elements (CREs) operate on their target genes. Determining the position of CREs and genes relative to TADs has become instrumental to trace gene expression changes during evolution and in diseases. Here we will review recent studies and discuss TADs as structural units with respect to their conservation and stability during genome reorganization. Furthermore, we describe how TAD restructuring contributed to morphological novelties during evolution but also their deleterious effects associated with disease. Despite considering TADs as structural units, the nested and dynamic scaffold within TADs contributes to tissue-specific gene expression, implying that such changes can also account for gene expression differences during evolution.

Enhancer Logic and Mechanics in Development and Disease

Enhancers are distally located genomic cis-regulatory elements that integrate spatiotemporal cues to coordinate gene expression in a tissue-specific manner during metazoan development. Enhancer function depends on a combination of bound transcription factors and cofactors that regulate local chromatin structure, as well as on the topological interactions that are necessary for their activity. Numerous genome-wide studies concur that the vast majority of disease-associated variations occur within non-coding genomic sequences, in other words the 'cis-regulome', and this underscores their relevance for human health. Advances in DNA sequencing and genome-editing technologies have dramatically expanded our ability to identify enhancers and investigate their properties in vivo, revealing an extraordinary level of interconnectivity underlying cis-regulatory networks. We discuss here these recently developed methodologies, as well as emerging trends and remaining questions in the field of enhancer biology, and how perturbation of enhancer activities/functions results in enhanceropathies.

Worldwide distribution of the DCDC2 READ1 regulatory element and its relationship with phoneme variation across languages

Languages evolve rapidly due to an interaction between sociocultural factors and underlying phonological processes that are influenced by genetic factors. DCDC2 has been strongly associated with core components of the phonological processing system in animal models and multiple independent studies of populations and languages. To characterize subtle language differences arising from genetic variants associated with phonological processes, we examined the relationship between READ1, a regulatory element in DCDC2, and phonemes in languages of 43 populations across five continents. Variation in READ1 was significantly correlated with the number of consonants. Our results suggest that subtle cognitive biases conferred by different READ1 alleles are amplified through cultural transmission that shape consonant use by populations over time.

DCDC2 is a gene strongly associated with components of the phonological processing system in animal models and in multiple independent studies of populations and languages. We propose that it may also influence population-level variation in language component usage. To test this hypothesis, we investigated the evolution and worldwide distribution of the READ1 regulatory element within DCDC2, and compared its distribution with variation in different language properties. The mutational history of READ1 was estimated by examining primate and archaic hominin sequences. This identified duplication and expansion events, which created a large number of polymorphic alleles based on internal repeat units (RU1 and RU2). Association of READ1 alleles was studied with respect to the numbers of consonants and vowels for languages in 43 human populations distributed across five continents. Using population-based approaches with multivariate ANCOVA and linear mixed effects analyses, we found that the RU1-1 allele group of READ1 is significantly associated with the number of consonants within languages independent of genetic relatedness, geographic proximity, and language family. We propose that allelic variation in READ1 helped create a subtle cognitive bias that was amplified by cultural transmission, and ultimately shaped consonant use by different populations over time.

Genomic positional conservation identifies topological anchor point RNAs linked to developmental loci

BACKGROUND

The mammalian genome is transcribed into large numbers of long noncoding RNAs (lncRNAs), but the definition of functional lncRNA groups has proven difficult, partly due to their low sequence conservation and lack of identified shared properties. Here we consider promoter conservation and positional conservation as indicators of functional commonality.

RESULTS

We identify 665 conserved lncRNA promoters in mouse and human that are preserved in genomic position relative to orthologous coding genes. These positionally conserved lncRNA genes are primarily associated with developmental transcription factor loci with which they are coexpressed in a tissue-specific manner. Over half of positionally conserved RNAs in this set are linked to chromatin organization structures, overlapping binding sites for the CTCF chromatin organiser and located at chromatin loop anchor points and borders of topologically associating domains (TADs). We define these RNAs as topological anchor point RNAs (tapRNAs). Characterization of these noncoding RNAs and their associated coding genes shows that they are functionally connected: they regulate each other's expression and influence the metastatic phenotype of cancer cells in vitro in a similar fashion. Furthermore, we find that tapRNAs contain conserved sequence domains that are enriched in motifs for zinc finger domain-containing RNA-binding proteins and transcription factors, whose binding sites are found mutated in cancers.

CONCLUSIONS

This work leverages positional conservation to identify lncRNAs with potential importance in genome organization, development and disease. The evidence that many developmental transcription factors are physically and functionally connected to lncRNAs represents an exciting stepping-stone to further our understanding of genome regulation.

Embryonic timing, axial stem cells, chromatin dynamics, and the Hox clock

Collinear regulation of Hox genes in space and time has been an outstanding question ever since the initial work of Ed Lewis in 1978. Here we discuss recent advances in our understanding of this phenomenon in relation to novel concepts associated with large-scale regulation and chromatin structure during the development of both axial and limb patterns. We further discuss how this sequential transcriptional activation marks embryonic stem cell-like axial progenitors in mammals and, consequently, how a temporal genetic system is further translated into spatial coordinates via the fate of these progenitors. In this context, we argue the benefit and necessity of implementing this unique mechanism as well as the difficulty in evolving an alternative strategy to deliver this critical positional information.

The HoxD cluster is a dynamic and resilient TAD boundary controlling the segregation of antagonistic regulatory landscapes

The mammalian HoxD cluster lies between two topologically associating domains (TADs) matching distinct enhancer-rich regulatory landscapes. During limb development, the telomeric TAD controls the early transcription of Hoxd genes in forearm cells, whereas the centromeric TAD subsequently regulates more posterior Hoxd genes in digit cells. Therefore, the TAD boundary prevents the terminal Hoxd13 gene from responding to forearm enhancers, thereby allowing proper limb patterning. To assess the nature and function of this CTCF-rich DNA region in embryos, we compared chromatin interaction profiles between proximal and distal limb bud cells isolated from mutant stocks where various parts of this boundary region were removed. The resulting progressive release in boundary effect triggered inter-TAD contacts, favored by the activity of the newly accessed enhancers. However, the boundary was highly resilient, and only a 400-kb deletion, including the whole-gene cluster, was eventually able to merge the neighboring TADs into a single structure. In this unified TAD, both proximal and distal limb enhancers nevertheless continued to work independently over a targeted transgenic reporter construct. We propose that the whole HoxD cluster is a dynamic TAD border and that the exact boundary position varies depending on both the transcriptional status and the developmental context.

The three-dimensional genome: regulating gene expression during pluripotency and development

The precise expression of genes in time and space during embryogenesis is largely influenced by communication between enhancers and promoters, which is propagated and governed by the physical proximity of these elements in the nucleus. Here, we review how chromatin domains organize the genome by guiding enhancers to their target genes thereby preventing non-specific interactions with other neighboring regions. We also discuss the dynamics of chromatin interactions between enhancers and promoters, as well as the consequent changes in gene expression, that occur in pluripotent cells and during development. Finally, we evaluate how genomic changes such as deletions, inversions and duplications affect 3D chromatin configuration overall and lead to ectopic enhancer-promoter contacts, and thus gene misexpression, which can contribute to abnormal development and disease.

Appendages and gene regulatory networks: Lessons from the limbless

Among squamate reptiles, dozens of lineages have independently evolved complete or partial limb reduction. This remarkable convergence of limbless and limb-reduced phenotypes provides multiple natural replicates of different ages to explore the evolution and development of the vertebrate limb and the gene regulatory network that controls its formation. The most successful and best known of the limb-reduced squamates are snakes, which evolved a limb-reduced body form more than 100 million years ago. Recent studies have revealed the unexpected finding that many ancient limb enhancers are conserved in the genomes of snakes. Analyses in limbed animals show that many of these limb enhancers are also active during development of the phallus, suggesting that these enhancers may have been retained in snakes due their importance in regulating transcription in the external genitalia. This hypothesis is substantiated by functional tests of snake enhancers, which demonstrate that snake enhancer elements have lost limb function while retaining genital enhancer function. The large degree of overlap in the gene regulatory networks deployed during limb and phallus development may act to constrain the divergence of shared gene network components and the evolution of appendage morphology. Future studies will reveal whether limb regulatory elements have undergone similar functional changes in other lineages of limb-reduced squamates.