RNA-sequencing analysis of TCDD-induced responses in zebrafish liver reveals high relatedness to in vivo mammalian models and conserved biological pathways

Whole-Transcriptome Analysis Reveals the RNA Profiles in Mouse Bone Marrow Mesenchymal Stem Cells or Zebrafish Embryos After Exposure to Environmental Level of Tributyltin

To understand the underlying molecular mechanisms, mouse bone marrow mesenchymal stem cells (BMSCs) and zebrafish embryos were exposed to the control group and Tributyltin (TBT) group (10 ng/L, environmental concentration) for 48 h, respectively. The expression profiles of RNAs were investigated using whole-transcriptome analysis in mouse BMSCs or zebrafish embryos after TBT exposure. For mouse BMSCs, the results showed 2,449 differentially expressed (DE) mRNAs, 59 DE miRNAs, 317 DE lncRNAs, and 15 circRNAs. Similarly, for zebrafish embryos, the results showed 1,511 DE mRNAs, 4 DE miRNAs, 272 DE lncRNAs, and 28 circRNAs. According to KEGG pathway analysis showed that DE RNAs were mainly associated with immune responses, signaling, and cellular interactions. Competing endogenous RNA (ceRNA) network analysis revealed that the regulatory network of miRNA-circRNA constructed in zebrafish embryos was more complex compared to that of mouse BMSCs.

Toxicity assessment of hexafluoropropylene oxide-dimer acid on morphology, heart physiology, and gene expression during zebrafish (Danio rerio) development

Hexafluoropropylene oxide-dimer acid (HFPO-DA) is one of the emerging replacements for the "forever" carcinogenic and toxic long-chain PFAS. HFPO-DA is a polymerization aid used for manufacturing fluoropolymers, whose global distribution and undetermined toxic properties are a concern regarding human and ecological health. To assess embryotoxic potential, zebrafish embryos were exposed to HFPO-DA at concentrations of 0.5-20,000 mg/L at 24-, 48-, and 72-h post-fertilization (hpf). Heart rate increased significantly in embryos exposed to 2 mg/L and 10 mg/L HFPO-DA across all time points. Spinal deformities and edema phenotypes were evident among embryos exposed to 1000-16,000 mg/L HFPO-DA at 72 hpf. A median lethal concentration (LC₅₀) was derived as 7651 mg/L at 72 hpf. Shallow RNA sequencing analysis of 9465 transcripts identified 38 consistently differentially expressed genes at 0.5 mg/L, 1 mg/L, 2 mg/L, and 10 mg/L HFPO-DA exposures. Notably, seven downregulated genes were associated with visual response, and seven upregulated genes were expressed in or regulated the cardiovascular system. This study identifies biological targets and molecular pathways affected during animal development by an emerging, potentially problematic, and ubiquitous industrial chemical.

A New In Vivo Zebrafish Bioassay Evaluating Liver Steatosis Identifies DDE as a Steatogenic Endocrine Disruptor, Partly through SCD1 Regulation

Non-alcoholic fatty liver disease (NAFLD), which starts with liver steatosis, is a growing worldwide epidemic responsible for chronic liver diseases. Among its risk factors, exposure to environmental contaminants, such as endocrine disrupting compounds (EDC), has been recently emphasized. Given this important public health concern, regulation agencies need novel simple and fast biological tests to evaluate chemical risks. In this context, we developed a new in vivo bioassay called StAZ (Steatogenic Assay on Zebrafish) using an alternative model to animal experimentation, the zebrafish larva, to screen EDCs for their steatogenic properties. Taking advantage of the transparency of zebrafish larvae, we established a method based on fluorescent staining with Nile red to estimate liver lipid content. Following testing of known steatogenic molecules, 10 EDCs suspected to induce metabolic disorders were screened and DDE, the main metabolite of the insecticide DDT, was identified as a potent inducer of steatosis. To confirm this and optimize the assay, we used it in a transgenic zebrafish line expressing a blue fluorescent liver protein reporter. To obtain insight into DDE's effect, the expression of several genes related to steatosis was analyzed; an up-regulation of scd1 expression, probably relying on PXR activation, was found, partly responsible for both membrane remodeling and steatosis.

NMR technique revealed the metabolic interference mechanism of the combined exposure to cadmium and tributyltin in grass carp larvae

Widespread human activity has resulted in the presence of different pollutants in the aquatic environment that does not exist in isolation. The study of the effects of contamination of aquatic organisms is of great significance. To assess the individual and combined toxicity of cadmium (Cd) and tributyltin (TBT) to aquatic organisms, juvenile grass carp (Ctenopharyngodon idella) were exposed to Cd (2.97 mg/L), TBT (7.5 μg/L), and their mixture MIX. The biological response was evaluated by nuclear magnetic resonance (NMR) analysis of plasma metabolites. Plasma samples at 1, 2, 4, 8, 16, 32, and 48 days post-exposure were analyzed using detection by NMR technique. The typical correlation analysis (CCA) analysis revealed that TBT had the greatest effect on plasma metabolism, followed by MIX and Cd. The interference pathway to grass carp was similar to that of TBT and MIX. Both Cd and TBT exposure alone or in combination can lead to metabolic abnormalities in TCA cycle-related pathways and interfere with energy metabolism. These results provide more detailed information for the metabolic study of pollutants and data for assessing the health risks of Cd, TBT, and MIX at the metabolic level.

