Similarities and differences in the regulation of HoxD genes during chick and mouse limb development

Leveraging chicken embryos for studying human enhancers

The dynamic activity of complex gene regulatory networks stands at the core of all cellular functions that define cell identity and behaviour. Gene regulatory networks comprise transcriptional enhancers, acted upon by cell-specific transcription factors to control gene expression in a spatial and temporal specific manner. Enhancers are found in the non-coding genome; pathogenic variants can disrupt enhancer activity and lead to disease. Correlating non-coding variants with aberrant enhancer activity remains a significant challenge. Due to their clinical significance, there is a longstanding interest in understanding enhancer function during early embryogenesis. With the onset of the omics era, it is now feasible to identify putative tissue-specific enhancers from epigenome data. However, such predictions in vivo require validation. The early stages of chick embryogenesis closely parallel those of human, offering an accessible in vivo model in which to assess the activity of putative human enhancer sequences. This review explores the unique advantages and recent advancements in employing chicken embryos to elucidate the activity of human transcriptional enhancers and the potential implications of these findings in human disease.

Splicing is dynamically regulated during limb development

Modifications to highly conserved developmental gene regulatory networks are thought to underlie morphological diversification in evolution and contribute to human congenital malformations. Relationships between gene expression and morphology have been extensively investigated in the limb, where most of the evidence for alterations to gene regulation in development consists of pre-transcriptional mechanisms that affect expression levels, such as epigenetic alterations to regulatory sequences and changes to cis-regulatory elements. Here we report evidence that alternative splicing (AS), a post-transcriptional process that modifies and diversifies mRNA transcripts, is dynamic during limb development in two mammalian species. We evaluated AS patterns in mouse (Mus musculus) and opossum (Monodelphis domestica) across the three key limb developmental stages: the ridge, bud, and paddle. Our data show that splicing patterns are dynamic over developmental time and suggest differences between the two mammalian taxa. Additionally, multiple key limb development genes, including Fgf8, are differentially spliced across the three stages in both species, with expression levels of the conserved splice variants, Fgf8a and Fgf8b, changing across developmental time. Our data demonstrates that AS is a critical mediator of mRNA diversity in limb development and provides an additional mechanism for evolutionary tweaking of gene dosage.

Expression Proteomics and Histone Analysis Reveal Extensive Chromatin Network Changes and a Role for Histone Tail Trimming during Cellular Differentiation

In order to understand the coordinated proteome changes associated with differentiation of a cultured cell pluripotency model, protein expression changes induced by treatment of NT2 embryonal carcinoma cells with retinoic acid were monitored by mass spectrometry. The relative levels of over 5000 proteins were mapped across distinct cell fractions. Analysis of the chromatin fraction revealed major abundance changes among chromatin proteins and epigenetic pathways between the pluripotent and differentiated states. Protein complexes associated with epigenetic regulation of gene expression, chromatin remodelling (e.g., SWI/SNF, NuRD) and histone-modifying enzymes (e.g., Polycomb, MLL) were found to be extensively regulated. We therefore investigated histone modifications before and after differentiation, observing changes in the global levels of lysine acetylation and methylation across the four canonical histone protein families, as well as among variant histones. We identified the set of proteins with affinity to peptides housing the histone marks H3K4me3 and H3K27me3, and found increased levels of chromatin-associated histone H3 tail trimming following differentiation that correlated with increased expression levels of cathepsin proteases. We further found that inhibition of cathepsins B and D reduces histone H3 clipping. Overall, the work reveals a global reorganization of the cell proteome congruent with differentiation, highlighting the key role of multiple epigenetic pathways, and demonstrating a direct link between cathepsin B and D activity and histone modification.

Imputation of 3D genome structure by genetic-epigenetic interaction modeling in mice

