Purpose: Cell-to-cell variability in apoptosis signaling contributes to heterogenic responses to cytotoxic stress in clinically heterogeneous neoplasia, such as acute myeloid leukemia (AML). The BCL-2 proteins BAX and BAK can commit mammalian cells to apoptosis and are inhibited by retrotranslocation from the mitochondria into the cytosol. The subcellular localization of BAX and BAK could determine the cellular predisposition to apoptotic death.

Experimental Design: The relative localization of BAX and BAK was determined by fractionation of AML cell lines and patient samples of a test cohort and a validation cohort.

Results: This study shows that relative BAX localization determines the predisposition of different AML cell lines to apoptosis. Human AML displays a surprising variety of relative BAX localizations. In a test cohort of 48 patients with AML, mitochondria-shifted BAX correlated with improved patient survival, FLT3-ITD status, and leukocytosis. Analysis of a validation cohort of 80 elderly patients treated with myelosuppressive chemotherapy confirmed that relative BAX localization correlates with probability of disease progression, FLT3-ITD status, and leukocytosis. Relative BAX localization could therefore be helpful to identify elderly or frail patients who may benefit from cytotoxic therapy.

Conclusions: In this retrospective analysis of two independent AML cohorts, our data suggest that Bax localization may predict prognosis of patients with AML and cellular predisposition to apoptosis, combining the actual contribution of known and unknown factors to a final “common path.” Clin Cancer Res; 23(16); 4805–16. ©2017 AACR.

Translational Relevance

Differential sensitivity of malignant cells to chemotoxic stress based on cell-to-cell variability in cell death signaling has become a central challenge in oncology. The proapoptotic BCL-2 proteins BAX and BAK can commit human cells to apoptosis and are regulated by dynamic shuttling between cytosol and mitochondria. The analysis of the relative BAX/BAK localization in human cells and pretreatment patient samples from two cohorts of patients with AML predict cellular disposition to apoptosis and treatment response. Our data suggest relative BAX/BAK localization as potential tumor-based biomarker to predict therapeutic response to cytotoxic stimuli in AML and perhaps other tumor entities.

Acute myeloid leukemia (AML) is a disease characterized by clinical and biological heterogeneity. Recently, genomics have facilitated insights into the genetic landscape of AML and identified novel mutations with functional implications (1–3). Genetic alterations have been described to be of major prognostic relevance (4–6). Genetic diversity, however, is not the sole determinant of relapse, drug resistance, and aggressiveness of leukemia biology (7). Large-scale genetic analysis approaches naturally lack insights on posttranscriptional protein expression, protein regulation, and localization, that is, the actual level of activation. Cell-to-cell variability in protein levels and in some cases subcellular localization induce large differences in the responsiveness to uniform stimuli between individual cells and descendants from different clones (8–10). Within a clonal cell population, some cells rapidly induce programmed cell death in response to death receptor–induced apoptosis, while other cells appear more resistant and survive (11). This “intratumoral” heterogeneity and the underlying differences in cell signaling critically influence cellular predisposition to apoptosis upon cytotoxic stress (e.g., induced by chemotherapy). Differential sensitivity of malignant cells to chemotoxic stress has become a central challenge in oncology. Therefore, the identification of prognostic factors and the analysis of individual cellular predisposition to apoptosis is a major focus of translational cancer research. Investigation of the apoptotic machinery may offer novel strategies to sensitize resistant cells to cytotoxic agents.

Analysis of signaling between prosurvival BCL-2 proteins, the proapoptotic BCL-2 proteins BAX and BAK, and BH3-only proteins (sharing only the BH3 domain with BCL-2) on the outer mitochondrial membrane has been proposed to predict cellular predisposition to chemotoxic cell death. Cellular cytotoxic response is suggested to depend on the presence of prosurvival BCL-2 proteins and the accumulated BH3-only proteins on the mitochondria, termed “mitochondrial priming,” resulting perhaps from prior sublethal cell stresses (12–15). Prosurvival BCL-2 proteins inhibit BAX and BAK via direct interactions or by sequestering “activator” BH3-only proteins, thereby preventing their interaction with BAX and BAK (16–21). In response to stress, BH3-only protein signaling is thought to activate BAX and BAK, which in turn permeabilize the outer mitochondrial membrane and release of cytochrome c (cyt c). Cytosolic cyt c results in mitochondrial dysfunction and initiates the caspase cascade that efficiently dismantles the cell (22, 23). The activation of BAX or BAK commits the cell to apoptosis necessitating a tight control of these proapoptotic BCL-2 proteins (24). In healthy cells, BAX and BAK are controlled by constant retrotranslocation of mitochondria-associated protein into the cytosol (25, 26). Permanent translocation of BAX and BAK to the mitochondria establishes an equilibrium between cytosolic and mitochondrial protein pools (25, 27). These processes are dependent on the porin VDAC2, which serves as a platform for BAX retrotranslocation (28). Differential rates of BAX and BAK shuttling determine localizations of BAX (predominantly in the cytosol) and BAK (largely on the mitochondria; ref. 26). Retrotranslocation of BAX and BAK depends on interactions with prosurvival BCL-2 proteins (25, 26). Importantly, BAX retrotranslocation rates determine the size of the mitochondrial BAX pool and cellular response to apoptosis stimulation (29). Moreover, BH3-only proteins reduce the rate of BAX retrotranslocation (25), suggesting a reflection of the BCL-2 protein signaling on the mitochondria (and “mitochondrial priming”) and perhaps other factors in the equilibrium between cytosolic and mitochondrial BAX.

Here, we provide first evidence for the correlation of the subcellular BAX (and BAK) localization with susceptibility to cytotoxic therapy in human AML cell lines and primary AML blasts. While BAX shift toward the mitochondria correlates with increased cellular sensitivity to chemotoxic stress and beneficial prognostic markers, BAX shift toward the cytoplasm rather correlates with cellular resistance to chemotherapeutic stress and adverse prognostic markers. Thus, relative BAX localization has the potential to combine all relevant factors to assess the cellular susceptibility to apoptosis induction as one of the “final effectors” of cell death.

Cell culture

