Printer Friendly

Comparison of gene expression patterns of key growth genes between different rate growths in zebrafish (Danio rerio) siblings.


Fish larvae stage is a developmental phase after the embryogenesis, in which the metamorphosis process generates a major morphological and physiological changes intrinsically associated with each species of fish (Dufour et al, 2012). From a point of view of cohort growth, at hatching the larvae showed similar body length, however within a short period of time is possible observed size heterogeneity among individuals; due a natural phenomenon called growth rate variability (Deangelis & Coutant, 1979; Kestemont et al., 2003). Growth rate variability is a natural phenomenon, which influencing individual performance either from an ecological perspective (Pepin et al., 2015) or in a fish aquaculture productivity (Goldan et al., 1997; Lekang, 2013). In fish farming, size heterogeneity is not ideal, because the subordinate fish have less access to the feed, have more stress and increase the possibility of cannibalism (Lekang, 2013). Size heterogeneity in the cohorts are expressed mathematically as the coefficient of variation (CV) (Weiner & Solbrig, 1984).

The physiological process of growing is mainly regulated by the growth hormone (GH)/insulin-like growth factor 1 (IGF-1) axis (Reinecke et al., 2005). Which is a pleiotropic physiological axis of hormones and cellular receptors that regulate: nutrients metabolism, protein synthesis in muscle and in general tissue growth and osmotic balance (Butler & Le Roith, 2001; Reinecke et al., 2005). The endocrine activity of GH has two pathways: the direct, in which the physiological effects are mediated by GH binding to its receptor (GHR), and indirect pathway, in which GH induces IGF-1 secretion and promotes biological activities (Canosa et al., 2007). Among insulin growth factors: iGf-1, IGF-2 and iGf-3; IGF-1 is the main hormone in the regulation of larvae growth, because IGF-2 expression is only observed at embryogenesis (Wood et al., 2005) and IGF-3 is gonad specific (Wang et al., 2008). How a counterpart of hormones are the growth axis receptors: the GHR is a single transmembrane glycoprotein that belongs to the class I cytokine receptor superfamily (Perez-Sanchez et al., 2002), and the IGF-1 receptor belongs to the tyrosine kinase superfamily of transmembrane receptor like the insulin (Wood et al., 2005). On the other hand, insulin-like growth factor binding protein (IGFBP) is a protein family of six members can bind IGF-1 and IGF-2, their roles are to increase the half-life of IGFs and distribution, which also are been described in fish (Daza et al., 2011; Reindl & Sheridan, 2012), the most abundantly in plasma are IGFBP-1, IGFBP-2 and IGFBP-3 (Wood et al., 2005).

Furthermore many peptides or hormones stimulating or inhibiting the GH/IGF-1 axis (Canosa et al., 2007; Chang & Wong, 2009), and affect the growth. The thyroid hormones triiodothyronine ([T.sub.3]) and thyroxine ([T.sub.4]) have stimulatory effects on GH and IGF-1, and they are the principal factors controlling metamorphosis in fish larvae (Wang & Zhang, 2011; McMenamin & Parichy, 2013). Ghrelin is a peptide hormone mainly secreted in the oxyntic mucosa of the stomach and is the ligand for the growth hormone secretagogue receptor (GHS-R); which is the other endocrine pathway that stimulates the secretion of GH by the pituitary gland in addition to the GH-releasing hormone (Dimaraki & Jaffe, 2006). Additionally, ghrelin is an orexigenic factor that increases food intake and plays an important role in energy and glucose homeostasis (Peter & Chang, 1999; Nakazato et al., 2001; Unniappan & Peter, 2005; Dimaraki & Jaffe, 2006; Arcamone et al., 2009; Pradhan et al., 2013). Cholecystokinin (CCK) is a peptide hormone secreted by the gastrointestinal tract, and its effects include gallbladder and pancreatic secretion, gastric and intestinal motor function, reduced food intake and stimulation of GH secretion (Canosa et al., 2007; Crespo et al., 2014; Micale et al., 2014; Dalmolin et al., 2015).

Many reports have observed that fasting or poor nutritional status in fish alters the mRNA expression of components of the GH/IGF-1 axis, such as: starving Oncorhynchus kisutch and Lates calcarifer (Duan & Plisetskaya, 1993; Matthews et al., 1997) or in fasting Anguilla japonica, Dicentrarchus labrax, Ictalurus punctatus and Oncorhynchus mykiss (Duan & Hirano, 1992; Norbeck et al., 2007; Terova et al., 2007; Peterson et al., 2009).

Growth rate variability has been attributed to both biotic and abiotic mechanisms, which have been categorized as either "imposed" or "inherent" (Huston & DeAngelis, 1987; Kestemont et al., 2003). Imposed mechanisms include: temperature, day length, food availability, and interactive factors such as food competition. On the other hand, inherent mechanisms have strong genetic influence, so high cohort genetic variability increases the growth rate variability (Nicieza et al., 1994; Hutchings & Jones, 1998; Ohlberger et al., 2013). Minimal information is available regarding gene expression patterns of growth factors due the growth rate variability process, as well is necessary define a main mechanism and tested in an isolated fashion. Food competition was the main mechanism proposed for this study, which may influence a poor nutritional status in some individuals; hence larvae growth rate may be associated with key growth gene expressions patters. Our aim was assess the key growth gene patterns observed at growth rate variability process. This study was carried out in environmentally controlled conditions (for the control inherent-non interactive mechanism) and to minimize genetic variability in the zebrafish larvae cohort, we used siblings (for the control imposed genetic mechanism).


Experimental animals

