Single-nucleotide polymorphism analysis of GH, GHR, and IGF-1 genes in minipigs

Abstract

Tibetan (TB) and Bama (BM) miniature pigs are two popular pig breeds that are used as experimental animals in China due to their small body size. Here, we analyzed single-nucleotide polymorphisms (SNPs) in gene fragments that are closely related to growth traits [growth hormone (GH), growth hormone receptor (GHR), and insulin-like growth factor (IGF)-1)] in these pig breeds and a large white (LW) control pig breed. On the basis of the analysis of 100 BMs, 108 TBs, and 50 LWs, the polymorphic distribution levels of GH, GHR, and IGF-1 were significantly different among these three pig breeds. According to correlation analyses between SNPs and five growth traits - body weight (BW), body length (BL), withers height (WH), chest circumference (CC), and abdomen circumference (AC) - three SNP loci in BMs and four SNP loci in TBs significantly affected growth traits. Three SNP sites in BMs and four SNP sites in TBs significantly affected growth traits. SNPs located in the GH gene fragment significantly affected BL and CC at locus 12 and BL at locus 45 in BMs, and also BW, WH, CC, and AC at locus 45 and WH and CC at locus 93 in TBs. One SNP at locus 85 in the BM GHR gene fragment significantly affected all growth traits. All indices were significantly reduced with a mixture of alleles at locus 85. These results provide more information regarding the genetic background of these minipig species and indicate useful selection markers for pig breeding programs.

minipigs; SNP; GH ; GHR ; IGF-1

Introduction

Dwarfism is caused by a genetic disease that results in a disproportionately short stature. Although the mechanisms of human Laro-type dwarfism, sex-linked dwarfism in chickens, and other dwarfism conditions have been well studied (¹), the underlying cause of dwarfism in miniature pigs remains elusive. Growth hormone (GH), growth hormone receptor (GHR), and insulin-like growth factor-1 (IGF-1) are candidate genes for evaluating pig growth traits (²). GH is a polypeptide produced by the adenohypophysis, and is a main factor that promotes somatic growth in vertebrates (³). GH is also involved in morphology, physiology, metabolism, immunology, reproduction, and behavior (⁴), and it promotes skeletal system growth and facilitates amino acid incorporation during protein synthesis (⁵). GHR is a member of the type I cytokine receptor family and is a transmembrane receptor for growth hormone, which is related to drip loss (⁶), body weight (BW) (⁷,⁸), and marbling score (⁹). GHR is expressed in many tissues and mediates the effect of GH. IGF-1 has complex biological functions and plays important roles in cell differentiation and proliferation, animal growth, and metabolism. IGF-1 function is particularly important for embryonic development (¹⁰). The porcine IGF-1 gene is also important for regulating body growth, development, and metabolism. Due to the wide and diverse geographic features of China, many minipig breeds are naturally distributed around the country. Over the last two decades, some of these minipig breeds have been adopted and used as laboratory animals, including Tibetan minipigs (TBs), Guizhou xiang minipigs, Guangxi Bama minipigs (BMs), Wuzhishan minipigs, and Banna minipigs (¹¹). The Laboratory Animal Center of Southern Medical University (China) first imported TBs from Tibet to Guangzhou for laboratory animal research in 2004. The acclimatization and experimental animalization of these minipigs have been completed (¹²). TB is a unique breed that lives in high altitude environments (¹³). These pigs, which grow slowly and have thin skin, high meat cutability, and extra-fine muscle fibers, exhibit strong adaptability and resistance to harsh environments (¹⁴).

Minipig breeding has been partially hindered by the paucity of studies regarding functional genes associated with growth and development (¹⁵). BMs and TBs are valuable, rare Chinese species that play a unique role in studying new pig breeds. We previously performed single-nucleotide polymorphism (SNP) analysis and molecular genetics research on MSTN, ESR, H-FABP, SLA, DQA, DQB, and PPARr-2 genes. In the present study, we analyzed GH, GHR, and IGF-1 SNPs by performing DNA sequencing in BMs and TBs to further explore the effects of these genes on growth traits. Large white pigs (LWs) from the United Kingdom, also known as Yorkshire pigs, were used as the control.

