A comprehensive coding and microRNA transcriptome of vertebral bone in postlarvae and juveniles of Senegalese sole (Solea senegalensis).

Understanding vertebral bone development is essential to prevent skeletal malformations in farmed fish related to genetic and environmental factors. This is an important issue in Solea senegalensis, with special impact of spinal anomalies in postlarval and juvenile stages. Vertebral bone transcriptomics in farmed fish mainly comes from coding genes, and barely on miRNA expression. Here, we used RNA-seq of spinal samples to obtain the first comprehensive coding and miRNA transcriptomic repertoire for postlarval and juvenile vertebral bone, covering different vertebral phenotypes and egg-incubation temperatures related to skeleton health in S. senegalensis. Coding genes, miRNA and pathways regulating bone development and growth were identified. Differential transcriptomic profiles and suggestive mRNA-miRNA interactions were found between postlarvae and juveniles. Bone-related genes and functions were associated with the extracellular matrix, development and regulatory processes, calcium binding, retinol and lipid metabolism or response to stimulus, including those revealed by the miRNA targets related to signaling, cellular and metabolic processes, growth, cell proliferation and biological adhesion. Pathway enrichment associated with fish skeleton were identified when comparing postlarvae and juveniles: growth and bone development functions in postlarvae, while actin cytoskeleton, focal adhesion and proteasome related to bone remodeling in juveniles. The transcriptome data disclosed candidate coding and miRNA gene markers related to bone cell processes, references for functional studies of the anosteocytic bone of S. senegalensis. This study establishes a broad transcriptomic foundation to study healthy and anomalous spines under early thermal conditions across life-stages in S. senegalensis, and for comparative analysis of skeleton homeostasis and pathology in fish and vertebrates.

De novo transcriptome assembly of Perkinsus olseni trophozoite stimulated in vitro with Manila clam (Ruditapes philippinarum) plasma.

The protistan parasite Perkinsus olseni is a deadly causative agent of perkinsosis, a molluscan disease affecting Manila clam (Ruditapes philippinarum), having a significant impact on world mollusc production. Deciphering the underlying molecular mechanisms in R. philippinarum-P. olseni interaction is crucial for controlling this parasitosis. The present study investigated the transcriptional expression in the parasite trophozoite using RNA-seq. Control and treatment (in vitro challenged with Manila clam-plasma) P. olseni trophozoite RNA were extracted and sequenced on the Illumina HiSeq 2000 instrument using a 100-bp paired-end sequencing strategy. Paired reads (64.7 million) were de novo assembled using Trinity, and the resultant transcripts were further clustered using CAP3. The re-constructed P. olseni transcriptome contains 47,590 unique transcripts of which 23,505 were annotated to 9764 unique proteins. A large number of genes were associated with Gene Ontology terms such as stress and immune-response, cell homeostasis, antioxidation, cell communication, signal transduction, signalling and proteolysis. Among annotated transcripts, a preliminary gene expression analysis detected 679 up-regulated and 478 down-regulated genes, linked to virulence factors, anti-oxidants, adhesion and immune-response molecules. Genes of several metabolic pathways such as DOXP/MEP, FAS II or folate biosynthesis, which are potential therapeutic targets, were identified. This study is the first description of the P. olseni transcriptome, and provides a substantial genomic resource for studying the molecular mechanisms of the host-parasite interaction in perkinsosis. In this sense, it is also the first evaluation of the parasite gene expression after challenge with clam extracellular products.

Expressed sequence tags (ESTs) from immune tissues of turbot (Scophthalmus maximus) challenged with pathogens.

BACKGROUND: The turbot (Scophthalmus maximus; Scophthalmidae; Pleuronectiformes) is a flatfish species of great relevance for marine aquaculture in Europe. In contrast to other cultured flatfish, very few genomic resources are available in this species. Aeromonas salmonicida and Philasterides dicentrarchi are two pathogens that affect turbot culture causing serious economic losses to the turbot industry. Little is known about the molecular mechanisms for disease resistance and host-pathogen interactions in this species. In this work, thousands of ESTs for functional genomic studies and potential markers linked to ESTs for mapping (microsatellites and single nucleotide polymorphisms (SNPs)) are provided. This information enabled us to obtain a preliminary view of regulated genes in response to these pathogens and it constitutes the basis for subsequent and more accurate microarray analysis. RESULTS: A total of 12584 cDNAs partially sequenced from three different cDNA libraries of turbot (Scophthalmus maximus) infected with Aeromonas salmonicida, Philasterides dicentrarchi and from healthy fish were analyzed. Three immune-relevant tissues (liver, spleen and head kidney) were sampled at several time points in the infection process for library construction. The sequences were processed into 9256 high-quality sequences, which constituted the source for the turbot EST database. Clustering and assembly of these sequences, revealed 3482 different putative transcripts, 1073 contigs and 2409 singletons. BLAST searches with public databases detected significant similarity (e-value < or = 1e-5) in 1766 (50.7%) sequences and 816 of them (23.4%) could be functionally annotated. Two hundred three of these genes (24.9%), encoding for defence/immune-related proteins, were mostly identified for the first time in turbot. Some ESTs showed significant differences in the number of transcripts when comparing the three libraries, suggesting regulation in response to these pathogens. A total of 191 microsatellites, with 104 having sufficient flanking sequences for primer design, and 1158 putative SNPs were identified from these EST resources in turbot. CONCLUSION: A collection of 9256 high-quality ESTs was generated representing 3482 unique turbot sequences. A large proportion of defence/immune-related genes were identified, many of them regulated in response to specific pathogens. Putative microsatellites and SNPs were identified. These genome resources constitute the basis to develop a microarray for functional genomics studies and marker validation for genetic linkage and QTL analysis in turbot.

Design and performance of a turbot (Scophthalmus maximus) oligo-microarray based on ESTs from immune tissues.

An expressed sequence tag database from immune tissues was used to design the first high-density turbot (Scophthalmus maximus) oligo-microarray with the aim of identifying candidate genes for tolerance to pathogens. Specific oligonucleotides (60 mers) were successfully designed for 2,716 out of 3,482 unique sequences of the database. An Agilent custom oligo-microarray 8 x 15 k (five replicates/gene; eight microarrays/slide) was constructed. The performance of the microarray and the sources of variation along microarray analysis were examined on spleen pools of controls and Aeromonas salmonicida-challenged fish at 3 days postinfection. Only 48 out of 2,716 probes did not show signal of hybridization on the 32 microarrays employed, thus demonstrating the consistency of the bioinformatic applications of our database. An asymmetric hierarchical design was employed to ascertain the noise associated with biological and technical (RNA extraction, labeling, hybridization, slide, and dye bias) factors using 1C and 2C labeling approaches. The high correlation coefficient between replicates at most factors tested demonstrated the high reproducibility of the signal. An analysis of random-effects variance revealed that technical variation was mostly negligible, and biological variation represented the main factor, even using pooled samples. One-color approach performed at least as well as 2C, suggesting their usefulness due to its higher design flexibility and lower cost. A relevant proportion of genes turn out to be differentially labeled depending on fluorophore, which alerts for the likely need of swapping replication in 2C experiments. A set of differentially expressed genes and enriched functions related to immune/defense response were detected at 3 days postchallenging.