A general approach for identifying distant regulatory elements applied to the Gdf6 gene

Transcriptional Interference Regulates the Evolutionary Development of Speech

The human capacity to speak is fundamental to our advanced intellectual, technological and social development. Yet so very little is known regarding the evolutionary genetics of speech or its relationship with the broader aspects of evolutionary development in primates. In this study, we describe a large family with evolutionary retrograde development of the larynx and wrist. The family presented with severe speech impairment and incremental retrograde elongations of the pisiform in the wrist that limited wrist rotation from 180° to 90° as in primitive primates. To our surprise, we found that a previously unknown primate-specific gene TOSPEAK had been disrupted in the family. TOSPEAK emerged de novo in an ancestor of extant primates across a 540 kb region of the genome with a pre-existing highly conserved long-range laryngeal enhancer for a neighbouring bone morphogenetic protein gene GDF6. We used transgenic mouse modelling to identify two additional GDF6 long-range enhancers within TOSPEAK that regulate GDF6 expression in the wrist. Disruption of TOSPEAK in the affected family blocked the transcription of TOSPEAK across the 3 GDF6 enhancers in association with a reduction in GDF6 expression and retrograde development of the larynx and wrist. Furthermore, we describe how TOSPEAK developed a human-specific promoter through the expansion of a penta-nucleotide direct repeat that first emerged de novo in the promoter of TOSPEAK in gibbon. This repeat subsequently expanded incrementally in higher hominids to form an overlapping series of Sp1/KLF transcription factor consensus binding sites in human that correlated with incremental increases in the promoter strength of TOSPEAK with human having the strongest promoter. Our research indicates a dual evolutionary role for the incremental increases in TOSPEAK transcriptional interference of GDF6 enhancers in the incremental evolutionary development of the wrist and larynx in hominids and the human capacity to speak and their retrogression with the reduction of TOSPEAK transcription in the affected family.

GDF6 Knockdown in a Family with Multiple Synostosis Syndrome and Speech Impairment

Multiple synostoses syndrome type 4 (SYNS4; MIM 617898) is an autosomal dominant disorder characterized by carpal-tarsal coalition and otosclerosis-associated hearing loss. SYSN4 has been associated with GDF6 gain-of-function mutations. Here we report a five-generation SYNS4 family with a reduction in GDF6 expression resulting from a chromosomal breakpoint 3' of GDF6. A 30-year medical history of the family indicated bilateral carpal-tarsal coalition in ~50% of affected family members and acquired otosclerosis-associated hearing loss in females only, whereas vertebral fusion was present in all affected family members, most of whom were speech impaired. All vertebral fusions were acquired postnatally in progressive fashion from a very early age. Thinning across the 2nd cervical vertebral interspace (C2-3) in the proband during infancy progressed to block fusion across C2-7 and T3-7 later in life. Carpal-tarsal coalition and pisiform expansion were bilaterally symmetrical within, but varied greatly between, affected family members. This is the first report of SYNS4 in a family with reduced GDF6 expression indicating a prenatal role for GDF6 in regulating development of the joints of the carpals and tarsals, the pisiform, ears, larynx, mouth and face and an overlapping postnatal role in suppression of aberrant ossification and synostosis of the joints of the inner ear (otosclerosis), larynx and vertebrae. RNAseq gene expression analysis indicated >10 fold knockdown of NOMO3, RBMXL1 and NEIL2 in both primary fibroblast cultures and fresh white blood cells. Together these results provide greater insight into the role of GDF6 in skeletal joint development.

siRNA Mediate RNA Interference Concordant with Early On-Target Transient Transcriptional Interference

Exogenous siRNAs are commonly used to regulate endogenous gene expression levels for gene function analysis, genotype–phenotype association studies and for gene therapy. Exogenous siRNAs can target mRNAs within the cytosol as well as nascent RNA transcripts within the nucleus, thus complicating siRNA targeting specificity. To highlight challenges in achieving siRNA target specificity, we targeted an overlapping gene set that we found associated with a familial form of multiple synostosis syndrome type 4 (SYSN4). In the affected family, we found that a previously unknown non-coding gene TOSPEAK/C8orf37AS1 was disrupted and the adjacent gene GDF6 was downregulated. Moreover, a conserved long-range enhancer for GDF6 was found located within TOSPEAK which in turn overlapped another gene which we named SMALLTALK/C8orf37. In fibroblast cell lines, SMALLTALK is transcribed at much higher levels in the opposite (convergent) direction to TOSPEAK. siRNA targeting of SMALLTALK resulted in post transcriptional gene silencing (PTGS/RNAi) of SMALLTALK that peaked at 72 h together with a rapid early increase in the level of both TOSPEAK and GDF6 that peaked and waned after 24 h. These findings indicated the following sequence of events: Firstly, the siRNA designed to target SMALLTALK mRNA for RNAi in the cytosol had also caused an early and transient transcriptional interference of SMALLTALK in the nucleus; Secondly, the resulting interference of SMALLTALK transcription increased the transcription of TOSPEAK; Thirdly, the increased transcription of TOSPEAK increased the transcription of GDF6. These findings have implications for the design and application of RNA and DNA targeting technologies including siRNA and CRISPR. For example, we used siRNA targeting of SMALLTALK to successfully restore GDF6 levels in the gene therapy of SYNS4 family fibroblasts in culture. To confidently apply gene targeting technologies, it is important to first determine the transcriptional interference effects of the targeting reagent and the targeted gene.

Differential transcriptome dynamics during the onset of conceptus elongation and between female and male porcine embryos

BACKGROUND

Porcine embryos undergo rapid differentiation and expansion between Days 8 and 12 before attaching to the maternal uterine epithelial surface after Day 13. It is known that maternal recognition of pregnancy and successful implantation are driven by mutual interactions between the elongated conceptus and the maternal endometrium. While most of the genes involved in regulation of embryo development are located on autosomal chromosomes, gene expression on sex chromosomes is modulating development through sex-specific transcription. To gain more insights into the dynamic transcriptome of preimplantation embryos at the onset of elongation and into X-linked gene expression, RNA-seq analyses were performed for single female and male porcine embryos collected on Days 8, 10, and 12 of pregnancy.

RESULTS

A high number of genes were differentially expressed across the developmental stages (2174 and 3275 for Days 8 vs 10, and 10 vs 12, respectively). The majority of differentially expressed genes (DEGs) were involved in embryo elongation, development, and embryo-maternal interaction. Interestingly, a number of DEGs was found with respect to embryo sex (137, 37, and 56 on Days 8, 10 and 12, respectively). At Day 8, most of these DEGs were X-linked (96). Strikingly, the number of DEGs encoded on the X chromosome dramatically decreased from Day 10 to Day 12.

CONCLUSIONS

The obtained results deepen the understanding about temporary transcriptomic changes in porcine embryos during the phase of conceptus elongation, meanwhile reveal dynamic compensation of X chromosome in the female and distinct transcriptional differences between female and male embryos.

An integrative genomic analysis of the Longshanks selection experiment for longer limbs in mice

Evolutionary studies are often limited by missing data that are critical to understanding the history of selection. Selection experiments, which reproduce rapid evolution under controlled conditions, are excellent tools to study how genomes evolve under selection. Here we present a genomic dissection of the Longshanks selection experiment, in which mice were selectively bred over 20 generations for longer tibiae relative to body mass, resulting in 13% longer tibiae in two replicates. We synthesized evolutionary theory, genome sequences and molecular genetics to understand the selection response and found that it involved both polygenic adaptation and discrete loci of major effect, with the strongest loci tending to be selected in parallel between replicates. We show that selection may favor de-repression of bone growth through inactivating two limb enhancers of an inhibitor, Nkx3-2. Our integrative genomic analyses thus show that it is possible to connect individual base-pair changes to the overall selection response.