A Mini-review of the Toxicity of Pollutants to Fish Under Different Salinities

In recent years, with the development of the global economy, water pollution has increased. Pollutants migrate, accumulate, and diffuse in aquatic environments. Most of the pollutants eventually enter aquatic organisms. The accumulation of pollutants affects the development and reproduction of organisms, and many pollutants have teratogenic, carcinogenic, and/or mutagenic effects. Aquatic organisms in estuaries and coastal areas are under pressure due to both salinity and pollutants. Among them, salinity, as an environmental factor, may affect the behavior of pollutants in the aquatic environment, causing changes in their toxic effects on fishes. Salinity also directly affects the growth and development of fishes. Therefore, this paper focuses on metals and organic pollutants and discusses the toxic effects of pollutants on fish under different salinities. This research is of great significance to environmental protection and ecological risk assessment of aquatic environments.

Transcriptomic and proteomic analysis of Chinese rare minnow (Gobiocypris rarus) larvae in response to acute waterborne cadmium or mercury stress

In this study, Chinese rare minnow (Gobiocypris rarus) larvae were exposed to the control group, Cd concentrations (0.5 and 2.5 mg/L), and Hg concentrations (0.1 and 0.3 mg/L) for 96 h. Transcriptome analysis showed that 816 and 1599 significantly differentially expressed genes (DEGs) were identified in response to 2.5 mg/L Cd²⁺ and 0.3 mg/L Hg²⁺, respectively. Functional enrichment analysis revealed that DEGs were mostly associated with immune responses after Cd exposure, such as antigen processing and presentation, phagosome, apoptosis, and lysosome. Similarly, functional enrichment analysis showed that many pathways were mostly involved in metabolism after Hg exposure, such as glutathione metabolism and starch and sucrose metabolism. Results of two-dimensional electrophoresis (2-DE) showed that the abundance of 10 protein spots was significantly altered in the Cd²⁺ treatments. The proteomic analysis demonstrated that Cd toxicity might impair cytoskeletal and cell motility-related protein activity in the liver of G. rarus. Similarly, the abundance of 24 protein spots was significantly altered in the Hg²⁺ treatments. Hg toxicity regulates the expression of proteins belonging to several functional categories, including cytoskeleton, oxidative stress, digestive system, and energy metabolism. This study provides valuable relevant insight into the molecular mechanisms in response to Cd or Hg toxicity in aquatic organisms and will help screen for potential biomarkers to respond to Cd and Hg pollutants.

Immune defence mechanisms presented in liver homogenates and bile of gilthead seabream (Sparus aurata)

Because the role of the liver of fishes in providing possible immunity remains largely unknown, the aim of this work was to identify and characterize different humoral defence mechanisms in the liver homogenates and bile of gilthead seabream (Sparus aurata) for the first time. Total protein levels and several immune parameters (complement activity, lysozyme and immunoglobulin M level) were studied. Furthermore, the activity of some lytic (proteases, antiproteases, esterase, alkaline phosphatase) and antioxidant (superoxide dismutase, catalase and peroxidase) enzymes was determined. Finally, bacteriostatic activity on three opportunist fish pathogens (Vibrio harveyi, Vibrio angillarum and Photobacterium damselae) was measured. Lysozyme and antiprotease activity were undetected in liver and bile, while natural haemolytic complement activity was only detected in bile, and immunoglobulin M was detected in both samples. The levels of proteases, esterase and antioxidant enzymes were greater in bile than in liver homogenates, while the level of alkaline phosphatase was very low in both samples. In addition, while no bacteriostatic activity was detected on liver homogenates, the bile revealed a very potent bacteriostatic activity against all the tested pathogenic bacteria. These results corroborate that fish liver - especially fish bile - contains many factors involved in innate immunity that could be useful for better understanding the role of the liver as an organ involved in fish immune functions as well as the possible contribution of bile to gut mucosal immunity.

Immunotoxicity of Xenobiotics in Fish: A Role for the Aryl Hydrocarbon Receptor (AhR)?

The impact of anthropogenic contaminants on the immune system of fishes is an issue of growing concern. An important xenobiotic receptor that mediates effects of chemicals, such as halogenated aromatic hydrocarbons (HAHs) and polyaromatic hydrocarbons (PAHs), is the aryl hydrocarbon receptor (AhR). Fish toxicological research has focused on the role of this receptor in xenobiotic biotransformation as well as in causing developmental, cardiac, and reproductive toxicity. However, biomedical research has unraveled an important physiological role of the AhR in the immune system, what suggests that this receptor could be involved in immunotoxic effects of environmental contaminants. The aims of the present review are to critically discuss the available knowledge on (i) the expression and possible function of the AhR in the immune systems of teleost fishes; and (ii) the impact of AhR-activating xenobiotics on the immune systems of fish at the levels of immune gene expression, immune cell proliferation and immune cell function, immune pathology, and resistance to infectious disease. The existing information indicates that the AhR is expressed in the fish immune system, but currently, we have little understanding of its physiological role. Exposure to AhR-activating contaminants results in the modulation of numerous immune structural and functional parameters of fish. Despite the diversity of fish species studied and the experimental conditions investigated, the published findings rather uniformly point to immunosuppressive actions of xenobiotic AhR ligands in fish. These effects are often associated with increased disease susceptibility. The fact that fish populations from HAH- and PAH-contaminated environments suffer immune disturbances and elevated disease susceptibility highlights that the immunotoxic effects of AhR-activating xenobiotics bear environmental relevance.

