Common mechanism underlies repeated evolution of extreme pollution tolerance

Independently evolved pollution resistance in four killifish populations is largely explained by few variants of large effect

The genetic architecture of phenotypic traits can affect the mode and tempo of trait evolution. Human-altered environments can impose strong natural selection, where successful evolutionary adaptation requires swift and large phenotypic shifts. In these scenarios, theory predicts that adaptation is due to a few adaptive variants of large effect, but empirical studies that have revealed the genetic architecture of rapidly evolved phenotypes are rare, especially for populations inhabiting polluted environments. Fundulus killifish have repeatedly evolved adaptive resistance to extreme pollution in urban estuaries. Prior studies, including genome scans for signatures of natural selection, have revealed some of the genes and pathways important for evolved pollution resistance, and provide context for the genotype-phenotype association studies reported here. We created multiple quantitative trait locus (QTL) mapping families using progenitors from four different resistant populations, and using RAD-seq genetically mapped variation in sensitivity (developmental perturbations) following embryonic exposure to a model toxicant PCB-126. We found that one to two large-effect QTL loci accounted for resistance to PCB-mediated developmental toxicity. QTLs harbored candidate genes that govern the regulation of aryl hydrocarbon receptor (AHR) signaling. One QTL locus was shared across all populations and another was shared across three populations. One QTL locus showed strong signatures of recent natural selection in the corresponding wild population but another QTL locus did not. Some candidate genes for PCB resistance inferred from genome scans in wild populations were identified as QTL, but some key candidate genes were not. We conclude that rapidly evolved resistance to the developmental defects normally caused by PCB-126 is governed by few genes of large effect. However, other aspects of resistance beyond developmental phenotypes may be governed by additional loci, such that comprehensive resistance to PCB-126, and to the mixtures of chemicals that distinguish urban estuaries more broadly, may be more genetically complex.

Bioenergetic Effects of Polycyclic Aromatic Hydrocarbon Resistance Manifest Later in Life in Offspring of Fundulus heteroclitus from the Elizabeth River

Shifts in key physiological processes can confer resistance to chemical pollutants. However, these adaptations may come with certain trade-offs, such as altered energy metabolic processes, as evident in Atlantic killifish (Fundulus heteroclitus) in Virginia's Elizabeth River (ER) that have evolved resistance to polycyclic aromatic hydrocarbons (PAHs). We seek to understand the bioenergetic costs of PAH resistance among subpopulations of Atlantic killifish with differing contamination levels in order to examine how these changes manifest across multiple life stages and how these costs might be exacerbated by additional stressors. Bioenergetics data revealed differences in metabolic rates between offspring of PAH-resistant fish and reference fish were absent or minimal in both the embryo and larval stages but pronounced at the juvenile life stage, suggesting that bioenergetic changes in pollution-adapted killifish manifest later in life. We also provide evidence that killifish from remediated sites are more sensitive to PAH exposure than killifish from nonremediated sites, suggesting loss of PAH tolerance following relaxed selection. Collectively, our data suggest that the fitness consequences associated with evolved resistance to anthropogenic stressors may manifest differently over time and depend on the magnitude of the selection pressure. This information can be valuable in effective risk and remediation assessments as well as in broadening our understanding of species responses to environmental change.

Connecting Suborganismal Data to Bioenergetic Processes: Killifish Embryos Exposed to a Dioxin-Like Compound

A core challenge for ecological risk assessment is to integrate molecular responses into a chain of causality to organismal or population-level outcomes. Bioenergetic theory may be a useful approach for integrating suborganismal responses to predict organismal responses that influence population dynamics. We describe a novel application of dynamic energy budget (DEB) theory in the context of a toxicity framework (adverse outcome pathways [AOPs]) to make quantitative predictions of chemical exposures to individuals, starting from suborganismal data. We use early-life stage exposure of Fundulus heteroclitus to dioxin-like chemicals (DLCs) and connect AOP key events to DEB processes through "damage" that is produced at a rate proportional to the internal toxicant concentration. We use transcriptomic data of fish embryos exposed to DLCs to translate molecular indicators of damage into changes in DEB parameters (damage increases somatic maintenance costs) and DEB models to predict sublethal and lethal effects on young fish. By changing a small subset of model parameters, we predict the evolved tolerance to DLCs in some wild F. heteroclitus populations, a data set not used in model parameterization. The differences in model parameters point to reduced sensitivity and altered damage repair dynamics as contributing to this evolved resistance. Our methodology has potential extrapolation to untested chemicals of ecological concern. Environ Toxicol Chem 2023;42:2040-2053. © 2023 Oak Ridge National Laboratory and The Authors. Environmental Toxicology and Chemistry published by Wiley Periodicals LLC on behalf of SETAC.

Urban evolutionary ecology brings exaptation back into focus

The contribution of pre-existing phenotypic variation to evolution in novel environments has long been appreciated. Nevertheless, evolutionary ecologists have struggled with communicating these aspects of the adaptive process. In 1982, Gould and Vrba proposed terminology to distinguish character states shaped via natural selection for the roles they currently serve ('adaptations') from those shaped under preceding selective regimes ('exaptations'), with the intention of replacing the inaccurate 'preadaptation'. Forty years later, we revisit Gould and Vrba's ideas which, while often controversial, continue to be widely debated and highly cited. We use the recent emergence of urban evolutionary ecology as a timely opportunity to reintroduce the ideas of Gould and Vrba as an integrated framework to understand contemporary evolution in novel environments.

Harnessing plant-microbiome interactions for bioremediation across a freshwater urbanization gradient

Urbanization impacts land, air, and water, creating environmental gradients between cities and rural areas. Urban stormwater delivers myriad co-occurring, understudied, and mostly unregulated contaminants to aquatic ecosystems, causing a pollution gradient. Recipient ecosystems host interacting species that can affect each others' growth and responses to these contaminants. For example, plants and their microbiomes often reciprocally increase growth and contaminant tolerance. Here, we identified ecological variables affecting contaminant fate across an urban-rural gradient using 50 sources of the aquatic plant Lemna minor (duckweed) and associated microbes, and two co-occurring winter contaminants of temperate cities, benzotriazole and salt. We conducted experiments totalling >2,500 independent host-microbe-contaminant microcosms. Benzotriazole and salt negatively affected duckweed growth, but not microbial growth, and duckweeds maintained faster growth with their local, rather than disrupted, microbiota. Benzotriazole transformation products of plant, microbial, and phototransformation pathways were linked to duckweed and microbial growth, and were affected by salt co-contamination, microbiome disruption, and source sites of duckweeds and microbes. Duckweeds from urban sites grew faster and enhanced phytotransformation, but supported less total transformation of benzotriazole. Increasing microbial community diversity correlated with greater removal of benzotriazole, but taxonomic groups may explain shifts across transformation pathways: the genus Aeromonas was linked to increasing phototransformation. Because benzotriazole toxicity could depend on amount and type of in situ transformation, this variation across duckweeds and microbes could be harnessed for better management of urban stormwater. Broadly, our results demonstrate that plant-microbiome interactions harbour manipulable variation for bioremediation applications.

Cardiac physiology and metabolic gene expression during late organogenesis among F. heteroclitus embryo families from crosses between pollution-sensitive and -resistant parents

BACKGROUND

