Signatures of selection in natural populations adapted to chronic pollution

Transcriptomic analysis provides insights into molecular mechanisms of thermal physiology

Physiological trait variation underlies health, responses to global climate change, and ecological performance. Yet, most physiological traits are complex, and we have little understanding of the genes and genomic architectures that define their variation. To provide insight into the genetic architecture of physiological processes, we related physiological traits to heart and brain mRNA expression using a weighted gene co-expression network analysis. mRNA expression was used to explain variation in six physiological traits (whole animal metabolism (WAM), critical thermal maximum (CT_max), and four substrate specific cardiac metabolic rates (CaM)) under 12 °C and 28 °C acclimation conditions. Notably, the physiological trait variations among the three geographically close (within 15 km) and genetically similar F. heteroclitus populations are similar to those found among 77 aquatic species spanning 15–20° of latitude (~ 2,000 km). These large physiological trait variations among genetically similar individuals provide a powerful approach to determine the relationship between mRNA expression and heritable fitness related traits unconfounded by interspecific differences. Expression patterns explained up to 82% of metabolic trait variation and were enriched for multiple signaling pathways known to impact metabolic and thermal tolerance (e.g., AMPK, PPAR, mTOR, FoxO, and MAPK) but also contained several unexpected pathways (e.g., apoptosis, cellular senescence), suggesting that physiological trait variation is affected by many diverse genes.

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.

Interindividual plasticity in metabolic and thermal tolerance traits from populations subjected to recent anthropogenic heating

To better understand temperature's role in the interaction between local evolutionary adaptation and physiological plasticity, we investigated acclimation effects on metabolic performance and thermal tolerance among natural Fundulus heteroclitus (small estuarine fish) populations from different thermal environments. Fundulus heteroclitus populations experience large daily and seasonal temperature variations, as well as local mean temperature differences across their large geographical cline. In this study, we use three populations: one locally heated (32°C) by thermal effluence (TE) from the Oyster Creek Nuclear Generating Station, NJ, and two nearby reference populations that do not experience local heating (28°C). After acclimation to 12 or 28°C, we quantified whole-animal metabolic (WAM) rate, critical thermal maximum (CTmax) and substrate-specific cardiac metabolic rate (CaM, substrates: glucose, fatty acids, lactate plus ketones plus ethanol, and endogenous (i.e. no added substrates)) in approximately 160 individuals from these three populations. Populations showed few significant differences due to large interindividual variation within populations. In general, for WAM and CTmax, the interindividual variation in acclimation response (log2 ratio 28/12°C) was a function of performance at 12°C and order of acclimation (12-28°C versus 28-12°C). CTmax and WAM were greater at 28°C than 12°C, although WAM had a small change (2.32-fold) compared with the expectation for a 16°C increase in temperature (expect 3- to 4.4-fold). By contrast, for CaM, the rates when acclimatized and assayed at 12 or 28°C were nearly identical. The small differences in CaM between 12 and 28°C temperature were partially explained by cardiac remodeling where individuals acclimatized to 12°C had larger hearts than individuals acclimatized to 28°C. Correlation among physiological traits was dependent on acclimation temperature. For example, WAM was negatively correlated with CTmax at 12°C but positively correlated at 28°C. Additionally, glucose substrate supported higher CaM than fatty acid, and fatty acid supported higher CaM than lactate, ketones and alcohol (LKA) or endogenous. However, these responses were highly variable with some individuals using much more FA than glucose. These findings suggest interindividual variation in physiological responses to temperature acclimation and indicate that additional research investigating interindividual may be relevant for global climate change responses in many species.

On biological evolution and environmental solutions

Drawing insights from multiple disciplines is essential for finding integrative solutions that are required to tackle complex environmental problems. Human activities are causing unprecedented influence on global ecosystems, culminating in the loss of species and fundamental changes in the selective environments of organisms across the tree of life. Our collective understanding about biological evolution can help identify and mitigate many of the environmental problems in the Anthropocene. To this end, we propose a stronger integration of environmental sciences with evolutionary biology.

Shared genomic outliers across two divergent population clusters of a highly threatened seagrass

The seagrass, Zostera capensis, occurs across a broad stretch of coastline and wide environmental gradients in estuaries and sheltered bays in southern and eastern Africa. Throughout its distribution, habitats are highly threatened and poorly protected, increasing the urgency of assessing the genomic variability of this keystone species. A pooled genomic approach was employed to obtain SNP data and examine neutral genomic variation and to identify potential outlier loci to assess differentiation across 12 populations across the ∼9,600 km distribution of Z. capensis. Results indicate high clonality and low genomic diversity within meadows, which combined with poor protection throughout its range, increases the vulnerability of this seagrass to further declines or local extinction. Shared variation at outlier loci potentially indicates local adaptation to temperature and precipitation gradients, with Isolation-by-Environment significantly contributing towards shaping spatial variation in Z. capensis. Our results indicate the presence of two population clusters, broadly corresponding to populations on the west and east coasts, with the two lineages shaped only by frequency differences of outlier loci. Notably, ensemble modelling of suitable seagrass habitat provides evidence that the clusters are linked to historical climate refugia around the Last Glacial Maxi-mum. Our work suggests a complex evolutionary history of Z. capensis in southern and eastern Africa that will require more effective protection in order to safeguard this important ecosystem engineer into the future.

Evolutionary Toxicogenomics of the Striped Killifish (Fundulus majalis) in the New Bedford Harbor (Massachusetts, USA)

In this paper, we used a Genotyping-by-Sequencing (GBS) approach to find and genotype more than 4000 genome-wide SNPs (Single Nucleotide Polymorphisms) from striped killifish exposed to a variety of polychlorinated biphenyls (PCBs) and other aromatic pollutants in New Bedford Harbor (NBH, Massachusetts, USA). The aims of this study were to identify the genetic consequences of exposure to aquatic pollutants and detect genes that may be under selection. Low genetic diversity (H_E and π) was found in the site exposed to the highest pollution level, but the pattern of genetic diversity did not match the pollution levels. Extensive connectivity was detected among sampling sites, which suggests that balanced gene flow may explain the lack of genetic variation in response to pollution levels. Tests for selection identified 539 candidate outliers, but many of the candidate outliers were not shared among tests. Differences among test results likely reflect different test assumptions and the complex pollutant mixture. Potentially, selectively important loci are associated with 151 SNPs, and enrichment analysis suggests a likely involvement of these genes with pollutants that occur in NBH. This result suggests that selective processes at genes targeted by pollutants may be occurring, even at a small geographical scale, and may allow the local striped killifish to resist the high pollution levels.

Population genomics of rapid evolution in natural populations: polygenic selection in response to power station thermal effluents

BACKGROUND

Examples of rapid evolution are common in nature but difficult to account for with the standard population genetic model of adaptation. Instead, selection from the standing genetic variation permits rapid adaptation via soft sweeps or polygenic adaptation. Empirical evidence of this process in nature is currently limited but accumulating.

