Is the Nutritional Value of Fish Fillet Related to Fish Maturation or Fish Age? Integrated Analysis of Transcriptomics and Metabolomics in Blunt Snout Bream (Megalobrama amblycephala).

BACKGROUND/AIMS: Fish is a protein-rich food and is increasingly favored by consumers. It has been well recognized that the flesh composition of fish is closely related to its maturation and growth stage, but few studies have explored these differences. Additionally, hormone residues in fish after artificial induction of reproduction also attract consumer concern. In this study, we attempt to address these concerns by using a combination of transcriptomics and metabolomics analyses to identify key regulated pathways, genes, and metabolites that may affect the flesh nutrition of one typical aquaculture species in China, blunt snout bream (Megalobrama amblycephala). METHODS: The four groups of fish were used for transcriptomics and metabolomics analyses, including one-year-old immature (group I), two-year-old immature (group II), two-year-old mature (group III) and successfully spawned (group IV) female M. amblycephala after artificial induction of reproduction. RESULTS: We identified a total of 1460 differential compounds and 1107 differentially expressed unigenes in muscle among the different groups. Differential metabolites related to fish age (group II vs group I, group III vs group I) were largely enriched in "Glycerophospholipid metabolism", "Linoleic acid metabolism", "α-Linolenic acid metabolism", and "Biosynthesis of unsaturated fatty acids". Between these two pairwise comparisons, metabolites that are beneficial to human health, such as docosapentaenoic acid, α-Linolenic acid, eicosapentaenoic acid, and docosahexaenoic acid were found to be significantly decreased in two-year-old (group II, group III) compared with one-year-old (group I) M. amblycephala. Only one differential metabolite related to fish maturation, a triglyceride, was detected between groups III and II. Transcriptomics data showed that differently expressed genes (between group III vs group II, group III vs group I) related to maturation were highly enriched in "Cell adhesion molecules (CAMs)", "Sphingolipid metabolism" and "Phagosome". DEGs (between group II vs group I, group III vs group I) relating to fish age were enriched in the "cGMP-PKG signaling pathway", "FoxO signaling pathway", and "AMPK signaling pathway". The gene-metabolite interaction network showed pivotal genes, including fumarate hydratase and GNPAT, which played a major role in the regulation of glycerphospholipid metabolism. The nutritional components were also measured, which verified the metabolomics results. Moreover, the metabolomics results showed that after 24 hours of artificial hormone injection, the drug was completely metabolized. CONCLUSION: Integrated analysis demonstrated that the nutrition value of fish fillet was much more related to fish age compared with maturation status in M. amblycephala females.

Characterization and spatiotemporal expression analysis of nine bone morphogenetic protein family genes during intermuscular bone development in blunt snout bream.

