Fundulus as the premier teleost model in environmental biology: opportunities for new insights using genomics

Method dependency of maximum oxygen uptake rate and its repeatability in the Gulf killifish, Fundulus grandis

The maximum rate at which fish can take up oxygen from their environment to fuel aerobic metabolism is an important feature of their physiology and ecology. Methods to quantify maximum oxygen uptake rate (ṀO2), therefore, should reliably and reproducibly estimate the highest possible ṀO2 by an individual or species under a given set of conditions (peak ṀO2). This study determined peak ṀO2 and its repeatability in Gulf killifish, Fundulus grandis, subjected to three methods to elevate metabolism: swimming at increasing water speeds, during recovery after an exhaustive chase, and after ingestion of a large meal. Estimates of peak ṀO2 during swimming and after an exhaustive chase were repeatable across two trials, whereas peak ṀO2 after feeding was not. Peak ṀO2 determined by the three methods was significantly different from one another, being highest during swimming, lowest after an exhaustive chase, and intermediate after feeding. In addition, peak ṀO2 during recovery from an exhaustive chase depended on the length of time of recovery: in nearly 60% of the trials, values within the first hour of the chase were lower than those measured later. A novel and important finding was that an individual's peak ṀO2 was not repeatable when compared across methods. Therefore, the peak ṀO2 estimated for a group of fish, as well as the ranking of individual ṀO2 within that group, depends on the method used to elevate aerobic metabolism.

Complete mitochondrial genome of the Starhead Topminnow Fundulus dispar (Cyprinodontiformes: Fundulidae)

Topminnows of the Teleost genus Fundulus serve as model organisms in ecotoxicology because of their broad physiological tolerance and propensity to breed in captivity. This research has been primarily limited to intraspecific comparisons, due to incomplete understanding of the evolutionary history of the genus, which is necessary for use of phylogenetic comparative methods. Interspecific relationships of topminnows remain unresolved, despite recent advances in mitochondrial and nuclear genome sequencing. Specifically, interrelationships of a group containing the starhead topminnows (Fundulus blairae, F. dispar, F. escambiae, F. lineolatus, and F. nottii) typically yield low node support values. Here, we present the first annotated mitochondrial genome of the Starhead Topminnow (F. dispar) and provide a phylogenetic hypothesis for starhead topminnows within the genus Fundulus. DNA was isolated from a specimen of F. dispar collected in Kentucky, USA. The circular genome is 16,564 bp long and contains 13 protein-coding genes, two ribosomal RNAs (rRNA), 22 transfer RNAs (tRNA), and one control region (D-loop). Our phylogenetic analysis supports a sister relationship between F. dispar and a group containing F. notatus and F. olivaceus. This data helps to resolve the phylogenetic placement of starhead topminnows.

Host ecology and biogeography drive parasite community composition in Atlantic killifishes

Understanding the mechanisms of parasite community assembly can be confounded by phylogenetic distance among host species. Addressing this requires focusing on parasite communities within closely related taxa. Thus, we took a macroecological approach to examining parasite community structure within Killifish species in the genus Fundulus to disentangle the effects of host phylogeny and ecological variables. We constructed a database of parasite communities within Fundulus species from 15 published and unpublished surveys covering the Atlantic coast of the US and Canada. The database was expanded by sampling sites in underrepresented provinces and states, totaling 10 Fundulus species from 57 unique geographic sites. Univariate analysis of observed parasite species richness among Fundulus populations in the dataset found that latitude, climate type, and salinity were the dominant factors determining parasite species richness. Multivariate analysis found that host species and landscape type were the most important factors in determining the similarity of parasite assemblages. Unexpectedly, parasite species richness decreased in low latitudes, and host phylogenetic distance was not found to be a significant factor in the similarity of parasite communities. These results indicate that commonly reported large-scale drivers of parasite community structure, such as latitude and phylogeny, could have diminished significance at the host genus level relative to host ecology, biogeography, and local landscape factors.

The effects of winter cold acclimation on acute and chronic cadmium bioaccumulation and toxicity in the banded killifish (Fundulus diaphanus)

Temperate freshwater fishes can experience large seasonal temperature fluctuations that could affect their exposure and sensitivity to trace metals. Yet, temperature effects are overlooked in ecotoxicology studies, especially for cold temperatures typical of the winter. In the present study, the effects of long-term cold acclimation on Cd bioaccumulation and toxicity were investigated in a freshwater fish, the banded killifish (Fundulus diaphanus). Killifish were acclimated to 14 °C or gradually cooled (2 °C/week) to 4 °C and cold acclimated for 6 weeks. Then, both acclimation groups were exposed to environmentally realistic waterborne Cd concentrations (0, 0.5 or 5 µg Cd L-1) for a further 28 d at their respective acclimation temperatures. Tissue metal bioaccumulation, fish survival, condition, and markers of oxidative and ionoregulation stress, were measured after 0, 2, 5 and 28 days of Cd exposure. Cadmium tissue accumulation increased over the exposure duration and was typically lower in cold-acclimated fish. In agreement with this lower bioaccumulation, fewer Cd toxic effects were observed in cold-acclimated fish. There was little evidence of a difference in intrinsic Cd sensitivity between 4 °C- and 14 °C-acclimated fish, as Cd toxicity appeared to closely follow Cd bioaccumulation. Our study suggests that current environmental water quality guidelines would be protective in the winter for the abundant and ecologically-important banded killifish.

A novel gnotobiotic experimental system for Atlantic salmon (Salmo salar L.) reveals a microbial influence on mucosal barrier function and adipose tissue accumulation during the yolk sac stage

Gnotobiotic models have had a crucial role in studying the effect that commensal microbiota has on the health of their animal hosts. Despite their physiological and ecological diversity, teleost fishes are still underrepresented in gnotobiotic research. Moreover, a better understanding of host-microbe interactions in farmed fish has the potential to contribute to sustainable global food supply. We have developed a novel gnotobiotic experimental system that includes the derivation of fertilized eggs of farmed and wild Atlantic salmon, and gnotobiotic husbandry of fry during the yolk sac stage. We used a microscopy-based approach to estimate the barrier function of the skin mucus layer and used this measurement to select the derivation procedure that minimized adverse effects on the skin mucosa. We also used this method to demonstrate that the mucus barrier was reduced in germ-free fry when compared to fry colonized with two different bacterial communities. This alteration in the mucus barrier was preceded by an increase in the number of cells containing neutral mucosubstances in the anterior segment of the body, but without changes in the number of cells containing acidic substances in any of the other segments studied along the body axis. In addition, we showed how the microbial status of the fry temporarily affected body size and the utilization of internal yolk stores during the yolk sac stage. Finally, we showed that the presence of bacterial communities associated with the fry, as well as their composition, affected the size of adipose tissue. Fry colonized with water from a lake had a larger visceral adipose tissue depot than both conventionally raised and germ-free fry. Together, our results show that this novel gnotobiotic experimental system is a useful tool for the study of host-microbe interactions in this species of aquacultural importance.

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

BACKGROUND

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

RESULTS

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

CONCLUSION

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

Endorsing Sustainable Enterprises Among Promising Entrepreneurs: A Comparative Study of Factor-Driven Economy and Efficiency-Driven Economy

