Fast skeletal muscle transcriptome of the gilthead sea bream (Sparus aurata) determined by next generation sequencing

Muscle Transcriptome Sequencing Revealed Thermal Stress-Responsive Regulatory Genes in Farmed Rohu, Labeo rohita (Hamilton, 1822)

Rohu, Labeo rohita, is one of the most important aquaculture species in the Indian subcontinent. Understanding the molecular-level physiological responses to thermal stress or climate change is essential. In the present work, transcriptome sequencing was carried out in the muscle tissue of the rohu in response to heat stress (35 °C) in comparison with the control (28 °C). A total of 125 Gb of sequence data was generated, and the raw-reads were filtered and trimmed, which resulted in 484 million quality reads. Reference-based assembly of reads was performed using L. rohita genome, and a total of 90.17% of reads were successfully mapped. A total of 37,462 contigs were assembled with an N50 value of 1854. The differential expression analysis revealed a total of 107 differentially expressed genes (DEGs) (15 up-, 37 down-, and 55 neutrally regulated) as compared to the control group (Log2FC > 2, P < 0.05). Gene enrichment analysis of DEGs indicates that transcripts were associated with molecular, biological, and cellular activities. The randomly selected differentially expressed transcripts were validated by RT-qPCR and found consistent expression patterns in line with the RNA-seq data. Several transcripts such as SERPINE1(HSP47), HSP70, HSP90alpha, Rano class II histocompatibility A beta, PGC-1 and ERR-induced regulator, proto-oncogene c-Fos, myozenin2, alpha-crystallin B chain-like protein, angiopoietin-like protein 8, and acetyl-CoA carboxylases have been identified in muscle tissue of rohu that are associated with stress/immunity. This study identified the key biomarker SERPINE1 (HSP47), which showed significant upregulation (~ 2- to threefold) in muscle tissue of rohu exposed to high temperature. This study can pave a path for the identification of stress-responsive biomarkers linked with thermal adaptations in the farmed carps.

The Emerging Role of Long Non-Coding RNAs in Development and Function of Gilthead Sea Bream (Sparus aurata) Fast Skeletal Muscle

Long non-coding RNAs (lncRNAs) are an emerging group of ncRNAs that can modulate gene expression at the transcriptional or translational levels. In the present work, previously published transcriptomic data were used to identify lncRNAs expressed in gilthead sea bream skeletal muscle, and their transcription levels were studied under different physiological conditions. Two hundred and ninety lncRNAs were identified and, based on transcriptomic differences between juveniles and adults, a total of seven lncRNAs showed potential to be important for muscle development. Our data suggest that the downregulation of most of the studied lncRNAs might be linked to increased myoblast proliferation, while their upregulation might be necessary for differentiation. However, with these data, as it is not possible to propose a formal mechanism to explain their effect, bioinformatic analysis suggests two possible mechanisms. First, the lncRNAs may act as sponges of myoblast proliferation inducers microRNAs (miRNAs) such as miR-206, miR-208, and miR-133 (binding energy MEF < -25.0 kcal). Secondly, lncRNA20194 had a strong predicted interaction towards the myod1 mRNA (ndG = -0.17) that, based on the positive correlation between the two genes, might promote its function. Our study represents the first characterization of lncRNAs in gilthead sea bream fast skeletal muscle and provides evidence regarding their involvement in muscle development.

Metabolic Remodeling in Skeletal Muscle Atrophy as a Therapeutic Target

Skeletal muscle is a highly responsive tissue, able to remodel its size and metabolism in response to external demand. Muscle fibers can vary from fast glycolytic to slow oxidative, and their frequency in a specific muscle is tightly regulated by fiber maturation, innervation, or external causes. Atrophic conditions, including aging, amyotrophic lateral sclerosis, and cancer-induced cachexia, differ in the causative factors and molecular signaling leading to muscle wasting; nevertheless, all of these conditions are characterized by metabolic remodeling, which contributes to the pathological progression of muscle atrophy. Here, we discuss how changes in muscle metabolism can be used as a therapeutic target and review the evidence in support of nutritional interventions and/or physical exercise as tools for counteracting muscle wasting in atrophic conditions.

The Estrogen-Responsive Transcriptome of Female Secondary Sexual Traits in the Gulf Pipefish

Sexual dimorphism often results from hormonally regulated trait differences between the sexes. In sex-role-reversed vertebrates, females often have ornaments used in mating competition that are expected to be under hormonal control. Males of the sex-role-reversed Gulf pipefish (Syngnathus scovelli) develop female-typical traits when they are exposed to estrogens. We aimed to identify genes whose expression levels changed during the development and maintenance of female-specific ornaments. We performed RNA-sequencing on skin and muscle tissue in male Gulf pipefish with and without exposure to estrogen to investigate the transcriptome of the sexually dimorphic ornament of vertical iridescent bands found in females and estrogen-exposed males. We further compared differential gene expression patterns between males and females to generate a list of genes putatively involved in the female secondary sex traits of bands and body depth. A detailed analysis of estrogen-receptor binding sites demonstrates that estrogen-regulated genes tend to have nearby cis-regulatory elements. Our results identified a number of genes that differed between the sexes and confirmed that many of these were estrogen-responsive. These estrogen-regulated genes may be involved in the arrangement of chromatophores for color patterning, as well as in the growth of muscles to achieve the greater body depth typical of females in this species. In addition, anaerobic respiration and adipose tissue could be involved in the rigors of female courtship and mating competition. Overall, this study generates a number of interesting hypotheses regarding the genetic basis of a female ornament in a sex-role-reversed pipefish.

Identification and characterization of parvalbumin-like protein in Trichophyton violaceum

