Opposing functions of the ETS factor family define Shh spatial expression in limb buds and underlie polydactyly

Mutagenesis Sensitivity Mapping of Human Enhancers In Vivo

Distant-acting enhancers are central to human development. However, our limited understanding of their functional sequence features prevents the interpretation of enhancer mutations in disease. Here, we determined the functional sensitivity to mutagenesis of human developmental enhancers in vivo. Focusing on seven enhancers active in the developing brain, heart, limb and face, we created over 1700 transgenic mice for over 260 mutagenized enhancer alleles. Systematic mutation of 12-basepair blocks collectively altered each sequence feature in each enhancer at least once. We show that 69% of all blocks are required for normal in vivo activity, with mutations more commonly resulting in loss (60%) than in gain (9%) of function. Using predictive modeling, we annotated critical nucleotides at base-pair resolution. The vast majority of motifs predicted by these machine learning models (88%) coincided with changes to in vivo function, and the models showed considerable sensitivity, identifying 59% of all functional blocks. Taken together, our results reveal that human enhancers contain a high density of sequence features required for their normal in vivo function and provide a rich resource for further exploration of human enhancer logic.

DNA-binding factor footprints and enhancer RNAs identify functional non-coding genetic variants

BACKGROUND

Genome-wide association studies (GWAS) have revealed a multitude of candidate genetic variants affecting the risk of developing complex traits and diseases. However, the highlighted regions are typically in the non-coding genome, and uncovering the functional causative single nucleotide variants (SNVs) is challenging. Prioritization of variants is commonly based on genomic annotation with markers of active regulatory elements, but current approaches still poorly predict functional variants. To address this, we systematically analyze six markers of active regulatory elements for their ability to identify functional variants.

RESULTS

We benchmark against molecular quantitative trait loci (molQTL) from assays of regulatory element activity that identify allelic effects on DNA-binding factor occupancy, reporter assay expression, and chromatin accessibility. We identify the combination of DNase footprints and divergent enhancer RNA (eRNA) as markers for functional variants. This signature provides high precision, but with a trade-off of low recall, thus substantially reducing candidate variant sets to prioritize variants for functional validation. We present this as a framework called FINDER-Functional SNV IdeNtification using DNase footprints and eRNA.

CONCLUSIONS

We demonstrate the utility to prioritize variants using leukocyte count trait and analyze variants in linkage disequilibrium with a lead variant to predict a functional variant in asthma. Our findings have implications for prioritizing variants from GWAS, in development of predictive scoring algorithms, and for functionally informed fine mapping approaches.

TBX3 is essential for establishment of the posterior boundary of anterior genes and upregulation of posterior genes together with HAND2 during the onset of limb bud development

During limb bud formation, axis polarities are established as evidenced by the spatially restricted expression of key regulator genes. In particular, the mutually antagonistic interaction between the GLI3 repressor and HAND2 results in distinct and non-overlapping anterior-distal Gli3 and posterior Hand2 expression domains. This is a hallmark of the establishment of antero-posterior limb axis polarity, together with spatially restricted expression of homeodomain and other transcriptional regulators. Here, we show that TBX3 is required for establishment of the posterior expression boundary of anterior genes in mouse limb buds. ChIP-seq and differential gene expression analysis of wild-type and mutant limb buds identifies TBX3-specific and shared TBX3-HAND2 target genes. High sensitivity fluorescent whole-mount in situ hybridisation shows that the posterior expression boundaries of anterior genes are positioned by TBX3-mediated repression, which excludes anterior genes such as Gli3, Alx4, Hand1 and Irx3/5 from the posterior limb bud mesenchyme. This exclusion delineates the posterior mesenchymal territory competent to establish the Shh-expressing limb bud organiser. In turn, HAND2 is required for Shh activation and cooperates with TBX3 to upregulate shared posterior identity target genes in early limb buds.

Conserved Cis-Acting Range Extender Element Mediates Extreme Long-Range Enhancer Activity in Mammals

While most mammalian enhancers regulate their cognate promoters over moderate distances of tens of kilobases (kb), some enhancers act over distances in the megabase range. The sequence features enabling such extreme-distance enhancer-promoter interactions remain elusive. Here, we used in vivo enhancer replacement experiments in mice to show that short- and medium-range enhancers cannot initiate gene expression at extreme-distance range. We uncover a novel conserved cis-acting element, Range EXtender (REX), that confers extreme-distance regulatory activity and is located next to a long-range enhancer of Sall1. The REX element itself has no endogenous enhancer activity. However, addition of the REX to other short- and mid-range enhancers substantially increases their genomic interaction range. In the most extreme example observed, addition of the REX increased the range of an enhancer by an order of magnitude, from its native 71kb to 840kb. The REX element contains highly conserved [C/T]AATTA homeodomain motifs. These motifs are enriched around long-range limb enhancers genome-wide, including the ZRS, a benchmark long-range limb enhancer of Shh. Mutating the [C/T]AATTA motifs within the ZRS does not affect its limb-specific enhancer activity at short range, but selectively abolishes its long-range activity, resulting in severe limb reduction in knock-in mice. In summary, we identify a sequence signature globally associated with long-range enhancer-promoter interactions and describe a prototypical REX element that is necessary and sufficient to confer extreme-distance gene activation by remote enhancers.

Proliferation-driven mechanical compression induces signalling centre formation during mammalian organ development

Localized sources of morphogens, called signalling centres, play a fundamental role in coordinating tissue growth and cell fate specification during organogenesis. However, how these signalling centres are established in tissues during embryonic development is still unclear. Here we show that the main signalling centre orchestrating development of rodent incisors, the enamel knot (EK), is specified by a cell proliferation-driven buildup in compressive stresses (mechanical pressure) in the tissue. Direct mechanical measurements indicate that the stresses generated by cell proliferation are resisted by the surrounding tissue, creating a circular pattern of mechanical anisotropy with a region of high compressive stress at its centre that becomes the EK. Pharmacological inhibition of proliferation reduces stresses and suppresses EK formation, and application of external pressure in proliferation-inhibited conditions rescues the formation of the EK. Mechanical information is relayed intracellularly through YAP protein localization, which is cytoplasmic in the region of compressive stress that establishes the EK and nuclear in the stretched anisotropic cells that resist the pressure buildup around the EK. Together, our data identify a new role for proliferation-driven mechanical compression in the specification of a model signalling centre during mammalian organ development.

Affinity-optimizing enhancer variants disrupt development

