Sexual selection drives asymmetric introgression in wall lizards

Sexual color ornamentation, microhabitat choice, and thermal physiology in the common wall lizard (Podarcis muralis)

Common wall lizards (Podarcis muralis) in Italy show a striking variation in body coloration across the landscape, with highly exaggerated black and green colors in hot and dry climates and brown and white colors in cool and wet climates. Males are more intensely colored than females, and previous work has suggested that the maintenance of variation in coloration across the landscape reflects climatic effects on the strength of male-male competition, and through this sexual selection. However climatic effects on the intensity of male-male competition would need to be exceptionally strong to fully explain the geographic patterns of color variation. Thus, additional processes may contribute to the maintenance of color variation. Here we test the hypothesis that selection for green and black ornamentation in the context of male-male competition is opposed by selection against ornamentation because the genes involved in the regulation of coloration have pleiotropic effects on thermal physiology, such that ornamentation is selected against in cool climates. Field observations revealed no association between body coloration and microhabitat use or field active body temperatures. Consistent with these field data, lizards at the extreme ends of the phenotypic distribution for body coloration did not show any differences in critical minimum temperature, preferred body temperature, temperature-dependent metabolic rate, or evaporative water loss when tested in the laboratory. Combined, these results provide no evidence that genes that underlie sexual ornamentation are selected against in cool climate because of pleiotropic effects on thermal biology.

Mitonuclear and phenotypic discordance in an Atlantic Forest frog hybrid zone

Discordance between mitochondrial and nuclear DNA is common among animals and can be the result of a number of evolutionary processes, including incomplete lineage sorting and introgression. Particularly relevant in contact zones, mitonuclear discordance is expected because the mitochondrial genome is haploid and primarily uniparentally inherited, whereas nuclear loci are evolving at slower rates. In addition, when closely related taxa come together in hybrid zones, the distribution of diagnostic phenotypic characters and their concordance with the mitochondrial or nuclear lineages can also inform on historical and ongoing dynamics within hybrid zones. Overall, genetic and phenotypic discordances provide evidence for evolutionary divergence and processes that maintain boundaries among sister species or lineages. In this study, we characterized patterns of genetic and phenotypic variation in a contact zone between Cycloramphus dubius and Cycloramphus boraceiensis, two sister species of frogs endemic to the Atlantic Coastal Forest of Brazil. We examined genomic-scale nuclear diversification across 19 populations, encompassing the two parental forms and a contact zone between them. We compared the distribution of genomic DNA variability with that of a mitochondrial locus (16S) and two morphological traits (dorsal tubercles and body size). Our results reveal multiple divergent lineages with ongoing admixture. We detected discordance in patterns of introgression across the three data types. Cycloramphus dubius males are significantly larger than C. boraceiensis males, and we posit that competition among males in the hybrid zone, coupled with mate choice by females, may be one mechanism leading to patterns of introgression observed between the species.

Adaptive introgression reveals the genetic basis of a sexually selected syndrome in wall lizards

The joint expression of particular colors, morphologies, and behaviors is a common feature of adaptation, but the genetic basis for such "phenotypic syndromes" remains poorly understood. Here, we identified a complex genetic architecture associated with a sexually selected syndrome in common wall lizards, by capitalizing on the adaptive introgression of coloration and morphology into a distantly related lineage. Consistent with the hypothesis that the evolution of phenotypic syndromes in vertebrates is facilitated by developmental linkage through neural crest cells, most of the genes associated with the syndrome are involved in neural crest cell regulation. A major locus was a ~400-kb region, characterized by standing structural genetic variation and previously implied in the evolutionary innovation of coloration and beak size in birds. We conclude that features of the developmental and genetic architecture contribute to maintaining trait integration, facilitating the extensive and rapid introgressive spread of suites of sexually selected characters.

The ecological importance of hybridization

Hybridization as an evolutionary process has been studied in depth over the past few decades. Research has focused on its role in shaping reproductive barriers, its adaptive value, and its genomic consequences. In contrast, our knowledge of ecological dimensions of hybridization is still in its infancy, despite hybridization being an inherently ecological interaction. Using examples from various organisms, we show that hybridization can affect and be affected by non-reproductive interactions, including predation, competition, parasitism, mutualism, commensalism, and organism-environment interactions, with significant implications for community structure and ecosystem functioning. However, since these dimensions of hybridization have mostly been revealed from studies designed to decipher other evolutionary processes, we argue that much of the eco-evolutionary importance of hybridization is yet to be discovered.

Signalling on islands: the case of Lilford’s wall lizard (Podarcis lilfordi gigliolii) from Dragonera

Studies of the effects of insularity on animal signals are scarce, particularly in lizards. Here, we use Lilford’s wall lizard from Dragonera (Podarcis lilfordi gigliolii) to ask how island conditions have affected its repertoire of social signals, focusing on two visual signals shared by many Podarcis species: ultraviolet (UV)–blue-reflecting ventrolateral colour patches and visual displays. We examined whether the number or spectral characteristics of the UV–blue patches are associated with traits related to individual quality. We also used visual models to assess visual conspicuousness and to measure sexual dichromatism. We did not observe foot shakes or any other visual displays usually found in continental Podarcis. We found that none of the UV–blue patch variables covaried with morphometric variables indicative of fighting ability or body condition in males, suggesting that this coloration does not signal individual quality. We also found very little sexual dichromatism. In particular, the UV–blue patches of females seem over-expressed and more similar to those of males than those of continental Podarcis. Ancestral state reconstruction reveals that the lack of sexual dimorphism in the UV–blue patches is a derived condition for P. lilfordi gigliolii and other Podarcis living on small islands. Our results thus show a pattern of reduced social signalling in P. lilfordi gigliolii relative to mainland Podarcis, with some signals being lost or under-expressed (visual displays) and others losing their signalling function (UV–blue patches). We hypothesize that these changes are attributable to the high population density of P. lilfordi gigliolii, which discourages territorial behaviour and promotes extreme social tolerance, making most social signals unnecessary. More work will be needed to determine whether this is a common pattern in lizards inhabiting small and densely populated islands.