Parvalbumins play crucial physiological roles in neuromuscular systems of vertebrates, such as cell-cycle, development of neurons, contraction of muscles, and regulation of intracellular calcium. To perform these neuromuscular functions, parvalbumin may be in associated with other proteins including calbindin, carbonic anhydrase, and cytochrome oxidase. Humans may show an IgE-specific hypersensitivity to parvalbumins after consumption of some distinct fish species. While this protein is abundant in fish muscles, literature review of publications related to fish parvalbumins, do not point to the presence of parvalbumins in eukaryotic microbes. In this study, we propose that distantly related parvalbumins may be found in some non-fish species. Bioinformatics studies such as multiple sequence alignment (MSA), phylogenetic analysis as well as molecular-based experiments indicate that, at least two parvalbumins sequences (UniProt IDs: A0A178F775 and A0A178F7E4) with EF-hand domains and Ca²⁺-binding sites could be identified in Trichophyton violaceum, a pathogenic fungal species. It was determined that both genes consisted of a single exon and encoded for parvalbumin proteins possessing conserved amino acid motifs. Antigenicity prediction revealed antigenic sites located in both sides of the Ca²⁺-binding site of the first EF-hand domain. Our phylogenetic analysis revealed that one of parvalbumins (UniProt ID: 0A178F775) can be evolved to other parvalbumins in T. violaceum (UniProt ID: A0A178F7E4) and fish species through evolutionary phenomenon. To confirm our in-silico findings, we designed three primer pairs to detect one of the T. violaceum parvalbumins (UniProt ID: A0A178F7E4) by polymerase chain reaction (PCR); one primer pair showed a strong and specific band in agarose gel electrophoresis. To evaluate the specificity of the method, the primers were tested on extracted DNA from Trichophyton rubrum and T. mentagrophytes. The results demonstrated that the evaluated parvalbumin gene (UniProt ID: A0A178F7E4) was T. violaceum-specific and this pathogenic fungus can be differentiated from T. rubrum and T. mentagrophytes through identification of parvalbumin genes. Further studies are necessary to unravel the biochemical and physiological functions of parvalbumins in T. violaceum.

Teleost Metamorphosis: The Role of Thyroid Hormone

In most teleosts, metamorphosis encompasses a dramatic post-natal developmental process where the free-swimming larvae undergo a series of morphological, cellular and physiological changes that enable the larvae to become a fully formed, albeit sexually immature, juvenile fish. In all teleosts studied to date thyroid hormones (TH) drive metamorphosis, being the necessary and sufficient factors behind this developmental transition. During metamorphosis, negative regulation of thyrotropin by thyroxine (T4) is relaxed allowing higher whole-body levels of T4 that enable specific responses at the tissue/cellular level. Higher local thyroid cellular signaling leads to cell-specific responses that bring about localized developmental events. TH orchestrate in a spatial-temporal manner all local developmental changes so that in the end a fully functional organism arises. In bilateral teleost species, the most evident metamorphic morphological change underlies a transition to a more streamlined body. In the pleuronectiform lineage (flatfishes), these metamorphic morphological changes are more dramatic. The most evident is the migration of one eye to the opposite side of the head and the symmetric pelagic larva development into an asymmetric benthic juvenile. This transition encompasses a dramatic loss of the embryonic derived dorsal-ventral and left-right axis. The embryonic dorsal-ventral axis becomes the left-right axis, whereas the embryonic left-right axis becomes, irrespectively, the dorsal-ventral axis of the juvenile animal. This event is an unparalleled morphological change in vertebrate development and a remarkable display of the capacity of TH-signaling in shaping adaptation and evolution in teleosts. Notwithstanding all this knowledge, there are still fundamental questions in teleost metamorphosis left unanswered: how the central regulation of metamorphosis is achieved and the neuroendocrine network involved is unclear; the detailed cellular and molecular events that give rise to the developmental processes occurring during teleost metamorphosis are still mostly unknown. Also in flatfish, comparatively little is still known about the developmental processes behind asymmetric development. This review summarizes the current knowledge on teleost metamorphosis and explores the gaps that still need to be challenged.

Dietary Creatine Supplementation in Gilthead Seabream (Sparus aurata) Increases Dorsal Muscle Area and the Expression of myod1 and capn1 Genes

Creatine (Cr) is an amino acid derivative with an important role in the cell as energy buffer that has been largely used as dietary supplement to increase muscle strength and lean body mass in healthy individuals and athletes. However, studies in fish are scarce. The aim of this work is to determine whether dietary Cr supplementation affects muscle growth in gilthead seabream (Sparus aurata) juveniles. Fish were fed ad libitum for 69 days with diets containing three increasing levels of creatine monohydrate (2, 5, and 8%) that were compared with a non-supplemented control (CTRL) diet. At the end of the trial, the fast-twist skeletal muscle growth dynamics (muscle cellularity) and the expression of muscle-related genes were evaluated. There was a general trend for Cr-fed fish to be larger and longer than those fed the CTRL, but no significant differences in daily growth index (DGI) were registered among dietary treatments. The dorsal cross-sectional muscle area (DMA) of fish fed Cr 5 and Cr 8% was significantly larger than that of fish fed CTRL. The groups supplemented with Cr systematically had a higher relative number of both small-sized (≤20 μm) and large-sized fibers (≥120 μm). Dorsal total fibers number was highest in fish fed 5% Cr. In fish supplemented with 5% Cr, the relative expression of myogenic differentiation 1 (myod1) increased almost four times compared to those fed the CTRL diet. The relative expression of calpain 3 (capn3) was highest in fish fed diets with 2% Cr supplementation, but did not differ significantly from those fed the CTRL or Cr 5%. The myod1 gene expression had a positive and significant correlation with that of capn1, capns1a, and capn3 expression. These results suggest that the observed modulation of gene expression was not enough to produce a significant alteration in muscle phenotype under the tested conditions, as a non-significant increase in muscle fiber diameter and higher total number of fiber was observed, but still resulted in increased DMA. Additional studies may be required in order to better clarify the effect of dietary Cr supplementation in fish, possibly in conjunction with induced resistance training.

Full genome survey and dynamics of gene expression in the greater amberjack Seriola dumerili

Background

Teleosts of the genus Seriola, commonly known as amberjacks, are of high commercial value in international markets due to their flesh quality and worldwide distribution. The Seriola species of interest to Mediterranean aquaculture is the greater amberjack (Seriola dumerili). This species holds great potential for the aquaculture industry, but in captivity, reproduction has proved to be challenging, and observed growth dysfunction hinders their domestication. Insights into molecular mechanisms may contribute to a better understanding of traits like growth and sex, but investigations to unravel the molecular background of amberjacks have begun only recently.

Findings

