- Research article
- Open Access
Transcriptome analysis reveals the molecular mechanisms underlying growth superiority in a novel grouper hybrid (Epinephelus fuscogutatus♀ × E. lanceolatus♂)
BMC Genetics volume 17, Article number: 24 (2016)
Groupers (Epinephelus spp.) have been widely cultivated in China and South-East Asian countries. As a novel hybrid offspring crossed between E. fuscogutatus♀ and E. lanceolatus♂, Hulong grouper exhibits significant growth superiority over its female parent, which made it a promising farmed species in grouper aquaculture industry in China. Hulong grouper present a good combination of beneficial traits from both parent species, but the molecular mechanisms of its heterosis still remain poorly understood.
Based on RNA sequencing and gene expression profiling, we conducted comparative transcriptome analyses between Hulong grouper and its parents E. fuscoguttatus & E. lanceolatus. Six hundred sixty-two and 5239 differentially expressed genes (DEGs) were identified in the brains and livers, respectively. GO enrichment analysis of these DEGs revealed that metabolic process and catalytic activity were the most enriched GO terms. Further analysis showed the expressions of GnRH1and GnRH3 in the brain, and GH/IGF axis related genes such as IGF-1, IGF-2b, IGFBP-1, IGFBP-2, IGFBP-4 and IGFBP-5a in the liver of the hybrid F1 were significantly up-regulated, which is in accordance with the growth superiority of hybrid grouper. Meanwhile, expressions of genes related to the protein and glycogen synthesis pathway, such as PI3KC, PI3KR, Raptor, EIF4E3, and PP1 were up-regulated, while PYG expression was down-regulated. These changes might contribute to increased protein and glycogen synthesis in the hybrid grouper.
We identified a number of differentially expressed genes such as GnRH1 and GnRH3, and genes involved in GH/IGF axis and its downstream signaling pathways for protein and glycogen synthesis in Hulong Grouper. These findings provided molecular basis underlying growth superiority of hybrid grouper, and comprehensive insights into better understanding the molecular mechanisms and regulative pathways regulating heterosis in fish.
As economically important fish species in marine aquaculture, groupers (Epinephelus spp.) are known for their delicious taste, tender flesh and rich nutrition . In the past decades, grouper industry has developed rapidly and many grouper species have been widely cultivated in China and South-East Asian countries [2, 3]. However, the sustainable development of grouper industry has been threatened by the degradation of germplasm resources and availability of grouper fries in hatcheries . To resolve these problems, hybridization technology has been introduced into grouper artificial breeding.
Hybridization is defined as a successful mating strategy for two species with one or more heritable traits, and is often used in artificial breeding to obtain potentially desirable traits in their offsprings [5, 6]. In grouper species, hybridization was firstly achieved between white-spotted green grouper (E. amblycephalus) and red grouper (E. akaara) . From then on, more efforts have been made in developing hybrid groupers and about eight hybrids with various favorable traits, such as faster growth and development, stronger immunity and higher survival rate were obtained [8–12].
Brown-marbled grouper (E. fuscoguttatus, Efu) and Giant grouper (E. lanceolatus, Ela) are both important breeding fishes with different growth rates and disease resistance. The former is a slow growing but long-lived species with high disease resistance [13, 14], while Giant grouper is popular breeding species for its rapid growth, reaching up to 3 kg in the first year . Previously, by fertilizing the eggs of Brown-marbled groupers with the sperms of Giant grouper, we obtained a novel hybrid grouper, Hulong grouper (Hyb). It was shown that Hulong grouper combines physiological features of both parental species , and its growth rates is 35.9 % faster than that of the maternal Brown-spotted grouper . According to our preliminary data of growth (Zhang Y. et al., unpublished data), at the age of 18 months, the hybrid Grouper (Hyb) reaches to 1174.7 ± 264.2 g (n = 10), while the paternal Grouper (Ela) is only 838.1 ± 168.2 g (n = 15). Although the sample size was not large enough, the overall result was statistically significant and demonstrated that the hybrid Grouper (Hyb) grows faster than the paternal Grouper (Ela) at the age of 18 months. In recent years, many researches have examined their biological and physiological features of the hybrid groupers, however, very little is known about its genetic mechanisms of heterosis.
High-throughput RNA sequencing technology (RNA-Seq) has been used as an effective tool for transcriptome analysis, aiming to discover, profile and quantify RNA transcripts . RNA-Seq combines the advantages of microarray and EST-sequencing, including single-base resolution, high throughput, low background noise, and high sensitivity, which make it feasible for mapping of transcribed regions, quantification of gene expression levels, and distinction of different isoforms and allelic expression [18–22]. In recently decades, RNA-Seq has been widely used to identify differential genes expressions related to heterosis in crop plants and cultured fishes [23–26].
In teleosts, growth hormone/insulin-like growth factor (GH/IGF) axis, including the effector hormones GH, IGFs (IGF-1 and IGF-2) and a regulatory feedback loop between them, plays an important role in regulating somatic growth [27, 28]. GH is a pituitary hormone regulating various physiological processes such as somatic growth, behavior, immune function, lipid and protein metabolism, osmoregulation, and feeding behavior in fishes [27, 29, 30]. It was widely acknowledged that GH stimulates hepatic and circulating IGF-1 levels in teleosts [31, 32], whereas IGF system, consisting of IGFs (IGF-1 and IGF-2), cell surface receptors (IGF-1R and IGF-2R) and IGF-binding proteins (IGFBP-1 to −6), plays a crucial role in cell growth, proliferation and differentiation . Among them, IGF-1 is positively correlated with growth rate and IGF-1 levels in serum is served as a growth index in fishes [34, 35], as demonstrated in Nile tilapia (Oreochromis niloticus) , Mozambique tilapia (Oreochromis mossambicus)  and Mud carp (Cirrhinus molitorella) .
In our study, comparative transcriptomic analysis for hybrid grouper and its parents Efu and Ela was conducted to explore the molecular basis underlying the growth superiority of hybrid grouper. The brain is an important organ that centralized control over the other organs of the body by driving secretion of various hormones to regulate the growth of soma. The liver is also critical for its major roles in metabolism with numerous functions in the fish, including regulation of glycogen storage, plasma protein synthesis, and hormone production. Therefore, the brain and liver were selected as the genetic mechanism targets. Our results provided comprehensive data with respect to differential gene expressions in the GH-IGF pathway and its downstream signaling pathways involved in the heterosis of the hybrid grouper. This will contribute towards improving our understanding of the molecular mechanisms and regulative pathways regulating heterosis in fish.
Transcriptome sequencing and alignment to the reference grouper genome
To provide transcriptome profiling of the hybrid F1 and its parents, brain and liver cDNA libraries from the three groupers (Hyb and its parents Efu & Ela) were separately prepared and subjected to RNA-Seq analysis in BGI-Shenzhen. In total, 323 million high-quality clean reads were obtained after data filtering (Table 1). Subsequently, the clean reads were aligned to the E. coioides genome (Zhang Y. et al., unpublished data) using SOAP aligner 2.0 . About 46–61 % of reads were mapped to the reference genome (≤5 base mismatches), in which 28–41 % of reads were mapped to the gene regions (Table 1). In addition, 28–40 % of those uniquely mapped reads were further used for gene quantification analysis (Table 1). The species of reference genome is orange-spotted grouper (E. coioides), which is different from the hybrid grouper and its two parents. The above-mentioned low mapping rates could mainly result from the species differences.
DEGs among the hybrid F1 and its parents Efu & Ela
RPKMs (Reads Per Kilobase Transcriptome per Million mapped reads) were used to quantify gene expression levels. The RPKM values of each gene in the brain and liver tissues of the hybrid grouper were compared with those in its parents Efu & Ela, respectively (Additional file 1). Through filtering by the criteria that false discovery rate (FDR) ≤0.05 and absolute log2 (ratio) ≥1, we identified DEGs between the hybrid grouper and its parents (Fig. 1a and b) (Additional files 2, 3, 4, 5). In the brain, there were 285 and 467 up-regulated genes in the Hyb compared to Ela and Efu, respectively, and the corresponding numbers of down-regulated genes were 467 and 257, and a total of 631 up-regulated genes and 992 down-regulated genes were identified in Efu relative to Ela (Fig. 1a). In the liver, the Hyb had 5264 and 2699 down-regulated genes compared to Ela and Efu, respectively, and Ela had 3497 down-regulated genes relative to Efu (Fig. 1b). These results suggest more divergent distribution of gene expression exists in the livers than in the brains. We also determined whether the genes differentially expressed between the parents overlap the genes that differ between each parent and the hybrid. In the brain, the number of overlap DGEs between ‘Ela vs Efu’ and ‘Ela vs Hyb’ was 393, and the number of overlap DGEs between ‘Ela vs Efu’ and ‘Ela vs Hyb’ was 470 (Fig. 1c). In the liver, there were 4629 DGEs between the parents, including 2789 and 1632 overlap DGEs when the Hyb was compared to Ela and Efu, respectively (Fig. 1d). The number of DGEs among three species in the brain and liver was 78 and 633, respectively (Fig. 1c and d). The hierarchical clustering maps of these DGEs were presented in Additional file 6: Figure S1 and Additional file 7: Figure S2. In the brain, 35 and 52 DE genes in the hybrid show overdominance and underdominance, respectively. In the liver, there were 595 and 1508 DE genes expressed more highly or lowly in the hybrid than in either parent, respectively.
GO enrichment of DEGs
Gene ontology (GO) annotation of genes was achieved by Blast2GO , and. 53,665 GO terms were found and then assigned to 9503 annotated genes in the reference genome. Among them, 13,244 terms is related to cellular components, 14,191 terms for molecular function and 26,230 terms for biological process (Additional file 8: Figure S3).
The identified DEGs were subsequently used for enrichment analysis by GO:Termfinder software using the hypergeometric test [41, 42], and P-values were corrected using Bonferroni method . Being selected significantly enriched GO terms with Q-value < 0.05, the liver DEGs between Hyb and Ela, and between Hyb and Efu were enriched into 40 and 22 GO terms respectively, which provides an overview of ontology content (Additional files 9 and 10). In the molecular function and biological process categories, metabolic process (GO:0008152) and catalytic activity (GO:0003824) were the most enriched GO terms (Fig. 2a and b), suggesting that extensive metabolic and catalytic activities in the liver might be associated with growth superiority of the hybrid grouper.
GO enrichment analysis was also performed for the brain DEGs between Hyb and its parents, and only 5 GO terms were enriched in the brain (Fig. 2c and d, Additional files 11 and 12), which is consistent with the fact that less number of DEGs were identified in brains than that in livers of these three grouper species.
Comprehensive analysis of DEGs
We further analyzed the DEGs in GH/IGF axis and its downstream signaling pathways, which is critical in the growth and development of fish. In brains, we found GH and the members of Gonadotropin-releasing hormone (GnRH), GnRH1 and GnRH3 were up-regulated in Hyb over the Efu (Fig. 3a). And genes in the GH/IGF axis, including IGF-1, IGF-2b, IGFBP-1, IGFBP-2a, IGFBP-4, IGFBP-5a and IGFBP-5b were differentially expressed in livers of the Hyb over its parents (Fig. 3b). In addition, the DEGs were identified in protein and glycogen synthesis signaling pathways, including PI3K regulatory subunit (PI3KR), PI3K catalytic subunit (PI3KC), regulatory associated protein of mTOR (Raptor), eukaryotic translation initiation factor 4E (EIF4E3), protein phosphatase 1 (PP1), and glycogen phosphorylase (PYG) (Fig. 3b).
Validation DEGs by quantitative real-time PCR (qRT-PCR)
A subset of night important DEGs was randomly selected for qRT-PCR validation in livers and brains of the three grouper species, and the qRT-PCR primers were designed based on the mapped sequences (see the primer and sequence list in Additional file 13). The results were further compared with those generated from the transcriptome sequencing. Our results showed that the data from the two different methods are consistent (Fig. 4).
In summary, our data demonstrated that differential expression of GnRH, and genes related to GH/IGF axis and its downstream signaling pathways may be related to growth superiority in the hybrid grouper by increasing protein and glycogen synthesis. Our finding provides insights into molecular basis underlying the heterosis of hybrids in fish.
Hybridization is one of the most effective ways to improve physiological properties of species . Over the years, many research efforts have been made in exploring the molecular basis of heterosis. The recent research progress includes the study on DEGs and its association with heterosis [25, 45]. For instance, 3488 DEGs were identified between the parental lines of rice and their F1 hybrids and further GO enrichment analysis of these DEGs revealed that metabolism categories were the most enriched terms, indicating metabolism pathway related DEGs might contribute to the heterosis of hybrid rice .
Here we conducted comparative transcriptomic study on hybrid Hulong and its parents Efu and Ela. 323 million high-quality reads were generated from the brains and livers of Hyb and its parental species. On average, 54.12 % of reads were mapped to the E. coioides genome (Zhang Y. et al., unpublished data), and 34.60 % of the reads were mapped to gene regions. Among the annotated transcripts, 15,107 DEGs in the liver and 2948 DEGs in the brain were identified, suggesting that the DEGs in the liver of Hulong grouper might play a more important role in its heterosis than those in its brain.
DEGs associated with the GH/IGF axis
The growth of vertebrates including fish is primarily regulated by the GH/IGF axis . In present study, we found several growth-associated genes differentially expressed in the Hyb over its parents, such as GnRH1, GnRH3, and GH in the brain, and IGF-1, IGF-2b, IGFBP-1, IGFBP-2a, IGFBP-4, IGFBP-5a and IGFBP-5b in the liver (Fig. 5). Among them, the expression levels of IGF-1 and IGF-2b in the hybrid F1 liver were significantly higher than those in their parents, indicating IGF-1 and IGF-2b might play a critical role in the growth of hybrid grouper. Although the expression of GH in the brain of hybrid grouper is only slightly higher than that in its paternal Ela, the possible explanation is that GH might be inhibited by the negative feedback from increased IGF-I level in the hybrid grouper .
Many neuroendocrinological factors such as GnRH1, GnRH3, pituitary adenylate cyclase activating polypeptide 1 (PACAP1), PACAP1 type 1a/b receptor (PACAP1ra/b1), thyrotropin releasing hormone (Trh), and dopamine receptor D1a (Drd1a) promote GH secretion in the brain [47–50]. Consistently, we found the gene expression level of GnRH1 and GnRH3 is higher in the hybrid grouper, compared with that in the maternal Efu , and the gene expression of PACAP1, PACAPra/b1, Trh, and Dra1a in the Hyb is slightly higher in Hyb than its maternal Efu.
IGFBPs were reported to bind IGFs with higher affinities than IGFR, leading to prevention of IGFs degradation, thereby prolonging the half-lives of IGFs in serum . In the present study, we found the significant higher expression levels of IGFBP-2a and IGFBP-4 in Hyb than those in its parents, indicating that IGFBP-4 and IGFBP-2 might be the dominant circulating binding proteins in groupers, as opposed to IGFBP3 in mammals . In addition, our results indicated that the accumulation of IGFBP-1, IGFBP-5a and IGFBP-5b mRNAs in Hyb were intermediated between the two parents, while members of IGFBP family exhibited distinct expression patterns among these three fish species, demonstrating the independent physiological roles of its individual IGFBP family members . Taken together, our results suggest that the differential gene expression related to GH/IGF axis might play an important role in the growth superiority of the hybrid grouper.
DEGs involved in protein and glycogen synthesis
IGF-1 activates a series of phosphorylation events including PI3K/AKT pathway by binding its receptor IGFR1 , thus lead to many anabolic effects including protein synthesis  and glycogen synthesis . Upon IGF-1 stimulation, PI3K/AKT activates mTOR , then Raptor is combined with mTOR to form mTORC1, which alters the activity of EIF4E-BP1 by phosphorylation . The phosphorylated EIF4E-BP1 favors the dissociation of EIF4E from the inhibitory complex of EIF4E and EIF4E-BP1 . Ultimately, the increased availability of EIF4E promotes the elevation of protein synthesis . In our study, the higher expression level of PI3KC, PI3KR, Raptor and EIF4E3 was observed in the hybrid F1 (Fig. 5), which suggests the strengthened capacity of protein synthesis in the hybrid grouper.
Another function of the PI3K/AKT pathway is to facilitate glycogen synthesis through generating activated PI3K and stimulating additional kinases, especially AKT (protein kinase B), and further activates protein PP1 (KEGG 04910) . Subsequently, the activated PP1 inhibits PYG activity and then accelerates the glycogenesis process . In line with this, our study showed that the expression of many genes involved in the PI3K/AKT pathway, such as PI3KR, PI3KC, and PP1 were up-regulated in the liver of Hyb compared with that in the maternal Efu, while the expression levels of PYG in Hyb were lower than those in Efu, which potentially leads to the maintenance of the glycogenesis in Hyb.
Using RNA-seq, we identified differentially expressed genes in Hulong grouper. GO enrichment results indicated that these genes were related to a variety of molecular functions, such as metabolic process and catalytic activity. Further analysis revealed that differential gene expression in GH/IGF axis and downstream signaling pathways might contribute to growth superiority of the hybrid grouper through enhancing protein synthesis and glycogen synthesis.
Fish and sample preparation
E. fuscoguttatus (♀), E. lanceolatus (♂) and their hybrid F1 (Hulong) were cultivated under same breeding conditions in Daya Bay Seawater Fish Farm in Huizhou, Guangdong Province, China. At the age of 18 months, three grouper species were obtained and accommodated under same natural conditions of food, water, light, and density for 3 days. Subsequently, three individuals from each group were randomly selected for further experiments. Fresh brain and liver tissues were respectively collected and immediately stored in liquid nitrogen. All experiments were performed in accordance with the guidelines of the Animal Ethics Committee and were approved by the Institutional Review Board on Bioethics and Biosafety of BGI.
RNA extraction, libraries construction and high-throughput sequencing
Total RNA from the brain and liver tissues was extracted with Trizol reagent (Invitrogen, Carlsbad, CA, USA) and purified using RNeasy Animal Mini Kit (Qiagen, Valencia, CA) according to the manufacturer’s instructions. RNA Integrity Number (RIN) values were measured for all samples. Equal amount of total RNA from three fish individuals in each group were pooled and used to construct cDNA libraries for RNA sequencing.
The cDNA libraries were prepared following the manufacturer’s instructions (Illumina, San Diego, CA). In brief, the poly-A containing mRNA molecules was purified using poly-T oligo-attached magnetic beads. Subsequently, the mRNAs were fragmented into small pieces. The cleaved RNA fragments were copied into first strand cDNA using reverse transcriptase and random primers, and the second strand cDNA was synthesized using DNA polymerase I and RNase H. These cDNA fragments were ligated with the adapters, and these products were then purified and enriched with PCR to create the final cDNA libraries. Finally, the cDNA libraries were sequenced through Illumina HiSeq2000 system at BGI-tech (Shenzhen, China). Raw reads were submitted to SRA database in NCBI under accession number SRP056564.
Processing of raw reads and quantification of differential gene expression levels
Raw reads were filtered by removing adaptors, reads with more than 5 % unknown nucleotides, and sequences shorter than 20 nt and low quality with Q <20. The resultant clean reads were mapped to the orange-spotted grouper (E. coioides) genome (Zhang Y. et al., unpublished data; using SOAP aligner 2.0 ), which has been de novo assembled at BGI-Shenzhen, China.
The quantification and differential expressed genes analyses were carried out by the Cufflink program 2.1.1. RPKMs (Reads Per Kilobase Transcriptome per Million mapped reads) were used to quantify the mapped whole gene expression levels. The differential expressed genes were determined through filtering by the criteria that false discovery rate (FDR) ≤0.05 and absolute log2 (ratio) ≥1.
Validation of deep sequencing
Nine genes with significant differential expression among Hyb and its parents were selected randomly to validate the deep sequencing results. After purification of total RNA using SV Total RNA Isolation System (Promega, USA), reverse transcription reactions were performed to synthesize the cDNAs. Subsequent real-time PCRs were carried out using gene-specific primers (Additional file 10). β-actin was used as an internal control and the 2-ΔΔCt method  was used to calculate relative expression amounts. All samples were examined in triplicate.
Availability of supporting data
The data sets supporting the results of this article are included within the article and its additional files. The transcriptome reads produced in this study have been deposited in the National Center for Biotechnology Information (NCBI) SRA database with accession number of SRX974310 & SRX 974314, SRX 969045 & SRX974311, and SRX974309 & SRX974313 for the brain & the liver of the Hyb, Ela and Efu respectively. Access to the data is available upon publication at http://www.ncbi.nlm.nih.gov/sra/.
α-melanocyte stimulating hormone
protein kinase B
differentially expressed genes
dopamine receptor D1a
eukaryotic translation initiation factor 4E
eukaryotic translation initiation factor 4E binding protein 1
expressed sequence tags
false discovery rate
insulin-like growth factor
IGF binding protein
pituitary adenylate cyclase activating polypeptide 1
phosphoinositide-3-kinase catalytic subunit
phosphoinositide-3-kinase regulatory subunit
protein phosphatase 1
regulatory associated protein of Mtor
reads per kilobase transcriptome per million mapped reads
thyrotropin releasing hormone
Yashiro R. Statue of grouper breeding and culture in Thailand; 1996.http://library.enaca.org/Grouper/Research/Breeding/2000/0803.htm. Accessed 3 August 1998.
Guo CY, Huang YH, Wei SN, Ouyang ZL, Yan Y, Huang XH, et al. Establishment of a new cell line from the heart of giant grouper, Epinephelus lanceolatus (Bloch), and its application in toxicology and virus susceptibility. J Fish Dis. 2013;38:175–86.
Pierre S, Gaillard S, Prévot‐D’Alvise N, Aubert J, Rostaing‐Capaillon O, Leung‐Tack D, et al. Grouper aquaculture: Asian success and Mediterranean trials. Aquat Conserv Mar Freshwat Ecosyst. 2008;18(3):297–308.
Kiriyakit A, Gallardo WG, Bart AN. Successful hybridization of groupers (< i > Epinephelus coioides</i > x < i > Epinephelus lanceolatus</i>) using cryopreserved sperm. Aquaculture. 2011;320(1):106–12.
Nieto Feliner G, Fuertes Aguilar J. Hybrids and hybrid zones. Trends Ecol Evol. 1998;13(7):282.
Urbanelli S, Porretta D, Mastrantonio V, Bellini R, Pieraccini G, Romoli R, et al. Hybridization, Natural Selection and Evolution of Reproductive Isolation: A 25-Years Survey of an Artificial Sympatric Area between Two Mosquito Sibling Species of the Aedes Mariae Complex. Evolution. 2014;68:3030–8.
Tseng W, Poon C. Hybridization of < i > Epinephelus</i > species. Aquaculture. 1983;34(1):177–82.
James C, Al-Thobaiti S, Rasem B, Carlos M. Potential of grouper hybrid (Epinephelus fuscoguttatus x E. polyphekadion) for Aquaculture. Naga. 1999;22(1):19–23.
Glamuzina B, Tutman P, Skaramuca B. Hybridization of Mediterranean groupers: Epinephelus marginatus♀ × E. aeneus♂ and early development. Aquac Res. 1999;30(8):625–8.
Glamuzina B, Glavić N, Skaramuca B, Kozul V, Tutman P. Early development of the hybrid < i > Epinephelus costae</i > ♀ × < i > E. marginatus</i > ♂. Aquaculture. 2001;198(1):55–61.
YZ Liufu HZ, Liu XC, Lin HR, Huang GG. Preliminary Study on the Hybrid Red-spotted Grouper(Epinephelus akaara)♂ × Orange-spotted Grouper(Epinephelus coioides)♀. Acta Sci Nat Univ Sunyatseni. 2007;46(3):72–5.
Yan Wang YZ, Haifa Z, Huang W, Haoran L. Comparision of morphology between two kinds of hybrid groupers and their parents. J Fish China. 2014;38(6):778–85.
Grandcourt E. Demographic Characteristics of Selected Epinepheline Groupers (family Serranidae, Subfamily Epinephelinae) from Aldabra Atoll, Seychelles. Washington: National Museum of Natural History, Smithsonian Institution; 2005.
Apines-Amar MJS, Amar EC, Faisan Jr JP, Pakingking Jr RV, Satoh S. Dietary onion and ginger enhance growth, hemato-immunological responses, and disease resistance in brown-marbled grouper, Epinephelus fuscoguttatus. AACL Bioflux. 2012;5(4):231–9.
Sadovy Y, Donaldson T, Graham T, McGilvray F, Muldoon G, Phillips M, et al. While stocks last: The live reef food fish trade. Manila: Asian Development Bank; 2003.
Center DA. Qinglong grouper and Hulong grouper. J Ocean Fish. 2010;9:37.
Zhou H: Genetic analysis for two grouper hybrids and their parents. Master degree thesis of Sun Yat-sen University. 2011;1:1–2.
Wang Z, Gerstein M, Snyder M. RNA-Seq: a revolutionary tool for transcriptomics. Nat Rev Genet. 2009;10(1):57–63.
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.
Vera JC, Wheat CW, Fescemyer HW, Frilander MJ, Crawford DL, Hanski I, et al. Rapid transcriptome characterization for a nonmodel organism using 454 pyrosequencing. Mol Ecol. 2008;17(7):1636–47.
Wilhelm BT, Landry J-R. RNA-Seq—quantitative measurement of expression through massively parallel RNA-sequencing. Methods. 2009;48(3):249–57.
Ozsolak F, Milos PM. RNA sequencing: advances, challenges and opportunities. Nat Rev Genet. 2011;12(2):87–98.
Song G, Guo Z, Liu Z, Qu X, Jiang D, Wang W, et al. The phenotypic predisposition of the parent in F1 hybrid is correlated with transcriptome preference of the positive general combining ability parent. BMC Genomics. 2014;15(1):297.
Bi Y-M, Meyer A, Downs GS, Shi X, El-Kereamy A, Lukens L, et al. High throughput RNA sequencing of a hybrid maize and its parents shows different mechanisms responsive to nitrogen limitation. BMC Genomics. 2014;15(1):77.
Zhai R, Feng Y, Wang H, Zhan X, Shen X, Wu W, et al. Transcriptome analysis of rice root heterosis by RNA-Seq. BMC Genomics. 2013;14:19.
Gao Y, Zhang H, Gao Q, Wang L, Zhang F, Siva VS, et al. Transcriptome Analysis of Artificial Hybrid Pufferfish Jiyan-1 and Its Parental Species: Implications for Pufferfish Heterosis. PLoS One. 2013;8(3):e58453.
Duan C. The insulin-like growth factor system and its biological actions in fish. Am Zool. 1997;37(6):491–503.
Duan C. Nutritional and developmental regulation of insulin-like growth factors in fish. J Nutr. 1998;128(2 Suppl):306S–14S.
Albalat A, Gomez-Requeni P, Rojas P, Medale F, Kaushik S, Vianen GJ, et al. Nutritional and hormonal control of lipolysis in isolated gilthead seabream (Sparus aurata) adipocytes. Am J Physiol Regul Integr Comp Physiol. 2005;289(1):R259–65.
Kawauchi H, Sower SA. The dawn and evolution of hormones in the adenohypophysis. Gen Comp Endocrinol. 2006;148(1):3–14.
Pierce AL, Dickey JT, Larsen DA, Fukada H, Swanson P, Dickhoff WW. A quantitative real-time RT-PCR assay for salmon IGF-I mRNA, and its application in the study of GH regulation of IGF-I gene expression in primary culture of salmon hepatocytes. Gen Comp Endocrinol. 2004;135(3):401–11.
Vong QP, Chan KM, Cheng CH. Quantification of common carp (Cyprinus carpio) IGF-I and IGF-II mRNA by real-time PCR: differential regulation of expression by GH. J Endocrinol. 2003;178(3):513–21.
Ferry Jr RJ, Cerri RW, Cohen P. Insulin-like growth factor binding proteins: new proteins, new functions. Horm Res. 1999;51(2):53–67.
Reinecke M. Influences of the environment on the endocrine and paracrine fish growth hormone-insulin-like growth factor-I system. J Fish Biol. 2010;76(6):1233–54.
Adams MD, Kelley JM, Gocayne JD, Dubnick M, Polymeropoulos MH, Xiao H, et al. Complementary DNA sequencing: expressed sequence tags and human genome project. Science. 1991;252(5013):1651–6.
Vera Cruz EM, Brown CL, Luckenbach JA, Picha ME, Bolivar RB, Borski RJ. Insulin-like growth factor-I cDNA cloning, gene expression and potential use as a growth rate indicator in Nile tilapia,< i > Oreochromis niloticus</i> Aquaculture. 2006;251(2):585–95.
Kajimura S, Uchida K, Yada T, Riley L, Byatt J, Collier R, et al. Stimulation of insulin-like growth factor-I production by recombinant bovine growth hormone in Mozambique tilapia, Oreochromis mossambicus. Fish Physiol Biochem. 2001;25(3):221–30.
Zhang D-C, Huang Y-Q, Shao Y-Q, Jiang S-G. Molecular cloning, recombinant expression, and growth–promoting effect of mud carp (< i > Cirrhinus molitorella</i>) insulin-like growth factor-I. Gen Comp Endocrinol. 2006;148(2):203–12.
Luo R, Liu B, Xie Y, Li Z, Huang W, Yuan J, et al. SOAPdenovo2: an empirically improved memory-efficient short-read de novo assembler. Gigascience. 2012;1(1):18.
Conesa A, Gotz S. Blast2GO: A comprehensive suite for functional analysis in plant genomics. Int J Plant Genomics. 2008;2008:619832.
Flight RM, Wentzell PD. Potential bias in GO::TermFinder. Brief Bioinform. 2009;10(3):289–94.
Boyle EI, Weng S, Gollub J, Jin H, Botstein D, Cherry JM, et al. GO::TermFinder--open source software for accessing Gene Ontology information and finding significantly enriched Gene Ontology terms associated with a list of genes. Bioinformatics. 2004;20(18):3710–5.
Abdi H. The Bonferonni and Šidák corrections for multiple comparisons. Encycl Meas Stat. 2007;3:103–7.
Shull GH. The composition of a field of maize. J Hered. 1908;1:296–301.
Zhang HY, He H, Chen LB, Li L, Liang MZ, Wang XF, et al. A genome-wide transcription analysis reveals a close correlation of promoter INDEL polymorphism and heterotic gene expression in rice hybrids. Mol Plant. 2008;1(5):720–31.
Clark RG, Robinson IC. Up and down the growth hormone cascade. Cytokine Growth Factor Rev. 1996;7(1):65–80.
Li WS, Lin HR, Wong AO. Effects of gonadotropin-releasing hormone on growth hormone secretion and gene expression in common carp pituitary. Comp Biochem Physiol B Biochem Mol Biol. 2002;132(2):335–41.
Wong AO, Leung MY, Shea WL, Tse LY, Chang JP, Chow BK. Hypophysiotropic action of pituitary adenylate cyclase-activating polypeptide (PACAP) in the goldfish: immunohistochemical demonstration of PACAP in the pituitary, PACAP stimulation of growth hormone release from pituitary cells, and molecular cloning of pituitary type I PACAP receptor. Endocrinology. 1998;139(8):3465–79.
Wong AO, Chang JP, Peter RE. Dopamine stimulates growth hormone release from the pituitary of goldfish, Carassius auratus, through the dopamine D1 receptors. Endocrinology. 1992;130(3):1201–10.
Lin XW, Lin HR, Peter RE. The regulatory effects of thyrotropin-releasing hormone on growth hormone secretion from the pituitary of common carp in vitro. Fish Physiol Biochem. 1993;11(1–6):71–6.
Hwa V, Oh Y, Rosenfeld RG. The insulin-like growth factor-binding protein (IGFBP) superfamily. Endocr Rev. 1999;20(6):761–87.
Firth SM, Baxter RC. Cellular actions of the insulin-like growth factor binding proteins. Endocr Rev. 2002;23(6):824–54.
Wood AW, Duan C, Bern HA. Insulin-like growth factor signaling in fish. Int Rev Cytol. 2005;243:215–85.
Stitt TN, Drujan D, Clarke BA, Panaro F, Timofeyva Y, Kline WO, et al. The IGF-1/PI3K/Akt pathway prevents expression of muscle atrophy-induced ubiquitin ligases by inhibiting FOXO transcription factors. Mol Cell. 2004;14(3):395–403.
Burgos SA, Cant JP. IGF-1 stimulates protein synthesis by enhanced signaling through mTORC1 in bovine mammary epithelial cells. Domest Anim Endocrinol. 2010;38(4):211–21.
Kalloo-Hosein HE, Whitehead JP, Soos M, Tavare JM, Siddle K, O’Rahilly S. Differential signaling to glycogen synthesis by the intracellular domain of the insulin versus the insulin-like growth factor-1 receptor. Evidence from studies of TrkC-chimeras. J Biol Chem. 1997;272(39):24325–32.
Bibollet-Bahena O, Almazan G. IGF-1-stimulated protein synthesis in oligodendrocyte progenitors requires PI3K/mTOR/Akt and MEK/ERK pathways. J Neurochem. 2009;109(5):1440–51.
Ma XM, Blenis J. Molecular mechanisms of mTOR-mediated translational control. Nat Rev Mol Cell Biol. 2009;10(5):307–18.
Kimball SR, Shantz LM, Horetsky RL, Jefferson LS. Leucine regulates translation of specific mRNAs in L6 myoblasts through mTOR-mediated changes in availability of eIF4E and phosphorylation of ribosomal protein S6. J Biol Chem. 1999;274(17):11647–52.
Kanehisa M, Goto S. KEGG: kyoto encyclopedia of genes and genomes. Nucleic Acids Res. 2000;28(1):27–30.
Zhang T, Wang S, Lin Y, Xu W, Ye D, Xiong Y, et al. Acetylation negatively regulates glycogen phosphorylase by recruiting protein phosphatase 1. Cell Metab. 2012;15(1):75–87.
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.
This work was financially supported by China 863 project (No. 2012AA10A407), Natural Science Foundation of China (No. 31572596; No. 31370047), Guangdong Provincial Natural Science Foundation (No. 2015A030313069), Special Fund for Fisheries-Scientific Research of Guangdong Province (No. A201400A01, A201501A03 and A201501A09), Fundamental Research Funds for the Central Universities (151gzs102, 151gzs121), Shenzhen and Hong Kong Innovation Circle Project (SGLH20131010105856414), Special Project on the Integration of Industry, Education and Research of Guangdong Province (2013B090800017), Shenzhen Special Program for Future Industrial Development (JSGG20141020113728803), and Shenzhen Dapeng Special Program for Industrial Development (KY20140104).
The authors declare that they have no competing interests.
YS contributed study design. XFL, CG, DW, XZ, CF and GH performed the bioinformatical analysis. DW and LX contributed to the fish and sample preparation. CG contributed to the RT-PCR experiments. YZ, QS, HL gave technical advice, and YS, CG, XFL, QS, XY and GH co-wrote the paper. All authors read and approved the final manuscript.
RPKMs of all genes detected in the brains and livers of three grouper species. (XLSX 4270 kb)
Differentially expressed genes between the brains of Hybrid F1 and maternal E. fuscoguttatus. (XLSX 87 kb)
Differentially expressed genes between the brains of Hybrid F1 and paternal E. lanceolatus. (XLSX 73 kb)
Differentially expressed genes between the livers of Hybrid F1 and maternal E. fuscoguttatus. (XLSX 460 kb)
Differentially expressed genes between the livers of Hybrid F1 and paternal E. lanceolatus. (XLSX 595 kb)
The hierarchical clustering map of DGEs among three species in the brain. Efu, Ela, and Hyb denote E. fuscoguttatus, E. lanceolatus and their hybrid F1, respectively. (PDF 40 kb)
The hierarchical clustering map of DGEs among three species in the liver. Efu, Ela, and Hyb denote E. fuscoguttatus, E. lanceolatus and their hybrid F1, respectively. (PDF 236 kb)
Functional annotation of grouper transcripts based on GO categorization. The left y-axis indicates the percentage of a specific category of genes in that main category. The right y-axis indicates the number of genes in a category. (TIF 521 kb)
GO enrichment of DEGs between the livers of Hybrid F1 and maternal E. fuscoguttatus. (XLSX 309 kb)
GO enrichment of DEGs between the livers of Hybrid F1 and paternal E. lanceolatus. (XLSX 464 kb)
GO enrichment of DEGs between the brains of Hybrid F1 and maternal E. fuscoguttatus. (XLSX 61 kb)
GO enrichment of DEGs between the brains of Hybrid F1 and paternal E. lanceolatus. (XLSX 51 kb)
Sequences of the primers used for real-time RT-PCR and related gene sequences. (XLSX 13 kb)
About this article
Cite this article
Sun, Y., Guo, C., Wang, D. et al. Transcriptome analysis reveals the molecular mechanisms underlying growth superiority in a novel grouper hybrid (Epinephelus fuscogutatus♀ × E. lanceolatus♂). BMC Genet 17, 24 (2016) doi:10.1186/s12863-016-0328-y
- Growth superiority
- GH/IGF axis
- Hybrid grouper