If you don't remember your password, you can reset it by entering your email address and clicking the Reset Password button. You will then receive an email that contains a secure link for resetting your password
If the address matches a valid account an email will be sent to __email__ with instructions for resetting your password
Department of Biochemistry, Tohoku University Graduate School of Medicine, Seiryo-machi 2-1, Sendai 980-8575, JapanCenter for Regulatory Epigenome and Diseases, Tohoku University Graduate School of Medicine, Seiryo-machi 2-1, Sendai 980-8575, Japan
To whom correspondence should be addressed: Dept. of Biochemistry, Tohoku University Graduate School of Medicine, Seiryo-machi 2-1, Sendai 980-8575, Japan. Tel.: 81-22-717-7596; Fax: 81-22-717-7598
Department of Biochemistry, Tohoku University Graduate School of Medicine, Seiryo-machi 2-1, Sendai 980-8575, JapanCenter for Regulatory Epigenome and Diseases, Tohoku University Graduate School of Medicine, Seiryo-machi 2-1, Sendai 980-8575, Japan
Ferroptosis is an iron-dependent programmed cell death event, whose regulation and physiological significance remain to be elucidated. Analyzing transcriptional responses of mouse embryonic fibroblasts exposed to the ferroptosis inducer erastin, here we found that a set of genes related to oxidative stress protection is induced upon ferroptosis. We considered that up-regulation of these genes attenuates ferroptosis induction and found that the transcription factor BTB domain and CNC homolog 1 (BACH1), a regulator in heme and iron metabolism, promotes ferroptosis by repressing the transcription of a subset of the erastin-induced protective genes. We noted that these genes are involved in the synthesis of GSH or metabolism of intracellular labile iron and include glutamate-cysteine ligase modifier subunit (Gclm), solute carrier family 7 member 11 (Slc7a11), ferritin heavy chain 1 (Fth1), ferritin light chain 1 (Ftl1), and solute carrier family 40 member 1 (Slc40a1). Ferroptosis has also been previously shown to induce cardiomyopathy, and here we observed that Bach1−/− mice are more resistant to myocardial infarction than WT mice and that the severity of ischemic injury is decreased by the iron-chelator deferasirox, which suppressed ferroptosis. Our findings suggest that BACH1 represses genes that combat labile iron-induced oxidative stress, and ferroptosis is stimulated at the transcriptional level by BACH1 upon disruption of the balance between the transcriptional induction of protective genes and accumulation of iron-mediated damage. We propose that BACH1 controls the threshold of ferroptosis induction and may represent a therapeutic target for alleviating ferroptosis-related diseases, including myocardial infarction.
Relationships among cortical glutathione levels, brain amyloidosis, and memory in healthy older adults investigated in vivo with 1H-MRS and Pittsburgh compound-B PET.
). Considering the involvement of lipid hydroperoxide, ferroptosis may be executed at the edge of the oxidative stress response. Therefore, ferroptosis may be a regulated process involving the oxidative stress response. However, the regulatory mechanism underlying ferroptosis has yet to be elucidated in full. In particular, it is not known at all how the regulatory genes of ferroptosis are transcriptionally controlled in induction of ferroptosis.
is a heme-binding transcription factor required for the proper regulation of the oxidative stress response and metabolic pathways related to heme and iron (
). BACH1 represses Hmox1 encoding heme oxygenase-1 (HO-1), Fth1 and Ftl1 encoding ferritin proteins, Gclm and Gclc encoding glutamate-cysteine ligase modifier and catalytic subunits, and other genes involved in the oxidative stress response (
Bach1 repression of ferritin and thioredoxin reductase1 is heme-sensitive in cells and in vitro and coordinates expression with heme oxygenase1, β-globin, and NADP(H) quinone (oxido) reductase 1.
). We hypothesized that BACH1 might regulate ferroptosis by inhibiting the expression of these genes. In addition, although BACH1 is involved in the exacerbation of various diseases involving oxidative stress, such as ischemic heart disease (
), it is unknown whether or not they are attributable to ferroptosis. We considered that BACH1 may exacerbate the severity of these diseases through ferroptosis.
To understand the regulatory mechanism underlying ferroptosis, we analyzed the transcriptome response in ferroptotic cells with RNA-Sequence (RNA-Seq). We also examined whether or not BACH1 was involved in the regulation of ferroptosis by comparing ferroptosis and the expression of ferroptosis-induced genes between wildtype (WT) and Bach1−/− mouse embryonic fibroblasts (MEFs). Furthermore, we assessed the influence of BACH1 and ferroptosis on the severity of acute myocardial infarction (AMI) in model mice. We found that BACH1 promoted ferroptosis by directly repressing genes involved in the synthesis of glutathione (GSH) and sequestration of free labile iron. BACH1 also increased the severity of AMI, which was mitigated by the iron chelator deferasirox (DFX). Our findings highlight the coordinated transcriptional response and its regulation by BACH1 upon ferroptosis.
Results
Regulatory genes of ferroptosis show compensatory up-regulation in ferroptotic cells
Changes in metabolic and biological processes occur in ferroptotic cells (
), and carried out RNA-seq analyses. Genes related to oxidative stress and iron metabolism showed significant induction in their expression (Fig. 1A). Many of these genes, for example, Slc7a11, Gclm, Gclc, are considered to possess inhibitory effects on ferroptosis (
). Therefore, ferroptosis accompanies the induction of genes that can restrict the execution of ferroptosis.
Figure 1Regulatory genes of ferroptosis are up-regulated with decreased BACH1 protein in the induction of ferroptosis.A, RNA-seq was performed in WT MEFs (9th passage: P9) with only DMSO (DMSO group) or DMSO + 3 μm erastin (Erastin group) for 24 h. A heat map of gene expression profiles shows the genes registered to map04216 (Ferroptosis pathway) of Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway map. The genes were arranged from the bottom in the order of fold-change of erastin group to DMSO group. n = 3 per group. B and C, Western blotting for BACH1, HO-1, and GAPDH of WT MEFs (7th passage: P7, 11th passage: P11) exposed to erastin for 12 h. Representative image (B) and the intensity of bands (C). D and E, Western blotting for BACH1 and GAPDH in WT MEFs (8th passage) exposed to erastin and 25 μm MG132 for 6 h. Representative image (D) and the intensity of bands (E). In A, p value by the differential expression analysis performed on edge R. In C and E, error bars represent S.D. p value by Tukey's test after two-way ANOVA.
). Indeed, the amount of BACH1 protein was decreased in MEFs exposed to erastin, which was accompanied by the induction of Hmox1 (Fig. 1, B and C, and Fig. S1). Although BACH1 protein accumulated after proteasome inhibition by MG132, the degradation of BACH1 in response to erastin was not completely blocked (Fig. 1D and E). These observations suggest that BACH1 protein is degraded by proteasome-dependent and -independent mechanisms during ferroptosis. With the reduction in BACH1 protein, the production of its mRNA was induced (Fig. S1), suggesting the presence of feedback regulation of BACH1.
These observations suggest that, when cells are exposed to erastin, the expression of genes that counteract ferroptosis is induced in part by a reduction in BACH1 protein and that the amount or activity of BACH1 and the kinetics of its feedback regulation may influence ferroptosis by suppressing this counteracting subprogram of ferroptosis.
BACH1 promotes ferroptosis
To clarify whether or not BACH1 regulates ferroptosis, we treated WT and Bach1−/− MEFs with erastin, stained them with propidium iodide (PI) and annexin V, and compared the cell death by a flow cytometry analysis (Fig. S2). Bach1−/− MEFs showed less cell death in response to erastin than WT cells (Fig. 2, A and B). When the erastin-treated MEFs were observed with a transmission electron microscope, mitochondrial spheroids, which reflect oxidative stress (
), were frequently observed in WT MEFs (Fig. 2C). Although mitochondrial spheroids were not clear in Bach1−/− MEFs exposed to erastin, condensed matrix structures, which are also observed as parts of mitochondrial spheroids (
), were observed (Fig. 2C). These results may reflect lower oxidative stress and higher resistance to erastin of Bach1−/− MEFs than WT MEFs. The cell death in our experiments was inhibited by the iron chelator deferoxamine (DFO) (Fig. 2, D and E), confirming that this death was ferroptosis. These results showed that BACH1 promoted ferroptosis in MEFs.
Figure 2BACH1 promotes ferroptosis.A and B, optical microscope image (A) and quantification of cell death by flow cytometer (B) of WT and Bach1−/− MEFs (11th passage: P11) exposed to erastin for 24 h. Scale bars in A represent 100 μm. C, transmission electron microscope image of WT and Bach1−/− MEFs (P9) exposed to erastin for 10 h. The portion surrounded by dashed line: a mitochondrial spheroid. Arrow: mitochondrial inner membrane-matrix condensation. Scale bars represent 200 nm. D and E, optical microscope image (D) and quantification of cell death by flow cytometer (E) of WT and Bach1−/− MEFs (P12) exposed to erastin and DFO for 24 h. Scale bars in D represents 200 μm. A, B, D, and E, representative of three independent experiments. Error bars of B represent S.D. The box and whisker plots of E show the 25th and 75th percentile quartiles and median values (center black line) and maximum and minimum values of the data. p value of B by t test. p value of E by three-way ANOVA.
It should be noted that Bach1−/− MEFs lost their resistance to ferroptosis under high doses of erastin (Fig. 2, A and B). Although decreased BACH1 by erastin is considered as a protective response to ferroptosis (Fig. 1, A–C), there is also its limitation to the inhibition of cell death. However, at least at the early stage, the reduction in BACH1 protein may be a part of the ferroptosis program, and BACH1 may set the threshold for ferroptosis. Execution of ferroptosis may be determined by the basal amount of BACH1 and how rapidly it is degraded in response to ferroptosis inducers.
BACH1 represses the expression of genes involved in the GSH synthesis pathway
) and the genes involved in the pathway of GSH synthesis were up-regulated as BACH1 protein decreases (Fig. 1, A–C). Therefore, BACH1 may decrease GSH by repressing the expression of these genes. To investigate this possibility, we measured the intracellular GSH concentrations in WT and Bach1−/− MEFs. Because the sensitivity to ferroptosis decreased as we passaged MEFs (Fig. S3, A and B), we performed the measurement of GSH and the subsequent analysis of gene expression at three different stages of primary culture (7th, 9th, and 11th passage; P7, P9, and P11). The amounts of GSH at P7 and P9 were significantly higher in Bach1−/− MEFs than in WT cells and a similar trend was observed at P11 (Fig. 3A), suggesting that BACH1 promoted ferroptosis by reducing GSH within cells.
Figure 3BACH1 decreases GSH by repressing the transcription of Slc7a11, Gclm, and Gclc.A, intracellular concentration of GSH in WT and Bach1−/− MEFs (7th, 9th, and 11th passage: P7, P9, and P11) by UHPLC/MS/MS. B, ChIP-seq analysis of the binding of BACH1, MAFK for gene regions of Slc7a11, Gclm, and Gclc in M1 cells. C, qPCR analysis for Slc7a11, Gclm, and Gclc mRNA relative to Actb mRNA in WT and Bach1−/− MEFs (P7, P9, P11). n = 3 of independent lots of MEFs per genotype. D and E, Western blotting for BACH1, SLC7A11, GCLM, GCLC, and GAPDH in WT and Bach1−/− MEFs (P10) exposed to 25 μm MG132. Representative image (D) and the intensity of bands (E). The box and whisker plots of A show the 25th and 75th percentile quartiles and median values (center black line) and maximum and minimum values of the data. Error bars of C and E represent S.D. p value of A, C, and E by t test.
The double knockout of Bach1 and Bach2 in mice reveals shared compensatory mechanisms in regulating alveolar macrophage function and lung surfactant homeostasis.
), we found peaks of BACH1 and its partner MAFK in the regulatory regions of genes encoding molecules for GSH synthesis, including Gclm, Gclc, and Slc7a11 (Fig. 3B). Furthermore, by comparing the expression of these genes in WT and Bach1−/− MEFs by quantitative PCR (qPCR), the expression of all of these genes was confirmed to be higher in Bach1−/− MEFs than in WT cells (Fig. 3C). These results suggested that BACH1 bound to the regulatory regions of these genes to repress their expression.
A comparison of the protein amounts of SLC7A11, GCLM, and GCLC in MEFs by Western blotting revealed that more GCLM protein was present in Bach1−/− MEFs than in WT cells (Fig. 3, D and E, and Fig. S4, A and B). Although the amounts of SLC7A11 protein were similar in WT and Bach1−/− MEFs (Fig. 3, D and E, and Fig. S4, A and B), more SLC7A11 protein was present in Bach1−/− MEFs than in WT cells when they were treated with proteasome inhibitor MG132 (Fig. 3, D and E). These observations suggest that the amount of SLC7A11 protein is further tuned by proteasomal-mediated degradation. There were no marked differences in the amount of GCLC protein with or without MG132 (Fig. 3, D and E, and Fig. S4, A and B). BACH1 may affect the expression of GCLC protein under certain circumstances. Given these results, we surmised that BACH1 decreased the amount of GSH by repressing the expression of Gclm and in part Slc7a11.
Knockdown of the genes in GSH synthesis pathway and Hmox1 increases the sensitivity of ferroptosis
We next examined whether or not the resistance of Bach1−/− MEFs against ferroptosis was actually dependent on the increased expression of the genes involved in the GSH synthesis pathway. Although it is not always statistically significant, knockdown of any of Slc7a11, Gclm, and Gclc resulted in slight but reproducible increases in ferroptosis in both WT and Bach1−/− MEFs (Figs. S5, A–C; S6, A–C; and S7, A–C). These results show that the genes involved in the GSH synthesis pathway have inhibitory effects against ferroptosis and suggest that BACH1 promotes ferroptosis by repressing their expression. We next examined the effect of knockdown of Hmox1. WT MEFs became more sensitive to ferroptosis by knockdown of Hmox1 than cells with control knockdown (Figs. S5D, S6D, and S7D). We thus concluded that HO-1 works as an inhibitor of ferroptosis under our experimental conditions. However, the effect of HO-1 to accelerate ferroptosis has also been reported (
). The function of HO-1 in ferroptosis might differ depending on the situations of cells. Importantly, knockdown of Slc7a11, Gclm, Gclc, or Hmox1 did not decrease the observed differences in ferroptosis between WT and Bach1−/− MEFs (Fig. S7). These results suggest that the role of BACH1 in promoting ferroptosis depends on the repression of multiple genes involved in ferroptosis.
BACH1 accelerates ferroptosis by suppressing labile iron metabolism
To explore other target genes of BACH1 in the regulation of ferroptosis, we examined genes involved in the regulation of iron metabolism (Fth1, Ftl1, Slc40a1, Tfrc, Mfn2, and Fxn), heavy metal stress (Mt1), and lipoperoxidation (Gpx4). Some of these genes were up-regulated in response to erastin (see Fig. 1A). Among these genes, ferritin genes (Fth1 and Ftl1) and the ferroportin gene (Slc40a1) were dramatically up-regulated in Bach1−/− MEFs (Fig. 4A), and binding peaks of BACH1 and MAFK were observed near their regulatory regions (Fig. 4B). In contrast, the expression of Tfrc, Mfn2, Fxn, Mt1, and Gpx4 was only mildly increased in Bach1−/− MEFs (Fig. S8A). There were no strong binding peaks of BACH1 or MAFK in the regulatory regions of these genes (Fig. S8B). Considering that both ferritin and ferroportin reduce the availability of free labile iron and are known to inhibit ferroptosis (
), these results suggest that BACH1 promotes ferroptosis by repressing the transcription of ferritin and ferroportin genes.
Figure 4BACH1 increases mitochondrial labile iron during ferroptosis by repressing the transcription of genes of ferritin and ferroportin.A, qPCR analysis for Fth1, Ftl1, and Slc40a1 mRNA relative to Actb mRNA in WT and Bach1−/− MEFs (7th, 9th, and 11th passage: P7, P9, and P11). n = 3 of independent lots of MEFs per genotype. B, ChIP-seq analysis of the binding of BACH1, MAFK for gene regions of Fth1, Ftl1, and Slc40a1 in M1 cells. C and D, WT and Bach1−/− MEFs (P8) were exposed to erastin for 24 h and afterward, stained by Mito-FerroGreen, a fluorophore. Representative flow cytometry images (C) showing the strategy that was implemented for the sorting of Mito-FerroGreen positive cells. The positive cells were judged as mitochondrial-labile iron positive cells. Quantification of mitochondrial-labile iron (D). MFI, mean fluorescence intensity. E, conceptual diagram. Error bars of A and D represent S.D. p value of A by t test. p value of D by Tukey's test after two-way ANOVA.
). Bach1−/− MEFs contained more mitochondrial labile iron than WT without erastin (Fig. 4, C and D). It is probably due to the effects of HO-1 overexpression. Although prominent increases in labile iron and cell death were observed in WT MEFs in response to erastin, these alterations were much less in Bach1−/− MEFs (Fig. 4, C and D, and Fig. S9, A and B). These results suggest BACH1 increases mitochondrial labile iron by repressing the transcription of ferritin and ferroportin during ferroptosis. These findings, along with the regulation of GSH synthesis pathway by BACH1, suggest that BACH1 accelerates ferroptosis by decreasing the intracellular activity of GSH and increasing the oxidative activity of labile iron (Fig. 4E).
BACH1 aggravates acute myocardial infarction by promoting ferroptosis
Finally, we tried to examine whether or not the promotion of ferroptosis by BACH1 is involved in pathological changes in vivo. It has long been known that iron and lipid peroxidation are involved in oxidative stress-related cell death, including ischemia-reperfusion injury (
) (Fig. 5A). In this model, Bach1−/− mice showed less severe injuries than WT mice as judged by the post-operative survival rate and an evaluation of the cardiac function with echocardiography (Fig. 5, B and C,Fig. S10, A–C, and Movie S1, A–D). The infarct area on pathological specimens was also smaller in Bach1−/− mice than in WT mice (Fig. 5, D and E). These results suggest that BACH1 exacerbates the pathology of AMI.
Figure 5BACH1 aggravates AMI.A, experimental process. B, Kaplan-Meier curve of each group. C, left ventricular fractional shortening (LVFS) on echocardiogram. D and E, mice was dissected after 9 weeks from operation. D, representative photographs of heart sections stained with elastica Masson staining. E, infarct size to left ventricular section. Scale bars in D represent 2 mm. F, mice was dissected next day from operation. Transmission electron microscope image of normal and infarct area of hearts of mice next day from operation. Orange arrow: normal mitochondria. Yellow arrow: mitochondrial spheroid-like abnormality. Scale bars represent 500 nm. Error bars of C represent S.D. The box and whisker plots of E show the 25th and 75th percentile quartiles and median values (center black line) and maximum and minimum values of the data. p value of B by log-rank test between WT (AMI) and Bach1−/− (AMI). p value of C by Tukey-Kramer method after two-way ANOVA. p value of E by t test.
To investigate whether or not ferroptosis is involved in the pathology, we observed the myocardial infarct regions using a transmission electron microscope. Mitochondrial spheroid-like abnormality was observed in both WT and Bach1−/− mice (Fig. 5F). This alteration was consistent with the findings of MEFs exposed to erastin (Fig. 2C). We then investigated whether or not the pathological changes could be improved by administering DFX (Fig. 6A), which is a clinically used iron chelator. First, we confirmed that it inhibited ferroptosis in MEFs (Fig. S10, D and E). Although there was no improvement in the survival rates in WT or Bach1−/− mice (Fig. 6B), an improvement in the cardiac function on echocardiography was observed in the DFX group, which was more prominent in the WT mice than Bach1−/− mice (Fig. 6, C and D, and Fig. S10, F–K). The DFX group of WT mice showed a reduction in the infarct area; however, no such effect was noted in Bach1−/− mice (Fig. 6, E and F). These results suggest that BACH1 exacerbates the pathology of AMI by promoting ferroptosis.
Figure 6An iron chelator DFX alleviates AMI.A, experimental process. B, Kaplan-Meier curve of each group. C and D, left ventricular fractional shortening (LVFS) of WT mice (C) and Bach1−/− mice (D) on echocardiogram. E and F, mice were dissected after 9 weeks from operation. E, representative photographs of heart sections stained with elastica Masson staining. F, infarct size to left ventricular section. Scale bars in E represent 2 mm. Error bars of C and D represent S.D. The box and whisker plots of F show the 25th and 75th percentile quartiles and median values (center black line) and maximum and minimum values of the data. p value of C, D, and F by t test.
), how their expression is regulated during ferroptosis remains unclear. In this study, we found that many of the inhibitory genes of ferroptosis were coordinately up-regulated upon induction of ferroptosis with erastin (Fig. 1A). Such a coordinated response may be a mechanism for restricting ferroptosis. We further showed that BACH1 directly counteracted this coordinated response of genes, including Hmox1, Slc7a11, Gclm, Gclc, Fth1, Ftl1, and Slc40a1 (Figs. 3, B and C, and 4, A and B), which are involved in the metabolism of GSH or labile iron. The protein amount of BACH1 was reduced upon the induction of ferroptosis (Fig. 1B). Our results suggest that BACH1 is degraded in part by proteasome during ferroptosis (Fig. 1, D and E). BACH1 is known to be degraded by proteasome under oxidative stress including heme (
). Therefore, reduction of the BACH1 protein level may trigger the coordinated induction of these genes as a subprogram of the initial phase of ferroptosis program. Cells can then integrate distinct signals leading to BACH1 degradation, and thus judge whether or not they should undergo ferroptosis. Thus, BACH1 sets the threshold for whether or not ferroptosis occurs in response to lipid peroxide synthesized.
NRF2 is known to activate some of the genes that are repressed by BACH1, including Hmox1, Slc7a11, Gclm, and Gclc (
Cobalt induces heme oxygenase-1 expression by a hypoxia-inducible factor-independent mechanism in Chinese hamster ovary cells: regulation by Nrf2 and MafG transcription factors.
). Therefore, ferroptosis execution may depend on the initial amounts and kinetics of the induction or reduction of these transcription factors. This mechanism may extend our understanding of the regulation of ferroptosis, wherein ferroptosis is a cell death programmed at the level of the gene regulatory network.
We showed that GSH was higher in Bach1−/− MEFs than WT cells (Fig. 3A). Our results strongly suggest that BACH1 decreases intracellular GSH by repressing the expression of Gclm, Gclc, and Slc7a11 (Fig. 3, B and C). Indeed, the protein amount of GCLM was higher in Bach1−/− MEFs than in WT cells (Fig. 3, D and E, and Fig. S4, A and B). However, the protein amounts of GCLC and SLC7A11 were similar between WT and Bach1−/− MEFs (Fig. 3, D and E, and Fig. S4, A and B). Cells may have additional mechanisms to tune strictly the protein amounts of GCLC and SLC7A11, managing the intracellular GSH amount and maintaining homeostasis. We found that SLC7A11 was further regulated by proteasomal degradation (Fig. 3, D and E). This observation suggests that the decision to undergo ferroptosis may be made based upon whether or not cells can induce efficiently inhibitory proteins like SLC7A11. Cells with higher amounts of SLC7A11 may likely be protected from ferroptosis. Gclc and Slc7a11 may be critical factors for cells, with the transcriptional regulation by BACH1 and additional layers of regulation, although these points will need to be explored in further studies.
Reports on the function of HO-1 are conflicting, with studies conversely describing it to promote or inhibit ferroptosis (
). These discrepant findings may be because HO-1 degrades prooxidant heme to produce not only the radical scavengers biliverdin and bilirubin but also free iron that mediates ferroptosis through Fenton reaction (
). Therefore, to allow HO-1 to function effectively as an anti-oxidative stress enzyme, it is essential to suppress the reactivity of labile iron derived from heme. We showed that BACH1 represses the expression of the genes of ferritin and ferroportin (Fig. 4, A and B), both of which reduce the intracellular availability of labile iron. By increasing the expression of not only HO-1 but also ferritin and ferroportin during the induction of ferroptosis (Fig. 1A), the prooxidant activities of heme and heme-derived free iron can be suppressed efficiently, thus protecting cells from ferroptosis. Conversely, BACH1 represses the expression of ferritin and ferroportin in addition to HO-1, thus effectively promoting ferroptosis (Fig. 4E). Indeed, we showed Bach1−/− MEFs were more resistant to the elevation of mitochondrial labile iron than WT MEFs during ferroptosis (Fig. 4, C and D). Based on the present and previous findings, we proposed a model in which BACH1 accelerates ferroptosis by suppressing two major intracellular counteracting mechanisms against ferroptosis: the GSH synthesis pathway and the system for the sequestration of labile iron (Fig. 4E).
In addition, we showed that ferroptosis was involved in the pathology of not only ischemia-reperfusion injury (
) but also AMI. The severity of AMI was improved by the iron chelator, DFX particularly in WT mice (Fig. 6, C–F). The peripheral areas of AMI are naturally reperfused by angiogenesis, where ferroptosis is likely induced. Our results here suggest that the therapeutic effect of DFX is expected in AMI and ischemia-reperfusion injury. Necroptosis is also reportedly involved in cardiac ischemic disease (
). Therefore, the double inhibition of ferroptosis and necroptosis may lead to a more effective treatment of AMI. In addition, this study suggests that Bach1−/− mice are more resistant to AMI than WT mice because of their lower rate of ferroptosis than in WT mice (Figure 5, Figure 6). BACH1 may be a potential therapeutic target of AMI in the future.
Ferroptosis is thought to play a major role in cancer suppression (
). Our results suggest that cancer cells may acquire resistance against ferroptosis by decreasing BACH1 protein, thus eluding elimination by ferroptosis. We previously reported that BACH1 promotes the proliferation of MEFs transformed with H-RasV12 and their tumor formation in a mouse transplantation model (
). BACH1 is therefore considered to have dual functions in cancers: promoting cell proliferation and cell death through ferroptosis. Cancer cells may adapt to their surrounding environment by changing the expression of BACH1; cancer cells may highly express BACH1 during stages of proliferation and metastasis but may reduce their levels of BACH1 under stress conditions, such as toxicity due to anti-cancer drugs. Such flexibility in the amount of BACH1 protein expressed may enhance the malignancy of cancer cells. Therapy that targets this flexibility, such as the down-regulation of BACH1 in response to erastin, may expand the field of potential cancer treatments.
Experimental procedures
Mice
The generation of Bach1−/− mice on the C57BL/6J background was described previously (
). Mice 13 weeks of age were analyzed for models of AMI. Animals were euthanized by cervical dislocation under anesthetic inhalation overdose with isoflurane before anatomy. These mice were bred at the animal facility of Tohoku University. Mice were housed under specific pathogen-free conditions. All experiments performed in this study were approved by the Institutional Animal Care and Use Committee of the Tohoku University Environmental and Safety Committee.
Mice models of AMI
Induction of AMI was performed as described previously (
). The mice were subjected to ligation of the proximal LAD to induce AMI. They were randomly assigned to sham or AMI (Fig. 5A), DMSO, or DFX groups (Fig. 6A). To follow up the time course of LV function after AMI, we performed transthoracic two-dimensional echocardiography. For histological analysis and analysis with transmission electron microscope, the heart was divided along the short axis at the center of the infarct.
Histopathological analysis
Excised hearts were fixed with 4% paraformaldehyde for histological and immunohistochemical examination. After 24–48 h of fixation and dehydration through increasing concentrations of ethanol, the tissue specimens were embedded in paraffin and sliced at 3 μm in thickness. The sections were used for Elastica-Masson staining. The extent of infarct area was calculated as a rate of fibrotic area using the following formula: fibrotic area/(LV free wall + interventricular septum) × 100 (%) with use of Photoshop software (Adobe).
Transmission electron microscopy
Cells and hearts were treated in 2.5% glutaraldehyde in 0.1 m cacodylate buffer (pH 7.4) for at least 24 h, and washed with 0.1 m cacodylate buffer 4 times and then treated with 1% OsO4 in 0.1 m cacodylate buffer for 90 min. After dehydration through an ethanol series (50–100% ethanol), cells were embedded in Epon resin. Thin sections were cut with a microtome (Leica EM UC-7), stained with 2% uranyl acetate and 0.4% lead citrate, and examined and photographed under a transmission electron microscope (Hitachi High-Technologies H-7600).
Isolation and culture of MEFs
MEFs were derived from 13.5-day-old embryos of WT or Bach1−/− mice. Following removal of the head and organs, embryos were rinsed with PBS (Nissui, Tokyo, Japan), minced and digested with trypsin (0.05% (v/v) solution containing 0.53 mm EDTA) (Gibco), and 1.8 mg/ml of DNase I (Roche, Basel, Switzerland) in PBS and incubated for 60 min at 37 °C. Trypsin was inactivated by addition of Dulbecco's modified Eagle's medium with high glucose (Gibco) containing 10% (v/v) fetal bovine serum (Sigma-Aldrich), 1× minimal essential medium nonessential amino acids (Gibco), and 0.1 mm 2-mercaptoethanol (Sigma-Aldrich). MEFs from a single embryo were plated into a 100-mm diameter culture dish and incubated at 37 °C in 3% oxygen (1st passage: P1). MEFs from embryos of homosexual littermates were mixed at 2nd passage (P2) and stocked.
MEFs were maintained at 37 °C in culture medium (Dulbecco's modified Eagle's medium with (Gibco) containing 10% fetal bovine serum (Sigma-Aldrich), 1× minimal essential medium nonessential amino acids (Gibco), penicillin/streptomycin (100 units/ml and 100 μg/ml each) (Gibco) and 0.1 mm 2-mercaptoethanol (Sigma-Aldrich)) in 3% oxygen for experiments. The number of passages were recorded for each lot of MEFs. From 5th to 11th passage MEFs were used for all experiments.
Reagents
Erastin, DMSO, and DFO were purchased from Sigma-Aldrich. MG132 was purchased from Calbiochem (San Diego, CA). DFX was transferred as raw material from Novartis Pfarma (Basel, Switzerland). S-Adenosylmethionine-13C5,15N, and GSH-13C2,15N were purchased from Taiyo Nissan Corp. (Tokyo, Japan) and used as internal standard for MS. Methanol, acetonitrile, and ammonium hydroxide for MS were purchased from Kanto Chemical (Tokyo, Japan). Ammonium bicarbonate (1 mol/liter) for MS was purchased from Cell Science and Technology Institute, Inc. (Miyagi, Japan). Formic acid for MS was purchased from Wako Pure Chemical Industries (Osaka, Japan).
Sample preparation for UHPLC/MS/MS
MEFs (3–8 × 106 cells for each lot) were suspended in 100 μl of methanol containing the internal standards (0.2 μg/ml of S-adenosylmethionine-13C5,15N for positive ion mode (Pos) and 1 μg/ml of GSH-13C2,15N for negative ion mode (Neg)), and were homogenized by mixing for 30 s followed by sonication for 10 min. After centrifugation at 16,400 × g for 20 min at 4 °C followed by deproteinization, 3 μl of each extract was analyzed by ultra HPLC triple quadrupole MS (UHPLC/MS/MS).
UHPLC/MS/MS analysis
The UHPLC/MS/MS analysis was performed on an AcquityTM Ultra Performance LC I-class system (Waters Corp., Milford, UK) interfaced to a Waters Xevo TQ-S MS/MS system equipped with electrospray ionization. The MS/MS was performed using the multiple reaction monitoring mode with a scan time of 1 ms for each compound. The transitions of the precursor ion to the product ion, cone voltage, and collision energy are listed in Table S1. The other settings are as follows: 3.5 kV (Pos) or 2.5 kV (Neg) capillary voltage, 30 V cone voltage, 50 V source offset, source temperature at 150 °C, 150 liter/h of cone gas (N2) flow rate, desolvation temperature at 450 °C, 1000 liter/h of desolvation gas flow, 0.15 min/ml of collision gas flow, 7.00 bar nebulization gas (N2) flow. LC separation was performed as described before (
), using a normal-phase column (ZIC-pHILIC; 100 × 2.1 mm inner diameter, 5-μm particle size; Sequant, Darmstadt, Germany) with a gradient elution using solvent A (10 mmol/liter of NH4HCO3, adjusted to pH 9.2 using ammonia solution) and B (acetonitrile) at 300 μl/min: 99 to 70% B from 0.5 to 4.0 min, 70 to 1% B from 4.0 to 6.5 min, 1% B for 2.5 min, and 99% B for 9 min until the end of the run. The oven temperature was 20 °C. The data were collected using the MassLynx version 4.1 software (Waters Corp.) and the ratio of the peak area of analyte to the internal standard was analyzed by Traverse MS (Reifycs Inc., Tokyo, Japan).
RNA interference
All siRNAs (siControl: Stealth RNAiTM siRNA Negative Control, Med GC; siGclm #1: MSS204722; siGclm #3: MSS204724; siSlc7a11 #1: MSS218649; siSlc7a11 #2: MSS218650; siGclc #2: MSS204720; siGclc #3: MSS204721; siHmox1 #1: MSS247281; siHmox1 #3: MSS274857) were obtained from Invitrogen. Sequences of the siRNAs are described in Table S2. 2 × 106 cells of MEFs were transfected with 1.2 nm siRNAs using Amaxa Nucleofector II (Lonza, Basel, Switzerland) and MF 1 Nucleofector kit (Lonza) according to the manufacturer's protocols. After transfection, MEFs were passaged to dishes or culture plate with culture medium.
Western blotting
Cells were trypsinized, pelleted, and washed twice in PBS. Cells were lysed beyond 5 min in SDS sample buffer (62.5 mm Tris-HCl (pH 6.8), 1% (v/v) 2-mercaptoethanol, 1% (w/v) SDS, 10% (w/v) glycerol, 0.02% (w/v) bromphenol blue). Lysates were resolved on 7.5–10% SDS-PAGE gels and transferred to polyvinylidene difluoride membranes (Millipore). The antibody for detection of β-actin (sc-1616) was purchased from Santa Cruz Biotechnology (Santa Cruz, CA). The antibody for detection of HO-1 (ADI-SPA-896) was purchased from Enzo Life Science (New York). The antibodies for GAPDH (ab8245), Gclc (ab53179), and Gclm (ab124827) were purchased from Abcam (Cambridge, UK). The antibody for Slc7a11 (119-11215) was purchased by RayBiotech (Norcross, GA). The antibody for BACH1 was described previously (
). All bands to be compared were quantified on the same image and were within the linear range of detection of the software.
Quantitative PCR with reverse transcription
Total RNA was purified with RNeasy plus micro kit or RNeasy plus mini kit (Qiagen, Hilden, Germany). Complementary DNA was synthesized by a SuperScript III First-Strand Synthesis System (Invitrogen). Quantitative PCR was performed using LightCycler Fast Start DNA Master SYBR Green I, and LightCycler nano (Roche) or LightCycler 96 (Roche). mRNA transcript abundance was normalized to that of Actb. Sequences of the qPCR primers are described in Table S3.
Administration of erastin and cell death assessment by flow cytometry
Before administration of erastin, the medium was exchanged to the experimental medium (culture medium without 2-mercaptoethanol and penicillin/streptomycin). Erastin was dissolved in DMSO and administered to experimental medium with DMSO. The concentration of DMSO was adjusted among each samples. Cell death was assessed 24 h after administration of erastin. PI and Annexin V staining were used for assessment of cell death. APC-Annexin V was purchased from BD Bioscience (Franklin Lakes, NJ). MEFs were stained by APC-Annexin V according to the manufacturer's protocols. PI was added to aliquot (1 μg/ml) before flow cytometry. The MEFs were sorted with a FACS Aria II (BD) and analyzed by FlowJo software (Tree Star, Ashland, OR). MEFs of positive, of whether at least Annexin V or PI, was assessed as dead cells. The gating strategy for assessing dead cells (Figs. 2, B and E; 4, C and D; and Figs. S3B, S7, A–D, S9B, and S10E) was shown in Fig. S2.
ChIP-Seq
We used ChIP-seq data of BACH1 and MAFK in the M1 cell line from GEO (Gene Expression Ominibus) data set GSE79139 that deposited for our previous report (
The double knockout of Bach1 and Bach2 in mice reveals shared compensatory mechanisms in regulating alveolar macrophage function and lung surfactant homeostasis.
) was used according to the manufacturer's protocol. MEFs on 12-well plates (Corning, NY) were washed three times with Hank's balanced salt solution (HBSS) (Gibco) to remove the residual medium. Then MEFs were treated with 5 μm Mito-FerroGreen with HBSS for 30 min at 37 °C. After incubation, Mito-FerroGreen was removed by washing twice with HBSS. MEFs were stained by APC-Annexin V and PI, and sorted with a FACS Aria II (BD Bioscience) and analyzed by FlowJo software (Tree Star). MEFs that were negative for both Annexin V and PI were sorted as alive cells (Fig. 4, C and D). The gating strategy was shown in Fig. S2.
RNA-Seq
Total RNA was purified using an RNeasy plus mini kit (Qiagen). To remove ribosomal RNA (rRNA), 4 μg of the total RNA was treated with a GeneRead rRNA Depletion kit (Qiagen) and then with an RNeasy MiniElute kit (Qiagen) for cleanup. For fragmentation, 100 ng of the rRNA-depleted RNA was incubated at 95 °C for 10 min and purified by a Magnetic Beads Cleanup Module (Thermo Fisher Scientific, Carlsbad, CA). The libraries were constructed with an RNA-seq library kit version 2 (Thermo Fisher Scientific) on ABI library builder (Thermo Fisher Scientific), and was barcoded with Ion Xpress RNA-seq BC primer (Thermo Fisher Scientific). The library fragments with a size range of 100–200 bp were selected with Agencourt AMPure XP beads (Beckman Coulter, Brea, CA). Templates were prepared on the Ion Chef system using an Ion PI Hi-Q Chef kit (Thermo Fisher Scientific) and sequencing was performed on an Ion Proton system using the Ion PI Hi-Q sequencing kit (Thermo Fisher Scientific) and the PI v3 chip (Thermo Fisher Scientific). The sequence data were obtained as fastq files. The sequence data were aligned to reference hg19 using the RNASeq Analysis plugin from Ion torrent suite software (Thermo Fisher Scientific). Mapped reads were counted for each gene using HTSeq version 0.9.1 htseq-count. The differential expression analysis was performed on edge R version 3.16.5 after removal of low count lead genes using three biological replicates for each condition (less than 5 leads per gene in the sample and counts per million mapped reads (cpm) of 1 or less).
Statistics
For all experiments, differences of data sets were considered statistically significant when p values were lower than 0.05. Statistical comparisons were performed using the t test in comparison between the two groups, and one-, two-, or three-way ANOVA followed by Tukey's test or Tukey-Kramer method in comparison among multiple groups. For the t test, Student's t test was used when the S.D. of the groups was not significantly different by f test. Welch's t test was used when the S.D. of the groups was significantly different by f test.
Author contributions
H. N., M. M., and K. I. conceptualization; H. N., H. K., K. S., M. S., and Y. I. resources; H. N. formal analysis; H. N. validation; H. N., M. M., and D. S. investigation; H. N., M. M., T. S., D. S., H. S., and K. I. methodology; H. N. writing-original draft; M. M. and K. I. funding acquisition; K. I. supervision; K. I. project administration; K. I. writing-review and editing.
Acknowledgments
We thank members of the Department of Biochemistry, Tohoku University Graduate School of Medicine, for discussions and support and the Biomedical Research Core of Tohoku University Graduate School of Medicine for technical support. We thank Novartis for the raw material transfer of DFX for this study.
Relationships among cortical glutathione levels, brain amyloidosis, and memory in healthy older adults investigated in vivo with 1H-MRS and Pittsburgh compound-B PET.
Bach1 repression of ferritin and thioredoxin reductase1 is heme-sensitive in cells and in vitro and coordinates expression with heme oxygenase1, β-globin, and NADP(H) quinone (oxido) reductase 1.
The double knockout of Bach1 and Bach2 in mice reveals shared compensatory mechanisms in regulating alveolar macrophage function and lung surfactant homeostasis.
Cobalt induces heme oxygenase-1 expression by a hypoxia-inducible factor-independent mechanism in Chinese hamster ovary cells: regulation by Nrf2 and MafG transcription factors.
This work was supported in part by Grants-in-Aid from the Japan Society for the Promotion of Science 19K07680 and 16K07108 (to M. M.) and 15H02506, 24390066, 21249014, and 18H04021 (to K. I.) and Agency for Medical Research and Development Grant JP16gm050001 (to K. I.). H. Nishizawa received DFX as raw material from Novartis Pharma for this study. The authors declare that they have no conflicts of interest with the contents of this article.
49550-0&id=gr1.jpg){kind=link}
49550-0&id=gr2.jpg){kind=link}
49550-0&id=gr3.jpg){kind=link}
49550-0&id=gr4.jpg){kind=link}
49550-0&id=gr5.jpg){kind=link}
49550-0&id=gr6.jpg){kind=link}