Understanding business trails among promising aspirants may contribute to an actual motive for diminishing ecological tracks and escort to developing devotion toward deciding intentions across various entrepreneurial types and tiers solely from the sustainability domains. Therefore, this study endeavors to comprehend and seek to employ the Theory of Planned Behavior (TPB) to inspect the relationship between antecedents on sustainable enterprise intention and sustainable value creation. In this study, we used the convenience sampling method and the quantitative data of 1,070 respondents from Pakistan and China and applied a SmartPLS structural equation model and partial least square path modeling by mediational and multigroup analyses. Findings divulge that R² (79.8%) value in the Pakistan sample of attitudes to sustainability, perceived entrepreneurial desirability (PED), and perceived entrepreneurial feasibility (PEF) was comparatively higher than in China. The R² (75.6%) variance value on sustainable entrepreneurial intentions (SEI) was recorded higher in the Pakistani sample. However, the relationship of environmental values, self-efficacy, and extrinsic and intrinsic rewards show positive and significant mediational effects on both the economies of SEI. The findings disclosed an inconsistent character of extrinsic rewards, general self-efficacy, and job security depict negative significant impacts of aspirations on sustainable entrepreneurship (SE) among promising entrepreneurs on sustainability enterprises in both Pakistan and China. This study extends on existing entrepreneurship literature. Results supported the designed hypothesis and played a significant role in shedding light on an individual trait underpinning a career in a sustainable business start-up. The study looks at the issue from the viewpoint of sustainability domains. It seeks to determine the individual PED, PEF, and attitude toward sustainable entrepreneurship (ATSE) as the mediational variables. The study highlights the importance of work values in pursuing sustainability-oriented entrepreneurship programs for promising aspirants to improve their entrepreneurial skills and knowledge podium, which will encourage them to become sustainable future entrepreneurs. Furthermore, the study provides understandings for ratifying sustainable openings and debates the potential paths for sustainable business growth and opportunities among nascent entrepreneurs in both economies.

Draft genome assemblies using sequencing reads from Oxford Nanopore Technology and Illumina platforms for four species of North American Fundulus killifish

BACKGROUND

Whole-genome sequencing data from wild-caught individuals of closely related North American killifish species (Fundulus xenicus, Fundulus catenatus, Fundulus nottii, and Fundulus olivaceus) were obtained using long-read Oxford Nanopore Technology (ONT) PromethION and short-read Illumina platforms.

FINDINGS

Draft de novo reference genome assemblies were generated using a combination of long and short sequencing reads. For each species, the PromethION platform was used to generate 30-45× sequence coverage, and the Illumina platform was used to generate 50-160× sequence coverage. Illumina-only assemblies were fragmented with high numbers of contigs, while ONT-only assemblies were error prone with low BUSCO scores. The highest N50 values, ranging from 0.4 to 2.7 Mb, were from assemblies generated using a combination of short- and long-read data. BUSCO scores were consistently >90% complete using the Eukaryota database.

CONCLUSIONS

High-quality genomes can be obtained from a combination of using short-read Illumina data to polish assemblies generated with long-read ONT data. Draft assemblies and raw sequencing data are available for public use. We encourage use and reuse of these data for assembly benchmarking and other analyses.

Exposure to wastewater effluent disrupts hypoxia responses in killifish (Fundulus heteroclitus)

Hypoxia (low oxygen) often occurs in aquatic ecosystems that receive effluent from municipal wastewater treatment plants (WWTP). The combination of hypoxia and WWTP effluent could impair fish health, because WWTP effluent contains multiple contaminants that could disrupt the physiological pathways fish use to cope with hypoxia, but the interactive effects of these stressors on fish physiology are poorly understood. We have examined this issue by exposing mummichog killifish (Fundulus heteroclitus) to hypoxia (5 and 2 kPa O₂) and/or 100% WWTP effluent for 21 days in a full factorial design. We then measured hypoxia tolerance, whole-animal metabolism, gill morphology, haematology, and tissue metabolites. In clean water, killifish responded to chronic hypoxia with improvements in hypoxia tolerance, as reflected by increases in time to loss of equilibrium at 0.5 kPa (t_LOE). These improvements occurred in association with increases in the exposed surface of gill lamellae that resulted from a regression of interlamellar cell mass (ILCM). Concurrent exposure to wastewater attenuated the increases in t_LOE and gill remodeling in chronic hypoxia, and nearly depleted brain glycogen stores. Therefore, exposure to WWTP effluent can disrupt the physiological mechanisms fish use to cope with chronic hypoxia and impair hypoxia tolerance. Our research suggests that the combination of stressors near WWTPs can have interactive effects on the physiology and health of fish.

Mummichog gill and operculum exhibit functionally consistent claudin-10 paralog profiles and Claudin-10c hypersaline response

Claudin (Cldn)-10 tight junction (TJ) proteins are hypothesized to form the paracellular Na⁺ secretion pathway of hyposmoregulating mummichog (Fundulus heteroclitus) branchial epithelia. Organ-specific expression profiles showed that only branchial organs [the gill and opercular epithelium (OE)] exhibited abundant cldn-10 paralog transcripts, which typically increased following seawater (SW) to hypersaline (2SW) challenge. Post-translational properties, protein abundance, and ionocyte localization of Cldn-10c, were then examined in gill and OE. Western blot analysis revealed two Cldn-10c immunoreactive bands in the mummichog gill and OE at ∼29 kDa and ∼40 kDa. The heavier protein could be eliminated by glycosidase treatment, demonstrating the novel presence of a glycosylated Cldn-10c. Protein abundance of Cldn-10c increased in gill and OE of 2SW-exposed fish. Cldn-10c localized to the sides of gill and OE ionocyte apical crypts and partially colocalized with cystic fibrosis transmembrane conductance regulator and F-actin, consistent with TJ complex localization. Cldn-10c immunofluorescent intensity increased but localization was unaltered by 2SW conditions. In support of our hypothesis, cldn-10/Cldn-10 TJ protein dynamics in gill and OE of mummichogs and TJ localization are functionally consistent with the creation and maintenance of salinity-responsive, cation-selective pores that facilitate Na⁺ secretion in hyperosmotic environments.

Summary: The role of claudin-10 tight junction proteins in paracellular salt secretion across fish branchial epithelia is indicated by organ-specific responses to hyperosmotic conditions and their association with salt secreting transcellular proteins

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

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

The chemical defensome of five model teleost fish

How an organism copes with chemicals is largely determined by the genes and proteins that collectively function to defend against, detoxify and eliminate chemical stressors. This integrative network includes receptors and transcription factors, biotransformation enzymes, transporters, antioxidants, and metal- and heat-responsive genes, and is collectively known as the chemical defensome. Teleost fish is the largest group of vertebrate species and can provide valuable insights into the evolution and functional diversity of defensome genes. We have previously shown that the xenosensing pregnane x receptor (pxr, nr1i2) is lost in many teleost species, including Atlantic cod (Gadus morhua) and three-spined stickleback (Gasterosteus aculeatus), but it is not known if compensatory mechanisms or signaling pathways have evolved in its absence. In this study, we compared the genes comprising the chemical defensome of five fish species that span the teleosteii evolutionary branch often used as model species in toxicological studies and environmental monitoring programs: zebrafish (Danio rerio), medaka (Oryzias latipes), Atlantic killifish (Fundulus heteroclitus), Atlantic cod, and three-spined stickleback. Genome mining revealed evolved differences in the number and composition of defensome genes that can have implication for how these species sense and respond to environmental pollutants, but we did not observe any candidates of compensatory mechanisms or pathways in cod and stickleback in the absence of pxr. The results indicate that knowledge regarding the diversity and function of the defensome will be important for toxicological testing and risk assessment studies.

Genomic Signature of Shifts in Selection and Alkaline Adaptation in Highland Fish

