Transcriptome profiling in fast versus slow-growing rainbow trout across seasonal gradients

Mustn1 in Skeletal Muscle: A Novel Regulator?

Skeletal muscle is a complex organ essential for locomotion, posture, and metabolic health. This review explores our current knowledge of Mustn1, particularly in the development and function of skeletal muscle. Mustn1 expression originates from Pax7-positive satellite cells in skeletal muscle, peaks during around the third postnatal month, and is crucial for muscle fiber differentiation, fusion, growth, and regeneration. Clinically, Mustn1 expression is potentially linked to muscle-wasting conditions such as muscular dystrophies. Studies have illustrated that Mustn1 responds dynamically to injury and exercise. Notably, ablation of Mustn1 in skeletal muscle affects a broad spectrum of physiological aspects, including glucose metabolism, grip strength, gait, peak contractile strength, and myofiber composition. This review summarizes our current knowledge of Mustn1's role in skeletal muscle and proposes future research directions, with a goal of elucidating the molecular function of this regulatory gene.

Comparative transcriptomic analyses of brain-liver-muscle in channel catfish (Ictalurus punctatus) with differential growth rate

Growth is an important economically trait for aquatic animals. The popularity of farmed channel catfish (Ictalurus punctatus) in China has recently surged, prompting a need for research into the genetic mechanisms that drive growth and development to expedite the selection of fast-growing variants. In this study, the brain, liver and muscle transcriptomes of channel catfish between fast-growing and slow-growing groups were analyzed using RNA-Seq. Totally, 63, 110 and 86 differentially expressed genes (DEGs) were from brain, liver and muscle tissues. DEGs are primarily involved in growth, development, metabolism and immunity, which are related to the growth regulation of channel catfish, such as growth hormone receptor b (ghrb), fibroblast growth factor receptor 4 (fgfr4), bone morphogenetic protein 1a (bmp1a), insulin-like growth factor 2a (igf2a), collagen, type I, alpha 1a (col1a1a), acyl-CoA synthetase long chain family member 2 (acsl2) and caveolin 1 (cav1). This study advances our knowledge of the genetic mechanisms accounting for differences in growth rate and offers crucial gene resources for future growth-related molecular breeding programs in channel catfish.

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.

Mustn1 ablation in skeletal muscle results in functional alterations

Mustn1, a gene expressed exclusively in the musculoskeletal system, was shown in previous in vitro studies to be a key regulator of myogenic differentiation and myofusion. Other studies also showed Mustn1 expression associated with skeletal muscle development and hypertrophy. However, its specific role in skeletal muscle function remains unclear. This study sought to investigate the effects of Mustn1 in a conditional knockout (KO) mouse model in Pax7 positive skeletal muscle satellite cells. Specifically, we investigated the potential effects of Mustn1 on myogenic gene expression, grip strength, alterations in gait, ex vivo investigations of isolated skeletal muscle isometric contractions, and potential changes in the composition of muscle fiber types. Results indicate that Mustn1 KO mice did not present any substantial phenotypic changes or significant variations in genes related to myogenic differentiation and fusion. However, an approximately 10% decrease in overall grip strength was observed in the 2-month-old KO mice in comparison to the control wild type (WT), but this decrease was not significant when normalized by weight. KO mice also generated approximately 8% higher vertical force than WT at 4 months in the hindlimb. Ex vivo experiments revealed decreases in about 20 to 50% in skeletal muscle contractions and about 10%-20% fatigue in soleus of both 2- and 4-month-old KO mice, respectively. Lastly, immunofluorescent analyses showed a persistent increase of Type IIb fibers up to 15-fold in the KO mice while Type I fibers decreased about 20% and 30% at both 2 and 4 months, respectively. These findings suggest a potential adaptive or compensatory mechanism following Mustn1 loss, as well as hinting at an association between Mustn1 and muscle fiber typing. Collectively, Mustn1's complex roles in skeletal muscle physiology requires further research, particularly in terms of understanding the potential role of Mustn1 in muscle repair and regeneration, as well as with influence of exercise. Collectively, these will offer valuable insights into Mustn1's key biological functions and regulatory pathways.

Identification of differentially expressed genes and SNPs linked to harvest body weight of genetically improved rohu carp, Labeo rohita

In most of the aquaculture selection programs, harvest body weight has been a preferred performance trait for improvement. Molecular interplay of genes linked to higher body weight is not elucidated in major carp species. The genetically improved rohu carp with 18% average genetic gain per generation with respect to harvest body weight is a promising candidate for studying genes' underlying performance traits. In the present study, muscle transcriptome sequencing of two groups of individuals, with significant difference in breeding value, belonging to the tenth generation of rohu carp was performed using the Illumina HiSeq 2000 platform. A total of 178 million paired-end raw reads were generated to give rise to 173 million reads after quality control and trimming. The genome-guided transcriptome assembly and differential gene expression produced 11,86,119 transcripts and 451 upregulated and 181 downregulated differentially expressed genes (DEGs) between high-breeding value and low-breeding value (HB & LB) groups, respectively. Similarly, 39,158 high-quality coding SNPs were identified with the Ts/Tv ratio of 1.23. Out of a total of 17 qPCR-validated transcripts, eight were associated with cellular growth and proliferation and harbored 13 SNPs. The gene expression pattern was observed to be positively correlated with RNA-seq data for genes such as myogenic factor 6, titin isoform X11, IGF-1 like, acetyl-CoA, and thyroid receptor hormone beta. A total of 26 miRNA target interactions were also identified to be associated with significant DETs (p-value < 0.05). Genes such as Myo6, IGF-1-like, and acetyl-CoA linked to higher harvest body weight may serve as candidate genes in marker-assisted breeding and SNP array construction for genome-wide association studies and genomic selection.

