Cellular Expression of Smarca4 (Brg1)-regulated Genes in Zebrafish Retinas
- Monica R Hensley†1,
- Farida Emran†2,
- Sylvia Bonilla1,
- Liyun Zhang1,
- Wenxuan Zhong3,
- Paul Grosu4,
- John E Dowling2 and
- Yuk Fai Leung1, 2Email author
© Hensley et al; licensee BioMed Central Ltd. 2011
Received: 24 February 2011
Accepted: 14 July 2011
Published: 14 July 2011
In a recent genomic study, Leung et al. used a factorial microarray analysis to identify Smarca4 (Brg1)-regulated genes in micro-dissected zebrafish retinas. Two hundred and fifty nine genes were grouped in three-way ANOVA models which carried the most specific retinal change. To validate the microarray results and to elucidate cellular expression patterns of the significant genes for further characterization, 32 known genes were randomly selected from this group. In situ hybridization of these genes was performed on the same types of samples (wild-type (WT) and smarca4 a50/a50 (yng) mutant) at the same stages (36 and 52 hours post-fertilization (hpf)) as in the microarray study.
Thirty out of 32 riboprobes showed a positive in situ staining signal. Twenty seven out of these 30 genes were originally further classified as Smarca4-regulated retinal genes, while the remaining three as retinal-specific expression independent of Smarca4 regulation. It was found that 90.32% of the significant microarray comparisons that were used to identify Smarca4-regulated retinal genes had a corresponding qualitative expression change in the in situ hybridization comparisons. This is highly concordant with the theoretical true discovery rate of 95%. Hierarchical clustering was used to investigate the similarity of the cellular expression patterns of 25 out of the 27 Smarca4-regulated retinal genes that had a sufficiently high expression signal for an unambiguous identification of retinal expression domains. Three broad groups of expression pattern were identified; including 1) photoreceptor layer/outer nuclear layer specific expression at 52 hpf, 2) ganglion cell layer (GCL) and/or inner nuclear layer (INL) specific expression at both 36 & 52 hpf, and 3) GCL and/or INL specific expression at 52 hpf only. Some of these genes have recently been demonstrated to play key roles in retinal cell-type specification, differentiation and lamination. For the remaining three retinal-specific genes that are independent of Smarca4 regulation, they all had a subtle expression difference between WT and smarca4 a50/a50 retinas as detected by in situ hybridization. This subtle expression difference was also detected by the original microarray analysis. However, the difference was lower than the fold change cut-off used in that study and hence these genes were not inferred as Smarca4-regulated retinal genes.
This study has successfully investigated the expression pattern of 32 genes identified from the original factorial microarray analysis. The results have demonstrated that the true discovery rate for identifying Smarca4-regulated retinal genes is 90.3%. Hence, the significant genes from the microarray study are good candidates for cell-type specific markers and will aid further investigation of retinal differentiation.
The retina contains six neuronal cell types and one major glial cell type, all of which originated from the same progenitor cell. In vertebrate retinas, the early proliferation of multipotent progenitors produces a sufficient number of cells that will be specified as different retinal cell types at around the time of their cell cycle withdrawal. The cell types that can be specified at a particular stage of development are controlled by both intrinsic and extrinsic signals; and this temporal restriction on cell fate specification is commonly referred to as the competence model . In general, ganglion cells (GCs) are the first cell type while Muller cells (MCs) are the last cell type to be specified. These reversibly specified cells then undergo irreversible determination after which their cell fate can not be changed by external signals . Specification and determination are collectively called commitment. Finally, these committed cells undergo terminal differentiation during which they send out neuronal processes and synapse with each other.
Many intrinsic factors and extrinsic signals that control retinal cell specification and differentiation have been identified. For example, Atoh7 and Crx are transcription factors that specify retinal ganglion cells  and photoreceptors [4, 5] respectively. Research on zebrafish retina has also identified key signalling molecules and processes that regulate terminal differentiation and lamination. These include Shh , cell polarity regulation , cell adhesion  and chromatin remodelling . For example, in a retinal terminal differentiation mutant smarca4 a50/a50 (young/yng), all retinal cell types could be specified but failed to fully differentiate due to a null mutation of a chromatin remodelling component smarca4 (brg1) [9, 10]. Nonetheless, the underlying genetic circuitry that controls terminal differentiation is still poorly understood .
In a recent microarray study of retinal development in zebrafish, a factorial design was used to identify specifically Smarca4-regulated retinal genes in micro-dissected wild-type (WT) and smarca4 a50/a50 retinas . Seven hundred and thirty one genes were identified as regulated by Smarca4 with different levels of retinal specificity, as defined by the factorial ANOVA model used. In particular, 259 genes were grouped in the three-way ANOVA models and considered to have the most specific retinal change. To study the regulation of these genes during retinal development, it is essential to validate their differential expression and elucidate their cellular expression patterns. The present study investigated the cellular expression patterns of 32 randomly selected known genes from this group of three-way ANOVA models by whole-mount in situ hybridization. A comparative analysis of the results with the original microarray findings was performed to evaluate the performance of the factorial microarray analysis.
In situhybridization analysis of Smarca4-regulated genes validates differential expression results obtained by factorial microarray analysis
In a recent factorial microarray analysis that aimed at identifying Smarca4-regulated retinal genes in micro-dissected retinas , 259 genes were categorized in three-way ANOVA models (see Methods for details). These genes have the most specific gene expression change in the smarca4 a50/a50 retinas at 52 hours post-fertilization (hpf). Since the differentiating retinal cells are beginning to organize themselves into three cellular layers and synapse with each other at this stage, the cellular expression analysis of these Smarca4-regulated retinal genes would potentially assist further functional characterization of retinal differentiation and lamination.
To facilitate performance evaluation of this microarray analysis approach, 32 known genes were randomly selected from these 259 genes and whole-mount in situ hybridization performed on embryos collected at the same stages that were used in the microarray analysis. These include WT and smarca4 a50/a50 embryos at 36 & 52 hpf. The corresponding retinal samples at the same stages are WR36 & WR52 for WT, and YR36 & YR52 for smarca4 a50/a50 respectively (former nomenclature yng will be used in the sample abbreviations for smarca4 a50/a50 to aid comparison with the published data in ). The factorial microarray analysis also further categorized these 259 genes into three functional groups: (i) Smarca4-regulated retinal differentiation genes (194 genes), (ii) retinal-specific genes independent of Smarca4 regulation (35 genes), and (iii) Smarca4-regulated genes outside retina (54 genes). Only genes in the first group have a significant differential retinal gene expression in smarca4 a50/a50 compared to WT (see Methods for the statistical criteria).
Among the 32 randomly selected known genes for in situ hybridization, 29 (olfm2, rlbp1l, guk1, ndrg1, glra4b, robo2, irx7, barhl2, vangl1, dtnbp1, lmo4l, cdh11, elovl4, ctbp2, irx4a, rho, id2a, sv2b, wnt11, calb2l, ckmt1, kal1a, aanat2, pbx1a, fzd8b, tfap2a, nme2l, rcv1 and gnat1) are in the functional group (i) as described above, three (foxn4, six3a and bhlhe22) are in group (ii) and none is in group (iii). Two (cdh11 and wnt11) out of the 29 genes from group (i) showed no discernable signal in all retinal samples in the in situ hybridization experiments. For cdh11, there was considerable expression signal detected outside the eye region (data not shown). This indicates that the probe could detect gene expression. However, there was no discernible expression signal in the eye. For wnt11, since there was no signal in the whole embryo (data not shown), an additional probe that is deposited in the public domain (ZFIN ID: ZDB-GENE-990603-12; probe cb748) was acquired and the in situ hybridization repeated with smarca4 a50/a50 and their WT siblings, as well as WT embryos collected from AB and TL strains. In all conditions, the staining in the other reported embryonic regions, including otic vesicles, lower jaw and myoseptum, was clearly observed. However, there was no positive staining in the retinas, even after an overnight staining (Additional file 1, Figure S1 and data not shown). Thus, the in situ hybridization results do not support the expression of wnt11 and cdh11 in the retinas at 36 and 52 hpf.
Cellular expression patterns of Smarca4-regulated retinal genes
I. Photoreceptor layer/ONL specific expression at 52 hpf
A total of eight genes, including rho, gnat1, guk1, nme2l, elovl4, rcv1, ndrg1 and aanat2, were found to express specifically in photoreceptor layer/ONL in the WT retinas (Figure 4). Their expression was suppressed in all smarca4 a50/a50 retinas. Two sub-clusters were also observed (Figure 4G). The genes in the first sub-cluster only expressed in a ventral patch of cells in the retina by 52 hpf (rho & gnat1, Figure 4A; nme2l & elovl4, Figure 4C). The genes in the second sub-cluster expressed to a different degree in the other parts of the ONL at the same stage in addition to the ventral patch of the retina (guk1, Figure 4B; rcv1 and aanat2, Figure 4D & 4F respectively), except for ndrg1 which did not express in this domain (Figure 4E). Also, the expression patterns of rho and gnat1, as well as nme2l and elovl4 were highly similar and were tightly clustered with each other (p-value <0.05; Figure 2D, red colour in the dendrogram). Interestingly, rho, gnat1 and rcv1, components of the visual cycle, were also weakly expressed in the ventral retina at 36 hpf (Figure 4A & 4D). This is earlier than the stage (~50-52 hpf) of the first rho expression previously reported , but is consistent with the public in situ hybridization data (ZFIN ID: ZDB-GENE-990415-271), which show the first expression of rho in the same location at around the same stage (Prim-15 - 25; ~30-36 hpf).
II. GCL and/or INL specific expression at both 36 & 52 hpf
A group of five genes, including irx4a, barhl2, pbx1a, rlbp1l and irx7, were expressed in WT retinas at both 36 and 52 hpf (Figure 2D). Their expressions were all suppressed in smarca4 a50/a50 retinas at the corresponding stages. These genes were segregated into two sub-clusters because of their diverse expression patterns at 36 hpf (Figure 5F). The first sub-cluster contains irx4a, barhl2 and pbx1a (Figure 5A, B & 5C). This group of genes expressed in the GC region at 36 hpf to a variable degree (Figure 5, arrows). For example, barhl2 was expressed in the VP, AV, AD and PD but not in the PV domain from the ventral view (Figure 5B, arrows), suggesting that it is a gene that may be regulated by the neurogenic wave in the retina . By 52 hpf, all of these three genes were expressing in either the GCL or both GCL and INL-b, which contains the presumptive ACs. The second sub-cluster, containing rlbp1l and irx7 (Figure 5D & 5E) , is clustered within a larger family of genes that only expressed at 52 hpf and will be discussed in the next section. Nonetheless, these two genes have a slightly unique expression pattern: whereas they were expressed on the ventral side of the 36 hpf retinas (Figure 5D & 5E, arrows), their expression domain was restricted to the whole INL by 52 hpf.
III. GCL and/or INL specific expression at 52 hpf only
A total of 14 genes (irx7, rlbp1l, sv2b, ctbp2, kal1a, ckmt1, robo2, vangl1, glra4b, calb2l, lmo4l, olfm2, id2a and tfap2a) were clustered in a group with a p-value < 0.1 (Figure 2D, green colour part of the dendrogram). Among them, only irx7 and rlbp1l started to express at 36 hpf and were discussed above (Figure 5D & 5E). For the remaining 12 genes in the cluster, they only expressed at 52 hpf (Figure 3). Interestingly, several genes inside this cluster have a highly similar and unique expression pattern in that they were tightly grouped in two sub-clusters with a p-value< 0.05 (Figure 2D, red colour part of the dendrogram). The first sub-cluster consists of four genes - tfap2a, id2a, olfm2 and lmo4l (Figure 3A) that exclusively expressed in the INL-b, and three genes - calb2l, glra4b and vangl1 (Figure 3B) that expressed in both GCL and INL-b. The second sub-cluster consists of two genes - robo2 (Figure 3C) and ckmt1 (Figure 3D) that both expressed strongly in GCL. Also, ckmt1 had an addition expression domain in the INL-a. There were three genes - kal1a (Figure 3E), ctbp2 (Figure 3F) and sv2b (Figure 3G) that did not belong to these two sub-clusters and formed a separate one. All of them had broad expression domains in INL that at least spanned the INL-m and INL-a. Sv2b further expressed in the GCL as well. Interestingly, kal1a was strongly expressed in a subset of cells in the AC region that the neuronal processes were intensely stained and highlighted the inner plexiform layer (IPL) (Figure 3E, arrow).
IV. Fzd8b only expressed in smarca4 a50/a50 retinas at 36 hpf
Cellular expression patterns of retinal-specific genes that were not inferred to be regulated by Smarca4
This study has successfully investigated the expression of 32 genes identified from a microarray study of Smarca4-regulated retinal genes from micro-dissected zebrafish retinas and elucidated the cellular expression patterns. The original microarray analysis has inferred 29 of these 32 genes as Smarca4-regulated retinal genes and the remaining three as retinal specific genes whose expression is independent of Smarca4 regulation. The in situ hybridization analysis has demonstrated that the sign of the expression changes matches with the microarray results well (Figure 1). In particular, 90.3% (28/31 cases) of the significant microarray expression comparisons have a corresponding differential expression between WT and smarca4 a50/a50 retinas in the in situ hybridization experiments. This suggests strongly that factorial microarray analysis, can efficiently detect true biological differences among different types of micro-dissected zebrafish retinas. Twenty seven out of the 29 Smarca4-regulated retinal genes have a positive in situ hybridization signal in the retinas, and 25 of them have the cellular expression domains that could be clearly discerned. Three general groups of expression patterns can be identified from these 25 genes (Figure 2D): 1) photoreceptor layer/ONL specific expression at 52 hpf, 2) GCL and/or INL specific expression at both 36 & 52 hpf and 3) GCL and/or INL specific expression at 52 hpf only.
The acquisition of the cellular expression patterns and the clustering of these genes can facilitate downstream characterization, because genes that are expressed in the same domain are likely to be regulated by the same upstream factor(s) and/or are used to control development and functions of specific cell type(s). For example, a group of photoreceptor/ONL genes were expressed on the ventral side of WT retinas at 52 hpf and were suppressed in the smarca4 a50/a50 retinas (Figure 4). This is the beginning stage of photoreceptor generation in the ventronasal patch of the retina , and many of these genes are related to the visual cycle. Interestingly, rcv1 and aanat2 are expressed in virtually the whole ONL at 52 hpf in addition to the ventronasal patch (Figure 4D & 4F), while ndrg1 is expressed in the whole ONL but not the ventronasal patch at this stage (Figure 4E). This extensive expression pattern is obviously different from the previously reported expression patterns of cones and rods differentiation, as indicated by the wave-like expression pattern of cone opsin and sporadic expression of rod opsin respectively starting from the ventronasal retina at around 52 hpf . Since rcv1 also participates in the visual cycle , the extensive expression of rcv1, aanat2 and ndrg1 in the ONL at 52 hpf suggests that there may be a different regulatory circuit that controls their expression. Together with the detection of rho, gnat1 and rcv1 in the ventral retina at 36 hpf (Figure 4A & 4D), these observations also suggest that the differentiation program of photoreceptors may begin earlier than 50-52 hpf. In the second example, a few genes including tfap2a, id2a, olfm2 and lmo4l (Figure 3A) , are expressed specifically on the basal side of the INL which is the presumptive region of ACs. Their expression domains were highly overlapped with each other that these genes were grouped together in a highly significant sub-cluster (p-value < 0.05; Figure 2D). This suggests that they may play a role in ACs development. Indeed, a recent perturbation experiment of id2a in zebrafish  have demonstrated such a role. In addition, the knockout of an orthologue of lmo4l, LMO4, in mice leads to a reduction of amacrine cells . In the third example, barhl2 is a gene that has been demonstrated to regulate the development of GCs and ACs in mice . Its expression pattern in the zebrafish correlates nicely with the development of GCs and ACs (Figure 5B), suggesting that barhl2 is likely to play a similar role in zebrafish retinal development.
One key objective of the original microarray study was to identify Smarca4-regulated retinal genes because their expression should either mediate terminal differentiation and/or be a consequence of terminal differentiation. The hallmarks of terminal differentiation of neurons include neurite outgrowth and synaptogenesis, which would be manifested as the formation of the retinal lamination. The current in situ hybridization analysis of the candidate genes, which are mostly differentially expressed at 52 hpf, during which the retinal lamination is first observed, may potentially provide new insights into the molecular control of this process. There are at least two non-exclusive possibilities: 1) these genes control the differentiation of a specific cell type which in turn orchestrates the formation of the final neuronal connection and lamination and/or 2) these genes control patterning and lamination of the retina directly. In the microarray study, it was speculated that tfap2, id2 and irx transcription factor families would play a key role in retinal differentiation and lamination , and this has been proven by studies of tfap2a knockout mice , and id2a  and irx1a [22, 23] knockdowns in zebrafish. In fact, id2a has also been demonstrated to play a role in the development of ACs as discussed above and the knockdown of tfap2c in tfap2a zebrafish mutant shows a very severe eye development phenotype, in which the eye size becomes significantly smaller . These results together suggest that these transcription factor families are integral components of retinal differentiation and lamination, and that the remaining members of these groups may play key roles in retinal development. Indeed, our investigation has demonstrated that irx7, a gene that has an expression pattern that spans the whole INL at 52 hpf, plays a role in the differentiation of both INL and photoreceptors, as well as lamination ( & manuscript in prep). Also, in the original microarray study, Cdk5 and activators were demonstrated to play a role in retinal differentiation . The execution of their functions likely involves a subset of the genes presented in this study. Thus, the knockdown of these key mediators of Smarca4-network and subsequent in situ hybridization analysis of the other Smarca4-regulated genes, which is currently in progress, will likely further the understanding of retinal differentiation.
While the microarray analysis accurately excluded genes that did not meet all criteria to be Smarca4-regulated retinal genes, some of these excluded genes might have a role in subtle developmental change in retinas that are regulated by Smarca4. For instance, the fold change of a key transcription factor in a network can be less than two, and yet it can play a substantial role in the developmental process that it controls. Foxn4 was excluded by the original microarray analysis as being a Smarca4-regulated retinal gene because its fold change was below the two-fold cut-off despite having a significant corresponding contrast. Nonetheless, it was obvious that there was a discernible difference in the expression level and cellular expression pattern of this gene between WR52 and YR52 (Figure 7A). The expression seemed to be relatively higher on the dorsal side of the WR52; while in YR52, the expression was generally higher but there was no such regional difference. In addition, the expression of foxn4 at 52 hpf is restricted to the proliferative marginal zone region; while at 36 hpf, it is primarily located in the proliferative part of the retinal neuro-epithelium (Figure 7A and data not shown). These results suggest that foxn4 may also play a role in the Smarca4-regulated retinal differentiation through the regulation of the retinal progenitors. Indeed, foxn4 has been demonstrated to express in retinal progenitor cells in mice and the null mutant has retinal dysplasia . While a small and consistent fold change of a gene can be biologically significant, it is necessary to implement a cut-off to ensure the false discovery rate is low. This is a necessary trade-off to maximize the number of true positive genes, especially in the first genomic study of micro-dissected zebrafish retinas. It is obvious that while the original factorial microarray analysis have efficiently identified true Smarca4-regulated retinal expression, it is possible that some genes that are subtly regulated by Smarca4 did not meet the statistical cut-off and were excluded from this category.
This study has investigated the expression of 32 genes identified from a factorial microarray analysis that was designed for searching Smarca4-regulated retinal genes. Twenty nine of them were grouped in this category while three were excluded by the original analysis. The in situ hybridization results have strongly supported that the original factorial microarray analysis can efficiently identify true Smarca4-regulated retinal genes with a true positive rate of 90.3%. Thus, this gene list is not only a rich source for cell-type specific markers, but also may assist in formulating new hypothesis to study the underlying genetic regulatory circuits that control retinal development.
Zebrafish (AB & TL wild-type (WT) and smarca4 a50/+ (yng)) were maintained according to standard procedures . All protocols were approved by the Harvard Standing Committee on the Use of Animals in Research and Teaching and Purdue Animal Care and Use Committee.
Egg collection, embryo staging and collection
Wild-type (WT) and smarca4 a50/a50 (yng) embryos were collected at the same stages (i.e. 36 and 52 postfertilization (hpf)) as in the microarray analysis . To ensure all embryos were at a similar stage during collection, embryos were collected every 20 minutes in E3 medium . Zebrafish embryo staging was performed as described , with the following modifications. First, embryos from each 20-minute collection were inspected at around 10-12 hpf, and any abnormal looking embryos were discarded. Then a sample of embryos, typically 5-7 for each 20-minute collection, was staged. This would be considered the average stage for all the embryos collected during the same period. The staged embryos were treated with 0.003% phenylthiourea (PTU) (Sigma; St. Louis, MO) in E3 medium starting at 23 hpf to prevent melanization. At least 10 embryos were set aside and kept in E3 as staging reference. The embryos treated with PTU were staged again at 36 and 52 hpf immediately before collection, using the staging reference embryos that were not treated with PTU. Standard staging criteria were used in conjunction with the eye size parameter for these stages . Smarca4 a50/a50 embryos were staged by using their WT siblings collected from the same clutch at the same time. Embryos were dechorionated by treating with 1 mg/mL Pronase (Sigma; St. Louis, MO) for 1 to 2 minutes until embryos started to come out from the chorion; then all embryos were rinsed extensively with E3 medium. All embryos were immediately fixed in 4% paraformaldhyde (PFA)(Sigma; St. Louis, MO) in phosphate buffered saline (PBS), dehydrated and stored in 100% methanol at -20°C as described .
Complementary DNA library preparation
Total RNA was extracted and purified from 2-day old zebrafish embryos as described . Complementary DNA (cDNA) library was prepared by reverse transcribing messenger RNAs from purified total RNA using Superscript II reverse transcriptase (Invitrogen, Carlsbad, CA) and an anchored primer 5'-dT20VN-3' (Integrated DNA Technologies, Coralville, IA).
Gene selection for in situhybridization analysis
indicating that the expression of the gene was not affected by any of the three biological changes.
Genes that were categorized in the three-way model group had the most specific regulation in the smarca4 a50/a50 retinas, because the three-way interaction term (T*M*R) indicates the presence of a mutation-specific effect on gene expression in the retina at 52 hpf. Since this is the stage when retinal lamination is first formed, the first cellular expression investigation was focused on this group of genes. Thirty two known genes were randomly selected from 259 three-way models for further in situ hybridization analysis. Also, in the same microarray analysis study, these 259 models were further categorized into three functional groups: (i) Smarca4-regulated retinal genes, (ii) Retina-specific genes independent of Smarca4 regulation, and (iii) Smarca4-regulated genes outside the retina. Note that a significant 3-way interaction term does not necessarily translate to a final overall significance, because that was inferred by two criteria : (1) a significant contrast (false discovery rate p-value < 0.05), which is the specific comparison of two biological conditions using the corresponding coefficients from the ANOVA model and (2) fold change of the corresponding comparison > = 2. A comparison can be insignificant because the T*M*R coefficient can cancel out the effect on expression imposed by other terms, even though the T*M*R coefficient is significant by itself. To infer significant Smarca4-regulated retinal genes for the functional group (i), the following contrasts were used:
while for YR52/WR52:
The goal of the primers design was to select primers that specifically flanked 500-850 bps of the candidate genes so that the resulting amplicon will overlap with the target sequence of the probesets of the Affymetrix zebrafish whole genome array (Affymetrix, Santa Clara, CA) as much as possible for riboprobe synthesis. Since this GeneChip was used in the microarray study , designing a riboprobe with maximal sequence overlap with the probeset sequence also allows for better expression validation of the original microarray results. The primers were designed using the following procedures until a pair of primers was chosen for that particular gene. First, a computer script (available upon request) was written to analyze and locate the consensus sequence of a probeset from the Affymetrix design file. Primers that flanked 500-850 bps including the target region were then selected by Primer3 , using default parameters except the option 'Max Poly-X' was set to 3. However, this approach was not always feasible, as Affymetrix often selected the target region for their probesets from the 3' end or the untranslated region (UTR) of a gene. The low complexity of the 3'UTR sometimes precluded an efficient primer design and/or PCR amplification. In this case the target region of the Primer3 would then be moved 5' in 50 bps increments until a pair of primers was picked. Sometime the consensus sequence was too short, i.e. shorter than 500 bps. In that case, the ReqSeq sequence would be obtained by the respective Genbank ID. The longer sequence would then be used for primer design as described above. The specificity of the selected primers was checked by blasting NCBI's zebrafish nr and htgs databases. All primers selected by this exercise are shown in Additional file 5, Table S4.
Whole-mount in situhybridization
Specific DNA fragments were amplified from a 2-day old zebrafish cDNA library and cloned into pGEM-T easy vector (Promega, Madison, WI) for probe synthesis. All constructs were sequenced verified. Riboprobe synthesis and whole-mount in situ hybridization were performed using procedures described . The stringency washes after probe hybridization and before signal detection were performed using the Biolane semi-automated in situ hybridization machine (INTAVIS Bioanalytical Instruments AG, Koeln, Germany). The specific parameters of the washing steps for the embryos are listed as follows:
Twice in 50% formamide/2X SSCT (Saline-Sodium Citrate with 0.1% Tween-20) for 20 minutes each at 65°C
Once in 2X SSCT for 15 minutes at 65°C
Twice in 0.2X SSCT for 20 minutes each at 65°C
Twice in PBST (PBS with 0.1% Tween-20) for five minutes each at room temperature
Once in Block buffer (2 mg/mL Bovine Serum Albumin, 2% normal sheep serum in PBST) for two hours at room temperature
Once in 1:3000 Anti-dioxgenin-alkaline phosphatase antibody (Roche Applied Science; Indianapolis, IN) in block for 12 hours at 4°C
◦ The antibody was pre-adsorpted at 1:1000 dilution in Block buffer with at least 400 fixed embryos for an overnight at 4°C before used
◦ The pre-adsoprtion step slightly increases signal clarity, but not very substantial, so this can be considered as an optional preparation
◦ The final working stock can be re-used for a few times and should be stored at 4°C
Six times in PBST at room temperature; the first one for five minutes, the next four for 30 minutes each and the final one is programmed for 12.5 hours, which can be stopped at any time after the first 30 minutes.
Three times in staining buffer (100 mM Tris pH9.5, 50mM MgCl2, 100 mM NaCl, 0.1% Tween-20 and 1 mM levamisol (Sigma; St. Louis, MO)) for five minutes each at room temperature
At least 20 embryos were analyzed for each antisense probe for each developmental stage and genotype. The same number of WT embryos at the same developmental stage was also collected and used for the sense probe control. All embryos used for the characterization of the same gene were stained for the same period of time to maximize comparability between conditions. The samples were destained by a 2:1 mixture of benzyl benzoate and benzyl alcohol and stored in 4% PFA at 4°C.
Embryos were mounted in 3% methylcellulose on a depression slide for observation and imaging. All images were taken by a Spot RT3 slider CCD camera (Diagnostic Instruments, Sterling Heights, MI) mounted on an Olympus SZX16 Stereomicroscope (Olympus, Center Valley, PA). The images were subsequently merged using Adobe Photoshop CS3 (Adobe Systems Incorporated, San Jose, CA).
Hierarchical clustering analysis of the expression domains of the genes in different samples was conducted as follows: First, the expression domain of the genes, as defined in Figure 2A, B & 2C, was visually scored as 0 (absent) or 1 (present) (Additional file 3, Table S2). Then the dissimilarity of these expression domains between each gene was calculated by a binary distance measure and hierarchical clustering conducted using average group linkage agglomeration. Finally, p-values of the clustering results were calculated by a multi-scale bootstrap resampling method . The number of bootstrap replications was 10000. All data analyses were conducted in the R statistical environment version 2.11.1 http://www.r-project.org. Heatmaps were generated using Multiexperiment Viewer (MeV) http://www.tm4.org.
We thank Doug Melton for granting us the access to a Biolane machine and Joe Qiao Zhou for his advice on optimizing the in situ hybridization procedures for the machine processing at the beginning stage of the project. This work was supported by grants from the Knights Templar Eye Foundation to LZ and YFL, a grant from the Office of Science (BER), U.S. Department of Energy and an NSF grant DMS 1120256 to WZ. NIH/NEI grants R01 EY000811 to JED and F32 EY018044 to FE, and awards from Hope for Vision & Showalter Research Trust to YFL.
- Livesey FJ, Cepko CL: Vertebrate neural cell-fate determination: lessons from the retina. Nat Rev Neurosci. 2001, 2: 109-18. 10.1038/35053522.View ArticlePubMedGoogle Scholar
- Gilbert SF: Developmental biology. 2006, Sunderland, Mass.: Sinauer Associates, Inc. Publishers, 8Google Scholar
- Kay JN, Finger-Baier KC, Roeser T, Staub W, Baier H: Retinal ganglion cell genesis requires lakritz, a Zebrafish atonal Homolog. Neuron. 2001, 30: 725-36. 10.1016/S0896-6273(01)00312-9.View ArticlePubMedGoogle Scholar
- Chen S, Wang QL, Nie Z, Sun H, Lennon G, Copeland NG, Gilbert DJ, Jenkins NA, Zack DJ: Crx, a novel Otx-like paired-homeodomain protein, binds to and transactivates photoreceptor cell-specific genes. Neuron. 1997, 19: 1017-30. 10.1016/S0896-6273(00)80394-3.View ArticlePubMedGoogle Scholar
- Furukawa T, Morrow EM, Cepko CL: Crx, a novel otx-like homeobox gene, shows photoreceptor-specific expression and regulates photoreceptor differentiation. Cell. 1997, 91: 531-41. 10.1016/S0092-8674(00)80439-0.View ArticlePubMedGoogle Scholar
- Neumann CJ, Nuesslein-Volhard C: Patterning of the zebrafish retina by a wave of sonic hedgehog activity. Science. 2000, 289: 2137-9. 10.1126/science.289.5487.2137.View ArticlePubMedGoogle Scholar
- Wei X, Malicki J: nagie oko, encoding a MAGUK-family protein, is essential for cellular patterning of the retina. Nat Genet. 2002, 31: 150-7. 10.1038/ng883.View ArticlePubMedGoogle Scholar
- Masai I, Lele Z, Yamaguchi M, Komori A, Nakata A, Nishiwaki Y, Wada H, Tanaka H, Nojima Y, Hammerschmidt M, et al: N-cadherin mediates retinal lamination, maintenance of forebrain compartments and patterning of retinal neurites. Development. 2003, 130: 2479-94. 10.1242/dev.00465.View ArticlePubMedGoogle Scholar
- Gregg RG, Willer GB, Fadool JM, Dowling JE, Link BA: Positional cloning of the young mutation identifies an essential role for the Brahma chromatin remodeling complex in mediating retinal cell differentiation. Proc Natl Acad Sci USA. 2003, 100: 6535-40. 10.1073/pnas.0631813100.PubMed CentralView ArticlePubMedGoogle Scholar
- Link BA, Fadool JM, Malicki J, Dowling JE: The zebrafish young mutation acts non-cell-autonomously to uncouple differentiation from specification for all retinal cells. Development. 2000, 127: 2177-88.PubMedGoogle Scholar
- Johnson KG, Shewan D, Holt CE: Visual System Development in Vertebrates. Book Visual System Development in Vertebrates. 2005, City: John Wiley & Sons, LtdGoogle Scholar
- Leung YF, Ma P, Link BA, Dowling JE: Factorial microarray analysis of zebrafish retinal development. Proc Natl Acad Sci USA. 2008, 105: 12909-14. 10.1073/pnas.0806038105.PubMed CentralView ArticlePubMedGoogle Scholar
- Hu M, Easter SS: Retinal neurogenesis: the formation of the initial central patch of postmitotic cells. Dev Biol. 1999, 207: 309-21. 10.1006/dbio.1998.9031.View ArticlePubMedGoogle Scholar
- Schmitt EA, Dowling JE: Comparison of topographical patterns of ganglion and photoreceptor cell differentiation in the retina of the zebrafish, Danio rerio. J Comp Neurol. 1996, 371: 222-34. 10.1002/(SICI)1096-9861(19960722)371:2<222::AID-CNE3>3.0.CO;2-4.View ArticlePubMedGoogle Scholar
- Schmitt EA, Dowling JE: Early retinal development in the zebrafish, Danio rerio: light and electron microscopic analyses. J Comp Neurol. 1999, 404: 515-36. 10.1002/(SICI)1096-9861(19990222)404:4<515::AID-CNE8>3.0.CO;2-A.View ArticlePubMedGoogle Scholar
- Raymond PA, Barthel LK, Curran GA: Developmental patterning of rod and cone photoreceptors in embryonic zebrafish. J Comp Neurol. 1995, 359: 537-50. 10.1002/cne.903590403.View ArticlePubMedGoogle Scholar
- Stryer L: Visual excitation and recovery. J Biol Chem. 1991, 266: 10711-4.PubMedGoogle Scholar
- Uribe RA, Gross JM: Id2a influences neuron and glia formation in the zebrafish retina by modulating retinoblast cell cycle kinetics. Development. 2010, 137: 3763-74. 10.1242/dev.050484.PubMed CentralView ArticlePubMedGoogle Scholar
- Duquette PM, Zhou X, Yap NL, MacLaren EJ, Lu JJ, Wallace VA, Chen HH: Loss of LMO4 in the retina leads to reduction of GABAergic amacrine cells and functional deficits. PLoS One. 2010, 5: e13232-10.1371/journal.pone.0013232.PubMed CentralView ArticlePubMedGoogle Scholar
- Ding Q, Chen H, Xie X, Libby RT, Tian N, Gan L: BARHL2 differentially regulates the development of retinal amacrine and ganglion neurons. J Neurosci. 2009, 29: 3992-4003. 10.1523/JNEUROSCI.5237-08.2009.PubMed CentralView ArticlePubMedGoogle Scholar
- Bassett EA, Williams T, Zacharias AL, Gage PJ, Fuhrmann S, West-Mays JA: AP-2alpha knockout mice exhibit optic cup patterning defects and failure of optic stalk morphogenesis. Hum Mol Genet. 2010, 19: 1791-804. 10.1093/hmg/ddq060.PubMed CentralView ArticlePubMedGoogle Scholar
- Cheng CW, Yan CH, Hui CC, Strahle U, Cheng SH: The homeobox gene irx1a is required for the propagation of the neurogenic waves in the zebrafish retina. Mech Dev. 2006, 123: 252-63. 10.1016/j.mod.2005.12.001.View ArticlePubMedGoogle Scholar
- Choy SW, Cheng CW, Lee ST, Li VW, Hui MN, Hui CC, Liu D, Cheng SH: A cascade of irx1a and irx2a controls shh expression during retinogenesis. Dev Dyn. 2010, 239: 3204-14. 10.1002/dvdy.22462.View ArticlePubMedGoogle Scholar
- Li W, Cornell RA: Redundant activities of Tfap2a and Tfap2c are required for neural crest induction and development of other non-neural ectoderm derivatives in zebrafish embryos. Dev Biol. 2007, 304: 338-54. 10.1016/j.ydbio.2006.12.042.PubMed CentralView ArticlePubMedGoogle Scholar
- Li S, Mo Z, Yang X, Price SM, Shen MM, Xiang M: Foxn4 controls the genesis of amacrine and horizontal cells by retinal progenitors. Neuron. 2004, 43: 795-807. 10.1016/j.neuron.2004.08.041.View ArticlePubMedGoogle Scholar
- Westerfield M: The zebrafish book: a guide for the laboratory use of zebrafish (Danio rerio). 2000, Eugene, OR: University of Oregon Press, 4Google Scholar
- Nusslein-Volhard C, Dahm R, (eds): Zebrafish: a practical approach. 2002, Oxford; New York: Oxford University Press, 1Google Scholar
- Kimmel CB, Ballard WW, Kimmel SR, Ullmann B, Schilling TF: Stages of embryonic development of the zebrafish. Dev Dyn. 1995, 203: 253-310. 10.1002/aja.1002030302.View ArticlePubMedGoogle Scholar
- Leung YF, Dowling JE: Gene expression profiling of zebrafish embryonic retina. Zebrafish. 2005, 2: 269-83. 10.1089/zeb.2005.2.269.View ArticlePubMedGoogle Scholar
- Rozen S, Skaletsky HJ: Primer3 on the WWW for general users and for biologist programmers. Bioinformatics Methods and Protocols: Methods in Molecular Biology. Edited by: Krawetz S, Misener S. 2000, Totowa, NJ Humana Press, 365-386.Google Scholar
- Suzuki R, Shimodaira H: Pvclust: an R package for assessing the uncertainty in hierarchical clustering. Bioinformatics. 2006, 22: 1540-2. 10.1093/bioinformatics/btl117.View 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.