Single-cell profiling of transcriptome and histone modifications with EpiDamID.

Recent advances in single-cell sequencing technologies have enabled simultaneous measurement of multiple cellular modalities, but the combined detection of histone post-translational modifications and transcription at single-cell resolution has remained limited. Here, we introduce EpiDamID, an experimental approach to target a diverse set of chromatin types by leveraging the binding specificities of single-chain variable fragment antibodies, engineered chromatin reader domains, and endogenous chromatin-binding proteins. Using these, we render the DamID technology compatible with the genome-wide identification of histone post-translational modifications. Importantly, this includes the possibility to jointly measure chromatin marks and transcription at the single-cell level. We use EpiDamID to profile single-cell Polycomb occupancy in mouse embryoid bodies and provide evidence for hierarchical gene regulatory networks. In addition, we map H3K9me3 in early zebrafish embryogenesis, and detect striking heterochromatic regions specific to notochord. Overall, EpiDamID is a new addition to a vast toolbox to study chromatin states during dynamic cellular processes.

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.

A pendulum of induction between the epiblast and extra-embryonic endoderm supports post-implantation progression.

Embryogenesis is supported by dynamic loops of cellular interactions. Here, we create a partial mouse embryo model to elucidate the principles of epiblast (Epi) and extra-embryonic endoderm co-development (XEn). We trigger naive mouse embryonic stem cells to form a blastocyst-stage niche of Epi-like cells and XEn-like cells (3D, hydrogel free and serum free). Once established, these two lineages autonomously progress in minimal medium to form an inner pro-amniotic-like cavity surrounded by polarized Epi-like cells covered with visceral endoderm (VE)-like cells. The progression occurs through reciprocal inductions by which the Epi supports the primitive endoderm (PrE) to produce a basal lamina that subsequently regulates Epi polarization and/or cavitation, which, in return, channels the transcriptomic progression to VE. This VE then contributes to Epi bifurcation into anterior- and posterior-like states. Similarly, boosting the formation of PrE-like cells within blastoids supports developmental progression. We argue that self-organization can arise from lineage bifurcation followed by a pendulum of induction that propagates over time.

Control of Hoxd gene transcription in the mammary bud by hijacking a preexisting regulatory landscape.

Vertebrate Hox genes encode transcription factors operating during the development of multiple organs and structures. However, the evolutionary mechanism underlying this remarkable pleiotropy remains to be fully understood. Here, we show that Hoxd8 and Hoxd9, two genes of the HoxD complex, are transcribed during mammary bud (MB) development. However, unlike in other developmental contexts, their coexpression does not rely on the same regulatory mechanism. Hoxd8 is regulated by the combined activity of closely located sequences and the most distant telomeric gene desert. On the other hand, Hoxd9 is controlled by an enhancer-rich region that is also located within the telomeric gene desert but has no impact on Hoxd8 transcription, thus constituting an exception to the global regulatory logic systematically observed at this locus. The latter DNA region is also involved in Hoxd gene regulation in other contexts and strongly interacts with Hoxd9 in all tissues analyzed thus far, indicating that its regulatory activity was already operational before the appearance of mammary glands. Within this DNA region and neighboring a strong limb enhancer, we identified a short sequence conserved in therian mammals and capable of enhancer activity in the MBs. We propose that Hoxd gene regulation in embryonic MBs evolved by hijacking a preexisting regulatory landscape that was already at work before the emergence of mammals in structures such as the limbs or the intestinal tract.

Role of a polymorphism in a Hox/Pax-responsive enhancer in the evolution of the vertebrate spine.

