Abstract: The soft-shelled turtle Pelodiscus sinensis is a high-profile turtle species because of its nutritional and medicinal value in Asian countries. However, little is known about the genes that are involved in formation of their nutritional quality traits, especially the molecular mechanisms responsible for unsaturated fatty acid and collagen biosynthesis. In the present study, the transcriptomes from six tissues from Pelodiscus sinensis were sequenced using an Illumina paired-end sequencing platform. We obtained more than 47 million sequencing reads and 73 954 unigenes with an average size of 754 bp by de-novo assembly. In total, 55.19% of the unigenes (40 814) had significant similarity with proteins in the National Center of Biotechnology Information (NCBI) non-redundant protein database and Swiss-Prot database (E-value <10⁻⁵). Of these annotated unigenes, 9 156 and 11 947 unigenes were assigned to 52 gene ontology categories (GO) and 25 clusters of orthologous groups (COG), respectively. In total, 26 496 (35.83%) unigenes were assigned to 242 pathways using the Kyoto Encyclopedia of Genes and Genomes pathway database (KEGG). In addition, we found a number of highly expressed genes involved in the regulation of P. sinensis unsaturated fatty acid biosynthesis and collagen formation, including desaturases, growth factors, transcription factors, and extracellular matrix components. Our data represent the most comprehensive sequence resource available for the Chinese soft-shelled turtle and could provide a basis for new research on this turtle as well as the molecular genetics and functional genomics of other terrapins. To our knowledge, we report for the first time, the large-scale RNA sequencing (RNA-Seq) of terrapin animals and would enrich the knowledge of turtles for future research.

[1]Altschul, S.F., Madden, T.L., Schaffer, A.A., Zhang, J.H., Zhang, Z., Miller, W., Lipman, D.J., 1997. Gapped BLAST and PSI-BLAST: a new generation of protein database search programs. Nucleic Acids Res., 25(17):3389-3402.

[2]Ameri, K., Harris, A.L., 2008. Activating transcription factor 4. Int. J. Biochem. Cell Biol., 40(1):14-21.

[3]Ashburner, M., Ball, C.A., Blake, J.A., Botstein, D., Butler, H., Cherry, J.M., 2000. Gene ontology: tool for the unification of biology. Nat. Genet., 25(1):25-29.

[4]Bhogal, R.K., Stoica, C.M., McGaha, T.L., Bona, C.A., 2005. Molecular aspects of regulation of collagen gene expression in fibrosis. J. Clin. Immunol., 25(6):592-603.

[5]Camacho, C., Coulouris, G., Avagyan, V., Ma, N., Papadopoulos, J., Bealer, K., Madden, T.L., 2009. BLAST⁺: architecture and applications. BMC Bioinf., 10(1):421.

[6]Cameron, M., Williams, H.E., Cannane, A., 2004. Improved gapped alignment in BLAST. IEEE/ACM Trans. Comput. Biol. Bioinf., 1(3):116-129.

[7]Cao, H., Gerhold, K., Mayers, J.R., Wiest, M.M., Watkins, S.M., Hotamisligil, G.S., 2008. Identification of a lipokine, a lipid hormone linking adipose tissue to systemic metabolism. Cell, 134(6):933-944.

[8]Chang, L., Chen, J.J., Xiao, Y.M., Xia, Y.P., 2011. De novo characterization of Lycoris sprengeri transcriptome using Illumina GA II. Afr. J. Biotechnol., 10(57):12147-12155.

[9]Collins, L.J., Biggs, P.J., Voelckel, C., Joly, S., 2008. An approach to transcriptome analysis of non-model organisms using short-read sequences. Genome Inf., 21:3-14.

[10]Conesa, A., Götz, S., García-Gómez, J.M., Terol, J., Talón, M., Robles, M., 2005. Blast2GO: a universal tool for annotation, visualization and analysis in functional genomics research. Bioinformatics, 21(18):3674-3676.

