Exploring the larval transcriptome of the common sole (Solea solea L.)

The multi-level regulation of clownfish metamorphosis by thyroid hormones

Most marine organisms have a biphasic life cycle during which pelagic larvae transform into radically different juveniles. In vertebrates, the role of thyroid hormones (THs) in triggering this transition is well known, but how the morphological and physiological changes are integrated in a coherent way with the ecological transition remains poorly explored. To gain insight into this question, we performed an integrated analysis of metamorphosis of a marine teleost, the false clownfish (Amphiprion ocellaris). We show how THs coordinate a change in color vision as well as a major metabolic shift in energy production, highlighting how it orchestrates this transformation. By manipulating the activity of liver X regulator (LXR), a major regulator of metabolism, we also identify a tight link between metabolic changes and metamorphosis progression. Strikingly, we observed that these regulations are at play in the wild, explaining how hormones coordinate energy needs with available resources during the life cycle.

Rapid physiological and transcriptomic changes associated with oxygen delivery in larval anemonefish suggest a role in adaptation to life on hypoxic coral reefs

Connectivity of coral reef fish populations relies on successful dispersal of a pelagic larval phase. Pelagic larvae must exhibit high swimming abilities to overcome ocean and reef currents, but once settling onto the reef, larvae transition to endure habitats that become hypoxic at night. Therefore, coral reef fish larvae must rapidly and dramatically shift their physiology over a short period of time. Taking an integrative, physiological approach, using swimming respirometry, and examining hypoxia tolerance and transcriptomics, we show that larvae of cinnamon anemonefish (Amphiprion melanopus) rapidly transition between "physiological extremes" at the end of their larval phase. Daily measurements of swimming larval anemonefish over their entire early development show that they initially have very high mass-specific oxygen uptake rates. However, oxygen uptake rates decrease midway through the larval phase. This occurs in conjunction with a switch in haemoglobin gene expression and increased expression of myoglobin, cytoglobin, and neuroglobin, which may all contribute to the observed increase in hypoxia tolerance. Our findings indicate that critical ontogenetic changes in the gene expression of oxygen-binding proteins may underpin the physiological mechanisms needed for successful larval recruitment to reefs.

Gene Expression Profiles of the Immuno-Transcriptome in Equine Asthma

BACKGROUND

Mild equine asthma (MEA) and severe equine asthma (SEA) are two of the most frequent equine airway inflammatory diseases, but knowledge about their pathogenesis is limited. The goal of this study was to investigate gene expression differences in the respiratory tract of MEA- and SEA-affected horses and their relationship with clinical signs.

METHODS

Clinical examination and endoscopy were performed in 8 SEA- and 10 MEA-affected horses and 7 healthy controls. Cytological and microbiological analyses of bronchoalveolar lavage (BAL) fluid were performed. Gene expression profiling of BAL fluid was performed by means of a custom oligo-DNA microarray.

RESULTS

In both MEA and SEA, genes involved in the genesis, length, and motility of respiratory epithelium cilia were downregulated. In MEA, a significant overexpression for genes encoding inflammatory mediators was observed. In SEA, transcripts involved in bronchoconstriction, apoptosis, and hypoxia pathways were significantly upregulated, while genes involved in the formation of the protective muco-protein film were underexpressed. The SEA group also showed enrichment of gene networks activated during human asthma.

CONCLUSIONS

The present study provides new insight into equine asthma pathogenesis, representing the first step in transcriptomic analysis to improve diagnostic and therapeutic approaches for this respiratory disease.

Transcriptomic responses and apoptosis in larval red drum (Sciaenops ocellatus) co-exposed to crude oil and ultraviolet (UV) radiation

Ultraviolet (UV) radiation can significantly increase the toxicity of polycyclic aromatic hydrocarbons (PAHs) in crude oil to early life stage (ELS) fishes through photo-induced /photo-enhanced toxicity. However, little is known about the sub-lethal effects and mechanisms of photo-induced PAH toxicity in ELS fishes. The present study investigated apoptosis and global transcriptomic effects in larval red drum (Sciaenops ocellatus) (24-72 h post-fertilization) following co-exposure to oil (0.29-0.30 μg/L ∑PAH₅₀) and UV. Apoptosis was quantified using the TUNEL assay, and transcriptomic effects were assessed using RNA sequencing analysis. Apoptotic fluorescence was greatest in the eyes and skin following 24 and 48 h co-exposure to oil and UV, indicating photo-induced toxicity. Consistent with these phenotypic responses, pathways associated with phototransduction, eye development, and dermatological disease were among the top predicted pathways impacted. The present study is the first to provide global transcriptomic analysis of UV and oil co-exposure in an ELS fish.

