PUMPING IONS: RAPID PARALLEL EVOLUTION OF IONIC REGULATION FOLLOWING HABITAT INVASIONS

The Genome Architecture of the Copepod Eurytemora carolleeae - the Highly Invasive Atlantic Clade of the Eurytemoraaffinis Species Complex

Copepods are among the most abundant organisms on the planet and play critical functions in aquatic ecosystems. Among copepods, populations of the Eurytemora affinis species complex are numerically dominant in many coastal habitats and serve as food sources for major fisheries. Intriguingly, certain populations possess the unusual capacity to invade novel salinities on rapid time scales. Despite their ecological importance, high-quality genomic resources have been absent for calanoid copepods, limiting our ability to comprehensively dissect the genome architecture underlying the highly invasive and adaptive capacity of certain populations. Here, we present the first chromosome-level genome of a calanoid copepod, from the Atlantic clade (Eurytemora carolleeae) of the E. affinis species complex. This genome was assembled using high-coverage PacBio long-read and Hi-C sequences of an inbred line, generated through 30 generations of full-sib mating. This genome, consisting of 529.3 Mb (contig N50 = 4.2 Mb, scaffold N50 = 140.6 Mb), was anchored onto four chromosomes. Genome annotation predicted 20,262 protein-coding genes, of which ion transport-related gene families were substantially expanded based on comparative analyses of 12 additional arthropod genomes. Also, we found genome-wide signatures of historical gene body methylation of the ion transport-related genes and the significant clustering of these genes on each chromosome. This genome represents one of the most contiguous copepod genomes to date and is among the highest quality marine invertebrate genomes. As such, this genome provides an invaluable resource to help yield fundamental insights into the ability of this copepod to adapt to rapidly changing environments.

Assessment of the physiological performance of the invasive oriental shrimp Palaemon macrodactylus from an atypical marine population

Since 2000, a well-established population of the invasive oriental shrimp Palaemon macrodactylus has been present in fully marine conditions in the southwestern Atlantic Ocean (~38° S). To assess the physiological performance of this atypical population restricted to fully marine conditions, we conducted a laboratory experiment in which individuals were transferred from 35 ‰S (local seawater) to 2 ‰S; 5 ‰S; 10 ‰S; 20 ‰S; 50 ‰S and 60‰ for short (6 h), medium (48 h), and long (>504 h) acclimation periods. We measured the time course response of relevant parameters in the shrimp's hemolymph; activity of Na+, K+-ATPase (NKA), and V-H+-ATPase (VHA); and muscle water content. Shrimp showed great osmoregulatory plasticity, being able to survive for long periods between 5 ‰S and 50 ‰S, whereas no individual survived after transfer to either 2 ‰S or 60 ‰S. Shrimp hyper-regulated hemolymph osmolality at 5 ‰S and 10 ‰S, hypo-regulated at 35 ‰S and 50 ‰S, and isosmoticity was close to 20 ‰S. Compared to 35 ‰S, prolonged acclimation to 5 ‰S caused a decrease in hemolymph osmolality (~34%) along with sodium and chloride concentrations (~24%); the NKA and VHA activities decreased by ~52% and ~88%, respectively, while muscle water content was tightly regulated. Our results showed that the atypical population of P. macrodactylus studied here lives in a chronic hypo-osmo-ion regulatory state and suggest that fully marine conditions contribute to its poor performance at the lower limit of salinity tolerance (<5 ‰S).

Ionic regulatory strategies of crabs: the transition from water to land

Terrestrial crabs (brachyurans and anomurans) have invaded land following a variety of pathways from marine and/or via freshwater environments. This transition from water to land requires physiological, ecological, and behavioral adaptations to allow the exploitation of these new environmental conditions. Arguably, the management of salt and water balance (e.g., osmoregulation) is integral for their survival and success in an environment where predominantly low-salinity aquatic (e.g., freshwater) water sources are found, sometimes in only minimal amounts. This requires a suite of morphological and biochemical modifications, especially at the branchial chamber of semi-terrestrial and terrestrial crabs to allow reprocessing of urine to maximize ion uptake. Using knowledge gained from electrophysiology, biochemistry, and more recent molecular biology techniques, we present summarized updated models for ion transport for all major taxonomic groups of terrestrial crabs. This is an exciting and fast-moving field of research, and we hope that this review will stimulate further study. Terrestrial crabs retain their crown as the ideal model group for studying the evolutionary pathways that facilitated terrestrial invasion.

Rethinking ecological niches and geographic distributions in face of pervasive human influence in the Anthropocene

Species are distributed in predictable ways in geographic spaces. The three principal factors that determine geographic distributions of species are biotic interactions (B), abiotic conditions (A), and dispersal ability or mobility (M). A species is expected to be present in areas that are accessible to it and that contain suitable sets of abiotic and biotic conditions for it to persist. A species' probability of presence can be quantified as a combination of responses to B, A, and M via ecological niche modeling (ENM; also frequently referred to as species distribution modeling or SDM). This analytical approach has been used broadly in ecology and biogeography, as well as in conservation planning and decision-making, but commonly in the context of 'natural' settings. However, it is increasingly recognized that human impacts, including changes in climate, land cover, and ecosystem function, greatly influence species' geographic ranges. In this light, historical distinctions between natural and anthropogenic factors have become blurred, and a coupled human-natural landscape is recognized as the new norm. Therefore, B, A, and M (BAM) factors need to be reconsidered to understand and quantify species' distributions in a world with a pervasive signature of human impacts. Here, we present a framework, termed human-influenced BAM (Hi-BAM, for distributional ecology that (i) conceptualizes human impacts in the form of six drivers, and (ii) synthesizes previous studies to show how each driver modifies the natural BAM and species' distributions. Given the importance and prevalence of human impacts on species distributions globally, we also discuss implications of this framework for ENM/SDM methods, and explore strategies by which to incorporate increasing human impacts in the methodology. Human impacts are redefining biogeographic patterns; as such, future studies should incorporate signals of human impacts integrally in modeling and forecasting species' distributions.

Paralog switching facilitates diadromy: ontogenetic, microevolutionary and macroevolutionary evidence

