Model-based genotype and ancestry estimation for potential hybrids with mixed-ploidy

Selection leads to remarkable variability in the outcomes of hybridisation across replicate hybrid zones

Hybrid zones have been viewed as an opportunity to see speciation in action. When hybrid zones are replicated, it is assumed that if the same genetic incompatibilities are maintaining reproductive isolation across all instances of secondary contact, those incompatibilities should be identifiable by consistent patterns in the genome. In contrast, changes in allele frequencies due to genetic drift should be idiosyncratic for each hybrid zone. To test this assumption, we simulated 20 replicates of each of 12 hybrid zone scenarios with varied genetic incompatibilities, rates of migration, selection and different starting population size ratios of parental species. We found remarkable variability in the outcomes of hybridisation in replicate hybrid zones, particularly with Bateson-Dobzhansky-Muller incompatibilities and strong selection. We found substantial differences among replicates in the overall genomic composition of individuals, including admixture proportions, inter-specific ancestry complement and number of ancestry junctions. Additionally, we found substantial variation in genomic clines among replicates at focal loci, regardless of locus-specific selection. We conclude that processes other than selection are responsible for some consistent outcomes of hybridisation, whereas selection on incompatibilities can lead to genomically widespread and highly variable outcomes. We highlight the challenge of mapping between pattern and process in hybrid zones and call attention to how selection against incompatibilities will commonly lead to variable outcomes. We hope that this study informs future research on replicate hybrid zones and encourages further development of statistical techniques, theoretical models and exploration of additional axes of variation to understand reproductive isolation.

Genetic and behavioral differences between above and below ground Culex pipiens bioforms

Efficiency of mosquito-borne disease transmission is dependent upon both the preference and fidelity of mosquitoes as they seek the blood of vertebrate hosts. While mosquitoes select their blood hosts through multi-modal integration of sensory cues, host-seeking is primarily an odor-guided behavior. Differences in mosquito responses to hosts and their odors have been demonstrated to have a genetic component, but the underlying genomic architecture of these responses has yet to be fully resolved. Here, we provide the first characterization of the genomic architecture of host preference in the polymorphic mosquito species, Culex pipiens. The species exists as two morphologically identical bioforms, each with distinct avian and mammalian host preferences. Cx. pipiens females with empirically measured host responses were prepared into reduced representation DNA libraries and sequenced to identify genomic regions associated with host preference. Multiple genomic regions associated with host preference were identified on all 3 Culex chromosomes, and these genomic regions contained clusters of chemosensory genes, as expected based on work in Anopheles gambiae complex mosquitoes and in Aedes aegypti. One odorant receptor and one odorant binding protein gene showed one-to-one orthologous relationships to differentially expressed genes in A. gambiae complex members with divergent host preferences. Overall, our work identifies a distinct set of odorant receptors and odorant binding proteins that may enable Cx. pipiens females to distinguish between their vertebrate blood host species, and opens avenues for future functional studies that could measure the unique contributions of each gene to host preference phenotypes.

The genetic consequences of range expansion and its influence on diploidization in polyploids

Despite newly formed polyploids being subjected to myriad fitness consequences, the relative prevalence of polyploidy both contemporarily and in ancestral branches of the tree of life suggests alternative advantages that outweigh these consequences. One proposed advantage is that polyploids may more easily colonize novel habitats such as deglaciated areas. However, previous research conducted in diploids suggests that range expansion comes with a fitness cost as deleterious mutations may fix rapidly on the expansion front. Here, we interrogate the potential consequences of expansion in polyploids by conducting spatially explicit forward-in-time simulations to investigate how ploidy and inheritance patterns impact the relative ability of polyploids to expand their range. We show that under realistic dominance models, autopolyploids suffer greater fitness reductions than diploids as a result of range expansion due to the fixation of increased mutational load that is masked in the range core. Alternatively, the disomic inheritance of allopolyploids provides a shield to this fixation resulting in minimal fitness consequences. In light of this advantage provided by disomy, we investigate how range expansion may influence cytogenetic diploidization through the reversion to disomy in autotetraploids. We show that under a wide range of parameters investigated for two models of diploidization, disomy frequently evolves more rapidly on the expansion front than in the range core, and that this dynamic inheritance model has additional effects on fitness. Together our results point to a complex interaction between dominance, ploidy, inheritance, and recombination on fitness as a population spreads across a geographic range.

