skip to main content

Title: Admixture mapping identifies introgressed genomic regions in North American canids

Hybrid zones typically contain novel gene combinations that can be tested by natural selection in a unique genetic context. Parental haplotypes that increase fitness can introgress beyond the hybrid zone, into the range of parental species. We used the Affymetrix canineSNPgenotyping array to identify genomic regions tagged by multiple ancestry informative markers that are more frequent in an admixed population than expected. We surveyed a hybrid zone formed in the last 100 years as coyotes expanded their range into eastern North America. Concomitant with expansion, coyotes hybridized with wolves and some populations became more wolflike, such that coyotes in the northeast have the largest body size of any coyote population. Using a set of 3102 ancestry informative markers, we identified 60 differentially introgressed regions in 44 canines across this admixture zone. These regions are characterized by an excess of exogenous ancestry and, in northeastern coyotes, are enriched for genes affecting body size and skeletal proportions. Further, introgressed wolf‐derived alleles have penetrated into SouthernUScoyote populations. Because no wolves currently exist in this area, these alleles are unlikely to have originated from recent hybridization. Instead, they probably originated from intraspecific gene flow or ancient admixture. We show that grey wolf and coyote admixture has far‐reaching effects and, in addition to phenotypically transforming admixed populations, allows for the differential movement of alleles from different parental species to be tested in new genomic backgrounds.

more » « less
Author(s) / Creator(s):
 ;  ;  ;  
