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

Multi-organ transcriptomics analysis of a slowly growing fish rock carp (Procypris rabaudi) reveals insights into mechanism of growth rate regulation

To explore the patterns of differentially expressed genes (DEGs) associated with different growth rates in rock carp (Procypris rabaudi), transcriptome sequencing was performed on the muscle, liver, and brain tissues of rock carp. Subsequently, bioinformatics analysis was conducted, and 2129, 1380, and 415 DEGs were identified in the muscle, liver, and brain tissues, respectively. GO enrichment and KEGG pathway analysis revealed that genes related to appetite regulation, protein degradation and digestion, lipid transport and metabolisms, and glycolysis/gluconeogenesis were upregulated in individuals with slower growth rates. Differential expression analysis identified 21 genes associated with feeding and metabolism across three tissues, including mc4r, npy, and npry in brain tissue; fatp, fabp, pparα, and apo in liver tissue; and prss, ctrl, and cela in muscle tissue. All these genes were upregulated in the slow-growing fish. Furthermore, weighted gene co-expression network analyses, including three modules (yellow, turquoise, and brown), significantly associated with growth. A network map that included these three modules enabled the identification of a series of hub genes, including rp13a, ube2o, h6pd, etc. These genes may be key candidate genes regulating the growth of rock carp. This study contributes to a deeper understanding of the growth control mechanism in rock carp and offers a scientific basis for efficient breeding and species improvement.

Transcriptome analysis of large yellow croaker (Larimichthys crocea) at different growth rates

The unsynchronized growth of the large yellow croaker (Larimichthys crocea), which impacts growth efficiency, poses a challenge for aquaculture practitioners. In our study, juvenile stocks of large yellow croaker were sorted by size after being cultured in offshore cages for 4 months. Subsequently, individuals from both the fast-growing (FG) and slow-growing (SG) groups were sampled for analysis. High-throughput RNA-Seq was employed to identify genes and pathways that are differentially expressed during varying growth rates, which could suggest potential physiological mechanisms that influence growth rate. Our transcriptome analysis identified 382 differentially expressed genes (DEGs), comprising 145 upregulated and 237 downregulated genes in comparison to the SG group. GO and KEGG enrichment analyses indicated that these DEGs are predominantly involved in signal transduction and biochemical metabolic pathways. Quantitative PCR (qPCR) results demonstrated that cat, fasn, idh1, pgd, fgf19, igf2, and fads2 exhibited higher expression levels, whereas gadd45b and gadd45g showed lower expression compared to the slow-growing group. In conclusion, the differential growth rates of large yellow croaker are intricately associated with cellular proliferation, metabolic rates of the organism, and immune regulation. These findings offer novel insights into the molecular mechanisms and regulatory aspects of growth in large yellow croaker and enhance our understanding of growth-related genes.

Eicosapentaenoic acid-mediated activation of PGAM2 regulates skeletal muscle growth and development via the PI3K/AKT pathway

Eicosapentaenoic acid regulates glucose uptake in skeletal muscle and significantly affects whole-body energy metabolism. However, the underlying molecular mechanism remains unclear. Here we report that eicosapentaenoic acid activates phosphoglycerate mutase 2, which mediates the conversion of 2-phosphoglycerate into 3-phosphoglycerate. This enzyme plays a pivotal role in glycerol degradation, thereby facilitating the proliferation and differentiation of satellite cells in skeletal muscle. Interestingly, phosphoglycerate mutase 2 inhibits mitochondrial metabolism, promoting the formation of fast-type muscle fibers. Treatment with eicosapentaenoic acid and phosphoglycerate mutase 2 knockdown induced opposite transcriptomic changes, most of which were enriched in the PI3K-AKT signaling pathway. Phosphoglycerate mutase 2 activated the PI3K-AKT signaling pathway, which inhibited the phosphorylation of FOXO1, and, in turn, inhibited mitochondrial function and promoted the formation of fast-type muscle fibers. Our results suggest that eicosapentaenoic acid promotes skeletal muscle growth and regulates glucose metabolism by targeting phosphoglycerate mutase 2 and activating the PI3K/AKT signaling pathway.

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.

