Open Journal of Animal Sciences
Vol.2 No.3(2012), Article ID:21371,8 pages DOI:10.4236/ojas.2012.23025

Evaluation of genetic diversity and conservation priorities for Egyptian chickens

Sherif Ramadan1,2, Boniface B. Kayang3, Eiji Inoue4, Keijiro Nirasawa5, Hiroshi Hayakawa6, Shin’ichi Ito7, Miho Inoue-Murayama1*

1Wildlife Research Center of Kyoto University, Kyoto, Japan; *Corresponding Author: mmurayama@wrc.kyoto-u.ac.jp

2Faculty of Veterinary Medicine, Benha University, Moshtohor, Egypt

3Department of Animal Science, College of Agriculture and Consumer Sciences, University of Ghana, Legon, Ghana

4Graduate School of Science, Kyoto University, Kyoto, Japan

5Animal Breeding and Reproduction Research Team, National Institute of Livestock and Grassland Science, Tsukuba, Japan

6Poultry Research Department, Gifu Prefectural Livestock Research Institute, Seki, Japan

7Faculty of Applied Biological Sciences, Gifu University, Gifu, Japan

Received 23 April 2012; revised 27 May 2012; accepted 7 June 2012

Keywords: Conservation; Egyptian Chickens; Genetic Diversity; Microsatellite; Population Structure

ABSTRACT

In this study, 21 microsatellite markers were used to genotype 196 Egyptian local chickens obtained from Fayoumi (n = 35), Dandarawy (n = 30), Baladi (n = 29), Sinai (n = 30), El-Salam (n = 36), and Golden Montazah (n = 36) strains. The results were compared to two pure commercial chicken populations reared in Japan-White Leghorn (n = 42) and Rhode Island Red (n = 43). A total of 162 alleles were observed, with an average of 7.7 alleles per locus. The average expected heterozygosity for the Egyptian chickens was 0.595. The closest pairwise Nei’s genetic distance was recorded between Sinai and Golden Montazah (0.038) and the smallest pairwise FST value (0.006) was observed between Baladi and Sinai. The most probable structure clustering of the eight studied populations was at K = 6. Baladi, Sinai and Golden Montazah strains were clustered together forming admixed mosaic cluster. Dandarawy ranked firstly and contributed the most to aggregate genetic diversity based on two prioritization methods. The information resulting from this study may be used as an initial guide to design further investigations for development of sustainable genetic improvement and conservation programs for the Egyptian chicken genetic resources.

1. INTRODUCTION

Conservation of genetic diversity is one of the main current issues in the conservation biology literature [1]. Conservation is not only about endangered breeds but also about those that are not being utilized efficiently [2]. More than 7500 different breeds of livestock are recognized globally [3]. Conservation of all livestock breeds is considered to be financially infeasible [4], so that priorities need to be set on which population/breed is to be conserved. Both genetic diversity and non-genetic criteria are important for prioritizing breeds for conservation. The non-genetic criteria include threat status and breed merit. The threat status includes risk of extinction and efficiency of the breed utilization, and breed merit includes economic or productive, ecological and sociocultural values of the breeds [5]. As a result of many years of domestication and breeding, a wide variety of chicken breeds exist today. However, an increasing number of local breeds are under threat of extinction and valuable genotypes and traits are at risk of being lost [6]. The genetic erosion of these local breeds may lead to the loss of valuable genetic variability in specific characteristics that are momentarily unimportant in commercial breeding strategies [7].

Egyptian local chickens are subdivided into three groups according to their external morphology [8]. The first group includes pure native breeds, as Fayoumi and Dandarawy. The second group includes mongrel fowl, such as the Baladi and Sinai strains, which originated from hybridization among exotic and Egyptian autochthonous chickens continued along with different times of old trade dispersal and colonization to Egypt. The third group includes improved local strains which originated from crossing between local and standardized exotic chicken strains accompanied by selection for fast growth, such as El-Salam strain [9] and for high egg production, such as Golden Montazah strain [10]. With regard to the commercial sector in Egypt, commercial broilers have contributed 63% of the total poultry production in 2005. This could reflect the substantial growing of commercial chicken industries in Egypt at the expense of native chicken resources, improvement and maintenance [8]. Egypt possesses versatile varieties of chickens including local types highly adapted to harsh conditions and thought to constitute genetic reservoirs. For instance, the Fayoumi breed has been demonstrated by several studies to possess increased resistance to coccidiosis [11] and Marek’s disease [12], and can thus be seen as a unique breed from the viewpoint of disease resistance [12]. Similarly, there is evidence for superiority in heat tolerance, of Sinai strain over White Leghorn and broiler chicks [13].

