Go to Top Go to Bottom
Maw, Kawabe, Shimogiri, Rerkamnuaychoke, Kawamoto, Masuda, and Okamoto: Genetic Diversity and Population Structure in Native Chicken Populations from Myanmar, Thailand and Laos by Using 102 Indels Markers

Abstract

The genetic diversity of native chicken populations from Myanmar, Thailand, and Laos was examined by using 102 insertion and/or deletion (indels) markers. Most of the indels loci were polymorphic (71% to 96%), and the genetic variability was similar in all populations. The average observed heterozygosities (HO) and expected heterozygosities (HE) ranged from 0.205 to 0.263 and 0.239 to 0.381, respectively. The coefficients of genetic differentiation (Gst) for all cumulated populations was 0.125, and the Thai native chickens showed higher Gst (0.088) than Myanmar (0.041) and Laotian (0.024) populations. The pairwise Fst distances ranged from 0.144 to 0.308 among populations. A neighbor-joining (NJ) tree, using Nei’s genetic distance, revealed that Thai and Laotian native chicken populations were genetically close, while Myanmar native chickens were distant from the others. The native chickens from these three countries were thought to be descended from three different origins (K = 3) from STRUCTURE analysis. Genetic admixture was observed in Thai and Laotian native chickens, while admixture was absent in Myanmar native chickens.

INTRODUCTION

Native chickens may be regarded as greatly diversified populations because of their long-term adaptation in response to varied agro-ecological zones. Native chickens possess unique adaptive traits that permit them to survive and reproduce under the harsh climatic, nutritional, and management conditions typically associated with low input-output production systems (Mwacharo et al., 2006). There are many kinds of native chickens in Asia. Most of them have not been improved and have lower productivity than the improved foreign breeds. Many kinds of native chickens can be found in all areas where farmers live. Farmers rear native chickens as important protein sources. Chickens were domesticated from the Red jungle fowl in Southeast Asia (Nishibori et al., 2005). Myanmar, Thailand, and Laos are three neighboring countries in Southeast Asia that share their borders on the Mekong River, which originates far away, in Quinghai Province, China. Native chickens have been domesticated in these countries for generations and have yielded a wide variety of chickens. It is very interesting to assess the genetic diversity of native chickens in those countries where native chickens are known to have been domesticated in the region of their ancestral species, Red jungle fowl.
Various tools offer a new approach to investigate both the genetic diversity and phylogenetic relationship among the sub-species of Gallus gallus and the domestic chicken. Some studies using blood protein variation in native chicken populations from Laos (Okamoto et al., 1999) and Myanmar (Okamoto et al., 2004) have been reported. Mitochondrial DNA sequencing is also a powerful tool for population genetic analysis. Liu et al. (2006) reported the geographical distribution of a wide range of domestic chickens in Eurasian regions. Kawabe et al. (2014) reported the genetic diversity of Laotian native fowls by mtDNA variation. The genetic characterization of Thai indigenous chickens and commercial chickens has been assessed by 20 microsatellite loci (Dorji et al., 2011). In addition, the genetic variation and phylogeographic analysis of native chicken populations from Myanmar and Thailand have been studied using 98 autosomal single nucleotide polymorphism (SNP) markers (Riztyan et al., 2012).
Insertion and/or deletion (indels) polymorphisms are diallelic markers with the potential for use in the study of genetic diversity and phylogenic relationships. Indels markers have many genetic advantages for analytical use, they are wide spread throughout the genome, all of the polymorphisms derive from a single mutation event, and they have reduced mutation rates (Natalle et al., 2010). In our previous study, we analyzed the genetic diversity of Myanmar and Indonesian native chickens together with two jungle fowl species using 102 indels polymorphisms (Maw et al., 2012). We found that genetic variability was higher among native chicken populations than in Red jungle fowl (Gallus gallus bankiva) and Green jungle fowl (Gallus varius). The native chickens from the two countries were genetically close to each other and distinct from Red and Green jungle fowl from Java Island, Indonesia.
In this study, we aimed to assess the genetic diversity and to clarify the population structure of native chicken populations from Myanmar, Thailand, and Laos by using 102 indels polymorphisms from previous studies as genetic markers.

