Local adaptation within a hybrid species

Ecological divergence among populations may be strongly influenced by their genetic background. For instance, genetic admixture through introgressive hybridization or hybrid speciation is likely to affect the genetic variation and evolvability of phenotypic traits. We studied geographic variation in two beak dimensions and three other phenotypic traits of the Italian sparrow (Passer italiae), a young hybrid species formed through interbreeding between house sparrows (P. domesticus) and Spanish sparrows (P. hispaniolensis). We found that beak morphology was strongly influenced by precipitation regimes and that it appeared to be the target of divergent selection within Italian sparrows. Interestingly, however, the degree of parental genetic contribution in the hybrid species had no effect on phenotypic beak variation. Moreover, beak height divergence may mediate genetic differentiation between populations, consistent with isolation-by-adaptation within this hybrid species. The study illustrates how hybrid species may be relatively unconstrained by their admixed genetic background, allowing them to adapt rapidly to environmental variation.     

North African hybrid sparrows (Passer domesticus,P. hispaniolensis) back from oblivion – ecological segregation and asymmetric mitochondrial introgression between parental species

A stabilized hybrid form of the house sparrow (Passer domesticus) and the Spanish sparrow (P. hispaniolensis) is known as Passer italiae from the Italian Peninsula and a few Mediterranean islands. The growing attention for the Italian hybrid sparrow and increasing knowledge on its biology and genetic constitution greatly contrast the complete lack of knowledge of the long‐known phenotypical hybrid sparrow populations from North Africa. Our study provides new data on the breeding biology and variation of mitochondrial DNA in three Algerian populations of house sparrows, Spanish sparrows, and phenotypical hybrids. In two field seasons, the two species occupied different breeding habitats: Spanish sparrows were only found in rural areas outside the cities and bred in open‐cup nests built in large jujube bushes. In contrast, house sparrows bred only in the town centers and occupied nesting holes in walls of buildings. Phenotypical hybrids were always associated with house sparrow populations. House sparrows and phenotypical hybrids started breeding mid of March, and most pairs had three successive clutches, whereas Spanish sparrows started breeding almost one month later and had only two successive clutches. Mitochondrial introgression is strongly asymmetric because about 75% of the rural Spanish sparrow population carried house sparrow haplotypes. In contrast, populations of the Italian hybrid form, P. italiae, were genetically least diverse among all study populations and showed a near‐fixation of house sparrow haplotypes that elsewhere were extremely rare or that were even unique for the Italian Peninsula. Such differences between mitochondrial gene pools of Italian and North African hybrid sparrow populations provide first evidence that different demographic histories have shaped the extant genetic diversity observed on both continents.     

Genetic admixture despite ecological segregation in a North African sparrow hybrid zone (Aves,Passeriformes, Passer domesticus × Passer hispaniolensis)

Under different environmental conditions, hybridization between the same species might result in different patterns of genetic admixture. Particularly, species pairs with large distribution ranges and long evolutionary history may have experienced several independent hybridization events over time in different zones of overlap. In birds, the diverse hybrid populations of the house sparrow (Passer domesticus) and the Spanish sparrow (Passer hispaniolensis) provide a striking example. Throughout their range of sympatry, these two species do not regularly interbreed; however, a stabilized hybrid form (Passer italiae) exists on the Italian Peninsula and on several Mediterranean islands. The spatial distribution pattern on the Eurasian continent strongly contrasts the situation in North Africa, where house sparrows and Spanish sparrows occur in close vicinity of phenotypically intermediate populations across a broad mosaic hybrid zone. In this study, we investigate patterns of divergence and admixture among the two parental species, stabilized and nonstabilized hybrid populations in Italy and Algeria based on a mitochondrial marker, a sex chromosomal marker, and 12 microsatellite loci. In Algeria, despite strong spatial and temporal separation of urban early‐breeding house sparrows and hybrids and rural late‐breeding Spanish sparrows, we found strong genetic admixture of mitochondrial and nuclear markers across all study populations and phenotypes. That pattern of admixture in the North African hybrid zone is strikingly different from i) the Iberian area of sympatry where we observed only weak asymmetrical introgression of Spanish sparrow nuclear alleles into local house sparrow populations and ii) the very homogenous Italian sparrow population where the mitogenome of one parent (P. domesticus) and the Z‐chromosomal marker of the other parent (P. hispaniolensis) are fixed. The North African sparrow hybrids provide a further example of enhanced hybridization along with recent urbanization and anthropogenic land‐use changes in a mosaic landscape.     