Enhancers control the location and timing of gene expression and contain the majority of variants associated with disease1-3. The ZRS is arguably the most well-studied vertebrate enhancer and mediates the expression of Shh in the developing limb4. Thirty-one human single-nucleotide variants (SNVs) within the ZRS are associated with polydactyly4-6. However, how this enhancer encodes tissue-specific activity, and the mechanisms by which SNVs alter the number of digits, are poorly understood. Here we show that the ETS sites within the ZRS are low affinity, and identify a functional ETS site, ETS-A, with extremely low affinity. Two human SNVs and a synthetic variant optimize the binding affinity of ETS-A subtly from 15% to around 25% relative to the strongest ETS binding sequence, and cause polydactyly with the same penetrance and severity. A greater increase in affinity results in phenotypes that are more penetrant and more severe. Affinity-optimizing SNVs in other ETS sites in the ZRS, as well as in ETS, interferon regulatory factor (IRF), HOX and activator protein 1 (AP-1) sites within a wide variety of enhancers, cause gain-of-function gene expression. The prevalence of binding sites with suboptimal affinity in enhancers creates a vulnerability in genomes whereby SNVs that optimize affinity, even slightly, can be pathogenic. Searching for affinity-optimizing SNVs in genomes could provide a mechanistic approach to identify causal variants that underlie enhanceropathies.

Single-nucleotide variants within heart enhancers increase binding affinity and disrupt heart development

Transcriptional enhancers direct precise gene expression patterns during development and harbor the majority of variants associated with phenotypic diversity, evolutionary adaptations, and disease. Pinpointing which enhancer variants contribute to changes in gene expression and phenotypes is a major challenge. Here, we find that suboptimal or low-affinity binding sites are necessary for precise gene expression during heart development. Single-nucleotide variants (SNVs) can optimize the affinity of ETS binding sites, causing gain-of-function (GOF) gene expression, cell migration defects, and phenotypes as severe as extra beating hearts in the marine chordate Ciona robusta. In human induced pluripotent stem cell (iPSC)-derived cardiomyocytes, a SNV within a human GATA4 enhancer increases ETS binding affinity and causes GOF enhancer activity. The prevalence of suboptimal-affinity sites within enhancers creates a vulnerability whereby affinity-optimizing SNVs can lead to GOF gene expression, changes in cellular identity, and organismal-level phenotypes that could contribute to the evolution of novel traits or diseases.

Widespread regulatory specificities between transcriptional co-repressors and enhancers in Drosophila

Gene expression is controlled by the precise activation and repression of transcription. Repression is mediated by specialized transcription factors (TFs) that recruit co-repressors (CoRs) to silence transcription, even in the presence of activating cues. However, whether CoRs can dominantly silence all enhancers or display distinct specificities is unclear. In this work, we report that most enhancers in Drosophila can be repressed by only a subset of CoRs, and enhancers classified by CoR sensitivity show distinct chromatin features, function, TF motifs, and binding. Distinct TF motifs render enhancers more resistant or sensitive to specific CoRs, as we demonstrate by motif mutagenesis and addition. These CoR-enhancer compatibilities constitute an additional layer of regulatory specificity that allows differential regulation at close genomic distances and is indicative of distinct mechanisms of transcriptional repression.

A spatio-temporally constrained gene regulatory network directed by PBX1/2 acquires limb patterning specificity via HAND2

A lingering question in developmental biology has centered on how transcription factors with widespread distribution in vertebrate embryos can perform tissue-specific functions. Here, using the murine hindlimb as a model, we investigate the elusive mechanisms whereby PBX TALE homeoproteins, viewed primarily as HOX cofactors, attain context-specific developmental roles despite ubiquitous presence in the embryo. We first demonstrate that mesenchymal-specific loss of PBX1/2 or the transcriptional regulator HAND2 generates similar limb phenotypes. By combining tissue-specific and temporally controlled mutagenesis with multi-omics approaches, we reconstruct a gene regulatory network (GRN) at organismal-level resolution that is collaboratively directed by PBX1/2 and HAND2 interactions in subsets of posterior hindlimb mesenchymal cells. Genome-wide profiling of PBX1 binding across multiple embryonic tissues further reveals that HAND2 interacts with subsets of PBX-bound regions to regulate limb-specific GRNs. Our research elucidates fundamental principles by which promiscuous transcription factors cooperate with cofactors that display domain-restricted localization to instruct tissue-specific developmental programs.

Recruitment of transcription factor ETS1 to activated accessible regions promotes the transcriptional program of cilia genes

Defects in cilia genes, which are critical for cilia formation and function, can cause complicated ciliopathy syndromes involving multiple organs and tissues; however, the underlying regulatory mechanisms of the networks of cilia genes in ciliopathies remain enigmatic. Herein, we have uncovered the genome-wide redistribution of accessible chromatin regions and extensive alterations of expression of cilia genes during Ellis-van Creveld syndrome (EVC) ciliopathy pathogenesis. Mechanistically, the distinct EVC ciliopathy-activated accessible regions (CAAs) are shown to positively regulate robust changes in flanking cilia genes, which are a key requirement for cilia transcription in response to developmental signals. Moreover, a single transcription factor, ETS1, can be recruited to CAAs, leading to prominent chromatin accessibility reconstruction in EVC ciliopathy patients. In zebrafish, the collapse of CAAs driven by ets1 suppression subsequently causes defective cilia proteins, resulting in body curvature and pericardial oedema. Our results depict a dynamic landscape of chromatin accessibility in EVC ciliopathy patients, and uncover an insightful role for ETS1 in controlling the global transcriptional program of cilia genes by reprogramming the widespread chromatin state.

The shh limb enhancer is activated in patterned limb regeneration but not in hypomorphic limb regeneration in Xenopus laevis

Xenopus young tadpoles regenerate a limb with the anteroposterior (AP) pattern, but metamorphosed froglets regenerate a hypomorphic limb after amputation. The key gene for AP patterning, shh, is expressed in a regenerating limb of the tadpole but not in that of the froglet. Genomic DNA in the shh limb-specific enhancer, MFCS1 (ZRS), is hypermethylated in froglets but hypomethylated in tadpoles: shh expression may be controlled by epigenetic regulation of MFCS1. Is MFCS1 specifically activated for regenerating the AP-patterned limb? We generated transgenic Xenopus laevis lines that visualize the MFCS1 enhancer activity with a GFP reporter. The transgenic tadpoles showed GFP expression in hoxd13-and shh-expressing domains of developing and regenerating limbs, whereas the froglets showed no GFP expression in the regenerating limbs despite having hoxd13 expression. Genome sequence analysis and co-transfection assays using cultured cells revealed that Hoxd13 can activate Xenopus MFCS1. These results suggest that MFCS1 activation correlates with regeneration of AP-patterned limbs and that re-activation of epigenetically inactivated MFCS1 would be crucial to confer the ability to non-regenerative animals for regenerating a properly patterned limb.

Genome-wide screening reveals the genetic basis of mammalian embryonic eye development

BACKGROUND