Publisher / Repository:
Date Published:
Journal Name:
Molecular Ecology
Page Range / eLocation ID:
p. 2443-2453
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. Abstract

    Incompletely isolated species provide an opportunity to investigate the genetic mechanisms and evolutionary forces that maintain distinct species in the face of ongoing gene flow. Here, we use field surveys and reduced representation sequencing to characterize the patterns of reproductive isolation, admixture and genomic divergence between populations of the outcrossing wildflowerMimulus guttatusand selfingM. nasutus. Focusing on a single site where these two species have come into secondary contact, we find that phenological isolation is strong, although incomplete, and is likely driven by divergence in response to photoperiod. In contrast to previous field studies, which have suggested that F1‐hybrid formation might be rare, we discover patterns of genomic variation consistent with ongoing introgression. Strikingly, admixed individuals vary continuously from highly admixed to nearly pureM. guttatus, demonstrating ongoing hybridization and asymmetric introgression fromM. nasutusintoM. guttatus. Patterns of admixture and divergence across the genome show that levels of introgression are more variable than expected by chance. Some genomic regions show a reduced introgression, including one region that overlaps a critical photoperiodQTL, whereas other regions show elevated levels of interspecific gene flow. In addition, we observe a genome‐wide negative relationship between absolute divergence and the local recombination rate, potentially indicating natural selection againstM. nasutusancestry inM. guttatusgenetic backgrounds. Together, our results suggest thatMimulusspeciation is both ongoing and dynamic and that a combination of divergence in phenology and mating system, as well as selection against interspecific alleles, likely maintains these sympatric species.

    more » « less
  2. Coyotes are ubiquitous on the North American landscape as a result of their recent expansion across the continent. They have been documented in the heart of some of the most urbanized cities, such as Chicago, Los Angeles, and New York City. Here, we explored the genomic composition of 16 coyotes in the New York metropolitan area to investigate genomic demography and admixture for urban-dwelling canids in Queens County, New York. We identified moderate-to-high estimates of relatedness among coyotes living in Queens (r = 0.0–0.5) and adjacent neighborhoods, suggestive of a relatively small population. Although we found low background levels of domestic-dog ancestry across most coyotes in our sample (5%), we identified a male suspected to be a first-generation coyote–dog hybrid with 46% dog ancestry, as well as his two putative backcrossed offspring that carried approximately 25% dog ancestry. The male coyote–dog hybrid and one backcrossed offspring each carried two transposable element insertions that are associated with human-directed hypersociability in dogs and gray wolves. An additional, unrelated coyote with little dog ancestry also carried two of these insertions. These genetic patterns suggest that gene flow from domestic dogs may become an increasingly important consideration as coyotes continue to inhabit metropolitan regions. 
    more » « less
  3. Abstract

    Hybrid zones are geographic regions where isolating barriers between divergent populations are challenged by admixture. Identifying factors that facilitate or inhibit hybridization in sympatry can illuminate the processes that maintain those reproductive barriers. We analysed patterns of hybridization and phenotypic variation across two newly discovered hybrid zones between three subspecies of barn swallow (Hirundo rustica). These subspecies differ in ventral coloration and wing length, traits that are targets of sexual and natural selection, respectively, and are associated with genome‐wide differentiation in allopatry. We tested the hypothesis that the degree of divergence in these traits is associated with the extent of hybridization in secondary contact. We applied measures of population structure based on >23,000SNPs to confirm that named subspecies correspond to distinct genomic clusters, and assessed coincidence between geographic clines for ancestry and phenotype. Although gene flow was ongoing across both hybrid zones and pairwiseFSTbetween subspecies was extremely low, we found striking differences in the extent of hybridization. In the more phenotypically differentiated subspecies pair, clines for ancestry, wing length and ventral coloration were steep and coincident, suggestive of strong isolation and, potentially, selection associated with phenotype. In the less phenotypically differentiated pair, gene flow and phenotypic variation occurred over a wide geographic span, indicative of weaker isolation. Traits associated with genome‐wide differentiation in allopatry may thus also contribute to isolation in sympatry. We discuss potentially important additional roles for evolutionary history and ecology in shaping variation in the extent hybridization between closely related pairs of subspecies.

    more » « less
  4. Abstract

    A major goal of speciation research is to reveal the genomic signatures that accompany the speciation process. Genome scans are routinely used to explore genome‐wide variation and identify highly differentiated loci that may contribute to ecological divergence, but they do not incorporate spatial, phenotypic or environmental data that might enhance outlier detection. Geographic cline analysis provides a potential framework for integrating diverse forms of data in a spatially explicit framework, but has not been used to study genome‐wide patterns of divergence. Aided by a first‐draft genome assembly, we combined anFCTscan and geographic cline analysis to characterize patterns of genome‐wide divergence between divergent pollination ecotypes ofMimulus aurantiacus.FCTanalysis of 58 872SNPs generated viaRAD‐seq revealed little ecotypic differentiation (meanFCT = 0.041), although a small number of loci were moderately‐to‐highly diverged. Consistent with our previous results from the geneMaMyb2, which contributes to differences in flower colour, 130 loci have cline shapes that recapitulate the spatial pattern of trait divergence, suggesting that they may reside in or near the genomic regions that contribute to pollinator isolation. In the narrow hybrid zone between the ecotypes, extensive admixture among individuals and low linkage disequilibrium between markers indicate that most outlier loci are scattered throughout the genome, rather than being restricted to one or a few divergent regions. In addition to revealing the genomic consequences of ecological divergence in this system, we discuss how geographic cline analysis is a powerful but under‐utilized framework for studying genome‐wide patterns of divergence.

    more » « less
  5. Abstract

    Many parasites infect multiple hosts, but estimating the transmission across host species remains a key challenge in disease ecology. We investigated the within and across host species dynamics of canine distemper virus (CDV) in grizzly bears (Ursus arctos) and wolves (Canis lupus) of the Greater Yellowstone Ecosystem (GYE). We hypothesized that grizzly bears may be more likely to be exposed toCDVduring outbreaks in the wolf population because grizzly bears often displace wolves while scavenging carcasses. We used serological data collected from 1984 to 2014 in conjunction with Bayesian state‐space models to infer the temporal dynamics ofCDV. These models accounted for the unknown timing of pathogen exposure, and we assessed how different testing thresholds and the potential for testing errors affected our conclusions. We identified three mainCDVoutbreaks (1999, 2005, and 2008) in wolves, which were more obvious when we used higher diagnostic thresholds to qualify as seropositive. There was some evidence for increased exposure rates in grizzly bears in 2005, but the magnitude of the wolf effect on bear exposures was poorly estimated and depended upon our prior distributions. Grizzly bears were exposed toCDVprior to wolf reintroduction and during time periods outside of known wolf outbreaks, thus wolves are only one of several potential routes for grizzly bear exposures. Our modeling approach accounts for several of the shortcomings of serological data and is applicable to many wildlife disease systems, but is most informative when testing intervals are short.CDVcirculates in a wide range of carnivore species, but it remains unclear whether the disease persists locally within theGYEcarnivore community or is periodically reintroduced from distant regions with larger host populations.

    more » « less