Hybrid speciation in sparrows II: a role for sex chromosomes?

Homoploid hybrid speciation in animals is poorly understood, mainly because of the scarcity of well‐documented cases. Here, we present the results of a multilocus sequence analysis on the house sparrow (Passer domesticus), Spanish sparrow (P. hispaniolensis) and their proposed hybrid descendant, the Italian sparrow (P. italiae). The Italian sparrow is shown to be genetically intermediate between the house sparrow and Spanish sparrow, exhibiting genealogical discordance and a mosaic pattern of alleles derived from either of the putative parental species. The average variation on the Z chromosome was significantly reduced compared with autosomal variation in the putative parental species, the house sparrow and Spanish sparrow. Additionally, divergence between the two species was elevated on the Z chromosome relative to the autosomes. This pattern of variation and divergence is consistent with reduced introgression of Z‐linked genes and/or a faster‐Z effect (increased rate of adaptive divergence on the Z). F_ST‐outlier tests were consistent with the faster‐Z hypothesis: two of five Z‐linked loci (CHD1Z and PLAA) were identified as candidates for being subject to positive, divergent selection in the putative parental species. Interestingly, the two latter genes showed a mosaic pattern in the (hybrid) Italian sparrow; that is, the Italian sparrow was found to be fixed for Spanish sparrow alleles at CHD1Z and to mainly have house sparrow alleles at PLAA. Preliminary evidence presented in this study thus suggests that sex chromosomes may play a significant role in this case of homoploid hybrid speciation.     

Hybrid speciation in sparrows I: phenotypic intermediacy,genetic admixture and barriers to gene flow

Homoploid hybrid speciation is thought to require unusual circumstances to yield reproductive isolation from the parental species, and few examples are known from nature. Here, we present genetic evidence for this mode of speciation in birds. Using Bayesian assignment analyses of 751 individuals genotyped for 14 unlinked, nuclear microsatellite loci, we show that the phenotypically intermediate Italian sparrow (Passer italiae) does not form a cluster of its own, but instead exhibits clear admixture (over its entire breeding range) between its putative parental species, the house sparrow (P. domesticus) and the Spanish sparrow (P. hispaniolensis). Further, the Italian sparrow possesses mitochondrial (mt) DNA haplotypes identical to both putative parental species (although mostly of house sparrow type), indicating a recent hybrid origin. Today, the Italian sparrow has a largely allopatric distribution on the Italian peninsula and some Mediterranean islands separated from its suggested parental species by the Alps and the Mediterranean Sea, but co‐occurs with the Spanish sparrow on the Gargano peninsula in southeast Italy. No evidence of interbreeding was found in this sympatric population. However, the Italian sparrow hybridizes with the house sparrow in a sparsely populated contact zone in the Alps. Yet, the contact zone is characterized by steep clines in species‐specific male plumage traits, suggesting that partial reproductive isolation may also have developed between these two taxa. Thus, geographic and reproductive barriers restrict gene flow into the nascent hybrid species. We propose that an origin of hybrid species where the hybrid lineage gets geographically isolated from its parental species, as seems to have happened in this system, might be more common in nature than previously assumed.     

Population divergence and morphometric integration in the greenfinch (Carduelis chloris) – evolution against the trajectory of least resistance?

Theory predicts that genetic and phenotypic correlations among traits may direct the process of short-term evolution by limiting the directions of variation available to natural selection to act on. We studied correlations between 14 skeletal traits in 10 geographically distinct and relatively young greenfinch (Carduelis chloris) populations to unravel whether the divergence among populations has occurred into directions predicted by the within-population correlations (cf. drift/correlated responses models), or whether it is better explained by ‘adaptive’ models, which predict no necessary association between within- and among-population correlations (allometries). We found that the within-population character correlations (or covariances) did not predict character divergence between populations. This was because the first eigenvector of the among-population correlation/covariance matrix, summarizing the major dimension of divergence, was a bipolar body:beak dimension, and distinct from the (≈ isometric) first eigenvector of within-population matrix. Hence, as the divergence among greenfinch populations cannot be satisfactorily accommodated by drift/correlated response models, an adaptive basis for divergence is suggested. The second major axis of within-population variation was a classical ‘group size’ factor revealing that beak size was more or less free to vary independently of body size. Consequently, even if the divergence among populations cannot be simply accommodated to expectations of drift and correlated response models, it is striking that the most pronounced size-independent (nonallometric) changes had occurred along the second largest dimension of variance. This could mean that selection pressures which shape integration within populations are the same as those that cause divergence among populations. A relaxed beak:body integration could also occur as a result of species level selection favouring taxa in which independent evolution of beak and body is made possible.     