Illumina HiSeq sequencing generated a high-coverage greater amberjack genome sequence comprising 45 909 scaffolds. Comparative mapping to the Japanese yellowtail (Seriola quinqueriadiata) and to the model species medaka (Oryzias latipes) allowed the generation of in silico groups. Additional gonad transcriptome sequencing identified sex-biased transcripts, including known sex-determining and differentiation genes. Investigation of the muscle transcriptome of slow-growing individuals showed that transcripts involved in oxygen and gas transport were differentially expressed compared with fast/normal-growing individuals. On the other hand, transcripts involved in muscle functions were found to be enriched in fast/normal-growing individuals.

Conclusion

The present study provides the first insights into the molecular background of male and female amberjacks and of fast- and slow-growing fish. Therefore, valuable molecular resources have been generated in the form of a first draft genome and a reference transcriptome. Sex-biased genes, which may also have roles in sex determination or differentiation, and genes that may be responsible for slow growth are suggested.

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.

Trans-omics approaches used to characterise fish nutritional biorhythms in leopard coral grouper (Plectropomus leopardus)

Aquaculture is now a major supplier of fish, and has the potential to be a major source of protein in the future. Leopard coral groupers are traded in Asian markets as superior fish, and production via aquaculture has commenced. As feeding efficiency is of great concern in aquaculture, we sought to examine the metabolism of leopard coral groupers using trans-omics approaches. Metabolic mechanisms were comprehensively analysed using transcriptomic and metabolomic techniques. This study focused on the dynamics of muscular metabolites and gene expression. The omics data were discussed in light of circadian rhythms and fasting/feeding. The obtained data suggest that branched-chain amino acids played a role in energy generation in the fish muscle tissues during fasting. Moreover, glycolysis, TCA cycles, and purine metabolic substances exhibited circadian patterns, and gene expression also varied. This study is the first step to understanding the metabolic mechanisms of the leopard coral grouper.

Gene Expression Profiling in Fish Toxicology: A Review

In this review, we present an overview of transcriptomic responses to chemical exposures in a variety of fish species. We have discussed the use of several molecular approaches such as northern blotting, differential display reverse transcription-polymerase chain reaction (DDRT-PCR), suppression subtractive hybridization (SSH), real time quantitative PCR (RT-qPCR), microarrays, and next-generation sequencing (NGS) for measuring gene expression. These techniques have been mainly used to measure the toxic effects of single compounds or simple mixtures in laboratory conditions. In addition, only few studies have been conducted to examine the biological significance of differentially expressed gene sets following chemical exposure. Therefore, future studies should focus more under field conditions using a multidisciplinary approach (genomics, proteomics and metabolomics) to understand the synergetic effects of multiple environmental stressors and to determine the functional significance of differentially expressed genes. Nevertheless, recent developments in NGS technologies and decreasing costs of sequencing holds the promise to uncover the complexity of anthropogenic impacts and biological effects in wild fish populations.

Skeletal muscle fiber type: using insights from muscle developmental biology to dissect targets for susceptibility and resistance to muscle disease

Skeletal muscle fibers are classified into fiber types, in particular, slow twitch versus fast twitch. Muscle fiber types are generally defined by the particular myosin heavy chain isoforms that they express, but many other components contribute to a fiber's physiological characteristics. Skeletal muscle fiber type can have a profound impact on muscle diseases, including certain muscular dystrophies and sarcopenia, the aging-induced loss of muscle mass and strength. These findings suggest that some muscle diseases may be treated by shifting fiber type characteristics either from slow to fast, or fast to slow phenotypes, depending on the disease. Recent studies have begun to address which components of muscle fiber types mediate their susceptibility or resistance to muscle disease. However, for many diseases it remains largely unclear why certain fiber types are affected. A substantial body of work has revealed molecular pathways that regulate muscle fiber type plasticity and early developmental muscle fiber identity. For instance, recent studies have revealed many factors that regulate muscle fiber type through modulating the activity of the muscle regulatory transcription factor MYOD1. Future studies of muscle fiber type development in animal models will continue to enhance our understanding of factors and pathways that may provide therapeutic targets to treat muscle diseases. WIREs Dev Biol 2016, 5:518-534. doi: 10.1002/wdev.230 For further resources related to this article, please visit the WIREs website.

Integrative Transcriptome, Genome and Quantitative Trait Loci Resources Identify Single Nucleotide Polymorphisms in Candidate Genes for Growth Traits in Turbot

Growth traits represent a main goal in aquaculture breeding programs and may be related to adaptive variation in wild fisheries. Integrating quantitative trait loci (QTL) mapping and next generation sequencing can greatly help to identify variation in candidate genes, which can result in marker-assisted selection and better genetic structure information. Turbot is a commercially important flatfish in Europe and China, with available genomic information on QTLs and genome mapping. Muscle and liver RNA-seq from 18 individuals was carried out to obtain gene sequences and markers functionally related to growth, resulting in a total of 20,447 genes and 85,344 single nucleotide polymorphisms (SNPs). Many growth-related genes and SNPs were identified and placed in the turbot genome and genetic map to explore their co-localization with growth-QTL markers. Forty-five SNPs on growth-related genes were selected based on QTL co-localization and relevant function for growth traits. Forty-three SNPs were technically feasible and validated in a wild Atlantic population, where 91% were polymorphic. The integration of functional and structural genomic resources in turbot provides a practical approach for QTL mining in this species. Validated SNPs represent a useful set of growth-related gene markers for future association, functional and population studies in this flatfish species.

Substantial Downregulation of Myogenic Transcripts in Skeletal Muscle of Atlantic Cod during the Spawning Period

Gonadal maturation is an extremely energy consuming process for batch spawners and it is associated with a significant decrease in growth and seasonal deterioration in flesh quality. Our knowledge about the molecular mechanisms linking sexual maturation and muscle growth is still limited. In the present study, we performed RNA-Seq using 454 GS-FLX pyrosequencing in fast skeletal muscle sampled from two-year-old Atlantic cod (Gadus morhua) at representative time points throughout the reproductive cycle (August, March and May). In total, 126,937 good quality reads were obtained, with 546 nucleotide length and 52% GC content on average. RNA-Seq analysis using the CLC Genomics Workbench with the Atlantic cod reference UniGene cDNA data revealed 59,581 (46.9%) uniquely annotated reads. Pairwise comparison for expression levels identified 153 differentially expressed UniGenes between time points. Notably, we found a significant suppression of myh13 and myofibrillar gene isoforms in fast skeletal muscle during the spawning season. This study uncovered a large number of differentially expressed genes that may be influenced by gonadal maturation, thus representing a significant contribution to our limited understanding of the molecular mechanisms regulating muscle wasting and regeneration in batch spawners during their reproductive cycle.