The teleost fish Fundulus heteroclitus inhabit estuaries heavily polluted with persistent and bioaccumulative chemicals. While embryos of parents from polluted sites are remarkably resistant to toxic sediment and develop normally, embryos of parents from relatively clean estuaries, when treated with polluted sediment extracts, are developmentally delayed, displaying deformities characteristic of pollution-induced embryotoxicity. To gain insight into parental effects on sensitive and resistant phenotypes during late organogenesis, we established sensitive, resistant, and crossed embryo families using five female and five male parents from relatively clean and predominantly PAH-polluted estuaries each, measured heart rates, and quantified individual embryo expression of 179 metabolic genes.

RESULTS

Pollution-induced embryotoxicity manifested as morphological deformities, significant developmental delays, and altered cardiac physiology was evident among sensitive embryos resulting from crosses between females and males from relatively clean estuaries. Significantly different heart rates among several geographically unrelated populations of sensitive, resistant, and crossed embryo families during late organogenesis and pre-hatching suggest site-specific adaptive cardiac physiology phenotypes relative to pollution exposure. Metabolic gene expression patterns (32 genes, 17.9%, at p < 0.05; 11 genes, 6.1%, at p < 0.01) among the embryo families indicate maternal pollutant deposition in the eggs and parental effects on gene expression and metabolic alterations.

CONCLUSION

Heart rate differences among sensitive, resistant, and crossed embryos is a reliable phenotype for further explorations of adaptive mechanisms. While metabolic gene expression patterns among embryo families are suggestive of parental effects on several differentially expressed genes, a definitive adaptive signature and metabolic cost of resistant phenotypes is unclear and shows unexpected sensitive-resistant crossed embryo expression profiles. Our study highlights physiological and metabolic gene expression differences during a critical embryonic stage among pollution sensitive, resistant, and crossed embryo families, which may contribute to underlying resistance mechanisms observed in natural F. heteroclitus populations living in heavily contaminated estuaries.

Long-lasting effects of chronic exposure to chemical pollution on the hologenome of the Manila clam

Chronic exposure to pollutants affects natural populations, creating specific molecular and biochemical signatures. In the present study, we tested the hypothesis that chronic exposure to pollutants might have substantial effects on the Manila clam hologenome long after removal from contaminated sites. To reach this goal, a highly integrative approach was implemented, combining transcriptome, genetic and microbiota analyses with the evaluation of biochemical and histological profiles of the edible Manila clam Ruditapes philippinarum, as it was transplanted for 6 months from the polluted area of Porto Marghera (PM) to the clean area of Chioggia (Venice lagoon, Italy). One month post-transplantation, PM clams showed several modifications to its resident microbiota, including an overrepresentation of the opportunistic pathogen Arcobacter spp. This may be related to the upregulation of several immune genes in the PM clams, potentially representing a host response to the increased abundance of deleterious bacteria. Six months after transplantation, PM clams demonstrated a lower ability to respond to environmental/physiological stressors related to the summer season, and the hepatopancreas-associated microbiota still showed different compositions among PM and CH clams. This study confirms that different stressors have predictable effects in clams at different biological levels and demonstrates that chronic exposure to pollutants leads to long-lasting effects on the animal hologenome. In addition, no genetic differentiation between samples from the two areas was detected, confirming that PM and CH clams belong to a single population. Overall, the obtained responses were largely reversible and potentially related to phenotypic plasticity rather than genetic adaptation. The results here presented will be functional for the assessment of the environmental risk imposed by chemicals on an economically important bivalve species.

Association between Chromosome 4 and mercury accumulation in muscle of the three-spined stickleback (Gasterosteus aculeatus)

Anthropogenic stressors, such as pollutants, act as selective factors that can leave measurable changes in allele frequencies in the genome. Metals are of particular concern among pollutants, because of interference with vital biological pathways. We use the three-spined stickleback as a model for adaptation to mercury pollution in natural populations. We collected sticklebacks from 21 locations in Flanders (Belgium), measured the accumulated levels of mercury in the skeletal muscle tissue, and genotyped the fish by sequencing (GBS). The spread of muscle mercury content across locations was considerable, ranging from 21.5 to 327 ng/g dry weight (DW). We then conducted a genome-wide association study (GWAS) between 28,450 single nucleotide polymorphisms (SNPs) and the accumulated levels of mercury, using different approaches. Based on a linear mixed model analysis, the GWAS yielded multiple hits with a single top hit on Chromosome 4, with eight more SNPs suggestive of association. A second approach, a latent factor mixed model analysis, highlighted one single SNP on Chromosome 11. Finally, an outlier test identified one additional SNP on Chromosome 4 that appeared under selection. Out of all ten SNPs we identified as associated with mercury in muscle, three SNPs all located on Chromosome 4 and positioned within a 2.5 kb distance of an annotated gene. Based on these results and the genome coverage of our SNPs, we conclude that the selective effect of mercury pollution in Flanders causes a significant association with at least one locus on Chromosome 4 in three-spined stickleback.

Genetic differentiation in pesticide resistance between urban and rural populations of a nontarget freshwater keystone interactor, Daphnia magna

There is growing evidence that urbanization drives adaptive evolution in response to thermal gradients. One such example is documented in the water flea Daphnia magna. However, organisms residing in urban lentic ecosystems are increasingly exposed to chemical pollutants such as pesticides through run-off and aerial transportation. The extent to which urbanization drives the evolution of pesticide resistance in aquatic organisms and whether this is impacted by warming and thermal adaptation remains limitedly studied. We performed a common garden rearing experiment using multiple clonal lineages originating from five replicated urban and rural D. magna populations, in which we implemented an acute toxicity test exposing neonates (<24h) to either a solvent control or the organophosphate pesticide chlorpyrifos. Pesticide exposures were performed at two temperatures (20°C vs. 24°C) to test for temperature-associated differences in urbanization-driven evolved pesticide resistance. We identified a strong overall effect of pesticide exposure on Daphnia survival probability (-72.8 percentage points). However, urban Daphnia genotypes showed higher survival probabilities compared to rural ones in the presence of chlorpyrifos (+29.7 percentage points). Our experiment did not reveal strong temperature x pesticide or temperature x pesticide x urbanization background effects on survival probability. The here observed evolution of resistance to an organophosphate pesticide is a first indication Daphnia likely also adapts to pesticide pollution in urban areas. Increased pesticide resistance could facilitate their population persistence in urban ponds, and feed back to ecosystem functions, such as top-down control of algae. In addition, adaptive evolution of nontarget organisms to pest control strategies and occupational pesticide use may modulate how pesticide applications affect genetic and species diversity in urban areas.

Genomics of Adaptations in Ungulates

Ungulates are a group of hoofed animals that have long interacted with humans as essential sources of food, labor, clothing, and transportation. These consist of domesticated, feral, and wild species raised in a wide range of habitats and biomes. Given the diverse and extreme environments inhabited by ungulates, unique adaptive traits are fundamental for fitness. The documentation of genes that underlie their genomic signatures of selection is crucial in this regard. The increasing availability of advanced sequencing technologies has seen the rapid growth of ungulate genomic resources, which offers an exceptional opportunity to understand their adaptive evolution. Here, we summarize the current knowledge on evolutionary genetic signatures underlying the adaptations of ungulates to different habitats.

CRISPR/Cas9 mutagenesis reveals a role for ABCB1 in gut immune responses to Vibrio diazotrophicus in sea urchin larvae