Mustn1 ablation in skeletal muscle results in increased glucose tolerance concomitant with upregulated GLUT expression in male mice

Glucose homeostasis is closely regulated to maintain energy requirements of vital organs and skeletal muscle plays a crucial role in this process. Mustn1 is expressed during embryonic and postnatal skeletal muscle development and its function has been implicated in myogenic differentiation and myofusion. Whether Mustn1 plays a role in glucose homeostasis in anyway remains largely unknown. As such, we deleted Mustn1 in skeletal muscle using a conditional knockout (KO) mouse approach. KO mice did not reveal any specific gross phenotypic alterations in skeletal muscle. However, intraperitoneal glucose tolerance testing (IPGTT) revealed that 2-month-old male KO mice had significantly lower glycemia than their littermate wild type (WT) controls. These findings coincided with mRNA changes in genes known to be involved in glucose metabolism, tolerance, and insulin sensitivity; 2-month-old male KO mice had significantly higher expression of GLUT1 and GLUT10 transporters, MUP-1 while OSTN expression was lower. These differences in glycemia and gene expression were statistically insignificant after 4 months. Identical experiments in female KO and WT control mice did not indicate any differences at any age. Our results suggest a link between Mustn1 expression and glucose homeostasis during a restricted period of skeletal muscle development/maturation. While this is an observational study, Mustn1's relationship to glucose homeostasis appears to be more complex with a possible connection to other key proteins such as GLUTs, MUP-1, and OSTN. Additionally, our data indicate temporal and sex differences. Lastly, our findings strengthen the notion that Mustn1 plays a role in the metabolic capacity of skeletal muscle.

Aquaculture omics: An update on the current status of research and data analysis

Aquaculture is the fast-growing agricultural sector and has the ability to meet the growing demand for protein nutritional security for future population. In future aquaculture is going to be the major source of fish proteins as capture fisheries reached at its maximum. However, several challenges need to overcome such as lack of genetically improved strains/varieties, lack of species-specific feed/functional feed, round the year availability of quality fish seed, pollution of ecosystems and increased frequencies of disease occurrence etc. In recent years, the continuous development of high throughput sequencing technology has revolutionized the biological sciences and provided necessary tools. Application of 'omics' in aquaculture research have been successfully used to resolve several productive and reproductive issues and thus ensure its sustainability and profitability. To date, high quality draft genomes of over fifty fish species have been generated and successfully used to develop large number of single nucleotide polymorphism markers (SNPs), marker panels and other genomic resources etc in several aquaculture species. Similarly, transcriptome profiling and miRNAs analysis have been used in aquaculture research to identify key transcripts and expression analysis of candidate genes/miRNAs involved in reproduction, immunity, growth, development, stress toxicology and disease. Metagenome analysis emerged as a promising scientific tool to analyze the complex genomes contained within microbial communities. Metagenomics has been successfully used in the aquaculture sector to identify novel and potential pathogens, antibiotic resistance genes, microbial roles in microcosms, microbial communities forming biofloc, probiotics etc. In the current review, we discussed application of high-throughput technologies (NGS) in the aquaculture sector.

The first transcriptome sequencing and data analysis of the Javan mahseer (Tor tambra)

The Javan mahseer (Tor tambra) is one of the most valuable freshwater fish found in Tor species. To date, other than mitogenomic data (BioProject: PRJNA422829), genomic and transcriptomic resources for this species are still lacking which is crucial to understand the molecular mechanisms associated with important traits such as growth, immune response, reproduction and sex determination. For the first time, we sequenced the transcriptome from a whole juvenile fish using Illumina NovaSEQ6000 generating raw paired-end reads. De novo transcriptome assembly generated a draft transcriptome (BUSCO5 completeness of 91.2% [Actinopterygii_odb10 database]) consisting of 259,403 putative transcripts with a total and N50 length of 333,881,215 bp and 2283 bp, respectively. A total count of 77,503 non-redundant protein coding sequences were predicted from the transcripts and used for functional annotation. We mapped the predicted proteins to 304 known KEGG pathways with signal transduction cluster having the highest representation followed by immune system and endocrine system. In addition, transcripts exhibiting significant similarity to previously published growth-and immune-related genes were identified which will facilitate future molecular breeding of Tor tambra.

Revelation of candidate genes and molecular mechanism of reproductive seasonality in female rohu (Labeo rohita Ham.) by RNA sequencing

BACKGROUND

Carp fish, rohu (Labeo rohita Ham.) is important freshwater aquaculture species of South-East Asia having seasonal reproductive rhythm. There is no holistic study at transcriptome level revealing key candidate genes involved in such circannual rhythm regulated by biological clock genes (BCGs). Seasonality manifestation has two contrasting phases of reproduction, i.e., post-spawning resting and initiation of gonadal activity appropriate for revealing the associated candidate genes. It can be deciphered by RNA sequencing of tissues involved in BPGL (Brain-Pituitary-Gonad-Liver) axis controlling seasonality. How far such BCGs of this fish are evolutionarily conserved across different phyla is unknown. Such study can be of further use to enhance fish productivity as seasonality restricts seed production beyond monsoon season.

RESULT