Microphthalmia, anophthalmia, and coloboma (MAC) spectrum disease encompasses a group of eye malformations which play a role in childhood visual impairment. Although the predominant cause of eye malformations is known to be heritable in nature, with 80% of cases displaying loss-of-function mutations in the ocular developmental genes OTX2 or SOX2, the genetic abnormalities underlying the remaining cases of MAC are incompletely understood. This study intended to identify the novel genes and pathways required for early eye development. Additionally, pathways involved in eye formation during embryogenesis are also incompletely understood. This study aims to identify the novel genes and pathways required for early eye development through systematic forward screening of the mammalian genome.

RESULTS

Query of the International Mouse Phenotyping Consortium (IMPC) database (data release 17.0, August 01, 2022) identified 74 unique knockout lines (genes) with genetically associated eye defects in mouse embryos. The vast majority of eye abnormalities were small or absent eyes, findings most relevant to MAC spectrum disease in humans. A literature search showed that 27 of the 74 lines had previously published knockout mouse models, of which only 15 had ocular defects identified in the original publications. These 12 previously published gene knockouts with no reported ocular abnormalities and the 47 unpublished knockouts with ocular abnormalities identified by the IMPC represent 59 genes not previously associated with early eye development in mice. Of these 59, we identified 19 genes with a reported human eye phenotype. Overall, mining of the IMPC data yielded 40 previously unimplicated genes linked to mammalian eye development. Bioinformatic analysis showed that several of the IMPC genes colocalized to several protein anabolic and pluripotency pathways in early eye development. Of note, our analysis suggests that the serine-glycine pathway producing glycine, a mitochondrial one-carbon donator to folate one-carbon metabolism (FOCM), is essential for eye formation.

CONCLUSIONS

Using genome-wide phenotype screening of single-gene knockout mouse lines, STRING analysis, and bioinformatic methods, this study identified genes heretofore unassociated with MAC phenotypes providing models to research novel molecular and cellular mechanisms involved in eye development. These findings have the potential to hasten the diagnosis and treatment of this congenital blinding disease.

A large, ten-generation family with autosomal dominant preaxial polydactyly/triphalangeal thumb: Historical, clinical, genealogical, and molecular studies

We present a large, ten-generation family of 273 individuals with 84 people having preaxial polydactyly/triphalangeal thumb due to a pathogenic variant in the zone of polarizing activity regulatory sequence (ZRS) within the exon 5 of LMBR1. The causative change maps to position 396 of the ZRS, located at position c.423 + 4909C > T (chr7:156791480; hg38; LMBR1 ENST00000353442.10; rs606231153 NG_009240.2) in the intron 5 of LMBR1. The first affected individual with the disorder was traced back to mid-1700, when some settlers and workers established in Cervera de Buitrago, a small village about 82 km North to Madrid. Clinical and radiological studies of most of the affected members have been performed for 42 years (follow-up of the family by LFGA). Molecular studies have confirmed a pathogenic variant in the ZRS that segregates in this family. To the best of our knowledge, this is the largest family with preaxial polydactyly/triphalangeal thumb reported so far.

TATTOO-seq delineates spatial and cell type-specific regulatory programs in the developing limb

The coordinated differentiation of progenitor cells into specialized cell types and their spatial organization into distinct domains is central to embryogenesis. Here, we developed and applied an unbiased spatially resolved single-cell transcriptomics method to identify the genetic programs underlying the emergence of specialized cell types during mouse limb development and their spatial integration. We identify multiple transcription factors whose expression patterns are predominantly associated with cell type specification or spatial position, suggesting two parallel yet highly interconnected regulatory systems. We demonstrate that the embryonic limb undergoes a complex multiscale reorganization upon perturbation of one of its spatial organizing centers, including the loss of specific cell populations, alterations of preexisting cell states' molecular identities, and changes in their relative spatial distribution. Our study shows how multidimensional single-cell, spatially resolved molecular atlases can allow the deconvolution of spatial identity and cell fate and reveal the interconnected genetic networks that regulate organogenesis and its reorganization upon genetic alterations.

Function and Constraint in Enhancer Sequences with Multiple Evolutionary Origins

Thousands of human gene regulatory enhancers are composed of sequences with multiple evolutionary origins. These evolutionarily "complex" enhancers consist of older "core" sequences and younger "derived" sequences. However, the functional relationship between the sequences of different evolutionary origins within complex enhancers is poorly understood. We evaluated the function, selective pressures, and sequence variation across core and derived components of human complex enhancers. We find that both components are older than expected from the genomic background, and complex enhancers are enriched for core and derived sequences of similar evolutionary ages. Both components show strong evidence of biochemical activity in massively parallel report assays. However, core and derived sequences have distinct transcription factor (TF)-binding preferences that are largely similar across evolutionary origins. As expected, given these signatures of function, both core and derived sequences have substantial evidence of purifying selection. Nonetheless, derived sequences exhibit weaker purifying selection than adjacent cores. Derived sequences also tolerate more common genetic variation and are enriched compared with cores for expression quantitative trait loci associated with gene expression variability in human populations. In conclusion, both core and derived sequences have strong evidence of gene regulatory function, but derived sequences have distinct constraint profiles, TF-binding preferences, and tolerance to variation compared with cores. We propose that the step-wise integration of younger derived with older core sequences has generated regulatory substrates with robust activity and the potential for functional variation. Our analyses demonstrate that synthesizing study of enhancer evolution and function can aid interpretation of regulatory sequence activity and functional variation across human populations.

Natural and Experimental Rewiring of Gene Regulatory Regions

The successful development and ongoing functioning of complex organisms depend on the faithful execution of the genetic code. A critical step in this process is the correct spatial and temporal expression of genes. The highly orchestrated transcription of genes is controlled primarily by cis-regulatory elements: promoters, enhancers, and insulators. The medical importance of this key biological process can be seen by the frequency with which mutations and inherited variants that alter cis-regulatory elements lead to monogenic and complex diseases and cancer. Here, we provide an overview of the methods available to characterize and perturb gene regulatory circuits. We then highlight mechanisms through which regulatory rewiring contributes to disease, and conclude with a perspective on how our understanding of gene regulation can be used to improve human health.

Pancreas agenesis mutations disrupt a lead enhancer controlling a developmental enhancer cluster

Sequence variants in cis-acting enhancers are important for polygenic disease, but their role in Mendelian disease is poorly understood. Redundancy between enhancers that regulate the same gene is thought to mitigate the pathogenic impact of enhancer mutations. Recent findings, however, have shown that loss-of-function mutations in a single enhancer near PTF1A cause pancreas agenesis and neonatal diabetes. Using mouse and human genetic models, we show that this enhancer activates an entire PTF1A enhancer cluster in early pancreatic multipotent progenitors. This leading role, therefore, precludes functional redundancy. We further demonstrate that transient expression of PTF1A in multipotent progenitors sets in motion an epigenetic cascade that is required for duct and endocrine differentiation. These findings shed insights into the genome regulatory mechanisms that drive pancreas differentiation. Furthermore, they reveal an enhancer that acts as a regulatory master key and is thus vulnerable to pathogenic loss-of-function mutations.

•