The ABC transporter ABCB1 plays an important role in the disposition of xenobiotics. Embryos of most species express high levels of this transporter in early development as a protective mechanism, but its native substrates are not known. Here, we used larvae of the sea urchin Strongylocentrotus purpuratus to characterize the early life expression and role of Sp-ABCB1a, a homolog of ABCB1. The results indicate that while Sp-ABCB1a is initially expressed ubiquitously, it becomes enriched in the developing gut. Using optimized CRISPR/Cas9 gene editing methods to achieve high editing efficiency in the F0 generation, we generated ABCB1a crispant embryos with significantly reduced transporter efflux activity. When infected with the opportunistic pathogen Vibrio diazotrophicus, Sp-ABCB1a crispant larvae demonstrated significantly stronger gut inflammation, immunocyte migration and cytokine Sp-IL-17 induction, as compared with infected control larvae. The results suggest an ancestral function of ABCB1 in host-microbial interactions, with implications for the survival of invertebrate larvae in the marine microbial environment.

The Importance of Eco-evolutionary Potential in the Anthropocene

Humans are dominant global drivers of ecological and evolutionary change, rearranging ecosystems and natural selection. In the present article, we show increasing evidence that human activity also plays a disproportionate role in shaping the eco-evolutionary potential of systems—the likelihood of ecological change generating evolutionary change and vice versa. We suggest that the net outcome of human influences on trait change, ecology, and the feedback loops that link them will often (but not always) be to increase eco-evolutionary potential, with important consequences for stability and resilience of populations, communities, and ecosystems. We also integrate existing ecological and evolutionary metrics to predict and manage the eco-evolutionary dynamics of human-affected systems. To support this framework, we use a simple eco–evo feedback model to show that factors affecting eco-evolutionary potential are major determinants of eco-evolutionary dynamics. Our framework suggests that proper management of anthropogenic effects requires a science of human effects on eco-evolutionary potential.

Pesticide exposure affects reproductive capacity of common toads (Bufo bufo) in a viticultural landscape

Amphibian populations are declining worldwide at alarming rates. Among the large variety of contributing stressors, chemical pollutants like pesticides have been identified as a major factor for this decline. Besides direct effects on aquatic and terrestrial amphibian stages, sublethal effects like impairments in reproduction can affect a population. Therefore, we investigated the reproductive capacity of common toads (Bufo bufo) in the pesticide-intensive viticultural landscape of Palatinate in Southwest Germany along a pesticide gradient. In a semi-field study, we captured reproductively active common toad pairs of five breeding ponds with different pesticide contamination level and kept them in a net cage until spawning. Toads from more contaminated ponds showed an increased fecundity (more eggs) but decreased fertilization rates (fewer hatching tadpoles) as well as lower survival rates and reduced size in Gosner stage 25, suggesting that the higher exposed populations suffer from long-term reproductive impairments. In combination with acute toxicity effects, the detected sublethal effects, which are mostly not addressed in the ecological risk assessment of pesticides, pose a serious threat on amphibian populations in agricultural landscapes.

Lower bioenergetic costs but similar immune responsiveness under a heat wave in urban compared to rural damselflies

There is mounting evidence that the widespread phenotypic changes in response to urbanization may reflect adaptations caused by rapid evolutionary processes driven by urban-related stressors. Compared to increased habitat fragmentation and pollution, adaptations towards another typical urban-related stressor, that is higher and longer lasting very high temperatures (heat waves), are much less studied. Notably, the sensitivities to heat waves of life-history traits and important fitness-related physiological traits such as immune responsiveness and bioenergetic variables (energy availability, energy consumption and their balance) have never been contrasted between urban and rural populations. By conducting a laboratory common-garden experiment, we compared effects of a simulated heat wave on life history (survival and growth rate), immune responsiveness and bioenergetic variables between three urban and three rural populations of the damselfly Coenagrion puella. Because energy-mediated trade-off patterns may only be detected under energetically costly manipulations, all larvae were immune-challenged by simulating ectoparasitism by water mites. As expected, the simulated heat wave caused negative effects on nearly all response variables. The immune responsiveness, on the other hand, increased under the heat wave, consistent with a trade-off pattern between immune function and growth, and this similarly between urban and rural populations. A key finding was that urban larvae suffered less from the simulated heat wave compared to the rural larvae in terms of a lower heat wave-induced depletion in energy availability. This suggests an adaptation of urban populations to better cope with the stronger and more frequent heat waves in cities. Notably, this urbanization-driven evolution in the bioenergetic variables was not apparent in the absence of a heat wave. Given that changes in energy budgets have strong fitness consequences, our findings suggest that the evolved higher ability to cope with heat waves is fundamental for the survival of urban damselfly populations.

Ecological and Evolutionary Consequences of Anticancer Adaptations

Cellular cheating leading to cancers exists in all branches of multicellular life, favoring the evolution of adaptations to avoid or suppress malignant progression, and/or to alleviate its fitness consequences. Ecologists have until recently largely neglected the importance of cancer cells for animal ecology, presumably because they did not consider either the potential ecological or evolutionary consequences of anticancer adaptations. Here, we review the diverse ways in which the evolution of anticancer adaptations has significantly constrained several aspects of the evolutionary ecology of multicellular organisms at the cell, individual, population, species, and ecosystem levels and suggest some avenues for future research.

Tributyltin disrupts fin development in Fundulus heteroclitus from both PCB-sensitive and resistant populations: Investigations of potential interactions between AHR and PPARγ

Tributyltin (TBT) and dioxin-like polychlorinated biphenyls (PCBs) are environmental contaminants that are highly toxic to fish and co-occur in New Bedford Harbor (NBH), an estuarine Superfund site located in Massachusetts, USA. Atlantic killifish (Fundulus heteroclitus) that reside in NBH (and other highly contaminated sites along the east coast of the United States) have developed resistance to activation of the aryl hydrocarbon receptor (AHR) pathway and the toxicity of dioxin-like chemicals, such as 3,3',4,4',5-pentachlorobiphenyl, PCB126. In many biological systems, TBT disregulates adipose and bone development via the PPARγ-RXR pathway; AHR activation also disrupts adipose and bone homeostasis, potentially through molecular crosstalk between AHR and PPARγ. However, little is known about how co-exposure and the interaction of these pathways modulate the toxicological effects of these contaminants. Here, we tested the hypotheses that TBT would induce teratogenesis in killifish via activation of PPARγ and that PCB126 co-exposure would suppress PPARγ pathway activation in PCB-sensitive killifish from a reference site (Scorton Creek, SC, PCB-sensitive) but not in PCB-tolerant NBH killifish. Killifish embryos from both populations exposed to TBT (50 and 100 nM) displayed caudal fin deformities. TBT did not change the expression of pparg or its target genes related to adipogenesis (fabp11a and fabp1b) in either population. However, expression of osx/sp7, an osteoblast marker gene, and col2a1b, a chondroblast marker gene, was significantly suppressed by TBT only in SC killifish. An RXR-specific agonist, but not a PPARγ-specific agonist, induced caudal fin deformities like those observed in TBT-treated embryos. PCB126 did not induce caudal fin deformities and did not exacerbate TBT-induced fin deformities. Further, PCB126 increased expression of pparg in SC embryos and not NBH embryos, but did not change the expression of fabp1b. Taken together, these results suggest that in killifish embryos the PPARγ pathway is regulated in part by AHR, but is minimally active at least in this early life stage. In killifish, RXR activation, rather than PPARγ activation, appears to be the mechanism by which TBT induces caudal fin teratogenicity, which is not modulated by AHR responsiveness.

Rapid resistance to pesticide control is predicted to evolve in an invasive fish

