Hybridization and mitochondrial genome introgression between Rana chensinensis and R. kukunoris

Mitochondrial genome (mito‐genome) introgression among metazoans is commonplace, and several biological processes may promote such introgression. We examined two proposed processes for the mito‐genome introgression between Rana chensinensis and R. kukunoris: natural hybridization and sex‐biased dispersal. We sampled 477 individuals from 28 sites in the potential hybrid zone in the western Tsinling Mountains. Mitochondrial gene (cyt‐b) trees were used to examine the introgression events. Microsatellite DNA loci, cyt‐b and morphological data were used to identify hybrids and to examine the extent of natural hybridization. We detected rampant bidirectional introgressions, both ancient and recent, between the two species. Furthermore, we found a wide hybrid zone, and frequent and asymmetric hybridization. The hybrid zone cline analysis revealed a clear mitochondrial–nuclear discordance; while most nuclear markers displayed similar and steep clines, cyt‐b had a displaced cline centre and a more gradual and wider cline. We also detected strong and asymmetric historical maternal gene flow across the hybrid zone. This widespread hybridization and detected low mito‐nuclear conflicts may, at least partially, explain the high frequency of introgression. Lastly, microsatellite data and population genetic methods were used to assess sex‐biased dispersal. A weak pattern of female‐biased dispersal was detected in both species, suggesting it may not play an important role in the observed introgression. Our data are consistent with the hybridization hypothesis, but support for the sex‐biased dispersal hypothesis is weak. We further suggest that selective advantages of the R. kukunoris‐type mito‐genome in thermal adaptation may also contribute to the introgression between the two species.     

Limited gene exchange between two sister species of leaf beetles within a hybrid zone in the Alps

Analysing genomic variation within and between sister species is a first step towards understanding species boundaries. We focused on two sister species of cold‐resistant leaf beetles, Gonioctena quinquepunctata and G. intermedia, whose ranges overlap in the Alps. A previous study of DNA sequence variation had revealed multiple instances of mitochondrial genome introgression in this region, suggesting recent hybridization between the two species. To evaluate the extent of gene exchange resulting from these hybridization events, we sampled individuals of both species inside and outside the hybrid zone and analysed genomic variation among them using RAD‐seq markers. Individual levels of introgression in the nuclear genome were estimated first by defining species‐specific SNPs (displaying a fixed difference between species) a priori and second by using model‐based methods. Both types of analyses indicated little gene exchange, if any, between species at the level of the nuclear genome. Whereas the first method suggested slightly more gene flow, we argue that it has likely overestimated introgression in the phylogeographic context of this study. We conclude that strong intrinsic barriers prevent genetic exchange at the level of the nuclear genome between the two species. The apparent discrepancy observed between introgression occurring in the nuclear and mitochondrial genomes could be explained by selection acting in favour of the latter. Also, these results have consequences for the phylogeographic study of each species, since we can assume that genetic diversity in the overlapping portion of their ranges is not the product of introgression.     

Spatial and temporal patterns of environmental DNA detection to inform sampling protocols in lentic and lotic systems

The development of efficient sampling protocols for the capture of environmental DNA (eDNA) could greatly help improve accuracy of occupancy monitoring for species that are difficult to detect. However, the process of developing a protocol in situ is complicated for rare species by the fact that animal locations are often unknown. We tested sampling designs in lake and stream systems to determine the most effective eDNA sampling protocols for two rare species: the Sierra Nevada yellow‐legged frog (Rana sierrae) and the foothill yellow‐legged frog (Rana boylii). We varied water volume, spatial sampling, and seasonal timing in lakes and streams; in lakes we also tested multiple filter types. We found that filtering 2 L versus 1 L increased the odds of detection in streams 5.42X (95% CI: 3.2–9.19X) in our protocol, from a probability of 0.51–0.85 per technical replicate. Lake sample volumes were limited by filter clogging, and we found no effect of volume or filter type. Sampling later in the season increased the odds of detection in streams by 1.96X for every 30 days (95% CI: 1.3–2.97X) but there was no effect for lakes. Spatial autocorrelation of the quantity of yellow‐legged frog eDNA captured in streams ceased between 100 and 200 m, indicating that sampling at close intervals is important.     

