Skip to main content
Download PDF
  • Research article
  • Open access
  • Published:

Transcriptome sequencing and metabolome analysis of food habits domestication from live prey fish to artificial diets in mandarin fish (Siniperca chuatsi)

BMC Genomics volume 22, Article number: 129 (2021) Cite this article

Abstract

Background

As economical traits, food habits domestication can reduce production cost in aquaculture. However, the molecular mechanism underlying food habits domestication has remained elusive. Mandarin fish (Siniperca chuatsi) only feed on live prey fish and refuse artificial diets. In the present study, we domesticated mandarin fish to feed on artificial diets. The two groups were obtained, the fish did not eat artificial diets or ate artificial diets during all of the three domestication processes, named Group W or X, respectively.

Results

Using transcriptome and metabolome analysis, we investigated the differentially expressed genes and metabolites between the two groups, and found three common pathways related to food habit domestication, including retinol metabolism, glycerolipid metabolism, and biosynthesis of unsaturated fatty acids pathways. Furthermore, the western blotting and bisulfite sequencing PCR analysis were performed. The gene expression of TFIIF and histone methyltransferase ezh1 were significantly increased and decreased in the fish of Group X, respectively. The total DNA methylation levels of TFIIF gene and tri-methylation of histone H3 at lysine 27 (H3K27me3) were significantly higher and lower in the fish of Group X, respectively.

Conclusion

It was speculated that mandarin fish which could feed on artificial diets, might be attributed to the lower expression of ezh1, resulting in the decreased level of H3K27me3 and increased level of DNA methylation of TFIIF gene. The high expression of TFIIF gene might up-regulate the expression of genes in retinol metabolism, glycerolipid metabolism and glycerophosphoric metabolism pathways. Our study indicated the relationship between the methylation of DNA and histone and food habits domestication, which might be a novel molecular mechanism of food habits domestication in animals.

Background

Food habits domestication can reduce production cost in animals. Mandarin fish, as an economic species, has very unique food preference. In the wild, as soon as they start to feed, they feed exclusively on live fry of other fish species [1]. Our previous study showed transcriptome determining of food preference (dead prey fish), and indicated that retinal photosensitivity, appetite control, circadian rhythm, learning and memory outputs might be involved in the food habit domestication of dead prey fish [2]. Compared to dead prey fish, the domestication of mandarin fish to accept artificial diets can provide more profitability. However, little studies investigate the molecular regulatory mechanisms of the domestication to accept artificial diets in mandarin fish.

Previous research showed that the hormones from central nervous systems play important roles in the food intake control, such as neuropeptide Y (NPY) and agouti-related protein (AgRP) [3, 4]. In giant panda, Tas1r1 pseudogenization reinforced the herbivorous life style because of the diminished attraction of returning to meat eating in the absence of Tas1r1 [5, 6]. The ion channels polycystic kidney disease 1-like 3 (PKD1L3) and PKD2L1 linked to sour taste, and the integral membrane protein CD36 is a putative “fat taste” receptor [7]. In a leaf-eating colobine monkey, metabolism genes pancreatic ribonuclease gene (RNASE1) was contributed to its food habits (leaves) [8]. However, little is known about the genetic and metabolic regulation on the food habits domestication of mandarin fish.

It has been noted that epigenetic status might be modified by environment and diets [9]. In mice, by feeding the diets with high levels of methyl donors (e.g. folic acid) to pregnant dams, it was possible to modify the expression of the agouti gene in the offspring with the high levels of DNA methylation [10, 11]. Histone modifications correlate with transcriptional activation and repression. The maternal undernutrition led to a decreased H3K27me3 level of the promotor region and increased expression of pomc gene in offspring mice [12]. Therefore, whether epigenetic regulation plays an important role in the food habits domestication is unknown.

In the present study, we domesticated the mandarin fish to accept artificial diets, and conducted the transcriptome sequencing and metabolome analysis to search the common pathways of transcriptome sequencing and metabolome. In addition, using western blotting and bisulfite sequencing PCR, we examined the methylation of histone and DNA, to investigate the molecular mechanism of food habits domestication in mandarin fish, which could promote the culture of mandarin fish with artificial diets.

Results

Pathway classification map of the differentially expressed genes based on transcriptome sequencing

