A novel tamanavirus (Flaviviridae) of the European common frog (Rana temporaria) from the UK.

Flavivirids are small, enveloped, positive-sense RNA viruses from the family Flaviviridae with genomes of ~9-13 kb. Metatranscriptomic analyses of metazoan organisms have revealed a diversity of flavivirus-like or flavivirid viral sequences in fish and marine invertebrate groups. However, no flavivirus-like virus has been identified in amphibians. To remedy this, we investigated the virome of the European common frog (Rana temporaria) in the UK, utilizing high-throughput sequencing at six catch locations. De novo assembly revealed a coding-complete virus contig of a novel flavivirid ~11.2 kb in length. The virus encodes a single ORF of 3456 aa and 5' and 3' untranslated regions (UTRs) of 227 and 666 nt, respectively. We named this virus Rana tamanavirus (RaTV), as BLASTp analysis of the polyprotein showed the closest relationships to Tamana bat virus (TABV) and Cyclopterus lumpus virus from Pteronotus parnellii and Cyclopterus lumpus, respectively. Phylogenetic analysis of the RaTV polyprotein compared to Flavivirus and Flavivirus-like members indicated that RaTV was sufficiently divergent and basal to the vertebrate Tamanavirus clade. In addition to the Mitcham strain, partial but divergent RaTV, sharing 95.64-97.39 % pairwise nucleotide identity, were also obtained from the Poole and Deal samples, indicating that RaTV is widespread in UK frog samples. Bioinformatic analyses of predicted secondary structures in the 3'UTR of RaTV showed the presence of an exoribonuclease-resistant RNA (xrRNA) structure standard in flaviviruses and TABV. To examine this biochemically, we conducted an in vitro Xrn1 digestion assay showing that RaTV probably forms a functional Xrn1-resistant xrRNA.

Gene Flow Limits Adaptation along Steep Environmental Gradients.

When environmental variation is spatially continuous, dispersing individuals move among nearby sites with similar habitat conditions. But as an environmental gradient becomes steeper, gene flow may connect more divergent habitats, and this is predicted to reduce the slope of the adaptive cline that evolves. We compared quantitative genetic divergence of Rana temporaria frog populations along a 2,000-m elevational gradient in eastern Switzerland (new experimental results) with divergence along a 1,550-km latitudinal gradient in Fennoscandia (previously published results). Both studies found significant countergradient variation in larval development rate (i.e., animals from cold climates developed more rapidly). The cline was weaker with elevation than with latitude. Animals collected on both gradients were genotyped at ∼2,000 single-nucleotide polymorphism markers, revealing that dispersal distance was 30% farther on the latitudinal gradient but 3.9 times greater with respect to environmental conditions on the elevational gradient. A meta-analysis of 19 experimental studies of anuran populations spanning temperature gradients revealed that countergradient variation in larval development, while significant overall, was weaker when measured on steeper gradients. These findings support the prediction that adaptive population divergence is less pronounced, and maladaptation more pervasive, on steep environmental gradients.

Are glacial refugia hotspots of speciation and cytonuclear discordances? Answers from the genomic phylogeography of Spanish common frogs.

Subdivided Pleistocene glacial refugia, best known as "refugia within refugia", provided opportunities for diverging populations to evolve into incipient species and/or to hybridize and merge following range shifts tracking the climatic fluctuations, potentially promoting extensive cytonuclear discordances and "ghost" mtDNA lineages. Here, we tested which of these opposing evolutionary outcomes prevails in northern Iberian areas hosting multiple historical refugia of common frogs (Rana cf. temporaria), based on a genomic phylogeography approach (mtDNA barcoding and RAD-sequencing). We found evidence for both incipient speciation events and massive cytonuclear discordances. On the one hand, populations from northwestern Spain (Galicia and Asturias, assigned to the regional endemic R. parvipalmata), are deeply-diverged at mitochondrial and nuclear genomes (~4 My of independent evolution), and barely admix with northeastern populations (assigned to R. temporaria sensu stricto) across a narrow hybrid zone (~25 km) located in the Cantabrian Mountains, suggesting that they represent distinct species. On the other hand, the most divergent mtDNA clade, widespread in Cantabria and the Basque country, shares its nuclear genome with other R. temporaria s. s. lineages. Patterns of population expansions and isolation-by-distance among these populations are consistent with past mitochondrial capture and/or drift in generating and maintaining this ghost mitochondrial lineage. This remarkable case study emphasizes the complex evolutionary history that shaped the present genetic diversity of refugial populations, and stresses the need to revisit their phylogeography by genomic approaches, in order to make informed taxonomic inferences.