Xenobiotic resistance is commonly found in species with short generation times such as bacteria, annual plants, and insects. Nevertheless, the fundamental evolutionary principles that govern the spread of resistance alleles hold true for species with longer generation times. One such example could occur with sea lamprey (Petromyzon marinus), a parasitic invasive species in the Laurentian Great Lakes that decimated native fish populations prior to its control with the pesticide 3-trifluoromethyl-4-nitrophenol (TFM). Since the 1950s, tributaries have been treated annually with TFM, where treatments effectively remove most, but not all, larval sea lamprey. We developed an eco-genetic model of sea lamprey to examine factors affecting the evolution of resistance and found that resistance alleles rapidly rise to fixation after 40-80 years of treatment, despite the species' relatively long generation time (4-7 years). The absence of natal homing allows resistant individuals to spread quickly throughout the entire system, but also makes the early detection of resistance challenging. High costs of resistance and density independent reproduction can delay, but not prevent, the onset of resistance. These results illustrate that sea lamprey have the potential to evolve resistance to their primary control agent in the near future, highlighting the urgent need for alternative controls.

Adaptive introgression enables evolutionary rescue from extreme environmental pollution

Radical environmental change that provokes population decline can impose constraints on the sources of genetic variation that may enable evolutionary rescue. Adaptive toxicant resistance has rapidly evolved in Gulf killifish (Fundulus grandis) that occupy polluted habitats. We show that resistance scales with pollution level and negatively correlates with inducibility of aryl hydrocarbon receptor (AHR) signaling. Loci with the strongest signatures of recent selection harbor genes regulating AHR signaling. Two of these loci introgressed recently (18 to 34 generations ago) from Atlantic killifish (F. heteroclitus). One introgressed locus contains a deletion in AHR that confers a large adaptive advantage [selection coefficient (s) = 0.8]. Given the limited migration of killifish, recent adaptive introgression was likely mediated by human-assisted transport. We suggest that interspecies connectivity may be an important source of adaptive variation during extreme environmental change.

The role of intraspecies variation in fish neurobehavioral and neuropharmacological phenotypes in aquatic models

Intraspecies variation is common in both clinical and animal research of various brain disorders. Relatively well-studied in mammals, intraspecies variation in aquatic fish models and its role in their behavioral and pharmacological responses remain poorly understood. Like humans and mammals, fishes show high variance of behavioral and drug-evoked responses, modulated both genetically and environmentally. The zebrafish (Danio rerio) has emerged as a particularly useful model organism tool to access neurobehavioral and drug-evoked responses. Here, we discuss recent findings and the role of the intraspecies variance in neurobehavioral, pharmacological and toxicological studies utilizing zebrafish and other fish models. We also critically evaluate common sources of intraspecies variation and outline potential strategies to improve data reproducibility and translatability.

Altered lipid homeostasis in a PCB-resistant Atlantic killifish (Fundulus heteroclitus) population from New Bedford Harbor, MA, U.S.A

Sentinel species such as the Atlantic killifish (Fundulus heteroclitus) living in urban waterways can be used as toxicological models to understand impacts of environmental metabolism disrupting compound (MDC) exposure on both wildlife and humans. Exposure to MDCs is associated with increased risk of metabolic syndrome, including impaired lipid and glucose homeostasis, adipogenesis, appetite control, and basal metabolism. MDCs are ubiquitous in the environment, including in aquatic environments. New Bedford Harbor (NBH), Massachusetts is polluted with polychlorinated biphenyls (PCBs), and, as we show for the first time, tin (Sn). PCBs and organotins are ligands for two receptor systems known to regulate lipid homeostasis, the aryl hydrocarbon receptor (AHR) and the peroxisome proliferator-activated receptors (PPARs), respectively. In the current study, we compared lipid homeostasis in laboratory-reared killifish from NBH (F2) and a reference location (Scorton Creek, Massachusetts; F1 and F2) to evaluate how adaptation to local conditions may influence responses to MDCs. Adult killifish from each population were exposed to 3,3',4,4',5-pentachlorobiphenyl (PCB126, dioxin-like), 2,2',4,4',5,5'-hexachlorobiphenyl (PCB153, non-dioxin-like), or tributyltin (TBT, a PPARγ ligand) by a single intraperitoneal injection and analyzed after 3 days. AHR activation was assessed by measuring cyp1a mRNA expression. Lipid homeostasis was evaluated phenotypically by measuring liver triglycerides and organosomatic indices, and at the molecular level by measuring the mRNA expression of pparg and ppara and a target gene for each receptor. Acute MDC exposure did not affect phenotypic outcomes. However, overall NBH killifish had higher liver triglycerides and adiposomatic indices than SC killifish. Both season and population were significant predictors of the lipid phenotype. Acute MDC exposure altered hepatic gene expression only in male killifish from SC. PCB126 exposure induced cyp1a and pparg, whereas PCB153 exposure induced ppara. TBT exposure did not induce ppar-dependent pathways. Comparison of lipid homeostasis in two killifish populations extends our understanding of how MDCs act on fish and provides a basis to infer adaptive benefits of these differences in the wild.

Characterization of the Fundulus heteroclitus embryo transcriptional response and development of a gene expression-based fingerprint of exposure for the alternative flame retardant, TBPH (bis (2-ethylhexyl)-tetrabromophthalate)

Although alternative Flame Retardant (FR) chemicals are expected to be safer than the legacy FRs they replace, their risks to human health and the environment are often poorly characterized. This study used a small volume, fish embryo system to reveal potential mechanisms of action and diagnostic exposure patterns for TBPH (bis (2-ethylhexyl)-tetrabromophthalate), a component of several widely-used commercial products. Two different concentration of TBPH were applied to sensitive early life stages of an ecologically important test species, Fundulus heteroclitus (Atlantic killifish), with a well-annotated genome. Exposed fish embryos were sampled for transcriptomics or chemical analysis of parent compound and primary metabolite or observed for development and survival through larval stage. Global transcript profiling using RNA-seq was conducted (n = 16 per treatment) to provide a non-targeted and statistically robust approach to characterize TBPH gene expression patterns. Transcriptomic analysis revealed a dose-response in the expression of genes associated with a surprisingly limited number of biological pathways, but included the aryl hydrocarbon receptor signal transduction pathway, which is known to respond to several toxicologically-important chemical classes. A transcriptional fingerprint using Random Forests was developed that was able to perfectly discriminate exposed vs. non-exposed individuals in test sets. These results suggest that TBPH has a relatively low potential for developmental toxicity (at least in fishes), despite concerns related to its structural similarities to endocrine disrupting chemicals and that the early life stage Fundulus system may provide a convenient test system for exposure characterization. More broadly, this study advances the usefulness of a biological testing and analysis system utilizing non-targeted transcriptomics profiling and early developmental endpoints that complements current screening methods to characterize chemicals of ecological and human health concern.

Toads phenotypically adjust their chemical defences to anthropogenic habitat change

Despite the well-documented effects of human-induced environmental changes on the morphology, physiology, behaviour and life history of wild animals, next to nothing is known about how anthropogenic habitats influence anti-predatory chemical defence, a crucial fitness component of many species. We investigated the amount and composition of defensive toxins in adult common toads (Bufo bufo) captured in natural, agricultural and urban habitats, and in their offspring raised in a common-garden experiment. We found that, compared to toads captured from natural habitats, adults from both types of anthropogenic habitats had larger toxin glands (parotoids) and their toxin secretion contained higher concentrations of bufagenins, the more potent class of bufadienolide toxins. Furthermore, urban toads had lower concentrations of bufotoxins, the compounds with lower toxicity. None of these differences were present in the captive-raised juveniles; instead, toadlets originating from agricultural habitats had smaller parotoids and lower bufotoxin concentrations. These results suggest that toads' chemical defences respond to the challenges of anthropogenic environments via phenotypic plasticity. These responses may constitute non-adaptive consequences of pollution by endocrine-disrupting chemicals as well as adaptive adjustments to the altered predator assemblages of urban and agricultural habitats.