Considerable genetic diversity and structure despite narrow endemism and limited ecological specialization in the Hayden's ringlet, Coenonympha haydenii

Understanding the processes that underlie the development of population genetic structure is central to the study of evolution. Patterns of genetic structure, in turn, can reveal signatures of isolation by distance (IBD), barriers to gene flow, or even the genesis of speciation. However, it is unclear how severe range restriction might impact the processes that dominate the development of genetic structure. In narrow endemic species, is population structure likely to be adaptive in nature, or rather the result of genetic drift? In this study, we investigated patterns of genetic diversity and structure in the narrow endemic Hayden's ringlet butterfly. Specifically, we asked to what degree genetic structure in the Hayden's ringlet can be explained by IBD, isolation by resistance (IBR) (in the form of geographic or ecological barriers to migration between populations), and isolation by environment (in the form of differences in host plant availability and preference). We employed a genotyping-by-sequencing (GBS) approach coupled with host preference assays, Bayesian modelling, and population genomic analyses to answer these questions. Our results suggest that despite their restricted range, levels of genetic diversity in the Hayden's ringlet are comparable to those seen in more widespread butterfly species. Hayden's ringlets showed a strong preference for feeding on grasses relative to sedges, but neither larval preference nor potential host availability at sampling sites correlated with genetic structure. We conclude that geography, in the form of IBR and simple IBD, was the major driver of contemporary patterns of differentiation in this narrow endemic species.

Botany and geogenomics: Constraining geological hypotheses in the neotropics with large-scale genetic data derived from plants

Decades of empirical research have revealed how the geological history of our planet shaped plant evolution by establishing well-known patterns (e.g., how mountain uplift resulted in high rates of diversification and replicate radiations in montane plant taxa). This follows a traditional approach where botanical data are interpreted in light of geological events. In this synthesis, I instead describe how by integrating natural history, phylogenetics, and population genetics, botanical research can be applied alongside geology and paleontology to inform our understanding of past geological and climatic processes. This conceptual shift aligns with the goals of the emerging field of geogenomics. In the neotropics, plant geogenomics is a powerful tool for the reciprocal exploration of two long standing questions in biology and geology: how the dynamic landscape of the region came to be and how it shaped the evolution of the richest flora. Current challenges that are specific to analytical approaches for plant geogenomics are discussed. I describe the scale at which various geological questions can be addressed from biological data and what makes some groups of plants excellent model systems for geogenomics research. Although plant geogenomics is discussed with reference to the neotropics, the recommendations given here for approaches to plant geogenomics can and should be expanded to exploring long-standing questions on how the earth evolved with the use of plant DNA.

Patterns of admixture and introgression in a mosaic Mytilus galloprovincialis and Mytilus edulis hybrid zone in SW England

The study of hybrid zones offers important insights into speciation. Earlier studies on hybrid populations of the marine mussel species Mytilus edulis and Mytilus galloprovincialis in SW England provided evidence of admixture but were constrained by the limited number of molecular markers available. We use 57 ancestry-informative SNPs, most of which have been mapped genetically, to provide evidence of distinctive differences between admixed populations in SW England and asymmetrical introgression from M. edulis to M. galloprovincialis. We combine the genetic study with analysis of phenotypic traits of potential ecological and adaptive significance. We demonstrate that hybrid individuals have brown mantle edges unlike the white or purple in the parental species, suggesting allelic or non-allelic genomic interactions. We report differences in gonad development stage between the species consistent with a prezygotic barrier between the species. By incorporating results from publications dating back to 1980, we confirm the long-term stability of the hybrid zone despite higher viability of M. galloprovincialis. This stability coincides with a dramatic change in temperature of UK coastal waters and suggests that these hybrid populations might be resisting the effects of global warming. However, a single SNP locus associated with the Notch transmembrane signalling protein shows a markedly different pattern of variation to the others and might be associated with adaptation of M. galloprovincialis to colder northern temperatures.

Estimating admixture pedigrees of recent hybrids without a contiguous reference genome

