Genetic Analysis of Nigerian Indigenous Goat Populations Using Microsattelite Markers

Document Type: Research Articles

Authors

1 National Agricultural Extension and Research Liaison Services, Ahmadu Bello University, Zaria, Nigeria

2 Department of Animal Science, Ahmadu Bello University, Zaria, Nigeria

3 Department of Animal Science, Federal University, Oye-Ekiti, Nigeria

4 National Animal Production Reseacrch Institute, Ahmadu Bello University, Zaria, Nigeria

5 Department of Animal Science, Kaduna State University, Kafanchan Campus, Nigeria

Abstract

This study was conducted to determine the genetic relationships among and population history of four Nigerian indigenous goat populations. A total of 200 goats from three breeds namely Sahel (60), Red Sokoto (60), West-African dwarf (60) and Kano Brown, a strain (20), were used for the study. Tissue samples were collected from the ear of animals using an Allflex® ear punch tissue sample collector, and aliquoted into plastic tubes containing the Allflex® tissue preservative. DNA extraction, amplification and sequencing were carried out at the International Livestock Research Institute (ILRI), Nairobi, Kenya. Genetic analysis of the DNA was carried out using 25 microsatellite markers proposed by Food and Agricultural Organization and International Society for Animal Genetics (FAO-ISAG). From the results of the study, the highest heterozygosity was observed in Kano Brown goats (0.68±0.04), which was followed by that of Red Sokoto goats (0.64±0.04) and that of West African dwarf goat (0.58±0.05) as the least. In the entire goat populations, low inbreeding was observed; mean Fis was 0.105 and Fit was 0.129. As expected, populations of Red Sokoto and Kano Brown showed higher genetic similarity as was seen in the genetic distance (0.022), confirming the notion that the Kano Brown is a strain from the Red Sokoto breed. Gene flow (Nm) played an important role in the genetic uniformity in populations of narrow geographical vicinity (14.868). The dendogram displayed a remarkable degree of consistency with the geographical origin of goats in this country. The information obtained in this study will aid in rational development, utilization and conservation of Nigerian indigenous goats.

Keywords


INTRODUCTION

Small ruminants play a very significant but often underestimated role in the livelihoods of rural farmers in Nigeria (Ajala et al. 2008). They are also an important source of animal protein and immediate cash, and have a variety of socio-cultural values (Odoi et al. 2000). Sahel (S), Red Sokoto (RS), Kano Brown (KB) and West African dwarf (WAD) goats are the most important goat populations found in Nigeria. The Sahel goats (Borno White or Kyalla) are common in the arid or Sahel region of the country. It is medium to large in size, long legged, and covered with coat of short fine hair and ears are long and droppy (Akpa et al. 1998). The Red Sokoto (Maradi) goat is the most numerous goat breed in Nigeria. It was originally confined to the Niger Republic and Sokoto province, but has now spread throughout the Savannah belt. The skin of Maradi goat is among the most valuable of all goat skins. It is used in the leather industry locally and internationally (Akpa et al. 1998). The Kano brown goat is an ecotype of Red Sokoto goat. It is light brown in colour with similar body weight as the Red Sokoto goat. While the West African dwarf (WAD) goat, which is trypanotolerant, is widely distributed across the rain-forest belt of Southern Nigeria. They are short-legged and small-bodied animals, weighing between 22 to 26 kg (Mourad et al. 2000). An assessment of the genetic variability in domestic goats is required to meet current production needs in various environments, to allow sustained genetic improvement, and to facilitate rapid adaptation to adverse weather conditions. Microsatellites are ideal molecular markers for characterization. Microsatellites are widely used as genetic markers for the analysis of genetic variability within and between breeds due to their high number and distribution throughout the genome and the efficacy of genotyping (Okpeku et al. 2011a). Microsatellite markers have been identified and used in several livestock species; in cattle, Sharma et al. (2015) studied the genetic diversity and relationship of Indian cattle inferred from microsatellite and mitochondrial DNA markers, similarly Karacaoren and Kadarmideen (2008) evaluated the principal component and clustering analysis of functional trait in Swiss dairy cattle while Sraphet et al. (2008) carried out a research on the use of cattle microsatellite markers to assess genetic diversity of Thai swamp Buffallo. Ibrahim et al. (2010) studied the genetic diversity in Balkhi, Hashtnagri and Michni sheep populations using SSR markers; Laval et al. (2000) studied the genetic diversity of eleven European pig breeds. Also, Mtileni et al. (2012) studied the genetic diversity and conservation of South African indigenous chicken populations. Several studies have been documented on the genetic diversity of goats based on microsatellite markers, such as Swiss breeds (Saitbekova et al. 1999), Chinese indigenous populations (Li et al. 2002; Li et al. 2004), Indian domestic goats (Rout et al. 2008), goats from Europe and Middle East (Canon et al. 2006). Scanty studies have been carried out on the genetic diversity of Nigerian goat (Adebambo et al. 2011; Okpeku et al. 2011a; Ojo, 2014; Yakubu et al. 2014). The diversity in livestock is rapidly changing due to changes in market demand, as farmers have preference for more prolific species which are more productive. It is therefore important to conserve the genetic diversity among indigenous livestock species. A Comprehensive knowledge of the existing genetic variability is the first step for conservation and exploitation of domestic animal biodiversity (Li et al. 2002). The four goat populations of Nigeria (Sahel, Red Sokoto, Kano Brown and West-African dwarf goats) have been phenotypically characterized but their genetic characterization has not been fully exploited. Hence, the present study, characterization of Nigerian goat populations at molecular level using microsatellite markers, is aimed at estimating allelic frequency of different gene loci and the genetic relationships between the Nigerian goat populations, as well as determining the time of divergence for all the four goats populations.

 

