Effect of sexual maturation on muscle gene expression of rainbow trout: RNA-Seq approach

Pikeperch muscle tissues: a comparative study of structure, enzymes, genes, and proteins in wild and farmed fish

Pikeperch (Sander lucioperca) is a freshwater species and an internationally highly demanded fish in aquaculture. Despite intensive research efforts on this species, fundamental knowledge of skeletal muscle biology and structural characteristics is missing. Therefore, we conducted a comprehensive analysis of skeletal muscle parameters in adult pikeperch from two different origins, wild-caught specimens from a lake and those reared in a recirculating aquaculture system. The analyses comprised the biochemical characteristics (nucleic acid, protein content), enzyme activities (creatine kinase, lactate dehydrogenase, NADP-dependent isocitrate dehydrogenase), muscle-specific gene and protein expression (related to myofibre formation, regeneration and permanent growth, muscle structure), and muscle fibre structure. The findings reveal distinct differences between the skeletal muscle of wild and farmed pikeperch. Specifically, nucleic acid content, enzyme activity, and protein expression varied significantly. The higher enzyme activity observed in wild pikeperch suggests greater metabolically activity in their muscles. Conversely, farmed pikeperch indicated a potential for pronounced muscle growth. As the data on pikeperch skeletal muscle characteristics is sparse, the purpose of our study is to gain fundamental insights into the characteristics of adult pikeperch muscle. The presented data serve as a foundation for further research on percids' muscle biology and have the potential to contribute to advancements and adaptations in aquaculture practices.

Bioinformatics approaches and big data analytics opportunities in improving fisheries and aquaculture

Aquaculture has witnessed an excellent growth rate during the last two decades and offers huge potential to provide nutritional as well as livelihood security. Genomic research has contributed significantly toward the development of beneficial technologies for aquaculture. The existing high throughput technologies like next-generation technologies generate oceanic data which requires extensive analysis using appropriate tools. Bioinformatics is a rapidly evolving science that involves integrating gene based information and computational technology to produce new knowledge for the benefit of aquaculture. Bioinformatics provides new opportunities as well as challenges for information and data processing in new generation aquaculture. Rapid technical advancements have opened up a world of possibilities for using current genomics to improve aquaculture performance. Understanding the genes that govern economically relevant characteristics, necessitates a significant amount of additional research. The various dimensions of data sources includes next-generation DNA sequencing, protein sequencing, RNA sequencing gene expression profiles, metabolic pathways, molecular markers, and so on. Appropriate bioinformatics tools are developed to mine the biologically relevant and commercially useful results. The purpose of this scoping review is to present various arms of diverse bioinformatics tools with special emphasis on practical translation to the aquaculture industry.

Integrated Analyses of DNA Methylation and Gene Expression of Rainbow Trout Muscle under Variable Ploidy and Muscle Atrophy Conditions

Rainbow trout, Oncorhynchus mykiss, is an important cool, freshwater aquaculture species used as a model for biological research. However, its genome reference has not been annotated for epigenetic markers affecting various biological processes, including muscle growth/atrophy. Increased energetic demands during gonadogenesis/reproduction provoke muscle atrophy in rainbow trout. We described DNA methylation and its associated gene expression in atrophying muscle by comparing gravid, diploid females to sterile, triploid females. Methyl Mini-seq and RNA-Seq were simultaneously used to characterize genome-wide DNA methylation and its association with gene expression in rainbow trout muscle. Genome-wide enrichment in the number of CpGs, accompanied by depleted methylation levels, was noticed around the gene transcription start site (TSS). Hypermethylation of CpG sites within ±1 kb on both sides of TSS (promoter and gene body) was weakly/moderately associated with reduced gene expression. Conversely, hypermethylation of the CpG sites in downstream regions of the gene body +2 to +10 kb was weakly associated with increased gene expression. Unlike mammalian genomes, rainbow trout gene promotors are poor in CpG islands, at <1% compared to 60%. No signs of genome-wide, differentially methylated (DM) CpGs were observed due to the polyploidy effect; only 1206 CpGs (0.03%) were differentially methylated, and these were primarily associated with muscle atrophy. Twenty-eight genes exhibited differential gene expression consistent with methylation levels of 31 DM CpGs. These 31 DM CpGs represent potential epigenetic markers of muscle atrophy in rainbow trout. The DM CpG-harboring genes are involved in apoptosis, epigenetic regulation, autophagy, collagen metabolism, cell membrane functions, and Homeobox proteins. Our study also identified genes explaining higher water content and modulated glycolysis previously shown as characteristic biochemical signs of rainbow trout muscle atrophy associated with sexual maturation. This study characterized DNA methylation in the rainbow trout genome and its correlation with gene expression. This work also identified novel epigenetic markers associated with muscle atrophy in fish/lower vertebrates.