Transcriptome Sequencing Analysis Reveals Dynamic Changes in Major Biological Functions during the Early Development of Clearhead Icefish, Protosalanx chinensis

Early development, when many important developmental events occur, is a critical period for fish. However, research on the early development of clearhead icefish is very limited, especially in molecular research. In this study, we aimed to explore the dynamic changes in the biological functions of five key periods in clearhead icefish early development, namely the YL (embryonic), PM (first day after hatching), KK (fourth day after hatching), LC (seventh day after hatching), and SL (tenth day after hatching) stages, through transcriptome sequencing and different analysis strategies. A trend expression analysis and an enrichment analysis revealed that the expression ofgenes encoding G protein-coupled receptors and their ligands, i.e., prss1_2_3, pomc, npy, npb, sst, rln3, crh, gh, and prl that are associated with digestion and feeding regulation gradually increased during early development. In addition, a weighted gene co-expression network analysis (WGCNA) showed that eleven modules were significantly associated with early development, among which nine modules were significantly positively correlated. Through the enrichment analysis and hub gene identification results of these nine modules, it was found that the pathways related to eye, bone, and heart development were significantly enriched in the YL stage, and the ccnd2, seh1l, kdm6a, arf4, and ankrd28 genes that are associated with cell proliferation and differentiation played important roles in these developmental processes; the pak3, dlx3, dgat2, and tas1r1 genes that are associated with jaw and tooth development, TG (triacylglycerol) synthesis, and umami amino acid receptors were identified as hub genes for the PM stage; the pathways associated with aerobic metabolism and unsaturated fatty acid synthesis were significantly enriched in the KK stage, with the foxk, slc13a2_3_5, ndufa5, and lsc2 genes playing important roles; the pathways related to visual perception were significantly enriched in the LC stage; and the bile acid biosynthetic and serine-type peptidase activity pathways were significantly enriched in the SL stage. These results provide a more detailed understanding of the processes of early development of clearhead icefish.

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.

Multi-organ transcriptomic landscape of Ambystoma velasci metamorphosis

Metamorphosis is a postembryonic developmental process that involves morphophysiological and behavioral changes, allowing organisms to adapt into a novel environment. In some amphibians, aquatic organisms undergo metamorphosis to adapt in a terrestrial environment. In this process, these organisms experience major changes in their circulatory, respiratory, digestive, excretory and reproductive systems. We performed a transcriptional global analysis of heart, lung and gills during diverse stages of Ambystoma velasci to investigate its metamorphosis. In our analyses, we identified eight gene clusters for each organ, according to the expression patterns of differentially expressed genes. We found 4064 differentially expressed genes in the heart, 4107 in the lung and 8265 in the gills. Among the differentially expressed genes in the heart, we observed genes involved in the differentiation of cardiomyocytes in the interatrial zone, vasculogenesis and in the maturation of coronary vessels. In the lung, we found genes differentially expressed related to angiogenesis, alveolarization and synthesis of the surfactant protein. In the case of the gills, the most prominent biological processes identified are degradation of extracellular matrix, apoptosis and keratin production. Our study sheds light on the transcriptional responses and the pathways modulation involved in the transformation of the facultative metamorphic salamander A. velasci in an organ-specific manner.

Ontogeny of Expression and Activity of Digestive Enzymes and Establishment of gh/igf1 Axis in the Omnivorous Fish Chelon labrosus

Simple Summary

Thick-lipped grey mullet (Chelon labrosus) feeds on the lowest trophic levels during adult stages, for which it is considered a viable candidate for an economically and environmentally sustainable aquaculture. Similar to most of marine fish species, C. labrosus produce a large number of eggs, leading to morphologically and anatomically larvae that are not completely mature and have to pass through substantial differentiation and development in their functional systems to acquire adult features. Therefore, the study of the development of digestive tract and of the growth regulation can provide useful information to adapt the feeding protocols and rearing conditions to the physiological requirements at each stage. This work aimed to evaluate the early ontogeny of key digestive enzymes and somatotropic factors at biochemical and/or transcriptional levels. Our results evidenced that maturation of the digestive system and acquisition of the adult mode of digestion occurs around 60 to 70 days post hatch (dph), when starch or other low-cost carbohydrate-based compounds could be used in formulated diets at increasing levels. Furthermore, our results implied an independent expression of the studied somatotropic genes during the first 40 dph and establishment of a functional growth hormone/insulin-like growth factor 1 axis from 50 dph onward.

Abstract

