Gis-based niche models reveal unifying climatic mechanisms that maintain the location of avian hybrid zones in a North American suture zone

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.

Local niches explain coexistence in environmentally-distinct contact zones between Western Mediterranean vipers

Species' ecological niches are frequently analysed to gain insights into how anthropogenic changes affect biodiversity. Coping with these changes often involves shifts in niche expression, which can disrupt local biotic interactions. Secondary contact zones, where competition and ecological segregation commonly occur, are ideal for studying the ecological factors influencing species' niches. In this study, we investigated the effect of climate and landscape factors on the ecological niches of two viper species, Vipera aspis and Vipera latastei, across three contact zones in northern Iberia, characterized by varying levels of landscape alteration. Using niche overlap tests, ecological niche models and spatial analyses we observed local variation in the expression of the species' niches across the three contact zones, resulting from the different abiotic and biotic conditions of each area. Rather than spatial niche segregation, we observed high niche overlap, suggesting niche convergence. A pattern of asymmetrical niche variation was identified in all contact zones, driven by species' climatic tolerances and the environmental conditions of each area. V. aspis generally exhibited a wider niche, except in the southernmost zone where the pure Mediterranean climate favored V. latastei. Human-induced landscape changes intensified niche asymmetry, by favoring the most generalist V. aspis over the specialist V. latastei, increasing habitat overlap, and likely competition. This study presents a comprehensive analysis of niche expression at range margins, anticipating a heightened impact of landscape changes in V. latastei. The methodological framework implemented here, and our findings, hold significant relevance for biodiversity management and conservation in human-impacted areas.

Can demographic histories explain long-term isolation and recent pulses of asymmetric gene flow between highly divergent grey fox lineages?

Secondary contact zones between deeply divergent, yet interfertile, lineages provide windows into the speciation process. North American grey foxes (Urocyon cinereoargenteus) are divided into western and eastern lineages that diverged approximately 1 million years ago. These ancient lineages currently hybridize in a relatively narrow zone of contact in the southern Great Plains, a pattern more commonly observed in smaller-bodied taxa, which suggests relatively recent contact after a long period of allopatry. Based on local ancestry inference with whole-genome sequencing (n = 43), we identified two distinct Holocene pulses of admixture. The older pulse (500-3500 YBP) reflected unidirectional gene flow from east to west, whereas the more recent pulse (70-200 YBP) of admixture was bi-directional. Augmented with genotyping-by-sequencing data from 216 additional foxes, demographic analyses indicated that the eastern lineage declined precipitously after divergence, remaining small throughout most of the late Pleistocene, and expanding only during the Holocene. Genetic diversity in the eastern lineage was highest in the southeast and lowest near the contact zone, consistent with a westward expansion. Concordantly, distribution modelling indicated that during their isolation, the most suitable habitat occurred far east of today's contact zone or west of the Great Plains. Thus, long-term isolation was likely caused by the small, distant location of the eastern refugium, with recent contact reflecting a large increase in suitable habitat and corresponding demographic expansion from the eastern refugium. Ultimately, long-term isolation in grey foxes may reflect their specialized bio-climatic niche. This system presents an opportunity for future investigation of potential pre- and post-zygotic isolating mechanisms.

Lower survival of hybrid grosbeaks, but not towhees, suggests a molt divide disfavors hybrids

Although avian hybrid zones in the Great Plains have been studied for almost 70 years, we know surprisingly little about the fitness costs to hybrids that keep these zones narrow. We compare age ratios in grosbeaks (Pheucticus ludovicianus and P. melanocephalus) and towhees (Pipilo erythropthalums and P. maculatus), two species pairs that differ in their life histories and molt schedules, to evaluate survival between hybrids and parentals. We then contrast molt and migratory divides as possible sources of selection against hybrids. Hybrid grosbeaks had 27%-33% lower survival relative to their parentals, whereas hybrid towhees had survival rates similar to parentals. Age ratio data for hybrid grosbeaks suggest high mortality in older birds, as expected if selection operates after the first year of life. This pattern is consistent with parental species of grosbeaks having contrasting molt schedules relative to migration, suggesting high mortality costs to hybrids driven by molt biology, which are expressed later in life. Contrasts in molt schedules are absent in towhees. While migratory divides may exist for towhees and grosbeaks, the low adult survival of hybrid grosbeaks suggest that molt may be an important and underappreciated source of selection maintaining this and other narrow avian hybrid zones.

Complex cycles of divergence and migration shape lineage structure in the common kingsnake species complex

Aim

The Nearctic is a complex patchwork of habitats and geologic features that form barriers to gene flow resulting in phylogeographic structure and speciation in many lineages. Habitats are rarely stable over geologic time, and the Nearctic has undergone major climatic changes in the past few million years. We use the common kingsnake species complex to study how climate, geography, and history influence lineage formation over a large, complex landscape.

Location

Nearctic/North America.

Taxon

Common kingsnake, Lampropeltis getula, species complex.

Methods

We analyzed genome-wide sequence data from 51 snakes spanning the majority of the species complex's range. We used population clustering, generalized dissimilarity modeling, and coalescent methods to identify the number of genetic clusters within the L. getula complex, infer the environmental correlates of genetic differentiation, and estimate models of divergence and gene flow among lineages.

Results

We identified three major lineages within the L. getula complex and further continuous spatial structure within lineages. The most important ecological correlates of genetic distance in the complex are related to aridity and precipitation, consistent with lineage breaks at the Great Plains/Desert ecotone and the Cochise Filter Barrier. Lineages are estimated to have undergone multiple rounds of isolation and secondary contact, with highly asymmetric migration occurring at present.

Main conclusions

