Evolutionary relationships and diversification of barhl genes within retinal cell lineages
© Schuhmacher et al; licensee BioMed Central Ltd. 2011
Received: 12 June 2011
Accepted: 21 November 2011
Published: 21 November 2011
Basic helix-loop-helix and homeodomain transcription factors have been shown to specify all different neuronal cell subtypes composing the vertebrate retina. The appearance of gene paralogs of such retina-specific transcription factors in lower vertebrates, with differently evolved function and/or conserved non-coding elements, might provide an important source for the generation of neuronal diversity within the vertebrate retinal architecture. In line with this hypothesis, we investigated the evolution of the homeobox Barhl family of transcription factors, barhl1 and barhl2, in the teleost and tetrapod lineages. In tetrapod barhl2, but not barhl1, is expressed in the retina and is important for amacrine cell specification. Zebrafish has three barhl paralogs: barhl1.1, barhl1.2 and barhl2, but their precise spatio-temporal retinal expression, as well as their function is yet unknown.
Here we performed a meticulous expression pattern comparison of all known barhl fish paralogs and described a novel barhl paralog in medaka. Our detailed analysis of zebrafish barhl gene expression in wild type and mutant retinas revealed that only barhl1.2 and barhl2 are present in the retina. We also showed that these two paralogs are expressed in distinct neuronal lineages and are differently regulated by Atoh7, a key retinal-specific transcription factor. Finally, we found that the two retained medaka fish barhl paralogs, barhl1 and barhl2, are both expressed in the retina, in a pattern reminiscent of zebrafish barhl1.2 and barhl2 respectively. By performing phylogenetic and synteny analysis, we provide evidence that barhl retinal expression domain is an ancestral feature, probably lost in tetrapods due to functional redundancy.
Functional differences among retained paralogs of key retina-specific transcription factors between teleosts and tetrapods might provide important clues for understanding their potential impact on the generation of retinal neuronal diversity. Intriguingly, within teleosts, retention of zebrafish barhl1.2 and its medaka ortholog barhl1 appears to correlate with the acquisition of distinct signalling mechanisms by the two genes within distinct retinal cell lineages. Our findings provide a starting point for the study of barhl gene evolution in relation to the generation of cell diversity in the vertebrate retina.
The vertebrate retina is organized into a complex network of cell layers, namely the ganglion cell layer (GCL) which contains retinal ganglion cells (RGCs) and displaced amacrine cells (ACs), the inner nuclear layer (INL) which consists of ACs, horizontal, bipolar and Müller glia cells, and the outer nuclear layer (ONL) which is made up of cone and rod photoreceptors. This strikingly complex architectural plan of the retina is extremely well conserved across vertebrate species, probably in direct correlation with the conservation of the key regulatory factors that govern retinal development. Several members of the basic helix-loop-helix (bHLH) and homeodomain family of transcription factors are known to play a role in the determination of retinal progenitor competence and cell fate, a function that is highly conserved from fish to mammals . Much less is known on the contribution of different functional paralogs of retina-specific transcription factors, which arose subsequently to rounds of whole genome duplication (WGD) during vertebrate evolution . Indeed, it has been proposed that after WGD, duplicated genes can either accumulate loss-of-function mutations and are functionally lost (non-functionalization [3, 4]) or acquire a new function (neo-functionalization), or split the ancestral function between the paralogs (sub-functionalization) ), therefore adding complexity to the developmental gene network that shapes organ formation. The genes of the barhl family encoding the homeobox transcription factors Barhl1 and Barhl2, have been shown to be expressed in more or less overlapping domains of the central nervous system and have partially redundant functions in neural subtype cell identity, migration and survival [5, 6]; however, barhl2 members appear to be uniquely expressed in the retina [7, 8]. In particular, Barhl2 is a pan-vertebrate regulator of the specification and survival of ACs and RGCs [9–11]. Forced expression of Barhl2 in the mouse retina promotes the differentiation of glycinergic amacrine cells at the expense of bipolar and Müller cells . Additionally, analysis of Barhl2-null retinas suggests that Barhl2 plays a critical role in both AC subtype determination and in RGC survival . The Xenopus Barhl2 ortholog (previously named Xbh1) has been shown to be expressed in RGCs and in presumptive AC precursors, and to promote RGC differentiation downstream of the bHLH transcription factor Atoh7 . While Xenopus, mouse, rat and human have one copy of barhl1 and barhl2 each, zebrafish has three barhl paralogs possibly due to a further genome duplication event that teleosts underwent during evolution after the split from the tetrapod lineage [12, 13]. On the basis of protein sequence alignment and phylogenetic analysis, it has been suggested that two of these orthologs belong to the barhl1 paralog group (nominated barhl1.1 and barhl1.2) while the third belongs to the barhl2 group [6, 12]. In contrast to mouse and Xenopus, all three barhl seem to be expressed both in the brain and in the retina . In medaka fish (Oryzias latipes), only one barhl has been described so far (olbarhl); based on phylogenetic analysis olbarhl has been clustered to the barhl1 group of paralogs [6, 12, 14]. To get more insights into the evolution of barhl paralogs with respect to retinal differentiation we took advantage of the zebrafish and medaka model systems to perform a comprehensive comparative analysis of barhl gene expression as compared to the one in tetrapod. By detailed in situ hybridization analysis we confirmed that barhl1.2 and barhl2 are expressed in the zebrafish retina but not barhl1.1. A meticulous inspection of barhl1.2 and barhl2 transcript distribution indicates that while barhl2 appears to recapitulate the expression of its mammalian and Xenopus counterpart, the spatio-temporal expression pattern of barhl1.2 is non redundant to that of barhl2, suggesting that barhl1.2 might have evolved non redundantly with respect to barhl2 in the retina. Interestingly, we have detected that barhl1.2 shows a very early onset of expression which is highly overlapping with the expression of the atoh7 gene demarcating the onset of RGC genesis. Furthermore, we also describe a new barhl paralog in medaka which, based on protein alignment, could be assigned to the Barhl2 group. By comparing the expression of medaka barhl1 and barhl2, we found that they are both expressed in the retina, in a pattern reminiscent of zebrafish barhl1.2 and barhl2 respectively. These results combined with phylogenetic and synteny analysis suggest that barhl retinal expression domain is an ancestral feature that has been specifically lost in tetrapod probably due to functional redundancy following the duplication-supplementation paradigm [3, 4].
Zebrafish barhl1.2 and barhl2, but not barhl1.1, are expressed in distinct spatio-temporal domains of the developing retina
Expression of barhl1.2 and barhl2 in the retina suggests distinct regulation by atoh7
Identification of a novel barhl paralog in medaka
Conserved gene synteny between barhl genes
According to this study, we propose that vertebrates have two homologs of barhl (barhl1 and barhl2) due to WGD events that occurred before the emergence of vertebrates. This is in contrast with Drosophila, where two BarH genes probably arose from tandem duplication in the same locus. It is a topic of speculation whether there have been one or two rounds of WGD in early vertebrate evolution before the split of the teleosts from the other jawed vertebrates, and currently the most widely accepted view is that there were two [29, 30]. Our findings are consistent with the latter view as the third paralog found in zebrafish most likely originated from the additional round of whole genome duplication that marked the rise of teleosts and did not affect other vertebrates . This hypothesis is in agreement with the fact that in basal deuterostomes such as sea squirt and amphioxus, only one barhl locus could be found by BLAST search. From this speculation one would assume that all teleost fish once had at least three, maybe even four paralogs of barhl, some of which were lost. We found that in zebrafish both barhl2 and barhl1.2 are expressed in the developing retina. Barhl1.1, like other vertebrate barhl1 homologs, is not expressed in the retina whereas the medaka barhl1 and barhl2 are both expressed in the retina in a similar fashion as barhl1.2 and barhl2 in zebrafish, thus suggesting that medaka barhl1 has a similar dynamic of transcriptional regulation as zebrafish barhl1.2. Given that the barhl expression in the retina is an ancestral feature that was already present in Drosophila, the most parsimonious explanation for barhl1 gene evolution is that its retinal expression was maintained within the teleost lineage but was lost in Xenopus and mammals. Thereafter within teleosts, the zebrafish barhl1.1 also lost its retinal expression probably due to a redundant function with barhl1.2 and relaxed evolutionary pressure in its locus. This can be further illustrated in the context of the Tetraodon barhl1 paralogs. These sequences do only resemble partial or split Barhl proteins and therefore the positions in the tree are excluded from the teleost barhl1 cluster. We have called the split sequences Tetraodon Barhl1.2 and Barhl1.3 and the protein that clusters within the teleost group Barhl1.1. In contrast, Barhl2 protein sequences are arranged in correspondence with the evolution of the species, while the retinal expression of Barhl2 has been conserved throughout the animal kingdom. This could be related to changes in function or regulation of barhl1 paralogs that led to diversification and, in some cases, retention rather than loss of a duplicate. Unfortunately, without access to expression information for Tetraodon and stickleback, no further hypothesis on the relationship of function and copy number can be made.
In line with these observations, our in situ hybridization analysis suggests that the zebrafish barhl2 expression pattern closely resembles that of barhl2 ortholog in other vertebrates, and therefore might play a role in the same neuronal lineages. The mouse Barhl2 takes part in the specification of the ACs, and later aspects of RGC maturation [9, 10]. Functional experiments will be necessary to test this hypothesis in zebrafish. It also remains to be demonstrated whether the expression of barhl2 within the GCL is restricted to displaced ACs or RGCs or both, and whether part of barhl2 expression might be dependent on atoh7, like in other vertebrates . As the expression of atoh7 and barhl2 are mostly complementary, our hypothesis is that barhl2 transcriptional activation in the AC lineage is independent from Atoh7. Other candidate factors might be at work in inducing barhl2, particularly those factors that have been shown to favour AC fate at the expense of RGCs in zebrafish . The striking complementary expression of atoh7 and barhl2 that we observed (see Figure 4) would rather supports a negative feedback between the two genes. Interestingly, similar mutually excluding expression domains have been found also in Drosophila between the homolog atonal and bar genes . In this study, a negative feedback has been shown to occur between the two Drosophila genes, thus suggesting the conservation of an ancestral feature. The expression of barhl2 in the zebrafish atoh7-/- retina at 40 hpf appears to be retained in time but its pattern of expression is shifted towards the lens (Figure 5B). Most ACs born at this stage become displaced in the presumptive GCL (which is devoid of RGCs, ); which might account for the localization of barhl2 expression in this domain. Thus, our current data doesn't provide, neither it excludes evidences for negative feedback between atoh7 and barhl2 in zebrafish. In contrast to barhl1.1, we found that barhl1.2 is expressed in the zebrafish retina. Surprisingly, this expression occurs nearly in synchrony with the one of atoh7 and the development of RGCs. Similarly, we found that the medaka barhl1 spatio-temporal expression is reminiscent of what has been previously described for atoh7 expression in the medaka retina, being first detected at stage 25 and becoming restricted to the CMZ at stage 30-35 . This observation further supports the hypothesis of conserved dynamics of retinal lineages specification in both fish "twins", at least with respect to atoh7 and barhl genes. Thus, due to this extreme similarity between barhl1.2 and atoh7 expression within the RGC lineage, it is also possible that Barhl1.2 function became redundant with respect to the one of Atoh7 and was therefore lost in the tetrapod lineage. In light of this intriguing hypothesis, it will be very interesting to ask what is the retained function of the zebrafish barhl1.2 in the RGC lineage and in relation to atoh7.
After branching of the metazoan clades into protostomes and deuterostomes, a tandem replication in Drosophila melanogaster led to the existence of two Bar factors in the fruit fly that are assumed to be redundant. In other insects, such as Drosophila ananassae, Anopheles or Apis melifera, only one Bar can be found (results of TBLASTN search). In the line of the deuterostomes, at least one round of WGD took place before the rise of the gnathostomes (jawed vertebrates). We assume that before this period of intense genome rearrangement, there was only one Barhl, as we can find in the genome sequences of Branchiostoma floridae (amphioxus) or Ciona savignyi (sea squirt, results of TBLASTN search). Going further up in the evolutionary tree, at least two Barhls could be found in every species we looked at. These observations raise the intriguing hypothesis that Barhl factors have evolved in complexity proportional to the evolving diversity of retinal arrangement and rhabdomeric photoreceptor-derived cell types . In summary, our study provides an example on how retained gene paralogs might have evolved in contributing to the specification of distinct cell-lineages in the vertebrate retina. A larger scale analysis of the functional implications of the differences in retinal key transcription factors between teleosts and other vertebrates will help to speculate further on the role of duplication retention in the evolution of vertebrates, and more specifically during eye evolution.
In summary, our teleost comparison of Barhl orthologs highlights differences in expression patterns within retinal cell populations and regarding Atoh7 regulation. By extensive cross-species analysis of the barhl loci, these differences could be linked to differential selective pressure: while the barhl2 locus remains under evolutionary constraint, we show that the barhl1 locus rapidly evolves, thereby leading to functional differences within barhl1 paralogs. Our cross-species analysis provides insights on how retained gene paralogs might have evolved in relation to distinct cell-lineages in the vertebrate retina.
Breeding and rising of zebrafish followed standard protocols . Zebrafish embryos were treated with 0.0045% 1-Phenyl-2-Thiourea (PTU) in medium after gastrulation to prevent pigment formation. Medaka embryos were kept in ERM medium containing 1 g/l NaCl, 30 mg/l KCl, 40 mg/l CaCl2 × 2 H2O and 163 mg/l MgSO4 × 7 H2O in deionized water.
All fish are housed in the fish facility of our laboratory, which was built according to the local animal welfare standards (Tierschutzgesetz 111, Abs. 1, Nr. 1) and in accordance with European Union animal welfare guidelines. The facility is under the supervision of the local representative of the animal welfare agency. No animal experiments were performed. Embryos of medaka (Oryzias latipes) and zebrafish (Danio rerio) were used exclusively at stages prior to hatching (not considered as animals according to German law and European union regulations). Zebrafish and medaka were raised and maintained as described previously . The following strains were used for zebrafish embryos: wild type WIK/AB and for medaka embryos: the Cab wild type strain."
Amplification and cloning of medaka Barhl2
An 800 bp fragment was polymerase chain reaction (PCR)-amplified from mixed stage medaka cDNA eye and brain using the forward primer (GAGATAGACACCGTGGGAACTGG) and reverse primer (CTGATGGAGTCCGGTACATGCTG) designed to bind in exons 1 and 4 of olbarhl2 (ENSORLT00000002844, EnsEMBL v58). Cycling conditions: five cycles 95°C, 10 sec, 65°C, 20 sec, 72°C, 4 min; followed by 28 cycles with annealing at 60°C. A Taq DNA-Polymerase was used for A-Tailing, incubation was for 30 min at 72°C. The PCR product was cloned into pCRII TOPO TA vector (Invitrogen) and sequenced using T7 and SP6 promoters. The sequence has been submitted to the EMBL [GenBank: JQ008931].
Whole-mount in situ hybridization
Single whole-mount in situ hybridization of barhl genes was carried out as previously described in , for the zebrafish embryos and in , for the medaka embryos. Riboprobes where labelled with digoxigenin-UTP (Roche Applied Science). Hybridization with the probe was carried on over night at 65°C/68°C. Anti-DIG primary antibody coupled to alkaline phosphatase (Roche Molecular Biochemicals) and NBT-BCIP (Roche Molecular Biochemicals) was used for signal detection. For the FISH on zebrafish embryos, modifications were applied to the method described in , as suggested by Stephan Kirchmaier (unpublished). Embryos were washed with 100 μl Tyramide Signal Amplification (TSA) solution and incubated with FITC in TSA. Incubation with the barhl2 probe was for 30 minutes, incubation with the atoh7 probe was for 40 minutes. Embryos were then kept in the dark for all following steps. For detection and staining of the antisense probes, embryos were washed 5 × 10 min with TNT (0.1M Tris pH7.5, 0.15M NaCl, 0.1% Tween20), incubated with 1% H2O2 in TNT for 20 min and washed again 5 × 10 min. Embryos were blocked in TNB (2% DIG Block in TNT) for 1 h at RT and afterwards incubated with Anti-Digoxigenin-POD Fab fragments diluted 1:100 in TNB. For signal detection, Fluorescein (FITC), Cyanine 3 (Cy3) or Cyanine 5 (Cy5) Fluorophore Tyramide by PerkinElmer was used. Embryos were then incubated in 1 × 4',6-Diamidin-2-phenylindol (DAPI) in TNT over night at 4°C and washed several times in TNT the next day. Embryos stained with NBT/BCIP were mounted in 87% Glycerol on microscope slides and imaged with a Leica DM5000B, 10x or 20x air objectives, Leica CD500 microscope camera and Leica FireCam 1.7.1 software. Double fluorescent embryos were mounted in 100 × 15 mm glass bottom dishes in 1.5% low melting agarose. Confocal stacks were taken using the Leica SP5 confocal microscope, 20x water immersion objective and Leica Application Suite (LAS) software. FITC was excited at 488 nm by Argon laser, Cy3 by the 568 nm Helium-Neon laser, Cy5 at 633 nm by Helium-Neon and DAPI by an UV laser. Emission was sensed at 500-550 nm for FITC, 650-700 nm for Cy3, 650-800 nm for Cy5 and 400-500 nm for DAPI. Emission channels were imaged sequentially to avoid bleed-through of the two fluorescent signals. Pictures were processed using the open source software ImageJ version 1.43 and Adobe Photoshop CS3.
Protein sequences were obtained from EnsEMBL Genome Browser (v58) after using TBLASTN to perform a search for the zebrafish protein sequence of Barhl2 and Barhl1.1, respectively, against the genomic DNA of the used species to check the integrity of the annotated protein sequences. The following protein sequences were used: Danio rerio Barhl1.1 (Accession number ENSDARP00000016114) Danio rerio Barhl1.2 (Acc. No. ENSDARP00000051473) Danio rerio Barhl2 (Acc. No. ENSDARP00000093436) Homo sapiens BARHL1 (Acc. No. ENSP00000263610) Homo sapiens BARHL2 (Acc. No. ENSP00000359474) Mus musculus BARHL1 (Acc. No. ENSMUSP00000053147) Mus musculus BARHL2 (Acc. No. ENSMUSP00000084005) Xenopus tropicalis barhl1 (Acc. No. ENSXETP00000013720) Xenopus tropicalis barhl2 (Acc. No. ENSXETP00000051744) Gasterosteus aculeatus Barhl1 (ENSGACP00000023974) Gasterosteus aculeatus Barhl2 (ENSGACP00000005730) Tetraodon nigroviridis Barhl2 (ENSTNIP00000016761) Tetraodon nigroviridis Barhl1.1 (ENSTNIP00000008103) Tetraodon nigroviridis Barhl1.2 (ENSTNIP00000008179) Tetraodon nigroviridis Barhl1.3 (ENSTNIP00000008783) Takifugu rubripes Barhl1 (ENSTRUP00000013477) Takifugu rubripes Barhl2 (ENSTRUP00000032575) Oryzias latipes Barhl1 (Acc. No. ENSORLP00000015485) Ciona savignyi Barhl (ENSCSAVP00000019219) Drosophila melanogaster Bar1 (Acc. No. FBpp0074204) Drosophila melanogaster Bar2 (Acc. No. FBpp00742043) Sequences obtained from National Center for Biotechnology's (NCBI) GenBank database (Release 178): Xenopus laevis Barhl1 (Acc. No. AAG14451.1) Xenopus laevis Barhl2 (Acc. No. NP_001082021.1) Salmo salar Barhl1 (Acc. No. NP_001167081.1) Branchiostoma floridae Barhl1 (XP_002596391.1). The protein sequence of Oryzias latipes Barhl2 was obtained from translating the sequenced cDNA clone using the ExPASy online translate tool from Swiss Institute of Bioinformatics (SIB). Alignment of the sequences was produced using MUSCLE online at the European Bioinformatics Institute (EBI) . A phylogenetic tree was assembled using BioNJ online at phylogeny.fr performing 1000 bootstraps and using Jones-Taylor-Thornton matrix. This distance-based algorithm is claimed to be well suited for comparison of sequences with high substitution rates [36, 37]. The tree was visualized using the open source software Dendroscope .
Chromosomal loci of barhl genes in human, zebrafish, stickleback, Tetraodon and medaka were compared by identifying all genes that occur in more than one of the loci. The position of each of these genes was then searched in all species using the EnsEMBL database search function (EnsEMBL release v58). Position and identity of genes were plotted schematically according to their order and orientation (Figure 9, not to scale). Exact position of each gene can be found in the Additional file 1: Table 1.
We would like to thank S. Kirchmaier for providing the double FISH protocol; M. Eichenlaub for providing the medaka cDNA and helping with the in situ hybridization; B. Wittbrodt, E. Leist, A. Saraceno, M. Majewski, B. Seiferling, T. Kellner, L. Schertel, C. Mueller for fish maintenance and technical assistance. We are grateful to L. Ettwiller, M. Eichenlaub, P. Jusuf, G. Lupo, M. Carl, W.A. Harris, S. Schneider and J. Wittbrodt for discussion and suggestions on the manuscript, and to J. Wittbrodt for generous support. SA is funded by LGFG (Funding program of the State of Baden-Württemberg). This work was supported by the DFG (Eigenen Stelle to LP).
- Wang JC, Harris WA: The role of combinational coding by homeodomain and bHLH transcription factors in retinal cell fate specification. Dev Biol. 2005, 285: 101-115. 10.1016/j.ydbio.2005.05.041.View ArticlePubMed
- Goode DL, Cooper GM, Schmutz J, Dickson M, Gonzales E, Tsai M, Karra K, Davydov E, Batzoglou S, Myers RM, Sidow A: Evolutionary constraint facilitates interpretation of genetic variation in resequenced human genomes. Genome Res. 2010, 20: 301-310. 10.1101/gr.102210.109.View ArticlePubMedPubMed Central
- Force A, Lynch M, Pickett FB, Amores A, Yan YL, Postlethwait J: Preservation of duplicate genes by complementary, degenerative mutations. Genetics. 1999, 151: 1531-1545.PubMedPubMed Central
- Prince VE, Pickett FB: Splitting pairs: the diverging fates of duplicated genes. Nat Rev Genet. 2002, 3: 827-837. 10.1038/nrg928.View ArticlePubMed
- Kawauchi D, Muroyama Y, Sato T, Saito T: Expression of major guidance receptors is differentially regulated in spinal commissural neurons transfated by mammalian Barh genes. Dev Biol. 2010, 344: 1026-1034. 10.1016/j.ydbio.2010.06.025.View ArticlePubMed
- Reig G, Cabrejos ME, Concha ML: Functions of BarH transcription factors during embryonic development. Dev Biol. 2007, 302: 367-375. 10.1016/j.ydbio.2006.10.008.View ArticlePubMed
- Patterson KD, Cleaver O, Gerber WV, White FG, Krieg PA: Distinct expression patterns for two Xenopus Bar homeobox genes. Dev Genes Evol. 2000, 210: 140-144. 10.1007/s004270050020.View ArticlePubMed
- Saito T, Sawamoto K, Okano H, Anderson DJ, Mikoshiba K: Mammalian BarH homolog is a potential regulator of neural bHLH genes. Dev Biol. 1998, 199: 216-225. 10.1006/dbio.1998.8889.View ArticlePubMed
- 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.View ArticlePubMedPubMed Central
- Mo Z, Li S, Yang X, Xiang M: Role of the Barhl2 homeobox gene in the specification of glycinergic amacrine cells. Development. 2004, 131: 1607-1618. 10.1242/dev.01071.View ArticlePubMed
- Poggi L, Vottari T, Barsacchi G, Wittbrodt J, Vignali R: The homeobox gene Xbh1 cooperates with proneural genes to specify ganglion cell fate within the Xenopus neural retina. Development. 2004, 131: 2305-2315. 10.1242/dev.01099.View ArticlePubMed
- Colombo A, Reig G, Mione M, Concha ML: Zebrafish BarH-like genes define discrete neural domains in the early embryo. Gene Expr Patterns. 2006, 6: 347-352. 10.1016/j.modgep.2005.09.011.View ArticlePubMed
- Furlong RF, Holland PW: Were vertebrates octoploid?. Philos Trans R Soc Lond B Biol Sci. 2002, 357: 531-544. 10.1098/rstb.2001.1035.View ArticlePubMedPubMed Central
- Poggi L, Carl M, Vignali R, Barsacchi G, Wittbrodt J: Expression of a medaka (Oryzias latipes) Bar homolog in the differentiating central nervous system and retina. Mech Dev. 2002, 114: 193-196. 10.1016/S0925-4773(02)00054-0.View ArticlePubMed
- Jusuf PR, Almeida AD, Randlett O, Joubin K, Poggi L, Harris WA: Origin and determination of inhibitory cell lineages in the vertebrate retina. J Neurosci. 2011, 31: 2549-2562. 10.1523/JNEUROSCI.4713-10.2011.View ArticlePubMedPubMed Central
- Masai I, Stemple DL, Okamoto H, Wilson SW: Midline signals regulate retinal neurogenesis in zebrafish. Neuron. 2000, 27: 251-263. 10.1016/S0896-6273(00)00034-9.View ArticlePubMed
- Brown NL, Patel S, Brzezinski J, Glaser T: Math5 is required for retinal ganglion cell and optic nerve formation. Development. 2001, 128: 24972508-
- Ghiasvand NM, Rudolph DD, Mashayekhi M, Brzezinski JAt, Goldman D, Glaser T: Deletion of a remote enhancer near ATOH7 disrupts retinal neurogenesis, causing NCRNA disease. Nat Neurosci. 2011, 14: 578-586. 10.1038/nn.2798.View ArticlePubMedPubMed Central
- 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-736. 10.1016/S0896-6273(01)00312-9.View ArticlePubMed
- Wang SW, Kim BS, Ding K, Wang H, Sun D, Johnson RL, Klein WH, Gan L: Requirement for math5 in the development of retinal ganglion cells. Genes Dev. 2001, 15: 24-29. 10.1101/gad.855301.View ArticlePubMedPubMed Central
- Feng L, Xie ZH, Ding Q, Xie X, Libby RT, Gan L: MATH5 controls the acquisition of multiple retinal cell fates. Mol Brain. 2010, 3: 36-10.1186/1756-6606-3-36.View ArticlePubMedPubMed Central
- Poggi L, Vitorino M, Masai I, Harris WA: Influences on neural lineage and mode of division in the zebrafish retina in vivo. J Cell Biol. 2005, 171: 991-999. 10.1083/jcb.200509098.View ArticlePubMedPubMed Central
- Furutani-Seiki M, Wittbrodt J: Medaka and Zebrafish, an evolutionary twin study. Mech Dev. 2004, 629-637. 7-8
- Peden E, Kimberly E, Gengyo-Ando K, Mitani S, Xue D: Control of sex-specific apoptosis in C. elegans by the BarH homeodomain protein CEH-30 and the transcriptional repressor UNC-37/Groucho. Genes Dev. 2007, 21: 3195-3207. 10.1101/gad.1607807.View ArticlePubMedPubMed Central
- Iwamatsu T: Stages of normal development in the medaka Oryzias latipes. Mech Dev. 2004, 605-618. 7-8
- Souren M, Martinez-Morales JR, Makri P, Wittbrodt B, Wittbrodt J: A global survey identifies novel upstream components of the Ath5 neurogenic network. Genome Biology. 2009, 10 (9): R92-10.1186/gb-2009-10-9-r92.View ArticlePubMedPubMed Central
- Kikuta H, Laplante M, Navratilova P, Komisarczuk AZ, Engstrom PG, Fredman D, Akalin A, Caccamo M, Sealy I, Howe K, et al: Genomic regulatory blocks encompass multiple neighboring genes and maintain conserved synteny in vertebrates. Genome Res. 2007, 17: 545-555. 10.1101/gr.6086307.View ArticlePubMedPubMed Central
- Ravi V, Venkatesh B: Rapidly evolving fish genomes and teleost diversity. Curr Opin Genet Dev. 2008, 18: 544-550. 10.1016/j.gde.2008.11.001.View ArticlePubMed
- Kuraku S, Meyer A: The evolution and maintenance of Hox gene clusters in vertebrates and the teleost-specific genome duplication. Int J Dev Biol. 2009, 53: 765-773. 10.1387/ijdb.072533km.View ArticlePubMed
- Panopoulou G, Poustka AJ: Timing and mechanism of ancient vertebrate genome duplications -- the adventure of a hypothesis. Trends Genet. 2005, 21: 559-567. 10.1016/j.tig.2005.08.004.View ArticlePubMed
- Lim J, Choi KW: Induction and autoregulation of the anti-proneural gene Bar during retinal neurogenesis in Drosophila. Development. 2004, 131: 5573-5580. 10.1242/dev.01426.View ArticlePubMed
- Lim J, Choi KW: Bar homeodomain proteins are anti-proneural in the Drosophila eye: transcriptional repression of atonal by Bar prevents ectopic retinal neurogenesis. Development. 2003, 130: 5965-5974. 10.1242/dev.00818.View ArticlePubMed
- Arendt D: Evolution of eyes and photoreceptor cell types. Int J Dev Biol. 2003, 47: 563-571.PubMed
- Westerfield M: The zebrafish book. 1994
- Edgar RC: MUSCLE: multiple sequence alignment with high accuracy and high throughput. Nucleic Acids Res. 2004, 32: 1792-1797. 10.1093/nar/gkh340.View ArticlePubMedPubMed Central
- Dereeper A, Guignon V, Blanc G, Audic S, Buffet S, Chevenet F, Dufayard JF, Guindon S, Lefort V, Lescot M, et al: Phylogeny.fr: robust phylogenetic analysis for the non-specialist. Nucleic Acids Res. 2008, 36: W465-469. 10.1093/nar/gkn180.View ArticlePubMedPubMed Central
- Gascuel O: BIONJ: an improved version of the NJ algorithm based on a simple model of sequence data. Mol Biol Evol. 1997, 14: 685-695.View ArticlePubMed
- Huson DH, Richter DC, Rausch C, Dezulian T, Franz M, Rupp R: Dendroscope: An interactive viewer for large phylogenetic trees. BMC Bioinformatics. 2007, 8: 460-10.1186/1471-2105-8-460.View ArticlePubMedPubMed Central
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.