RESULTS

We provide genome-wide analyses of rapid evolution in Fundulus heteroclitus populations subjected to recently elevated temperatures due to coastal power station thermal effluents using 5449 SNPs across two effluent-affected and four reference populations. Bayesian and multivariate analyses of population genomic structure reveal a substantial portion of genetic variation that is most parsimoniously explained by selection at the site of thermal effluents. An FST outlier approach in conjunction with additional conservative requirements identify significant allele frequency differentiation that exceeds neutral expectations among exposed and closely related reference populations. Genomic variation patterns near these candidate loci reveal that individuals living near thermal effluents have rapidly evolved from the standing genetic variation through small allele frequency changes at many loci in a pattern consistent with polygenic selection on the standing genetic variation.

CONCLUSIONS

While the ultimate trajectory of selection in these populations is unknown and we survey only a minority of genomic loci, our findings suggest that polygenic models of adaptation may play important roles in large, natural populations experiencing recent selection due to environmental changes that cause broad physiological impacts.

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.

Toxicology tailored low density oligonucleotide microarray for the thicklip grey mullets (Chelon labrosus): Biomarker gene transcription profile after caging in a polluted harbour

In aquatic organisms inhabiting polluted waters genes are activated to build an adaptive/compensatory defence against the possible effects of pollutants. Such responses can be used as biomarkers of exposure to chemical compounds, outlining the molecular mechanisms activated under specific pollution scenarios. With the aim of exploiting such approach in environmental health assessment, toxicologically relevant gene fragments were sequenced in the thicklip grey mullet (Chelon labrosus) and a toxicologically tailored low-density (160 genes) oligonucleotide microarray was customised. The tool was validated comparing organ/sex specific gene expression profiles and characterising responses under laboratory exposure to model chemicals. Finally, juvenile mullets were caged in a polluted harbour and hepatic gene expression profiles analysed after 5 and 21 days of deployment. Cages were deployed in the inner (IH) and outer (OH) Pasaia harbour, Bay of Biscay. Mussels (Mytilus galloprovincialis) were also caged as biological matrix for chemical bioaccumulation analysis and stress biomarkers measurements. Slightly higher concentrations of chemicals (metals, tributyltin, PAHs, phthalates) were quantified in IH than in OH, fish bile metabolites also revealing higher availability of PAHs in IH. Lysosome membrane stability in mussels was reduced, indicating stress condition in both sites. The developed microarray discriminated mullets showing distinctive expression profiles depending on site and deployment time. Genes related to immune and hypoxia responses were regulated comparing IH and OH at day 5. Phase I and II biotransformation genes, such as cyp2, cyp3 and ugt, were up-regulated in IH, together with the aryl hydrocarbon receptor 2 (ahr2) and the ahr repressor. Similarly, TBT-binding proteins and genes involved in lipid metabolism (pparγ, cyp7) were up-regulated with deployment time. Even if nowadays higher throughput approaches for gene expression analyses are available, the developed mullet tool constitutes a comprehensive tool to assess molecular responses of mullets exposed to pollutants, although it remains to be explored whether it can be applied to assess pollutant exposure in active pollution monitorings and in environmental health assessment.

The role of contamination history and gender on the genotoxic responses of the crayfish Procambarus clarkii to a penoxsulam-based herbicide

The responses of non-target organisms to pesticide exposure are still poorly explored in what concerns the development of adjustments favouring population success. Owing to the vital role of DNA integrity, it is important to identify genome-maintenance skills and their determinant factors. Thus, the major aims of the present study were: (i) to assess the genotoxicity of the penoxsulam-based herbicide (Viper^®) to the crayfish Procambarus clarkii; (ii) to understand the influence of gender and contamination history in the genotoxic responses following exposure to this herbicide; (iii) to investigate the damage mechanisms involved in putative adjustments shown by P. clarkii. Two populations were tested, one from a reference site and the other from a historically contaminated site. Specimens from both populations were exposed to Viper^®, considering environmentally relevant penoxsulam concentrations (20 and 40 µg L^-1) and to a model genotoxicant (EMS). Comet assay was adopted to assess the genetic damage in gills. The results disclosed the genotoxicity of the herbicide to crayfish (a non-target organism). Additionally, organisms exposed to the highest concentration of penoxsulam signalized the influence of factor "population" towards the genotoxic pressure (measured as effective DNA breaks): P2 males from the historically impacted population displayed a significantly higher susceptibly (by up to 53.98%) when compared to control, while the homologous group from the reference population presented levels similar to its respective control. When DNA lesion-repair enzymes were considered, DNA oxidation patterns suggested an increased ability of this gender (39.75% lower than negative control) to deal with this particular type of damage, namely considering pyrimidines oxidation. It is worth remarking that the influence of the exposure history on the protection/vulnerability to the penoxsulam-based herbicide was only evident in males, despite depending on the type of DNA damage: when the non-specific damage was considered, organisms from the impacted population seemed to be more vulnerable while regarding to the oxidative damage, males from the impacted population appeared to be more protected than organisms that have never been exposed to penoxsulam. Overall, the influence of factors "gender" and "contamination history" was demonstrated as well as its dependence on DNA damage type was evident. EMS groups did not present the differences between populations, reinforcing the agent-specific adjustment hypothesis.These findings highlighted the importance of considering differential physiological backgrounds in ecogenotoxicological analysis, hence favouring the elaboration of more plausible and holistic approaches integrating the environmental risk assessment of pesticides.

Genome-wide scan reveals signatures of selection related to pollution adaptation in non-model estuarine Atlantic killifish (Fundulus heteroclitus)

