Transcriptome Analysis of the Liver and Muscle Tissues of Black Carp (Mylopharyngodon piceus) of Different Growth Rates

Hepatopancreas Transcriptome Analysis of Spinibarbus sinensis to Reveal Different Growth-Related Genes

Spinibarbus sinensis, also known as Qingbo, is an important economic fish in China. However, the detailed mechanisms underlying its growth are still unknown. To excavate the genes and signaling pathways related to its growth, we compared the transcriptome profiles of the hepatopancreas tissues of S. sinensis, with two groups of growth rate for evaluation. An average of 66,304,909 and 68,739,585 clean reads were obtained in the fast growth (FG) and slow growth (SG) group, respectively. The differential gene expression analysis results showed that 272 differentially expressed genes (DEGs) were screened between the FG and SG groups, including 101 up-regulated genes and 171 down-regulated genes. Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analysis results showed that GO terms related to metabolic process, organic substance metabolic process, and catalytic activity were enriched, pathway signals related to steroid biosynthesis and protein digestion and absorption were also detected. Meanwhile, the potential key regulatory genes sst2, fndc4, and cckra related to the growth of S. sinensis were screened. Reverse transcript fluorescence quantitative PCR (RT-qPCR) validation of 18 DEGs associated with growth differences showed that the RT-qPCR results were consistent with RNA-seq analysis, and nine genes, stk31, gpr149, angptl1, fstl1, sik1, ror2, nlrc3, pdlim2, and nav2 were significantly expressed in the FG group. bmp1, stc1, gpatch8, sstrt2, s100a1, ktf6, cckar6, sync1, bhlha15, a total of nine genes were significantly expressed in the SG group. This study provides basic information for improving the growth characteristics of S. sinensis and the functional research of candidate genes.

Transcriptome Analysis of Different Aquaculture Substrates on the Immune Response of Babylonia areolata

To assess the impact of different substrates in a recirculating water system on the immune response and antioxidant capacity of Babylonia areolata, we conducted a comparative analysis of the transcriptomes and antioxidant performance of the digestive glands in three substrate environments (sand-S group, ceramic granules-C group, and PVC breeding nest-P group). Transcriptome results revealed that the S group and P group exhibited the highest number of differentially expressed genes (DEGs), with a total of 2218 DEGs, including 928 upregulated and 1290 downregulated DEGs. The C group and P group had 1055 DEGs in common, with 316 upregulated and 739 downregulated DEGs. The C group and S group had the fewest DEGs, with 521 in total, including 303 upregulated and 218 downregulated DEGs. GO enrichment analysis showed that in the S vs P group, terms such as catalytic activity, membrane part, and cellular process were enriched with 287, 262, and 180 DEGs, respectively. In the C vs P group, binding, cellular process, and cell part were enriched with 146, 135, and 127 DEGs, respectively. In the C vs S group, catalytic activity, membrane part, and metabolic process were enriched with 90, 83, and 59 DEGs, respectively. Kegg enrichment analysis revealed significant changes in immune-related pathways in the S vs P group, including lysosome, phagosome, and leukocyte transendothelial migration, with 30, 13, and 10 enriched DEGs, respectively. In the C vs P group, phagosome, drug metabolism-other enzymes, and N-Glycan biosynthesis showed significant changes in immune-related pathways, with 9, 6, and 4 enriched DEGs, respectively. In the C vs S group, lysosome, PPAR signaling pathway, and fatty acid degradation exhibited significant changes in immune-related pathways, with 8, 4, and 3 enriched DEGs, respectively. Regarding antioxidant capacity, the S group showed significantly higher total T-AOC than the other experimental groups, while CAT, SOD, POD, and AKP were lower than in the C and P groups. The ACP level in the Sand group was not significantly different from the P group but significantly lower than the C group. In conclusion, substrate environments significantly influence the immune-related genes and key antioxidant enzyme activities in B. areolata.

Single nucleotide polymorphism SNP19140160 A > C is a potential breeding locus for fast-growth largemouth bass (Micropterus salmoides)

