Bioinvasion threatens the genetic integrity of native diversity and a natural hybrid zone: smooth‐shelled blue mussels (Mytilus spp.) in the Strait of Magellan

Smooth‐shelled blue mussels of the Mytilus edulis species complex are widely distributed bivalve molluscs whose introductions threaten native marine biodiversity (non‐indigenous species – NIS). The aim of the present study was to identify the species and hybrids of Mytilus present in the Magellan Region (southern Chile). Results indicate that three mussel species of the Mytilus edulis complex are found in the region – M. edulis, M. chilensis (or the Southern Hemisphere lineage of Mytilus galloprovincialis), and M. galloprovincialis of Northern Hemisphere origin. For the first time, alleles of the introduced M. trossulus are reported from the Southern Hemisphere. In the Strait of Magellan the native Pacific blue mussel, Mytilus chilensis and the native Atlantic blue mussel, Mytilus edulis, meet and mix at a natural hybrid zone (about 125 km in length). This is the first record of a natural Mytilus hybrid zone in the Southern Hemisphere and is also the first record of the co‐occurrence of genes from all four Mytilus species in any one region. These results contribute to the knowledge of the biodiversity and delimitation of mussel species in southern South America, and highlight how introduced species may threaten the genetic integrity of native species through hybridization and introgression.     

Historical changes in the distributions of invasive and endemic marine invertebrates are contrary to global warming predictions: the effects of decadal climate oscillations

Aim We tested whether a hybrid zone that has formed between an endemic and an invasive species of marine mussel has shifted poleward as expected under a general hypothesis of global warming or has responded instead to decadal climate oscillations. Location We sampled 15 locations on the coast of California, USA, that span the distributions of the two species of marine mussels and their hybrids. Methods Mussels were sampled in 2005–08 and analysed at three nuclear gene loci using methods identical to those used in a study a decade earlier in order to document the genetic architecture of this system. Change in the system was determined by comparing the frequency of species‐specific alleles and multi‐locus genotypes over the intervening decade. Climate variation over the same period was examined by comparing the Pacific Decadal Oscillation (PDO), El Niño/Southern Oscillation (ENSO), upwelling indices and sea surface temperature (SST) during and prior to the study period. Results Contrary to the general expectations of global warming we show that the highly invasive warm‐water mussel Mytilus galloprovincialis and the hybrid zone formed with the endemic species Mytilus trossulus has rapidly contracted southwards. Mytilus galloprovincialis declined in abundance over the northern third of its geographic range (c. 540 km) and has become rare or absent across the northern 200 km of the range it previously colonized during its initial invasion. The distribution of the native species M. trossulus has remained unchanged over the same time period. Main conclusions The large‐scale range shift in the warm‐water invasive species M. galloprovincialis and the hybrid zone it forms with M. trossulus has been exceptionally rapid and is in the opposite direction to that predicted by the global warming hypotheses. This shift, however, is consistent with decadal climate variation associated with the ENSO and the PDO. Since the biogeography of this system was first described in 1999, the PDO has shifted from a warm phase, dominated by frequent and large El Niño events, to a cold‐phase period, with minimal ENSO activity. Thus recent decadal climate variation can oppose global trends in average temperature and this study illustrates the need to integrate the effects of climate change across multiple time‐scales.     

Differential introgression from a sister species explains high FST outlier loci within a mussel species

