Comparative genomics of the SOX9 region in human and Fugu rubripes: conservation of short regulatory sequence elements within large intergenic regions

Conserved Noncoding Elements Evolve Around the Same Genes Throughout Metazoan Evolution

Conserved noncoding elements (CNEs) are DNA sequences located outside of protein-coding genes that can remain under purifying selection for up to hundreds of millions of years. Studies in vertebrate genomes have revealed that most CNEs carry out regulatory functions. Notably, many of them are enhancers that control the expression of homeodomain transcription factors and other genes that play crucial roles in embryonic development. To further our knowledge of CNEs in other parts of the animal tree, we conducted a large-scale characterization of CNEs in more than 50 genomes from three of the main branches of the metazoan tree: Cnidaria, Mollusca, and Arthropoda. We identified hundreds of thousands of CNEs and reconstructed the temporal dynamics of their appearance in each lineage, as well as determining their spatial distribution across genomes. We show that CNEs evolve repeatedly around the same genes across the Metazoa, including around homeodomain genes and other transcription factors; they also evolve repeatedly around genes involved in neural development. We also show that transposons are a major source of CNEs, confirming previous observations from vertebrates and suggesting that they have played a major role in wiring developmental gene regulatory mechanisms since the dawn of animal evolution.

Impact of Zinc Transport Mechanisms on Embryonic and Brain Development

The trace element zinc (Zn) binds to over ten percent of proteins in eukaryotic cells. Zn flexible chemistry allows it to regulate the activity of hundreds of enzymes and influence scores of metabolic processes in cells throughout the body. Deficiency of Zn in humans has a profound effect on development and in adults later in life, particularly in the brain, where Zn deficiency is linked to several neurological disorders. In this review, we will summarize the importance of Zn during development through a description of the outcomes of both genetic and early dietary Zn deficiency, focusing on the pathological consequences on the whole body and brain. The epidemiology and the symptomology of Zn deficiency in humans will be described, including the most studied inherited Zn deficiency disease, Acrodermatitis enteropathica. In addition, we will give an overview of the different forms and animal models of Zn deficiency, as well as the 24 Zn transporters, distributed into two families: the ZIPs and the ZnTs, which control the balance of Zn throughout the body. Lastly, we will describe the TRPM7 ion channel, which was recently shown to contribute to intestinal Zn absorption and has its own significant impact on early embryonic development.

Evolutionary conservation in noncoding genomic regions

Humans may share more genomic commonalities with other species than previously thought. According to current estimates, ~5% of the human genome is functionally constrained, which is a much larger fraction than the ~1.5% occupied by annotated protein-coding genes. Hence, ~3.5% of the human genome comprises likely functional conserved noncoding elements (CNEs) preserved among organisms, whose common ancestors existed throughout hundreds of millions of years of evolution. As whole-genome sequencing emerges as a standard procedure in genetic analyses, interpretation of variations in CNEs, including the elucidation of mechanistic and functional roles, becomes a necessity. Here, we discuss the phenomenon of noncoding conservation via four dimensions (sequence, regulatory conservation, spatiotemporal expression, and structure) and the potential significance of CNEs in phenotype variation and disease.

Sequence alteration in the enhancer contributes to the heterochronic Sox9 expression in marsupial cranial neural crest

Neonates of marsupial mammals are altricial at birth, because their gestation period is relatively short compared to placental mammals. Yet, as they need to travel to the teat from the birth canal, and suckle on the mother's milk, forelimbs and jaws develop significantly early. Previous studies in opossum (Monodelphis domestica), an experimental marsupial model, have revealed that cranial neural crest cells are generated significantly early compared to those in placental mammals, such as mouse, leading to an early development of jaw primordia. We have previously found that Sox9, an important neural crest-specifier gene, is expressed in the future cranial neural crest of the opossum embryonic ectoderm significantly earlier than that in mouse or quail embryos. As Sox9 is essential for neural crest formation in various vertebrates, it seems likely that the heterochronic expression of Sox9 is critical for the early cranial neural crest formation in the marsupial embryos. In this study, we show a marsupial-specific sequence in the Sox9 neural crest enhancer E3. We also reveal that the mouse E3 enhancer is activated in the cranial neural crest cells of quail embryos, that the E3 enhancer with marsupial-specific sequence is activated earlier in the Pax7-expressing neural border prior to the onset of endogenous Sox9 expression, and that a misexpression of cMyb, which is also a transcriptional activator of Pax7, in the neural border can ectopically activate the "marsupialized" enhancer. Thus, we suggest that the modification of the E3 enhancer sequence in the marsupial ancestor would have promoted the early expression of Sox9 in the neural border, facilitating the early formation of the cranial neural crest cells and the subsequent heterochronic development of the jaw primordia.

Non-coding variation in disorders of sex development

Genetic studies in Disorders of Sex Development (DSD), representing a wide spectrum of developmental or functional conditions of the gonad, have mainly been oriented towards the coding genome. Application of genomic technologies, such as whole-exome sequencing, result in a molecular genetic diagnosis in ∼50% of cases with DSD. Many of the genes mutated in DSD encode transcription factors such as SRY, SOX9, NR5A1, and FOXL2, characterized by a strictly regulated spatiotemporal expression. Hence, it can be hypothesized that at least part of the missing genetic variation in DSD can be explained by non-coding mutations in regulatory elements that alter gene expression, either by reduced, mis- or overexpression of their target genes. In addition, structural variations such as translocations, deletions, duplications or inversions can affect the normal chromatin conformation by different mechanisms. Here, we review non-coding defects in human DSD phenotypes and in animal models. The wide variety of non-coding defects found in DSD emphasizes that the regulatory landscape of known and to be discovered DSD genes has to be taken into consideration when investigating the molecular pathogenesis of DSD.

RNA binding proteins regulate anabolic and catabolic gene expression in chondrocytes

OBJECTIVE

Regulation of anabolic and catabolic factors is considered essential in maintaining the homoeostasis of healthy articular cartilage. In this study we investigated the influence of RNA binding proteins (RNABPs) in this process.

DESIGN

Using small interfering RNA (siRNA), RNABP expression was knocked down in SW1353 chondrosarcoma cells and human articular chondrocytes. Gene expression and messenger RNA (mRNA) decay of anabolic (SOX9, Aggrecan) and catabolic (matrix metalloproteinase (MMP)13) factors were analysed using reverse transcription quantitative polymerase chain reaction (RT-qPCR). RNA-electromobility shift assays (EMSAs) were used to investigate RNABP interactions with the SOX9 mRNA 3' untranslated region (UTR). Immunohistochemical localisation of MMP13 and the RNABP human antigen R (HuR) was performed in E13.5 and E16.5 mouse embryo sections.

RESULTS

SOX9 mRNA, mRNA half-life and protein expression were increased with siRNA targeting the RNABP tristetraprolin (TTP) in both HACs and SW1353s. TTP knockdown also stimulated aggrecan mRNA expression but did not affect its stability. RNA-EMSAs demonstrated that adenine uracil (AU)-rich elements in the SOX9 mRNA 3'UTR interacted with chondrocyte proteins with three specific elements interacting with TTP. HuR knockdown significantly increased MMP13 expression and also regulated the expression of a number of known transcriptional repressors of MMP13. HuR was ubiquitously expressed within mouse embryos yet displayed regional down-regulation within developing skeletal structures.

CONCLUSION

This study demonstrates for the first time how RNABPs are able to affect the balance of anabolic and catabolic gene expression in human chondrocytes. The post-transcriptional mechanisms controlled by RNABPs present novel avenues of regulation and potential points of intervention for controlling the expression of SOX9 and MMP13 in chondrocytes.

Identification of Sex-Linked SNPs and Sex-Determining Regions in the Yellowtail Genome

Unlike the conservation of sex-determining (SD) modes seen in most mammals and birds, teleost fishes exhibit a wide variety of SD systems and genes. Hence, the study of SD genes and sex chromosome turnover in fish is one of the most interesting topics in evolutionary biology. To increase resolution of the SD gene evolutionary trajectory in fish, identification of the SD gene in more fish species is necessary. In this study, we focused on the yellowtail, a species widely cultivated in Japan. It is a member of family Carangidae in which no heteromorphic sex chromosome has been observed, and no SD gene has been identified to date. By performing linkage analysis and BAC walking, we identified a genomic region and SNPs with complete linkage to yellowtail sex. Comparative genome analysis revealed the yellowtail SD region ancestral chromosome structure as medaka-fugu. Two inversions occurred in the yellowtail linage after it diverged from the yellowtail-medaka ancestor. An association study using wild yellowtails and the SNPs developed from BAC ends identified two SNPs that can reasonably distinguish the sexes. Therefore, these will be useful genetic markers for yellowtail breeding. Based on a comparative study, it was suggested that a PDZ domain containing the GIPC protein might be involved in yellowtail sex determination. The homomorphic sex chromosomes widely observed in the Carangidae suggest that this family could be a suitable marine fish model to investigate the early stages of sex chromosome evolution, for which our results provide a good starting point.

Otic placode cell specification and proliferation are regulated by Notch signaling in avian development

BACKGROUND

The entire inner ear including the cochlear-vestibular ganglion arises from a simple epithelium, the otic placode. Precursors for the placode originate from a pool of progenitors located in ectoderm next to the future hindbrain, the pre-otic field, where they are intermingled with future epibranchial and epidermal cells. While the importance of secreted proteins, such as FGFs and Wnts, in imparting otic identity has been well studied, how precursors for these different fates segregate locally is less well understood.

RESULTS

(1) The Notch ligand Delta1 and the Notch target Hes5-2 are expressed in a part of pre-otic field before otic commitment, indicative of active Notch signaling, and this is confirmed using a Notch reporter. (2) Loss and gain-of-function approaches reveal that Notch signaling regulates both proliferation and specification of pre-otic progenitors.

CONCLUSIONS

Our results identify a novel function of Notch signaling in cell fate determination in the pre-otic field of avian embryos.

Deep conservation of cis-regulatory elements in metazoans

Despite the vast morphological variation observed across phyla, animals share multiple basic developmental processes orchestrated by a common ancestral gene toolkit. These genes interact with each other building complex gene regulatory networks (GRNs), which are encoded in the genome by cis-regulatory elements (CREs) that serve as computational units of the network. Although GRN subcircuits involved in ancient developmental processes are expected to be at least partially conserved, identification of CREs that are conserved across phyla has remained elusive. Here, we review recent studies that revealed such deeply conserved CREs do exist, discuss the difficulties associated with their identification and describe new approaches that will facilitate this search.

Mapping association between long-range cis-regulatory regions and their target genes using synteny

In chordates, long-range cis-regulatory regions are involved in the control of transcription initiation (either as repressors or enhancers). Their main characteristics are that (i) they can be located as far as 1 Mb away from the transcription start site of the target gene, (ii) they can regulate more than one gene, and (iii) they are usually orientation-independent. Therefore, proper characterization of functional interactions between long-range cis-regulatory regions and their target genes remains problematic. We present a novel method to predict such interactions based on the analysis of rearrangements between the human and 16 other vertebrate genomes. Our method is based on the assumption that genome rearrangements that would disrupt the functional interaction between a cis-regulatory region and its target gene are likely to be deleterious. Therefore, conservation of synteny through evolution would be an indication of a functional interaction. We use our algorithm to predict the association between a set of 123,905 human candidate regulatory regions to their target gene(s). This genome-wide map of interactions has many potential applications, including the selection of candidate regions prior to in vivo experimental characterization, a better characterization of regulatory regions involved in position effect diseases, and an improved understanding of the mechanisms and importance of long-range regulation.

Copy number variation in patients with disorders of sex development due to 46,XY gonadal dysgenesis

Disorders of sex development (DSD), ranging in severity from mild genital abnormalities to complete sex reversal, represent a major concern for patients and their families. DSD are often due to disruption of the genetic programs that regulate gonad development. Although some genes have been identified in these developmental pathways, the causative mutations have not been identified in more than 50% 46,XY DSD cases. We used the Affymetrix Genome-Wide Human SNP Array 6.0 to analyse copy number variation in 23 individuals with unexplained 46,XY DSD due to gonadal dysgenesis (GD). Here we describe three discrete changes in copy number that are the likely cause of the GD. Firstly, we identified a large duplication on the X chromosome that included DAX1 (NR0B1). Secondly, we identified a rearrangement that appears to affect a novel gonad-specific regulatory region in a known testis gene, SOX9. Surprisingly this patient lacked any signs of campomelic dysplasia, suggesting that the deletion affected expression of SOX9 only in the gonad. Functional analysis of potential SRY binding sites within this deleted region identified five putative enhancers, suggesting that sequences additional to the known SRY-binding TES enhancer influence human testis-specific SOX9 expression. Thirdly, we identified a small deletion immediately downstream of GATA4, supporting a role for GATA4 in gonad development in humans. These CNV analyses give new insights into the pathways involved in human gonad development and dysfunction, and suggest that rearrangements of non-coding sequences disturbing gene regulation may account for significant proportion of DSD cases.

Disruption of long-distance highly conserved noncoding elements in neurocristopathies

One of the key discoveries of vertebrate genome sequencing projects has been the identification of highly conserved noncoding elements (CNEs). Some characteristics of CNEs include their high frequency in mammalian genomes, their potential regulatory role in gene expression, and their enrichment in gene deserts nearby master developmental genes. The abnormal development of neural crest cells (NCCs) leads to a broad spectrum of congenital malformation(s), termed neurocristopathies, and/or tumor predisposition. Here we review recent findings that disruptions of CNEs, within or at long distance from the coding sequences of key genes involved in NCC development, result in neurocristopathies via the alteration of tissue- or stage-specific long-distance regulation of gene expression. While most studies on human genetic disorders have focused on protein-coding sequences, these examples suggest that investigation of genomic alterations of CNEs will provide a broader understanding of the molecular etiology of both rare and common human congenital malformations.

Conserved noncoding elements and the evolution of animal body plans

The genomes of vertebrates, flies, and nematodes contain highly conserved noncoding elements (CNEs). CNEs cluster around genes that regulate development, and where tested, they can act as transcriptional enhancers. Within an animal group CNEs are the most conserved sequences but between groups they are normally diverged beyond recognition. Alternative CNEs are, however, associated with an overlapping set of genes that control development in all animals. Here, we discuss the evidence that CNEs are part of the core gene regulatory networks (GRNs) that specify alternative animal body plans. The major animal groups arose >550 million years ago. We propose that the cis-regulatory inputs identified by CNEs arose during the "re-wiring" of regulatory interactions that occurred during early animal evolution. Consequently, different animal groups, with different core GRNs, contain alternative sets of CNEs. Due to the subsequent stability of animal body plans, these core regulatory sequences have been evolving in parallel under strong purifying selection in different animal groups.

Highly conserved non-coding elements on either side of SOX9 associated with Pierre Robin sequence

Pierre Robin sequence (PRS) is an important subgroup of cleft palate. We report several lines of evidence for the existence of a 17q24 locus underlying PRS, including linkage analysis results, a clustering of translocation breakpoints 1.06-1.23 Mb upstream of SOX9, and microdeletions both approximately 1.5 Mb centromeric and approximately 1.5 Mb telomeric of SOX9. We have also identified a heterozygous point mutation in an evolutionarily conserved region of DNA with in vitro and in vivo features of a developmental enhancer. This enhancer is centromeric to the breakpoint cluster and maps within one of the microdeletion regions. The mutation abrogates the in vitro enhancer function and alters binding of the transcription factor MSX1 as compared to the wild-type sequence. In the developing mouse mandible, the 3-Mb region bounded by the microdeletions shows a regionally specific chromatin decompaction in cells expressing Sox9. Some cases of PRS may thus result from developmental misexpression of SOX9 due to disruption of very-long-range cis-regulatory elements.

