Comprehensive characterization of the cis-regulatory code responsible for the spatio-temporal expression of olSix3.2in the developing medaka forebrain
© Conte and Bovolenta.; licensee BioMed Central Ltd. 2007
Received: 23 February 2007
Accepted: 6 July 2007
Published: 06 July 2007
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.
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.
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.
Embryonic development is coordinated by networks of evolutionary conserved regulatory genes that encode transcription factors and components of cell signaling pathways, which in many instances are repetitively exploited in space and time to generate appropriate outcomes in target cells.
Progressive specification of the vertebrate prosencephalon indeed follows this rule [1, 2] and requires, among other factors, recurrent use of Six3, which is a member of the Six/sine oculis family of homeobox transcription factors . In all vertebrates, Six3 is expressed from the neurula stage in the anteriormost neural plate and then in its derivatives: the developing eyes and olfactory placodes, the hypothalamic pituitary regions, and the ventral telencephalon. In mouse and chick, this distribution overlaps with that of its closely related homolog, namely Six6 . However, with time Six3 and Six6 expressions progressively segregate to different brain regions, and Six3 - but not Six6 - is additionally expressed in the olfactory bulb, cerebral cortex, hippocampus, midbrain, and cerebellum . Consistent with this expression, Six3-null mice die at birth, lacking most of the head structures anterior to the midbrain, including eyes , and mutations in SIX3 have been found in humans affected by holoprosencephaly and aprosencephaly/atelencephaly [6, 7]. During mammalian lens induction, Six3 is essential in the presumptive lens ectoderm to activate Pax6 and possibly Sox2 expression . In addition, morpholino-based knockdown of the medaka fish Six3 demonstrates the concentration-dependent need for the function of this transcription factor for proximo-distal patterning of the optic vesicles . Biochemical and functional studies have also shown that Six3, as well as Six6, can induce ectopic retinal tissues and control retinal neuroblast proliferation, acting as transcriptional repressors through the interaction with members of the groucho family of transcriptional co-repressors [10–15]. Furthermore, Six3, but not Six6, functionally interacts with the DNA replication inhibitor Geminin, controlling the balance between cell proliferation and differentiation with a mechanism that is independent of transcriptional regulation .
How the activity of Six3 - or that of any other gene with multiple functions during embryo development - is diversified remains to be elucidated. This could be facilitated by defining the precise gene regulatory network that controls its spatio-temporal expression. It is now well established that control of gene expression is executed through sets of cis-regulatory regions within the noncoding DNA of animal genomes. These cis-regulatory modules have variable length and contain clusters of DNA-binding sites for different transcription factors. These modules work as promoter enhancers or silencers and collectively constitute a unique code for the switching on and off of gene activity [17–19].
The experimental definition of the organization of these specific cis-regulatory elements has progressed substantially in both Drosophila and sea urchin . In contrast, our understanding of how these modules are combined to generate precise gene expression patterns in vertebrates is still rather limited. Possible causes of this are the increased genome complexity and the slow and laborious process of testing the functional significance of identified elements in mammals . Recently, however, computational approaches based on multispecies genomic sequence alignments, combining both closely related and highly divergent organisms, have facilitated identification of highly conserved noncoding sequences, which in many cases appear to coincide with the regulatory modules of genes that play critical roles in development. Analyses of the complex regulation of genes such as Sox2, Sox9, Otx2, Shh, and Irx provide some illustrative examples [21–27]. Functional testing of 'enhancer' activity has also progressed, thanks to the use of alternative and relatively faster 'transgenic' approaches based on the use of nonmammalian vertebrate model systems [20, 25].
Here, we have taken advantage of both the power of computational analysis and the particular compact genome and high transgenesis efficiency of the medaka fish (Oryzia latipes)  to dissect the regulatory control of one of the two Six3 medaka homologs, olSix3.2, that we identified during the course of this study. olSix3.2 is more closely related to the mammalian Six3 than the previously described medaka homolog  (hereafter referred to as 'olSix3.1'). Similar to other related studies [23–25], we identified and functionally characterized sets of cis-regulatory modules that control the olSix3.2 promoter, showing that at least some of these cis-regulatory elements are conserved in other vertebrates, although they are dispersed over a greater stretch of DNA. Going a step further, we have also used combinations and deletions of the identified cis-regulatory modules to elucidate the regulatory code of olSix3.2, which is composed of two enhancers, two silencers, and two 'silencer blockers' used in a combinatorial manner. This comprehensive description of the olSix3.2 cis-regulatory code provides a unique framework for defining the network of trans-acting factors that control the evolutionary conserved activity of Six3 during forebrain development.
Isolation, characterization, and expression of olSix3.2
In order to identify the elements that regulate Six3 expression using the medaka fish (Oryzia latipes) as a model, we used the available olSix3.1 coding sequence (AJ000937) as a query to search public databases (see Materials and methods, below) for the ortholog genomic loci of the closely related species Fugu rubripes, Tetraodon nigroviridis, and Danio rerio (zebrafish). This search retrieved four different loci, one for the fugu and the tetraodon, and two for the zebrafish (six3a and six3b). Alignment of about 20 kilobases (kb) of the retrieved sequences upstream of the Six3 translational start sites identified a cluster of conserved noncoding blocks roughly contained within the first 4.5 kb (data not shown). In the case of the zebrafish, alignment of the six3a or six3b loci yielded comparable results. This information was used to amplify from genomic DNA a fragment of the medaka Six3 locus that contains the corresponding conserved noncoding blocks and the entire first exon.
Interestingly, nucleotide and amino acid sequence alignment of the partially amplified olSix3 coding region did not completely overlap with that reported for the previously identified olSix3.1  but identified - as in zebrafish and Xenopus [30, 31] - a second Six3-related gene in the medaka genome, namely olSix3.2 (AM494407).
Cloning and sequencing of the entire olSix3.2 coding region revealed a two-exon structure, similar to that of olSix3.1 and the mouse Six3, in which the first exon encodes the Six and homeobox domains. olSix3.1 and olSix3.2 exhibited 76% and 63% identity at the nucleotide and amino acid levels, respectively. Interestingly, comparison of the amino acid sequence (81% versus 59%; Additional data file 1) and genomic organization, together with phylogenetic analysis (Additional data file 2), demonstrated that olSix3.2 was more closely related to the mammalian Six3 than the previously identified olSix3.1, which instead falls in between the Six3 and Six6 branches of the family (Additional data file 2).
In conclusion, the distribution of olSix3.2 appeared closely related to that reported for the chick and mouse Six3 [4, 32, 33], whereas the combined expression patterns of olSix3.1 and olSix6 resembled that reported for Six6 [34, 35].
The cis-regulatory elements responsible for olSix3.2expression are contained in a 4.5 kb genomic region ending with a distal 'silencer'
Owing to selective pressure, functional elements in genomes evolve at a slower pace than nonfunctional regions [36–39]. A number of recent studies have functionally demonstrated that a proportion of the highly conserved noncoding regions present in vertebrate genomes correspond to regulatory elements with enhancer activity [21, 39]. We therefore asked whether the region containing the cluster of ten highly conserved noncoding elements was necessary and sufficient to control the entire expression of olSix3.2.
To this end we fused this 4.5 kb genomic region, including the first nine nucleotides of the coding sequence, in frame with a nuclear EGFP (enhanced green fluorescent protein) reporter (Figure 2b). This construct, containing the ten conserved noncoding blocks (termed A-L; Figure 2b), was used to generate three independent stable transgenic medaka lines, which all exhibited a spatio-temporal distribution of the reporter virtually identical to that observed for the endogenous olSix3.2 both at embryonic (compare Figure 1 with Figure 2c-n) and adult stages (not shown). We thus concluded that this region was sufficient to control the entire expression of olSix3.2.
In addition to regulatory elements, sequence conservation could reflect the existence of natural anti-sense mRNAs  or of alternative and yet uncharacterized exons of Six3. However, reverse transcription polymerase chain reaction (RT-PCR) analysis and in situ hybridization studies excluded these possibilities (data not shown). We thus assumed that the ten modules, identified on the basis of their conservation among teleosts (the precise nucleotide sequence of each module is provided in Additional data file 3), could all potentially contain elements that are involved in the regulation of olSix3.2. To test whether this assumption was correct, we generated a series of constructs (named cI to cXXVII) carrying different combinations of the A-L modules, which were then functionally assayed by generating and analyzing three independent stable transgenic lines for the vast majority of the constructs. In each case, the pattern of expression of the EGFP reporter was compared with that observed with construct I (cI), containing the full 4.5 kb sequence (Figure 2i-n) and was always consistent with that observed in F0 injected embryos.
Early expression of olSix3.2 in the anterior neural structures depends on one enhancer, whereas that in the lens placode requires the additional activity of four cis-regulatory modules
Notably, modules D to H (cXVIII; Figure 4a) were also sufficient to induce caudal expansion of reporter expression, with a pattern identical to that observed in the absence of the A module (Figure 3e-g), indicating that modules I and L do not contribute to this expansion or to early expression of the gene.
During organogenesis, appropriate expression of olSix3.2requires the combined activity of two silencers, one enhancer, and two putative 'silencer blockers'
To determine the minimal region of module I involved in the control of this expression, we engineered five 5' to 3' stepwise deletions covering the entire module (cXXIII to cXXVII; Figure 5a). Notably, deletions two, three, and four resulted in progressive abrogation of EGFP expression in the thalamic, hypothalamic (Figyre 5b-d), and telencephalic regions (not shown). This strongly suggests that module I contains a 5' to 3' organized succession of cis-regulatory elements that control the posterior to anterior spatio-temporal organization of olSix3.2 expression in the developing brain. This interpretation was further supported by the injection of two internal deletion constructs (cXXVIII and cXXIX) in which the stretches of nucleotides apparently responsible for hypothalamic and telencephalic expression were removed from cXXIII (Figure 5a). Indeed, in 11% (close to transgenic efficiency) of the embryos analyzed in F0, EGFP fluorescence was not detected in the telencephalon (cXXVIII; Figure 5f) or in the hypothalamus and telencephalon (cXXIX; Figure 5g), clearly indicating that deleted elements are the main driver of olSix3.2 expression in these regions.
When tested alone, block I did not drive EGFP expression in the differentiating retina, whereas activity of the D block was sufficient to maintain reporter expression only in the prospective neural retina (Figure 4a,d,e). Thus, olSix3.2 expression in the differentiating retina appeared to depend on a combination of modules different from those tested thus far. The search for this code demonstrated that only the combined activity of the E to I modules (cVII; Figure 4a) was effective in supporting EGFP expression in the late developing retina.
Identification and characterization of conserved regions among vertebrate
Altogether these data provide a detailed picture of the regulatory code that governs olSix3.2 expression during eye and brain development in medaka. As summarized in Figure 6, this spatio-temporal code is provided by the combined use of at least seven different modules, all conserved among fishes, with distinct enhancer, silencer, or silencer blocker activities. The next logical question was whether this regulatory organisation was conserved in the Six3 locus of vertebrates other than fishes.
Although the human construct (h-cI) we injected drove EGFP expression only in the late olSix3.2 expression domain, according to what expected given the sole presence of modules G to L, we could not exclude that this human region contained regulatory information not readable in fish. Thus, to rule out possible cross-species interferences, we amplified from genomic DNA the equivalent Xenopus region, in which the G to L elements are organized as in humans (Figure 7a). Transgenesis analysis in Xenopus embryos using a construct containing this fragment (X-cI; Figure 7c) yielded results equivalent to those observed with the human fragment; EGFP reporter expression was detected only at later stages of brain development in the expected domain of Xsix3.2 expression (Figure 7j,k). This supports the idea that the regulatory information for early Six3 expression in vertebrates other than fishes reside in as yet unidentified genomic regions.
Six3 is an important regulator of vertebrate forebrain development. Gene regulatory network models predict that the precise spatio-temporal expression pattern of genes fundamental for embryo development must be orchestrated by the interaction of various regulatory regions . Supporting the model, we functionally demonstrated that the entire expression of the newly identified olSix3.2 is orchestrated by the combined use of seven different cis-regulatory modules (Figure 6) and that at least part of this regulation is conserved in the Six3 locus of vertebrates other than fishes. Two main 'enhancer' modules (D and I) are responsible for olSix3.2 expression at early and late stages of brain development, respectively. Their activity is spatially refined by the function of two 'silencers' and two 'silencer blockers'. In addition, olSix3.2 expression in the lens ectoderm and in the differentiating retina requires the combined activity of five different cis-regulatory modules. This apparently simple regulation may hide additional organization, as we have demonstrated for the I enhancer, in which an organized sequence of cis-regulatory elements control the posterior to anterior expression of olSix3.2 in the brain.
The availability of different genome sequences and the development of analytical bioinformatic tools have facilitated study of cis regulation of a number of genes with evolutionary conserved roles in vertebrate embryonic development. Some of these studies have focused, as has ours, on a specific gene or a gene cluster, identifying enhancers that are involved in the control of specific expression domains [21, 23, 25, 45–49]. However, possibly because of the size of the genomic regions that are involved, or to the laborious and time consuming use of mice, or the limitations of chick electroporation in validating regulatory activities, these studies have mostly focused on each enhancer as a separate entity, thus missing the effects of possible cooperative activities. Other recent and extremely informative studies, based on medium or small throughput screens in zebrafish, have instead systematically tested the autonomously enhancing function of large numbers of highly conserved noncoding elements positioned in areas surrounding developmentally important genes, with positive identification only of a fraction of them [21, 39]. Because each element is tested in an unconstrained context, negative regulators as well as modulatory functions of surrounding endogenous elements are also undetected using these approaches . In contrast, possibly benefiting from the high transgenesis efficiency of the medaka fish  and its compact genome, we were able to assign enhancer, silencer, and modulatory functions to the majority of the highly conserved noncoding elements surrounding the Six3 gene in fishes. Testing different combinations of these elements, we have also established their required interactions for proper expression of the gene. Thus, to our knowledge, we provide the first description of the regulatory code necessary for the expression of a vertebrate gene and offer a unique framework to define the entire interplay of trans-acting factors that control the evolutionary conserved use of Six3 during forebrain development.
Teleosts are the most diverse class of vertebrates with a huge variety of different species; they are characterized by broad size range and dynamic organization of genomes, which are the result of an initial genome duplication followed by subsequent independent evolution of the different lineages [51, 52]. Comparison of divergent teleost genomes largely separated in the phylogenetic tree, such as the medaka and zebrafish genomes (approximately 115 to 200 million years ), is thus a powerful tool with which to study gene regulatory mechanisms. Adopting this strategy, we identified a cluster of potential regulatory modules in the Six3 locus, which were barely identifiable in a comparison among mammalian genomes (compare Figure 2a with Figure 7a). VISTA analysis of the available genomic sequences flanking the homologous vertebrate Six3 genes revealed several blocks of highly conserved noncoding sequences in the gene surroundings. Although a few of these blocks were located downstream of the coding sequence (data not shown), we demonstrated that the pattern of olSix3.2 expression could be recapitulated by 4.5 kb of genomic sequence flanking the 5' end of the gene. This conclusion is based on a relatively efficient (roughly 20% of injected embryos) and highly reproducible (basically 100%, albeit with different EGFP intensity, thus excluding chromosomal position effects) transgenic analysis using three independent and stable medaka lines generated for all of the constructs we tested. Thus, we are fairly confident that we identified the main regulatory region for olSix3.2, although we cannot entirely exclude the possibility that additional or duplicated regulatory elements positioned in untested regions may contribute to a refinement of the main expression domain. Indeed, redundant cis-regulatory elements have been reported to control specific expression domains in different genes, including Otx2, Shh, and Sox2 [22, 23, 25].
According to our analysis, the regulatory region of olSix3.2 is relatively compact as compared with those reported for other genes that are involved in neural development, such as Sox2, Sox9, Otx2, Pax6, and Shh, for which enhancers located in the region of 10, 100, and even 1,000 kb away from their promoters have been reported [23–27, 53, 54], even in the compact Fugu genome . Genes with complex patterns of expression are predicted to have more regulatory elements and occupy significantly more space in the genome than those with simpler expressions that are restricted to populations of cells with similarities or shared identity . It is thus possible that the compactness of the olSix3.2 regulatory region might reflect the association that exists among the main territories in which the gene is expressed. Indeed, the specification of telencephalic and eye fields appears to be closely linked , and the initial expression of olSix3.2 in both regions appears to depend on the activity of a single enhancer element (D) and a distal silencer (A), which constrains the expression domain to the anteriormost neural tube. This hypothesis could also explain why the combined activities of five different modules (D to H) are instead needed to control expression in the lens placode, which is the only non-neural domain of olSix3.2 expression. Nevertheless, compactness does not appear to be, at least in this case, a reflection of simplicity, because each of the conserved modules may include additional regulatory organization. This is the case of module I, which is the main enhancer involved in the late embryonic expression of the gene. Stepwise and internal deletions of this module have revealed a peculiar organization, in a 5' to 3' direction, of a series of cis-regulatory elements that are required for the posterior to anterior spatio-temporal expression of olSix3.2 in the thalamus, hypothalamus, and telencephalon. The activity of the I module is refined by a silencer, G, the activity of which is modulated by two silencer blockers that act in a temporal sequence, thus establishing an elaborate control code. Furthermore, although the L module per se has no activity, we cannot totally exclude the possibility that this module might contribute, together with modules E to H, to the regulation of I, because it was present in the constructs used for this analysis.
Alternatively, the short-range regulation of olSix3.2 may be linked to the chromosomal localization of the Six genes, which are organized in two evolutionarily conserved clusters . Although the expression of the other Six family members (Six1, Six2, Six4, and Six5) is mostly associated with tissues of mesodermal and ectodermal origin , it is possible that genes within the same cluster (Six4, Six1, and Six6) will share a few regulatory elements, which might have imposed constrains against rearrangement during evolution .
In silico comparison identified ten conserved modules in the teleost Six3 locus. Transgenic analysis in medaka demonstrated clear regulatory activity for seven of them, whereas modules B, C, and L did not influence EGFP reporter expression. Although these modules might have subtle regulatory activities below the resolution of our analysis, their conservation could reflect other important roles in gene transcription control, such as regulation of chromatin structure or - in the particular case of module L - they may contribute to minimal promoter functions.
The regulatory region we have studied belongs to a newly identified medaka Six3 gene, namely olSix3.2. Genomic organization and phylogenetic analysis suggests that olSix3.2 is more closely related to the mammalian Six3 than the previously identified olSix3.1 . Like its mammalian homolog , olSix3.2 is strongly expressed in various forebrain regions where its paralog is not expressed. Our comparative expression study suggests that the combination of expression domains of olSix3.1, olSix3.2, and the related olSix6 correspond to the combined tissue distribution observed for the mouse and chick Six3 and Six6 [32–34], with a preponderant expression of olSix3.1 in the eye, of olSix3.2 in the telencephalic and thalamic regions, and of olSix6 in the hypothalamus. Genetic abrogation studies in mice demonstrated that Six3 is necessary for the formation of forebrain, which is absent in homozygous embryos . Genetic deletion of Six6 instead is associated with pituitary defects, absence or hypoplasia of the optic nerves, and chiasm and alteration in neural retina proliferation . How the functions of olSix3.1, olSix3.2, and olSix6 relate to those described in the mouse for Six3 and Six6 is still unresolved and knock-down analysis of all three genes in medaka will be necessary to address this issue. Thus far, morpholino-based knock-down of olSix3.1 results in forebrain and eye defects, including loss of optic stalk markers , whereas preliminary analysis indicates that olSix3.2 morphants are characterized by strong ventral forebrain defects with minor eye malformations (De la Torre A, Conte I, Bovolenta P, unpublished observations), suggesting that the two olSix3 paralogs may cover Six3 as well as part of the mouse Six6 functions, a possibility that is also supported by the phylogenetic position of olSix3.1, which falls almost in between the Six3 and Six6 branches of the Six gene family (Additional data file 2).
Comparative analysis of the regulatory code of the three medaka genes currently ongoing in our laboratory might be useful in complementing these studies by providing insights into the sub-functionalization or neo-functionalization of olSix3.1, olSix3.2, and olSix6 as compared with their mammalian counterparts. Furthermore, they will help to elucidate whether Six3 and Six6 have arisen from the duplication of a common ancestor, as previously proposed , possibly duplicating at least part of their regulatory region. This is an important point because, with the comparison parameters used, we were unable to identify in other vertebrate species the conservation and distribution of the A to F regulatory modules characterized in fishes. This is particularly important for the A and D modules, which are the main regulators of early Six3 expression in fishes. Informatics searches of corresponding regions in mammalian genomes yielded no clear information, suggesting that these modules might be present outside the regions that we analyzed or they might have evolved differently in other vertebrate genomes, making their search even more difficult than that of the H and I modules. Alternatively, these modules may represent a new acquisition of olSix3.2 caused by teleost genome duplication.
In our study, we demonstrated strong functional conservation between fishes and other vertebrates only for the G, H, and I modules, which control late expression of olSix3.2. The sequences that compose the H and I modules in fishes were intermixed and differently arranged in other vertebrate genomes, although their function was strongly conserved when assayed in medaka and Xenopus transgenesis. This suggests that sequences from different vertebrates are activated by common transcription factors, although the binding sites for these factors might be distributed, oriented, or represented in different numbers among species. An additional explanation for the different arrangement of the H and I modules might be species-specific nucleotide modifications, which have been proposed to contribute to gene transcriptional evolution [59–61].
Conservation of regulatory function between human and fish in the absence of clear sequence conservation has previously been reported also for the RET gene. In this case, lack of correlation between the two events was even more marked, and different in silico analysis designed to detect shorter stretches of sequence similarities or the existence of inversion and rearrangement failed to detect alignable sequences . Thus, our data, together with few additional observations [63–65], strongly support the idea proposed by Fisher and colleagues  that some relevant regulatory information might be conserved among species at a level that is not detectable using genomic sequence alignment.
Our study established the cis-regulatory code required for the proper expression of olSix3.2 and demonstrates that there is a need to test different combinations of highly conserved putative cis-regulatory regions to elucidate how each conserved element contributes to the spatio-temporal control of gene expression. In fact, one limitation of previous studies that have used transgenic analysis to test the function of highly conserved noncoding sequences is the identification of single enhancers uprooted from possible interactions with the remaining regulatory elements. Our comprehensive description of the olSix3.2 regulatory code is now a powerful starting point from which to define the entire interplay of trans-acting factors that control the evolutionarily conserved use of Six3 during forebrain development. From a broader perspective, this type of information will be necessary to elucidate the composition and evolution of vertebrate gene regulatory networks, as compared with those of invertebrates such as Drosophila and sea urchin, in which this type of information is accumulating at a much faster pace .
Materials and methods
Microinjection and establishment of transgenic lines
Adult and embryonic medaka fishes (Oryzia latipes) from the Cab inbred strain were used throughout the study. Fertilized eggs were collected immediately and incubated at 4 to 10°C in Yamamoto's embryo rearing medium to suppress further development . DNA was prepared using a High Pure Plasmid Isolation Kit (Roche, Basel, Switzerland). DNA injections (10 ng/μl DNA in ISceI enzyme reaction) were performed as previously described . Embryos were staged according to the method proposed by Iwamatsu , raised to sexual maturity, and transgenic founder fishes were identified by out-crossing to wild-type fishes. Transcriptional activation of the constructs was monitored by EGFP expression observed in living embryos under UV fluorescent stereo-microscopy (Leica Microsystems, Wetzlar, Germany). Xenopus laevis embryos were obtained and raised as described previously . Xenopus transgenesis was performed following the same procedures as used for the medaka embryos.
Whole-mount in situhybridization
Whole-mount in situ hybridizations were performed as previously described using digoxigenin labelled riboprobes . Anti-sense and sense riboprobes for medaka olSix3.1, olSix3.2, and olSix6 and the Xenopus Xsix3.2 were used. A minimum of 40 embryos were hybridized for each marker and condition. In toto hybridized embryos were photographed, embedded in gelatine/albumine block, and further sectioned using a vibratome (Leica Microsystems, Wetzlar, Germany).
The vertebrate Six3 genomic sequences were retrieved from public databases: Genome Browser UCSC  and JGI . The genomic sequence of olSix3.2 was isolated from medaka genomic DNA using the following primers: olSix3 forward CCTCATTAAATGTCGCTAAC, and olSix3 reverse cgcctaatgacac cagcctc. Sequence alignments were performed using the VISTA  and Multalign programs , which are available at the corresponding websites [69, 70]. The criterion used for comparisons was a minimum 75% nucleotide identity with a window size of over 100 bp. Phylogenetic analysis was performed using the PHYLIP package . The results were plotted using the Tree-view software package . olSix3.2 protein sequences were scanned for motifs using online software available at HGMP  and NCBI .
Isolation of olSix3.2cDNA
Total RNAs from medaka embryos at different stages were isolated by RNAzol B (Campro Scientific, Berlin, Germany) and treated with Dnase I (Invitrogen, Carlsbad, CA). RT-PCR reactions were performed using SUPERSCRIPT II (Invitrogen, Carlsbad, CA), as described previously . PCR using olSix3.2 specific primers was performed using 2 μl of the reverse transcription reaction as a template with the High Fidelity PCR system (Roche, Basel, Switzerland). Oligonucleotide primers used to isolate olSix3.2 cDNA are listed in Additional data file 4.
A 4.5 kb region of olSix3.2 genomic sequence containing nine nucleotides (corresponding to the first three amino acids) of the coding region was cloned in frame with EGFP reporter gene into the pSKII-ISceI-EGFP vector , to create the cI construct. Xenopus and human sequences were amplified from corresponding genomic DNA and cloned in the pSKII-ISceI-EGFP vector with the same strategy. The medaka deleted constructs pSix3.2ΔXhoI (cII), pSix3.2ΔXhoI-NsiI (cVIII), and pSix3.2ΔXbaI-HindIII (cIX) were obtained by digesting the pSix3.2-4.5 kb construct using the indicated enzymes. All the other deleted constructs (pSix3.2Δel1, pSix3.2Δel2, pSix3.2Δel3, pSix3.2Δel4, pSix3.2Δel5, pSix3.2Δel6, and pSix3.2Δel7; cXXIII to cXXIX) were obtained by PCR amplification from pSix3.2-4.5 kb and then cloned into pSKII-ISceI-EGFP vector. The A, B, and C modules were deleted by restriction enzyme digestion (A, NarI/KpnI; B, BtsI/BglII; and C, BamHI/ClaI) and inserted (in sense and anti-sense orientations) into the polylinker of pSix3.2ΔXhoI (cIII to cV), pSix3.2Δel1 and pSKII-ISceI-Tk-EGFP (containing the tyrosine kinase minimal promoter) vectors to test their potential regulatory activity. All of the other modules were amplified and cloned (in sense and anti-sense orientations) into the polylinker of pSix3.2Δel1 (cVI to cVII, cX, and cXII to cXVI) and pSKII-ISceI-Tk-EGFP (cXI and cXVII to cXXII). The primer sequences used to generate these constructs are shown in Supplementary Table I. All constructs were verified by automated sequencing.
Additional data files
The following additional data are available with the online version of this manuscript. Additional data file 1 is a Figure reporting the amino acid sequence alignment of Six3 genes from different vertebrate species. Additional data file 2 is a figure illustrating the phylogenetic tree of the SIX family. Additional data file 3 provides the precise nucleotide sequences of modules A to L described in the report. Additional data file 4 is a table listing the sequences of the primers used to amplify the DNA fragments, which were used to design the different constructs described in the report.
We are grateful to Drs WA Harris, H Kondoh, and M Manzanares for critical reading of the manuscript, and to Dr J Wittbrodt and members of our laboratory for many helpful suggestions. We wish to thank I Dompablo for excellent technical assistance. This study was supported by grants from Spanish Ministerio de Educación y Ciencia (BFU-2004-01585) and in part by the EU (QLG3-CT-2001-01460) and the HFSPO (RGP0040/2001-M) to PB. A Telethon Foundation (GFP03007) and MEC (SB2003-0182) fellowships supported the postdoctoral work of IC.
- Wilson SW, Houart C: Early steps in the development of the forebrain. Dev Cell. 2004, 6: 167-181. 10.1016/S1534-5807(04)00027-9.PubMedPubMed CentralView ArticleGoogle Scholar
- Esteve P, Bovolenta P: Secreted inducers in vertebrate eye development: more functions for old morphogens. Curr Opin Neurobiol. 2006, 16: 13-19. 10.1016/j.conb.2006.01.001.PubMedView ArticleGoogle Scholar
- Rodríguez de Córdoba S, Gallardo ME, Lopez-Rios J, Bovolenta P: The human Six family of homeobox genes. Curr Genomics. 2001, 2: 231-242. 10.2174/1389202013350832.View ArticleGoogle Scholar
- Conte I, Morcillo J, Bovolenta P: Comparative analysis of Six 3 and Six 6 distribution in the developing and adult mouse brain. Dev Dyn. 2005, 234: 718-725. 10.1002/dvdy.20463.PubMedView ArticleGoogle Scholar
- Lagutin OV, Zhu CC, Kobayashi D, Topczewski J, Shimamura K, Puelles L, Russell HR, McKinnon PJ, Solnica-Krezel L, Oliver G: Six3 repression of Wnt signaling in the anterior neuroectoderm is essential for vertebrate forebrain development. Genes Dev. 2003, 17: 368-379. 10.1101/gad.1059403.PubMedPubMed CentralView ArticleGoogle Scholar
- Wallis DE, Muenke M: Molecular mechanisms of holoprosencephaly. Mol Genet Metab. 1999, 68: 126-138. 10.1006/mgme.1999.2895.PubMedView ArticleGoogle Scholar
- Pasquier L, Dubourg C, Gonzales M, Lazaro L, David V, Odent S, Encha-Razavi F: First occurrence of aprosencephaly/atelencephaly and holoprosencephaly in a family with a SIX3 gene mutation and phenotype/genotype correlation in our series of SIX3 mutations. J Med Genet. 2005, 42: e4-10.1136/jmg.2004.023416.PubMedPubMed CentralView ArticleGoogle Scholar
- Liu W, Lagutin OV, Mende M, Streit A, Oliver G: Six3 activation of Pax6 expression is essential for mammalian lens induction and specification. EMBO J. 2006, 25: 5383-5395. 10.1038/sj.emboj.7601398.PubMedPubMed CentralView ArticleGoogle Scholar
- Carl M, Loosli F, Wittbrodt J: Six3 inactivation reveals its essential role for the formation and patterning of the vertebrate eye. Development. 2002, 129: 4057-4063.PubMedGoogle Scholar
- Bernier G, Panitz F, Zhou X, Hollemann T, Gruss P, Pieler T: Expanded retina territory by midbrain transformation upon overexpression of Six6 (Optx2) in Xenopus embryos. Mech Dev. 2000, 93: 59-69. 10.1016/S0925-4773(00)00271-9.PubMedView ArticleGoogle Scholar
- Kobayashi M, Nishikawa K, Suzuki T, Yamamoto M: The homeobox protein Six3 interacts with the Groucho corepressor and acts as a transcriptional repressor in eye and forebrain formation. Dev Biol. 2001, 232: 315-326. 10.1006/dbio.2001.0185.PubMedView ArticleGoogle Scholar
- Loosli F, Winkler S, Wittbrodt J: Six3 overexpression initiates the formation of ectopic retina. Genes Dev. 1999, 13: 649-654.PubMedPubMed CentralView ArticleGoogle Scholar
- Lopez-Rios J, Tessmar K, Loosli F, Wittbrodt J, Bovolenta P: Six3 and Six6 activity is modulated by members of the groucho family. Development. 2003, 130: 185-195. 10.1242/dev.00185.PubMedView ArticleGoogle Scholar
- Zhu CC, Dyer MA, Uchikawa M, Kondoh H, Lagutin OV, Oliver G: Six3-mediated auto repression and eye development requires its interaction with members of the Groucho-related family of co-repressors. Development. 2002, 129: 2835-2849.PubMedGoogle Scholar
- Zuber ME, Perron M, Philpott A, Bang A, Harris WA: Giant eyes in Xenopus laevis by overexpression of XOptx2. Cell. 1999, 98: 341-352. 10.1016/S0092-8674(00)81963-7.PubMedView ArticleGoogle Scholar
- Del Bene F, Tessmar-Raible K, Wittbrodt J: Direct interaction of geminin and Six3 in eye development. Nature. 2004, 427: 745-749. 10.1038/nature02292.PubMedView ArticleGoogle Scholar
- Davidson EH, Erwin DH: Gene regulatory networks and the evolution of animal body plans. Science. 2006, 311: 796-800. 10.1126/science.1113832.PubMedView ArticleGoogle Scholar
- Levine M, Davidson EH: Gene regulatory networks for development. Proc Natl Acad Sci USA. 2005, 102: 4936-4942. 10.1073/pnas.0408031102.PubMedPubMed CentralView ArticleGoogle Scholar
- Michelson AM: Deciphering genetic regulatory codes: a challenge for functional genomics. Proc Natl Acad Sci USA. 2002, 99: 546-548. 10.1073/pnas.032685999.PubMedPubMed CentralView ArticleGoogle Scholar
- Allende ML, Manzanares M, Tena JJ, Feijoo CG, Gomez-Skarmeta JL: Cracking the genome's second code: enhancer detection by combined phylogenetic footprinting and transgenic fish and frog embryos. Methods. 2006, 39: 212-219. 10.1016/j.ymeth.2005.12.005.PubMedView ArticleGoogle Scholar
- de la Calle-Mustienes E, Feijoo CG, Manzanares M, Tena JJ, Rodriguez-Seguel E, Letizia A, Allende ML, Gomez-Skarmeta JL: A functional survey of the enhancer activity of conserved non-coding sequences from vertebrate Iroquois cluster gene deserts. Genome Res. 2005, 15: 1061-1072. 10.1101/gr.4004805.PubMedPubMed CentralView ArticleGoogle Scholar
- Jeong Y, El-Jaick K, Roessler E, Muenke M, Epstein DJ: A functional screen for sonic hedgehog regulatory elements across a 1 Mb interval identifies long-range ventral forebrain enhancers. Development. 2006, 133: 761-772. 10.1242/dev.02239.PubMedView ArticleGoogle Scholar
- Kurokawa D, Kiyonari H, Nakayama R, Kimura-Yoshida C, Matsuo I, Aizawa S: Regulation of Otx2 expression and its functions in mouse forebrain and midbrain. Development. 2004, 131: 3319-3331. 10.1242/dev.01220.PubMedView ArticleGoogle Scholar
- Kurokawa D, Takasaki N, Kiyonari H, Nakayama R, Kimura-Yoshida C, Matsuo I, Aizawa S: Regulation of Otx2 expression and its functions in mouse epiblast and anterior neuroectoderm. Development. 2004, 131: 3307-3317. 10.1242/dev.01219.PubMedView ArticleGoogle Scholar
- Uchikawa M, Ishida Y, Takemoto T, Kamachi Y, Kondoh H: Functional analysis of chicken Sox2 enhancers highlights an array of diverse regulatory elements that are conserved in mammals. Dev Cell. 2003, 4: 509-519. 10.1016/S1534-5807(03)00088-1.PubMedView ArticleGoogle Scholar
- Bagheri-Fam S, Barrionuevo F, Dohrmann U, Gunther T, Schule R, Kemler R, Mallo M, Kanzler B, Scherer G: Long-range upstream and downstream enhancers control distinct subsets of the complex spatiotemporal Sox9 expression pattern. Dev Biol. 2006, 291: 382-397. 10.1016/j.ydbio.2005.11.013.PubMedView ArticleGoogle Scholar
- Bagheri-Fam S, Ferraz C, Demaille J, Scherer G, Pfeifer D: Comparative genomics of the SOX9 region in human and Fugu rubripes: conservation of short regulatory sequence elements within large intergenic regions. Genomics. 2001, 78: 73-82. 10.1006/geno.2001.6648.PubMedView ArticleGoogle Scholar
- Wittbrodt J, Shima A, Schartl M: Medaka: a model organism from the far East. Nat Rev Genet. 2002, 3: 53-64. 10.1038/nrg704.PubMedView ArticleGoogle Scholar
- Loosli F, Koster RW, Carl M, Krone A, Wittbrodt J: Six3, a medaka homologue of the Drosophila homeobox gene sine oculis is expressed in the anterior embryonic shield and the developing eye. Mech Dev. 1998, 74: 159-164. 10.1016/S0925-4773(98)00055-0.PubMedView ArticleGoogle Scholar
- Ghanbari H, Seo HC, Fjose A, Brandli AW: Molecular cloning and embryonic expression of Xenopus Six homeobox genes. Mech Dev. 2001, 101: 271-277. 10.1016/S0925-4773(00)00572-4.PubMedView ArticleGoogle Scholar
- Seo HC, Drivenes , Ellingsen S, Fjose A: Expression of two zebrafish homologues of the murine Six3 gene demarcates the initial eye primordia. Mech Dev. 1998, 73: 45-57. 10.1016/S0925-4773(98)00028-8.PubMedView ArticleGoogle Scholar
- Bovolenta P, Mallamaci A, Puelles L, Boncinelli E: Expression pattern of cSix3, a member of the Six/sine oculis family of transcription factors. Mech Dev. 1998, 70: 201-203. 10.1016/S0925-4773(97)00183-4.PubMedView ArticleGoogle Scholar
- Oliver G, Mailhos A, Wehr R, Copeland NG, Jenkins NA, Gruss P: Six3, a murine homologue of the sine oculis gene, demarcates the most anterior border of the developing neural plate and is expressed during eye development. Development. 1995, 121: 4045-4055.PubMedGoogle Scholar
- Lopez-Rios J, Gallardo ME, Rodriguez de Cordoba S, Bovolenta P: Six9 (Optx2), a new member of the six gene family of transcription factors, is expressed at early stages of vertebrate ocular and pituitary development. Mech Dev. 1999, 83: 155-159. 10.1016/S0925-4773(99)00017-9.PubMedView ArticleGoogle Scholar
- Toy J, Yang JM, Leppert GS, Sundin OH: The optx2 homeobox gene is expressed in early precursors of the eye and activates retina-specific genes. Proc Natl Acad Sci USA. 1998, 95: 10643-10648. 10.1073/pnas.95.18.10643.PubMedPubMed CentralView ArticleGoogle Scholar
- Hardison RC: Conserved noncoding sequences are reliable guides to regulatory elements. Trends Genet. 2000, 16: 369-372. 10.1016/S0168-9525(00)02081-3.PubMedView ArticleGoogle Scholar
- Thomas JW, Touchman JW, Blakesley RW, Bouffard GG, Beckstrom-Sternberg SM, Margulies EH, Blanchette M, Siepel AC, Thomas PJ, McDowell JC, et al: Comparative analyses of multi-species sequences from targeted genomic regions. Nature. 2003, 424: 788-793. 10.1038/nature01858.PubMedView ArticleGoogle Scholar
- Nobrega MA, Pennacchio LA: Comparative genomic analysis as a tool for biological discovery. J Physiol. 2004, 554: 31-39. 10.1113/jphysiol.2003.050948.PubMedPubMed CentralView ArticleGoogle Scholar
- Woolfe A, Goodson M, Goode DK, Snell P, McEwen GK, Vavouri T, Smith SF, North P, Callaway H, Kelly K, et al: Highly conserved non-coding sequences are associated with vertebrate development. PLoS Biol. 2005, 3: e7-10.1371/journal.pbio.0030007.PubMedPubMed CentralView ArticleGoogle Scholar
- Alfano G, Vitiello C, Caccioppoli C, Caramico T, Carola A, Szego MJ, McInnes RR, Auricchio A, Banfi S: Natural antisense transcripts associated with genes involved in eye development. Hum Mol Genet. 2005, 14: 913-923. 10.1093/hmg/ddi084.PubMedView ArticleGoogle Scholar
- Oki M, Kamakaka RT: Blockers and barriers to transcription: competing activities?. Curr Opin Cell Biol. 2002, 14: 299-304. 10.1016/S0955-0674(02)00327-7.PubMedView ArticleGoogle Scholar
- West AG, Fraser P: Remote control of gene transcription. Hum Mol Genet. 2005, 14 (Spec No 1): R101-R111. 10.1093/hmg/ddi104.PubMedView ArticleGoogle Scholar
- Frazer KA, Pachter L, Poliakov A, Rubin EM, Dubchak I: VISTA: computational tools for comparative genomics. Nucleic Acids Res. 2004, W273-W279. 10.1093/nar/gkh458. 32 Web ServerGoogle Scholar
- Corpet F: Multiple sequence alignment with hierarchical clustering. Nucleic Acids Res. 1988, 16: 10881-10890. 10.1093/nar/16.22.10881.PubMedPubMed CentralView ArticleGoogle Scholar
- Kimura C, Takeda N, Suzuki M, Oshimura M, Aizawa S, Matsuo I: Cis-acting elements conserved between mouse and pufferfish Otx2 genes govern the expression in mesencephalic neural crest cells. Development. 1997, 124: 3929-3941.PubMedGoogle Scholar
- Kimura-Yoshida C, Kitajima K, Oda-Ishii I, Tian E, Suzuki M, Yamamoto M, Suzuki T, Kobayashi M, Aizawa S, Matsuo I: Characterization of the pufferfish Otx2 cis-regulators reveals evolutionarily conserved genetic mechanisms for vertebrate head specification. Development. 2004, 131: 57-71. 10.1242/dev.00877.PubMedView ArticleGoogle Scholar
- Takemoto T, Uchikawa M, Kamachi Y, Kondoh H: Convergence of Wnt and FGF signals in the genesis of posterior neural plate through activation of the Sox2 enhancer N-1. Development. 2006, 133: 297-306. 10.1242/dev.02196.PubMedView ArticleGoogle Scholar
- Hutcheson DA, Hanson MI, Moore KB, Le TT, Brown NL, Vetter ML: bHLH-dependent and -independent modes of Ath5 gene regulation during retinal development. Development. 2005, 132: 829-839. 10.1242/dev.01653.PubMedView ArticleGoogle Scholar
- Lee AP, Koh EG, Tay A, Brenner S, Venkatesh B: Highly conserved syntenic blocks at the vertebrate Hox loci and conserved regulatory elements within and outside Hox gene clusters. Proc Natl Acad Sci USA. 2006, 103: 6994-6999. 10.1073/pnas.0601492103.PubMedPubMed CentralView ArticleGoogle Scholar
- Thermes V, Grabher C, Ristoratore F, Bourrat F, Choulika A, Wittbrodt J, Joly JS: I-SceI meganuclease mediates highly efficient transgenesis in fish. Mech Dev. 2002, 118: 91-98. 10.1016/S0925-4773(02)00218-6.PubMedView ArticleGoogle Scholar
- Furutani-Seiki M, Wittbrodt J: Medaka and zebrafish, an evolutionary twin study. Mech Dev. 2004, 121: 629-637. 10.1016/j.mod.2004.05.010.PubMedView ArticleGoogle Scholar
- Muller F, Blader P, Strahle U: Search for enhancers: teleost models in comparative genomic and transgenic analysis of cis regulatory elements. Bioessays. 2002, 24: 564-572. 10.1002/bies.10096.PubMedView ArticleGoogle Scholar
- Kleinjan DA, Seawright A, Mella S, Carr CB, Tyas DA, Simpson TI, Mason JO, Price DJ, van Heyningen V: Long-range downstream enhancers are essential for Pax6 expression. Dev Biol. 2006, 299: 563-581. 10.1016/j.ydbio.2006.08.060.PubMedPubMed CentralView ArticleGoogle Scholar
- Lettice LA, Heaney SJ, Purdie LA, Li L, de Beer P, Oostra BA, Goode D, Elgar G, Hill RE, de Graaff E: A long-range Shh enhancer regulates expression in the developing limb and fin and is associated with preaxial polydactyly. Hum Mol Genet. 2003, 12: 1725-1735. 10.1093/hmg/ddg180.PubMedView ArticleGoogle Scholar
- Nelson CE, Hersh BM, Carroll SB: The regulatory content of intergenic DNA shapes genome architecture. Genome Biol. 2004, 5: R25-10.1186/gb-2004-5-4-r25.PubMedPubMed CentralView ArticleGoogle Scholar
- Gallardo ME, Lopez-Rios J, Fernaud-Espinosa I, Granadino B, Sanz R, Ramos C, Ayuso C, Seller MJ, Brunner HG, Bovolenta P, et al: Genomic cloning and characterization of the human homeobox gene SIX6 reveals a cluster of SIX genes in chromosome 14 and associates SIX6 hemizygosity with bilateral anophthalmia and pituitary anomalies. Genomics. 1999, 61: 82-91. 10.1006/geno.1999.5916.PubMedView ArticleGoogle Scholar
- Duboule D, Wilkins AS: The evolution of 'bricolage'. Trends Genet. 1998, 14: 54-59. 10.1016/S0168-9525(97)01358-9.PubMedView ArticleGoogle Scholar
- Li X, Perissi V, Liu F, Rose DW, Rosenfeld MG: Tissue-specific regulation of retinal and pituitary precursor cell proliferation. Science. 2002, 297: 1180-1183.PubMedGoogle Scholar
- Cowles CR, Hirschhorn JN, Altshuler D, Lander ES: Detection of regulatory variation in mouse genes. Nat Genet. 2002, 32: 432-437. 10.1038/ng992.PubMedView ArticleGoogle Scholar
- Wray GA: Transcriptional regulation and the evolution of development. Int J Dev Biol. 2003, 47: 675-684.PubMedGoogle Scholar
- Teichmann SA, Babu MM: Gene regulatory network growth by duplication. Nat Genet. 2004, 36: 492-496. 10.1038/ng1340.PubMedView ArticleGoogle Scholar
- Fisher S, Grice EA, Vinton RM, Bessling SL, McCallion AS: Conservation of RET regulatory function from human to zebrafish without sequence similarity. Science. 2006, 312: 276-279. 10.1126/science.1124070.PubMedView ArticleGoogle Scholar
- Oda-Ishii I, Bertrand V, Matsuo I, Lemaire P, Saiga H: Making very similar embryos with divergent genomes: conservation of regulatory mechanisms of Otx between the ascidians Halocynthia roretzi and Ciona intestinalis. Development. 2005, 132: 1663-1674. 10.1242/dev.01707.PubMedView ArticleGoogle Scholar
- Pennacchio LA, Rubin EM: Genomic strategies to identify mammalian regulatory sequences. Nat Rev Genet. 2001, 2: 100-109. 10.1038/35052548.PubMedView ArticleGoogle Scholar
- Gottgens B, Barton LM, Gilbert JG, Bench AJ, Sanchez MJ, Bahn S, Mistry S, Grafham D, McMurray A, Vaudin M, et al: Analysis of vertebrate SCL loci identifies conserved enhancers. Nat Biotechnol. 2000, 18: 181-186. 10.1038/72635.PubMedView ArticleGoogle Scholar
- Iwamatsu T: Stages of normal development in the medaka Oryzias latipes. Mech Dev. 2004, 121: 605-618. 10.1016/j.mod.2004.03.012.PubMedView ArticleGoogle Scholar
- UCSC Genome Bioinformatics. [http://genome.ucsc.edu]
- Eukaryotic Genomics. [http://genome.jgi-psf.org]
- VISTA. [http://genome.lbl.gov/vista/index.shtml]
- MultAlin. [http://bioinfo.genopole-toulouse.prd.fr/multalin/multalin.html]
- Retief JD: Phylogenetic analysis using PHYLIP. Methods Mol Biol. 2000, 132: 243-258.PubMedGoogle Scholar
- Page RD: TreeView: an application to display phylogenetic trees on personal computers. Comput Appl Biosci. 1996, 12: 357-358.PubMedGoogle Scholar
- HGMP/RFCGR. [http://fugu.biology.qmul.ac.uk]
- BLAST. [http://0-www.ncbi.nlm.nih.gov.brum.beds.ac.uk/BLAST]
- Conte I, Lestingi M, den Hollander A, Alfano G, Ziviello C, Pugliese M, Circolo D, Caccioppoli C, Ciccodicola A, Banfi S: Identification and characterisation of the retinitis pigmentosa 1-like1 gene (RP1L1): a novel candidate for retinal degenerations. Eur J Hum Genet. 2003, 11: 155-162. 10.1038/sj.ejhg.5200942.PubMedView ArticleGoogle Scholar
This article is published under license to BioMed Central Ltd. This is an open access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.