Scanning genomes for loci with high levels of population differentiation has become a standard of population genetics. F_ST outlier loci are most often interpreted as signatures of local selection, but outliers might arise for many other reasons too often left unexplored. Here, we tried to identify further the history and genetic basis underlying strong differentiation at F_ST outlier loci in a marine mussel. A genome scan of genetic differentiation has been conducted between Atlantic and Mediterranean populations of Mytilus galloprovincialis. The differentiation was low overall (F_ST = 0.03), but seven loci (2%) were strong F_ST outliers. We then analysed DNA sequence polymorphism at two outlier loci. The genetic structure proved to be the consequence of differential introgression of alleles from the sister‐hybridizing species Mytilus edulis. Surprisingly, the Mediterranean population was the most introgressed at these two loci, although the contact zone between the two species is nowadays localized along the Atlantic coasts of France and the British Isles. A historical contact between M. edulis and Mediterranean M. galloprovincialis should have happened during glacial periods. It proved difficult to disentangle two hypotheses: (i) introgression was adaptive, implying edulis alleles have been favoured in Mediterranean populations, or (ii) the genetic architecture of the barrier to edulis gene flow is different between the two M. galloprovincialis backgrounds. Five of the seven outliers between M. galloprovincialis populations were also outliers between M. edulis and Atlantic M. galloprovincialis, which would support the latter hypothesis. Differential introgression across semi‐permeable barriers to gene flow is a neglected scenario to interpret outlying loci that may prove more widespread than anticipated.     

Development of a real-time PCR assay for detection of Mytilus species specific alleles: Application to a sampling survey in Scotland

Shellfish aquaculture is a growing industry in Scotland, dominated by the production of the mussel Mytilus edulis, the native species. Recently the discovery of Mytilus galloprovincialis and Mytilus trossulus together with M. edulis and all 3 hybrids in cultivation in some Scottish sea lochs led to questions regarding the distribution of mussel species in Scotland. The establishment of an extensive sampling survey, involving the collection of mussels at 34 intertidal sites and 10 marinas around Scotland, motivated the development of a high-throughput method for identification of Mytilus alleles from samples. Three Taqman®-MGB probes and one set of primers were designed, based on the previously described Me 15/16 primers targeting the adhesive protein gene sequence, and samples were screened for the presence of M. edulis, M. galloprovincialis and M. trossulus alleles using real-time PCR. Mytilus edulis alleles were identified in samples from all 44 sites. Mytilus galloprovincialis alleles were found together with M. edulis alleles extensively in northern parts of the west and east coasts. Mytilus trossulus alleles were identified in samples from 6 sites in the west and south-west of Scotland. Because M. trossulus is generally undesirable in cultivation and therefore preventing the geographical spread of this species across Scotland is considered beneficial by the shellfish aquaculture industry, these 6 samples were further analysed for genotype frequencies using conventional PCR. Although distribution of the non-native species M. galloprovincialis and M. trossulus have proven to be more widespread than previously thought, there is no evidence from our study of either M. trossulus or M. galloprovincialis acting as an invasive species in Scotland. The real-time PCR method developed in this study has proven to be a rapid and effective tool for the identification of M. edulis, M. galloprovincialis and M. trossulus alleles from samples and should prove useful in future surveys, ecological or aquaculture management related studies in both unispecific and mixed species areas of these species.     

Formation of a recent hybrid zone offers insight into the geographic puzzle and maintenance of species boundaries in musk turtles

Speciation is the result of an accumulation of reproductive barriers between populations, but pinpointing the factors involved is often difficult. However, hybrid zones can form when these barriers are not complete, especially when lineages come into contact in intermediate or modified habitats. We examine a hybrid zone between two closely related riverine turtle species, Sternotherus depressus and S. peltifer, and use dual‐digest RAD sequencing to understand how this hybrid zone formed and elucidate genomic patterns of reproductive isolation. First, the geographical extent and timing of formation of the hybrid zone is established to provide context for understanding the role of extrinsic and intrinsic reproductive isolating mechanisms in this system. The strength of selection on taxon‐specific contributions to maintenance of the hybrid zone is then inferred using a Bayesian genomic cline model. These analyses identify a role for selection inhibiting introgression in some genomic regions at one end of the hybrid zone and promoting introgression in many loci at the other. When selective pressures necessary to generate outliers to the genomic cline are considered with the geographical and temporal context of this hybrid zone, we conclude that habitat‐specific selection probably limits introgression from S. depressus to S. peltifer in the direction of river flow. However, selection is mediating rapid, unidirectional introgression from S. peltifer to S. depressus, which is probably facilitated by anthropogenic habitat alteration. These findings indicate a potentially imminent threat of population‐level genomic extinction for an already imperiled species due to ongoing human‐caused habitat alteration.     

