Advances in dairy Research

Advances in dairy Research
Open Access

ISSN: 2329-888X

+44 1300 500008

Research Article - (2015) Volume 3, Issue 1

Effect of SNPs in Prolactin Promoter on Milk Traits in Egyptian Buffalo

Mohammad Abu El- Magd1*, Abeer Gamal Zaki2, Khaled Abdelalem Kahilo3, Mohamed El-Sayd Barakat2 and Ibrahim Fatouh Hassan3
1Department of Anatomy and Embryology, Faculty of Veterinary Medicine, Kafrelsheikh University, Egypt
2Department of Biochemistry, Animal Health Research Institute, Kafrelsheikh, Egypt
3Department of Biochemistry, Faculty of Veterinary Medicine, Kafrelsheikh University, Egypt
*Corresponding Author: Dr. Mohammad Abu El- Magd, Department of Anatomy and Embryology, Faculty of Veterinary Medicine, Kafrelsheikh University, Kafrelsheikh, Egypt, Tel: +201068160301, Fax: +20473215497 Email:

Abstract

Prolactin (PRL) plays a crucial role in the initiation and maintenance of lactation in animals. It is a potential quantitative trait locus and genetic marker of production traits in dairy cattle. However, no SNPs was detected in, and no association studies was conducted on, PRL gene in Egyptian buffalo. Therefore, this study aims to detect PRL gene polymorphisms and to study their effect on in Egyptian water buffalo. A PCR product of 245 bp, including part of the promoter of PRL was amplified using PCR and subsequently, SSCP and nucleotide sequencing were conducted to identify its different allelic patterns in Egyptian buffaloes. Our results revealed presence of a novel T72C SNP in the PRL promoter with two genotypes (TC, CC). Chi-square test indicated presence of genetic disequilibrium with respect to studied variations and revealed significant difference in distribution of variants among buffaloes. Genetic diversity parameter showed moderate degree of polymorphism. The two genotypes (TC, CC) did not show association with milk traits (milk yield, fat%, protein%, and lactose%, solid %) in Egyptian buffaloes. We conclude that no significant difference in milk traits between Egyptian buffaloes have TC or CC genotypes and so breeders should not depend on this SNP (alone) in their selection for dairy animals.

Keywords: Prolactin; Egyptian buffalo; PCR; Sequencing; SNPs

Introduction

As a result to the actual increase in the world population and shortage of food supply, it becomes necessary for developing countries to maximize the production of their native animals to guarantee a sustainable source of food of animal origin. Among domestic animals, the water buffalo (Bubalus bubalis ), particularly the river buffalo, which is the main buffalo breed in Egypt, holds great promise and potential for animal production. Indeed, the Egyptian buffalo plays a crucial role in Egypt economy by their livestock productivity from milk and meat. It contributes about 2.7% and 8.4% to the world buffalo's milk and meat, respectively [1]. The great importance of the Egyptian buffalo (Bubalus bubalis ) and its superiority to the domesticated cattle is owed to their great adaptive capacity to the tropical climate, excellent nutritional benefits, resistance to the diseases and also the Egyptian consumers prefer buffalo milk to cattle milk because its white colour and good flavor. In addition, buffalo milk is characterized by high milk fat (7%), solid not fat (16%), lower levels of phospholipids and cholesterol and the fat has a higher proportion of saturated fatty acids, so that Egyptian buffalo consider the main dairy animal in Egypt [2].

To improve animal traits and productivity, the animals should be selected periodically on the basis of the genotype which will reflect on their productivity and on the national economics. There is a considerable interest in the application of molecular genetics technologies in the form of specific DNA markers that are associated with various productivity traits to promote more efficient and relatively easy selection and breeding of farm animals with an advantage for inheritable traits of meat and milk productivity. Many candidate genes have been identified and selected for analysis based on a known relationship with productivity traits [3-6]. The candidate genes that related to milk traits are called lactogenic genes. Among different candidates, the prolactin (PRL) gene seems to be promising, because it plays a crucial role in mammary gland development and in the initiation, maintenance of lactation and expression of milk protein genes [7]. Prolactin is a single chian polypeptide hormone secreted from specialized lactotroph cells of the anterior pituitary gland, furthermore, it is produced by numerous other cells and tissues, including the mammary gland [8,9]. The main functions of PRL are regulation of reproduction, promotion of lactation in mammals, synthesis of milk (lactogenesis) and maintenance of milk secretion (galactopoiesis) [10].

