Skip to main content
  • Research article
  • Open access
  • Published:

Molecular characterization and expression analysis of five different elongation factor 1 alpha genes in the flatfish Senegalese sole (Solea senegalensis Kaup): Differential gene expression and thyroid hormones dependence during metamorphosis

Abstract

Background

Eukaryotic elongation factor 1 alpha (eEF1A) is one of the four subunits composing eukaryotic translation elongation factor 1. It catalyzes the binding of aminoacyl-tRNA to the A-site of the ribosome in a GTP-dependent manner during protein synthesis, although it also seems to play a role in other non-translational processes. Currently, little information is still available about its expression profile and regulation during flatfish metamorphosis. With regard to this, Senegalese sole (Solea senegalensis) is a commercially important flatfish in which eEF1A gene remains to be characterized.

Results

The development of large-scale genomics of Senegalese sole has facilitated the identification of five different eEF1A genes, referred to as SseEF1A1, SseEF1A2, SseEF1A3, SseEF1A4, and Sse42Sp50. Main characteristics and sequence identities with other fish and mammalian eEF1As are described. Phylogenetic and tissue expression analyses allowed for the identification of SseEF1A1 and SseEF1A2 as the Senegalese sole counterparts of mammalian eEF1A1 and eEF1A2, respectively, and of Sse42Sp50 as the ortholog of Xenopus laevis and teleost 42Sp50 gene. The other two elongation factors, SseEF1A3 and SseEF1A4, represent novel genes that are mainly expressed in gills and skin. The expression profile of the five genes was also studied during larval development, revealing different behaviours. To study the possible regulation of SseEF1A gene expressions by thyroid hormones (THs), larvae were exposed to the goitrogen thiourea (TU). TU-treated larvae exhibited lower SseEF1A4 mRNA levels than untreated controls at both 11 and 15 days after treatment, whereas transcripts of the other four genes remained relatively unchanged. Moreover, addition of exogenous T4 hormone to TU-treated larvae increased significantly the steady-state levels of SseEF1A4 with respect to untreated controls, demonstrating that its expression is up-regulated by THs.

Conclusion

We have identified five different eEF1A genes in the Senegalese sole, referred to as SseEF1A1, SseEF1A2, SseEF1A3, SseEF1A4, and Sse42Sp50. The five genes exhibit different expression patterns in tissues and during larval development. TU and T4 treatments demonstrate that SseEF1A4 is up-regulated by THs, suggesting a role in the translational regulation of the factors involved in the dramatic changes that occurs during Senegalese sole metamorphosis.

Background

G proteins constitute a large superfamily of regulatory proteins that show a high similarity in primary structure and organization of the functional center [1–3]. They all share a common GTPase cycle, being turned on by binding GTP and off by hydrolyzing GTP to GDP. Additional regulatory proteins are often required to induce the conformational changes that occur during this cycle: guanine nucleotide exchange factors, which catalyze release of bound GDP and promote its replacement by GTP, and GTPase-activating proteins, which accelerate GTP hydrolysis. Such transition between active and inactive forms allows for them to serve as molecular switches and to make G proteins suitable for the regulation of a wide range of cellular processes such as signal transduction, cytoskeletal reorganizations, vesicular transport and protein synthesis [4]. The superfamily of G proteins includes three main classes: Ras-like GTPases, Gα subunits of heterotrimeric G proteins, and the translation elongation factors.

The eukaryotic translation elongation factor 1 alpha, currently termed eEF1A, is a member of the G protein family, and one of the four subunits that compose the eukaryotic elongation factor 1 [5, 6]. It is the second-most abundant protein after actin, comprising 1–3% of the total protein content in normal growing cells [7, 8]. In its GTP-liganded form, eEF1A catalyzes the first step of the elongation cycle during peptide synthesis, interacting with aminoacyl-tRNA to bring it to the acceptor site of the ribosome [6, 9]. Following hydrolysis of the GTP, the complex eEF1A·GDP is released from the ribosome and further converted to the active GTP-bound form by elongation factor 1 beta, which acts as a guanine nucleotide exchange factor for eEF1A. Apart from its key role in protein translation, other functions have been reported as well. eEF1A has been characterized as an actin binding protein with the ability to bundle F actin [10, 11] displaying microtubule severing [12] and stabilizing [13] activities, although the biological significance of these interactions is still unclear. In addition, eEF1A is thought to play a role in mediating signal transduction, as it has been shown to interact with PLCγ [14], and to be an activator of phosphoinositol 4-kinase [11, 15], an enzyme belonging to the signal transduction cascade activated by growth factors. Furthermore, an altered expression of eEF1A has been found in several kinds of tumors [16–20]. It has been also demonstrated that eEF1A is required for the degradation of Nα-acetylated proteins probably acting as a ubiquitin C-terminal hydrolase, thus rendering the substrate more easily degraded by the proteasome or by reducing the correct folding of the protein/ubiquitins complex for the subsequent proteolytic degradation [21, 22]. Finally, upregulation of eEF1A mediated either by p53 [23, 24] or in response to oxidative stress [25, 26] seems to be involved in apoptosis. It is likely that this broad diversity of functions may explain why eEF1A is such a well-conserved protein in eukaryotic organisms.

Eukaryotes possess a variable number of eEF1A genes. In plants from 10 to 15 eEF1A genes have been described in maize [27], and nine in cotton [28]. In fungi, the yeast Saccharomyces cerevisiae contains two eEF1A genes [29, 30], and Mucor racemosus has three [31]. In animals, the fruit fly Drosophila melanogaster contains two genes, one of which is found to be expressed only in certain stages of development [32]. Three eEF1A genes in Xenopus laevis are expressed at different developmental stages: the oocyte form (eEF-1αO), present in male and female germ cells, the somatic form (eEF-1αS), hardly detectable in oocytes but present in embryos and adult cells, and the 42Sp50 form, detected only in oocytes [33–35]. In the mammalian genome many eEF1A-like loci exist, most of which appear to be pseudogenes [36–38]. In fact, only two of these sequences are currently known to be actively expressed in human [39, 40], mouse [41, 42], rat [43–46], and rabbit [47, 48]. The two eEF1A protein isoforms encoded by these genes are currently termed eEF1A1 and eEF1A2, and display a quite different expression pattern in tissues. Whereas eEF1A1 is expressed in all tissues, eEF1A2 is only present in brain, heart, and skeletal muscle, which are tissues composed of cells locked in a state of nonproliferation [40, 45, 49, 50]. In teleosts, only one actively transcribed eEF1A gene has been reported in Danio rerio [51], Sparus aurata [52], Oreochromis niloticus [53], and Oryzias latipes [54]. In flatfish, partial sequences of one expressed eEF1A gene have been described in Scophthalmus maximus [55], Paralichthys olivaceus [56], and Paralichthys lethostigma [57].

Senegalese sole, Solea senegalensis (Pleuronectiformes: Soleidae), is a commercially important flatfish. This species undergoes metamorphosis from 12 to 19 days post hatch (DPH) during larval development. The process involves jaw and head restructuring, eye migration from the left to the right side, and a change from a symmetrical to an asymmetrical body shape [58]. These drastic morphological changes have been shown to be regulated by thyroid hormones (THs) in flatfish [59, 60]. Nevertheless, few data are available about the biological role of eEF1A in flatfish metamorphosis. In fact, only in Scophthalmus maximus this question has been somehow addressed. In this species, coincident T4 and eEF1A expression peaks are detected at the climax of metamorphosis [55]. Nevertheless, the gene expression pattern of Senegalese sole eEF1A during larval development and in tissues, and its dependence on THs, remains to be elucidated.

Large-scale genomics of Senegalese sole has allowed for the availability of a high number of EST sequences. In this study, we have identified and characterized five different eEF1A genes, referred to as SseEF1A1, SseEF1A2, SseEF1A3, SseEF1A4, and Sse42Sp50. The main sequence features are described hereinafter. A phylogenetic analysis was carried out to identify putative ortholog sequences. Gene expression profiles during larval development and in tissues from juvenile soles were explored using real-time PCR. Additionally, thiourea (TU) and T4 treatments were carried out in order to reveal the possible dependence of all SseEF1A expressions on THs.

Results

Molecular characterization of Senegalese sole eEF1A genes

Five Senegalese sole eEF1A genes, referred to as SseEF1A1, SseEF1A2, SseEF1A3, SseEF1A4, and Sse42Sp50, were identified after EST analysis of a normalized cDNA library constructed from different larval stages (pre-, meta- and post-metamorphosis), undifferentiated gonads, and six adult tissues (testis, ovary, stomach, intestine, liver and brain). A total of 69, 10, 3, 1 and 4 clones were identified as SseEF1A1, SseEF1A2, SseEF1A3, SseEF1A4 and Sse42Sp50, respectively [DDBJ:AB326302 to AB326306]. Sequencing of specific PCR products amplified from a pre-metamorphic library using Senegalese sole eEF1A-specific and the universal T3 and T7 primers confirmed sequences for all SseEF1A loci.