Identifying how the demands of migration are met at the level of gene expression is critical for understanding migratory physiology and can potentially reveal how migratory forms evolve from nonmigratory forms and vice versa. Among fishes, migration between freshwater and seawater (diadromy) requires considerable osmoregulatory adjustments, powered by the ion pump Na+, K+-ATPase (NKA) in the gills. Paralogs of the catalytic α-subunit of the pump (NKA α1a and α1b) are reciprocally upregulated in fresh- and seawater, a response known as paralog-switching, in gills of some diadromous species. We tested ontogenetic changes in NKA α-subunit paralog expression patterns, comparing pre-migrant and migrant alewife (Alosa pseudoharengus) sampled in their natal freshwater environment and after 24 h in seawater. In comparison to pre-migrants, juvenile out-migrants exhibited stronger paralog switching via greater downregulation of NKA α1a in seawater. We also tested microevolutionary changes in the response, exposing juvenile diadromous and landlocked alewife to freshwater (0 ppt) and seawater (30 ppt) for 2, 5, and 15 days. Diadromous and landlocked alewife exhibited salinity-dependent paralog switching, but levels of NKA α1b transcription were higher and the decrease in NKA α1a was greater after seawater exposure in diadromous alewife. Finally, we placed alewife α-subunit NKA paralogs in a macroevolutionary context. Molecular phylogenies show alewife paralogs originated independently of paralogs in salmonids and other teleosts. This study demonstrated that NKA paralog switching is tied to halohabitat profile and that duplications of the NKA gene provided the substrate for multiple, independent molecular solutions that support a diadromous life history.

Genome structural variants shape adaptive success of an invasive urban malaria vector Anopheles stephensi

Global changes are associated with the emergence of several invasive species. However, the genomic determinants of the adaptive success of an invasive species in a new environment remain poorly understood. Genomic structural variants (SVs), consisting of copy number variants, play an important role in adaptation. SVs often cause large adaptive shifts in ecologically important traits, which makes SVs compelling candidates for driving rapid adaptations to environmental changes, which is critical to invasive success. To address this problem, we investigated the role SVs play in the adaptive success of Anopheles stephensi , a primary vector of urban malaria in South Asia and an invasive malaria vector in several South Asian islands and Africa. We collected whole genome sequencing data from 115 mosquitoes from invasive island populations and four locations from mainland India, an ancestral range for the species. We identified 2,988 duplication copy number variants and 16,038 deletions in these strains, with ∼50% overlapping genes. SVs are enriched in genomic regions with signatures of selective sweeps in the mainland and invasive island populations, implying a putative adaptive role of SVs. Nearly all high-frequency SVs, including the candidate adaptive variants, in the invasive island populations are present on the mainland, suggesting a major contribution of existing variation to the success of the island populations. Among the candidate adaptive SVs, three duplications involving toxin-resistance genes evolved, likely due to the widespread application of insecticides in India since the 1950s. We also identify two SVs associated with the adaptation of An. stephensi larvae to brackish water in the island and two coastal mainland populations, where the mutations likely originated. Our results suggest that existing SVs play a vital role in the evolutionary success of An. stephensi in new environmental conditions.

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.

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.

Early onset of urea synthesis and ammonia detoxification pathways in three terrestrially developing frogs

Frogs evolved terrestrial development multiple times, necessitating mechanisms to avoid ammonia toxicity at early stages. Urea synthesis from ammonia is a key adaptation that reduces water dependence after metamorphosis. We tested for early expression and plasticity of enzymatic mechanisms of ammonia detoxification in three terrestrial-breeding frogs: foam-nest-dwelling larvae of Leptodactylus fragilis (Lf) and arboreal embryos of Hyalinobatrachium fleischmanni (Hf) and Agalychnis callidryas (Ac). Activity of two ornithine-urea cycle (OUC) enzymes, arginase and CPSase, and levels of their products urea and CP in tissues were high in Lf regardless of nest hydration, but reduced in experimental low- vs. high-ammonia environments. High OUC activity in wet and dry nests, comparable to that under experimental high ammonia, suggests terrestrial Lf larvae maintain high capacity for urea excretion regardless of their immediate risk of ammonia toxicity. This may aid survival through unpredictably long waiting periods before rain enables their transition to water. Moderate levels of urea and CP were present in Hf and Ac tissues and enzymatic activities were lower than in Lf. In both species, embryos in drying clutches can hatch and enter the water early, behaviorally avoiding ammonia toxicity. Moreover, glutamine synthetase was active in early stages of all three species, condensing ammonia and glutamate to glutamine as another mechanism of detoxification. Enzyme activity appeared highest in Lf, although substrate and product levels were higher in Ac and Lf. Our results reveal that multiple biochemical mechanisms of ammonia detoxification occur in early life stages of anuran lineages that evolved terrestrial development.

Mechanisms of Na+ uptake from freshwater habitats in animals

Life in fresh water is osmotically and energetically challenging for living organisms, requiring increases in ion uptake from dilute environments. However, mechanisms of ion uptake from freshwater environments are still poorly understood and controversial, especially in arthropods, for which several hypothetical models have been proposed based on incomplete data. One compelling model involves the proton pump V-type H+ ATPase (VHA), which energizes the apical membrane, enabling the uptake of Na+ (and other cations) via an unknown Na+ transporter (referred to as the "Wieczorek Exchanger" in insects). What evidence exists for this model of ion uptake and what is this mystery exchanger or channel that cooperates with VHA? We present results from studies that explore this question in crustaceans, insects, and teleost fish. We argue that the Na+/H+ antiporter (NHA) is a likely candidate for the Wieczorek Exchanger in many crustaceans and insects; although, there is no evidence that this is the case for fish. NHA was discovered relatively recently in animals and its functions have not been well characterized. Teleost fish exhibit redundancy of Na+ uptake pathways at the gill level, performed by different ion transporter paralogs in diverse cell types, apparently enabling tolerance of low environmental salinity and various pH levels. We argue that much more research is needed on overall mechanisms of ion uptake from freshwater habitats, especially on NHA and other potential Wieczorek Exchangers. Such insights gained would contribute greatly to our general understanding of ionic regulation in diverse species across habitats.

Experimental evolution reveals the synergistic genomic mechanisms of adaptation to ocean warming and acidification in a marine copepod

Metazoan adaptation to global change relies on selection of standing genetic variation. Determining the extent to which this variation exists in natural populations, particularly for responses to simultaneous stressors, is essential to make accurate predictions for persistence in future conditions. Here, we identified the genetic variation enabling the copepod Acartia tonsa to adapt to experimental ocean warming, acidification, and combined ocean warming and acidification (OWA) over 25 generations of continual selection. Replicate populations showed a consistent polygenic response to each condition, targeting an array of adaptive mechanisms including cellular homeostasis, development, and stress response. We used a genome-wide covariance approach to partition the allelic changes into three categories: selection, drift and replicate-specific selection, and laboratory adaptation responses. The majority of allele frequency change in warming (57%) and OWA (63%) was driven by shared selection pressures across replicates, but this effect was weaker under acidification alone (20%). OWA and warming shared 37% of their response to selection but OWA and acidification shared just 1%, indicating that warming is the dominant driver of selection in OWA. Despite the dominance of warming, the interaction with acidification was still critical as the OWA selection response was highly synergistic with 47% of the allelic selection response unique from either individual treatment. These results disentangle how genomic targets of selection differ between single and multiple stressors and demonstrate the complexity that nonadditive multiple stressors will contribute to predictions of adaptation to complex environmental shifts caused by global change.