De novo assembly and characterization of tissue-specific transcriptome in the endangered golden mahseer, Tor putitora

The golden mahseer (Tor putitora) graces most of the Himalayan Rivers of India and neighboring South Asian countries. Despite its several importance as a research model, as food, and in sport fishing, knowledge on transcriptome database is nil. Therefore, it was targeted to develop reference transcriptome databases of the species using next-generation sequencing. In the present study, 100,540,130 high-quality paired-end reads were obtained from six cDNA libraries of spleen, liver, gill, kidney, muscle, and brain with 28.4 GB data using Illumina paired-end sequencing technology. Tissue-specific transcriptomes as well as complete transcriptome assembly were analyzed for concise representation of the study. In brief, the de novo assembly of individual tissue resulted in an average of 31,829 (18,512–46,348) contigs per sample, while combined transcriptome comprised 77,907 unique transcript fragments (unigenes) assembled from reads of six tissues. Approximately 75,407 (96.8%) unigenes could be annotated according to their homology matches in the nr, SwisseProt, GO, or KEGG databases. Comparative analysis showed that 84% of the unigenes have significant similarity to zebra fish RefSeq proteins. Tissue-specific-dominated genes were also identified to hypothesize their localization and expression in individual tissue. In addition, 2485 simple sequence repeats (SSRs) were detected from 77,907 transcripts in the combined transcriptome of the golden mahseer. This study has generated organ-specific transcriptome profiles, which will be helpful to understand the local adaptation, genome evolution, and also future functional studies on immune system of the golden mahseer.

•

Organ specific and concatenated transcriptome of 6 organs of golden mahseer was generated using illumina sequencing.

•

A total of 100,540,130 high-quality paired-ends reads with 28.4 GB data that assembled into 77,907 contigs were obtained.

•

Contigs were annotated using GO and KEGG pathway and organ-specific few genes were identified.

•

Thirty four usable SSR markers were extracted from 2485 identified.

Transcriptomic Analysis of Metabolic Pathways in Milkfish That Respond to Salinity and Temperature Changes

Milkfish (Chanos chanos), an important marine aquaculture species in southern Taiwan, show considerable euryhalinity but have low tolerance to sudden drops in water temperatures in winter. Here, we used high throughput next-generation sequencing (NGS) to identify molecular and biological processes involved in the responses to environmental changes. Preliminary tests revealed that seawater (SW)-acclimated milkfish tolerated lower temperatures than the fresh water (FW)-acclimated group. Although FW- and SW-acclimated milkfish have different levels of tolerance for hypothermal stress, to date, the molecular physiological basis of this difference has not been elucidated. Here, we performed a next-generation sequence analysis of mRNAs from four groups of milkfish. We obtained 29669 unigenes with an average length of approximately 1936 base pairs. Gene ontology (GO) analysis was performed after gene annotation. A large number of genes for molecular regulation were identified through a transcriptomic comparison in a KEGG analysis. Basal metabolic pathways involved in hypothermal tolerance, such as glycolysis, fatty acid metabolism, amino acid catabolism and oxidative phosphorylation, were analyzed using PathVisio and Cytoscape software. Our results indicate that in response to hypothermal stress, genes for oxidative phosphorylation, e.g., succinate dehydrogenase, were more highly up-regulated in SW than FW fish. Moreover, SW and FW milkfish used different strategies when exposed to hypothermal stress: SW milkfish up-regulated oxidative phosphorylation and catabolism genes to produce more energy budget, whereas FW milkfish down-regulated genes related to basal metabolism to reduce energy loss.

RNAseq analysis of fast skeletal muscle in restriction-fed transgenic coho salmon (Oncorhynchus kisutch): an experimental model uncoupling the growth hormone and nutritional signals regulating growth

BACKGROUND

Coho salmon (Oncorhynchus kisutch) transgenic for growth hormone (Gh) express Gh in multiple tissues which results in increased appetite and continuous high growth with satiation feeding. Restricting Gh-transgenics to the same lower ration (TR) as wild-type fish (WT) results in similar growth, but with the recruitment of fewer, larger diameter, muscle skeletal fibres to reach a given body size. In order to better understand the genetic mechanisms behind these different patterns of muscle growth and to investigate how the decoupling of Gh and nutritional signals affects gene regulation we used RNA-seq to compare the fast skeletal muscle transcriptome in TR and WT coho salmon.

RESULTS

Illumina sequencing of individually barcoded libraries from 6 WT and 6 TR coho salmon yielded 704,550,985 paired end reads which were used to construct 323,115 contigs containing 19,093 unique genes of which >10,000 contained >90 % of the coding sequence. Transcripts coding for 31 genes required for myoblast fusion were identified with 22 significantly downregulated in TR relative to WT fish, including 10 (vaspa, cdh15, graf1, crk, crkl, dock1, trio, plekho1a, cdc42a and dock5) associated with signaling through the cell surface protein cadherin. Nineteen out of 44 (43 %) translation initiation factors and 14 of 47 (30 %) protein chaperones were upregulated in TR relative to WT fish.

CONCLUSIONS

TR coho salmon showed increased growth hormone transcripts and gene expression associated with protein synthesis and folding than WT fish even though net rates of protein accretion were similar. The uncoupling of Gh and amino acid signals likely results in additional costs of transcription associated with protein turnover in TR fish. The predicted reduction in the ionic costs of homeostasis in TR fish associated with increased fibre size were shown to involve multiple pathways regulating myotube fusion, particularly cadherin signaling.

Dynamics of gene expression patterns during early development of the European seabass (Dicentrarchus labrax)