The cDNA libraries were constructed from W and X groups of mandarin fish, and sequenced using the Illumina Hiseq2000 system. High quality reads were assembled. After removing the partial overlapping sequences, a total of 77,312 distinct sequences were obtained (All-Unigene, mean size: 1138 bp, N50: 2334 bp). In these unigene, 49.06% (37,927) were less than 500 bp, 50.94% (39,385) were longer than 500 bp, in which 34.38% (26,578) were longer than 1000 bp. We found 54 genes to be differential expressed among the two groups, 29 and 25 genes are up-regulated and down-regulated in mandarin fish of Group X, respectively. The metabolic pathway showed the most differential expressed genes (Fig. 1a and b), in which lipid metabolism, signal transduction and global overview maps showed 10, 6 and 13 genes to be differentially expressed, respectively (Fig. 1a). And the rich factor of steroid biosynthesis and glycerolipid metabolism is largest of all (Fig. 1b). The details of the differential expressed genes between the two groups were presented in Table 1. The sequencing data in this study have been deposited in the Sequence Read Archive (SRA) database (accession number: PRJNA613186).

Fig. 1
figure 1

a Pathway classification map of the differentially expressed genes. b Rich factor of the differentially expressed genes of different pathway based on transcriptome sequencing

Full size image
Table 1 Identification of differentially expressed gene based on transcriptome sequencing
Full size table

Analysis of differential metabolites of two groups

We analyzed the metabolic profiles of the two groups by LC-MS in positive (ESI+) and negative (ESI−) scan modes, and selected 9249 irons for subsequent analyses (4155 irons in ESI+ mode and 5094 irons in ESI− mode).

The normalized data were analyzed by PCA and PLS-DA with multivariate analysis. The PCA result showed the positive and negative ions from the different groups were in the two clusters, and were separated clearly by the first two components (Fig. 2a). PLS-DA result showed the clear separation of the two groups (Fig. 2b), suggesting the significant biochemical changes. The hierarchical clustering analysis (HCA) of the differential metabolites showed that Group X and W showed significant difference (Fig. 2c). The information of these metabolomic biomarkers was listed in Table 2.

Fig. 2
figure 2

a PCA scores scatter plot in positive ion (left) and in negative ion (right) scan modes for the two groups. b PLS-DA scores scatter plot in positive ion (left) and in negative ion (right) scan modes for the two groups. c The heat map of differential metabolites from the related pathways between the two groups in both positive and negative mode. Each line represents a differential metabolite and each cross represents a plasma sample group. Different colors represent different abundance intensity, and the higher abundance intensity shows a gradual increase from dark color to red color

Full size image
Table 2 Potential metabolomics biomarkers identified between the two groups
Full size table

To identify the metabolites, we used the freely accessible database of Kyoto Encyclopedia of Genes and Genomes (KEGG) to elucidate the putative function of the metabolites. 44 and 20 irons were identified by MS1 and MS2 level in positive mode respectively, and 24 and 11 irons in MS1 and MS2 level in negative mode respectively. The details of differential ions between the two groups were presented in Table 3.

Table 3 Identification of differential ions based on metabolome
Full size table

The common pathways of differential metabolites and genes

In retinol metabolism pathway, retinol, 9-cis-retinol and 11-cis-retinol metabolites were higher in mandarin fish of Group X than those of Group W, RDH (retinol dehydrogenase) gene expression was consistently higher in Group X (Fig. 3a). In glycerolipid metabolism pathway, triacylglycerol lipase gene expression was higher in mandarin fish of Group X, and glycerophosphoric metabolites was also higher in Group X (Fig. 3b). In biosynthesis of unsaturated fatty acids pathway, stearoyl-CoA gene expression and DPA (docosapentaenoic acid) metabolites were higher in fish of Group X than those in Group W (Fig. 3c).

Fig. 3
figure 3

Pathways of the differentially expressed genes and metabolites based on transcriptome and metabolome. a Retinol metabolism; b Glycerolipid metabolism; c Biosynthesis of unsaturated fatty

Full size image

TFIIF gene expression and DNA methylation

As is shown in Fig. 4a, General transcription factor IIF (TFIIF) gene expression was higher in the mandarin fish of Group X than that of Group W. We then analyzed the CpG islands at − 5000 bp upstream from the transcription initiation site (designated as 0) of TFIIF by methylation analysis software. As shown in Fig. 4b, one CpG islands containing 9 CpG sites existed in − 3619 to − 3574 bp of TFIIF gene. The total DNA methylation level was significantly higher in the fish of Group X than that of Group W (Table 4).

