Transcriptome analysis of grass carp (Ctenopharyngodon idella) between fast- and slow-growing fish

Muscle Fiber Characteristics and Transcriptome Analysis in Slow- and Fast-Growing Megalobrama amblycephala

Growth is an important trait in aquaculture that is influenced by various factors, among which genetic regulation plays a crucial role. Megalobrama amblycephala, one of the most important freshwater species in China, exhibits wide variations in body mass among individuals of the same age within the same pool. But the molecular mechanisms underlying wide variation in body mass remain unclear. Here, we performed muscle histological and transcriptome analysis of muscle tissues from Fast-Growing (FG) and Slow-Growing (SG) M. amblycephala at the age of 4 months old (4 mo) and 10 months old (10 mo) to elucidate its muscle development and growth mechanism. The muscle histological analysis showed smaller diameter and higher total number of muscle fibers in FG compared to SG at 4 mo, while larger diameter and total number of muscle fibers were detected in FG at 10 mo. The transcriptome analysis of muscle tissue detected 1171 differentially expressed genes (DEGs) between FG and SG at 4 mo, and 718 DEGs between FG and SG at 10 mo. Furthermore, 44 DEGs were consistently up-regulated in FG at both 4 mo and 10 mo. Up-regulated DEGs in FG at 4 mo were mainly enriched in the pathways related to cell proliferation, while down-regulated DEGs were significantly enriched in cell fusion and muscle contraction. Up-regulated DEGs in FG at 10 mo were mainly enriched in the pathways related to cell proliferation and protein synthesis. Therefore, these results provide novel insights into the molecular mechanism of M. amblycephala muscle growth at different stages, and will be of great guiding significance to promote the fast growth of M. amblycephala.

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.

Transcriptome-based nutrigenomics analysis reveals the roles of dietary taurine in the muscle growth of juvenile turbot (Scophthalmus maximus)

The present study explored transcriptomics and gene regulation variations in the muscle of turbot fed with dietary taurine. A 70-day feeding trial was conducted using turbot (initial body weight: 3.66 ± 0.02 g) fed with different levels of dietary taurine: 0 % (C), 0.4 % (T2), 1.2 % (T4) and 2.0 % (T6). Two methods were used to analyze and verify the taurine effects on muscle growth: (1) real-time quantitative PCR (qRT-PCR) for the key muscle growth-related genes and (2) transcriptomic analysis by next-generation sequencing (NGS). The results showed that 1.2 % of dietary taurine supplementation significantly increased the expression of muscle growth stimulatory genes, including TauT, myoD, Myf5, myogenin and follistatin. And also, the 1.2 % level significantly decreased the expression of the muscle growth-restricting gene (myostatin). Meanwhile, transcriptomics analysis found that 1.2 % dietary taurine supplementation significantly increased the number of up-regulated genes linked to metabolic pathways. In contrast, taurine significantly enriched the actin cytoskeleton and metabolic pathways in the T4 and T2 groups, respectively. These findings align with the gene ontology (GO) analysis, which indicated a higher number of cellular component (CC) gene expressions at a 1.2 % of dietary taurine compared to a 0.4 % of dietary taurine supplementation. In conclusion, dietary taurine had positive impacts on the growth-stimulatory genes. Moreover, 1.2 % of dietary taurine supplementation is important to the metabolic pathway enrichment.

Integrated analysis of mRNAs and lncRNAs reveals candidate marker genes and potential hub lncRNAs associated with growth regulation of the Pacific Oyster, Crassostrea gigas

BACKGROUND

The Pacific oyster, Crassostrea gigas, is an economically important shellfish around the world. Great efforts have been made to improve its growth rate through genetic breeding. However, the candidate marker genes, pathways, and potential lncRNAs involved in oyster growth regulation remain largely unknown. To identify genes, lncRNAs, and pathways involved in growth regulation, C. gigas spat was cultured at a low temperature (15 ℃) to yield a growth-inhibited model, which was used to conduct comparative transcriptome analysis with spat cultured at normal temperature (25 ℃).

RESULTS

In total, 8627 differentially expressed genes (DEGs) and 1072 differentially expressed lncRNAs (DELs) were identified between the normal-growth oysters (cultured at 25 ℃, hereinafter referred to as NG) and slow-growth oysters (cultured at 15 ℃, hereinafter referred to as SG). Functional enrichment analysis showed that these DEGs were mostly enriched in the AMPK signaling pathway, MAPK signaling pathway, insulin signaling pathway, autophagy, apoptosis, calcium signaling pathway, and endocytosis process. LncRNAs analysis identified 265 cis-acting pairs and 618 trans-acting pairs that might participate in oyster growth regulation. The expression levels of LNC_001270, LNC_003322, LNC_011563, LNC_006260, and LNC_012905 were inducible to the culture temperature and food abundance. These lncRNAs were located at the antisense, upstream, or downstream of the SREBP1/p62, CDC42, CaM, FAS, and PIK3CA genes, respectively. Furthermore, the expression of the trans-acting lncRNAs, including XR_9000022.2, LNC_008019, LNC_015817, LNC_000838, LNC_00839, LNC_011859, LNC_007294, LNC_006429, XR_002198885.1, and XR_902224.2 was also significantly associated with the expression of genes enriched in AMPK signaling pathway, insulin signaling pathway, autophagy, apoptosis, calcium signaling pathway, and endocytosis process.