Changing climates combined with a large and geologically complex landscape have resulted in a mosaic of discrete and spatially continuous genetic structure. Multiple rounds of isolation and secondary contact as climate fluctuated over the past ~4.4 My have likely driven the evolution of discrete lineages that maintain high levels of gene flow. Continuous structure is strongly shaped by aridity and precipitation, suggesting roles for major precipitation gradients in helping to maintain lineage identity in the face of gene flow when lineages are in geographic contact.

Bird Communities in a Changing World: The Role of Interspecific Competition

Significant changes in the environment have the potential to affect bird species abundance and distribution, both directly, through a modification of the landscape, habitats, and climate, and indirectly, through a modification of biotic interactions such as competitive interactions. Predicting and mitigating the consequences of global change thus requires not only a sound understanding of the role played by biotic interactions in current ecosystems, but also the recognition and study of the complex and intricate effects that result from the perturbation of these ecosystems. In this review, we emphasize the role of interspecific competition in bird communities by focusing on three main predictions derived from theoretical and empirical considerations. We provide numerous examples of population decline and displacement that appeared to be, at least in part, driven by competition, and were amplified by environmental changes associated with human activities. Beyond a shift in relative species abundance, we show that interspecific competition may have a negative impact on species richness, ecosystem services, and endangered species. Despite these findings, we argue that, in general, the role played by interspecific competition in current communities remains poorly understood due to methodological issues and the complexity of natural communities. Predicting the consequences of global change in these communities is further complicated by uncertainty regarding future environmental conditions and the speed and efficacy of plastic and evolutionary responses to fast-changing environments. Possible directions of future research are highlighted.

Revisiting a classic hybrid zone: Movement of the northern flicker hybrid zone in contemporary times

Natural hybrid zones have provided important insights into the evolutionary process, and their geographic dynamics over time can help to disentangle the underlying biological processes that maintain them. Here, we leverage replicated sampling of an identical transect across the hybrid zone between yellow-shafted and red-shafted flickers in the Great Plains to assess its stability over ∼60 years (1955-1957 to 2016-2018). We identify a ∼73-km westward shift in the hybrid zone center toward the range of the red-shafted flicker, but find no associated changes in width over our sampling period. In fact, the hybrid zone remains remarkably narrow, suggesting some kind of selective pressure maintains the zone. By comparing to previous work in the same geographic region, it appears likely that the movement in the hybrid zone has occurred in the years since the early 1980s. This recent movement may be related to changes in climate or land management practices that have allowed westward movement of yellow-shafted flickers into the Great Plains.

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.

Ecological Divergence and the History of Gene Flow in the Nearctic Milksnakes (Lampropeltis triangulum Complex)

Many phylogeographic studies on species with large ranges have found genetic-geographic structure associated with changes in habitat and physical barriers preventing or reducing gene flow. These interactions with geographic space, contemporary and historical climate, and biogeographic barriers have complex effects on contemporary population genetic structure and processes of speciation. While allopatric speciation at biogeographic barriers is considered the primary mechanism for generating species, more recently it has been shown that parapatric modes of divergence may be equally or even more common. With genomic data and better modeling capabilities, we can more clearly define causes of speciation in relation to biogeography and migration between lineages, the location of hybrid zones with respect to the ecology of parental lineages, and differential introgression of genes between taxa. Here, we examine the origins of three Nearctic milksnakes (Lampropeltis elapsoides, Lampropeltis triangulum and Lampropeltis gentilis) using genome-scale data to better understand species diversification. Results from artificial neural networks show that a mix of a strong biogeographic barrier, environmental changes, and physical space has affected genetic structure in these taxa. These results underscore conspicuous environmental changes that occur as the sister taxa L. triangulum and L. gentilis diverged near the Great Plains into the forested regions of the Eastern Nearctic. This area has been recognized as a region for turnover for many vertebrate species, but as we show here the contemporary boundary does not isolate these sister species. These two species likely formed in the mid-Pleistocene and have remained partially reproductively isolated over much of this time, showing differential introgression of loci. We also demonstrate that when L. triangulum and L. gentilis are each in contact with the much older L. elapsoides, some limited gene flow has occurred. Given the strong agreement between nuclear and mtDNA genomes, along with estimates of ecological niche, we suggest that all three lineages should continue to be recognized as unique species. Furthermore, this work emphasizes the importance of considering complex modes of divergence and differential allelic introgression over a complex landscape when testing mechanisms of speciation. [Cline; delimitation; Eastern Nearctic; Great Plains; hybrids; introgression; speciation.].

Estimating hybridization in the wild using citizen science data: a path forward

Genomic evidence of introgression in natural populations has reinvigorated the study of hybridization in recent years. Still, it is largely unknown how frequently individual organisms mate across species lines. Recently, Justyn et al. 2020 suggested that eBird, one of the world's largest citizen science databases, may supply adequate data for estimating hybridization rates. Here, we compare Justyn et al.'s estimates-and their conclusions that hybridization is rare-with estimates from museum and molecular data. We also estimate hybridization using eBird observations from areas and times when hybridization is possible, namely, in contact zones during the breeding season. These estimates are all considerably higher than those reported in Justyn et al., emphasizing that inferences from multiple datasets can differ radically. Finally, we demonstrate an approach for predicting the location of hybrid zones using eBird data, which can be done with high confidence and with unprecedented resolution. We show that citizen science data, far from settling the question of how frequently bird species hybridize, instead offer a promising step toward more focused study of hybrid zones. This article is protected by copyright. All rights reserved.

Genetic, bioacoustic and morphological analyses reveal cryptic speciation in the warbling vireo complex (Vireo gilvus: Vireonidae: Passeriformes)