Degradation of Muscle Quality in Hybrid Grouper (♀ Epinephelus fuscoguttatus × ♂ Epinephelus lanceolatu) Due to Oxidative Damage Caused by Ingestion of Oxidized Fish Oil

The aim of the study was to investigate the effects of fresh fish oil (FFO) and oxidized fish oil (OFO) diets on the muscle quality of hybrid grouper (♀ Epinephelus fuscoguttatus × ♂ E. lanceolatu). Hybrid grouper were fed with diets containing 9% FFO or OFO for 60 days. Muscle sample were collected at 0, 30, and 60 days and the selected indexes of muscle were measured. Malondialdehyde (MDA), reactive oxygen species (ROS), total cholesterol (TC) and triglycerides (TG) in grouper muscle accumulated gradually with prolonged ingestion time, especially OFO group. Total saturated fatty acids (ΣSAFA) was significantly reduced and total polyunsaturated fatty acids (ΣPUFA) was significantly increased of muscle in FFO group; meanwhile, the muscle ΣSAFA and monounsaturated fatty acids (ΣMUFA) contents in the OFO group were significantly higher than those in the FFO group and the ΣPUFA (especially C22:5n3, C22:6n3) contents was significantly lower than that in the FFO group at 60 days. Consumption of OFO diet for 60 days reduced the diversity of volatile compounds, significantly reduced the content of total esters and increased the content of total aldehydes and total aromatics in grouper muscle. Furthermore, ingestion of OFO diet significantly reduced the mRNA expression of fraction growth factors and antioxidant genes in the muscle of grouper. In conclusion, the increasing MDA content in FO and the oxidative rancidity of PUFA can cause the deterioration of grouper quality and flavor due to oxidative muscle damage.

Gene expression profiling of trout muscle during flesh quality recovery following spawning

Background

Sexual maturation causes loss of fish muscle mass and deterioration of fillet quality attributes that prevent market success. We recently showed that fillet yield and flesh quality recover in female trout after spawning. To gain insight into the molecular mechanisms regulating flesh quality recovery, we used an Agilent-based microarray platform to conduct a large-scale time course analysis of gene expression in female trout white muscle from spawning to 33 weeks post-spawning.

Results

In sharp contrast to the situation at spawning, muscle transcriptome of female trout at 33 weeks after spawning was highly similar to that of female trout of the same cohort that did not spawn, which is consistent with the post-spawning flesh quality recovery. Large-scale time course analysis of gene expression in trout muscle during flesh quality recovery following spawning led to the identification of approximately 3340 unique differentially expressed genes that segregated into four major clusters with distinct temporal expression profiles and functional categories. The first cluster contained approximately 1350 genes with high expression at spawning and downregulation after spawning and was enriched with genes linked to mitochondrial ATP synthesis, fatty acid catabolism and proteolysis. A second cluster of approximately 540 genes with transient upregulation 2 to 8 weeks after spawning was enriched with genes involved in transcription, RNA processing, translation, ribosome biogenesis and protein folding. A third cluster containing approximately 300 genes upregulated 4 to 13 weeks after spawning was enriched with genes encoding ribosomal subunits or regulating protein folding. Finally, a fourth cluster that contained approximately 940 genes with upregulation 8 to 24 weeks after spawning, was dominated by genes encoding myofibrillar proteins and extracellular matrix components and genes involved in glycolysis.