In Egypt, microsatellites marker analyses were involved in some recent studies to assess genetic diversity within and between local chicken strains [14,15]. In this study, we evaluated the genetic diversity and the breed contribution to aggregate genetic diversity as an important criterion for its conservation by utilizing three different prioritization methods in order to set the priorities for conservation of Egyptian chickens based on microsatellite genetic markers.

2. MATERIALS AND METHODS

2.1. Sample Collection and DNA Extraction

Feather samples were obtained from a total of 196 birds from six Egyptian local chicken strains: Fayoumi (n = 35), Dandarawy (n = 30), Baladi (n = 29), Sinai (n = 30), El-Salam (n = 36), and Golden Montazah (n = 36). For comparative purpose, samples were also obtained from two exotic pure chicken breeds: White Leghorn (WL, n = 42) and Rhode Island Red (RIR, n = 43). The Egyptian samples were collected from Al-Azzab poultry farms belonging to the Poultry Integrated Project at the Fayoum governorate, Egypt. Flock sizes for each breed were about 5000 birds, with sex ratios of one rooster per ten hens. Samples of White Leghorn were collected from the National Institute of Livestock and Grassland Science, Tsukuba, Japan and those of Rhode Island Red from Gifu Prefectural Livestock Research Institute, Gifu, Japan. DNAs were extracted from feather samples using the QIAGEN DNeasy Tissue Kit (QIAGEN, Valencia, CA, USA).

2.2. Microsatellite Genotyping

