Daniel M. Ibrahim

Formation of new chromatin domains determines pathogenicity of genomic duplications

Chromosome conformation capture methods have identified subchromosomal structures of higher-order chromatin interactions called topologically associated domains (TADs) that are separated from each other by boundary regions. By subdividing the genome into discrete regulatory units, TADs restrict the contacts that enhancers establish with their target genes. However, the mechanisms that underlie partitioning of the genome into TADs remain poorly understood. Here we show by chromosome conformation capture (capture Hi-C and 4C-seq methods) that genomic duplications in patient cells and genetically modified mice can result in the formation of new chromatin domains (neo-TADs) and that this process determines their molecular pathology. Duplications of non-coding DNA within the mouse Sox9 TAD (intra-TAD) that cause female to male sex reversal in humans, showed increased contact of the duplicated regions within the TAD, but no change in the overall TAD structure. In contrast, overlapping duplications that extended over the next boundary into the neighbouring TAD (inter-TAD), resulted in the formation of a new chromatin domain (neo-TAD) that was isolated from the rest of the genome. As a consequence of this insulation, inter-TAD duplications had no phenotypic effect. However, incorporation of the next flanking gene, Kcnj2, in the neo-TAD resulted in ectopic contacts of Kcnj2 with the duplicated part of the Sox9 regulatory region, consecutive misexpression of Kcnj2, and a limb malformation phenotype. Our findings provide evidence that TADs are genomic regulatory units with a high degree of internal stability that can be sculptured by structural genomic variations. This process is important for the interpretation of copy number variations, as these variations are routinely detected in diagnostic tests for genetic disease and cancer. This finding also has relevance in an evolutionary setting because copy-number differences are thought to have a crucial role in the evolution of genome complexity.

Regulation of cell polarity in the cartilage growth plate and perichondrium of metacarpal elements by HOXD13 and WNT5A

The morphology of bones is genetically determined, but the molecular mechanisms that control shape, size and the overall gestalt of bones remain unclear. We previously showed that metacarpals in the synpolydactyly homolog (spdh) mouse, which carries a mutation in Hoxd13 similar to the human condition synpolydactyly (SPD), were transformed to carpal-like bones with cuboid shape that lack cortical bone and a perichondrium and are surrounded by a joint surface. Here we provide evidence that spdh metacarpal growth plates have a defect in cell polarization with a random instead of linear orientation. In parallel prospective perichondral cells failed to adopt the characteristic flattened cell shape. We observed a similar cell polarity defect in metacarpals of Wnt5a(-/-) mice. Wnt5a and the closely related Wnt5b were downregulated in spdh handplates, and HOXD13 induced expression of both genes in vitro. Concomitant we observed mislocalization of core planar cell polarity (PCP) components DVL2 and PRICKLE1 in spdh metacarpals indicating a defect in the WNT/PCP pathway. Conversely the WNT/β-CATENIN pathway, a hallmark of joint cells lining carpal bones, was upregulated in the perichondral region. Finally, providing spdh limb explant cultures with cells expressing either HOXD13 or WNT5A led to a non-cell autonomous partial rescue of cell polarity the perichondral region and restored the expression of perichondral markers. This study provides a so far unrecognized link between HOX proteins and cell polarity in the perichondrium and the growth plate, a failure of which leads to transformation of metacarpals to carpal-like structures.

Characterization of hundreds of regulatory landscapes in developing limbs reveals two regimes of chromatin folding

Complex regulatory landscapes control the pleiotropic transcriptional activities of developmental genes. For most genes, the number, location, and dynamics of their associated regulatory elements are unknown. In this work, we characterized the three-dimensional chromatin microarchitecture and regulatory landscape of 446 limb-associated gene loci in mouse using Capture-C, ChIP-seq, and RNA-seq in forelimb, hindlimb at three developmental stages, and midbrain. The fine mapping of chromatin interactions revealed a strong preference for functional genomic regions such as repressed or active domains. By combining chromatin marks and interaction peaks, we annotated more than 1000 putative limb enhancers and their associated genes. Moreover, the analysis of chromatin interactions revealed two regimes of chromatin folding, one producing interactions stable across tissues and stages and another one associated with tissue and/or stage-specific interactions. Whereas stable interactions associate strongly with CTCF/RAD21 binding, the intensity of variable interactions correlates with changes in underlying chromatin modifications, specifically at the viewpoint and at the interaction site. In conclusion, this comprehensive data set provides a resource for the characterization of hundreds of limb-associated regulatory landscapes and a framework to interpret the chromatin folding dynamics observed during embryogenesis.

