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

Evidence for Fgf and Wnt regulation of Lhx2 during limb development via two limb-specific Lhx2-associated cis-regulatory modules

Introduction

In vertebrate limb morphogenesis, wingless-related integration site (Wnt) proteins and fibroblast growth factors (Fgfs) secreted from the apical ectodermal ridge (AER) coordinate proximodistal outgrowth. Fgfs also sustain sonic hedgehog (Shh) in the zone of polarizing activity (ZPA). Shh directs anteroposterior patterning and expansion and regulates AER-Fgfs, establishing a positive regulatory feedback loop that is vital in sustaining limb outgrowth. The transcription factor LIM homeodomain 2 (Lhx2) is expressed in the distal mesoderm and coordinates AER and ZPA signals that control cellular proliferation, differentiation, and shaping of the developing limb. Yet how Lhx2 is transcriptionally regulated to support such functions has only been partially characterized.

Methods/Results

We have identified two limb-specific cis-regulatory modules (CRMs) active within the Lhx2 expression domain in the limb. Chromatin conformation analysis of the Lhx2 locus in mouse embryonic limb bud cells predicted CRMs-Lhx2 promoter interactions. Single-cell RNA-sequencing analysis of limb bud cells revealed co-expression of several Fgf-related Ets and Wnt-related Tcf/Lef transcripts in Lhx2-expressing cells. Additionally, disruption of Ets and Tcf/Lef binding sites resulted in loss of reporter-driven CRM activity. Finally, binding of β-catenin to both Lhx2-associated CRMs supports the associated binding of Tcf/Lef transcription factors.

Discussion

These results suggest a role for Ets and Tcf/Lef transcription factors in the regulation of Lhx2 expression through these limb-specific Lhx2-associated CRMs. Moreover, these CRMs provide a mechanism for Fgf and Wnt signaling to localize and maintain distal Lhx2 expression during vertebrate limb development.

Anterior-posterior constraint on Hedgehog signaling by hhip in teleost fin elaboration

Pectoral fins, the anterior paired fins in fish, have enhanced maneuvering abilities due to morphological changes. Teleosts have fewer radial bones in their pectoral fins than basal species, resulting in more-elaborate fins. The mechanism behind this radial constraint change in teleosts is unclear. Here, we found that mutations in hhip, which encodes an antagonist of Hedgehog signaling, led to an increase in radial bones in a localized region. Expression of the Shh genes, encoding ligands of Hedgehog signaling, coincided with notable hhip expression specifically during early development. We suggest that a negative feedback effect of Hedgehog signaling by hhip regulates the constraint of the pectoral fin in zebrafish. Additionally, re-analysis of hhip-related gene expression data in zebrafish and basal species revealed that the notable hhip expression during early development is a characteristic of zebrafish that is not observed in basal species. Region-specific expression of Hox13 genes in the zebrafish pectoral fin indicated that the median region, analogous to the region with abundant radials in basal species, is expanded in hhip-/- zebrafish. These data underscore potential morphological evolution through constrained diversity.

Development of a ZRS Reporter System for the Newt (Cynops pyrrhogaster) During Terrestrial Limb Regeneration

BACKGROUND

Newts, a type of urodele amphibian, offer remarkable insights into regenerative medicine due to their extraordinary tissue regeneration capabilities-a challenging feat in humans. During limb regeneration of adult newts, fascinating cellular and molecular processes are revealed, including scarless healing, de-differentiation of mature cells, and regeneration of limbs and digits. Sonic hedgehog (Shh), crucial for vertebrate limb development, is regulated by the zone of polarizing activity regulatory sequence (ZRS) in the limb bud zone of polarizing activity (ZPA). The metamorphosed (terrestrial) newt can reactivate Shh during regeneration, facilitating proper limb patterning. Cell types capable of regulating the ZRS in metamorphosed newts remain unknown. The identification of such cell types provides invaluable insight into novel regenerative mechanisms.