Thick-lipped grey mullet (Chelon labrosus) is a candidate for sustainable aquaculture due to its omnivorous/detritivorous feeding habit. This work aimed to evaluate its digestive and growth potentials from larval to early juvenile stages. To attain these objectives the activity of key digestive enzymes was measured from three until 90 days post hatch (dph). Expression of genes involved in digestion of proteins (try2, ctr, pga2, and atp4a), carbohydrates (amy2a), and lipids (cel and pla2g1b), together with two somatotropic factors (gh and igf1) were also quantified. No chymotrypsin or pepsin activities were detected. While specific activity of trypsin and lipase were high during the first 30 dph and declined afterward, amylase activity was low until 57 dph and increased significantly beyond that point. Expression of try2, ctr, amy2a, and cel increased continuously along development, and showed a peak at the end of metamorphosis. Expression of pla2g1b, pga2 and atp4a increased until the middle of metamorphosis and decreased afterwars. Most of these trends contrast the usual patterns in carnivorous species and highlight the transition from larvae, with high protein requirements, to post-larvae/juvenile stages, with omnivorous/detritivorous feeding preferences. Somatotropic genes, gh and igf1, showed approximately inverse expression patterns, suggesting the establishment of the Gh/Igf1 axis from 50 dph.

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.

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.

Dynamic transcriptome sequencing and analysis during early development in the bighead carp (Hypophthalmichthys nobilis)

BACKGROUND

Early development is a key process of the life history of fish. However, the relationship between the transcriptome and the dynamic regulation of early development is still uncharacterized in the bighead carp (Hypophthalmichthys nobilis). In the present study, we performed transcriptome analysis of six development stages in H. nobilis, aiming to understand the dynamic molecular regulation of early development in this fish.

RESULTS

A total of 76,573 unigenes were assembled from clean sequence reads, with an average length of 1768 base. Among which, 41,742 (54.54%) unigenes were annotated to public protein databases, and an additional 59,014 simple sequence repeat (SSR) loci were identified among the unigenes. Furthermore, 30,199 differentially expressed transcripts (DETs) (fold change > 4 or < 0.25, and the false discovery rate FDR < 0.01) were observed in comparisons between the adjacent developmental stages, and nine expression patterns (profiles) were simulated using series-cluster analysis across six developmental stages. The unigenes expression level markedly increased after the DS1 stage (early blastula), and the numbers of DETs gradually decreased during subsequent development. The largest transcriptomic change (up- or down-regulated) was detected during the period from DS1 to DS2 (6-somite stage), which was enriched for many biological processes and metabolic pathways related to maternal to zygotic transition (MZT). Distinctly protein-protein interaction (PPI) networks were plotted for DETs during the period from DS1 to DS2. The genes (or proteins) from the same pathways were integrated together, and showed with obvious co-regulation patterns. In the series-cluster analysis, a remarkable profile of gene expression (profile_48) was identified that is probably related to the hatching during H. nobilis development, and the strict co-expression of a hatching enzyme gene (hce1) with 33 other annotated genes was identified from this profile.

CONCLUSIONS

The results indicated that strict dynamic regulation occurs during the early development in H. nobilis, especially in embryogenesis before hatching. This study provides valuable new information and transcriptomic resources related to H. nobilis early development, and for certain events such as MZT and hatching.

Experimental infection model for vibriosis in Dover sole (Solea solea) larvae as an aid in studying its pathogenesis and alternative treatments

Severe economic losses due to diseases in marine larviculture may be linked to vibriosis. To better understand the pathogenesis of vibriosis and evaluate new ways to prevent and combat this important disease, there is a great need for reliable and reproducible experimental infection models. The present study aimed at developing a challenge model for vibriosis in Dover sole larvae and testing its applicability to study the effect of the probiotic treatment. For that purpose, larvae were challenged at 10 days post hatching with Vibrio anguillarum WT, V. anguillarum HI610 or V. harveyi WT. Following administration of V. anguillarum WT via immersion at 1 × 10⁷ colony forming units/mL, a larval mortality of 50% was observed at 17 days post-inoculation. In a next step, the probiotic potential of 371 isolates retrieved from Dover sole was assessed by screening for their inhibitory effects against Vibrio spp. and absence of haemolytic activity. One remaining isolate (V. proteolyticus) and V. lentus, known for its protective characteristics in seabass larvae, were further tested in vivo by means of the pinpointed experimental infection model. Neither isolate provided via the water or feed proved to be protective for the Dover sole larvae against challenge with V. anguillarum WT. This developed challenge model constitutes a firm basis to expedite basic and applied research regarding the pathogenesis and treatment of vibriosis as well as for studying the impact of (a)biotic components on larval health.

Larval Red Drum (Sciaenops ocellatus) Sublethal Exposure to Weathered Deepwater Horizon Crude Oil: Developmental and Transcriptomic Consequences