Main features of the cDNAs are summarized in Table 1. The cDNA lengths ranged between 1,523 and 1,743 nucleotides (nt) for SseEF1A3 and SseEF1A1, respectively. All of them contained a short 5'-untranslated region (between 31–73 nt) followed by an open reading frame of 1,365 nt (455 codons; positions 74–1,438 for Sse42Sp50), 1,386 nt (462 codons; positions 48–1,433, 35–1,420, and 32–1,417 for SseEF1A1, SseEF1A3, and SseEF1A4, respectively), and 1,389 nt (463 codons; positions 40–1,428 for SseEF1A2). The 3'-untranslated region varied between 103 (SseEF1A3) and 310 (SseEF1A1) nt long and included one (SseEF1A2, SseEF1A4 and Sse42Sp50) or two (SseEF1A1 and SseEF1A3) canonical polyadenylation signals (AATAAA), and a short oligo-A tail. The nucleotide sequences surrounding the ATG initiator codon fit well to the consensus sequence of Kozak [61], including A (SseEF1A1, SseEF1A2, SseEF1A3 and SseEF1A4) or G (Sse42Sp50) at position -3 and G at position +4 (considering position +1 the A from the ATG start codon). The predicted molecular weights and isoelectric points of the translated amino acid sequences of Senegalese sole eEF1A genes ranged between 49.86-50.53 kDa and 9.13-9.69, respectively.

Table 1 Main features of the cDNAs encoding for SseEF1As.

All SseEF1As possessed the characteristic regions G1 to G4 critical in GDP/GTP exchange, GTP-induced conformational change and GTP hydrolysis (Figure 1) [1]. The three consensus elements GXXXXGK (G14-K20), DXXG (D91-G94) and NKXD (N153-D156) [2] were identified in the GTP-binding domain of the Senegalese sole polypeptides. In addition, the predicted amino acid sequences were scanned against the PROSITE database of sequence motifs. The conserved motif DKLKAERERGITIDI(A/S) was identified at positions 61 to 76 in each of SseEF1As as the GTP-binding elongation factor signature (Figure 1).

Figure 1
figure 1

Comparison of the primary structure of SseEF1As and various others eEF1A proteins (see Table 3). The alignment was performed using MegAlign software. Dots represent identity with SseEF1A1, and dashes represent gaps. G1 to G4 indicate the critical regions involved in GDP/GTP exchange and GTP hydrolysis. The consensus sequence composed of the three consensus elements GXXXXGK (G14-K20), DXXG (D91-G94), and NKXD (N153-D156) present in the GTP-binding domain is shaded in grey. The GTP-binding elongation factor signature corresponding to amino acids 61 to 76 is boxed.

SseEF1A genes shared high sequence homology at nucleotide level in the coding region (62.9–88.6% identity) and at amino acid level (66.7–89.8% identity). A multiple alignment was performed with the DNA (Table 2; below diagonal) and amino acid sequences (Table 2; above diagonal) of eEF1A from other vertebrates (see Table 3). As observed for Senegalese sole sequences, the identities obtained in all the pairwise comparisons were lower at DNA than at amino acid level, indicating that the highest sequence variation was at the third codon positions. Identities of SseEF1A1 and SseEF1A2 amino acid sequences were higher than those determined for SseEF1A3 and SseEF1A4 with regard to eEF1A from the other species. By the other side, 42Sp50 sequences exhibited high identities, with values above 87% in the comparisons performed.

Table 2 Percentage of amino acid (above diagonal) and nucleotide (below diagonal) sequence identity among Senegalese sole elongation factors and several other eEF1As (see Table 3) as calculated using MegAlign.
Table 3 Vertebrate eEF1A sequences used in this study. The abbreviated name, description, accession number and species are indicated in each case.

The phylogenetic relationships of the five SseEF1As as well as the predicted protein sequence of other vertebrate eEF1A genes were determined using the maximum likelihood (ML) and neighbor-joining (NJ) methods. The ProtTest analysis [62] determined the JTT + G model as being the bestfit model of amino acid sequence evolution (-lnL = 8,642.63) with a gamma shape value (four rate categories) of 0.69. These settings were also employed for the NJ analysis. The topology of the trees obtained with these parameters are depicted in Figure 2 (ML) and Figure 3 (NJ), with Drosophila melanogaster [GeneBank:NM_079872) and Artemia [EMBL:X03349] eEF1A proteins used as outgroups. Complex evolutionary relatedness was obtained. SseEF1A1 was linked together with other teleost eEF1As in a consistent clade (63 and 56% of bootstrap value for ML and NJ, respectively). SseEF1A2 was clustered in a highly supported branch (92 and 79% of bootstrap value for ML and NJ, respectively) with turbot eEF1A as well as with two other teleost sequences from Oryzias latipes and Gasterosteus aculeatus. In addition, all 42Sp50 proteins were grouped together in a highly consistent ramification (100% bootstrap support in ML and NJ trees). Interestingly, both SseEF1A3 and SseEF1A4 appeared in a more basal position of both trees, and no consistent clustering with any of the sequences added was observed.

Figure 2
figure 2

Phylogenetic relationships among SseEF1As and a wide range of vertebrate eEF1As (see Table 3) using the maximum likelihood method. D. melanogaster and Artemia sp elongation factor 1 alpha proteins were used as outgroups to root tree. Only bootstrap values higher than 50% are indicated on each branch. The scale for branch length (0.1 substitutions/site) is shown below the tree.

Figure 3
figure 3

Phylogenetic relationships among SseEF1As and a wide range of vertebrate eEF1As (see Table 3) using the neighbor-joining method. D. melanogaster and Artemia sp elongation factor 1 alpha proteins were used as outgroups to root tree. Only bootstrap values higher than 50% are indicated on each branch. The scale for branch length (0.1 substitutions/site) is shown below the tree.

Expression levels of SseEF1A genes in tissues

Steady-state levels of the five SseEF1A transcripts were quantitated in liver, spleen, intestine, stomach, head-kidney, gills, skeletal muscle, brain, heart, and skin from juvenile soles (Figure 4A). Relative gene expression levels were normalized by measuring Ubiquitin gene and expressed relative to liver. All SseEF1A genes were found in detectable amounts in the tissues examined. SseEF1A1 transcripts were quite similar in all tissues examined except in muscle (60-fold lower than in liver; P < 0.05). SseEF1A2 reached the highest expression levels in skeletal muscle, heart, and brain (28, 17, and 9-fold higher than liver, respectively; P < 0.05). SseEF1A3 and SseEF1A4 showed very similar expression patterns, as they were strongly expressed in gills (750 and 13,000-fold higher than liver, respectively; P < 0.001) and skin (500 and 6,000-fold higher than liver, respectively; P < 0.001). Finally, Sse42Sp50 was expressed at a relatively high level in brain (32-fold higher than liver; P < 0.001).

Figure 4
figure 4

A) Relative expression levels in tissues of the five SseEF1A genes. Expression values were normalized to those of Ubiquitin. Data were expressed as the mean fold change (mean ± SEM, n = 3) from the calibrator group (liver). Values marked with an asterisk are significantly different from liver values (P < 0.05). L: liver, Sp: spleen, I: intestine, St: stomach, Hk: head-kidney, G: gills, M: skeletal muscle, B: brain, H: heart, Sk: skin. B) Comparison of the relative levels of SseEF1A transcripts in tissues. Data are expressed as the ratio (calculated using 2-(ΔCt)) of target mRNA to Ubiquitin mRNA.

Although they exhibited differential expression profiles, we calculated the relative amount of the five SseEF1A mRNA levels in the 10 tissues examined (Figure 4B). As a whole, SseEF1A1 transcripts were the most abundant with 60, 2,000, 62,000, and 35,000-fold higher overall mean expression ratios than SseEF1A2, SseEF1A3, SseEF1A4, and Sse42Sp50, respectively. Nevertheless, SseEF1A2 showed the highest values in muscle (20-fold higher than SseEF1A1), and it was only 4 and 3-fold lower expressed than SseEF1A1 in brain and heart, respectively. SseEF1A3 and SseEF1A4 reached relatively important expression levels in gills (7 and 3-fold lower than SseEF1A1, respectively) and skin (14 and 8-fold lower than SseEF1A1, respectively). Finally, Sse42Sp50 was expressed at very low amounts in all tissues.

Expression levels and regulation during larval development