Transcriptome analysis of the muscle of fast- and slow-growing phoenix barb (Spinibarbus denticulatus denticulatus)

Growth rate is a commercial trait in aquaculture that is influenced by multiple factors, among which genetic composition plays a fundamental role in the growth rate of species. The phoenix barb (Spinibarbus denticulatus denticulatus) is a widely distributed freshwater fish species in South China. Although S. d. denticulatus is reared in South China, the molecular mechanisms underlying the growth rate of the species remain unclear. Here, the authors performed transcriptome analysis of muscle tissues from fast-growing (FG) and slow-growing (SG) S. d. denticulatus at 90, 150, and 300 days after hatch (DAH) to elucidate its growth mechanism. Kyoto Encyclopedia of Genes and Genomes (KEGG) analysis revealed that differentially expressed genes (DEGs) between the two groups were enriched in pathways related to muscle growth, glycolysis, and energy and lipid metabolism. Nonetheless, a higher number of DEGs were identified in the FG vs. SG groups at 90 and 300 DAH compared with 150 DAH. DEGs identified at 90 DAH were mainly enriched in the GH/IGF axis, PI3K-Akt signalling pathway, AMPK signalling pathway and lipid metabolism highly expressed in FG individuals. DEGs identified at 300 DAH were mainly enriched in PI3K-Akt signalling pathway, glycolysis/gluconeogenesis, gene translation and lipid metabolism. In addition, some genes were expressed during the early growth stage in FG individuals but expressed during the late stage in SG individuals, indicating considerable variations in the expression profiles of growth-related genes at different developmental stages. Overall, these findings contribute to the understanding of the growth mechanism of S. d. denticulatus, which would be useful for the propagation of fast-growing breeds.

Growth regulation in brook charr Salvelinus fontinalis

Fish growth can be modulated through genetic selection. However, it is not known whether growth regulatory mechanisms modulated by genetic selection can provide information about phenotypic growth variations among families or populations. Following a five-generation breeding program that selected for the absence of early sexual maturity and increased growth in brook charr we aimed to understand how the genetic selection process modifies the growth regulatory pathway of brook charr at the molecular level. To achieve this, we studied the regulation of growth traits at three different levels: 1) between lines-one under selection, the other not, 2) among-families expressing differences in average growth phenotypes, which we termed family performance, and 3) among individuals within families that expressed extreme growth phenotypes, which we termed slow- and fast-growing. At age 1+, individuals from four of the highest performing and four of the lowest performing families in terms of growth were sampled in both the control and selected lines. The gene expression levels of three reference and ten target genes were analyzed by real-time PCR. Results showed that better growth performance (in terms of weight and length at age) in the selected line was associated with an upregulation in the expression of genes involved in the growth hormone (GH)/insulin growth factor-1 (IGF-1) axis, including the igf-1 receptor in pituitary; the gh-1 receptor and igf-1 in liver; and ghr and igf-1r in white muscle. When looking at gene expression within families, family performance and individual phenotypes were associated with upregulations of the leptin receptor and neuropeptid Y-genes related to appetite regulation-in the slower-growing phenotypes. However, other genes related to appetite (ghrelin, somatostatin) or involved in muscle growth (myosin heavy chain, myogenin) were not differentially expressed. This study highlights how transcriptomics may improve our understanding of the roles of different key endocrine steps that regulate physiological performance. Large variations in growth still exist in the selected line, indicating that the full genetic selection potential has not been reached.

Integrative analysis of miRNA and mRNA expression associated with the immune response in the intestine of rainbow trout (Oncorhynchus mykiss) infected with infectious hematopoietic necrosis virus

Rainbow trout (Oncorhynchus mykiss), an economically important cold-water fish cultured worldwide, suffers from infectious hematopoietic necrosis virus (IHNV) infection, resulting in huge financial losses. In order to understand the immune response of rainbow trout during virus infection, we explored trout intestine transcriptome profiles following IHNV challenge, and identified 3355 differentially expressed genes (DEGs) and 80 differentially expressed miRNAs (DEMs). Transcriptome analysis revealed numerous DEGs involved in immune responses, such as toll-like receptor 3 (TLR3), toll-like receptor 7/8 (TLR7/8), tripartite motif-containing 25 (TRIM25), DExH-Box helicase 58 (DHX58), interferon-induced with helicase C domain 1 (IFIH1), interferon regulatory factor 3 (IRF3/7), signal transducer and activator of transcription 1 (STAT1) and heat shock protein 90-alpha 1 (HSP90A1). Integrated analysis identified five key miRNAs (miR-19-y, miR-181-z, miR-203-y, miR-143-z and miR-206-y) targeting at least two important immune genes (TRIM25, DHX58, STAT1, TLR7/8 and HSP90A1). Gene Ontology and Kyoto Encyclopedia of Genes and Genomes analyses showed that DEGs and target genes were significantly enriched in various immune-related terms including immune system process, binding, cell part and pathways of Toll-like receptor signalling, RIG-I-like receptor signalling, NOD-like receptor signalling, JAK-STAT signalling, PI3K-Akt signalling, NF-kappa B signalling, IL-17 signalling and AGE-RAGE signalling. In addition, protein-protein interaction networks (PPI) was used to display highly interactive DEG networks involving eight immune-related pathways. The expression trends of 12 DEGs and 10 DEMs were further verified by quantitative real-time PCR, which confirmed the reliability of the transcriptome sequencing results. This study expands our understanding of the immune response of rainbow trout infected with IHNV, and provides valuable resources for future studies on the immune molecular mechanism and disease resistance breeding.

Efficacy of Utilization of All-Plant-Based and Commercial Low-Fishmeal Feeds in Two Divergently Selected Strains of Rainbow Trout (Oncorhynchus mykiss): Focus on Growth Performance, Whole-Body Proximate Composition, and Intestinal Microbiome

The present study aimed to investigate the growth performance, whole-body proximate composition, and intestinal microbiome of rainbow trout strains when selected and non-selected for weight gain on all-plant protein diets. A 2x2 factorial design was applied, where a selected (United States) and a non-selected (ITA) rainbow trout strain were fed using either an all-plant protein (PP) or a commercial low-FM diet (C). Diets were fed to five replicates of 20 (PP) or 25 (C) fish for 105 days. At the end of the trial, growth parameters were assessed, and whole fish (15 pools of three fish/diet) and gut samples (six fish/diet) were collected for whole-body proximate composition and gut microbiome analyses, respectively. Independent of the administered diet, the United States strain showed higher survival, final body weight, weight gain, and specific growth rate when compared to the ITA fish (p < 0.001). Furthermore, decreased whole-body ether extract content was identified in the PP-fed United States rainbow trout when compared to the ITA strain fed the same diet (p < 0.001). Gut microbiome analysis revealed the Cetobacterium probiotic-like genus as clearly associated with the United States rainbow trout, along with the up-regulation of the pathway involved in starch and sucrose metabolism. In summary, the overall improvement in growth performance and, to a lesser extent, whole-body proximate composition observed in the selected rainbow trout strain was accompanied by specific, positive modulation of the intestinal microbiome.

Increased accuracy of genomic predictions for growth under chronic thermal stress in rainbow trout by prioritizing variants from GWAS using imputed sequence data

