- Research article
- Open Access
Genetic evidence from Indian red jungle fowl corroborates multiple domestication of modern day chicken
© Kanginakudru et al; licensee BioMed Central Ltd. 2008
- Received: 19 October 2006
- Accepted: 10 June 2008
- Published: 10 June 2008
Domestication of chicken is believed to have occurred in Southeast Asia, especially in Indus valley. However, non-inclusion of Indian red jungle fowl (RJF), Gallus gallus murghi in previous studies has left a big gap in understanding the relationship of this major group of birds. In the present study, we addressed this issue by analyzing 76 Indian birds that included 56 G. g. murghi (RJF), 16 G. g. domesticus (domestic chicken) and 4 G. sonneratii (Grey JF) using both microsatellite markers and mitochondrial D-loop sequences. We also compared the D-loop sequences of Indian birds with those of 779 birds obtained from GenBank.
Microsatellite marker analyses of Indian birds indicated an average FST of 0.126 within G. g. murghi, and 0.154 within G. g. domesticus while it was more than 0.2 between the two groups. The microsatellite-based phylogenetic trees showed a clear separation of G. g. domesticus from G. g. murghi, and G. sonneratii. Mitochondrial DNA based mismatch distribution analyses showed a lower Harpending's raggedness index in both G. g. murghi (0.001515) and in Indian G. g. domesticus (0.0149) birds indicating population expansion. When meta analysis of global populations of 855 birds was carried out using median joining haplotype network, 43 Indian birds of G. g. domesticus (19 haplotypes) were distributed throughout the network sharing haplotypes with the RJFs of different origins.
Our results suggest that the domestication of chicken has occurred independently in different locations of Asia including India. We found evidence for domestication of Indian birds from G. g. spadiceus and G. g. gallus as well as from G. g. murghi, corroborating multiple domestication of Indian and other domestic chicken. In contrast to the commonly held view that RJF and domestic birds hybridize in nature, the present study shows that G. g. murghi is relatively pure. Further, the study also suggested that the chicken populations have undergone population expansion, especially in the Indus valley.
- Mismatch Distribution
- Domestic Chicken
- Domestication Event
- Domestic Bird
- Chicken Population
Archeological findings have indicated that the 'mother of all poultry' is the Southeast (SE) Asian Red jungle fowl (RJF) (Gallus gallus). Since domestication of chicken has been observed at the Indus valley as early as 3,200 BC, it is believed to be the epicenter of chicken domestication . However, later day excavations in Peiligan Neolithic sites of China have raised questions about the exclusive domestication at Indus valley, suggesting alternate and possibly earlier domestication centers . It is proposed that G. gallus, the wild RJF found in the forests of SE Asia and India, spread to other parts of the world when people domesticated the chicken, resulting in many chicken breeds [3, 4]. Subsequent to domestication, the extensive breeding programmes have resulted in sixty or so breeds of chicken representing four distinct lineages: egg-type, game, meat-type and bantam . While some authors suggest monophyletic origin of domestic chicken [6, 7], others provide evidence for multiple and independent domestication events . Such inconsistent observations are attributable to the fact that the initial studies were done with relatively small set of samples. In all these reported studies the native RJFs of Indian sub-continent, G. g. murghi were not represented in the analyses due to lack of sequence or molecular marker information on this group of birds.
Taxonomically, genus Gallus is composed of four species, G. gallus (RJF), G. lafayettei (Lafayette's JF), G. varius (Green JF) and G. sonneratii (Grey JF – GJF). Presently there are 5 sub-species of RJF, G. g. gallus (SE Asian RJF), and G. g. spadiceus, G. g. bankiva, G. g. murghi (Indian RJF) and G. g. jabouillei . These classifications are mainly based on phenotypic traits and geographic distribution of the populations. In literature, wild and domesticated birds are often referred to as 'fowls' and 'chicken', respectively. The domestic chicken is considered either as a sub-species of RJF (G. g. domesticus) or as a separate species, G. domesticus. However, tight clustering of the different sub-species discounts this existing taxonomical hierarchy  rendering sub-species status within RJF redundant.
Besides the taxonomical intricacies, the researchers are also concerned about the genetic integrity and conservation status of the RJF in the wild and those held in avicultural collections. It is suspected that the domestic chicken is hybridizing with the wild RJF resulting in erosion of genetic purity of the wild birds [4, 10, 11]. Most of these earlier studies are based on either phenotypic characters or DNA analyses confined to small samples. Phylogenetic analyses of mitochondrial D-loop sequence and nuclear genes have indicated possible hybridization between GJF-RJF/domestic birds . In the light of these reports it is important to assess the genetic uniqueness of Indian RJFs not only from conservation point of view, but also for using them in population studies.
It is well known that the patterns of genetic differences can reveal the demographic history of the population under study. Since changes in population size leave characteristic molecular signatures, by measuring such changes one can reconstruct the population history. Mismatch distributions, also known as 'pairwise differences' give information on genetic differences between pairs of subjects and can be used to understand the population history . Mitochondrial DNA (mtDNA) is inherited maternally  and is often employed in population genetic analyses due to its high copy number, haploid nature and absence of/rare recombination events . A recent study involving the analysis of chicken mitochondrial DNA sequence from pre-historical samples pointed pre-columbian Polynesian origin of American breeds , thus suggesting the importance of mitochondrial D-loop sequence in determining the history of chicken domestication. Microsatellite markers, on the other hand, are nuclear markers and are used extensively in population genetic analyses because they are highly polymorphic, ubiquitously distributed throughout the genome, are having high mutation rates, co-dominant in nature, selectively neutral and are amenable to PCR-based high through-put analysis . Hilel et al.  characterized 52 chicken breeds using 22 microsatellite markers and concluded that the origin of domestic chicken to be from RJFs, as supported by mt DNA analyses .
In all the earlier studies these two marker systems have been used independently to study chicken populations [7, 17–19]. In the present study we have combined the results emanating from these two informative marker systems to address the questions relating to (i) evolutionary status of Indian RJF and chicken and (ii) extent of gene flow between Indian RJF and chicken, RJF-GJF in comparison to the world population.
Genetic identity of Indian fowls
In the present study, phylogenetic and demographic profiling analyses were carried out using 76 Indian birds that belong to seven populations. We employed 11 microsatellite markers and also sequenced 650 bp of hyper variable region (D-loop) of mitochondrial genomes of two species of fowls namely G. sonneratii (n = 4) and G. gallus, which includes two subspecies G. g. murghi, (n = 56) and G. g. domesticus (n = 16).
Microsatellite loci and PCR conditions used in the present study, and population-wise number of alleles obtained for Indian chicken. Original references are listed in the table.
Allele Size range
Total no. of alleles*
G. g. murghi
G. g. domesticus
Population pairwise FST values of microsatellite (below diagonal) and mitochondrial D-loop sequence (above diagonal) for 76 Indian chicken classified into 7 groups.
Population-wise Nei's genetic distance calculated using GenAlEx program showed a higher average distance between G. sonneratii and G. g. domesticus (2.099) than between G. sonneratii – G. g. murghi (0.758) and G. g. murghi – G. g. domesticus (1.695) combinations (Additional file 1). Intra-population average distances were lower for both domestic and RJF groups than across the population distances, which is consistent with the observation that among population variation was more than within population variation. The results suggest very rare genetic exchange between the RJF and domestic chicken populations, at least in recent history.
We also constructed a genetic distance based neighbor-joining (NJ) tree to obtain the genetic relationship among Indian birds. The result clearly points to the fact that hybridization between Indian chicken and Indian RJF G. g. murghi in the wild is extremely rare (Fig. 1b).
We also sequenced and analyzed 650 bp of the D-loop region of 76 Indian birds to derive the matrilineal population history by using coalescent-based models.
The pairwise FST values calculated using Arlequin were very low within G. g. murghi when compared to G. g. domesticus (Table 2). The average FST values were more for G. sonneratii – G. g. domesticus combination than for G. sonneratii-G. g. murghi or G. g. murghi – G. g. domesticus combinations. Contrary to microsatellite based analysis, mitochondrial analysis showed a lower mean pairwise FST value for G. sonneratii – G. g. domesticus combination. These results again point out that there is hardly any genetic exchange across the three categories of birds, namely GJF, RJF and domestic birds.
Phylogeny and median-joining network profiles of Indian birds in comparison to global populations
Standard diversity indices and mismatch distribution analyses obtained from d-loop sequence of 855 birds belonging to 13 populations. Coturnix japonica that was used as an outgroup is not included.
G. gallus domesticus (Domestic chicken)
G. gallus (RJF)
Other species of Gallus
Haplotype diversity (Hd)
Sum of square freqs.
Number of observed transitions
Number of observed transversions
Number of substitutions
Number of polymorphic sites
Nucleotide diversity Standard deviation (±)
Mismatch observed mean
Harpending's Raggedness index
P (Sim. Rag. > = Obs. Rag.)
P(D random < D obs)
Prob(sim_Fs < = obs_Fs)
The sequenced regions of D-loop were well conserved in most of the birds. The 62 bp insertion element found in G. sonneratii was the only indel of considerable size present in the whole sequence. The salient mutation observed in G. g. murghi was the presence of 'T' in most of the birds (54 out of 56 i.e. 96.4%) (Fig. 4) that was present in only 33.3% of G. g. gallus, 24.4% of G. g. spadiceus and 13.1% of all the G. g. domesticus birds. Only three Indian G. g. domesticus (6.9%) had 'T' at this position. This position corresponds to the nucleotide number 360 (where there is 'C') in the complete mitochondrial genome sequence of G. gallus (Acc. No. NC_001323) available in GenBank. Across the portion of the sequenced mitochondrial genome we observed that transition to transversion ratio was lower in case of G. g. murghi (45/47 = 0.99) than in G. g. gallus (31/2 = 15.5) or G. g. spadiceus (23/4 = 5.75 – Table 3).
Mismatch distribution analyses carried out under sudden expansion model showed that mean pairwise differences were highest in G. g. murghi among G. gallus sub-species. However, another species of Gallus, G. sonneratii had the highest mean pairwise differences (Table 3). Both Tajima's D and Fu's Fs were significantly (P < 0.05) negative in G. g. murghi, suggesting the departure from neutrality (Table 3). Tajima's D was significant only in G. g. spadiceus and G. g. murghi. Fu's FS, a better indicator for estimating the departures from neutral theory, showed a significant negative value (P of simFS < = obsFS is < 0.01) for G. g. murghi and G. g. domesticus (Indonesia). Nucleotide diversity was also high (0.030) in case of G. g. murghi (except for G. sonneratii) amongst all the groups studied (Table 3).
The FST values showed a very high differentiation between the out group C. japonica and the other 13 populations, with a value above 0.9. Most of the pairwise FST values were significant with a P value < 0.05. FST value  based NJ tree showed the divergence of G. g. murghi from G. g. gallus, G. g. spadiceus and G. g. domesticus (Additional file 3). AMOVA calculations showed that the majority of variation found within G. gallus subspecies was between the domestic and jungle fowl populations (76%), while 'among the group' variation was 7% with an overall FST value of 0.234 (P = 0).
In the present study we attempted to understand the contribution of Indian red jungle fowl, G. g. murghi to the domestication event. As of now, no sequence information is available from this group of birds which, in all likelihood, was contributor of one of the earliest known chicken domestication event, i.e. in Mohanjo-Daro. We studied Indian birds belonging to two species of Gallus and compared them with the worldwide bird populations. Since microsatellite markers and the D-loop sequence of mitochondrial DNA have a high mutation rates, they provide information about recent evolutionary history as compared to slow mutating genes that provide data about ancient history . To reconstruct the recent past, we used both these marker systems and also addressed the issue of genetic purity of wild birds.
From the mtDNA analysis, we observed that a group of G. g. domesticus birds had the G. g. murghi haplotype (H_6), while a few others shared haplotypes with G. g. gallus (H_21, H_59, H_61 etc.) and with G. g. spadiceus (H_39, H_42, H_58 etc.). This is also true for Indian chicken that have originated by independent domestication from G. g. murghi as well as possibly from other G. g. subspecies. Interestingly, sharing of different haplotypes by Indian domestic chicken clarify that the present day Asiatic chicken might have originated from different progenitors by multiple domestication events and such multi-origin breeds could still be observed in a single geographical location. This is consistent with the observation of Oka et al. , who showed that the present day native Japanese chicken are having multiple origin.
The multilocus microsatellite data as well as the D-loop sequence of Indian chicken showed departure from neutrality as indicated by significant negative value of Tajima's D and Fu's Fs for G. g. murghi suggesting the possible population expansion of Indian birds. These results are also consistent with mismatch distribution analyses and significant value for Harpending's raggedness index (Table 3). Taken together with the genetic diversity analyses, we surmise that G. g. murghi and G. g. domesticus (India) must have undergone population expansion. It is believed that population expansion follows a domestication event. This fact holds true even in case of chicken, where we did observe population expansion, as did previous studies.
Our analyses revealed that the sampled Indian birds are relatively pure with very rare hybridization between G. g. murghi and G. g. domesticus (India). Nishibori et al.  suggested the hybridization of RJF with domestic and grey jungle fowl. In the present study of Indian birds, we did not come across noticeable hybridization at least in the recent past, as indicated by very low FST values for mtDNA (Table 2) and microsatellite markers and also a clear separation of RJF clades from domestic chicken in microsatellite based phylogeny. All these results indicate the genetic integrity of the G. g. murghi.
In the present study, we observed predominant occurrence of a characteristic 'T' nucleotide in 96.4% of G. g. murghi birds that is absent in most of the Indian domestic chicken further supporting the occurrence of negligible hybridization between G. g. murghi and G. g. domesticus (India). All the sampled G. sonneratii had 'A' nucleotide in this position. If the frequent hybridization is occurring then it is expected that at least a few of G. g. murghi to have 'A' at this position. Since we found only one RJF with 'A' at this site, and also due to clear separation of the clades, it is unlikely that G. g. murghi and G. sonneratii hybridize in the wild contrary to the observations made by Nishibori et al.  who suggested the possible hybridization between the RJF and GJF. Such contrasting observations may be due to the limited number of samples (3 RJF and 3 GF) used by Nisibori et al. . Taken together with the Fumihito et. al.'s [6, 7] and Liu et al.'s  observations our results prompted us to question the sub-species status for G. g. gallus, G. g. spadiceus, G. g. murghi and G. g. domesticus. In the light of these findings, we recommend that G. gallus should be classified as G. g. gallus that should include all RJFs and G. g. domesticus birds. At the same time, after confirmation of reproductive isolation, G. g. bankiva could be placed into a separate species, Gallus bankiva.
For the first time by analyzing hitherto unreported samples of G. g. murghi and also including the G. g. domesticus (India), we confirm that the domestication of chicken has occurred independently from G. g. murghi. We also provide evidence that there is little genetic exchange between G. g. murghi and G. g. domesticus (India) and minimal hybridization between G. sonneratii and G. g. murghi. Comparison of Indian RJF and domestic bird to that of world population also supports the previous studies of obsoleteness of the sub-species status given to RJFs and domestic chicken.
Samples, DNA isolation
About 0.5 ml of blood was collected from the wing vein of live birds into vials containing 5 mM EDTA and genomic DNA was isolated as per standard protocols .
PCR amplification of microsatellite loci
The microsatellite loci used in the study, the primer details and PCR conditions are available in Table 1. The PCR products were separated on 3.5% Metaphor agarose gel along with pUC/MspI digest (MBI Fermentas) ladder. For GeneScan analysis, the PCR products generated using fluorescent dUTPs were dissolved in 2 μl of formamide gel loading buffer with 0.3 μl of ROX-500™ GeneScan ruler (Perkin Elmer) and separated on a 5% polyacrylamide-7 M Urea gel.
PCR amplification and Sequencing of D-loop of mitochondrial DNA
The D-loop hypervariable region was PCR amplified using primers described elsewhere . The sequencing was carried out on both the strands. We also used an internal primer (5'GTGGAATATAGGTTAATGCC 3') to obtain the sequence information from 5' region without any ambiguity. 50 ng PCR product was used in a sequencing reaction that contained 8 μl of Ready reaction mix (BDT v 3.0, Applied Biosystems, Foster City, CA) and 5 picomoles of primer. The sequencing was carried out in ABI Prism 3100 Genetic Analyzer (Applied Biosystems).
The individuals were genotyped based on allele size data. Allele frequency and heterozygosity were calculated using GenAlEx . F-statistics were used as a measure of diversity within and between populations and were estimated using GenAlEx. For PCA, genetic distance was calculated from the allele data and the genetic distance was plotted as PCA using GenAlEx. Population history parameters were calculated using Arlequin . The microsatellite allele frequency data from different populations was bootstrapped using seqboot of Phylip and the output file was used for construction of maximum likelihood (ML) tree using Phylip 'contml' program. For distance based NJ tree, the genetic distance was calculated using GenAlEx program and the resulting distance matrix was used to construct NJ tree with Kimura-2-parameter (K2P) option in MEGA.
Mitochondrial D-loop sequences of 779 Jungle Fowls and G. g. domesticus birds as well as other species of Gallus were obtained from GenBank  (available as Additional file 4). The sequences of 855 birds, including the 76 samples sequenced in the present study were aligned using ClustalX program , manually edited using GeneDoc  and the region conserved in all the birds was used for analysis. With gaps there were a total of 482 bp sequence, with highest number of nucleotides (460–462) coming from the G. sonneratii and approximately 400 bp from RJFs. Similar to previous phylogenetic studies, we used Coturnix japonica as an outgroup in our study. The phylogenetic trees were constructed using Phylip 3.5  or MEGA . The best fit model was selected using the program MODELTEST  as implemented in HyPhy . The alpha value obtained from best-fit model was used for gamma correction in haplotype NJ tree, which was constructed using Tajima and Nei's model in MEGA3. Haplotype data was obtained using DnaSP  and was used to construct haplotype network using Network program . Default parameters were used for obtaining the median joining network tree. Population genetic structure was measured using AMOVA as implemented in Arlequin with 1000 permutations. Tajima's D , Fu's Fs  and other population genetics parameters were also calculated using Arlequin for which the significance was tested after 1000 simulation steps . For such analyses, populations were defined depending on the data being analyzed – e.g. in case of Indian chicken, as 2G – 2 populations viz. G. g. domesticus and G. g. murghi, 6G-3 sub populations each based on the sampling location in both G. g. domesticus (C, J, M) and G. g. murghi populations (M-RJF, B-RJF, K-RJF) or 7G (6G and G. sonneratii). Such a classification was carried out to study the population genetic structure within the subdivided populations. Whenever necessary, as in case of subdivided populations, realignment of the sequence was carried out using clustalX.
We investigated the demographic profiles of chicken populations based on coalescence theory and analyzed pairwise mismatch distribution to confirm the population expansion , using Arlequin. The parameter of demographic expansion τ was estimated with a generalized nonlinear least squares approach and approximate confidence intervals were obtained with 1000 parametric bootstrap replicates. The goodness-of-fit of the observed data to a simulated model of expansion was tested with the sum of squared deviations and Harpending's raggedness index was estimated . For bootstrap phylogenetic NJ tree, the aligned mtDNA sequence was run in MEGA with 1000 replicates with a 50% cutoff option.
We acknowledge the financial support received from the Department of Forests, Govt. of Haryana, India. We thank Mr. JPL Srivastava, Principal Chief Conservator of Forests, Govt. of Haryana for his support throughout the research.
- Zeuner FE: A History of Domesticated Animals. 1963, Harper and Row, New YorkGoogle Scholar
- West B, Zhou B-X: Did chicken go north? New evidence for domestication. J Archaeol Sci. 1988, 15: 515-533. 10.1016/0305-4403(88)90080-5.View ArticleGoogle Scholar
- Stevens L: Genetic and evolution of the domestic fowl. 1991, New York , Cambridge University PressView ArticleGoogle Scholar
- Peterson AT, Brisbin IL: Genetic endangerment of wild red jungle fowl (Gallus gallus)?. Bird Conservation International. 1999, 9: 387-394.View ArticleGoogle Scholar
- Moiseyeva IG, Romanov MN, Nikiforov AA, Sevastyanova AA, Semyenova SK: Evolutionary relationships of Red Jungle Fowl and chicken breeds. Genet Sel Evol. 2003, 35 (4): 403-423. 10.1051/gse:2003031.PubMed CentralView ArticlePubMedGoogle Scholar
- Fumihito A, Miyake T, Takada M, Shingu R, Endo T, Gojobori T, Kondo N, Ohno S: Monophyletic origin and unique dispersal patterns of domestic fowls. Proc Natl Acad Sci U S A. 1996, 93 (13): 6792-6795. 10.1073/pnas.93.13.6792.PubMed CentralView ArticlePubMedGoogle Scholar
- Fumihito A, Miyake T, Sumi S, Takada M, Ohno S, Kondo N: One subspecies of the red junglefowl (Gallus gallus gallus) suffices as the matriarchic ancestor of all domestic breeds. Proc Natl Acad Sci U S A. 1994, 91 (26): 12505-12509. 10.1073/pnas.91.26.12505.PubMed CentralView ArticlePubMedGoogle Scholar
- Liu YP, Wu GS, Yao YG, Miao YW, Luikart G, Baig M, Beja-Pereira A, Ding ZL, Palanichamy MG, Zhang YP: Multiple maternal origins of chickens: out of the Asian jungles. Mol Phylogenet Evol. 2006, 38 (1): 12-19. 10.1016/j.ympev.2005.09.014.View ArticlePubMedGoogle Scholar
- Niu D, Fu Y, Luo J, Ruan H, Yu XP, Chen G, Zhang YP: The origin and genetic diversity of Chinese native chicken breeds. Biochem Genet. 2002, 40 (5-6): 163-174. 10.1023/A:1015832108669.View ArticlePubMedGoogle Scholar
- Brisbin IL, Peterson AT, Okimoto R, Amato RG: Characterization of the genetic status of populations of red jungle fowl. J Bombay Natural History Society. 2002, 99: 217-223.Google Scholar
- Nishibori M, Shimogiri T, Hayashi T, Yasue H: Molecular evidence for hybridization of species in the genus Gallus except for Gallus varius. Anim Genet. 2005, 36 (5): 367-375. 10.1111/j.1365-2052.2005.01318.x.View ArticlePubMedGoogle Scholar
- Harpending HC, Batzer MA, Gurven M, Jorde LB, Rogers AR, Sherry ST: Genetic traces of ancient demography. Proc Natl Acad Sci U S A. 1998, 95 (4): 1961-1967. 10.1073/pnas.95.4.1961.PubMed CentralView ArticlePubMedGoogle Scholar
- Berlin S, Ellegren H: Evolutionary genetics. Clonal inheritance of avian mitochondrial DNA. Nature. 2001, 413 (6851): 37-38. 10.1038/35092623.View ArticlePubMedGoogle Scholar
- Savolainen P, Zhang YP, Luo J, Lundeberg J, Leitner T: Genetic evidence for an East Asian origin of domestic dogs. Science. 2002, 298 (5598): 1610-1613. 10.1126/science.1073906.View ArticlePubMedGoogle Scholar
- Storey AA, Ramirez JM, Quiroz D, Burley DV, Addison DJ, Walter R, Anderson AJ, Hunt TL, Athens JS, Huynen L, Matisoo-Smith EA: Radiocarbon and DNA evidence for a pre-Columbian introduction of Polynesian chickens to Chile. Proc Natl Acad Sci U S A. 2007, 104 (25): 10335-10339. 10.1073/pnas.0703993104.PubMed CentralView ArticlePubMedGoogle Scholar
- Schlotterer C: The evolution of molecular markers--just a matter of fashion?. Nat Rev Genet. 2004, 5 (1): 63-69. 10.1038/nrg1249.View ArticlePubMedGoogle Scholar
- Hillel J, Groenen MA, Tixier-Boichard M, Korol AB, David L, Kirzhner VM, Burke T, Barre-Dirie A, Crooijmans RP, Elo K, Feldman MW, Freidlin PJ, Maki-Tanila A, Oortwijn M, Thomson P, Vignal A, Wimmers K, Weigend S: Biodiversity of 52 chicken populations assessed by microsatellite typing of DNA pools. Genet Sel Evol. 2003, 35 (5): 533-557. 10.1051/gse:2003038.PubMed CentralView ArticlePubMedGoogle Scholar
- Crooijmans RP, van Kampen AJ, van der Poel JJ, Groenen MA: Highly polymorphic microsatellite markers in poultry. Anim Genet. 1993, 24 (6): 441-443.View ArticlePubMedGoogle Scholar
- Fumihito A, Miyake T, Takada M, Ohno S, Kondo N: The genetic link between the Chinese bamboo partridge (Bambusicola thoracica) and the chicken and junglefowls of the genus Gallus. Proc Natl Acad Sci U S A. 1995, 92 (24): 11053-11056. 10.1073/pnas.92.24.11053.PubMed CentralView ArticlePubMedGoogle Scholar
- Hudson RR, Slatkin M, Maddison WP: Estimation of levels of gene flow from DNA sequence data. Genetics. 1992, 132 (2): 583-589.PubMed CentralPubMedGoogle Scholar
- Jorde LB, Bamshad M, Rogers AR: Using mitochondrial and nuclear DNA markers to reconstruct human evolution. Bioessays. 1998, 20 (2): 126-136. 10.1002/(SICI)1521-1878(199802)20:2<126::AID-BIES5>3.0.CO;2-R.View ArticlePubMedGoogle Scholar
- Oka T, Ino Y, Nomura K, Kawashima S, Kuwayama T, Hanada H, Amano T, Takada M, Takahata N, Hayashi Y, Akishinonomiya F: Analysis of mtDNA sequences shows Japanese native chickens have multiple origins. Anim Genet. 2007, 38 (3): 287-293. 10.1111/j.1365-2052.2007.01604.x.View ArticlePubMedGoogle Scholar
- Loftus RT, MacHugh DE, Bradley DG, Sharp PM, Cunningham P: Evidence for two independent domestications of cattle. Proc Natl Acad Sci U S A. 1994, 91 (7): 2757-2761. 10.1073/pnas.91.7.2757.PubMed CentralView ArticlePubMedGoogle Scholar
- Giuffra E, Kijas JM, Amarger V, Carlborg O, Jeon JT, Andersson L: The origin of the domestic pig: independent domestication and subsequent introgression. Genetics. 2000, 154 (4): 1785-1791.PubMed CentralPubMedGoogle Scholar
- Sambrook J, W. RD: Molecular Cloning, A laboratory manual. 2001, New York. , Cold Spring Harbor Laboratory Press.Google Scholar
- Peakall R, Smouse PE: GENALEX 6: genetic analysis in Excel. Population genetic software for teaching and research. Molecular Ecology Notes. 2006, 6: 288-295. 10.1111/j.1471-8286.2005.01155.x.View ArticleGoogle Scholar
- Schneider S D. Roessli, L. Excoffier.: A software for population genetics data analysis. 2000, Geneva , Genetics and Biometry Lab, Dept. of Anthropology, University of Geneva, Ver. 2.000.Google Scholar
- GenBank. [http://www.ncbi.nlm.nih.gov]
- Thompson JD, Gibson TJ, Plewniak F, Jeanmougin F, Higgins DG: The CLUSTAL_X windows interface: flexible strategies for multiple sequence alignment aided by quality analysis tools. Nucleic Acids Res. 1997, 25 (24): 4876-4882. 10.1093/nar/25.24.4876.PubMed CentralView ArticlePubMedGoogle Scholar
- Nicholas KB, Nicholas HBJ: Genedoc, a Tool for Editing and Annotating Multiple Sequence Alignments. 1997, Distributed by the authorGoogle Scholar
- Felsenstein J: PHYLIP. Phylogeny interface package (version 3.2). Cladistics. 1989, 5: 164–166-Google Scholar
- Kumar S, Tamura K, Nei M: MEGA3: Integrated software for Molecular Evolutionary Genetics Analysis and sequence alignment. Brief Bioinform. 2004, 5 (2): 150-163. 10.1093/bib/5.2.150.View ArticlePubMedGoogle Scholar
- Posada D, Crandall KA: Modeltest: testing the model of DNA substitution. Bioinformatics. 1998, 14: 817-818. 10.1093/bioinformatics/14.9.817.View ArticlePubMedGoogle Scholar
- Pond SLK, Frost SDW, Muse SV: HyPhy: hypothesis testing using phylogenies Bioinformatics. Bioinformatics. 2005, 21: 676-679. 10.1093/bioinformatics/bti079.View ArticlePubMedGoogle Scholar
- Rozas J, Sanchez-DelBarrio JC, Messeguer X, Rozas R: DnaSP, DNA polymorphism analyses by the coalescent and other methods. Bioinformatics. 2003, 19 (18): 2496-2497. 10.1093/bioinformatics/btg359.View ArticlePubMedGoogle Scholar
- Network program. [http://www.fluxus-technology.com/]
- Tajima F: Statistical method for testing the neutral mutation hypothesis by DNA polymorphism. Genetics. 1989, 123 (3): 585-595.PubMed CentralPubMedGoogle Scholar
- Fu YX: Statistical tests of neutrality of mutations against population growth, hitchhiking and background selection. Genetics. 1997, 147 (2): 915-925.PubMed CentralPubMedGoogle Scholar
- Rogers AR, Harpending H: Population growth makes waves in the distribution of pairwise genetic differences. Mol Biol Evol. 1992, 9 (3): 552-569.PubMedGoogle Scholar
- Harpending HC: Signature of ancient population growth in a low-resolution mitochondrial DNA mismatch distribution. Hum Biol. 1994, 66 (4): 591-600.PubMedGoogle Scholar
- Cheng HH, Levin I, Vallejo RL, Khatib H, Dodgson JB, Crittenden LB, Hillel J: Development of a genetic map of the chicken with markers of high utility. Poult Sci. 1995, 74 (11): 1855-1874.View ArticlePubMedGoogle Scholar
- Crooijmans RP, van der Poel JJ, Groenen MA: Functional genes mapped on the chicken genome. Anim Genet. 1995, 26 (2): 73-78.View ArticlePubMedGoogle Scholar
- Crooijmans RP, Dijkhof RJ, van der Poel JJ, Groenen MA: New microsatellite markers in chicken optimized for automated fluorescent genotyping. Anim Genet. 1997, 28 (6): 427-437. 10.1111/j.1365-2052.1997.00205.x.View ArticlePubMedGoogle Scholar
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.