Individual variation in aquatic toxicology: Not only unwanted noise

The mean value of any parameter and its changes are usually discussed, when ecotoxicological studies are carried out. However, also the variation of any parameter and its changes can be important components of the responses to environmental contamination. Although the homogeneity of variances is commonly tested, testing is done for the use of correct statistical methods, not because of exploring the possibility that variability and its changes could be important components of environmental responses. We evaluated recent aquatic toxicological literature and found that in the majority of articles indicating that homogeneity of variances was tested and giving the result of testing, the assumption of homogeneity was not fulfilled. Further, it was observed that in some studies experimental treatment clearly affected the variability. In this commentary we discuss the reasons for variability: measurement errors, experimental design, genetic heterogeneity and phenotypic plasticity, and conclude that even after accounting for experimental design and genetic makeup significant variability remains. This plasticity may change in environmental responses as suggested by a hypothetical example, and as confirmed by experimental data. As a consequence, the changes of variability can be significant, even when the means do not differ. Because of this, variability and its changes should always be analysed and reported. This will be easy, since the datasets are exactly the same for comparing the variances and means, and as normally variances are tested for homogeneity. It is likely that much new information about the responses of organisms to environmental contamination will be obtained. However, the present journal practises tend to discourage one from concentrating on anything but the mean. In contrast, we think it is imperative that variability is always included as an endpoint in data analysis in the future.

Urbanization threaten the pollination of Gentiana dahurica

With rapid spread of the urbanization, many environmental factors, such as climate, soil pH and nutrients have been changed. However, the plant pollination affected by urbanization was seldom conducted. Here, we studied the flower visitation rates, seed production, pollen limitation and flower morphological characters of Gentiana dahurica at 3 populations along an urban-peri-urban gradient around Xi'ning over 4 consecutive years, aiming to test the effects of urbanization on plant pollination service. Our results showed that the pollinator visit frequencies, interannual stability of pollinator assemblages and visit frequencies declined with the intensification of urbanization. As urbanization intensified, plant borne more flowers and the flower morphological sizes became "longer" (the length of flowers, filaments and styles were increased, but the width of flowers kept stable at the 3 populations); the flower duration, especially the female phase duration prolonged. The seed-set ratio of G. dahurica in natural condition decreased and more severe pollen limitation occurred in more urbanized populations. Also, an interannual variation of seed-set ratio and index of pollen limitation (IPL), which related with the variation of pollinator visit frequencies, were found in this study. These results suggest that the pollination service can be threatened by urbanization over a long-time interval for G. dahurica. This finding highlights the importance of pollinator affections acting on plant pollination system. Additionally, as pollinator assemblages and visit frequencies interannually changed, a long-time scale observation is needed to understand the plant-pollinator relationships.

Memory effect of arsenic-induced cellular response and its influences on toxicity of titanium dioxide nanoparticle

Toxicity of arsenic (As) has been widely characterized. However, few studies focus on whether cell responses induced by As at nontoxic concentration could be inherited and further change cell tolerance to another pollutant. In this study, human A549 and HeLa cells were exposed to As at nontoxic concentrations for 10 or 15 passages, then the cells were recovered in the cell medium without As. At 25th passage, residual As in both type of cells was completely removed through the recovery process. And no abnormity in cell viability was identified in both type of cells between control and As-treated groups. Above results indicated that As exposure-recovery treatment had limited influence on phenotype of the cells. However, gene expression profiles determined by high-throughput sequencing showed that As exposure-recovery treatment induced similar expression modification of genes related to inflammation, oxidative stress and epigenetic modulation in the A549 and HeLa cells after recovery of 10 or 15 passages, indicating that As-induced cellular responses have been partially memorized at transcriptional level. The memory effect might play key roles in increased tolerance of the A549 and HeLa cells to adverse effects (cell viability, intracellular reactive oxygen species (ROS) generation and plasma membrane damage) induced by titanium dioxide nanoparticles (as representative pollutant). This study shed new lights on toxic effects induced by As at nontoxic concentration, which is useful for risk assessment of combined effects of As and other pollutants.

Hepatic metabolite profiling of polychlorinated biphenyl (PCB)-resistant and sensitive populations of Atlantic killifish (Fundulus heteroclitus)

Atlantic killifish inhabiting polluted sites along the east coast of the U.S. have evolved resistance to toxic effects of contaminants. One such contaminated site is the Acushnet River estuary, near New Bedford Harbor (NBH), Massachusetts, which is characterized by very high PCB concentrations in the sediments and in the tissues of resident killifish. Though killifish at this site appear to be thriving, the metabolic costs of survival in a highly contaminated environment are not well understood. In this study we compared the hepatic metabolite profiles of resistant (NBH) and sensitive populations (Scorton Creek (SC), Sandwich, MA) using a targeted metabolomics approach in which polar metabolites were extracted from adult fish livers and quantified. Our results revealed differences in the levels of several metabolites between fish from the two sites. The majority of these metabolites are associated with one-carbon metabolism, an important pathway that supports multiple physiological processes including DNA and protein methylation, nucleic acid biosynthesis and amino acid metabolism. We measured the gene expression of DNA methylation (DNA methyltransferase 1, dnmt1) and demethylation genes (Ten-Eleven Translocation (TET) genes) in the two populations, and observed lower levels of dnmt1 and higher levels of TET gene expression in the NBH livers, suggesting possible differences in DNA methylation profiles. Consistent with this, the two populations differed significantly in the levels of 5-methylcytosine and 5-hydroxymethylcytosine nucleotides. Overall, our results suggest that the unique hepatic metabolite signatures observed in NBH and SC reflect the adaptive mechanisms for survival in their respective habitats.

Variable gene transcription underlies phenotypic convergence of hypoxia tolerance in sculpins

Background

The degree by which mechanisms underlying phenotypic convergence are similar among taxa depends on the number of evolutionary paths available for selection to act upon. Likelihood of convergence will be influenced by an interplay of factors such as genetic architecture, phylogenetic history and population demography. To determine if there is convergence or divergence in mechanisms underlying phenotypic similarity, we assessed whether gene transcription patterns differed among species with similar levels of hypoxia tolerance.

Results

Three species of marine fish from the superfamily Cottoidea (smoothhead sculpin [Artedius lateralis], sailfin sculpin [Nautichthys oculofasciatus] and Pacific staghorn sculpin [Leptocottus armatus]), all of which have previously been shown to share the same level of hypoxia tolerance, were exposed to short-(8 h) and longer-term (72 h) hypoxia and mRNA transcripts were assessed using a custom microarray. We examined hypoxia-induced transcription patterns in metabolic and protein production pathways and found that a high proportion of genes associated with these biological processes showed significant differences among the species. Specifically, the data suggest that the smoothhead sculpin, unlike the sailfin sculpin and the Pacific staghorn sculpin, relied on amino acid degradation rather than glycolysis or fatty acid oxidation to generate ATP during hypoxia exposure. There was also variation across the species in the transcription of genes involved in protein production (e.g. mRNA processing and protein translation), such that it increased in the smoothhead sculpin, decreased in the sailfin sculpin and was variable in the Pacific staghorn sculpin.

Conclusions

Changes in metabolic and protein production pathways are part of the key responses of fishes to exposures to environmental hypoxia. Yet, species with similar overall hypoxia tolerance exhibited different transcriptional responses in these pathways, indicating flexibility and complexity of interactions in the evolution of the mechanisms underlying the hypoxia tolerance phenotype. The variation in the hypoxia-induced transcription of genes across species with similar hypoxia tolerance suggests that similar whole-animal phenotypes can emerge from divergent evolutionary paths that may affect metabolically important functions.