The words of the regulatory code are arranged in a variable manner in highly conserved enhancers

The cis-regulatory regions of many developmental regulators and transcription factors are believed to be highly conserved in the genomes of vertebrate species, suggesting specific regulatory mechanisms for these gene classes. We functionally characterized five notochord enhancers, whose sequence is highly conserved, and systematically mutated two of them. Two subregions were identified to be essential for expression in the notochord of the zebrafish embryo. Synthetic enhancers containing the two essential regions in front of a TATA-box drive expression in the notochord while concatemerization of the subregions alone is not sufficient, indicating that the combination of the two sequence elements is required for notochord expression. Both regions are present in the five functionally characterized notochord enhancers. However, the position, the distance and relative orientation of the two sequence motifs can vary substantially within the enhancer sequences. This suggests that the regulatory grammar itself does not dictate the high evolutionary conservation between these orthologous cis-regulatory sequences. Rather, it represents a less well-conserved layer of sequence organization within these sequences.

Ectopic SOX9 mediates extracellular matrix deposition characteristic of organ fibrosis

Appropriate temporospatial expression of the transcription factor SOX9 is important for normal development of a wide range of organs. Here, we show that when SOX9 is expressed ectopically, target genes become expressed that are associated with disease. Histone deacetylase inhibitors in clinical trials for cancer therapy induced SOX9 expression via enhanced recruitment of nuclear factor Y (NF-Y) to CCAAT elements in the SOX9 proximal promoter. The effect of histone deacetylase inhibitors could be elicited in cells that normally lack SOX9, such as hepatocytes. In human fetal hepatocytes, this aberrant induction of SOX9 protein caused ectopic expression of COL2A1 and COMP1 that encode extracellular matrix (ECM) components normally associated with chondrogenesis. Previously, ectopic expression of this "chondrogenic" profile has been implicated in vascular calcification. More broadly, inappropriate ECM deposition is a hallmark of fibrosis. We demonstrated that induction of SOX9 expression also occurred during activation of fibrogenic cells from the adult liver when the transcription factor was responsible for expression of the major component of fibrotic ECM, type 1 collagen. These combined data identify new aspects in the regulation of SOX9 expression. They support a role for SOX9 beyond normal development as a transcriptional regulator in the pathology of fibrosis.

Comparative genomics for detecting human disease genes

Originally, comparative genomics was geared toward defining the synteny of genes between species. As the human genome project accelerated, there was an increase in the number of tools and means to make comparisons culminating in having the genomic sequence for a large number of organisms spanning the evolutionary tree. With this level of resolution and a long history of comparative biology and comparative genetics, it is now possible to use comparative genomics to build or select better animal models and to facilitate gene discovery. Comparative genomics takes advantage of the functional genetic information from other organisms, (vertebrates and invertebrates), to apply it to the study of human physiology and disease. It allows for the identification of genes and regulatory regions, and for acquiring knowledge about gene function. In this chapter, the current state of comparative genomics and the available tools are discussed in the context of developing animal model systems that reflect the clinical picture.

Organization of conserved elements near key developmental regulators in vertebrate genomes

Sequence conservation has traditionally been used as a means to target functional regions of complex genomes. In addition to its use in identifying coding regions of genes, the recent availability of whole genome data for a number of vertebrates has permitted high-resolution analyses of the noncoding "dark matter" of the genome. This has resulted in the identification of a large number of highly conserved sequence elements that appear to be preserved in all bony vertebrates. Further positional analysis of these conserved noncoding elements (CNEs) in the genome demonstrates that they cluster around genes involved in developmental regulation. This chapter describes the identification and characterization of these elements, with particular reference to their composition and organization.

Long-range gene control and genetic disease

The past two decades have seen great progress in the elucidation of the genetic basis of human genetic disease. Many clinical phenotypes have been linked with mutations or deletions in specific causative genes. However, it is often less recognized that in addition to the integrity of the protein-coding sequences, human health critically also depends on the spatially, temporally, and quantitatively correct expression of those genes. Genetic disease can therefore equally be caused by disruption of the regulatory mechanisms that ensure proper gene expression. The term "position effect" is used in those situations where the expression level of a gene is deleteriously affected by an alteration in its chromosomal environment, while maintaining an intact transcription unit. Here, we review recent advances in our understanding of the possible mechanisms of a number of "position effect" disease cases and discuss the findings with respect to current models for genome organization and long-range control of gene expression.

Retroviral enhancer detection insertions in zebrafish combined with comparative genomics reveal genomic regulatory blocks - a fundamental feature of vertebrate genomes

A large-scale enhancer detection screen was performed in the zebrafish using a retroviral vector carrying a basal promoter and a fluorescent protein reporter cassette. Analysis of insertional hotspots uncovered areas around developmental regulatory genes in which an insertion results in the same global expression pattern, irrespective of exact position. These areas coincide with vertebrate chromosomal segments containing identical gene order; a phenomenon known as conserved synteny and thought to be a vestige of evolution. Genomic comparative studies have found large numbers of highly conserved noncoding elements (HCNEs) spanning these and other loci. HCNEs are thought to act as transcriptional enhancers based on the finding that many of those that have been tested direct tissue specific expression in transient or transgenic assays. Although gene order in hox and other gene clusters has long been known to be conserved because of shared regulatory sequences or overlapping transcriptional units, the chromosomal areas found through insertional hotspots contain only one or a few developmental regulatory genes as well as phylogenetically unrelated genes. We have termed these regions genomic regulatory blocks (GRBs), and show that they underlie the phenomenon of conserved synteny through all sequenced vertebrate genomes. After teleost whole genome duplication, a subset of GRBs were retained in two copies, underwent degenerative changes compared with tetrapod loci that exist as single copy, and that therefore can be viewed as representing the ancestral form. We discuss these findings in light of evolution of vertebrate chromosomal architecture and the identification of human disease mutations.

Hypoxia-inducible factor 2alpha is essential for hypoxic induction of the human articular chondrocyte phenotype

OBJECTIVE

To uncover the mechanism by which hypoxia enhances cartilage matrix synthesis by human articular chondrocytes.

METHODS

The hypoxic response was investigated by exposing normal (nonarthritic) human articular chondrocyte cultures to 20% oxygen and 1% oxygen. Induction of the differentiated phenotype was confirmed at the gene and protein levels. In its first reported application in human articular chondrocytes, the RNA interference method was used to directly investigate the role of specific transcription factors in this process. Small interfering RNA directed against hypoxia-inducible factor 1alpha (HIF-1alpha), HIF-2alpha, and SOX9 were delivered by lipid-based transfection of primary and passaged human articular chondrocytes. The effect of each knockdown on hypoxic induction of the chondrocyte phenotype was assessed.

RESULTS

Hypoxia enhanced matrix synthesis and SOX9 expression of human articular chondrocytes at both the gene and protein levels. Although HIF-1alpha knockdown had no effect, depletion of HIF-2alpha abolished this hypoxic induction. Thus, we provide the first evidence that HIF-2alpha, but not HIF-1alpha, is essential for hypoxic induction of the human articular chondrocyte phenotype. In addition, depletion of SOX9 prevented hypoxic induction of matrix genes, indicating that the latter are not direct HIF targets but are up-regulated by hypoxia via SOX9.

CONCLUSION

