MSTRG3207 promotes apoptosis in zebrafish ZF4 cells via sponging dre-miR-736/bbc3/LOC101885512 axis during cold acclimation.

Long non-coding RNAs (lncRNAs) play essential roles in a variety of biological processes. It has been recently reported that lncRNAs can regulate mRNA expression by binding to microRNAs (miRNAs) as competing endogenous RNAs (ceRNAs). However, the involvement of this regulatory mechanism during cold acclimation in fish remains unclear. In this study, we constructed a ceRNA network mediated by lncRNAs in cold-acclimated zebrafish ZF4 cells through bioinformatic analysis of the mRNA, miRNA, and lncRNA profiles obtained from ZF4 cells cultured at 18 °C for 30 days. A previously uncharacterized lncRNA, MSTRG3207, was selected for further analysis. MSTRG3207 was upregulated and dre-miR-736 was downregulated during cold acclimation. MSTRG3207 was cloned by rapid amplification of cDNA ends (RACE) and functionally characterized. The binding of MSTRG3207 to dre-miR-736 was validated by dual-luciferase reporter assay. Under cold acclimation, MSTRG3207 promoted apoptosis by sponging dre-miR-736 and upregulating bbc3 and LOC101885512, two apoptotic genes targeted by dre-miR-736. Taken together, our findings indicate that MSTRG3207 upregulation promotes apoptosis by sponging dre-miR-736 during cold acclimation in fish.

MiR-202-3p determines embryo viability during mid-blastula transition.

Developmental growth is an intricate process involving the coordinated regulation of the expression of various genes, and microRNAs (miRNAs) play crucial roles in diverse processes throughout animal development. The mid-blastula transition (MBT) is a developmental milestone when maternal RNAs are cleared and the zygotic genome programmed asynchronous cell division begins to drive embryogenesis. While mechanisms underlying MBT have been intensively revealed, factors regulating cell proliferation at the transition remain largely unknown. We report here a microRNA, miR-202-3p to be a key factor that determines embryonic fate during MBT in zebrafish. A miR-202-3p antagomir specifically terminated embryo development at the mid-blastula stage. In vivo deletion of the miR-202 locus recapitulated the fatal phenotypes, which were rescued only by miR-202-3p or its precursor. Transcriptome comparison revealed >250 RNAs including both maternal and zygotic origins were dysregulated at MBT in the miR-202-/- embryos, corresponding with arrays of homeostatic disorders leading to massive apoptosis. A trio of genes: nfkbiaa, perp and mgll, known to be intimately involved with cell proliferation and survival, were identified as direct targets of miR-202-3p. Importantly, over- or under-expression of any of the trio led to developmental delay or termination at the blastula or gastrula stages. Furthermore, nfkbiaa and perp were shown to inter-regulate each other. Thus, miR-202-3p mediates a regulatory network whose components interact closely during MBT to determine embryonic viability and development.

Regulation of MYB by distal enhancer elements in human myeloid leukemia.

MYB plays vital roles in regulating proliferation and differentiation of hematopoietic progenitor cells, dysregulation of MYB has been implicated in the pathogenesis of leukemia. Although the transcription of MYB has been well studied, its detailed underlying regulatory mechanisms still remain elusive. Here, we detected the long-range interaction between the upstream regions, -34k and -88k, and the MYB promoter in K562, U937, and HL-60 cells using circularized chromosome conformation capture (4C) assay, which declined when MYB was downregulated during chemical-induced differentiation. The enrichment of enhancer markers, H3K4me1 and H3K27ac, and enhancer activity at the -34k and -88k regions were confirmed by ChIP-qPCR and luciferase assay respectively. ChIP-qPCR showed the dynamic binding of GATA1, TAL1, and CCAAT/enhancer-binding protein (C/EBPß) at -34k and -88k during differentiation of K562 cells. Epigenome editing by a CRISPR-Cas9-based method showed that H3K27ac at -34k enhanced TF binding and MYB expression, while DNA methylation inhibited MYB expression. Taken together, our data revealed that enhancer elements at -34k are required for MYB expression, TF binding, and epigenetic modification are closely involved in this process in human myeloid leukemia cells.

Risk SNP-induced lncRNA-SLCC1 drives colorectal cancer through activating glycolysis signaling.