Larval and embryonic stages are the most critical period in the life cycle of marine fish. Key developmental events occur early in development and are influenced by external parameters like stress, temperature, salinity, and photoperiodism. Any failure may cause malformations, developmental delays, poor growth, and massive mortalities. Advanced understanding of molecular processes underlying marine larval development may lead to superior larval rearing conditions. Today, the new sequencing and bioinformatic methods allow transcriptome screens comprising messenger (mRNA) and microRNA (miRNA) with the scope of detecting differential expression for any species of interest. In the present study, we applied Illumina technology to investigate the transcriptome of early developmental stages of the European seabass ( Dicentrarchus labrax). The European seabass, in its natural environment, is a euryhaline species and has shown high adaptation processes in early life phases. During its embryonic and larval phases the European seabass lives in a marine environment and as a juvenile it migrates to coastal zones, estuaries, and lagoons. Investigating the dynamics of gene expression in its early development may shed light on factors promoting phenotypic plasticity and may also contribute to the improvement and advancement of rearing methods of the European seabass, a species of high economic importance in European and Mediterranean aquaculture. We present the identification, characterization, and expression of mRNA and miRNA, comprising paralogous genes and differentially spliced transcripts from early developmental stages of the European seabass. We further investigated the detection of possible interactions of miRNA with mRNA.

Characterization of the transcriptome of fast and slow muscle myotomal fibres in the pacu (Piaractus mesopotamicus)

BACKGROUND

The Pacu (Piaractus mesopotamicus) is a member of the Characiform family native to the Prata Basin (South America) and a target for the aquaculture industry. A limitation for the development of a selective breeding program for this species is a lack of available genetic information. The primary objectives of the present study were 1) to increase the genetic resources available for the species, 2) to exploit the anatomical separation of myotomal fibres types to compare the transcriptomes of slow and fast muscle phenotypes and 3) to systematically investigate the expression of Ubiquitin Specific Protease (USP) family members in fast and slow muscle in response to fasting and refeeding.

RESULTS

We generated 0.6 Tb of pair-end reads from slow and fast skeletal muscle libraries. Over 665 million reads were assembled into 504,065 contigs with an average length of 1,334 bp and N50 = 2,772 bp. We successfully annotated nearly 47% of the transcriptome and identified around 15,000 unique genes and over 8000 complete coding sequences. 319 KEGG metabolic pathways were also annotated and 380 putative microsatellites were identified. 956 and 604 genes were differentially expressed between slow and fast skeletal muscle, respectively. 442 paralogues pairs arising from the teleost-specific whole genome duplication were identified, with the majority showing different expression patterns between fibres types (301 in slow and 245 in fast skeletal muscle). 45 members of the USP family were identified in the transcriptome. Transcript levels were quantified by qPCR in a separate fasting and refeeding experiment. USP genes in fast muscle showed a similar transient increase in expression with fasting as the better characterized E3 ubiquitin ligases.

CONCLUSION

We have generated a 53-fold coverage transcriptome for fast and slow myotomal muscle in the pacu (Piaractus mesopotamicus) significantly increasing the genetic resources available for this important aquaculture species. We describe significant differences in gene expression between muscle fibre types for fundamental components of general metabolism, the Pi3k/Akt/mTor network and myogenesis, including detailed analysis of paralogue expression. We also provide a comprehensive description of USP family member expression between muscle fibre types and with changing nutritional status.

Systematic variation in the pattern of gene paralog retention between the teleost superorders Ostariophysi and Acanthopterygii

Teleost fish underwent whole-genome duplication around 450 Ma followed by diploidization and loss of 80-85% of the duplicated genes. To identify a deep signature of this teleost-specific whole-genome duplication (TSGD), we searched for duplicated genes that were systematically and uniquely retained in one or other of the superorders Ostariophysi and Acanthopterygii. TSGD paralogs comprised 17-21% of total gene content. Some 2.6% (510) of TSGD paralogs were present as pairs in the Ostariophysi genomes of Danio rerio (Cypriniformes) and Astyanax mexicanus (Characiformes) but not in species from four orders of Acanthopterygii (Gasterosteiformes, Gasterosteus aculeatus; Tetraodontiformes, Tetraodon nigroviridis; Perciformes, Oreochromis niloticus; and Beloniformes, Oryzias latipes) where a single copy was identified. Similarly, 1.3% (418) of total gene number represented cases where TSGD paralogs pairs were systematically retained in the Acanthopterygian but conserved as a single copy in Ostariophysi genomes. We confirmed the generality of these results by phylogenetic and synteny analysis of 40 randomly selected linage-specific paralogs (LSPs) from each superorder and completed with the transcriptomes of three additional Ostariophysi species (Ictalurus punctatus [Siluriformes], Sinocyclocheilus species [Cypriniformes], and Piaractus mesopotamicus [Characiformes]). No chromosome bias was detected in TSGD paralog retention. Gene ontology (GO) analysis revealed significant enrichment of GO terms relative to the human GO SLIM database for "growth," "Cell differentiation," and "Embryo development" in Ostariophysi and for "Transport," "Signal Transduction," and "Vesicle mediated transport" in Acanthopterygii. The observed patterns of paralog retention are consistent with different diploidization outcomes having contributed to the evolution/diversification of each superorder.

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.

IGF-I and amino acids effects through TOR signaling on proliferation and differentiation of gilthead sea bream cultured myocytes

Skeletal muscle growth and development is controlled by nutritional (amino acids, AA) as well as hormonal factors (insulin-like growth factor, IGF-I); however, how its interaction modulates muscle mass in fish is not clearly elucidated. The purpose of this study was to analyze the development of gilthead sea bream cultured myocytes to describe the effects of AA and IGF-I on proliferating cell nuclear antigen (PCNA) and myogenic regulatory factors (MRFs) expression, as well as on the transduction pathways involved in its signaling (TOR/AKT). Our results showed that AA and IGF-I separately increased the number of PCNA-positive cells and, together produced a synergistic effect. Furthermore, AA and IGF-I, combined or separately, increased significantly Myogenin protein expression, whereas MyoD was not affected. These results indicate a role for these factors in myocyte proliferation and differentiation. At the mRNA level, AA significantly enhanced PCNA expression, but no effects were observed on the expression of the MRFs or AKT2 and FOXO3 upon treatment. Nonetheless, we demonstrated for the first time in gilthead sea bream that AA significantly increased the gene expression of TOR and its downstream effectors 4EBP1 and 70S6K, with IGF-I having a supporting role on 4EBP1 up-regulation. Moreover, AA and IGF-I also activated TOR and AKT by phosphorylation, respectively, being this activation decreased by specific inhibitors. In summary, the present study demonstrates the importance of TOR signaling on the stimulatory role of AA and IGF-I in gilthead sea bream myogenesis and contributes to better understand the potential regulation of muscle growth and development in fish.

