Next Article in Journal
Cadmium Caused Different Toxicity to Photosystem I and Photosystem II of Freshwater Unicellular Algae Chlorella pyrenoidosa (Chlorophyta)
Next Article in Special Issue
Development and Applications of a Zebrafish (Danio rerio) CYP1A-Targeted Monoclonal Antibody (CRC4) with Reactivity across Vertebrate Taxa: Evidence for a Conserved CYP1A Epitope
Previous Article in Journal
Ambient Air Pollution and Risk for Stroke Hospitalization: Impact on Susceptible Groups
Previous Article in Special Issue
Developmental, Behavioral and Transcriptomic Changes in Zebrafish Embryos after Smoke Dye Exposure
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
Toxics
Volume 10
Issue 7
10.3390/toxics10070351
toxics-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

share Share announcement Help format_quote Cite question_answer Discuss in SciProfiles thumb_up
...
Endorse textsms
...
Comment
first_page
settings
Order Article Reprints
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Evaluating Phenotypic and Transcriptomic Responses Induced by Low-Level VOCs in Zebrafish: Benzene as an Example

by
Chia-Chen Wu
1,
Jessica R. Blount
2,
Alex Haimbaugh
2,3,
Samantha Heldman
2,3,
Jeremiah N. Shields
2 and
Tracie R. Baker
1,2,3,*
1
Department of Global and Environmental Health, University of Florida, Gainesville, FL 32610, USA
2
Institute of Environmental Health Sciences, Wayne State University, Detroit, MI 48202, USA
3
Department of Pharmacology, School of Medicine, Wayne State University, Detroit, MI 48202, USA
*
Author to whom correspondence should be addressed.
Toxics 2022, 10(7), 351; https://doi.org/10.3390/toxics10070351
Submission received: 15 May 2022 / Revised: 23 June 2022 / Accepted: 24 June 2022 / Published: 27 June 2022
(This article belongs to the Special Issue Fish Models for Human Toxicology)
Browse Figures
Versions Notes

Abstract

:
Urban environments are plagued by complex mixtures of anthropogenic volatile organic compounds (VOCs), such as mixtures of benzene, toluene, ethylene, and xylene (BTEX). Sources of BTEX that drive human exposure include vehicle exhaust, industrial emissions, off-gassing of building material, as well as oil spillage and leakage. Among the BTEX mixture, benzene is the most volatile compound and has been linked to numerous adverse health outcomes. However, few studies have focused on the effects of low-level benzene on exposure during early development, which is a susceptible window when hematological, immune, metabolic, and detoxification systems are immature. In this study, we used zebrafish to conduct a VOC exposure model and evaluated phenotypic and transcriptomic responses following 0.1 and 1 ppm benzene exposure during the first five days of embryogenesis (n = 740 per treatment). The benzene body burden was 2 mg/kg in 1 ppm-exposed larval zebrafish pools and under the detection limit in 0.1 ppm-exposed fish. No observable phenotypic changes were found in both larvae except for significant skeletal deformities in 0.1 ppm-exposed fish (p = 0.01) compared with unexposed fish. Based on transcriptomic responses, 1 ppm benzene dysregulated genes that were implicated with the development of hematological system, and the regulation of oxidative stress response, fatty acid metabolism, immune system, and inflammatory response, including apob, nfkbiaa, serpinf1, foxa1, cyp2k6, and cyp2n13 from the cytochrome P450 gene family. Key genes including pik3c2b, pltp, and chia.2 were differentially expressed in both 1 and 0.1 ppm exposures. However, fewer transcriptomic changes were induced by 0.1 ppm compared with 1 ppm. Future studies are needed to determine if these transcriptomic responses during embryogenesis have long-term consequences at levels equal to or lower than 1 ppm.

1. Introduction

Benzene, toluene, ethylbenzene, and xylene (BTEX) is a common volatile organic compound (VOC) mixture associated with the petrochemical industry. They are common constituents of crude oil and thus are important raw materials or intermediates in the manufacture of many chemicals. They are also byproducts of incomplete combustion and can be released into the air during oil and gas operations. Despite legislation and enforcement, anthropogenic sources such as vehicle exhaust, gas emissions from manufacturing and chemical facilities, and spills or leaks involving gasoline and diesel fuel still occur and can lead to ambient levels ranging from 0.1 to 100 µg/m3 in industrialized urban regions [1,2,3,4]. Environmental BTEX contribute significantly to poor air quality and pose a public health hazard across the globe [5,6].
Among the BTEX mixture compounds, benzene is the most volatile compound and has been linked to numerous adverse health outcomes that involve neurological, respiratory, cardiovascular, endocrine, renal, hepatic, immunological, and hematological systems [7]. Developmental exposure to benzene, particularly, needs more attention since it has been associated with abnormal birth outcomes in humans at non-occupational settings [8,9]. Emerging epidemiological studies indicate that children exposed to environmental benzene have elevated abnormalities in blood cells formation and hematologic functions compared with unexposed children [10,11,12]. The mechanisms of benzene-induced toxicity are not completely understood, but are believed to be associated with the metabolites of benzene, such as benzene oxide, phenol, hydroquinone, and catechol [13]. These metabolites cause damage in hematopoietic cells, lead to chromosomal aberration, produce oxidative stress, damage the immune system, and alter gene expression and epigenetic regulation [14,15,16]. Children and fetuses are more susceptible to low-level benzene exposure than adults due to ongoing development of host defenses, detoxification mechanisms, and hematopoietic and metabolic systems [17].
The lack of information on effects of developmental exposure to low-level benzene necessitates the use of a high-throughput animal model to address exposure concerns for this pervasive chemical. Due to the short generation times, high fecundity, and ease of husbandry, zebrafish is an excellent model for developmental toxicity screening for VOC. Since zebrafish share approximately 70% homologous genes with humans, results from toxicology assays using this vertebrate model are translatable to human health and disease outcomes [18]. In addition, zebrafish hematological and immune systems are similar to mammals [19,20], and the molecular pathways governing developmental hematopoiesis are highly conserved between zebrafish and mammals [21]. The innate immune system matures during zebrafish embryonic development and the full adaptive immunity is functional after 4–6 weeks postfertilization. Larval zebrafish life stages are suitable for studying the development of innate immunity and the regulation of the inflammatory response [22].
In this study, we evaluate low levels of benzene exposure (0.1 and 1 ppm) during the first five days of embryogenesis to identify the potential-induced phenotypic and transcriptomic changes and their associations with diseases and pathways implicated in human benzene exposure. To maintain consistent and rigorous control over the benzene concentrations during the exposure period, chemical exposures were conducted in sealed containers with daily water changes including renewal of the exposure chemicals. The benzene levels in water and whole larvae body were monitored. The results of this study will provide molecular signatures of benzene exposure in zebrafish that will guide and add dimension to the interpretation of VOC exposure studies in mammalian and human tissues and populations.

2. Methods

2.1. Fish Husbandry