Long non-coding RNAs (lncRNAs) play key roles in colorectal carcinogenesis. Here, we aimed to identify the risk SNP-induced lncRNAs and to investigate their roles in colorectal carcinogenesis. First, we identified rs6695584 as the causative SNP in 1q41 locus. The A>G mutation of rs6695584 created a protein-binding motif of BATF, altered the enhancer activity, and subsequently activated lncSLCC1 expression. Further validation in two independent CRC cohorts confirmed the upregulation of lncSLCC1 in CRC tissues, and revealed that increased lncSLCC1 expression was associated with poor survival in CRC patients. Mechanistically, lncRNA-SLCC1 interacted with AHR and transcriptionally activated HK2 expression, the crucial enzyme in glucose metabolism, thereby driving the glycolysis pathway and accelerating CRC tumor growth. The functional assays revealed that lncSLCC1 induced glycolysis activation and tumor growth in CRC mediated by HK2. In addition, HK2 was upregulated in colorectal cancer tissues and positively correlated with lncSLCC1 expression and patient survival. Taken together, our findings reveal a risk SNP-mediated oncogene lncRNA-SLCC1 promotes CRC through activating the glycolysis pathway.

Inhibition of Embryonic HSP 90 Function Promotes Variation of Cold Tolerance in Zebrafish.

Accumulating evidence indicates that heat shock protein 90 (HSP90) plays essential roles in modulation of phenotypic plasticity in vertebrate development, however, the roles of HSP90 in modulation of cold tolerance capacity in fish are still unclear. In the present study, we showed that transient inhibition of embryonic HSP90 function by a chemical inhibitor or low conductivity stress promoted variation of cold tolerance capacity in adult zebrafish. Further work showed that embryonic HSP90 inhibition enhanced cold tolerance in adult zebrafish could be transmitted to their offspring. RNA-seq data showed that embryonic HSP90 inhibition enhanced cold tolerance involves variation of gene expression related to proteasome, lysosome, autophagy, and ribosome. Experiments with zebrafish ZF4 cells showed that two differentially expressed genes atg9b and psmd12 were up-regulated by radicicol treatment and provided protective roles for cells under cold stress, indicating that up-regulation of autophagy and proteasome function contributes to enhanced cold tolerance. The present work sheds a light on the roles of HSP90 in regulation of phenotypic plasticity associated with thermal adaptation in fish.

Identification and characterization of miRNAs involved in cold acclimation of zebrafish ZF4 cells.

MicroRNAs (miRNAs) play vital roles in various biological processes under multiple stress conditions by leading to mRNA cleavage or translational repression. However, the detailed roles of miRNAs in cold acclimation in fish are still unclear. In the present study, high-throughput sequencing was performed to identify miRNAs from 6 small RNA libraries from the zebrafish embryonic fibroblast ZF4 cells under control (28°C, 30 days) and cold-acclimation (18°C, 30 days) conditions. A total of 414 miRNAs, 349 known and 65 novel, were identified. Among those miRNAs, 24 (19 known and 5 novel) were up-regulated, and 23 (9 known and 14 novel) were down-regulated in cold acclimated cells. The Gene Ontology (GO) and Kyoto Encyclopaedia of Genes and Genomes (KEGG) enrichment analyses indicated that the target genes of known differentially expressed miRNAs (DE-miRNA) are involved in cold acclimation by regulation of phosphorylation, cell junction, intracellular signal transduction, ECM-receptor interaction and so on. Moreover, both miR-100-3p inhibitor and miR-16b mimics could protect ZF4 cells under cold stress, indicating the involvement of miRNA in cold acclimation. Further study showed that miR-100-3p and miR-16b could regulate inversely the expression of their target gene (atad5a, cyp2ae1, lamp1, rilp, atxn7, tnika, btbd9), and that overexpression of miR-100-3p disturbed the early embryonic development of zebrafish. In summary, the present data show that miRNAs are closely involved in cold acclimation in zebrafish ZF4 cells and provide information for further understanding of the roles of miRNAs in cold acclimation in fish.

Characterization of lncRNAs involved in cold acclimation of zebrafish ZF4 cells.