MATERIALS AND METHODS

A total of 200 randomly sampled indigenous goats from three (3) breeds namely: Sahel, Red-Sokoto and West African dwarf goat and Kano-Brown were used for the study. For each of the first 3 breeds, a maximum of sixty (60) animals were sampled, and only twenty (20) were sampled from the Kano Brown.

 

Sample collection and DNA extraction

Microsatellite markers and genotyping

Twenty five microsatellite markers were used for the genotyping. They were randomly chosen from FAO recommended list (FAO, 2011) (Table 1). The annealing temperature of the primers ranged between 50 ˚C to 65˚C. A master mix was prepared according to the manufacturers specifications containing water, buffer, dNTPs, primers and TaqDNA Polymerase in a single tube and latter aliquoted into individual tubes. Magnesium chloride (MgCl2) and template DNA solutions were then added. The sample was gently vortexed and briefly centrifuged to collect all drops from the walls of the tube and the samples were placed in a thermo cycler using the following setup: an initial denaturation at 95 ˚C (4 minutes), annealing at 50 ˚C (45 seconds) and extension at 72 ˚C (15 minutes) for 35 cycles (Ojo, 2014). After amplification, the PCR products were ran on a 1.5% agarose gel and stained with PET dye and visualized by ultra-violet light exposure. Bands of the correct size were excised from the gel, documented and sequenced. The study was conducted in accordance with institutional guidelines on the care and use of animals for scientific research, and in compliance with generally accepted rules of best practice worldwide.

 

Data analysis

The alleles and genotypic frequencies were directly identified from the gel, while genetic distance and polymorphic information content were computed using GenAlEX 6.5 genetic analysis statistical package of Peakall and Smouse (2012).

 

Table 1 Characteristics of microsatellite markers used for the study

 

 

The deviations of the locus/population combinations from Hardy Weinberg equilibrium (HWE) were determined using the statistical program of GenAlex 6.5 (Peakall and Smouse, 2012). Phylogenetic trees were constructed from Nei’s standard genetic distance (DS) and Cavalli- Sforza and Edwards genetic distance measure (DA) generated using the Neighbour Joining (NJ) and the Unpaired Group with Arithmetic Mean (UPGMA) methods in the TREEVIEW program, GENEPOP version 3.4 (Raymond and Rousset, 1995).

 

 RESULTS AND DISCUSSION

