Background Despite their species abundance and primary economic importance, genomic information about copepods is still limited. In particular, genomic resources are lacking for the copepod Calanus sinicus, which is a dominant species in the coastal waters of East Asia. In this study, we performed de novo transcriptome sequencing to produce a large number of expressed sequence tags for the copepod C. sinicus. Results Copepodid larvae and adults were used as the basic material for transcriptome sequencing. Using 454 pyrosequencing, a total of 1,470,799 reads were obtained, which were assembled into 56,809 high quality expressed sequence tags. Based on their sequence similarity to known proteins, about 14,000 different genes were identified, including members of all major conserved signaling pathways. Transcripts that were putatively involved with growth, lipid metabolism, molting, and diapause were also identified among these genes. Differentially expressed genes related to several processes were found in C. sinicus copepodid larvae and adults. We detected 284,154 single nucleotide polymorphisms (SNPs) that provide a resource for gene function studies. Conclusion Our data provide the most comprehensive transcriptome resource available for C. sinicus. This resource allowed us to identify genes associated with primary physiological processes and SNPs in coding regions, which facilitated the quantitative analysis of differential gene expression. These data should provide foundation for future genetic and genomic studies of this and related species.
References
[1]
Humes AG (1994) How many copepods? Hydrobiologia 292: 1–7.
[2]
Bron JE, Frisch D, Goetze E, Johnson SC, Lee CE, et al. (2011) Observing copepods through a genomic lens. Frontiers in Zoology 8: 22.
[3]
Beaugrand G, Brander KM, Lindley JA, Souissi S, Reid PC (2003) Plankton effect on cod recruitment in the North Sea. Nature 426: 661–664.
[4]
Frangoulis C, Christou E, Hecq J (2004) Comparison of marine copepod outfluxes: nature, rate, fate and role in the carbon and nitrogen cycles. Advances in Marine Biology 47: 253–309.
[5]
Richardson AJ (2008) In hot water: zooplankton and climate change. ICES Journal of Marine Science: Journal du Conseil 65: 279–295.
[6]
Minxiao W, Song S, Chaolun L, Xin S (2011) Distinctive mitochondrial genome of Calanoid copepod Calanus sinicus with multiple large non-coding regions and reshuffled gene order: Useful molecular markers for phylogenetic and population studies. BMC Genomics 12: 73.
[7]
Hudson ME (2008) Sequencing breakthroughs for genomic ecology and evolutionary biology. Molecular Ecology Resources 8: 3–17.
[8]
Collins LJ, Biggs PJ, Voelckel C, Joly S (2008) An approach to transcriptome analysis of non-model organisms using short-read sequences. Genome Informatics 21: 3–14.
[9]
Rothberg JM, Leamon JH (2008) The development and impact of 454 sequencing. Nature Biotechnology 26: 1117–1124.
[10]
Marguerat S, B？hler J (2009) RNA-seq: from technology to biology. Cellular and Molecular Life Sciences 67: 569–579.
[11]
Morozova O, Hirst M, Marra MA (2009) Applications of new sequencing technologies for transcriptome analysis. Annual Review of Genomics and Human Genetics 10: 135–151.
[12]
Emrich SJ, Barbazuk WB, Li L, Schnable PS (2007) Gene discovery and annotation using LCM-454 transcriptome sequencing. Genome Research 17: 69–73.
[13]
Vera JC, Wheat CW, Fescemyer HW, Frilander MJ, Crawford DL, et al. (2008) Rapid transcriptome characterization for a nonmodel organism using 454 pyrosequencing. Molecular Ecology 17: 1636–1647.
[14]
Parchman TL, Geist KS, Grahnen JA, Benkman CW, Buerkle CA (2010) Transcriptome sequencing in an ecologically important tree species: assembly, annotation, and marker discovery. BMC Genomics 11: 180.
[15]
Fraser B, Weadick C, Janowitz I, Rodd H, Hughes K (2011) Sequencing and characterization of the guppy (Poecilia reticulata) transcriptome. BMC Genomics 12: 202.
[16]
Ben EC, Nathan S, Kristen P, Yuichiro S, Siegfried R, et al. (2011) The maternal and early embryonic transcriptome of the milkweed bug Oncopeltus fasciatus. BMC Genomics 12: 61.
[17]
Gregory R, Darby AC, Irving H, Coulibaly MB, Hughes M, et al. (2011) A De Novo Expression Profiling of Anopheles funestus, Malaria Vector in Africa, Using 454 Pyrosequencing. PLoS One 6: e17418.
[18]
Hou R, Bao Z, Wang S, Su H, Li Y, et al. (2011) Transcriptome Sequencing and De Novo Analysis for Yesso Scallop (Patinopecten yessoensis) Using 454 GS FLX. PLoS One 6: e21560.
[19]
Barreto FS, Moy GW, Burton RS (2011) Interpopulation patterns of divergence and selection across the transcriptome of the copepod Tigriopus californicus. Molecular Ecology 20: 560–572.
[20]
Meyer E, Aglyamova G, Wang S, Buchanan-Carter J, Abrego D, et al. (2009) Sequencing and de novo analysis of a coral larval transcriptome using 454 GSFlx. BMC Genomics 10: 219.
[21]
Hulsemann K (1994) Calanus sinicus Brodsky and C. jashmovi, nom. nov.(Copepoda: Calanoida) of the north-western Pacific Ocean: a comparison, with notes on the integumental pore pattern in Calanus s. str. Invertebrate Systematics 8: 1461–1482.
[22]
Chen Q (1964) Study on the reproduction, sex ratio and body size of Calanus sinicus (in Chinese with English abstract). Oceanologia et Limnologia Sinica 6: 272–287.
[23]
Uye S (2000) Why does Calanus sinicus prosper in the shelf ecosystem of the Northwest Pacific Ocean? ICES Journal of Marine Science: Journal du Conseil 57: 1850–1855.
[24]
Kang HK, Lee CR, Choi KH (2011) Egg production rate of the copepod Calanus sinicus off the Korean coast of the Yellow Sea during spring. Ocean Science Journal 46: 133–143.
[25]
Wang S, Li C, Sun S, Ning X, Zhang W (2009) Spring and autumn reproduction of Calanus sinicus in the Yellow Sea. Marine Ecology Progress Series 379: 123–133.
[26]
Sun S, Zhang G (2005) Over-summering strategy of Calanus sinicus. GLOBEC Int Newsl 11: 34.
[27]
Pu XM, Sun S, Yang B, Zhang GT, Zhang F (2004) Life history strategies of Calanus sinicus in the southern Yellow Sea in summer. Journal of Plankton Research 26: 1059–1068.
[28]
Bi H, Sun S, Gao S, Zhang G (2001) The ecological characteristics of zooplankton community in the Bo hai Sea II. The distribution of copepoda abundance and seasonal dynamics. Acta Ecologica Einica 21: 177–185.
[29]
Zhang W, Tang D, Yang B, Gao S, Sun J, et al. (2009) Onshore–offshore variations of copepod community in northern South China Sea. Hydrobiologia 636: 257–269.
[30]
O'Neil S, Dzurisin J, Carmichael R, Lobo N, Emrich S, et al. (2010) Population-level transcriptome sequencing of nonmodel organisms Erynnis propertius and Papilio zelicaon. BMC Genomics 11: 310.
[31]
Zeng V, Villanueva K, Ewen-Campen B, Alwes F, Browne W, et al. (2011) De novo assembly and characterization of a maternal and developmental transcriptome for the emerging model crustacean Parhyale hawaiensis. BMC Genomics 12: 581.
[32]
Brockman W, Alvarez P, Young S, Garber M, Giannoukos G, et al. (2008) Quality scores and SNP detection in sequencing-by-synthesis systems. Genome Research 18: 763–770.
[33]
Irigoien X (2004) Some ideas about the role of lipids in the life cycle of Calanus finmarchicus. Journal of Plankton Research 26: 259–263.
[34]
Dubrovsky EB (2005) Hormonal cross talk in insect development. Trends in Endocrinology and Metabolism 16: 6–11.
[35]
Kidokoro K, Iwata K, Fujiwara Y, Takeda M (2006) Effects of juvenile hormone analogs and 20-hydroxyecdysone on diapause termination in eggs of Locusta migratoria and Oxya yezoensis. Journal of Insect Physiology 52: 473–479.
[36]
Pierceall WE, Li C, Biran A, Miura K, Raikhel AS, et al. (1999) E75 expression in mosquito ovary and fat body suggests reiterative use of ecdysone-regulated hierarchies in development and reproduction. Molecular and Cellular Endocrinology 150: 73–89.
[37]
Henrich VC, Livingston L, Gilbert LI (1993) Developmental requirements for the ecdysoneless (ecd) locus in Drosophila melanogaster. Developmental Genetics 14: 369–377.
[38]
Li T, Bender M (2000) A conditional rescue system reveals essential functions for the ecdysone receptor (EcR) gene during molting and metamorphosis in Drosophila. Development 127: 2897–2905.
[39]
Riddiford LM (2008) Juvenile hormone action: a 2007 perspective. Journal of Insect Physiology 54: 895–901.
[40]
Helvig C, Koener J, Unnithan G, Feyereisen R (2004) CYP15A1, the cytochrome P450 that catalyzes epoxidation of methyl farnesoate to juvenile hormone III in cockroach corpora allata. Proceedings of the National Academy of Sciences of the United States of America 101: 4024–4029.
[41]
Vannini L, Ciolfi S, Dallai R, Frati F, Hoffmann KH, et al. (2010) Putative-farnesoic acid O-methyltransferase (FAMeT) in medfly reproduction. Archives of Insect Biochemistry and Physiology 75: 92–106.
[42]
Tarrant AM, Baumgartner MF, Verslycke T, Johnson CL (2008) Differential gene expression in diapausing and active Calanus finmarchicus (Copepoda). Marine Ecology Progress Series 355: 193–207.
[43]
Denlinger DL (2002) Regulation of diapause. Annual Review of Entomology 47: 93–122.
[44]
Aruda AM, Baumgartner MF, Reitzel AM, Tarrant AM (2011) Heat shock protein expression during stress and diapause in the marine copepod Calanus finmarchicus. Journal of Insect Physiology 57: 665–675.
[45]
Jakobsson A, Westerberg R, Jacobsson A (2006) Fatty acid elongases in mammals: their regulation and roles in metabolism. Progress in Lipid Research 45: 237–249.
[46]
Tillman JA, Seybold SJ, Jurenka RA, Blomquist GJ (1999) Insect pheromones – an overview of biosynthesis and endocrine regulation. Insect Biochemistry and Molecular biology 29: 481–514.
[47]
Wang S, (2009) Reproduction, population recruitment and life history of Calanus sinicus in the Yellow Sea (Dissertation, in Chinese with English abstract). Qingdao: Institute of Oceanology, Chinese Academy of Sciences.
[48]
Liu Q, Yuan Y, Lin J, Zhong Y (2010) Advance of researches on insect cuticular proteins and the regulation mechanism of their gene expression (in Chinese with English abstract). Chinese Journal of Applied Entomology 002: 247–255.
[49]
Eichner C, Frost P, Dysvik B, Jonassen I, Kristiansen B, et al. (2008) Salmon louse (Lepeophtheirus salmonis) transcriptomes during post molting maturation and egg production, revealed using EST-sequencing and microarray analysis. BMC Genomics 9: 126.
[50]
Harrison PM, Arosio P (1996) The ferritins: molecular properties, iron storage function and cellular regulation. Biochimica et Biophysica Acta 1275: 161–203.
[51]
Reif DW (1992) Ferritin as a source of iron for oxidative damage. Free Radical Biology & Medicine 12: 417–427.
[52]
MacRae TH (2010) Gene expression, metabolic regulation and stress tolerance during diapause. Cellular and Molecular Life Sciences 67: 2405–2424.
[53]
Vos MJ, Hageman J, Carra S, Kampinga HH (2008) Structural and Functional Diversities between Members of the Human HSPB, HSPH, HSPA, and DNAJ Chaperone Families. Biochemistry 47: 7001–7011.
[54]
Wang M (2010) Application of molecular markers to the researches on pelagic copepods in the Chinese coastal regions (Dissertation, in Chinese with English abstract). Qingdao: Institute of Oceanology, Chinese Academy of Sciences.
[55]
Chen Z, Xue C, Zhu S, Zhou F, Ling XB, et al. (2005) GoPipe: streamlined gene ontology annotation for batch anonymous sequences with statistics. Progress in Biochemistry and Biophysics 32: 187–191.
[56]
Colbourne JK, Pfrender ME, Gilbert D, Thomas WK, Tucker A, et al. (2011) The Ecoresponsive Genome of Daphnia pulex. Science 331: 555–561.
[57]
Kanehisa M, Goto S, Furumichi M, Tanabe M, Hirakawa M (2010) KEGG for representation and analysis of molecular networks involving diseases and drugs. Nucleic Acids Research 38: 355–360.
[58]
Mortazavi A, Williams BA, McCue K, Schaeffer L, Wold B (2008) Mapping and quantifying mammalian transcriptomes by RNA-Seq. Nature Methods 5: 621–628.
[59]
Wang L, Feng Z, Wang X, Zhang X (2010) DEGseq: an R package for identifying differentially expressed genes from RNA-seq data. Bioinformatics 26: 136–138.