Long non-coding RNAs (lncRNAs) are increasingly regarded as a key role in regulating diverse biological processes in various tissues and species. Although the cold responsive lncRNAs have been reported in plants, no data is available on screening and functional prediction of lncRNAs in cold acclimation in fish so far. Here we compared the expression profile of lncRNAs in cold acclimated zebrafish embryonic fibroblast cells (ZF4) cultured at 18°C for 30 days with that of cells cultured at 28°C as control by high-throughput sequencing. Totally 8,363 novel lncRNAs were identified. Including known and novel lncRNAs, there are 347 lncRNAs up-regulated and 342 lncRNAs down-regulated in cold acclimated cells. Among the differentially expressed lncRNAs, 74 and 61 were detected only in control cells or cold-acclimated cells, respectively. The Gene Ontology (GO) and Kyoto Encyclopaedia of Genes and Genomes (KEGG) enrichment analyses of adjacent genes to the differentially expressed lncRNAs showed that the enriched genes are involved in electron transport, cell adhesion, oxidation-reduction process, and so on. We also predicted the target genes of the differentially expressed lncRNAs by looking for interactions between lncRNAs and mRNAs, and constructed an interaction network. In summary, our genome-wide systematic identification and functional prediction of cold responsive lncRNAs in zebrafish cells suggests a crucial role of lincRNAs in cold acclimation in fish.

Characterization of Selenoprotein P cDNA of the Antarctic toothfish Dissostichus mawsoni and its role under cold pressure.

Our previous study using comparative genome analysis revealed a significant gene copy number gain of Dissostichus mawsoni selenoprotein P (Dm-SEPP) during the evolutionary radiation of Antarctic notothenioids, suggesting that Dm-SEPP contribute to this process, but the detailed structure and function of this gene product remain unclear. In the present study, the Dm-SEPP cDNA was cloned and characterized. The Dm-SEPP cDNA contains 17 selenocysteines (Sec) encoded by TGA codons and 2 typical SECIS elements located in the 3'-UTR. Evolutionary analysis of the Dm-SEPP gene revealed that it's closely related to the SEPP gene of zebrafish (Danio rerio), showing 51% amino acid similarity. Over-expression of Dm-SEPP could protect mammalian cells under cold pressure, probably via eliminating ROS. Further study showed an increase of endogenous SEPP in zebrafish ZF4 cells under cold pressure, and knockdown of SEPP decreased cell viability, accompanied with increased ROS. Our results suggested a protective role of Dm-SEPP in cold adaptation in Antarctic notothenioids.

Distal regulation of c-myb expression during IL-6-induced differentiation in murine myeloid progenitor M1 cells.

The c-Myb transcription factor is a major regulator that controls differentiation and proliferation of hematopoietic progenitor cells, which is frequently deregulated in hematological diseases, such as lymphoma and leukemia. Understanding of the mechanisms regulating the transcription of c-myb gene is challenging as it lacks a typical promoter and multiple factors are involved. Our previous studies identified some distal regulatory elements in the upstream regions of c-myb gene in murine myeloid progenitor M1 cells, but the detailed mechanisms still remain unclear. In the present study, we found that a cell differentiation-related DNase1 hypersensitive site is located at a -28k region upstream of c-myb gene and that transcription factors Hoxa9, Meis1 and PU.1 bind to the -28k region. Circular chromosome conformation capture (4C) assay confirmed the interaction between the -28k region and the c-myb promoter, which is supported by the enrichment of CTCF and Cohesin. Our analysis also points to a critical role for Hoxa9 and PU.1 in distal regulation of c-myb expression in murine myeloid cells and cell differentiation. Overexpression of Hoxa9 disrupted the IL-6-induced differentiation of M1 cells and upregulated c-myb expression through binding of the -28k region. Taken together, our results provide an evidence for critical role of the -28k region in distal regulatory mechanism for c-myb gene expression during differentiation of myeloid progenitor M1 cells.

Variation of DNA Methylome of Zebrafish Cells under Cold Pressure.

DNA methylation is an essential epigenetic mechanism involved in multiple biological processes. However, the relationship between DNA methylation and cold acclimation remains poorly understood. In this study, Methylated DNA Immunoprecipitation Sequencing (MeDIP-seq) was performed to reveal a genome-wide methylation profile of zebrafish (Danio rerio) embryonic fibroblast cells (ZF4) and its variation under cold pressure. MeDIP-seq assay was conducted with ZF4 cells cultured at appropriate temperature of 28°C and at low temperature of 18°C for 5 (short-term) and 30 (long-term) days, respectively. Our data showed that DNA methylation level of whole genome increased after a short-term cold exposure and decreased after a long-term cold exposure. It is interesting that metabolism of folate pathway is significantly hypomethylated after short-term cold exposure, which is consistent with the increased DNA methylation level. 21% of methylation peaks were significantly altered after cold treatment. About 8% of altered DNA methylation peaks are located in promoter regions, while the majority of them are located in non-coding regions. Methylation of genes involved in multiple cold responsive biological processes were significantly affected, such as anti-oxidant system, apoptosis, development, chromatin modifying and immune system suggesting that those processes are responsive to cold stress through regulation of DNA methylation. Our data indicate the involvement of DNA methylation in cellular response to cold pressure, and put a new insight into the genome-wide epigenetic regulation under cold pressure.