Cryptic species are closely related taxa that are difficult to separate morphologically, but are reproductively isolated. Here we examine the warbling vireo complex (Vireo gilvus), a widespread songbird speculated to be comprised of more than one cryptic species. We included three taxa within the complex: two of the western (Vireo gilvus swainsonii and Vireo gilvus brewsteri) subspecies and the single eastern (Vireo gilvus gilvus) subspecies. We used mtDNA and microsatellite loci to assess the congruence of genetic data to the current subspecies boundaries. We then incorporated bioacoustic, morphometric and ecological niche modelling analyses to further examine differences. We found two genetic groups with mtDNA analysis, splitting eastern and western warbling vireos. Microsatellite analyses revealed four genetic groups: an eastern group, a Black Hills group and two western groups that do not agree with current western subspecies boundaries based on phenotypic data. Our results suggest that eastern and western warbling vireos have been reproductively isolated for a long period of time and therefore may be best treated as separate species. However, more research into areas of contact to examine the presence of hybridization is advised before making a taxonomic revision. Differences between the two western genetic groups appear less clear, requiring additional research.

Mitochondrial Genomes of the United States Distribution of Gray Fox (Urocyon cinereoargenteus) Reveal a Major Phylogeographic Break at the Great Plains Suture Zone

We examined phylogeographic structure in gray fox (Urocyon cinereoargenteus) across the United States to identify the location of secondary contact zone(s) between eastern and western lineages and investigate the possibility of additional cryptic intraspecific divergences. We generated and analyzed complete mitochondrial genome sequence data from 75 samples and partial control region mitochondrial DNA sequences from 378 samples to investigate levels of genetic diversity and structure through population- and individual-based analyses including estimates of divergence (FST and SAMOVA), median joining networks, and phylogenies. We used complete mitochondrial genomes to infer phylogenetic relationships and date divergence times of major lineages of Urocyon in the United States. Despite broad-scale sampling, we did not recover additional major lineages of Urocyon within the United States, but identified a deep east-west split (∼0.8 million years) with secondary contact at the Great Plains Suture Zone and confirmed the Channel Island fox (Urocyon littoralis) is nested within U. cinereoargenteus. Genetic diversity declined at northern latitudes in the eastern United States, a pattern concordant with post-glacial recolonization and range expansion. Beyond the east-west divergence, morphologically-based subspecies did not form monophyletic groups, though unique haplotypes were often geographically limited. Gray foxes in the United States displayed a deep, cryptic divergence suggesting taxonomic revision is needed. Secondary contact at a common phylogeographic break, the Great Plains Suture Zone, where environmental variables show a sharp cline, suggests ongoing evolutionary processes may reinforce this divergence. Follow-up study with nuclear markers should investigate whether hybridization is occurring along the suture zone and characterize contemporary population structure to help identify conservation units. Comparative work on other wide-ranging carnivores in the region should test whether similar evolutionary patterns and processes are occurring.

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.

Phylogeography and the Evolutionary History of Sunflower (Helianthus annuus L.): Wild Diversity and the Dynamics of Domestication

Patterns of genetic variation in crops are the result of selection and demographic changes that occurred during their domestication and improvement. In many cases, we have an incomplete picture of the origin of crops in the context of their wild progenitors, particularly with regard to the processes producing observed levels of standing genetic variation. Here, we analyzed sequence diversity in cultivated sunflower (Helianthus annuus L.) and its wild progenitor (common sunflower, also H. annuus) to reconstruct phylogeographic relationships and population genetic/demographic patterns across sunflower. In common sunflower, south-north patterns in the distribution of nucleotide diversity and lineage splitting indicate a history of rapid postglacial range expansion from southern refugia. Cultivated sunflower accessions formed a clade, nested among wild populations from the Great Plains, confirming a single domestication event in central North America. Furthermore, cultivated accessions sorted by market type (i.e., oilseed vs. confectionery) rather than breeding pool, recapitulating the secondary development of oil-rich cultivars during its breeding history. Across sunflower, estimates of nucleotide diversity and effective population sizes suggest that cultivated sunflower underwent significant population bottlenecks following its establishment ~5000 years ago. The patterns inferred here corroborate those from previous studies of sunflower domestication, and provide a comprehensive overview of its evolutionary history.

Population genomics of two congeneric Palaearctic shorebirds reveals differential impacts of Quaternary climate oscillations across habitats types

Intracontinental biotic divisions across the vast Palaearctic region are not well-characterized. Past research has revealed patterns ranging from a lack of population structure to deep divergences along varied lines of separation. Here we compared biogeographic patterns of two Palaearctic shorebirds with different habitat preferences, Whimbrel (Numenius phaeopus) and Eurasian curlew (N. arquata). Using genome-wide markers from populations across the Palaearctic, we applied a multitude of population genomic and phylogenomic approaches to elucidate population structure. Most importantly, we tested for isolation by distance and visualized barriers and corridors to gene flow. We found shallow Palaearctic population structure in subpolar bog and tundra-breeding whimbrels, consistent with other species breeding at a similarly high latitude, indicating connectivity across the tundra belt, both presently and during southward shifts in periods of global cooling. In contrast, the temperate grassland-breeding Eurasian curlew emerged in three distinct clades corresponding to glacial refugia. Barriers to gene flow coincided with areas of topographic relief in the central Palaearctic for whimbrels and further east for Eurasian curlews. Our findings highlight the interplay of historic and ecological factors in influencing present-day population structure of Palaearctic biota.

Modelling the potential effects of climate factors on Dubas bug (Ommatissus lybicus) presence/absence and its infestation rate: A case study from Oman

BACKGROUND

Ommatissus lybicus de Bergevin (Hemiptera: Tropiduchidae) (Dubas Bug, DB) is an insect pest attacking date palms. It occurs in Arab countries including Oman. In this paper, the logistic, ordinary least square, and geographical weighted regressions were applied to model the absence/presence and density of DB against climate factors. A method is proposed for modelling spatially correlated prorations annually over the study period, based on annual and seasonal outbreaks. The historical 2006-2015 climate data were obtained from weather stations located in nine governorates in northern Oman, while dataloggers collected the 2017 microclimate data in eight of these nine governorates.

