Genetic tools for the study of the mangrove killifish, Kryptolebias marmoratus, an emerging vertebrate model for phenotypic plasticity

Kryptolebias marmoratus (Kmar), a teleost fish of the order Cyprinodontiformes, has a suite of unique phenotypes and behaviors not observed in other fishes. Many of these phenotypes are discrete and highly plastic-varying over time within an individual, and in some cases reversible. Kmar and its interfertile sister species, K. hermaphroditus, are the only known self-fertile vertebrates. This unusual sexual mode has the potential to provide unique insights into the regulation of vertebrate sexual development, and also lends itself to genetics. Kmar is easily adapted to the lab and requires little maintenance. However, its internal fertilization and small clutch size limits its experimental use. To support Kmar as a genetic model, we compared alternative husbandry techniques to maximize recovery of early cleavage-stage embryos. We find that frequent egg collection enhances yield, and that protease treatment promotes the greatest hatching success. We completed a forward mutagenesis screen and recovered several mutant lines that serve as important tools for genetics in this model. Several will serve as useful viable recessive markers for marking crosses. Importantly, the mutant kissylips lays embryos at twice the rate of wild-type. Combining frequent egg collection with the kissylips mutant background allows for a substantial enhancement of early embryo yield. These improvements were sufficient to allow experimental analysis of early development and the successful mono- and bi-allelic targeted knockout of an endogenous tyrosinase gene with CRISPR/Cas9 nucleases. Collectively, these tools will facilitate modern developmental genetics in this fascinating fish, leading to future insights into the regulation of plasticity.

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

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

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

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

Turquoise killifish: A natural model of age-dependent brain degeneration

Turquoise killifish (Nothobranchius furzeri) are naturally short-lived vertebrates that display a wide range of spontaneous age-related changes, including onset of cancer, reduced mobility, and cognitive decline. Here, we focus on describing the phenotypic spectrum of the aging killifish brain. As turquoise killifish age, their dopaminergic and noradrenergic neurons undergo a significant decline in number. Furthermore, brain aging in turquoise killifish is associated with several glial-specific changes, such as an increase in the astrocyte-covered surface area and an increase in the numbers of microglial cells, i.e. the brain-specific macrophage population. Killifish brains undergo age-dependent reduced proteasome activity and increased protein aggregation, including the aggregation of the Parkinson's disease marker α-synuclein. Parallel to brain degeneration, turquoise killifish develop spontaneous age-related gut dysbiosis, which has been proposed to affect human neurodegenerative disease. Finally, aged turquoise killifish display declined learning capacity. We argue that, taken together, the molecular, cellular and functional changes that spontaneously take place during aging in killifish brains are consistent with a robust degenerative process that shares remarkable similarities with human neurodegenerative diseases. Hence, we propose that turquoise killifish represent a powerful model of spontaneous brain degeneration which can be effectively used to explore the causal mechanisms underlying neurodegenerative diseases.

The holy grail of longevity research

A new technology to study physiology and cognition elevates African turquoise killifish as a model organism for studies of aging in vertebrates.

The roles of calcium and salinity in protecting against physiological symptoms of waterborne zinc toxicity in the euryhaline killifish (Fundulus heteroclitus)

In fresh water, environmental Ca ameliorates Zn toxicity because Ca²⁺ and Zn²⁺ compete for uptake at the gills. Zn toxicity is also lower in sea water, but it is unclear whether this is due to increased Ca²⁺ concentration, and/or to the other ions present at higher salinity. Using the euryhaline killifish, we evaluated the relative roles of Ca²⁺ (as CaNO₃) versus the other ions contributing to salinity in protecting against physiological symptoms of Zn²⁺ toxicity. Killifish were exposed to a sublethal level of Zn (500 μg/L, as ZnSO₄) for 96 h in either fresh water (0 % salinity) at low (1 mmol/L) and High Ca (10 mmol/L) or 35 ppt sea water (100 % salinity) at low (1 mmol/L) and High Ca (10 mmol/L). At 0 % salinity, High Ca partly or completely protected against the following effects of Zn seen at Low Ca: elevated plasma Zn, hypocalcaemia, inhibited unidirectional Ca²⁺ influx, inhibited branchial Na⁺/K⁺ATPase and Ca²⁺ATPase activities, and oxidative stress in gills, liver, intestine, and muscle. At 100 % salinity, in the presence of 1 mmol/L (Low Ca), Zn caused no disturbances in most of these same parameters, showing that the "non-Ca" components of sea water alone provided complete protection. However, for a few endpoints (inhibited intestinal Ca²⁺ATPase activity, oxidative stress in gill and liver), High Ca (10 mmol/L) was needed to provide full protection against Zn in 100 % salinity. There was no instance where the combination of 100 % salinity and High Ca failed to provide complete protection against Zn-induced disturbances in sea water.