The Deepwater Horizon (DWH) incident resulted in extensive oiling of the pelagic zone and shoreline habitats of many commercially important fish species. Exposure to the water-accommodated fraction (WAF) of oil from the spill causes developmental toxicity through cardiac defects in pelagic fish species. However, few studies have evaluated the effects of the oil on near-shore estuarine fish species such as red drum (Sciaenops ocellatus). Following exposure to a certified weathered slick oil (4.74 μg/L ∑PAH₅₀) from the DWH event, significant sublethal impacts were observed ranging from impaired nervous system development [average 17 and 22% reductions in brain and eye area at 48 h postfertilization (hpf), respectively] to abnormal cardiac morphology (100% incidence at 24, 48, and 72 hpf) in red drum larvae. Consistent with the phenotypic responses, significantly differentially expressed transcripts, enriched gene ontology, and altered functions and canonical pathways predicted adverse outcomes in nervous and cardiovascular systems, with more pronounced changes at later larval stages. Our study demonstrated that the WAF of weathered slick oil of DWH caused morphological abnormalities predicted by a suite of advanced bioinformatic tools in early developing red drum and also provided the basis for a better understanding of molecular mechanisms of crude oil toxicity in fish.

Transcriptome Analysis Reveals Increases in Visceral Lipogenesis and Storage and Activation of the Antigen Processing and Presentation Pathway during the Mouth-Opening Stage in Zebrafish Larvae

The larval phase of the fish life cycle has the highest mortality, particularly during the transition from endogenous to exogenous feeding. However, the transcriptional events underlying these processes have not been fully characterized. To understand the molecular mechanisms underlying mouth-opening acclimation, RNA-seq was used to investigate the transcriptional profiles of the endogenous feeding, mixed feeding and exogenous feeding stages of zebrafish larvae. Differential expression analysis showed 2172 up-regulated and 2313 down-regulated genes during this stage. Genes associated with the assimilation of exogenous nutrients such as the arachidonic acid metabolism, linoleic acid metabolism, fat digestion and absorption, and lipogenesis were activated significantly, whereas dissimilation including the cell cycle, homologous recombination, and fatty acid metabolism were inhibited, indicating a physiological switch for energy storage occurred during the mouth-opening stage. Moreover, the immune recognition involved in the antigen processing and presentation pathway was activated and nutritional supply seemed to be required in this event confirmed by qPCR. These results suggested the energy utilization during the mouth-opening stage is more tended to be reserved or used for some important demands, such as activity regulation, immune defense, and lipid deposition, instead of rapid growth. The findings of this study are important for understanding the physiological switches during the mouth-opening stage.

Developmental transcriptomic analyses for mechanistic insights into critical pathways involved in embryogenesis of pelagic mahi-mahi (Coryphaena hippurus)

Mahi-mahi (Coryphaena hippurus) is a commercially and ecologically important species of fish occurring in tropical and temperate waters worldwide. Understanding early life events is crucial for predicting effects of environmental stress, which is largely restricted by a lack of genetic resources regarding expression of early developmental genes and regulation of pathways. The need for anchoring developmental stages to transcriptional activities is highlighted by increasing evidence on the impacts of recurrent worldwide oil spills in this sensitive species during early development. By means of high throughput sequencing, we characterized the developmental transcriptome of mahi-mahi at three critical developmental stages, from pharyngula embryonic stage (24 hpf) to 48 hpf yolk-sac larva (transition 1), and to 96 hpf free-swimming larva (transition 2). With comparative analysis by multiple bioinformatic tools, a larger number of significantly altered genes and more diverse gene ontology terms were observed during transition 2 than transition 1. Cellular and tissue development terms were more significantly enriched in transition 1, while metabolism related terms were more enriched in transition 2, indicating a switch progressing from general embryonic development to metabolism during the two transitions. Special focus was given on the most significant common canonical pathways (e.g. calcium signaling, glutamate receptor signaling, cAMP response element-binding protein signaling, cardiac β-adrenergic signaling, etc.) and expression of developmental genes (e.g. collagens, myosin, notch, glutamate metabotropic receptor etc.), which were associated with morphological changes of nervous, muscular, and cardiovascular system. These data will provide an important basis for understanding embryonic development and identifying molecular mechanisms of abnormal development in fish species.

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.

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.

Developmental transcriptome analysis and identification of genes involved in formation of intestinal air-breathing function of Dojo loach, Misgurnus anguillicaudatus

