Polymorphisms of the Ovine Leptin Gene and its Association with Growth and Carcass Traits in Three Iranian Sheep Breeds

Document Type: Research Paper

Authors

1 Department of Animal and Poultry Science, College of Aboureyhan, University of Tehran, P.O. Box 11365-4117, Tehran, I.R. Iran

2 Department of Clinical Genetic, National Institute of Genetic Engineering and Biotechnology, P.O. Box 14965/161,Tehran, I.R. Iran

Abstract

Leptin (LEP), the expression product of the obese gene produced primarily in the adipose tissue, is related to feed intake, growth and lipid metabolism. The aim of this study was to study the possible association between polymorphism of the LEP gene with growth and carcass traits among three Iranian sheep breeds of the Shal, Zandi and Zel varieties. A total of 180 purebred animals of the Shal, Zandi and Zel, breeds were chosen for this study. Three flocks, each comprising of 60 ewes of three breeds, were derived from Aboureyhan sheep populations. The Shal (n=18), Zandi (n=24) and Zel (n=17) lambs were screened for polymorphism of the LEP gene. Following genomic DNA extraction from whole blood samples, polymerase chain reaction (PCR) was carried out in order to amplify a 260 bp fragment of the target gene. Polymorphisms were detected using the single strand conformational polymorphism (SSCP) technique. Two genotypes of AA and AG with frequencies of 0.53 and 0.47 in the Shal 0.70 and 0.30 in the Zandi and 0.65, 0.35 in the Zel breed were observed, respectively. The frequencies of alleles A and G were 0.74, 0.26 in the Shal breed, 0.85, 0.15 in the Zandi breed and 0.82, 0.18 in the Zel breed, respectively. Chi-Square test (c2) confirmed the Hardy-Weinberg equilibrium for the LEP loci. Average heterozygosity (30%) of the LEP locus for the three breeds was slightly low. Comparison of the sequence of the target gene available in the GeneBank with the results of the present study showed two single nucleotide polymorphisms (SNP) A®G and T®C transitions at 113 and 165 bp positions, respectively. In the Shal breed, the A113G SNP associated with an increase in cold carcass weight (p< 0.01), fat-tail percent (p< 0.05) and total body fat weight (p< 0.05). In the Zel breed, the A113G SNP was associated with an increase in fat-tail percent (p< 0.05) and reduction in slaughter weight (p< 0.05), cold carcass weight (p< 0.01) and lean meat weight (p< 0.01). Therefore, a significant association between SNP within the LEP gene and certain carcass traits in the Shal and Zel breeds is proposed. In the Zandi breed the A113G SNP was not associated with carcass traits but showed a reduction in weaning weight (P< 0.05).

Keywords


INTRODUCTION

The recently discovered hormone LEP is the main factor in regulation of appetite (Thomas et al., 2001). The LEP gene was discovered in 1994 by positional cloning techniques. The gene consists of two introns and three exons that are involved in the synthesis of a protein with 167 amino acids, which is secreted into blood after cleavage of the 21-amino-acid signal peptide (Dela et al., 1996; Zhang et al., 1994). The adipose tissue is the main source of LEP. In fact, LEP’s receptors are found in many organs including the brain, adipose tissue, stomach, gonads, muscles, hypothalamus and hypophysis (Houseknecht et al., 1998; Dyer et al., 1997). LEP acts as an adipostat or circulating signal that carries information from fat depots or energy reserves to regulatory appetite centers and is thought to play an important role in many physiological functions; lipogenesis, angiogenesis, thermogenesis, feed intake regulation, fertility, fetal growth, immunological processes, body weight, hypothalamus-hypophysis-adrenal and hypophysis-gonads axis (Barb et al., 2001; Houseknecht et al., 1998; Schneider et al., 2000). Genetic differences in the LEP gene were first reported in mice and humans (Ohshiro., 2000; Halaas et al., 1995). In recent years studies have been performed on the association between  LEP gene polymorphisms and production traits in dairy cattle (Fiona et al., 2002), sheep (Boucher et al., 2006) and poultry (Taouis et al., 1998). The Shal and Zandi breeds are known as fat-tailed breeds. These two breeds differ in body size and breeding conditions. The Shal is large in size while the Zandi is medium in size. The only tailed sheep breed of Iran is the Zel. The aim of this study was to identify polymorphisms of the ovine LEP gene and its possible association with growth and carcass traits in three Iranian sheep breeds.