The pancreas agenesis enhancer (EnhP) activates PTF1A in early pancreatic progenitors

•

EnhP also activates other progenitor PTF1A enhancers

•

This master key function explains why EnhP is vulnerable to loss-of-function mutations

•

Transient PTF1A expression in progenitors controls pancreas growth and endocrinogenesis

A Ctnnb1 enhancer regulates neocortical neurogenesis by controlling the abundance of intermediate progenitors

β-catenin-dependent canonical Wnt signaling plays a plethora of roles in neocortex (Ncx) development, but its function in regulating the abundance of intermediate progenitors (IPs) is elusive. Here we identified neCtnnb1, an evolutionarily conserved cis-regulatory element with typical enhancer features in developing Ncx. neCtnnb1 locates 55 kilobase upstream of and spatially close to the promoter of Ctnnb1, the gene encoding β-catenin. CRISPR/Cas9-mediated activation or interference of the neCtnnb1 locus enhanced or inhibited transcription of Ctnnb1. neCtnnb1 drove transcription predominantly in the subventricular zone of developing Ncx. Knock-out of neCtnnb1 in mice resulted in compromised expression of Ctnnb1 and the Wnt reporter in developing Ncx. Importantly, knock-out of neCtnnb1 lead to reduced production and transit-amplification of IPs, which subsequently generated fewer upper-layer Ncx projection neurons (PNs). In contrast, enhancing the canonical Wnt signaling by stabilizing β-catenin in neCtnnb1-active cells promoted the production of IPs and upper-layer Ncx PNs. ASH2L was identified as the key trans-acting factor that associates with neCtnnb1 and Ctnnb1’s promoter to maintain Ctnnb1’s transcription in both mouse and human Ncx progenitors. These findings advance understanding of transcriptional regulation of Ctnnb1, and provide insights into mechanisms underlying Ncx expansion during development.

Etv2 regulates enhancer chromatin status to initiate Shh expression in the limb bud

Sonic hedgehog (Shh) is essential for limb development, and the mechanisms that govern the propagation and maintenance of its expression has been well studied; however, the mechanisms that govern the initiation of Shh expression are incomplete. Here we report that ETV2 initiates Shh expression by changing the chromatin status of the developmental limb enhancer, ZRS. Etv2 expression precedes Shh in limb buds, and Etv2 inactivation prevents the opening of limb chromatin, including the ZRS, resulting in an absence of Shh expression. Etv2 overexpression in limb buds causes nucleosomal displacement at the ZRS, ectopic Shh expression, and polydactyly. Areas of nucleosome displacement coincide with ETS binding site clusters. ETV2 also functions as a transcriptional activator of ZRS and is antagonized by ETV4/5 repressors. Known human polydactyl mutations introduce novel ETV2 binding sites in the ZRS, suggesting that ETV2 dosage regulates ZRS activation. These studies identify ETV2 as a pioneer transcription factor (TF) regulating the onset of Shh expression, having both a chromatin regulatory role and a transcriptional activation role.

The pZRS non-coding regulatory mutation resulting in triphalangeal thumb-polysyndactyly syndrome changes the pattern of local interactions

Herein, we report on a large Polish family presenting with a classical triphalangeal thumb-polysyndactyly syndrome (TPT-PS). This rare congenital limb anomaly is generally caused by microduplications encompassing the Sonic Hedgehog (SHH) limb enhancer, termed the zone of polarizing activity (ZPA) regulatory sequence (ZRS). Recently, a pathogenic variant in the pre-ZRS (pZRS), a conserved sequence located near the ZRS, has been described in a TPT-PS Dutch family. We performed targeted ZRS sequencing, array comparative genomic hybridization, and whole-exome sequencing. Next, we sequenced the recently described pZRS region. Finally, we performed a circular chromatin conformation capture-sequencing (4C-seq) assay on skin fibroblasts of one affected family member and control samples to examine potential alterations in the SHH regulatory domain and functionally characterize the identified variant. We found that all affected individuals shared a recently identified pathogenic point mutation in the pZRS region: NC_000007.14:g.156792782C>G (GRCh38/hg38), which is the same as in the Dutch family. The results of 4C-seq experiments revealed increased interactions within the whole SHH regulatory domain (SHH-LMBR1 TAD) in the patient compared to controls. Our study expands the number of TPT-PS families carrying a pathogenic alteration of the pZRS and underlines the importance of routine pZRS sequencing in the genetic diagnostics of patients with TPT-PS or similar phenotypes. The pathogenic mutation causative for TPT-PS in our patient gave rise to increased interactions within the SHH regulatory domain in yet unknown mechanism.

The importance of considering regulatory domains in genome-wide analyses - the nearest gene is often wrong!

The expression of a large number of genes is regulated by regulatory elements that are located far away from their promoters. Identifying which gene is the target of a specific regulatory element or is affected by a non-coding mutation is often accomplished by assigning these regions to the nearest gene in the genome. However, this heuristic ignores key features of genome organisation and gene regulation; in that the genome is partitioned into regulatory domains, which at some loci directly coincide with the span of topologically associated domains (TADs), and that genes are regulated by enhancers located throughout these regions, even across intervening genes. In this review, we examine the results from genome-wide studies using chromosome conformation capture technologies and from those dissecting individual gene regulatory domains, to highlight that the phenomenon of enhancer skipping is pervasive and affects multiple types of genes. We discuss how simply assigning a genomic region of interest to its nearest gene is problematic and often leads to incorrect predictions and highlight that where possible information on both the conservation and topological organisation of the genome should be used to generate better hypotheses.

The article has an associated Future Leader to Watch interview.

Summary: Identifying which gene is the target of an enhancer is often accomplished by assigning it to the nearest gene, here we discuss how this heuristic can lead to incorrect predictions.

Perfect and imperfect views of ultraconserved sequences

Across the human genome, there are nearly 500 'ultraconserved' elements: regions of at least 200 contiguous nucleotides that are perfectly conserved in both the mouse and rat genomes. Remarkably, the majority of these sequences are non-coding, and many can function as enhancers that activate tissue-specific gene expression during embryonic development. From their first description more than 15 years ago, their extreme conservation has both fascinated and perplexed researchers in genomics and evolutionary biology. The intrigue around ultraconserved elements only grew with the observation that they are dispensable for viability. Here, we review recent progress towards understanding the general importance and the specific functions of ultraconserved sequences in mammalian development and human disease and discuss possible explanations for their extreme conservation.

Modeling uniquely human gene regulatory function via targeted humanization of the mouse genome