A total of ~ 150 Gb of transcriptomic data of four tissues viz., BPGL were generated using Illumina TruSeq. De-novo assembled BPGL tissues revealed 75,554 differentially expressed transcripts, 115,534 SSRs, 65,584 SNPs, 514 pathways, 5379 transcription factors, 187 mature miRNA which regulates candidate genes represented by 1576 differentially expressed transcripts are available in the form of web-genomic resources. Findings were validated by qPCR. This is the first report in carp fish having 32 BCGs, found widely conserved in fish, amphibian, reptile, birds, prototheria, marsupials and placental mammals. This is due to universal mechanism of rhythmicity in response to environment and earth rotation having adaptive and reproductive significance.

CONCLUSION

This study elucidates evolutionary conserved mechanism of photo-periodism sensing, neuroendocrine secretion, metabolism and yolk synthesis in liver, gonadal maturation, muscular growth with sensory and auditory perception in this fish. Study reveals fish as a good model for research on biological clock besides its relevance in reproductive efficiency enhancement.

Comparative Transcriptome Analysis Revealed Genes Involved in Sexual and Polyploid Growth Dimorphisms in Loach (Misgurnus anguillicaudatus)

Simple Summary

Misgurnus anguillicaudatus not only exhibits sexual size dimorphism, but also shows polyploid size dimorphism. Here, we performed comparative transcriptome integration analysis of multiple tissues of diploid and tetraploid M. anguillicaudatus of both sexes. We found that differences in energy metabolism and steroid hormone synthesis levels may be the main causes of sexual and polyploidy growth dimorphisms of M. anguillicaudatus. Fast-growing M. anguillicaudatus (tetraploids, females) have higher levels of energy metabolism and lower steroid hormone synthesis and fatty acid degradation abilities than slow-growing M. anguillicaudatus (diploids, males).

Abstract

Sexual and polyploidy size dimorphisms are widespread phenomena in fish, but the molecular mechanisms remain unclear. Loach (Misgurnus anguillicaudatus) displays both sexual and polyploid growth dimorphism phenomena, and are therefore ideal models to study these two phenomena. In this study, RNA-seq was used for the first time to explore the differentially expressed genes (DEGs) between both sexes of diploid and tetraploid loaches in four tissues (brain, gonad, liver, and muscle). Results showed that 21,003, 17, and 1 DEGs were identified in gonad, liver, and muscle tissues, respectively, between females and males in both diploids and tetraploids. Regarding the ploidy levels, 4956, 1496, 2187, and 1726 DEGs were identified in the brain, gonad, liver, and muscle tissues, respectively, between tetraploids and diploids of the same sex. When both sexual and polyploid size dimorphisms were considered simultaneously in the four tissues, only 424 DEGs were found in the gonads, indicating that these gonadal DEGs may play an important regulatory role in regulating sexual and polyploid size dimorphisms. Regardless of the sex or ploidy comparison, the significant DEGs involved in glycolysis/gluconeogenesis and oxidative phosphorylation pathways were upregulated in faster-growing individuals, while steroid hormone biosynthesis-related genes and fatty acid degradation and elongation-related genes were downregulated. This suggests that fast-growing loaches (tetraploids, females) have higher energy metabolism levels and lower steroid hormone synthesis and fatty acid degradation abilities than slow-growing loaches (diploids, males). Our findings provide an archive for future systematic research on fish sexual and polyploid dimorphisms.

Cross-species Exon Capture and Whole Exome Sequencing: Application, Utility and Challenges for Genomic Resource Development in Non-model Species

Comprehending the genetic architecture of complex traits has many applications in evolution, ecology, conservation biology and plant and animal production systems. Underlying research questions in these fields are diverse species that often have limited genetic information available. In aquaculture, for example, genetic progress has been slow in many species due to a lack in such genetic information. In this study, zebrafish (as a well-studied model species) was used in cross-species transfer to develop genomic resources and identify candidate genes underling growth differentials in dusky kob. Dusky kob is a Sciaenid finfish and an emerging aquaculture species. The zebrafish All Exon Predesigned Probe-set capture protocol was used to enrich fractionated DNA samples from kob, classified as either large or small, before massive parallel sequencing on the Ion Torrent platform. Although vast quantities of sequence data were generated, only about 30% of contigs could be identified as zebrafish homologues. There were numerous species-specific sequences and inconsistent coverage of sequencing products across samples, likely due to non-specific binding of the probe-set as a result of the evolutionary divergence between zebrafish and kob. Nonetheless, more than 55,000 SNPs could be reliably identified and genotyped to the individual level. Using SNP genotypic divergence estimates, between large and small cohorts, a number of candidate genes associated with growth was also identified for future investigation. These findings contribute to the growing body of evidence demonstrating the utility of a cross-species capture approach in the development of important genomic resources for understanding traits of interest in species without reference genomes.

Differences in DNA methylation between slow and fast muscle in Takifugu rubripes

Fish skeletal muscle is comprised of fast muscle (FM) and slow muscle (SM), which constitutes 60% of total the body mass. Fish skeletal muscle can affect fish swimming activity, which is important for aquaculture due to its growth-potentiating effects. DNA methylation can influence gene expression level. We previously identified multiple differentially expressed genes (DEGs) between FM and SM in Takifugu rubripes. However, it is unknown if the expression levels of these DEGs are influenced by DNA methylation. In the present study, we used DNA methylation sequencing to study the DNA methylation profiles of FM and SM in T. rubripes. SM had higher overall methylation levels than FM. A total of 8479 differentially methylated genes (DMGs) and 3407 DMGs containing differentially methylated regions (DMRs) in the promoter regions between FM and SM were identified. After enrichment analysis, we found functionally relevant DMGs between FM and SM, including Kapca, Plcd3a, Plcd1, Pi3k, Tsp4b and Pgfrb in the hedgehog signaling pathway and phosphatidylinositol (PI)-related pathways. Due to the different methylation levels of these genes between FM and SM, the expression levels of Kapca, Plcd3a, Plcd1, Pi3k, and Tsp4b were higher in FM and Pgfrb was higher in SM. There were differences in the hedgehog signaling pathway and PI-related pathways between FM and SM. In SM, the cytokine-cytokine receptor interaction promoted focal adhesion, while ECM-receptor interactions promoted focal adhesion in FM. These results provide information regarding the difference between FM and SM in T. rubripes.

