Conserved and non-conserved enhancers direct tissue specific transcription in ancient germ layer specific developmental control genes
© Chatterjee et al; licensee BioMed Central Ltd. 2011
Received: 31 August 2011
Accepted: 20 October 2011
Published: 20 October 2011
Identifying DNA sequences (enhancers) that direct the precise spatial and temporal expression of developmental control genes remains a significant challenge in the annotation of vertebrate genomes. Locating these sequences, which in many cases lie at a great distance from the transcription start site, has been a major obstacle in deciphering gene regulation. Coupling of comparative genomics with functional validation to locate such regulatory elements has been a successful method in locating many such regulatory elements. But most of these studies looked either at a single gene only or the whole genome without focusing on any particular process. The pressing need is to integrate the tools of comparative genomics with knowledge of developmental biology to validate enhancers for developmental transcription factors in greater detail
Our results show that near four different genes (nkx3.2, pax9, otx1b and foxa2) in zebrafish, only 20-30% of highly conserved DNA sequences can act as developmental enhancers irrespective of the tissue the gene expresses in. We find that some genes also have multiple conserved enhancers expressing in the same tissue at the same or different time points in development. We also located non-conserved enhancers for two of the genes (pax9 and otx1b). Our modified Bacterial artificial chromosome (BACs) studies for these 4 genes revealed that many of these enhancers work in a synergistic fashion, which cannot be captured by individual DNA constructs and are not conserved at the sequence level. Our detailed biochemical and transgenic analysis revealed Foxa1 binds to the otx1b non-conserved enhancer to direct its activity in forebrain and otic vesicle of zebrafish at 24 hpf.
Our results clearly indicate that high level of functional conservation of genes is not necessarily associated with sequence conservation of its regulatory elements. Moreover certain non conserved DNA elements might have role in gene regulation. The need is to bring together multiple approaches to bear upon individual genes to decipher all its regulatory elements.
One of the paradigms of development is the regulation of the genome in a precise and synchronized manner to form a highly complex embryo with diverse and specialized cell types. Though the major cell types in the embryo contain the same genetic material, they are very different from each other in both morphology as well as function. The generation of this cellular diversity by the genome is controlled by cis-regulatory elements. Cis-regulatory elements are DNA elements that are a key component of the genome's non-coding functional sequences and consist of promoters, enhancers, silencers, insulators and locus control regions (LCR). The idea that animal development is regulated by cis-regulatory DNA elements is well established and has been elegantly described in invertebrates [1–3]. Of the many cis-regulatory elements, enhancers are critical in modulating tissue specific and time dependent gene expression during embryonic development [4–7]. Enhancers are thought to consist of clustered target sites for a number of transcription factors and collectively form the genomic instructions for developmental gene regulatory networks. Hence any approach to elucidate such networks necessitates the discovery of all constituent enhancer elements and their genomic locations.
The identification and characterization of such cis-regulatory regions within the non-coding region of vertebrate genomes remains a challenge for the post genomics era. Traditionally enhancers have been identified through deletion assays [8, 9] and in vitro footprinting [10, 11]. But recent advances in sequencing coupled with improved alignment tools have made comparative genomics one of the favored methods for enhancer detection [12–15].
It is suggested that 1% of the human non-coding genome is at least 70% conserved in the mouse genome over a region of 100 bp or longer [16, 17]. Many such conserved non-coding elements (CNEs) lie next to critical developmental control genes and have been shown to be developmental enhancers when selected and tested in mice and zebrafish [14, 18, 19]. It has also been demonstrated that genomic regions that have no sequence level conservation among different species can act as robust developmental enhancers [13, 20]. Recently chromatin immunoprecipitation followed by massive sequencing (ChIP-Seq) for P300 and ChIP followed by microarray (ChIP-Chip) on histone modifying mark H3K4me3 (Histone 3 Lysine 4 trimethylation) has been shown to be useful in identifying tissue specific enhancers on a global scale, many of which have very weak sequence conservation [21–24]. Though such studies are a valuable resource they are limited by both the amount of material generated from a particular tissue as well as the availability of a specific antibody. It is also interesting to note that all enhancers are not marked by specific histone marks or co-activator proteins.
Our study takes a "gene centric" approach to locate and validate CNEs as transcriptional enhancers, thus helping to link an enhancer to gene expression and helping to understand transcriptional control of specific genes in greater depth. We looked for CNEs to validate as enhancers around four critical developmental control genes that are expressed in tissues that are ancient in origin and conserved in evolution. We sought to investigate if genes, irrespective of which germ layer they expressed in, followed the general trend of having both sequence conserved and non-conserved enhancers or if there was a bias depending on the germ layer. We picked one gene expressing in ectoderm (otx1b), one expressing in endoderm (foxa2) and two from mesoderm (nkx3.2 and pax9). As mesoderm is considered to be of more recent origin [25, 26], we chose 2 genes expressing in this particular germ layer to assess whether there was any bias in enhancer conservation due to evolutionary time. We hypothesized that using stringent conservation criteria over large phylogenetic distances will lead to a better filtering out of the functional enhancers from non-functional CNEs.
Using modified Bacterial artificial chromosomes (BACs) to detect regulatory elements is highly beneficial when testing for multiple regulatory domains from noncontiguous DNA, which act in concert to regulate expression of a gene [27, 28]. It has also helped in cases where when the non-coding regulatory DNA is not conserved across species and thus not recognizable prior to testing. The success of this approach lies in fact that for genes with multiple expression domains, most of these domains can be observed in a single transgenic embryo as opposed to screening for multiple embryos for individual DNA constructs. Hence, we also modified the BACs containing the 4 genes of interest to locate and validate cis-regulatory elements, which are not conserved at the sequence level.
CNEs in and around nkx3.2
UCSC genome browser's alignment of sequences near Nkx3.2 in zebrafish revealed four conserved elements around the gene matching our criterion (≥ 100 bp length, ≥ 60% conservation) named Bx1 to Bx4 and varied in length from 318 to 990 bp (Figure 1A and Additional File 1). These pieces were co-injected with the minimal promoter-EGFP construct into 1-cell zebrafish embryos and observed for EGFP expression in the domains of expression of nkx3.2 in zebrafish at different time points. Only one DNA element, Bx4 drove EGFP expression that was detected along the antero-posterior axis in the region of the emerging sclerotome and myotome at 48 hpf (Figure 1D and Additional File 2). Because transient transgenics in zebrafish tend to produce mosaic expression we did immunohistochemistry with an EGFP antibody to detect the exact expression domain of this enhancer as previously described . Comparison with the RNA in situ data at the time point revealed the expression of the enhancer matched that of the endogenous gene (Additional File 3). There was no EGFP expression from the remaining CNEs at the developmental time points under observation. For this gene we also cloned all the four enhancer elements upstream of the promoter in the same order as they are found in the genome and this "mega" construct was also tested. We did not see any significant increase in expression level or domain from all the four enhancers versus only the single functional CNE (data not shown). Thus only 1 out of the 4 CNEs functioned as a developmental enhancer at all time points under observation.
A 186 kb BAC (CH73-353E16) containing the zebrafish nkx3.2gene was identified from the UCSC genome browser and modified by inserting the EGFP ORF and a selection cassette was inserted just upstream of the translational start site of the gene (Additional File 4). This ensured that the all the cis-regulatory elements in the 186 kb genomic region around Nkx3.2 that was cloned into the BAC were captured in the assay. This modified BAC was injected into 1-cell zebrafish embryos and the embryos were followed over its development to detect EGFP in domains of expression of nkx3.2. 70% of the embryos gave distinct expression along the antero-posterior axis as marked by sclerotome/myotome, at 24 and 48 hpf (Figure 1E), which correlated with the nkx3.2 expression as detected by the RNA in situ (Figure 1C, D). We were however unable to detect any EGFP expression in the branchial arch, which is also a domain of expression for the gene. This led us to conclude that the regulatory elements for this domain lay outside the 186 kb genomic DNA present in the BAC.
CNEs in and around pax9
Alignment of the sequences around pax9 revealed 11 conserved elements matching our criteria ranging in size from 311 bp to 1.2 kb (Figure 2A and Additional File 1) until there was a synteny break in one or both species that were used in the alignment. Six of these pieces were at the 5' end of the gene while five were from the 3' end, with the furthest being 178 kb away. Three out of these CNEs faithfully recaptured EGFP expression along the A-P axis as marked by the sclerotome/myotome when co-injected with the minimal promoter-reporter construct. All the three CNEs drove expression of EGFP along the presumptive sclerotome and myotome at different time points in development. Elements Px2 and Px4 which were 23 and 27 kb away from the gene on the 5' end drove expression in the sclerotome and myotome at 36 hpf (Figure 2D, E and Additional File 2), whereas Px7 which was 17 kb away from the gene at the 3' end drove EGFP expression in the same domain at 48 hpf (Figure 2F and Additional File 2). Thus 27% of the CNEs around pax9 were functional in our assay for enhancers.
Since all the three functional CNEs were expressing in the cells along the antero-posterior axis which contain both the sclerotome and myotome, we did immunohistochemistry for EGFP as described previously on the transgenic fish to ascertain the exact tissue which these enhancers were active in. Ours results show that all the three enhancers are expressing in the subset of sclerotomal cells which also expresses pax9 as seen by RNA in situ hybridization. (Additional File 5)
Branchial Arch enhancer of pax9
A 188 kb BAC (CH211-62I6) containing 188 Kb of genomic DNA flanking Pax9 was modified with EGFP and injected into 1-cell zebrafish embryo and followed over development. The modified BAC recapitulated the gene expression data in the sclerotome and myotome at 24 to 48 hpf. 70% of the embryos at 48 hpf reproduced this expression pattern in dorsal sclerotome consistently (Figure 2G).
CNEs in and around otx1b
Genome alignment revealed 16 CNEs in the syntenic block around otx1b (Figure 4A and Additional File 1). All the CNEs were present on the 5' end of the gene with the farthest being 172 kb away. There was no synteny at the 3' end of the gene with mouse and human genomes. Ox4 strongly recapitulated the gene expression at 6 hpf at the neural plate (Figure 4E and Additional File 2). Both Ox2 and Ox3 that are adjacent to each other and about 8-10 kb away from the gene drove EGFP expression in the forebrain at 48 hpf (Figure 4F, G and Additional File 2). We did not detect any EGFP expression from the other CNEs when they were assayed individually with the reporter construct. Thus 19% of the conserved elements assayed in our system were positive as tissue specific enhancers.
The 111 kb BAC (CH73-220O18) containing the Otx1 gene was modified with a reporter and was injected into the 1-cell fertilized zebrafish embryos as previously described. The EGFP expression could be detected as early as 6 hpf in the neural plate (data not shown). An individual CNE also drove reporter gene expression in the same tissue at the same time point and was present in the genomic region cloned in the BAC under testing. 69% of the injected embryos also showed EGFP expression in the forebrain midbrain and otic vesicle at 48 hpf, which overlapped with the endogenous gene expressions at the time point observed by RNA in situ (Figure 4D). We detected a much more extensive expression domain for the BAC as compared to the individual CNEs, thus leading us to believe that there exist more regulatory elements in that genomic region than detected by sequence level constrain.
CNEs in and around foxa2
Using our criteria, we located six CNEs around Foxa2, all of which were located at the 3' end of the gene, the closest being 2 kb away and the farthest being 60 kb away and ranging in size from 425 bp to 1220 bp (Figure 5A and Additional File 1). Each of these CNEs was individually ligated with the reporter construct and the injected embryos were visualized at different time points in development. Two of the CNEs, Fx2 and Fx6 could drive expression in the forebrain at 48 hpf (Figure 5G, H and Additional File 2). There was no detectable expression of EGFP from the other individual CNEs at the time points under observation. So only 33% of the CNEs assayed were functional as developmental enhancers in zebrafish.
The 179 kb BAC (CH211-38M14) containing the genomic insert in and around zebrafish foxa2 was modified and injected into 1-cell zebrafish embryos. The BAC recapitulated the endogenous gene expression particularly in the rostral end of the embryo in the forebrain and pharyngeal endoderm at 48 hpf (Figure 5E-F). We failed to detect EGFP expression in any other foxa2 expression domains.
Bac dissection to detect non-conserved enhancers
Biochemical and transgenic assays to detect a core binding site
We made "sliding" 5 bp mutations (nucleotide transversions e.g. A to C) along the 40 bp (EM4) region to narrow down the binding site (Additional File 6). Mutations in nucleotide 21 to 25 (GGAGG) completely abrogated binding of the protein with weakened binding observed for mutations in nucleotides 26-30(GAGGG) (Figure 7B). Injection of the EGFP reporter construct harboring the mutant 5 bp binding sites independently in the zebrafish showed weakened EGFP expression but not complete loss of activity (Figure 7E, F). This led us to believe that the nucleotide 21-30 together harbored the functional site in the enhancer. We went ahead and designed mutant probes for the whole 10 bp (21-30 nucleotides) region and EMSA and transgenic assays showed a complete abrogation of binding and markedly reduced EGFP expression in the zebrafish embryos (Figure 7C, G). This gave conclusive proof that the actual binding site and the core functional domain of the non-conserved enhancer lay within these 10 nucleotides.
Enhancer binding transcription factor
This study focused on examining CNEs directing expression of four developmental control genes that are ancient in origin and conserved during evolution and expressed in the three different germ layers. This allowed us to ask two very important questions regarding the nature and role of these conserved elements: (1) are all the conserved elements found in the syntenic region around a gene functional? (2) Can we find all the enhancers for a particular gene by looking only at the CNEs? The CNEs selected for our study were selected based on a criteria of ≥ 60% sequence identity and ≥ 100 bp in length across three phylogenetically distant species (human, mouse, zebrafish). As the teleost fishes and mammals last shared a common ancestor about 450 million years ago, predating the emergence of the majority of all extant vertebrates, it strongly implied that any non-coding sequences conserved between these two groups are likely to be fundamental to vertebrate life and hence functional. The design of our study has also allowed us to address some major arguments which have been put forward before to explain the absence of functionality in these CNEs in human-mouse comparative studies. By being able to look at potentially all developmental time points we had an advantage over most of the studies where the functional testing was done in mouse and was restricted to one developmental time point . Yet our data clearly indicates that even DNA elements that are conserved very stringently over long evolutionary time might not have functional roles as enhancers during development. The fact that many of these conserved non-coding elements are non-functional as enhancers and non-conserved DNA can be functional is not totally unexpected. Previous studies [9, 20] by other groups looking at single genes also revealed similar trends in other genomic loci in the zebrafish. Conversely it has also been shown that absence of sequence level homology does not preclude a fragment of DNA from being functional as a developmental enhancer [13, 43]. Thus our data supports the hypothesis that in zebrafish many of the enhancers for important developmental control genes are not constrained at the sequence level.
We have also demonstrated that there is still a need to combine high throughput comparative genomics with detailed biochemical assays to locate and validate cis-regulatory elements in the vertebrate genome. By doing such detailed assays on one of the non-conserved enhancers we have demonstrated that putative transcription factor binding sites and the actual transcription factor binding to an enhancer can be determined with great accuracy and can lead to a more complete analysis of an enhancer. Though a recent study discovered motifs in zebrafish enhancers, they only predicted what putative transcription factors could bind them , we went beyond only prediction and actually validated the responsible transcription factor. We believe our study will help in opening the door for future studies looking in finer detail at individual enhancers and how they control gene expression.
We have also shown the utility and efficiency of using BAC transgenics in zebrafish to locate and validate distal cis-regulatory elements for different genes. BACs have been successfully used to validate cis-regulatory elements for numerous genes in mouse [28, 45] and in certain cases multiple overlapping BACs have been used to locate all the regulatory elements of a gene . BAC modification for transgenics has also been employed for a few zebrafish genes in the past [46–49]. We have done an extensive spatio-temporal gene expression study of BAC transgenics for four critical developmental control genes. Thus we have demonstrated that circular modified BACs can be used without linearization for large-scale transgenic studies in zebrafish making it an attractive model system for similar future studies.
Finally it has not escaped our attention that only about 20-30% of the CNEs tested were functional as developmental enhancers. Most of these fragments could indeed be non-functional but alternately other possible explanations exist. One possibility is that these elements are acting in a synergistic manner and require other elements together to function. This was also demonstrated by our BAC studies, which show much more extensive expression patterns as compared to the individual CNEs. Also the possibility exists that these CNEs could be negative regulatory elements and hence our assay fails to capture this activity. Some recent results from other groups have indicated the same for the human genome . Also other functions such as chromatin attachment sites, miRNA genes or splice regulatory regions may also reside in such highly conserved non-coding sequences [50–54].
Thus it is becoming increasing clear that the degree of sequence conservation is not a measure of how 'important' or 'crucial' a functional element might be, but its relevance is shaped by the underlying molecular mechanisms that determine its particular function. Hence as we demonstrated here, future needs will be to combine evolutionary, functional and bioinformatic approaches to understand how these sequences function at the molecular level to determine the nature of these interactions and finally decipher what 98% of our genome encodes.
RNA in situ hybridization
DNA clones for making in situ hybridization probes were obtained from the Expressed Sequence Tag (EST) clone collection at the Genome Institute of Singapore and whole mount in situ hybridization using digoxigenin (Roche) labeled RNA probes was performed as described previously 
Fish maintenance and embryo staging
Fish were maintained in the GIS zebrafish facility according to standard procedures. Crosses were set up in the evening and the barrier was lifted the next morning. After half an hour, the fertilized embryos were collected and maintained at 28.5°C in egg water supplemented with methylene blue. The stages of embryos were indicated as hpf (hours post-fertilization) or dpf (days post-fertilization) . From 24 hours onward the embryos were raised in egg water containing Phenylthiourea. This helps to prevent pigmentation and allows better visualization of the Enhanced Green Fluorescent Protein (EGFP) signal. The GFP2 filter of stereomicroscope (Leica Microsystems, Germany) was used to observe the EGFP signals in transgenic zebrafish. The Leica DFC camera was used to take photographs and was captured using the IM50 software from Leica.
Genome alignments and Conserved Non-coding Elements (CNEs)
The University of California Santa Cruz (UCSC) genome browser's (http://genome.ucsc.edu) zebrafish (Danio rerio) Zv6 genome assembly was used as a base genome to align the genome sequences with mouse (mm8 assembly) and human (hg18 assembly) for regions around genes of interest. This assembly was used initially as the comparative genomics tracks were best annotated. Aligning the sequences with a recent assembly (Zv8) has shown that all the CNEs tested have still remained conserved and non-coding in the new annotation. For the positive CNEs the new coordinates from the Zv8 assembly are listed in Additional File 1. The closest gene to the start of each CNE was determined by the looking at known gene annotations on the UCSC genome browser. The UCSC browser's conservation tracks (phastcon conserved elements) were used to determine conserved elements (≥ 60% identity, ≥ 100 bp length). Transcribed sequences in the conserved set were filtered out using known genes, spliced ESTs and mRNA annotations (intronic conservation was allowed). We confined our search for CNEs up to 1 Mb or till wherever there was a break in synteny around the gene of interest between the three species (zebrafish-mouse-human).
We assayed for enhancer activity in embryos by either co-injecting the candidate enhancer element with a minimal promoter-reporter vector in a method previously described [19, 57] or by ligating individual CNEs to the minimal vector. Though there has been suggestions that using transposon based vectors might be a more efficient way to assay for enhancers, our previous study and also studies from other labs have shown that circular vectors containing the enhancer-reporter construct as well as co-injection of the reporter construct with the DNA are equally efficient in transient transgenic assays in zebrafish [34, 58]. In a series of control experiments we injected just the reporter construct alone (negative control), which was found to be transcriptionally inactive on its own. For preparation of DNA for microinjections, CNEs were PCR amplified from zebrafish genomic DNA or bacterial artificial chromosomes (BACs) containing the region of interest (50-100 bp sequences immediately flanking the core conserved region on either side was included in the PCR). DNA was purified using the QIAquick PCR purification kit (Qiagen, USA). The minimal vector was constructed by cloning in the EGFP ORF (Clontech, USA) under the control of a minimal promoter from the mouse β-globin gene. This fragment was released from the vector backbone using SmaI restriction enzyme and gel purified using QIAquick gel extraction kit (Qiagen, USA). Element DNA (at 20-50 pg/μl), reporter fragment (50 pg/μl) and phenol red (0.1% as tracer) were combined and injected into at least 200 1-cell zebrafish embryos using a Pico-Injector (PLI-100 Harvard Instrument, USA). Three independent injections were done for all constructs (Additional File 2).
For locating the pax9 non-conserved branchial arch enhancer, we did 500 bp overlapping PCRs spanning the 15 kb genomic region (chr17:47,911,703-47,927,438) from the final exon of pax9 to the functional CNE Px7 which is located in the third intron of the neighboring gene slc25a21. Each fragment was ligated to the reporter construct individually tested. We did not PCR out the three exons of the slc25a21. The embryos were observed at various time points in development (from 6 hpf to 3 dfp) for GFP expression. On basis of the GFP expression a CNE was scored as 1) positive, 2) negative or 3) positive but not in "expected" domain of the gene. The domains of expression were annotated by referring to the RNA in situ hybridization results for the gene and also by cross-referring to existing gene expression data at ZFIN (The Zebrafish Model Organism Database) (http://www.zfin.org).
BAC targeting construct and modification
The UCSC genome browser was used to locate 150-250 kb BACs containing the gene of interest. A BAC was selected which contained sufficient flanking genomic sequences around the gene. This enabled us to retain the regulatory elements in their proper genomic context and hence increase the possibility of capturing the maximum number of distal cis-regulatory elements. A BAC targeting vector was constructed by PCR amplifying a FRT-PGK-gb2-NEO-FRT cassette from a commercially available vector (Gene Bridges, GmbH). This was cloned into the NotI site of a pBlueScript vector SK+ (Stratagene, USA) containing the EGFP ORF between its XhoI-NotI sites. This fragment containing the EGFP ORF and selection cassette was PCR amplified with 50 bp homology arms on either side of the translational start site of each gene of interest and homologously recombined using the GeneBridges BAC modification kit (Gene Bridges, GmbH) as per the manufacturer's instructions (Additional File 4) . Circular modified BAC DNA was purified using Qiagen Maxiprep kit and injected into one-cell zebrafish embryo and EGFP expression followed over development.
The otx1b BAC was digested with XbaI restriction enzyme and the fragments generated were co-injected with the reporter vector into embryos. The fragments that drove EGFP expression were mapped back to the zebrafish genome assembly to determine its position with respect to the otx1b gene. The 6.4 kb non-conserved fragment was further subdivided into 1 kb and 200 bp fragments by generating overlapping PCRs. The 200 bp PCR generated genomic fragment that drove expression of EGFP was used in subsequent EMSAs.
Electrophoretic mobility shift assays (EMSAs)
Nuclear proteins were extracted from 24 hpf zebrafish embryos using the NE-PER nuclear and cytoplasmic extraction kit (Thermo Scientific, USA). EMSAs were carried out using DNA probes modified with 5' Cy5 labels (Sigma Proligo). The 200 bp non-conserved functional element was sub divided into six 40 bp and one 30 bp overlapping fragments having a 10 bp overlap. Equimolar amounts of complementary strands were mixed and heated to 95°C followed by gradual cooling to ambient temperature over at least 5 h to anneal the probes. For binding studies, double-stranded DNA probes at 10 nM were mixed with 10 ug of nuclear proteins and 500 ng of Poly dI-dC (Sigma) in a buffer containing 40 mM Tris-HCl (pH 8.0), 0.4 mg/ml BSA, 200 uM ZnCl2 , 400 mM KCl, 40% glycerol and 0.4% IGEPAL and incubated at 4°C in the dark for one hour. The bound and unbound probes were subsequently run on a pre-run 8% 1X TBE polyacrylamide gel for approximately 30 min at 200 V. The fluorescence was detected using a Typhoon 9140 PhosphorImager (Amersham Biosciences). At least three independent experiments for each binding site were performed to ascertain binding.
Mutation studies on the 40 bp genomic DNA fragment (EM4) that showed a shift was carried out by sequentially mutating 5 and 10 bases each along the sequence by transversion and carrying out EMSA on each mutant probe. For supershift experiments, 10 ug of 24 hpf zebrafish nuclear extract was incubated with 10 nM EM4 probe and 1 ug of antibody against the three predicted transcription factors (FOXA1: ab23738, PBX2: ab 66942, LEF1: ab 52017, all from Abcam plc), incubated for 1 hour and run on a gel as described above.
The following two morpholinos against foxa1 were injected into about 200 1 to 2-cell embryos at a concentration of 0.6-1.2 pico molar: 5'-CGCCCAACATTATGGAGGAAATCC-3' (M1) and 5'-CTTCCATTTTCACTGCGCCCAACAT-3' (M2). The specificity of foxa1 morpholinos was confirmed by using a standard scrambled control morpholino (CMO) 5'CCTCTTACCTCAGTTACAATTTATA-3' (Gene Tools).
Western blot was carried out as previously described  with some modifications. 10 morphants and wild type 24 hpf zebrafish was used to extract protein using the NE-PER nuclear and cytoplasmic extraction kit (Thermo Scientific, USA) and 25 ug of nuclear extract was loaded on the gel. Foxa1 antibody (ab23738, Abcam Plc) was applied in 1:500 dilution, followed by anti mouse IgG at 1:2500 dilution.
Total RNA was extracted from wild type and morphant 24 hpf zebrafish using TRIzol (Invitrogen, USA) as per manufacture's instruction. Reverse transcription was performed using Superscript III (Invitogen, USA) and oligo DT primers. PCR was performed on the cDNA using the primers 5'-AATCTCCATCCGTCTACATT-3' and 5'-CAGGCCGTTCATGGCGTAGG-3' for otx1b and 5'-CGGTGACATCAAGGAGCT-3' and 5'-TCGTGGATACCGCAAGATTCC-3' for b-actin.
Acknowledgements and fundings
This work is supported by the Agency for Science, Technology and Research (Singapore).
- Goto T, Macdonald P, Maniatis T: Early and late periodic patterns of even skipped expression are controlled by distinct regulatory elements that respond to different spatial cues. Cell. 1989, 57 (3): 413-422. 10.1016/0092-8674(89)90916-1.View ArticlePubMedGoogle Scholar
- Frasch M, Levine M: Complementary patterns of even-skipped and fushi tarazu expression involve their differential regulation by a common set of segmentation genes in Drosophila. Genes Dev. 1987, 1 (9): 981-995. 10.1101/gad.1.9.981.View ArticlePubMedGoogle Scholar
- Davidson EH, Rast JP, Oliveri P, Ransick A, Calestani C, Yuh CH, Minokawa T, Amore G, Hinman V, Arenas-Mena C, et al: A genomic regulatory network for development. Science. 2002, 295 (5560): 1669-1678. 10.1126/science.1069883.View ArticlePubMedGoogle Scholar
- Banerji J, Rusconi S, Schaffner W: Expression of a beta-globin gene is enhanced by remote SV40 DNA sequences. Cell. 1981, 27 (2 Pt 1): 299-308.View ArticlePubMedGoogle Scholar
- Blackwood EM, Kadonaga JT: Going the distance: a current view of enhancer action. Science. 1998, 281 (5373): 60-63.View ArticlePubMedGoogle Scholar
- Khoury G, Gruss P: Enhancer elements. Cell. 1983, 33 (2): 313-314. 10.1016/0092-8674(83)90410-5.View ArticlePubMedGoogle Scholar
- Levine M, Tjian R: Transcription regulation and animal diversity. Nature. 2003, 424 (6945): 147-151. 10.1038/nature01763.View ArticlePubMedGoogle 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 (2): 382-397. 10.1016/j.ydbio.2005.11.013.View ArticlePubMedGoogle Scholar
- Antonellis A, Huynh JL, Lee-Lin SQ, Vinton RM, Renaud G, Loftus SK, Elliot G, Wolfsberg TG, Green ED, McCallion AS, et al: Identification of neural crest and glial enhancers at the mouse Sox10 locus through transgenesis in zebrafish. PLoS Genet. 2008, 4 (9): e1000174-10.1371/journal.pgen.1000174.PubMed CentralView ArticlePubMedGoogle Scholar
- Keene MA, Corces V, Lowenhaupt K, Elgin SC: DNase I hypersensitive sites in Drosophila chromatin occur at the 5' ends of regions of transcription. Proc Natl Acad Sci USA. 1981, 78 (1): 143-146. 10.1073/pnas.78.1.143.PubMed CentralView ArticlePubMedGoogle Scholar
- McGhee JD, Wood WI, Dolan M, Engel JD, Felsenfeld G: A 200 base pair region at the 5' end of the chicken adult beta-globin gene is accessible to nuclease digestion. Cell. 1981, 27 (1 Pt 2): 45-55.View ArticlePubMedGoogle 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 (3): 212-219. 10.1016/j.ymeth.2005.12.005.View ArticlePubMedGoogle 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 (5771): 276-279. 10.1126/science.1124070.View ArticlePubMedGoogle Scholar
- Pennacchio LA, Ahituv N, Moses AM, Prabhakar S, Nobrega MA, Shoukry M, Minovitsky S, Dubchak I, Holt A, Lewis KD, et al: In vivo enhancer analysis of human conserved non-coding sequences. Nature. 2006, 444 (7118): 499-502. 10.1038/nature05295.View ArticlePubMedGoogle Scholar
- Petrykowska HM, Vockley CM, Elnitski L: Detection and characterization of silencers and enhancer-blockers in the greater CFTR locus. Genome Res. 2008, 18 (8): 1238-1246. 10.1101/gr.073817.107.PubMed CentralView ArticlePubMedGoogle Scholar
- Cooper GM, Brudno M, Stone EA, Dubchak I, Batzoglou S, Sidow A: Characterization of evolutionary rates and constraints in three Mammalian genomes. Genome Res. 2004, 14 (4): 539-548. 10.1101/gr.2034704.PubMed CentralView ArticlePubMedGoogle Scholar
- DeSilva U, Elnitski L, Idol JR, Doyle JL, Gan W, Thomas JW, Schwartz S, Dietrich NL, Beckstrom-Sternberg SM, McDowell JC, et al: Generation and comparative analysis of approximately 3.3 Mb of mouse genomic sequence orthologous to the region of human chromosome 7q11.23 implicated in Williams syndrome. Genome Res. 2002, 12 (1): 3-15. 10.1101/gr.214802.PubMed CentralView ArticlePubMedGoogle Scholar
- Sandelin A, Bailey P, Bruce S, Engstrom PG, Klos JM, Wasserman WW, Ericson J, Lenhard B: Arrays of ultraconserved non-coding regions span the loci of key developmental genes in vertebrate genomes. BMC Genomics. 2004, 5 (1): 99-10.1186/1471-2164-5-99.PubMed CentralView ArticlePubMedGoogle 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 (1): e7-10.1371/journal.pbio.0030007.PubMed CentralView ArticlePubMedGoogle Scholar
- McGaughey DM, Vinton RM, Huynh J, Al-Saif A, Beer MA, McCallion AS: Metrics of sequence constraint overlook regulatory sequences in an exhaustive analysis at phox2b. Genome Res. 2008, 18 (2): 252-260. 10.1101/gr.6929408.PubMed CentralView ArticlePubMedGoogle Scholar
- Blow MJ, McCulley DJ, Li Z, Zhang T, Akiyama JA, Holt A, Plajzer-Frick I, Shoukry M, Wright C, Chen F, et al: ChIP-Seq identification of weakly conserved heart enhancers. Nat Genet. 42 (9): 806-810.Google Scholar
- Visel A, Blow MJ, Li Z, Zhang T, Akiyama JA, Holt A, Plajzer-Frick I, Shoukry M, Wright C, Chen F, et al: ChIP-seq accurately predicts tissue-specific activity of enhancers. Nature. 2009, 457 (7231): 854-858. 10.1038/nature07730.PubMed CentralView ArticlePubMedGoogle Scholar
- Heintzman ND, Stuart RK, Hon G, Fu Y, Ching CW, Hawkins RD, Barrera LO, Van Calcar S, Qu C, Ching KA, et al: Distinct and predictive chromatin signatures of transcriptional promoters and enhancers in the human genome. Nat Genet. 2007, 39 (3): 311-318. 10.1038/ng1966.View ArticlePubMedGoogle Scholar
- Heintzman ND, Hon GC, Hawkins RD, Kheradpour P, Stark A, Harp LF, Ye Z, Lee LK, Stuart RK, Ching CW, et al: Histone modifications at human enhancers reflect global cell-type-specific gene expression. Nature. 2009, 459 (7243): 108-112. 10.1038/nature07829.PubMed CentralView ArticlePubMedGoogle Scholar
- Martindale MQ, Pang K, Finnerty JR: Investigating the origins of triploblasty: 'mesodermal' gene expression in a diploblastic animal, the sea anemone Nematostella vectensis (phylum, Cnidaria; class, Anthozoa). Development. 2004, 131 (10): 2463-2474. 10.1242/dev.01119.View ArticlePubMedGoogle Scholar
- Rodaway A, Patient R: Mesendoderm. an ancient germ layer?. Cell. 2001, 105 (2): 169-172. 10.1016/S0092-8674(01)00307-5.View ArticlePubMedGoogle Scholar
- Khandekar M, Suzuki N, Lewton J, Yamamoto M, Engel JD: Multiple, distant Gata2 enhancers specify temporally and tissue-specific patterning in the developing urogenital system. Mol Cell Biol. 2004, 24 (23): 10263-10276. 10.1128/MCB.24.23.10263-10276.2004.PubMed CentralView ArticlePubMedGoogle Scholar
- Mortlock DP, Guenther C, Kingsley DM: A general approach for identifying distant regulatory elements applied to the Gdf6 gene. Genome Res. 2003, 13 (9): 2069-2081. 10.1101/gr.1306003.PubMed CentralView ArticlePubMedGoogle Scholar
- Tribioli C, Lufkin T: The murine Bapx1 homeobox gene plays a critical role in embryonic development of the axial skeleton and spleen. Development. 1999, 126 (24): 5699-5711.PubMedGoogle Scholar
- Lettice L, Hecksher-Sorensen J, Hill R: The role of Bapx1 (Nkx3.2) in the development and evolution of the axial skeleton. J Anat. 2001, 199 (Pt 1-2): 181-187.PubMed CentralView ArticlePubMedGoogle Scholar
- Miller CT, Yelon D, Stainier DY, Kimmel CB: Two endothelin 1 effectors, hand2 and bapx1, pattern ventral pharyngeal cartilage and the jaw joint. Development. 2003, 130 (7): 1353-1365. 10.1242/dev.00339.View ArticlePubMedGoogle Scholar
- Wilson J, Tucker AS: Fgf and Bmp signals repress the expression of Bapx1 in the mandibular mesenchyme and control the position of the developing jaw joint. Dev Biol. 2004, 266 (1): 138-150. 10.1016/j.ydbio.2003.10.012.View ArticlePubMedGoogle Scholar
- Crotwell PL, Mabee PM: Gene expression patterns underlying proximal-distal skeletal segmentation in late-stage zebrafish, Danio rerio. Dev Dyn. 2007, 236 (11): 3111-3128. 10.1002/dvdy.21352.View ArticlePubMedGoogle Scholar
- Chatterjee S, Min L, Karuturi RK, Lufkin T: The role of post-transcriptional RNA processing and plasmid vector sequences on transient transgene expression in zebrafish. Transgenic Res. 19 (2): 299-304.Google Scholar
- Mise T, Iijima M, Inohaya K, Kudo A, Wada H: Function of Pax1 and Pax9 in the sclerotome of medaka fish. Genesis. 2008, 46 (4): 185-192. 10.1002/dvg.20381.View ArticlePubMedGoogle Scholar
- Hammond KL, Whitfield TT: The developing lamprey ear closely resembles the zebrafish otic vesicle: otx1 expression can account for all major patterning differences. Development. 2006, 133 (7): 1347-1357. 10.1242/dev.02306.View ArticlePubMedGoogle Scholar
- Foucher I, Mione M, Simeone A, Acampora D, Bally-Cuif L, Houart C: Differentiation of cerebellar cell identities in absence of Fgf signalling in zebrafish Otx morphants. Development. 2006, 133 (10): 1891-1900. 10.1242/dev.02352.View ArticlePubMedGoogle Scholar
- Sasaki H, Hogan BL: Differential expression of multiple fork head related genes during gastrulation and axial pattern formation in the mouse embryo. Development. 1993, 118 (1): 47-59.PubMedGoogle Scholar
- Ang SL, Wierda A, Wong D, Stevens KA, Cascio S, Rossant J, Zaret KS: The formation and maintenance of the definitive endoderm lineage in the mouse: involvement of HNF3/forkhead proteins. Development. 1993, 119 (4): 1301-1315.PubMedGoogle Scholar
- Monaghan AP, Kaestner KH, Grau E, Schutz G: Postimplantation expression patterns indicate a role for the mouse forkhead/HNF-3 alpha, beta and gamma genes in determination of the definitive endoderm, chordamesoderm and neuroectoderm. Development. 1993, 119 (3): 567-578.PubMedGoogle Scholar
- Norton WH, Mangoli M, Lele Z, Pogoda HM, Diamond B, Mercurio S, Russell C, Teraoka H, Stickney HL, Rauch GJ, et al: Monorail/Foxa2 regulates floorplate differentiation and specification of oligodendrocytes, serotonergic raphe neurones and cranial motoneurones. Development. 2005, 132 (4): 645-658. 10.1242/dev.01611.PubMed CentralView ArticlePubMedGoogle Scholar
- Wingender E, Dietze P, Karas H, Knuppel R: TRANSFAC: a database on transcription factors and their DNA binding sites. Nucleic Acids Res. 1996, 24 (1): 238-241. 10.1093/nar/24.1.238.PubMed CentralView ArticlePubMedGoogle Scholar
- Hare EE, Peterson BK, Iyer VN, Meier R, Eisen MB: Sepsid even-skipped enhancers are functionally conserved in Drosophila despite lack of sequence conservation. PLoS Genet. 2008, 4 (6): e1000106-10.1371/journal.pgen.1000106.PubMed CentralView ArticlePubMedGoogle Scholar
- Li Q, Ritter D, Yang N, Dong Z, Li H, Chuang JH, Guo S: A systematic approach to identify functional motifs within vertebrate developmental enhancers. Dev Biol. 337 (2): 484-495.Google Scholar
- Gong S, Zheng C, Doughty ML, Losos K, Didkovsky N, Schambra UB, Nowak NJ, Joyner A, Leblanc G, Hatten ME, et al: A gene expression atlas of the central nervous system based on bacterial artificial chromosomes. Nature. 2003, 425 (6961): 917-925. 10.1038/nature02033.View ArticlePubMedGoogle Scholar
- Sato Y, Miyasaka N, Yoshihara Y: Hierarchical regulation of odorant receptor gene choice and subsequent axonal projection of olfactory sensory neurons in zebrafish. J Neurosci. 2007, 27 (7): 1606-1615. 10.1523/JNEUROSCI.4218-06.2007.View ArticlePubMedGoogle Scholar
- Jessen JR, Meng A, McFarlane RJ, Paw BH, Zon LI, Smith GR, Lin S: Modification of bacterial artificial chromosomes through chi-stimulated homologous recombination and its application in zebrafish transgenesis. Proc Natl Acad Sci USA. 1998, 95 (9): 5121-5126. 10.1073/pnas.95.9.5121.PubMed CentralView ArticlePubMedGoogle Scholar
- Yang Z, Jiang H, Chachainasakul T, Gong S, Yang XW, Heintz N, Lin S: Modified bacterial artificial chromosomes for zebrafish transgenesis. Methods. 2006, 39 (3): 183-188. 10.1016/j.ymeth.2006.04.011.View ArticlePubMedGoogle Scholar
- Long Q, Meng A, Wang H, Jessen JR, Farrell MJ, Lin S: GATA-1 expression pattern can be recapitulated in living transgenic zebrafish using GFP reporter gene. Development. 1997, 124 (20): 4105-4111.PubMedGoogle Scholar
- Xie X, Lu J, Kulbokas EJ, Golub TR, Mootha V, Lindblad-Toh K, Lander ES, Kellis M: Systematic discovery of regulatory motifs in human promoters and 3' UTRs by comparison of several mammals. Nature. 2005, 434 (7031): 338-345. 10.1038/nature03441.PubMed CentralView ArticlePubMedGoogle Scholar
- Bejerano G, Pheasant M, Makunin I, Stephen S, Kent WJ, Mattick JS, Haussler D: Ultraconserved elements in the human genome. Science. 2004, 304 (5675): 1321-1325. 10.1126/science.1098119.View ArticlePubMedGoogle Scholar
- Glazko GV, Koonin EV, Rogozin IB, Shabalina SA: A significant fraction of conserved noncoding DNA in human and mouse consists of predicted matrix attachment regions. Trends Genet. 2003, 19 (3): 119-124. 10.1016/S0168-9525(03)00016-7.View ArticlePubMedGoogle Scholar
- Glazov EA, Pheasant M, Nahkuri S, Mattick JS: Evidence for control of splicing by alternative RNA secondary structures in Dipteran homothorax pre-mRNA. RNA Biol. 2006, 3 (1): 36-39. 10.4161/rna.3.1.2719.View ArticlePubMedGoogle Scholar
- Bejerano G, Lowe CB, Ahituv N, King B, Siepel A, Salama SR, Rubin EM, Kent WJ, Haussler D: A distal enhancer and an ultraconserved exon are derived from a novel retroposon. Nature. 2006, 441 (7089): 87-90. 10.1038/nature04696.View ArticlePubMedGoogle Scholar
- Qiu X, Xu H, Haddon C, Lewis J, Jiang YJ: Sequence and embryonic expression of three zebrafish fringe genes: lunatic fringe, radical fringe, and manic fringe. Dev Dyn. 2004, 231 (3): 621-630. 10.1002/dvdy.20155.View ArticlePubMedGoogle Scholar
- Kimmel CB, Ballard WW, Kimmel SR, Ullmann B, Schilling TF: Stages of embryonic development of the zebrafish. Dev Dyn. 1995, 203 (3): 253-310. 10.1002/aja.1002030302.View ArticlePubMedGoogle Scholar
- Muller F, Williams DW, Kobolak J, Gauvry L, Goldspink G, Orban L, Maclean N: Activator effect of coinjected enhancers on the muscle-specific expression of promoters in zebrafish embryos. Mol Reprod Dev. 1997, 47 (4): 404-412. 10.1002/(SICI)1098-2795(199708)47:4<404::AID-MRD6>3.0.CO;2-O.View ArticlePubMedGoogle Scholar
- Goode DK, Callaway HA, Cerda GA, Lewis KE, Elgar G: Minor change, major difference: divergent functions of highly conserved cis-regulatory elements subsequent to whole genome duplication events. Development. 138 (5): 879-884.Google Scholar
- Chatterjee S, Lufkin T: Fishing for function: zebrafish BAC transgenics for functional genomics. Mol Biosyst. 2011, 7 (8): 2345-2351. 10.1039/c1mb05116d.PubMed CentralView ArticlePubMedGoogle 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.