Molecular genotyping of the samples was carried out with a set of 21 autosomal (CA)n di-nucleotide microsatellite markers that are as uniformly distributed as possible throughout the chicken genome. These markers are from the revised set of microsatellites originally recommended by the FAO MoDAD project (http://www.fao.org/AG/AGAInfo/programmes/en/genetics/documents/ITWG3_Inf3.pdf) for diversity studies in chicken. These markers were used in multiplex PCR reactions employing the QIAGEN Multiplex PCR Kit (QIAGEN, Valencia, CA, USA). PCR was carried out in 10 µl reactions containing 20 ng of DNA template, 0.2 µM of each primer, of which the forward ones were fluorescently labelled (6-FAM, NED, and HEX) and 2x QIAGEN Multiplex PCR Master Mix. After an initial incubation of 95˚C for 15 min, PCR amplification was performed for 35 cycles consisting of 94˚C for 30 sec, 60˚C - 63˚C annealing for 90 sec, 72˚C for 60 sec, followed by a final extension of 60˚C for 30 min. Subsequently, the PCR products were electrophoresed on an ABI 3130 × l DNA Sequencer (Applied Biosystems) and the size of fragments was estimated based on 400 HD Rox size marker using the GENEMAPPER software (Applied Biosystems).

2.3. Data Analysis

Genetic diversity was assessed by calculating the observed and effective number of alleles (NA and Ne), mean number of alleles (MNA), observed heterozygosity (HO) and expected heterozygosity (HE) by using GENALEX version 6.0 [16]. Polymorphism information content (PIC) was calculated by using Molkin version 2.0 [17].  F-statistics [fixation coefficient of an individual within a subpopulation (FIS), fixation coefficient of an individual within the total population (FIT), and fixation coefficient of a subpopulation within the total population (FST)] per locus, in addition to pairwise FST [18] across the eight studied populations were calculated using GENEPOP version 3.4 [19]. Genetic distances among the eight populations were evaluated by Nei’s genetic distance [20]. A phylogenetic tree was constructed based on the Nei’s genetic distance (DA) by using the neighbor-joining (NJ) method [21]. The robustness of tree topologies was evaluated with a bootstrap test of 1000 resampling across loci. These processes were conducted using POPULATIONS version 1.2.30 software (http://bioinformatics-org/~tryphon/populations/).

We investigated the genetic structure of the sampled populations using a Bayesian clustering procedure implemented in STRUCTURE with the admixture method. We analyzed the clustering of the eight studied populations by using independent allele frequencies model [22]. We did 50 runs for each different value of K (2 ≤ K ≤ 8) with 60,000 iterations following a burn-in period of 100,000. Pairwise comparisons of the 50 solutions of each K value were run along with 50 permutations using CLUMPP software [23]. The software calculated the highest pairwise similarity index (H). CLUMPP software also outputs a mean of the permuted matrices across replicates after aligning the cluster membership coefficients of these replicate. Finally, the clustering pattern with the highest H value was graphically displayed for the selected K value using DISTRUCT software [24]. The most probable clustering numbers (best DK value) was assessed according to the equation [25]:

STRUCTURE software identified the migrants and admixed individuals. If the membership coefficient of an individual was more than 0.80, it was assigned to the cluster completely. If the value was lower than 0.80, it indicated that the individual was admixed and assigned it to two or more population clusters. The membership coefficient could also identify the migrants who had infiltrated into other chicken population clusters [26].

Different prioritization methods were utilized through measuring the breed contribution to aggregate genetic diversity as in the following:

1) According to Ollivier and Foulley [27], the contribution to between-breed diversity (CB) was computed by estimation of Weitzman values [28] based on the Nei’s genetic distance [20] with WEITZPRO [29]. Within breed contributions to diversity (CW) were calculated using the average values of within-breed expected heterozygosity as in the formula:, where Hk is the contribution to within-breed diversity (CW) of breed k, H(S) is the average internal heterozygosity of the whole set S and the average internal heterozygosity of the set excluding breed k. The aggregate diversity (D1) was obtained after weighting CB by FST and CW by 1 – FST according to the following equation: D1 = FSTCB + (1 – FST)CW. Positive contributions to diversity from a given population using the Ollivier and Foulley [27] method means that the remaining dataset decreases the overall diversity; consequently, the assessed population would be preferred for conservation.

2) According to Petit et al. [30], the rarefacted number of alleles per locus (k) was used to assess the contribution of the ith population to the total allelic richness as CT = CS + CD where CS is the contribution to the total allelic richness due to the allelic richness of the ith population and CD is the contribution due to its divergence. Positive contributions to diversity from ith population using the method of Petit et al. [30] mean that the remaining set has a lower number of alleles than the original set; consequently, the ith population would be preferred for conservation. This procedure was computed using Molkin version 2.0 [17].

3) According to Caballero and Toro [31], the partitions of the total gene diversity was calculated as in the following equation: (1 – ƒ) = (1 – t) + D, where f is the average global coancestry; ƭ is the average coancestry between populations; (1 – ƒ) = GDT representing the total gene diversity; (1 – t) = GDw, representing the within population component; and D is the Nei genetic distance between populations, representing the between population component. Positive contributions to diversity from a given population using the method of Caballero and Toro [31] mean that the remaining dataset increases the overall diversity; consequently, the assessed population would not be preferred for conservation. This procedure was computed using Molkin version 2.0 [17].

3. RESULTS AND DISCUSSION

3.1. Marker Polymorphisms and Population Diversity

A total of 162 alleles were observed across all the eight populations, out of which 144 alleles (144/162, 88.9%), including 18 unique ones (18/144, 12.5%), were observed in the six Egyptian populations (Tables 1 and 2). In this study, across the six Egyptian populations the estimated means of NA (6.9), Ne (3.0) and HE (0.595) are relatively lower than those of Eltanany et al. [15] who reported values of 7.34, 3.00 and 0.653, respectively, across ten Egyptian chicken strains using 29 microsatellites loci. The FST value across the 21 studied loci showed a relatively high mean (0.082) indicating that there is genetic differentiation among the six Egyptian local strains. The estimated FST value was lower than that measured between pure-bred commercial chicken lines in a study in Zimbabwe which showed 0.357 of total genetic variation owing to line differences [32]. However, it was slightly higher than the 0.068 previously reported across ten Egyptian chicken strains [15]. In this study, FST recorded a high value (0.222) after adding the two exotic pure populations (WL and RIR) indicating that, there is high genetic differentiation between these two exotic pure breeds and Egyptian chicken populations (Table 1). The relatively low but positive FIS average (0.051), in addition to the eleven loci showing a deficit of heterozygote might indicate non-random mating and also these loci might be under morphological or productive traits of selective interest. Moreover, FIS is used to obtain a deeper insight to appraise the degree of inbreeding and endangerment potentiality and is considered as an important tool to judge the conservation priority [33]. Accordingly, when FIS is less than 0.05, the breeds are not in danger; between 0.05 - 0.15, they are potentially endangered; between 0.15 - 0.25, they are minimally endangered; between 0.25 - 0.40, they are endangered; and more than 0.40, they are critically endangered. In this study, Fayoumi, Dandarawy, El-Salam and RIR populations showed high levels of inbreeding (0.110, 0.053, 0.095 and 0.083, respectively), posing their potential endangerment [33].

Table 1. Observed (NA) and effective (Ne) number of alleles, polymorphism information content (PIC), observed (HO) and expected (HE) heterozygosities, and F-statistics (FIS, FST, and FIT) across the six Egyptian strains.

Table 2. Mean observed (MNA) and effective (MNe) number of alleles, unique alleles, observed (HO) and expected (HE) heterozygosities, and fixation coefficient of an individual within a subpopulation (FIS) per breed.

In respect to the within population genetic diversity, the studied eight chicken populations could be categorized into a low diversity class (Fayoumi, WL, and RIR) and a high diversity class which includes the remaining five populations. This is in agreement with breed history and management. These populations (WL and RIR) had undergone selection for high growth rate (RIR) and high egg production (WL). Moreover, the Fayoumi strain recorded the highest value of FIS (0.110) and complete allele fixation (monomorphic) of the MCW014 locus (Table 2). This might be attributed to its narrow genetic base as it is an ancient native chicken bred as a closed population. The two Egyptian mongrel strains (Baladi and Sinai) recorded the highest genetic diversity (MNA = 5.9; Ne = 3.4; HO = 0.622, and HE = 0.645 for Baladi and MNA = 5.4; Ne = 3.4; HO = 0.648, and HE = 0.660 for Sinai) among the eight studied populations and this might be attributed to their wide genetic bases due to hybridization among exotic and Egyptian autochthonous chickens continued along with different times of old trade dispersal and colonization to Egypt [8].

3.2. Genetic Relationship

The Nei’s genetic distance (DA) and pairwise FST statistic were estimated for the eight studied chicken populations across the 21 microsatellite loci (Table 3). The closest pairwise Nei’s genetic distance was recorded between the Sinai and Golden Montazah strains (0.038) and this was supported by clustering in the neighbor-joining phylogenetic tree (Figure 1). Similarly, the lowest pairwise FST value was recorded between the Baladi and Sinai strains (0.006). The close relation between Sinai and Golden Montazah and also between Baladi and Sinai can be attributed to the mongrel nature of Baladi and Sinai strains which originated from hybridization among exotic and Egyptian autochthonous chickens.

3.3. Population Structure and Individual’s Assignment

The most probable structure clustering of the eight studied chicken populations was at K = 6 (Figure 2). The pure breeds (WL, RIR, Fayoumi and Dandarawy) in addition to El-Salam were assigned independently into their respective clusters while the remaining three populations (Baladi, Sinai and Golden Montazah) were clustered together forming admixed mosaic cluster. A probable explanation for the separation of El-Salam to form its own cluster is that it might have experienced high inbreeding and low gene flow from other strains under this study. This is in line with the relatively high FIS (0.095) value within the El-Salam strain. The high genetic admixture and migrations between Baladi, Sinai, and Golden Montazah strains could contribute to gather them forming the admixed mosaic cluster.

