ddRAD‐seq data reveal significant genome‐wide population structure and divergent genomic regions that distinguish the mallard and close relatives in North America
- 3 April 2019
- journal article
- research article
- Published by Wiley in Molecular Ecology
- Vol. 28 (10) , 2594-2609
Abstract
Recently evolved species typically share genetic variation across their genomes due to incomplete lineage sorting and/or ongoing gene flow. Given only subtle allele frequency differences at most loci and the expectation that divergent selection may affect only a tiny fraction of the genome, distinguishing closely related species based on multi‐locus data requires substantial genomic coverage. In this study, we used ddRAD‐seq to sample the genomes of five recently diverged, New World “mallards” (Anas spp.), a group of dabbling duck species characterized by diagnosable phenotypic differences but minimal genetic differentiation. With increased genomic sampling, we aimed to characterize population structure within this group and identify genomic regions that may have experienced divergent selection during speciation. We analyzed 3,017 autosomal ddRAD‐seq loci and 177 loci from the Z‐chromosome. In contrast to previous studies, the ddRAD‐seq data were sufficient to assign individuals to their respective species or subspecies and to generate estimates of gene flow in a phylogenetic framework. We find limited evidence of contemporary gene flow between the dichromatic mallard and several monochromatic taxa, but find evidence for historical gene flow between some monochromatic species pairs. We conclude that the overall genetic similarity of these taxa likely reflects retained ancestral polymorphism rather than recent and extensive gene flow. Thus, despite recurring cases of hybridization in this group, our results challenge the current dogma predicting the genetic extinction of the New World monochromatic dabbling ducks via introgressive hybridization with mallards. Moreover, ddRAD‐seq data were sufficient to identify previously unknown outlier regions across the Z‐chromosome and several autosomal chromosomes that may have been involved in the diversification of species in this recent radiation.Keywords
This publication has 90 references indexed in Scilit:
- RADSeq: next-generation population geneticsBriefings in Functional Genomics, 2010
- Speciation genetics: current status and evolving approachesPhilosophical Transactions of the Royal Society B: Biological Sciences, 2010
- A quasi-Newton acceleration for high-dimensional optimization algorithmsStatistics and Computing, 2009
- Reconstructing Indian population historyNature, 2009
- Fast model-based estimation of ancestry in unrelated individualsGenome Research, 2009
- The evolution of restricted recombination in sex chromosomesTrends in Ecology & Evolution, 2009
- Measurement of the human allele frequency spectrum demonstrates greater genetic drift in East Asians than in EuropeansNature Genetics, 2007
- PLINK: A Tool Set for Whole-Genome Association and Population-Based Linkage AnalysesThe American Journal of Human Genetics, 2007
- Rapid and cost-effective polymorphism identification and genotyping using restriction site associated DNA (RAD) markersGenome Research, 2006
- Mathematical model for studying genetic variation in terms of restriction endonucleases.Proceedings of the National Academy of Sciences, 1979