Email updates

Keep up to date with the latest news and content from BMC Evolutionary Biology and BioMed Central.

Open Access Research article

A homogenous nature of native Chinese duck matrilineal pool

Da-Qian He12, Qing Zhu1, Shi-Yi Chen1, Hui-Ying Wang2, Yi-Ping Liu1* and Yong-Gang Yao3*

Author Affiliations

1 College of Animal Science and Technology, Sichuan Agriculture University, Ya'an, Sichuan, PR China

2 Institute of Animal Sciences, Shanghai Academy of Agricultural Sciences, Shanghai, PR China

3 State Key Laboratory of Genetic Resource and Evolution, Kunming Institute of Zoology, Chinese Academy of Sciences, Kunming, Yunnan, PR China

For all author emails, please log on.

BMC Evolutionary Biology 2008, 8:298  doi:10.1186/1471-2148-8-298

The electronic version of this article is the complete one and can be found online at: http://www.biomedcentral.com/1471-2148/8/298


Received:13 May 2008
Accepted:29 October 2008
Published:29 October 2008

© 2008 He et al; licensee 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.

Abstract

Background

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).

Results

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.

Conclusion

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.

Background

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 [1]. 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 [1]. 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 webcite 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. [2]. 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) [3]. 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 [4]. Analysis of the nuclear DNA sequence variation in representative Chinese duck breeds, wild mallard, and spot-billed duck by RAPD (random-amplified polymorphic DNA) [5], AFLP (amplified fragment length polymorphism) [6] and SSR (simple sequence repeat, microsatellite) markers [7] 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 [8]. After the initial domestication in China, the Chinese domestic ducks were subsequently introduced into Japan via Taiwan [9].

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.

Results

mtDNA control region and Cyt b sequence 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 EF126702EF126706, EF126708, EF126711EF126718, EF126721, EF126723EF126728, EF126731), five domestic duck samples (Anas platyrhynchos) were from Thailand (GenBank accession nos. EU013948EU013952), 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.

Additional file 1. Sequence variation of 41 mtDNA control region haplotypes identified in 449 domestic ducks and wild mallards. Variable sites were scored relative to the reference sequence (abbreviated as RS, GenBank accession number NC_009684). Dots (·) denote identity with the reference sequence, and the number of samples sharing the same haplotype is listed in the right column (under the capital N).

Format: DOC Size: 58KB Download file

This file can be viewed with: Microsoft Word ViewerOpen Data

Additional file 2. Haplotype distribution among different domestic duck breeds/populations and wild mallard populations.

Format: XLS Size: 33KB Download file

This file can be viewed with: Microsoft Excel ViewerOpen Data

Table 1. Sample information and genetic diversity of Chinese domestic ducks and mallards

thumbnailFigure 1. Sample locations in this study. Samples marked in black and white circles refer to the domestic ducks and the wild mallards, respectively. The abbreviations of different breeds/populations were defined in Table 1.

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.

thumbnailFigure 2. Rooted neighbour-joining tree of Chinese domestic ducks and mallards based on the mtDNA control region haplotypes. The tree was rooted by the Muscovy duck (Cairina moschata). The values on the branch are bootstrap support based on 1000 replications and those values lower than 50% were omitted.

thumbnailFigure 3. Network profile of the mtDNA control region haplotypes in Chinese domestic ducks and mallards. The links are labelled by the nucleotide positions to designate transitions; the single transversion at site 434 are further specified by adding suffix A. Recurrent mutations are underlined. The order of the mutations on a branch is arbitrary. The black colour denotes the wild mallards, while the white represents the domestic ducks. Circle area is proportional to haplotype frequency.

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 [31] and the mismatch distribution analysis [32] to estimate the demographic history of duck. A statistically significant rejection of the neutrality hypothesis revealed by the Fs test [31] and a unimodal or Poisson-like mismatch distribution [32] 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 [31] 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 [32] and was consistent with an explanation for population fragmentation if it was not caused by the small sample size [33] 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.

Table 2. Fu's Fs test for domestic duck and mallard samples

thumbnailFigure 4. Mismatch distributions for Chinese domestic ducks and mallards.

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.

Table 3. Main matrilineal components in Chinese domestic animals

thumbnailFigure 5. Mismatch distributions for domestic yak, goat, cattle, and chicken in China. The data were from our previous studies and others [21-23,39]. The sample sizes for yak, goat, cattle, and chicken were 52, 795, 209, and 900, respectively.

Discussion

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 [14] and the mtDNA control region sequence variation [3].

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 [9]. 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. [9] 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 [36]. Among the two haplotypes in the Northeast Asian group [9], 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. [9], the emerging pattern from their duck Cyt b sequences is consistent with our results. In fact, Hitosugi et al. [9] 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 [38]. 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 [29]. Study of the paternal markers in wild spot-billed duck, mallard, and domestic ducks will help to clarify this important issue.

Conclusion

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.

