Hybrid Zone Research Articles

Structure of the hybrid zone between allied species of the common vole, Microtus arvalis and M. obscurus: Influence of genetic factors and landscape-geographic conditions

Patterns of introgression of several genetic markers across the hybrid zone between allied species of the common vole Microtus arvalis s. str. и M. obscurus were studied in four its sections: in northwest Nizhny Novgorod region, east Vladimir and southwest of Nizhny Novgorod regions, south Lipetsk region, and northwest Voronezh region. Analysis of the clinal variability for three molecular-genetic markers (cytb, tp53, SMCY11) and for karyotypes showed a structural similarity between the “Vladimir – Nizhny Novgorod”, “Nizhny Novgorod”, and “Voronezh” sections. The maximal width was shown for the cytb cline, the minimal width – for the SMCY11 cline; the tp53 cline and chromosomal cline occupy intermediate position for this parameter. Furthermore, in these transects the centre of the cline for the cytb is shifted southeastward (into the distribution range of M. obscurus) from the centres of three other clines. The revealed asymmetric introgression of mitochondrial genome from M. arvalis to M. obscurus may be explained by the fact that the hybrid zone was formed as a result of invasion of M. obscurus into the range of M. arvalis. The “Lipetsk” transect differs from three above-mentioned transects by very narrow clines with nearly coinciding centres. Such characteristics of the “Lipetsk” transect are obviously caused by localization of the hybrid zone in this section along river Voronezh. The obtained results led us to suppose that the structure of the studied hybrid zone is determined mainly by coinciding (or non-coinciding) of its centre with local physical barriers.

Genome-wide SNP assessment of contemporary European red deer genetic structure highlights the distinction of peripheral populations and the main admixture zones in Europe.

Genome-wide technologies open up new possibilities to clarify questions on genetic structure and phylogeographic history of taxa previously studied with microsatellite loci and mitochondrial sequences. Here, we used 736 individual red deer (Cervus elaphus) samples genotyped at 35,701 single nucleotide polymorphism loci (SNPs) to assess the population structure of the species throughout Europe. The results identified 28 populations, with higher degrees of genetic distinction in peripheral compared to mainland populations. Iberian red deer show high genetic differentiation, with lineages in Western and Central Iberia maintaining their distinctiveness, which supports separate refugial ranges within Iberia along with little recent connection between Iberian and the remaining Western European populations. The Norwegian population exhibited the lowest variability and the largest allele frequency differences from mainland European populations, compatible with a history of bottlenecks and drift during post-glacial colonization from southern refugia. Scottish populations showed high genetic distance from the mainland but high levels of diversity. Hybrid zones were found between Eastern and Western European lineages in Central Europe as well as in the Pyrenees, where red deer from France are in close contact with Iberian red deer. Anthropogenic restocking has promoted the Pyrenean contact zone, admixture events in populations on the Isle of Rum and in the Netherlands, and at least partly the admixture of the two main lineages in central-eastern Europe. Our analysis enabled detailed resolution of population structure of a large mammal widely distributed throughout Europe and contributes to resolving the evolutionary history, which can also inform conservation and management policies.

Reproductive character displacement: insights from genital morphometrics in damselfly hybrid zones.

Reproductive Character Displacement (RCD) refers to the phenomenon of greater differences in reproductive characters between two species when they occur in sympatry compared to when they occur in allopatry to prevent maladaptive hybridization. We explored whether reinforcement of a mechanical barrier involved in the first contact point between male and female genital traits during copulation in the cross between Ischnura graellsii males and Ischnura elegans females has led to RCD, and whether it supports the lock-and-key hypothesis of genital evolution. We employed geometric morphometrics to analyze the shape and size of male and female genital traits, controlling for environmental and geographic factors. Consistent with an increase in mechanical isolation via reinforcement, we detected larger divergence in genital traits between the species in sympatry than in allopatry, and also stronger signal in females than in males. In the Northwest (NW) hybrid zone, we detected RCD in I. graellsii males and I. elegans females, while in the Northcentral (NC) hybrid zone we detected RCD only in I. elegans females and I. elegans males. The detection of RCD in both sexes of I. elegans was consistent with the lock-and-key hypothesis of genital evolution via female choice for conspecific males in this species. Our study highlights the importance of using geometric morphometrics to deal with the complexity of female reproductive structures while controlling for environmental and geographic factors to investigate RCD. This study contributes valuable insights into the dynamics of reproductive isolation mechanisms and genital coevolution.

Recurrent hybridization and gene flow shaped Norway and Siberian spruce evolutionary history over multiple glacial cycles.

Most tree species underwent cycles of contraction and expansion during the Quaternary. These cycles led to an ancient and complex genetic structure that has since been affected by extensive gene flow and by strong local adaptation. The extent to which hybridization played a role in this multi-layered genetic structure is important to be investigated. To study the effect of hybridization on the joint population genetic structure of two dominant species of the Eurasian boreal forest, Picea abies and P. obovata, we used targeted resequencing and obtained around 480 K nuclear SNPs and 87 chloroplast SNPs in 542 individuals sampled across most of their distribution ranges. Despite extensive gene flow and a clear pattern of Isolation-by-Distance, distinct genetic clusters emerged, indicating the presence of barriers and corridors to migration. Two cryptic refugia located in the large hybrid zone between the two species played a critical role in shaping their current distributions. The two species repeatedly hybridized during the Pleistocene and the direction of introgression depended on latitude. Our study suggests that hybridization helped both species to overcome main shifts in their distribution ranges during glacial cycles and highlights the importance of considering whole species complex instead of separate entities to retrieve complex demographic histories.

Asymmetric mate preference and reproductive interference mediate climate-induced changes in mate availability in a small mammal hybrid zone.

Range expansion and contraction are among the most common biotic responses to changing environmental conditions, yet much is to be learned about the mechanisms that underlie range-edge population dynamics, especially when those areas are points of secondary contact between closely related species. Here, we present field-measured parentage data that document the reproductive outcomes of changes in mate availability at a secondary contact zone between two species of woodrat in the genus Neotoma. Changes in mate availability resulted from drought-driven differential survival between the species and their hybrids. As the availability of conspecific mates declined, rates of hybridization increased, leading to the accumulation of admixed individuals in the zone of contact. Patterns of reproductive success in the wild appear to be the result of a combination of both pre-mating isolation and post-zygotic selection resulting from genomic incompatibilities between the parental lineages. Evidence of asymmetric mate preference between the parental lineages came from both skewed reproductive output in the field and laboratory preference trials. Moreover, partial genomic incompatibility was evident from the near-zero reproductive success of F1 males and because nearly all surviving hybrids had one pure parent. Nonetheless, high reproductive success of F1 females and backcrossing in both parental directions allow for introgression between the parental species. These findings reveal how climate change may alter evolutionary outcomes for species at the edge of their ranges through an interplay of behavioral, demographic, and genetic mechanisms.

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.

Higher locomotor costs of pregnancy in viviparous compared to oviparous common lizards (Zootoca vivipara).

Pregnancy is a physiological cost of reproduction for animals that rely on fleeing to avoid predators. Costs of reproduction are predicted to differ between alternative reproductive strategies or modes, such as egg-laying (oviparity) or live-bearing (viviparity). However, disentangling the factors that comprise this cost and how it differs for oviparous or viviparous females is challenging due to myriad environmental, biological, and evolutionary confounds. Here, we tested the hypothesis that costs of pregnancy differ between oviparous and viviparous common lizards (Zootoca vivipara). We predicted that the degree of locomotor impairment during pregnancy and therefore the cost of reproduction would be higher in viviparous females. We conducted our experiment in a hybrid zone containing oviparous and viviparous common lizards. Due to the common environment and inclusion of hybrid individuals, we could infer whether differences were inherent to parity mode. We found that the average and maximum running speed of pregnant females was slower than after they had given birth or laid eggs. Viviparous females experienced an additional pregnancy weight burden and for a longer time period, but were not slower at running than pregnant oviparous females. In addition, we found a parity mode-specific effect of reproductive investment; producing larger clutches was costlier for the locomotor performance of viviparous females for reasons other than the mass increase. Locomotor costs were found to be intermediate in hybrid females, indicating that they are specific to each reproductive mode. Our study shows that viviparous females experience an additional physical and physiological cost of pregnancy and reproductive investment. This two-fold cost implies that viviparous females modulate resource allocation decisions and/or adjust their behavioural responses that result from locomotor impairment.