Freshwater adaptation in prickly sculpin (Pisces: Cottidae): intraspecific comparisons reveal evidence for water pH and Na+ concentration driving diversity in gill H+-ATPase and ionoregulation

Phenotypic divergence is a hallmark of adaptive radiation. One example involves differentiation in physiological traits involved in ion regulation among species with contrasting life-styles and living in distinct environments. Differentiation in ion regulation and its ecological implications among populations within species are, however, less well understood. To address this knowledge gap, we collected prickly sculpin (Cottus asper) from distinct habitat types including coastal rivers connected to estuaries, coastal lakes, and interior lakes, all from British Columbia, Canada. We tested for differences in plasma Na+ and Cl-, gill Na+/K+-ATPase and H+-ATPase activities and protein abundance as well as changes in body mass, and arterial blood pH in fish sampled from the field and acclimated to two different freshwater conditions in the laboratory including artificial lake water (ALW) and ion-poor water (IPW). We also tested for associations between environmental water chemistry and the physiological characteristics associated with ion regulation. Transfer to IPW resulted upregulation in gill Na+/K+-ATPase and H+-ATPase activities as well as increases in gill H+-ATPase protein expression level in each habitat compared with the common ALW treatment. Despite the presence of population-within-habitat type differences, significant habitat-type effects were revealed in most of the ion regulation characteristics examined under different acclimation conditions. Significantly lower plasma Cl- was detected in fish from coastal rivers compared to fish from the other two habitat types during the IPW treatment, which was also significantly lower compared with ALW. Similarly, gill Na+/K+-ATPase activity was lower in the coastal river populations in IPW than in fish from coastal and interior lakes, which was not in accordance with the protein expressed in the gill. For gill H+-ATPase, fish from interior lake populations had the highest level of activity across all habitat types under all conditions, which was related to the protein levels in the gill. The activity of gill H+-ATPase was positively correlated with the combined effect of water Na+ and pH under the ALW treatment. Our results suggest that variation in habitat may be an important factor driving differences in gill Na+/K+-ATPase and H+-ATPase activities across populations of C. asper. Further, the combined effect of water Na+ and pH may have played a key role in physiological adaptation in C. asper during post-glacial freshwater colonization and dispersal.

The genome-wide identification and adaptive evolution of slc9 genes in Leuciscus waleckii under extremely alkaline conditions

The solute carrier family 9 (slc9) genes, especially slc9a isoform coding proteins contribute to electroneutral countertransport of H⁺ for Na⁺ across the plasmalemmal and organellar membranes, intracellular pH and cellular volume regulation as well as the electrolyte, acid-base, and fluid volume homeostasis at the systemic level. These functional properties determine a potential basis for organisms to challenge stressful conditions. However, these well-done researches have been reported more in mammals. Thus, in this study, a total of eleven slc9 genes were identified from the latest version genome of L. waleckii, a cyprinid fish that could tolerate extremely alkaline environments (pH 9.6). The evolutionary footprint of slc9 genes was uncovered via the analysis of copy numbers, gene structure, motif composition, chromosome location and phylogenetic relationship. More importantly, there were two SNPs located on 5' UTR and three non-synonymous mutations in the coding region of the slc9a3.2 gene by comparing freshwater with alkaline water populations attached to resequencing technology. Slc9a3.2 gene was a statistically significant low expression in gill tissue with extremely alkaline pressure. Generally, slc9 gene family in L. waleckii was highly conserved. Several important SNPs with high Fst values were identified where non-synonymous mutations occurred between freshwater and alkaline water populations, and they may play an important role in specific functional differentiation. Slc9 genes had clear tissue expression preferences and were involved in abiotic stress response, indicating their roles in physiological function and strong self-regulating capacity. Our insight into the genetic variations that take place in the individual genes under extreme conditions could provide a feasible example for studying specific molecular mechanisms based on genomic data with increasing environmental stress.

Genome-wide signatures of synergistic epistasis during parallel adaptation in a Baltic Sea copepod

The role of epistasis in driving adaptation has remained an unresolved problem dating back to the Evolutionary Synthesis. In particular, whether epistatic interactions among genes could promote parallel evolution remains unexplored. To address this problem, we employ an Evolve and Resequence (E&R) experiment, using the copepod Eurytemora affinis, to elucidate the evolutionary genomic response to rapid salinity decline. Rapid declines in coastal salinity at high latitudes are a predicted consequence of global climate change. Based on time-resolved pooled whole-genome sequencing, we uncover a remarkably parallel, polygenic response across ten replicate selection lines, with 79.4% of selected alleles shared between lines by the tenth generation of natural selection. Using extensive computer simulations of our experiment conditions, we find that this polygenic parallelism is consistent with positive synergistic epistasis among alleles, far more so than other mechanisms tested. Our study provides experimental and theoretical support for a novel mechanism promoting repeatable polygenic adaptation, a phenomenon that may be common for selection on complex physiological traits.

Repeated Genetic Targets of Natural Selection Underlying Adaptation of Fishes to Changing Salinity

Ecological transitions across salinity boundaries have led to some of the most important diversification events in the animal kingdom, especially among fishes. Adaptations accompanying such transitions include changes in morphology, diet, whole-organism performance, and osmoregulatory function, which may be particularly prominent since divergent salinity regimes make opposing demands on systems that maintain ion and water balance. Research in the last decade has focused on the genetic targets underlying such adaptations, most notably by comparing populations of species that are distributed across salinity boundaries. Here, we synthesize research on the targets of natural selection using whole-genome approaches, with a particular emphasis on the osmoregulatory system. Given the complex, integrated and polygenic nature of this system, we expected that signatures of natural selection would span numerous genes across functional levels of osmoregulation, especially salinity sensing, hormonal control, and cellular ion exchange mechanisms. We find support for this prediction: genes coding for V-type, Ca2+, and Na+/K+-ATPases, which are key cellular ion exchange enzymes, are especially common targets of selection in species from six orders of fishes. This indicates that while polygenic selection contributes to adaptation across salinity boundaries, changes in ATPase enzymes may be of particular importance in supporting such transitions.

A Multi-Species Comparison and Evolutionary Perspectives on Ion Regulation in the Antennal Gland of Brachyurans