The evolution of uniquely human traits likely entailed changes in developmental gene regulation. Human Accelerated Regions (HARs), which include transcriptional enhancers harboring a significant excess of human-specific sequence changes, are leading candidates for driving gene regulatory modifications in human development. However, insight into whether HARs alter the level, distribution, and timing of endogenous gene expression remains limited. We examined the role of the HAR HACNS1 (HAR2) in human evolution by interrogating its molecular functions in a genetically humanized mouse model. We find that HACNS1 maintains its human-specific enhancer activity in the mouse embryo and modifies expression of Gbx2, which encodes a transcription factor, during limb development. Using single-cell RNA-sequencing, we demonstrate that Gbx2 is upregulated in the limb chondrogenic mesenchyme of HACNS1 homozygous embryos, supporting that HACNS1 alters gene expression in cell types involved in skeletal patterning. Our findings illustrate that humanized mouse models provide mechanistic insight into how HARs modified gene expression in human evolution.

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.

Dynamic Interactions of Transcription Factors and Enhancer Reprogramming in Cancer Progression

While improved tumor treatment has significantly reduced the overall mortality rates, invasive progression including recurrence, therapy resistance and metastasis contributes to the majority of deaths caused by cancer. Enhancers are essential distal DNA regulatory elements that control temporal- or spatial-specific gene expression patterns during development and other biological processes. Genome-wide sequencing has revealed frequent alterations of enhancers in cancers and reprogramming of distal enhancers has emerged as one of the important features for tumors. In this review, we will discuss tumor progression-associated enhancer dynamics, its transcription factor (TF) drivers and how enhancer reprogramming modulates gene expression during cancer invasive progression. Additionally, we will explore recent advancements in contemporary technology including single-cell sequencing, spatial transcriptomics and CUT&RUN, which have permitted integrated studies of enhancer reprogramming in vivo. Given the essential roles of enhancer dynamics and its drivers in controlling cancer progression and treatment outcome, understanding these changes will be paramount in mitigating invasive events and discovering novel therapeutic targets.

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.

Amplicon genome fishing (AGF): a rapid and efficient method for sequencing target cis-regulatory regions in nonmodel organisms

Cis-regulatory sequences play a crucial role in regulating gene expression and are evolutionary hot spots that drive phenotypic divergence among organisms. Sequencing some cis-regulatory regions of interest in many different species is common in comparative genetic studies. For nonmodel organisms lacking genomic data, genome walking is often the preferred method for this type of application. However, applying genome walking will be laborious and time-consuming when the number of cis-regulatory regions and species to be analyzed is large. In this study, we propose a novel method called amplicon genome fishing (AGF), which can isolate and sequence cis-regulatory regions of interest for any organism. The main idea of the AGF method is to use fragments amplified from the target cis-regulatory regions as enrichment baits to capture and sequence the whole target cis-regulatory regions from genomic library pools. Unlike genome walking, the AGF method is based on hybridization capture and high-throughput sequencing, which makes this method rapid and efficient for projects where some cis-regulatory regions have to be sequenced for many species. We used human amplicons as capture baits and successfully sequenced five target enhancer regions of Homo sapiens, Mus musculus, Gallus gallus, and Xenopus tropicalis, proving the feasibility and repeatability of AGF. To show the utility of the AGF method in real studies, we used it to sequence the ZRS enhancer, a cis-regulatory region associated with the limb loss of snakes, for twenty-three vertebrate species (includes many limbless species never sequenced before). The newly obtained ZRS sequences provide new perspectives into the relationship between the ZRS enhancer's evolution and limb loss in major tetrapod lineages.

Apical ectodermal ridge regulates three principal axes of the developing limb

Understanding limb development not only gives insights into the outgrowth and differentiation of the limb, but also has clinical relevance. Limb development begins with two paired limb buds (forelimb and hindlimb buds), which are initially undifferentiated mesenchymal cells tipped with a thickening of the ectoderm, termed the apical ectodermal ridge (AER). As a transitional embryonic structure, the AER undergoes four stages and contributes to multiple axes of limb development through the coordination of signalling centres, feedback loops, and other cell activities by secretory signalling and the activation of gene expression. Within the scope of proximodistal patterning, it is understood that while fibroblast growth factors (FGFs) function sequentially over time as primary components of the AER signalling process, there is still no consensus on models that would explain proximodistal patterning itself. In anteroposterior patterning, the AER has a dual-direction regulation by which it promotes the sonic hedgehog (Shh) gene expression in the zone of polarizing activity (ZPA) for proliferation, and inhibits Shh expression in the anterior mesenchyme. In dorsoventral patterning, the AER activates Engrailed-1 (En1) expression, and thus represses Wnt family member 7a (Wnt7a) expression in the ventral ectoderm by the expression of Fgfs, Sp6/8, and bone morphogenetic protein (Bmp) genes. The AER also plays a vital role in shaping the individual digits, since levels of Fgf4/8 and Bmps expressed in the AER affect digit patterning by controlling apoptosis. In summary, the knowledge of crosstalk within AER among the three main axes is essential to understand limb growth and pattern formation, as the development of its areas proceeds simultaneously.

Establishing the pattern of the vertebrate limb

The vertebrate limb continues to serve as an influential model of growth, morphogenesis and pattern formation. With this Review, we aim to give an up-to-date picture of how a population of undifferentiated cells develops into the complex pattern of the limb. Focussing largely on mouse and chick studies, we concentrate on the positioning of the limbs, the formation of the limb bud, the establishment of the principal limb axes, the specification of pattern, the integration of pattern formation with growth and the determination of digit number. We also discuss the important, but little understood, topic of how gene expression is interpreted into morphology.

Gene regulatory landscape of the sonic hedgehog locus in embryonic development

The organs of vertebrate species display a wide variety of morphology. A remaining challenge in evolutionary developmental biology is to elucidate how vertebrate lineages acquire distinct morphological features. Developmental programs are driven by spatiotemporal regulation of gene expression controlled by hundreds of thousands of cis-regulatory elements. Changes in the regulatory elements caused by the introduction of genetic variants can confer regulatory innovation that may underlie morphological novelties. Recent advances in sequencing technology have revealed a number of potential regulatory variants that can alter gene expression patterns. However, a limited number of studies demonstrate causal dependence between genetic and morphological changes. Regulation of Shh expression is a good model to understand how multiple regulatory elements organize tissue-specific gene expression patterns. This model also provides insights into how evolution of molecular traits, such as gene regulatory networks, lead to phenotypic novelty.

Novel IQCE variations confirm its role in postaxial polydactyly and cause ciliary defect phenotype in zebrafish

Polydactyly is one of the most frequent inherited defects of the limbs characterized by supernumerary digits and high-genetic heterogeneity. Among the many genes involved, either in isolated or syndromic forms, eight have been implicated in postaxial polydactyly (PAP). Among those, IQCE has been recently identified in a single consanguineous family. Using whole-exome sequencing in patients with uncharacterized ciliopathies, including PAP, we identified three families with biallelic pathogenic variations in IQCE. Interestingly, the c.895_904del (p.Val301Serfs*8) was found in all families without sharing a common haplotype, suggesting a recurrent mechanism. Moreover, in two families, the systemic phenotype could be explained by additional pathogenic variants in known genes (TULP1, ATP6V1B1). RNA expression analysis on patients' fibroblasts confirms that the dysfunction of IQCE leads to the dysregulation of genes associated with the hedgehog-signaling pathway, and zebrafish experiments demonstrate a full spectrum of phenotypes linked to defective cilia: Body curvature, kidney cysts, left-right asymmetry, misdirected cilia in the pronephric duct, and retinal defects. In conclusion, we identified three additional families confirming IQCE as a nonsyndromic PAP gene. Our data emphasize the importance of taking into account the complete set of variations of each individual, as each clinical presentation could finally be explained by multiple genes.