The means and standard errors of the effective alleles, observed and expected heterozygosities for each population are shown in Table 2. The mean number of different alleles, effective alleles, observed and expected heterozygosities as well as fixation index for the Sahel goat population were 10.00 ± 0.85, 4.37 ± 0.40, 0.61 ± 0.05, 0.70 ± 0.04 and 0.14 ± 0.044 respectively. Whereas for the Red Sokoto population, the mean number of different alleles, effective alleles, observed and expected heterozygosities, and the fixation index was 9.96 ± 0.80, 4.31± 0.38, 0.64 ± 0.04, 0.71 ± 0.03 and 0.09 ± 0.044, respectively. In Kano Brown population, the mean number of different alleles was 7.60 ± 0.61, while that of the effective alleles was 4.67 ± 0.43. The mean observed and expected heterozygosities for this population were 0.68 ± 0.04 and 0.72 ± 0.04, respectively and the fixation index was 0.04 ± 0.04. Also, the mean number of alleles for the West African dwarf goats was 9.52 ± 0.80, while the mean effective allele was 4.18 ± 0.42. The values obtained for the observed and expected heterozygosities for this population were 0.58 ± 0.05 and 0.68 ± 0.04, and the fixation index was 0.17 ± 0.06. The overall means of the different alleles and effective alleles were 9.27 ± 0.39 and 4.38 ± 0.20, respectively. The overall means for the observed and expected heterozygosities were 0.63 ± 0.02 and 0.70 ± 0.02, respectively, while the fixation index was 0.11 ± 0.024. In all cases, the average observed heterozygosities were lower than the expected heterozygosities. The standard genetic distances calculated for the four goat populations (Table 3) showed that the shortest distance (0.022) was observed between Red Sokoto and Kano Brown, followed by the distance between Red Sokoto and Sahel goats (0.035). However, the farthest distance was obtained between West-African dwarf and Kano Brown (0.121). The average observed and expected heterozygosities of the different goat populations as well as their Hardy-Weinberg equilibrium values are presented in Table 4. The average expected heterozygosity was lowest in West African dwarf goats (0.679), while the Kano Brown goats had the highest expected heterozygosity (0.721). Similarly, Kano Brown goats had the highest observed heterozygosity than all the other three goat breeds. In all cases, the average observed heterozygosities were lower than the expected heterozygosities and all populations showed significant deviations (P<0.01) from Hardy-Weinberg expectations. The F-statistics (Fis, Fit and Fst) and estimates of gene-flow across the goat populations are shown in Table 5. The Fis value for the entire population using the 25 markers ranged from -0.043 to 0.592. The Fis values were negative for five (5) loci (ILSTS029, ILSTS087, INRABERN 185, SPS 113 and SRYM 18) while the remaining loci revealed positive Fis values ranging from 0.012 at ILSTS11 to 0.592 at INRABERN172. The mean Fis value was 0.105 ± 0.040. The Fit value ranged from -0.016 to 0.601 and a mean of 0.129 ± 0.041, whereas the Fst value across the population ranged from 0.007 at INRA023 to 0.169 at MAF65 with a mean value of 0.03 ± 0.07. Gene flow (Nm) ranged from 1.226 for MAF65 to 35.598 for INRA023 with a mean value of 14.868 ± 1.799. The dendogram or phylogenetic tree of the Nigerian goat mitochondria DNA (mtDNA) sequences is shown in Figure 1. The tree shows that all the breeds originated from a common source or ancestry. However, West African dwarf goat diverged first and is more genetically distinct. This was followed by Sahel and lastly Red Sokoto and Kano Brown goats. Although separated, they are more genetically closely related. In the present study twenty five micro-satellite loci were used to evaluate the genetic relationships within and between the Nigerian goat populations. The twenty five micro-satellite markers were all polymorphic in the goat populations with the exception of ILST005. The different measures of genetic variation such as Na, Ne and I, were highly polymorphic across loci. Comparatively, allelic frequency in Nigerian goats of 7.6 to 10.0 was higher than those of Iranian goats (7.3) and Canary Island goats (5.9) as reported by Mahmoudi et al. (2010) and Martinez et al. (2006) respectively, but lower than Spanish Guadarrama goats with allelic frequency of 9-36 (Serrano et al. 2009). The mean observed heterozygosity was lower than those of the expected as such, suggesting a deficit in heterozygosity as observed in the F values as well. The relative low heterozygosity suggests the presence of inbreeding, probably resulting from minimal reproductive isolation and larger effects of migration and genetic drift across the entire Nigerian goat population. The close relationship observed between Red Sokoto and Kano Brown may have resulted from gene flow as these animals are in the same geographical location. The genetic distance between Red Sokoto and West African dwarf goats (0.29) was lower than the 0.39 reported by Adebambo (2003) but higher than 0.27 reported by Okpeku  et al. (2011a) for the same breeds. This may have resulted from wider geographical locations from which these animals were sampled.

 