CONCLUSIONS

In this study, we identified the critical growth-related genes and lncRNAs that could be utilized as candidate markers to illustrate the molecular mechanisms underlying the growth regulation of Pacific oysters.

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.

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.

Starvation Affects the Muscular Morphology, Antioxidant Enzyme Activity, Expression of Lipid Metabolism-Related Genes, and Transcriptomic Profile of Javelin Goby (Synechogobius hasta)

Fish in natural and cultivated environments can be challenged by starvation. However, inducing starvation in a controlled manner cannot only reduce feed consumption but also reduces aquatic eutrophication and even improves farmed fish quality. This study investigated the effects of starvation on the muscular function, morphology, and regulatory signaling in javelin goby (Synechogobius hasta) by evaluating the biochemical, histological, antioxidant, and transcriptional changes in the musculature of S. hasta subjected to 3, 7, and 14 days fasting. The muscle glycogen and triglyceride levels in S. hasta were gradually reduced under starvation, reaching their lowest at the end of the trial (P < 0.05). The levels of glutathione and superoxide dismutase were significantly elevated after 3-7 days of starvation (P < 0.05), but later returned to the level of the control group. The muscle of starved S. hasta developed structural abnormalities in some areas after 7 days of food deprivation, and more vacuolation and more atrophic myofibers were observed in 14-day fasted fish. The transcript levels of stearoyl-CoA desaturase 1 (scd1), the key gene involved in the biosynthesis of monounsaturated fatty acids, were markedly lower in the groups starved for 7 or more days (P < 0.05). However, the relative expressions of genes associated with lipolysis were decreased in the fasting experiment (P < 0.05). Similar declines in the transcriptional response to starvation were found in muscle fatp1 and ppar γ abundance (P < 0.05). Furthermore, the de novo transcriptome of muscle tissue from the control, 3-day and 14-day starved S. hasta generated 79,255 unigenes. The numbers of differentially expressed genes (DEGs) identified by pairwise comparisons among three groups were 3276, 7354, and 542, respectively. The enrichment analysis revealed that the DEGs were primarily involved in metabolism-related pathways, including ribosome, TCA pathway, and pyruvate metabolism. Moreover, the qRT-PCR results of 12 DEGs validated the expression trends observed in the RNA-seq data. Taken together, these findings demonstrated the specific phenotypical and molecular responses of muscular function and morphology in starved S. hasta, which may offer preliminary reference data for optimizing operational strategies incorporating fasting/refeeding cycles in aquaculture.

Splenic protection network revealed by transcriptome analysis in inactivated vaccine-immunized flounder (Paralichthys olivaceus) against Edwardsiella tarda infection

The protective immune response produced by fish after vaccination is crucial for vaccine effectiveness. Our previous studies have shown inactivated vaccine against Edwardsiella tarda can induce immune response in flounder (Paralichthys olivaceus). To elucidate the protective immune response at the genetic level, in this study, flounder was immunized with inactivated E. tarda for 5 weeks, and then they were challenged with E. tarda. The spleen was dissected at 7th day post immunization, 1st and 7th day post challenge, respectively. Transcriptome analysis showed that average of 46 million clean reads were obtained per library, while percentage of clean reads being mapped to reference genome was more than 89% in all cases, which suggested good quality of samples. As for differentially expressed genes (DEGs) identification in inactivated E. tarda groups, at 7th day post immunization, 1422 DEGs were identified and significantly enriched in innate immune-related pathways, such as Phagosome, Cell adhesion molecules and NF-kappa B signaling pathway; At 1st post challenge day, 1210 DEGs were identified and enriched to Antigen processing and presentation and Cell adhesion molecules, indicating that the pathogen was rapidly recognized and delivered; At 7th post challenge day, 1929 DEGs were identified, belonged to Toll-like receptor signaling pathway, Antigen processing and presentation, Th1 and Th2 cell differentiation and Th17 cell differentiation. Compared to 7th post immunization day, 73 immune-associated DEGs were identified at 1st post challenge day. Protein-protein interaction networks analysis revealed 11 hub genes (TLR7, TLR3, CXCR4, IFIH1, TLR8 etc), associated with recognition of pathogens and activation of innate immunity; while for 7th post challenge day, 141 immune-associated DEGs were identified. 30 hub genes (IL6, STAT1, HSP90A.1, TLR7, IL12β etc) were associated with stimulation of lymphocyte differentiation and activation of cellular immunity. Ten immune-related genes were randomly selected for RT-qPCR validation at each time point. In conclusion, data revealed protection of flounder against E. tarda infection by inactivated vaccine is mediated via immediate recognition of pathogen and subsequently activation of cellular immunity. Results give new aspect for vaccine protection cascades, is good references for vaccine evaluation.