A novel ZRS variant causes preaxial polydactyly type I by increased sonic hedgehog expression in the developing limb bud

Purpose

Preaxial polydactyly (PPD) is a common congenital hand malformation classified into four subtypes (PPD I–IV). Variants in the zone of polarizing activity regulatory sequence (ZRS) within intron 5 of the LMBR1 gene are linked to most PPD types. However, the genes responsible for PPD I and the underlying mechanisms are unknown.

Methods

A rare large four-generation family with isolated PPD I was subjected to genome-wide genotyping and sequence analysis. In vitro and in vivo functional studies were performed in Caco-2 cells, 293T cells, and a knockin transgenic mouse model.

Results

A novel g.101779T>A (reference sequence: NG_009240.2; position 446 of the ZRS) variant segregates with all PPD I–affected individuals. The knockin mouse with this ZRS variant exhibited PPD I phenotype accompanying ectopic and excess expression of Shh. We confirmed that HnRNP K can bind the ZRS and SHH promoters. The ZRS mutant enhanced the binding affinity for HnRNP K and upregulated SHH expression.

Conclusion

Our results identify the first PPD I disease-causing variant. The variant leading to PPD I may be associated with enhancing SHH expression mediated by HnRNP K. This study adds to the ZRS-associated syndromes classification system for PPD and clarifies the underlying molecular mechanisms.

Genome Sequencing of the Japanese Eel (Anguilla japonica) for Comparative Genomic Studies on tbx4 and a tbx4 Gene Cluster in Teleost Fishes

Limbs originated from paired fish fins are an important innovation in Gnathostomata. Many studies have focused on limb development-related genes, of which the T-box transcription factor 4 gene (tbx4) has been considered as one of the most essential factors in the regulation of the hindlimb development. We previously confirmed pelvic fin loss in tbx4-knockout zebrafish. Here, we report a high-quality genome assembly of the Japanese eel (Anguilla japonica), which is an economically important fish without pelvic fins. The assembled genome is 1.13 Gb in size, with a scaffold N50 of 1.03 Mb. In addition, we collected 24 tbx4 sequences from 22 teleost fishes to explore the correlation between tbx4 and pelvic fin evolution. However, we observed complete exon structures of tbx4 in several pelvic-fin-loss species such as Ocean sunfish (Mola mola) and ricefield eel (Monopterus albus). More interestingly, an inversion of a special tbx4 gene cluster (brip1-tbx4-tbx2b- bcas3) occurred twice independently, which coincides with the presence of fin spines. A nonsynonymous mutation (M82L) was identified in the nuclear localization sequence (NLS) of the Japanese eel tbx4. We also examined variation and loss of hindlimb enhancer B (HLEB), which may account for pelvic fin loss in Tetraodontidae and Diodontidae. In summary, we generated a genome assembly of the Japanese eel, which provides a valuable genomic resource to study the evolution of fish tbx4 and helps elucidate the mechanism of pelvic fin loss in teleost fishes. Our comparative genomic studies, revealed for the first time a potential correlation between the tbx4 gene cluster and the evolutionary development of toxic fin spines. Because fin spines in teleosts are usually venoms, this tbx4 gene cluster may facilitate the genetic engineering of toxin-related marine drugs.

Structural variation in the 3D genome

Structural and quantitative chromosomal rearrangements, collectively referred to as structural variation (SV), contribute to a large extent to the genetic diversity of the human genome and thus are of high relevance for cancer genetics, rare diseases and evolutionary genetics. Recent studies have shown that SVs can not only affect gene dosage but also modulate basic mechanisms of gene regulation. SVs can alter the copy number of regulatory elements or modify the 3D genome by disrupting higher-order chromatin organization such as topologically associating domains. As a result of these position effects, SVs can influence the expression of genes distant from the SV breakpoints, thereby causing disease. The impact of SVs on the 3D genome and on gene expression regulation has to be considered when interpreting the pathogenic potential of these variant types.

Pigeon foot feathering reveals conserved limb identity networks

The tetrapod limb is a stunning example of evolutionary diversity, with dramatic variation not only among distantly related species, but also between the serially homologous forelimbs (FLs) and hindlimbs (HLs) within species. Despite this variation, highly conserved genetic and developmental programs underlie limb development and identity in all tetrapods, raising the question of how limb diversification is generated from a conserved toolkit. In some breeds of domestic pigeon, shifts in the expression of two conserved limb identity transcription factors, PITX1 and TBX5, are associated with the formation of feathered HLs with partial FL identity. To determine how modulation of PITX1 and TBX5 expression affects downstream gene expression, we compared the transcriptomes of embryonic limb buds from pigeons with scaled and feathered HLs. We identified a set of differentially expressed genes enriched for genes encoding transcription factors, extracellular matrix proteins, and components of developmental signaling pathways with important roles in limb development. A subset of the genes that distinguish scaled and feathered HLs are also differentially expressed between FL and scaled HL buds in pigeons, pinpointing a set of gene expression changes downstream of PITX1 and TBX5 in the partial transformation from HL to FL identity. We extended our analyses by comparing pigeon limb bud transcriptomes to chicken, anole lizard, and mammalian datasets to identify deeply conserved PITX1- and TBX5-responsive components of the limb identity program. Our analyses reveal a suite of predominantly low-level gene expression changes that are conserved across amniotes to regulate the identity of morphologically distinct limbs.

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.

11q24.2q24.3 microdeletion in two families presenting features of Jacobsen syndrome, without intellectual disability: Role of FLI1, ETS1, and SENCR long noncoding RNA

This report presents two families with interstitial 11q24.2q24.3 deletion, associated with malformations, hematologic features, and typical facial dysmorphism, observed in Jacobsen syndrome (JS), except for intellectual disability (ID). The smallest 700 Kb deletion contains only two genes: FLI1 and ETS1, and a long noncoding RNA, SENCR, narrowing the minimal critical region for some features of JS. Consistent with recent literature, it adds supplemental data to confirm the crucial role of FLI1 and ETS1 in JS, namely FLI1 in thrombocytopenia and ETS1 in cardiopathy and immune deficiency. It also supports that combined ETS1 and FLI1 haploinsufficiency explains dysmorphic features, notably ears, and nose anomalies. Moreover, it raises the possibility that SENCR, a long noncoding RNA, could be responsible for limb defects, because of its early role in endothelial cell commitment and function. Considering ID and autism spectrum disorder, which are some of the main features of JS, a participation of ETS1, FLI1, or SENCR cannot be excluded. But, considering the normal neurodevelopment of our patients, their role would be either minor or with an important variability in penetrance. Furthermore, according to literature, ARHGAP32 and KIRREL3 seem to be the strongest candidate genes in the 11q24 region for other Jacobsen patients.