From a spawning with one pair of adult zebrafish we obtained 100 viable eggs, which were incubated at 26[degrees]C; only those larvae that hatched between 48 to 72 h after spawning were included. In the experimental design proposed, the use hatching siblings sought to reduce the cohort genetic variability. Larvae, were maintained in glass flasks with 2 L of E2 methylene blue media (Westerfield, 2000) under controlled light/dark conditions (14L/10D), with a 30% of water change every day. At 6 days post-hatching (dph) or 156 accumulated thermal units (ATU), the body lengths of co-hatched zebrafish larvae were measured. The larvae standard length measurement was conducted under stereoscopic microscope with a Motic[R] Images Plus 2.0ML software according to the proceeding proposed by Parichy et al. (2009); the larvae were previously anesthetized by benzocaine 20% (0.2 mL [L.sup.-1]). Subsequently, total RNA was isolated from 30 zebrafish larvae (4 mm body length) to establish the initial state of gene expression (reference group); these larvae were grouped into five pools or biological replicates of six larvae each. The other co-hatching siblings larvae were reared for 14 days and were fed with rotifers (Brachionus plicatilis), at the rate of 200 rotifers per larva per day (Lawrence et al., 2012). At 20 dph, the body length of each remaining larva was measured, and then 60 larvae were classified into two groups: large larvae and small larvae. The small larvae group was composed of 30 individuals with a body length of 4 mm, and they were divided into five pools or biological replicates with 6 larvae each. The large larvae group was also composed of 30 individuals divided into five pools or biological replicates of 6 larvae each. The large larvae groups were organized in its body length as follow: 4 larvae groups of 5 mm and one larvae group of 6 mm. This study was conducted in strict accordance with the recommendations in the Guide for the Care and Use of Laboratory Animals of the National Institutes of Health. The protocol was approved by the Committee on the Ethics of Animal Experiments of the INTA Universidad de Chile.

Total RNA isolation and reverse transcription

All zebrafish larvae pools (n = 15, small-large reference) were placed in 1.7 mL microcentrifuge tubes; excess liquid was removed, and euthanasia was performed by freezing at -80[degrees]C in liquid nitrogen. Total RNA of all pools was isolated using 800 [micro]L of Tripure[R] reagent (Roche) (Lan et al., 2009), according to the manufacturer's instructions. RNA was quantified with a spectrophotometry at 260 and 280 nm (Nano-Drop[R]) and the RNA quality was assessed with 1% agarose gel electrophoresis. They were treated with RQ1 RNase-Free DNase (cat. M6101, Promega[R]) to avoid genomic DNA amplification, the absence of genomic DNA was confirmed by PCR on the treated RNA. The first-strand cDNA synthesis was performed using the ImProm-II[TM] Reverse Transcription System (Promega[R]). Total RNA was combined with 0.5 [micro]g [reaction.sup.-1] Oligo[(dT).sub.15] Primer (cat. C1101, Promega[R]) for a final volume of 5 pL and incubated at 70[degrees]C for 5 min. Next, 15 [micro]L of the transcription mix (ImProm-II[TM] 5X Reaction Buffer 4.6 pL, 2.25 mM of Mg[Cl.sub.2], 0.5 mM each dNTP, Recombinant RNasin[R] Ribonuclease Inhibitor (Promega[R]) 20 [micro]L and 1 [micro]L ImProm-II[TM] Reverse Transcriptase (Promega[R])) was added. Following the addition of transcription mix, the reaction was maintained at 25[degrees]C for 5 min and then transferred to 42[degrees]C for 60 min. The reverse transcription reactions were stopped by heating the mixture at 70[degrees]C for 15 min.

qPCR analysis

The gene-specific oligonucleotide primers for growth hormone (gh), growth hormone receptor (ghra), insulin-like growth factor 1a (igf-1a), insulin-like growth factor receptor a and b (igf-1r a and b), insulin-like growth factor binding protein 1 (igfbp-1), ghrelin (ghrl), cholecystokinin a (cck) and glucokinase (gck) were developed using Primer-BLAST (NCBI) (Ye et al., 2012). To test the modulation of T3 and T4 in fish larvae, we assessed the gene expression of their precursor protein, thyroglobulin (tg). For normalization of cDNA loading, all samples were run in parallel using the housekeeping gene elongation factor I-alpha (efl[alpha]) as the reference gene (McCurley & Callard, 2008). All primers are listed in Table 1. The relative mRNA expression levels of target genes and the reference gene (efl[alpha]) were quantified using real-time PCR analysis with AriaMx Real-Time PCR (Agilent Technologies).

Amplification of specific PCR products was detected using the FastStart Essential DNA Green Master[R] (Roche), according to the manufacturer's instructions. All cDNA examples were analyzed in duplicate. The amplification protocol used was as follows: one initial step of 10 min at 95[degrees]C (denaturation and enzyme activation), followed by 45 cycles of 95[degrees]C for 10 s, 60[degrees]C for 5 s and 72[degrees]C for 15 s. After the amplification, melting curve analysis was performed over a range of 50-95[degrees]C to verify that a single PCR product was generated at the end of the assay.

Data & statistical analysis

The cohort's coefficients of variation were calculated based on the formula proposed by Sokal & Rohlf (1995):

Cv = (1 + 1/4xn) x (Sx100/x)

where n is the number of observations; 5 is the sample standard deviation, and x is the sample mean. Density histograms were made using the program R-3.1.2 for Windows (32/64 bit) (R Core Team, 2014).