Therapeutic potential of growth differentiation factors in the treatment of degenerative disc diseases

Intervertebral disc (IVD) degeneration is a major contributing factor to chronic low back pain and disability, leading to imbalance between anabolic and catabolic processes, altered extracellular matrix composition, loss of tissue hydration, inflammation, and impaired mechanical functionality. Current treatments aim to manage symptoms rather than treat underlying pathology. Therefore, IVD degeneration is a target for regenerative medicine strategies. Research has focused on understanding the molecular process of degeneration and the identification of various factors that may have the ability to halt and even reverse the degenerative process. One such family of growth factors, the growth differentiation factor (GDF) family, have shown particular promise for disc regeneration in in vitro and in vivo models of IVD degeneration. This review outlines our current understanding of IVD degeneration, and in this context, aims to discuss recent advancements in the use of GDF family members as anabolic factors for disc regeneration. An increasing body of evidence indicates that GDF family members are central to IVD homeostatic processes and are able to upregulate healthy nucleus pulposus cell marker genes in degenerative cells, induce mesenchymal stem cells to differentiate into nucleus pulposus cells and even act as chemotactic signals mobilizing resident cell populations during disc injury repair. The understanding of GDF signaling and its interplay with inflammatory and catabolic processes may be critical for the future development of effective IVD regeneration therapies.

Genetics of scapula and pelvis development: An evolutionary perspective

In tetrapods, the scapular and pelvic girdles perform the important function of anchoring the limbs to the trunk of the body and facilitating the movement of each appendage. This shared function, however, is one of relatively few similarities between the scapula and pelvis, which have significantly different morphologies, evolutionary histories, embryonic origins, and underlying genetic pathways. The scapula evolved in jawless fish prior to the pelvis, and its embryonic development is unique among bones in that it is derived from multiple progenitor cell populations, including the dermomyotome, somatopleure, and neural crest. Conversely, the pelvis evolved several million years later in jawed fish, and it develops from an embryonic somatopleuric cell population. The genetic networks controlling the formation of the pelvis and scapula also share similarities and differences, with a number of genes shaping only one or the other, while other gene products such as PBX transcription factors act as hierarchical developmental regulators of both girdle structures. Here, we provide a detailed review of the cellular processes and genetic networks underlying pelvis and scapula formation in tetrapods, while also highlighting unanswered questions about girdle evolution and development.

A novel enhancer near the Pitx1 gene influences development and evolution of pelvic appendages in vertebrates

Vertebrate pelvic reduction is a classic example of repeated evolution. Recurrent loss of pelvic appendages in sticklebacks has previously been linked to natural mutations in a pelvic enhancer that maps upstream of Pitx1. The sequence of this upstream PelA enhancer is not conserved to mammals, so we have surveyed a large region surrounding the mouse Pitx1 gene for other possible hind limb control sequences. Here we identify a new pelvic enhancer, PelB, that maps downstream rather than upstream of Pitx1. PelB drives expression in the posterior portion of the developing hind limb, and deleting the sequence from mice alters the size of several hind limb structures. PelB sequences are broadly conserved from fish to mammals. A wild stickleback population lacking the pelvis has an insertion/deletion mutation that disrupts the structure and function of PelB, suggesting that changes in this ancient enhancer contribute to evolutionary modification of pelvic appendages in nature.

Epigenetic profiling of growth plate chondrocytes sheds insight into regulatory genetic variation influencing height

GWAS have identified hundreds of height-associated loci. However, determining causal mechanisms is challenging, especially since height-relevant tissues (e.g. growth plates) are difficult to study. To uncover mechanisms by which height GWAS variants function, we performed epigenetic profiling of murine femoral growth plates. The profiled open chromatin regions recapitulate known chondrocyte and skeletal biology, are enriched at height GWAS loci, particularly near differentially expressed growth plate genes, and enriched for binding motifs of transcription factors with roles in chondrocyte biology. At specific loci, our analyses identified compelling mechanisms for GWAS variants. For example, at CHSY1, we identified a candidate causal variant (rs9920291) overlapping an open chromatin region. Reporter assays demonstrated that rs9920291 shows allelic regulatory activity, and CRISPR/Cas9 targeting of human chondrocytes demonstrates that the region regulates CHSY1 expression. Thus, integrating biologically relevant epigenetic information (here, from growth plates) with genetic association results can identify biological mechanisms important for human growth.

Ancient selection for derived alleles at a GDF5 enhancer influencing human growth and osteoarthritis risk

Variants in GDF5 are associated with human arthritis and decreased height, but the causal mutations are still unknown. We surveyed the Gdf5 locus for regulatory regions in transgenic mice and fine-mapped separate enhancers controlling expression in joints versus growing ends of long bones. A large downstream regulatory region contains a novel growth enhancer (GROW1), which is required for normal Gdf5 expression at ends of developing bones and for normal bone lengths in vivo. Human GROW1 contains a common base-pair change that decreases enhancer activity and colocalizes with peaks of positive selection in humans. The derived allele is rare in Africa but common in Eurasia and is found in Neandertals and Denisovans. Our results suggest that an ancient regulatory variant in GROW1 has been repeatedly selected in northern environments and that past selection on growth phenotypes explains the high frequency of a GDF5 haplotype that also increases arthritis susceptibility in many human populations.

Heads, Shoulders, Elbows, Knees, and Toes: Modular Gdf5 Enhancers Control Different Joints in the Vertebrate Skeleton

Synovial joints are crucial for support and locomotion in vertebrates, and are the frequent site of serious skeletal defects and degenerative diseases in humans. Growth and differentiation factor 5 (Gdf5) is one of the earliest markers of joint formation, is required for normal joint development in both mice and humans, and has been genetically linked to risk of common osteoarthritis in Eurasian populations. Here, we systematically survey the mouse Gdf5 gene for regulatory elements controlling expression in synovial joints. We identify separate regions of the locus that control expression in axial tissues, in proximal versus distal joints in the limbs, and in remarkably specific sub-sets of composite joints like the elbow. Predicted transcription factor binding sites within Gdf5 regulatory enhancers are required for expression in particular joints. The multiple enhancers that control Gdf5 expression in different joints are distributed over a hundred kilobases of DNA, including regions both upstream and downstream of Gdf5 coding exons. Functional rescue tests in mice confirm that the large flanking regions are required to restore normal joint formation and patterning. Orthologs of these enhancers are located throughout the large genomic region previously associated with common osteoarthritis risk in humans. The large array of modular enhancers for Gdf5 provide a new foundation for studying the spatial specificity of joint patterning in vertebrates, as well as new candidates for regulatory regions that may also influence osteoarthritis risk in human populations.

Joints, such as the hip and knee, are crucial for support and locomotion in animals, and are the frequent sites of serious human diseases such as arthritis. The Growth and differentiation factor 5 (Gdf5) gene is required for normal joint formation, and has been linked to risk of common arthritis in Eurasians. Here, we surveyed the mouse gene for the regulatory information that controls Gdf5's expression pattern in stripes at sites of joint formation. The gene does not have a single regulatory sequence that drives expression in all joints. Instead, Gdf5 has multiple different control sequences that show striking specificity for joints in the head, vertebral column, shoulder, elbow, wrist, hip, knee, and digits. Rescue experiments show that multiple control sequences are required to restore normal joint formation in Gdf5 mutants. The joint control sequences originally found in mice are also present in humans, where they are marked as active regions during fetal development and post-natal life, and map to a large region associated with arthritis risk in human populations. Regulatory variants in the human GDF5 control sequences can now be studied for their potential role in altering joint development or disease risk at particular locations in the skeleton.

Ror2 signaling is required for local upregulation of GDF6 and activation of BMP signaling at the neural plate border