Gene expression is known to be affected by interactions between local genetic variation and DNA accessibility, with the latter organized into three-dimensional chromatin structures. Analyses of these interactions have previously been limited, obscuring their regulatory context, and the extent to which they occur throughout the genome. Here, we undertake a genome-scale analysis of these interactions in a genetically diverse population to systematically identify global genetic-epigenetic interaction, and reveal constraints imposed by chromatin structure. We establish the extent and structure of genotype-by-epigenotype interaction using embryonic stem cells derived from Diversity Outbred mice. This mouse population segregates millions of variants from eight inbred founders, enabling precision genetic mapping with extensive genotypic and phenotypic diversity. With 176 samples profiled for genotype, gene expression, and open chromatin, we used regression modeling to infer genetic-epigenetic interactions on a genome-wide scale. Our results demonstrate that statistical interactions between genetic variants and chromatin accessibility are common throughout the genome. We found that these interactions occur within the local area of the affected gene, and that this locality corresponds to topologically associated domains (TADs). The likelihood of interaction was most strongly defined by the three-dimensional (3D) domain structure rather than linear DNA sequence. We show that stable 3D genome structure is an effective tool to guide searches for regulatory elements and, conversely, that regulatory elements in genetically diverse populations provide a means to infer 3D genome structure. We confirmed this finding with CTCF ChIP-seq that revealed strain-specific binding in the inbred founder mice. In stem cells, open chromatin participating in the most significant regression models demonstrated an enrichment for developmental genes and the TAD-forming CTCF-binding complex, providing an opportunity for statistical inference of shifting TAD boundaries operating during early development. These findings provide evidence that genetic and epigenetic factors operate within the context of 3D chromatin structure.

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.

Transcriptome-based comparison reveals key genes regulating allometry growth of forelimb and hindlimb bone in duck embryos

Allometric growth of the forelimb and hindlimb is a widespread phenomenon observed in vertebrates. As a typical precocial bird, ducks exhibit more advanced development of their hindlimbs compared to their forelimbs, enabling them to walk shortly after hatching. This phenomenon is closely associated with the development of long bones in the embryonic stage. However, the molecular mechanism governing the allometric growth of duck forelimb and hindlimb bones is remains elusive. In this study, we employed phenotypic, histological, and gene expression analyses to investigate developmental differences between the humerus (forelimb bone) and tibia/femur (hindlimb bones) in duck embryos. Our results revealed a gradual increase in weight and length disparity between the tibia and humerus from E12 to E28 (embryo age). At E12, endochondral ossification was observed solely in the tibia but not in the humerus. The number of differentially expressed genes (DEGs) gradually increased at H12 vs. T12, H20 vs. T20, and H28 vs. T28 stages consistent with phenotypic variations. A total of 38 DEGs were found across all 3 stages. Protein-protein interaction network analysis demonstrated strong interactions among members of HOXD gene family (HOXD3/8/9/10/11/12), HOXB gene family (HOXB8/9), TBX gene family (TBX4/5/20), HOXA11, SHOX2, and MEIS2. Gene expression profiling indicated higher expression levels for all HOXD genes in the humerus compared to tibia while opposite trends were observed for HOXA/HOXB genes with low or no expression detected in the humerus. These findings suggest distinct roles played by different clusters within HOX gene family during skeletal development regulation of duck embryo's forelimbs versus hind limbs. Notably, TBX4 exhibited high expression levels specifically in tibia whereas TBX5 showed similar patterns exclusively within humerus as seen previously across other species' studies. In summary, this study identified key regulatory genes involved in allometric growth of duck forelimb and hindlimb bones during embryonic development. Skeletal development is a complex physiological process, and further research is needed to elucidate the regulatory role of candidate genes in endochondral ossification.

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.

Hnrnpk is essential for embryonic limb bud development as a transcription activator and a collaborator of insulator protein Ctcf

Proper development of the limb bud relies on the concordance of various signals, but its molecular mechanisms have not yet been fully illustrated. Here we report that heterogeneous nuclear ribonucleoprotein K (hnRNPK) is essential for limb bud development. Its ablation in the limb bud results in limbless forelimbs and severe deformities of the hindlimbs. In terms of mechanism, hnRNPK functions as a transcription activator for the vital genes involved in the three regulatory axes of limb bud development. Simultaneously, for the first time we elucidate that hnRNPK binds to and coordinates with the insulator protein CCCTC binding factor (CTCF) to maintain a three-dimensional chromatin architecture. Ablation of hnRNPK weakens the binding strength of CTCF to topologically associating domain (TAD) boundaries, then leading to the loose TADs, and decreased interactions between promoters and enhancers, and further decreased transcription of developmental genes. Our study establishes a fundamental and novel role of hnRNPK in regulating limb bud development.