The aryl hydrocarbon receptor activates ceramide biosynthesis in mice contributing to hepatic lipogenesis

Aryl hydrocarbon receptor (AHR) activation via 2,3,7,8-tetrachlorodibenzofuran (TCDF) induces the accumulation of hepatic lipids. Here we report that AHR activation by TCDF (24  μg/kg body weight given orally for five days) induced significant elevation of hepatic lipids including ceramides in mice, was associated with increased expression of key ceramide biosynthetic genes, and increased activity of their respective enzymes. Results from chromatin immunoprecipitation (ChIP), electrophoretic mobility shift assay (EMSA) and cell-based reporter luciferase assays indicated that AHR directly activated the serine palmitoyltransferase long chain base subunit 2 (Sptlc2, encodes serine palmitoyltransferase 2 (SPT2)) gene whose product catalyzes the initial rate-limiting step in de novo sphingolipid biosynthesis. Hepatic ceramide accumulation was further confirmed by mass spectrometry-based lipidomics. Taken together, our results revealed that AHR activation results in the up-regulation of Sptlc2, leading to ceramide accumulation, thus promoting lipogenesis, which can induce hepatic lipid accumulation.

Identification of molecular toxicity pathways across early life-stages of zebrafish exposed to PCB126 using a whole transcriptomics approach

Although withdrawn from the market in the 1980s, polychlorinated biphenyls (PCBs) are still found ubiquitously in the aquatic environment and pose a serious risk to biota due to their teratogenic potential. In fish, early life-stages are often considered most sensitive with regard to their exposure to PCBs and other dioxin-like compounds. However, little is known about the molecular drivers of the frequently observed teratogenic effects. Therefore, the aims of our study were to: (1) characterize the baseline transcriptome profiles at different embryonic life-stages in zebrafish (Danio rerio); and (2) to identify the molecular response to PCB exposure and life-stage specific-effects of the chemical on associated processes. For both objectives, embryos were sampled at 12, 48, and 96 h post-fertilization (hpf) and subjected to Illumina sequence-by-synthesis and RNAseq analysis. Results revealed that with increasing age more genes and related pathways were upregulated both in terms of number and magnitude. Yet, other transcripts followed an opposite pattern with greater transcript abundance at the earlier time points. Additionally, embryos were exposed to PCB126, a potent agonist of the aryl hydrocarbon receptor (AHR). ClueGO network analysis revealed significant enrichment of genes associated with basic cell metabolism, communication, and homeostasis as well as eye development, muscle formation, and skeletal formation. We selected eight genes involved in the affected pathways for an in-depth characterization of their regulation throughout normal embryogenesis and after exposure to PCB126 by quantification of transcript abundances every 12 h until 118 hpf. Among these, fgf7 and c9 stood out because of their strong upregulation by PCB126 exposure at 48 and 96 hpf, respectively. Cyp2aa12 was upregulated from 84 hpf on. Fabp10ab, myhz1.1, col8a1a, sulf1, and opn1sw1 displayed specific regulation depending on the developmental stage. Overall, we demonstrate that (1) the developmental transcriptome of zebrafish is highly dynamic, and (2) dysregulation of gene expression by exposure to PCB126 was significant and in several cases not directly connected to AHR-signaling. Hence, this study improves the understanding of linkages between molecular events and apical outcomes that are of regulatory relevance.

Genomic regions underlying metabolic and neuronal signaling pathways are temporally consistent in a moving avian hybrid zone

The study of hybrid zones can provide insight into the genetic basis of species differences that are relevant for the maintenance of reproductive isolation. Hybrid zones can also provide insight into climate change, species distributions, and evolution. The hybrid zone between black-capped chickadees (Poecile atricapillus) and Carolina chickadees (Poecile carolinensis) is shifting northward in response to increasing winter temperatures but is not increasing in width. This pattern indicates strong selection against chickadees with admixed genomes. Using high-resolution genomic data, we identified regions of the genomes that are outliers in both time points and do not introgress between the species; these regions may be involved in the maintenance of reproductive isolation. Genes involved in metabolic regulation processes were overrepresented in this dataset. Several gene ontology categories were also temporally consistent-including glutamate signaling, synaptic transmission, and catabolic processes-but the nucleotide variants leading to this pattern were not. Our results support recent findings that hybrids between black-capped and Carolina chickadees have higher basal metabolic rates than either parental species and suffer spatial memory and problem-solving deficits. Metabolic breakdown, as well as spatial memory and problem-solving, in hybrid chickadees may act as strong postzygotic isolation mechanisms in this moving hybrid zone.

Aryl hydrocarbon receptor ligands increase ABC transporter activity and protein expression in killifish (Fundulus heteroclitus) renal proximal tubules