Increased genetic divergence between two closely related fir species in areas of range overlap

Because of introgressive hybridization, closely related species can be more similar to each other in areas of range overlap (parapatry or sympatry) than in areas where they are geographically isolated from each other (allopatry). Here, we report the reverse situation based on nuclear genetic divergence between two fir species, Abies chensiensis and Abies fargesii, in China, at sites where they are parapatric relative to where they are allopatric. We examined genetic divergence across 126 amplified fragment length polymorphism (AFLP) markers in a set of 172 individuals sampled from both allopatric and parapatric populations of the two species. Our analyses demonstrated that AFLP divergence was much greater between the species when comparisons were made between parapatric populations than between allopatric populations. We suggest that selection in parapatry may have largely contributed to this increased divergence.     

Contrasting genetic patterns between two coexisting Eleutherococcus species in northern China

Climate oscillations are the key factors to understand the patterns in modern biodiversity. East Asia harbors the most diverse temperate flora, largely because an extensive terrestrial ice cap was absent during repeated Pleistocene glaciation–interglacial cycles. Comparing the demographic histories of species that are codistributed and are close relatives may provide insight into how the process of climate change influences species ranges. In this study, we compared the spatial genetic structure and demographic histories of two coexisting Eleutherococcus species, Eleutherococcus senticosus and E. sessiliflorus. Both species are distributed in northern China, regions that are generally considered to be sensitive to climatic fluctuations. These regions once hosted temperate forest, but this temperate forest was replaced by tundra and taiga forest during the Last Glacial Maximum (LGM), according to pollen records. Using three chloroplast DNA fragments, we assessed the genetic structure of 20 and 9 natural populations of E. senticosus and E. sessiliflorus, respectively. Extremely contrasting genetic patterns were found between the two species; E. sessiliflorus had little genetic variation, whereas E. senticosus had considerably higher levels of genetic variation (15 haplotypes). We speculated that a recent severe bottleneck may have resulted in the extremely low genetic diversity in E. sessiliflorus. In E. senticosus, populations in Northeast China (NEC) harbored all of the haplotypes found in this species and included private haplotypes. The populations in NEC had higher levels of genetic diversity than did those from North China (NC). Therefore, we suggest that both the NC and NEC regions can sustain LGM refugia and that lineage admixture from multiple refugia took place after the LGM elevated the local genetic diversity in NEC. In NEC, multiple genetic hot spots were found in the Changbai Mountains and the Xiaoxing'an Range, which implied that multiple locations in NEC may sustain LGM refugia, even in the Xiaoxing'an Range.     

Hybridization and introgression among three Aconitum (Ranunculaceae) species of different ploidy levels in the Tatra Mountains (Western Carpathians)

Hybridization among species of A conitum effects their morphology and ecology. In this study the hybridization between the diploid 2n(2x) = 16 ( A . lasiocarpum and A . variegatum) and tetraploid 2n(4x) = 32 ( A . firmum) species was documented in the Tatra Mountains (Western Carpathians) in a small, local population in a semi‐natural site. The hybrid species were: homoploid A . × pawlowskii ( A . lasiocarpum × A . variegatum), and triploid A . × berdaui ( A . firmum × A . variegatum, 2n(3x) = 24). Chloroplast DNA (cpDNA) alleles formed two distinct haplotypes, one typical for the tetraploid and another for diploid lines, shared between the tetraploid, triploid and diploid groups, indicating introgressive hybridization. The presumed gene flow was from the tetraploid to diploid species via the triploid bridge. The only two specimens of A . × pawlowskii that harbored tetraploid ( A . firmum) type cpDNA possessed bracteoles of A . firmum‐type. The remaining introgressed (cpDNA and Inter Simple Sequence Repeats (ISSR)) specimens ( A . variegatum) were morphologically pure, implying cryptic introgression. ISSR loci shared between the tetraploid A . firmum and diploid A . variegatum support the hypothesis of an ancient allopolyploid origin of A . firmum and the diploid species of A . variegatum‐type as one of its parent.     

Nuclear introgression without mitochondrial introgression in two turtle species exhibiting sex‐specific trophic differentiation

Despite the presence of reproductive barriers between species, interspecific gene introgression has been documented in a range of natural systems. Comparing patterns of genetic introgression in biparental versus matrilineal markers can potentially reveal sex‐specific barriers to interspecific gene flow. Hybridization has been documented in the freshwater turtles Graptemys geographica and G. pseudogeographica, whose ranges are largely sympatric. Morphological differentiation between the species is restricted to females, with female G. geographica possessing large heads and jaws compared to the narrow heads of G. pseudogeographica females. If hybrid females are morphologically intermediate, they may be less successful at exploiting parental feeding niches, thereby limiting the introgression of maternally inherited, but not biparental, molecular markers. We paired sequence data with stable isotope analysis and examined sex‐specific genetic introgression and trophic differentiation in sympatric populations of G. geographica and G. pseudogeographica. We observed introgression from G. pseudogeographica into G. geographica at three nuclear loci, but not at the mitochondrial locus. Analysis of ?¹⁵N and ?¹³C was consistent with species differences in trophic positioning in females, but not males. These results suggest that ecological divergence in females may reduce the opportunity for gene flow in this system.     

Hybridization between Felis silvestris silvestris and Felis silvestris catus in two contrasted environments in France

European wildcat (Felis silvestris silvestris) populations are fragmented throughout most of the whole range of the subspecies and may be threatened by hybridization with the domestic cat F.s. catus. The underlying ecological processes promoting hybridization remain largely unknown. In France, wildcats are mainly present in the northeast and signs of their presence in the Pyrenees have been recently provided. However, no studies have been carried out in the French Pyrenees to assess their exposure to hybridization. We compared two local populations of wildcats, one living in a continuous forest habitat in the French Pyrenees, the other living in a highly fragmented forest‐agricultural landscape in northeastern France to get insights into the variability of hybridization rates. Strong evidence of hybridization was detected in northeastern France and not in the Pyrenees. Close kin in the Pyrenees were not found in the same geographic location contrary to what was previously reported for females in the northeastern wildcat population. The two wildcat populations were significantly differentiated (F_ST = 0.072) to an extent close to what has been reported (F_ST = 0.103) between the Iberian population, from which the Pyrenean population may originate, and the German population, which is connected to the northeastern population. The genetic diversity of the Pyrenean wildcats was lower than that of northeastern wildcat populations in France and in other parts of Europe. The lower hybridization in the Pyrenees may result from the continuity of natural forest habitats. Further investigations should focus on linking landscape features to hybridization rates working on local populations.     

Persistence at distributional edges: Columbia spotted frog habitat in the arid Great Basin,USA

