Phylogeography and systematics of zebra mussels and related species

Evolution of ion transporter Na+/K+-ATPase expression in the osmoregulatory maxillary glands of an invasive copepod

While many freshwater invaders originate from saline habitats, the physiological mechanisms involved are poorly understood. We investigated the evolution of ion transporter Na+/K+-ATPase (NKA) protein expression between ancestral saline and freshwater invading populations of the copepod Eurytemora carolleae (Atlantic clade of the E. affinis complex). We compared in situ NKA expression between populations under common-garden conditions at three salinities in the maxillary glands. We found the evolution of reduced NKA expression in the freshwater population under freshwater conditions and reduced plasticity (canalization) across salinities, relative to the saline population. Our results support the hypothesis that maxillary glands are involved in ion reabsorption from excretory fluids at low-salinity conditions in the saline population. However, mechanisms of freshwater adaptation, such as increased ion uptake from the environment, might reduce the need for ion reabsorption in the freshwater population. These patterns of ion transporter expression contribute insights into the evolution of ionic regulation during habitat change.

Biochemical basis of resistance to multiple contaminations in the native and invasive populations of Dreissena polymorpha

The zebra mussel Dreissena polymorpha (Pallas, 1771) is an invasive species and a valuable bioindicator in the inland waters. Nevertheless, the biochemical reasons for the unique competitiveness of zebra mussels are not clear. This study aimed to compare the native and invasive populations of D. polymorpha in their ability to withstand the same multiple environmental challenges (i.e. chemical: microplastics and caffeine; physical: temperatures). The specimens from the invasive population in west Ukraine (Tn) and native population at lower streams of river Dnipro (south Ukraine, Kh) were sampled in the August of 2021 y. Molluscs from both populations were treated simultaneously with microplastics (MP, 1 mg L-1, pore size 2 μm); caffeine (Caf, 20.0 μg L-1) at 18 °C, elevated temperature (25 °C) and MP and Caf combinations at 18 °C and 25 °C for 14 days. Untreated molluscs exposed at 18 °C represented control groups. A set of the 20 markers of oxidative stress, biotransformation, detoxification and apoptosis were assayed in the total soft tissues. From the two controls, Kh-group indicated lower stressful impact. However, both populations increased caspase-3 and GST activities and lysosomal instability in most exposures, and cholinesterase and phenoloxidase activities under the heating and combine exposures, indicating the remarkable properties to respond to new challenges and synergistic effect of mixtures. Inter-population differences were related to the metallothionein, cathepsin D, cytochrome P540 and oxidative stress responses that reflect population-dependent adverse outcome pathways. The discriminant analysis separated two populations with a substantially lesser magnitude of responses in the invasive population as a sign of higher resistance.

Genome architecture underlying salinity adaptation in the invasive copepod Eurytemora affinis species complex: A review

With climate change, habitat salinity is shifting rapidly throughout the globe. In addition, many destructive freshwater invaders are recent immigrants from saline habitats. Recently, populations of the copepod Eurytemora affinis species complex have invaded freshwater habitats multiple times independently from saline estuaries on three continents. This review discusses features of this species complex that could enhance their evolutionary potential during rapid environmental change. Remarkably, across independent freshwater invasions, natural selection has repeatedly favored the same alleles far more than expected. This high degree of parallelism is surprising, given the expectation of nonparallel evolution for polygenic adaptation. Factors such as population structure and the genome architecture underlying critical traits under selection might help drive rapid adaptation and parallel evolution. Given the preponderance of saline-to-freshwater invasions and climate-induced salinity change, the principles found here could provide invaluable insights into mechanisms operating in other systems and the potential for adaptation in a changing planet.

What we know and don't know about the invasive zebra (Dreissena polymorpha) and quagga (Dreissena rostriformis bugensis) mussels

We summarized existing knowledge on Dreissena polymorpha (the zebra mussel) and D. r. bugensis (the quagga mussel), including data on their taxonomy, systematics, evolution, life cycle, reproduction, feeding, growth and longevity, population dynamics, interspecific competition, habitat requirements, and distribution within and among waterbodies. We analyzed the history of spread of both species and the major pathways and vectors of their spread in Europe and North America. Special consideration was given to their ecological and economic impacts and their natural enemies, like waterfowl, fishes, and parasites, as well as the prevention of their introduction, early detection, control, and containment. We also outlined the most salient ecosystem services provided by zebra and quagga mussels, including water purification, nutrient recycling, culling the effects of eutrophication, biomonitoring, and their role as a food resource for fish and waterfowl. Finally, we identified major knowledge gaps and key studies needed to better understand the biology, ecology, and impacts of these aggressive freshwater invaders. Our review indicates that much crucial information on the quagga mussel is still missing, including key life history parameters, like spawning cues, fecundity, and longevity, particularly for the profundal zone of deep lakes.

Recognizing salinity threats in the climate crisis

Climate change is causing habitat salinity to transform at unprecedented rates across the globe. While much of the research on climate change has focused on rapid shifts in temperature, far less attention has focused on the effects of changes in environmental salinity. Consequently, predictive studies on the physiological, evolutionary, and migratory responses of organisms and populations to the threats of salinity change are relatively lacking. This omission represents a major oversight, given that salinity is among the most important factors that define biogeographic boundaries in aquatic habitats. In this perspective, we briefly touch on responses of organisms and populations to rapid changes in salinity occurring on contemporary time scales. We then discuss factors that might confer resilience to certain taxa, enabling them to survive rapid salinity shifts. Next, we consider approaches for predicting how geographic distributions will shift in response to salinity change. Finally, we identify additional data that are needed to make better predictions in the future. Future studies on climate change should account for the multiple environmental factors that are rapidly changing, especially habitat salinity.

The genome of the zebra mussel, Dreissena polymorpha: a resource for comparative genomics, invasion genetics, and biocontrol

The zebra mussel, Dreissena polymorpha, continues to spread from its native range in Eurasia to Europe and North America, causing billions of dollars in damage and dramatically altering invaded aquatic ecosystems. Despite these impacts, there are few genomic resources for Dreissena or related bivalves. Although the D. polymorpha genome is highly repetitive, we have used a combination of long-read sequencing and Hi-C-based scaffolding to generate a high-quality chromosome-scale genome assembly. Through comparative analysis and transcriptomics experiments, we have gained insights into processes that likely control the invasive success of zebra mussels, including shell formation, synthesis of byssal threads, and thermal tolerance. We identified multiple intact steamer-like elements, a retrotransposon that has been linked to transmissible cancer in marine clams. We also found that D. polymorpha have an unusual 67 kb mitochondrial genome containing numerous tandem repeats, making it the largest observed in Eumetazoa. Together these findings create a rich resource for invasive species research and control efforts.

Advanced molecular-based surveillance of quagga and zebra mussels: A review of environmental DNA/RNA (eDNA/eRNA) studies and considerations for future directions

Sensitive methods, capable of rapidly and accurately detecting aquatic invasive species, are in demand. Molecular-based approaches, such as environmental DNA (eDNA) surveys, satisfy these requirements and have grown in popularity. As such, eDNA surveys could aid the effort to combat the colonisation and spread of two notoriously invasive freshwater mussel species, the quagga mussel (Dreissena rostriformis bugensis) and zebra mussel (D. polymorpha), through improved surveillance ability. Here, we provide a review of dreissenid eDNA literature (both grey and published), summarising efforts involved in the development of various assays for use in multiple different technologies (e.g. quantitative PCR, high-throughput sequencing and loop-mediated isothermal amplification) and sampling scenarios. We discuss important discoveries made along the way, including novel revelations involving environmental RNA (eRNA), as well as the advantages and limitations of available methods and instrumentation. In closing, we highlight critical remaining gaps, where further investigation could lead to advancements in dreissenid monitoring capacity.

Rapid evolution of genome-wide gene expression and plasticity during saline to freshwater invasions by the copepod Eurytemora affinis species complex

Saline migrants into freshwater habitats constitute among the most destructive invaders in aquatic ecosystems throughout the globe. However, the evolutionary and physiological mechanisms underlying such habitat transitions remain poorly understood. To explore the mechanisms of freshwater adaptation and distinguish between adaptive (evolutionary) and acclimatory (plastic) responses to salinity change, we examined genome-wide patterns of gene expression between ancestral saline and derived freshwater populations of the Eurytemora affinis species complex, reared under two different common-garden conditions (0 versus 15 PSU). We found that evolutionary shifts in gene expression (between saline and freshwater inbred lines) showed far greater changes and were more widespread than acclimatory responses to salinity (0 versus 15 PSU). Most notably, 30-40 genes showing evolutionary shifts in gene expression across the salinity boundary were associated with ion transport function, with inorganic cation transmembrane transport forming the largest Gene Ontology category. Of particular interest was the sodium transporter, the Na⁺ /H⁺ antiporter (NHA) gene family, which was discovered in animals relatively recently. Thirty key ion regulatory genes, such as NHA paralogue #7, demonstrated concordant evolutionary and plastic shifts in gene expression, suggesting the evolution of ion transporter function and plasticity during rapid invasions into novel salinities. Moreover, freshwater invasions were associated with the evolution of reduced plasticity in the freshwater population, again for the same key ion transporters, consistent with the predicted evolution of canalization following adaptation to stressful conditions. Our results have important implications for understanding evolutionary and physiological mechanisms of range expansions by some of the most widespread invaders in aquatic habitats.

Differentiation of sympatric zebra and quagga mussels in ecotoxicological studies: A comparison of morphometric data, gene expression, and body metal concentrations

The zebra mussel is among the best studied freshwater molluscs in ecotoxicology, but information on the quagga mussel is lacking. Considering its potential spread, we selected a river in France in which zebra and quagga mussels coexisted, and then we used genetic markers to differentiate the two species and compared morphological parameters. cDNA sequencing assays of ten genes already used in zebra mussels were performed on quagga mussels to obtain functional specific primers. Then we analyzed the expression of genes involved in cellular metabolic activities (Cytochrome-c-oxidase - cox, and ATP synthase - atp), detoxification processes (Glutathione-S-Transferase - gst), oxidative stress (Catalase - cat), and digestive functions (Amylase - amy) on the two species. Whereas morphometric analysis underlined similarities in shape between the two species, relative gene expression profiles and metal concentrations evidenced strong differences. Quagga mussels notably presented half as high concentrations in Cd and Pb, two particularly toxic elements, as zebra mussels. These results imply that i) particular attention should be paid to properly distinguish the two species considering their similar external appearance, and ii) zebra mussels cannot be replaced by quagga mussels in ecotoxicological studies without preliminary investigations on biomarker response patterns. To our knowledge, this study is the first to have undertaken such an approach in gene expression analysis in quagga mussels, and more generally to have compared such biomarker responses of zebra and quagga mussels in the field.

Phylogeography and taxonomic status of trout and salmon from the Ponto-Caspian drainages, with inferences on European Brown Trout evolution and taxonomy

Current taxonomy of western Eurasian trout leaves a number of questions open; it is not clear to what extent some species are distinct genetically and morphologically. The purpose of this paper was to explore phylogeography and species boundaries in freshwater and anadromous trout from the drainages of the Black and the Caspian Seas (Ponto-Caspian). We studied morphology and mitochondrial phylogeny, combining samples from the western Caucasus within the potential range of five nominal species of trout that are thought to inhabit this region, and using the sequences available from GenBank. Our results suggest that the genetic diversity of trout in the Ponto-Caspian region is best explained with the fragmentation of catchments. (1) All trout species from Ponto-Caspian belong to the same mitochondrial clade, separated from the other trout since the Pleistocene; (2) the southeastern Black Sea area is the most likely place of diversification of this clade, which is closely related to the clades from Anatolia; (3) The species from the Black Sea and the Caspian Sea drainages are monophyletic; (4) except for the basal lineage of the Ponto-Caspian clade, Salmo rizeensis, all the lineages produce anadromous forms; (5) genetic diversification within the Ponto-Caspian clade is related to Pleistocene glacial waves; (6) the described morphological differences between the species are not fully diagnostic, and some earlier described differences depend on body size; the differences between freshwater and marine forms exceed those between the different lineages. We suggest a conservative taxonomic approach, using the names S. rizeensis and Salmo labrax for trout from the Black Sea basin and Salmo caspius and Salmo ciscaucasicus for the fish from the Caspian basin.

Evolutionary mechanisms of habitat invasions, using the copepod Eurytemora affinis as a model system