Many widespread and persistent organic pollutants, for example, 2,3,7,8-tetrachlorodibenzo-p-dioxin (TCDD) and some polychlorinated biphenyls, activate the aryl hydrocarbon receptor (AhR) causing it to translocate to the cell nucleus where it transactivates target genes, increasing expression of a number of xenobiotic metabolizing enzymes as well as some transporters. AhR's ability to target transporters within the kidney is essentially unexplored. We show here that exposing isolated killifish (Fundulus heteroclitus) renal proximal tubules to micromolar β-naphthoflavone (BNF) or nanomolar TCDD roughly doubled the transport activity of Multidrug resistance-associated proteins Mrp2 and Mrp4, P-glycoprotein (P-gp) and Breast cancer resistance protein (Bcrp), all ATP-driven xenobiotic efflux pumps and critical determinants of renal xenobiotic excretion. These effects were abolished by actinomycin D and cycloheximide and by the AhR antagonist, α-naphthoflavone, indicating that increased transport activity was dependent on transcription and translation as well as ligand binding to AhR. Quantitative immunostaining of renal tubules exposed to BNF and TCDD showed increased luminal membrane expression of Mrp2, Mrp4, P-gp and Bcrp. Thus, in these renal tubules, the four ABC transporters are targets of AhR action.

Toxic effects of a methanolic coal dust extract on fish early life stage

Coal dust is a contaminant that impacts the terrestrial and aquatic environment with a complex mixture of chemicals, including PAHs and metals. This study aims to evaluate the toxic effect of a methanolic coal dust extract on a fish early life stage by analyzing phenotypic alterations, transcriptome changes, and mortality in zebrafish (ZF) embryos. ZF embryos were exposed to methanolic coal dust extract at 1-5000 mg·L-1 and monitored using bright field microscopy 24 and 48 hpf to determine malformations and mortality. In situ hybridization, RNA sequencing, and qRT-PCR were employed to identify transcriptome changes in malformed embryos. Three malformed phenotypes were generated in a dose-dependent manner. In situ hybridization analysis revealed brain, somite, dorsal cord, and heart tube development biomarker alterations. Gene expression profile analysis identified changes in genes related to structural constituent of muscle, calcium ion binding, actin binding, melanin metabolic process, muscle contraction, sarcomere organization, cardiac myofibril assembly, oxidation-reduction process, pore complex, supramolecular fiber, striated muscle thin filament, Z disc, and intermediate filament. This study shows, for the first time, the malformations generated by a mixture of pollutants from a methanolic coal dust extract on a fish early life stage, constituting a potential risk for normal embryonic development of other aquatic vertebrate organisms. Furthermore, we establish that phenotypes and changes in gene expression induced by the extract constitute a target for future studies about mechanical toxicity and their utility as sensitive tools in environmental risk assessments for biota and humans exposed to coal mining activities.

Comparative analysis of the liver transcriptome in the red-eared slider Trachemys scripta elegans under chronic salinity stress

The red-eared slider (Trachemys scripta elegans), identified as one of the 100 most invasive species in the world, is a freshwater turtle originally from the eastern United States and northeastern Mexico. Field investigations have shown that T. s. elegans can survive and lay eggs in saline habitats. In order to understand the molecular mechanisms of salinity adaptation, high-throughput RNA-Seq was utilized to identify the changes in gene expression profiles in the liver of T. s. elegans in response to elevated salinity. We exposed individuals to 0, 5, or 15 psu (practical salinity units) for 30 days. A total of 157.21 million reads were obtained and assembled into 205138 unigenes with an average length of 620 bp and N50 of 964 bp. Of these, 1019 DEGs (differentially expressed genes) were found in the comparison of 0 vs. 5 psu, 1194 DEGs in 0 vs. 15 psu and 1180 DEGs in 5 vs. 15 psu, which are mainly related to macromolecule metabolic process, ion transport, oxidoreductase activity and generation of precursor metabolites and energy by GO (Gene Ontology) enrichment analyses. T. s. elegans can adapt itself into salinity by balancing the entry of sodium and chloride ions via the up-regulation expression genes of ion transport (potassium voltage-gated channel subfamily H member 5, KCNH5; erine/threonine-protein kinase 32, STK32; salt-inducible kinase 1, SIK1; adiponectin, ACDC), and by accumulating plasma urea and free amino acid via the up-regulation expression genes of amino acid metabolism (ornithine decarboxylase antizyme 3, OAZ3; glutamine synthetase, GLUL; asparaginase-like protein 1b, ASRGL; L-amino-acid oxidase-like, LAAO; sodium-dependent neutral amino acid transporter B, SLC6A15s; amino acid permease, SLC7A9) in response to osmotic regulation. An investment of energy to maintain their homeostatic balance is required to salinity adaptation, therefore, the genes related to energy production and conversion (F-ATPase protein 6, ATP6; cytochrome c oxidase subunit I, COX1; cytochrome c oxidase subunit III, COX3; cytochrome b, CYTb; cytochrome P450 17A1, CYP17A1) were up-regulated with the increase of gene expression associated with lipid metabolism (apolipoprotein E precursor, APoE; coenzyme Q-binding protein, CoQ10; high-density lipoprotein particle, SAA) and carbohydrate metabolism (HK, MIP). These findings improve our understanding of the underlying molecular mechanisms involved in salinity adaptation and provide general guidance to illuminate the invasion potential of T. s. elegans into saline environments.

Proteomic Analysis of Zebrafish (Danio rerio) After Chemical Exposure

Traditional toxicological screens based on the zebrafish model use observable phenotypic endpoints during their development to determine the toxicity of teratogens. Yet toxicity does not always translate to obvious phenotypic changes and the criteria used to score the toxicity of a teratogen are frequently subjected to human perception. The advancement in omics-based technologies has allowed us to quantitatively and objectively determine the toxicity of a teratogen based on biomolecular changes. The field of proteomics has been gaining popularity as a valuable tool in toxicology. Hence, in this chapter, we described a protocol for both label-free and label-based proteomic methods to analyse proteomic changes in both embryos and adult livers of zebrafish exposed to the teratogen TCDD (tetrachlorodibenzo-p-dioxin) as an example.