Expression patterns of SseEF1A genes during larval development (from 2 to 22 DPH) were also determined. Data were normalized to the housekeeping gene glyceraldehyde-3-phosphate dehydrogenase (GAPDH2) [63], and further expressed relative to 2 DPH. All transcripts were detected very early at 2 DPH (Figure 5A). Nevertheless, they displayed different expression profiles during development. SseEF1A1 transcripts remained relatively constant with no significant changes between pre-metamorphosis and metamorphosis. A similar expression pattern was observed for SseEF1A3, although a significant peak in mRNA levels was observed at first feeding (1.89-fold; P < 0.05). The expression profile of SseEF1A2 and SseEF1A4 was quite similar. Both transcripts were constant until 15 DPH, when they increased significantly (7 and 22-fold, respectively; P < 0.001), and thereafter they rose progressively until the end of metamorphosis. In relation to Sse42Sp50, the highest expression levels were detected at first developmental stages with no significant changes from 6 to 22 DPH. As observed in tissues, SseEF1A1 was the most abundantly expressed of all genes during larval development (Figure 5B).

Figure 5
figure 5

A) Relative SseEF1A expression levels during larval development (from 2 to 22 DPH) in Senegalese sole. Expression values were normalized to those of GAPDH2. Data are expressed as the mean fold change (mean ± SEM, n = 3) from the calibrator group (2 DPH). Different letters denote days that are significantly different (P < 0.05) analyzed by ANOVA followed by a Tukey test. The interval for the metamorphic process is shaded. B) Comparison of the relative levels of SseEF1A transcripts during larval development. Data are expressed as the ratio (calculated using 2-(ΔCt)) of target mRNA to GAPDH2 mRNA.

To study the involvement of THs on the expression of SseEF1A genes, 7 DPH larvae were exposed to the goitrogen TU. As a consequence of the TU treatment, the metamorphic process was blocked at S1-S3 stages as determined by the degree of eye migration. No differences in survivability were observed with respect to the untreated control (not shown). mRNA levels for SseEF1A genes were quantified in whole larvae pools collected at 8 hours, and 6 days, 11 days, and 15 days after treatment (dat). Untreated control larvae exhibited expression profiles similar to those described above in all cases (Figure 6). No significant differences in gene expression were observed for SseEF1A1, SseEF1A2, SseEF1A3 and Sse42sp50 between untreated control and TU-treated larvae. However, TU-treated larvae showed 3 and 4-fold lower (P < 0.05) SseEF1A4 mRNA levels than untreated controls at both 11 and 15 dat, respectively.

Figure 6
figure 6

Relative SseEF1A expression levels as determined by real-time quantitative PCR in the untreated (grey) and TU-treated (black) groups. Untreated and TU-treated samples were collected for RNA isolation at four different times: 8 hours, and 6 days, 11 days, and 15 days after treatment starting on 7 DPH. Expression values were normalized to those of GAPDH2. Data are expressed as the mean fold change (mean ± SEM, n = 3) from the calibrator group (control 8 hours). Values marked with an asterisk are significantly different (P < 0.05 or less) from the corresponding untreated control group values.

At the sight of these results, we carried out a rescue assay to determine the ability of T4 to revert the TU effect on the expression of SseEF1A4. For this purpose, 4 DPH larvae were exposed to TU to ensure the lack of endogenously synthesized T4 at day 7, when the exogenous T4 (100 ppb) was supplied. mRNA levels of SseEF1A4 gene were determined in larvae sampled at 8 and 13 days after T4 treatment (datt). The metamorphosis blocking induced by TU was not observed in TU+T4 treated larvae. The expression patterns of untreated and TU-treated larvae were similar to those described above (Figures 5 and 6). Interestingly, TU+T4 treated larvae exhibited significantly higher SseEF1A4 transcript levels than untreated (10.5 and 2-fold; P < 0.05) and TU-treated (16.8 and 8.7-fold; P < 0.05) larvae both at 8 and 13 datt, respectively (Figure 7).

Figure 7
figure 7

Relative SseEF1A4 expression levels as determined by real-time quantitative PCR in the untreated control (grey), and in TU (black) and TU+T4 (white) treated groups. Larvae samples were collected for rRNA isolation at 8 and 13 days after T4 treatment starting 7 DPH. Expression values were normalized to those of GAPDH2. Data are expressed as the mean fold change (mean ± SEM, n = 3) from the calibrator group (untreated control day 8). Values marked with an asterisk are significantly different (P < 0.05 or less) from the corresponding untreated control group values.

Discussion

Three rounds of large-scale gene duplications, referred to as 1R, 2R, and 3R or fish-specific genome duplication, have been proposed to take place in the fish evolutionary lineage [64, 65]. As a result, several gene copies have been described in teleosts for some group of genes including glycolytic enzymes [63, 64], Hox genes [66, 67], hormones and their receptors [68, 69], and keratins [70]. The majority of gene duplicates have been lost or silenced during evolution but in some instances paralogue genes have picked up new functions (neofunctionalization), or divided the ancestral function between paralogs (subfunctionalization) [71, 72]. Such a selective advantage would have avoided loss of functionally active duplicates. In this study, we have obtained the complete cDNA sequence of five distinct Senegalese sole eEF1A genes, referred to as SseEF1A1, SseEF1A2, SseEF1A3, SseEF1A4, and Sse42Sp50. All of them share high identity at both nucleotide (62.9-88.6%) and amino acid (66.7-89.8%) levels. Amino acid sequence analysis of the predicted polypeptides revealed the presence of the consensus amino acid sequence that identifies a GTP-binding protein [1, 2] in all of them with spacing between the characteristic sequence elements GXXXXGK (G14-K20), DXXG (D91-G94) and NKXD (N153-D156) that fell into the expected size (40–80 residues). Furthermore, the five Senegalese sole proteins possessed the characteristic motif that provides a signature for GTP-binding elongation factor. Altogether, these results indicate that SseEF1A genes actually encode for functional eEF1A proteins.

Phylogenetic analyses were carried out using a wide range of vertebrate eEF1A predicted proteins. A number of sequences found in the complete genome of several teleosts and mammals belonging to eEF1A family were also included, although it is not known if such sequences are actively transcribed or represent pseudogenes. Phylogeny gives strong support to the classification of SseEF1A1 as the Senegalese sole ortholog of teleost eEF1A. Furthermore, we can consider that teleost eEF1A, including SseEF1A1, is the fish counterpart of tetrapod eEF1A1. Such a conclusion is based on two main observations. Firstly, the high sequence identities of teleost with tetrapod eEF1A proteins as obtained with the sequences employed in this study, ranging between 86.6–89.9% (not shown). And secondly, gene expression analyses indicate that teleost eEF1A [52–54] as well as SseEF1A is a highly abundant mRNA transcript that exhibits a relatively uniform distribution across many tissues, as observed in mammalian eEF1A1 [40, 45, 48, 49, 73]. In contrast, strikingly enough, turbot SmaEF1A as well as the predicted polypeptides of one O. latipes and G. aculeatus sequences appeared closely linked to SseEF1A2 in a different clade to teleost eEF1A, being even more closely related to mammalian (and avian) eEF1A isoforms. According to the expression pattern in tissues of SseEF1A2, we suggest that this fish paralog represents the counterpart to mammalian eEF1A2. In this sense, northern blot and immunoblotting analyses have revealed that expression of mammalian eEF1A2 is limited to the terminally differentiated cells of the skeletal muscle, brain, and heart [40, 45, 48, 49, 73], the same tissues where SseEF1A2 reached the highest expression levels. If this gene is only expressed in the species analyzed here or in a wide range of flatfish or teleosts is a problem that needs to be resolved in future studies. In any case, the presence of these two genes in fish and mammals suggests that such gene duplication was an early evolutionary event. An interesting issue is that, in contrast to teleost, a close phylogenetic relatedness between mammalian eEF1A1 and eEF1A2 isoforms was observed. One possible explanation for this important difference is the increased evolutionary rate that is known to occur in fish genes [74] making sequence divergence between mammalian and fish counterparts comparable to a pair of fish paralogous genes.

The biological roles of mammalian eEF1A isoforms have been widely studied [6, 8, 75]. At present, it is unclear why there are two isoforms of eEF1A. In vitro, the translational elongation activity is similar but there is a difference in the off-rate for GDP as eEF1A2 binds GDP more strongly than GTP, whereas the opposite is true for eEF1A1 [48]. Given that eEF1A2 is expressed specifically in cells that have permanently withdrawn from the cell cycle such as neurons, cardiomyocites and myocites [40, 41, 48], it has been suggested that eEF1A2 may be involved in pathways of protein synthesis specifically for nonproliferating cells or in the non-translational roles suggested for eEF1A. It seems that eEF1A2 is anti-apoptotic [76], which is closely related with its possible role in tumourigenesis [17, 20, 77]. If we consider the high degree of structural and functional relatedness existing between eEF1A from different organisms, it is highly probable that fish orthologous genes carry out the same functions as mammalian isoforms.