The Effect of Continuous Light on Growth and Muscle-Specific Gene Expression in Atlantic Salmon (Salmo salar L.) Yearlings

Photoperiod is associated to phenotypic plasticity of somatic growth in several teleost species, however, the molecular mechanisms underlying this phenomenon are currently unknown. The effect of a continuous lighting (LD 24:0), compared with the usual hatchery lighting (HL) regime, on the growth rate and gene expression of myogenic regulatory factors (MRFs: MyoD1 paralogs, Myf5, and MyoG) myosin heavy chain (MyHC), and MSTN paralogs in the white muscles of hatchery-reared Atlantic salmon yearlings was evaluated over a 6-month period (May–October). The levels of gene expression were determined using real-time PCR. Continuous lighting was shown to have a positive effect on weight gain. MyHC, MyoD1c, MyoD1b, and MSTN1a/b mRNA expression was influenced by the light regime applied. In all the studied groups, a significant positive correlation was observed between the expression levels of MRFs and MSTN paralogs throughout the experiment. The study demonstrated seasonal patterns regarding the simultaneous expression of several MRFs. MyoD1a, MyoG, and MyHC mRNA expression levels were elevated in the mid-October, but MyoD1b/c, and Myf5 mRNA levels decreased by the end of this month. In general, the findings showed that constant lighting affected the regulatory mechanisms of muscle growth processes in salmon.

Revealing Alteration in the Hepatic Glucose Metabolism of Genetically Improved Carp, Jayanti Rohu Labeo rohita Fed a High Carbohydrate Diet Using Transcriptome Sequencing

Although feed cost is the greatest concern in aquaculture, the inclusion of carbohydrates in the fish diet, and their assimilation, are still not well understood in aquaculture species. We identified molecular events that occur due to the inclusion of high carbohydrate levels in the diets of genetically improved 'Jayanti rohu' Labeo rohita. To reveal transcriptional changes in the liver of rohu, a feeding experiment was conducted with three doses of gelatinized starch (20% (control), 40%, and 60%). Transcriptome sequencing revealed totals of 15,232 (4464 up- and 4343 down-regulated) and 15,360 (4478 up- and 4171 down-regulated) differentially expressed genes. Up-regulated transcripts associated with glucose metabolisms, such as hexokinase, PHK, glycogen synthase and PGK, were found in fish fed diets with high starch levels. Interestingly, a de novo lipogenesis mechanism was found to be enriched in the livers of treated fish due to up-regulated transcripts such as FAS, ACCα, and PPARγ. The insulin signaling pathways with enriched PPAR and mTOR were identified by Kyoto Encyclopedia of Genes and Genome (KEGG) as a result of high carbohydrates. This work revealed for the first time the atypical regulation transcripts associated with glucose metabolism and lipogenesis in the livers of Jayanti rohu due to the inclusion of high carbohydrate levels in the diet. This study also encourages the exploration of early nutritional programming for enhancing glucose efficiency in carp species, for sustainable and cost-effective aquaculture production.

Transcriptomic Response to Selective Breeding for Fast Growth in Rainbow Trout (Oncorhynchus mykiss)

Genetic improvement for faster growth is a conventional approach to increase growth rates in aquaculture species; however, the genetic and physiological factors regulating growth performance in fish are not fully characterized. The objective of this study was to identify physiological mechanisms associated with faster growth rates by comparing the liver and muscle transcriptome of a rainbow trout line selectively bred for fast growth (growth line, GL) and a contemporary randomly mated control line (synthetic control, SC) from the same selective breeding program. A third genetic line from a commercial egg supplier (commercial A, CA) was also included to characterize differences in gene expression profiles between populations. Body weight of the GL at harvest was approximately 20% and 8% heavier (p < 0.05) than SC and CA, respectively. There were 145 and 36 differentially expressed genes (DEG) in liver and white muscle, respectively, between the GL and SC that were enriched for the growth hormone/insulin-like growth factor axis (GH/IGF) and PI3K-Akt, JAK-STAT, MAPK, and cAMP signal transduction pathways. A greater concentration of plasma IGF-I was detected in the GL compared with SC (p < 0.05). A unique gene profile was detected in CA, with 11 and 210 DEG in liver and white muscle; these genes associated with innate immunity, complement systems, and metabolic pathways. Collectively, these findings provide a more extensive characterization of the fast-growth phenotype in fish that furthers knowledge of the physiological basis for genetic variation in growth performance in selectively bred rainbow trout.

Comparative Analysis of the Transcriptome and Distribution of Putative SNPs in Two Rainbow Trout (Oncorhynchus mykiss) Breeding Strains by Using Next-Generation Sequencing