Fig. 4
figure 4

TFIIF gene expression and DNA methylation. a TFIIF gene expression. b Illustration of the region of CpG islands sites, which includes 9 CpG sites, DNA methylation patterns of the two groups (X and W) analyzed by BSP. Each line represents one individual bacterial clone, and each circle represents one single CpG dinucleotide. Open circles show unmethylated CpGs and black circles show methylated CpGs

Full size image
Table 4 Methylation of each CpG island of TFIIF gene between the two groups
Full size table

Ezh1 gene expression and histone methylation

The mRNA expression of histone methyltransferase ezh1 gene was lower in the mandarin fish of Group X (Fig. 5a). As histone methyltransferase Ezh1 could methylate ‘Lys-27’ of histone H3, we analyzed the H3K27me3 levels of the two groups. The results showed that H3K27me3 level was also lower in the mandarin fish of Group X than that of Group W (Fig. 5b).

Fig. 5
figure 5

a Validation of ezh1 mRNA expression. b The H3K27me3 protein level of between Group X and W. Data are mean ± SEM (n = 6), significant difference is marked with an asterisk (P < 0.05)

Full size image

Discussion

In rearing conditions, mandarin fish accept only live prey fish, refusing dead prey fish or artificial diets [13]. Although pervious research showed the methods of mandarin fish domestication [14], little is known about the mechanism of food habits domestication. In the present study, we domesticated the mandarin fish to feed on artificial diets, and found a part of mandarin fish could accept artificial diets easily (Group X), but another part could not accept completely (Group W). To uncover the molecular mechanism why mandarin fish refuses artificial diets, we conducted the transcriptome sequencing and metabolome analysis. The results showed that the differentially expressed gene between the two groups were enriched in metabolism, in which the global and overview maps and lipid metabolism were the most enriched. And the rich factors of steroid biosynthesis and glycerolipid metabolism were the highest. The metabolome results showed that the pathways with different metabolites were mostly enriched in the metabolic pathways, which were consistent with transcriptome sequencing results. Previous research has shown that the most important pathways related to the domestication of dead prey fish in mandarin fish included the retinal photosensitivity, circadian rhythm, appetite control, learning and memory pathway [2]. Our results showed that metabolism, especially lipid metabolism, might contribute to the domestication of artificial diets, which was different from the domestication of dead prey fish, as the different constituents between dead prey fish and artificial diets.

We then analyzed the pathways in which the differential genes or metabolites were involved, the common pathways which showed the most enriched differential genes and metabolites, were retinol metabolism, glycerolipid metabolism and biosynthesis of unsaturated fatty acids. For retinol metabolism, retinol, 9-cis-retinol and 11-cis-retinol metabolites were higher in the Group X, consistently the RDH (retinoldehydrogenase) gene expression was higher in the Group X, suggesting a better visual acuity in the mandarin fish which could be easy to accept artificial diets. Animals make food choices on the basis of the nutritional, physiological, environmental, and sociocultural factors [7], sensory system is of significance to food choices. It is critical for mandarin fish to catch prey fish though the perception of shape and motion with well-developed scotopic vision [13]. Salmo spp. shows the same motion and shape of food, they have high visual acuity, thus can feed swiftly by darting, the offered food pellet can be captured immediately before it falls down to the bottom of the tank [15,16,17]. Because of the low visual acuity and feeding only by stalking, mandarin fish can not recognize the prey before the time when food pellet fall to the bottom of tank, thus makes it difficult to feed mandarin fish with artificial diets [13]. The retinol metabolism dysfunction might be contributed to the lower visual ability in the mandarin fish which refused artificial diets.

In glycerolipid metabolism pathway, the gene expression of triacylglycerol lipase was higher in the mandarin fish of Group X, and the glycerophosphoric acid metabolite was also higher in the Group X. In the biosynthesis of unsaturated fatty acids pathway, stearoyl-CoA gene expression and docosapentaenoic acid (DPA) metabolite were higher in the Group X. These results suggested that mandarin fish which could accept artificial diets well, might be attributed to the better capacity of glycerolipid metabolism and unsaturated fatty acids biosynthesis. Live food diets (such as zooplankton) and dry formulated diets have different fat levels and influences in European grayling [18]. Artificial diets might have more fat and energy than live prey fish, suggesting that mandarin fish which accept artificial diets could make good use of fat, while the fish which refuse artificial diets could not.