In many human-altered ecosystems, organisms are increasingly faced with more diverse and complex environmental stressors and pollutant mixtures, to which the adaptations necessary to survive exposure are likely to be numerous and varied. Improving our understanding of the molecular mechanisms that underlie complex polygenic adaptations in natural settings requires significant toxicological, biochemical, physiological, and genomic data rarely available for non-model organisms. Here, we build upon two decades of study of adaptation to anthropogenic pollutants in a population of Atlantic killifish (Fundulus heteroclitus) that inhabits the creosote-contaminated Atlantic Wood Industries Superfund (AW) site on the Elizabeth River, Virginia in the United States. To better understand the genotypes that underlie previously characterized resistance to PCBs and PAHs, we performed Restriction site-Associated DNA sequencing (RADseq) on killifish from AW and two relatively clean reference sites (King's Creek-KC, and Mains Creek-MC). Across the genome, we analyzed over 83,000 loci and 12,000 single nucleotide polymorphisms (SNPs). Shared across both comparisons of killifish from polluted (AW) and relatively unpolluted (KC and MC) sites, we found eight genomic regions with smoothed FST values significantly (p < 0.001) elevated above background. Using the recently published F. heteroclitus reference genome, we identified candidate genes in these significant regions involved in the AHR pathway (e.g. AIP, ARNT1c), as well as genes relating to cardiac structure and function. These genes represent both previously characterized and potentially novel molecular adaptations involved with various aspects of resistance to these environmental toxins.

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.

Genetic and epigenetic variation in Spartina alterniflora following the Deepwater Horizon oil spill

Catastrophic events offer unique opportunities to study rapid population response to stress in natural settings. In concert with genetic variation, epigenetic mechanisms may allow populations to persist through severe environmental challenges. In 2010, the Deepwater Horizon oil spill devastated large portions of the coastline along the Gulf of Mexico. However, the foundational salt marsh grass, Spartina alterniflora, showed high resilience to this strong environmental disturbance. Following the spill, we simultaneously examined the genetic and epigenetic structure of recovering populations of S. alterniflora to oil exposure. We quantified genetic and DNA methylation variation using amplified fragment length polymorphism and methylation sensitive fragment length polymorphism (MS‐AFLP) to test the hypothesis that response to oil exposure in S. alterniflora resulted in genetically and epigenetically based population differentiation. We found high genetic and epigenetic variation within and among sites and found significant genetic differentiation between contaminated and uncontaminated sites, which may reflect nonrandom mortality in response to oil exposure. Additionally, despite a lack of genomewide patterns in DNA methylation between contaminated and uncontaminated sites, we found five MS‐AFLP loci (12% of polymorphic MS‐AFLP loci) that were correlated with oil exposure. Overall, our findings support genetically based differentiation correlated with exposure to the oil spill in this system, but also suggest a potential role for epigenetic mechanisms in population differentiation.

On the maintenance of genetic variation and adaptation to environmental change: considerations from population genomics in fishes

The first goal of this paper was to overview modern approaches to local adaptation, with a focus on the use of population genomics data to detect signals of natural selection in fishes. Several mechanisms are discussed that may enhance the maintenance of genetic variation and evolutionary potential, which have been overlooked and should be considered in future theoretical development and predictive models: the prevalence of soft sweeps, polygenic basis of adaptation, balancing selection and transient polymorphisms, parallel evolution, as well as epigenetic variation. Research on fish population genomics has provided ample evidence for local adaptation at the genome level. Pervasive adaptive evolution, however, seems to almost never involve the fixation of beneficial alleles. Instead, adaptation apparently proceeds most commonly by soft sweeps entailing shifts in frequencies of alleles being shared between differentially adapted populations. One obvious factor contributing to the maintenance of standing genetic variation in the face of selective pressures is that adaptive phenotypic traits are most often highly polygenic, and consequently the response to selection should derive mostly from allelic co-variances among causative loci rather than pronounced allele frequency changes. Balancing selection in its various forms may also play an important role in maintaining adaptive genetic variation and the evolutionary potential of species to cope with environmental change. A large body of literature on fishes also shows that repeated evolution of adaptive phenotypes is a ubiquitous evolutionary phenomenon that seems to occur most often via different genetic solutions, further adding to the potential options of species to cope with a changing environment. Moreover, a paradox is emerging from recent fish studies whereby populations of highly reduced effective population sizes and impoverished genetic diversity can apparently retain their adaptive potential in some circumstances. Although more empirical support is needed, several recent studies suggest that epigenetic variation could account for this apparent paradox. Therefore, epigenetic variation should be fully integrated with considerations pertaining to role of soft sweeps, polygenic and balancing selection, as well as repeated adaptation involving different genetic basis towards improving models predicting the evolutionary potential of species to cope with a changing world.

Sewage treatment plant associated genetic differentiation in the blue mussel from the Baltic Sea and Swedish west coast

Human-derived environmental pollutants and nutrients that reach the aquatic environment through sewage effluents, agricultural and industrial processes are constantly contributing to environmental changes that serve as drivers for adaptive responses and evolutionary changes in many taxa. In this study, we examined how two types of point sources of aquatic environmental pollution, harbors and sewage treatment plants, affect gene diversity and genetic differentiation in the blue mussel in the Baltic Sea area and off the Swedish west coast (Skagerrak). Reference sites (REF) were geographically paired with sites from sewage treatments plant (STP) and harbors (HAR) with a nested sampling scheme, and genetic differentiation was evaluated using a high-resolution marker amplified fragment length polymorphism (AFLP). This study showed that genetic composition in the Baltic Sea blue mussel was associated with exposure to sewage treatment plant effluents. In addition, mussel populations from harbors were genetically divergent, in contrast to the sewage treatment plant populations, suggesting that there is an effect of pollution from harbors but that the direction is divergent and site specific, while the pollution effect from sewage treatment plants on the genetic composition of blue mussel populations acts in the same direction in the investigated sites.

A Cost-Effective Approach to Sequence Hundreds of Complete Mitochondrial Genomes

We present a cost-effective approach to sequence whole mitochondrial genomes for hundreds of individuals. Our approach uses small reaction volumes and unmodified (non-phosphorylated) barcoded adaptors to minimize reagent costs. We demonstrate our approach by sequencing 383 Fundulus sp. mitochondrial genomes (192 F. heteroclitus and 191 F. majalis). Prior to sequencing, we amplified the mitochondrial genomes using 4–5 custom-made, overlapping primer pairs, and sequencing was performed on an Illumina HiSeq 2500 platform. After removing low quality and short sequences, 2.9 million and 2.8 million reads were generated for F. heteroclitus and F. majalis respectively. Individual genomes were assembled for each species by mapping barcoded reads to a reference genome. For F. majalis, the reference genome was built de novo. On average, individual consensus sequences had high coverage: 61-fold for F. heteroclitus and 57-fold for F. majalis. The approach discussed in this paper is optimized for sequencing mitochondrial genomes on an Illumina platform. However, with the proper modifications, this approach could be easily applied to other small genomes and sequencing platforms.

Cross-resistance in Gulf killifish (Fundulus grandis) populations resistant to dioxin-like compounds

The Houston Ship Channel (HSC) in Houston, Texas is an aquatic environment with a long history of contamination, including polychlorinated dibenzodioxins (PCDD), polychlorinated dibenzofurans (PCDF), polychlorinated biphenyls (PCBs), polycyclic aromatic hydrocarbons (PAHs), and heavy metals. Populations of Gulf killifish (Fundulus grandis) from the HSC have adapted to resist developmental cardiac deformities caused by dioxin-like compounds (DLCs). Contaminants in the HSC have acted as a strong selective pressure on resident Gulf killifish populations. Rapid adaptation can lead to fitness costs, some as a direct result of the mechanisms involved in the adaptive process, whereas other adaptations may be more general. To explore potential fitness costs, we evaluated two Gulf killifish populations with documented resistance to DLC-induced cardiac teratogenesis (Patrick Bayou and Vince Bayou), and one previously characterized reference population (Gangs Bayou). We also characterized a previously unstudied population from Galveston Bay as an additional reference population (Smith Point). We tested the sensitivity of F1 larvae from these four populations to two classes of pesticides (pyrethroid (permethrin) and carbamate (carbaryl)) and two model pro-oxidants (tert-butyl hydroquinone (tBHQ) and tert-butyl hydroperoxide (tBOOH)). In addition, we explored their responses to hypoxia and measured resting metabolic rates (M.O2). Both adapted populations were cross-resistant to the toxicity of carbaryl and both pro-oxidants tested. There were no population differences in sensitivity to permethrin. On the other hand, one reference population (Gangs Bayou) was less sensitive to hypoxia, and maintained a lower M.O2 . However, there were no differences in hypoxia tolerance or resting metabolic rate between the second reference and the two adapted populations. This investigation emphasizes the importance of including multiple reference populations to clearly link fitness costs or cross-resistance to pollution adaptation, rather than to unrelated environmental or ecological differences. When compared to previous literature on adapted populations of Fundulus heteroclitus, we see a mixture of similarities and differences, suggesting that F. grandis adapted phenotypes likely involve multiple mechanisms, which may not be completely consistent among adapted populations.

Disruption of the stress response in wastewater treatment works effluent-exposed three-spined sticklebacks persists after translocation to an unpolluted environment

The hypothalamic-pituitary-adrenal/interrenal (HPA/I) axis plays a key role in responding to biotic and abiotic challenges in all vertebrates. Recent studies have shown that the apical response of the HPI axis to stressors in three-spined sticklebacks varies in proportion to the concentration of wastewater treatment works (WWTW) effluent to which the fish are exposed. This study was conducted to determine whether between-site variation in stress responsiveness among WWTW effluent-exposed sticklebacks is persistent or reversible. Sticklebacks from eight sites in north-west England affected by WWTW effluent and exhibiting between-population variation in HPI axis reactivity, were moved to a clean-water aquarium environment. After five months in the contaminant-free environment the responsiveness of these fish to a standardised stressor was determined, by measuring the rate of stress-induced cortisol release across the gills, and compared with the responses of fish newly sampled from the eight original capture sites. Inter-site differences in the reactivity of the HPI axis, proportional to the effluent concentration at each site, persisted among the translocated female sticklebacks for at least 5 months. In male fish however, the direct relationship between stress responsiveness and site-specific effluent was not evident 5 months post-translocation. These results support previous observations that the HPA/I axis, a non-reproductive endocrine system, is vulnerable to modulation by anthropogenic factors in fish and show for the first time that, in female fish at least, this modulation is not transient. The mechanisms underlying these observations, and the implications for the fitness and resilience of affected populations, requires investigation.

Evaluating the impact of a fluoropolymer plant on a river macrobenthic community by a combined chemical, ecological and genetic approach

Effect-based monitoring is a recommended approach suggested in European Guidelines to assess the response of ecosystem affected by a pollution source, considering the effects at community, population, individual but also at suborganism level. A combined chemical, ecological and genetic approach was applied in order to assess the impact of a fluoropolymer plant on the macrobenthic community of the Northern Italian river Bormida (Piedmont region). The macrobenthic community living downstream of the industrial discharge was chronically exposed to a mixture of perfluoroalkyl substances (PFAS), with perfluorooctanoic acid as the main compound, at concentrations up to several μgL(-1). Ecological assessment proved that the downstream community was not substantially different from that living upstream of the pollution source. The impact on community is not quantifiable with the traditional monitoring methods used for ecological classification under European regulation because macrobenthic communities showed only slight differences in their structure. In order to highlight effects on genetic variability of the native population, a subcellular analysis by using the AFLP (Amplified Fragment Length Polymorphism) genetic technique was applied to genotype of individuals of a selected species (Hydropsyche modesta, Trichoptera) collected in the two sampling sites. Percentage of variation between the two populations was 6.8%, a threshold compatible with a genetic drift induced in the downstream population. The genetic study carried out in field identified a significant divergence between exposed and non-exposed populations, but at present it is not possible to associate this divergence to a specific effect induced by PFAS.

Do stressful conditions make adaptation difficult? Guppies in the oil-polluted environments of southern Trinidad

The ability of populations to rapidly adapt to new environments will determine their future in an increasingly human-modified world. Although meta-analyses do frequently uncover signatures of local adaptation, they also reveal many exceptions. We suggest that particular constraints on local adaptation might arise when organisms are exposed to novel stressors, such as anthropogenic pollution. To inform this possibility, we studied the extent to which guppies (Poecilia reticulata) show local adaptation to oil pollution in southern Trinidad. Neutral genetic markers revealed that paired populations in oil-polluted versus not-polluted habitats diverged independently in two different watersheds. Morphometrics revealed some divergence (particularly in head shape) between these environments, some of which was parallel between rivers. Reciprocal transplant experiments in nature, however, found little evidence of local adaptation based on survival and growth. Moreover, subsequent laboratory experiments showed that the two populations from oil-polluted sites showed only weak local adaptation even when compared to guppies from oil-free northern Trinidad. We conclude that guppies show little local adaptation to oil pollution, which might result from the challenges associated with adaptation to particularly stressful environments. It might also reflect genetic drift owing to small population sizes and/or high gene flow between environments.

Effects of Anthropogenic Pollution on the Oxidative Phosphorylation Pathway of Hepatocytes from Natural Populations of Fundulus heteroclitus

Persistent organic pollutants (POPs), including polycyclic aromatic hydrocarbons (PAHs) and polychlorinated biphenyls (PCBs), potentially target mitochondria and cause toxicity. We compared the effects of POPs on mitochondrial respiration by measuring oxidative phosphorylation (OxPhos) metabolism in hepatocytes isolated from lab-depurated Fundulus heteroclitus from a Superfund site contaminated with PAHs (Elizabeth River VA, USA) relative to OxPhos metabolism in individuals from a relatively clean, reference population (King's Creek VA, USA). In individuals from the polluted Elizabeth River population, OxPhos metabolism displayed lower LEAK and lower activities in complex III, complex IV, and E State, but higher activity in complex I compared to individuals from the reference King's Creek population. To test the supposition that these differences were due to or related to the chronic PAH contamination history of the Elizabeth River population, we compared the OxPhos functions of undosed individuals from the polluted and reference populations to individuals from these populations dosed with a PAH {benzo [α] pyrene (BaP)} or a PCB {PCB126 (3,3',4,4',5-pentachlorobiphenyl)}, respectively. Exposure to PAH or PCB affected OxPhos in the reference King's Creek population but had no detectable effects on the polluted Elizabeth River population. Thus, PAH exposure significantly increased LEAK, and exposure to PCB126 significantly decreased State 3, E state and complex I activity in the reference King's Creek population. These data strongly implicate an evolved tolerance in the Elizabeth River fish where dosed fish are not affected by PAH exposure and undosed fish show decreased LEAK and increased State 3 and E state.

Expression and function of ryanodine receptor related pathways in PCB tolerant Atlantic killifish (Fundulus heteroclitus) from New Bedford Harbor, MA, USA

Atlantic killifish (Fundulus heteroclitus) thrive in New Bedford Harbor (NBH), MA, highly contaminated with polychlorinated biphenyls (PCBs). Resident killifish have evolved tolerance to dioxin-like (DL) PCBs, whose toxic effects through the aryl hydrocarbon receptor (AhR) are well studied. In NBH, non-dioxin like PCBs (NDL PCBs), which lack activity toward the AhR, vastly exceed levels of DL congeners yet how killifish counter NDL toxic effects has not been explored. In mammals and fish, NDL PCBs are potent activators of ryanodine receptors (RyR), Ca(2+) release channels necessary for a vast array of physiological processes. In the current study we compared the expression and function of RyR related pathways in NBH killifish with killifish from the reference site at Scorton Creek (SC, MA). Relative to the SC fish, adults from NBH displayed increased levels of skeletal muscle RyR1 protein, and increased levels of FK506-binding protein 12 kDa (FKBP12) an accessory protein essential for NDL PCB-triggered changes in RyR channel function. In accordance with increased RyR1 levels, NBH killifish displayed increased maximal ligand binding, increased maximal response to Ca(2+) activation and increased maximal response to activation by the NDL PCB congener PCB 95. Compared to SC, NBH embryos and larvae had increased levels of mtor and ryr2 transcripts at multiple stages of development, and generations, while levels of serca2 were decreased at 9 days post-fertilization in the F1 and F2 generations. These findings suggest that there are compensatory and heritable changes in RyR mediated Ca(2+) signaling proteins or potential signaling partners in NBH killifish.

Ten years of transcriptomics in wild populations: what have we learned about their ecology and evolution?

Molecular ecology has moved beyond the use of a relatively small number of markers, often noncoding, and it is now possible to use whole-genome measures of gene expression with microarrays and RNAseq (i.e. transcriptomics) to capture molecular response to environmental challenges. While transcriptome studies are shedding light on the mechanistic basis of traits as complex as personality or physiological response to catastrophic events, these approaches are still challenging because of the required technical expertise, difficulties with analysis and cost. Still, we found that in the last 10 years, 575 studies used microarrays or RNAseq in ecology. These studies broadly address three questions that reflect the progression of the field: (i) How much variation in gene expression is there and how is it structured? (ii) How do environmental stimuli affect gene expression? (iii) How does gene expression affect phenotype? We discuss technical aspects of RNAseq and microarray technology, and a framework that leverages the advantages of both. Further, we highlight future directions of research, particularly related to moving beyond correlation and the development of additional annotation resources. Measuring gene expression across an array of taxa in ecological settings promises to enrich our understanding of ecology and genome function.

The Elizabeth River Story: A Case Study in Evolutionary Toxicology

The Elizabeth River system is an estuary in southeastern Virginia, surrounded by the towns of Chesapeake, Norfolk, Portsmouth, and Virginia Beach. The river has played important roles in U.S. history and has been the location of various military and industrial activities. These activities have been the source of chemical contamination in this aquatic system. Important industries, until the 1990s, included wood treatment plants that used creosote, an oil-derived product that is rich in polycyclic aromatic hydrocarbons (PAH). These plants left a legacy of PAH pollution in the river, and in particular Atlantic Wood Industries is a designated Superfund site now undergoing remediation. Numerous studies examined the distribution of PAH in the river and impacts on resident fauna. This review focuses on how a small estuarine fish with a limited home range, Fundulus heteroclitus (Atlantic killifish or mummichog), has responded to this pollution. While in certain areas of the river this species has clearly been impacted, as evidenced by elevated rates of liver cancer, some subpopulations, notably the one associated with the Atlantic Wood Industries site, displayed a remarkable ability to resist the marked effects PAH have on the embryonic development of fish. This review provides evidence of how pollutants have acted as evolutionary agents, causing changes in ecosystems potentially lasting longer than the pollutants themselves. Mechanisms underlying this evolved resistance, as well as mechanisms underlying the effects of PAH on embryonic development, are also described. The review concludes with a description of ongoing and promising efforts to restore this historic American river.

Salinity-dependent copper accumulation in the guppy Poecilia vivipara is associated with CTR1 and ATP7B transcriptional regulation

Copper (Cu) accumulation and regulation of key-genes involved in Cu homeostasis were evaluated in freshwater- and saltwater-acclimated guppies Poecilia vivipara. Fish were exposed (96h) to environmentally relevant concentrations of dissolved Cu (0, 5.0, 9.0 and 20.0μg/L). In freshwater guppies, gill and liver Cu accumulation was dependent on Cu concentration in the exposure medium. In saltwater guppies, this dependence was observed only in the gut. These findings indicate that Cu accumulation was salinity- and tissue-dependent. Key genes involved in Cu metabolism were sequenced for the first time in P. vivipara. Transcripts coding for the high-affinity copper transporter (CTR1) and copper-transporting ATPase (ATP7B) were identified using polymerase chain reaction (PCR) and gene sequencing. The full-length CTR1 open reading frame (1560bp) and a partial ATP7B (690bp) were discovered. Predicted amino acid sequences shared high identities with the CTR1 of Fundulus heteroclitus (81%) and the ATP7B of Sparus aurata (87%). Basal transcriptional levels addressed by RT-qPCR in control fish indicate that CTR1 and ATP7B was highly transcribed in liver of freshwater guppies while CTR1 was highly transcribed in gut of saltwater guppies. This could explain the higher Cu accumulation observed in liver of freshwater guppies and in gut of saltwater guppies, because CTR1 is involved in Cu uptake. Reduced gill mRNA expression of CTR1 was observed in freshwater guppies exposed to 20.0μg/L Cu and in saltwater guppies exposed to 5.0μg/L Cu. In turn, reduced mRNA expression of gut ATP7B was observed in freshwater and salt water guppies exposed to 9.0 and 20.0μg/L Cu. Liver CTR1 and ATP7B transcription were not affected by Cu exposure. These findings suggest that gill CTR1 and gut ATP7B are down-regulated to limit Cu absorption after exposure to dissolved Cu, while liver CTR1 and ATP7B levels are maintained to allow Cu storage and detoxification. In conclusion, findings reported here indicate that Cu accumulation in the euryhaline guppy P. vivipara is tissue specific and dependent on water salinity. They also suggest that Cu homeostasis involves a differential transcriptional regulation of the newly identified Cu transporters, CTR1 and ATP7B.

Targeted approach to identify genetic loci associated with evolved dioxin tolerance in Atlantic killifish (Fundulus heteroclitus)

Background

The most toxic aromatic hydrocarbon pollutants are categorized as dioxin-like compounds (DLCs) to which extreme tolerance has evolved independently and contemporaneously in (at least) four populations of Atlantic killifish (Fundulus heteroclitus). Surprisingly, the magnitude and phenotype of DLC tolerance is similar among these killifish populations that have adapted to varied, but highly aromatic hydrocarbon-contaminated urban/industrialized estuaries of the US Atlantic coast. Multiple tolerant and neighboring sensitive killifish populations were compared with the expectation that genetic loci associated with DLC tolerance would be revealed.

Results

Since the aryl hydrocarbon receptor (AHR) pathway partly or fully mediates DLC toxicity in vertebrates, single nucleotide polymorphisms (SNPs) from 42 genes associated with the AHR pathway were identified to serve as targeted markers. Wild fish (N = 36/37) from four highly tolerant killifish populations and four nearby sensitive populations were genotyped using 59 SNP markers. Similar to other killifish population genetic analyses, strong genetic differentiation among populations was detected, consistent with isolation by distance models. When DLC-sensitive populations were pooled and compared to pooled DLC-tolerant populations, multi-locus analyses did not distinguish the two groups. However, pairwise comparisons of nearby tolerant and sensitive populations revealed high differentiation among sensitive and tolerant populations at these specific loci: AHR 1 and 2, cathepsin Z, the cytochrome P450s (CYP1A and 3A30), and the NADH dehydrogenase subunits. In addition, significant shifts in minor allele frequency were observed at AHR2 and CYP1A loci across most sensitive/tolerant pairs, but only AHR2 exhibited shifts in the same direction across all pairs.

Conclusions

The observed differences in allelic composition at the AHR2 and CYP1A SNP loci were identified as significant among paired sensitive/tolerant populations of Atlantic killifish with multiple statistical tests. The genetic patterns reported here lend support to the argument that AHR2 and CYP1A play a role in the adaptive response to extreme DLC contamination. Additional functional assays are required to isolate the exact mechanism of DLC tolerance.

Genomic and physiological responses to strong selective pressure during late organogenesis: few gene expression changes found despite striking morphological differences

Background

Adaptations to a new environment, such as a polluted one, often involve large modifications of the existing phenotypes. Changes in gene expression and regulation during critical developmental stages may explain these phenotypic changes. Embryos from a population of the teleost fish, Fundulus heteroclitus, inhabiting a clean estuary do not survive when exposed to sediment extract from a site highly contaminated with polycyclic aromatic hydrocarbons (PAHs) while embryos derived from a population inhabiting a PAH polluted estuary are remarkably resistant to the polluted sediment extract. We exposed embryos from these two populations to surrogate model PAHs and analyzed changes in gene expression, morphology, and cardiac physiology in order to better understand sensitivity and adaptive resistance mechanisms mediating PAH exposure during development.

Results

The synergistic effects of two model PAHs, an aryl hydrocarbon receptor (AHR) agonist (β-naphthoflavone) and a cytochrome P4501A (CYP1A) inhibitor (α-naphthoflavone), caused significant developmental delays, impaired cardiac function, severe morphological alterations and failure to hatch, leading to the deaths of reference embryos; resistant embryos were mostly unaffected. Unexpectedly, patterns of gene expression among normal and moderately deformed embryos were similar, and only severely deformed embryos showed a contrasting pattern of gene expression. Given the drastic morphological differences between reference and resistant embryos, a surprisingly low percentage of genes, 2.24% of 6,754 analyzed, show statistically significant differences in transcript levels during late organogenesis between the two embryo populations.

Conclusions

Our study demonstrates important contrasts in responses between reference and resistant natural embryo populations to synergistic effects of surrogate model PAHs that may be important in adaptive mechanisms mediating PAH effects during fish embryo development. These results suggest that statistically significant changes in gene expression of relatively few genes contribute to the phenotypic changes and large morphological differences exhibited by reference and resistant populations upon exposure to PAH pollutants. By correlating cardiac physiology and morphology with changes in gene expression patterns of reference and resistant embryos, we provide additional evidence for acquired resistance among embryos whose parents live at heavily contaminated sites.

Population bottlenecks, genetic diversity and breeding ability of the three-spined stickleback (Gasterosteus aculeatus) from three polluted English Rivers

Pollution is a significant environmental pressure on fish populations in both freshwater and marine environments. Populations subjected to chronic exposure to pollutants can experience impacts ranging from altered reproductive capacity to changes in population genetic structure. Few studies, however, have examined the reproductive vigor of individuals within populations inhabiting environments characterized by chronic pollution. In this study we undertook an analysis of populations of three-spined sticklebacks (Gasterosteus aculeatus) from polluted sites, to determine levels of genetic diversity, assess for evidence of historic population genetic bottlenecks and determine the reproductive competitiveness of males from these locations. The sites chosen included locations in the River Aire, the River Tees and the River Birket, English rivers that have been impacted by pollution from industrial and/or domestic effluents for over 100 years. Male reproductive competitiveness was determined via competitive breeding experiments with males and females derived from a clean water site, employing DNA microsatellites to determine parentage outcome. Populations of stickleback collected from the three historically polluted sites showed evidence of recent population bottlenecks, although only the River Aire population showed low genetic diversity. In contrast, fish collected from two relatively unpolluted sites within the River Gowy and Houghton Springs showed weak, or no evidence of such bottlenecks. Nevertheless, males derived from polluted sites were able to reproduce successfully in competition with males derived from clean water exposures, indicating that these bottlenecks have not resulted in any substantial loss of reproductive fitness in males.

Differential tolerance to copper, but no evidence of population-level genetic differences in a widely-dispersing native barnacle

Despite many estuaries having high levels of metal pollution, species are found to persist in these stressful environments. Populations of estuarine invertebrates exposed to toxic concentrations of such metals may be under selection. However, in species with a wide-dispersal potential, any short-term results of localized selection may be counteracted by external recruitment from populations not under selection. The barnacle Amphibalanus variegatus is found in nearshore coastal environments as well as sheltered embayments and estuaries, including metal-impacted estuaries, from New South Wales, Australia to Western Australia. The fertilised eggs of A. variegatus are brooded internally and released as larvae (nauplii), which remain in the water-column for ~14 days before settling. Hence the species has a considerable dispersal capacity. The purpose of this study was to examine whether populations of A. variegatus from metal-impacted sites, displayed a greater tolerance to a toxicant (copper) than reference populations. Adult barnacles where collected from twenty sites within two metal-impacted and fourteen sites within two reference estuaries. Within 24 h, adults were induced to spawn and the offspring were exposed to copper in a laboratory assay. Larvae collected from the metal-impacted estuaries demonstrated a greater tolerance to copper compared to those from reference sites. To determine if selection/localised in the metal impacted sites was occurring, the genetic structure of populations at three sites was examined using an AFLP methodology. No evidence of unique population identity and or selection (outlier loci) was detected suggesting that: (1) the tolerance displayed in the assay was derived from acclimation during development; and/or (2) that the populations are open preventing the fixation of any unique alleles.

Population genetics provides new insights into biomarker prevalence in dab (Limanda limanda L.): a key marine biomonitoring species

Bioindicators are species for which some quantifiable aspect of its biology, a biomarker, is assumed to be sensitive to ecosystem health. However, there is frequently a lack of information on the underlying genetic and environmental drivers shaping the spatiotemporal variance in prevalence of the biomarkers employed. Here, we explore the relative role of potential variables influencing the spatiotemporal prevalence of biomarkers in dab, Limanda limanda, a species used as a bioindicator of marine contaminants. Firstly, the spatiotemporal genetic structure of dab around UK waters (39 samples across 15 sites for four years: 2005–2008) is evaluated with 16 microsatellites. Two temporally stable groups are identified corresponding to the North and Irish Seas (average between basin = 0.007; = 0.022). Secondly, we examine the association between biomarker prevalence and several variables, including genetic structuring, age and contaminant exposure. Genetic structure had significant interactive effects, together with age and some contaminants, in the prevalence of some of the biomarkers considered, namely hyperpigmentation and liver lesions. The integration of these data sets enhanced our understanding of the relationship between biomarker prevalence, exposure to contaminants and population-specific response, thereby yielding more informative predictive models of response and prospects for environmental remediation.

Going where traditional markers have not gone before: utility of and promise for RAD sequencing in marine invertebrate phylogeography and population genomics

Characterization of large numbers of single-nucleotide polymorphisms (SNPs) throughout a genome has the power to refine the understanding of population demographic history and to identify genomic regions under selection in natural populations. To this end, population genomic approaches that harness the power of next-generation sequencing to understand the ecology and evolution of marine invertebrates represent a boon to test long-standing questions in marine biology and conservation. We employed restriction-site-associated DNA sequencing (RAD-seq) to identify SNPs in natural populations of the sea anemone Nematostella vectensis, an emerging cnidarian model with a broad geographic range in estuarine habitats in North and South America, and portions of England. We identified hundreds of SNP-containing tags in thousands of RAD loci from 30 barcoded individuals inhabiting four locations from Nova Scotia to South Carolina. Population genomic analyses using high-confidence SNPs resulted in a highly-resolved phylogeography, a result not achieved in previous studies using traditional markers. Plots of locus-specific FST against heterozygosity suggest that a majority of polymorphic sites are neutral, with a smaller proportion suggesting evidence for balancing selection. Loci inferred to be under balancing selection were mapped to the genome, where 90% were located in gene bodies, indicating potential targets of selection. The results from analyses with and without a reference genome supported similar conclusions, further highlighting RAD-seq as a method that can be efficiently applied to species lacking existing genomic resources. We discuss the utility of RAD-seq approaches in burgeoning Nematostella research as well as in other cnidarian species, particularly corals and jellyfishes, to determine phylogeographic relationships of populations and identify regions of the genome undergoing selection.

Interactive effects of inbreeding and endocrine disruption on reproduction in a model laboratory fish

Inbreeding depression is expected to be more severe in stressful environments. However, the extent to which inbreeding affects the vulnerability of populations to environmental stressors, such as chemical exposure, remains unresolved. Here we report on the combined impacts of inbreeding and exposure to an endocrine disrupting chemical (the fungicide clotrimazole) on zebrafish (Danio rerio). We show that whilst inbreeding can negatively affect reproductive traits, not all traits are affected equally. Inbreeding depression frequently only became apparent when fish were additionally stressed by chemical exposure. Embryo viability was significantly reduced in inbred exposed fish and there was a tendency for inbred males to sire fewer offspring when in direct competition with outbred individuals. Levels of plasma 11-ketotestosterone, a key male sex hormone, showed substantial inbreeding depression that was unaffected by addition of the fungicide. In contrast, there was no effect of inbreeding or clotrimazole exposure on egg production. Overall, our data provide evidence that stress may amplify the effects of inbreeding on key reproductive traits, particularly those associated with male fitness. This may have important implications when considering the consequences of exposure to chemical pollutants on the fitness of wild populations.

Outlier Loci and Selection Signatures of Simple Sequence Repeats (SSRs) in Flax (Linum usitatissimum L.)

Genomic microsatellites (gSSRs) and expressed sequence tag-derived SSRs (EST-SSRs) have gained wide application for elucidating genetic diversity and population structure in plants. Both marker systems are assumed to be selectively neutral when making demographic inferences, but this assumption is rarely tested. In this study, three neutrality tests were assessed for identifying outlier loci among 150 SSRs (85 gSSRs and 65 EST-SSRs) that likely influence estimates of population structure in three differentiated flax sub-populations (F_ST = 0.19). Moreover, the utility of gSSRs, EST-SSRs, and the combined sets of SSRs was also evaluated in assessing genetic diversity and population structure in flax. Six outlier loci were identified by at least two neutrality tests showing footprints of balancing selection. After removing the outlier loci, the STRUCTURE analysis and the dendrogram topology of EST-SSRs improved. Conversely, gSSRs and combined SSRs results did not change significantly, possibly as a consequence of the higher number of neutral loci assessed. Taken together, the genetic structure analyses established the superiority of gSSRs to determine the genetic relationships among flax accessions, although the combined SSRs produced the best results. Genetic diversity parameters did not differ statistically (P > 0.05) between gSSRs and EST-SSRs, an observation partially explained by the similar number of repeat motifs. Our study provides new insights into the ability of gSSRs and EST-SSRs to measure genetic diversity and structure in flax and confirms the importance of testing for the occurrence of outlier loci to properly assess natural and breeding populations, particularly in studies considering only few loci.

AFLP genome scan to detect genetic structure and candidate loci under selection for local adaptation of the invasive weed Mikania micrantha

Why some species become successful invaders is an important issue in invasive biology. However, limited genomic resources make it very difficult for identifying candidate genes involved in invasiveness. Mikania micrantha H.B.K. (Asteraceae), one of the world's most invasive weeds, has adapted rapidly in response to novel environments since its introduction to southern China. In its genome, we expect to find outlier loci under selection for local adaptation, critical to dissecting the molecular mechanisms of invasiveness. An explorative amplified fragment length polymorphism (AFLP) genome scan was used to detect candidate loci under selection in 28 M. micrantha populations across its entire introduced range in southern China. We also estimated population genetic parameters, bottleneck signatures, and linkage disequilibrium. In binary characters, such as presence or absence of AFLP bands, if all four character combinations are present, it is referred to as a character incompatibility. Since character incompatibility is deemed to be rare in populations with extensive asexual reproduction, a character incompatibility analysis was also performed in order to infer the predominant mating system in the introduced M. micrantha populations. Out of 483 AFLP loci examined using stringent significance criteria, 14 highly credible outlier loci were identified by Dfdist and Bayescan. Moreover, remarkable genetic variation, multiple introductions, substantial bottlenecks and character compatibility were found to occur in M. micrantha. Thus local adaptation at the genome level indeed exists in M. micrantha, and may represent a major evolutionary mechanism of successful invasion. Interactions between genetic diversity, multiple introductions, and reproductive modes contribute to increase the capacity of adaptive evolution.

Biomarkers of organic contamination in the South American fish Poecilia vivipara and Jenynsia multidentata

South American cyprinodontiform fish are potential candidates to be used as model biomarker species of exposure in environmental toxicology. The aim of this study was to identify molecular and biochemical biomarkers of pollution using Poecilia vivipara (Poecilidae) and Jenynsia multidentata (Anablepidae). Partial nucleotide sequences for cytochrome P-450 1A (CYP1A), a classical biomarker of exposure to organic contaminants in fish, were identified in P. vivipara and J. multidentata (approximately 650 nucleotides) using degenerated primers and polymerase chain reaction (PCR). These sequences shared approximately 90% identity in the predicted amino acid sequence with the corresponding CYP1A region of Fundulus heteroclitus. Real-time quantitative PCR (RT-qPCR) analysis confirmed that CYP1A transcription was markedly induced in the liver and gills of J. multidentata (approximately185-fold and 20-fold, respectively) and P. vivipara (122-fold and 739-fold, respectively) 24 h after exposure to 1 μM synthetic CYP1A inducer β-naphthoflavone (BNF). At 24 h after injection with 1 μg/g environmental carcinogenic contaminant benzo[a]pyrene (BaP), a decreased total antioxidant capacity against peroxyl radicals was observed both in liver of J. multidentata and gills of P. vivipara. BaP injection in both fish did not produce changes in lipid peroxide (thiobarbituric acid-reactive substances, TBARS) levels, suggesting an absence of an oxidative stress condition. The newly identified CYP1A may thus serve as general biomarker of exposure to organic contaminant in future studies using P. vivipara and J. multidentata. Data also indicate the importance of species-specific differences in biomarker responses in these South American cyprinodontiform fish, suggesting distinct resistance/susceptibility properties to polycyclic aromatic hydrocarbons.

AFLP genome scan to detect genetic structure and candidate loci under selection for local adaptation of the invasive weed Mikania micrantha

Why some species become successful invaders is an important issue in invasive biology. However, limited genomic resources make it very difficult for identifying candidate genes involved in invasiveness. Mikania micrantha H.B.K. (Asteraceae), one of the world's most invasive weeds, has adapted rapidly in response to novel environments since its introduction to southern China. In its genome, we expect to find outlier loci under selection for local adaptation, critical to dissecting the molecular mechanisms of invasiveness. An explorative amplified fragment length polymorphism (AFLP) genome scan was used to detect candidate loci under selection in 28 M. micrantha populations across its entire introduced range in southern China. We also estimated population genetic parameters, bottleneck signatures, and linkage disequilibrium. In binary characters, such as presence or absence of AFLP bands, if all four character combinations are present, it is referred to as a character incompatibility. Since character incompatibility is deemed to be rare in populations with extensive asexual reproduction, a character incompatibility analysis was also performed in order to infer the predominant mating system in the introduced M. micrantha populations. Out of 483 AFLP loci examined using stringent significance criteria, 14 highly credible outlier loci were identified by Dfdist and Bayescan. Moreover, remarkable genetic variation, multiple introductions, substantial bottlenecks and character compatibility were found to occur in M. micrantha. Thus local adaptation at the genome level indeed exists in M. micrantha, and may represent a major evolutionary mechanism of successful invasion. Interactions between genetic diversity, multiple introductions, and reproductive modes contribute to increase the capacity of adaptive evolution.

Evolutionary and functional analyses of cytochrome P4501A promoter polymorphisms in natural populations

The functional importance of variable, transcriptional regulatory sequences within and among natural populations is largely unexplored. We analysed the cytochrome P4501A (CYP1A) promoter in three populations of the minnow, Fundulus heteroclitus, because two SNPs in the promoter and first intron of CYP1A are under selection in populations adapted to pollutants. To define the importance of these SNPs, 1630 bp of the CYP1A promoter and first intron and exon were sequenced in eight individuals from three populations: a population from a polluted environment resistant to some aromatic pollutants and two flanking reference populations. CYP1A is induced by many aromatic pollutants, but in populations adapted to pollutants, CYP1A has been shown to be refractory to induction. We were interested in understanding whether variation in the CYP1A promoter explains mechanism(s) of adaptation to these aromatic pollutants. The CYP1A promoter was extremely variable (an average of 9.3% of the promoter nucleotides varied among all populations) and exhibited no fixed differences between populations. As CYP1A is poorly inducible in adapted fish, we hypothesized that CYP1A promoter regions might vary functionally between populations. Unexpectedly, in vitro analysis showed significantly greater transcription from CYP1A promoters found in the population from the polluted environment relative to promoters found in both reference populations. Thus, despite extensive variation among populations and lack of fixed differences between populations, individuals from a polluted environment have significantly enhanced promoter activity. These data demonstrate that intraspecific variation, which provides the raw material for natural selection to act on, can occur while maintaining promoter function.

Common mechanism underlies repeated evolution of extreme pollution tolerance

Human alterations to the environment can exert strong evolutionary pressures, yet contemporary adaptation to human-mediated stressors is rarely documented in wildlife populations. A common-garden experimental design was coupled with comparative transcriptomics to discover evolved mechanisms enabling three populations of killifish resident in urban estuaries to survive normally lethal pollution exposure during development, and to test whether mechanisms are unique or common across populations. We show that killifish populations from these polluted sites have independently converged on a common adaptive mechanism, despite variation in contaminant profiles among sites. These populations are united by a similarly profound desensitization of aryl-hydrocarbon receptor-mediated transcriptional activation, which is associated with extreme tolerance to the lethal effects of toxic dioxin-like pollutants. The rapid, repeated, heritable and convergent nature of evolved tolerance suggests that ancestral killifish populations harboured genotypes that enabled adaptation to twentieth-century industrial pollutants.

Applications and implications of neutral versus non-neutral markers in molecular ecology

The field of molecular ecology has expanded enormously in the past two decades, largely because of the growing ease with which neutral molecular genetic data can be obtained from virtually any taxonomic group. However, there is also a growing awareness that neutral molecular data can provide only partial insight into parameters such as genetic diversity, local adaptation, evolutionary potential, effective population size, and taxonomic designations. Here we review some of the applications of neutral versus adaptive markers in molecular ecology, discuss some of the advantages that can be obtained by supplementing studies of molecular ecology with data from non-neutral molecular markers, and summarize new methods that are enabling researchers to generate data from genes that are under selection.