Developmental exposure to pesticides that disrupt retinoic acid signaling causes persistent retinoid and behavioral dysfunction in zebrafish

Early developmental exposure to environmental toxicants may play a role in the risk for developing autism. A variety of pesticides have direct effects on retinoic acid (RA) signaling and as RA signaling has important roles in neurodevelopment, such compounds may cause developmental neurotoxicity through an overlapping adverse outcome pathway. It is hypothesized that a pesticide's embryonic effects on retinoid function may correspond with neurobehavioral disruption later in development. In the current studies, we determined the effects of RA-acting pesticides on neurobehavioral development in zebrafish. Buprofezin and imazalil caused generalized hypoactivity in the larval motility test, whereas chlorothalonil and endosulfan I led to selective hypoactivity and hyperactivity, respectively. With buprofezin, chlorothalonil, and imazalil, hypoactivity and/or novel anxiety-like behaviors persisted in adulthood and buprofezin additionally decreased social attraction responses in adulthood. Endosulfan I did not produce significant adult behavioral effects. Using qPCR analyses of adult brain tissue, we observed treatment-induced alterations in RA synthesis or catabolic genes, indicating persistent changes in RA homeostasis. These changes were compound-specific, with respect to expression directionality, and potential patterns of homeostatic disruption. Results suggest the likely persistence of disruptions in RA signaling well into adulthood and may represent compensatory mechanisms following early life stage exposures. This study demonstrates that early developmental exposure to environmental toxicants that interfere with RA signaling causes short as well as long-term behavioral disruption in a well-established zebrafish behavioral model and expand upon the meaning of the RA adverse outcome pathway, indicating that observed effects likely correspond with the nature of underlying homeostatic effects.

Legacy and emerging per- and polyfluoroalkyl substances suppress the neutrophil respiratory burst

Per- and polyfluoroalkyl substances (PFASs) are used in a multitude of processes and products, including nonstick coatings, food wrappers, and fire-fighting foams. These chemicals are environmentally-persistent, ubiquitous, and can be detected in the serum of 98% of Americans. Despite evidence that PFASs alter adaptive immunity, few studies have investigated their effects on innate immunity. The report here presents results of studies that investigated the impact of nine environmentally-relevant PFASs [e.g. perfluorooctanoic acid (PFOA), perfluorooctane sulfonic acid potassium salt (PFOS-K), perfluorononanoic acid (PFNA), perfluorohexanoic acid (PFHxA), perfluorohexane sulfonic acid (PFHxS), perfluorobutane sulfonic acid (PFBS), ammonium perfluoro(2-methyl-3-oxahexanoate) (GenX), 7H-perfluoro-4-methyl-3,6-dioxa-octane sulfonic acid (Nafion byproduct 2), and perfluoromethoxyacetic acid sodium salt (PFMOAA-Na)] on one component of the innate immune response, the neutrophil respiratory burst. The respiratory burst is a key innate immune process by which microbicidal reactive oxygen species (ROS) are rapidly induced by neutrophils in response to pathogens; defects in the respiratory burst can increase susceptibility to infection. The study here utilized larval zebrafish, a human neutrophil-like cell line, and primary human neutrophils to ascertain whether PFAS exposure inhibits ROS production in the respiratory burst. It was observed that exposure to PFHxA and GenX suppresses the respiratory burst in zebrafish larvae and a human neutrophil-like cell line. GenX also suppressed the respiratory burst in primary human neutrophils. This report is the first to demonstrate that these PFASs suppress neutrophil function and support the utility of employing zebrafish larvae and a human cell line as screening tools to identify chemicals that may suppress human immune function.

Comparative analysis of sucrose-embedding for whole-body zebrafish MSI by IR-MALDESI

Infrared matrix-assisted laser desorption electrospray ionization mass spectrometry imaging (IR-MALDESI) conventionally utilizes fresh-frozen biological tissues with an ice matrix to improve the detection of analytes. Sucrose-embedding with paraformaldehyde fixation has demonstrated feasibility as an alternative matrix for analysis by IR-MALDESI by preserving tissue features and enhancing ionization of lipids. However, investigating multi-organ systems provides broader context for a biological study and can elucidate more information about a disease state as opposed to a single organ. Danio rerio, or zebrafish, are model organisms for various disease states and can be imaged as a multi-organ sample to analyze morphological and metabolomic preservation as a result of sample preparation. Herein, whole-body zebrafish were imaged to compare sucrose-embedding with paraformaldehyde fixation against conventional fresh-frozen sample preparation. Serial sections were analyzed with and without an ice matrix to evaluate if sucrose functions as an alternative energy-absorbing matrix for IR-MALDESI applications across whole-body tissues. The resulting four conditions were compared in terms of total putative lipid annotations and category diversity, coverage across the entire m/z range, and ion abundance. Ultimately, sucrose-embedded zebrafish had an increase in putative lipid annotations for the combination of putative annotations with and without the application of an ice matrix relative to fresh-frozen tissues which require the application of an ice matrix. Upon the use of an ice matrix, a greater number of high mass putative lipid annotations (e.g., glycerophospholipids, glycerolipids, and sphingolipids) were identified. Conversely, without an ice matrix, sucrose-embedded sections elucidated more putative annotations in lower molecular weight lipids, including fatty acyls and sterol lipids. Similar to the mouse brain model, sucrose-embedding increased putative lipid annotation and abundance for whole-body zebrafish.

Maximized Spatial Information and Minimized Acquisition Time of Top-Hat IR-MALDESI-MSI of Zebrafish Using Nested Regions of Interest (nROIs)

Increasing the spatial resolution of a mass spectrometry imaging (MSI) method results in a more defined heatmap of the spatial distribution of molecules across a sample, but it is also associated with the disadvantage of increased acquisition time. Decreasing the area of the region of interest to achieve shorter durations results in the loss of potentially valuable information in larger specimens. This work presents a novel MSI method to reduce the time of MSI data acquisition with variable step size imaging: nested regions of interest (nROIs). Using nROIs, a small ROI may be imaged at a higher spatial resolution while nested inside a lower-spatial-resolution peripheral ROI. This conserves the maximal spatial and chemical information generated from target regions while also decreasing the necessary acquisition time. In this work, the nROI method was characterized on mouse liver and applied to top-hat MSI of zebrafish using a novel optical train, which resulted in a significant improvement in both acquisition time and spatial detail of the zebrafish. The nROI method can be employed with any step size pairing and adapted to any method in which the acquisition time of larger high-resolution ROIs poses a practical challenge.

TCDD alters essential transcriptional regulators of osteogenic differentiation in multipotent mesenchymal stem cells

Differentiation of multipotent mesenchymal stem cells (MSCs) into bone-forming osteoblasts requires strict coordination of transcriptional pathways. Aryl hydrocarbon receptor ligands, such as 2,3,7,8-tetrachlorodibenzo-p-dioxin (TCDD), have been shown to alter osteoblast differentiation in vitro and bone formation in multiple developmental in vivo models. The goal of the present study was to establish a global transcriptomic landscape during early, intermediate, and apical stages of osteogenic differentiation in vitro in response to TCDD exposure. Human bone-derived mesenchymal stem cells (hBMSCs) were cultured in growth media (GM), osteogenic differentiation media (ODM), or ODM containing 10 nM TCDD (ODM + TCDD), thus enabling a comparison of the transcriptomic profiles of undifferentiated, differentiated, and differentiated-TCDD-exposed hBMSCs, respectively. In this test system, exposure to TCDD attenuated the differentiation of hBMSCs into osteoblasts as evidenced by reduced alkaline phosphatase activity and mineralization. At various timepoints, we observed altered expression of genes that play a role in the Wnt, fibroblast growth factor, bone morphogenetic protein/transforming growth factor beta developmental pathways, as well as pathways related to extracellular matrix organization and deposition. Reconstruction of gene regulatory networks with the interactive dynamic regulatory event miner (iDREM) analysis revealed modulation of transcription factors (TFs) including POLR3G, NR4A1, RDBP, GTF2B, POU2F2, and ZEB1, which may putatively influence osteoblast differentiation and the requisite deposition and mineralization of bone extracellular matrix. We demonstrate that the combination of RNA-Seq data in conjunction with the iDREM regulatory model captures the transcriptional dynamics underlying MSC differentiation under different conditions in vitro. Model predictions are consistent with existing knowledge and provide a new tool to identify novel pathways and TFs that may facilitate a better understanding of the osteoblast differentiation process, perturbation by exogenous agents, and potential intervention strategies targeting those specific pathways.

Nonylphenol Polyethoxylates Enhance Adipose Deposition in Developmentally Exposed Zebrafish