Sex-Chromosome Recombination in Common Frogs Brings Water to the Fountain-of-Youth.

According to the canonical model of sex-chromosome evolution, the degeneration of Y or W chromosomes (as observed in mammals and birds, respectively) results from an arrest of recombination in the heterogametic sex, driven by the fixation of sexually antagonistic mutations. However, sex chromosomes have remained homomorphic in many lineages of fishes, amphibians, and nonavian reptiles. According to the "fountain-of-youth" model, this homomorphy results from occasional events of sex reversal. If recombination arrest in males is controlled by maleness per se (and not by genotype), then Y chromosomes are expected to recombine in XY females, preventing their long-term degeneration. Here, we provide field support for the fountain-of-youth, by showing that sex-chromosome recombination in Rana temporaria only depends on phenotypic sex: naturally occurring XX males show the same restriction of recombination as XY males (average map length ∼2 cM), while XY females recombine as much as XX females (average map length ∼150 cM). Our results challenge several common assumptions regarding the evolution of sex chromosomes, including the role of sexually antagonistic genes as drivers of recombination arrest, and that of chromosomal inversions as underlying mechanisms.

A reciprocal translocation radically reshapes sex-linked inheritance in the common frog.

X and Y chromosomes can diverge when rearrangements block recombination between them. Here we present the first genomic view of a reciprocal translocation that causes two physically unconnected pairs of chromosomes to be coinherited as sex chromosomes. In a population of the common frog (Rana temporaria), both pairs of X and Y chromosomes show extensive sequence differentiation, but not degeneration of the Y chromosomes. A new method based on gene trees shows both chromosomes are sex-linked. Furthermore, the gene trees from the two Y chromosomes have identical topologies, showing they have been coinherited since the reciprocal translocation occurred. Reciprocal translocations can thus reshape sex linkage on a much greater scale compared with inversions, the type of rearrangement that is much better known in sex chromosome evolution, and they can greatly amplify the power of sexually antagonistic selection to drive genomic rearrangement. Two more populations show evidence of other rearrangements, suggesting that this species has unprecedented structural polymorphism in its sex chromosomes.

Causes and consequences of reciprocal translocations on sex chromosomes.

Under XY sex determination, the Y chromosome is only inherited via males, whereas the X chromosome is predominantly found in females. Thus, it is favourable when alleles with high male fitness become associated with the Y chromosome and when alleles with high female fitness become associated with the X chromosome. These favourable associations can be strengthened through linkage. Rearrangements, such as inversions and sex chromosome-autosome fusions, can increase linkage and thereby become favoured (Charlesworth, 2017). In a From the Cover article in this issue of Molecular Ecology, Toups, Rodrigues, Perrin, and Kirkpatrick (2019) present the first genomic analysis of a sex chromosome reciprocal translocation, a particularly dramatic chromosomal rearrangement that modifies linkage with the sex chromosome. As a result of reciprocal translocation, one studied population of the common frog (Rana temporaria, Figure 1) displays a remarkable sex-determining system in which there are two physically unlinked sex chromosomes that are exclusively cotransmitted (Figure 2a).

Mapping of quantitative trait loci for life history traits segregating within common frog populations.

The evolution of complex traits is often shaped by adaptive divergence. However, very little is known about the number, effect size, and location of the genomic regions influencing the variation of these traits in natural populations. Based on a dense linkage map of the common frog, Rana temporaria, we have localized, for the first time in amphibians, three significant and nine suggestive quantitative trait loci (QTLs) for metabolic rate, growth rate, development time, and weight at metamorphosis, explaining 5.6-18.9% of the overall phenotypic variation in each trait. We also found a potential pleiotropic QTL between development time and size at metamorphosis that, if confirmed, might underlie the previously reported genetic correlation between these traits. Furthermore, we demonstrate that the genetic variation linked to fitness-related larval traits segregates within Rana temporaria populations. This study provides the first insight into the genomic regions that affect larval life history traits in anurans, providing a valuable resource to delve further into the genomic basis of evolutionary change in amphibians.