All 42Sp50 proteins grouped together in a highly consistent cluster both in ML and NJ trees. Although orthologous genes have been described in other teleosts, we could no identify a mammalian counterpart in any of the complete genomes available in databases. Thus, it seems that this gene was lost early in the evolution of tetrapods. Little information is found concerning its possible function. In X. laevis, 42Sp50 transcripts have only been detected in oocytes, and the protein interacts with other components of the 42S RNP particle, which serves as storage of 5S rRNA and various tRNAs [34]. Senegalese sole Sse42Sp50 was expressed in all tissues examined, but at a highly significant lesser amount than SseEF1A1. Phylogeny of SseEF1A3 and SseEF1A4 is more difficult to interpret because it was poorly resolved, mainly due to the fact that no teleost (and no mammalian) sequence significantly clustered with them. Furthermore, we have failed to identify possible counterparts in the complete genome of other teleosts and mammals available in databases. Nevertheless, we argue that both SseEF1A3 and SseEF1A4 are probably the result of recent gene duplications. Both genes expressed mainly in gills and skin, two tissues involved in respiratory and osmoregulatory processes [78]. It is tempting to assume that they could be involved in the translational regulation of factors playing a role in such important biological activities in the Senegalese sole. And in the case of SseEF1A4, an additional role in metamorphosis will be discussed below.

Metamorphosis in flatfish is characterized by the change from a symmetrical, pelagic larva to an asymmetric, benthic juvenile, and by the eye migration from one side of the head to the other. THs have shown to play a key role in flatfish metamorphosis. Important physiological and molecular switches in erythropoiesis, myosin light chains, troponin T, and gill mitochondria-rich cells as well as differentiation of gastric glands and epidermal adult cells, and re-organization of the neural retina skin have been reported to occur during this dramatic process in a THs dependent manner [79–85]. More recently, it has also been reported that THs down-regulate GAPDH1 transcripts after metamorphosis climax [63] demonstrating the dependence of the glycolytic pathway on THs. With regard to eEF1A an expression peak at the climax of metamorphosis, coincident with a T4 peak, has been detected in turbot [55]. Interestingly, regulation by THs of the germ cell specific EF-1αO gene has been demonstrated during metamorphosis in X. laevis [86]. The start of metamorphosis has been associated to a surge of THs that increase their levels until the metamorphic climax, and reduce towards post-climax [59, 60]. In this study, the possible regulation of SseEF1A genes by THs was analyzed using TU and combined TU+T4 treatments. TU is a blocking agent of THs synthesis that can reduce T4 levels by about 95% [82]. Exogenous treatments with TU as well as other types of thyroid inhibitors have proved useful to study the involvement of THs on metamorphosis in flatfish [59, 82, 84] and the transcriptional regulation of genes involved in the pituitary-thyroid axis [87, 88]. TU-treated larvae exhibited significantly lower SseEF1A4 transcripts levels than untreated larvae at late metamorphosis, whereas mRNA levels of the other four genes remained relatively unchanged. In contrast, addition of T4 increased the amount of SseEF1A4 transcripts to significantly higher levels than those detected in untreated and TU-treated larvae. These results demonstrate that THs up-regulate either directly or indirectly SseEF1A4 expression, and constitute the first report of the dependence on THs of an elongation factor gene during flatfish metamorphosis. Although further research will be required to elucidate the precise role of SseEF1A4 in the drastic changes that occur during metamorphosis, it can be hypothesized that SseEF1A4 is involved in the translational regulation of the factors that promote an appropriate transition from larval to juvenile stage in response to THs. In fact, elongation factors are directly involved in the regulation of peptide synthesis [89, 90]. Yet, other non-translational functions during metamorphosis cannot be ruled out.

Conclusion

In conclusion, in this work we describe the sequence and main features of five eEF1A genes in the Senegalese sole, referred to as SseEF1A1, SseEF1A2, SseEF1A3, SseEF1A4, and Sse42Sp50. Combined phylogenetic and tissue expression analyses revealed SseEF1A1 and SseEF1A2 as the putative counterparts of mammalian eEF1A1 and eEF1A2, respectively, and Sse42Sp50 as the ortholog of X. laevis and teleost 42Sp50. In contrast, SseEF1A3 and SseEF1A4 represent novel genes, as orthologous genes have not been identified neither in fish nor mammals. The five genes exhibited different expression patterns during larval development. TU and T4 treatments demonstrated that SseEF1A4 was up-regulated by THs, suggesting a role in the translational regulation of the factors involved in the dramatic changes that occur during Senegalese sole metamorphosis. The regulation and interaction of Senegalese sole elongation factors with other genes are questions that need to be clarified in the future and will help to improve the aquaculture conditions for this species. Wih regard to this, present work constitutes undoubtly a useful framework for further studies.

Methods

Source of fish and experimental rearing conditions

All experimental sole larvae were obtained from fertilized eggs collected from breeding tanks, where breeders spawned naturally under environmental conditions. Eggs were incubated at a density of 2,000 eggs L-1 in 300 L cylinder conical tanks with gentle aeration, and water exchange every two hours. Temperature and salinity during all experiments were 20°C and 38 psu, respectively. Newly hatched larvae were transferred to a 400 L tank at an initial density of 45 to 50 larvae L-1 with a 16L:8D photoperiod, and a light intensity of 600–800 lux. Larvae were fed rotifers (Brachionus plicatilis) 3 DPH until 9 DPH. Enriched Artemia metanauplii were fed to them from 7 DPH until the end of the experiment. Pools of larvae from 2 to 22 DPH (n = 3) were collected, washed with DEPC water, frozen in liquid nitrogen, and stored at -80°C until analysis. Metamorphic stages (S0-S4) were classified according to [58].

Chronic exposure to 30 ppm TU was achieved in 200 L round tanks. TU was added on day 7, and water was kept stagnant for 24 hours. After this exposure, 20% of the water was exchanged daily, with the subsequent addition of the eliminated TU. A second tank with the same characteristics was used to rear untreated control sole larvae. Larvae were initially stocked at a density of 100 individual L-1, and lights were kept off until the onset of external feeding on 3 DPH. Fluorescent lamps supplied an illumination of 800 lux on the water surface, while a 16L:8D light cycle was used. Larvae were fed rotifers (Brachionus plicatilis) and 48 hours Artemia metanauplii, according to the experimental design. The Haptophyceae Isochrysis galbana (T-ISO strain) was employed for the enrichment of both live preys, and was also added (2 mg dry weight L-1 d-1) to larval tanks during rotifer feed. Daily larval samples were removed for metamorphosis control. Pools of larvae (n = 3) were also collected, washed with DEPC water, frozen in liquid nitrogen, and stored at -80°C until analysis.

For the rescue experiments, culture conditions were the same as described above, but 30 ppm TU was added on day 4 to ensure lack of endogenously synthesized T4 at the commencement of exogenous T4 treatment. At day 7, untreated control larvae were transferred to a pair of 16 L cylinder conical tubes at an initial density of 45 larvae L-1. TU-treated larvae were transferred to 4 tubes, two of which were supplemented with 100 ppb T4. After this new exposure, 20% of the water was exchanged every two days, with the subsequent addition of the eliminated TU and TU+T4. Pools of larvae (n = 3) were collected 8 and 13 days after the commencement of T4 treatment, washed with DEPC water, frozen in liquid nitrogen, and stored at -80°C until analysis.

Juvenile Senegalese sole individuals (average weight = 169.60 ± 24.3 g; n = 3) were obtained from IFAPA Centro El Toruño facilities (El Puerto Santa María, Cádiz, Spain). They were sacrificed by immersion in tricaine methanesulfonate (MS-222) according to the guidelines on the care and use of fish in research, teaching and testing from the Canadian Council on Animal Care (2005). Head-kidney, liver, testis, brain, heart, skin, skeletal muscle, spleen, intestine, stomach, and gills were rapidly dissected, frozen in liquid nitrogen, and stored at -80°C until use.

Identification and sequence analyses of elongation factor 1 alpha cDNAs in Senegalese sole

Ten cDNA libraries were constructed from different larval stages and adult tissues of Senegalese sole using the ZAP Express® cDNA Syntesis kit and Zap Express cDNA Gigapack® III Gold Cloning kit (Stratagene) following the manufacturer's protocol (Cerdà et al., in preparation). Libraries were pooled and normalized, and approximately 11,000 randomly selected clones were sequenced from the 3'-end. Expressed sequence tags (ESTs) encoding eEF1A genes were identified after EST annotation. Sequences have been deposited with accession numbers [DDBJ:AB326302 to AB326306].