Starter feed for carnivorous species as a practical replacement of bloodworms for a vertebrate model organism in ageing, the turquoise killifish Nothobranchius furzeri

The absence of a controlled diet is unfortunate in a promising model organism for ageing, the turquoise killifish (Nothobranchius furzeri Jubb, 1971). Currently captive N. furzeri are fed bloodworms but it is not known whether this is an optimal diet. Replacing bloodworms with a practical dry feed would reduce diet variability. In the present study, we estimated the nutritional value of the diet ingested by wild fish and determined the fish-body amino acid profile as a proxy for their nutritional requirements. We compared the performance of fish fed four commercial feeds containing 46%-64% protein to that achieved with bloodworms and that of wild fish. Wild fish target a high-protein (60%) diet and this is supported by their superior performance on high-protein diets in captivity. In contrast, feeds for omnivores led to slower growth, lower fecundity and unnatural liver size. In comparison to wild fish, a bloodworm diet led to lower body condition, overfeeding and male liver enlargement. Out of the four dry feeds tested, the fish fed Aller matched wild fish in body condition and liver size, and was comparable to bloodworms in terms of growth and fecundity. A starter feed for carnivorous species appears to be a practical replacement for bloodworms for N. furzeri. The use of dry feeds improved performance in comparison to bloodworms and thus may contribute to reducing response variability and improving research reproducibility in N. furzeri research.

Heart development in two populations of Atlantic killifish (Fundulus heteroclitus) following exposure to a polycyclic aromatic hydrocarbon mixture

Historic industrial pollution of the Elizabeth River, Virginia resulted in polycyclic aromatic hydrocarbon (PAH) contamination in sediments. Atlantic killifish (Fundulus heteroclitus) inhabiting the Atlantic Wood (AW) industrial site adapted to complex PAH mixture at this Superfund site. Their embryos have proved highly resistant to cardiac abnormalities indicative of PAH toxicity. In this study, embryos spawned from adults collected at AW and King's Creek (KC), a reference site, were exposed at 24 h post fertilization (hpf) to Elizabeth River Sediment Extract (ERSE), a complex PAH mixture, in a range of concentrations (0, 5.04, 50.45, 100.90, 151.35, or 252.25 µg/L total PAHs). Embryos were processed for histology at 144 hpf to enable evaluations of hearts at tissue and cellular levels. Morphometry and severity scoring were used to evaluate the extent of alterations. Unexposed embryos were similar in both populations. ERSE exposure resulted in multiple changes to hearts of KC embryos but not AW. Alterations were particularly evident in KC embryos exposed to concentrations above 1% ERSE (50.45 µg/L), which had thinner ventricular walls and larger pericardial edema. Individuals with moderate pericardial edema maintained arrangement and proximity of heart chambers, but changes were seen in ventricular myocytes. Severe pericardial edema was prevalent in exposed KC embryos and typically resulted in tube heart formation. Ventricles of tube hearts had very thin walls composed of small, basophilic cells and lacked trabeculae. Edematous pericardial fluid contained small amounts of proteinaceous material, as did controls, and was free of cells. This fluid was primarily unstained, suggesting water influx due to increased permeability. The use of histological approaches provided more specific detail for tissue and cellular effects in hearts of embryos exposed to PAHs and enabled understanding of potential links to later life effects of early life exposure.

Conspecific density and environmental complexity impact behaviour of turquoise killifish (Nothobranchius furzeri)

