Shared and unique patterns of phenotypic diversification along a stream gradient in two congeneric species

Rapid development of increased neonicotinoid tolerance in non-target freshwater amphipods

The comprehensive assessment of the long-term impacts of constant exposure to pollutants on wildlife populations remains a relatively unexplored area of ecological risk assessment. Empirical evidence to suggest that multigenerational exposure affects the susceptibility of organisms is scarce, and the underlying mechanisms in the natural environment have yet to be fully understood. In this study, we first examined the arthropod candidate species, Gammarus roeselii that - unlike closely related species - commonly occurs in many contaminated river systems of Central Europe. This makes it a suitable study organism to investigate the development of tolerances and phenotypic adaptations along pollution gradients. In a 96-h acute toxicity assay with the neonicotinoid thiacloprid, we indeed observed a successive increase in tolerance in populations coming from contaminated regions. This was accompanied by a certain phenotypic change, with increased investment into reproduction. To address the question of whether these changes are plastic or emerged from longer lasting evolutionary processes, we conducted a multigeneration experiment in the second part of our study. Here, we used closely-related Hyalella azteca and pre-exposed them for multiple generations to sublethal concentrations of thiacloprid in a semi-static design (one week renewal of media containing 0.1 or 1.0 µg/L thiacloprid). The pre-exposed individuals were then used in acute toxicity assays to see how quickly such adaptive responses can develop. Over only two generations, the tolerance to the neonicotinoid almost doubled, suggesting developmental plasticity as a plausible mechanism for the rapid adaptive response to strong selection factors such as neonicotinoid insecticides. It remains to be discovered whether the plasticity of rapidly developed tolerance is species-specific and explains why closely related species - which may not have comparable adaptive response capabilities - disappear in polluted habitats. Overall, our findings highlight the neglected role of developmental plasticity during short- and long-term exposure of natural populations to pollution. Moreover, our results show that even pollutant levels seven times lower than concentrations found in the study region have a clear impact on the developmental trajectories of non-target species.

Cryptic Species in Ecotoxicology

The advent of genetic methods has led to the discovery of an increasing number of species that previously could not be distinguished from each other on the basis of morphological characteristics. Even though there has been an exponential growth of publications on cryptic species, such species are rarely considered in ecotoxicology. Thus, the particular question of ecological differentiation and the sensitivity of closely related cryptic species is rarely addressed. Tackling this question, however, is of key importance for evolutionary ecology, conservation biology, and, in particular, regulatory ecotoxicology. At the same time, the use of species with (known or unknown) cryptic diversity might be a reason for the lack of reproducibility of ecotoxicological experiments and implies a false extrapolation of the findings. Our critical review includes a database and literature search through which we investigated how many of the species most frequently used in ecotoxicological assessments show evidence of cryptic diversity. We found a high proportion of reports indicating overlooked species diversity, especially in invertebrates. In terrestrial and aquatic realms, at least 67% and 54% of commonly used species, respectively, were identified as cryptic species complexes. The issue is less prominent in vertebrates, in which we found evidence for cryptic species complexes in 27% of aquatic and 6.7% of terrestrial vertebrates. We further exemplified why different evolutionary histories may significantly determine cryptic species' ecology and sensitivity to pollutants. This in turn may have a major impact on the results of ecotoxicological tests and, consequently, the outcome of environmental risk assessments. Finally, we provide a brief guideline on how to deal practically with cryptic diversity in ecotoxicological studies in general and its implementation in risk assessment procedures in particular. Environ Toxicol Chem 2023;42:1889-1914. © 2023 The Authors. Environmental Toxicology and Chemistry published by Wiley Periodicals LLC on behalf of SETAC.

Infection with acanthocephalans increases tolerance of Gammarus roeselii (Crustacea: Amphipoda) to pyrethroid insecticide deltamethrin

