Physiological adaptation along environmental gradients and replicated hybrid zone structure in swordtails (Teleostei: Xiphophorus)

Hybridization in the Anthropocene - how pollution and climate change disrupt mate selection in freshwater fish

Chemical pollutants and/or climate change have the potential to break down reproductive barriers between species and facilitate hybridization. Hybrid zones may arise in response to environmental gradients and secondary contact between formerly allopatric populations, or due to the introduction of non-native species. In freshwater ecosystems, field observations indicate that changes in water quality and chemistry, due to pollution and climate change, are correlated with an increased frequency of hybridization. Physical and chemical disturbances of water quality can alter the sensory environment, thereby affecting chemical and visual communication among fish. Moreover, multiple chemical compounds (e.g. pharmaceuticals, metals, pesticides, and industrial contaminants) may impair fish physiology, potentially affecting phenotypic traits relevant for mate selection (e.g. pheromone production, courtship, and coloration). Although warming waters have led to documented range shifts, and chemical pollution is ubiquitous in freshwater ecosystems, few studies have tested hypotheses about how these stressors may facilitate hybridization and what this means for biodiversity and species conservation. Through a systematic literature review across disciplines (i.e. ecotoxicology and evolutionary biology), we evaluate the biological interactions, toxic mechanisms, and roles of physical and chemical environmental stressors (i.e. chemical pollution and climate change) in disrupting mate preferences and inducing interspecific hybridization in freshwater fish. Our study indicates that climate change-driven changes in water quality and chemical pollution may impact visual and chemical communication crucial for mate choice and thus could facilitate hybridization among fishes in freshwater ecosystems. To inform future studies and conservation management, we emphasize the importance of further research to identify the chemical and physical stressors affecting mate choice, understand the mechanisms behind these interactions, determine the concentrations at which they occur, and assess their impact on individuals, populations, species, and biological diversity in the Anthropocene.

Ecological niche contributes to the persistence of the western × glaucous-winged gull hybrid zone

Hybrid zones occur in nature when populations with limited reproductive barriers overlap in space. Many hybrid zones persist over time, and different models have been proposed to explain how selection can maintain hybrid zone stability. More empirical studies are needed to elucidate the role of ecological adaptation in maintaining stable hybrid zones. Here, we investigated the role of exogenous factors in maintaining a hybrid zone between western gulls (Larus occidentalis) and glaucous-winged gulls (L. glaucescens). We used ecological niche models (ENMs) and niche similarity tests to quantify and examine the ecological niches of western gulls, glaucous-winged gulls, and their hybrids. We found evidence of niche divergence between all three groups. Our results support the bounded superiority model, providing further evidence that exogenous selection favoring hybrids may be an important factor in maintaining this stable hybrid zone.

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.

Validity of Xiphophorus fish as models for human disease

Platyfish and swordtails of the genus Xiphophorus provide a well-established model for melanoma research and have become well known for this feature. Recently, modelling approaches for other human diseases in Xiphophorus have been developed or are emerging. This Review provides a comprehensive summary of these models and discusses how findings from basic biological and molecular studies and their translation to medical research demonstrate that Xiphophorus models have face, construct and predictive validity for studying a broad array of human diseases. These models can thus improve our understanding of disease mechanisms to benefit patients.

Ecological niche contributes to the persistence of the Western x Glaucous-winged Gull hybrid zone

Hybrid zones occur in nature when populations with limited reproductive barriers overlap in space. Many hybrid zones persist over time, and different models have been proposed to explain how selection can maintain hybrid zone stability. More empirical studies are needed to elucidate the role of ecological adaptation in maintaining stable hybrid zones. Here, we investigated the role of exogenous factors in maintaining a hybrid zone between western gulls (Larus occidentalis) and glaucous-winged gulls (L. glaucescens). We used ecological niche models (ENMs) and niche similarity tests to quantify and examine the ecological niches of western gulls, glaucous-winged gulls, and their hybrids. We found evidence of niche divergence between all three groups. Our results best support the bounded superiority model, providing further evidence that exogenous selection favoring hybrids may be an important factor in maintaining this stable hybrid zone.

Investigating the utility of Anchored Hybrid Enrichment data to investigate the relationships among the Killifishes (Actinopterygii: Cyprinodontiformes), a globally distributed group of fishes

The Killifishes (Cyprinodontiformes) are a diverse and well-known group of fishes that contains sixteen families inclusive of Anablepidae, Aphaniidae Aplocheilidae, Cubanichthyidae, Cyprinodontidae, Fluviphylacidae, Fundulidae, Goodeidae, Nothobranchiidae, Orestiidae, Pantanodontidae, Poeciliidae, Procatopodidae, Profundulidae, Rivulidae, and Valenciidae and more than 1,200 species that are globally distributed in tropical and temperate, freshwater and estuarine habitats. The evolutionary relationships among the families within the group, based on different molecular and morphological data sets, have remained uncertain. Therefore, the objective of this study was to use a targeted approach, anchored hybrid enrichment, to investigate the phylogenetic relationships among the families within the Cyprindontiformes. This study included more than 100 individuals, representing all sixteen families within the Cyprinodontiformes, including many recently diagnosed families. We recovered an average of 244 loci per individual. These data were submitted to phylogenetic analyses (RaxML and ASTRAL) and although we recovered many of the same relationships as in previous studies of the group, several novel sets of relationships for other families also were recovered. In addition, two well-established clades (Suborders Cyprinodontoidei and Aplocheilodei) were recovered as monophyletic and are in agreement with most previous studies. We also assessed the degree of gene tree discordance in our dataset to evaluate support for alternative topological hypotheses for interfamilial relationships within the Cyprinodontiformes using a variety of different analyses. The results from this study will provide a robust, historical framework needed to investigate a plethora of biogeographic, taxonomic, ecological, and physiological questions for this group of fishes.

Do habitat and elevation promote hybridization during secondary contact between three genetically distinct groups of warbling vireo (Vireo gilvus)?