The relative expression levels of the genes were calculated by the method of Pfaffl (2001), using the reference group as a control group in the equation. The primer PCR efficiency (E) was calculated for each gene fluorescence curve with LigRegPCR 12.18 software (Udvardi et al., 2008), and the efficiency rates for the transcripts were as follows: 1.96 for gh, 1.87 for igf-1, 1.91 for tg, 1.8 for igfbp, 1.82 for igf-1r(a/b), 1.7 for gck, 1.88 cck, 1.96 ghrelin and 1,84 for efl[alpha] over the entire quantification range. The differences in the gene expression levels were analyzed by a Wilcox-Mann-Whitney test (Derveaux et al., 2010) between the small and large pools using R-3.1.2 for Windows (32/64 bit) (R Core Team, 2014), P-values <0.05 were considered significant. In addition, gene expression was analyzed by principal component analysis (PCA) (Abdi & Williams, 2010). The principal component analyses were made using the FactoMineR packages and the biplot by Factoextra and ggplot2 packages in R-3.1.2 (Ringner, 2008; R Core Team, 2014).


Body size heterogeneity

The distributions of larval body length at 6 dph and 20 dph are presented in the (Figs. 1a-1b). The larvae body length distribution observed in the beginning of the study (6 dph), showed a mode associated with 4.0 mm in body length. This mode represented nearly 72% of the measured larvae, and the remaining larvae were shorter than 4 mm. At 20 dph, body length distribution evolved to show two modes, 4.0 mm (46%) and 5.0 mm (44%), and 10% of the larvae were 6 mm long, representing the longest larvae. The coefficient of variation between the cohorts changed from 5.79 at 6 dph to 14.22 at 20 dph. The mortality during the rearing period was 8%.

Modulation of genes related to growth

In the GH/IGF-1 axis, differences in gene expression between small and large larvae are presented in Fig. 2. Only were statistically significant for the igf-1 (a/b) receptor (Wilcox-Mann-Whitney test, P = 0.02) and igflbp (Wilcox-Mann-Whitney test, P = 0.01), the other genes of growth axis were not statistically significant (P > 0.05) and showed similar levels of gene expression between the groups. The expression levels of the other genes evaluated (tg, ghrelin, cck, gck) were not statistically significant (Fig. 3), although for gck, the Wilcox-Mann-Whitney test yielded P = 0.07.

Principal components analysis (PCA)

This descriptive analysis revealed other aspects of variability among individuals (the larvae pools) and variables (genes), enabling us to elucidate the main components of the variability. The PCA results revealed that approximately 84.4% of the inertia was explained by four components or dimensions (PCs): PC1 = 26.6%, PC2 = 23.6%, PC3 = 21% and PC4 = 13.2%. The PCA with PC1 and PC2 (Figs. 4a-4b). Figure 4a shows the correlation among variables (gene expression). There is a positive correlation between gck and igf-1r, as well as between tg and igfbp-1. A projection in the first component (PC1) is apparent in both groups and the two projections are negatively correlated. Figure 4b shows the individuals (larvae pools) and circumscribed groups associated with larvae body length. This aggrupation is projected onto PC1, the large larvae are on the negative side and the small larvae on the positive side. The correlations among components and variables are presented in Table 2, which shows that the most important genes correlated with PC1 were igfbp-1 (0.86, P = 0.001), gck (-0.67, P = 0.03) and tg (0.65, P = 0.03). Therefore, gck was significantly correlated with PC1 and its projection was associated with large larvae, while tg and igfbp-1 were significantly correlated with PC1, although their projections were associated with small larvae (Figs. 4a-4b). In PC2, the most important correlations were cck (0.92, P = 0.0001) and ghr (0.66, P = 0.03).


The present study confirms that the growth rate variability phenomenon modulate the gene expression patterns of growth endocrine control genes. The siblings zebrafish cohort presented a growth rate variability after the rearing period, because the larvae body length distribution began with low CV (5.8%) at 6 dph and was raised to 14.22% at 20 dph. The size length variation observed in the study is according to the size variation observed in other studies in Sparus aurata (Goldan et al., 1997) or in Sciaenops ocellatus (Smith & Fuiman, 2003) and this situation is consistent with the report of Kestemont et al. (2003) in the growth rate variability phenomenon.

Larvae body length groups not showed significant differences in gh and igf-1 mRNA levels. However, the mean of gh mRNA levels was slightly higher in smaller larvae than in large larvae, this trend conforms to the expectations associated with the poor nutritional status (Wood et al., 2005; Norbeck et al., 2007; Savage, 2013). Conversely, was observed uniformity in igf-1 mRNA levels between the larvae groups, in contrast to the most fasting reports in different fish species (Wood et al., 2005; Norbeck et al., 2007; Peterson & Waldbieser, 2009; Reinecke, 2010; Kawaguchi et al., 2013; Tian et al., 2015; Taniyama et al., 2016). However, Wen-Ying et al. (2012) in Carassius auratus gibelio, Breves et al. (2014) and Fox et al. (2010) in Mozambique tilapia (Oreochromis mossambicus), and Hevroy et al. (2011) in Atlantic salmon (Salmon salar) only observed significant differences at the protein levels but not at the mRNA levels by fasting challenge. Less information is available about the ghr and igf-1r on teleost fish, regarding to the nutritional regulation. Regard ghr mRNA levels modulation by fasting challenge, the fish studies showed variable results. In channel catfish (Ictalurus punctatus) by increasing feeding levels, ghr mRNA levels did not change significantly between the groups (Peterson et al., 2008). However in rainbow trout (Oncorhynchus mykiss) the fasting reduce ghr the mRNA levels (Walock et al., 2014); conversely fasting increase ghra and ghrb mRNA levels in zebrafish (Danio rerio) (Tian et al., 2015). Regard to igf-1r mRNA levels the fish studies also have been showed variable results. As well, the in Ictalurus punctatus and Oncorhynchus mykiss, it did not show mRNA levels modulation by fasting (Peterson et al., 2009, Gabillard et al., 2003); though Norbeck et al. (2007) found igf-1r up-regulation in gill but not in skeletal muscle in fasting Oncorhynchus mykiss. Nevertheless, this studies analyses igf-1r mRNA levels by specific tissue, nonetheless this receptor is expressed in all body tissues (Wood et al., 2005; Nimptsch & Giovannucci, 2012); therefore, these results give a partial interpretation of the igf-1r expression. Conversely, our study used the complete larvae and the igf-1r gene expression adds all body tissues. Hence, the significant difference between the body length groups in igf-1r suggests that the gene expression of igf-1r may be more sensible by the larvae nutritional status than its main ligand igf-1.