Effects of Dietary Supplementation with Aurantiochytrium sp. on Zebrafish Growth as Determined by Transcriptomics

The marine protist Aurantiochytrium produces several bioactive chemicals, including EPA (eicosapentaenoic acid), DHA (docosahexaenoic acid), and other critical fish fatty acids. It has the potential to improve growth and fatty acid profiles in aquatic taxa. This study evaluated zebrafish growth performance in response to diets containing 1% to 3% Aurantiochytrium sp. crude extract (TE) and single extract for 56 days. Growth performance was best in the 1% TE group, and therefore, this concentration was used for further analyses of the influence of Aurantiochytrium sp. Levels of hepatic lipase, glucose-6-phosphate dehydrogenase, acetyl-CoA oxidase, glutathione peroxidase, and superoxide dismutase increased significantly in response to 1% TE, while malic enzyme activity, carnitine lipid acylase, acetyl-CoA carboxylase, fatty acid synthase, and malondialdehyde levels decreased. These findings suggest that Aurantiochytrium sp. extract can modulate lipase activity, improve lipid synthesis, and decrease oxidative damage caused by lipid peroxidation. Transcriptome analysis revealed 310 genes that were differentially expressed between the 1% TE group and the control group, including 185 up-regulated genes and 125 down-regulated genes. Kyoto Encyclopedia of Genes and Genomes (KEGG) and Gene Ontology (GO) pathway analyses of the differentially expressed genes revealed that Aurantiochytrium sp. extracts may influence liver metabolism, cell proliferation, motility, and signal transduction in zebrafish.

Dynamic Transcriptional Landscape of Grass Carp (Ctenopharyngodon idella) Reveals Key Transcriptional Features Involved in Fish Development

A high-quality baseline transcriptome is a valuable resource for developmental research as well as a useful reference for other studies. We gathered 41 samples representing 11 tissues/organs from 22 important developmental time points within 197 days of fertilization of grass carp eggs in order to systematically examine the role of lncRNAs and alternative splicing in fish development. We created a high-quality grass carp baseline transcriptome with a completeness of up to 93.98 percent by combining strand-specific RNA sequencing and single-molecule real-time RNA sequencing technologies, and we obtained temporal expression profiles of 33,055 genes and 77,582 transcripts during development and tissue differentiation. A family of short interspersed elements was preferentially expressed at the early stage of zygotic activation in grass carp, and its possible regulatory components were discovered through analysis. Additionally, after thoroughly analyzing alternative splicing events, we discovered that retained intron (RI) alternative splicing events change significantly in both zygotic activation and tissue differentiation. During zygotic activation, we also revealed the precise regulatory characteristics of the underlying functional RI events.

Dynamic transcriptome analysis of the muscles in high-fat diet-induced obese zebrafish (Danio rerio) under 5-HT treatment

Peripheral 5-hydroxytryptamine (5-HT, also called serotonin) is reportedly a potential therapeutic target in obesity-related metabolic diseases due to its regulatory role in energy homeostasis in mammals. However, information on the detailed effect of peripheral 5-HT on the energy metabolism in fishes, especially the lipid metabolism, and the underlying mechanism remains elusive. In this study, a diet-induced obesity model was developed in the zebrafish (Danio rerio), a prototypical animal model for metabolic disorders. The zebrafish were fed a high-fat diet for 8 weeks and were simultaneously injected with PBS, 0.1 mM and 10 mM 5-HT, intraperitoneally. The body weight was significantly lower in the zebrafish injected with 0.1 mM 5-HT (P < 0.05), however, there was no change in body length (P > 0.05) at the end of the 8-week treatment. The muscle tissues from the zebrafish treated with PBS and 5-HT were collected for transcriptomic analysis and the RNA-seq revealed 1134, 3713, and 2535 genes were screened out compared to the muscular DEGs among three groups. The enrichment analysis revealed DEGs to be significantly associated with multiple metabolic pathways, including ribosome, oxidative phosphorylation, proteasome, PPAR signaling pathway, and ferroptosis. Additionally, the qRT-PCR validated 12 DEGs out of which 10 genes exhibited consistent trends. Taken together, this data provided useful information on the transcriptional characteristics of the muscle tissue in the obese zebrafish exposed to 5-HT, offering important insights into the regulatory effect of peripheral 5-HT in teleosts, as well as novel approaches for preventing and treating obesity-related metabolic dysfunction.

