cis-regulatory mutations are a genetic cause of human limb malformations

Cis-regulatory landscapes in the evolution and development of the mammalian skull

Extensive morphological variation found in mammals reflects the wide spectrum of their ecological adaptations. The highest morphological diversity is present in the craniofacial region, where geometry is mainly dictated by the bony skull. Mammalian craniofacial development represents complex multistep processes governed by numerous conserved genes that require precise spatio-temporal control. A central question in contemporary evolutionary biology is how a defined set of conserved genes can orchestrate formation of fundamentally different structures, and therefore how morphological variability arises. In principle, differential gene expression patterns during development are the source of morphological variation. With the emergence of multicellular organisms, precise regulation of gene expression in time and space is attributed to cis-regulatory elements. These elements contribute to higher-order chromatin structure and together with trans-acting factors control transcriptional landscapes that underlie intricate morphogenetic processes. Consequently, divergence in cis-regulation is believed to rewire existing gene regulatory networks and form the core of morphological evolution. This review outlines the fundamental principles of the genetic code and genomic regulation interplay during development. Recent work that deepened our comprehension of cis-regulatory element origin, divergence and function is presented here to illustrate the state-of-the-art research that uncovered the principles of morphological novelty. This article is part of the theme issue 'The mammalian skull: development, structure and function'.

Disruption of T-box transcription factor eomesa results in abnormal development of median fins in Oujiang color common carp Cyprinus carpio

Median fins are thought to be ancestors of paired fins which in turn give rise to limbs in tetrapods. However, the developmental mechanisms of median fins remain largely unknown. Nonsense mutation of the T-box transcription factor eomesa in zebrafish results in a phenotype without dorsal fin. Compared to zebrafish, the common carp undergo an additional round of whole genome duplication, acquiring an extra copy of protein-coding genes. To verify the function of eomesa genes in common carp, we established a biallelic gene editing technology in this tetraploidy fish through simultaneous disruption of two homologous genes, eomesa1 and eomesa2. We targeted four sites located upstream or within the sequences encoding the T-box domain. Sanger sequencing data indicated the average knockout efficiency was around 40% at T1-T3 sites and 10% at T4 site in embryos at 24 hours post fertilization. The individual editing efficiency was high to about 80% at T1-T3 sites and low to 13.3% at T4 site in larvae at 7 days post fertilization. Among 145 mosaic F0 examined at four months old, three individuals (Mutant 1-3) showed varying degrees of maldevelopment in the dorsal fin and loss of anal fin. Genotyping showed the genomes of all three mutants were disrupted at T3 sites. The null mutation rates on the eomesa1 and eomesa2 loci were 0% and 60% in Mutant 1, 66.7% and 100% in Mutant 2, and 90% and 77.8% in Mutant 3, respectively. In conclusion, we demonstrated a role of eomesa in the formation and development of median fins in Oujiang color common carp and established an method that simultaneously disrupt two homologous genes with one gRNA, which would be useful in genome editing in other polyploidy fishes.

Recommendations for clinical interpretation of variants found in non-coding regions of the genome

BACKGROUND

The majority of clinical genetic testing focuses almost exclusively on regions of the genome that directly encode proteins. The important role of variants in non-coding regions in penetrant disease is, however, increasingly being demonstrated, and the use of whole genome sequencing in clinical diagnostic settings is rising across a large range of genetic disorders. Despite this, there is no existing guidance on how current guidelines designed primarily for variants in protein-coding regions should be adapted for variants identified in other genomic contexts.

METHODS

We convened a panel of nine clinical and research scientists with wide-ranging expertise in clinical variant interpretation, with specific experience in variants within non-coding regions. This panel discussed and refined an initial draft of the guidelines which were then extensively tested and reviewed by external groups.

RESULTS

We discuss considerations specifically for variants in non-coding regions of the genome. We outline how to define candidate regulatory elements, highlight examples of mechanisms through which non-coding region variants can lead to penetrant monogenic disease, and outline how existing guidelines can be adapted for the interpretation of these variants.

CONCLUSIONS

These recommendations aim to increase the number and range of non-coding region variants that can be clinically interpreted, which, together with a compatible phenotype, can lead to new diagnoses and catalyse the discovery of novel disease mechanisms.

Enhancer-gene specificity in development and disease

Enhancers control the establishment of spatiotemporal gene expression patterns throughout development. Over the past decade, the development of new technologies has improved our capacity to link enhancers with their target genes based on their colocalization within the same topological domains. However, the mechanisms that regulate how enhancers specifically activate some genes but not others within a given domain remain unclear. In this Review, we discuss recent insights into the factors controlling enhancer specificity, including the genetic composition of enhancers and promoters, the linear and 3D distance between enhancers and their target genes, and cell-type specific chromatin landscapes. We also discuss how elucidating the molecular principles of enhancer specificity might help us to better understand and predict the pathological consequences of human genetic, epigenetic and structural variants.

CNPY4 inhibits the Hedgehog pathway by modulating membrane sterol lipids

The Hedgehog (HH) pathway is critical for development and adult tissue homeostasis. Aberrant HH signaling can lead to congenital malformations and diseases including cancer. Although cholesterol and several oxysterol lipids have been shown to play crucial roles in HH activation, the molecular mechanisms governing their regulation remain unresolved. Here, we identify Canopy4 (CNPY4), a Saposin-like protein, as a regulator of the HH pathway that modulates levels of membrane sterol lipids. Cnpy4-/- embryos exhibit multiple defects consistent with HH signaling perturbations, most notably changes in digit number. Knockdown of Cnpy4 hyperactivates the HH pathway in vitro and elevates membrane levels of accessible sterol lipids, such as cholesterol, an endogenous ligand involved in HH activation. Our data demonstrate that CNPY4 is a negative regulator that fine-tunes HH signal transduction, revealing a previously undescribed facet of HH pathway regulation that operates through control of membrane composition.

North Carolina Macular Dystrophy: Phenotypic Variability and Computational Analysis of Disease-Associated Noncoding Variants

Purpose

North Carolina macular dystrophy (NCMD) is an autosomal dominant, congenital disorder affecting the central retina. Here, we report clinical and genetic findings in three families segregating NCMD and use epigenomic datasets from human tissues to gain insights into the effect of NCMD-implicated variants.

Methods

Clinical assessment and genetic testing were performed. Publicly available transcriptomic and epigenomic datasets were analyzed and the activity-by-contact method for scoring enhancer elements and linking them to target genes was used.

Results