In fish fasting challenge influence the modulation of gastrointestinal peptides as cholecystokinin and ghrelin. The cck mRNA levels decrease by fasting in different fish species (Murashita et al., 2006; Feng et al., 2012; Ji et al., 2015), conversely ghrelin increase mRNA levels by fasting (Amole & Unniappan, 2009; Zhou & Xue, 2009; Tian et al., 2015; Volkoff, 2015; Blanco et al., 2016). These gene expression modulations were not observed in our results.

The most likely explanation of the key growth gene expression patterns results could be associated to the larvae nutritional status. The igfbp-1 and gck gene expressions are regulated by nutritional status or the glucose levels. The igfbp-1 is regulated by insulin levels and in consequently with glucose levels (Lee et al., 1993) and fasting increase its mRNA levels in fish (Shimizu et al., 2006; Hevroy et al., 2011; Kawaguchi et al., 2013; Breves et al., 2014), according to the our study results. On the other hand, gck is a liver enzyme that catalyzes the phosphorylation of glucose to glucose-6-phosphate (Enes et al., 2009), is associated with individual nutritional status and showing an upregulation by feed intake (Caseras et al., 2000; Gonzalez-Alvarez et al., 2009; Panserat et al., 2014). Likewise, the main significant correlations in PCA analysis were according with this interpretation. The large larvae could be associated with a high food intake or better nutritional status and this increase gck mRNA levels; this is according with the high significant correlation observed between these factors. On the other hand, the igfbp-1 was correlated with small larvae; this outcome was according with the igfbp-1 modulation by poor nutritional status. Hence, the results suggest that the growth rate variability was associated with the larvae nutritional status, and possibly influenced by food competition (Ruzzante, 1994). However, intensity level of the nutrient restriction was not similar than fasting challege, because the food competition does not prevent small larvae from feeding; rather, they have less access to food. As well, this light fasting condition could prevent observed significant differences in gck or igf-1 mRNA levels between the body length groups.

In conclusion, our results suggest that growth rate variability affect the gene expression of igfbp-1 and igf-1r genes, the increment in igfbp-1 mRNA levels observed in the small larvae suggest that nutritional status is associated to their growth rate. Future researches have to include protein assess and different feed levels to understand the growth rate variability influence in larvae zebrafish growth rate.

DOI: 10.3856/vol45-issue4-fulltext-12


This investigation was supported by a grant (FONDE CYT Post-Doctorado No. 3130518) from CONICYTChile and Concurso Nacional de Insercion de Capital Humano Avanzado en la Academia No. 79110002 from CONICYT-Chile.


Abdi, H. & L.J. Williams. 2010. Principal component analysis. Wiley interdisciplinary reviews: computational statistics, 2: 433-459.

Amole, N. & S. Unniappan. 2009. Fasting induces preproghrelin mRNA expression in the brain and gut of zebrafish, Danio rerio. Gen. Comp. Endocr., 161: 133-137.

Arcamone, N., S. Neglia, G. Gargiulo, V. Esposito, E. Varricchio, P. Battaglini, P. De Girolamo & F. Russo. 2009. Distribution of ghrelin peptide in the gastrointestinal tract of stomachless and stomach-containing teleosts. Microsc. Res. Tech., 72: 525-533.

Blanco, A.M., M. Gomez-Boronat, I. Redondo, A.I. Valenciano & M.J. Delgado. 2016. Periprandial changes and effects of short- and long-term fasting on ghrelin, GOAT, and ghrelin receptors in goldfish (Carassius auratus). J. Comp. Physiol. B, 186: 727-738.

Breves, J.P., C.K. Tipsmark, B.A. Stough, A.P. Seale, B.R. Flack, B.P. Moorman, D.T. Lerner & E.G Grau. 2014. Nutritional status and growth hormone regulate insulin-like growth factor binding protein (igfbp) transcripts in Mozambique tilapia. Gen. Comp. Endocr., 207: 66-73.

Butler, A.A. & D. Le Roith. 2001. Control of growth by the somatropic axis: growth hormone and the insulin-like growth factors have related and independent roles. Ann. Rev. Physiol., 63: 141-164.

Canosa, L.F., J.P. Chang & R.E. Peter. 2007. Neuroendocrine control of growth hormone in fish. Gen. Comp. Endocrin., 151: 1-26.

Caseras, A., I. Meton, F. Fernandez & I.V. Baanante. 2000. Glucokinase gene expression is nutritionally regulated in liver of gilthead sea bream (Sparus aurata). BBA-Gene Struct. Express., 1493: 135-141.

Crespo, C.S., A.P. Cachero, P. Jimenez, V. Barrios & E. Arilla. 2014. Peptides and food intake. Front. Endocrinol., 5(58): 1-13.

Chang, J.P. & A.O.L Wong. 2009. Growth hormone regulation in fish: a multifactorial model with hypothalamic, peripheral and local autocrine/paracrine signals. In: J. Nicholas, D.G. Bernier, N. Bernier, G. Van der Kraak, A. Farrell & C. Brauner (eds.). Fish physiology. Academic Press, New York, pp. 151-195.