MATERIALS AND METHODS

Birds and DNA extraction

A total of 278 blood samples were collected from nine native chicken populations, comprising three populations from each of three countries. The blood samples of Myanmar native chickens were collected from Yangon (YG) and Pegu (PG, Lower Myanmar) and from Mandalay (MD, Central Myanmar). From Thailand, the blood samples were collected from Trat and Chantaburi (TC, Eastern area), Kanchahaburi (KC, Western area), and from Nan, Chiang Rai and Lampang (NCL, Northern area). From Laos, the blood samples were collected from Viangchan and Louangphrabang (VC and LP, Central area) and from Pakxe (PK, Southern area). The sampled populations, abbreviations, and number of samples in each population are presented in Table 1. Genomic DNA was extracted from blood samples using standard phenol-chloroform extraction protocols (Sambrook, 1989). Detailed information and genotyping methods of the 102 indels markers used in this study (m1 to m102) can be found in our previous study (Maw et al., 2012). The genotypes of the indels polymorphisms were determined by size differences among the PCR fragments.

Statistical analysis

The genetic variability of each population was assessed by calculating the proportion of polymorphic loci (Ppoly: Lewontin and Hubby, 1966), and average observed heterozygosity (HO) and expected heterozygosity (HE) (Nei, 1978). A chi square (x2) approximation was used to test for Hardy-Weinberg equilibrium (HWE; Weir, 1996) and calculate the degree of genetic differentiation (Gst; Nei, 1973). Using Fstat version 2.9.3.2, pair-wise distance matrix (Fst) values (Weir and Cockerham, 1984) were estimated. A neighbor-joining (NJ) tree was constructed using MEGA 4 (Tamura et al., 2007), and we computed the pairwise standard genetic distance (Ds) (Nei, 1972) using PHYLIP ver. 3.69 (Felsenstein, 2009). Population structure was investigated using the Bayesian clustering algorithm implemented in the STRUCTURE 2.3.3 application (Pritchard et al., 2000). We applied a burn-in period of 20,000 followed by 20,000 iterations of the Markov chain Monte Carlo (MCMC) inference (Pritchard et al., 2000). The probable number of genetic clusters (K) was set to the range of 1≤K≤15, with 20 independent runs for each value of K. For the inference of the K value, we used the estimated likelihood of the probability of data Ln Pr (X/K) and the second order rate of change of the likelihood function with respect to KK). The ΔK value was calculated using the Structure Harvester 0.56.3 application (Riztyan et al., 2011). The optimal K value was determined by the greatest ΔK value as described in Evanno et al. (2005). The DISTRUCT 1.1 application (Rosenberg, 2004) was used to display the probabilities for population subdivisions. Subdivisions were represented as colors, and individuals were depicted as bars divided into color segments corresponding to the membership coefficients for the populations.

RESULTS

Indels polymorphisms and genetic variability

The mean (±standard error) Ppoly, HO, and HE values for each population are listed in Table 2. The native chickens from the three countries showed polymorphism in most of the indels loci. The highest (0.967±0.018) and the lowest (0.713±0.045) Ppoly were found in the KC and PG populations. The average HO ranged from 0.205±0.005 in the NCL population to 0.263±0.004 in the VC population. The average HE ranged from 0.239±0.004 in the PG population to 0.381±0.004 in the KC population. Tests of HWE for all cumulated populations showed HWE at p<0.05. In addition, each populations of Myanmar, Thailand, and Laos were in HWE (p<0.05). The degree of Gst was 0.041 in Myanmar, 0.088 in Thailand, and 0.024 in Laos (Table 2). The Gst among the nine populations was calculated as 0.125.

Pairwise Fst distance and phylogenic analysis