A previously described, heterozygous, noncoding variant upstream of the PRDM13 gene was detected in all six affected study participants (chr6:100,040,987G>C [GRCh37/hg19]). Interfamilial and intrafamilial variability were observed; the visual acuity ranged from 0.0 to 1.6 LogMAR and fundoscopic findings ranged from visually insignificant, confluent, drusen-like macular deposits to coloboma-like macular lesions. Variable degrees of peripheral retinal spots (which were easily detected on widefield retinal imaging) were observed in all study subjects. Notably, a 6-year-old patient developed choroidal neovascularization and required treatment with intravitreal bevacizumab injections. Computational analysis of the five single nucleotide variants that have been implicated in NCMD revealed that these noncoding changes lie within two putative enhancer elements; these elements are predicted to interact with PRDM13 in the developing human retina. PRDM13 was found to be expressed in the fetal retina, with greatest expression in the amacrine precursor cell population.

Conclusions

We provide further evidence supporting the role of PRDM13 dysregulation in the pathogenesis of NCMD and highlight the usefulness of widefield retinal imaging in individuals suspected to have this condition.

Deletion of CTCF sites in the SHH locus alters enhancer-promoter interactions and leads to acheiropodia

Acheiropodia, congenital limb truncation, is associated with homozygous deletions in the LMBR1 gene around ZRS, an enhancer regulating SHH during limb development. How these deletions lead to this phenotype is unknown. Using whole-genome sequencing, we fine-mapped the acheiropodia-associated region to 12 kb and show that it does not function as an enhancer. CTCF and RAD21 ChIP-seq together with 4C-seq and DNA FISH identify three CTCF sites within the acheiropodia-deleted region that mediate the interaction between the ZRS and the SHH promoter. This interaction is substituted with other CTCF sites centromeric to the ZRS in the disease state. Mouse knockouts of the orthologous 12 kb sequence have no apparent abnormalities, showcasing the challenges in modelling CTCF alterations in animal models due to inherent motif differences between species. Our results show that alterations in CTCF motifs can lead to a Mendelian condition due to altered enhancer-promoter interactions.

Atlas of Musculoskeletal Stem Cells with the Soft and Hard Tissue Differentiation Architecture

Although being of utmost importance for human health and mobility, stem cell identity and hierarchical organization of musculoskeletal progenitors remain largely unexplored. Here, cells from E10.5, E12.5, and E15.5 murine limbs are analyzed by high throughput single-cell RNA sequencing to illustrate the cellular architecture during limb development. Single-cell transcriptional profiling demonstrates the identity and differentiation architecture of musculoskeletal stem cells (MSSC), soft and hard tissue progenitors through expression pattern of musculoskeletal markers (scleraxis [Scx], Hoxd13, Sox9, and Col1a1). This is confirmed by genetic in vivo lineage tracing. Moreover, single-cell analyses of Scx knockout mice tissues illustrates that Scx regulates MSSC self-renewal and proliferation potential. A high-throughput and low-cost multi-tissues RNA sequencing strategy further provides evidence that musculoskeletal system tissues, including muscle, bone, meniscus, and cartilage, are all abnormally developed in Scx knockout mice. These results establish the presence of an indispensable limb Scx+Hoxd13+ MSSC population and their differentiation into soft tissue progenitors (Scx+Col1a1+) and hard tissue progenitors (Scx+Sox9+). Collectively, this study paves the way for systematically decoding the complex molecular mechanisms and cellular programs of musculoskeletal tissues morphogenesis in limb development and regeneration.

Modulating gene regulation to treat genetic disorders

Over a thousand diseases are caused by mutations that alter gene expression levels. The potential of nuclease-deficient zinc fingers, TALEs or CRISPR fusion systems to treat these diseases by modulating gene expression has recently emerged. These systems can be applied to modify the activity of gene-regulatory elements - promoters, enhancers, silencers and insulators, subsequently changing their target gene expression levels to achieve therapeutic benefits - an approach termed cis-regulation therapy (CRT). Here, we review emerging CRT technologies and assess their therapeutic potential for treating a wide range of diseases caused by abnormal gene dosage. The challenges facing the translation of CRT into the clinic are discussed.

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.

Comprehensive In Vivo Interrogation Reveals Phenotypic Impact of Human Enhancer Variants

Establishing causal links between non-coding variants and human phenotypes is an increasing challenge. Here, we introduce a high-throughput mouse reporter assay for assessing the pathogenic potential of human enhancer variants in vivo and examine nearly a thousand variants in an enhancer repeatedly linked to polydactyly. We show that 71% of all rare non-coding variants previously proposed as causal lead to reporter gene expression in a pattern consistent with their pathogenic role. Variants observed to alter enhancer activity were further confirmed to cause polydactyly in knockin mice. We also used combinatorial and single-nucleotide mutagenesis to evaluate the in vivo impact of mutations affecting all positions of the enhancer and identified additional functional substitutions, including potentially pathogenic variants hitherto not observed in humans. Our results uncover the functional consequences of hundreds of mutations in a phenotype-associated enhancer and establish a widely applicable strategy for systematic in vivo evaluation of human enhancer variants.

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.

Predicting the effects of SNPs on transcription factor binding affinity

MOTIVATION

Genome-wide association studies have revealed that 88% of disease-associated single-nucleotide polymorphisms (SNPs) reside in noncoding regions. However, noncoding SNPs remain understudied, partly because they are challenging to prioritize for experimental validation. To address this deficiency, we developed the SNP effect matrix pipeline (SEMpl).

RESULTS

SEMpl estimates transcription factor-binding affinity by observing differences in chromatin immunoprecipitation followed by deep sequencing signal intensity for SNPs within functional transcription factor-binding sites (TFBSs) genome-wide. By cataloging the effects of every possible mutation within the TFBS motif, SEMpl can predict the consequences of SNPs to transcription factor binding. This knowledge can be used to identify potential disease-causing regulatory loci.

AVAILABILITY AND IMPLEMENTATION

SEMpl is available from https://github.com/Boyle-Lab/SEM_CPP.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Holt-Oram syndrome: clinical and molecular description of 78 patients with TBX5 variants