Conclusion

Overall, our study indicates that white muscle tissue restoration and flesh quality recovery after spawning are associated with transcriptional changes promoting anaerobic ATP production, muscle fibre hypertrophic growth and extracellular matrix remodelling. The generation of the first database of genes associated with post-spawning muscle recovery may provide insights into the molecular and cellular mechanisms controlling muscle yield and fillet quality in fish and provide a useful list of potential genetic markers for these traits.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-021-08228-3.

Integrated metabolomic and gene expression analyses to study the effects of glycerol monolaurate on flesh quality in large yellow croaker (Larimichthys crocea)

To improve the quality of cultured large yellow croaker (Larimichthys crocea), this study was performed to study the impacts of glycerol monolaurate (GML) on the nutritional value, growth performance, muscle texture, and taste intensity of L. crocea. The results showed that GML as a feed additive significantly increased the crude lipid content and reduced the diameters of muscle fibers, which in turn markedly altered the flesh texture in terms of cohesiveness. Moreover, the taste indicators (umami and richness) and flavor-related amino acid (glutamic acid, glycine, and proline) contents of L. crocea muscle were significantly higher in the GML group. Metabolomic and gene expression analyses showed that GML supplementation could significantly improve amino acid biosynthesis and metabolism, promote protein and lipid synthesis, and activate myogenic-related signaling pathways of L. crocea. Consequently, adding an appropriate amount of GML to fish feed would be conducive to providing healthy, nutrient-rich and acceptably flavored aquatic-products.

Genome-wide scan for common variants associated with intramuscular fat and moisture content in rainbow trout

BACKGROUND

Genetic improvement of fillet quality attributes is a priority of the aquaculture industry. Muscle composition impacts quality attributes such as flavor, appearance, texture, and juiciness. Fat and moisture make up about ~ 80% of the tissue weight. The genetic architecture underlying the fat and moisture content of the muscle is still to be fully explored in fish. A 50 K gene transcribed SNP chip was used for genotyping 789 fish with available phenotypic data for fat and moisture content. Genotyped fish were obtained from two consecutive generations produced in the National Center for Cool and Cold Water Aquaculture (NCCCWA) growth-selective breeding program. Estimates of SNP effects from weighted single-step GBLUP (WssGBLUP) were used to perform genome-wide association (GWA) analysis to identify quantitative trait loci (QTL) associated with the studied traits.

RESULTS

Using genomic sliding windows of 50 adjacent SNPs, 137 and 178 SNPs were identified as associated with fat and moisture content, respectively. Chromosomes 19 and 29 harbored the highest number of SNPs explaining at least 2% of the genetic variation in fat and moisture content. A total of 61 common SNPs on chromosomes 19 and 29 affected the aforementioned traits; this association suggests common mechanisms underlying intramuscular fat and moisture content. Additionally, based on single-marker GWA analyses, 8 and 24 SNPs were identified in association with fat and moisture content, respectively.

CONCLUSION

SNP-harboring genes were primarily involved in lipid metabolism, cytoskeleton remodeling, and protein turnover. This work provides putative SNP markers that could be prioritized and used for genomic selection in breeding programs.

Modification of the content of n-3 highly unsaturated fatty acid, chemical composition, and lipid nutritional indices in the meat of grass carp (Ctenopharyngodon idella) fed alfalfa (Medicago sativa) pellets

