Acessibilidade / Reportar erro

Association of IGF-I gene polymorphisms with milk yield and body size in Chinese dairy goats

Abstract

The association of IGF-I gene polymorphisms with certain traits in 708 individuals of two Chinese dairy-goat breeds (Guanzhong and Xinong Saanen) was investigated. Polymerase chain reaction-single strand conformation polymorphism (PCR-SSCP) and DNA sequencing methods were employed in screening for genetic variation. Two novel mutations were detected in the 5'-flanking region and in intron 4 of IGF-I gene, viz., g.1617 G > A and g.5752 G > C (accession D26119.2), respectively. The associations of the g.1617 G > A mutation with milk yield and the body size were not significant (p > 0.05). However, in the case of g.5752 G > C, Xinong Saanen dairy goats with the CG genotype presented longer bodies (p < 0.05). Chest circumference (p < 0.05) was larger in Guanzhong goats with the GG genotype. In Xinong Saanen dairy goats with the CC genotype, milk yields were significantly higher during the first and second lactations (p < 0.05). Hence, the g.5752 G > C mutation could facilitate association analysis and serve as a genetic marker for Chinese dairy-goat breeding and genetics.

IGF-I gene; dairy goat; milk yield; PCR- SSCP


ANIMAL GENETICS

SHORT COMMUNICATION

Association of IGF-I gene polymorphisms with milk yield and body size in Chinese dairy goats

Chanjuan DengI; Rongnuan MaI; Xiangpeng YueI; Xianyong LanI; Hong ChenI,II; Chuzhao LeiI

IShaanxi Key Laboratory of Molecular Biology for Agriculture, College of Animal Science and Technology, Northwest A&F University, Yangling, Shaanxi, China

IIInstitute of Cellular and Molecular Biology, Xuzhou Normal University, Xuzhou, Jiangsu, China

Send correspondence to Send correspondence to: Chuzhao Lei Shaanxi Key Laboratory of Molecular Biology for Agriculture College of Animal Science and Technology, Northwest A & F University Yangling, 712100 Shaanxi, China E-mail: leichuzhao1118@126.com

ABSTRACT

The association of IGF-I gene polymorphisms with certain traits in 708 individuals of two Chinese dairy-goat breeds (Guanzhong and Xinong Saanen) was investigated. Polymerase chain reaction-single strand conformation polymorphism (PCR-SSCP) and DNA sequencing methods were employed in screening for genetic variation. Two novel mutations were detected in the 5'-flanking region and in intron 4 of IGF-I gene, viz., g.1617 G > A and g.5752 G > C (accession D26119.2), respectively. The associations of the g.1617 G > A mutation with milk yield and the body size were not significant (p > 0.05). However, in the case of g.5752 G > C, Xinong Saanen dairy goats with the CG genotype presented longer bodies (p < 0.05). Chest circumference (p < 0.05) was larger in Guanzhong goats with the GG genotype. In Xinong Saanen dairy goats with the CC genotype, milk yields were significantly higher during the first and second lactations (p < 0.05). Hence, the g.5752 G > C mutation could facilitate association analysis and serve as a genetic marker for Chinese dairy-goat breeding and genetics.

Key words:IGF-I gene, dairy goat, milk yield, PCR- SSCP.

Animal growth is controlled by a complex system, in which the somatotropic axis plays a key role. Genes that operate in this axis are responsible for postnatal growth, mainly by way of growth hormones acting on bones and muscles mediated by insulin-like growth factor I (IGF-I) (Sellier, 2000). Thus, IGF-I is an important growth factor involved in a variety of physiological processes including cell differentiation, embryogenesis, and regulation of metabolism, reproduction, fetal development, and growth (Adam et al., 2000; Shen et al., 2003). The IGF-I gene was mapped on chromosome 5 in bovine (Miller et al., 1991). The IGF-I gene is expressed in most tissues, but mainly in liver (Wang et al., 2003). It encodes a polypeptide of molecular weight of 7.5 kDa with 70 amino acids (Daughaday and Rotwein, 1989). The amino acid sequence of IGF-I is identical in human, cattle, dog, and pig (Nixon et al., 1999).

In the nineties, several lines of evidence revealed an association of the high level of IGF-I in the blood with fast growth in cattle (Barash et al., 1998; Sirotkin et al., 2000). In the following decades, cattle IGF-I gene mutations (C/T mutation in the 5'-flanking region) were found to be associated with the percentage of fat and protein in milk (Ge et al., 2001; Eulalia et al., 2006). However, these results were quite different from Hines' discovery whence no association was found between C/T mutation in 5'-flanking region and dairy production traits in Holstein cattle (Hines et al., 1998). The association of IGF-I gene polymorphism (short random repeat in the 5'-flanking region) with body weight, both at birth and weaning, has also been described (Li et al., 2004; Moody et al., 1996), although such an association was not found by Curi et al. (2005). Recently, an IGF-I intron 2 polymorphism was found to correlate positively with the twinning rate (Kim et al., 2009).

Dairy goats are numerous in northwest China, among which, the Xinong Saanen (SN) and Guanzhong (GZ) breeds are the best-known, due not only to their tolerance of crushed feed and the local harsh weather, but also to their high milk yield. As it is difficult to culture high quality domestic breeds over a short period by the traditional breeding and genetics methods, the major breeders recently focused on using DNA markers for developing breeds through marker-assisted selection (MAS). Thus, it is important to identify significant associations between polymorphisms within crucial candidate genes and production traits. The aim was to find IGF-I gene SNPs in two Chinese dairy goat breeds and search for the possible associations of these with milk yield and body size.

708 blood samples were collected from Shaanxi province (440 for GZ breed from Sanyuan Guanzhong dairy goat breeding base in Shaanxi province; 202 for SN breed from Qianyang Saanen dairy goat breeding base in Shaanxi province; 66 for SN breed, also from animal farm in Northwest A&F University in Shaanxi province). All goats were 2 to 3 years old adults. Measurements were taken for all the individuals (BH = body height; BL = body length; ChC = chest circumference). Based on a 305-day lactation period, the first and second lactation data from each of 202 SN dairy goats from the Qianyang Saanen dairy goat breeding base, were recorded for every goat.

Blood samples were drawn from the jugular vein into curettes containing ACD (anticoagulant citrate dextrose). The samples were stored at -80 °C for several weeks, whereupon genomic DNA was extracted with the use of phenol/chloroform (Sambrock et al., 2001).

Based on the goat IGF-I gene sequence published in the GenBank (D26119.2), five pairs of primers were designed. Primer sequences, annealing temperatures and amplicons are listed in Table 1. The P5 primers for the 5' -flanking end in Table 1 referred to those used in cattle (Ge et al., 2001), whereas the remainder were designed in this study.

The final volume of 25 µL contained 50 ng of genomic DNA, 0.5 µL of each primer (10 pM), 2.5 µL of a 10xbuffer (including 20 mM Mg2+), 0.5 µL of dNTPs (10 mM) and 0.25 µL of Taq DNA polymerase (5 U/µL Dongsheng. Guangzhou, China). Amplification was carried out using PCR, under the following conditions: 5 min at 96 °C for initial DNA denaturation, 35 cycles of denaturation at 94 °C for 40 s; annealing at X °C (shown in Table 1) for 40 s for each pair of primer; extension at 72 °C for 40 s. The final cycle was followed by a 10 min extension at 72 °C.

Aliquots of 4 µL of PCR products were mixed with 8 µL of denaturing solution (95% formamide, 25 mM EDTA, 0.025% xylene-cyanole and 0.025% bromophenol blue) and heated for 10 min at 98 °C and chilled on ice. Denatured DNA was subjected to 10% PAGE (80 mm x 73 mm x 0.75 mm) in 1x TBE buffer and constant voltage 250V for 15 min, then 160V for 2 h. The gel was stained with 0.1% silver nitrate (Lan et al., 2007). After the polymorphisms were detected, the PCR products of different electrophoresis patterns were sequenced by the Huada Genetic Sequence Company (Beijing, China). The sequences were analyzed with DNASTAR software (Version 7.10, DNASTAR, Madison, USA).

Based on the individual number of different genotypes in analyzed breeds, Hardy-Weinberg equilibriums, genotypic frequencies and allelic frequencies were directly calculated. The statistical software SPSS (Version 16.0) was used to analyze the relationships between the genotypes and traits in goats. The adjusted Linear Model I with fixed effects was used to evaluate the relationships between genotypes and milk yield of Xinong Saanen dairy goats. The model applied was: Yijk = µ + Ai + Gj + (AG)ij + eijk, where Yijk was the trait measured in each ijkth animal, m the overall population mean, Ai the fixed effect due to the ith age, Gj the fixed effect associated with the jth genotype, (AG)ij the interaction between the ith age and the jth genotype, and eijk the random error. Effects associated with farm, sex and season of birth (spring versus fall), and the ages of dams and sires, were not considered in the linear model, since preliminary statistical analysis indicated that these effects did not have significant influence on the variability of traits in SN female dairy goats. The adjusted Linear Model II with fixed effects was used to analyze the relationships between genotypes and body size in these two breeds. The model applied was: Yijk = µ + Bi + Gj + (BG)ij + eijk, where Yijk was the trait measured in each ijkth animal, m the overall population mean, Bi the fixed effect due to the ith breed, Gj the fixed effect associated with jth genotype, (BG)ij the interaction between the ith breed and the jth genotype, and finally eijk the random error. Effects associated with farm, sex and season of birth (spring versus fall), as well as the ages of both dams and sires, were not used in the linear model, as the preliminary statistical analyses indicated that these effects did not have significant influences on variability of traits in this two breeds. All the goats were adults, so the age did not influence on the body size. The least square means (LSMs), with standard errors and multiple range tests for IGF-I genotypes and traits, were calculated.

Using PCR-SSCP method with five pairs of primers of the goat IGF-I gene, only two (P1 and P4) of them demonstrated polymorphic patterns in both breeds. In the P1 locus (fragment amplified by the P1 pair of primers), there were only two kinds of SSCP polymorphic patterns detected in all the 708 samples. We named them GG, with three bands, and AG, with five bands (Figure 1a). In the P4 locus (fragment amplified by the P4 pair of primers), three SSCP polymorphic patterns were found in the 708 samples, which were named GG pattern, with four bands, CG pattern, with five bands, and CC pattern also with four bands (Figure 1b).



In order to clarify the nucleotide mutations in these two loci, nine PCR products, each with a different pattern, were sequenced at each locus. The three remaining products (P2, P3, and P5) were also sequenced. On comparing with GenBank sequence D26119.2, 249 bp PCR products, including the 5'-flanking region upstream from the start codon (ATG), were amplified by the PI pair of primers. There was a G to A mutation, namely g.1617 G > A (accession D26119.2). The P4 primers were used to amplify the whole of exon 4 and part of introns 3 and 4. There was a G to C variation in both breeds, viz., g.5752 G > C (accession D26119.2), located in intron 4. No mutation was detected in the remaining loci (P2, P3, and P5).

The allelic and genotypic frequencies of each locus obtained for the different genetic groups, can be seen in Table 2. For the g.1617 G > A mutation (accession D26119.2), frequencies of the IGF-I-A allele in the analyzed populations were 0.08 and 0.11 for GZ and SN dairy goats, respectively, with genotype distribution of both breeds in agreement with Hardy-Weinberg equilibrium (p > 0.05). For the g.5752 G > C mutation (accession D26119.2), the GZ breed was not in Hardy-Weinberg equilibrium (p < 0.01), whereas the SN breed was in Hardy-Weinberg equilibrium (p > 0.05). The allelic frequencies for C were 0.19 and 0.30 for GZ and SN dairy goats, respectively (Table 2).