Following postglacial expansion, secondary contact can occur between genetically distinct lineages. These genetic lineages may be associated with specific habitat or environmental variables and therefore, their distributions in secondary contact could reflect such conditions within these areas. Here we used mtDNA, microsatellite, and morphological data to study three genetically distinct groups of warbling vireo (Vireo gilvus) and investigate the role that elevation and habitat play in their distributions. We studied two main contact zones and within each contact zone, we examined two separate transects. Across the Great Plains contact zone, we found that hybridization between eastern and western groups occurs along a habitat and elevational gradient, whereas hybridization across the Rocky Mountain contact zone was not as closely associated with habitat or elevation. Hybrids in the Great Plains contact zone were more common in transitional areas between deciduous and mixed-wood forests, and at lower elevations (<1000 m). Hybridization patterns were similar along both Great Plains transects indicating that habitat and elevation play a role in hybridization between distinct eastern and western genetic groups. The observed patterns suggest adaptation to different habitats, perhaps originating during isolation in multiple Pleistocene refugia, is facilitating hybridization in areas where habitat types overlap.

Cognitive Phenotypic Plasticity: Environmental Enrichment Affects Learning but Not Executive Functions in a Teleost Fish, Poecilia reticulata

Many aspects of animal cognition are plastically adjusted in response to the environment through individual experience. A remarkable example of this cognitive phenotypic plasticity is often observed when comparing individuals raised in a barren environment to individuals raised in an enriched environment. Evidence of enrichment-driven cognitive plasticity in teleost fish continues to grow, but it remains restricted to a few cognitive traits. The purpose of this study was to investigate how environmental enrichment affects multiple cognitive traits (learning, cognitive flexibility, and inhibitory control) in the guppy, Poecilia reticulata. To reach this goal, we exposed new-born guppies to different treatments: an enrichment environment with social companions, natural substrate, vegetation, and live prey or a barren environment with none of the above. After a month of treatment, we tested the subjects in a battery of three cognitive tasks. Guppies from the enriched environment learned a color discrimination faster compared to guppies from the environment with no enrichments. We observed no difference between guppies of the two treatments in the cognitive flexibility task, requiring selection of a previously unrewarded stimulus, nor in the inhibitory control task, requiring the inhibition of the attack response toward live prey. Overall, the results indicated that environmental enrichment had an influence on guppies' learning ability, but not on the remaining cognitive functions investigated.

Genomic data and multi-species demographic modelling uncover past hybridization between currently allopatric freshwater species

Evidence for ancient interspecific gene flow through hybridization has been reported in many animal and plant taxa based on genetic markers. The study of genomic patterns of closely related species with allopatric distributions allows the assessment of the relative importance of vicariant isolating events and past gene flow. Here, we investigated the role of gene flow in the evolutionary history of four closely related freshwater fish species with currently allopatric distributions in western Iberian rivers-Squalius carolitertii, S. pyrenaicus, S. torgalensis and S. aradensis-using a population genomics dataset of 23,562 SNPs from 48 individuals, obtained through genotyping by sequencing (GBS). We uncovered a species tree with two well-differentiated clades: (i) S. carolitertii and S. pyrenaicus; and (ii) S. torgalensis and S. aradensis. By using D-statistics and demographic modelling based on the site frequency spectrum, comparing alternative demographic scenarios of hybrid origin, secondary contact and isolation, we found that the S. pyrenaicus North lineage is likely the result of an ancient hybridization event between S. carolitertii (contributing ~84%) and S. pyrenaicus South lineage (contributing ~16%), consistent with a hybrid speciation scenario. Furthermore, in the hybrid lineage, we identify outlier loci potentially affected by selection favouring genes from each parental lineage at different genomic regions. Our results suggest that ancient hybridization can affect speciation and that freshwater fish species currently in allopatry are useful to study these processes.

Differences in Thermal Tolerance between Parental Species Could Fuel Thermal Adaptation in Hybrid Wood Ants

AbstractGenetic variability is essential for adaptation and could be acquired via hybridization with a closely related lineage. We use ants to investigate thermal adaptation and the link between temperature and genetic variation arising from hybridization. We test for differences in cold and heat tolerance between Finnish Formica polyctena and Formica aquilonia wood ants and their naturally occurring hybrids. Using workers, we find that the parental individuals differ in both cold and heat tolerances and express thermal limits that reflect their global distributions. Hybrids, however, cannot combine thermal tolerance of parental species as they have the same heat tolerance as F. polyctena but not the same cold tolerance as F. aquilonia. We then focus on a single hybrid population to investigate the relationship between temperature variation and genetic variation across 16 years using reproductive individuals. On the basis of the thermal tolerance results, we expected the frequency of putative F. polyctena alleles to increase in warm years and F. aquilonia alleles to increase in cold years. We find support for this in hybrid males but not in hybrid females. These results contribute to understanding the outcomes of hybridization, which may be sex specific or depend on the environment. Furthermore, genetic variability resulting from hybridization could help hybrid wood ants cope with changing thermal conditions.

Population genetic structure of sharpbelly Hemiculter leucisculus (Basilesky, 1855) and morphological diversification along climate gradients in China

Sharpbelly Hemiculter leucisculus (Basilewski, 1855) is a small, widespread, and native cyprinid fish with prominent habitat suitability and high invasive potential and is becoming the dominant species in freshwater ecosystems under intensified environmental disturbances. But how H. leucisculus acclimates to extremely heterogeneous environments remains unclear. In current study, the genetic structure of H. leucisculus was analyzed using Bayesian phylogenetic inference, haplotype network, and STRUCTURE base on cytb gene across 18 populations spanning 20 degrees of latitude and 18 degrees of longitude in China. The morphological diversification of body size and shape for H. leucisculus along the climate gradient was studied. The results showed that the 18 H. leucisculus populations were divided into 3 clusters: one cluster mainly from Huanghe River Basin, another cluster mainly from Yangzi River Basin, and H cluster containing Hainan and Beihai populations. The fish from southern populations were deeper bodied while individuals from northern populations were more slender. Inland individuals were more streamlined while coastal individuals were of deeper body. The partial Mantel test predicts that the potential mechanism underlining the intraspecies morphological diversification along climate gradients is primarily the divergent selection pressures among different environments, while genetic variation had less contribution to morphological differentiation. The formation of the Nanling Mountain Range could drive genetic differentiation between Beihai population and those from Yangzi River Basin. The present results highlight strong selective pressures of climate on widespread species and enrich morphological differentiation basis of acclimation for species with high habitat suitability and invasive potential.

Hybridization alters growth and migratory life-history expression of native trout