The receptor tyrosine kinase Ror2 is a major Wnt receptor that activates β-catenin-independent signaling and plays a conserved role in the regulation of convergent extension movements and planar cell polarity in vertebrates. Mutations in the ROR2 gene cause recessive Robinow syndrome in humans, a short-limbed dwarfism associated with craniofacial malformations. Here, we show that Ror2 is required for local upregulation of gdf6 at the neural plate border in Xenopus embryos. Ror2 morphant embryos fail to upregulate neural plate border genes and show defects in the induction of neural crest cell fate. These embryos lack the spatially restricted activation of BMP signaling at the neural plate border at early neurula stages, which is required for neural crest induction. Ror2-dependent planar cell polarity signaling is required in the dorsolateral marginal zone during gastrulation indirectly to upregulate the BMP ligand Gdf6 at the neural plate border and Gdf6 is sufficient to rescue neural plate border specification in Ror2 morphant embryos. Thereby, Ror2 links Wnt/planar cell polarity signaling to BMP signaling in neural plate border specification and neural crest induction.

Microinjection for Transgenesis and Genome Editing in Threespine Sticklebacks

The threespine stickleback fish has emerged as a powerful system to study the genetic basis of a wide variety of morphological, physiological, and behavioral phenotypes. The remarkably diverse phenotypes that have evolved as marine populations adapt to countless freshwater environments, combined with the ability to cross marine and freshwater forms, provide a rare vertebrate system in which genetics can be used to map genomic regions controlling evolved traits. Excellent genomic resources are now available, facilitating molecular genetic dissection of evolved changes. While mapping experiments generate lists of interesting candidate genes, functional genetic manipulations are required to test the roles of these genes. Gene regulation can be studied with transgenic reporter plasmids and BACs integrated into the genome using the Tol2 transposase system. Functions of specific candidate genes and cis-regulatory elements can be assessed by inducing targeted mutations with TALEN and CRISPR/Cas9 genome editing reagents. All methods require introducing nucleic acids into fertilized one-cell stickleback embryos, a task made challenging by the thick chorion of stickleback embryos and the relatively small and thin blastomere. Here, a detailed protocol for microinjection of nucleic acids into stickleback embryos is described for transgenic and genome editing applications to study gene expression and function, as well as techniques to assess the success of transgenesis and recover stable lines.

Evolving New Skeletal Traits by cis-Regulatory Changes in Bone Morphogenetic Proteins

Changes in bone size and shape are defining features of many vertebrates. Here we use genetic crosses and comparative genomics to identify specific regulatory DNA alterations controlling skeletal evolution. Armor bone-size differences in sticklebacks map to a major effect locus overlapping BMP family member GDF6. Freshwater fish express more GDF6 due in part to a transposon insertion, and transgenic overexpression of GDF6 phenocopies evolutionary changes in armor-plate size. The human GDF6 locus also has undergone distinctive regulatory evolution, including complete loss of an enhancer that is otherwise highly conserved between chimps and other mammals. Functional tests show that the ancestral enhancer drives expression in hindlimbs but not forelimbs, in locations that have been specifically modified during the human transition to bipedalism. Both gain and loss of regulatory elements can localize BMP changes to specific anatomical locations, providing a flexible regulatory basis for evolving species-specific changes in skeletal form.

A distinct regulatory region of the Bmp5 locus activates gene expression following adult bone fracture or soft tissue injury

Bone morphogenetic proteins (BMPs) are key signaling molecules required for normal development of bones and other tissues. Previous studies have shown that null mutations in the mouse Bmp5 gene alter the size, shape and number of multiple bone and cartilage structures during development. Bmp5 mutations also delay healing of rib fractures in adult mutants, suggesting that the same signals used to pattern embryonic bone and cartilage are also reused during skeletal regeneration and repair. Despite intense interest in BMPs as agents for stimulating bone formation in clinical applications, little is known about the regulatory elements that control developmental or injury-induced BMP expression. To compare the DNA sequences that activate gene expression during embryonic bone formation and following acute injuries in adult animals, we assayed regions surrounding the Bmp5 gene for their ability to stimulate lacZ reporter gene expression in transgenic mice. Multiple genomic fragments, distributed across the Bmp5 locus, collectively coordinate expression in discrete anatomic domains during normal development, including in embryonic ribs. In contrast, a distinct regulatory region activated expression following rib fracture in adult animals. The same injury control region triggered gene expression in mesenchymal cells following tibia fracture, in migrating keratinocytes following dorsal skin wounding, and in regenerating epithelial cells following lung injury. The Bmp5 gene thus contains an "injury response" control region that is distinct from embryonic enhancers, and that is activated by multiple types of injury in adult animals.

Harnessing mobile genetic elements to explore gene regulation

Sequences that regulate expression of a gene in cis but are located at large distances along the DNA from the gene, as found with most developmentally regulated genes in higher vertebrates, are difficult to identify if those sequences are not conserved across species. Mutating suspected gene-regulatory sequences to alter expression then becomes a hit-or-miss affair. The relaxed specificity of transposon insertions offers an opportunity to develop alternate strategies, to scan in an unbiased manner, pieces of chromosomal DNA cloned in BACs for transcription enhancing elements. This article illustrates how insertions of Tn10 with enhancer-traps into BAC DNA containing the gene, and its germ-line expression in zebrafish, have identified distal regulatory elements functionally. Transposition of Tn10 first introduces the enhancer-trap with a loxP site randomly into BAC DNA. Cre-recombination between the inserted loxP and the loxP endogenous to a BAC-end positions the enhancer-trap to the newly created truncated end of BAC DNA. The procedure generates a library of integration-ready enhancer-trap BACs with progressive truncations from an end in a single experiment. Individual enhancer-trap BACs from the library can be evaluated functionally in zebrafish or mice. Furthermore, the ability to readily alter sequences in a small transposon plasmid containing a regulatory domain of the gene allows re-introduction of altered parts of a BAC back into itself. It serves as a useful strategy to functionally dissect multiple discontinuous regulatory domains of a gene quickly. These methodologies have been successfully used in identifying novel regulatory domains of the Amyloid Precursor Protein (appb) gene in zebrafish, and provided important clues for regulation of the gene in humans.

A molecular basis for classic blond hair color in Europeans

Hair color differences are among the most obvious examples of phenotypic variation in humans. Although genome-wide association studies (GWAS) have implicated multiple loci in human pigment variation, the causative base-pair changes are still largely unknown. Here we dissect a regulatory region of the KITLG gene (encoding KIT ligand) that is significantly associated with common blond hair color in northern Europeans. Functional tests demonstrate that the region contains a regulatory enhancer that drives expression in developing hair follicles. This enhancer contains a common SNP (rs12821256) that alters a binding site for the lymphoid enhancer-binding factor 1 (LEF1) transcription factor, reducing LEF1 responsiveness and enhancer activity in cultured human keratinocytes. Mice carrying ancestral or derived variants of the human KITLG enhancer exhibit significant differences in hair pigmentation, confirming that altered regulation of an essential growth factor contributes to the classic blond hair phenotype found in northern Europeans.

Hox11 genes are required for regional patterning and integration of muscle, tendon and bone

Development of the musculoskeletal system requires precise integration of muscles, tendons and bones. The molecular mechanisms involved in the differentiation of each of these tissues have been the focus of significant research; however, much less is known about how these tissues are integrated into a functional unit appropriate for each body position and role. Previous reports have demonstrated crucial roles for Hox genes in patterning the axial and limb skeleton. Loss of Hox11 paralogous gene function results in dramatic malformation of limb zeugopod skeletal elements, the radius/ulna and tibia/fibula, as well as transformation of the sacral region to a lumbar phenotype. Utilizing a Hoxa11eGFP knock-in allele, we show that Hox11 genes are expressed in the connective tissue fibroblasts of the outer perichondrium, tendons and muscle connective tissue of the zeugopod region throughout all stages of development. Hox11 genes are not expressed in differentiated cartilage or bone, or in vascular or muscle cells in these regions. Loss of Hox11 genes disrupts regional muscle and tendon patterning of the limb in addition to affecting skeletal patterning. The tendon and muscle defects in Hox11 mutants are independent of skeletal patterning events as disruption of tendon and muscle patterning is observed in Hox11 compound mutants that do not have a skeletal phenotype. Thus, Hox genes are not simply regulators of skeletal morphology as previously thought, but are key factors that regulate regional patterning and integration of the musculoskeletal system.