The genome of recently admixed individuals or hybrids has characteristic genetic patterns that can be used to learn about their recent admixture history. One of these are patterns of interancestry heterozygosity, which can be inferred from SNP data from either called genotypes or genotype likelihoods, without the need for information on genomic location. This makes them applicable to a wide range of data that are often used in evolutionary and conservation genomic studies, such as low-depth sequencing mapped to scaffolds and reduced representation sequencing. Here we implement maximum likelihood estimation of interancestry heterozygosity patterns using two complementary models. We furthermore develop apoh (Admixture Pedigrees of Hybrids), a software that uses estimates of paired ancestry proportions to detect recently admixed individuals or hybrids, and to suggest possible admixture pedigrees. It furthermore calculates several hybrid indices that make it easier to identify and rank possible admixture pedigrees that could give rise to the estimated patterns. We implemented apoh both as a command line tool and as a Graphical User Interface that allows the user to automatically and interactively explore, rank and visualize compatible recent admixture pedigrees, and calculate the different summary indices. We validate the performance of the method using admixed family trios from the 1000 Genomes Project. In addition, we show its applicability on identifying recent hybrids from RAD-seq data of Grant's gazelle (Nanger granti and Nanger petersii) and whole genome low-depth data of waterbuck (Kobus ellipsiprymnus) which shows complex admixture of up to four populations.

Variable rates of hybridization among contact zones between a pair of topminnow species, Fundulus notatus and F. olivaceus

Pairs of species that exhibit broadly overlapping distributions, and multiple geographically isolated contact zones, provide opportunities to investigate the mechanisms of reproductive isolation. Such naturally replicated systems have demonstrated that hybridization rates can vary substantially among populations, raising important questions about the genetic basis of reproductive isolation. The topminnows, Fundulus notatus and F. olivaceus, are reciprocally monophyletic, and co-occur in drainages throughout much of the central and southern United States. Hybridization rates vary substantially among populations in isolated drainage systems. We employed genome-wide sampling to investigate geographic variation in hybridization, and to assess the possible importance of chromosome fusions to reproductive isolation among nine separate contact zones. The species differ by chromosomal rearrangements resulting from Robertsonian (Rb) fusions, so we hypothesized that Rb fusion chromosomes would serve as reproductive barriers, exhibiting steeper genomic clines than the rest of the genome. We observed variation in hybridization dynamics among drainages that ranged from nearly random mating to complete absence of hybridization. Contrary to predictions, our use of genomic cline analyses on mapped species-diagnostic SNP markers did not indicate consistent patterns of variable introgression across linkage groups, or an association between Rb fusions and genomic clines that would be indicative of reproductive isolation. We did observe a relationship between hybridization rates and population phylogeography, with the lowest rates of hybridization tending to be found in populations inferred to have had the longest histories of drainage sympatry. Our results, combined with previous studies of contact zones between the species, support population history as an important factor in explaining variation in hybridization rates.

Population genomic evidence of selection on structural variants in a natural hybrid zone

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 two Lycaeides butterfly taxa (L. melissa and Jackson Hole Lycaeides) 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 Hole Lycaeides ancestry 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 Hole Lycaeides. Excess Jackson Hole Lycaeides ancestry 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.

Population Genomic Analysis of Diploid-Autopolyploid Species

This chapter outlines an empirical analysis of genome-wide single-nucleotide polymorphism (SNP) variation and its underlying drivers among multiple natural populations within a diploid-autopolyploid species. The aim is to reconstruct the genetic structure among natural populations of varying ploidy and infer footprints of selection in these populations, framed around specific questions that are typically encountered when analyzing a mixed-ploidy data set,e.g., addressing the relevance of natural whole-genome duplication for speciation and adaptation. We briefly review the options for the analysis of polyploid population genomic data involving variant calling, population structure, demographic history inference, and selection scanning approaches. Further, we provide suggestions for methods and associated software, possible caveats, and examples of their application to mixed-ploidy and autopolyploid data sets.

Limited evidence of a genetic basis for sex determination in the common creek chub, Semotilus atromaculatus

Sexual reproduction is almost universal in vertebrates; therefore, each animal species which uses it must have a mechanism for designating sex as male or female. Fish, especially, have a wide range of sex determining systems. In the present study, we aimed to identify a genetic basis for sex determination in the common creek chub (Semotilus atromaculatus) using genotyping-by-sequencing data. No sex-associated markers were found by RADSex or a GWAS using GEMMA; however, Weir and Cockerham locus-specific F_ST analysis and discriminant analysis of principal components revealed genetic differentiation between the sexes at several loci. While no explicit sex determination mechanism has been yet discovered in creek chub, these loci are potential candidates for future studies. Incompatible systems are thought to increase reproductive isolation but interspecific hybridization is common among groups such as cyprinid minnows; thus, studies such as ours can provide insight into hybridization and evolutionary diversification of this clade. We also highlight technical challenges involved in studying sex determination in evolutionary groups with extremely variable mechanisms and without heteromorphic sex chromosomes.