Human-mediated hybridization threatens many native species, but the effects of introgressive hybridization on life-history expression are rarely quantified, especially in vertebrates. We quantified the effects of non-native rainbow trout admixture on important life-history traits including growth and partial migration behavior in three populations of westslope cutthroat trout over five years. Rainbow trout admixture was associated with increased summer growth rates in all populations and decreased spring growth rates in two populations with cooler spring temperatures. These results indicate that non-native admixture may increase growth under warmer conditions, but cutthroat trout have higher growth rates during cooler periods. Non-native admixture consistently increased expression of migratory behavior, suggesting that there is a genomic basis for life-history differences between these species. Our results show that effects of interspecific hybridization on fitness traits can be the product of genotype-by-environment interactions even when there are minor differences in environmental optima between hybridizing species. These results also indicate that while environmentally mediated traits like growth may play a role in population-level consequences of admixture, strong genetic influences on migratory life-history differences between these species likely explains the continued spread of non-native hybridization at the landscape-level, despite selection against hybrids at the population-level.

The role of the environment in the spatial dynamics of an extensive hybrid zone between two neotropical cats

Identifying factors that create and maintain a hybrid zone is of great interest to ecology, evolution and, more recently, conservation biology. Here, we investigated the role of environmental features in shaping the spatial dynamics of a hybrid zone between the southern tigrina, Leopardus guttulus, and Geoffroy's cat, L. geoffroyi, testing for exogenous selection as the main force acting on its maintenance. These Neotropical felid species are mainly allopatric, with a restricted area of sympatry in the ecotone between the Atlantic Forest and Pampa biomes. As both biomes have experienced high rates of anthropogenic habitat alteration, we also analysed the influence of habitat conversion on the hybrid zone structure. To do this, we used 13 microsatellite loci to identify potential hybrids and generated ecological niche models for them and their parental species. We compared the influence of variables on parental species and hybrid occurrence and calculated the amount of niche overlap among them. Parental species showed different habitat requirements and predicted co-occurrence was restricted to the forest-grassland mosaic of the ecotone. However, hybrids were found beyond this area, mainly in the range of L. geoffroyi. Hybrids demonstrated higher tolerance to habitat alteration than parental types, with a probability of occurrence that was positively related with mosaics of cropland areas and remnants of natural vegetation. These results indicate that exogenous selection alone does not drive the dynamics of the hybrid zone, and that habitat conversion influences its structure, potentially favouring hybrids over parental species.

Patterns of introgression vary within an avian hybrid zone

Background

Exploring hybrid zone dynamics at different spatial scales allows for better understanding of local factors that influence hybrid zone structure. In this study, we tested hypotheses about drivers of introgression at two spatial scales within the Saltmarsh Sparrow (Ammospiza caudacuta) and Nelson’s Sparrow (A. nelsoni) hybrid zone. Specifically, we evaluated the influence of neutral demographic processes (relative species abundance), natural selection (exogenous environmental factors and genetic incompatibilities), and sexual selection (assortative mating) in this mosaic hybrid zone. By intensively sampling adults (n = 218) and chicks (n = 326) at two geographically proximate locations in the center of the hybrid zone, we determined patterns of introgression on a fine scale across sites of differing habitat. We made broadscale comparisons of patterns from the center with those of prior studies in the southern edge of the hybrid zone.

Results

A panel of fixed SNPs (135) identified from ddRAD sequencing was used to calculate a hybrid index and determine genotypic composition/admixture level of the populations. Another panel of polymorphic SNPs (589) was used to assign paternity and reconstruct mating pairs to test for sexual selection. On a broad-scale, patterns of introgression were not explained by random mating within marshes. We found high rates of back-crossing and similarly low rates of recent-generation (F1/F2) hybrids in the center and south of the zone. Offspring genotypic proportions did not meet those expected from random mating within the parental genotypic distribution. Additionally, we observed half as many F1/F2 hybrid female adults than nestlings, while respective male groups showed no difference, in support of Haldane’s Rule. The observed proportion of interspecific mating was lower than expected when accounting for mate availability, indicating assortative mating was limiting widespread hybridization. On a fine spatial scale, we found variation in the relative influence of neutral and selective forces between inland and coastal habitats, with the smaller, inland marsh influenced primarily by neutral demographic processes, and the expansive, coastal marsh experiencing higher selective pressures in the form of natural (exogenous and endogenous) and sexual selection.

Conclusions

Multiple drivers of introgression, including neutral and selective pressures (exogenous, endogenous, and sexual selection), are structuring this hybrid zone, and their relative influence is site and context-dependent.

Salinity, not genetic incompatibilities, limits the establishment of the invasive hybrid cattail Typha × glauca in coastal wetlands

Hybrids of a single pair of parent species can be much more common in some geographical regions than in others. The reasons for this are not well understood, but could help explain processes such as species diversification or the range expansion of invasive hybrids. The widespread cattails Typha latifolia and T. angustifolia seldom hybridize in some parts of their range, but in other areas produce the dominant hybrid T. × glauca. We used a combination of field and greenhouse experiments to investigate why T. × glauca has invaded wetlands in the Laurentian Great Lakes region of southern Ontario, Canada, but is much less common in the coastal wetlands of Nova Scotia (NS) in eastern Canada. One potentially important environmental difference between these two regions is salinity. We therefore tested three hypotheses: (1) T. latifolia and T. angustifolia in NS are genetically incompatible; (2) the germination or growth of T. × glauca is reduced by salinity; and (3) T. latifolia, a main competitor of T. × glauca, is locally adapted to saline conditions in NS. Our experiments showed that NS T. latifolia and T. angustifolia are genetically compatible, and that saline conditions do not impede growth of hybrid plants. However, we also found that under conditions of high salinity, germination rates of hybrid seeds were substantially lower than those of NS T. latifolia. In addition, germination rates of NS T. latifolia were higher than those of Ontario T. latifolia, suggesting local adaptation to salinity in coastal wetlands. This study adds to the growing body of literature which identifies the important roles that local habitat and adaptation can play in the distributions and characteristics of hybrid zones.

Temperature effects on performance and physiology of two prairie stream minnows

