Determining evolutionary origin and phylogenetic relationships of mallard-like ducks of Oceania, greater Indonesia, and the Philippines with ddRAD-seq data

AIM

We aim to determine the evolutionary origins and population genetics of mallard-like ducks of Oceania, greater Indonesia, and the Philippines.

LOCATION

Oceania, greater Indonesia, and the Philippines.

TAXON

Mallard (Anas platyrhynchos), Pacific black duck (A. superciliosa spp.), and Philippine duck (A. luzonica) METHODS: Thousands of nuclear ddRAD-seq loci and the mitochondrial DNA control region were assayed across individuals representative of each species' range. We assessed population structure and phylogenetic relationships, as well as estimated demographic histories to reconstruct the biogeographical history of each species.

RESULTS

Philippine and Pacific black ducks represent unique genetic lineages that diverged from the mallard 1-2 million years ago. We find no support for the Philippine duck representing a hybrid species as once posited; however, their low levels of genetic diversity requires further attention. We find a lack of substructure among Philippine ducks. However, we found pronounced differentiation between subspecies of Pacific black ducks, especially between A. s. superciliosa from New Zealand and A. s. rogersi from Australia, Papua New Guinea, and Timor-Leste, Indonesia. Anas superciliosa pelewensis gave mixed results; individuals from the Solomon Islands were differentiated from the other subspecies, but those from the island of Aunu'u, American Samoa, were genetically more similar to A. s. rogersi than A. s. pelewensis samples from the Solomon Islands. Finally, we find limited evidence of interspecific gene flow at evolutionary scales, and mallard introgression among contemporary samples.

MAIN CONCLUSIONS

Mallard-like ducks radiated across Oceania, greater Indonesia, and the Philippines within the last 2 million years. Only the Pacific black duck showed unique sub-structuring that largely followed known sub-species ranges, except for A. s. pelewensis. We posit that the high interrelatedness among Solomon Island samples suggests that their genetic distinctiveness may simply be the result of high levels of genetic drift. In contrast, we conclude that mainland Australian Pacific black ducks were the most likely source for the recent colonization of American Samoa. As a result, our findings suggest that either the A. s. pelewensis subspecies designations and/or its geographical range may require re-evaluation. Continued re-evaluation of evolutionary and taxonomic relationships is necessary when attempting to reconstruct and understand biogeographical histories, with important implications towards any attempts to implement conservation strategies.

simRestore: A decision-making tool for adaptive management of the native genetic status of wild populations

Anthropogenic hybridization, or higher and non-natural rates of gene flow directly and indirectly induced by human activities, is considered a significant threat to biodiversity. The primary concern for conservation is the potential for genomic extinction and loss of adaptiveness for native species due to the extensive introgression of non-native genes. To alleviate or reverse trends for such scenarios requires the direct integration of genomic data within a model framework for effective management. Towards this end, we developed the simRestore R program as a decision-making tool that integrates ecological and genomic information to simulate ancestry outcomes from optimized conservation strategies. In short, the program optimizes supplementation and removal strategies across generations until a set native genetic threshold is reached within the studied population. Importantly, in addition to helping with initial decision-making, simulations can be updated with the outcomes of ongoing efforts, allowing for the adaptive management of populations. After demonstrating functionality, we apply and optimize among actionable management strategies for the endangered Hawaiian duck for which the current primary threat is genetic extinction through ongoing anthropogenic hybridization with feral mallards. Simulations demonstrate that supplemental and removal efforts can be strategically tailored to move the genetic ancestry of Hawaii's hybrid populations towards Hawaiian duck without the need to completely start over. Further, we discuss ecological parameter sensitivity, including which factors are most important to ensure genetic outcomes (i.e. number of offspring). Finally, to facilitate use, the program is also available online as a Shiny Web application.

Population genetics and geographic origins of mallards harvested in northwestern Ohio