Brachyurans inhabit a variety of habitats and have evolved diverse osmoregulatory patterns. Gills, antennal glands and a lung-like structure are important organs of crabs that maintain their homeostasis in different habitats. Species use different processes to regulate ions in the antennal gland, especially those with high terrestriality such as Grapsoidea and Ocypodoidea. Our phylogenetic generalized least square (PGLS) result also suggested that there is a correlation between antennal gland NKA activity and urine-hemolymph ratio for Na⁺ concentration in hypo-osmotic environments among crabs. Species with higher antennal gland NKA activity showed a lower urine-hemolymph ratio for Na⁺ concentration under hypo-osmotic stress. These phenomenon may correlate to the structural and functional differences in gills and lung-like structure among crabs. However, a limited number of studies have focused on the structural and functional differences in the antennal gland among brachyurans. Integrative and systemic methods like next generation sequencing and proteomics method can be useful for investigating the differences in multi-gene expression and sequences among species. These perspectives can be combined to further elucidate the phylogenetic history of crab antennal glands.

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.

Expression levels and activities of energy-yielding ATPases in the oligohaline neritid snail Theodoxus fluviatilis under changing environmental salinities

The aquatic gastropod Theodoxus fluviatilis occurs in Europe and adjacent areas of Asia. The snail species has formed two genetically closely related subgroups, the freshwater ecotype (FW) and the brackish water ecotype (BW). Other than individuals of the FW ecotype, those of the BW ecotype survive in salinities of up to 28‰. Coastal aquatic ecosystems may be affected by climate change due to salinization. Thus, we investigated how the two Theodoxus ecotypes adjust to changes in environmental salinity, focusing on the question whether Na+/K+-ATPase or V-ATPase are regulated on the transcriptional, the translational or at the activity level under changing external salinities. Animals were gradually adjusted to extreme salinities in containers under long-day conditions and constant temperature. Whole body RNA- or protein extracts were prepared. Semi-quantitative PCR- and western blot-analyses did not reveal major changes in transcript or protein abundances for the two transporters under low or high salinity conditions. No significant changes in ATPase activities in whole body extracts of animals adjusted to high or low salinity conditions were detected. We conclude that constitutive expression of ATPases is sufficient to support osmotic and ion regulation in this species under changing salinities given the high level of tolerance with respect to changes in body fluid volume.

Ion Transporter Gene Families as Physiological Targets of Natural Selection During Salinity Transitions in a Copepod

Salinity is a key factor that structures biodiversity on the planet. With anthropogenic change, such as climate change and species invasions, many populations are facing rapid and dramatic changes in salinity throughout the globe. Studies on the copepod Eurytemora affinis species complex have implicated ion transporter gene families as major loci contributing to salinity adaptation during freshwater invasions. Laboratory experiments and population genomic surveys of wild populations have revealed evolutionary shifts in genome-wide gene expression and parallel genomic signatures of natural selection during independent salinity transitions. Our results suggest that balancing selection in the native range and epistatic interactions among specific ion transporter paralogs could contribute to parallel freshwater adaptation. Overall, these studies provide unprecedented insights into evolutionary mechanisms underlying physiological adaptation during rapid salinity change.

Effects of Bicarbonate Stress on Serum Ions and Gill Transporters in Alkali and Freshwater Forms of Amur Ide (Leuciscus waleckii)

The Amur ide (Leuciscus waleckii) is a fish in the Cyprinidae family. Compared with other Amur ide living in freshwater ecosystems, the Amur ide population in Lake Dali Nor of China is famous for its high tolerance to the alkaline conditions of 54 mM (pH 9.6). Yet, surprisingly, the ionoregulatory mechanism responsible for this remarkable alkaline adaptation remains unclear. Therefore, this study sought to investigate how bicarbonate affects the acid-base balancing and ionoregulatory responses of this animal. Here, using a comparative approach, the alkali form of Amur ide and its ancestral freshwater form living in other freshwater basins were each exposed to 50 mM (pH 9.59 ± 0.09), a level close to the alkalinity of Lake Dali Nor, and their physiological (AE1) adjustment of ions and acid-base regulation were investigated. This study highlighted differences in blood pH and serum ions (e.g., Na⁺, K⁺, Cl⁻, and Ca²⁺), Na⁺/K⁺ ATPase (NKA) activity and its mRNA level, and mRNA expression of gill transporters (Na⁺/H⁺ exchanger member 2 and/or 3, Na⁺/HCO3- cotransporter (NBC1), Cl⁻/HCO3- exchanger, Na⁺/Cl⁻ cotransporter (NCC), Na⁺/K⁺/2Cl⁻ (NKCC1), SLC26A5, and SLC26A6) for alkalinity adaptation between the two forms of Amur ide differing in alkalinity tolerance. Specifically, close relationships among the serum Na⁺ and mRNA levels of NCC, NKCC1, and NHE, and also NKA and NBC1, in addition to serum Cl⁻ and bicarbonate transporters (e.g., SLC26A5 and SLC26A6), characterized the alkali form of Amur ide. We propose that this ecotype can ensure its transepithelial Cl⁻ and Na⁺ uptake/base secretions are highly functional, by its basolateral NKA with NBC1 and apical ionic transporters, and especially NCC incorporated with other transporters (e.g., SLC26). This suggests an evolved strong ability to maintain an ion osmotic and acid-base balance for more effectively facilitating its adaptability to the high alkaline environment. This study provides new insights into the physiological responses of the alkaline form of the Amur ide fish for adapting to extreme alkaline conditions. This information could be used as a reference to cultivating alkaline-tolerant fish species in abandoned alkaline waters.

Balanced polymorphism fuels rapid selection in an invasive crab despite high gene flow and low genetic diversity

Adaptation across environmental gradients has been demonstrated in numerous systems with extensive dispersal, despite high gene flow and consequently low genetic structure. The speed and mechanisms by which such adaptation occurs remain poorly resolved, but are critical to understanding species spread and persistence in a changing world. Here, we investigate these mechanisms in the European green crab Carcinus maenas, a globally distributed invader. We focus on a northwestern Pacific population that spread across >12 degrees of latitude in 10 years from a single source, following its introduction <35 years ago. Using six locations spanning >1500 km, we examine genetic structure using 9376 single nucleotide polymorphisms (SNPs). We find high connectivity among five locations, with significant structure between these locations and an enclosed lagoon with limited connectivity to the coast. Among the five highly connected locations, the only structure observed was a cline driven by a handful of SNPs strongly associated with latitude and winter temperature. These SNPs are almost exclusively found in a large cluster of genes in strong linkage disequilibrium that was previously identified as a candidate for cold tolerance adaptation in this species. This region may represent a balanced polymorphism that evolved to promote rapid adaptation in variable environments despite high gene flow, and which now contributes to successful invasion and spread in a novel environment. This research suggests an answer to the paradox of genetically depauperate yet successful invaders: populations may be able to adapt via a few variants of large effect despite low overall diversity.

Negative relationship between thermal tolerance and plasticity in tolerance emerges during experimental evolution in a widespread marine invertebrate