Through imputation of genotypes, genome-wide association study (GWAS) and genomic prediction (GP) using whole-genome sequencing (WGS) data are cost-efficient and feasible in aquaculture breeding schemes. The objective was to dissect the genetic architecture of growth traits under chronic heat stress in rainbow trout (Oncorhynchus mykiss) and to assess the accuracy of GP based on imputed WGS and different preselected single nucleotide polymorphism (SNP) arrays. A total of 192 and 764 fish challenged to a heat stress experiment for 62 days were genotyped using a customized 1 K and 26 K SNP panels, respectively, and then, genotype imputation was performed from a low-density chip to WGS using 102 parents (36 males and 66 females) as the reference population. Imputed WGS data were used to perform GWAS and test GP accuracy under different preselected SNP scenarios. Heritability was estimated for body weight (BW), body length (BL) and average daily gain (ADG). Estimates using imputed WGS data ranged from 0.33 ± 0.05 to 0.55 ± 0.05 for growth traits under chronic heat stress. GWAS revealed that the top five cumulatively SNPs explained a maximum of 0.94%, 0.86% and 0.51% of genetic variance for BW, BL and ADG, respectively. Some important functional candidate genes associated with growth-related traits were found among the most important SNPs, including signal transducer and activator of transcription 5B and 3 (STAT5B and STAT3, respectively) and cytokine-inducible SH2-containing protein (CISH). WGS data resulted in a slight increase in prediction accuracy compared with pedigree-based method, whereas preselected SNPs based on the top GWAS hits improved prediction accuracies, with values ranging from 1.2 to 13.3%. Our results support the evidence of the polygenic nature of growth traits when measured under heat stress. The accuracies of GP can be improved using preselected variants from GWAS, and the use of WGS marginally increases prediction accuracy.

Comparative Physiological and Transcriptomic Profiling Offers Insight into the Sexual Dimorphism of Hepatic Metabolism in Size-Dimorphic Spotted Scat (Scatophagus argus)

The spotted scat (Scatophagus argus) is an economically important cultured marine fish that exhibits a typical sexual size dimorphism (SSD). SSD has captivated considerable curiosity for farmed fish production; however, up till now the exact underlying mechanism remains largely unclear. As an important digestive and metabolic organ, the liver plays key roles in the regulation of fish growth. It is necessary to elucidate its significance as a downstream component of the hypothalamic-pituitary-liver axis in the formation of SSD. In this study, the liver physiological differences between the sexes were evaluated in S. argus, and the activity of several digestive and metabolic enzymes were affected by sex. Females had higher amylase, protease, and glucose-6-phosphate dehydrogenase activities, while males exhibited markedly higher hepatic lipase and antioxidant enzymes activities. A comparative transcriptomics was then performed to characterize the responsive genes. Illumina sequencing generated 272.6 million clean reads, which were assembled into 79,115 unigenes. A total of 259 differentially expressed genes were identified and a few growth-controlling genes such as igf1 and igfbp1 exhibited female-biased expression. Further analyses showed that several GO terms and pathways associated with metabolic process, particularly lipid and energy metabolisms, were significantly enriched. The male liver showed a more active mitochondrial energy metabolism, implicating an increased energy expenditure associated with reproduction. Collectively, the female-biased growth dimorphism of S. argus may be partially attributed to sexually dimorphic metabolism in the liver. These findings would facilitate further understanding of the nature of SSD in teleost fish.

Alternative migratory tactics in brown trout (Salmo trutta) are underpinned by divergent regulation of metabolic but not neurological genes

The occurrence of alternative morphs within populations is common, but the underlying molecular mechanisms remain poorly understood. Many animals, for example, exhibit facultative migration, where two or more alternative migratory tactics (AMTs) coexist within populations. In certain salmonid species, some individuals remain in natal rivers all their lives, while others (in particular, females) migrate to sea for a period of marine growth. Here, we performed transcriptional profiling ("RNA-seq") of the brain and liver of male and female brown trout to understand the genes and processes that differentiate between migratory and residency morphs (AMT-associated genes) and how they may differ in expression between the sexes. We found tissue-specific differences with a greater number of genes expressed differentially in the liver (n = 867 genes) compared with the brain (n = 10) between the morphs. Genes with increased expression in resident livers were enriched for Gene Ontology terms associated with metabolic processes, highlighting key molecular-genetic pathways underlying the energetic requirements associated with divergent migratory tactics. In contrast, smolt-biased genes were enriched for biological processes such as response to cytokines, suggestive of possible immune function differences between smolts and residents. Finally, we identified evidence of sex-biased gene expression for AMT-associated genes in the liver (n = 12) but not the brain. Collectively, our results provide insights into tissue-specific gene expression underlying the production of alternative life histories within and between the sexes, and point toward a key role for metabolic processes in the liver in mediating divergent physiological trajectories of migrants versus residents.