A 120-d feeding trial was conducted to determine the effect of alfalfa (Medicago sativa) feeding on growth and chemical composition, fatty acid content, and nutritional and lipid indices of the meat of grass carp (Ctenopharyngodon idella). Two experimental diets were used: alfalfa pellet (AP) diet and artificial grain diet (GD). Final weight, feed conversion rate, and protein efficiency ratio were significantly greater in the GD group (P < 0.05). However, no differences in the length and condition factor were observed. The composition of the meat differed between treatments. The protein content was significantly greater in the AP group (P < 0.05), while the lipid and cholesterol contents were significantly greater in the GD group (P < 0.05). A greater proportion of saturated, n-6 polyunsaturated, and n-6 highly unsaturated fatty acids was obtained in the GD group. The AP group accumulated a greater concentration of eicosapentaenoic (EPA), docosapentaenoic (DPA), and docosahexaenoic (DHA) acids (P < 0.05). The fatty acid composition of the meat determined a significant decrease in the thrombogenicity index and saturation index (S/P) in the AP group (P < 0.05). The Elongase index was greater in the GD group (P < 0.05). In contrast, the AP group had a greater index of Δ9 Desaturase and Δ5 + Δ6 Desaturase for n-3 and n-6 fatty acids (P < 0.05). These results suggest that alfalfa feeding decreases the growth of C. idella but improves the quality of meat by increasing the protein, EPA, and DHA contents. It also reduces cholesterol content and improves nutritional indices.

Characterization and localization of cyclin B3 transcript in both oocyte and spermatocyte of the rainbow trout (Oncorhynchus mykiss)

B-type cyclins are regulatory subunits with distinct roles in the cell cycle. To date, at least three subtypes of B-type cyclins (B1, B2, and B3) have been identified in vertebrates. Previously, we reported the characterization and expression profiles of cyclin B1 and B2 during gametogenesis in the rainbow trout (Oncorhynchus mykiss). In this paper, we isolated another subtype of cyclin B, cyclin B3 (CB3), from a cDNA library of the rainbow trout oocyte. The full-length CB3 cDNA (2,093 bp) has an open reading frame (1,248 bp) that encodes a protein of 416 amino acid residues. The CB3 transcript was widely distributed in all the examined tissues, namely, eye, gill, spleen, brain, heart, kidney, stomach, skin, muscle, and, especially, gonad. Northern blot analysis indicated only one form of the CB3 transcript in the testis and ovary. In situ hybridization revealed that, in contrast to cyclin B1 and B2 transcripts, CB3 transcripts were localized in the oocytes, spermatocytes, and spermatogonia. These findings strongly suggest that CB3 plays a role not only as a mitotic cyclin in spermatogonial proliferation during early spermatogenesis but also during meiotic maturation of the spermatocyte and oocyte in the rainbow trout.

Genome-Wide Association Study Identifies Genomic Loci Affecting Filet Firmness and Protein Content in Rainbow Trout

Filet quality traits determine consumer satisfaction and affect profitability of the aquaculture industry. Soft flesh is a criterion for fish filet downgrades, resulting in loss of value. Filet firmness is influenced by many factors, including rate of protein turnover. A 50K transcribed gene SNP chip was used to genotype 789 rainbow trout, from two consecutive generations, produced in the USDA/NCCCWA selective breeding program. Weighted single-step GBLUP (WssGBLUP) was used to perform genome-wide association (GWA) analyses to identify quantitative trait loci affecting filet firmness and protein content. Applying genomic sliding windows of 50 adjacent SNPs, 212 and 225 SNPs were associated with genetic variation in filet shear force and protein content, respectively. Four common SNPs in the ryanodine receptor 3 gene (RYR3) affected the aforementioned filet traits; this association suggests common mechanisms underlying filet shear force and protein content. Genes harboring SNPs were mostly involved in calcium homeostasis, proteolytic activities, transcriptional regulation, chromatin remodeling, and apoptotic processes. RYR3 harbored the highest number of SNPs (n = 32) affecting genetic variation in shear force (2.29%) and protein content (4.97%). Additionally, based on single-marker analysis, a SNP in RYR3 ranked at the top of all SNPs associated with variation in shear force. Our data suggest a role for RYR3 in muscle firmness that may be considered for genomic- and marker-assisted selection in breeding programs of rainbow trout.

Genome-Wide Association Analysis With a 50K Transcribed Gene SNP-Chip Identifies QTL Affecting Muscle Yield in Rainbow Trout

Detection of coding/functional SNPs that change the biological function of a gene may lead to identification of putative causative alleles within QTL regions and discovery of genetic markers with large effects on phenotypes. This study has two-fold objectives, first to develop, and validate a 50K transcribed gene SNP-chip using RNA-Seq data. To achieve this objective, two bioinformatics pipelines, GATK and SAMtools, were used to identify ~21K transcribed SNPs with allelic imbalances associated with important aquaculture production traits including body weight, muscle yield, muscle fat content, shear force, and whiteness in addition to resistance/susceptibility to bacterial cold-water disease (BCWD). SNPs ere identified from pooled RNA-Seq data collected from ~620 fish, representing 98 families from growth- and 54 families from BCWD-selected lines with divergent phenotypes. In addition, ~29K transcribed SNPs without allelic-imbalances were strategically added to build a 50K Affymetrix SNP-chip. SNPs selected included two SNPs per gene from 14K genes and ~5K non-synonymous SNPs. The SNP-chip was used to genotype 1728 fish. The average SNP calling-rate for samples passing quality control (QC; 1,641 fish) was ≥ 98.5%. The second objective of this study was to test the feasibility of using the new SNP-chip in GWA (Genome-wide association) analysis to identify QTL explaining muscle yield variance. GWA study on 878 fish (representing 197 families from 2 consecutive generations) with muscle yield phenotypes and genotyped for 35K polymorphic markers (passing QC) identified several QTL regions explaining together up to 28.40% of the additive genetic variance for muscle yield in this rainbow trout population. The most significant QTLs were on chromosomes 14 and 16 with 12.71 and 10.49% of the genetic variance, respectively. Many of the annotated genes in the QTL regions were previously reported as important regulators of muscle development and cell signaling. No major QTLs were identified in a previous GWA study using a 57K genomic SNP chip on the same fish population. These results indicate improved detection power of the transcribed gene SNP-chip in the target trait and population, allowing identification of large-effect QTLs for important traits in rainbow trout.

Muscle transcriptome resource for growth, lipid metabolism and immune system in Hilsa shad, Tenualosa ilisha

The information on the genes involved in muscle growth, lipid metabolism and immune systems would help to understand the mechanisms during the spawning migration in Hilsa shad, which in turn would be useful in its future domestication process. The primary objective of this study was to generate the transcriptome profile of its muscle through RNA seq. The total RNA was isolated and library was prepared from muscle tissue of Tenualosa ilisha, which was collected from Padma River at Farakka, India. The prepared library was then sequenced by Illumina HiSeq platform, HiSeq 2000, using paired-end strategy. A total of 8.68 GB of pair-end reads of muscle transcriptome was generated, and 43,384,267 pair-end reads were assembled into 3,04,233 contigs, of which 23.99% of assembled contigs has length ≥ 150 bp. The total GO terms were categorised into cellular component, molecular function and biological process through PANTHER database. Fifty-three genes related to muscle growth were identified and genes in different pathways were: 75 in PI3/AKT, 46 in mTOR, 76 in MAPK signalling, 24 in Janus kinase-signal transducer and activator of transcription, 45 in AMPK and 27 in cGMP pathways. This study also mined the genes involved in lipid metabolism, in which glycerophospholipid metabolism contained highest number of genes (32) and four were found to be involved in fatty acid biosynthesis. There were 58 immune related genes found, in which 31 were under innate and 27 under adaptive immunity. The present study included a large genomic resource of T. ilisha muscle generated through RNAseq, which revealed the essential dataset for our understanding of regulatory processes, specifically during the seasonal spawning migration. As Hilsa is a slow growing fish, the genes identified for muscle growth provided the basic information to study myogenesis. In addition, genes identified for lipid metabolism and immune system would provide resources for lipid synthesis and understanding of Hilsa defense mechanisms, respectively.

Integrated analysis of lncRNA and mRNA expression in rainbow trout families showing variation in muscle growth and fillet quality traits

Muscle yield and quality traits are important for the aquaculture industry and consumers. Genetic selection for these traits is difficult because they are polygenic and result from multifactorial interactions. To study the genetic architecture of these traits, phenotypic characterization of whole body weight (WBW), muscle yield, fat content, shear force and whiteness were measured in ~500 fish representing 98 families from a growth-selected line. RNA-Seq was used to sequence the muscle transcriptome of different families exhibiting divergent phenotypes for each trait. We have identified 240 and 1,280 differentially expressed (DE) protein-coding genes and long noncoding RNAs (lncRNAs), respectively, in fish families exhibiting contrasting phenotypes. Expression of many DE lncRNAs (n = 229) was positively correlated with overlapping, neighboring or distantly located protein-coding genes (n = 1,030), resulting in 3,392 interactions. Three DE antisense lncRNAs were co-expressed with sense genes known to impact muscle quality traits. Forty-four DE lncRNAs had potential sponge functions to miRNAs that affect muscle quality traits. This study (1) defines muscle quality associated protein-coding and noncoding genes and (2) provides insight into non-coding RNAs involvement in regulating growth and fillet quality traits in rainbow trout.

Crosstalk among lncRNAs, microRNAs and mRNAs in the muscle 'degradome' of rainbow trout

In fish, protein-coding and noncoding genes involved in muscle atrophy are not fully characterized. In this study, we characterized coding and noncoding genes involved in gonadogenesis-associated muscle atrophy, and investigated the potential functional interplay between these genes. Using RNA-Seq, we compared expression pattern of mRNAs, long noncoding RNAs (lncRNAs) and microRNAs of atrophying skeletal muscle from gravid females and control skeletal muscle from age-matched sterile individuals. A total of 852 mRNAs, 1,160 lncRNAs and 28 microRNAs were differentially expressed (DE) between the two groups. Muscle atrophy appears to be mediated by many genes encoding ubiquitin-proteasome system, autophagy related proteases, lysosomal proteases and transcription factors. Transcripts encoding atrogin-1 and mir-29 showed exceptional high expression in atrophying muscle, suggesting an important role in bulk muscle proteolysis. DE genes were co-localized in the genome with strong expression correlation, and they exhibited extensive 'lncRNA-mRNA', 'lncRNA-microRNA', 'mRNA-microRNA' and 'lncRNA-protein' physical interactions. DE genes exhibiting potential functional interactions comprised the highly correlated 'lncRNA-mRNA-microRNA' gene network described as 'degradome'. This study pinpoints extensive coding and noncoding RNA interactions during muscle atrophy in fish, and provides valuable resources for future mechanistic studies.

Diet Affects Muscle Quality and Growth Traits of Grass Carp (Ctenopharyngodon idellus): A Comparison Between Grass and Artificial Feed

Fish muscle, the main edible parts with high protein level and low fat level, is consumed worldwide. Diet contributes greatly to fish growth performance and muscle quality. In order to elucidate the correlation between diet and muscle quality, the same batch of juvenile grass carp (Ctenopharyngodon idellus) were divided into two groups and fed with either grass (Lolium perenne, Euphrasia pectinata and Sorghum sudanense) or artificial feed, respectively. However, the different two diets didn't result in significant differences in all the detected water quality parameters (e.g., Tm, pH, DO, NH3/[Formula: see text]-N, [Formula: see text]-N, [Formula: see text], TN, TP, and TOC) between the two experimental groups. After a 4-month culture period, various indexes and expression of myogenic regulatory factor (MRFs) and their related genes were tested. The weight gain of the fish fed with artificial feed (AFG) was nearly 40% higher than the fish fed with grass (GFG). Significantly higher alkaline phosphatase, total cholestrol, high density cholestrol and total protein were detected in GFG as compared to AFG. GFG also showed increased hardness, resilience and shear force in texture profile analysis, with significantly bigger and compact muscle fibers in histologic slices. The fat accumulation was most serious in the abdomen muscle of AFG. Additionally, the expression levels of MyoG, MyoD, IGF-1, and MSTNs were higher, whereas Myf-5, MRF4, and IGF-2 were lower in most positional muscles of GFG as compared to AFG. Overall, these results suggested that feeding grass could promote muscle growth and development by stimulating muscle fiber hypertrophy, as well as significantly enhance the expression of CoL1As. Feeding C. idellus with grass could also improve flesh quality by improving muscle characteristics, enhancing the production of collagen, meanthile, reducing fat accumulation and moisture in muscle, but at the cost of a slower growth.