OBJECTIVE

In this study, we developed the first newt ZRS reporter.

METHODS

We isolated and characterized the newt ZRS enhancer (nZRS), identifying conserved DNA binding sites. Several binding sites with medical relevance were conserved in the newt ZRS. In functional analysis, we developed a system composed of a transgenic nZRS reporter newt and a new newt anti-Shh antibody, which allowed Shh monitoring during limb regeneration.

RESULTS

We identified a group of Schwann cells capable of ZRS reporter and Shh protein expression during terrestrial limb regeneration.

CONCLUSIONS

This system provides a valuable in vivo approach for future genetic studies of patterning during limb regeneration.

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.

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.

Emergence of a left-right symmetric body plan in vertebrate embryos

External bilateral symmetry is a prevalent feature in vertebrates, which emerges during early embryonic development. To begin with, vertebrate embryos are largely radially symmetric before transitioning to bilaterally symmetry, after which, morphogenesis of various bilateral tissues (e.g somites, otic vesicle, limb bud), and structures (e.g palate, jaw) ensue. While a significant amount of work has probed the mechanisms behind symmetry breaking in the left-right axis leading to asymmetric positioning of internal organs, little is known about how bilateral tissues emerge at the same time with the same shape and size and at the same position on the two sides of the embryo. By discussing emergence of symmetry in many bilateral tissues and structures across vertebrate model systems, we highlight that understanding symmetry establishment is largely an open field, which will provide deep insights into fundamental problems in developmental biology for decades to come.

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.

Rapid and Quantitative Functional Interrogation of Human Enhancer Variant Activity in Live Mice

Functional analysis of non-coding variants associated with human congenital disorders remains challenging due to the lack of efficient in vivo models. Here we introduce dual-enSERT, a robust Cas9-based two-color fluorescent reporter system which enables rapid, quantitative comparison of enhancer allele activities in live mice of any genetic background. We use this new technology to examine and measure the gain- and loss-of-function effects of enhancer variants linked to limb polydactyly, autism, and craniofacial malformation. By combining dual-enSERT with single-cell transcriptomics, we characterize variant enhancer alleles at cellular resolution, thereby implicating candidate molecular pathways in pathogenic enhancer misregulation. We further show that independent, polydactyly-linked enhancer variants lead to ectopic expression in the same cell populations, indicating shared genetic mechanisms underlying non-coding variant pathogenesis. Finally, we streamline dual-enSERT for analysis in F0 animals by placing both reporters on the same transgene separated by a synthetic insulator. Dual-enSERT allows researchers to go from identifying candidate enhancer variants to analysis of comparative enhancer activity in live embryos in under two weeks.

What is an enhancer?

Tight control of the transcription process is essential for the correct spatial and temporal gene expression pattern during development and in homeostasis. Enhancers are at the core of correct transcriptional activation. The original definition of an enhancer is straightforward: a DNA sequence that activates transcription independent of orientation and direction. Dissection of numerous enhancer loci has shown that many enhancer-like elements might not conform to the original definition, suggesting that enhancers and enhancer-like elements might use multiple different mechanisms to contribute to transcriptional activation. Here, we review methodologies to identify enhancers and enhancer-like elements and discuss pitfalls and consequences for our understanding of transcriptional regulation.

Developmental regulation of conserved non-coding element evolution provides insights into limb loss in squamates

Limb loss shows recurrent phenotypic evolution across squamate lineages. Here, based on three de novo-assembled genomes of limbless lizards from different lineages, we showed that divergence of conserved non-coding elements (CNEs) played an important role in limb development. These CNEs were associated with genes required for limb initiation and outgrowth, and with regulatory signals in the early stage of limb development. Importantly, we identified the extensive existence of insertions and deletions (InDels) in the CNEs, with the numbers ranging from 111 to 756. Most of these CNEs with InDels were lineage-specific in the limbless squamates. Nearby genes of these InDel CNEs were important to early limb formation, such as Tbx4, Fgf10, and Gli3. Based on functional experiments, we found that nucleotide mutations and InDels both affected the regulatory function of the CNEs. Our study provides molecular evidence underlying limb loss in squamate reptiles from a developmental perspective and sheds light on the importance of regulatory element InDels in phenotypic evolution.

