Research Article

 

Molecular Characterization of GHR1 Gene and Expression Analysis of Some Growth-Related Genes in Oreochromis niloticus

 

Rania Said Aboukila1, Shabaan Abd Elatif Hemeda2, Abeer Fikry El Nahas2, Walaa Slouma Hamouda Abd El Naby2*

1Veterinarian in the Veterinary Medicine Directorate- Kafer Elsheikh, Egypt; 2Genetics and Genetic Engineering-Department of Animal Husbandry and Animal Wealth Development, Faculty of Veterinary Medicine, Alexandria University, Egypt.

 

Received | March 12, 2021; Accepted | March 17, 2021; Published | June 15, 2021

*Correspondence | Walaa S.H. Abd El Naby, Genetics and Genetic Engineering-Department of Animal Husbandry and Animal Wealth Development, Faculty of Veterinary Medicine, Alexandria University, Egypt; Email: [email protected]

Citation | Aboukila RS, Hemeda SAE, El Nahas AF, Abd El Naby WSH (2021). Molecular characterization of GHR1 gene and expression analysis of some growth-related genes in oreochromis niloticus. Adv. Anim. Vet. Sci. 9(7): 1025-1033.

DOI | http://dx.doi.org/10.17582/journal.aavs/2021/9.7.1025.1033

ISSN (Online) | 2307-8316; ISSN (Print) | 2309-3331

Copyright © 2021 Abd El Naby et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

 

INTRODUCTION

 

Fish has much economic importance as a source of animal protein as well as many valuable products (Yan and Wang, 2010). The growth rate of cultured fish is one of the most important factors for the success of aquaculture. Body mass is depending on skeletal muscles because they occupy 70% of body weight. The growth of fish skeletal muscle is controlled by group of genes as growth hormone (GH), growth hormone receptor (GHR), insulin-like growth factor (IGF1), myogenic regulatory factors (MRFs), myostatine (MSTN) and others (De-Santis et al., 2007; Seiliez et al., 2012; Fuentes et al., 2013; Dias et al., 2019; Teng et al., 2020). The growth variation that related to extrinsic factors such as feeding protocol, duration of the feeding, and others, influence on the juvenile Nile tilapia body weight gain through their effect on the expression of muscle growth-related genes (Nebo et al., 2013). However, this variation caused by an extrinsic factor has a genetic base caused by the polymorphism at the coding regions of some of these genes are associated with variation in body size (Wringe et al., 2010; Tong and Sun, 2015). Selection for growth traits in farm animals and fish based on the polymorphism in the growth-related genes can be screened using many techniques (Yang et al., 2013). Genetic polymorphism is a difference in DNA sequence among individuals, groups, or populations. It includes single nucleotide polymorphisms (SNPs), sequence repeats, recombination, insertions, and deletions. SNPs reflect a particular position at which more than one nucleotide is identified and only two alleles at a SNP locus. The SNPs may lead to a change in the phenotypic appearance of the individual. The variation can be utilized as a marker for good traits and culling individuals having weak traits (Vignal et al., 2002).

 

Growth hormone receptor (GHR) is a critical regulator of growth and metabolism. It induces a phosphorylation cascade for signaling and gene expression by binding with GH (Ma et al., 2007). The GH belongs to the GH/prolactin (PRL)/ somatolactin (SL) family. It is a pluripotent hormone produced by the pituitary gland which has an important role in many physiological functions such as growth, metabolism, cellular differentiation, immunity, and reproduction in vertebrate (Rienecke et al., 2005). The IGFs are the mediators of the growth hormone that promotes mitosis, induces cell differentiation, and regulates body growth IGF-1 is structurally related to pro-insulin (Laron, 2001). Moreover, IGF-1 is a potent mitogenic factor that induces growth and stimulates differentiated functions of several cell types (Rienecke et al., 2005). The Myogenine gene (MyoG) is one of the members of the myogenic regulatory factors (MRFs) which include also MyoD, Myf5, and MRF4. The MRFs can convert a variety of non-muscle cells into myoblasts and myotubes (Johnston et al., 2011). MSTN gene is an important negative regulator gene of the skeletal muscle mass; it is known as growth and differentiation factor-8 (GDF-8), it inhibits muscle growth (Baylan et al., 2015).

 

Little data are available on such comparison of gene expression profile of growth-related genes among the small and large weight Nile tilapia, also, no data are available on the polymorphism of GHR1 among different body weight fish. So, this study was conducted to investigate GHR1 gene polymorphisms in large and small body weight using DNA sequencing. Also, to quantify the expression levels of some growth-related genes (GH, GHR1, IGF-1, MyoG, and MSTN) among the variant sizes’ mono-sex Nile tilapia.

 

