De novo assembly, characterization and functional annotation of Senegalese sole (Solea senegalensis) and common sole (Solea solea) transcriptomes: integration in a database and design of a microarray

Assessing the role of ulvan as immunonutrient in Solea senegalensis

Immunonutrition is a promising and viable strategy for the development of prophylactic measures in aquaculture. Ulvan, a sulphated marine polysaccharide from green seaweeds, has many biological activities including the immunomodulatory ones. The aim of this study was to assess the short and long-term effects of an ulvan-rich extract obtained from U. ohnoi as immunonutrient in Senegalese sole juveniles. In this work, an ulvan-rich extract from Ulva ohnoi has been obtained by the hot water method and isolated by ethanol precipitation. The FTIR analysis revealed that the ulvan-rich extact had very similar characteristics to previously published ulvan spectra. The total sulfate and protein content was 24.85 ± 3.98 and 0.91 ± 0.04 %, respectively. In vitro assays performed in Senegalese sole (Solea senegalensis) macrophages showed that the ulvan obtained in this study did not compromise the cell viability at concentrations up to 1 mg ml-1 and expression levels of lyg, irf1, il6, il10, c7, tf and txn were significantly upregulated in a concentration dependent-manner. Finally, S. senegalensis juveniles were fed basal diets and diets supplemented with the ulvan-rich extract at ratios 1 and 2 % for 30 days and then, challenged with Photobacterium damselae subsp. piscicida (Phdp). Thereafter, ulvan was withdrawn from the diet and all juveniles were fed the basal diet for 30 days. At 30 days post withdrawal (dpw), juveniles were challenged with Phdp. The expression profiles of a set of genes related to the immune system in spleen were evaluated as well as the lysozyme, peroxidase and bactericidal activity in plasma. Dietary effects of 1 % ulvan resulted in a boost of the immune response and increased disease resistance at short-term whereas juveniles fed diets supplemented with 2 % ulvan showed a significant decrease in the bactericidal activity and lack of protection against Phdp. At long-term (30 days after the withdrawal of ulvan), an improved response was observed in juveniles previously fed 1 % ulvan.

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.

Photoreceptor distributions, visual pigments and the opsin repertoire of Atlantic halibut (Hippoglossus hippoglossus)

Fishes often have cone photoreceptors organized in lattice-like mosaic formations. In flatfishes, these lattices undergo dramatic changes during metamorphosis whereby a honeycomb mosaic of single cones in the larva is replaced by a square mosaic of single and double cones in the adult. The spatio-temporal dynamics of this transition are not well understood. Here, we describe the photoreceptors and mosaic formations that occur during the larva to juvenile transition of Atlantic halibut from the beginning of eye migration to its completion. To gauge the possibility of colour vision, visual pigments in juveniles were measured by microspectrophotometry and the opsin repertoire explored using bioinformatics. At the start of eye migration, the larva had a heterogeneous retina with honeycomb mosaic in the dorsonasal and ventrotemporal quadrants and a square mosaic in the ventronasal and dorsotemporal quadrants. By the end of metamorphosis, the square mosaic was present throughout the retina except in a centrodorsotemporal area where single, double and triple cones occurred randomly. Six cone visual pigments were found with maximum absorbance (λ_max, in nm) in the short [S(431) and S(457)], middle [M(500), M(514) and M(527)], and long [L(550)] wavelengths, and a rod visual pigment with λ_max at 491 nm. These pigments only partially matched the opsin repertoire detected by query of the Atlantic halibut genome. We conclude that the Atlantic halibut undergoes a complex re-organization of photoreceptors at metamorphosis resulting in a multi-mosaic retina adapted for a demersal life style.

Overview and New Insights Into the Diversity, Evolution, Role, and Regulation of Kisspeptins and Their Receptors in Teleost Fish

In the last two decades, kisspeptin (Kiss) has been identified as an important player in the regulation of reproduction and other physiological functions in vertebrates, including several fish species. To date, two ligands (Kiss1, Kiss2) and three kisspeptin receptors (Kissr1, Kissr2, Kissr3) have been identified in teleosts, likely due to whole-genome duplication and loss of genes that occurred early in teleost evolution. Recent results in zebrafish and medaka mutants have challenged the notion that the kisspeptin system is essential for reproduction in fish, in marked contrast to the situation in mammals. In this context, this review focuses on the role of kisspeptins at three levels of the reproductive, brain-pituitary-gonadal (BPG) axis in fish. In addition, this review compiled information on factors controlling the Kiss/Kissr system, such as photoperiod, temperature, nutritional status, sex steroids, neuropeptides, and others. In this article, we summarize the available information on the molecular diversity and evolution, tissue expression and neuroanatomical distribution, functional significance, signaling pathways, and gene regulation of Kiss and Kissr in teleost fishes. Of particular note are recent advances in understanding flatfish kisspeptin systems, which require further study to reveal their structural and functional diversity.

Evaluation of an oral DNA nanovaccine against photobacteriosis in Solea senegalensis

Infectious diseases are one of the main causes of social and economical losses in world aquaculture. Senegalese sole (Solea senegalensis) is an important species for aquaculture in southern Europe, whose production is affected by the appearance of bacterial diseases such as photobacteriosis, a septicemia caused by Photobacterium damselae subsp. piscicida (Phdp). The aim of this study was to obtain an oral DNA nanovaccine and to evaluate its efficacy against Phdp in S. senegalensis juveniles. For this purpose, the amplified product corresponding to the protein inosine-5'-monophophate dehydrogenase (IMPDH) from Phdp, was cloned into the expression vector pcDNA™6.2/C-EmGFP-GW obtaining the DNA vaccine named as pPDPimpdh. The correct transcription and protein expression was verified at 48 h post tansfection in HEK293 cells. Chitosan nanoparticles (CS-TPP NPs) were prepared by ionotropic gelation and their features were appropriate for use as oral delivery system. Therefore, pPDPimpdh was protected with chitosan CS-TPP NPs throughout complex coacervation method giving as a result a DNA nanovaccine referred as CS-TPP+pPDPimpdh NPs. Sole juveniles were vaccinated orally with CS-TPP NPs, pPDPimpdh and CS-TPP+pPDPimpdh NPs followed by a challenge with Phdp at 30 days post vaccination (dpv). The relative percentage survival (RPS) for pPDPimpdh vaccinated groups was 6.25%, probably due to its degradation in the digestive tract. RPS value obtained for CS-TPP NPs and CS-TPP+pPDPimpdh NPs was 40% and antibodies were observed in both cases. However, a delay in mortality was observed in sole juveniles vaccinated orally with CS-TPP+pPDPimpdh NPs. In fact, an upregulation of tf, mhcII, cd8a and igm in the posterior gut and c3, hamp1, tf and cd4 in spleen was observed in juveniles vaccinated with CS-TPP+pPDPimpdh NPs. After challenge, a modulation of cd8a and cd4 expression levels in the posterior gut and c3, tf, lyg, cd4, igm and igt expression levels in spleen was observed. Moreover, the concentration of lysozyme in skin mucus significantly increased in fish vaccinated orally with CS-TPP+pPDPimpdh NPs at 11 dpc. These data indicate that oral vaccination with CS-TPP+pPDPimpdh NPs could be acting through the non-specific immune responses as well as the specific humoral and cell mediated immunity and provide the first step toward a development of an oral DNA nanovaccine against Phdp in sole.

A Comprehensive Integrated Genetic Map of the Complete Karyotype of Solea senegalensis (Kaup 1858)

Solea senegalensis aquaculture production has experienced a great increase in the last decade and, consequently, the genome knowledge of the species is gaining attention. In this sense, obtaining a high-density genome mapping of the species could offer clues to the aquaculture improvement in those aspects not resolved so far. In the present article, a review and new processed data have allowed to obtain a high-density BAC-based cytogenetic map of S. senegalensis beside the analysis of the sequences of such BAC clones to achieve integrative data. A total of 93 BAC clones were used to localize the chromosome complement of the species and 588 genes were annotated, thus almost reaching the 2.5% of the S. senegalensis genome sequences. As a result, important data about its genome organization and evolution were obtained, such as the lesser gene density of the large metacentric pair compared with the other metacentric chromosomes, which supports the theory of a sex proto-chromosome pair. In addition, chromosomes with a high number of linked genes that are conserved, even in distant species, were detected. This kind of result widens the knowledge of this species’ chromosome dynamics and evolution.

Development of whole-genome multiplex assays and construction of an integrated genetic map using SSR markers in Senegalese sole

