Genetic Diversity of Indigenous Chicken (Gallus Gallus domesticus) from Ecozones of Egypt and Kingdom of Saudi Arabia ()
1. Introduction
Chickens’ ecotypes are being utilized in rural and backyard farming of developing countries all over the world. The significance of these ecotypes in the rustic economy is tremendous in various counties [1] - [7]. This significance is due to affordable low production cost, scavenging competency as well as adaptability to harsh and stressful environmental conditions [8] [9]. Moreover, chicken ecotypes supply work security to the family as well as getting the accessibility of food. Jobless youth and ladies can likewise acquire pay through poultry cultivating [10]. Still, local ecotypes have low egg productivity and light mature body size [11] [12].
The poultry industry is globally dominated by commercial chicken breeds, where breeding strategy is based on few chicken genotypes. Such a strategy has the shortcoming of erosion of valuable genetic resources (ecotypes) [13]. Their low productivity does not meet the global demands of intensive production systems, as such ecotypes could barely compete with commercial breeds. Therefore, ecotypes are negatively selected regardless of their good quality of egg and meat, disease resistance as well as adaptation to the local environment [13]. The establishment of a framework for preserving these genetic resources is of need [8]. Investigating the genetic makeup of indigenous chicken breeds was carried out in several studies [11] [14] [15] [16] and has been used as a sampling guide in genetic diversity research [17].
At Taif governorate (~1.7 km above sea) local chickens are acclimated to the rough environment of high altitudes, such as low oxygen tension, as well as daily temperature fluctuations [18] [19] [20] [21].
Microsatellites (MS) markers are normally used because they are plentiful, codominant, randomly spread throughout the genome, and very polymorphic, also because of their reproducibility [22] [23]. This coherence has made MS markers a dependable tool for genetic diversity evaluation and assessment in many commercial chicken breeds and ecotypes [13] [18] [19] [20]. Quite a few MS markers are available in chickens and have been utilized to develop linkage maps in numerous chicken breeds [22] [24]. These markers introduced efficient machinery to QTL research and have additionally been effectively utilized to consider the hereditary associations among and within chicken populations [23].
Osaman et al., [25] utilized the complete sequence of mitochondrial DNA D-loop to explain the genetic descent of Egyptian indigenous chicken and Asian chicken. Results of this study revealed that both Egyptian native chicken and West and Central Asian chicken are sharing the same common ancestor as they branched together in the same clade. However, this work did not include native chicken strains from Saudi Arabia. Therefore, the current study was carried out as an attempt to introduce a better understanding of the genetic characterization of 2 ecotypes sampled from two different ecozones at Taif (Saudi Arabia) and Ismailia (Egypt). A dense microsatellites panel of 39 markers was genotyped for two ecotypes.
2. Materials and Methods
2.1. Sample Collection
Chicken samples were collected from individual framers in local village either in Ismailia (ISM) or Taif (TA) with same age as possible. From each location 25 chicken were collected and blood sample collected from chicken using EDTA contain collected tube.
2.2. DNA Extraction
Collected blood samples were used for DNA isolation using DNA extraction kit (QIAGEN) according to manufacture manual. Extracted DNA quality checked by electrophoresis in a minigel while the quantity checked spectrophotometrically using NanoDrop 2000C (Spectronic Genesys, Thermo Electron Corporation). DNA samples were stored at −20˚C until microsatellites analysis.
2.3. Microsatellite Genotyping
Thirty nine MS markers were utilized to assess DNA polymorphism of these 2 ecotypes. These MS markers were chosen based on reports by the International Society of Animal Genetics (ISAG)-FAO to think about the hereditary assorted variety of chickens [26] and additionally 150 microsatellites were investigated for the segregation of five thoroughbred KNC lines [20]. Names of each of the 39 MS markers, chromosome number and respective allelic size are appeared in Table 1.
![]()
Table 1. Marker’s name, Chromosomal number (Chr) and allele size (bp) of 39 genotyped Microsatellites markers.
Polymerase Chain Reaction (PCR) amplifications were in absolute volume of 20 µL, 50 ng of genomic DNA, 10 pmol of fluorescent dye (FAM, VIC, NED, PET) labeled modified forward primer and normal reverse primer (Applied Biosystems, USA), 2.5 mM of each dNTPs (GeNet Bio, Korea), 10 X reaction buff (GeNet Bio, Korea), 2.5 unit of prime Taq DNA polymerase (GeNet Bio, Korea). The PCR was performed in an initial denaturation at 95C for 10 min followed by 35 cycles of 30 sec of denaturation at 95C, 30 sec of annealing at 60C, 30 sec of extension at 72C and final extension at 72C for 10 min using My-Genie 96 Thermal Cycler (Bioneer, Korea). The PCR products were initially electrophoresis on 3% agarose gel with ethidium bromide (EtBr) and confirmed whether they gave single PCR DNA band under the UV light. When the bands were clearly appeared, further genotyping will be performed. For the microsatellite genotyping, more than 20 times diluted PCR products will be used. The genotyping reaction contained 1 µL of diluted PCR products, 10 µl of Hi-Di Formamide (Applied Biosystems, USA) and 0.1 µL of GeneScan-500 LIZ size standard marker (Applied Biosystems, USA). After dilution, genotyping reaction mixture was denatured for 2 min at 95C and fragment analysis was performed using capillary array in Genetic analyzer 3130xl (Applied Biosystems, USA). The MS genotypes will be identified using GeneMapper ver.3.7 (Applied Biosystems, USA).
3. Data Analyses
Exploratory data analysis and graphical representations of the results were carried out under R statistical environment R Core Team [27]. Fundamental measures of genetic diversity, such as total number of alleles, allele frequencies, mean number of alleles, observed, and expected heterozygosity (Nei, 1987) were computed using adegenet R statistical package [28]. Deviation from Hardy-Weinberg equilibrium (HWE) using chi square test, and phylogenic analysis using Nei’s distance [29] as well as population subdivision was examined using Weir and Cockerham [30] unbiased estimator of Wrights fixation indices (FIT, FST and FIS) [31] was carried out using “pegas” package R Population and Evolutionary Genetics Analysis System [32]. The null hypothesis was that the estimates were not significantly different from zero, and the level of significance (P < 0.01) was adjusted using Bonferroni correction. Allelic richness number was estimated using hierfstat package [33]. Phylogenetic analysis was carried out using APE package [34].
4. Results and Discussion
Figure 1 & Figure 2 show a descriptive overview of number of alleles per locus, Polymorphic Information Content (PIC) observed and expected Heterozygosity (Hobs and Hexp) for both ISM and TA ecotypes.
For the number of alleles per locus a considerable difference was observed in both ecotypes. The ISM ecotype has a sum of 157 alleles, number of alleles per
![]()
Figure 1. (a) Number of Alleles, (b) Polymorphic Information Content “PIC”, (c) Observed heterozygosity “Hobs” and (d) Expected heterozygosity “Hexp” of ISM ecotype.
locus ranged from 1 to 10 alleles, averaged 4.2 ± 2.2. The TA ecotype exhibited a lower number of alleles sum 138, the number of alleles per locus ranged from 1 to 7 averaged 3.6 ± 1.6. These estimates are lower than the reported on other studies, e.g. Abebe et al., [13] on Swedish native chicken (4.7), Muchadeyi et al., [16] on Zimbabwe ecotypes (9.7), Van-Marle-Koster et al., [35] on South Africa native chicken (6.1), and Ramadan et al., [36] on native Egyptian and commercial chicken breeds (7.7). These lower estimates of the present study could be attributed to utilizing a larger number of loci. The number of alleles per locus could also be represented as “allelic richness”. The number of observed alleles is basically a function of sample size. Allelic richness (Rt) was estimated across all loci (Table 2). One definition of Rt is the mean number of alleles per locus [37]. The significance of Rt is that decrease in Rt might result in reducing ability of the population to acclimate to future environmental changes, as such an assorted variety is the rough material for evolutionary forces, for example, common determination [38].
The mean of polymorphic information content (PIC) for ISM and TA ecotype were 0.47 ± 0.21, and 0.41 ± 0.2. This means that this penal of markers are
![]()
Figure 2. (a) Number of Alleles, (b) Polymorphic Information Content “PIC”, (c) observed heterozygosity “Hobs” and (d) Expected heterozygosity “Hexp” of TA ecotype.
reasonably informative as PIC exceeds 0.25 Botstein et al., (1980). For ISM ecotype 54% of markers are highly informative as PIC was higher than 50%, where only 38% of TA’s markers were highly informative. The mean value of PIC in this study was smaller than the estimates of other studies on different native breeds e.g., Seo et al., [20] estimate of PIC was 0.77 on five native Korean lines, Ramadan et al., [36] reported 0.65 on six Egyptian native breeds, and Choi et al., [39] on commercial Korean native breed reported mean PIC value of 0.68, and Soltan et al., [40] estimation of PIC was 0.84 on Sinai and Norfa Egyptian native chicken. However, lower estimate of 0.55 was reported by Abebe et al., [13] on five local Swedish breeds. Again, these differences in mean value of PIC could be ascribed to the dense panel of markers were used in the present study.
The presence of private alleles was also investigated (Table 3). Only 4 markers did show private alleles between the two ecotypes, namely ADL278, MCW087, MCW222 and MCW295. For ISM ecotype the highest number of private alleles was 9, where the highest number of private alleles for TA ecotype was only 5. Existence of private alleles might be indicating of diversifying selection of these specific alleles. In general, the percentage of shared alleles between the two
![]()
Table 2. Allelic Richness (Rt) of genotyped Microsatellites markers across ISM and TA ecotypes.
![]()
Table 3. Number of private alleles of genotyped Microsatellites markers across ISM and TA ecotypes.
ecotypes was 45%. The incidence of heterozygous individuals in population is measured by the observed proportion of heterozygosity (Hobs). A similar average values of (Hobs) for both ecotypes were attained although narrower range of values were found for TA ecotype. For the ISM ecotype the Hobs ranged from 0.0 to 0.8 with average of 0.27 ± 0.24. For the TA ecotype the Hobs ranged from 0.0 to 0.64 with average of 0.27 ± 0.23. The expected heterozygosity (Hexp) is an estimate of the likelihood that a pair of gametes arbitrary picked from genetic stock are of various alleles. In the contrary to Hobs the averages of Hexp differed between the two ecotypes, 0.52 ± 0.23 vs. 0.45 ± 0.21 for ISM and TA. This variability of the two parameters might be result of variation in evolutionary process on each of the two ecotypes.
Test for deviation from Hardy Weinberg Equilibrium was also carried out. For the ISM ecotypes only 8 loci were significantly deviated from Hardy Weinberg expectation, where 12 loci of TA ecotype.
Breaking up the genetic diversity within and among populations is best described through Wright’s F-statistics (Wright, 1951). Table 4 shows Wright’s fixation indices (FST, Fit, and FIS) across ISM and TA ecotypes for the genotyped microsatellites markers. Figure 3 shows a boxplot for the estimates of Wright’s
![]()
Table 4. Wright’s fixation indices (FST, Fit, and FIS) across ISM and TA ecotypes for the genotyped Microsatellites markers.
![]()
Figure 3. Wright’s fixation indices (FST, FIT, and FIS) across ISM and TA ecotypes across 39 microsatellites loci.
FST, F1T and F1S. The average values of the three indices were all positive. The average of FST was small (0.06) which indicates similarity of allele frequencies within each breed, which in turn indicate that genetic variation among ecotypes accounted for only 6%, which shows homology of allele frequencies of both ecotypes [41], thus indicate that genetic variation between ecotypes represents 6% of entire genetic variation. This small estimate of FST could be considered as indicative that none of the 39 MS loci were under selection [42] [43]. The estimate of FST ranged from −0.02 to 0.32. This estimate was lower than the estimates of Abebe et al. [13] on Swedish local breed, Seo et al. [20] on Korean chicken as well as Ramadan et al., [36], Soltan et al., [40] on Egyptian native breeds.
Average value of global heterozygosity index (FIT) was 0.102, ranged from −0.04 to 0.48. This index could be considered as correlation between alleles within individuals with respect to combined sample [41]. The inbreeding coefficient (FIS) across the 39 MS loci averaged 0.65 ranged from 0.14 to 1.00. This estimate is higher than the other studies on indigenous chicken breeds e.g. 0.187 on Swedish breeds [13], 0.0093 on Korean indigenous breeds [20] and 0.018 on Egyptian indigenous breeds [36], but Soltan et al., [40] reported somewhat closer estimate on Sinai and Norfa indigenous Egyptian chicken breeds. FIS is also taken as an indicative of endangerment potentiality and devised to assess the conservation priorities [44]. Thus, when FIS is lower than 0.05, breed is not in danger; the range from 0.05 to 0.15, breed is potentially endangered; between 0.15 - 0.25, they are slightly endangered; between 0.25 - 0.40, they are endangered; and more than 0.40, breed is seriously endangered [36] [40].
Diversity assessment between population is normally quantified by estimation of genetic distance. Both Nei’s genetic distance [29] as well as Nei’s pairwise FST index between all pairs of populations [45] was estimated for the two ecotypes over the 39 MS loci. The estimate of Nei’s genetic distance was 0.2 where Nei’s pairwise FST was 0.29. This small value is comfortable with the percentage of shared alleles of the two ecotypes 45%.
5. Conclusion
The results of the present study are, presumably, the first to recount the genetic diversity between two ecozones Egyptian and Saudi Arabia ecotypes. To be concluded, a small proportion of genetic diversity due to allele frequency differences was between these two ecotypes was only 6%. This result is also supported by small estimate of Nei’s genetic distance (0.2) as well as small Nei’s pairwise FST (0.29). The percentage of shared alleles of the two ecotypes was 45%.