Understanding how organisms adapt to aquatic life at high altitude is fundamental in evolutionary biology. This objective has been addressed primarily related to hypoxia adaptation by recent comparative studies, whereas highland fish has also long suffered extreme alkaline environment, insight into the genomic basis of alkaline adaptation has rarely been provided. Here, we compared the genomes or transcriptomes of 15 fish species, including two alkaline tolerant highland fish species and their six alkaline intolerant relatives, three alkaline tolerant lowland fish species, and four alkaline intolerant species. We found putatively consistent patterns of molecular evolution in alkaline tolerant species in a large number of shared orthologs within highland and lowland fish taxa. Remarkably, we identified consistent signatures of accelerated evolution and positive selection in a set of shared genes associated with ion transport, apoptosis, immune response, and energy metabolisms in alkaline tolerant species within both highland and lowland fish taxa. This is one of the first comparative studies that began to elucidate the consistent genomic signature of alkaline adaptation shared by highland and lowland fish. This finding also highlights the adaptive molecular evolution changes that support fish adapting to extreme environments at high altitude.

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.

20th Pollutant Responses in Marine Organisms (PRIMO 20): Global issues and fundamental mechanisms caused by pollutant stress in marine and freshwater organisms

The 20th Pollutant Responses in Marine Organisms (PRIMO 20) conference provided a forum for scientists from around the world to communicate novel toxicological research findings specifically focused on aquatic organisms, by combining applied and basic research at the intersection of environmental and mechanistic toxicology. The work highlighted in this special issue of Aquatic Toxicology, a special issue of Marine Environmental Research, and presented through posters and presentations, encompass important and emerging topics in freshwater and marine toxicology. This includes multiple types of emerging contaminants including microplastics and UV filtering chemicals. Other studies aimed to further our understanding of the effects of endocrine disrupting chemicals, pharmaceuticals, and personal care products. Further research presented in this virtual issue examined the interactive effects of chemicals and pathogens, while the final set of manuscripts demonstrates continuing efforts to combine traditional biomonitoring, data from -omic technologies, and modeling for use in risk assessment and management. An additional goal of PRIMO meetings is to address the link between environmental and human health. Several articles in this issue of Aquatic Toxicology describe the appropriateness of using aquatic organisms as models for human health, while the keynote speakers, as described in the editorial below, presented research that highlighted bioaccumulation of contaminants such as PFOS and mercury from fish to marine mammals and coastal human populations such as the Gullah/GeeChee near Charleston, South Carolina, USA.

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.

Plasticity, repeatability and phenotypic correlations of aerobic metabolic traits in a small estuarine fish

Standard metabolic rate (SMR), maximum metabolic rate (MMR), absolute aerobic scope (AAS) and critical oxygen tension (P_crit) were determined for the Gulf killifish, Fundulus grandis, an ecologically dominant estuarine fish, acclimated to lowered salinity, elevated temperature and lowered oxygen concentration. Acclimation to low salinity resulted in a small, but significant, elevation of P_crit (suggesting lower tolerance of hypoxia); acclimation to elevated temperature increased SMR, MMR, AAS and P_crit; acclimation to low oxygen led to a small increase in SMR, but substantial decreases in MMR, AAS and P_crit Variation in these metabolic traits among individuals was consistent and repeatable when measured during multiple control exposures over 7 months. Trait repeatability was unaffected by acclimation condition, suggesting that repeatability of these traits is not context dependent. There were significant phenotypic correlations between specific metabolic traits: SMR was positively correlated with MMR and P_crit; MMR was positively correlated with AAS; and AAS was negatively correlated with P_crit In general, within-individual variation contributed more than among-individual variation to these phenotypic correlations. The effects of acclimation on these traits demonstrate that aerobic metabolism is plastic and influenced by the conditions experienced by these fish in the dynamic habitats in which they occur; however, the repeatability of these traits and the correlations among them suggest that these traits change in ways that maintain the rank order of performance among individuals across a range of environmental variation.

Plasticity, repeatability and phenotypic correlations of aerobic metabolic traits in a small estuarine fish

Standard metabolic rate (SMR), maximum metabolic rate (MMR), absolute aerobic scope (AAS) and critical oxygen tension (Pcrit) were determined for the Gulf killifish, Fundulus grandis, an ecologically dominant estuarine fish, acclimated to lowered salinity, elevated temperature and lowered oxygen concentration. Acclimation to low salinity resulted in a small, but significant, elevation of Pcrit (suggesting lower tolerance of hypoxia); acclimation to elevated temperature increased SMR, MMR, AAS and Pcrit; acclimation to low oxygen led to a small increase in SMR, but substantial decreases in MMR, AAS and Pcrit Variation in these metabolic traits among individuals was consistent and repeatable when measured during multiple control exposures over 7 months. Trait repeatability was unaffected by acclimation condition, suggesting that repeatability of these traits is not context dependent. There were significant phenotypic correlations between specific metabolic traits: SMR was positively correlated with MMR and Pcrit; MMR was positively correlated with AAS; and AAS was negatively correlated with Pcrit In general, within-individual variation contributed more than among-individual variation to these phenotypic correlations. The effects of acclimation on these traits demonstrate that aerobic metabolism is plastic and influenced by the conditions experienced by these fish in the dynamic habitats in which they occur; however, the repeatability of these traits and the correlations among them suggest that these traits change in ways that maintain the rank order of performance among individuals across a range of environmental variation.

Hypoxia acclimation alters reactive oxygen species homeostasis and oxidative status in estuarine killifish (Fundulus heteroclitus)

Hypoxia is common in aquatic environments, and exposure to hypoxia followed by re-oxygenation is often believed to induce oxidative stress. However, there have been relatively few studies of reactive oxygen species (ROS) homeostasis and oxidative status in fish that experience natural hypoxia-re-oxygenation cycles. We examined how exposure to acute hypoxia (2 kPa O2) and subsequent re-oxygenation (to 20 kPa O2) affects redox status, oxidative damage and anti-oxidant defenses in estuarine killifish (Fundulus heteroclitus), and whether these effects were ameliorated or potentiated by prolonged (28 days) acclimation to either constant hypoxia or intermittent cycles of nocturnal hypoxia (12 h:12 h normoxia:hypoxia). Acute hypoxia and re-oxygenation led to some modest and transient changes in redox status, increases in oxidized glutathione, depletion of scavenging capacity and oxidative damage to lipids in skeletal muscle. The liver had greater scavenging capacity, total glutathione concentrations and activities of anti-oxidant enzymes (catalase, glutathione peroxidase) than muscle, and generally experienced less variation in glutathiones and lipid peroxidation. Unexpectedly, acclimation to constant hypoxia or intermittent hypoxia led to a more oxidizing redox status (muscle and liver) and it increased oxidized glutathione (muscle). However, hypoxia-acclimated fish exhibited little to no oxidative damage (as reflected by lipid peroxidation and aconitase activity), in association with improvements in scavenging capacity and catalase activity in muscle. We conclude that hypoxia acclimation leads to adjustments in ROS homeostasis and oxidative status that do not reflect oxidative stress, but may instead be part of the suite of responses that killifish use to cope with chronic hypoxia.

Using asexual vertebrates to study genome evolution and animal physiology: Banded (Fundulus diaphanus) x Common Killifish (F. heteroclitus) hybrid lineages as a model system

Wild, asexual, vertebrate hybrids have many characteristics that make them good model systems for studying how genomes evolve and epigenetic modifications influence animal physiology. In particular, the formation of asexual hybrid lineages is a form of reproductive incompatibility, but we know little about the genetic and genomic mechanisms by which this mode of reproductive isolation proceeds in animals. Asexual lineages also provide researchers with the ability to produce genetically identical individuals, enabling the study of autonomous epigenetic modifications without the confounds of genetic variation. Here, we briefly review the cellular and molecular mechanisms leading to asexual reproduction in vertebrates and the known genetic and epigenetic consequences of the loss of sex. We then specifically discuss what is known about asexual lineages of Fundulus diaphanus x F. heteroclitus to highlight gaps in our knowledge of the biology of these clones. Our preliminary studies of F. diaphanus and F. heteroclitus karyotypes from Porter's Lake (Nova Scotia, Canada) agree with data from other populations, suggesting a conserved interspecific chromosomal arrangement. In addition, genetic analyses suggest that: (a) the same major clonal lineage (Clone A) of F. diaphanus x F. heteroclitus has remained dominant over the past decade, (b) some minor clones have also persisted, (c) new clones may have recently formed, and iv) wild clones still mainly descend from F. diaphanus ♀ x F. heteroclitus ♂ crosses (96% in 2017-2018). These data suggest that clone formation may be a relatively rare, but continuous process, and there are persistent environmental or genetic factors causing a bias in cross direction. We end by describing our current research on the genomic causes and consequences of a transition to asexuality and the potential physiological consequences of epigenetic variation.

Potential osmoprotective roles of branchial heat shock proteins towards Na+, K+-ATPase in milkfish (Chanos chanos) exposed to hypotonic stress

In euryhaline teleosts, osmoregulatory mechanisms vary with osmotic stresses, and heat shock proteins (HSPs) play a central role in maintaining cellular homeostasis. The present study aimed to investigate the expression and potential roles of HSP70 and HSP90 in the gills of seawater (SW)- and freshwater (FW)-acclimated milkfish (Chanos chanos). Four HSP genes, including Cchsc70 (heat shock cognate 70), Cchsp70, Cchsp90α, and Cchsp90β, were analyzed in milkfish gills. Among these genes, only the mRNA abundance of branchial Cchsp90α was significantly lower in the FW-acclimated than in the SW-acclimated milkfish. Immunoblotting showed no significant difference in the relative protein abundance of branchial HSP70 and HSP90 between the two groups. The time-course experiments (from SW to FW) showed that the protein abundance of HSP70 and HSP90 at the 3 h and 6 h post-transfer and then declined gradually. To further illustrate the potential osmoregulatory roles of HSP70 and HSP90, their interaction with Na⁺, K⁺-ATPase (NKA, the primary driving force for osmoregulation) was analyzed using co-immunoprecipitation. The results showed the interaction between HSP70, HSP90 and NKA after acclimation to SW or FW increased within 3 h; and then returned to normal levels within 7 days. To our knowledge, the present study was the first to demonstrate that the interaction between HSP70, HSP90 and NKA changes with hypotonic stress in euryhaline teleosts. Before the transfer, no interaction was detected. When transferred to FW from SW, the interaction of HSP70 and HSP90 with NKA were detected. The results suggested that HSP70 and HSP90 participated in the acute responses of osmoregulatory mechanisms to protect branchial NKA from hypotonic stress in milkfish.

Long-Term Ecological Impacts from Oil Spills: Comparison of Exxon Valdez, Hebei Spirit, and Deepwater Horizon

The long-term ecological impacts of the Exxon Valdez oil spill (EVOS) are compared to two extensively studied and more recent large spills: Deepwater Horizon (DWH) and the Hebei Spirit oil spill (HSOS). Each of the three spills differed in magnitude and duration of oil released, environmental conditions, ecological communities, response and clean up measures, and ecological recovery. The EVOS began on March 24, 1989, and released 40.8 million liters of Alaska North Slope crude oil into the cold, nearly pristine environment of Prince William Sound, Alaska. EVOS oiled wildlife and rocky intertidal shorelines and exposed early life stages of fish to embryotoxic levels of polycyclic aromatic hydrocarbons (PAH). Long-term impacts following EVOS were observed on seabirds, sea otters, killer whales, and subtidal communities. The DWH spill began on April 20, 2010, and released 507 million liters of light Louisiana crude oil from 1600 m on the ocean floor into the Gulf of Mexico over an 87-day period. The DWH spill exposed a diversity of complex aquatic communities in the deep ocean, offshore pelagic areas, and coastal environments to petroleum hydrocarbons. Large-scale persistent ecological effects included impacts to deep ocean corals, failed recruitment of oysters over multiple years, damage to coastal wetlands, and reduced dolphin, sea turtle, and seabird populations. The HSOS began on December 7, 2007, and released approximately 13 million liters of Middle East crude oils into ecologically sensitive areas of the Taean area of western Korea. Environmental conditions and the extensive initial cleanup of HSOS oil stranded on shorelines limited the long-term impacts to changes in composition and abundance of intertidal benthic communities. Comparisons of EVOS, DWH, and HSOS show the importance and complexity of the interactions among the environment, oil spill dynamics, affected ecological systems, and response actions.

Mummichog (Fundulus heteroclitus) are less sensitive to 17α-ethinylestradiol (EE2) than other common model teleosts: A comparative review of reproductive effects

The environmental estrogen 17α-ethinylestradiol (EE₂) will depress or completely inhibit egg production in many common model teleosts at low concentrations (≤0.5 ng/L; Runnalls et al., 2015). This inhibition is not seen in the estuarine killifish, or mummichog (Fundulus heteroclitus), even when exposed to 100 ng/L EE₂. This relative insensitivity to EE₂ exposure indicates species-specific mechanisms for compensating for exogenous estrogenic exposure. This review compares various reproductive responses elicited by EE₂ in mummichog to other common model teleosts, such as zebrafish (Danio rerio) and fathead minnow (Pimephales promelas), identifying key endpoints where mummichog differ from other studied fish. For example, EE₂ accumulates primarily in the liver/gall bladder of mummichog, which is different than zebrafish and fathead minnow in which accumulation is predominantly in the carcass. Despite causing species-specific differences in fecundity, EE₂ has been shown to consistently induce hepatic vitellogenin in males and cause feminization/sex reversal during gonadal differentiation in larval mummichog, similar to other species. In addition, while gonadal steroidogenesis and plasma steroid levels respond to exogenous EE_2, it is generally at higher concentrations than observed in other species. In mummichog, production of 17β-estradiol (E₂) by full grown ovarian follicles remains high; unlike other teleost models where E₂ synthesis decreases as 17α,20β-dihydroxy-4-prenen-3-on levels increase to induce oocyte maturation. New evidence in mummichog indicates some dissimilarity in gonadal steroidogenic gene expression responses compared to gene expression responses in zebrafish and fathead minnow exposed to EE₂. The role of ovarian physiology continues to warrant investigation regarding the tolerance of mummichog to exogenous EE₂ exposure. Here we present a comprehensive review, highlighting key biological differences in response to EE₂ exposure between mummichog and other commonly used model teleosts.

Establishment of a brain cell line (FuB-1) from mummichog (Fundulus heteroclitus) and its application to fish virology, immunity and nanoplastics toxicology

The marine fish mummichog (Fundulus heteroclitus), extensively used as research model, including in ecotoxicology, for over a century has been surpassed by other fish species. This fact may be associated with the lack of cell lines from this species, excellent models for the comprehension of fish physiology, immunology, toxicology and virology, that contribute to the reduction in the number of animals used in research. We have generated, for the first time, a brain-derived cell line from mummichog, FuB-1, and evaluated its application to the fields of fish virology, immunity and toxicology. First, FuB-1 cells show epithelial morphology and neural stem/astroglial origin. Secondly, FuB-1 cells effectively supports the replication of both spring viremia carp (SVCV) and infectious pancreatic necrosis (IPNV) viruses, but not nodavirus (NNV), indicating its potential use for fish virology. Related to this, FuB-1 cells infected with NNV up-regulate the transcription of genes related to the antiviral immune response, leading to cell resistance; while they are unaltered when infected with IPNV and SVCV, facilitating viral replication. Finally, FuB-1 cells were used for toxicological purposes and we demonstrated that exposure to either polystyrene nanoplastics (PS-100) or several human-usage pharmaceuticals are cytotoxic. Additionally, PS-100 particles increase the antioxidant catalase and glutathione S-transferase activities and decrease the total non-protein thiols in FuB-1 cells. However, PS-100 particles are able to reduce the cytotoxic effects induced by the pharmaceuticals. In conclusion, we have generated a cell line from mummichog, which might represent a valuable model for fish studies in the fields of virology, immunology and toxicology.

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.

Salinity-dependent expression of ncc2 in opercular epithelium and gill of mummichog (Fundulus heteroclitus)

Mummichogs (Fundulus heteroclitus) can tolerate abrupt changes in environmental salinity because of their ability to rapidly adjust the activities of ionocytes in branchial and opercular epithelia. In turn, the concerted expression of sub-cellular effectors of ion transport underlies adaptive responses to fluctuating salinities. Exposure to seawater (SW) stimulates the expression of Na⁺/K⁺/2Cl^- cotransporter 1 (nkcc1) and cystic fibrosis transmembrane regulator (cftr) mRNAs in support of ion extrusion by SW-type ionocytes. Given the incomplete understanding of how freshwater (FW)-type ionocytes actually operate in mummichogs, the transcriptional responses essential for ion absorption in FW environments remain unresolved. In a subset of species, a 'fish-specific' Na⁺/Cl^- cotransporter denoted Ncc2 (Slc12a10) is responsible for the uptake of Na⁺ and Cl^- across the apical surface of FW-type ionocytes. In the current study, we identified an ncc2 transcript that is highly expressed in gill filaments and opercular epithelium of FW-acclimated mummichogs. Within 1 day of transfer from SW to FW, ncc2 levels in both tissues increased in parallel with reductions in nkcc1 and cftr. Conversely, mummichogs transferred from FW to SW exhibited marked reductions in ncc2 concurrent with increases in nkcc1 and cftr. Immunohistochemical analyses employing a homologous antibody revealed apical Ncc2-immunoreactivity in Na⁺/K⁺-ATPase-immunoreactive ionocytes of FW-acclimated animals. Our combined observations suggest that Ncc2/ncc2-expressing ionocytes support the capacity of mummichogs to inhabit FW environments.

Focal adhesion kinase and osmotic responses in ionocytes of Fundulus heteroclitus, a euryhaline teleost fish

Cystic Fibrosis Transmembrane conductance Regulator (CFTR) anion channels are the regulated exit pathway in Cl^- secretion by teleost salt secreting ionocytes of the gill and opercular epithelia of euryhaline teleosts. By confocal light immunocytochemistry using regular and phospho-antibodies directed against conserved sites, we found that killifish CFTR (kfCFTR) and the tyrosine kinase Focal Adhesion Kinase (FAK) phosphorylated at Y407 (FAKpY407) and FAKpY397 are colocalized at the apical membrane and in subjacent membrane vesicles of ionocytes. Hypotonic shock and the α-2 adrenergic agonist clonidine rapidly and reversibly inhibit Cl^- secretion by isolated opercular epithelia, simultaneous with dephosphorylation of FAKpY407 and increased FAKpY397, located in the apical membrane of ionocytes in the opercular epithelium. FAKpY407 is re-phosphorylated at the apical membrane of ionocytes and Cl^- secretion rapidly restored by hypertonic shock, detectable at 2 min., maximum at 5 min and still elevated at 30 min. In isolated opercular epithelia, the FAK phosphorylation inhibitor Y15 and p38MAP kinase inhibitor SB203580 significantly blunted the recovery of short-circuit current (I_sc, equal to Cl^- secretion rate) after hypertonic shock. The cSRC inhibitor saracatinib dephosphorylated FAKpY861 seen near tight junctions of pavement cells, and reduced the increase in epithelial resistance normally seen with clonidine inhibition of ion transport, while FAKpY397 was unaffected. The results show rapid osmosensitive responses in teleost fish ionocytes involve phosphorylation of CFTR by FAKpY407, an opposing role for FAKpY397 and a possible role for FAKpY861 in tight junction dynamics.

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.

The osmorespiratory compromise: physiological responses and tolerance to hypoxia are affected by salinity acclimation in the euryhaline Atlantic killifish (Fundulus heteroclitus)

The characteristics of the fish gill that maximize gas exchange are the same that promote diffusion of ions and water to and from the environment; therefore, physiological trade-offs are likely to occur. Here, we investigated how salinity acclimation affects whole-animal respiratory gas exchange during hypoxia using Fundulus heteroclitus, a fish that inhabits salt marshes where salinity and oxygen levels vary greatly. Salinity had marked effects on hypoxia tolerance, with fish acclimated to 11 and 35 ppt showing much longer time to loss of equilibrium (LOE) in hypoxia than 0 ppt-acclimated fish. Fish acclimated to 11 ppt (isosmotic salinity) exhibited the greatest capacity to regulate oxygen consumption rate (Ṁ _O2 ) under hypoxia, as measured through the regulation index (RI) and P _crit At 35 ppt, fish had a higher routine metabolic rate (RMR) but a lower RI than fish at 11 ppt, but there were no differences in gill morphology, ventilation or blood O₂ transport properties between these groups. In contrast, 0 ppt-acclimated fish had the highest ventilation and lowest O₂ extraction efficiency in normoxia and hypoxia, indicating a higher ventilatory workload in order to maintain similar levels of Ṁ _O2 These differences were related to alterations in gill morphology, where 0 ppt-acclimated fish had the smallest lamellar surface area with the greatest epithelial cell coverage (i.e. thicker lamellae, longer diffusion distance) and a larger interlamellar cell mass, contrasting with 11 ppt-acclimated fish, which had overall the highest respiratory surface area. The alteration of an array of physiological parameters provides evidence for a compromise between salinity and hypoxia tolerance in killifish acclimated to freshwater.

Mitochondrial Ecophysiology: Assessing the Evolutionary Forces That Shape Mitochondrial Variation

The mitonuclear species concept hypothesizes that incompatibilities between interacting gene products of the nuclear and mitochondrial genomes are a major factor establishing and maintaining species boundaries. However, most of the data available to test this concept come from studies of genetic variation in mitochondrial DNA, and clines in the mitochondrial genome across contact zones can be produced by a variety of forces. Here, we show that using a combination of population genomic analyses of the nuclear and mitochondrial genomes and studies of mitochondrial function can provide insight into the relative roles of neutral processes, adaptive evolution, and mitonuclear incompatibility in establishing and maintaining mitochondrial clines, using Atlantic killifish (Fundulus heteroclitus) as a case study. There is strong evidence for a role of secondary contact following the last glaciation in shaping a steep mitochondrial cline across a contact zone between northern and southern subspecies of killifish, but there is also evidence for a role of adaptive evolution in driving differentiation between the subspecies in a variety of traits from the level of the whole organism to the level of mitochondrial function. In addition, studies are beginning to address the potential for mitonuclear incompatibilities in admixed populations. However, population genomic studies have failed to detect evidence for a strong and pervasive influence of mitonuclear incompatibilities, and we suggest that polygenic selection may be responsible for the complex patterns observed. This case study demonstrates that multiple forces can act together in shaping mitochondrial clines, and illustrates the challenge of disentangling their relative roles.

The osmorespiratory compromise in the euryhaline killifish: water regulation during hypoxia