Holt-Oram syndrome (HOS) is an autosomal dominant condition characterised by the association of congenital heart defect (CHD), with or without rhythm disturbances and radial defects, due to TBX5 variants. The diagnosis is challenged by the variability of expression and the large phenotypic overlap with other conditions, like Okihiro syndrome, TAR syndrome or Fanconi disease. We retrospectively reviewed 212 patients referred for suspicion of HOS between 2002 and 2014, who underwent TBX5 screening. A TBX5 variant has been identified in 78 patients, representing the largest molecular series ever described. In the cohort, 61 met the previously described diagnostic criteria and 17 have been considered with an uncertain HOS diagnosis. A CHD was present in 91% of the patients with a TBX5 variant, atrial septal defects being the most common (61.5%). The genotype-phenotype study highlights the importance of some critical features in HOS: the septal characteristic of the CHD, the bilateral and asymmetric characteristics of the radial defect and the presence of shoulder or elbow mobility defect. Besides, 21 patients presented with an overlapping condition. Among them, 13 had a typical HOS presentation. We discuss the strategies that could be adopted to improve the molecular delineation of the remaining typical patients.

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.

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.

The developmental and genetic basis of 'clubfoot' in the peroneal muscular atrophy mutant mouse

Genetic factors underlying the human limb abnormality congenital talipes equinovarus ('clubfoot') remain incompletely understood. The spontaneous autosomal recessive mouse 'peroneal muscular atrophy' mutant (PMA) is a faithful morphological model of human clubfoot. In PMA mice, the dorsal (peroneal) branches of the sciatic nerves are absent. In this study, the primary developmental defect was identified as a reduced growth of sciatic nerve lateral motor column (LMC) neurons leading to failure to project to dorsal (peroneal) lower limb muscle blocks. The pma mutation was mapped and a candidate gene encoding LIM-domain kinase 1 (Limk1) identified, which is upregulated in mutant lateral LMC motor neurons. Genetic and molecular analyses showed that the mutation acts in the EphA4-Limk1-Cfl1/cofilin-actin pathway to modulate growth cone extension/collapse. In the chicken, both experimental upregulation of Limk1 by electroporation and pharmacological inhibition of actin turnover led to defects in hindlimb spinal motor neuron growth and pathfinding, and mimicked the clubfoot phenotype. The data support a neuromuscular aetiology for clubfoot and provide a mechanistic framework to understand clubfoot in humans.

Two Types of Etiological Mutation in the Limb-Specific Enhancer of Shh

An enhancer named MFCS1 regulates Sonic hedgehog (Shh) expression in the posterior mesenchyme of limb buds. Several mutations in MFCS1 induce ectopic Shh expression in the anterior limb bud, and these result in preaxial polydactyly (PPD). However, the molecular basis of ectopic Shh expression remains elusive, although some mutations are known to disrupt the negative regulation of Shh expression in the anterior limb bud. Here, we analyzed the molecular mechanism of ectopic Shh expression in PPD including in a mouse mutation-hemimelic extra toes (Hx)-and in other MFCS1 mutations in different species. First, we generated transgenic mouse lines with a LacZ reporter cassette flanked with tandem repeats of 40 bp MFCS1 fragments harboring a mutation. The transgenic mouse line with the Hx-type fragment showed reporter expression exclusively in the anterior, but not in the posterior margins of limb buds. In contrast, no specific LacZ expression was observed in lines carrying the MFCS1 fragment with other mutations. Yeast one-hybrid assays revealed that the msh-like homeodomain protein, MSX1, bound specifically to the Hx sequence of MFCS1. Thus, PPD caused by mutations in MFCS1 has two major types of molecular etiology: loss of a cis-motif for negative regulation of Shh, and acquisition of a new cis-motif binding to a preexisting transcription factor, as represented by the Hx mutation.

Variant Interpretation: Functional Assays to the Rescue

Classical genetic approaches for interpreting variants, such as case-control or co-segregation studies, require finding many individuals with each variant. Because the overwhelming majority of variants are present in only a few living humans, this strategy has clear limits. Fully realizing the clinical potential of genetics requires that we accurately infer pathogenicity even for rare or private variation. Many computational approaches to predicting variant effects have been developed, but they can identify only a small fraction of pathogenic variants with the high confidence that is required in the clinic. Experimentally measuring a variant's functional consequences can provide clearer guidance, but individual assays performed only after the discovery of the variant are both time and resource intensive. Here, we discuss how multiplex assays of variant effect (MAVEs) can be used to measure the functional consequences of all possible variants in disease-relevant loci for a variety of molecular and cellular phenotypes. The resulting large-scale functional data can be combined with machine learning and clinical knowledge for the development of "lookup tables" of accurate pathogenicity predictions. A coordinated effort to produce, analyze, and disseminate large-scale functional data generated by multiplex assays could be essential to addressing the variant-interpretation crisis.

Whole-genome sequencing approaches for conservation biology: Advantages, limitations and practical recommendations

Whole-genome resequencing (WGR) is a powerful method for addressing fundamental evolutionary biology questions that have not been fully resolved using traditional methods. WGR includes four approaches: the sequencing of individuals to a high depth of coverage with either unresolved or resolved haplotypes, the sequencing of population genomes to a high depth by mixing equimolar amounts of unlabelled-individual DNA (Pool-seq) and the sequencing of multiple individuals from a population to a low depth (lcWGR). These techniques require the availability of a reference genome. This, along with the still high cost of shotgun sequencing and the large demand for computing resources and storage, has limited their implementation in nonmodel species with scarce genomic resources and in fields such as conservation biology. Our goal here is to describe the various WGR methods, their pros and cons and potential applications in conservation biology. WGR offers an unprecedented marker density and surveys a wide diversity of genetic variations not limited to single nucleotide polymorphisms (e.g., structural variants and mutations in regulatory elements), increasing their power for the detection of signatures of selection and local adaptation as well as for the identification of the genetic basis of phenotypic traits and diseases. Currently, though, no single WGR approach fulfils all requirements of conservation genetics, and each method has its own limitations and sources of potential bias. We discuss proposed ways to minimize such biases. We envision a not distant future where the analysis of whole genomes becomes a routine task in many nonmodel species and fields including conservation biology.

Gene Regulatory Elements, Major Drivers of Human Disease

Gene expression changes, the driving forces for cellular diversity in multicellular organisms, are regulated by a diverse set of gene regulatory elements that direct transcription in specific cells. Mutations in these elements, ranging from chromosomal aberrations to single-nucleotide polymorphisms, are a major cause of human disease. However, we currently have a very limited understanding of how regulatory element genotypes lead to specific phenotypes. In this review, we discuss the various methods of regulatory element identification, the different types of mutations they harbor, and their impact on human disease. We highlight how these variations can affect transcription of multiple genes in gene regulatory networks. In addition, we describe how novel technologies, such as massively parallel reporter assays and CRISPR/Cas9 genome editing, are beginning to provide a better understanding of the functional roles that these elements have and how their alteration can lead to specific phenotypes.