Electronic supplementary material

The online version of this article (10.1186/s12862-018-1275-1) contains supplementary material, which is available to authorized users.

The developmental and genetic trajectory of coloration in the guppy (Poecilia reticulata)

Examining the association between trait variation and development is crucial for understanding the evolution of phenotypic differences. Male guppy ornamental caudal fin coloration is one trait that shows a striking degree of variation within and between guppy populations. Males initially have no caudal fin coloration, then gradually develop it as they reach sexual maturity. For males, there is a trade-off between female preference for caudal fin coloration and increased visibility to predators. This trade-off may reach unique endpoints in males from different predation regimes. Caudal fin coloration includes black melanin, orange/yellow pteridines or carotenoids, and shimmering iridescence. This study examined the phenotypic trajectory and genetics associated with color development. We found that black coloration always developed first, followed by orange/yellow, then iridescence. The ordering and timing of color appearance was the same regardless of predation regime. The increased expression of melanin synthesis genes correlated well with the visual appearance of black coloration, but there was no correlation between carotenoids or pteridine synthesis gene expression and the appearance of orange/yellow. The lack of orange/yellow coloration in earlier male caudal fin developmental stages may be due to reduced expression of genes underlying the development of orange/yellow xanthophores.

Endocrine disruptors in breeding ponds and reproductive health of toads in agricultural, urban and natural landscapes

Many chemical pollutants have endocrine disrupting effects which can cause lifelong reproductive abnormalities in animals. Amphibians are the most threatened group of vertebrates, but there is little information on the nature and quantity of pollutants occurring in typical amphibian breeding habitats and on the reproductive capacities of amphibian populations inhabiting polluted areas. In this study we investigated the occurrence and concentrations of endocrine disrupting chemicals in the water and sediment of under-studied amphibian breeding habitats in natural, agricultural and urbanized landscapes. Also, we captured reproductively active common toads (Bufo bufo) from these habitats and let them spawn in a 'common garden' to assess among-population differences in reproductive capacity. Across 12 ponds, we detected 41 out of the 133 contaminants we screened for, with unusually high concentrations of glyphosate and carbamazepine. Levels of polycyclic aromatic hydrocarbons, nonylphenol and bisphenol-A increased with urban land use, whereas levels of organochlorine and triazine pesticides and sex hormones increased with agricultural land use. Toads from all habitats had high fecundity, fertilization rate and offspring viability, but the F1 generation originating from agricultural and urban ponds had reduced development rates and lower body mass both as larvae and as juveniles. Females with small clutch mass produced thicker jelly coat around their eggs if they originated from agricultural and urban ponds compared with natural ponds. These results suggest that the observed pollution levels did not compromise reproductive potential in toads, but individual fitness and population viability may be reduced in anthropogenically influenced habitats, perhaps due to transgenerational effects and/or costs of tolerance to chemical contaminants.

Phenotypic selection on floral traits in an urban landscape

Native species are increasingly living in urban landscapes associated with abiotic and biotic changes that may influence patterns of phenotypic selection. However, measures of selection in urban and non-urban environments, and exploration of the mechanisms associated with such changes, are uncommon. Plant-animal interactions have played a central role in the evolution of flowering plants and are sensitive to changes in the urban landscape, and thus provide opportunities to explore how urban environments modify selection. We evaluated patterns of phenotypic selection on the floral and resistance traits of Gelsemium sempervirens in urban and non-urban sites. The urban landscape had increased florivory and decreased pollen receipt, but showed only modest differences in patterns of selection. Directional selection for one trait, larger floral display size, was stronger in urban compared to non-urban sites. Neither quadratic nor correlational selection significantly differed between urban and non-urban sites. Pollination was associated with selection for larger floral display size in urban compared to non-urban sites, due to the differences in the translation of pollination into seeds rather than pollinator selectivity. Thus, our data suggest that urban landscapes may not result in sweeping differences in phenotypic selection but rather modest differences for some traits, potentially mediated by species interactions.

Stable genetic structure and connectivity in pollution-adapted and nearby pollution-sensitive populations of Fundulus heteroclitus

Populations of the non-migratory estuarine fish Fundulus heteroclitus inhabiting the heavily polluted New Bedford Harbour (NBH) estuary have shown inherited tolerance to local pollutants introduced to their habitats in the past 100 years. Here we examine two questions: (i) Is there pollution-driven selection on the mitochondrial genome across a fine geographical scale? and (ii) What is the pattern of migration among sites spanning a strong pollution gradient? Whole mitochondrial genomes were analysed for 133 F. heteroclitus from seven nearby collection sites: four sites along the NBH pollution cline (approx. 5 km distance), which had pollution-adapted fish, as well as one site adjacent to the pollution cline and two relatively unpolluted sites about 30 km away, which had pollution-sensitive fish. Additionally, we used microsatellite analyses to quantify genetic variation over three F. heteroclitus generations in both pollution-adapted and sensitive individuals collected from two sites at two different time points (1999/2000 and 2007/2008). Our results show no evidence for a selective sweep of mtDNA in the polluted sites. Moreover, mtDNA analyses revealed that both pollution-adapted and sensitive populations harbour similar levels of genetic diversity. We observed a high level of non-synonymous mutations in the most polluted site. This is probably associated with a reduction in Ne and concomitant weakening of purifying selection, a demographic expansion following a pollution-related bottleneck or increased mutation rates. Our demographic analyses suggest that isolation by distance influences the distribution of mtDNA genetic variation between the pollution cline and the clean populations at broad spatial scales. At finer scales, population structure is patchy, and neither spatial distance, pollution concentration or pollution tolerance is a good predictor of mtDNA variation. Lastly, microsatellite analyses revealed stable population structure over the last decade.

Gene expression signatures of mating system evolution

The diversity of mating systems among animals is astounding. Importantly, similar mating systems have evolved even across distantly related taxa. However, our understanding of the mechanisms underlying these convergently evolved phenotypes is limited. Here, we examine on a genomic scale the neuromolecular basis of social organization in cichlids of the tribe Ectodini from Lake Tanganyika. Using field-collected males and females of four closely related species representing two independent evolutionary transitions from polygyny to monogamy, we take a comparative transcriptomic approach to test the hypothesis that these independent transitions have recruited similar gene sets. Our results demonstrate that while lineage and species exert a strong influence on neural gene expression profiles, social phenotype can also drive gene expression evolution. Specifically, 331 genes (∼6% of those assayed) were associated with monogamous mating systems independent of species or sex. Among these genes, we find a strong bias (4:1 ratio) toward genes with increased expression in monogamous individuals. A highly conserved nonapeptide system known to be involved in the regulation of social behavior across animals was not associated with mating system in our analysis. Overall, our findings suggest deep molecular homologies underlying the convergent or parallel evolution of monogamy in different cichlid lineages of Ectodini.

Evolution of life in urban environments

Our planet is an increasingly urbanized landscape, with over half of the human population residing in cities. Despite advances in urban ecology, we do not adequately understand how urbanization affects the evolution of organisms, nor how this evolution may affect ecosystems and human health. Here, we review evidence for the effects of urbanization on the evolution of microbes, plants, and animals that inhabit cities. Urbanization affects adaptive and nonadaptive evolutionary processes that shape the genetic diversity within and between populations. Rapid adaptation has facilitated the success of some native species in urban areas, but it has also allowed human pests and disease to spread more rapidly. The nascent field of urban evolution brings together efforts to understand evolution in response to environmental change while developing new hypotheses concerning adaptation to urban infrastructure and human socioeconomic activity. The next generation of research on urban evolution will provide critical insight into the importance of evolution for sustainable interactions between humans and our city environments.

