Genetic mechanisms involved in the evolution of the cephalopod camera eye revealed by transcriptomic and developmental studies
© Yoshida and Ogura; licensee BioMed Central Ltd. 2011
Received: 12 February 2011
Accepted: 24 June 2011
Published: 24 June 2011
Coleoid cephalopods (squids and octopuses) have evolved a camera eye, the structure of which is very similar to that found in vertebrates and which is considered a classic example of convergent evolution. Other molluscs, however, possess mirror, pin-hole, or compound eyes, all of which differ from the camera eye in the degree of complexity of the eye structures and neurons participating in the visual circuit. Therefore, genes expressed in the cephalopod eye after divergence from the common molluscan ancestor could be involved in eye evolution through association with the acquisition of new structural components. To clarify the genetic mechanisms that contributed to the evolution of the cephalopod camera eye, we applied comprehensive transcriptomic analysis and conducted developmental validation of candidate genes involved in coleoid cephalopod eye evolution.
We compared gene expression in the eyes of 6 molluscan (3 cephalopod and 3 non-cephalopod) species and selected 5,707 genes as cephalopod camera eye-specific candidate genes on the basis of homology searches against 3 molluscan species without camera eyes. First, we confirmed the expression of these 5,707 genes in the cephalopod camera eye formation processes by developmental array analysis. Second, using molecular evolutionary (dN/dS) analysis to detect positive selection in the cephalopod lineage, we identified 156 of these genes in which functions appeared to have changed after the divergence of cephalopods from the molluscan ancestor and which contributed to structural and functional diversification. Third, we selected 1,571 genes, expressed in the camera eyes of both cephalopods and vertebrates, which could have independently acquired a function related to eye development at the expression level. Finally, as experimental validation, we identified three functionally novel cephalopod camera eye genes related to optic lobe formation in cephalopods by in situ hybridization analysis of embryonic pygmy squid.
We identified 156 genes positively selected in the cephalopod lineage and 1,571 genes commonly found in the cephalopod and vertebrate camera eyes from the analysis of cephalopod camera eye specificity at the expression level. Experimental validation showed that the cephalopod camera eye-specific candidate genes include those expressed in the outer part of the optic lobes, which unique to coleoid cephalopods. The results of this study suggest that changes in gene expression and in the primary structure of proteins (through positive selection) from those in the common molluscan ancestor could have contributed, at least in part, to cephalopod camera eye acquisition.
Animal eyes have long been considered a classic example of convergent evolution. In recent decades, this view has changed due to the discovery of shared developmental regulatory genes for eye formation. Several genes, such as Pax-6/eyeless (ey) , eyes absent , dachshund , and sine oculis , together with their orthologs in metazoan animals, are able to induce the formation of ectopic eyes in flies and have been regarded as essential eye regulator genes among metazoan animals [5, 6]. Most of the genes involved in eye development had already existed in the last common ancestors of cnidarians and bilaterians . Such evidence suggests that some conserved genes have similarly contributed to eye development across a wide range of animals.
The question remains as to what genetic background could have contributed to the evolution of the coleoid cephalopod camera eye. We hypothesized that changes in expresssion patterns and functions of pre-existing genes as well as the gain and loss of genes have played important roles in the evolution of the camera eye in coleoid cephalopods. Previously, Ogura et al.  revealed that more than 60% of octopus eye ESTs were commonly observed in the human eye, indicating that changes of expression patterns and functions have been involved in the evolution of the camera eyes in the octopus and humans.
Identifying genes specifically expressed in the coleoid cephalopod camera eye
To collect cephalopod camera eye-specific candidate genes, we performed an in vitro homology search using an array designed to estimate gene expression among species on the basis of inter-species hybridization signals. In this in vitro homology search array, probes were designed, using the probe design method described in Ogura et al. , to estimate distantly related homologs of molluscan eye genes. Microarray probes from Octopus vulgaris (octopus), Loligo bleekeri (squid), Idiosepius paradoxus (pygmy squid), Mizuhopecten (Patinopecten) yessoensis (scallop) and Nautilus pompilius (nautilus) were derived from our eye EST libraries and the NCBI EST database. We applied total RNA samples from the pygmy squid, scallop and nautilus to a custom Agilent microarray designed to detect specific gene expression in molluscan eyes. The microarray results for camera eye-specific genes allowed us to eliminate 173 and 162 genes found to be significantly expressed, using the Feature Extraction output (see methods), in the nautilus and scallop eye RNA, respectively. In addition to elimination on the basis of the microarray results, we also eliminated genes having homolog(s) among the nautilus or scallop eye ESTs, based on BLASTN and TBLASTX searches (<1e-4). As a result, a total of 5,707 genes (2,607, 1,014 and 2,086 genes from the pygmy squid, squid, and octopus, respectively) were selected as cephalopod camera eye-specific candidate genes (Figure 2). Although this array-based estimation may include false-positive genes (expressed in nautilus or scallop eyes but not detected in the array), as described in Gilad et al. , we expect highly conserved genes (common eye genes and housekeeping genes) to be eliminated using this analysis.
Confirmation of cephalopod camera eye-specific genes expression in the developmental stages of the pygmy squid by developmental array
Next, we performed microarray experiments using pygmy squid embryonic RNA to estimate genes contributing to the camera eye developmental process in coleoid cephalopods. We extracted total RNAs from whole embryos at stage 20 and stage 30 (hatchlings), and from the eyes at stage 25. The reason for the selection of these three stages is that essential developmental events occur in each stage: stage 20 sees the establishment of eye placodes, stage 25 sees eye vesicle closure and the onset of lens formation, and stage 30 sees the completion of eye development (though genes used for the growth and maintenance of eyes are still active). Among the cephalopod camera eye-specific candidate genes selected above, 53.8% (3,075/5,707 genes) of the genes were found to have a positive signal against pygmy squid embryonic RNA.
Estimation of lineage-specific functional changes in coleoid cephalopod camera eye-specific genes by molecular evolutionary analysis
We hypothesized that taxon-specific gene modification contributed to the evolutionary steps in coleoid cephalopod camera eye development. To clarify the involvement of gene recruitment and subsequent positive selection in camera eye acquisition, we estimated the dN/dS ratio of camera eye-specific candidate genes after the divergence of cephalopods from the common molluscan ancestor. The cephalopod camera eye-specific data were concatenated together with homologous genes estimated by BLASTN and TBLASTX searches against the molluscan NCBI nucleotide collection and molluscan ESTs (an e-value of less than 1e-20). Phylogenetic trees based on the gene sets were analyzed using the Codeml program in PAML ver. 4.4 b  under the F3x4 model. Finally, we obtained 1,391 sets of multiple alignments including cephalopods and other molluscs, and dN/dS values of the camera eye-specific genes. Of the camera eye-specific genes, 156 (68, 43, 24 from the pygmy squid, octopus, and squid, respectively) were identified as positively selected genes in the cephalopod lineage (the dN/dS ratio was significantly higher than 2, Additional file 1, Table S1).
Homologs of the cornichon gene were found among the octopus and the pygmy squid ESTs (5primeCluster0328, DB919298), and both homologs were found to be positively selected genes in the cephalopod lineage. Cornichon-like genes in vertebrates have been shown to be located directly downstream of Pax-6 . It is of great interest that the region downstream of Pax-6 is included among the camera eye-specific genes, although it is still unclear whether the squid cornichon is controlled by Pax-6 or not.
Centaurin gamma is a small GTP family member with NTPase activity  (Additional file 2, Figure S1, Additional file 3, Supplementary data S1). It is known to be expressed in human cancer cells is known; however, its function in vertebrate development is unclear. Localization of the centaurin gamma homolog (OctEye_3919F_082) is described below.
Determination of common and unique genes by comparative gene expression analysis of cephalopods and vertebrates
Over-represented GO terms in the camera eye-specific genes.
protein binding (F)
nucleotide binding (F)
structural constituent of ribosome (F)
intracellular protein transport (P)
ATP binding (F)
translation initiation factor activity (F)
GTP binding (F)
signal transduction (P)
zinc ion binding (F)
small GTPase mediated signal transduction (P)
protein transport (P)
metal ion binding (F)
catalytic activity (F)
protein homodimerization activity (F)
protein modification process (P)
transferase activity (F)
Functional analysis of four candidate genes selected among the coleoid cephalopod camera eye-expressed genes
To clarify the function of the coleoid cephalopod camera eye-specific candidate genes in the camera eye formation process, we cloned and analyzed the localization of transcripts of three candidate genes in the pygmy squid embryo. We selected the Ets-4, HMGb3 and Hla-b associated transcript homologs from the transcription factor genes commonly expressed in vertebrates and cephalopods due to their high expression signals and the existence of human orthologs for further functional analysis (Additional file 1, Table S1). The centaurin gamma homolog was also selected from the positively selected gene candidate genes in the cephalopod lineage as it showed a high dN/dS ratio and had a human ortholog (Additional file 1, Table S1, Additional file 2, Figure S1, Additional file 3, Supplementary data S1).
Genetic mechanisms involved in camera eye evolution of the coleoid cephalopods
The coleoid cephalopod camera eye-specific candidate genes selected by our analysis might be divided into two groups. The first group consists of genes that commonly exist in animals, but whose expression level has been changed through the structural changes of the eyes after the divergence of the cephalopods. The Ets-4 homolog, for example, was expressed in the granular cell layers specifically found in the optic lobes in coleoid cephalopods. The fact that a large proportion of the genes expressed in the octopus eye was found in the genomes of other animals  supports the notion that the role of the genes might have changed in relation to the structural changes in the various eye types. The second group comprises genes that commonly exist, even in molluscan eyes, but whose functions have been changed by positive selection in relation to camera eye development in cephalopods. We identified more than a hundred genes that appear to have undergone positive selection (dN/dS values of more than 2) in the coleoid cephalopod lineage. The genetic mechanism represented by these two groups seems to have been contributed to acquisition of the camera eye.
Because the in vitro homology search array could be used to detect highly expressed genes in humans and squids , we thought it possible to remove highly expressed molluscan eye genes, including housekeeping genes, from the list of candidate genes to estimate cephalopod specificity. The cephalopod camera eye-specific candidate gene set also contain genes common among animals with eyes, all of which should, ideally, be removed along with the nautilus and scallop gene pool but which nevertheless remained due to the use of incomplete gene sets.
Functions of the coleoid cephalopod camera eye-specific genes found in this study
One of the most intriguing events in coleoid cephalopod camera eye evolution is the transition from the simple eye to a more complex visual eye. The transition was accompanied by the evolution of the eye components including the optic lobes, visual center of the brain. In this transition, expression changes and functional chagnges of the camera eye specific candidate genes might have been contributed to the acquisition of the eye compnents. For example, Ets-4 localization in the five-layered optic lobes of the pygmy squid appears to be involved in the organization of outer granular layer. As Erg and Fli-1, orthologs of Ets-4 in vertebretes, are major inducers of blood vessel development  and D-ets-3 and D-ets-6, orthologs of Ets-4 in Drosophila, are known to be expressed in the ventral nerve cord [24, 25], it is speculated that Ets-4 of the pygmy squid was recruited for the development of different organs, particularly in the conversion of neuron responses in cephalopod optic lobes.
The cornichon and centaurin gamma homologs identified as camera eye specific genes in cephalopods are also known to be expressed in the vertebrate camera eye. Even though these transduction-related genes are completely function-unknown, together with their upstream and downstream, their observed localization in the optic lobe suggests their involvement in the neuronal development of the coleoid cephalopod camera eye. The HMGb3 and Hla-b associated transcript homologs may be involved in visual circuit formation through chromatin remodeling and initiation of subsequent gene expressions. As a consequence of our study, a number of common genes from the common molluscan ancestor have been shown to be involved in coleoid cephalopod camera eye evolution through changes in their expression patterns and functions.
In this study, we selected 5,707 genes as coleoid cephalopod camera eye-specific candidates by homology searches and comparative expression against the nautilus and scallop eyes. We applied a molecular evolutionary approach to the cephalopod camera eye candidates to identify genes that had functionally changed after the divergence of cephalopods from the common ancestor. Our developmental approach, which included developmental arrays and in situ hybridization using the pygmy squid embryos, indicated the contribution of camera eye-specific genes to the eye developmental process. In particular, some camera eye specific genes, such as the Ets-4, cornichon, and centaurin gamma homologs, was shown to have contributed to camera eye evolution through the development of a complex visual circuit.
Pygmy squid cultures and RNA isolation
Japanese pygmy squid, Idiosepius paradoxus, specimens were obtained from Chita Peninsula, Nagoya, Japan and kept in our laboratory, as previously described , where the embryos were staged . The embryonic eyes of stage 25 embryos were dissected from the body using forceps and used for RNA extraction. A specimen of Nautilus pompilius was purchased from a retail aquarium in Kanagawa, Japan. A specimen of the Japanese scallop, Mizuhopecten (Patinopecten) yessoensis, was purchased from the Tokyo Central Wholesale Fish Market (caught in Aomori). Total RNA was extracted from the embryos or adult tissues using an E. Z. N. A. Mollusc RNA kit (Omega Bio-Tek, Inc.), according to the manufacturer's protocol.
Probe design for the in vitrohomology search of the molluscan eyes
Molluscan EST data for Loligo bleekeri (squid), Mizuhopecten (Patinopecten) yessoensis (scallop) and Nautilus pompilius (nautilus) were obtained from cDNA libraries constructed using the Creator™ SMART™ Library Construction Kit (Takara Clontech) with the TRIMMER-DIRECT cDNA normalization kit (Evrogen) and introduced into DH5α electro competent cells (Takara). The 5' regions of the cDNA libraries were sequenced using the M13 forward primer (5-TGT AAA ACG ACG GCC AGT-3) and Big Dye Terminator v3.1 (Applied Biosystems) on an ABI3030 sequencer (Applied Biosystems). Raw data for these ESTs are available under the following accession numbers: [DDBJ:FY298839-FY302524]. The sequence clean-up program Lucy (The Institute for Genomic Research, Rockville, MD; http://www.tigr.org/softlab) was used to trim vector sequences and low-quality sequences from the raw data. These ESTs were assembled using the MIRA assembler 3.2.1 , and then used for the subsequent analyses. The octopus ESTs were obtained as previously described  and from the following website http://www.cib.nig.ac.jp/dda/database/octopus.htm. We adopted ESTs of I. paradoxus from the NCBI EST database under the following accession numbers: [DDBJ, DB910977-DB920055]. Redundant sequences in the I. paradoxus ESTs were found using BLASTCLUST and removed. These data were concatenated and used as molluscan eye ESTs in the subsequent analyses, and are available upon request. We added the pubmed data sets for vertebrate eye genes, obtained from a public database, to identify eye gene expression across distantly related species. Finally, we obtained 12,128 molluscan probes (4,657, 1,483, 3,377, 1,493 and 1,118 from the pygmy squid, squid, octopus, scallop and nautilus, respectively) and designed a custom microarray with the Agilent 8 × 15K format custom microarray service.
To perform the comparative analysis of various developmental stages, we applied total RNA samples from the molluscan species onto the custom microarray. Total RNA from the embryonic stages were labeled and applied separately. Briefly, 0.5 μg of total RNA from each sample was used to synthesize fluorescent-labelled cRNA using Cyanine 3-CTP according to the manufacturer's protocol (Agilent Quick Amp Labeling Kit, one-color). Labeled DNA was hybridized for 16 hr at 65°C on the custom array. After hybridization, the microarray slides were washed according to the standard protocol (Agilent Technologies) and scanned on an Agilent microarray scanner. Data were analyzed using the Agilent Feature Extraction Software (v10.7, Additional file 7, Table S2). We adopted the microarray data under the following accession number [CIBEX, CBX187]. We selected wells (spots) exhibited a signal "well above the background" and "positive and siginificant" as siginificantly expressed genes using the Feature Extraction output.
Databases for homology searches among molluscs and vertebrates
Molluscan eye ESTs were obtained as described above. For comparison with the vertebrate eye EST data, we used data from Bodymap, UniGene, and NEI Bank (Accession numbers are listed as Additional file 8, Supplement data S3, Additional file 9, Supplement data S4, Additional file 10, Supplement data S5, respectively). For the comparative analysis of gene expression, we obtained human-eye ESTs (929 contigs) of the retina, corneal endothelium, cornea, and iris from BodyMap . NEI Bank, maintained by the National Eye Institute, a division of the NIH (USA), contains data (680,045 contigs) for several human-eye cDNA libraries , including tissue-separated databases of the ciliary body, cornea, fovea, iris, lens, optic nerve, retina, RPE choroids, trabecular meshwork, and other tissues/organs. We also used data sets from the vertebrate lens, retinal and eye libraries (484, 202 contigs from humans, mice, rats, rabbits, dogs, pigs, sticklebacks, zebrafish, salmon, and pufferfish) in the NCBI UniGene .
Functional annotation was conducted as follows. Gene ontology, which is defined by the Gene Ontology Consortium [31, 32], and used to categorize genes by (1) Molecular Function, (2) Biological Process, or (3) Cellular Component, was used to categorize cephalopod eye ESTs for the octopus, squid and pygmy squid. We adopted all GO terms for categorizing cephalopod eye genes using a blast2GO . Enrichment analysis of the camera eye-specific genes against total cephalopod ESTs excluding the camera eye-specific genes, was carried out with the blast2GO program (using Fisher's exact test, an e-value of less than 0.01, Table 1). The Pfam database  was used to estimate the domain structure of the molluscan eye ESTs. Domains in the translated amino acid sequences of the molluscan eye ESTs were estimated using batch search from the Pfam website.
We determined transcriptional factors in the cephalopod camera eye-specific candidate genes manually using the sequence description obtained by the blast2GO, and the HMM-name obtained by the Pfam programs. Homologs of six homeobox 2 (00672_Oc_096), ets-related isoform 4 (DB913089), lim domain-containing (06182_Oc_5_043), lim and sh3 domain protein (OctEye_4576F_074) and high mobility group b3b (01791_Oc_5_056) were identified as they have homeobox, ETS, LIM, LIM and high mobility group domains. We also identified genes involved in signal transduction cascades, homologs of rab5 protein (02945_Oc_5_065), ras-related protein rab-2 (DB916866), member ras oncogene family (ika1224-No19_E06_001, DB913730 and DB919285), ran gtpase-activating protein (DB912649), and dishevelled associated activator of morphogenesis 1 like (03249_Oc_5_081), which are possibly related to camera eye specification. Indeed, the Ras superfamily of small GTPases regulate many cellular regulatory and developmental pathways, including diurnal regulation of the octopus photoreceptors . As previously reported , retinal arrestin (08322_Oc_5_073), retinal dehydrogenase (DB913426), neuron-specific enolase (5primeCluster0558), and gelsolin (07345_Oc_5_049, 08462_Oc_5_015) are commonly expressed and are thought to be involved in camera eye formation.
Molecular evolutionary (dN/dS) analysis
We calculated the dN/dS ratio in the cephalopod lineage (squids or octopus) after divergence from other molluscs to estimate genes changed through positive selection in the cephalopod camera eye. The nucleotide sequences and deduced amino-acid sequences of the cephalopod eye cDNAs isolated in the present study were aligned together with estimated homologs of other molluscs using homology searches. The homologs were selected using BLASTN and TBLASTX against the non-cephalopod molluscan eye ESTs (scallop and nautilus) and the molluscan gene sets obtained from the NCBI nucleotide and EST collections (Taxonomy ID: 6447). Homology searches of the GenBank non-redundant database using BLAST with a cutoff value of < 1e-20 produced matches against 1,391 of the contigs. The homologs were concatenated and aligned using the Clustal W program (default options) . After removing gap positions, the dN/dS ratio (Kn divergence per Ks divergence) between each homologous gene set was estimated using the maximum likelihood method of Goldman and Yang , as implemented in the codeml program of the PAML package  under the F3 × 4 model. The statistical significance of positively selected genes was determined using a one-tailed student's t-test. The above-mentioned protocol was systematically performed using both Perl and shell scripts.
Targeted-gene cloning in the pygmy squid and expression analysis
We selected the cephalopod camera eye candidate genes using the above-mentioned criteria. Orthologous candidates in the pygmy squid were then cloned to validate their expression pattern. Fragments of the Ets-4, HMGb3, Hla-b associated transcript, cornichon and centaurin gamma homologs were obtained by PCR reaction using pygmy squid eye cDNA. Reverse transcription was carried out using a PrimeScript™ RT Reagent kit (TaKaRa), and the primers used to amplify the target genes are listed in Additional file 11, Table S3.
Longer length cDNAs of the Ets-4 and centaurin gamma homologs were obtained using a BD SMARTer RACE (rapid amplification of cDNA ends) cDNA amplification kit (Clontech). Ready-to-use first strand cDNA was synthesized according to the manufacturer's protocol using pygmy squid embryonic eye RNA. Total RNA was prepared from a number of pygmy squid hatchling specimens using an RNeasy mini RNA extraction kit (QIAGEN). The primers used to amplify the cDNA extremities were Ets4-5RACE, AAGTTGTATCACCTGGGAGGGCCG, Ets4-3RACE, CTTTTTGCAGCGCCGGTTCCTATT, centaurin gamma-5RACE, ACGCGAGGATTGCTTTCACTGATGG, and centaurin gamma-3RACE, GGTCCGCCGGAGATGCAGTTTACTC. The RACE reaction products were cloned into a plasmid using a pGEM T-vector system (Promega). Plasmid DNA from transformant colonies was purified with a GenElute™ Plasmid Miniprep kit (Sigma-Aldrich). Both strands of the plasmid DNA were fully sequenced downstream of an insert site by the dideoxy chain-termination method using a BigDye® Terminator v. 3.1. The sequences are available under the following accession numbers: Ets-4 homolog, AB586703; centaurin gamma homolog, AB586704.
Short fragments of these genes were amplified by PCR and subcloned into the pGEM T-vector. These were then used as a template to generate digoxigenin (DIG)-labeled antisense and sense probes by in vitro transcription with a DIG-RNA labeling kit (Roche) using SP6 RNA polymerase (Roche) and T7 RNA polymerase (Roche). Whole-mount in situ hybridization for stage 20, 22, and 25 of the pygmy squid embryos was carried out using established procedures .
Acknowledgements and Funding
We are truly grateful to Ms. C. Satoh and Ms. S. Takagi at Ochanimizu University for their assistance with the analyses. We would like to thank Mr. T. Kasugai at the Nagoya Port Foundation for his help in the collection of the adult pygmy squid specimens, and Dr. J. Sese and Ms. M. Fukuzaki for their help in constructing and analyzing the in vitro homology search array.
This study was supported by the grant from the Japanese Science and Technology Agency to AO.
- Halder G, Callaerts P, Gehring WJ: Induction of ectopic eyes by targeted expression of the eyeless gene in Drosophila. Science. 1995, 267: 1788-1792. 10.1126/science.7892602.View ArticlePubMedGoogle Scholar
- Bonini NM, Bui QT, Gray-Board GL, Warrick JM: The Drosophila eyes absent gene directs ectopic eye formation in a pathway conserved between flies and vertebrates. Development. 1997, 124: 4819-4826.PubMedGoogle Scholar
- Shen W, Mardon G: Ectopic eye development in Drosophila induced by directed dachshund expression. Development. 1997, 124: 45-52.PubMedGoogle Scholar
- Pignoni F, Hu B, Zavitz KH, Xiao J, Garrity PA, Zipursky SL: The eye-specification proteins So and Eya form a complex and regulate multiple steps in Drosophila eye development. Cell. 1997, 91: 881-891. 10.1016/S0092-8674(00)80480-8.View ArticlePubMedGoogle Scholar
- Gehring WJ, Ikeo K: Pax6: Mastering eye morphogenesis and eye evolution. Trends Genet. 1999, 15: 371-377. 10.1016/S0168-9525(99)01776-X.View ArticlePubMedGoogle Scholar
- Donner AL, Maas RL: Conservation and non-conservation of genetic pathways in eye specification. Int J Dev Biol. 2004, 48: 743-753. 10.1387/ijdb.041877ad.View ArticlePubMedGoogle Scholar
- Koyanagi M, Takano K, Tsukamoto H, Ohtsu K, Tokunaga F, Terakita A: Jellyfish vision starts with cAMP signaling mediated by opsin-G(s) cascade. Proc Natl Acad Sci USA. 2008, 105: 15576-15580. 10.1073/pnas.0806215105.View ArticlePubMedPubMed CentralGoogle Scholar
- Oakley TH: The eye as a replicating and diverging, modular developmental unit. Trends Ecol Evol. 2003, 18: 623-627. 10.1016/j.tree.2003.09.005.View ArticleGoogle Scholar
- Serb JM, Oakley TH: Hierarchical phylogenetics as a quantitative analytical framework for Evolutionary Developmental Biology. BioEssays. 2005, 27: 1158-1166. 10.1002/bies.20291.View ArticlePubMedGoogle Scholar
- Serb JM: Toward developing models to study the disease, ecology, and evolution of the eye in Mollusca. Amer Malac Bull. 2008, 26: 3-18. 10.4003/006.026.0202.View ArticleGoogle Scholar
- Hochner B, Shomrat T, Fiorito G: The octopus: a model for a comparative analysis of the evolution of learning and memory mechanisms. Biol Bull. 2006, 210: 308-317. 10.2307/4134567.View ArticlePubMedGoogle Scholar
- Young JZ: The central nervous system of Loligo. I. The optic lobe. Philos Trans R Soc Lond B Biol Sci. 1974, 267: 263-302. 10.1098/rstb.1974.0002.View ArticlePubMedGoogle Scholar
- Kriegstein AR: Development of the nervous system of Aplysia californica. PNAS. 1977, 74: 375-378. 10.1073/pnas.74.1.375.View ArticlePubMedPubMed CentralGoogle Scholar
- Ogura A, Ikeo K, Gojobori T: Comparative analysis of gene expression for convergent evolution of camera eye between octopus and human. Genome Res. 2004, 14: 1555-1561. 10.1101/gr.2268104.View ArticlePubMedPubMed CentralGoogle Scholar
- Ogura A, Yoshida MA, Fukuzaki M, Sese J: In vitro homology search array comprehensively reveals highly conserved genes and their functional characteristics in non-genome sequenced species. BMC Bioinfomatics. 2010, 11 (Suppl 4): S9-10.1186/1471-2105-11-S4-P9.View ArticleGoogle Scholar
- Gilad Y, Rifkin SA, Bertone P, Gerstein M, White KP: Multi-species microarrays reveal the effect of sequence divergence on gene expression profiles. Genome Res. 2005, 15 (5): 674-80. 10.1101/gr.3335705.View ArticlePubMedPubMed CentralGoogle Scholar
- Yang Z: PAML 4: Phylogenetic analysis by maximum likelihood. Mol Biol Evol. 2007, 24: 1586-1591. 10.1093/molbev/msm088.View ArticlePubMedGoogle Scholar
- Wolf LV, Yang Y, Wang J, Xie Q, Braunger B, Tamm ER, Zavadil J, Cvekl A: Identification of Pax6-dependent gene regulatory networks in the mouse lens. PLoS ONE. 2009, 4: e4159-10.1371/journal.pone.0004159.View ArticlePubMedPubMed CentralGoogle Scholar
- Soundararajan M, Yang X, Elkins JM, Sobott F, Doyle DA: The centaurin γ -1 GTPase-like domain functions as an NTPase. Biochem J. 2007, 401: 679-688. 10.1042/BJ20060555.View ArticlePubMedPubMed CentralGoogle Scholar
- Randi AM, Sperone A, Dryden NH, Birdsey GM: Regulation of angiogenesis by ETS transcription factors. Biochem Soc Trans. 2009, 37: 1248-1253. 10.1042/BST0371248.View ArticlePubMedGoogle Scholar
- Budelmann BU, Schipp R, Boletzky SV: Cephalopoda. Microscopic Anatomy of Invertebrates, (Mollusca II). Edited by: Harrison FW, Ruppert EE. 1997, New York: Wiley-Liss, 6A: 119-414.Google Scholar
- Yamamoto M: Normal embryonic stages of the pygmy cuttlefish, Idiosepius pygmaeus paradoxus Ortmann. Zool Sci. 1988, 5: 989-998.Google Scholar
- Yamamoto M, Shimazaki Y, Shigeno S: Atlas of the embryonic brain in the pygmy squid, Idiosepius paradoxus. Zool Sci. 2003, 20: 163-179. 10.2108/zsj.20.163.View ArticlePubMedGoogle Scholar
- Chen T, Bunting M, Karim FD, Thummel CS: Isolation and characterization of five Drosophila genes that encode an ets-related DNA binding domain. Dev Biol. 1992, 151: 176-191. 10.1016/0012-1606(92)90225-6.View ArticlePubMedGoogle Scholar
- Hsu T, Schulz RA: Sequence and functional properties of Ets genes in the model organism Drosophila. Oncogene. 2000, 19: 6409-6416. 10.1038/sj.onc.1204033.View ArticlePubMedGoogle Scholar
- Yoshida MA, Shigeno S, Tsuneki K, Furuya H: Squid vascular endothelial growth factor receptor: a shared molecular signature in the convergent evolution of closed circulatory systems. Evol Dev. 2009, 12: 25-33.View ArticleGoogle Scholar
- Chevreux B, Pfisterer T, Drescher B, Driesel AJ, Müller WE, Wetter T, Suhai S: Using the miraEST Assembler for Reliable and Automated mRNA Transcript Assembly and SNP Detection in Sequenced ESTs. Genome Research. 2004, 14: 1147-1159. 10.1101/gr.1917404.View ArticlePubMedPubMed CentralGoogle Scholar
- Hishiki T, Kawamoto S, Morishita S, Okubo K: BodyMap: a human and mouse gene expression database. Nucleic Acids Res. 2000, 28: 136-138. 10.1093/nar/28.1.136.View ArticlePubMedPubMed CentralGoogle Scholar
- Wistow G, Peterson K, Gao J, Buchoff P, Jaworski C, Bowes-Rickman C, Ebright JN, Hauser MA, Hoover D: NEIBank: genomics and bioinformatics resources for vision research. Mol Vis. 2008, 14: 1327-1337.PubMedPubMed CentralGoogle Scholar
- Pontius JU, Wagner L, Schuler GD: UniGene: a unified view of the transcriptome. The NCBI Handbook. 2003, Information NCfB. Bethesda (MD)Google Scholar
- Ashburner M, Ball CA, Blake JA: Gene ontology: Tool for the unification of biology. The Gene Ontology Consortium. Nat Genet. 2000, 25: 25-29. 10.1038/75556.View ArticlePubMedPubMed CentralGoogle Scholar
- The Gene Ontology Consortium: Creating the gene ontology resource: Design and implementation. Genome Res. 2001, 11: 1425-1433. 10.1101/gr.180801.View ArticleGoogle Scholar
- Conesa A, Götz S, García-Gómez JM, Terol J, Talón M, Robles M: Blast2GO: a universal tool for annotation, visualization and analysis in functional genomics research. Bioinformatics. 2005, 21: 3674-3676. 10.1093/bioinformatics/bti610.View ArticlePubMedGoogle Scholar
- Finn RD, Mistry J, Tate J, Coggill P, Heger A, Pollington JE, Gavin OL, Gunasekaran P, Ceric G, Forslund K, Holm L, Sonnhammer EL, Eddy SR, Bateman A: The Pfam protein families database. Nucleic Acids Res. 2010, 38: D211-222. 10.1093/nar/gkp985.View ArticlePubMedGoogle Scholar
- Gray SM, Kelly S, Robles LJ: Rho signaling mediates cytoskeletal re-arrangements in octopus photoreceptors. Amer Malac Bull. 2008, 26: 19-26. 10.4003/006.026.0203.View ArticleGoogle Scholar
- Thompson JD, Higgins DG, Gibson TJ: CLUSTAL W: improving the sensitivity of progressive multiple sequence alignment through sequence weighting, position-specific gap penalties and weight matrix choice. Nucleic Acids Res. 1994, 22: 4673-4680. 10.1093/nar/22.22.4673.View ArticlePubMedPubMed CentralGoogle Scholar
- Goldman N, Yang Z: A codon-based model of nucleotide substitution for protein-coding DNA sequences. Mol Biol Evol. 1994, 11 (5): 725-736.PubMedGoogle 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.