Alkylphenol polyethoxylates (APEOs), such as nonylphenol ethoxylates (NPEOs), are high-production-volume surfactants used in laundry detergents, hard-surface cleaners, pesticide formulations, textile production, oils, paints, and other products. NPEOs comprise −80% of the total production of APEOs and are widely reported across diverse environmental matrices. Despite a growing push for replacement products, APEOs continue to be released into the environment through wastewater at significant levels. Research into related nonionic surfactants from varying sources has reported metabolic health impacts, and we have previously demonstrated that diverse APEOs and alcohol polyethoxylates promote adipogenesis in the murine 3T3-L1 pre-adipocyte model. These effects appeared to be independent of the base alkylphenol and related to the ethoxylate chain length, though limited research has evaluated NPEO exposures in animal models. The goals of this study were to assess the potential of NPEOs to promote adiposity (Nile red fluorescence quantification) and alter growth and/or development (toxicity, length, weight, and energy expenditure) of developmentally exposed zebrafish (Danio rerio). We also sought to expand our understanding of the ability to promote adiposity through evaluation in human mesenchymal stem cells. Herein, we demonstrated consistent adipogenic effects in two separate human bone-marrow-derived mesenchymal stem cell models, and that nonylphenol and its ethoxylates promoted weight gain and increased adipose deposition in developmentally exposed zebrafish. Notably, across both cell and zebrafish models we report increasing adipogenic/obesogenic activity with increasing ethoxylate chain lengths up to maximums around NPEO-6 and then decreasing activity with the longest ethoxylate chain lengths. This research suggests metabolic health concerns for these common obesogens, suggesting further need to assess molecular mechanisms and better characterize environmental concentrations for human health risk assessments.

Vitamin D deficiency promotes accumulation of bioactive lipids and increased endocannabinoid tone in zebrafish

Vitamin D is well known for its traditional role in bone mineral homeostasis; however, recent evidence suggests that vitamin D also plays a significant role in metabolic control. This study served to investigate putative linkages between vitamin D deficiency (VDD) and metabolic disruption of bioactive lipids by MS imaging. Our approach employed infrared-matrix-assisted laser desorption electrospray ionization MS imaging for lipid metabolite profiling in 6-month-old zebrafish fed either a VDD or a vitamin D-sufficient (VDS) diet. Using a lipidomics pipeline, we found that VDD zebrafish had a greater abundance of bioactive lipids (N-acyls, endocannabinoids [ECs], diacylglycerols/triacylglycerols, bile acids/bile alcohols, and vitamin D derivatives) suggestive of increased EC tone compared with VDS zebrafish. Tandem MS was performed on several differentially expressed metabolites with sufficient ion abundances to aid in structural elucidation and provide additional support for MS annotations. To confirm activation of the EC pathways, we subsequently examined expression of genes involved in EC biosynthesis, metabolism, and receptor signaling in adipose tissue and liver from VDD and VDS zebrafish. Gene expression changes were congruent with increased EC tone, with VDD zebrafish demonstrating increased synthesis and metabolism of anandamide compared with VDS zebrafish. Taken together, our data suggest that VDD may promote accumulation of bioactive lipids and increased EC tone in zebrafish.

Thyroid Receptor Antagonism of Chemicals Extracted from Personal Silicone Wristbands within a Papillary Thyroid Cancer Pilot Study

Research suggests that thyroid cancer incidence rates are increasing, and environmental exposures have been postulated to be playing a role. To explore this possibility, we conducted a pilot study to investigate the thyroid disrupting bioactivity of chemical mixtures isolated from personal silicone wristband samplers within a thyroid cancer cohort. Specifically, we evaluated TRβ antagonism of chemical mixtures extracted from wristbands (n = 72) worn by adults in central North Carolina participating in a case-control study on papillary thyroid cancer. Sections of wristbands were solvent-extracted and analyzed via mass spectrometry to quantify a suite of semivolatile chemicals. A second extract from each wristband was used in a bioassay to quantify TRβ antagonism in human embryonic kidney cells (HEK293/17) at concentrations ranging from 0.1 to 10% of the original extract (by volume). Approximately 70% of the sample extracts tested at a 1% extract concentration exhibited significant TRβ antagonism, with a mean of 30% and a range of 0-100%. Inhibited cell viability was noted in >20% of samples that were tested at 5 and 10% concentrations. Antagonism was positively associated with wristband concentrations of several phthalates, organophosphate esters, and brominated flame retardants. These results suggest that personal passive samplers may be useful in evaluating the bioactivities of mixtures that people contact on a daily basis. We also report tentative associations between thyroid receptor antagonism, chemical concentrations, and papillary thyroid cancer case status. Future research utilizing larger sample sizes, prospective data collection, and measurement of serum thyroid hormone levels (which were not possible in this study) should be utilized to more comprehensively evaluate these associations.

Zebrafish show long-term behavioral impairments resulting from developmental vitamin D deficiency

Vitamin D has been shown in a wide variety of species to play critical roles in neurodevelopment. Vitamin D deficiency disrupts development of the brain and can cause lasting behavioral dysfunction. Zebrafish have become an important model for the study of development in general and neurodevelopment in particular. Zebrafish were used in the current study to characterize the effects of developmental vitamin D deficiency on behavioral function. Adult zebrafish that had been chronically fed a vitamin D deficient or replete diets were bred and the offspring were continued on those diets. The offspring were behaviorally tested as adults. In the novel tank diving test the vitamin D deficient diet significantly lowered the vertical position of fish indicative of more anxiety-like behavior. In the novel tank diving test swimming activity was also significantly decreased by vitamin D deficiency. Startle response was increased by developmental vitamin D deficiency during the early part of the test. No significant effects of vitamin D deficiency were seen with social affiliation and predatory stimulus avoidance tests. These results indicate a phenotype of vitamin D deficiency characterized by more anxiety-like behavior. This result was relatively specific inasmuch as few or no behavioral effects were seen in other behavioral tests.

Developmental exposure of zebrafish to vitamin D receptor acting drugs and environmental toxicants disrupts behavioral function

Vitamin D receptor (VDR) signaling is important for optimal neurobehavioral development. Disruption of VDR signaling by environmental toxicants during early development might contribute to the etiology of behavioral dysfunction. In the current set of studies, we examined ten compounds known to affect VDR function in vitro for neurobehavioral effects in vivo in zebrafish. Zebrafish embryos were exposed to concentrations of the compounds in their water during the first 5 days post-fertilization. On day 5, the embryos were tested in an alternating light-dark locomotor assay using a computerized video tracking system. We found that most of the compounds produced significant changes in locomotor behavior in exposed zebrafish larvae, although the direction of the effect (i.e., hypo- or hyperactivity) and the sensitivity of the effect to changes in illumination condition varied across the compounds. The nature of the behavioral effects generally corresponded to the effects these compounds have been shown to exert on VDR. These studies lay a foundation for further investigation to determine whether behavioral dysfunction persists into adulthood and if so which behavioral functions are affected. Zebrafish can be useful for screening compounds identified in high throughput in vitro assays to provide an initial test for how those compounds would affect construction and behavioral function of a complex nervous system, helping to bridge the gap between in vitro neurotoxicity assays and mammalian models for risk assessment in humans.

Methods for Cryosectioning and Mass Spectrometry Imaging of Whole-Body Zebrafish

The zebrafish (Danio rerio) is an ideal model for whole animal studies of lipid metabolism and lipid-related disease. In this work, infrared matrix-assisted laser desorption electrospray ionization (IR-MALDESI) mass spectrometry imaging (MSI) was applied for direct visualization of lipid and metabolite distributions across various organs in whole-body zebrafish tissue sections. Detailed methods for overcoming the challenges of cryosectioning adult male zebrafish for MSI and complementary histological imaging are described. Representative two-dimensional ion maps demonstrated organ specific localization of lipid analytes allowing for visualization of areas of interest including the brain, liver, intestines, and skeletal muscle. A high resolving power mass spectrometer was utilized for accurate mass measurements, which permitted the use of open-source, web-based tools for MS¹ annotations including METASPACE and METLIN. Whole-body MSI with IR-MALDESI allowed for broad lipid coverage with high spatial resolution, illustrating the potential of this technique for studying lipid-related diseases using zebrafish as a model organism.

Maternal vitamin D deficiency and developmental origins of health and disease (DOHaD)