Turning natural adaptations to oncogenic factors into an ally in the war against cancer

Both field and experimental evolution studies have demonstrated that organisms naturally or artificially exposed to environmental oncogenic factors can, sometimes rapidly, evolve specific adaptations to cope with pollutants and their adverse effects on fitness. Although numerous pollutants are mutagenic and carcinogenic, little attention has been given to exploring the extent to which adaptations displayed by organisms living in oncogenic environments could inspire novel cancer treatments, through mimicking the processes allowing these organisms to prevent or limit malignant progression. Building on a substantial knowledge base from the literature, we here present and discuss this progressive and promising research direction, advocating closer collaboration between the fields of medicine, ecology, and evolution in the war against cancer.

Urban driven phenotypic changes: empirical observations and theoretical implications for eco-evolutionary feedback

Emerging evidence that cities drive micro-evolution raises the question of whether rapid urbanization of Earth might impact ecosystems by causing systemic changes in functional traits that regulate urban ecosystems' productivity and stability. Intraspecific trait variation-variation in organisms' morphological, physiological or behavioural characteristics stemming from genetic variability and phenotypic plasticity-has significant implications for ecological functions such as nutrient cycling and primary productivity. While it is well established that changes in ecological conditions can drive evolutionary change in species' traits that, in turn, can alter ecosystem function, an understanding of the reciprocal and simultaneous processes associated with such interactions is only beginning to emerge. In urban settings, the potential for rapid trait change may be exacerbated by multiple selection pressures operating simultaneously. This paper reviews evidence on mechanisms linking urban development patterns to rapid phenotypic changes, and differentiates phenotypic changes for which there is evidence of micro-evolution versus phenotypic changes which may represent plasticity. Studying how humans mediate phenotypic trait changes through urbanization could shed light on fundamental concepts in ecological and evolutionary theory. It can also contribute to our understanding of eco-evolutionary feedback and provide insights for maintaining ecosystem function over the long term.This article is part of the themed issue 'Human influences on evolution, and the ecological and societal consequences'.

Parallel habitat acclimatization is realized by the expression of different genes in two closely related salamander species (genus Salamandra)

The utilization of similar habitats by different species provides an ideal opportunity to identify genes underlying adaptation and acclimatization. Here, we analysed the gene expression of two closely related salamander species: Salamandra salamandra in Central Europe and Salamandra infraimmaculata in the Near East. These species inhabit similar habitat types: 'temporary ponds' and 'permanent streams' during larval development. We developed two species-specific gene expression microarrays, each targeting over 12 000 transcripts, including an overlapping subset of 8331 orthologues. Gene expression was examined for systematic differences between temporary ponds and permanent streams in larvae from both salamander species to establish gene sets and functions associated with these two habitat types. Only 20 orthologues were associated with a habitat in both species, but these orthologues did not show parallel expression patterns across species more than expected by chance. Functional annotation of a set of 106 genes with the highest effect size for a habitat suggested four putative gene function categories associated with a habitat in both species: cell proliferation, neural development, oxygen responses and muscle capacity. Among these high effect size genes was a single orthologue (14-3-3 protein zeta/YWHAZ) that was downregulated in temporary ponds in both species. The emergence of four gene function categories combined with a lack of parallel expression of orthologues (except 14-3-3 protein zeta) suggests that parallel habitat adaptation or acclimatization by larvae from S. salamandra and S. infraimmaculata to temporary ponds and permanent streams is mainly realized by different genes with a converging functionality.

Incorporating evolutionary insights to improve ecotoxicology for freshwater species

Ecotoxicological studies have provided extensive insights into the lethal and sublethal effects of environmental contaminants. These insights are critical for environmental regulatory frameworks, which rely on knowledge of toxicity for developing policies to manage contaminants. While varied approaches have been applied to ecotoxicological questions, perspectives related to the evolutionary history of focal species or populations have received little consideration. Here, we evaluate chloride toxicity from the perspectives of both macroevolution and contemporary evolution. First, by mapping chloride toxicity values derived from the literature onto a phylogeny of macroinvertebrates, fish, and amphibians, we tested whether macroevolutionary relationships across species and taxa are predictive of chloride tolerance. Next, we conducted chloride exposure tests for two amphibian species to assess whether potential contemporary evolutionary change associated with environmental chloride contamination influences chloride tolerance across local populations. We show that explicitly evaluating both macroevolution and contemporary evolution can provide important and even qualitatively different insights from those obtained via traditional ecotoxicological studies. While macroevolutionary perspectives can help forecast toxicological end points for species with untested sensitivities, contemporary evolutionary perspectives demonstrate the need to consider the environmental context of exposed populations when measuring toxicity. Accounting for divergence among populations of interest can provide more accurate and relevant information related to the sensitivity of populations that may be evolving in response to selection from contaminant exposure. Our data show that approaches accounting for and specifically examining variation among natural populations should become standard practice in ecotoxicology.

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.

Complexities of gene expression patterns in natural populations of an extremophile fish (Poecilia mexicana, Poeciliidae)

Variation in gene expression can provide insights into organismal responses to environmental stress and physiological mechanisms mediating adaptation to habitats with contrasting environmental conditions. We performed an RNA-sequencing experiment to quantify gene expression patterns in fish adapted to habitats with different combinations of environmental stressors, including the presence of toxic hydrogen sulphide (H2 S) and the absence of light in caves. We specifically asked how gene expression varies among populations living in different habitats, whether population differences were consistent among organs, and whether there is evidence for shared expression responses in populations exposed to the same stressors. We analysed organ-specific transcriptome-wide data from four ecotypes of Poecilia mexicana (nonsulphidic surface, sulphidic surface, nonsulphidic cave and sulphidic cave). The majority of variation in gene expression was correlated with organ type, and the presence of specific environmental stressors elicited unique expression differences among organs. Shared patterns of gene expression between populations exposed to the same environmental stressors increased with levels of organismal organization (from transcript to gene to physiological pathway). In addition, shared patterns of gene expression were more common between populations from sulphidic than populations from cave habitats, potentially indicating that physiochemical stressors with clear biochemical consequences can constrain the diversity of adaptive solutions that mitigate their adverse effects. Overall, our analyses provided insights into transcriptional variation in a unique system, in which adaptation to H2 S and darkness coincide. Functional annotations of differentially expressed genes provide a springboard for investigating physiological mechanisms putatively underlying adaptation to extreme environments.