Crustacean amphipods serve as intermediate hosts for parasites and are at the same time sensitive indicators of environmental pollution in aquatic ecosystems. The extent to which interaction with the parasite influences their persistence in polluted ecosystems is poorly understood. Here, we compared infections of Gammarus roeselii with two species of Acanthocephala, Pomphorhynchus laevis, and Polymorphus minutus, along a pollution gradient in the Rhine-Main metropolitan region of Frankfurt am Main, Germany. Prevalence of P. laevis was very low at the unpolluted upstream reaches (P ≤ 3%), while higher prevalence (P ≤ 73%) and intensities of up to 9 individuals were found further downstream-close to an effluent of a large wastewater treatment plant (WWTP). Co-infections of P. minutus and P. laevis occurred in 11 individuals. Highest prevalence of P. minutus was P ≤ 9% and one parasite per amphipod host was the maximum intensity recorded. In order to assess whether the infection affects survival in the polluted habitats, we tested the sensitivity of infected and uninfected amphipods towards the pyrethroide insecticide deltamethrin. We found an infection-dependent difference in sensitivity within the first 72 h, with an effect concentration (24 h EC₅₀) of 49.8 ng/l and 26.6 ng/l for infected and uninfected G. roeselii, respectively. Whereas final host abundance might partially explain the high prevalence of P. laevis in G. roeselii, the results of the acute toxicity test suggest a beneficial effect of acanthocephalan infection for G. roeselii at polluted sites. A strong accumulation of pollutants in the parasite could serve as a sink for pesticide exposure of the host. Due to the lack of a co-evolutionary history between parasite and host and a lack of behavioral manipulation (unlike in co-evolved gammarids), the predation risk by fish remains the same, explaining high local prevalence. Thus, our study exemplifies how organismic interaction can favor the persistence of a species under chemical pollution.

Non-parallel morphological divergence following colonization of a new host plant

Adaptation to new ecological niches is known to spur population diversification and may lead to speciation if gene flow is ceased. While adaptation to the same ecological niche is expected to be parallel, it is more difficult to predict whether selection against maladaptive hybridization in secondary sympatry results in parallel divergence also in traits that are not directly related to the ecological niches. Such parallelisms in response to selection for reproductive isolation can be identified through estimating parallelism in reproductive character displacement across different zones of secondary contact. Here, we use a host shift in the phytophagous peacock fly Tephritis conura, with both host races represented in two geographically separate areas East and West of the Baltic Sea to investigate convergence in morphological adaptations. We asked (i) if there are consistent morphological adaptations to a host plant shift and (ii) if the response to secondary sympatry with the alternate host race is parallel across contact zones. We found surprisingly low and variable, albeit significant, divergence between host races. Only one trait, the length of the female ovipositor, which serves an important function in the interaction with the hosts, was consistently different between host races. Instead, co-existence with the other host race significantly affected the degree of morphological divergence, but the divergence was largely driven by different traits in different contact zones. Thus, local stochastic fixation or reinforcement could generate trait divergence, and additional evidence is needed to conclude whether divergence is locally adaptive.

Predator Environment Does Not Predict Life History in the Morphologically Constrained Fish Alfaro cultratus (Cyprinodontiformes: Poeciliidae)

Predation is known to have a significant effect on life history diversification in a variety of species. However, physical constraints of body shape and size can sometimes limit life history divergence. We test this idea in the Costa Rican livebearing fishAlfaro cultratus. Individuals in this species have a narrow body and keeled ventral surface, and females do not develop a distended abdomen when pregnant like other livebearing fishes. Here, we describe the life history ofA. cultratusfrom 20 different populations across both high-predation and low-predation environments. We found significantly lower reproductive allotment in females from high-predation environments than in females from low-predation environments, but no significant difference in female or male size at maturity, number of offspring produced by females, or size of offspring. We found thatA. cultratusexhibit isometric patterns of allocation for clutch dry mass in relation to female dry mass in high-predation and low-predation environments. Our results suggest that body shape constraints in this species limit the life history divergence we typically see between populations from high-predation and low-predation environments in other species.

A century later: Adaptive plasticity and rapid evolution contribute to geographic variation in invasive mosquitofish

One century after their introduction to Europe, eastern mosquitofish (Gambusia holbrooki) represent a natural experiment to determine the relative contributions of adaptive plasticity and rapid evolutionary change in creating large-scale geographic variation in phenotypes. We evaluated the population-genetic structure and invasion history based on allele length polymorphisms of 15 nuclear microsatellites, which we quantified for N = 660 individuals from 23 populations sampled in 2013 across the invasive range of G. holbrooki in Europe. We analysed body-shape and life-history variation in N = 1331 individuals from 36 populations, sampled in 2013 and 2017, and tested heritability of phenotypic differences in a subset of four populations using a common-garden experiment. The genetic structure of wild-caught individuals suggested a single introduction for all European mosquitofish, which were genetically impoverished compared to their native counterparts. We found some convergent patterns of phenotypic divergence across native and invasive climatic gradients (e.g., increased body size in colder/more northern populations); however, several phenotypic responses were not consistent between sampling years, pointing towards plastic phenotypes. Our analysis of common-garden reared individuals uncovered moderate heritability estimates only for two measures of male body size (intraclass correlation coefficient, ICC = 0.628 and 0.556) and offspring fat content (ICC = 0.734), while suggesting high levels of plasticity in most other phenotypic traits (ICC ≤ 0.407). Our results highlight the importance of phenotypic plasticity in invasive species during range expansions and demonstrate that strong selective pressures-in this case towards increased body size in colder environments-simultaneously promote rapid evolutionary divergence.

