Highly variable reproductive isolation among pairs of Catostomus species

Phylogeography and paleoclimatic range dynamics explain variable outcomes to contact across a species' range

Replicability of divergence after contact is a poorly characterized process, particularly in the contexts of phylogeography and postglacial range dynamics within species. Using contact zones located at the leading-, mid- and rear-edges of a species' range, we examined variation in outcomes to contact between divergent lineages of Campanula americana. We investigated whether contact zones vary in quantity and directionality of gene flow, how phylogeographic structure differs between contact zones, and how historic range dynamics may affect outcomes to contact. We found that all contact zones formed at similar times via primary contact yet detected significant admixture in only the rear-edge (RE) contact zone. In the northern leading-edge contact zone and the mid-range Virginia contact zone, gene flow was minimal and asymmetric. In the southern RE contact zone, gene flow was strong and symmetric. Asymmetric admixture in the leading-edge and Virginia contact zones matches the directionality of a known cosmopolitan cytonuclear incompatibility between lineages of C. americana. Our results emphasize the dependence of speciation processes on phylogeographic structure, evolutionary history and range dynamics.

The extent of introgression between incipient Clarkia species is determined by temporal environmental variation and mating system

Introgression is pervasive across the tree of life but varies across taxa, geography, and genomic regions. However, the factors modulating this variation and how they may be affected by global change are not well understood. Here, we used 200 genomes and a 15-y site-specific environmental dataset to investigate the effects of environmental variation and mating system divergence on the magnitude of introgression between a recently diverged outcrosser-selfer pair of annual plants in the genus Clarkia. These sister taxa diverged very recently and subsequently came into secondary sympatry where they form replicated contact zones. Consistent with observations of other outcrosser-selfer pairs, we found that introgression was asymmetric between taxa, with substantially more introgression from the selfer to the outcrosser. This asymmetry was caused by a bias in the direction of initial F1 hybrid formation and subsequent backcrossing. We also found extensive variation in the outcrosser's admixture proportion among contact zones, which was predicted nearly entirely by interannual variance in spring precipitation. Greater fluctuations in spring precipitation resulted in higher admixture proportions, likely mediated by the effects of spring precipitation on the expression of traits that determine premating reproductive isolation. Climate-driven hybridization dynamics may be particularly affected by global change, potentially reshaping species boundaries and adaptation to novel environments.

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.

Two transects reveal remarkable variation in gene flow on opposite ends of a European toad hybrid zone

Speciation entails a reduction in gene flow between lineages. The rates at which genomic regions become isolated varies across space and time. Barrier markers are linked to putative genes involved in (processes of) reproductive isolation, and, when observed over two transects, indicate species-wide processes. In contrast, transect-specific putative barrier markers suggest local processes. We studied two widely separated transects along the 900 km hybrid zone between Bufo bufo and B. spinosus, in northern and southern France, for ~1200 RADseq markers. We used genomic and geographic cline analyses to identify barrier markers based on their restricted introgression, and found that some markers are transect-specific, while others are shared between transects. Twenty-six barrier markers were shared across both transects, of which some are clustered in the same chromosomal region, suggesting that their associated genes are involved in reduced gene flow across the entire hybrid zone. Transect-specific barrier markers were twice as numerous in the southern than in the northern transect, suggesting that the overall barrier effect is weaker in northern France. We hypothesize that this is consistent with a longer period of secondary contact in southern France. The smaller number of introgressed genes in the northern transect shows considerably more gene flow towards the southern (B. spinosus) than the northern species (B. bufo). We hypothesize that hybrid zone movement in northern France and hybrid zone stability in southern France explain this pattern. The Bufo hybrid zone provides an excellent opportunity to separate a general barrier effect from localized gene flow-reducing conditions.

The population genetics of speciation by cascade reinforcement