A common challenge in the conservation of broadly distributed, yet imperiled species is understanding which factors facilitate persistence at distributional edges, locations where populations are often vulnerable to extirpation due to changes in climate, land use, or distributions of other species. For Columbia spotted frogs (Rana luteiventris) in the Great Basin (USA), a genetically distinct population segment of conservation concern, we approached this problem by examining (1) landscape‐scale habitat availability and distribution, (2) water body‐scale habitat associations, and (3) resource management‐identified threats to persistence. We found that areas with perennial aquatic habitat and suitable climate are extremely limited in the southern portion of the species’ range. Within these suitable areas, native and non‐native predators (trout and American bullfrogs [Lithobates catesbeianus]) are widespread and may further limit habitat availability in upper‐ and lower‐elevation areas, respectively. At the water body scale, spotted frog occupancy was associated with deeper sites containing abundant emergent vegetation and nontrout fish species. Streams with American beaver (Castor canadensis) frequently had these structural characteristics and were significantly more likely to be occupied than ponds, lakes, streams without beaver, or streams with inactive beaver ponds, highlighting the importance of active manipulation of stream environments by beaver. Native and non‐native trout reduced the likelihood of spotted frog occupancy, especially where emergent vegetation cover was sparse. Intensive livestock grazing, low aquatic connectivity, and ephemeral hydroperiods were also negatively associated with spotted frog occupancy. We conclude that persistence of this species at the arid end of its range has been largely facilitated by habitat stability (i.e., permanent hydroperiod), connectivity, predator‐free refugia, and a commensalistic interaction with an ecosystem engineer. Beaver‐induced changes to habitat quality, stability, and connectivity may increase spotted frog population resistance and resilience to seasonal drought, grazing, non‐native predators, and climate change, factors which threaten local or regional persistence.     

Maintaining their genetic distance: Little evidence for introgression between widely hybridizing species of Geum with contrasting mating systems

Within the plant kingdom, many genera contain sister lineages with contrasting outcrossing and inbreeding mating systems that are known to hybridize. The evolutionary fate of these sister lineages is likely to be influenced by the extent to which they exchange genes. We measured gene flow between outcrossing Geum rivale and selfing Geum urbanum, sister species that hybridize in contemporary populations. We generated and used a draft genome of G. urbanum to develop dd‐RAD data scorable in both species. Coalescent analysis of RAD data from allopatric populations indicated that the species diverged 2–3 Mya, and that historical gene flow between them was extremely low (1 migrant every 25 generations). Comparison of genetic divergence between species in sympatry and allopatry, together with an analysis of allele frequencies in potential parental and hybrid populations, provided no evidence of contemporary introgression in sympatric populations. Cluster‐ and species‐specific marker analyses revealed that, apart from four early‐generation hybrids, individuals in sympatric populations fell into two genetically distinct groups that corresponded exactly to their morphological species classification with maximum individual admixture estimates of only 1–3%. However, we did observe joint segregation of four putatively introgressed SNPs across two scaffolds in the G. urbanum population that was associated with significant morphological variation, interpreted as tentative evidence for rare, recent interspecific gene flow. Overall, our results indicate that despite the presence of hybrids in contemporary populations, genetic exchange between G. rivale and G. urbanum has been extremely limited throughout their evolutionary history.     

From hybridization to introgression between two closely related sympatric ant species

Interspecific hybridization is becoming more frequent worldwide due to increasing global changes and translocations of organisms. For individual organisms, the most significant negative consequences are sterility or inviability of hybrid offspring. However, hybridization sometimes leads to fertile offspring, promoting introgression from one species into another. In such situations, hybridization can play a key role in evolution and speciation. Combining hypervariable DNA (microsatellites) and mitochondrial DNA markers with the use of several modeling methods allow an efficient detection of hybridization processes. The present study therefore investigates hybridization between two ant species, Tetramorium immigrans and T. caespitum, using multiple methods, and systematically comparing results with simulated data to ensure accurate identification of hybrids. Introgression was revealed both by backcross detection based on 14 nuclear microsatellite loci and by mitochondrial‐nuclear discordance based on comparison with mitochondrial DNA (cytochrome c oxidase subunit I). Results were spatially consistent, with hybrids located at latitudes where parental species are sympatric. The causes and consequences of hybridization and introgression between T. caespitum and T. immigrans remain to be further investigated, especially because T. immigrans could be an invasive species in France.     

Fungal infection alters the selection,dispersal and drift processes structuring the amphibian skin microbiome