MATERIALS AND METHODS

Animals: A total of 180 animals belonging to the Shal, Zandi and Zel breeds were selected for this study. Three flocks, each consisting of 60 ewes of the three breeds, were derived randomly. In each breed, ewes were mated with rams of the same breed. Population sets of randomly selected purebred Shal (n=18), Zandi (n=24) and Zel (n=17) lambs were selected. The choice of these three breeds was based on their opposite phenotypes for fat deposition and muscle accretion. The Shal and Zandi are known as fat-tailed native breeds of the Tehran province. These two breeds differed in size and living conditions. The Shal is large in size while the Zandi is medium in size. The only tailed sheep breed of Iran is the Zel, living as a native of two Northern provinces of the country. This breed is small in size and early-maturing (Kashan et al., 2005). Lambs were reared with their dams until an average age of 90 days. At the time of weaning (approximately 90 days old), only the male lambs were selected in each breed. All lambs were then raised under the same controlled-environmental conditions. Animals were fed a ration without limitation to finishing weight. All lambs were slaughtered according to their targeted fattening period (90 days) after a 12 h fasting period.

Growth and Carcass data: From birth to weaning and during the 90-days fattening period, the lambs were weighed once a week. At the end of the fattening period all lambs were weighed. The weights of carcasses were also recorded. The carcasses were then chilled at 4ºC for 24 h and weighed again. Subsequently, the carcasses were split longitudinally into two parts. The left side of the carcass was dissected into lean meat, bone, subcutaneous fat (SCF), intramuscular fat (IMF), trimmings and weighed separately according to Kashan et al. (2005). Internal fat was obtained from the fat surrounding the intestine and kidney. The whole soft tissue (lean meat and fat) of the left side carcasses was grounded and passed twice through a plate with a 4mm orifice. Dry matter (DM), crude protein (CP), ash and chemical fat (ether extract) contents were detected according to the Association of Official Analytical Chemist (AOAC) standardized methods (1995). 
       Blood samples and PCR amplification: Blood samples were obtained from three pure breeds of Shal, Zandi and Zel lambs. Each sample contains 3-5 ml blood taken from jugular vein. Genomic DNA was purified from 100 ml of blood samples using the Diatom DNA prep100 (DNA extraction kit) (gene Fanavaran Co., Tehran) The Primers were designed based on available ovine genomic sequences (GeneBank accession numbers: U84247 and AY831682). The primer sequences are as follows: forward 5´-CGCAAGGTCCAGGATGACACC-3´; and reverse 5´-GTCTGGGAGGGAGGAGAGTGA-3´. The amplicons were 260 bp in length that covered exon 2 and part of intron 2 of the ovine LEP gene. The amplification reactions were performed in a total volume of 25 ml containing 200 ng of genomic DNA, 1 ml of each primer (10 pM), 0.5 ml of dNTPs (10 mM), 0.8 ml of MgCl2 (1.6 mM), 1.2 U of Taq DNA polymerase and 2.5 ml of 10X PCR buffer in 35 Cycles (initial denaturation: 94°C, 2 min;  denaturation: 94°C, 30s:  annealing:  62.5°C, 25s; extention: 72°C, 45s)  followed by a final extension step at 72°C for 7 min.