BACKGROUND

Largemouth bass (Micropterus salmoides) has significant economic value as a high-yielding fish species in China's freshwater aquaculture industry. Determining the major genes related to growth traits and identifying molecular markers associated with these traits serve as the foundation for breeding strategies involving gene pyramiding. In this study, we screened restriction-site associated DNA sequencing (RAD-seq) data to identify single nucleotide polymorphism (SNP) loci potentially associated with extreme growth differences between fast-growth and slow-growth groups in the F1 generation of a largemouth bass population.

RESULTS

We subsequently identified associations between these loci and specific candidate genes related to four key growth traits (body weight, body length, body height, and body thickness) based on SNP genotyping. In total, 4,196,486 high-quality SNPs were distributed across 23 chromosomes. Using a population-specific genotype frequency threshold of 0.7, we identified 30 potential SNPs associated with growth traits. Among the 30 SNPs, SNP19140160, SNP9639603, SNP9639605, and SNP23355498 showed significant associations; three of them (SNP9639603, SNP9639605, and SNP23355498) were significantly associated with one trait, body length, in the F1 generation, and one (SNP19140160) was significantly linked with four traits (body weight, height, length, and thickness) in the F1 generation. The markers SNP19140160 and SNP23355498 were located near two growth candidate genes, fam174b and ppip5k1b, respectively, and these candidate genes were closely linked with growth, development, and feeding. The average body weight of the group with four dominant genotypes at these SNP loci in the F1 generation population (703.86 g) was 19.63% higher than that of the group without dominant genotypes at these loci (588.36 g).

CONCLUSIONS

Thus, these four markers could be used to construct a population with dominant genotypes at loci related to fast growth. These findings demonstrate how markers can be used to identify genes related to fast growth, and will be useful for molecular marker-assisted selection in the breeding of high-quality largemouth bass.

Multiple Tissues Transcriptome of Zig-Zag Eel (Mastacembelus armatus) with Different Growth Rates

In order to explore the main regulatory genes and related pathways of growth traits, transcriptome sequencing was first performed on the brain, liver, and muscle tissues of 3-month-old M. armatus with different growth rates. By comparative transcriptome analysis of fast-growing and slow-growing groups of M. armatus, a total of 2887 DEGs were screened, of which 59 up-regulated genes and 105 down-regulated genes were detected in the brain, 146 up-regulated genes and 202 down-regulated genes were detected in the liver, and 529 up-regulated genes and 1846 down-regulated genes were detected in muscle, including insulin-like growth factor binding protein 1a (IGFBP1A), insulin-like growth factor binding protein 1b (IGFBP1B), myosin, light chain 1 (MYL1), and myoglobin (MB). Through Gene Ontology (GO) and the Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analysis, we identified a total of 288 significantly enriched GO entries and 68 significantly enriched KEGG pathways related to growth, such as skeletal muscle tissue development, insulin-like growth factor binding, and the mitotic cell cycle. These key genes and signaling pathways may play a key role in regulating the growth of M. armatus. Digging into the regulatory mechanisms of these key genes will provide a theoretical basis for further exploration of the molecular mechanisms related to the growth and development of M. armatus, and help to breed new varieties of M. armatus with rapid growth.

Muscle transcriptome analysis provides new insights into the growth gap between fast- and slow-growing Sinocyclocheilus grahami

