A common flame-retardant and plasticizer, triphenyl phosphate (TPhP) is an aryl phosphate ester found in many aquatic environments at nM concentrations. Yet, most studies interrogating its toxicity have used µM concentrations. In this study, we used the model organism zebrafish (Danio rerio) to uncover the developmental impact of nM exposures to TPhP at the phenotypic and molecular levels. At concentrations of 1.5–15 nM (0.5 µg/L–5 µg/L), chronically dosed 5dpf larvae were shorter in length and had pericardial edema phenotypes that had been previously reported for exposures in the µM range. Cardiotoxicity was observed but did not present as cardiac looping defects as previously reported for µM concentrations. The RXR pathway does not seem to be involved at nM concentrations, but the tbx5a transcription factor cascade including natriuretic peptides (nppa and nppb) and bone morphogenetic protein 4 (bmp4) were dysregulated and could be contributing to the cardiac phenotypes. We also demonstrate that TPhP is a weak pro-oxidant, as it increases the oxidative stress response within hours of exposure. Overall, our data indicate that TPhP can affect animal development at environmentally relevant concentrations and its mode of action involves multiple pathways.

Freshwater animals are exposed to anthropogenic contaminants and are biomonitors of water quality and models of the deleterious impacts of exposure. Sponges, such as Ephydatia muelleri, constantly pump water and are effective indicators of water-soluble contaminants. Zebrafish (Danio rerio), native to Southeast Asia, live in the water column and feed at the water-sediment interface and are exposed to both water-soluble and insoluble contaminants. While sponges and zebrafish diverged ∼700 million years ago, they share common genetic elements, and their response to contaminants can be predictive to a wide-range of animals. An emerging contaminant, bisphenol S, was tested to evaluate its toxicity during development. The toxicity and mechanism(s) of action of BPS is not well known. Water-borne exposures to BPS caused differing hatching rates, morphological changes, and shared gene expression changes of toxicologically-relevant genes. This study shows that BPS causes similarly adverse developmental impacts pointing to some overlapping mechanisms of action.

Mono-2-ethylhexyl phthalate (MEHP) is the primary metabolite of the ubiquitous plasticizer and toxicant, di-2-ethylhexyl phthalate. MEHP exposure has been linked to abnormal development, increased oxidative stress, and metabolic syndrome in vertebrates. Nuclear factor, Erythroid 2 Like 2 (Nrf2), is a transcription factor that regulates gene expression in response to oxidative stress. We investigated the role of Nrf2a in larval steatosis following embryonic exposure to MEHP. Wild-type and nrf2a mutant (m) zebrafish embryos were exposed to 0 or 200 μg/l MEHP from 6 to either 96 (histology) or 120 hours post fertilization (hpf). At 120 hpf, exposures were ceased and fish were maintained in clean conditions until 15 days post fertilization (dpf). At 15 dpf, fish lengths and lipid content were examined, and the expression of genes involved in the antioxidant response and lipid processing was quantified. At 96 hpf, a subset of animals treated with MEHP had vacuolization in the liver. At 15 dpf, deficient Nrf2a signaling attenuated fish length by 7.7%. MEHP exposure increased hepatic steatosis and increased expression of peroxisome proliferator-activated receptor alpha target fabp1a1 . Cumulatively, these data indicate that developmental exposure alone to MEHP may increase risk for hepatic steatosis and that Nrf2a does not play a major role in this phenotype.

Interactions between regulatory pathways allow organisms to adapt to their environment and respond to stress. One interaction that has been recently identified occurs between the aryl hydrocarbon receptor (AHR) and the nuclear factor erythroid-2 related factor (NRF) family. Each transcription factor regulates numerous downstream genes involved in the cellular response to toxicants and oxidative stress; they are also implicated in normal developmental pathways. The zebrafish model was used to explore the role of AHR regulation of nrf genes during development and in response to toxicant exposure. To determine if AHR1b is responsible for transcriptional regulation of 6 nrf genes during development, a loss-of-function experiment using morpholino-modified oligonucleotides was conducted followed by a chromatin immunoprecipitation study at the beginning of the pharyngula period (24 h postfertilization). The expression of nrf1a was AHR1b dependent and its expression was directly regulated through specific XREs in its cis-promoter. However, nrf1a expression was not altered by exposure to 2, 3, 7, 8-tetrachlorodibenzo-p-dioxin (TCDD), a toxicant and prototypic AHR agonist. The expression of nrf1b, nrf2a, and nfe2 was induced by TCDD, and AHR1b directly regulated their expression by binding to cis-XRE promoter elements. Last, nrf2b and nrf3 were neither induced by TCDD nor regulated by AHR1b. These results show that AHR1b transcriptionally regulates nrf genes under toxicant modulation via binding to specific XREs. These data provide a better understanding of how combinatorial molecular signaling potentially protects embryos from embryotoxic events following toxicant exposure.

