RNA sequencing and de novo assembly of the digestive gland transcriptome in Mytilus galloprovincialis fed with toxinogenic and non-toxic strains of Alexandrium minutum

BMC Research Notes

Alberto Pallavicini²  Paola Venier⁴  Alfred Beran³  Elena Riccardi¹  Anna Milandri¹  Chiara Manfrin²  Gianluca De Moro²  Marco Gerdol² 

[1] Fondazione Centro Ricerche Marine, viale Amerigo Vespucci 2, Cesenatico, Forlì-Cesena 47042, Italy;Laboratory of Genetics, Department of Life Sciences, University of Trieste, Via Licio Giorgeri 5, Trieste 34126, Italy;Istituto nazionale di Oceanografia e di Geofisica sperimentale, Dipartimento di Oceanografia Biologica, via Auguste Piccard 54, Santa Croce, Trieste 34151, Italy;Department of Biology, University of Padua, Via Bassi 58 / B 35121 Padova (PD) Padua, Italy

关键词: transcriptome;    RNA-seq;    Alexandrium minutum;    De novo assembly;    Paralytic shellfish poisoning;    Mytilus galloprovincialis;   

Others  :  1127263
DOI  :  10.1186/1756-0500-7-722

 received in 2013-11-14, accepted in 2014-10-02,  发布年份 2014

【 摘 要 】

Background

The Mediterranean mussel Mytilus galloprovincialis is marine bivalve with a relevant commercial importance as well as a key sentinel organism for the biomonitoring of environmental pollution. Here we report the RNA sequencing of the mussel digestive gland, performed with the aim: a) to produce a high quality de novo transcriptome assembly, thus improving the genetic and molecular knowledge of this organism b) to provide an initial assessment of the response to paralytic shellfish poisoning (PSP) on a molecular level, in order to identify possible molecular markers of toxin accumulation.

Results

The comprehensive de novo assembly and annotation of the transcriptome yielded a collection of 12,079 non-redundant consensus sequences with an average length of 958 bp, with a high percentage of full-length transcripts. The whole-transcriptome gene expression study indicated that the accumulation of paralytic toxins produced by the dinoflagellate Alexandrium minutum over a time span of 5 days scarcely affected gene expression, but the results need further validation with a greater number of biological samples and naturally contaminated specimens.

Conclusion

The digestive gland reference transcriptome we produced significantly improves the data collected from previous sequencing efforts and provides a basic resource for expanding functional genomics investigations in M. galloprovincialis. Although not conclusive, the results of the RNA-seq gene expression analysis support the classification of mussels as bivalves refractory to paralytic shellfish poisoning and point out that the identification molecular biomarkers of PSP in the digestive gland of this organism is problematic.

【 授权许可】

   
2014 Gerdol et al.; licensee BioMed Central Ltd.

【 图 表 】

Figure 1.

Figure 2.

Figure 3.

【 参考文献 】