Combined data related to body-size and milk-yield at the first and second lactations are shown in Table 3. For g.1617 G > A (accession D26119.2), there was no significant association with body size and the milk yield in the two breeds (p > 0.05). For g.5752 G > C (accession D26119.2), significant statistical results were encountered in milk yield and BL (body length) of individuals with the CC and CG genotypes (p < 0.05 and p < 0.05, respectively) in Saanen dairy goats. The BL of SN breed with CG genotype was significantly longer than in the others. Individuals with genotype CC of the IGF-I gene had a superior milk yield at the first and second location when compared to those of individuals with genotype CG and GG (p < 0.05). In the Guanzhong population, chest circumference (ChC) of GZ goats with the GG genotype was the largest of the three genotypes, whereas ChC for individuals with CG genotype was larger than that for CC carriers (p < 0.05).

When focusing on the IGF-I gene, an attempt was made to detect nucleotide mutations in primer amplified fragments, and to analyze their associations with body size and milk yield in Chinese dairy goats. The analyzed length of the IGF-I gene was 1588 bp, including part of the 5'- flanking, and almost the entire coding region except for exon 2 which is only 15 bp long. In the g.1617 G > A mutation (accession D26119.2), the lack of AA homozygotes indicated that the mutation in the goat IGF-I gene might give rise to diminished production in the population under analysis. These two dairy goat breeds are major milk producers, originating from a varied dairy-goat breeding basis in Shaanxi province (P.R. China). As the two populations were over one year old, individuals with low production had possibly already been eliminated. Hence, it is reasonable to surmise that AA homozygotes is almost completely lost in quite large local populations. Obviously, the reason of the lack of genotype AA in these breeding farms requires further investigation.

Significant statistical results were found as regards body size and milk yield among different genotypes in g.5752 G > C (accession D26119.2) (p < 0.05) (Table 3). Being located in one of the intron regions, the SNP studied here may not be a causal mutation. Thus, we suggest that it could be linked to another mutation in the coding or regulatory regions of a gene which is a causal mutation for production traits. Introns have been shown to affect the transcriptional efficiency of numerous genes in a variety of organisms (Greenwood and Kelsoe, 2003; LeHir et al., 2003). Further analysis is required to validate both the association found in this study and the physiological significance of the intron 4 mutation in IGF-I gene.

In conclusion, our study revealed the first known significant associations of IGF-I gene polymorphisms with milk yield and body size in dairy goats. The CC genotype at the g.5752 G > C (accession D26119.2) locus could be used as a molecular marker for superior milk yield. Furthermore, this study will be of practical use for the improvement of Chinese native dairy goats and the local breeding of those with higher milk production.

Acknowledgments

This study was supported by the National 863 Program of China (No. 2008AA10Z138), the Young Topnotch Researcher Support Project of Northwest A&F University (No.QNGG-2009-007), the "13115" Sci-Tech Innovation Program of Shaanxi Province (No. 2008ZDKG-11), the Research Fund for the Doctor Program of Higher Education of China (No.20080712001), and the Natural Science Foundation of Jiangsu Province (No. BK2008120).

Received: August 19, 2009; Accepted: December 9, 2009.