Whether populations can adapt to predicted climate change conditions, and how rapidly, are critical questions for the management of natural systems. Experimental evolution has become an important tool to answer these questions. In order to provide useful, realistic insights into the adaptive response of populations to climate change, there needs to be careful consideration of how genetic differentiation and phenotypic plasticity interact to generate observed phenotypic changes. We exposed three populations of the widespread copepod Acartia tonsa (Crustacea) to chronic, sublethal temperature selection for 15 generations. We generated thermal survivorship curves at regular intervals both during and after this period of selection to track the evolution of thermal tolerance. Using reciprocal transplants between ambient and warming conditions, we also tracked changes in the strength of phenotypic plasticity in thermal tolerance. We observed significant increases in thermal tolerance in the Warming lineages, while plasticity in thermal tolerance was strongly reduced. We suggest these changes are driven by a negative relationship between thermal tolerance and plasticity in thermal tolerance. Our results indicate that adaptation to warming through an increase in thermal tolerance might not reduce vulnerability to climate change if the increase comes at the expense of tolerance plasticity. These results illustrate the importance of considering changes in both a trait of interest and the trait plasticity during experimental evolution.

Invasion of the body snatchers: the role of parasite introduction in host distribution and response to salinity in invaded estuaries

In dynamic systems, organisms are faced with variable selective forces that may impose trade-offs. In estuaries, salinity is a strong driver of organismal diversity, while parasites shape species distributions and demography. We tested for trade-offs between low-salinity stress and parasitism in an invasive castrating parasite and its mud crab host along salinity gradients of two North Carolina rivers. We performed field surveys every six to eight weeks over 3 years to determine factors influencing parasite prevalence, host abundance, and associated taxa diversity. We also looked for signatures of low-salinity stress in the host by examining its response (time-to-right and gene expression) to salinity. We found salinity and temperature significantly affected parasite prevalence, with low-salinity sites (less than 10 practical salinity units (PSU)) lacking infection, and populations in moderate salinities at warmer temperatures reaching prevalence as high as 60%. Host abundance was negatively associated with parasite prevalence. Host gene expression was plastic to acclimation salinity, but several osmoregulatory and immune-related genes demonstrated source-dependent salinity response. We identified a genetic marker that was strongly associated with salinity against a backdrop of no neutral genetic structure, suggesting possible selection on standing variation. Our study illuminates how selective trade-offs in naturally dynamic systems may shape host evolutionary ecology.

Using seasonal genomic changes to understand historical adaptation to new environments: Parallel selection on stickleback in highly-variable estuaries

Parallel evolution is considered strong evidence for natural selection. However, few studies have investigated the process of parallel selection as it plays out in real time. The common approach is to study historical signatures of selection in populations already well adapted to different environments. Here, to document selection under natural conditions, we study six populations of threespine stickleback (Gasterosteus aculeatus) inhabiting bar-built estuaries that undergo seasonal cycles of environmental changes. Estuaries are periodically isolated from the ocean due to sandbar formation during dry summer months, with concurrent environmental shifts that resemble the long-term changes associated with postglacial colonization of freshwater habitats by marine populations. We used pooled whole-genome sequencing to track seasonal allele frequency changes in six of these populations and search for signatures of natural selection. We found consistent changes in allele frequency across estuaries, suggesting a potential role for parallel selection. Functional enrichment among candidate genes included transmembrane ion transport and calcium binding, which are important for osmoregulation and ion balance. The genomic changes that occur in threespine stickleback from bar-built estuaries could provide a glimpse into the early stages of adaptation that have occurred in many historical marine to freshwater transitions.

Molecular mechanisms of local adaptation for salt-tolerance in a treefrog

Salinization is a global phenomenon affecting ecosystems and forcing freshwater organisms to deal with increasing levels of ionic stress. However, our understanding of mechanisms that permit salt tolerance in amphibians is limited. This study investigates mechanisms of salt tolerance in locally adapted, coastal populations of a treefrog, Hyla cinerea. Using a common garden experiment, we (i) determine the extent that environment (i.e., embryonic and larval saltwater exposure) or genotype (i.e., coastal vs. inland) affects developmental benchmarks and transcriptome expression, and (ii) identify genes that may underpin differences in saltwater tolerance. Differences in gene expression, survival, and plasma osmolality were most strongly associated with genotype. Population genetic analyses on expressed genes also delineated coastal and inland groups based on genetic similarity. Coastal populations differentially expressed osmoregulatory genes including ion transporters (atp1b1, atp6V1g2, slc26a), cellular adhesion components (cdh26, cldn1, gjb3, ocln), and cytoskeletal components (odc1-a, tgm3). Several of these genes are the same genes expressed by euryhaline fish after exposure to freshwater, which is a novel finding for North American amphibians and suggests that these genes may be associated with local salinity adaptation. Coastal populations also highly expressed glycerol-3-phosphate dehydrogenase 1 (gpd1), which indicates they use glycerol as a compatible osmolyte to reduce water loss - another mechanism of saltwater tolerance previously unknown in frogs. These data signify that Hyla cinerea inhabiting coastal, brackish wetlands have evolved a salt-tolerant ecotype, and highlights novel candidate pathways that can lead to salt tolerance in freshwater organisms facing habitat salinization.

Tolerance Patterns and Transcriptomic Response to Extreme and Fluctuating Salinities across Populations of the Intertidal Copepod Tigriopus californicus

AbstractPopulations that tolerate extreme environmental conditions with frequent fluctuations can give valuable insights into physiological limits and adaptation. In some estuarine and marine ecosystems, organisms must adapt to extreme and fluctuating salinities, but not much is known about how varying salinities impact local adaptation across a wide geographic range. We used eight geographically and genetically divergent populations of the intertidal copepod Tigriopus californicus to test whether northern populations have greater tolerance to low salinity stresses, as they experience greater precipitation and less evaporation. We used a common-garden experiment approach and exposed all populations to acute low (1 and 3 ppt) and high (110 and 130 ppt) salinities for 24 h and to a fluctuation between baseline salinity and moderate low (7 ppt) and high (80 ppt) salinities for over 49 h. We also performed RNA sequencing at several time points during the fluctuation between baseline and salinity of 7 ppt to understand the molecular basis of divergence between two populations with differing physiological responses. We present these novel findings: (1) acute low salinity conditions caused more deaths than high salinity; (2) molecular processes that elevate proline levels increased in salinity of 7 ppt, which contrasts with other physiological studies in T. californicus that mainly associated accumulation of proline with hyperosmotic stress; and (3) tolerance to a salinity fluctuation did not follow a latitudinal trend but was instead governed by a complex interplay of factors, including population and duration of salinity stress. This highlights the importance of including a wider variety of environmental conditions in empirical studies to understand local adaptation.