Freshwater- and seawater-acclimated Fundulus heteroclitus were exposed to acute hypoxia (10% air saturation, 3 h), followed by normoxic recovery (3 h). In both salinities, ventilation increased and heart rate fell in the classic manner, while Ṁ _O2  initially declined by ∼50%, with partial restoration by 3 h of hypoxia, and no O₂ debt repayment during recovery. Gill paracellular permeability (measured with [¹⁴C] PEG-4000) was 1.4-fold higher in seawater, and declined by 50% during hypoxia with post-exposure overshoot to 188%. A similar pattern with smaller changes occurred in freshwater. Drinking rate (also measured with [¹⁴C] PEG-4000) was 8-fold higher in seawater fish, but declined by ∼90% during hypoxia in both groups, with post-exposure overshoots to ∼270%. Gill diffusive water flux (measured with ³H₂O) was 1.9-fold higher in freshwater fish, and exhibited a ∼35% decrease during hypoxia, which persisted throughout recovery, but was unchanged during hypoxia in seawater fish. Nevertheless, freshwater killifish gained mass while seawater fish lost mass during hypoxia, and these changes were not corrected during normoxic recovery. We conclude that this hypoxia-tolerant teleost beneficially reduces gill water permeability in a salinity-dependent fashion during hypoxia, despite attempting to simultaneously improve Ṁ _O2 , but nevertheless incurs a net water balance penalty in both freshwater and seawater.

Developing a High-Quality Linkage Map for the Atlantic Killifish Fundulus heteroclitus

Killifish (Fundulus heteroclitus) are widely distributed among different aquatic environments where they demonstrate an impressive range of highly-plastic and locally adaptive phenotypes. High-throughput sequencing has begun to unravel the mechanisms and evolutionary history of these interesting features by establishing relationships in the genotype-phenotype map. However, some genotype-phenotype analyses require a higher order of contiguity than what initial scaffolded (fragmented genome assembly where contigs have been assemble into scaffolds) genome assemblies can provide. Here, we used 5,685 high-quality RAD-Seq markers from a single mapping family to order 84% of the scaffolded genome assembly to 24 chromosomes. This serves to: 1) expand the killifish genomic toolkit, 2) estimate genome-wide recombination rates, and 3) compare genome synteny to humans and other fishes. After initially building our map, we found that the selection of thresholds for sequence data filtration highly impacted scaffold placement in the map. We outline each step of the approach that dramatically improved our map to help guide others toward more effective linkage mapping for genome assembly. Our final map supports strong conservation of genomic synteny among closely related fish species and reveals previously described chromosomal rearrangements between more distantly related clades. However, we also commonly found minor scaffold misorientations in F. heteroclitus and in other assemblies, suggesting that further mapping (such as optical mapping) is necessary for finer scale resolution of genome structure. Lastly, we discuss the problems that would be expected from misoriented/unplaced scaffolds and stress the importance of a quality mapped genome as a key feature for further investigating population and comparative genomic questions with F. heteroclitus and other taxa.

Exploring the advantages and limitations of sampling methods commonly used in research facilities for zebrafish health inspections

Examining zebrafish populations for the presence of disease is an integral component of managing fish health in research facilities. Currently, many different strategies are used for zebrafish fish health inspections, which is a scenario that may result in subjective and biased diagnostic evaluations. The goal of this study was to compare the success of pathogen detection between a sample size of randomly selected fish (n = 60) that provides 95% confidence in pathogen detection based on a presumed pathogen prevalence level ≥5%, and other subpopulations and sample numbers commonly submitted for diagnostic testing within a 1000 tank, 30,000 fish, recirculating research system. This included fish collected from a sump tank (n = 53), sentinel fish (n = 11), and fish that were found moribund or freshly dead (n = 18). Additionally, five fish from each subpopulation were collected for histopathologic examination. A second study used retrospective data to examine pathogen distribution between systems (n = 2-5) in multi-system facilities (n = 5) using a sample size of 60 fish per system. For the pathogens detected, results supported the use of representative sample numbers rather than smaller numbers of populations considered more at risk. The exception to this is for the moribund/mortality group, which may be a resource for targeted surveillance of select pathogens. Each system within multi-system facilities should be considered separate units in terms of fish health inspections and biosecurity. Development of these evidence-based standards for fish health inspections in zebrafish systems enhances fish welfare, provides identification of potentially zoonotic pathogens, and ensures scientific integrity and reproducibility of research results.

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

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

Exposure to arsenic during embryogenesis impairs olfactory sensory neuron differentiation and function into adulthood

Arsenic is a contaminant of food and drinking water. Epidemiological studies have reported correlations between arsenic exposure and neurodevelopmental abnormalities, such as reduced sensory functioning, while in vitro studies have shown that arsenic reduces neurogenesis and alters stem cell differentiation. The goal of this study was assess whether arsenic exposure during embryogenesis reduced olfactory stem cell function and/or numbers, and if so, whether those changes persist into adulthood. Killifish (Fundulus heteroclitus) embryos were exposed to 0, 10, 50 or 200 ppb arsenite (As^III) until hatching, and juvenile fish were raised in clean water. At 0, 2, 4, 8, 16, 28 and 40 weeks of age, odorant response tests were performed to assess specific olfactory sensory neuron (OSN) function. Olfactory epithelia were then collected for immunohistochemical analysis of stem cell (Sox2) and proliferating cell numbers (PCNA), as well as the number and expression of ciliated (calretinin) and microvillus OSNs (G_αi3) at 0, 4, 16 and 28 weeks. Odorant tests indicated that arsenic exposure during embryogenesis increased the start time of killifish responding to pheromones, and this altered start time persisted to 40 weeks post-exposure. Response to the odorant taurocholic acid (TCA) was also reduced through week 28, while responses to amino acids were not consistently altered. Immunohistochemistry was used to determine whether changes in odorant responses were correlated to altered cell numbers in the olfactory epithelium, using markers of proliferating cells, progenitor cells, and specific OSNs. Comparisons between response to pheromones and PCNA + cells indicated that, at week 0, both parameters in exposed fish were significantly reduced from the control group. At week 28, all exposure are still significantly different than control fish, but now with higher PCNA expression coupled with reduced pheromone responses. A similar trend was seen in the comparisons between Sox2-expressing progenitor cells and response to pheromones, although Sox2 expression in the 28 week-old fish only recovers back to the level of control fish rather than being significantly higher. Comparisons between calretinin expression (ciliated OSNs) and response to TCA demonstrated that both parameters were reduced in the 200 ppb arsenic-exposed fish in at weeks 4, 16, and 28. Correlations between TCA response and the number of PCNA + cells revealed that, at 28 weeks of age, all arsenic exposure groups had reductions in response to TCA, but higher PCNA expression, similar to that seen with the pheromones. Few changes in G_αi3 (microvillus OSNs) were seen. Thus, it appears that embryonic-only exposure to arsenic has long-term reductions in proliferation and differentiation of olfactory sensory neurons, leading to persistent effects in their function.

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.

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

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

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.

Mechanisms of Ca2+ uptake in freshwater and seawater-acclimated killifish, Fundulus heteroclitus, and their response to acute salinity transfer