HeLa cells were cultured in DMEM medium supplemented with 10 mmol/L HEPES and 10% heat-inactivated FBS in 5% CO2 at 37°C. All AML cell lines and primary AML cells were cultured in RPMI1640 medium supplemented with 10 mmol/L Hepes, 10% heat-inactivated FBS, 4.5 g/L glucose, and 1 mmol/L sodium pyruvate in 5% CO2 at 37°C. Cells were regularly (in 4-week intervals) tested for potential mycoplasma infection using the Venor GeM kit (Biochrom).

Whole-cell lysis and subcellular fractionation

For the standard curve whole-cell lysates, HeLa cell lysate was used to establish a standardized mix of all analyzed proteins, ensuring similar standardization of all patient samples. HeLa cells were harvested, washed with PBS, and subsequently resuspended in cell lysis buffer (20 mmol/L Tris, 100 mmol/L NaCl, 1 mmol/L EDTA, 0.5% Triton X-100, pH 7.5) supplemented with protease inhibitors. Upon incubation on ice for 15 minutes, the samples were centrifuged at 15,000 × g for 10 minutes at 4°C. The supernatants were subjected to acetone precipitation, followed by resuspension in SDS sample buffer and storage at −80°C. To obtain mitochondrial and cytosolic fractions of patient bone marrow samples, mononuclear cells were isolated by Ficoll-Paque density gradient centrifugation. Cells were washed with PBS, resuspended in RPMI supplemented with 10% DMSO and 20% FCS, and stored in liquid nitrogen. Aspirates were then thawed and centrifuged for 5 minutes at 1,500 × g at 4°C. The mean percentage of blasts in the bone marrow aspirates (n = 44) was 64% (95% CI: 54%–74%). Upon washing with PBS, the cell pellet was resuspended in SEM buffer (10 mmol/L HEPES, 250 mmol/L sucrose, pH 7.2) supplemented with protease inhibitors and homogenized using the MINILYS (PeqLab) system. Subsequently, samples were centrifuged at 1,500 × g for 5 minutes at 4°C. The supernatant was transferred to a new tube and centrifuged for 30 minutes at 13,000 × g at 4°C. While sedimented mitochondria were washed two times, the supernatant of this step, the cytosolic fraction, was ultracentrifuged at 150,000 × g for 1 hour at 4°C. Finally, the samples were separated by SDS-PAGE and subjected to Western blot analysis.

Quantification of Western blot data

Band intensities were captured using the Fujifilm LAS-4000 imager and quantified using ImageJ. On the basis of a whole-cell lysate standard curve within each individual blot, BAX (E63, Abcam) and BAK (Millipore) were then quantified using SigmaPlot. The amounts of mitochondrial and cytosolic protein were then determined as relative to COX IV (Invitrogen) and β-actin (Millipore), respectively. Finally, ratios were built from the mitochondrial and cytosolic values to obtain the relative protein localization. Accumulation of the cleaved form of the caspase substrate PARP was analyzed by Western blot (pAB, Cell Signaling Technology) to determine apoptosis progression.

Patients

