Widespread horizontal genomic exchange does not erode species barriers among sympatric ducks
© Kraus et al; licensee BioMed Central Ltd. 2012
Received: 8 November 2011
Accepted: 2 April 2012
Published: 2 April 2012
The study of speciation and maintenance of species barriers is at the core of evolutionary biology. During speciation the genome of one population becomes separated from other populations of the same species, which may lead to genomic incompatibility with time. This separation is complete when no fertile offspring is produced from inter-population matings, which is the basis of the biological species concept. Birds, in particular ducks, are recognised as a challenging and illustrative group of higher vertebrates for speciation studies. There are many sympatric and ecologically similar duck species, among which fertile hybrids occur relatively frequently in nature, yet these species remain distinct.
We show that the degree of shared single nucleotide polymorphisms (SNPs) between five species of dabbling ducks (genus Anas) is an order of magnitude higher than that previously reported between any pair of eukaryotic species with comparable evolutionary distances. We demonstrate that hybridisation has led to sustained exchange of genetic material between duck species on an evolutionary time scale without disintegrating species boundaries. Even though behavioural, genetic and ecological factors uphold species boundaries in ducks, we detect opposing forces allowing for viable interspecific hybrids, with long-term evolutionary implications. Based on the superspecies concept we here introduce the novel term "supra-population" to explain the persistence of SNPs identical by descent within the studied ducks despite their history as distinct species dating back millions of years.
By reviewing evidence from speciation theory, palaeogeography and palaeontology we propose a fundamentally new model of speciation to accommodate our genetic findings in dabbling ducks. This model, we argue, may also shed light on longstanding unresolved general speciation and hybridisation patterns in higher organisms, e.g. in other bird groups with unusually high hybridisation rates. Observed parallels to horizontal gene transfer in bacteria facilitate the understanding of why ducks have been such an evolutionarily successful group of animals. There is large evolutionary potential in the ability to exchange genes among species and the resulting dramatic increase of effective population size to counter selective constraints.
Biology has seen the proposition of several species concepts. Of these, the biological species concept  is historically the most influential; according to it all individuals belong to the same species if they produce viable and fertile offspring in nature, i.e., they share a common gene pool. To account for inherent difficulties to test this concept in practice, especially in allopatric populations that never encounter each other, biologists tend to supplement it by elements of the morphospecies concept (which is as old as the study of nature). With the advance of molecular genetic data over the past decades many researchers now define species by genetic characteristics rather than morphological ones because genetics provides a means of actually measuring recent or ongoing genetic connectivity between species . Species boundaries are strengthened by accumulation of genomic incompatibilities preventing formation of zygotes, so called Dobzhansky-Muller incompatibilities [3–5]. Once evolved, post-zygotic isolation is irreversible, in contrast to pre-zygotic barriers such as mate recognition. There is much evidence that post-zygotic barriers evolve slowly in birds [5, 6], potentially contributing to the high rates of hybridisation observed in this group  and explaining why genetic distances can be low in spite of large morphological differences .
When populations diverge into species their gene pools become disconnected, and even in the absence of ecological differentiation stochastic effects, i.e., genetic drift, will drive each new species towards increased differentiation. If introgression of genetic material of one species into another occurs regularly enough in the absence of genomic incompatibility, one would expect that these events oppose genetic drift by exchange of alleles that the two subsequently will have in common. Such potential sharing of alleles at genetic loci through genetic admixture can directly be observed by the study of genetic markers. One type of genetic marker that has recently received a lot of attention is the 'single nucleotide polymorphism' (SNP) . Due to the abundance of SNPs in genomes and suitability for high automation in genotyping, SNPs can be characterised in large numbers, yielding a representative image of an entire genome. With SNP data from multiple species, one can study the sharing of genetic material at the same loci, providing a new means of studying species divergence by the speed of loss of genetic coherence.
Results and Discussion
Genotypic differentiation between Anas platyrhynchosand other duck species
We screened 364 SNPs developed for the mallard, Anas platyrhynchos,  in the genomes of six duck species, five of genus Anas and one of Aythya, the latter mainly for outgroup comparison: Anas platyrhynchos (N = 197), Anas acuta (northern pintail, N = 7), Anas crecca, (common teal, N = 9), Anas penelope (Eurasian wigeon, N = 14), Anas strepera (gadwall, N = 10) and Aythya fuligula (tufted duck, N = 17). The SNPs were evaluated for minor allele frequency (MAF) spectrum, Hardy-Weinberg equilibrium and linkage disequilibrium in Anas platyrhynchos from nine localities on three continents. The great majority of SNPs does not significantly deviate from neutrality and are unlinked.
SNP sharing among duck species is unexpectedly high
Interaction between population size and persistence time
Census size Nc
Effective size Ne
Persistence time (mean)
p = 0.5
p = 0.1
Increased population size by ongoing interspecific hybridisation
What can then explain the high level of shared polymorphisms? We argue that these (and other closely related) duck species are part of a superspecies complex, here defined as a group of distinct species that frequently hybridise, with fertile offspring as the result. The superspecies concept was put forward by Mayr in 1931 , as a translation of the German expression Artenkreis, based on the work of Rensch . Initially, it was used to assign species status to allopatric "races" that were too distinct to be lumped into the same species [27–29] (superspecies sensu stricto). Later, the definition was widened by Kiriakoff  and Mayr and Short  to be no longer exclusive to allopatric populations. For the Anas platyrhynchos complex this concept has previously been used by Scherer . Being aware that "superspecies" is not an official taxonomic category we here choose to use the term superspecies (sensu lato) to embrace the sympatric distribution of interbreeding duck species. In doing so, we do not attempt to redefine nomenclatural classification schemes, nor do we propose to change current nomenclature. The term superspecies is clearly "an evolutionary taxonomy category but not nomenclatural rank" , thus to be preferred when studying biological systems rather than nomenclature.
There is longstanding anecdotal, morphological and experimental evidence for high hybridisation rates in ducks [7, 12, 22], but molecular proof has been limited thus far. Two studies using mitochondrial DNA in the Anas rubripes/platyrhynchos  and Anas zonorhyncha/platyrhynchos  complexes confirm hybridisation between these species. These findings were corroborated by studies investigating one to two nuclear markers [35, 36]. Our study, using shared polymorphisms at hundreds of independent loci across the entire genome provides a more powerful means of analysing gene pool connectivity between closely related species and our results are consistent with a high level of genetic transfer between species via hybrid production and backcrossing.
A STRUCTURE  analysis identified several cases where genetic admixture from other species seems supported by their genotypes. When all six duck species were analysed jointly with the genetic clustering software STRUCTURE, all non-Anas platyrhynchos individuals were assigned to the same cluster (Additional file 1). Anas acuta individuals in particular showed partial Anas platyrhynchos genome admixture, and many Anas platyrhynchos individuals displayed some admixture from other species. When Anas platyrhynchos individuals were excluded, STRUCTURE assigned Anas penelope, Anas strepera and Aythya fuligula individuals to their species specific clusters, although one Anas strepera individual (ANST001) was almost fully assigned to Anas penelope. Anas acuta and Anas crecca were lumped into one cluster, and the hybrid was correctly assigned to that cluster by only 50% of its genome (Additional file 2). Excluding the hybrid from analysis did not alter the assignment of these two species to the same cluster. The same data sets were analysed with comparable settings in the software InStruct , which does not assume Hardy-Weinberg equilibrium in the inferred populations, and yields qualitatively similar results as the STRUCTURE analysis. This may be direct evidence of partial gene pool sharing between species, hence the establishment of a superspecies complex.
For example, a superspecies complex comprising Anas platyrhynchos, Anas acuta and Anas crecca would have a joint census population size of 31 million individuals and hence an Ne of 3.1 million (see methods for sources and assumptions), although sub-division of this possible superspecies due to assortative mating makes this an over-estimate. However, an Ne of 3.1 million results in a mean persistence time of almost 9 million years (for initial MAF = 0.5). With an estimated most recent common ancestor at 6.4 Mya, these species could have on average retained even SNPs of lower MAF = 0.2. We refer to this analysis as 'persistence time analysis'.
Species status and the supra-population concept
The ducks studied here have not only remained morphologically distinct, their genetic cluster species designation  is strongly supported by principal component analysis of SNP genotypes: we find clear genetic differentiation between Anas platyrhynchos and the other duck species, as well as among these (Figure 2c). Even though all these species live in sympatry, such a combined population is highly structured by assortative mating. While geographical substructure would be indicated by the term "meta-population", the situation in ducks leads us to define a new term that does not have a geographical connotation: "supra-population". We define a supra-population as a group of individuals that are part of the same sympatric superspecies complex and within which natural hybridisation occurs. Individuals of a superspecies complex are genetically-connected hybridising species, in which species barriers are primarily maintained by pre-zygotic factors.
A new model of speciation in ducks
Genomic incompatibilities usually lead to irreversible post-zygotic isolation of populations, but other, reversible, barriers can also be strong drivers of speciation. Visual cues have been identified as drivers of speciation in sexually dimorphic bird species [8, 40] while sexual imprinting alone can explain assortative mating in modeling studies . An empirical example from another Anatid species, the snow goose Anser caerulescens, which has two wide ranging colour morphs, nicely illustrates the case . At any rate, a model for speciation in ducks must be able to explain the observed pattern of genetic and morphological differentiation in spite of the high degree horizontal gene exchange.
Paleogeographic and paleoclimatic evidence suggest that ecological conditions have been favourable for a duck radiation 6-12 Mya. This late Miocene period was warm and humid [43, 44], but in transition towards a colder climate. Precipitation remained relatively high [45–47], making wetlands abundant and turning large inland salt water bodies brackish or even freshwater (e.g., Lake Pannon in Eurasia [48–50]). Globally, during this transition towards a colder, wet climate tropical forests were largely replaced by open grasslands [51–53], a habitat well suited for ducks. The fossil record of ducks beyond the Pleistocene is still very poor  but the few studies on the subject suggest that morphological change in respective duck species has been very limited over the last few million years [55, 56], after a larger waterfowl species turn-over 15-23 million years ago . The first fossil that resembles Anas platyrhynchos is thought to be from the late Pliocene, about 5 Mya . This is close to the suggested lower bound of divergence times of some Anas species in the latest phylogeny of Anatidae . We propose that an Anas-like duck split into multiple sister morphs sympatrically and simultaneously at that time, subsequently diverging by assortative mating. Our results indicate that the resulting cluster of species still exchanges portions of their genomes. We argue that since branching off of the Anas clade at least 6 Mya these mostly sympatric species remain separate by isolating mechanisms other than genetic incompatibilities, mostly by assortative mating. Though we acknowledge that this speciation scenario rests on the assumption of widespread sympatry for millions of years, we feel comfortable in making this claim. Although we only sampled five species for the present study, our model system sensu lato is the specious genus Anas, and even though species distributions change over time there certainly have always been several Anas species living in sympatry.
Theoretical studies suggest that sexual imprinting can drive speciation even in sympatry . Moreover, experimental manipulations clearly demonstrate that individuals of Anas platyrhynchos can be imprinted on nearly any species of waterfowl but when raised in isolation they recognise conspecifics as mates . This suggests that imprinting is important but incomplete in ducks; genetic factors also contribute to mate recognition. The presence of assortative mating and recognition mechanisms are prerequisites for sympatric speciation leading to a superspecies complex around Anas platyrhynchos.
The amount of shared polymorphism between the studied duck species cannot be explained by large population sizes of the respective species only. We suggest extraordinary and evolutionarily sustained hybridisation rates as drivers of ongoing gene pool mixing. Gene flow continues and will allow the transfer of genetic material among duck species. At present, extensive hybridisation still occurs. The genetic compatibility of different duck species, combined with mixed effects of genetically determined and imprinted mate choice leads to speciation reversals  despite genotypically and morphologically defined species boundaries. Present-day occurrence of Anas platyrhynchos in large numbers and wide geographical extent may even drive some of their close relatives to extinction by hybridisation . This is a major concern in many parts of the world, especially where Anas platyrhynchos is not indigenous . Many species of the genus Anas are hard to fit into the biological species concept because their evolution has rather led to a superspecies complex with discernable lineages. Besides the five dabbling duck species studied here, it is likely that many more of the ca. 40 Anas species are part of the global supra-population.
Besides conservation implications, this creates large evolutionary potential, comparable to bacteria, which are able to exchange genes among different species by horizontal gene transfer. Further, increasing effective population sizes into the millions may allow non-adaptive evolutionary processes to act, opening up additional degrees of evolutionary freedom . SNP-based analysis at hundreds of independent loci across the entire genome, as done here, may serve to re-evaluate long-standing puzzling patterns of speciation and hybridisation in several bird groups, such as other waterfowl, galliforms, hummingbirds and woodpeckers , as well as in many other organisms where species pairs exhibit unusually high levels of hybridisation.
In total, 212 individuals of Anas platyrhynchos obtained from nine localities representing Eurasian and North American populations were sampled and their blood was stored on FTA cards  at room temperature until DNA isolation. Numbers of samples with abbreviation codes in brackets: Eurasian samples were from Austria (25, ATHO), Estonia (22, EETA), Portugal (32, PTDJ), and three Russian localities: Yaroslavl (25, RUYA), Omsk (12, RUOM) and Tomsk (32, RUTO). North America was represented by Ontario (7, CALM), Manitoba (20, CARM) and Alaska (22, USMF). Preliminary multivariate clustering of SNP genotypes (see below) positioned 15 of these individuals far outside the Anas platyrhynchos species cluster, sometimes well within the clusters of other duck species (Additional file 3). We discarded these 15 individuals as mislabelled because they showed obvious deviation from their putative genotypic species cluster. Details are available in Additional file 4.
A set of 67 samples from other duck species were obtained world wide from various sources (hunting bags, live-trapped, zoos) and localities. Most often blood on FTA cards was used, sometimes other tissues stored in ethanol, and also previously isolated DNA from collections. The cross species testing was applied to ducks of the following Anas and Aythya species (numbers of samples and abbreviation code in brackets): Anas acuta (7; ANAC), Anas crecca (9; ANCR), Anas penelope (14; ANPE), Anas strepera (10; ANST), Aythya fuligula (17; AYFU) and one F1 hybrid between Anas acuta and Anas platyrhynchos (ANACPLA). Using the same procedure as with the Anas platyrhynchos set, we identified nine of these samples as apparently mislabelled. These were excluded from all subsequent analyses (Additional file 5). Details are available in Additional file 6.
DNA isolation was done using the Gentra Systems Puregene DNA purification Kit according to the manufacturer's instructions, with modifications when handling of FTA cards. Appropriate amounts of tissue or blood on FTA cards were digested with 9 μg Proteinase K (Sigma) in Cell Lysis Solution (Gentra Systems) at 65°C over night, or longer in case of some tissues. Proteins were subsequently precipitated with Protein Precipitation Solution (Gentra Systems) and spun down together with the FTA card material. DNA from the supernatant was precipitated with isopropanol and washed twice with 70% ethanol. DNA quantity and purity were measured using the Nanodrop ND1000. Samples with 260/280 nm absorption ratios less than 1.8 were purified again.
We used Illumina's GoldenGate Genotyping assay, on the Illumina BeadXpress. The marker set consisted of 384 SNPs  ("mallard 384 SNP set"), which are numbered according to their dbSNP accession numbers from ss263068950 (SNP 0) to ss263069333 (SNP 383). Raw genotyping results were analysed in GenomeStudio (Illumina), and SNP clusters adjusted by hand. The respective OPA (oligo pooled assay) and cluster files can be found online with this paper (Additional file 7 and Additional file 8).
SNP set evaluation
We assessed technical and biological properties of the SNP set in Anas platyrhynchos:
i) Minor allele frequencies and heterozygosity
For each of the nine localities we counted the occurrences of each of the two alleles. The count of the allele occurring less frequently (minor allele) was divided by the total number of alleles, giving the population wide frequency of minor alleles per locus (minor allele frequency, MAF). Additionally, we counted heterozygote individuals as a fraction of all individuals (observed heterozygosity, Hobs).
ii) Hardy-Weinberg equilibrium
Each locus was tested for deviation from Hardy-Weinberg equilibrium in each locality with the software Arlequin 188.8.131.52  using the analog to Fisher's exact test for arbitrary table size  (1,000,000 Markov chain steps, 100,000 dememorisation steps).
iii) Linkage disequilibrium
Per locality, pairs of SNP loci were tested for presence of linkage disequilibrium (LD) in Arlequin. The implemented likelihood-ratio test  employs the EM algorithm  to infer haplotypes from unphased genotypic data to test for statistical significance of LD. Repeated use of a SNP in multiple statistical tests requires a correction of the significance level α. In our 364 SNP data set each SNP is involved in 66066 pairwise tests, significance levels for LD are thus Bonferroni corrected.
iv) Physical SNP locations inferred from chicken genome
We searched the SNP positions on the chicken genome (from Kraus et al. ) in Ensembl  for chicken gene information using Bioconductor  with biomaRt in R . We recorded if a SNP was situated in a gene, or even intron.
Persistence times of SNPs
The equation for mean persistence time t(p) is a combination of the time to loss and to fixation [72, 73]. It can be written as -4N e [(1-p) ln(1-p) + pln(p)] where p denotes the initial MAF and Ne the effective population size (for derivation see ref. , page 112, eqn. 3.10). To calculate the persistence time t(p) of a SNP, an estimate of the effective population size (Ne) from the census population size (Nc) is thus required. Estimates of current census population sizes (Nc) of the investigated duck species were taken from the BirdLife species fact sheets . Upper estimates were used when a range was given. The ratio between Ne and Nc for species of dabbling ducks has to our knowledge not been studied, but it is probably fairly low as most census estimates are based on winter counts made several months before the breeding season starts and most mortality may occur before breeding . Further, dabbling ducks are generally r-selected and their population sizes fluctuate greatly by swift responses to benign and detrimental conditions [76, 77], with Ne being dominated by the smallest values . Estimated Ne: Nc ratios from white-winged wood ducks (Asarcornis scutulata, formerly Cairina scutulata) range between 0.052 (genetic measurements) and 0.094 (demographic measurements) . Thus, we use a ratio of 0.1 as a conservative estimate (on the high side) for the ratio of Ne to Nc.
The generation time has been set to one year for clarity. As mentioned above, many individuals do not reproduce at all, and those that do are in the vast majority of the cases first-years . The actual generation time value should lie in the range between 1.1 and 1.2 years, and this has no effect on our interpretations.
Interspecific, genetic admixture
A Bayesian genetic clustering algorithm as implemented in the software STRUCTURE  (version 2.3.3) was used to test for genetic admixture, i.e., the incorporation of genes from one discrete population/species into another. Two datasets were analysed: i) all Anas platyrhynchos and other duck species (the same individuals as analysed by PCA, see Figure 2c); ii) only the other species (cf. Figure 2b) plus the putative hybrid between Anas platyrhynchos and Anas acuta. A value of K = 6 simulated clusters (as many clusters as species) was chosen in the analysis of all ducks (i), and consequently K = 5 when Anas platyrhynchos was excluded (ii). Default settings were used with the admixture model of STRUCTURE, run for 300,000 steps (the first 100,000 discarded as burn-in). Additionally, we compared the results of the STRUCTURE analysis with those of the program InStruct  which is designed to perform the same analysis as STRUCTURE but not depending on Hardy-Weinberg equilibrium. The same datasets and settings were used, including the default settings, with the same values for K. Mode 1 - "infer population structure only with admixture" - in InStruct was chosen because it is most comparable to the program STRUCTURE as explained in its manual. The dataset containing only non-Anas platyrhynchos ducks (K = 5) was also run for the same amount of iteration steps. The larger dataset, all ducks combined (K = 6), was run substantially longer because the Markov chain converged very slowly (2,000,000 steps, of which 1,000,000 were discarded as burn-in).
Multivariate genetic clustering of genotypes
We tested for genetic similarity of individuals using principal component analysis (PCA) on their genotypes with the program smartpca from the Eigenstrat package  with default settings, but outlier removal switched off. The analysis was repeated for every new subset of the data.
We are thankful to the following persons and institutions for providing Anas platyrhynchos samples: Ernst Niedermayer, Hans Jörg Damm (Stiftung Fürst Liechtenstein, Austria), David Rodrigues (Escola Superior Agrária de Coimbra, Portugal), Brandt Meixell, Danielle Mondloch, Jonathan Runstadler (University of Alaska Fairbanks, USA), V.N. Stepanov, O. Tutenkov, V.I. Zalogin, Sergey Gashkov, Sergey Fokin (State Informational-Analytical Centre of Game Animal and Environment of Hunting Department of Russia), Urmas Võro, David Lamble, Garry Grigg, Aaron Everingham, Holly Middleton (Simon Fraser University, Vancouver, Canada). Rolik Grzegorz (Zoo Opole, Poland), Magnus Hellström (Ottenby Bird Observatory, Sweden), Michael Wink, Javier Gonzales (University of Heidelberg, Germany), Dirk Ullrich (Alpenzoo Innsbruck, Austria), Kamil Čihák (Zoo Dvur Kralove, Czech Republic), Marina Euler (Tierpark Lange Erlen, Switzerland), Sascha Knauf (Opel Zoo, Germany), Yang Liu (University of Bern, Switzerland), Mathieu Boos (CNRS Strasbourg, France), Crystal Matthews (Virginia Aquarium, USA), Timm Spretke (Zoologischer Garten Halle, Germany) and Valery Buzun provided samples from duck species other than Anas platyrhynchos. Technical assistance with genotyping was provided by Bert Dibbits. Daniël Goedbloed helped with the software package Eigenstrat. We thank Michael Turelli and Carlo Dietl for discussions. Javier Gonzales provided unpublished data on divergence times of duck species, and Brian Cade helped with statistics. The WWT, Slimbridge, UK, provided drawings for Figure 1. This work was financially supported by the KNJV (Royal Netherlands Hunters Association), the Dutch Ministry of Agriculture, the Faunafonds and the Stichting de Eik trusts (both in The Netherlands) and the Swedish Environmental Protection Agency, grants V-220-08 and V-205-09.
- Mayr E: Systematics and the origin of Species. 1942, New York: Columbia University PressGoogle Scholar
- Mallet J: A species definition for the modern synthesis. Trends Ecol Evol. 1995, 10: 294-299. 10.1016/0169-5347(95)90031-4.PubMedView ArticleGoogle Scholar
- Dobzhansky T: Studies on hybrid sterility. II Localization of sterility factors in Drosophila pseudoobscura hybrids. Genetics. 1936, 21: 113-135.PubMedGoogle Scholar
- Muller HJ: Isolating mechanisms, evolution and temperature. Temperature, Evolution, Development. Edited by: Dobzhansky T. 1942, Lancaster, PA: JaquesCattell Press, 6: 71-125. [Dobzhansky T (Series Editor): Biological Symposia: A Series of Volumes Devoted to Current Symposia in the Field of Biology]Google Scholar
- Gourbiére S, Mallet J: Are species real? The shape of the species boundary with exponential failure, reinforcement, and the "missing snowball". Evolution. 2010, 64: 1-24. 10.1111/j.1558-5646.2009.00844.x.PubMedView ArticleGoogle Scholar
- Prager EM, Wilson AC: Slow evolutionary loss of the potential for interspecific hybridization in birds: a manifestation of slow regulatory evolution. Proc Natl Acad Sci USA. 1975, 72: 200-204. 10.1073/pnas.72.1.200.PubMedPubMed CentralView ArticleGoogle Scholar
- Grant PR, Grant BR: Hybridization of bird species. Science. 1992, 256: 193-197. 10.1126/science.256.5054.193.PubMedView ArticleGoogle Scholar
- Grant PR, Grant BR: Genetics and the origin of bird species. Proc Natl Acad Sci USA. 1997, 94: 7768-7775. 10.1073/pnas.94.15.7768.PubMedPubMed CentralView ArticleGoogle Scholar
- Morin PA, Luikart G, Wayne RK: SNPs in ecology, evolution and conservation. Trends Ecol Evol. 2004, 19: 208-216. 10.1016/j.tree.2004.01.009.View ArticleGoogle Scholar
- Seehausen O: Conservation: Losing biodiversity by reverse speciation. Curr Biol. 2006, 16: R334-R337.PubMedView ArticleGoogle Scholar
- Seehausen O, Takimoto G, Roy D, Jokela J: Speciation reversal and biodiversity dynamics with hybridization in changing environments. Mol Ecol. 2008, 17: 30-44. 10.1111/j.1365-294X.2007.03529.x.PubMedView ArticleGoogle Scholar
- Tubaro PL, Lijtmaer DA: Hybridization patterns and the evolution of reproductive isolation in ducks. Biol Jour Linn Soc. 2002, 77: 193-200. 10.1046/j.1095-8312.2002.00096.x.View ArticleGoogle Scholar
- Scherer S, Hilsberg T: Hybridisierung und Verwandtschaftsgrade innerhalb der Anatidae - eine systematische und evolutionstheoretische Betrachtung [in German with English summary]. J Ornith. 1982, 123: 357-380. 10.1007/BF01643271.View ArticleGoogle Scholar
- Kraus RHS, Kerstens HHD, van Hooft P, Crooijmans RPMA, van Der Poel JJ, Elmberg J, Vignal A, Huang Y, Li N, Prins HHT, Groenen MAM: Genome wide SNP discovery, analysis and evaluation in mallard (Anas platyrhynchos). BMC Genomics. 2011, 12: 150-10.1186/1471-2164-12-150.PubMedPubMed CentralView ArticleGoogle Scholar
- Kraus RHS, Zeddeman A, Sartakov D, Soloviev SA, Ydenberg RC, Prins HHT: Evolution and connectivity in the world-wide migration system of the mallard: Inferences from mitochondrial DNA. BMC Genet. 2011, 12: 99-PubMedPubMed CentralView ArticleGoogle Scholar
- Michelizzi VN, Wu X, Dodson MV, Michal JJ, Zambrano-Varon J, McLean DJ, Jiang Z: A global view of 54,001 single nucleotide polymorphisms (SNPs) on the Illumina BovineSNP50 BeadChip and their transferability to Water Buffalo. Int J Biol Sci. 2011, 7: 18-27.PubMed CentralView ArticleGoogle Scholar
- Megens H-J, Crooijmans RPMA, Bastiaansen JWM, Kerstens HHD, Coster A, Jalving R, Vereijken A, Silva P, Muir WM, Cheng HH, et al: Comparison of linkage disequilibrium and haplotype diversity on macro- and microchromosomes in chicken. BMC Genet. 2009, 10: 86-PubMedPubMed CentralView ArticleGoogle Scholar
- Miller JM, Poissant J, Kijas JW, Coltman DW: A genome-wide set of SNPs detects population substructure and long range linkage disequilibrium in wild sheep. Mol Ecol Res. 2011, 11: 314-10.1111/j.1755-0998.2010.02918.x.View ArticleGoogle Scholar
- Charlesworth B, Bartolomé C, Noël V: The detection of shared and ancestral polymorphisms. Genet Res. 2005, 86: 149-157. 10.1017/S0016672305007743.PubMedView ArticleGoogle Scholar
- Hedges SB, Dudley J, Kumar S: TimeTree: A public knowledge-base of divergence times among organisms. Bioinformatics. 2006, 22: 2971-2972. 10.1093/bioinformatics/btl505.PubMedView ArticleGoogle Scholar
- Ramos-Onsins SE, Stranger BE, Mitchell-Olds T, Aguadé M: Multilocus analysis of variation and speciation in the closely related species Arabidopsis halleri and A. lyrata. Genetics. 2004, 166: 373-388. 10.1534/genetics.166.1.373.PubMedPubMed CentralView ArticleGoogle Scholar
- Gonzalez J, Düttmann H, Wink M: Phylogenetic relationships based on two mitochondrial genes and hybridization patterns in Anatidae. J Zool. 2009, 279: 310-318. 10.1111/j.1469-7998.2009.00622.x.View ArticleGoogle Scholar
- Hartl DL, Clark AG: Principles of Population Genetics. 2007, Sunderland, MA, USA: Sinauer Associates, 4Google Scholar
- Clark AG: Neutral behavior of shared polymorphism. Proc Natl Acad Sci USA. 1997, 94: 7730-7734. 10.1073/pnas.94.15.7730.PubMedPubMed CentralView ArticleGoogle Scholar
- Mayr E: Birds collected during the Whitney South Sea Expedition. XII. Notes on Halcyon chloris and some of its subspecies. Am Mus Novit. 1931, 469: 1-10.Google Scholar
- Rensch B: Grenzfälle von Rasse und Art [in German]. Journ f Ornith. 1928, 76: 222-231. 10.1007/BF01923570.View ArticleGoogle Scholar
- Urban EK, Fry CH, Keith S: Introduction to "Birds of Africa". Birds of Africa, Vol II. Edited by: Urban EK, Fry CH, Keith S. 1986, London: Academic Press, Inc., II: xi-xvi.Google Scholar
- Amadon D: The superspecies concept. Syst Zool. 1966, 15: 245-249.View ArticleGoogle Scholar
- Haffer J: Superspecies and species limits in vertebrates. Z Zool Syst Evol. 1986, 24: 169-190.View ArticleGoogle Scholar
- Kiriakoff SG: On the nomenclature of the superspecies. Syst Zool. 1967, 16: 281-282.View ArticleGoogle Scholar
- Mayr E, Short LL: Species Taxa of North-American Birds A Contribution to Comparative Systematics. 1970, Cambridge, Mass.: Publications of the Nuttall Ornithological Club, No 9Google Scholar
- Dubois A: New proposals for naming lower-ranked taxa within the frame of the International Code of Zoological Nomenclature. C R Biol. 2006, 329: 823-840. 10.1016/j.crvi.2006.07.003.PubMedView ArticleGoogle Scholar
- Avise JC, Ankney CD, Nelson WS: Mitochondrial gene trees and the evolutionary relationship of Mallard and Black Ducks. Evolution. 1990, 44: 1109-1119. 10.2307/2409570.View ArticleGoogle Scholar
- Kulikova IV, Zhuravlev YN, McCracken KG: Asymmetric hybridization and sex-biased gene flow between Eastern Spot-billed Ducks (Anas zonorhyncha) and Mallards (A. platyrhynchos) in the Russian Far East. Auk. 2004, 121: 930-949. 10.1642/0004-8038(2004)121[0930:AHASGF]2.0.CO;2.View ArticleGoogle Scholar
- Peters JL, McCracken KG, Zhuravlev YN, Lu Y, Wilson RE, Johnson KP, Omland KE: Phylogenetics of wigeons and allies (Anatidae: Anas): The importance of sampling multiple loci and multiple individuals. Mol Phylogenet Evol. 2005, 35: 209-224.PubMedView ArticleGoogle Scholar
- Peters JL, Zhuravlev Y, Fefelov I, Logie A, Omland KE: Nuclear loci and coalescent methods support ancient hybridization as cause of mitochondrial paraphyly between gadwall and falcated duck (Anas spp.). Evolution. 2007, 61: 1992-2006. 10.1111/j.1558-5646.2007.00149.x.PubMedView ArticleGoogle Scholar
- Pritchard JK, Stephens M, Donnelly P: Inference of population structure using multilocus genotype data. Genetics. 2000, 155: 945-959.PubMedPubMed CentralGoogle Scholar
- Gao H, Williamson S, Bustamante CD: A Markov chain Monte Carlo approach for joint inference of population structure and inbreeding rates from multilocus genotype data. Genetics. 2007, 176: 1635-1651. 10.1534/genetics.107.072371.PubMedPubMed CentralView ArticleGoogle Scholar
- Whitlock MC, Barton NH: The effective size of a subdivided population. Genetics. 1997, 146: 427-441.PubMedPubMed CentralGoogle Scholar
- Brodin A, Haas F: Speciation by perception. Anim Behav. 2006, 72: 139-146.View ArticleGoogle Scholar
- Immelmann K: Ecological significance of imprinting and early learning. Annu Rev Ecol Syst. 1975, 6: 15-37. 10.1146/annurev.es.06.110175.000311.View ArticleGoogle Scholar
- Cooke F, Cooch FG: The genetics of polymorphism in the goose Anser caerulescens. Evolution. 1968, 22: 289-300. 10.2307/2406528.View ArticleGoogle Scholar
- Zachos J, Pagani H, Sloan L, Thomas E, Billups K: Trends, rhythms, and aberrations in global climate 65 Ma to present. Science. 2001, 292: 686-693. 10.1126/science.1059412.PubMedView ArticleGoogle Scholar
- Jiménez-Moreno G, Fauquette S, Suc JP: Vegetation, climate and palaeoaltitude reconstructions of the Eastern Alps during the Miocene based on pollen records from Austria, Central Europe. J Biogeography. 2008, 35: 1638-1649. 10.1111/j.1365-2699.2008.01911.x.View ArticleGoogle Scholar
- Böhme M, Ilg A, Winklhofer M: Late Miocene "washhouse" climate in Europe. Earth Planet Sc Lett. 2008, 275: 393-401. 10.1016/j.epsl.2008.09.011.View ArticleGoogle Scholar
- Utescher T, Mosbrugger V, Ashraf AR: Terrestrial climate evolution in northwest Germany over the last 25 million years. Palaios. 2000, 15: 430-449.View ArticleGoogle Scholar
- Bruch AA, Uhl D, Mosbrugger V: Miocene climate in Europe - Patterns and evolution. A first synthesis of NECLIME. Palaeogeogr Palaeocl. 2007, 253: 1-7. 10.1016/j.palaeo.2007.03.030.View ArticleGoogle Scholar
- Harzhauser M, Piller WE: Benchmark data of a changing sea - Palaeogeography, Palaeobiogeography and events in the Central Paratethys during the Miocene. Palaeogeogr Palaeocl. 2007, 253: 8-31. 10.1016/j.palaeo.2007.03.031.View ArticleGoogle Scholar
- Harzhauser M, Latal C, Piller WE: The stable isotope archive of Lake Pannon as a mirror of Late Miocene climate change. Palaeogeogr Palaeocl. 2007, 249: 335-350. 10.1016/j.palaeo.2007.02.006.View ArticleGoogle Scholar
- Magyar I, Geary DH, Müller P: Paleogeographic evolution of the Late Miocene Lake Pannon in Central Europe. Palaeogeogr Palaeocl. 1999, 147: 151-167. 10.1016/S0031-0182(98)00155-2.View ArticleGoogle Scholar
- Cerling TE, Harris JM, MacFadden BJ, Leakey MG, Quade J, Eisenmann V, Ehleringer JR: Global vegetation change through the Miocene/Pliocene boundary. Nature. 1997, 389: 153-10.1038/38229.View ArticleGoogle Scholar
- Agustí J: Sanz de Siria A, Garcés M: Explaining the end of the hominoid experiment in Europe. J Hum Evol. 2003, 45: 145-153. 10.1016/S0047-2484(03)00091-5.PubMedView ArticleGoogle Scholar
- Zhisheng A, Kutzbach JE, Prell WL, Porter SC: Evolution of Asian monsoons and phased uplift of the Himalaya - Tibetan plateau since Late Miocene times. Nature. 2001, 411: 62-66. 10.1038/35075035.PubMedView ArticleGoogle Scholar
- Manegold A, Brink JS: Descriptions and palaeoecological implications of bird remains from the Middle Pleistocene of Florisbad, South Africa. Paläontol Z. 2010, 85: 19-32.View ArticleGoogle Scholar
- Worthy TH: Pliocene waterfowl (Aves: Anseriformes) from South Australia and a new genus and species. Emu. 2008, 108: 153-165. 10.1071/MU07063.View ArticleGoogle Scholar
- Olson SL: The identity of the fossil ducks described from Australia by C.W. De Vis. Emu. 1977, 77: 127-131. 10.1071/MU9770127.View ArticleGoogle Scholar
- Worthy TH: Descriptions and phylogenetic relationships of two new genera and four new species of Oligo-Miocene waterfowl (Aves: Anatidae) from Australia. Zool J Linn Soc. 2009, 156: 411-454. 10.1111/j.1096-3642.2008.00483.x.View ArticleGoogle Scholar
- Mlíkovský J: Cenozoic birds of the world Part 1: Europe. 2002, Prague, Czech Republic: Ninox PressGoogle Scholar
- Seiger MB: A computer simulation study of influence of imprinting on population structure. Am Nat. 1967, 101: 47-57. 10.1086/282468.View ArticleGoogle Scholar
- Schutz F: Objektfixierung geschlechtlicher Reaktionen bei Anatiden und Hühnern [in German]. Naturwissenschaften. 1963, 50: 624-625. 10.1007/BF00632393.View ArticleGoogle Scholar
- Mank JE, Carlson JE, Brittingham MC: A century of hybridization: Decreasing genetic distance between American black ducks and mallards. Conserv Genet. 2004, 5: 395-403.View ArticleGoogle Scholar
- Rhymer JM: Extinction by hybridization and introgression in anatine ducks. Acta Zool Sin. 2006, 52 (Supplement): 583-585.Google Scholar
- Lynch M: The frailty of adaptive hypotheses for the origins of organismal complexity. Proc Natl Acad Sci USA. 2007, 104: 8597-8604. 10.1073/pnas.0702207104.PubMedPubMed CentralView ArticleGoogle Scholar
- Smith LM, Burgoyne LA: Collecting, archiving and processing DNA from wildlife samples using FTA® databasing paper. BMC Ecol. 2004, 4: 4-10.1186/1472-6785-4-4.PubMedPubMed CentralView ArticleGoogle Scholar
- Excoffier L, Lischer HEL: Arlequin suite ver 3.5: A new series of programs to perform population genetics analyses under Linux and Windows. Mol Ecol Res. 2010, 10: 564-567. 10.1111/j.1755-0998.2010.02847.x.View ArticleGoogle Scholar
- Guo SW, Thompson EA: Performing the exact Test of Hardy-Weinberg proportion for multiple alleles. Biometrics. 1992, 48: 361-372. 10.2307/2532296.PubMedView ArticleGoogle Scholar
- Excoffier L, Slatkin M: Incorporating genotypes of relatives into a test of linkage disequilibrium. Am J Hum Genet. 1998, 62: 171-180. 10.1086/301674.PubMedPubMed CentralView ArticleGoogle Scholar
- Excoffier L, Slatkin M: Maximum-likelihood estimation of molecular haplotype frequencies in a diploid population. Mol Biol Evol. 1995, 12: 921-927.PubMedGoogle Scholar
- Birney E, Andrews TD, Bevan P, Caccamo M, Chen Y, Clarke L, Coates G, Cuff J, Curwen V, Cutts T, et al: An overview of Ensembl. Genome Res. 2004, 14: 925-928. 10.1101/gr.1860604.PubMedPubMed CentralView ArticleGoogle Scholar
- Gentleman RC, Carey VJ, Bates DM, Bolstad B, Dettling M, Dudoit S, Ellis B, Gautier L, Ge Y, Gentry J, et al: Bioconductor: open software development for computational biology and bioinformatics. Genome Biol. 2004, 5: R80-10.1186/gb-2004-5-10-r80.PubMedPubMed CentralView ArticleGoogle Scholar
- R Development Core Team: R: A language and environment for statistical computing. 2009, Vienna, Austria: R Foundation for Statistical Computing, [http://www.R-project.org]Google Scholar
- Kimura M, Ohta T: Average number of generations until extinction of an individual mutant gene in a finite population. Genetics. 1969, 63: 701-709.PubMedPubMed CentralGoogle Scholar
- Kimura M, Ohta T: Average number of generations until fixation of a mutant gene in a finite population. Genetics. 1969, 61: 763-771.PubMedPubMed CentralGoogle Scholar
- BirdLife International: Species factsheets. 2010, Downloaded from http://www.birdlife.org on 20/10/2010Google Scholar
- Gunnarsson G, Elmberg J, Dessborn L, Jonzén N, Pöysä H, Valkama J: Survival estimates, mortality patterns, and population growth of Fennoscandian mallards Anas platyrhynchos. Ann Zool Fennici. 2008, 45: 483-495.View ArticleGoogle Scholar
- Nudds TD: Niche dynamics and organization of waterfowl guilds in variable environments. Ecology. 1983, 64: 319-330. 10.2307/1937079.View ArticleGoogle Scholar
- Patterson JH: Can ducks be managed by regulation? Experiences in Canada. Trans North Am Wildl and Nat Resour Conf. 1979, 44: 130-139.Google Scholar
- Tomlinson C, Mace GM, Black JM, Hewston N: Improving the management of a highly inbred species: the case of the white-winged wood duck Cairina scutulata in captivity. Wildfowl. 1991, 42: 123-133.Google Scholar
- Price AL, Patterson NJ, Plenge RM, Weinblatt ME, Shadick NA, Reich D: Principal components analysis corrects for stratification in genome-wide association studies. Nat Genet. 2006, 38: 904-909. 10.1038/ng1847.PubMedView ArticleGoogle 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.