Earth's atmosphere has warmed by ~1°C over the past century and continues to warm at an increasing rate. Effects of atmospheric warming are already visible in most major ecosystems and are evident across all levels of biological organization. Linking functional responses of individuals to temperature is critical for predicting responses of populations and communities to global climate change. The southern redbelly dace Chrosomus erythrogaster and the central stoneroller Campostoma anomalum are two minnows (Cyprinidae) that commonly occur in the Flint Hills region of the USA but show different patterns of occurrence, with dace largely occupying headwater reaches and stonerollers persisting in both headwater and intermediate-sized streams. We tested for differences between species in critical thermal maximum, energy metabolism, sustained swimming and activity over an ecologically relevant temperature gradient of acclimation temperatures. Typically, metrics increased with acclimation temperature for both species, although stoneroller activity decreased with temperature. We observed a significant interaction between species and temperature for critical thermal maxima, where stonerollers only had higher critical thermal maxima at the coldest temperature and at warm temperatures compared to the dace. We did not find evidence suggesting differences in the energy metabolism of dace and stonerollers. We detected interspecific differences in sustained swimming performance, with dace having higher swimming speed than stonerollers regardless of acclimation temperature. Finally, there was a significant interaction between temperature and species for activity; dace activity was higher at intermediate and warm temperatures compared to stonerollers. We observed subtle interspecific differences in how performance metrics responded to temperature that did not always align with observed patterns of distribution for these species. Thus, other ecological factors likely are important drivers of distributional patterns in these species.

Natural and sexual selection drive multivariate phenotypic divergence along climatic gradients in an invasive fish