Alignments of the predicted polypeptide sequences were carried out, and the sequence identities calculated by the MegAlign program from the LASERGENE software suite. For this purpose as well as for the ML and NJ phylogenetic analyses, predicted polypeptides of a range of vertebrate eEF1A sequences retrieved from GenBank/EMBL/DDBJ were used (Table 3). Additional sequences belonging to the eEF1A family were retrieved from ENSEMBL (Table 3). Drosophila melanogaster [GeneBank:NM_079872) and Artemia [EMBL:X03349] elongation factor 1 alpha sequences were used as outgroups to root ML and NJ trees. The bestfit model of sequence evolution was determined by using the ProtTest v 1.3 [62]. Phylogenetic analysis was carried out using the PHYLIP package [91] as follows. A bootstrap analysis was performed using SEQBOOT (200 replicates for ML; 1000 replicates for NJ). Data were then analyzed by ML using the software PHYML that generated 200 trees. The consensus phylogenetic tree was then obtained (CONSENSE). For NJ, data were analyzed with PROTDIST that generated 1000 distance matrices. The program NEIGHBOR was then employed to generate 1000 trees. Finally, the consensus phylogenetic tree was obtained using CONSENSE. Trees were drawn using the TreeViewX program v 0.5.0. Presence of known domains in the predicted protein sequence of Senegalese sole eEF1As was determined by scanning in the PROSITE database [92].

RNA isolation and gene expression analysis

Homogenization of juvenile tissues and larvae was carried out using the Lysing Matrix D (Q-BioGene) for 40 s at speed setting 6 in the Fastprep FG120 instrument (Bio101). Total RNA was isolated from 50 mg of pooled larvae using the RNeasy Mini Kit (Qiagen). All RNA isolation procedures were performed in accordance with the manufacturer's protocol. In all cases, total RNA was treated twice with DNase I using the RNase-Free DNase kit (Qiagen) for 30 min in order to avoid amplification of genomic DNA. RNA sample quality was checked using Experion (Bio-Rad), and quantification was performed spectrophotometrically. Total RNA (1 μg) from each sample was reverse-transcribed using the iScript™ cDNA Synthesis kit (Bio-Rad). Lack of genomic DNA contamination was confirmed by PCR amplification of RNA samples in the absence of cDNA synthesis.

Real-time analysis was carried out on an iCycler (Bio-Rad). Matching oligonucleotide primers were designed using the Oligo v 6.89 software (Medprobe). To confer specificity, each Senegalese sole eEF1A primer pair (Table 4) was designed using as template the corresponding 3'-UTR, where the high sequence variability made impossible to correctly align them. In addition, amplification of specific SseEF1A encoding cDNA fragments was verified by direct sequencing of PCR products in a previous assay. Real-time reactions were accomplished in a 25 μl volume containing cDNA generated from 10 ng of original RNA template, 300 nM each of specific forward (F) and reverse (R) primers, and 12.5 μl of iQ™ SYBR Green Supermix (Bio-Rad). The amplification protocol used was as follows: initial 7 min denaturation and enzyme activation at 95°C, 40 cycles of 95°C for 15 s, and 70°C for 30 s. Each assay was done in duplicate. For normalization of cDNA loading, all samples were run in parallel using the housekeeping genes GAPDH2 [63] or Ubiquitin as loading standards for larval development and tissues, respectively (Table 4). To estimate efficiencies, a standard curve was generated for each primer pair based on known quantities of cDNA (10-fold serial dilutions corresponding to cDNA transcribed from 100 to 0.01 ng of total RNA). All calibration curves exhibited correlation coefficients higher than 0.99, and the corresponding real-time PCR efficiencies were in the range 0.92–0.98. Relative mRNA expression was determined using the 2-(ΔΔCt) method [93]. To compare the amount of each of the Senegalese sole transcripts, the relative mRNA expression was determined by subtracting the Ct value of GAPDH2 or Ubiquitin mRNA from the Ct value of the target mRNA. Data were then expressed as the ratio (calculated using 2-(ΔCt)) of target mRNA to reference (GAPDH2 or Ubiquitin) mRNA. Owing to the small differences among amplicons both in size (62–142) and %GC (52.8, 58.1, 52.4, 46.7 and 67.8% for SseEF1A1, SseEF1A2, SseEF1A3, SseEF1A4 and Sse42Sp50, respectively), the relative sensitivity factor KRS was assumed to be one [94].

Table 4 Primers used for real-time gene expression analysis. The amplicon size generated by each primer pair is shown.

Results were expressed as mean ± SEM. Comparisons among groups were carried out with one-way analysis of variance (ANOVA), followed by a Tukey test for identification of the statistically distinct groups. Significance was accepted for P < 0.05.

Abbreviations

dat:

days after TU treatment

datt:

days after T4 treatment

DPH:

days post hatch

ML:

maximum likelihood

NJ:

neighbor-joining

SEM:

standard error of the mean

THs:

thyroid hormones

TU:

thiourea