[11]Eble, J.A., Kassner, A., Niland, S., Mörgelin, M., Grifka, J., Grässel, S., 2006. Collagen XVI harbors an integrin α1β1 recognition site in its C-terminal domains. J. Biol. Chem., 281(35):25745-25756.

[12]Feldmeyer, B., Wheat, C.H., Krezdorn, N., Rotter, B., Pfenninger, M., 2011. Short read Illumina data for the de novo assembly of a non-model snail species transcriptome (Radix balthica, Basommatophora, Pulmonata), and a comparison of assembler performance. BMC Genomics, 12(1):317.

[13]Gao, X.G., Han, J.B., Lu, Z.C., Li, Y.F., He, C.B., 2012. Characterization of the spotted seal Phoca largha transcriptome using Illumina paired-end sequencing and development of SSR markers. Comp. Biochem. Physiol. Part D: Genom. Proteom., 7(3):277-284.

[14]Garg, R., Patel, R.K., Tyagi, A.K., Jain, M., 2011. De novo assembly of chickpea transcriptome using short reads for gene discovery and marker identification. DNA Res., 18(1):53-63.

[15]Gong, X., Niu, C.J., Zhang, Z.B., 2011. cDNA cloning and tissue expression for l-gulonolactone oxidase gene in soft-shelled turtle Pelodiscus sinensis a species with the ability to synthesize ascorbic acid. Fish Sci., 77(4):547-555.

[16]Grabherr, M.G., Haas, B.J., Yassour, M., Levin, J.Z., Thompson, D.A., Amit, I., Adiconis, X., Fan, L., Raychowdhury, R., Zeng, Q.D., 2011. Full-length transcriptome assembly from RNA-Seq data without a reference genome. Nat. Biotech., 29(7):644-652.

[17]Harris, M.A., Clark, J., Ireland, A., Lomax, J., Ashburner, M., Foulger, R., Eilbeck, K., Lewis, S., Marshall, B., Mungall, C., et al., 2004. The gene ontology (GO) database and informatics resource. Nucleic Acids Res., 32:D258-D261.

[18]Hegedűs, Z., Zakrzewska, A., Ágoston, V.C., Ordas, A., Ŕacz, P., Mink, M., Spaink, H.P., Meijer, A.H., 2009. Deep sequencing of the zebrafish transcriptome response to mycobacterium infection. Mol. Immunol., 46(15):2918-2930.

[19]Heine, U.I., Munoz, E.F., Flanders, K.C., Roberts, A.B., Sporn, M.B., 1990. Colocalization of TGF-β1 and collagen I and III, fibronectin and glycosaminoglycans during lung branching morphogenesis. Development, 109(1):29-36.

[20]Hillier, L.W., Reinke, V., Green, P., Hirst, M., Marra, M.A., Waterson, R.H., 2009. Massively parallel sequencing of the polyadenylated transcriptome of C. elegans. Genome Res., 19(4):657-666.

[21]Hou, R., Bao, Z.M., Wang, S., Su, H.L., Li, Y., Du, H., Hu, J., Wang, S., Hu, X., 2011. Transcriptome sequencing and de novo analysis for Yesso Scallop (Patinopecten yessoensis) using 454 GS FLX. PLoS One, 6(6):e21560.

[22]Huang, C.H., Lin, W.Y., Chu, J.H., 2005. Dietary lipid level influences fatty acid profiles, tissue composition, and lipid peroxidation of soft-shelled turtle, Pelodiscus sinensis. Comp. Biochem. Physiol. A: Mol. Integr. Physiol., 142(3):383-388.

[23]Kanehisa, M., Goto, S., Kawashima, S., Okuno, Y., Hattori, M., 2004. The KEGG resource for deciphering the genome. Nucleic Acids Res., 32(S1):D277-D280.

[24]Li, R., Zhu, H., Ruan, J., Qian, W., Fang, X., Shi, Z., Li, Y., Li, S., Shan, G., Kristiansen, K., et al., 2010. De novo assembly of human genomes with massively parallel short read sequencing. Genome Res., 20(2):265-272.