Identifying Distal cis-acting Gene-Regulatory Sequences by Expressing BACs Functionalized with loxP-Tn10 Transposons in Zebrafish

Bacterial Artificial Chromosomes (BACs) are large pieces of DNA from the chromosomes of organisms propagated faithfully in bacteria as large extra-chromosomal plasmids. Expression of genes contained in BACs can be monitored after functionalizing the BAC DNA with reporter genes and other sequences that allow stable maintenance and propagation of the DNA in the new host organism. The DNA in BACs can be altered within its bacterial host in several ways. Here we discuss one such approach, using Tn10 mini-transposons, to introduce exogenous sequences into BACs for a variety of purposes. The largely random insertions of Tn10 transposons carrying lox sites have been used to position mammalian cell-selectable antibiotic resistance genes, enhancer-traps and inverted repeat ends of the vertebrate transposon Tol2 precisely at the ends of the genomic DNA insert in BACs. These modified BACs are suitable for expression in zebrafish or mouse, and have been used to functionally identify important long-range gene regulatory sequences in both species. Enhancer-trapping using BACs should prove uniquely useful in analyzing multiple discontinuous DNA domains that act in concert to regulate expression of a gene, and is not limited by genome accessibility issues of traditional enhancer-trapping methods.

Stratified whole genome linkage analysis of Chiari type I malformation implicates known Klippel-Feil syndrome genes as putative disease candidates

Chiari Type I Malformation (CMI) is characterized by displacement of the cerebellar tonsils below the base of the skull, resulting in significant neurologic morbidity. Although multiple lines of evidence support a genetic contribution to disease, no genes have been identified. We therefore conducted the largest whole genome linkage screen to date using 367 individuals from 66 families with at least two individuals presenting with nonsyndromic CMI with or without syringomyelia. Initial findings across all 66 families showed minimal evidence for linkage due to suspected genetic heterogeneity. In order to improve power to localize susceptibility genes, stratified linkage analyses were performed using clinical criteria to differentiate families based on etiologic factors. Families were stratified on the presence or absence of clinical features associated with connective tissue disorders (CTDs) since CMI and CTDs frequently co-occur and it has been proposed that CMI patients with CTDs represent a distinct class of patients with a different underlying disease mechanism. Stratified linkage analyses resulted in a marked increase in evidence of linkage to multiple genomic regions consistent with reduced genetic heterogeneity. Of particular interest were two regions (Chr8, Max LOD = 3.04; Chr12, Max LOD = 2.09) identified within the subset of "CTD-negative" families, both of which harbor growth differentiation factors (GDF6, GDF3) implicated in the development of Klippel-Feil syndrome (KFS). Interestingly, roughly 3-5% of CMI patients are diagnosed with KFS. In order to investigate the possibility that CMI and KFS are allelic, GDF3 and GDF6 were sequenced leading to the identification of a previously known KFS missense mutation and potential regulatory variants in GDF6. This study has demonstrated the value of reducing genetic heterogeneity by clinical stratification implicating several convincing biological candidates and further supporting the hypothesis that multiple, distinct mechanisms are responsible for CMI.

Manipulating the Mouse Genome Using Recombineering

Genetically engineered mouse models are indispensable for understanding the biological function of genes, understanding the genetic basis of human diseases and for preclinical testing of novel therapies. Generation of such mouse models has been possible because of our ability to manipulate the mouse genome. Recombineering is a highly efficient recombination-based method of genetic engineering that has revolutionized our ability to generate mouse models. Since recombineering technology is not dependent on the availability of restriction enzyme recognition sites, it allows us to modify the genome with great precision. It requires homology arms as short as 40 bases for recombination, which makes it relatively easy to generate targeting constructs to insert, change or delete either a single nucleotide or a DNA fragment several kb in size; insert selectable markers, reporter genes or add epitope tags to any gene of interest. In this review, we focus on the development of recombineering technology and its application in the generation of transgenic and knockout or knock-in mouse models. High throughput generation of gene targeting vectors, used to construct knockout alleles in mouse embryonic stem cells, is now feasible because of this technology. The challenge now is to use the "designer" mice to develop novel therapies to prevent, cure or effectively manage some the most debilitating human diseases.

A misplaced lncRNA causes brachydactyly in humans

Translocations are chromosomal rearrangements that are frequently associated with a variety of disease states and developmental disorders. We identified 2 families with brachydactyly type E (BDE) resulting from different translocations affecting chromosome 12p. Both translocations caused downregulation of the parathyroid hormone-like hormone (PTHLH) gene by disrupting the cis-regulatory landscape. Using chromosome conformation capturing, we identified a regulator on chromosome 12q that interacts in cis with PTHLH over a 24.4-megabase distance and in trans with the sex-determining region Y-box 9 (SOX9) gene on chromosome 17q. The element also harbored a long noncoding RNA (lncRNA). Silencing of the lncRNA, PTHLH, or SOX9 revealed a feedback mechanism involving an expression-dependent network in humans. In the BDE patients, the human lncRNA was upregulated by the disrupted chromosomal association. Moreover, the lncRNA occupancy at the PTHLH locus was reduced. Our results document what we believe to be a novel in cis- and in trans-acting DNA and lncRNA regulatory feedback element that is reciprocally regulated by coding genes. Furthermore, our findings provide a systematic and combinatorial view of how enhancers encoding lncRNAs may affect gene expression in normal development.

Vascular tissues are a primary source of BMP2 expression during bone formation induced by distraction osteogenesis

Prior studies showed that bone regeneration during distraction osteogenesis (DO) was dependent on vascular tissue development and that inhibition of VEGFR signaling diminished the expression of BMP2. A combination of micro-computed tomography (μCT) analysis of vascular and skeletal tissues, immunohistological and histological analysis of transgenic mice containing a BAC transgene in which β-galactosidase had been inserted into the coding region of BMP2 and qRT-PCR analysis, was used to examine how the spatial temporal expression of the morphogenetic signals that drive skeletal and vascular tissue development is coordinated during DO. These results showed that BMP2 expression was induced in smooth muscle and vascular endothelial cells of arteries and veins, capillary endothelial cells, hypertrophic chondrocytes and osteocytes. BMP2 was not expressed by lymphatic vessels or macrophages. Separate peaks of BMP2 mRNA expression were induced in the surrounding muscular tissues and the distraction gap and corresponded first with large vessel collateralization and arteriole remodeling followed by periods of angiogenesis in the gap region. Immunohistological and qRT-PCR analysis of VEGF receptors and ligands showed that mesenchymal cells, lining cells and chondrocytes, expressed VEGFA, although PlGF expression was only seen in mesenchymal cells within the gap region. On the other hand VEGFR2 appeared to be predominantly expressed by vascular endothelial and hematopoietic cells. These results suggest that bone and vascular tissue formation is coordinated via a mutually supporting set of paracrine loops in which blood vessels primarily synthesize the morphogens that promote bone formation while mesenchymal cells primarily synthesize the morphogens that promote vascular tissue formation.

The BMP ligand Gdf6 prevents differentiation of coronal suture mesenchyme in early cranial development

