Familiarity and sex modulate size-dependent following behaviour in the Mediterranean killifish

Living in a group offers the chance to follow the choices and the behaviours of other individuals. Following a group mate might confer fitness advantages if the group mate knows about resources such as food or shelters. Shoaling fish often follow larger group mates which, in most species, are generally older and therefore more experienced. Yet, the effect of individuals' characteristics other than size on following behaviour remains to be understood. For example, familiar fish and female shoals have been reported as more cohesive, which might be due to a differential tendency to follow in relation to familiarity and sex. Here, we investigated whether size, familiarity, sex, and the interaction between these factors affect fish following behaviour. We observed pairs of differently sized Mediterranean killifish, Aphanius fasciatus, exploring a new environment, and we recorded whether the rear fish followed the front fish when the latter changed swimming direction. In female and male pairs, and in unfamiliar pairs, smaller fish were more likely to follow the directional change of the larger fish than vice versa. In mixed-sex pairs and in familiar pairs, however, size did not affect following behaviour and larger fish followed as much as smaller fish did. Our results revealed that killifish's following decisions are determined by the size of the individuals, their level of familiarity, and their sex. These characteristics may have a notable impact on the behaviour of fish groups in nature.

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.

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.

Embryonic arsenic exposure reduces intestinal cell proliferation and alters hepatic IGF mRNA expression in killifish (Fundulus heteroclitus)

Arsenic (As) is a toxicant found in food and water throughout the world, and studies suggested that exposure early in life reduces growth. Thus, the goal of this study was to examine mechanisms by which As impacted organismal growth. Killifish (Fundulus heteroclitus) were exposed to 0, 10, 50, or 200 ppb As as embryos and, after hatching, were reared in clean water for up to 40 weeks. Metabolism studies revealed that killifish biotransform As such that monomethylated and dimethylated arsenicals account for 15-17% and 45-61%, respectively, of the total metal. Growth, as measured by condition factor (CF), was significantly and dose-dependently reduced at 8 weeks of age but was similar to controls by 40 weeks. To determine mechanisms underlying the observed initial decrease, intestinal proliferation and morphology were examined. Arsenic-exposed fish exhibited significant 1.3- to 1.5-fold reduction in intestinal villus height and 1.4- to 1.6-fold decrease in proliferating cell nuclear antigen (PCNA+) intestinal cells at all weeks examined. In addition, there were significant correlations between CF, PCNA+ cells, and intestinal villus height. Upon examining whether fish might compensate for the intestinal changes, it was found that hepatic mRNA expression of insulin-like growth factor 1 (IGF-1) and its binding protein (IGFBP-1) were dose-dependently increased. These results indicate that embryonic exposure initially diminished growth, and while intestinal cell proliferation remained reduced, fish appear to compensate by enhancing transcript levels of hepatic IGF-1 and IGFBP-1.

Probabilistic analytical modelling of predator-prey interactions in fishes

Predation is a fundamental interaction between species, yet it is largely unclear what tactics are successful for the survival or capture of prey. One challenge in this area comes with how to test theoretical ideas about strategy with experimental measurements of features such as speed, flush distance and escape angles. Tactics may be articulated with an analytical model that predicts the motion of predator or prey as they interact. However, it may be difficult to recognize how the predictions of such models relate to behavioural measurements that are inherently variable. Here, we present an alternative approach for modelling predator-prey interactions that uses deterministic dynamics, yet incorporates experimental kinematic measurements of natural variation to predict the outcome of biological events. This technique, called probabilistic analytical modelling (PAM), is illustrated by the interactions between predator and prey fish in two case studies that draw on recent experiments. In the first case, we use PAM to model the tactics of predatory bluefish ( Pomatomus saltatrix) as they prey upon smaller fish ( Fundulus heteroclitus). We find that bluefish perform deviated pure pursuit with a variable pursuit angle that is suboptimal for the time to capture. In the second case, we model the escape tactics of zebrafish larvae ( Danio rerio) when approached by adult predators of the same species. Our model successfully predicts the measured patterns of survivorship using measured probability density functions as parameters. As these results demonstrate, PAM is a data-driven modelling approach that can be predictive, offers analytical transparency, and does not require numerical simulations of system dynamics. Though predator-prey interactions demonstrate the use of this technique, PAM is not limited to studying biological systems and has broad utility that may be applied towards understanding a wide variety of natural and engineered dynamical systems where data-driven modelling is beneficial.