Transcriptional assessment by microarray analysis and large-scale meta-analysis of the metabolic capacity of cardiac and skeletal muscle tissues to cope with reduced nutrient availability in Gilthead Sea Bream (Sparus aurata L.)

The effects of nutrient availability on the transcriptome of cardiac and skeletal muscle tissues were assessed in juvenile gilthead sea bream fed with a standard diet at two feeding levels: (1) full ration size and (2) 70 % satiation followed by a finishing phase at the maintenance ration. Microarray analysis evidenced a characteristic transcriptomic profile for each muscle tissue following changes in oxidative capacity (heart > red skeletal muscle > white skeletal muscle). The transcriptome of heart and secondly that of red skeletal muscle were highly responsive to nutritional changes, whereas that of glycolytic white skeletal muscle showed less ability to respond. The highly expressed and nutritionally regulated genes of heart were mainly related to signal transduction and transcriptional regulation. In contrast, those of white muscle were enriched in gene ontology (GO) terms related to proteolysis and protein ubiquitination. Microarray meta-analysis using the bioinformatic tool Fish and Chips ( http://fishandchips.genouest.org/index.php ) showed the close association of a representative cluster of white skeletal muscle with some of cardiac and red skeletal muscle, and many GO terms related to mitochondrial function appeared to be common links between them. A second round of cluster comparisons revealed that mitochondria-related GOs also linked differentially expressed genes of heart with those of liver from cortisol-treated gilthead sea bream. These results show that mitochondria are among the first responders to environmental and nutritional stress stimuli in gilthead sea bream, and functional phenotyping of this cellular organelle is highly promising to obtain reliable markers of growth performance and well-being in this fish species.

The Atlantic Bonito (Sarda sarda, Bloch 1793) transcriptome and detection of differential expression during larvae development

The Atlantic bonito (Sarda sarda, Bloch 1793) belongs to the important marine fish species with a wide geographical distribution covering the Atlantic Ocean, the Mediterranean and its bordering seas. Aquaculture practices for this species are still in their infancies and scientific studies are seldom undertaken, mainly because of difficulties in sampling. Thus for small tuna species like the Atlantic bonito only little is known about its biology and regarding the molecular background even less information is available. In the production of marine fish it is known that the most critical period is the larval stages, as high growth rates as well as significant developmental changes take place. In this study we have investigated the transcriptome of the Atlantic bonito of five larvae stages applying Illumina sequencing technology. For non-model species like aquaculture species, transcriptome analysis of RNA samples from individuals using Illumina sequencing technology is technically efficient and cost effective. In the present study a total number of 169,326,711 paired-end reads with a read length of 100 base pairs were generated resulting in a reference transcriptome of 68,220 contigs with an average length of 2054 base pairs. For differential expression analyses single end reads were obtained from different developmental stages and mapped to the constructed reference transcriptome. Differential expression analyses revealed in total 18,657 differentially expressed transcripts and were assigned to five distinguished groups. Each of the five clusters shows stage specific gene expression. We present for the first time in the Atlantic bonito an extensive RNA-Seq based characterization of its transcriptome as well as significant information on differential expression among five developmental larvae stages. The generated transcripts, including SNP and microsatellite information for candidate molecular markers and gene expression information will be a valuable resource for future genetic and molecular studies.

Differential gene expression in small and large rainbow trout derived from two seasonal spawning groups

Background

Growth in fishes is regulated via many environmental and physiological factors and is shaped by the genetic background of each individual. Previous microarray studies of salmonid growth have examined fish experiencing either muscle wastage or accelerated growth patterns following refeeding, or the influence of growth hormone and transgenesis. This study determines the gene expression profiles of genetically unmanipulated large and small fish from a domesticated salmonid strain reared on a typical feeding regime. Gene expression profiles of white muscle and liver from rainbow trout (Oncorhynchus mykiss) from two seasonal spawning groups (September and December lots) within a single strain were examined when the fish were 15 months of age to assess the influence of season (late fall vs. onset of spring) and body size (large vs. small).

Results

Although IGFBP1 gene expression was up-regulated in the livers of small fish in both seasonal lots, few expression differences were detected in the liver overall. Faster growing Dec. fish showed a greater number of differences in white muscle expression compared to Sept. fish. Significant differences in the GO Generic Level 3 categories ‘response to external stimulus’, ‘establishment of localization’, and ‘response to stress’ were detected in white muscle tissue between large and small fish. Larger fish showed up-regulation of cytoskeletal component genes while many genes related to myofibril components of muscle tissue were up-regulated in small fish. Most of the genes up-regulated in large fish within the ‘response to stress’ category are involved in immunity while in small fish most of these gene functions are related to apoptosis.

Conclusions

A higher proportion of genes in white muscle compared to liver showed similar patterns of up- or down-regulation within the same size class across seasons supporting their utility as biomarkers for growth in rainbow trout. Differences between large and small Sept. fish in the ‘response to stress’ and ‘response to external stimulus’ categories for white muscle tissue, suggests that smaller fish have a greater inability to handle stress compared to the large fish. Sampling season had a significant impact on the expression of genes related to the growth process in rainbow trout.

Transcriptome of Atlantic cod (Gadus morhua L.) early embryos from farmed and wild broodstocks

Significant efforts have been made to elucidate factors affecting egg quality in fish. Recently, we have shown that eggs originating from wild broodstock (WB) of Atlantic cod (Gadus morhua L.) are of superior quality to those derived from farmed broodstock (FB), and this is associated with differences in the chemical composition of egg yolk. However, maternal transcripts, accumulated during oogenesis, have not been studied extensively in fish. The aim of the present study was to characterize putative maternal mRNA transcriptome in fertilized eggs of Atlantic cod and to compare transcript pools between WB and FB in order to investigate the relation between egg developmental potential and putative maternal mRNA deposits. We performed high-throughput 454 pyrosequencing. For each WB and FB group, five cDNA libraries were individually tagged and sequenced, resulting in 98,687 (WB) and 119,333 (FB) average reads per library. Sequencing reads were de novo assembled, annotated, and mapped. Out of 13,726 identified isotigs, 238 were differentially expressed between WB and FB, with 155 isotigs significantly upregulated in WB. The sequence reads were mapped to 11,340 different Atlantic cod transcripts and 158 sequences were differentially expressed between the 2 groups. Important transcripts involved in fructose metabolism, fatty acid metabolism, glycerophospholipid metabolism, and oxidative phosphorylation were differentially represented between the two broodstock groups, showing potential as biomarkers of egg quality in teleosts. Our findings contribute to the hypothesis that maternal mRNAs affect egg quality and, consequently, the early development of fish.

Transient inactivation of myostatin induces muscle hypertrophy and overcompensatory growth in zebrafish via inactivation of the SMAD signaling pathway

Myostatin (MSTN) is the main negative regulator of muscle growth and development in vertebrates. In fish, little is known about the molecular mechanisms behind how MSTN inactivation triggers skeletal muscle enhancement, particularly regarding the signaling pathways involved in this process. Moreover, there have not been reports on the biotechnological applications of MSTN and its signal transduction. In this context, zebrafish underwent compensatory growth using fasting and refeeding trials, and MSTN activity was inactivated with dominant negative LAPD76A recombinant proteins during the refeeding period, when a rapid, compensatory muscle growth was observed. Treated fish displayed an overcompensation of growth characterized by higher muscle hypertrophy and growth performance than constantly fed, control fish. Treatment with LAPD76A recombinant proteins triggered inactivation of the SMAD signaling pathway in skeletal muscle, the main signal transduction used by MSTN to achieve its biological actions. Therefore, transient inactivation of MSTN during the compensatory growth of zebrafish led to a decrease in the SMAD signaling pathway in muscle, triggering muscle hypertrophy and finally improving growth performance, thus, zebrafish achieved an overcompensation of growth. The present study shows an attractive strategy for improving muscle growth in a fish species by mixing a classical strategy, such as compensatory growth, and a biotechnological approach, such as the use of recombinant proteins for inhibiting the biological actions of MSTN. The mix of both strategies may represent a method that could be applied in order to improve growth in commercial fish of interest for aquaculture.

Characterisation and expression of calpain family members in relation to nutritional status, diet composition and flesh texture in gilthead sea bream (Sparus aurata)

Calpains are non-lysosomal calcium-activated neutral proteases involved in a wide range of cellular processes including muscle proteolysis linked to post-mortem flesh softening. The aims of this study were (a) to characterise several members of the calpain system in gilthead sea bream and (b) to examine their expression in relation to nutritional status and muscle tenderisation. We identified the complete open reading frame of gilthead sea bream calpains1-3, sacapn1, sacapn2, sacapn3, and two paralogs of the calpain small subunit1, sacapns1a and sacapns1b. Proteins showed 63–90% sequence identity compared with sequences from mammals and other teleost fishes, and the characteristic domain structure of vertebrate calpains. Transcripts of sacapn1, sacapn2, sacapns1a and sacapns1b had a wide tissue distribution, whereas sacapn3 was almost exclusively detected in skeletal muscle. Next, we assessed transcript expression in skeletal muscle following alteration of nutritional status by (a) fasting and re-feeding or (b) feeding four experimental diets with different carbohydrate-to-protein ratios. Fasting significantly reduced plasma glucose and increased free fatty acids and triglycerides, together with a significant increase in sacapns1b expression. Following 7 days of re-feeding, plasma parameters returned to fed values and sacapn1, sacapn2, sacapns1a and sacapns1b expression was significantly reduced. Furthermore, an increase in dietary carbohydrate content (11 to 39%) diminished growth but increased muscle texture, which showed a significant correlation with decreased sacapn1 and sacapns1a expression, whilst the other calpains remained unaffected. This study has demonstrated that calpain expression is modulated by nutritional status and diet composition in gilthead sea bream, and that the expression of several calpain members is correlated with muscle texture, indicating their potential use as molecular markers for flesh quality in aquaculture production.

Comparative analysis of a teleost skeleton transcriptome provides insight into its regulation

An articulated endoskeleton that is calcified is a unifying innovation of the vertebrates, however the molecular basis of the structural divergence between terrestrial and aquatic vertebrates, such as teleost fish, has not been determined. In the present study long-read next generation sequencing (NGS, Roche 454 platform) was used to characterize acellular perichondral bone (vertebrae) and chondroid bone (gill arch) in the gilthead sea bream (Sparus auratus). A total of 15.97 and 14.53Mb were produced, respectively from vertebrae and gill arch cDNA libraries and yielded 32,374 and 28,371 contigs (consensus sequences) respectively. 10,455 contigs from vertebrae and 10,625 contigs from gill arches were annotated with gene ontology terms. Comparative analysis of the global transcriptome revealed 4249 unique transcripts in vertebrae, 4201 unique transcripts in the gill arches and 3700 common transcripts. Several core gene networks were conserved between the gilthead sea bream and mammalian skeleton. Transcripts for putative endocrine factors were identified in acellular gilthead sea bream bone suggesting that in common with mammalian bone it can act as an endocrine tissue. The acellular bone of the vertebra, in contrast to current opinion based on histological analysis, was responsive to a short fast and significant (p<0.05) down-regulation of several transcripts identified by NGS, osteonectin, osteocalcin, cathepsin K and IGFI occurred. In gill arches fasting caused a significant (p<0.05) down-regulation of osteocalcin and up-regulation of MMP9.

A combined strategy involving Sanger and 454 pyrosequencing increases genomic resources to aid in the management of reproduction, disease control and genetic selection in the turbot (Scophthalmus maximus)

Background

Genomic resources for plant and animal species that are under exploitation primarily for human consumption are increasingly important, among other things, for understanding physiological processes and for establishing adequate genetic selection programs. Current available techniques for high-throughput sequencing have been implemented in a number of species, including fish, to obtain a proper description of the transcriptome. The objective of this study was to generate a comprehensive transcriptomic database in turbot, a highly priced farmed fish species in Europe, with potential expansion to other areas of the world, for which there are unsolved production bottlenecks, to understand better reproductive- and immune-related functions. This information is essential to implement marker assisted selection programs useful for the turbot industry.

Results

Expressed sequence tags were generated by Sanger sequencing of cDNA libraries from different immune-related tissues after several parasitic challenges. The resulting database (“Turbot 2 database”) was enlarged with sequences generated from a 454 sequencing run of brain-hypophysis-gonadal axis-derived RNA obtained from turbot at different development stages. The assembly of Sanger and 454 sequences generated 52,427 consensus sequences (“Turbot 3 database”), of which 23,661 were successfully annotated. A total of 1,410 sequences were confirmed to be related to reproduction and key genes involved in sex differentiation and maturation were identified for the first time in turbot (AR, AMH, SRY-related genes, CYP19A, ZPGs, STAR FSHR, etc.). Similarly, 2,241 sequences were related to the immune system and several novel key immune genes were identified (BCL, TRAF, NCK, CD28 and TOLLIP, among others). The number of genes of many relevant reproduction- and immune-related pathways present in the database was 50–90% of the total gene count of each pathway. In addition, 1,237 microsatellites and 7,362 single nucleotide polymorphisms (SNPs) were also compiled. Further, 2,976 putative natural antisense transcripts (NATs) including microRNAs were also identified.

Conclusions

The combined sequencing strategies employed here significantly increased the turbot genomic resources available, including 34,400 novel sequences. The generated database contains a larger number of genes relevant for reproduction- and immune-associated studies, with an excellent coverage of most genes present in many relevant physiological pathways. This database also allowed the identification of many microsatellites and SNP markers that will be very useful for population and genome screening and a valuable aid in marker assisted selection programs.

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

Background

The gilthead sea bream (Sparus aurata) is the main fish species cultured in the Mediterranean area and constitutes an interesting model of research. Nevertheless, transcriptomic and genomic data are still scarce for this highly valuable species. A transcriptome database was constructed by de novo assembly of gilthead sea bream sequences derived from public repositories of mRNA and collections of expressed sequence tags together with new high-quality reads from five cDNA 454 normalized libraries of skeletal muscle (1), intestine (1), head kidney (2) and blood (1).

Results

Sequencing of the new 454 normalized libraries produced 2,945,914 high-quality reads and the de novo global assembly yielded 125,263 unique sequences with an average length of 727 nt. Blast analysis directed to protein and nucleotide databases annotated 63,880 sequences encoding for 21,384 gene descriptions, that were curated for redundancies and frameshifting at the homopolymer regions of open reading frames, and hosted at http://www.nutrigroup-iats.org/seabreamdb. Among the annotated gene descriptions, 16,177 were mapped in the Ingenuity Pathway Analysis (IPA) database, and 10,899 were eligible for functional analysis with a representation in 341 out of 372 IPA canonical pathways. The high representation of randomly selected stickleback transcripts by Blast search in the nucleotide gilthead sea bream database evidenced its high coverage of protein-coding transcripts.

Conclusions

The newly assembled gilthead sea bream transcriptome represents a progress in genomic resources for this species, as it probably contains more than 75% of actively transcribed genes, constituting a valuable tool to assist studies on functional genomics and future genome projects.

Sequencing of the first ayu (Plecoglossus altivelis) macrophage transcriptome and microarray development for investigation the effect of LECT2 on macrophages

Macrophages play an important role in first-line host defense of innate immune in fishes. However, it is difficult to investigate cellular mechanism of immune response in fish species with little genomic information available. Here we present the first use of RNA-Sequencing to study the macrophage transcriptome of ayu, Plecoglossus altivelis, which is an economically important fish in East Asia. De novo assembly generated 49,808 non-redundant consensus sequences, among which 23,490 transcripts found respective coding sequences. 15,707 transcripts are predicted to be involved in known metabolic or signaling pathways. The sequences were then used to develop a microarray for measurement the effect of recombinant LECT2 on ayu macrophages. LECT2 altered expression of a variety of genes mainly implicated in actin cytoskeleton, pattern recognition receptors and cytokines. Meanwhile, LECT2 enhanced phagocytosis, bacterial killing, and respiratory burst in ayu macrophages, which supported the thought derived from the microarray data that LECT2 activates macrophages. In conclusion, our results contribute to understanding the specific regulation mechanism of LECT2 in macrophage activation, and the combination of transcriptome analysis and microarray assay is a good method for screening a special tissue or cell response to a stimulus or pathogen in non-model fish species.

Development temperature has persistent effects on muscle growth responses in gilthead sea bream

Initially we characterised growth responses to altered nutritional input at the transcriptional and tissue levels in the fast skeletal muscle of juvenile gilthead sea bream. Fish reared at 21-22°C (range) were fed a commercial diet at 3% body mass d(-1) (non-satiation feeding, NSF) for 4 weeks, fasted for 4d (F) and then fed to satiation (SF) for 21d. 13 out of 34 genes investigated showed consistent patterns of regulation between nutritional states. Fasting was associated with a 20-fold increase in MAFbx, and a 5-fold increase in Six1 and WASp expression, which returned to NSF levels within 16h of SF. Refeeding to satiation was associated with a rapid (<24 h) 12 to 17-fold increase in UNC45, Hsp70 and Hsp90α transcripts coding for molecular chaperones associated with unfolded protein response pathways. The growth factors FGF6 and IGF1 increased 6.0 and 4.5-fold within 16 h and 24 h of refeeding respectively. The average growth in diameter of fast muscle fibres was checked with fasting and significant fibre hypertrophy was only observed after 13d and 21d SF. To investigate developmental plasticity in growth responses we used the same experimental protocol with fish reared at either 17.5-18.5°C (range) (LT) or 21-22°C (range) (HT) to metamorphosis and then transferred to 21-22°C. There were persistent effects of development temperature on muscle growth patterns with 20% more fibres of lower average diameter in LT than HT group of similar body size. Altering the nutritional input to the muscle to stimulate growth revealed cryptic changes in the expression of UNC45 and Hsp90α with higher transcript abundance in the LT than HT groups, whereas there were no differences in the expression of MAFbx and Six1. It was concluded that myogenesis and gene expression patterns during growth are not fixed, but can be modified by temperature during the early stages of the life cycle.