Fish models are essential for research in many biological and medical disciplines. With a typical lifespan of only 6 months, the Turquoise killifish (Nothobranchius furzeri) was recently established as a time- and cost-efficient model to facilitate whole-life and multigenerational studies in several research fields, including behavioural ecotoxicology. Essential information on the behavioural norm and on how laboratory conditions affect behaviour, however, is deficient. In the current study, we examined the impact of the social and structural environment on a broad spectrum of behavioural endpoints in N. furzeri. While structural enrichment affected only fish boldness and exploratory behaviour, fish rearing density affected the total body length, locomotor activity, boldness, aggressiveness and feeding behaviour of N. furzeri individuals. Overall, these results contribute to compiling a behavioural baseline for N. furzeri that increases the applicability of this new model species. Furthermore, our findings will fuel the development of improved husbandry protocols to maximize the welfare of N. furzeri in a laboratory setting.

Sex-Specific Differences of Steroid Receptors Following Exposure to Environmentally Relevant Concentrations of Phenothiazine in Fundulus heteroclitus

Phenothiazine (PTZ) is a heterocyclic thiazine compound used for industrial and medical purposes. Through environmental surveillance studies, PTZ was found being discharged into a local river in Connecticut. Phenothiazine has been shown to act similarly to endocrine disrupting chemicals. This study sought to identify sex specific hormone receptor changes in Fundulus heteroclitus in response to PTZ exposure. Fundulus heteroclitus, also known as mummichog, are small fish native to the Atlantic coast of the United States and Canada. They reside in brackish waters and can survive harsh toxic environments. This model organism is native to the polluted waters found in Connecticut. In this study, fish were exposed to PTZ concentrations of 0.5 ppm, 1.0 ppm, and 2.0 ppm for 1 week. Following exposure, brain, liver, and gonad tissues were harvested; cDNA was synthesized; and mRNA expression was assessed for 6 different hormone receptors. Compared with vehicle control (ethanol) differences in mRNA expression, levels of hormone receptors were observed in various tissues from male and female fish. Many of the tissues assessed showed changes in expression level, while only female liver and testis showed no change. These results implicate PTZ as a potential endocrine disrupting compound to mummichog at environmentally relevant concentrations.

Toxicological effects of micronized tire crumb rubber on mummichog (Fundulus heteroclitus) and fathead minnow (Pimephales promelas)

Recent studies on the distribution of microplastics in the Charleston Harbor, SC, USA revealed that a large part of the microplastic particles that are found in the intertidal sediments are tire wear particles. These particles originate from the wear of tire treads on roadways, and wash into the estuary during rain events. The abundance of these particles has raised questions about potential toxicity to aquatic organisms that ingest these particles. The synthetic rubber in car tires consists of a large variety of chemicals, which can vary between manufacturers, but usually contains styrene-butadiene rubber, carbon black and zinc. To investigate the potential toxicity of tire wear particles, both mummichogs (Fundulus heteroclitus) and fathead minnow (Pimephales promelas) were exposed to different concentrations of crumb rubber particles (38-355 µm) in a 7-day static renewal exposure. Dissection of the fish revealed that crumb rubber was ingested and accumulated in the intestinal tract. At the highest concentration tested (6 g/L) partial mortality was observed in the fathead minnow, which came close to the assumed LC₅₀. To investigate if polynuclear aromatic hydrocarbons (PAHs) were leaching from the particles, bile fluorescence was measured, together with potential induction of cytochrome P450-1A through the ethoxyresorufin-O-deethylase (EROD) assay. Elevated levels of 2-, 4-, and 5-, ring structures representative of PAHs were detected in the bile of exposed animals. Bile fluorescence indicated that 4-ring PAH compounds were the most bioavailable from the crumb rubber particles. Induction of EROD activity was observed in exposed animals at environmentally relevant concentrations of the crumb rubber particles (<1-2 g/L), and this elevated EROD activity indicated that PAH compounds from the crumb rubber particles were being metabolized in both mummichogs and fathead minnow.

The effect of population density on growth and survival of a Neotropical annual killifish

We used a field experiment to test the effects of population density on the growth rate and survival of Austrolebias bellottii, a Neotropical annual killifish. Effects differed between the sexes: males at high densities achieved a smaller final size and experienced higher mortality while no such effects were observed in females. This sex-specific effect could be an indirect consequence of mate competition.

Fast-start escape performance across temperature and salinity gradients in mummichog Fundulus heteroclitus