Selective breeding can significantly improve the establishment of sustainable and profitable aquaculture fish farming. For rainbow trout (Oncorhynchus mykiss), one of the main aquaculture coldwater species in Europe, a variety of selected hatchery strains are commercially available. In this study, we investigated the genetic variation between the local Born strain, selected for survival, and the commercially available Silver Steelhead strain, selected for growth. We sequenced the transcriptome of six tissues (gills, head kidney, heart, liver, spleen, and white muscle) from eight healthy individuals per strain, using RNA-seq technology to identify strain-specific gene-expression patterns and single nucleotide polymorphisms (SNPs). In total, 1760 annotated genes were differentially expressed across all tissues. Pathway analysis assigned them to different gene networks. We also identified a set of SNPs, which are heterozygous for one of the two breeding strains: 1229 of which represent polymorphisms over all tissues and individuals. Our data indicate a strong genetic differentiation between Born and Silver Steelhead trout, despite the relatively short time of evolutionary separation of the two breeding strains. The results most likely reflect their specifically adapted genotypes and might contribute to the understanding of differences regarding their robustness toward high stress and pathogenic challenge described in former studies.

Transcriptome Analysis Identified Genes for Growth and Omega-3/-6 Ratio in Saline Tilapia

Growth and omega-3/-6 ratio are important traits in aquaculture. The mechanisms underlying quick growth and high omega-3/-6 ratio in fish are not fully understood. The consumption of the meat of tilapia suffers a bad reputation due to its low omega-3/-6 ratio. To facilitate the improvement of these traits and to understand more about the mechanisms underlying quick growth and high omega-3/-6 ratio, we conducted transcriptome analysis in the muscle and liver of fast- and slow-growing hybrid saline tilapia generated by crossing Mozambique tilapia and red tilapia. A transcriptome with an average length of 963 bp was generated by using 486.65 million clean 100 bp paired-end reads. A total of 42,699 annotated unique sequences with an average length of 3.4 kb were obtained. Differentially expressed genes (DEGs) in the muscle and liver were identified between fast- and slow-growing tilapia. Pathway analysis classified these genes into many pathways. Ten genes, including foxK1, sparc, smad3, usp38, crot, fadps, sqlea, cyp7b1, impa1, and gss, from the DEGs were located within QTL for growth and omega-3, which were previously detected content in tilapia, suggesting that these ten genes could be important candidate genes for growth and omega-3 fatty acid content. Analysis of SNPs in introns 1 and 2 of foxK1 revealed that the SNPs were significantly associated with growth and omega-3/-6 ratio. This study lays the groundwork for further investigation of the molecular mechanisms underlying the phenotypic variation of these two traits and provides SNPs for selecting these traits at fingerling stage.

Genome-wide association analysis for body weight identifies candidate genes related to development and metabolism in rainbow trout (Oncorhynchus mykiss)

Growth is one of the most important traits from both a physiological and economic perspective in aquaculture species. Thus, identifying the genomic regions and genes underpinning genetic variation for this trait is of particular interest in several fish species, including rainbow trout. In this work, we perform a genome-wide association study (GWAS) to identify the genomic regions associated with body weight at tagging (BWT) and at 18 months (BW18M) using a dense SNP panel (57 k) and 4596 genotyped rainbow trout from 105 full-sib families belonging to a Chilean breeding population. Analysis was performed by means of single-step GBLUP approach. Genetic variance explained by 20 adjacent SNP windows across the whole genome is reported. To further explore candidate genes, we focused on windows that explained the highest proportion of genetic variance in the top 10 chromosomes for each trait. The main window from the top 10 chromosomes was explored by BLAST using the first and last SNP position of each window to determine the target nucleotide sequence. As expected, the percentage of genetic variance explained by windows was relatively low, due to the polygenic nature of body weight. The most important genomic region for BWT and BW18M were located on chromosomes 15 and 24 and they explained 2.14% and 3.02% of the genetic variance for each trait, respectively. Candidate genes including several growth factors, genes involved in development of skeletal muscle and bone tissue and nutrient metabolism were identified within the associated regions for both traits BWT and BW18M. These results indicate that body weight is polygenic in nature in rainbow trout, with the most important loci explaining as much as 3% of the genetic variance for the trait. The genes identified here represent good candidates for further functional validation to uncover biological mechanisms underlying variation for growth in rainbow trout.

Transcriptomic profiling reveals molecular regulation of seasonal reproduction in Tibetan highland fish, Gymnocypris przewalskii

BACKGROUND

The Tibetan highland fish, Gymnocypris przewalskii, migrates from Lake Qinghai to its spawning grounds every summer. This seasonal reproduction is critically regulated by intrinsic and extrinsic signals. However, the molecular mechanisms that process environmental oscillations to initiate the seasonal mating are largely unknown.

RESULTS

A transcriptomic analysis was conducted on the brain and gonad of male and female G. przewalskii in reproductive and nonreproductive seasons. We obtained 2034, 760, 1158 and 17,856 differentially expressed genes between the reproductively active and dormant female brain, male brain, ovary and testis. Among these genes, DIO2 was upregulated in the reproductively active brain and gonad of both males and females. Neuroactive ligand-receptor genes were activated in male and female brain. Functional enrichment analysis suggested that retinol metabolism was uniquely stimulated in reproductively active males. Genes involved in GnRH signaling and sex hormone synthesis exhibited higher expression levels in brain and gonad during the reproductive season. A co-expression network classified all the genes into 9 modules. The network pinpointed CDC42 as the hub gene that connected the pathways in responsible for modulating reproduction in G. przewalskii. Meanwhile, the sex pheromone receptor gene prostaglandin receptor was identified to link to multiple endocrine receptors, such as GnRHR2 in the network.

CONCLUSIONS

The current study profiled transcriptomic variations between reproductively active and dormant fish, highlighting the potential regulatory mechanisms of seasonal reproduction in G. przewalskii. Our data suggested that the seasonal regulation of reproduction in G. przewalskii was controlled by the external stimulation of photoperiodic variations. The activated transcription of neuroendocrine and sex hormone synthesis genes contributed to seasonal reproduction regulation in G. przewalskii, which was presumably influenced by the increased day-length during the breeding season.

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.

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.

Co-expression Analysis of Sirtuins and Related Metabolic Biomarkers in Juveniles of Gilthead Sea Bream (Sparus aurata) With Differences in Growth Performance

Sirtuins (SIRTs) represent a conserved protein family of deacetylases that act as master regulators of metabolism, but little is known about their roles in fish and livestock animals in general. The present study aimed to assess the value of SIRTs for the metabolic phenotyping of fish by assessing their co-expression with a wide-representation of markers of energy and lipid metabolism and intestinal function and health in two genetically different gilthead sea bream strains with differences in growth performance. Fish from the fast-growing strain exhibited higher feed intake, feed efficiency and plasma IGF-I levels, along with higher hepatosomatic index and lower mesenteric fat (lean phenotype). These observations suggest differences in tissue energy partitioning with an increased flux of fatty acids from adipose tissue toward the liver. The resulting increased risk of hepatic steatosis may be counteracted in the liver by reduced lipogenesis and enhanced triglyceride catabolism, in combination with a higher and more efficient oxidative metabolism in white skeletal muscle. These effects were supported by co-regulated changes in the expression profile of SIRTs (liver, sirt1; skeletal muscle, sirt2; adipose tissue, sirt5-6) and markers of oxidative metabolism (pgc1α, cpt1a, cs, nd2, cox1), mitochondrial respiration uncoupling (ucp3) and fatty acid and triglyceride metabolism (pparα, pparγ, elovl5, scd1a, lpl, atgl) that were specific to each strain and tissue. The anterior intestine of the fast-growing strain was better suited to cope with improved growth by increased expression of markers of nutrient absorption (fabp2), epithelial barrier integrity (cdh1, cdh17) and immunity (il1β, cd8b, lgals1, lgals8, sIgT, mIgT), which were correlated with low expression levels of sirt4 and markers of fatty acid oxidation (cpt1a). In the posterior intestine, the fast-growing strain showed a consistent up-regulation of sirt2, sirt3, sirt5 and sirt7 concurrently with increased expression levels of markers of cell proliferation (pcna), oxidative metabolism (nd2) and immunity (sIgT, mIgT). Together, these findings indicate that SIRTs may play different roles in the regulation of metabolism, inflammatory tone and growth in farmed fish, arising as powerful biomarkers for a reliable metabolic phenotyping of fish at the tissue-specific level.

Differential gene expression during early development in recently evolved and sympatric Arctic charr morphs

Phenotypic differences between closely related taxa or populations can arise through genetic variation or be environmentally induced, leading to altered transcription of genes during development. Comparative developmental studies of closely related species or variable populations within species can help to elucidate the molecular mechanisms related to evolutionary divergence and speciation. Studies of Arctic charr (Salvelinus alpinus) and related salmonids have revealed considerable phenotypic variation among populations and in Arctic charr many cases of extensive variation within lakes (resource polymorphism) have been recorded. One example is the four Arctic charr morphs in the ∼10,000 year old Lake Thingvallavatn, which differ in numerous morphological and life history traits. We set out to investigate the molecular and developmental roots of this polymorphism by studying gene expression in embryos of three of the morphs reared in a common garden set-up. We performed RNA-sequencing, de-novo transcriptome assembly and compared gene expression among morphs during an important timeframe in early development, i.e., preceding the formation of key trophic structures. Expectedly, developmental time was the predominant explanatory variable. As the data were affected by some form of RNA-degradation even though all samples passed quality control testing, an estimate of 3'-bias was the second most common explanatory variable. Importantly, morph, both as an independent variable and as interaction with developmental time, affected the expression of numerous transcripts. Transcripts with morph effect, separated the three morphs at the expression level, with the two benthic morphs being more similar. However, Gene Ontology analyses did not reveal clear functional enrichment of transcripts between groups. Verification via qPCR confirmed differential expression of several genes between the morphs, including regulatory genes such as AT-Rich Interaction Domain 4A (arid4a) and translin (tsn). The data are consistent with a scenario where genetic divergence has contributed to differential expression of multiple genes and systems during early development of these sympatric Arctic charr morphs.

Somatotropic Axis Regulation Unravels the Differential Effects of Nutritional and Environmental Factors in Growth Performance of Marine Farmed Fishes