Breaking enhancers to gain insights into developmental defects

Despite ground-breaking genetic studies that have identified thousands of risk variants for developmental diseases, how these variants lead to molecular and cellular phenotypes remains a gap in knowledge. Many of these variants are non-coding and occur at enhancers, which orchestrate key regulatory programs during development. The prevailing paradigm is that non-coding variants alter the activity of enhancers, impacting gene expression programs, and ultimately contributing to disease risk. A key obstacle to progress is the systematic functional characterization of non-coding variants at scale, especially since enhancer activity is highly specific to cell type and developmental stage. Here, we review the foundational studies of enhancers in developmental disease and current genomic approaches to functionally characterize developmental enhancers and their variants at scale. In the coming decade, we anticipate systematic enhancer perturbation studies to link non-coding variants to molecular mechanisms, changes in cell state, and disease phenotypes.

Noncoding variants alter GATA2 expression in rhombomere 4 motor neurons and cause dominant hereditary congenital facial paresis

Hereditary congenital facial paresis type 1 (HCFP1) is an autosomal dominant disorder of absent or limited facial movement that maps to chromosome 3q21-q22 and is hypothesized to result from facial branchial motor neuron (FBMN) maldevelopment. In the present study, we report that HCFP1 results from heterozygous duplications within a neuron-specific GATA2 regulatory region that includes two enhancers and one silencer, and from noncoding single-nucleotide variants (SNVs) within the silencer. Some SNVs impair binding of NR2F1 to the silencer in vitro and in vivo and attenuate in vivo enhancer reporter expression in FBMNs. Gata2 and its effector Gata3 are essential for inner-ear efferent neuron (IEE) but not FBMN development. A humanized HCFP1 mouse model extends Gata2 expression, favors the formation of IEEs over FBMNs and is rescued by conditional loss of Gata3. These findings highlight the importance of temporal gene regulation in development and of noncoding variation in rare mendelian disease.

Regulatory de novo mutations underlying intellectual disability

The genetic aetiology of a major fraction of patients with intellectual disability (ID) remains unknown. De novo mutations (DNMs) in protein-coding genes explain up to 40% of cases, but the potential role of regulatory DNMs is still poorly understood. We sequenced 63 whole genomes from 21 ID probands and their unaffected parents. In addition, we analysed 30 previously sequenced genomes from exome-negative ID probands. We found that regulatory DNMs were selectively enriched in fetal brain-specific enhancers as compared with adult brain enhancers. DNM-containing enhancers were associated with genes that show preferential expression in the prefrontal cortex. Furthermore, we identified recurrently mutated enhancer clusters that regulate genes involved in nervous system development (CSMD1, OLFM1, and POU3F3). Most of the DNMs from ID probands showed allele-specific enhancer activity when tested using luciferase assay. Using CRISPR-mediated mutation and editing of epigenomic marks, we show that DNMs at regulatory elements affect the expression of putative target genes. Our results, therefore, provide new evidence to indicate that DNMs in fetal brain-specific enhancers play an essential role in the aetiology of ID.

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.

Gene expression changes during the evolution of the tetrapod limb

