NiVec database (2011-11-21 release, http://www.ncbi.nlm.nih. gov/VecScreen/UniVec.html). For all contigs longer than 250 bp the open reading frames most likely to encode proteins were identified using the transcripts_to_best_scoring_ORFs.pl script distributed with the 2011-10-29 release of Trinity. The 20 best BLASTP matches for each predicted protein in the NCBI nr database (downloaded 2011-1004) were identified using a local installation of Blast2 [24]. The Blast2 output was used as the input for Blast2GO [25] to assign gene ontology and IEC enzyme codes to proteins, to map enzyme code assignments onto KEGG maps, and to identify the organismal distribution of the best Blast2 hits.RT-PCRRT-PCR was performed using a cDNA pool generated from RNA isolated from a stage 17 T. scripta L cells. Moreover, there was no evidence of any inflammatory cellular embryo. Genes were amplified from the cDNA pool using Taq polymerase (NEB) for 35 cycles with a 60uC annealing temperature and a 1 minute extension time. Primers for each gene (Table 1) were designed to generate a 500?50 bp PCR product and have 65uC annealing temperatures using Primer3 [29].In Situ HybridizationA BMP5 probe was amplified using primers tBmp5NotIR (59TTTGCGGCCGCTGGCTAAGGGAGGACTCT-39) and tBmp5SalF (59TTTGTCGACAGGGGAGAATCACCAAAGA-39). Whole mount stage 15 embryos were hybridized according to [30]. Briefly, embryos were fixed in 4 paraformal-Accession NumbersThe RNA-seq sequences have been deposited in the NCBI Ied kidney origin proteins with previously identified human candidate biomarkers of Sequence Read Archive as accession SRX121294 and theTable 2. Similarity between existing and new T. scripta sequences.length of existing Genbank sequence EF524559.1| Trachemys scripta paired-box protein 1 (Pax1) mRNA, partial cds EF524561.1| Trachemys scripta paired-box protein 3 (Pax3) mRNA, partial cds EF524562.1| Trachemys scripta twist1-like protein mRNA, partial cds EF524563.1| Trachemys scripta dermo-1 (Dermo1) mRNA, partial cds EF524564.1| Trachemys scripta engrailed 1 (En1) mRNA, partial cds EF524565.1| Trachemys scripta gremlin 1 mRNA, partial cds EF524567.1| Trachemys scripta SRY sex determining region Y-box 9 (Sox9) mRNA, partial cds EF527274.1| Trachemys scripta bone morphogenetic protein 4 precursor, mRNA, partial cds EF527276.1| Trachemys scripta homeobox-containing Msx2-like protein (MSX2) mRNA, partial cds AY327846.2|Trachemys scripta bone morphogenetic protein 23148522 2 precursor (BMP-2) mRNA, partial cds. Total length Average identity 614 465 397 614 717 402 340 488 396 1342BLASTN HSP sizes (identical/total length) 576/578 464/465 393/396 447/474, 87/94 717/717 402/402 340/340 488/488 395/396 1273/Length of embryonic transcriptome assembly sequence identity 921 3309 2476 1023 1548 928 3556 1775 735 2789 19060 99.2 99.7 99.8 99.2 94.0 100.0 100.0 100.0 100.0 99.7 99.Existing T. scripta sequences in Genbank were used as queries in a BLASTN search of our assembled sequences. The BLAST HSP sizes represent the sizes of the sequence matches between existing sequences and new T. scripta transcriptome assembly sequences. doi:10.1371/journal.pone.0066357.tTable 3. Top protein hits by species.Species 8,620 5,517 4,651 4,336 2,010 1,087 1,398 627 391 28,637 5,517 67,980 85,348 76.1 79.7 17,368 23.9 20.3 1,095,781 1.2 1.0 261,907 1.4 23.9 675,684 2.2 61.7 34,431 4.9 3.1 17,735 3.8 1.6 2.3 1.6 0.0 0.1 20,676 7.0 1.9 3.7 13,291 15.1 1.2 12.5 17,704 16.2 1.6 10.1 17,368 19.3 1.6 12.2 36,985 30.1 3.4 8.Common nameNumber of top BLAST Number of sequences in NCBI hits vs. transcriptome protein database of sequence.NiVec database (2011-11-21 release, http://www.ncbi.nlm.nih. gov/VecScreen/UniVec.html). For all contigs longer than 250 bp the open reading frames most likely to encode proteins were identified using the transcripts_to_best_scoring_ORFs.pl script distributed with the 2011-10-29 release of Trinity. The 20 best BLASTP matches for each predicted protein in the NCBI nr database (downloaded 2011-1004) were identified using a local installation of Blast2 [24]. The Blast2 output was used as the input for Blast2GO [25] to assign gene ontology and IEC enzyme codes to proteins, to map enzyme code assignments onto KEGG maps, and to identify the organismal distribution of the best Blast2 hits.RT-PCRRT-PCR was performed using a cDNA pool generated from RNA isolated from a stage 17 T. scripta embryo. Genes were amplified from the cDNA pool using Taq polymerase (NEB) for 35 cycles with a 60uC annealing temperature and a 1 minute extension time. Primers for each gene (Table 1) were designed to generate a 500?50 bp PCR product and have 65uC annealing temperatures using Primer3 [29].In Situ HybridizationA BMP5 probe was amplified using primers tBmp5NotIR (59TTTGCGGCCGCTGGCTAAGGGAGGACTCT-39) and tBmp5SalF (59TTTGTCGACAGGGGAGAATCACCAAAGA-39). Whole mount stage 15 embryos were hybridized according to [30]. Briefly, embryos were fixed in 4 paraformal-Accession NumbersThe RNA-seq sequences have been deposited in the NCBI Sequence Read Archive as accession SRX121294 and theTable 2. Similarity between existing and new T. scripta sequences.length of existing Genbank sequence EF524559.1| Trachemys scripta paired-box protein 1 (Pax1) mRNA, partial cds EF524561.1| Trachemys scripta paired-box protein 3 (Pax3) mRNA, partial cds EF524562.1| Trachemys scripta twist1-like protein mRNA, partial cds EF524563.1| Trachemys scripta dermo-1 (Dermo1) mRNA, partial cds EF524564.1| Trachemys scripta engrailed 1 (En1) mRNA, partial cds EF524565.1| Trachemys scripta gremlin 1 mRNA, partial cds EF524567.1| Trachemys scripta SRY sex determining region Y-box 9 (Sox9) mRNA, partial cds EF527274.1| Trachemys scripta bone morphogenetic protein 4 precursor, mRNA, partial cds EF527276.1| Trachemys scripta homeobox-containing Msx2-like protein (MSX2) mRNA, partial cds AY327846.2|Trachemys scripta bone morphogenetic protein 23148522 2 precursor (BMP-2) mRNA, partial cds. Total length Average identity 614 465 397 614 717 402 340 488 396 1342BLASTN HSP sizes (identical/total length) 576/578 464/465 393/396 447/474, 87/94 717/717 402/402 340/340 488/488 395/396 1273/Length of embryonic transcriptome assembly sequence identity 921 3309 2476 1023 1548 928 3556 1775 735 2789 19060 99.2 99.7 99.8 99.2 94.0 100.0 100.0 100.0 100.0 99.7 99.Existing T. scripta sequences in Genbank were used as queries in a BLASTN search of our assembled sequences. The BLAST HSP sizes represent the sizes of the sequence matches between existing sequences and new T. scripta transcriptome assembly sequences. doi:10.1371/journal.pone.0066357.tTable 3. Top protein hits by species.Species 8,620 5,517 4,651 4,336 2,010 1,087 1,398 627 391 28,637 5,517 67,980 85,348 76.1 79.7 17,368 23.9 20.3 1,095,781 1.2 1.0 261,907 1.4 23.9 675,684 2.2 61.7 34,431 4.9 3.1 17,735 3.8 1.6 2.3 1.6 0.0 0.1 20,676 7.0 1.9 3.7 13,291 15.1 1.2 12.5 17,704 16.2 1.6 10.1 17,368 19.3 1.6 12.2 36,985 30.1 3.4 8.Common nameNumber of top BLAST Number of sequences in NCBI hits vs. transcriptome protein database of sequence.