The Gh/Prl/Sl family has evolved differentially through evolution, resulting in varying relationships between the somatotropic axis and growth rates within and across fish species. This is due to a wide range of endogenous and exogenous factors that make this association variable throughout season and life cycle, and the present minireview aims to better define the nutritional and environmental regulation of the endocrine growth cascade over precisely defined groups of fishes, focusing on Mediterranean farmed fishes. As a result, circulating Gh and Igf-i are revitalized as reliable growth markers, with a close association with growth rates of gilthead sea bream juveniles with deficiency signs in both macro- or micro-nutrients. This, together with other regulated responses, promotes the use of Gh and Igf-i as key performance indicators of growth, aerobic scope, and nutritional condition in gilthead sea bream. Moreover, the sirtuin-energy sensors might modulate the growth-promoting action of somatotropic axis. In this scenario, transcripts of igf-i and gh receptors mirror changes in plasma Gh and Igf-i levels, with the ghr-i/ghr-ii expression ratio mostly unaltered over season. However, this ratio is nutritionally regulated, and enriched plant-based diets or diets with specific nutrient deficiencies downregulate hepatic ghr-i, decreasing the ghr-i/ghr-ii ratio. The same trend, due to a ghr-ii increase, is found in skeletal muscle, whereas impaired growth during overwintering is related to increase in the ghr-i/ghr-ii and igf-ii/igf-i ratios in liver and skeletal muscle, respectively. Overall, expression of insulin receptors and igf receptors is less regulated, though the expression quotient is especially high in the liver and muscle of sea bream. Nutritional and environmental regulation of the full Igf binding protein 1-6 repertoire remains to be understood. However, tissue-specific expression profiling highlights an enhanced and nutritionally regulated expression of the igfbp-1/-2/-4 clade in liver, whereas the igfbp-3/-5/-6 clade is overexpressed and regulated in skeletal muscle. The somatotropic axis is, therefore, highly informative of a wide-range of growth-disturbing and stressful stimuli, and multivariate analysis supports its use as a reliable toolset for the assessment of growth potentiality and nutrient deficiencies and requirements, especially in combination with selected panels of other nutritionally regulated metabolic biomarkers.

G0S2a1 (G0/G1 switch gene 2a1) is downregulated by TNF-α in grass carp (Ctenopharyngodon idellus) hepatocytes through PPARα inhibition

G0/G1 switch gene 2 plays an important role in the regulation of lipolysis in mammals, but little is known about its gene (G0S2) structure and function in fish. In the present study, two genes, G0S2a and G0S2b were isolated and characterized from grass carp Ctenopharyngodon idella, which encode peptides of 111 and 84 amino acids, respectively. Moreover, alternative multiple exon usage resulted in a significant variation in the 5'-region of G0S2a transcripts yielding two isoforms (G0S2a1 and G0S2a2). Phylogenetic and synteny analyses indicated that G0S2a and G0S2b could have originated from the teleost-specific genome duplication event. Analysis of the exon-intron structures clarified that G0S2a contained an extra intron compared with G0S2b. G0S2a1, G0S2a2 and G0S2b mRNAs were highly expressed in adipose tissue and liver. G0S2a was localized to the cytoplasm and nucleus, while G0S2b was mainly localized in cytoplasm, suggesting that G0S2a and G0S2b may have different functions in grass carp. PPARα agonist caused an increase in G0S2a1 and G0S2b expression, revealing that they are subject to transcriptional control by PPARα-mediated signals. TNF-α treatment decreased G0S2a1 and G0S2a2 transcripts that paralleled TNF-α downregulation of PPARα; however, only the effects of TNF-α on G0S2a1 were attenuated by treatment with PPARα agonist. Our findings identify G0S2a, not G0S2b, as a target gene for TNF-α and reveal that TNF-α suppresses G0S2a1 gene expression through a PPARα-dependent pathway in grass carp hepatocytes.

Chronic stress inhibits growth and induces proteolytic mechanisms through two different nonoverlapping pathways in the skeletal muscle of a teleost fish

Chronic stress detrimentally affects animal health and homeostasis, with somatic growth, and thus skeletal muscle, being particularly affected. A detailed understanding of the underlying endocrine and molecular mechanisms of how chronic stress affects skeletal muscle growth remains lacking. To address this issue, the present study assessed primary (plasma cortisol), secondary (key components of the GH/IGF system, muscular proteolytic pathways, and apoptosis), and tertiary (growth performance) stress responses in fine flounder ( Paralichthys adspersus) exposed to crowding chronic stress. Levels of plasma cortisol, glucocorticoid receptor 2 ( gr2), and its target genes ( klf15 and redd1) mRNA increased significantly only at 4 wk of crowding ( P < 0.05). The components of the GH/IGF system, including ligands, receptors, and their signaling pathways, were significantly downregulated at 7 wk of crowding ( P < 0.05). Interestingly, chronic stress upregulated the ubiquitin-proteasome pathway and the intrinsic apoptosis pathways at 4wk ( P < 0.01), whereas autophagy was only significantly activated at 7 wk ( P < 0.05), and meanwhile the ubiquitin-proteasome and the apoptosis pathways returned to control levels. Overall growth was inhibited in fish in the 7-wk chronic stress trial ( P < 0.05). In conclusion, chronic stress directly affects muscle growth and downregulates the GH/IGF system, an action through which muscular catabolic mechanisms are promoted by two different and nonoverlapping proteolytic pathways. These findings provide new information on molecular mechanisms involved in the negative effects that chronic stress has on muscle anabolic/catabolic signaling balance.

Differential gene expression and SNP association between fast- and slow-growing turbot (Scophthalmus maximus)

Growth is among the most important traits for animal breeding. Understanding the mechanisms underlying growth differences between individuals can contribute to improving growth rates through more efficient breeding schemes. Here, we report a transcriptomic study in muscle and brain of fast- and slow-growing turbot (Scophthalmus maximus), a relevant flatfish in European and Asian aquaculture. Gene expression and allelic association between the two groups were explored. Up-regulation of the anaerobic glycolytic pathway in the muscle of fast-growing fish was observed, indicating a higher metabolic rate of white muscle. Brain expression differences were smaller and not associated with major growth-related genes, but with regulation of feeding-related sensory pathways. Further, SNP variants showing frequency differences between fast- and slow-growing fish pointed to genomic regions likely involved in growth regulation, and three of them were individually validated through SNP typing. Although different mechanisms appear to explain growth differences among families, general mechanisms seem also to be involved, and thus, results provide a set of useful candidate genes and markers to be evaluated for more efficient growth breeding programs and to perform comparative genomic studies of growth in fish and vertebrates.

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.