Large scale phenotypic characterisation of Hierophis viridiflavus (Squamata: Serpentes): climatic and environmental drivers suggest the role of evolutionary processes in a polymorphic species

Colour variability is largely widespread in the animal world as it is tightly associated with fitness and survivorship. Therefore, the drivers and implications of such variability have been of great interest for zoologists in the past decades. Reptiles are excellent models to investigate colour variations and expression under different conditions. Here, we focused on melanism occurrence in the two main lineages of Hierophis viridiflavus at the scale of the species distribution, by extracting available data from iNaturalist, a citizen science network, with the aim of detecting any pure effect of climate or local habitat on colour expression. Our analyses highlighted that habitat does not explain differences in phenotypes, whereas marked effects of geographic and climatic variables were detected. However, the observed climatic effects could be a proxy of the geographical distribution of the two groups, and thus the high occurrence of bright colourations in western populations of the eastern lineage could be addressed to an ongoing event of asymmetric gene flow in contact zones. The current distribution of phenotypes could be the outcome of the evolutionary history of the species combined with the geological history of the Mediterranean region. This investigation, though, is only preliminary and molecular analyses on highly variable regions of the genome are mandatory to address this issue.

Reproductive isolating mechanisms contributing to asymmetric hybridization in Killifishes (Fundulus spp.)

When species hybridize, one F1 hybrid cross type often predominates. Such asymmetry can arise from differences in a variety of reproductive barriers, but the relative roles and concordance of pre-mating, post-mating prezygotic, and post-zygotic barriers in producing these biases in natural animal populations have not been widely investigated. Here, we study a population of predominantly F1 hybrids between two killifish species (Fundulus heteroclitus and F. diaphanus) in which >95% of F1 hybrids have F. diaphanus mothers and F. heteroclitus fathers (D♀ × H♂). To determine why F. heteroclitus × F. diaphanus F1 hybrids (H♀ × D♂) are so rare, we tested for asymmetry in pre-mating reproductive barriers (female preference and male aggression) at a common salinity (10 ppt) and post-mating, pre-zygotic (fertilization success) and post-zygotic (embryonic development time and hatching success) reproductive barriers at a range of ecologically relevant salinities (0, 5, 10, and 15 ppt). We found that F. heteroclitus females preferred conspecific males, whereas F. diaphanus females did not, matching the observed cross bias in the wild. Naturally rare H♀ × D♂ crosses also had lower fertilization success than all other cross types, and a lower hatching success than the prevalent D♀ × H♂ crosses at the salinity found in the hybrid zone centre (10 ppt). Furthermore, the naturally predominant D♀ × H♂ crosses had a higher hatching success than F. diaphanus crosses at 10 ppt, which may further increase their relative abundance. The present study suggests that a combination of incomplete mating, post-mating pre-zygotic and post-zygotic reproductive isolating mechanisms act in concert to produce hybrid asymmetry in this system.

Experimentally testing mate preference in an avian system with unidirectional bill color introgression

Mating behavior can play a key role in speciation by inhibiting or facilitating gene flow between closely related taxa. Hybrid zones facilitate a direct examination of mating behavior and the traits involved in establishing species barriers. The long-tailed finch (Poephila acuticauda) has two hybridizing subspecies that differ in bill color (red and yellow), and the yellow bill phenotype appears to have introgressed ~350 km eastward following secondary contact. To examine the role of mate choice on bill color introgression, we performed behavioral assays using natural and manipulated bill colors. We found an assortative female mating preference for males of their own subspecies when bill color was not manipulated. However, we did not find this assortative preference in trials based on artificially manipulated bill color. This could suggest that assortative preference is not fixed entirely on bill color and instead may be based on a different trait (e.g., song) or a combination of traits, or alternatively may be due to lower statistical power alongside the bill manipulations being unconvincing to the female choosers. Intriguingly, we find a bias in the inheritance of bill color in captive bred F1 hybrid females. Previous modeling suggests that assortative mate preference and this kind of partial dominance in the underlying genes may together contribute to introgression, making the genetic architecture of bill color in this system a priority for future research.

Conflict over fertilization underlies the transient evolution of reinforcement

When two species meet in secondary contact, the production of low fitness hybrids may be prevented by the adaptive evolution of increased prezygotic isolation, a process known as reinforcement. Theoretical challenges to the evolution of reinforcement are generally cast as a coordination problem, i.e., "how can statistical associations between traits and preferences be maintained in the face of recombination?" However, the evolution of reinforcement also poses a potential conflict between mates. For example, the opportunity costs to hybridization may differ between the sexes or species. This is particularly likely for reinforcement based on postmating prezygotic (PMPZ) incompatibilities, as the ability to fertilize both conspecific and heterospecific eggs is beneficial to male gametes, but heterospecific mating may incur a cost for female gametes. We develop a population genetic model of interspecific conflict over reinforcement inspired by "gametophytic factors", which act as PMPZ barriers among Zea mays subspecies. We demonstrate that this conflict results in the transient evolution of reinforcement-after females adaptively evolve to reject gametes lacking a signal common in conspecific gametes, this gamete signal adaptively introgresses into the other population. Ultimately, the male gamete signal fixes in both species, and isolation returns to pre-reinforcement levels. We interpret geographic patterns of isolation among Z. mays subspecies considering these findings and suggest when and how this conflict can be resolved. Our results suggest that sexual conflict over fertilization may pose an understudied obstacle to the evolution of reinforcement.