Growth Differentiation Factor-6 (Gdf6) is a member of the Bone Morphogenetic Protein (BMP) family of secreted signaling molecules. Previous studies have shown that Gdf6 plays a role in formation of a diverse subset of skeletal joints. In mice, loss of Gdf6 results in fusion of the coronal suture, the intramembranous joint that separates the frontal and parietal bones. Although the role of GDFs in the development of cartilaginous limb joints has been studied, limb joints are developmentally quite distinct from cranial sutures and how Gdf6 controls suture formation has remained unclear. In this study we show that coronal suture fusion in the Gdf6-/- mouse is due to accelerated differentiation of suture mesenchyme, prior to the onset of calvarial ossification. Gdf6 is expressed in the mouse frontal bone primordia from embryonic day (E) 10.5 through 12.5. In the Gdf6-/- embryo, the coronal suture fuses prematurely and concurrently with the initiation of osteogenesis in the cranial bones. Alkaline phosphatase (ALP) activity and Runx2 expression assays both showed that the suture width is reduced in Gdf6+/- embryos and is completely absent in Gdf6-/- embryos by E12.5. ALP activity is also increased in the suture mesenchyme of Gdf6+/- embryos compared to wild-type. This suggests Gdf6 delays differentiation of the mesenchyme occupying the suture, prior to the onset of ossification. Therefore, although BMPs are known to promote bone formation, Gdf6 plays an inhibitory role to prevent the osteogenic differentiation of the coronal suture mesenchyme.

Mutation of Rubie, a novel long non-coding RNA located upstream of Bmp4, causes vestibular malformation in mice

Background

The vestibular apparatus of the vertebrate inner ear uses three fluid-filled semicircular canals to sense angular acceleration of the head. Malformation of these canals disrupts the sense of balance and frequently causes circling behavior in mice. The Epistatic circler (Ecl) is a complex mutant derived from wildtype SWR/J and C57L/J mice. Ecl circling has been shown to result from the epistatic interaction of an SWR-derived locus on chromosome 14 and a C57L-derived locus on chromosome 4, but the causative genes have not been previously identified.

Methodology/Principal Findings

We developed a mouse chromosome substitution strain (CSS-14) that carries an SWR/J chromosome 14 on a C57BL/10J genetic background and, like Ecl, exhibits circling behavior due to lateral semicircular canal malformation. We utilized CSS-14 to identify the chromosome 14 Ecl gene by positional cloning. Our candidate interval is located upstream of bone morphogenetic protein 4 (Bmp4) and contains an inner ear-specific, long non-coding RNA that we have designated Rubie (RNA upstream of Bmp4 expressed in inner ear). Rubie is spliced and polyadenylated, and is expressed in developing semicircular canals. However, we discovered that the SWR/J allele of Rubie is disrupted by an intronic endogenous retrovirus that causes aberrant splicing and premature polyadenylation of the transcript. Rubie lies in the conserved gene desert upstream of Bmp4, within a region previously shown to be important for inner ear expression of Bmp4. We found that the expression patterns of Bmp4 and Rubie are nearly identical in developing inner ears.

Conclusions/Significance

Based on these results and previous studies showing that Bmp4 is essential for proper vestibular development, we propose that Rubie is the gene mutated in Ecl mice, that it is involved in regulating inner ear expression of Bmp4, and that aberrant Bmp4 expression contributes to the Ecl phenotype.

High temperature requirement factor A1 (HTRA1) gene regulates angiogenesis through transforming growth factor-β family member growth differentiation factor 6

Genome-wide association study (GWAS) has identified genetic variants in the promoter region of the high temperature requirement factor A1 (HTRA1) gene associated with age-related macular degeneration (AMD). As a secreted serine protease, HTRA1 has been reported to interact with members of the transforming growth factor-β (TGF-β) family and regulate their signaling pathways. Growth differentiation factor 6 (GDF6), a member of the TGF-β family, is involved in ectoderm patterning and eye development. Mutations in GDF6 have been associated with abnormal eye development that may result in microphthalmia and anophthalmia. In this report, we identified a single nucleotide polymorphism (SNP) rs6982567 A/G near the GDF6 gene that is significantly associated with AMD (p value = 3.54 × 10(-8)). We demonstrated that the GDF6 AMD risk allele (rs6982567 A) is associated with decreased expression of the GDF6 and increased expression of HTRA1. Similarly, the HTRA1 AMD risk allele (rs10490924 T) is associated with decreased GDF6 and increased HTRA1 expression. We observed decreased vascular development in the retina and significant up-regulation of GDF6 gene in the RPE layer, retinal and brain tissues in HTRA1 knock-out (htra1(-/-)) mice as compared with the wild-type counterparts. Furthermore, we showed enhanced SMAD signaling in htra1(-/-) mice. Our data suggests a critical role of HTRA1 in the regulation of angiogenesis via TGF-β signaling and identified GDF6 as a novel disease gene for AMD.

Conserved and non-conserved enhancers direct tissue specific transcription in ancient germ layer specific developmental control genes

BACKGROUND

Identifying DNA sequences (enhancers) that direct the precise spatial and temporal expression of developmental control genes remains a significant challenge in the annotation of vertebrate genomes. Locating these sequences, which in many cases lie at a great distance from the transcription start site, has been a major obstacle in deciphering gene regulation. Coupling of comparative genomics with functional validation to locate such regulatory elements has been a successful method in locating many such regulatory elements. But most of these studies looked either at a single gene only or the whole genome without focusing on any particular process. The pressing need is to integrate the tools of comparative genomics with knowledge of developmental biology to validate enhancers for developmental transcription factors in greater detail

RESULTS

Our results show that near four different genes (nkx3.2, pax9, otx1b and foxa2) in zebrafish, only 20-30% of highly conserved DNA sequences can act as developmental enhancers irrespective of the tissue the gene expresses in. We find that some genes also have multiple conserved enhancers expressing in the same tissue at the same or different time points in development. We also located non-conserved enhancers for two of the genes (pax9 and otx1b). Our modified Bacterial artificial chromosome (BACs) studies for these 4 genes revealed that many of these enhancers work in a synergistic fashion, which cannot be captured by individual DNA constructs and are not conserved at the sequence level. Our detailed biochemical and transgenic analysis revealed Foxa1 binds to the otx1b non-conserved enhancer to direct its activity in forebrain and otic vesicle of zebrafish at 24 hpf.

CONCLUSION

Our results clearly indicate that high level of functional conservation of genes is not necessarily associated with sequence conservation of its regulatory elements. Moreover certain non conserved DNA elements might have role in gene regulation. The need is to bring together multiple approaches to bear upon individual genes to decipher all its regulatory elements.

Fishing for function: zebrafish BAC transgenics for functional genomics

Transgenics using bacterial artificial chromosomes (BACs) offers a great opportunity to look at gene regulation in a developing embryo. The modified BAC containing a reporter inserted just before the translational start site of the gene of interest allows for the visualization of spatio-temporal gene expression. Though this method has been used in the mouse model extensively, its utility in zebrafish studies is relatively new. This review aims to look at the utility of making BAC transgenics in zebrafish and its applications in functional genomics. We look at the various methods to modify the BAC, some limitations and what the future holds.

Identification of a distant cis-regulatory element controlling pharyngeal arch-specific expression of zebrafish gdf6a/radar

Skeletal formation is an essential and intricately regulated part of vertebrate development. Humans and mice deficient in growth and differentiation factor 6 (Gdf6) have numerous skeletal abnormalities, including joint fusions and cartilage reductions. The expression of Gdf6 is dynamic and in part regulated by distant evolutionarily conserved cis-regulatory elements. radar/gdf6a is a zebrafish ortholog of Gdf6 and has an essential role in embryonic patterning. Here, we show that radar is transcribed in the cells surrounding and between the developing cartilages of the ventral pharyngeal arches, similar to mouse Gdf6. A 312 bp evolutionarily conserved region (ECR5), 122 kilobases downstream, drives expression in a pharyngeal arch-specific manner similar to endogenous radar/gdf6a. Deletion analysis identified a 78 bp region within ECR5 that is essential for transgene activity. This work illustrates that radar is regulated in the pharyngeal arches by a distant conserved element and suggests radar has similar functions in skeletal development in fish and mammals.