Sequential and directional insulation by conserved CTCF sites underlies the Hox timer in stembryos

During development, Hox genes are temporally activated according to their relative positions on their clusters, contributing to the proper identities of structures along the rostrocaudal axis. To understand the mechanism underlying this Hox timer, we used mouse embryonic stem cell-derived stembryos. Following Wnt signaling, the process involves transcriptional initiation at the anterior part of the cluster and a concomitant loading of cohesin complexes enriched on the transcribed DNA segments, that is, with an asymmetric distribution favoring the anterior part of the cluster. Chromatin extrusion then occurs with successively more posterior CTCF sites acting as transient insulators, thus generating a progressive time delay in the activation of more posterior-located genes due to long-range contacts with a flanking topologically associating domain. Mutant stembryos support this model and reveal that the presence of evolutionary conserved and regularly spaced intergenic CTCF sites controls the precision and the pace of this temporal mechanism.

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.

Developmental and evolutionary comparative analysis of a regulatory landscape in mouse and chicken

Modifications in gene regulation are driving forces in the evolution of organisms. Part of these changes involve cis-regulatory elements (CREs), which contact their target genes through higher-order chromatin structures. However, how such architectures and variations in CREs contribute to transcriptional evolvability remains elusive. We use Hoxd genes as a paradigm for the emergence of regulatory innovations, as many relevant enhancers are located in a regulatory landscape highly conserved in amniotes. Here, we analysed their regulation in murine vibrissae and chicken feather primordia, two skin appendages expressing different Hoxd gene subsets, and compared the regulation of these genes in these appendages with that in the elongation of the posterior trunk. In the two former structures, distinct subsets of Hoxd genes are contacted by different lineage-specific enhancers, probably as a result of using an ancestral chromatin topology as an evolutionary playground, whereas the gene regulation that occurs in the mouse and chicken embryonic trunk partially relies on conserved CREs. A high proportion of these non-coding sequences active in the trunk have functionally diverged between species, suggesting that transcriptional robustness is maintained, despite considerable divergence in enhancer sequences.

Summary: Analyses of the relationships between chromatin architecture and regulatory activities at the HoxD locus show that ancestral transcription patterns can be maintained while new regulations evolve.

Promoter and enhancer RNAs regulate chromatin reorganization and activation of miR-10b/HOXD locus, and neoplastic transformation in glioma

miR-10b is silenced in normal neuroglial cells of the brain but commonly activated in glioma, where it assumes an essential tumor-promoting role. We demonstrate that the entire miR-10b-hosting HOXD locus is activated in glioma via the cis-acting mechanism involving 3D chromatin reorganization and CTCF-cohesin-mediated looping. This mechanism requires two interacting lncRNAs, HOXD-AS2 and LINC01116, one associated with HOXD3/HOXD4/miR-10b promoter and another with the remote enhancer. Knockdown of either lncRNA in glioma cells alters CTCF and cohesin binding, abolishes chromatin looping, inhibits the expression of all genes within HOXD locus, and leads to glioma cell death. Conversely, in cortical astrocytes, enhancer activation is sufficient for HOXD/miR-10b locus reorganization, gene derepression, and neoplastic cell transformation. LINC01116 RNA is essential for this process. Our results demonstrate the interplay of two lncRNAs in the chromatin folding and concordant regulation of miR-10b and multiple HOXD genes normally silenced in astrocytes and triggering the neoplastic glial transformation.

Sequential in cis mutagenesis in vivo reveals various functions for CTCF sites at the mouse HoxD cluster