Identification of SNPs associated with muscle yield and quality traits using allelic-imbalance analyses of pooled RNA-Seq samples in rainbow trout

Background

Coding/functional SNPs change the biological function of a gene and, therefore, could serve as “large-effect” genetic markers. In this study, we used two bioinformatics pipelines, GATK and SAMtools, for discovering coding/functional SNPs with allelic-imbalances associated with total body weight, muscle yield, muscle fat content, shear force, and whiteness. Phenotypic data were collected for approximately 500 fish, representing 98 families (5 fish/family), from a growth-selected line, and the muscle transcriptome was sequenced from 22 families with divergent phenotypes (4 low- versus 4 high-ranked families per trait).

Results

GATK detected 59,112 putative SNPs; of these SNPs, 4798 showed allelic imbalances (>2.0 as an amplification and <0.5 as loss of heterozygosity). SAMtools detected 87,066 putative SNPs; and of them, 4962 had allelic imbalances between the low- and high-ranked families. Only 1829 SNPs with allelic imbalances were common between the two datasets, indicating significant differences in algorithms. The two datasets contained 7930 non-redundant SNPs of which 4439 mapped to 1498 protein-coding genes (with 6.4% non-synonymous SNPs) and 684 mapped to 295 lncRNAs. Validation of a subset of 92 SNPs revealed 1) 86.7–93.8% success rate in calling polymorphic SNPs and 2) 95.4% consistent matching between DNA and cDNA genotypes indicating a high rate of identifying SNPs with allelic imbalances. In addition, 4.64% SNPs revealed random monoallelic expression. Genome distribution of the SNPs with allelic imbalances exhibited high density for all five traits in several chromosomes, especially chromosome 9, 20 and 28. Most of the SNP-harboring genes were assigned to important growth-related metabolic pathways.

Conclusion

These results demonstrate utility of RNA-Seq in assessing phenotype-associated allelic imbalances in pooled RNA-Seq samples. The SNPs identified in this study were included in a new SNP-Chip design (available from Affymetrix) for genomic and genetic analyses in rainbow trout.

Electronic supplementary material

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

Triploid Atlantic salmon shows similar performance, fatty acid composition and proteome response to diploids during early freshwater rearing

There is currently renewed interest in farming triploid Atlantic salmon. Improving farming requires identifying triploid specific phenotypic and physiological traits that are uniquely derived from ploidy per se and developed under optimal growing conditions. This study investigated firstly, the impact of ploidy on growth performance and whole body composition of Atlantic salmon at different early freshwater stages [34dph (days post-hatching) alevin, 109dph fry, and 162dph parr] and secondly, whether phenotypic differences at these stages were reflected in protein samples collected from whole fish, white muscle or liver tissue. Female diploid and triploid Atlantic salmon (n=3) were first fed at 35dph and then maintained by feeding to satiation on commercial feeds. Triploids were significantly lower in weight at the late alevin and fry stages but matched diploid weight at the parr stage. The whole-body lipid content was significantly higher for triploids at the parr stage, while the whole-body lipid class profile was broadly similar and was largely not affected by ploidy. Comparative label-free shotgun proteomic analysis did not detect significant alterations in protein expression between diploids and triploids at any growth stage. The present results indicate that ploidy under optimal growing conditions and during early freshwater stages only result in small phenotypic differences in weight and whole body lipid content that were not reflected at the proteome level. These findings suggest that optimal husbandry conditions for freshwater Atlantic salmon are similar between ploidies, at least for all-female populations.

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.