Based on our data, we propose a novel mechanism whereby hypoxia promotes cartilage matrix synthesis specifically through HIF-2alpha-mediated SOX9 induction of key cartilage genes. These findings have potential application for the development of cartilage repair therapies.

Multiple conserved regulatory elements with overlapping functions determine Sox10 expression in mouse embryogenesis

Expression and function of the transcription factor Sox10 is predominant in neural crest cells, its derivatives and in oligodendrocytes. To understand how Sox10 expression is regulated during development, we analysed the potential of evolutionary conserved non-coding sequences in the Sox10 genomic region to function as enhancers. By linking these sequences to a β-galactosidase marker gene under the control of a minimal promoter, five regulatory regions were identified that direct marker gene expression in transgenic mice to Sox10 expressing cell types and tissues in a defined temporal pattern. These possible enhancers of the Sox10 gene mediate Sox10 expression in the otic vesicle, in oligodendrocytes and in several neural crest derivatives including the developing peripheral nervous system and the adrenal gland. They furthermore exhibit overlapping activities and share binding sites for Sox, Lef/Tcf, Pax and AP2 transcription factors. This may explain high level and robustness of Sox10 expression during embryonic development.

Correct dosage of Fog2 and Gata4 transcription factors is critical for fetal testis development in mice

Previous reports suggested that humans and mice differ in their sensitivity to the genetic dosage of transcription factors that play a role in early testicular development. This difference implies that testis determination might be somewhat different in these two species. We report that the Fog2 and Gata4 transcription factors are haploinsufficient for testis determination in mice. Whether gonadal sex reversal occurs depends on genetic background (i.e., modifier genes). For example, C57BL/6J (B6) XY mice develop testes if they are heterozygous for a mutant Fog2 (Fog2-) or Gata4 (Gata4(ki)) allele. However, if the B6 Y chromosome (Y(B6)) is replaced by the AKR Y chromosome (Y(AKR)), B6 Fog2-/+ XY(AKR) mice develop ovaries, and B6 Gata4(ki)/+ XY(AKR) mice develop ovaries and ovotestes (gonads containing both ovarian and testicular tissue). Furthermore, DBA/2J (D2) Fog2-/+ XY(AKR) mice and (B6 x D2)F1 hybrid Gata4(ki)/+ XY(AKR) mice develop testes. Sry is expressed in the mutant XY gonads, indicating that the lack of Sry expression is not the cause of ovarian tissue development in B6 Fog2-/+ or Gata4(ki)/+ XY(AKR) mice. However, up-regulation of Sox9 expression, which is critical for normal testicular development, does not occur in mutant XY gonads that develop as ovaries. We conclude that under certain genetic conditions, Sox9 up-regulation depends on the proper dosage of Fog2 and Gata4. We propose that in humans the FOG2 and/or GATA4 genes might be haploinsufficient for normal testis determination and thus could be the cause of some previously unassigned cases of XY gonadal sex reversal.

Comparative genomics reveals functional transcriptional control sequences in the Prop1 gene

Mutations in PROP1 are a common genetic cause of multiple pituitary hormone deficiency (MPHD). We used a comparative genomics approach to predict the transcriptional regulatory domains of Prop1 and tested them in cell culture and mice. A BAC transgene containing Prop1 completely rescues the Prop1 mutant phenotype, demonstrating that the regulatory elements necessary for proper PROP1 transcription are contained within the BAC. We generated DNA sequences from the PROP1 genes in lemur, pig, and five different primate species. Comparison of these with available human and mouse PROP1 sequences identified three putative regulatory sequences that are highly conserved. These are located in the PROP1 promoter proximal region, within the first intron of PROP1, and downstream of PROP1. Each of the conserved elements elicited orientation-specific enhancer activity in the context of the Drosophila alcohol dehydrogenase minimal promoter in both heterologous and pituitary-derived cells lines. The intronic element is sufficient to confer dorsal expansion of the pituitary expression domain of a transgene, suggesting that this element is important for the normal spatial expression of endogenous Prop1 during pituitary development. This study illustrates the usefulness of a comparative genomics approach in the identification of regulatory elements that may be the site of mutations responsible for some cases of MPHD.

SOX9cre1, a cis-acting regulatory element located 1.1 Mb upstream of SOX9, mediates its enhancement through the SHH pathway

SOX9 is a temporal and tissue-specific transcription factor involved in male sexual development and bone formation. Haploinsufficiency of SOX9 is known to cause campomelic dysplasia (CD). CD cases without SOX9 coding region mutations have been described in association with translocations that have breakpoints mapping as far as 932 kb upstream from the gene. These rearrangements suggest that position effects acting from a great distance regulate SOX9 gene expression. Studies of one such case (900 kb upstream to SOX9) have led to the delineation of a potential 2.1 kb cis-acting regulatory element 1.1 Mb upstream of SOX9, termed SOX9cre1. We investigated the role of this putative regulator in SOX9 expression. SOX9cre1 increases the activity of a minimal SOX9 promoter in reporter constructs in a dose-dependent and tissue-specific manner, consistent with an enhancer role. In silico studies identify a putative binding site within SOX9cre1 for GLI1, a downstream mediator of sonic hedgehog (SHH). Furthermore, the stimulation of primary human chondrocyte cells in culture with SHH increases endogenous SOX9 expression 3-fold. Electrophoresis mobility shift assay (EMSA) studies that demonstrate physical interactions between the GLI1 transcription factor and a putative binding site within SOX9cre1, as well as experiments in which reporter constructs are co-transfected with GLI1, suggest a direct interaction between GLI1 and SOX9cre1. GLI1-SOX9cre1 interactions are verified in chromatin immunoprecipitation experiments. These data support a direct molecular link between the Hh signaling pathway and SOX9 regulation, wherein SHH stimulates SOX9 through its mediator GLI1, and are consistent with a mechanism of SOX9 regulation through distal chromatin interactions.

Two novel translocation breakpoints upstream of SOX9 define borders of the proximal and distal breakpoint cluster region in campomelic dysplasia

The semilethal skeletal malformation syndrome campomelic dysplasia (CD) with or without XY sex reversal is caused by mutations within the SOX9 gene on 17q24.3 or by chromosomal aberrations (translocations, inversions or deletions) with breakpoints outside the SOX9 coding region. The previously published CD translocation breakpoints upstream of SOX9 fall into two clusters: a proximal cluster with breakpoints between 50-300 kb and a distal cluster with breakpoints between 899-932 kb. Here, we present clinical, cytogenetic and molecular data from two novel CD translocation cases. Case 1 with karyotype 46,XY,t(1;17)(q42.1;q24.3) has characteristic symptoms of CD, including mild tibial bowing, cryptorchidism and hypospadias. By standard fluorescence in situ hybridization (FISH) and by high-resolution fiber FISH, the 17q breakpoint was mapped 375 kb from SOX9, defining the centromeric border of the proximal breakpoint cluster region. Case 2 with karyotype 46,X,t(Y;17)(q11.2;q24.3) has the acampomelic form of CD and complete XY sex reversal. By FISH and somatic cell hybrid analysis, the 17q breakpoint was mapped 789 kb from SOX9, defining the telomeric border of the distal breakpoint cluster region. We discuss the structure of the 1 Mb cis-control region upstream of SOX9 and the correlation between the position of the 14 mapped translocation breakpoints with respect to disease severity and XY sex reversal.

Sex determination and gonadal development in mammals