Distal regions of the human IFNG locus direct cell type-specific expression

Genes, such as IFNG, which are expressed in multiple cell lineages of the immune system, may employ a common set of regulatory elements to direct transcription in multiple cell types or individual regulatory elements to direct expression in individual cell lineages. By employing a bacterial artificial chromosome transgenic system, we demonstrate that IFNG employs unique regulatory elements to achieve lineage-specific transcriptional control. Specifically, a one 1-kb element 30 kb upstream of IFNG activates transcription in T cells and NKT cells but not in NK cells. This distal regulatory element is a Runx3 binding site in Th1 cells and is needed for RNA polymerase II recruitment to IFNG, but it is not absolutely required for histone acetylation of the IFNG locus. These results support a model whereby IFNG uses cis-regulatory elements with cell type-restricted function.

Control of BMP gene expression by long-range regulatory elements

Much evidence suggests that "developmental regulator" genes, like those encoding transcription factors and signaling molecules, are typically controlled by many modular, tissue-specific cis-regulatory elements that function during embryogenesis. These elements are often far from gene coding regions and promoters. Bone morphogenetic proteins (BMPs) drive many processes in development relating to organogenesis and differentiation. Four BMP family members, Bmp2, Bmp4, Bmp5, and Gdf6, are now known to be under the control of distant cis-regulatory elements. BMPs are thus firmly placed in the category of genes prone to this phenomenon. The analysis of distant BMP regulatory elements has provided insight into the many pleiotropic effects of BMP genes, and underscores the biological importance of non-coding genomic DNA elements.

BMP2 is essential for post natal osteogenesis but not for recruitment of osteogenic stem cells

The effects of BMP2 on bone marrow stromal cell differentiation and bone formation after bone marrow ablation were determined using C57 BL/6J (B6) mice. Inhibition of BMP2 expression with lentiviral BMP2 shRNA prevented both mineralized nodule formation in vitro and bone formation in vivo, and blocked the expression of Runx2 and osterix, transcriptional determinants of terminal osteogenic differentiation. No effect was observed on the expression of Sox9, a transcription factor, which is the one of the first transcriptional determinant to be expressed in committed chondroprogenitor and osteoprogenitor cells. In vitro studies showed that exogenously added BMP7 rescued the expression of osterix and enhanced the expression of Sox9, but had no effect on the expression of Runx2, while it only partially recovered the development of mineral deposition in the cultures. On the other hand, the exogenous addition of BMP2 rescued both Runx2 and osterix expression, did not enhance the expression of Sox9, but fully recovered the inhibition of mineral deposition in the cultures. Using antibodies against CD146 and Sox9, immunohistological examination of the cell populations found in the medullary space three days after bone marrow ablation, showed qualitatively equal numbers of cells expressing these skeletal progenitor and stem cell markers in control and BMP2 shRNA treated animals. Fluorescence Activated Cell Sorting (FACS) analysis of the cells found with the marrow cavities at three days after marrow ablation using CD146 antibody showed near equal numbers of immunopositive cells in both control and shRNA treated animals. In summary, the differences observed in vitro for BMP2 and BMP7 on osteogenic gene expression and mineralization suggest that they have differing effects on bone cell differentiation. These results further demonstrate that in vivo BMP2 is a central morphogenetic regulator of post natal osteoprogenitor differentiation, but does not affect recruitment of progenitors to the osteoblastic lineage.

Duplications involving a conserved regulatory element downstream of BMP2 are associated with brachydactyly type A2

Autosomal-dominant brachydactyly type A2 (BDA2), a limb malformation characterized by hypoplastic middle phalanges of the second and fifth fingers, has been shown to be due to mutations in the Bone morphogenetic protein receptor 1B (BMPR1B) or in its ligand Growth and differentiation factor 5 (GDF5). A linkage analysis performed in a mutation-negative family identified a novel locus for BDA2 on chromosome 20p12.3 that incorporates the gene for Bone morphogenetic protein 2 (BMP2). No point mutation was identified in BMP2, so a high-density array CGH analysis covering the critical interval of approximately 1.3 Mb was performed. A microduplication of approximately 5.5 kb in a noncoding sequence approximately 110 kb downstream of BMP2 was detected. Screening of other patients by qPCR revealed a similar duplication in a second family. The duplicated region contains evolutionary highly conserved sequences suggestive of a long-range regulator. By using a transgenic mouse model we can show that this sequence is able to drive expression of a X-Gal reporter construct in the limbs. The almost complete overlap with endogenous Bmp2 expression indicates that a limb-specific enhancer of Bmp2 is located within the identified duplication. Our results reveal an additional functional mechanism for the pathogenesis of BDA2, which is duplication of a regulatory element that affects the expression of BMP2 in the developing limb.

Identification of an ancient Bmp4 mesoderm enhancer located 46 kb from the promoter

Bone morphogenetic protein 4 (Bmp4) is a multi-functional, developmentally regulated gene that is essential for mouse development, as most Bmp4-null mouse embryos die at the onset of gastrulation and fail to develop mesoderm. Little is known about the transcriptional regulation of Bmp4. To identify potential long-range cis-regulatory elements that direct its complex spatiotemporal expression patterns, we surveyed the mouse Bmp4 locus using two overlapping bacterial artificial chromosome (BAC) reporter transgenes. Our findings indicate that tissue-specific cis-regulatory elements reside greater than 28 kb 5' or 3' to the mouse Bmp4 transcription unit. In addition, comparative analyses identified three noncoding evolutionarily conserved regions (ECRs), spaced around the gene and conserved from mammals to fish, that are maintained in a syntenic group across vertebrates. Deletion of one of these conserved sequences (ECR2) from a BAC transgene revealed a tissue-specific requirement for ECR2 in driving Bmp4 expression in extraembryonic and embryonic mesoderm. Furthermore, a 467 bp mouse sequence containing ECR2 reproducibly directed lacZ minigene expression in mesoderm. Taken together, this shows that an ancient, mesoderm-specific cis-regulatory element resides nearly 50 kb 5' to mouse Bmp4.

Shaping skeletal growth by modular regulatory elements in the Bmp5 gene

Cartilage and bone are formed into a remarkable range of shapes and sizes that underlie many anatomical adaptations to different lifestyles in vertebrates. Although the morphological blueprints for individual cartilage and bony structures must somehow be encoded in the genome, we currently know little about the detailed genomic mechanisms that direct precise growth patterns for particular bones. We have carried out large-scale enhancer surveys to identify the regulatory architecture controlling developmental expression of the mouse Bmp5 gene, which encodes a secreted signaling molecule required for normal morphology of specific skeletal features. Although Bmp5 is expressed in many skeletal precursors, different enhancers control expression in individual bones. Remarkably, we show here that different enhancers also exist for highly restricted spatial subdomains along the surface of individual skeletal structures, including ribs and nasal cartilages. Transgenic, null, and regulatory mutations confirm that these anatomy-specific sequences are sufficient to trigger local changes in skeletal morphology and are required for establishing normal growth rates on separate bone surfaces. Our findings suggest that individual bones are composite structures whose detailed growth patterns are built from many smaller lineage and gene expression domains. Individual enhancers in BMP genes provide a genomic mechanism for controlling precise growth domains in particular cartilages and bones, making it possible to separately regulate skeletal anatomy at highly specific locations in the body.

