Novel single nucleotide polymorphisms in candidate genes for growth in tilapia ( Oreochromis niloticus )

ABSTRACT

The objective of the present work was to identify and validate single nucleotide variations located in candidate genes to growth traits in tilapia (Oreochromis niloticus). Two transitions were identified in the promoter region of the growth hormone gene (GH); eight nucleotide changes were identified in introns and promoter region of the IGF-I gene; and a transition (T/C) was identified in the Myogenin gene (MyoG). The highest genotypic frequency (0.8) for GHpA1 and MyoG was found in the GG and TT homozygous individuals, while the highest frequency (0.9) for GHpB1 was observed in the CT heterozygous fish. There was no genotypic frequency in the CC homozygous tilapia for the GHpB1 and MyoG markers. Based on their allelic frequencies, validation as novel single nucleotide polymorphisms (SNP) of those variations located at O. niloticus GH and MyoG genes was possible. These new markers will allow their association with growth traits in tilapia to be exploited in order to determine their potential use as assisted-selection markers.

Key Words:
assisted selection markers; genotypic frequency; growth hormone gene; myogenin gene

Introduction

Determination of the genetic variability is an essential step for the implementation of genetic improvement programs focused on the selection of faster growing fish with lower feed conversion rates and resistant to diseases (Lupchinski Jr. et al., 2011). Studies of genetic diversity at DNA level represent an expansion field in aquaculture aimed at finding out those DNA variations associated with productive phenotypes, in order to use them as tools for assisting the offspring selection at any early stage and predict their productive performance (Na-Nakorn and Moeikum, 2009). This strategy is known as gene-assisted selection (GAS) (De-Santis and Jerry, 2007). Application of GAS requires molecular markers. Currently, many aquatic species are improving their productive traits SNP as molecular markers from GAS (De-Santis and Jerry, 2007). Fish growth is an important productive trait and identification of SNP for growth performance located in the somatotropic axis genes has been achieved for different species (De-Santis and Jerry, 2007). There are some candidate genes associated with growth in different species of fish such as the growth hormone (GH) (Blanck et al., 2009; De-Santis and Jerry, 2007), insulin-growth factor I (IGF-1) (De-Santis and Jerry, 2007; Thompson et al., 2010), and myogenin (MyoG) (Ulloa et al., 2013; Zou and Liu, 2015) genes, which were selected for this work. The growth hormone plays an important role in the stimulation of somatic growth in teleost (Tao and Boulding, 2003), while IGF-I is a mediator of the endocrine action of hormonal growth and is able to stimulate the synthesis of DNA and subsequent cell division in many cell types (Wood et al., 2005). Myogenin belongs to the family of myogenic regulatory factors: MRF (MyoD, Myf-5, Myogenin, and MRF- 4); this gene is a specific transcription factor in muscle basic-helix-loop-helix involved in coordination of skeletal muscle development and also participates in the myogenesis and its repair (Wu et al., 2012).

Among the most exploited fish species, tilapia is considered the one with best conditions for freshwater aquaculture (Lupchinski Jr. et al., 2011). However, tilapia production depends on the improvements in the farming systems that allow the reduction of production costs increasing productivity. To date, there are no reports on molecular markers allowing the GAS application for growth traits in tilapia. Therefore, the aim of the present work was to search genetic variation and single nucleotide polymorphisms (SNP) in three of the main genes related to growth in tilapia (O. niloticus), and their validation.

Material and Methods

Fifteen females from three commercial hatcheries in Mexico (H1 = 5, H2 = 5, and H3 = 5) were used to evaluate the experimental discovery population to identify genetic variability in O. niloticus (H1: O. niloticus, Stirling; H2: O. niloticus, Chitralada, and H3: O. niloticus, Egypt). Blood samples were extracted from the caudal fin, and DNA was isolated using the Genelute Mammalia Genomic DNA commercial kit (Cat.G1N350, Sigma - Aldrich). The specific primers were designed to amplify coded regions of three genes and promoter regions for GH and IGF-I genes using the Primer Select V7.0.0 program from Lasergene Suite, DNASTAR. We used sequences available from public database to develop tilapia-specific primer (GH, GenBank accession M84774.1; IGF-I, GenBank accession AF033797; and MyoG, GenBank accession GU246717) which was re-sequenced. DNA amplification was carried out by polymerase chain reaction (PCR), using 50 ng DNA, 1.25U Taq DNA polymerase (Promega Com, Madison, WI, USA), and 0.4 mM dNTPs. Different concetrations of MgCl₂ and alignment temperatures were used for each primer-pair (Table 1). Purified PCR products were bidirectional-sequenced in an automatic sequencer ABI PRISM 3100 Genetic Analyzer (Applied Biosystems^(r)). Once nucleotide sequences were obtained, sequence alignments and SNP identification were achieved using PolyPhred software (Nickerson et al., 1997).