SSCP and sequencing: PCR products were purified using DNA purification kit (Gene Fanavaran Co., Tehran) before being sequenced. For genotyping of the LEP locus, SSCP was carried out on the PCR product according to the standard protocol (Sambrook et al., 1989). One microliter aliquots of the PCR products were mixed with 4  ml of dye  95% formamide,  25 mM EDTA, 0.025% xylene-cyanol  and 0.025% bromophenol blue), incubated at 95ºC for 10 min and then chilled on ice. Denatured DNA was loaded on a 6% acrylamid gel in 1×TBE buffer and electrophoresed at a constant voltage of 150 V for 15 h. The gel was stained with 0.1% silver nitrate solution. (Xianyong et al., 2007). The PCR fragments from different SSCP patterns in different breeds were cleaned up and sequenced by the ABI PRIZM 377 DNA sequencer (Perkin-Elmer Corporation).

Statistical analysis: The genotypic and allelic frequencies, expected mean, observed and Nei’s heterozygosities were calculated using the Pop Gene 32 (version 1.31) software. The Hardy-Weinberg equilibrium in the populations was also tested. The sequences were aligned using the Sequencer Software (version. 4.2). The General Linear Model (GLM) procedure of Statistical Analysis System SAS (SAS Institute Inc, version 9.2, USA) was used to test the association of SNP marker genotypes of the LEP gene with growth and carcass traits. The initial and final live weights of the lambs were used as co-variate for body and carcass weight, while empty body weight was used as a co-variate for the carcass components.
The linear model used was as follows:
 
Yijk = m + Gi + b (xi+x) + eijk                                 (1)

Where Yijk is the considered dependent variable (growth and carcass traits); m, the overall mean, Gi, the fixed effect of the ith genotype at the LEP locus (I = 1, 2, 3 for AA, AG and GG, respectively); b (xi-xij¯), the covariate and eijk is the residual effect. Analyses were performed for each breed, separately.


RESULTS

PCR-SSCP and sequencing: The 260 bp PCR product consisting of exon 2 and part of intron 2 of the LEP gene was amplified. Following PCR-SSCP and sequencing, two alleles (A and G) were found and two genotypes (AA and AG) were observed. (Fig. 1). Two SNPs in intron 2 were noticed by A®G transition at the 113 bp position and T®C transition at the 165 bp position of the ovine LEP gene (Fig. 2). According to the results of sequencing, C®T SNP for all samples was similar. 

Genetic variability: The LEP allelic and genotypic frequencies and chi-square (χ2) values are summarized in Table 1. The genotypic frequencies of AA and AG were 0.53, 0.47 in the Shal lambs and 0.70, 0.30 in the Zandi lambs and 0.65, 0.35 in the Zel lambs, respectively. The GG genotype was not observed in these populations. Allelic frequencies for A and G were 0.76, 0.24 in the Shal lambs and 0.85, 0.15 in the Zandi lambs and 0.82, 0.18 in the Zel lambs, respectively. The χ2 test confirmed the Hardy-Weinberg equilibrium in the studied populations. Observed and expected heterozygosities, Nei’s and average heterozygosities of the LEP gene for the three Iranian sheep populations were slightly low (Table 2).

Association analysis: The GLM procedure was performed to determine possible association between the LEP gene and growth and carcass traits. Significant associations were found for the A®G SNP (Table 3). In the Shal breed, the presence of the A®G SNP is associated with an increase in cold carcass weight (2473 g; p< 0.01), fat-tail (3.3%; p< 0.05) and total body fat (1074.1 g; p< 0.05) of the carcass traits. In the Zel breed, the presence of the A®G SNP is associated with an increase in fat-tail (1.6%; p< 0.05) and a reduction in slaughter weight (-577.9 g; p < 0.05), cold carcass weight (577.4 g; p< 0.01) and lean meat weight (303.3 g; p< 0.05). This study did not indicate any obvious association between the intron 2 SNP and growth traits in the Shal breed.
DISCUSSION