Transcriptome comparison reveals a genetic network regulating the lower temperature limit in fish.

Transcriptional plasticity is a major driver of phenotypic differences between species. The lower temperature limit (LTL), namely the lower end of survival temperature, is an important trait delimiting the geographical distribution of a species, however, the genetic mechanisms are poorly understood. We investigated the inter-species transcriptional diversification in cold responses between zebrafish Danio rerio and tilapia Oreochromis niloticus, which were reared at a common temperature (28 °C) but have distinct LTLs. We identified significant expressional divergence between the two species in the orthologous genes from gills when the temperature cooled to the LTL of tilapia (8 °C). Five KEGG pathways were found sequentially over-represented in the zebrafish/tilapia divergently expressed genes in the duration (12 hour) of 8 °C exposure, forming a signaling cascade from metabolic regulation to apoptosis via FoxO signaling. Consistently, we found differential progression of apoptosis in the gills of the two species in which zebrafish manifested a delayed and milder apoptotic phenotype than tilapia, corresponding with a lower LTL of zebrafish. We identified diverged expression in 25 apoptosis-related transcription factors between the two species which forms an interacting network with diverged factors involving the FoxO signaling and metabolic regulation. We propose a genetic network which regulates LTL in fishes.

Global identification of the genetic networks and cis-regulatory elements of the cold response in zebrafish.

The transcriptional programs of ectothermic teleosts are directly influenced by water temperature. However, the cis- and trans-factors governing cold responses are not well characterized. We profiled transcriptional changes in eight zebrafish tissues exposed to mildly and severely cold temperatures using RNA-Seq. A total of 1943 differentially expressed genes (DEGs) were identified, from which 34 clusters representing distinct tissue and temperature response expression patterns were derived using the k-means fuzzy clustering algorithm. The promoter regions of the clustered DEGs that demonstrated strong co-regulation were analysed for enriched cis-regulatory elements with a motif discovery program, DREME. Seventeen motifs, ten known and seven novel, were identified, which covered 23% of the DEGs. Two motifs predicted to be the binding sites for the transcription factors Bcl6 and Jun, respectively, were chosen for experimental verification, and they demonstrated the expected cold-induced and cold-repressed patterns of gene regulation. Protein interaction modeling of the network components followed by experimental validation suggested that Jun physically interacts with Bcl6 and might be a hub factor that orchestrates the cold response in zebrafish. Thus, the methodology used and the regulatory networks uncovered in this study provide a foundation for exploring the mechanisms of cold adaptation in teleosts.

Small-Molecule Bcl2 BH4 Antagonist for Lung Cancer Therapy.

The BH4 domain of Bcl2 is required for its antiapoptotic function, thus constituting a promising anticancer target. We identified a small-molecule Bcl2-BH4 domain antagonist, BDA-366, that binds BH4 with high affinity and selectivity. BDA-366-Bcl2 binding induces conformational change in Bcl2 that abrogates its antiapoptotic function, converting it from a survival molecule to a cell death inducer. BDA-366 suppresses growth of lung cancer xenografts derived from cell lines and patient without significant normal tissue toxicity at effective doses. mTOR inhibition upregulates Bcl2 in lung cancer cells and tumor tissues from clinical trial patients. Combined BDA-366 and RAD001 treatment exhibits strong synergy against lung cancer in vivo. Development of this Bcl2-BH4 antagonist may provide a strategy to improve lung cancer outcome.

The Bmi-1 helix-turn and ring finger domains are required for Bmi-1 antagonism of (-) epigallocatechin-3-gallate suppression of skin cancer cell survival.