The genetic composition of mallards in eastern North America has been changed by release of domestically-raised, game-farm mallards to supplement wild populations for hunting. We sampled 296 hatch-year mallards harvested in northwestern Ohio, October-December 2019. The aim was to determine their genetic ancestry and geographic origin to understand the geographic extent of game-farm mallard introgression into wild populations in more westward regions of North America. We used molecular analysis to detect that 35% of samples were pure wild mallard, 12% were early generation hybrids between wild and game-farm mallards (i.e., F1-F3), and the remaining 53% of samples were assigned as part of a hybrid swarm. Percentage of individuals in our study with some form of hybridization with game-farm mallard (65%) was greater than previously detected farther south in the mid-continent (~4%), but less than the Atlantic coast of North America (~ 92%). Stable isotope analysis using δ2Hf suggested that pure wild mallards originated from areas farther north and west than hybrid mallards. More specifically, 17% of all Ohio samples had δ2Hf consistent with more western origins in the prairies, parkland, or boreal regions of the mid-continent of North America, with 55%, 35%, and 10% of these being genetically wild, hybrid swarm, and F3, respectively. We conclude that continued game-farm introgression into wild mallards is not isolated to the eastern population of mallards in North America, and may be increasing and more widespread than previously detected. Mallards in our study had greater incidence of game-farm hybridization than other locales in the mid-continent but less than eastern North American regions suggesting further need to understand game-farm mallard genetic variation and movement across the continent.

Genomic and morphological data shed light on the complexities of shared ancestry between closely related duck species

Causes for genomic and morphological similarities among recently radiated species are often multifaceted and are further convoluted among species that readily interbreed. Here, we couple genomic and morphological trait comparisons to test the extent that ancestry and gene flow explain the retention of mallard-like traits within a sister species, the Mexican duck. First, we confirm that these taxa remain genetically structured, and that Mexican ducks exhibit an isolation-by-distance pattern. Despite the assumption of wide-spread hybridization, we found only a few late-stage hybrids, all from the southwestern USA. Next, assessing 23 morphological traits, we developed a genetically-vetted morphological key that is > 97% accurate in distinguishing across sex-age cohorts of Mexican ducks, mallards, and hybrids. During key development, we determined that 25% of genetically pure, immature male Mexican ducks of the northern population naturally displayed mallard-like traits in their formative plumage. In fact, applying this key to 55 museum specimens, we identified that only four of the 14 specimens originally classified as phenotypic hybrids were truly hybrids. We discuss how genomic and morphological comparisons shed light into the mechanism(s) underlying the evolution of complex phenotypic traits in recent radiations, and how misunderstanding the true morphological diversity within Mexican ducks resulted in taxonomic revisions that hindered conservation efforts.

Phylogenomics, introgression, and demographic history of South American true toads (Rhinella)

The effects of genetic introgression on species boundaries and how they affect species' integrity and persistence over evolutionary time have received increased attention. The increasing availability of genomic data has revealed contrasting patterns of gene flow across genomic regions, which impose challenges to inferences of evolutionary relationships and of patterns of genetic admixture across lineages. By characterizing patterns of variation across thousands of genomic loci in a widespread complex of true toads (Rhinella), we assess the true extent of genetic introgression across species thought to hybridize to extreme degrees based on natural history observations and multi-locus analyses. Comprehensive geographic sampling of five large-ranged Neotropical taxa revealed multiple distinct evolutionary lineages that span large geographic areas and, at times, distinct biomes. The inferred major clades and genetic clusters largely correspond to currently recognized taxa; however, we also found evidence of cryptic diversity within taxa. While previous phylogenetic studies revealed extensive mito-nuclear discordance, our genetic clustering analyses uncovered several admixed individuals within major genetic groups. Accordingly, historical demographic analyses supported that the evolutionary history of these toads involved cross-taxon gene flow both at ancient and recent times. Lastly, ABBA-BABA tests revealed widespread allele sharing across species boundaries, a pattern that can be confidently attributed to genetic introgression as opposed to incomplete lineage sorting. These results confirm previous assertions that the evolutionary history of Rhinella was characterized by various levels of hybridization even across environmentally heterogeneous regions, posing exciting questions about what factors prevent complete fusion of diverging yet highly interdependent evolutionary lineages.

Spontaneous hybridization and introgression between walleye (Sander vitreus) and sauger (Sander canadensis) in two large reservoirs: Insights from genotyping by sequencing

