Phylogenetic analysis of the vertebrate Excitatory/Neutral Amino Acid Transporter (SLC1/EAAT) family reveals lineage specific subfamilies
© Gesemann et al; licensee BioMed Central Ltd. 2010
Received: 13 January 2010
Accepted: 29 April 2010
Published: 29 April 2010
The composition and expression of vertebrate gene families is shaped by species specific gene loss in combination with a number of gene and genome duplication events (R1, R2 in all vertebrates, R3 in teleosts) and depends on the ecological and evolutionary context. In this study we analyzed the evolutionary history of the solute carrier 1 (SLC1) gene family. These genes are supposed to be under strong selective pressure (purifying selection) due to their important role in the timely removal of glutamate at the synapse.
In a genomic survey where we manually annotated and analyzing sequences from more than 300 SLC1 genes (from more than 40 vertebrate species), we found evidence for an interesting evolutionary history of this gene family. While human and mouse genomes contain 7 SLC1 genes, in prototheria, sauropsida, and amphibia genomes up to 9 and in actinopterygii up to 13 SLC1 genes are present. While some of the additional slc1 genes in ray-finned fishes originated from R3, the increased number of SLC1 genes in prototheria, sauropsida, and amphibia genomes originates from specific genes retained in these lineages.
Phylogenetic comparison and microsynteny analyses of the SLC1 genes indicate, that theria genomes evidently lost several SLC1 genes still present in the other lineage. The genes lost in theria group into two new subfamilies of the slc1 gene family which we named slc1a8/eaat6 and slc1a9/eaat7.
The phylogeny of the SLC1/EAAT gene family demonstrates how multiple genome reorganization and duplication events can influence the number of active genes. Inactivation and preservation of specific SLC1 genes led to the complete loss of two subfamilies in extant theria, while other vertebrates have retained at least one member of two newly identified SLC1 subfamilies.
Genomes of extant species are shaped by extensive gene loss and duplication events. In the radiation of vertebrate species, two whole genome duplication events at the base of the vertebrate lineage (approximately 500 million years (Mya) ago) [1–3] and a third round of genome duplication at the base of the teleost lineage (about 350 Mya ago [4–7]) are proposed to play a crucial role [1–3]. Duplicated genes are expected to be functionally redundant and therefore released of selective pressure, consequently leading to gene loss. Nevertheless many duplicated paralogs are retained in modern genomes, for instance an estimated 15 - 24% of duplicated paralogs are present in extant teleost genomes [6, 8–10]. The duplication-complementation-degeneration model seeks to explain the retention of duplicated genes by subfunctionalization, where the function of the essential ancestral gene is distributed to two genes, each fulfilling part of the ancestral gene's function due to regulatory mutations [11, 12]. Another most interesting possible event following gene duplication is neofunctionalization, in which one of the two paralogous genes is free to acquire a new function that differs from the ancestral gene, due to the essential (ancestral) function being carried out by the unchanged paralog. Such functionalization events may pave the way for speciation. An alternative driving force in speciation might be divergent resolution, where random paralog losses in two allopatric populations can lead to diversity [13–15]. While whole genome duplications are one way to increase the gene repertoire within a genome, tandem duplication events as well as lineage specific gene loss are other means to change the number of active genes (e.g. [16, 17].
In the context of lineage specific gene loss and duplication events, we studied members of the vertebrate slc1 gene family of neutral and excitatory amino acid transporters (EAATs), whose presence on glia cells or neurons is indispensable for precise and sustained synaptic activity [18–22]. EAATs are involved in the removal of glutamate from excitatory synapses, a process not only essential to ensure precise termination of synaptic transmission but also to avoid neurotoxic accumulation implicated in a number of diseases . Moreover, transport of glutamate is associated with an increased chloride conductance across the membrane, which hyperpolarizes cells . Based on these important functions high-affinity glutamate transporters in the nervous system are supposed to be under strong purifying selection (pressure to stay the same).
Gene Aliases in Different Nomenclatures.
In teleosts, given the evidence for a third round of genome duplication (R3), an increased number of SLC1 genes are expected in their genomes. We indeed identified up to 13 SLC1s in fish genomes. While some of these additional SLC1 genes clearly are the consequence of the teleost specific R3, others are evidently members of two additional subfamilies (now called SLC1A8/EAAT6, SLC1A9/EAAT7), that were lost in the lineage leading to human and mouse. Interestingly members of these two subfamilies can also be found in prototherian, sauropsidan, and amphibian genomes, suggesting that these genes are specifically retained in these lineages.
The zebrafish genome contains 13 slc1genes on separate chromosomal locations
As a basis to study the evolutionary history of a given gene family, we analyzed the abundance of slc1 genes in the zebrafish genome. While in the human and mouse 7 different members of the slc1 gene family have been described, we identified and annotated 13 slc1 family genes in the zebrafish genome. Cloning from whole embryo cDNA indicated that all of the 13 annotated slc1 sequences are indeed transcribed. Sequencing of the amplified cDNA fragments revealed no significant deviation from our predicted sequences, except for some sites of single nucleotide polymorphisms.
The fact that the slc1 gene family in zebrafish consists of nearly twice as many genes as present in genomes of human and mouse suggests that at least some of the zebrafish slc1 orthologs have originated from the teleost specific whole-genome duplication. Physical and virtual mapping indicated that none of the sequences within a subgroup of slc1 are located within the same chromosomal cluster (Additional File 1), supporting the hypothesis that the genes originated from whole genome duplications rather than from individual tandem duplications.
Gene duplicates originating from whole genome duplications should be parts of large blocks of duplicated gene pairs, called paralogons. This holds true for slc1a3a/eaat1a and slc1a3b/eaat1b, which are located on chromosome 10 and 5 respectively. Moreover, these chromosomes are derived from a common protochromosome, further supporting a whole genome duplication event [7, 8].
Evidence for the existence of additional SLC1 subgroups
SLC1 gene retention and loss across vertebrate lineages
Database analyses indicated that all therian species have lost SLC1A8 and SLC1A9. In contrast, the members of egg laying vertebrates (prototheria, amphibia, and sauropsida) display various gene losses and retentions (Fig. 5B). Monotremates (prototheria) and birds have retained SLC1A8 and lost SLC1A9, whereas the reverse pattern holds true for lepidosauromorphes (lizards) that have retained SLC1A9 and lost SLC1A8. Xenopus tropicalis, as a representative of the amphibians on the other hand has retained members of both subfamilies. Interestingly, additional gene losses occurred in a variety of lineages. For instance, chicken (Gallus gallus) has lost SLC1A5 and the green anole lizard (Anolis carolinensis) has lost SLC1A6 (Fig. 5B).
Phylogenetic reconstruction including now SLC1 sequences from teleosts, amphibia, sauropsida and mammalia clearly indicates that the newly assigned SLC1A8 and SLC1A9 each form a separate clade (Fig. 5C).
In order to get an additional glimpse at the evolutionary history of SLC1 genes, we analyzed the genomes of two basal vertebrates, namely the elephant shark (Callorhinchus milii) and the sea lamprey (Petromyzon marinus). The elephant shark, belonging to the cartilaginous fish, represents an independent evolutionary branch of the vertebrates, whose genomes are thought to have undergone the same number of genome duplications as the mammalian branch . The low (1.4×) coverage of the elephant shark genome prevented the assembly of complete slc1 cDNA sequences. Nevertheless, our comparison revealed only exons for one ortholog per subfamily, with the notable exception of the SLC1A2/A9 subgroup, where we detected a retained SLC1A9 ortholog. Interestingly we also identified two paralogs of the SLC1A2/A9 subgroup in the lamprey, a jawless vertebrate at the base of the vertebrate tree (for a phylogenetic tree see Additional File 6).
Chromosomal location of putative human SLC1A8 and SLC1A9remnants
Furthermore synteny analysis in Xenopus and chicken also revealed ELAVL1 in close proximity to SLC1A8 within the same chromosomal fragment. Similarly we found the NUCB1 gene in the vicinity of SLC1A9 in the Xenopus and Anolis genome (data not shown).
The expression patterns of zebrafish slc1a3/eaat1 and slc1a2/eaat2paralogs suggest subfunctionalization
As an example to compare the spatial and temporal expression pattern of duplicated genes, we investigated the expression patterns of the zebrafish slc1a3 and slc1a2 subfamily genes.
We found non-overlapping expression for the slc1a3b/eaat1b paralog in neuromasts of the lateral line organ (Fig. 7C and 7D). Neuromasts consist of sensory hair cells and support cells that are homologous to mammalian inner ear hair cells and their adjacent support cells. The observed expression is consistent with a comparable function of this transporter in neuromasts, as observed for the ribbon synapse of the mammalian inner ear .
The slc1a2/eaat2 paralogs are both expressed in glia cells of the larval zebrafish brain (Fig. 7E-H) in a similar pattern to slc1a3a/eaat1a, which is consistent with SLC1A2/EAAT2 expression in mammalian astrocytes. In the retina, however, we found a complementary expression pattern with slc1a2a/eaat2a transcripts in the inner (Fig. 7E and 7F, arrow) and slc1a2b/eaat2b transcripts in the outer retina (Fig. 7H, arrow).
While the expression of zebrafish slc1a2b/eaat2b paralogs is in line with mammalian slc1a2/eaat2 expression in the outer retina, the slc1a2a/eaat2a expression in the inner retina has not been reported in mammals, suggestive of neofunctionalization.
Taken together our expression data for the slc1a3/eaat1 and slc1a2/eaat2 paralogs reveal subfunctionalization events in the zebrafish, since each paralog fulfills presumably parts of the ancestral function as predicted by the Duplication Degeneration-Complementation model [11, 12].
Modern genomes are shaped by frequent duplication and deletion events. Most dramatic are whole genome duplications that were prominently stated by Susumo Ohno as the most important driving force in the evolution of metazoans . This idea has immediate appeal since it provides an elegant route to the generation of new functions without compromising essential ancestral gene functions. Following the split from invertebrates, modern vertebrate genomes harbour remnants of two whole genome duplications [2, 3], and there is now overwhelming evidence that ray-finned fishes (teleosts) had an additional whole genome duplication dating back about 320-350 million years ago (Mya) , leading to many unique teleost specific genes. While analyzing solute carrier family 1 (slc1) genes in zebrafish, we found as expected evidence for many duplicated genes. However, some zebrafish slc1 genes turned out to be not the result of the teleost specific R3. Further analyses ruled out the possibility that these SLC1 genes originated from the teleost specific genome duplication. One possible explanation for these additional genes is lineage and species specific gene loss. This has long been recognized in unicellular organisms [27, 28]. However more recent reports demonstrate that such events are also observed in vertebrate species [16, 29–31], where the impact of such species specific gene loss remains unclear. In the case of the SLC1 family, where five of the seven human and mouse SLC1 genes are grouped in the functionally defined excitatory amino acid transporter (EAAT) family, in teleosts we found members of two additional subfamilies, now called SLC1A8/EAAT6 and SLC1A9/EAAT7, in teleosts. These two subfamilies have been lost in all therian mammals. This is quite surprising, as these proteins are involved in the removal of glutamate from vertebrate synapses and are therefore likely under immediate selective pressure.
In order to reconstruct the phylogenetic history of this gene family we studied a number of vertebrate species from all major vertebrate lineages. We found representatives of both new subfamilies in all analyzed teleost species, showing that the existence of these genes is common to teleost species and is not restricted to zebrafish. Analysis of amphibian genomes reveal the existence of one member of each SLC1A8/EAAT6 and SLC1A9/EAAT7 subfamily, arguing that the tetrapod ancestral species still contained the full complement of SLC1 subfamilies. In this respect it is interesting to note that a study of glutamate transporters in the tiger salamander Ambystoma tigrinum retina reports the cloning of a SLC1A8/EAAT6, as well as a SLC1A9/EAAT7 gene (referred to as EAAT2b and EAAT5b in the report, respectively) [32, 33], which indeed group into the two new subfamilies (data not shown). These findings provide additional evidence for a scenario where the genes of these subfamilies have been lost after the amphibia - amnionta split.
A clue to when the first subfamily was lost is provided by the sauropsida. Here the Lepidsauromorpha (represented by the green anole lizard Anolis carolinensis) have retained SLC1A9 but lost SLC1A8. The situation is reversed in the archosauromorpha (represented by the bird species zebrafinch and chicken), where SLC1A9 was lost but SLC1A8 was retained. This suggests that the genome of the last common ancestor of the sauropsida still contained both subfamilies, followed by independent gene loss in the two branches of the sauropsida. Following this logic the genome of the amniote ancestors should have contained members of both subfamilies.
In mammals only the protherian species (represented by the duck-billed Platypus Ornithorhynchus anatinus) have retained the member of SLC1A8, while all other remaining mammals have lost members of both subfamilies. In the most parsimonious scenario the last common mammalian ancestor has already lost SLC1A9 but has retained SLC1A8. Subsequently SLC1A8 was lost at the base of the therian lineage, whereas a member of this subfamily was retained in the branched off protheria species. Hence the SLC1A8 and SLC1A9 subfamilies have been lost independently in a number of lineages. While for instance SLC1A8 has been lost at least twice independently in the archosauromorpha and therian lineage, SLC1A9 has been lost independently in the therian and lepidosauromorpha lineage. In order to infer the ancestral situation we assessed the SLC1/EAAT repository in species that are at the base of the gnathosomes.
Genome information for two species of interest, the elephant shark (Callorhinchus milii) and the sea lamprey (Petromyzon marinus), is available. While the former represents the cartilaginous fish, an independent branch of the vertebrates that has split before the division of the tetrapod and the teleost lineage, the later belongs to a clade of skulled chordate animals lacking jaws. Since both lineages branched off before R3, cartilaginous fish as well as lampreys are expected to have undergone the same number of genome duplications as the mammalian branch. Consequently the gene repertoire between lamprey, sharks, and mammals should be comparable [25, 34]. While we found for each therian SLC1 gene a corresponding elephant shark ortholog, we identified several exons of a retained SLC1A9 ortholog, demonstrating that this basal species have retained members of additional SLC1 gene families. Interestingly we found only one exon sequence for a putative SLC1A8/EAAT6 ortholog in the lamprey and no trace of this gene in the elephant shark. This suggests that either this gene has been lost independently in both of these basal species (assuming that the putative lamprey exon is a false positive one) or that limited genomic information prevented us from identifying these corresponding genes.
Up to now, we have focused our discussion mainly on the SLC1A8 and SLC1A9 subfamilies. An interesting case is also the loss of SLC1A6/EAAT4 in the green anole lizard. This gene is not duplicated in any species analyzed and is prominently expressed in the cerebellum. We can only speculate that another SLC1 gene in lizards must have taken over the function of SLC1A6 in glutamate removal from cerebellar synapses. Support for such a scenario comes from a study comparing the similarity of the glutamatergic sytem in the turtle Chrysemys picta and rodents .
We have documented multiple gene loss in vertebrate lineages, raising the question if traces of these lost genes can still be found in extant genomes. Most of these gene losses must have happened many Mya ago, so finding defunct gene sequences is rather unlikely. However by microsynteny analyses, we found the putative genomic location of both lost SLC1 genes in the human genome. Intriguingly both lost genes were most likely located on human chromosome 19. While the putative SLC1A8 location is on the short arm of chromosome 19 (19p13.2/3), the most likely position of the inactivated SLC1A9 is on the long arm at position 19q13.3. For this human chromosome extensive rearrangements within the short arm and many intrachromosomal breaks in the long arm have been reported, even since the time of primate rodent divergence . Therefore it seems well conceivable that many genes, including SLC1A8 and SLC1A9, may have been rendered non-functional during these rearrangement events. Our microsynteny analyses indicated that the ELAVL1 gene is in close proximity of both zebrafish SLC1A8 genes, and that other genes in this region map like ELAVL1 to the human chromosome 19. Interestingly this synteny is conserved in chicken and Xenopus, where ELAVL1 and SLC1A8 are neighbors on chromosome 28 (890K-930K) and scaffold 112 (300K-400K) respectively (data not shown).
The phylogenetic analysis of a gene family spanning multiple taxa can yield important information about their evolutionary history. In the case of the solute carrier family 1, consistent with the duplication-complementation-degeneration model, most of the 7 human or mouse orthologs have one or two teleost orthologs. However, for two teleost SLC1 subfamilies we found no corresponding therian counterpart. This indicates that therians must have lost these two subfamilies, now named SLC1A8/EAAT6 and SLC1A9/EAAT7.
These two subfamilies are not teleost specific, but can be found throughout different vertebrates. An analysis of the major vertebrate lineages revealed an intriguing pattern of lineage specific gene losses, shaping the phylogenetic history of SLC1 genes. In all non-therian vertebrate species we found at least one member of the SLC1A8/9 subfamilies. Amphibians have neither lost SLC1A8 nor SLC1A9, lepidosauria (extent lizards, snakes and tuatara) on the one hand have lost the SLC1A8 ortholog, while birds on the other hand have lost SLC1A9. Interestingly egg laying mammals (platypus) are the only mammals that have also retained a SLC1A8 ortholog.
By studying one gene family across a number of vertebrate taxa, inferences can be made about the evolutionary history of these genes. Such an analysis is needed to guide the interpretation of functional data obtained for those genes and provides a fascinating scenario to study the evolution of vertebrate gene families.
Zebrafish maintenance and breeding
Wild-type fish of the Tübingen (Tü) inbred strain were bred and crossed as previously described . Embryos were raised in E3 medium (5 mM NaCl, 0.17 mM KCl, 0.33 mM CaCl2, 0.33 mM MgSO4) at 28°C. Stages refer to the development in days post fertilization (dpf).
Annotation of zebrafish SLC1 cDNAs
As many genes predicted within GenBank are produced by automated processes and have been shown to contain numerous errors such as missing exons and wrong predictions of start and termination sites, all cDNA sequences for the described slc1 genes were annotated manually. Sequences from the species listed in Additional File 3 were used for our analysis. Sequences were identified and annotated using combined information from EST, genome (GeneBank, http://www.ncbi.nlm.nih.gov; Ensembl, http://www.ensembl.org/index.html; version 50/51, 2008) and WGS http://blast.ncbi.nlm.nih.gov/Blast.cgi?PAGE=Nucleotides&PROGRAM=blastn&BLAST_SPEC=TraceArchive&BLAST_PROGRAMS=megaBlast&PAGE_TYPE=BlastSearch databases. Human and mouse sequences were used as initial query. Genomic regions coding for SLC1 family genes were identified by using the tblastx alignment algorithm (for details on intron/exon sizes see Additional File 1). Corresponding genomic fragments were analyzed using the GENSCAN gene prediction program http://genes.mit.edu/GENSCAN.html. The obtained putative cDNA and protein sequences were compared to the corresponding human and mouse orthologs and GENSCAN prediction errors were corrected by manual inspection of the intron/exon boundaries in false predicted regions (small exons and minor class splice sites are usually not recognized ). Gaps in the assembled sequences due to inaccurate or incomplete genome sequencing were wherever possible filled by corresponding EST or WGS sequences. Sequences found by EST and WGS searches were again blasted against the nr database and only if the top hits corresponded to the original query sequences were used in the annotation process. Note that in the case of duplicated genes annotation errors using WGS sequences could not be excluded due to the fact that individual exons could not be linked together. Sequence alignment of cDNA fragments and overlapping EST sequences was done using the SeqMan software (Lasergene, DNASTAR, Madison WI).
Phylogenetic tree analysis
Coding sequences of slc1 genes (except for the exon 234 analysis which was done on nucleotide sequences) were translated into proteins using the EditSeq software (Lasergene, DNASTAR, Madison WI) and obtained protein sequences were used to generate a combined sequence file in FASTA format (see Additional File 7 (protein FASTA) and Additional File 8 (exon 234 FASTA). Sequence alignment and phylogenetic analysis was performed on the Phylogeny.fr platform http://www.phylogeny.fr/version2_cgi/phylogeny.cgi. Sequences were aligned using MUSCLE (v3.7, Edgar 2004) configured for highest accuracy (MUSCLE with default settings). After alignment, ambiguous regions (i.e. containing gaps and/or are poorly aligned) were removed with Gblocks (v0.91b, ) using the following parameters: minimum length of a block after gap cleaning = 10; no gap positions were allowed in the final alignment; all segments with contiguous nonconserved positions bigger than 8 were rejected; minimum number of sequences for a flank position = 85%. The phylogenetic tree was reconstructed using the maximum likelihood method implemented in the PhyML program (v3.0 aLRT, . The default substitution model was selected assuming an estimated proportion of invariant sites and 4 gamma-distributed rate categories to account for rate heterogeneity across sites. The gamma shape parameter was estimated directly from the data. Reliability for internal branch was assessed using the bootstrapping method (100 bootstrap replicates) and bootstrap values above 50% (0.5) are shown. Graphical representation and editing of the phylogenetic tree was done using TreeDyn (v198.3,  and the obtained svg files were colored using the CorelDraw program. Animal icons in Figure 5 were purchased from iStockphoto http://www.istockphoto.com.
Alternative alignments and trees were calculated using the alignment programs ClustalW  and T-Coffee , and the neighbour joining (PHYLIP package 3.66 ; distances calculated using ProtDist  or MRBAYES (v3.1.2,  phylogeny programs. Using these alternative methods revealed some potential differences in phylogeny and a Bayesian calculated tree including human, mouse and zebrafish SLC1 sequences is shown in Additional File 2.
Approximately 50 5 dpf zebrafish larvae were homogenized to prepare total RNA using the Qiagen RNeasy kit (Qiagen), according to the manufacturer's instructions. Oligo dT primers were used to reverse-transcribe the total RNA with reverse transcriptase (First Strand Kit, Invitrogen, Carlsbad, CA). Polymerase chain reaction (PCR) was performed using sequence-specific oligonucleotide primers. Amplified fragments were ligated into the pCR-II vector (TA Cloning Kit Dual Promoter, Invitrogen, Carlsbad, CA, USA) and at least three independently amplified cDNA fragments per gene were sequenced to confirm our previously annotated zebrafish sequences. Accession numbers for the zebrafish SLC1/EAAT genes are: SLC1A1/EAAT3 (HM138690), SLC1A2a/EAAT2a (HM138691), SLC1A2b/EAAT2b (HM138692), SLC1A3a/EAAT1a (HM138693), SLC1A3b/EAAT1b (HM138694), SLC1A4 (HM138695), SLC1A5 (HM138696), SLC1A6/EAAT4 (HM138697), SLC1A7a/EAAT5a (HM138698), SLC1A7b/EAAT5b (HM138699), SLC1A8a/EAAT6a (HM138700), SLC1A8b/EAAT6b (HM138701), SLC1A9/EAAT7 (HM138702).
All zebrafish slc1 genes except slc1a6/eaat4 were localized in the zebrafish genome by radiation hybrid mapping  using the LN54 panel (Loeb/NIH/5000/4000, M. Ekker, Ottawa, Canada) and specific primer pairs for the genomic sequences (see Additional File 9).
For microsynteny analyses, genes flanking the slc1 locus of human or zebrafish were used as initial queries for a tblastx  search against the zebrafish or human database (ncbi nr/nt database) limited to the desired species (homo sapiens, danio rerio). The different hits were compared with each other and the hit with the highest conservation (identical aa and conserved aa) was selected (note that due to frequently occurring alignment gaps, hits with the lowest expect scores (E values) not necessarily representing the ones with highest identity). Only hits displaying a significant homology in length as well as conservation were used for further analyses. The hit with the highest conservation was used in a reciprocal tblastx search against the corresponding database and only genes which identified the initially used query are depicted.
Whole-mount in situbybridization
In vitro transcription of DIG-labeled probes was performed using the Roche RNA Labeling Kit (Roche Diagnostics, Rotkreuz, Switzerland). RNA probes were hydrolyzed to obtain fragments of about 300 bp length. Whole-mount in situ hybridization on 3 dpf zebrafish larvae was performed as follows. PTU treated larvae were anesthetized on ice and immediately fixed in 4% paraformaldehyde in 0.2 M phosphate buffer (pH 7.4) for 40 min at ambient temperature (RT). The in situ hybridization was performed according to the Zebrafish Book protocol  using an automated in situ hybridization machine (BioLane HTI, Hölle & Hüttner, Tübingen, Germany).
We thank A. Bohnert for the physical mapping of the zebrafish slc1 genes and Corinne Hodel and Dr. Edda Kastenhuber for critical reading of the manuscript. This work was supported by an ETH internal grant (TH-0-20150-3) (AL, SCFN) and by a grant of the Swiss National Science Foundation (PP00A-68868) (SCFN).
- Ohno S: Evolution by Gene Duplication. 1970, Springer Verlag, New YorkView ArticleGoogle Scholar
- Dehal P, Boore JL: Two rounds of whole genome duplication in the ancestral vertebrate. PLoS Biol. 2005, 3 (10): e314-10.1371/journal.pbio.0030314.PubMed CentralView ArticlePubMedGoogle Scholar
- Putnam NH, Butts T, Ferrier DE, Furlong RF, Hellsten U, Kawashima T, Robinson-Rechavi M, Shoguchi E, Terry A, Yu JK, et al: The amphioxus genome and the evolution of the chordate karyotype. Nature. 2008, 453 (7198): 1064-1071. 10.1038/nature06967.View ArticlePubMedGoogle Scholar
- Taylor JS, Braasch I, Frickey T, Meyer A, Peer Van de Y: Genome duplication, a trait shared by 22000 species of ray-finned fish. Genome Res. 2003, 13 (3): 382-390. 10.1101/gr.640303.PubMed CentralView ArticlePubMedGoogle Scholar
- Meyer A, Peer Van de Y: From 2R to 3R: evidence for a fish-specific genome duplication (FSGD). Bioessays. 2005, 27 (9): 937-945. 10.1002/bies.20293.View ArticlePubMedGoogle Scholar
- Jaillon O, Aury JM, Brunet F, Petit JL, Stange-Thomann N, Mauceli E, Bouneau L, Fischer C, Ozouf-Costaz C, Bernot A, et al: Genome duplication in the teleost fish Tetraodon nigroviridis reveals the early vertebrate proto-karyotype. Nature. 2004, 431 (7011): 946-957. 10.1038/nature03025.View ArticlePubMedGoogle Scholar
- Kasahara M, Naruse K, Sasaki S, Nakatani Y, Qu W, Ahsan B, Yamada T, Nagayasu Y, Doi K, Kasai Y, et al: The medaka draft genome and insights into vertebrate genome evolution. Nature. 2007, 447 (7145): 714-719. 10.1038/nature05846.View ArticlePubMedGoogle Scholar
- Woods IG, Wilson C, Friedlander B, Chang P, Reyes DK, Nix R, Kelly PD, Chu F, Postlethwait JH, Talbot WS: The zebrafish gene map defines ancestral vertebrate chromosomes. Genome Res. 2005, 15 (9): 1307-1314. 10.1101/gr.4134305.PubMed CentralView ArticlePubMedGoogle Scholar
- Postlethwait JH, Woods IG, Ngo-Hazelett P, Yan YL, Kelly PD, Chu F, Huang H, Hill-Force A, Talbot WS: Zebrafish comparative genomics and the origins of vertebrate chromosomes. Genome Res. 2000, 10 (12): 1890-1902. 10.1101/gr.164800.View ArticlePubMedGoogle Scholar
- Brunet FG, Crollius HR, Paris M, Aury JM, Gibert P, Jaillon O, Laudet V, Robinson-Rechavi M: Gene loss and evolutionary rates following whole-genome duplication in teleost fishes. Mol Biol Evol. 2006, 23 (9): 1808-1816. 10.1093/molbev/msl049.View ArticlePubMedGoogle Scholar
- Force A, Lynch M, Pickett FB, Amores A, Yan YL, Postlethwait J: Preservation of duplicate genes by complementary, degenerative mutations. Genetics. 1999, 151 (4): 1531-1545.PubMed CentralPubMedGoogle Scholar
- Lynch M, Force A: The probability of duplicate gene preservation by subfunctionalization. Genetics. 2000, 154 (1): 459-473.PubMed CentralPubMedGoogle Scholar
- Lynch M, Conery JS: The evolutionary fate and consequences of duplicate genes. Science. 2000, 290 (5494): 1151-1155. 10.1126/science.290.5494.1151.View ArticlePubMedGoogle Scholar
- Taylor JS, Peer Van de Y, Meyer A: Genome duplication, divergent resolution and speciation. Trends Genet. 2001, 17 (6): 299-301. 10.1016/S0168-9525(01)02318-6.View ArticlePubMedGoogle Scholar
- Evans BJ: Genome evolution and speciation genetics of clawed frogs (Xenopus and Silurana). Front Biosci. 2008, 13: 4687-4706. 10.2741/3033.View ArticlePubMedGoogle Scholar
- Feiner N, Begemann G, Renz AJ, Meyer A, Kuraku S: The origin of bmp16, a novel Bmp2/4 relative, retained in teleost fish genomes. BMC Evol Biol. 2009, 9 (1): 277-10.1186/1471-2148-9-277.PubMed CentralView ArticlePubMedGoogle Scholar
- 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 (5-6): 765-773. 10.1387/ijdb.072533km.View ArticlePubMedGoogle Scholar
- Amara SG, Fontana AC: Excitatory amino acid transporters: keeping up with glutamate. Neurochem Int. 2002, 41 (5): 313-318. 10.1016/S0197-0186(02)00018-9.View ArticlePubMedGoogle Scholar
- Arriza JL, Fairman WA, Wadiche JI, Murdoch GH, Kavanaugh MP, Amara SG: Functional comparisons of three glutamate transporter subtypes cloned from human motor cortex. J Neurosci. 1994, 14 (9): 5559-5569.PubMedGoogle Scholar
- Grewer C, Rauen T: Electrogenic glutamate transporters in the CNS: molecular mechanism, pre-steady-state kinetics, and their impact on synaptic signaling. J Membr Biol. 2005, 203 (1): 1-20. 10.1007/s00232-004-0731-6.PubMed CentralView ArticlePubMedGoogle Scholar
- Huang YH, Bergles DE: Glutamate transporters bring competition to the synapse. Curr Opin Neurobiol. 2004, 14 (3): 346-352. 10.1016/j.conb.2004.05.007.View ArticlePubMedGoogle Scholar
- Kanai Y, Hediger MA: The glutamate/neutral amino acid transporter family SLC1: molecular, physiological and pharmacological aspects. Pflugers Arch. 2004, 447 (5): 469-479. 10.1007/s00424-003-1146-4.View ArticlePubMedGoogle Scholar
- Beart PM, O'Shea RD: Transporters for L-glutamate: an update on their molecular pharmacology and pathological involvement. Br J Pharmacol. 2007, 150 (1): 5-17. 10.1038/sj.bjp.0706949.PubMed CentralView ArticlePubMedGoogle Scholar
- Ryan RM, Vandenberg RJ: A channel in a transporter. Clin Exp Pharmacol Physiol. 2005, 32 (1-2): 1-6. 10.1111/j.1440-1681.2005.04164.x.View ArticlePubMedGoogle Scholar
- Venkatesh B, Kirkness EF, Loh YH, Halpern AL, Lee AP, Johnson J, Dandona N, Viswanathan LD, Tay A, Venter JC, et al: Survey sequencing and comparative analysis of the elephant shark (Callorhinchus milii) genome. PLoS Biol. 2007, 5 (4): e101-10.1371/journal.pbio.0050101.PubMed CentralView ArticlePubMedGoogle Scholar
- Glowatzki E, Cheng N, Hiel H, Yi E, Tanaka K, Ellis-Davies GC, Rothstein JD, Bergles DE: The glutamate-aspartate transporter GLAST mediates glutamate uptake at inner hair cell afferent synapses in the mammalian cochlea. J Neurosci. 2006, 26 (29): 7659-7664. 10.1523/JNEUROSCI.1545-06.2006.View ArticlePubMedGoogle Scholar
- Kawahara Y, Imanishi T: A genome-wide survey of changes in protein evolutionary rates across four closely related species of Saccharomyces sensu stricto group. BMC Evol Biol. 2007, 7: 9-PubMed CentralPubMedGoogle Scholar
- Aravind L, Watanabe H, Lipman DJ, Koonin EV: Lineage-specific loss and divergence of functionally linked genes in eukaryotes. Proc Natl Acad Sci USA. 2000, 97 (21): 11319-11324. 10.1073/pnas.200346997.PubMed CentralView ArticlePubMedGoogle Scholar
- Mauti O, Sadhu R, Gemayel J, Gesemann M, Stoeckli ET: Expression patterns of plexins and neuropilins are consistent with cooperative and separate functions during neural development. BMC Dev Biol. 2006, 6: 32-10.1186/1471-213X-6-32.PubMed CentralView ArticlePubMedGoogle Scholar
- Tian X, Pascal G, Fouchecourt S, Pontarotti P, Monget P: Gene birth, death, and divergence: the different scenarios of reproduction-related gene evolution. Biol Reprod. 2009, 80 (4): 616-621. 10.1095/biolreprod.108.073684.View ArticlePubMedGoogle Scholar
- Demuth JP, Hahn MW: The life and death of gene families. Bioessays. 2009, 31 (1): 29-39. 10.1002/bies.080085.View ArticlePubMedGoogle Scholar
- Eliasof S, Arriza JL, Leighton BH, Kavanaugh MP, Amara SG: Excitatory amino acid transporters of the salamander retina: identification, localization, and function. J Neurosci. 1998, 18 (2): 698-712.PubMedGoogle Scholar
- Eliasof S, Arriza JL, Leighton BH, Amara SG, Kavanaugh MP: Localization and function of five glutamate transporters cloned from the salamander retina. Vision Res. 1998, 38 (10): 1443-1454. 10.1016/S0042-6989(97)00452-5.View ArticlePubMedGoogle Scholar
- Kuraku S, Meyer A, Kuratani S: Timing of genome duplications relative to the origin of the vertebrates: did cyclostomes diverge before or after?. Mol Biol Evol. 2009, 26 (1): 47-59. 10.1093/molbev/msn222.View ArticlePubMedGoogle Scholar
- Keifer J, Carr MT: Immunocytochemical localization of glutamate receptor subunits in the brain stem and cerebellum of the turtle Chrysemys picta. J Comp Neurol. 2000, 427 (3): 455-468. 10.1002/1096-9861(20001120)427:3<455::AID-CNE11>3.0.CO;2-X.View ArticlePubMedGoogle Scholar
- Grimwood J, Gordon LA, Olsen A, Terry A, Schmutz J, Lamerdin J, Hellsten U, Goodstein D, Couronne O, Tran-Gyamfi M, et al: The DNA sequence and biology of human chromosome 19. Nature. 2004, 428 (6982): 529-535. 10.1038/nature02399.View ArticlePubMedGoogle Scholar
- Mullins MC, Hammerschmidt M, Haffter P, Nusslein-Volhard C: Large-scale mutagenesis in the zebrafish: in search of genes controlling development in a vertebrate. Curr Biol. 1994, 4 (3): 189-202. 10.1016/S0960-9822(00)00048-8.View ArticlePubMedGoogle Scholar
- Patel AA, Steitz JA: Splicing double: insights from the second spliceosome. Nat Rev Mol Cell Biol. 2003, 4 (12): 960-970. 10.1038/nrm1259.View ArticlePubMedGoogle Scholar
- Castresana J: Selection of conserved blocks from multiple alignments for their use in phylogenetic analysis. Mol Biol Evol. 2000, 17 (4): 540-552.View ArticlePubMedGoogle Scholar
- Guindon S, Gascuel O: A simple, fast, and accurate algorithm to estimate large phylogenies by maximum likelihood. Syst Biol. 2003, 52 (5): 696-704. 10.1080/10635150390235520.View ArticlePubMedGoogle Scholar
- Chevenet F, Brun C, Banuls AL, Jacq B, Christen R: TreeDyn: towards dynamic graphics and annotations for analyses of trees. BMC Bioinformatics. 2006, 7: 439-10.1186/1471-2105-7-439.PubMed CentralView ArticlePubMedGoogle 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 (22): 4673-4680. 10.1093/nar/22.22.4673.PubMed CentralView ArticlePubMedGoogle Scholar
- Notredame C, Higgins DG, Heringa J: T-Coffee: A novel method for fast and accurate multiple sequence alignment. J Mol Biol. 2000, 302 (1): 205-217. 10.1006/jmbi.2000.4042.View ArticlePubMedGoogle Scholar
- Felsenstein J: Phylogeny interface package. Cladistics. 1989, 5: 164-166.Google Scholar
- Elias I, Lagergren J: Fast computation of distance estimators. BMC Bioinformatics. 2007, 8: 89-10.1186/1471-2105-8-89.PubMed CentralView ArticlePubMedGoogle Scholar
- Huelsenbeck JP, Ronquist F: MRBAYES: Bayesian inference of phylogenetic trees. Bioinformatics. 2001, 17 (8): 754-755. 10.1093/bioinformatics/17.8.754.View ArticlePubMedGoogle Scholar
- Walter MA, Spillett DJ, Thomas P, Weissenbach J, Goodfellow PN: A method for constructing radiation hybrid maps of whole genomes. Nat Genet. 1994, 7 (1): 22-28. 10.1038/ng0594-22.View ArticlePubMedGoogle Scholar
- Altschul SF, Gish W, Miller W, Myers EW, Lipman DJ: Basic local alignment search tool. J Mol Biol. 1990, 215 (3): 403-410.View ArticlePubMedGoogle Scholar
- Karlsson J, von Hofsten J, Olsson PE: Generating transparent zebrafish: a refined method to improve detection of gene expression during embryonic development. Mar Biotechnol (NY). 2001, 3 (6): 522-527. 10.1007/s1012601-0053-4.View ArticleGoogle Scholar
- Westerfield M: THE ZEBRAFISH BOOK: A guide for the laboratory use of zebrafish (Danio rerio). 2000, Eugene: University of Oregon PressGoogle 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.