Unravelling the origin of the common wall lizards (Podarcismuralis) in south-eastern Europe using mitochondrial evidence

The origin of the common wall lizards (Podarcismuralis) populations in south-eastern Europe (namely in Bulgaria and Romania), representing the north-eastern range border of this species, was addressed using mitochondrial DNA. We compared cytochrome b sequences from Bulgaria and Romania with those from the contiguous range in Central Europe that are available from previous studies. We recorded five main haplogroups in Bulgaria and Romania, belonging to the Central Balkan clade. However, haplogroup III was recorded in more localities than previously found. Additionally, signs of haplotype admixture were identified in several populations along the Danube River. The presence of the Southern Alps haplotype in one population from Otopeni, Bucharest (Romania) and its close phylogenetic relationships to north Italy populations suggests human-mediated introductions of this wall lizard clade in Romania. Our results confirm that P.muralis can have non-native lineages and admixture through active human-mediated transport.

Intraspecific hybridisation of an invasive lizard on Lord Howe Island

Human-mediated dispersal of animals often acts to bring populations that have been separated for substantial periods of evolutionary time (e.g. millions of years) in their native range into contact in their introduced range. Whether these taxa successfully interbreed in the introduced range provides information on the strength of reproductive isolation amongst them. The invasive delicate skink (Lampropholis delicata) has been accidentally introduced to Lord Howe Island from four genetically divergent (>2 million years) regions of the species’ native range in eastern Australia. We used mitochondrial DNA and microsatellite data to investigate whether the individuals from four of the native-range source regions are interbreeding on Lord Howe Island. Our analyses indicate that intraspecific hybridisation among individuals from all four native-range source regions is occurring. Although there is little evidence for hybrids in the northern end of Lord Howe Island (proportion of hybrids: 0–0.02; n = 31), there is a high proportion of hybrids in the central (0.33–0.69; n = 59) and southern regions (0.38–0.75; n = 8) of the island. Given the strong evidence for interbreeding among all four native-range source regions examined, and the relatively minor morphological, life-history and phenotypic variation among them, we suggest that the delicate skink should continue to be treated as a single, widespread, but variable species.

Population genetic differentiation and genomic signatures of adaptation to climate in an abundant lizard

Species distributed across climatic gradients will typically experience spatial variation in selection, but gene flow can prevent such selection from causing population genetic differentiation and local adaptation. Here, we studied genomic variation of 415 individuals across 34 populations of the common wall lizard (Podarcis muralis) in central Italy. This species is highly abundant throughout this region and populations belong to a single genetic lineage, yet there is extensive phenotypic variation across climatic regimes. We used redundancy analysis to, first, quantify the effect of climate and geography on population genomic variation in this region and, second, to test if climate consistently sorts specific alleles across the landscape. Climate explained 5% of the population genomic variation across the landscape, about half of which was collinear with geography. Linear models and redundancy analyses identified loci that were significantly differentiated across climatic regimes. These loci were distributed across the genome and physically associated with genes putatively involved in thermal tolerance, regulation of temperature-dependent metabolism and reproductive activity, and body colouration. Together, these findings suggest that climate can exercise sufficient selection in lizards to promote genetic differentiation across the landscape in spite of high gene flow.

Asymmetric allelic introgression across a hybrid zone of the coal tit (Periparus ater) in the central Himalayas

In the Himalayas, a number of secondary contact zones have been described for vicariant vertebrate taxa. However, analyses of genetic divergence and admixture are missing for most of these examples. In this study, we provide a population genetic analysis for the coal tit (Periparus ater) hybrid zone in Nepal. Intermediate phenotypes between the distinctive western "spot-winged tit" (P. a. melanolophus) and Eastern Himalayan coal tits (P. a. aemodius) occur across a narrow range of <100 km in western Nepal. As a peculiarity, another distinctive cinnamon-bellied form is known from a single population so far. Genetic admixture of western and eastern mitochondrial lineages was restricted to the narrow zone of phenotypically intermediate populations. The cline width was estimated 46 km only with a center close to the population of the cinnamon-bellied phenotype. In contrast, allelic introgression of microsatellite loci was asymmetrical from eastern P. a. aemodius into far western populations of phenotypic P. a. melanolophus but not vice versa. Accordingly, the microsatellite cline was about 3.7 times wider than the mitochondrial one.

Population genomics of wall lizards reflects the dynamic history of the Mediterranean Basin

The Mediterranean Basin has experienced extensive change in geology and climate over the past six million years. Yet, the relative importance of key geological events for the distribution and genetic structure of the Mediterranean fauna remains poorly understood. Here, we use population genomic and phylogenomic analyses to establish the evolutionary history and genetic structure of common wall lizards (Podarcis muralis). This species is particularly informative because, in contrast to other Mediterranean lizards, it is widespread across the Iberian, Italian, and Balkan Peninsulas, and in extra-Mediterranean regions. We found strong support for six major lineages within P. muralis, which were largely discordant with the phylogenetic relationship of mitochondrial DNA. The most recent common ancestor of extant P. muralis was likely distributed in the Italian Peninsula, and experienced an “Out-of-Italy” expansion following the Messinian salinity crisis (∼5 Mya), resulting in the differentiation into the extant lineages on the Iberian, Italian, and Balkan Peninsulas. Introgression analysis revealed that both inter- and intraspecific gene flows have been pervasive throughout the evolutionary history of P. muralis. For example, the Southern Italy lineage has a hybrid origin, formed through admixture between the Central Italy lineage and an ancient lineage that was the sister to all other P. muralis. More recent genetic differentiation is associated with the onset of the Quaternary glaciations, which influenced population dynamics and genetic diversity of contemporary lineages. These results demonstrate the pervasive role of Mediterranean geology and climate for the evolutionary history and population genetic structure of extant species.