Multi-omic approach provides insights into osmoregulation and osmoconformation of the crab Scylla paramamosain

Osmoregulation and osmoconformation are two mechanisms through which aquatic animals adapt to salinity fluctuations. The euryhaline crab Scylla paramamosain, being both an osmoconformer and osmoregulator, is an excellent model organism to investigate salinity adaptation mechanisms in brachyurans. In the present study, we used transcriptomic and proteomic approaches to investigate the response of S. paramamosain to salinity stress. Crabs were transferred from a salinity of 25 ppt to salinities of 5 ppt or 33 ppt for 6 h and 10 days. Data from both approaches revealed that exposure to 5 ppt resulted in upregulation of ion transport and energy metabolism associated genes. Notably, acclimation to low salinity was associated with early changes in gene expression for signal transduction and stress response. In contrast, exposure to 33 ppt resulted in upregulation of genes related to amino acid metabolism, and amino acid transport genes were upregulated only at the early stage of acclimation to this salinity. Our study reveals contrasting mechanisms underlying osmoregulation and osmoconformation within the salinity range of 5–33 ppt in the mud crab, and provides novel candidate genes for osmotic signal transduction, thereby providing insights on understanding the salinity adaptation mechanisms of brachyuran crabs.

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.

Spiralian Genomes Reveal Gene Family Expansions Associated with Adaptation to Freshwater

The colonization of freshwater habitats by marine-adapted organisms represents a major transition that has only occurred a few times in the evolution of animals. Only around half of the extant animal phyla have representatives in both marine and freshwater environments and even within those phyla some major clades are restricted to marine environments. Moving from marine to freshwater environments can create severe osmotic and ionic stresses and the mechanisms that animals have used to adapt to those stresses are still not well understood. In this study, we downloaded amino acid sequence data from 11 spiralian animal genomes (four freshwater taxa representing four different phyla as well as 7 marine taxa) and identified a number of gene family expansions that have occurred exclusively in the freshwater lineages. Further investigation of these gene families and the timing and nature of their expansions will illuminate one of the major evolutionary transitions in the history of life on Earth.

Geographical origin determines responses to salinity of Mediterranean caddisflies

Many freshwater ecosystems worldwide, and particularly Mediterranean ones, show increasing levels of salinity. These changes in water conditions could affect abundance and distribution of inhabiting species as well as the provision of ecosystem services. In this study we conduct laboratory experiments using the macroinvertebrate Smicridea annulicornis as a model organism. Our factorial experiments were designed to evaluate the effects of geographical origin of organisms and salinity levels on survival and behavioral responses of caddisflies. The experimental organisms were captured from rivers belonging to three hydrological basins along a 450 Km latitudinal gradient in the Mediterranean region of Chile. Animals were exposed to three conductivity levels, from 180 to 1400 μS/cm, close to the historical averages of the source rivers. We measured the behavioral responses to experimental stimuli and the survival time. Our results showed that geographical origin shaped the behavioral and survival responses to salinity. In particular, survival and activity decreased more strongly with increasing salinity in organisms coming from more dilute waters. This suggests local adaptation to be determinant for salinity responses in this benthic invertebrate species. In the current scenario of fast temporal and spatial changes in water levels and salt concentration, the conservation of geographic intra-specific variation of aquatic species is crucial for lowering the risk of salinity-driven biodiversity loss.

Epigenetic patterns associated with an ascidian invasion: a comparison of closely related clades in their native and introduced ranges

Environmentally induced epigenetic modifications have been proposed as one mechanism underlying rapid adaptive evolution of invasive species. Didemnum vexillum is an invasive colonial ascidian that has established in many coastal waters worldwide. Phylogenetic analyses have revealed that D. vexillum populations consist of two distinct clades; clade B appears to be restricted to the native range (Japan), whereas clade A is found in many regions throughout the world, including New Zealand. The spread of D. vexillum clade A suggests that it might be intrinsically more invasive than clade B, despite low levels of genetic diversity compared to populations from the native region. This study investigated whether D. vexillum clade A exhibits epigenetic signatures (specifically differences in DNA methylation) associated with invasiveness. Global DNA methylation patterns were significantly different between introduced clade A colonies, and both clades A and B in the native range. Introduced colonies also showed a significant reduction in DNA methylation levels, which could be a mechanism for increasing phenotypic plasticity. High levels of DNA methylation diversity were maintained in the introduced population, despite reduced levels of genetic diversity, which may allow invasive populations to respond quickly to changes in new environments. Epigenetic changes induced during the invasion process could provide a means for rapid adaptation despite low levels of genetic variation in introduced populations.

The Toxicogenome of Hyalella azteca: A Model for Sediment Ecotoxicology and Evolutionary Toxicology

Hyalella azteca is a cryptic species complex of epibenthic amphipods of interest to ecotoxicology and evolutionary biology. It is the primary crustacean used in North America for sediment toxicity testing and an emerging model for molecular ecotoxicology. To provide molecular resources for sediment quality assessments and evolutionary studies, we sequenced, assembled, and annotated the genome of the H. azteca U.S. Lab Strain. The genome quality and completeness is comparable with other ecotoxicological model species. Through targeted investigation and use of gene expression data sets of H. azteca exposed to pesticides, metals, and other emerging contaminants, we annotated and characterized the major gene families involved in sequestration, detoxification, oxidative stress, and toxicant response. Our results revealed gene loss related to light sensing, but a large expansion in chemoreceptors, likely underlying sensory shifts necessary in their low light habitats. Gene family expansions were also noted for cytochrome P450 genes, cuticle proteins, ion transporters, and include recent gene duplications in the metal sequestration protein, metallothionein. Mapping of differentially expressed transcripts to the genome significantly increased the ability to functionally annotate toxicant responsive genes. The H. azteca genome will greatly facilitate development of genomic tools for environmental assessments and promote an understanding of how evolution shapes toxicological pathways with implications for environmental and human health.

The Toxicogenome of Hyalella azteca: A Model for Sediment Ecotoxicology and Evolutionary Toxicology

Hyalella azteca is a cryptic species complex of epibenthic amphipods of interest to ecotoxicology and evolutionary biology. It is the primary crustacean used in North America for sediment toxicity testing and an emerging model for molecular ecotoxicology. To provide molecular resources for sediment quality assessments and evolutionary studies, we sequenced, assembled, and annotated the genome of the H. azteca U.S. Lab Strain. The genome quality and completeness is comparable with other ecotoxicological model species. Through targeted investigation and use of gene expression data sets of H. azteca exposed to pesticides, metals, and other emerging contaminants, we annotated and characterized the major gene families involved in sequestration, detoxification, oxidative stress, and toxicant response. Our results revealed gene loss related to light sensing, but a large expansion in chemoreceptors, likely underlying sensory shifts necessary in their low light habitats. Gene family expansions were also noted for cytochrome P450 genes, cuticle proteins, ion transporters, and include recent gene duplications in the metal sequestration protein, metallothionein. Mapping of differentially expressed transcripts to the genome significantly increased the ability to functionally annotate toxicant responsive genes. The H. azteca genome will greatly facilitate development of genomic tools for environmental assessments and promote an understanding of how evolution shapes toxicological pathways with implications for environmental and human health.