RESULTS

Logistic regression model showed the percentages of correctly predicted values using a cut-off point of 0.5 were 90%, 88% and 84%, indicating good classification accuracy. OLS and GWR models showed an overall trend of strong linear correlation between DB infestation levels and short- and long-term climate factors. The three models suggested that precipitation, elevation, temperature, humidity, wind direction and wind speed are important in influencing the spatial distribution and the presence/absence of dense DB populations.

CONCLUSION

The results provide an improved understanding of climate factors that impact DB's spread and is considered useful for managing DB infestations in date palm plantations. © 2019 Society of Chemical Industry.

Phylogeography and niche modelling: reciprocal enlightenment

Phylogeography examines the spatial genetic structure of species. Environmental niche modelling (or ecological niche modelling; ENM) examines the environmental limits of a species’ ecological niche. These two fields have great potential to be used together. ENM can shed light on how phylogeographical patterns develop and help identify possible drivers of spatial structure that need to be further investigated. Specifically, ENM can be used to test for niche differentiation among clades, identify factors limiting individual clades and identify barriers and contact zones. It can also be used to test hypotheses regarding the effects of historical and future climate change on spatial genetic patterns by projecting niches using palaeoclimate or future climate data. Conversely, phylogeographical information can populate ENM with within-species genetic diversity. Where adaptive variation exists among clades within a species, modelling their niches separately can improve predictions of historical distribution patterns and future responses to climate change. Awareness of patterns of genetic diversity in niche modelling can also alert conservationists to the potential loss of genetically diverse areas in a species’ range. Here, we provide a simplistic overview of both fields, and focus on their potential for integration, encouraging researchers on both sides to take advantage of the opportunities available.

Genomic divergence in allopatric Northern Cardinals of the North American warm deserts is linked to behavioral differentiation

Biogeographic barriers are considered important in initiating speciation through geographic isolation, but they rarely indiscriminately and completely reduce gene flow across entire communities. Explicitly demonstrating which factors are associated with gene-flow levels across barriers would help elucidate how speciation is initiated and isolation maintained. Here, we investigated the association of behavioral isolation on population differentiation in Northern Cardinals (Cardinalis cardinalis) distributed across the Cochise Filter Barrier, a region of transitional habitat which separates the Sonoran and Chihuahuan deserts of North America. Using genomewide markers, we modeled demographic history by fitting the data to isolation and isolation-with-migration models. The best-fit model indicated that desert populations diverged in the Pleistocene with low, historic, and asymmetric gene flow across the barrier. We then tested behavioral isolation using reciprocal call-broadcast experiments to compare song recognition between deserts, controlling for song dialect changes within deserts. We found that male Northern Cardinals in both deserts were most aggressive to local songs and failed to recognize across-barrier songs. A correlation of genomic differentiation and strong song discrimination is consistent with a model where speciation is initiated across a barrier and maintained by behavioral isolation.

Species Distributions, Quantum Theory, and the Enhancement of Biodiversity Measures

Species distributions are typically represented by records of their observed occurrence at a given spatial and temporal scale. Such records are inevitably incomplete and contingent on the spatial-temporal circumstances under which the observations were made. Moreover, organisms may respond differently to similar environmental conditions at different places or moments, so their distribution is, in principle, not completely predictable. We argue that this uncertainty exists, and warrants considering species distributions as analogous to coherent quantum objects, whose distributions are better described by a wavefunction rather than by a set of locations. We use this to extend the existing concept of "dark diversity", which incorporates into biodiversity metrics those species that could, but which have not yet been observed to, inhabit a region-thereby developing the idea of "potential biodiversity". We show how conceptualizing species' distributions in this way could help overcome important weaknesses in current biodiversity metrics, both in theory and by using a worked case study of mammal distributions in Spain over the last decade. We propose that considerable theoretical advances could eventually be gained through interdisciplinary collaboration between biogeographers and quantum physicists. [Biogeography; favorability; physics; predictability; probability; species occurrence; uncertainty; wavefunction.

Origin and cross-century dynamics of an avian hybrid zone

Background

Characterizations of the dynamics of hybrid zones in space and time can give insights about traits and processes important in population divergence and speciation. We characterized a hybrid zone between tanagers in the genus Ramphocelus (Aves, Thraupidae) located in southwestern Colombia. We evaluated whether this hybrid zone originated as a result of secondary contact or of primary differentiation, and described its dynamics across time using spatial analyses of molecular, morphological, and coloration data in combination with paleodistribution modeling.

Results

Models of potential historical distributions based on climatic data and genetic signatures of demographic expansion suggested that the hybrid zone likely originated following secondary contact between populations that expanded their ranges out of isolated areas in the Quaternary. Concordant patterns of variation in phenotypic characters across the hybrid zone and its narrow extent are suggestive of a tension zone, maintained by a balance between dispersal and selection against hybrids. Estimates of phenotypic cline parameters obtained using specimens collected over nearly a century revealed that, in recent decades, the zone appears to have moved to the east and to higher elevations, and may have become narrower. Genetic variation was not clearly structured along the hybrid zone, but comparisons between historical and contemporary specimens suggested that temporal changes in its genetic makeup may also have occurred.

Conclusions

Our data suggest that the hybrid zone likey resulted from secondary contact between populations. The observed changes in the hybrid zone may be a result of sexual selection, asymmetric gene flow, or environmental change.

Electronic supplementary material

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

Environmental (in)dependence of a hybrid zone: Insights from molecular markers and ecological niche modeling in a hybrid zone of Origanum (Lamiaceae) on the island of Crete