A novel approach to wildlife transcriptomics provides evidence of disease-mediated differential expression and changes to the microbiome of amphibian populations.

Ranaviruses are responsible for a lethal, emerging infectious disease in amphibians and threaten their populations throughout the world. Despite this, little is known about how amphibian populations respond to ranaviral infection. In the United Kingdom, ranaviruses impact the common frog (Rana temporaria). Extensive public engagement in the study of ranaviruses in the UK has led to the formation of a unique system of field sites containing frog populations of known ranaviral disease history. Within this unique natural field system, we used RNA sequencing (RNA-Seq) to compare the gene expression profiles of R. temporaria populations with a history of ranaviral disease and those without. We have applied a RNA read-filtering protocol that incorporates Bloom filters, previously used in clinical settings, to limit the potential for contamination that comes with the use of RNA-Seq in nonlaboratory systems. We have identified a suite of 407 transcripts that are differentially expressed between populations of different ranaviral disease history. This suite contains genes with functions related to immunity, development, protein transport and olfactory reception among others. A large proportion of potential noncoding RNA transcripts present in our differentially expressed set provide first evidence of a possible role for long noncoding RNA (lncRNA) in amphibian response to viruses. Our read-filtering approach also removed significantly more bacterial reads from libraries generated from positive disease history populations. Subsequent analysis revealed these bacterial read sets to represent distinct communities of bacterial species, which is suggestive of an interaction between ranavirus and the host microbiome in the wild.

Dmrt1 polymorphism and sex-chromosome differentiation in Rana temporaria.

Sex-determination mechanisms vary both within and among populations of common frogs, opening opportunities to investigate the molecular pathways and ultimate causes shaping their evolution. We investigated the association between sex-chromosome differentiation (as assayed from microsatellites) and polymorphism at the candidate sex-determining gene Dmrt1 in two Alpine populations. Both populations harboured a diversity of X-linked and Y-linked Dmrt1 haplotypes. Some males had fixed male-specific alleles at all markers ("differentiated" Y chromosomes), others only at Dmrt1 ("proto-" Y chromosomes), while still others were genetically indistinguishable from females (undifferentiated X chromosomes). Besides these XX males, we also found rare XY females. The several Dmrt1 Y haplotypes differed in the probability of association with a differentiated Y chromosome, which we interpret as a result of differences in the masculinizing effects of alleles at the sex-determining locus. From our results, the polymorphism in sex-chromosome differentiation and its association with Dmrt1, previously inferred from Swedish populations, are not just idiosyncratic features of peripheral populations, but also characterize highly diverged populations in the central range. This implies that an apparently unstable pattern has been maintained over long evolutionary times.

Molecular phenotyping of maternally mediated parallel adaptive divergence within Rana arvalis and Rana temporaria.

When similar selection acts on the same traits in multiple species or populations, parallel evolution can result in similar phenotypic changes, yet the underlying molecular architecture of parallel phenotypic divergence can be variable. Maternal effects can influence evolution at ecological timescales and facilitate local adaptation, but their contribution to parallel adaptive divergence is unclear. In this study, we (i) tested for variation in embryonic acid tolerance in a common garden experiment and (ii) used molecular phenotyping of egg coats to investigate the molecular basis of maternally mediated parallel adaptive divergence in two amphibian species (Rana arvalis and Rana temporaria). Our results on three R. arvalis and two R. temporaria populations show that adaptive divergence in embryonic acid tolerance is mediated via maternally derived egg coats in both species. We find extensive polymorphism in egg jelly coat glycoproteins within both species and that acid-tolerant clutches have more negatively charged egg jelly - indicating that the glycosylation status of the jelly coat proteins is under divergent selection in acidified environments, likely due to its impact on jelly water balance. Overall, these data provide evidence for parallel mechanisms of adaptive divergence in two species. Our study highlights the importance of studying intraspecific molecular variation in egg coats and, specifically, their glycoproteins, to increase understanding of underlying forces maintaining variation in jelly coats.

High-density linkage maps fail to detect any genetic component to sex determination in a Rana temporaria family.