Morphological Variation in Wild Marmosets (Callithrix penicillata and C. geoffroyi) and Their Hybrids

Evolutionary theory and observation predict wider phenotypic variation in hybrids than parental species. Emergent phenotypic novelty in hybrids may in turn drive new adaptations or speciation by breaking parental phenotypic constraints. Primate hybridization is often documented through genetic evidence, but knowledge about the primate hybrid phenotype remains limited due to a small number of available studies on hybrid primate morphology. Here, we examine pelage and morphometric variation in two Brazilian marmoset species (Callithrix penicillata and C. geoffroyi) and their hybrids. Hybrids were sampled in an anthropogenic hybrid zone in the municipality of Viçosa, Minas Gerais state, Brazil. We analyzed hybrid facial and body pelage color variation, and compared 13 morphometric measures between hybrids and parental species. Five different hybrid facial morphotypes were observed, varying from intermediate to parental-like. Hybrid facial morphotypes were biased towards C. penicillata, suggesting that the pelage of this species may be dominant to that of C. geoffroyi in this context, and indicating that mate preference, and therefore gene flow/introgression, may be biased towards C. penicillata within the hybrid zone. Hybrid morphometric features were on average intermediate to parental species traits, but transgressive hybrids were also observed, suggesting that morphometric variation for the studied traits is consistent with Rieseberg’s complementary allele model. Finally, we observed a decoupling of facial patterning and size/shape in hybrids, relative to parent phenotypes, suggesting that an important factor driving phenotypic novelty within the Viçosa marmoset hybrid zone might be the loosening of evolutionary constraints on phenotypic trait integration.     

Fourteen polymorphic microsatellite loci characterized in the house sparrow Passer domesticus (Passeridae,Aves)

We characterized 14 polymorphic microsatellite loci in the house sparrow Passer domesticus. Four loci were isolated from house sparrow genomic libraries and 10 loci were identified by testing 100 loci that had been originally isolated in other passerine species. Loci were characterized in 37–54 unrelated sparrows from British and Norwegian populations. Each locus displayed between two and 31 alleles, with the observed heterozygosity ranging between 0.30 and 0.91.     

Phenotypic Changes Across a Geographic Gradient: The Case of Three Sympatric Dolphin Species

Phenotypic changes in the mammalian mandible can occur at different spatial and temporal scales. We investigated mandibular size and shape variation in three extant closely related dolphins (Cetacea, Odontoceti): Tursiops truncatus, Stenella coeruleoalba and Delphinus delphis in order to test the hypothesis that similar phenotypic changes occur across the same geographical gradient. Our data included 219 specimens representative of the following geographic locations: the Mediterranean Sea, the eastern north Atlantic and the North Sea. Each mandibula was photographed laterally and spatial positioning of eight homologous 2D landmarks was recorded. After applying generalised Procrustes analysis (GPA), intraspecific variation was first investigated between sexes and among populations to allow further pooling of samples. Size and shape differences among populations and species were investigated through multivariate ordination techniques (PCA), Procrustes ANOVA and allometric analyses. In all three species, Mediterranean populations clearly differed in mandible shape from the extra-Mediterranean ones. Among the three, the direction of geographic phenotypic changes was significantly similar in the striped and common dolphin, while the bottlenose dolphin was the most divergent species, differing both in size and allometric trajectory. Shape variation of the two former species highlighted a morphological convergence in the Atlantic, and a phenotypic divergence in the Mediterranean. Shape differences among the three dolphin species were interpreted in the light of different prey preferences, feeding strategies and habitat partitioning to avoid direct competition.     

Predictors of genomic differentiation within a hybrid taxon

Hybridization is increasingly recognized as an important evolutionary force. Novel genetic methods now enable us to address how the genomes of parental species are combined in hybrid lineages. However, we still do not know the relative importance of admixed proportions, genome architecture and local selection in shaping hybrid genomes. Here, we take advantage of the genetically divergent island populations of Italian sparrow on Crete, Corsica and Sicily to investigate the predictors of genomic variation within a hybrid taxon. We test if differentiation is affected by recombination rate, selection, or variation in ancestry proportions. We find that the relationship between recombination rate and differentiation is less pronounced within hybrid lineages than between the parent species, as expected if purging of minor parent ancestry in low recombination regions reduces the variation available for differentiation. In addition, we find that differentiation between islands is correlated with differences in signatures of selection in two out of three comparisons. Signatures of selection within islands are correlated across all islands, suggesting that shared selection may mould genomic differentiation. The best predictor of strong differentiation within islands is the degree of differentiation from house sparrow, and hence loci with Spanish sparrow ancestry may vary more freely. Jointly, this suggests that constraints and selection interact in shaping the genomic landscape of differentiation in this hybrid species.     

Does divergence from normal patterns of integration increase as chromosomal fusions increase in number? A test on a house mouse hybrid zone

Chromosomal evolution is widely considered an important driver of speciation because it can promote the establishment of reproductive barriers. Karyotypic reorganization is also expected to affect the mean phenotype, as well as its development and patterns of phenotypic integration, through processes such as variation in genetic linkage between quantitative trait loci or between regulatory regions and their targets. Here we explore the relationship between chromosomal evolution and phenotypic integration by analyzing a well-known house mouse parapatric contact zone between a highly derived Robertsonian (Rb) race (2n = 22) and populations with standard karyotype (2n = 40). Populations with hybrid karyotypes are scattered throughout the hybrid zone connecting the two parental races. Using mandible shape data and geometric morphometrics, we test the hypothesis that patterns of integration progressively diverge from the “normal” integration pattern observed in the standard race as they accumulate Rb fusions. We find that the main pattern of integration observed between the posterior and anterior part of the mandible can be largely attributed to allometry. We find no support for a gradual increase in divergence from normal patterns of integration as fusions accumulate. Surprisingly, however, we find that the derived Rb race (2n = 22) has a distinct allometric trajectory compared with the standard race. Our results suggest that either individual fusions disproportionately affect patterns of integration or that there are mechanisms which “purge” extreme variants in hybrids (e.g. reduced fitness of hybrid shape).     

Contrasting adaptive immune defenses and blood parasite prevalence in closely related Passer sparrows

Immune system components differ in their functions and costs, and immune defense profiles are likely to vary among species with differing ecologies. We compared adaptive immune defenses in two closely related species that have contrasting inflammatory immune responses, the widespread and abundant house sparrow (Passer domesticus) and the less abundant tree sparrow (Passer montanus). We found that the house sparrow, which we have previously shown mounts weaker inflammatory responses, exhibits stronger adaptive immune defenses, including antibody responses, natural antibody titers, and specific T-cell memory, than the tree sparrow. Conversely, tree sparrows, which mount strong inflammatory responses, also mount stronger nonspecific inflammatory T-cell responses but weaker specific adaptive responses. Prevalence of avian malaria parasite infections, which are controlled by adaptive immune defenses, was higher in the geographically restricted tree sparrow than in the ubiquitous house sparrow. Together these data describe distinct immune defense profiles between two closely related species that differ greatly in numbers and distributions. We suggest that these immunological differences could affect fitness in ways that contribute to the contrasting abundances of the two species in North American and Western Europe.     

Genetic and phenotypic correlations in a natural population of song sparrows

We estimated heritabilities, and genetic and phenotypic correlations between beak and body traits in the song sparrow ( Melospiza melodia ). We compared these estimates to values for the same traits in the Galápagos finches, Geospiza (Boag, 1983; Grant, 1983). Morphological variance is low in the song sparrow, and our results show that genetic and phenotypic correlations are considerably lower than correlations in the morphologically more variable Geospiza. Comparison using a larger sample of Galapagos populations confirms the existence of an association between variance and correlation for phenotypic values. We suggest two possible explanations for this association. First, most traits studied are functionally related, and the joint evolution of variance and correlation may have resulted from stabilizing selection about a line of optimal allometry between traits. Alternatively, introgression between populations and species could have caused correlation and variance to evolve jointly. Both selection and introgression were probably influential in producing the observed pattern, but it is not possible to estimate their relative importance with current data. Genetic and phenotypic correlations were correlated in the song sparrow, but heritabilities of traits varied greatly. As a result, the genetic variance-covariance matrix for traits is not simply a constant multiple of the phenotypic matrix. Evolutionary response to natural selection cannot, therefore, be predicted from the measurement of phenotypic characteristics alone.     