Dalmolin, C., D. Almeida, M. Figueiredo & L. Marins. 2015. Food intake and appetite control in a GHtransgenic zebrafish. Fish Physiol. Biochem., 5: 1131-1141.

Daza, D.O., G. Sundstrom, C.A. Bergqvist, C. Duan & D. Larhammar. 2011. Evolution of the insulin-like growth factor binding protein (IGFBP) family. Endocrinology, 152: 2278-2289.

Deangelis, D.L. & C.C. Coutant. 1979. Growth rates and size distributions of first-year smallmouth bass populations: some conclusions from experiments and a model. T. Am. Fish. Soc., 108: 137-141.

Derveaux, S., J. Vandesompele & J. Hellemans. 2010. How to do successful gene expression analysis using real-time PCR. Methods, 50: 227-230.

Dimaraki, E.V. & C.A. Jaffe. 2006. Role of endogenous ghrelin in growth hormone secretion, appetite regulation and metabolism. Rev. Endocr. Metab. Disord., 7: 237-249.

Duan, C. & T. Hirano. 1992. Effects of insulin-like growth factor-I and insulin on the in-vitro uptake of sulphate by eel branchial cartilage: evidence for the presence of independent hepatic and pancreatic sulphation factors. J. Endocrinol., 133: 211-219.

Duan, C. & E.M. Plisetskaya. 1993. Nutritional regulation of insulin-like growth factor-I mRNA expression in salmon tissues. J. Endocrinol., 139: 243-252.

Dufour, S., K. Rousseau & B.G. Kapoor. 2012. Metamorphosis in fish. Taylor & Francis Group, New York, 268 pp.

Enes, P., S. Panserat, S. Kaushik & A. Oliva-Teles. 2009. Nutritional regulation of hepatic glucose metabolism in fish. Fish Physiol. Biochem., 35: 519-539.

Feng, K., G.R. Zhang, K.J. Wei, B.X. Xiong, T. Liang & H.C. Ping. 2012. Molecular characterization of cholecystokinin in grass carp (Ctenopharyngodon idellus): cloning, localization, developmental profile, and effect of fasting and refeeding on expression in the brain and intestine. Fish Physiol. Biochem., 38: 1825-1834.

Fox, B.K., J.P. Breves, L.K. Davis, A.L. Pierce, T. Hirano & E.G. Grau. 2010. Tissue-specific regulation of the growth hormone/insulin-like growth factor axis during fasting and re-feeding: Importance of muscle expression of IGF-I and IGF-II mRNA in the tilapia. Gen. Comp. Endocr., 166: 573-580.

Gabillard, J.C., C. Weil, P.Y. Rescan, I. Navarro, J. Gutierrez & P.Y. Le Bail. 2003. Effects of environmental temperature on IGF1, IGF2, and IGF type I receptor expression in rainbow trout (Oncorhynchus mykiss). Gen. Comp. Endocrinol., 133: 233-242.

Goldan, O., D. Popper & I. Karplus. 1997. Management of size variation in juvenile gilthead sea bream (Sparus aurata). 1. Particle size and frequency of feeding dry and live food. Aquaculture, 152: 181-190.

Gonzalez-Alvarez, R., D. Ortega-Cuellar, A. Hernandez-Mendoza, E. Moreno-Arriola, K. Villasenor-Mendoza, A. Galvez-Mariscal, M.E. Perez-Cruz, I. Morales-Salas & A. Velazquez-Arellano. 2009. The hexokinase gene family in the zebrafish: structure, expression, functional and phylogenetic analysis. Comp. Biochem. Physiol. B, 152: 189-195.

Hevroy, E.M., C. Azpeleta, M. Shimizu, A. Lanzen, H. Kaiya, M. Espe & P.A. Olsvik. 2011. Effects of short-term starvation on ghrelin, GH-IGF system, and IGF-binding proteins in Atlantic salmon. Fish Physiol. Biochem., 37: 217-232.

Huston, M.A. & D.L. DeAngelis. 1987. Size bimodality in monospecific populations: a critical review of potential mechanisms. Am. Natur., 129: 678-707.

Hutchings, J.A. & M.E.B. Jones. 1998. Life history variation and growth rate thresholds for maturity in Atlantic salmon, Salmo salar. Can. J. Fish. Aquat. Sci., 55: 22-47.

Ji, W., H.-C. Ping, K.-J. Wei, G.-R. Zhang, Z.-C. Shi, R. B. Yang, G.-W. Zou & W.-M. Wang. 2015. Ghrelin, neuropeptide Y (NPY) and cholecystokinin (CCK) in blunt snout bream (Megalobrama amblycephala): cDNA cloning, tissue distribution and mRNA expression changes responding to fasting and refeeding. Gen. Comp. Endocr., 223: 108-119.

Kawaguchi, K., N. Kaneko, M. Fukuda, Y. Nakano, S. Kimura, A, Hara & M. Shimizu. 2013. Responses of insulin-like growth factor (IGF)-I and two IGF-binding protein-1 subtypes to fasting and re-feeding, and their relationships with individual growth rates in yearling masu salmon (Oncorhynchus masou). Comp. Biochem. Physiol. A, 165: 191-198.

Kestemont, P., S. Jourdan, M. Houbart, C. Melard, M. Paspatis, P. Fontaine, A. Cuvier, M. Kentouri & E. Baras. 2003. Size heterogeneity, cannibalism and competition in cultured predatory fish larvae: biotic and abiotic influences. Aquaculture, 227: 333-356.

Lan, C.-C., R. Tang, I. Un San Leong & D.R. Love. 2009. Quantitative real-time RT-PCR (qRT-PCR) of zebrafish transcripts: optimization of RNA extraction, quality control considerations, and data analysis. Cold Spring Harbor Protocols, 4: 1-12.