Distinct global shifts in genomic binding profiles of limb malformation-associated HOXD13 mutations

Gene regulation by transcription factors (TFs) determines developmental programs and cell identity. Consequently, mutations in TFs can lead to dramatic phenotypes in humans by disrupting gene regulation. To date, the molecular mechanisms that actually cause these phenotypes have been difficult to address experimentally. ChIP-seq, which couples chromatin immunoprecipitation with high-throughput sequencing, allows TF function to be investigated on a genome-wide scale, enabling new approaches for the investigation of gene regulation. Here, we present the application of ChIP-seq to explore the effect of missense mutations in TFs on their genome-wide binding profile. Using a retroviral expression system in chicken mesenchymal stem cells, we elucidated the mechanism underlying a novel missense mutation in HOXD13 (Q317K) associated with a complex hand and foot malformation phenotype. The mutated glutamine (Q) is conserved in most homeodomains, a notable exception being bicoid-type homeodomains that have lysine (K) at this position. Our results show that the mutation results in a shift in the binding profile of the mutant toward a bicoid/PITX1 motif. Gene expression analysis and functional assays using in vivo overexpression studies confirm that the mutation results in a partial conversion of HOXD13 into a TF with bicoid/PITX1 properties. A similar shift was not observed with another mutation, Q317R, which is associated with brachysyndactyly, suggesting that the bicoid/PITX1-shift observed for Q317K might be related to the severe clinical phenotype. The methodology described can be used to investigate a wide spectrum of TFs and mutations that have not previously been amenable to ChIP-seq experiments.

Genome-wide binding of posterior HOXA/D transcription factors reveals subgrouping and association with CTCF

Homeotic genes code for key transcription factors (HOX-TFs) that pattern the animal body plan. During embryonic development, Hox genes are expressed in overlapping patterns and function in a partially redundant manner. In vitro biochemical screens probing the HOX-TF sequence specificity revealed largely overlapping sequence preferences, indicating that co-factors might modulate the biological function of HOX-TFs. However, due to their overlapping expression pattern, high protein homology, and insufficiently specific antibodies, little is known about their genome-wide binding preferences. In order to overcome this problem, we virally expressed tagged versions of limb-expressed posterior HOX genes (HOXA9-13, and HOXD9-13) in primary chicken mesenchymal limb progenitor cells (micromass). We determined the effect of each HOX-TF on cellular differentiation (chondrogenesis) and gene expression and found that groups of HOX-TFs induce distinct regulatory programs. We used ChIP-seq to determine their individual genome-wide binding profiles and identified between 12,721 and 28,572 binding sites for each of the nine HOX-TFs. Principal Component Analysis (PCA) of binding profiles revealed that the HOX-TFs are clustered in two subgroups (Group 1: HOXA/D9, HOXA/D10, HOXD12, and HOXA13 and Group 2: HOXA/D11 and HOXD13), which are characterized by differences in their sequence specificity and by the presence of cofactor motifs. Specifically, we identified CTCF binding sites in Group 1, indicating that this subgroup of HOX-proteins cooperates with CTCF. We confirmed this interaction by an independent biological assay (Proximity Ligation Assay) and demonstrated that CTCF is a novel HOX cofactor that specifically associates with Group 1 HOX-TFs, pointing towards a possible interplay between HOX-TFs and chromatin architecture.

Odd skipped-related 1 identifies a population of embryonic fibro-adipogenic progenitors regulating myogenesis during limb development