To elucidate the regulatory mechanism of up-regulated gene expression in the mandarin fish of Group X, we analyzed the differentially expressed genes based on transcriptome sequencing. The results showed TFIIF gene expression was significantly increased in the Group X. TFIIF communicates with a number of factors to regulate gene transcription. It has been reported that TFIIF directly binds to basal factors of TFIID, TFIIE and TFIIB [19]. TFIIF has been shown to be necessary for most, if not all, preinitiation complex formation and gene transcription [20, 21]. It suggested that in the mandarin fish which accepted artificial diets, the up-regulated transcription of genes, involved in retinol metabolism, glycerolipid metabolism and biosynthesis of unsaturated fatty acids, might be contributed to the increased TFIIF expression.

To uncover why TFIIF was upregulated in the mandarin fish of Group X, the expressions of histone methyltransferases were analyzed based on transcriptome. The expression of histone-lysine N-methyltransferase ezh1 was significantly decreased in the Group X. Histone methyltransferases EZH1 and EZH2 catalyze the tri-methylation of H3K27, which serves as an epigenetic signal for chromatin condensation and transcriptional repression [22]. In mice, Ezh1 was required for neonatal heart regeneration after myocardial infarction and overexpression of Ezh1 promoted heart regeneration by upregulating cardiac muscle growth genes [23]. Furthermore, we observed the protein level of tri-methylation of histone H3 at lysine 27 was lower in the Group X, suggesting an active function of gene expression. The decreased abundance of histone H3K27me3 was also found in FOXO1 (forkhead box protein O1) in HFD (high fat diets) fed rats, which persisted even after 8 weeks of diet reversal [24]. In addition, the total DNA methylation level of TFIIF was significantly higher in the mandarin fish of Group X than those of Group W. The mRNA level of TFIIF was higher in fish of Group X, our results showed a positive effect of DNA methylation on gene expression. In soybean, the genome-wide methylation profiles showed that hyper-methylated genes had higher gene expression [25], in ventricular septal defect patients, genome-wide DNA methylation data showed 12 hypermethylated genes had a higher gene expression [26]. In Arabidopsis, upon loss of CpG methylation, there was target-specific enrichment of H3K27me3 in heterochromatin that correlated with transcriptional reactivation, it is suggested that there was an antagonistic effect between CpG methylation and H3K27me3 [27]. It is hypothesized that the lower ezh1 expression in the mandarin fish of Group X, could be contributed to the decreased methylation at ‘Lys-27’ of histone H3, and then up-regulating the expression and methylation of TFIIF gene.

Conclusions

Our research indicated the individual differences of acceptance on artificial diets in mandarin fish and the potential molecular mechanism. The mandarin fish which feed on artificial diets, could down-regulate the ezh1 expression, repressing the tri-methylation level of histone H3 at lysine 27, and then resulting in the increased DNA methylation and mRNA expression of TFIIF gene. TFIIF as an important transcription factor, might regulate the expression of genes involved in retinol metabolism, glycerolipid metabolism and glycerophosphoric metabolism, and modify the acceptance on artificial diets of mandarin fish. These results suggested the potential effect of histone methylation on food habits domestication in mandarin fish.

Methods

Fish domestication and sampling

Mandarin fish (69.9 ± 10.2 g) were obtained from Chinese Perch Research Center of Huazhong Agricultural University (Wuhan, Hubei Province, China) and maintained in the aquarium (12 tanks, 50 fish per tank) at constant temperature (25 ± 0.5 °C). Mandarin fish were and domesticated and fed with artificial diets (Table 5) and divided into two groups: fish did not eat artificial diets and fish ate artificial diets. The fish did not eat artificial diets during the first domestication process was then fed with live fish prey for three days, starved for two days and fed with artificial diets for one day, and then we selected the fish did not eat artificial diets during the second domestication process and repeated the domestication process for one more time. The fish ate artificial diets during the first domestication process was fed with live fish prey for one days and fed with artificial diets for three days, then we selected the fish ate artificial diets during the second domestication process and repeated the domestication process for one more time. Finally, the two groups were obtained, the fish did not eat artificial diets or ate artificial diets during all of the three domestication processes, named Group W (n = 56) or X (n = 24), respectively. Six fish were used for real-time quantitative PCR. Six fish were used for western blotting. Ten fish were used for metabolome, and three fish were used for transcriptome sequencing.

Table 5 Composition of artificial diets
Full size table