A critical review of histopathological findings associated with endocrine and non-endocrine hepatic toxicity in fish models

Although frequently examined as a target organ for non-endocrine toxicity, histopathological evaluation of the liver is becoming a routine component of endocrine disruption studies that utilize various fish species as test subjects. However, the interpretation of microscopic liver findings can be challenging, especially when attempting to distinguish adverse changes associated with endocrine disrupting substances from those caused by systemic or direct hepatic toxicity. The purpose of this project was to conduct a critical assessment of the available peer-reviewed and grey literature concerning the histopathologic effects of reproductive endocrine active substances (EAS) and non-endocrine acting substances in the livers of fish models, and to determine if liver histopathology can be used to reliably distinguish endocrine from non-endocrine etiologies. The results of this review suggest that few compound-specific histopathologic liver effects have been identified, among which are estrogen agonist-induced increases in hepatocyte basophilia and proteinaceous intravascular fluid in adult male teleosts, and potentially, decreased hepatocyte basophilia in female fish exposed to substances that possess androgenic, anti-estrogenic, or aromatase inhibitory activity. This review also used published standardized methodology to assess the credibility of the histopathology data in each of the 117 articles that reported liver effects of treatment, and consequently it was determined that in only 37% of those papers were the data considered either highly credible or credible. The outcome of this work highlights the value of histopathologic liver evaluation as an investigative tool for EAS studies, and provides information that may have implications for EAS hazard assessment.

Hepatoprotective and antioxidant effects of dietary Glycyrrhiza polysaccharide against TCDD-induced hepatic injury and RT-PCR quantification of AHR2, ARNT2, CYP1A mRNA in Jian Carp (Cyprinus carpio var. Jian)

To evaluate the protective effects of Glycyrrhiza polysaccharide (GPS) against 2,3,7,8-tetrachlorodibenzo-p-dioxin (TCDD)-induced hepatotoxicity in Jian carp, the fish were fed diets containing GPS at doses of 0.1, 0.5 and 1.0g/kg for 60days before an intraperitoneal injection of 0.6μg/kg TCDD at a volume of 0.05mL/10g body weight. At 72hr post-injection, blood and liver samples were taken for biochemical analysis and the fish liver samples were used for the preparation of pathological slices. The results showed that increases in alanine aminotransferase (GPT), aspartate aminotransferase (GOT), lactate dehydrogenase (LDH), and alkaline phosphatase (AKP) in serum induced by TCDD were significantly inhibited by pre-treatment with 1.0g/kg GPS. Following the 1.0g/kg GPS pre-treatment, total protein (TP), albumin (Alb), catalase (CAT), glutathione peroxidase (GPx), total antioxidant capacity (T-AOC) and superoxide dismutase (SOD) activities in liver tissue increased significantly, malondialdehyde (MDA) formation (P<0.05 or P<0.01) was significantly inhibited, and the expression of cytochrome P4501A (CYP1A), aryl hydrocarbon receptor 2 (AHR2) and aryl hydrocarbon receptor nuclear translocator 2 (ARNT2) mRNA (P<0.05) was significantly enhanced. Histological observations on fish liver were obtained by preparing paraffin tissue sections via HE staining, and the results showed that histological changes were obviously reduced by 0.5 and 1.0g/kg GPS. GPS significantly reduced liver tissue damage caused by TCDD. Overall, these results proved the hepatoprotective effect of GPS in protecting against fish liver injury induced by TCDD, and supported the use of GPS (1.0g/kg) as a hepatoprotective and antioxidant agent in fish.

ZEBRAFISH AS AN IN VIVO MODEL FOR SUSTAINABLE CHEMICAL DESIGN

Heightened public awareness about the many thousands of chemicals in use and present as persistent contaminants in the environment has increased the demand for safer chemicals and more rigorous toxicity testing. There is a growing recognition that the use of traditional test models and empirical approaches is impractical for screening for toxicity the many thousands of chemicals in the environment and the hundreds of new chemistries introduced each year. These realities coupled with the green chemistry movement have prompted efforts to implement more predictive-based approaches to evaluate chemical toxicity early in product development. While used for many years in environmental toxicology and biomedicine, zebrafish use has accelerated more recently in genetic toxicology, high throughput screening (HTS), and behavioral testing. This review describes major advances in these testing methods that have positioned the zebrafish as a highly applicable model in chemical safety evaluations and sustainable chemistry efforts. Many toxic responses have been shown to be shared among fish and mammals owing to their generally well-conserved development, cellular networks, and organ systems. These shared responses have been observed for chemicals that impair endocrine functioning, development, and reproduction, as well as those that elicit cardiotoxicity and carcinogenicity, among other diseases. HTS technologies with zebrafish enable screening large chemical libraries for bioactivity that provide opportunities for testing early in product development. A compelling attribute of the zebrafish centers on being able to characterize toxicity mechanisms across multiple levels of biological organization from the genome to receptor interactions and cellular processes leading to phenotypic changes such as developmental malformations. Finally, there is a growing recognition of the links between human and wildlife health and the need for approaches that allow for assessment of real world multi-chemical exposures. The zebrafish is poised to be an important model in bridging these two conventionally separate areas of toxicology and characterizing the biological effects of chemical mixtures that could augment its role in sustainable chemistry.