References

  1. Bourne HR, Sanders DA, McCormick F: The GTPase superfamily: conserved structure and molecular mechanism. Nature 1991, 349: 117-127. 10.1038/349117a0.

    Article  CAS  PubMed  Google Scholar 

  2. Dever TE, Glynias MJ, Merrick WC: GTP-binding domain: three consensus sequence elements with distinct spacing. Proc Natl Acad Sci USA 1987, 84: 1814-1818. 10.1073/pnas.84.7.1814.

    Article  CAS  PubMed  Google Scholar 

  3. Kubarenko AV, Sergiev PV, Rodnina MV: GTPases of translational apparatus. Mol Biol 2005, 39: 746-761. 10.1007/s11008-005-0080-2.

    Article  CAS  Google Scholar 

  4. Bourne HR, Sanders DA, McCormick F: The GTPase superfamily: a conserved switch for diverse cellular functions. Nature 1990, 348: 125-132. 10.1038/348125a0.

    Article  CAS  PubMed  Google Scholar 

  5. Browne GJ, Proud CG: Regulation of peptide-chain elongation in mammalian cells. Eur J Biochem 2002, 269: 5360-5368. 10.1046/j.1432-1033.2002.03290.x.

    Article  CAS  PubMed  Google Scholar 

  6. Ejiri S: Moonlighting functions of polypeptide elongation factor 1: from actin bundling to zinc finger protein R1-associated nuclear localization. Biosci Biotechnol Biochem 2002, 66: 1-21. 10.1271/bbb.66.1.

    Article  CAS  PubMed  Google Scholar 

  7. Slobin LI: The role of eucaryotic factor Tu in protein synthesis. The measurement of the elongation factor Tu content of rabbit reticulocytes and other mammalian cells by a sensitive radioimmunoassay. Eur J Biochem 1980, 110: 555-563. 10.1111/j.1432-1033.1980.tb04898.x.

    Article  CAS  PubMed  Google Scholar 

  8. Condeelis J: Elongation factor 1 alpha, translation and the cytoskeleton. Trends Biochem Sci 1995, 20: 169-170. 10.1016/S0968-0004(00)88998-7.

    Article  CAS  PubMed  Google Scholar 

  9. Moldave K: Eukaryotic protein synthesis. Annu Rev Biochem 1985, 54: 1109-1149. 10.1146/annurev.bi.54.070185.005333.

    Article  CAS  PubMed  Google Scholar 

  10. Yang F, Demma M, Warren V, Dharmawardhane S, Condeelis J: Identification of an actin-binding protein from Dictyostelium as elongation factor 1a. Nature 1990, 347: 494-496. 10.1038/347494a0.

    Article  CAS  PubMed  Google Scholar 

  11. Yang W, Boss WF: Regulation of phosphatidylinositol 4-kinase by the protein activator PIK-A49. Activation requires phosphorylation of PIK-A49. J Biol Chem 1994, 269: 3852-3857.

    CAS  PubMed  Google Scholar 

  12. Shiina N, Gotoh Y, Kubomura N, Iwamatsu A, Nishida E: Microtubule severing by elongation factor 1 alpha. Science 1994, 266: 282-285. 10.1126/science.7939665.

    Article  CAS  PubMed  Google Scholar 

  13. Durso NA, Cyr RJ: A calmodulin-sensitive interaction between microtubules and a higher plant homolog of elongation factor-1 alpha. Plant Cell 1994, 6: 893-905. 10.1105/tpc.6.6.893.

    Article  CAS  PubMed  Google Scholar 

  14. Kim MJ, Si F, Kim SJ, Hong SB, Hwang JI, Lee HJ, Lee SJ, Chang JS, Lee YH, Ryu SH, Suh PG: The SH2-SH2-SH3 domain of phospholipase C-gamma1 directly binds to translational elongation factor-1alpha. Mol Cells 1999, 9: 631-637.

    CAS  PubMed  Google Scholar 

  15. Yang W, Burkhart W, Cavallius J, Merrick WC, Boss WF: Purification and characterization of a phosphatidylinositol 4-kinase activator in carrot cells. J Biol Chem 1993, 268: 392-398.

    CAS  PubMed  Google Scholar 

  16. Adams SM, Sharp MG, Walker RA, Brammar WJ, Varley JM: Differential expression of translation-associated genes in benign and malignant human breast tumours. Br J Cancer 1992, 65: 65-71.

    Article  CAS  PubMed  Google Scholar 

  17. Anand N, Murthy S, Amann G, Wernick M, Porter LA, Cukier IH, Collins C, Gray JW, Diebold J, Demetrick DJ, Lee JM: Protein elongation factor EEF1A2 is a putative oncogene in ovarian cancer. Nat Genet 2002, 31: 301-305.

    CAS  PubMed  Google Scholar 

  18. Edmonds BT, Wyckoff J, Yeung YG, Wang Y, Stanley ER, Jones J, Segall J, Condeelis J: Elongation factor-1 alpha is an overexpressed actin binding protein in metastatic rat mammary adenocarcinoma. J Cell Sci 1996, 109: 2705-2714.

    CAS  PubMed  Google Scholar 

  19. Grant AG, Flomen RM, Tizard ML, Grant DA: Differential screening of a human pancreatic adenocarcinoma lambda gt11 expression library has identified increased transcription of elongation factor EF-1 alpha in tumour cells. Int J Cancer 1992, 50: 740-745. 10.1002/ijc.2910500513.

    Article  CAS  PubMed  Google Scholar 

  20. Tomlinson VA, Newbery HJ, Wray NR, Jackson J, Larionov A, Miller WR, Dixon JM, Abbott CM: Translation elongation factor eEF1A2 is a potential oncoprotein that is overexpressed in two-thirds of breast tumours. BMC Cancer 2005, 5: 113. 10.1186/1471-2407-5-113.

    Article  PubMed  Google Scholar 

  21. Gonen H, Dickman D, Schwartz AL, Ciechanover A: Protein synthesis elongation factor EF-1 alpha is an isopeptidase essential for ubiquitin-dependent degradation of certain proteolytic substrates. Adv Exp Med Biol 1996, 389: 209-219.

    Article  CAS  PubMed  Google Scholar 

  22. Gonen H, Smith CE, Siegel NR, Kahana C, Merrick WC, Chakraburtty K, Schwartz AL, Ciechanover A: Protein synthesis elongation factor EF-1 alpha is essential for ubiquitin-dependent degradation of certain N alpha-acetylated proteins and may be substituted for by the bacterial elongation factor EF-Tu. Proc Natl Acad Sci USA 1994, 91: 7648-7652. 10.1073/pnas.91.16.7648.

    Article  CAS  PubMed  Google Scholar 

  23. Kato MV: The mechanisms of death of an erythroleukemic cell line by p53: involvement of the microtubule and mitochondria. Leuk Lymphoma 1999, 33: 181-186.

    CAS  PubMed  Google Scholar 

  24. Kato MV, Sato H, Nagayoshi M, Ikawa Y: Upregulation of the elongation factor-1alpha gene by p53 in association with death of an erythroleukemic cell line. Blood 1997, 90: 1373-1378.

    CAS  PubMed  Google Scholar 

  25. Chen E, Proestou G, Bourbeau D, Wang E: Rapid up-regulation of peptide elongation factor EF-1alpha protein levels is an immediate early event during oxidative stress-induced apoptosis. Exp Cell Res 2000, 259: 140-148. 10.1006/excr.2000.4952.

    Article  CAS  PubMed  Google Scholar 

  26. Duttaroy A, Bourbeau D, Wang XL, Wang E: Apoptosis rate can be accelerated or decelerated by overexpression or reduction of the level of elongation factor-1 alpha. Exp Cell Res 1998, 238: 168-176. 10.1006/excr.1997.3819.

    Article  CAS  PubMed  Google Scholar 

  27. Carneiro NP, Hughes PA, Larkins BA: The eEFIA gene family is differentially expressed in maize endosperm. Plant Mol Biol 1999, 41: 801-813. 10.1023/A:1006391207980.

    Article  CAS  PubMed  Google Scholar 

  28. Xu WL, Wang XL, Wang H, Li XB: Molecular characterization and expression analysis of nine cotton GhEF1A genes encoding translation elongation factor 1A. Gene 2007, 389: 27-35. 10.1016/j.gene.2006.09.014.

    Article  CAS  PubMed  Google Scholar 

  29. Nagashima K, Kasai M, Nagata S, Kaziro Y: Structure of the two genes coding for polypeptide chain elongation factor 1 alpha (EF-1 alpha) from Saccharomyces cerevisiae . Gene 1986, 45: 265-273. 10.1016/0378-1119(86)90024-7.

    Article  CAS  PubMed  Google Scholar 

  30. Nagata S, Nagashima K, Tsunetsugu-Yokota Y, Fujimura K, Miyazaki M, Kaziro Y: Polypeptide chain elongation factor 1 alpha (EF-1 alpha) from yeast: nucleotide sequence of one of the two genes for EF-1 alpha from Saccharomyces cerevisiae . EMBO J 1984, 3: 1825-1830.

    CAS  PubMed  Google Scholar 

  31. Linz JE, Katayama C, Sypherd PS: Three genes for the elongation factor EF-1 alpha in Mucor racemosus . Mol Cell Biol 1986, 6: 593-600.

    Article  CAS  PubMed  Google Scholar 

  32. Hovemann B, Richter S, Walldorf U, Cziepluch C: Two genes encode related cytoplasmic elongation factors 1 alpha (EF-1 alpha) in Drosophila melanogaster with continuous and stage specific expression. Nucleic Acids Res 1988, 16: 3175-3194. 10.1093/nar/16.8.3175.

    Article  CAS  PubMed  Google Scholar 

  33. Abdallah B, Hourdry J, Krieg PA, Denis H, Mazabraud A: Germ cell-specific expression of a gene encoding eukaryotic translation elongation factor 1 alpha (eEF-1 alpha) and generation of eEF-1 alpha retropseudogenes in Xenopus laevis . Proc Natl Acad Sci USA 1991, 88: 9277-9281. 10.1073/pnas.88.20.9277.

    Article  CAS  PubMed  Google Scholar 

  34. Dje MK, Mazabraud A, Viel A, le Maire M, Denis H, Crawford E, Brown DD: Three genes under different developmental control encode elongation factor 1-alpha in Xenopus laevis . Nucleic Acids Res 1990, 18: 3489-3493. 10.1093/nar/18.12.3489.

    Article  CAS  PubMed  Google Scholar 

  35. Krieg PA, Varnum SM, Wormington WM, Melton DA: The mRNA encoding elongation factor 1-alpha (EF-1 alpha) is a major transcript at the midblastula transition in Xenopus . Dev Biol 1989, 133: 93-100. 10.1016/0012-1606(89)90300-X.

    Article  CAS  PubMed  Google Scholar 

  36. Lund A, Knudsen SM, Vissing H, Clark B, Tommerup N: Assignment of human elongation factor 1alpha genes: EEF1A maps to chromosome 6q14 and EEF1A2 to 20q13.3. Genomics 1996, 36: 359-361. 10.1006/geno.1996.0475.

    Article  CAS  PubMed  Google Scholar 

  37. Madsen HO, Poulsen K, Dahl O, Clark BF, Hjorth JP: Retropseudogenes constitute the major part of the human elongation factor 1 alpha gene family. Nucleic Acids Res 1990, 18: 1513-1516. 10.1093/nar/18.6.1513.

    Article  CAS  PubMed  Google Scholar 

  38. Uetsuki T, Naito A, Nagata S, Kaziro Y: Isolation and characterization of the human chromosomal gene for polypeptide chain elongation factor-1 alpha. J Biol Chem 1989, 264: 5791-5798.

    CAS  PubMed  Google Scholar 

  39. Brands JH, Maassen JA, van Hemert FJ, Amons R, Moller W: The primary structure of the alpha subunit of human elongation factor 1. Structural aspects of guanine-nucleotide-binding sites. Eur J Biochem 1986, 155: 167-171. 10.1111/j.1432-1033.1986.tb09472.x.

    Article  CAS  PubMed  Google Scholar 

  40. Knudsen SM, Frydenberg J, Clark BF, Leffers H: Tissue-dependent variation in the expression of elongation factor-1 alpha isoforms: isolation and characterisation of a cDNA encoding a novel variant of human elongation-factor 1 alpha. Eur J Biochem 1993, 215: 549-554. 10.1111/j.1432-1033.1993.tb18064.x.

    Article  CAS  PubMed  Google Scholar 

  41. Lee S, Ann DK, Wang E: Cloning of human and mouse brain cDNAs coding for S1, the second member of the mammalian elongation factor-1 alpha gene family: analysis of a possible evolutionary pathway. Biochem Biophys Res Commun 1994, 203: 1371-1377. 10.1006/bbrc.1994.2336.

    Article  CAS  PubMed  Google Scholar 

  42. Lu XA, Werner D: The complete cDNA sequence of mouse elongation factor 1 alpha (EF 1 alpha) mRNA. Nucleic Acids Res 1989, 17: 442. 10.1093/nar/17.1.442.

    Article  CAS  PubMed  Google Scholar 

  43. Ann DK, Lin HH, Lee S, Tu ZJ, Wang E: Characterization of the statin-like S1 and rat elongation factor 1 alpha as two distinctly expressed messages in rat. J Biol Chem 1992, 267: 699-702.

    CAS  PubMed  Google Scholar 

  44. Ann DK, Moutsatsos IK, Nakamura T, Lin HH, Mao PL, Lee MJ, Chin S, Liem RK, Wang E: Isolation and characterization of the rat chromosomal gene for a polypeptide (pS1) antigenically related to statin. J Biol Chem 1991, 266: 10429-10437.

    CAS  PubMed  Google Scholar 

  45. Lee S, Wolfraim LA, Wang E: Differential expression of S1 and elongation factor-1 alpha during rat development. J Biol Chem 1993, 268: 24453-24459.

    CAS  PubMed  Google Scholar 

  46. Shirasawa T, Sakamoto K, Akashi T, Takahashi H, Kawashima A: Nucleotide sequence of rat elongation factor-1 alpha cDNA. Nucleic Acids Res 1992, 20: 909. 10.1093/nar/20.4.909.

    Article  CAS  PubMed  Google Scholar 

  47. Cavallius J, Merrick WC: Nucleotide sequence of rabbit elongation factor 1 alpha cDNA. Nucleic Acids Res 1992, 20: 1422. 10.1093/nar/20.6.1422.

    Article  CAS  PubMed  Google Scholar 

  48. Kahns S, Lund A, Kristensen P, Knudsen CR, Clark BF, Cavallius J, Merrick WC: The elongation factor 1 A-2 isoform from rabbit: cloning of the cDNA and characterization of the protein. Nucleic Acids Res 1998, 26: 1884-1890. 10.1093/nar/26.8.1884.

    Article  CAS  PubMed  Google Scholar 

  49. Lee S, Francoeur AM, Liu S, Wang E: Tissue-specific expression in mammalian brain, heart, and muscle of S1, a member of the elongation factor-1 alpha gene family. J Biol Chem 1992, 267: 24064-24068.

    CAS  PubMed  Google Scholar 

  50. Lee S, LeBlanc A, Duttaroy A, Wang E: Terminal differentiation-dependent alteration in the expression of translation elongation factor-1 alpha and its sister gene, S1, in neurons. Exp Cell Res 1995, 219: 589-597. 10.1006/excr.1995.1268.

    Article  CAS  PubMed  Google Scholar 

  51. Gao D, Li Z, Murphy T, Sauerbier W: Structure and transcription of the gene for translation elongation factor 1 subunit alpha of zebrafish ( Danio rerio ). Biochim Biophys Acta 1997, 1350: 1-5.

    Article  CAS  PubMed  Google Scholar 

  52. Nowell MA, Power DM, Guerreiro PM, Llewellyn L, Ramsurn VV, Wigham T, Sweeney GE: Cloning and Expression of an Elongation Factor-1alpha in Sea Bream ( Sparus aurata ) Larvae and Adult Tissue. Mar Biotechnol 2000, 2: 173-179.

    CAS  PubMed  Google Scholar 

  53. Mochida K, Matsubara T: Molecular cloning of an elongaton factor 1α and its mRNA localization in testis of the Nile tilapia Oreochromis niloticus . Fish Sci 2002, 68: 830-837. 10.1046/j.1444-2906.2002.00499.x.

    Article  CAS  Google Scholar 

  54. Kinoshita M, Nakata T, Adachi K, Yabe T, Yokoyama Y, Toyohara H, Hirata T, Takayama E, Mikawa S, Kioka N, Sakaguchi M: cDNA cloning of polypeptide chain elongation factor 1α from medaka Oryzias latipes . Fish Sci 1999, 65: 133-137.

    CAS  Google Scholar 

  55. Marchand O, Duffraisse M, Triqueneaux G, Safi R, Laudet V: Molecular cloning and developmental expression patterns of thyroid hormone receptors and T3 target genes in the turbot ( Scophtalmus maximus ) during post-embryonic development. Gen Comp Endocrinol 2004, 135: 345-357. 10.1016/j.ygcen.2003.10.012.

    Article  CAS  PubMed  Google Scholar 

  56. Kitano T, Takamune K, Kobayashi T, Nagahama Y, Abe SI: Suppression of P450 aromatase gene expression in sex-reversed males produced by rearing genetically female larvae at a high water temperature during a period of sex differentiation in the Japanese flounder ( Paralichthys olivaceus ). J Mol Endocrinol 1999, 23: 167-176. 10.1677/jme.0.0230167.

    Article  CAS  PubMed  Google Scholar 

  57. Luckenbach JA, Early LW, Rowe AH, Borski RJ, Daniels HV, Godwin J: Aromatase cytochrome P450: cloning, intron variation, and ontogeny of gene expression in southern flounder ( Paralichthys lethostigma ). J Exp Zoolog 2005, 303: 643-656. 10.1002/jez.a.198.

    Article  Google Scholar 

  58. Fernández-Díaz C, Yúfera M, Cañavate JP, Moyano FJ, Alarcón FJ, Díaz M: Growth and physiological changes during metamorphosis of Senegal sole reared in the laboratory. J Fish Biol 2001, 58: 1-13. 10.1111/j.1095-8649.2001.tb00557.x.

    Article  Google Scholar 

  59. Miwa S, Tagawa M, Inui Y, Hirano T: Thyroxine surge in metamorphosing flounder larvae. Gen Comp Endocrinol 1988, 70: 158-163. 10.1016/0016-6480(88)90105-0.

    Article  CAS  PubMed  Google Scholar 

  60. Power DM, Llewellyn L, Faustino M, Nowell MA, Bjornsson BT, Einarsdottir IE, Canario AV, Sweeney GE: Thyroid hormones in growth and development of fish. Comp Biochem Physiol C Toxicol Pharmacol 2001, 130: 447-459. 10.1016/S1532-0456(01)00271-X.

    Article  CAS  PubMed  Google Scholar 

  61. Kozak M: Point mutations define a sequence flanking the AUG initiator codon that modulates translation by eukaryotic ribosomes. Cell 1986, 44: 283-292. 10.1016/0092-8674(86)90762-2.

    Article  CAS  PubMed  Google Scholar 

  62. Abascal F, Zardoya R, Posada D: ProtTest: selection of best-fit models of protein evolution. Bioinformatics 2005, 21: 2104-2105. 10.1093/bioinformatics/bti263.

    Article  CAS  PubMed  Google Scholar 

  63. Manchado M, Infante C, Asensio E, Cañavate JP: Differential gene expression and dependence on thyroid hormones of two glyceraldehyde-3-phosphate dehydrogenases in the flatfish Senegalese sole ( Solea senegalensis Kaup). Gene 2007, 400: 1-8. 10.1016/j.gene.2007.05.008.

    Article  CAS  PubMed  Google Scholar 

  64. Steinke D, Hoegg S, Brinkmann H, Meyer A: Three rounds (1R/2R/3R) of genome duplications and the evolution of the glycolytic pathway in vertebrates. BMC Biology 2006, 4: 16. 10.1186/1741-7007-4-16.

    Article  PubMed  Google Scholar 

  65. Volff J: Genome evolution and biodiversity in teleost fish. Heredity 2005, 94: 280-294. 10.1038/sj.hdy.6800635.

    Article  CAS  PubMed  Google Scholar 

  66. Hoegg S, Brinkmann H, Taylor JS, Meyer A: Phylogenetic timing of the fish-specific genome duplication correlates with the diversification of teleost fish. J Mol Evol 2004, 59: 190-203. 10.1007/s00239-004-2613-z.

    Article  CAS  PubMed  Google Scholar 

  67. Hoegg S, Meyer A: Hox clusters as models for vertebrate genome evolution. Trends Genet 2005, 21: 421-424. 10.1016/j.tig.2005.06.004.

    Article  CAS  PubMed  Google Scholar 

  68. Moncaut N, Somoza G, Power DM, Canario AV: Five gonadotrophin-releasing hormone receptors in a teleost fish: isolation, tissue distribution and phylogenetic relationships. J Mol Endocrinol 2005, 34: 767-779. 10.1677/jme.1.01757.

    Article  CAS  PubMed  Google Scholar 

  69. Sherwood NM, Adams BA: Gonadotroping-Releasing hormone in fish: Evolution, expression and regulation of the GnRH gene. In Hormones and Their Receptors in Fish Reproduction. Volume 4. Edited by: Melamed P, Sherwood NM. Singapore: World Scientific Publishing Company; 2005:1-39.

    Chapter  Google Scholar 

  70. Krushna Padhi B, Akimenko MA, Ekker M: Independent expansion of the keratin gene family in teleostean fish and mammals: an insight from phylogenetic analysis and radiation hybrid mapping of keratin genes in zebrafish. Gene 2006, 368: 37-45. 10.1016/j.gene.2005.09.016.

    Article  CAS  PubMed  Google Scholar 

  71. 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.

    CAS  PubMed  Google Scholar 

  72. Sidow A: Gen(om)e duplications in the evolution of early vertebrates. Curr Opin Genet Dev 1996, 6: 715-722. 10.1016/S0959-437X(96)80026-8.

    Article  CAS  PubMed  Google Scholar 

  73. Khalyfa A, Bourbeau D, Chen E, Petroulakis E, Pan J, Xu S, Wang E: Characterization of elongation factor-1A (eEF1A-1) and eEF1A-2/S1 protein expression in normal and wasted mice. J Biol Chem 2001, 276: 22915-22922. 10.1074/jbc.M101011200.

    Article  CAS  PubMed  Google Scholar 

  74. Robinson-Rechavi M, Laudet V: Evolutionary rates of duplicate genes in fish and mammals. Mol Biol Evol 2001, 18: 681-683.

    Article  CAS  PubMed  Google Scholar 

  75. Lamberti A, Caraglia M, Longo O, Marra M, Abbruzzese A, Arcari P: The translation elongation factor 1A in tumorigenesis, signal transduction and apoptosis: review article. Amino Acids 2004, 26: 443-448. 10.1007/s00726-004-0088-2.

    Article  CAS  PubMed  Google Scholar 

  76. Ruest LB, Marcotte R, Wang E: Peptide elongation factor eEF1A-2/S1 expression in cultured differentiated myotubes and its protective effect against caspase-3-mediated apoptosis. J Biol Chem 2002, 277: 5418-5425. 10.1074/jbc.M110685200.

    Article  CAS  PubMed  Google Scholar 

  77. Tomlinson VA, Newbery HJ, Bergmann JH, Boyd J, Scott D, Wray NR, Sellar GC, Gabra H, Graham A, Williams AR, Abbott CM: Expression of eEF1A2 is associated with clear cell histology in ovarian carcinomas: overexpression of the gene is not dependent on modifications at the EEF1A2 locus. Br J Cancer 2007, 96: 1613-1620. 10.1038/sj.bjc.6603748.

    Article  CAS  PubMed  Google Scholar 

  78. Shephard KL: Functions for fish mucus. Rev Fish Biol 1994, 4: 401-429. 10.1007/BF00042888.

    Article  Google Scholar 

  79. Evans BI, Fernald RD: Retinal transformation at metamorphosis in the winter flounder ( Pseudopleuronectes americanus ). Vis Neurosci 1993, 10: 1055-1064.

    Article  CAS  PubMed  Google Scholar 

  80. Miwa S, Inui Y: Thyroid hormone stimulates the shift of erythrocyte populations during metamorphosis of the flounder. J Exp Biol 1991, 259: 222-228.

    CAS  Google Scholar 

  81. Miwa S, Yamano K, Inui Y: Thyroid hormone stimulates gastric development in flounder larvae during metamorphosis. J Exp Biol 1991, 261: 424-430.

    Google Scholar 

  82. Schreiber AM, Specker JL: Metamorphosis in the summer flounder ( Paralichthys dentatus ): stage-specific developmental response to altered thyroid status. Gen Comp Endocrinol 1998, 111: 156-166. 10.1006/gcen.1998.7095.

    Article  CAS  PubMed  Google Scholar 

  83. Tagawa M, Aritaki M: Production of symmetrical flatfish by controlling the timing of thyroid hormone treatment in spotted halibut Verasper variegatus . Gen Comp Endocrinol 2005, 141: 184-189. 10.1016/j.ygcen.2005.01.013.

    Article  CAS  PubMed  Google Scholar 

  84. Yamano K, Takano-Ohmuro H, Obinata T, Inui Y: Effect of thyroid hormone on developmental transition of myosin light chains during flounder metamorphosis. Gen Comp Endocrinol 1994, 93: 321-326. 10.1006/gcen.1994.1036.

    Article  CAS  PubMed  Google Scholar 

  85. Yamano K, Miwa S, Obinata T, Inui Y: Thyroid hormone regulates developmental changes in muscle during flounder metamorphosis. Gen Comp Endocrinol 1991, 81: 464-472. 10.1016/0016-6480(91)90174-5.

    Article  CAS  PubMed  Google Scholar 

  86. Abdallah B, Sachs L, Hourdry J, Wegnez M, Denis H, Demeneix B, Mazabraud A: Thyroid hormone regulation of germ cell-specific EF-1 alpha expression during metamorphosis of Xenopus laevis . Int J Dev Biol 1996, 40: 507-514.

    CAS  PubMed  Google Scholar 

  87. Larsen DA, Dickey JT, Dickhoff WW: Quantification of salmon alpha- and thyrotropin (TSH) beta-subunit messenger RNA by an RNase protection assay: regulation by thyroid hormones. Gen Comp Endocrinol 1997, 107: 98-108. 10.1006/gcen.1997.6900.

    Article  CAS  PubMed  Google Scholar 

  88. Pradet-Balade B, Schmitz M, Salmon C, Dufour S, Querat B: Down-regulation of TSH subunit mRNA levels by thyroid hormones in the European eel. Gen Comp Endocrinol 1997, 108: 191-198. 10.1006/gcen.1997.6960.

    Article  CAS  PubMed  Google Scholar 

  89. Riis B, Rattan SI, Clark BF, Merrick WC: Eukaryotic protein elongation factors. Trends Biochem Sci 1990, 15: 420-424. 10.1016/0968-0004(90)90279-K.

    Article  PubMed  Google Scholar 

  90. Ryazanov AG, Rudkin BB, Spirin AS: Regulation of protein synthesis at the elongation stage. New insights into the control of gene expression in eukaryotes. FEBS Lett 1991, 285: 170-175. 10.1016/0014-5793(91)80798-8.

    Article  CAS  PubMed  Google Scholar 

  91. Felsenstein J: PHYLIP-phylogeny inference package (version 3.66). Cladistics 1989, 5: 164-166.

    Google Scholar 

  92. PROSITE Database of protein domains, families and functional sites[http://www.expasy.org/prosite/]

  93. Livak KJ, Schmittgen TD: Analysis of relative gene expression data using real-time quantitative PCR and the 2-ΔΔCtMethod. Methods 2001, 25: 402-408. 10.1006/meth.2001.1262

    Article  CAS  PubMed  Google Scholar 

  94. Manchado M, Infante C, Asensio E, Canavate J, Douglas S: Comparative sequence analysis of the complete set of 40S ribosomal proteins in the Senegalese sole ( Solea senegalensis Kaup) and Atlantic halibut ( Hippoglossus hippoglossus L.) (Teleostei: Pleuronectiformes): phylogeny and tissue- and development-specific expression. BMC Evol Biol 2007, 7: 107. 10.1186/1471-2148-7-107.

    Article  PubMed  Google Scholar 

Download references

Acknowledgements

This study was supported by Genoma España through the Pleurogene Project. CI is grateful to Gonzalo Martínez (ICMAN, CSIC) for his help in the elaboration of this manuscript.

Author information

Authors and Affiliations

Authors

Corresponding authors

Correspondence to Carlos Infante, Esther Asensio, José Pedro Cañavate or Manuel Manchado.

Additional information

Authors' contributions

EA performed the Senegalese sole cultures and hormone exposures, as well as larval and tissue samplings. CI and MM performed the relative quantitation assays and wrote the manuscript. JP-C and MM designed the experiments, and participated in manuscript preparation. All authors read and approved the final manuscript.

Authors’ original submitted files for images

Rights and permissions

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.

Reprints and permissions

About this article

Cite this article

Infante, C., Asensio, E., Cañavate, J.P. et al. Molecular characterization and expression analysis of five different elongation factor 1 alpha genes in the flatfish Senegalese sole (Solea senegalensis Kaup): Differential gene expression and thyroid hormones dependence during metamorphosis. BMC Molecular Biol 9, 19 (2008). https://doi.org/10.1186/1471-2199-9-19

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/1471-2199-9-19

Keywords