LHX2 Mediates the FGF-to-SHH Regulatory Loop during Limb Development

During limb development, fibroblast growth factors (Fgfs) govern proximal⁻distal outgrowth and patterning. FGFs also synchronize developmental patterning between the proximal⁻distal and anterior⁻posterior axes by maintaining Sonic hedgehog (Shh) expression in cells of the zone of polarizing activity (ZPA) in the distal posterior mesoderm. Shh, in turn, maintains Fgfs in the apical ectodermal ridge (AER) that caps the distal tip of the limb bud. Crosstalk between Fgf and Shh signaling is critical for patterned limb development, but the mechanisms underlying this feedback loop are not well-characterized. Implantation of Fgf beads in the proximal posterior limb bud can maintain SHH expression in the former ZPA domain (evident 3 h after application), while prolonged exposure (24 h) can induce SHH outside of this domain. Although temporally and spatially disparate, comparative analysis of transcriptome data from these different populations accentuated genes involved in SHH regulation. Comparative analysis identified 25 candidates common to both treatments, with eight linked to SHH expression or function. Furthermore, we demonstrated that LHX2, a LIM Homeodomain transcription factor, is an intermediate in the FGF-mediated regulation of SHH. Our data suggest that LHX2 acts as a competency factor maintaining distal posterior SHH expression subjacent to the AER.

Genome-wide association study identifies nine novel loci for 2D:4D finger ratio, a putative retrospective biomarker of testosterone exposure in utero

The ratio of the length of the index finger to that of the ring finger (2D:4D) is sexually dimorphic and is commonly used as a non-invasive biomarker of prenatal androgen exposure. Most association studies of 2D:4D ratio with a diverse range of sex-specific traits have typically involved small sample sizes and have been difficult to replicate, raising questions around the utility and precise meaning of the measure. In the largest genome-wide association meta-analysis of 2D:4D ratio to date (N = 15 661, with replication N = 75 821), we identified 11 loci (9 novel) explaining 3.8% of the variance in mean 2D:4D ratio. We also found weak evidence for association (β = 0.06; P = 0.02) between 2D:4D ratio and sensitivity to testosterone [length of the CAG microsatellite repeat in the androgen receptor (AR) gene] in females only. Furthermore, genetic variants associated with (adult) testosterone levels and/or sex hormone-binding globulin were not associated with 2D:4D ratio in our sample. Although we were unable to find strong evidence from our genetic study to support the hypothesis that 2D:4D ratio is a direct biomarker of prenatal exposure to androgens in healthy individuals, our findings do not explicitly exclude this possibility, and pathways involving testosterone may become apparent as the size of the discovery sample increases further. Our findings provide new insight into the underlying biology shaping 2D:4D variation in the general population.

Genome-wide association study identifies nine novel loci for 2D:4D finger ratio, a putative retrospective biomarker of testosterone exposure in utero

The ratio of the length of the index finger to that of the ring finger (2D:4D) is sexually dimorphic and is commonly used as a non-invasive biomarker of prenatal androgen exposure. Most association studies of 2D:4D ratio with a diverse range of sex-specific traits have typically involved small sample sizes and have been difficult to replicate, raising questions around the utility and precise meaning of the measure. In the largest genome-wide association meta-analysis of 2D:4D ratio to date (N = 15 661, with replication N = 75 821), we identified 11 loci (9 novel) explaining 3.8% of the variance in mean 2D:4D ratio. We also found weak evidence for association (β = 0.06; P = 0.02) between 2D:4D ratio and sensitivity to testosterone [length of the CAG microsatellite repeat in the androgen receptor (AR) gene] in females only. Furthermore, genetic variants associated with (adult) testosterone levels and/or sex hormone-binding globulin were not associated with 2D:4D ratio in our sample. Although we were unable to find strong evidence from our genetic study to support the hypothesis that 2D:4D ratio is a direct biomarker of prenatal exposure to androgens in healthy individuals, our findings do not explicitly exclude this possibility, and pathways involving testosterone may become apparent as the size of the discovery sample increases further. Our findings provide new insight into the underlying biology shaping 2D:4D variation in the general population.

The Conserved Sonic Hedgehog Limb Enhancer Consists of Discrete Functional Elements that Regulate Precise Spatial Expression

Expression of sonic hedgehog (Shh) in the limb bud is regulated by an enhancer called the zone of polarizing activity regulatory sequence (ZRS), which, in evolution, belongs to an ancient group of highly conserved cis regulators found in all classes of vertebrates. Here, we examined the endogenous ZRS in mice, using genome editing to establish the relationship between enhancer composition and embryonic phenotype. We show that enhancer activity is a consolidation of distinct activity domains. Spatial restriction of Shh expression is mediated by a discrete repressor module, whereas levels of gene expression are controlled by large overlapping domains containing varying numbers of HOXD binding sites. The number of HOXD binding sites regulates expression levels incrementally. Substantial portions of conserved sequence are dispensable, indicating the presence of sequence redundancy. We propose a collective model for enhancer activity in which function is an integration of discrete expression activities and redundant components that drive robust expression.

•

The ancient vertebrate enhancer, the ZRS, shows sequence plasticity

•

Discrete regulatory activities are assigned to specific sites in the enhancer

•

The number of HOXD binding sites determines the level of Shh expression

•

Robust expression is a collective of regulatory and redundant information

Noncanonical Hox, Etv4, and Gli3 gene activities give insight into unique limb patterning in salamanders

Limb development in salamanders is unique among tetrapods in significant ways. Not only can salamanders regenerate lost limbs repeatedly and throughout their lives, but also the preaxial zeugopodial element and digits form before the postaxial ones and, hence, with a reversed polarity compared to all other tetrapods. Moreover, in salamanders with free-swimming larval stages, as exemplified by the axolotl (Ambystoma mexicanum), each digit buds independently, instead of undergoing a paddle stage. Here, we report gene expression patterns of Hoxa and d clusters, and other crucial transcription factors during axolotl limb development. During early phases of limb development, expression patterns are mostly similar to those reported for amniotes and frogs. Likewise, Hoxd and Shh regulatory landscapes are largely conserved. However, during late digit-budding phases, remarkable differences are present: (i) the Hoxd13 expression domain excludes developing digits I and IV, (ii) we expand upon previous observation that Hoxa11 expression, which traditionally marks the zeugopodium, extends distally into the developing digits, and (iii) Gli3 and Etv4 show prolonged expression in developing digits. Our findings identify derived patterns in the expression of key transcription factors during late phases of salamander limb development, and provide the basis for a better understanding of the unique patterning of salamander limbs.