Material and Methods

Animals

Animal experiments were performed according to the Guidelines on Animal Care and Use established by the Southern Medical University Animal Care and Use Committee. A total of 100 blood samples from 3- to 8-month-old BMs and 108 blood samples from 3- to 8-month-old TBs were collected at the Laboratory Animal Center, Southern Medical University, Guangzhou. Fifty ear samples from LWs were randomly collected at the Dalingshan Food Company, Dongguan, China. This study was carried out in strict accordance with the recommendations in the Guide for the Care and Use of Laboratory Animals of the National Institutes of Health (USA). The animal use protocol was approved by the Institutional Animal Care and Use Committee of Southern Medical University.

Genomic DNA was extracted using Tiangen genomic DNA extraction kits (Tiangen Biotech, China). The purity and concentration of the DNA samples were determined using agarose gel electrophoresis and ultraviolet spectrophotometry. The samples were diluted to 100 ng/μL and stored at −20°C.

GH, GHR, and IGF-1 gene fragment amplification

Primers for GH (Gene ID: 347618782), GHR (Gene ID: 347618778), and IGF-1 (Gene ID: 347618789) gene fragments were designed using the Primer 5.0 software (PREMIER Biosoft, Canada). Polymerase chain reaction (PCR) was carried out in a 50-μL volume containing 50 ng template DNA, 2 μL primers (10 μM each), 25 μL of 2×MasterMix (0.05 units/μL Taq DNA polymerase, 4 mM MgCl₂, 4 mM dNTPs), and double-distilled (dd)H₂O (Aidlab Biotechnologies Co., Ltd., China). PCR reactions were performed under the following conditions: a 3-min hot start at 95°C, 35 cycles of denaturation at 94°C for 30 s, annealing (GH at 60°C, GHR at 58°C, IGF-1 at 55°C) for 30 s, extension at 72°C for 90 s, and a final extension at 72°C for 10 min. PCR products were loaded onto a 1% agarose gel and visualized using a gel imaging system (Bio-Rad, USA).

DNA sequencing

The PCR products from pig samples were sequenced by the Invitrogen Trading Company (China).

Data analysis

Polymorphism results were analyzed using DNA STAR Lasergene v7.1, Chromas, Variant Reporter, Popgene32, polymorphism information content (PIC) Calc 0.6, SPSS13.0 (ANOVA, chi-square test), and other biological software by searching SNP loci sites and identifying genotypes.

Results

SNP loci of GH, GHR, and IGF-1 gene fragments in three pig breeds

A 221-bp fragment of the GH gene, a 254-bp fragment of the GHR gene, and a 486-bp fragment of the IGF-1 gene were amplified by PCR. GH, GHR, and IGF-1 gene fragments from the three pig breeds were of the expected sizes.

GH gene fragments revealed five SNP loci (A12G, T45C, G84A, G93A, and C133T) in BMs and TBs, but three SNP loci (A12G, T45C, and G93A) in LWs. Only one SNP locus (A85C) was found in the GHR gene fragment of BMs, and no loci were found in the other two pig breeds. In IGF-1 gene fragments, three SNP loci (A100T, A285G, T300C) were observed in BMs and TBs, whereas two SNP loci (A285G and T300C) were found in LWs. On the basis of the number of SNP loci in these three gene fragments among the three pig breeds, the most diversity was observed in the BM breed, and the least diversity was observed in LWs. The differences were also closely related to their geographical distribution and breeding methods. In general, BMs came from the Guangxi Provinces of China, and they were bred from an individual family with very small numbers and without any selection. TBs have a relatively smaller geographical distribution than BMs. In contrast, LWs are bred on large commercial pig farms with specific selection. The LW samples collected for this experiment came from a limited number of farmers.

