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

Campomelic dysplasia (CD), a human skeletal malformation syndrome with XY sex reversal, is caused by heterozygous mutations in and around the gene SOX9. SOX9 has an extended 5' control region, as indicated by CD translocation breakpoints scattered over 1 Mb proximal to SOX9 and by expression data from mice transgenic for human SOX9-spanning yeast artificial chromosomes. To identify long-range regulatory elements within the SOX9 5' control region, we compared approximately 3.7 Mb and 195 kb of sequence around human and Fugu rubripes SOX9, respectively. We identified only seven and five protein-coding genes in the human and F. rubripes sequences, respectively. Four of the F. rubripes genes have been mapped in humans; all reside on chromosome 17 but show extensive intrachromosomal gene shuffling compared with the gene order in F. rubripes. In both species, very large intergenic distances separate SOX9 from its directly flanking genes: 2 Mb and 500 kb on either side of SOX9 in humans, and 68 and 97 kb on either side of SOX9 in F. rubripes. Comparative sequence analysis of the intergenic regions revealed five conserved elements, E1-E5, up to 290 kb 5' to human SOX9 and up to 18 kb 5' to F. rubripes SOX9, and three such elements, E6-E8, 3' to SOX9. Where available, mouse sequences confirm conservation of the elements. From the yeast artificial chromosome transgenic data, elements E3-E5 are candidate enhancers for SOX9 expression in limb and vertebral column, and 8 of 10 CD translocation breakpoints separate these elements from SOX9.     

Regulation of growth differentiation factor 9 expression in oocytes in vivo: a key role of the E-box

Growth differentiation factor 9 (GDF9) is preferentially expressed in oocytes and is essential for female fertility. To identify regulatory elements that confer high-level expression of GDF9 in the ovary but repression in other tissues, we generated transgenic mice in which regions of the Gdf9 locus were fused to reporter genes. Two transgenes (-10.7/+5.6mGdf9-GFP) and (-3.3/+5.6mGdf9-GFP) that contained sequences either 10.7 or 3.3 kb upstream and 5.6 kb downstream of the Gdf9 initiation codon demonstrated expression specifically in oocytes, thereby mimicking endogenous Gdf9 expression. In contrast, transgenes -10.7mGdf9-Luc and -3.3mGdf9-Luc, which lacked the downstream 5.6-kb region, demonstrated reporter expression not only in oocytes but also high expression in male germ cells. This suggests that the downstream 5.6-kb sequence contains a testis-specific repressor element and that 3.3 kb of 5'-flanking sequence contains all the cis-acting elements for directing high expression of Gdf9 to female (and male) germ cells. To define sequences responsible for oocyte expression of Gdf9, we analyzed sequences of Gdf9 genes from 16 mammalian species. The approximately 400 proximal base pairs upstream of these Gdf9 genes are highly conserved and contain a perfectly conserved E-box (CAGCTG) sequence. When this 400-bp region was placed upstream of a luciferase reporter (-0.4mGdf9-Luc), oocyte-specific expression was observed. However, a similar transgene construct (-0.4MUT-mGdf9-Luc) with a mutation in the E-box abolished oocyte expression. Likewise, the presence of an E-box mutation in a longer construct (-3.3MUT-mGdf9-Luc) abolished expression in the ovary but not in the testis. These observations indicate that the E-box is a key regulatory sequence for Gdf9 expression in the ovary.     

Detection of potential GDF6 regulatory elements by multispecies sequence comparisons and identification of a skeletal joint enhancer

The identification of noncoding functional elements within vertebrate genomes, such as those that regulate gene expression, is a major challenge. Comparisons of orthologous sequences from multiple species are effective at detecting highly conserved regions and can reveal potential regulatory sequences. The GDF6 gene controls developmental patterning of skeletal joints and is associated with numerous, distant cis-acting regulatory elements. Using sequence data from 14 vertebrate species, we performed novel multispecies comparative analyses to detect highly conserved sequences flanking GDF6. The complementary tools WebMCS and ExactPlus identified a series of multispecies conserved sequences (MCSs). Of particular interest are MCSs within noncoding regions previously shown to contain GDF6 regulatory elements. A previously reported conserved sequence at -64 kb was also detected by both WebMCS and ExactPlus. Analysis of LacZ-reporter transgenic mice revealed that a 440-bp segment from this region contains an enhancer for Gdf6 expression in developing proximal limb joints. Several other MCSs represent candidate GDF6 regulatory elements; many of these are not conserved in fish or frog, but are strongly conserved in mammals.     