Gene-Flow in a Mosaic Hybrid Zone: Is Local Introgression Adaptive?

Genome-wide scans of genetic differentiation between hybridizing taxa can identify genome regions with unusual rates of introgression. Regions of high differentiation might represent barriers to gene flow, while regions of low differentiation might indicate adaptive introgression—the spread of selectively beneficial alleles between reproductively isolated genetic backgrounds. Here we conduct a scan for unusual patterns of differentiation in a mosaic hybrid zone between two mussel species, Mytilus edulis and M. galloprovincialis. One outlying locus, mac-1, showed a characteristic footprint of local introgression, with abnormally high frequency of edulis-derived alleles in a patch of M. galloprovincialis enclosed within the mosaic zone, but low frequencies outside of the zone. Further analysis of DNA sequences showed that almost all of the edulis allelic diversity had introgressed into the M. galloprovincialis background in this patch. We then used a variety of approaches to test the hypothesis that there had been adaptive introgression at mac-1. Simulations and model fitting with maximum-likelihood and approximate Bayesian computation approaches suggested that adaptive introgression could generate a “soft sweep,” which was qualitatively consistent with our data. Although the migration rate required was high, it was compatible with the functioning of an effective barrier to gene flow as revealed by demographic inferences. As such, adaptive introgression could explain both the reduced intraspecific differentiation around mac-1 and the high diversity of introgressed alleles, although a localized change in barrier strength may also be invoked. Together, our results emphasize the need to account for the complex history of secondary contacts in interpreting outlier loci.     

Mito‐nuclear discordance across a recent contact zone for California voles

To examine the processes that maintain genetic diversity among closely related taxa, we investigated the dynamics of introgression across a contact zone between two lineages of California voles (Microtus californicus). We tested the prediction that introgression of nuclear loci would be greater than that for mitochondrial loci, assuming ongoing gene flow across the contact zone. We also predicted that genomic markers would show a mosaic pattern of differentiation across this zone, consistent with genomes that are semi‐permeable. Using mitochondrial cytochrome b sequences and genome‐wide loci developed via ddRAD‐seq, we analyzed genetic variation for 10 vole populations distributed along the central California coast; this transect included populations from within the distributions of both parental lineages as well as the putative contact zone. Our analyses revealed that (1) the two lineages examined are relatively young, having diverged ca. 8.5–54 kya, (2) voles from the contact zone in Santa Barbara County did not include F1 or early generation backcrossed individuals, and (3) there appeared to be little to no recurrent gene flow across the contact zone. Introgression patterns for mitochondrial and nuclear markers were not concordant; only mitochondrial markers revealed evidence of introgression, putatively due to historical hybridization. These differences in genetic signatures are intriguing given that the contact zone occurs in a region of continuous vole habitat, with no evidence of past or present physical barriers. Future studies that examine specific isolating mechanisms, such as microhabitat use and mate choice, will facilitate our understanding of how genetic boundaries are maintained in this system.     

Influence of genotype and geography on shell shape and morphometric trait variation among North Atlantic blue mussel (Mytilus spp.) populations

The influence of geography and genotype on shell shape (outline) and trait (morphometric) variation among North Atlantic blue mussels and their hybrids has been examined. Shape differences among reference taxa (Mytilus trossulus, Mytilus edulis and Mytilus galloprovincialis) were consistent with an association between taxon‐specific genes and shape genes. Newfoundland M. edulis × M. trossulus populations and northern Quebec M. trossulus populations exhibited an uncoupling of taxon‐specific genes from shape genes, whereas Nova Scotia M. trossulus populations and SW England M. edulis × M. galloprovincialis populations exhibited an association between taxon‐specific genes and shape genes. We found no evidence of a geographic effect (NE versus NW Atlantic) for shape variation, indicating that the genotype effect is stronger than any geographic effect at macrogeographic scales. Pronounced differences were observed in trait variability consistent with an association between taxon‐specific genes and trait genes in European populations, and trait divergence of New York M. edulis from all European mussels. Trait variability in mussels from Newfoundland, Nova Scotia and northern Quebec indicated an uncoupling of taxon genes from trait genes, whereas trait variability in SW England M. edulis × M. galloprovincialis populations was consistent with background genotype, indicating a strong association between taxon genes and trait genes. A pronounced macrogeographic split (NE versus NW Atlantic) regardless of taxonomy was observed, indicating that geography exerts a greater influence than genotype on trait variation at the macrogeographic scale. This is consistent with pronounced within‐taxon genetic divergence, indicative of different selection regimes or more likely of different evolutionary histories of mussels on either side of the North Atlantic. © 2009 The Linnean Society of London, Biological Journal of the Linnean Society, 2009, 96 , 875–897.     

ASYMMETRIC INTROGRESSION OF MITOCHONDRIAL DNA AMONG EUROPEAN POPULATIONS OF BLUE MUSSELS (MYTILUS SPP.)

Abstract.—Mytilus edulis and M. galloprovincialis are two blue mussel species that coexist in western Europe. Previously, we reported that M. galloprovincialis populations contain female and male haplotypes that are fixed in M. edulis populations as well as unique haplotypes. This study assesses whether paraphyly for these species is due to introgression or incomplete lineage extinction. The lineage extinction hypothesis predicts that the shared mtDNA haplotypes in M. galloprovincialis will be significantly diverged from those in M. edulis and form distinct sequence clades. In contrast, the introgression hypothesis proposes that M. edulis haplotypes have only recently been introduced into M. galloprovincialis through hybridization with relatively little divergence accumulating between the shared RFLP haplotypes. We examined 80 mtl6S gene sequences for both the maternal and paternal mtDNA lineages from mussels sampled from various European populations and found strong support for the introgression hypothesis. In addition, we found that M. edulis mtDNA haplotypes appear to be introgressing into mussel populations in the Baltic Sea, which have predominantly M. trossulus nuclear genotypes, indicating that introgressive hybridization is prevalent among European mussel populations.     

Combining hydrodynamic modelling with genetics: can passive larval drift shape the genetic structure of Baltic Mytilus populations?

While secondary contact between Mytilus edulis and Mytilus trossulus in North America results in mosaic hybrid zone formation, both species form a hybrid swarm in the Baltic. Despite pervasive gene flow, Baltic Mytilus species maintain substantial genetic and phenotypic differentiation. Exploring mechanisms underlying the contrasting genetic composition in Baltic Mytilus species will allow insights into processes such as speciation or adaptation to extremely low salinity. Previous studies in the Baltic indicated that only weak interspecific reproductive barriers exist and discussed the putative role of adaptation to environmental conditions. Using a combination of hydrodynamic modelling and multilocus genotyping, we investigate how oceanographic conditions influence passive larval dispersal and hybrid swarm formation in the Baltic. By combining our analyses with previous knowledge, we show a genetic transition of Baltic Mytilus species along longitude 12°‐13°E, that is a virtual line between Malmö (Sweden) and Stralsund (Germany). Although larval transport only occurs over short distances (10–30 km), limited larval dispersal could not explain the position of this genetic transition zone. Instead, the genetic transition zone is located at the area of maximum salinity change (15–10 psu). Thus, we argue that selection results in weak reproductive barriers and local adaptation. This scenario could maintain genetic and phenotypic differences between Baltic Mytilus species despite pervasive introgressive hybridization.     

Larval settlement can explain the adult distribution of Mytilus californianus Conrad but not of M. galloprovincialis Lamarck or M. trossulus Gould in Moss Landing, central California: Evidence from genetic identification of spat

