Livestock Research for Rural Development 32 (5) 2020 | LRRD Search | LRRD Misssion | Guide for preparation of papers | LRRD Newsletter | Citation of this paper |
Local sheep breeds play an important role for the Albanian farmers as the mains source of milk, meat and wool. The aim of this study was to detect the genetic polymorphism of exon 17 in alpha -casein gene and exon 4, in kappa -casein gene in three Albanian local sheep breeds: Bardhoka, Ruda and Shkodrane, by the use of PCR-SSCP and PCR-RFLP respectively. The sequence analysis showed a single nucleotide polymorphism (SNP) at position 137 (C>T) and at 168 (T>C) of exon 17 CSN1S1 and exon 4 CSN3, respectively. Analysis of molecular variance (AMOVA) showed that 4% of differentiation was between populations. Genetic differentiation among populations (FST) is low, only 4.2%. Gene flow between populations was 5.97.
Keywords: genetic polymorphism, gene flow, AMOVA
Local sheep breeds are the main livestock species reared by Albanian farmers, which are managed in extensive and semi extensive conditions and are the main source of milk, meat and wool mainly for family consumption (Hoda and Marsan 2012).
Milk proteins include caseins and whey proteins (Farrell Jr et al 2004). Caseins are a family of phosphoproteins synthesized in the mammary gland in response to lactogenic hormones and secreted as micelles (Feligini et al 2005). Casein genetic polymorphisms are important and well known due to their effects on quantitative traits and technological properties of milk (Ceriotti et al 2004).
The genetic polymorphism of casein genes is studied in different Egyptian breeds (Othman et al 2013a,b), Czech (Sztankoova et al 2012), Italian (Sacchi et al 2005) etc. Feligini et al (2005) have studied SNP in kappa casein coding gene of Pag, Sarda and Pramenka sheep breeds. Genetic polymorphism of CSN1S1 casein gene is studied also in Bulgarian sheep breeds. (Hristova at al 2011, Gencheva and Georgieva 2019).
Previously different markers are used to characterize the local sheep breeds in Albania like microsatellite (Hoda and Marsan 2012, Hoda et al 2010b), AFLP (Hoda et al 2010a), mtDNA (Hoda et al 2014), SNP (Hoda et al 2011), which have shown that these breeds face a high level of admixture. Milk protein polymorphisms have been considered as potential tools for selection of dairy sheep breeds (Ceriotti et al 2004).
The aim of this study was to investigate allele frequencies at casein loci of milk protein genes (exon 17 CSN1S1and exon 4 CSN3) in three Albanian local sheep breeds, Bardhoka, Ruda and Shkodrane.
A total of 93 individuals from three local sheep breeds: Bardhoka, Ruda and Shkodrane were included in the analysis. A description of the breed is provided previously (Hoda and Marsan 2012). Blood samples were collected from 31 unrelated individuals, per each breed, according to the recommendations of ISAG/FAO. DNA was isolated from the whole blood using the standard procedure of phenol-chloroform extraction and Proteinase K digestion. The exon 17 from ovine CSN1S1 gene was analyzed using PCR - SSCP method, meanwhile exon 4 of CSN3 gene was determined by PCR- RFLP analysis as described by Pariset et al (2006).
GenAlEx 6.5 program (Peakall and Smouse 2012) was used to compute the allele frequency per locus, observed and expected heterozygosity, Information index. ARLEQUIN software (Excoffier et al 2005) was used for the analysis of molecular variance (AMOVA). GENEPOP software (Rousset 2008) was used to calculate fixation index F (Weir and Cockerham 1984).
Sequence analysis of CSN1S1exon 17 in three local sheep breeds from Albania revealed transition C to T, at position 137. The C SNP has the highest frequency at Bardhoka breed (78.6%; Figure 1). The analysis of Exon 4 of CSN3gene detected a SNP at position 168 (T>C).The highest frequency of T was found at Bardhoka breed (92.9%; Figure 1). Ruda breed displayed the highest values of genetic diversity at each locus, while Bardhoka showed the lowest value for the CSN3 locus.
Figure 1. Allelic frequencies per locus and breed |
Table 1. Information index, diversity and unbiased diversity by locus and populations |
||||||||||
Bardhoka |
Ruda |
Shkodrane |
||||||||
CSN1S1 |
CSN3 |
MLP |
CSN1S1 |
CSN3 |
MLP |
CSN1S1 |
CSN3 |
MLP |
MLP |
|
I |
0.340 |
0.154 |
0.247 |
0.596 |
0.534 |
0.565 |
0.510 |
0.401 |
0.456 |
0.423 |
Ho |
0.214 |
0.071 |
0.143 |
0.433 |
0.194 |
0.313 |
0.207 |
0.207 |
0.207 |
0.221 |
He |
0.191 |
0.069 |
0.130 |
0.406 |
0.350 |
0.378 |
0.328 |
0.238 |
0.283 |
0.264 |
uHe |
0.195 |
0.070 |
0.132 |
0.413 |
0.355 |
0.384 |
0.334 |
0.242 |
0.288 |
0.268 |
MLP = Mean over loci, per pop |
According to Takezaki and Nei (1996), a marker can be used in the selection program if the average heterozygosity in a breed is in the range of 0.3 to 0.8. Based on the results obtained here, both loci cannot be considered suitable for selection, however, they could be used for estimation of breed genetic diversity.
Table 2. F-statistics and Estimates of Nm over All Pops for each Locus |
||||
Locus |
Fis |
Fit |
Fst |
Nm |
CSN1S1EX17 |
0.077 |
0.107 |
0.033 |
7.411 |
CSN3_EX4 |
0.281 |
0.319 |
0.052 |
4.532 |
Mean |
0.179 |
0.213 |
0.042 |
5.972 |
F statistics and gene flow (Nm), per each locus, are shown in Table 2. The data show that the global heterozygote deficit across populations (FIT) is 21.3%. The heterozygote deficit due to inbreeding within populations (FIS) is 17.9%. Genetic differentiation among populations (FST) is low, only 4.2%. The low level of genetic differentiation is supported by gene flow across breeds over two loci, which is 5.97.
Table 3. AMOVA results over 2 loci |
|||
Source |
Sum of |
Variance |
Variation |
Among Populations |
1.965 |
0.022 |
3.98 |
Among animals |
27.22 |
0.049 |
17.44% |
Within animals |
19.500 |
0.223 |
78.65% |
Total |
48.687 |
0.284 |
|
The analysis of molecular variance (AMOVA) displayed 3.98% differentiation among populations, similar as explained by FST values(Table 3). AMOVA shows that only 4% of total variation is found among the population.
The substitution of nucleotide between variant C to T at position 137 of exon 17 change the aminoacid Thr137 with Ile137. At exon 4 of CSN3 gene, thymine was substituted for a cytosine while the corresponding amino acid remains a tyrosine. The same finding has shown Feligini et al (2005), for Pag, Pramenka and Sarda breeds, from Croatia, Serbia and Italy, respectively.
The use of genetic polymorphism of casein genes in molecular marker-assisted selection (MAS) is of great interest to improve milk productivity in farm animals (Kumar et al 2006). The genetic polymorphisms of milk proteins are often associated with quantitative and qualitative parameters in milk being potential candidate markers that should be included in breeding strategies (Selvaggi et al 2014). Identification of polymorphisms in ovine milk is more important because sheep milk is commonly used solely for cheesemaking (Feligini et al 2005). Therefore the results obtained here can be used also for an association between genetic variants and ovine milk traits, as previously done by Corral et al (2010) and to consider the molecular markers in selection schemes for milk production in sheep breeds.
Ceriotti G, Chessa S, Bolla P, Budelli E, Bianchi L, Duranti E and Caroli A 2004 Single nucleotide polymorphisms in the ovine casein genes detected by polymerase chain reaction-single strand conformation polymorphism, Journal of dairy science. 87 (8): 2606-2613.
Corral J M, Padilla J A and Izquierdo M 2010 Associations between milk protein genetic polymorphisms and milk production traits in Merino sheep breed, Livestock Science. 129 (1-3): 73-79.
Excoffier L, Laval G and Schneider S 2005 Arlequin (version 3.0): an integrated software package for population genetics data analysis, Evolutionary bioinformatics. 1, p. 117693430500100003.
Farrell Jr H M, Jimenez-Flores R, Bleck G T, Brown E M, Butler J E, Creamer L K, Hicks C L, Hollar C M, Ng-Kwai-Hang K F and Swaisgood H E 2004 Nomenclature of the proteins of cows milk—Sixth revision, Journal of dairy science. 87 (6): 1641-1674.
Feligini M, Vlaco S, Curik V C, Parma P, Greppi G and Enne G 2005 A single nucleotide polymorphism in the sheep $$-casein coding region, Journal of dairy research. 72 (3): 317-321.
Gencheva D and Georgieva S 2019 Genetic diversity and population structure of Bulgarian autochthonous sheep breeds based on nucleotide variation in Alpha S1-casein gene, Bulgarian Journal of Agricultural Science. 25, p 3.
Hoda A and Marsan P A 2012 Genetic characterization of Albanian sheep breeds by microsatellite markers, Analysis of genetic variation in animals 3-27.
Hoda A, Ajmone-Marsan P, Hykaj G and Consotrium E 2010a Genetic diversity in albanian sheep breeds estimated by AFLP markers, Albanian Journal of Agricultural Sciences. 9 (2): 23-29.
Hoda A, Bytyqi H and Bozgo V 2014 Genetic diversity of sheep breeds from Albania and Kosovo revealed by mtDNA, Journal of Endocytobiosis and Cell Research. 66, p. 69.
Hoda A, Hykaj G and consortium E 2010b Genetic diversity and population structure in albanian local sheep breeds analyzed by microsatellite markers, Albanian Journal of Agricultural Sciences. 9 (1): 55-62.
Hoda A, Hykaj G, Sena L and Delia E 2011 Population structure in three Albanian sheep breeds using 36 single nucleotide polymorphisms, Acta Agriculturae Scand Section A. 61 (1): 12-20.
Hristova D and others 2011 Genetic polymorphism of alpha S1-casein gene in Bulgarian sheep breeds, Agricultural Science and Technology. 3 (1): 8-12.
Kumar D, Gupta N, Ahlawat S, Satyanarayana R, Sunder S and Gupta S C 2006 Single strand confirmation polymorphism (SSCP) detection in exon I of the a-lactalbumin gene of Indian Jamunapari milk goats (Capra hircus), Genetics and Molecular Biology. 29 (2): 287-289.
Othman O E, El-Fiky S A, Hassan N A, Mahfouz E R and Balabel E A 2013a Genetic polymorphism detection of two alpha-Casein genes in three Egyptian sheep breeds, Journal of Genetic Engineering and Biotechnology. 11 (2): 129-134.
Othman O E, El-Fiky S A, Hassan N A, Mahfouz E R and Balabel E A 2013b Genetic variations of B beta and K-casein genes in Egyptian sheep breeds, Journal of Applied Biosciences. 64 (1): 4858-4866.
Pariset L, Cappuccio I, Ajmone-Marsan P, Bruford M, Dunner S, Cortes O, Erhardt G, Prinzenberg E-M, Gutscher K, Joost S and others 2006 Characterization of 37 breed-specific single-nucleotide polymorphisms in sheep, Journal of Heredity. 97 (5): 531-534.
Peakall R and Smouse P E 2012 GenAlEx 6.5: genetic analysis in Excel. Population genetic software for teaching and researchdan update. Bioinformatics 28, 2537e2539, .
Rousset F 2008 genepop007: a complete re-implementation of the genepop software for Windows and Linux, Molecular ecology resources. 8 (1): 103-106.
Sacchi P, Chessa S, Budelli E, Bolla P, Ceriotti G, Soglia D, Rasero R, Cauvin E and Caroli A 2005 Casein haplotype structure in five Italian goat breeds, Journal of Dairy Science. 88 (4): 1561-1568.
Selvaggi M, Laudadio V, Dario C and Tufarelli V 2014 Investigating the genetic polymorphism of sheep milk proteins: a useful tool for dairy production, Journal of the Science of Food and Agriculture. 94 (15): 3090-3099.
Sztankoova Z, Matlova V and Rychtarova J 2012 Genetic polymorphism at the milk protein genes (CSN1S1, CSN2 and CSN3) in the Czech Sumava and Walachian sheep breeds, Grant Journal. 1 (2): 63-66.
Takezaki N and Nei M 1996 Genetic distances and reconstruction of phylogenetic trees from microsatellite DNA, Genetics. 144 (1): 389-399.
Weir B S and Cockerham C C 1984 Estimating F-statistics for the analysis of population structure, evolution. 38 (6): 1358-1370.
Received 19 February 2020; Accepted 30 March 2020; Published 1 May 2020