Table 2 Mean and standard errors for Shannon index, observed and expected heterozygosities and fixation index for each population

 

N: sample size; Na: no. of alleles; Ne: no. of effective alleles; I: Shannon index; Ho: observed heterozygosity; He: expected heterozygosity and Fit: fixation index.

SE: standard error.

 

Table 3 Pair wise population matrix of Nei unbiased genetic distance

 

 

Table 4 Hardy-Weinberg equilibrium for the Nigerian goat populations

 

Ho: observed heterozygosity; He: expected heterozygosity and HWE: Hardy Weinberg equilibrium.

** (P<0.01).

 

The significant deviations from the Hardy-Weinberg Equilibrium for the Nigerian goat populations, implies that there was a deficit of heterozygotes resulting from the non-random union of gametes in the population i.e. the mating among individuals in the population which are more related than the average relationship (Ojo, 2014). Similarly, migration, mutation and selection might have caused these deviations. Deviations from Hardy-Weinberg equilibrium at microsatellite level have been reported in various studies (Barker et al. 2001; Hassan et al. 2003; Laval et al. 2000; Luikart et al. 1999). The presence of negative Fis values at five (5) loci (ILSTS029, ILSTS087, INRABERN 185, SPS 113 and SRYM 18) suggests heterozygote deficiencies which have also been reported in other studies on goats (Barker et al. 1997; Luikart et al. 1999; Agha et al. 2008; Rout et al. 2008; Dixit et al. 2009). This heterozygote deficiency may have arisen due to population sub-structure from pooling together different populations (admixture) in the analysis (Cerda-Flores et al. 2002; Muema et al. 2009). Additional factors include population subdivision owing to genetic drift, null alleles and inbreeding (Hoarau et al. 2005). The average Fis inbreeding coefficient for the entire goat population (0.11±0.04) was higher (0.04±0.13) than that obtained by Okpeku et al. (2011b) and could be a consequence of management systems (like tethering of the animals in the case of small flock size during crop season, herding of the flock by children and / or women, where the flock size is larger and some localized transhumance during the dry season). Tethering minimizes the number of different bucks that any particular female (doe) can mate, leading to service by only males preferred by the owner. While women and children herders usually avoid the mixing of different flocks or exchange of bucks from different areas even at the watering points. Localized transhumance also leads to minimal encounters between different flocks. All these practices and events tend to increase the rate of inbreeding (Ojo, 2014). The relatively low total inbreeding (Fit=0.129) observed for the entire goat population may merely be an indication of gene flow and migration between the Nigerian goat populations. This implies a gradual erosion of desirable traits within the Nigerian goat population and if not controlled, deleterious genes would begin to express themselves with time.

 

Table 5 F-statistics and estimates of gene flow (Nm) over all populations

 

Fis: within population inbreeding; Fit: total inbreeding estimate; Fst: measure of population differentiation and Nm: gene flow.

SE: standard error.

 

Figure 1 Dendogram representing the genetic relationship amongst four Nigerian populations

 

 