A point mutation in the pre-ZRS disrupts sonic hedgehog expression in the limb bud and results in triphalangeal thumb-polysyndactyly syndrome

PURPOSE

The zone of polarizing activity regulatory sequence (ZRS) is an enhancer that regulates sonic hedgehog during embryonic limb development. Recently, mutations in a noncoding evolutionary conserved sequence 500 bp upstream of the ZRS, termed the pre-ZRS (pZRS), have been associated with polydactyly in dogs and humans. Here, we report the first case of triphalangeal thumb-polysyndactyly syndrome (TPT-PS) to be associated with mutations in this region and show via mouse enhancer assays how this mutation leads to ectopic expression throughout the developing limb bud.

METHODS

We used linkage analysis, whole-exome sequencing, Sanger sequencing, fluorescence in situ hybridization, multiplex ligation-dependent probe amplification, single-nucleotide polymorphism array, and a mouse transgenic enhancer assay.

RESULTS

Ten members of a TPT-PS family were included in this study. The mutation was linked to chromosome 7q36 (LOD score 3.0). No aberrations in the ZRS could be identified. A point mutation in the pZRS (chr7:156585476G>C; GRCh37/hg19) was detected in all affected family members. Functional characterization using a mouse transgenic enhancer essay showed extended ectopic expression dispersed throughout the entire limb bud (E11.5).

CONCLUSION

Our work describes the first mutation in the pZRS to be associated with TPT-PS and provides functional evidence that this mutation leads to ectopic expression of this enhancer within the developing limb.

Direct Interactions Between Gli3, Wnt8b, and Fgfs Underlie Patterning of the Dorsal Telencephalon

A key step in the development of the cerebral cortex is a patterning process, which subdivides the telencephalon into several molecularly distinct domains and is critical for cortical arealization. This process is dependent on a complex network of interactions between signaling molecules of the Fgf and Wnt gene families and the Gli3 transcription factor gene, but a better knowledge of the molecular basis of the interplay between these factors is required to gain a deeper understanding of the genetic circuitry underlying telencephalic patterning. Using DNA-binding and reporter gene assays, we here investigate the possibility that Gli3 and these signaling molecules interact by directly regulating each other's expression. We show that Fgf signaling is required for Wnt8b enhancer activity in the cortical hem, whereas Wnt/β-catenin signaling represses Fgf17 forebrain enhancer activity. In contrast, Fgf and Wnt/β-catenin signaling cooperate to regulate Gli3 expression. Taken together, these findings indicate that mutual interactions between Gli3, Wnt8b, and Fgf17 are crucial elements of the balance between these factors thereby conferring robustness to the patterning process. Hence, our study provides a framework for understanding the genetic circuitry underlying telencephalic patterning and how defects in this process can affect the formation of cortical areas.

Zone of Polarizing Activity Regulatory Sequence Mutations/Duplications with Preaxial Polydactyly and Longitudinal Preaxial Ray Deficiency in the Phenotype: A Review of Human Cases, Animal Models, and Insights Regarding the Pathogenesis

Clinicians and scientists interested in developmental biology have viewed preaxial polydactyly (PPD) and longitudinal preaxial ray deficiency (LPAD) as two different entities. Point mutations and duplications in the zone of polarizing activity regulatory sequence (ZRS) are associated with anterior ectopic expression of Sonic Hedgehog (SHH) in the limb bud and usually result in a PPD phenotype. However, some of these mutations/duplications also have LPAD in the phenotype. This unusual PPD-LPAD association in ZRS mutations/duplications has not been specifically reviewed in the literature. The author reviews this unusual entity and gives insights regarding its pathogenesis.

Evolution: Enhanced Footing for Snake Limb Development

Two groups have studied the loss of limbs in snake evolution by focusing on a long-distance cis-acting enhancer of Sonic Hedgehog. They find a progressive degeneration of binding sites for key transcription factors, mirroring the progressive limblessness occurring in these reptiles.

Fibroblast growth factors (FGFs) prime the limb specific Shh enhancer for chromatin changes that balance histone acetylation mediated by E26 transformation-specific (ETS) factors

Sonic hedgehog (Shh) expression in the limb bud organizing centre called the zone of polarizing activity is regulated by the ZRS enhancer. Here, we examine in mouse and in a mouse limb-derived cell line the dynamic events that activate and restrict the spatial activity of the ZRS. Fibroblast growth factor (FGF) signalling in the distal limb primes the ZRS at early embryonic stages maintaining a poised, but inactive state broadly across the distal limb mesenchyme. The E26 transformation-specific transcription factor, ETV4, which is induced by FGF signalling and acts as a repressor of ZRS activity, interacts with the histone deacetylase HDAC2 and ensures that the poised ZRS remains transcriptionally inactive. Conversely, GABPα, an activator of the ZRS, recruits p300, which is associated with histone acetylation (H3K27ac) indicative of an active enhancer. Hence, the primed but inactive state of the ZRS is induced by FGF signalling and in combination with balanced histone modification events establishes the restricted, active enhancer responsible for patterning the limb bud during development.

As an animal embryo develops, specific genes need to be switched on and off at the right time and place to ensure that the embryo’s tissues and organs form properly. Proteins called transcription factors control the activity of individual genes by binding to regions of DNA known as enhancers. Changes in the way DNA is packaged inside cells can affect the ability of transcription factors to access the enhancers, and therefore also influence when particular genes are switched on or off.

Sonic hedgehog (or Shh for short) is a gene that helps to control various aspects of development including the formation of the limbs and brain. The limb forms from a structure in the embryo referred to as the limb bud. An enhancer called ZRS regulates the precise position within the limb bud where the Shh gene is active in a region designated as the “zone of polarizing activity”. Yet, it was not known how the enhancer is controlled to ensure this pattern is achieved. Peluso et al. investigated the events that lead to ZRS becoming active in mice embryos.

The experiments show that the ZRS enhancer exists in three different states in cells across the limb bud: poised, active and inactive. The enhancer is poised in a broad region of the limb bud in cells that are potentially able to switch on the Shh gene. Proteins called fibroblast growth factors drive the enhancer to enter this poised state by altering the way the DNA containing the enhancer is packaged in the cell. Specific transcription factors are able to bind to the poised enhancer and it is the balance between these different transcription factors that activates the enhancer in the zone of polarizing activity. Furthermore in the region of the limb bud where the fibroblast growth factors are not present the ZRS is inactive.

These findings show that fibroblast growth factors, in combination with other changes to the ZRS enhancer, restrict the area in which the enhancer is active to a particular region of the limb bud. Differences in enhancer elements are known to underlie a range of inherited characteristics and may influence whether an individual develops many common diseases. In the future, investigating how cells control the activity of enhancers may provide clues to identifying new targets for drugs to treat some of these diseases.