Vitamin D is an essential nutrient that is metabolized in the body to generate an active metabolite (1,25(OH)2D) with hormone-like activity and highly diverse roles in cellular function. Vitamin D deficiency (VDD) is a prevalent but easily preventable nutritional disturbance. Emerging evidence demonstrates the importance of sufficient vitamin D concentrations during fetal life with deficiencies leading to long-term effects into adulthood. Here, we provide a detailed review and perspective of evidence for the role of maternal VDD in offspring long term health, particularly as it relates to Developmental Origins of Health and Disease (DOHaD). We focus on roles in neurobehavioral and cardiometabolic disorders in humans and highlight recent findings from zebrafish and rodent models that probe potential mechanisms linking early life VDD to later life health outcomes. Moreover, we explore evidence implicating epigenetic mechanisms as a mediator of this link. Gaps in our current understanding of how maternal VDD might result in deleterious offspring outcomes later in life are also addressed.

Evidence for Aryl hydrocarbon Receptor-Mediated Inhibition of Osteoblast Differentiation in Human Mesenchymal Stem Cells

Multipotent mesenchymal stem cells (MSCs) maintain the ability to differentiate into adipogenic, chondrogenic, or osteogenic cell lineages. There is increasing concern that exposure to environmental agents such as aryl hydrocarbon receptor (AhR) ligands, may perturb the osteogenic pathways responsible for normal bone formation. The objective of the current study was to evaluate the potential of the prototypic AhR ligand 2,3,7,8-tetrachlorodibenzo-p-dioxin (TCDD) to disrupt osteogenic differentiation of human bone-derived MSCs (hBMSCs) in vitro. Primary hBMSCs from three donors were exposed to 10 nM TCDD and differentiation was interrogated using select histological, biochemical, and transcriptional markers of osteogenesis. Exposure to 10 nM TCDD resulted in an overall consistent attenuation of alkaline phosphatase (ALP) activity and matrix mineralization at terminal stages of differentiation in primary hBMSCs. At the transcriptional level, the transcriptional regulator DLX5 and additional osteogenic markers (ALP, OPN, and IBSP) displayed attenuated expression; conversely, FGF9 and FGF18 were consistently upregulated in each donor. Expression of stem cell potency markers SOX2, NANOG, and SALL4 decreased in the osteogenic controls, whereas expression in TCDD-treated cells resembled that of undifferentiated cells. Coexposure with the AhR antagonist GNF351 blocked TCDD-mediated attenuation of matrix mineralization, and either fully or partially rescued expression of genes associated with osteogenic regulation, extracellular matrix, and/or maintenance of multipotency. Thus, experimental evidence from this study suggests that AhR transactivation likely attenuates osteoblast differentiation in multipotent hBMSCs. This study also underscores the use of primary human MSCs to evaluate osteoinductive or osteotoxic potential of chemical and pharmacologic agents in vitro.

Relation of contaminants to fish intersex in riverine sport fishes

Endocrine active compounds (EACs) are pollutants that have been recognized as an emerging and widespread threat to aquatic ecosystems globally. Intersex, the presence of female germ cells within a predominantly male gonad, is considered a biomarker of endocrine disruption caused by EACs. We measured a suite of EACs and assessed their associated impacts on fish intersex occurrence and severity in a large, regulated river system in North Carolina and South Carolina, USA. Our specific objective was to determine the relationship of contaminants in water, sediment, and fish tissue with the occurrence and severity of the intersex condition in wild, adult black bass (Micropterus), sunfish (Lepomis), and catfish (Ictaluridae) species at 11 sites located on the Yadkin-Pee Dee River. Polycyclic aromatic hydrocarbons (PAHs), ethinylestradiol (EE2), and heavy metals were the most prevalent contaminants that exceeded effect levels for the protection of aquatic organisms. Fish intersex condition was most frequently observed and most severe in black basses and was less frequently detected and less severe in sunfishes and catfishes. The occurrence of the intersex condition in fish showed site-related effects, rather than increasing longitudinal trends from upstream to downstream. Mean black bass and catfish tissue contaminant concentrations were higher than that of sunfish, likely because of the latter's lower trophic position in the food web. Principal component analysis identified waterborne PAHs as the most correlated environmental contaminant with intersex occurrence and severity in black bass and sunfish. As indicated by the intersex condition, EACs have adverse but often variable effects on the health of wild sport fishes in this river, likely due to fluctuations in EAC inputs and the dynamic nature of the riverine system. These findings enhance the understanding of the relationship between contaminants and fish health and provide information to guide ecologically comprehensive conservation and management decisions.

An object-oriented Bayesian network approach for establishing swine manure-borne natural estrogenic compounds budget

A facility-wide estrogen budget model was developed to assess the excretion of natural estrogens by swine in a commercial swine farrowing concentrated animal feeding operations (CAFO) in North Carolina, using an object-oriented Bayesian network (OOBN) approach. The OOBN model is the combination of twelve objects of Bayesian network models, which characterize the estrogen budget flows based on the sow reproductive cycle (i.e., pre-estrus, estrus, and lactation) for the three natural estrogen types [estrone (E1), 17β-estradiol (E2), and estriol (E3)] within each barn. This OOBN model provides a mechanism to quantify the levels of the natural estrogens and their probabilistic distributions with regard to estrogen type, waste sources such as urine, feces, and recycling lagoon slurry, and animal reproductive status. Moreover, the OOBN model allows us to assess the overall contribution of natural estrogen compounds from each operational unit of the CAFO, while accounting for the uncertainties. Results from the OOBN model indicate a rank order of lactating sows > gestating sows > breeding sows in terms of contribution of estrogen loads to the total natural estrogen budget. As to estrogen type, E1 was found as the major estrogen metabolite with the summed concentrations of urine, feces, and flushing slurry wastes exceeding 3000 ng/l > 90% of the time. As to waste sources, the flushing slurry waste was found to be a major contributor of the estrogen budget compared with urine and feces wastes from barn animals.

Confirmation of high-throughput screening data and novel mechanistic insights into VDR-xenobiotic interactions by orthogonal assays

High throughput screening (HTS) programs have demonstrated that the Vitamin D receptor (VDR) is activated and/or antagonized by a wide range of structurally diverse chemicals. In this study, we examined the Tox21 qHTS data set generated against VDR for reproducibility and concordance and elucidated functional insights into VDR-xenobiotic interactions. Twenty-one potential VDR agonists and 19 VDR antagonists were identified from a subset of >400 compounds with putative VDR activity and examined for VDR functionality utilizing select orthogonal assays. Transient transactivation assay (TT) using a human VDR plasmid and Cyp24 luciferase reporter construct revealed 20/21 active VDR agonists and 18/19 active VDR antagonists. Mammalian-2-hybrid assay (M2H) was then used to evaluate VDR interactions with co-activators and co-regulators. With the exception of a select few compounds, VDR agonists exhibited significant recruitment of co-regulators and co-activators whereas antagonists exhibited considerable attenuation of recruitment by VDR. A unique set of compounds exhibiting synergistic activity in antagonist mode and no activity in agonist mode was identified. Cheminformatics modeling of VDR-ligand interactions were conducted and revealed selective ligand VDR interaction. Overall, data emphasizes the molecular complexity of ligand-mediated interactions with VDR and suggest that VDR transactivation may be a target site of action for diverse xenobiotics.

Windows of Susceptibility and Consequences of Early Life Exposures to 17β-estradiol on Medaka (Oryzias latipes) Reproductive Success

Estrogens and estrogen mimics are commonly found in surface waters and are associated with deleterious effects in fish populations. Impaired fertility and fecundity in fish following chronic exposures to estrogens and estrogen mimics during critical windows in development are well documented. However, information regarding differential reproductive effects of exposure within defined developmental stages remains sparse. In this study, reproductive capacity was assessed in Japanese medaka (Oryzias latipes) after exposure to two concentrations of 17β-estradiol (E2β; 2 ng/L and 50 ng/L) during four distinct stages of development: gonad development, gonad differentiation, development of secondary sex characteristics (SSC) and gametogenesis. Exposure to E2β did not adversely impact survival, hatch success, growth, or genotypic ratios. In contrast, exposure to 50 ng/L E2β during SSC development altered phenotypic ratios and SSC. Exposure to both E2β treatments reduced reproductive capacity (fertility, fecundity) by 7.3-57.4% in adult medaka breeding pairs, with hindrance of SSC development resulting in the largest disruption in breeding capacity (51.6-57.4% decrease) in the high concentration. This study documents differential effects among four critical stages of development and provides insight into factors (window of exposure, exposure concentration and duration of exposure period) contributing to reproductive disruption in fish.

Endocrine active contaminants in aquatic systems and intersex in common sport fishes