Table 3. Nei’s genetic distance (DA: above diagonal) and pairwise FST (below diagonal) estimates for the 21 microsatellite loci between the eight studied chicken strains.

Figure 1. Neighbor-joining tree of the six Egyptian and the two exotic pure chicken populations based on the 21 microsatellite loci. The consensus tree was generated with 1000 bootstraps over loci and bootstrap values lower than 50 are not shown in the diagram.

The distribution of 35 admixed and five migrant individuals in the inferred six clusters according to their membership coefficients were evaluated (Table 4). All the studied populations recorded more than 0.80 membership coefficients in their inferred clusters except Baladi strain (0.69). Despite that Sinai and Golden Montazah strains had more than 0.80 membership coefficient, they were clustered together in the same cluster with Baladi strain. This might be explained by the genetically similar nature of these strains as confirmed by their pairwise genetic distances (Table 3) and the mosaic admixture in the STRUCTURE dendrogram (Figure 2). Thus, based on the preceding results, when the Baladi, Sinai, and Golden Montazah chicken strains are to be used for future research work it would be necessary to consider, in addition to random selection, the admixed individuals and migrants.

3.4. Conservation Prioritization of the Studied Strains

In the current study different prioritization methods were utilized to measure the breed contribution to aggregate genetic diversity as an important criterion for its conservation (Table 5). All such methods revealed that Fayoumi strain contributed negatively to aggregate genetic diversity (CW = –4.20, D1 = –1.15, D2 = –1.89 and GD = 1.72). Therefore, Fayoumi according to such a determined criterion may be ranked last for conservation.

On the contrary, Dandarawy contributed the most (CB = 34.92, D2 = 2.49, GD = –1.40) to aggregate genetic diversity according to [27,30,31]. The two Egyptian mongrel strains ranked second (D1 = 2.90, GD = –1.23 for Sinai and D1 = 1.73, D2 = 2.27 for Baladi), while Egyptian synthetic strains (Golden Montazah then El-Salam) came in the third position according to their contribution to aggregate genetic diversity. The preceding prioritization of the breeds for conservation is based only on molecular genetic marker information, but when we combine other non-genetic criteria the ranking may become different. Thus, according to the preceding prioritization methods, Fayoumi ranked the last, but after considering its high level of inbreeding and breed merit in term of disease resistance ability (Marek’s disease and coccidiosis), it may get advanced ranking. Similarly, Sinai strain ranked second, but after considering its breed merit (superiority in heat tolerance) it may get a different ranking. 

Figure 2. Structure clustering of the six Egyptian and two exotic pure chicken populations obtained for K = 6. The percentage inside the parenthesis is the average pairwise similarity index (H) of the individuals Q matrix, while K is the cluster number.

Table 4. Number of admixed and migrant individuals in the inferred clusters.

Table 5. Contribution of each strain to aggregate genetic diversity.

In conclusion, the results from this study confirm the applicability and efficiency of this microsatellite panel for assessing genetic variation and setting the conservation priorities for Egyptian local chickens. Consideration of breed merits and threat status, in addition to genetic diversity, enabled us to balance the trade-offs between conserving diversity as insurance against future uncertainties and current sustainable utilization. More detailed information about non-genetic aspect (threat status and breed merits) and a conceptual framework for a maximum utility through a weighted summation of measures of neutral diversity, breed merits and threat status of Egyptian chickens merits consideration.

4. ACKNOWLEDGEMENTS