Species interactions drive diverse evolutionary outcomes. Speciation by cascade reinforcement represents one example of how species interactions can contribute to the proliferation of species. This process occurs when the divergence of mating traits in response to selection against interspecific hybridization incidentally leads to reproductive isolation among populations of the same species. Here, we investigated the population genetic outcomes of cascade reinforcement in North American chorus frogs (Hylidae: Pseudacris). Specifically, we estimated the frequency of hybridization among three taxa, assessed genetic structure within the focal species, P. feriarum, and ascertained the directionality of gene flow within P. feriarum across replicated contact zones via coalescent modeling. Through field observations and preliminary experimental crosses, we assessed whether hybridization is possible under natural and laboratory conditions. We found that hybridization occurs among P. feriarum and two conspecifics at a low rate in multiple contact zones, and that gene flow within the former species is unidirectional from allopatry into sympatry with these other species in three of four contact zones studied. We found evidence of substantial genetic structuring within P. feriarum including a divergent western allopatric cluster, a behaviorally-distinct sympatric South Carolina cluster, and several genetically-overlapping clusters from the remainder of the distribution. Furthermore, we found sub-structuring between reinforced and nonreinforced populations in the two most intensely-sampled contact zones. Our literature review indicated that P. feriarum hybridizes with at least five heterospecifics at the periphery of its range providing a mechanism for further intraspecific diversification. This work strengthens the evidence for cascade reinforcement in this clade, revealing the geographic and genetic landscape upon which this process can contribute to the proliferation of species.

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.

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.

Reinforcement and the Proliferation of Species

Adaptive radiations are characterized by the rapid proliferation of species. Explaining how adaptive radiations occur therefore depends, in part, on identifying how populations become reproductively isolated-and ultimately become different species. Such reproductive isolation could arise when populations adapting to novel niches experience selection to avoid interbreeding and, consequently, evolve mating traits that minimize such hybridization via the process of reinforcement. Here, we highlight that a downstream consequence of reinforcement is divergence of conspecific populations, and this further divergence can instigate species proliferation. Moreover, we evaluate when reinforcement will-and will not-promote species proliferation. Finally, we discuss empirical approaches to test what role, if any, reinforcement plays in species proliferation and, consequently, in adaptive radiation. To date, reinforcement's downstream effects on species proliferation remain largely unknown and speculative. Because the ecological and evolutionary contexts in which adaptive radiations occur are conducive to reinforcement and its downstream consequences, adaptive radiations provide an ideal framework in which to evaluate reinforcement's role in diversification.

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

Non-random mating among individuals can lead to spatial clustering of genetically similar individuals and population stratification. This deviation from panmixia is commonly observed in natural populations. Consequently, individuals can have parentage in single populations or involving hybridization between differentiated populations. Accounting for this mixture and structure is important when mapping the genetics of traits and learning about the formative evolutionary processes that shape genetic variation among individuals and populations. Stratified genetic relatedness among individuals is commonly quantified using estimates of ancestry that are derived from a statistical model. Development of these models for polyploid and mixed-ploidy individuals and populations has lagged behind those for diploids. Here, we extend and test a hierarchical Bayesian model, called entropy, which can use low-depth sequence data to estimate genotype and ancestry parameters in autopolyploid and mixed-ploidy individuals (including sex chromosomes and autosomes within individuals). Our analysis of simulated data illustrated the trade-off between sequencing depth and genome coverage and found lower error associated with low-depth sequencing across a larger fraction of the genome than with high-depth sequencing across a smaller fraction of the genome. The model has high accuracy and sensitivity as verified with simulated data and through analysis of admixture among populations of diploid and tetraploid Arabidopsis arenosa.

Reticulate evolution as a management challenge: Patterns of admixture with phylogenetic distance in endemic fishes of western North America