Every bone in the skeleton has a specific shape and size. These characteristic features must be under separate genetic control, because individual bones can undergo striking morphological changes in different species. Researchers have long postulated that the morphology of individual bones arises from the local activity of many separate growth domains around each bone's surface. Differential growth within such domains could modify size, curvature, and formation of specific processes. Here, we show that local growth domains around individual bones are controlled by independent regulatory sequences in bone morphogenetic protein (BMP) genes. We identify multiple regulatory sequences in the Bmp5 gene that control expression in particular bones, rather than all bones. We show that some of these elements are remarkably specific for individual subdomains around the surface of individual bones. Finally, we show that local BMP signaling is necessary and sufficient to trigger highly localized growth patterns in ribs and nasal cartilages. These results suggest that the detailed pattern of growth of individual skeletal structures is encoded in part by multiple regulatory sequences in BMP genes. Gain and loss of anatomy-specific sequences in BMP genes may provide a flexible genomic mechanism for modifying local skeletal anatomy during vertebrate evolution.

Assaying the regulatory potential of mammalian conserved non-coding sequences in human cells

The fraction of experimentally active conserved non-coding sequences within any given cell type is low, so classical assays are unlikely to expose their potential.

Background

Conserved non-coding sequences in the human genome are approximately tenfold more abundant than known genes, and have been hypothesized to mark the locations of cis-regulatory elements. However, the global contribution of conserved non-coding sequences to the transcriptional regulation of human genes is currently unknown. Deeply conserved elements shared between humans and teleost fish predominantly flank genes active during morphogenesis and are enriched for positive transcriptional regulatory elements. However, such deeply conserved elements account for <1% of the conserved non-coding sequences in the human genome, which are predominantly mammalian.

Results

We explored the regulatory potential of a large sample of these 'common' conserved non-coding sequences using a variety of classic assays, including chromatin remodeling, and enhancer/repressor and promoter activity. When tested across diverse human model cell types, we find that the fraction of experimentally active conserved non-coding sequences within any given cell type is low (approximately 5%), and that this proportion increases only modestly when considered collectively across cell types.

Conclusions

The results suggest that classic assays of cis-regulatory potential are unlikely to expose the functional potential of the substantial majority of mammalian conserved non-coding sequences in the human genome.

Germline translocations in mice: unique tools for analyzing gene function and long-distance regulatory mechanisms

Translocations have provided invaluable tools for identifying both cancer-linked genes and loci associated with heritable human diseases, but heritable human translocations are rare and few mouse models exist. Here we report progress on analysis of a collection of heritable translocations generated by treatment of mice with specific chemicals or radiation during late spermatogenic stages. The translocation mutants exhibit a range of visible phenotypes reflecting the disruption of coding sequences or the separation of genes from essential regulatory elements. The breakpoints of both radiation-induced and chemically induced mutations in these mice are remarkably clean, with very short deletions, duplications, or inversions in some cases, and ligation mediated by microhomology, suggesting nonhomologous end joining as the major path of repair. These mutations provide new tools for the discovery of novel genes and regulatory elements linked to human developmental disorders and new clues to the molecular basis of human genetic disease.

BMP-2 vs. BMP-4 expression and activity in glucocorticoid-arrested MC3T3-E1 osteoblasts: Smad signaling, not alkaline phosphatase activity, predicts rescue of mineralization

Pharmacological glucocorticoids (GCs) inhibit bone formation, leading to osteoporosis. GCs inhibit bone morphogenetic protein-2 (Bmp2) expression, and rhBMP-2 restores mineralization in GC-arrested osteoblast cultures. To better understand how GCs regulate BMPs, we investigated Bmp transcription, as well as rhBMP-induced Smad and alkaline phosphatase (ALP) activity. Bmp2 cis-regulatory regions were analyzed by reporter plasmids and LacZ-containing bacterial artificial chromosomes. We found that GCs inhibited Bmp2 via a domain > 50 kb downstream of the coding sequence. Bmp expression was evaluated by RT-PCR; whereas GCs strongly inhibited Bmp2, Bmp4 was abundantly expressed and resistant to GCs. Both rhBMP-2 and rhBMP-4 restored mineralization in GC-arrested cultures; rhBMP-2 was 5-fold more effective when dosing was based on ALP activation, however, the rhBMPs were equipotent when dosing was based on Smad transactivation. In conclusion, GCs regulate Bmp2 via a far-downstream domain, and activation of Smad, not ALP, best predicts the pro-mineralization potential of rhBMPs.

Dual hindlimb control elements in the Tbx4 gene and region-specific control of bone size in vertebrate limbs

The Tbx4 transcription factor is crucial for normal hindlimb and vascular development, yet little is known about how its highly conserved expression patterns are generated. We have used comparative genomics and functional scanning in transgenic mice to identify a dispersed group of enhancers controlling Tbx4 expression in different tissues. Two independent enhancers control hindlimb expression, one located upstream and one downstream of the Tbx4 coding exons. These two enhancers, hindlimb enhancer A and hindlimb enhancer B (HLEA and HLEB), differ in their primary sequence, in their precise patterns of activity within the hindlimb, and in their degree of sequence conservation across animals. HLEB is highly conserved from fish to mammals. Although Tbx4 expression and hindlimb development occur at different axial levels in fish and mammals, HLEB cloned from either fish or mouse is capable of driving expression at the appropriate position of hindlimb development in mouse embryos. HLEA is highly conserved only in mammals. Deletion of HLEA from the endogenous mouse locus reduces expression of Tbx4 in the hindlimb during embryogenesis, bypasses the embryonic lethality of Tbx4-null mutations, and produces viable, fertile mice with characteristic changes in the size of bones in the hindlimb but not the forelimb. We speculate that dual hindlimb enhancers provide a flexible genomic mechanism for altering the strength and location of Tbx4 expression during normal development, making it possible to separately modify the size of forelimb and hindlimb bones during vertebrate evolution.

Misexpression of Six2 is associated with heritable frontonasal dysplasia and renal hypoplasia in 3H1 Br mice

A radiation-induced mouse mutant, Brachyrrhine (Br), exhibits frontonasal dysplasia and renal hypoplasia, two malformations associated with deficiencies in mesenchymal condensation. The purpose of this study was to resolve the Br locus, evaluate possible candidate genes, and identify developmental defects in the mutant chondrocranium. Linkage analysis mapped the Br mutation to a critical region distal to D17Mit76, which contains only one gene, the transcription factor Six2. Sequence analysis of the Six2 gene, including 1.5 kb of the promoter, failed to reveal the Br mutation. However, homozygous Br/Br embryos showed almost complete absence of Six2 mRNA and protein in craniofacial and renal tissues while heterozygous Br/+ embryos displayed intermediate Six2 levels. Mutant embryos displayed malformations of neural crest-derived structures of the anterior cranium where Six2 is normally expressed. These data suggest a mutation in a novel cis-acting regulatory region inhibits Six2 expression and is associated with frontonasal dysplasia and renal hypoplasia.

Searching for functional genetic variants in non-coding DNA

1. The search for DNA sequence variants for complex human polygenic conditions has been a strong focus of recent genetic research. While gene loci have been identified, few variants in the coding sequences of these genes have been found, suggesting that non-coding sequence variation may underlie many complex conditions. 2. Non-coding DNA harbours regulatory elements capable of making changes to gene expression. However, regulatory DNA sequences are currently difficult to recognize and their function is poorly understood, complicating the task of assigning potential functional significance to non-coding variation. 3. Comparative genomics, the study of evolutionary DNA conservation, has enabled the emergent field of non-coding DNA identification in human disease analysis. 4. This brief review will focus on the potential of a relatively high throughput technique based on comparative genomics, that may aid in the identification of functionally important non-coding sequence variation in complex diseases.

Approaches to comparative sequence analysis: towards a functional view of vertebrate genomes

The comparison of genomic sequences is now a common approach to identifying and characterizing functional regions in vertebrate genomes. However, for theoretical reasons and because of practical issues, the generation of these data sets is non-trivial and can have many pitfalls. We are currently seeing an explosion of comparative sequence data, the benefits and limitations of which need to be disseminated to the scientific community. This Review provides a critical overview of the different types of sequence data that are available for analysis and of contemporary comparative sequence analysis methods, highlighting both their strengths and limitations. Approaches to determining the biological significance of constrained sequence are also explored.