Major changes in the vertebrate anatomy have preceded the conquest of land by the members of this taxon, and continuous changes in limb shape and use have occurred during the later radiation of tetrapods. While the main, conserved mechanisms of limb development have been discerned over the past century using a combination of classical embryological and molecular methods, only recent advances made it possible to identify and study the regulatory changes that have contributed to the evolution of the tetrapod appendage. These advances include the expansion of the model repertoire from traditional genetic model species to non-conventional ones, a proliferation of predictive mathematical models that describe gene interactions, an explosion in genomic data and the development of high-throughput methodologies. These revolutionary innovations make it possible to identify specific mutations that are behind specific transitions in limb evolution. Also, as we continue to apply them to more and more extant species, we can expect to gain a fine-grained view of this evolutionary transition that has been so consequential for our species as well.

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.

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.

Cancer Evo-Dev: A Theory of Inflammation-Induced Oncogenesis

Chronic inflammation is a prerequisite for the development of cancers. Here, we present the framework of a novel theory termed as Cancer Evolution-Development (Cancer Evo-Dev) based on the current understanding of inflammation-related carcinogenesis, especially hepatocarcinogenesis induced by chronic infection with hepatitis B virus. The interaction between genetic predispositions and environmental exposures, such as viral infection, maintains chronic non-resolving inflammation. Pollution, metabolic syndrome, physical inactivity, ageing, and adverse psychosocial exposure also increase the risk of cancer via inducing chronic low-grade smoldering inflammation. Under the microenvironment of non-resolving inflammation, pro-inflammatory factors facilitate the generation of somatic mutations and viral mutations by inducing the imbalance between the mutagenic forces such as cytidine deaminases and mutation-correcting forces including uracil-DNA glycosylase. Most cells with somatic mutations and mutated viruses are eliminated in survival competition. Only a small percentage of mutated cells survive, adapt to the hostile environment, retro-differentiate, and function as cancer-initiating cells via altering signaling pathways. These cancer-initiating cells acquire stem-ness, reprogram metabolic patterns, and affect the microenvironment. The carcinogenic process follows the law of "mutation-selection-adaptation". Chronic physical activity reduces the levels of inflammation via upregulating the activity and numbers of NK cells and lymphocytes and lengthening leukocyte telomere; downregulating proinflammatory cytokines including interleukin-6 and senescent lymphocytes especially in aged population. Anti-inflammation medication reduces the occurrence and recurrence of cancers. Targeting cancer stemness signaling pathways might lead to cancer eradication. Cancer Evo-Dev not only helps understand the mechanisms by which inflammation promotes the development of cancers, but also lays the foundation for effective prophylaxis and targeted therapy of various cancers.

Genetic basis for an evolutionary shift from ancestral preaxial to postaxial limb polarity in non-urodele vertebrates

In most tetrapod vertebrates, limb skeletal progenitors condense with postaxial dominance. Posterior elements (such as ulna and fibula) appear prior to their anterior counterparts (radius and tibia), followed by digit-appearance order with continuing postaxial polarity. The only exceptions are urodele amphibians (salamanders), whose limb elements develop with preaxial polarity and who are also notable for their unique ability to regenerate complete limbs as adults. The mechanistic basis for this preaxial dominance has remained an enigma and has even been proposed to relate to the acquisition of novel genes involved in regeneration. However, recent fossil evidence suggests that preaxial polarity represents an ancestral rather than derived state. Here, we report that 5'Hoxd (Hoxd11-d13) gene deletion in mouse is atavistic and uncovers an underlying preaxial polarity in mammalian limb formation. We demonstrate this shift from postaxial to preaxial dominance in mouse results from excess Gli3 repressor (Gli3R) activity due to the loss of 5'Hoxd-Gli3 antagonism and is associated with cell-cycle changes promoting precocious cell-cycle exit in the anterior limb bud. We further show that Gli3 knockdown in axolotl results in a shift to postaxial dominant limb skeleton formation, as well as expanded paddle-shaped limb-bud morphology and ensuing polydactyly. Evolutionary changes in Gli3R activity level, which also played a key role in the fin-to-limb transition, appear to be fundamental to the shift from preaxial to postaxial polarity in formation of the tetrapod limb skeleton.

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.