The study of the copepod Eurytemora affinis has provided unprecedented insights into mechanisms of invasive success. In this invited review, I summarize a subset of work from my laboratory to highlight key insights gained from studying E. affinis as a model system. Invasive species with brackish origins are overrepresented in freshwater habitats. The copepod E. affinis is an example of such a brackish invader, and has invaded freshwater habitats multiple times independently in recent years. These invasions were accompanied by the evolution of physiological tolerance and plasticity, increased body fluid regulation, and evolutionary shifts in ion transporter (V‐type H⁺ATPase, Na⁺, K⁺‐ATPase) activity and expression. These evolutionary changes occurred in parallel across independent invasions in nature and in laboratory selection experiments. Selection appears to act on standing genetic variation during invasions, and maintenance of this variation is likely facilitated through ‘beneficial reversal of dominance’ in salinity tolerance across habitats. Expression of critical ion transporters is localized in newly discovered Crusalis leg organs. Increased freshwater tolerance is accompanied by costs to development time and greater requirements for food. High‐food concentration increases low‐salinity tolerance, allowing saline populations to invade freshwater habitats. Mechanisms observed here likely have relevance for other taxa undergoing fundamental niche expansions.

Genetic diversity in introduced golden mussel populations corresponds to vector activity

We explored possible links between vector activity and genetic diversity in introduced populations of Limnoperna fortunei by characterizing the genetic structure in native and introduced ranges in Asia and South America. We surveyed 24 populations: ten in Asia and 14 in South America using the mitochondrial cytochrome c oxidase subunit I (COI) gene, as well as eight polymorphic microsatellite markers. We performed population genetics and phylogenetic analyses to investigate population genetic structure across native and introduced regions. Introduced populations in Asia exhibit higher genetic diversity (H(E) = 0.667-0.746) than those in South America (H(E) =  0.519-0.575), suggesting higher introduction effort for the former populations. We observed pronounced geographical structuring in introduced regions, as indicated by both mitochondrial and nuclear markers based on multiple genetic analyses including pairwise Ф(ST), F(ST), bayesian clustering method, and three-dimensional factorial correspondence analyses. Pairwise F(ST) values within both Asia (F(ST) = 0.017-0.126, P = 0.000-0.009) and South America (F(ST) =0.004-0.107, P = 0.000-0.721) were lower than those between continents (F(ST) = 0.180-0.319, P = 0.000). Fine-scale genetic structuring was also apparent among introduced populations in both Asia and South America, suggesting either multiple introductions of distinct propagules or strong post-introduction selection and demographic stochasticity. Higher genetic diversity in Asia as compared to South America is likely due to more frequent propagule transfers associated with higher shipping activities between source and donor regions within Asia. This study suggests that the intensity of human-mediated introduction vectors influences patterns of genetic diversity in non-indigenous species.

Feasting in fresh water: impacts of food concentration on freshwater tolerance and the evolution of food × salinity response during the expansion from saline into fresh water habitats

Saline to freshwater invasions have become increasingly common in recent years. A key hypothesis is that rates of freshwater invasions have been amplified in recent years by increased food concentration, yet this hypothesis has remained unexplored. We examined whether elevated food concentration could enhance freshwater tolerance, and whether this effect evolves following saline to freshwater invasions. We examined physiological response to salinity and food concentration in a 2 × 2 factorial design, using ancestral brackish and freshwater invading populations of the copepod Eurytemora affinis. We found that high food concentration significantly increases low-salinity tolerance. This effect was reduced in the freshwater population, indicating evolution following the freshwater invasion. Thus, ample food could enable freshwater invasions, allowing subsequent evolution of low-salinity tolerance even under food-poor conditions. We also compared effects of food concentration on freshwater survival between two brackish populations from the native range. Impacts of food concentration on freshwater survival differed between the brackish populations, suggesting variation in functional properties affecting their propensity to invade freshwater habitats. The key implication is that high food concentration could profoundly extend range expansions of brackishwater species into freshwater habitats, potentially allowing for condition-specific competition between saline invaders and resident freshwater species.

Effects of shell morphology on mechanics of zebra and quagga mussel locomotion