The Fst estimate which is commonly called the coefficient of gene differentiation obtained across all the loci implies that genotypic frequencies among the populations were not randomly distributed, meaning that the Nigerian goat populations have some genetic similarity amongst the sampled populations. The mean Fst value across all loci was low (0.03) compared to previous results of Nigerian goats (Fst=0.11, Okpeku et al. 2011b), Asian goats (Fst=0.14, Barker et al. 2001) and Chinese goats (Fst= 0.105, Li et al. 2002). However, Mujibi (2005) reported a moderate Fst of 0.058 for West African dwarf goats in Kenya. This variation may have resulted from population and breed differences. According to Laval et al. (2000), migration, interbreeding and genetic drift may exert a greater effect on the reduction in genetic differentiation between populations. The high mean gene flow (14.868) across the populations suggests mobility and considerable exchange of genetic material among these goats which may have resulted from minimal reproductive isolation and larger effects of migration resulting in low genetic differentiation and inbreeding. These could be attributed to the fact that most of these animals originated from Northern Nigeria where normadic pastoralism is the dominant livestock management system and to the extensive system of management, allowing the animals to roam freely and fend for themselves in most rural households and communities in the south. This enables and reinforces the ability of related animals to meet on pasture to breed or for neighbours to exchange related animals for upkeep or breeding (Okpeku et al. 2011b). The dendogram separated the four Nigerian goat breeds according to their geographical location in the country. Evidenced from the dendogram, supports first divergence of West African dwarf goat, followed by Sahel and latter Red Sokoto and Kano Brown goat. From the dendogram, geographically adjacent populations were observed to be more genetically related as observed in Red Sokoto, Kano Brown and Sahel goats which are Northern breeds, though Sahel goat is located at the extreme North. However, West African dwarf is a Southern breed and was observed to be more genetically distinct.

 

CONCLUSION

The close kinship between Red Sokoto and Kano Brown goats suggest some past crossing between these two populations. The result of the dendogram was consistent with the background of the origin, history and geographical location of these goat populations. The dendogram showed that West African dwarf goat population is more distinct from the other goat populations. The information obtained in this study will aid rational development, utilization and conservation of the Nigerian indigenous goats (Ojo, 2014).

 

ACKNOWLEDGEMENT

Authors thank the staff of DNA laboratory, International Livestock Research Institute (ILRI), Kenya, for the molecular analysis. Also we wish to thank Dr C. Mukasa for his valuable contribution in the collection of the samples.

Adebambo O.A. (2003). Animal breeds:A Nation’s Heritage. An Inaugural Lecture Delivered at University of Agriculture, Abeokuta, Nigeria.

Adebambo A.O., Adebambo O.A., Williams J.L., Blott S. and Urquart B. (2011). Genetic distance between two popular Nigerian goat breeds used for milk production. Livest. Res. Rural. Dev. Available at: http://www.lrrd.org/lrrd23/2/adeb23026.htm.

Agha S.H., Pilla F., Galal S., Shaat I., D’Andrea M., Reale S., Abdelsalam A.Z.A. and Li M.H. (2008). Genetic diversity in Egyptian and Italian goat breeds measured with microsatellite polymorphism. J. Anim. Breed. Genet. 125, 194-200.

Ajala M.K., Lamidi O.S. and Otaru S.M. (2008). Peri-Urban ruminant production in Northern Guinea Savannah, Nigeria. J. Anim. Vet. Adv. 3(3), 138-146.

Akpa G.N., Duru S. and Amos T.T. (1998). Influence of strain and sex on estimation of within-age-group body weight of Nigerian Maradi goats from their linear body measurements. Trop. Agric. Trinidad. 75, 462-467.

Barker J.S.F., Moore S.S., Hetzel D.J.S., Tan S.G. and Byrne K. (1997). Genetic diversity of Asian water buffalo (Bubalusbubalis): Microsatellite variation and a comparison with protein coding loci. Anim. Genet. 28, 103-115.