Combined effects of flow speed and sub-lethal insecticide exposure on predator-prey interactions between the California killifish and an infaunal polychaete

Hydrodynamics and pollution affect estuarine populations, but their ecological effects have rarely been studied in combination. We conducted two laboratory experiments to quantify whether predator-prey interactions between California killifish, Fundulus parvipinnis, and the polychaete Polydora cornuta vary with flow speed and chlorpyrifos exposure. In one experiment, only F. parvipinnis was exposed to chlorpyrifos; in the other, only P. cornuta was exposed. The flume included a 300-cm² area of sediment with 24 P. cornuta in a central patch (98 cm²). We videotaped groups of three killifish for 50 min at one of four flow speeds (6, 9, 12, or 15 cm/s) and recorded the proportion of bites directed at the prey patch. Unexposed killifish directed 70% of their bites at the prey patch at 6 cm/s, and prey-patch selection decreased as flow increased. Killifish exposed to chlorpyrifos directed 41% of their bites at the prey patch at 6 cm/s with reduced prey-patch selection relative to unexposed fish at 9 and 12 cm/s. At 15 cm/s, both exposed and unexposed fish displayed non-selective biting. Worms were videotaped to quantify their deposit- and suspension-feeding activities. Exposing worms to chlorpyrifos reduced total feeding activity by ~30%. Suspension feeding was more common at faster flow speeds, but the time worms spent suspension feeding relative to deposit feeding was unaffected by chlorpyrifos. No behavioral changes were noted in either species when the other was exposed to chlorpyrifos. This study highlights how hydrodynamic conditions can alter the relative importance of a toxicant's effects on predator-prey interactions.

Profiling microRNA expression in Atlantic killifish (Fundulus heteroclitus) gill and responses to arsenic and hyperosmotic stress

The Atlantic killifish (Fundulus heteroclitus), native to estuarine areas of the Atlantic coast of the United States, has become a valuable ecotoxicological model as a result of its ability to acclimate to rapid environmental changes and adapt to polluted habitats. MicroRNAs (miRNAs) are highly conserved small RNAs that regulate gene expression and play critical roles in stress responses in a variety of organisms. Global miRNA expression in killifish and the potential roles miRNA have in environmental acclimation have yet to be characterized. Accordingly, we profiled miRNA expression in killifish gill for the first time and identified a small group of highly expressed, well-conserved miRNAs as well as 16 novel miRNAs not yet identified in other organisms. Killifish respond to large fluctuations in salinity with rapid changes in gene expression and protein trafficking to maintain osmotic balance, followed by a secondary phase of gene and protein expression changes that enable remodeling of the gills. Arsenic, a major environmental toxicant, was previously shown to inhibit gene expression responses in killifish gill, as well the ability of killifish to acclimate to a rapid increase in salinity. Thus, we examined the individual and combined effects of salinity and arsenic on miRNA expression in killifish gill. Using small RNA sequencing, we identified 270 miRNAs expressed in killifish, and found that miR-135b was differentially expressed in response to arsenic and at 24 h following transfer to salt water. Predicted targets of miR-135b are involved in ion transport, cell motility and migration, GTPase mediated signal transduction and organelle assembly. Consistent with previous studies of these two environmental stressors, we found a significant interaction (i.e., arsenic dependent salinity effect), whereby killifish exposed to arsenic exhibited an opposite response in miR-135b expression at 24 h post hyperosmotic challenge compared to controls. By examining mRNA expression of predicted miRNA targets during salinity acclimation and arsenic exposure, we found that miR-135b targets were significantly more likely to decrease during salinity acclimation than non-targets. Our identification of a significant interaction effect of arsenic and salinity on miR-135b expression supports the hypothesis that arsenic alters upstream regulators of stress response networks, which may adversely affect the killifish response to osmotic stress. The characterization of miRNAs in this ecotoxicological model will be a valuable resource for future studies investigating the role of miRNAs in response to environmental stress.