Dojo loach, Misgurnus anguillicaudatus is a freshwater fish species of the loach family Cobitidae, using its posterior intestine as an accessory air-breathing organ. Little is known about the molecular regulatory mechanisms in the formation of intestinal air-breathing function of M. anguillicaudatus. Here high-throughput sequencing of mRNAs was performed from six developmental stages of posterior intestine of M. anguillicaudatus: 4-Dph (days post hatch) group, 8-Dph group, 12-Dph group, 20-Dph group, 40-Dph group and Oyd (one-year-old) group. These six libraries were assembled into 81300 unigenes. Totally 40757 unigenes were annotated. Subsequently, 35291 differentially expressed genes (DEGs) were scanned among different developmental stages and clustered into 20 gene expression profiles. Finally, 15 key pathways and 25 key genes were mined, providing potential targets for candidate gene selection involved in formation of intestinal air-breathing function in M. anguillicaudatus. This is the first report of developmental transcriptome of posterior intestine in M. anguillicaudatus, offering a substantial contribution to the sequence resources for this species and providing a deep insight into the formation mechanism of its intestinal air-breathing function. This report demonstrates that M. anguillicaudatus is a good model for studies to identify and characterize the molecular basis of accessory air-breathing organ development in fish.

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.

Sequencing and characterization of a multi-organ Arctic charr transcriptome: A toolbox for investigating polymorphism and seasonal life in a high Arctic fish

The Arctic charr (Salvelinus alpinus L.) inhabits fresh water ecosystems of the high North. The species has developed a strong phenotypic plasticity and variability in life history characteristics which has made this species an attractive model for investigations on phenotype plasticity, morph formation and ecological speciation. Further, the extreme seasonal variations in environmental conditions (e.g. food availability) in the high North induce seasonal changes in phenotype, which require precise timing mechanisms and physiological preparations. Individual gating of life-history strategies (e.g. formation of resident and sea-migrating morphs) and transitions (e.g. maturation) depends on conditional traits (size/energy status) at specific assessment time windows, and complex neuroendocrine regulation, which so far is poorly understood. In the absence of a reference genome, and in order to facilitate the investigation of the complex biological mechanisms of this unique fish model, the present study reveals a reference transcriptome for the Arctic charr. Using Roche 454 GS FLX+, we targeted various organs being either at the crossroads of many key pathways (neuroendocrine, metabolic, behavioral), of different ontological origins or displaying complementary physiological functions. The assemblage yielded 34,690 contigs greater than 1000bp with an average length (1690bp) and annotation rate (52%) within the range, or even higher, than what has been previously obtained with other teleost de novo transcriptomes. We dramatically improve the publically available transcript data on this species that may indeed be useful for various disciplines, from basic research to applied aspects related to conservation issues and aquaculture.

Transcriptome profiling reveals that feeding wild zooplankton to larval Atlantic cod (Gadus morhua) influences suites of genes involved in oxidation-reduction, mitosis, and selenium homeostasis

BACKGROUND

Larval nutrition and growth are key issues for wild and cultured cod. While it was shown previously that larval cod fed wild zooplankton grow faster than those fed only rotifers, the mechanisms involved in this enhanced growth are not completely understood. We used microarrays to identify larval cod transcripts that respond to feeding with small amounts of wild zooplankton (5-10 % of live prey items). The larval transcriptome was compared between 3 treatment groups [fed rotifers (RA), rotifers with protein hydrolysate (RA-PH), or rotifers with zooplankton (RA-Zoo)] at 9-10 mm length [26-30 days post-hatch (dph)] to identify a robust suite of zooplankton-responsive genes (i.e. differentially expressed between RA-Zoo and both other groups).

RESULTS

The microarray experiment identified 147 significantly up-regulated and 156 significantly down-regulated features in RA-Zoo compared with both RA and RA-PH. Gene ontology terms overrepresented in the RA-Zoo responsive gene set included "response to selenium ion" and several related to cell division and oxidation-reduction. Ten selenoprotein-encoding genes, and 2 genes involved in thyroid hormone generation, were up-regulated in RA-Zoo compared with both other groups. Hierarchical clustering of RA-Zoo responsive genes involved in oxidation-reduction and selenium homeostasis demonstrated that only the zooplankton treatment had a considerable and consistent impact on the expression of these genes. Fourteen microarray-identified genes were selected for QPCR involving 9-13 mm larvae, and 13 of these were validated as differentially expressed between RA-Zoo and both other groups at ~9 mm. In contrast, in age-matched (34-35 dph; ~11 mm RA and RA-PH, ~13 mm RA-Zoo) and size-matched (~13 mm) older larvae, only 2 and 3 genes, respectively, showed the same direction of RA-Zoo-responsive change as in ~9 mm larvae.

CONCLUSIONS

The modulation of genes involved in selenium binding, redox homeostasis, and thyroid hormone generation in ~9 mm RA-Zoo larvae in this study may be in response to the relatively high levels of selenium, iodine, and LC-PUFA (potentially causing oxidative stress) in zooplankton. Nonetheless, only a subset of zooplankton-responsive genes in ~9 mm larvae remained so in older larvae, suggesting that the observed transcriptome changes are largely involved in initiating the period of growth enhancement.