Fish gut-liver immunity during homeostasis or inflammation revealed by integrative transcriptome and proteome studies

The gut-associated lymphoid tissue, connected with liver via bile and blood, constructs a local immune environment of both defense and tolerance. The gut-liver immunity has been well-studied in mammals, yet in fish remains largely unknown, even though enteritis as well as liver and gallbladder syndrome emerged as a limitation in aquaculture. In this study, we performed integrative bioinformatic analysis for both transcriptomic (gut and liver) and proteomic (intestinal mucus and bile) data, in both healthy and infected tilapias. We found more categories of immune transcripts in gut than liver, as well as more adaptive immune in gut meanwhile more innate in liver. Interestingly reduced differential immune transcripts between gut and liver upon inflammation were also revealed. In addition, more immune proteins in bile than intestinal mucus were identified. And bile probably providing immune effectors to intestinal mucus upon inflammation was deduced. Specifically, many key immune transcripts in gut or liver as well as key immune proteins in mucus or bile were demonstrated. Accordingly, we proposed a hypothesized profile of fish gut-liver immunity, during either homeostasis or inflammation. Current data suggested that fish gut and liver may collaborate immunologically while keep homeostasis using own strategies, including potential unique mechanisms.

TBBPA exposure during a sensitive developmental window produces neurobehavioral changes in larval zebrafish

Tetrabromobisphenol A (TBBPA), one of the most widely used brominated flame retardants (BFRs), is a ubiquitous contaminant in the environment and in the human body. This study demonstrated that zebrafish embryos exposed to TBBPA during a sensitive window of 8-48 h post-fertilization (hpf) displayed morphological malformations and mortality. Zebrafish exposed exclusively between 48 and 96 hpf were phenotypically normal. TBBPA was efficiently absorbed and accumulated in zebrafish embryos, but was eliminated quickly when the exposure solution was removed. Larval behavior assays conducted at 120 hpf indicated that exposure to 5 μM TBBPA from 8 to 48 hpf produced larvae with significantly lower average activity and speed of movement in the normal condition than in those exposed from 48 to 96 hpf. Specifically, 8-48 hpf-exposed larvae spent significantly less time in both activity bursts and gross movements compared to control or 48-96 hpf exposed larvae. Consistent with the motor deficits, TBBPA induced apoptotic cell death, delayed cranial motor neuron development, inhibited primary motor neuron development and loosed muscle fiber during the early developmental stages. To further explore TBBPA-induced developmental and neurobehavioral toxicity, RNA-Seq analysis was used to identify early transcriptional changes following TBBPA exposure. In total, 1969 transcripts were significantly differentially expressed (P < 0.05, FDR < 0.05, 1.5-FC) upon TBBPA exposure. Functional and pathway analysis of the TBBPA transcriptional profile identified biological processes involved in nerve development, muscle filament sliding and contraction, and extracellular matrix disassembly and organization changed significantly. In addition, TBBPA also led to an elevation in the expression of genes encoding uridine diphosphate glucuronyl transferases (ugt), which could affect thyroxine (T4) metabolism and subsequently lead to neurobehavioral changes. In summary, TBBPA exposure during a narrow, sensitive developmental window perturbs various molecular pathways and results in neurobehavioral deficits in zebrafish.

Identification of aryl hydrocarbon receptor signaling pathways altered in TCDD-treated red seabream embryos by transcriptome analysis

2,3,7,8-Tetrachlorodibenzo-p-dioxin (TCDD) induces a broad spectrum of toxic effects including craniofacial malformation and neural damage in fish embryos. These effects are mainly mediated by the aryl hydrocarbon receptor (AHR). However, the mode of action between TCDD-induced AHR activation and adverse outcomes is not yet understood. To provide a comprehensive picture of the AHR signaling pathway in fish embryos exposed to TCDD, red seabream (Pagrus major) embryos were treated with graded concentrations of TCDD (0.3-37nM) in seawater, or with a mixture of TCDD and 500nM CH223191, an AHR-specific antagonist. The transcriptome of red seabream embryos was analyzed using a custom-made microarray with 6000 probes specifically prepared for this species. A Jonckheere-Terpstra test was performed to screen for genes that demonstrated altered mRNA expression levels following TCDD exposure. The signals of 1217 genes (as human homologs) were significantly altered in a TCDD concentration-dependent manner (q-value<0.2). Notably, the TCDD-induced alteration in mRNA expression was alleviated by co-exposure to CH223191, suggesting that the mRNA expression level of these genes was regulated by AHR. To identify TCDD-activated pathways, the microarray data were further subjected to gene set enrichment analysis (GSEA) and functional protein-protein interaction (PPI) network analysis. GSEA demonstrated that the effects of TCDD on sets of genes involved calcium, mitogen-activated protein kinase (MAPK), actin cytoskeleton, chemokine, T cell receptor, melanoma, vascular endothelial growth factor (VEGF), axon guidance, and renal cell carcinoma signaling pathways. These results suggest the hypotheses that TCDD induces immunosuppression via the calcium, MAPK, chemokine, and T cell receptor signaling pathways, neurotoxicity via VEGF signaling, and axon guidance alterations and teratogenicity via the dysregulation of the actin cytoskeleton and melanoma and renal cell carcinoma signaling pathways. Furthermore, the PPI network analysis indicated that the adverse outcome pathways of TCDD in the embryos might be propagated through several hub genes such as cell division control protein 42, phosphoinositide-3-kinase regulatory subunit 1, and guanine nucleotide-binding proteins. Understanding these pathways potentially allows for exploring the adverse outcome pathway of the effects of TCDD on the red seabream embryos.

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.