Associate Editor: Bertram Brenig

  • Adam CL, Gadd TS, Findlay PA and Wathes DC (2000) IGF-I stimulation of luteinizing hormone secretion, IGF-binding proteins (IGFBPs) and expression of mRNAs for IGFs, IGF receptors and IGFBPs in the bovine pituitary gland. J Endocrinol 166:247-254.
  • Barash H, Aharoni Y, Brosh A and Holzer Z (1998) Effect of low energy diets followed by a compensatory diet on body weight gain and plasma hormone concentrations in bull calves. J Dairy Sci 81:250-254.
  • Curi RA, Oliveira HN, Silveira AC and Lopes CR (2005) Effects of polymorphic microsatellites in the regulatory region of IGF-I and GHR on growth and carcass traits in beef cattle. Livest Prod Sci 94:159-167.
  • Daughaday WH and Rotwein P (1989) Insulin-like growth factors I and II. Peptide messenger ribonucleic acid and gene structures, serum and tissue concentrations. Endocr Rev 10:68-91.
  • Eulalia S, Lech, Z, Jolanta O, Nina S, Emilia B and Krzyzewski J (2006) Effect of polymorphism in IGF-I gene on production traits in Polish Holstein-Friesian cattle. Anim Sci Pap Rep 24:225-237.
  • Ge W, Davis ME, Hines HC, Irvin KM and Simmen RC (2001) Association of a genetic marker with blood serum insulin-like growth factor-I concentration and growth traits in Angus cattle. J Anim Sci 79:1757-1762.
  • Greenwood TA and Kelsoe JR (2003) Promoter and intronic variants affect the transcriptional regulation of the human dopamine transporter gene. Genome 82:511-520.
  • Hines HC, Ge W, Zhao Q and Davis ME (1998) Association of genetic markers in growth hormone and insulin-like growth factor I loci with lactation traits in Holsteins. Anim Genet 29:69.
  • Kim ES, Shi X, Cobanoglu O, Weigel K, Berger PJ and Kirkpatrick BW (2009) Refined mapping of twinning rate quantitative trait loci on bovine chromosome 5 and analysis of insulin-like growth factor-I as a positional candidate gene. J Anim Sci 87:835-843.
  • Lan XY, Pan CY, Chen H, Zhang CL, Li JY and Zhao M (2007) An AluI PCR-RFLP detecting a silent allele at the goat POUIFI locus and its association with production traits. Small Ruminant Res 73:8-12.
  • LeHir H, Nott A and Moore M (2003) How introns influence and enhance eukaryotic gene expression. Trends Biochem Sci 28:215-220.
  • Li C, Basarab J, Snelling WM, Benkel B, Murdoch B, Hansen C and Moore SS (2004) Assessment of positional candidate genes myf5 and IGF-I for growth on bovine chromosome 5 in commercial lines of Bos taurus Anim Sci 82:1-7.
  • Miller JR, Thomsen PD, Dixon SC, Tucker EM, Konfortov BA and Harbitz I (1991) Synteny mapping of the bovine IGHG2, CRC and IGF-I genes. Anim Genet 23:51-58.
  • Moody DE, Pomp D, Newman S and Macheil MD (1996) Characterization of DNA polymorphisms in three populations of Hereford cattle and their associations with growth and maternal EPD in line 1 Herefords. J Anim Sci 74:1784-1793.
  • Nixon AJ, Brower-Toland BD and Sanell LJ (1999) Primary nucleotide structure of predominant and alternate splice forms of equine insulin-like growth factor I and their gene expression patterns in tissues. Am J Vet Res 60:1234-1241.
  • Sambrock J, Fritsch EF and Maniatis T (2001) Molecular Cloning: A Laboratory Manual. 3rd edition. Cold Spring Harbor, New York.
  • Sellier P (2000) Genetically caused retarded growth in animals. Domest Anim Endocrinol 19:105-119.
  • Shen WH, Wisniowski P, Ahmed L, Boyle DW, Denne SC and Liechty EA (2003) Protein anabolic effects of insulin and IGF-I in the ovine fetus. Am J Physiol Endocrinol Metab 284:E748-E756.
  • Sirotkin AV, Chrenek P, Makarevich AV, Huba J and Bulla J (2000) Interrelationships between breed, growth hormone, plasma IGF-I level and meat performance in bulls of different ages. ArchTierzucht 6:591-596.
  • Wang Y, Price SE and Jiang H (2003) Cloning and characterization of the bovine class 1 and class 2 insulin-like growth factor-I mRNAs. Domest Anim Endocrinol 25:315-328.
  • Send correspondence to:

    Chuzhao Lei
    Shaanxi Key Laboratory of Molecular Biology for Agriculture
    College of Animal Science and Technology, Northwest A & F University
    Yangling, 712100 Shaanxi, China
    E-mail:
  • Publication Dates

    • Publication in this collection
      30 Apr 2010
    • Date of issue
      2010

    History

    • Received
      19 Aug 2009
    • Accepted
      09 Dec 2009
    Sociedade Brasileira de Genética Rua Cap. Adelmio Norberto da Silva, 736, 14025-670 Ribeirão Preto SP Brazil, Tel.: (55 16) 3911-4130 / Fax.: (55 16) 3621-3552 - Ribeirão Preto - SP - Brazil
    E-mail: editor@gmb.org.br