Mono(2-ethylhexyl) phthalate (MEHP) is the bioactive metabolite of di(2-ethylhexyl) phthalate, a plasticizing agent and persistent environmental contaminant associated with obesity, developmental abnormalities, and oxidative stress. Nrf2 (Nfe2l2) is a transcription factor that regulates cytoprotective genes as part of the adaptive antioxidant response. We previously identified the pancreas as a sensitive target of oxidative stress during embryonic development. The goals of this study were to 1) characterize the effects of MEHP exposure on pancreatic development, and 2) determine whether oxidative stress contributes to MEHP embryotoxicity. Zebrafish (Danio rerio) embryos from AB wildtype and Tg(ins:GFP;nrf2afh318/fh318) were exposed to 0 or 200 μg/L MEHP at 3 h post fertilization (hpf) through 168 hpf to assess pancreatic organogenesis. MEHP exposure significantly decreased β-cell area at all timepoints (48, 72, 96, 168 hpf), but Nrf2a did not significantly protect against islet hypomorphism. Tg(gcga:GFP) embryos exposed to MEHP showed a decrease in α-cell area in the islet across the same timepoints. Tg(ptf1a:GFP) embryos were assessed at 80 and 168 hpf for exocrine pancreas length. MEHP exposure decreased growth of the exocrine pancreas. Expression of pancreas genes insa, sst2 and ptf1a was significantly reduced by MEHP exposure compared to controls. Glutathione (GSH) concentrations and redox potentials were quantified at 72 hpf by HPLC, but no significant changes were observed. However, expression of the GSH-related genes gstp1 and gsr were significantly altered by MEHP exposure. These data indicate that the developing pancreas is a sensitive target tissue of embryonic exposure to MEHP.

Redox signaling is important for embryogenesis, guiding pathways that govern processes crucial for embryo patterning, including cell polarization, proliferation, and apoptosis. Exposure to pro-oxidants during this period can be deleterious, resulting in altered physiology, teratogenesis, later-life diseases, or lethality. We previously reported that the glutathione antioxidant defense system becomes increasingly robust, including a doubling of total glutathione and dynamic shifts in the glutathione redox potential at specific stages during embryonic development in the zebrafish, Danio rerio. However, the mechanisms underlying these changes are unclear, as is the effectiveness of the glutathione system in ameliorating oxidative insults to the embryo at different stages. Here, we examine how the glutathione system responds to the model pro-oxidants tert-butylhydroperoxide and tert-butylhydroquinone at different developmental stages, and the role of Nuclear factor erythroid 2-related factor (Nrf) proteins in regulating developmental glutathione redox status. Embryos became increasingly sensitive to pro-oxidants after 72h post-fertilization (hpf), after which the duration of the recovery period for the glutathione redox potential was increased. To determine whether the doubling of glutathione or the dynamic changes in glutathione redox potential are mediated by zebrafish paralogs of Nrf transcription factors, morpholino oligonucleotides were used to knock down translation of Nrf1 and Nrf2 (nrf1a, nrf1b, nrf2a, nrf2b). Knockdown of Nrf1a or Nrf1b perturbed glutathione redox state until 72 hpf. Knockdown of Nrf2 paralogs also perturbed glutathione redox state but did not significantly affect the response of glutathione to pro-oxidants. Nrf1b morphants had decreased gene expression of glutathione synthesis enzymes, while hsp70 increased in Nrf2b morphants. This work demonstrates that despite having a more robust glutathione system, embryos become more sensitive to oxidative stress later in development, and that neither Nrf1 nor Nrf2 alone appear to be essential for the response and recovery of glutathione to oxidative insults.