The zebrafish used in this study were wild-type AB strains kept in 27–30 °C reverse osmosis water buffered with 60 mg/L Instant Ocean Salts in a recirculating system (Aquaneering Inc., San Diego, CA, USA) on a 14 h:10 h light/dark cycle. Twice daily, adult fish were fed Aquatox Fish Diet (Zeigler, PA, USA) and spirulina flakes (Zeigler, PA, USA). To obtain embryos for VOC exposure experiment, adult fish were housed in spawning tanks with a sex ratio of two females: one male. Four hours post fertilization (hpf), embryos were collected, cleaned with 58 ppm bleach (diluted from Chlorox bleach, the Clorox Co., Oakland, CA, USA) for 5 min, and rinsed with egg water that was made up of 600 mg/L Instant Ocean salts (Spectrum Brands, Madison, WI, USA) in RO water. The clean eggs were used for the exposure experiment. The Institutional Animal Care and Use Committee at Wayne State University approved these zebrafish protocols on 14 October 2020 (protocol #IACUC-19-02-0938) in accordance with the Guide for the Care and Use of Laboratory Animals [23].

2.2. Chemical Exposure

2.2.1. The Stock Chemicals and Dilutions

The 99.8% benzene stock (Sigma Aldrich, St. Louis, MO, USA) was serially diluted in a septum-sealed amber vial with 40 mL egg water to create 998 ppm and 99.8 ppm benzene stock dilutions. To minimize benzene loss through volatilization, benzene stock dilutions were mixed fresh daily, and caps were never removed from septum vials during or after benzene injection into egg water.

2.2.2. Exposures

Embryonic benzene exposures were conducted in 30 mL glass vials with septum caps. Each vial contained 20 embryos in 20 mL of egg water per exposure level until 5 days post fertilization (dpf). To inject benzene stock dilution into the glass vials, 20 μL of egg water, the amount equal to the stock dilution about to be added, was taken out to avoid pressure imbalance. Then, 20 μL, 998 ppm, and 99.8 ppm benzene stock dilution, or egg water, was injected through the sealed septum cap directly into the egg water with a gas tight syringe to reach 1, 0.1, and 0 ppm concentrations. Vials were gently inverted several times to mix and then incubated at 28 °C.
Over the course of five days, daily deaths and hatch rates were recorded. Dead embryos were removed from the vial. Water change was conducted in each vial daily by removing 20 mL exposure solution, replenishing 20 mL fresh egg water, and injecting benzene stock solutions or egg water as on day 0. To end chemical exposure on 5 dpf, larval fish were collected, rinsed three times with egg water, and divided into 24-well plates for subsequent behavioral analysis, abnormality, and mortality screening. The exposure experiment was repeated three times with the same replicate number of embryos per concentration. Embryos of each exposure experiment were obtained from 13 to 18 spawns of different parental stocks in our fish husbandry. A total of 740 embryos were used per treatment.

2.3. Measurements of Body Burden

To measure the benzene uptake, 2 mL water samples, 50 embryos, and approximately 50 larval fish were collected and processed to measure benzene concentration on 1 and 5 dpf. The water samples from all three exposure experiments were collected before water change. The embryos and larval fish were sampled from two exposure experiments. Embryos or fish were placed into glass vials and euthanized by cooling the vials on ice for 30 min. Egg water was then removed completely, and methanol was added to limit benzene volatilization. Fish or eggs were weighed and quickly transferred into methanol to 2 mL screw-top vials. These vials were then sonicated for 30 min in a cup horn sonicator with the frequency of 2 s on and 1 s off at 90% amplitude in 4 °C water bath. The resultant samples were then sent to Ann Arbor Technical Services (Ann Arbor, MI, USA) and analyzed by purge and trap GC/MS according to the USEPA method 8260C [24].

2.4. Abnormality and Mortality Screening

All zebrafish embryos were screened for mortality daily throughout the 5-day exposure period. On day 5, the following morphological abnormalities were assessed for the live embryos: numbers of unhatched embryos, skeletal deformities, improperly inflated swim bladder, yolk sac edema, and heart edema. Chi-square analysis and Bonferroni pairwise adjustment were conducted to compare if the morphological abnormalities differ from the control treatment.

2.5. Behavioral Screening

On day 8, fish with swim bladders and without morphological abnormalities were selected for behavioral testing. For each exposure concentration, 32 larval fish were divided into 24-well plates, with each well housing a single fish and 2 mL of fish water. The plated fish were placed into a Noldus DanioVision Observation Chamber (Noldus Information Technology, Wageningen, The Netherlands) with a constant temperature of 28.5 ± 0.5 °C. In the observation chamber, larval movements were tracked over 24 min, with alternating light-dark cycles lasting 3 min each. Movement data, including distance traveled, velocity, turn angle, and angular velocity were integrated every 30 s. All behavioral assays took place between 14:00 and 22:00.

2.6. RNA Isolation

Larval fish were quickly euthanized with tricaine methanesulfonate solution and then pooled in groups of five in RNALater. RNA was isolated using the RNeasy Lipid Tissue Mini Kit (Qiagen, Hilden, Germany) according to manufacturer’s instructions. Briefly, the frozen larval pool was homogenized in 300 μL of Qiazol lysis reagent (Qiagen, Hilden, Germany) with a sterile pestle and motor (Kimble-Kontes, Vineland, NJ, USA) until no tissues were seen. Then 700 μL of Qiazol reagent (Qiagen, Hilden, Germany) and subsequently 200 μL of chloroform (Fisher Scientific, Waltham, MA, USA) were added to yield RNA partitions at the upper aqueous phase. The aqueous phase was collected and mixed with the same amount of ethanol for appropriate binding conditions. The solutions were applied to the RNeasy Mini spin column (Qiagen, Hilden, Germany), purified, and eluted in 30 μL of nuclease-free water following the manufacturer’s protocol. The concentrations of eluted RNA were measured using a Qubit 3.0 Fluorometer (Invitrogen, CA, USA). In total, 15–20 larvae per treatment were used for RNA isolation.

2.7. Transcriptomic Analysis

From isolated RNA, 3′ mRNA-seq libraries were prepared using QuantSeq 3′ mRNA-Seq Library Prep Kit FWD for Illumina (Lexogen, Vienna, Austria). Samples were normalized to 40 ng/μL (total input of 200 ng in 5 μL) and amplified at 17 cycles. The sample concentrations and qualities of each sample were determined using a Qubit® 3.0 Fluorometer and Qubit® dsDNA Broad Range Assay Kit (Invitrogen, Carlsbad, CA, USA). The libraries were run on an Agilent TapeStation 2200 (Agilent Technologies, Santa Clara, CA, USA) for quality control and sequenced on a NovaSeq 6000 (Illumina, CA, USA) with single-end 75 bp reads. Reads were aligned to D. rerio (Build danRer10) using the BlueBee Genomics Platform (BlueBee, Rijswijk, The Netherlands). Differential gene expression between the control and exposure treatment was evaluated using DEseq2 (available through GenePattern; Broad Institute, Cambridge, MA, USA). Genes with significant changes in expression, as defined by absolute log2 fold change value ≥ 0.75 and p value < 0.05, were uploaded into the Ingenuity Pathway Analysis software (IPA; QIAGEN Bioinformatics, Redwood City, CA, USA) for analysis using Ensembl gene ID as identifiers. The significant biological process and canonical pathways were defined as those with p values < 0.5, (or −log (p-value) ≥ 1.3) based on the right-tailed Fisher’s Exact Test. The z-score was generated through an IPA software to predict the activation (≥2) or inhibition (≤−2) of a pathway.

3. Results

We measured the benzene levels in water and larval tissues during the exposure period. At the end of day-1 exposure (before the water change), mean benzene levels in the water phase were 0.74 ± 0.3 and 0.14 ± 0.05 ppm for 1 and 0.1 ppm benzene exposures, respectively. After 5-day exposure, the mean benzene body burden was 2.0 mg/kg in the 1 ppm-exposed larval fish pool, and below the detection limit (<0.064 mg/kg) in 0.1 ppm-exposed fishes. For phenotypic responses, neither 0.1 ppm nor 1 ppm benzene exposures affected the hatch rate or mortality compared with controls (Figure 1). On day 5, surviving fish were assessed for morphological abnormalities including skeletal deformities, underinflated swim bladders, and yolk sac and heart edemas. Only the 0.1 ppm exposure group showed twice as many skeletal deformities (26%, p = 0.01) and heart edema (21%, p = 0.06) than the control (11% and 10%, respectively). No other significant morphological differences were observed between 1 ppm-exposed and unexposed fish (Figure 1). In order to investigate the potential neurological effects of benzene exposure during development, we performed behavioral assays that track larval zebrafish movement during periods of light/dark stimulus. We observed no differences between control fish and those exposed to 1 ppm or 0.1 ppm benzene (Figure 1).
The transcriptomic analysis uncovered 239 and 49 differentially expressed genes (DEGs) in fish exposed to 1 and 0.1 ppm benzene, respectively, with 15 transcripts overlapping between the two exposure groups (Figure 2, Table S1). All overlapping transcripts expressed in the same direction. For 1 and 0.1 ppm exposures 152 and 37 transcripts were upregulated, while 86 and 12 transcripts were downregulated, respectively. These DEGs were implicated in disorders in gastrointestinal, endocrine, reproductive, and hematological systems (Table S2). The top associated networks and canonical pathways of these dysregulated genes were involved with stress response, inflammation, and the development of the hematological and immunological system (Figure 2b). Only 1 ppm exposure was predicted to inactivate or activate some biological processes (Table S3), including decreased cognition (z = −2.4), and learning ability (z = −2.2).
Table 1 shows the DEGs of interest that were either dysregulated in both benzene concentrations or involved with multiple pathways. The DEGs shared between two benzene exposure levels included pik3c2b, mep1a.2, chrnb2, and pltp. We found the upregulation of pik3c2b involved most of the implicated canonical pathways. Pik3c2b encodes phosphatidylinositol kinases and is responsible for phosphoinositide metabolism in vertebrate tissues, which is essential for membrane trafficking and intracellular signaling cascades in numerous biological processes [25,26]. Another overexpressed gene, mep1a.2 is orthologous to human mep1a and encodes alpha meprin A that has metalloendopeptidase activity. In our study, both 1 and 0.1 ppm benzene upregulated mep1a.2 1.7-fold higher than the control. The dysregulation of mep1a.2 was predicted to involve intestinal inflammatory response [27] and angiogenesis effect [28]. Chrnb2, encoding for beta 2 nicotinic cholinergic receptor, has been implicated with impaired cognition and learning ability following 1 ppm exposure. For both benzene exposures, pathway analysis implies chrnb2 can be involved with AMPK signaling and inflammatory response. Another shared gene that is involved with innate immune response is pltp, encoding phospholipid transfer protein (PLTP), which can modulate the production of proinflammatory cytokines [29]. Some shared genes were not involved with top canonical pathways but play important roles in embryogenesis, such as chia.2 and taf4a. Chia.2, one of the genes from chitinase gene family that is essential for zebrafish development [30], was upregulated approximately two-fold higher than control following both exposures. Their encoding proteins also regulate type 2 immune responses and can be triggered during allergic rhinitis in mammals [31]. Another gene that is associated with the production of blood cells and epithelial cells in liver during embryogenesis, taf4a, benzene downregulated the expression compared with the control [32,33]. For 0.1 ppm exposure, mapk11 was associated with most canonical pathways, which is involved in MAPK signaling in various cellular processes including embryogenesis and innate immunity [34].
Transcriptomic changes following 1 ppm benzene exposure were associated with detoxification processes, hematological development, and inflammation response. We detected upregulation of two genes from the cytochrome P450 (CYP) family commonly involved in the biotransformation of xenobiotics, cyp2k6 (orthologous to human cyp2w1) and cyp2n13 [35,36]. Many genes from the solute carrier (slc) family were also dysregulated, including downregulation of slc2a1a and slc9a5, as well as upregulation of slc23a3, slc27a2b, and slc26a3.2. These genes are implicated with hematological diseases or inflammatory responses. Several DEGs are implicated with hematopoietic neoplasm, such as the downregulation of antxr; it encodes the anthrax toxin receptor and can inhibit vascular development in zebrafish [37]. The Apob expression also can alter angiogenesis by mediating the expression of the receptors of the vascular endothelial growth factor [38]. Among the genes of interest that were exclusively differentially expressed following 1 ppm exposure, c3a.3 had the highest fold change. This gene is orthologous to human complement C3 and is essential for innate immunity [39]. Along with c3a.3, another complement component, c2 (si:ch1073-280e3.1), was also upregulated following 1 ppm exposure. Two cytokine-related genes, il15 and ifih1, encoding interleukin and interferon, respectively, are implicated with innate immune response and downregulated following 1 ppm exposure. A negative regulator of immunity responses [40], nfkbiaa, encoding the inhibitor of nuclear factor κB (NF-κB) transcription factors, was upregulated following 1 ppm benzene exposure. Some 1 ppm upregulated genes are responsible for multiple functions, such as foxa2 and serpinf1. Foxa2 encodes for transcription factors of the Forkhead box A family and is responsible for glucose and lipid homeostasis in liver [41] and adipose cells [42]. In addition, it is also involved in the development of embryonic structures [43], maturation of intestinal goblet cells [44], and modulation of inflammation response in skeletal and muscle cells [45]. The other gene, serpinf1, is also involved in hepatic lipid metabolism [46], and can be mediated by interleukin and vascular endothelial growth factors to inhibit angiogenesis [47].

4. Discussion

This is the first study to use zebrafish as a toxicological model for benzene exposure. Our results show that embryonic exposure to benzene levels as low as 0.1 ppm exerted skeletal abnormalities, and both 0.1 ppm and 1 ppm induced transcriptomic changes. Following a five-day exposure, up to 2 mg/kg benzene was detected in 1 ppm-exposed larval zebrafish pools, which can be viewed as the adsorbed dose of benzene in entire larvae tissues after metabolism. Cells in the embryonic zebrafish liver perform xenobiotic metabolism [48], and benzene metabolism can occur in the liver and bone marrow under exposure levels ranging from nanogram to milligram levels per kilogram body weight [14,49,50,51,52].
Among phenotypic responses, the only statistically significant parameter in our study was skeletal deformities following 0.1 ppm exposure. This morphological abnormality might be related to the dysregulation of mapk11, which is essential for controlling the initiation of dorsalization signals in zebrafish embryogenesis [53]. Among previous studies, skeletal malformations were found across species after prenatal exposure at levels larger than 47 ppm [54]. While skeletal abnormalities were not found in 1 ppm-exposed fish, more studies are needed to investigate if low-level benzene causes non-monotonic skeletal malformation during embryogenic stage.
The number of dysregulated genes implicated in the top canonical pathways indicates a monotonic dose–response relationship in transcriptomic changes. Among all, pik3c2b can be regarded as a key gene since it was induced by both exposure concentrations and was implicated in the most canonical pathways, including xenobiotic metabolism signaling, HIF1a signaling, AMPK signaling, FAK signaling, PI3k/Akt signaling, glucocorticoid receptor signaling, NRF2-mediated oxidative stress response, and G-protein coupled receptor signaling. Other genes from the PIK gene family, pik3cg, pik3r1, and pik3r2, were also differentially expressed in workers with chorionic benzene exposures [55]. For 0.1 ppm exposure specifically, mapk11 was one of the few DEGs that was implicated in more than one pathway, including xenobiotic metabolism signaling, AMPK signaling, antioxidant action of vitamin C, and glucocorticoid receptor signaling. For 1 ppm exposure, hsp90aa1, igf1, slc2a1, nfkbia, cyp2k6 and cyp2n13 were implicated in multiple pathways related to oxidative stress and detoxification process. The CYP protein family is known to oxidize benzene [56]. Common CYP genes, such as cyp2e1 and cyp2f1 were upregulated in workers exposed to low-level airborne benzene [57]. One of the CYP enzymes, CYP2B1, was found to be destructed by benzene metabolites, phenobarbital, and then be rescued by ascorbate (vitamin C) [56]. In our study, several genes relevant to antioxidant action of vitamin C were altered by 1 ppm exposure, such as the overexpression of nfkbia, plaat1, and slc23a3, as well as the downregulation of slc2a1.
We also found benzene exposure induced transcriptomic changes that were closely relevant to fatty acid oxidation, including two shared genes, pltp and mecr, and 1 ppm exclusive DEGs, apob, serpinf1, foxa2, c3, and hmgcra. Among the top canonical pathways, LXR/RXR function showed the most statistical significance (−log(p-value) = 3.1) and a trend toward activation (Z = 1.0). Zebrafish and mammals have conserved functions in LXR/RXR activation, which regulate the biosynthesis of fatty acid and cholesterol during embryogenesis [58]. In previous studies, benzene exposures increased the expression of fatty acid transport- and β-oxidation-related enzymes in mice [59]. Changes in fatty acid oxidation metabolism were detected in human plasma [60] and mice bone marrow cells following low-level occupational and high-level exposure, respectively. Since fatty acid oxidation is essential for the differentiation [61] and self-renewal [62] of hematopoietic stem cells, as well as the growth and maintenance of leukemic cells [63,64], they may be associated with the molecular mechanisms of benzene hematotoxicity.
Transcriptomic changes of 1 ppm benzene exposure were also related to the perturbation of the immune system during embryogenesis. Previous studies found benzene can dysregulate the cytokine network in vitro and in vivo, potentially through alteration of gene expression, or direct damage of hematopoietic progenitor cells [65,66,67]. In our study, one of the top networks is the crosstalk between PI3k/Akt signaling and PLTP modulation, which is implicated with the alteration of immune cell trafficking and inflammatory activation. The key gene of the oxidative stress response [68], pik3c2b, is also involved with PI3k/Akt signaling that regulates immune cell development, differentiation, and function [69]. Several metabolic-related DEGs, including the upregulation of apob and downregulation of pltp, are implicated with the activation of inflammatory response through PLTP modulation. The apob encoding protein is predicted to have low-density lipoprotein particle receptor binding activity and can transfer cholesterol, lipoprotein, and triglycerides. Elevated levels of lipoproteins can trigger the adhesion and migration of immune cells to the endothelial surface and the production of several plasma inflammatory markers, such as the tumor necrosis factor and interleukin [70]. The inflammatory response is strongly affected by circulation and disaggregation of lipoprotein-associated proteins, which can be modulated by PLTP [29]. Another mediator of cytokine homeostasis [45], foxA2, was upregulated following benzene exposure in our study. FoxA2 expression can also be upregulated by fatty acids and peroxisome proliferator-activated receptors in skeletal muscle in the mice model [45]. Nfkbia is another detoxification gene that is also involved with inflammatory response. It encodes the inhibitor of NF-κB signaling and thus can alter the expression of proinflammatory genes [71]. In our study, nfkbia was upregulated following 1 ppm exposure, and we observed that a cytokine gene, IL-15, was downregulated simultaneously.
Our study shows that benzene exposure of 1 ppm can induce transcriptomic alterations that are associated with oxidative stress response, fatty acid metabolism, immune system and inflammatory responses, and hematological system development. In addition, the number of dysregulated genes implicated for the top canonical pathways has a monotonic dose–response relationship with benzene levels. Previous studies in humans also show that benzene causes hematotoxic effects following occupational exposure at levels lower than 1 ppm [72,73,74]. Since the disruption of hematological system development is interlinked with benzene metabolism and associated inflammatory responses, our study indicates that zebrafish are great models to evaluate the effects of benzene exposure on transcriptomic changes. Future studies should be conducted to investigate if benzene-induced transcriptomic changes during embryogenesis have long-term consequences at levels equal to or lower than 1 ppm. Our results also demonstrate that the use of experimental zebrafish is reproducible and well-controlled model for VOC exposure studies. Since VOC contaminants often occur in the environment as complex mixtures, the appropriate integration of in vivo, in vitro, and in silico data with epidemiological studies is paramount for robust health risk assessments [75,76].

Supplementary Materials

The following supporting information can be downloaded at: https://www.mdpi.com/article/10.3390/toxics10070351/s1. Table S1: Differentially expressed genes due to benzene exposure. Table S2: Top IPA pathways affected by 1 or 0.1 ppm benzene exposure. Table S3: Estimated Z-score of biological process following 1 ppm benzene exposure.

Author Contributions

Conceptualization, T.R.B. and J.N.S.; methodology, J.R.B., J.N.S. and S.H.; validation, J.R.B., A.H. and S.H.; formal analysis, C.-C.W., J.R.B. and A.H.; investigation, C.-C.W., J.R.B. and A.H.; data curation, C.-C.W. and J.R.B.; writing—original draft preparation, C.-C.W. and J.R.B.; writing—review and editing, C.-C.W., J.R.B. and T.R.B.; visualization, T.R.B.; supervision, T.R.B.; funding acquisition, T.R.B. All authors have read and agreed to the published version of the manuscript.

Funding

Funding was provided by the National Institute of Environmental Health Sciences [R01 ES030722 to TRB, AH; P30 ES020957 to TRB].

Institutional Review Board Statement

All zebrafish use protocols were approved by the Institutional Animal Care and Use Committee at Wayne State University, according to the National Institutes Health Guide to the Care and Use of Laboratory Animals (Protocol 16-03-054; approved 4 August 2016).

Data Availability Statement

All data is published in this manuscript and Supplementary Materials.

Acknowledgments

We acknowledge Emily Crofts, Kim Bauman, and all members of the Warrior Aquatic, Translational, and Environmental Research (WATER) lab at Wayne State University for help with zebrafish care and husbandry. We would like to acknowledge the Wayne State University Applied Genomics Technology Center, including Katherine Gurdziel, for providing sequencing services and the use of Ingenuity Pathway Analysis Software.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. McKenzie, L.M.; Witter, R.Z.; Newman, L.S.; Adgate, J.L. Human Health Risk Assessment of Air Emissions from Development of Unconventional Natural Gas Resources. Sci. Total Environ. 2012, 424, 79–87. [Google Scholar] [CrossRef] [PubMed]
  2. Miller, L.; Xu, X.; Grgicak-Mannion, A.; Brook, J.; Wheeler, A. Multi-Season, Multi-Year Concentrations and Correlations amongst the BTEX Group of VOCs in an Urbanized Industrial City. Atmos. Environ. 2012, 61, 305–315. [Google Scholar] [CrossRef]
  3. Miri, M.; Rostami Aghdam Shendi, M.; Ghaffari, H.R.; Ebrahimi Aval, H.; Ahmadi, E.; Taban, E.; Gholizadeh, A.; Yazdani Aval, M.; Mohammadi, A.; Azari, A. Investigation of Outdoor BTEX: Concentration, Variations, Sources, Spatial Distribution, and Risk Assessment. Chemosphere 2016, 163, 601–609. [Google Scholar] [CrossRef] [PubMed]
  4. Duan, X.; Li, Y. Sources and Fates of BTEX in the General Environment and Its Distribution in Coastal Cities of China. J. Environ. Sci. Public Health 2017, 1, 86–106. [Google Scholar] [CrossRef] [Green Version]
  5. Miller, C.J.; Runge-Morris, M.; Cassidy-Bushrow, A.E.; Straughen, J.K.; Dittrich, T.M.; Baker, T.R.; Petriello, M.C.; Mor, G.; Ruden, D.M.; O’leary, B.F.; et al. A Review of Volatile Organic Compound Contamination in Post-Industrial Urban Centers: Reproductive Health Implications Using a Detroit Lens. Int. J. Environ. Res. Public Health 2020, 17, 8755. [Google Scholar] [CrossRef] [PubMed]
  6. Davidson, C.J.; Hannigan, J.H.; Bowen, S.E. Effects of Inhaled Combined Benzene, Toluene, Ethylbenzene, and Xylenes (BTEX): Toward an Environmental Exposure Model. Environ. Toxicol. Pharmacol. 2021, 81, 103518. [Google Scholar] [CrossRef]
  7. Bahadar, H.; Mostafalou, S.; Abdollahi, M. Current Understandings and Perspectives on Non-Cancer Health Effects of Benzene: A Global Concern. Toxicol. Appl. Pharmacol. 2014, 276, 83–94. [Google Scholar] [CrossRef] [Green Version]
  8. Slama, R.; Thiebaugeorges, O.; Goua, V.; Aussel, L.; Sacco, P.; Bohet, A.; Forhan, A.; Ducot, B.; Annesi-Maesano, I.; Heinrich, J. Maternal Personal Exposure to Airborne Benzene and Intrauterine Growth. Environ. Health Perspect. 2009, 117, 1313–1321. [Google Scholar] [CrossRef]
  9. Chen, D.; Cho, S.-I.; Chen, C.; Wang, X.; Damokosh, A.I.; Ryan, L.; Smith, T.J.; Christiani, D.C.; Xu, X. Exposure to Benzene, Occupational Stress, and Reduced Birth Weight. Occup. Environ. Med. 2000, 57, 661–667. [Google Scholar] [CrossRef] [Green Version]
  10. Lee, C.R.; Yoo, C.I.; Lee, J.H.; Kim, S.-R.; Kim, Y. Hematological Changes of Children Exposed to Volatile Organic Compounds Containing Low Levels of Benzene. Sci. Total Environ. 2002, 299, 237–245. [Google Scholar] [CrossRef]
  11. D’Andrea, M.A.; Reddy, G.K. Health Effects of Benzene Exposure among Children Following a Flaring Incident at the British Petroleum Refinery in Texas City. Pediatr. Hematol. Oncol. 2014, 31, 1–10. [Google Scholar] [CrossRef] [PubMed]
  12. D’Andrea, M.A.; Reddy, G.K. Adverse Health Effects of Benzene Exposure among Children Following a Flaring Incident at the British Petroleum Refinery in Texas City. Clin. Pediatr. 2016, 55, 219–227. [Google Scholar] [CrossRef] [PubMed]
  13. Rappaport, S.M.; Kim, S.; Lan, Q.; Vermeulen, R.; Waidyanatha, S.; Zhang, L.; Li, G.; Yin, S.; Hayes, R.B.; Rothman, N.; et al. Evidence That Humans Metabolize Benzene via Two Pathways. Environ. Health Perspect. 2009, 117, 946–952. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  14. Cooper, K.R.; And, V.K.; Snyder, R. Correlation of Benzene Metabolism and Histological Lesions in Rainbow Trout (Salmo Gairdneri). Drug Metab. Rev. 1984, 15, 673–696. [Google Scholar] [CrossRef] [PubMed]
  15. Khan, H.A. Benzene’s Toxicity: A Consolidated Short Review of Human and Animal Studies by HA Khan. Hum. Exp. Toxicol. 2007, 26, 677–685. [Google Scholar] [CrossRef]
  16. Spatari, G.; Allegra, A.; Carrieri, M.; Pioggia, G.; Gangemi, S. Epigenetic Effects of Benzene in Hematologic Neoplasms: The Altered Gene Expression. Cancers 2021, 13, 2392. [Google Scholar] [CrossRef]
  17. D’Andrea, M.A.; Reddy, G.K. Health Risks Associated with Benzene Exposure in Children: A Systematic Review. Glob. Pediatr. Health 2018, 5, 2333794X18789275. [Google Scholar] [CrossRef]
  18. Howe, K.; Clark, M.D.; Torroja, C.F.; Torrance, J.; Berthelot, C.; Muffato, M.; Collins, J.E.; Humphray, S.; McLaren, K.; Matthews, L.; et al. The Zebrafish Reference Genome Sequence and Its Relationship to the Human Genome. Nature 2013, 496, 498–503. [Google Scholar] [CrossRef] [Green Version]
  19. Trede, N.S.; Langenau, D.M.; Traver, D.; Look, A.T.; Zon, L.I. The Use of Zebrafish to Understand Immunity. Immunity 2004, 20, 367–379. [Google Scholar] [CrossRef] [Green Version]
  20. Ciau-Uitz, A.; Monteiro, R.; Kirmizitas, A.; Patient, R. Developmental Hematopoiesis: Ontogeny, Genetic Programming and Conservation. Exp. Hematol. 2014, 42, 669–683. [Google Scholar] [CrossRef]
  21. Robertson, A.L.; Avagyan, S.; Gansner, J.M.; Zon, L.I. Understanding the Regulation of Vertebrate Hematopoiesis and Blood Disorders—Big Lessons from a Small Fish. FEBS Lett. 2016, 590, 4016–4033. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  22. Stockhammer, O.W.; Zakrzewska, A.; Hegedûs, Z.; Spaink, H.P.; Meijer, A.H. Transcriptome Profiling and Functional Analyses of the Zebrafish Embryonic Innate Immune Response to Salmonella Infection. J. Immunol. 2009, 182, 5641–5653. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  23. National Research Council. Guide for the Care and Use of Laboratory Animals, 8th ed.; The National Academies Press: Washington, DC, USA, 2010. [Google Scholar]
  24. U.S. Environmental Protection Agency. Method 8260C Volatile Organic Compounds by Gas Chromatography/Mass Spectrometry (GC/MS); U.S. Environmental Protection Agency: Washington, DC, USA, 2006.
  25. Krauß, M.; Haucke, V. Phosphoinositide-metabolizing Enzymes at the Interface between Membrane Traffic and Cell Signalling. EMBO Rep. 2007, 8, 241–246. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  26. Di Paolo, G.; De Camilli, P. Phosphoinositides in Cell Regulation and Membrane Dynamics. Nature 2006, 443, 651–657. [Google Scholar] [CrossRef]
  27. Banerjee, S.; Oneda, B.; Yap, L.; Jewell, D.; Matters, G.; Fitzpatrick, L.; Seibold, F.; Sterchi, E.; Ahmad, T.; Lottaz, D.; et al. MEP1A Allele for Meprin A Metalloprotease Is a Susceptibility Gene for Inflammatory Bowel Disease. Mucosal. Immunol. 2009, 2, 220–231. [Google Scholar] [CrossRef] [Green Version]
  28. Schütte, A.; Hedrich, J.; Stöcker, W.; Becker-Pauly, C. Let It Flow: Morpholino Knockdown in Zebrafish Embryos Reveals a pro-Angiogenic Effect of the Metalloprotease Meprin A2. PLoS ONE 2010, 5, e8835. [Google Scholar] [CrossRef]
  29. Gautier, T.; Paul, C.; Deckert, V.; Desrumaux, C.; Klein, A.; Labbé, J.; Le Guern, N.; Athias, A.; Monier, S.; Hammann, A.; et al. Innate Immune Response Triggered by Triacyl Lipid a is Dependent on Phospholipid Transfer Protein (PLTP) Gene Expression. FASEB J. 2010, 24, 3544–3554. [Google Scholar] [CrossRef]
  30. Koch, B.E.V.; Stougaard, J.; Spaink, H.P. Spatial and Temporal Expression Patterns of Chitinase Genes in Developing Zebrafish Embryos. Gene Expr. Patterns 2014, 14, 69–77. [Google Scholar] [CrossRef] [Green Version]
  31. Cho, W.S.; Kim, T.H.; Lee, H.M.; Lee, S.H.; Lee, S.H.; Yoo, J.H.; Kim, Y.S.; Lee, S.H. Increased Expression of Acidic Mammalian Chitinase and Chitotriosidase in the Nasal Mucosa of Patients with Allergic Rhinitis. Laryngoscope 2010, 120, 870–875. [Google Scholar] [CrossRef]
  32. Alpern, D.; Langer, D.; Ballester, B.; Le Gras, S.; Romier, C.; Mengus, G.; Davidson, I. TAF4, a Subunit of Transcription Factor II D, Directs Promoter Occupancy of Nuclear Receptor HNF4A during Post-Natal Hepatocyte Differentiation. Elife 2014, 3, e03613. [Google Scholar] [CrossRef]
  33. Langer, D.; Martianov, I.; Alpern, D.; Rhinn, M.; Keime, C.; Dollé, P.; Mengus, G.; Davidson, I. Essential Role of the TFIID Subunit TAF4 in Murine Embryogenesis and Embryonic Stem Cell Differentiation. Nat. Commun. 2016, 7, 11063. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  34. Krens, S.F.G.; Spaink, H.P.; Snaar-Jagalska, B.E. Functions of the MAPK Family in Vertebrate-Development. FEBS Lett. 2006, 580, 4984–4990. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  35. Verbueken, E.; Bars, C.; Ball, J.S.; Periz-Stanacev, J.; Marei, W.F.A.; Tochwin, A.; Gabriëls, I.J.; Michiels, E.D.G.; Stinckens, E.; Vergauwen, L.; et al. From MRNA Expression of Drug Disposition Genes to in Vivo Assessment of CYP-Mediated Biotransformation during Zebrafish Embryonic and Larval Development. Int. J. Mol. Sci. 2018, 19, 3976. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  36. Nawaji, T.; Yamashita, N.; Umeda, H.; Zhang, S.; Mizoguchi, N.; Seki, M.; Kitazawa, T.; Teraoka, H. Cytochrome P450 Expression and Chemical Metabolic Activity before Full Liver Development in Zebrafish. Pharmaceuticals 2020, 13, 456. [Google Scholar] [CrossRef]
  37. Bolcome, R.E.; Sullivan, S.E.; Zeller, R.; Barker, A.P.; Collier, R.J.; Chan, J. Anthrax Lethal Toxin Induces Cell Death-Independent Permeability in Zebrafish Vasculature. Proc. Natl. Acad. Sci. USA 2008, 105, 2439–2444. [Google Scholar] [CrossRef] [Green Version]
  38. Avraham-Davidi, I.; Ely, Y.; Pham, V.N.; Castranova, D.; Grunspan, M.; Malkinson, G.; Gibbs-Bar, L.; Mayseless, O.; Allmog, G.; Lo, B. ApoB-Containing Lipoproteins Regulate Angiogenesis by Modulating Expression of VEGF Receptor 1. Nat. Med. 2012, 18, 967–973. [Google Scholar] [CrossRef] [Green Version]
  39. Zhang, S.; Cui, P. Complement System in Zebrafish. Dev. Comp. Immunol. 2014, 46, 3–10. [Google Scholar] [CrossRef]
  40. Chatterjee, B.; Banoth, B.; Mukherjee, T.; Taye, N.; Vijayaragavan, B.; Chattopadhyay, S.; Gomes, J.; Basak, S. Late-Phase Synthesis of IκBα Insulates the TLR4-Activated Canonical NF-ΚB Pathway from Noncanonical NF-ΚB Signaling in Macrophages. Sci. Signal. 2016, 9, ra120. [Google Scholar] [CrossRef] [Green Version]
  41. Zhang, L.; Rubins, N.E.; Ahima, R.S.; Greenbaum, L.E.; Kaestner, K.H. Foxa2 Integrates the Transcriptional Response of the Hepatocyte to Fasting. Cell Metab. 2005, 2, 141–148. [Google Scholar] [CrossRef] [Green Version]
  42. Wolfrum, C.; Shih, D.Q.; Kuwajima, S.; Norris, A.W.; Kahn, C.R.; Stoffel, M. Role of Foxa-2 in Adipocyte Metabolism and Differentiation. J. Clin. Investig. 2003, 112, 345–356. [Google Scholar] [CrossRef] [Green Version]
  43. Dal-Pra, S.; Thisse, C.; Thisse, B. FoxA Transcription Factors Are Essential for the Development of Dorsal Axial Structures. Dev. Biol. 2011, 350, 484–495. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  44. Lai, Y.-R.; Lu, Y.-F.; Lien, H.-W.; Huang, C.-J.; Hwang, S.-P.L. Foxa2 and Hif1ab Regulate Maturation of Intestinal Goblet Cells by Modulating Agr2 Expression in Zebrafish Embryos. Biochem. J. 2016, 473, 2205–2218. [Google Scholar] [CrossRef] [PubMed]
  45. Phua, W.W.T.; Tan, W.R.; Yip, Y.S.; Hew, I.D.; Wee, J.W.K.; Cheng, H.S.; Leow, M.K.S.; Wahli, W.; Tan, N.S. PPARβ/δ Agonism Upregulates Forkhead Box A2 to Reduce Inflammation in C2C12 Myoblasts and in Skeletal Muscle. Int. J. Mol. Sci. 2020, 21, 1747. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  46. Borg, M.L.; Andrews, Z.B.; Duh, E.J.; Zechner, R.; Meikle, P.J.; Watt, M.J. Pigment Epithelium-Derived Factor Regulates Lipid Metabolism via Adipose Triglyceride Lipase. Diabetes 2011, 60, 1458–1466. [Google Scholar] [CrossRef] [Green Version]
  47. Shellenberger, T.D.; Mazumdar, A.; Henderson, Y.; Briggs, K.; Wang, M.; Chattopadhyay, C.; Jayakumar, A.; Frederick, M.; Clayman, G.L. Headpin: A Serpin with Endogenous and Exogenous Suppression of Angiogenesis. Cancer Res. 2005, 65, 11501–11509. [Google Scholar] [CrossRef] [Green Version]
  48. Goessling, W.; Sadler, K.C. Zebrafish: An Important Tool for Liver Disease Research. Gastroenterology 2015, 149, 1361–1377. [Google Scholar] [CrossRef] [Green Version]
  49. Turteltaub, K.W.; Mani, C. Benzene Metabolism in Rodents at Doses Relevant to Human Exposure from Urban Air. Res. Rep. Health Eff. Inst. 2003, 113, 1–26. [Google Scholar]
  50. Valcke, M.; Krishnan, K. Assessing the Impact of the Duration and Intensity of Inhalation Exposure on the Magnitude of the Variability of Internal Dose Metrics in Children and Adults. Inhal. Toxicol. 2011, 23, 863–877. [Google Scholar] [CrossRef]
  51. Sammett, D.; Lee, E.W.; Kocsis, J.J.; Snyder, R. Partial Hepatectomy Reduces Both Metabolism and Toxicity of Benzene. J. Toxicol. Environ. Health Part A Curr. Issues 1979, 5, 785–792. [Google Scholar] [CrossRef]
  52. Sheets, P.L.; Carlson, G.P. Kinetic Factors Involved in the Metabolism of Benzene in Mouse Lung and Liver. J. Toxicol. Environ. Health Part A 2004, 67, 421–430. [Google Scholar] [CrossRef]
  53. Fujii, R.; Yamashita, S.; Hibi, M.; Hirano, T. Asymmetric P38 Activation in Zebrafish: Its Possible Role in Symmetric and Synchronous Cleavage. J. Cell Biol. 2000, 150, 1335–1347. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  54. Wilbur, S.B.; Keith, S.; Faroon, O.; Wohlers, D. Toxicological Profile for Benzene; U.S. Department of Health & Human Services: Washington, DC, USA, 2007. [Google Scholar] [CrossRef]
  55. Gao, A.; Yang, J.; Yang, G.; Niu, P.; Tian, L. Differential Gene Expression Profiling Analysis in Workers Occupationally Exposed to Benzene. Sci. Total Environ. 2014, 472, 872–879. [Google Scholar] [CrossRef] [PubMed]
  56. Gut, I.; Nedelcheva, V.; Soucek, P.; Stopka, P.; Tichavská, B. Cytochromes P450 in Benzene Metabolism and Involvement of Their Metabolites and Reactive Oxygen Species in Toxicity. Environ. Health Perspect. 1996, 104 (Suppl. S6), 1211–1218. [Google Scholar] [PubMed] [Green Version]
  57. Thomas, R.; Hubbard, A.E.; McHale, C.M.; Zhang, L.; Rappaport, S.M.; Lan, Q.; Rothman, N.; Vermeulen, R.; Guyton, K.Z.; Jinot, J.; et al. Characterization of Changes in Gene Expression and Biochemical Pathways at Low Levels of Benzene Exposure. PLoS ONE 2014, 9, e91828. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  58. Pinto, C.L.; Kalasekar, S.M.; McCollum, C.W.; Riu, A.; Jonsson, P.; Lopez, J.; Swindell, E.C.; Bouhlatouf, A.; Balaguer, P.; Bondesson, M.; et al. Lxr Regulates Lipid Metabolic and Visual Perception Pathways during Zebrafish Development. Mol. Cell. Endocrinol. 2016, 419, 29–43. [Google Scholar] [CrossRef] [Green Version]
  59. Sun, R.; Cao, M.; Zhang, J.; Yang, W.; Wei, H.; Meng, X.; Yin, L.; Pu, Y. Benzene Exposure Alters Expression of Enzymes Involved in Fatty Acid β-Oxidation in Male C3H/He Mice. Int. J. Environ. Res. Public Health 2016, 13, 1068. [Google Scholar] [CrossRef] [Green Version]
  60. Sun, R.; Xu, K.; Zhang, Q.; Jiang, X.; Man, Z.; Yin, L.; Zhang, J.; Pu, Y. Plasma Metabonomics Investigation Reveals Involvement of Fatty Acid Oxidation in Hematotoxicity in Chinese Benzene-Exposed Workers with Low White Blood Cell Count. Environ. Sci. Pollut. Res. 2018, 25, 32506–32514. [Google Scholar] [CrossRef]
  61. Ito, K.; Carracedo, A.; Weiss, D.; Arai, F.; Ala, U.; Avigan, D.E.; Schafer, Z.T.; Evans, R.M.; Suda, T.; Lee, C.-H. A PML–PPAR-δ Pathway for Fatty Acid Oxidation Regulates Hematopoietic Stem Cell Maintenance. Nat. Med. 2012, 18, 1350–1358. [Google Scholar] [CrossRef] [Green Version]
  62. Ito, K.; Suda, T. Metabolic Requirements for the Maintenance of Self-Renewing Stem Cells. Nat. Rev. Mol. Cell Biol. 2014, 15, 243–256. [Google Scholar] [CrossRef] [Green Version]
  63. Ricciardi, M.R.; Mirabilii, S.; Allegretti, M.; Licchetta, R.; Calarco, A.; Torrisi, M.R.; Foa, R.; Nicolai, R.; Peluso, G.; Tafuri, A. Targeting the Leukemia Cell Metabolism by the CPT1a Inhibition: Functional Preclinical Effects in Leukemias. Blood J. Am. Soc. Hematol. 2015, 126, 1925–1929. [Google Scholar] [CrossRef] [Green Version]
  64. Samudio, I.; Harmancey, R.; Fiegl, M.; Kantarjian, H.; Konopleva, M.; Korchin, B.; Kaluarachchi, K.; Bornmann, W.; Duvvuri, S.; Taegtmeyer, H. Pharmacologic Inhibition of Fatty Acid Oxidation Sensitizes Human Leukemia Cells to Apoptosis Induction. J. Clin. Investig. 2010, 120, 142–156. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  65. Wang, J.; Guo, X.; Chen, Y.; Zhang, W.; Ren, J.; Gao, A. Association between Benzene Exposure, Serum Levels of Cytokines and Hematological Measures in Chinese Workers: A Cross-Sectional Study. Ecotoxicol. Environ. Saf. 2021, 207, 111562. [Google Scholar] [CrossRef] [PubMed]
  66. Cassidy-Bushrow, A.E.; Burmeister, C.; Birbeck, J.; Chen, Y.; Lamerato, L.; Lemke, L.D.; Li, J.; Mor, G.; O’Leary, B.F.; Peters, R.M.; et al. Ambient BTEX Exposure and Mid-Pregnancy Inflammatory Biomarkers in Pregnant African American Women. J. Reprod. Immunol. 2021, 145, 103305. [Google Scholar] [CrossRef]
  67. Minciullo, P.L.; Navarra, M.; Calapai, G.; Gangemi, S. Cytokine Network Involvement in Subjects Exposed to Benzene. J. Immunol. Res. 2014, 2014, 937987. Available online: https://www.hindawi.com/journals/jir/2014/937987/ (accessed on 10 May 2022). [CrossRef] [PubMed] [Green Version]
  68. Koundouros, N.; Poulogiannis, G. Phosphoinositide 3-Kinase/Akt Signaling and Redox Metabolism in Cancer. Front. Oncol. 2018, 8, 160. [Google Scholar] [CrossRef]
  69. Okkenhaug, K. Signalling by the Phosphoinositide 3-Kinase Family in Immune Cells. Annu. Rev. Immunol. 2013, 31, 675–704. [Google Scholar] [CrossRef] [Green Version]
  70. Berbée, J.F.P.; Havekes, L.M.; Rensen, P.C.N. Apolipoproteins Modulate the Inflammatory Response to Lipopolysaccharide. J. Endotoxin Res. 2005, 11, 97–103. [Google Scholar] [CrossRef]
  71. Lawrence, T. The Nuclear Factor NF-KappaB Pathway in Inflammation. Cold Spring Harb. Perspect. Biol. 2009, 1, a001651. [Google Scholar] [CrossRef] [Green Version]
  72. Swaen, G.M.H.; van Amelsvoort, L.; Twisk, J.J.; Verstraeten, E.; Slootweg, R.; Collins, J.J.; Burns, C.J. Low Level Occupational Benzene Exposure and Hematological Parameters. Chem. Biol. Interact. 2010, 184, 94–100. [Google Scholar] [CrossRef]
  73. Schnatter, A.R.; Glass, D.C.; Tang, G.; Irons, R.D.; Rushton, L. Myelodysplastic Syndrome and Benzene Exposure among Petroleum Workers: An International Pooled Analysis. J. Natl. Cancer Inst. 2012, 104, 1724–1737. [Google Scholar] [CrossRef] [Green Version]
  74. Lan, Q.; Zhang, L.; Li, G.; Vermeulen, R.; Weinberg, R.S.; Dosemeci, M.; Rappaport, S.M.; Shen, M.; Alter, B.P.; Wu, Y.; et al. Hematotoxicity in Workers Exposed to Low Levels of Benzene. Science 2004, 306, 1774–1776. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  75. Batterman, S.; Su, F.-C.; Li, S.; Mukherjee, B.; Jia, C. Personal Exposure to Mixtures of Volatile Organic Compounds: Modeling and Further Analysis of the RIOPA Data. Res. Rep. Health. Eff. Inst. 2014, 181, 3. [Google Scholar]
  76. Hernández, A.F.; Tsatsakis, A.M. Human Exposure to Chemical Mixtures: Challenges for the Integration of Toxicology with Epidemiology Data in Risk Assessment. Food Chem. Toxicol. 2017, 103, 188–193. [Google Scholar] [CrossRef] [PubMed]
Toxics 10 00351 g001 550
Figure 1. Phenotypic responses following benzene exposure from 4 h post-fertilization (hpf) to 5 days post-fertilization (dpf). (a) Average mortality (left) and hatch rate (right) of live embryos over the course of exposure period (n = 740 per treatment); (b) percentage of embryos with morphological abnormalities at 5 dpf per treatment (n = 134, 149, 146 for 0, 0.1, and 1 ppm benzene); or (c) average distance moved by larval zebrafish during light and dark cycles (n = 32 per treatment). * Condition is significantly different from control (0 ppm; p < 0.05).
Figure 1. Phenotypic responses following benzene exposure from 4 h post-fertilization (hpf) to 5 days post-fertilization (dpf). (a) Average mortality (left) and hatch rate (right) of live embryos over the course of exposure period (n = 740 per treatment); (b) percentage of embryos with morphological abnormalities at 5 dpf per treatment (n = 134, 149, 146 for 0, 0.1, and 1 ppm benzene); or (c) average distance moved by larval zebrafish during light and dark cycles (n = 32 per treatment). * Condition is significantly different from control (0 ppm; p < 0.05).
Toxics 10 00351 g001
Toxics 10 00351 g002 550
Figure 2. Transcriptomic responses following benzene exposure from 4 h post-fertilization to 5 days post-fertilization (dpf). (a) Venn diagram of the number of differentially expressed genes (DEGs) upregulated or downregulated at 0.1, 1 ppm, and both; and (b) number of DEGs implicated in top-rank canonical pathways.
Figure 2. Transcriptomic responses following benzene exposure from 4 h post-fertilization to 5 days post-fertilization (dpf). (a) Venn diagram of the number of differentially expressed genes (DEGs) upregulated or downregulated at 0.1, 1 ppm, and both; and (b) number of DEGs implicated in top-rank canonical pathways.
Toxics 10 00351 g002
Table
Table 1. Dysregulated genes of benzene exposure based on gene expression fold changes relative to unexposed controls.
Table 1. Dysregulated genes of benzene exposure based on gene expression fold changes relative to unexposed controls.
Gene SymbolGene Name1 ppm0.1 ppm
Inflammatory response
pik3c2bPhosphatidylinositol-4-phosphate 3-kinase, catalytic subunit type 2 beta2.81.85
mep1a.2Meprin A, alpha (PABA peptide hydrolase), tandem duplicate 21.71.73
chia.2Acidic chitinase 22.011.90
chrnb2Cholinergic receptor, nicotinic, beta 20.590.54
pltpPhospholipid transfer protein0.570.58
nfkbiaaNuclear factor of kappa light polypeptide gene enhancer in B-cells inhibitor, alpha a1.70
apobApolipoprotein B1.70
si:ch1073-280e3.1Orthologous to human c2, complement c21.70
c3a.3Complement C3a, tandem duplicate 32.16
il15Interleukin 150.49
ifih1Interferon induced with helicase C domain 10.54
mapk11Mitogen-activated protein kinase 11 0.58
Hematological system development and function
acad11Acyl-CoA dehydrogenase family, member 113.20
antxr2ANTXR cell adhesion molecule 20.49
gabrb1Gamma-aminobutyric acid type A receptor subunit beta10.58
hic1HIC ZBTB transcriptional repressor 11.72
hsp90aa1.1Heat shock protein 90, alpha (cytosolic), class A member 1, tandem duplicate 11.82
kif21bKinesin family member 21B2.12
slc27a2bSolute carrier family 27 member 21.72
slc9a5Solute carrier family 9 member A50.56
lancl2Lanthionine synthetase C-like 20.59
Lipid/fatty acid metabolism
mecrMitochondrial trans-2-enoyl-CoA reductase 2.141.97
cyp2k6Cytochrome P450, family 2, subfamily K, polypeptide 61.72
cyp2n13Cytochrome P450, family 2, subfamily N, polypeptide 131.84
hmgcra3-hydroxy-3-methylglutaryl-CoA reductase a1.94
serpinf1Serpin family F member 11.84
Embryogenesis
taf4aTAF4A RNA polymerase II, TATA box binding protein (TBP)-associated factor0.500.52
foxa2Forkhead box A21.74
tacr1aTachykinin receptor 1a0.58
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Wu, C.-C.; Blount, J.R.; Haimbaugh, A.; Heldman, S.; Shields, J.N.; Baker, T.R. Evaluating Phenotypic and Transcriptomic Responses Induced by Low-Level VOCs in Zebrafish: Benzene as an Example. Toxics 2022, 10, 351. https://doi.org/10.3390/toxics10070351

AMA Style

Wu C-C, Blount JR, Haimbaugh A, Heldman S, Shields JN, Baker TR. Evaluating Phenotypic and Transcriptomic Responses Induced by Low-Level VOCs in Zebrafish: Benzene as an Example. Toxics. 2022; 10(7):351. https://doi.org/10.3390/toxics10070351

Chicago/Turabian Style

Wu, Chia-Chen, Jessica R. Blount, Alex Haimbaugh, Samantha Heldman, Jeremiah N. Shields, and Tracie R. Baker. 2022. "Evaluating Phenotypic and Transcriptomic Responses Induced by Low-Level VOCs in Zebrafish: Benzene as an Example" Toxics 10, no. 7: 351. https://doi.org/10.3390/toxics10070351

APA Style

Wu, C. -C., Blount, J. R., Haimbaugh, A., Heldman, S., Shields, J. N., & Baker, T. R. (2022). Evaluating Phenotypic and Transcriptomic Responses Induced by Low-Level VOCs in Zebrafish: Benzene as an Example. Toxics, 10(7), 351. https://doi.org/10.3390/toxics10070351

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop