Investigation of genetic relationship among populations has been
traditionally based on the analysis of allele frequencies at different loci.
The prime objective of this research was to measure the genetic polymorphism
of five microsatellite markers (McMA2, BM6444, McMA26, HSC, and OarHH35) and
study genetic diversity of 14 sheep types in Iran. Genomic DNA was extracted
from blood samples of 565 individuals using an optimized salting-out DNA
extraction procedure. The polymerase chain reaction (PCR) was successfully
performed with the specific primers. Some locus–population combinations were
not at Hardy–Weinberg equilibrium (
Small ruminants, especially native breed types, play an important role in the livelihoods of a considerable part of human population in the tropics from socioeconomic aspects. Therefore, an integrated attempt in terms of management and genetic improvement to enhance production is of crucial importance (Mohammadabadi and Sattayimokhtari, 2013). Economical and biological efficiency of sheep production enterprises generally improves by increasing productivity and reproductive performance of ewes (Mohammadabadi and Sattayimokhtari, 2013). There are more than 50 million sheep in Iran, of 27 breeds and ecotypes (Khodabakhshzadeh et al., 2016) that have not defined well as distinct breeds. However, they are considered as geographically defined populations. The need to maintain and improve local genetic resources has been recognized as a priority at the world level. Biodiversity studies depicting a deep picture of the genetic variability of the available sheep breeds provide favorable opportunities for both genetic conservation programs and enhancing production efficiency by means of controlled and well-designed crossbreeding systems exploiting breed diversities, heterosis and breed complementarity (Esmailizadeh et al., 2012).
Genetic diversity in indigenous breeds is a major concern considering the necessity of preserving what may be a precious and irreplaceable richness with regard to new productive demands (Khodabakhshzadeh et al., 2016). Conservation should be based on a deep knowledge of the genetic resources of the specific breed (Zamani et al., 2015). Therefore, it is important to try to characterize genetically indigenous breeds. Genes affecting polygenic traits and characterizing milk or meat productions are difficult to identify (Soufy et al., 2009; Shojaei et al., 2011). The maintenance of genetic diversity in livestock species requires the adequate implementation of conservation priorities and sustainable management programs, which should be based on comprehensive information regarding the structure of the populations, including sources of genetic variability among and within breeds. Genetic diversity is an essential component for population survival, evolution, genetic improvement and adaptation to changing environmental conditions (Kumar et al., 2006). Molecular methods based on molecular markers, such as random amplification of polymorphic DNA (RAPD), restriction fragment length polymorphism (RFLP) and microsatellites, are useful tools to study the genetic variations. Short tandem repeats known as microsatellites are widely used as molecular markers of choice for genetic studies. Advantages of microsatellites are a high degree of polymorphism due to existence of several alleles at each locus, their large number, distribution throughout the genome, a high level of polymorphism, neutrality with respect to selection, codominant inheritance and easy automation of analytical procedures (Canon et al., 2006; Mohammadifar et al., 2009; Mohammadabadi et al., 2010). Several studies have investigated the genetic diversity in sheep using microsatellites (Buchanan and Thue, 1998; Esmaeilkhanian and Banabazi, 2006; Bhatia and Arora, 2007; Nanekarani et al., 2010; Sun et al., 2010; Jakaria et al., 2012; Musthafa Muneeb et al., 2012; Hepsibha et al., 2013; Crispim et al., 2014), but a study of all Iranian sheep together has not been performed until now. Hence, the aim of the present study was to evaluate the genetic diversity within and between 14 sheep types in Iran (Kermani, Pakistani, Lori, Arabi, Dalagh, Baluchi, Iran-Black, Gharegol, Arman, Lori–Bakhtiari, Kermani–Pakistani, Kermani–Romanov, Lori–Bakhtiari–Romanov and Lori–Bakhtiari–Pakistani) using five microsatellite markers and to measure the distance among these breeds.
In this study, 565 blood samples were collected from different individuals of
14 sheep types in Iran (Kermani (KER),
Characteristics of selected microsatellite in present study.
In this study, five microsatellite markers across the sheep genome were
used. Gradient PCR was used to optimize the annealing temperature for each
marker. The PCR products were tested in agarose gel (0.8 %) to estimate the
best annealing temperature for each primer. The studied microsatellite
markers, their primer sequences, detected annealing temperature and their
allele size ranges are shown in Table 1. The selected microsatellites were
amplified with PCR using genomic DNA extracted from individual animals (in
total 25
Genotypes were assigned for each animal based on allele size data. Frequencies and number of alleles for each locus, observed and expected heterozygosity were estimated using FSTAT (version 2.9.3.2) (Goudet, 2002). The polymorphic information content (PIC) value was calculated according to Buchanan and Thue (1998). Nei's standard genetic distances (DS) among populations were computed by POPGENE (Yeh et al., 1999). This software was also used to construct the dendrogram of unweighted pair group with arithmetic mean (UPGMA).
The number of actual alleles (Na) and effective alleles (Ne) for different combinations of locus–population and for each population.
The PCR reactions were successfully performed with all primers. All the
microsatellite loci were found to be highly polymorphic. In total, 65
alleles were detected; the HSC marker in PAK sheep breed and overall showed
the highest number of alleles per locus (14 and 15 alleles respectively)
while the OarHH35 marker in KER-PAK Sheep tape showed the lowest number of
alleles (5 alleles) (Table 2) with a mean of
Polymorphic information content (PIC), Shannon information index
(I) and expected heterozygosity (
PIC and the Shannon information index are another measure of genetic variability indicating the informativeness of the assessed loci. PIC values ranged from 0.84 (HSC marker) to 0.92 (BM6444 marker) with a mean PIC value of 0.87 (Table 3), indicating that all loci were highly polymorphic.
Mean polymorphic information content (PIC), mean Shannon
information index (I), mean expected heterozygosity (
The mean expected heterozygosity over all breeds in the present study was
0.85 (Table 3), so it does not likely encounter problems that result from
inbreeding depression. Mean estimates of expected heterozygosity overall
loci and types were 0.86 (Table 4). ARB had the highest expected
heterozygosity for all the loci (
When the Hardy–Weinberg testing was performed for the loci, deviations from
the Hardy–Weinberg equilibrium were found to be significant (
UPGMA phylogenetic tree based on Nei genetic distance. Location of these 14 studied breeds is shown in Fig. 2. Subspecies are indicated as Kermani (KER), Pakistani (PAK), Lori (LOR), Arabi (ARB), Dalagh (DAL), Baluchi (BAL), Iran-Black (IRB), Gharegol (GHA), Arman (ARM), Lori–Bakhtiari (LRB), Kermani–Pakistani (KER-PAK), Kermani–Romanov (KER-ROM), Lori–Bakhtiari–Romanov (LRB-ROM) and Lori–Bakhtiari–Pakistani (LRB-PAK).
Estimators of F statistics at each locus across the 14 sheep types.
Genetic distance matrices based on
Population differentiation examined by fixation indices such as
Genetic variation is a basic requirement for animal breeding, whereas a high genetic variation is needed for genetic improvement of domestic animals (Askari et al., 2011). The number of alleles at different marker loci serves as a measure of the genetic variability having direct impact on differentiation of breeds within a species. Since 80 % of the markers exhibited four or more alleles, the microsatellite loci screened in this study were appropriate in expressing the molecular characteristics and/or genetic variation in the population. The effective number of alleles at each locus provides information on predominant alleles. Since allelic diversity in the studied sheep breed populations was high, it can be concluded that, in these animals, genetic diversity is sufficiently high and they have a good gene pool for breeding programs. Because of the limitations in studied population number and microsatellite loci, further research with more population and loci numbers needs to be done to be able correctly evaluate the genetic relationship among the populations. However, a higher number of alleles for each locus showed that all the markers used were appropriate to analyze diversity in studied breeds. The level of variation depicted by number of alleles at each locus was similar to earlier reports on sheep breeds (Arora and Bhatia, 2006; Jakaria et al., 2012; Hepsibha et al., 2013; Crispim et al., 2014).
Among the types, the mean of PIC ranged from 0.72 in LRB-PAK to 0.88 in ARB, and Shannon information index ranged from 1.58 in KER-PAK to 2.27 in ARB. . However, PIC in the range of 0.69 to 0.92 was reported for this locus in some Iranian sheep breeds in a previous study (Esmaeilkhanian and Banabazi, 2006; Banabazi et al., 2007; Nanekarani et al., 2010). Our estimate of Shannon information index was in agreement with that reported by Nanekarani et al. (2010) on three types of Iranian sheep populations and Musthafa Muneeb et al. (2012) on Najdi sheep populations. This index estimate was higher than that of Mehraban sheep (Zamani et al., 2011) and Indian Bellary sheep (Kumar et al., 2006).
A high value of average expected heterozygosity within a breed could be attributed to the large number of alleles detected in the tested loci (Kalinwski, 2002). The high heterozygosity is attributed to the wider distribution of population and to presumably the larger flock size. It also proves that no controlled breeding is followed as males and females are allowed to graze as a large flock. Heterozygosity values were in concordance with those of Najdi sheep populations (Musthafa Muneeb et al., 2012), Indian sheep breeds (Bhatia and Arora, 2007) and five Iranian indigenous sheep populations (Esmaeilkhanian and Banabazi, 2006). The expected heterozygosity in this study was higher than that of some of loci (OarHH35) used in Hu sheep in China (Sun et al., 2010). Thus, differences of heterozygosity values obtained could be due to type of markers used, the sequence of simple sequence repeat (SSR) markers, and sampling and sample size.
In total, 58.6 % of the loci in some breeds deviated from HWE (
Location of the 14 study sites in Iran. The putative
subspecies are indicated as Kermani (KER), Pakistani (PAK), Lori (LOR),
Arabi (ARB), Dalagh (DAL), Baluchi (BAL), Iran-Black (IRB), Gharegol (GHA),
Arman (ARM), Lori–Bakhtiari (LRB), Kermani–Pakistani (KER-PAK),
Kermani–Romanov (KER-ROM), Lori–Bakhtiari–Romanov (LRB-ROM) and
Lori–Bakhtiari–Pakistani (LRB-PAK). Source of original map:
OCHA/ReliefWeb
(
The within-breed inbreeding estimate (
There are three separate clusters on the dendrogram. One includes 12 breeds
and another consists of Arabi and then Pakistani in a separate branch. Since
The breeds included in this study are of great importance to the sheep holders in Iran. This study presents an investigation of variability at the DNA level within and between some Iranian sheep breeds. The results indicated that all studied breeds exhibited considerable genetic variation, based on their high mean number of alleles and gene diversity. According to the selective standard of microsatellite loci, microsatellite loci ought to have at least four alleles per locus to be considered useful for the evaluation of genetic diversity. Based on this criterion, the five microsatellite loci used in the present study can be considered useful for the evaluation of genetic diversity within and among populations and for the selection of breeding animals from divergent groups maximizing genetic variation and consequently fitness. A future direction to our study can be studying all of the Iranian indigenous sheep breeds to better evaluate the level of inbreeding and establish appropriate conservation strategies with the aim to avoid losses of genetic diversity.
The original data of the paper are available upon request from the corresponding author.
MM and AEK designed this project, MTVE collected blood from animals and performed experiments in the laboratory, and MM and MTVE wrote the paper.
The authors declare that they have no conflict of interest.
We would like to thank the sheep farmers who provided sheep blood samples for this study. Edited by: Steffen Maak Reviewed by: Sonia Zakizadeh and Masood Asadi Fozi