Hemigrapsus sanguineus, the Asian shore crab, has rapidly replaced Carcinus maenas, the green crab, as the most abundant crab on rocky shores in the northwest Atlantic since its introduction to the United States (USA) in 1988. The northern edge of this progressing invasion is the Gulf of Maine, where Asian shore crabs are only abundant in the south. We compared H. sanguineus population densities to those from published 2005 surveys and quantified genetic variation using the cytochrome c oxidase subunit I gene. We found that the range of H. sanguineus had extended northward since 2005, that population density had increased substantially (at least 10-fold at all sites), and that Asian shore crabs had become the dominant intertidal crab species in New Hampshire and southern Maine. Despite the significant increase in population density of H. sanguineus, populations only increased by a factor of 14 in Maine compared to 70 in southern New England, possibly due to cooler temperatures in the Gulf of Maine. Genetically, populations were predominantly composed of a single haplotype of Japanese, Korean, or Taiwanese origin, although an additional seven haplotypes were found. Six of these haplotypes were of Asian origin, while two are newly described. Large increases in population sizes of genetically diverse individuals in Maine will likely have a large ecological impact, causing a reduction in populations of mussels, barnacles, snails, and other crabs, similar to what has occurred at southern sites with large populations of this invasive crab species.

Development is a complex and well-defined process characterized by rapid cell proliferation and apoptosis. At this stage in life, a developmentally young organism is more sensitive to toxicants as compared to an adult. In response to pro-oxidant exposure, members of the Cap'n'Collar (CNC) basic leucine zipper (b-ZIP) transcription factor family (including Nfe2 and Nfe2-related factors, Nrfs) activate the expression of genes whose protein products contribute to reduced toxicity. Here, we studied the role of the CNC protein, Nfe2, in the developmental response to pro-oxidant exposure in the zebrafish (Danio rerio). Following acute waterborne exposures to diquat or tert-buytlhydroperoxide (tBOOH) at one of three developmental stages, wildtype (WT) and nfe2 knockout (KO) embryos and larvae were morphologically scored and their transcriptomes sequenced. Early in development, KO animals suffered from hypochromia that was made more severe through exposure to pro-oxidants; this phenotype in the KO may be linked to decreased expression of alas2, a gene involved in heme synthesis. WT and KO eleutheroembryos and larvae were phenotypically equally affected by exposure to pro-oxidants, where tBOOH caused more pronounced phenotypes as compared to diquat. Comparing diquat and tBOOH exposed embryos relative to the WT untreated control, a greater number of genes were up-regulated in the tBOOH condition as compared to diquat (tBOOH: 304 vs diquat: 148), including those commonly found to be differentially regulated in the vertebrate oxidative stress response (OSR) (e.g. hsp70.2, txn1, and gsr). When comparing WT and KO across all treatments and times, there were 1170 genes that were differentially expressed, of which 33 are known targets of the Nrf proteins Nrf1 and Nrf2. More specifically, in animals exposed to pro-oxidants a total of 968 genes were differentially expressed between WT and KO across developmental time, representing pathways involved in coagulation, embryonic organ development, body fluid level regulation, erythrocyte differentiation, and oxidation-reduction, amongst others. The greatest number of genes that changed in expression between WT and KO occurred in animals exposed to diquat at 2h post fertilization (hpf). Across time and treatment, there were six genes (dhx40, cfap70, dnajb9b, slc35f4, spi-c, and gpr19) that were significantly up-regulated in KO compared to WT and four genes (fhad1, cyp4v7, nlrp12, and slc16a6a) that were significantly down-regulated. None of these genes have been previously identified as targets of Nfe2 or the Nrf family. These results demonstrate that the zebrafish Nfe2 may be a regulator of both primitive erythropoiesis and the OSR during development.

Introduced over 200 yr ago to the east coast of North America, Carcinus maenas now ranges from New York to Newfoundland. In the 1980s, a secondary invasion of European lineages, termed northern haplotypes, occurred in Nova Scotia. Young-of-the-year sampled in 2007 revealed that northern haplotypes were present in low frequencies at several northwestern Atlantic sites as far south as New York; a model predicted an increase in their range and frequency over time. We collected samples in 2013 and 2014 to determine the haplotypes of adult crabs from New York to Nova Scotia. Six haplotypes, encompassing previously identified northern and southern haplotypes, 1 novel southern haplotype, and 1 Scandinavian haplotype, were identified in 275 crabs sampled at 11 sites. Northern haplotypes were only found in Nova Scotia, Beals Island (Maine), and Mount Desert Island (Maine) at a frequency of 60, 8, and 24%, respectively; remaining sites were predominantly composed of a previously identified southern haplotype. Northern haplotypes are limited in adult crabs to Mount Desert Island and north, indicating that the southern haplotype is selectively favored at some point during their life history, recruitment of northern larvae is limited south of Mount Desert Island, or entire year-classes post-2007 were lost. Our results do not support the predictions of an increase in the range and frequency of northern haplotypes, at least among adults, and indicate that a more complete knowledge of factors affecting C. maenas life stages is necessary to understand the current distribution of haplotypes.