Limb development: a paradigm of gene regulation

The limb is a commonly used model system for developmental biology. Given the need for precise control of complex signalling pathways to achieve proper patterning, the limb is also becoming a model system for gene regulation studies. Recent developments in genomic technologies have enabled the genome-wide identification of regulatory elements that control limb development, yielding insights into the determination of limb morphology and forelimb versus hindlimb identity. The modulation of regulatory interactions - for example, through the modification of regulatory sequences or chromatin architecture - can lead to morphological evolution, acquired regeneration capacity or limb malformations in diverse species, including humans.

Progressive Loss of Function in a Limb Enhancer during Snake Evolution

The evolution of body shape is thought to be tightly coupled to changes in regulatory sequences, but specific molecular events associated with major morphological transitions in vertebrates have remained elusive. We identified snake-specific sequence changes within an otherwise highly conserved long-range limb enhancer of Sonic hedgehog (Shh). Transgenic mouse reporter assays revealed that the in vivo activity pattern of the enhancer is conserved across a wide range of vertebrates, including fish, but not in snakes. Genomic substitution of the mouse enhancer with its human or fish ortholog results in normal limb development. In contrast, replacement with snake orthologs caused severe limb reduction. Synthetic restoration of a single transcription factor binding site lost in the snake lineage reinstated full in vivo function to the snake enhancer. Our results demonstrate changes in a regulatory sequence associated with a major body plan transition and highlight the role of enhancers in morphological evolution. PAPERCLIP.

An increased duplication of ZRS region that caused more than one supernumerary digits preaxial polydactyly in a large Chinese family

Preaxial polydactyly (PPD) is inherited in an autosomal dominant fashion and characterized by the presence of one or more supernumerary digits on the thumb side. It had been identified that point mutation or genomic duplications of the long-range limb-specific cis-regulator - zone of polarizing activity regulatory sequence (ZRS) cause PPD or other limb deformities such as syndactyly type IV (SD4) and Triphalangeal thumb-polysyndactyly syndrome (TPTPS). Most previously reported cases involved with no more than one extra finger; however, the role of the point mutation or genomic duplications of ZRS in the case of more than one redundant finger polydactyly remains unclear. In this article, we reported a family case of more than one redundant finger polydactyly on the thumb side for bilateral hands with a pedigree chart of the family. Results of quantitative PCR (qPCR) and sequence analysis suggested that the relative copy number (RCN) of ZRS but not point mutation (including insertion and deletion) was involved in all affected individuals.

A Novel ZRS Mutation in a Chinese Patient with Preaxial Polydactyly and Triphalangeal Thumb

Preaxial polydactyly (PPD; OMIM 603596), which is characterized as having supernumerary fingers, is an unusual congenital hand abnormality. Triphalangeal thumb (TPT; OMIM 190600) is identified by an extra phalangeal bone and is often found in association with PPD. When in combination, the disease is referred to as PPD type II (PPD II; OMIM 174500). Previous studies have demonstrated that variations in the zone of polarizing activity regulatory sequence (ZRS; chr7:156,583,796-156,584,569; hg19) region are associated with PPD II. In this study, our patient was diagnosed with PPD II, having bilateral thumb duplication and unilateral TPT (on the right hand). Further investigation of possible causative genes identified a de novo heterozygous ZRS mutation (ZRS 428T>A). This novel mutation was neither found in 200 normal controls nor reported in online databases. Moreover, the bioinformatics program Genomic Evolutionary Rate Profiling (GERP) revealed this site (ZRS428) to be evolutionarily highly conserved, and the 428T>A point mutation was predicted to be deleterious by MutationTaster. In conclusion, the affected individual shows bilateral thumb duplication, but unilateral TPT making this case special. Thus, our findings not only further support the important role of ZRS in limb morphogenesis and expand the spectrum of ZRS mutations, but also emphasize the significance of genetic diagnosis and counseling of families with digit number and identity alterations as well.

Genomic features of human limb specific enhancers

To elucidate important cellular and molecular interactions that regulate patterning and skeletal development, vertebrate limbs served as a model organ. A growing body of evidence from detailed studies on a subset of limb regulators like the HOXD cluster or SHH, reveals the importance of enhancers in limb related developmental and disease processes. Exploiting the recent genome-wide availability of functionally confirmed enhancer dataset, this study establishes regulatory interactions for dozens of human limb developmental genes. From these data, it appears that the long-range regulatory interactions are fairly common during limb development. This observation highlights the significance of chromosomal breaks/translocations in human limb deformities. Transcriptional factor (TF) analysis predicts that the differentiation of early nascent limb-bud into future territories entail distinct TF interaction networks. Conclusively, an important motivation for annotating the human limb specific regulatory networks is to pave way for the systematic exploration of their role in disease and evolution.

Integrative Tissue-Specific Functional Annotations in the Human Genome Provide Novel Insights on Many Complex Traits and Improve Signal Prioritization in Genome Wide Association Studies

Extensive efforts have been made to understand genomic function through both experimental and computational approaches, yet proper annotation still remains challenging, especially in non-coding regions. In this manuscript, we introduce GenoSkyline, an unsupervised learning framework to predict tissue-specific functional regions through integrating high-throughput epigenetic annotations. GenoSkyline successfully identified a variety of non-coding regulatory machinery including enhancers, regulatory miRNA, and hypomethylated transposable elements in extensive case studies. Integrative analysis of GenoSkyline annotations and results from genome-wide association studies (GWAS) led to novel biological insights on the etiologies of a number of human complex traits. We also explored using tissue-specific functional annotations to prioritize GWAS signals and predict relevant tissue types for each risk locus. Brain and blood-specific annotations led to better prioritization performance for schizophrenia than standard GWAS p-values and non-tissue-specific annotations. As for coronary artery disease, heart-specific functional regions was highly enriched of GWAS signals, but previously identified risk loci were found to be most functional in other tissues, suggesting a substantial proportion of still undetected heart-related loci. In summary, GenoSkyline annotations can guide genetic studies at multiple resolutions and provide valuable insights in understanding complex diseases. GenoSkyline is available at http://genocanyon.med.yale.edu/GenoSkyline.

Bat Accelerated Regions Identify a Bat Forelimb Specific Enhancer in the HoxD Locus