Genetic variability: We observed two SNP in inron 2 in these three Iranian sheep populations whereas, three SNPs, A®G, C®T and C®G in intron 2 of the LEP gene in the Suffolk and Dorset breeds have been detected (Boucher et al., 2006).
    Two genotypes have been identified (AA and AG) with frequencies ranging from 0.53 to 0.70 for AA and 0.30 to 0.47 for AG genotype in Shal, Zandi, and Zel lambs. Results of this study are to some extent almost similar to published data on other sheep breeds. The genotypic frequencies of AA and AG in the ovine LEP gene of Dorset and Suffolk breeds have been found to be 0.75, 0.25 and 0.87, 0.13, respectively (Boucher et al., 2006).
    The frequencies of alleles were (0.76 and 0.24), (0.85 and 0.15) and (0.82 and 0.18) for A and G, in the Shal, Zandi and Zel lambs, respectively. Boucher et al. (2006) found allelic frequencies of 0.87 and 0.13 for A and G in Dorset sheep, respectively, which is in agreement with our results. Allelic frequencies of A and T for E2JW SNP in the bovine LEP gene of the Charolay cattle have been found to be 0.91 and 0.09, respectively (Schenkel et al., 2005). However, our results of sequencing suggested that C®T mutation is fixed in all studied population.
    The studied populations show a low degree of genetic variability for the LEP gene. Our explanation is that a few number of ram which have been used for breeding may lower the effective number of population with conservation so that heterozygosity decreases as a result of inbreeding.

Association analysis: We observed significant association in the presence of A®G SNP with some carcass traits in Shal and Zel breeds. Published data indicated a similar association between lean meat and the LEP locus (Boucher et al., 2006). The UASMS1 SNP is associated with fat yield, grade fat and lean meat yield (Schenkel et al., 2005). It has previously been reported that the A®G SNP in the ovine LEP gene has no significant association with total dissected fat (Boucher et al., 2006). Although no significant association was observed between A®G SNP and carcass traits in the Zandi breed, but a significant reduction (1255.3g) in weaning weight was detected (p< 0.01). It has been reported that intronic mutations, such as the A®G variant, may affect gene regulation and transcription levels (Tokuhiro et al., 2003). Intronic mutations can also result in splicing abnormalities, which often change the structure of mature protein (Faustino et al., 2003).


CONCLUSION
 
This study demonstrated the polymorphism of the LEP gene and its association with growth and carcass traits in three Iranian sheep breeds. The results obtained from this study suggest a significant association between the A113G variant of the LEP gene and carcass traits in the Shal and Zel breeds only. Generally, the lambs with the G allele had a higher carcass fat. However, the A®G SNP in intron 2 of the LEP gene did not have such an effect on growth traits and only in the Zandi breed a significant reduction in weaning weight was observed. To date, this was the first study that attempted to detect allele variation in the ovine LEP gene in Iranian sheep breeds. Further studies will be needed before application of this SNP to the breeding industry.

Acknowledgments

This study was conducted at the Research Station of Aboureyhan College of Agriculture, University of Tehran and the National Institute of Genetic Engineering and Biotechnology (NIGEB), Tehran, Iran. Authors gratefully acknowledge their excellent cooperation in this project.