The role of environment and the relative significance of endogenous versus exogenous selection in shaping hybrid zones have been crucial issues in the studies of hybridization. Recent advances in ecological niche modeling (ENM) offer new methodological tools, especially in combination with the genotyping of individuals in the hybrid zone. Here, we study the hybrid zone between the widely known spices Origanum onites and Origanum vulgare ssp. hirtum in Crete. We analyze the genetic structure of both parental taxa and their hybrid Origanum × intercendens using AFLP markers on 15 sympatric and 12 allopatric populations and employ ecological niche modeling and niche similarity tests to study their niche patterns. We complement these analyses with seed viability measurements. Our study revealed that the hybridizing taxa O. onites and O. vulgare ssp. hirtum and the resulting genotypic classes showed geographical and environmental niche similarities based on the predictions of ENMs and the subsequent similarity tests. The occurrence of the hybrid zone is not directly dependent on environmental factors which favor the fitness of the hybrid compared to the parental taxa, but rather on aspects such as historical factors and management practices, which may contribute to the localization and maintenance of the contact zone between parental species. Our results suggest that if a minimum required niche differentiation between genotypic classes is not achieved, environmental dependence might not have a prominent role on the outcome of the hybridization.

Changing climate mediates sapsucker (Aves: Sphyrapicus) hybrid zone movement

Hybrid zones, where two divergent taxa meet and interbreed, offer unique opportunities to investigate how climate contributes to reproductive isolation between closely related taxa and how these taxa may respond to climatic changes. Red-naped (Sphyrapicus nuchalis) and Red-breasted (Sphyrapicus ruber) sapsuckers (Aves: Picidae) hybridize along a narrow contact zone that stretches from northern California to British Columbia. The hybrid zone between these species has been studied extensively for more than 100 years and represents an excellent system for investigations of the evolution of reproductive isolation. Shifts in the proportions of phenotypes at hybrid localities since 1910 that were inferred using specimens from museum collections were confirmed using species distribution models. We predicted the historical, current, and future distributions of parental and hybrid sapsuckers using Random Forests models to quantify how climate change is affecting hybrid zone movement in the Pacific Northwest. We found observed distribution shifts of parental sapsuckers were likely the result of climate change over the past 100 years, with these shifts predicted to continue for both sapsuckers over the next 80 years. We found Red-breasted Sapsuckers are predicted to continue to expand, while Red-naped Sapsuckers are predicted to contract substantially under future climate scenarios. As a result of the predicted changes, the amount of overlap in the distribution of these sapsuckers may decrease. Using hybrid phenotypes, we found the climate niche occupied by the hybrid zone is predicted to disappear under future conditions. The disappearance of this climate niche where the two parental species come into contact and hybridize may lead to a substantial reduction in genetic introgression. Understanding the impacts of global climate change on hybrid zones may help us to better understand how speciation has been shaped by climate in the past, as well as how evolution may respond to climate change in the future.

Variation in complex mating signals in an "island" hybrid zone between Stenobothrus grasshopper species

Two grasshopper species Stenobothrus rubicundus and S. clavatus were previously shown to meet in a narrow hybrid zone on Mount Tomaros in northern Greece. The species are remarkable for their complex courtship songs accompanied by conspicuous movements of antennae and wings. We analyzed variations in forewing morphology, antenna shape, and courtship song across the hybrid zone using a geographic information system, and we documented three contact zones on Mount Tomaros. All male traits and female wings show abrupt transitions across the contact zones, suggesting that these traits are driven by selection rather than by drift. Male clines in antennae are displaced toward S. clavatus, whereas all clines in wings are displaced toward S. rubicundus. We explain cline discordance as depending on sexual selection via female choice. The high covariance between wings and antennae found in the centers of all contact zones results from high levels of linkage disequilibria among the underlying loci, which in turn more likely results from assortative mating than from selection against hybrids. The covariance is found to be higher in clavatus-like than rubicundus-like populations, which implies asymmetric assortative mating in parental-like sites of the hybrid zone and a movement of the hybrid zone in favor of S. clavatus.

Potential Implications of Climate Change on Aegilops Species Distribution: Sympatry of These Crop Wild Relatives with the Major European Crop Triticum aestivum and Conservation Issues

Gene flow from crop to wild relatives is a common phenomenon which can lead to reduced adaptation of the wild relatives to natural ecosystems and/or increased adaptation to agrosystems (weediness). With global warming, wild relative distributions will likely change, thus modifying the width and/or location of co-occurrence zones where crop-wild hybridization events could occur (sympatry). This study investigates current and 2050 projected changes in sympatry levels between cultivated wheat and six of the most common Aegilops species in Europe. Projections were generated using MaxEnt on presence-only data, bioclimatic variables, and considering two migration hypotheses and two 2050 climate scenarios (RCP4.5 and RCP8.5). Overall, a general decline in suitable climatic conditions for Aegilops species outside the European zone and a parallel increase in Europe were predicted. If no migration could occur, the decline was predicted to be more acute outside than within the European zone. The potential sympatry level in Europe by 2050 was predicted to increase at a higher rate than species richness, and most expansions were predicted to occur in three countries, which are currently among the top four wheat producers in Europe: Russia, France and Ukraine. The results are also discussed with regard to conservation issues of these crop wild relatives.

Differential effects of climate and species interactions on range limits at a hybrid zone: potential direct and indirect impacts of climate change