Improving the reliability and ecological validity of pharmaceutical risk assessment: Turquoise killifish (Nothobranchius furzeri) as a model in behavioral ecotoxicology

Pharmaceuticals are essential for human well-being, but their increasing and continuous use pollutes the environment. Although behavioral ecotoxicology is increasingly advocated to assess the effects of pharmaceutical pollution on wildlife and ecosystems, a consensus on the actual environmental risks is lacking for most compounds. The main limitation is the lack of standardized reproducible tests that are based on sensitive behavioral endpoints and that accommodate a high ecological relevance. In the present study, we assessed the impact of a 3-wk exposure to the antidepressant fluoxetine on multiple behavioral traits in the promising new model organism Nothobranchius furzeri (turquoise killifish). Overall, our study shows that fluoxetine can impact feeding behavior, habitat choice in a novel environment, and antipredator response of N. furzeri individuals; effects on spontaneous activity and exploration tendency were less pronounced. However, effects became only apparent when individuals were exposed to fluoxetine concentrations that were 10 times higher than typical concentrations in natural aquatic environments. Ecotoxicologists are challenged to maximize both the reliability and ecological validity of risk assessments of pollutants. Our study contributes to the development of a time- and cost-efficient, standardized ecotoxicological test based on sensitive, ecologically relevant behavioral endpoints in N. furzeri. Environ Toxicol Chem 2019;38:262-270. © 2018 SETAC.

Toxicity of the biocide polycarbamate, used for aquaculture nets, to some marine fish species

We investigated toxic effects of the antifouling biocide polycarbamate (PC) on marine fish by conducting acute, early-life stage toxicity (ELS), and embryo toxicity tests. Mummichog (Fundulus heteroclitus) 96-h LC₅₀ values for hatched larvae (body weight about 2.0 mg) and juveniles (660 ± 36 mg) were about 12 and 630 μg/L, respectively. The ELS test using mummichog embryos yielded a lowest-observed-effect concentration of 3.9 μg/L and a no-observed-effect concentration of 2.1 μg/L with growth as the most sensitive endpoint. The embryo toxicity test for spotted halibut (Verasper variegatus) revealed a 10-d EC₅₀ of 8.1 μg/L with abnormality as an endpoint. During the ELS and embryo toxicity tests, morphological abnormalities (notochord undulation) were induced in the embryos. Biochemical and gene-expression analysis suggest that PC-induced morphological abnormalities involve disruption of lysyl oxidase-mediated collagen fiber organization, essential for notochord formation, and inhibition of gene expression related to notochord formation_.

Pharmacological evidence that DAPI inhibits NHE2 in Fundulus heteroclitus acclimated to freshwater

Ionoregulation in the euryhaline killifish Fundulus heteroclitus has been intensively studied over the last two decades using a variety of techniques. However, there has been limited use of pharmacological inhibitors to identify proteins involved in ion transport for this species. In this study, we used a range of pharmacological inhibitors (EIPA, DAPI, ethoxzolamide, bumetanide, bafilomycin, phenamil, hydrochlorothiazide) to investigate the proteins involved in Na⁺ transport in freshwater (1 mM Na⁺) acclimated F. heteroclitus. Our results indicate that Na⁺ uptake under these conditions is sensitive to both EIPA (NHE-specific inhibitor) and DAPI (putative ASIC-specific inhibitor), but not to any of the other inhibitors. Results for EIPA are consistent with previous studies indicating F. heteroclitus relies solely on NHE2 for Na⁺ transport across the apical membrane of ionocytes. In contrast, results for DAPI are surprising given previous studies that have indicated the H⁺-ATPase is basolaterally located in F. heteroclitus and so cannot contribute to Na⁺ uptake via ASIC. The lack of bafilomycin sensitivity in the current study is consistent with a basolaterally located H⁺-ATPase. This suggests that DAPI is not an ASIC-specific inhibitor as has been previously assumed, and that it may also inhibit NHE2. Finally, we did not observe Na⁺ uptake to be sensitive to ethoxzolamide, suggesting that carbonic anhydrase may not be involved in generating the H⁺ needed to maintain NHE activity in freshwater as has been previously proposed.