[25]Li, X.L., Zhang, C.L., Fang, W.H., Lin, F.C., 2008. White-spot disease of Chinese soft-shelled turtles (Trionyx sinensis) caused by Paecilomyces lilacinus. J. Zhejiang Univ. Sci. B, 9(7):578-581.

[26]Liu, C.C., Liu, Y., Jin, Y.Z., 2010. Extraction and antioxidant activity of collagen from the Chinese soft-shelled turtle (Pelodiscus sinensis). Adv. Mater. Res., 152-153:1788-1792.

[27]Lu, J.F., Wan, Q., Yin, Z.M., Lin, L., Weng, S.B., Ye, Y.W., Jiang, S.T., 2010. Extraction and characterization of collagen from calipash of Chinese soft-shelled turtle (Pelodiscus sinensis). J. Fish China, 34:981-988 (in Chinese).

[28]Metzker, M.L., 2010. Sequencing technologies—the next generation. Nat. Rev. Genet., 11(1):31-46.

[29]Ministry of Agriculture and Fisheries Bureau of China, 2011. China Fishery Statistical Yearbook. China Agriculture Press, Beijing, China, p.41 (in Chinese).

[30]Mortazavi, A., Williams, B.A., McCue, K., Schaeffer, L., Wold, B., 2008. Mapping and quantifying mammalian transcriptomes by RNA-Seq. Nat. Methods, 5(7):621-628.

[31]Myers, J.C., Yang, H.Y., D′Ippolito, J.A., Presente, A., Miller, M.K., Dion, A.S., 1994. The triple-helical region of human type XIX collagen consists of multiple collagenous subdomains and exhibits limited sequence homology to alpha 1(XVI). J. Biol. Chem., 269(28):18549-18557.

[32]Nobuhiro, N., Hatsumi, K., Shizuka, K., Masanobu, M., 2009. Preparation and characterization of collagen from soft-shelled turtle (Pelodiscus Sinensis) skin for biomaterial applications. J. Biomat. Sci., Polymer Ed., 20(5-6):567-576.

[33]Parchman, T.L., Geist, K.S., Grahnen, J.A., Benkman, C.W., Buerkle, C.A., 2010. Transcriptome sequencing in an ecologically important tree species: assembly, annotation, and marker discovery. BMC Genomics, 11(1):180.

[34]Pereira, S.L., Leonard, A.E., Mukerji, P., 2003. Recent advances in the study of fatty acid desaturases from animals and lower eukaryotes. Prostaglandins Leukot. Essent. Fatty Acids, 68(2):97-106.

[35]Schroeder, A., Mueller, O., Stocker, S., Salowsky, R., Leiber, M., Gassmann, M., Lightfoot, S., Menzel, W., Granzow, M., Ragg, T., 2006. The RIN: an RNA integrity number for assigning integrity values to RNA measurements. BMC Mol. Biol., 7(1):3.

[36]Schuster, S.C., 2008. Next-generation sequencing transforms today’s biology. Nat. Methods, 5(1):16-18.

[37]Storch, J., McDermott, L., 2009. Structural and functional analysis of fatty acid-binding proteins. J. Lipid Res., 50(Suppl.):S126-S131.

[38]Tatusov, R.L., Natale, D.A., Garkavtsev, I.V., Tatusova, T.A., Shankavaram, U.T., Rao, B.S., Kiryutin, B., Galperin, M.Y., Fedorova, N.D., Koonin, E.V., 2001. The COG database: new developments in phylogenetic classification of proteins from complete genomes. Nucleic Acids Res., 29(1):22-28.

[39]Wang, D.Z., Tang, Z.R., Tan, Y.J., 1997. Biochemical compositions of Chinese soft-shelled turtle (Trionyx sinensis). I. Contents of normal nutrients and composition of muscle fatty acids. Acta Hydrobiol. Sin., 21:199-305 (in Chinese).