Ecological outcomes of hybridization vary extensively in Catostomus fishes

Hybridization outcomes vary geographically and can depend on the environment. Hybridization can also reshape biotic interactions, leading to ecological shifts. If hybrids function differently ecologically in ways that enhance or reduce fitness, and those ecological roles vary geographically, ecological factors might explain variation in hybridization outcomes. However, relatively few studies have focused on ecological traits of hybrids. We compared the feeding ecology of Catostomus fish species and hybrids by using stable isotopes (δ13 C and δ15 N) as a proxy for diet and habitat use, and compared two native species, an introduced species, and three interspecific hybrid crosses. We included hybrids and parental species from seven rivers where hybridization outcomes vary. Relative isotopic niches of native species varied geographically, but native species did not fully overlap in isotopic space in any river sampled, suggesting little overlap of resource use between historically sympatric species. The introduced species overlapped with one or both native species in every river, suggesting similar resource use and potential competition. Hybrids occupied intermediate, matching, or more transgressive isotopic niches, and varied within and among rivers. Ecological outcomes of hybridization varied across locations, implying that hybridization might have unpredictable, idiosyncratic ecological effects.

Additive genetic effects in interacting species jointly determine the outcome of caterpillar herbivory

Plant-insect interactions are common and important in basic and applied biology. Trait and genetic variation can affect the outcome and evolution of these interactions, but the relative contributions of plant and insect genetic variation and how these interact remain unclear and are rarely subject to assessment in the same experimental context. Here, we address this knowledge gap using a recent host-range expansion onto alfalfa by the Melissa blue butterfly. Common garden rearing experiments and genomic data show that caterpillar performance depends on plant and insect genetic variation, with insect genetics contributing to performance earlier in development and plant genetics later. Our models of performance based on caterpillar genetics retained predictive power when applied to a second common garden. Much of the plant genetic effect could be explained by heritable variation in plant phytochemicals, especially saponins, peptides, and phosphatidyl cholines, providing a possible mechanistic understanding of variation in the species interaction. We find evidence of polygenic, mostly additive effects within and between species, with consistent effects of plant genotype on growth and development across multiple butterfly species. Our results inform theories of plant-insect coevolution and the evolution of diet breadth in herbivorous insects and other host-specific parasites.

Variable hybridization between two Lake Tanganyikan cichlid species in recent secondary contact

Closely related taxa frequently exist in sympatry before the evolution of robust reproductive barriers, which can lead to substantial gene flow. Post-divergence gene flow can promote several disparate trajectories of divergence ranging from the erosion of distinctiveness and eventual collapse of the taxa to the strengthening of reproductive isolation. Among many relevant factors, understanding the demographic history of divergence (e.g. divergence time, extent of historical gene flow) can be particularly informative when examining contemporary gene flow between closely related taxa because this history can influence gene flow's prevalence and consequences. Here, we used genotyping-by-sequencing data to investigate speciation and contemporary hybridization in two closely related and sympatrically distributed Lake Tanganyikan cichlid species in the genus Petrochromis. Demographic modeling supported a speciation scenario involving divergence in isolation followed by secondary contact with bidirectional gene flow. Further investigation of this recent gene flow found evidence of ongoing hybridization between the species that varied in extent between different co-occurring populations. Relationships between abundance and the degree of admixture across populations suggest that the availability of conspecific mates may influence patterns of hybridization. These results, together with the observation that sets of recently diverged cichlid taxa are generally geographically separated in the lake, suggest that ongoing speciation in Lake Tanganyikan cichlids relies on initial spatial isolation. Additionally, the spatially heterogeneous patterns of admixture between the Petrochromis species illustrates the complexities of hybridization when species are in recent secondary contact.

Hybridization decreases native cutthroat trout reproductive fitness