Biomarker response and hypothalamus-pituitary-interrenal axis functioning in Arctic charr from Bjørnøya (74°30' N), Norway, with high levels of organohalogenated compounds

The populations of Arctic charr (Salvelinus alpinus) residing in Lake Ellasjøen at Bjørnøya Island in the Norwegian Arctic (74° 30'N, 19° 00'E) possess substantially higher levels of organohalogenated compounds (strongly dominated by polychlorinated biphenyls, PCBs) than conspecifics residing in other, proximate lakes on the island. In the present study we sampled large (<400g), immature charr from Lake Ellasjøen (high PCB levels) and Lake Laksvatn (reference lake, low PCB levels) by hook and line for an immediate blood sampling, and blood and tissue sampling after a 1h confinement stressor. This was done in order to investigate possible effects of pollutants on an acute stress performance in a high-latitude fish species by comparing muscle PCB levels, hepatic cytochrome P4501A (CYP1A) biomarker activation and functioning of the hypothalamus-pituitary-interrenal (HPI) axis between these two populations of Arctic charr. As expected sum PCB muscle levels were 8-fold higher on a wet weigh basis, and 19-fold higher on a lipid weight basis, in charr from Ellasjøen than in charr from Laksvatn. This was accompanied by a 3.5-fold higher liver cyp1a mRNA abundance in the Ellasjøen charr compared to Laksvatn charr. Brain transcript levels encoding glucocorticoid receptor 1 and 2 (GR2) and corticotropin-releasing factor, and pituitary transcript levels encoding GR2 and proopiomelanocortin A1 and A2 were higher in Ellasjøen charr than in Laksvatn charr, while interrenal transcript levels encoding melanocortin 2 receptor and steroidogenic acute regulatory protein were lower. There were no differences in plasma cortisol concentration between the two charr populations immediately after capture and one hour after confinement. The strong biomarker response to OHCs and altered mRNA abundances of key genes related to HPI axis functioning in the Ellasjøen charr suggest endocrine disruptive effects of OHCs in this charr population. Possible ecological implications are not known, but it cannot be excluded that a slower growth rate in Ellasjøen charr compared to Laksvatn charr due to an increased metabolic demand associated with the activation of xenobiotic defense and detoxification systems may have contributed to the lower body mass of Ellasjøen charr compared to Laksvatn charr.

When evolution is the solution to pollution: Key principles, and lessons from rapid repeated adaptation of killifish (Fundulus heteroclitus) populations

For most species, evolutionary adaptation is not expected to be sufficiently rapid to buffer the effects of human‐mediated environmental changes, including environmental pollution. Here we review how key features of populations, the characteristics of environmental pollution, and the genetic architecture underlying adaptive traits, may interact to shape the likelihood of evolutionary rescue from pollution. Large populations of Atlantic killifish (Fundulus heteroclitus) persist in some of the most contaminated estuaries of the United States, and killifish studies have provided some of the first insights into the types of genomic changes that enable rapid evolutionary rescue from complexly degraded environments. We describe how selection by industrial pollutants and other stressors has acted on multiple populations of killifish and posit that extreme nucleotide diversity uniquely positions this species for successful evolutionary adaptation. Mechanistic studies have identified some of the genetic underpinnings of adaptation to a well‐studied class of toxic pollutants; however, multiple genetic regions under selection in wild populations seem to reflect more complex responses to diverse native stressors and/or compensatory responses to primary adaptation. The discovery of these pollution‐adapted killifish populations suggests that the evolutionary influence of anthropogenic stressors as selective agents occurs widely. Yet adaptation to chemical pollution in terrestrial and aquatic vertebrate wildlife may rarely be a successful “solution to pollution” because potentially adaptive phenotypes may be complex and incur fitness costs, and therefore be unlikely to evolve quickly enough, especially in species with small population sizes.

Applying evolutionary genetics to developmental toxicology and risk assessment

Evolutionary thinking continues to challenge our views on health and disease. Yet, there is a communication gap between evolutionary biologists and toxicologists in recognizing the connections among developmental pathways, high-throughput screening, and birth defects in humans. To increase our capability in identifying potential developmental toxicants in humans, we propose to apply evolutionary genetics to improve the experimental design and data interpretation with various in vitro and whole-organism models. We review five molecular systems of stress response and update 18 consensual cell-cell signaling pathways that are the hallmark for early development, organogenesis, and differentiation; and revisit the principles of teratology in light of recent advances in high-throughput screening, big data techniques, and systems toxicology. Multiscale systems modeling plays an integral role in the evolutionary approach to cross-species extrapolation. Phylogenetic analysis and comparative bioinformatics are both valuable tools in identifying and validating the molecular initiating events that account for adverse developmental outcomes in humans. The discordance of susceptibility between test species and humans (ontogeny) reflects their differences in evolutionary history (phylogeny). This synthesis not only can lead to novel applications in developmental toxicity and risk assessment, but also can pave the way for applying an evo-devo perspective to the study of developmental origins of health and disease.

Hepatic transcriptional responses to copper in the three-spined stickleback are affected by their pollution exposure history

Some fish populations inhabiting contaminated environments show evidence of increased chemical tolerance, however the mechanisms contributing to this tolerance, and whether this is heritable, are poorly understood. We investigated the responses of two populations of wild three-spined stickleback (Gasterosteus aculeatus) with different histories of contaminant exposure to an oestrogen and copper, two widespread aquatic pollutants. Male stickleback originating from two sites, the River Aire, with a history of complex pollution discharges, and Siblyback Lake, with a history of metal contamination, were depurated and then exposed to copper (46μg/L) and the synthetic oestrogen ethinyloestradiol (22ng/L). The hepatic transcriptomic response was compared between the two populations and to a reference population with no known history of exposure (Houghton Springs, Dorset). Gene responses included those typical for both copper and oestrogen, with no discernable difference in response to oestrogen between populations. There was, however, some difference in the magnitude of response to copper between populations. Siblyback fish showed an elevated baseline transcription of genes encoding metallothioneins and a lower level of metallothionein induction following copper exposure, compared to those from the River Aire. Similarly, a further experiment with an F1 generation of Siblyback fish bred in the laboratory found evidence for elevated transcription of genes encoding metallothioneins in unexposed fish, together with an altered transcriptional response to 125μg/L copper, compared with F1 fish originating from the clean reference population exposed to the same copper concentration. These data suggest that the stickleback from Siblyback Lake have a differential response to copper, which is inherited by the F1 generation in laboratory conditions, and for which the underlying mechanism may include an elevation of baseline transcription of genes encoding metallothioneins. The genetic and/or epigenetic mechanisms contributing to this inherited alteration of metallothionein transcription have yet to be established.

Evolutionary toxicology in an omics world

Evolutionary toxicology is a young field that has grown rapidly in the past two decades. The potential of this field comes from the ability to link chemical contamination to multigenerational and population-wide effects in various species. The advancements and rapidly decreasing costs of -omic tools are improving the power and resolution of evolutionary toxicology studies. In this manuscript, we aim to address the trajectories and perspectives for conducting evolutionary toxicology studies with -omic approaches. We discuss the complementarity of using multiple -omic tools (genomics, eDNA, transcriptomics, proteomics, and metabolomics) for utility in understanding the toxicological relevance of adaptive responses in populations. In addition, we discuss phenotypic plasticity and its relevance to transcriptomic studies in toxicology. As evolutionary toxicology grows and expands its capacity to link toxicology with population-wide end points, we emphasize the applications of such studies in answering questions about ecological and population health, as well as future applicability to regulation. Thus, we aim to emphasize the enormous potential for evolutionary toxicology in an -omics world and give perspectives on the directions of future investigations.

The genomic landscape of rapid repeated evolutionary adaptation to toxic pollution in wild fish

Atlantic killifish populations have rapidly adapted to normally lethal levels of pollution in four urban estuaries. Through analysis of 384 whole killifish genome sequences and comparative transcriptomics in four pairs of sensitive and tolerant populations, we identify the aryl hydrocarbon receptor-based signaling pathway as a shared target of selection. This suggests evolutionary constraint on adaptive solutions to complex toxicant mixtures at each site. However, distinct molecular variants apparently contribute to adaptive pathway modification among tolerant populations. Selection also targets other toxicity-mediating genes and genes of connected signaling pathways; this indicates complex tolerance phenotypes and potentially compensatory adaptations. Molecular changes are consistent with selection on standing genetic variation. In killifish, high nucleotide diversity has likely been a crucial substrate for selective sweeps to propel rapid adaptation.