Open chromatin in grapevine marks candidate CREs and with other chromatin features correlates with gene expression

Vitis vinifera is an economically important crop and a useful model in which to study chromatin dynamics. In contrast to the small and relatively simple genome of Arabidopsis thaliana, grapevine contains a complex genome of 487 Mb that exhibits extensive colonization by transposable elements. We used Hi-C, ChIP-seq and ATAC-seq to measure how chromatin features correlate to the expression of 31 845 grapevine genes. ATAC-seq revealed the presence of more than 16 000 open chromatin regions, of which we characterize nearly 5000 as possible distal enhancer candidates that occur in intergenic space > 2 kb from the nearest transcription start site (TSS). A motif search identified more than 480 transcription factor (TF) binding sites in these regions, with those for TCP family proteins in greatest abundance. These open chromatin regions are typically within 15 kb from their nearest promoter, and a gene ontology analysis indicated that their nearest genes are significantly enriched for TF activity. The presence of a candidate cis-regulatory element (cCRE) > 2 kb upstream of the TSS, location in the active nuclear compartment as determined by Hi-C, and the enrichment of H3K4me3, H3K4me1 and H3K27ac at the gene are correlated with gene expression. Taken together, these results suggest that regions of intergenic open chromatin identified by ATAC-seq can be considered potential candidates for cis-regulatory regions in V. vinifera. Our findings enhance the characterization of a valuable agricultural crop, and help to clarify the understanding of unique plant biology.

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.

Ultraconserved enhancer function does not require perfect sequence conservation

Ultraconserved enhancer sequences show perfect conservation between human and rodent genomes, suggesting that their functions are highly sensitive to mutation. However, current models of enhancer function do not sufficiently explain this extreme evolutionary constraint. We subjected 23 ultraconserved enhancers to different levels of mutagenesis, collectively introducing 1,547 mutations, and examined their activities in transgenic mouse reporter assays. Overall, we find that the regulatory properties of ultraconserved enhancers are robust to mutation. Upon mutagenesis, nearly all (19/23, 83%) still functioned as enhancers at one developmental stage, as did most of those tested again later in development (5/9, 56%). Replacement of endogenous enhancers with mutated alleles in mice corroborated results of transgenic assays, including the functional resilience of ultraconserved enhancers to mutation. Our findings show that the currently known activities of ultraconserved enhancers do not necessarily require the perfect conservation observed in evolution and suggest that additional regulatory or other functions contribute to their sequence constraint.

Super-Enhancers and CTCF in Early Embryonic Cell Fate Decisions

Cell fate decisions are the backbone of many developmental and disease processes. In early mammalian development, precise gene expression changes underly the rapid division of a single cell that leads to the embryo and are critically dependent on autonomous cell changes in gene expression. To understand how these lineage specifications events are mediated, scientists have had to look past protein coding genes to the cis regulatory elements (CREs), including enhancers and insulators, that modulate gene expression. One class of enhancers, termed super-enhancers, is highly active and cell-type specific, implying their critical role in modulating cell-type specific gene expression. Deletion or mutations within these CREs adversely affect gene expression and development and can cause disease. In this mini-review we discuss recent studies describing the potential roles of two CREs, enhancers and binding sites for CTCF, in early mammalian development.

Congenital Malformations in Sea Turtles: Puzzling Interplay between Genes and Environment

Simple Summary

Congenital malformations can lead to embryonic mortality in many species, and sea turtles are no exception. Genetic and/or environmental alterations occur during early development in the embryo, and may produce aberrant phenotypes, many of which are incompatible with life. Causes of malformations are multifactorial; genetic factors may include mutations, chromosomal aberrations, and inbreeding effects, whereas non-genetic factors may include nutrition, hyperthermia, low moisture, radiation, and contamination. It is possible to monitor and control some of these factors (such as temperature and humidity) in nesting beaches, and toxic compounds in feeding areas, which can be transferred to the embryo through their lipophilic properties. In this review, we describe possible causes of different types of malformations observed in sea turtle embryos, as well as some actions that may help reduce embryonic mortality.

Abstract

The completion of embryonic development depends, in part, on the interplay between genetic factors and environmental conditions, and any alteration during development may affect embryonic genetic and epigenetic regulatory pathways leading to congenital malformations, which are mostly incompatible with life. Oviparous reptiles, such as sea turtles, that produce numerous eggs in a clutch that is buried on the beach provide an opportunity to study embryonic mortality associated with malformations that occur at different times during development, or that prevent the hatchling from emerging from the nest. In sea turtles, the presence of congenital malformations frequently leads to mortality. A few years ago, a detailed study was performed on external congenital malformations in three species of sea turtles from the Mexican Pacific and Caribbean coasts, the hawksbill turtle, Eretmochelys imbricata (n = 23,559 eggs), the green turtle, Chelonia mydas (n = 17,690 eggs), and the olive ridley, Lepidochelys olivacea (n = 20,257 eggs), finding 63 types of congenital malformations, of which 38 were new reports. Of the three species, the olive ridley showed a higher incidence of severe anomalies in the craniofacial region (49%), indicating alterations of early developmental pathways; however, several malformations were also observed in the body, including defects in the carapace (45%) and limbs (33%), as well as pigmentation disorders (20%), indicating that deviations occurred during the middle and later stages of development. Although intrinsic factors (i.e., genetic mutations or epigenetic modifications) are difficult to monitor in the field, some environmental factors (such as the incubation temperature, humidity, and probably the status of feeding areas) are, to some extent, less difficult to monitor and/or control. In this review, we describe the aetiology of different malformations observed in sea turtle embryos, and provide some actions that can reduce embryonic mortality.

Dissection of the Fgf8 regulatory landscape by in vivo CRISPR-editing reveals extensive intra- and inter-enhancer redundancy

Developmental genes are often regulated by multiple elements with overlapping activity. Yet, in most cases, the relative function of those elements and their contribution to endogenous gene expression remain poorly characterized. An example of this phenomenon is that distinct sets of enhancers have been proposed to direct Fgf8 in the limb apical ectodermal ridge and the midbrain-hindbrain boundary. Using in vivo CRISPR/Cas9 genome engineering, we functionally dissect this complex regulatory ensemble and demonstrate two distinct regulatory logics. In the apical ectodermal ridge, the control of Fgf8 expression appears distributed between different enhancers. In contrast, we find that in the midbrain-hindbrain boundary, one of the three active enhancers is essential while the other two are dispensable. We further dissect the essential midbrain-hindbrain boundary enhancer to reveal that it is also composed by a mixture of essential and dispensable modules. Cross-species transgenic analysis of this enhancer suggests that its composition may have changed in the vertebrate lineage.

Chromatin topology and the timing of enhancer function at the HoxD locus

The HoxD gene cluster is critical for proper limb formation in tetrapods. In the emerging limb buds, different subgroups of Hoxd genes respond first to a proximal regulatory signal, then to a distal signal that organizes digits. These two regulations are exclusive from one another and emanate from two distinct topologically associating domains (TADs) flanking HoxD, both containing a range of appropriate enhancer sequences. The telomeric TAD (T-DOM) contains several enhancers active in presumptive forearm cells and is divided into two sub-TADs separated by a CTCF-rich boundary, which defines two regulatory submodules. To understand the importance of this particular regulatory topology to control Hoxd gene transcription in time and space, we either deleted or inverted this sub-TAD boundary, eliminated the CTCF binding sites, or inverted the entire T-DOM to exchange the respective positions of the two sub-TADs. The effects of such perturbations on the transcriptional regulation of Hoxd genes illustrate the requirement of this regulatory topology for the precise timing of gene activation. However, the spatial distribution of transcripts was eventually resumed, showing that the presence of enhancer sequences, rather than either their exact topology or a particular chromatin architecture, is the key factor. We also show that the affinity of enhancers to find their natural target genes can overcome the presence of both a strong TAD border and an unfavorable orientation of CTCF sites.

Perspective: Quality Versus Quantity; Is It Important to Assess the Role of Enhancers in Complex Disease from an In Vivo Perspective?

Sequencing of the human genome has permitted the development of genome-wide association studies (GWAS) to analyze the genetics of a number of complex disorders such as depression, anxiety and substance abuse. Thanks to their ability to analyze huge cohort sizes, these studies have successfully identified thousands of loci associated with a broad spectrum of complex diseases. Disconcertingly, the majority of these GWAS hits occur in non-coding regions of the genome, much of which controls the cell-type-specific expression of genes essential to health. In contrast to gene coding sequences, it is a challenge to understand the function of this non-coding regulatory genome using conventional biochemical techniques in cell lines. The current commentary scrutinizes the field of complex genetics from the standpoint of the large-scale whole-genome functional analysis of the promoters and cis-regulatory elements using chromatin markers. We contrast these large scale quantitative techniques against comparative genomics and in vivo analyses including CRISPR/CAS9 genome editing to determine the functional characteristics of these elements and to understand how polymorphic variation and epigenetic changes within these elements might contribute to complex disease and drug response. Most importantly, we suggest that, although the role of chromatin markers will continue to be important in identifying and characterizing enhancers, more emphasis must be placed on their analysis in relevant in-vivo models that take account of the appropriate cell-type-specific roles of these elements. It is hoped that offering these insights might refocus progress in analyzing the data tsunami of non-coding GWAS and whole-genome sequencing "hits" that threatens to overwhelm progress in the field.

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.

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.

Developmentally regulated Shh expression is robust to TAD perturbations

Topologically associating domains (TADs) have been proposed to both guide and constrain enhancer activity. Shh is located within a TAD known to contain all its enhancers. To investigate the importance of chromatin conformation and TAD integrity on developmental gene regulation, we have manipulated the Shh TAD - creating internal deletions, deleting CTCF sites, and deleting and inverting sequences at TAD boundaries. Chromosome conformation capture and fluorescence in situ hybridisation assays were used to investigate the changes in chromatin conformation that result from these manipulations. Our data suggest that these substantial alterations in TAD structure have no readily detectable effect on Shh expression patterns or levels of Shh expression during development - except where enhancers are deleted - and result in no detectable phenotypes. Only in the case of a larger deletion at one TAD boundary could ectopic influence of the Shh limb enhancer be detected on a gene (Mnx1) in the neighbouring TAD. Our data suggests that, contrary to expectations, the developmental regulation of Shh expression is remarkably robust to TAD perturbations.

A multidisciplinary review of triphalangeal thumb

Despite being a rare congenital limb anomaly, triphalangeal thumb is a subject of research in various scientific fields, providing new insights in clinical research and evolutionary biology. The findings of triphalangeal thumb can be predictive for other congenital anomalies as part of an underlying syndrome. Furthermore, triphalangeal thumb is still being used as a model in molecular genetics to study gene regulation by long-range regulatory elements. We present a review that summarizes a number of scientifically relevant topics that involve the triphalangeal thumb phenotype. Future initiatives involving multidisciplinary teams collaborating in the field of triphalangeal thumb research can lead to a better understanding of the pathogenesis and molecular mechanisms of this condition as well as other congenital upper limb anomalies.

FGF, BMP, and RA signaling are sufficient for the induction of complete limb regeneration from non-regenerating wounds on Ambystoma mexicanum limbs