The molecular events leading to the development of the bat wing remain largely unknown, and are thought to be caused, in part, by changes in gene expression during limb development. These expression changes could be instigated by variations in gene regulatory enhancers. Here, we used a comparative genomics approach to identify regions that evolved rapidly in the bat ancestor, but are highly conserved in other vertebrates. We discovered 166 bat accelerated regions (BARs) that overlap H3K27ac and p300 ChIP-seq peaks in developing mouse limbs. Using a mouse enhancer assay, we show that five Myotis lucifugus BARs drive gene expression in the developing mouse limb, with the majority showing differential enhancer activity compared to the mouse orthologous BAR sequences. These include BAR116, which is located telomeric to the HoxD cluster and had robust forelimb expression for the M. lucifugus sequence and no activity for the mouse sequence at embryonic day 12.5. Developing limb expression analysis of Hoxd10-Hoxd13 in Miniopterus natalensis bats showed a high-forelimb weak-hindlimb expression for Hoxd10-Hoxd11, similar to the expression trend observed for M. lucifugus BAR116 in mice, suggesting that it could be involved in the regulation of the bat HoxD complex. Combined, our results highlight novel regulatory regions that could be instrumental for the morphological differences leading to the development of the bat wing.

Congenital Anomalies of the Hand--Principles of Management

Physicians who specialize in pediatric orthopedics and hand surgery frequently encounter congenital hand abnormalities, despite their relative rarity. The treating physician should be aware of the associated syndromes and malformations that may, in some cases, be fatal if not recognized and treated appropriately. Although these congenital disorders have a wide variability, their treatment principles are similar in that the physician should promote functional use and cosmesis for the hand. This article discusses syndactyly, preaxial polydactyly and post-axial polydactyly, and the hypoplastic thumb.

Minor Loops in Major Folds: Enhancer-Promoter Looping, Chromatin Restructuring, and Their Association with Transcriptional Regulation and Disease

The organization and folding of chromatin within the nucleus can determine the outcome of gene expression. Recent technological advancements have enabled us to study chromatin interactions in a genome-wide manner at high resolution. These studies have increased our understanding of the hierarchy and dynamics of chromatin domains that facilitate cognate enhancer–promoter looping, defining the transcriptional program of different cell types. In this review, we focus on vertebrate chromatin long-range interactions as they relate to transcriptional regulation. In addition, we describe how the alteration of boundaries that mark discrete regions in the genome with high interaction frequencies within them, called topological associated domains (TADs), could lead to various phenotypes, including human diseases, which we term as “TADopathies.”

The Evolution of Bony Vertebrate Enhancers at Odds with Their Coding Sequence Landscape

Enhancers lie at the heart of transcriptional and developmental gene regulation. Therefore, changes in enhancer sequences usually disrupt the target gene expression and result in disease phenotypes. Despite the well-established role of enhancers in development and disease, evolutionary sequence studies are lacking. The current study attempts to unravel the puzzle of bony vertebrates’ conserved noncoding elements (CNE) enhancer evolution. Bayesian phylogenetics of enhancer sequences spotlights promising interordinal relationships among placental mammals, proposing a closer relationship between humans and laurasiatherians while placing rodents at the basal position. Clock-based estimates of enhancer evolution provided a dynamic picture of interspecific rate changes across the bony vertebrate lineage. Moreover, coelacanth in the study augmented our appreciation of the vertebrate cis-regulatory evolution during water–land transition. Intriguingly, we observed a pronounced upsurge in enhancer evolution in land-dwelling vertebrates. These novel findings triggered us to further investigate the evolutionary trend of coding as well as CNE nonenhancer repertoires, to highlight the relative evolutionary dynamics of diverse genomic landscapes. Surprisingly, the evolutionary rates of enhancer sequences were clearly at odds with those of the coding and the CNE nonenhancer sequences during vertebrate adaptation to land, with land vertebrates exhibiting significantly reduced rates of coding sequence evolution in comparison to their fast evolving regulatory landscape. The observed variation in tetrapod cis-regulatory elements caused the fine-tuning of associated gene regulatory networks. Therefore, the increased evolutionary rate of tetrapods’ enhancer sequences might be responsible for the variation in developmental regulatory circuits during the process of vertebrate adaptation to land.

A statistical framework to predict functional non-coding regions in the human genome through integrated analysis of annotation data

Identifying functional regions in the human genome is a major goal in human genetics. Great efforts have been made to functionally annotate the human genome either through computational predictions, such as genomic conservation, or high-throughput experiments, such as the ENCODE project. These efforts have resulted in a rich collection of functional annotation data of diverse types that need to be jointly analyzed for integrated interpretation and annotation. Here we present GenoCanyon, a whole-genome annotation method that performs unsupervised statistical learning using 22 computational and experimental annotations thereby inferring the functional potential of each position in the human genome. With GenoCanyon, we are able to predict many of the known functional regions. The ability of predicting functional regions as well as its generalizable statistical framework makes GenoCanyon a unique and powerful tool for whole-genome annotation. The GenoCanyon web server is available at http://genocanyon.med.yale.edu

Cis-regulatory somatic mutations and gene-expression alteration in B-cell lymphomas

BACKGROUND

With the rapid increase of whole-genome sequencing of human cancers, an important opportunity to analyze and characterize somatic mutations lying within cis-regulatory regions has emerged. A focus on protein-coding regions to identify nonsense or missense mutations disruptive to protein structure and/or function has led to important insights; however, the impact on gene expression of mutations lying within cis-regulatory regions remains under-explored. We analyzed somatic mutations from 84 matched tumor-normal whole genomes from B-cell lymphomas with accompanying gene expression measurements to elucidate the extent to which these cancers are disrupted by cis-regulatory mutations.

RESULTS

We characterize mutations overlapping a high quality set of well-annotated transcription factor binding sites (TFBSs), covering a similar portion of the genome as protein-coding exons. Our results indicate that cis-regulatory mutations overlapping predicted TFBSs are enriched in promoter regions of genes involved in apoptosis or growth/proliferation. By integrating gene expression data with mutation data, our computational approach culminates with identification of cis-regulatory mutations most likely to participate in dysregulation of the gene expression program. The impact can be measured along with protein-coding mutations to highlight key mutations disrupting gene expression and pathways in cancer.

CONCLUSIONS

Our study yields specific genes with disrupted expression triggered by genomic mutations in either the coding or the regulatory space. It implies that mutated regulatory components of the genome contribute substantially to cancer pathways. Our analyses demonstrate that identifying genomically altered cis-regulatory elements coupled with analysis of gene expression data will augment biological interpretation of mutational landscapes of cancers.