Mammalian Hox gene clusters contain a range of CTCF binding sites. In addition to their importance in organizing a TAD border, which isolates the most posterior genes from the rest of the cluster, the positions and orientations of these sites suggest that CTCF may be instrumental in the selection of various subsets of contiguous genes, which are targets of distinct remote enhancers located in the flanking regulatory landscapes. We examined this possibility by producing an allelic series of cumulative in cis mutations in these sites, up to the abrogation of CTCF binding in the five sites located on one side of the TAD border. In the most impactful alleles, the global chromatin architecture of the locus was modified, yet not drastically, illustrating that CTCF sites located on one side of a strong TAD border are sufficient to organize at least part of this insulation. Spatial colinearity in the expression of these genes along the major body axis was nevertheless maintained, despite abnormal expression boundaries. In contrast, strong effects were scored in the selection of target genes responding to particular enhancers, leading to the misregulation of Hoxd genes in specific structures. Altogether, while most enhancer-promoter interactions can occur in the absence of this series of CTCF sites, the binding of CTCF in the Hox cluster is required to properly transform a rather unprecise process into a highly discriminative mechanism of interactions, which is translated into various patterns of transcription accompanied by the distinctive chromatin topology found at this locus. Our allelic series also allowed us to reveal the distinct functional contributions for CTCF sites within this Hox cluster, some acting as insulator elements, others being necessary to anchor or stabilize enhancer-promoter interactions, and some doing both, whereas they all together contribute to the formation of a TAD border. This variety of tasks may explain the amazing evolutionary conservation in the distribution of these sites among paralogous Hox clusters or between various vertebrates.

Conserved and species-specific chromatin remodeling and regulatory dynamics during mouse and chicken limb bud development

Chromatin remodeling and genomic alterations impact spatio-temporal regulation of gene expression, which is central to embryonic development. The analysis of mouse and chicken limb development provides important insights into the morphoregulatory mechanisms, however little is known about the regulatory differences underlying their morphological divergence. Here, we identify the underlying shared and species-specific epigenomic and genomic variations. In mouse forelimb buds, we observe striking synchrony between the temporal dynamics of chromatin accessibility and gene expression, while their divergence in chicken wing buds uncovers species-specific regulatory heterochrony. In silico mapping of transcription factor binding sites and computational footprinting establishes the developmental time-restricted transcription factor-DNA interactions. Finally, the construction of target gene networks for HAND2 and GLI3 transcriptional regulators reveals both conserved and species-specific interactions. Our analysis reveals the impact of genome evolution on the regulatory interactions orchestrating vertebrate limb bud morphogenesis and provides a molecular framework for comparative Evo-Devo studies.

The vertebrate limb bud is a paradigm to uncover the fundamental mechanisms that govern embryogenesis and evolutionary diversification. Here the authors compare mouse and chicken limb bud development to study the impact of genome evolution on conserved and divergent gene regulatory interactions.

Mesomelic dysplasias associated with the HOXD locus are caused by regulatory reallocations

Human families with chromosomal rearrangements at 2q31, where the human HOXD locus maps, display mesomelic dysplasia, a severe shortening and bending of the limb. In mice, the dominant Ulnaless inversion of the HoxD cluster produces a similar phenotype suggesting the same origin for these malformations in humans and mice. Here we engineer 1 Mb inversion including the HoxD gene cluster, which positioned Hoxd13 close to proximal limb enhancers. Using this model, we show that these enhancers contact and activate Hoxd13 in proximal cells, inducing the formation of mesomelic dysplasia. We show that a secondary Hoxd13 null mutation in-cis with the inversion completely rescues the alterations, demonstrating that ectopic HOXD13 is directly responsible for this bone anomaly. Single-cell expression analysis and evaluation of HOXD13 binding sites suggests that the phenotype arises primarily by acting through genes normally controlled by HOXD13 in distal limb cells. Altogether, these results provide a conceptual and mechanistic framework to understand and unify the molecular origins of human mesomelic dysplasia associated with 2q31.

Induction of a chromatin boundary in vivo upon insertion of a TAD border