Symbiotic microbial communities are important for host health, but the processes shaping these communities are poorly understood. Understanding how community assembly processes jointly affect microbial community composition is limited because inflexible community models rely on rejecting dispersal and drift before considering selection. We developed a flexible community assembly model based on neutral theory to ask: How do dispersal, drift and selection concurrently affect the microbiome across environmental gradients? We applied this approach to examine how a fungal pathogen affected the assembly processes structuring the amphibian skin microbiome. We found that the rejection of neutrality for the amphibian microbiome across a fungal gradient was not strictly due to selection processes, but was also a result of species‐specific changes in dispersal and drift. Our modelling framework brings the qualitative recognition that niche and neutral processes jointly structure microbiomes into quantitative focus, allowing for improved predictions of microbial community turnover across environmental gradients.     

Hybridization promotes color polymorphism in the aposematic harlequin poison frog,Oophaga histrionica

Whether hybridization can be a mechanism that drives phenotypic diversity is a widely debated topic in evolutionary biology. In poison frogs (Dendrobatidae), assortative mating has been invoked to explain how new color morphs persist despite the expected homogenizing effects of natural selection. Here, we tested the complementary hypothesis that new morphs arise through hybridization between different color morphs. Specifically, we (1) reconstructed the phylogenetic relationships among the studied populations of a dart‐poison frog to provide an evolutionary framework, (2) tested whether microsatellite allele frequencies of one putative hybrid population of the polymorphic frog O. histrionica are intermediate between O. histrionica and O. lehmanni, and (3) conducted mate‐choice experiments to test whether putatively intermediate females prefer homotypic males over males from the other two populations. Our findings are compatible with a hybrid origin for the new morph and emphasize the possibility of hybridization as a mechanism generating variation in polymorphic species. Moreover, because coloration in poison frogs is aposematic and should be heavily constrained, our findings suggest that hybridization can produce phenotypic novelty even in systems where phenotypes are subject to strong stabilizing selection.     

Climatic influences on fire regimes in the northern Sierra Nevada mountains, Lake Tahoe Basin, Nevada, USA

Aim The goal of this study was to understand better the role of interannual and interdecadal climatic variation on local pre‐EuroAmerican settlement fire regimes in fire‐prone Jeffrey pine (Pinus jeffreyi Grev. & Balf.) dominated forests in the northern Sierra Nevada Mountains. Location Our study was conducted in a 6000‐ha area of contiguous mixed Jeffrey pine‐white fir (Abies concolor Gordon & Glend.) forest on the western slope of the Carson Range on the eastern shore of Lake Tahoe, Nevada. Methods Pre‐EuroAmerican settlement fire regimes (i.e. frequency, return interval, extent, season) were reconstructed in eight contiguous watersheds for a 200‐year period (1650–1850) from fire scars preserved in the annual growth rings of nineteenth century cut stumps and recently dead pre‐settlement Jeffrey pine trees. Superposed epoch analysis (SEA) and correlation analysis were used to examine relationships between tree ring‐based reconstructions of the Palmer Drought Severity Index (PDSI), Southern Oscillation Index (SOI), Pacific Decadal Oscillation (PDO) and pre‐EuroAmerican fire regimes in order to assess the influence of drought and equatorial and north Pacific teleconnections on fire occurrence and fire extent. Results For the entire period of record (1650–1850), wet conditions were characteristic of years without fires. In contrast, fire years were associated with drought. Drought intensity also influenced fire extent and the most widespread fires occurred in the driest years. Years with widespread fires were also preceded by wet conditions 3 years before the fire. Widespread fires were also associated with phase changes of the PDO, with the most widespread burns occurring when the phase changed from warm (positive) to cold (negative) conditions. Annual SOI and fire frequency or extent were not associated in our study. At decadal time scales, burning was more widespread during decades that were dryer and characterized by La Niña and negative PDO conditions. Interannual and interdecadal fire–climate relationships were not stable over time. From 1700 to 1775 there was no interannual relationship between drought, PDO, and fire frequency or extent. However, from 1775 to 1850, widespread fires were associated with dry years preceded by wet years. This period also had the strongest association between fire extent and the PDO. In contrast, fire–climate associations at interdecadal time scales were stronger in the earlier period than in the later period. The change from strong interdecadal to strong interannual climate influence was associated with a breakdown in decadal scale constructive relationships between PDO and SOI. Main conclusions Climate strongly influenced pre‐settlement pine forest fire regimes in northern Sierra Nevada. Both interannual and interdecadal climatic variation regulated conditions conducive to fire activity, and longer term changes in fire frequency and extent correspond with climate‐mediated changes observed in both the northern and southern hemispheres. The sensitivity of fire regimes to shifts in modes of climatic variability suggests that climate was a key regulator of pine forest ecosystem structure and dynamics before EuroAmerican settlement. An understanding of pre‐EuroAmerican fire–climate relationships may provide useful insights into how fire activity in contemporary forests may respond to future climatic variation.     