Male fish are susceptible to developing intersex, a condition characterized by the presence of testicular oocytes. In the present study, the relationship between intersex and exposure to estrogenic endocrine active contaminants (EACs) was assessed for 2 genera of sport fish, Micropterus and Lepomis, at 20 riverine sites. Seasonal trends and relationships between EACs and intersex (prevalence and severity) were examined at varying putative sources of EACs throughout North Carolina, identified as point sources, nonpoint sources, and reference sites. Intersex was identified in both genera, which was documented for the first time in wild-caught Lepomis. Intersex was more prevalent (59.8%) and more severe (1.6 mean rank) in Micropterus, which was highly correlation to EACs in sediment. In contrast, intersex was less common (9.9%) and less severe (0.2 mean rank) in Lepomis and was highly correlated to EACs in the water column. The authors found that concentrations of polycyclic aromatic hydrocarbons, polychlorinated biphenyls, industrial EACs, and estrogens were highest at point source sites; however, no source type variation was identified in the prevalence or severity of intersex, nor were there seasonal trends in intersex or EAC concentrations. The authors' results associate genus-specific prevalence of intersex with specific EAC classes in common sport fishes having biological, ecological, and conservation implications. Environ Toxicol Chem 2017;36:959-968. © 2016 SETAC.

Exposure to a PBDE/OH-BDE mixture alters juvenile zebrafish (Danio rerio) development

Polybrominated diphenyl ethers (PBDEs) and their metabolites (e.g., hydroxylated BDEs [OH-BDEs]) are contaminants frequently detected together in human tissues and are structurally similar to thyroid hormones. Thyroid hormones partially mediate metamorphic transitions between life stages in zebrafish, making this a critical developmental window that may be vulnerable to chemicals disrupting thyroid signaling. In the present study, zebrafish were exposed to 6-OH-BDE-47 (30 nM; 15 μg/L) alone, or to a low-dose (30 μg/L) or high-dose (600 μg/L) mixture of PentaBDEs, 6-OH-BDE-47 (0.5-6 μg/L), and 2,4,6-tribromophenol (5-100 μg/L) during juvenile development (9-23 d postfertilization) and evaluated for developmental endpoints mediated by thyroid hormone signaling. Fish were sampled at 3 time points and examined for developmental and skeletal morphology, apical thyroid and skeletal gene markers, and modifications in swimming behavior (as adults). Exposure to the high-dose mixture resulted in >85% mortality within 1 wk of exposure, despite being below reported acute toxicity thresholds for individual congeners. The low-dose mixture and 6-OH-BDE-47 groups exhibited reductions in body length and delayed maturation, specifically relating to swim bladder, fin, and pigmentation development. Reduced skeletal ossification was also observed in 6-OH-BDE-47-treated fish. Assessment of thyroid and osteochondral gene regulatory networks demonstrated significantly increased expression of genes that regulate skeletal development and thyroid hormones. Overall, these results indicate that exposures to PBDE/OH-BDE mixtures adversely impact zebrafish maturation during metamorphosis. Environ Toxicol Chem 2017;36:36-48. © 2016 SETAC.

From the Cover: Embryonic Exposure to TCDD Impacts Osteogenesis of the Axial Skeleton in Japanese medaka, Oryzias latipes

Recent studies from mammalian, fish, and in vitro models have identified bone and cartilage development as sensitive targets for dioxins and other aryl hydrocarbon receptor ligands. In this study, we assess how embryonic 2,3,7,8-tetrachlorochlorodibenzo-p-dioxin (TCDD) exposure impacts axial osteogenesis in Japanese medaka (Oryzias latipes), a vertebrate model of human bone development. Embryos from inbred wild-type Orange-red Hd-dR and 3 transgenic medaka lines (twist:EGFP, osx/sp7:mCherry, col10a1:nlGFP) were exposed to 0.15 nM and 0.3 nM TCDD and reared until 20 dpf. Individuals were stained for mineralized bone and imaged using confocal microscopy to assess skeletal alterations in medial vertebrae in combination with a qualitative spatial analysis of osteoblast and osteoblast progenitor cell populations. Exposure to TCDD resulted in an overall attenuation of vertebral ossification characterized by truncated centra, and reduced neural and hemal arch lengths. Effects on mineralization were consistent with modifications in cell number and cell localization of transgene-labeled osteoblast and osteoblast progenitor cells. Endogenous expression of osteogenic regulators runt-related transcription factor 2 (runx2) and osterix (osx/sp7), and extracellular matrix genes osteopontin (spp1), collagen type I alpha I (col1), collagen type X alpha I (col10a1), and osteocalcin (bglap/osc) was significantly diminished at 20 dpf following TCDD exposure as compared with controls. Through global transcriptomic analysis more than 590 differentially expressed genes were identified and mapped to select pathological states including inflammatory disease, connective tissue disorders, and skeletal and muscular disorders. Taken together, results from this study suggest that TCDD exposure inhibits axial bone formation through dysregulation of osteoblast differentiation. This approach highlights the advantages and sensitivity of using small fish models to investigate how xenobiotic exposure may impact skeletal development.

Sharing the Roles: An Assessment of Japanese Medaka Estrogen Receptors in Vitellogenin Induction

Teleost fish express at least three estrogen receptor (ER) subtypes. To date, however, the individual role of these ER subtypes in regulating expression of estrogen responsive genes remains ambiguous. Here, we investigate putative roles of three ER subtypes in Japanese medaka (Oryzias latipes), using vitellogenin (VTG) I and II as model genes. We identify specific ligand/receptor/promoter dynamics, using transient transactivation assays that incorporate luciferase reporters comprising 3kb promoter/enhancer regions of medaka VTGI and VTGII genes. Four steroidal estrogens (17β-estradiol, estrone, estriol, and 17α-estradiol) were tested in these assays. Results indicate that all three medaka ERs (mERs) are capable of initiating transactivation of both VTG I and II, with ERβ2 exhibiting greatest activity. Promoter deletion analysis suggests that ligand-specific receptor transactivation and utilization of regional-specific estrogen response elements may be associated with differential activities of each medaka ER. Further, cluster analysis of in vivo gene expression and in vitro transactivation suggests that all three ER subtypes putatively play a role in up-regulation of VTG. Results illustrate that preferential ligand/receptor/promoter interactions may have direct implications for VTG gene expression and other ER-mediated regulatory functions that are relevant to the risk assessment of estrogenic compounds.

Advancing toxicology research using in vivo high throughput toxicology with small fish models

Small freshwater fish models, especially zebrafish, offer advantages over traditional rodent models, including low maintenance and husbandry costs, high fecundity, genetic diversity, physiology similar to that of traditional biomedical models, and reduced animal welfare concerns. The Collaborative Workshop on Aquatic Models and 21st Century Toxicology was held at North Carolina State University on May 5-6, 2014, in Raleigh, North Carolina, USA. Participants discussed the ways in which small fish are being used as models to screen toxicants and understand mechanisms of toxicity. Workshop participants agreed that the lack of standardized protocols is an impediment to broader acceptance of these models, whereas development of standardized protocols, validation, and subsequent regulatory acceptance would facilitate greater usage. Given the advantages and increasing application of small fish models, there was widespread interest in follow-up workshops to review and discuss developments in their use. In this article, we summarize the recommendations formulated by workshop participants to enhance the utility of small fish species in toxicology studies, as well as many of the advances in the field of toxicology that resulted from using small fish species, including advances in developmental toxicology, cardiovascular toxicology, neurotoxicology, and immunotoxicology. We alsoreview many emerging issues that will benefit from using small fish species, especially zebrafish, and new technologies that will enable using these organisms to yield results unprecedented in their information content to better understand how toxicants affect development and health.

Predicting characteristics of rainfall driven estrogen runoff and transport from swine AFO spray fields

Animal feeding operations (AFOs) have been implicated as potentially major sources of estrogenic contaminants into the aquatic environment due to the relatively minimal treatment of waste and potential mobilization and transport of waste components from spray fields. In this study a Bayesian network (BN) model was developed to inform management decisions and better predict the transport and fate of natural steroidal estrogens from these sites. The developed BN model integrates processes of surface runoff and sediment loss with the modified universal soil loss equation (MUSLE) and the soil conservation service curve number (SCS-CN) runoff model. What-if scenario simulations of lagoon slurry wastes to the spray fields were conducted for the most abundant natural estrogen estrone (E1) observed in the system. It was found that E1 attenuated significantly after 2 months following waste slurry application in both spring and summer seasons, with the overall attenuation rate predicted to be higher in the summer compared to the spring. Using simulations of rainfall events in conjunction with waste slurry application rates, it was predicted that the magnitude of E1 runoff loss is significantly higher in the spring as compared to the summer months, primarily due to spray field crop management plans. Our what-if scenario analyses suggest that planting Bermuda grass in the spray fields is likely to reduce runoff losses of natural estrogens near the water bodies and ecosystems, as compared to planting of soybeans.