Barker J.S.F., Tan S.G., Moore S.S., Mukherjee T.K., Matheson J.L. and Selvaraj O.S. (2001). Genetic variation within and relationships among populations of Asian goats (Capra hircus). J. Anim. Breed. Genet. 118(4), 213-233.

Cañon J., Garcia D., Garcia-Atance M.A., Obexer-Ruff G., Lenstra J.A., Ajmone-Marsan P., Dunner S. and Econogene Consortium (2006). Geographical partitioning of goat diversity in Europe and the Middle East. Anim. Genet. 37, 327-334.

Cerda-Flores R.M., Villalobos-Torres M.C., Barrera-Saldaña H.A., Cortés-Prieto L.M., Barajas L.O., Rivas F., Carracedo A., Zhong Y., Barton S.A. and Chakraborty R. (2002). Genetic admixture in three Mexican Mestizo populations based on D1S80 and HLA-DQA1 loci. American J. Hum. Biol. 14(2), 257-263.

Dixit S.P., Verma N.K., Aggarwal R.A.K., Kumar S., Chander R.,Vyas M.K. and Singh K.P. (2009). Genetic structure and differentiation of three Indian goat breeds. Asian-Australasian J. Anim. Sci. 22(9), 1234-1240.

FAO. (2011). Food and Agriculture Organization of the United Nations (FAO). Rome, Italy.

Hassan A.A., Abou M.H.A., Oraby H.A., De Hondt E.L. and Nahas S.M. (2003). Genetic diversity of three sheep breeds in Egypt based on micro-satellite Analysis. J. Engin. Biotechnol. 1(1), 141-150.

Hoarau G., Boon E., Jongma D.N., Ferber S., Palsson J., Van Der Veer H.W., Rijnsdorp A.D., Stam W.T. and Olsen J.L (2005). Low effective population size and evidence for inbreeding in an overexploited flatfish, plaice (Pleuronectes platessa). Proc. Soc. Biol. Sci. 272, 497-503.

Ibrahim M., Ahmad S., Swati Z.A. and Khan M.S. (2010). Genetic diversity in Balkhi, Hashtnagri and Michni sheep populations using SSR markers. African J. Biotechnol. 9(45), 7617-7628.

Karacaoren B. and Kadarmideen H.N. (2008). Principal component and clustering analysis of functional traits in Swiss dairy cattle. Turkey J. Vet. Anim. Sci. 32, 163-171.

Laval G., Iannuccelli N., Legault C., Milan D., Groenen M.A., Giuffra E., Andersson L., Nissen P.H., Jorgensen C.B., Beeckmann P., Geldermann H., Foulley J.L., Chevalet C. and Ollivier L. (2000). Genetic diversity of eleven European pig breeds. Genet. Select. Evol. 32(2), 187-203.

Li M.H., Zhao S.H., Bian C., Wang H.S., Wei H., Liu B., Yu M., Fan B., Chen S.L., Zhu M.J., Li S.J., Xiong T.A. and Li K. (2002). Genetic relationships among twelve Chinese indigenous goat populations based on microsatellite analysis. Genet. Select. Evol. 34(6), 729-744.

Li X.L. and Valentini A. (2004). Genetic diversity of Chinese goat breeds based on microsatellite markers. J. Anim. Breed. Genet. 121, 350-355.

Luikart G., Biju-Duval M.P., Ertugrul O., Zagdsuren Y., Maudet C. and Taberlet P. (1999). Power of 22 microsatellite markers in fluorescent multiplexes for parentage testing in goats (Capra hircus). Anim. Genet. 30, 431-438.

Mahmoudi B., Bayat M., Sadeghi R., Babayev M. and Abdollahi H. (2010). Genetic diversity among three goat populations assessed by microsatellite DNA markers in Iran. Global Vet. 4(2), 118-124.

Martinez A.M., Acosta J., Vega-pla J.L. and Delgado J.V. (2006). Analysis of the genetic structure of the canary goat populations using microsatellites. Livest. Sci. 102, 140-145.

Mourad M., Gbanamou G. and Balde I.B. (2000). Performance of West African dwarf goats under the extensive system of production in Faranah Guinea. Pp. 227-230 in Proc. 7th Int. Conf. Goats. Tours, France.