Genomic identification of regulatory elements by evolutionary sequence comparison and functional analysis

Despite remarkable recent advances in genomics that have enabled us to identify most of the genes in the human genome, comparable efforts to define transcriptional cis-regulatory elements that control gene expression are lagging behind. The difficulty of this task stems from two equally important problems: our knowledge of how regulatory elements are encoded in genomes remains elementary, and there is a vast genomic search space for regulatory elements, since most of mammalian genomes are noncoding. Comparative genomic approaches are having a remarkable impact on the study of transcriptional regulation in eukaryotes and currently represent the most efficient and reliable methods of predicting noncoding sequences likely to control the patterns of gene expression. By subjecting eukaryotic genomic sequences to computational comparisons and subsequent experimentation, we are inching our way toward a more comprehensive catalog of common regulatory motifs that lie behind fundamental biological processes. We are still far from comprehending how the transcriptional regulatory code is encrypted in the human genome and providing an initial global view of regulatory gene networks, but collectively, the continued development of comparative and experimental approaches will rapidly expand our knowledge of the transcriptional regulome.

Analysis of mouse Cdh6 gene regulation by transgenesis of modified bacterial artificial chromosomes

Classic cadherins are cell adhesion molecules whose expression patterns are dynamically modulated in association with their diverse functions during morphogenesis. The large size and complexity of cadherin loci have made it a challenge to investigate the organization of cis-regulatory modules that control their spatiotemporal patterns of expression. Towards this end, we utilized bacterial artificial chromosomes (BACs) containing the Cdh6 gene, a mouse type II classic cadherin, to systematically identify cis-regulatory modules that govern its expression. By inserting a lacZ reporter gene into the Cdh6 BAC and generating a series of modified variants via homologous recombination or transposon insertions that have been examined in transgenic mice, we identified an array of genomic regions that contribute to specific regulation of the gene. These regions span approximately 350 kb of the locus between 161-kb upstream and 186-kb downstream of the Cdh6 transcription start site. Distinct modules independently regulate compartmental expression (i.e. forebrain, hindbrain rhombomeres, and spinal cord) and/or cell lineage-specific expression patterns (i.e. neural crest subpopulations such as Schwann cells) of Cdh6 at the early developmental stages. With respect to regulation of expression in neural crest cells, we have found that distinct regions contribute to different aspects of expression and have identified a short 79-bp region that is implicated in regulating expression in cells once they have emigrated from the neural tube. These results build a picture of the complex organization of Cdh6 cis-regulatory modules and highlight the diverse inputs that contribute to its dynamic expression during early mouse embryonic development.

A mouse transgene drives embryonic dorsal posterior commissure expression

In this report we generated mice co-transgenic for a minimal promoter LacZ construct and a mouse BAC from the gene poor region upstream of the Hoxd cluster. In addition to expression in the distal limb, genital bud, and spinal cord, we show that this BAC transgene also reproducibly drives unique bilateral, dorsal posterior commissure expression. The ability of this BAC to direct posterior commissure expression makes it worthy of further study as a valuable tool in transgenic/targeting experiments.

Characterization of a familial t(16;22) balanced translocation associated with congenital cataract leads to identification of a novel gene, TMEM114, expressed in the lens and disrupted by the translocation

Molecular characterization of chromosomal rearrangements is a powerful resource in identification of genes associated with monogenic disorders. We describe the molecular characterization of a balanced familial chromosomal translocation, t(16;22)(p13.3;q11.2), segregating with congenital lamellar cataract. This led to the discovery of a cluster of lens-derived expressed sequence tags (ESTs) close to the 16p13.3 breakpoint. This region harbors a locus associated with cataract and microphthalmia. Long-range PCR and 16p13.3 breakpoint sequencing identified genomic sequence in a human genome sequence gap, and allowed identification of a novel four-exon gene, designated TMEM114, which encodes a predicted protein of 223 amino acids. The breakpoint lies in the promoter region of TMEM114 and separates the gene from predicted eye-specific upstream transcription factor binding sites. There is sequence conservation among orthologs down to zebrafish. The protein is predicted to contain four transmembrane domains with homology to the lens intrinsic membrane protein, LIM2 (also known as MP20), in the PMP-22/EMP/MP20 family. TMEM114 mutation screening in 130 congenital cataract patients revealed missense mutations leading to the exchange of highly-conserved amino acids in the first extracellular domain of the protein (p.I35T, p.F106L) in two separate patients and their reportedly healthy sibling and mother, respectively. In the lens, Tmem114 shows expression in the lens epithelial cells extending into the transitional zone where early fiber differentiation occurs. Our findings implicate dysregulation of expression of this novel human gene, TMEM114, in mammalian cataract formation.

Constraints on utilization of the EDA-signaling pathway in threespine stickleback evolution

Many traits evolve in parallel in widely separated populations. The evolutionary radiation of threespine sticklebacks provides a powerful model for testing the molecular basis of parallel evolution in vertebrates. Although marine sticklebacks are completely covered with bony armor plates, most freshwater populations have dramatic reductions in plates. Recent genetic studies have shown that major changes in armor patterning are likely due to regulatory alterations in the gene encoding the secreted signaling molecule ectodysplasin (EDA). In mammals, mutations in many different components of the EDA-signaling pathway produce similar changes in hair, teeth, sweat glands, and dermal bones. To test whether other genes in the EDA pathway also control natural variation in armor plates, we identified and mapped stickleback EDA Receptor (EDAR), the EDAR-Associated Death Domain adaptor, Tumor Necrosis Factor Receptor (TNFR) SuperFamily member 19, its adaptor TNFR-Associated Factor 6, and the downstream regulator nuclear factor kappa B Essential Modulator (NEMO). In contrast to the diversity of genes underlying ectodermal dysplasia disease phenotypes in humans, none of these EDA pathway components map to chromosomes previously shown to modify armor plates in natural populations, though EDAR showed a small but significant effect on plate number. We further investigated whether these genes exhibit differences in copy number, target size, or genomic organization that might make them less suitable targets for evolutionary change. In comparison with EDA, all these genes have smaller surrounding noncoding (putative regulatory) regions, with fewer evolutionarily conserved regions. We suggest that the presence of highly modular cis-acting control sequences may be a key factor influencing the likelihood that particular genes will serve as the basis of major phenotypic changes in nature.

Relevance of BAC transgene copy number in mice: transgene copy number variation across multiple transgenic lines and correlations with transgene integrity and expression

Bacterial artificial chromosomes (BACs) are excellent tools for manipulating large DNA fragments and, as a result, are increasingly utilized to engineer transgenic mice by pronuclear injection. The demand for BAC transgenic mice underscores the need for careful inspection of BAC integrity and fidelity following transgenesis, which may be crucial for interpreting transgene function. Thus, it is imperative that reliable methods for assessing these parameters are available. However, there are limited data regarding whether BAC transgenes routinely integrate in the mouse genome as intact molecules, how BAC transgenes behave as they are passed through the germline across successive generations, and how variation in BAC transgene copy number relates to transgene expression. To address these questions, we used TaqMan real-time PCR to estimate BAC transgene copy number in BAC transgenic embryos and lines. Here we demonstrate the reproducibility of copy number quantification with this method and describe the variation in copy number across independent transgenic lines. In addition, polymorphic marker analysis suggests that the majority of BAC transgenic lines contain intact molecules. Notably, all lines containing multiple BAC copies also contain all BAC-specific markers. Three of 23 founders analyzed contained BAC transgenes integrated into more than one genomic location. Finally, we show increased BAC transgene copy number correlates with increased BAC transgene expression. In sum, our efforts have provided a reliable method for assaying BAC transgene integrity and fidelity, and data that should be useful for researchers using BACs as transgenic vectors.