We thank Dr. Ahmed Radwan, Director of Fowl Integrated Project, Fayoum Governorate, Egypt, for the great facilities he kindly offered during the feather sampling of Egyptian chickens. We are grateful to Dr. Azusa Hayano, Wildlife Research Center of Kyoto University for her assistance in data analysis. This study was supported financially in part by the Ministry of Education, Culture, Sports, Science and Technology (MEXT) with a Grant-in-aid for Science Research (#21310150 to MI-M) and the Global Center of Excellence Program “Formation of a Strategic Base for Biodiversity and Evolutionary Research: from Genome to Ecosystem”.

REFERENCES

  1. Frankham, R. (1995) Conservation genetics. Annual Review of Genetics, 29, 305-327. doi:10.1146/-annurev.ge.29.120195.001513
  2. Barker, J.S. (2001) Conservation and management of genetic diversity: A domestic animal perspective. Canadian Journal of Forest Research, 31, 588-595. doi:10.1139/x00-180
  3. FAO (2007) The state of the world’s animal genetic resources for food and agriculture. FAO, Rome. http://www.fao.org/docrep/010/a1250e/a1250e00.htm
  4. Bennewitz, J., Eding, H., Ruane, J. and Simianer, H. (2007) Selection of breeds for conservation. In: Oldenbrook, K., Ed., Utilization and Conservation of Farm Animal Genetic Resources, Wageningen Academic Publishers, Wageningen, 131-146.
  5. Ruane, J. (2000) A framework for prioritizing domestic animal breeds for Conservation purposes at the national level: A Norwegian case study. Conservation Biology, 14, 1385-1393. doi:10.1046/j.1523-1739.2000.99276.x
  6. Blackburn, H. (2006) The National Animal Germplasm Program: Challenge and opportunities for poultry genetic resources. Poultry Science, 85, 210-215.
  7. Weigend, S., Vef, E., Wesch, G., Meckenstock, E., Seibold, R. and Ellendorff, F. (1995) Concept for conserving genetic resources in poultry in Germany. Archiv für Geflügelkunde, 59, 327-334.
  8. Hosny, F.A. (2006) The structure and importance of the commercial and village based poultry systems in Egypt. In: Food and Agriculture Organization of the United Nations, Ed., Poultry Sector Country Review, 1-39.
  9. Abd El-Gawad, E.M., Balat, M.M., AbouElElla, N.Y., Ali, M.M. and Omran, Kh.M. (1983) “El-Salam” a new locally developed strain of chickens. Agricultural Research Revolution, 61, 147-157.
  10. Mahmoud, T.H., Sayed, I.F. and Madkour, Y. H. (1974) “The Golden Montazah” a new variety of chickens. Agricultural Research Revolution Cairo, 52, 97-105.
  11. Pinard-Van Der Laan, M.H., Monvoisin, J.L., Pery, P., Hamet, N. and Thomas, M. (1998) Comparison of outbred lines of chickens for resistance to experimental infection with coccidiosis (Eimeria tenella). Poultry Science, 77, 185-191.
  12. Tixier-Boichard, M., Bordas, A. and Rognon, X. (2009) Characterisation and monitoring of poultry genetic resources. World’s Poultry Science Journal, 65, 272-285. doi:10.1017/S0043933909000233
  13. Arad, Z., Marder, J. and Soller, M. (1981) Effect of gradual acclimation to temperatures up to 44˚C on productive performance of the desert Bedouin fowl, the commercial White Leghorn and the two reciprocal crossbreds. British Poultry Science, 22, 511-520. doi:10.1080/00071688108447918
  14. Roushdy, Kh., Zein El-Dein, A., Fathi, M.M., Ali, U.M. and Assy, H.M. (2008) Microsatellite genetic differentiation analysis of two local chicken breeds compared with foreign Hy-line strain. International Journal of Poultry Science, 7, 1045-1053. doi:10.3923/ijps.2008.1045.1053
  15. Eltanany, M., Philipp, U., Weigend, S. and Distl, O. (2011) Genetic diversity of ten Egyptian chicken strains using 29 microsatellite markers. Animal Genetics, 42, 666-669. doi:10.1111/j.1365-2052.2011.02185.x
  16. Peakall, R. and Smouse, P.E. (2006) GENALEX 6: Genetic analysis in Excel. Population genetic software for teaching and research. Molecular Ecology Notes, 6, 288- 295. doi:10.1111/j.1471-8286.2005.01155.x
  17. Gutiérrez, J.P., Royo, L.J., Álvarez, I. and Goyache, F. (2005) Molkin v2.0: A computer program for genetic analysis of populations using molecular coancestry information. Journal of Heredity, 96, 718-721. doi:10. 1093/jhered/esi118
  18. Weir, B.S. and Cockerham, C.C. (1984) Estimating F-statistics for the analysis of population structure. Evolution, 38, 1358-1370. doi:10.2307/ 2408641
  19. Raymond, M. and Rousset, F. (1995) GENEPOP (version 1.2): Population genetics software for exact tests and ecumenicism. Journal of Heredity, 86, 248-249.
  20. Nei, M., Tajima, F. and Tateno, Y. (1983) Accuracy of estimated phylogenetic trees from molecular data. II. Gene frequency data. Journal of Molecular Evolution, 19, 153-170. doi:10.1007/BF02300753
  21. Saitou, N. and Nei, M. (1987) The neighbor-joining method: A new method for reconstructing phylogenetic trees. Molecular Biology and Evolution, 4, 406-425.
  22. Pritchard, J.K., Stephens, M. and Donnelly, P. (2000) Inference of population structure using multilocus genotype data. Genetics, 155, 945-959.
  23. Jakobsson, M. and Rosenberg, N.A. (2007) CLUMPP: A cluster matching and permutation program for dealing with label switching and multimodality in analysis of population structure. Bioinformatics, 23, 1801-1806. doi:10.1093/bioinformatics/btm233
  24. Rosenberg, N.A. (2004) DISTRUCT: A program for the graphical display of population structure. Molecular Ecology Notes, 4, 137-138. doi:10.1046/j.1 471 -8286.2003.00566.x
  25. Evanno, G., Regnaut, S. and Goudet, J. (2005) Detecting the number of clusters of individuals using the software STRUCTURE, a simulation study. Molecular Ecology, 14, 2611-2620. doi:10.1111/j.1365-294X. 2005.02553.x
  26. Li, H.F., Han, W., Zhu, Y.F., Shu, J.T., Zhang, X.Y. and Chen, K.W. (2009) Analysis of genetic structure and relationship among nine indigenous Chinese chicken populations by the Structure program. Journal of Genetics, 88, 197-203. doi:10.1007/s12041-009-0028-8
  27. Ollivier, L. and Foulley, J.L. (2005) Aggregate diversity: new approach combining within and between-breed genetic diversity. Livestock Production Science, 95, 247-254. doi:10.1016/j.livprodsci.2005.01.005
  28. Weitzman, M.L. (1993) What to preserve? An application of diversity theory to Crane conservation. Quarterly Journal of Economics, 108, 157-183. doi: 10.2307/2118499
  29. Derban, S., Foulley, J.L. and Ollivier, L. (2002) WEITZPRO: A software for analyzing genetic diversity. Institut National de la Recherche Agronomique, Paris. http://www23sgqa.jouy.inra.fr/article.php3?id_article=3
  30. Petit, R.J., El Mousadik, A. and Pons, O. (1998) Identifying populations for conservation on the basis of genetic markers. Conservation Biology, 12, 844-855. doi:10.1046/j.1523-1739.1998.96489.x
  31. Caballero, A. and Toro, M.A. (2002) Analysis of genetic diversity for the management of conserved subdivided populations. Conservation Genetics, 3, 289-299. doi:10.1023/A:1019956205473
  32. Muchadeyi, F.C., Eding, H. Wollny, C.B., Groeneveld, E., Makuza, S.M., Shamseldin, R., Simianer, H. and Weigend, S. (2007) Absence of population substructuring in Zimbabwe chicken ecotypes inferred using microsatellite analysis. Animal Genetics, 38, 332-339. doi:10.1111/j.1365-2052.2007.01606.x
  33. Simon, D.L. and Buchenauer, D. (Eds.) (1993) Genetic diversity of European livestock breeds. Results of Monitoring by the EAAP Working Group on Animal Genetic Resources; EAAP Animal Genetic Data Bank, Institute of Animal Breeding and Genetics, School of Veterinary Medicine, Hannover, EAAP Publication, No. 66, Wageningen Pers, Wageningen, 591 p.