Predicting evolutionary rescue via evolving plasticity in stochastic environments

Phenotypic plasticity and its evolution may help evolutionary rescue in a novel and stressful environment, especially if environmental novelty reveals cryptic genetic variation that enables the evolution of increased plasticity. However, the environmental stochasticity ubiquitous in natural systems may alter these predictions, because high plasticity may amplify phenotype-environment mismatches. Although previous studies have highlighted this potential detrimental effect of plasticity in stochastic environments, they have not investigated how it affects extinction risk in the context of evolutionary rescue and with evolving plasticity. We investigate this question here by integrating stochastic demography with quantitative genetic theory in a model with simultaneous change in the mean and predictability (temporal autocorrelation) of the environment. We develop an approximate prediction of long-term persistence under the new pattern of environmental fluctuations, and compare it with numerical simulations for short- and long-term extinction risk. We find that reduced predictability increases extinction risk and reduces persistence because it increases stochastic load during rescue. This understanding of how stochastic demography, phenotypic plasticity, and evolution interact when evolution acts on cryptic genetic variation revealed in a novel environment can inform expectations for invasions, extinctions, or the emergence of chemical resistance in pests.

The chicken or the egg? Adaptation to desiccation and salinity tolerance in a lineage of water beetles

Transitions from fresh to saline habitats are restricted to a handful of insect lineages, as the colonization of saline waters requires specialized mechanisms to deal with osmotic stress. Previous studies have suggested that tolerance to salinity and desiccation could be mechanistically and evolutionarily linked, but the temporal sequence of these adaptations is not well established for individual lineages. We combined molecular, physiological and ecological data to explore the evolution of desiccation resistance, hyporegulation ability (i.e., the ability to osmoregulate in hyperosmotic media) and habitat transitions in the water beetle genus Enochrus subgenus Lumetus (Hydrophilidae). We tested whether enhanced desiccation resistance evolved before increases in hyporegulation ability or vice versa, or whether the two mechanisms evolved in parallel. The most recent ancestor of Lumetus was inferred to have high desiccation resistance and moderate hyporegulation ability. There were repeated shifts between habitats with differing levels of salinity in the radiation of the group, those to the most saline habitats generally occurring more rapidly than those to less saline ones. Significant and accelerated changes in hyporegulation ability evolved in parallel with smaller and more progressive increases in desiccation resistance across the phylogeny, associated with the colonization of meso- and hypersaline waters during global aridification events. All species with high hyporegulation ability were also desiccation-resistant, but not vice versa. Overall, results are consistent with the hypothesis that desiccation resistance mechanisms evolved first and provided the physiological basis for the development of hyporegulation ability, allowing these insects to colonize and diversify across meso- and hypersaline habitats.

RNA sequencing analysis of transcriptional change in the freshwater mussel Elliptio complanata after environmentally relevant sodium chloride exposure

To identify potential biomarkers of salt stress in a freshwater sentinel species, we examined transcriptional responses of the common mussel Elliptio complanata to controlled sodium chloride (NaCl) exposures. Ribonucleic acid sequencing (RNA-Seq) of mantle tissue identified 481 transcripts differentially expressed in adult mussels exposed to 2 ppt NaCl (1.2 ppt chloride) for 7 d, of which 290 had nonoverlapping intervals. Differentially expressed gene categories included ion and transmembrane transport, oxidoreductase activity, maintenance of protein folding, and amino acid metabolism. The rate-limiting enzyme for synthesis of taurine, an amino acid frequently linked to osmotic stress in aquatic species, was upregulated, as was the transmembrane ion pump sodium/potassium adenosine 5'-triphosphatase. These patterns confirm a primary transcriptional response to the experimental dose, albeit likely overlapping with nonspecific secondary stress responses. Substantial involvement of the heat shock protein 70 chaperone family and the water-transporting aquaporin family was not detected, however, in contrast to some studies in other bivalves. A subset of the most significantly regulated genes was confirmed by quantitative polymerase chain reaction in an independent sample. Cluster analysis showed separation of mussels exposed to 2 ppt NaCl from control mussels in multivariate space, but mussels exposed to 1 ppt NaCl were largely indistinguishable from controls. Transcriptome-scale analysis of salt exposure under laboratory conditions efficiently identified candidate biomarkers for further functional analysis and field validation. Environ Toxicol Chem 2017;36:2352-2366. © Published 2017 Wiley Periodicals Inc. on behalf of SETAC. This article is a US government work and, as such, is in the public domain in the United States of America.

Comparative molecular analyses of select pH- and osmoregulatory genes in three freshwater crayfish Cherax quadricarinatus, C. destructor and C. cainii

Systemic acid-base balance and osmotic/ionic regulation in decapod crustaceans are in part maintained by a set of transport-related enzymes such as carbonic anhydrase (CA), Na⁺/K⁺-ATPase (NKA), H⁺-ATPase (HAT), Na⁺/K⁺/2Cl⁻ cotransporter (NKCC), Na⁺/Cl⁻/HCO\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{upgreek} \usepackage{mathrsfs} \setlength{\oddsidemargin}{-69pt} \begin{document} }{}${}_{3}^{-}$\end{document}3− cotransporter (NBC), Na⁺/H⁺ exchanger (NHE), Arginine kinase (AK), Sarcoplasmic Ca⁺²-ATPase (SERCA) and Calreticulin (CRT). We carried out a comparative molecular analysis of these genes in three commercially important yet eco-physiologically distinct freshwater crayfish, Cherax quadricarinatus, C. destructor and C. cainii, with the aim to identify mutations in these genes and determine if observed patterns of mutations were consistent with the action of natural selection. We also conducted a tissue-specific expression analysis of these genes across seven different organs, including gills, hepatopancreas, heart, kidney, liver, nerve and testes using NGS transcriptome data. The molecular analysis of the candidate genes revealed a high level of sequence conservation across the three Cherax sp. Hyphy analysis revealed that all candidate genes showed patterns of molecular variation consistent with neutral evolution. The tissue-specific expression analysis showed that 46% of candidate genes were expressed in all tissue types examined, while approximately 10% of candidate genes were only expressed in a single tissue type. The largest number of genes was observed in nerve (84%) and gills (78%) and the lowest in testes (66%). The tissue-specific expression analysis also revealed that most of the master genes regulating pH and osmoregulation (CA, NKA, HAT, NKCC, NBC, NHE) were expressed in all tissue types indicating an important physiological role for these genes outside of osmoregulation in other tissue types. The high level of sequence conservation observed in the candidate genes may be explained by the important role of these genes as well as potentially having a number of other basic physiological functions in different tissue types.

