Identification of potential pathways and microRNA-mRNA networks associated with benzene metabolite hydroquinone-induced hematotoxicity in human leukemia K562 cells
BMC Pharmacology and Toxicology volume 23, Article number: 20 (2022)
Hydroquinone (HQ) is a phenolic metabolite of benzene with a potential risk for hematological disorders and hematotoxicity in humans. In the present study, an integrative analysis of microRNA (miRNA) and mRNA expressions was performed to identify potential pathways and miRNA-mRNA network associated with benzene metabolite hydroquinone-induced hematotoxicity.
K562 cells were treated with 40 μM HQ for 72 h, mRNA and miRNA expression changes were examined using transcriptomic profiles and miRNA microarray, and then bioinformatics analysis was performed.
Out of all the differentially expressed genes (DEGs) and differentially expressed miRNAs (DEMs) induced by HQ, 1482 DEGs and 10 DEMs were up-regulated, and 1594 DEGs and 42 DEMs were down-regulated. HQ-induced DEGs were involved in oxidative stress, apoptosis, DNA methylation, histone acetylation and cellular response to leukemia inhibitory factor GO terms, as well as metabolic, Wnt/β-catenin, NF-κB, and leukemia-related pathways. The regulatory network of mRNAs and miRNAs includes 23 miRNAs, 1108 target genes, and 2304 potential miRNAs-mRNAs pairs. MiR-1246 and miR-224 had the potential to be major regulators in HQ-exposed K562 cells based on the miRNAs-mRNAs network.
This study reinforces the use of in vitro model of HQ exposure and bioinformatic approaches to advance our knowledge on molecular mechanisms of benzene hematotoxicity at the RNA level.
Hydroquinone (HQ) is a phenolic metabolite of benzene representing potential risks for hematological disorders and hematotoxicity [1,2,3]. Extensive research has shown that HQ may contribute to benzene-induced leukemia by oxidative stress, DNA damage, cell cycle regulation, and apoptosis [4,5,6,7,8,9,10]. HQ can inhibit erythroid differentiation in HD3 chicken erythroblast cells, K562 human leukemia cells, and U937 human leukemia cells [4, 11,12,13,14]. Our previous studies suggest that benzene metabolite phenol, 1,2,4-benzenetriol, and HQ considerably inhibit hemin-induced erythroid differentiation in K562 cells [14, 15]. The underlying molecular mechanisms involved in HQ toxicity are not fully understood.
MicroRNAs (miRNAs) are endogenous, highly conserved small non-coding RNAs that modulate various biological processes like cell survival, proliferation, metabolism, differentiation, and apoptosis [16,17,18,19]. The miRNA-mediated gene silencing can involve translational repression, co-translational protein degradation, competition for the cap structure, inhibition of ribosomal subunit joining, inhibition of mRNA circularization through deadenylation, or deadenylation-dependent decay [20, 21]. But recently, miRNAs have been found to mediate gene activation through the bidirectional transcription of the human genome, binding to the enhancer, or recruiting a protein complex with transcriptional activators to the gene promoter [22, 23]. Repression is more common in eukaryotes whereas posttranscriptional upregulation has been observed in specific cell types with distinct transcripts or conditions . MiRNAs have been identified as exerting a powerful effect upon acute myeloid leukemia (AML) development and can be used as potential biomarkers for leukemia diagnosis and prognosis [24,25,26]. Various studies have assessed the role of miRNAs in particulate matter-exposed human pulmonary epithelial cells, sulfur mustard-resistance of the keratinocyte cell line, alternariol and altertoxin II-treated HepG2 cells, and environmental toxicants-exposed aquatic organisms to explore the corresponding toxicity mechanisms [27,28,29,30]. As Liang et al. point out, miRNA-451a and miRNA-486-5p expression are notably lower in HQ-treated CD34+ hematopoietic progenitor cells and K562 cells . HQ may activate apoptotic signals via inhibiting the tumor-suppressive effects of miR-7-5p in TK6 lymphoblastoid cells . However, the role of miRNAs and their target mRNAs expression in benzene hematotoxicity has not been fully addressed yet and needs further study.
Transcriptome sequencing with high sensitivity and good reproducibility is a bargain for detecting low-expression genes . The integration of genomic tools contributes to investigating the molecular mechanism of toxicity . Human leukemia K562 cells were derived from a patient with chronic myeloid leukemia . Here we performed transcriptomic profiles and miRNA microarray to identify mRNAs and miRNAs changes, constructed the mRNAs and miRNAs regulatory network, and performed an in-depth bioinformatic analysis in HQ-induced K562 cells. The differentially expressed mRNAs and miRNAs participated in the metabolic pathways, Wnt/β-catenin pathway, NF-κB pathway, etc. by regulating oxidative stress, apoptosis, DNA methylation, histone acetylation, resulting in hematotoxicity including leukemia. These findings provided a theoretical basis for understanding the molecular mechanisms of benzene hematotoxicity, laying the foundation for future validation of in vivo models as well as therapeutic targets and prognostic factors.
Human leukemia K562 cells were purchased from Cell Resource Center, PekingUnion Medical College (CRC/PUMC, China), and were cultured as described previously [15, 34]. After K562 cells were treated with 40 μM HQ (Sigma-Aldrich) for 72 h, the cells were harvested for further study.
The genome-wide transcriptome analysis was analyzed as described previously . Differentially expressed genes (DEGs) analysis was indicated with the absolute log2 (fold change of HQ/C) values and the adjusted P-value <0.05 was considered differentially expressed. Data are representative of three independent experiments. The volcano plot of DEGs was performed by creating scatter plots in Excel software. Select the data of downregulated DEGs including log2 (fold change of HQ/C) values and the adjusted P-value; choose the scatter plot to represent the relationship between the data sets; Add green color for the dots. The data of upregulated DEGs and genes without differential changes were performed in the same plot and added red and blue colors, respectively. Gene Ontology (GO) enrichment analysis of DEGs was conducted by over-representation analysis (ORA) using the online database WebGestalt (http://www.webgestalt.org). Significantly enriched GO terms in DEGs compared to the genome background were defined by Wallenius’ non-central hypergeometric distribution adjusting for gene length bias . The pathway analysis of DEGs was conducted using the Kyoto Encyclopedia of Genes and Genomes (KEGG) database (http://www.kegg.jp/) .
miRNA microarray analysis
The miRNA microarray was analyzed as described previously . Total RNA was extracted with Ribozol™ RNA Extraction Reagent (Invitrogen, USA) according to the manufacturer’s protocol. The expression of microRNAs contained in the miRBase was analyzed by microarray using the μParaflo™ Microfluidic Biochip Technology of LC sciences (Huston, TX, USA). A total of 1 μg of miRNA-enriched RNA was labeled with Cy3 using the ULS™ microRNA labeling kit (Kreatech, USA) and hybridized on the microarray. Hybridization images were collected using a laser scanner and digitized using Array-Pro image analysis software. The signals were normalized using a locally-weighted regression (LOWESS) filter.
Differentially expressed miRNAs (DEMs) analysis was indicated with the absolute log2 (fold change of HQ/C) values and P-value <0.05 was considered differentially expressed. Data are representative of three independent experiments. Hierarchical clustering analysis of the DEMs was performed and graphs were generated using the ggplot2 package in R software.
miRNA-mRNA regulatory network construction
The potential target mRNAs of miRNAs were screened from miRWalk (http://mirwalk.umm.uni-heidelberg.de), miRanda (https://microbit.org), TargetScan (http://www.targetscan.org), miRDB (http://mirdb.org) and RNA22 (https://cm.jefferson.edu). The miRNA target genes that appeared in at least three prediction databases were screened and quantified from the results of transcriptome analysis. Cytoscape software (version 3.8.0) was employed to construct and analyze the miRNA-mRNA regulatory network.
The effects of HQ on DEGs and DEMs in K562 cells
The volcano plot of overall gene expression showed that 3076 genes were differentially expressed in K562 cells treated with 40 μM HQ for 72 h, including 1482 upregulated and 1594 downregulated DEGs (Fig. 1a). Moreover, there were 84 upregulated and 307 downregulated DEGs over a 2-fold change after HQ exposure. The top 10 upregulated DEGs and top 10 downregulated DEGs based on the absolute log2 (fold change of HQ/C) values after HQ exposure for 72 h were set out in Table 1. What stands out in the table is that LILRA6 (log2 (HQ/C) = 4.17) was the most upregulated DEG and HBB (log2 (HQ/C) = −4.95) was the most downregulated DEG.
The heatmap and hierarchical clustering analysis showed that 52 miRNAs were differentially expressed after HQ exposure, including 10 upregulated and 42 downregulated DEMs (Fig. 1b). In addition, 7 upregulated and 23 downregulated DEMs were over 2-fold change after HQ exposure. These results indicated that majority of genes and miRNAs were downregulated after HQ exposure. The top 10 upregulated DEMs and top 10 downregulated DEMs based on the absolute log2 (fold change of HQ/C) values after HQ exposure for 72 h were set out in Table 2. It is apparent from this table that miR-1246 (log2 (HQ/C) = 4.49) was the most upregulated DEM and miR-4423-3p was the most downregulated DEM.
DEGs involved in the oxidative stress, apoptosis, DNA methylation, histone acetylation and cellular response to leukemia inhibitory factor GO terms in HQ-induced K562 cells
The top 30 enriched GO terms of target DEGs were presented in Fig. 2. GO analysis showed that HQ-upregulated DEGs enriched in biological processes for RNA processing, cellular nitrogen compound metabolic process, and ribonucleoprotein complex biogenesis, whereas HQ-downregulated DEGs enriched in biological processes for single-organism cellular process, localization, and response to stimulus. Our previous studies have demonstrated that DNA methylation, histone acetylation, and ROS have major influences on the transcription of erythroid-specific genes [15, 39,40,41]. Here, we focused on the GO terms of oxidative stress, apoptosis, DNA methylation, and histone acetylation in HQ-induced K562 cells.
HQ treatment considerably increases the relative ROS levels in HD3 chicken erythroblast cells, HL-60 promyelocytic leukemia cells (HL-60 cells), Jurkat T-lymphoblastic leukemia cells (Jurkat cells), and K562 cells [4, 40, 42, 43]. In the present study, there were 38 upregulated DEGs and 37 downregulated DEGs in the response to oxidative stress term (GO:0006979), among which there were 27 upregulated DEGs and 20 downregulated DEGs in the cellular response to oxidative stress term (GO:0034599) in HQ-induced K562 cells. Taken together, these indicated the crucial role of oxidative stress in HQ exposure in K562 cells.
Activation of CASP3 and CASP8 both have pivotal roles in the execution phase of cell apoptosis. As mentioned, HQ induces apoptosis in HL-60 cells, Jurkat cells, human bone marrow mononuclear cells (BMMNC), and K562 cells accompanied by CASP3 or CASP8 activation [7, 44, 45]. Consistently, HQ exposure markedly upregulated the mRNA levels of CASP3 and CASP8 to 1.53-fold and 1.38-fold of that in the control K562 cells in the present study. Moreover, there were 17 upregulated DEGs and 14 downregulated DEGs in the positive regulation of apoptotic signaling pathway (GO:2001235), contributing to HQ-induced apoptosis in K562 cells.
Exposure to benzene metabolites phenol and HQ cause an increase in DNA methylation levels at some erythroid-specific genes clusters, including α-globin, β-globin, and GATA1 gene clusters in K562 cells [39, 46]. There were 5 upregulated DEGs (GATAD2A, PRMT5, BRCA1, BEND3, and MTRR) and 2 downregulated DEGs (TDRKH and SPI1) in the DNA methylation term (GO:0006306) in HQ-induced K562 cells in the present study. To date, we have viewed no data on the effects of HQ on histone acetylation, but the effects of environmental chemicals on histone modifications have been explored recently [47, 48]. Presently, there were 17 upregulated DEGs (LIF, BEND3, ATG5, NAA50, HAT1, TAF5L, POLE3, POLE4, GTF3C4, NOC2L, ING3, BRCA1, KANSL2, WDR5, RUVBL1, ACTL6A, SET, and BEND3) and 5 downregulated DEGs (SPI1, MAPK3, WBP2, FLCN, and YEATS2) in the histone acetylation term (GO:0016573) in HQ-induced K562 cells, indicating their potential roles in HQ exposure.
Exposure to HQ can markedly change the mRNA levels and DNA methylation levels of erythroid-specific genes, as well as reactive oxygen species (ROS) levels in K562 cells [14, 39,40,41]. Transcriptomic analysis contributed to high-throughput data for HQ-induced DEGs in K562 cells, which further support the role of oxidative stress, apoptosis, DNA methylation, histone acetylation in benzene metabolite HQ-induced hematotoxicity. In addition, there were 20 upregulated DEGs and 10 downregulated DEGs in the cellular response to leukemia inhibitory factor term (GO:1990830) in HQ-induced K562 cells (Table 3). Different biological processes in response to HQ exposure might result in these genes expression changes and take part in the mechanism of benzene-induced leukemia.
DEGs involved in the metabolic, Wnt/β-catenin, NF-κB, and leukemia-related pathways in HQ-induced K562 cells
The HQ-upregulated DEGs were annotated with 243 KEGG pathway whereas the downregulated DEGs included 262 KEGG pathway. The top 20 enriched pathways were presented in Fig. 3. KEGG pathway analysis demonstrated that HQ-upregulated DEGs were the most significantly enriched in the ribosome biogenesis in the eukaryotes pathway (annotated 50 DEGs, P = 4.03E-24), spliceosome, RNA transport, and proteasome whereas HQ-downregulated DEGs enriched in the axon guidance pathway (annotated 28 DEGs, P = 4.96E-05), biosynthesis of amino acids, ribosome and carbon metabolism (Fig. 3).
In addition, most genes were annotated as the metabolic pathways, attaining 121genes in both the upregulated and the downregulated DEGs. The cytochrome P450 family 1 subfamily A member 1 (CYP1A1) participates in metabolic pathways. The frizzled class receptor 2 (FZD2) gene encodes a protein in the beta-catenin canonical signaling pathway. HQ has been reported to considerably upregulate the mRNA level of CYP1A1, whereas downregulate the mRNA level of FZD2 in human epidermal keratinocytes cells (HEK cells) . In our results, the mRNA levels of CYP1A1 and FZD2 in K562 cells showed similar trends to the HEK cells.
The canonical Wnt/β-catenin pathway is a key factor in multiple biological processes like cell proliferation, apoptosis, motility, and differentiation. HQ treatment reduces the level of β-catenin, a dominant part in the canonical Wnt signaling pathway in BMMNC cells . Presently, HQ considerably downregulated the mRNA levels of WNT6 and WNT11, as well as other 12 genes (FZD2, PLCB2, VANGL2, NFATC4, FZD4, FRAT1, PPARD, CTNNBIP1, RAC2, LRP5, FRAT2, and NFATC2) in the Wnt signaling pathway (hsa04310) (Supplementary Fig. 1).
RNA-seq has identified that hydroquinone promotes DNA homologous recombination repair via activating the NF-κB pathway in the human osteosarcoma cell line (U2OS/DR-GFP) . In our study, the mRNA level of NFKB1 was upregulated after HQ exposure, indicating the NF-κB pathway might be activated in HQ-induced K562 cells.
Furthermore, there were 10, 12, 13, and 39 DEGs annotated as the KEGG terms of chronic myeloid leukemia, acute myeloid leukemia, hematopoietic cell lineage, and human T-cell leukemia virus 1 infection in HQ-induced K562 cells, respectively (Table 4) . The related DEGs were highlighted in the KEGG maps in Supplementary Fig. 2–5. HQ exposure might change these genes’ expression to activate the leukemia-related pathways, resulting in potential benzene hematotoxicity.
Identification of miR-1246 and miR-224 as major regulators in the miRNA-mRNA network in HQ-induced K562 cells
The mutual regulation of DEMs and DEGs was analyzed to construct an interaction network after HQ exposure in K562 cells. It is a widely held view that miRNAs can induce the degradation of their target mRNAs, so we focused on those miRNAs and mRNAs with opposite expression trends. There were 10 upregulated DEMs and 1594 downregulated DEGs, as well as 42 downregulated DEMs and 1482 upregulated DEGs in HQ-induced K562 cells. The entire analysis led to the identification of 23 miRNAs, 1108 target genes, and 2304 potential miRNAs-mRNAs pairs in the network. Moreover, the miRNA-mRNA network of DEGs and DEMs over a 2-fold change was constructed by a cystoscope (Fig. 4).
The miRNA-mRNA network was constructed with 109 potential pairs containing 7 upregulated DEMs and 75 downregulated target DEGs in Fig. 4a; the network had 30 potential pairs containing 8 upregulated DEMs and 15 downregulated target DEGs in Fig. 4b. Per miRNA had 2–29 target genes and per gene was connected to 1–5 miRNAs. The lists of top miRNAs and target genes in the network were shown in Table 5 and Table 6. On the other hand, we found that majority of DEGs and DEMs were downregulated after HQ exposure.
HMDD (the Human microRNA Disease Database) is a database for experiment-based miRNA and human disease associations (http://www.cuilab.cn/hmdd) . When the DEMs over a 2-fold change were screened from the HMDD 3.2, 5 miRNAs (miR-27a, miR-224, miR-1246, miR-18a, and miR-18b) were associated with leukemia, indicating potential roles in benzene-induced leukemia. According to the miRNA-mRNA regulatory network, there were more interactions between upregulated miRNAs and downregulated genes than the other way around. In our study, miR-1246 and miR-224 were the most upregulated DEMs with 17 target DEGs whereas miR-18a and miR-18b were downregulated in HQ-induced K562 cells. However, it remains unknown why MiR-27a was upregulated in HQ-induced K562 cells but downregulated in acute leukemia cell lines and primary samples compared to hematopoietic stem-progenitor cells . Therefore, miR-1246 and miR-224 were identified to be major regulators for HQ exposure in K562 cells. These results indicated the crucial role of miR-1246 and miR-224 in benzene hematotoxicity.
Our previous studies have demonstrated that hemin-induced erythroid differentiation is concentration-dependently and time-dependently inhibited by benzene metabolites exposure (phenol, 1,2,4-benzenetriol, and hydroquinone) in K562 cells [14, 15, 39, 46]. In the present study, the concentration of 40 μM HQ was selected to correspond to no obvious cytotoxicity but markedly inhibiting erythroid differentiation in K562 cells when exposed for 72 h . Hemoglobin subunit beta (HBB) loci are associated with beta-thalassemia, sickle cell anemia, and heinz body anemias [53,54,55] . HBB was the most dramatically downregulated gene after HQ exposure, which might play important roles in benzene hematotoxicity.
It has been reported that the NF-κB pathway is activated in the development of chronic myeloid leukemia and acute myeloid leukemia [56,57,58,59]. In our study, HQ exposure upregulated the mRNA level of NFKB1 in chronic myeloid leukemia and acute myeloid leukemia pathway, which partly supported that HQ might promote leukemia development by activating the NF-κB pathway. HQ exposure downregulated the mRNA levels of MAP 2 K2 and MAPK3 in chronic myeloid leukemia and acute myeloid leukemia pathway (Table 4). This may inhibit the MAPK signaling pathway and thus inhibit cell proliferation (Supplementary Fig. 2, 3). These results were consistent with our previous study that HQ induced a concentration-dependent decrease in the viabilities in K562 cells . The roles of metabolic pathways, Wnt/β-catenin pathway, and NF-κB pathway in benzene hematotoxicity need further study.
The downregulated miR-1246/1248 are key nodes that reveal the possible relapse mechanisms for pediatric T cell acute lymphoblastic leukemia. As reported, miR-1246 is one of the most highly enriched miRNAs in AML derived extracellular vesicles [16, 60, 61]. MiR-1246, a hundreds-fold alteration in microvesicles from three different leukemia cell lines (K562, Nalm-6, and Jurkat), is upregulated to activate the expression of C6orf211 and C19orf10 to promote tumor progression in patients diagnosed as chronic myeloid leukemia . MiR-224 expression is considerably upregulated in the bone marrow of pediatric AML patients and can be used as noninvasive biomarkers for the early prediction of hepatocellular carcinoma development . These results indicated the crucial role of miR-1246 and miR-224 in hematotoxicity. MiR-29a is highly expressed in arsenic-induced peripheral neuropathy, which is consistent with higher expression after HQ exposure in this study . Low expression of miR-18a distinguishes pediatric and adult acute lymphoblastic leukemia from each other . PML/RARα-regulated miR-181a/b cluster targets the tumor suppressor RASSF1A in acute promyelocytic leukemia .
In further research, the expression, target genes, and biological function of miR-1246 and miR-224 will be confirmed by more techniques in HQ-induced K562 cells, CD34+ hematopoietic progenitor cells, U937 human leukemia cells, and human peripheral blood mononuclear cells (PBMCs). Inhibition of miR-1246 and miR-224 will contribute to their regulation for HQ exposure. On the other hand, hydroquinone has been reported to inhibit PRV infection in mouse neuroblastoma N2a cells, protect neurons from transient cerebral ischemia, and reduce gliosis in a gerbil model of transient cerebral ischemia [67, 68]. The study on the effect of HQ on neural systems will contribute to a full understanding of benzene toxicity. Further assessment in suitable animal models or any data of miRNAs expression of benzene-exposed patients will be immensely beneficial to further studies.
In summary, the present study identified differentially expressed genes and miRNAs in HQ-induced K562 cells using transcriptomic profiles and miRNA microarray. The miRNA-mRNA network can help us better understand the molecular mechanisms between miRNAs and their target genes. MiR-1246 and miR-224 had the potential to be major regulators for HQ exposure in K562 cells based on the miRNAs-mRNAs network and were reported to be associated with leukemia, suggesting potential biomarkers for the evaluation of benzene hematotoxicity. Apart from the miRNAs-mRNAs regulation, how miRNAs regulate protein expression remains elusive. Further proteomics study is to be performed to elucidate the underlying mechanism. The GO and KEGG pathways provide a framework for further studies in suitable in vitro and animal models, which will contribute to developing new strategies for the prevention and rapid diagnosis of benzene hematotoxicity.
Availability of data and materials
The datasets used and/or analyzed during the current study are available from the corresponding author on reasonable request.
Aksoy M. Hematotoxicity and carcinogenicity of benzene. Environ Health Perspect. 1989;82:193–7.
Huff J. Benzene-induced cancers: Abridged history and occupational health impact. Int J Occup Environ Health. 2007;13:213–21.
Rothman N, Smith MT, Hayes RB, Li GL, Irons RD, Dosemeci M, et al. An epidemiologic study of early biologic effects of benzene in Chinese workers. Environ Health Perspect. 1996;104:1365–70.
Badham HJ, Winn LM. In utero and in vitro effects of benzene and its metabolites on erythroid differentiation and the role of reactive oxygen species. Toxicol Appl Pharmacol. 2010;224:273–9.
Gaskell M, McLuckie KI, Farmer PB. Genotoxicity of the benzene metabolites para-benzoquinone and hydroquinone. Chem Biol Interact. 2005;153:267–70.
Hatzi VI, Terzoudi GI, Pantelias GE, Spiliopoulou C, Makropoulos V. The benzene metabolite hydroquinone enhances G2-chromosomal radiosensitivity by inducing a less-efficient G2-M-checkpoint in irradiated lymphocytes. Int J Oncol. 2007;31:145–52.
Inayat-Hussain SH, Winski SL, Ross D. Differential involvement of caspases in hydroquinone-induced apoptosis in human leukemic hl-60 and jurkat cells. Toxicol Appl Pharmacol. 2001;175:95–103.
McDonald TA, Holland NT, Skibola C, Duramad P, Smith MT. Hypothesis: Phenol and hydroquinone derived mainly from diet and gastrointestinal flora activity are causal factors in leukemia. Leukemia. 2001;15:10–20.
Smith MT. The mechanism of benzene-induced leukemia: A hypothesis and speculations on the causes of leukemia. Environ Health Perspect. 1996;104:1219–25.
Tan Q, Li J, Peng J, Liu Z, Liu J, Zhang H, et al. E4F1 silencing inhibits the cell growth through cell-cycle arrest in malignant transformed cells induced by droquinone. J Biochem Mol Toxicol. 2019;33:e22269.
Chen YJ, Liu WH, Chang LS. Hydroquinone-induced FOXP3-ADAM17-Lyn-Akt-p21 signaling axis promotes malignant progression of human leukemia U937 cells. Arch Toxicol. 2017;91:983–97.
Liang B, Chen Y, Yuan W, Qin F, Zhang Q, Deng N, et al. Down-regulation of miRNA-451a and miRNA-486-5p involved in benzene-induced inhibition on erythroid cell differentiation in vitro and in vivo. Arch Toxicol. 2018;92:259–72.
Ma X, Zhang X, Luo J, Liang B, Peng J, Chen C, et al. MiR-486-5p-directed MAGI1/Rap1/RASSF5 signaling pathway contributes to hydroquinone-induced inhibition of erythroid differentiation in K562 cells. Toxicol in Vitro. 2020;66:104830.
Wu XR, Xue M, Li XF, Wang Y, Wang J, Han QL, et al. Phenolic metabolites of benzene inhibited the erythroid differentiation of K562 cells. Toxicol Lett. 2011;203:190–9.
Yu CH, Li Y, Zhao X, Yang SQ, Li L, Cui NX, et al. Benzene metabolite 1,2,4-benzenetriol changes DNA methylation and histone acetylation of erythroid-specific genes in K562 cells. Arch Toxicol. 2019;93:137–47.
Butler JT, Abdelhamed S, Gao L, Lim J, Horton TM, Kurre P. Leukemic stress targets the mTOR pathway to suppress residual HSC in the BM microenvironment. Blood. 2019;134:3730.
Guo H, Ingolia NT, Weissman JS, Bartel DP. Mammalian microRNAs predominantly act to decrease target mRNA levels. Nature. 2010;466:835–40.
Siomi H, Siomi MC. Posttranscriptional regulation of microRNA biogenesis in animals. Mol Cell. 2010;38:323–32.
Valencia-Sanchez MA, Liu J, Hannon GJ, Parker R. Control of translation and mRNA degradation by miRNAs and siRNAs. Genes Dev. 2006;20:515–24.
Eulalio A, Huntzinger E, Izaurralde E. Getting to the root of miRNA-mediated gene silencing. Cell. 2008;132(1):9–14.
Behm-Ansmant I, Rehwinkel J, Izaurralde E. MicroRNAs silence gene expression by repressing protein expression and/or by promoting mRNA decay. Cold Spring Harb Symp Quant Biol. 2006;71:523–30.
Pu M, Chen J, Tao Z, Miao L, Qi X, Wang Y, et al. Regulatory network of miRNA on its target: coordination between transcriptional and post-transcriptional regulation of gene expression. Cell Mol Life Sci. 2019;76(3):441–51.
Vasudevan S. Posttranscriptional upregulation by microRNAs. Wiley Interdiscip Rev RNA. 2012;3:311–30.
Marcucci G, Maharry KS, Metzeler KH, Volinia S, Wu YZ, Mrózek K, et al. Clinical role of microRNAs in cytogenetically normal acute myeloid leukemia: miR-155 upregulation independently identifies high-risk patients. J Clin Oncol. 2013;31:2086–93.
Lopotova T, Zackova M, Klamova H, Moravcova J. MicroRNA-451 in chronic myeloid leukemia: miR-451-BCR-ABL regulatory loop? Leuk Res. 2011;35:974–7.
Wallace JA, O'Connell RM. MicroRNAs and acute myeloid leukemia: therapeutic implications and emerging concepts. Blood. 2017;130:1290–301.
Bertucci A, Pierron F, Ye T, Gonzalez P, Couture P, Baudrimont M. Identification and expression of microRNAs in european eels Anguilla anguilla from two natural sites with different pollution levels. Environ Pollut. 2019;250:274–83.
Rothmiller S, Wolf M, Worek F, Steinritz D, Thiermann H, Schmidt A. Alteration of miRNA expression in a sulfur mustard resistant cell line. Toxicol Lett. 2018;293:38–44.
Vejdovszky K, Sack M, Jarolim K, Aichinger G, Somoza MM, Marko D. In vitro combinatory effects of the Alternaria mycotoxins alternariol and altertoxin II and potentially involved miRNAs. Toxicol Lett. 2017;267:45–52.
Xu Z, Wang N, Xu Y, Hua L, Zhou D, Zheng M, et al. Effects of chronic PM2.5 exposure on pulmonary epithelia: Transcriptome analysis of mRNA-exosomal miRNA interactions. Toxicol Lett. 2019;316:49–59.
Luo H, Liang H, Chen Y, Chen S, Xu Y, Xu L, et al. miR-7-5p overexpression suppresses cell proliferation and promotes apoptosis through inhibiting the ability of DNA damage repair of PARP-1 and BRCA1 in TK6 cells exposed to hydroquinone. Chem Biol Interact. 2018;283:84–90.
Yadav NK, Shukla P, Omer A, Pareek S, Singh RK. Next generation sequencing: potential and application in drug discovery. Sci World J. 2014;2014:802437.
Robbens J, van der Ven K, Maras M, Blust R, De Coen W. Ecotoxicological risk assessment using DNA chips and cellular reporters. Trends Biotechnol. 2007;25:460–6.
Lozzio CB, Lozzio BB. Human chronic myelogenous leukemia cell-line with positive Philadelphia chromosome. Blood. 1975;45:321–34.
Gao F, Kataoka M, Liu N, Liang T, Huang ZP, Gu F, et al. Therapeutic role of miR-19a/19b in cardiac regeneration and protection from myocardial infarction. Nat Commun. 2019;10:1802–15.
Young MD, Wakefield MJ, Smyth GK, Oshlack A. Gene ontology analysis for RNA-seq: accounting for selection bias. Genome Biol. 2010;11:R14.
Kanehisa M, Goto S. KEGG: Kyoto encyclopedia of genes and genomes. Nucleic Acids Res. 2000;28:27–30.
Guo J, Miao Y, Xiao B, Huan R, Jiang Z, Meng D, et al. Differential expression of microRNA species in human gastric cancer versus non-tumorous tissues. J Gastroenterol Hepatol. 2009;24:652–7.
Li Y, Wu XR, Li XF, Suriguga YCH, Li YR, Yi ZC. Changes in DNA methylation of erythroid-specific genes in K562 cells exposed to phenol and hydroquinone. Toxicology. 2013;312:108–14.
Yu CH, Suriguga LY, Li YR, Tang KY, Jiang L, Yi ZC. The role of ROS in hydroquinone-induced inhibition of K562 cell erythroid differentiation. Biomed Environ Sci. 2014;27:212–4.
Yu CH, Jiang L, Wang Y, Cui NX, Zhao X, Yi ZC. Inhibition of erythroid differentiation of human leukemia K562 cells by N-acetylcysteine and ascorbic acid through downregulation of ROS. Biomed Environ Sci. 2018;31:247–51.
Ishihama M, Toyooka T, Ibuki Y. Toxicology in Vitro Generation of phosphorylated histone H2AX by benzene metabolites. Toxicol in Vitro. 2008;22:1861–8.
Kim YJ, Woo HD, Kim BM, Lee YJ, Kang SJ, Cho YH, et al. Risk assessment of hydroquinone: differential responses of cell growth and lethality correlated to hydroquinone concentration. J Toxicol Environ Health A. 2009;72:1272–8.
Wang Y, Zhang GY, Han QL, Wang J, Suriguga LY, Yu CH, et al. Phenolic metabolites of benzene induced caspase-dependent cytotoxicities to K562 cells accompanied with decrease in cell surface sialic acids. Environ Toxicol. 2014;29:1437–51.
Li J, Jiang S, Chen Y, Ma R, Chen J, Qian S, et al. Benzene metabolite hydroquinone induces apoptosis of bone marrow mononuclear cells through inhibition of β-catenin signaling. Toxicol in Vitro. 2018;46:361–9.
Tang KY, Yu CH, Jiang L, Gong M, Liu WJ, Wang Y, et al. Long-term exposure of K562 cells to benzene metabolites inhibited erythroid differentiation and elevated methylation in erythroid specific genes. Toxicol Res. 2016;5:1284–97.
Baccarelli A, Bollati V. Epigenetics and environmental chemicals. Curr Opin Pediatr. 2009;21:243–51.
Philbrook NA, Winn LM. Investigating the effects of in utero benzene exposure on epigenetic modifications in maternal and fetal CD-1 mice. Toxicol Appl Pharmacol. 2015;289:12–9.
Pyo JJ, Ahn S, Jin SH, An S, Lee E, Choi J, et al. Keratinocyte-derived IL-36γ plays a role in hydroquinone-induced chemical leukoderma through inhibition of melanogenesis in human epidermal melanocytes. Arch Toxicol. 2019;93:2307–20.
Yang X, Lu Y, He F, Hou F, Xing C, Xu P, et al. Benzene metabolite hydroquinone promotes DNA homologous recombination repair via the NF-κB pathway. Carcinogenesis. 2019;40:1021–30.
Huang Z, Shi J, Gao Y, Cui C, Zhang S, Li J, et al. HMDD v3.0: a database for experimentally supported human microRNA-disease associations. Nucleic Acids Res. 2019;47:D1013–7.
Scheibner KA, Teaboldt B, Hauer MC, Chen X, Cherukuri S, Guo Y, et al. MiR-27a functions as a tumor suppressor in acute leukemia by regulating 14-3-3θ. PLoS One. 2012;7:e50895.
Akinsheye I, Alsultan A, Solovieff N, Ngo D, Baldwin CT, Sebastiani P, et al. Fetal hemoglobin in sickle cell anemia. Blood. 2011;118:19–27.
Galanello R, Origa R. Beta-thalassemia. Orphanet J Rare Dis. 2010;5:11.
Poyart C, Wajcman H. Hemolytic anemias due to hemoglobinopathies. Mol Asp Med. 1996;17:129–42.
Wang J, Li Q, Wang C, Xiong Q, Lin Y, Sun Q, et al. Knock-down of CIAPIN1 sensitizes K562 chronic myeloid leukemia cells to Imatinib by regulation of cell cycle and apoptosis-associated members via NF-κB and ERK5 signaling pathway. Biochem Pharmacol. 2016;99:132–45.
Chorzalska A, Ahsan N, Rao RSP, Roder K, Yu X, Morgan J, et al. Overexpression of Tpl2 is linked to imatinib resistance and activation of MEK-ERK and NF-κB pathways in a model of chronic myeloid leukemia. Mol Oncol. 2018;12(5):630–47.
Shen MY, Wang Y, Cui SY, Wu XL, Guo Y, Xu RR. MicroRNA-125a regulates proliferation and apoptosis of acute myeloid leukemia through targeting NF-κB pathway. Eur Rev Med Pharmacol Sci. 2019;23(9):3594–601.
Fisher DAC, Malkova O, Engle EK, Miner CA, Fulbright MC, Behbehani GK, et al. Mass cytometry analysis reveals hyperactive NF Kappa B signaling in myelofibrosis and secondary acute myeloid leukemia. Leukemia. 2017;31(9):1962–74.
Abdelhamed S, Butler JT, Doron B, Halse A, Nemecek E, Wilmarth PA, et al. Extracellular vesicles impose quiescence on residual hematopoietic stem cells in the leukemic niche. EMBO Rep. 2019;20:e47546.
Luo M, Zhang Q, Xia M, Hu F, Ma Z, Chen Z, et al. Differential co-expression and regulatory network analysis uncover the relapse factor and mechanism of T cell acute leukemia. Mol Ther Nucleic Acids. 2018;12:184–94.
Lu L, Li H, Chen X, Xiong W, Huang S. The chromosome open reading frame genes targeted by abnormal micrornas in microvesicles from chronic myeloid leukemia. Blood. 2014;124:5509.
Amr KS, Elmawgoud Atia HA, Elazeem Elbnhawy RA, Ezzat WM. Early diagnostic evaluation of miR-122 and miR-224 as biomarkers for hepatocellular carcinoma. Genes Dis. 2017;4:215–21.
Chatterjee D, Bandyopadhyay A, Sarma N, Basu S, Roychowdhury T, Roy SS, et al. Role of microRNAs in senescence and its contribution to peripheral neuropathy in the arsenic exposed population of West Bengal, India. Environ Pollut. 2018;233:596–603.
Mosakhani N, Missiry ME, Vakkila E, Knuutila S, Vakkila J. Low expression of miR-18a as a characteristic of pediatric acute lymphoblastic leukemia. J Pediatr Hematol Oncol. 2017;39:585–8.
Bräuer-Hartmann D, Hartmann JU, Wurm AA, Gerloff D, Katzerke C, Verga Falzacappa MV, et al. PML/RARα-Regulated miR-181a/b Cluster Targets the Tumor Suppressor RASSF1A in Acute Promyelocytic Leukemia. Cancer Res. 2015;75:3411–24.
Fang L, Gao Y, Lan M, Jiang P, Bai J, Li Y, et al. Hydroquinone inhibits PRV infection in neurons in vitro and in vivo. Vet Microbiol. 2020;250:108864.
Park JH, Park CW, Ahn JH, Choi SY, Shin MC, Cho JH, et al. Neuroprotection and reduced gliosis by pre- and post-treatments of hydroquinone in a gerbil model of transient cerebral ischemia. Chem Biol Interact. 2017;278:230–8.
This work was supported by grants from the National Natural Science Foundation of China (Project No. 81070091 and 81573192).
Ethics approval and consent to participate
Consent for publication
The authors declare that there is no conflict of interest.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
. HQ-regulated DEGs of K562 cells in the Wnt signaling pathway. The Wnt signaling pathway was downloaded from the KEGG database (https://www.kegg.jp/kegg-bin/show_pathway?map04310). In K562 cells, HQ-upregulated DEGs were colored pink and HQ-downregulated DEGs were colored green. Supplementary Figure 2. HQ-regulated DEGs of K562 cells in chronic myeloid leukemia pathway. The chronic myeloid leukemia pathway was downloaded from the KEGG database (https://www.kegg.jp/kegg-bin/show_pathway?map05220). In K562 cells, HQ-upregulated DEGs were colored pink and HQ-downregulated DEGs were colored green. Supplementary Figure 3. HQ-regulated DEGs of K562 cells in acute myeloid leukemia pathway. The acute myeloid leukemia pathway was downloaded from the KEGG database (https://www.kegg.jp/kegg-bin/show_pathway?map05221). In K562 cells, HQ-upregulated DEGs were colored pink and HQ-downregulated DEGs were colored green. Supplementary Figure 4. HQ-regulated DEGs of K562 cells in hematopoietic cell lineage pathway. The hematopoietic cell lineage pathway was downloaded from the KEGG database (https://www.kegg.jp/kegg-bin/show_pathway?map04640). In K562 cells, HQ-upregulated DEGs were colored pink and HQ-downregulated DEGs were colored green. Supplementary Figure 5. HQ-regulated DEGs of K562 cells in human T-cell leukemia virus 1 infection pathway. The human T-cell leukemia virus 1 infection pathway was downloaded from the KEGG database (https://www.kegg.jp/kegg-bin/show_pathway?map05166). In K562 cells, HQ-upregulated DEGs were colored pink and HQ-downregulated DEGs were colored green.
About this article
Cite this article
Yu, CH., Yang, SQ., Li, L. et al. Identification of potential pathways and microRNA-mRNA networks associated with benzene metabolite hydroquinone-induced hematotoxicity in human leukemia K562 cells. BMC Pharmacol Toxicol 23, 20 (2022). https://doi.org/10.1186/s40360-022-00556-8