Efficient identification of regulatory sequences in the chicken genome by a powerful combination of embryo electroporation and genome comparison

Recently expanded knowledge of gene regulation clearly indicates that the regulatory sequences of a gene, usually identified as enhancers, are widely distributed in the gene locus, revising the classical view that they are clustered in the vicinity of genes. To identify regulatory sequences for Sox2 expression governing early neurogenesis, we scanned the 50-kb region of the chicken Sox2 locus for enhancer activity utilizing embryo electroporation, resulting in identification of a number of enhancers scattered throughout the analyzed genomic span. The 'pan-neural' Sox2 expression in early embryos is actually brought about by the composite activities of five separate enhancers with distinct spatio-temporal specificities. These and other functionally defined enhancers exactly correspond to extragenic sequence blocks that are conspicuously conserved between the chicken and mammalian genomes and that are embedded in sequences with a wide range of sequence conservation between humans and mice. The sequences conserved between amniotes and teleosts correspond to subregions of the enhancer subsets which presumably represent core motifs of the enhancers, and the limited conservation partly reflects divergent expression patterns of the gene. The phylogenic distance between the chicken and mammals appears optimal for identifying a battery of genetic regulatory elements as conserved sequence blocks, and chicken embryo electroporation facilitates functional characterization of these elements.     

斑马鱼肌肉高表达基因近端非编码元件分析及功能检测

为鉴定鱼类肌肉组织特异性顺式调控元件,通过分析斑马鱼多个组织的转录组数据,筛选出肌肉高表达基因及低表达基因.通过MEME对肌肉高表达基因和低表达基因非编码区序列特征进行分析,在5个肌肉高表达基因的转录起始位点上游发现了序列保守的DNA区域,包含6个排列顺序一致的DNA基序.将其中一段目标片段插入具有Tol2转座子元件的...     

The evolutionary origin of nodal-related genes in teleosts

Because of an extra whole-genome duplication, zebrafish and other teleosts have two copies of genes that are present in a single copy in tetrapod genomes. Some zebrafish genes, however, are present in triplicate. For example, the nodal-related genes encode secreted proteins of the transforming growth factor beta superfamily that are required in all vertebrates to induce the mesoderm and endoderm, pattern all three germ layers, and establish the left-right axis. Zebrafish have three nodal-related genes, called ndr1/squint, ndr2/cyclops, and ndr3/southpaw. As part of an analysis of enhancer elements controlling zebrafish nodal-related gene expression, we analyzed the nodal loci in the sequenced genomes of five teleost species and four tetrapod species. Each teleost genome contains three nodal-related genes, indicating that squint, cyclops, and southpaw orthologues were present early in the teleost lineage. The genes flanking the nodal-related genes are also conserved, demonstrating a high degree of conserved synteny. Although we found little homology outside of the coding sequences in this region, pufferfish enhancer sequences work in zebrafish embryos to drive reporter gene expression in the squint expression pattern. This indicates a high degree of functional conservation of enhancer elements within the teleosts. We conclude that the ancestral squint and cyclops genes arose during the teleost-specific whole-genome duplication event and that southpaw emerged from a subsequent duplication event involving ancestral squint.     

Structure, organization, and expression of the rat cardiac myosin light chain-2 gene. Identification of a 250-base pair fragment which confers cardiac-specific expression

The present study characterized the structure, organization, and expression of the rat cardiac myosin light chain (MLC) -2 gene. The rat cardiac MLC-2 gene has seven exons which display complete conservation with the exon structure of the rat fast twitch skeletal MLC-2 gene. A 250-base pair (bp) sequence of the 5'-flanking region contains CArG motifs and additional cis elements, each greater than 10 bp in length, which were conserved in sequence and relative position with the chick cardiac MLC-2 gene. A series of MLC-2/luciferase fusion genes consisting of nested 5' deletions of the MLC-2 5'-flanking region were constructed and transfected into primary neonatal rat myocardial cells and a non-myocardial cell line (CV-1), demonstrating that this 250 bp of the MLC-2 5'-flanking region was sufficient to confer cardiac specific expression on a luciferase reporter gene. This study suggests the presence of important proximal regulatory sequences in the MLC-2 5'-flanking region which are capable of directing the cardiac specific expression of the rat cardiac myosin light chain-2 gene.     

Comparison of diverged Hoxc8 early enhancer activities reveals modification of regulatory interactions at conserved cis-acting elements

The Hoxc8 early enhancer that controls the initiation and establishment of Hoxc8 expression in the developing mouse embryo is found in different vertebrate lineages including mammals, birds and fish. Mouse and Fugu Hoxc8 early enhancers (200 bp) have diverged in the composition of elements located towards the 3' region. However, they share cis-acting elements A-E located in the 5' region. Mutations at these elements in the context of the mouse Hoxc8 early enhancer affect reporter gene expression in the posterior neural tube, somites and lateral plate mesoderm of day 9.5 mouse embryos. Here, we demonstrate that mutations introduced at the same elements but in the context of the Fugu Hoxc8 early enhancer had different consequences on the reporter gene expression in transgenic mouse embryos. Furthermore, in contrast to the mouse enhancer the Fugu enhancer does not utilize elements D and E in achieving posterior neural tube and somite expression. These results suggest that the diverged sequences prevent regulatory interactions at conserved cis-acting elements. We propose that divergent sequences modify regulatory interactions at conserved elements by providing a "contextual change". Our finding that the enhancer elements do not act in a unitary fashion but function in the context of the surrounding sequence brings a new dimension to the study of cis-regulatory evolution.     

2S storage protein gene of Douglas-fir: characterization and activity of promoter in transgenic tobacco seeds.

To date a few sequences regulating expression of conifer seed-specific genes have been reported. To characterize Douglas-fir (Pseudotsuga menziesii [Mirb] Franco) 2S albumin storage protein genes, a genomic DNA sequence containing upstream promoter sequences was isolated by screening a Douglas-fir genomic library. Sequence analysis of the Douglas-fir gPm2S1 promoter revealed the presence of RY-repeated elements (GCATGC), and multiple E-box motifs (CANNTG) and ACGT-core elements, features characteristic of 2S storage protein genes in angiosperms. When fused to the GUS reporter gene, the 1.16 kb Douglas-fir 2S promoter sequence was sufficient to direct transient expression in both developing Douglas-fir embryos and maternally derived haploid megagametophytes. Analysis of this promoter construct in transgenic tobacco showed that expression was restricted to embryo and endosperm in developing seeds and was not detected in vegetative tissues of two-week-old seedlings. These results strongly suggest that both structural and regulatory elements as well as upstream signaling components controlling the expression of 2S albumin genes are highly conserved during evolution.     

Surveying phylogenetic footprints in large gene clusters: applications to Hox cluster duplications

Evolutionarily conserved non-coding genomic sequences represent a potentially rich source for the discovery of gene regulatory regions. Since these elements are subject to stabilizing selection they evolve much more slowly than adjacent non-functional DNA. These so-called phylogenetic footprints can be detected by comparison of the sequences surrounding orthologous genes in different species. Therefore the loss of phylogenetic footprints as well as the acquisition of conserved non-coding sequences in some lineages, but not in others, can provide evidence for the evolutionary modification of cis-regulatory elements. We introduce here a statistical model of footprint evolution that allows us to estimate the loss of sequence conservation that can be attributed to gene loss and other structural reasons. This approach to studying the pattern of cis-regulatory element evolution, however, requires the comparison of relatively long sequences from many species. We have therefore developed an efficient software tool for the identification of corresponding footprints in long sequences from multiple species. We apply this novel method to the published sequences of HoxA clusters of shark, human, and the duplicated zebrafish and Takifugu clusters as well as the published HoxB cluster sequences. We find that there is a massive loss of sequence conservation in the intergenic region of the HoxA clusters, consistent with the finding in [Chiu et al., PNAS 99 (2002) 5492]. The loss of conservation after cluster duplication is more extensive than expected from structural reasons. This suggests that binding site turnover and/or adaptive modification may also contribute to the loss of sequence conservation.