Sex chromosome differentiation in Rana temporaria varies strikingly among populations or families: whereas some males display well-differentiated Y haplotypes at microsatellite markers on linkage group 2 (LG2), others are genetically undistinguishable from females. We analysed with RADseq markers one family from a Swiss lowland population with no differentiated sex chromosomes, and where sibship analyses had failed to detect any association between the phenotypic sex of progeny and parental haplotypes. Offspring were reared in a common tank in outdoor conditions and sexed at the froglet stage. We could map a total of 2177 SNPs (1123 in the mother, 1054 in the father), recovering in both adults 13 linkage groups (= chromosome pairs) that were strongly syntenic to Xenopus tropicalis despite > 200 My divergence. Sexes differed strikingly in the localization of crossovers, which were uniformly distributed in the female but limited to chromosome ends in the male. None of the 2177 markers showed significant association with offspring sex. Considering the very high power of our analysis, we conclude that sex determination was not genetic in this family; which factors determined sex remain to be investigated.

The genetic contribution to sex determination and number of sex chromosomes vary among populations of common frogs (Rana temporaria).

The patterns of sex determination and sex differentiation have been shown to differ among geographic populations of common frogs. Notably, the association between phenotypic sex and linkage group 2 (LG2) has been found to be perfect in a northern Swedish population, but weak and variable among families in a southern one. By analyzing these populations with markers from other linkage groups, we bring two new insights: (1) the variance in phenotypic sex not accounted for by LG2 in the southern population could not be assigned to genetic factors on other linkage groups, suggesting an epigenetic component to sex determination; (2) a second linkage group (LG7) was found to co-segregate with sex and LG2 in the northern population. Given the very short timeframe since post-glacial colonization (in the order of 1000 generations) and its seemingly localized distribution, this neo-sex chromosome system might be the youngest one described so far. It does not result from a fusion, but more likely from a reciprocal translocation between the original Y chromosome (LG2) and an autosome (LG7), causing their co-segregation during male meiosis. By generating a strict linkage between several important genes from the sex-determination cascade (Dmrt1, Amh and Amhr2), this neo-sex chromosome possibly contributes to the 'differentiated sex race' syndrome (strictly genetic sex determination and early gonadal development) that characterizes this northern population.

Sex-chromosome differentiation and 'sex races' in the common frog (Rana temporaria).

Sex-chromosome differentiation was recently shown to vary among common frog populations in Fennoscandia, suggesting a trend of increased differentiation with latitude. By rearing families from two contrasted populations (respectively, from northern and southern Sweden), we show this disparity to stem from differences in sex-determination mechanisms rather than in XY-recombination patterns. Offspring from the northern population display equal sex ratios at metamorphosis, with phenotypic sexes that correlate strongly with paternal LG2 haplotypes (the sex chromosome); accordingly, Y haplotypes are markedly differentiated, with male-specific alleles and depressed diversity testifying to their smaller effective population size. In the southern population, by contrast, a majority of juveniles present ovaries at metamorphosis; only later in development do sex ratios return to equilibrium. Even at these later stages, phenotypic sexes correlate only mildly with paternal LG2 haplotypes; accordingly, there are no recognizable Y haplotypes. These distinct patterns of gonadal development fit the concept of 'sex races' proposed in the 1930s, with our two populations assigned to the 'differentiated' and 'semi-differentiated' races, respectively. Our results support the suggestion that 'sex races' differ in the genetic versus epigenetic components of sex determination. Analysing populations from the 'undifferentiated race' with high-density genetic maps should help to further test this hypothesis.

Local adaptation with high gene flow: temperature parameters drive adaptation to altitude in the common frog (Rana temporaria).

Both environmental and genetic influences can result in phenotypic variation. Quantifying the relative contributions of local adaptation and phenotypic plasticity to phenotypes is key to understanding the effect of environmental variation on populations. Identifying the selective pressures that drive divergence is an important, but often lacking, next step. High gene flow between high- and low-altitude common frog (Rana temporaria) breeding sites has previously been demonstrated in Scotland. The aim of this study was to assess whether local adaptation occurs in the face of high gene flow and to identify potential environmental selection pressures that drive adaptation. Phenotypic variation in larval traits was quantified in R. temporaria from paired high- and low-altitude sites using three common temperature treatments. Local adaptation was assessed using Q(ST)-F(ST) analyses, and quantitative phenotypic divergence was related to environmental parameters using Mantel tests. Although evidence of local adaptation was found for all traits measured, only variation in larval period and growth rate was consistent with adaptation to altitude. Moreover, this was only evident in the three mountains with the highest high-altitude sites. This variation was correlated with mean summer and winter temperatures, suggesting that temperature parameters are potentially strong selective pressures maintaining local adaptation, despite high gene flow.