The Bmi-1 Polycomb group (PcG) protein is an important epigenetic regulator of chromatin status. Elevated Bmi-1 expression is observed in skin cancer and contributes to cancer cell survival. (-) Epigallocatechin-3-gallate (EGCG), an important green tea-derived cancer prevention agent, reduces Bmi-1 level resulting in reduced skin cancer cell survival. This is associated with increased p21(Cip1) and p27(Kip1) expression, reduced cyclin, and cyclin dependent kinase expression, and increased cleavage of apoptotic markers. These EGCG-dependent changes are attenuated by vector-mediated maintenance of Bmi-1 expression. In the present study, we identify Bmi-1 functional domains that are required for this response. Bmi-1 expression reverses the EGCG-dependent reduction in SCC-13 cell survival, but Bmi-1 mutants lacking the helix-turn-helix-turn-helix-turn (Bmi-1ΔHT) or ring finger (Bmi-1ΔRF) domains do not reverse the EGCG impact. The reduction in Ring1B ubiquitin ligase activity, observed in the presence of mutant Bmi-1, is associated with reduced ability of these mutants to interact with and activate Ring1B ubiquitin ligase, the major ligase responsible for the ubiquitination of histone H2A during chromatin condensation. This results in less chromatin condensation leading to increased tumor suppressor gene expression and reduced cell survival; thereby making the cells more susceptible to the anti-survival action of EGCG. We further show that these mutants act in a dominant-negative manner to inhibit the action of endogenous Bmi-1. Our results suggest that the HT and RF domains are required for Bmi-1 ability to maintain skin cancer cell survival in response to cancer preventive agents.

A proteasome inhibitor-stimulated Nrf1 protein-dependent compensatory increase in proteasome subunit gene expression reduces polycomb group protein level.

The polycomb group (PcG) proteins, Bmi-1 and Ezh2, are important epigenetic regulators that enhance skin cancer cell survival. We recently showed that Bmi-1 and Ezh2 protein level is reduced by treatment with the dietary chemopreventive agents, sulforaphane and green tea polyphenol, and that this reduction involves ubiquitination of Bmi-1 and Ezh2, suggesting a key role of the proteasome. In the present study, we observe a surprising outcome that Bmi-1 and Ezh2 levels are reduced by treatment with the proteasome inhibitor, MG132. We show that this is associated with a compensatory increase in the level of mRNA encoding proteasome protein subunits in response to MG132 treatment and an increase in proteasome activity. The increase in proteasome subunit level is associated with increased Nrf1 and Nrf2 level. Moreover, knockdown of Nrf1 attenuates the MG132-dependent increase in proteasome subunit expression and restores Bmi-1 and Ezh2 expression. The MG132-dependent loss of Bmi-1 and Ezh2 is associated with reduced cell proliferation, accumulation of cells in G(2), and increased apoptosis. These effects are attenuated by forced expression of Bmi-1, suggesting that PcG proteins, consistent with a prosurvival action, may antagonize the action of MG132. These studies describe a compensatory Nrf1-dependent, and to a lesser extent Nrf2-dependent, increase in proteasome subunit level in proteasome inhibitor-treated cells and confirm that PcG protein levels are regulated by proteasome activity.

Suppression of AP1 transcription factor function in keratinocyte suppresses differentiation.

Our previous study shows that inhibiting activator protein one (AP1) transcription factor function in murine epidermis, using dominant-negative c-jun (TAM67), increases cell proliferation and delays differentiation. To understand the mechanism of action, we compare TAM67 impact in mouse epidermis and in cultured normal human keratinocytes. We show that TAM67 localizes in the nucleus where it forms TAM67 homodimers that competitively interact with AP1 transcription factor DNA binding sites to reduce endogenous jun and fos factor binding. Involucrin is a marker of keratinocyte differentiation that is expressed in the suprabasal epidermis and this expression requires AP1 factor interaction at the AP1-5 site in the promoter. TAM67 interacts competitively at this site to reduce involucrin expression. TAM67 also reduces endogenous c-jun, junB and junD mRNA and protein level. Studies with c-jun promoter suggest that this is due to reduced transcription of the c-jun gene. We propose that TAM67 suppresses keratinocyte differentiation by interfering with endogenous AP1 factor binding to regulator elements in differentiation-associated target genes, and by reducing endogenous c-jun factor expression.

(-)-Epigallocatechin-3-gallate and DZNep reduce polycomb protein level via a proteasome-dependent mechanism in skin cancer cells.