Climate Shapes the Geographic Distribution and Introgressive Spread of Color Ornamentation in Common Wall Lizards

AbstractClimate can exert an effect on the strength of sexual selection, but empirical evidence is limited. Here, we tested whether climate predicts the geographic distribution and introgressive spread of sexually selected male color ornamentation across 114 populations of the common wall lizard, Podarcis muralis. Coloration was highly structured across the landscape and did not reflect genetic differentiation. Instead, color ornamentation was consistently exaggerated in hot and dry environments, suggesting that climate-driven selection maintains geographic variation in spite of gene flow. Introgression of color ornamentation into a distantly related lineage appears to be ongoing and was particularly pronounced in warm climates with wet winters and dry summers. Combined, these results suggest that sexual ornamentation is consistently favored in climates that allow a prolonged reproductive season and high and reliable opportunities for lizard activity. This pattern corroborates theoretical predictions that such climatic conditions reduce the temporal clustering of receptive females and increase male-male competition, resulting in strong sexual selection. In summary, we provide compelling evidence for the importance of climate for the evolution of color ornamentation, and we demonstrate that geographic variation in the strength of sexual selection influences introgression of this phenotype.

Anthropogenic evolution in an insect wing polymorphism following widespread deforestation

Anthropogenic environmental change can underpin major shifts in natural selective regimes, and can thus alter the evolutionary trajectories of wild populations. However, little is known about the evolutionary impacts of deforestation-one of the most pervasive human-driven changes to terrestrial ecosystems globally. Absence of forest cover (i.e. exposure) has been suggested to play a role in selecting for insect flightlessness in montane ecosystems. Here, we capitalize on human-driven variation in alpine treeline elevation in New Zealand to test whether anthropogenic deforestation has caused shifts in the distributions of flight-capable and flightless phenotypes in a wing-polymorphic lineage of stoneflies from the Zelandoperla fenestrata species complex. Transect sampling revealed sharp transitions from flight-capable to flightless populations with increasing elevation. However, these phenotypic transitions were consistently delineated by the elevation of local treelines, rather than by absolute elevation, providing a novel example of human-driven evolution in response to recent deforestation. The inferred rapid shifts to flightlessness in newly deforested regions have implications for the evolution and conservation of invertebrate biodiversity.

Fossil-calibrated time tree of Podarcis wall lizards provides limited support for biogeographic calibration models

Podarcis wall lizards are endemic to the Mediterranean Basin where they represent the predominant reptile group. Despite being extensively used as model organisms in evolutionary and ecological studies their phylogeny and historical biogeography are still incompletely resolved. Moreover, molecular clock calibrations used in wall lizard phylogeography are based on the assumption of vicariant speciation triggered by the abrupt Mediterranean Sea level rise at the end of the Messinian salinity crisis (MSC). However, the validity of this biogeographic calibration remains untested. In this study we inferred a robust time tree based on multilocus data and fossil calibrations using both gene concatenation and species-tree approaches and including models with gene-flow. We found five deeply divergent, geographically coherent, and well-supported clades comprising species from i) Iberian Peninsula and North Africa; ii) Western Mediterranean islands, iii) Sicilian and Maltese islands; and iv-v) Balkan region and Aegean islands. The mitochondrial tree shows some inconsistencies with the species tree that warrant future investigation. Diversification of main clades is estimated in a short time frame during the Middle Miocene and might have been associated with a period of global climate cooling with the establishment of a marked climatic zonation in Europe. Cladogenetic events within the main clades are scattered throughout the time tree, from the Late Miocene to the Early Pleistocene, suggesting that speciation events in wall lizards reflect a complex interplay between regional topography, climate and geological history rather than a shared major climatic or paleogeographic event. Our absolute time estimates, as well as a relative dating approach, demonstrate that the assumption of a causal link between sea-level rise at the end of the MSC and the diversification of many island endemics is not justified. This study reinforces the notion that multiple dispersal and vicariant events, at different time frames, are required to explain current allopatric distributions and to account for the historical assembly of Mediterranean biota, and cautions against the use of biogeographic calibrations based on the assumption of vicariance.

Biased Hybridization and Its Impact on Adaptive Introgression

Gene exchange between species can influence ecological and evolutionary processes ranging from population rescue to adaptive radiation. Genomic tools have provided new insights into the prevalence and nature of gene exchange between species. However, much remains unknown of how ecological, behavioral, and evolutionary factors determine what genetic variation moves between species in the first place. In particular, more research is needed that evaluates whether such factors bias gene flow from one species to another, and whether any such biases affect how genetic variation from another species is ultimately retained in the genome of a given species. Addressing this issue is crucial in a changing world where hybridization and introgression might determine which species succeed and which become extinct.

Maternal reproductive output and F1 hybrid fitness may influence contact zone dynamics

