- Research article
- Open Access
Global gene expression patterns of grass carp following compensatory growth
BMC Genomics volume 16, Article number: 184 (2015)
Compensatory growth is accelerated compared with normal growth and occurs when growth-limiting conditions are overcome. Most animals, especially fish, are capable of compensatory growth, but the mechanisms remain unclear. Further investigation of the mechanism of compensatory growth in fish is needed to improve feeding efficiency, reduce cost, and explore growth-related genes.
In the study, grass carp, an important farmed fish in China, were subjected to a compensatory growth experiment followed by transcriptome analysis by RNA-sequencing. Samples of fish from starved and re-feeding conditions were compared with the control. Under starved conditions, 4061 and 1988 differentially expressed genes (DEGs) were detected in muscle and liver tissue when compared the experimental group with control group, respectively. After re-feeding, 349 and 247 DEGs were identified in muscle and liver when the two groups were compared. Moreover, when samples from experimental group in starved and re-feeding conditions were compared, 4903 and 2444 DEGs were found in muscle and liver. Most of these DEGs were involved in metabolic processes, or encoded enzymes or proteins with catalytic activity or binding functions, or involved in metabolic and biosynthetic pathways. A number of the more significant DEGs were subjected to further analysis. Under fasting conditions, many up-regulated genes were associated with protein ubiquitination or degradation, whereas many down-regulated genes were involved in the metabolism of glucose and fatty acids. Under re-feeding conditions, genes participating in muscle synthesis and fatty acid metabolism were up-regulated significantly, and genes related to protein ubiquitination or degradation were down-regulated. Moreover, Several DEGs were random selected for confirmation by real-time quantitative PCR.
Global gene expression patterns of grass carp during compensatory growth were determined. To our knowledge, this is a first reported for a teleost fish. The results will enhance our understanding of the mechanism of compensatory growth in teleost fish.
Growth is an important trait in fish farming and is one of the primary targets in breeding programs aimed at improving productivity . Growth of fish is governed by multiple genes and is also influenced by various environmental factors . Several important growth-related genes have been identified in various domesticated vertebrates, including growth hormone (GH), growth hormone receptor (GHR), insulin-like growth factor (IGF) I and II, growth hormone-releasing hormone (GHRH), leptins, growth hormone inhibiting hormone (GHIH), myostatin (MSTN), myogenic regulatory factors (MRFs), and many others [3,4]. However, reports on other growth-related genes in fish are relatively limited. Quantitative trait loci (QTLs) have been successfully used to locate growth-associated genes to particular regions of the fish genome, but identifying individual genes has proved difficult [5-8]. Other methods such as compensatory growth experiments proved fruitful for understanding regulatory mechanisms connected with growth in fish [9-12].
Compensatory growth is a period of accelerated growth that follows growth-limiting conditions once non-limiting conditions are renewed . Characteristic features of compensatory growth include increased food-intake, accelerated mitosis and enhanced rate of food utilization. Compensatory growth was first reported a century ago and has been widely studied in vertebrates . Most animals, especially fish, are capable of compensatory growth . However, the regulatory mechanisms and global gene expression patterns of compensatory growth in fish remain poorly understood. Expanding knowledge in this area is important for identifying growth-associated genes, increasing the efficiency of feeding and reducing the cost of fish farming.
Transcriptome or RNA-sequencing (RNA-seq) is a technology based on next-generation sequencing that is being successfully applied to transcriptome analysis . RNA-seq has proved advantageous for characterizing the gene expression profiles of both model and non-model species, despite only being available for a short time [17,18]. Moreover, RNA-seq has strengthened our understanding of the breadth and depth of eukaryotic transcriptomes. In fish, transcriptome profiles can be mapped and annotated by RNA-seq, and numerous biological processes such as development, host immune response, stress response, and adaptive evolution are now better understood due to this technique . RNA-seq has been applied to zebrafish, channel catfish, European sea bass, rainbow trout, and grass carp [20-24].
Grass carp (Ctenopharyngodon idellus), an important aquaculture species in China, accounts for more than 18% of total freshwater aquaculture production in this country. Production of grass carp reached 478.2 million tons in 2012, which making it the most highly consumed freshwater fish worldwide . However, much of our knowledge on grass carp genes is restricted to immunity-related or disease resistance-associated genes, and growth–related genes are not well understood [26-29]. Therefore, it is important to explore growth-related genes to inform grass carp breeding programs aimed at improving growth traits. To this end, grass carp were subjected to compensatory growth and transcriptome analysis by RNA-seq in this study. Global gene expression patterns during compensatory growth were investigated, and some significant differentially expressed genes (DEGs) were identified and annotated. These results enhance our understanding of the mechanism of compensatory growth in grass carp and will be of benefit to future grass carp breeding programs.
Animal welfare and experimental procedures were carried out in accordance with the Guide for the Care and Use of Laboratory Animals (Ministry of Science and Technology of China, 2006), and the protocol was approved by the committee of the Institute of Hydrobiology, Chinese Academy of Sciences (CAS). All surgery was performed under eugenol anesthesia, and all efforts were made to minimize suffering.
Healthy grass carp at three months old were used in the study. Grass carp weighing 2–3 g and an average length of 5 cm were obtained from the Guan Qiao Experimental Station, Institute of Hydrobiology, CAS, and acclimatized in aerated fresh water at 26-28°C for one week before processing. Fish were fed with commercial feed (The component of the feed that used in the study was provided in Additional file 1) twice a day and water was exchanged daily. If no abnormal symptoms were observed, grass carp were selected for further study. Fish were divided into control and experimental groups that each included three repeated subgroups. Each subgroup contained 65 grass carp in a separate tank. Five fish were sampled in each repeated subgroups at each time point.
Compensatory growth experiment
The compensatory growth experiment was carried out after no abnormal symptom were observed, and lasted for five weeks. At the beginning of the experiment, fish in the two groups were weighed. During the first week, fish in the control group were fed twice a day, whereas no feed was given to the experimental group. After the first week, five fish from each subgroup were weighed and muscle and liver were sampled. These samples were named as C-1-M (control, first week, muscle), C-l-L (control, first week, Liver), E-1-M (experimental, first week, muscle), and E-1-L (experimental, first week, Liver). In the second week, fish in both groups were fed twice a day, and fish were weighed and sampled at the end of the second week as described. These were named as C-2-M (control, second week, muscle), C-2-L (control, second week, Liver), E-2-M (experimental, second week, muscle), and E-2-L (experimental, second week, Liver). Tissue samples from the same subgroups were mixed equivalently for RNA isolation. The remaining fish were fed twice a day and weighed every week until the end of the experiment, and weights were subjected to statistical analysis. The Specific growth rate (SGR) was calculated as described previously . Briefly, SGR means ((In W2-W1)/(T2-T1) × 100), where W2 is the weight at the end of the growth interval and W1 is the weight at the beginning of the growth interval, while T2-T1 represents the duration (days) of the growing interval. In this study, SGRs were calculated for control and experimental group during the three time intervals: 0 ~ 1 weeks, 1 ~ 2 weeks, and 2 ~ 5 weeks. In addition, an independent repeated experiment was carried out by the similar method that described above. The repeated experiment was lasted for three weeks and samples were used for qPCR analysis.
RNA isolation, library construction and sequencing
RNA was isolated using Trizol reagent (Invitrogen, USA) according the manufacturer’s protocol, concentration was measured by the Qubit RNA assay kit (Life Technologies, USA), and integrity was assessed with the RNA nano 6000 assay kit (Agilent Technologies, USA). RNA of sufficient quality was used in library construction. Sequencing libraries were generated using the NEBNext Ultra RNA library prep kit for Illumina (New England Biolabs, USA) following the manufacturer’s protocol. Briefly, mRNA was purified from total RNA using poly-T oligo-attached magnetic beads and fragmented by NEBNext first strand synthesis reaction buffer (New England Biolabs, USA). First strand cDNA was synthesized using a random hexamer primer and M-MuLV reverse transcriptase. Second strand cDNA synthesis was subsequently performed using DNA polymerase I and RNase H. After adenylation of the 3’ end of DNA fragments, NEBNext adaptors with hairpin loop structure were ligated in preparation for hybridization. 3 μl USER enzyme (New England Biolabs, USA) was used with size-selected, adaptor-ligated cDNA at 37°C for 15 min followed by 5 min at 95°C prior to PCR using phusion high-fidelity DNA polymerase, universal PCR primers and index (X) primer. Finally, PCR products were purified using an AMPure XP system and library quality was assessed using an Agilent Bioanalyzer 2100 system. Libraries were sequenced on an Illumina Hiseq 2000 platform and 100 bp single-end reads were generated.
Raw data reads in fastq format were initially processed using in-house perl scripts. In this step, clean data (clean reads) were obtained by removing adapter, poly-N and poor quality data. The Q20, Q30, and GC content of the clean data were calculated, and all downstream analysis was performed the clean high quality data.
Clean data were mapped to the grass carp reference genome (Bioproject: PRJNA39737, unpublished data) using TopHat2 software . Two base mismatches were allowed in the mapping process, total mapped reads were calculated, and the mapped regions (exon, intron, and intergenic) were counted.
HTSeq software was used to count the number of reads mapped to each gene , and the reads per kilobase of the exon model per million mapped reads (RPKM) were calculated for each gene based on the length of the gene and the number of reads mapped to the gene .
Differential expression analysis
Differential expression analysis of two groups/conditions was performed using the DESeq package . The resulting p-values were adjusted using the Benjamini and Hochberg’s approach for controlling the false discovery rate. Genes with an adjusted p-value <0.05 (padj <0.05) found by DESeq were assigned as differentially expressed.
Gene Ontology (GO) enrichment analysis of DEGs was implemented by the GOseq R package , in which gene length bias was corrected. GO terms with corrected p-values less than 0.05 were considered significantly enriched by DEGs.
The Kyoto Encyclopedia of Genes and Genomes (KEGG) database is used for understanding high-level functional information in biological systems from molecules, cells, organisms and ecosystems, and is particularly powerful for large-scale molecular datasets generated by genome sequencing and other high-throughput experimental approaches . In this study, KOBAS software was employed to test the statistical enrichment of DEGs in KEGG pathways . KEGG terms with corrected p-values less than 0.05 were considered significant.
Validation of DEGs by qPCR
In order to confirm the reliability of data obtained by RNA-seq, twelve DEGs were random selected for validation by qPCR. The primers were listed in Additional file 2 and the cDNA sequences (completely or partially) of these genes were shown in Additional file 3. The RNA samples from an independent repeated study and were used for reverse transcription. First strand cDNAs were obtained using a random hexamer primer and the ReverTra Ace kit (Toyobo, Japan). qPCR was carried out in a Bio-rad fluorescence quantitative PCR instrument (Bio-rad, USA). Each qPCR mixture contained 0.8 μL sense and reverse primers, 1 μL template, 10 μL 2 × SYBR mix (TOYOBO, Japan), and 7.4 μL ddH2O. Three replicates were conducted for each sample and β-actin gene was used as an internal control to normalize. Only the primer with efficiency of 90% ~ 110% was used for qPCR. The program for qPCR was as follows: 95°C for 10 s, 40 cycles of 95°C for 5 s and 60°C for 20 s. Relative expression level was calculated using the 2-△△Ct method . All data are given as mean ± standard deviation of three replicates.
Changes in body weight and SGR during compensatory growth
The weight of fish in two groups was recorded at six time points and curves were drawn (Figure 1A). For the control group, a total increase of 1.16 g in body weight and a growth rate of 39.8% was obtained. For the experimental group, a total increase of 1.26 g in body weight and a growth rate of 49.5% was acquired. Moreover, the SGR in different time intervals was calculated (Figure 1B). In the first week, the weight of the experimental group decreased sharply following the induced starvation (12.2% decrease in body weight), indicated by a negative SGR (−1.84 ± 0.52). During the following week of re-feeding, the weight of the experimental group increased rapidly and resulted in a positive of SGR (3.95 ± 0.36), which is significantly higher than (P < 0.01) that in the control group (SGR = 1.64 ± 0.42). The elevated SGR that characterizes compensatory growth subsequently declined back to low level during the 2 ~ 5 weeks of realimentation, whereas the SGR of experimental group (SGR = 0.91 ± 0.14) was still significantly higher than (P < 0.01) that in the control group (SGR = 0.34 ± 0.19) (Figure 1B).
Preliminary analysis of RNA-seq data
At different time points, muscle and liver samples from control and experimental groups were used for library construction. Duplicates were performed for each class, therefore 24 libraries were constructed in total. These libraries were sequenced using the Illumina Hiseq. 2000 platform, and raw reads, clean reads, Q20, total mapped reads, and unique mapped reads for each library were recorded (Table 1). All libraries gave Q20 ≥ 95%, total mapped reads ≥89%, and unique mapped reads ≥85%. Moreover, the percentage of the total mapped reads that mapped to the genome region was calculated and ≥83% mapped to the exon for all libraries (data not shown). This confirmed the high quality of the sequencing data and suitability for further analysis. The sequencing data in this study have been deposited in the Sequence Read Archive (SRA) at the National Center for Biotechnology Information (NCBI) (accession number: SRP055685).
Analysis of gene expression level
The total reads number that mapped to the genome region of each gene was calculated using HTSeq software, and the expression level of each gene was calculated according to the method of RPKM . RPKM intervals and gene numbers in each RPKM interval of all 24 libraries are listed in Additional file 4. RPKM intervals were similar for duplicate samples. However, RPKM intervals of samples from muscle were different from liver samples. The number of genes with a low RPKM interval (1 ~ 3) was greater in liver, whereas the number of genes with a high RPKM interval (15 ~ 60) was greater in muscle. The RPKM of each gene was calculated from the average of all replicates. Moreover, the R2 value of the Pearson product–moment correlation coefficient was ≥0.92 for all replicates (data not shown), which further confirmed the repeatability of the sampling.
Identification of DEGs
From integration of the replicates, the 24 libraries were condensed into 8 sample groups named C-1-M, C-l-L, E-1-M, E-1-L, C-2-M, C-2-L, E-2-M, and E-2-L. To identify DEGs, samples were subjected to a series of paired-comparisons. Muscle and liver samples from the experimental group that experienced fasting conditions were compared with the appropriate control group (E-1-M/C-1-M and E-1-L/C-1-L). Muscle and liver samples from the experimental group following re-feeding were compared with the appropriate control group (E-2-M/C-2-M and E-2-L/C-2-L). In addition, muscle and liver samples from the experimental group following re-feeding were compared with experimental groups following fasting conditions (E-2-M/E-1-M and E-2-L/E-1-L). The number of DEGs identified from the different paired-comparisons is listed in Table 2. In fasting conditions, 4061 DEGs were detected in muscle (2124 up-regulated and 1937 down-regulated) and 1988 DEGs were identified in liver (761 up-regulated and 1227 down-regulated). Following re-feeding, 349 and 27 DEGs were discovered in muscle (281 up-regulated and 68 down-regulated) and liver (148 up-regulated and 99 down-regulated), respectively. Moreover, when fasting and re-feeding experimental groups were compared, 4903 DEGs were identified in muscle (2668 up-regulated and 2235 down-regulated) and 2444 DEGs were detected in liver (1512 up-regulated and 932 down-regulated). Detailed information of all DEGs is shown in Additional file 5. These DEGs were subjected to Venn diagram analysis (Figure 2), which identified 892 genes in both E-1-M/C-1-M and E-1-L/C-1-L (Figure 2A), 29 genes in both E-2-M/C-2-M and E-2-L/C-2-L (Figure 2B), and 1205 genes in both E-2-M/E-1-M andE-2-L/E-1-L (Figure 2C).
GO enrichment analysis
GO enrichment analysis was performed to investigate the possible roles of DEGs. For all six paired-comparisons, annotated genes were categorized into three main categories, namely biological process, molecular function, and cellular component (Figure 3, top 30 most enriched terms). The biological process category included high representation for genes involved in single-organism metabolic processes, organonitrogen compound metabolism, small molecule metabolism, oxidation-reduction, general metabolic and organic acid metabolic processes. Catalytic activity, oxidoreductase activity, cofactor binding, coenzyme binding, and other binding terms were significant enriched in the molecular function category. In the cellular component category, intracellular, myosin complex, extracellular matrix, actin cytoskeleton, and non-membrane-bound organelle terms were abundant. In the comparison of E-2-L/C-2-L, no cellular component term was enriched, and the number of terms for biological process and molecular function was also low. Detailed information of enriched terms is listed in Additional file 6.
In order to identify possible biochemical pathways in which DEGs operate, KEGG pathway analysis was carried out for all six paired-comparisons, and significantly enriched pathways are listed in Table 3. The most of enriched pathways were metabolism- or biosynthesis-related pathways such as metabolic pathways, carbon metabolism, fatty acid metabolism, glycine, serine and threonine metabolism, citrate cycle (TCA cycle), and biosynthesis of amino acids. In addition, some pathways involved in genetic information processing were also enriched, such as protein processing in endoplasmic reticulum, DNA replication, aminoacyl-tRNA biosynthesis, ribosome biogenesis in eukaryotes, and RNA transport. In all the significant enriched pathways, metabolic pathway was the top category that included the greatest number of annotated genes. In the comparison of E-2-M/C-2-M and E-2-L/C-2-L, fewer pathways were enriched, indicating less DEGs in these comparisons.
Identification of significant DEGs between experimental and control groups
The more significant DEGs may play an important role in response to changes in the environment , therefore these were identified and annotated. The top five significant DEGs (up-regulated and down-regulated) are listed in Table 4. In fasting conditions, DEGs exhibited a log2FoldChange ≥5 for both up- and down-regulated genes, indicating a significant change in expression level. Specifically, calcium-binding and coiled-coil domain-containing protein 1 (CALCOCO1), kelch-like protein 38 (KLH38), methyltransferase-like protein 21C (METTL21C), F-box only protein 32 (FBOX32), SPRY domain-containing SOCS box protein 3 (SPSB3), hepcidin-1, nociceptin receptor, insulin-like growth factor-binding protein (IGFBP) 1 and 4, and krueppel-like factor 9 (KLF-9) were all up-regulated significantly. Type-4 ice-structuring protein (AFP4), major facilitator superfamily domain-containing protein 2A (MFSD2A) and genes involved in the metabolism of glucose and fatty acids such as group 3 secretory phospholipase A2 (PLA2G3), glucokinase, fatty acid binding protein 6 (FABP6), elongation of very long chain fatty acids protein-4 (ELOVL-4), apolipoprotein A-I, and long-chain-fatty-acid-CoA ligase 1 (ACSL-1) were all down-regulated significantly. Following re-feeding, changes in expression level were less pronounced than those under fasting conditions. However, it was apparent that all up-regulated DEGs were muscle-related (myosin heavy chain (MYL), myosin-13, parvalbumin beta, parvalbumin-2, myosin regulatory light chain 2 (MYL2), myosin light chain 3 (MYL-3), and troponin C). Genes such as lovastatin nonaketide synthase, kelch domain-containing protein 1 (KLHDC-1), ankyrin repeat domain-containing protein 29 (ANKRD-29), kyphoscoliosis peptidase, ankyrin repeat and SOCS box protein 2 (ASB-2), cytosolic phospholipase A2 gamma (CPLA2-γ), growth/differentiation factor 15 (GDF-15), ankyrin repeat domain- containing protein 37 (ANKRD-37), nuclear receptor coactivator 7 (NCOA7), and ferritin (middle subunit) were all down-regulated significantly.
Identification of significant DGEs in experimental group between fasting and re-feeding conditions
To further investigate the mechanism of compensatory growth, significant DEGs in experimental group between re-feeding and fasting conditions were also identified and annotated. All listed DEGs showed a log2FoldChange ≥5, indicating a marked change in expression level between re-feeding and fasting (Table 5). LA2G3, Actin, ELOVL-4, glycine amidinotransferase (GATM), MYL, AFP4, elongation of very long chain fatty acids protein 6 (ELOVL-6), ACSL-1, and zinc finger FYVE domain-containing protein 9 (ZFYVE9)) were all up-regulated. CALCOCO1, alpha-2-HS-glycoprotein (AHSG), protein-glutamine gamma- glutamyltransferase (TGM), SPSB3, inactive dual specificity phosphatase 27 (DUSP27), calcium- independent phospholipase A2 (iPLA2), heme oxygenase (HO), solute carrier organic anion transporter family member 1C1 (SLCO1C1), solute carrier family 13 member 2 (SLC13A2, and suppressor of cytokine signaling 2 (SOCS2) were significantly down-regulated DEGs. These significant DEGs may play an important role in the response to re-feeding and fasting in teleost fish. The cDNA sequences (completely or partially) of significant DGEs could be available in Additional file 3.
Confirmation of DEGs by qPCR
To confirm the RNA-seq data, twelve DEGs were random selected for qPCR analysis. The RNA samples that form an independent repeated study and were used for reverse transcription and qPCR analysis. For each of paired-comparison, two genes were random selected. The random selected DEGs were macrophage migration inhibitory factor (MIF),peroxiredoxin 3 (PRDX3), apolipoprotein Eb (APOEb), elongation factor 1-alpha (EF-1a), apolipoprotein A-I-1 (APOA-I-1), poly (A) binding protein interacting protein 2B (PAIP2B), pleckstrin and Sec7 domain containing 2 (PSD2), fructose-bisphosphate aldolase b (ALDOb), fructose-bisphosphate aldolase a (ALDOa), complement factor D (CFD), eukaryotic translation elongation factor 1 alpha 1-like 2 (EFF1a1L2), glyceraldehyde-3-phosphate dehydrogenase (GAPDH). As shown in Figure 4, the expression patterns of all twelve DEGs that obtained by qPCR were similar to that in RNA-seq, although the relative expression level was not completely consistent. The results confirmed the reliability and accuracy of the RNA-seq data (Figure 4).
Due to changes in season, variation of temperature, unbalanced food availability, and other growth-stunting conditions, the growth and development of fish can often be restricted. However, when conditions returned to normal, fish can undergo obvious compensatory growth . Despite numerous reports on compensatory growth in fish [40-42], the mechanism remains unknown, and global gene expression patterns resulting from compensatory growth are unclear. In order to reveal the mechanism and explore growth-related genes, a compensatory growth experiment was carried out on grass carp and global gene expression patterns were determined using an RNA-seq approach. As reported previously, compensatory growth can be classified into four types: over-compensation, full compensation, partial compensation, and no compensation . In our study, the total increased weight of the experimental group was 1.26 g in five weeks, which was more than that in the control group. Moreover, the SGR of experimental group was significant higher than that in the control group during the week of re-feeding. A phenomenon of hyperphagia was observed in the experimental group during re-feeding (data not shown), which is a characteristic of compensatory growth. In addition, the RNA-seq data obtained in the study also implied that the compensation growth happened. Thus, although the initial body weight of fish in experimental group was slightly lower than that in control group, the results still demonstrated that at least partial compensation occurred.
Under fasting conditions, 4061 and 1988 DEGs were identified in muscle and liver, respectively, many of which were involved in metabolic processes, catalytic activity, binding functions, and participated in metabolic and biosynthetic pathways, according to GO and KEGG annotation results. Up- and down-regulation of metabolic and biosynthetic pathways may therefore be a strategy undertaken by grass carp in response to fasting conditions. Following re-feeding, 349 and 247 DEGs were identified in muscle and liver, respectively. Although less than that in under fasting conditions, most DEGs were up-regulated. These results suggest differences between the experimental group and control groups were reduced after re-feeding, but the differences that were apparent may contribute greatly to the fast increase in body weight observed. In addition, when samples from experimental groups in different conditions were compared, 4093 and 2444 DEGs were identified in muscle and liver, respectively, indicating extensive differences in gene expression between starved and re-feed fish. These DEGs may be particularly important for compensatory growth in grass carp.
The more significant DEGs that showed the largest changes in expression level were annotated. Under starvation, most of the significant down-regulated genes were involved in the metabolism of glucose and fatty acids (PLA2G3, Glucokinase, FABP6, ELOVL-4, ELOVL-6, Apolipoprotein A-I, and ACSL-1) [44-49]. MFSD2A and AFP4 are known to be important for growth and development [50,51], and down-regulation of these genes indicates that metabolism of glucose and fatty acids is decreased sharply during fasting conditions, which prevented growth and was consistent with the observed reduced weight under these conditions. Of the significant up-regulated genes, CALCOCO1, KLH38, FBOX32, SPSB3, and METTL21C are associated with modification, ubiquitination, or degradation of proteins [52-56], whereas the nociceptin receptor plays an important role in response to stimulation . Up-regulation of these genes suggests muscle proteins may be used as the major energy source to maintain basic metabolism in response to fasting, which may also contribute to the decreased weight. In addition, two insulin-like growth factor binding proteins, IGFBP-1 and 4, were down-regulated significantly. IGFBPs are reported to bind to insulin-like growth factors (IGFs), which prevents binding between IGFs and their cognate receptors, thereby inhibiting the activities of IGFs [9,12,58-61]. Thus, significant up-regulated of IGFBP-1 and 4 may be another reason for the decreased weight of fish under fasting conditions.
After re-feeding, the most significantly up-regulated genes were myosin-related, such as MYL, Myosin-13, Parvalbumin beta, Parvalbumin-2, MYL-2, MYL-3, and Troponin C [62-67]. This strongly indicates that biosynthesis of myosin or muscle occurred, which resulted in increased weight following re-feeding. Of the significantly down-regulated genes, Lovastatin nonaketide synthase is known to participate in the methylation of proteins . KLHDC-1, ANKRD-29, ANKRD-37, and ASB-2 are also related to the ubiquitination or degradation of proteins [69-71], and Kyphoscoliosis peptidase hydrolyzes muscle-specific proteins . GDF-15 is associated with appetite, and high expression level of GDF-15 may reduce appetite and weight in mice . Ferritin (middle subunit) is a protein involved in response to cellular emergencies . Down-regulation of these genes suggests that following re-feeding, fish exit from the emergency response, and the presence of food stimulates appetite. Subsequently, muscle proteins are no longer used as the major energy source, and rapid compensatory growth occurred.
Significant DEGs in experimental group between fasting and re-feeding conditions were also identified. PLA2G3, ELOVL-4, ELOVL-6, and ACSL-1, involved in the metabolism of glucose and fatty acids and identified above, were all significantly up-regulated. Actin, GATM, and MYL encode muscle-related proteins , and AFP4 and ZFYVE9 are important for growth and development [51,76]. Up-regulation of these genes indicates stimulation of glucose and fatty acid metabolism, and enhanced biosynthesis of muscle, which may explain the increased body weight. Of the significantly down-regulated genes, CALCOCO1 and SPSB3, associated with the ubiquitination and degradation of proteins and identified above, were down-regulated. AHSG is able to make mice insensitive to insulin, and inhibits growth . DUSP27 encodes a protein that catalyzes the hydrolysis of amino acids , while HO is important in the response to cellular emergencies . SLCO1C1 encodes a receptor that mediates uptake of thyroid hormones  and SOCS2 negatively regulates growth hormone and IGFs . Down-regulation of these genes indicates that muscle proteins are no longer used as the major energy source following re-feeding. In addition, growth hormone and IGFs may be positively regulated, which may accelerate the growth and development of grass carp.
In conclusion, grass carp were subjected to compensatory growth and global gene expression patterns were determined by RNA-seq. Numerous DEGs were identified and several significant DEGs were annotated. This study expands our understanding of the mechanism of compensatory growth, and will provide a reference for growth-related genes in grass carp.
Availability of supporting data
The sequencing data in this study could be available form the Sequence Read Archive (SRA) at the National Center for Biotechnology Information (NCBI) (accession number: SRP055685). Other supporting data are included as additional files.
Differentially expressed genes
Reads per kilobase of exon model per million mapped reads
Kyoto encyclopedia of genes and genomes
Calcium-binding and coiled-coil domain-containing protein 1
Kelch-like protein 38
Methyltransferase-like protein 21C
F-box only protein 32
SPRY domain-containing SOCS box protein 3
Insulin-like growth factor-binding protein 1
Insulin-like growth factor-binding protein 4
Krueppel-like factor 9
Group 3 secretory phospholipase A2
Fatty acid binding protein 6
Elongation of very long chain fatty acids protein-4
Major facilitator superfamily domain-containing protein 2A
Myosin heavy chain
Myosin regulatory light chain 2
Myosin light chain 3
Kelch domain-containing protein 1
Ankyrin repeat domain-containing protein 29, ASB-2, Ankyrin repeat and SOCS box protein 2
Cytosolic phospholipase A2 gamma
Growth/differentiation factor 15
Ankyrin repeat domain-containing protein 37
Nuclear receptor coactivator 7
Elongation of very long chain fatty acids protein 6
Zinc finger FYVE domain-containing protein 9
Protein-glutamine gamma-glutamyltransferase, DUSP27, Inactive dual specificity phosphatase 27
Calcium-independent phospholipase A2
Solute carrier organic anion transporter family member 1C1
Solute carrier family 13 member 2
Suppressor of cytokine signaling 2. MIF, macrophage migration inhibitory factor
Elongation factor 1-alpha
Poly (A) binding protein interacting protein 2B
Pleckstrin and Sec7 domain containing 2
Fructose-bisphosphate aldolase b
Fructose-bisphosphate aldolase a
Complement factor D
Eukaryotic translation elongation factor 1 alpha 1-like 2
Glyceraldehyde-3 -phosphate dehydrogenase
Gui JF, Zhu ZY. Molecular basis and genetic improvement of economically important traits in aquaculture animals. Chin Sci Bull. 2012;57:1751–60.
Johnston IA, Bower NI, Macqueen DJ. Growth and the regulation of myotomal muscle mass in teleost fish. J Exp Biol. 2011;214(10):1617–28.
De-santis C, Jerry DR. Candidate growth genes in finfish—where should we be looking? Aquaculture. 2007;272:22–38.
Devlin RH, Sakhrani D, Tymchuk WE, Rise ML, Goh B. Domestication and growth hormone transgenesis cause similar changes in gene expression in coho salmon (Oncorhynchus kisutch). Proc Natl Acad Sci U S A. 2009;106(9):3047–52.
Laine VN, Shikano T, Herczeg G, Vilkki J, Merilä J. Quantitative trait loci for growth and body size in the nine-spined stickleback Pungitius pungitius L. Mol Ecol. 2013;22(23):58615876.
Loukovitis D, Sarropoulou E, Tsigenopoulos CS, Batargias C, Magoulas A, Apostolidis AP, et al. Quantitative trait loci involved in sex determination and body growth in the gilthead sea bream (Sparus aurata L.) through targeted genome scan. PLoS One. 2011;6(1):e16599.
Sánchez-Molano E, Cerna A, Toro MA, Bouza C, Hermida M, Pardo BG, et al. Detection of growth-related QTL in turbot (Scophthalmus maximus). BMC Genomics. 2011;12:473.
Hutson AM, Liu Z, Kucuktas H, Umali-Maceina G, Su B, Dunham RA. Quantitative trait loci map for growth and morphometric traits using a channel catfish x blue catfish interspecific hybrid system. J Anim Sci. 2014;92(5):1850–65.
Chauvigné F, Gabillard JC, Weil C, Rescan PY. Effect of refeeding on IGFI, IGFII, IGF receptors, FGF2, FGF6, and myostatin mRNA expression in rainbow trout myotomal muscle. Gen Comp Endocrinol. 2003;132(2):209–15.
Montserrat N, Gabillard JC, Capilla E, Navarro MI, Gutiérrez J. Role of insulin, insulin-like growth factors, and muscle regulatory factors in the compensatorygrowth of the trout (Oncorhynchus mykiss). Gen Comp Endocrinol. 2007;150(3):462–72.
Nebo C, Portella MC, Carani FR, de Almeida FL, Padovani CR, Carvalho RF, et al. Short periods of fasting followed by refeeding change the expression of muscle growth-related genes in juvenile Nile tilapia (Oreochromis niloticus). Comp Biochem Physiol B Biochem Mol Biol. 2013;164(4):268–74.
Picha ME, Biga PR, Galt N, McGinty AS, Gross K, Hedgpeth VS, et al. Overcompensation of circulating and local insulin-like growth factor-1 during catch-up growth in hybrid striped bass (Morone chrysops × Morone saxatilis) following temperature and feeding manipulations. Aquaculture. 2014;428:174–83.
Won ET, Borski RJ. Endocrine regulation of compensatory growth in fish. Front Endocrinol (Lausanne). 2013;4:74.
Osborne TB, Mendel LB. Acceleration of growth after retardation. Am J Physiol. 1916;40:16–20.
Ali M, Nicieza A, Wootton RJ. Compensatory growth in fishes: a response to growth depression. Fish Fisheries. 2003;4(2):147–90.
Morozova O, Hirst M, Marra MA. Applications of new sequencing technologies for transcriptome analysis. Annu Rev Genomics Hum Genet. 2009;10:135–51.
Wang Z, Gerstein M, Snyder M. RNA-Seq: a revolutionary tool for transcriptomics. Nat Rev Genet. 2009;10:57–63.
Costa V, Angelini C, De Feis I, Ciccodicola A. Uncovering the complexity of transcriptomes with RNA-Seq. J Biomed Biotechnol. 2010;2010:853916.
Qian X, Ba Y, Zhuang Q, Zhong G. RNA-Seq technology and its application in fish transcriptomics. OMICS. 2014;18(2):98–110.
Collins JE, White S, Searle SM, Stemple DL. Incorporating RNA-seq data into the zebrafish Ensembl genebuild. Genome Res. 2012;22:2067–78.
Liu S, Zhang Y, Zhou Z, Waldbieser G, Sun F, Lu J, et al. Efficient assembly andannotation of the transcriptome of catfish by RNA-Seq analysis of a doubled haploid homozygote. BMC Genomics. 2012;13:595.
Sarropoulou E, Galindo-Villegas J, Garcıá-Alcázar A, Kasapidis P, Mulero V. Characterization of European seabass transcripts by RNA seq after oral vaccine against V. anguillarum. Mar Biotechnol. 2012;14:634–42.
Palstra AP, Beltran S, Burgerhout E, Brittijn SA, Magnoni LJ, Henkel CV, et al. Deep RNA sequencing of the skeletal muscle transcriptome in swimming fish. PLoS One. 2013;8:e53171.
Shi M, Huang R, Du F, Pei Y, Liao L, Zhu Z, et al. RNA-seq profiles from grass carp tissues after reovirus (GCRV) infection based on singular and modular enrichment analyses. Mol Immunol. 2014;61(1):44–53.
FAO. Fishery and aquaculture statistics yearbook. Rome: Food and Agriculture Oranization of the United Nations; 2013.
Chen WQ, Xu QQ, Chang MX, Zou J, Secombes CJ, Peng KM, et al. Molecular characterization and expression analysis of the IFN-gamma related gene (IFN-gammarel) in grass carp Ctenopharyngodon idella. Vet Immunol Immunopathol. 2010;134(3–4):199–207.
Du F, Su J, Huang R, Liao L, Zhu Z, Wang Y. Cloning and preliminary functional studies of the JAM-A gene in grass carp (Ctenopharyngodon idellus). Fish Shellfish Immunol. 2013;34(6):1476–684.
Cai J, Yang L, Wang B, Huang Y, Tang J, Lu Y, et al. Identification of a novel N4BP1-like gene from grass carp (Ctenopharyngodon idella) in response to GCRV infection. Fish Shellfish Immunol. 2014;36(1):223–8.
Feng X, Su J, Yang C, Yan N, Rao Y, Chen X. Molecular characterizations of grass carp (Ctenopharyngodon idella) TBK1 gene and its roles in regulating IFN-I pathway. Dev Comp Immunol. 2014;45(2):278–90.
Kim D, Pertea G, Trapnell C, Pimentel H, Kelley R, Salzberg SL. TopHat2: accurate alignment of transcriptomes in the presence of insertions, deletions andgene fusions. Genome Biol. 2013;14(4):R36.
Simon Anders, Paul Theodor Pyl, Wolfgang Huber. HTSeq-A Python framework to work with high-throughput sequencing data. BioRxiv. 2014, doi: 10.1101/002824.
Mortazavi A, Williams BA, McCue K, Schaeffer L, Wold B. Mapping and quantifying mammalian transcriptomes by RNA-Seq. Nat Methods. 2008;5(7):621–8.
Anders S, Huber W. Differential expression analysis for sequence count data. Genome Bio. 2010;11(10):R106.
Young MD, Wakefield MJ, Smyth GK, Oshlack A. Gene ontology analysis for RNA-seq: accounting for selection bias. Genome Biol. 2010;11(2):R14.
Kanehisa M, Araki M, Goto S, Hattori M, Hirakawa M, Itoh M, et al. KEGG for linking genomes to life and the environment. Nucleic Acids Res. 2008;36:D480–4.
Xie C, Mao X, Huang J, Ding Y, Wu J, Dong S, et al. KOBAS 2.0:web server for annotation and identification of enriched pathways and diseases. Nucleic Acids Res. 2011;39:W316–22.
Livak KJ, Schmittgen TD. Analysis of relative gene expression data using real-time quantitative PCR and the 2 (−Delta Delta C (T)) Method. Methods. 2001;25(4):402–8.
Townsend JP, Cavalieri D, Hartl DL. Population genetic variation in genome-wide gene expression. Mol Biol Evol. 2003;20(6):955–63.
Mommsen TP. Paradigms of growth in fish. Comp Biochem Physiol B Biochem Mol Biol. 2001;129(2–3):207–19.
Abolfathi M, Hajimoradloo A, Ghorbani R, Zamani A. Compensatory growth in juvenile roach Rutilus caspicus: effect of starvation and re-feeding on growth and digestive surface area. J Fish Biol. 2012;81(6):1880–90.
Sundström LF, Kaspersson R, Näslund J, Johnsson JI. Density-dependent compensatory growth in brown trout (Salmo trutta) in nature. PLoS One. 2013;8(5):e63287.
Zhu K, Chen L, Zhao J, Wang H, Wang W, Li Z, et al. Molecular characterization and expression patterns of myogenin in compensatory growth of Megalobrama amblycephala. Comp Biochem Physiol B Biochem Mol Biol. 2014;170:10–7.
Jobling M. Fish Bioenergetics. London: Chapman & Hall; 1994.
Breslow JL, Ross D, McPherson J, Williams H, Kurnit D, Nussbaum AL, et al. Isolation and characterization of cDNA clones for human apolipoprotein A-I. Proc Natl Acad Sci U S A. 1982;79(22):6861–5.
Dennis EA. Diversity of group types, regulation, and function of phospholipase A2. J Biol Chem. 1994;269(18):13057–60.
Börchers T, Hohoff C, Buhlmann C, Spener F. Heart-type fatty acid binding protein - involvement in growth inhibition and differentiation. Prostaglandins Leukot Essent Fatty Acids. 1997;57(1):77–84.
Malhotra KT, Malhotra K, Lubin BH, Kuypers FA. Identification and molecular characterization of acyl-CoA synthetase in human erythrocytes and erythroid precursors. Biochem J. 1999;344:135–43.
Iynedjian PB. Molecular physiology of mammalian glucokinase. Cell Mol Life Sci. 2009;66(1):27–42.
Ohno Y, Suto S, Yamanaka M, Mizutani Y, Mitsutake S, Igarashi Y, et al. ELOVL1 production of C24 acyl-CoAs is linked to C24 sphingolipid synthesis. Proc Natl Acad Sci U S A. 2010;107(43):18439–44.
Berger JH, Charron MJ, Silver DL. Major facilitator superfamily domain-containing protein 2a (MFSD2A) has roles in body growth, motor function, and lipid metabolism. PLoS One. 2012;7(11):e50629.
Xiao Q, Xia JH, Zhang XJ, Li Z, Wang Y, Zhou L, et al. Type-IV antifreeze proteins are essential for epiboly and convergence in gastrulation of zebrafish embryos. Int J Biol Sci. 2014;10(7):715–32.
Takahashi K, Inuzuka M, Ingi T. Cellular signaling mediated by calphoglin-induced activation of IPP and PGM. Biochem Biophys Res Commun. 2004;325(1):203–14.
Tintignac LA, Lagirand J, Batonnet S, Sirri V, Leibovitch MP, Leibovitch SA. Degradation of MyoD mediated by the SCF (MAFbx) ubiquitin ligase. J Biol Chem. 2005;280(4):2847–56.
Kuang Z, Lewis RS, Curtis JM, Zhan Y, Saunders BM, Babon JJ, et al. The SPRY domain-containing SOCS box protein SPSB2 targets iNOS for proteasomal degradation. J Cell Biol. 2010;190(1):129–41.
Dhanoa BS, Cogliati T, Satish AG, Bruford EA, Friedman JS. Update on the Kelch-like (KLHL) gene family. Hum Genomics. 2013;7(1):13.
Cloutier P, Lavallée-Adam M, Faubert D, Blanchette M, Coulombe B. A newly uncovered group of distantly related lysine methyltransferases preferentially interact with molecular chaperones to regulate their activity. PLoS Genet. 2013;9(1):e1003210.
Spampinato S, Di Toro R, Alessandri M, Murari G. Agonist-induced internalization and desensitization of the human nociceptin receptor expressed in CHO cells. Cell Mol Life Sci. 2002;59(12):2172–83.
Tazuke SI, Mazure NM, Sugawara J, Carland G, Faessen GH, Suen LF, et al. Hypoxia stimulates insulin-like growth factor binding protein 1 (IGFBP-1) gene expression inHepG2 cells: possible model for IGFBP-1 expression in fetal hypoxia. Proc Natl Acad Sci U S A. 1998;95(17):10188–93.
Clemmons DR. Use of mutagenesis to probe IGF-binding protein structure/function relationships. Endocr Rev. 2001;22(6):800–17.
Zhou R, Diehl D, Hoeflich A, Lahm H, Wolf E. IGF-binding protein-4: biochemical characteristics and functional consequences. J Endocrinol. 2003;178(2):177–93.
Li M, Li Y, Lu L, Wang X, Gong Q, Duan C. Structural, gene expression, and functional analysis of the fugu (Takifugu rubripes) insulin-like growth factor binding protein-4 gene. Am J Physiol Regul Integr Comp Physiol. 2009;296(3):R558–66.
Strehler EE, Periasamy M, Strehler-Page MA, Nadal-Ginard B. Myosin light-chain 1 and 3 gene has two structurally distinct and differentially regulated promoters evolving at different rates. Mol Cell Biol. 1985;5(11):3168–82.
Iwami Y, Ojima T, Inoue A, Nishita K. Primary structure of myosin heavy chain from fast skeletal muscle of Chum salmon Oncorhynchus keta. Comp Biochem Physiol B Biochem Mol Biol. 2002;133(2):257–67.
Sachdev S, Raychowdhury MK, Sarkar S. Human fast skeletal myosin light chain 2 cDNA: isolation, tissue specific expression of the single copy gene, comparative sequence analysis of isoforms and evolutionary relationships. DNA Seq. 2003;14(5):339–50.
Gillis TE, Marshall CR, Tibbits GF. Functional and evolutionary relationships of troponin C. Physiol Genomics. 2007;32(1):16–27.
Bloemink MJ, Deacon JC, Resnicow DI, Leinwand LA, Geeves MA. The superfast human extraocular myosin is kinetically distinct from the fast skeletal IIa, IIb, and IId isoforms. J Biol Chem. 2013;288(38):27469–79.
Van Do T, Hordvik I, Endresen C, Elsayed S. The major allergen (parvalbumin) of codfish is encoded by at least two isotypic genes: cDNAcloning, expression and antibody binding of the recombinant allergens. Mol Immunol. 2003;39(10):595–602.
Heever JP, Hutchinson CR, Vederas JC. Lovastatin nonaketide synthase catalyses an intramolecular Diels-Alder reaction of a substrate analogue. J Am Chem Soc. 2000;122(46):11519–20.
Heuzé ML, Guibal FC, Banks CA, Conaway JW, Conaway RC, Cayre YE, et al. ASB2 is an Elongin BC-interacting protein that can assemble with Cullin 5 and Rbx1 toreconstitute an E3 ubiquitin ligase complex. J Biol Chem. 2005;280(7):5468–74.
Huang J, Chen F, Del Casino C, Autino A, Shen M, Yuan S, et al. An ankyrin repeat-containing protein, characterized as a ubiquitin ligase, is closely associated with membrane-enclosed organelles and required for pollen germination and pollen tube growth in lily. Plant Physiol. 2006;140(4):1374–83.
Sekine Y, Hatanaka R, Watanabe T, Sono N, Iemura S, Natsume T, et al. The Kelch repeat protein KLHDC10 regulates oxidative stress-induced ASK1 activation by suppressing PP5. Mol Cell. 2012;48(5):692–704.
Blanco G, Coulton GR, Biggin A, Grainge C, Moss J, Barrett M, et al. The kyphoscoliosis (ky) mouse is deficient in hypertrophic responses and is caused by a mutation in a novel muscle- specific protein. Hum Mol Genet. 2001;10(1):9–16.
Macia L, Tsai VW, Nguyen AD, Johnen H, Kuffner T, Shi YC, et al. Macrophage inhibitory cytokine 1 (MIC-1/GDF15) decreases food intake, body weight and improves glucose tolerance in mice on normal & obesogenic diets. PLoS One. 2012;7(4):e34868.
Lee JH, Pooley NJ, Mohd-Adnan A, Martin SA. Cloning and Characterisation of Multiple Ferritin Isoforms in the Atlantic Salmon (Salmo salar). PLoS One. 2014;9(7):e103729.
Cullen ME, Yuen AH, Felkin LE, Smolenski RT, Hall JL, Grindle S, et al. Myocardial expression of the arginine: glycine amidinotransferase gene is elevated in heart failure and normalized after recovery: potential implications for local creatine synthesis. Circulation. 2006;114(1 Suppl):I16–20.
Tsukazaki T, Chiang TA, Davison AF, Attisano L, Wrana JL. SARA, a FYVE domain protein that recruits Smad2 to the TGF beta receptor. Cell. 1998;95(6):779–91.
Mathews ST, Singh GP, Ranalletta M, Cintron VJ, Qiang X, Goustin AS, et al. Improved insulin sensitivity and resistance to weight gain in mice null for the Ahsg gene. Diabetes. 2002;51(8):2450–8.
Lountos GT, Tropea JE, Waugh DS. Structure of human dual-specificity phosphatase 27 at 2.38 A resolution. Acta Crystallogr D Biol Crystallogr. 2011;67(5):471–9.
Piantadosi CA, Withers CM, Bartz RR, MacGarvey NC, Fu P, Sweeney TE, et al. Heme oxygenase-1 couples activation of mitochondrial biogenesis to anti-inflammatory cytokine expression. J Biol Chem. 2011;286(18):16374–85.
van der Deure WM, Hansen PS, Peeters RP, Kyvik KO, Friesema EC, Hegedüs L, et al. Thyroid hormone transport and metabolism by organic anion transporter 1C1 and consequences of genetic variation. Endocrinology. 2008;149(10):5307–14.
Dey BR, Spence SL, Nissley P, Furlanetto RW. Interaction of human suppressor of cytokine signaling (SOCS)-2 with the insulin-like growth factor-I receptor. J Biol Chem. 1998;273(37):24095–101.
This work was funded by the National Natural Science Foundation of China (No. 31130055), National High Technology Research and Development Program (No. 2011AA100403), and Direction Program of Chinese Academy of Sciences (No. KSCX2-EW-N-004-3).
The authors declare that they have no competing interests.
LBH contributed to bioinformatics analysis and wrote the manuscript. YYP participated in primers design and qPCR analysis. YJ, YML, and LJL performed sample collection and RNA isolation. ZYZ and YPW contributed to the design of the research and manuscript editing. All authors read and approved the final manuscript.
The components of the feed that used in the study.
Sequences and efficiencies of primers that used in qPCR analysis.
The cDNA sequences (completely or partially) of DEGs that mentioned in the study.
Summary of RPKM level in different samples.
Detail information of DGEs in different comparisons.
Top 30 GO enriched terms in different comparisons.
About this article
Cite this article
He, L., Pei, Y., Jiang, Y. et al. Global gene expression patterns of grass carp following compensatory growth. BMC Genomics 16, 184 (2015). https://doi.org/10.1186/s12864-015-1427-2
- Grass carp
- Compensatory growth
- Transcriptome analysis
- RNA sequencing
- Differentially expressed genes