Intermuscular bones (IBs) only exist in the myosepta of lower teleosts and its molecular mechanism remains to be clarified. Bone morphogenetic proteins (BMPs) have been demonstrated to be involved in various physiological processes, including bone and cartilage formation. In this study, we firstly obtained and characterized nine bmp genes for Megalobrama amblycephala, which belongs to Cyprinidae and have a certain amount of IBs. Sequence alignment and phylogenetic analysis both documented that the mature proteins of M. amblycephala bmp genes were highly conserved with other corresponding homologs, respectively, indicating that the function of each bmp gene has been conserved throughout evolution. As a step to characterize potential involvement of bmp genes in IB formation and development, spatiotemporal expressions of nine bmp genes (bmp2a, bmp2b, bmp3, bmp4, bmp5, bmp7b, bmp8a, bmp14 and bmp16) were investigated during the key development stages of IBs. During the ossification process from stage I (the IBs haven't emerged) to stage IV (all of the IBs ossified in the tail with the mature morphology), the expression profiles revealed that bmp16 was the most abundant transcript while bmp4 had the lowest abundance. The mRNA levels of bmp3, bmp4, bmp5 and bmp8a increased significantly at stage II, suggesting their roles in stimulating IB formation. The expression of bmp7b reached the highest level at stage III (the rapid period of IB development), suggesting potential involvement of bmp7b in promoting osteoblast differentiation. With the exception of bmp7b and bmp16, most bmp genes appeared a significant increase at IB maturation phase (stage IV), which means that they may play important roles in maintenance of IB morphogenesis. Spatial tissue distribution of bmp genes showed that most bmp genes were observed at the highest level in developing IBs at one year old fish. Spatiotemporal expression patterns suggest the potential key roles of these bmp genes in IBs formation and maintenance in fish, being as possible promoters or inhibitors.

Dynamic mRNA and miRNA expression analysis in response to intermuscular bone development of blunt snout bream (Megalobrama amblycephala).

Intermuscular bone (IB), which occurs only in the myosepta of lower teleosts, is attracting more attention because they are difficult to remove and make the fish unpleasant to eat. By gaining a better understanding of the genetic regulation of IB development, an integrated analysis of miRNAs and mRNAs expression profiling was performed on Megalobrama amblycephala. Four key development stages were selected for transcriptome and small RNA sequencing. A number of significantly differentially expressed miRNAs/genes associated with bone formation and differentiation were identified and the functional characteristics of these miRNAs/genes were revealed by GO function and KEGG pathway analysis. These were involved in TGF-ß, ERK and osteoclast differentiation pathways known in the literature to affect bone formation and differentiation. MiRNA-mRNA interaction pairs were detected from comparison of expression between different stages. The function annotation results also showed that many miRNA-mRNA interaction pairs were likely to be involved in regulating bone development and differentiation. A negative regulation effect of two miRNAs was verified through dual luciferase reporter assay. As a unique public resource for gene expression and regulation during the IB development, this study is expected to provide forwards ideas and resources for further biological researches to understand the IBs' development.

Mitochondrial Genome Variation after Hybridization and Differences in the First and Second Generation Hybrids of Bream Fishes.

Hybridization plays an important role in fish breeding. Bream fishes contribute a lot to aquaculture in China due to their economically valuable characteristics and the present study included five bream species, Megalobrama amblycephala, Megalobrama skolkovii, Megalobrama pellegrini, Megalobrama terminalis and Parabramis pekinensis. As maternal inheritance of mitochondrial genome (mitogenome) involves species specific regulation, we aimed to investigate in which way the inheritance of mitogenome is affected by hybridization in these fish species. With complete mitogenomes of 7 hybrid groups of bream species being firstly reported in the present study, a comparative analysis of 17 mitogenomes was conducted, including representatives of these 5 bream species, 6 first generation hybrids and 6 second generation hybrids. The results showed that these 17 mitogenomes shared the same gene arrangement, and had similar gene size and base composition. According to the phylogenetic analyses, all mitogenomes of the hybrids were consistent with a maternal inheritance. However, a certain number of variable sites were detected in all F1 hybrid groups compared to their female parents, especially in the group of M. terminalis (♀) × M. amblycephala (♂) (MT×MA), with a total of 86 variable sites between MT×MA and its female parent. Among the mitogenomes genes, the protein-coding gene nd5 displayed the highest variability. The number of variation sites was found to be related to phylogenetic relationship of the parents: the closer they are, the lower amount of variation sites their hybrids have. The second generation hybrids showed less mitogenome variation than that of first generation hybrids. The non-synonymous and synonymous substitution rates (dN/dS) were calculated between all the hybrids with their own female parents and the results indicated that most PCGs were under negative selection.

Identification of MicroRNA for Intermuscular Bone Development in Blunt Snout Bream (Megalobrama amblycephala).

Intermuscular bone (IB), which occurs only in the myosepta of the lower teleosts, is attracting more attention of researchers due to its particular development and lack of genetic information. MicroRNAs (miRNAs) are emerging as important regulators for biological processes. In the present study, miRNAs from IBs and connective tissue (CT; encircled IBs) from six-month-old Megalobrama amblycephala were characterized and compared. The results revealed the sequences and expression levels of 218 known miRNA genes (belonging to 97 families). Of these miRNAs, 44 known microRNA sequences exhibited significant expression differences between the two libraries, with 24 and 20 differentially-expressed miRNAs exhibiting higher expression in the CT and IBs libraries, respectively. The expressions of 11 miRNAs were selected to validate in nine tissues. Among the high-ranked predicted gene targets, differentiation, cell cycle, metabolism, signal transduction and transcriptional regulation were implicated. The pathway analysis of differentially-expressed miRNAs indicated that they were abundantly involved in regulating the development and differentiation of IBs and CT. This study characterized the miRNA for IBs of teleosts for the first time, which provides an opportunity for further understanding of miRNA function in the regulation of IB development.