Forty-eight bone marrow samples of the test cohort were obtained from the Hematology Tumor Bank Magdeburg (HTM) at the Department of Hematology and Oncology, Medical Center, Otto-von-Guericke University (Magdeburg, Germany). These samples were obtained after written informed consent of patients and in accordance with the Declaration of Helsinki. This study was approved by the institutional review board of the University Hospital Magdeburg and the local ethics committee (file #15/2008). The majority of patients were diagnosed with AML (n = 45; 93.8%). Other myeloid diagnoses included blast phase of chronic myeloid leukemia (n = 1), blastic plasmacytoid dendritic cell neoplasm (BPDCN; n = 1), and one patient was diagnosed with acute biphenotypic leukemia. Three of all samples (6.3%) were from relapsed patients. Analyzing the AML samples (n = 45) only, three were taken in relapse (6.7%), and 74.4% (32/43) were de novo AML. Of these AML samples, 20.9% (9/43) were FLT3-ITD (FMS-like tyrosine kinase 3 internal tandem duplication) positive, and 18.6% (8/43) were at good, 51.2% (22/43) were at intermediate, and 30.2% (13/43) were at poor risk according ELN criteria (6), respectively.

SNV analyses

DNA from 40 available samples (low n = 8; medL n = 9; medH n = 11; high n = 12) were isolated and genotyped using the Illumina OmniExpressExome 1.2 Chip. The data were analyzed using Illumina Genome Studio 2011.1 Genotyping Module to obtain normalized genotype data. CNV analysis was carried out using the CNVpartition 3.2.0 algorithm within Genome Studio using default parameters. Somatic copy-number alterations were determined using GISTIC 2.0. Identification of differential single nucleotide variants (SNV) was performed using PLINK (P < 0.001).

Statistical analysis

For the main statistical analysis, Sigma Plot was used. The logistic regression and ROC calculation was done using SAS Software 9.4. A significance level of α = 5% was used for testing; however, all P values have explorative character only. For the ROS, curve area under curve and 95% confidence intervals were calculated using the Mann–Whitney method.

Caspase-3/7 assay

For caspase assays, cells were washed with ice-cold 1× PBS and resuspended in ice-cold Cell Lysis Buffer (20 mmol/L TRIS pH 7.4, 100 mmol/L NaCl, 1 mmol/L EDTA, including protease inhibitor cocktail and 0.5 % Triton X-100). Whole-cell lysate was incubated with Caspase-3/7 substrate (BD Pharmingen) for 60 minutes at 37°C, and protein concentration was determined by a Bradford Assay (Roth). Substrate cleavage was measured for 50 cycles with 10-second delay (excitation at 380 nm, emission at 430–460 nm). Kinetics were determined and calculated to the amount of protein per sample.

Lactate dehydrogenase assay

Cells were seeded in 96-well plate format with mock-treated samples and lysis control on the same plate. Cells were treated with 1 or 20 μmol/L daunorubicin, 1 μmol/L actinomycin D, 100 μmol/L etoposide, 5 μmol/L doxorubicin, 10 μmol/L ABT-737, 10 μmol/L UMI-77, or 1 μmol/L staurosporin for 24 hours. Then, 50 μL supernatant was transferred to fresh a 96-well plate and mixed with substrate (Promega). Lactate dehydrogenase (LDH) activity was monitored over 25 cycles with a 25-second delay, and kinetics were calculated.

High mitochondrial BAX levels correlate with improved AML patient survival

The regulation of BAX and BAK by steady translocation to the mitochondria and retrotranslocation back into the cytosol of nonapoptotic cells raises the question whether BAX/BAK localization could predict cellular predisposition to apoptosis. Therefore, the relative distribution of BAX and BAK between cytosolic and mitochondrial pools in primary bone marrow samples of human AML was analyzed in a test cohort by fractionating 48 patient samples into cytosol and mitochondria (Fig. 1A). To allow comparison of samples not analyzed on the same Western blot, relative mitochondrial and relative cytosolic protein was determined using the fractionation loading controls COX IV and β-actin, respectively. All detected proteins were normalized using a titration of the same HeLa cell extract and thus the same amounts of all detected proteins for the evaluation of every analyzed sample. Both relative protein pools were then combined to calculate the relative protein localization (relative BAX localization = mitochondrial BAX/cytosolic BAX), with high relative protein localization values pointing to increased mitochondrial protein and low values to predominantly cytosolic protein independent of cellular protein levels (Fig. 1B). Surprisingly, the analysis of patient bone marrow samples revealed a high interindividual variability of relative BAX and BAK localizations (Fig. 1C; Supplementary Table S1; Supplementary Fig. S1). Nonetheless, relative BAX localizations show a general trend towards cytosolic localizations, whereas BAK is present in larger mitochondrial pools (Supplementary Fig. S2A). Comparison of cytosolic and mitochondrial pools of BAX or BAK shows no apparent correlation (Fig. 1D; Supplementary Fig. S2B). The absence of correlations between the levels of mitochondrial and cytosolic BAX suggests that BAX localization does not depend on cellular BAX levels. Also the maintenance of a constant ratio between mitochondrial and cytosolic BCL-2 protein pools is not apparent.

The analysis of test cohort quartiles [low, mostly cytosolic; high, largely mitochondrial, medium high (medH) and medium (medL)] of relative BAX localization reveals a strikingly better survival probability among patients with high relative BAX localization (Fig. 1E). Low relative BAX localization values increase the odds of death during treatment by 11.5 times. The same difference is apparent analyzing only patients from the test cohort with de novo AML (Supplementary Fig. S3). On the basis of the similar response of patients with intermediate or low relative BAX localization, the use of relative BAX localization as continuous value is not appropriate to distinguish responder and nonresponder (Supplementary Fig. S4). However, categorizing patients with AML according to their relative BAX localization could predict treatment response.

Cytosolic BAX localization is associated with decreased predisposition to apoptosis and FLT3-ITD

Among test cohort patients, relative BAX localization shows a strong positive correlation with relative BAK localization (Fig. 2A). A strong correlation between cytosolic BAX and relative BAX localization suggests high variability of cytosolic BAX pools compared with similar mitochondrial levels detected among the patients with AML (Supplementary Fig. S5A and S5B). These results indicate that cells tolerate a larger variety of BAX concentrations in the cytosol, probably explained by the fact that cytosolic BAX is monomeric and lacks apparent protein interactions (30). Of note, the strong positive correlation between cytosolic BAX and BAK suggest that cells also tolerate a large variety of cytosolic BAK despite a greater tendency of this BCL-2 protein to localize on the mitochondria (Supplementary Fig. S5C). The subcellular distribution of BAX and BAK is clearly regulated by the same retrotranslocation process in pretreatment bone marrow samples of patients with AML. Therefore, the localization analysis of BAX is sufficient to assess the cellular BAX/BAK regulation.

The potential of relative BAX localization as a diagnostic tool was analyzed in the test cohort of 48 patients with AML (Supplementary Table S2.) and correlated with clinical data. Exemplarily, FMS-like tyrosine kinase 3 internal tandem duplications (FLT3-ITD), a bona fide mutational marker of negative prognostic impact in AML (31, 32), were investigated. Strikingly, FLT3-ITD–positive patients show low relative BAX localization values, suggesting increased cytosolic BAX and increased protection against apoptosis (Fig. 2B). Accordingly, FLT3-ITD is associated with a nonfavorable prognosis (Supplementary Fig. S6A). FLT3-ITD has been associated with activation of both prosurvival BCL-2 proteins and BH3-only proteins and could therefore induce changes in the relative BAX localization (33–36). Analysis of patient collective quartiles [low, mostly cytosolic; high, largely mitochondrial, medium high (medH) and medium (medL)] of relative BAX localization confirms predominant occurrence of FLT3-ITD in samples with low relative BAX localization values (Fig. 2C). The relative risk of FLT3-ITD in this group is at least 6.25 times of that in the groups with medium and high relative BAX localization values. The same coherence is apparent based on the analysis of the relative BAK localization, again suggesting a BAX localization–based analysis is sufficient for the regulatory status of both prosurvival BCL-2 proteins (Supplementary Fig. S6B). On the other hand, the analysis based on BAX levels instead of BAX localization shows FLT3-ITD occurrence in patients with high cellular BAX levels (Supplementary Fig. S7). Remarkably, patients with intermediate BAX localization values show inferior survival compared with high relative BAX localization patients despite a similar chance of FLT3-ITD (Figs. 1E and 2B). Of note, one FLT3-ITD–positive individual with high relative BAX localization value shows prolonged survival following chemotherapy. FLT3-ITD thus may influence relative BAX localization, but many other factors (e.g., genetic or epigenetic) also determine BAX localization and the predisposition to apoptosis.

Consistently, single-nucleotide variants (SNVs) occurring in patient samples of high or other relative BAX localization were analyzed. Overall, distribution and frequency of SNVs was similar across the different BAX localizations reflecting a quite homogenous pattern of genetic alterations and very little genomic diversity among the different patient groups (Fig. 2D). Notably, patient samples with high relative BAX localization displayed a significantly lower number of genetic losses compared with all other samples. GISTIC analyses further confirmed a high concordance of driving copy-number alterations across the different groups and very few specific alterations (Supplementary Table S3). Notably, while no recurrent losses were detected in high BAX localization, both intermediate and low relative BAX localization quartiles showed losses on 5q31.3. Nevertheless, comparison of the differences between SNVs in BAX high versus the other subgroups showed little overlap among the different subgroups, indicating the different BAX localizations might be determined by diverse molecular mechanisms and cell stress rather than common genetic alterations (Supplementary Fig. S8).

In an attempt to identify further factors that are reflected by the relative BAX localization, we analyzed potential and known predictors of response to therapy and survival in AML. In the test cohort, patient age seems to have no influence on relative BAX localization (Supplementary Fig. S9A). Surprisingly, the established overall risk rating shows no striking congruence with relative BAX localization (Supplementary Fig. S9B). However, patients with therapy-associated AML present with elevated mitochondrial BAX levels (Supplementary Fig. S10). This association is conceivable, as prior nontoxic cellular stress could be reflected in the relative BAX localization (29). Strikingly, the initial white blood count is significantly increased in patients with low relative BAX localization compared with patients with BAX shifted to their mitochondria (Fig. 2E). Relative BAX localization shows a negative correlation with initial leukocytosis (Fig. 2F). A shift of the BAX population towards the mitochondria and high predisposition to commitment to mitochondrial apoptosis thus is associated with lower initial leukocyte counts.

Treatment of primary AML samples of seven patients with six different apoptotic stimuli shows reduced apoptotic response by AML cells with low relative BAX localization (Fig. 2G). Patient samples with high relative Bax localization show increased apoptotic response, resulting in a correlation of relative BAX localization and apoptotic response, particularly for doxorubicin- (and daunorubicin-) treated cells (Fig. 2H). Relative Bax localization can identify cellular predisposition to apoptosis.

Relative BAX localization determines predisposition of AML cell lines to apoptosis

The link between BAX localization and predisposition to apoptosis in AML was further tested by determining the relative BAX localization in eight different AML cell lines (Fig. 3A; Supplementary Fig. S11). Consistent with the AML patient data, the analysis shows increased pools of mitochondrial BAK compared with BAX (Supplementary Fig. S12A). However, in AML cell lines, this difference is much more pronounced, while relative cytosolic BAX and relative mitochondrial BAX also lack an apparent correlation (Supplementary Fig. S12B).

The effects of relative BAX localization in AML cell lines on the induction of cell death were investigated in response to cytotoxic stress induced by daunorubicin. Application of 1 μmol/L daunorubicin for 24 hours resulted 2%–16% cell death depending on the cell line. The signal-to-noise ratio was increased using 20 μmol/L daunorubicin (37), resulting in robust cell death induction in seven of the eight tested AML cell lines (Fig. 3B). However, HEL cells did not induce cell death even at 20 μmol/L daunorubicin, revealing remarkably reduced predisposition to cell death compared with other cell lines. Noteworthy, HEL cells show a particular low relative BAX localization (Fig. 3C). Treatment with cytarabine (ara-C) for 24 hours did not induce significant apoptosis in HEL or OCI-AML3 cells (Supplementary Fig. S13). Comparing the apoptotic response of HEL cells and OCI-AML3 cells with increased mitochondrial BAX pools in response to daunorubicin treatment demonstrates a robust activation of caspases-3/7 in OCI-AML3 cells, whereas HEL cells show low levels of active caspase after 24 hours (Fig. 3D). This difference is corroborated by accumulation of the cleaved form of the caspase substrate PARP in OCI-AML3 cells but not HEL cells (Fig. 3E). These results corroborate the correlation between BAX localization and apoptotic response, resulting from altered BAX retrotranslocation rates, as previously shown in HCT116 and MEF cells (29). In parallel, HEL cells show increased resistance against different apoptotic stimuli and the BH3 mimetic UMI-77 (Fig. 3F; Supplementary Fig. S14). Among the 8 AML cell lines, MV4-11 cells contain the FLT3-ITD. FLT3 inhibition results in a shift of BAX towards the mitochondria suggesting an increase of apoptosis predisposition of MV4-11 cells but not THP-1 cells (Fig. 3G). However, MV4-11 cells are not particularly protected against apoptosis induction, suggesting an influence on relative BAX localization by FLT3 activity among other factors. The responses of AML cell lines to apoptotic stimuli suggest that relative BAX localization could predict predisposition to apoptosis.

Relative BAX localization is associated with FLT3-ITD in AML of the elderly

Next, we focused on the patients with an age of 60 years and older who are frequently not eligible for intensive chemotherapy or hematopoietic stem cell transplantation. Association of the susceptibility to cytotoxic therapy with relative BAX localization could facilitate treatment decisions toward intensive treatment versus nonintensive approaches such as targeted- or immune-therapies and best supportive care. Among elderly patients (n = 25), a pronounced difference in relative BAX localization became evident when comparing FLT3-ITD–positive to FLT3-ITD–negative patients (Fig. 4A). While 60% of elderly patients with low relative BAX localization levels harbor FLT3-ITD mutations, no elderly patient with high relative BAX localization values is FLT3-ITD–positive (Fig. 4B). These observations are supported by the analysis based on relative BAK localization (Supplementary Fig. S15). In addition, different blast counts between samples with low and high relative BAX localization values were detected (Fig. 4C). Therefore, it is tempting to speculate on a potential use of relative BAX localization to guide clinical decisions as diagnostic or prognostic tool for elderly patients with AML.

BAX and BAK localization is highly variable and correlates with clinical markers in human AML

To evaluate our previous insights into the localization of BAX and BAK in human AML and correlation of BAX localization with clinical markers, we analyzed a validation cohort of 80 elderly patients with AML treated with myelosuppressive chemotherapy within the AMLSG. Consistently with the test cohort, correlation between mitochondrial and cytosolic protein levels is not apparent (Fig. 5A; Supplementary Table S4; Supplementary Fig. S16). As observed in the test cohort, a wide variety of BAX/BAK localizations were detected with a general trend toward cytosolic BAX localization and mitochondrial BAK (Fig. 5B). A strong correlation also exists between the relative localizations of BAX and BAK in the validation cohort (Fig. 5C), demonstrating similar regulation of both proapoptotic BCL-2 proteins in human AML.

Different relative BAX/BAK localizations were detectable in the validation cohort depending on the mutational status of FLT3, although the difference was less pronounced compared with the test cohort (Fig. 5D). The same effect could be observed for the subgroup of intermediate risk patients according to ELN risk criteria (Supplementary Fig. S17). When applying the test cohort criteria and categorizing patients according to low (mostly cytosolic), medium low (medL) medium high (medH), and high (largely mitochondrial) BAX localization, decreased risk of patients with high mitochondrial BAX for FLT3-ITD positivity can be confirmed (Fig. 5E). Relative BAX localization also correlates with leukocytosis in the validation cohort (Fig. 5F; Supplementary Fig. S18), suggesting a similar correlation between relative BAX localization and predisposition to apoptosis in both patient cohorts.

Relative BAX localization correlates with therapeutic response in AML of the elderly

Next, we analyzed the potential link between relative BAX localization and therapeutic response. Of note, 16 patients from the validation cohort that underwent allogeneic stem cell transplantation were excluded, as their therapeutic response is rather related to immunogenic control than to susceptibility to cytotoxic therapy. The analysis shows a tendency of mitochondrial BAX localization towards positive response to chemotherapy (Supplementary Fig. S19). Moreover, patients with high relative BAX localization are less likely refractory to chemotherapy (Fig. 6A). This group of patients also shows a trend toward improved survival (Fig. 6B). However, elderly patients from the test cohort with high relative BAX localization values show strikingly improved survival (Fig. 6C). The effects observed in the intensively treated validation cohort (Fig. 6B) are less pronounced compared with the test cohort, where patients received nonmyelosuppressive regimens (cytoreductive therapies or best supportive care; Figs. 1E and 6C). Again, the analysis of relative BAX localization as a continuous value provides no significant difference in receiver operating curve for validation cohort patients (Supplementary Fig. S20). The analysis basing on categorized BAX levels is showing a clear trend in both patient cohorts but results are still suffering from the small number of participants. Logistic regression on survival as a function of relative BAX localization, sex, FLT3-ITD, ELN-risk, and total-risk showed a significant influence of increased mitochondrial BAX levels OR = 1.61 (95% confidence interval, 1.1–2.4), while all other showed trends on survival influence but failed to reach significance (Supplementary Table S5).

BAX and BAK share common regulation of their localization in primary AML patient samples and AML cell lines. In healthy cells, both proapoptotic BCL-2 proteins are retrotranslocated from the mitochondria into the cytosol by the same prosurvival BCL-2 protein-dependent mechanism (25, 26). Shuttling at similar rates, BAX and BAK share similar apoptotic activity and subcellular localization (26). Therefore, the analysis of the relative BAX localization is sufficient to determine the status of BAX/BAK signaling. The total cellular protein level is not linked to BAX regulation. Increased BAX level could even be associated with more apoptosis-resistant scenarios (Supplementary Fig. S6). This argues against previous attempts to predict cellular survival or chemotherapy response solely by analyzing cellular BAX levels (38, 39). Also, the determination of the BAX/BCL-2 ratio, as previously suggested (40, 41), considers a potentially important interaction for BAX retrotranslocation, but neglects BAK-dependent apoptosis and potential cell protection from BAX/BAK activity by other prosurvival proteins, such as Bcl-xL and Mcl-1 (25). The analysis of relative BAX/BAK localization includes the combined activities of prosurvival BCL-2 proteins and potential influence of BH3-proteins, like Bim, on BAX/BAK retrotranslocation (25, 26, 29). The equilibrium between cytosolic and mitochondrial BAX pools displays a surprising diversity in AML patient samples that reflects the heterogeneous response to chemotoxic stress. Significant differences mainly in the cytosolic BAX pool hold the promise of discriminating between clinically different AML.

Cellular survival requires fast BAX retrotranslocation from the mitochondria (26, 27). Consequently, a shift of the equilibrium between cytosolic and mitochondrial BAX toward the mitochondria results in increased apoptosis induction in response to apoptotic stimuli (29). In parallel to other cell cultures, reduced predisposition to apoptosis is observed in HEL cells with cytosol-shifted BAX compared with other AML cell lines. The other tested cell lines show a similar response to apoptotic stimuli, corroborating previous results from BH3 profiling (42). While BH3 profiling has the potential to dissect contributions of different key proteins, cell permeabilization, required for this type of analysis, removes cytosolic protein pools and triggers reequilibration potentially interfering with the determination of cellular predisposition to apoptosis. Relative BAX localization reveals the status of cellular BAX/BAK regulation, including contributions of known and unknown factors but cannot dissect their contributions.

Relative BAX localization itself is influenced by a plethora of factors, probably including genetic aberrations, for example, FLT3-ITD. However, genetic determinants of BAX regulation are heterogeneous in the analyzed patients' samples. The use of genomic data for risk stratification is probably prone to underscore the effect of mutational events on transcriptional, translational, and posttranslational consequences. They may be further influenced by the microenvironment or epigenetic alterations. Presence of complex aberrant cytogenetic changes has been described to correlate with resistance to low-dose cytarabine treatment before (43). However, even patients without any cytogenetic aberrations do not uniformly respond to chemotherapy. Most recently, an integrative genomic and clinical analysis of 1,540 AML patients enrolled in three prospective clinical trials and treated with myelosuppressive chemotherapy identified a set of 5,234 driver mutations with more than one driver mutation occurring in the vast majority of patients (86%; ref. 44). Moreover, cooccurring mutations influenced outcome and survival also beyond the previously described categories (6), highlighting the complexity of risk-assessment in AML. Here, determination of cellular susceptibility toward induction of apoptosis could potentially serve as a rapid and suitable assessment to predict therapeutic response and to stratify patients towards myelosuppressive chemotherapy versus rather nonintensive therapy such as targeted agents, epigenetic treatment, or immunotherapy.

Clinically, the predictive value of certain prognostic markers has a distinct impact, depending on the patient's fitness and comorbidities. AML of the elderly or frail patient requires detailed assessment of physical status and past medical history before patients are stratified toward intensive treatment versus nonintensive therapies or best supportive care. The majority of patients diagnosed with AML is above the age of 60 years and many of them do not have any curative options due to comorbidities or lack of response. Achievement of remission leads to prolonged survival and improved quality of life of the patient; however, which patients will respond is difficult to predict. When focusing on these patient-oriented aspects, prediction of response by using one major regulator of cellular resistance could facilitate stratification of elderly and frail patients to chemotherapeutic approaches versus targeted therapies or best supportive care.

In principle, BAX localization can be analyzed within a short period of time (less than 72 hours) and requires low amounts of cells (<1 × 107 cells) derived from the peripheral blood or bone marrow. Our analysis has shown that the use of relative BAX localization as categorical but not continuous value identifies patients that are likely to have favorable therapeutic response. The strength of the analysis of the relative BAX localization in comparison with other available markers is the direct detection of the cellular predisposition to apoptosis in response to cytotoxic stress or specific inhibitors (e.g., BH3 mimetics). This final stage of cell fate decision could incorporate various effectors of leukemia development and maintenance: genomic and epigenetic aberrations, metabolic changes, cell-intrinsic and cell-extrinsic signaling, and modifications/homeostasis of the proteome.

In our analysis, determination of relative BAX localization shows that BAX shift toward the mitochondria is associated with predisposition to apoptosis and improved patient survival. These findings indicate increased susceptibility of the predominant AML clone to cytotoxic therapies. BAX localization shift toward the cytosol, on the other hand, leads to a resistance phenotype and is associated with a high tumor burden and FLT3-ITD, described to coincide with poor prognostic outcome. These correlations are apparent, although not homogeneous, in both cohorts of patients with AML analyzed in this article. Of note, small sample size and heterogeneity of consolidation (e.g., allogeneic stem cell transplantation) or second-line therapies (demethylating agents, targeted therapies) may significantly influence analysis of survival as readout for response. Still, BAX localization may predict response and outcome and its predictive value is even better detectable in analysis of rather small but more homogeneous subgroups (Fig. 6B and C). Prospective analysis of response rates to first-line treatment especially in a higher number of older patients with AML is necessary to validate its predictive value in the future.

Our findings indicate the potential use of BAX localization to assess for response to cytotoxic therapy in AML but also in the treatment of other tumor entities, as the underlying BAX regulation is a general mechanism that prevents BAX activity in human cells. Nevertheless, this analysis identifies also the limits of our analysis, as we assess biologic response of the dominant AML clone at diagnosis. Subclones may be selected or emerge during chemotherapy. Overall, we provide first evidence that BAX localization is associated with response to chemotherapy.

No potential conflicts of interest were disclosed.

Conception and design: F. Reichenbach, C. Wiedenmann, F. Edlich

Development of methodology: F. Reichenbach, C. Wiedenmann, F. Edlich

Acquisition of data (provided animals, acquired and managed patients, provided facilities, etc.): F. Reichenbach, C. Wiedenmann, E. Schalk, D. Wolleschak, K. Döhner, J.U. Marquardt, F. Heidel

Analysis and interpretation of data (e.g., statistical analysis, biostatistics, computational analysis): F. Reichenbach, C. Wiedenmann, K. Funk, P. Scholz-Kreisel, F. Todt, K. Döhner, J.U. Marquardt, F. Heidel, F. Edlich, D. Becker

Writing, review, and/or revision of the manuscript: F. Reichenbach, C. Wiedenmann, E. Schalk, P. Scholz-Kreisel, K. Döhner, J.U. Marquardt, F. Heidel, F. Edlich

Administrative, technical, or material support (i.e., reporting or organizing data, constructing databases): F. Reichenbach, F. Heidel

Study supervision: F. Edlich

This work was supported by the DFG Emmy Noether program, the Collaborative Research Cluster (CRC) 746, the Else Kröner-Fresenius-Stiftung, the Wilhelm Sander-Stiftung, the Spemann Graduate School of Biology and Medicine (SGBM, GSC-4) and the Centre for Biological Signalling Studies (BIOSS, EXC-294) funded by the Excellence Initiative of the German Federal and State Governments (F. Reichenbach, C. Wiedenmann, F. Todt, F. Edlich). F. Heidel was supported by a grant of the German Research Foundation, Collaborative Research Cluster (CRC) 854 (Project A20), the ProExcellence Research Initiative “RegenerAging” (State of Thuringia), the Thuringian country programme ProExzellenz (RegenerAging - FSU-I-03/14) of the Thuringian Ministry for Research (TMWWDG), and the German Jose-Carreras Leukemia Society (DJCLS SP12/08). K. Döhner was supported by a grant of the German Research Foundation, Collaborative Research Cluster (CRC) 1074 (Project B3) and the AMLSG study group. J.U. Marquardt is supported by the German Cancer Aid (DKH 110989), the Wilhelm Sander-Stiftung, and the Volkswagen Foundation (Lichtenberg program).

The costs of publication of this article were defrayed in part by the payment of page charges. This article must therefore be hereby marked advertisement in accordance with 18 U.S.C. Section 1734 solely to indicate this fact.

1.
The Cancer Genome Atlas Research Network
. 
Genomic and epigenomic landscapes of adult de novo acute myeloid leukemia
.
N Engl J Med
2013
;
368
:
2059
74
.
2.
Ley
TJ
,
Ding
L
,
Walter
MJ
,
McLellan
MD
,
Lamprecht
T
,
Larson
DE
, et al
DNMT3A mutations in acute myeloid leukemia
.
N Engl J Med
2010
;
363
:
2424
33
.
3.
Mardis
ER
,
Ding
L
,
Dooling
DJ
,
Larson
DE
,
McLellan
MD
,
Chen
K
, et al
Recurring mutations found by sequencing an acute myeloid leukemia genome
.
N Engl J Med
2009
;
361
:
1058
66
.
4.
Bullinger
L
,
Döhner
K
,
Bair
E
,
Fröhling
S
,
Schlenk
RF
,
Tibshirani
R
, et al
Use of gene-expression profiling to identify prognostic subclasses in adult acute myeloid leukemia
.
N Engl J Med
2004
;
350
:
1605
16
.
5.
Byrd
JC
,
Mrózek
K
,
Dodge
RK
,
Carroll
AJ
,
Edwards
CG
,
Arthur
DC
, et al
Pretreatment cytogenetic abnormalities are predictive of induction success, cumulative incidence of relapse, and overall survival in adult patients with de novo acute myeloid leukemia: results from Cancer and Leukemia Group B (CALGB 8461)
.
Blood
2002
;
100
:
4325
36
.
6.
Döhner
H
,
Estey
EH
,
Amadori
S
,
Appelbaum
FR
,
Büchner
T
,
Burnett
AK
, et al
Diagnosis and management of acute myeloid leukemia in adults: recommendations from an international expert panel, on behalf of the European LeukemiaNet
.
Blood
2010
;
115
:
453
74
.
7.
Heidel
FH
,
Arreba-Tutusaus
P
,
Armstrong
SA
,
Fischer
T
. 
Evolutionarily conserved signaling pathways: ACTINg in the shadows of acute myelogenous leukemia's genetic diversity
.
Clin Cancer Res
2015
;
21
:
240
8
.
8.
Chang
HH
,
Hemberg
M
,
Barahona
M
,
Ingber
DE
,
Huang
S
. 
Transcriptome-wide noise controls lineage choice in mammalian progenitor cells
.
Nature
2008
;
453
:
544
7
.
9.
Feinerman
O
,
Veiga
J
,
Dorfman
JR
,
Germain
RN
,
Altan-Bonnet
G
. 
Variability and robustness in T cell activation from regulated heterogeneity in protein levels
.
Science
2008
;
321
:
1081
4
.
10.
Sigal
A
,
Milo
R
,
Cohen
A
,
Geva-Zatorsky
N
,
Klein
Y
,
Liron
Y
, et al
Variability and memory of protein levels in human cells
.
Nature
2006
;
444
:
643
6
.
11.
Spencer
SL
,
Gaudet
S
,
Albeck
JG
,
Burke
JM
,
Sorger
PK
. 
Non-genetic origins of cell-to-cell variability in TRAIL-induced apoptosis
.
Nature
2009
;
459
:
428
32
.
12.
Ryan
JA
,
Brunelle
JK
,
Letai
A
. 
Heightened mitochondrial priming is the basis for apoptotic hypersensitivity of CD4+ CD8+ thymocytes
.
Proc Natl Acad Sci U S A
2010
;
107
:
12895
900
.
13.
Certo
M
,
Del Gaizo Moore
V
,
Nishino
M
,
Wei
G
,
Korsmeyer
S
,
Armstrong
SA
, et al
Mitochondria primed by death signals determine cellular addiction to antiapoptotic BCL-2 family members
.
Cancer Cell
2006
;
9
:
351
65
.
14.
Ni Chonghaile
T
,
Sarosiek
KA
,
Vo
T-T
,
Ryan
JA
,
Tammareddi
A
,
Moore
VDG
, et al
Pretreatment mitochondrial priming correlates with clinical response to cytotoxic chemotherapy
.
Science
2011
;
334
:
1129
33
.
15.
Brunelle
JK
,
Ryan
J
,
Yecies
D
,
Opferman
JT
,
Letai
A
. 
MCL-1-dependent leukemia cells are more sensitive to chemotherapy than BCL-2-dependent counterparts
.
J Cell Biol
2009
;
187
:
429
42
.
16.
Kim
H
,
Rafiuddin-Shah
M
,
Tu
H-C
,
Jeffers
JR
,
Zambetti
GP
,
Hsieh
JJ-D
, et al
Hierarchical regulation of mitochondrion-dependent apoptosis by BCL-2 subfamilies
.
Nat Cell Biol
2006
;
8
:
1348
58
.
17.
Kuwana
T
,
Bouchier-Hayes
L
,
Chipuk
JE
,
Bonzon
C
,
Sullivan
BA
,
Green
DR
, et al
BH3 domains of BH3-only proteins differentially regulate BAX-mediated mitochondrial membrane permeabilization both directly and indirectly
.
Mol Cell
2005
;
17
:
525
35
.
18.
Letai
A
,
Bassik
MC
,
Walensky
LD
,
Sorcinelli
MD
,
Weiler
S
,
Korsmeyer
SJ
. 
Distinct BH3 domains either sensitize or activate mitochondrial apoptosis, serving as prototype cancer therapeutics
.
Cancer Cell
2002
;
2
:
183
92
.
19.
Llambi
F
,
Moldoveanu
T
,
Tait
SWG
,
Bouchier-Hayes
L
,
Temirov
J
,
McCormick
LL
, et al
A unified model of mammalian BCL-2 protein family interactions at the mitochondria
.
Mol Cell
2011
;
44
:
517
31
.
20.
Willis
SN
,
Chen
L
,
Dewson
G
,
Wei
A
,
Naik
E
,
Fletcher
JI
, et al
Proapoptotic BAK is sequestered by Mcl-1 and Bcl-xL, but not BCL-2, until displaced by BH3-only proteins
.
Genes Dev
2005
;
19
:
1294
305
.
21.
Willis
SN
,
Fletcher
JI
,
Kaufmann
T
,
van Delft
MF
,
Chen
L
,
Czabotar
PE
, et al
Apoptosis initiated when BH3 ligands engage multiple BCL-2 homologs, not BAX or BAK
.
Science
2007
;
315
:
856
9
.
22.
Bratton
SB
,
Cohen
GM
. 
Caspase cascades in chemically-induced apoptosis
.
Adv Exp Med Biol
2001
;
500
:
407
20
.
23.
Green
DR
,
Kroemer
G
. 
The pathophysiology of mitochondrial cell death
.
Science
2004
;
305
:
626
9
.
24.
Lindsten
T
,
Ross
AJ
,
King
A
,
Zong
WX
,
Rathmell
JC
,
Shiels
HA
, et al
The combined functions of proapoptotic BCL-2 family members BAK and BAX are essential for normal development of multiple tissues
.
Mol Cell
2000
;
6
:
1389
99
.
25.
Edlich
F
,
Banerjee
S
,
Suzuki
M
,
Cleland
MM
,
Arnoult
D
,
Wang
C
, et al
Bcl-x(L) retrotranslocates BAX from the mitochondria into the cytosol
.
Cell
2011
;
145
:
104
16
.
26.
Todt
F
,
Cakir
Z
,
Reichenbach
F
,
Emschermann
F
,
Lauterwasser
J
,
Kaiser
A
, et al
Differential retrotranslocation of mitochondrial BAX and BAK
.
EMBO J
2015
;
34
:
67
80
.
27.
Edlich
F
. 
The great migration of BAX and BAK
.
Mol Cell Oncol
2015
;
2
:
e995029
.
28.
Lauterwasser
J
,
Todt
F
,
Zerbes
RM
,
Nguyen
TN
,
Craigen
W
,
Lazarou
M
, et al
The porin VDAC2 is the mitochondrial platform for BAX retrotranslocation
.
Sci Rep
2016
;
6
:
32994
.
29.
Todt
F
,
Cakir
Z
,
Reichenbach
F
,
Youle
RJ
,
Edlich
F
. 
The C-terminal helix of Bcl-x(L) mediates BAX retrotranslocation from the mitochondria
.
Cell Death Differ
2013
;
20
:
333
42
.
30.
Hsu
YT
,
Youle
RJ
. 
BAX in murine thymus is a soluble monomeric protein that displays differential detergent-induced conformations
.
J Biol Chem
1998
;
273
:
10777
83
.
31.
Fröhling
S
,
Schlenk
RF
,
Breitruck
J
,
Benner
A
,
Kreitmeier
S
,
Tobis
K
, et al
Prognostic significance of activating FLT3 mutations in younger adults (16 to 60 years) with acute myeloid leukemia and normal cytogenetics: a study of the AML Study Group Ulm
.
Blood
2002
;
100
:
4372
80
.
32.
Gale
RE
,
Hills
R
,
Pizzey
AR
,
Kottaridis
PD
,
Swirsky
D
,
Gilkes
AF
, et al
Relationship between FLT3 mutation status, biologic characteristics, and response to targeted therapy in acute promyelocytic leukemia
.
Blood
2005
;
106
:
3768
76
.
33.
Bagrintseva
K
,
Geisenhof
S
,
Kern
R
,
Eichenlaub
S
,
Reindl
C
,
Ellwart
JW
, et al
FLT3-ITD-TKD dual mutants associated with AML confer resistance to FLT3 PTK inhibitors and cytotoxic agents by overexpression of Bcl-x(L)
.
Blood
2005
;
105
:
3679
85
.
34.
Nordigården
A
,
Kraft
M
,
Eliasson
P
,
Labi
V
,
Lam
EW-F
,
Villunger
A
, et al
BH3-only protein Bim more critical than Puma in tyrosine kinase inhibitor-induced apoptosis of human leukemic cells and transduced hematopoietic progenitors carrying oncogenic FLT3
.
Blood
2009
;
113
:
2302
11
.
35.
Yang
X
,
Liu
L
,
Sternberg
D
,
Tang
L
,
Galinsky
I
,
DeAngelo
D
, et al
The FLT3 internal tandem duplication mutation prevents apoptosis in interleukin-3-deprived BaF3 cells due to protein kinase A and ribosomal S6 kinase 1-mediated BAD phosphorylation at serine 112
.
Cancer Res
2005
;
65
:
7338
47
.
36.
Yoshimoto
G
,
Miyamoto
T
,
Jabbarzadeh-Tabrizi
S
,
Iino
T
,
Rocnik
JL
,
Kikushige
Y
, et al
FLT3-ITD up-regulates MCL-1 to promote survival of stem cells in acute myeloid leukemia via FLT3-ITD-specific STAT5 activation
.
Blood
2009
;
114
:
5034
43
.
37.
McCollum
AK
,
TenEyck
CJ
,
Stensgard
B
,
Morlan
BW
,
Ballman
KV
,
Jenkins
RB
, et al
P-Glycoprotein-mediated resistance to Hsp90-directed therapy is eclipsed by the heat shock response
.
Cancer Res
2008
;
68
:
7419
27
.
38.
Gascoyne
RD
,
Krajewska
M
,
Krajewski
S
,
Connors
JM
,
Reed
JC
. 
Prognostic significance of BAX protein expression in diffuse aggressive non-Hodgkin's lymphoma
.
Blood
1997
;
90
:
3173
8
.
39.
Krajewski
S
,
Blomqvist
C
,
Franssila
K
,
Krajewska
M
,
Wasenius
VM
,
Niskanen
E
, et al
Reduced expression of proapoptotic gene BAX is associated with poor response rates to combination chemotherapy and shorter survival in women with metastatic breast adenocarcinoma
.
Cancer Res
1995
;
55
:
4471
8
.
40.
Del Poeta
G
,
Venditti
A
,
Del Principe
MI
,
Maurillo
L
,
Buccisano
F
,
Tamburini
A
, et al
Amount of spontaneous apoptosis detected by BAX/BCL-2 ratio predicts outcome in acute myeloid leukemia (AML)
.
Blood
2003
;
101
:
2125
31
.
41.
Mackey
TJ
,
Borkowski
A
,
Amin
P
,
Jacobs
SC
,
Kyprianou
N
. 
BCL-2/BAX ratio as a predictive marker for therapeutic response to radiotherapy in patients with prostate cancer
.
Urology
1998
;
52
:
1085
90
.
42.
Vo
T-T
,
Ryan
J
,
Carrasco
R
,
Neuberg
D
,
Rossi
DJ
,
Stone
RM
, et al
Relative mitochondrial priming of myeloblasts and normal HSCs determines chemotherapeutic success in AML
.
Cell
2012
;
151
:
344
55
.
43.
Dennis
M
,
Russell
N
,
Hills
RK
,
Hemmaway
C
,
Panoskaltsis
N
,
McMullin
M-F
, et al
Vosaroxin and vosaroxin plus low-dose Ara-C (LDAC) vs. low-dose Ara-C alone in older patients with acute myeloid leukemia
.
Blood
2015
;
125
:
2923
32
.
44.
Papaemmanuil
E
,
Gerstung
M
,
Bullinger
L
,
Gaidzik
VI
,
Paschka
P
,
Roberts
ND
, et al
Genomic classification and prognosis in acute myeloid leukemia
.
N Engl J Med
2016
;
374
:
2209
21
.