Using native and invasive livebearing fishes (Poeciliidae, Teleostei) for the integrated biological assessment of pollution in urban streams

Invasive species are increasingly replacing native species, especially in anthropogenically transformed or polluted habitats. This opens the possibility to use invasive species as indicator taxa for the biological assessment of pollution. Integrated biological assessment, however, additionally relies on the application of multiple approaches to quantify physiological or cytogenetic responses to pollution within the same focal species. This is challenging when species are restricted to either polluted or unpolluted sites. Here, we make use of a small group of neotropical livebearing fishes (family Poeciliidae) for the integrated biological assessment of water quality. Comparing urban and suburban stream sections that receive varying degrees of pollution from industrial and domestic waste waters in and around the Brazilian city of Uberlândia, we demonstrate that two members of this family may indeed serve as indicators of water pollution levels. The native species Phalloceros caudimaculatus appears to be replaced by invasive guppies (Poecilia reticulata) at heavily polluted sites. Nevertheless, we demonstrate that both species could be used for the assessment of bioaccumulation of heavy metals (Pb, Cu, and Cr). Ambient (sediment) concentrations predicted concentrations in somatic tissue across species (R²-values between 0.74 and 0.96). Moreover, we used cytogenetic methods to provide an estimate of genotoxic effects of water pollution and found pollution levels (multiple variables, condensed into principal components) to predict the occurrence of nuclear abnormalities (e.g., frequencies of micro-nucleated cells) across species (R² between 0.69 and 0.83). The occurrence of poeciliid fishes in urban and polluted environments renders this family a prime group of focal organisms for biological water quality monitoring and assessment. Both species could be used interchangeably to assess genotoxic effects of water pollution, which may facilitate future comparative analyses over extensive geographic scales, as members of the family Poeciliidae have become invasive in tropical and subtropical regions worldwide.

Small-scale phenotypic differentiation along complex stream gradients in a non-native amphipod

Background

Selective landscapes in rivers are made up by an array of selective forces that vary from source to downstream regions or between seasons, and local/temporal variation in fitness maxima can result in gradual spatio-temporal variation of phenotypic traits. This study aimed at establishing freshwater amphipods as future model organisms to study adaptive phenotypic diversification (evolutionary divergence and/or adaptive plasticity) along stream gradients.

Methods

We collected Gammarus roeselii from 16 sampling sites in the Rhine catchment during two consecutive seasons (summer and winter). Altogether, we dissected n = 1648 individuals and quantified key parameters related to morphological and life-history diversification, including naturally selected (e.g., gill surface areas) as well as primarily sexually selected traits (e.g., male antennae). Acknowledging the complexity of selective regimes in streams and the interrelated nature of selection factors, we assessed several abiotic (e.g., temperature, flow velocity) and biotic ecological parameters (e.g., conspecific densities, sex ratios) and condensed them into four principal components (PCs).

Results

Generalized least squares models revealed pronounced phenotypic differentiation in most of the traits investigated herein, and components of the stream gradient (PCs) explained parts of the observed differences. Depending on the trait under investigation, phenotypic differentiation could be ascribed to variation in abiotic conditions, anthropogenic disturbance (influx of thermally polluted water), or population parameters. For example, female fecundity showed altitudinal variation and decreased with increasing conspecific densities, while sexual dimorphism in the length of male antennae—used for mate finding and assessment—increased with increasing population densities and towards female-biased sex ratios.

Conclusions

We provide a comprehensive protocol for comparative analyses of intraspecific variation in life history traits in amphipods. Whether the observed phenotypic differentiation over small geographical distances reflects evolutionary divergence or plasticity (or both) remains to be investigated in future studies. Independent of the mechanisms involved, variation in several traits is likely to have consequences for ecosystem functions. For example, leaf-shredding in G. roeselii strongly depends on body size, which varied in dependence of several ecological parameters.

Electronic supplementary material

The online version of this article (10.1186/s12983-019-0327-8) contains supplementary material, which is available to authorized users.

Natural and sexual selection drive multivariate phenotypic divergence along climatic gradients in an invasive fish