Patterning of the vertebrate skeleton requires the coordinated activity of Hox genes. In particular, Hox10 proteins are essential to set the transition from thoracic to lumbar vertebrae because of their rib-repressing activity. In snakes, however, the thoracic region extends well into Hox10-expressing areas of the embryo, suggesting that these proteins are unable to block rib formation. Here, we show that this is not a result of the loss of rib-repressing properties by the snake proteins, but rather to a single base pair change in a Hox/Paired box (Pax)-responsive enhancer, which prevents the binding of Hox proteins. This polymorphism is also found in Paenungulata, such as elephants and manatees, which have extended rib cages. In vivo, this modified enhancer failed to respond to Hox10 activity, supporting its role in the extension of rib cages. In contrast, the enhancer could still interact with Hoxb6 and Pax3 to promote rib formation. These results suggest that a polymorphism in the Hox/Pax-responsive enhancer may have played a role in the evolution of the vertebrate spine by differently modulating its response to rib-suppressing and rib-promoting Hox proteins.

The king cobra genome reveals dynamic gene evolution and adaptation in the snake venom system.

Snakes are limbless predators, and many species use venom to help overpower relatively large, agile prey. Snake venoms are complex protein mixtures encoded by several multilocus gene families that function synergistically to cause incapacitation. To examine venom evolution, we sequenced and interrogated the genome of a venomous snake, the king cobra (Ophiophagus hannah), and compared it, together with our unique transcriptome, microRNA, and proteome datasets from this species, with data from other vertebrates. In contrast to the platypus, the only other venomous vertebrate with a sequenced genome, we find that snake toxin genes evolve through several distinct co-option mechanisms and exhibit surprisingly variable levels of gene duplication and directional selection that correlate with their functional importance in prey capture. The enigmatic accessory venom gland shows a very different pattern of toxin gene expression from the main venom gland and seems to have recruited toxin-like lectin genes repeatedly for new nontoxic functions. In addition, tissue-specific microRNA analyses suggested the co-option of core genetic regulatory components of the venom secretory system from a pancreatic origin. Although the king cobra is limbless, we recovered coding sequences for all Hox genes involved in amniote limb development, with the exception of Hoxd12. Our results provide a unique view of the origin and evolution of snake venom and reveal multiple genome-level adaptive responses to natural selection in this complex biological weapon system. More generally, they provide insight into mechanisms of protein evolution under strong selection.

Regulatory role for a conserved motif adjacent to the homeodomain of Hox10 proteins.

Development of the vertebrate axial skeleton requires the concerted activity of several Hox genes. Among them, Hox genes belonging to the paralog group 10 are essential for the formation of the lumbar region of the vertebral column, owing to their capacity to block rib formation. In this work, we explored the basis for the rib-repressing activity of Hox10 proteins. Because genetic experiments in mice demonstrated that Hox10 proteins are strongly redundant in this function, we first searched for common motifs among the group members. We identified the presence of two small sequences flanking the homeodomain that are phylogenetically conserved among Hox10 proteins and that seem to be specific for this group. We show here that one of these motifs is required but not sufficient for the rib-repressing activity of Hox10 proteins. This motif includes two potential phosphorylation sites, which are essential for protein activity as their mutation to alanines resulted in a total loss of rib-repressing properties. Our data indicates that this motif has a significant regulatory function, modulating interactions with more N-terminal parts of the Hox protein, eventually triggering the rib-repressing program. In addition, this motif might also regulate protein activity by alteration of the protein's DNA-binding affinity through changes in the phosphorylation state of two conserved tyrosine residues within the homeodomain.

Antagonism between H3K27me3 and genome-lamina association drives atypical spatial genome organization in the totipotent embryo.

In mammals, early embryonic development exhibits highly unusual spatial positioning of genomic regions at the nuclear lamina, but the mechanisms underpinning this atypical genome organization remain elusive. Here, we generated single-cell profiles of lamina-associated domains (LADs) coupled with transcriptomics, which revealed a striking overlap between preimplantation-specific LAD dissociation and noncanonical broad domains of H3K27me3. Loss of H3K27me3 resulted in a restoration of canonical LAD profiles, suggesting an antagonistic relationship between lamina association and H3K27me3. Tethering of H3K27me3 to the nuclear periphery showed that the resultant relocalization is partially dependent on the underlying DNA sequence. Collectively, our results suggest that the atypical organization of LADs in early developmental stages is the result of a tug-of-war between intrinsic affinity for the nuclear lamina and H3K27me3, constrained by the available space at the nuclear periphery. This study provides detailed insight into the molecular mechanisms regulating nuclear organization during early mammalian development.