The outcome of secondary contact between divergent lineages or species may be influenced by both the reproductive traits of parental species and the fitness of offspring; however, their relative contributions have rarely been evaluated, particularly in longer-lived vertebrate species. We performed pure and reciprocal laboratory crosses between Ctenophorus decresii (tawny dragon) and C. modestus (swift dragon) to examine how parental reproductive traits and ecologically relevant offspring fitness traits may explain contact zone dynamics in the wild. The two species meet in a contact zone of post-F1 hybrids with asymmetric backcrossing and predominantly C. modestus mtDNA haplotypes. We found no evidence for reduced parental fecundity or offspring fitness for F1 hybrid crosses. However, maternal reproductive strategy differed between species, irrespective of the species of their mate. Ctenophorus modestus females had higher fecundity and produced more and larger clutches with lower embryonic mortality. Parental species also influenced sex ratios and offspring traits, with C. modestus ♀ × C. decresii ♂ hybrids exhibiting higher trait values for more fitness measures (growth rate, sprint speed, bite force) than offspring from all other pairings. Together, these patterns are consistent with the prevalence of C. modestus mtDNA in the contact zone, and asymmetric backcrossing likely reflects fitness effects that manifest in the F2 generation. Our results highlight how parental species can influence multiple offspring traits in different ways, which together may combine to influence offspring fitness and shape contact zone dynamics.

Male–male contest limits the expression of assortative mate preferences in a polymorphic poison frog

Codivergence of sexual traits and mate preferences can lead to assortative mating and subsequently reproductive isolation. However, mate choice rarely operates without intrasexual competition, and the effects of the latter on speciation are often overlooked. Maintaining trait polymorphisms despite gene flow and limiting assortative female preferences for less-competitive male phenotypes are two important roles that male–male competition may play in the speciation process. Both roles rely on the assumption that male–male competition limits the expression of divergent female preferences. We tested this assumption in the highly color-polymorphic strawberry poison frog (Oophaga pumilio). Females prefer males of the local color, suggesting that reproductive isolation may be evolving among color morphs. However, this inference does not account for male–male competition, which is also color-mediated. We housed females with two differently colored males, and compared reproductive patterns when the more attractive male was the territory holder versus when he was the nonterritorial male. Females mated primarily with the territory winner, regardless of coloration, suggesting that when a choice must be made between the two, male territoriality overrides female preferences for male coloration. Our results highlight the importance of considering the combined effects of mate choice and intrasexual competition in shaping phenotypic divergence and speciation.

Genetic and phenotypic evidence of a contact zone between divergent colour morphs of the iconic red-eyed treefrog

Hybrid zones act as natural laboratories where divergent genomes interact, providing powerful systems for examining the evolutionary processes underlying biological diversity. In this study, we characterized patterns of genomic and phenotypic variation resulting from hybridization between divergent intraspecific lineages of the Neotropical red-eyed treefrog (Agalychnis callidryas). We found genetic evidence of a newly discovered contact zone and phenotypic novelty in leg colour-a trait suspected to play a role in mediating assortative mating in this species. Analysis of hybrid ancestry revealed an abundance of later-generation F_n individuals, suggesting persistence of hybrids in the contact zone. Hybrids are predominantly of southern ancestry but are phenotypically more similar to northern populations. Genome-wide association mapping revealed QTL with measurable effects on leg-colour variation, but further work is required to dissect the architecture of this trait and establish causal links. Further, genomic cline analyses indicated substantial variation in patterns of introgression across the genome. Directional introgression of loci associated with different aspects of leg colour are inherited from each parental lineage, creating a distinct hybrid colour pattern. We show that hybridization can generate new phenotypes, revealing the evolutionary processes that potentially underlie patterns of phenotypic diversity in this iconic polytypic frog. Our study is consistent with a role of hybridization and sexual selection in lineage diversification, evolutionary processes that have been implicated in accelerating divergence in the most phenotypically diverse species.

No evidence for differential sociosexual behavior and space use in the color morphs of the European common wall lizard (Podarcis muralis)

Explaining the evolutionary origin and maintenance of color polymorphisms is a major challenge in evolutionary biology. Such polymorphisms are commonly thought to reflect the existence of alternative behavioral or life-history strategies under negative frequency-dependent selection. The European common wall lizard Podarcis muralis exhibits a striking ventral color polymorphism that has been intensely studied and is often assumed to reflect alternative reproductive strategies, similar to the iconic "rock-paper-scissors" system described in the North American lizard Uta stansburiana. However, available studies so far have ignored central aspects in the behavioral ecology of this species that are crucial to assess the existence of alternative reproductive strategies. Here, we try to fill this gap by studying the social behavior, space use, and reproductive performance of lizards showing different color morphs, both in a free-ranging population from the eastern Pyrenees and in ten experimental mesocosm enclosures. In the natural population, we found no differences between morphs in site fidelity, space use, or male-female spatial overlap. Likewise, color morph was irrelevant to sociosexual behavior, space use, and reproductive success within experimental enclosures. Our results contradict the commonly held hypothesis that P. muralis morphs reflect alternative behavioral strategies, and suggest that we should instead turn our attention to alternative functional explanations.

The Spatial Signature of Introgression After a Biological Invasion With Hybridization

The accumulation of genome-wide molecular data has emphasized the important role of hybridization in the evolution of many organisms, which may carry introgressed genomic segments resulting from past admixture events with other taxa. Despite a number of examples of hybridization occurring during biological invasions, the resulting spatial patterns of genomic introgression remain poorly understood. Preliminary simulation studies have suggested a heterogeneous spatial level of introgression for invasive taxa after range expansion. We investigated in detail the robustness of this pattern and its persistence over time for both invasive and local organisms. Using spatially explicit simulations, we explored the spatial distribution of introgression across the area of colonization of an invasive taxon hybridizing with a local taxon. The general pattern for neutral loci supported by our results is an increasing introgression of local genes into the invasive taxon with the increase in the distance from the source of the invasion and a decreasing introgression of invasive genes into the local taxon. However, we also show there is some variation in this general trend depending on the scenario investigated. Spatial heterogeneity of introgression within a given taxon is thus an expected neutral pattern in structured populations after a biological invasion with a low to moderate amount of hybridization. We further show that this pattern is consistent with published empirical observations. Using additional simulations, we argue that the spatial pattern of Neanderthal introgression in modern humans, which has been documented to be higher in Asia than in Europe, can be explained by a model of hybridization with Neanderthals in Eurasia during the range expansion of modern humans from Africa. Our results support the view that weak hybridization during range expansion may explain spatially heterogeneous introgression patterns without the need to invoke selection.