Evolutionary responses to crude oil from the Deepwater Horizon oil spill by the copepod Eurytemora affinis

The BP Deepwater Horizon Oil Disaster was the most catastrophic offshore oil spill in U.S. history, yet we still have a poor understanding of how organisms could evolve in response to the toxic effects of crude oil. This study offers a rare analysis of how fitness-related traits could evolve rapidly in response to crude oil toxicity. We examined evolutionary responses of populations of the common copepod Eurytemora affinis residing in the Gulf of Mexico, by comparing crude oil tolerance of populations collected before versus after the Deepwater Horizon oil spill of 2010. In addition, we imposed laboratory selection for crude oil tolerance for ~8 generations, using an E. affinis population collected from before the oil spill. We found evolutionary increases in crude oil tolerance in the wild population following the oil spill, relative to the population collected before the oil spill. The post-oil spill population showed increased survival and rapid development time in the presence of crude oil. In contrast, evolutionary responses following laboratory selection were less clear; though, development time from metamorphosis to adult in the presence of crude oil did become more rapid after selection. We did find that the wild population, used in both experiments, harbored significant genetic variation in crude oil tolerance, upon which selection could act. Thus, our study indicated that crude oil tolerance could evolve, but perhaps not on the relatively short time scale of the laboratory selection experiment. This study contributes novel insights into evolutionary responses to crude oil, in directly examining fitness-related traits before and after an oil spill, and in observing evolutionary responses following laboratory selection.

Regulation of gene expression is associated with tolerance of the Arctic copepod Calanus glacialis to CO2-acidified sea water

Ocean acidification is the increase in seawater pCO ₂ due to the uptake of atmospheric anthropogenic CO ₂, with the largest changes predicted to occur in the Arctic seas. For some marine organisms, this change in pCO ₂, and associated decrease in pH, represents a climate change-related stressor. In this study, we investigated the gene expression patterns of nauplii of the Arctic copepod Calanus glacialis cultured at low pH levels. We have previously shown that organismal-level performance (development, growth, respiration) of C. glacialis nauplii is unaffected by low pH. Here, we investigated the molecular-level response to lowered pH in order to elucidate the physiological processes involved in this tolerance. Nauplii from wild-caught C. glacialis were cultured at four pH levels (8.05, 7.9, 7.7, 7.5). At stage N6, mRNA was extracted and sequenced using RNA-seq. The physiological functionality of the proteins identified was categorized using Gene Ontology and KEGG pathways. We found that the expression of 151 contigs varied significantly with pH on a continuous scale (93% downregulated with decreasing pH). Gene set enrichment analysis revealed that, of the processes downregulated, many were components of the universal cellular stress response, including DNA repair, redox regulation, protein folding, and proteolysis. Sodium:proton antiporters were among the processes significantly upregulated, indicating that these ion pumps were involved in maintaining cellular pH homeostasis. C. glacialis significantly alters its gene expression at low pH, although they maintain normal larval development. Understanding what confers tolerance to some species will support our ability to predict the effects of future ocean acidification on marine organisms.

Selection by higher-order effects of salinity and bacteria on early life-stages of Western Baltic spring-spawning herring

Habitat stratification by abiotic and biotic factors initiates divergence of populations and leads to ecological speciation. In contrast to fully marine waters, the Baltic Sea is stratified by a salinity gradient that strongly affects fish physiology, distribution, diversity and virulence of important marine pathogens. Animals thus face the challenge to simultaneously adapt to the concurrent salinity and cope with the selection imposed by the changing pathogenic virulence. Western Baltic spring‐spawning herring (Clupea harengus) migrate to spawning grounds characterized by different salinities to which herring are supposedly adapted. We hypothesized that herring populations do not only have to cope with different salinity levels but that they are simultaneously exposed to higher‐order effects that accompany the shifts in salinity, that is induced pathogenicity of Vibrio bacteria in lower saline waters. To experimentally evaluate this, adults of two populations were caught in their spawning grounds and fully reciprocally crossed within and between populations. Larvae were reared at three salinity levels, representing the spawning ground salinity of each of the two populations, or Atlantic salinity conditions resembling the phylogenetic origin of Clupea harengus. In addition, larvae were exposed to a Vibrio spp. infection. Life‐history traits and gene expression analysis served as response variables. Herring seem adapted to Baltic Sea conditions and cope better with low saline waters. However, upon a bacterial infection, herring larvae suffer more when kept at lower salinities implying reduced resistance against Vibrio or higher Vibrio virulence. In the context of recent climate change with less saline marine waters in the Baltic Sea, such interactions may constitute key future stressors.

Toxicological perspective on the osmoregulation and ionoregulation physiology of major ions by freshwater animals: Teleost fish, crustacea, aquatic insects, and Mollusca

Anthropogenic sources increase freshwater salinity and produce differences in constituent ions compared with natural waters. Moreover, ions differ in physiological roles and concentrations in intracellular and extracellular fluids. Four freshwater taxa groups are compared, to investigate similarities and differences in ion transport processes and what ion transport mechanisms suggest about the toxicity of these or other ions in freshwater. Although differences exist, many ion transporters are functionally similar and may belong to evolutionarily conserved protein families. For example, the Na+ /H+ -exchanger in teleost fish differs from the H+ /2Na+ (or Ca2+ )-exchanger in crustaceans. In osmoregulation, Na+ and Cl- predominate. Stenohaline freshwater animals hyperregulate until they are no longer able to maintain hypertonic extracellular Na+ and Cl- concentrations with increasing salinity and become isotonic. Toxic effects of K+ are related to ionoregulation and volume regulation. The ionic balance between intracellular and extracellular fluids is maintained by Na+ /K+ -adenosine triphosphatase (ATPase), but details are lacking on apical K+ transporters. Elevated H+ affects the maintenance of internal Na+ by Na+ /H+ exchange; elevated HCO3- inhibits Cl- uptake. The uptake of Mg2+ occurs by the gills or intestine, but details are lacking on Mg2+ transporters. In unionid gills, SO42- is actively transported, but most epithelia are generally impermeant to SO42- . Transporters of Ca2+ maintain homeostasis of dissolved Ca2+ . More integration of physiology with toxicology is needed to fully understand freshwater ion effects. Environ Toxicol Chem 2017;36:576-600. Published 2016 Wiley Periodicals Inc. on behalf of SETAC. This article is a US government work and, as such, is in the public domain in the United States of America.