Molecular cloning, functional characterization, and evolutionary analysis of vitamin D receptors isolated from basal vertebrates

The vertebrate genome is a result of two rapid and successive rounds of whole genome duplication, referred to as 1R and 2R. Furthermore, teleost fish have undergone a third whole genome duplication (3R) specific to their lineage, resulting in the retention of multiple gene paralogs. The more recent 3R event in teleosts provides a unique opportunity to gain insight into how genes evolve through specific evolutionary processes. In this study we compare molecular activities of vitamin D receptors (VDR) from basal species that diverged at key points in vertebrate evolution in order to infer derived and ancestral VDR functions of teleost paralogs. Species include the sea lamprey (Petromyzon marinus), a 1R jawless fish; the little skate (Leucoraja erinacea), a cartilaginous fish that diverged after the 2R event; and the Senegal bichir (Polypterus senegalus), a primitive 2R ray-finned fish. Saturation binding assays and gel mobility shift assays demonstrate high affinity ligand binding and classic DNA binding characteristics of VDR has been conserved across vertebrate evolution. Concentration response curves in transient transfection assays reveal EC50 values in the low nanomolar range, however maximum transactivational efficacy varies significantly between receptor orthologs. Protein-protein interactions were investigated using co-transfection, mammalian 2-hybrid assays, and mutations of coregulator activation domains. We then combined these results with our previous study of VDR paralogs from 3R teleosts into a bioinformatics analysis. Our results suggest that 1, 25D3 acts as a partial agonist in basal species. Furthermore, our bioinformatics analysis suggests that functional differences between VDR orthologs and paralogs are influenced by differential protein interactions with essential coregulator proteins. We speculate that we may be observing a change in the pharmacodynamics relationship between VDR and 1, 25D3 throughout vertebrate evolution that may have been driven by changes in protein-protein interactions between VDR and essential coregulators.

Functional diversification of vitamin D receptor paralogs in teleost fish after a whole genome duplication event

The diversity and success of teleost fishes (Actinopterygii) has been attributed to three successive rounds of whole-genome duplication (WGD). WGDs provide a source of raw genetic material for evolutionary forces to act upon, resulting in the divergence of genes with altered or novel functions. The retention of multiple gene pairs (paralogs) in teleosts provides a unique opportunity to study how genes diversify and evolve after a WGD. This study examines the hypothesis that vitamin D receptor (VDR) paralogs (VDRα and VDRβ) from two distantly related teleost orders have undergone functional divergence subsequent to the teleost-specific WGD. VDRα and VDRβ paralogs were cloned from the Japanese medaka (Beloniformes) and the zebrafish (Cypriniformes). Initial transactivation studies using 1α, 25-dihydroxyvitamin D3 revealed that although VDRα and VDRβ maintain similar ligand potency, the maximum efficacy of VDRβ was significantly attenuated compared with VDRα in both species. Subsequent analyses revealed that VDRα and VDRβ maintain highly similar ligand affinities; however, VDRα demonstrated preferential DNA binding compared with VDRβ. Protein-protein interactions between the VDR paralogs and essential nuclear receptor coactivators were investigated using transactivation and mammalian two-hybrid assays. Our results imply that functional differences between VDRα and VDRβ occurred early in teleost evolution because they are conserved between distantly related species. Our results further suggest that the observed differences may be associated with differential protein-protein interactions between the VDR paralogs and coactivators. We speculate that the observed functional differences are due to subtle ligand-induced conformational differences between the two paralogs, leading to divergent downstream functions.

Transport of steroid hormones, phytoestrogens, and estrogenic activity across a swine lagoon/sprayfield system

The inflow, transformation, and attenuation of natural steroid hormones and phytoestrogens and estrogenic activity were assessed across the lagoon/sprayfield system of a prototypical commercial swine sow operation. Free and conjugated steroid hormones (estrogens, androgens, and progesterone) were detected in urine and feces of sows across reproductive stages, with progesterone being the most abundant steroid hormone. Excreta also contained phytoestrogens indicative of a soy-based diet, particularly, daidzein, genistein, and equol. During storage in barn pits and the anaerobic lagoon, conjugated hormones dissipated, and androgens and progesterone were attenuated. Estrone and equol persisted along the waste disposal route. Following application of lagoon slurry to agricultural soils, all analytes exhibited attenuation within 2 days. However, analytes including estrone, androstenedione, progesterone, and equol remained detectable in soil at 2 months postapplication. Estrogenic activity in the yeast estrogen screen and T47D-KBluc in vitro bioassays generally tracked well with analyte concentrations. Estrone was found to be the greatest contributor to estrogenic activity across all sample types. This investigation encompasses the most comprehensive suite of natural hormone and phytoestrogen analytes examined to date across a livestock lagoon/sprayfield and provides global insight into the fate of these analytes in this widely used waste management system.

A Bayesian network model for assessing natural estrogen fate and transport in a swine waste lagoon

Commercial swine waste lagoons are regarded as a major reservoir of natural estrogens, which have the potential to produce adverse physiological effects on exposed aquatic organisms and wildlife. However, there remains limited understanding of the complex mechanisms of physical, chemical, and biological processes that govern the fate and transport of natural estrogens within an anaerobic swine lagoon. To improve lagoon management and ultimately help control the offsite transport of these compounds from swine operations, a probabilistic Bayesian network model was developed to assess natural estrogen fate and budget and then compared against data collected from a commercial swine field site. In general, the model was able to describe the estrogen fate and budget in both the slurry and sludge stores within the swine lagoon. Sensitivity analysis within the model demonstrated that the estrogen input loading from the associated barn facility was the most important factor in controlling estrogen concentrations within the lagoon slurry storage, whereas the settling rate was the most significant factor in the lagoon sludge storage. The degradation reactions were shown to be minor in both stores based on prediction of average total estrogen concentrations. Management scenario evaluations demonstrated that the best possible management options to reduce estrogen levels in the lagoon are either to adjust the estrogen input loading from swine barn facilities or to effectively enhance estrogen bonding with suspended solids through the use of organic polymers or inorganic coagulants.

Bridging the gap from screening assays to estrogenic effects in fish: potential roles of multiple estrogen receptor subtypes

This study seeks to delineate the ligand interactions that drive biomarker induction in fish exposed to estrogenic pollutants and provide a case study on the capacity of human (h) estrogen receptor (ER)-based in vitro screening assays to predict estrogenic effects in aquatic species. Adult male Japanese medaka (Oryzias latipes) were exposed to solutions of singular steroidal estrogens or to the estrogenic extract of an anaerobic swine waste lagoon. All exposure concentrations were calibrated to be equipotent based on the yeast estrogen screen (YES), which reports activation of hERα. These exposures elicited significantly different magnitudes of hepatic vitellogenin and choriogenin gene induction in the male medaka. Effects of the same YES-calibrated solutions in the T47D-KBluc assay, which reports activation of hERα and hERβ, generally recapitulated observations in medaka. Using competitive ligand binding assays, it was found that the magnitude of vitellogenin/choriogenin induction by different estrogenic ligands correlated positively with preferential binding affinity for medaka ERβ subtypes, which are highly expressed in male medaka liver prior to estrogen exposure. Results support emerging evidence that ERβ subtypes are critically involved in the teleost estrogenic response, with the ERα:ERβ ratio being of particular importance. Accordingly, incorporation of multiple ER subtypes into estrogen screening protocols may increase predictive value for the risk assessment of aquatic systems, including complex estrogenic mixtures.

Comprehensive assessment of hormones, phytoestrogens, and estrogenic activity in an anaerobic swine waste lagoon

In this study, the distribution of steroid hormones, phytoestrogens, and estrogenic activity was thoroughly characterized within the anaerobic waste lagoon of a typical commercial swine sow operation. Three independent rounds of sampling were conducted in June 2009, April 2010, and February 2011. Thirty-seven analytes in lagoon slurry and sludge were assessed using LC/MS-MS, and yeast estrogen screen was used to determine estrogenic activity. Of the hormone analytes, steroidal estrogens were more abundant than androgens or progesterone, with estrone being the predominant estrogen species. Conjugated hormones were detected only at low levels. The isoflavone metabolite equol was by far the predominant phytoestrogen species, with daidzein, genistein, formononetin, and coumestrol present at lower levels. Phytoestrogens were often more abundant than steroidal estrogens, but contributed minimally toward total estrogenic activity. Analytes were significantly elevated in the solid phases of the lagoon; although low observed log KOC values suggest enhanced solubility in the aqueous phase, perhaps due to dissolved or colloidal organic carbon. The association with the solid phase, as well as recalcitrance of analytes to anaerobic degradation, results in a markedly elevated load of analytes and estrogenic activity within lagoon sludge. Overall, findings emphasize the importance of adsorption and transformation processes in governing the fate of these compounds in lagoon waste, which is ultimately used for broadcast application as a fertilizer.