Invasive species that rapidly spread throughout novel distribution ranges are prime models to investigate climate-driven phenotypic diversification on a contemporary scale. Previous studies on adaptive diversification along latitudinal gradients in fish have mainly considered body size and reported either increased or decreased body size towards higher latitudes (i.e. Bergmann's rule). Our study is the first to investigate phenotypic divergence in multiple traits, including sexually selected traits (size and shape of the male copulatory organ, the gonopodium) of invasive Gambusia affinis in China. We studied body size, life history traits and morphological variation across populations spanning 17 degrees of latitude and 16 degrees of longitude. Even though we found phenotypic variation along climatic gradients to be strongest in naturally selected traits, some sexually selected traits also showed systematic gradual divergence. For example, males from southern populations possessed wider gonopodia with increased armament. Generally, males and females diverged in response to different components of climatic gradients (latitudinal or longitudinal variation) and in different trait suites. We discuss that not only temperature regimes, but also indirect effects of increased resource and mate competition (as a function of different extrinsic overwinter mortality rates) alter the selective landscape along climatic gradients.

Gradients in predation risk in a tropical river system

The importance of predation risk as a key driver of evolutionary change is exemplified by the Northern Range in Trinidad, where research on guppies living in multiple parallel streams has provided invaluable insights into the process of evolution by natural selection. Although Trinidadian guppies are now a textbook example of evolution in action, studies have generally categorized predation as a dichotomous variable, representing high or low risk. Yet, ecologists appreciate that community structure and the attendant predation risk vary substantially over space and time. Here, we use data from a longitudinal study of fish assemblages at 16 different sites in the Northern Range to quantify temporal and spatial variation in predation risk. Specifically we ask: 1) Is there evidence for a gradient in predation risk? 2) Does the ranking of sites (by risk) change with the definition of the predator community (in terms of species composition and abundance currency), and 3) Are site rankings consistent over time? We find compelling evidence that sites lie along a continuum of risk. However, site rankings along this gradient depend on how predation is quantified in terms of the species considered to be predators and the abundance currency is used. Nonetheless, for a given categorization and currency, rankings are relatively consistent over time. Our study suggests that consideration of predation gradients will lead to a more nuanced understanding of the role of predation risk in behavioral and evolutionary ecology. It also emphasizes the need to justify and report the definition of predation risk being used.

Molecular evolution and expression of oxygen transport genes in livebearing fishes (Poeciliidae) from hydrogen sulfide rich springs

Hydrogen sulfide (H2S) is a natural toxicant in some aquatic environments that has diverse molecular targets. It binds to oxygen transport proteins, rendering them non-functional by reducing oxygen-binding affinity. Hence, organisms permanently inhabiting H2S-rich environments are predicted to exhibit adaptive modifications to compensate for the reduced capacity to transport oxygen. We investigated 10 lineages of fish of the family Poeciliidae that have colonized freshwater springs rich in H2S—along with related lineages from non-sulfidic environments—to test hypotheses about the expression and evolution of oxygen transport genes in a phylogenetic context. We predicted shifts in the expression of and signatures of positive selection on oxygen transport genes upon colonization of H2S-rich habitats. Our analyses indicated significant shifts in gene expression for multiple hemoglobin genes in lineages that have colonized H2S-rich environments, and three hemoglobin genes exhibited relaxed selection in sulfidic compared to non-sulfidic lineages. However, neither changes in gene expression nor signatures of selection were consistent among all lineages in H2S-rich environments. Oxygen transport genes may consequently be predictable targets of selection during adaptation to sulfidic environments, but changes in gene expression and molecular evolution of oxygen transport genes in H2S-rich environments are not necessarily repeatable across replicated lineages.

Predation drives morphological convergence in the Gambusia panuco species group among lotic and lentic habitats

Fish morphology is often constrained by a trade-off between optimizing steady vs. unsteady swimming performance due to opposing effects of caudal peduncle size. Lotic environments tend to select for steady swimming performance, leading to smaller caudal peduncles, whereas predators tend to select for unsteady swimming performance, leading to larger caudal peduncles. However, it is unclear which aspect of performance should be optimized across heterogeneous flow and predation environments and how this heterogeneity may affect parallel phenotypic evolution. We investigated this question among four Gambusia species in north-eastern Mexico, specifically the riverine G. panuco, the spring endemics G. alvarezi and G. hurtadoi, and a fourth species, G. marshi, found in a variety of habitats with varying predation pressure in the Cuatro Ciénegas Basin and Río Salado de Nadadores. We employed a geometric morphometric analysis to examine how body shapes of both male and female fish differ among species and habitats and with piscivore presence. We found that high-predation and low-predation species diverged morphologically, with G. marshi exhibiting a variable, intermediate body shape. Within G. marshi, body morphology converged in high-predation environments regardless of flow velocity, and fish from high-predation sites had larger relative caudal peduncle areas. However, we found that G. marshi from low-predation environments diverged in morphology between sub-basins of Cuatro Ciénegas, indicating other differences among these basins that merit further study. Our results suggest that a morphological trade-off promotes parallel evolution of body shape in fishes colonizing high-predation environments and that changing predation pressure can strongly impact morphological evolution in these species.