Anthropogenic activities may facilitate undesirable hybridization and genomic introgression between fish species. Walleye (Sander vitreus) and sauger (Sander canadensis) are economically valuable freshwater species that can spontaneously hybridize in areas of sympatry. Levels of genomic introgression between walleye and sauger may be increased by modifications to waterbodies (e.g., reservoir development) and inadvertent propagation of hybrids in stocking programs. We used genotyping by sequencing (GBS) to examine 217 fish from two large reservoirs with mixed populations of walleye and sauger in Saskatchewan, Canada (Lake Diefenbaker, Tobin Lake). Analyses with 20,038 (r90) and 478 (r100) single nucleotide polymorphisms clearly resolved walleye and sauger, and classified hybrids with high confidence. F1, F2, and multigeneration hybrids were detected in Lake Diefenbaker, indicating potentially high levels of genomic introgression. In contrast, only F1 hybrids were detected in Tobin Lake. Field classification of fish was unreliable; 7% of fish were misidentified based on broad species categories. Important for activities such as brood stock selection, 12 of 173 (7%) fish field identified as pure walleye, and one of 24 (4%) identified as pure sauger were actually hybrids. In addition, two of 15 (13%) field-identified hybrids were actually pure walleye or sauger. We conclude that hybridization and introgression are occurring in Saskatchewan reservoirs and that caution is warranted when using these populations in stocking programs. GBS offers a powerful and flexible tool for examining hybridization without preidentification of informative loci, eliminating some of the key challenges associated with other marker types.

Phylogenomics reveals ancient and contemporary gene flow contributing to the evolutionary history of sea ducks (Tribe Mergini)

Insight into complex evolutionary histories continues to build through broad comparative phylogenomic and population genomic studies. In particular, there is a need to understand the extent and scale that gene flow contributes to standing genomic diversity and the role introgression has played in evolutionary processes such as hybrid speciation. Here, we investigate the evolutionary history of the Mergini tribe (sea ducks) by coupling multi-species comparisons with phylogenomic analyses of thousands of nuclear ddRAD-seq loci, including Z-sex chromosome and autosomal linked loci, and the mitogenome assayed across all extant sea duck species in North America. All sea duck species are strongly structured across all sampled marker types (pair-wise species Φ_ST > 0.2), with clear genetic assignments of individuals to their respective species, and phylogenetic relationships recapitulate known relationships. Despite strong species integrity, we identify at least 18 putative hybrids; with all but one being late generational backcrosses. Most interesting, we provide the first evidence that an ancestral gene flow event between long-tailed ducks (Clangula hyemalis) and true Eiders (Somateria spp.) not only moved genetic material into the former species, but likely generated a novel species - the Steller's eider (Polysticta stelleri) - via hybrid speciation. Despite generally low contemporary levels of gene flow, we conclude that hybridization has and continues to be an important process that shifts novel genetic variation between species within the tribe Mergini. Finally, we outline methods that permit researchers to contrast genomic patterns of contemporary versus ancestral gene flow when attempting to reconstruct potentially complex evolutionary histories.

Lousy grouse: Comparing evolutionary patterns in Alaska galliform lice to understand host evolution and host-parasite interactions

Understanding both sides of host-parasite relationships can provide more complete insights into host and parasite biology in natural systems. For example, phylogenetic and population genetic comparisons between a group of hosts and their closely associated parasites can reveal patterns of host dispersal, interspecies interactions, and population structure that might not be evident from host data alone. These comparisons are also useful for understanding factors that drive host-parasite coevolutionary patterns (e.g., codivergence or host switching) over different periods of time. However, few studies have compared the evolutionary histories between multiple groups of parasites from the same group of hosts at a regional geographic scale. Here, we used genomic data to compare phylogenomic and population genomic patterns of Alaska ptarmigan and grouse species (Aves: Tetraoninae) and two genera of their associated feather lice: Lagopoecus and Goniodes. We used whole-genome sequencing to obtain hundreds of genes and thousands of single-nucleotide polymorphisms (SNPs) for the lice and double-digest restriction-associated DNA sequences to obtain SNPs from Alaska populations of two species of ptarmigan. We found that both genera of lice have some codivergence with their galliform hosts, but these relationships are primarily characterized by host switching and phylogenetic incongruence. Population structure was also uncorrelated between the hosts and lice. These patterns suggest that grouse, and ptarmigan in particular, share habitats and have likely had historical and ongoing dispersal within Alaska. However, the two genera of lice also have sufficient dissimilarities in the relationships with their hosts to suggest there are other factors, such as differences in louse dispersal ability, that shape the evolutionary patterns with their hosts.

Assessing changes in genomic divergence following a century of human-mediated secondary contact among wild and captive-bred ducks