The disruption of a novel limb cis-regulatory element of SHH is associated with autosomal dominant preaxial polydactyly-hypertrichosis

The expression gradient of the morphogen Sonic Hedgehog (SHH) is crucial in establishing the number and the identity of the digits during anteroposterior patterning of the limb. Its anterior ectopic expression is responsible for preaxial polydactyly (PPD). Most of these malformations are due to the gain-of-function of the Zone of Polarizing Activity Regulatory Sequence, the only limb-specific enhancer of SHH known to date. We report a family affected with a novel condition associating PPD and hypertrichosis of the upper back, following an autosomal dominant mode of inheritance. This phenotype is consistent with deregulation of SHH expression during limb and follicle development. In affected members, we identified a 2 kb deletion located ~240 kb upstream from the SHH promoter. The deleted sequence is capable of repressing the transcriptional activity of the SHH promoter in vitro, consistent with a silencer activity. We hypothesize that the deletion of this silencer could be responsible for SHH deregulation during development, leading to a PPD-hypertrichosis phenotype.

Genetic manipulation of a high-affinity PHR1 target cis-element to improve phosphorous uptake in Oryza sativa L

Phosphorus (P) is an essential macronutrient for crop development and production. Phosphate starvation response 1 (PHR1) acts as the central regulator for Pi-signaling and Pi-homeostasis in plants by binding to the cis-element PHR1 binding sequence (P1BS; GNATATNC). However, how phosphate starvation-induced gene expression is regulated remains obscure. In this work, we investigated the DNA binding affinity of the PHR1 ortholog OsPHR2 to its downstream target genes in Oryza sativa (rice). We confirmed that a combination of P1BS and P1BS-like motifs are essential for stable binding by OsPHR2. Furthermore, we report that variations in P1BS motif bases affected the binding affinity of OsPHR2 and that the highest affinity motif was GaATATtC (designated the A-T-type P1BS). We also found that a combination of two A-T-type P1BS elements in tandem, namely HA-P1BS, was very efficient for binding of OsPHR2. Using the cis-regulator HA-P1BS, we modified the promoters of Transporter Traffic Facilitator 1 (PHF1), a key factor controlling endoplasmic reticulum-exit of phosphate transporters to the plasma membrane, for efficient uptake of phosphorous in an energetically neutral way. Transgenic plants with the modified promoters showed significantly enhanced tolerance to low phosphate stress in both solution and soil conditions, which provides a new strategy for crop improvement to enhance tolerance of nutrient deficiency.

Intrafamilial variability of ZRS-associated syndrome: characterization of a mosaic ZRS mutation by pyrosequencing

During limb development, the spatio-temporal expression of sonic hedgehog (SHH) is driven by the Zone of polarizing activity Regulatory Sequence (ZRS), located 1 megabase upstream from SHH. Gain-of-function mutations of this enhancer, which cause ectopic expression of SHH, are known to be responsible for congenital limb malformations with variable expressivity, ranging from preaxial polydactyly or triphalangeal thumbs to polysyndactyly, which may also be associated with mesomelic deficiency. In this report, we describe a patient affected with mirror-image polydactyly of the four extremities and bilateral tibial deficiency. The proband's father had isolated preaxial polydactyly type II (PPD2). Using Sanger sequencing, a ZRS point mutation (NC_000007.14, g.156584153A>G, UCSC, Build hg.19) was only identified in the patient. However, pyrosequencing analysis enabled the detection of a 10% somatic mosaic in the blood and saliva from the father. To our knowledge, this is the first description of a ZRS mosaic mutation. This report highlights the complexity of genotype-phenotype correlation in ZRS-associated syndromes and the importance of detecting somatic mosaicism for accurate genetic counselling.

Genome-wide identification of signaling center enhancers in the developing limb

The limb is widely used as a model developmental system and changes to gene expression patterns in its signaling centers, notably the zone of polarizing activity (ZPA) and the apical ectodermal ridge (AER), are known to cause limb malformations and evolutionary differences in limb morphology. Although several genes that define these limb signaling centers have been described, the identification of regulatory elements that are active within these centers has been limited. By dissecting mouse E11.5 limbs that fluorescently mark the ZPA or AER, followed by fluorescence-activated cell sorting and low-cell H3K27ac ChIP-seq, we identified thousands of specific signaling-center enhancers. Our ChIP-seq datasets show strong correlation with ZPA- and AER-expressed genes, previously characterized functional ZPA and AER enhancers and enrichment for relevant biological terms related to limb development and malformation for the neighboring genes. Using transgenic assays, we show that several of these sequences function as ZPA and AER enhancers. Our results identify novel ZPA and AER enhancers that could be important regulators of genes involved in the establishment of these specialized regions and the patterning of tetrapod limbs.

A novel ZRS mutation leads to preaxial polydactyly type 2 in a heterozygous form and Werner mesomelic syndrome in a homozygous form

Point mutations in the zone of polarizing activity regulatory sequence (ZRS) are known to cause human limb malformations. Although most mutations cause preaxial polydactyly (PPD), triphalangeal thumb (TPT) or both, a mutation in position 404 of the ZRS causes more severe Werner mesomelic syndrome (WMS) for which malformations include the distal arm or leg bones in addition to the hands and/or feet. Of more than 15 reported families with ZRS mutations, only one homozygous individual has been reported, with no change in phenotype compared with heterozygotes. Here, we describe a novel point mutation in the ZRS, 402C>T (AC007097.4:g.105548C>T), that is transmitted through two Mexican families with one homozygous individual. The homozygous phenotype for this mutation, WMS, is more severe than the numerous heterozygous individuals genotyped from both families who have TPT and PPD. A mouse transgenic enhancer assay shows that this mutation causes an expansion of the enhancer's expression domain in the developing mouse limb, confirming its pathogenicity. Combined, our results identify a novel ZRS mutation in the Mexican population, 402C>T, and suggest that a dosage effect exists for this ZRS mutation.

Epidemiology of congenital upper limb anomalies in Stockholm, Sweden, 1997 to 2007: application of the Oberg, Manske, and Tonkin classification

PURPOSE

To investigate the epidemiology of congenital upper limb anomalies (CULA) based on the newly proposed Oberg, Manske, and Tonkin (OMT) classification, to compare this classification with the International Federation of Societies for Surgery of the Hand (IFSSH) classification, and to provide incidence rates of the different CULA.

METHODS