Wolbachia have made it twice: Hybrid introgression between two sister species of Eurema butterflies

Wolbachia, cytoplasmically inherited endosymbionts of arthropods, are known to hijack their host reproduction in various ways to increase their own vertical transmission. This may lead to the selective sweep of associated mitochondria, which can have a large impact on the evolution of mitochondrial lineages. In Japan, two different Wolbacahia strains (wCI and wFem) are found in two sister species of pierid butterflies, Eurema mandarina and Eurema hecabe. In both species, females infected with wCI (C females) produce offspring with a nearly 1:1 sex ratio, while females infected with both wCI and wFem (CF females) produce all‐female offspring. Previous studies have suggested the historical occurrence of hybrid introgression in C individuals between the two species. Furthermore, hybrid introgression in CF individuals is suggested by the distinct mitochondrial lineages between C females and CF females of E. mandarina. In this study, we performed phylogenetic analyses based on nuclear DNA and mitochondrial DNA markers of E. hecabe with previously published data on E. mandarina. We found that the nuclear DNA of this species significantly diverged from that of E. mandarina. By contrast, mitochondrial DNA haplotypes comprised two clades, mostly reflecting Wolbachia infection status rather than the individual species. Collectively, our results support the previously suggested occurrence of two independent historical events wherein the cytoplasms of CF females and C females moved between E. hecabe and E. mandarina through hybrid introgression.     

Extensive hybridization and past introgression between divergent lineages in a quasi‐clonal hermaphroditic fish: Ramifications for species concepts and taxonomy

Extreme inbreeding is expected to reduce the incidence of hybridization, serving as a prezygotic barrier. Mangrove rivulus is a small killifish that reproduces predominantly by self‐fertilization, producing highly homozygous lines throughout its geographic range. The Bahamas and Caribbean are inhabited by two highly diverged phylogeographic lineages of mangrove rivulus, Kryptolebias marmoratus and a ‘Central clade’ closely related to K. hermaphroditus from Brazil. The two lineages are largely allopatric, but recently were found in syntopy on San Salvador, Bahamas, where a single hybrid was reported. To better characterize the degree of hybridization and the possibility of secondary introgression, here we conducted a detailed genetic analysis of the contact zone on San Salvador. Two mixed populations were identified, one of which contained sexually mature hybrids. The distribution of heterozygosity at diagnostic microsatellite loci in hybrids showed that one of these hybrids was an immediate offspring from the K. marmoratus x Central clade cross, whereas the remaining five hybrids were products of reproduction by self‐fertilization for 1–3 generations following the initial cross. Two hybrids had mitochondrial haplotypes of K. marmoratus and the remaining four hybrids had a haplotype of the Central clade, indicating that crosses go in both directions. In hybrids, alleles of parental lineages were represented in equal proportions suggesting lack of recent backcrossing to either of the parental lineages. However, sympatric populations of two lineages were less diverged than allopatric populations, consistent with introgression. Results are discussed in terms of applicability of the biological species concept for isogenic, effectively clonal, organisms.     