Sinocyclocheilus grahami is an economically valuable and famous fish in Yunnan Province, China. However, given its slow growth (40 g/2 years) and large growth differences among individuals, its growth performance needs to be improved for sustainable future use, in which molecular breeding technology can play an important role. In the current study, we conducted muscle transcriptomic analysis to investigate the growth gaps among individuals and the mechanism underlying growth within 14 fast- and 14 slow-growth S. grahami. In total, 1,647 differentially expressed genes (DEGs) were obtained, including 947 up-regulated and 700 down-regulated DEGs in fast-growth group. Most DEGs were significantly enriched in ECM-receptor interaction, starch and sucrose metabolism, glycolysis/gluconeogenesis, pyruvate metabolism, amino acids biosynthesis and metabolism, peroxisome, and PPAR signaling pathway. Some genes related to glycogen degradation, glucose transport, and glycolysis (e.g., adipoq, prkag1, slc2a1, agl, pygm, pgm1, pfkm, gapdh, aldoa, pgk1, pgam2, bpgm, and eno3) were up-regulated, while some genes related to fatty acid degradation and transport (e.g., acox1, acaa1, fabp1b.1, slc27a1, and slc27a2) and amino acid metabolism (e.g., agxt, shmt1, glula, and cth) were down-regulated in the fast-growth group. Weighted gene co-expression network analysis identified col1a1, col1a2, col5a1, col6a2, col10a1, col26a1, bglap, and krt15 as crucial genes for S. grahami growth. Several genes related to bone and muscle growth (e.g., bmp2, bmp3, tgfb1, tgfb2, gdf10, and myog) were also up-regulated in the fast-growth group. These results suggest that fast-growth fish may uptake adequate energy (e.g., glucose, fatty acid, and amino acids) from fodder, with excess energy substances used to synthesize collagen to accelerate bone and muscle growth after normal life activities are maintained. Moreover, energy uptake may be the root cause, while collagen synthesis may be the direct reason for the growth gap between fast- and slow-growth fish. Hence, improving food intake and collagen synthesis may be crucial for accelerating S. grahami growth, and further research is required to fully understand and confirm these associations.

Impact of Titanium Dioxide-Graphene Oxide (TiO2-GO) Composite Nanoparticle on the Juveniles of the Giant River Prawn, Macrobrachium rosenbergii: Physio-Biochemistry and Transcriptional Response

Nanomaterials are used in many fields, resulting in inevitably releasing into the aquatic environment. The presence of nanomaterials, including TiO₂-GO in the aquatic environment, can be toxic to aquatic organisms. However, few studies have focused on the effects of TiO₂-GO composite nanoparticle on crustaceans. In the present study, the giant river prawn Macrobrachium rosenbergii juveniles were exposed to two concentrations of TiO₂-GO composite nanoparticle (0.1 and 0.5 mg/L). The effects of TiO₂-GO composite exposure on activities of digestive and antioxidant-related enzymes and expressions of growth and immune-related genes at the transcriptome were studied. The results showed that the survival rate and growth performance were not negatively affected by TiO₂-GO composite at the two exposure levels. Nevertheless, exposure to TiO₂-GO composite causes an effect on the activities of digestive and antioxidant enzymes in the juvenile prawns. The enzyme activities of CAT, SOD, GSH-Px, AMS, TPS, and LPS in the 0.1 mg/L TiO₂-GO composite experimental group were markedly reduced than those in the control group. Additionally, the expression level of genes involved in growth and immunity was significantly affected by TiO₂-GO composite. After exposure to the 0.1 mg/L TiO₂-GO composite, the mRNA expression level of MSTN was significantly increased, but the level of EcR, Raptor, and CaBP was significantly decreased. However, the mRNA levels of the CTL, TLR, JAK, and STAT were significantly increased after exposure to the 0.5 mg/L concentration of TiO₂-GO composite. Furthermore, to understand the molecular mechanism of M. rosenbergii under TiO₂-GO composite exposure, RNA-Seq was employed to analyze the changes of the muscle and hepatopancreas transcriptome. Compared with the control group, we identified 5166 and 4784 differentially expressed genes (DEGs) in the muscle and hepatopancreas, respectively (p < 0.05). Based on gene ontology and KEGG analysis, significant differences were observed in the DEGs involved in activity and binding, metabolism, immune response, and environmental information processing. These results showed that exposure to TiO₂-GO composite nanoparticle led to the changes of enzyme activity and gene expression, suggesting that TiO₂-GO composite existing in aquatic environments would disrupt the physiology of M. rosenbergii. This study will serve as a foundation for subsequent research into the evaluation of nanomaterial toxicity on crustacean species.

Transcriptome analysis reveals the involvement of ubiquitin-proteasome pathway in the regulation of muscle growth of rice flower carp