The relative contributions of climate versus interspecific interactions in shaping species distributions have important implications for closely related species at contact zones. When hybridization occurs within a contact zone, these factors regulate hybrid zone location and movement. While a hybrid zone's position may depend on both climate and interactions between the hybridizing species, little is known about how these factors interact to affect hybrid zone dynamics. Here, we utilize SDM (species distribution modeling) both to characterize the factors affecting the current location of a moving North American avian hybrid zone and to predict potential direct and indirect effects of climate change on future distributions. We focus on two passerine species that hybridize where their ranges meet, the Black‐capped (Poecile atricapillus) and Carolina (P. carolinensis) chickadee. Our contemporary climate models predict the occurrence of climatically suitable habitat extending beyond the hybrid zone for P. atricapillus only, suggesting that interspecific interactions primarily regulate this range boundary in P. atricapillus, while climatic factors regulate P. carolinensis. Year 2050 climate models predict a drastic northward shift in suitable habitat for P. carolinensis. Because of the greater importance of interspecific interactions for regulating the southern range limit of P. atricapillus, these climate‐mediated shifts in the distribution of P. carolinensis may indirectly lead to a range retraction in P. atricapillus. Together, our results highlight the ways climate change can both directly and indirectly affect species distributions and hybrid zone location. In addition, our study lends support to the longstanding hypothesis that abiotic factors regulate species' poleward range limits, while biotic factors shape equatorial range limits.

Environmental factors contribute to the formation and maintenance of the contact zone observed in deciduous broad-leaved tree species in Japan

Contact zones are defined as areas where populations from different refugia meet during a postglacial expansion and distinct DNA lineages are mixedly distributed. In Japan, contact zones of various plants and animals were reported from the Kinki-Chugoku region. These contact zones appear to be maintained without any drastic topographic barriers such as those observed in the Alps and Pyrenees Mountains. In this study, the mechanisms underlying the formation and/or maintenance of these contact zones were investigated using six deciduous broad-leaved tree species (Carpinus laxiflora, C. tschonoskii, C. japonica, Magnolia obovata, Padus grayana, and Euonymus oxyphyllus). First, the precise location of the contact zones was examined by intensive genetic analysis of the six species. Second, the relationships between the geographic location of the contact zone and various environmental factors, including climate and topography, were investigated by generalized additive models to reveal the mechanisms of the formation and maintenance of the contact zones. As a result, four of the six examined plant species clearly showed a geographically common contact zone in Hyogo Prefecture and its adjacent areas. The results of the generalized additive models indicate that the pattern of low habitat suitability estimated by ecological niche modeling was the most important factor for determining the location of the common contact zone. These results suggest that areas with low habitat suitability in Hyogo Prefecture restrict the migration and gene flow of the four species in this region, and thus, they maintain the pattern of the contact zones. This study suggests that there are major effects of habitat suitability on the formation and maintenance of the contact zones.

Munroa argentina, a Grass of the South American Transition Zone, Survived the Andean Uplift, Aridification and Glaciations of the Quaternary

The South American Transition Zone (SATZ) is a biogeographic area in which not only orogeny (Andes uplift) and climate events (aridification) since the mid-Miocene, but also Quaternary glaciation cycles had an important impact on the evolutionary history of the local flora. To study this effect, we selected Munroa argentina, an annual grass distributed in the biogeographic provinces of Puna, Prepuna and Monte. We collected 152 individuals from 20 localities throughout the species' range, ran genetic and demographic analyses, and applied ecological niche modeling. Phylogenetic and population genetic analyses based on cpDNA and AFLP data identified three phylogroups that correspond to the previously identified subregions within the SATZ. Molecular dating suggests that M. argentina has inhabited the SATZ since approximately 3.4 (4.2-1.2) Ma and paleomodels predict suitable climate in these areas during the Interglacial period and the Last Glacial Maximum. We conclude that the current distribution of M. argentina resulted from the fragmentation of its once continuous range and that climate oscillations promoted ecological differences that favored isolation by creating habitat discontinuity.

Contrasting the effects of environment, dispersal and biotic interactions to explain the distribution of invasive plants in alpine communities

Despite considerable efforts devoted to investigate the community assembly processes driving plant invasions, few general conclusions have been drawn so far. Three main processes, generally acting as successive filters, are thought to be of prime importance. The invader has to disperse (1st filter) into a suitable environment (2nd filter) and succeed in establishing in recipient communities through competitive interactions (3rd filter) using two strategies: competition avoidance by the use of different resources (resource opportunity), or competitive exclusion of native species. Surprisingly, despite the general consensus on the importance of investigating these three processes and their interplay, they are usually studied independently. Here we aim to analyse these three filters together, by including them all: abiotic environment, dispersal and biotic interactions, into models of invasive species distributions. We first propose a suite of indices (based on species functional dissimilarities) supposed to reflect the two competitive strategies (resource opportunity and competition exclusion). Then, we use a set of generalised linear models to explain the distribution of seven herbaceous invaders in natural communities (using a large vegetation database for the French Alps containing 5,000 community-plots). Finally, we measure the relative importance of competitive interaction indices, identify the type of coexistence mechanism involved and how this varies along environmental gradients. Adding competition indices significantly improved model's performance, but neither resource opportunity nor competitive exclusion were common strategies among the seven species. Overall, we show that combining environmental, dispersal and biotic information to model invasions has excellent potential for improving our understanding of invader success.

Niche-habitat mechanisms and biotic interactions explain the coexistence and abundance of congeneric sandgrouse species

Ascertaining which niche processes allow coexistence between closely related species is of special interest in ecology. We quantified variations in the environmental niches and densities of two congeneric species, the pin-tailed and the black-bellied sandgrouse (Pterocles alchata and Pterocles orientalis) in allopatry and sympatry under similar abiotic, habitat and dispersal contexts to understand their coexistence. Using principal component analysis, we defined environmental gradients (niche dimensions) including abiotic, habitat and anthropogenic variables, and calculated niche breadth, position and overlap of both species in sympatry and allopatry. Additionally, sandgrouse density was modelled as a function of the niche dimensions and the density of the other species. We found evidence that each species occupies distinct environmental niches in sympatry and in allopatry. The black-bellied sandgrouse exploits a broader range of environmental conditions (wider niche breadth) while the pin-tailed sandgrouse reaches high densities where conditions seem to match its optimum. In sympatry, both species shift their niches to intermediate positions, indicating the importance of abiotic factors in setting coexistence areas. Environmental conditions determine regional densities of pin-tailed sandgrouse whereas biotic interactions explain the density of the black-bellied sandgrouse in areas with abiotic conditions similarly conducive for both species. Highly suitable areas for the pin-tailed sandgrouse fall beyond the upper thermal limit of the black-bellied sandgrouse, leading to niche segregation and low densities for the latter. Finally, local niche shift and expansion plus possible heterospecific aggregation allow the pin-tailed sandgrouse to thrive in a priori less favourable environments. This work provides insight into how different mechanisms allow species coexistence and how species densities vary in sympatry compared to allopatry as a result of environmental filtering and biotic interactions.

