Hox genes are required during the morphogenesis of both vertebrate digits and external genitals. We investigated whether transcription in such distinct contexts involves a shared enhancer-containing landscape. We show that the same regulatory topology is used, yet with some tissue-specific enhancer-promoter interactions, suggesting the hijacking of a regulatory backbone from one context to the other. In addition, comparable organizations are observed at both HoxA and HoxD clusters, which separated through genome duplication in an ancestral invertebrate animal. We propose that this convergent regulatory evolution was triggered by the preexistence of some chromatin architecture, thus facilitating the subsequent recruitment of the appropriate transcription factors. Such regulatory topologies may have both favored and constrained the evolution of pleiotropic developmental loci in vertebrates.

No abstract

Polycomb group (PcG) proteins are essential for the repression of key factors during early development. In Drosophila, the polycomb repressive complexes (PRC) associate with defined polycomb response DNA elements (PREs). In mammals, however, the mechanisms underlying polycomb recruitment at targeted loci are poorly understood. We have used an in vivo approach to identify DNA sequences of importance for the proper recruitment of polycomb proteins at the HoxD locus. We report that various genomic re-arrangements of the gene cluster do not strongly affect PRC2 recruitment and that relatively small polycomb interacting sequences appear necessary and sufficient to confer polycomb recognition and targeting to ectopic loci. In addition, a high GC content, while not sufficient to recruit PRC2, may help its local spreading. We discuss the importance of PRC2 recruitment over Hox gene clusters in embryonic stem cells, for their subsequent coordinated transcriptional activation during development.

In many animal species with a bilateral symmetry, Hox genes are clustered either at one or at several genomic loci. This organization has a functional relevance, as the transcriptional control applied to each gene depends upon its relative position within the gene cluster. It was previously noted that vertebrate Hox clusters display a much higher level of genomic organization than their invertebrate counterparts. The former are always more compact than the latter, they are generally devoid of repeats and of interspersed genes, and all genes are transcribed by the same DNA strand, suggesting that particular factors constrained these clusters toward a tighter structure during the evolution of the vertebrate lineage. Here, we investigate the importance of uniform transcriptional orientation by engineering several alleles within the HoxD cluster, such as to invert one or several transcription units, with or without a neighboring CTCF site. We observe that the association between the tight structure of mammalian Hox clusters and their regulation makes inversions likely detrimental to the proper implementation of this complex genetic system. We propose that the consolidation of Hox clusters in vertebrates, including transcriptional polarity, evolved in conjunction with the emergence of global gene regulation via the flanking regulatory landscapes, to optimize a coordinated response of selected subsets of target genes in cis.

Genome-wide association studies for non-syndromic orofacial clefting (OFC) have identified single nucleotide polymorphisms (SNPs) at loci where the presumed risk-relevant gene is expressed in oral periderm. The functional subsets of such SNPs are difficult to predict because the sequence underpinnings of periderm enhancers are unknown. We applied ATAC-seq to models of human palate periderm, including zebrafish periderm, mouse embryonic palate epithelia, and a human oral epithelium cell line, and to complementary mesenchymal cell types. We identified sets of enhancers specific to the epithelial cells and trained gapped-kmer support-vector-machine classifiers on these sets. We used the classifiers to predict the effects of 14 OFC-associated SNPs at 12q13 near KRT18. All the classifiers picked the same SNP as having the strongest effect, but the significance was highest with the classifier trained on zebrafish periderm. Reporter and deletion analyses support this SNP as lying within a periderm enhancer regulating KRT18/KRT8 expression.