Discordance between taxonomy and population genomic data: An avian example relevant to the United States Endangered Species Act.

Population genomics can reveal cryptic biological diversity that may impact fitness while simultaneously serving to delineate relevant conservation units. Here, we leverage the power of whole-genome resequencing for conservation by studying 433 individual lesser prairie-chicken (Tympanuchus pallidicinctus; LEPC, a federally endangered species of conservation concern in the United States) and greater prairie-chicken (Tympanuchus cupido; GRPC, a legally huntable species throughout much of its range). The genomic diversity of two formally recognized distinct population segments (DPSs) of LEPCs is similar, but they are genetically distinct. Neither DPS is depleted of its genomic diversity, neither is especially inbred, and temporal diversity is relatively stable in both conservation units. Interspecific differentiation between the two species was only slightly higher than that observed between LEPC DPSs, due largely to bidirectional introgression. The high resolution provided by our dataset identified a genomic continuum between the two species such that individuals sampled from the hybrid zone were imperfectly assigned to their presumptive species when considering only their physical characteristics. The admixture between the two species is reflected in the spectrum of individual ancestry coefficients, which has legal implications for the "take" of individuals under the Endangered Species Act. Overall, our data highlight the recurring dissonance between static policies and dynamic species boundaries that are increasingly obvious in the population genomic era.

Climate change differentially alters distribution of two marten species in a hybrid zone.

Species' ranges are shifting rapidly with climate change, altering the composition of biological communities and interactions within and among species. Hybridization is among the species interactions that may change markedly with climate change, yet it is understudied relative to others. We used non-invasive genetic detections to build a maximum entropy species distribution model and investigate the factors that delimit the present and future ranges of American marten (Martes americana) and Pacific marten (Martes caurina) in a contact zone in the Northern Rockies. We found that climate change will decrease the suitable habitat predicted for both species, as well as the amount of overlap in predicted suitable habitat between the species. Interestingly, predicted suitable habitat for Pacific marten extended further north in the study region than our genetic detections for the species, suggesting that biotic factors, such as interactions with American marten, may affect the realized range of this species. Our results suggest that future work investigating the interactions among biotic and abiotic factors that influence hybrid zone dynamics is important for predicting the futures of these two species in this area under climate change.

Waves of Colonization and Gene Flow in a Great Speciator.

Secondary contact between previously allopatric lineages offers a test of reproductive isolating mechanisms that may have accrued in isolation. Such instances of contact can produce stable hybrid zones-where reproductive isolation can further develop via reinforcement or phenotypic displacement-or result in the lineages merging. Ongoing secondary contact is most visible in continental systems, where steady input from parental taxa can occur readily. In oceanic island systems, however, secondary contact between closely related species of birds is relatively rare. When observed on sufficiently small islands, relative to population size, secondary contact likely represents a recent phenomenon. Here, we examine the dynamics of a group of birds whose apparent widespread hybridization influenced Ernst Mayr's foundational work on allopatric speciation: the whistlers of Fiji (Aves: Pachycephala). We demonstrate two clear instances of secondary contact within the Fijian archipelago, one resulting in a hybrid zone on a larger island, and the other resulting in a wholly admixed population on a smaller, adjacent island. We leveraged low genome-wide divergence in the hybrid zone to pinpoint a single genomic region associated with observed phenotypic differences. We use genomic data to present a new hypothesis that emphasizes rapid plumage evolution and post-divergence gene flow.

Accelerated mitochondrial evolution and asymmetric fitness of hybrids contribute to the persistence of Helix thessalica in the Helix pomatia range.