The experimental fish were anesthetized with MS-222 (200 mg/L) (Redmond, WA, USA) and sacrificed by decapitation according to the ethical guidelines of Huazhong Agricultural University. Immediately after the surgical resection, the liver tissue was frozen in liquid nitrogen and stored at − 80 °C until used. The blood was drawn from the tail vein. The plasma was obtained from whole blood sample separated with 4000 r/centrifuge for 10 min. The animal protocol was approved by the Institutional Animal Care and Use Ethics Committee of Huazhong Agricultural University (Wuhan, China) (HZAUFI-2017-015).

RNA isolation and reverse transcription

Total RNA was extracted using Trizol reagent, and one microgram of total RNA was synthesized to complementary DNA (cDNA) by Revert Aid™ Reverse Transcriptase (TaKaRa, Tokyo, Japan). The protocols were following the manufacturer’s instructions.

Transcriptome sequencing

For transcriptome analysis, the equal amount of total RNA was used to construct the libraries for each group (n = 3) using MGIEasy RNA kit (BGI, Wuhan, China). The paired-end cDNA libraries were constructed and sequenced using BGISEQ-500 system (BGI, Wuhan, China). The processes of image deconvolution, base calling, unigene assembly, annotation and expression level estimation were carried out as described by You et al. [28]. The differentially expressed genes were identified as DEGseq method described before [29], and the significance of gene expression difference was judged with Fold Change ≥2.00 and FDR (False Discovery Rate) ≤ 0.001 as the threshold. The analysis of GO function and KEGG pathway of differentially expressed genes were carried out.

Metabolome

Serum samples (40 μl, including QC samples) were added to new Eppendorf tubes with ice-cold methanol (120 μl), vortex mixed for 1 min, placed in holding for 30 min at − 20 °C, and centrifuged at 4000 g for 20 min at 4 °C. 25 μl of supernatant and 225 μl 50% methanol were mixed. Then 20 μl of mixture from each sample were mixed as quality control samples, 60 μl of mixture was conducted as samples. All samples were stored at − 80 °C (ten biological replicates for each group).