Genetic diversity and population structure in contemporary house sparrow populations along an urbanization gradient

House sparrow (Passer domesticus) populations have suffered major declines in urban as well as rural areas, while remaining relatively stable in suburban ones. Yet, to date no exhaustive attempt has been made to examine how, and to what extent, spatial variation in population demography is reflected in genetic population structuring along contemporary urbanization gradients. Here we use putatively neutral microsatellite loci to study if and how genetic variation can be partitioned in a hierarchical way among different urbanization classes. Principal coordinate analyses did not support the hypothesis that urban/suburban and rural populations comprise two distinct genetic clusters. Comparison of F(ST) values at different hierarchical scales revealed drift as an important force of population differentiation. Redundancy analyses revealed that genetic structure was strongly affected by both spatial variation and level of urbanization. The results shown here can be used as baseline information for future genetic monitoring programmes and provide additional insights into contemporary house sparrow dynamics along urbanization gradients.     

Genetic and Morphometric Divergence of an Invasive Bird: The Introduced House Sparrow (Passer domesticus) in Brazil

Introduced species are interesting systems for the study of contemporary evolution in new environments because of their spatial and temporal scales. For this study we had three aims: (i) to determine how genetic diversity and genetic differentiation of introduced populations of the house sparrow (Passer domesticus) in Brazil varies with range expansion, (ii) to determine how genetic diversity and differentiation in Brazil compares to ancestral European populations; and (iii) to determine whether selection or genetic drift has been more influential on phenotypic divergence. We used six microsatellite markers to genotype six populations from Brazil and four populations from Europe. We found slightly reduced levels of genetic diversity in Brazilian compared to native European populations. However, among introduced populations of Brazil, we found no association between genetic diversity and time since introduction. Moreover, overall genetic differentiation among introduced populations was low indicating that the expansion took place from large populations in which genetic drift effects would likely have been weak. We found significant phenotypic divergence among sites in Brazil. Given the absence of a spatial genetic pattern, divergent selection and not genetic drift seems to be the main force behind most of the phenotypic divergence encountered. Unravelling whether microevolution (e.g., allele frequency change), phenotypic plasticity, or both mediated phenotypic divergence is challenging and will require experimental work (e.g., common garden experiments or breeding programs).     

Patterns of morphological variation among populations of the widespread annual killifish Nothobranchius orthonotus are independent of genetic divergence and biogeography

Populations of annual killifish of the genus Nothobranchius occur in patchily distributed temporary pools in the East African savannah. Their fragmented distribution and low dispersal ability result in highly structured genetic clustering of their populations. In this study, we examined body shape variation in a widely distributed species, Nothobranchius orthonotus with known phylogeographic structure. We tested whether genetic divergence of major mitochondrial lineages forming two candidate species is congruent with phenotypic diversification, using linear and geometric morphometry analyses of body shape in 23 wild populations. We also conducted a common‐garden experiment with two wild‐derived populations to control for the effect of local environmental conditions on body shape. We identified different allometric trajectories for different mitochondrial lineages and candidate species in both sexes. However, in a principal components analysis of population‐level body shape, the separation among mitochondrial lineages was incomplete. Higher similarity of mitochondrial lineages belonging to different candidate species than that of same candidate species prevented distinction of the two candidate species on the basis of body shape. Analysis at the individual level demonstrated that N. orthonotus express high intrapopulation variability, with major overlap among individuals from all populations. In conclusion, we suggest that N. orthonotus be considered as a single species with an extensive geographic range, strong population genetic structure and high morphological variability.     

Surveillance for microbes and range expansion in house sparrows