Genome-wide identification of suitable zebrafish Danio rerio reference genes for normalization of gene expression data by RT-qPCR

In this study, to systematically identify the most stably expressed genes for internal reference in zebrafish Danio rerio investigations, 37 D. rerio transcriptomic datasets (both RNA sequencing and microarray data) were collected from gene expression omnibus (GEO) database and unpublished data, and gene expression variations were analysed under three experimental conditions: tissue types, developmental stages and chemical treatments. Forty-four putative candidate genes were identified with the c.v. <0·2 from all datasets. Following clustering into different functional groups, 21 genes, in addition to four conventional housekeeping genes (eef1a1l1, b2m, hrpt1l and actb1), were selected from different functional groups for further quantitative real-time (qrt-)PCR validation using 25 RNA samples from different adult tissues, developmental stages and chemical treatments. The qrt-PCR data were then analysed using the statistical algorithm refFinder for gene expression stability. Several new candidate genes showed better expression stability than the conventional housekeeping genes in all three categories. It was found that sep15 and metap1 were the top two stable genes for tissue types, ube2a and tmem50a the top two for different developmental stages, and rpl13a and rp1p0 the top two for chemical treatments. Thus, based on the extensive transcriptomic analyses and qrt-PCR validation, these new reference genes are recommended for normalization of D. rerio qrt-PCR data respectively for the three different experimental conditions.

High Conservation in Transcriptomic and Proteomic Response of White Sturgeon to Equipotent Concentrations of 2,3,7,8-TCDD, PCB 77, and Benzo[a]pyrene

Adverse effects associated with exposure to dioxin-like compounds (DLCs) are mediated primarily through activation of the aryl hydrocarbon receptor (AHR). However, little is known about the cascades of events that link activation of the AHR to apical adverse effects. Therefore, this study used high-throughput, next-generation molecular tools to investigate similarities and differences in whole transcriptome and whole proteome responses to equipotent concentrations of three agonists of the AHR, 2,3,7,8-TCDD, PCB 77, and benzo[a]pyrene, in livers of a nonmodel fish, the white sturgeon (Acipenser transmontanus). A total of 926 and 658 unique transcripts were up- and down-regulated, respectively, by one or more of the three chemicals. Of the transcripts shared by responses to all three chemicals, 85% of up-regulated transcripts and 75% of down-regulated transcripts had the same magnitude of response. A total of 290 and 110 unique proteins were up- and down-regulated, respectively, by one or more of the three chemicals. Of the proteins shared by responses to all three chemicals, 70% of up-regulated proteins and 48% of down-regulated proteins had the same magnitude of response. Among treatments there was 68% similarity between the global transcriptome and global proteome. Pathway analysis revealed that perturbed physiological processes were indistinguishable between equipotent concentrations of the three chemicals. The results of this study contribute toward more completely describing adverse outcome pathways associated with activation of the AHR.

Retene causes multifunctional transcriptomic changes in the heart of rainbow trout (Oncorhynchus mykiss) embryos

Fish are particularly sensitive to aryl hydrocarbon receptor (AhR)-mediated developmental toxicity. The molecular mechanisms behind these adverse effects have remained largely unresolved in salmonids, and for AhR-agonistic polycyclic aromatic hydrocarbons (PAHs). This study explored the cardiac transcriptome of rainbow trout (Oncorhynchus mykiss) eleuteroembryos exposed to retene, an AhR-agonistic PAH. The embryos were exposed to retene (nominal concentration 32 μg/L) and control, their hearts were collected before, at and after the onset of the visible signs of developmental toxicity, and transcriptomic changes were studied by microarray analysis. Retene up- or down-regulated 122 genes. The largest Gene Ontology groups were signal transduction, transcription, apoptosis, cell growth, cytoskeleton, cell adhesion/mobility, cardiovascular development, xenobiotic metabolism, protein metabolism, lipid metabolism and transport, and amino acid metabolism. Together these findings suggest that retene affects multiple signaling cascades in the heart of rainbow trout embryos, and potentially disturbs processes related to cardiovascular development and function.

Chronic TCDD exposure results in the dysregulation of gene expression in splenic B-lymphocytes and in the impairments in T-cell and B-cell differentiation in mouse model