Digest: Biased male-male competition drives asymmetric introgression in lizards

Introgression, gene flow from one population into another, can be asymmetric. Yang et al. suggest that reduction of gene flow in one direction, rather than elevated gene flow in the opposite direction, explains the pattern of asymmetric introgression between two lizard lineages. The authors propose that a dominant male phenotype in one lineage blocks a submissive male phenotype from another lineage in mating with females of the opposite lineage. This case underscores just how capricious introgression can be.

Spatial variation in gene flow across a hybrid zone reveals causes of reproductive isolation and asymmetric introgression in wall lizards

Hybrid zones provide insights into the evolution of reproductive isolation. Sexual selection can contribute to the evolution of reproductive barriers, but it remains poorly understood how sexual traits impact gene flow in secondary contact. Here, we show that a recently evolved suite of sexual traits that function in male-male competition mediates gene flow between two lineages of wall lizards (Podarcis muralis). Gene flow was relatively low and asymmetric in the presence of exaggerated male morphology and coloration compared to when the lineages share the ancestral phenotype. Putative barrier loci were enriched in genomic regions that were highly differentiated between the two lineages and showed low concordance between the transects. The exception was a consistently low genetic exchange around ATXN1, a gene that modulates social behavior. We suggest that this gene may contribute to the male mate preferences that are known to cause lineage-assortative mating in this species. Although female choice modulates the degree of reproductive isolation in a variety of taxa, wall lizards demonstrate that both male-male competition and male mate choice can contribute to the extent of gene flow between lineages.

Dispersal Predicts Hybrid Zone Widths across Animal Diversity: Implications for Species Borders under Incomplete Reproductive Isolation

Hybrid zones occur as range boundaries for many animal taxa. One model for how hybrid zones form and stabilize is the tension zone model, a version of which predicts that hybrid zone widths are determined by a balance between random dispersal into hybrid zones and selection against hybrids. Here, we examine whether random dispersal and proxies for selection against hybrids (genetic distances between hybridizing pairs) can explain variation in hybrid zone widths across 131 hybridizing pairs of animals. We show that these factors alone can explain ∼40% of the variation in zone width among animal hybrid zones, with dispersal explaining far more of the variation than genetic distances. Patterns within clades were idiosyncratic. Genetic distances predicted hybrid zone widths particularly well for reptiles, while this relationship was opposite tension zone predictions in birds. Last, the data suggest that dispersal and molecular divergence set lower bounds on hybrid zone widths in animals, indicating that there are geographic restrictions on hybrid zone formation. Overall, our analyses reinforce the fundamental importance of dispersal in hybrid zone formation and more generally in the ecology of range boundaries.

Divergent male and female mate preferences do not explain incipient speciation between lizard lineages

Diversification in sexual signals is often taken as evidence for the importance of sexual selection in speciation. However, in order for sexual selection to generate reproductive isolation between populations, both signals and mate preferences must diverge together. Furthermore, assortative mating may result from multiple behavioral mechanisms, including female mate preferences, male mate preferences, and male–male competition; yet their relative contributions are rarely evaluated. Here, we explored the role of mate preferences and male competitive ability as potential barriers to gene flow between 2 divergent lineages of the tawny dragon lizard, Ctenophorus decresii, which differ in male throat coloration. We found stronger behavioral barriers to pairings between southern lineage males and northern lineage females than between northern males and southern females, indicating incomplete and asymmetric behavioral isolating barriers. These results were driven by both male and female mate preferences rather than lineage differences in male competitive ability. Intrasexual selection is therefore unlikely to drive the outcome of secondary contact in C. decresii, despite its widely acknowledged importance in lizards. Our results are consistent with the emerging view that although both male and female mate preferences can diverge alongside sexual signals, speciation is rarely driven by divergent sexual selection alone.

Testing an hypothesis of hybrid zone movement for toads in France

Hybrid zone movement may result in substantial unidirectional introgression of selectively neutral material from the local to the advancing species, leaving a genetic footprint. This genetic footprint is represented by a trail of asymmetric tails and displaced cline centres in the wake of the moving hybrid zone. A peak of admixture linkage disequilibrium is predicted to exist ahead of the centre of the moving hybrid zone. We test these predictions of the movement hypothesis in a hybrid zone between common (Bufo bufo) and spined toads (B. spinosus), using 31 nuclear and one mtDNA SNPs along a transect in the northwest of France. Average effective selection in Bufo hybrids is low and clines vary in shape and centre. A weak pattern of asymmetric introgression is inferred from cline discordance of seven nuclear markers. The dominant direction of gene flow is from B. spinosus to B. bufo and is in support of southward movement of the hybrid zone. Conversely, a peak of admixture linkage disequilibrium north of the hybrid zone suggests northward movement. These contrasting results can be explained by reproductive isolation of the B. spinosus and B. bufo gene pools at the southern (B. spinosus) side of the hybrid zone. The joint occurrence of asymmetric introgression and admixture linkage disequilibrium can also be explained by the combination of low dispersal and random genetic drift due to low effective population sizes.