Geographic variation in sex-chromosome differentiation in the common frog (Rana temporaria).

In sharp contrast with birds and mammals, sex-determination systems in ectothermic vertebrates are often highly dynamic and sometimes multifactorial. Both environmental and genetic effects have been documented in common frogs (Rana temporaria). One genetic linkage group, mapping to the largest pair of chromosomes and harbouring the candidate sex-determining gene Dmrt1, associates with sex in several populations throughout Europe, but association varies both within and among populations. Here, we show that sex association at this linkage group differs among populations along a 1500-km transect across Sweden. Genetic differentiation between sexes is strongest (FST  = 0.152) in a northern-boreal population, where male-specific alleles and heterozygote excesses (FIS  = -0.418 in males, +0.025 in females) testify to a male-heterogametic system and lack of X-Y recombination. In the southernmost population (nemoral climate), in contrast, sexes share the same alleles at the same frequencies (FST  = 0.007 between sexes), suggesting unrestricted recombination. Other populations show intermediate levels of sex differentiation, with males falling in two categories: some cluster with females, while others display male-specific Y haplotypes. This polymorphism may result from differences between populations in the patterns of X-Y recombination, co-option of an alternative sex-chromosome pair, or a mixed sex-determination system where maleness is controlled either by genes or by environment depending on populations or families. We propose approaches to test among these alternative models, to disentangle the effects of climate and phylogeography on the latitudinal trend, and to sort out how this polymorphism relates to the 'sexual races' described in common frogs in the 1930s.

Using genetic variation to infer associations with climate in the common frog, Rana temporaria.

Recent and historical species' associations with climate can be inferred using molecular markers. This knowledge of population and species-level responses to climatic variables can then be used to predict the potential consequences of ongoing climate change. The aim of this study was to predict responses of Rana temporaria to environmental change in Scotland by inferring historical and contemporary patterns of gene flow in relation to current variation in local thermal conditions. We first inferred colonization patterns within Europe following the last glacial maximum by combining new and previously published mitochondrial DNA sequences. We found that sequences from our Scottish samples were identical to (92%), or clustered with, the common haplotype previously identified from Western Europe. This clade showed very low mitochondrial variation, which did not allow inference of historical colonization routes but did allow interpretation of patterns of current fine-scale population structure without consideration of confounding historical variation. Second, we assessed fine-scale microsatellite-based patterns of genetic variation in relation to current altitudinal temperature gradients. No population structure was found within altitudinal gradients (average FST=0.02), despite a mean annual temperature difference of 4.5 °C between low- and high-altitude sites. Levels of genetic diversity were considerable and did not vary between sites. The panmictic population structure observed, even along temperature gradients, is a potentially positive sign for R. temporaria persistence in Scotland in the face of a changing climate. This study demonstrates that within taxonomic groups, thought to be at high risk from environmental change, levels of vulnerability can vary, even within species.

Phenotypic plasticity in the hepatic transcriptome of the European common frog (Rana temporaria): the interplay between environmental induction and geographical lineage on developmental response.

Phenotypic plasticity might facilitate adaptation to new environmental conditions through the enhancement of initial survival of organisms. Once a population is established, further adaptation and diversification may occur through adaptive trait evolution. While several studies have found evidence for this mechanism using phenotypic traits, much less is known at the level of gene expression. Here, we use an islands system of frog populations that show local adaptation and phenotypic plasticity to pool drying conditions in development time until metamorphoses. We examined gene expression differences in Rana temporaria tadpole livers with respect to pool drying at the source population and in response to simulated pool drying in the laboratory. Using a MAGEX cDNA microarray and quantitative real-time polymerase chain reaction (qPCR), we identified an increase in several gene transcripts in response to artificial pool drying including thyroid hormone receptor alpha and beta, carbamoyl phosphate synthetase 1, ornithine transcarbamylase and catalase. In addition, these gene transcripts also showed greater abundance in island populations that developed faster. Hence, the gene transcripts were related to both constitutive response (higher levels in island populations that developed faster) and plastic response (increased abundance under decreasing water levels). This pattern is in accordance with genetic accommodation, which predicts similarities between plastic gene expression and constitutive expression in locally adapted populations.

Within-population polymorphism of sex-determination systems in the common frog (Rana temporaria).