Admixture in natural populations is a long-standing management challenge, with population genomic approaches offering means for adjudication. We now more clearly understand the permeability of species boundaries and the potential of admixture for promoting adaptive evolution. These issues particularly resonate in western North America, where tectonism and aridity have fragmented and reshuffled rivers over millennia, in turn promoting reticulation among endemic fishes, a situation compounded by anthropogenic habitat modifications and non-native introductions. The melding of historic and contemporary admixture has both confused and stymied management. We underscore this situation with a case study that quantifies basin-wide admixture among a group of native and introduced fishes by employing double-digest restriction site-associated DNA (ddRAD) sequencing. Our approach: (a) quantifies the admixed history of 343 suckers (10 species of Catostomidae) across the Colorado River Basin; (b) gauges admixture within the context of phylogenetic distance and "ecological specialization"; and (c) extrapolates potential drivers of introgression across hybrid crosses that involve endemic as well as invasive species. Our study extends across an entire freshwater basin and expands previous studies more limited in scope both geographically and taxonomically. Our results detected admixture involving all 10 species, with habitat alterations not only accelerating the breakdown of reproductive isolation, but also promoting introgression. Hybridization occurred across the genus despite phylogenetic distance, whereas introgression was only detected within subgenera, implicating phylogenetic distance and/or ecological specialization as drivers of reproductive isolation. Understanding the extent of admixture and reproductive isolation across multiple species serves to disentangle their reticulate evolutionary histories and provides a broadscale perspective for basin-wide conservation and management.

Gene flow and species delimitation in fishes of Western North America: Flannelmouth (Catostomus latipinnis) and Bluehead sucker (C. Pantosteus discobolus)

The delimitation of species boundaries, particularly those obscured by reticulation, is a critical step in contemporary biodiversity assessment. It is especially relevant for conservation and management of indigenous fishes in western North America, represented herein by two species with dissimilar life histories codistributed in the highly modified Colorado River (i.e., flannelmouth sucker, Catostomus latipinnis; bluehead sucker, C. (Pantosteus) discobolus). To quantify phylogenomic patterns and examine proposed taxonomic revisions, we first employed double-digest restriction site-associated DNA sequencing (ddRAD), yielding 39,755 unlinked SNPs across 139 samples. These were subsequently evaluated with multiple analytical approaches and by contrasting life history data. Three phylogenetic methods and a Bayesian assignment test highlighted similar phylogenomic patterns in each, but with considerable difference in presumed times of divergence. Three lineages were detected in bluehead sucker, supporting elevation of C. (P.) virescens to species status and recognizing C. (P.) discobolus yarrowi (Zuni bluehead sucker) as a discrete entity. Admixture in the latter necessitated a reevaluation of its contemporary and historic distributions, underscoring how biodiversity identification can be confounded by complex evolutionary histories. In addition, we defined three separate flannelmouth sucker lineages as ESUs (evolutionarily significant units), given limited phenotypic and genetic differentiation, contemporary isolation, and lack of concordance (per the genealogical concordance component of the phylogenetic species concept). Introgression was diagnosed in both species, with the Little Colorado and Virgin rivers in particular. Our diagnostic methods, and the agreement of our SNPs with previous morphological, enzymatic, and mitochondrial work, allowed us to partition complex evolutionary histories into requisite components, such as isolation versus secondary contact.

Coevolution of male and female mate choice can destabilize reproductive isolation

Sexual interactions play an important role in the evolution of reproductive isolation, with important consequences for speciation. Theoretical studies have focused on the evolution of mate preferences in each sex separately. However, mounting empirical evidence suggests that premating isolation often involves mutual mate choice. Here, using a population genetic model, we investigate how female and male mate choice coevolve under a phenotype matching rule and how this affects reproductive isolation. We show that the evolution of female preferences increases the mating success of males with reciprocal preferences, favouring mutual mate choice. However, the evolution of male preferences weakens indirect selection on female preferences and, with weak genetic drift, the coevolution of female and male mate choice leads to periodic episodes of random mating with increased hybridization (deterministic 'preference cycling' triggered by stochasticity). Thus, counterintuitively, the process of establishing premating isolation proves rather fragile if both male and female mate choice contribute to assortative mating.

