Multidrug resistance-1 (MDR1) acts as a chemotherapeutic drug efflux pump in tumor cells, although its physiological functions remain enigmatic. Using a recently developed MDR1-knockin reporter allele (Abcb1aAME), we found that constitutive MDR1 expression among hematopoietic cells was observed in cytolytic lymphocytes—including CD8+ cytotoxic T lymphocytes (CTLs) and natural killer cells—and regulated by Runt-related (Runx) transcription factors. Whereas MDR1 was dispensable for naive CD8+ T cell development, it was required for both the normal accumulation of effector CTLs following acute viral infection and the protective function of memory CTLs following challenge with an intracellular bacterium. MDR1 acted early after naive CD8+ T cell activation to suppress oxidative stress, enforce survival, and safeguard mitochondrial function in nascent CTLs. These data highlight an important endogenous function of MDR1 in cell-mediated immune responses and suggest that ongoing efforts to intentionally inhibit MDR1 in cancer patients could be counterproductive.

Multidrug resistance-1 (MDR1; encoded by ABCB1 in humans and Abcb1a and Abcb1b in mice) is a membrane-associated, ATP-dependent efflux pump recognized and named for removing cytostatic drugs from tumor cells. Accordingly, MDR1 has been viewed as a “dedicated drug handler” in mammalian cells (Borst and Schinkel, 2013; Zhou, 2008), and numerous MDR1 inhibitors have been developed and tested in clinical cancer trials, with little success (Bauer et al., 2005; Kelly et al., 2012; O’Brien et al., 2010; Peck et al., 2001; Seiden et al., 2002). However, MDR1 is also expressed in a number of normal cell types and tissues (Schinkel et al., 1995; Sugawara et al., 1988; Thiebaut et al., 1987), and both the presence of MDR1 orthologues in prokaryotes and a growing body of literature suggest that this transporter has conserved endogenous functions in eukaryotes that extend beyond interacting with synthetic medicines.

In the immune system, MDR1 expression has been reported in skin dendritic cells, CD4+-induced T regulatory and T effector (Teff) cells, CD8+ CTLs, and natural killer (NK) cells (Chaudhary et al., 1992; Chaudhary and Roninson, 1991; Egashira et al., 1999; Randolph et al., 1998). MDR1 has been suggested to regulate egress of skin dendritic cells into lymphatic vessels, promote induced T regulatory development, and protect IFN-γ–producing (T helper [Th]1) and IL-17–secreting (Th17) CD4+ T cells from bile acid–driven oxidative stress in the small intestine (Cao et al., 2017; Randolph et al., 1998; Tanner et al., 2013). By contrast, the function of MDR1 in CTLs and NK cells has remained controversial (Egashira et al., 1999; Gupta et al., 1992), but has important implications in the design and delivery of vaccines and immunotherapies.

A paucity of genetic mouse models and specific antibodies has hampered a more robust understanding of MDR1 expression and function in vivo. Mice lacking one (Abcb1a−/−; Schinkel et al., 1994) or both (Abcb1a−/−Abcb1b−/−; hereafter referred to as Abcb1a/1b−/−; Schinkel et al., 1997) MDR1 genes are only available on the FVB background, in contrast to most C57BL/6J (B6)–derived transgenic and knockout alleles. In addition, no currently available antibodies detect MDR1 proteins in live mouse cells by flow cytometry. Thus, MDR1 expression is most commonly inferred indirectly, based on the efflux of fluorescent transport substrates, such as rhodamine 123 (Rh123; Ludescher et al., 1992). However useful, these dyes are also substrates for related multidrug transporters (e.g., ABCC1, ABCC3, ABCG2) and are, thus, nonspecific (Strouse et al., 2013; Zhou et al., 2001). To begin addressing these limitations, we previously generated an MDR1-reporter allele (Abcb1aAME) on the B6 background, in which the Abcb1a stop codon was replaced with a bicistronic reporter cassette containing a P2A peptide and a fluorescent transgene, ametrine, to reflect endogenous Abcb1a mRNA levels (Cao et al., 2017).

Using Abcb1aAME/+ reporter mice here, we found that cytolytic lymphocytes, including CD8+ CTLs and NK cells, constitutively express Abcb1a. Runt-related (Runx) transcription factors, which instruct CTL and NK cell development (Cruz-Guilloty et al., 2009; Lotem et al., 2013; Rapp et al., 2017; Taniuchi et al., 2002), bound and maintained chromatin accessibility at the MDR1 (Abcb1a-Abcb1b) locus. Most importantly, endogenous MDR1 expression was required for both effector and memory CTL responses to intracellular pathogens in vivo and acted early after CD8+ T cell activation to suppress oxidative stress, promote mitochondrial function, and enforce survival. These results unmasked a previously unrecognized function of MDR1 in cell-mediated immunity, which could be leveraged to identify endogenously generated MDR1 transport substrates and to inform the safe and rationale use of MDR1 inhibitors in human clinical trials.

Endogenous MDR1 expression across the hematopoietic system

An integrated view of endogenous MDR1 expression throughout the hematopoietic system is lacking and could provide insight into protein function. Therefore, we used heterozygous Abcb1aAME/+ reporter mice to quantify steady-state Abcb1a expression in >100 immune cell types and developmental stages from five major lymphoid and nonlymphoid tissues: bone marrow, thymus, spleen, lung, and small intestine lamina propria (siLP; Fig. 1, A and B; and Table S1). This analysis incorporated 11 high-content (10–13 color) flow cytometry panels and used parallel gating of Abcb1aAME/+ reporter and wild-type B6 subsets (Table S1), to account for variable auto-fluorescence between cell types and to quantify normalized Abcb1a expression (Fig. 1 A).

Abcb1a expression was low or absent throughout most stages of bone marrow hematopoiesis and thymic T cell development, whereas it was expressed in several mature myeloid and lymphoid subsets in the spleen and upregulated in nearly all lymphocytes in siLP (Fig. 1, C and D; and Fig. S1). Elevated Abcb1a expression in the intestine was not a generic feature of all mucosal tissues, as Abcb1a expression in these same cell types was significantly lower in the lung vs. gut (mean Abcb1a expression in lung, 442.8; mean Abcb1a expression in siLP, 844.3; P = 0.0051, paired Student’s t test; Fig. 1, C and D; and Fig. S1 A). Thus, preferential expression of Abcb1a in intestinal lymphocytes is compatible with current models that consider MDR1 a transporter of xenobiotic substances in the gastrointestinal tract (Borst and Schinkel, 2013; Cohen, 1989; Schinkel, 1997).