The Senegalese sole (Solea senegalensis) is an economically important flatfish species. In this study, a genome draft was analyzed to identify microsatellite (SSR) markers for whole-genome genotyping. A subset of 224 contigs containing SSRs were preselected and validated by using a de novo female hybrid assembly. Overall, the SSR density in the genome was 886.7 markers per megabase of genomic sequences and the dinucleotide motif was the most abundant (52.4%). In silico comparison identified a set of 108 SSRs (with di-, tetra- or pentanucleotide motifs) widely distributed in the genome and suitable for primer design. A total of 106 markers were structured in thirteen multiplex PCR assays (with up to 10-plex) and the amplification conditions were optimized with a high-quality score. Main genetic diversity statistics and genotyping reliability were assessed. A subset of 40 high polymorphic markers were selected to optimize four supermultiplex PCRs (with up to 11-plex) for pedigree analysis. Theoretical exclusion probabilities and real parentage allocation tests using parent–offspring information confirmed their robustness and effectiveness for parental assignment. These new SSR markers were combined with previously published SSRs (in total 229 makers) to construct a new and improved integrated genetic map containing 21 linkage groups that matched with the expected number of chromosomes. Synteny analysis with respect to C. semilaevis provided new clues on chromosome evolution in flatfish and the formation of metacentric and submetacentric chromosomes in Senegalese sole.

An improved de novo assembling and polishing of Solea senegalensis transcriptome shed light on retinoic acid signalling in larvae

Senegalese sole is an economically important flatfish species in aquaculture and an attractive model to decipher the molecular mechanisms governing the severe transformations occurring during metamorphosis, where retinoic acid seems to play a key role in tissue remodeling. In this study, a robust sole transcriptome was envisaged by reducing the number of assembled libraries (27 out of 111 available), fine-tuning a new automated and reproducible set of workflows for de novo assembling based on several assemblers, and removing low confidence transcripts after mapping onto a sole female genome draft. From a total of 96 resulting assemblies, two "raw" transcriptomes, one containing only Illumina reads and another with Illumina and GS-FLX reads, were selected to provide SOLSEv5.0, the most informative transcriptome with low redundancy and devoid of most single-exon transcripts. It included both Illumina and GS-FLX reads and consisted of 51,348 transcripts of which 22,684 code for 17,429 different proteins described in databases, where 9527 were predicted as complete proteins. SOLSEv5.0 was used as reference for the study of retinoic acid (RA) signalling in sole larvae using drug treatments (DEAB, a RA synthesis blocker, and TTNPB, a RA-receptor agonist) for 24 and 48 h. Differential expression and functional interpretation were facilitated by an updated version of DEGenes Hunter. Acute exposure of both drugs triggered an intense, specific and transient response at 24 h but with hardly observable differences after 48 h at least in the DEAB treatments. Activation of RA signalling by TTNPB specifically increased the expression of genes in pathways related to RA degradation, retinol storage, carotenoid metabolism, homeostatic response and visual cycle, and also modified the expression of transcripts related to morphogenesis and collagen fibril organisation. In contrast, DEAB mainly decreased genes related to retinal production, impairing phototransduction signalling in the retina. A total of 755 transcripts mainly related to lipid metabolism, lipid transport and lipid homeostasis were altered in response to both treatments, indicating non-specific drug responses associated with intestinal absorption. These results indicate that a new assembling and transcript sieving were both necessary to provide a reliable transcriptome to identify the many aspects of RA action during sole development that are of relevance for sole aquaculture.

Dietary Antioxidant Supplementation Promotes Growth in Senegalese Sole Postlarvae

Somatic growth is a balance between protein synthesis and degradation, and it is largely influenced by nutritional clues. Antioxidants levels play a key role in protein turnover by reducing the oxidative damage in the skeletal muscle, and hence promoting growth performance in the long-term. In the present study, Senegalese sole postlarvae (45 days after hatching, DAH) were fed with three experimental diets, a control (CTRL) and two supplemented with natural antioxidants: curcumin (CC) and grape seed (GS). Trial spanned for 25 days and growth performance, muscle cellularity and the expression of muscle growth related genes were assessed at the end of the experiment (70 DAH). The diets CC and GS significantly improved growth performance of fish compared to the CTRL diet. This enhanced growth was associated with larger muscle cross sectional area, with fish fed CC being significantly different from those fed the CTRL. Sole fed the CC diet had the highest number of muscle fibers, indicating that this diet promoted muscle hyperplastic growth. Although the mean fiber diameter did not differ significantly amongst treatments, the proportion of large-sized fibers (>25 μm) was also higher in fish fed the CC diet suggesting increased hypertrophic growth. Such differences in the phenotype were associated with a significant up-regulation of the myogenic differentiation 2 (myod2) and the myomaker (mymk) transcripts involved in myocyte differentiation and fusion, respectively, during larval development. The inclusion of grape seed extract (GS diet) resulted in a significant increase in the expression of myostatin1. These results demonstrate that both diets (CC and GS) can positively modulate muscle development and promote growth in sole postlarvae. This effect is more prominent in CC fed fish, where increased hyperplastic and hypertrophic growth of the muscle was associated with an upregulation of myod2 and mymk genes.

Effects of the sulfated polysaccharide ulvan from Ulva ohnoi on the modulation of the immune response in Senegalese sole (Solea senegalensis)

Sulfated polysaccharides derived from green seaweeds exhibit many beneficial biological activities and have great potential to be used as nutraceutical in aquaculture. In this work, we evaluated the effects of the sulfated polysaccharide ulvan from Ulva ohnoi on Senegalese sole (Solea senegalensis) juveniles at the transcriptomic level. Cytotoxicity assay performed in liver primary cell cultures from sole determined that the different ulvan concentrations assayed did not impair cell viability. Juveniles were intraperitoneally (IP) injected with ulvan (0.5 mg/fish) followed by a challenge with Photobacterium damselae subsp. piscicida (Phdp) at 7 days. RNASeq analyses at 2 days post injection (dpi) revealed that 402 transcripts were differentially expressed in liver between ulvan IP injected and control groups before the challenge. Genes related to bacterial and antiviral defence, complement system, chemokines, proteasomes and antigen presentation were upregulated in ulvan treated groups. A detailed expression analysis of sixteen genes related to innate and adaptive immune system was performed in two systemic tissues: liver and spleen. Ulvan injection provoked the upregulation of tlr22 and a transient inflammatory response was initiated in both liver and spleen at 2 dpi. As consequence, expression of acute phase proteins, antimicrobial peptides and complement genes was induced. Moreover, expression of mhcI, mhcII, psmb10 and bcl6 was also induced 2 dpi. At 2 dpi with Phdp, inflammatory cytokines and genes related to bacterial and antiviral defense, iron metabolism, complement system and antigen presentation were differentially modulated in survival juveniles previously IP injected with ulvan. Moreover, mortality was retarded in ulvan treated juveniles. These results provide new evidence about the role of ulvan as a bioactive compound with immunomodulatory activity in Senegalese sole as well as its possible use as vaccine adjuvant against Phdp. This is the first published study that evaluates the transcriptomic response of Senegalese sole IP injected with ulvan.

Nanostructured recombinant protein particles raise specific antibodies against the nodavirus NNV coat protein in sole

Nervous necrosis virus (NNV) reassortant strains RGNNV/SJNNV have emerged as a potent threat to the Mediterranean marine aquaculture industry, causing viral encephalopathy and retinopathy (VER) in Senegalese sole (Solea senegalensis). In this study, a cheap and practical vaccine strategy using bacterial inclusion bodies made of the coat protein of a virulent reassortant strain of this betanodavirus was devised. The nanostructured recombinant protein nanoparticles, VNNV-C^NP, were administered without adjuvant to two groups of juvenile sole, one by intraperitoneal injection and the other by oral intubation. Specific antibodies were raised in vivo against the NNV coat protein via both routes, with a substantial specific antibody expansion in the injected group 30 days post homologous prime boost. Expression levels of five adaptive immune-related genes, cd8a, cd4, igm, igt and arg2, were also quantified in intestine, spleen and head kidney. Results showed cd4 and igm were upregulated in the head kidney of injected fish, indicating activation of an adaptive systemic response, while intubated fish exhibited a mucosal response in the intestine. Neither route showed significant differential expression of cd8a. The specific antibody response elicited in vivo and the lack of any signs of toxicity over the 6-week study period in young fish (n = 100), evidences the potential of the nanoparticle as a vaccine candidate.