Phenotypic characterization of transgenic Japanese medaka (Oryzias latipes) that express a red fluorescent protein in hepatocytes

Transgenic organisms that express fluorescent proteins are used frequently for in vivo visualization of proteins and cells. The phenotype of a transgenic medaka (Oryzias latipes) strain that expresses a red fluorescent protein (RFP) in hepatocytes was characterized using light and fluorescence microscopy, immunohistochemistry, and transmission electron microscopy (TEM). Expression of RFP was first detected by confocal fluorescence microscopy in the location of the liver bud of live medaka embryos at 60 hr postfertilization (developmental stage 27). Subsequently, RFP signal was observed exclusively in hepatocytes throughout life using fluorescence microscopy in live fish and immunohistochemistry in formalin-fixed, paraffin-embedded liver sections. As the fish aged, prominent intracytoplasmic eosinophilic inclusions immunoreactive for RFP were observed by light microscopy and were correlated with membrane-bound electron dense inclusions on TEM. These results define the onset and location of RFP expression in the Tg(zf.L-fabp:DsRed) medaka and characterize a histologic phenotype that results from RFP expression in hepatocytes.

Immunohistochemical characterization of the hepatic progenitor cell compartment in medaka (Oryzias latipes) following hepatic injury

Laboratory fish species are used increasingly in biomedical research and are considered robust models for the study of regenerative processes. Studies investigating the response of the fish liver to injury have demonstrated the presence of a ductular reaction and oval-like cells in injured and regenerating liver. To date, however, it is unclear if this cell population is the piscine equivalent of oval cells (OCs) or intermediate hepatobiliary cells (IHBCs) identified in rodents and man, respectively. The present study defines the process of OC differentiation in fish liver using histopathology, immunohistochemistry and transmission electron microscopy. To generate OC proliferation in Japanese medaka (Oryzias latipes), hepatic injury was induced by exposure of adult fish to either microcystin LR or dimethylnitrosamine. A transgenic strain of medaka expressing a red fluorescent protein (RFP) exclusively in hepatocytes was used. The morphological response to injury was characterized by a ductular reaction comprised of cytokeratin (CK) AE1/AE3(+) OCs progressing to IHBCs variably positive for CK and RFP and finally mature RFP(+) hepatocytes and CK(+) cholangiocytes. These observations support a bipotential differentiation pathway of fish OCs towards hepatocytes and cholangiocytes. Ultrastructural morphology confirmed the presence of OCs and differentiation towards hepatocytes. These results demonstrated clear similarities between patterns of reaction to injury in fish and mammalian livers. They also confirm the presence of, and support the putative bipotential lineage capabilities of, the fish OC.

Anchoring hepatic gene expression with development of fibrosis and neoplasia in a toxicant-induced fish model of liver injury

Fish have been used as laboratory models to study hepatic development and carcinogenesis but not for pathogenesis of hepatic fibrosis. In this study, a dimethylnitrosamine-induced fish model of hepatic injury was developed in Japanese medaka (Oryzias latipes) and gene expression was anchored with the development of hepatic fibrosis and neoplasia. Exposed livers exhibited mild hepatocellular degenerative changes 2 weeks' postexposure. Within 6 weeks, hepatic fibrosis/cirrhosis was evident with development of neoplasia by 10 weeks. Stellate cell activation and development of fibrosis was associated with upregulation of transforming growth factor beta 1 (tgfb1), tgfb receptor 2, mothers against decapentaplegic homolog 3 (smad3a), smad3b, beta-catenin (ctnnb1), myc, matrix metalloproteinase (mmp2), mmp14a, mmp14b, tissue inhibitors of metalloproteinase (timp) 2a, timp2b, timp3, collagen type I alpha 1a (col1a1a), and col1a1b and a less pronounced increase in mmp13 and col4a1 expression. Tgfb receptor I expression was unchanged. Immunohistochemistry suggested that biliary epithelial cells and stellate cells were the main producers of TGF-β1. This study identified a group of candidate genes likely to be involved in the development of hepatic fibrosis and demonstrated that the TGF-β pathway likely plays a major role in the pathogenesis. These results support the medaka as a viable fish model of hepatic fibrosis.

TCDD disrupts hypural skeletogenesis during medaka embryonic development

Defective bone and cartilage development account for a large number of human birth defects annually. Normal skeletogenesis involves cartilage development in early morphogenesis through a highly coordinated and orchestrated series of events involving commitment and differentiation of mesenchymal cells to chondrocytes followed by a highly programmed process of structural maturation. Recent developmental studies with laboratory model fish demonstrate that exposure to 2,3,7,8-tetrachlorodibenzo-p-dioxin (TCDD) results in cartilage and skeletal abnormalities. In this study, we exposed embryonic medaka to TCDD to induce developmental modification(s) of both cartilage and bone formation. Emphasis is placed on cell-rich hyaline cartilage of the hypural plate where both chondrogenesis and osteogenesis are impaired by TCDD exposure. In this model, TCDD exposure results in a concentration-dependent impairment of mesenchymal cell recruitment, chondrocyte cell proliferation, differentiation, and progression to hypertrophy. Gene expression of ColA2, a marker of chondrocyte terminal differentiation in hypural structures, is markedly attenuated consistent with hypural dysmorphogenesis. Assessment of hypural structure using a transgenic medaka expressing mCherry under control of the osterix promoter illustrated significant attenuation in expression of the osteoblast gene marker and lack of formation of a calcified perichondral sheath surrounding hypural anlage. Overall, these studies illustrate that TCDD impacts terminal differentiation and growth of cartilage and bone in axial structures not likely derived from neural crest progenitors in medaka hypurals.

The evolution of farnesoid X, vitamin D, and pregnane X receptors: insights from the green-spotted pufferfish (Tetraodon nigriviridis) and other non-mammalian species

Background

The farnesoid X receptor (FXR), pregnane X receptor (PXR), and vitamin D receptor (VDR) are three closely related nuclear hormone receptors in the NR1H and 1I subfamilies that share the property of being activated by bile salts. Bile salts vary significantly in structure across vertebrate species, suggesting that receptors binding these molecules may show adaptive evolutionary changes in response. We have previously shown that FXRs from the sea lamprey (Petromyzon marinus) and zebrafish (Danio rerio) are activated by planar bile alcohols found in these two species. In this report, we characterize FXR, PXR, and VDR from the green-spotted pufferfish (Tetraodon nigriviridis), an actinopterygian fish that unlike the zebrafish has a bile salt profile similar to humans. We utilize homology modelling, docking, and pharmacophore studies to understand the structural features of the Tetraodon receptors.

Results

Tetraodon FXR has a ligand selectivity profile very similar to human FXR, with strong activation by the synthetic ligand GW4064 and by the primary bile acid chenodeoxycholic acid. Homology modelling and docking studies suggest a ligand-binding pocket architecture more similar to human and rat FXRs than to lamprey or zebrafish FXRs. Tetraodon PXR was activated by a variety of bile acids and steroids, although not by the larger synthetic ligands that activate human PXR such as rifampicin. Homology modelling predicts a larger ligand-binding cavity than zebrafish PXR. We also demonstrate that VDRs from the pufferfish and Japanese medaka were activated by small secondary bile acids such as lithocholic acid, whereas the African clawed frog VDR was not.

Conclusions

Our studies provide further evidence of the relationship between both FXR, PXR, and VDR ligand selectivity and cross-species variation in bile salt profiles. Zebrafish and green-spotted pufferfish provide a clear contrast in having markedly different primary bile salt profiles (planar bile alcohols for zebrafish and sterically bent bile acids for the pufferfish) and receptor selectivity that matches these differences in endogenous ligands. Our observations to date present an integrated picture of the co-evolution of bile salt structure and changes in the binding pockets of three nuclear hormone receptors across the species studied.

Use of medaka in toxicity testing

Small aquarium fishes are increasingly used as animal models, and one of these, the Japanese Medaka (Oryzias latipes), is frequently utilized for toxicity testing. While these vertebrates have many similarities with their terrestrial counterparts, there are differences that must be considered if these organisms are to be used to their highest potential. Commonly, testing may employ either the developing embryo or adults; both are easy to use and work with. To illustrate the utility and breadth of toxicity testing possible using medaka fish, we present protocols for assessing neurotoxicity in developing embryos, evaluating toxicant effects on sexual phenotype after treatment with endocrine-disrupting chemicals by sexual genotyping, and measuring hepatotoxicity in adult fish after treatment with a model hepatotoxicant. The methods run the gamut from immunohistology through PCR to basic histological techniques.

