The aim of this study was to investigate the genetic diversity and structure in three Iranian native goat breeds (Markhoz, Mahabadi and Lori) and the Beetal imported breed using inter-simple sequence repeat (ISSR) markers and also to investigate ISSR markers' potential in order to genetically separate single (S) and twin-birth (T) subpopulations. Blood samples were collected from 210 animals for this purpose. In total, 16 primers were used, and finally 5 primers were selected based on the number of clear bands and the level of polymorphisms. The result of this study showed that 76 of 86 observed fragments were polymorphic. Genetic diversity for each breed ranged from 0.23 in the Beetal breed to 0.26 in the Markhoz breed; this represents a relatively similar genetic diversity in these breeds. An unweighted pair group method with arithmetic mean (UPGMA) dendrogram based on the Nei's standard genetic distance between the breeds studied showed that three Iranian goat breeds (Mahabadi, Lori and Markhoz) were clustered closer together, while the Beetal breed formed a separate cluster. In the constructed dendrogram of the subpopulations, the S and T subpopulations of each breed were clustered together. The constructed dendrogram of the Beetal breed and the S and T subpopulations of all breeds studied showed a separate cluster for the Beetal breed as an imported breed and another cluster for the S and T subpopulations as Iranian native breeds. The current study showed that the ISSR markers studied had no potential to genetically separate S and T subpopulations. On the other hand, these ISSR markers can be used for the clustering of distinct populations.
List of primers used for ISSR amplification and their annealing temperature.
Goats are probably one of the earliest ruminants domesticated around 10 000 years ago (Zeder and Hesse, 2000). Today, goats are raised all over the world due to their tolerance to various environments, adaptability to nutrient-poor diets and manageable size. Goats are multipurpose animals that produce milk, meat and wool and contribute to the economy of farmers living in arid and semiarid regions, including south of Iran (Vahidi et al., 2014). The Iranian goat population is 21 million, which is 2.57 % of the world's goat population (861.90 million) and 4.08 % of the Asian goat population (514.4 million).
Genetic diversity is an essential component for population survival, evolution, genetic improvement and adaptation to changing environmental conditions (Kumar et al., 2009). The conservation of genetic variability has arisen in the last years because the existence of a large gene pool is important for the potential future breeding preservation and for sustainable animal production system development. Comprehensive knowledge of the existing genetic variability is the first step for the conservation and exploitation of domestic animal biodiversity.
Molecular tools permit us to characterize genetic resources at the DNA level. Inter-simple sequence repeat (ISSR) regions in the genome are currently useful markers for a wide range of molecular genetic studies such as establishing population structure, population differentiation and the reconstruction of phylogenetic relationships between populations, especially in plants (Wang et al., 2008). The ISSR markers are created using a single primer complementary to the microsatellite repeat (from 4 to 12 repeat units) and carrying the 3' sequence consisting of one or two arbitrary nucleotides (the so-called “anchor”). Such primers enable the amplification of the DNA fragments which are situated between two rather closely located inverted microsatellite sequences (Sulimova, 2004). The ISSR marker has a mild technical difficulty, good reproducibility and a reasonable cost, permitting its use for genetic studies of populations (Chen et al., 2005; Fagen and Nianhe, 2005)
The aim of this study was to investigate the genetic diversity and structure in three Iranian native goat breeds (Markhoz, Mahabadi and Lori) and the Beetal imported breed using ISSR markers and also to investigate ISSR markers' potential to separate single (S) and twin-birth (T) subpopulations genetically.
All protocols used in this study were approved by the Animal Care and Use
Committee of Razi University. Blood samples were collected from a total of
210 individuals belonging to three Iranian native goat breeds and the Beetal
imported breed: 60 Markhoz does (30 with an LS (litter size) of 1 and 30 with
an LS of 2) from Sanandaj, 60 Mahabadi does (30 with an LS of 1 and 30 with an LS of 2)
from Mahābād, 60 Lori does (30 with an LS of 1 and 30 with an LS of 2) from Lorestān
and 30 Beetal does (LS
DNA was extracted from whole blood using a DIAtom DNA Prep 100 kit
(GenFanavaran Co, Iran). First, 16 ISSR primers were used, but 5 primers
were selected due to the number of clear bands and polymorphic properties.
Annealing temperature and sequences for all primers are listed in Table 1.
Polymerase chain reaction was performed in a 15
Various measures of genetic variation at different loci across studied goat breeds.
Na: observed number of alleles; Ne: effective number of alleles; H: Nei's genetic diversity; I: Shannon's information index.
Various measures of genetic variability among different goat breeds.
Na: observed number of alleles; Ne: effective number of alleles; H: Nei's genetic diversity; I: Shannon's information index.
The measurement of diversity including gene diversity (H), observed number of alleles (Na), effective number of alleles (Ne) and Shannon's information index was done by the POPGEN 1.31 software (Yeh et al., 1999). The software was also employed to determine Nei's unbiased genetic distances (Nei, 1978). A dendrogram based on Nei's unbiased genetic distances, using the unweighted pair group method with arithmetic mean (UPGMA), was generated to show the genetic distances of the populations or subpopulations investigated in this study.
In the current study, the repeatability of the PCR products was initially
examined by the 16 primers on agarose gel for 21 samples. The test indicated
that the patterns of ISSRs are polymorphic and reproducible in five primers
including (ATG)
The genetic characterization of the studied goat breads is shown in Table 3. Among the breeds, the mean number of alleles ranged from 1.71 in Beetal to 1.80 in the Markhoz breed. Nei's gene diversity and Shannon's information index were relatively similar in the breeds studied (Table 3). It seemed that the breeds studied have relatively similar genetic diversity. Among the breeds, total heterozygosity (Ht), the average heterozygosity within populations (Hs) and diversity between the four populations (Dst) were 0.26, 0.25 and 0.01, respectively. The average genetic differentiation coefficient (Gst) value (Weir and Cockerham, 1984), another measure of genetic differentiation, over all loci was 0.056 indicating that 5.6 % of total genetic variation corresponded to differences between populations, whereas 94.4 % was explained by differences between individuals.
Nei's (1978) genetic identity (I, above diagonal) and genetic distance (D, below diagonal) among the four goat breeds.
Nei's (1978) genetic identity (I, above diagonal) and genetic distance (D, below diagonal) between single (S) and twin-birth (T) subpopulations of goat breeds.
The results of the estimated Nei's unbiased measures of genetic distances (Nei, 1978) for the breeds studied are presented in Table 4. According to this, the highest genetic distance between the breeds was 0.0264, which indicates a high genetic similarity between the breeds studied. A cluster analysis based on Nei's standard distance matrix (Nei, 1978) (Table 4; Fig. 1) showed two main groups: one consisting of the three Iranian goat breeds, which exhibited the Mahabadi and Lori breeds in one cluster, followed by Markhoz, and the other consisting of the Beetal breed. Nei's genetic distance between the S and T subpopulations of the breeds studied is shown in Table 5. Although the genetic distance between all subpopulations is low, this value between the S and T subpopulations of each breed is less than that of different breeds.
The phylogenetic dendrogram showing the genetic relationships between goat populations using the Nei (1978) standard genetic distance obtained from five ISSR markers.
The phylogenetic dendrogram showing the genetic relationships between the single (S) and twin-birth (T) subpopulations of the goat breeds studied using the Nei (1978) standard genetic distance obtained from five ISSR markers.
The phylogenetic dendrogram of the S and T subpopulations of the breeds studied based on Nei's genetic distance analysis (Fig. 2) placed the S and T subpopulations of each breed in one cluster, and, overall, they were divided into two groups: one includes S Mahabadi and T Mahabadi in one branch with S Lori and T Lori subpopulations. Another consists of Beetal and then the S Markhoz and T Markhoz subpopulations.
Nei's (1978) genetic identity and genetic distance between the Beetal breed and the S and T populations of the other breeds studied are shown in Table 6. According to this, the S and T populations are closer together and have a greater distance to the Beetal breed. The phylogenetic relationship between the Beetal breed and the S and T populations of the other breeds studied was constructed on the basis of Nei's standard distance matrix (Nei, 1978) (Fig. 3). As seen in Fig. 3, the S and T populations were grouped in one cluster and the Beetal breed in another one.
Genetic variation is a basic requirement for animal breeding, whereas high genetic variation is needed for the genetic improvement of domestic animals (Askari et al., 2011). The analysis of genetic diversity is a method to estimate the variation present in populations (Weir, 1990). Therefore, we used ISSR markers to assess genetic variation based on the mean allelic frequency of the 86 ISSR fragments obtained in the breeds studied. The total number of alleles per locus in the present study ranged from 14 to 19. The higher number of alleles for each locus showed that all the markers used were appropriate to analyze diversity in the breeds studied.
Nei's (1978) genetic identity (I, above diagonal) and genetic distance (D, below diagonal) between the Beetal breed and the single (S) and twin-birth (T) subpopulations of the other breeds studied.
The phylogenetic dendrogram showing the genetic relationships between the Beetal breed and the S and T populations of the other breeds studied using the Nei (1978) standard genetic distance obtained from five ISSR markers.
Genetic diversity for each breed ranged from 0.23 in Beetal to 0.26 in Markhoz. The measures of genetic variation in the present study indicated a moderate level (0.27) of Nei's genetic diversity in the four goat breeds studied. No significant differences in the genetic diversity were observed within populations, and they had a relatively similar genetic diversity. This result was in accordance with those of Italian goat populations (Ajmone-Marsan et al., 2001), Moxoto goat breed (Oliveira et al., 2005) and six Iranian indigenous goat populations (Mirhoseini and Badrbarin, 2013). The genetic diversity in this study was higher than that of some Iranian breeds (Askari et al., 2011; Zamani et al., 2011), and Najdi, Harri and Aradi Saudi goat breeds (Sabir et al., 2013); however, the diversity was noticeably lower than that of Korean and Chinese goats (Kim et al., 2002), twelve Chinese indigenous goat populations (Li et al., 2002), southern Italian goat populations (Iamartino et al., 2005), Egyptian and Italian goat breeds (Agha et al., 2008), and Indian goat breeds (Dixit et al., 2012) and was also slightly lower in comparison to the color types of the Markhoz goat populations (Moradii et al., 2014). Zamani et al. (2011), Askari et al. (2011) and Moradii et al. (2014) used two ISSR markers for their studies; the sequence of ISSR markers in their studies was also different from those in the present study. Sabir et al. (2013) used 30 goats in their study, which is much lower than the sample size in our study (210 goats); also, the sequence of some ISSR markers in their study was different from those in the present study. On the other hand, Kim et al. (2002), Li et al. (2002), Imartino et al. (2005), Agha et al. (2008) and Dixit et al. (2012) used microsatellite markers in their studies. Thus, the differences in genetic diversity values obtained could be due to the type of markers used, the sequence of ISSR markers, sampling and sample size.
The Shannon information index was 0.41, indicating that 41 % of effective alleles (a mean of 1.44 alleles per locus) were on average identified per primer. Our estimate of the Shannon information index was in accordance with what is reported by Moradii et al. (2014) on the color types of Markhoz goat populations and by Mirhoseini and Badbarin (2013) on six Iranian indigenous goat populations. This estimated index was higher than that for Mehraban sheep (Zamani et al., 2011), Rayini goats (Askari et al., 2011), and Najdi, Harri and Aradi Saudi goat breeds (Sabir et al., 2013). However, the Shannon information index in this study was lower than that of Kermani sheep (Esfandyarpour et al., 2008).
An overall mean of Gst
The phylogenetic dendrogram of the breeds (Fig. 1) showed that the Mahabadi and Lori breeds were grouped together, with the Markhoz breed deviating from them. This type of clustering was completely consistent with the fact that Markhoz is a Mohair-producing breed but Mahabadi and Lori are known as meat and milk-producing breeds. Furthermore, the Beetal breed formed a separate cluster from Iranian goats. It is noteworthy that this breed has been recently imported from Pakistan into Sīstān va Balūchestān province in Iran. Therefore, this observed difference can be due to geographical distance, which is exacerbated by the official borders between the two countries. Similar population clustering according to geographic location was previously observed in 12 Chinese indigenous goat breeds (Li et al., 2008, 2002), Egyptian and Italian goat breeds (Agha et al., 2008), Indian goat breeds (Dixit et al., 2012), and 6 Iranian indigenous goat populations (Mirhoseini and Badbarin, 2013).
In the dendrogram constructed of the subpopulations of the breeds studied (Fig. 2), the S and T subpopulations per breed were grouped in the same cluster. In fact, the ISSR markers used in this study showed a low potential for separating the S and T subpopulations as different clusters. The ISSR markers amplify intergenic and non-coding regions of the genome, and these regions may not have any effect on the number of kids (twining trait). It seems that a study on major genes known for the twining trait can distinguish S subpopulations from T subpopulations. For example, Rasouli et al. (2016) reported that the twinning rate in Markhoz goats was significantly influenced by the IGF-I and IGFBP-3 genotypes; also, the association of two single nucleotide polymorphisms (SNPs) in the growth hormone (GH) gene with the litter size and superovulation response in goat breeds has been studied, and the results have shown that the two loci of the GH gene are highly associated with abundant prolificacy and a superovulation response in goat breeds (Zhang et al., 2011).
The phylogenetic tree constructed between the Beetal breed and the S and T subpopulations for all the breeds studied (Fig. 3) showed a separate cluster for the Beetal breed as an imported breed and another cluster for the S and T subpopulations as Iranian native breeds. This finding confirmed that geographic distance and different agroclimatic conditions and probably breeding objectives have an important effect on genetic diversity.
Moradii et al. (2014) showed that the phylogenetic dendrogram of color types using ISSR markers was in agreement with their phenotypic similarities, in which saddle brown and sandy brown color types were grouped in one cluster and black types in another one. However, the current study is the first report using ISSR marker analysis to investigate the possibility of different clustering for the twining trait in several goat populations (Agha et al., 2008; Dixit et al., 2012; Javanrouh Aliabad et al., 2004; Li et al., 2008; Mirhoseini and Badbarin, 2013).
However, the low number of ISSR markers (five markers) and the sample size (30 goats per subpopulation) in the current study may be other reasons for the low potential of the ISSR markers studied to show S and T subpopulations in different clusters. It seems that this would be clearer using many more ISSR markers in larger goat populations.
In conclusion, this study has shown the moderate genetic diversity in some native goat breeds. This research has also shown that the ISSR markers studied had no potential to genetically separate S and T subpopulations, while ISSR markers can be used to cluster distinct populations.
The original data of the paper will be available upon request to the corresponding author.
A. Abdolmohammadi and S. Foroutanifar designed this project, L. Simaei Soltani carried it out and wrote the paper, and A. Zebarjadi helped us in the molecular analysis of ISSR markers.
The authors declare that they have no conflict of interest.
We would like to thank to H. M. Moradi Shahr Babak and M. Khaldari for help in collecting blood samples and also Razi University for financially supporting this project. Edited by: S. Maak Reviewed by: T. Gruszecki and two anonymous referees