Straying from the flatfish retinal plan: Cone photoreceptor patterning in the common sole (Solea solea) and the Senegalese sole (Solea senegalensis)

The retinas of nonmammalian vertebrates have cone photoreceptor mosaics that are often organized as highly patterned lattice-like distributions. In fishes, the two main lattice-like patterns are composed of double cones and single cones that are either assembled as interdigitized squares or as alternating rows. The functional significance of such orderly patterning is unknown. Here, the cone mosaics in two species of Soleidae flatfishes, the common sole and the Senegalese sole, were characterized and compared to those from other fishes to explore variability in cone patterning and how it may relate to visual function. The cone mosaics of the common sole and the Senegalese sole consisted of single, double, and triple cones in formations that differed from the traditional square mosaic pattern reported for other flatfishes in that no evidence of higher order periodicity was present. Furthermore, mean regularity indices for single and double cones were conspicuously lower than those of other fishes with "typical" square and row mosaics, but comparable to those of goldfish, a species with lattice-like periodicity in its cone mosaic. Opsin transcripts detected by quantitative polymerase chain reaction (sws1, sws2, rh2.3, rh2.4, lws, and rh1) were uniformly expressed across the retina of the common sole but, in the Senegalese sole, sws2, rh2.4, and rh1 were more prevalent in the dorsal retina. Microspectrophotometry revealed five visual pigments in the retina of the common sole [S(472), M(523), M(536), L(559), and rod(511)] corresponding to the repertoire of transcripts quantified except for sws1. Overall, these results indicate a loss of cone mosaic patterning in species that are primarily nocturnal or dwell in low light environments as is the case for the common sole and the Senegalese sole. The corollary is that lattice-like patterning of the cone mosaic may improve visual acuity. Ecological and physiological correlates derived from observations across multiple fish taxa that live in low light environments and do not possess lattice-like cone mosaics are congruent with this claim.

Hox Genes Polymorphism Depicts Developmental Disruption of Common Sole Eggs

In sole aquaculture production, consistency in the quality of produced eggs throughout the year is unpredictable. Hox genes have a crucial role in controlling embryonic development and their genetic variation could alter the phenotype dramatically. In teleosts genome duplication led paralog hox genes to become diverged. Direct association of polymorphism in hoxa1a, hoxa2a & hoxa2b of Solea solea with egg viability indicates hoxa2b as a potential genetic marker. High Resolution Melt (HRM) analysis was carried out in 52 viable and 61 non-viable eggs collected at 54±6 hours post fertilization (hpf). Allelic and genotypic frequencies of polymorphism were analyzed and results illustrated a significantly increased risk for non-viability for minor alleles and their homozygous genotypes. Haplotype analysis demonstrated a significant recessive effect on the risk of non-viability, by increasing the odds of disrupting embryonic development up to three-fold. Phylogenetic analysis showed that the paralog genes hoxa2a and hoxa2b, are separated distinctly in two clades and presented a significant ω variation, revealing their diverged evolutionary rate.

Transcriptional profiling of skeletal muscle reveals starvation response and compensatory growth in Spinibarbus hollandi

BACKGROUND

Spinibarbus hollandi is an economically important fish species in southern China. This fish is known to have nutritional and medicinal properties; however, its farming is limited by its slow growth rate. In the present study, we observed that a compensatory growth phenomenon could be induced by adequate refeeding following 7 days of fasting in S. hollandi. To understand the starvation response and compensatory growth mechanisms in this fish, the muscle transcriptomes of S. hollandi under control, fasting, and refeeding conditions were profiled using next-generation sequencing (NGS) techniques.

RESULTS

More than 4.45 × 108 quality-filtered 150-base-pair Illumina reads were obtained from all nine muscle samples. De novo assemblies yielded a total of 156,735 unigenes, among which 142,918 (91.18%) could be annotated in at least one available database. After 7 days of fasting, 2422 differentially expressed genes were detected, including 1510 up-regulated genes and 912 down-regulated genes. Genes involved in fat, protein, and carbohydrate metabolism were significantly up-regulated, and genes associated with the cell cycle, DNA replication, and immune and cellular structures were inhibited during fasting. After refeeding, 84 up-regulated genes and 16 down-regulated genes were identified. Many genes encoding the components of myofibers were significantly up-regulated. Histological analysis of muscle verified the important role of muscle hypertrophy in compensatory growth.

CONCLUSION

In the present work, we reported the transcriptome profiles of S. hollandi muscle under different conditions. During fasting, the genes involved in the mobilization of stored energy were up-regulated, while the genes associated with growth were down-regulated. After refeeding, muscle hypertrophy contributed to the recovery of growth. The results of this study may help to elucidate the mechanisms underlying the starvation response and compensatory growth.

Genomic and phylogenetic analysis of choriolysins, and biological activity of hatching liquid in the flatfish Senegalese sole

The hatching enzymes or choriolysins are key proteases in fish life cycle controlling the release of larvae to surrounding environment that have been suggested as target for novel biotechnological uses. Due to the large amounts of eggs released by the flatfish Solea senegalensis, during the spawning season, the hatching liquid properties and choriolysin-encoding genes were investigated in this species. A genomic analysis identified four putative genes referred to as SseHCEa, SseHCEb, SseLCE and SseHE. The phylogenetic analysis classified these paralogs into two clades, the clade I containing SseHCE paralogs and the clade II containing two well-supported subclades named as HE and LCE. The two SseHCE paralogs were intron-less and both genes were tandemly arrayed very close in the genome. The synteny and gene rearrangement identified in the flatfish lineage indicated that the duplication of these two paralogs occurred recently and they are under divergent evolution. The genes SseHE and SseLCE were structured in 8 exons and 7 introns and the synteny was conserved in teleosts. Expression studies confirmed that the four genes were expressed in the hatching gland cells and they migrate co-ordinately from the head to around the yolk sac close to the hatch with specific temporal and intensity expression profiles. Although the mRNA levels of the four genes peaked in the hours previous to larval hatching, the SseHCE and SseLCE paralogs kept a longer expression than SseHE after hatching. These expression patterns were consistent even when larvae were incubated at different temperatures that modified hatching times. The analysis of hatching-liquid using SDS-PAGE and zymography analyses of hatching liquid identified a major band of expected choriolysin size. The optimal pH for protease activity was 8.5 and inhibition assays using EDTA demonstrated that most of the activity in the hatching liquid was due to metalloproteases with Ca²⁺ ions acting as the most effective metal to restore the activity. All these data provide new clues about the choriolysin evolution and function in flatfish with impact in the aquaculture and the blue cosmetic industry.

Integrative genetic map of repetitive DNA in the sole Solea senegalensis genome shows a Rex transposon located in a proto-sex chromosome

Repetitive sequences play an essential role in the structural and functional evolution of the genome, particularly in the sexual chromosomes. The Senegalese sole (Solea senegalensis) is a valuable flatfish in aquaculture albeit few studies have addressed the mapping and characterization of repetitive DNA families. Here we analyzed the Simple Sequence Repeats (SSRs) and Transposable elements (TEs) content from fifty-seven BAC clones (spanning 7.9 Mb) of this species, located in chromosomes by multiple fluorescence in situ hybridization (m-BAC-FISH) technique. The SSR analysis revealed an average density of 675.1 loci per Mb and a high abundance (59.69%) of dinucleotide coverage was observed, being ‘AC’ the most abundant. An SSR-FISH analysis using eleven probes was also carried out and seven of the 11 probes yielded positive signals. ‘AC’ probes were present as large clusters in almost all chromosomes, supporting the bioinformatic analysis. Regarding TEs, DNA transposons (Class II) were the most abundant. In Class I, LINE elements were the most abundant and the hAT family was the most represented in Class II. Rex/Babar subfamily, observed in two BAC clones mapping to chromosome pair 1, showed the longest match. This chromosome pair has been recently reported as a putative sexual proto-chromosome in this species, highlighting the possible role of the Rex element in the evolution of this chromosome. In the Rex1 phylogenetic tree, the Senegalese sole Rex1 retrotransposon could be associated with one of the four major ancient lineages in fish genomes, in which it is included O. latipes.

Understanding pseudo-albinism in sole (Solea senegalensis): a transcriptomics and metagenomics approach

