An official website of the United States government
Abstract Structural variants (SVs) can promote speciation by directly causing reproductive isolation or by suppressing recombination across large genomic regions. Whereas examples of each mechanism have been documented, systematic tests of the role of SVs in speciation are lacking. Here, we take advantage of long‐read (Oxford nanopore) whole‐genome sequencing and a hybrid zone between twoLycaeidesbutterfly taxa (L.melissaand Jackson HoleLycaeides) to comprehensively evaluate genome‐wide patterns of introgression for SVs and relate these patterns to hypotheses about speciation. We found >100,000 SVs segregating within or between the two hybridizing species. SVs and SNPs exhibited similar levels of genetic differentiation between species, with the exception of inversions, which were more differentiated. We detected credible variation in patterns of introgression among SV loci in the hybrid zone, with 562 of 1419 ancestry‐informative SVs exhibiting genomic clines that deviated from null expectations based on genome‐average ancestry. Overall, hybrids exhibited a directional shift towards Jackson HoleLycaeidesancestry at SV loci, consistent with the hypothesis that these loci experienced more selection on average than SNP loci. Surprisingly, we found that deletions, rather than inversions, showed the highest skew towards excess ancestry from Jackson HoleLycaeides. Excess Jackson HoleLycaeidesancestry in hybrids was also especially pronounced for Z‐linked SVs and inversions containing many genes. In conclusion, our results show that SVs are ubiquitous and suggest that SVs in general, but especially deletions, might disproportionately affect hybrid fitness and thus contribute to reproductive isolation.
more »
« less
This content will become publicly available on November 1, 2025
A Next Generation of Hierarchical Bayesian Analyses of Hybrid Zones Enables Model‐Based Quantification of Variation in Introgression in R
ABSTRACT Hybrid zones, where genetically distinct groups of organisms meet and interbreed, offer valuable insights into the nature of species and speciation. Here, we present a new R package,bgchm, for population genomic analyses of hybrid zones. This R package extends and updates the existingbgcsoftware and combines Bayesian analyses of hierarchical genomic clines with Bayesian methods for estimating hybrid indexes, interpopulation ancestry proportions, and geographic clines. Compared to existing software,bgchmoffers enhanced efficiency through Hamiltonian Monte Carlo sampling and the ability to work with genotype likelihoods combined with a hierarchical Bayesian approach, enabling inference for diverse types of genetic data sets. The package also facilitates the quantification of introgression patterns across genomes, which is crucial for understanding reproductive isolation and speciation genetics. We first describe the models underlyingbgchmand then provide an overview of the R package and illustrate its use through the analysis of simulated and empirical data sets. We show thatbgchmgenerates accurate estimates of model parameters under a variety of conditions, especially when the genetic loci analyzed are highly ancestry informative. This includes relatively robust estimates of genome‐wide variability in clines, which has not been the focus of previous models and methods. We also illustrate how both selection and genetic drift contribute to variability in introgression among loci and how additional information can be used to help distinguish these contributions. We conclude by describing the promises and limitations ofbgchm, comparingbgchmto other software for genomic cline analyses, and identifying areas for fruitful future development.
more »
« less
- Award ID(s):
- 1844941
- PAR ID:
- 10584198
- Publisher / Repository:
- John Wiley & Sons Ltd.
- Date Published:
- Journal Name:
- Ecology and Evolution
- Volume:
- 14
- Issue:
- 11
- ISSN:
- 2045-7758
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
-
Abstract When organisms experience secondary contact after allopatric divergence, genomic regions can introgress differentially depending on their relationships with adaptation, reproductive isolation, recombination, and drift. Analyses of genome‐wide patterns of divergence and introgression could provide insight into the outcomes of hybridization and the potential relationship between allopatric divergence and reproductive isolation. Here, we generate population genetic data (26,262 SNPs; 353 individuals) using a reduced‐representation sequencing approach to quantify patterns of ancestry, differentiation, and introgression between a pair of ecologically distinct mammals—the desert woodrat (N.lepida) and Bryant's woodrat (N.bryanti)—that hybridize at a sharp ecotone in southern California. Individual ancestry estimates confirmed that hybrids were rare in this bimodal hybrid zone, and entirely consisted of a few F1individuals and a broad range of multigenerational backcrosses. Genomic cline analyses indicated more than half of loci had elevated introgression from one genomic background into the other. However, introgression was not associated with relative or absolute measures of divergence, and loci with extreme values for both were not typically found near detoxification enzymes previously implicated in dietary specialization for woodrats. The decoupling of differentiation and introgression suggests that processes other than adaptation, such as drift, may underlie the extreme clines at this contact zone.more » « less
-
Abstract BackgroundPatterns of multi-locus differentiation (i.e., genomic clines) often extend broadly across hybrid zones and their quantification can help diagnose how species boundaries are shaped by adaptive processes, both intrinsic and extrinsic. In this sense, the transitioning of loci across admixed individuals can be contrasted as a function of the genome-wide trend, in turn allowing an expansion of clinal theory across a much wider array of biodiversity. However, computational tools that serve to interpret and consequently visualize ‘genomic clines’ are limited, and users must often write custom, relatively complex code to do so. ResultsHere, we introduce the ClineHelpR R-package for visualizing genomic clines and detecting outlier loci using output generated by two popular software packages, bgc and Introgress. ClineHelpR bundles both input generation (i.e., filtering datasets and creating specialized file formats) and output processing (e.g., MCMC thinning and burn-in) with functions that directly facilitate interpretation and hypothesis testing. Tools are also provided for post-hoc analyses that interface with external packages such as ENMeval and RIdeogram. ConclusionsOur package increases the reproducibility and accessibility of genomic cline methods, thus allowing an expanded user base and promoting these methods as mechanisms to address diverse evolutionary questions in both model and non-model organisms. Furthermore, the ClineHelpR extended functionality can evaluate genomic clines in the context of spatial and environmental features, allowing users to explore underlying processes potentially contributing to the observed patterns and helping facilitate effective conservation management strategies.more » « less
-
Abstract Invasive species can impact native populations through competition, predation, habitat alteration, and disease transmission, but also genetically through hybridization. Potential outcomes of hybridization span the continuum from extinction to hybrid speciation and can be further complicated by anthropogenic habitat disturbance. Hybridization between the native green anole lizard (Anolis carolinensis) and a morphologically similar invader (A. porcatus) in south Florida provides an ideal opportunity to study interspecific admixture across a heterogeneous landscape. We used reduced‐representation sequencing to describe introgression in this hybrid system and to test for a relationship between urbanization and non‐native ancestry. Our findings indicate that hybridization between green anole lineages was probably a limited, historic event, producing a hybrid population characterized by a diverse continuum of ancestry proportions. Genomic cline analyses revealed rapid introgression and disproportionate representation of non‐native alleles at many loci and no evidence for reproductive isolation between parental species. Three loci were associated with urban habitat characteristics; urbanization and non‐native ancestry were positively correlated, although this relationship did not remain significant when accounting for spatial nonindependence. Ultimately, our study demonstrates the persistence of non‐native genetic material even in the absence of ongoing immigration, indicating that selection favouring non‐native alleles can override the demographic limitation of low propagule pressure. We also note that not all outcomes of admixture between native and non‐native species should be considered intrinsically negative. Hybridization with ecologically robust invaders can lead to adaptive introgression, which may facilitate the long‐term survival of native populations otherwise unable to adapt to anthropogenically mediated global change.more » « less
-
Abstract The outcomes of speciation across organismal dimensions (e.g., ecological, genetic, phenotypic) are often assessed using phylogeographic methods. At one extreme, reproductively isolated lineages represent easily delimitable species differing in many or all dimensions, and at the other, geographically distinct genetic segments introgress across broad environmental gradients with limited phenotypic disparity. In the ambiguous gray zone of speciation, where lineages are genetically delimitable but still interacting ecologically, it is expected that these lineages represent species in the context of ontology and the evolutionary species concept when they are maintained over time with geographically well‐defined hybrid zones, particularly at the intersection of distinct environments. As a result, genetic structure is correlated with environmental differences and not space alone, and a subset of genes fail to introgress across these zones as underlying genomic differences accumulate. We present a set of tests that synthesize species delimitation with the speciation process. We can thereby assess historical demographics and diversification processes while understanding how lineages are maintained through space and time by exploring spatial and genome clines, genotype‐environment interactions, and genome scans for selected loci. Employing these tests in eight lineage‐pairs of snakes in North America, we show that six pairs represent 12 “good” species and that two pairs represent local adaptation and regional population structure. The distinct species pairs all have the signature of divergence before or near the mid‐Pleistocene, often with low migration, stable hybrid zones of varying size, and a subset of loci showing selection on alleles at the hybrid zone corresponding to transitions between distinct ecoregions. Locally adapted populations are younger, exhibit higher migration, and less ecological differentiation. Our results demonstrate that interacting lineages can be delimited using phylogeographic and population genetic methods that properly integrate spatial, temporal, and environmental data.more » « less
