Livestock Research for Rural Development 21 (6) 2009 Guide for preparation of papers LRRD News

Citation of this paper

Genetic diversity and distances of Albanian local sheep breeds using microsatellite markers

A Hoda, P Dobi and G Hyka

Department of Animal Production, Faculty of Agriculture and Environment, Agricultural University of Tirana, (AUT), Albania
hodanila@yahoo.com

Abstract

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


Introduction

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.

 

Materials and methods 

Sample collection and 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.



Figure 1.  The distribution of sampling locations of four local sheep breeds in Albania



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


Statistical analyses

 

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)

           

Results 

Microsatellite loci and genetic variability

 

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
4 autochthones sheep breeds

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).

 

Genetic variation within and among breeds

 

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.

 

Genetic distances

 

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.



Figure 2.  Consensus “unrooted” tree, based on DS genetic distance, through 6 microsatellite marker


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.



Figure 3.  PCA based on the allele frequencies


From both analyses Bardhoka and Shkodrane breeds formed the closest group.

 

Discussion 

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.

 

Literature 

Á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

Go to top