Arguably the most defining moment in our lives is fertilization, the point at which we inherit either an X or a Y chromosome from our father. The profoundly different journeys of male and female life are thus decided by a genetic coin toss. These differences begin to unfold during fetal development, when the Y-chromosomal Sry ("sex-determining region Y") gene is activated in males and acts as a switch that diverts the fate of the undifferentiated gonadal primordia, the genital ridges, towards testis development. This sex-determining event sets in train a cascade of morphological changes, gene regulation, and molecular interactions that directs the differentiation of male characteristics. If this does not occur, alternative molecular cascades and cellular events drive the genital ridges toward ovary development. Once testis or ovary differentiation has occurred, our sexual fate is further sealed through the action of sex-specific gonadal hormones. We review here the molecular and cellular events (differentiation, migration, proliferation, and communication) that distinguish testis and ovary during fetal development, and the changes in gene regulation that underpin these two alternate pathways. The growing body of knowledge relating to testis development, and the beginnings of a picture of ovary development, together illustrate the complex mechanisms by which these organ systems develop, inform the etiology, diagnosis, and management of disorders of sexual development, and help define what it is to be male or female.

Comprehensive characterization of the cis-regulatory code responsible for the spatio-temporal expression of olSix3.2 in the developing medaka forebrain

BACKGROUND

Embryonic development is coordinated by sets of cis-regulatory elements that are collectively responsible for the precise spatio-temporal organization of regulatory gene networks. There is little information on how these elements, which are often associated with highly conserved noncoding sequences, are combined to generate precise gene expression patterns in vertebrates. To address this issue, we have focused on Six3, an important regulator of vertebrate forebrain development.

RESULTS

Using computational analysis and exploiting the diversity of teleost genomes, we identified a cluster of highly conserved noncoding sequences surrounding the Six3 gene. Transgenesis in medaka fish demonstrates that these sequences have enhancer, silencer, and silencer blocker activities that are differentially combined to control the entire distribution of Six3.

CONCLUSION

This report provides the first example of the precise regulatory code necessary for the expression of a vertebrate gene, and offers a unique framework for defining the interplay of trans-acting factors that control the evolutionary conserved use of Six3.

Distant regulatory elements in a Sox10-beta GEO BAC transgene are required for expression of Sox10 in the enteric nervous system and other neural crest-derived tissues

Sox10 is an essential transcription factor required for development of neural crest-derived melanocytes, peripheral glia, and enteric ganglia. Multiple transcriptional targets regulated by Sox10 have been identified; however, little is known regarding regulation of Sox10. High sequence conservation surrounding 5' exons 1 through 3 suggests these regions might contain functional regulatory elements. However, we observed that these Sox10 genomic sequences do not confer appropriate cell-specific transcription in vitro when linked to a heterologous reporter. To identify elements required for expression of Sox10 in vivo, we modified bacterial artificial chromosomes (BACs) to generate a Sox10betaGeoBAC transgene. Our approach leaves endogenous Sox10 loci unaltered, circumventing haploinsufficiency issues that arise from gene targeting. Sox10betaGeoBAC expression closely approximates Sox10 expression in vivo, resulting in expression in anterior dorsal neural tube at embryonic day (E) 8.5 and in cranial ganglia, otic vesicle, and developing dorsal root ganglia at E10.5. Characterization of Sox10betaGeoBAC expression confirms the presence of essential regulatory regions and additionally identifies previously unreported expression in thyroid parafollicular cells, thymus, salivary, adrenal, and lacrimal glands. Fortuitous deletions in independent Sox10betaGeoBAC lines result in loss of transgene expression in peripheral nervous system lineages and coincide with evolutionarily conserved regions. Our analysis indicates that Sox10 expression requires the presence of distant cis-acting regulatory elements. The Sox10betaGeoBAC transgene offers one avenue for specifically testing the role of individual conserved regions in regulation of Sox10 and makes possible analysis of Sox10+ derivatives in the context of normal neural crest development.

Misexpression of Sox9 in mouse limb bud mesenchyme induces polydactyly and rescues hypodactyly mice

Our previous studies have demonstrated the essential roles of the transcription factor Sox9 in the commitment of mesenchymal cells to a chondrogenic cell lineage and in overt chondrogenesis during limb bud development. However, it remains unknown if Sox9 induces chondrogenesis in mesenchyme ectopically in vivo as a master regulator of chondrogenesis. In this study, we first generated mutant mice in which Sox9 was misexpressed in the limb bud mesenchyme. The mutant mouse embryos exhibited polydactyly in limb buds in association with ectopic expression of Sox5 and Sox6 although markers for the different axes of limb bud development showed a normal pattern of expression. Misexpression of Sox9 stimulated cell proliferation in limb bud mesenchyme, suggesting that Sox9 has a role in recruiting mesenchymal cells to mesenchymal condensation. Second, despite the facts that misexpression of Sonic hedgehog (Shh) induces polydactyly in a number of mutant mice and Shh-null mutants have severely defective cartilage elements in limb buds, misexpression of Sox9 did not restore limb bud phenotypes in Shh-null mutants. Rather, there was no expression of Sox9 in digit I of Hoxa13Hd mutant embryos, and Sox9 partially rescued hypodactyly in Hoxa13Hd mutant embryos. These results provide evidence that Sox9 induces ectopic chondrogenesis in mesenchymal cells and strongly suggest that its expression may be regulated by Hox genes during limb bud development.

Regulation of SOX9 mRNA in human articular chondrocytes involving p38 MAPK activation and mRNA stabilization

Human articular chondrocytes rapidly lose their phenotype in monolayer culture. Recently we have shown that overexpression of the transcription factor SOX9 greatly enhanced re-expression of the phenotype in three-dimensional aggregate cultures. Here we show that endogenous SOX9 mRNA can be rapidly up-regulated in subcultured human articular chondrocytes if grown in alginate, in monolayer with cytochalasin D, or with specific inhibition of the RhoA effector kinases ROCK1 and -2, which all prevent actin stress fiber formation. Disruption of actin stress fibers using any of these redifferentiation stimuli also supported the superinduction of SOX9 by cycloheximide. The superinduction was blocked by inhibitors of the p38 MAPK signaling pathway and involved the stabilization of SOX9 mRNA. Furthermore stimulation of chondrocyte p38 MAPK activity with interleukin-1beta resulted in increased levels of SOX9 mRNA, and this was again dependent on the absence of actin stress fibers in the cells. In this study of chondrocyte redifferentiation we have provided further evidence of the early involvement of SOX9 and have discovered a novel post-transcriptional regulatory mechanism activated by p38 MAPK, which stabilized SOX9 mRNA.

The role of TGFbetas and Sox9 during limb chondrogenesis

The majority of the skeletal elements, except the flat bones of the skull, are formed by endochondral ossification, in which cartilage is replaced by bone. The formation of cartilage is a multi-step process termed chondrogenesis, during which undifferentiated mesenchymal cells condense and undergo differentiation towards chondrocytes. Notwithstanding recent advances, our knowledge of the detailed mechanisms implicated in cartilage and bone formation is still scarce. Recent genetic, cellular and biochemical studies have highlighted the importance of TGFbeta signaling and the activity of the transcription factor Sox9 during the early stages of vertebrate limb chondrogenesis.

Comparative evolutionary genomics of the HADH2 gene encoding Abeta-binding alcohol dehydrogenase/17beta-hydroxysteroid dehydrogenase type 10 (ABAD/HSD10)

Background

The Aβ-binding alcohol dehydrogenase/17β-hydroxysteroid dehydrogenase type 10 (ABAD/HSD10) is an enzyme involved in pivotal metabolic processes and in the mitochondrial dysfunction seen in the Alzheimer's disease. Here we use comparative genomic analyses to study the evolution of the HADH2 gene encoding ABAD/HSD10 across several eukaryotic species.

Results

Both vertebrate and nematode HADH2 genes showed a six-exon/five-intron organization while those of the insects had a reduced and varied number of exons (two to three). Eutherian mammal HADH2 genes revealed some highly conserved noncoding regions, which may indicate the presence of functional elements, namely in the upstream region about 1 kb of the transcription start site and in the first part of intron 1. These regions were also conserved between Tetraodon and Fugu fishes. We identified a conserved alternative splicing event between human and dog, which have a nine amino acid deletion, causing the removal of the strand β_F. This strand is one of the seven strands that compose the core β-sheet of the Rossman fold dinucleotide-binding motif characteristic of the short chain dehydrogenase/reductase (SDR) family members. However, the fact that the substrate binding cleft residues are retained and the existence of a shared variant between human and dog suggest that it might be functional. Molecular adaptation analyses across eutherian mammal orthologues revealed the existence of sites under positive selection, some of which being localized in the substrate-binding cleft and in the insertion 1 region on loop D (an important region for the Aβ-binding to the enzyme). Interestingly, a higher than expected number of nonsynonymous substitutions were observed between human/chimpanzee and orangutan, with six out of the seven amino acid replacements being under molecular adaptation (including three in loop D and one in the substrate binding loop).

Conclusion

Our study revealed that HADH2 genes maintained a reasonable conserved organization across a large evolutionary distance. The conserved noncoding regions identified among mammals and between pufferfishes, the evidence of an alternative splicing variant conserved between human and dog, and the detection of positive selection across eutherian mammals, may be of importance for further research on ABAD/HSD10 function and its implication in the Alzheimer's disease.

The Wilms tumor gene, Wt1, is required for Sox9 expression and maintenance of tubular architecture in the developing testis

Mutation of the transcription factor and tumor suppressor gene WT1 results in a range of genitourinary anomalies in humans, including 46,XY gonadal dysgenesis, indicating that WT1 plays a critical role in sex determination. However, because knockout of Wt1 in mice results in apoptosis of the genital ridge, it is unknown whether WT1 is required for testis development after the initial steps of sex determination. To address this question, we generated a mouse strain carrying a Wt1 conditional knockout allele and ablated Wt1 function specifically in Sertoli cells by embryonic day 14.5, several days after testis determination. Wt1 knockout resulted in disruption of developing seminiferous tubules and subsequent progressive loss of Sertoli cells and germ cells such that postnatal mutant testes were almost completely devoid of these cell types and were severely hypoplastic. Thus, Wt1 is essential for the maintenance of Sertoli cells and seminiferous tubules in the developing testes. Of particular note, expression of the testis-determining gene Sox9 in mutant Sertoli cells was turned off at embryonic day 14.5 after Wt1 ablation, suggesting that WT1 regulates Sox9, either directly or indirectly, after Sry expression ceases. Our data, along with previous work demonstrating the role of Wt1 at early stages of gonadal development, thus indicate that Wt1 is essential at multiple steps in testicular development.

Conservation of RET regulatory function from human to zebrafish without sequence similarity

Evolutionary sequence conservation is an accepted criterion to identify noncoding regulatory sequences. We have used a transposon-based transgenic assay in zebrafish to evaluate noncoding sequences at the zebrafish ret locus, conserved among teleosts, and at the human RET locus, conserved among mammals. Most teleost sequences directed ret-specific reporter gene expression, with many displaying overlapping regulatory control. The majority of human RET noncoding sequences also directed ret-specific expression in zebrafish. Thus, vast amounts of functional sequence information may exist that would not be detected by sequence similarity approaches.

Long-range upstream and downstream enhancers control distinct subsets of the complex spatiotemporal Sox9 expression pattern

SOX9 is an evolutionary conserved transcription factor that is expressed in a variety of tissues, with essential functions in cartilage, testis, heart, glial cell, inner ear and neural crest development. By comparing human and pufferfish genomic sequences, we previously identified eight highly conserved sequence elements between 290 kb 5' and 450 kb 3' to human SOX9. In this study, we assayed the regulatory potential of elements E1 to E7 in transgenic mice using a lacZ reporter gene driven by a 529 bp minimal mouse Sox9 promoter. We found that three of these elements and the Sox9 promoter control distinct subsets of the tissue-specific expression pattern of Sox9. E3, located 251 kb 5' to SOX9, directs lacZ expression to cranial neural crest cells and to the inner ear. E1 is located 28 kb 5' to SOX9 and controls expression in the node, notochord, gut, bronchial epithelium and pancreas. Transgene expression in the neuroectoderm is mediated by E7, located 95 kb 3' to SOX9, which regulates expression in the telencephalon and midbrain, and by the Sox9 minimal promoter which controls expression in the ventral spinal cord and hindbrain. We show that E3-directed reporter gene expression in neural crest cells of the first but not of the second and third pharyngeal arch is dependent on beta-catenin, revealing a complex regulation of Sox9 in cranial neural crest cells. Moreover, we identify and discuss highly conserved transcription factor binding sites within enhancer E3 that are in good agreement with current models for neural crest and inner ear development. Finally, we identify enhancer E1 as a cis-regulatory element conserved between vertebrates and invertebrates, indicating that some cis-regulatory sequences that control developmental genes in vertebrates might be phylogenetically ancient.

A novel approach to identifying regulatory motifs in distantly related genomes

A two-step procedure for identifying regulatory motifs in distantly related organisms is described that combines the advantages of sequence alignment and motif detection approaches.

Although proven successful in the identification of regulatory motifs, phylogenetic footprinting methods still show some shortcomings. To assess these difficulties, most apparent when applying phylogenetic footprinting to distantly related organisms, we developed a two-step procedure that combines the advantages of sequence alignment and motif detection approaches. The results on well-studied benchmark datasets indicate that the presented method outperforms other methods when the sequences become either too long or too heterogeneous in size.

Trade-offs in detecting evolutionarily constrained sequence by comparative genomics

As whole-genome sequencing efforts extend beyond more traditional model organisms to include a deep diversity of species, comparative genomic analyses will be further empowered to reveal insights into the human genome and its evolution. The discovery and annotation of functional genomic elements is a necessary step toward a detailed understanding of our biology, and sequence comparisons have proven to be an integral tool for that task. This review is structured to broadly reflect the statistical challenges in discriminating these functional elements from the bulk of the genome that has evolved neutrally. Specifically, we review the comparative genomics literature in terms of specificity, sensitivity, and phylogenetic scope, as well as the trade-offs that relate these factors in standard analyses. We consider the impact of an expanding diversity of orthologous sequences on our ability to resolve functional elements. This impact is assessed through both recent comparative analyses of deep alignments and mathematical modeling.

Deletion of long-range sequences at Sox10 compromises developmental expression in a mouse model of Waardenburg-Shah (WS4) syndrome

The transcription factor SOX10 is mutated in the human neurocristopathy Waardenburg-Shah syndrome (WS4), which is characterized by enteric aganglionosis and pigmentation defects. SOX10 directly regulates genes expressed in neural crest lineages, including the enteric ganglia and melanocytes. Although some SOX10 target genes have been reported, the mechanisms by which SOX10 expression is regulated remain elusive. Here, we describe a transgene-insertion mutant mouse line (Hry) that displays partial enteric aganglionosis, a loss of melanocytes, and decreased Sox10 expression in homozygous embryos. Mutation analysis of Sox10 coding sequences was negative, suggesting that non-coding regulatory sequences are disrupted. To isolate the Hry molecular defect, Sox10 genomic sequences were collected from multiple species, comparative sequence analysis was performed and software was designed (ExactPlus) to identify identical sequences shared among species. Mutation analysis of conserved sequences revealed a 15.9 kb deletion located 47.3 kb upstream of Sox10 in Hry mice. ExactPlus revealed three clusters of highly conserved sequences within the deletion, one of which shows strong enhancer potential in cultured melanocytes. These studies: (i) present a novel hypomorphic Sox10 mutation that results in a WS4-like phenotype in mice; (ii) demonstrate that a 15.9 kb deletion underlies the observed phenotype and likely removes sequences essential for Sox10 expression; (iii) combine a novel in silico method for comparative sequence analysis with in vitro functional assays to identify candidate regulatory sequences deleted in this strain. These studies will direct further analyses of Sox10 regulation and provide candidate sequences for mutation detection in WS4 patients lacking a SOX10-coding mutation.