Transcription factors in the CNC-bZIP family (NFE2, NRF1, NRF2 and NRF3) regulate genes with a wide range of functions in response to both physiological and exogenous signals, including those indicating changes in cellular redox status. Given their role in helping to maintain cellular homeostasis, it is imperative to understand the expression, regulation, and function of CNC-bZIP genes during embryonic development. We explored the expression and function of six nrf genes (nfe2, nrf1a, nrf1b, nrf2a, nrf2b, and nrf3) using zebrafish embryos as a model system. Analysis by microarray and quantitative RT-PCR showed that genes in the nrf family were expressed throughout development from oocytes to larvae. The spatial expression of nrf3 suggested a role in regulating the development of the brain, brachia and pectoral fins. Knock-down by morpholino anti-sense oligonucleotides suggested that none of the genes were necessary for embryonic viability, but nfe2 was required for proper cellular organization in the pneumatic duct and subsequent swim bladder function, as well as for proper formation of the otic vesicles. nrf genes were induced by the oxidant tert-butylhydroperoxide, and some of this response was regulated through family members Nrf2a and Nrf2b. Our results provide a foundation for understanding the role of nrf genes in normal development and in regulating the response to oxidative stress in vertebrate embryos.

A common geographical pattern of genetic variation is the one-dimensional cline. Clines may be maintained by diversifying selection across a geographical gradient but can also reflect historical processes such as allopatry followed by secondary contact. To identify loci that may be undergoing diversifying selection, we examined the distribution of geographical variation patterns across the range of the killifish (Fundulus heteroclitus) in 310 loci, including microsatellites, allozymes, and single nucleotide polymorphisms. We employed two approaches to detect loci under strong diversifying selection. First, we developed an automated method to identify clinal variation on a per-locus basis and examined the distribution of clines to detect those that exhibited signifcantly steeper slopes. Second, we employed a classic [Formula: see text]-outlier method as a complementary approach. We also assessed performance of these techniques using simulations. Overall, latitudinal clines were detected in nearly half of all loci genotyped (i.e., all eight microsatellite loci, 12 of 16 allozyme loci and 44% of the 285 SNPs). With the exception of few outlier loci (notably mtDNA and malate dehydrogenase), the positions and slopes of Fundulus clines were statistically indistinguishable. The high frequency of latitudinal clines across the genome indicates that secondary contact plays a central role in the historical demography of this species. Our simulation results indicate that accurately detecting diversifying selection using genome scans is extremely difficult in species with a strong signal of secondary contact; neutral evolution under this history produces clines as steep as those expected under selection. Based on these results, we propose that demographic history can explain all clinal patterns observed in F. heteroclitus without invoking natural selection to either establish or maintain the pattern we observe today.

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

Evolution by natural selection acts on natural populations amidst migration, gene-by-environmental interactions, constraints, and tradeoffs, which affect the rate and frequency of adaptive change. We asked how many and how rapidly loci change in populations subject to severe, recent environmental changes. To address these questions, we used genomic approaches to identify randomly selected single nucleotide polymorphisms (SNPs) with evolutionarily significant patterns in three natural populations of Fundulus heteroclitus that inhabit and have adapted to highly polluted Superfund sites. Three statistical tests identified 1.4-2.5% of SNPs that were significantly different from the neutral model in each polluted population. These nonneutral patterns in populations adapted to highly polluted environments suggest that these loci or closely linked loci are evolving by natural selection. One SNP identified in all polluted populations using all tests is in the gene for the xenobiotic metabolizing enzyme, cytochrome P4501A (CYP1A), which has been identified previously as being refractory to induction in the three highly polluted populations. Extrapolating across the genome, these data suggest that rapid evolutionary change in natural populations can involve hundreds of loci, a few of which will be shared in independent events.