Interactions between hosts and parasites influence the success of host introductions and range expansions post-introduction. However, the physiological mechanisms mediating these outcomes are little known. In some vertebrates, variation in the regulation of inflammation has been implicated, perhaps because inflammation imparts excessive costs, including high resource demands and collateral damage upon encounter with novel parasites. Here, we tested the hypothesis that variation in the regulation of inflammation contributed to the spread of house sparrows (Passer domesticus) across Kenya, one of the world''s most recent invasions of this species. Specifically, we asked whether inflammatory gene expression declines with population age (i.e. distance from Mombasa (dfM), the site of introduction around 1950). We compared expression of two microbe surveillance molecules (Toll-like receptors, TLRs-2 and 4) and a proinflammatory cytokine (interleukin-6, IL-6) before and after an injection of an immunogenic component of Gram-negative bacteria (lipopolysaccharide, LPS) among six sparrow populations. We then used a best-subset model selection approach to determine whether population age (dfM) or other factors (e.g. malaria or coccidian infection, sparrow density or genetic group membership) best-explained gene expression. For baseline expression of TLR-2 and TLR-4, population age tended to be the best predictor with expression decreasing with population age, although other factors were also important. Induced expression of TLRs was affected by LPS treatment alone. For induced IL-6, only LPS treatment reliably predicted expression; baseline expression was not explained by any factor. These data suggest that changes in microbe surveillance, more so than downstream control of inflammation via cytokines, might have been important to the house sparrow invasion of Kenya.     

When Anthropogenic River Disturbance Decreases Hybridisation between Non-Native and Endemic Cyprinids and Drives an Ecomorphological Displacement towards Juvenile State in Both Species

Understanding the impact of non-native species on native species is a major challenge in molecular ecology, particularly for genetically compatible fish species. Invasions are generally difficult to study because their effects may be confused with those of environmental or human disturbances. Colonized ecosystems are differently impacted by human activities, resulting in diverse responses and interactions between native and non-native species. We studied the dynamics between two Cyprinids species (invasive Chondrostoma nasus and endemic Parachondrostoma toxostoma) and their hybrids in 16 populations (from allopatric to sympatric situations and from little to highly fragmented areas) corresponding to 2,256 specimens. Each specimen was assigned to a particular species or to a hybrid pool using molecular identification (cytochrome b and 41 microsatellites). We carried out an ecomorphological analysis based on size, age, body shape, and diet (gut vacuity and molecular fecal contents). Our results contradicted our initial assumptions on the pattern of invasion and the rate of introgression. There was no sign of underperformance for the endemic species in areas where hybridisation occurred. In the unfragmented zone, the introduced species was found mostly downstream, with body shapes similar to those in allopatric populations while both species were found to be more insectivorous than the reference populations. However, high level of hybridisation was detected, suggesting interactions between the two species during spawning and/or the existence of hybrid swarm. In the disturbed zone, introgression was less frequent and slender body shape was associated with diatomivorous behaviour, smaller size (juvenile characteristics) and greater gut vacuity. Results suggested that habitat degradation induced similar ecomorphological trait changes in the two species and their hybrids (i.e. a transition towards a pedomorphic state) where the invasive species is more affected than the native species. Therefore, this study reveals a diversity of relationships between two genetically compatible species and emphasizes constraints on the invasion process in disturbed areas.     

Geographic structure of genetic and phenotypic variation in the hybrid zone between Quercus affinis and Q. laurina in Mexico

Analyzing the structure of hybrid zones is important for inferring their origin, dynamics and evolutionary significance. We examined the geographic structure of phenotypic and genetic variation in the contact zone between two Mexican red oaks, Quercus affinis and Q. laurina. A total of 105 individuals from seven populations were sampled along a 600‐km latitudinal gradient representing the distribution area of the two species and their contact zone. Individuals were genotyped for nine nuclear and four chloroplast DNA microsatellite loci (ncSSR and cpSSR, respectively), and characterized for several leaf and acorn traits. The cpSSR data revealed extensive haplotype sharing among populations of the two species, while a Bayesian assignment analysis based on ncSSRs identified two main genetic groups, each corresponding to one of the species, and two populations in the contact zone showing evidence of admixture. The proportion of genetic ancestry in the populations was strongly associated with latitude and showed a pattern of variation with the shape of a narrow sigmoidal cline. The variation in three of the seven phenotypic traits was partially congruent with molecular variation, while the other traits did not conform to a geographic cline but instead were correlated with environmental variables. In conclusion, the hybrid zone between the two oak species has some of the characteristics of a tension zone, but heterogeneous variation across traits suggests differential introgression and the action of extrinsic selection.