Interbreeding and introgression between recently diverged species is common. However, the processes that prevent these species from merging where they co-occur are not well understood. We studied the mechanisms that allowed an isolated group of populations of the snail Helix thessalica to persist within the range of the related Helix pomatia despite high gene flow. Using genomic cline analysis, we found that the nuclear gene flow between the two taxa across the mosaic hybrid zone was not different from that expected under neutral admixture, but that the exchange of mtDNA was asymmetric. Tests showed that there is relaxed selection in the mitochondrial genome of H. thessalica and that the substitution rate is elevated compared to that of H. pomatia. A lack of hybrids that combine the mtDNA of H. thessalica with a mainly (>46%) H. pomatia genomic background indicates that the nuclear-encoded mitochondrial proteins of H. pomatia are not well adapted to the more rapidly evolving proteins and RNAs encoded by the mitochondrion of H. thessalica. The presumed reduction of fitness of hybrids with the fast-evolving mtDNA of H. thessalica and a high H. pomatia ancestry, similar to 'Darwin's Corollary to Haldane's rule', resulted in a relative loss of H. pomatia nuclear ancestry compared to H. thessalica ancestry in the hybrid zone. This probably prevents the H. thessalica populations from merging quickly with the surrounding H. pomatia populations and supports the hypothesis that incompatibilities between rapidly evolving mitochondrial genes and nuclear genes contribute to speciation.

Hybridization and introgression in deeply differentiated salamander species – molecular genetics and a reappraisal of Dr. Louis Vallée’s osteological data

Abstract Two morphologically and genetically distinct salamander species (the crested newt, Triturus cristatus and the marbled newt, T. marmoratus) engage over a large area in the west of France where they hybridize at ca. 4%. The species interaction is characterized by ecological differentiation and limited gene flow beyond the F1 hybrid generation. Incompletely isolated species like these allow to investigate the genetic mechanisms and evolutionary forces that maintain their identity in the face of ongoing gene flow, to which a large mosaic hybrid zone provides excellent opportunities. A reanalysis of published morphological data supports the partial breakdown of the species barrier whereas extensive genetic data show that introgression is low yet asymmetric, in line with the dynamics of species replacement inferred earlier. The current work, with seven diagnostic nuclear markers studied for a large sample, revises the estimates of introgression in both T. cristatus (to 0.24%) and T. marmoratus (to 0.11%). Difficulties remain in the recognition of potential triploid hybrids versus backcross hybrids. Haldane’s rule is partially supported, but deeper analyses require the use of a molecular marker for sex that is not yet available.

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.

Contrasting effect of hybridization on genetic differentiation in three rockfish species with similar life history.

Hybridization can provide evolutionary benefits (e.g., population resilience to climate change) through the introduction of adaptive alleles and increase of genetic diversity. Nevertheless, management strategies may be designed based only on the parental species within a hybrid zone, without considering the hybrids. This can lead to ineffective spatial management of species, which can directly harm population diversity and negatively impact food webs. Three species of rockfish (Brown Rockfish (Sebastes caurinus), Copper Rockfish (S. auriculatus), and Quillback Rockfish (S. maliger)) are known to hybridize within Puget Sound, Washington, but genetic data from these species are used to infer population structure in the entire genus, including in species that do not hybridize. The goal of this project was to estimate the hybridization rates within the region and determine the effect of hybridization on geographic patterns of genetic structure. We sequenced 290 Brown, Copper, and Quillback rockfish using restriction-site associated DNA sequencing (RADseq) from four regions within and outside Puget Sound, Washington. We show that (i) hybridization within Puget Sound was asymmetrical, not recent, widespread among individuals, and relatively low level within the genome, (ii) hybridization affected population structure in Copper and Brown rockfish, but not in Quillback Rockfish and (iii) after taking hybridization into account we found limited directional dispersal in Brown and Copper rockfish, and evidence for two isolated populations in Quillback Rockfish. Our results suggest that rockfish population structure is species-specific, dependent on the extent of hybridization, and cannot be inferred from one species to another despite similar life history.

Unraveling the genomic landscape of Campylorhynchus wrens along western Ecuador's precipitation gradient: Insights into hybridization, isolation by distance, and isolation by the environment.