Evolutionary solutions to the physiological challenges of life in highly variable habitats can span the continuum from evolution of a cosmopolitan plastic phenotype to the evolution of locally adapted phenotypes. Killifish (Fundulus sp.) have evolved both highly plastic and locally adapted phenotypes within different selective contexts, providing a comparative system in which to explore the genomic underpinnings of physiological plasticity and adaptive variation. Importantly, extensive variation exists among populations and species for tolerance to a variety of stressors, and we exploit this variation in comparative studies to yield insights into the genomic basis of evolved phenotypic variation. Notably, species of Fundulus occupy the continuum of osmotic habitats from freshwater to marine and populations within Fundulus heteroclitus span far greater variation in pollution tolerance than across all species of fish. Here, we explore how transcriptome regulation underpins extreme physiological plasticity on osmotic shock and how genomic and transcriptomic variation is associated with locally evolved pollution tolerance. We show that F. heteroclitus quickly acclimate to extreme osmotic shock by mounting a dramatic rapid transcriptomic response including an early crisis control phase followed by a tissue remodeling phase involving many regulatory pathways. We also show that convergent evolution of locally adapted pollution tolerance involves complex patterns of gene expression and genome sequence variation, which is confounded with body-weight dependence for some genes. Similarly, exploiting the natural phenotypic variation associated with other established and emerging model organisms is likely to greatly accelerate the pace of discovery of the genomic basis of phenotypic variation.

Background By targeting SNPs contained in both coding and non-coding areas of the genome, we are able to identify genetic differences and characterize genome-wide patterns of variation among individuals, populations and species. We investigated the utility of 454 sequencing and MassARRAY genotyping for population genetics in natural populations of the teleost, Fundulus heteroclitus as well as closely related Fundulus species ( F. grandis , F. majalis and F. similis ). Results We used 454 pyrosequencing and MassARRAY genotyping technology to identify and type 458 genome-wide SNPs and determine genetic differentiation within and between populations and species of Fundulus . Specifically, pyrosequencing identified 96 putative SNPs across coding and non-coding regions of the F. heteroclitus genome: 88.8% were verified as true SNPs with MassARRAY. Additionally, putative SNPs identified in F. heteroclitus EST sequences were verified in most (86.5%) F. heteroclitus individuals; fewer were genotyped in F. grandis (74.4%), F. majalis (72.9%), and F. similis (60.7%) individuals. SNPs were polymorphic and showed latitudinal clinal variation separating northern and southern populations and established isolation by distance in F. heteroclitus populations. In F. grandis , SNPs were less polymorphic but still established isolation by distance. Markers differentiated species and populations. Conclusions In total, these approaches were used to quickly determine differences within the Fundulus genome and provide markers for population genetic studies.

BACKGROUND:Populations of the teleost fish Fundulus heteroclitus appear to flourish in heavily polluted and geographically separated Superfund sites. Populations from three Superfund sites (New Bedford Harbor, MA, Newark Bay, NJ, and Elizabeth River, VA) have independently evolved adaptive resistance to chemical pollutants. In these polluted populations, natural selection likely has altered allele frequencies of loci that affect fitness or that are linked to these loci. The aim of this study was to identify loci that exhibit non-neutral behavior in the F. heteroclitus genome in polluted populations versus clean reference populations.RESULTS:To detect signatures of natural selection and thus identify genetic bases for adaptation to anthropogenic stressors, we examined allele frequencies for many hundreds of amplified fragment length polymorphism markers among populations of F. heteroclitus. Specifically, we contrasted populations from three Superfund sites (New Bedford Harbor, MA, Newark Bay, NJ, and Elizabeth River, VA) to clean reference populations flanking the polluted sites. When empirical FST values were compared to a simulated distribution of FST values, 24 distinct outlier loci were identified among pairwise comparisons of pollutant impacted F. heteroclitus populations and both surrounding reference populations. Upon removal of all outlier loci, there was a strong correlation (R2=0.79, p<0.0001) between genetic and geographical distance. This apparently neutral evolutionary pattern was not evident when outlier loci were included (R2=0.092, p=0.0721). Two outlier loci were shared between New Bedford Harbor and Elizabeth River populations, and two different loci were shared between Newark Bay and Elizabeth River populations.CONCLUSION:In total, 1% to 6% of loci are implicated as being under selection or linked to areas of the genome under selection in three F. heteroclitus populations that reside in polluted estuaries. Shared loci among polluted sites indicate that selection may be acting on multiple loci involved in adaptation, and loci shared between polluted sites potentially are involved in a generalized adaptive response.