All samples were acquired by the LC-MS system followed machine orders. Chromatographic separations were performed using an ultra performance liquid chromatography (UPLC) system (Waters, USA) and a high-resolution tandem mass spectrometer SYNAPT G2 XS QTOF (Waters, USA) was used to detect metabolites as the methods described by Huang et al. [30]. The mass spectrometry data were acquired in Centroid MSE mode. Statistical analysis was performed as previous [31]. Putative metabolites were first derived by searching the exact molecular mass data from redundant m/z peaks against the online HMDB (http://www.hmdb.ca/), METLIN (http://metlin.scripps.edu/) and KEGG (www.genome.jp/kegg/) databases.

Real-time quantitative PCR

We searched the cDNA sequences from transcriptome data of mandarin fish, designed the primers with Primer 5.0 software (Table 6). The potential housekeeping gene according to the literature [32], beta-actin, b2m, rpl13a, and hmbs were examined, and Rpl13a gene was more stable and selected as the internal control. Real-time quantitative PCR was performed with MyiQ™ 2 Two-Color Real-Time PCR Detection System (Bio-Rad, Hercules, USA) as the methods described by Liang et al. [33]. The target gene expression relative to rpl13a expression were calculated by the optimized comparative Ct (2-ΔΔCt) value method [34]. Data with six biological replicates and three technical replicates were presented as mean ± S.E.M.

Table 6 Nucleotide sequences of the primers
Full size table

DNA methylation analysis

TIANamp Genomic DNA Kit (Tiangen, Beijing, China) was used for Genomic DNA extraction, and then the DNA was treated with sodium bisulfite by EZ DNA Methylation Kit (Zymo Research, USA). The BSP primers were designed by the online MethPrimer software 14 and Primer 5.0 (Table 6). As the methods described by Cai et al. [35], PCR products were subcloned and sequenced. Six samples from the mandarin fish of Group W or X were analyzed with five technical replicates.

Western blotting

The proteins from liver tissue were separated on 10% SDS-PAGE gel and transferred onto PVDF membrane. Tri-Methyl-Histone H3 (lys27) (C36B11) Rabbit mAb was obtained from Cell Signaling Technology (Danvers, MA). The protein level H3K27me3 were detected by western blotting with the antibody (1:1000–1:4000) according to manufacturer’s instructions. Blots were probed by second antibody labeled with IR-Dye 680 or 800 cw (1:2000–1:4000, Licor, Lincoln, NE, USA) and the membranes were visualized and quantified as the methods described by You et al. [28]. (six biological replicates for each group).

Statistical analysis

SPSS 19.0 software was used for the statistical analyses. The normality and homogeneity of variances were tested using the Shapiro-Wilk’s test and Levene’s test, respectively. The significant differences were found using one-way analysis of variance (ANOVA), followed by Duncan’s multiple range tests and Fisher’s least significant difference post hoc test. Differences with P < 0.05 were considered to be significant.

Availability of data and materials

The sequencing data in this study have been deposited in the Sequence Read Archive (SRA) database (accession number: PRJNA613186).

Abbreviations

EZH1:

Histone H3-lysine27 N-trimethyltransferase EZH1

H3K27me3:

Histone H3-lysine27 tri-methylation

TFIIF:

General transcription factor IIF

NPY:

Neuropeptide Y

AgRP:

Agouti-related protein

Tas1r1:

Taste receptor type 1 member 1

PKD1L3:

Polycystic kidney disease 1-like 3

PKD2L1:

Polycystic kidney disease 2-like 1

RNASE1:

Pancreatic ribonuclease

B2m:

Beta-2-microglobulin

Rpl13a:

Ribosomal protein L13a

Hmbs:

Hydroxymethylbilane synthase

QC:

Quality control

PCA:

Principal component analysis

LC-MS:

Liquid chromatography mass spectrometry

KEGG:

Kyoto Encyclopedia of Genes and Genomes

RDH:

Retinol dehydrogenase

DPA:

Docosapentaenoic acid

GO:

Gene Ontology

FOXO1:

Forkhead box protein O1

HFD:

High fat diets

UPLC:

Ultra performance liquid chromatography

FDR:

False Discovery Rate

ANOVA:

One-way analysis of variance

References

  1. Liang XF, Lin X, Li S, Liu JK. Impact of environmental and innate factors on the food habit of Chinese perch Siniperca chuatsi (Basilewsky) (Percichthyidae). Aquac Res. 2008;39:150–7.

    Article  Google Scholar 

  2. He S, Liang XF, Sun J, Li L, Yu Y, HuangW Q. CM, Cao L, Bai XL, Tao YX. Insights into food preference in hybrid f1 of siniperca chuatsi (♀) × siniperca scherzeri (♂) mandarin fish through transcriptome analysis. BMC Genomics. 2013;14.

  3. Morton GJ, Cummings DE, Baskin DG, Barsh GS, Schwartz MW. Central nervous system control of food intake and body weight. Nature. 2006;443:289–95.

    Article  CAS  PubMed  Google Scholar 

  4. Cone RD. Anatomy and regulation of the central melanocortin system. Nat Neurosci. 2005;8:571–8.

    Article  CAS  PubMed  Google Scholar 

  5. Zhao H, Yang JR, Xu H, Zhang J. Pseudogenization of the umami taste receptor gene Tas1r1 in the giant panda coincided with its dietary switch to bamboo. Mol Biol Evol. 2010;27:2669–73.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  6. Li R, Fan W, Tian G, Zhu H, He L, Cai J, Huang Q, Cai Q, Li B, Bai Y, et al. The sequence and de novo assembly of the giant panda genome. Nature. 2010;463:311–7.

    Article  CAS  PubMed  Google Scholar 

  7. Garcia-Bailo B, Toguri C, Eny KM, El-Sohemy A. Genetic variation in taste and its influence on food selection. OMICS. 2009;13:69–80.

    Article  CAS  PubMed  Google Scholar 

  8. Zhang J, Zhang YP, Rosenberg HF. Adaptive evolution of a duplicated pancreatic ribonuclease gene in a leaf-eating monkey. Nat Genet. 2002;30:411–5.

    Article  CAS  PubMed  Google Scholar 

  9. Johnson IT, Belshaw NJ. Environment, diet and cpg island methylation: epigenetic signals in gastrointestinal neoplasia. Food Chem Toxicol Int J Published Br Industrial Biol Res Assoc. 2008;46(4):0–1359.

    Google Scholar 

  10. Wolff GL, Kodell RL, Moore SR, Cooney CA. Maternal epigenetics and methyl supplements affect agouti gene expression in Avy/a mice. FASEB J. 1998;12:949–57.

    Article  CAS  PubMed  Google Scholar 

  11. Cooney CA, Dave AA, Wolff GL. Maternal methyl supplements in mice affect epigenetic variation and DNA methylation of offspring. J Nutr. 2002;132:2393S–400S.

    Article  CAS  PubMed  Google Scholar 

  12. Begum G, Stevens A, Smith EB, Connor K, Challis JRG, Bloomfield F, White A. Epigenetic changes in fetal hypothalamic energy regulating pathways are associated with maternal undernutrition and twinning. FASEB J. 2012;26(4):1694–703.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  13. Liang XF, Liu JK, Huang BY. The role of sense organs in the feeding behaviour of Chinese perch. J Fish Biol. 1998;52:1058–67.

    Article  Google Scholar 

  14. Liang X, Oku H, Ogata H, Liu J, He X. Weaning Chinese perch Siniperca chuatsi (Basilewsky) onto artificial diets based upon its specific sensory modality infeeding. Aquac Res. 2001;32:76–82.

    Article  Google Scholar 

  15. Irvine JR, Northcote TG. Selection by young rainbow trout (Salmo gairdneri) in simulated stream environments for live and dead prey of different sizes. Can J Fish Aquat Sci. 1983;40:1745–9.

    Article  Google Scholar 

  16. Stradmeyer L, Metcalfe NB, Thorpe JE. Effect of food pellet shape and texture on the feeding response of juvenile Atlantic salmon. Aquaculture. 1988;73:217–28.

    Article  Google Scholar 

  17. Stradmeyer L. A behavioural method to test feeding responses of fish to pelleted diets. Aquaculture. 1989;79:303–10.

    Article  Google Scholar 

  18. Carlstein M. Natural food and artificial, dry starter diets: effects on growth and survival in intensively reared European grayling. Aquac Int. 1993;1(2):112–23.

    Article  Google Scholar 

  19. McEwan IJ, Gustafsson J. Interaction of the human androgen receptor transactivation function with the general transcription factor TFIIF. Proc Natl Acad Sci U S A. 1997;94:8485–90.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  20. Tan S, Conaway RC, Conaway JW. Dissection of transcription factor TFIIF functional domains required for initiation and elongation. Proc Natl Acad Sci U S A. 1995;92:6042–6.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  21. Killeen MT, Greenblatt JF. The general transcription factor RAP30 binds to RNA polymerase II and prevents it from binding nonspecifically to DNA. Mol Cell Biol. 1992;12:30–7.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  22. Lui JC, Garrison P, Nguyen Q, Ad M, Keembiyehetty C, Chen W, et al. Ezh1 and ezh2 promote skeletal growth by repressing inhibitors of chondrocyte proliferation and hypertrophy. Nat Commun. 2016;7:13685.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  23. Ai S, Yu X, Li Y, Peng Y, Li C, Yue Y, et al. Divergent requirements for ezh1 in heart development versus regeneration. Circ Res. 2017;121:106–12.

  24. Kumar S, Pamulapati H, Tikoo K. Fatty acid induced metabolic memory involves alterations in renal histone h3k36me2 and h3k27me3. Mol Cell Endocrinol. 2015;422:233–42.

    Article  PubMed  Google Scholar 

  25. Yanwei L, Xianlong D, Xuan W, Tingting H, Hao Z, Longshu Y, et al. Genome-wide comparative analysis of dna methylation between soybean cytoplasmic male-sterile line njcms5a and its maintainer njcms5b. BMC Genomics. 2017;18(1):596.

    Article  Google Scholar 

  26. Marcel G, Cornelia D, Huanhuan C, Ilona D, Kerstin S, Sophia S, et al. Comparative dna methylation and gene expression analysis identifies novel genes for structural congenital heart diseases. Cardiovasc Res. 2016;112:464–77.

  27. Mathieu O, Probst AV, Paszkowski J. Distinct regulation of histone h3 methylation at lysines 27 and 9 by cpg methylation in arabidopsis. EMBO J. 2005;24(15):2783–91.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  28. You JJ, Ren P, He S, Liang XF, Zhang YP. Histone methylation of h3k4 involved in the anorexia of carnivorous mandarin fish (siniperca chuatsi) after feeding on a carbohydrate-rich diet. Front Endocrinol. 2020;11:323.

  29. Wang L, Feng Z, Wang X, Wang X, Zhang X. DEGseq, an R package for identifying differentially expressed genes from RNA-seq data. Bioinformatics. 2010;26(1):136–8.

    Article  PubMed  Google Scholar 

  30. Huang JF, Shen ZY, Mao QL, Zhang XM, Zhang B, Wu JS, et al. Systematic analysis of bottlenecks in a multibranched and multilevel regulated pathway: the molecular fundamentals of l-methionine biosynthesis in Escherichia coli. ACS Synth Biol. 2018;7:2577–89.

  31. Zhou LF, Zhao BW, Guan NN, et al. Plasma metabolomics profiling for fish maturation in blunt snout bream. Metabolomics. 2017;13(4):40.

    Article  Google Scholar 

  32. Vandesompele J, Preter KD, Pattyn F, Poppe B, Roy NV, Paepe AD. Accurate normalization of real-time quantitative RT-PCR data by geometric averaging of multiple internal control genes. Genome Biol. 2018;3:1–12.

    Google Scholar 

  33. Liang H, He S, Liang XF, Lu HL, Chen K. Feeding habit transition induced by social learning through CAMKII signaling in chinese perch (siniperca chuatsi). Aquaculture. 2020;533:736211.

  34. Livak KJ, Schmittgen TD. Analysis of relative gene expression data using real-time quantitative PCR and the 2−ΔΔCT method. Methods. 2001;25:402–8.

    Article  CAS  PubMed  Google Scholar 

  35. Cai WJ, He S, Liang XF, Yuan XC. DNA methylation of T1R1 gene in the vegetarian adaptation of grass carp ctenopharyngodon idella. Sci Rep. 2018;8:6934.

    Article  PubMed  PubMed Central  Google Scholar 

Download references

Acknowledgments

Not applicable.

Funding

This work was financially supported by the China Agriculture Research System (CARS-46), the National Natural Science Foundation of China (31772822), and the Open Project of Anhui Province Key Laboratory of Aquaculture & Stock Enhancement (AHSC202002).

Author information

Authors and Affiliations

  1. College of Fisheries, Chinese Perch Research Center, Huazhong Agricultural University, 1 Shizishan Street, Wuhan, 430070, Hubei, China

    Shan He, Jun-Jie You, Xu-Fang Liang, Zhi-Lu Zhang & Yan-Peng Zhang

  2. Innovation Base for Chinese Perch Breeding, Key Lab of Freshwater Animal Breeding, Ministry of Agriculture, Wuhan, 430070, China

    Shan He, Jun-Jie You, Xu-Fang Liang, Zhi-Lu Zhang & Yan-Peng Zhang

  3. Anhui Province Key Laboratory of Aquaculture & Stock Enhancement, Fishery Institute of Anhui Academy of Agricultural Sciences, Hefei, 230031, China

    Shan He & Jun-Jie You

Authors
  1. Shan He
    View author publications

    You can also search for this author in PubMed Google Scholar

  2. Jun-Jie You
    View author publications

    You can also search for this author in PubMed Google Scholar

  3. Xu-Fang Liang
    View author publications

    You can also search for this author in PubMed Google Scholar

  4. Zhi-Lu Zhang
    View author publications

    You can also search for this author in PubMed Google Scholar

  5. Yan-Peng Zhang
    View author publications

    You can also search for this author in PubMed Google Scholar

Contributions

S.H., J.J.Y., Z.L.Z. and Y.P.Z. contributed to the sample preparation and examination. S.H. and L.L. performed the assembly and bioinformatical analysis. S.H. and J.J.Y. performed the data analysis. S.H. and X.F.L. gave technical advice and contributed to the study design. S.H. and J.J.Y. wrote the paper. All authors read and approved the final manuscript.

Corresponding author

Correspondence to Xu-Fang Liang.

Ethics declarations

Ethics approval and consent to participate

The animal protocol was approved by the Institutional Animal Care and Use Ethics Committee of Huazhong Agricultural University (Wuhan, China) (HZAUFI-2017-015).

Consent for publication

Not applicable.

Competing interests

The authors declare that they have no conflict of interest.

Additional information

Publisher’s Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated in a credit line to the data.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

He, S., You, JJ., Liang, XF. et al. Transcriptome sequencing and metabolome analysis of food habits domestication from live prey fish to artificial diets in mandarin fish (Siniperca chuatsi). BMC Genomics 22, 129 (2021). https://doi.org/10.1186/s12864-021-07403-w

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s12864-021-07403-w

Keywords

BMC Genomics

ISSN: 1471-2164

Contact us