TCDD induced pericardial edema and relative COX-2 expression in medaka (Oryzias Latipes) embryos

Exposure to dioxin and other aryl hydrocarbon receptor (AhR) ligands results in multiple, specific developmental cardiovascular phenotypes including pericardial edema and circulatory failure in small aquarium fish models. Although phenotypes are well described, mechanistic underpinnings for such toxicities remain elusive. Here we suggest that AhR activation results in stimulation of inflammation and "eicosanoid" pathways, which contribute to the observed developmental, cardiovascular phenotypes. We demonstrate that medaka embryos exposed to 2,3,7,8-tetrachlorodibenzo-p-dioxin (TCDD) (0.05-1 ppb) during early development result in a dose-related increase in the prevalence of pericardial edema and that this phenotype correlates with an increase in cyclooxygenase-2 (COX-2) gene expression. Those individuals exhibiting the edema phenotype had significantly greater COX-2 mRNA than their nonedematous cohort. Selective pharmacological inhibition of COX-2, with NS-398, and genetic knock down of COX-2 with a translation initiation morpholino significantly attenuated prevalence and severity of edema phenotype. Subsequently, exposures of medaka embryos to arachidonic acid (AA) resulted in recapitulation of the pericardial edema phenotype and significantly increased COX-2 expression only in those individuals exhibiting the edema phenotype compared with their nonedematous cohort. AA exposure does not result in significant induction of cytochrome P450 1A expression, suggesting that pericardial edema can be induced independent of AhR/aryl hydrocarbon receptor nuclear translocator/dioxin response element interactions. Results from this study demonstrate that developmental exposure to TCDD results in an induction of inflammatory mediators including COX-2, which contribute to the onset, and progression of heart dysmorphogenesis in the medaka model.

Two farnesoid X receptor alpha isoforms in Japanese medaka (Oryzias latipes) are differentially activated in vitro

The nuclear receptor farnesoid X receptor alpha (FXRalpha, NR1H4) is activated by bile acids in multiple species including mouse, rat, and human and in this study we have identified two isoforms of Fxralpha in Japanese medaka (Oryzias latipes), a small freshwater teleost. Both isoforms share a high amino acid sequence identity to mammalian FXRalpha (approximately 70% in the ligand-binding domain). Fxralpha1 and Fxralpha2 differ within the AF1 domain due to alternative splicing at the fourth intron-exon boundary. This process results in Fxralpha1 having an extended N-terminus compared to Fxralpha2. A Gal4DBD-FxralphaLBD fusion construct was activated by chenodeoxycholic, cholic, deoxycholic and lithocholic acids, and the synthetic agonist GW4064 in transient transactivation assays. Activation of the Gal4DBD-FxralphaLBD fusion construct was enhanced by addition of PGC-1alpha, as demonstrated through titration assays. Surprisingly, when the full-length versions of the two Fxralpha isoforms were compared in transient transfection assays, Fxralpha2 was activated by C(24) bile acids and GW4064, while Fxralpha1 was not significantly activated by any of the compounds tested. Since the only significant difference between the full-length constructs was sequence in the AF1 domain, these experiments highlight a key functional region in the Fxralpha AF1 domain. Furthermore, mammalian two-hybrid studies demonstrated the ability of Fxralpha2, but not Fxralpha1, to interact with PGC-1alpha and SRC-1, and supported our results from the transient transfection reporter gene activation assays. These data demonstrate that both mammalian and teleost FXR (Fxralpha2 isoform) are activated by primary and secondary bile acids.

Exposure to the synthetic FXR agonist GW4064 causes alterations in gene expression and sublethal hepatotoxicity in eleutheroembryo medaka (Oryzias latipes)

The small freshwater teleost, medaka (Oryzias latipes), has a history of usage in studies of chronic toxicity of liver and biliary system. Recent progress with this model has focused on defining the medaka hepatobiliary system. Here we investigate critical liver function and toxicity by examining the in vivo role and function of the farnesoid X receptor alpha (FXRalpha, NR1H4), a member of the nuclear receptor superfamily that plays an essential role in the regulation of bile acid homeostasis. Quantitative mRNA analysis of medaka FXRalpha demonstrates differential expression of two FXRalpha isoforms designated Fxralpha1 and Fxralpha2, in both free swimming medaka embryos with remaining yolk (eleutheroembryos, EEs) and adults. Activation of medaka Fxralpha in vivo with GW4064 (a strong FXRalpha agonist) resulted in modification of gene expression for defined FXRalpha gene targets including the bile salt export protein, small heterodimer partner, and cytochrome P450 7A1. Histological examination of medaka liver subsequent to GW4064 exposure demonstrated significant lipid accumulation, cellular and organelle alterations in both hepatocytes and biliary epithelial cells of the liver. This report of hepatobiliary injury following GW4064 exposure extends previous investigations of the intrahepatic biliary system in medaka, reveals sensitivity to toxicant exposure, and illustrates the need for added resolution in detection and interpretation of toxic responses in this vertebrate.

Non invasive in vivo investigation of hepatobiliary structure and function in STII medaka (Oryzias latipes): methodology and applications

Background

A novel transparent stock of medaka (Oryzias latipes; STII), recessive for all pigments found in chromatophores, permits transcutaneous imaging of internal organs and tissues in living individuals. Findings presented describe the development of methodologies for non invasive in vivo investigation in STII medaka, and the successful application of these methodologies to in vivo study of hepatobiliary structure, function, and xenobiotic response, in both 2 and 3 dimensions.

Results

Using brightfield, and widefield and confocal fluorescence microscopy, coupled with the in vivo application of fluorescent probes, structural and functional features of the hepatobiliary system, and xenobiotic induced toxicity, were imaged at the cellular level, with high resolution (< 1 μm), in living individuals. The findings presented demonstrate; (1) phenotypic response to xenobiotic exposure can be investigated/imaged in vivo with high resolution (< 1 μm), (2) hepatobiliary transport of solutes from blood to bile can be qualitatively and quantitatively studied/imaged in vivo, (3) hepatobiliary architecture in this lower vertebrate liver can be studied in 3 dimensions, and (4) non invasive in vivo imaging/description of hepatobiliary development in this model can be investigated.

Conclusion

The non-invasive in vivo methodologies described are a unique means by which to investigate biological structure, function and xenobiotic response with high resolution in STII medaka. In vivo methodologies also provide the future opportunity to integrate molecular mechanisms (e.g., genomic, proteomic) of disease and toxicity with phenotypic changes at the cellular and system levels of biological organization. While our focus has been the hepatobiliary system, other organ systems are equally amenable to in vivo study, and we consider the potential for discovery, within the context of in vivo investigation in STII medaka, as significant.

Paralogous vitamin D receptors in teleosts: transition of nuclear receptor function

The availability of multiple teleost (bony fish) genomes is providing unprecedented opportunities to understand the diversity and function of gene duplication events using comparative genomics. Here we describe the cloning and functional characterization of two novel vitamin D receptor (VDR) paralogs from the freshwater teleost medaka (Oryzias latipes). VDR sequences were identified through mining of the medaka genome database in which gene organization and structure was determined. Two distinct VDR genes were identified in the medaka genome and mapped to defined loci. Each VDR sequence exhibits unique intronic organization and dissimilar 5' untranslated regions, suggesting they are not isoforms of the same gene locus. Phylogenetic comparison with additional teleosts and mammalian VDR sequences illustrate that two distinct clusters are formed separating aquatic and terrestrial species. Nested within the teleost cluster are two separate clades for VDRalpha and VDRbeta. The topology of teleost VDR sequences is consistent with the notion of paralogous genes arising from a whole genome duplication event prior to teleost radiation. Functional characterization was conducted through the development of VDR expression vectors including Gal4 chimeras containing the yeast Gal4 DNA binding domain fused to the medaka VDR ligand binding domain and full-length protein. The common VDR ligand 1alpha,25-dihydroxyvitamin D3 [1alpha,25(OH)(2)D(3)] resulted in significant transactivation activity with both the Gal4 and full-length constructs of medaka (m) VDRbeta. Comparatively, transactivation of mVDRalpha with 1alpha,25(OH)(2)D(3) was highly attenuated, suggesting a functional divergence between these two nuclear receptor paralogs. We additionally demonstrate through coactivator studies that mVDRalpha is still functional; however, it exhibits a different sensitivity to 1alpha,25(OH)(2)D(3), compared with VDRbeta. These results suggest that in mVDRalpha and VDRbeta have undergone a functional divergence through a process of sub- and/or neofunctionalization of VDR nuclear receptor gene pairs.