GH, GHR, and IGF-1 gene fragments and their genotype frequencies

We statistically analyzed the genotype frequencies of GH, GHR, and IGF-1 gene fragments (Tables 1- 3, K independent samples test), and the PICs are shown in Tables 1- 3. The following PICs were identified at each locus of the GH gene fragments from BMs, TBs, and LWs, respectively: 0.372, 0.360, and 0.310 at A12G; 0.373, 0.375, and 0.340 at T45C; 0.122, 0.256, and 0 at G84A; 0.029, 0.211, and 0.374 at G93A; and 0.143, 0.369, and 0 at C133T. The G84A and C133T SNP loci were not found in the LW breed. In the GH gene fragment GG was only observed in TBs at locus 84. The dominant allele was G, and no AA was present at locus 93 in BMs. Ninety-four percent GA was found in LWs at locus 93. The percentages of T at locus 133 were 42.6% in TBs and 0% in LWs. The genotype distribution among the three groups was significantly different (P<0.05, chi-square test) for four of the five SNP sites.

For the GHR gene, the SNP at A85C was found only in BMs. Both TBs and LWs contain the AA allele in this locus. The genotype distribution among the three groups was significantly different (P<0.05, chi-square test).

For the IGF-1 gene fragment, three SNPs were observed (A100T, A285G, and T300C) in the two minipig breeds, but two SNPs (A285G and T300C) were observed in LWs. AA and TT were present at 98% at locus 285 and locus 300 in LWs, respectively. The percentages of A were 70.5% in BMs and 29.6% in TBs at locus 100, whereas the percentages of A were 24.5% in BMs and 68.1% in TBs at locus 285. The genotype distribution among the three groups was significantly different (P<0.05, chi-square test).

Among the three gene fragments, the average observed and expected heterozygosities of the GH gene in the TB populations were higher than those of BMs. However, it was lower for the IGF-1 gene, and no SNP sites were observed in the GHR gene of TBs (Table 4).

Relationship between growth traits and GH, GHR, and IGF-1 gene polymorphisms

SNPs in the GH gene fragment significantly affected body length (BL) and chest circumference at locus 12 and BL at locus 45 in BMs; BW, withers height (WH), CC, and abdomen circumference (AC) at locus 45; and WH and CC at locus 93 in TBs. One SNP was observed at locus 85 in BMs. The GHR gene fragment significantly affected all of the growth trait characteristics. All indices were significantly reduced with a mixture of alleles in this locus. These results indicate that this locus has the potential to be used as a breeding selection marker for BMs. SNPs in the IGF-1 gene fragment did not affect growth trait characteristics in BMs. However, the SNPs at loci 285 and 300 significantly affected BW, BL, and CC in TBs. These three indices were significantly reduced with the mixture of alleles at these two loci. These two SNPs may be used as breeding selection markers for TBs.

Discussion

The GH gene fragment displayed the most diversity among the three gene fragments and pig breeds as demonstrated by the number of SNP loci and PIC values. Only one SNP locus was found in the GHR gene fragment of BMs. Although the IGF-1 gene fragment was selected to evaluate the genetic diversity among pig breeds, our results suggest that GH is the best candidate for this purpose.

Polymorphism levels were different among the SNP sites. The average PIC is an ideal index to assess gene fragment polymorphisms. A PIC value >0.5 indicates a highly polymorphic region, 0.25<PIC<0.5 indicates a moderately polymorphic site, and a PIC value <0.25 signifies a low polymorphic site (¹⁶). In the GH gene fragment, moderately polymorphic sites were observed at A12G and T45C in BMs; A12G, T45C, G84A, and C133T in TBs, and A12G, T45C, and G93A in LWs. The only moderately polymorphic site found in BMs was at A85C in the GHR gene fragment. The three SNP sites in the IGF-1 gene fragment were moderately polymorphic for BMs and TBs but showed low polymorphism for LWs. However, only low polymorphism was observed in TBs in two previous studies (¹⁷-¹⁹) using the same IGF-1 gene fragment. The difference in these studies may be due to the different sample sizes. LW body size is much bigger than that of TBs and BMs (²⁰), and only three moderately polymorphic sites were found in LWs. These results indicate that the values of each SNP site vary among the three pig breeds.