In- and outdoor reproduction of first generation common sole Solea solea under a natural photothermal regime: Temporal progression of sexual maturation assessed by monitoring plasma steroids and gonadotropin mRNA expression

Reproduction of many temperate fishes is seasonal and maturation and spawning of gametes are under photothermal control. Reproductive success of first generation (G1) common sole Solea solea in captivity has been low. In this study, the sexual maturation status has been assessed during the prespawning months in G1 sole that were housed (a) outdoor under the natural photoperiod and temperature, or (b) indoor under artificial photothermal induction. Maturation was assessed in male and female G1 broodstock in November as controls, after which the remaining population was divided over two outdoor flow-through tanks placed in a pond and two indoor recirculating aquaculture system (RAS) tanks. Subsequently, maturation status (gonadosomatic index GSI and plasma levels of testosterone T and 17β-estradiol E2) was assessed in one tank for each condition in January, February and during spawning in early April, while fish in the other tank were not disturbed in achieving reproductive success. Quantitative real-time PCR was performed to determine species-specific gonadotropin mRNA expression in females. Successful G1 spawning and egg fertilisation occurred in all experimental tanks. Gonadal development was similar under both conditions. Higher E2 and T levels were found in indoor housed females. Gonadotropin expression revealed similar profiles between outdoor and indoor housed females. G1 sole could be reproduced in the outdoor tanks under the natural photoperiod and in the indoor tanks under artificial simulation of this regime that includes a potentially crucial chilling period of 2-3 months at 5-7 °C.

Orchestrating change: The thyroid hormones and GI-tract development in flatfish metamorphosis

Metamorphosis in flatfish (Pleuronectiformes) is a late post-embryonic developmental event that prepares the organism for the larval-to-juvenile transition. Thyroid hormones (THs) play a central role in flatfish metamorphosis and the basic elements that constitute the thyroid axis in vertebrates are all present at this stage. The advantage of using flatfish to study the larval-to-juvenile transition is the profound change in external morphology that accompanies metamorphosis making it easy to track progression to climax. This important lifecycle transition is underpinned by molecular, cellular, structural and functional modifications of organs and tissues that prepare larvae for a successful transition to the adult habitat and lifestyle. Understanding the role of THs in the maturation of organs and tissues with diverse functions during metamorphosis is a major challenge. The change in diet that accompanies the transition from a pelagic larvae to a benthic juvenile in flatfish is associated with structural and functional modifications in the gastrointestinal tract (GI-tract). The present review will focus on the maturation of the GI-tract during metamorphosis giving particular attention to organogenesis of the stomach a TH triggered event. Gene transcripts and biological processes that are associated with GI-tract maturation during Atlantic halibut metamorphosis are identified. Gene ontology analysis reveals core biological functions and putative TH-responsive genes that underpin TH-driven metamorphosis of the GI-tract in Atlantic halibut. Deciphering the specific role remains a challenge. Recent advances in characterizing the molecular, structural and functional modifications that accompany the appearance of a functional stomach in Atlantic halibut are considered and future research challenges identified.

Microarray: a global analysis of biomineralization-related gene expression profiles during larval development in the pearl oyster, Pinctada fucata

BACKGROUND

The molluscan Pinctada fucata is an important pearl-culturing organism to study biomineralization mechanisms. Several biomineralization-related genes play important roles regulating shell formation, but most previous work has focused only on their functions in adult oysters. Few studies have investigated biomineralization during larval development, when the shell is initially constructed and formed until the juvenile stage in dissoconch shells. Here, we report, for the first time, a global gene analysis during larval development of P. fucata based on a microarray and reveal the relationships between biomineralization-related genes and the shell formation process.

RESULTS

Based on the P. fucata mantle transcriptome, 58,940 probes (60 nt), representing 58,623 transcripts, were synthesized. The gene expression profiles of the fertilized egg, trochophore, D-shaped, and umbonal stage larvae, as well as juveniles were analyzed by microarray performance. The expression patterns of the biomineralization-related genes changed corresponding to their regulatory function during shell formation. Matrix proteins chitin synthase and PFMG2 were highly expressed at the D-shaped stage, whereas PFMG6, PFMG8 and PfN23 were significantly up-regulated at the umbonal stage, indicating different roles regulating the formation of either periostracum, Prodissoconch I or Prodissoconch II shells. However, the majority of matrix proteins were expressed at high levels at the juvenile stage, and the shells comprised both an aragonitic nacreous layer and a calcitic prismatic layer as adults. We also identified five new genes that were significantly up-regulated in juveniles. These genes were expressed particularly in the mantle and coded for secreted proteins with tandem-arranged repeat units, as most matrix proteins. RNAi knockdown resulted in disrupted nacreous and prismatic shell layers, indicating their potential roles in shell formation.