Fast-start predator-escape performance of mummichogs Fundulus heteroclitus was tested across field-informed variation in temperature (24, 30 and 36°C) and salinity (2, 12 and 32 ppt). Performance was similar across temperatures and salinities when fish were allowed to acclimate to these conditions. However, when mummichogs experienced acute temperature changes, performance exhibited thermal dependence in two contrasting ways. Fast-start turning rates and linear speeds varied directly with the temperature at which the manoeuvre was executed, but these aspects of performance varied inversely with acclimation temperature, with cool-acclimated fish exhibiting faster starts across test temperatures. Temperature effects were consistent across salinities. These results suggest that while mummichogs increase performance with acute temperature increases, long-term rises in sea temperature may cause these fish to become more susceptible to predation during abrupt cooling events, such as when storm events flood shallow water estuaries with cool rainwater.

Interspecific variation in hypoxia tolerance and hypoxia acclimation responses in killifish from the family Fundulidae

Hypoxia is a pervasive stressor in aquatic environments, and both phenotypic plasticity and evolutionary adaptation could shape the ability to cope with hypoxia. We investigated evolved variation in hypoxia tolerance and the hypoxia acclimation response across fundulid killifishes that naturally experience different patterns of hypoxia exposure. We compared resting O2 consumption rate (ṀO2 ), and various indices of hypoxia tolerance [critical O2 tension (Pcrit), regulation index (RI), O2 tension (PO2 ) at loss of equilibrium (PLOE) and time to LOE (tLOE) at 0.6 kPa O2] in Fundulus confluentus, Fundulus diaphanus, Fundulus heteroclitus, Fundulus rathbuni, Lucania goodei and Lucania parva We examined the effects of chronic (28 days) exposure to constant hypoxia (2 kPa) or nocturnal intermittent hypoxia (12 h normoxia:12 h hypoxia) in a subset of species. Some species exhibited a two-breakpoint model in ṀO2  caused by early, modest declines in ṀO2  in moderate hypoxia. We found that hypoxia tolerance varied appreciably across species: F. confluentus was the most tolerant (lowest PLOE and Pcrit, longest tLOE), whereas F. rathbuni and F. diaphanus were the least tolerant. However, there was not a consistent pattern of interspecific variation for different indices of hypoxia tolerance, with or without taking phylogenetic relatedness into account, probably because these different indices are underlain by partially distinct mechanisms. Hypoxia acclimation generally improved hypoxia tolerance, but the magnitude of plasticity and responsiveness to different hypoxia patterns varied interspecifically. Our results therefore suggest that hypoxia tolerance is a complex trait that is best appreciated by considering multiple indices of tolerance.

Influence of 96h sub-lethal copper exposure on aerobic scope and recovery from exhaustive exercise in killifish (Fundulus heteroclitus)

Production of industrial effluents have led to increased copper (Cu) pollution of aquatic ecosystems, impacting the physiology of aquatic vertebrates. Past work has shown that Cu exerts its toxicity by disruption ion regulation and/ or increasing oxidative stress. However, it remains unclear how Cu may influence aerobic metabolism and hypoxia tolerance, two possible targets of its toxicity. To address this issue, we exposed freshwater acclimated killifish (F. heteroclitus) to a 96 h Cu exposure at a target concentration of 100 μg L^-1. We determined resting oxygen consumption (ṀO₂), ṀO_2max after exhaustive exercise, and followed ṀO₂ for 3 h in post-exercise recovery in water with either no Cu or 100 μg L^-1 Cu. We assessed hypoxia tolerance by determining the critical oxygen tension (P_crit). It was found that killifish exposed to combined 96 h Cu exposure and Cu present during metabolic measurements, showed a significant decrease in ṀO_2max and in aerobic scope (ṀO_2max - ṀO_2rest), compared to control fish. However, changes in blood and muscle lactate and muscle glycogen were not consistent with an upregulation of anaerobic metabolism as compensation for reduced aerobic performance in Cu exposed fish. Hypoxia tolerance was not influenced by the 96 h Cu exposure or by presence or absence of Cu during the P_crit test. This study suggests that Cu differentially influences responses to changes in oxygen demand and oxygen availability.

Predation drives the evolution of brain cell proliferation and brain allometry in male Trinidadian killifish, Rivulus hartii