In sharp contrast with birds and mammals, the sex chromosomes of ectothermic vertebrates are often undifferentiated, for reasons that remain debated. A linkage map was recently published for Rana temporaria (Linnaeus, 1758) from Fennoscandia (Eastern European lineage), with a proposed sex-determining role for linkage group 2 (LG2). We analysed linkage patterns in lowland and highland populations from Switzerland (Western European lineage), with special focus on LG2. Sibship analyses showed large differences from the Fennoscandian map in terms of recombination rates and loci order, pointing to large-scale inversions or translocations. All linkage groups displayed extreme heterochiasmy (total map length was 12.2 cM in males, versus 869.8 cM in females). Sex determination was polymorphic within populations: a majority of families (with equal sex ratios) showed a strong correlation between offspring phenotypic sex and LG2 paternal haplotypes, whereas other families (some of which with female-biased sex ratios) did not show any correlation. The factors determining sex in the latter could not be identified. This coexistence of several sex-determination systems should induce frequent recombination of X and Y haplotypes, even in the absence of male recombination. Accordingly, we found no sex differences in allelic frequencies on LG2 markers among wild-caught male and female adults, except in one high-altitude population, where nonrecombinant Y haplotypes suggest sex to be entirely determined by LG2. Multifactorial sex determination certainly contributes to the lack of sex-chromosome differentiation in amphibians.

Refugia within refugia as a key to disentangle the genetic pattern of a highly variable species: the case of Rana temporaria Linnaeus, 1758 (Anura, Ranidae).

Two distinct lineages of Rana temporaria are known in the Palaearctic region, but it is uncertain whether this species persisted in one or more Pleistocene refugia. We resolved the phylogeographic history and genetic variability of R. temporaria in the Italian peninsula, a 'traditional' Pleistocene refugium, and related our findings to patterns described for other European populations. We sequenced the mitochondrial markers Cox I and cytochrome b. Phylogenetic reconstruction only indicated the presence of haplotypes belonging to the Western lineage in the Italian peninsula. Overall, the genetic variability of Italian populations was higher than other European populations, which shared haplotypes with the Alpine populations. We demonstrated subdivision into five main Italian sublineages, which was associated with a geographical structure of populations in two divergent groups. In particular, one Apennine group might have resulted from bottlenecks during the last interglacials ages. In contrast, Alpine populations were recently diverged and showed incomplete lineage sorting. Our data indicate that the Italian peninsula served as refugium for the Western lineage of R. temporaria. Dispersion towards Central Europe probably started only from the western slope of the Alps via a rapid leading edge expansion. The identified structure is partially congruent with traditional peripheral refugia identified for plants. This evolutionary scenario does not support any taxonomic distinction at the subspecific level for R. temporaria.

Discordant patterns of nuclear and mitochondrial introgression in Iberian populations of the European common frog (Rana temporaria).

Amphibians often show complex histories of intraspecific and interspecific genetic introgression, which might differ in mitochondrial and nuclear genes. In our study of the genetic differentiation of the European common frog, Rana temporaria (159 specimens from 23 populations were analyzed for 24 presumptive allozyme loci; 82 specimens were sequenced for a 540-bp fragment of the mitochondrial 16S rRNA gene), multilocus correspondence analysis (CA) and Bayesian assignment tests of the nuclear data were concordant in identifying 2 population groups corresponding to 1) the Pyrenees in the east and 2) the Galicia and Asturias regions in the west, the latter corresponding to the subspecies R. temporaria parvipalmata. Geographically intermediate populations were genetically intermediate in the allozyme CA and, less clearly in the Bayesian assignment, with mitochondrial haplotypes exclusively belonging to the parvipalmata group. This indicates different degrees of introgression in the mitochondrial and nuclear genomes. Although Pyrenean high-altitude populations are morphologically distinct from low-altitude populations, these 2 groups were not separate clusters in any analysis. This suggests that the morphological differences may be due to fast adaptation to elevational gradients, likely under maintenance of gene flow, and that the underlying genetic changes are not detectable by the analyzed markers. We argue that a parsimonious explanation for the observed pattern along the east-west axis in northern Spain may be competition between invading and resident populations, with no need to invoke selection. However, in order to conclusively rule out selective processes, additional and finer scale data are required to test for asymmetric mating preference/behaviour, sex-biased gene flow, or sex-biased survival of potential hybrids.