Past lake shore dynamics explain present pattern of unidirectional introgression across a habitat barrier

Introgression patterns between divergent lineages are often characterized by asymmetry in the direction and among-marker variation in the extent of gene flow, and therefore inform on the mechanisms involved in differentiation and speciation. In the present study, we test the hypothesis that unidirectional introgression between two phenotypically and genetically distinct lineages of the littoral, rock-dwelling cichlid fish Tropheus moorii across a wide sandy bay is linked to observed differences in mate preferences between the two lineages. This hypothesis predicts bi-directional nuclear gene flow and was rejected by congruent patterns of introgression in mtDNA, AFLP and microsatellite markers, with admixture confined to the populations west of the bay. This pattern can be explained on the basis of habitat changes in the course of lake level fluctuations, which first facilitated the development of a symmetric admixture zone including the area corresponding to the present sand bay and then shaped asymmetry by causing local extinctions and cessation of gene flow when this area became once more inhabitable. This conforms with previous assumptions that habitat dynamics are a primary determinant of population-level evolution in Tropheus. In this respect, Tropheus may be representative of species whose preferred habitat is subject to frequent re-structuring.

Assessing Behavioral Associations in a Hybrid Zone through Social Network Analysis: Complex Assortative Behaviors Structure Associations in a Hybrid Quail Population

Behavior can strongly influence rates and patterns of hybridization between animal populations and species. Yet few studies have examined reproductive behaviors in natural hybrid zones within the fine-scale social context in which they naturally occur. We use radio-frequency identification tags with social network analyses to test whether phenotypic similarity in plumage and mass correlate with social behavior throughout a breeding season in a California and Gambel's quail hybrid zone. We use a novel approach to partition phenotypic variation in a way that does not confound differences between sexes and species, and we illustrate the complex ways that phenotype and behavior structure the social environment, mating opportunities, and male-male associations. Associations within the admixed population were random with respect to species-specific plumage but showed strong patterns of assortment based on sexually dimorphic plumage, monomorphic plumage, and mass. Weak behavioral reproductive isolation in this admixed population may be the result of complex patterns of phenotypic assortment based on multiple traits rather than a lack of phenotypic discrimination. More generally, our results support the utility of social network analyses for analyzing behavioral factors affecting genetic exchange between populations and species.

Hybridization and introgression during density-dependent range expansion: European wildcats as a case study

Interbreeding between historically allopatric species with incomplete reproductive barriers may result when species expand their range. The genetic consequences of such hybridization depend critically on the dynamics of the range expansion. Hybridization models during range expansion have been developed but assume dispersal to be independent from neighboring population densities. However, organisms may disperse because they are attracted by conspecifics or because they prefer depopulated areas. Here, through spatially explicit simulations, we assess the effect of various density-dependent dispersal modes on the introgression between two species. We find huge introgression from the local species into the invasive one with all dispersal modes investigated, even when the hybridization rate is relatively low. This represents a general expectation for neutral genes even if the dispersal modes differ in colonization times and amount of introgression. Invasive individuals attracted by conspecifics need more time to colonize the whole area and are more introgressed by local genes, whereas the opposite is found for solitary individuals. We applied our approach to a recent expansion of European wildcats in the Jura Mountains and the hybridization with domestic cats. We show that the simulations explained better the observed level of introgression at nuclear, mtDNA, and Y chromosome markers, when using solitary dispersal for wildcats instead of random or gregarious dispersal, in accordance with ecological knowledge. Thus, use of density-dependent dispersal models increases the predictive power of the approach.

Differential introgression of a female competitive trait in a hybrid zone between sex-role reversed species

Mating behavior between recently diverged species in secondary contact can impede or promote reproductive isolation. Traditionally, researchers focus on the importance of female mate choice and male-male competition in maintaining or eroding species barriers. Although female-female competition is widespread, little is known about its role in the speciation process. Here, we investigate a case of interspecific female competition and its influence on patterns of phenotypic and genetic introgression between species. We examine a hybrid zone between sex-role reversed, Neotropical shorebird species, the northern jacana (Jacana spinosa) and wattled jacana (J. jacana), in which female-female competition is a major determinant of reproductive success. Previous work found that females of the more aggressive and larger species, J. spinosa, disproportionately mother hybrid offspring, potentially by monopolizing breeding territories in sympatry with J. jacana. We find a cline shift of female body mass relative to the genetic center of the hybrid zone, consistent with asymmetric introgression of this competitive trait. We suggest that divergence in sexual characteristics between sex-role reversed females can influence patterns of gene flow upon secondary contact, similar to males in systems with more typical sex roles.

Maladaptive secondary sexual characteristics reduce the reproductive success of hybrids between native and non-native salmonids

Human-mediated hybridization between introduced and native species is one of the most serious threats to native taxa. Although field studies have attempted to quantify the relative fitness or reproductive success of parental species and their hybrids, only a few studies have unraveled the factors determining the fitness of hybrids. Here, we hypothesized that maladaptive secondary sexual characteristics may reduce fitness of hybrids between two fish species. To test this, we evaluated the reproductive success of introduced brook trout (BT: Salvelinus fontinalis), native white-spotted charr (WSC: S. leucomaenis) and their hybrids in a natural stream in Hokkaido, Japan, where the two parental species show remarkably different male secondary sexual characteristics, such as elongated jaws and deeper bodies. We predicted that introgression from WSC is maladaptive for BT males because the BT male has more prominent secondary sexual characteristics. Our results suggest that both sexual selection and outbreeding depression in males and females significantly influence an individual's reproductive success. Our results also suggest that asymmetric introgression may increase the risks to persistence in the recipient species.

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.