The external environment influences brain cell proliferation, and this might contribute to brain plasticity underlying adaptive behavioural changes. Additionally, internal genetic factors influence the brain cell proliferation rate. However, to date, researchers have not examined the importance of environmental versus genetic factors in causing natural variation in brain cell proliferation. Here, we examine brain cell proliferation and brain growth trajectories in free-living populations of Trinidadian killifish, Rivulus hartii, exposed to contrasting predation environments. Compared to populations without predators, populations in high predation (HP) environments exhibited higher rates of brain cell proliferation and a steeper brain growth trajectory (relative to body size). To test whether these differences in the wild persist in a common garden environment, we reared first-generation fish originating from both predation environments in uniform laboratory conditions. Just as in the wild, brain cell proliferation and brain growth in the common garden were greater in HP populations than in no predation populations. The differences in cell proliferation observed across the brain in both the field and common garden studies indicate that the differences are probably genetically based and are mediated by evolutionary shifts in overall brain growth and life-history traits.

The effect of acute warming and thermal acclimation on maximum heart rate of the common killifish Fundulus heteroclitus

Common killifish Fundulus heteroclitus were acclimated to ecologically relevant temperatures (5, 15 and 33°C) and their maximum heart rate (f_Hmax ) was measured at each acclimation temperature during an acute warming protocol. Acclimation to 33°C increased peak f_Hmax by up to 32% and allowed the heart to beat rhythmically at a temperature 10°C higher when compared with acclimation to 5°C. Independent of acclimation temperature, peak f_Hmax occurred about 3°C cooler than the temperature that first produced cardiac arrhythmias. Thus, when compared with previously published values for the critical thermal maximum of F. heteroclitus, the temperature for peak f_Hmax was cooler and the temperature that first produced cardiac arrhythmias was similar to these critical thermal maxima. The considerable thermal plasticity of f_Hmax demonstrated in the present study is entirely consistent with eurythermal ecology of killifish, as shown previously for another eurythermal fish Gillichthys mirabilis.

Physiological effects of 5α-dihydrotestosterone in male mummichog (Fundulus heteroclitus) are dose and time dependent

Numerous anthropogenic sources, such as pulp mill and sewage treatment effluents, contain androgenic endocrine disrupting compounds that alter the reproductive status of aquatic organisms. The current study injected adult male mummichog (Fundulus heteroclitus) with 0 (control), 1 pg/g, 1 ng/g or 1 μg/g body weight of the model androgen 5α-dihydrotestosterone (DHT) with the intent to induce a period of plasma sex hormone depression, a previously-observed effect of DHT in fish. A suite of gonadal steroidogenic genes were assessed during sex hormone depression and recovery. Fish were sampled 6, 12, 16, 18, 24, 30 and 36 h post-injection, and sections of testis tissue were either snap frozen immediately or incubated for 24 h at 18 °C to determine in vitro gonadal hormone production and then frozen. Plasma testosterone (T) and 11-ketotestosterone (11KT) were depressed beginning 24 h post-injection. At 36 h post-injection plasma T remained depressed while plasma 11KT had recovered. In snap frozen tissue there was a correlation between plasma sex hormone depression and downregulation of key steroidogenic genes including steroidogenic acute regulatory protein (star), cytochrome P450 17a1 (cyp17a1), 3β-hydroxysteroid dehydrogenase (3βhsd), 11β-hydroxysteroid dehydrogenase (11βhsd) and 17β-hydroxysteroid dehydrogenase (17βhsd). Similar to previous studies, 3βhsd was the first and most responsive gene during DHT exposure. Gene responses from in vitro tissue were more variable and included the upregulation of 3βhsd, 11βhsd and star during the period of hormone depression. The differential expression of steroidogenic genes from the in vitro testes compared to the snap frozen tissues may be due to the lack of regulators from the hypothalamo-pituitary-gonadal axis present in whole-animal systems. Due to these findings it is recommended to use snap frozen tissue, not post-incubation tissue from in vitro analysis, for gonadal steroidogenic gene expression to more accurately reflect in vivo responses.

Embryonic Fundulus heteroclitus responses to sediment extracts from differentially contaminated sites in the Elizabeth River, VA