The pairwise Fst distances are shown in Table 3. The Fst distances between the native chicken populations ranged from 0.045 to 0.094 in Myanmar, from 0.116 to 0.154 in Thai, and from 0.018 to 0.037 in Laos. The pairwise Fst distances between Myanmar and Thailand native chicken populations ranged from 0.194 to 0.308; between Myanmar and Laotian, 0.144 to 0.249; and between Thai and Laotian, 0.178 to 0.292. The smallest distance (0.144) was between the YG and PK population, and the largest (0.308) was between the PG and NCL populations.
The NJ tree constructed from the pairwise standard Ds between the nine populations is shown in Figure 1. The populations from the three countries were assembled as respective clusters, and the tree is composed of two major clades. The first clade was composed of the Thai and Laotian native chicken populations. The second clade was composed of the Myanmar native chicken populations, which was located outside the first clade.

Population substructure by using structure analysis

The STRUCTURE analysis revealed that the Ln Pr (X/K) at K = 3 with the maximum likelihood value of –23,407.74. The maximum value for the ad hoc statistic, ΔK, was 171.107. Therefore, the probable number of clusters (K) was thought to be three. The contribution of the three ancestral populations (K = 3) to the nine populations is presented in Figure 2. The proportion of membership suggested that the native chicken populations from the three countries were descended from three different origins. The presence of admixture was observed among the Thai and Laotian native chicken populations, whereas the Myanmar native chickens lacked admixture.

DISCUSSION

The results revealed that all of the native chicken populations showed polymorphism in most of the indels loci (71% to 96%). This is larger than that observed in a previous study by Väli et al. (2008). Additionally, their study stated that 81 (86.2%) and 76 (80.9%) out of 94 indels markers could be validated as polymorphic loci in dogs (n = 7) and wolves (n = 18), respectively (Väli et al., 2008). The HE values of chickens native to Myanmar in the present study were similar to those obtained by SNPs in an earlier report by Riztyan et al. (2012). However, for Thai and Laotian native chicken populations, the HE values in the present study were higher than those obtained by SNPs and blood protein polymorphisms in earlier studies (Okamoto et al., 1999; Riztyan et al., 2012). The HE values of Myanmar native chickens were lower than those of the Thai and Laotian populations. Therefore, the genetic variability of the Thai and Laotian native chickens was higher than that of Myanmar populations. This may be linked to the practice of crossbreeding between native chickens and captured Red jungle fowl for cockfighting and as a source of game birds for sport in Thailand and the Philippines (Nishida et al., 2000). Free-scavenging chickens can easily reproduce with wild ones, which may lead to increased genetic variability. There was no HWE deviation among any of the populations or within any populations, suggesting that random mating is predominant among native chickens in those three countries.
The Gst value of Thai native chicken populations was 0.088, indicating that the degree of genetic differentiation is higher in Thai native chickens than Myanmar and Laotian native chickens. The Gst values in the present study were higher than those of the chickens native to Myanmar and Laos (0.024 and 0.025), and lower than those of the Thai native chickens (0.127), estimated from the genetic constitution of endogenous protease gene polymorphisms (Okumura et al., 2006). The Gst of Thai native chickens in the present study was also lower than that seen in the study using of 20 microsatellite loci (0.183) of Dorji et al. (2011). The potential source of these variations could be the sample populations, sample sizes, characteristics of different markers and experimental conditions. However, the Gst of the Myanmar and Laotian native chickens was close to the range of 0.001 to 0.039 (Kinoshita et al., 2002), estimated from egg white protein polymorphisms of local populations in Asian countries (Myanmar, Indonesia, China, Nepal, Vietnam, and Laos). The Gst for all the populations was 0.125, indicating that the genetic differentiation among the native chickens from the three countries was not large. Therefore, Myanmar, Thai, and Laotian native chickens can be regarded as genetically close populations.
The lower pairwise Fst distances observed within each country (0.018 to 0.154) and higher pairwise Fst distances observed between countries (0.144 to 0.308), reflect the fact that the populations within the countries were genetically closer populations than those between countries. The pairwise Fst distances between the Thai and Laotian populations were closer than those between the Myanmar and Thai populations. This is consistent with the topology of the NJ tree constructed from Ds distances. In the NJ tree, the populations of native chickens from Thailand and Laos were genetically close to one another, while the Myanmar native chickens were genetically distant from them. A similar genetic relationship was oserved by Okumura et al. (2006).
According to the STRUCTURE analysis, native chicken populations from each country were found to have distinct origins because of the long history of domestication in each country. Chickens became domesticated from Red jungle fowl in Southeast Asia nearly 10,000 years ago (Sawai et al., 2010). Interestingly, a small component of admixture from Myanmar populations was found in both Thai and Laotian populations. This was in agreement with the report of Riztyan et al. (2012). In addition, a small component of admixture from Laos was found in Thai populations. However, the admixture from Thai and Laotian populations was not found in Myanmar native chicken populations. The probable causes of admixture from Myanmar to Thai populations may be linked to the influence of Myanmar upon Thailand during the Ayutthaya era, 400 years ago. Myanmar native chickens seem to have been carried along by Myanmar’s military (Cocks, 1919; Hmannan Yazawin 2, 1829) and then spread to Laos from Thailand.
The present study examined the genetic diversity and population structure of native chickens from Myanmar, Thailand, and Laos using 102 indels polymorphisms as genetic markers. The native chickens from these three countries demonstrated polymorphisms in most indels loci and are genetically rich diversified populations. The genetic differentiation was higher in Thai native chickens than Myanmar and Laotian populations. According to the NJ tree, Thai and Laotian native chicken populations were genetically closer than Myanmar native chickens. The population structure of the nine chicken populations from the three countries was thought to be based on three different origins (K = 3). It was found that genetic admixture was observed in Thai and Laotian native chicken populations. The genetic composition of Myanmar native chickens is pure because of the lack of admixture from the other two countries.