Table 1
Primers used for resequencing of GH, IGF-I, and MyoG genes in tilapia

• ¹
  Annealing temperature during PCR.
• ²
  MgCl₂ concentration (mM).

The progeny from the experimental discovery population was used for validation of the observed genetic variations in the GH and MyoG genes. Eighty-six juveniles of tilapia from the discovery population were used for gene validation. The population was genotyped using allelic discrimination assays designed and synthesized in Applied Biosystems^(r) (Table 2). Because of the location of the variations found in the promoter region of the IGF-I gene, it was not possible to design assays of allelic discrimination that would allow their validation. All samples were standardized using 0.25 µL DNA (25 ng), 6.25 μL Taqman PCR master mix (Applied Biosystems), and 0.3125 μL mixture of probes. Each assay was conducted on an individual basis in optical 96 well-plates in ABI PRISM^(r) 7000 equipment (Sequence Detection System), under the following conditions: one cycle at 50 °C/2 min and 95 °C/10 min, and 40 cycles of two steps at 92 °C/15 s and 60 °C/1 min. ABI PRISM 7000 program (Real-Time Sequence Detection Software) was used for the analysis of each genotype. Results were analyzed with ABI PRISM^(r) 7000 (Real-Time Sequence Detection Software). The genotype of each sample was confirmed by visual inspection. Allelic and genotypic frequencies for three markers (GHpA1 A/G, GHpB1 C/T, and MyoG C/T) were generated according to Falconer and Mackay (2001; Figure 1).

Table 2
Primers and probes used for validation of markers GHpA1 A/G, GHpB1 C/T, and MyoG C/T

Figure 1
Procedures for selected discovery population and validation for candidate growth genes in O. niloticus.

Results and Discussion

A candidate genes approach was used to search and validate genetic variation with potential to be associated with growth traits in tilapia. Some studies support this approach as an effective tool for the discovery of SNP in genes that regulate specific traits (Primmer et al., 2002; De-Santis and Jerry, 2007; Hemmer-Hansen et al., 2011).

In GH gene, a 1.563 pb region was re-sequenced (100% promoter region, 66% exons, and 42% introns), and two transitions in the promoter region were found: one located at −284 (A/G), and other at −271 (C/T). Also, the 2,944 pb of the complete IGF-I gene was re-sequenced. A total of eight nucleotide changes were identified in discovery population (six transversions and two transitions), four located in the introns (T/A, G/A, C/G, and T/G), and four in the promoter region (G/C, T/A, C/G, and C/T). The average number of polymorphisms present in introns was higher compared with exons and promoter regions. A 1,146 bp region from MyoG gene was re-sequenced and included 66% of the gene coded region. This was achieved in six individuals where a transition in the non-coding region (C/T) was identified. A total of 11 nucleotide variations were identified in the discovery population. Insulin growth factor I gene showed the highest genetic variability, followed by GH and MyoG genes. For GH and MyoG genes, observed variations were validated with minor allelic frequencies greater than 1% (Table 3).

Table 3
Genotypic and allelic frequencies for markers GHpA1 A/G, GHpB1 C/T, and MyoG C/T in tilapia

Several studies in fish species have supported the identification, validation, and association of SNP located in candidate genes for growth, maturation, and reproduction traits. The validated markers found in this study for O. niloticus will allow a better understanding of heritability and molecular basis of this complex trait (Diopere et al., 2013).

Conclusions

Novel and validated SNP obtained in this study will allow the evaluation of their association in tilapia growth traits in order to determine their potential use as assisted-selection markers.

Acknowledgments