Sex-specific differences in transcriptome profiles of brain and muscle tissue of the tropical gar

Background

The tropical gar (Atractosteus tropicus) is the southernmost species of the seven extant species of gar fishes in the world. In Mexico and Central America, the species is an important food source due to its nutritional quality and low price. Despite its regional importance and increasing concerns about overexploitation and habitat degradation, basic genetic information on the tropical gar is lacking. Determining genetic information on the tropical gar is important for the sustainable management of wild populations, implementation of best practices in aquaculture settings, evolutionary studies of ancient lineages, and an understanding of sex-specific gene expression. In this study, the transcriptome of the tropical gar was sequenced and assembled de novo using tissues from three males and three females using Illumina sequencing technology. Sex-specific and highly differentially expressed transcripts in brain and muscle tissues between adult males and females were subsequently identified.

Results

The transcriptome was assembled de novo resulting in 80,611 transcripts with a contig N50 of 3,355 base pairs and over 168 kilobases in total length. Male muscle, brain, and gonad as well as female muscle and brain were included in the assembly. The assembled transcriptome was annotated to identify the putative function of expressed transcripts using Trinotate and SwissProt, a database of well-annotated proteins. The brain and muscle datasets were then aligned to the assembled transcriptome to identify transcripts that were differentially expressed between males and females. The contrast between male and female brain identified 109 transcripts from 106 genes that were significantly differentially expressed. In the muscle comparison, 82 transcripts from 80 genes were identified with evidence for significant differential expression. Almost all genes identified as differentially expressed were sex-specific. The differentially expressed transcripts were enriched for genes involved in cellular functioning, signaling, immune response, and tissue-specific functions.

Conclusions

This study identified differentially expressed transcripts between male and female gar in muscle and brain tissue. The majority of differentially expressed transcripts had sex-specific expression. Expanding on these findings to other developmental stages, populations, and species may lead to the identification of genetic factors contributing to the skewed sex ratio seen in the tropical gar and of sex-specific differences in expression in other species. Finally, the transcriptome assembly will open future research avenues on tropical gar development, cell function, environmental resistance, and evolution in the context of other early vertebrates.

Electronic supplementary material

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

Transcriptome assembly and identification of genes and SNPs associated with growth traits in largemouth bass (Micropterus salmoides)

Growth is one of the most crucial economic traits of all aquaculture species, but the molecular mechanisms involved in growth of largemouth bass (Micropterus salmoides) are poorly understood. The objective of this study was to screen growth-related genes of M. salmoides by RNA sequencing and identify growth-related single-nucleotide polymorphism (SNP) markers through a growth association study. The muscle transcriptomes of fast- and slow-growing largemouth bass were obtained using the RNA-Seq technique. A total of 54,058,178 and 54,742,444 qualified Illumina read pairs were obtained for the fast-growing and slow-growing groups, respectively, giving rise to 4,865,236,020 and 4,926,819,960 total clean bases, respectively. Gene expression profiling showed that 3,530 unigenes were differentially expressed between the fast-growing and slow-growing phenotypes (false discovery rate ≤0.001, the absolute value of log₂ (fold change) ≥1), including 1,441 up-regulated and 2,889 down-regulated unigenes in the fast-growing largemouth bass. Analysis of these genes revealed that several signalling pathways, including the growth hormone-insulin-like growth factor 1 axis and signalling pathway, the glycolysis pathway, and the myostatin/transforming growth factor beta signalling pathway, as well as heat shock protein, cytoskeleton, and myofibril component genes might be associated with muscle growth. From these genes, 10 genes with putative SNPs were selected, and 17 SNPs were genotyped successfully. Marker-trait analysis in 340 individuals of Youlu No. 1 largemouth bass revealed three SNPs associated with growth in key genes (phosphoenolpyruvate carboxykinase 1, FOXO3b, and heat shock protein beta-1). This research provides information about key genes and SNPs related to growth, providing new clues to understanding the molecular basis of largemouth bass growth.

Identification of differentially expressed genes associated with differential body size in mandarin fish (Siniperca chuatsi)

Body size is an obvious and important characteristic of fish. Mandarin fish Siniperca chuatsi (Basilewsky) is one of the most valuable perciform species widely cultured in China. Individual differences in body size are common in mandarin fish and significantly influence the aquaculture production. However, little is currently known about its genetic control. In this study, digital gene expression profiling and transcriptome sequencing were performed in mandarin fish with differential body size at 30 and 180 days post-hatch (dph), respectively. Body weight, total length and body length of fish with big-size were significantly higher than those with small-size at both 30 and 180 dph (P < 0.05). 2171 and 2014 differentially expressed genes were identified between small-size and big-size fish at 30 and 180 dph, respectively. RT quantitative PCR (qPCR) analysis showed that the differential expression of 10 selected genes in mandarin fish that went through the same training procedure. The genes were involved in the growth hormone-insulin-like growth factor axis, cell proliferation and differentiation, appetite control, glucose metabolism, reproduction and sexual size dimorphism pathways. This study will help toward a comprehensive understanding of the complexity of regulation of body size in mandarin fish individuals and provide valuable information for future research.