Deep sequencing for de novo construction of a marine fish (Sparus aurata) transcriptome database with a large coverage of protein-coding transcripts

[Background] The gilthead sea bream (Sparus aurata) is the main fish species cultured in the Mediterranean area and constitutes an interesting model of research. Nevertheless, transcriptomic and genomic data are still scarce for this highly valuable species. A transcriptome database was constructed by de novo assembly of gilthead sea bream sequences derived from public repositories of mRNA and collections of expressed sequence tags together with new high-quality reads from five cDNA 454 normalized libraries of skeletal muscle (1), intestine (1), head kidney (2) and blood (1).[Results] Sequencing of the new 454 normalized libraries produced 2,945,914 high-quality reads and the de novo global assembly yielded 125,263 unique sequences with an average length of 727 nt. Blast analysis directed to protein and nucleotide databases annotated 63,880 sequences encoding for 21,384 gene descriptions, that were curated for redundancies and frameshifting at the homopolymer regions of open reading frames, and hosted at http://www.nutrigroup-iats.org/seabreamdb. Among the annotated gene descriptions, 16,177 were mapped in the Ingenuity Pathway Analysis (IPA) database, and 10,899 were eligible for functional analysis with a representation in 341 out of 372 IPA canonical pathways. The high representation of randomly selected stickleback transcripts by Blast search in the nucleotide gilthead sea bream database evidenced its high coverage of protein-coding transcripts.[Conclusions] The newly assembled gilthead sea bream transcriptome represents a progress in genomic resources for this species, as it probably contains more than 75% of actively transcribed genes, constituting a valuable tool to assist studies on functional genomics and future genome projects.This work was funded under EU seventh Framework Programme by projects ARRAINA (Advanced Research Initiatives for Nutrition & Aquaculture, FP7/2007-2013; grant agreement nº 288925) and AQUAEXCEL (Aquaculture Infrastructures for Excellence in European Fish Research, FP7/2007-2012; grant agreement nº 262336). Additional funding was obtained by Generalitat Valenciana (research grant PROMETEO 2010/006) and Spanish Government through AQUAGENOMICS (Ingenio-2010 Programme), AQUAFAT (AGL2009-07797) and ENTEROMYXCONTROL (AGL2009-13282-C02-01) projects.We also acknowledge institutional support from the Unit of Information Resources for Research at the “Consejo Superior de Investigaciones Científicas” (CSIC) for the article-processing charges contribution.Peer Reviewe