We investigated the spatial distribution of adult and newly settled mussels (Mytilus galloprovincialis Lamarck, Mytilus trossulus Gould and Mytilus californianus Conrad) on the shore at Moss Landing, California to test the hypothesis that adult distributions are a result of settlement patterns. Adult M. californianus were most abundant on a wave-exposed rocky jetty and adults of Blue mussels (M. trossulus and M. galloprovincialis) were more abundant inside the protected Moss Landing harbor. Using taxon-specific polymerase chain reactions, we monitored recruitment during continuous 1-2 week intervals on fibrous scrubbing pads for 12 months in 2002-2003. All mussel species settled in greatest numbers on the exposed jetty, and Blue mussels settled in greater numbers there than did M. californianus. Because Blue mussels settled abundantly where their adults were rare, post-settlement mortality appeared to be the strongest influence on adult distribution. In contrast, M. californianus settled mostly in their adult habitat.     

Glacial history of the European marine mussels Mytilus,inferred from distribution of mitochondrial DNA lineages

Mussels of the genus Mytilus have been used to assess the circumglacial phylogeography of the intertidal zone. These mussels are representative components of the intertidal zone and have rapidly evolving mitochondrial DNA, suitable for high resolution phylogeographic analyses. In Europe, the three Mytilus species currently share mitochondrial haplotypes, owing to the cases of extensive genetic introgression. Genetic diversity of Mytilus edulis, Mytilus trossulus and Mytilus galloprovincialis was studied using a 900-bp long part of the most variable fragment of the control region from one of their two mitochondrial genomes. To this end, 985 specimens were sampled along the European coasts, at sites ranging from the Black Sea to the White Sea. The relevant DNA fragments were amplified, sequenced and analyzed. Contrary to the earlier findings, our coalescence and nested cladistics results show that only a single M. edulis glacial refugium existed in the Atlantic. Despite that, the species survived the glaciation retaining much of its diversity. Unsurprisingly, M. galloprovincialis survived in the Mediterranean Sea. In a relatively short time period, around the climatic optimum at 10 ky ago, the species underwent rapid expansion coupled with population differentiation. Following the expansion, further contemporary gene flow between populations was limited.     

A population genetic study of the hybrid zone of Mytilus trossulus Gould, 1850 and an introduced species,M. galloprovincialis Lamarck, 1819, (Bivalvia: Mytilidae) in peter the great bay in the Sea of Japan

This paper examines the genetic variability of the Pacific mussel Mytilus trossulus and an introduced Atlantic species, M. galloprovincialis, in the northwestern Sea of Japan (Peter the Great Bay and Kievka Bay). The genotyping of individuals from eight populations was carried out using eight polymorphic enzyme loci and two nuclear DNA markers (Me-5 and ITS-1,2); the occurrence frequency of parent species and their hybrids was determined. The enzyme and nuclear markers demonstrated concordant genetic variation. The genotypes of the native species M. trossulus were predominant in the samples studied. The frequency of the introduced species M. galloprovincialis in the total material was relatively low; however, it reached 42 ± 2% in samples that were collected in Possjet Bay near the town of Zarubino in a zone of active international navigation. In this area the greatest number of hybrids was found as well. It is concluded that the invasion of M. galloprovincialis in the northwestern Sea of Japan is continuing; permanent populations of this mussel appeared in Possjet Bay that were not recorded here previously.     

Study of Two Species of Mussels, Mytilus trossulusand Mytilus galloprovincialis(Bivalvia, Mytilidae), and Their Hybrids in Peter the Great Bay of the Sea of Japan with the Use of PCR Markers