CONCLUSIONS

Our results add a global perspective on larval expression patterns of P. fucata genes and propose a mechanism of how biomineralization-related genes regulate the larval shell formation process. These results increase knowledge about biomineralization-related genes and highlight new aspects of shell formation mechanisms.

In-depth transcriptome analysis of Coilia ectenes, an important fish resource in the Yangtze River: de novo assembly, gene annotation

Coilia ectenes is an important teleost species in the Yangtze River and a model organism that can be used to study the protection of fish resources. In this report, we performed de novo transcriptome sequencing of ten cDNA libraries from the brain, gill, heart, intestine, kidney, liver, muscle, stomach, ovary, and testis tissues. A total of 352 million raw reads of 100 base pairs were generated, and 130,113 transcripts, corresponding to 65,350 non-redundant transcripts, with a mean length of 1520 bp, were assembled. BLASTx-based gene annotation (E-value<1 × 10(-5)) allowed the identification of 73,900 transcripts against at least one of four databases, including the NCBI non-redundant database, the GO database, the COG database, and the KEGG database. Our study provides a valuable resource for C. ectenes genomic and transcriptomic data that will facilitate future functional studies of C. ectenes.

Fatty acid composition of eggs and its relationships to egg and larval viability from domesticated common sole (Solea solea) breeders

The study of lipids and fatty acids (FAs) has been used in the assessment of egg quality because their composition can influence the fertilization rate, hatching, survival and growth of marine fish larvae. For these reasons, the lipid content (TL) and fatty acid composition of common sole (Solea solea) eggs were measured and correlated to egg and larval viability parameters throughout an entire reproductive season. Seventeen batches of fertile eggs obtained from natural spawning of captive breeders were characterized for the TL, FA profile, hatching rate (HR) and survival rate of larvae (SR) at 0-6 days post-hatching (dph). The egg FA composition reflected the composition of the feed supplied to the broodstock during summer and autumn (before and during vitellogenesis) rather than that supplied during the spawning season. In general, the egg FA profile showed minimal differences among the early-, mid- and late-spawning periods (possibly due to the change of the diet and/or water temperature) indicating that it is possible to obtain a similar egg quality in terms of egg FA profile over 2 months of spawning. Saturated FAs and monounsaturated FAs (MUFA) were positively correlated with HR, while TL, 22 : 6n-3 (DHA), 20 : 4n-6 (ARA), polyunsaturated FAs of the (n-3) series (n-3 PUFA) and polyunsaturated FAs of the (n-6) series were negatively correlated (p ≤ 0.05). MUFA, 20 : 5n-3 (EPA), n-6/n-3 were positively correlated with SR, while DHA, n-3 PUFA, DHA/EPA were negatively correlated (p ≤ 0.05). In conclusion, the feed supplied before and during vitellogenesis has a major role in determining the egg FA profile in common sole. The relationships found between TL and FAs with egg and larval viability parameters differ from many other farmed marine fish species, which may suggest the need for a specific broodstock feed for this species.

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

Background

Senegalese sole (Solea senegalensis) and common sole (S. solea) are two economically and evolutionary important flatfish species both in fisheries and aquaculture. Although some genomic resources and tools were recently described in these species, further sequencing efforts are required to establish a complete transcriptome, and to identify new molecular markers. Moreover, the comparative analysis of transcriptomes will be useful to understand flatfish evolution.

Results

A comprehensive characterization of the transcriptome for each species was carried out using a large set of Illumina data (more than 1,800 millions reads for each sole species) and 454 reads (more than 5 millions reads only in S. senegalensis), providing coverages ranging from 1,384x to 2,543x. After a de novo assembly, 45,063 and 38,402 different transcripts were obtained, comprising 18,738 and 22,683 full-length cDNAs in S. senegalensis and S. solea, respectively. A reference transcriptome with the longest unique transcripts and putative non-redundant new transcripts was established for each species. A subset of 11,953 reference transcripts was qualified as highly reliable orthologs (>97% identity) between both species. A small subset of putative species-specific, lineage-specific and flatfish-specific transcripts were also identified. Furthermore, transcriptome data permitted the identification of single nucleotide polymorphisms and simple-sequence repeats confirmed by FISH to be used in further genetic and expression studies. Moreover, evidences on the retention of crystallins crybb1, crybb1-like and crybb3 in the two species of soles are also presented. Transcriptome information was applied to the design of a microarray tool in S. senegalensis that was successfully tested and validated by qPCR. Finally, transcriptomic data were hosted and structured at SoleaDB.

Conclusions

