<title>Introduction</title> <p>Selenoproteins are a diverse family of proteins containing Selenium (Se) in the form of the non-canonical amino acid selenocysteine (Sec). Selenocysteine, the 21st amino acid, is similar to cysteine (Cys) but with Se replacing Sulphur. In many cases the homologous gene of a known selenoprotein is present with cysteine in the place of Sec in a different genome. Selenocysteine is coded by the opal STOP codon (TGA). Since this codon normally signifies an end to translation, a number of factors combine to achieve the co-translational recoding of TGA to Sec (##FIG##0##Figure 1##). The 3′ UTRs of selenoprotein transcripts contain a stem-loop structure called a <bold>SE</bold>leno<bold>C</bold>ysteine <bold>I</bold>nsertion <bold>S</bold>equence (SECIS) element. This is recognised by the <bold>S</bold>ECIS <bold>B</bold>inding <bold>P</bold>rotein 2 (SBP2), which binds to both the SECIS element and the ribosome. SBP2, in turn, recruits the Sec-specific Elongation Factor EFsec, and the selenocysteine transfer RNA, tRNA<sup>Sec</sup>. SBP2 and tRNA<sup>Sec</sup> form a complex with the tRNA <bold>Se</bold>leno<bold>c</bold>ysteine associated <bold>p</bold>rotein, secp43, which is believed to be involved in the regulation of selenoprotein translation ##REF##16508009##[1]##. Ribosomal protein L30 has recently been shown to interact with the SECIS element and compete with SBP2 for SECIS binding in a Magnesium dependent manner ##REF##15821744##[2]##.</p> <p>Unlike most amino acids, which are aminoacylated onto their cognate tRNAs, Sec is synthesized from serine in a multi-step reaction while bound to its unique tRNA<sup>[Ser]Sec</sup>
##UREF##0##[3]##. Although this reaction is well understood in prokaryotes (e.g. ##UREF##1##[4]##), the details of the eukaryotic pathway remain elusive. It has recently been demonstrated that the protein previously known as <bold>S</bold>oluble <bold>L</bold>iver <bold>A</bold>ntigen/<bold>L</bold>iver <bold>P</bold>ancreas antigen (SLA/LP) is the eukaryotic homolog of bacterial Sec synthetase (SecS), and converts the seryl-tRNA<sup>[Ser]Sec</sup> to selenocysteil-tRNA<sup>[Ser]Sec</sup>
##REF##17194211##[5]##. A <bold>P</bold>hospho<bold>s</bold>eryl <bold>t</bold>RNA <bold>K</bold>inase (PSTK) has also been identified and shown to convert seryl-tRNA to phosphoseryl-tRNA, a likely intermediate to selenocysteil-tRNA ##REF##15317934##[6]##. Finally, Selenophosphate Synthetase 1 and 2 (SPS1 and SPS2), which exhibit sequence similarity, catalyse the formation of mono-selenophosphate, the donor compound of Selenium necessary for the synthesis of selenocysteine. A summary of the selenocysteine biosynthesis and selenoprotein transcription pathways can be seen in ##FIG##0##Figure 1##. Interestingly, SPS2 is itself a selenoprotein. Since selenocysteine, and therefore mono-selenophosphate, is necessary for the expression of SPS2, it has been suggested that SPS1 manufactures basal levels of this compound and the more reactive SPS2 takes over under stimulatory conditions ##REF##16737768##[7]##.</p> <p>Selenoproteins exist in all domains of life, Eukarya, Eubacteria and Archaea. However, no selenoproteins have been found in higher plants (one has been identified in the green alga <italic>Chlamydomonas reinhardtii</italic>) ##REF##12110581##[8]## or fungi. Vertebrate genomes encode up to 25 selenoprotein genes ##REF##12775843##[9]##, ##REF##14710190##[10]##, while invertebrate genomes encode fewer ##REF##15843685##[11]##, ##REF##11493597##[12]##. Three selenoprotein genes have been found in <italic>D. melanogaster</italic>
##REF##11493597##[12]##, <italic>SPS2</italic>, <italic>SelH</italic> and <italic>SelK</italic>. SPS2 is involved in selenoprotein biosynthesis (see above), while SelH and SelK are poorly characterized functionally, but they seem to play an antioxidant role ##REF##17503775##[13]##, ##REF##16962588##[14]##. It has been reported that inhibiting either SelK or SelH expression significantly reduces viability in embryos ##REF##14750951##[15]##. Both SelK and SelH have Cys-paralogs in the <italic>D. melanogaster</italic> genome.</p> <p>Remarkably, only one selenoprotein (Thioredoxin reductase) has been identified in the <italic>C. elegans</italic> genome ##REF##15843685##[11]##. That the entire machinery of selenoprotein synthesis has been conserved in <italic>C. elegans</italic> for synthesizing a single protein had been taken, until very recently, as an indication that selenoproteins are essential for animal life. Indeed, mouse tRNA<sup>Sec</sup> knock-outs have been shown to be lethal in-utero ##REF##9159106##[16]##. Similarly, mutant flies for SPS1 do not contain selenoproteins and are lethal at third instar larvae ##REF##10444382##[17]##. In contrast, Hirosawa-Takamori et al ##REF##14978508##[18]## have reported that mutant flies for EFsec also fail to decode TGA as Sec but are viable and fertile.</p> <p>Recently, we have shown ##REF##17994087##[19]## that one fly, <italic>Drosophila willistoni</italic>, lacks selenoprotein genes, being the first animal reported to lack these proteins. More recently Lobanov et al. have reported that other insect genomes also appear to lack selenoproteins ##REF##18156471##[20]##. In this paper, we extend these results by performing an exhaustive analysis of all available arthropod genomic sequences searching for selenoproteins and selenoprotein factors.</p> <p>First, we analyzed the genomes of the 12 Drosophila species recently sequenced ##REF##17994087##[19]##. In addition to the fact that in <italic>D. willistoni</italic> two of the known insect selenoproteins (SelH and SelK) are Cys-homologs, while the third (SPS2) appears to have been lost ##REF##17994087##[19]##, we have found that many of the genes involved in selenoprotein synthesis have been lost in <italic>D. willistoni</italic>, including the tRNA specific for Sec (tRNA<sup>Sec</sup>). This is strongly indicative that <italic>D. willistoni</italic> not only lacks the <italic>D. melanogaster</italic> selenoprotein reference complement but that it has lost the ability to synthesize selenoproteins altogether. However, other genes thought to be involved in selenoprotein synthesis are as conserved in <italic>D. willistoni</italic> as in the other Drosophila genomes, suggesting that these proteins are involved in additional pathways other than selenoprotein synthesis. Overall, our analyses show selenoprotein evolution in Drosophila to be a very dynamic process; other deviations from the reference selenoprotein complement include the loss of SelK as a selenoprotein in <italic>D. persimilis</italic> and the duplication of SelH in <italic>D. grimshawi</italic>.</p> <p>We have also analyzed the sequences of all other available insect genomes (the mosquitoes <italic>Anopheles gambiae</italic> and <italic>Aedes aegypti</italic>, the honey bee <italic>Apis mellifera</italic>, the wasp <italic>Nasonia vitripennis</italic>, the beetle <italic>Tribolium castaneum</italic> and the silkworm <italic>Bombyx mori</italic>), and found that, while mosquitoes share the selenoprotein complement of <italic>D. melanogaster</italic>, selenoproteins have been lost in the wasp, the honey bee, the silkworm and the beetle. Analysis of available sequence data from other arthropoda (including cDNA, EST, protein and genomic data) suggests that the loss of selenoproteins has been confined to the infraclass Endopterygota, affecting species of all orders investigated (Hymenoptera, Lepidoptera, Diptera and Coleoptera). Interestingly, however, it is not possible to identify a single evolutionary event leading to the loss of selenoproteins in all these species. That most known Diptera still conserve selenoproteins and the mosaic pattern of selenoprotein loss in the other species suggest, instead, multiple independent events of selenoprotein loss in insects. This pattern of gene loss is consistent with a relaxation of the selective constraints acting on insects to maintain selenoproteins, which could be related to the differences in antioxidant defense systems between insects and other metazoans.</p>
<title>Methods</title> <title>Accession Numbers</title> <p>The accession numbers for each of the <italic>D. melanogaster</italic> genes used in this study are as follows: SelK : [GenBank:AAF48111.2]; SelH : [GenBank:AAF48293.3]; SPS2 : [GenBank:AAN10746.2]; SPS1 :[GenBank:AAM70998.1]; SBP2 : [GenBank:AAF50448.2]; EFsec: [GenBank:AAF46721.1]; Secp43 : [GenBank:AAL90383.1]; SecS : [GenBank:AAS65099.1]; PSTK : [GenBank:AAF48985.2]; and tRNA<sup>Sec</sup> : [FLYBASE:FBgn0011987]. The sequences of the other eukaryotic selenoproteins used can be found at <ext-link ext-link-type="uri" xlink:href="http://genome.imim.es/datasets/2007selenoinsects/#1">http://genome.imim.es/datasets/2007selenoinsects/#1</ext-link>.</p> <title>Genome Sequence Data</title> <p>The genomes of the Drosophila species were downloaded from the Drosophila Sequencing Consortium wiki (<ext-link ext-link-type="uri" xlink:href="http://rana.lbl.gov/drosophila/caf1.html">http://rana.lbl.gov/drosophila/caf1.html</ext-link>), we used the Comparative Analysis Freeze 1 (CAF1). The species are: <italic>Drosophila ananassae, Drosophila erecta, Drosophila grimshawi, Drosophila melanogaster, Drosophila mojavensis, Drosophila persimilis, Drosophila pseudoobscura, Drosophila sechellia, Drosophila simulans, Drosophila virilis, Drosophila willistoni</italic> and <italic>Drosophila yakuba.</italic>
</p> <p>The <italic>A. mellifera</italic> genome ##REF##17073008##[21]## (apiMel2, January 2005) was downloaded from UCSC, <ext-link ext-link-type="ftp" xlink:href="ftp://hgdownload.cse.ucsc.edu/goldenPath/apiMel2/">ftp://hgdownload.cse.ucsc.edu/goldenPath/apiMel2/</ext-link>
</p> <p>The <italic>T. castaneum</italic>
##REF##18362917##[22]## (release 1.1, April 2006) sequences were downloaded from NCBI, <ext-link ext-link-type="ftp" xlink:href="ftp://ftp.ncbi.nih.gov/genomes/Tribolium_castaneum/">ftp://ftp.ncbi.nih.gov/genomes/Tribolium_castaneum/</ext-link>
</p> <p>The <italic>A. gambiae</italic> sequences ##REF##12364791##[23]## (anoGam1, February 2003) were downloaded from UCSC, <ext-link ext-link-type="ftp" xlink:href="ftp://hgdownload.cse.ucsc.edu/goldenPath/anoGam1/">ftp://hgdownload.cse.ucsc.edu/goldenPath/anoGam1/</ext-link>
</p> <p>The <italic>A. aegypti</italic> sequences ##REF##17510324##[24]## (AaegL1, March 2006) were downloaded from VectorBase, <ext-link ext-link-type="ftp" xlink:href="ftp://ftp.vectorbase.org/public_data/organism_data/aaegypti/">ftp://ftp.vectorbase.org/public_data/organism_data/aaegypti/</ext-link>
</p> <p>The <italic>N. vitripennis</italic> sequences (Nas1.0, March 8, 2007) were downloaded from the Human Genome Sequencing Center at the Baylor College of Medicine, <ext-link ext-link-type="uri" xlink:href="http://www.hgsc.bcm.tmc.edu/projects/nasonia/">http://www.hgsc.bcm.tmc.edu/projects/nasonia/</ext-link>
</p> <p>The <italic>B. mori</italic> sequences ##UREF##2##[25]## (release 1, October 2003) were downloaded from SilkDB, <ext-link ext-link-type="uri" xlink:href="http://silkworm.genomics.org.cn/silkworm/">http://silkworm.genomics.org.cn/silkworm/</ext-link>
</p> <p>The <italic>D. pulex</italic> sequences were produced by the US Department of Energy Joint Genome Institute (<ext-link ext-link-type="uri" xlink:href="http://www.jgi.doe.gov/">http://www.jgi.doe.gov/</ext-link>) in collaboration with the Daphnia Genomics Consortium (<ext-link ext-link-type="uri" xlink:href="http://daphnia.cgb.indiana.edu">http://daphnia.cgb.indiana.edu</ext-link>).</p> <p>The sequences of all other species investigated were downloaded using the NCBI ENTREZ data retrieval service, <ext-link ext-link-type="uri" xlink:href="http://www.ncbi.nlm.nih.gov/gquery/gquery.fcgi">http://www.ncbi.nlm.nih.gov/gquery/gquery.fcgi</ext-link>.</p> <title>Selenoprotein search in insect genomes</title> <p>The sequences of known selenoproteins— <italic>D. melanogaster</italic> when available, human when not— (see <ext-link ext-link-type="uri" xlink:href="http://genome.imim.es/datasets/2008selenoinsects/#1">http://genome.imim.es/datasets/2008selenoinsects/#1</ext-link>) were searched using the program TBLASTN ##REF##9254694##[26]## against the genomic sequences of each investigated species. The resulting regions of high similarity (see <ext-link ext-link-type="uri" xlink:href="http://genome.imim.es/datasets/2008selenoinsects/#2">http://genome.imim.es/datasets/2008selenoinsects/#2</ext-link>) were then extracted from the target genome, and specifically aligned to the query selenoprotein amino acid sequence using the genewise ##REF##15123596##[27]## and exonerate ##REF##15713233##[28]## programs with default parameters. The output of these programs was manually analyzed to build the exonic structure and the amino acid sequence of the predicted selenoprotein in the target genome. Selenoprotein genes were also searched in cDNA, EST and protein sequences when available for the investigated species.</p> <p>We also investigated all arthropods with sufficient sequence data available (at least 100 genomic, EST, or peptide sequences) for the presence of the known eukaryotic selenoproteins. These were: Amblyomma americanum, Anoplophora glabripennis, Antheraea pernyi, Bactrocera dorsalis, Bactrocera oleae, Bemisia tabaci, Bombyx mandarina, Ceratitis capitata, Daphnia Pulex, Haematobia irritans, Laupala kohalensis, Acyrthosiphon pisum, Homalodisca coagulata, Ixodes scapularis, Locusta migratoria, Nasonia giraulti, Ostrinia furnacalis, Ostrinia nubilalis, Pediculus humanus, Pyrocoelia rufa, Reticulitermes flavipes, Schizaphis graminum, Thermobia domestica, and Triatoma dimidiata.</p> <title>Prediction of tRNA<sup>Sec</sup>
</title> <p>tRNAScanSe ##REF##9023104##[29]## was used to scan each genome for the presence of a tRNA<sup>Sec</sup> gene first with default parameters and then, if no selenocysteine tRNA was found, using only Cove analysis (-C option) which increases the sensitivity. tRNAScanSe uses three models for tRNA<sup>Sec</sup>: Sec(e), Sec(p) and Sec. Sec(e) matches a selenocysteine model based specifically on eukaryotic tRNAs, Sec(p) matches a selenocysteine model based specifically on prokaryotic tRNAs and SeC means that the anticodon identified is UCA, but the predicted tRNA does not match specific SeC models (Lowe T.M., pers. comm.).</p> <p>It must be stressed that tRNAScanSe models for tRNA<sup>Sec</sup> are not as trustworthy as those for other tRNAs due to the small number of tRNA<sup>Sec</sup> sequences available. tRNAScanSe fails to predict a tRNA<sup>Sec</sup> in at least one species (<italic>Takifugu rubripes)</italic> known to code for selenoproteins (Chapple C.E. unpublished data). Therefore the lack of a tRNA<sup>Sec</sup> prediction, although indicative, is not conclusive evidence for the absence of said gene in a given genome.</p> <title>Multiple Alignments of selenoprotein genes</title> <p>The alignments of the amino acid sequences of selenoproteins, selenoprotein cys-homologs and selenoprotein factors were obtained using a combination of the programs clustalw ##REF##12824352##[30]##, t_coffee ##REF##10964570##[31]## and mafft ##REF##12136088##[32]##. Where necessary, the alignments were manually edited using SEAVIEW ##REF##9021275##[33]##. The alignment images presented here were created by jalview ##REF##14960472##[34]##.</p> <title>Phylogenetic trees</title> <p>Phylogenetic trees were built using the online service phylogeny.fr ##REF##18424797##[35]## (“Advanced” Mode, no multiple alignment, all else default) which implements PhyML ##REF##14530136##[36]## for the construction of phylogenetic trees and treedyn ##REF##17032440##[37]## for producing the images presented here.</p> <title>Syntenic Alignments</title> <p>For the Drosophila species, we built the syntenic regions for each of the three selenoprotein genes and each selenoprotein factor. For this we used the annotations produced by the Drosophila Sequencing Consortium ##REF##17994087##[19]##. We selected the 20 surrounding genes of each selenoprotein gene in <italic>D. melanogaster</italic>. We then checked the position of each of these on the target genome. If a gene was annotated as being on the same sequence (scaffold or chromosome depending on the genome) as the target selenoprotein, it was designated “found” and if not, “missing”. For some <italic>D. melanogaster</italic> genes, the target genome had no annotated homolog. In these cases the <italic>D. melanogaster</italic> gene was searched against the target genome using TBLASTN, and the resulting HSP was extended to the full-length protein by genewise and/or exonerate. The distance between the genes was not taken into account, only the order in which they were found. We also built syntenic regions using the whole genome multi-species alignments produced by Lior Pachter's group at UC Berkeley ##REF##17994087##[19]##. For each of the insect selenoproteins and selenoprotein factors, we extracted the region containing the gene in question and the immediately adjacent genes both upstream and downstream.</p> <p>Although we attempted to do the same for the other insects investigated, synteny between them was not sufficient and we were unable to build the necessary alignments.</p> <title>Search for novel selenoproteins in <italic>D. willistoni</italic>
</title> <p>A modified version of the gene prediction software geneid ##REF##10779490##[38]## capable of predicting selenoprotein genes was run on the <italic>D. willistoni</italic> genome. This method has already been described in ##REF##11493597##[12]##. Briefly, it consists of predicting all possible SECIS elements in the target genome then running geneid with the position of these elements given as external information. Geneid will only predict a TGA-containing gene if a SECIS element is found at a suitable distance downstream. The resulting predictions are usually compared against the protein and EST non-redundant sequence databases, as well as against other genome sequences, in search of supporting evidence in the form significant alignments including the aligned Sec-Sec or Sec-Cys.</p> <p>We also searched for all possible TGA-containing open reading frames (ORFs) in the <italic>D. willistoni</italic> genome. This approach is described in full in Taskov et al ##REF##15843685##[11]##. In summary, all TGA-containing ORFs, defined as genomic sequences between two non-TGA stop codons with at least one in-frame TGA and no other in-frame stop codons, are searched in the genome of interest; the resulting sequences are translated in the appropriate frame and compared against the non-redundant protein and EST databases, as well as against other genomes (of insects, in this case). Query sequences where the in-frame TGA is shown to align to either another TGA in the target sequence or to a cysteine residue, and which show conservation extending past the TGA are kept as candidates and further analyzed for the presence of SECIS elements.</p> <title>SECIS prediction</title> <p>The SECIS elements in this paper were predicted using SECISearch ##REF##12775843##[9]##, which can predict potential SECIS elements as well as assess their thermodynamic stability. Three different patterns of decreasing strictness were used allowing us to find both standard and non-standard SECIS elements (see <ext-link ext-link-type="uri" xlink:href="http://genome.imim.es/datasets/2008selenoinsects/#3">http://genome.imim.es/datasets/2008selenoinsects/#3</ext-link>.)</p>
<title>Results</title> <title>Loss of selenoproteins in <italic>D. willistoni</italic>
</title> <p>The three known <italic>D. melanogaster</italic> selenoproteins (SelK, SelH, and SPS2) are found as selenoproteins in all Drosophila genomes except <italic>Drosophila persimilis</italic> and <italic>D. willistoni</italic>. SelK is not a selenoprotein in <italic>D. persimilis</italic>, while in <italic>D. willistoni</italic> SelK and SelH are Cys-homologs, and SPS2 appears to have been lost.</p> <title>SelH</title> <p>As can be seen in ##FIG##1##Figure 2##, SelH appears to be as conserved in <italic>D. willistoni</italic> as in the other Drosophila genomes. However, a number of residues around the Cys/Sec, conserved across all Drosophila (and other Diptera) are mutated in <italic>D. willistoni</italic>, suggesting adaptive changes to compensate the change from Sec to Cys, thereby maintaining the function of SelH. Such compensatory changes have been reported for thioredoxin reductases ##REF##14569031##[39]##. The <italic>SelH</italic> SECIS element, strongly conserved across Drosophila species (##SUPPL##0##Figure S1##), cannot be found in <italic>D. willistoni</italic>, nor can any alternative SECIS element. Interestingly, SelH has been duplicated in <italic>Drosophila grimshawi</italic>, where we found two distinct <italic>SelH</italic> selenoprotein genes (##FIG##1##Figure 2##).</p> <title>SelK</title> <p>With the exception of the Cys to Sec change, SelK is as conserved in <italic>D. willistoni</italic> as in the other Drosophila species (##FIG##2##Figure 3##). The <italic>D. melanogaster</italic> SelK Cys-paralog (CG1840) is only present in the melanogaster group (<italic>D. simulans, D. sechellia, D. melanogaster, D. yakuba, D. erecta</italic> and <italic>D. ananassae</italic>) and so is missing in <italic>D. willistoni</italic>. Indeed, the phylogenetic tree including the SelK and SelK Cys-paralogs in the 12 Drosophila species clearly shows that despite the absence of a Sec residue, <italic>D. willistoni</italic> SelK clusters with the selenoproteins and not the cysteine paralogs (##SUPPL##1##Figure S2##). Interestingly, the <italic>SelK</italic> SECIS element, strongly conserved across selenoprotein containing Drosophila, can still be recognized, although degenerate, in the genome of <italic>D. willistoni</italic> (##SUPPL##2##Figure S3##).</p> <p>SelK is not a selenoprotein in <italic>D. persimilis</italic> either (##FIG##2##Figure 3##). In a previously unreported selenoprotein disabling event, the insertion of a T nucleotide has caused a frameshift, eliminating the in-frame TGA and the subsequent STOP codon, adding nine codons downstream to the next STOP (##SUPPL##3##Figure S4##). Consistent with the disabling mutation, the <italic>SelK</italic> SECIS is degenerate in <italic>D. persimilis</italic> (##SUPPL##2##Figure S3##)</p> <title>SPS2</title> <p>SPS2 appears to have been lost in <italic>D. willistoni</italic>. Indeed, the <italic>D. melanogaster SPS1</italic> and <italic>SPS2</italic> map to the same location in the <italic>D. willistoni</italic> genome, but analysis of the sequence alignments (##FIG##3##Figure 4##) clearly reveals that the <italic>D. willistoni</italic> gene is the SPS1 homolog, as confirmed by the tree built from the multiple alignment of insect SPS1 and SPS2 proteins (##FIG##4##Figure 5##). We could not find a secondary match for the <italic>D. melanogaster</italic> SPS2 in <italic>D. willistoni</italic>, suggesting that this protein is lost in this species.</p> <p>The above analyses strongly indicate that none of the known <italic>D. melanogaster</italic> selenoproteins is a selenoprotein in <italic>D. willistoni</italic>. From these analyses, however, we cannot conclude that <italic>D. willistoni</italic> lacks selenoprotein genes, since other selenoproteins not in <italic>D. melanogaster</italic> could be present in <italic>D. willistoni</italic>. However, we think this is highly unlikely. First, we have compared all known eukaryotic selenoproteins against the <italic>D. willistoni</italic> genome and have found the Cys-homologs typically found in <italic>D. melanogaster</italic> (15-kDa, Glutathione peroxidase (GPx), <italic>thioredoxin</italic> reductase (TR) and SelR). Second, we ran a modified version of the gene predictor geneid ##REF##10779490##[38]## capable of predicting selenoprotein genes and, in addition, we searched for all possible TGA-containing exons in the genome of <italic>D. willistoni</italic> (see <xref ref-type="sec" rid="s2">Methods</xref>). However, after screening the predictions made by these two methods for conservation across the predicted Sec-encoding TGA and potential SECIS elements, all candidates were discarded. The strongest evidence that <italic>D. willistoni</italic> not only lacks selenoprotein genes, but also the capacity to synthesize selenoproteins, comes however from the analysis of the genes involved in selenoprotein biosynthesis. Indeed, that <italic>SPS2</italic> is lost in <italic>D. willistoni</italic> already indicates that selenoprotein synthesis is strongly compromised. Arguably, <italic>SPS1</italic>, present in <italic>D. willistoni</italic> (see below), could rescue SPS2 function. It has been demonstrated however, that selenoprotein biosynthesis is severely impaired in SPS2 knockdown NIH3T3 cells, and that transfection of SPS1 does not restore selenoprotein biosynthesis, suggesting that SPS1 does not complement SPS2 function ##REF##17346238##[40]##. Our analyses indicate, moreover, that not only SPS2, but also other crucial components of the selenoprotein biosynthesis machinery have also been lost in <italic>D. willistoni</italic>.</p> <p>Below we describe our results for each of the factors known to be involved in selenoprotein biosynthesis. We will not focus on ribosomal protein L30 because, as an ubiquitous component of the ribosome, it was present in all species investigated and SECIS binding does not appear to be its primary function.</p> <title>tRNA<sup>Sec</sup>
</title> <p>We used tRNAScanSe to predict tRNA<sup>Sec</sup> genes in each Drosophila genome (see <xref ref-type="sec" rid="s2">Methods</xref>). No suitable tRNA<sup>Sec</sup> could be found in the genome of <italic>D. willistoni</italic>, but high scoring candidates were found in all other Drosophila species (the score of the only tRNA<sup>Sec</sup> prediction in <italic>D. willistoni</italic> was 23.16, while those of the other drosophila ranged from 50.88 to 56.88, see ##TAB##0##Table 1##). Moreover, the <italic>D. willistoni</italic> prediction was clearly less conserved than those in the other Diptera (##FIG##5##Figure 6##), and it did not map to the syntenic region of this gene in the other Drosophila genomes. Instead, the syntenic region in <italic>D. willistoni</italic> shows a deletion spanning the tRNA<sup>Sec</sup> locus (see <ext-link ext-link-type="uri" xlink:href="http://genome.imim.es/datasets/2008selenoinsects/#4">http://genome.imim.es/datasets/2008selenoinsects/#4</ext-link>). The upstream (CG7754) and the downstream (CG12384) immediately adjacent genes are present. These data strongly suggest that <italic>D. willistoni</italic> has lost tRNA<sup>Sec</sup>.</p> <title>EfSec</title> <p>
<italic>EFsec</italic> was found to be highly conserved in the genomes of all the Drosophila species (##SUPPL##4##Figure S5##), but absent in <italic>D. willistoni</italic>. The best candidate found in <italic>D. willistoni</italic> was in fact the gene <italic>EFtau</italic>. No residual (pseudogenised) sequence was found when investigating the syntenic region. Instead, <italic>D. willistoni</italic> shows a gap spanning the <italic>EFsec</italic> locus. Both the upstream (CG10795) and the downstream (CG9707) immediately adjacent genes are present (see <ext-link ext-link-type="uri" xlink:href="http://genome.imim.es/datasets/2008selenoinsects/#4">http://genome.imim.es/datasets/2008selenoinsects/#4</ext-link>).</p> <title>SecS</title> <p>Although conserved in the other Drosophila (##SUPPL##5##Figure S6##), no <italic>SecS</italic> homolog was found in the <italic>D. willistoni</italic> genome. Analysis of the syntenic region across the Drosophila genomes reveals, however, that while the gene upstream of <italic>SecS</italic> (CG2922) is strongly conserved in the Drosophila genomes (including that of <italic>D. willistoni</italic>), a huge deletion, eliminating <italic>SecS</italic> and the gene downstream (CG2919), is present in the genome of <italic>D. willistoni</italic> (see <ext-link ext-link-type="uri" xlink:href="http://genome.imim.es/datasets/2008selenoinsects/#4">http://genome.imim.es/datasets/2008selenoinsects/#4</ext-link>).</p> <title>pstk</title> <p>
<italic>pstk</italic> was also found to be missing in the <italic>D. willistoni</italic> genome and present in the other Drosophila (##SUPPL##6##Figure S7##). However, detailed analysis of the syntenic region in <italic>D. willistoni</italic> of the <italic>D. melanogaster pstk</italic> locus revealed a sequence that could be considered a pseudogenised pstk (see <ext-link ext-link-type="uri" xlink:href="http://genome.imim.es/datasets/2008selenoinsects/#4">http://genome.imim.es/datasets/2008selenoinsects/#4</ext-link>).</p> <title>SBP2</title> <p>An SBP2 homolog can be found in <italic>D. willistoni</italic> (##FIG##6##Figure 7##), in the expected syntenic region. The C-terminal region, strongly conserved across the drosophila, is also conserved in <italic>D. willistoni</italic>indicating that the lack of selenoproteins has not relaxed the selective constraints acting on this region of the protein. This region, however, contains the SBP2 SECIS binding domain. Within this domain a region of 19 amino acids bounded by two Glutamic Acid (E) or Aspartic Acid (D) residues has been shown to be essential for SECIS recognition. Indeed, the specific distance between these two amino acids seems to be a defining feature of the SECIS binding capacity of SBP2 (##REF##12095701##[41]##, Krol A., pers. comm.). Interestingly, this region in <italic>D. willistoni</italic> has an insertion, which could impair SECIS binding capacity (##FIG##6##Figure 7##). The N-terminal region is less conserved overall, but is particularly degenerated in <italic>D. willistoni</italic>.</p> <title>SPS1</title> <p>Although SPS1 is present and highly conserved in <italic>D. willistoni</italic> (##FIG##3##Figure 4##), the phylogenetic tree derived from this protein's amino acid sequence places <italic>D. willistoni</italic> as sister group to the rest of the Drosophila (##FIG##4##Figure 5##), suggesting a relative acceleration of the rate of evolutionary change in this gene.</p> <title>Secp43</title> <p>
<italic>Secp43</italic> was found to be present and highly conserved in all 12 fly genomes, including <italic>D. willistoni</italic> (##SUPPL##7##Figure S8##).</p> <p>In summary, our analyses reveal three different modes of evolution for the proteins involved in selenoprotein biosynthesis in the willistoni branch. A number of genes seem to have been evolving free of selective constraints, and they cannot be found in the genome of <italic>D. willistoni</italic> (tRNA<sup>Sec</sup> , <italic>SecS, EFsec</italic>) or they can only be found as residual pseudogenes (<italic>pstk</italic>). These proteins are probably directly involved in selenoprotein biosynthesis, and unlikely to be involved in unrelated functions. A second group of genes (<italic>SPS1</italic> and <italic>Secp43</italic>), in contrast, are as conserved (or almost as conserved) in <italic>D. willistoni</italic> as in the other Drosophila. Selenoprotein loss does not seem to have greatly influenced their evolutionary rate, and they are therefore likely to have additional (and perhaps more important) functional roles not directly related to selenoprotein biosynthesis. SBP2 exhibits a third, intermediate behavior, with the N-terminal region of the protein showing little similarity to the sequence conserved across the Drosophila genomes, but the C-terminal region strongly conserved. This suggests that, while SBP2 plays an important role in selenoprotein biosynthesis, it may also be involved in other unrelated functions.</p> <title>Loss of selenoproteins in other insects</title> <p>In order to get a clearer understanding of the evolution of selenoprotein genes in the class Insecta, we have analyzed the sequences of all published insect genomes: the Diptera <italic>A. gambiae</italic>
##REF##12364791##[23]## and <italic>A. aegypti</italic>
##REF##17510324##[24]##, the Hymenoptera <italic>A. mellifera</italic> (honey bee) ##REF##17073008##[21]##, and <italic>N. vitripennis</italic> (wasp), the Coleopteron <italic>T. castaneum</italic> (beetle) ##REF##18362917##[22]##, and the Lepidopteron <italic>B. mori</italic> (silkworm) ##UREF##2##[25]##.</p> <p>Lobanov et al. ##REF##18156471##[20]## have recently reported that the species <italic>T. castaneum</italic> and <italic>B. mori</italic> lack selenoproteins. Our results, summarized in ##TAB##1##Table 2##, indicate that in addition to these species and to <italic>D. willistoni</italic>, the Hymenoptera <italic>N. vitripennis</italic> and <italic>A. mellifera</italic> have also lost both selenoproteins and the ability to synthesize them, while mosquitoes maintain the selenoprotein complement of <italic>D. melanogaster</italic>. Like in <italic>D. willistoni</italic>, key components of the selenoprotein machinery have been lost in insects lacking selenoproteins, but not in insects containing them (##TAB##1##Table 2##). Thus, tRNA<sup>Sec</sup>, <italic>EFsec, pstk</italic>, and <italic>SecS</italic>, which could not be found in <italic>D. willistoni</italic> –and we, therefore, speculated were exclusively involved in selenoprotein biosynthesis–also cannot be found in the genomes of other selenoprotein lacking insects (or they can only be found as pseudogenes: <italic>SecS</italic> in <italic>T. castaneum</italic> and <italic>pstk</italic> in <italic>D. willistoni</italic>). They are, however, present in the selenoprotein coding genomes of the mosquitoes. <italic>SecS</italic>, found in the genome of <italic>A. mellifera</italic> is the only exception to this trend (##SUPPL##5##Figure S6##).</p> <p>As with <italic>D. willistoni</italic>, the matches found in these genomes for <italic>EFsec</italic> actually correspond to <italic>EFtau</italic>. The results of our search for tRNA<sup>Sec</sup> (##TAB##0##Table 1##) are at first glance puzzling: consistent with the observed pattern of presence/absence of selenoproteins, no eukaryotic tRNA<sup>Sec</sup> could be predicted in the genomes of <italic>A. mellifera</italic>, and <italic>N. vitripennis</italic>, and only a poor one in the genome of <italic>D. willistoni</italic>, while a very strong candidate can be identified in the selenoprotein containing genome of <italic>A. aegypti</italic>. However, relatively strong tRNA<sup>Sec</sup> predictions are obtained in the genomes of <italic>B. mori</italic> and <italic>T. castaneum</italic>, which lack selenoproteins, while in contrast only a relatively poor prediction is obtained in the selenoprotein containing genome of <italic>A. gambiae</italic>. Close inspection of the predicted tRNA<sup>Sec</sup> genes, however, shows that the tRNA<sup>Sec</sup> predicted in <italic>A. gambiae</italic> strongly resembles that of the other selenoprotein containing Diptera, while the sequence of the tRNA<sup>Sec</sup> predicted in <italic>B. mori, T. castaneum</italic> and <italic>D. willistoni</italic> are very divergent (##FIG##5##Figure 6##).</p> <p>
<italic>SPS1</italic> and <italic>Secp43</italic>, which were as conserved in <italic>D. willistoni</italic> as in the other Drosophila species–and which we therefore speculated are involved in additional functions not related to selenoprotein biosynthesis–can be found in all investigated species, irrespective of selenoprotein coding capacity. Although we found no good <italic>secp43</italic> candidate in <italic>A. gambiae</italic>, we did find a chimeric match against chromosome 2R and chromosome 3L, suggesting an assembly error. More intriguing is the case of <italic>SPS1</italic> in <italic>A. mellifera</italic>. As with <italic>D. willistoni</italic>, the <italic>D. melanogaster SPS1</italic> and <italic>SPS2</italic> map to the same location in the <italic>A. mellifera</italic> genome and analysis of the sequence alignments strongly suggests that the <italic>A. mellifera</italic> gene is the SPS1 homolog (##FIG##3##Figure 4##). In contrast with all other insect <italic>SPS1s</italic>, the <italic>A. mellifera SPS1</italic> contains a TGA codon at the position of the conserved Arginine (R) residue, which is orthologous to the Sec encoding TGA in <italic>SPS2</italic>. In addition, <italic>A. mellifera</italic> is the only species of those which we claim lack selenoproteins that retains the selenoprotein specific factor <italic>SecS</italic>. However, only an unstable SECIS (Free energy: -1.20), exhibiting moreover a bulge immediately after the conserved core (##SUPPL##8##Figure S9##) can be found in the region immediately 3′ of this gene. In addition, as we have already seen, not only can none of the known fly selenoproteins be found in <italic>A. mellifera</italic> but three of the factors involved in selenoprotein synthesis (<italic>pstk, Efsec</italic> and tRNA<sup>Sec</sup>) are also absent. We believe therefore that the <italic>A. mellifera SPS1</italic> is not a selenoprotein, but that the TGA codon is readthrough by an alternative mechanism to produce the full length SPS1 protein. Stop codon readthrough is not uncommon in <italic>D. melanogaster</italic>
##REF##12874049##[42]##, ##REF##17994088##[43]##. TGA is known to be the most “leaky” STOP codon (e.g. ##REF##11972340##[44]##) and has been shown to direct incorporation of Arginine ##REF##11726686##[45]##. Structural analysis of the putative <italic>SPS1</italic> mRNA using mfold ##REF##10329189##[46]##, ##REF##12824337##[47]## showed that the TGA is in a region of high structural stability and forms part of a stem-loop (data not shown) both of which have been shown, together with a Guanine residue observed at the position +3 downstream of the TGA codon, to enhance readthrough in <italic>D. melanogaster</italic>
##REF##12874049##[42]##.</p> <p>
<italic>SBP2</italic>, which was partially conserved in <italic>D. willistoni</italic>, is more confusing when analyzed across all insect genomes. <italic>SBP2</italic> appears to be absent in the selenoprotein lacking genomes of <italic>B. mori</italic> and <italic>T. castaneum</italic>. The conserved C-terminal region is recognizable, however, albeit quite divergent in the genomes of the selenoprotein lacking Hymenoptera. Within Diptera, <italic>SBP2</italic> is present both in mosquitoes and flies. Although no disabling insertion, as in <italic>D. willistoni</italic>, can be found in the SECIS binding domain of other selenoprotein lacking species, overall this domain is slightly more conserved within selenoprotein containing genomes (##FIG##6##Figure 7##).</p> <p>We have extended our analysis by searching genomic, EST, cDNA and peptide sequences available for other arthropod species (including the Arachnids <italic>Ixodes scapularis</italic> and <italic>Amblyomma americanum</italic>, and the Crustacean <italic>Daphnia pulex</italic> ). Results are summarized in ##FIG##7##Figure 8##. Since none of these genomes (except <italic>D. pulex</italic>) is complete, lack of evidence for selenoproteins cannot be taken to indicate total loss of selenoproteins in a given species. Results are interesting, notwithstanding. We have found no evidence of selenoprotein loss in the genomes of any species outside the Endopterygota. Within this infraclass, species of the Hymenoptera, Lepidoptera and Coleoptera orders whose genomes have been completely sequenced do not code for selenoprotein genes, and no evidence of selenoproteins can be found in the partially sequenced species from these orders. In contrast, within Diptera, we found both selenoprotein coding and non-coding species. In summary, existing data do not support selenoprotein loss outside the Endopterygota, but within this infraclass, the loss appears to be generalized across the orders investigated, with the exception of the Diptera.</p>
<title>Discussion</title> <p>Until very recently, selenoproteins were believed to be essential for animal life. The discovery that <italic>D. willistoni</italic>
##REF##17994087##[19]##, <italic>T. castaneum</italic> and <italic>B. mori</italic>
##REF##18156471##[20]## have lost the ability to encode selenoproteins has shaken this belief. Here we present a comprehensive analysis of the selenoproteomes of all available insect genomes, and identify two other insect species, <italic>N. vitripennis</italic> and <italic>A. mellifera</italic>, which have also lost selenoproteins. Through our analysis we have been able to reconstruct a broader picture of selenoprotein extinction in insects.</p> <p>Strong evidence in support of our conclusions comes from the concomitant loss of the factors required for selenoprotein biosynthesis, that we find systematically associated with the loss of selenoproteins (or with their conversion into Cys-homologs, see ##TAB##1##Table 2##). Contradicting our conclusions, however, selenoprotein genes have been reported, after computational analyses, in the genomes of <italic>B. mori</italic>
##UREF##2##[25]## and <italic>A. mellifera</italic>
##REF##17073008##[21]##, both of which we claim do not code for selenoproteins. The case of <italic>B. mori</italic> is, in our opinion, unquestionable. The genome of this species lacks all factors specifically associated with selenoprotein synthesis: tRNA<sup>Sec</sup>, <italic>EFsec</italic>, <italic>SecS, pstk</italic> and <italic>SBP2</italic>. Moreover, the known fly selenoproteins are either Cys-homologs (<italic>SelK</italic>) or lost (<italic>SelH</italic> and <italic>SPS2</italic> ). Finally, the evidence provided in support of novel selenoprotein genes is rather weak, and lacks experimental verification. We therefore confidently conclude that <italic>B. mori</italic> does not contain selenoprotein genes. The case of <italic>A. mellifera</italic> is more controversial (see <xref ref-type="sec" rid="s3">Results</xref>). However, the absence of many essential components of the selenoprotein biosynthesis machinery and of all known fly selenoproteins in <italic>A. mellifera</italic>, as well as the weak SECIS element predicted in <italic>SPS1</italic>, strongly suggest that SPS1 is not a selenoprotein in <italic>A. mellifera</italic>, and that this species lacks the capacity to synthesize selenoproteins.</p> <p>The selective loss of selenoproteins in some insect species allows us to investigate the specificity of the functional role attributed to the known selenoprotein factors and even pinpoint the protein regions responsible for such functional specificity (see <xref ref-type="sec" rid="s3">Results</xref>). We have found that three such factors: tRNA<sup>Sec</sup>, <italic>EFsec</italic>, and <italic>pstk</italic>, are either lost or very degenerate in the genomes of species lacking selenoproteins. We conclude, therefore, that the major role of these genes is indeed related to selenoprotein biosynthesis. We believe that SecS also belongs to this group, although this protein is present in <italic>A. mellifera</italic> (##SUPPL##5##Figure S6##). <italic>SPS1</italic> and <italic>Secp43</italic>, on the other hand, are present and conserved in all genomes investigated irrespective of whether or not they code for selenoproteins. We believe, therefore, that these genes are likely to play a very important role, unrelated to selenoproteins. This does not necessarily mean that they are not involved in selenoprotein biosynthesis. The relative acceleration of the rate of sequence change observed in <italic>D. willistoni</italic> suggests that this may be the case at least for SPS1 (##FIG##4##Figure 5##). Finally, our results indicate that, while SBP2's major role is probably related to selenoprotein biosynthesis, it may also have secondary roles that slow or prevent its elimination from the genome of selenoprotein lacking species. Conversely, selenoprotein loss can also be used to identify novel selenoprotein factors. We have systematically searched for genes that are consistently absent in genomes that do not code for selenoproteins but present in those that do. However, the incomplete status of many of the genomes analyzed confounds analysis of the results of such a search, and we have not been able to confidently identify novel candidate selenoprotein factors.</p> <p>With very few exceptions (<italic>pstk D. willistoni</italic>, and <italic>SecS</italic> in <italic>T. castaneum</italic>), we have not found pseudogenes for most of the lost selenoproteins and selenoprotein factors. Thus, selenoproteins and selenoprotein factors appear to be either present in a given genome as functional proteins or totally absent (without even a recognizable fossil sequence relic). Certainly, the large phylogenetic distance separating many of the species we have investigated confounds the identification of orthologous genomic regions evolving for extended evolutionary times free of selective constraints. The availability of the genomes of the twelve drosophila species ##REF##17994087##[19]##, however, allows us to pinpoint the syntenic regions in <italic>D. willistoni</italic> corresponding to the regions that, in the other Drosophila genomes, contain the selenoprotein factors. Thus, for tRNA<sup>Sec</sup>, <italic>EFsec, pstk</italic> and <italic>SecS</italic> —all lost in <italic>D. willistoni</italic> —we have been able to identify the syntenic regions in the genome of this species (see <ext-link ext-link-type="uri" xlink:href="http://genome.imim.es/datasets/2008selenoinsects/#4">http://genome.imim.es/datasets/2008selenoinsects/#4</ext-link>). For all factors but <italic>pstk</italic>, a deletion in the genome of <italic>D. willistoni</italic> has eliminated the syntenic region, but the upstream and downstream genes can still be found. In the case of <italic>SecS</italic>, the deletion also includes the gene downstream. The fate of SPS2 is harder to determine, since we have not been able to find the syntenic region in <italic>D. willistoni</italic>. It appears therefore, that during Drosophila evolution the entire deletion of non-functional regions and genes is more common—at least in selenoprotein associated genes—than sequence degradation and consequent pseudogenisation. Such behavior is consistent with the dynamic nature of genome micro-structure and the high rate of genome rearrangements observed in Drosophila ##REF##17994087##[19]##.</p> <p>Our data strongly suggest that selenoprotein loss has been confined to the Endopterygota within a general trend of reduction in selenoprotein number in this group compared with other insects (##FIG##7##Figure 8##). Endopterygota are among the most diverse group of insects, comprising eleven orders. Complete genome sequences are only available for species from the orders Coleoptera, Hymenoptera, Lepidoptera, and Diptera. While certainly insufficient for definitive conclusions, available data suggest that selenoprotein loss may be general within Hymenoptera, Coleoptera, and Lepidoptera (all species from these orders with completely sequenced genomes lack selenoproteins and many of the necessary selenoprotein factors, and no evidence of selenoproteins or selenoprotein factors can be found in the species with only partial sequence data available, see ##FIG##7##Figure 8##). In the Diptera, in contrast, only the genome of <italic>D. willistoni</italic>, out of the 14 so far sequenced, has lost selenoproteins. While the phylogenetic relationship between the different orders of this group remain controversial (e.g. ##UREF##4##[49]##–##UREF##5##[51]##), the selective selenoprotein loss that occurred within the Diptera rules out a single evolutionary event as the origin of the pattern of selenoprotein extinction observed in this group. We believe that this pattern is more consistent with a relaxation in insects— accentuated in the Endopterygota— of the selective constraints to maintain selenoproteins that appear to be acting across metazoans (S. Castellano, pers. comm.). In this scenario, different Endopterygota species would be losing selenoproteins independently. Within this general trend of selenoprotein extinction, punctual events of selenoprotein expansion are still possible, such as the duplication of <italic>SelH</italic> observed in <italic>D. grimshawi.</italic>
</p> <p>Lack of intermediate pseudogene evidence confounds the investigation of the evolutionary events that lead to selenoprotein extinction within the Endopterygota. Nevertheless, two contrasting hypotheses can be formulated. One possibility is that selenoprotein genes are lost (or converted to Cys-homologs) first and this triggers the loss of the selenoprotein factors—by rendering these genes functionally irrelevant and easing the selective constraints acting upon then. The alternative hypothesis is that a disabling mutation in a selenoprotein factor occurs first. This renders the mutated species incapable of selenoprotein biosynthesis, and triggers the conversion of selenoproteins into Cys-homologs or simply their elimination from the genome. Primary loss of selenoproteins appears to be more consistent with the existing Drosophila data. Indeed, <italic>D. persimilis</italic> has lost <italic>SelK</italic> as a selenoprotein and, with the exception of the selenoprotein factor <italic>SPS2</italic>, retains only <italic>SelH</italic> as a selenoprotein. Since <italic>SelH</italic> has a Cys-homolog in <italic>D. persimilis</italic>, the selective constraints acting to maintain the selenoprotein factors, which may already be weak in the Drosophila, may be almost non-existent in this species. The alternative hypothesis cannot be discarded, however, because experimental data exists showing that mutant flies for <italic>EFsec</italic>, which fail to decode TGA as Sec, are viable and fertile ##REF##14978508##[18]##. The initial mutation of a selenoprotein factor would therefore not necessarily affect the survival of a species survival, but it would trigger selenoprotein loss. In fact, independent events of selenoprotein extinction may have occurred through different evolutionary paths.</p> <p>Whatever the path leading to selenoprotein loss, our data indicate that, in contrast to all other known metazoans, selenoproteins and selenoprotein synthesis are dispensable in Drosophila and apparently in the entire Endopterygota infraclass. This agrees with the aforementioned data showing viability of mutant <italic>EFsec</italic> flies ##REF##14978508##[18]##. We also think our hypothesis is consistent with the data by Alsina et al. showing that mutant flies for <italic>SPS1</italic> do not contain selenoproteins and are lethal at the third larval instar ##REF##10444382##[17]##. Indeed, our data indicates that one of the major roles of SPS1 is likely to be unrelated to selenoprotein biosynthesis (see <xref ref-type="sec" rid="s3">Results</xref>). Therefore, lethality in mutant <italic>SPS1</italic> flies is probably being induced through the impairment of this major role, and not through its effect on selenoprotein biosynthesis. More difficult to reconcile with our hypothesis of dispensability of selenoproteins in Endopterygota are the experiments in which inhibition of either SelK or SelH expression significantly reduces viability in <italic>D. melanogaster</italic> embryos ##REF##16962588##[14]##, ##REF##14612395##[52]##. The possibility exists, however, that the phenotypes observed in the affected flies result from off-target effects of the RNAi molecules used to interfere with <italic>SelH</italic> and <italic>SelK</italic>.</p> <p>Why the selective constraints acting on selenoproteins have relaxed in Endopterygota remains unknown. However, it is tempting to speculate that they are related to the differences in antioxidant defense systems between Drosophila (and likely other insects) and other metazoans ##REF##11158675##[53]##–##REF##17098255##[56]##. Among other differences, flies do not utilize glutathione peroxidases and have replaced glutathione reductase with non-selenoprotein thioredoxin reductases to reduce glutathione ##REF##11158675##[53]##. The components of the oxidative stress defense system in insects may have thus become independent of selenoproteins, rendering selenoproteins (which in vertebrates have well-established anti-oxidant and cellular redox functions) dispensable. This hypothesis is supported by a recent comparison of the antioxidant defense systems of <italic>A. mellifera, A. gambiae</italic> and <italic>D. melanogaster</italic>
##REF##17069640##[57]## that showed a differential expansion of antioxidant gene families between the Sec-lacking <italic>A. mellifera and</italic> the two Sec-encoding Diptera.</p> <p>Selenoproteins, on the other hand, may also play a functional role in metazoans by sequestering selenium. Sequestration of selenium, whose excess in diet is highly toxic, would thus be compromised in the selenoprotein lacking Endopterygota. We have investigated other non-selenoprotein selenium binding proteins, but have found no signature specific to selenoprotein-lacking insects. How these animals deal with excess selenium represents, therefore, an avenue of future research. In general, that some animals can live without selenoproteins should contribute to a better understanding of the functional role and evolutionary history of this intriguing family of proteins, the most striking exception to the universality of the genetic code.</p>
<p>Conceived and designed the experiments: CEC RG. Performed the experiments: CEC. Analyzed the data: CEC RG. Wrote the paper: CEC RG.</p> <title>Background</title> <p>Selenoproteins are a diverse family of proteins notable for the presence of the 21st amino acid, selenocysteine. Until very recently, all metazoan genomes investigated encoded selenoproteins, and these proteins had therefore been believed to be essential for animal life. Challenging this assumption, recent comparative analyses of insect genomes have revealed that some insect genomes appear to have lost selenoprotein genes.</p> <title>Methodology/Principal Findings</title> <p>In this paper we investigate in detail the fate of selenoproteins, and that of selenoprotein factors, in all available arthropod genomes. We use a variety of <italic>in silico</italic> comparative genomics approaches to look for known selenoprotein genes and factors involved in selenoprotein biosynthesis. We have found that five insect species have completely lost the ability to encode selenoproteins and that selenoprotein loss in these species, although so far confined to the Endopterygota infraclass, cannot be attributed to a single evolutionary event, but rather to multiple, independent events. Loss of selenoproteins and selenoprotein factors is usually coupled to the deletion of the entire no-longer functional genomic region, rather than to sequence degradation and consequent pseudogenisation. Such dynamics of gene extinction are consistent with the high rate of genome rearrangements observed in Drosophila. We have also found that, while many selenoprotein factors are concomitantly lost with the selenoproteins, others are present and conserved in all investigated genomes, irrespective of whether they code for selenoproteins or not, suggesting that they are involved in additional, non-selenoprotein related functions.</p> <title>Conclusions/Significance</title> <p>Selenoproteins have been independently lost in several insect species, possibly as a consequence of the relaxation in insects of the selective constraints acting across metazoans to maintain selenoproteins. The dispensability of selenoproteins in insects may be related to the fundamental differences in antioxidant defense between these animals and other metazoans.</p>
<title>Supporting Information</title>
<p>We thank Sergi Castellano, from Janelia Farm, and Alain Krol, from the Institut de Biologie Molculaire et Cellulaire for helpful discussions, Montserrat Corominas, from Universitat de Barcelona, for critical reading of the manuscript and John Chapple for help in editing. We also thank the sequencing centers for providing the raw genome sequence data necessary for this investigation.</p>
<fig id="pone-0002968-g001" position="float"><object-id pub-id-type="doi">10.1371/journal.pone.0002968.g001</object-id><label>Figure 1</label><caption><title>Selenocysteine biosynthesis and selenoprotein translation pathways.</title><p>Selenoproteins incorporate the amino acid Selenocysteine (Sec) which is coded by the codon UGA, normally a stop codon. The recoding of UGA as a Sec codon is mediated by a structural element on the 3′ Untranslated Region (UTR) of selenoprotein mRNAs, the SElenoCysteine Insertion Sequence (SECIS). This is recognised by the SECIS Binding Protein 2 (SBP2), which binds to both the SECIS element and the ribosome. SBP2, in turn, recruits the Sec-specific Elongation Factor EFsec, and the selenocysteine transfer RNA, tRNA<sup>Sec</sup>. SBP2 and tRNA<sup>Sec</sup> form a complex with the tRNA Selenocysteine associated protein, secp43. Sec is synthesized from serine in a multi-step reaction: Ser-tRNA<sup>[Sec]</sup> is phosphorylated by A Phosphoseryl tRNA Kinase (PSTK) and converted to Sec-tRNA<sup>[Sec]</sup> by Sec synthetase (SecS). Secp43 is also known to be involved in the conversion from seryl to selenocysteyl but its exact role is unclear. Finally, Selenophosphate Synthetase 2 (SPS2), catalyses the formation of mono-selenophosphate, the donor compound of Selenium necessary for the synthesis of selenocysteine, from either selenite (SeO<sub>3</sub>) or from an unstable selenide compound depicted as (Se<sup>2−</sup>). The exact role of SPS1 is still not clear. This figure was partially adapted from ##REF##16737768##[7]##.</p></caption></fig> <fig id="pone-0002968-g002" position="float"><object-id pub-id-type="doi">10.1371/journal.pone.0002968.g002</object-id><label>Figure 2</label><caption><title>Alignment of insect SelH proteins.</title><p>The black arrow shows the position of the selenocysteine (U) residue (cysteine in <italic>D. willistoni</italic>). Here, as in the other alignments, only insect species encoding SelH are shown.</p></caption></fig> <fig id="pone-0002968-g003" position="float"><object-id pub-id-type="doi">10.1371/journal.pone.0002968.g003</object-id><label>Figure 3</label><caption><title>Alignment of insect SelK and SelK cysteine paralogs.</title><p>SelK has been duplicated, producing a Cys-paralog, in species of the melanogaster group. These paralogs are shown with a “C” after the name of the species. The black arrow shows the position of the selenocysteine (U) residue (cysteine in <italic>D. Willistoni</italic> SelK and in the SelK Cys-paralogs).</p></caption></fig> <fig id="pone-0002968-g004" position="float"><object-id pub-id-type="doi">10.1371/journal.pone.0002968.g004</object-id><label>Figure 4</label><caption><title>Alignment of insect SPS1 and SPS2 proteins.</title><p>The black arrow shows the position of the selenocysteine residue in SPS2 and Arganine or Cysteine in SPS1. In <italic>A. mellifera</italic>, we use “?” to denote the codon TGA. Although we believe that in this case, TGA is being readthrough to incorporate Arginine (R), it actually aligns with the Sec-incorporating codon in SPS2.</p></caption></fig> <fig id="pone-0002968-g005" position="float"><object-id pub-id-type="doi">10.1371/journal.pone.0002968.g005</object-id><label>Figure 5</label><caption><title>Phylogenetic tree for the insect SPS1 and SPS2 proteins.</title><p>This tree was built from the alignment of all insects sequences in ##FIG##3##Figure 4##. Note that the D. willistoni sequence (in magenta) clusters with the other SPS1 sequences. This is also the case of the sequence from <italic>A. mellifera</italic> (in blue), in spite of the fact that the in-frame UGA codon in this sequence aligns with the Sec codon in the insect SPS2 sequences.</p></caption></fig> <fig id="pone-0002968-g006" position="float"><object-id pub-id-type="doi">10.1371/journal.pone.0002968.g006</object-id><label>Figure 6</label><caption><title>Alignment of insect tRNA<sup>Sec</sup> sequences.</title><p>The black arrow points to the position of the TCA anticodon.</p></caption></fig> <fig id="pone-0002968-g007" position="float"><object-id pub-id-type="doi">10.1371/journal.pone.0002968.g007</object-id><label>Figure 7</label><caption><title>Alignment of insect SBP2 proteins.</title><p>Alignment of insect SBP2 proteins. The conserved region containing the insertion in <italic>D. willistoni</italic> is bound by a magenta box.</p></caption></fig> <fig id="pone-0002968-g008" position="float"><object-id pub-id-type="doi">10.1371/journal.pone.0002968.g008</object-id><label>Figure 8</label><caption><title>Selenoprotein extinction in arthropoda.</title><p>Species whose genomes do not code for selenoprotein genes are shown in red. Sec-encoding species are shown in green with the number of selenoproteins found in each genome in parentheses next to its name. Species for which the available data was inconclusive are shown in white. The phylogenetic relationships have been taken from the ncbi's Taxonomy database (<ext-link ext-link-type="uri" xlink:href="http://www.ncbi.nlm.nih.gov/Taxonomy/">http://www.ncbi.nlm.nih.gov/Taxonomy/</ext-link>) and the Tree of Life project (<ext-link ext-link-type="uri" xlink:href="http://www.tolweb.org/tree/">http://www.tolweb.org/tree/</ext-link>). The Drosphilidae tree was taken from the Drosophila Sequencing Consortium wiki (<ext-link ext-link-type="uri" xlink:href="http://rana.lbl.gov/drosophila/caf1.html">http://rana.lbl.gov/drosophila/caf1.html</ext-link>).</p></caption></fig>
<table-wrap id="pone-0002968-t001" position="float"><object-id pub-id-type="doi">10.1371/journal.pone.0002968.t001</object-id><label>Table 1</label><caption><title>tRNA<sup>Sec</sup> predictions and their scores in each species.</title></caption><alternatives><table frame="hsides" rules="groups"><colgroup span="1"><col align="left" span="1"/><col align="center" span="1"/></colgroup><thead><tr><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1">tRNA<sup>Sec</sup>
</td></tr></thead><tbody><tr><td align="left" rowspan="1" colspan="1">
<italic>Drosophila</italic>
</td><td align="left" rowspan="1" colspan="1">50.88–56.88, Sec(e)</td></tr><tr><td align="left" rowspan="1" colspan="1">
<bold><italic>D. willistoni</italic></bold>
</td><td align="left" rowspan="1" colspan="1">
<bold>23.16 (Sec)</bold>
</td></tr><tr><td align="left" rowspan="1" colspan="1">
<italic>A. gambiae</italic>
</td><td align="left" rowspan="1" colspan="1">32.53, Sec(p)</td></tr><tr><td align="left" rowspan="1" colspan="1">
<italic>A. aegypti</italic>
</td><td align="left" rowspan="1" colspan="1">53.89, Sec(e)</td></tr><tr><td align="left" rowspan="1" colspan="1">
<bold><italic>A. mellifera</italic></bold>
</td><td align="left" rowspan="1" colspan="1">NA</td></tr><tr><td align="left" rowspan="1" colspan="1">
<bold><italic>N. vitripennis</italic></bold>
</td><td align="left" rowspan="1" colspan="1">NA</td></tr><tr><td align="left" rowspan="1" colspan="1">
<bold><italic>B. mori</italic></bold>
</td><td align="left" rowspan="1" colspan="1">
<bold>55.07, Sec</bold>
</td></tr><tr><td align="left" rowspan="1" colspan="1">
<bold><italic>T. castaneum</italic></bold>
</td><td align="left" rowspan="1" colspan="1">
<bold>45.28, Sec(e)</bold>
</td></tr></tbody></table></alternatives></table-wrap> <table-wrap id="pone-0002968-t002" position="float"><object-id pub-id-type="doi">10.1371/journal.pone.0002968.t002</object-id><label>Table 2</label><caption><title>A summary of the results for each selenoprotein and selenoprotein factor in all completely sequenced insect genomes.</title></caption><alternatives><table frame="hsides" rules="groups"><colgroup span="1"><col align="left" span="1"/><col align="center" span="1"/><col align="center" span="1"/><col align="center" span="1"/><col align="center" span="1"/><col align="center" span="1"/><col align="center" span="1"/><col align="center" span="1"/><col align="center" span="1"/><col align="center" span="1"/><col align="center" span="1"/></colgroup><thead><tr><td align="left" rowspan="1" colspan="1"/><td align="left" rowspan="1" colspan="1">SelH</td><td align="left" rowspan="1" colspan="1">SelK</td><td align="left" rowspan="1" colspan="1">SPS2</td><td align="left" rowspan="1" colspan="1">SPS1</td><td align="left" rowspan="1" colspan="1">SBP2</td><td align="left" rowspan="1" colspan="1">EfSec</td><td align="left" rowspan="1" colspan="1">tRNA<sup>Sec</sup>
</td><td align="left" rowspan="1" colspan="1">Secp43</td><td align="left" rowspan="1" colspan="1">SecS</td><td align="left" rowspan="1" colspan="1">pstk</td></tr></thead><tbody><tr><td align="left" rowspan="1" colspan="1">
<italic>Drosophila</italic>
</td><td align="left" rowspan="1" colspan="1">+</td><td align="left" rowspan="1" colspan="1">+<xref ref-type="table-fn" rid="nt103">*</xref>
</td><td align="left" rowspan="1" colspan="1">+</td><td align="left" rowspan="1" colspan="1">+</td><td align="left" rowspan="1" colspan="1">+</td><td align="left" rowspan="1" colspan="1">+</td><td align="left" rowspan="1" colspan="1">+</td><td align="left" rowspan="1" colspan="1">+</td><td align="left" rowspan="1" colspan="1">+</td><td align="left" rowspan="1" colspan="1">+</td></tr><tr><td align="left" rowspan="1" colspan="1">
<bold><italic>D. willistoni</italic></bold>
</td><td align="left" rowspan="1" colspan="1">Cys</td><td align="left" rowspan="1" colspan="1">Cys</td><td align="left" rowspan="1" colspan="1">−</td><td align="left" rowspan="1" colspan="1">+</td><td align="left" rowspan="1" colspan="1">+</td><td align="left" rowspan="1" colspan="1">−</td><td align="left" rowspan="1" colspan="1">-(?)</td><td align="left" rowspan="1" colspan="1">+</td><td align="left" rowspan="1" colspan="1">−</td><td align="left" rowspan="1" colspan="1">-(?)</td></tr><tr><td align="left" rowspan="1" colspan="1">
<italic>A. gambiae</italic>
</td><td align="left" rowspan="1" colspan="1">+</td><td align="left" rowspan="1" colspan="1">+</td><td align="left" rowspan="1" colspan="1">+</td><td align="left" rowspan="1" colspan="1">+</td><td align="left" rowspan="1" colspan="1">+</td><td align="left" rowspan="1" colspan="1">+</td><td align="left" rowspan="1" colspan="1">+</td><td align="left" rowspan="1" colspan="1">-(?)</td><td align="left" rowspan="1" colspan="1">+</td><td align="left" rowspan="1" colspan="1">+</td></tr><tr><td align="left" rowspan="1" colspan="1">
<italic>A. aegypti</italic>
</td><td align="left" rowspan="1" colspan="1">+</td><td align="left" rowspan="1" colspan="1">+</td><td align="left" rowspan="1" colspan="1">+</td><td align="left" rowspan="1" colspan="1">+</td><td align="left" rowspan="1" colspan="1">+</td><td align="left" rowspan="1" colspan="1">+</td><td align="left" rowspan="1" colspan="1">+</td><td align="left" rowspan="1" colspan="1">+</td><td align="left" rowspan="1" colspan="1">+</td><td align="left" rowspan="1" colspan="1">+</td></tr><tr><td align="left" rowspan="1" colspan="1">
<bold><italic>A. mellifera</italic></bold>
</td><td align="left" rowspan="1" colspan="1">−</td><td align="left" rowspan="1" colspan="1">−</td><td align="left" rowspan="1" colspan="1">−</td><td align="left" rowspan="1" colspan="1">+</td><td align="left" rowspan="1" colspan="1">+</td><td align="left" rowspan="1" colspan="1">−</td><td align="left" rowspan="1" colspan="1">−</td><td align="left" rowspan="1" colspan="1">+</td><td align="left" rowspan="1" colspan="1">+</td><td align="left" rowspan="1" colspan="1">−</td></tr><tr><td align="left" rowspan="1" colspan="1">
<bold><italic>N. vitripennis</italic></bold>
</td><td align="left" rowspan="1" colspan="1">−</td><td align="left" rowspan="1" colspan="1">−</td><td align="left" rowspan="1" colspan="1">−</td><td align="left" rowspan="1" colspan="1">+</td><td align="left" rowspan="1" colspan="1">+</td><td align="left" rowspan="1" colspan="1">−</td><td align="left" rowspan="1" colspan="1">−</td><td align="left" rowspan="1" colspan="1">+</td><td align="left" rowspan="1" colspan="1">−</td><td align="left" rowspan="1" colspan="1">−</td></tr><tr><td align="left" rowspan="1" colspan="1">
<bold><italic>B. mori</italic></bold>
</td><td align="left" rowspan="1" colspan="1">−</td><td align="left" rowspan="1" colspan="1">−</td><td align="left" rowspan="1" colspan="1">−</td><td align="left" rowspan="1" colspan="1">+</td><td align="left" rowspan="1" colspan="1">−</td><td align="left" rowspan="1" colspan="1">−</td><td align="left" rowspan="1" colspan="1">-(?)</td><td align="left" rowspan="1" colspan="1">+</td><td align="left" rowspan="1" colspan="1">−</td><td align="left" rowspan="1" colspan="1">−</td></tr><tr><td align="left" rowspan="1" colspan="1">
<bold><italic>T. castaneum</italic></bold>
</td><td align="left" rowspan="1" colspan="1">−</td><td align="left" rowspan="1" colspan="1">Cys</td><td align="left" rowspan="1" colspan="1">−</td><td align="left" rowspan="1" colspan="1">+</td><td align="left" rowspan="1" colspan="1">−</td><td align="left" rowspan="1" colspan="1">−</td><td align="left" rowspan="1" colspan="1">-(?)</td><td align="left" rowspan="1" colspan="1">+</td><td align="left" rowspan="1" colspan="1">-(?)</td><td align="left" rowspan="1" colspan="1">−</td></tr></tbody></table></alternatives></table-wrap>
<supplementary-material content-type="local-data" id="pone.0002968.s001"><label>Figure S1</label><caption><p>Alignment of insect SelH SECIS elements. Alignment of predicted SelH SECIS elements from each of the insects investigated in which the gene was found. Magenta boxes bound the conserved regions of the SECIS element.</p><p>(0.64 MB DOC)</p></caption></supplementary-material> <supplementary-material content-type="local-data" id="pone.0002968.s002"><label>Figure S2</label><caption><p>Phylogenetic tree built from the alignment of SelK and SelK cysteine paralogs. Phylogenetic tree built from an alignment of SelK and SelK cysteine paralogs (identified with a “C” after the name of the species) across the 12 Drosophila and <italic>A. Gambiae</italic> (used as an outgroup to root the tree). <italic>D. willistoni</italic> is shown in magenta.</p><p>(0.05 MB TIF)</p></caption></supplementary-material> <supplementary-material content-type="local-data" id="pone.0002968.s003"><label>Figure S3</label><caption><p>Alignment of insect SelK SECIS elements. Alignment of predicted SelK SECIS elements from each of the each of the insects investigated in which the gene was found. Magenta boxes bound the conserved regions of the SECIS element. Note the loss of conservation in the <italic>D. persimilis</italic> and <italic>D.willistoni</italic> fossil SECIS elements.</p><p>(0.65 MB DOC)</p></caption></supplementary-material> <supplementary-material content-type="local-data" id="pone.0002968.s004"><label>Figure S4</label><caption><p>SelK cDNA alignment across the 12 Drosophila. Only the terminal region upstream of the stop codon (the last codon in the alignment) is shown. The arrows point to the inserted “T” which has caused a frameshift in <italic>D. persimilis</italic> and the selenocysteine codon (TGA). See ##FIG##2##Figure 3## for the effect of the frameshift on the protein sequence of SelK in <italic>D. persimilis.</italic>
</p><p>(0.55 MB DOC)</p></caption></supplementary-material> <supplementary-material content-type="local-data" id="pone.0002968.s005"><label>Figure S5</label><caption><p>EFsec alignment across all insects investigated.</p><p>(2.41 MB TIF)</p></caption></supplementary-material> <supplementary-material content-type="local-data" id="pone.0002968.s006"><label>Figure S6</label><caption><p>SecS alignment across all insects investigated.</p><p>(2.93 MB TIF)</p></caption></supplementary-material> <supplementary-material content-type="local-data" id="pone.0002968.s007"><label>Figure S7</label><caption><p>PSTK alignment across all insects investigated.</p><p>(1.72 MB TIF)</p></caption></supplementary-material> <supplementary-material content-type="local-data" id="pone.0002968.s008"><label>Figure S8</label><caption><p>secp43 alignment across all insects investigated.</p><p>(3.27 MB DOC)</p></caption></supplementary-material> <supplementary-material content-type="local-data" id="pone.0002968.s009"><label>Figure S9</label><caption><p>
<italic>A. mellifera</italic> SPS1 potential SECIS element. Potential SECIS element 3′ of <italic>A. mellifera</italic> SPS1 gene using the program SECISearch (free energy: -1.20).</p><p>(0.06 MB DOC)</p></caption></supplementary-material>
<table-wrap-foot><fn id="nt101"><p>Species lacking selenoproteins are shown in <bold>bold.</bold> “Sec(e)” matches a model based specifically on eukarotic tRNAs, “Sec(p)” matches a model based specifically on prokaryotic tRNAs and “SeC” means the anticodon identified is TCA, but does not match any specific SeC models (these are much less certain, and could be due to problems in correctly locating the tRNA anticodon)</p></fn></table-wrap-foot> <table-wrap-foot><fn id="nt102"><p>Species lacking selenoproteins are shown in bold. “+” means the gene was present and conserved, “Cys” means the gene was found as a cysteine homolog, “−” means the gene was absent and “-(?)” means that a candidate could be found but the conservation, or score in the case of tRNASec was too low for it to be considered a bona fide homolog.</p></fn><fn id="nt103"><label>*</label><p>Except for <italic>D.persimilis,</italic> see text.</p></fn></table-wrap-foot> <fn-group><fn fn-type="COI-statement"><p><bold>Competing Interests: </bold>The authors have declared that no competing interests exist.</p></fn><fn fn-type="financial-disclosure"><p><bold>Funding: </bold>The work described here is funded by grants from the Spanish Ministery of Education and Science and from the BioSapiens European Network of Excellence to RG. CEC is reciepient of a pre-doctoral fellowship from the Spanish Ministery of Education and Science.</p></fn></fn-group>
<graphic xlink:href="pone.0002968.g001"/> <graphic xlink:href="pone.0002968.g002"/> <graphic xlink:href="pone.0002968.g003"/> <graphic xlink:href="pone.0002968.g004"/> <graphic xlink:href="pone.0002968.g005"/> <graphic xlink:href="pone.0002968.g006"/> <graphic id="pone-0002968-t001-1" xlink:href="pone.0002968.t001"/> <graphic xlink:href="pone.0002968.g007"/> <graphic id="pone-0002968-t002-2" xlink:href="pone.0002968.t002"/> <graphic xlink:href="pone.0002968.g008"/>
<media xlink:href="pone.0002968.s001.doc"><caption><p>Click here for additional data file.</p></caption></media> <media xlink:href="pone.0002968.s002.tif"><caption><p>Click here for additional data file.</p></caption></media> <media xlink:href="pone.0002968.s003.doc"><caption><p>Click here for additional data file.</p></caption></media> <media xlink:href="pone.0002968.s004.doc"><caption><p>Click here for additional data file.</p></caption></media> <media xlink:href="pone.0002968.s005.tif"><caption><p>Click here for additional data file.</p></caption></media> <media xlink:href="pone.0002968.s006.tif"><caption><p>Click here for additional data file.</p></caption></media> <media xlink:href="pone.0002968.s007.tif"><caption><p>Click here for additional data file.</p></caption></media> <media xlink:href="pone.0002968.s008.doc"><caption><p>Click here for additional data file.</p></caption></media> <media xlink:href="pone.0002968.s009.doc"><caption><p>Click here for additional data file.</p></caption></media>
[{"label": ["3"], "element-citation": ["\n"], "surname": ["Small-Howard", "Berry"], "given-names": ["A", "MJ"], "year": ["2006"], "article-title": ["Selenium: Its molecular biology and role in human health, 2nd edition chapter Selenocysteine biosynthesis and incorporation may require supramolecular complexes."], "publisher-loc": ["New York, New York"], "publisher-name": ["Springer Science+Business Media"], "fpage": ["81"], "lpage": ["95"]}, {"label": ["4"], "element-citation": ["\n"], "surname": ["Bock", "Rother", "Leibundgut", "Ban"], "given-names": ["A", "M", "M", "N"], "year": ["2006"], "article-title": ["Selenium: Its molecular biology and role in human health, 2nd edition, chapter Selenium metabolism in prokaryotes,"], "publisher-loc": ["New York, New York"], "publisher-name": ["Springer Science+Business Media"], "fpage": ["9"], "lpage": ["28"]}, {"label": ["25"], "element-citation": ["\n"], "surname": ["Mita", "Kasahara", "Sasaki", "Nagayasu", "Yamada"], "given-names": ["K", "M", "S", "Y", "T"], "year": ["Feb. 2002"], "article-title": ["The Genome Sequence of Silkworm, "], "italic": ["Bombyx mori"], "source": ["DNA Res"], "volume": ["11"], "fpage": ["27"], "lpage": ["35"]}, {"label": ["48"], "element-citation": ["\n"], "surname": ["Chen", "Duan", "Jiang", "Liu", "Zhao"], "given-names": ["P", "J", "L", "Q", "P"], "year": ["2006"], "article-title": ["In silico identification of silkworm selenoproteomes."], "source": ["Chinese Science Bulletin"], "volume": ["51"], "fpage": ["2860"], "lpage": ["2867"]}, {"label": ["49"], "element-citation": ["\n"], "surname": ["Whiting"], "given-names": ["MF"], "year": ["2002"], "article-title": ["Phylogeny of the holometabolous insect orders based on 18s ribosomal DNA: when bad things happen to good data."], "source": ["EXS"], "volume": ["(92)"], "fpage": ["69"], "lpage": ["83"]}, {"label": ["51"], "element-citation": ["\n"], "surname": ["Cameron", "Beckenbach", "Dowton", "Whiting"], "given-names": ["SL", "AT", "MA", "MF"], "year": ["2006"], "article-title": ["Evidence from mitochondrial genomics on interordinal relationships in insects."], "source": ["Arthropod Systematics & Phylogeny"], "volume": ["64"], "fpage": ["27"], "lpage": ["34"]}]
PLoS One. 2008 Aug 13; 3(8):e2968