Examining natural selection in wild populations is challenging, but crucial to understanding many ecological and evolutionary processes. Additionally, in hybridizing populations, natural selection may be an important determinant of the eventual outcome of hybridization. We characterized several components of relative fitness in hybridizing populations of Yellowstone cutthroat trout and rainbow trout in an effort to better understand the prolonged persistence of both parental species despite predictions of extirpation. Thousands of genomic loci enabled precise quantification of hybrid status in adult and subsequent juvenile generations; a subset of those data also identified parent-offspring relationships. We used linear models and simulations to assess the effects of ancestry on reproductive output and mate choice decisions. We found a relatively low number of late-stage (F3+) hybrids and an excess of F2 juveniles relative to the adult generation in one location, which suggests the presence of hybrid breakdown decreasing the fitness of F2+ hybrids later in life. Assessments of reproductive output showed that Yellowstone cutthroat trout are more likely to successfully reproduce and produce slightly more offspring than their rainbow trout and hybrid counterparts. Mate choice appeared to be largely random, though we did find statistical support for slight female preference for males of similar ancestry. Together, these results show that native Yellowstone cutthroat trout are able to outperform rainbow trout in terms of reproduction and suggests that management action to exclude rainbow trout from spawning locations may bolster the now-rare Yellowstone cutthroat trout.

The Genetic Population Structure of Lake Tanganyika's Lates Species Flock, an Endemic Radiation of Pelagic Top Predators

Understanding genetic connectivity plays a crucial role in species conservation decisions, and genetic connectivity is an important component of modern fisheries management. In this study, we investigated the population genetics of four endemic Lates species of Lake Tanganyika (Lates stappersii, L. microlepis, L. mariae, and L. angustifrons) using reduced-representation genomic sequencing methods. We find the four species to be strongly differentiated from one another (mean interspecific FST = 0.665), with no evidence for contemporary admixture. We also find evidence for strong genetic structure within L. mariae, with the majority of individuals from the most southern sampling site forming a genetic group that is distinct from the individuals at other sampling sites. We find evidence for much weaker structure within the other three species (L. stappersii, L. microlepis, and L. angustifrons). Our ability to detect this weak structure despite small and unbalanced sample sizes and imprecise geographic sampling locations suggests the possibility for further structure undetected in our study. We call for further research into the origins of the genetic differentiation in these four species-particularly that of L. mariae-which may be important for conservation and management of this culturally and economically important clade of fishes.

When do autopolyploids need poly-sequencing data?

The sequencing depth required to genotype autopolyploid populations is a very controversial topic. Different studies have adopted variable depth values without a clear guide on the optimal sequencing depth value. Many studies suggest high depth thresholds for different ploidies that may not be practical and substantially increase the overall genotyping cost for different projects. However, such conservative thresholds may not be required to achieve the most common research goals. In fact, some recent reports in the field of quantitative genetics found that much lower sequencing depth thresholds could achieve the same accuracy as high depth thresholds. In this manuscript, I discuss when researchers need to use stringent sequencing depth thresholds and when they can use more relaxed ones. I support my argument by calculating the probabilities of sampling different homologues at a given sequencing depth. I also discuss the uses and the uncertainty in calculating a continuous allelic dosage as the proportion of sequencing reads that hold the alternative allele, which is becoming a common method now in quantitative genetics to replace discrete dosage estimation.

Genomic and common garden approaches yield complementary results for quantifying environmental drivers of local adaptation in rubber rabbitbrush, a foundational Great Basin shrub

The spatial structure of genomic and phenotypic variation across populations reflects historical and demographic processes as well as evolution via natural selection. Characterizing such variation can provide an important perspective for understanding the evolutionary consequences of changing climate and for guiding ecological restoration. While evidence for local adaptation has been traditionally evaluated using phenotypic data, modern methods for generating and analyzing landscape genomic data can directly quantify local adaptation by associating allelic variation with environmental variation. Here, we analyze both genomic and phenotypic variation of rubber rabbitbrush (Ericameria nauseosa), a foundational shrub species of western North America. To quantify landscape genomic structure and provide perspective on patterns of local adaptation, we generated reduced representation sequencing data for 17 wild populations (222 individuals; 38,615 loci) spanning a range of environmental conditions. Population genetic analyses illustrated pronounced landscape genomic structure jointly shaped by geography and environment. Genetic-environment association (GEA) analyses using both redundancy analysis (RDA) and a machine-learning approach (Gradient Forest) indicated environmental variables (precipitation seasonality, slope, aspect, elevation, and annual precipitation) influenced spatial genomic structure and were correlated with allele frequency shifts indicative of local adaptation at a consistent set of genomic regions. We compared our GEA-based inference of local adaptation with phenotypic data collected by growing seeds from each population in a greenhouse common garden. Population differentiation in seed weight, emergence, and seedling traits was associated with environmental variables (e.g., precipitation seasonality) that were also implicated in GEA analyses, suggesting complementary conclusions about the drivers of local adaptation across different methods and data sources. Our results provide a baseline understanding of spatial genomic structure for E. nauseosa across the western Great Basin and illustrate the utility of GEA analyses for detecting the environmental causes and genetic signatures of local adaptation in a widely distributed plant species of restoration significance.