Comparative Efficacy of Pimephales promelas, Fundulus diaphanus, and Gambusia affinis and Influence of Prey Density for Biological Control of Culex pipiens molestus Larvae

Larval survival times and density-dependent feeding behavior were evaluated with the use of 2 species of fish native to the northeastern USA (Pimephales promelas and Fundulus diaphanus), and the potentially invasive Gambusia affinis. Each species was provided 10, 20, 30, 45, and 60 4th-stage larvae of Culex pipiens molestus/fish in the laboratory and digital images were recorded to quantify the number of surviving larvae at various intervals. Daily feeding rates were greatest at the highest larval density. These were 49.69 ± 4.07 larvae for P. promelas, 60 larvae for F. diaphanus, and 36.44 ± 6.6 larvae for G. affinis. Survival analysis was used to compare efficacy of each fish species over time. All fish species consumed larvae at similar rates at lower densities, but significant differences occurred at densities of 30-60 larvae/fish. Survival times of larvae at the highest density were 44 ± 7.9 h for P. promelas, 15 ± 3.4 h for F. diaphanus, and 70.6 ±13 h for G. affinis. In order to evaluate feeding rate as a function of prey density, we compared consumption rates 1.5 h after feeding with the use of a 4-parameter logistic model. Fundulus diaphanus and G. affinis feeding aligned with the 4-parameter model, indicating that initial feeding rates for these species increased with prey density to an upper limit (satiation). Pimephales promelas feeding within 1.5 h did not align with this model, suggesting that early feeding rates for this species are not heavily influenced by prey density.

Embryonic-only arsenic exposure alters skeletal muscle satellite cell function in killifish (Fundulus heteroclitus)

Arsenic is a contaminant found worldwide in drinking water and food. Epidemiological studies have correlated arsenic exposure with reduced weight gain and improper muscular development, while in vitro studies show that arsenic exposure impairs myogenic differentiation. The purpose of this study was to use Fundulus heteroclitus or killifish as a model organism to determine if embryonic-only arsenic exposure permanently reduces the number or function of muscle satellite cells. Killifish embryos were exposed to 0, 50, 200, or 800 ppb arsenite (As^III) until hatching, and then juvenile fish were raised in clean water. At 28, 40, and 52 weeks after hatching, skeletal muscle injuries were induced by injecting cardiotoxin into the trunk of the fish just posterior to the dorsal fin. Muscle sections were collected at 0, 3 and 10 days post-injury. Collagen levels were used to assess muscle tissue damage and recovery, while levels of proliferating cell nuclear antigen (PCNA) and myogenin were quantified to compare proliferating cells and newly formed myoblasts. At 28 weeks of age, baseline collagen levels were 105% and 112% greater in 200 and 800 ppb groups, respectively, and at 52 weeks of age, were 58% higher than controls in the 200 ppb fish. After cardiotoxin injury, collagen levels tend to increase to a greater extent and take longer to resolve in the arsenic exposed fish. The number of baseline PCNA(+) cells were 48-216% greater in 800 ppb exposed fish compared to controls, depending on the week examined. However, following cardiotoxin injury, PCNA is reduced at 28 weeks in 200 and 800 ppb fish at day 3 during the recovery period. By 52 weeks, there are significant reductions in PCNA in all exposure groups at day 3 of the recovery period. Based on these results, embryonic arsenic exposure increases baseline collagen levels and PCNA(+) cells in skeletal muscle. However, when these fish are challenged with a muscle injury, the proliferation and differentiation of satellite cells into myogenic precursors is impaired and instead, the fish appear to be favoring a fibrotic resolution to the injury.

Later life swimming performance and persistent heart damage following subteratogenic PAH mixture exposure in the Atlantic killifish (Fundulus heteroclitus)