Pseudo-albinism is a pigmentation disorder observed in flatfish aquaculture with a complex, multi-factor aetiology. We tested the hypothesis that pigmentation abnormalities are an overt signal of more generalised modifications in tissue structure and function, using as a model the Senegalese sole and two important innate immune barriers, the skin and intestine, and their microbiomes. Stereological analyses in pseudo-albino sole revealed a significantly increased mucous cell number in skin (P < 0.001) and a significantly thicker muscle layer and lamina propria in gut (P < 0.001). RNA-seq transcriptome analysis of the skin and gut identified 573 differentially expressed transcripts (DETs, FDR < 0.05) between pseudo-albino and pigmented soles (one pool/tissue from 4 individuals/phenotype). DETs were mainly linked to pigment production, skin structure and regeneration and smooth muscle contraction. The microbiome (16 S rRNA analysis) was highly diverse in pigmented and pseudo-albino skin but in gut had low complexity and diverged between the two pigmentation phenotypes. Quantitative PCR revealed significantly lower loads of Mycoplasma (P < 0.05) and Vibrio bacteria (P < 0.01) in pseudo-albino compared to the control. The study revealed that pseudo-albinism in addition to pigmentation changes was associated with generalised changes in the skin and gut structure and a modification in the gut microbiome.

Nitric oxide-dependent regulation of sweet pepper fruit ripening

Ripening is a complex physiological process that involves changes in reactive nitrogen and oxygen species that govern the shelf-life and quality of fruits. Nitric oxide (NO)-dependent changes in the sweet pepper fruit transcriptome were determined by treating fruits at the initial breaking point stage with NO gas. Fruits were also harvested at the immature (green) and ripe (red) stages. Fruit ripening in the absence of NO resulted in changes in the abundance of 8805 transcripts whose function could be identified. Among these, functional clusters associated with reactive oxygen/nitrogen species and lipid metabolism were significantly modified. NO treatment resulted in the differential expression of 498 genes framed within these functional categories. Biochemical analysis revealed that NO treatment resulted in changes in fatty acid profiling, glutathione and proline contents, and the extent of lipid peroxidation, as well as increases in the activity of ascorbate peroxidase and lipoxygenase. These data provide supporting evidence for the crucial role of NO in the ripening of pepper fruit.

Feed and immersion challenges with lymphocystis disease virus (LCDV) reveals specific mechanisms for horizontal transmission and immune response in senegalese sole post-larvae

The horizontal transmission of lymphocystis disease virus (LCDV) through contaminated water and feed (using artemia as vehicle) and the associated immune gene expression profiles in Senegalese sole post-larvae were investigated. All specimens analyzed were positive for LCDV DNA detection at 1-day post-challenge (1 dpc) with the highest viral levels in specimens infected through the immersion route. However, the percentage of LCDV-positive animals and number of viral DNA copies dropped progressively at 2 and 7 dpc. The histological analysis identified structural changes in the skin, muscle and gills of sole post-larvae LCDV-challenged by immersion. In situ hybridization confirmed a wide distribution of LCDV in the skin, gut, surrounding vessels in trunk muscle and head kidney in the immersion route, while the signals were restricted to the liver and lamina propria in the feeding treatment. Expression analysis using a set of 22 genes related to innate immune defense system demonstrated clear differences in the time-course response to LCDV as function of the infection route. Most antiviral defense genes, the proinflammatory cytokines, the complement c3, g-type lysozyme and T-cell markers cd4 and cd8a were rapidly induced in the feeding-infected post-larvae, and they were remained activated at 2 dpc. In contrast, in the immersion-infected post-larvae the induction of most defensive genes was delayed, with a low intensity at 2 dpc. All these data demonstrate that LCDV can horizontally infect Senegalese sole post-larvae through the water or feed although with different patterns of histopathological disorders, virus distribution and route-specific expression profiles.

Sole head transcriptomics reveals a coordinated developmental program during metamorphosis

Most teleosts undergo a thyroid hormone (TH) regulated larval to juvenile transition known as metamorphosis. In Pleuronectiformes (flatfish), metamorphosis is most dramatic, and one eye of the symmetric pelagic larvae migrates to the opposite side of the head, giving rise to an asymmetric benthic juvenile with both eyes on the same side of the body. Asymmetric development occurs mostly in the head. To understand the genetic mechanisms underlying this developmental change we have generated a Solea senegalensis metamorphosing flatfish head transcriptome. Our results reveal that THs acting as integrative factors direct a stepwise genetic program that initiates a specific organismal level response followed by cell specific responses that lead to the long-term changes that characterise the post-metamorphic identity and physiology of the head. Flatfish head asymmetric development during metamorphosis and its TH dependency is conserved thus we consider the findings in sole most likely representative of all flatfish species.

TransFlow: a modular framework for assembling and assessing accurate de novo transcriptomes in non-model organisms

BACKGROUND

The advances in high-throughput sequencing technologies are allowing more and more de novo assembling of transcriptomes from many new organisms. Some degree of automation and evaluation is required to warrant reproducibility, repetitivity and the selection of the best possible transcriptome. Workflows and pipelines are becoming an absolute requirement for such a purpose, but the issue of assembling evaluation for de novo transcriptomes in organisms lacking a sequenced genome remains unsolved. An automated, reproducible and flexible framework called TransFlow to accomplish this task is described.

RESULTS

TransFlow with its five independent modules was designed to build different workflows depending on the nature of the original reads. This architecture enables different combinations of Illumina and Roche/454 sequencing data, and can be extended to other sequencing platforms. Its capabilities are illustrated with the selection of reliable plant reference transcriptomes and the assembling six transcriptomes (three case studies for grapevine leaves, olive tree pollen, and chestnut stem, and other three for haustorium, epiphytic structures and their combination for the phytopathogenic fungus Podosphaera xanthii). Arabidopsis and poplar transcriptomes revealed to be the best references. A common result regarding de novo assemblies is that Illumina paired-end reads of 100 nt in length assembled with OASES can provide reliable transcriptomes, while the contribution of longer reads is noticeable only when they complement a set of short, single-reads.

CONCLUSIONS

TransFlow can handle up to 181 different assembling strategies. Evaluation based on principal component analyses allows its self-adaptation to different sets of reads to provide a suitable transcriptome for each combination of reads and assemblers. As a result, each case study has its own behaviour, prioritises evaluation parameters, and gives an objective and automated way for detecting the best transcriptome within a pool of them. Sequencing data type and quantity (preferably several hundred millions of 2×100 nt or longer), assemblers (OASES for Illumina, MIRA4 and EULER-SR reconciled with CAP3 for Roche/454) and strategy (preferably scaffolding with OASES, and probably merging with Roche/454 when available) arise as the most impacting factors.

Evidence for a Robertsonian fusion in Solea senegalensis (Kaup, 1858) revealed by zoo-FISH and comparative genome analysis

Background

Solea senegalensis (Kaup, 1858) is a commercially important flatfish species, belonging to the Pleuronectiformes order. The taxonomy of this group has long been controversial, and the karyotype of the order presents a high degree of variability in diploid number, derived from chromosomal rearrangements such as Robertsonian fusions. Previously it has been proposed that the large metacentric chromosome of S. senegalensis arises from this kind of chromosome rearrangement and that this is a proto-sex chromosome.

Results

In this work, the Robertsonian origin of the large metacentric chromosome of S. senegalensis has been tested by the Zoo-FISH technique applied to two species of the Soleidae family (Dicologlossa cuneata and Dagetichthys lusitanica), and by comparative genome analysis with Cynoglossus semilaevis. From the karyotypic analysis we were able to determine a chromosome complement comprising 2n = 50 (FN = 54) in D. cuneata and 2n = 42 (FN = 50) in D. lusitanica. The large metacentric painting probe gave consistent signals in four acrocentric chromosomes of the two Soleidae species; and the genome analysis proved a common origin with four chromosome pairs of C. semilaevis. As a result of the genomic analysis, up to 61 genes were annotated within the thirteen Bacterial Artificial Chromosome clones analysed.

Conclusions

These results confirm that the large metacentric chromosome of S. senegalensis originated from a Robertsonian fusion and provide new data about the chromosome evolution of S. senegalensis in particular, and of Pleuronectiformes in general.

Electronic supplementary material

The online version of this article (10.1186/s12864-018-5216-6) contains supplementary material, which is available to authorized users.

Antioxidant capacity and immunomodulatory effects of a chrysolaminarin-enriched extract in Senegalese sole