Hierarchical genetic structure and implications for conservation of the world's largest salmonid, Hucho taimen

Population genetic analyses can evaluate how evolutionary processes shape diversity and inform conservation and management of imperiled species. Taimen (Hucho taimen), the world’s largest freshwater salmonid, is threatened, endangered, or extirpated across much of its range due to anthropogenic activity including overfishing and habitat degradation. We generated genetic data using high throughput sequencing of reduced representation libraries for taimen from multiple drainages in Mongolia and Russia. Nucleotide diversity estimates were within the range documented in other salmonids, suggesting moderate diversity despite widespread population declines. Similar to other recent studies, our analyses revealed pronounced differentiation among the Arctic (Selenge) and Pacific (Amur and Tugur) drainages, suggesting historical isolation among these systems. However, we found evidence for finer-scale structure within the Pacific drainages, including unexpected differentiation between tributaries and the mainstem of the Tugur River. Differentiation across the Amur and Tugur basins together with coalescent-based demographic modeling suggests the ancestors of Tugur tributary taimen likely diverged in the eastern Amur basin, prior to eventual colonization of the Tugur basin. Our results suggest the potential for differentiation of taimen at different geographic scales, and suggest more thorough geographic and genomic sampling may be needed to inform conservation and management of this iconic salmonid.

Genomic time-series data show that gene flow maintains high genetic diversity despite substantial genetic drift in a butterfly species

Effective population size affects the efficacy of selection, rate of evolution by drift, and neutral diversity levels. When species are subdivided into multiple populations connected by gene flow, evolutionary processes can depend on global or local effective population sizes. Theory predicts that high levels of diversity might be maintained by gene flow, even very low levels of gene flow, consistent with species long-term effective population size, but tests of this idea are mostly lacking. Here, we show that Lycaeides buttery populations maintain low contemporary (variance) effective population sizes (e.g., ~200 individuals) and thus evolve rapidly by genetic drift. In contrast, populations harbored high levels of genetic diversity consistent with an effective population size several orders of magnitude larger. We hypothesized that the differences in the magnitude and variability of contemporary versus long-term effective population sizes were caused by gene flow of sufficient magnitude to maintain diversity but only subtly affect evolution on generational time scales. Consistent with this hypothesis, we detected low but non-trivial gene flow among populations. Furthermore, using short-term population-genomic time-series data, we documented patterns consistent with predictions from this hypothesis, including a weak but detectable excess of evolutionary change in the direction of the mean (migrant gene pool) allele frequencies across populations, and consistency in the direction of allele frequency change over time. The documented decoupling of diversity levels and short-term change by drift in Lycaeides has implications for our understanding of contemporary evolution and the maintenance of genetic variation in the wild.

Multigenerational backcrossing and introgression between two woodrat species at an abrupt ecological transition

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 F₁ individuals 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.

Tracing evolutionary history and admixture in mixed-ploidy systems

Polyploidy, resulting from whole genome duplication, is a widespread phenomenon throughout Eukaryotic kingdoms. It is estimated that 15% of speciation events in plants is due to polyploidization and 16% of plant species encompass ploidy variation. In spite of the evolutionary and economic significance of polyploidy, there is a limited set of tools that would allow routine population genetic and genomic analysis of polyploid systems, in particular for the inference of population diversity and differentiation from large genome-wide data sets (Dufresne et al., 2014; Meirmans, 2020). Such a shortage is striking especially when compared to the rapid development of such tools in diploid systems over the last decade. Consequently, population genomic research in polyploids is still lagging behind diploids, especially in autopolyploids, for example, polyploids with multiple (>2) similar genome copies. In this issue of Molecular Ecology Resources, Shastry et al. (2021) develops the first software specifically designed to infer population structure and ancestry in large genome-wide mixed-ploidy data sets, opening novel opportunities for polyploid population genomic analysis.