Variable hybridization outcomes in trout are predicted by historical fish stocking and environmental context

Hybridization can profoundly affect the genomic composition and phenotypes of closely related species, and provides an opportunity to identify mechanisms that maintain reproductive isolation between species. Recent evidence suggests that hybridization outcomes within a species pair can vary across locations. However, we still do not know how variable outcomes of hybridization are across geographic replicates, and what mechanisms drive that variation. In this study, we described hybridization outcomes across 27 locations in the North Fork Shoshone River basin (Wyoming, USA) where native Yellowstone cutthroat trout and introduced rainbow trout co-occur. We used genomic data and hierarchical Bayesian models to precisely identify ancestry of hybrid individuals. Hybridization outcomes varied across locations. In some locations, only rainbow trout and advanced backcrossed hybrids towards rainbow trout were present, while trout in other locations had a broader range of ancestry, including both parental species and first-generation hybrids. Later-generation intermediate hybrids were rare relative to backcrossed hybrids and rainbow trout individuals. Using an individual-based simulation, we found that outcomes of hybridization in the North Fork Shoshone River basin deviate substantially from what we would expect under null expectations of random mating and no selection against hybrids. Since this deviation implies that some mechanisms of reproductive isolation function to maintain parental taxa and a diversity of hybrid types, we then modelled hybridization outcomes as a function of environmental variables and stocking history that are likely to affect prezygotic barriers to hybridization. Variables associated with history of fish stocking were the strongest predictors of hybridization outcomes, followed by environmental variables that might affect overlap in spawning time and location.

Asymmetric introgression between fishes in the Red River basin of Texas is associated with variation in water quality

When ecologically divergent taxa encounter one another, hybrid zones can form when reproductive isolation is incomplete. The location of such hybrid zones can be influenced by environmental variables, and an ecological context can provide unique insights into the mechanisms by which species diverge and are maintained. Two ecologically differentiated species of small benthic fishes, the endemic and imperiled prairie chub, Macrhybopsis australis, and the shoal chub, Macrhybopsis hyostoma, are locally sympatric within the upper Red River Basin of Texas. We integrated population genomic data and environmental data to investigate species divergence and the maintenance of species boundaries in these two species. We found evidence of advanced-generation asymmetric hybridization and introgression, with shoal chub alleles introgressing more frequently into prairie chubs than the reciprocal. Using a Bayesian Genomic Cline framework, patterns of genomic introgression were revealed to be quite heterogeneous, yet shoal chub alleles were found to have likely selectively introgressed across species boundaries significantly more often than prairie chub alleles, potentially explaining some of the observed asymmetry in hybridization. These patterns were remarkably consistent across two sampled geographic regions of hybridization. Several environmental variables were found to significantly predict individual admixture, suggesting ecological isolation might maintain species boundaries.

The genetic legacy of 50 years of desert bighorn sheep translocations

Conservation biologists have increasingly used translocations to mitigate population declines and restore locally extirpated populations. Genetic data can guide the selection of source populations for translocations and help evaluate restoration success. Bighorn sheep (Ovis canadensis) are a managed big game species that suffered widespread population extirpations across western North America throughout the early 1900s. Subsequent translocation programs have successfully re-established many formally extirpated bighorn herds, but most of these programs pre-date genetically informed management practices. The state of Nevada presents a particularly well-documented case of decline followed by restoration of extirpated herds. Desert bighorn sheep (O. c. nelsoni) populations declined to less than 3,000 individuals restricted to remnant herds in the Mojave Desert and a few locations in the Great Basin Desert. Beginning in 1968, the Nevada Department of Wildlife translocated ~2,000 individuals from remnant populations to restore previously extirpated areas, possibly establishing herds with mixed ancestries. Here, we examined genetic diversity and structure among remnant herds and the genetic consequences of translocation from these herds using a genotyping-by-sequencing approach to genotype 17,095 loci in 303 desert bighorn sheep. We found a signal of population genetic structure among remnant Mojave Desert populations, even across geographically proximate mountain ranges. Further, we found evidence of a genetically distinct, potential relict herd from a previously hypothesized Great Basin lineage of desert bighorn sheep. The genetic structure of source herds was clearly reflected in translocated populations. In most cases, herds retained genetic evidence of multiple translocation events and subsequent admixture when founded from multiple remnant source herds. Our results add to a growing literature on how population genomic data can be used to guide and monitor restoration programs.

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.