Genetic background of HSH in three Polish families and a patient with an X;9 translocation

Hypomagnesemia with secondary hypocalcemia (HSH) is a rare inherited disease, characterised by neurological symptoms, such as tetany, muscle spasms and seizures, due to hypocalcemia. It has been suggested that HSH is genetically heterogeneous, but only one causative gene, TRPM6, on chromosome 9 has so far been isolated. We have now studied the genetic background of HSH in four Polish patients belonging to three families, and a HSH patient carrying an apparently balanced X;9 translocation. The translocation patient has long been considered as an example of the X-linked form of HSH. We identified six TRPM6 gene mutations, of which five were novel, in the Polish patients. All the alterations were either nonsense/splicing or missense mutations. The clinical picture of the patients was similar to the HSH patients reported earlier. No genotype-phenotype correlation could be detected. Sequencing did not reveal any TRPM6 or TRPM7 gene mutations in the female HSH patient with an X;9 translocation. Isolation of the translocation breakpoint showed that the chromosome 9 specific breakpoint mapped within satellite III repeat sequence. The X-chromosomal breakpoint was localised to the first intron of the vascular endothelial growth factor gene, VEGFD. No other sequence alterations were observed within the VEGFD gene. Even though the VEGFD gene was interrupted by the X;9 translocation, it seems unlikely that VEGFD is causing the translocation patient's HSH-like phenotype. Furthermore, re-evaluation of patient's clinical symptoms suggests that she did not have a typical HSH.

A functional survey of the enhancer activity of conserved non-coding sequences from vertebrate Iroquois cluster gene deserts

Recent studies of the genome architecture of vertebrates have uncovered two unforeseen aspects of its organization. First, large regions of the genome, called gene deserts, are devoid of protein-coding sequences and have no obvious biological role. Second, comparative genomics has highlighted the existence of an array of highly conserved non-coding regions (HCNRs) in all vertebrates. Most surprisingly, these structural features are strongly associated with genes that have essential functions during development. Among these, the vertebrate Iroquois (Irx) genes stand out on both fronts. Mammalian Irx genes are organized in two clusters (IrxA and IrxB) that span >1 Mb each with no other genes interspersed. Additionally, a large number of HCNRs exist within Irx clusters. We have systematically examined the enhancer activity of HCNRs from the IrxB cluster using transgenic Xenopus and zebrafish embryos. Most of these HCNRs are active in subdomains of endogenous Irx expression, and some are candidates to contain shared enhancers of neighboring genes, which could explain the evolutionary conservation of Irx clusters. Furthermore, HCNRs present in tetrapod IrxB but not in fish may be responsible for novel Irx expression domains that appeared after their divergence. Finally, we have performed a more detailed analysis on two IrxB ultraconserved non-coding regions (UCRs) duplicated in IrxA clusters in similar relative positions. These four regions share a core region highly conserved among all of them and drive expression in similar domains. However, inter-species conserved sequences surrounding the core, specific for each of these UCRs, are able to modulate their expression.

Secrets to a healthy Sox life: lessons for melanocytes

Sox proteins are transcriptional regulators with a high-mobility-group domain as sequence-specific DNA-binding domain. For function, they generally require other transcription factors as partner proteins. Sox proteins furthermore affect DNA topology and may shape the conformation of enhancer-bound multiprotein complexes as architectural proteins. Recent studies suggest that Sox proteins are tightly regulated in their expression by many signalling pathways, and that their transcriptional activity is subject to post-translational modification and sequestration mechanisms. Sox proteins are thus ideally suited to perform their many different functions as transcriptional regulators throughout mammalian development. Their unique properties also cause Sox proteins to escape detection in many standard transcription assays. In melanocytes, studies have so far focused on the Sox10 protein which functions both during melanocyte specification and at later times in the melanocyte lineage. During specification, Sox10 activates the Mitf gene as the key regulator of melanocyte development. At later stages, it ensures cell-type specific expression of melanocyte genes such as Dopachrome tautomerase. Both activities require cooperation with transcriptional partner proteins such as Pax-3, CREB and eventually Mitf. If predictions can be made from other cell lineages, further functions of Sox proteins in melanocytes may still lie ahead.

Identification of neuronal enhancers of the proopiomelanocortin gene by transgenic mouse analysis and phylogenetic footprinting

The proopiomelanocortin (POMC) gene is expressed in the pituitary and arcuate neurons of the hypothalamus. POMC arcuate neurons play a central role in the control of energy homeostasis, and rare loss-of-function mutations in POMC cause obesity. Moreover, POMC is the prime candidate gene within a highly significant quantitative trait locus on chromosome 2 associated with obesity traits in several human populations. Here, we identify two phylogenetically conserved neuronal POMC enhancers designated nPE1 (600 bp) and nPE2 (150 bp) located approximately 10 to 12 kb upstream of mammalian POMC transcriptional units. We show that mouse or human genomic regions containing these enhancers are able to direct reporter gene expression to POMC hypothalamic neurons, but not the pituitary of transgenic mice. Conversely, deletion of nPE1 and nPE2 in the context of the entire transcriptional unit of POMC abolishes transgene expression in the hypothalamus without affecting pituitary expression. Our results indicate that the nPEs are necessary and sufficient for hypothalamic POMC expression and that POMC expression in the brain and pituitary is controlled by independent sets of enhancers. Our study advances the understanding of the molecular nature of hypothalamic POMC neurons and will be useful to determine whether polymorphisms in POMC regulatory regions play a role in the predisposition to obesity.

Neuropeptides and the control of food intake in fish

The brain, particularly the hypothalamus, integrates input from factors that stimulate (orexigenic) and inhibit (anorexigenic) food intake. In fish, the identification of appetite regulators has been achieved by the use of both peptide injections followed by measurements of food intake, and by molecular cloning combined with gene expression studies. Neuropeptide Y (NPY) is the most potent orexigenic factor in fish. Other orexigenic peptides, orexin A and B and galanin, have been found to interact with NPY in the control of food intake in an interdependent and coordinated manner. On the other hand cholecystokinin (CCK), cocaine and amphetamine-regulated transcript (CART), and corticotropin-releasing factor (CRF) are potent anorexigenic factors in fish, the latter being involved in stress-related anorexia. CCK and CART have synergistic effects on food intake and modulate the actions of NPY and orexins. Although leptin has not yet been identified in fish, administration of mammalian leptin inhibits food intake in goldfish. Moreover, leptin induces CCK gene expression in the hypothalamus and its actions are mediated at least in part by CCK. Other orexigenic factors have been identified in teleost fish, including the agouti-related protein (AgRP) and ghrelin. Additional anorexigenic factors include bombesin (or gastrin-releasing peptide), alpha-melanocyte-stimulating hormone (alpha-MSH), tachykinins, and urotensin I. In goldfish, nutritional status can modify the expression of mRNAs encoding a number of these peptides, which provides further evidence for their roles as appetite regulators: (1) brain mRNA expression of CCK, CART, tachykinins, galanin, ghrelin, and NPY undergo peri-prandial variations; and (2) fasting increases the brain mRNA expression of NPY, AgRP, and ghrelin as well as serum ghrelin levels, and decreases the brain mRNA expression of tachykinins, CART, and CCK. This review will provide an overview of recent findings in this field.

Position effects due to chromosome breakpoints that map approximately 900 Kb upstream and approximately 1.3 Mb downstream of SOX9 in two patients with campomelic dysplasia