Environmental gradients have the potential to influence genetic differentiation among populations ultimately leading to allopatric speciation. However, environmental gradients can also facilitate hybridization between closely related taxa. We investigated a putative hybrid zone in western Ecuador, involving two polytypic wren species (Aves: Troglodytidae), Campylorhynchus zonatus and C. fasciatus. Our study addressed two primary questions: (1) Is there evidence of population structure and genetic admixture between these taxa in western Ecuador? and (2) What are the relative contributions of isolation by distance and isolation by the environment to the observed genetic differentiation along the environmental gradient in this region? We analyzed 4409 single-nucleotide polymorphisms (SNPs) from 112 blood samples sequenced using ddRadSeq and a de novo assembly. The optimum number of genetic clusters ranged from 2 to 4, aligning with geographic origins, known phylogenetics, and physical or ecological constraints. We observed notable transitions in admixture proportions along the environmental gradient in western Ecuador between C. z. brevirostris and the northern and southern genetic clusters of C. f. pallescens. Genetic differentiation between the two C. f. pallescens populations could be attributed to an unreported potential physical barrier in central western Ecuador, where the proximity of the Andes to the coastline restricts lowland habitats, limiting dispersal and gene flow, especially among dry-habitat specialists. The observed admixture in C. f. pallescens suggests that this subspecies may be a hybrid between C. z. brevirostris and C. fasciatus, with varying degrees of admixture in western Ecuador and northwestern Peru. We found evidence of isolation by distance, while isolation by the environment was less pronounced but still significant for annual mean precipitation and precipitation seasonality. This study enhances our understanding of avian population genomics in tropical regions.

Ecological displacement in a Rocky Mountain hybrid zone informs management of North American martens (Martes)

ContextParapatric sister species are ideal for tests of ecological interactions. Pacific (Martes caurina) and American pine (M. americana) martens are economically and culturally valuable furbearers that hybridize in the north-central Rocky Mountains. Despite preliminary evidence of biased introgression, the hybrid zone has been geographically stable for 70 years, but interspecific ecological interactions have yet to be examined in detail.ObjectivesWe test whether ecological interactions may influence the outcome of hybridization in this system. To that end, we estimate the fundamental niche of each species and gauge how suitability landscapes change when the two species are in contact.MethodsWe genotyped > 400 martens from the Rocky Mountain hybrid zone to diagnose individuals to species-level and identify putative hybrids. We then built range-wide ecological niche models for each species, excluding individuals in the hybrid zone, to approximate their respective fundamental niches. Those models were projected into the hybrid zone and compared with niche models trained on individuals within the hybrid zone to assess how niche dynamics change when the species are in sympatry.ResultsThe fundamental niche of each species differed significantly, while the hybrid zone was equally suitable for both. Niches of each species based on models built within the hybrid zone showed that Pacific martens utilized significantly less suitable habitat than expected based on their range-wide fundamental niche, suggesting that species interactions shape local hybridization. We detected few admixed individuals (12%), with no evidence of directional (sex or species) biases. Interstate-90 further acts as a major dispersal barrier.ConclusionsNorth American martens are currently managed as a single species by some state agencies, yet significant ecological and genetic differences indicate they should be managed separately. The observed ecological displacement of Pacific martens by American pine martens may partially explain the mixed success of historical, mixed-species wildlife translocations and cautions such translocations in the future. Landscape-scale consideration of ecological dynamics, in addition to molecular compatibility, will be essential to the success of future translocations.

Genetic identification, morphology and distribution of Natrixhelvetica subspecies in southern and western Switzerland (Reptilia, Squamata, Serpentes).