[40]Wang, W., Lian, N., Li, L., Moss, H.E., Wang, W., Perrien, D.S., Elefteriou F., Yang, X., 2009. Atf4 regulates chondrocyte proliferation and differentiation during endochondral ossification by activating Ihh transcription. Development, 136(24):4143-4153.

[41]Wang, X.W., Luan, J.B., Li, J.M., Bao, Y.Y., Zhang, C.X., Liu, S.S., 2010b. De novo characterization of a whitefly transcriptome and analysis of its gene expression during development. BMC Genomics, 11(1):400.

[42]Wang, Z., Gerstein, M., Snyder, M., 2009. RNA-Seq: a revolutionary tool for transcriptomics. Nat. Rev. Genet., 10(1):57-63.

[43]Wang, Z.Y., Fang, B.P., Chen, J.Y., Zhang, X.J., Luo, Z.X., Huang, L.F., Chen, X.L., 2010a. De novo assembly and characterization of root transcriptome using Illumina paired-end sequencing and development of cSSR markers in sweetpotato (Ipomoea batatas). BMC Genomics, 11(1):726.

[44]Watanabe, R., Fujii, H., Odani, S., Sakakibara, J., Yamamoto, A., Ito, M., Ono, T., 1994. Molecular cloning of a cDNA encoding a novel fatty acid-binding protein from rat skin. Biochem. Biophys. Res. Commun., 200(1):253-259.

[45]Wei, W.L., Qi, X.Q., Wang, L.H., Zhang, Y.X., Hua, W., Li, D.H., Lv, H.X., Zhang, X.R., 2011. Characterization of the sesame (Sesamum indicum L.) global transcriptome using Illumina paired-end sequencing and development of EST-SSR markers. BMC Genomics, 12(1):451.

[46]Xiang, L.X., He, D., Dong, W.R., Zhang, Y.W., Shao, J.Z., 2010. Deep sequencing-based transcriptome profiling analysis of bacteria-challenged Lateolabrax japonicus reveals insight into the immune relevant genes in marine fish. BMC Genomics, 11(1):472.

[47]Xie, F.L., Burklew, C.E., Yang, Y.F., Liu, M., Xiao, P., Zhang, B.H., Qiu, D.Y., 2012. De novo sequencing and a comprehensive analysis of purple sweet potato (Impomoea batatas L.) transcriptome. Planta, 236(1):101-113.

[48]Xue, J., Bao, Y.Y., Li, B.L., Cheng, Y.B., Peng, Z.Y., Liu, H., Xu, H.J., Zhu, Z.R., Lou, Y.G., Cheng, J.A., Zhang, C.X., 2010. Transcriptome analysis of the brown planthopper Nilaparvata lugens. PloS One, 5(12):e14233.

[49]Yao, B.J., Zhao, Y., Wang, Q., Zhang, M., Liu, M.C., Liu, H.L., Li, J., 2012a. De novo characterization of the antler tip of Chinese Sika deer transcriptome and analysis of gene expression related to rapid growth. Mol. Cell Biochem., 364(1-2):93-100.

[50]Yao, B.J., Zhao, Y., Zhang, H.S., Zhang, M., Liu, M.C., Liu, H.L., Li, J., 2012b. Sequencing and de novo analysis of the Chinese Sika deer antler-tip transcriptome during the ossification stage using Illumina RNA-Seq technology. Biotechnol. Lett., 34(5):813-822.

[51]Ye, J., Fang, L., Zheng, H., Zhang, Y., Chen, J., Zhang, Z., Wang, J., Li, S., Li, R., Bolund, L., Wang, J., 2006. WEGO: a web tool for plotting GO annotations. Nucleic Acids Res., 34:W293-W297.

[52]Zhan, X.A., Xu, Z.R., Qian, L.C., 2000. Muscle and fat quality of Chinese soft-shelled turtle. J. Zhejiang Univ. (Agric. & Life Sci.), 26(4):457-460 (in Chinese).

[53]List of electronic supplementary materials

[54]Table S1 KEGG pathways for Pelodiscus sinensis