Bivalve mollusks of the genus Mytilus(M. trossulusand M. galloprovincialis) occurring in Peter the Great Bay of the Sea of Japan were first studied in Russia. A region of nonrepetitive sequences of the gene encoding the polyphenolic adhesive protein bissus was used as a species-specific genetic marker. After amplification using specific primers, a 126-bp fragment was found to amplify in all representatives ofM. galloprovincialiscollected from driftwood in the gulf Posset (the southwestern part of Peter the Great Bay). M. trossulusspecimens from the same region were shown to have a 168-bp fragment. In Vostok Gulf (the eastern part of Peter the Great Bay), both artificially grown mussels and those from natural habitats contained a 168-bp fragment or two fragments (126- and 168-bp) that corresponded to a hybrid form between the above species. The possibility of using this genetic marker to identify closely related Mytilusstrains and their hybrids in similar habitats, near the Primorye coast in particular, was demonstrated. The presence of approximately 9% of hybrid specimens confirms that a zone of hybridization between M. trossulusand M. galloprovincialismay exist in this region.     

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.     

Genetic diversity of Southern hemisphere blue mussels (Bivalvia: Mytilidae) and the identification of non‐indigenous taxa

The taxonomic and evolutionary affinities of Southern hemisphere smooth‐shelled blue mussels are unclear, with studies using different marker types having identified different relationships among various geographic regions. Using an existing and a new molecular assay, the present study builds on previous work to test the distribution of blue mussels native to and introduced to the Southern hemisphere. Populations of Mytilus were sampled from New Zealand, Australia, and Chile. The nuclear‐DNA marker Me 15/16 was used to identify the taxonomic status of 484 individuals. A new restriction fragment length polymorphism (RFLP) assay was used to identify the hemisphere of origin for a subset of Mytilus galloprovincialis. The Me15/16 marker identified 478 pure M. galloprovincialis from Southern hemisphere sites and six Mytilus edulis/M. galloprovincialis hybrids from the Auckland Islands (New Zealand) and Chile. A cytoplasmic RFLP identified Northern hemisphere M. galloprovincialis in almost every Southern hemisphere region. The presence of native M. galloprovincialis at high latitudes (up to 52°S) has implications for our understanding of environmentally induced selective constraints considered to determine species distributions. Widespread occurrence of invasive Northern hemisphere blue mussels in the Southern hemisphere is documented for the first time. Identification of inter‐specific hybrids (M. edulis × M. galloprovincialis) in Chile and in the Auckland Islands (subantarctic New Zealand) illustrates that environments ranging from international ports to remote protected locations are vulnerable to bioinvasion. © 2010 The Linnean Society of London, Biological Journal of the Linnean Society, 2010, 101 , 898–909.     

Barriers to hybridisation and their conservation implications for a highly threatened Australian fish species

Hybridisation and introgression are natural phenomena that may lead to the transfer of adaptive alleles from one species to another and increased species diversity. At the same time, hybridisation and subsequent introgression threaten many species world‐wide through the loss of genetic and species diversity. In Australia, introgressive hybridisation between native and alien species has not typically been considered a significant threat to native biodiversity because of the taxonomic distance between native and alien biota. However, many native fish have been introduced outside their natural range. Recently, four taxa in the genus Melanotaenia have been nationally listed as threatened due to introgressive hybridisation with introduced Melanotaenia splendida. We examined pre‐ and post‐zygotic barriers to hybridisation between M. splendida and one of these threatened taxa—Running River rainbowfish (RRR)—to assess the potential for hybridisation to occur. We used dichotomous mate choice experiments to examine pre‐zygotic barriers and mating experiments to examine post‐zygotic barriers. Size was not a significant predictor of the proportion of time subjects spent with a potential mate, nor was there any significant difference in the amount of time subjects spent with potential mates of their own or the opposite species. Eggs from hybrid pairings with female RRR had a slightly higher hatching rate than those from hybrid pairings with female M. splendida, but neither were significantly different from intraspecies crosses. We could not identify any definite barriers to hybridisation, demonstrating that the introduction of “native” fish species outside their natural range poses a higher risk of hybridisation than previously thought. We call for better education around the consequences of moving “native” fish and the development of rapid response plans to deal with recently established alien populations of Australian fish species in order to prevent future extinctions due to introgressive hybridisation.