Admixture of hybrid swarms of native and introduced lizards in cities is determined by the cityscape structure and invasion history

Introductions of non-native lineages increase opportunities for hybridization. Non-native lineages of the common wall lizard, Podarcis muralis, are frequently introduced in cities where they hybridize with native populations. We aimed at unravelling the invasion history and admixture of native and non-native wall lizards in four German cities using citywide, comprehensive sampling. We barcoded and genotyped 826 lizards and tested if gene flow in populations composed of admixed native and introduced lineages is facilitated by similar environmental factors to those in native populations by comparing fine-scale landscape genetic patterns. In cities with non-native lineages, lizards commonly occurred in numerous clusters of hybrid swarms, which showed variable lineage composition, consisting of up to four distinct evolutionary lineages. Hybrid swarms held vast genetic diversity and showed recent admixture with other hybrid swarms. Landscape genetic analyses showed differential effects of cityscape structures across cities, but identified water bodies as strong barriers to gene flow in both native and admixed populations. By contrast, railway tracks facilitated gene flow of admixed populations only. Our study shows that cities represent unique settings for hybridization, caused by multiple introductions of non-native taxa. Cityscape structure and invasion histories of cities will determine future evolutionary pathways at these novel hybrid zones.

Interspecific competition, hybridization, and reproductive isolation in secondary contact: missing perspectives on males and females

Research on sexual selection and hybridization has focused on female mate choice and male-male competition. While the evolutionary outcomes of interspecific female preference have been well explored, we are now gaining a better understanding of the processes by which male-male competition between species in secondary contact promotes reproductive isolation versus hybridization. What is relatively unexplored is the interaction between female choice and male competition, as they can oppose one another or align with similar outcomes for reproductive isolation. The role of female-female competition in hybridization is also not well understood, but could operate similarly to male-male competition in polyandrous and other systems where costs to heterospecific mating are low for females. Reproductive competition between either sex of sympatric species can cause the divergence and/or convergence of sexual signals and recognition, which in turn influences the likelihood for interspecific mating. Future work on species interactions in secondary contact should test the relative influences of both mate choice and competition for mates on hybridization outcomes, and should not ignore the possibilities that females can compete over mating resources, and males can exercise mate choice.

Color and behavior differently predict competitive outcomes for divergent stickleback color morphs

Our knowledge of how male competition contributes to speciation is dominated by investigations of competition between within-species morphs or closely related species that differ in conspicuous traits expressed during the breeding season (e.g. color, song). In such studies, it is important to consider the manner in which putatively sexually selected traits influence the outcome of competitive interactions within and between types because these traits can communicate information about competitor quality and may not be utilized by homotypic and heterotypic receivers in the same way. We studied the roles of breeding color and aggressive behaviors in competition within and between two divergent threespine stickleback Gasterosteus aculeatus color types. Our previous work in this system showed that the switch from red to black breeding coloration is associated with changes in male competition biases. Here, we find that red and black males also use different currencies in competition. Winners of both color types performed more aggressive behaviors than losers, regardless of whether the competitor was of the same or opposite color type. But breeding color differently predicted competitive outcomes for red and black males. Males who were redder at the start of competition were more likely to win when paired with homotypic competitors and less likely to win when paired with heterotypic competitors. In contrast, black color, though expressed in the breeding season and condition dependent, was unrelated to competitive outcomes. Placing questions about the role of male competition in speciation in a sexual signal evolution framework may provide insight into the "why and how" of aggression biases and asymmetries in competitive ability between closely related morphs and species.

The Role of Sexual Selection in Local Adaptation and Speciation

Sexual selection plays several intricate and complex roles in the related processes of local adaptation and speciation. In some cases sexual selection can promote these processes, but in others it can be inhibitory. We present theoretical and empirical evidence supporting these dual effects of sexual selection during local adaptation, allopatric speciation, and speciation with gene flow. Much of the empirical evidence for sexual selection promoting speciation is suggestive rather than conclusive; we present what would constitute strong evidence for sexual selection driving speciation. We conclude that although there is ample evidence that sexual selection contributes to the speciation process, it is very likely to do so only in concert with natural selection.

Females drive asymmetrical introgression from rare to common species in Darwin's tree finches

The consequences of hybridization for biodiversity depend on the specific ecological and evolutionary context in which it occurs. Understanding patterns of gene flow among hybridizing species is crucial for determining the evolutionary trajectories of species assemblages. The recently discovered hybridization between two species of Darwin's tree finches (Camarhynchus parvulus and C. pauper) on Floreana Island, Galápagos, presents an exciting opportunity to investigate the mechanisms causing hybridization and its potential evolutionary consequences under conditions of recent habitat disturbance and the introduction of invasive pathogens. In this study, we combine morphological and genetic analysis with pairing observations to explore the extent, direction and drivers of hybridization and to test whether hybridization patterns are a result of asymmetrical pairing preference driven by females of the rarer species (C. pauper). We found asymmetrical introgression from the critically endangered, larger-bodied C. pauper to the common, smaller-bodied C. parvulus, which was associated with a lack of selection against heterospecific males by C. pauper females. Examination of pairing data showed that C. parvulus females paired assortatively, whereas C. pauper females showed no such pattern. This study shows how sex-specific drivers can determine the direction of gene flow in hybridizing species. Furthermore, our results suggest the existence of a hybrid swarm comprised of C. parvulus and hybrid birds. We discuss the influence of interspecific abundance differences and susceptibility to the invasive parasite Philornis downsi on the observed hybridization and recommend that the conservation of this iconic species group should be managed jointly rather than species-specific.