Invasive species that rapidly spread throughout novel distribution ranges are prime models to investigate climate-driven phenotypic diversification on a contemporary scale. Previous studies on adaptive diversification along latitudinal gradients in fish have mainly considered body size and reported either increased or decreased body size towards higher latitudes (i.e. Bergmann's rule). Our study is the first to investigate phenotypic divergence in multiple traits, including sexually selected traits (size and shape of the male copulatory organ, the gonopodium) of invasive Gambusia affinis in China. We studied body size, life history traits and morphological variation across populations spanning 17 degrees of latitude and 16 degrees of longitude. Even though we found phenotypic variation along climatic gradients to be strongest in naturally selected traits, some sexually selected traits also showed systematic gradual divergence. For example, males from southern populations possessed wider gonopodia with increased armament. Generally, males and females diverged in response to different components of climatic gradients (latitudinal or longitudinal variation) and in different trait suites. We discuss that not only temperature regimes, but also indirect effects of increased resource and mate competition (as a function of different extrinsic overwinter mortality rates) alter the selective landscape along climatic gradients.

Natural selection interacts with recombination to shape the evolution of hybrid genomes

To investigate the consequences of hybridization between species, we studied three replicate hybrid populations that formed naturally between two swordtail fish species, estimating their fine-scale genetic map and inferring ancestry along the genomes of 690 individuals. In all three populations, ancestry from the "minor" parental species is more common in regions of high recombination and where there is linkage to fewer putative targets of selection. The same patterns are apparent in a reanalysis of human and archaic admixture. These results support models in which ancestry from the minor parental species is more likely to persist when rapidly uncoupled from alleles that are deleterious in hybrids. Our analyses further indicate that selection on swordtail hybrids stems predominantly from deleterious combinations of epistatically interacting alleles.

Population genetic dynamics of Himalayan-Hengduan tree peonies, Paeonia subsect. Delavayanae

According to the present taxonomical treatment, Paeonia subsect. Delavayanae consists of only two species (P. delavayi and P. ludlowii) endemic to the Himalayan-Hengduan Mountains. Although P. ludlowii can be distinguished from P. delavayi on the basis of a series of morphological characters, the species delimitation remains controversial because the more widespread one, P. delavayi, exhibits considerable morphological diversity. Both chloroplast DNA markers and nuclear microsatellites or simple sequence repeats (nSSR) are used herein to reveal genetic diversity and relationships of the two taxa included in this subsection, and ecological niche modeling (ENM) is employed to get insights into their paleodistribution. Our results show that genetic boundaries between the two currently recognized species are unclear, probably due to recent divergence. Paeonia ludlowii is budding from P. delavayi, probably by genetic isolation but also by shifting its niche to the harsher upland Tibetan conditions. Paeonia delavayi itself would be, however, under active speciation, showing significant genetic differentiation and morphological diversity. Whereas P. ludlowii would have endured the Pleistocene glacial periods by in situ persistence in local, small refugia, a 'dual' model seems to apply for P. delavayi (in situ persistence and retreat to refugia). The rarity of P. ludlowii and high evolutionary potential of P. delavayi imply high priority for in situ conservation of both taxa. The Himalayan-Hengduan Mountains are an ideal arena for differentiation within subsect. Delavayanae of Paeonia, by means of expansions/contractions/displacements, vertical migrations, and local survival/extinctions in response to the Neogene climate fluctuations and geological changes.

Low temperature and low salinity drive putatively adaptive growth differences in populations of threespine stickleback

Colonisation can expose organisms to novel combinations of abiotic and biotic factors and drive adaptive divergence. Yet, studies investigating the interactive effects of multiple abiotic factors on the evolution of physiological traits remain rare. Here we examine the effects of low salinity, low temperature, and their interaction on the growth of three North American populations of threespine stickleback (Gasterosteus aculeatus). In north-temperate freshwater habitats, stickleback populations experience a combination of low salinity and low winter temperatures that are not experienced by the ancestral marine and anadromous populations. Here we show that both salinity and temperature, and their interaction, have stronger negative effects on marine and anadromous populations than a freshwater population. Freshwater stickleback showed only a ~20% reduction in specific growth rate when exposed to 4 °C, while marine and anadromous stickleback showed sharp declines (82% and 74% respectively) under these conditions. The modest decreases in growth in freshwater stickleback in fresh water in the cold strongly suggest that this population has the capacity for physiological compensation to offset the negative thermodynamic effects of low temperature on growth. These results are suggestive of adaptive evolution in response to the interactive effects of low salinity and low temperature during freshwater colonisation.

Climate, Demography, and Zoogeography Predict Introgression Thresholds in Salmonid Hybrid Zones in Rocky Mountain Streams

Among the many threats posed by invasions of nonnative species is introgressive hybridization, which can lead to the genomic extinction of native taxa. This phenomenon is regarded as common and perhaps inevitable among native cutthroat trout and introduced rainbow trout in western North America, despite that these taxa naturally co-occur in some locations. We conducted a synthetic analysis of 13,315 genotyped fish from 558 sites by building logistic regression models using data from geospatial stream databases and from 12 published studies of hybridization to assess whether environmental covariates could explain levels of introgression between westslope cutthroat trout and rainbow trout in the U.S. northern Rocky Mountains. A consensus model performed well (AUC, 0.78–0.86; classification success, 72–82%; 10-fold cross validation, 70–82%) and predicted that rainbow trout introgression was significantly associated with warmer water temperatures, larger streams, proximity to warmer habitats and to recent sources of rainbow trout propagules, presence within the historical range of rainbow trout, and locations further east. Assuming that water temperatures will continue to rise in response to climate change and that levels of introgression outside the historical range of rainbow trout will equilibrate with those inside that range, we applied six scenarios across a 55,234-km stream network that forecast 9.5–74.7% declines in the amount of habitat occupied by westslope cutthroat trout populations of conservation value, but not the wholesale loss of such populations. We conclude that introgression between these taxa is predictably related to environmental conditions, many of which can be manipulated to foster largely genetically intact populations of westslope cutthroat trout and help managers prioritize conservation activities.

Unique evolutionary trajectories in repeated adaptation to hydrogen sulphide-toxic habitats of a neotropical fish (Poecilia mexicana)

Replicated ecological gradients are prime systems to study processes of molecular evolution underlying ecological divergence. Here, we investigated the repeated adaptation of the neotropical fish Poecilia mexicana to habitats containing toxic hydrogen sulphide (H2 S) and compared two population pairs of sulphide-adapted and ancestral fish by sequencing population pools of >200 individuals (Pool-Seq). We inferred the evolutionary processes shaping divergence and tested the hypothesis of increase of parallelism from SNPs to molecular pathways. Coalescence analyses showed that the divergence occurred in the face of substantial bidirectional gene flow. Population divergence involved many short, widely dispersed regions across the genome. Analyses of allele frequency spectra suggest that differentiation at most loci was driven by divergent selection, followed by a selection-mediated reduction of gene flow. Reconstructing allelic state changes suggested that selection acted mainly upon de novo mutations in the sulphide-adapted populations. Using a corrected Jaccard index to quantify parallel evolution, we found a negligible proportion of statistically significant parallel evolution of Jcorr  = 0.0032 at the level of SNPs, divergent genome regions (Jcorr  = 0.0061) and genes therein (Jcorr  = 0.0091). At the level of metabolic pathways, the overlap was Jcorr  = 0.2545, indicating increasing parallelism with increasing level of biological integration. The majority of pathways contained positively selected genes in both sulphide populations. Hence, adaptation to sulphidic habitats necessitated adjustments throughout the genome. The largely unique evolutionary trajectories may be explained by a high proportion of de novo mutations driving the divergence. Our findings favour Gould's view that evolution is often the unrepeatable result of stochastic events with highly contingent effects.

Slow climate velocities of mountain streams portend their role as refugia for cold-water biodiversity

The imminent demise of montane species is a recurrent theme in the climate change literature, particularly for aquatic species that are constrained to networks and elevational rather than latitudinal retreat as temperatures increase. Predictions of widespread species losses, however, have yet to be fulfilled despite decades of climate change, suggesting that trends are much weaker than anticipated and may be too subtle for detection given the widespread use of sparse water temperature datasets or imprecise surrogates like elevation and air temperature. Through application of large water-temperature databases evaluated for sensitivity to historical air-temperature variability and computationally interpolated to provide high-resolution thermal habitat information for a 222,000-km network, we estimate a less dire thermal plight for cold-water species within mountains of the northwestern United States. Stream warming rates and climate velocities were both relatively low for 1968-2011 (average warming rate = 0.101 °C/decade; median velocity = 1.07 km/decade) when air temperatures warmed at 0.21 °C/decade. Many cold-water vertebrate species occurred in a subset of the network characterized by low climate velocities, and three native species of conservation concern occurred in extremely cold, slow velocity environments (0.33-0.48 km/decade). Examination of aggressive warming scenarios indicated that although network climate velocities could increase, they remain low in headwaters because of strong local temperature gradients associated with topographic controls. Better information about changing hydrology and disturbance regimes is needed to complement these results, but rather than being climatic cul-de-sacs, many mountain streams appear poised to be redoubts for cold-water biodiversity this century.

Ecological divergence and conservatism: spatiotemporal patterns of niche evolution in a genus of livebearing fishes (Poeciliidae: Xiphophorus)

Background

Ecological factors often have a strong impact on spatiotemporal patterns of biodiversity. The integration of spatial ecology and phylogenetics allows for rigorous tests of whether speciation is associated with niche conservatism (constraints on ecological divergence) or niche divergence. We address this question in a genus of livebearing fishes for which the role of sexual selection in speciation has long been studied, but in which the potential role of ecological divergence during speciation has not been tested.

Results

By combining reconstruction of ancestral climate tolerances and disparity indices, we show that the earliest evolutionary split in Xiphophorus was associated with significant divergence for temperature variables. Niche evolution and present day niches were most closely associated with each species’ geographic distribution relative to a biogeographic barrier, the Trans-Mexican Volcanic Belt. Tests for similarity of the environmental backgrounds of closely related species suggested that the relative importance of niche conservatism and divergence during speciation varied among the primary clades of Xiphophorus. Closely related species in the two swordtail clades exhibited higher levels of niche overlap than expected given environmental background similarity indicative of niche conservatism. In contrast, almost all species of platyfish had significantly divergent niches compared to environmental backgrounds, which is indicative of niche divergence.

Conclusion

The results suggest that the relative importance of niche conservatism and divergence differed among the clades of Xiphophorus and that traits associated with niche evolution may be more evolutionarily labile in the platyfishes. Our results ultimately suggest that the taxonomic scale of tests for conservatism and divergence could greatly influence inferences of their relative importance in the speciation process.

Electronic supplementary material

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

Ecological differentiation and habitat unsuitability maintaining a ground beetle hybrid zone

Exogenous selection via interactions between organisms and environments may influence the dynamics of hybrid zones between species in multiple ways. Two major models of a hybrid zone allowed us to hypothesize that environmental conditions influence hybrid zone dynamics in two ways. In the first model, an environmental gradient determines the mosaic distribution at the boundary between ecologically differentiated species (mosaic hybrid zone model). In the second model, a patch of unsuitable habitat traps a hybrid zone between species whose hybrids are unfit (tension zone model). To test these, we examined the environmental factors influencing the spatial structure of a hybrid zone between the ground beetles Carabus maiyasanus and C. iwawakianus using GIS‐based quantification of environmental factors and a statistical comparison of species distribution models (SDMs). We determined that both of the hypothetical processes can be important in the hybrid zone. We detected interspecific differences in the environmental factors in presence localities and their relative contribution in SDMs. SDMs were not identical between species even within contact areas, but tended to be similar within the range of each species. These results suggest an association between environments and species, and provide evidence that ecological differentiation between species plays a role in the maintenance of the hybrid zone. Contact areas were characterized by a relatively high temperature, low precipitation, and high topological wetness. Thus, the contact areas were regarded as being located in an unsuitable habitat with a drier climate, where those populations are likely to occur in patches with limited precipitation concentrated. A comparison of spatial scales suggests that exogenous selection via environmental factors may be weaker than endogenous selection via genitalic incompatibility.

Patterns of hybridization among cutthroat trout and rainbow trout in northern Rocky Mountain streams

Introgressive hybridization between native and introduced species is a growing conservation concern. For native cutthroat trout and introduced rainbow trout in western North America, this process is thought to lead to the formation of hybrid swarms and the loss of monophyletic evolutionary lineages. Previous studies of this phenomenon, however, indicated that hybrid swarms were rare except when native and introduced forms of cutthroat trout co‐occurred. We used a panel of 86 diagnostic, single nucleotide polymorphisms to evaluate the genetic composition of 3865 fish captured in 188 locations on 129 streams distributed across western Montana and northern Idaho. Although introgression was common and only 37% of the sites were occupied solely by parental westslope cutthroat trout, levels of hybridization were generally low. Of the 188 sites sampled, 73% contained ≤5% rainbow trout alleles and 58% had ≤1% rainbow trout alleles. Overall, 72% of specimens were nonadmixed westslope cutthroat trout, and an additional 3.5% were nonadmixed rainbow trout. Samples from seven sites met our criteria for hybrid swarms, that is, an absence of nonadmixed individuals and a random distribution of alleles within the sample; most (6/7) were associated with introgression by Yellowstone cutthroat trout. In streams with multiple sites, upstream locations exhibited less introgression than downstream locations. We conclude that although the widespread introduction of nonnative trout within the historical range of westslope cutthroat trout has increased the incidence of introgression, sites containing nonadmixed populations of this taxon are common and broadly distributed.

Spatiotemporal environmental heterogeneity and the maintenance of the tailspot polymorphism in the variable platyfish (Xiphophorus variatus)

Genetic variation is critical for adaptive evolution. Despite its importance, there is still limited evidence in support of some prominent theoretical models explaining the maintenance of genetic polymorphism within populations. We examined 84 populations of Xiphophorus variatus, a livebearing fish with a genetic polymorphism associated with physiological performance, to test: (1) whether niche differentiation explains broad-scale maintenance of polymorphism, (2) whether polymorphism is maintained among populations by local adaptation and migration, or (3) whether heterogeneity in explicit environmental variables could be linked to levels of polymorphism within populations. We found no evidence of climatic niche differentiation that could generate or maintain broad geographic variation in polymorphism. Subsequently, hierarchical partitioning of genetic richness and partial mantel tests revealed that 76% of the observed genetic richness was partitioned within populations with no effect of geographic distance on polymorphism. These results strongly suggest a lack of migration-selection balance in the maintenance of polymorphism, and model selection confirmed a significant relationship between environmental heterogeneity and genetic richness within populations. Few studies have demonstrated such effects at this scale, and additional studies in other taxa should examine the generality of gene-by-environment interactions across populations to better understand the dynamics and scale of balancing selection.

Genotype-environment associations support a mosaic hybrid zone between two tidal marsh birds

Local environmental features can shape hybrid zone dynamics when hybrids are bounded by ecotones or when patchily distributed habitat types lead to a corresponding mosaic of genotypes. We investigated the role of marsh-level characteristics in shaping a hybrid zone between two recently diverged avian taxa - Saltmarsh (Ammodramus caudacutus) and Nelson's (A. nelsoni) sparrows. These species occupy different niches where allopatric, with caudacutus restricted to coastal marshes and nelsoni found in a broader array of wetland and grassland habitats and co-occur in tidal marshes in sympatry. We determined the influence of habitat types on the distribution of pure and hybrid sparrows and assessed the degree of overlap in the ecological niche of each taxon. To do this, we sampled and genotyped 305 sparrows from 34 marshes across the hybrid zone and from adjacent regions. We used linear regression to test for associations between marsh characteristics and the distribution of pure and admixed sparrows. We found a positive correlation between genotype and environmental variables with a patchy distribution of genotypes and habitats across the hybrid zone. Ecological niche models suggest that the hybrid niche was more similar to that of A. nelsoni and habitat suitability was influenced strongly by distance from coastline. Our results support a mosaic model of hybrid zone maintenance, suggesting a role for local environmental features in shaping the distribution and frequency of pure species and hybrids across space.

Environmental heterogeneity generates opposite gene-by-environment interactions for two fitness-related traits within a population

Theory predicts that environmental heterogeneity offers a potential solution to the maintenance of genetic variation within populations, but empirical evidence remains sparse. The live-bearing fish Xiphophorus variatus exhibits polymorphism at a single locus, with different alleles resulting in up to five distinct melanistic "tailspot" patterns within populations. We investigated the effects of heterogeneity in two ubiquitous environmental variables (temperature and food availability) on two fitness-related traits (upper thermal limits and body condition) in two different tailspot types (wild-type and upper cut crescent). We found gene-by-environment (G × E) interactions between tailspot type and food level affecting upper thermal limits (UTL), as well as between tailspot type and thermal environment affecting body condition. Exploring mechanistic bases underlying these G × E patterns, we found no differences between tailspot types in hsp70 gene expression despite significant overall increases in expression under both thermal and food stress. Similarly, there was no difference in routine metabolic rates between the tailspot types. The reversal of relative performance of the two tailspot types under different environmental conditions revealed a mechanism by which environmental heterogeneity can balance polymorphism within populations through selection on different fitness-related traits.

Multiple mating and reproductive skew in parental and introgressed females of the live-bearing fish Xiphophorus birchmanni

Just as mating patterns can promote speciation or hybridization, the presence of hybridization can shape mating patterns within a population. In this study, we characterized patterns of multiple mating and reproductive skew in a naturally hybridizing swordtail fish species, Xiphophorus birchmanni. We quantified multiple mating using microsatellite markers to genotype embryos from 43 females collected from 2 wild populations. We also used a suite of single-nucleotide polymorphism markers to categorize females and their inferred mates as either parental X. birchmanni or as introgressed individuals, which carried alleles from a sister species, X. malinche. We found that parental and introgressed X. birchmanni females mated multiply with both parental and introgressed males. We found no difference in mating patterns or reproductive skew between parental and introgressed X. birchmanni females. However, nonintrogressed X. birchmanni males mated more often with large, fecund females. These females also had the greatest levels of skew in fertilization success of males. Thus, our results show that X. birchmanni has a polygynandrous mating system and that introgression of X. malinche alleles has only subtle effects on mating patterns in this species.

Adaptive radiation along a thermal gradient: preliminary results of habitat use and respiration rate divergence among whitefish morphs

Adaptive radiation is considered an important mechanism for the development of new species, but very little is known about the role of thermal adaptation during this process. Such adaptation should be especially important in poikilothermic animals that are often subjected to pronounced seasonal temperature variation that directly affects metabolic function. We conducted a preliminary study of individual lifetime thermal habitat use and respiration rates of four whitefish (Coregonus lavaretus (L.)) morphs (two pelagic, one littoral and one profundal) using stable carbon and oxygen isotope values of otolith carbonate. These morphs, two of which utilized pelagic habitats, one littoral and one profundal recently diverged via adaptive radiation to exploit different major niches in a deep and thermally stratified subarctic lake. We found evidence that the morphs used different thermal niches. The profundal morph had the most distinct thermal niche and consistently occupied the coldest thermal habitat of the lake, whereas differences were less pronounced among the shallow water pelagic and littoral morphs. Our results indicated ontogenetic shifts in thermal niches: juveniles of all whitefish morphs inhabited warmer ambient temperatures than adults. According to sampling of the otolith nucleus, hatching temperatures were higher for benthic compared to pelagic morphs. Estimated respiration rate was the lowest for benthivorous profundal morph, contrasting with the higher values estimated for the other morphs that inhabited shallower and warmer water. These preliminary results suggest that physiological adaptation to different thermal habitats shown by the sympatric morphs may play a significant role in maintaining or strengthening niche segregation and divergence in life-history traits, potentially contributing to reproductive isolation and incipient speciation.

Sexual ornaments, body morphology, and swimming performance in naturally hybridizing swordtails (teleostei: xiphophorus)

Determining the costs of sexual ornaments is complicated by the fact that ornaments are often integrated with other, non-sexual traits, making it difficult to dissect the effect of ornaments independent of other aspects of the phenotype. Hybridization can produce reduced phenotypic integration, allowing one to evaluate performance across a broad range of multivariate trait values. Here we assess the relationship between morphology and performance in the swordtails Xiphophorus malinche and X. birchmanni, two naturally-hybridizing fish species that differ extensively in non-sexual as well as sexual traits. We took advantage of novel trait variation in hybrids to determine if sexual ornaments incur a cost in terms of locomotor ability. For both fast-start and endurance swimming, hybrids performed at least as well as the two parental species. The sexually-dimorphic sword did not impair swimming performance per se. Rather, the sword negatively affected performance only when paired with a sub-optimal body shape. Studies seeking to quantify the costs of ornaments should consider that covariance with non-sexual traits may create the spurious appearance of costs.

Cold tolerance performance of westslope cutthroat trout (Oncorhynchus clarkii lewisi) and rainbow trout (Oncorhynchus mykiss) and its potential role in influencing interspecific hybridization

Hybridization between rainbow trout (Oncorhynchus mykiss (Walbaum, 1792)) and westslope cutthroat trout (Oncorhynchus clarkii lewisi (Girard, 1856)) occurs commonly when rainbow trout are introduced into the range of westslope cutthroat trout. Typically, hybridization is most common in warmer, lower elevation habitats, but much less common in colder, higher elevation habitats. We assessed the tolerance to cold water temperature (i.e., critical thermal minimum, CTMin) in juvenile rainbow trout and westslope cutthroat trout to test the hypothesis that westslope cutthroat trout better tolerate low water temperature, which may explain the lower prevalence of rainbow trout and interspecific hybrids in higher elevation, cold-water habitats (i.e., the “elevation refuge hypothesis”). All fish had significantly lower CTMin values (i.e., were better able to tolerate low temperatures) when they were acclimated to 15 °C (mean CTMin = 1.37 °C) versus 18 °C (mean CTMin = 1.91 °C; p < 0.001). Westslope cutthroat trout tended to have lower CTMin than rainbow trout from two populations, second–generation (F2) hybrids between two rainbow trout populations, and backcrossed rainbow trout at 15 °C (cross type × acclimation temperature interaction; p = 0.018). Differential adaptation to cold water temperatures may play a role in influencing the spatial distribution of hybridization between sympatric species of trout.

High-resolution mapping reveals hundreds of genetic incompatibilities in hybridizing fish species

Hybridization is increasingly being recognized as a common process in both animal and plant species. Negative epistatic interactions between genes from different parental genomes decrease the fitness of hybrids and can limit gene flow between species. However, little is known about the number and genome-wide distribution of genetic incompatibilities separating species. To detect interacting genes, we perform a high-resolution genome scan for linkage disequilibrium between unlinked genomic regions in naturally occurring hybrid populations of swordtail fish. We estimate that hundreds of pairs of genomic regions contribute to reproductive isolation between these species, despite them being recently diverged. Many of these incompatibilities are likely the result of natural or sexual selection on hybrids, since intrinsic isolation is known to be weak. Patterns of genomic divergence at these regions imply that genetic incompatibilities play a significant role in limiting gene flow even in young species.

Is there metabolic cold adaptation in terrestrial ectotherms? Exploring latitudinal compensation in the invasive snail Cornu aspersum

Lower temperatures, extreme seasonality and shorter growing seasons at higher latitudes are expected to cause a decline in metabolic rates and annual growth rates of ectotherms. If a reduction in the rates of these biological processes involves a reduction in fitness, then organisms may evolve compensatory responses for the constraints imposed by high-latitude habitats. To test the existence of a latitudinal compensation in ectotherms, we used a common-garden experiment to investigate the extent to which the level of energy turnover (measured as standard metabolic rate, SMR) and the energy budget (energy allocation to growth) are affected by climatic constraints in three populations of the land snail Cornu aspersum, distributed across a latitudinal gradient of 1300 km in Chile. Our results did not support the existence of a latitudinal compensation in metabolic rates (metabolic cold adaptation). However, there was a countergradient variation (CnGV) for growth rate in which the highest latitudinal population exhibited greater growth rates than their counterparts from lower latitudes. Surprisingly, this CnGV pattern was accompanied by a lower apparent dry-matter digestibility, which could highlight a differential assimilation of ingested nutrients into somatic tissue, revealing enhanced growth efficiency in snails from the highest latitudinal habitat. Our evidence highlights that adjustments in energy allocation to the digestive machinery and to protein storage could act as a latitudinal compensation for enhanced growth efficiency in snails from the highest latitudinal population.

Hybridization at an ecotone: ecological and genetic barriers between three Iberian vipers

The formation of stable genetic boundaries between emerging species is often diagnosed by reduced hybrid fitness relative to parental taxa. This reduced fitness can arise from endogenous and/or exogenous barriers to gene flow. Although detecting exogenous barriers in nature is difficult, we can estimate the role of ecological divergence in driving species boundaries by integrating molecular and ecological niche modelling tools. Here, we focus on a three-way secondary contact zone between three viper species (Vipera aspis, V. latastei and V. seoanei) to test for the contribution of ecological divergence to the development of reproductive barriers at several species traits (morphology, nuclear DNA and mitochondrial DNA). Both the nuclear and mitochondrial data show that all taxa are genetically distinct and that the sister species V. aspis and V. latastei hybridize frequently and backcross over several generations. We find that the three taxa have diverged ecologically and meet at a hybrid zone coincident with a steep ecotone between the Atlantic and Mediterranean biogeographical provinces. Integrating landscape and genetic approaches, we show that hybridization is spatially restricted to habitats that are suboptimal for parental taxa. Together, these results suggest that niche separation and adaptation to an ecological gradient confer an important barrier to gene flow among taxa that have not achieved complete reproductive isolation.

Adaptations to "Thermal Time" Constraints in Papilio: Latitudinal and Local Size Clines Differ in Response to Regional Climate Change

Adaptations to "thermal time" (=Degree-day) constraints on developmental rates and voltinism for North American tiger swallowtail butterflies involve most life stages, and at higher latitudes include: smaller pupae/adults; larger eggs; oviposition on most nutritious larval host plants; earlier spring adult emergences; faster larval growth and shorter molting durations at lower temperatures. Here we report on forewing sizes through 30 years for both the northern univoltine P. canadensis (with obligate diapause) from the Great Lakes historical hybrid zone northward to central Alaska (65° N latitude), and the multivoltine, P. glaucus from this hybrid zone southward to central Florida (27° N latitude). Despite recent climate warming, no increases in mean forewing lengths of P. glaucus were observed at any major collection location (FL to MI) from the 1980s to 2013 across this long latitudinal transect (which reflects the "converse of Bergmann's size Rule", with smaller females at higher latitudes). Unlike lower latitudes, the Alaska, Ontonogon, and Chippewa/Mackinac locations (for P. canadensis) showed no significant increases in D-day accumulations, which could explain lack of size change in these northernmost locations. As a result of 3-4 decades of empirical data from major collection sites across these latitudinal clines of North America, a general "voltinism/size/D-day" model is presented, which more closely predicts female size based on D-day accumulations, than does latitude. However, local "climatic cold pockets" in northern Michigan and Wisconsin historically appeared to exert especially strong size constraints on female forewing lengths, but forewing lengths quickly increased with local summer warming during the recent decade, especially near the warming edges of the cold pockets. Results of fine-scale analyses of these "cold pockets" are in contrast to non-significant changes for other Papilio populations seen across the latitudinal transect for P. glaucus and P. canadensis in general, highlighting the importance of scale in adaptations to climate change. Furthermore, we also show that rapid size increases in cold pocket P. canadensis females with recent summer warming are more likely to result from phenotypic plasticity than genotypic introgression from P. glaucus, which does increase size in late-flight hybrids and P. appalachiensis.

Climate-Driven Reshuffling of Species and Genes: Potential Conservation Roles for Species Translocations and Recombinant Hybrid Genotypes

Comprising 50%-75% of the world's fauna, insects are a prominent part of biodiversity in communities and ecosystems globally. Biodiversity across all levels of biological classifications is fundamentally based on genetic diversity. However, the integration of genomics and phylogenetics into conservation management may not be as rapid as climate change. The genetics of hybrid introgression as a source of novel variation for ecological divergence and evolutionary speciation (and resilience) may generate adaptive potential and diversity fast enough to respond to locally-altered environmental conditions. Major plant and herbivore hybrid zones with associated communities deserve conservation consideration. This review addresses functional genetics across multi-trophic-level interactions including "invasive species" in various ecosystems as they may become disrupted in different ways by rapid climate change. "Invasive genes" (into new species and populations) need to be recognized for their positive creative potential and addressed in conservation programs. "Genetic rescue" via hybrid translocations may provide needed adaptive flexibility for rapid adaptation to environmental change. While concerns persist for some conservationists, this review emphasizes the positive aspects of hybrids and hybridization. Specific implications of natural genetic introgression are addressed with a few examples from butterflies, including transgressive phenotypes and climate-driven homoploid recombinant hybrid speciation. Some specific examples illustrate these points using the swallowtail butterflies (Papilionidae) with their long-term historical data base (phylogeographical diversity changes) and recent (3-decade) climate-driven temporal and genetic divergence in recombinant homoploid hybrids and relatively recent hybrid speciation of Papilio appalachiensis in North America. Climate-induced "reshuffling" (recombinations) of species composition, genotypes, and genomes may become increasingly ecologically and evolutionarily predictable, but future conservation management programs are more likely to remain constrained by human behavior than by lack of academic knowledge.