Some organisms, such as the Mexican axolotl, have the capacity to regenerate complicated biological structures throughout their lives. Which molecular pathways are sufficient to induce a complete endogenous regenerative response in injured tissue is an important question that remains unanswered. Using a gain-of-function regeneration assay, known as the Accessory Limb Model (ALM), we and others have begun to identify the molecular underpinnings of the three essential requirements for limb regeneration; wounding, neurotrophic signaling, and the induction of pattern from cells that retain positional memory. We have previously shown that treatment of Mexican axolotls with exogenous retinoic acid (RA) is sufficient to induce the formation of complete limb structures from blastemas that were generated by deviating a nerve bundle into an anterior-located wound site on the limb. Here we show that these ectopic structures are capable of regenerating and inducing new pattern to form when grafted into new anterior-located wounds. We additionally found that the expression of Alx4 decreases, and Shh expression increases in these anterior located blastemas, but not in the mature anterior tissues, supporting the hypothesis that RA treatment posteriorizes blastema tissue. Based on these and previous observations, we used the ALM assay to test the hypothesis that a complete regenerative response can be generated by treating anterior-located superficial limb wounds with a specific combination of growth factors at defined developmental stages. Our data shows that limb wounds that are first treated with a combination of FGF-2, FGF-8, and BMP-2, followed by RA treatment of the resultant mid-bud stage blastema, will result in the generation of limbs with complete proximal/distal and anterior/posterior limb axes. Thus, the minimal signaling requirements from the nerve and a positional disparity are achieved with the application of this specific combination of signaling molecules.

Shh signaling influences the phenotype of Pitx1-/- hindlimbs

Forelimbs (FLs) and hindlimbs (HLs) develop under the instructive and integrated guidance of signaling centers and transcription factor (TF) action. The development of structures specific to each limb type depends on the limb-specific modulation of these integrated components. Pitx1 is a transcription factor gene expressed in HL, absent in FL, and required for HL-specific patterning and development, in particular for formation of anterior HL skeletal elements. Pitx1 achieves this function by direct TF action on the core limb program, which is largely shared between FL and HL. Shh signaling plays a crucial role in anterior-posterior (AP) patterning in both FL and HL. The present work assessed the relationship between Shh signaling and Pitx1 action for AP patterning. We found that reducing the gene dosage of Shh in the context of the Pitx1-/- HL decreases the severity of the Pitx1-/- phenotype, in particular, the loss of anterior limb structures and the shortening of femur length. However, this did not rescue HL-specific patterning features. Thus, Pitx1 action integrates Shh signaling but not for limb-type-specific patterning.

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.

The vertebrate limb: An evolving complex of self-organizing systems

The paired appendages (fins or limbs) of jawed vertebrates contain an endoskeleton consisting of nodules, bars and, in some groups, plates of cartilage, or bone arising from replacement of cartilaginous templates. The generation of the endoskeletal elements occurs by processes involving production and diffusion of morphogens, with, variously, positive and negative feedback circuits, adhesion, and receptor dynamics with similarities to the mechanism for chemical pattern formation proposed by Alan Turing. This review presents a unified interpretation of the evolution and functioning of these mechanisms. Studies are described indicating that protocondensations, compacted mesenchymal cell aggregates that prefigure the appendicular skeleton, arise through the adhesive activity of galectin-1, a matricellular protein with skeletogenic homologs in all jawed vertebrates. In the cartilaginous and lobe-finned fishes (and to a variable extent in ray-finned fishes) it additionally cooperates with an isoform of galectin-8 to constitute a self-organizing network capable of generating arrays of preskeletal nodules, bars and plates. Further, in the tetrapods, a putative galectin-8 control module was acquired that may have enabled proximodistal increase in the number of protocondensations. In parallel to this, other self-organizing networks emerged that acted, via Bmp, Wnt, Sox9 and Runx2, as well as transforming factor-β and fibronectin, to convert protocondensations into skeletal tissues. The progressive appearance and integration of these skeletogenic networks over evolution occurred in the context of an independently evolved system of Hox protein and Shh gradients that interfaced with them to tune the spatial wavelengths and refine the identities of the resulting arrays of elements.