Transcriptomes and molecular markers identified in this study represent a valuable source for future genomic studies in these economically important species. Orthology analysis provided new clues regarding sole genome evolution indicating a divergent evolution of crystallins in flatfish. The design of a microarray and establishment of a reference transcriptome will be useful for large-scale gene expression studies. Moreover, the integration of transcriptomic data in the SoleaDB will facilitate the management of genomic information in these important species.

Electronic supplementary material

The online version of this article (doi:10.1186/1471-2164-15-952) contains supplementary material, which is available to authorized users.

Population genomics of marine fishes: next-generation prospects and challenges

Over the past few years, technological advances have facilitated giant leaps forward in our ability to generate genome-wide molecular data, offering exciting opportunities for gaining new insights into the ecology and evolution of species where genomic information is still limited. Marine fishes are valuable organisms for advancing our understanding of evolution on historical and contemporary time scales, and here we highlight areas in which research on these species is likely to be particularly important in the near future. These include possibilities for gaining insights into processes on ecological time scales, identifying genomic signatures associated with population divergence under gene flow, and determining the genetic basis of phenotypic traits. We also consider future challenges pertaining to the implementation of genome-wide coverage through next-generation sequencing and genotyping methods in marine fishes. Complications associated with fast decay of linkage disequilibrium, as expected for species with large effective population sizes, and the possibility that adaptation is associated with both soft selective sweeps and polygenic selection, leaving complex genomic signatures in natural populations, are likely to challenge future studies. However, the combination of high genome coverage and new statistical developments offers promising solutions. Thus, the next generation of studies is likely to truly facilitate the transition from population genetics to population genomics in marine fishes. This transition will advance our understanding of basic evolutionary processes and will offer new possibilities for conservation and management of valuable marine resources.

Identification of hypoxia-regulated genes in the liver of common sole (Solea solea) fed different dietary lipid contents

Coastal systems could be affected by hypoxic events brought about by global change. These areas are essential nursery habitats for several fish species including the common sole (Solea solea L.). Tolerance of fish to hypoxia depends on species and also on their physiological condition and nutritional status. Indeed, high dietary lipid content has been recently shown to negatively impact the resistance of sole to a severe hypoxic challenge. In order to study the molecular mechanisms involved in the early response to hypoxic stress, the present work examined the hepatic transcriptome in common sole fed diets with low and high lipid content, exposed to severe hypoxia. The activity of AMP-activated protein kinase (AMPK) was also investigated through the quantification of threonine-172 phosphorylation in the alpha subunit. The results show that hypoxia consistently regulates several actors involved in energy metabolism pathways and particularly AMPKα, as well as some involved in cell growth and maintenance or unfolded protein response. Our findings reveal that (1) the expression of genes involved in biological processes with high energy cost or implicated in aerobic ATP synthesis was down-regulated by hypoxia, contrary to genes involved in neoglucogenesis or in angiogenesis, (2) the consumption of high lipid induced regulation of metabolic pathways going against this energy saving, and (3) this control was fine-tuned by the regulation of several transcriptomic factors. These results provide insight into the biological processes involved in the hepatic response to hypoxic stress and underline the negative impact of high lipid consumption on the tolerance of common sole to hypoxia.

Generation and characterization of the sea bass Dicentrarchus labrax brain and liver transcriptomes

The sea bass Dicentrarchus labrax is the center of interest of an increasing number of basic or applied research investigations, even though few genomic or transcriptomic data is available. Current public data only represent a very partial view of its transcriptome. To fill this need, we characterized brain and liver transcriptomes in a generalist manner that would benefit the entire scientific community. We also tackled some bioinformatics questions, related to the effect of RNA fragment size on the assembly quality. Using Illumina RNA-seq, we sequenced organ pools from both wild and farmed Atlantic and Mediterranean fishes. We built two distinct cDNA libraries per organ that only differed by the length of the selected mRNA fragments. Efficiency of assemblies performed on either or both fragments size differed depending on the organ, but remained very close reflecting the quality of the technical replication. We generated more than 19,538Mbp of data. Over 193million reads were assembled into 35,073 contigs (average length=2374bp; N50=3257). 59% contigs were annotated with SwissProt, which corresponded to 12,517 unique genes. We compared the Gene Ontology (GO) contig distribution between the sea bass and the tilapia. We also looked for brain and liver GO specific signatures as well as KEGG pathway coverage. 23,050 putative micro-satellites and 134,890 putative SNPs were identified. Our sampling strategy and assembly pipeline provided a reliable and broad reference transcriptome for the sea bass. It constitutes an indisputable quantitative and qualitative improvement of the public data, as it provides 5 times more base pairs with fewer and longer contigs. Both organs present unique signatures consistent with their specific physiological functions. The discrepancy in fragment size effect on assembly quality between organs lies in their difference in complexity and thus does not allow prescribing any general strategy. This information on two key organs will facilitate further functional approaches.