Along with manipulating habitat, the direct release of domesticated individuals into the wild is a practice used worldwide to augment wildlife populations. We test between possible outcomes of human-mediated secondary contact using genomic techniques at both historical and contemporary timescales for two iconic duck species. First, we sequence several thousand ddRAD-seq loci for contemporary mallards (Anas platyrhynchos) throughout North America and two domestic mallard types (i.e., known game-farm mallards and feral Khaki Campbell's). We show that North American mallards may well be becoming a hybrid swarm due to interbreeding with domesticated game-farm mallards released for hunting. Next, to attain a historical perspective, we applied a bait-capture array targeting thousands of loci in century-old (1842-1915) and contemporary (2009-2010) mallard and American black duck (Anas rubripes) specimens. We conclude that American black ducks and mallards have always been closely related, with a divergence time of ~600,000 years before present, and likely evolved through prolonged isolation followed by limited bouts of gene flow (i.e., secondary contact). They continue to maintain genetic separation, a finding that overturns decades of prior research and speculation suggesting the genetic extinction of the American black duck due to contemporary interbreeding with mallards. Thus, despite having high rates of hybridization, actual gene flow is limited between mallards and American black ducks. Conversely, our historical and contemporary data confirm that the intensive stocking of game-farm mallards during the last ~100 years has fundamentally changed the genetic integrity of North America's wild mallard population, especially in the east. It thus becomes of great interest to ask whether the iconic North American mallard is declining in the wild due to introgression of maladaptive traits from domesticated forms. Moreover, we hypothesize that differential gene flow from domestic game-farm mallards into the wild mallard population may explain the overall temporal increase in differentiation between wild black ducks and mallards, as well as the uncoupling of genetic diversity and effective population size estimates across time in our results. Finally, our findings highlight how genomic methods can recover complex population histories by capturing DNA preserved in traditional museum specimens.

Old divergence and restricted gene flow between torrent duck (Merganetta armata) subspecies in the Central and Southern Andes

AIM

To investigate the structure and rate of gene flow among populations of habitat-specialized species to understand the ecological and evolutionary processes underpinning their population dynamics and historical demography, including speciation and extinction.

LOCATION

Peruvian and Argentine Andes.

TAXON

Two subspecies of torrent duck (Merganetta armata).

METHODS

We sampled 156 individuals in Peru (M. a. leucogenis; Chillón River, n = 57 and Pachachaca River, n = 49) and Argentina (M. a. armata; Arroyo Grande River, n = 33 and Malargüe River, n = 17), and sequenced the mitochondrial DNA (mtDNA) control region to conduct coarse and fine-scale demographic analyses of population structure. Additionally, to test for differences between subspecies, and across genetic markers with distinct inheritance patterns, a subset of individuals (Peru, n = 10 and Argentina, n = 9) was subjected to partial genome resequencing, obtaining 4,027 autosomal and 189 Z-linked double-digest restriction-associated DNA sequences.

RESULTS

Haplotype and nucleotide diversities were higher in Peru than Argentina across all markers. Peruvian and Argentine subspecies showed concordant species-level differences (ΦST mtDNA = 0.82; ΦST autosomal = 0.30; ΦST Z chromosome = 0.45), including no shared mtDNA haplotypes. Demographic parameters estimated for mtDNA using IM and IMa2 analyses, and for autosomal markers using ∂a∂i (isolation-with-migration model), supported an old divergence (mtDNA = 600,000 years before present (ybp), 95% HPD range = 1.2 Mya to 200,000 ybp; and autosomal ∂a∂i = 782,490 ybp), between the two subspecies, characteristic of deeply diverged lineages. The populations were well-differentiated in Argentina but moderately differentiated in Peru, with low unidirectional gene flow in each country.

MAIN CONCLUSIONS

We suggest that the South American Arid Diagonal was preexisting and remains a current phylogeographic barrier between the ranges of the two torrent duck subspecies, and the adult territoriality and breeding site fidelity to the rivers define their population structure.

Coast to coast: High genomic connectivity in North American scoters

Dispersal shapes demographic processes and therefore is fundamental to understanding biological, ecological, and evolutionary processes acting within populations. However, assessing population connectivity in scoters (Melanitta sp.) is challenging as these species have large spatial distributions that span remote landscapes, have varying nesting distributions (disjunct vs. continuous), exhibit unknown levels of dispersal, and vary in the timing of the formation of pair bonds (winter vs. fall/spring migration) that may influence the distribution of genetic diversity. Here, we used double-digest restriction-associated DNA sequence (ddRAD) and microsatellite genotype data to assess population structure within the three North American species of scoter (black scoter, M. americana; white-winged scoter, M. deglandi; surf scoter, M. perspicillata), and between their European congeners (common scoter, M. nigra; velvet scoter, M. fusca). We uncovered no or weak genomic structure (ddRAD Φ ST < 0.019; microsatellite F ST < 0.004) within North America but high levels of structure among European congeners (ddRAD Φ ST > 0.155, microsatellite F ST > 0.086). The pattern of limited genomic structure within North America is shared with other sea duck species and is often attributed to male-biased dispersal. Further, migratory tendencies (east vs. west) of female surf and white-winged scoters in central Canada are known to vary across years, providing additional opportunities for intracontinental dispersal and a mechanism for the maintenance of genomic connectivity across North America. In contrast, the black scoter had relatively elevated levels of divergence between Alaska and Atlantic sites and a second genetic cluster found in Alaska at ddRAD loci was concordant with its disjunct breeding distribution suggestive of a dispersal barrier (behavioral or physical). Although scoter populations appear to be connected through a dispersal network, a small percentage (<4%) of ddRAD loci had elevated divergence which may be useful in linking areas (nesting, molting, staging, and wintering) throughout the annual cycle.

Identifying hybrids & the genomics of hybridization: Mallards & American black ducks of Eastern North America

Resolving evolutionary relationships and establishing population structure depends on molecular diagnosability that is often limited for closely related taxa. Here, we use 3,200 ddRAD-seq loci across 290 mallards, American black ducks, and putative hybrids to establish population structure and estimate hybridization rates. We test between traditional assignment probability and accumulated recombination events based analyses to assign hybrids to generational classes. For hybrid identification, we report the distribution of recombination events complements ADMIXTURE simulation by extending resolution past F4 hybrid status; however, caution against hybrid assignment based on accumulated recombination events due to an inability to resolve F1 hybrids. Nevertheless, both analyses suggest that there are relatively few backcrossed stages before a lineage's hybrid ancestry is lost and the offspring are effectively parental again. We conclude that despite high rates of observed interspecific hybridization between mallards and black ducks in the middle part of the 20th century, our results do not support the predicted hybrid swarm. Conversely, we report that mallard samples genetically assigned to western and non-western clusters. We indicate that these non-western mallards likely originated from game-farm stock, suggesting landscape level gene flow between domestic and wild conspecifics.

Ultraconserved elements (UCEs) illuminate the population genomics of a recent, high-latitude avian speciation event

Using a large, consistent set of loci shared by descent (orthologous) to study relationships among taxa would revolutionize among-lineage comparisons of divergence and speciation processes. Ultraconserved elements (UCEs), highly conserved regions of the genome, offer such genomic markers. The utility of UCEs for deep phylogenetics is clearly established and there are mature analytical frameworks available, but fewer studies apply UCEs to recent evolutionary events, creating a need for additional example datasets and analytical approaches. We used UCEs to study population genomics in snow and McKay's buntings (Plectrophenax nivalis and P. hyperboreus). Prior work suggested divergence of these sister species during the last glacial maximum (∼18-74 Kya). With a sequencing depth of ∼30× from four individuals of each species, we used a series of analysis tools to genotype both alleles, obtaining a complete dataset of 2,635 variable loci (∼3.6 single nucleotide polymorphisms/locus) and 796 invariable loci. We found no fixed allelic differences between the lineages, and few loci had large allele frequency differences. Nevertheless, individuals were 100% diagnosable to species, and the two taxa were different genetically (F ST = 0.034; P = 0.03). The demographic model best fitting the data was one of divergence with gene flow. Estimates of demographic parameters differed from published mtDNA research, with UCE data suggesting lower effective population sizes (∼92,500-240,500 individuals), a deeper divergence time (∼241,000 years), and lower gene flow (2.8-5.2 individuals per generation). Our methods provide a framework for future population studies using UCEs, and our results provide additional evidence that UCEs are useful for answering questions at shallow evolutionary depths.

Dissection by genomic and plumage variation of a geographically complex hybrid zone between two Australian non-sister parrot species, Platycercus adscitus and Platycercus eximius

The study of hybrid zones advances understanding of the speciation process, and approaches incorporating genomic data are increasingly used to draw significant conclusions about the impact of hybridisation. Despite the progress made, the complex interplay of factors that can lead to substantially variable hybridisation outcomes are still not well understood, and many systems and/or groups remain comparatively poorly studied. Our study aims to broaden the literature on avian hybrid zones, investigating a potentially geographically and temporally complex putative hybrid zone between two native Australian non-sister parrot species, the pale-headed and eastern rosellas (Platycercus adscitus and Platycercus eximius, respectively). We analysed six plumage traits and >1400 RADseq loci and detected hybrid individuals and an unexpectedly complex geographic structure. The hybrid zone is larger than previously described due to either observer bias or its movement over recent decades. It comprises different subregions where genetic and plumage signals of admixture vary markedly in their concordance. Evidence of contemporary hybridisation (later generation and backcrossed individuals) both within and beyond the previously defined zone, when coupled with a lack of F1 hybrids and differential patterns of introgression among potentially diagnostic loci, indicates a lack of post-zygotic barriers to gene flow between species. Despite ongoing gene flow, species boundaries are likely maintained largely by strong pre-mating barriers. These findings are discussed in detail and future avenues for research into this system are proposed, which would be of benefit to the speciation and hybrid zone literature.

Migration-Selection Balance Drives Genetic Differentiation in Genes Associated with High-Altitude Function in the Speckled Teal (Anas flavirostris) in the Andes

Local adaptation frequently occurs across populations as a result of migration-selection balance between divergent selective pressures and gene flow associated with life in heterogeneous landscapes. Studying the effects of selection and gene flow on the adaptation process can be achieved in systems that have recently colonized extreme environments. This study utilizes an endemic South American duck species, the speckled teal (Anas flavirostris), which has both high- and low-altitude populations. High-altitude speckled teal (A. f. oxyptera) are locally adapted to the Andean environment and mostly allopatric from low-altitude birds (A. f. flavirostris); however, there is occasional gene flow across altitudinal gradients. In this study, we used next-generation sequencing to explore genetic patterns associated with high-altitude adaptation in speckled teal populations, as well as the extent to which the balance between selection and migration have affected genetic architecture. We identified a set of loci with allele frequencies strongly correlated with altitude, including those involved in the insulin-like signaling pathway, bone morphogenesis, oxidative phosphorylation, responders to hypoxia-induced DNA damage, and feedback loops to the hypoxia-inducible factor pathway. These same outlier loci were found to have depressed gene flow estimates, as well as being highly concentrated on the Z-chromosome. Our results suggest a multifactorial response to life at high altitudes through an array of interconnected pathways that are likely under positive selection and whose genetic components seem to be providing an effective genomic barrier to interbreeding, potentially functioning as an avenue for population divergence and speciation.

A history of hybrids? Genomic patterns of introgression in the True Geese

Background

The impacts of hybridization on the process of speciation are manifold, leading to distinct patterns across the genome. Genetic differentiation accumulates in certain genomic regions, while divergence is hampered in other regions by homogenizing gene flow, resulting in a heterogeneous genomic landscape. A consequence of this heterogeneity is that genomes are mosaics of different gene histories that can be compared to unravel complex speciation and hybridization events. However, incomplete lineage sorting (often the outcome of rapid speciation) can result in similar patterns. New statistical techniques, such as the D-statistic and hybridization networks, can be applied to disentangle the contributions of hybridization and incomplete lineage sorting. We unravel patterns of hybridization and incomplete lineage sorting during and after the diversification of the True Geese (family Anatidae, tribe Anserini, genera Anser and Branta) using an exon-based hybridization network approach and taking advantage of discordant gene tree histories by re-sequencing all taxa of this clade. In addition, we determine the timing of introgression and reconstruct historical effective population sizes for all goose species to infer which demographic or biogeographic factors might explain the observed patterns of introgression.

Results

We find indications for ancient interspecific gene flow during the diversification of the True Geese and were able to pinpoint several putative hybridization events. Specifically, in the genus Branta, both the ancestor of the White-cheeked Geese (Hawaiian Goose, Canada Goose, Cackling Goose and Barnacle Goose) and the ancestor of the Brent Goose hybridized with Red-breasted Goose. One hybridization network suggests a hybrid origin for the Red-breasted Goose, but this scenario seems unlikely and it not supported by the D-statistic analysis. The complex, highly reticulated evolutionary history of the genus Anser hampered the estimation of ancient hybridization events by means of hybridization networks. The reconstruction of historical effective population sizes shows that most species showed a steady increase during the Pliocene and Pleistocene. These large effective population sizes might have facilitated contact between diverging goose species, resulting in the establishment of hybrid zones and consequent gene flow.

Conclusions

Our analyses suggest that the evolutionary history of the True Geese is influenced by introgressive hybridization. The approach that we have used, based on genome-wide phylogenetic incongruence and network analyses, will be a useful procedure to reconstruct the complex evolutionary histories of many naturally hybridizing species groups.

Electronic supplementary material

The online version of this article (doi:10.1186/s12862-017-1048-2) contains supplementary material, which is available to authorized users.

Hybridization between two cryptic filamentous brown seaweeds along the shore: analysing pre- and postzygotic barriers in populations of individuals with varying ploidy levels

We aimed to study the importance of hybridization between two cryptic species of the genus Ectocarpus, a group of filamentous algae with haploid-diploid life cycles that include the principal genetic model organism for the brown algae. In haploid-diploid species, the genetic structure of the two phases of the life cycle can be analysed separately in natural populations. Such life cycles provide a unique opportunity to estimate the frequency of hybrid genotypes in diploid sporophytes and meiotic recombinant genotypes in haploid gametophytes allowing the effects of reproductive barriers preventing fertilization or preventing meiosis to be untangle. The level of hybridization between E. siliculosus and E. crouaniorum was quantified along the European coast. Clonal cultures (568 diploid, 336 haploid) isolated from field samples were genotyped using cytoplasmic and nuclear markers to estimate the frequency of hybrid genotypes in diploids and recombinant haploids. We identified admixed individuals using microsatellite loci, classical assignment methods and a newly developed Bayesian method (XPloidAssignment), which allows the analysis of populations that exhibit variations in ploidy level. Over all populations, the level of hybridization was estimated at 8.7%. Hybrids were exclusively observed in sympatric populations. More than 98% of hybrids were diploids (40% of which showed signs of aneuploidy) with a high frequency of rare alleles. The near absence of haploid recombinant hybrids demonstrates that the reproductive barriers are mostly postzygotic and suggests that abnormal chromosome segregation during meiosis following hybridization of species with different genome sizes could be a major cause of interspecific incompatibility in this system.

Identifying hybridization and admixture using SNPs: application of the DArTseq platform in phylogeographic research on vertebrates

Next-generation sequencing (NGS) approaches are increasingly being used to generate multi-locus data for phylogeographic and evolutionary genetics research. We detail the applicability of a restriction enzyme-mediated genome complexity reduction approach with subsequent NGS (DArTseq) in vertebrate study systems at different evolutionary and geographical scales. We present two case studies using SNP data from the DArTseq molecular marker platform. First, we used DArTseq in a large phylogeographic study of the agamid lizard Ctenophorus caudicinctus, including 91 individuals and spanning the geographical range of this species across arid Australia. A low-density DArTseq assay resulted in 28 960 SNPs, with low density referring to a comparably reduced set of identified and sequenced markers as a cost-effective approach. Second, we applied this approach to an evolutionary genetics study of a classic frog hybrid zone (Litoria ewingii-Litoria paraewingi) across 93 individuals, which resulted in 48 117 and 67 060 SNPs for a low- and high-density assay, respectively. We provide a docker-based workflow to facilitate data preparation and analysis, then analyse SNP data using multiple methods including Bayesian model-based clustering and conditional likelihood approaches. Based on comparison of results from the DArTseq platform and traditional molecular approaches, we conclude that DArTseq can be used successfully in vertebrates and will be of particular interest to researchers working at the interface between population genetics and phylogenetics, exploring species boundaries, gene exchange and hybridization.

Naturally occurring hybrids of coral reef butterflyfishes have similar fitness compared to parental species

Hybridisation can produce evolutionary novelty by increasing fitness and adaptive capacity. Heterosis, or hybrid vigour, has been documented in many plant and animal taxa, and is a notable consequence of hybridisation that has been exploited for decades in agriculture and aquaculture. On the contrary, loss of fitness in naturally occurring hybrid taxa has been observed in many cases. This can have negative consequences for the parental species involved (wasted reproductive effort), and has raised concerns for species conservation. This study evaluates the relative fitness of previously documented butterflyfish hybrids of the genus Chaetodon from the Indo-Pacific suture zone at Christmas Island. Histological examination confirmed the reproductive viability of Chaetodon hybrids. Examination of liver lipid content showed that hybrid body condition was not significantly different from parent species body condition. Lastly, size at age data revealed no difference in growth rates and asymptotic length between hybrids and parent species. Based on the traits measured in this study, naturally occurring hybrids of Chaetodon butterflyfishes have similar fitness to their parental species, and are unlikely to supplant parental species under current environmental conditions at the suture zone. However, given sufficient fitness and ongoing genetic exchange between the respective parental species, hybrids are likely to persist within the suture zone.

Contrasting molecular and morphological evidence for the identification of an anomalous Buteo: a cautionary tale for hybrid diagnosis

An adult Buteo was found dead as a road-kill south of Sacramento, California, and was thought to represent the first state record of the eastern Red-shouldered Hawk (B. lineatus lineatus;). It is now a specimen in the Museum of Wildlife and Fisheries Biology (WFB 4816) at the University of California, Davis. We examined this specimen and found that many of its plumage characters differed from all other adult Red-shouldered Hawks examined, including nominate adults. Plumage markings and measurements were intermediate between Red-tailed Hawk (Buteo jamaicensis, ssp calurus) and Red-shouldered Hawk (ssp elegans), leading us to hypothesize that the bird was a hybrid. However, mtDNA sequences and nuDNA microsatellites proved definitively that the bird was a Red-shouldered Hawk, most likely of eastern origin. This case illustrates that apparent hybrids or apparent vagrants could be individuals with anomalous phenotypes caused by rare genetic variation or novel epigenetic effects.

An examination of species limits in the Aulacorhynchus "prasinus" toucanet complex (Aves: Ramphastidae)

The number of species recognized in Aulacorhynchus toucanets has varied tremendously over the past century. Revisors seem to disagree on whether head and bill coloration are useful indicators of species limits, especially in the A. "prasinus" complex. Using morphometrics, I tested the hypothesis that the major color-based subspecific groups of A. "prasinus" sensu lato are simply "cookie-cutter" (i.e., morphologically nearly identical) toucanets with different head and bill colorations. Univariate and multivariate analyses show that they are not simply morphological replicates of different colors: a complex array of morphometric similarities and dissimilarities occur between the major subspecific groups, and these variations differ between the sexes. Latitude and longitude had a small but significant association with female (but not male) PC1 and PC2. Hybridization and intergradation were also considered using plumage and bill characters as a surrogate to infer gene flow. Hybridization as indicated by phenotype appears to be substantial between A. "p." cyanolaemus and A. "p." atrogularis and nonexistent between other major groups, although from genetic evidence it is likely rare between A. "p." albivitta and A. "p." cyanolaemus. The congruence and complexities of the morphological and color changes occurring among these groups suggest that ecological adaptation (through natural selection) and social selection have co-occurred among these groups and that species limits are involved. Further, hybridization is not evident at key places, despite in many cases (hypothetical) opportunity for gene flow. Consequently, I recommend that this complex be recognized as comprising five biological species: A. wagleri, prasinus, caeruleogularis, albivitta, and atrogularis. Four of these also have valid subspecies within them, and additional work may eventually support elevation of some of these subspecies to full species. Species limits in South America especially need more study.

Recurrent hybridization and recent origin obscure phylogenetic relationships within the 'white-headed' gull (Larus sp.) complex

Species complexes that have undergone recent radiations are often characterized by extensive allele sharing due to recent ancestry and (or) introgressive hybridization. This can result in discordant evolutionary histories of genes and heterogeneous genomes, making delineating species limits difficult. Here we examine the phylogenetic relationships among a complex group of birds, the white-headed gulls (Aves: Laridae), which offer a unique window into the speciation process due to their recent evolutionary history and propensity to hybridize. Relationships were examined among 17 species (61 populations) using a multilocus approach, including mitochondrial and nuclear intron DNA sequences and microsatellite genotype information. Analyses of microsatellite and intron data resulted in some species-based groupings, although most species were not represented by a single cluster. Considerable allele and haplotype sharing among white-headed gull species was observed; no locus contained a species-specific clade. Despite this, our multilocus approach provided better resolution among some species than previous studies. Interestingly, most clades appear to correspond to geographic locality: our BEAST analysis recovered strong support for a northern European/Icelandic clade, a southern European/Russian clade, and a western North American/canus clade, with weak evidence for a high latitude clade spanning North America and northwestern Europe. This geographical structuring is concordant with behavioral observations of pervasive hybridization in areas of secondary contact. The extent of allele and haplotype sharing indicates that ecological and sexual selection are likely not strong enough to complete reproductive isolation within several species in the white-headed gull complex. This suggests that just a few genes are driving the speciation process.