Lawrence, C., E. Sanders & E. Henry. 2012. Methods for culturing saltwater rotifers (Brachionus plicatilis) for rearing larval zebrafish. Zebrafish, 9: 140-146.

Lee, P.D.K., C.A. Conover & D.R. Powell. 1993. Regulation and function of insulin-like growth factor-binding protein-1. Proc. Soc. Exp. Biol. Med., 204: 429.

Lekang, O.I. 2013. Internal transport and size grading. In: O.-I. Lekang (ed.). Aquaculture engineering. John Wiley & Sons, Oxford, 403 pp.

Matthews, S.J., A.K.K. Kinhult, P. Hoeben, V.R. Sara & T.A. Anderson. 1997. Nutritional regulation of insulin-like growth factor-I mRNA expression in barramundi, Lates calcarifer. J. Molec. Endocrinol., 18: 273-276.

McCurley, A.T. & G.V. Callard. 2008. Characterization of housekeeping genes in zebrafish: male-female differences and effects of tissue type, developmental stage and chemical treatment. BMC Molec. Biol., 9: 1-12.

McMenamin, S.K. & D.M. Parichy. 2013. Metamorphosis in teleosts. In: S. Yun-Bo, (ed.). Current topics in developmental biology. Academic Press, New York, pp. 127-165.

Micale, V., S. Campo, A. D'Ascola, M.C. Guerrera, M.B. Levanti, A. Germana & U. Muglia. 2014. Cholecystokinin: how many functions? Observations in seabreams. Gen. Comp. Endocrinol., 205: 166-167.

Murashita, K., H. Fukada, H. Hosokawa & T. Masumoto. 2006. Cholecystokinin and peptide Y in yellowtail (Seriola quinqueradiata): Molecular cloning, real-time quantitative RT-PCR, and response to feeding and fasting. Gen. Comp. Endocr., 145: 287-297.

Nakazato, M., N. Murakami, Y. Date, M. Kojima, H. Matsuo, K. Kangawa & S. Matsukura. 2001. A role for ghrelin in the central regulation of feeding. Nature, 409: 194-198.

Nicieza, A.G., F.G. Reyes-Gavilan & F. Brana. 1994. Differentiation in juvenile growth and bimodality patterns between northern and southern populations of Atlantic salmon (Salmo salar L.). Can. J. Zool., 72: 1603-1610.

Nimptsch, K. & E. Giovannucci. 2012. Epidemiology of IGF-1 and Cancer. In: D. LeRoith (ed.). Insulin-like growth factors and cancer. Springer, Neew York, pp. 1-24.

Norbeck, L.A., J.D. Kittilson & M.A. Sheridan. 2007. Resolving the growth-promoting and metabolic effects of growth hormone: differential regulation of GH-IGF-I system components. Gen. Comp. Endocrin., 151: 332-341.

Ohlberger, J., J. Otero, E. Edeline, I.J. Winfield, N.C. Stenseth & L.A. Vollestad. 2013. Biotic and abiotic effects on cohort size distributions in fish. Oikos, 122: 835-844.

Panserat, S., N. Rideau & S. Polakof. 2014. Nutritional regulation of glucokinase: a cross-species story. Nutr. Res. Rev., 27: 21-47.

Parichy, D.M., M.R. Elizondo, M.G. Mills, T.N. Gordon & R.E. Engeszer. 2009. Normal table of postembryonic zebrafish development: staging by externally visible anatomy of the living fish. Dev. Dynam., 238: 2975-3015.

Pepin, P., D. Robert, C. Bouchard, J.F. Dower, M. Falardeau, L. Fortier, G.P. Jenkins, V. Leclerc, K. Levesque, J.K. Llopiz, M.G. Meekan, H.M. Murphy, M., Ringuette, P. Sirois & S. Sponaugle. 2015. Once upon a larva: revisiting the relationship between feeding success and growth in fish larvae. ICES J. Mar. Sci., 72: 359-373.

Perez-Sanchez, J., J.A. Calduch-Giner, M. Mingarro, S. Vega-Rubin de Celis, P. Gomez-Requeni, A. Saera-Vila, A. Astola & M. Valdivia. 2002. Overview of fish growth hormone family. New insights in genomic organization and heterogeneity of growth hormone receptors. Fish Physiol. Biochem., 27: 243-258.

Peter, R.E. & J.P. Chang. 1999. Brain regulation of growth hormone secretion and food intake in fish. In: P.D. Prasada-Rao & R.E. Peter (eds.). Neural regulation in the vertebrate endocrine system. Springer Business Media, New York, pp. 55-67.

Peterson, B.C., A.L. Bilodeau-Bourgeois & B.C. Small. 2009. Response of the somatotropic axis to alterations in feed intake of channel catfish (Ictalurus punctatus). Comp. Biochem. Physiol. A, 153: 457-463.

Peterson, B.C., B.C. Small, G.C. Waldbieser & B.G. Bosworth. 2008. Endocrine responses of fast- and slow-growing families of channel catfish. N. Am. J. Aquacult., 70: 240-250.

Peterson, B.C. & G.C. Waldbieser. 2009. Effects of fasting on IGF-I, IGF-II, and IGF-binding protein mRNA concentrations in channel catfish (Ictalurus punctatus). Dom. Anim. Endocrinol., 37: 74-83.

Pfaffl, M.W. 2001. A new mathematical model for relative quantification in real-time RT-PCR. Nucleic Acids Res., 29: 2002-2007.

Pradhan, G., S.L. Samson & Y.X. Sun. 2013. Ghrelin: much more than a hunger hormone. Curr. Opin. Clin. Nutr., 16: 619-624.