Data analysis

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 EU013948EU013952) 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 AY506945AY506947) 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 [45]. 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. [46] using program Network 4.1 http://www.fluxus-engineering.com/sharenet.htm webcite. We estimated the haplotype diversity (h) and nucleotide diversity (π) of the breeds/populations using DnaSP 4.10 [47]. The Fs values of the neutrality test [31] and the mismatch distributions [32] for the domestic ducks and wild mallards, as well as the total samples were analyzed using DnaSP 4.10 [47]. 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 EU677846EU678193 and EU678194EU678245, respectively.

Authors' contributions

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.

Acknowledgements

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.

References

  1. Chen GH, Wang KH, Wang JY, Ding C, Yang N: The genetic resources of poultry in China. Shanghai: Shanghai Scientific and Technical Publishers; 2004.

    (in Chinese)

  2. 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)

    OpenURL

  3. 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. PubMed Abstract | Publisher Full Text OpenURL

  4. 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)

    OpenURL

  5. 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)

    OpenURL

  6. 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)

    OpenURL

  7. 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)

    OpenURL

  8. Xie C-X: The breeding history of poultry in China. Beijing: Chinese Agriculture Publishers; 1995.

    (in Chinese)

  9. 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. Publisher Full Text OpenURL

  10. 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)

    OpenURL

  11. 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)

    OpenURL

  12. 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)

    OpenURL

  13. 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)

    OpenURL

  14. 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. Publisher Full Text OpenURL

  15. 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)

    OpenURL

  16. Maak S, Neumann K, von Lengerken G, Gattermann R: First seven microsatellites developed for the Peking duck (Anas platyrhynchos).

    Anim Genet 2000, 31:233. PubMed Abstract | Publisher Full Text OpenURL

  17. 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. Publisher Full Text OpenURL

  18. 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. PubMed Abstract | BioMed Central Full Text | PubMed Central Full Text OpenURL

  19. 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. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  20. 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. PubMed Abstract | Publisher Full Text OpenURL

  21. 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. PubMed Abstract | Publisher Full Text OpenURL

  22. 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. PubMed Abstract | Publisher Full Text OpenURL

  23. 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. PubMed Abstract | Publisher Full Text OpenURL

  24. Bruford MW, Bradley DG, Luikart G: DNA markers reveal the complexity of livestock domestication.

    Nat Rev Genet 2003, 4:900-910. PubMed Abstract | Publisher Full Text OpenURL

  25. 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. PubMed Abstract | Publisher Full Text OpenURL

  26. 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. PubMed Abstract | Publisher Full Text OpenURL

  27. 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. PubMed Abstract | BioMed Central Full Text | PubMed Central Full Text OpenURL

  28. 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. PubMed Abstract | Publisher Full Text OpenURL

  29. 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. Publisher Full Text OpenURL

  30. 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. Publisher Full Text OpenURL

  31. Fu Y-X: Statistical tests of neutrality of mutations against population growth, hitchhiking and background selection.

    Genetics 1997, 147:915-925. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  32. Rogers AR, Harpending H: Population growth makes waves in the distribution of pairwise genetic differences.

    Mol Biol Evol 1992, 9:552-569. PubMed Abstract | Publisher Full Text OpenURL

  33. 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. PubMed Abstract | Publisher Full Text OpenURL

  34. 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. PubMed Abstract | Publisher Full Text OpenURL

  35. Sorenson MD, Quinn TW: Numts: A challenge for avian systematics and population biology.

    Auk 1998, 115:214-221. OpenURL

  36. Ruokonen M, KVIST L, Lumme J: Close relatedness between mitochondrial DNA from seven Anser goose species.

    J Evol Biol 2000, 13:532-540. Publisher Full Text OpenURL

  37. 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. PubMed Abstract | BioMed Central Full Text | PubMed Central Full Text OpenURL

  38. Yao Y-G, Kong Q-P, Salas A, Bandelt H-J: Pseudo-mitochondrial genome haunts disease studies.

    J Med Genet 2008, in press.

    doi:10.1136/jmg.2008.059782

    PubMed Abstract | Publisher Full Text OpenURL

  39. 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, in press. OpenURL

  40. Ho SY, Larson G: Molecular clocks: when times are a-changin'.

    Trends Genet 2006, 22:79-83. PubMed Abstract | Publisher Full Text OpenURL

  41. 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. PubMed Abstract | Publisher Full Text OpenURL

  42. Graur D, Martin W: Reading the entrails of chickens: molecular timescales of evolution and the illusion of precision.

    Trends Genet 2004, 20:80-86. PubMed Abstract | Publisher Full Text OpenURL

  43. 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. PubMed Abstract | Publisher Full Text OpenURL

  44. 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. PubMed Abstract | Publisher Full Text OpenURL

  45. 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. PubMed Abstract | Publisher Full Text OpenURL

  46. Bandelt H-J, Forster P, Röhl A: Median-joining networks for inferring intraspecific phylogenies.

    Mol Biol Evol 1999, 16:37-48. PubMed Abstract | Publisher Full Text OpenURL

  47. Rozas J, Sanchez-DelBarrio JC, Messeguer X, Rozas R: DnaSP, DNA polymorphism analyses by the coalescent and other methods.

    Bioinformatics 2003, 19:2496-2497. PubMed Abstract | Publisher Full Text OpenURL

  48. 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. PubMed Abstract | Publisher Full Text OpenURL