This study was supported by the projects SIP-IPN (20120542; 20131517; 20141467). H. Rodríguez-González acknowledges COFAA and EDI for their economic support. B. L. Cuevas-Rodríguez acknowledges financial support during their Ph.D. studies to CONACYT (Reg. No. 269508) and Instituto Politécnico Nacional for PIFI fellowship. The authors also thank Dr. Carlos L. Calderón Vásquez and the staff of Genómica Funcional Laboratory, CIIDIR- SINALOA-IPN.

References

• Blanck, D. V.; Gasparino, E.; Ribeiro, R. P. and Marques, D. S. 2009. Polymorphism in the GH1-PstI gene associated to corporal characteristics in Nile tilapia strains. Pesquisa Agropecuária Brasileira 44:599-604.
• De-Santis, C. and Jerry, D. R. 2007. Candidate growth genes in finfish - where should we be looking? Aquaculture 272:22-38.
• Diopere, E.; Hellemans, B.; Volckaert, F. A. and Maes, G. E. 2013. Identification and validation of single nucleotide polymorphisms in growth-and maturation-related candidate genes in sole (Solea solea L.). Marine Genomics 9:33-38.
• Falconer, D. S. and Mackay, T. F. C. 2001. Introducción a la genética cuantitativa. 1.ed. Acribia, Barcelona, España.
• Hemmer-Hansen, J. A. K. O. B.; Nielsen, E. E.; Meldrup, D. and Mittelholzer, C. 2011. Identification of single nucleotide polymorphisms in candidate genes for growth and reproduction in a nonmodel organism; the Atlantic cod, Gadus morhua Molecular Ecology Resources 11:71-80.
• Lupchinski Jr., E.; Vargas, L.; Lopera-Barrero, N. M.; Ribeiro, R. P.; Povh, J. A.; Gasparino, E.; and Braccini, G. L. 2011. Genetic characterization of three Nile tilapia (Oreochromis niloticus) strains. Archivos de Zootecnia 60:985-995.
• Na-Nakorn, U. and Moeikum, T. 2009. Genetic diversity of domesticated stocks of striped catfish, Pangasianodon hypophthalmus (Sauvage, 1878), in Thailand: Relevance to broodstock management regimes. Aquaculture 297:70-77.
• Nickerson, D. A. O.; Vincent, L. T. and Scott, L. T. 1997. PolyPhred: automating the detection and genotyping of single nucleotide substitutions using fluorescence-based resequencing. Nucleic Acids Research 25:2745-2751.
• Primmer, C. R.; Borge, T.; Lindell, J. and Saetre, G. P. 2002. Single-nucleotide polymorphism characterization in species with limited available sequence information: high nucleotide diversity revealed in the avian genome. Molecular Ecology 11:603-612.
• Tao, W. J. and Boulding, E. G. 2003. Associations between single nucleotide polymorphisms in candidate genes and growth rate in Arctic charr (Salvelinus alpinus L.). Heredity 91:60-69.
• Thomson, M. J.; Ismail, A. M.; McCouch, S. R. and Mackill, M. J. 2010. Marker assisted breeding. p.451-469. In: Abiotic stress adaptation in plants: physiological, molecular and genomic foundation. Pareek, A.; Sopory, S. K.; Bohnert, H. J. and Govindjee, eds. Springer, New York.
• Ulloa, P. E.; Peña, A. A.; Lizama, C. D.; Araneda, C.; Iturra, P.; Neira, R. and Medrano, J. F. 2013. Growth response and expression of muscle growth-related candidate genes in adult zebrafish fed plant and fishmeal protein-based diets. Zebrafish 10:99-109.
• Wood, A. W.; Duan, C. and Bern, H. A. 2005. Insulin-like growth factor signaling in fish. International Review of Cytology 243:215-285.
• Wu, Q.; Yao, H. D.; Zhang, Z. W.; Zhang, B.; Meng, F. Y.; Xu, S. W. and Wang, X. L. 2012. Possible correlation between selenoprotein W and myogenic regulatory factors in chicken embryonic myoblasts. Biological Trace Element Research 150:166-172.
• Zou, L. and Liu, B. 2015. Identification of a Serum amyloid A gene and the association of SNPs with Vibrio-resistance and growth traits in the clam Meretrix meretrix. Fish & Shellfish Immunology 43:301-309.

Publication Dates

• Publication in this collection
  June 2016

History

• Received
  22 June 2015
• Accepted
  10 Apr 2016