Livestock Research for Rural Development 21 (6) 2009 | Guide for preparation of papers | LRRD News | Citation of this paper |
Sheep are considered as an important livestock species in Albania. An analysis of six microsatellite loci in 124 unrelated individuals was performed to define the genetic structure and variability of four local sheep breeds: “Bardhoka”, “Ruda”, “Recka” and “Shkodrane”. 31 individuals were selected per each breed. Allele diversity, observed heterozygosity, expected heterozygosit, GST estimates, DS genetic distance were calculated. Wright F-statistics, (FIS, FIT and FST) values were estimated.
A total of 98 alleles were found. The average number of alleles per locus was 8.15. Within breeds, the mean number of alleles ranged from 7.7 in Bardhoka to 8.8 in Shkodrane. Mean expected heterozygosity (He) ranged from 0.755 in Bardhoka to 0.781 in Shkodrane. Allelic richness varied from 7.63 in Bardhoka to 8.71 in Shkodrane. Inbreeding for all population is rather high FIS =0.061 ranging from -0.011 in Ruda to 0.087 in Shkodrane. The genetic differentiation, mean FsT was 0.106 and demonstrated that 89.4% of total genetic variation is due to genetic differentiation within each population. In the neighbor –joining tree based on DS distance, Bardhoka and Shkodrane were grouped together, than Ruda and finally Recka.
The Principal Component Analysis (PCA) analysis supported the conclusion suggested by dendrogram based on DS distance
Key words: genetic distance, genetic diversity, indigenous breeds, microsatellite, sheep
Sheep is one of the most important livestock species in Albania. There are four indigenous sheep breeds: Bardhoka, Ruda, Shkodrane and Recka. Bardhoka breed is classified in long tail group. Its origin is from north/northeast of Albania, as well as in Western part of Kosovo. Ruda is triple production breed with half-fine wool. This breed is part of Cigaya family, regarding the wool quality. It is originated and widespread in north-eastern part of Albania. The Shkodrane sheep breed belongs to the long tail group. The area of origin is north Albania. Recka breed is part of Cakel group. It is located all over the Albanian territory, especially on mountain areas of Central and South Albania. There were applied crossing between Bardhoka and Shkodrane, known as Baca, in order to improve the milk production, but no crossbreeding with Recka. Also, more than 30 years ago are introduced Cigaya rams for breeding Recka, in order to improve this breed. Selection programs in sheep are less advanced and systematic then in other livestock. Lack of parentage control and breed purity has facilitated a continuous gene flow. Because artificial insemination with frozen sperm is rarely practiced, this gene flow is limited by distance and geography. On the other hand, local management may have led to genetic isolation, which reduces the effective population size.
Polymorphic DNA markers are very useful in assessment of genetic diversity within and between breeds. Microsatellite is widely used as genetic markers for the analysis of genetic variability within and between breeds due to their high number, distribution throughout the genome and the efficacy of genotyping.
Recently, they are used for the study of diversity in European sheep breeds (Álvarez et al 2004, Arranz et al 1998, Diez-Tascón et al 2000, Farid et al 2000, Pariset et al 2003, Rendo et al 2004, Stahlberger-Saitbekova et al 2001, Tapio et al 2005, Tapio et al 2003, Peter et al 2007, Cinkulov et al 2008).
For long time the genetic diversity of indigenous sheep breed has been evaluated on the basis of biochemical methods. In this study that diversity has been characterized four local sheep breeds using 6 microsatellite markers. The objective was to characterize the relationships among the breeds, by estimating genetic distances from microsatellite markers.
Blood samples of 124 unrelated animals representing 4 different Albanian sheep breeds were analyzed. For each breed were sampled 31 individuals. There were selected maximum three unrelated individuals (two females and one male) per flock based on the information provided by the farmer. Sampling was carried out from an average of 11 flocks per breed, in the traditional rearing area of each breed. There were selected only the farms with pure breed animals. The breeds are marginally farmed and indigenous: Bardhoka (Bar), Ruda (Rud), Shkodrane (Sko) and Recka (Rck) (Figure 1, Table 1). A total of 6 microsatellite markers were analyzed: BM8125, MAF65, OarCP34, OarFCB304, OarHH47, and OarVH72.
|
|
Table 1. Summary of sheep breeds sampled |
||||
Breed |
Acronym |
Location |
Utility |
DNA sample |
Bardhoka |
Bar |
North-West Albania |
Milk* |
31 |
Ruda |
Rud |
West-north west Albania |
Meat-Wool* |
31 |
Shkodrane |
Sko |
North-East Albania |
Milk* |
31 |
Recka |
Rck |
Southern part of Albania |
Milk* |
31 |
* This is main utility of these breeds that are triple purpose breeds |
Allele frequencies and tests of genotype frequencies for deviation from Hardy-Weinberg equilibrium (HWE) were carried out using exact tests of the GENEPOP V.1.2 package (Raymond and Rousset 2001), performing a probability test using Markov chain Monte Carlo simulation (dememorization 10,000, batches 1000, iterations per batch 5,000). Significant levels were calculated per locus, per population, and over all loci and populations combined. Genetix version 4.05.2 (http://www.genetix.univ-montp2.fr/genetix/genetix.htm) was used for measuring genetic diversity within population, by the calculation of observed heterozygosity (HO) and mean unbiased estimates of gene diversity (He), (Nei 1978), mean number of alleles (MNA) per locus, the number of private alleles (PA, alleles found in only one breed). The program FSTAT (http://www2.unil.ch/popgen/softwares/fstat.htm) was used for the calculation of corrected allele diversity (allelic richness). F-statistics is widely used to characterize population genetic structure. These statistics allow the partition of genetic diversity/heterozygosity within and among populations. FIS measures the heterozygote deficit within populations, FST among populations, and FIT the global deficit of heterozygotes. Therefore, F-statistic parameters, (FIT, FIS, and FST) were estimated in the form of F, q and f respective estimators of these parameters proposed by Weir and Cockerham (1984). These were computed using FSTAT program (Goudet 1995). The significance was tested by 1000 permutations (Goudet 1995).
There were calculated two genetic distances from allele frequencies using GENETIX software, Nei genetic distance (DA) that is useful in case of closely related populations, where the genetic drift is the primary factor of genetic differentiation. It was calculated also standard genetic distance of Nei (1972) (DS). DS distance is used for the construction of Neighbor-Joining consensus tree (Saitou and Nei 1987) with Phylip package (http//evolution.genetics.washington.edu/phylip.html). Bootstrap (1000 replicates) resampling was performed to test the robustness of the dendrogram topology.
PCA using microsatellite data was performed using XLSTAT program (Angresti 1990; Saporta 1991)
The allele and genotype frequencies of six microsatellite loci were determined in 4 Albanian sheep breeds. In table 2 are summarized the estimates of various measures of diversity at the breed level. All the markers were polymorphic.
Table 2. Measures of genetic variability – Observed allele number (NA), allelic richness AR), observed heterozygosity (HO), Gene diversity (HE), FIS across the four breeds |
||||||
|
NA |
AR |
FIS |
HO |
HE |
|
BM8125 |
Bardhoka |
6 |
5.966 |
-0.117 |
0.733 |
0.656 |
|
Ruda |
7 |
6.932 |
-0.057 |
0.806 |
0.763 |
|
Shkodrane |
6 |
6 |
0.028 |
0.677 |
0.697 |
|
Recka |
7 |
6.899 |
-0.009 |
0.813 |
0.805 |
MAF65 |
Bardhoka |
8 |
7.871 |
-0.107 |
0.871 |
0.787 |
|
Ruda |
7 |
6.806 |
0.056 |
0.677 |
0.718 |
|
Shkodrane |
7 |
6.99 |
0.119 |
0.677 |
0.769 |
|
Recka |
10 |
9.71 |
-0.083 |
0.844 |
0.779 |
OARCP3 |
Bardhoka |
6 |
5.967 |
0.016 |
0.800 |
0.813 |
|
Ruda |
6 |
5.935 |
-0.148 |
0.903 |
0.787 |
|
Shkodrane |
7 |
6.933 |
0.083 |
0.733 |
0.8 |
|
Recka |
7 |
6.812 |
0.105 |
0.688 |
0.768 |
OARFCB |
Bardhoka |
10 |
10 |
0.322 |
0.483 |
0.712 |
|
Ruda |
9 |
8.742 |
-0.09 |
0.806 |
0.74 |
|
Shkodrane |
13 |
12.668 |
0.06 |
0.742 |
0.789 |
|
Recka |
8 |
7.884 |
0.127 |
0.625 |
0.716 |
OarHH4 |
Bardhoka |
8 |
8 |
0.191 |
0.655 |
0.81 |
|
Ruda |
9 |
8.994 |
0.154 |
0.710 |
0.839 |
|
Shkodrane |
12 |
11.674 |
0.144 |
0.645 |
0.754 |
|
Recka |
10 |
9.812 |
-0.024 |
0.875 |
0.855 |
OarVH7 |
Bardhoka |
8 |
8 |
0.054 |
0.793 |
0.839 |
|
Ruda |
9 |
8.997 |
0.007 |
0.871 |
0.877 |
|
Shkodrane |
8 |
7.999 |
0.082 |
0.800 |
0.871 |
|
Recka |
8 |
7.812 |
-0.044 |
0.875 |
0.838 |
N: sample size; NA: Number of alleles; AR: allelic richness; FIS: within population inbreeding estimates; HO: observed heterozygosity; HE: Expected heterozygosity. |
The total number of alleles and allele size for each locus are presented in Table 3.
Table 3. Alleles
number, allelic range, observed heterozygosity (HO),
expected heterozygosity (HE), PIC values, fixation
indices (FIS, FIT, FST) for each
marker in |
|||||||||||
Loci |
NA |
Allelic range (bp) |
HO |
HE |
HS |
HT |
PIC |
FIS |
FIT |
FST |
GST |
BM8125 |
9 |
108-126 |
0.757 |
0.719 |
0.73 |
0.769 |
0.743 |
-0.036 |
0.033 |
0.067*** |
0.066 |
MAF65 |
19 |
117-138 |
0.767 |
0.751 |
0.763 |
0.848 |
0.946 |
-0.006 |
0.125*** |
0.13*** |
0.129 |
OARCP3 |
14 |
111-135 |
0.781 |
0.779 |
0.792 |
0.876 |
0.949 |
0.014 |
0.139*** |
0.127*** |
0.124 |
OARFCB |
24 |
149-189 |
0.664 |
0.726 |
0.739 |
0.846 |
0.946 |
0.098* |
0.246*** |
0.164*** |
0.161 |
OarHH4 |
23 |
121-155 |
0.721 |
0.800 |
0.815 |
0.89 |
0.956 |
0.11** |
0.212*** |
0.112*** |
0.111 |
OarVH7 |
9 |
125-141 |
0.835 |
0.842 |
0.856 |
0.877 |
0.951 |
0.024 |
0.054 |
0.031*** |
0.03 |
Overall |
16 |
|
0.754 |
0.769 |
0.782 |
0.851 |
0.915 |
0.034* |
0.135*** |
0.106*** |
0.104 |
*p<0.05; **p<0.01; ***p<0.001 |
A total of 98 alleles were detected over all loci in 124 individuals. The average number of alleles varied from 9 (BM8125 and OARVH7) to 24 (OARFCB). All markers showed more than 5 alleles. Observed heterozygosity (HO) per locus ranged from 0.664 (OARFCB) to 0.835 (OARVH7). There were found some breed specific alleles. Frequencies of private alleles for Bardhoka and Shkodrane were <5% and OARFCB304179 for Ruda had the frequency of 8.1%. In Recka the frequency of breed specific alleles was higher than 5%. In total 1 from 24 locus-population comparisons revealed significant departures (p>0.05) from Hardy-Weinberg proportions.
The polymorphism information content (PIC) values for all loci ranged from 0.743 (BM8125) to 0.956 (OarHH4). The PIC values are very high, indicating that these loci are suitable for genetic studies of sheep breeds.
At individual loci, FIS estimates over all sheep breeds ranged between -0.036 (locus BM8125) and 0.11 (locus OarHH4), (Table 4).
Table 4. Measures of genetic variability in autochthones sheep breed |
|||||||||
Breed |
Acronym |
N |
Ne |
AR |
He |
Ho |
MNA |
ENA |
FIS |
Bardhoka |
Bar |
31 |
4.34 |
7.63 |
0.755 |
0.723 |
7.7 |
3.9 |
0.061 |
Ruda |
Rud |
31 |
4.76 |
7.73 |
0.775 |
0.796 |
7.8 |
4.1 |
-0.011 |
Shkodrane |
Sko |
31 |
4.53 |
8.71 |
0.766 |
0.713 |
8.8 |
3.9 |
0.087 |
Recka |
Rck |
31 |
4.77 |
8.15 |
0.781 |
0.787 |
8.3 |
4 |
0.009** |
All breeds |
124 |
4.6 |
11.90 |
0.769 |
0.755 |
8.15 |
3.8 |
0.061 |
On average, breeds had a 3.4% significant deficit of heterozygotes (p<0.001), whereas the total population had a 13.5% significant deficit of heterozygotes (p<0.001). To the observed deficit contribute only the markers Oarfcb and OarHH4.
The different estimates of genetic differentiation (FST and GST) with FIS and FIT are shown in table 2. FST values of genetic differentiation and GST values of breed differentiation were similar. The average genetic differentiation between all breeds was 10.6% significantly different from zero (p<0.001). FST values indicated that about 11% (p<0.001) of the total genetic variation was explained by a population difference and 89% correspond to the differences among individuals. The overall GST value was 0.104. That indicates that 10.4% of the total genetic diversity was observed among populations while 99.6% was within population.
The mean genetic differentiation within population (HS) was 0.783 and in total populations (HT) was 0.851. The lowest value showed BM8125 (0.73 and 0.769 respectively) and highest at OarVH7 (0.856 and 0.877 respectively).
Allelic richness (AR) ranged from 7.63 in Bardhoka to 8.71 in Shkodrane. The mean expected heterozygosity (HE) ranged from 0.755 in Bardhoka to 0.781 in Recka. Observed heterozygosity (HO) varied from 0.713 in Shkodrane to 0.796 in Ruda. Mean values of HE and HO, overall loci and breeds were 0.769 and 0.755 respectively. The mean number of alleles per locus (MNA) varied from 7.7 in Bardhoka to 8.8 in Shkodrane. Positive values of FIS per locus and per breed indicate inbreeding in all breeds, except of Ruda. FIS values ranged from -0.011 in Ruda to 0.087 in Shkodrane.
Nei’ genetic distance DA ranged from 0.067 between Ruda and Shkodrane to 1.807 between Shkodrane and Recka (Table 5).
Table 5. DA (above) and DS (below) genetic distance |
||||
|
Bardhoka |
Ruda |
Shkodrane |
Recka |
Bardhoka |
|
0.118 |
0.111 |
1.665 |
Ruda |
0.024 |
|
0.067 |
1.714 |
Shkodrane |
0.016 |
0.018 |
|
1.807 |
Recka |
0.199 |
0.195 |
0.204 |
|
The highest DS value was observed between Shkodrane and Recka (0.204) and the lowest was between Bardhoka and Shkodrane. A NJ topology tree based on DS genetic distance relating 4 sheep populations studied is presented in figure 2.
|
|
The number at the nodes is values for 1000 bootstrap resampling of the typed loci. The bootstrap values obtained at the nodes are high (>50%) and suggest that the robustness of the tree is high. The phylogenetic tree showed the grouping of Bardhoka and Shkodrane. The most separated, appeared Recka.
In the scattergram (figure 3), Bardhoka and Shkodrane exhibited a close relationship, forming the same group, thus supporting the conclusions suggested by the dendrograme based on Nei’s distance.
|
|
From both analyses Bardhoka and Shkodrane breeds formed the closest group.
This study investigated the allele frequencies and genotype distributions of 6 microsatellite loci in 124 individuals of local sheep breeds. The genetic analysis of 4 Albanian local sheep breeds with 6 microsatellite markers showed high gene diversity. The microsatellite loci for genetic diversity studies should have more than four alleles. The microsatellite markers showed a high genetic polymorphism. The total number of alleles per locus in the present study ranged from 9 to 24. The high number of alleles for each locus, as well as the high PIC values, suggested that all the markers used were appropriate to analyze diversity in Albanian local sheep breeds. The mean number of alleles per breed is a good indicator of genetic variation within populations. Gene diversity for each breed ranged from 0.75 in Bardhoka to 0.78 in Recka. According to Takezaki and Nei (1996), a marker can be useful in measuring genetic variation if the average heterozygosity ranges from 0.3 to 0.8 in the population. This confirmed that these markers were appropriate for measuring genetic variation. These breeds are characterized by a local distribution, small population size, living in extensive condition, grazing in the poor pastures, linked to traditional farming system. This may be the explanation for the high degree of genetic variability.
The high value of heterozygosity may be result of application of empiric selection and different breeding objectives, as well as application of uncontrolled natural mating. The heterozygosity values observed in Albanian sheep breeds are higher to other sheep breeds: 0.68 in French Mutton Merino by Diez Tascon et al (2000); 0.70 in Muzzafaranagri sheep by Arora and Bhatia (2004); 0.71 in Baltic sheep breeds by Tapio et al (2005); 0.67 in Nali and Chokla sheep breeds by Sodhi et al (2006); 0.69 in Magra sheep by Arora and Bhatia (2006).
The genetic differentiation (FST) among the local goat breeds in this study is 10.6%. The FST values obtained in this study suggested a moderate genetic differentiation between Albanian sheep breeds. The values of genetic differentiation are higher than those reported by other authors. The genetic differentiation were 5.7% between 57 European and Middle Eastern sheep breeds (Peter et al 2007), 5.2% between West Balkan pramenka sheep types (Cinkulov et al 2008), 8% between finish sheep breed, (Tapio et al 2003), 8.8% between Baltic sheep breeds (Tapio et al 2005), or 8.5% (Forbes et al 1995). In all these cases the number of markers was much higher than in our case. It is clear that 89.4% of genetic variation corresponds to differences among individuals and 10.6% is result of differences among breeds.
In the present study FIS=0.035, FIT=0.137 and FST=0.106. The values of FIS are positive except of Ruda. The cause for the loss of heterosigotes may be the mating between the relatives. Only Bardhoka breed showed 2 of 6 loci with a significant deficit of heterozygotes, therefore the mating between relatives is not likely. The excess of homozygotes could be due to locus under selection, or presence of “null alleles” that may lead to a falce observation of homozygotes. A very feasible explanation for the defficit of heterozygotes is the presence of Wahlund effect (population substructure within the breed). Since sampling is carried out in 11 flock per breed, therefore it is possible the presence of genetic structure.
The genetic differentiation (10.6%) may be as result of inbreeding and of random genetic drift. The moderate values of FST, is most likely attributed to genetic drift between the investigated breed. There is no herd book since 1990 and maybe there is some admixture of breeds. In Albania are applied crossing between Bardhoka and Shkodrane, known as Baca, in order to improve the milk production. This may be one explanation for the smallest distance between Bardhoka and Shkodrane and their grouping, revealed by phylogenetic analysis and PCA. Recka showed always the largest genetic distance with other breeds. There is no crossbreeding with Recka. Also, more than 30 years ago are introduced Cigaya rams for breeding Recka, in order to improve this breed. These can explain why Recka behave like a separate breed.
Álvarez I, Royo, L J, Fernández I, Gutiérrez J P, Gómez E and Goyache F 2004 Genetic relationships and admixture among sheep breeds from Northern Spain assessed using microsatellites. Journal of Animal Science 82: 2246-2252 http://jas.fass.org/cgi/reprint/82/8/2246
Angresti A 1990 Categorical data analysis. New York Willey Interscience.
Arora R and Bhatia S 2004 Genetic structure of Muzzafarnagri sheep based on microsatellite analysis. Small Ruminant Research 54: 227-230
Arora R and Bhatia S 2006 Genetic diversity of Magra sheep from India using microsatellite analysis. Asian-Australian Journal of Animal Science 19: 938-942
Arranz J J, Bayon Y and San Primitivo F 1998 Genetic relationships among Spanish sheep using microsatellites. Animal Genetics 29: 435-440
Cinkulov M, Popovski Z, Porcu K, Tanaskovska B, Hodzic A, Bytyqi H, Mehmeti H, Margeta V, Djedovic R, Hoda A, Trailovic R, Brka M, Markovic B, Vazic B, Vegara M, Oslaker I and Kantanen J 2008 Genetic diversity and structure of the West Balkan pramenka sheep types as revealed by microsatellite and mitochondrial DNA anaylsis, Journal of Animal Breeding and Genetics 125: 417-426
Diez-Tascon C, Littlejohn R P, Almeida P A and Crawford A M 2000 Genetic variation within the Merino sheep breeds analysis of closely related populations using microsatellites. Animal Genetics 31: 243–251
Farid A, O'reilly E, Dollard C and Kelsey, C R 2000 Genetic analysis of ten sheep breeds using microsatellite markers. Canadian Journal of Animal Science 80: 9-17
Forbes S H, Hogg J T, Buchanan F C, Crawford A M and Allendorf F W 1995 Microsatellite evolution in congeneric mammals domestic and bighorn sheep. Molecular Biology and Ecology 12: 1106-1113
Goudet J 1995 FSTAT V 2.9.3 a computer programme to calculate F-statistics http//www.unil.ch/izea/software/fstat.html. Journal of Heredity 8: 485-486
Nei M 1972 Genetic distance between populations. AAmerican Naturalist. 106, 283-292.
Nei M 1978 Estimation of average heterozygosity and genetic distance from a small number of individuals. Genetics 89, 583–90.
Pariset L, Savarese M C, Cappucio I and Valentini A 2003 Use of microsatellites for genetic variation and inbreeding analysis in Sarda sheep flocks of central Italy. Journal of Animal Breedind and Genetics 120: 425-432
Peter C, Bruford M, Perez T, Dalamitra S, Hewitt G, Erhardt G and The Econogene Consortium 2007 Genetic diversity and subdivision of 57 European and Middle Eastern sheep breeds. Animal Genetics 38: 37 – 44
Raymond M and Rousset F 2001 GENEPOP a population genetic software for exact test and ecunemism. http://jhered.oxfordjournals.org/cgi/reprint/86/3/248
Rendo F,
Iriondo M, Jugo B M, Mazón L I, Aguirre A, Vicario A and Estonba A 2004
Tracking diversity and differentiation in six sheep breeds from the North
Iberian Peninsula through DNA variation. Small Ruminant
Research 52 : 196-202
Saitou N and Nei M 1987 The neighbor-joining method: a new method for reconstructing phylogenetic trees. Mol. Biol. Ecol. 4, 406-425.
Saporta G 1991 Probability analyse de donnis et statistiques. Paris. Technip.
Sodhi M, Mukesh M and Bhatia S 2006 Characterizing Nali and Chokla sheep differentiation with microsatellite markers, Small Ruminant Research 65: 185-192
Stahlberger-Saitbekova N, Schläpfer J, Dolf G and Gaillard C 2001 Genetic relationships in Swiss sheep breeds based on microsatellite analysis. Journal of Animal Breeding and Genetics 118: 379-387
Takezaki N and Nei M 1996 Genetic distances and reconstruction of phylogenetic trees from microsatellite DNA. Genetics 144: 389–3999
Tapio M, Miceikiene I, Vilkki J and Kantanen J 2003 Comparison of microsatellite and blood protein diversity in sheep inconsistencies in fragmented breeds. Molecular Ecology 12: 2045-2056
Tapio I, Tapio M, Grislis Z, Holm L E, Jeppsson S, Kantanen J, Miceikiene I, Olsaker I, Viinalass H and Eythorsdottir E 2005 Unfolding of population structure in Baltic sheep breeds using microsatellite analysis. Heredity 94: 448-456
Weir B S and Cockerham C 1984 Estimating F-statistics for the analysis of population structure. Evolution 38: 1358-1369
Received 5 March 2009; Accepted 15 May 2009; Published 1 June 2009