2,3,7,8-Tetrachlorodibenzo-p-dioxin (TCDD) exposure in humans is associated with marked immune suppressions and increased incidence of lymphoblastic diseases. To elucidate mechanisms of impairments in humoral immune responses, we used a murine model. Following a 20-week administration of low doses of TCDD, we observed severely reduced antibody titers, dramatically decreased number of splenic Th1 and Th2 cells and an increase in CD19(+) B cells. Transcriptional profiling of CD19(+) B cells showed that markers of pre-B cells were significantly elevated, indicating delayed B cell maturation. These changes in B cells were accompanied by decreases of T helper cell numbers and reduced IgM and IgG titers. A transcriptome analysis of splenic B cells followed by Ingenuity Pathway Analysis (IPA) revealed a set of differentially expressed genes known to play roles in tumorigenesis, cell-proliferation and cell-migration. The most up-regulated transcript gene was Eph receptor A2 (EphA2), a known oncogene, and the most down-regulated transcript was ZBTB16 that codes for a negative transcriptional regulator important in epigenetic chromatin remodeling. IPA identified cAMP-responsive element modulator (CREM) and cAMP-responsive element binding protein 1 (CREB1) as top upstream regulators. Consistently, a MAPPER promoter database analysis showed that all top dysregulated genes had CREM and/or CREB1 binding sites in their promoter regions. In summary, our data showed that chronic TCDD exposure in mice caused suppressed humoral immunity accompanied with profound dysregulation of gene expression in splenic B-lymphocytes, likely through cAMP-dependent pathways. This dysregulation resulted in impairments in T-cell and B-cell differentiation and activation of the tumorigenic transcription program.

Identification of Putative Nuclear Receptors and Steroidogenic Enzymes in Murray-Darling Rainbowfish (Melanotaenia fluviatilis) Using RNA-Seq and De Novo Transcriptome Assembly

Murray-Darling rainbowfish (Melanotaenia fluviatilis [Castelnau, 1878]; Atheriniformes: Melanotaeniidae) is a small-bodied teleost currently under development in Australasia as a test species for aquatic toxicological studies. To date, efforts towards the development of molecular biomarkers of contaminant exposure have been hindered by the lack of available sequence data. To address this, we sequenced messenger RNA from brain, liver and gonads of mature male and female fish and generated a high-quality draft transcriptome using a de novo assembly approach. 149,742 clusters of putative transcripts were obtained, encompassing 43,841 non-redundant protein-coding regions. Deduced amino acid sequences were annotated by functional inference based on similarity with sequences from manually curated protein sequence databases. The draft assembly contained protein-coding regions homologous to 95.7% of the complete cohort of predicted proteins from the taxonomically related species, Oryzias latipes (Japanese medaka). The mean length of rainbowfish protein-coding sequences relative to their medaka homologues was 92.1%, indicating that despite the limited number of tissues sampled a large proportion of the total expected number of protein-coding genes was captured in the study. Because of our interest in the effects of environmental contaminants on endocrine pathways, we manually curated subsets of coding regions for putative nuclear receptors and steroidogenic enzymes in the rainbowfish transcriptome, revealing 61 candidate nuclear receptors encompassing all known subfamilies, and 41 putative steroidogenic enzymes representing all major steroidogenic enzymes occurring in teleosts. The transcriptome presented here will be a valuable resource for researchers interested in biomarker development, protein structure and function, and contaminant-response genomics in Murray-Darling rainbowfish.

Associations between DNA methylation and schizophrenia-related intermediate phenotypes - a gene set enrichment analysis

Multiple genetic approaches have identified microRNAs as key effectors in psychiatric disorders as they post-transcriptionally regulate expression of thousands of target genes. However, their role in specific psychiatric diseases remains poorly understood. In addition, epigenetic mechanisms such as DNA methylation, which affect the expression of both microRNAs and coding genes, are critical for our understanding of molecular mechanisms in schizophrenia. Using clinical, imaging, genetic, and epigenetic data of 103 patients with schizophrenia and 111 healthy controls of the Mind Clinical Imaging Consortium (MCIC) study of schizophrenia, we conducted gene set enrichment analysis to identify markers for schizophrenia-associated intermediate phenotypes. Genes were ranked based on the correlation between DNA methylation patterns and each phenotype, and then searched for enrichment in 221 predicted microRNA target gene sets. We found the predicted hsa-miR-219a-5p target gene set to be significantly enriched for genes (EPHA4, PKNOX1, ESR1, among others) whose methylation status is correlated with hippocampal volume independent of disease status. Our results were strengthened by significant associations between hsa-miR-219a-5p target gene methylation patterns and hippocampus-related neuropsychological variables. IPA pathway analysis of the respective predicted hsa-miR-219a-5p target genes revealed associated network functions in behavior and developmental disorders. Altered methylation patterns of predicted hsa-miR-219a-5p target genes are associated with a structural aberration of the brain that has been proposed as a possible biomarker for schizophrenia. The (dys)regulation of microRNA target genes by epigenetic mechanisms may confer additional risk for developing psychiatric symptoms. Further study is needed to understand possible interactions between microRNAs and epigenetic changes and their impact on risk for brain-based disorders such as schizophrenia.

Zebrafish as animal model for aquaculture nutrition research

The aquaculture industry continues to promote the diversification of ingredients used in aquafeed in order to achieve a more sustainable aquaculture production system. The evaluation of large numbers of diets in aquaculture species is costly and requires time-consuming trials in some species. In contrast, zebrafish (Danio rerio) can solve these drawbacks as an experimental model, and represents an ideal organism to carry out preliminary evaluation of diets. In addition, zebrafish has a sequenced genome allowing the efficient utilization of new technologies, such as RNA-sequencing and genotyping platforms to study the molecular mechanisms that underlie the organism’s response to nutrients. Also, biotechnological tools like transgenic lines with fluorescently labeled neutrophils that allow the evaluation of the immune response in vivo, are readily available in this species. Thus, zebrafish provides an attractive platform for testing many ingredients to select those with the highest potential of success in aquaculture. In this perspective article aspects related to diet evaluation in which zebrafish can make important contributions to nutritional genomics and nutritional immunity are discussed.