Most of Switzerland is inhabited by the nominotypical subspecies of the barred grass snake (Natrixhelveticahelvetica), which is characterized by mitochondrial DNA lineage E. Only in the northeast of the country, the common grass snake (N.natrix) occurs and hybridizes with N.h.helvetica in a narrow contact zone. However, we discovered that in southern and western Switzerland barred grass snakes representing another mtDNA lineage (lineage C) are widely distributed. Lineage C is typical for Alpine populations of the southern subspecies N.h.sicula. Our microsatellite analyses of the Swiss samples revealed differences between the two subspecies and also a substructure with two clusters in each subspecies. Furthermore, we discovered a contact and hybrid zone of N.h.helvetica and N.h.sicula along the northern shore of Lake Geneva and also confirm that interbreeding with alien common grass snakes (N.n.moreotica, mtDNA lineage 7) occurs there. This finding is of concern for nature conservation and measures should be taken to prevent further genetic pollution. Using morphometrics, we found no differences between the two subspecies of N.helvetica, while N.natrix was slightly distinct from N.helvetica.

Digest: Asymmetrical introgression of belly coloration across a manakin hybrid zone in Panama.

How does sexual selection impact introgression dynamics across a hybrid zone? Long et al. (2024) used historical (1989-1994) and contemporary (2017-2020) samples to quantify the stability of a Panamanian hybrid zone between golden-collared manakins (Manacus vitellinus) and white-collared manakins (M. candei). Their analyses revealed a spatially stable hybrid zone, except for one trait, belly plumage coloration, which has introgressed into the distribution of the white-collared manakin. This finding suggests possible sexual selection for this trait.

Fluctuating selection in a Monkeyflower hybrid zone.

While hybridization was viewed as a hindrance to adaptation and speciation by early evolutionary biologists, recent studies have demonstrated the importance of hybridization in facilitating evolutionary processes. However, it is still not well-known what role spatial and temporal variation in natural selection play in the maintenance of naturally occurring hybrid zones. To identify whether hybridization is adaptive between two closely related monkeyflower species, Mimulus guttatus and Mimulus laciniatus, we performed repeated reciprocal transplants between natural hybrid and pure species' populations. We planted parental genotypes along with multiple experimental hybrid generations in a dry (2021) and extremely wet (2023) year in the Sierra Nevada, CA. By taking fine scale environmental measurements, we found that the environment of the hybrid zone is more similar to M. laciniatus's seasonally dry rocky outcrop habitat than M. guttatus's moist meadows. In our transplants hybridization does not appear to be maintained by a consistent fitness advantage of hybrids over parental species in hybrid zones, but rather a lack of strong selection against hybrids. We also found higher fitness of the drought adapted species, M. laciniatus, than M. guttatus in both species' habitats, as well as phenotypic selection for M. laciniatus-like traits in the hybrid habitat in the dry year of our experiment. These findings suggest that in this system hybridization might function to introduce drought-adapted traits and genes from M. laciniatus into M. guttatus, specifically in years with limited soil moisture. However, we also find evidence of genetic incompatibilities in second generation hybrids in the wetter year, which may balance a selective advantage of M. laciniatus introgression. Therefore, we find that hybridization in this system is both potentially adaptive and costly, and that the interaction of positive and negative selection likely determines patterns of gene flow between these Mimulus species.

Genomic insights into hybrid zone formation: The role of climate, landscape, and demography in the emergence of a novel hybrid lineage.

Population demographic changes, alongside landscape, geographic and climate heterogeneity, can influence the timing, stability and extent of introgression where species hybridise. Thus, quantifying interactions across diverged lineages, and the relative contributions of interspecific genetic exchange and selection to divergence at the genome-wide level is needed to better understand the drivers of hybrid zone formation and maintenance. We used seven latitudinally arrayed transects to quantify the contributions of climate, geography and landscape features to broad patterns of genetic structure across the hybrid zone of Populus trichocarpa and P. balsamifera and evaluated the demographic context of hybridisation over time. We found genetic structure differed among the seven transects. While ancestry was structured by climate, landscape features influenced gene flow dynamics. Demographic models indicated a secondary contact event may have influenced contemporary hybrid zone formation with the origin of a putative hybrid lineage that inhabits regions with higher aridity than either of the ancestral groups. Phylogenetic relationships based on chloroplast genomes support the origin of this hybrid lineage inferred from demographic models based on the nuclear data. Our results point towards the importance of climate and landscape patterns in structuring the contact zones between P. trichocarpa and P. balsamifera and emphasise the value whole genome sequencing can have to advancing our understanding of how neutral processes influence divergence across space and time.