Assembly-free rapid differential gene expression analysis in non-model organisms using DNA-protein alignment

Background

RNA-seq is being increasingly adopted for gene expression studies in a panoply of non-model organisms, with applications spanning the fields of agriculture, aquaculture, ecology, and environment. For organisms that lack a well-annotated reference genome or transcriptome, a conventional RNA-seq data analysis workflow requires constructing a de-novo transcriptome assembly and annotating it against a high-confidence protein database. The assembly serves as a reference for read mapping, and the annotation is necessary for functional analysis of genes found to be differentially expressed. However, assembly is computationally expensive. It is also prone to errors that impact expression analysis, especially since sequencing depth is typically much lower for expression studies than for transcript discovery.

Results

We propose a shortcut, in which we obtain counts for differential expression analysis by directly aligning RNA-seq reads to the high-confidence proteome that would have been otherwise used for annotation. By avoiding assembly, we drastically cut down computational costs – the running time on a typical dataset improves from the order of tens of hours to under half an hour, and the memory requirement is reduced from the order of tens of Gbytes to tens of Mbytes. We show through experiments on simulated and real data that our pipeline not only reduces computational costs, but has higher sensitivity and precision than a typical assembly-based pipeline. A Snakemake implementation of our workflow is available at: https://bitbucket.org/project_samar/samar.

Conclusions

The flip side of RNA-seq becoming accessible to even modestly resourced labs has been that the time, labor, and infrastructure cost of bioinformatics analysis has become a bottleneck. Assembly is one such resource-hungry process, and we show here that it can be avoided for quick and easy, yet more sensitive and precise, differential gene expression analysis in non-model organisms.

Supplementary Information

The online version contains supplementary material available at (10.1186/s12864-021-08278-7).

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.

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.

Intestinal Transcriptome Analysis Reveals Enrichment of Genes Associated with Immune and Lipid Mechanisms, Favoring Soybean Meal Tolerance in High-Growth Zebrafish (Danio Rerio)

The molecular mechanisms underlying fish tolerance to soybean meal (SBM) remain unclear. Identifying these mechanisms would be beneficial, as this trait favors growth. Two fish replicates from 19 experimental families were fed fishmeal-(100FM) or SBM-based diets supplemented with saponin (50SBM + 2SPN) from juvenile to adult stages. Individuals were selected from families with a genotype-by-environment interaction higher (HG-50SBM + 2SPN, 170 ± 18 mg) or lower (LG-50SBM + 2SPN, 76 ± 10 mg) weight gain on 50SBM + 2SPN for intestinal transcriptomic analysis. A histological evaluation confirmed middle intestinal inflammation in the LG- vs. HG-50SBM + 2SPN group. Enrichment analysis of 665 differentially expressed genes (DEGs) identified pathways associated with immunity and lipid metabolism. Genes linked to intestinal immunity were downregulated in HG fish (mpx, cxcr3.2, cftr, irg1l, itln2, sgk1, nup61l, il22), likely dampening inflammatory responses. Conversely, genes involved in retinol signaling were upregulated (rbp4, stra6, nr2f5), potentially favoring growth by suppressing insulin responses. Genes associated with lipid metabolism were upregulated, including key components of the SREBP (mbtps1, elov5l, elov6l) and cholesterol catabolism (cyp46a1), as well as the downregulation of cyp7a1. These results strongly suggest that transcriptomic changes in lipid metabolism mediate SBM tolerance. Genotypic variations in DEGs may become biomarkers for improving early selection of fish tolerant to SMB or others plant-based diets.

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.

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

In this study, we used RNA-seq to analyze the muscle and liver tissues of black carps (Mylopharyngodon piceus) of different growth rates from the same batch to evaluate their growth traits. We have two groups; they are the black carp group with fast-growth rate and the slow-growth rate. A total of 23,132 genes were enriched in the Gene Ontology analysis, and 285 related pathways were found in the Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analysis. The KEGG pathway analysis showed significant differences in the expression of some genes involved in growth- and development-related metabolic pathways such as the FoxO signaling pathway, p53 signaling pathway, PI3K-Akt signaling pathway, apoptosis, TGF-β signaling pathway, and insulin signaling pathway. The numbers of differentially expressed genes in muscle and liver are 1913 and 1775. Nine of the differently expressed genes involved in the different growth traits and accuracy of the transcriptome data were validated using quantitative real-time PCR. We found that the expression levels of some growth-related genes were significantly higher in the fast-growth rate black carps than in the slow-growth rate black carps. The large number of transcriptome sequences obtained in this study has enriched the black carp gene resources, and the obtained differentially expressed genes and related pathway analysis provide valuable information for understanding the growth traits of the black carp.