Differential developmental toxicity of naphthoic acid isomers in medaka (Oryzias latipes) embryos

Polycyclic aromatic hydrocarbons (PAHs) are widespread persistent pollutants that readily undergo biotic and abiotic conversion to numerous transformation products in rivers, lakes and estuarine sediments. Here we characterize the developmental toxicity of four PAH transformation products each structural isomers of hydroxynaphthoic acid: 1H2NA, 2H1NA, 2H3NA, and 6H2NA. Medaka fish (Oryzias latipes) embryos and eleutheroembryos were used to determine toxicity. A 96-well micro-plate format was used to establish a robust, statistically significant platform for assessment of early life stages. Individual naphthoic acid isomers demonstrated a rank order of toxicity with 1H2NA>2H1NA>2H3NA>6H2NA being more toxic. Abnormalities of circulatory system were most pronounced including pericardial edema and tube heart. To determine if HNA isomers were AhR ligands, spatial-temporal expression and activity of CYP1A was measured via in vivo EROD assessments. qPCR measurement of CYP1A induction proved different between isomers dosed at respective concentrations affecting 50% of exposed individuals (EC50s). In vitro, all ANH isomers transactivated mouse AhR using a medaka CYP1A promoter specific reporter assay. Circulatory abnormalities followed P450 induction and response was consistent with PAH toxicity. A 96-well micro-plates proved suitable as exposure chambers and provided statistically sound evaluations as well as efficient toxicity screens. Our results demonstrate the use of medaka embryos for toxicity analysis thereby achieving REACH objectives for the reduction of adult animal testing in toxicity evaluations.

An in vivo look at vertebrate liver architecture: three-dimensional reconstructions from medaka (Oryzias latipes)

Understanding three-dimensional (3D) hepatobiliary architecture is fundamental to elucidating structure/function relationships relevant to hepatobiliary metabolism, transport, and toxicity. To date, factual information on vertebrate liver architecture in 3 dimensions has remained limited. Applying noninvasive in vivo imaging to a living small fish animal model we elucidated, and present here, the 3D architecture of this lower vertebrate liver. Our investigations show that hepatobiliary architecture in medaka is based on a polyhedral (hexagonal) structural motif, that the intrahepatic biliary system is an interconnected network of canaliculi and bile-preductules, and that parenchymal architecture in this lower vertebrate is more related to that of the mammalian liver than previously believed. The in vivo findings presented advance our comparative 3D understanding of vertebrate liver structure/function, help clarify previous discrepancies among vertebrate liver conceptual models, and pose interesting questions regarding the "functional unit" of the vertebrate liver.

Comparisons of polychromatic and monochromatic UV-based treatments of bisphenol-A in water via toxicity assessments

Polychromatic ultraviolet irradiation, such as from medium pressure (MP) Hg lamps may enhance the UV degradation of environmental pollutants as compared to low pressure (LP) Hg UV sources emitting monochromatic irradiation. Typically, studies involving destruction of environmental pollutants such as endocrine disrupting compounds (EDCs) are based on measurement of the parent compound decay using analytical chemistry, but such information is insufficient to determine an effective treatment endpoint because the identity and biological activity of many transformation products remain unknown. Bioanalytical methods to assess residual biological activity of a treated water offers one means to compare removal efficiency of EDC activity between MP- and LP-UV lamps under photolysis and UV/H2O2 oxidation. In this study, changes in estrogenic activity of bisphenol-A (BPA) as a function of UV treatment were evaluated using both an in vitro yeast estrogen screen and in vivo vitellogenin assay with Japanese medaka (Oryzias latipes) fish. Decay of BPA parent compound and formation of degradation products were followed using HPLC analysis. Results demonstrated that MP-UV direct photolysis more effectively removed BPA and associated estrogenic activity compared to LP-UV lamps. UV in combination with H2O2 significantly removed estrogenic activity in vitro and in vivo compared to direct photolysis; however, no significant difference in removal rates was found between the two lamps under UV/H2O2 oxidation. Furthermore, the UV/H2O2 process was effective for reducing embryo toxicity of BPA, but resulted in the production of acidic intermediates, causing acute toxicity and delayed hatching in some medaka embryos.

Biological assessments of a mixture of endocrine disruptors at environmentally relevant concentrations in water following UV/H2O2 oxidation

Numerous studies have investigated degradation of individual endocrine disrupting compounds (EDCs) in lab or natural waters. However, natural variations in water matrices and mixtures of EDCs in the environment may confound analysis of the treatment efficiency. Because chemical based analytical methods cannot represent the combined or synergistic activities between water quality parameters and/or the EDC mixtures at environmentally relevant concentrations (microg L(-1)-ng L(-1)), bioanalytical assessments of residual estrogenic activity in treated water were used to evaluate the performance of the UV based advanced oxidation process for estrogenic contaminants in water. Four EDCs including estradiol (E(2)), ethinyl estradiol (EE(2)), bisphenol-A (BPA) and nonylphenol (NP) were spiked individually or as a mixture at mug L(-1)-ng L(-1) in laboratory or natural river water. The removal rates of estrogenic activity were quantitatively evaluated by in vitro yeast estrogen screen (YES) and in vivo Vitellogenin (VTG) assays with Japanese medaka fish (Oryzias latipes). UV in combination with 10 ppm H(2)O(2) as an oxidation process was capable of decreasing in vitro and in vivo estrogenic activity, however, in vivo estrogenic activity of the EDC mixture in natural water was not completely removed at UV fluence up to 2000 mJ cm(-2). The removal rates of in vitro estrogenic activity of the EDC mixtures were lower than those observed for single compounds, and slower in natural waters, likely due to lower steady-state concentrations of hydroxyl radicals (*OH) in the presence of *OH scavengers from the water matrix and EDC mixture.

Estrogenic modulation of CYP3A38, CYP3A40, and CYP19 in mature male medaka (Oryzias latipes)

We examined cytochrome P450 production and activity and circulating hormone concentrations in male medaka exposed to 17beta-estradiol (E2) or 17alpha-ethinylestradiol (EE2). Intraperitoneal injection of E2 at 1, 10, or 100 microg/g-fish completely suppressed CYP3A38 protein production and suppressed CYP3A40 protein levels by 89%, 52%, or 47%, respectively. CYP3A38 and CYP3A40 mRNA expression was unaltered, and CYP3A enzymatic activity initially increased and then decreased with increasing E2 dose. Males co-cultured with females were exposed to a markedly high concentration (43 ng/L) of E2 secreted by females. CYP3A protein levels in co-cultured males were suppressed. Serum testosterone (TE) and 11keto-testosterone levels in co-cultured males were downregulated to 40% of pre-exposure levels. Serum E2 levels increased in co-cultured males or males exposed to EE2. Testicular CYP19, which converts TE to E2, increased by 9.5 times in males exposed to 50 ng/L EE2 and by 21.5 times in those exposed to 100 ng/L EE2. Male medaka exposed to EE2 showed increased serum Vtg levels. Estrogenic exposure induced Vtg production, suppressed CYP3A protein production, downregulated TE metabolism, and enhanced CYP19 activity. Serum E2 endogenously induced by CYP19 could contribute to Vtg induction in male medaka.

UV/H2O2 degradation of endocrine-disrupting chemicals in water evaluated via toxicity assays

Due to rising concern regarding the presence of endocrine-disrupting chemicals (EDCs) in surface water and groundwater throughout the United States, Asia and Europe, treatment of these chemicals in drinking water and wastewater to protect human health and the environment is an area of great interest. Many conventional treatment schemes are relatively ineffective in removing EDCs from water and wastewater. This is concerning because these chemicals are biologically active at very low concentrations and effects of mixtures are relatively unknown. 17-alpha-oestradiol (E2) and 17-beta-ethinyl-oestradiol (EE2), suspected EDCs, were degraded significantly by the UV/H2O2 AOP. The UV/H2O2 processes using either low or medium pressure lamps were degraded EDCs by between 80 and 99.3% at a 15 ppm H2O2 concentration and a UV dose of 1,000 mJ/cm2. Significantly greater removal was noted when the removal was based on total oestrogenic activity using a yeast oestrogen screen (YES) assay. These data indicated that a dose of less than 200 mJ/cm2 completely removed oestrogenic activity in lab water. Values for natural waters were slightly higher. A steady state model was developed to determine EDC destruction efficiency in waters of differing quality. The model effectively predicted destruction in water, where concentrations of all scavenging species were known. Based on these results it was concluded than complete destruction of oestrogenic activity was possible under practical advanced oxidation conditions for a variety of water qualities.