Sites along the Elizabeth River are contaminated with polycyclic aromatic hydrocarbons (PAHs) from historical creosote production and other industrial processes. Previous studies have demonstrated that Atlantic killifish collected from sites throughout the Elizabeth River display resistance to the teratogenic effects of PAH-exposure in a manner commensurate with sediment PAH concentrations. The current study characterized various chemical pollutants in sediment and investigated the effects of aqueous sediment extracts from sites along the Elizabeth River to the cardiac development of Atlantic killifish embryos from fish collected from an uncontaminated reference site. Embryonic cardiac deformities were more prevalent after exposure to extracts from sites with high PAH loads. However, activation of cytochrome P4501A, a gene up-regulated by PAH-induction of the aryl hydrocarbon receptor and measured using an in ovo EROD assay, did not consistently increase with PAH concentrations. This work further characterizes sediments in the Elizabeth River, as well as provides insight into the evolutionary pressures at each ER site.

Salinity affects growth but not thermal preference of adult mummichogs Fundulus heteroclitus

The effects of an ecologically relevant range of salinities (2, 12, 22, 32) on thermal preferences and growth of adult mummichogs Fundulus heteroclitus were determined for fish from a southern Chesapeake Bay population. Salinity did not affect the mean temperature selected by F. heteroclitus in a thermal gradient, which was identified as 26.6°C based on observations of 240 individuals. Salinity and temperature had significant and interacting effects on growth rates of F. heteroclitus measured over 12 weeks. Growth rates were highest overall and remained high over a broader range of temperatures at moderate salinities (12 and 22), while high growth rates were shifted toward lower temperatures for fish grown at a salinity of 2 and higher temperatures at a salinity of 32. Significant reductions in growth relative to the optimal conditions (28.6°C, salinity of 22) were observed at the coolest (19.6°C) and warmest (33.6°C) temperature tested at all salinities, as well as temperatures ≥ 26.6°C at a salinity of 2, ≥ 28.6°C at a salinity of 12 and ≤ 26.6°C at a salinity of 32. Growth rates provide a long-term, organismal measure of performance and results of this study indicate that performance may be reduced under conditions that the highly euryhaline F. heteroclitus can otherwise easily tolerate. The combination of reduced salinity and increased temperature that is predicted for temperate estuaries as a result of climate change may have negative effects on growth of this ecologically important species.

Rapid growth and large body size in annual fish populations are compromised by density-dependent regulation

We tested the effect of population density on maximum body size in three sympatric species of annual killifishes Nothobranchius spp. from African ephemeral pools. We found a clear negative effect of population density on body size, limiting their capacity for extremely fast development and rapid growth. This suggests that density-dependent population regulation and the ephemeral character of their habitat impose contrasting selective pressures on the life history of annual killifishes.

Functionalized silver nanoparticles depress aerobic metabolism in the absence of overt toxicity in brackish water killifish, Fundulus heteroclitus

Engineered nanomaterials (ENMs) tend to precipitate in saline waters so the majority of aquatic toxicity studies have focused on freshwaters, where bioavailability is presumed to be higher. Recent studies have illustrated that some ENM formulations are bioavailable and bioactive in salt water and that their effects are more pronounced at the physiological than biochemical level. These findings raise concerns regarding the effects of ENMs on marine organisms. Therefore, our goal was to characterize the effects of polyvinylpyrolidone-functionalized silver ENMs (nAg) on aerobic performance in the killifish (Fundulus heteroclitus), a common euryhaline teleost. Fish were exposed to 80 μg L^-1 of 5 nm nAg for 48 h in brackish water (12 ppt) and routine (ṀO_2min) and maximum (ṀO_2max) rates of oxygen consumption were quantified. Silver dissolution was minimal and nAg remained well dispersed in brackish water, with a hydrodynamic diameter of 21.0 nm, compared to 19.3 in freshwater. Both ṀO_2min and ṀO_2max were significantly lower (by 53 and 30%, respectively) in killifish exposed to nAg and a reduction in ṀO₂ variability suggested spontaneous activity was suppressed. Neither gill Na⁺/K⁺-ATPase activity, nor various other biochemical markers were affected by nAg exposure. The results illustrate that a common ENM formulation is bioactive in salt water and, as in previous studies on functionalized copper ENMs, that effects are more pronounced at the whole animal than the biochemical level.

Interactive effects of 3,4-DCA and temperature on the annual killifish Nothobranchius furzeri