MATERIAL AND METHODS

 

Experimental design and samples collection

Thirty mono-sex Nile Tilapia (Oreochromis niloticus) were obtained from a private mono-sex fish farm in El-Riad, Kafr El-sheikh, Egypt from the same bond during the fishing time. The collected fish were divided into two groups according to their body weight. The large body weight group (n=15) ranged from 394 to 485 gm, and small body weight group ranged from 105 to 290 gm. Whole blood was collected from the caudal vein on anti-coagulant and was stored at -20 ºC until DNA extraction for screening of the polymorphism of GHR1gene.

 

Other selected fish were categorized according to their weights into large (412, 427, and 475 gm), medium (229, 235, and 242 gm) and small (132, 142, and 183 gm) and were used for gene expression study. Liver tissue samples were collected and stored at -80ºC for RNA Extraction. The experimental design is illustrated in Fig. (1). All experimental procedures used for this study were approved by the local ethics committee of animal use, Faculty of Veterinary Medicine, Alexandria University, Egypt.

 

 

DNA Extraction, PCR amplification of GHR1 gene, DNA sequencing and data analysis

DNA was extracted using Gene JET Whole Blood Genomic DNA Purification Mini Kit (Thermo Scientific, Lithuania), according to manufacture instructions. The primers used for amplification of GHR1 were shown in Table (1). A segment of GHR1 from 1430-2779 bp was amplified using two primers, they were designed using NCBI tool to amplify region one (coding region; 709 bp, from 1430 bp to 2138 bp), and primer 2 to amplify region two (700 bp, from 2080 to 2779) (represent the terminal part of the coding region and 3`-untranslated region of GHR1 gene (2237-2779 bp). Amplification of GHR1 gene fragment was performed in 50 μl reaction volume containing 25 μl master mix, 5 μl genomic DNA, 1 μl of each primer and 18 μl dH2O. The final reaction mixture was placed in a thermal cycler (Techne, TC-3000, USA). The PCR program was carried out by initial denaturation

 

Table 1: Primers sequences used for amplification of the coding region of GHR1 gene.

 

Primers	Sequence (5`→3`)	Annealing temperature (ºC)	Amplicon size (bp)
Primer 1	F: TCCCAACCTGTCAGTCAACAC R: AAGGCGGGTTAAGGAGCAG	60	709
Primer 2	F: AGACTACACAGTGGTTCAGGAG R: TGAGTTAAGAAAGATAGGAACGGT	57	700

 

Table 2: List of primers used in qRT-PCR.

 

Primers	Sequence (5`→3`)	Accession No.	Amplicon size (bp)
GH	F: GGAGTTTCCCAGTCGCTCTC R: TCAGCAGCAAGATTCCCGTT	XM_003442542	95
GHR1	F: ACAGCCATCACAGTCTGCTC R: TGTACCCCACAGGCATAGGT	XM_019365174	99
IGF-1	F: CTGTGGAGAGCGAGGCTTTT R: CAGCTTTGGAAGCAGCACTC	NM_001279503	99
MyoG	F: CAGCCACACTGAGGGAGAAG R: TCTGGTTGGGGTTCATCAGC	NM_001279526	90
MSTN	F: GCATCTGTCTCAGATCGTGCT R: TGCCATCATTACAATTGTCTCCG	XM_003458832	95
β-actin	F: TGGCAATGAGAGGTTCCG R: TGCTGTTGTAGGTGGTTTCG	XM_003443127	136

GH = Growth Hormone; GHR1 = Growth Hormone receptor1; IGF-1=Insulin like growth factor1; MyoG= Myogenine; MSTN = myostatin and β-actin=beta actin

 

at 95 ºC for 5 min followed by 35 cycles of 94 ºC for 1 min for DNA denaturation, annealing temperature as shown in Table (1) for 1 min, extension at 72 ºC for 1 min and final extension at 72 ºC for 10 min according to manufacture instructions of 2X PCR Master mix Solution, i-TaqTM, Intron Biotechnology, Korea with few modifications. PCR products of each sample and 100 bp DNA ladder (GeneDirex Taiwan) were loaded in 2% agarose gel stained with ethidium bromide. The electrophoresis gel was visualized and photographed using the gel documentation system (Gel Doc. Alpha-chem. Imager, USA).

 

PCR products were purified by using MEGAquick-spintm total fragment DNA purification kit (Intron Biotechnology, Korea) according to manufacture instruction. The sequence was performed for five fish (two large and three small weights) for each primer. The purified PCR products were sequenced in one direction (forward primer) by DNA sequencer (LGC Company, Germany). The sequence results were analyzed using Chromas 1.45 (http://www.technelysium.com.au), then the sequences’ comparisons were done using the BLAST database (https://blast.ncbi.nlm.nih.gov/Blast.cgi) The alignment of obtained sequences was done by Clustalw program version 1.8 (Thompson et al., 1994). The amino acid translation was done using the ExPASy-Translate tool (https://web.expasy.org/translate)

 

Total RNA Extraction, and cDNA Synthesis

The RNA was extracted from the liver tissue of nine fish (three large, three medium, and three small body weights) by Trizol (BioFlux, Japan) according to manufacture instructions. The extracted RNA from each group was used for cDNA synthesis by Hicen Script™ cDNA Synthesis Kit (Intron Biotechnology, United Kingdom) according to manufacture instruction.

 

Quantitative real-time PCR of selected target genes and data analysis

The qRT-PCR reaction was carried out to determine GHR1, GH, IGF-1, MyoG, MSTN and β-actin expression levels using primers listed in (Table 2). The reaction mixture is 20 μl, consisting of 2 μl cDNA, 0.8 μl of each primer (50 nm), 10 μl of Syber green (Bioline, United Kingdom), and 6.4 μl of RNase-free H2O. The qRT PCR program was an initial denaturation at 95°C for 10 min followed by 40 cycles of 95°C for 15 sec and annealing temperature 60°C.

 

The target genes’ expressions of GHR1, GH, IGF-1, MyoG and MSTN were calculated using a comparative threshold cycle (Ct) method (Rao et al., 2013), and the results were reported as fold change as compared to the calibrator (medium-size fish) after normalization of the transcript amount to β-actin gene as an endogenous control. The statistical significance analysis was done using one-way analysis of variance (ANOVA). Data are expressed as means±SD and P ≤0.05 are considered significant.

 

Table 3: Comparison of the SNPs in region 1 of GHR1 among the large and small size fish.

 

Large size fish			Small size fish
Fish 1	Fish 2	Fish 3	Fish 4	Fish 5	Fish 6
-	-	-	-	1663 C>T	-
-	-	-	-	1710 C>G	-
-	-	-	-	2016 C>T	-
2161 C>T	-	-	-		2161 C>T
-	-	-	-	2200 G>A	-

 

Table 4: Comparison of the SNPs in the coding and non-coding region 2 of GHR1 among the large and small size fish.

 

Large size fish			Small size fish
Fish 1	Fish 2	Fish 3	Fish 4	Fish 5	Fish 6
-	-	2116 C>A	-	-	-
-	-	2117 A>G	-	-	-
-	2161 C>T	-	-	-	-
-	-	-	-	2123 C>G	-
-	-	-	2200 G>A	2200 G>A	2200 G>A
-	-	-	2301 G>A	-	-
-	-	-	-	2309 C>T	-
-	-	-	2388 A>T	2388 A>T	2388 A>T
-	2501 A>G	2501 A>G	2501 A>G	2501 A>G	2501 A>G
-	2543 G>A	2543 G>A	2543 G>A	2543 G>A	2543 G>A
-	2625 G>A	2625 G>A	2625 G>A	2625 G>A	2625 G>A

Gray color of the tables represents SNPs at the 3`-untrnslated region of GHR1

 

RESULT

 

GHR1 gene polymorphism

The sequences of region 1 of GHR1 of Nile tilapia in two large body weights and three small body weights were submitted to gene bank and received the accession No# MW509674, MW509675, MW509676, MW509677, and MW509678. Comparison of the aligned sequences’ region 1 of GHR1 revealed a presence of a repeated SNP (2161 C>T) which was observed in both large and small body weight fish (Table 3).

 

The sequences of the region 2 of GHR1 of Nile tilapia in two large body weights and three small body weights were submitted to gene bank and received the accession No MW509679, MW509680, MW509681, MW509682, and MW509683. Sequence alignment between the studied fish revealed the presence of common SNPs at the coding region of GHR1 (2200 G>A), and the other at the 3`-untranslated region of the gene (2388 A>T) among the small body weight fish. Repeated three SNPs were also observed in both small and large body weights’ fish at the 3`-untranslated region of the gene (2501 A>G), (2543 G>A), and (2625 G>A) (Table 4).

 

The translation of region 1 of Nile tilapia GHR1 occurred in two large body weight, and three small body weight are illustrated in Figure (2). A change in amino acid number 460 from threonine to serine (460 T >S) and amino acid number 562 from threonine to isoleucine (562 T >I) was detected in fish number 5 (one of the small body weight group).

 

Amino acid sequence alignment of the coding region 2 of GHR1 is illustrated in Fig. (3). No amino acid difference was observed among the large and small body weight fish in this region.

 

 

Expression levels of selected muscle growth related genes in large, medium, and small body weight of Nile Tilapia

The expression profiles of selected growth-related genes (GH, GHR1, IGF-1, MyoG and MSTN) are variable between the small and large body weights’ fish. Expression of GH gene was significantly up regulated (P≤0.05) in small and large body weight fish to 20.49±1.04 and 6.54 ± 0.03-fold relative to control as shown in Fig. (4a). However, a significant increase in its expression was observed in small body weight fish. The GHR1 expression profile revealed a significant increase in both large and small body weight (21.7±0.46 and 8.86±2.66 fold) (Fig. 4b). The highest significant expression was observed in fish of large body weight. Meanwhile, IGF-1 expression level decreased in both large and small body weight fish (0.84 ± 0.97 and 0.10 ± 0.1) as shown in fig. (4c).

 

MyoG gene expression showed a significant upregulation (P≤0.05) in both large and small body weight fish relative to control (Fig. 4d). More increased expression of MyoG was observed in large body weight to 16.70±1.04 fold than small body weight Nile tilapia (6.80±0.60 fold). The expression level of MSTN revealed significant (P≤0.05) increased in small body weight fish (8.05±1.06 fold) and its increase in large body weight (5.99±2.06) (Fig. 4e).

 

DISCUSSION

 

In this study, the sequencing of GHR1 gene in Nile tilapia was applied on two large body weight fish and three small body weight fish. Many SNPs were detected in large and small body weight fish groups in GHR1 gene. The SNPs were detected in both regions 1 and 2. In the region one, the SNPs are presented individually in fish. However, in the region two repeated SNPs were detected among the three small-size fish and can be used as a strong marker in GHR1 gene for culling of small size fish. On the other hand, the SNPs at 2116 C>A and 2117 A>G reported in large body weight can be used as a marker of selection for increasing body weight in Tilapia fish. Tao et al. (2011) found 38 SNPs in GHRs of Cyprinus carpio, but there were only five SNPs that can be used in molecular breeding. The presence of SNPs in GHR1 revealed changes of some amino acid sequences in region 1 and 2. Likewise, Prinzio et al. (2010) detected three nucleotides substitutions in GHR1 by Southern blot in Zebrafish, they also detected shorter GHR1 isoforms, which have deletions that render changes in the intracellular domain, involved in different signal transduction pathways. Diopere et al. (2012) found that there were eight SNPs in GHR1 in Solea solea, which were associated with amino acid substitutions.

 

We demonstrate that the expressions of growth-related genes (GH, GHR1, IGF-1, MyoG, and MSTN) are variable between the small and large body weight fish. Similarly, many authors reported variable gene expression of growth-related genes and myogenesis between small and large body weight fish and animals (Carnevali et al., 2006; Filby and Tyler, 2007; Opazo et al., 2017).

 

Atlantic, rainbow trout and Chinook salmon showed the presence of different forms of GH, and GHR (Very et al., 2005; Von Schalburg et al., 2008). The different forms of GH and GHR in wild and domestic salmon and trout can show altered expression that could indicate functional differences. Changes in the relative ratio of expression of the corresponding genes may alter function, even though the combined overall expression remains the same (Tymchuk et al., 2009). Perez-Sanchez and Le Bail, (1999) thought that, ideal growth performance achieved with low plasma GH in combination with a high concentration of hepatic GHRs and high circulating levels of IGFs. In the present study, the GH gene expression showed a significant difference between the large and the small body weight fish. Its expression was more in the small body weight group than the large-body weight group. Similarly, Opazo et al. (2017) proved that the GH mRNA level was slightly higher in smaller larvae of zebrafish than in large larvae. However, Huang et al. (2012) reported that fast-growing group of Nile Tilapia had up regulation in GH mRNA expression level than the slow-growing group. Additionally, Ma et al. (2012) found that the extra-large Cynoglossus semilaevis male fish has increased GH mRNA level than normal male. Likewise, the GH mRNA is more expressed in fast growing juvenile group than the slow and moderate growing groups (Mun et al., 2019). The GH produced mainly by the pituitary gland (Reinecke et al., 2005). In our study, the low expression of GH in liver of large fish may be compensated by expression in other organs in the body. Also, GH may also act directly on target tissues like skeletal muscle to stimulate the production of IGF-1, which sequentially can act in a paracrine/autocrine way to stimulate tissue growth (Chauvigné et al., 2003).

 

GHR1 gene expression was found throughout the intestine, muscle, liver, ovary, spleen, and kidney but is not found in the tests of adult zebrafish. The liver GHR1 expression was reported in all species of fish (Prinzio et al., 2010). In the present study, the GHR1 gene expression showed more expression in the large body weight fish than the small-group fish. Similarly, Calduch-Giner et al. (2003) observed that, in gilthead sea bream, the levels of hepatic GHR were positively correlated with the body growth rate. Moreover, Huang et al. (2012) compared between fast and slow-growing Nile Tilapia, and they found more GHR gene expression in the fast-growing group than the slow-growing group. Mun et al. (2019) also found that the expression levels of GHR mRNA in the juvenile, red-spotted grouper liver was slightly higher in fast growing group than in moderate and slow-growing groups.

 

Segers et al. (2012) approved that GHR mRNA levels in large-size egg of mouthbrooding cichlid is more than that its level in small one. While, after hatching the young originated from a small egg had high levels of GHR mRNA levels in the yolk sac more than the young originated from a large egg. So, we could suggest that the high levels of GHR mRNA accelerated the growth rate of fish. Otherwise, Kocmarek et al. (2014) showed that GHR1 in white muscle and liver of the small-size rainbow trout group showed down regulation than the large-size group. In addition, Opazo et al. (2017) observed that no difference in mRNA level of GHR1 expression between the large and small larvae of zebrafish.

 

No significant difference in IGF-1 gene expression was observed between the large and small body weight fish, and its expression is lower than the control group in our study. Peterson et al. (2004) observed similarity in IGF-1 expression level between faster and slower growing groups of USDA103 channel catfish (Ictalurus punctatus). Similarly, Opazo et al. (2017) reported that IGF-I showed similarity in mRNA levels between the large and small larvae of zebrafish. While Carnevali et al. (2006) reported that the larvae of European sea bass have higher IGF-1 mRNA levels in the large size group more than in the small-size group in response to lactic acid bacteria probiotic. Also, the large body size zebrafish and fast-growing Nile tilapia has increased IGF-1 expression than the control (Amaral and Johnston, 2012; Huang et al., 2012). Prinzio et al. (2010) proved that the liver is the main organ contributing to IGF-I circulating levels in response to GH and observed that, the IGF-1 expressed in the slow-growing group more than its expression in fast and moderate-growing groups by 2.3-fold. However, Kocmarek et al. (2014) found that IGFs were down regulated in the liver and white muscle of small fish compared to large fish in rainbow trout.

 

In this study, the MyoG expression level showed a significant difference in large and small body weight fish than the control group. The large body weight fish group had more MyoG mRNA levels than the small body weight one. Sengupta, (2014) observed the steady expression of MyoG in different developmental stages of Labeo rohita fish between 6 hours and 12 months. While Chu et al. (2014) observed that MyoG appeared in early embryonic stage of Chinese Mandarin fish. The MyoG presented the highest levels of expression in the brain compared with the kidney, spleen, liver, heart, and white muscles in Chinese Mandarin fish (Chu et al., 2014). The expression level of MyoG in brown trout showed high level in zero days and decreased in developed stages (Churova et al., 2017).

 

The expression of liver MSTN in Nile tilapia in this study showed a high level in the small body weight fish than its level in the large body weight group. Likewise, in catfish showed increased expression of MSTN in association with a lower growth rate (Weber and Bosworth, 2005). Furthermore, Carnevali et al. (2006) showed that MSTN mRNA level in the small body weight group was more than that level in the large body weight group in the larvae of European sea bass. Furthermore, Biga and Meyer, (2009) found that there was a relation between rapidly growing GH-treated giant danio, and down-regulation of MSTN mRNA level. MSTN mRNA expression in zebrafish in the embryo was less than its expression in adult fish (Xu et al., 2003). Rodgers et al. (2003) found that in larval tilapia, prolonged fasting decreased MSTN mRNA levels and somatic growth. Whereas, fasting in adult’s tilapia does not impact myostatin mRNA levels. While Gabillard et al. (2013) observed that in many fish species muscle growth rate was not unambiguously associated to the level of MSTN mRNA expression.

 

CONCLUSION

 

Many SNPs have been detected in GHR1 gene among the large and small body weight Nile tilapia. The two SNPs reported in large body weight can be used for selection to increase body weight in tilapia fish. Moreover, two repeated SNPs were detected among the small body weight fish can be used as a marker-assisted selection for culling of small size fish. The expression profiles of selected genes (GH, GHR1, IGF-1, MyoG, and MSTN) are variable between the small and large body weight fish. The expression of GH and MSTN genes increased in small body weight fish than large body weight. GHR1 and MyoG expression increased in large body weight fish, indicate their important roles for the Nile tilapia growth traits.

 

SNPs at GHR1 and increased its expression together with the MyoG gene in large size fish proposed them as candidate genes to be more evaluated for efficient growth breeding programs in Nile tilapia.

 

ACKNOWLEDGMENTS

 

We are grateful for the members and assistants of the Department of Animal Husbandry and Animal Wealth Development, Faculty of Veterinary Medicine, Alexandria University, Egypt, for providing the necessary facilities to achieve the present study.

 

CONFLICT OF INTEREST

 

There is no conflict of interest.

 

AUTHORS CONTRIBUTIONS

 

All authors contributed equally at study design, sampling, methodology, interpretation of results, and writing of the manuscript.

 

REFERENCES

 

• Amaral IPG, Johnston IA (2012). Experimental selection for body size at age modifies early life-history traits and muscle gene expression in adult zebrafish. Exp. Biol. 215: 3895–3904. https://doi.org/10.1242/jeb.068908

• Baylan M, Mazi G, Ozcan N, Devrim BO, Güzel A I (2015). Investigation of the polymorphism of the myostatin gene by using PCR-RFLP in two tilapia species (Oreochromis aureus and Oreochromis niloticus). Food Agric Environ. 13 (2): 228-231.

• Biga PR, Meyer J (2009). Growth hormone differentially regulates growth and growth-related gene expression in closely related fish species. Comp Biochem Physiol B Biochem Mol Biol. 154: 465–473. https://doi.org/10.1016/j.cbpa.2009.07.023

• Calduch-Giner JA, Mingarro M, de Celis SV Rn, Boujard D, Pérez-Sánchez J (2003). Molecular cloning and characterization of gilthead sea bream (Sparus aurata) growth hormone receptor (GHR). Assessment of alternative splicing. Comp Biochem Physiol B Biochem Mol Biol. 136 (1): 1–13.  https://doi.org/10.1016/S1096-4959(03)00150-7

• Carnevali O, Vivo L, Sulpizio R, Gioacchini G, Olivotto I, Silvi S, Cresci A (2006). Growth improvement by probiotic in European sea bass juveniles (Dicentrar-chus labrax, L.), with particular attention to IGF-1, myostatin and cortisol gene expression. Aquac. 258: 430-438. https://doi.org/10.1016/j.aquaculture.2006.04.025

• Chu W, Li Y, Wu P, Chen D, Chen J (2014). Characterization and Expression Analysis of Myogenin Gene in White Muscle of Chinese Mandarin Fish, Siniperca chuatsi. J. Proteom. Bioinform. 7: 071-076.

• Chauvigné F, Gabillard JC, Weil C, Rescan PY (2003). Effect of refeeding on IGFI, IGFII, IGF receptors, FGF2, FGF6, and myostatin mRNA in rainbow trout myotomal muscle. Gen. Comp. Endocrinol. 132: 209–215. https://doi.org/10.1016/S0016-6480(03)00081-9

• Churova MV, Meshcheryakova OV, Ruchev M, Nemova NN (2017). Age- and stage-dependent variations of muscle-specific gene expression in brown trout Salmo trutta L. Comp Biochem Physiol B Biochem Mol Biol. 211: 16-21. https://doi.org/10.1016/j.cbpb.2017.04.001

• De-Santis C, Jerry DR (2007). Candidate growth genes in finfish -where should we be looking?. Aquaculture. 272(1-4): 22-38. https://doi.org/10.1016/j.aquaculture.2007.08.036

• Dias MAD, Neto RVR, Filho JSB, Jaser SKK, Lago AA, Hilsdorf AWS (2019). Growth hormone gene polymorphism associated with grow-out performance of Oreochromis niloticus strains. Aquaculture. 503: 105-110. https://doi.org/10.1016/j.aquaculture.2018.12.069

• Diopere E, Hellemans B, Volckaert FA, Maes GE (2012). Identification and validation of single nucleotide polymorphisms in growth- and maturation-related candidate genes in sole (Solea solea L). Mar. Genomics. 9: 33–38. https://doi.org/10.1016/j.margen.2012.09.001

• Filby AL, Tyler CR (2007). Cloning and characterization of cDNAs for hormones and/or receptors of growth hormone, insulin-like growth factor-I, thyroid hormone, and corticosteroid and the gender-, tissue-, and developmental- specific expression of their mRNA transcripts in fathead minnow (Pimephales promelas). Gen. and Comp. Endocrinol. 150 (1): 151–163. https://doi.org/10.1016/j.ygcen.2006.07.014

• Fuentes EN, Valdés JA, Molina F, Björnsson BT (2013). Regulation of skeletal muscle growth in fish by the growth hormone Insulin-like growth factor system. Gen. and Comp. Endocrinol.192:136-148. https://doi.org/10.1016/j.ygcen.2013.06.009

• Gabillard JC, Biga BR, Rescan BY, Seiliez I (2013). Revisiting the paradigm of myostatin in vertebrates: Insights from fishes. Gen. and Comp. Endocrinol. 194: 45–54. https://doi.org/10.1016/j.ygcen.2013.08.012

• Huang CW, Li YH, Hu SY, Chi JR, Lin GH, Lin CC et al. (2012). Differential expression patterns of growth-related microRNAs in the skeletal muscle of Nile tilapia (Oreochromis niloticus). J. Ani. Sci. 90: 4266–4279. https://doi.org/10.2527/jas.2012-5142

• Johnston IA, Bower NI, Macqueen DJ (2011). Growth and the regulation of myotomal muscle mass in teleost fish. J. Exp. Biol. 214: 1617-1628. https://doi.org/10.1242/jeb.038620

• Kocmarek AL, Ferguson MM, Danzmann RG (2014). Differential gene expression in small and large rainbow trout derived from two seasonal spawning groups. BMC Genomics. 15: 57-75. https://doi.org/10.1186/1471-2164-15-57

• Laron Z (2001). Insulin- like growth factor 1 (IGF-1): a growth hormone. Mol. Pathol. 54(5): 311–316. https://doi.org/10.1136/mp.54.5.311

• Ma X, Liu L, Zhang Y, Lin H (2007). Two growth hormone receptors in Nile tilapia (Oreochromis niloticus): molecular characterization, tissue distribution and expression profiles in gonad during reproductive cycle. Institute of Aquatic Economic Animals and Guangdong Provincial Key Laboratory for Aquatic Economic Animals, School of Life Sciences, Sun Yat-Sen University, Guangzhou. https://doi.org/10.1016/j.cbpb.2007.01.021

• Ma Q, Liu SH, Zhuang Z, Lin L, Sun Z, Liu CH, Ma H, Su Y, Tang Q (2012). Genomic structure, polymorphism and expression analysis of the growth hormone (GH) gene in female and male Half-smooth tongue sole (Cynoglossus semilaevis). Gene, 493: 92–104. https://doi.org/10.1016/j.gene.2011.11.015

• Mun SH, You JH, Oh HJ, Lee CH, Baek HJ, Lee YD, Kwon JY (2019). Expression Patterns of growth-related genes in Juvenile Red Spotted Grouper (Epinephelus akaara) with Different Growth Performance after Size Grading. Dev. Reprod. 23 (1): 35-42. https://doi.org/10.12717/DR.2019.23.1.035

• Nebo C, Portella MC, Carani FR, Alves AFL, Padovani CR, Carvalho RF et al (2013). 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. 164 (4): 268-74. https://doi.org/10.1016/j.cbpb.2013.02.003

• Opazo R, Valladares L, Romero J (2017). Comparison of gene expression patterns of key growth genes between different rate growths in zebrafish (Danio rerio) siblings. Lat. Am. J. Aquat. Res. 45 (4): 766-775. https://doi.org/10.3856/vol45-issue4-fulltext-12

• Perez-Sanchez J, Le bail PY (1999). Growth hormone axis as marker of nutritional status and growth performance in fish. Aquac. 177 (1-4), 117–128. https://doi.org/10.1016/S0044-8486(99)00073-3

• Peterson BC, Waldbieser GC, Bilodeau L (2004). IGF-I and IGF-II mRNA expression in slow and fast-growing families of USDA103 channel catfish (Ictalurus punctatus). Comp. Biochem. Physiol. A Mol. Integr. Physiol. 139: 317–323. https://doi.org/10.1016/j.cbpb.2004.09.015

• Prinzio M C, Botta BE, Barriga EH, Rios EA, Reyes AE, Arranz SE (2010). Growth hormone receptors in zebrafish (Danio rerio): Adult and embryonic expression patterns. Gene Express. Patterns. 10 (4-5): 214-225. https://doi.org/10.1016/j.gep.2010.03.001

• Rao X, Huang X, Zhou Z, Lin X (2013). An improvement of the 2ˆ(–delta delta CT) method for quantitative real-time polymerase chain reaction data analysis. Biostat Bioinforma Biomath. 3(3): 71–85.

• Reinecke M, Björnsson BT, Walton W, Dickho V, Stephen DM et al (2005). Growth hormone and insulin-like growth factors in fish. Where we are and where to go. Gen. Comp. Endocrinol. 142: 20–24. https://doi.org/10.1016/j.ygcen.2005.01.016

• Rodgers BD, Gregory MW, Kevin MK, Michael A (2003). Prolonged fasting and cortisol reduce myostatin mRNA levels in tilapia larvae; short-term fasting elevates. Am. J. Physiol. Reg. I. 284: R1277–R1286. https://doi.org/10.1152/ajpregu.00644.2002

• Segers FH, Berishvili G, Taborsky B (2012). Egg size-dependent expression of growth hormone receptor accompanies compensatory growth in fish. Roy. Soc., B. 279: 592–600. https://doi.org/10.1098/rspb.2011.1104

• Sengupta A, Mukherjee S, Bhattacharya S, Saha S K, Chattopadhyay A (2014). Expression Pattern of Myogenic Regulatory Transcription Factor mRNAs in the Embryo and Adult Labeo rohita (Hamilton, 1822). Int. J. Zool. 9. https://doi.org/10.1155/2014/259685

• Seiliez I, Sabin N, Gabillard J (2012). Myostatin inhibits proliferation but not differentiation of trout myoblasts. Mol. Cell Endocrinol. 351 (2): 220-226. https://doi.org/10.1016/j.mce.2011.12.011

• Tao WJ, Ma LJ, Ruan RX, Yu JH (2011). SNP loci associated with weight gain on growth hormone receptor genes in Cyprinus carpio var. jian. Acta Hydrobiol Sinica. 35(4): 622–629.

• Teng T, Zhao X, Li C (2020). Cloning and expression of IGF-I, IGF-II, and GHR genes and the role of their single-nucleotide polymorphisms in the growth of pikeperch (Sander lucioperca). Aquac. Int. 28(4): 1547-1561. https://doi.org/10.1007/s10499-020-00542-z

• Thompson JD, Higgins DG, Gibson TJ, Clustal W (1994). Improving sensitivity and progressive multiple sequence alignment through sequence weighting, position specific gap penalites and weight matrix choice. Nucleic Acids Res. 22: 4673-4680. https://doi.org/10.1093/nar/22.22.4673

• Tong J, Sun X (2015). Genetic and genomic analyses for economically important traits and their applications in molecular breeding of cultured fish. Sci. China Life Sci. 58 (2):178‐186. https://doi.org/10.1007/s11427-015-4804-9

• Tymchuk WE, Beckman B, Devlin RH (2009). Altered expression of growth hormone/Insulin-Like Growth Factor I axis hormones in domesticated fish. Endocrinol. 150 (4): 1809–1816. https://doi.org/10.1210/en.2008-0797

• Very NM, Kittilson JD, Norbeck LA, Sheridan MA (2005). Isolation, characterization, and distribution of two cDNAs encoding for growth hormone receptor in rainbow trout (Oncorhynchus mykiss). Comp Biochem Physiol B Biochem Mol Biol. 140 (4): 615–628. https://doi.org/10.1016/j.cbpc.2004.12.008

• Vignal A, Milan D, Sancristobal M, Eggenb A (2002). A review on SNP and other types of molecular markers and their use in animal genetics. Genet. Sel. Evol. 34: 275-305 https://doi.org/10.1186/1297-9686-34-3-275

• Von Schalburg KR, Yazawa R, Boer J, Lubieniecki K P, Goh B, Straub, Ch. A. et al., (2008). Isolation, characterization and comparison of Atlantic and Chinook salmon growth hormone 1 and 2. BMC genomics. 9: 522-533. https://doi.org/10.1186/1471-2164-9-522

• Weber TE, Bosworth BG (2005). Effects of 28-day exposure to cold temperature or feed restriction on growth, body composition, and expression of genes related to muscle growth and metabolism in channel catfish. Aquaculture. 246: 483–492. https://doi.org/10.1016/j.aquaculture.2005.02.032

• Wringe BF, Devlin RH, Ferguson MM, Moghadam HK, Sakhrani D, Danzmann RG (2010). Growth-related quantitative trait loci in domestic and wild rainbow trout (Oncorhynchus mykiss). BMC Genet. 11: 463–476. https://doi.org/10.1186/1471-2156-11-63

• Xu C, Wu G, Zohar Y, Du S (2003). Analysis of myostatin gene structure, expression and function in zebrafish. J. Exp. Biol. 206: 4067-4079. https://doi.org/10.1242/jeb.00635

• Yan B, Wang Z (2010). Growth salinity tolerance and microsatellite analysis of the F2 reciprocal hyprids of Oreochromis niloticus, Sarotherdon galilaeus at different salinities. Aquac. Res. 41: e336 -e344. https://doi.org/10.1111/j.1365-2109.2010.02542.x

• Yang W, Kang X, Yang Q, Lin Y, Fang M (2013). Review on the development of genotyping methods for assessing farm animal diversity. Anim. Sci. Bio. 4:2. https://doi.org/10.1186/2049-1891-4-2