Identifying Candidate Genes Involved in the Regulation of Early Growth Using Full-Length Transcriptome and RNA-Seq Analyses of Frontal and Parietal Bones and Vertebral Bones in Bighead Carp (Hypophthalmichthys nobilis)

Growth, one of the most important traits monitored in domestic animals, is essentially associated with bone development. To date, no large-scale transcriptome studies investigating bone development in bighead carp have been reported. In this study, we applied Isoform-sequencing technology to uncover the entire transcriptomic landscape of the bighead carp (Hypophthalmichthys nobilis) in early growth stage, and obtained 63,873 non-redundant transcripts, 20,907 long non-coding RNAs, and 1,579 transcription factors. A total of 381 alternative splicing events were seen in the frontal and parietal bones with another 784 events simultaneously observed in the vertebral bones. Coupling this to RNA sequencing (RNA-seq) data, we identified 27 differentially expressed unigenes (DEGs) in the frontal and parietal bones and 45 DEGs in the vertebral bones in the fast-growing group of fish, when compared to the slow-growing group of fish. Finally, 15 key pathways and 20 key DEGs were identified and found to be involved in regulation of early growth such as energy metabolism, immune function, and cytoskeleton function and important cellular pathways such as the arginine and proline metabolic pathway (p4ha1), FoxO signaling pathway (sgk1), cell adhesion molecules (b2m, ptprc, and mhcII), and peroxisome proliferator-activated receptor signaling pathway (scd). We established a novel full-length transcriptome resource and combined it with RNA-seq to elucidate the mechanism of genetic regulation of differential growth in bighead carp. The key DEGs identified in this study could fuel further studies investigating associations between growth and bone development and serve as a source of potential candidate genes for marker-assisted breeding programs.

Replacing fish oil and astaxanthin by microalgal sources produced different metabolic responses in juvenile rainbow trout fed 2 types of practical diets

Dietary fish oil supplementation provides n-3 long-chained polyunsaturated fatty acids for supporting fish growth and metabolism and enriching fillet with eicosapentaenoic acid (EPA; 20:5n-3) and docosahexaenoic acid (DHA; c22:6n-3). Two experiments were performed as a 3 × 2 factorial arrangement of dietary treatments for 16 wk to determine effects and mechanisms of replacing 0%, 50%, and 100% fish oil with DHA-rich microalgae in combination with synthetic vs. microalgal source of astaxanthin in plant protein meal (PM)- or fishmeal (FM)- based diets for juvenile rainbow trout (Oncorhynchus mykiss). Fish (22 ± 0.26 g) were stocked at 17/tank and 3 tanks/diet. The 100% fish oil replacement impaired (P < 0.0001) growth performance, dietary protein and energy utilization, body indices, and tissue accumulation of DHA and EPA in both diet series. The impairments were associated (P < 0.05) with upregulation of hepatic gene expression related to growth (ghr1and igf1) and biosynthesis of DHA and EPA (fads6 and evol5) that was more dramatic in the FM than PM diet-fed fish, and more pronounced on tissue EPA than DHA concentrations. The source of astaxanthin exerted interaction effects with the fish oil replacement on several measures including muscle total cholesterol concentrations. In conclusion, replacing fish oil by the DHA-rich microalgae produced more negative metabolic responses than the substitution of synthetic astaxanthin by the microalgal source in juvenile rainbow trout fed 2 types of practical diets.