The different SNP sites might differentially affect pig growth traits (²¹-²⁵). Among the nine SNP sites in these three gene fragments, three SNP sites in BMs and four SNP sites in TBs significantly affected growth traits. SNPs in the GH gene fragment significantly affected BL and CC at locus 12 and BL at locus 45 in BMs, whereas they affected BW, WH, CC, and AC at locus 45 and WH and CC at locus 93 in TBs. The SNP at locus 85 in the GHR gene fragment in BMs significantly affected all of the examined growth traits. All indices were reduced significantly with the mixture of alleles at this locus. These results indicate that this locus has the potential to be used as a breeding selection marker in BMs. SNPs in the IGF-1 gene fragment did not affect growth characteristics in BMs. However, the SNPs at loci 285 and 300 significantly affected BW, BL, and CC in TBs. These three indices were significantly reduced with the mixture of alleles at these two loci. Therefore, these two SNPs may be used as breeding selection markers for TBs.

Genetic heterozygosity indicates the proportion of the group with site heterozygotes at some loci. The average heterozygosity of loci reflects the level of variation in the genetic structure, and variability is directly correlated with heterozygosity and the ability to adapt to the environment (²⁰). Among the three gene fragments, the average observed and expected heterozygosities of the TB population was higher than those of BMs for the GH gene, lower for the IGF-1 gene, and no SNP sites were observed in the GHR gene of TBs.

In summary, the GH gene fragment represents the best candidate for SNP analysis among these three gene fragments, and some SNP sites are closely related to growth traits.

Acknowledgments