Mtileni B.J., Muchadeyi F.C., Maiwashe A., Groeneveld E., Groeneveld L.F., Dzama K. and Weigend S. (2012). Genetic diversity and conservation of South African indigenous chicken populations. J. Anim. Breed. Genet. 128, 209-221.

Muema E.K., Wakhungu J.W., Hanotte O. and Jianlin H. (2009). Genetic diversity and relationship of indigenous goats of Sub-Saharan Africa using microsatellite DNA markers. Livest. Res. Rural Dev.  Available at: http://www.lrrd.org/lrrd21/2/muem21028.htm.

Mujibi N. (2005). Genetic characterization of West African dwarf (WAD) goats using microsatellite markers. Kenyatta Univ., Nairobi, Kenya.

Odoi A., Gathuma J.M., Gachuiri C.K., Omore A., Gitau I., Wanyangu S.W. and Murage E. (2000). Constraints to sheep and goat production in small holder mixed farms in Kenya’s central highlands. Pp. 25-30 in Proc. 9th Int. Symp. Vet. Epidemiol. Econ. Breckenridge, Colorado.

Ojo O.A. (2014). Genetic diversity of Nigerian indigenous goat breeds using microsatellite markers. Ph D. Thesis. Ahmadu Bello Univ., Zaria, Nigeria.

Okpeku M., Peters S.O., Ozoje M.O., Adebambo O.A., Agaviezor B.O., O’Neill M.J. and Imumorin I.G. (2011a). Preliminary analysis of microsatellite-based genetic diversity of goats in Southern Nigeria. Anim. Genet. Res. 49, 33-41.

Okpeku M., Yakubu A., Peters S.O., Ozoje M.O., Ikeobi C.O.N., Adebambo O.A. and Imumorin I.G. (2011b). Application of multivariate principal component analysis to morphological characterization of indigenous goats in Southern Nigeria. Acta Agric. Slov. 98(2), 101-109.

Peakall R. and Smouse P.E. (2012). GenAlex 6.5 genetic analysis in excel, genetic software for teaching and research an update. Bioinformatics. 28, 2537-2539.

Raymond M. and Rousset F. (1995). GENEPOP (version 1.2): Population genetics software for exact test and ecumenicsm. J. Hered. 86, 248-249.

Rout P.K., Joshi M.B., Mandal A., Laloe D., Singh L. and Thangaraj K. (2008). Microsatellite based phylogeny of Indian domestic goats. BMC Genet. 9, 11-16.

Saitbekova N., Gaillard C., Obexer-Ruff G. and Dolf G. (1999). Genetic diversity in Swiss goat breeds based on microsatellite analysis. Anim. Genet. 30(1), 36-41.

Serrano M., Calvo J.H., Martinez M., Marcos-Carcavilla A., Cuevas J., Gonzalez C., Jurrado J.J. and Tejada P.D. (2009). Microsatellite based genetic diversity and population structure of endangered Spanish Guadarrama goat breed. BMC Genet. 29, 61-67.

Sharma R., Kishore A., Mukesh M., Ahlawat S., Maitra A., Pandey A.K. and Tantia M.S. (2015). Genetic diversity and relationship of Indian cattle inferred from microsatellite and mitochondrial DNA markers. BMC Genet. 16, 73-80.

Sraphet S., Moolmuang B., Na-Chiangmai A., Panyim S., Smith D.R. and Triwitayakorn K. (2008). Use of cattle microsatellite markers to assess genetic diversity of Thai Swamp buffalo (Bubalus bubalis). Asian Australasian J. Anim. Sci. 21(2), 177-180.

Yakubu A., Salako A.E., Akinyemi M.O. and Abdullah A.R. (2014). Genetic  characterization of non- descript goats in Nasarawa State, Nigeria using microsatellite markers: A preliminary investigation. Pp. 25-31 in Proc. 39th Conf. Nigerian Soc. Anim. Prod.  Ibadan, Nigeria.