R. Core Team. 2014. R: a language and environment for statistical computing. R Foundation for Statistical Computing, Vienna, Austria. []. Reviewed: 12 November 2016.

Reindl, K.M. & M.A. Sheridan. 2012. Peripheral regulation of the growth hormone-insulin-like growth factor system in fish and other vertebrates. Comp. Biochem. Physiol. A, 163: 231-245.

Reinecke, M. 2010. Influences of the environment on the endocrine and paracrine fish growth hormone-insulin-like growth factor-I system. J. Fish. Biol., 76: 1233-1254.

Reinecke, M., B.T. Bjornsson, W.W. Dickhoff, S.D. McCormick, I. Navarro, D.M. Power & J. Gutierrez. 2005. Growth hormone and insulin-like growth factors in fish: where we are and where to go. Gen. Comp. Endocrinol., 142: 20-24.

Ringner, M. 2008. What is principal component analysis? Nat. Biotech., 26: 303-304.

Ruzzante, D.E. 1994. Domestication effects on aggressive and schooling behavior in fish. Aquaculture, 120: 1-24.

Savage, M.O. 2013. Insulin-like growth factors, nutrition and growth. World Rev. Nutr. Diet., 106: 52-59. Shimizu, M., B.R. Beckman, A. Hara & W.W. Dickhoff.

2006. Measurement of circulating salmon IGF binding protein-1: assay development, response to feeding ration and temperature, and relation to growth parameters. J. Endocrinol., 188: 101-110.

Smith, M.E. & L.A. Fuiman. 2003. Causes of growth depensation in red drum, Sciaenops ocellatus, larvae. Environ. Biol. Fishes, 66: 49-60.

Sokal, R.R. & F.J. Rohlf. 1995. Biometry: the principles and practice of statistics in biological research. W.H. Freeman and Co, New York, 887 pp.

Taniyama, N., N. Kaneko, Y. Inatani, Y. Miyakoshi & M. Shimizu. 2016. Effects of seawater transfer and fasting on the endocrine and biochemical growth indices in juvenile chum salmon (Oncorhynchus keta). Gen. Comp. Endocr., 236: 146-156.

Terova, G., S. Rimoldi, V. Chini, R. Gornati, G. Bernardini & M. Saroglia. 2007. Cloning and expression analysis of insulin-like growth factor I and II in liver and muscle of sea bass (Dicentrarchus labrax, L.) during long-term fasting and refeeding. J. Fish Biol., 70: 219-233.

Tian, J., G. He, K.S. Mai & C.D. Liu. 2015. Effects of postprandial starvation on mRNA expression of endocrine-, amino acid and peptide transporter-, and metabolic enzyme-related genes in zebrafish (Danio rerio). Fish Physiol. Biochem., 41: 773-787.

Udvardi, M.K., T. Czechowski & W.R. Scheible. 2008. Eleven golden rules of quantitative RT-PCR. Plant Cell, 20: 1736-1737.

Unniappan, S. & R.E. Peter. 2005. Structure, distribution and physiological functions of ghrelin in fish. Comp. Biochem. Physiol. A, 140: 396-408.

Volkoff, H. 2015. Cloning, tissue distribution and effects of fasting on mRNA expression levels of leptin and ghrelin in red-bellied piranha (Pygocentrus nattereri). Gen. Comp. Endocr., 217: 20-27.

Walock, C.N., J.D. Kittilson & M.A. Sheridan. 2014. Characterization of a novel growth hormone receptor-encoding cDNA in rainbow trout and regulation of its expression by nutritional state. Gene, 533: 286-294.

Wang, Y. & S. Zhang. 2011. Expression and regulation by thyroid hormone (TH) of zebrafish IGF-I gene and amphioxus IGF-l gene with implication of the origin of TH/IGF signaling pathway. Comp. Biochem. Physiol. A, 160: 474-479.

Wang, D.-S., B. Jiao, C. Hu, X. Huang, Z. Liu & C.H.K. Cheng. 2008. Discovery of a gonad-specific IGF subtype in teleost. Biochem. Biophys. Res. Comm., 367: 336-341.

Weiner, J. & O. Solbrig. 1984. The meaning and measurement of size hierarchies in plant populations. Oecologia, 61: 334-336.

Wen-Ying., S., R. Gang & Z. Yao-Rong. 2012. Effects of compensatory growth on the levels of IGF-1, IGFBP-1 and expressions of IGF-1 mRNA, IGF-1R mRNA in Carassius auratus gibelio. Zool. Res., 33: 298-303.

Westerfield, M. 2000. The zebrafish book. A guide for the laboratory use of zebrafish (Danio rerio). University of Oregon Press, Eugene, 252 pp.

Wood, A.W., C.M. Duan & H.A. Bern. 2005. Insulin-like growth factor signaling in fish. Int. Rev. Cytol., 243: 215-285.

Ye, J., G. Coulouris, I. Zaretskaya, I. Cutcutache, S. Rozen & T.L. Madden. 2012. Primer-BLAST: A tool to design target-specific primers for polymerase chain reaction. BMC Bioinformatics, 13: 1-11.

Zhou, X.L. & C.R. Xue. 2009. Ghrelin inhibits the development of acute pancreatitis and nuclear factor kappa b activation in pancreas and liver. Pancreas, 38: 752-757.

Received: 18 August 2016; Accepted: 22 March 2017

Rafael Opazo (1), Luis Valladares (1) & Jaime Romero (1)

(1) Instituto de Nutricion y Tecnologia de los Alimentos (INTA), Universidad de Chile Corresponding author: Rafael Opazo (

Corresponding editor: Jesus Ponce-Palafox

Caption: Figure 1. Frequency histograms of body length (mm) zebrafish larvae cohort distribution a) at 6 dph and b) at 20 dph.