A second set of results was less consistent with the xenobiotic transport model and revealed that Abcb1a was also constitutively expressed in cytotoxic lymphocytes. First, Abcb1a was upregulated in the bone marrow as common lymphoid progenitors committed to the NK lineage (Fig. 1, C and D). Second, whereas Abcb1a expression was high in early double-negative (e.g., DN1, DN2) thymocytes, it was downregulated in later stages of thymic development and only reemerged late in CD8+, but not CD4+, T cell development (Fig. 1, C and D). This same pattern of Abcb1a expression in developing thymocytes was noted on the Immunological Genome Project’s website (available at: http://www.immgen.org), using RNA-sequencing (RNA-seq) as a readout. Third, naive CD8+ T cells displayed constitutive Abcb1a expression in all peripheral tissues, unlike naive CD4+ T cells (Fig. 1, C and D; and Fig. S1 A). Other lymphocytes with cytotoxic potential—invariant NK T cells and NK1.1+ γδT cells—also displayed constitutive Abcb1a expression (Fig. 1, C and D; and Fig. S1 A). Therefore, Abcb1a expression in mice appears to be both a locally acquired trait of lymphocytes in the intestine and a developmentally imprinted feature in cytotoxic immune cells.

Runx transcription factors promote gene expression at the MDR1 (Abcb1a-Abcb1b) locus

Abcb1a expression was higher in polyclonal CD8+ T cells with central memory or effector memory phenotypes compared with naive cells (Fig. 1 C), suggesting that MDR1 may play a role in the generation or maintenance of memory CTLs. To explore the regulation of Abcb1a expression in antigen-specific CD8+ T cells responding to infection, we crossed Abcb1aAME/+ reporter mice with Thy1.1 congenic/P14 TCR transgenic mice (hereafter referred to as P14 mice) that express a TCR specific for the lymphocytic choriomeningitis virus (LCMV) glycoprotein (gp)33-41 peptide presented on MHC class I (H2-Db; Pircher et al., 1989). Thy1.1+CD8+ naive P14 cells with or without the Abcb1aAME reporter allele were transferred into Thy1.2+ wild-type B6 hosts, which were subsequently infected with the Armstrong strain of LCMV (LCMVArm) to produce an acute infection or were left uninfected (Fig. 2 A). Both control and Abcb1aAME/+ reporter P14 cells underwent marked expansion by day 8 after infection, and both differentiated into expected mixtures of Klrg1+IL-7Rαlo terminal effector (TE) and Klrg1IL-7Rαhi memory precursor (MP) cells (Joshi et al., 2007; Fig. 2 A). Abcb1a expression was maintained in TE and Klrg1IL-7Rαlo early effector cells, relative to naive cells, but was upregulated in both MP and Klrg1+IL-7Rαhi double-positive (DP) effector cells from LCMV-infected recipients (Fig. 2, B and C). Thus, Abcb1a expression is highest in effector CD8+ T cells that develop the potential to give rise to both effector memory and central memory T cells.

Runx transcription factors, particularly Runx1 and Runx3, are essential for the development and maintenance of CD8+ T cells and NK cells (Naito et al., 2011). Runx3 is upregulated during positive selection of CD8+ thymocytes (Liu et al., 2005), consistent with the emergence of Abcb1a expression late in CD8+ T cell thymic development (Fig. 1, C and D). In the periphery, Runx3 promotes naive CD8+ T cell differentiation into both DP and MP CTLs, as well as tissue-resident memory cells, following acute infections (Milner et al., 2017; Wang et al., 2018). As these Runx3-dependent subsets all displayed elevated Abcb1a expression, relative to either naive CD8+ T cells or to many other hematopoietic cell lineages (Fig. 1, C and D), we used retroviral expression of previously validated short hairpin RNAs expressed in a microRNA context (shRNAmirs; Wang et al., 2018) to test whether Runx transcriptional activity is required to maintain CTL-intrinsic expression of Abcb1a in vivo (Fig. 2 D). Consistent with this possibility, depletion of Runx1, Runx3, or the obligate Runx DNA-binding cofactor Cbfb in P14 cells responding to LCMVArm infection each impaired endogenous Abcb1a expression in TE or MP cells (Fig. 2, D–F), although Runx1 depletion had the strongest effect (Fig. 2, D–F). Because Runx1 expression is increased in CTLs following Runx3 depletion (Wang et al., 2018), it seems less likely that Runx1 is sufficient to sustain Abcb1a expression in the absence of Runx3, and more likely that Runx1 and Runx3 act in concert to promote endogenous Abcb1a expression during CTL responses to infection.

The mouse MDR1 locus encodes two functionally redundant transporter proteins (Abcb1a, Abcb1b), with Abcb1a reflecting the mouse orthologue of human ABCB1. Unlike its stable expression in vivo, Abcb1a was lost upon in vitro TCR stimulation of naive CD8+ T cells that were subsequently cultured in high or low IL-2–containing media (Pipkin et al., 2010; Fig. S2, A and B). However, these cells retained MDR1 transport (Rh123 efflux) activity, apparently because of Abcb1b upregulation (Fig. S2, B and C). Loss of Abcb1a expression occurred within 48 h of TCR stimulation and paralleled Runx1 downregulation (Fig. S2 C), further suggesting that both Runx1 and Runx3 are required to maintain Abcb1a expression in CTLs. By contrast, Runx3 was sufficient in the absence of Runx1 to maintain, and even increase, Abcb1b expression in cultured CTLs; Abcb1b upregulation correlated with increased Runx3 protein expression, despite stable Runx3 mRNA expression (Fig. S2 C; data not shown), and shRNAmir-mediated depletion of either Runx3 or Cbfb reduced MDR1-dependent Rh123 efflux and abolished Abcb1b mRNA upregulation upon in vitro expansion (Fig. S2, D and E). These results suggest that interplay among Runx1, Runx3, and T cell activation signals underlies dynamic regulation of MDR1 (Abcb1a, Abcb1b) gene expression in CTLs.

To gain molecular insight into MDR1 gene regulation, we used the published assay for transposable-accessible chromatin (ATAC)-seq and chromatin immunoprecipitation (ChIP)–seq data sets (Kurachi et al., 2014; Lotem et al., 2013; Scott-Browne et al., 2016; Wang et al., 2018) to examine regions of chromatin accessibility and transcription factor binding that develop upon naive CD8 T cell activation in vitro and that are present in mature LCMV-specific CTL subsets in vivo. Multiple regions of chromatin accessibility were observed across the MDR1 (Abcb1a-Abcb1b) locus in naive CD8+ T cells; these were most accessible in memory CTLs (isolated 35 d after LCMVArm infection), much less accessible in “exhausted” CTLs (isolated 20 d after LCMV-clone 13 infection), and reduced in Runx3−/− naive CD8+ T cells after 12 h of TCR stimulation in vitro (Fig. 2 G). Runx3 occupied an accessible region in Abcb1b intron 19 in stimulated, but not resting, CD8+ T cells, together with IRF4 and T-bet—transcription factors that act concertedly with Runx3 to program effector and memory CTL differentiation. IRF4 and T-bet also bound additional regions of Runx3-dependent accessible chromatin in Abcb1a intron 3 (Fig. 2 G). Thus, Runx3 binds and maintains chromatin accessibility at multiple cis-regulatory regions across the MDR1 locus in activated CTLs, which is consistent with the preferential expression of Abcb1a in cytotoxic lymphocytes in vivo and may inform broader mechanisms of MDR1 gene regulation in other normal and malignant cell types.

MDR1 is required for CTL function in vivo

To decipher the function of MDR1 during CTL responses to infection, we backcrossed FVB mice lacking both MDR1 genes (Abcb1a/1b−/−) to pure B6 (12 generations) and then bred B6-derived Abcb1a/1b−/− mice with P14 mice to obtain animals in which gp33-specific CD8+ T cells are null for MDR1 function. We produced 1:1 congenic mixtures of Thy1.1+Thy1.2+ MDR1-sufficient and Thy1.1+Thy1.2 MDR1-deficient naive P14 cells, transferred 5000 mixed cells into wild-type B6 (Thy1.1Thy1.2+) mice, infected recipients with LCMVArm to induce an acute infection, and monitored the accumulation and differentiation of virus-specific CTLs over time (Fig. 3 A). P14 CTLs lacking MDR1 failed to accumulate to wild-type levels in the blood of LCMV-infected mice, with the relative abundance of MDR1-deficient vs. MDR1-sufficient P14 cells decreasing markedly within the first 1–2 weeks of LCMV infection before stabilizing after day 35 (Fig. 3 B). Thus, MDR1 function may be particularly important during the early response of antigen-primed CTLs. Frequencies and absolute numbers of MDR1-deficient P14 cells were also reduced in spleens of infected mice, relative to MDR1-sufficient counterparts; this was evident at the peak of the effector response (day 8), as well as at both early (day 35) and late (day 90) memory time points (Fig. 3, C and D). Loss of MDR1 impacted both Klrg1+IL-7Rαlo and Klrg1IL-7Rαhi CTLs equally (Fig. 3, C and D), further suggesting that MDR1 regulates early events in the activation or expansion of primed CTLs, as opposed to selectively regulating effector or memory CTL differentiation per se.

We also asked whether long-lived memory CTLs that develop in the absence of MDR1 are functional. For this, we FACS-purified Thy1.1+Thy1.2+ MDR1-sufficient or Thy1.1+Thy1.2 MDR1-deficient memory P14 cells from spleens of congenically transferred mice 90 d after LCMV infection, transferred 30,000 cells of each genotype into separate groups of naive B6 recipients, and then challenged these animals with a high dose (80,000 CFUs) of recombinant Listeria monocytogenes expressing gp33 (LM-gp33), which requires exogenous P14 memory cells for protection (Olson et al., 2013); unmanipulated naive B6 mice (receiving no P14 cells) were infected in parallel to confirm memory CTL-dependent protection (Fig. 3 E). Indeed, memory CTL-dependent protection against LM-gp33 challenge also required MDR1, as both naive B6 mice and mice receiving MDR1-deficient P14 memory cells displayed increased morbidity (weight loss) and mortality compared with recipients of MDR1-sufficient P14 memory cells (Fig. 3 F). Among animals surviving until day 5, significantly more L.monocytogenes was recovered from spleens of MDR1-deficient vs. MDR1-sufficient P14 recipients, whereas both the frequency and absolute numbers of MDR1-deficient P14 cells were reduced relative to MDR1-sufficient controls (Fig. 3, G–I). These data suggest that both primary and secondary CTL responses to intracellular pathogens require endogenous MDR1 function. By contrast, we found no evidence that MDR1 regulates the development or homeostasis of naive CD8+ T cells as (1) neither FVB- nor B6-derived MDR1-null mice showed altered frequencies or numbers of naive CD8+ T cells, relative to wild-type littermates; and (2) loss of MDR1 affected neither the engraftment nor the persistence (over 8 d) of naive P14 cells transferred into uninfected recipient mice (data not shown).

MDR1 enforces survival and mitochondrial fitness in early mitotic CTLs

To explore the timing and mechanisms by which endogenous MDR1 function might promote CTL responses to infection, we performed single-cell RNA-seq (scRNA-seq) analysis on congenically cotransferred control and MDR1-deficient P14 cells 8 d after LCMVArm infection. After normalizing gene expression and deconvoluting cell genotypes (based on antibody hashtags), we reduced the dimensionality of the data with Uniform Manifold Approximation and Projection (UMAP; Becht et al., 2018) and performed Louvain clustering. Both control and MDR1-deficient P14 cells segregated into five cell clusters, including early mitotic cells, Klrg1hi TE cells, and Il7rhi MP cells (Fig. 4 A). The expected distribution of control and MDR1-deficient CTLs was observed in all but one cluster (cluster 3), which corresponded to early mitotic CTLs (expressing high levels of Cdk1, Pcna, and Mki67) and which contained significantly fewer than expected MDR1-null cells (chi-square P = 0.396E-06; Fig. 4, B and C). Among this actively dividing CTL cluster, cells lacking MDR1 displayed increased gene expression associated with autophagy and apoptosis (Fig. 4 D); many of these genes were also increased in all MDR1-deficient vs. control CTLs, irrespective of cluster (data not shown). As autophagy is a major pathway through which damaged mitochondria are cleared (Um and Yun, 2017) and cell-intrinsic apoptosis is initiated by the release of cytochrome c from depolarized mitochondrial membranes (Green, 2019), these transcriptional signatures suggest that a primary function of MDR1 may be to promote mitochondrial fitness in recently activated CTLs. Consistent with this possibility, day 8 LCMV-specific CTLs produced more mitochondrial superoxide (mitoSOX) and less ATP in the absence of MDR1 (Fig. 4, E–G).

MDR1 mitigates activation-induced oxidative stress

Elucidating early activation events in small/physiological numbers of antigen-specific CD8+ T cells in vivo is not feasible. Therefore, we used in vitro cell culture to more precisely map the early events in CD8+ T cell activation and proliferation that may be impacted by MDR1 deficiency. Consistent with our in vivo results, naive CD8+ T cells lacking MDR1 displayed reduced accumulation vs. wild-type counterparts upon in vitro differentiation into either effector- or memory-like CTLs (Fig. 5 A). Blocking MDR1 transport activity with elacridar (Hyafil et al., 1993) also reduced the accumulation of wild-type CTLs—to a similar extent as Abcb1a/1b ablation—but had no impact on the already blunted expansion of MDR1-null CTLs (Fig. 5 A). In addition, MDR1 promoted CTL accumulation in a cell-intrinsic manner, which we confirmed using congenic coculture experiments where CD45.1 wild-type and CD45.2 MDR1-deficient CTLs were activated and expanded together in the same tissue culture wells (Fig. S2 F). Reduced accumulation of MDR1-deficient CTLs was associated with increased cell death (Fig. 5 B), whereas MDR1-deficient cells that remained viable displayed wild-type level activation (CD69 upregulation) and proliferation (carboxyfluorescein diacetate succinimidyl ester dilution) kinetics (Fig. S2 G). Also in line with our in vivo results, loss or inhibition of MDR1 in cultured CTLs led to (1) increased gene expression associated with autophagy and apoptosis (as well as hypoxia; Fig. 5 C and Fig. S2 H); and (2) increased mitoSOX and decreased ATP production (Fig. 5, D–F). These same metabolic parameters were unaffected by MDR1 deficiency in resting naive CD8+ T cells (Fig. 5, D–F). These in vitro results corroborate our in vivo findings and reinforce a model in which MDR1 acts early and selectively in activated CTL to enforce metabolic fitness.

Increased mitoSOX production (i.e., oxidative stress) can be both a cause and a consequence of defective ATP synthesis (Martínez-Reyes and Cuezva, 2014; Zorov et al., 2014). Thus, we finally sought to understand whether endogenous MDR1 function more directly regulates oxidation reduction (redox) or bioenergetic pathways. For this, we analyzed the kinetics of oxidative stress and ATP synthesis in recently activated naive CD8+ T cells that have or lack MDR1, and we focused on the synthesis and redox status of glutathione—an endogenous antioxidant cofactor produced by activated T cells to buffer oxidized free radicals (Wang et al., 2011)—to obtain a more comprehensive view of oxidative stress. As expected, TCR stimulation diminished reduced glutathione (GSH) stores in both wild-type and MDR1-deficient CD8+ T cells and within 12 h of TCR stimulation (Fig. 5 G). By contrast, only wild-type cells were capable of replenishing reduced GSH levels over the next 12–36 h (Fig. 5 G). Decreased levels of reduced GSH in activated MDR1-null CD8+ T cells were associated specifically with accumulation of oxidized glutathione dipeptides (GSSG), whereas both total and reduced GSH synthesis remained intact in the absence of MDR1 (Fig. 5 G, data not shown). In addition, alterations in the GSH redox balance preceded deficits in ATP synthesis in MDR1-deficient CTLs, the latter of which did not manifest until 48 h after stimulation (Fig. 5 H). Together, these data suggest that the primary endogenous function of MDR1 in CTLs may be to limit activation-induced oxidative stress.

General conclusions

By establishing contemporary reporter and loss-of-function alleles in mice, we have shown that the MDR1 transporter is highly expressed in CTLs, where it is regulated by interplay between Runx transcription factors and T cell activation signals and necessary for both primary and secondary responses to infection. These data highlight a previously unrecognized endogenous function of MDR1 in cell-mediated immunity and raise questions about the safety and utility of systemic MDR1 inhibitors in cancer therapy. At a molecular level, both in vivo scRNA-seq and in vitro mechanistic studies support a model in which MDR1 function is summoned rapidly after antigen-dependent CTL priming to buffer TCR-driven oxidized free radical production. It is interesting in this regard that production of reactive oxygen/nitrogen species is considered a general response of adaptive lymphocytes to antigen receptor signaling, whereas most noncytotoxic lymphocytes lack constitutive MDR1 expression (Fig. 1). Therefore, it is possible that the metabolic programs of activated CTLs are unique among adaptive lymphocytes in their generation of oxidized free radicals or endogenous MDR1 transport substrates. Finally, decreased mitochondrial metabolism and increased oxidative stress shown here in MDR1-deficient CTLs resembles that previously described in MDR1-deficient enterocytes, tumor cells, and CD4+ T cells exposed to bile acids in the intestine (Cao et al., 2017; Ho et al., 2018; Hwang et al., 2019). Thus, it remains possible that a unifying molecular function of MDR1 may eventually be traced to the transport of a common class of endogenous metabolites that is produced by all cells in a manner commensurate with responses to oxidative stress.

Mice

C57BL/6J (B6)–derived Abcb1aAME/+ reporter mice were generated as previously described (Cao et al., 2017). B6-derived wild-type (stock no. 000664) and Rag1−/− (stock no. 002216) mice were purchased from The Jackson Laboratory. FVB/N (FVB)-derived Abcb1a/1b−/− mice (model no. 1487) were purchased from Taconic and backcrossed for 10 generations under a breeding license to obtain B6-derived Abcb1a/1b−/− mice. B6-derived Abcb1aAME/+ reporter and Abcb1a/1b−/− mice were bred to Thy1.1 congenic/P14 TCR transgenic mice (Wang et al., 2018). All mice were bred and used in experiments in accordance with protocols approved by the Scripps Florida institutional animal care and use committee.

T cell isolation and culture

Naive CD8+ T cells were isolated from spleen and peripheral lymph nodes by magnetic negative isolation (Stem Cell Technologies, Inc.). Purified CD8+ T cells were cultured in DMEM supplemented with 10% heat-inactivated FBS, essential amino acids, nonessential amino acids, sodium pyruvate, arginine/asparagine/folic acid, 10 mM Hepes, 2 mM l-glutamine, 50 µM 2-mercaptoethanol, penicillin-streptomycin, and gentamycin (all from Thermo Fisher Scientific). For in vitro CD8+ T cell activation, magnetically isolated CD8+ T cells were activated at 4 × 105 cells/cm2 and 106 cells/ml in 24-well flat bottom plates coated with anti-CD3 (1 µg/ml) and anti-CD28 (1 µg/ml; BioLegend) after precoating with 200 µg/ml goat anti-hamster IgG (Thermo Fisher Scientific). Activated T cells were removed from antibody-coated wells after 48 h and recultured at 5 × 105 cells/ml in the indicated concentrations of recombinant human IL-2 (rhIL-2; National Institutes of Health Biorepository). For in vitro coculture experiments, 1:1 mixtures of magnetically enriched CD45.1 wild-type and CD45.2 Abcb1a/1b−/− naive CD8+ T cells were activated and cultured as above. In some experiments, MDR1 was blocked by adding elacridar (100 nM; Sigma-Aldrich) immediately before T cell activation and again every 36–48 h.

Retroviral plasmids and transductions

GFP- or ametrine-expressing murine retroviral vectors containing shRNAmirs against Cd4, Runx1, Runx3, and Cbfb have been described and validated previously (Wang et al., 2018). Generation of retroviral particles and T cell transductions were performed as previously described (Cao et al., 2017) using TransIT-LT1 transfection reagent (Mirus). Briefly, CD8+ T cells were transduced 24 h after activation; culture media was replaced with viral supernatant containing 10 µg/ml polybrene and centrifuged at 2,000 rpm for 1 h at room temperature. For in vitro culture experiments, transduced cells were kept in a 37°C incubator until 48 h after activation, after which time retroviral supernatants were replaced with fresh media containing 10 or 100 U/ml of rhIL-2, as indicated. For adoptive transfer of transduced P14 cells, T cell cultures were incubated for 3 h in a 37°C CO2 incubator after retroviral infections and centrifugation, and 5 × 105 cells were transferred into naive B6 mice for LCMV infection as above.

Flow cytometry

All surface and intracellular FACS stains were performed at 4°C for 30 min and then washed with PBS before data acquisition. For Rh123 efflux, cells were labeled with 0.1 µM Rh123 (Sigma-Aldrich) for 30 min at 4°C in ice-cold complete media. Labeled cells were then washed with ice-cold fresh media, resuspended in fresh room temperature media, and incubated for 30–60 min in a 37°C incubator. An aliquot of labeled cells was treated with elacridar (100 nM) to monitor background efflux activity. Any surface stains performed after Rh123 efflux were conducted at 4°C as above. Mitochondrial superoxide was quantified via FACS analysis of mitoSOX Red staining (Thermo Fisher Scientific; excitation/emission = 510/580 nm). Briefly, cells were labeled concomitantly with mitoSOX Red (5 µM) and cell surface FACS antibodies in PBS and incubated for 15 min in a 37°C incubator protected from light. To confirm specificity, an aliquot of cells was treated for 15 min with the mitochondrial complex III inhibitor, Antimycin A (1 µM; Abcam), at 37°C before mitoSOX staining. Labeled cells were washed twice and resuspended with warm PBS before acquisition on the cytometer. Anti-mouse antibodies used for FACS analysis included Alexa700-CD34, BV605-CD11b, BV711-CD127, PE-CF594-CD135, BV711-CD138, FITC-CD1d, Percp-Cy5.5-CD21/CD35, PE-CF594-CD25, PE-Cy7-CD27, BV421-CD5, Alexa700-CD8a, BUV395-CD8a, PE-CF594-CD93, PE-Cy7-CD95, BV421-Gr-1, Percp-Cy5.5-IgD, BV650-IgM, BUV395-Ly-6G, BV650-MHCII, BUV395-NK1.1, BV605-Sca1, PE-CF594-Siglec-F (from BD Biosciences); Alexa700-B220, FITC-B220, BV421-c-Kit, PE-c-Kit, BV421-CD103, Percp-Cy5.5-CD115, FITC-CD11b, BV605-CD11b, BV650-CD150, BUV395-CD19, FITC-CD19, FITC-CD25, BV650-CD3, FITC-CD3, BV711-CD4, Percp-Cy5.5-CD4, PE-Cy7-CD44, BV650-CD44, Percp-Cy5.5-CD48, Percp-Cy5.5-CD62L, PE-CD8A, PE-CX3CR1, FITC-Gr-1, BV421-Lin Cocktail, BV711-Ly-6C, Percp-Cy5.5-TCRβ, PE-Cy7-TCRβ, BV605-TCRβ, FITC-CD69, BV421-CD90.1, APC-KLRG1 (from BioLegend); and PE-Cy7-CD11c, PE-Cy7-CD122, FITC-CD16/CD32, FITC-CD24, and PE-CD38 (from eBioscience). Vital dyes used in this study included fixable viability eFluor 506 and eFluor 660 (eBioscience). All FACS data were acquired on LSRII or FACSCanto II instruments (BD Biosciences), and analyzed using FlowJo 9 or 10 software (TreeStar, Inc.). Primary FACS data of endogenous Abcb1aAme reporter expression are available at ImmPort (available at: https://www.immport.org) under the accession no. SDY1595.

Fluorescent microscopy

Spleens from B6 or Abcb1aAME/+ reporter mice were washed with PBS, embedded in Tissue-Tek O.C.T. Compound (Thermo Fisher Scientific) and flash frozen in liquid nitrogen. Tissues were cut on a CM1950 cryostat (Leica) at 5 µm thickness. Slides were postfixed in cold acetone (Sigma-Aldrich), blocked with 5% normal donkey serum (Sigma-Aldrich), and stained with anti-CD11c BOND RTU Primary (Leica Biosystems) with Alexa594 F(ab′)2 and donkey anti-rabbit IgG (H+L) secondary (1:500; Jackson Immunoresearch). Sections were mounted with ProLong Diamond Mountant (Thermo Fisher Scientific). Images were captured on an Eclipse microscope (Nikon).

LCMV/L.monocytogenes infections

5 × 103 naive CD8+ T cells isolated from spleens of wild-type P14, P14.Abcb1aAME/+, or P14.Abcb1a/1b−/− mice were injected (i.v.) into 6–8-wk-old wild-type B6 recipients. The next day, recipient animals were either left uninfected or infected (i.p.) with 2 × 105 PFU LCMVArm. In some experiments, naive CD8+ T cells from spleens of P14.Abcb1aAME/+ reporter mice were activated and transduced in vitro (as above) with GFP-expressing retroviral vectors containing control (shCD4), or Runx-targeting (shRunx1, shRunx3, shCbfb) shRNAmirs before in vivo transfer. 5 × 105 transduced P14 cells were injected (i.v.) into 6–8-wk-old wild-type B6 recipients, which were infected (i.p.) 1 h later with 2 × 105 PFU LCMVArm. For congenic T cell transfer and infection experiments, naive CD8+ T cells isolated from spleens of wild-type P14 (Thy1.1/Th1.2 heterozygous) or P14.Abcb1a/1b−/− (Thy1.1 homozygous) mice were mixed in a 1:1: ratio, and 5 × 103 mixed cells were coadoptively transferred (i.v.) into 6–8-wk-old wild-type B6 recipients before infection (i.p.) with 2 × 105 PFU LCMVArm 24 h later. For LM-gp33 challenge experiments, viable CD8+CD44hi P14 cells with (Thy1.1/Th1.2 heterozygous) or without (Thy1.1 homozygous) MDR1 were FACS-purified from spleens of LCMV-infected mice on day 90 after infection. 3 × 104 cells of each genotype were injected (i.v.) into separate groups of naive (6–8-wk-old) wild-type B6 recipients, which were then infected (i.v.) with 8 × 104 CFU of LM-gp33 24 h later. Mice were sacrificed 5 d after infection. Ex vivo T cell analysis was performed by FACS; bacterial titers were quantified by plating serial dilutions of spleen lysates on brain heart infusion agar as described previously (Olson et al., 2013).

ATP assays

ATP concentrations were quantified using the CellTiter-Glo Luminescent Cell Viability kit (Promega), in accordance with manufacturer’s instructions. Briefly, cells were resuspended in complete media, and 25 µl of titrating cells numbers (104, 5 × 103, 2.5 × 103, and 1.25 × 103) were pipetted, in triplicate, into standard white opaque 384-well microplates (Thermo Fisher Scientific). Control wells were prepared containing complete media without cells to determine background luminescence. Tenfold serial dilutions of ATP (100 mM–0.1 mM) were prepared fresh in complete media to generate ATP standard curves for each plate. An equal volume of 25 µl CellTiter-Glo reagent was added to each sample and shaken for 2 min on an orbital shaker to induce cell lysis. The plate was incubated for an additional 10 min at room temperature to stabilize the luminescent signal. Luminescence was recorded, in duplicate, using a SpectraMax M5 microplate reader (Molecular Devices) with an integration time of 0.5 s per well. ATP concentrations (per 5 × 103 cells) were calculated using the mean corrected luminescence signals, after background subtraction and normalization for cell viability (determined by flow cytometry).

GSH assays

Total, reduced, and oxidized (GSSG) GSH levels were quantified using the GSH/GSSG Ratio Detection Assay Kit (Abcam) following the manufacturer’s instructions. Briefly, 104 cells were lysed and deproteinized using a Deproteinizing (TCA) Sample Preparation Kit (Abcam). 25 µl of the deproteinized lysates were pipetted, in duplicate, into black side/clear bottom 384-well microplates (Greiner Bio-One). Twofold serial dilutions of reduced GSH (10 µM–0.1563 µM) and GSSG (5 µM–0.0781 µM) were prepared fresh for standard curves. For detection of reduced GSH, 25 µl of GSH assay mixture was added into each reduced GSH standard and sample wells. For detection of total GSH, 25 µl of total GSH assay mixture was added into each oxidized GSSG standard and sample wells. Plates were incubated for 30 min in the dark, and fluorescence was recorded at excitation/emission = 490/520 nm on the Molecular Devices SpectraMax M5 microplate reader. Reduced GSH and total GSH concentrations were calculated by fitting corrected (background subtracted) mean fluorescence intensities (MFIs) to standard curves; oxidized GSSG concentrations were calculated based on the difference between total GSH and reduced GSH.

Quantitative PCR (qPCR)

RNA was isolated from cultured or ex vivo isolated cells using RNeasy columns (Qiagen) as per the manufacturer’s instructions using on-column DNase treatment. RNA was used to synthesize cDNA via a high capacity cDNA reverse transcription kit (Life Technologies). Taqman qPCR was performed on a StepOnePlus real time PCR instrument (Life Technologies) using the following Taqman primer/probe sets (Life Technologies): Abcb1a (Mm00440761_m1), Abcb1b (Mm00440736_m1), Runx1 (Mm01213404_m1), Runx3 (Mm01213404_m1), and Gapdh (mm99999915_g1).

Bioinformatics analyses

ChIP-seq

Raw sequencing reads for Runx3 (GSE50131; Lotem et al., 2013), Tbx21 (GSE72408; Dominguez et al., 2015), and Irf4 (GSE54191; Kurachi et al., 2014) were downloaded from the Sequence Read Archive database. Reads were aligned to UCSC mm10 with Bowtie2 (Langmead and Salzberg, 2012), and peaks were called in MACS (Zhang et al., 2008) using base settings. Bigwig files of binding peaks were generated and visualized on the UCSC genome browser. Peaks were filtered in R-studio (using GenomeInfoDb and ChIP blacklist GenomicRanges packages) to remove reads with alternative annotations, mitochondrial DNA, or blacklist sites (Carroll et al., 2014).

ATAC-seq

ATAC-seq data were downloaded from Gene Expression Omnibus (GEO; GSE111149, GSE88987), aligned to UCSC mm10 with Bowtie2 (Langmead and Salzberg, 2012) and visualized on the UCSC genome browser.

RNA-Seq

Next-generation RNA-seq was performed on three biologically independent replicates of B6 wild-type or Abcb1a/1b−/− CD8+ T cells activated and expanded for 6 d in 10 U/ml of rhIL-2–containing media with or without elacridar. High quality total RNA (RNA integrity number >9 by Bioanalyzer; Agilent Technologies) was isolated using RNeasy columns with on-column DNase treatment (Qiagen) and processed with the NEBNext Ultra II ribosomal RNA-depletion module human/mouse/rat (New England BioLabs Inc.) using 400 ng total RNA as input. The ribosomal RNA-depleted RNA was processed using the NEBNext Ultra II directional RNA kit for Illumina (New England BioLabs Inc.) to generate strand-specific libraries. Final libraries were amplified by PCR (10 cycles), size selected using beads to enrich for fragments >400 bp, and purified using Ampure XP beads (Beckman Coulter). Final libraries were validated on the Agilent 2100 Bioanalyzer DNA chips and quantified on the Qubit (Thermo Fisher Scientific). Libraries were pooled at equimolar ratios and loaded onto the Illumina NextSeq500 flow cell and sequenced using 40-bp paired-end chemistry. Approximately 25 million reads pass filter (base quality score >Q30; suggesting <1 error in 1,000 bp) were generated per library. Raw sequencing reads (fastq files) were mapped to the mm10 genome, and normalized transcript abundance, displayed as transcripts per million (TPM), were quantified using Salmon (Patro et al., 2017). Differential gene expression was analyzed using DESeq2 and visualized using the MultiPlot studio module in GenePattern (https://cloud.genepattern.org/gp). Principle component analysis(PCA)of normalized TPM data were computed in R script. Gene ontology (GO) term enrichment analysis (using differentially expressed gene sets defined by DESeq2) was performed using the R package ClusterProfiler (Yu et al., 2012); heatmaps of representative genes within significantly enriched GO pathways were visualized using Morpheus (https://software.broadinstitute.org/morpheus/).

scRNA-seq

scRNA-seq analyses was performed on FACS-purified Thy1.1+Thy1.2+ MDR1-sufficient and Thy1.1+Thy1.2 MDR1-deficient P14 CD8+ T cells from 10 congenically transferred and LCMV-infected wild-type B6 mice at day 8 after infection. P14 cells with or without MDR1 were barcoded with TotalSeq-A anti-mouse hashtag antibodies (BioLegend) to enable sample multiplexing and then counted on the Countess II FL automated cell counter (Thermo Fisher Scientific) using Trypan blue to determine viability and cell number (Stoeckius et al., 2018). Samples were pooled equally and adjusted to a final concentration between 700–1,200 cells/µl; 35,000 total cells were loaded into one lane of the single-cell A chip kit (P/N 1000009). Single-cell gel beads in emulsion (GEM) were generated using the 10X Chromium Controller (10X Genomics). Single-cell GEM and sequencing libraries were generated using the Single-Cell 3′ Library and Gel Bead Kit V2 (P/N 120267) according to the manufacturer’s recommendations. Briefly, GEM-RT was performed on the recovered GEM in a C1000 Touch Thermal Cycler with the following program: 53°C for 45 min, 85°C for 5 min, and held at 4°C. After RT, single-strand cDNA was recovered with the DynaBeads MyOne Silane Beads (Thermo Fisher Scientific; P/N 2000048). The cDNA was amplified with the hashtag oligonucleotide (HTO) additive primer 5′-GTG​ACT​GGA​GTT​CAG​ACG​TGT​GC*T*C-3′ (asterisk indicates a phosphorothioate bond) added at 0.2 µM using the C1000 Touch Thermal Cycler with the 96-Deep Well Reaction Module: 98°C for 3 min; cycled 8×: 98°C for 15 s, 67°C for 20 s, and 72°C for 1 min; 72°C for 1 min; and held at 4°C. Amplified cDNA products were subjected to cleanup, using the SPRIselect Reagent Kit with 0.6× SPRI (Beckman Coulter; P/N B23317) to recover full-length cDNA, and supernatant was purified with 2× SPRI to recover the HTO product. cDNA products were then quantified by Qubit (Thermo Fisher Scientific) and visualized on the Agilent 2100 Bioanalyzer DNA high sensitivity chip. The HTO fraction was amplified with 2× Kapa HiFi HotStart ready mix (Kapa Biosystems, P/N KK2601) and the following oligos: P5 5′-AAT​GAT​ACG​GCG​ACC​ACC​GAG​ATC​TAC​ACT​CTT​TCC​CTA​CAC​GAC​GC*T*C-3′ and P7 index primer 5′-CAA​GCA​GAA​GAC​GGC​ATA​CGA​GAT​CGA​GTA​ATG​TGA​CTG​GAG​TTC​AGA​CGT​GT*G*C-3′. The following cycling conditions were used in the C1000 Touch Thermal Cycler to amplify the hashtag library: 95°C for 3 min; cycled 12×: 95°C for 20 s, 64°C for 30 s, and 72°C for 20 s; 72°C for 5 min; and held at 4°C. HTO library was purified with 1.6× SPRI beads. The cDNA was converted to indexed sequencing libraries according to the protocol in the Single-Cell 3′ Reagent Kit v2 User Guide (10X Genomics; CG00052 Rev E) that included the following steps: (1) enzymatic fragmentation, end-repair, and A-tailing; (2) double-sided size selection (0.6× SPRI and 0.8× SPRI to supernatant); (3) adaptor ligation; (4) postligation cleanup (0.8X SPRI); (5) index PCR with the Chromium i7 multiplex kit (P/N 120262) on the C1000 Touch Thermal cycler with the 96-Deep Well Reaction Module: 98°C for 45 s; cycled 15×: 98°C for 20 s, 54°C for 30 s, and 72°C for 20 s; 72°C for 1 min; and held at 4°C; and (6) post-PCR cleanup with double-sided SPRI: 0.6× and 0.8×. The final libraries were run on the Bioanalyzer and quantified with the NEBNext library quantification kit (P/N E7630). The cDNA and HTO libraries were pooled at 10:1 and sequenced for a sequencing depth of 50,000 reads per cell for the cDNA and 2000 reads per cell for the HTO library. Libraries were loaded (1.8 pM) on the Illumina NextSeq500 and sequenced with Read 1 25-bp, index i7 8-bp, and Read 2 98-bp. The Cell Ranger Software Suite pipeline was used to demultiplex samples, align reads to the mm10 reference genome, and generate barcoded data matrices (https://support.10xgenomics.com/single-cell-gene-expression/software/pipelines/latest/what-is-cell-ranger). Single-cell transcriptomic data were analyzed using the Python-based SCANPY framework (Wolf et al., 2018). Cells expressing <500 genes and genes detected in <3 cells were removed from analysis. Normalization of each cell by total counts over all genes was performed so that every cell had the same total of 10,000 reads, allowing counts to be comparable among cells yielding a total of 4,201 P14 CD8+ T cells—1,503 MDR1-sufficient and 2,698 MDR1-deficient—averaging 1,682 gene counts/cell. Normalized counts were then log-transformed and used for UMAP to visualize phenotypic clusters identified by the Louvain algorithm (Becht et al., 2018). Differential gene expression was analyzed using DESeq2 with the ZINBWAVE R package to account for the excess of zero read counts and better fit the data to a zero inflated negative binomial distribution (Van den Berge et al., 2018) and visualized with ggplot2. GO term enrichment analysis was performed using ClusterProfiler (Yu et al., 2012), as above. The scRNA-seq and bulk RNA-seq data sets presented in this paper are publicly available in the GEO database under accession no. GSE143141.

Quantification and statistical analyses

Statistical analyses were performed using Prism (GraphPad). P values were determined by paired or unpaired Student’s t tests or one-way or two-way ANOVA, as appropriate, and indicated throughout the figure legends. Differences were considered significant when P ≤ 0.05 (*P ≤ 0.05, **P ≤ 0.01, ***P ≤ 0.001, ****P ≤ 0.0001). Significance levels are specified throughout the figure legends. Data are shown as mean values ± SEM throughout.

Online supplemental material

Fig. S1 shows endogenous Abcb1aAme reporter expression (related to Fig. 1) in hematopoietic cells from different tissues (bone marrow, spleen, thymus, lung, siLP). Fig. S2 shows (1) regulation of MDR1 (Abcb1a, Abcb1b) gene expression in cultured CTL by Runx transcription factors (related to Fig. 2); (2) activation, proliferation, and accumulation of CD45.1 MDR1-sufficient and CD45.2 MDR1-deficient CTL in congenic coculture experiments (related to Fig. 5); and (3) PCA and differential gene expression, determined by RNA-seq in day 6 in vitro–expanded or MDR1-deficient memory-like CTL treated with or without elacridar (related to Fig. 5). Table S1 describes the FACS staining panels used to analyze Abcb1aAme reporter expression, including the specific cell types analyzed in each panel and associated gating strategies/phenotypes (related to Fig. 1 and Fig. S1).

The authors are grateful for technical support from the The Scripps Research Institute (TSRI)–Florida Flow Cytometry and Genomics Core Facilities, and for the employees of the TSRI-Florida Animal Resource Center for assistance with animal care. We thank Drs. Ananda Goldrath, Gustavo Martinez, and Joseph G. Crompton for critical discussions and review of the manuscript.

This work was supported by TSRI-Florida via the State of Florida, the National Institutes of Health (grants R01AI118931-01 and R21AI119728 to M.S. Sundrud; grants 5R01AI095634-06 and 5U19 AI109976-05 to M.E. Pipkin; grant R01HL125816 to S.B. Koralov), and the Crohn’s and Colitis Foundation (Senior Research Award 422515 to M.S. Sundrud). A. Sun received medical scientist training grants from the National Institutes of Health (5T32GM007308-36 and 5T32CA009161-39).

Author contributions: M.L. Chen, A. Eliason, K.M. Mendez, A.J. Getzler, C. Mukori, and S. Tsuda performed in vitro and in vivo experiments, analyzed data, and revised the manuscript. H. Diao performed bioinformatic analysis of ChIP-seq and ATAC-seq data. A. Sun, W. Cao and S.Y. Kim generated Abcb1aAME/+ reporter mice. N.E. Bruno provided critical discussions and revised the manuscript. M.E. Pipkin provided key reagents and discussions, supervised in vivo infection experiments, analyzed data, and revised the manuscript. S.B. Koralov supervised the generation of Abcb1aAME/+ reporter mice, provided critical discussions, and revised the manuscript. M.S. Sundrud conceived and supervised the study, designed experiments, analyzed data, and wrote the manuscript.

Bauer
,
B.
,
A.M.
Hartz
,
G.
Fricker
, and
D.S.
Miller
.
2005
.
Modulation of p-glycoprotein transport function at the blood-brain barrier
.
Exp. Biol. Med. (Maywood).
230
:
118
127
.
Becht
,
E.
,
L.
McInnes
,
J.
Healy
,
C.A.
Dutertre
,
I.W.H.
Kwok
,
L.G.
Ng
,
F.
Ginhoux
, and
E.W.
Newell
.
2018
.
Dimensionality reduction for visualizing single-cell data using UMAP
.
Nat. Biotechnol.
Borst
,
P.
, and
A.H.
Schinkel
.
2013
.
P-glycoprotein ABCB1: a major player in drug handling by mammals
.
J. Clin. Invest.
123
:
4131
4133
.
Cao
,
W.
,
H.
Kayama
,
M.L.
Chen
,
A.
Delmas
,
A.
Sun
,
S.Y.
Kim
,
E.S.
Rangarajan
,
K.
McKevitt
,
A.P.
Beck
,
C.B.
Jackson
, et al
.
2017
.
The xenobiotic transporter Mdr1 enforces T cell homeostasis in the presence of intestinal bile acids
.
Immunity.
47
:
1182
1196.e10
.
Carroll
,
T.S.
,
Z.
Liang
,
R.
Salama
,
R.
Stark
, and
I.
de Santiago
.
2014
.
Impact of artifact removal on ChIP quality metrics in ChIP-seq and ChIP-exo data
.
Front. Genet.
5
:
75
.
Chaudhary
,
P.M.
, and
I.B.
Roninson
.
1991
.
Expression and activity of P-glycoprotein, a multidrug efflux pump, in human hematopoietic stem cells
.
Cell.
66
:
85
94
.
Chaudhary
,
P.M.
,
E.B.
Mechetner
, and
I.B.
Roninson
.
1992
.
Expression and activity of the multidrug resistance P-glycoprotein in human peripheral blood lymphocytes
.
Blood.
80
:
2735
2739
.
Cohen
,
I.
1989
.
Teaching medical students to use simple language when talking with patients
.
J. La. State Med. Soc.
141
:
33
36
.
Cruz-Guilloty
,
F.
,
M.E.
Pipkin
,
I.M.
Djuretic
,
D.
Levanon
,
J.
Lotem
,
M.G.
Lichtenheld
,
Y.
Groner
, and
A.
Rao
.
2009
.
Runx3 and T-box proteins cooperate to establish the transcriptional program of effector CTLs
.
J. Exp. Med.
206
:
51
59
.
Dominguez
,
C.X.
,
R.A.
Amezquita
,
T.
Guan
,
H.D.
Marshall
,
N.S.
Joshi
,
S.H.
Kleinstein
, and
S.M.
Kaech
.
2015
.
The transcription factors ZEB2 and T-bet cooperate to program cytotoxic T cell terminal differentiation in response to LCMV viral infection
.
J. Exp. Med.
212
:
2041
2056
.
Egashira
,
M.
,
N.
Kawamata
,
K.
Sugimoto
,
T.
Kaneko
, and
K.
Oshimi
.
1999
.
P-glycoprotein expression on normal and abnormally expanded natural killer cells and inhibition of P-glycoprotein function by cyclosporin A and its analogue, PSC833
.
Blood.
93
:
599
606
.
Green
,
D.R.
2019
.
The coming decade of cell death research: five riddles
.
Cell.
177
:
1094
1107
.
Gupta
,
S.
,
C.H.
Kim
,
T.
Tsuruo
, and
S.
Gollapudi
.
1992
.
Preferential expression and activity of multidrug resistance gene 1 product (P-glycoprotein), a functionally active efflux pump, in human CD8+ T cells: a role in cytotoxic effector function
.
J. Clin. Immunol.
12
:
451
458
.
Ho
,
G.T.
,
R.E.
Aird
,
B.
Liu
,
R.K.
Boyapati
,
N.A.
Kennedy
,
D.A.
Dorward
,
C.L.
Noble
,
T.
Shimizu
,
R.N.
Carter
,
E.T.S.
Chew
, et al
.
2018
.
MDR1 deficiency impairs mitochondrial homeostasis and promotes intestinal inflammation
.
Mucosal Immunol.
11
:
120
130
.
Hwang
,
S.H.
,
M.C.
Kim
,
S.
Ji
,
Y.
Yang
,
Y.
Jeong
, and
Y.
Kim
.
2019
.
Glucose starvation induces resistance to metformin through the elevation of mitochondrial multidrug resistance protein 1
.
Cancer Sci.
110
:
1256
1267
.
Hyafil
,
F.
,
C.
Vergely
,
P.
Du Vignaud
, and
T.
Grand-Perret
.
1993
.
In vitro and in vivo reversal of multidrug resistance by GF120918, an acridonecarboxamide derivative
.
Cancer Res.
53
:
4595
4602
.
Joshi
,
N.S.
,
W.
Cui
,
A.
Chandele
,
H.K.
Lee
,
D.R.
Urso
,
J.
Hagman
,
L.
Gapin
, and
S.M.
Kaech
.
2007
.
Inflammation directs memory precursor and short-lived effector CD8(+) T cell fates via the graded expression of T-bet transcription factor
.
Immunity.
27
:
281
295
.
Kelly
,
R.J.
,
R.W.
Robey
,
C.C.
Chen
,
D.
Draper
,
V.
Luchenko
,
D.
Barnett
,
R.K.
Oldham
,
Z.
Caluag
,
A.R.
Frye
,
S.M.
Steinberg
, et al
.
2012
.
A pharmacodynamic study of the P-glycoprotein antagonist CBT-1® in combination with paclitaxel in solid tumors
.
Oncologist.
17
:
512
e523
.
Kurachi
,
M.
,
R.A.
Barnitz
,
N.
Yosef
,
P.M.
Odorizzi
,
M.A.
DiIorio
,
M.E.
Lemieux
,
K.
Yates
,
J.
Godec
,
M.G.
Klatt
,
A.
Regev
, et al
.
2014
.
The transcription factor BATF operates as an essential differentiation checkpoint in early effector CD8+ T cells
.
Nat. Immunol.
15
:
373
383
.
Langmead
,
B.
, and
S.L.
Salzberg
.
2012
.
Fast gapped-read alignment with Bowtie 2
.
Nat. Methods.
9
:
357
359
.
Liu
,
X.
,
B.J.
Taylor
,
G.
Sun
, and
R.
Bosselut
.
2005
.
Analyzing expression of perforin, Runx3, and Thpok genes during positive selection reveals activation of CD8-differentiation programs by MHC II-signaled thymocytes
.
J. Immunol.
175
:
4465
4474
.
Lotem
,
J.
,
D.
Levanon
,
V.
Negreanu
,
D.
Leshkowitz
,
G.
Friedlander
, and
Y.
Groner
.
2013
.
Runx3-mediated transcriptional program in cytotoxic lymphocytes
.
PLoS One.
8
:e80467.
Ludescher
,
C.
,
J.
Thaler
,
D.
Drach
,
J.
Drach
,
M.
Spitaler
,
C.
Gattringer
,
H.
Huber
, and
J.
Hofmann
.
1992
.
Detection of activity of P-glycoprotein in human tumour samples using rhodamine 123
.
Br. J. Haematol.
82
:
161
168
.
Martínez-Reyes
,
I.
, and
J.M.
Cuezva
.
2014
.
The H(+)-ATP synthase: a gate to ROS-mediated cell death or cell survival
.
Biochim. Biophys. Acta.
1837
:
1099
1112
.
Milner
,
J.J.
,
C.
Toma
,
B.
Yu
,
K.
Zhang
,
K.
Omilusik
,
A.T.
Phan
,
D.
Wang
,
A.J.
Getzler
,
T.
Nguyen
,
S.
Crotty
, et al
.
2017
.
Runx3 programs CD8+ T cell residency in non-lymphoid tissues and tumours
.
Nature.
552
:
253
257
.
Naito
,
T.
,
H.
Tanaka
,
Y.
Naoe
, and
I.
Taniuchi
.
2011
.
Transcriptional control of T-cell development
.
Int. Immunol.
23
:
661
668
.
O’Brien
,
M.M.
,
N.J.
Lacayo
,
B.L.
Lum
,
S.
Kshirsagar
,
S.
Buck
,
Y.
Ravindranath
,
M.
Bernstein
,
H.
Weinstein
,
M.N.
Chang
,
R.J.
Arceci
, et al
.
2010
.
Phase I study of valspodar (PSC-833) with mitoxantrone and etoposide in refractory and relapsed pediatric acute leukemia: a report from the Children’s Oncology Group
.
Pediatr. Blood Cancer.
54
:
694
702
.
Olson
,
J.A.
,
C.
McDonald-Hyman
,
S.C.
Jameson
, and
S.E.
Hamilton
.
2013
.
Effector-like CD8+ T cells in the memory population mediate potent protective immunity
.
Immunity.
38
:
1250
1260
.
Patro
,
R.
,
G.
Duggal
,
M.I.
Love
,
R.A.
Irizarry
, and
C.
Kingsford
.
2017
.
Salmon provides fast and bias-aware quantification of transcript expression
.
Nat. Methods.
14
:
417
419
.
Peck
,
R.A.
,
J.
Hewett
,
M.W.
Harding
,
Y.M.
Wang
,
P.R.
Chaturvedi
,
A.
Bhatnagar
,
H.
Ziessman
,
F.
Atkins
, and
M.J.
Hawkins
.
2001
.
Phase I and pharmacokinetic study of the novel MDR1 and MRP1 inhibitor biricodar administered alone and in combination with doxorubicin
.
J. Clin. Oncol.
19
:
3130
3141
.
Pipkin
,
M.E.
,
J.A.
Sacks
,
F.
Cruz-Guilloty
,
M.G.
Lichtenheld
,
M.J.
Bevan
, and
A.
Rao
.
2010
.
Interleukin-2 and inflammation induce distinct transcriptional programs that promote the differentiation of effector cytolytic T cells
.
Immunity
.
32
(
1
):
79
90
.
Pircher
,
H.
,
K.
Bürki
,
R.
Lang
,
H.
Hengartner
, and
R.M.
Zinkernagel
.
1989
.
Tolerance induction in double specific T-cell receptor transgenic mice varies with antigen
.
Nature.
342
:
559
561
.
Randolph
,
G.J.
,
S.
Beaulieu
,
M.
Pope
,
I.
Sugawara
,
L.
Hoffman
,
R.M.
Steinman
, and
W.A.
Muller
.
1998
.
A physiologic function for p-glycoprotein (MDR-1) during the migration of dendritic cells from skin via afferent lymphatic vessels
.
Proc. Natl. Acad. Sci. USA.
95
:
6924
6929
.
Rapp
,
M.
,
C.M.
Lau
,
N.M.
Adams
,
O.E.
Weizman
,
T.E.
O’Sullivan
,
C.D.
Geary
, and
J.C.
Sun
.
2017
.
Core-binding factor β and Runx transcription factors promote adaptive natural killer cell responses
.
Sci. Immunol.
2
:
2
.
Schinkel
,
A.H.
1997
.
The physiological function of drug-transporting P-glycoproteins
.
Semin. Cancer Biol.
8
:
161
170
.
Schinkel
,
A.H.
,
J.J.
Smit
,
O.
van Tellingen
,
J.H.
Beijnen
,
E.
Wagenaar
,
L.
van Deemter
,
C.A.
Mol
,
M.A.
van der Valk
,
E.C.
Robanus-Maandag
,
H.P.
te Riele
, et al
.
1994
.
Disruption of the mouse mdr1a P-glycoprotein gene leads to a deficiency in the blood-brain barrier and to increased sensitivity to drugs
.
Cell.
77
:
491
502
.
Schinkel
,
A.H.
,
E.
Wagenaar
,
L.
van Deemter
,
C.A.
Mol
, and
P.
Borst
.
1995
.
Absence of the mdr1a P-Glycoprotein in mice affects tissue distribution and pharmacokinetics of dexamethasone, digoxin, and cyclosporin A
.
J. Clin. Invest.
96
:
1698
1705
.
Schinkel
,
A.H.
,
U.
Mayer
,
E.
Wagenaar
,
C.A.
Mol
,
L.
van Deemter
,
J.J.
Smit
,
M.A.
van der Valk
,
A.C.
Voordouw
,
H.
Spits
,
O.
van Tellingen
, et al
.
1997
.
Normal viability and altered pharmacokinetics in mice lacking mdr1-type (drug-transporting) P-glycoproteins
.
Proc. Natl. Acad. Sci. USA.
94
:
4028
4033
.
Scott-Browne
,
J.P.
,
I.F.
López-Moyado
,
S.
Trifari
,
V.
Wong
,
L.
Chavez
,
A.
Rao
, and
R.M.
Pereira
.
2016
.
Dynamic changes in chromatin accessibility occur in CD8+ T cells responding to viral infection
.
Immunity.
45
:
1327
1340
.
Seiden
,
M.V.
,
K.D.
Swenerton
,
U.
Matulonis
,
S.
Campos
,
P.
Rose
,
G.
Batist
,
E.
Ette
,
V.
Garg
,
A.
Fuller
,
M.W.
Harding
, and
D.
Charpentier
.
2002
.
A phase II study of the MDR inhibitor biricodar (INCEL, VX-710) and paclitaxel in women with advanced ovarian cancer refractory to paclitaxel therapy
.
Gynecol. Oncol.
86
:
302
310
.
Stoeckius
,
M.
,
S.
Zheng
,
B.
Houck-Loomis
,
S.
Hao
,
B.Z.
Yeung
,
W.M.
Mauck
III
,
P.
Smibert
, and
R.
Satija
.
2018
.
Cell Hashing with barcoded antibodies enables multiplexing and doublet detection for single cell genomics
.
Genome Biol.
19
:
224
.
Strouse
,
J.J.
,
I.
Ivnitski-Steele
,
A.
Waller
,
S.M.
Young
,
D.
Perez
,
A.M.
Evangelisti
,
O.
Ursu
,
C.G.
Bologa
,
M.B.
Carter
,
V.M.
Salas
, et al
.
2013
.
Fluorescent substrates for flow cytometric evaluation of efflux inhibition in ABCB1, ABCC1, and ABCG2 transporters
.
Anal. Biochem.
437
:
77
87
.
Sugawara
,
I.
,
I.
Kataoka
,
Y.
Morishita
,
H.
Hamada
,
T.
Tsuruo
,
S.
Itoyama
, and
S.
Mori
.
1988
.
Tissue distribution of P-glycoprotein encoded by a multidrug-resistant gene as revealed by a monoclonal antibody, MRK 16
.
Cancer Res.
48
:
1926
1929
.
Taniuchi
,
I.
,
M.
Osato
,
T.
Egawa
,
M.J.
Sunshine
,
S.C.
Bae
,
T.
Komori
,
Y.
Ito
, and
D.R.
Littman
.
2002
.
Differential requirements for Runx proteins in CD4 repression and epigenetic silencing during T lymphocyte development
.
Cell.
111
:
621
633
.
Tanner
,
S.M.
,
E.M.
Staley
, and
R.G.
Lorenz
.
2013
.
Altered generation of induced regulatory T cells in the FVB.mdr1a-/- mouse model of colitis
.
Mucosal Immunol.
6
:
309
323
.
Thiebaut
,
F.
,
T.
Tsuruo
,
H.
Hamada
,
M.M.
Gottesman
,
I.
Pastan
, and
M.C.
Willingham
.
1987
.
Cellular localization of the multidrug-resistance gene product P-glycoprotein in normal human tissues
.
Proc. Natl. Acad. Sci. USA.
84
:
7735
7738
.
Um
,
J.H.
, and
J.
Yun
.
2017
.
Emerging role of mitophagy in human diseases and physiology
.
BMB Rep.
50
:
299
307
.
Van den Berge
,
K.
,
F.
Perraudeau
,
C.
Soneson
,
M.I.
Love
,
D.
Risso
,
J.P.
Vert
,
M.D.
Robinson
,
S.
Dudoit
, and
L.
Clement
.
2018
.
Observation weights unlock bulk RNA-seq tools for zero inflation and single-cell applications
.
Genome Biol.
19
:
24
.
Wang
,
R.
,
C.P.
Dillon
,
L.Z.
Shi
,
S.
Milasta
,
R.
Carter
,
D.
Finkelstein
,
L.L.
McCormick
,
P.
Fitzgerald
,
H.
Chi
,
J.
Munger
, and
D.R.
Green
.
2011
.
The transcription factor Myc controls metabolic reprogramming upon T lymphocyte activation
.
Immunity.
35
:
871
882
.
Wang
,
D.
,
H.
Diao
,
A.J.
Getzler
,
W.
Rogal
,
M.A.
Frederick
,
J.
Milner
,
B.
Yu
,
S.
Crotty
,
A.W.
Goldrath
, and
M.E.
Pipkin
.
2018
.
The transcription factor Runx3 establishes chromatin accessibility of cis-regulatory landscapes that drive memory cytotoxic T lymphocyte formation
.
Immunity.
48
:
659
674.e6
.
Wolf
,
F.A.
,
P.
Angerer
, and
F.J.
Theis
.
2018
.
SCANPY: large-scale single-cell gene expression data analysis
.
Genome Biol.
19
:
15
.
Yu
,
G.
,
L.G.
Wang
,
Y.
Han
, and
Q.Y.
He
.
2012
.
clusterProfiler: an R package for comparing biological themes among gene clusters
.
OMICS.
16
:
284
287
.
Zhang
,
Y.
,
T.
Liu
,
C.A.
Meyer
,
J.
Eeckhoute
,
D.S.
Johnson
,
B.E.
Bernstein
,
C.
Nusbaum
,
R.M.
Myers
,
M.
Brown
,
W.
Li
, and
X.S.
Liu
.
2008
.
Model-based analysis of ChIP-Seq (MACS)
.
Genome Biol.
9
:
R137
.
Zhou
,
S.F.
2008
.
Structure, function and regulation of P-glycoprotein and its clinical relevance in drug disposition
.
Xenobiotica.
38
:
802
832
.
Zhou
,
S.
,
J.D.
Schuetz
,
K.D.
Bunting
,
A.M.
Colapietro
,
J.
Sampath
,
J.J.
Morris
,
I.
Lagutina
,
G.C.
Grosveld
,
M.
Osawa
,
H.
Nakauchi
, and
B.P.
Sorrentino
.
2001
.
The ABC transporter Bcrp1/ABCG2 is expressed in a wide variety of stem cells and is a molecular determinant of the side-population phenotype
.
Nat. Med.
7
:
1028
1034
.
Zorov
,
D.B.
,
M.
Juhaszova
, and
S.J.
Sollott
.
2014
.
Mitochondrial reactive oxygen species (ROS) and ROS-induced ROS release
.
Physiol. Rev.
94
:
909
950
.

Competing Interests

Disclosures: The authors declare no competing interests exist.

Author notes

*

M.L. Chen and A. Sun contributed equally to this paper.

W. Cao’s present address is Atreca Pharmaceuticals, Redwood City, CA.

This article is distributed under the terms of an Attribution–Noncommercial–Share Alike–No Mirror Sites license for the first six months after the publication date (see http://www.rupress.org/terms/). After six months it is available under a Creative Commons License (Attribution–Noncommercial–Share Alike 4.0 International license, as described at https://creativecommons.org/licenses/by-nc-sa/4.0/).