Mammalian genomes are partitioned into sub-megabase to megabase-sized units of preferential interactions called topologically associating domains or TADs, which are likely important for the proper implementation of gene regulatory processes. These domains provide structural scaffolds for distant cis regulatory elements to interact with their target genes within the three-dimensional nuclear space and architectural proteins such as CTCF as well as the cohesin complex participate in the formation of the boundaries between them. However, the importance of the genomic context in providing a given DNA sequence the capacity to act as a boundary element remains to be fully investigated. To address this question, we randomly relocated a topological boundary functionally associated with the mouse HoxD gene cluster and show that it can indeed act similarly outside its initial genomic context. In particular, the relocated DNA segment recruited the required architectural proteins and induced a significant depletion of contacts between genomic regions located across the integration site. The host chromatin landscape was re-organized, with the splitting of the TAD wherein the boundary had integrated. These results provide evidence that topological boundaries can function independently of their site of origin, under physiological conditions during mouse development.

During development, enhancer sequences tightly regulate the spatio-temporal expression of target genes often located hundreds of kilobases away. This complex process is made possible by the folding of chromatin into domains, which are separated from one another by specific genomic regions referred to as boundaries. In order to understand whether such boundary sequences require their particular genomic contexts to achieve their isolating effect, we analyzed the impact of introducing one such boundary, taken from the HoxD locus, into a distinct topological domain. We show that this ectopic boundary splits the host domain into two sub-domains and affects the expression levels of a neighboring gene. We conclude that this sequence can work independently from its genomic context and thus carries all the information necessary to act as a boundary element.

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.

Developmental hourglass and heterochronic shifts in fin and limb development

How genetic changes are linked to morphological novelties and developmental constraints remains elusive. Here, we investigate genetic apparatuses that distinguish fish fins from tetrapod limbs by analyzing transcriptomes and open-chromatin regions (OCRs). Specifically, we compared mouse forelimb buds with the pectoral fin buds of an elasmobranch, the brown-banded bamboo shark (Chiloscyllium punctatum). A transcriptomic comparison with an accurate orthology map revealed both a mass heterochrony and hourglass-shaped conservation of gene expression between fins and limbs. Furthermore, open-chromatin analysis suggested that access to conserved regulatory sequences is transiently increased during mid-stage limb development. During this stage, stage-specific and tissue-specific OCRs were also enriched. Together, early and late stages of fin/limb development are more permissive to mutations than middle stages, which may have contributed to major morphological changes during the fin-to-limb evolution. We hypothesize that the middle stages are constrained by regulatory complexity that results from dynamic and tissue-specific transcriptional controls.

Animals come in all shapes and sizes. This diversity arose through genetic mutations during evolution, but it is unclear exactly how these variations led to the formation of new shapes. There is increasing evidence to suggest that not all shapes are possible and that variability between animals is limited by a phenomenon known as “developmental constraint”. These limitations direct parts of the body towards a specific shape as they develop in the embryo. Therefore, understanding the mechanisms underlying these developmental constraints could help explain how different body shapes evolved.

The limbs of humans and other mammals evolved from the fins of fish, and this transition is often used to study the role developmental constraints play in evolution. This is an ideal model as there is already a detailed fossil record mapping this evolutionary event, and data pinpointing some of the genes involved in the development of limbs and fins. But this data is incomplete, and a full comparison between the genes activated in the fin and the limb during embryonic development had not been achieved. This is because most fish used for research have undergone recent genetic changes, making it hard to spot which genetic differences are linked to the evolution of the limb.

To overcome this barrier, Onimaru et al. compared genetic data from the developing limbs of mice to the developing fins of the brown-banded bamboo shark, which evolves much slower than other fish. This revealed that although many genes commonly played a role in the development of the fin and the limb in the embryo, the activity of these shared genes was not the same. For example, genes that switched on in the late stages of limb development, switched off in the late stages of fin development. But in the middle of development, those differences were relatively small and both species activated very similar sets of genes. Many of these genes were pleiotropic, which means they have important roles in other tissues and therefore mutate less often. This suggests that the mid-stage of limb development is under the strongest level of constraint.

Darwin’s theory of natural selection explains that mutations drive evolution. But the theory cannot predict what kinds of new body shapes new mutations will produce. Understanding how the activity levels of different genes affect development could help to fill this knowledge gap. This has potential medical applications, for example, understanding why some genetic changes cause more serious problems than others. This work suggests that mutations in genes that are active during the mid-stage of limb development may have the most serious impact.