Caption: Figure 2. Relative gene expression of: a) Growth hormone (gh), b) Insulin-like growth factor 1 (igf-1), c) growth hormone receptor (ghra) and d) Insulin-like growth factor 1 receptor (igf-1ra/b), insulin-like growth factor binding protein 1 (igfbp-1) in the zebrafish larvae small and large body length groups assessment by qPCR. Bars represent mean [+ or -] SD, n = 5. The statistical significance was determined using the paired sample Wilcoxon Mann-Whitney U signed rank test (P < 0.05).

Caption: Figure 3. Relative gene expression of: a) Thyroglobulin (tg), b) Cholecystokinin (cck), c) Ghrelin and d) Glucokinase (gck) in the zebrafish larvae small and large body length groups assessment by RT-PCR. Bars represent mean [+ or -] and line the SD, n = 5. The statistical significance was determined using the paired sample Wilcoxon Mann-Whitney U signed rank test (P < 0.05).

Caption: Figure 4. PCA analysis a) PCA-Correlation loadings plot of the variables (gene expression) in the principal components PC1 and PC2, b) PCA-Score plot in the PC1 and PC2 principal components. growth hormone (gh), growth hormone receptor (ghr), insulin-like growth factor 1 (igf-1), insulin-like growth factor receptor (igf-1r), insulin-like growth factor binding protein 1 (igfbp-1), ghrelin (ghrl), cholecystokinin (cck) and glucokinase (gck), thyroglobulin (tg). The black individuals are the large body length zebrafish larvae pool and the red individuals are the small body length zebrafish larvae pool.
Table 1. Primers used for the quantification of the mRNA
expression by qPCR.

Target gene           Gene symbol      Genbank       Position
                                    accession no.

Cholecystokinin a        ccka       XM_001346104.4    222-302
Elongation factor     efl[alpha]     NM_131263.1     1414-1516
  1 alpha
Ghrelin                  ghrl       NM_001083872.1    23-158
Glucokinase               gck       NM_001045385.2   952-1119
Growth hormone            gh         NM_001123676     69-165
Growth hormone           ghra        NM_001083578     675-929
  receptor (a)
Insulin-like            igf-la       NM 131825.2      250-405
  growth factor 1
Insulin-like growth     Igfbp-1       NM_173283       575-716
  protein 1
Insulin-like growth   igf-1R(a/b)   NM_152968 and    3468-3619
  factor receptor                     NM_152969
  (a and b)
Thyroglobulin             tg         XM_689200.5     4274-4455
Target gene             Product     Sequence of primers
                      length (bp)   (5' [right arrow] 3')

Cholecystokinin a         80        (F) CGCCTGCTGGACAAATCAAC
                                    (R) GGCCAGTAGTTCGGTTAGGC
Elongation factor         103       (F) GTGCTGGCAAGGTCACAAAG
1 alpha                             (R) AGAGGTTGGGAAGAACACGC
Ghrelin                   136       (F) GCAGCATGTTTCTGCTCCTG
                                    (R) TCAGCAGCTTCTCTTCTGCC
Glucokinase               168       (F) ACGAGAAGCTGATTGGTGGG
                                    (R) TGTCCCCTGTGTCACTCTCA
Growth hormone            97        (F) CTGTTGCAGTTGGTGGTGGT
                                    (R) GGTGTTGCACACGGATGACT
Growth hormone            255       (F) TGAGTCGTTCAGGGTTGCACTT
  receptor (a)                      (R) CGCTGTCGCTGAATTCACCAAA
Insulin-like              156       (F) AGTGTACCATGCGCTGTCTC
  growth factor 1                   (R) AAAAGCCCCTGTCTCCACAC
Insulin-like growth       142       (F) AGTCAACGCGATACGCAAGAA
  factor-binding                    (R) TGTTTGTCGCAGTTTGGCAG
  protein 1
Insulin-like growth       152       (F) AGGCAAAGGGCTGCTGCCGGTG CGCTGG
  factor receptor                   (R) GCTCGTTGGACATGCCCTGGTA GGGCTG
  (a and b)
Thyroglobulin             182       (F) CTCCGACCATTCTCTCGCTC
                                    (R) GAGAGCAAAAGACCTGCCCT

Table 2. Correlation between components and variables.

Variable   PC1     PC2     PC3     PC4

gck        -0.67    0.35   -0.33    0.15
igf-1r     -0.49    0.22    0.11    0.74
ghrelin    -0.36   -0.44    0.75    0.04
ghr        -0.04    0.66    0.29   -0.57
cck         0.02    0.92    0.22   -0.09
igf-1       0.21     0.4    0.77    0.25
gh          0.56    0.51   -0.44    0.39
tg          0.65   -0.13    0.49    0.22
igfbp-1     0.86   -0.09   -0.14    0.05
COPYRIGHT 2017 Pontificia Universidad Catolica de Valparaiso, Escuela de Ciencias del Mar
No portion of this article can be reproduced without the express written permission from the copyright holder.
Copyright 2017 Gale, Cengage Learning. All rights reserved.

Article Details
Printer friendly Cite/link Email Feedback
Title Annotation:Research Article
Author:Opazo, Rafael; Valladares, Luis; Romero, Jaime
Publication:Latin American Journal of Aquatic Research
Date:Sep 1, 2017
Previous Article:Modelling fishery-induced impacts on the food web of the continental shelf off central-south Chile using a size-based network approach.
Next Article:Reproductive biology of Macrobrachium amazonicum (Heller, 1862) (Decapoda: Palaemonidae) in a reservoir situated in Minas Gerais State, southeastern...

Terms of use | Privacy policy | Copyright © 2019 Farlex, Inc. | Feedback | For webmasters