High-level, acute exposures to individual polycyclic aromatic hydrocarbons (PAHs) and complex PAH mixtures result in cardiac abnormalities in developing fish embryos. Whereas acute PAH exposures can be developmentally lethal, little is known about the later life consequences of early life, lower level PAH exposures in survivors. A population of PAH-adapted Fundulus heteroclitus from the PAH-contaminated Superfund site, Atlantic Wood Industries, Elizabeth River, Portsmouth, Virginia, United States, is highly resistant to acute PAH cardiac teratogenicity. We sought to determine and characterize long-term swimming performance and cardiac histological alterations of a subteratogenic PAH mixture exposure in both reference killifish and PAH-adapted Atlantic Wood killifish embryos. Killifish from a relatively uncontaminated reference site, King's Creek, Virginia, United States, and Atlantic Wood killifish were treated with dilutions of Elizabeth River sediment extract at 24 h post fertilization (hpf). Two proven subteratogenic dilutions, 0.1 and 1.0% Elizabeth River sediment extract (total PAH 5.04 and 50.4 µg/L, respectively), were used for embryo exposures. Then, at 5-mo post hatching, killifish were subjected to a swim performance test. A separate subset of these individuals was processed for cardiac histological analysis. Unexposed King's Creek killifish significantly outperformed the unexposed Atlantic Wood killifish in swimming performance as measured by Ucrit (i.e., critical swimming speed). However, King's Creek killifish exposed to Elizabeth River sediment extract (both 0.1 and 1.0%) showed significant declines in Ucrit. Histological analysis revealed the presence of blood in the pericardium of King's Creek killifish. Although Atlantic Wood killifish showed baseline performance deficits relative to King's Creek killifish, their pericardial cavities were nearly free of blood and atrial and ventricular alterations. These findings may explain, in part, the diminished swimming performance of King's Creek fish. Environ Toxicol Chem 2017;36:3246-3253. © 2017 SETAC.

A non-destructive BFCOD assay for in vivo measurement of cytochrome P450 3A (CYP3A) enzyme activity in fish embryos and larvae

There is increasing interest in quantifying the exposure and effects of anthropogenic contaminants in fish. Determination of exposures in wild fish is routinely performed, but methods to investigate potential effects are less established. One of the most relevant approaches would be the use of in vivo assays, but existing assays are often limited to in vitro determination of enzyme activity. Many pharmaceuticals and some persistent pollutants activate, and are metabolized by cytochrome P4503A (CYP3A), which make it a relevant and desirable target for biomarker research. We altered the established 7-benzyloxy-4-trifluoromethylcoumarin-O-debenzylation (BFCOD) in vitro protocol for CYP3A activity determination, developing a rapid and inexpensive method to measure in vivo (and in ovo) CYP3A activity in two fish systems: Gulf killifish (Fundulus grandis) and zebrafish (Danio rerio) early life stages. Even with very low concentrations of 7-benzyloxy-4-trifluoromethyl coumarin (BFC, 0.06 µM or 20 µg/L), we were able to detect significant induction in CYP3A activity in embryos of F. grandis, as well as in larvae of D. rerio in response to benzo[a]pyrene (BaP) and fluoranthene (FL) exposures. Because of concerns regarding the possible contribution of CYP1A to BFCOD activity from previous research, we have used a CYP1A post-translational inhibitor (FL) in order to calculate the contribution of CYP1A to the BFCOD assay. We also dosed with benzo[k]fluoranthene (BkF) and showed significant induction of CYP1A activity, with no concurrent increase in CYP3A activity. In this paper, we have taken an established in vitro CYP3A activity assay, and utilized the reaction in a novel way to allow for the non-destructive determination of CYP3A. In summary, we describe a sensitive, cheap, fast and easy modified BFCOD assay for in ovo and in vivo determination of CYP3A activity for use in moderate throughput early-life-stage fish experiments.

Cost of Tolerance: Physiological Consequences of Evolved Resistance to Inhabit a Polluted Environment in Teleost Fish Fundulus heteroclitus