The relationship between postmating reproductive isolation and reinforcement in Phlox

The process of speciation involves the accumulation of reproductive isolation (RI) between diverging lineages. Selection can favor increased RI via the process of reinforcement, whereby costs to hybridization impose selection for increased prezygotic RI. Reinforcement results in phenotypic divergence within at least one taxon, as a result of costly hybridization between sympatric taxa. The strength of selection driving reinforcement is determined by the cost of hybridization and the frequency of hybridization. We investigated the cost of hybridization by quantifying postmating RI barriers among Phlox species that comprise one of the best-studied cases of reinforcement. We determined if the strength of RI differs among lineages that have and have not undergone reinforcement, how much variability there is within species in RI, and whether RI is associated with phylogenetic relatedness. We found high RI for the species that underwent phenotypic divergence due to reinforcement; however, RI was also high between other species pairs. We found extensive variability in RI among individuals within species, and no evidence that the strength of RI was associated with phylogenetic relatedness. We suggest that phenotypic divergence due to reinforcement is associated with the frequency of hybridization and introgression, and not the cost of hybridization in this clade.

Admixture on the northern front: population genomics of range expansion in the white-footed mouse (Peromyscus leucopus) and secondary contact with the deer mouse (Peromyscus maniculatus)

Range expansion has genetic consequences expected to result in differentiated wave-front populations with low genetic variation and potentially introgression from a local species. The northern expansion of Peromyscus leucopus in southern Quebec provides an opportunity to test these predictions using population genomic tools. Our results show evidence of recent and post-glacial expansion. Genome-wide variation in P. leucopus indicates two post-glacial lineages are separated by the St. Lawrence River, with a more recent divergence of populations isolated by the Richelieu River. In two of three transects we documented northern populations with low diversity in at least one genetic measure, although most relationships were not significant. Consistent with bottlenecks and allele surfing during northward expansion, we document a northern-most population with low nucleotide diversity, divergent allele frequencies and the most private alleles, and observed heterozygosity indicates outcrossing. Ancestry proportions revealed putative hybrids of P. leucopus and P. maniculatus. A formal test for gene flow confirmed secondary contact, showing that a reticulate population phylogeny between P. maniculatus and P. leucopus was a better fit to the data than a bifurcating model without gene flow. Thus, we provide the first genomic evidence of gene flow between this pair of species in natural populations. Understanding the evolutionary consequences of secondary contact is an important conservation concern as climate-induced range expansions are expected to result in new hybrid zones between closely related species.

Inconsistent reproductive isolation revealed by interactions between Catostomus fish species

Interactions between species are central to evolution and ecology, but we do not know enough about how outcomes of interactions between species vary across geographic locations, in heterogeneous environments, or over time. Ecological dimensions of interactions between species are known to vary, but evolutionary interactions such as the establishment and maintenance of reproductive isolation are often assumed to be consistent across instances of an interaction between species. Hybridization among Catostomus fish species occurs over a large and heterogeneous geographic area and across taxa with distinct evolutionary histories, which allows us to assess consistency in species interactions. We analyzed hybridization among six Catostomus species across the Upper Colorado River basin (US mountain west) and found extreme variation in hybridization across locations. Different hybrid crosses were present in different locations, despite similar species assemblages. Within hybrid crosses, hybridization varied from only first generation hybrids to extensive hybridization with backcrossing. Variation in hybridization outcomes might result from uneven fitness of hybrids across locations, polymorphism in genetic incompatibilities, chance, unidentified historical contingencies, or some combination thereof. Our results suggest caution in assuming that one or a few instances of hybridization represent all interactions between the focal species, as species interactions vary substantially across locations.