Polycomb group (PcG) protein-dependent histone methylation and ubiquitination drives chromatin compaction leading to reduced tumor suppressor expression and increased cancer cell survival. Green tea polyphenols and S-adenosylhomocysteine (AdoHcy) hydrolase inhibitors are important candidate chemopreventive agents. Previous studies indicate that (-)-epigallocatechin-3-gallate (EGCG), a potent green tea polyphenol, suppresses PcG protein level and skin cancer cell survival. Inhibition of AdoHcy hydrolase with 3-deazaneplanocin A (DZNep) inhibits methyltransferases by reducing methyl group availability. In the present study, we examine the impact of EGCG and DZNep cotreatment on skin cancer cell function. EGCG and DZNep, independently and in combination, reduce the level of PcG proteins including Ezh2, eed, Suz12, Mel18 and Bmi-1. This is associated with reduced H3K27me3 and H2AK119ub formation, histone modifications associated with closed chromatin. Histone deacetylase 1 level is also reduced and acetylated H3 formation is increased. These changes are associated with increased tumor suppressor expression and reduced cell survival and are partially reversed by vector-mediated maintenance of Bmi-1 level. The reduction in PcG protein level is associated with increased ubiquitination and is reversed by proteasome inhibitors, suggesting proteasome-associated degradation.

Sodium selenite induces apoptosis in acute promyelocytic leukemia-derived NB4 cells through mitochondria-dependent pathway.

Our previous study has shown that sodium selenite can cause apoptosis in acute promyelocytic leukemia-derived NB4 cells in a caspase-dependent manner involving Deltapsim disruption and cleavage of Bcl-2, but more detailed mechanism(s) remain unclear. Here we showed that mitochondrial apoptosis signaling pathway played a vital role in apoptosis induced by sodium selenite based on the following findings: 1) cytochrome c release, activation of caspase 9, mitochondrial targeting, and oligermerization of Bax; 2) caspase 9, but not caspase 8, inhibitor could attenuate apoptosis; 3) downregulation of Bax and Bad by siRNA could delay sodium selenite-induced apoptosis. Further investigation showed that ROS was an essential inducer of deltapsim disruption and apoptosis by sodium selenite. Our findings here demonstrate that sodium selenite can induce apoptosis in NB4 cells through a mechanism involving ROS, activation of proapoptotic proteins Bad and Bax, Deltapsim disruption, release of cytochrome c, and consequent initiation of caspase cascade.

Exposure of human leukemia NB4 cells to increasing concentrations of selenite switches the signaling from pro-survival to pro-apoptosis.

Selenium at low concentrations has a chemopreventive role against cancer, while at high concentrations, selenite exerts a direct antitumor effect. However, the mechanisms behind these effects remain elusive. In this study, we found that different concentrations of selenite triggered different signal pathways in human leukemia NB4 cells. Low concentrations of selenite elicited mild endoplasmic reticulum (ER) stress and mediated cell survival by activating unfolded protein response signaling, whereas high concentrations of selenite induced severe ER stress and caused cell death by activation of the pro-apoptotic transcription factors GADD153. In addition, selenite at low concentrations activated other anti-apoptotic pathways, such as AKT and ERK, whereas high concentrations of selenite induced activation of p53 and oxidative stress, which mediated the antitumor activity of selenite by causing mitochondrial dysfunction and caspase activation. These findings uncover the molecular mechanisms of the chemopreventive and antitumor effects of different concentrations of selenite.

P53 transcription-independent activity mediates selenite-induced acute promyelocytic leukemia NB4 cell apoptosis.

Selenium, an essential trace element possessing anti-carcinogenic properties, can induce apoptosis in cancer cells. We have previously shown that sodium selenite can induce apoptosis by activating the mitochondrial apoptosis pathway in NB4 cells. However, the detailed mechanism remains unclear. Presently, we demonstrate that p53 contributes to apoptosis by directing signaling at the mitochondria. Immunofluorescent and Western blot procedures revealed selenite-induced p53 translocation to mitochondria. Inhibition of p53 blocked accumulation of reactive oxygen species (ROS) and loss of mitochondrial membrane potential, suggesting that mitochondrial p53 acts as an upstream signal of ROS and activates the mitochondrial apoptosis pathway. Selenite also disrupted cellular calcium ion homeostasis in a ROS-dependent manner and increased mitochondrial calcium ion concentration. p38 kinase mediated phosphorylation and mitochondrial translocation of p53. Taken together, these results indicate that p53 involves selenite-induced NB4 cell apoptosis by translocation to mitochondria and activation mitochondrial apoptosis pathway in a transcription-independent manner.