Evaluation of potential candidate genes involved in salinity tolerance in striped catfish (Pangasianodon hypophthalmus) using an RNA-Seq approach

Increasing salinity levels in freshwater and coastal environments caused by sea level rise linked to climate change is now recognized to be a major factor that can impact fish growth negatively, especially for freshwater teleost species. Striped catfish (Pangasianodon hypophthalmus) is an important freshwater teleost that is now widely farmed across the Mekong River Delta in Vietnam. Understanding the basis for tolerance and adaptation to raised environmental salinity conditions can assist the regional culture industry to mitigate predicted impacts of climate change across this region. Attempt of next generation sequencing using the ion proton platform results in more than 174 million raw reads from three tissue libraries (gill, kidney and intestine). Reads were filtered and de novo assembled using a variety of assemblers and then clustered together to generate a combined reference transcriptome. Downstream analysis resulted in a final reference transcriptome that contained 60,585 transcripts with an N50 of 683 bp. This resource was further annotated using a variety of bioinformatics databases, followed by differential gene expression analysis that resulted in 3062 transcripts that were differentially expressed in catfish samples raised under two experimental conditions (0 and 15 ppt). A number of transcripts with a potential role in salinity tolerance were then classified into six different functional gene categories based on their gene ontology assignments. These included; energy metabolism, ion transportation, detoxification, signal transduction, structural organization and detoxification. Finally, we combined the data on functional salinity tolerance genes into a hypothetical schematic model that attempted to describe potential relationships and interactions among target genes to explain the molecular pathways that control adaptive salinity responses in P. hypophthalmus. Our results indicate that P. hypophthalmus exhibit predictable plastic regulatory responses to elevated salinity by means of characteristic gene expression patterns, providing numerous candidate genes for future investigations.

mRNA-seq reveals skeletal muscle atrophy in response to handling stress in a marine teleost, the red cusk-eel (Genypterus chilensis)

Background

Fish reared under intensive conditions are repeatedly exposed to stress, which negatively impacts growth. Although most fish follow a conserved pattern of stress response, with increased concentrations of cortisol, each species presents specificities in the cell response and stress tolerance. Therefore, culturing new species requires a detailed knowledge of these specific responses. The red cusk-eel (Genypterus chilensis) is a new economically important marine species for the Chilean aquaculture industry. However, there is no information on the stress- and cortisol-induced mechanisms that decrease skeletal muscle growth in this teleost.

Results

Using Illumina RNA-seq technology, skeletal muscle sequence reads for G. chilensis were generated under control and handling stress conditions. Reads were mapped onto a reference transcriptome, resulting in the in silico identification of 785 up-regulated and 167 down-regulated transcripts. Gene ontology enrichment analysis revealed a significant up-regulation of catabolic genes associated with skeletal muscle atrophy. These results were validated by RT-qPCR analysis for ten candidates genes involved in ubiquitin-mediated proteolysis, autophagy and skeletal muscle growth. Additionally, using a primary culture of fish skeletal muscle cells, the effect of cortisol was evaluated in relation to red cusk-eel skeletal muscle atrophy.

Conclusions

The present data demonstrated that handling stress promotes skeletal muscle atrophy in the marine teleost G. chilensis through the expression of components of the ubiquitin-proteasome and autophagy-lysosome systems. Furthermore, cortisol was a powerful inductor of skeletal muscle atrophy in fish myotubes. This study is an important step towards understanding the atrophy system in non-model teleost species and provides novel insights on the cellular and molecular mechanisms that control skeletal muscle growth in early vertebrates.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-015-2232-7) contains supplementary material, which is available to authorized users.