Although aquatic organisms are increasingly exposed to pollutants and abnormally high temperatures as a consequence of climate change, interactive effects between those stressors remain poorly assessed. Especially in ectotherms, such as fish, increases in ambient temperature are expected to affect fitness-related traits and physiology. We used the turquoise killifish Nothobranchius furzeri to study the effects of a range of 3,4-dichloroaniline concentrations (0, 50, 100 μg/L) in combination with two temperature conditions (control and control +4 °C) during four months of exposure. As part of an integrated multi-level approach, we quantified effects on classic life history traits (size, maturation time, body mass, fecundity), critical thermal maximum and physiology (energy reserves and stress-associated enzymatic activity). While no interactive effects of 3,4-DCA exposure and increased temperature emerged, our results do show a negative effect of 3,4-DCA on thermal tolerance. This finding is of particular relevance in light of increasing temperatures under climate change. Due to increases in pest species and faster degradation of 3,4-DCA under higher temperatures, increased use of the pesticide is expected under climate change which, in turn, could result in a decreased tolerance of aquatic organisms to high temperatures.

The influence of salinity on the toxicity of Corexit at multiple life stages of Gulf killifish

Following the Deepwater Horizon oil spill, approximately 7 million liters of the dispersant Corexit 9500A were released to promote oil biodegradation by breaking up surface oil slick formation. This process is accomplished via amphipathic anionic surfactants within dispersants that facilitate the mixing of aqueous and lipid phases. However, the amphipathicity of Corexit may also cause it to interact with biological membranes like the gill, impairing gill function and ultimately disrupting physiological processes mediated by it, such as osmoregulation. The goal of this study was to investigate the osmoregulatory effects and toxicity of Corexit in Gulf killifish. Killifish at the embryonic, larval, juvenile, and adult life stages were exposed to Corexit in water of different salinities to assess the interactive effects of ontogeny and salinity on Corexit toxicity. Corexit was not toxic to embryos except when exposed in hyperosmotic water where it had negligible effects; however, its toxicity to killifish increased dramatically following hatch, showing its greatest deleterious effects in adults. Corexit tended to increase sodium and chloride burdens in killifish when exposed in hyperosmotic waters and reduced whole-body and plasma ion concentrations in fish exposed to hypoosmotic waters. However, Corexit exposure at hyperosmotic salinities resulted in an increased differential accumulation of sodium over chloride as killifish matured. These findings suggest that Corexit may impair gill structure or alter specific components of osmoregulatory function, thus impacting osmoregulation in hypersosmotic and hypoosmotic waters, potentially impairing survival during osmotic challenges. Furthermore, the magnitude of these impacts continues to increase concomitant with gill ontogeny.

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.

A novel variant of aquaporin 3 is expressed in killifish (Fundulus heteroclitus) intestine

Killifish (Fundulus heteroclitus) are euryhaline teleosts that are widely used in environmental and toxicological studies, and they are tolerant to arsenic, in part due to very low assimilation of arsenic from the environment. The mechanism of arsenic uptake by the intestine, a major route of arsenic uptake in humans is unknown. Thus, the goal of this study was to determine if aquaglyceroporins (AQPs), which transport water and other small molecules including arsenite across cell membranes, are expressed in the killifish intestine, and whether AQP expression is affected by osmotic stress. Through RT-PCR and sequence analysis of PCR amplicons, we demonstrated that the intestine expresses kfAQP3a and kfAQP3b, two previously identified variants, and also identified a novel variant of killifish AQP3 (kfAQP3c) in the intestine. The variants likely represent alternate splice forms. A BLAST search of the F. heteroclitus reference genome revealed that the AQP3 gene resides on a single locus, while an alignment of the AQP3 sequence among 384 individuals from eight population ranging from Rhode Island to North Carolina revealed that its coding sequence was remarkably conserved with no fixed polymorphism residing in the region that distinguishes these variants. We further demonstrate that the novel variant transports arsenite into HEK293T cells. Whereas kfAQP3a, which does not transport arsenite, was expressed in both freshwater (FW) and saltwater (SW) acclimated fish, kfAQP3b, an arsenic transporter, was expressed only in FW acclimated fish, and kfAQP3c was expressed only in SW acclimated fish. Thus, we have identified a novel, putative splice variant of kfAQP3, kfAQP3c, which transports arsenic and is expressed only in SW acclimated fish.