Fibro-adipogenic progenitors (FAPs) are an interstitial cell population in adult skeletal muscle that support muscle regeneration. During development, interstitial muscle connective tissue (MCT) cells support proper muscle patterning, however the underlying molecular mechanisms are not well understood and it remains unclear whether adult FAPs and embryonic MCT cells share a common lineage. We show here that mouse embryonic limb MCT cells expressing the transcription factor Osr1, differentiate into fibrogenic and adipogenic cells in vivo and in vitro defining an embryonic FAP-like population. Genetic lineage tracing shows that developmental Osr1+ cells give rise to a subset of adult FAPs. Loss of Osr1 function leads to a reduction of myogenic progenitor proliferation and survival resulting in limb muscle patterning defects. Transcriptome and functional analyses reveal that Osr1+ cells provide a critical pro-myogenic niche via the production of MCT specific extracellular matrix components and secreted signaling factors.Fibro-adipogenic progenitors (FAPs) form part of interstitial muscle connective tissue (MCT) in adults but the origin of this non-myogenic lineage is unclear. Here, the authors show that Odd skipped related 1 (Osr1) in mice marks embryonic MCT, giving rise to FAPs, and loss of Osr1 in the limb causes muscle defects.

A homozygous HOXD13 missense mutation causes a severe form of synpolydactyly with metacarpal to carpal transformation.

Synpolydactyly (SPD) is a rare congenital limb disorder characterized by syndactyly between the third and fourth fingers and an additional digit in the syndactylous web. In most cases SPD is caused by heterozygous mutations in HOXD13 resulting in the expansion of a N‐terminal polyalanine tract. If homozygous, the mutation results in severe shortening of all metacarpals and phalanges with a morphological transformation of metacarpals to carpals. Here, we describe a novel homozygous missense mutation in a family with unaffected consanguineous parents and severe brachydactyly and metacarpal‐to‐carpal transformation in the affected child. We performed whole exome sequencing on the index patient, followed by Sanger sequencing of parents and patient to investigate cosegregation. The DNA‐binding ability of the mutant protein was tested with electrophoretic mobility shift assays. We demonstrate that the c.938C>G (p.313T>R) mutation in the DNA‐binding domain of HOXD13 prevents binding to DNA in vitro. Our results show to our knowledge for the first time that a missense mutation in HOXD13 underlies severe brachydactyly with metacarpal‐to‐carpal transformation. The mutation is non‐penetrant in heterozygous carriers. In conjunction with the literature we propose the possibility that the metacarpal‐to‐carpal transformation results from a homozygous loss of functional HOXD13 protein in humans in combination with an accumulation of non‐functional HOXD13 that might be able to interact with other transcription factors in the developing limb.

Missense-Mutationen in Transkriptionsfaktoren

ZusammenfassungTranskriptionsfaktoren sind entscheidende Regulatoren der Embryonalentwicklung, da sie die Genexpression in jeder Zelle kontrollieren. Mutationen in Transkriptionsfaktoren liegen häufig angeborenen Entwicklungsdefekten zugrunde, jedoch ist die funktionelle Einschätzung der Pathogenität einzelner Transkriptionsfaktorvarianten anspruchsvoll, da die molekulare Funktionsweise von Transkriptionsfaktoren nicht vollkommen verstanden ist. Besonders Gain-of-Function-Mutationen führen häufig zu neuen, unerwarteten Phänotypen, deren funktionelle Charakterisierung eine Herausforderung darstellt. Die im letzten Jahrzehnt entwickelte ChIP-seq-Technologie ermöglicht es, die molekularen Mechanismen zu unterscheiden, welche Transkriptionsfaktor-assoziierten Krankheiten zugrunde liegen. Dieser Artikel fasst die molekularen Pathomechanismen diverser Transkriptionsfaktormutationen zusammen und versucht einen molekularbiologischen Rahmen für die Bewertung neuer Transkriptionsfaktormutationen zu geben.AbstractTranscription factors (TF) are key regulators that control the cell-type-specific gene expression in each individual cell and are crucial for the coordination of embryonic development. Mutations in TFs frequently underlie heritable developmental defects; however, the functional characterization of a variant TF is challenging, since the exact molecular mechanisms by which TFs exert their function are not fully understood. In particular, gain-of-function mutations can lead to novel phenotypes that are difficult to characterize functionally. Recent technological advances, in particular ChIP-seq, have enabled experimental approaches that can distinguish between distinct molecular mechanisms underlying TF-associated diseases. This article reviews the molecular pathomechanisms underlying various TF mutations and proposes approaches to previously unknown TF mutations.