Killifish (Fundulus heteroclitus) has been extensively used as a model for ion regulation by euryhaline fishes. Na⁺ and Cl^- dynamics have been well studied in killifish, but few studies have addressed that of Ca²⁺. Therefore, this study aimed to characterize Ca²⁺ fluxes in freshwater (FW) and seawater (SW)-acclimated killifish, their response to salinity transfer, and to elucidate the mechanisms of Ca²⁺ influx in FW and SW. SW killifish displayed a significantly higher Ca²⁺ influx rate than that of FW fish, while Ca²⁺ efflux rates were comparable in both salinities. Ca²⁺ influx was saturable in FW (K_m = 78 ± 19 µmol/L; J_max = 53 ± 3 nmol/g/h) and influx by SW killifish was linear up to 7 mmol/L Ca²⁺. In SW-acclimated fish, 36% of Ca²⁺ influx was attributed to "intestinal Ca²⁺ intake", likely caused by drinking, whereas intestinal Ca²⁺ intake in FW contributed to < 2% of total. Throughout the study, results suggested that "cation competition" in SW modulates Ca²⁺ influx. Therefore, we hypothesized that SW-acclimated fish actually have a higher affinity Ca²⁺ influx system than FW-acclimated fish but that it is competitively inhibited by competing SW cations. In agreement with this cation competition hypothesis, we demonstrated for the first time that "extra-intestinal" Ca²⁺ influx was inhibited by Mg²⁺ in both FW and SW-acclimated killifish. Following acute salinity transfer, extra-intestinal Ca²⁺ influx was rapidly regulated within 12-24 h, similar to Na⁺ and Cl^-. Ca²⁺ influx in FW was inhibited by La³⁺, an epithelial Ca²⁺ channel blocker, whereas La³⁺ had no significant effect in SW.

Distinct metabolic adjustments arise from acclimation to constant hypoxia and intermittent hypoxia in estuarine killifish (Fundulus heteroclitus)

Many fish experience daily cycles of hypoxia in the wild, but the physiological strategies for coping with intermittent hypoxia are poorly understood. We examined how killifish adjust O2 supply and demand during acute hypoxia, and how these responses are altered after prolonged acclimation to constant or intermittent patterns of hypoxia exposure. We acclimated killifish to normoxia (∼20 kPa O2), constant hypoxia (2 kPa) or intermittent cycles of nocturnal hypoxia (12 h:12 h normoxia:hypoxia) for 28 days, and then compared whole-animal O2 consumption rates (ṀO2) and tissue metabolites during exposure to 12 h of hypoxia followed by reoxygenation in normoxia. Normoxia-acclimated fish experienced a pronounced 27% drop in ṀO2 during acute hypoxia, and modestly increased ṀO2 upon reoxygenation. They strongly recruited anaerobic metabolism during acute hypoxia, indicated by lactate accumulation in plasma, muscle, liver, brain, heart and digestive tract, as well as a transient drop in intracellular pH, and they increased hypoxia inducible factor (HIF)-1α protein abundance in muscle. Glycogen, glucose and glucose-6-phosphate levels suggested that glycogen supported brain metabolism in hypoxia, while the muscle used circulating glucose. Acclimation to constant hypoxia caused a stable ∼50% decrease in ṀO2 that persisted after reoxygenation, with minimal recruitment of anaerobic metabolism, suggestive of metabolic depression. By contrast, fish acclimated to intermittent hypoxia maintained sufficient O2 transport to support normoxic ṀO2, modestly recruited lactate metabolism and increased ṀO2 dramatically upon reoxygenation. Both groups of hypoxia-acclimated fish had similar glycogen, ATP, intracellular pH and HIF-1α levels as normoxic controls. We conclude that different patterns of hypoxia exposure favour distinct strategies for matching O2 supply and O2 demand.

The fallacy of the P crit - are there more useful alternatives?

P _crit - generally defined as the P _O2  below which the animal can no longer maintain a stable rate of O₂ consumption (Ṁ _O2 ), such that Ṁ _O2  becomes dependent upon P _O2  - provides a single number into which a vast amount of experimental effort has been invested. Here, with specific reference to water-breathers, I argue that this focus on the P _crit is not useful for six reasons: (1) calculation of P _crit usually involves selective data editing; (2) the value of P _crit depends greatly on the way it is determined; (3) there is no good theoretical justification for the concept; (4) P _crit is not the transition point from aerobic to anaerobic metabolism, and it disguises what is really going on; (5) P _crit is not a reliable index of hypoxia tolerance; and (6) P _crit carries minimal information content. Preferable alternatives are loss of equilibrium (LOE) tests for hypoxia tolerance, and experimental description of full Ṁ _O2  versus P _O2  profiles accompanied by measurements of ventilation, lactate appearance and metabolic rate by calorimetry. If the goal is to assess the ability of the animal to regulate Ṁ _O2  from this profile in a mathematical fashion, promising, more informative alternatives to P _crit are the regulation index and Michaelis-Menten or sigmoidal allosteric analyses.

Can variation among hypoxic environments explain why different fish species use different hypoxic survival strategies?

In aquatic environments, hypoxia is a multi-dimensional stressor that can vary in O₂ level (partial pressure of O₂ in water, Pw_O2 ), rate of induction and duration. Natural hypoxic environments can therefore be very different from one another. For the many fish species that have evolved to cope with these different hypoxic environments, survival requires adjusting energy supply and demand pathways to maintain energy balance. The literature describes innumerable ways that fishes combine aerobic metabolism, anaerobic metabolism and metabolic rate depression (MRD) to accomplish this, but it is unknown whether the evolutionary paths leading to these different strategies are determined primarily by species' phylogenetic histories, genetic constraint or their native hypoxic environments. We explored this idea by devising a four-quadrant matrix that bins different aquatic hypoxic environments according to their duration and Pw_O2  characteristics. We then systematically mined the literature for well-studied species native to environments within each quadrant, and, for each of 10 case studies, described the species' total hypoxic response (THR), defined as its hypoxia-induced combination of sustained aerobic metabolism, enhanced anaerobic metabolism and MRD, encompassing also the mechanisms underlying these metabolic modes. Our analysis revealed that fishes use a wide range of THRs, but that distantly related species from environments within the same matrix quadrant have converged on similar THRs. For example, environments of moderately hypoxic Pw_O2  favoured predominantly aerobic THRs, whereas environments of severely hypoxic Pw_O2  favoured MRD. Capacity for aerial emergence as well as predation pressure (aquatic and aerial) also contributed to these responses, in addition to other biotic and abiotic factors. Generally, it appears that the particular type of hypoxia experienced by a fish plays a major role in shaping its particular THR.

Transcriptomic responses to extreme low salinity among locally adapted populations of Olympia oyster (Ostrea lurida)

The Olympia oyster (Ostrea lurida) is a foundation species inhabiting estuaries along the North American west coast. In California estuaries, O. lurida is adapted to local salinity regimes and populations differ in low salinity tolerance. In this study, oysters from three California populations were reared for two generations in a laboratory common garden and subsequently exposed to low salinity seawater. Comparative transcriptomics was then used to understand species-level responses to hyposmotic stress and population-level mechanisms underlying divergent salinity tolerances. Gene expression patterns indicate Olympia oysters are sensitive to hyposmotic stress: All populations respond to low salinity by up-regulating transcripts indicative of protein unfolding, DNA damage and cell cycle arrest after sub-lethal exposure. Among O. lurida populations, transcriptomic profiles differed constitutively and in response to low salinity. Despite two generations in common-garden conditions, transcripts encoding apoptosis modulators were constitutively expressed at significantly different levels in the most tolerant population. Expression of cell death regulators may facilitate cell fate decisions when salinity declines. Following low salinity exposure, oysters from the more tolerant population expressed a small number of mRNAs at significantly higher levels than less tolerant populations. Proteins encoded by these transcripts regulate ciliary activity within the mantle cavity and may function to prolong valve closure and reduce mortality in low salinity seawater. Collectively, gene expression patterns suggest sub-lethal impacts of hyposmotic stress in Olympia oysters are considerable and that even oysters with greater low salinity tolerance may be vulnerable to future freshwater flooding events.