The role of heterochromatin in 3D genome organization during preimplantation development.

During the early stages of mammalian development, the epigenetic state of the parental genome is completely reprogrammed to give rise to the totipotent embryo. An important aspect of this remodeling concerns the heterochromatin and the spatial organization of the genome. While heterochromatin and genome organization are intricately linked in pluripotent and somatic systems, little is known about their relationship in the totipotent embryo. In this review, we summarize the current knowledge on the reprogramming of both regulatory layers. In addition, we discuss available evidence on their relationship and put this in the context of findings in other systems.

Spatial chromatin organization and gene regulation at the nuclear lamina.

The nuclear lamina (NL) consists of a thin meshwork of lamins and associated proteins that lines the inner nuclear membrane (INM). In metazoan nuclei, a large proportion of the genome contacts the NL in broad lamina-associated domains (LADs). Contacts of the NL with the genome are believed to aid the spatial organization of chromosomes and contribute to the regulation of transcription. Here, we will focus on recent insights in the structural organization of the genome at the NL and the role of this organization in the regulation of gene expression.

DamC reveals principles of chromatin folding in vivo without crosslinking and ligation.

Current understanding of chromosome folding is largely reliant on chromosome conformation capture (3C)-based experiments, where chromosomal interactions are detected as ligation products after chromatin crosslinking. To measure chromosome structure in vivo, quantitatively and without crosslinking and ligation, we implemented a modified version of DNA adenine methyltransferase identification (DamID) named DamC, which combines DNA methylation-based detection of chromosomal interactions with next-generation sequencing and biophysical modeling of methylation kinetics. DamC performed in mouse embryonic stem cells provides the first in vivo validation of the existence of topologically associating domains (TADs), CTCF loops and confirms 3C-based measurements of the scaling of contact probabilities. Combining DamC with transposon-mediated genomic engineering shows that new loops can be formed between ectopic and endogenous CTCF sites, which redistributes physical interactions within TADs. DamC provides the first crosslinking- and ligation-free demonstration of the existence of key structural features of chromosomes and provides novel insights into how chromosome structure within TADs can be manipulated.

Reorganisation of Hoxd regulatory landscapes during the evolution of a snake-like body plan.

Within land vertebrate species, snakes display extreme variations in their body plan, characterized by the absence of limbs and an elongated morphology. Such a particular interpretation of the basic vertebrate body architecture has often been associated with changes in the function or regulation of Hox genes. Here, we use an interspecies comparative approach to investigate different regulatory aspects at the snake HoxD locus. We report that, unlike in other vertebrates, snake mesoderm-specific enhancers are mostly located within the HoxD cluster itself rather than outside. In addition, despite both the absence of limbs and an altered Hoxd gene regulation in external genitalia, the limb-associated bimodal HoxD chromatin structure is maintained at the snake locus. Finally, we show that snake and mouse orthologous enhancer sequences can display distinct expression specificities. These results show that vertebrate morphological evolution likely involved extensive reorganisation at Hox loci, yet within a generally conserved regulatory framework.

Snakes: hatching of a model system for Evo-Devo?

Evo-Devo studies rely on a collection of animal model systems belonging to different phylogenetic branches to try and understand how organisms carrying a similar set of genes and pathways can develop into such a variety of shapes and sizes. The squamate clade, however, has only recently started to receive the attention it deserves in particular due to extreme morphological and metabolic aspects and, consequently, the important insights that it could bring in different fields. The recent sequencing of several squamate genomes as well as the generation of high quality trancriptomes for different snake tissues now provide the necessary tools to complement biological studies. Here, we briefly report on recent work involving developing snake embryos to illustrate their interest to assess vertebrate developmental mechanisms. We also discuss the relevance to use snake species as Evo-Devo model systems and potential ways to cross the important limitations intrinsically associated with developmental and genetic studies of these fascinating animals.