Hybridization in a warmer world

Climate change is profoundly affecting the evolutionary trajectory of individual species and ecological communities, in part through the creation of novel species assemblages. How climate change will influence competitive interactions has been an active area of research. Far less attention, however, has been given to altered reproductive interactions. Yet, reproductive interactions between formerly isolated species are inevitable as populations shift geographically and temporally as a result of climate change, potentially resulting in introgression, speciation, or even extinction. The susceptibility of hybridization rates to anthropogenic disturbance was first recognized in the 1930s. To date, work on anthropogenically mediated hybridization has focused primarily on either physical habitat disturbance or species invasion. Here, I review recent literature on hybridization to identify how ecological responses to climate change will increase the likelihood of hybridization via the dissolution of species barriers maintained by habitat, time, or behavior. Using this literature, I identify several cases where novel hybrid zones have recently formed, likely as a result of changing climate. Future research should focus on identifying areas and taxonomic groups where reproductive species interactions are most likely to be influenced by climate change. Furthermore, a better understanding of the evolutionary consequences of climate-mediated secondary contact is urgently needed. Paradoxically, hybridization is both a major conservation concern and an important source of novel genetic and phenotypic variation. Hybridization may therefore both contribute to increasing rates of extinction and stimulate the creation of novel phenotypes that will speed adaptation to novel climates. Predicting which result will occur following secondary contact will be an important contribution to conservation for many species.

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.

Species distribution models contribute to determine the effect of climate and interspecific interactions in moving hybrid zones

Climate is a major factor delimiting species' distributions. However, biotic interactions may also be prominent in shaping geographical ranges, especially for parapatric species forming hybrid zones. Determining the relative effect of each factor and their interaction of the contact zone location has been difficult due to the lack of broad scale environmental data. Recent developments in species distribution modelling (SDM) now allow disentangling the relative contributions of climate and species' interactions in hybrid zones and their responses to future climate change. We investigated the moving hybrid zone between the breeding ranges of two parapatric passerines in Europe. We conducted SDMs representing the climatic conditions during the breeding season. Our results show a large mismatch between the realized and potential distributions of the two species, suggesting that interspecific interactions, not climate, account for the present location of the contact zone. The SDM scenarios show that the southerly distributed species, Hippolais polyglotta, might lose large parts of its southern distribution under climate change, but a similar gain of novel habitat along the hybrid zone seems unlikely, because interactions with the other species (H. icterina) constrain its range expansion. Thus, whenever biotic interactions limit range expansion, species may become 'trapped' if range loss due to climate change is faster than the movement of the contact zone. An increasing number of moving hybrid zones are being reported, but the proximate causes of movement often remain unclear. In a global context of climate change, we call for more interest in their interactions with climate change.

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

Local adaptation is often invoked to explain hybrid zone structure, but empirical evidence of this is generally rare. Hybrid zones between two poeciliid fishes, Xiphophorus birchmanni and X. malinche, occur in multiple tributaries with independent replication of upstream-to-downstream gradients in morphology and allele frequencies. Ecological niche modelling revealed that temperature is a central predictive factor in the spatial distribution of pure parental species and their hybrids and explains spatial and temporal variation in the frequency of neutral genetic markers in hybrid populations. Among populations of parentals and hybrids, both thermal tolerance and heat-shock protein expression vary strongly, indicating that spatial and temporal structure is likely driven by adaptation to local thermal environments. Therefore, hybrid zone structure is strongly influenced by interspecific differences in physiological mechanisms for coping with the thermal environment.

Characterizing the physical and genetic structure of the lodgepole pine × jack pine hybrid zone: mosaic structure and differential introgression

Understanding the physical and genetic structure of hybrid zones can illuminate factors affecting their formation and stability. In north-central Alberta, lodgepole pine (Pinus contorta Dougl. ex Loud. var. latifolia) and jack pine (Pinus banksiana Lamb) form a complex and poorly defined hybrid zone. Better knowledge of this zone is relevant, given the recent host expansion of mountain pine beetle into jack pine. We characterized the zone by genotyping 1998 lodgepole, jack pine, and hybrids from British Columbia, Alberta, Saskatchewan, Ontario, and Minnesota at 11 microsatellites. Using Bayesian algorithms, we calculated genetic ancestry and used this to model the relationship between species occurrence and environment. In addition, we analyzed the ancestry of hybrids to calculate the genetic contribution of lodgepole and jack pine. Finally, we measured the amount of gene flow between the pure species. We found the distribution of the pine classes is explained by environmental variables, and these distributions differ from classic distribution maps. Hybrid ancestry was biased toward lodgepole pine; however, gene flow between the two species was equal. The results of this study suggest that the hybrid zone is complex and influenced by environmental constraints. As a result of this analysis, range limits should be redefined.

Pleistocene and ecological effects on continental-scale genetic differentiation in the bobcat (Lynx rufus)