Figure 1
Neighbor-joining (NJ) tree constructed by a genetic distance matrix from nine chicken populations from Myanmar, Thailand, and Laos. YG, Yangon; PG, Pegu; MD, Mandalay; TC, Trat and Chantaburi, KC, Kanchahaburi; NCL, Nan, Chiang Rai, Lampang; VC, Viangchan; LP, Louangphrabang; PK, Pakxe.
ajas-28-1-14f1.gif
Figure 2
Individual assignment according to the estimated membership coefficient at K = 3. YG, Yangon; MD, Mandalay; PG, Pegu; TC, Trat and Chantaburi, KC, Kanchahaburi; NCL, Nan, Chiang Rai, Lampang; VC, Viangchan; LP, Louangphrabang; PK, Pakxe.
ajas-28-1-14f2.gif
Table 1
Sampled populations, abbreviations and number of samples in each population
Country Population Abbreviation No. of Sample
Myanmar Yangon YG 27
Mandalay MD 40
Pegu PG 13
Thailand Trat and Chantaburi TC 15
Kanchahaburi KC 35
Nan, Chiang Rai, Lampang NCL 19
Laos Viangchan VC 50
Louangphrabang LP 44
Pakxe PK 35
Total 278
Table 2
The genetic variability of nine native chicken populations from Myanmar, Thailand, and Laos
Population PPoly ±SE O ±SE E ±SE Gst
YG 0.901 0.033 0.216 0.004 0.266 0.004 0.125
MD 0.871 0.033 0.218 0.005 0.239 0.004
PG 0.713 0.045 0.229 0.005 0.263 0.005
TC 0.935 0.026 0.221 0.005 0.369 0.004
KC 0.967 0.018 0.212 0.004 0.381 0.004
NCL 0.946 0.024 0.205 0.005 0.340 0.004
VC 0.931 0.025 0.263 0.004 0.317 0.004
LP 0.902 0.030 0.244 0.004 0.288 0.004
PK 0.882 0.032 0.261 0.004 0.301 0.004
Myanmar native chickens 0.041
Thailand native chickens 0.088
Laos native chickens 0.024

