- Research article
- Open Access
A homogenous nature of native Chinese duck matrilineal pool
BMC Evolutionary Biologyvolume 8, Article number: 298 (2008)
China, with around 30 unique breeds, has a diverse duck genetic pool. Currently, there is no systematic report which investigates the genetic diversity, phylogenetic relationship, and matrilineal genetic structure of these domestic breeds and wild mallards (Anas platyrhynchos).
In this study, we sequenced the mitochondrial DNA (mtDNA) control region segments in 278 domestic ducks (Anas platyrhynchos domestica) from 19 indigenous breeds/populations and 70 wild mallard samples and analyzed them together with the 101 control region sequences from published sources. Fifty-two samples were then sequenced for a cytochrome b (Cyt b) gene fragment to solidify the pattern emerged from the control region sequences. All domestic duck and wild mallard haplotypes were essentially indistinguishable and were clustered together in the phylogenetic tree. There was no geographic differentiation and breed/population-specific distribution of duck lineages.
Our results showed that unlike other domesticated farm animals in China such as chicken, cattle, goat, and yak with multiple matrilineal components, the matrilineal pool of Chinese ducks was homogenous.
Domestic ducks (Anas platyrhynchos domestica) play a key role in the agricultural and economic sectors of Asia. There are numerous domestic duck breeds in China including 27 indigenous breeds, two introduced breeds, and a few developing breeds according to the 2004 systematic field investigation . About 70% of the recorded breeds are distributed along the mid to downstream regions of the Yangtze and Pearl rivers, as well as, in the coastal districts . With the increasing demand for duck products, including meat, eggs, and down feathers, the duck breeding industry has flourished in China. In comparison to other domestic animals in China, the present conservation situation for local duck genetic resources is above average. Despite this, the overall population size of eight renowned breeds (Beijing duck, Youxian Sheldrake, Liancheng white duck, Jianchang duck, Jinding duck, Shaoxing duck, Putian black duck, and Gaoyou duck) has seen a sharp decline during the last decade. Fortunately, these local breeds are included in the National Genetic Resources Protection list http://www.agri.gov.cn/blgg/t20060609_626418.htm which ensures essential measurements are in place to sustainably manage these genetic resources.
The domestic duck is thought to have been domesticated independently from the wild mallard (Anas platyrhynchos), see a review by Li et al. . A recent study of the mtDNA D-loop sequence variation among nine duck breeds along the Yangtze-Huai River also suggested a single matrilineal origin from the wild mallard (Anas platyrhynchos) . However, comparisons of mtDNA polymorphisms among two local domestic ducks, the wild mallard and one wild spot-billed duck (Anas zonorhyncha), using RFLP (restriction fragment length polymorphism) showed that both wild duck species contributed to the gene pool of domestic ducks in Fujian Province, Southeast China . Analysis of the nuclear DNA sequence variation in representative Chinese duck breeds, wild mallard, and spot-billed duck by RAPD (random-amplified polymorphic DNA) , AFLP (amplified fragment length polymorphism)  and SSR (simple sequence repeat, microsatellite) markers  further supported this notion. Based on the archaeological assemblages, Xie proposed that Chinese native ducks were domesticated in the middle and lower region of the Yangtze River more than 2,000 years ago . After the initial domestication in China, the Chinese domestic ducks were subsequently introduced into Japan via Taiwan .
Overall, molecular studies investigating the duck breeds and populations in China are limited, focusing primarily on blood protein polymorphisms [10–12], random amplified polymorphic DNA [5, 13], and microsatellite markers [7, 14–17]. Much of the research is overshadowed by the superficiality of the research work (especially those published in Chinese) and the limited breeds and populations investigated. mtDNA sequence variation has been used extensively to study the genetic structure and matrilineal origin of farm animals [18–28]. However, mtDNA reports about domestic duck and wild mallard are relatively sparse [3, 9, 29, 30]. In this study, we analysed the genetic diversity and phylogeographic profiles of 374 specimens from 26 Chinese domestic duck breeds/populations and five Thailand native duck samples based on mtDNA sequence variation. To better understand the matrilineal origin of domestic ducks in China, we also sequenced 70 individuals from three wild mallard populations and compared them with the domestic samples. Our results revealed that the Chinese native duck gene pool is homogeneous, without clear geographic differentiation of the regional breed/population pools.
mtDNA control region and Cyt bsequence variation
A total of 449 domestic duck and wild mallard samples were analyzed in this study (Table 1 and Figure 1). Among these samples, 374 individuals were from 26 Chinese domestic duck breeds/populations (including 96 reported sequences, which were shared by 22 haplotypes with GenBank accession numbers EF126702–EF126706, EF126708, EF126711–EF126718, EF126721, EF126723–EF126728, EF126731), five domestic duck samples (Anas platyrhynchos) were from Thailand (GenBank accession nos. EU013948–EU013952), and 70 individuals were from three wild mallard populations in China. Forty-one haplotypes were identified and were defined by 30 polymorphic sites in the 481 bp fragment (Additional file 1). Among these haplotypes, H6 was predominant and occurred in 282 samples from all 30 breeds/populations; haplotypes H22, H1, and H5 were distributed in 32, 21, and 12 individuals, respectively. Nineteen haplotypes were shared by two to ten samples, while the remaining eighteen haplotypes were found in only a single individual. The five Thailand native ducks shared the most common haplotype H6. Haplotype H14, a one-step mutation derive of haplotype H6, was only found in three wild mallard samples, whereas the other wild mallard samples shared haplotypes with the domestic samples. The detailed haplotype distribution pattern among different breeds/populations is shown in Additional file 2.
Based on the haplotype distribution pattern of the mtDNA control region sequences, we further selected 52 samples and analyzed a Cyt b fragment (Table 1), to solidify the pattern emerged from the control region sequences. Three haplotypes (defined by two variable sites 15274 A>T and 15574 A>G) were discerned in 52 Cyt b sequences analyzed in this study. Specifically, variant A15274T was detected in four Jianchang ducks and one Jinding duck and caused an amino acid change from Leucine to Phenylalanine. Variant A15574G was synonymous and appeared in one Jinding duck.
Phylogenetic profile of domestic ducks
In the rooted neighbour-joining (NJ) tree based on the mtDNA control region haplotypes, all domestic duck and wild mallard samples were clustered together, whereas the spot-billed ducks formed another cluster (Figure 2). The wild mallard samples were intermingled with the domestic samples. There was no breed/population-specific clustering pattern. The network based on the mtDNA control region haplotypes revealed similar pattern as the NJ tree, with no essential distinction between the wild and domestic samples (Figure 3). Moreover, the network presented a star-like profile, consistent with a pattern of population expansion in the past. Twenty-three haplotypes differed from haplotype H6 by one mutation, and the remaining 17 haplotypes diverged from H6 by no more than three-mutation distance (Figure 3). The 70 wild mallard samples shared the predominant haplotype and eight other haplotypes with the domestic samples. The paucity of sequence variation in the Cyt b sequences was consistent with the less resolved tree based on the mtDNA control region sequences. We did not perform phylogenetic analysis for Cyt b sequences due to the lack of sequence variation.
Genetic diversity of duck breeds
Among the 30 duck breeds/populations analyzed in this study, the genetic diversity of each population/breed varied substantially (Table 1). Only one haplotype was detected in each sample of the crossbred populations of the Cherry Valley duck, the Gaoyou duck (CT), the Zhenyi duck and the Gaoyou duck (BY), as well as, the Thailand domestic duck (TH). The Jingjiang Sheldrake (JJ) had the highest haplotype diversity, while the Beijing duck (BD) had the highest nucleotide diversity. The haplotype diversity and nucleotide diversity in all the domestic ducks was 0.619 ± 0.029 and 0.00218 ± 0.00016, respectively, which was higher than those of the wild mallard (haplotype diversity, 0.461 ± 0.072; nucleotide diversity, 0.00190 ± 0.00036). The average haplotype diversity and nucleotide diversity in all 449 samples was 0.598 ± 0.027 and 0.00217 ± 0.00015, respectively. The small sample size of some breeds/populations and their unique breeding histories may account for the high variability of genetic diversity (Table 1). In addition, potentially biased sampling, e.g. from one cultured population that is hard to avoid the potentially close affinity of individuals, could reduce the haplotype and nucleotide diversities.
Population demographic history
We performed the Fs test of Fu  and the mismatch distribution analysis  to estimate the demographic history of duck. A statistically significant rejection of the neutrality hypothesis revealed by the Fs test  and a unimodal or Poisson-like mismatch distribution  would characterize a population expansion in the past. In contrast, a population with a constant size in the past would have a multimodal mismatch distribution and insignificant Fs test result [31, 32]. As shown in Table 2, the Fs tests  for the domestic duck samples, wild mallard samples, and the total samples were all statistically significant (P < 0.05). The mismatch distributions for all samples (including domestic ducks and wild mallards) showed a (biased) unimodal shape (Figure 4). Similar distribution shape was detected for the domestic ducks. However, the wild mallard samples showed a seemingly bumpy mismatch distribution, which was characteristic of a population in equilibrium  and was consistent with an explanation for population fragmentation if it was not caused by the small sample size  or other demographic effects. Combined with the star-like shape of the network profile (Figure 3), we speculated that the domestic duck in China might have undergone population expansion in the past.
To further compare the results with other farm animals in China, we estimated the mismatch distributions for animals, such as domestic yak, goat, cattle, and chicken (Figure 5). All these farm animals showed multimodal mismatch distribution patterns, especially yak, goat, and cattle, which were apparently caused by multiple divergent matrilineal components in each pool (Table 3), instead of indicating a population in equilibrium. Indeed, the multimodal mismatch distributions of Chinese cattle were derived from the existence of both Bos taurus and Bos indicus types [22, 34]. An analysis of the mismatch distributions based on Bos taurus or Bos indicus types only showed a unimodal distribution (data not shown). Overall, the mismatch distribution patterns of goat, yak, cattle, and chicken were dramatically different when compared to that of the domestic duck. Note that the comparison between the domestic duck and other farm animals was restricted to their current matrilineal genetic structure and the domestication history of each animal was not considered.
In this study we discerned the genetic diversity and origin of Chinese native ducks. The weak phylogeographic structure of Chinese domestic ducks and wild mallards is clearly discerned in the NJ tree and the network profile based on the mtDNA control region haplotypes. The results suggest that the matrilineal pool of Chinese domestic ducks is homogeneous and all samples can be grouped into one haplogroup, which is analogous to a single matrilineal origin. The Cyt b gene sequences in 52 domestic duck and wild mallard samples that were selected based on their control region sequence information contain only two variable sites, which further supports the homogeneous nature of the duck matrilineal pool. The overall diversity of Chinese domestic duck estimated in our study was similar to the recent reports based on the microsatellite markers  and the mtDNA control region sequence variation .
In a recent study, Hitosugi and colleagues analyzed the Cyt b sequence variation of duck samples from East and Southeast Asia, and they found that the samples from Taiwan and Japan were grouped into the Northeast Asian group, whereas those from Indonesia, Vietnam, and west Malaysia were clustered into the Southeast Asian group . This paraphyletic pattern of duck mtDNAs directly contrasts the clustering pattern of Chinese samples. Because we analyzed up to 26 domestic duck breeds/populations (which covers about 80% of the recorded native duck breeds/populations in nearly all current habitats across China; Table 1 and Figure 1) and three wild mallard populations, we believed that our sampling could represent the entire matrilineal pool of Chinese duck (at least for the Mainland region) and we should identify lineages belonging to the two groups. An audit for the Cyt b sequences that were grouped into the Southeast Asian group reported by Hitosugi et al.  showed that all these sequences should be nuclear mitochondrial pseudogenes (NUMTs), as these Cyt b fragments contained a premature termination codon. The invasion of NUMT is not infrequent in birds [35–37] and it had been identified in several goose species . Among the two haplotypes in the Northeast Asian group , one shared sequence with the duck reference sequence (GenBank accession number NC_009684) and the other was a close derivative of the main haplotype in Chinese duck. Therefore, after eliminating the NUMTs in the data of Hitosugi et al. , the emerging pattern from their duck Cyt b sequences is consistent with our results. In fact, Hitosugi et al.  also sequenced the mtDNA control region sequences for these East Asian duck samples and observed no variation; a result that concurs with our findings.
An important issue which needs clarification concerns our belief that the mtDNA sequence data generated in the present study were free of NUMTs and DNA contamination, although we did not perform specific experiments for confirmation. Three lines of evidence justify our claim. First, the mutation spectrum left by spurious amplification of a NUMT could be detected because it usually is considerably different from the authentic natural spectrum . Among all blood/tissue samples analyzed in this study, we did not detect any sequences showing a remarkable divergence (which would be resulted from the different evolutionary rate between authentic mtDNA and the NUMTs). Second, we obtained consistent patterns between two different mtDNA fragments and all Cyt b sequences could be fully translated. Third, DNA contamination between different samples and co-amplification of NUMTs and authentic mtDNA would result in many heteroplasmic mutations in the sequence, which was not observed in this study.
Interestingly, the matrilineal pool homogeneity of Chinese native ducks differed substantially with other domesticated Chinese animals, as judged by the numbers of main matrilineal components (Table 3) and the mismatch distribution patterns (Figures 4 and 5). In contrast to the unimodal shape for Chinese native ducks, the mismatch distributions for Chinese domestic yak [23, 27], goat [26, 39], cattle [22, 25, 34], and chicken [20, 21] were bumpy or multimodal, which reflects the highly divergent lineages within the population (Table 3). Multiple domestication events of these latter farm animals could account for the difference in their respective matrilineal pools. We did not attempt to discern the domestication time for Chinese domestic ducks in this study because of concerns with the molecular clock [40–42]. Similarly, the indistinguishable matrilineal components among different breeds/regional populations of both domestic ducks and wild mallards have deterred any attempts to elucidate the initial domestication site(s) of Chinese domestic ducks.
Kulikova and colleagues detected two mtDNA haplogroups (A and B) of mallards that were collected from western Russia to mainland Alaska, and found that the haplogroup A was predominantly distributed in western Russia and North Asia territories [29, 30]. All Chinese domestic duck and wild mallard mtDNAs in this study could be grouped into haplogroup A, suggesting a close relationship between the mallard matrilineal pools of China and North Asia. The fact that five Thailand domestic duck samples shared haplotype H6 with Chinese domestic ducks and wild mallards is also evidence for a genetic continuality among these regional pools. Nonetheless, due to the small sample size of the Thailand domestic duck, such a speculation should be interpreted with caution.
During the past decade, many molecular studies (mainly based on nuclear genetic markers) showed that the wild mallard and spot-billed duck have contributed to the genetic pool of Chinese domestic ducks, with the mallard contributing a greater proportion [1, 4–7]. However, in the phylogenetic tree based on the mtDNA control region sequences, the three spot-billed ducks were robustly separated and none of the Chinese domestic ducks and wild mallards (including the Thailand sample) were clustered with the spot-billed ducks. A possible explanation for this unique pattern is that the genetic contribution of spot-billed duck to Chinese domestic duck pool might be via a male-biased gene introgression, as have been described between the eastern spot-billed duck and the mallard . Study of the paternal markers in wild spot-billed duck, mallard, and domestic ducks will help to clarify this important issue.
In this study, we analyzed Chinese domestic duck and wild mallard mtDNA control regions and partial Cyt b sequences variation. The results showed that all domestic duck and wild mallard haplotypes were essentially indistinguishable and were clustered together in the phylogenetic tree. The matrilineal pool of Chinese domestic ducks and mallards was rather homogeneous, without clear geographic differentiation of the regional breed/population pools. This pattern substantially differs from other farm animals in China.
Materials and methods
Sampling, DNA amplification and sequencing
In total, we collected 348 blood/tissue samples from 19 Chinese domestic duck breeds/populations (N = 278) and three wild mallard populations (N = 70) in East, Central and South China (Table 1 and Figure 1).
Genomic DNA was extracted from blood/tissue samples using a standard phenol/chloroform method. The mtDNA control region fragments (481 bp) were amplified and sequenced using primer pair L194: 5'-CCTACCTATCGGACTACCCTC-3'/H716: 5'-GCAGGTGTGTCCAGGCTTAGA-3', which was designed in this study. PCR amplification was performed in a 50-μl reaction mixture containing 100 ng of DNA, 10 mM Tris – HCl (pH 8.3), 2.5 mM MgCl2, 50 mM KCl, 10 pM of each primer, and 1 unit of Taq polymerase following 33 cycles of 50 s at 94°C, 40 s at 58°C, and 90 s at 72°C. PCR products were purified on spin columns and were directly sequenced for both strands using Big Dye Terminator v3.1 Cycle Sequencing Kit (Applied Biosystems, USA) on an ABI PRISM® 3100 DNA sequencer according to the manufacturer's manual.
The Cyt b fragment (508 bp) was amplified and sequenced using primer pair designed in this study (L15230: 5'-ACCCTGACCCGATTCTTC-3'/H15774: 5'-GATGCGAGTTGCCCGATGA-3') under the same condition as for mtDNA control region analysis but with a modification of annealing temperature (55°C) during the PCR amplification.
The domestic duck and wild mallard mtDNA control region and Cyt b sequences were edited and aligned using the DNAstar program (DNAS Inc, Madison, WI, USA). The sequences obtained were aligned based on the complete mitochondrial genome of Anas platyrhynchos (GenBank accession number NC_009684). We retrieved 106 Chinese domestic duck mtDNA control region sequences from nine native breeds/populations (including two breeds that were overlapped with our sampling) and five Thailand native ducks (GenBank accession numbers EU013948–EU013952) from GenBank. All these reported control region sequences were aligned with the new sequences generated in this study and were truncated to 481 bp fragments in the following analyses. Ten reported mtDNA sequences (shared by nine haplotypes; GenBank accession numbers EF126707, EF126709, EF126710, EF126719, EF126720, EF126722, EF126729, EF126730, and EF126732) were excluded for further analysis because of apparent sequencing and reading errors by applying an error-pinpointing approach described in our previous study and others [43, 44]. Three mtDNA control region sequences of the spot-billed duck (Anas zonorhyncha) were retrieved from GenBank (accession numbers AY506945–AY506947) and were included in the NJ tree analysis to discern their phylogenetic relationship to domestic samples.
The NJ tree was constructed based on the haplotypes identified in domestic ducks, mallards, and spot-billed ducks and the Kimura-2-parameters model using MEGA 4.0 . The Muscovy duck (Cairina moschata; GenBank accession number L16769) was chosen as the outgroup to root the tree. Median-joining network of the mtDNA control region sequence haplotypes was constructed according to Bandelt et al.  using program Network 4.1 http://www.fluxus-engineering.com/sharenet.htm. We estimated the haplotype diversity (h) and nucleotide diversity (π) of the breeds/populations using DnaSP 4.10 . The Fs values of the neutrality test  and the mismatch distributions  for the domestic ducks and wild mallards, as well as the total samples were analyzed using DnaSP 4.10 . We also estimated the mismatch distributions for some reported Chinese domestic animals, such as chicken, cattle, yak, and goat [21–23, 39], in comparison to that of the duck. The new duck mtDNA control region and Cyt b sequences generated in this study were deposited in GenBank under accession numbers EU677846–EU678193 and EU678194–EU678245, respectively.
Chen GH, Wang KH, Wang JY, Ding C, Yang N: The genetic resources of poultry in China. 2004, Shanghai: Shanghai Scientific and Technical Publishers, (in Chinese)
Li H-F, Qian K, Xu W-J, Gu R, Zhao R-G: The diversity of worldwide domestic ducks. China Poult. 2006, 28: 58-61. (in Chinese)
Zhang T-J, Li H-F, Chen K-W, Chang H, Tang Q-P, Zhang J-X: Genetic diversity and systematic evolution of Chinese domestic ducks along the Yangtze-Huai River. Biochem Genet. 2007, 45: 823-837. 10.1007/s10528-007-9121-y.
Chen Y, Zhao Y, Lu L, Lai Y: Comparative study of mitochondrial DNA restriction maps of genus of Anas ducks in Fujian Province. J Xiamen Univ Nat Sci. 1999, 38: 108-111. (in Chinese)
Chen Y-X, Zuo Z-H: Study on genetic relationship and diversity in ducks by RAPD. J Xiamen Univ Nat Sci. 2001, 40: 141-145. (in Chinese)
Yan F-H, Zuo Z-H, Chen M, Song Y-Q, Lv L-J, Chen Y-X: Analysis of genetic diversity and relationship among some Chinese domestic ducks and wild ducks using AFLP. J Xiamen Univ Nat Sci. 2005, 44: 729-733. (in Chinese)
Chen Q, Huang Z-Y, Yan F-H, Lu L-J: Analysis of genetic relationship among some local domestic ducks in Fujian Province. J Quanzhou Norm Univ Nat Sci Ed. 2006, 24: 90-94. (in Chinese)
Xie C-X: The breeding history of poultry in China. 1995, Beijing: Chinese Agriculture Publishers, (in Chinese)
Hitosugi S, Tsuda K, Okabayashi H, Tanabe Y: Phylogenetic relationships of mitochondrial DNA cytochrome b gene in East Asian ducks. J Poult Sci. 2007, 44: 141-145. 10.2141/jpsa.44.141.
Lin Y-Q, Wang L-Q, Zeng F-T, Qiu X-P: Electrophoresis patterns of blood plasma protein polymorphisms and the association with quantitative trait loci of duck. J Sichuan Agric Univ. 1996, 14: 85-91. (in Chinese)
Yang G-Q, Han L-M, Liu M, Tian H: Studies on the serum prealbumin polymorphisms of various strains of the Peking duck. J Shenyang Agric Univ. 2001, 32: 436-438. (in Chinese)
Zhang T, Lu L, Shen J, Chen J, Zhao R: Comparative study on serum esterase polymorphism between two lines of Shaoxing ducks. China Poult. 2004, 26: 10-12. (in Chinese)
Zuo Z-H, Chen Y-X, Lv L-J: Studies on genetic diversity and molecular marker of the Jinding duck. J Xiamen Univ Nat Sci. 2004, 43: 256-259. (in Chinese)
Li H, Yang N, Chen K, Chen G, Tang Q, Tu Y, Yu Y, Ma Y: Study on molecular genetic diversity of native duck breeds in China. World's Poult Sci J. 2006, 62: 603-611. 10.1017/S0043933906001152.
Ji W, Qian K, Li H, Ge H, Zhao XT, Duan X, Chen K: Analysis of genetic diversity of 7 domestic duck populations conserved in national waterfowl gene pool. China Poult. 2006, 28: 72-74. (in Chinese)
Maak S, Neumann K, von Lengerken G, Gattermann R: First seven microsatellites developed for the Peking duck (Anas platyrhynchos). Anim Genet. 2000, 31: 233-
Maak S, Wimmers K, Weigend S, Neumann K: Isolation and characterization of 18 microsatellites in the Peking duck (Anas platyrhynchos) and their application in other waterfowl species. Mol Ecol Notes. 2003, 3: 224-227. 10.1046/j.1471-8286.2003.00405.x.
Wu G-S, Yao Y-G, Qu K-X, Ding Z-L, Li H, Palanichamy Mg, Duan Z-Y, Li N, Chen Y-S, Zhang Y-P: Population phylogenomic analysis of mitochondrial DNA in wild boars and domestic pigs revealed multiple domestication events in East Asia. Genome Biol. 2007, 8: R245-10.1186/gb-2007-8-11-r245.
Naderi S, Rezaei HR, Taberlet P, Zundel S, Rafat SA, Naghash HR, El-Barody MA, Ertugrul O, Pompanon F: Large-scale mitochondrial DNA analysis of the domestic goat reveals six haplogroups with high diversity. PLoS ONE. 2007, 2: e1012-10.1371/journal.pone.0001012.
Liu Y-P, Zhu Q, Yao Y-G: Genetic relationship of Chinese and Japanese gamecocks revealed by mtDNA sequence variation. Biochem Genet. 2006, 44: 19-29.
Liu Y-P, Wu G-S, Yao Y-G, Miao Y-W, Luikart G, Baig M, Beja-Pereira A, Ding Z-L, Palanichamy MG, Zhang Y-P: Multiple maternal origins of chickens: out of the Asian jungles. Mol Phylogenet Evol. 2006, 38: 12-19. 10.1016/j.ympev.2005.09.014.
Lai S-J, Liu Y-P, Liu Y-X, Li X-W, Yao Y-G: Genetic diversity and origin of Chinese cattle revealed by mtDNA D-loop sequence variation. Mol Phylogenet Evol. 2006, 38: 146-154. 10.1016/j.ympev.2005.06.013.
Lai S-J, Chen S-Y, Liu Y-P, Yao Y-G: Mitochondrial DNA sequence diversity and origin of Chinese domestic yak. Anim Genet. 2007, 38: 77-80. 10.1111/j.1365-2052.2007.01555.x.
Bruford MW, Bradley DG, Luikart G: DNA markers reveal the complexity of livestock domestication. Nat Rev Genet. 2003, 4: 900-910. 10.1038/nrg1203.
Chen S-Y, Liu Y-P, Wang W, Gao C-Z, Yao Y-G, Lai S-J: Dissecting the matrilineal components of Tongjiang cattle from southwest China. Biochem Genet. 2008, 46: 206-215. 10.1007/s10528-008-9144-z.
Chen S-Y, Su Y-H, Wu S-F, Sha T, Zhang Y-P: Mitochondrial diversity and phylogeographic structure of Chinese domestic goats. Mol Phylogenet Evol. 2005, 37: 804-814. 10.1016/j.ympev.2005.06.014.
Guo S, Savolainen P, Su J, Zhang Q, Qi D, Zhou J, Zhong Y, Zhao X, Liu J: Origin of mitochondrial DNA diversity of domestic yaks. BMC Evol Biol. 2006, 6: 73-10.1186/1471-2148-6-73.
Lei CZ, Zhang W, Chen H, Lu F, Liu RY, Yang XY, Zhang HC, Liu ZG, Yao LB, Lu ZF, Zhao ZL: Independent maternal origin of Chinese swamp buffalo (Bubalus bubalis). Anim Genet. 2007, 38: 97-102. 10.1111/j.1365-2052.2007.01567.x.
Kulikova IV, Zhuravlev YN, McCracken KG: Asymmetric hybridization and sex-biased gene flow between eastern spot-billed ducks (Anas zonorhyncha) and maliards (A. platyrhynchos) in the Russia far east. Auk. 2004, 121: 930-949. 10.1642/0004-8038(2004)121[0930:AHASGF]2.0.CO;2.
Kulikova IV, Drovetski SV, Gibson DD, Harrigan RJ, Rohwer S, Sorenson MD, Winker K, Zhuravlev YN, McCracken KG: Phylogeography of the mallard (Anas platyrhynchos): hybridization, dispersal, and lineage sorting contribute to complex geographic structure. Auk. 2005, 122: 949-965. 10.1642/0004-8038(2005)122[0949:POTMAP]2.0.CO;2.
Fu Y-X: Statistical tests of neutrality of mutations against population growth, hitchhiking and background selection. Genetics. 1997, 147: 915-925.
Rogers AR, Harpending H: Population growth makes waves in the distribution of pairwise genetic differences. Mol Biol Evol. 1992, 9: 552-569.
Yao Y-G, Nie L, Harpending H, Fu Y-X, Yuan Z-G, Zhang Y-P: Genetic relationship of Chinese ethnic populations revealed by mtDNA sequence diversity. Am J Phys Anthropol. 2002, 118: 63-76. 10.1002/ajpa.10052.
Lei CZ, Chen H, Zhang HC, Cai X, Liu RY, Luo LY, Wang CF, Zhang W, Ge QL, Zhang RF, Lan XY, Sun WB: Origin and phylogeographical structure of Chinese cattle. Anim Genet. 2006, 37: 579-582. 10.1111/j.1365-2052.2006.01524.x.
Sorenson MD, Quinn TW: Numts: A challenge for avian systematics and population biology. Auk. 1998, 115: 214-221.
Ruokonen M, KVIST L, Lumme J: Close relatedness between mitochondrial DNA from seven Anser goose species. J Evol Biol. 2000, 13: 532-540. 10.1046/j.1420-9101.2000.00184.x.
Pereira SL, Baker AJ: Low number of mitochondrial pseudogenes in the chicken (Gallus gallus) nuclear genome: implications for molecular inference of population history and phylogenetics. BMC Evol Biol. 2004, 4: 17-24. 10.1186/1471-2148-4-17.
Yao Y-G, Kong Q-P, Salas A, Bandelt H-J: Pseudo-mitochondrial genome haunts disease studies. J Med Genet. 2008, doi:10.1136/jmg.2008.059782,
Liu Y-P, Cao S-X, Chen S-Y, Yao Y-G, Liu T-Z: Genetic diversity of Chinese domestic goat based on the mitochondrial DNA sequence variation. J Anim Breed Genet. 2008,
Ho SY, Larson G: Molecular clocks: when times are a-changin'. Trends Genet. 2006, 22: 79-83. 10.1016/j.tig.2005.11.006.
Ho SY, Phillips MJ, Cooper A, Drummond AJ: Time dependency of molecular rate estimates and systematic overestimation of recent divergence times. Mol Biol Evol. 2005, 22: 1561-1568. 10.1093/molbev/msi145.
Graur D, Martin W: Reading the entrails of chickens: molecular timescales of evolution and the illusion of precision. Trends Genet. 2004, 20: 80-86. 10.1016/j.tig.2003.12.003.
Bandelt H-J, Lahermo P, Richards M, Macaulay V: Detecting errors in mtDNA data by phylogenetic analysis. Int J Legal Med. 2001, 115: 64-69. 10.1007/s004140100228.
Yao Y-G, Bravi CM, Bandelt H-J: A call for mtDNA data quality control in forensic science. Forensic Sci Int. 2004, 141: 1-6. 10.1016/j.forsciint.2003.12.004.
Tamura K, Dudley J, Nei M, Kumar S: MEGA4: Molecular Evolutionary Genetics Analysis (MEGA) software version 4.0. Mol Biol Evol. 2007, 24: 1596-1599. 10.1093/molbev/msm092.
Bandelt H-J, Forster P, Röhl A: Median-joining networks for inferring intraspecific phylogenies. Mol Biol Evol. 1999, 16: 37-48.
Rozas J, Sanchez-DelBarrio JC, Messeguer X, Rozas R: DnaSP, DNA polymorphism analyses by the coalescent and other methods. Bioinformatics. 2003, 19: 2496-2497. 10.1093/bioinformatics/btg359.
Chen S-Y, Zhou F, Xiao H, Sha T, Wu S-F, Zhang Y-P: Mitochondrial DNA diversity and population structure of four Chinese donkey breeds. Anim Genet. 2006, 37: 427-429. 10.1111/j.1365-2052.2006.01486.x.
We thank the reviewers for helpful comments during the revision of the manuscript and Ms. Aleah Smith for help in editing the language. This work was financially supported by Chinese national key technologies R & D programme (Grant No: 2006BAD14B06), the project of educational bureau of Sichuan Province (Grant No: 07ZA070), and the National 863 Project of China (No. 2008AA101001). Y-GY was supported by the "Century Program" (or Hundreds-Talent Program) of the Chinese Academy of Sciences.
DQH collected samples, carried out the molecular genetic studies, participated in the sequence alignment, analyzed the data, and wrote the manuscript. QZ participated in the design of the study and helped to draft the manuscript. SYC collected the samples, carried out the molecular genetic studies, participated in the sequence alignment, analyzed the data, and wrote the manuscript. HYW collected samples. YPL conceived of the study, and participated in its design and coordination and wrote the manuscript. YGY conceived of the study, supervised the data analyses and interpretation, and wrote the manuscript. All authors read and approved the final manuscript.