The bovine PRL gene was mapped on chromosome 23 [11] and consists of 5 exons and 4 introns with 10 kb in size [12,13], encoding a 229-amino-acid prolactin precursor. The signal peptide contains 30 amino acids; thus the mature bovine PRL is composed of 199 amino acids [14]. Several polymorphic sites have been detected within PRL gene and statistically significant associations between PRL variants and milk production traits have been described in dairy cattle. A previous study has reported the associations between polymorphisms in the coding regions of the bovine PRL and milk production traits [15]. Scanty researches have been conducted exploring the genetic diversity on molecular genetic basis in buffalo all over the world in comparison with other farm animals. Moreover, only very few studies were conducted on the Egyptian buffalo. Therefore, the aim of this study is to identify PRL SNPs in Egyptian water buffalo and to study their association with milk traits.

Materials and Methods

Sampling and DNA extraction

This study involved 200 pure Beheiry Egyptian water buffalo (Bubalus bubalis ) kept on a farm located in Nucleus Herd, Nataff Gedeed station of MahaletMousa farm, Kafrelsheikh governorate. All records of milk production traits in Beheiry buffaloes were collected from farm records to be used for statistical analysis. Records covered the period from 2012 till 2014. Animals were kept under semi-open sheds. Lactating buffaloes were milked by hand twice daily at 7.00 a.m. and 5.00 p.m. throughout the lactation period, and milk production was recorded daily. Buffaloes were maintained under the same system of feeding in the farm. The animals were grazed on Egyptian clover (Trifolium Alexandrinum ) during December to May with concentrate mixture and rice straw. While, during June to November, animals were fed on concentrate mixture, rice straw and limited amount of clover hay. The concentrate feed mixture was given twice daily before milking, while rice straw was offered once daily at 9.00 a. m, whereas clover hay or (silage) in the summer was offered at 11.00 a.m. Multi mineral licking blocks were available for animals in the stalls.

Blood samples were collected by jugular vein puncture into vacutainer tubes contain an anticoagulant (disodium EDTA) and kept in ice box then kept at -20 for further use. The genomic DNA was extracted using Gene JET genomic DNA extraction kit following the manufacturer protocol (Jena Bioscience, Jena/Germany ).

Polymerase chain reaction (PCR)