This research was supported by the Guangdong Natural Science Foundation (#S2013010014720); China Postdoctoral Science Foundation (#2014M550439); the National High Technology Research and Development Program of China (863 Program: #2011AA100304); Science and Technology Planning Project (#2010U1E00811) of Guangzhou City, Guangdong Province, China; Conditions for scientific research project (#2010B060500001, #2011B060300014) of Guangdong Province, China; Key Program for International S&T Cooperation Projects of China (#2011DFA33290); Torch Program (#SQ2012GGHE000003A08), and Development Foundation of Special High-tech Zone (#SQ2011GX07E02680, #SQ2011SF11CO1543).

References

• ¹
  Godowski PJ, Leung DW, Meacham LR, Galgani JP, Hellmiss R, Keret R, et al. Characterization of the human growth hormone receptor gene and demonstration of a partial gene deletion in two patients with Laron-type dwarfism. Proc Natl Acad Sci U S A 1989; 86: 8083-8087, doi: 10.1073/pnas.86.20.8083.
  » https://doi.org/10.1073/pnas.86.20.8083
• ²
  Suzuki K, Nakagawa M, Katoh K, Kadowaki H, Shibata T, Uchida H, et al. Genetic correlation between serum insulin-like growth factor-1 concentration and performance and meat quality traits in Duroc pigs. J Anim Sci 2004; 82: 994-999.
• ³
  Figueiredo MA, Fernandes RV, Studzinski AL, Rosa CE, Corcini CD, Junior AS, et al. GH overexpression decreases spermatic parameters and reproductive success in two-years-old transgenic zebrafish males. Anim Reprod Sci 2013; 139: 162-167, doi: 10.1016/j.anireprosci.2013.03.012.
  » https://doi.org/10.1016/j.anireprosci.2013.03.012
• ⁴
  Devlin RH, Yesaki TY, Biagi CA, Donaldson EM, Swanson P, Chan W. Extraordinary salmon growth. Nature 1994; 371: 209-210, doi: 10.1038/371209a0.
  » https://doi.org/10.1038/371209a0
• ⁵
  Bockamp E, Sprengel R, Eshkind L, Lehmann T, Braun JM, Emmrich F, et al. Conditional transgenic mouse models: from the basics to genome-wide sets of knockouts and current studies of tissue regeneration. Regen Med 2008; 3: 217-235, doi: 10.2217/17460751.3.2.217.
  » https://doi.org/10.2217/17460751.3.2.217
• ⁶
  Di Stasio L, Destefanis G, Brugiapaglia A, Albera A, Rolando A. Polymorphism of the GHR gene in cattle and relationships with meat production and quality. Anim Genet 2005; 36: 138-140, doi: 10.1111/j.1365-2052.2005.01244.x.
  » https://doi.org/10.1111/j.1365-2052.2005.01244.x
• ⁷
  Curi RA, Oliveira HN, Silveira AC, Lopes CR. Effects of polymorphic microsatellites in the regulatory region of IGF1 and GHR on growth and carcass traits in beef cattle. Anim Genet 2005; 36: 58-62, doi: 10.1111/j.1365-2052.2004.01226.x.
  » https://doi.org/10.1111/j.1365-2052.2004.01226.x
• ⁸
  Sherman EL, Nkrumah JD, Murdoch BM, Li C, Wang Z, Fu A, et al. Polymorphisms and haplotypes in the bovine neuropeptide Y, growth hormone receptor, ghrelin, insulin-like growth factor 2, and uncoupling proteins 2 and 3 genes and their associations with measures of growth, performance, feed efficiency, and carcass merit in beef cattle. J Anim Sci 2008; 86: 1-16, doi: 10.2527/jas.2006-799.
  » https://doi.org/10.2527/jas.2006-799
• ⁹
  Hale CS, Herring WO, Shibuya H, Lucy MC, Lubahn DB, Keisler DH, et al. Decreased growth in angus steers with a short TG-microsatellite allele in the P1 promoter of the growth hormone receptor gene. J Anim Sci 2000; 78: 2099-2104.
• ¹⁰
  Zhang J, Yang R, Sun S, Sun L, Liu Y, Zhang Y, et al. Cloning and characterization of new transcript variants of insulin-like growth factor-I in Sika deer (Cervus elaphus). Growth Horm IGF Res 2013; 23: 120-127, doi: 10.1016/j.ghir.2013.04.003.
  » https://doi.org/10.1016/j.ghir.2013.04.003
• ¹¹
  Wang GH, Yin XM, Sun X, Zhao DM. General situation of miniature pig resources at home and abroad. Chin J Comp Med 2009; 19: 71-73.
• ¹²
  Wang YJ, Liu W, Chen C, Yan LM, Song J, Guo KY, et al. Irradiation induced injury reduces energy metabolism in small intestine of Tibet minipigs. PLoS One 2013; 8: e58970, doi: 10.1371/journal.pone.0058970.
  » https://doi.org/10.1371/journal.pone.0058970
• ¹³
  Wang YJ, Wu SJ, Guo KY, Chen C, Xie Q, Gu WW, et al. (1)(8)F-FDG uptake by spleen helps rapidly predict the dose level after total body irradiation in a Tibetan minipig model. Eur Radiol 2012; 22: 1844-1851, doi: 10.1007/s00330-012-2451-0.
  » https://doi.org/10.1007/s00330-012-2451-0
• ¹⁴
  Ryan MT, Hamill RM, O'Halloran AM, Davey GC, McBryan J, Mullen AM, et al. SNP variation in the promoter of the PRKAG3 gene and association with meat quality traits in pig. BMC Genet 2012; 13: 66, doi: 10.1186/1471-2156-13-66.
  » https://doi.org/10.1186/1471-2156-13-66
• ¹⁵
  Silio L, Rodriguez MC, Fernandez A, Barragan C, Benitez R, Ovilo C, et al. Measuring inbreeding and inbreeding depression on pig growth from pedigree or SNP-derived metrics. J Anim Breed Genet 2013; 130: 349-360.
• ¹⁶
  Hao LL, Yu H, Zhang Y, Sun SC, Liu SC, Zeng YZ, et al. Single nucleotide polymorphism analysis of exons 3 and 4 of IGF-1 gene in pigs. Genet Mol Res 2011; 10: 1689-1695, doi: 10.4238/vol10-3gmr1328.
  » https://doi.org/10.4238/vol10-3gmr1328
• ¹⁷
  Lang S, Deng JZ, Hao LL, Zeng YZ, Wang J, Li MQ, et al. Single nucleotide polymorphism analysis of IGF-1 gene in some miniature pigs and large breeds. Chin J Vet Sci 2011; 31: 750-758.
• ¹⁸
  Zambonelli P, Davoli R, Bigi M, Braglia S, De Paolis LF, Buttazzoni L, et al. SNPs detection in DHPS-WDR83 overlapping genes mapping on porcine chromosome 2 in a QTL region for meat pH. BMC Genet 2013; 14: 99, doi: 10.1186/1471-2156-14-99.
  » https://doi.org/10.1186/1471-2156-14-99
• ¹⁹
  Zhang GW, Gao L, Chen SY, Zhao XB, Tian YF, Wang X, et al. Single nucleotide polymorphisms in the FTO gene and their association with growth and meat quality traits in rabbits. Gene 2013; 527: 553-557, doi: 10.1016/j.gene.2013.06.024.
  » https://doi.org/10.1016/j.gene.2013.06.024
• ²⁰
  Xiao W, Yun P, Li GB, Hang GY, Lu YQ, Wang CD. A comparative study of different sources of growth curve of large white pigs. J Anhui Agri Sci 2007; 35: 2912-2913.
• ²¹
  An X, Bai L, Hou J, Zhao H, Peng J, Song Y, et al. Molecular cloning, tissue expression and SNP analysis in the goat nerve growth factor gene. Mol Biol Rep 2013; 40: 857-863, doi: 10.1007/s11033-012-2126-6.
  » https://doi.org/10.1007/s11033-012-2126-6
• ²²
  Aslam ML, Bastiaansen JW, Elferink MG, Megens HJ, Crooijmans RP, Blomberg lA, et al. Whole genome SNP discovery and analysis of genetic diversity in Turkey (Meleagris gallopavo). BMC Genomics 2012; 13: 391, doi: 10.1186/1471-2164-13-391.
  » https://doi.org/10.1186/1471-2164-13-391
• ²³
  Badke YM, Bates RO, Ernst CW, Schwab C, Fix J, Van Tassell CP, et al. Methods of tagSNP selection and other variables affecting imputation accuracy in swine. BMC Genet 2013; 14: 8, doi: 10.1186/1471-2156-14-8.
  » https://doi.org/10.1186/1471-2156-14-8
• ²⁴
  Wang J, Wang H, Jiang J, Kang H, Feng X, Zhang Q, et al. Identification of genome-wide copy number variations among diverse pig breeds using SNP genotyping arrays. PLoS One 2013; 8: e68683, doi: 10.1371/journal.pone.0068683.
  » https://doi.org/10.1371/journal.pone.0068683
• ²⁵
  Wang L, Liu X, Zhang L, Yan H, Luo W, Liang J, et al. Genome-wide copy number variations inferred from SNP genotyping arrays using a Large White and Minzhu intercross population. PLoS One 2013; 8: e74879, doi: 10.1371/journal.pone.0074879.
  » https://doi.org/10.1371/journal.pone.0074879

• First published online August 1, 2014.

Publication Dates

• Publication in this collection
  01 Aug 2014
• Date of issue
  Sept 2014

History

• Received
  25 Feb 2014
• Accepted
  12 May 2014