In this study, the same 562 individuals with a CULA who were analyzed in a previous epidemiologic study based on the IFSSH classification were reclassified according to the OMT classification. All children identified with CULA and born in Stockholm County between January 1, 1997 and December 31, 2007 were included in the study. During the period there were 261,914 live births in Stockholm County, and the population of Stockholm County was 1,949,516 inhabitants at the end of the period. From medical records and available radiographs, all cases were analyzed regarding type of CULA, sex, affected side, associated nonhand anomalies, and occurrence among relatives. Individuals with right and left side anomalies belonging to different OMT subgroups were counted as 2 anomalies; thus, the material consisted of 577 CULA in 562 children.

RESULTS

It was possible to organize all CULA into the OMT classification. The largest main category was malformations (429 cases), followed by deformations (124 cases), dysplasias (10 cases), and syndromes (14 cases). We present the relation between the IFSSH and OMT classifications, elucidate difficulties within the OMT classification, and propose additions to the classification.

CONCLUSIONS

This study confirms that the OMT classification is useful and accurate, but also points out difficulties. With further refinements, we regard the OMT classification as a needed and appropriate replacement for the IFSSH classification.

TYPE OF STUDY/LEVEL OF EVIDENCE

Diagnostic III.

Confirmation of genetic homogeneity of syndactyly type IV and triphalangeal thumb-polysyndactyly syndrome in a Chinese family and review of the literature

Syndactyly type IV (SD4) is inherited in an autosomal dominant fashion and characterized by complete cutaneous syndactyly of all fingers accompanied with polydactyly. Triphalangeal thumb-polysyndactyly syndrome (TPTPS) consists of a triphalangeal thumb, polydactyly, and syndactyly and is transmitted in an autosomal dominant manner with variable expression. Genomic duplications of the long-range limb-specific cis-regulator (ZRS) cause SD4 and TPTPS. Here, we report two individuals from a Chinese family with syndactyly. One individual had overlapping clinical symptoms of TPTPS and SD4, while the other had a typical SD4 with postaxial polydactyly of the toe. Results of quantitative PCR suggested that the duplication of ZRS involved all affected individuals, and array comparative genomic hybridization detected its size as 115.3 kb.

CONCLUSION

This work confirms the genetic homogeneity of SD4 and TPTPS. Our result expands the spectrum of ZRS duplications. TPTPS and SD4 should be considered as a continuum of phenotypes.

Currarino Syndrome and HPE Microform Associated with a 2.7-Mb Deletion in 7q36.3 Excluding SHH Gene

Holoprosencephaly (HPE) is the most common forebrain defect in humans. It results from incomplete midline cleavage of the prosencephalon and can be caused by environmental and genetic factors. HPE is usually described as a continuum of brain malformations from the most severe alobar HPE to the middle interhemispheric fusion variant or syntelencephaly. A microform of HPE is limited to craniofacial features such as congenital nasal pyriform aperture stenosis and single central maxillary incisor, without brain malformation. Among the heterogeneous causes of HPE, point mutations and deletions in the SHH gene at 7q36 have been identified as well as extremely rare chromosomal rearrangements in the long-range enhancers of this gene. Here, we report a boy with an HPE microform associated with a Currarino syndrome. Array CGH detected a de novo 2.7-Mb deletion in the 7q36.3 region including the MNX1 gene, usually responsible for the Currarino triad but excluding SHH, which is just outside the deletion. This new case provides further evidence of the importance of the SHH long-range enhancers in the HPE spectrum.

Identification of biologically relevant enhancers in human erythroid cells

Identification of cell type-specific enhancers is important for understanding the regulation of programs controlling cellular development and differentiation. Enhancers are typically marked by the co-transcriptional activator protein p300 or by groups of cell-expressed transcription factors. We hypothesized that a unique set of enhancers regulates gene expression in human erythroid cells, a highly specialized cell type evolved to provide adequate amounts of oxygen throughout the body. Using chromatin immunoprecipitation followed by massively parallel sequencing, genome-wide maps of candidate enhancers were constructed for p300 and four transcription factors, GATA1, NF-E2, KLF1, and SCL, using primary human erythroid cells. These data were combined with gene expression analyses, and candidate enhancers were identified. Consistent with their predicted function as candidate enhancers, there was statistically significant enrichment of p300 and combinations of co-localizing erythroid transcription factors within 1-50 kb of the transcriptional start site (TSS) of genes highly expressed in erythroid cells. Candidate enhancers were also enriched near genes with known erythroid cell function or phenotype. Candidate enhancers exhibited moderate conservation with mouse and minimal conservation with nonplacental vertebrates. Candidate enhancers were mapped to a set of erythroid-associated, biologically relevant, SNPs from the genome-wide association studies (GWAS) catalogue of NHGRI, National Institutes of Health. Fourteen candidate enhancers, representing 10 genetic loci, mapped to sites associated with biologically relevant erythroid traits. Fragments from these loci directed statistically significant expression in reporter gene assays. Identification of enhancers in human erythroid cells will allow a better understanding of erythroid cell development, differentiation, structure, and function and provide insights into inherited and acquired hematologic disease.

Functional characterization of tissue-specific enhancers in the DLX5/6 locus

Disruption of distaless homeobox 5 and 6 (Dlx5/6) in mice results in brain, craniofacial, genital, ear and limb defects. In humans, chromosomal aberrations in the DLX5/6 region, some of which do not encompass DLX5/6, are associated with split hand/foot malformation 1 (SHFM1) as well as intellectual disability, craniofacial anomalies and hearing loss, suggesting that the disruption of DLX5/6 regulatory elements could lead to these abnormalities. Here, we characterized enhancers in the DLX5/6 locus whose tissue-specific expression and genomic location along with previously characterized enhancers correlate with phenotypes observed in individuals with chromosomal abnormalities. By analyzing chromosomal aberrations at 7q21, we refined the minimal SHFM1 critical region and used comparative genomics to select 26 evolutionary conserved non-coding sequences in this critical region for zebrafish enhancer assays. Eight of these sequences were shown to function as brain, olfactory bulb, branchial arch, otic vesicle and fin enhancers, recapitulating dlx5a/6a expression. Using a mouse enhancer assay, several of these zebrafish enhancers showed comparable expression patterns in the branchial arch, otic vesicle, forebrain and/or limb at embryonic day 11.5. Examination of the coordinates of various chromosomal rearrangements in conjunction with the genomic location of these tissue-specific enhancers showed a correlation with the observed clinical abnormalities. Our findings suggest that chromosomal abnormalities that disrupt the function of these tissue-specific enhancers could be the cause of SHFM1 and its associated phenotypes. In addition, they highlight specific enhancers in which mutations could lead to non-syndromic hearing loss, craniofacial defects or limb malformations.