• [1]Pérez-Enciso M, Ferretti L: Massive parallel sequencing in animal genetics: Wherefroms and wheretos. Anim Genet 2010, 41:561-569.
• [2]Suárez-Ulloa V, Fernández-Tajes J, Manfrin C, Gerdol M, Venier P, Eirín-López JM: Bivalve omics: State of the art and potential applications for the biomonitoring of harmful marine compounds. Mar Drugs 2013, 11:4370-4389.
• [3]Smith SA, Wilson NG, Goetz FE, Feehery C, Andrade SCS, Rouse GW, Giribet G, Dunn CW: Resolving the evolutionary relationships of molluscs with phylogenomic tools. Nature 2011, 480:364-367.
• [4]Kocot KM, Cannon JT, Todt C, Citarella MR, Kohn AB, Meyer A, Santos SR, Schander C, Moroz LL, Lieb B, Halanych KM: Phylogenomics reveals deep molluscan relationships. Nature 2011, 477:452-456.
• [5]Takeuchi T, Kawashima T, Koyanagi R, Gyoja F, Tanaka M, Ikuta T, Shoguchi E, Fujiwara M, Shinzato C, Hisata K, Fujie M, Usami T, Nagai K, Maeyama K, Okamoto K, Aoki H, Ishikawa T, Masaoka T, Fujiwara A, Endo K, Endo H, Nagasawa H, Kinoshita S, Asakawa S, Watabe S, Satoh N: Draft genome of the pearl oyster Pinctada fucata: a platform for understanding bivalve biology. DNA Res 2012, 19:117-130.
• [6]Zhang G, Fang X, Guo X, Li L, Luo R, Xu F, Yang P, Zhang L, Wang X, Qi H, Xiong Z, Que H, Xie Y, Holland PWH, Paps J, Zhu Y, Wu F, Chen Y, Wang J, Peng C, Meng J, Yang L, Liu J, Wen B, Zhang N, Huang Z, Zhu Q, Feng Y, Mount A, Hedgecock D, et al.: The oyster genome reveals stress adaptation and complexity of shell formation. Nature 2012, 490:49-54.
• [7]Venier P, De Pitta C, Bernante F, Varotto L, De Nardi B, Bovo G, Roch P, Novoa B, Figueras A, Pallavicini A, Lanfranchi G: MytiBase: a knowledgebase of mussel (M. galloprovincialis) transcribed sequences. BMC Genomics 2009, 10:72. BioMed Central Full Text
• [8]Craft JA, Gilbert JA, Temperton B, Dempsey KE, Ashelford K, Tiwari B, Hutchinson TH, Chipman JK: Pyrosequencing of Mytilus galloprovincialis cDNAs: tissue-specific expression patterns. PLoS One 2010, 5:e8875.
• [9]Philipp EER, Kraemer L, Melzner F, Poustka AJ, Thieme S, Findeisen U, Schreiber S, Rosenstiel P: Massively parallel RNA sequencing identifies a complex immune gene repertoire in the lophotrochozoan Mytilus edulis. PLoS One 2012, 7:e33091.
• [10]Suárez-Ulloa V, Fernández-Tajes J, Aguiar-Pulido V, Rivera-Casas C, González-Romero R, Ausio J, Méndez J, Dorado J, Eirín-López J: The CHROMEVALOA database: a resource for the vvaluation of okadaic acid contamination in the marine environment based on the chromatin-associated transcriptome of the mussel Mytilus galloprovincialis. Mar Drugs 2013, 11:830-841.
• [11]Tanguy M, McKenna P, Gauthier-Clerc S, Pellerin J, Danger J-M, Siah A: Sequence analysis of a normalized cDNA library of Mytilus edulis hemocytes exposed to Vibrio splendidus LGP32 strain. Results Immunol 2013, 3:40-50.
• [12]Freer A, Bridgett S, Jiang J, Cusack M: Biomineral proteins from Mytilus edulis mantle tissue transcriptome. Marine Biotechnol 2013, 16:1-12.
• [13]Wang Z, Gerstein M, Snyder M: RNA-Seq: a revolutionary tool for transcriptomics. Nat Rev Genet 2009, 10:57-63.
• [14]Domeneghetti S, Manfrin C, Varotto L, Rosani U, Gerdol M, De Moro G, Pallavicini A, Venier P: How gene expression profiles disclose vital processes and immune responses in Mytilus spp. Invertebrate Surviv J 2011, 8:179-189.
• [15]Moreira R, Balseiro P, Planas JV, Fuste B, Beltran S, Novoa B, Figueras A: Transcriptomics of in vitro immune-stimulated hemocytes from the Manila clam Ruditapes philippinarum using high-throughput sequencing. PLoS One 2012, 7:e35009.
• [16]Pauletto M, Milan M, Moreira R, Novoa B, Figueras A, Babbucci M, Patarnello T, Bargelloni L: Deep transcriptome sequencing of Pecten maximus hemocytes: A genomic resource for bivalve immunology. Fish Shellfish Immunol 2014, 37:154-165.
• [17]Ghiselli F, Milani L, Chang PL, Hedgecock D, Davis JP, Nuzhdin SV, Passamonti M: De novo assembly of the Manila clam Ruditapes philippinarum transcriptome provides new insights into expression bias, mitochondrial doubly uniparental inheritance and sex determination. Mol Biol Evol 2012, 29:771-786.
• [18]de Sousa JT, Milan M, Bargelloni L, Pauletto M, Matias D, Joaquim S, Matias AM, Quillien V, Leitão A, Huvet A: A microarray-based analysis of gametogenesis in two Portuguese populations of the European clam Ruditapes decussatus. PLoS One 2014, 9:e92202.
• [19]Hou R, Bao Z, Wang S, Su H, Li Y, Du H, Hu J, Wang S, Hu X: Transcriptome sequencing and de novo analysis for Yesso scallop (Patinopecten yessoensis) using 454 GS FLX. PLoS One 2011, 6:e21560.
• [20]Clark M, Thorne M, Vieira F, Cardoso J, Power D, Peck L: Insights into shell deposition in the Antarctic bivalve Laternula elliptica: gene discovery in the mantle transcriptome using 454 pyrosequencing. BMC Genomics 2010, 11:362. BioMed Central Full Text
• [21]Hallegraeff GM, Steffensen DA, Wetherbee R: Three estuarine Australian dinoflagellates that can produce paralytic shellfish toxins. J Plankton Res 1988, 10:533.
• [22]James KJ, Carey B, O'Halloran J, Van Pelt FNAM, Škrabáková Z: Shellfish toxicity: Human health implications of marine algal toxins. Epidemiol Infect 2010, 138:927-940.
• [23]Terlau H, Heinemann SH, Stuhmer W, Pusch M, Conti F, Imoto K, Numa S: Mapping the site of block by tetrodotoxin and saxitoxin of sodium channel II. FEBS Lett 1991, 293:93-96.
• [24]Manfrin C, De Moro G, Torboli V, Venier P, Pallavicini A, Gerdol M: Physiological and molecular responses of bivalves to toxic dinoflagellates. Invertebrate Surviv J 2012, 9:184-199.
• [25]Conte FS: Economic impact of paralytic shellfish poison on the oyster industry in the Pacific United States. Aquaculture 1984, 39:331-343.
• [26]Anderson DM, Sullivan JJ, Reguera B: Paralytic shellfish poisoning in northwest Spain: The toxicity of the dinoflagellate Gymnodinium catenatum. Toxicon 1989, 27:665-674.
• [27]Coulson JC, Potts GR, Deans IR, Fraser SM: Dinoflagellate crop in the North Sea: Mortality of shags and other Sea Birds caused by paralytic shellfish poison. Nature 1968, 220:23-24.
• [28]Geraci JR: Humpback whales (Megaptera novaeangliae) fatally poisoned by dinoflagellate toxin. Can J Fish Aquat Sci 1989, 46:1895-1898.
• [29]Kvitek RG, Degange AR, Beitler MK: Paralytic shellfish poisoning toxins mediate feeding behavior of sea otters. Limnology & Oceanography 1991, 36:393-404.
• [30]Cembella AD, Quilliam MA, Lewis NI, Bauder AG, Dell'Aversano C, Thomas K, Jellett J, Cusack RR: The toxigenic marine dinoflagellate Alexandrium tamarense as the probable cause of mortality of caged salmon in Nova Scotia. Harmful Algae 2002, 1:313-325.
• [31]Bricelj VM, Lee JH, Cembella AD: Influence of dinoflagellate cell toxicity on uptake and loss of paralytic shellfish toxins in the northern quahog Mercenaria mercenaria. Mar Ecol Prog Ser 1991, 74:33-46.
• [32]MacQuarrie SP, Bricelj VM: Behavioral and physiological responses to PSP toxins in Mya arenaria populations in relation to previous exposure to red tides. Mar Ecol Prog Ser 2008, 366:59-74.
• [33]Twarog BM, Hidaka T, Yamaguchi H: Resistance to tetrodotoxin and saxitoxin in nerves of bivalve molluscs.A possible correlation with paralytic shellfish poisoning. Toxicon 1972, 10:273-278.
• [34]Twarog BM: Immunity to paralytic shellfish toxin in bivalve molluscs. In Proc Int Symp Coral Reefs. 2nd edition. Edited by Cameron AM, Campbell BM, Cribb AB. Brisbane, Australia: Great Barrier Reef Comm; 1974:505-512.
• [35]Bricelj VM, Lee JH, Cembella AD, Anderson DM: Uptake kinetics of paralytic shellfish toxins from the dinoflagellate Alexandrium fundyense in the mussel Mytilus edulis. Mar Ecol Prog Ser 1990, 63:117-188.
• [36]Shumway SE, Cucci TL: The effects of the toxic dinoflagellate Protogonyaulax tamarensis on the feeding and behaviour of bivalve molluscs. Aquat Toxicol 1987, 10:9-27.
• [37]Shumway SE, Pierce FC, Knowlton K: The effect of Protogonyaulax tamarensis on byssus production in Mytilus edulis l., Modiolus modiolus linnaeus, 1758 and Geukensia demissa dillwyn. Comp Biochem Physiol A Physiol 1987, 87:1021-1023.
• [38]Galimany E, Sunila I, Hégaret H, Ramón M, Wikfors GH: Experimental exposure of the blue mussel (Mytilus edulis, L.) to the toxic dinoflagellate Alexandrium fundyense: Histopathology, immune responses, and recovery. Harmful Algae 2008, 7:702-711.
• [39]Núñez-Acuña G, Aballay AE, Hégaret H, Astuya AP, Gallardo-Escárate C: Transcriptional responses of Mytilus chilensis exposed in vivo to Saxitoxin (STX). J Molluscan Stud 2013, 79:323-331.
• [40]Mat AM, Haberkorn H, Bourdineaud J-P, Massabuau J-C, Tran D: Genetic and genotoxic impacts in the oyster Crassostrea gigas exposed to the harmful alga Alexandrium minutum. Aquat Toxicol 2013, 140–141:458-465.
• [41]Manfrin C, Dreos R, Battistella S, Beran A, Gerdol M, Varotto L, Lanfranchi G, Venier P, Pallavicini A: Mediterranean mussel gene expression profile induced by okadaic acid exposure. Environ Sci Tech 2010, 44:8276-8283.
• [42]Ewen-Campen B, Shaner N, Panfilio K, Suzuki Y, Roth S, Extavour C: The maternal and early embryonic transcriptome of the milkweed bug Oncopeltus fasciatus. BMC Genomics 2011, 12:61. BioMed Central Full Text
• [43]Stüken A, Orr RJS, Kellmann R, Murray SA, Neilan BA, Jakobsen KS: Discovery of nuclear-encoded genes for the neurotoxin saxitoxin in dinoflagellates. PLoS One 2011, 6:e20096.
• [44]Negri AP, Jones GJ: Bioaccumulation of paralytic shellfish poisoning (PSP) toxins from the cyanobacterium Anabaena circinalis by the freshwater mussel Alathyria condola. Toxicon 1995, 33:667-678.
• [45]Pereira P, Dias E, Franca S, Pereira E, Carolino M, Vasconcelos V: Accumulation and depuration of cyanobacterial paralytic shellfish toxins by the freshwater mussel Anodonta cygnea. Aquat Toxicol 2004, 68:339-350.
• [46]Estrada N, Lagos N, García C, Maeda-Martínez A, Ascencio F: Effects of the toxic dinoflagellate Gymnodinium catenatum on uptake and fate of paralytic shellfish poisons in the Pacific giant lions-paw scallop Nodipecten subnodosus. Mar Biol 2007, 151:1205-1214.
• [47]Rodríguez-Juíz AM, Torrado M, Méndez J: Genome-size variation in bivalve molluscs determined by flow cytometry. Mar Biol 1996, 126:489-497.
• [48]Ieyama H, Kameoka O, Tan T, Yamasaki J: Chromosomes and nuclear DNA contents of some species of Mytilidae. Japanese Journal of Malacology 1994, 53:327-331.
• [49]Pallavicini A, del Mar CM, Gestal C, Dreos R, Figueras A, Venier P, Novoa B: High sequence variability of myticin transcripts in hemocytes of immune-stimulated mussels suggests ancient host-pathogen interactions. Dev Comp Immunol 2008, 32:213-226.
• [50]Rosani U, Varotto L, Rossi A, Roch P, Novoa B, Figueras A, Pallavicini A, Venier P: Massively parallel amplicon sequencing reveals isotype-specific variability of antimicrobial peptide transcripts in Mytilus galloprovincialis. PLoS One 2011, 6:e26680.
• [51]Mosquera E, Lopez JL, Alvarez G: Genetic variability of the marine mussel Mytilus galloprovincialis assessed using two-dimensional electrophoresis. Heredity 2003, 90:432-442.
• [52]Jackson D, Ellemor N, Degnan B: Correlating gene expression with larval competence, and the effect of age and parentage on metamorphosis in the tropical abalone Haliotis asinina. Mar Biol 2005, 147:681-697.
• [53]Gerdol M, Manfrin C, De Moro G, Figueras A, Novoa B, Venier P, Pallavicini A: The C1q domain containing proteins of the Mediterranean mussel Mytilus galloprovincialis: A widespread and diverse family of immune-related molecules. Dev Comp Immunol 2011, 35:635-643.
• [54]Gorbushin AM, Iakovleva NV: A new gene family of single fibrinogen domain lectins in Mytilus. Fish Shellfish Immunol 2011, 30:434-438.
• [55]Kong L, Zhang Y, Ye Z-Q, Liu X-Q, Zhao S-Q, Wei L, Gao G: CPC: assess the protein-coding potential of transcripts using sequence features and support vector machine. Nucleic Acids Res 2007, 35:W345-W349.
• [56]Ilott NE, Ponting CP: Predicting long non-coding RNAs using RNA sequencing. Methods 2013, 63:50-59.
• [57]Milan M, Coppe A, Reinhardt R, Cancela L, Leite R, Saavedra C, Ciofi C, Chelazzi G, Patarnello T, Bortoluzzi S, Bargelloni L: Transcriptome sequencing and microarray development for the Manila clam. Ruditapes philippinarum: genomic tools for environmental monitoring. BMC Genomics 2011, 12:234. BioMed Central Full Text
• [58]Delgado M, Estrada M, Camp J, Fernández JV, Santmartí M, Lletí C: Development of a toxic Alexandrium minutum Halim (Dinophyceae) bloom in the harbour of Sant Carles de la Ràpita (Ebro Delta, northwestern Mediterranean). Scientia Marina 1990, 54:1-7.
• [59]Maguer JF, Wafar M, Madec C, Morin P, Denn EEL: Nitrogen and phosphorus requirements of an Alexandrium minutum bloom in the Penzé Estuary, France. Limnol Oceanogr 2004, 49:1108-1114.
• [60]Garcés E, Bravo I, Vila M, Figueroa RI, Masó M, Sampedro N: Relationship between vegetative cells and cyst production during Alexandrium minutum bloom in Arenys de Mar harbour (NW Mediterranean). J Plankton Res 2004, 26:637-645.
• [61]Galluzzi L, Penna A, Bertozzini E, Vila M, Garcés E, Magnani M: Development of a real-time PCR assay for rapid detection and quantification of Alexandrium minutum (a dinoflagellate). Appl Environ Microbiol 2004, 70:1199-1206.
• [62]Van Lenning K, Vila M, Masó M, Garcés E, Anglès S, Sampedro N, Morales-Blake A, Camp J: Short-term variations in development of a recurrent toxic Alexandrium minutum-dominated dinoflagellate bloom induced by meteorological conditions. J Phycol 2007, 43:892-907.
• [63]Blanco J, Morono A, Franco J, Reyero MI: PSP detoxification kinetics in the mussel Mytilus galloprovincialis. One- and two-compartment models and the effect of some environmental variables. Mar Ecol Prog Ser 1997, 158:165-175.
• [64]Blanco J, Reyero MI, Franco J: Kinetics of accumulation and transformation of paralytic shellfish toxins in the blue mussel Mytilus galloprovincialis. Toxicon 2003, 42:777-784.
• [65]Fernández-Reiriz MJ, Navarro JM, Contreras AM, Labarta U: Trophic interactions between the toxic dinoflagellate Alexandrium catenella and Mytilus chilensis: Feeding and digestive behaviour to long-term exposure. Aquat Toxicol 2008, 87:245-251.
• [66]Botelho MJ, Vale C, Mota AM, Simões S, Gonalves MDL: Depuration kinetics of paralytic shellfish toxins in Mytilus galloprovincialis exposed to Gymnodinium catenatum: Laboratory and field experiments. J Environ Monit 2010, 12:2269-2275.
• [67]Fistarol GO, Legrand C, Selander E, Hummert C, Stolte W, Granéli E: Allelopathy in Alexandrium spp.: Effect on a natural plankton community and on algal monocultures. Aquat Microb Ecol 2004, 35:45-56.
• [68]Lelong A, Haberkorn H, Le Goïc N, Hégaret H, Soudant P: A new insight into allelopathic effects of Alexandrium minutum on photosynthesis and respiration of the diatom Chaetoceros neogracile revealed by photosynthetic-performance analysis and flow cytometry. Microb Ecol 2011, 62:919-930.
• [69]Tillmann U, John U, Cembella A: On the allelochemical potency of the marine dinoflagellate Alexandrium ostenfeldii against heterotrophic and autotrophic protists. J Plankton Res 2007, 29:527-543.
• [70]Juhl AR, Martins CA, Anderson DM: Toxicity of Alexandrium lusitanicum to gastropod larvae is not caused by paralytic-shellfish-poisoning toxins. Harmful Algae 2008, 7:567-573.
• [71]Mello DF, Silva PM, Barracco MA, Soudant P, Hégaret H: Effects of the dinoflagellate Alexandrium minutum and its toxin (saxitoxin) on the functional activity and gene expression of Crassostrea gigas hemocytes. Harmful Algae 2013, 26:45-51.
• [72]Ford SE, Bricelj VM, Lambert C, Paillard C: Deleterious effects of a nonPST bioactive compound(s) from Alexandrium tamarense on bivalve hemocytes. Mar Biol 2008, 154:241-253.
• [73]Gerdol M, De Moro G, Manfrin C, Venier P, Pallavicini A: Big defensins and mytimacins, new AMP families of the Mediterranean mussel Mytilus galloprovincialis. Dev Comp Immunol 2012, 36:390-399.
• [74]Toubiana M, Gerdol M, Rosani U, Pallavicini A, Venier P, Roch P: Toll-like receptors and MyD88 adaptors in Mytilus: Complete cds and gene expression levels. Dev Comp Immunol 2013, 40:158-166.
• [75]Connell LB, MacQuarrie SP, Twarog BM, Iszard M, Bricelj VM: Population differences in nerve resistance to paralytic shellfish toxins in softshell clam, Mya arenaria, associated with sodium channel mutations. Mar Biol 2007, 150:1227-1236.
• [76]Medhioub W, Ramondenc S, Vanhove AS, Vergnes A, Masseret E, Savar V, Amzil Z, Laabir M, Rolland JL: Exposure to the neurotoxic dinoflagellate, Alexandrium catenella, induces apoptosis of the hemocytes of the oyster, Crassostrea gigas. Mar Drugs 2013, 11:4799-4814.
• [77]EFSA: Scientific opinion of the panel on contaminants in the food chain on a request from the European Commission on marine biotoxins in shellfish – Saxitoxin group. EFSA Journal 2009, 1019:1-76.
• [78]Desbiens M, Cembella AD: Minimization of PSP toxin accumulation in cultured blue mussels (Mytilus edulis) by vertical displacement in the water column. In Toxic phytoplankton blooms in the sea. Edited by Smayda TJ, Shimizu YT. Amsterdam: Elsevier; 1993:395-400.
• [79]Giron Perez MI: Relationships between innate immunity in bivalve molluscs and environmental pollution. Invertebrate Surviv J 2010, 7:149-156.
• [80]Agatha S, StrÜDer-Kypke MC, Beran A: Morphologic and genetic variability in the marine planktonic ciliate Laboea strobila Lohmann, 1908 (Ciliophora, Oligotrichia), with notes on its ontogenesis. J Eukaryot Microbiol 2004, 51:267-281.
• [81]Lawrence JF, Niedzwiadek B, Menard C: Quantitative determination of paralytic shellfish poisoning toxins in shellfish using prechromatographic oxidation and liquid chromatography with fluorescence detection: interlaboratory study. J AOAC Int 2004, 87:83-100.
• [82]Quilliam MA, Janeček M, Lawrence JF: Characterization of the oxidation products of paralytic shellfish poisoning toxins by liquid chromatography/mass spectrometry. Rapid Commun Mass Spectrom 1993, 7:482-487.
• [83]Grabherr MG, Haas BJ, Yassour M, Levin JZ, Thompson DA, Amit I, Adiconis X, Fan L, Raychowdhury R, Zeng Q, Chen Z, Mauceli E, Hacohen N, Gnirke A, Rhind N, di Palma F, Birren BW, Nusbaum C, Lindblad-Toh K, Friedman N, Regev A: Full-length transcriptome assembly from RNA-Seq data without a reference genome. Nat Biotechnol 2011, 29:644-652.
• [84]Altschul SF, Madden TL, Schäffer AA, Zhang J, Zhang Z, Miller W, Lipman DJ: Gapped BLAST and PSI-BLAST: a new generation of protein database search programs. Nucleic Acids Res 1997, 25:3389-3402.
• [85]McIntyre L, Lopiano K, Morse A, Amin V, Oberg A, Young L, Nuzhdin S: RNA-seq: technical variability and sampling. BMC Genomics 2011, 12:293. BioMed Central Full Text
• [86]Altschul SF, Gish W, Miller W, Myers EW, Lipman DJ: Basic local alignment search tool. J Mol Biol 1990, 215:403-410.
• [87]Finn R, Clements J, Eddy SR: HMMER web server: interactive sequence similarity searching. Nucleic Acids Res 2011, 39:29-37.
• [88]Punta M, Coggill PC, Eberhardt RY, Mistry J, Tate J, Boursnell C, Pang N, Forslund K, Ceric G, Clements J, Heger A, Holm L, Sonnhammer ELL, Eddy SR, Bateman A, Finn RD: The Pfam protein families database. Nucleic Acids Res 2012, 40:D290-D301.
• [89]Hunter S, Apweiler R, Attwood TK, Bairoch A, Bateman A, Binns D, Bork P, Das U, Daugherty L, Duquenne L, Finn RD, Gough J, Haft D, Hulo N, Kahn D, Kelly E, Laugraud A, Letunic I, Lonsdale D, Lopez R, Madera M, Maslen J, McAnulla C, McDowall J, Mistry J, Mitchell A, Mulder N, Natale D, Orengo C, Quinn AF, et al.: InterPro: the integrative protein signature database. Nucleic Acids Res 2009, 37:D211-D215.
• [90]Powell S, Szklarczyk D, Trachana K, Roth A, Kuhn M, Muller J, Arnold R, Rattei T, Letunic I, Doerks T, Jensen LJ, von Mering C, Bork P: eggNOG v3.0: Orthologous groups covering 1133 organisms at 41 different taxonomic ranges. Nucleic Acids Res 2012, 40:D284-D289.
• [91]Ashburner M, Ball CA, Blake JA, Botstein D, Butler H, Cherry JM, Davis AP, Dolinski K, Dwight SS, Eppig JT, Harris MA, Hill DP, Issel-Tarver L, Kasarskis A, Lewis S, Matese JC, Richardson JE, Ringwald M, Rubin GM, Sherlock G: Gene ontology: Tool for the unification of biology. Nat Genet 2000, 25:25-29.
• [92]Kal AJ, Van Zonneveld AJ, Benes V, Van Den Berg M, Koerkamp MG, Albermann K, Strack N, Ruijter JM, Richter A, Dujon B, Ansorge W, Tabak HF: Dynamics of gene expression revealed by comparison of serial analysis of gene expression transcript profiles from yeast grown on two different carbon sources. Mol Biol Cell 1999, 10:1859-1872.