Microsatellite DNA markers for the Chinese wood frog (Rana chensinensis) and tests for their cross-utility in 15 ranid frog species

We developed 22 microsatellite markers for the Chinese wood frog (Rana chensinensis) to study the impact of landscape features on its population structure. Thirty‐four individuals from one breeding site were examined and 14 loci were polymorphic. The number of alleles, expected heterozygosity and observed heterozygosity varied from two to 14, from 0.0833 to 0.9118, and from 0.1376 to 0.8667, respectively. Cross‐species amplification was tested for 15 ranid frog species. The Plateau brown frog, Rana kukunoris (n = 23), was successfully amplified at 18 loci, and 15 were polymorphic with number of alleles varying from two to 18. Ten other species were also amplified at a limited number of loci.     

Hybrid zones between two genetically differentiated forms of the pond frog Rana lessonae in southern Italy

Two distinct population groups of the pond frog Rana lessonae were detected in peninsular Italy and Sicily by multilocus electrophoresis: one group inhabits the peninsula down to northern Calabria, the second occurs in southern Calabria and on Sicily. Fixed alternative alleles distinguish the two groups at 5 of the 25 loci examined; marked allele frequency differences were observed at two additional loci. On average, the two groups differ by a Nei's standard genetic distance of 0.4. A wide hybrid zone (about 120 km) occurs between the two groups, with high genotypic diversity and absence of pure parental genotypes in central Calabria. Patterns of allozyme variation suggest that at least two distinct contact and hybridization events occurred, one in the Catanzaro, the other in the Crati-Sibari plains, about 70 km to the north. Geological evidence indicates that these areas correspond to two main marine-flooded grabens that would have repeatedly interrupted or reduced genetic exchange during Plio-Pleistocene times. The finding of a fixed difference at the Mdhp-1 locus between Sicilian and Calabrian R. lessonae witnesses their continuing differentiation following their last separation by definitive opening of the Strait of Messina, about 50 000 years ago. The wide hybrid zone, the diversity of genotypes and the agreement with Hardy-Weinberg expectations suggest complete hybrid fertility. Different patterns of introgression were observed at the various loci. The pattern of allelic variation at loci in R. lessonae is paralleled by the pattern of variation in lessonae genomes of the sympatric hemiclonal hybrid Rana esculenta, into which lessonae genomes are introduced each generation as a result of hybridogenesis.     

Introgression between invasive and native blue mussels (genus Mytilus) in the central California hybrid zone

The ecological and genetic factors determining the extent of introgression between species in secondary contact zones remain poorly understood. Here, we investigate the relative importance of isolating barriers and the demographic expansion of invasive Mytilus galloprovincialis on the magnitude and the direction of introgression with the native Mytilus trossulus in a hybrid zone in central California. We use double‐digest restriction‐site‐associated DNA sequencing (ddRADseq) to genotype 1337 randomly selected single nucleotide polymorphisms and accurately distinguish early and advanced generation hybrids for the first time in the central California Mytilus spp. hybrid zone. Weak levels of introgression were observed in both directions but were slightly more prevalent from the native M. trossulus into the invasive M. galloprovincialis. Few early and advanced backcrossed individuals were observed across the hybrid zone confirming the presence of strong barriers to interbreeding. Heterogeneous patterns of admixture across the zone of contact were consistent with the colonization history of M. galloprovincialis with more extensive introgression in northern localities furthest away from the putative site of introduction in southern California. These observations reinforce the importance of dynamic spatial and demographic expansions in determining patterns of introgression between close congeners, even in those with high dispersal potential and well‐developed reproductive barriers. Our results suggest that the threat posed by invasive M. galloprovincialis is more ecological than genetic as it has displaced, and continues to displace the native M. trossulus from much of central and southern California.