Making a bat: The developmental basis of bat evolution

Bats are incredibly diverse, both morphologically and taxonomically. Bats are the only mammalian group to have achieved powered flight, an adaptation that is hypothesized to have allowed them to colonize various and diverse ecological niches. However, the lack of fossils capturing the transition from terrestrial mammal to volant chiropteran has obscured much of our understanding of bat evolution. Over the last 20 years, the emergence of evo-devo in non-model species has started to fill this gap by uncovering some developmental mechanisms at the origin of bat diversification. In this review, we highlight key aspects of studies that have used bats as a model for morphological adaptations, diversification during adaptive radiations, and morphological novelty. To do so, we review current and ongoing studies on bat evolution. We first investigate morphological specialization by reviewing current knowledge about wing and face evolution. Then, we explore the mechanisms behind adaptive diversification in various ecological contexts using vision and dentition. Finally, we highlight the emerging work into morphological novelties using bat wing membranes.

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.

Establishment and evolution of heterochromatin

The eukaryotic genome is packaged into transcriptionally active euchromatin and silent heterochromatin, with most studies focused on the former encompassing the majority of protein-coding genes. The recent development of various sequencing techniques has refined this classic dichromatic partition and has better illuminated the composition, establishment, and evolution of this genomic and epigenomic "dark matter" in the context of topologically associated domains and phase-separated droplets. Heterochromatin includes genomic regions that can be densely stained by chemical dyes, which have been shown to be enriched for repetitive elements and epigenetic marks, including H3K9me2/3 and H3K27me3. Heterochromatin is usually replicated late, concentrated at the nuclear periphery or around nucleoli, and usually lacks highly expressed genes; and now it is considered to be as neither genetically inert nor developmentally static. Heterochromatin guards genome integrity against transposon activities and exerts important regulatory functions by targeting beyond its contained genes. Both its nucleotide sequences and regulatory proteins exhibit rapid coevolution between species. In addition, there are dynamic transitions between euchromatin and heterochromatin during developmental and evolutionary processes. We summarize here the ever-changing characteristics of heterochromatin and propose models and principles for the evolutionary transitions of heterochromatin that have been mainly learned from studies of Drosophila and yeast. Finally, we highlight the role of sex chromosomes in studying heterochromatin evolution.

Impact of genome architecture on the functional activation and repression of Hox regulatory landscapes

BACKGROUND

The spatial organization of the mammalian genome relies upon the formation of chromatin domains of various scales. At the level of gene regulation in cis, collections of enhancer sequences define large regulatory landscapes that usually match with the presence of topologically associating domains (TADs). These domains often contain ranges of enhancers displaying similar or related tissue specificity, suggesting that in some cases, such domains may act as coherent regulatory units, with a global on or off state. By using the HoxD gene cluster, which specifies the topology of the developing limbs via highly orchestrated regulation of gene expression, as a paradigm, we investigated how the arrangement of regulatory domains determines their activity and function.

RESULTS

Proximal and distal cells in the developing limb express different levels of Hoxd genes, regulated by flanking 3' and 5' TADs, respectively. We characterized the effect of large genomic rearrangements affecting these two TADs, including their fusion into a single chromatin domain. We show that, within a single hybrid TAD, the activation of both proximal and distal limb enhancers globally occurred as when both TADs are intact. However, the activity of the 3' TAD in distal cells is generally increased in the fused TAD, when compared to wild type where it is silenced. Also, target gene activity in distal cells depends on whether or not these genes had previously responded to proximal enhancers, which determines the presence or absence of H3K27me3 marks. We also show that the polycomb repressive complex 2 is mainly recruited at the Hox gene cluster and can extend its coverage to far-cis regulatory sequences as long as confined to the neighboring TAD structure.

CONCLUSIONS

We conclude that antagonistic limb proximal and distal enhancers can exert their specific effects when positioned into the same TAD and in the absence of their genuine target genes. We also conclude that removing these target genes reduced the coverage of a regulatory landscape by chromatin marks associated with silencing, which correlates with its prolonged activity in time.