The microalgae are an important source of bioactive molecules including β-glucans that can be used as immunostimulants in aquaculture. In the present study, the antioxidant capacity, cytotoxicity and immunomodulatory activity of a chrysolaminarin-enriched extract obtained from the diatom Phaeodactylum tricornutum was evaluated. The extract showed a higher total antioxidant activity as determined by ORAC and FRAP assays and a lower DPPH scavenging activity than particulate yeast-β-glucan. The cytotoxicity test indicated that extract concentrations higher than 0.01% w/v could impair cell viability of human dermal fibroblasts. To evaluate the immunomodulatory activity, juvenile soles were intraperitoneally injected with the chrysolaminarin-enriched extract suspended in coconut oil (1 mg/fish) followed by a reinjection at 7 days. A sham group injected with the carrier solution was maintained as a negative control. Cumulated mortality of fish injected with the chrysolaminarin-enriched extract was 29.4% after six days and no mortality was recorded after extract reinjection. Expression analyses of fifteen genes related to the innate immune system in kidney, spleen and intestine showed temporal and organ-specific responses. A rapid (2 days post-injection; dpi) and strong induction of the pro-inflammatory il1b and the antimicrobial peptide hamp1 in the three immunological organs, the hsp90aa in kidney and spleen, irf3 in intestine and c3 in spleen was observed indicating a potent inflammatory response. The recovery of steady-state levels for all activated genes at 5 dpi, and the down-regulation of c-lectin receptor as well as some interferon-related genes (ifn1, irf1, irf3, irf8, irf9 and mx) in kidney and cxc10 in spleen indicated that the soles were able to activate a homeostatic response against the β-glucan insult. The reinjection of the chrysolaminarin-enriched extract did not activate a new inflammatory response but reduced the mRNA levels of hsp90aa and irf3 indicating that soles developed some resistance to β-glucans. Overall, these results reveal this enriched extract as a novel and potent source of β-glucans with antioxidant and immunomodulatory capacity suitable for immunostimulation in aquaculture.

Transcriptomic Profiles of Senegalese Sole Infected With Nervous Necrosis Virus Reassortants Presenting Different Degree of Virulence

Betanodaviruses [nervous necrosis virus (NNV)] are the causative agent of the viral encephalopathy and retinopathy, a disease that affects cultured Senegalese sole (Solea senegalensis). NNV reassortants, combining genomic segments from redspotted grouper nervous necrosis virus (RGNNV) and striped jack nervous necrosis virus (SJNNV) genotypes, have been previously isolated from several fish species. The wild-type reassortant wSs160.03, isolated from Senegalese sole, has been proven to be more virulent to sole than the parental genotypes (RGNNV and SJNNV), causing 100% mortality. Mutations at amino acids 247 (serine to alanine) and 270 (serine to asparagine) in the wSs160.03 capsid protein have allowed us to obtain a mutant reassortant (rSs160.03_247+270), which provokes a 40% mortality decrease. In this study, the RNA-Seq technology has been used to comparatively analyze Senegalese sole transcriptomes in two organs (head kidney and eye/brain) after infection with wild-type and mutant strains. A total of 633 genes were differentially expressed (DEGs) in animals infected with the wild-type isolate (with higher virulence), whereas 393 genes were differentially expressed in animals infected with the mutant strain (37.9% decrease in the number of DEGs). To study the biological functions of detected DEGs involved in NNV infection, a gene ontology (GO) enrichment analysis was performed. Different GO profiles were obtained in the following subclasses: (i) biological process; (ii) cellular component; and (iii) molecular function, for each viral strain tested. Immune response and proteolysis have been the predominant biological process after the infection with the wild-type isolate, whereas the infection with the mutant strain induces proteolysis in head kidney and inhibition of vasculogenesis in nervous tissue. Regarding the immune response, genes coding for proteins acting as mediators of type I IFN expression (DHX58, IRF3, IRF7) and IFN-stimulated genes (ISG15, Mx, PKR, Gig1, ISG12, IFI44, IFIT-1, to name a few) were upregulated in animals infected with the wild-type isolate, whereas no-differential expression of these genes was observed in samples inoculated with the mutant strain. The different transcriptomic profiles obtained could help to better understand the NNV pathogenesis in Senegalese sole, setting up the importance as virulence determinants of amino acids at positions 247 and 270 within the RNA2 segment.

Sma3s: A universal tool for easy functional annotation of proteomes and transcriptomes

The current cheapening of next-generation sequencing has led to an enormous growth in the number of sequenced genomes and transcriptomes, allowing wet labs to get the sequences from their organisms of study. To make the most of these data, one of the first things that should be done is the functional annotation of the protein-coding genes. But it used to be a slow and tedious step that can involve the characterization of thousands of sequences. Sma3s is an accurate computational tool for annotating proteins in an unattended way. Now, we have developed a completely new version, which includes functionalities that will be of utility for fundamental and applied science. Currently, the results provide functional categories such as biological processes, which become useful for both characterizing particular sequence datasets and comparing results from different projects. But one of the most important implemented innovations is that it has now low computational requirements, and the complete annotation of a simple proteome or transcriptome usually takes around 24 hours in a personal computer. Sma3s has been tested with a large amount of complete proteomes and transcriptomes, and it has demonstrated its potential in health science and other specific projects.

A transcriptomic approach to study the effect of long-term starvation and diet composition on the expression of mitochondrial oxidative phosphorylation genes in gilthead sea bream (Sparus aurata)

Background

The impact of nutritional status and diet composition on mitochondrial oxidative phosphorylation (OXPHOS) in fish remains largely unknown. To identify biomarkers of interest in nutritional studies, herein we obtained a deep-coverage transcriptome by 454 pyrosequencing of liver and skeletal muscle cDNA normalised libraries from long-term starved gilthead sea bream (Sparus aurata) and fish fed different diets.

Results

After clean-up of high-throughput deep sequencing reads, 699,991 and 555,031 high-quality reads allowed de novo assembly of liver and skeletal muscle sequences, respectively (average length: 374 and 441 bp; total megabases: 262 and 245 Mbp). An additional incremental assembly was completed by integrating data from both tissues (hybrid assembly). Assembly of hybrid, liver and skeletal muscle transcriptomes yielded, respectively, 19,530, 11,545 and 10,599 isotigs (average length: 1330, 1208 and 1390 bp, respectively) that were grouped into 15,954, 10,033 and 9189 isogroups. Following annotation, hybrid transcriptomic data were used to construct an oligonucleotide microarray to analyse nutritional regulation of the expression of 129 genes involved in OXPHOS in S. aurata. Starvation upregulated cytochrome c oxidase components and other key OXPHOS genes in the liver, which exhibited higher sensitive to food deprivation than the skeletal muscle. However, diet composition affected OXPHOS in the skeletal muscle to a greater extent than in the liver: most of genes upregulated under starvation presented higher expression among fish fed a high carbohydrate/low protein diet.

Conclusions

Our findings indicate that the expression of coenzyme Q-binding protein (COQ10), cytochrome c oxidase subunit 6A2 (COX6A2) and ADP/ATP translocase 3 (SLC25A6) in the liver, and cytochrome c oxidase subunit 5B isoform 1 (COX5B1) in the liver and the skeletal muscle, are sensitive markers of the nutritional condition that may be relevant to assess the effect of changes in the feeding regime and diet composition on fish farming.

Electronic supplementary material

The online version of this article (10.1186/s12864-017-4148-x) contains supplementary material, which is available to authorized users.

Phylogeny, expression patterns and regulation of DNA Methyltransferases in early development of the flatfish, Solea senegalensis

Background

The identification of DNA methyltransferases (Dnmt) expression patterns during development and their regulation is important to understand the epigenetic mechanisms that modulate larval plasticity in marine fish. In this study, dnmt1 and dnmt3 paralogs were identified in the flatfish Solea senegalensis and expression patterns in early developmental stages and juveniles were determined. Additionally, the regulation of Dnmt transcription by a specific inhibitor (5-aza-2′-deoxycytidine) and temperature was evaluated.

Results

Five paralog genes of dnmt3, namely dnmt3aa, dnmt3ab, dnmt3ba, dnmt3bb.1 and dnmt3bb.2 and one gene for dnmt1 were identified. Phylogenetic analysis revealed that the dnmt gene family was highly conserved in teleosts and three fish-specific genes, dnmt3aa, dnmt3ba and dnmt3bb.2 have evolved. The spatio-temporal expression patterns of four dnmts (dnmt1, dnmt3aa, dnmt3ab and dnmt3bb.1) were different in early larval stages although all of them reduced expression with the age and were detected in neural organs and dnmt3aa appeared specific to somites. In juveniles, the four dnmt genes were expressed in brain and hematopoietic tissues such as kidney, spleen and gills. Treatment of sole embryos with 5-aza-2′-deoxycytidine down-regulated dntm1 and up-regulated dntm3aa. Moreover, in lecithotrophic larval stages, dnmt3aa and dnmt3ab were temperature sensitive and their expression was higher in larvae incubated at 16 °C relative to 20 °C.