Although zebra mussels (Dreissena polymorpha) initially colonized shallow habitats within the North American Great Lakes, quagga mussels (Dreissena bugensis) are becoming dominant in both shallow- and deep-water habitats. Shell morphology differs among zebra, shallow quagga and deep quagga mussels but functional consequences of such differences are unknown. We examined effects of shell morphology on locomotion for the three morphotypes on hard (typical of shallow habitats) and soft (characteristic of deep habitats) sedimentary substrates. We quantified morphology using the polar moment of inertia, a parameter used in calculating kinetic energy that describes shell area distribution and resistance to rotation. We quantified mussel locomotion by determining the ratio of rotational (K(rot)) to translational kinetic energy (K(trans)). On hard substrate, K(rot):K(trans) of deep quagga mussels was fourfold greater than for the other morphotypes, indicating greater energy expenditure in rotation relative to translation. On soft substrate, K(rot):K(trans) of deep quagga mussels was approximately one-third of that on hard substrate, indicating lower energy expenditure in rotation on soft substrate. Overall, our study demonstrates that shell morphology correlates with differences in locomotion (i.e. K(rot):K(trans)) among morphotypes. Although deep quagga mussels were similar to zebra and shallow quagga mussels in terms of energy expenditure on sedimentary substrate, their morphology was energetically maladaptive for linear movement on hard substrate. As quagga mussels can possess two distinct morphotypes (i.e. shallow and deep morphs), they might more effectively utilize a broader range of substrates than zebra mussels, potentially enhancing their ability to colonize a wider range of habitats.

Zebra mussels anchor byssal threads faster and tighter than quagga mussels in flow

While the invasive zebra mussel Dreissena polymorpha has rapidly spread throughout the Great Lakes and inland waterways, it is being displaced by the quagga mussel Dreissena bugensis in shallow water habitats. However, zebra mussels remain dominant in areas with higher water velocity. We hypothesized that the persistence of zebra over quagga mussels in habitats with higher water velocity might result from greater rate and strength of byssal thread attachment. We examined whether zebra mussels relative to quagga mussels have: (1) higher byssal thread synthesis rate, (2) lower dislodgment in flow and (3) greater mechanical force required for detachment from substrate. Specifically, we examined byssal thread synthesis rate and dislodgment of both species in response to water velocities of 0, 50, 100 and 180 cm s–1. Byssal thread synthesis rate was significantly higher for zebra than for quagga mussels at all velocities. Dislodgment from the substrate increased for both species with increasing velocity but was significantly lower for zebra than for quagga mussels. We also tested the mechanical force to detach mussels after short (32 h) and long (two and three months) periods of attachment on hard substrate. Detachment force was significantly higher for zebra than for quagga mussels only after short-term attachment. Higher byssal thread synthesis rate in zebra mussels was a likely factor that minimized their dislodgment in flow and increased short-term attachment strength. Differences in byssal thread synthesis rate between the two species might partly account for the ability of zebra mussels to maintain dominance over quagga mussels in habitats with high velocities.

Evolutionary origins of invasive populations

What factors shape the evolution of invasive populations? Recent theoretical and empirical studies suggest that an evolutionary history of disturbance might be an important factor. This perspective presents hypotheses regarding the impact of disturbance on the evolution of invasive populations, based on a synthesis of the existing literature. Disturbance might select for life-history traits that are favorable for colonizing novel habitats, such as rapid population growth and persistence. Theoretical results suggest that disturbance in the form of fluctuating environments might select for organismal flexibility, or alternatively, the evolution of evolvability. Rapidly fluctuating environments might favor organismal flexibility, such as broad tolerance or plasticity. Alternatively, longer fluctuations or environmental stress might lead to the evolution of evolvability by acting on features of the mutation matrix. Once genetic variance is generated via mutations, temporally fluctuating selection across generations might promote the accumulation and maintenance of genetic variation. Deeper insights into how disturbance in native habitats affects evolutionary and physiological responses of populations would give us greater capacity to predict the populations that are most likely to tolerate or adapt to novel environments during habitat invasions. Moreover, we would gain fundamental insights into the evolutionary origins of invasive populations.

Ancient divisions, recent expansions: phylogeography and population genetics of the round goby Apollonia melanostoma