AOAC (1995). Official methods of analysis AOAC, Association of Official Analytical Chemists. Washington, DC, USA.
Barb CR, Hausman GJ, Houseknecht KL (2001). Biology of leptin in the pig. Domestic Animal Endocrinology. 21: 297-317.
Boucher D, Palin MF, Castonguay F, Gariépy C, Pothier F (2006). Detection of polymorphisms in the ovine leptin (LEP) gene: Association of a single nucleotide polymorphism with muscle growth and meat quality traits. Can J Anim Sci. 86: 31-35.
Bouchanan FC, Fitzsimmons CJ, van Kessel AG, Thue TD, Winkelman-Sim DC, Schmutz SM (2002). Association of a missense mutation in the bovine leptin gene with carcass fat content and leptin mRNA levels. Genet Sel Evol. 34: 105-116.
Dela BFC, Shan B, Chen JL (1996). Identification of the promoter of the mouse obese gene. Proc Natl Acad Sci USA. 93: 4096-4101.
Dyer CJ, Simmons JM, Matteri RL, Keisler DH (1997). Leptin receptor mRNA is expressed in ewe anterior pituitary and adipose tissues and is differentially expressed in hypothalamic regions of well- fed and feed-restricted ewes. Domest Anim Endocrinol. 14: 119-128.
Faustino NA, Cooper TA (2003). Pre-mRNA splicing and human disease. Genes Dev. 17: 419-437.
Fiona CB, Carolyn J., Andrew GV, Tracey DT, Dianne CW, Sheila MS (2002). Association of a missense mutation in the bovine leptin gene with carcass fat content and leptin mRNA levels. Genet Sel Evol. 34: 105-116.
Halaas JL, Gajiwala KS, Maffei M, Cohen SL, Chait BT, Rabinowitz D, Lallone RL, Burley SK, Friedman JM (1995). Weight-reducing effects of the plasma protein encoded by the obese gene. Science 269: 543-546.
Houseknecht KL, Baile CA, Matteri RL, Spurlock ME (1998). The biology of leptin: a review. J Anim Sci. 76: 1405-1420.
Kashan NEJ, Manafi-Azar GH, Afzalzadeh A, Salehi A (2005). Growth performance and carcass quality of fattening lambs from fat-tailed and tailed sheep breeds. Small Rumin Res. 60: 267-271.
Nei M (1973). Analysis of gene diversity in subdivided populations. Proc Natl Acad Sci. 70: 3321-3323.
Ohshiro Y, Ueda K, Nishi M, Ishigame M, Wakasaki H, Kawashima H, Furuta H, Sasaki H, Sanke T, Takasu N, Nanjo K (2000). A polymorphic marker in the leptin gene associated with Japanese morbid obesity. J Mol Med. 78: 516-520.
Ozaki K, Ohnishi Y, Iida A, Sekine A, Yamada R, Tsunoda T, Sato H, Hori M, Nakamura Y, Tanaka T (2002). Functional SNPs in the lymphotoxin-alpha gene that are associated with susceptibility to myocardial infarction. Nat Genet. 32: 650-654.
Sambrook J, Fritsch EF, Maniatis T (1989). Molecular cloning: A laboratory manual, 2nd, p. B.23. Cold Spring Harbor Laboratory Press, Cold Spring Harbor, New York. PP. 1350-1356.
Schenkel FS, Miller SP, Ye X, Moore SS, Nkrumah JD, Li C, Yu J, Mandell IB, Wilton JW, Williams JL (2005). Association of single nucleotide polymorphisms in the leptin gene with carcass and meat quality traits of beef cattle. J Anim Sci. 83: 2009-2020.
Schneider JE, Zhou D, Blum RM (2000). Leptin and metabolic control of reproduction. Horm Behav. 37: 306-326.
Taouis M, Chen JW, Davi aud C, Dupont J, Derouet M, Simon J  (1998). Cloning the chicken leptin gene. Gene. 208: 239-242.
Thomas L, Wallace JM, Aitken RP, Mercer JG, Trayhurn P, Hoggard N (2001). Circulating leptin during ovine pregnancy in relation to maternal nutrition, body composition and pregnancy outcome. J Endocrinol. 169: 465-476.
Tokuhiro S, Yamada R, Chang X, Suzuki A, Kochi Y, Sawada T, Suzuki M, Nagasaki M, Ohtsuki M, Ono M, Furukawa H, Nagashima M, Yoshino S, Mabuchi A, Sekine A, Saito S, Takahashi A, Tsunoda T, Nakamura Y, Yamamoto K (2003). An intronic SNP in a RUNX1 binding site of SLC22A4, encoding an organic cation transporter is associated with rheumatoid arthritis. Nat Genet. 35: 341-348.
Zhang Y, Proenca R, Maffei M, Barone M, Leopold L, Friedman JM (1994). Positional cloning of the mouse obese gene and its human homologue. Nature 372: 425-432