A partial sequence of prolactin gene promoter was amplified by PCR using forward primer 5\- CAGCCAGCAATTTTGATGA-3\ and reverse primer 5\- GATGTTCATTTCTGGTCAGTATG-3\ designed by Primer 3.0 software based on the published sequences of Indian buffalo (GenBank accession number, AF426315). The PCR was carried out in a reaction volume of 50 μL, containing 4.0 μL DNA template (approximately 100 ng), 1.0 μL (0.20 mM) dNTP, 5.0 μL buffer, 3.0 μL (2.5 mM) MgCl2, 2.0 μL 10 μmol/L forward primer, 2.0 μL 10 μmol/L reverse primer, 1.0 μL 10X Taq DNA polymerase (5 U/μL, Fermentas, #K1071, European Union), and 32 μL nuclease free water. Thermal cycling parameters were as follows: initial denaturation at 94°C for 5 min, 35 cycles of amplification (94°C for 40 s for DNA denaturation, annealing 56°C for 1 min, extension at 72°C for 1 min) and final extension at 72°C for 10 min. The samples were held at 4°C. PCR products (245bp) were resolved by electrophoresis on 2% agarose gel in 1X TBE, stained with ethidium bromide and visualized with UV light of gel documentation system (Biometra Biomedizinische Analytik, GmbH).

Single stranded conformational polymorphism (SSCP)

Five μL of PCR products were mixed with 5μL of denaturating dye (95% formamide, 25 mM EDTA, 0.025% xylene–cyanole and 0.025% bromophenol blue), heated in 95°C for 10min and chilled in ice then subjected in 12% PAGE (37.5:1 acrylamide: bis-acrylamide) in 1x TBE buffer at 200V for 6h at 4°C. The gel was stained with ethidium bromide as was previously described by [4].

DNA Sequencing

PCR products with expected size were purified using PCR purification kit following the manufacturer protocol (Jena Bioscience, # pp-201×s) to remove primer dimmers, primers, nucleotides, proteins, salt, agarose, ethidium bromide and other impurities. The purified PCR products were sent to MacroGen Company (South Korea) to sequence in both directions using ABI 3730XL DNA sequencer (Applied Biosystem, USA) and the identity of the sequenced PCR product was examined using Blast search against the published sequences of Indian buffalo on GenBank database. The alignments and assembly of the sequences were performed using Geneious 4.8.4 software.

Statistical Analysis

Allele and genotype frequencies were calculated by direct counting. Hardy–Weinberg equilibrium (HWE) assessed by the χ2 test and population genetic indices, i.e. He (gene heterozygosity) and Ne (effective allele numbers; reciprocal of homozygosity) were computed by POPGENE software [16]. Polymorphism information content (PIC) was calculated and classified into low polymorphism (PIC value <0.25), median polymorphism (0.250.5) according to Botstein et al. [17]. GLM was used to analyse the relationship between variations of the PRL promoter and milk traits using software SAS 6.12 and the following model: Ykljm= μ + Ak + Gl +Dj+ Eklm, where Ykljm was the trait measured on each of the klmth animals μ was the overall population mean Ak was fixed effect due to the kth age (in days), Gl was the fixed effect associated with lth genotype, Dj = covariate of the days in milking at the test day, and Ekljm was the random error. All data was described as least square means ± standard error of means (SEM). The significance of least square means was compared for significance using Duncan’s test.

Results and Discussion

Detection of novel T72C SNP in PRL promoter

Previous studies have considered the bovine PRL gene as an excellent candidate for linkage analysis with quantitative trait loci (QTL) affecting milk production traits [18-20]. Allelic variation in the structural or regulatory sequences of the PRL gene would be of interest because of the possible direct and indirect effect on milk production. In the current study, a locus containing a partial sequence of buffalo PRL promoter determined by PCR was examined to search for any prospective polymorphisms (Figure 1).

advances-dairy-research-ethidium-bromide

Figure 1: Ethidium bromide stained agarose gel of PCR products representing amplification of 245bp PRL locus (a part of the promoter) in 12 Egyptian buffaloes. M represents 100bp ladder.

Nucleotide sequences of this locus showed a novel T72C SNP (Figures 2A and 2B). In this SNP, the cytosine (C) nucleotide number 72 in PRL locus was replaced by thymine (T) nucleotide. SSCP patterns of PRL showed two different SSCP patterns; TC (three bands) and CC (one band) (Figure 3).

advances-dairy-research-nucleotide-sequences

Figure 2a: Whole nucleotide sequences of PRL.1 locus showed T72C SNP (small colored box) in one Egyptian buffalo.

advances-dairy-research-egyptian-buffaloes

Figure 2b: Nucleotide sequences of PRL showed a T72C SNP in 2 Egyptian buffaloes. The arrow indicates the position of the SNP.

advances-dairy-research-three-bands-pattern

Figure 3: SSCP patterns of PRL (promoter) showed two different bands patterns: three bands pattern in TC genotype (lanes 1, 2, 4-6) and one band pattern in CC genotype (lane 3) in 6 Egyptian buffaloes.

To our knowledge, this is the first association study applied on PRL gene of Egyptian buffalo. This SNP was novel and so was not overlapped with mutations reported in buffaloes or other animals. In agreement, some other SNPs were determined in the promoter regions of PRL in other animals: A1043G, A402G, A767C, G485T, C247A SNPs [19], and T175G and A446G SNPs [21] in cattle in addition to C499T SNP in pig [20].

Analysis of allele, genotype frequencies, Chi-Square (χ2) and genetic indices

The genotype frequencies of T72C were 0.37 (74) for CC and 0.63 (126) for CT, with allele frequencies of 0.3150 for C and 0.6850 for T. Chi-Square (χ2) equals 41.9 (P<0.05) among the 200 animals that were genotyped for this polymorphism. This means that all buffaloes were deviated from HWE. This deviation from HWE suggests a change in the distribution of alleles from one generation to the next.

The value of the difference between expected and observed He was 0.431, while the value of Ne was 1.759. The PIC values were medium (0.338). The moderate He and Ne of T72C SNP shows the moderate polymorphism, moderate genetic diversity, moderate ability to maintain allelic stability during selection or mutation. Moreover, the moderate PIC (moderate polymorphism) is suitable as marker in molecular breeding [22] and indicates absence of genetic drift, no isolation or population substructure, or that they had not passed through artificial selections with a greater degree in history. This is in agreement with the results from Hardy-Weinberg which also show deviation from HWE.

Association between T72C SNP with milk trait

SNPs occurring within the PRL gene may influence the chemical composition of milk or at least be an effective DNA marker of a subregion of dairy cattle genome [18]. Although the SNPs in this study were detected in non-coding regions, they may be associated with production traits as reported by previous studies [15,19,21]. One of these studies has identified two SNPs in bovine PRL promoter, A1043G and A402G that statistically associated with milk yield and fat content [15]. In addition, three other SNPs (A767C, G485T, and C247A) were identified within the promoter of the bovine PRL gene that significantly associated with milk yield, fat yield, and protein yield and protein percentage of dairy cattle [19]. Another two SNPs, T175G and A446G, also associated with milk traits were detected in the 5'UTR of the PRL promoter [21].

Therefore, we used the GLM to study the statistical association between the T72C SNP and some milk traits; including milk yield, fat %, protein %, lactose % and solid %, in Egyptian dairy buffaloes. The results of this association analysis were detailed in Table 1. The T72C SNPs in PRL promoter showed no significant association with any of the examined milk traits (P>0.05).

Traits  TC  CC
Lactation length 259.35±6.43 261.63±6.14
Milk yield 2042.52±68.29 2063.62±58.31
Fat percentage 6.20±0.02 6.46±0.03
Protein percentage 4.21±0.03 4.10±0.03
Lactose percentage 5.45±02 5.19±02
Solid percentage 17.54±0.31 17.36±0.23

Table 1: Least squares mean (±SE) of different production traits for SSCP pattern of buffalo PRL. Lactation length is given in days and milk yield in Kg.

Although this SNP is not associated with milk traits, it may be linked with another SNP in the same gene or in another QTL related to milk traits. The phenotype effect of a single SNP may be affected by other SNPs in other loci in the same gene or other linked genes especially in polygenic traits such as milk traits. Thus, analyses of single SNP may be no effect on the phenotype, however when in combination, two SNPs or more displayed more profound impacts than in separation.

This agrees with the conclusion of that the inheritance of genotype combinations was more effective than that of a single SNP [23]. Therefore, further investigations including screening of the whole PRL gene sequences a long with statistical analysis of a large population size can shed the light on the significance of the detected SNP in this study in combination with other correlated SNPs and whether they associated with milk traits or not.

Conclusion

In this study, we detected a novel T72C SNP in the promoter region of the PRL gene, however, the statistical analysis showed no association between this SNP and some milk traits (milk yield, fat%, protein %, lactose % and solid %) in Egyptian dairy water buffaloes. Before excluding this SNP from breeding program for dairy buffaloes, it is necessary to extend our study to involve a large population size and whole sequences of the PRL gene. Considering the economic importance of the milk to the livestock industry, it appears clearly essential to further research on PRL in the buffalo.

Acknowledgements

This study was funded by Science Technology Development Fund (STDF), Ministry of Scientific Research, Egypt. All experiments were conducted in the Biotechnology Lab, Faculty of Veterinary Medicine, Kafrelsheikh University, Egypt.

References

  1. Abdel-Salam SA, Sayed AI, Elsayed M, Abou-Bakr S (2010) Genetic gain in open nucleus breeding scheme to improve milk production in Egyptian Buffalo. LivestSci 131: 162-167.
  2. Ahmad S, Anjum FM, Huma N, Sameen A, Zahoor T (2013) Composition and physico-chemical characteristics of Buffalo milk with particular emphasis on lipids, proteins,minerals, enzymess and vitamins. . The Journal of Animal and Plant Sciences 23: 62-74.
  3. Spelman RJ, Bovenhuis H (1998) Moving from QTL experimental results to the utilization of QTL in breeding programmes. Anim Genet 29: 77-84.
  4. Abo-Al-Ela HG, El-Magd MA, El-Nahas AF, Mansour AA (2014) Association of a novel SNP in exon 10 of the IGF2 gene with growth traits in Egyptian water buffalo (Bubalusbubalis). Trop Anim Health Prod 46: 947-952.
  5. El-Magd MA, Abbas HE, El-kattawy AM, Mokhbatly A (2013) Novel polymorphisms of the IGF1R gene and their association with average daily gain in Egyptian buffalo (Bubalusbubalis). DomestAnimEndocrinol 45: 105-110.
  6. El-Magd MA, Abo-Al-Ela HG, El-Nahas A, Saleh AA, Mansour AA (2014) Effects of a novel SNP of IGF2R gene on growth traits and expression rate of IGF2R and IGF2 genes in gluteus medius muscle of Egyptian buffalo. Gene 540: 133-139.
  7. Othman OE, Zayed FA, Ali A. El Gawead b, El-Rahman MRA (2011) Genetic polymorphism of three genes associated with milk trait in Egyptian buffalo. Journal of Genetic Engineering and Biotechnology 9: 97–102.
  8. Bole-Feysot C, Goffin V, Edery M, Binart N, Kelly PA (1998) Prolactin (PRL) and its receptor: actions, signal transduction pathways and phenotypes observed in PRL receptor knockout mice. Endocr Rev 19: 225-268.
  9. Perrot-Applanat M, Gualillo O, Pezet A, Vincent V, Edery M, et al. (1997) Dominant negative and cooperative effects of mutant forms of prolactin receptor. MolEndocrinol 11: 1020-1032.
  10. Freeman ME, Kanyicska B, Lerant A, Nagy G (2000) Prolactin: structure, function, and regulation of secretion. Physiol Rev 80: 1523-1631.
  11. Hallerman EM, Theilmann JL, Beckmann JS, Soller M, Womack JE (1988) Mapping of bovine prolactin and rhodopsin genes in hybrid somatic cells. Anim Genet 19: 123-131.
  12. Cooke NE, Coit D, Shine J, Baxter JD, Martial JA (1981) Human prolactin. cDNA structural analysis and evolutionary comparisons. J BiolChem 256: 4007-4016.
  13. Truong AT, Duez C, Belayew A, Renard A, Pictet R, et al. (1984) Isolation and characterization of the human prolactin gene. EMBO J 3: 429-437.
  14. Cao X, Wang Q, Yan JB, Yang FK, Huang SZ, et al. (2002) [Molecular cloning and analysis of bovine prolactin full-long genomic as well as cDNA sequences]. Yi ChuanXueBao 29: 768-773.
  15. Lü A, Hu X, Chen H, Jiang J, Zhang C, et al. (2010) Single nucleotide polymorphisms in bovine PRL gene and their associations with milk production traits in Chinese Holsteins. MolBiol Rep 37: 547-551.
  16. Yeh FC, Yang R, Boyle T (1999) POPGENE Version 3.1, Microsoft Window-based Freeware for Population Genetic Analysis. University of Alberta, Edmonton, AB, Canada.
  17. Botstein D, White RL, Skolnick M, Davis RW (1980) Construction of a genetic linkage map in man using restriction fragment length polymorphisms. Am J Hum Genet 32: 314-331.
  18. Brym P, Kami Åski S, WÃ jcik E (2005) Nucleotide sequence polymorphism within exon 4 of the bovine prolactin gene and its associations with milk performance traits. J Appl Genet 46: 179-185.
  19. He F, Sun DX, Yu Y, Wang YC, Zhang Y (2006) Association between SNPs within prolactin gene and milk performance traits in Holstein dairy cattle. Asian Austral J Anim 19: 1384-1389.
  20. Korwin-Kossakowska A, Wyszyńska-Koko J, Sender G (2006) A novel polymorphism of the porcine prolactin gene (PRL). NeuroEndocrinolLett 27: 241-246.
  21. Li JT, Wang AH, Chen P, Li HB, Zhang CS, et al. (2006) Relationship between the polymorphisms of 5 ' regulation region of prolactin gene and milk traits in Chinese Holstein dairy cows. Asian Austral J Anim 19: 459-462.
  22. Vaiman D, Mercier D, Moazami-Goudarzi K, Eggen A, Ciampolini R, et al. (1994) A set of 99 cattle microsatellites: characterization, synteny mapping, and polymorphism. Mammalian genome: official journal of the International Mammalian Genome Society 5: 288-297.
  23. Fallin D, Cohen A, Essioux L, Chumakov I, Blumenfeld M, et al. (2001) Genetic analysis of case/control data using estimated haplotype frequencies: application to APOE locus variation and Alzheimer's disease. Genome Res 11: 143-151.
Citation: Magd MAE, Zaki AG, Kahilo KA, Barakat MES, Hassan IF (2015) Effect of SNPs in Prolactin Promoter on Milk Traits in Egyptian Buffalo. J Adv Dairy Res 3:128.

Copyright: © 2015 Magd MAE, et al. This is an open-access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.
Top