Anthropogenic stressors, including pollutants, are key evolutionary drivers. It is hypothesized that rapid evolution to anthropogenic changes may alter fundamental physiological processes (e.g., energy metabolism), compromising an organism's capacity to respond to additional stressors. The Elizabeth River (ER) Superfund site represents a "natural-experiment" to explore this hypothesis in several subpopulations of Atlantic killifish that have evolved a gradation of resistance to a ubiquitous pollutant-polycyclic aromatic hydrocarbons (PAH). We examined bioenergetic shifts and associated consequences in PAH-resistant killifish by integrating genomic, physiological, and modeling approaches. Population genomics data revealed that genomic regions encoding bioenergetic processes are under selection in PAH-adapted fish from the most contaminated ER site and ex vivo studies confirmed altered mitochondrial function in these fish. Further analyses extending to differentially PAH-resistant subpopulations showed organismal level bioenergetic shifts in ER fish that are associated with increased cost of living, decreased performance, and altered metabolic response to temperature stress-an indication of reduced thermal plasticity. A movement model predicted a higher energetic cost for PAH-resistant subpopulations when seeking an optimum habitat. Collectively, we demonstrate that pollution adaption and inhabiting contaminated environments may result in physiological shifts leading to compromised organismal capacity to respond to additional stressors.

Influence of euthanasia method on blood and gill variables in normoxic and hypoxic Gulf killifish Fundulus grandis

In many experiments, euthanasia, or humane killing, of animals is necessary. Some methods of euthanasia cause death through cessation of respiratory or cardiovascular systems, causing oxygen levels of blood and tissues to drop. For experiments where the goal is to measure the effects of environmental low oxygen (hypoxia), the choice of euthanasia technique, therefore, may confound the results. This study examined the effects of four euthanasia methods commonly used in fish biology (overdose of MS-222, overdose of clove oil, rapid cooling and blunt trauma to the head) on variables known to be altered during hypoxia (haematocrit, plasma cortisol, blood lactate and blood glucose) or reflecting gill damage (trypan blue exclusion) and energetic status (ATP, ADP and ATP:ADP) in Gulf killifish Fundulus grandis after 24 h exposure to well-aerated conditions (normoxia, 7·93 mg O₂  l^-1 , c. 150 mm Hg or c. 20 kPa) or reduced oxygen levels (0·86 mg O₂  l^-1 , c. 17 mm Hg or c. 2·2 kPa). Regardless of oxygen treatment, fish euthanized by an overdose of MS-222 had higher haematocrit and lower gill ATP:ADP than fish euthanized by other methods. The effects of 24 h hypoxic exposure on these and other variables, however, were equivalent among methods of euthanasia (i.e. there were no significant interactions between euthanasia method and oxygen treatment). The choice of an appropriate euthanasia method, therefore, will depend upon the magnitude of the treatment effects (e.g. hypoxia) relative to potential artefacts caused by euthanasia on the variables of interest.

Resistance to polycyclic aromatic hydrocarbon toxicity and associated bioenergetic consequences in a population of Fundulus heteroclitus

Several locations in the Elizabeth River, VA, USA are highly contaminated with polycyclic aromatic hydrocarbons (PAHs) due to the release of creosote mixtures from wood treatment facilities. Interestingly, some populations of Atlantic killifish (Fundulus heteroclitus) inhabiting the Elizabeth River (ER) are resistant to PAH-induced teratogenesis. However, evolutionary resistance to PAHs due to chronic PAH exposure is associated with reduced fitness and increased susceptibility to other environmental stressors in at least one PAH-resistant ER killifish population. More specifically, wild-caught and first generation PAH-resistant juvenile killifish have altered metabolic demands when compared to non-resistant fish. Herein, we investigated this association further by examining a previously under-studied population captured from the creosote-contaminated site Republic Creosoting (Rep). We assessed PAH toxicity and effects on energy metabolism in Rep killifish in comparison with killifish from the reference site Kings Creek (KC). Following exposures to simple and complex PAH mixtures, Rep killifish exhibited several phenotypes associated with PAH resistance including decreased incidences of developmental cardiovascular deformities and recalcitrant cytochrome P450 1A (CYP1A) activity. We evaluated bioenergetics in killifish embryos throughout development and found elevated basal oxygen consumption rates in Rep embryos relative to KC embryos. Furthermore, juvenile F1 Rep fish had significantly lower maximal metabolic rates and aerobic scopes than KC juveniles. These results suggest that populations of killifish that have adapted or evolved to withstand the toxicity associated with PAHs consequently have altered energetic metabolism or demands. Such consequences could result in an enhanced vulnerability to other environmental and anthropogenic stressors in PAH-resistant killifish.