Conclusion

Five dnmt3 and one dnmt1 paralog were identified in sole and their distinct developmental and tissue-specific expression patterns indicate that they may have different roles during development. The inhibitor 5-aza-2′-deoxycytidine modified the transcript abundance of dntm1 and dntm3aa in embryos, which suggests that a regulatory feedback mechanism exists for these genes. The impact of thermal regime on expression levels of dnmt3aa and dnmt3ab in lecithotrophic larval stages suggests that these paralogs might be involved in thermal programing.

Electronic supplementary material

The online version of this article (doi:10.1186/s12861-017-0154-0) contains supplementary material, which is available to authorized users.

Molecular characterization and expression analyses of the Solea senegalensis interferon-stimulated gene 15 (isg15) following NNV infections

Interferons are essential in fish resistance to viral infections. They induce interferon-stimulated genes, such as isg15. In this study, the Senegalese sole isg15 gene (ssisg15) has been characterized. As other isg15, ssisg15 contains a 402-bp intron sited in the 5'-UTR, and the full length cDNA is 1492-bp, including a 480-bp ORF. The expression analyses revealed basal levels of isg15 transcripts, and a clear induction after poly I:C injection, that reached maximum values in brain, head kidney and gills. The ssisg15 induction patterns were similar in RGNNV- and SJNNV-inoculated fish, whereas the reassortant (RG/SJ) isolate, which has higher replication fitness, triggered delayed but higher transcript levels. Furthermore, RG/SJ infection after poly I:C treatment reduced the induction of ssisg15 transcripts, suggesting an antagonistic mechanism against interferon type I system, that might allow an efficient viral replication at the initial steps of the infective process.

Gene expression profiles associated with lymphocystis disease virus (LCDV) in experimentally infected Senegalese sole (Solea senegalensis)

In the present study, the pathogenesis of lymphocystis disease virus (LCDV) and the immune gene expression patterns associated with this viral infection were determined in the flatfish Senegalese sole. The results indicate that LCDV spreads rapidly from the peritoneal cavity through the bloodstream to reach target organs such as kidney, gut, liver, and skin/fin. The viral load was highest in kidney and reduced progressively thorough the experiment in spite of the viral major capsid protein gene was transcribed. The LCDV injection activated a similar set of differentially expressed transcripts in kidney and intestine although with some differences in the intensity and time-course response. This set included antiviral-related transcripts (including the mx and interferon-related factors irf1, irf2, irf3, irf7, irf8, irf9, irf10), cytokines (il1b, il6, il8, il12 and tnfa) and their receptors (il1r, il8r, il10r, il15ra, il17r), chemokines (CXC-type, CC-type and IL-8), prostaglandins (cox-2), g-type lysozymes, hepcidin, complement fractions (c2, c4-1 and c4-2) and the antigen differentiation factors cd4, cd8a, and cd8b. The expression profile observed indicated that the host triggered a systemic defensive response including inflammation able to cope with the viral challenge.

Integrating genomic resources of flatfish (Pleuronectiformes) to boost aquaculture production

Flatfish have a high market acceptance thus representing a profitable aquaculture production. The main farmed species is the turbot (Scophthalmus maximus) followed by Japanese flounder (Paralichthys olivaceous) and tongue sole (Cynoglossus semilaevis), but other species like Atlantic halibut (Hippoglossus hippoglossus), Senegalese sole (Solea senegalensis) and common sole (Solea solea) also register an important production and are very promising for farming. Important genomic resources are available for most of these species including whole genome sequencing projects, genetic maps and transcriptomes. In this work, we integrate all available genomic information of these species within a common framework, taking as reference the whole assembled genomes of turbot and tongue sole (>210× coverage). New insights related to the genetic basis of productive traits and new data useful to understand the evolutionary origin and diversification of this group were obtained. Despite a general 1:1 chromosome syntenic relationship between species, the comparison of turbot and tongue sole genomes showed huge intrachromosomic reorganizations. The integration of available mapping information supported specific chromosome fusions along flatfish evolution and facilitated the comparison between species of previously reported genetic associations for productive traits. When comparing transcriptomic resources of the six species, a common set of ~2500 othologues and ~150 common miRNAs were identified, and specific sets of putative missing genes were detected in flatfish transcriptomes, likely reflecting their evolutionary diversification.

Analysis of the histone cluster in Senegalese sole (Solea senegalensis): evidence for a divergent evolution of two canonical histone clusters

The Senegalese sole (Solea senegalensis) is commercially very important and a priority species for aquaculture product diversification. The main histone cluster was identified within two BAC clones. However, two replacement histones (H1.0 and H3.3) were found in another BAC clone. Different types of canonical histones H2A and H2B were found within the same species for the first time. Phylogenetic analysis demonstrated that the different types of H1, H2A, and H2B histones were all more similar to each other than to canonical histones from other species. The canonical histone H3 of S. senegalensis differs from subtypes H3.1 and H3.2 in humans at the site of residue 96, where a serine is found instead of an alanine. This same polymorphism has been found only in Danio rerio. The karyotype of S. senegalensis comprises 21 pairs of chromosomes, distributed in 3 metacentric pairs, 2 submetacentric pairs, 4 subtelocentric pairs, and 12 acrocentric pairs. The two BAC clones that contain the clusters of canonical histones were both mapped on the largest metacentric pair, and mFISH analysis confirmed the co-location with the dmrt1 gene in that pair. Three chromosome markers have been identified which, in addition to those previously described, account for 18 chromosome pairs in S. senegalensis.

Ontogenetic onset of immune-relevant genes in the common sole (Solea solea)

Fish are free-living organisms since initial stages of development and are exposed to numerous pathogens before their lymphoid organs have matured and adaptive immunity has developed. Susceptibility to diseases and juvenile mortality represent key critical factors for aquaculture. In this context, the characterization of the appearance kinetics of the immune system key members will be useful in understanding the ability of a particular species in generating immune protection against invading pathogens at different developmental stages. The present study characterized, for the first time, the transcriptional onset of un-explored relevant genes of both innate and adaptive immune system during the Solea solea ontogenesis. Gene expression profiles of immune relevant genes was investigated, by means of DNA microarray, in ten developmental stages, from hatching (1 day post-hatching, dph) to accomplishment of the juvenile form (33 dph). The obtained results revealed that transcripts encoding relevant members of innate immune repertoire, such as lysozyme, AMPs (hepcidin, β-defensin), PPRs and complement components are generally characterized by high expression levels at first stages (i.e. hatch and first feeding) indicating protection from environmental pathogens even at early development. Transcription of adaptive immune genes (i.e. Class I and class II MHC, TCRs) differs from that of the innate immune system. Their onset coincides with metamorphosis and larvae-to-juvenile transition, and likely overlaps with the appearance and maturation of the main lymphoid organs. Finally, data collected suggest that at the end of metamorphosis S. solea cell-mediated immune system hasn't still undergone full maturation.

Transcriptomic profiles of the upper olfactory rosette in cultured and wild Senegalese sole (Solea senegalensis) males

The aims of this study were the characterization of the upper olfactory epithelium of cultured and wild Senegalese sole mature males at histological and transcriptomic (using RNA-Seq) level. No significant differences in tissue structure, cell types and cellular distribution pattern (olfactory sensory neurons) were identified between cultured and wild specimens. Deep transcriptomic analysis showed 2387 transcripts were differentially expressed between cultured and wild groups. A detailed analysis identified the differentially expressed transcripts included some olfactory receptors (OR, TAAR and V2R-like) and transcripts related with the control of reproduction such as the brain aromatase cytochrome P450 and tachykinin-3. Also a wide set of genes related with lipid sensing, metabolism and transport were differentially expressed and these transcripts were often down-regulated in cultured fish. Furthermore, cultured males presented a higher expression of genes related with goblet cells and mucin production that modulates innate and adaptive immune responses. All these changes in gene expression could be explained by different nutritional status and diet preference. The different expression of transcripts related to olfaction, reproduction, nutrient sensing and immune system demonstrate distinct differences in functionalities between cultured and wild soles providing new clues about the sexual dysfunction in this species.