Growth mechanism of economically important aquaculture species has aroused widespread interest among scholars. Rice flower carp (Cyprinus carpio), commonly cultured in rice-fish farming systems, shows wide variation in body mass at the same age, which limits the development of commercial aquaculture. In this study, muscle tissues from 20-month-old fish of different sizes were used for transcriptome analysis and muscle histological studies. The muscle histological analysis showed the muscle growth in rice flower carp main depends on the hypertrophic growth of muscle fibers. A total of 30,590 unigenes were generated by muscle trancriptome analysis, including 403 differentially expressed genes (DEGs). Of these, 157 DEGs were upregulated and 246 DEGs were downregulated. Nine unigenes related to the ubiquitin-proteasome pathway were identified using differential expression analysis. This study initially revealed that the differences in growth of rice flower carp could be due to hypertrophic growth of muscle fibers caused by higher protein deposition, and the ubiquitin-proteasome pathway was an important factor affecting the growth rate of rice flower carp. E3 ubiquitin-protein ligase ari7, g2e3, Neurl1 and rnf144ab were upregulated in the slow-growing fish, indicating the binding of ubiquitin to target protein was enhanced. Foxo3 was upregulated in the slow-growing fish, which could promote the muscle loss. Eif4a2 was upregulated in the fast-growing fish, increasing protein translation efficiency. Some genes related to active muscle contraction such as actb, actg, camk2a, and camk2b were upregulated in the fast-growing rice flower carp muscle. In summary, these results provide valuable information about the key genes for use as biomarkers of growth in selective breeding programs for rice flower carp and provide novel insights into the regulatory mechanisms of muscle growth.

Comprehensive analysis of hub mRNA, lncRNA and miRNA, and associated ceRNA networks implicated in grass carp (Ctenopharyngodon idella) growth traits

Grass carp (Ctenopharyngodon idella) is the most productive freshwater aquaculture fish in worldwide. However, the molecular mechanism of its growth traits has not been fully elucidated. Whole transcriptome analysis of the brain and hepatopancreas of 29 six-month-old grass carp with different growth rates was performed. Weighted gene co-expression network analysis (WGCNA) was used to construct a weighted gene co-expression network of mRNA, miRNA and lncRNA separately. A total of 35 hub mRNAs, 47 hub lncRNAs and 33 hub miRNAs were identified from the brain, 37 hub mRNAs, 110 hub lncRNAs and 36 hub miRNAs were identified from the hepatopancreas. The ceRNA networks in the brain and hepatopancreas were involved in brain development and nutrient metabolism, respectively. Overall, this is the first investigation of the growth-related transcriptomic characteristics in the brain and hepatopancreas of grass carp, thus will help us to further explore the molecular mechanism of grass carp growth rate.

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.

RNA-Seq Transcriptome Analysis of the Liver and Brain of the Black Carp (Mylopharyngodon piceus) During Fasting

The black carp (Mylopharyngodon piceus) is an important carnivorous freshwater-cultured species. To understand the molecular basis underlying the response of black carp to fasting, we used RNA-Seq to analyze the liver and brain transcriptome of fasting fish. Annotation to the NCBI database identified 66,609 unigenes, of which 22,841 were classified into the Gene Ontology database and 15,925 were identified in the Kyoto Encyclopedia of Genes and Genomes (KEGG) database. Comparative analysis of the expression profile between fasting and normal feeding fish revealed 13,737 differentially expressed genes (P < 0.05), of which 12,480 were found in liver tissue and 1257 were found in brain tissue. The KEGG pathway analysis showed significant differences in expression of genes involved in metabolic and immune pathways, such as the insulin signaling pathway, PI3K-Akt signaling pathway, cAMP signaling pathway, FoxO signaling pathway, AMPK signaling pathway, endocytosis, and apoptosis. Quantitative real-time PCR analysis confirmed that expression of the genes encoding the factors involved in those pathways differed between fasting and feeding fish. These results provide valuable information about the molecular response mechanism of black carp under fasting conditions.