Freshwater influences on embryos, hatching and larval survival of euryhaline Gulf killifish Fundulus grandis and potential constraints on habitat distribution

The influence of fresh water on potential habitat occupancy of early life-history stages of euryhaline Gulf killifish Fundulus grandis was determined by evaluating fertilization of freshwater-spawned eggs and subsequent survival of embryos and larvae in comparison with saline water (salinity 7). Overall per cent fertilization of eggs was low (mean ± s.e. = 20·21 ± 0·03%). Embryo survival was greater in saline water, but hatching rate (mean ± s.e. = 81·6 ± 0·1%) and post-hatch survival of larvae in fresh water (mean ± s.e. = 74·5 ± 0·1%) was relatively high. Therefore, the relative limitation of fresh water on habitat distribution of F. grandis changes with development, stimulating further questions on factors that may constrain habitat distribution of euryhaline fishes.

AHR-related activities in a creosote-adapted population of adult atlantic killifish, Fundulus heteroclitus, two decades post-EPA superfund status at the Atlantic Wood Site, Portsmouth, VA USA

Atlantic killifish, Fundulus heteroclitus, are adapted to creosote-based PAHs at the US EPA Superfund site known as Atlantic Wood (AW) on the southern branch of the Elizabeth River, VA USA. Subsequent to the discovery of the AW population in the early 1990s, these fish were shown to be recalcitrant to CYP1A induction by PAHs under experimental conditions, and even to the time of this study, killifish embryos collected from the AW site are resistant to developmental deformities typically associated with exposure to PAHs in reference fish. Historically, however, 90 +% of the adult killifish at this site have proliferative hepatic lesions including cancer of varying severity. Several PAHs at this site are known to be ligands for the aryl hydrocarbon receptor (AHR). In this study, AHR-related activities in AW fish collected between 2011 and 2013 were re-examined nearly 2 decades after first discovery. This study shows that CYP1A mRNA expression is three-fold higher in intestines of AW killifish compared to a reference population. Using immunohistochemistry, CYP1A staining in intestines was uniformly positive compared to negative staining in reference fish. Livers of AW killifish were examined by IHC to show that CYP1A and AHR2 protein expression reflect lesions-specific patterns, probably representing differences in intrinsic cellular physiology of the spectrum of proliferative lesions comprising the hepatocarcinogenic process. We also found that COX2 mRNA expression levels were higher in AW fish livers compared to those in the reference population, suggesting a state of chronic inflammation. Overall, these findings suggest that adult AW fish are responsive to AHR signaling, and do express CYP1A and AHR2 proteins in intestines at a level above what was observed in the reference population.

Silver toxicity across salinity gradients: the role of dissolved silver chloride species (AgCl x ) in Atlantic killifish (Fundulus heteroclitus) and medaka (Oryzias latipes) early life-stage toxicity

The influence of salinity on Ag toxicity was investigated in Atlantic killifish (Fundulus heteroclitus) early life-stages. Embryo mortality was significantly reduced as salinity increased and Ag(+) was converted to AgCl(solid). However, as salinity continued to rise (>5 ‰), toxicity increased to a level at least as high as observed for Ag(+) in deionized water. Rather than correlating with Ag(+), Fundulus embryo toxicity was better explained (R(2) = 0.96) by total dissolved Ag (Ag(+), AgCl2 (-), AgCl3 (2-), AgCl4 (3-)). Complementary experiments were conducted with medaka (Oryzias latipes) embryos to determine if this pattern was consistent among evolutionarily divergent euryhaline species. Contrary to Fundulus data, medaka toxicity data were best explained by Ag(+) concentrations (R(2) = 0.94), suggesting that differing ionoregulatory physiology may drive observed differences. Fundulus larvae were also tested, and toxicity did increase at higher salinities, but did not track predicted silver speciation. Alternatively, toxicity began to increase only at salinities above the isosmotic point, suggesting that shifts in osmoregulatory strategy at higher salinities might be an important factor. Na(+) dysregulation was confirmed as the mechanism of toxicity in Ag-exposed Fundulus larvae at both low and high salinities. While Ag uptake was highest at low salinities for both Fundulus embryos and larvae, uptake was not predictive of toxicity.

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.