Campomelic dysplasia (CD) is a semilethal skeletal malformation syndrome with or without XY sex reversal. In addition to the multiple mutations found within the sex-determining region Y-related high-mobility group box gene (SOX9) on 17q24.3, several chromosome anomalies (translocations, inversions, and deletions) with breakpoints scattered over 1 Mb upstream of SOX9 have been described. Here, we present a balanced translocation, t(4;17)(q28.3;q24.3), segregating in a family with a mild acampomelic CD with Robin sequence. Both chromosome breakpoints have been identified by fluorescence in situ hybridization and have been sequenced using a somatic cell hybrid. The 17q24.3 breakpoint maps approximately 900 kb upstream of SOX9, which is within the same bacterial artificial chromosome clone as the breakpoints of two other reported patients with mild CD. We also report a prenatal identification of acampomelic CD with male-to-female sex reversal in a fetus with a de novo balanced complex karyotype, 46,XY,t(4;7;8;17)(4qter-->4p15.1::17q25.1-->17qter;7qter-->7p15.3::4p15.1-->4pter;8pter-->8q12.1::7p15.3-->7pter;17pter-->17q25.1::8q12.1-->8qter). Surprisingly, the 17q breakpoint maps approximately 1.3 Mb downstream of SOX9, making this the longest-range position effect found in the field of human genetics and the first report of a patient with CD with the chromosome breakpoint mapping 3' of SOX9. By using the Regulatory Potential score in conjunction with analysis of the rearrangement breakpoints, we identified a candidate upstream cis-regulatory element, SOX9cre1. We provide evidence that this 1.1-kb evolutionarily conserved element and the downstream breakpoint region colocalize with SOX9 in the interphase nucleus, despite being located 1.1 Mb upstream and 1.3 Mb downstream of it, respectively. The potential molecular mechanism responsible for the position effect is discussed.

Fine mapping of chromosome 17 translocation breakpoints > or = 900 Kb upstream of SOX9 in acampomelic campomelic dysplasia and a mild, familial skeletal dysplasia

Previously, our group reported a five-generation family in which a balanced t(13;17) translocation is associated with a spectrum of skeletal abnormalities, including Robin sequence, hypoplastic scapulae, and a missing pair of ribs. Using polymerase chain reaction (PCR) with chromosome-specific markers to analyze DBA from somatic cell hybrids containing the derivative translocation chromosomes, we narrowed the breakpoint on each chromosome. Subsequent sequencing of PCR products spanning the breakpoints identified the breaks precisely. The chromosome 17 breakpoint maps approximately 932 kb upstream of the sex-determining region Y (SRY)-related high-mobility group box gene (SOX) within a noncoding transcript represented by two IMAGE cDNA clones. A growing number of reports have implicated chromosome 17 breakpoints at a distance of up to 1 Mb from SOX9 in some cases of campomelic dysplasia (CD). Although this multigeneration family has a disorder that shares some features with CD, their phenotype is significantly milder than any reported cases of (nonmosaic) CD. Therefore, this case may represent an etiologically distinct skeletal dysplasia or may be an extremely mild familial example of CD, caused by the most proximal translocation breakpoint from SOX9 reported to date. In addition, we have refined the breakpoint in a acampomelic CD case described elsewhere and have found that it lies approximately 900 kb upstream of SOX9.

Comparative genomics using fugu: a tool for the identification of conserved vertebrate cis-regulatory elements

With the imminent completion of the whole genome sequence of humans, increasing attention is being focused on the annotation of cis-regulatory elements in the human genome. Comparative genomics approaches based on evolutionary conservation have proved useful in the detection of conserved cis-regulatory elements. The pufferfish, Fugu rubripes, is an attractive vertebrate model for comparative genomics, by virtue of its compact genome and maximal phylogenetic distance from mammals. Fugu has lost a large proportion of nonessential DNA, and retained single orthologs for many duplicate genes that arose in the fish lineage. Non-coding sequences conserved between fugu and mammals have been shown to be functional cis-regulatory elements. Thus, fugu is a model fish genome of choice for discovering evolutionarily conserved regulatory elements in the human genome. Such evolutionarily conserved elements are likely to be shared by all vertebrates, and related to regulatory interactions fundamental to all vertebrates. The functions of these conserved vertebrate elements can be rapidly assayed in mammalian cell lines or in transgenic systems such as zebrafish/medaka and Xenopus, followed by validation of crucial elements in transgenic rodents.

A pair of Sox: distinct and overlapping functions of zebrafish sox9 co-orthologs in craniofacial and pectoral fin development

Understanding how developmental systems evolve after genome amplification is important for discerning the origins of vertebrate novelties, including neural crest, placodes, cartilage and bone. Sox9 is important for the development of these features, and zebrafish has two co-orthologs of tetrapod SOX9 stemming from an ancient genome duplication event in the lineage of ray-fin fish. We have used a genotype-driven screen to isolate a mutation deleting sox9b function, and investigated its phenotype and genetic interactions with a sox9a null mutation. Analysis of mutant phenotypes strongly supports the interpretation that ancestral gene functions partitioned spatially and temporally between Sox9 co-orthologs. Distinct subsets of the craniofacial skeleton, otic placode and pectoral appendage express each gene, and are defective in each single mutant. The double mutant phenotype is additive or synergistic. Ears are somewhat reduced in each single mutant but are mostly absent in the double mutant. Loss-of-function animals from mutations and morpholino injections, and gain-of-function animals injected with sox9a and sox9b mRNAs showed that sox9 helps regulate other early crest genes, including foxd3, sox10, snai1b and crestin, as well as the cartilage gene col2a1 and the bone gene runx2a; however, tfap2a was nearly unchanged in mutants. Chondrocytes failed to stack in sox9a mutants, failed to attain proper numbers in sox9b mutants and failed in both morphogenetic processes in double mutants. Pleiotropy can cause mutations in single copy tetrapod genes, such as Sox9, to block development early and obscure later gene functions. By contrast, subfunction partitioning between zebrafish co-orthologs of tetrapod genes, such as sox9a and sox9b, can relax pleiotropy and reveal both early and late developmental gene functions.

Control of RUNX2 isoform expression: The role of promoters and enhancers

The three mammalian RUNX genes constitute the family of runt domain transcription factors that are involved in the regulation of a number of developmental processes such as haematopoiesis, osteogenesis and neuronal differentiation. All three genes show a complex temporo-spatial pattern of expression. Since the three proteins are probably mutually interchangeable with regard to function, most of the specificity of each family member seems to be based on a tightly controlled regulation of expression. While RUNX gene expression is driven by two promoters for each gene, the promoter sequence alone does not seem to suffice for a proper expressional control. This review focuses on the available evidence for the existence of such control mechanisms and studies aiming at discovering cis-acting regulatory sequences of the RUNX2 gene.

Long-range control of gene expression: emerging mechanisms and disruption in disease

Transcriptional control is a major mechanism for regulating gene expression. The complex machinery required to effect this control is still emerging from functional and evolutionary analysis of genomic architecture. In addition to the promoter, many other regulatory elements are required for spatiotemporally and quantitatively correct gene expression. Enhancer and repressor elements may reside in introns or up- and downstream of the transcription unit. For some genes with highly complex expression patterns--often those that function as key developmental control genes--the cis-regulatory domain can extend long distances outside the transcription unit. Some of the earliest hints of this came from disease-associated chromosomal breaks positioned well outside the relevant gene. With the availability of wide-ranging genome sequence comparisons, strong conservation of many noncoding regions became obvious. Functional studies have shown many of these conserved sites to be transcriptional regulatory elements that sometimes reside inside unrelated neighboring genes. Such sequence-conserved elements generally harbor sites for tissue-specific DNA-binding proteins. Developmentally variable chromatin conformation can control protein access to these sites and can regulate transcription. Disruption of these finely tuned mechanisms can cause disease. Some regulatory element mutations will be associated with phenotypes distinct from any identified for coding-region mutations.