A novel mutation in the SHH long-range regulator (ZRS) is associated with preaxial polydactyly, triphalangeal thumb, and severe radial ray deficiency

Sonic Hedgehog (SHH) within the posteriorly located zone of polarizing activity is the main controller of the antero-posterior axis of limb development. The ZRS (zone of polarizing activity regulatory sequence) is a long-range limb-specific SHH enhancer. Several point mutations in the ZRS have been described in humans. These mutations cause enhanced SHH activity and ectopic anterior expression of SHH and a variable phenotype of preaxial polydactyly and triphalangeal thumb. Absent thumb or radius has not been reported with ZRS mutations. Here, we report on a family with a variable phenotype of preaxial polydactyly as well as absent thumb and radius, with kidney and cardiac defects. The family was screened for SALL1, SALL4, and TBX5 mutations, but all were normal. Finally, they were screened for ZRS mutations, which showed a novel point mutation within the ZRS, NG_009240.1: g.106954C>T (traditional nomenclature: ZRS619C>T) in the five affected members. This mutation was not previously reported in any public domain database, and was not found in our healthy and ethnically matched control individuals or unaffected family members. We hypothesize that interactions of SHH and SALL1 explain the overlapping features of the family described here and patients with Townes-Brocks syndrome.

Conserved cis-regulatory regions in a large genomic landscape control SHH and BMP-regulated Gremlin1 expression in mouse limb buds

Background

Mouse limb bud is a prime model to study the regulatory interactions that control vertebrate organogenesis. Major aspects of limb bud development are controlled by feedback loops that define a self-regulatory signalling system. The SHH/GREM1/AER-FGF feedback loop forms the core of this signalling system that operates between the posterior mesenchymal organiser and the ectodermal signalling centre. The BMP antagonist Gremlin1 (GREM1) is a critical node in this system, whose dynamic expression is controlled by BMP, SHH, and FGF signalling and key to normal progression of limb bud development. Previous analysis identified a distant cis-regulatory landscape within the neighbouring Formin1 (Fmn1) locus that is required for Grem1 expression, reminiscent of the genomic landscapes controlling HoxD and Shh expression in limb buds.

Results

Three highly conserved regions (HMCO1-3) were identified within the previously defined critical genomic region and tested for their ability to regulate Grem1 expression in mouse limb buds. Using a combination of BAC and conventional transgenic approaches, a 9 kb region located ~70 kb downstream of the Grem1 transcription unit was identified. This region, termed Grem1 Regulatory Sequence 1 (GRS1), is able to recapitulate major aspects of Grem1 expression, as it drives expression of a LacZ reporter into the posterior and, to a lesser extent, in the distal-anterior mesenchyme. Crossing the GRS1 transgene into embryos with alterations in the SHH and BMP pathways established that GRS1 depends on SHH and is modulated by BMP signalling, i.e. integrates inputs from these pathways. Chromatin immunoprecipitation revealed interaction of endogenous GLI3 proteins with the core cis-regulatory elements in the GRS1 region. As GLI3 is a mediator of SHH signal transduction, these results indicated that SHH directly controls Grem1 expression through the GRS1 region. Finally, all cis-regulatory regions within the Grem1 genomic landscape locate to the DNAse I hypersensitive sites identified in this genomic region by the ENCODE consortium.

Conclusions

This study establishes that distant cis-regulatory regions scattered through a larger genomic landscape control the highly dynamic expression of Grem1, which is key to normal progression of mouse limb bud development.

A novel ZRS mutation in a Balochi tribal family with triphalangeal thumb, pre-axial polydactyly, post-axial polydactyly, and syndactyly

Limb malformations are one of the most common types of human congenital malformations. Mutations in the ZRS enhancer of Sonic Hedgehog are thought to be responsible for pre-axial polydactyly in multiple independent families. Here, we describe a large Balochi tribal family from Southern Punjab, Pakistan, with a variable set of limb malformations and a novel ZRS mutation. The family has a limb phenotype characterized by triphalangeal thumb, pre-axial polydactyly, and post-axial polydactyly. There is also a high degree of phenotypic heterogeneity with less common clinical findings in the affected family members that include osseous syndactyly of forth-fifth fingers, clinodactyly, hypoplasia of mesoaxial fingers, and bifid halluces. The presentation in most of the affected patients was bilateral and symmetrical. A heterozygous C>A mutation at position 287 of the ZRS enhancer (chr7:156,584,283; hg19) was detected in all affected subjects and is absent from four unaffected family members, 42 unrelated samples, and multiple databases of human variation. Combined, these results identify a novel ZRS287 C>A mutation which leads to a variable spectrum of limb phenotypes.

Studies of TBX4 and chromosome 17q23.1q23.2: an uncommon cause of nonsyndromic clubfoot

Clubfoot is a common birth defect characterized by inward posturing and rigid downward displacement of one or both feet. The etiology of syndromic forms of clubfoot is varied and the causes of isolated clubfoot are not well understood. A microduplication of 2.2 Mb on chromosome 17q23.1q23.2 which includes T-box 4 (TBX4), a hindlimb-specific gene, and 16 other genes was recently identified in 3 of 66 families reported as nonsyndromic clubfoot, but additional non-foot malformations place them in the syndromic clubfoot category. Our study assesses whether variation in or around TBX4 contributes to nonsyndromic clubfoot. To determine whether this microduplication was a common cause of nonsyndromic clubfoot, 605 probands (from 148 multiplex and 457 simplex families) with nonsyndromic clubfoot were evaluated by copy number and oligonucleotide array CGH testing modalities. Only one multiplex family (0.68%) that had 16 with clubfoot and 9 with other foot anomalies, had a 350 kb microduplication, which included the complete duplication of TBX4 and NACA2 and partial duplication of BRIP1. The microduplication was transmitted in an autosomal dominant pattern and all with the microduplication had a range of phenotypes from short wide feet and toes to bilateral clubfoot. Minimal evidence was found for an association between TBX4 and clubfoot and no pathogenic sequence variants were identified in the two known TBX4 hindlimb enhancer elements. Altogether, these results demonstrate that variation in and around the TBX4 gene and the 17q23.1q23.2 microduplication are not a frequent cause of this common orthopedic birth defect and narrows the 17q23.1q23.2 nonsyndromic clubfoot-associated region.