Advances in ecological genomics in forest trees and applications to genetic resources conservation and breeding

Forest trees are an unparalleled group of organisms in their combined ecological, economic and societal importance. With widespread distributions, predominantly random mating systems and large population sizes, most tree species harbour extensive genetic variation both within and among populations. At the same time, demographic processes associated with Pleistocene climate oscillations and land-use change have affected contemporary range-wide diversity and may impinge on the potential for future adaptation. Understanding how these adaptive and neutral processes have shaped the genomes of trees species is therefore central to their management and conservation. As for many other taxa, the advent of high-throughput sequencing methods is expected to yield an understanding of the interplay between the genome and environment at a level of detail and depth not possible only a few years ago. An international conference entitled 'Genomics and Forest Tree Genetics' was held in May 2016, in Arcachon (France), and brought together forest geneticists with a wide range of research interests to disseminate recent efforts that leverage contemporary genomic tools to probe the population, quantitative and evolutionary genomics of trees. An important goal of the conference was to discuss how such data can be applied to both genome-enabled breeding and the conservation of forest genetic resources under land use and climate change. Here, we report discoveries presented at the meeting and discuss how the ecological genomic toolkit can be used to address both basic and applied questions in tree biology.

Studying the genetic basis of speciation in high gene flow marine invertebrates

A growing number of genes responsible for reproductive incompatibilities between species (barrier loci) exhibit the signals of positive selection. However, the possibility that genes experiencing positive selection diverge early in speciation and commonly cause reproductive incompatibilities has not been systematically investigated on a genome-wide scale. Here, I outline a research program for studying the genetic basis of speciation in broadcast spawning marine invertebrates that uses a priori genome-wide information on a large, unbiased sample of genes tested for positive selection. A targeted sequence capture approach is proposed that scores single-nucleotide polymorphisms (SNPs) in widely separated species populations at an early stage of allopatric divergence. The targeted capture of both coding and non-coding sequences enables SNPs to be characterized at known locations across the genome and at genes with known selective or neutral histories. The neutral coding and non-coding SNPs provide robust background distributions for identifying FST-outliers within genes that can, in principle, identify specific mutations experiencing diversifying selection. If natural hybridization occurs between species, the neutral coding and non-coding SNPs can provide a neutral admixture model for genomic clines analyses aimed at finding genes exhibiting strong blocks to introgression. Strongylocentrotid sea urchins are used as a model system to outline the approach but it can be used for any group that has a complete reference genome available.

What, if anything, are hybrids: enduring truths and challenges associated with population structure and gene flow

Hybridization is a potent evolutionary process that can affect the origin, maintenance, and loss of biodiversity. Because of its ecological and evolutionary consequences, an understanding of hybridization is important for basic and applied sciences, including conservation biology and agriculture. Herein, we review and discuss ideas that are relevant to the recognition of hybrids and hybridization. We supplement this discussion with simulations. The ideas we present have a long history, particularly in botany, and clarifying them should have practical consequences for managing hybridization and gene flow in plants. One of our primary goals is to illustrate what we can and cannot infer about hybrids and hybridization from molecular data; in other words, we ask when genetic analyses commonly used to study hybridization might mislead us about the history or nature of gene flow and selection. We focus on patterns of variation when hybridization is recent and populations are polymorphic, which are particularly informative for applied issues, such as contemporary hybridization following recent ecological change. We show that hybridization is not a singular process, but instead a collection of related processes with variable outcomes and consequences. Thus, it will often be inappropriate to generalize about the threats or benefits of hybridization from individual studies, and at minimum, it will be important to avoid categorical thinking about what hybridization and hybrids are. We recommend potential sampling and analytical approaches that should help us confront these complexities of hybridization.