Bioaccumulation and effects of dietary exposure to the alternative flame retardant, bis(2-ethylhexyl) tetrabromophthalate (TBPH), in the Atlantic killifish, Fundulus heteroclitus

Bis(2-ethylhexyl) tetrabromophthalate (TBPH), a high production volume flame retardant chemical used as a replacement for banned flame retardants, has been detected in media and human and wildlife tissues globally. We describe bioaccumulation and biological effects from dietary exposure of TBPH to an estuarine fish, Atlantic killifish, Fundulus heteroclitus. Briefly, adult fish were fed carrier control or chemically amended diets for 28 d, followed by 14 d of control diet feeding. Diets were amended with TBPH (TBPH_LO diet, 139 μg/g dry wt, or TBPH_HI diet, 4360 μg/g dry wt) or a polychlorinated biphenyl congener (PCB153 diet, 13 μg/g dry wt), which was included as a positive control for bioaccumulation. Although bioaccumulation of either chemical correlated with fish size, only a small proportion of the TBPH offered (<0.5% total TBPH) had bioaccumulated into TBPH-treated fish by 28 d. In contrast, 24.5% of the PCB153 offered was accounted for in 28-d PCB-treated fish. Although 28-d bioaccumulated concentrations of TBPH differed by sex and treatment, sexes did not differ in their rates of TBPH bioaccumulation, and the time to achieve 50% of 28 d concentration (T1/2 ) was estimated to be 13 d. Depuration rates of TBPH did not differ by sex or treatment, and the time after exposure to achieve T1/2 was estimated to be 22 d. Independent of treatment, male fish grew faster than female fish, but for both sexes reproductive condition (gonadal somatic index) declined unexpectedly over the experimental period. Across treatments, only the TBPH_LO treatment affected growth, reducing male but increasing female growth rates by small amounts relative to respective controls. In summary, our study used very high concentrations of dietary TBPH to contaminate fish tissues above the highest levels reported to date in wild biota, yet we observed few adverse biological effects. Environ Toxicol Chem 2018;37:2350-2360. © 2018 SETAC. This article is a US government work and, as such, is in the public domain in the United States of America.

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

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

Early ecotoxic effects of ZnO nanoparticle chronic exposure in Mytilus galloprovincialis revealed by transcription of apoptosis and antioxidant-related genes

Recently, China became one of the largest nanomaterial markets in the world. The wide use of ZnO nanoparticles in a number of products implies an increasing release in marine environment and consequently the evaluation of the potential effects upon marine organisms largely cultured in China for commercial purposes, such as invertebrate bivalves is a current need. To this aim, survival, bioaccumulation, and transcription pattern of key genes, p53, PDRP, SOD, CAT, and GST, involved in DNA damage/repair and antioxidation, in Mytilus galloprovincialis digestive gland, exposed to ZnO NPs (<100 nm) and ZnO bulk (150-200 nm) for 4 weeks, were evaluated. ZnSO₄ was also assessed to appraise the role of zinc ions. Starting from 72 h, increasing mortality values along the exposure time were observed for all ZnO compounds. The highest difference was evident after 28 d when NPs resulted three times more toxic than bulk, (LC₅₀) = 0.78 mg Zn/L (confidence limits: 0.64, 1.00) and 2.62 mg Zn/L (confidence limits: 1.00, 4.00), respectively. For ZnSO₄ the (LC₅₀) was always the lowest reaching the minimum value at 28 d 0.25 mg Zn/L (confidence limits: 0.10-0.40). Digestive gland showed higher uptake rate of ionic Zn respect to ZnO NPs and bulk during the first three days of exposure. In particular at the end of the exposure time (28 d) at 1 mg Zn/L the rank of Zn uptake rate was Zinc ion > ZnO NPs > ZnO bulk. The relative expression of investigated genes evidenced that distinct actions of apoptosis and antioxidation occurred in M. galloprovincialis exposed to ZnO NPs with a peculiar pattern dependent on exposure time and concentration. Application of the qRT-PCR technique revealed evidence of sensitivity to the nanomaterial since the first time of exposure.

Vertebrate Genome Evolution in the Light of Fish Cytogenomics and rDNAomics

To understand the cytogenomic evolution of vertebrates, we must first unravel the complex genomes of fishes, which were the first vertebrates to evolve and were ancestors to all other vertebrates. We must not forget the immense time span during which the fish genomes had to evolve. Fish cytogenomics is endowed with unique features which offer irreplaceable insights into the evolution of the vertebrate genome. Due to the general DNA base compositional homogeneity of fish genomes, fish cytogenomics is largely based on mapping DNA repeats that still represent serious obstacles in genome sequencing and assembling, even in model species. Localization of repeats on chromosomes of hundreds of fish species and populations originating from diversified environments have revealed the biological importance of this genomic fraction. Ribosomal genes (rDNA) belong to the most informative repeats and in fish, they are subject to a more relaxed regulation than in higher vertebrates. This can result in formation of a literal 'rDNAome' consisting of more than 20,000 copies with their high proportion employed in extra-coding functions. Because rDNA has high rates of transcription and recombination, it contributes to genome diversification and can form reproductive barrier. Our overall knowledge of fish cytogenomics grows rapidly by a continuously increasing number of fish genomes sequenced and by use of novel sequencing methods improving genome assembly. The recently revealed exceptional compositional heterogeneity in an ancient fish lineage (gars) sheds new light on the compositional genome evolution in vertebrates generally. We highlight the power of synergy of cytogenetics and genomics in fish cytogenomics, its potential to understand the complexity of genome evolution in vertebrates, which is also linked to clinical applications and the chromosomal backgrounds of speciation. We also summarize the current knowledge on fish cytogenomics and outline its main future avenues.

claudin-10 isoform expression and cation selectivity change with salinity in salt-secreting epithelia of Fundulusheteroclitus

To provide insight into claudin (Cldn) tight junction (TJ) protein contributions to branchial salt secretion in marine teleost fishes, this study examined cldn-10 TJ protein isoforms of a euryhaline teleost (mummichog; Fundulus heteroclitus) in association with salinity change and measurements of transepithelial cation selectivity. Mummichogs were transferred from freshwater (FW) to seawater (SW, 35‰) and from SW to hypersaline SW (2SW, 60‰) in a time course with transfer control groups (FW to FW, and SW to SW). FW to SW transfer increased mRNA abundance of cldn-10d and cldn-10e twofold, whilst cldn-10c and cldn-10f transcripts were unchanged. Transfer from SW to 2SW did not alter cldn-10d, and transiently altered cldn-10e abundance, but increased cldn-10c and cldn-10f fourfold. This was coincident with an increased number of single-stranded junctions (observed by transmission electron microscopy). For both salinity transfers, (1) cldn-10e mRNA was acutely responsive (i.e. after 24 h), (2) other responsive cldn-10 isoforms increased later (3-7 days), and (3) cystic fibrosis transmembrane conductance regulator (cftr) mRNA was elevated in accordance with established changes in transcellular Cl^- movement. Changes in mRNA encoding cldn-10c and -10f appeared linked, consistent with the tandem repeat locus in the Fundulus genome, whereas mRNA for tandem cldn-10d and cldn-10e seemed independent of each other. Cation selectivity sequence measured by voltage and conductance responses to artificial SW revealed Eisenman sequence VII: Na⁺>K⁺>Rb⁺∼Cs⁺>Li⁺ Collectively, these data support the idea that Cldn-10 TJ proteins create and maintain cation-selective pore junctions in salt-secreting tissues of teleost fishes.