Mummichog (Fundulus heteroclitus) continue to successfully produce eggs after exposure to high levels of 17α-ethinylestradiol

17α-Ethinylestradiol (EE2) is a potent estrogen used in birth-control pills. Previous laboratory and field studies have shown negative impacts in a variety of fish species after exposure to low levels of EE2, most notably a nearly complete shutdown of egg production. The present study demonstrates that mummichog (Fundulus heteroclitus), a small-bodied estuarine species, is able to continue to produce eggs after exposure for 28 d to 100 ng of EE2/L. No effect of EE2 on egg production was observed, whereas a >35-fold increase in vitellogenin (vtg 1) gene expression in males was found. The lack of response in egg production in fish exposed to high levels of EE2 warrants further investigations on species-specific responses to estrogens and endocrine disruptors in general.

Investigating the mechanisms of Ni uptake and sub-lethal toxicity in the Atlantic killifish Fundulus heteroclitus in relation to salinity

The Atlantic killifish (Fundulus heteroclitus) is a resilient estuarine species that may be subjected to anthropogenic contamination of its natural habitat, by toxicants such as nickel (Ni). We investigated Ni accumulation and potential modes of Ni toxicity, in killifish, as a function of environmental salinity. Killifish were acclimated to 4 different salinities [0 freshwater (FW), 10, 30 and 100% seawater (SW)] and exposed to 5 mg/L of Ni for 96 h. Tissue Ni accumulation, whole body ions, critical swim speed and oxidative stress parameters were examined. SW was protective against Ni accumulation in the gills and kidney. Addition of Mg and Ca to FW protected against gill Ni accumulation, suggesting competition with Ni for uptake. Concentration-dependent Ni accumulation in the gill exhibited saturable relationships in both FW- and SW-acclimated fish. However SW fish displayed a lower Bmax (i.e. lower number of Ni binding sites) and a lower Km (i.e. higher affinity for Ni binding). No effect of Ni exposure was observed on critical swim speed (Ucrit) or maximum rate of oxygen consumption (MO2max). Markers of oxidative stress showed either no effect (e.g. protein carbonyl formation), or variable effects that appeared to depend more on salinity than on Ni exposure. These data indicate that the killifish is very tolerant to Ni toxicity, a characteristic that may facilitate the use of this species as a site-specific biomonitor of contaminated estuaries.

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.

The effects of model androgen 5α-dihydrotestosterone on mummichog (Fundulus heteroclitus) reproduction under different salinities

Endocrine disrupting substances (EDSs) have the potential to disturb sensitive hormone pathways, particularly those involved in development and reproduction. Both fresh and estuarine water bodies receive inputs of EDSs from a variety of sources, including sewage effluent, industrial effluent and agricultural runoff. Based on current literature, freshwater species appear to respond to lower levels of EDSs than estuarine or marine species. Therefore, effects elicited by EDSs in freshwater teleosts may not be an accurate representation of how EDSs affect teleosts in estuarine and marine environments. To address this potential difference, a short-term reproductive bioassay was conducted under conditions of low and high salinity using mummichog (Fundulus heteroclitus), a euryhaline species that is native to the east coast of North America. The goals of this study were to determine the response of mummichog when exposed to an androgenic EDS and whether salinity affected the response. A model androgen, 5α-dihydrotestosterone (DHT), was selected for this experiment. Impacts on reproduction were evaluated at multiple biological levels, including physiological (sex steroid levels), organismal (gonad size and gonad morphology), and functional (egg production) endpoints. Under conditions of high salinity, egg production was significantly reduced at all exposure concentrations. Under conditions of low salinity, there were no significant differences based on DHT treatment; however, egg production in all treatment groups including the control were significantly reduced relative to the high salinity control group. Other reproductive endpoints, such as sex steroid production, showed stronger correlation to fecundity in females than males. This study demonstrates that mummichog fecundity is sensitive to androgenic endocrine disruption while also underscoring the importance of how changes in salinity, an environmental variable, can impact reproduction.