The potential for widespread, mobile species to exhibit genetic structure without clear geographic barriers is a topic of growing interest. Yet the patterns and mechanisms of structure--particularly over broad spatial scales--remain largely unexplored for these species. Bobcats occur across North America and possess many characteristics expected to promote gene flow. To test whether historical, topographic or ecological factors have influenced genetic differentiation in this species, we analysed 1 kb mtDNA sequence and 15 microsatellite loci from over 1700 samples collected across its range. The primary signature in both marker types involved a longitudinal cline with a sharp transition, or suture zone, occurring along the Great Plains. Thus, the data distinguished bobcats in the eastern USA from those in the western half, with no obvious physical barrier to gene flow. Demographic analyses supported a scenario of expansion from separate Pleistocene refugia, with the Great Plains representing a zone of secondary contact. Substructure within the two main lineages likely reflected founder effects, ecological factors, anthropogenic/topographic effects or a combination of these forces. Two prominent topographic features, the Mississippi River and Rocky Mountains, were not supported as significant genetic barriers. Ecological regions and environmental correlates explained a small but significant proportion of genetic variation. Overall, results implicate historical processes as the primary cause of broad-scale genetic differentiation, but contemporary forces seem to also play a role in promoting and maintaining structure. Despite the bobcat's mobility and broad niche, large-scale landscape changes have contributed to significant and complex patterns of genetic structure.

Why do species co-occur? A test of alternative hypotheses describing abiotic differences in sympatry versus allopatry using spadefoot toads

Areas of co-occurrence between two species (sympatry) are often thought to arise in regions where abiotic conditions are conducive to both species and are therefore intermediate between regions where either species occurs alone (allopatry). Depending on historical factors or interactions between species, however, sympatry might not differ from allopatry, or, alternatively, sympatry might actually be more extreme in abiotic conditions relative to allopatry. Here, we evaluate these three hypothesized patterns for how sympatry compares to allopatry in abiotic conditions. We use two species of congeneric spadefoot toads, Spea multiplicata and S. bombifrons, as our study system. To test these hypotheses, we created ecological niche models (specifically using Maxent) for both species to create a map of the joint probability of occurrence of both species. Using the results of these models, we identified three types of locations: two where either species was predicted to occur alone (i.e., allopatry for S. multiplicata and allopatry for S. bombifrons) and one where both species were predicted to co-occur (i.e., sympatry). We then compared the abiotic environment between these three location types and found that sympatry was significantly hotter and drier than the allopatric regions. Thus, sympatry was not intermediate between the alternative allopatric sites. Instead, sympatry occurred at one extreme of the conditions occupied by both species. We hypothesize that biotic interactions in these extreme environments facilitate co-occurrence. Specifically, hybridization between S. bombifrons females and S. multiplicata males may facilitate co-occurrence by decreasing development time of tadpoles. Additionally, the presence of alternative food resources in more extreme conditions may preclude competitive exclusion of one species by the other. This work has implications for predicting how interacting species will respond to climate change, because species interactions may facilitate survival in extreme habitats.

Relationships between spatio-temporal environmental and genetic variation reveal an important influence of exogenous selection in a pupfish hybrid zone

The importance of exogenous selection in a natural hybrid zone between the pupfishes Cyprinodon atrorus and Cyprinodon bifasciatus was tested via spatio-temporal analyses of environmental and genetic change over winter, spring and summer for three consecutive years. A critical influence of exogenous selection on hybrid zone regulation was demonstrated by a significant relationship between environmental (salinity and temperature) and genetic (three diagnostic nuDNA loci) variation over space and time (seasons) in the Rio Churince system, Cuatro Ciénegas, Mexico. At sites environmentally more similar to parental habitats, the genetic composition of hybrids was stable and similar to the resident parental species, whereas complex admixtures of parental and hybrid genotypic classes characterized intermediate environments, as did the greatest change in allelic and genotypic frequencies across seasons. Within hybrids across the entire Rio Churince system, seasonal changes in allelic and genotypic frequencies were consistent with results from previous reciprocal transplant experiments, which showed C. bifasciatus to suffer high mortality (75%) when exposed to the habitat of C. atrorus in winter (extreme temperature lows and variability) and summer (abrupt salinity change and extreme temperature highs and variability). Although unconfirmed, the distributional limits of C. atrorus and C. atrorus-like hybrids appear to be governed by similar constraints (predation or competition). The argument favouring evolutionary significance of hybridization in animals is bolstered by the results of this study, which links the importance of exogenous selection in a contemporary hybrid zone between C. atrorus and C. bifasciatus to previous demonstration of the long-term evolutionary significance of environmental variation and introgression on the phenotypic diversification Cuatro Ciénegas Cyprinodon.

The Role of Environmental Heterogeneity in Maintaining Reproductive Isolation between HybridizingPasserina(Aves: Cardinalidae) Buntings

Hybrid zones are useful systems in which to investigate processes important in creating and maintaining biological diversity. As they are often located in ecotones, patterns of environmental heterogeneity may influence hybridization, and may also influence the maintenance of reproductive isolation between hybridizing species. Focusing on the hybrid zone betweenPasserina amoena(Lazuli Bunting) andPasserina cyanea(Indigo Bunting), located in the eastern Rocky Mountain/western Great Plains ecotone, we examined the relationship between population-pairwise differences in the proportion of hybrids and environmental variation. Models including environmental variables explained more of the variation in hybridization rates across the ecotone than did models that only included the geographic distance between sampling localities as predictor variables (63.9% and 58.9% versus 38.8% and 39.9%, depending on how hybridization was quantified). In the models including environmental variables, the amount of rainfall during the warmest quarter had the greatest explanatory power, consistent with a hypothesis thatP. cyaneais better adapted to the mesic environments of eastern North America andP. amoenais better adapted to the xeric habitats of western North America. These results suggest that continued reproductive isolation between these species is mediated, at least partially, by differential adaptations to local environmental conditions.