YG, Yangon; MD, Mandalay; PG, Pegu; TC, Trat and Chantaburi, KC, Kanchahaburi; NCL, Nan, Chiang Rai, Lampang; VC, Viangchan; LP, Louangphrabang; PK, Pakxe.

Table 3
Pairwise Fst distances among nine native chicken populations from Myanmar, Thailand, and Laos
YG MD PG TC KC NCL VC LP PK
YG 0.045 0.094 0.194 0.239 0.286 0.162 0.189 0.144
MD 0.073 0.202 0.244 0.291 0.168 0.196 0.165
PG 0.246 0.243 0.308 0.205 0.249 0.196
TC 0.127 0.154 0.178 0.218 0.201
KC 0.116 0.213 0.252 0.228
NCL 0.260 0.292 0.284
VC 0.028 0.018
LP 0.037
PK

YG, Yangon; MD, Mandalay; PG, Pegu; TC, Trat and Chantaburi, KC, Kanchahaburi; NCL, Nan, Chiang Rai, Lampang; VC, Viangchan; LP, Louangphrabang; PK, Pakxe.

REFERENCES

Cocks SW. 1919. A short history of Burma. 2nd editionMacmillan and Co. Limited; London, UK:

Dorji N, Daungjinda M, Phasuk Y. 2011. Genetic characterization of Thai indigenous chickens compared with commercial lines. Trop Anim Health Prod 43:779–785.
crossref pmid
Evanno G, Regnaut S, Goudet J. 2005. Detecting the number of clusters of individuals using the software STRUCTURE: A simulation study. Mol Ecol 14:2611–2620.
crossref pmid
Felsenstein J. 2009. Phylogeny Inference Package Version 3.69.

Hmannan Yazawin 2. 1829. Tabinshwehti’s Ayutthaya campaign (in Burmese). 2003 editionMinistry of Information; Yangon, Myanmar:

Kawabe K, Worawut R, Taura S, Shimogiri T, Nishida T, Okamoto S. 2014. Genetic diversity of mtDNA D-loop polymorphisms in Laotian native fowl populations. Asian Australas J Anim Sci 27:19–23.
crossref pmid pmc
Kinoshita K, Okamoto S, Shimogiri T, Kawabe K, Nishida T, Kakizawa R, Yamamoto Y, Maeda Y. 2002. Gene constitution of egg white proteins of native chicken in Asian countries. Asian Australas J Anim Sci 15:157–165.
crossref
Lewontin RC, Hubby JL. 1966. A molecular approach to the study of genetic heterozygosity in natural populations of Drosophila pseduobscura. Genetics 54:595–609.
pmid pmc
Liu YP, Wu GS, Yao YG, Miao YW, Luikart G, Baig M, Beja-Pereira A, Ding ZL, Palanichamy MG, Zhang YP. 2006. Multiple maternal origins of chickens: Out of the Asian jungles. Mol Phylogenet Evol 38:12–19.
crossref pmid
Maw AA, Shimogiri T, Riztyan , Kawabe K, Kawamoto Y, Okamoto S. 2012. Genetic diversity of Myanmar and Indonesia native chickens together with two jungle fowl species by using 102 Indels polymorphisms. Asian Australas J Anim Sci. 25:p. 927–934.
crossref pmid pmc
Mwacharo JM, Jianlin H, Amano T. 2006. Native African chicken: valuable genetic resources for future breeding improvement. J Anim Genet Jpn 34:63–69.
crossref
Nishibori M, Shimogiri T, Hayashi T, Yasue H. 2005. Molecular evidence of hybridization of species in the genus Gallus except of Gallus varius. Anim Genet 36:367–375.
crossref pmid
Natalle SCF, Resque RL, Ribeiro-Rodrigues EM, Guerreiro JF, Santos NPC, Santos AR, Santos S. 2010. X-linked insertion/deletion polymorphisms: Forensic applications of a 33-markers panel. Int J Legal Med 124:589–593.
crossref pmid
Nei M. 1973. Analysis of gene diversity in subdivided populations. Proc Natl Acad Sci USA 70:3321–3323.
crossref pmid pmc
Nei M. 1978. Estimation of average heterozygosity and genetic distance from a small number of individuals. Genetics 89:583–590.
crossref pmid pmc pdf
Nishida T, Rerkamnuaychoke W, Tung DG, Saignaleus S, Okamoto S, Kawamoto Y, Kimura J, Kawabe K, Tsunekawa N, Otaka H, Hatashi Y. 2000. Morphological identification and ecology of Red Jungle Fowl in Thailand, Laos and Vietnam. Anim Sci J 71:470–480.
crossref
Okamoto S, Tsunekawa N, Kawamoto Y, Worawut R, Kawabe K, Maeda Y, Nishida T. 1999. Blood protein polymorphisms of native fowls in Laos. Asian Australas J Anim Sci 12:1011–1014.
crossref
Okamoto S, Kinoshita K, Yamamoto Y, Tsunoda K, Okabayashi H, Mannen H, Tanaka K, Yamagata T, Nozawa K, Nishibori M, Kurosawa Y, Maung Maung, Than Nyunt, Than Daing, Nay Win Hla, Maeda Y. 2004. Blood protein variations of native fowls in Myanmar. Researches Report of the Society for Native Livestock 21:237–242.