During the past two decades, the round goby Apollonia melanostoma (=Neogobius melanostomus) has expanded its range via shipping transport and canals, extending north and west from the Ponto-Caspian region of Eurasia and to the North American Great Lakes. Exotic populations of the round goby have been very successful in the Baltic Sea and the Great Lakes regions, exerting significant ecological changes. Our study evaluates the population genetic and biogeographical structure of the round goby across its native and nonindigenous ranges, in light of geological history and its expansion pathways. We analyzed seven new nuclear microsatellite loci and mitochondrial DNA cytochrome b gene sequences from 432 individuals in 22 locations. Population structure was tested using F(ST)-analogs, phylogenetic trees, clustering diagrams, Bayesian assignment tests and nested clade analyses. Results show that native populations in the Black vs. the Caspian Sea basins diverge by 1.4% and c. 350,000 years, corresponding to closure of their prior connections and supporting the taxonomic separation of the Black Sea A. m. melanostoma from the Caspian Sea A. m. affinis. Their within-basin populations diverge by approximately 0.4% and 100,000 years. Nonindigenous populations in the Baltic Sea and Danube and Dnieper Rivers trace to separate northern Black Sea origins, whereas the upper Volga River system houses mixed populations of A. m. melanostoma and A. m. affinis. Native populations average twice the genetic diversity of most exotic sites; however, sites in the Volga River system have high diversity due to mixing of the two taxa. Our results highlight how vicariance and anthropogenic disturbances have shaped a rapidly expanding species' genetic heritage.

Divergence with gene flow between Ponto-Caspian refugia in an anadromous cyprinid Rutilus frisii revealed by multiple gene phylogeography

The Black and Caspian Seas have experienced alternating periods of isolation and interconnection over many Milankovitch climate oscillations and most recently became separated when the meltwater overflow from the Caspian Sea ceased at the end of the last glaciation. Climate-induced habitat changes have indisputably had profound impacts on distribution and demography of aquatic species, yet uncertainties remain about the relative roles of isolation and dispersal in the response of species shared between the Black and Caspian Sea basins. We examined these issues using phylogeographical analysis of an anadromous cyprinid fish Rutilus frisii. Bayesian coalescence analyses of sequence variation at two nuclear and one mitochondrial genes suggest that the Black and Caspian Seas supported separate populations of R. frisii during the last glaciation. Parameter estimates from the fitted isolation-with-migration model showed that their separation was not complete, however, and that the two populations continued to exchange genes in both directions. These analyses also suggested that majority of migrations occurred during the Pleistocene, showing that the variation shared between the Black and Caspian Seas is the result of ancient dispersal along the temporary natural connections between the basins, rather than of incomplete lineage sorting or recent human-mediated dispersal. Gene flow between the refugial populations was therefore an important source of genetic variation, and we suggest that it facilitated the evolutionary response of the populations to changing climate.

Heterogeneity within the native range: population genetic analyses of sympatric invasive and noninvasive clades of the freshwater invading copepod Eurytemora affinis

Invasive species are often composed of highly differentiated populations or sibling species distributed across their native ranges. This study analysed patterns of distribution and the evolutionary and demographic histories of populations within the native range of the copepod species complex Eurytemora affinis. Genetic structure was analysed for samples from 17 locations from both the invaded and native ranges in the St Lawrence River drainage basin, using 652 base pairs of the mitochondrial cytochrome oxidase subunit I gene. This study revealed a high degree of heterogeneity in genetic structure and habitat type in the native range, as well as a bias in the sources of invasive populations. Two genetically distinct clades showed a pattern of niche partitioning within the St Lawrence basin. The noninvasive North Atlantic clade primarily occupied the central portion of the St Lawrence Middle Estuary, whereas the invasive Atlantic clade was more prevalent along the margins, in the upstream reaches of the estuary and downstream salt marshes. Habitat partitioning and genetic subdivision was also present within the Atlantic clade. The freshwater populations were genetically more proximate to the Atlantic clade populations in the estuary than to those in the salt marsh, suggesting the estuary as the source of the invasive populations. The freshwater invading populations showed evidence of a modest population bottleneck. Populations from both clades showed genetic signatures of demographic population expansions that preceded the timing of the last glacial maximum, supporting the St Lawrence as a secondary contact zone between the two clades. Additional analyses on physiological and evolutionary properties of populations in the native range, along with analysis of the selection regime within native habitats, might yield insights into the evolutionary potential to invade.