Transcriptional regulation of genes involved in retinoic acid metabolism in Senegalese sole larvae

The aim of this study was the characterization of transcriptional regulatory pathways mediated by retinoic acid (RA) in Senegalese sole larvae. For this purpose, pre-metamorphic larvae were treated with a low concentration of DEAB, an inhibitor of RALDH enzyme, until the end of metamorphosis. No differences in growth, eye migration or survival were observed. Nevertheless, gene expression analysis revealed a total of 20 transcripts differentially expressed during larval development and only six related with DEAB treatments directly involved in RA metabolism and actions (rdh10a, aldh1a2, crbp1, igf2r, rarg and cyp26a1) to adapt to a low-RA environment. In a second experiment, post-metamorphic larvae were exposed to the all-trans RA (atRA) observing an opposite regulation for those genes involved in RA synthesis and degradation (rdh10a, aldh1a2, crbp1 and cyp26a1) as well as other related with thyroid- (dio2) and IGF-axes (igfbp1, igf2r and igfbp5) to balance RA levels. In a third experiment, DEAB-pretreated post-metamorphic larvae were exposed to atRA and TTNPB (a specific RAR agonist). Both drugs down-regulated rdh10a and aldh1a2 and up-regulated cyp26a1 expression demonstrating their important role in RA homeostasis. Moreover, five retinoic receptors that mediate RA actions, the thyroid receptor thrb, and five IGF binding proteins changed differentially their expression. Overall, this study demonstrates that exogenous RA modulates the expression of some genes involved in the RA synthesis, degradation and cellular transport through RAR-mediated regulatory pathways establishing a negative feedback regulatory mechanism necessary to balance endogenous RA levels and gradients.

Molecular characterization and developmental expression patterns of apolipoprotein A-I in Senegalese sole (Solea senegalensis Kaup)

The apolipoprotein A-I (ApoA-I) is an essential component of the high density lipoproteins (HDL). In this study, the cDNA and genomic sequences of this apolipoprotein were characterized for first time in Solea senegalensis. The predicted polypeptide revealed conserved structural features including ten repeats in the lipid-binding domain and some residues involved in cholesterol interaction and binding. The gene structure analysis identified four exons and three introns. Moreover, the synteny analysis revealed that apoA-I did not localize with other apolipoproteins indicating a divergent evolution with respect to the apoA-IV and apoE cluster. The phylogenetic analyses identified two distinct apoA-I paralogs in Ostariophysi (referred to as Ia and Ib) and only one (Ib) in Acanthopterygii. Whole-mount in situ hybridization located the apoA-I signal mainly in the yolk syncytial layer in lecitotrophic larval stages. Later at mouth opening, the mRNA signals were detected mainly in liver and intestine compatible with its role in the HDL formation. Moreover, a clear signal was detected in some regions of the brain, retina and neural cord suggesting a role in local regulation of cholesterol homeostasis. After metamorphosis, apoA-I was also detected in other tissues such as gills, head kidney and spleen suggesting a putative role in immunity. Expression analyses in larvae fed two diets with different triacylglycerol levels indicated that apoA-I mRNA levels were more associated to larval size and development than dietary lipid levels. Finally, qPCR analyses of immature and mature transcripts revealed distinct expression profiles suggesting a posttranscriptional regulatory mechanism.

The transcriptome of metamorphosing flatfish

BACKGROUND

Flatfish metamorphosis denotes the extraordinary transformation of a symmetric pelagic larva into an asymmetric benthic juvenile. Metamorphosis in vertebrates is driven by thyroid hormones (THs), but how they orchestrate the cellular, morphological and functional modifications associated with maturation to juvenile/adult states in flatfish is an enigma. Since THs act via thyroid receptors that are ligand activated transcription factors, we hypothesized that the maturation of tissues during metamorphosis should be preceded by significant modifications in the transcriptome. Targeting the unique metamorphosis of flatfish and taking advantage of the large size of Atlantic halibut (Hippoglossus hippoglossus) larvae, we determined the molecular basis of TH action using RNA sequencing.

RESULTS

De novo assembly of sequences for larval head, skin and gastrointestinal tract (GI-tract) yielded 90,676, 65,530 and 38,426 contigs, respectively. More than 57 % of the assembled sequences were successfully annotated using a multi-step Blast approach. A unique set of biological processes and candidate genes were identified specifically associated with changes in morphology and function of the head, skin and GI-tract. Transcriptome dynamics during metamorphosis were mapped with SOLiD sequencing of whole larvae and revealed greater than 8,000 differentially expressed (DE) genes significantly (p < 0.05) up- or down-regulated in comparison with the juvenile stage. Candidate transcripts quantified by SOLiD and qPCR analysis were significantly (r = 0.843; p < 0.05) correlated. The majority (98 %) of DE genes during metamorphosis were not TH-responsive. TH-responsive transcripts clustered into 6 groups based on their expression pattern during metamorphosis and the majority of the 145 DE TH-responsive genes were down-regulated.

CONCLUSIONS

A transcriptome resource has been generated for metamorphosing Atlantic halibut and over 8,000 DE transcripts per stage were identified. Unique sets of biological processes and candidate genes were associated with changes in the head, skin and GI-tract during metamorphosis. A small proportion of DE transcripts were TH-responsive, suggesting that they trigger gene networks, signalling cascades and transcription factors, leading to the overt changes in tissue occurring during metamorphosis.

Integrated gene mapping and synteny studies give insights into the evolution of a sex proto-chromosome in Solea senegalensis

The evolution of genes related to sex and reproduction in fish shows high plasticity and, to date, the sex determination system has only been identified in a few species. Solea senegalensis has 42 chromosomes and an XX/XY chromosome system for sex determination, while related species show the ZZ/ZW system. Next-generation sequencing (NGS), multi-color fluorescence in situ hybridization (mFISH) techniques, and bioinformatics analysis have been carried out, with the objective of revealing new information about sex determination and reproduction in S. senegalensis. To that end, several bacterial artificial chromosome (BAC) clones that contain candidate genes involved in such processes (dmrt1, dmrt2, dmrt3, dmrt4, sox3, sox6, sox8, sox9, lh, cyp19a1a, amh, vasa, aqp3, and nanos3) were analyzed and compared with the same region in other related species. Synteny studies showed that the co-localization of dmrt1-dmrt2-drmt3 in the largest metacentric chromosome of S. senegalensis is coincident with that found in the Z chromosome of Cynoglossus semilaevis, which would potentially make this a sex proto-chromosome. Phylogenetic studies show the close proximity of S. senegalensis to Oryzias latipes, a species with an XX/XY system and a sex master gene. Comparative mapping provides evidence of the preferential association of these candidate genes in particular chromosome pairs. By using the NGS and mFISH techniques, it has been possible to obtain an integrated genetic map, which shows that 15 out of 21 chromosome pairs of S. senegalensis have at least one BAC clone. This result is important for distinguishing those chromosome pairs of S. senegalensis that are similar in shape and size. The mFISH analysis shows the following co-localizations in the same chromosomes: dmrt1-dmrt2-dmrt3, dmrt4-sox9-thrb, aqp3-sox8, cyp19a1a-fshb, igsf9b-sox3, and lysg-sox6.

Dietary fatty acid composition affects food intake and gut-brain satiety signaling in Senegalese sole (Solea senegalensis, Kaup 1858) larvae and post-larvae

Little is known how dietary lipids affect food intake during larval development of fish, especially with regard to fatty acid (FA) composition. In fact, very little work has been done on appetite regulation and food intake in fish larvae in general, due to biological and technical difficulties associated with this type of studies. A new method using fluorescent microspheres as markers was developed in this study to evaluate food intake and prey selectivity of Senegalese sole larvae and post-larvae. Food intake was quantified in fish fed Artemia metanauplii enriched with oils differing in FA profile: cod liver oil (CLO), linseed oil (LSO), soybean oil (SBO) or olive oil (OO). The fish did not preferentially ingest a specific diet when presented with a choice. However, pre-metamorphic larvae from the CLO treatment ingested more metanauplii per g body weight, while differences in post-larvae were not significant. These findings were developed further by analyzing mRNA levels of a range of putative anorexigenic (pyya, pyyb, glp1, cckl, cart1a, cart1b, cart2a, cart4, pomca, pomcb, crf) and orexigenic (gal, npy, agrp2) genes, to identify those which are significantly affected by feeding and/or dietary FA composition. The variety of expression patterns observed highlighted the complexity of appetite regulatory mechanisms. In general, fish fed the CLO diet tended to show gene expression patterns most dissimilar to the remaining treatments. Expression in pre-metamorphic larvae was generally less in accordance with the putative function of the genes than in post-larvae, which could suggest a yet underdeveloped regulatory system.