Okumura F, Shimogiri T, Kawabe K, Okamoto S, Nishibori M, Yamamoto Y, Maeda Y. 2006. Gene constitution of South-East Asian native Chickens, commercial chickens and jungle fowl using polymorphisms of four Calpain genes. Anim Sci J 77:188–195.
crossref
Pritchard JK, Stephens M, Donnelly P. 2000. Inference of population structure using multilocus genotype data. Genetics 155:945–959.
crossref pmid pmc pdf
Riztyan , Katano T, Shimogiri T, Kawabe K, Okamoto S. 2011. Genetic diversity and population structure of Indonesian native chickens based on single nucleotide polymorphism makers. Poult Sci 90:2471–2478.
crossref pmid
Riztyan , Nishida T, Worawut R, Shimogiri T, Maw AA, Kawabe K, Kawamoto Y, Okamoto S. 2012. Genetic variation and phylogeographic analysis of native chicken populations in Myanmar and Thailand. J Poult Sci 49:68–73.
crossref
Rosenberg NA. 2004. DISTRUCT: a program for the graphical display of population structure. Mol. Ecol Notes 4:137–138.
crossref
Sawai H, Kim HL, Kuno K, Suzuki S, Gotoh H, Takada M, Takahata N, Satta Y, Akishinonomiya F. 2010. The origin and genetic variation of domestic chickens with special reference to jungle fowls Gallus g. gallus and G. varius. PLoS One 5:e10639
crossref pmid pmc
Sambrook J, Fritsch EF, Maniatis T. 1989. Molecular cloning: A laboratory manual. 2nd editionCold Spring Laboratory Press; Cold Spring Harbor, New York, USA:

Tamura K, Dudley J, Nei M, Kumar S. 2007. MEGA4: Molecular evolutionary genetics analysis (MEGA) software version 4.0. Mol Biol Evol 24:1596–1599.
crossref pmid
Väli U, Brandström M, Johansson M, Ellegren H. 2008. Insertion-deletion polymorphisms (indels) as genetic markers in natural populations. BMC Genet 9:8
crossref pmid pmc
Weir BS. 1996. Genetic Data Analysis II. Sinauer; Sunderland, MA, USA:

Weir BS, Cockerham CC. 1984. Estimating F-statistic for the analysis of population structure. Evolution 38:1358–1370.
crossref pmid pdf


Editorial Office
Asian-Australasian Association of Animal Production Societies(AAAP)
Room 708 Sammo Sporex, 23, Sillim-ro 59-gil, Gwanak-gu, Seoul 08776, Korea   
TEL : +82-2-888-6558    FAX : +82-2-888-6559   
E-mail : editor@animbiosci.org               

Copyright © 2024 by Asian-Australasian Association of Animal Production Societies.

Developed in M2PI

Close layer
prev next