Role of the IFN I system against the VHSV infection in juvenile Senegalese sole (Solea senegalensis)

Senegalese sole is susceptible to marine VHSV isolates but is not affected by freshwater isolates, which may indicate differences regarding virus-host immune system interaction. IFN I induces an antiviral state in fish, stimulating the expression of genes encoding antiviral proteins (ISG). In this study, the stimulation of the Senegalese sole IFN I by VHSV infections has been evaluated by the relative quantification of the transcription of several ISG (Mx, Isg15 and Pkr) after inoculation with marine (pathogenic) and freshwater (non-pathogenic) VHSV isolates. Compared to marine VHSV, lower levels of RNA of the freshwater VHSV induced transcription of ISG to similar levels, with the Isg15 showing the highest fold induction. The protective role of the IFN I system was evaluated in poly I:C-inoculated animals subsequently challenged with VHSV isolates. The cumulative mortality caused by the marine isolate in the control group was 68%, whereas in the poly I:C-stimulated group was 5%. The freshwater VHSV isolate did not cause any mortality. Furthermore, viral RNA fold change and viral titers were lower in animals from the poly I:C + VHSV groups than in the controls. The implication of the IFN I system in the protection observed was confirmed by the transcription of the ISG in animals from the poly I:C + VHSV groups. However, the marine VHSV isolate exerts a negative effect on the ISG transcription at 3 and 6 h post-inoculation (hpi), which is not observed for the freshwater isolate. This difference might be partly responsible for the virulence shown by the marine isolate.

Hypothalamic fatty acid sensing in Senegalese sole (Solea senegalensis): response to long-chain saturated, monounsaturated, and polyunsaturated (n-3) fatty acids

We assessed the presence of fatty acid (FA)-sensing mechanisms in hypothalamus of Senegalese sole (Solea senegalensis) and investigated their sensitivity to FA chain length and/or level of unsaturation. Stearate (SA, saturated FA), oleate (OA, monounsaturated FA of the same chain length), α-linolenate [ALA, a n-3 polyunsaturated fatty acid (PUFA) of the same chain length], and eicosapentanoate (EPA, a n-3 PUFA of a larger chain length) were injected intraperitoneally. Parameters related to FA sensing and neuropeptide expression in the hypothalamus were assessed after 3 h and changes in accumulated food intake after 4, 24, and 48 h. Three FA sensing systems characterized in rainbow trout were also found in Senegalese sole and were activated by OA in a way similar to that previously characterized in rainbow trout and mammals. These hypothalamic FA sensing systems were also activated by ALA, differing from mammals, where n-3 PUFAs do not seem to activate FA sensors. This might suggest additional roles and highlights the importance of n-3 PUFA in fish diets, especially in marine species. The activation of FA sensing seems to be partially dependent on acyl chain length and degree of saturation, as no major changes were observed after treating fish with SA or EPA. The activation of FA sensing systems by OA and ALA, but not SA or EPA, is further reflected in the expression of hypothalamic neuropeptides involved in the control of food intake. Both OA and ALA enhanced anorexigenic capacity compatible with the activation of FA sensing systems.

Genomic characterization and expression analysis of four apolipoprotein A-IV paralogs in Senegalese sole (Solea senegalensis Kaup)

The apolipoprotein A-IV (ApoA-IV) plays a key role in lipid transport and feed intake regulation. In this work, four cDNA sequences encoding ApoA-IV paralogs were identified. Sequence analysis revealed conserved structural features including the common 33-codon block and nine repeated motifs. Gene structure analysis identified four exons and three introns except for apoA-IVAa1 (with only 3 exons). Synteny analysis showed that the four paralogs were structured into two clusters (cluster A containing apoA-IVAa1 and apoA-IVAa2 and cluster B with apoA-IVBa3 and apoA-IVBa4) linked to an apolipoprotein E. Phylogenetic analysis clearly separated the paralogs according to their cluster organization as well as revealed four subclades highly conserved in Acanthopterygii. Whole-mount analyses (WISH) in early larvae (0 and 1day post-hatch (dph)) showed that the four paralogs were mainly expressed in yolk syncytial layer surrounding the oil globules. Later, at 3 and 5dph, the four paralogs were mainly expressed in liver and intestine although with differences in their relative abundance and temporal expression patterns. Diet supply triggered the intensity of WISH signals in the intestine of the four paralogs. Quantification of mRNA abundance by qPCR using whole larvae only detected the induction by diet at 5dph. Moreover, transcript levels increased progressively with age except for apoA-IVAa2, which appeared as a low-expressed isoform. Expression analysis in juvenile tissues confirmed that the four paralogs were mainly expressed in liver and intestine and secondary in other tissues. The role of these ApoA-IV genes in lipid transport and the possible role of apoA-IVAa2 as a regulatory form are discussed.

Dietary vegetable oils: effects on the expression of immune-related genes in Senegalese sole (Solea senegalensis) intestine

The decreased availability of fish oil, traditionally used as oil source in marine aquafeeds, has lead to the search for alternatives oils. Vegetable oils (VO) are being extensively used as lipid sources in marine fish diets, inducing an imbalance on certain dietary fatty acids. Alteration on the dietary ratio of w-6/w-3 has been described to have detrimental effects on fish immunity. Senegalese sole has high susceptibility to stress and diseases, and little is known on the effects of dietary VO on its immunity. In this study, Senegalese sole juveniles were fed diets (56% crude protein, 12% crude lipid) containing linseed (100LO), soybean (100SO) or fish (100FO) oils as unique oil source. Growth, cortisol and intestinal fatty acid composition were determined after 90 days. Moreover, at the final of the experiment a stress test (5 min of net chasing) was carried out. To evaluate the effect of diets and stress on intestine immunology, expression profiles of a set of 53 immune-related genes using RT-qPCR was also performed. The use of VO did not induced changes in fish growth, but affected fatty acid profile of intestine and expression of immune-related genes. The use of SO (rich in n-6 fatty acids) induced an over-expression of those genes related to complement pathway, recognizing pathogen associated to molecular patterns, defensive response against bacteria, defensive response against viruses, antigen differentiation, cytokines and their receptors. This general over-expression could indicate an activation of inflammatory processes in fish gut. When a stress was applied, a decrease of mRNA levels of different immune-related genes with respect to the unstressed control could be observed in fish fed 100FO. However, fish fed 100LO, with a higher ALA/LA ratio, seemed to ameliorate the effects of combined effects of FO substitution plus stressful situation whereas fish fed 100SO did not show this type of response.

ReprOlive: a database with linked data for the olive tree (Olea europaea L.) reproductive transcriptome

Plant reproductive transcriptomes have been analyzed in different species due to the agronomical and biotechnological importance of plant reproduction. Here we presented an olive tree reproductive transcriptome database with samples from pollen and pistil at different developmental stages, and leaf and root as control vegetative tissues http://reprolive.eez.csic.es). It was developed from 2,077,309 raw reads to 1,549 Sanger sequences. Using a pre-defined workflow based on open-source tools, sequences were pre-processed, assembled, mapped, and annotated with expression data, descriptions, GO terms, InterPro signatures, EC numbers, KEGG pathways, ORFs, and SSRs. Tentative transcripts (TTs) were also annotated with the corresponding orthologs in Arabidopsis thaliana from TAIR and RefSeq databases to enable Linked Data integration. It results in a reproductive transcriptome comprising 72,846 contigs with average length of 686 bp, of which 63,965 (87.8%) included at least one functional annotation, and 55,356 (75.9%) had an ortholog. A minimum of 23,568 different TTs was identified and 5,835 of them contain a complete ORF. The representative reproductive transcriptome can be reduced to 28,972 TTs for further gene expression studies. Partial transcriptomes from pollen, pistil, and vegetative tissues as control were also constructed. ReprOlive provides free access and download capability to these results. Retrieval mechanisms for sequences and transcript annotations are provided. Graphical localization of annotated enzymes into KEGG pathways is also possible. Finally, ReprOlive has included a semantic conceptualisation by means of a Resource Description Framework (RDF) allowing a Linked Data search for extracting the most updated information related to enzymes, interactions, allergens, structures, and reactive oxygen species.