Radiation therapy induces immunosenescence mediated by p90RSK

Introduction

Radiation therapy (RT) is an important component of managing primary thoracic malignancies, e.g., lung cancer, breast cancer and lymphoma. The delayed effects of RT on the cardiovascular system causes substantial morbidity and mortality in cancer survivors (1). Cardiovascular disease (CVD) is a leading cause of premature morbidity and mortality among breast cancer (2–4), Hodgkin or non-Hodgkin lymphoma (5–8), and lung cancer (9, 10) survivors more than 5 years after their diagnosis and treatment. RT to the chest in Hodgkin lymphoma increases the patients’ risk of CVD by about sevenfold, based on the results of multivariate analyses controlling for other risk factors (11). The incidence of major coronary events increases linearly by 7.4% per Gy in breast cancer survivors (1). Modern RT techniques significantly decrease the doses of radiation to the heart, but in breast cancer patients, the heart may still receive doses of 1–5 Gray (Gy) (1). A more complete understanding of the mechanisms by which RT induces CVD could lead to specific preventive therapeutic approaches.

Risk factors for CVD are generally greater for cancer survivors than the general population, with higher values for mean body mass index, C-reactive protein and low density lipoprotein cholesterol (12). One of the common phenotypes observed in cancer survivors is accelerated senescence (12–14). It is becoming evident that genotoxic chemotherapy agents and radiation induce mitochondrial dysfunction and cellular senescence, i.e., stress-induced premature senescence (SIPS). DNA damaging agents can cause SIPS in a relatively short period (usually 3–10 days) with or without significant telomere shortening (15). Most DNA damage is repaired by DNA damage response (DDR) mechanisms within 24 h after stress (16), but, in contrast, telomeric DNA damage persists for months (17). Therefore, telomeric DNA damage-induced SIPS may be a feasible explanation for the late or delayed CVD effects triggered by cancer treatments.

Although cellular senescence can lead to cell cycle arrest, senescent cells remain metabolically active and secrete a variety of factors, including inflammatory cytokines, chemokines, growth factors, angiogenic factors, nitric oxide (NO), prostaglandin E2 (PGE2), and reactive oxygen species (ROS), constituting a senescence-associated secretory phenotype (SASP). To maintain tissue homeostasis, removing senescent cells promptly is crucial. The SASP cytokines can recruit immune cells to clear the senescent cells by efferocytosis. On the other hand, “immunosenescence” refers to the decline in the immune response with aging (18), including impaired clearance of senescent cells. The critical role of immunosenescence in regulating both tumorigenesis (19) and CVD (20) has been suggested, but the exact regulatory mechanism remains unclear.

As cells age or are exposed to mutagens, they may acquire mutations that lead to clonal expansion without malignant transformation. In the hematopoietic system, this process is known as clonal hematopoiesis of indeterminant potential (CHIP). In over 70% of people with CHIP, somatic mutations of TET2, DNMT3A, ASXL1, and JAK2 are commonly observed, and all but one (JAK2) are loss of function mutations (21, 22). CHIP is associated with a more than twofold increase in risk of CVD (22, 23), which is likely mediated through inflammatory pathways (22, 24). Among survivors of solid malignancies, CHIP-associated mutations are common and are detectable in around 30% of individuals in their 7th decade of age, and may adversely impact overall survival (25). Oncologic therapies, such as RT, are frequently dose-limited based on toxicity to normal organs—organs that may already have impaired function in the presence of CHIP-associated mutations (22, 26–28). It is undefined whether the presence of CHIP in patients undergoing RT for thoracic malignancies contributes to an increased risk of radiotherapy-induced CVD.

The critical role of NAD⁺ depletion during aging in regulating metabolic dysfunction is well established (29). Previously, we reported that NAD⁺ was depleted after IR via PARP activation, triggering persistent PISP in myeloid cells (30). In addition to PARP activation, CD38 may be involved in NAD⁺ depletion during aging (31). CD38 is a NAD glycohydrolase, and the levels and activity of CD38 were significantly increased by aging, leading to NAD⁺ depletion (31). CD38 is expressed mainly in immune cells, and recently, Chini et al. have reported that senescence-induced inflammation promoted CD38 accumulation and caused NAD⁺ depletion (32). Taken together, these reports suggested that both PARP activation and CD38 expression promote NAD⁺ depletion during the aging process.

T-bet (official gene name: TBX21) is an immune cell-specific member of the T-box family of transcription factors. In B cells, antigen stimulation activates B cell proliferation and induces affinity maturation and class switch recombination to produce high-affinity antibodies, which show different immune effector functions (33). T-bet induces Iγ2A transcripts and promotes IFNγ-mediated class-switching to the IgG2a isotype, which plays a critical role in the pathogenesis of humoral autoimmunity and protection against pathogens (34). Furthermore, T-bet promotes the survival of antigen-specific IgG2a⁺/CD38^hi memory B cells by upregulating mature B cell receptor transcription (35). Also, T-bet drives the migration of memory IgG2a⁺ B cell to sites of inflammation by increasing CXCR3 expression (36). Importantly, the age-associated B cell (ABC) is characterized by T-bet expression and increases overall inflammatory state (37). Although systemic T-bet depletion inhibited atherosclerosis formation (38), the role of T-bet in B cells on atherosclerosis formation remains unclear.

Previously, we reported the role of p90RSK activation in SASP after RT. IR increases p90RSK activation in myeloid cells and forms a nucleus-mitochondria positive feedback loop with p90RSK-mediated ERK5 S496 phosphorylation, leading to persistent mitochondrial (mt) ROS production and subsequent telomere DNA damage. We found that poly (ADP-ribose) polymerase activation induced by telomere DNA damage promoted nicotinamide mononucleotide (NAD⁺) depletion and consequent mitochondrial dysfunction, which elicited mtROS production, activated redox-sensitive kinase of p90RSK, and formed a positive feedback loop. Therefore, p90RSK activation played a crucial role in inducing SASP.

Although RT may induce immunosenescence (39), the mechanisms remain incompletely defined. Therefore, in this study, we use multiparameter mass cytometry (CyTOF) to explore the association of RT with various senescence-related factors and determinants of senescence including DDR -related markers, efferocytosis, and CHIP-associated drivers in PBMCs from thoracic cancer patients before and 3 months after RT.

Materials and methods

Patients selection

The study was approved by the Institutional Review Board of the University of Texas MD Anderson Cancer Center (#PA16-0971). The inclusion criteria of this study were as follows: (1) ≥ 18 years of age, (2) Karnofsky Performance Scale ≥ 70, (3) To receive radiation therapy (RT) with computed tomography (CT)-based treatment planning that will involve delivery of dose to the heart, (3) An estimated mean heart dose of at least 6 Gy and mean left ventricular dose of at least 5 Gy, as estimated by the treating radiation oncologist at the time of simulation.

The exclusion criteria of this study were as follows: (1) Unable or unwilling to give written informed consent, (2) Previous history of RT to the thorax or breast, (3) Allergy to gadolinium, (4) Implanted device that is non-MRI compatible, including cardiac devices and neurostimulators, (5) Pregnant or breast-feeding, (6) Atrial fibrillation or frequent ventricular or atrial premature beats, (7) GFR < 60 mL/min according to the Modification of Diet in Renal Disease equation, (8) Personal history of coronary artery disease or myocardial disease, (9) Personal history of hypertension, requiring > 1 antihypertensive agent to maintain blood pressure < 140/90, (10) Valvular stenosis or regurgitation of > moderate severity, (11) Heart failure (baseline NYHA > 2), (12) Systolic BP < 90 mmHg (15) Pulse < 50/min, (16) History of pulmonary hypertension or elevated right ventricular systolic pressures by echocardiogram, (17) Renal failure necessitating dialysis, (18) Suspicion or diagnosis of amyloidosis, (19) Suspicion or diagnosis of hemochromatosis, (20) Unable to obtain an MRI scan with gadolinium for any other reason that is not listed above.

Panel design and metal conjugation of antibodies

We generated a CyTOF panel which includes antibodies against CHIP drivers, SASP-related proteins, and various cell surface markers for PBMCs such as the markers described in the previous report (40). We used 37 antibodies, Ir DNA-Intercalator (Cell-ID Intercalator-Ir; Fluidigm #201192A) for nucleus staining and Rh DNA-intercalator (Cell-ID Intercalator-103Rh; Fluidigm #201103A) for viability staining, as listed in Supplementary Table 1. We designed the panel using Maxpar Panel Designer software (Fluidigm). The antibodies listed in Supplementary Table 1, except for the metal-conjugated antibodies that were purchased from Fluidigm, BioLegend, BD Biosciences, and Miltenyi, were conjugated to lanthanides using the Maxpar X8 Multimetal Labeling Kit (Fluidigm) or cadmium using the Maxpar MCP9 Antibody Labeling Kit (Fluidigm) according to the manufacturer’s protocol. To determine optimal antibody concentrations, three different concentrations of each metal-labeled antibody were tested by CyTOF before being used in analyses.

The 90 kDa ribosomal S6 kinases (p90RSKs) are a group of serine/threonine kinases consisting of 4 RSK isoforms (RSK1-4). The expressions of p90RSK1-4 isoforms in human are ubiquitous, and all of them express in immune cells.¹ The p-p90RSK antibody recognizes the phosphorylation of p90RSK1 S380 site, which has a homology between p90RSK1-3 at LFRGFSpFVA, and p90RSK4 at SpFVA. Therefore, at least this antibody can recognize p90RSK1-3 isoforms.

Cytometry sample preparation and staining

Whole blood was collected at two different time points, before (baseline) and 3 months after RT, from patients (n = 16) with thoracic cancer treated at The University of Texas MD Anderson Cancer Center (Supplementary Table 2). PBMCs were isolated using Ficoll-Paque Plus (Fisher #45-001-750) according to the manufacturer’s protocol and were cryopreserved in liquid nitrogen until analysis. Cryopreserved PBMC samples were thawed in RPMI containing 5% pooled human serum (PHS), 10 mM HEPES, 0.2% gentamicin, and 2.5 μ/mL benzonase nuclease, washed once with the same medium and repelleted by centrifugation at 300 g at 4^°C). One million cells were applied to the viability staining with Cell-ID Intercalator-103Rh and washed with Maxpar Cell staining buffer (Fluidigm #201068). For cell fixation, the cells were resuspended in 1x Fix I Buffer (Fluidigm #201065) and start the barcoding step using Cell-ID 20 Plex Pd Barcoding Kit (Fluidigm #201060). All barcoded samples were combined, and the staining procedure was performed according to protocol, i.e., Maxpar Phosphoprotein Staining with Fresh Fix (Fluidigm PN400278A4 protocol). The cells were blocked by Human TruStain FcX (BioLegend #422302), followed by cell surface staining for cell-type identification. After the wash with Maxpar Cell staining buffer, the cells were chilled on ice and incubated at 4^oC with methanol for permeabilization, followed by two washes with Maxpar Cell staining buffer. The cells were incubated with metal-labeled antibodies against cytoplasmic proteins followed by two washes. For the final fixation, the cells were incubated with 1.6% formaldehyde, followed by the incubation with Cell-ID Intercalator-Ir for nucleus staining. After washing twice with Maxpar Cell Staining Buffer and once with Maxpar water (Fluidigm #201069), the cells were filtered to make a single cell suspension and centrifuged. After adding EQ beads (Fluidigm #201078), CyTOF data were acquired on a Helios mass cytometer (Fluidigm).

Data acquisition for cytometry

CyTOF data were acquired on a Helios mass cytometer (Fluidigm) and were analyzed using the Astrolabe Cytometry Platform (Astrolabe Diagnostics, Inc.). Single-cell data have been clustered using the FlowSOM R package (41) and labeled using the Ek’Balam algorithm (42). Cell subset definitions were followed as previously described (43, 44). Cluster labeling, method implementation, and visualization were done through the Astrolabe Cytometry Platform (Astrolabe Diagnostics, Inc.).

Irradiation of human peripheral blood mononuclear cells followed by flow cytometry

Human peripheral blood mononuclear cells (hPBMC) were seeded into a 48-well plate in culture medium (RPMI 1640 with 10%FBS, 10 mM HEPES, 1 × non-essential amino acid mix, 1 mM sodium pyruvate, 100 mg/mL streptomycin, 0.55 mM 2-mercaptoethanol) at a density of 2 million cells per well. Treatments were performed in triplicate for each group. We pre-treated the cells with FMK (RSK inhibitor Fmk, Soluble in DMSO, Axon Medchem Cat#1848) or DMSO for 1 h before irradiation. The cells were irradiated at 2 Gy while the non-irradiation control was set aside. Twenty-four hours after irradiation, we harvested the cells of each group and analyzed them on an LSRII instrument (Becton Dickinson). Flow CyTOF data were analyzed using FlowJo software. The following Cell surface markers were used: human CD3 (Biotin, clone: UCHT1, eBioscience), CD14 (Per-Cy5.5, clone: 63D3, BioLegend), CD19 (PE-Cy7, clone: HIB19, BioLegend), CD27 (PE, clone: O323, Tonbo), IgD (FITC, clone: MEM-59, BioLegend), CD38 (BUV395, clone: HB7, BD), and T-bet (Alexa Fluor 647, clone:4B10, Biolegend). Mouse IgG1 Isotype Control APC (clone:MOPC-21, Tonbo) was used as a negative control.

Mitochondria-specific reactive oxygen species measurements

PBMC were washed with PBS and resuspended in RPMI 1640 supplemented with 10% FBS and 1% each of L-glutamine, sodium pyruvate, non-essential amino acids, 100 mg/ml streptomycin, 0.1 M HEPES, 55 μM β-mercaptoethanol. PBMC were seeded at a density of 2 × 10⁵/ml and cultured 6 h at 37^°C and 5%CO2. After resting, cells were pretreated with and without FMK-MEA (10 μM)/1 h and later treated with IR (2Gy) for 24 h. After that, cell were washed with a medium two times and incubated with Mito sox (5 μM) (Molecular Probes, Invitrogen, Carlsbad, CA) for 20 min at 37°C. Then cells were incubated with anti-CD19 APC and anti-CD14 PerCp Cy5.5 at 4°C for 30 min. Autofluorescence of unstained cells was used as a control for each sample. Approximately 100,000 gated events were acquired for each sample on BD LSR II (Becton-Dickinson, San Jose, CA). Dead cells and debris were excluded based on forward scatter and side scatter measurements. The percentage of MitoSOX expression positive cells were quantified on CD19 APC -positive gated Cells or CD14 PerCp Cy5.5-positive gated cells in the PE channel analyzed using FlowJo software (TreeStar, Ashland, OR).

Senescence associated β-galactosidase activity in cell cultures

After 24 h of ionizing radiation, PBMCs were incubated with 33 mM C12FDG (Cayman chemicals) for 30 min in a 37^°C water bath. Cells were centrifuged at 500 g for 5 min, washed with Cell Staining buffer (BD Biosciences), and stained with anti-CD19-APC and anti-CD14-PerCPCy5.5 at 4^°C for 30 min. Cells were then washed in Cell Staining buffer and resuspended in Cell Staining buffer and analyzed by a BD LSR II Flow cytometer. Gates were set as follows: Dead cells and debris were excluded based on forward scatter and side scatter measurements, B cells (CD19⁺), Monocytes (CD14⁺) and analyzed for C 12 FDG fluorescence. Data were analyzed with FlowJo software (v.10).

Statistical analysis

The multidimensional scaling (MDS) map was generated using the cmdscale R package (45). Differential abundance analysis was done using the edgeR R package (46, 47) as previously described (48). Differential expression analysis was done using the limma R package (49), as described previously (50). The frequency data within –log(FDR) > 0.5 and log [fold change (FC)] < –1 or 1 < log(FC) were considered as significant. The expression data within adjusted P-value (adj. P-Val) < 0.05 were considered as significant. For the adj. P-Val, limma R package is using the Benjamini-Hochberg correction.

We determined differences between 2 independent groups using the Student’s t-test (2-tailed) and, when applicable, 1-way analysis of variance followed by Bonferroni post-hoc testing for inter- group comparisons using GraphPad Prism (GraphPad Software, San Diego, CA, USA). We used Welch’s analysis for evaluating the equal variance of each group. P-values < 0.05 were considered statistically significant.

Results

The responses of immune cells after radiation treatment in cancer patients

Previous studies showed that ionizing radiation and various cancer treatments induced a pro-inflammatory senescent phenotype in myeloid cells (30). To characterize how RT can change immune cell phenotypes in cancer patients using multiparameter mass CyTOF, we enrolled 16 patients with esophagus or lung cancer who were about to start RT (Supplementary Table 2). Blood was collected at the time of enrollment and 3 months after the initial radiation treatment, and peripheral blood mononuclear cells (PBMCs) were collected. We then evaluated the phenotypic changes of immune cells by multiparameter mass CyTOF with general cell surface markers and various senescent markers as described in the methods. Analysis using the Astrolabe platform (40) automatically labeled canonical immune cell subsets (Figure 1A and Supplementary Table 2). Astrolabe identified eleven T cell subsets [including CD4⁺ and CD8⁺ T cells, naïve T_EMRA (effector memory cells re-expressing CD45RA), effector memory, and central memory cells], three B cell subsets, five myeloid cell subsets, two NK cell subsets, and granulocytes. We found that only B cells significantly decreased after RT (Figure 1B). However, when we looked at each immune cell subset, we found four immune cell subsets: CD4⁺ T_EMRA cell, B cell (Memory), T cell (unassigned), and naïve B cell (CD27^–) were significantly decreased after RT (Figure 1C). Our data suggest that RT may play a significant role in regulating B cell functions in both memory and naïve B cells, including antigen producing and presentation. Furthermore, CD4⁺ T_EMRA cells can show a cytotoxic phenotype with upregulation of CX3CR1 expression (51–54), which is related to CD4 and CD8 T cell cytotoxic potentials. Therefore, these data suggest the potential role of RT in regulating cytotoxic abilities of T cells.

FIGURE 1

Figure 1. Cell subset using the Astrolabe platform. The FlowSOM algorithm clustered cells and labeled subsets according to the canonical gating hierarchy. (A) Cell subset heat map. The rows represent cell subsets, and the columns represent cell protein expression. Tile intensity showed the median value in the subset according to the color key. (B,C) Multidimensional scaling (MDS) map for each cell type (B) and cell subsets clustered by Astrolabe platform as shown in (A). Each bubble represents a cell type (B) or subset (C), and the size of that bubble is determined by the median frequency of cells contained in that bubble across all samples analyzed. Each bubble is colored based on the magnitude of the fold change post-radiation therapy (RT) over baseline.

Additional subsets revealed using functional antibodies

We further included the staining data with functional antibodies, and the Astrolabe Platform derived 138 functional profiling subsets with unsupervised clustering by FlowSOM algorithm, which were labeled by the markers that provided the most significant separation among them (Figure 2 and Supplementary Figure 1). We found that the frequencies of 5 functional profiling subsets were significantly upregulated, and seven functional profiling subsets were downregulated after RT (Figures 2B,C, 3A and Supplementary Figure 2). Among 12 functional profiling subsets, we found that only three subsets, CD27^–/Ki67^lo/CD38^lo/DNMT3a^hi B cells, CD4⁺/CD25^lo/IL-7Ra^lo/TET2^lo T_Naive cells and CD4⁺/Il-7Ra^lo/PD-1^lo/CD25^lo T_EMRA cells, showed > 1% frequency among total cells and were decreased after RT. These data suggest that the effects of RT on immune cells were relatively specific to a few unique functional subsets of naïve B cell and CD4⁺ T cell at 3 months after RT.

FIGURE 2

Figure 2. Profiling subsets using the Astrolabe platform. (A) Profiling subset heat map for B cells only. Refer to Supplementary Figure 1 for the complete heat map. Rows are cell subsets; columns are channels. Tile intensity is the median value measured in each channel in each profiling. (B) Profiling of every cell cluster (138 clusters) was visualized using MDS. Each cluster is colored based on the logarithm of the ratio of intensity in the post-RT sample to pre-RT (baseline) sample. (C) A volcano plot of identified cell clusters profiled in (B), showing those with –log10 [the false discovery rate (FDR)] > 0.5, log(FC) < –1 or 1 < log(FC) and p < 0.05 (student’s t-test) after RT (pre-RT vs. post-RT as described in Figure 4A and Supplementary Figure 2) in light green. Once the frequency data within –log(FDR) > 0.5 and log[fold change (FC)] < –1 or 1 < log(FC) were considered as significant, 31 profiling subsets were changed by RT. Then 12 profiling subsets were picked out from the 31 subsets by student’s t-test.

FIGURE 3

Figure 3. B cell profiling subset and T-bet expression. (A) The frequency of each B cell profiling subset for all the cells in an experiment at pre-RT and post-RT. *p < 0.05 (or **p < 0.01) by Student’s t-test, –log(FDR) > 0.5 and log(FC) < –1 or 1 < log(FC). (B) T-bet expression in each B cell profiling subset. The significant increase of T-bet expression after RT was observed only in the Ki67^–DNMT3a⁺ naïve B cell subset. Mean ± SD, *adj. P-Val < 0.05.

FIGURE 4

Figure 4. The relationship between T-bet and p90RSK in Ki67^–DNMT3a⁺ naïve B cell subset. (A) p-p90RSK expression in each cell types. Mean ± SD, *p < 0.05, **p < 0.01, ***p < 0.001. Statistical significance was assessed by Dunn’s multiple comparisons test. (B) p-p90RSK expression in Ki67^–DNMT3a⁺ naïve B cell subset showed no difference between pre-RT and post-RT. (C) The relationship between p-p90RSK and T-bet expression in Ki67^–DNMT3a⁺ naïve B cell subset at pre-RT and post-RT. (D) The relationship between total p90RSK and T-bet expression in Ki67^–DNMT3a⁺ naïve B cell subset at pre-RT and post-RT. Statistical significance and correlation were assessed by Pearson correlation.

Radiation therapy uniquely induced T-bet and CD38 expression in Ki67^–DNMT3a⁺ naïve B cell subset in cancer patients

Astrolabe revealed five functional subsets of B cells (CD27^–). Among all the B cell subsets, CD27^–/Ki67^lo/CD38^lo/DNMT3a^hi (Ki67^–DNMT3a⁺ naïve B cell) was the largest, and CD27^–/Ki67^lo/CD38^lo/DNMT3a^lo was the second-largest subset (Figure 3A). Although B cells did not highly express T-bet compared to other cell types in the baseline (Supplementary Figure 3) and the T-bet expression in Ki67^–DNMT3a⁺ naïve B cell subset was not exceptionally high compared to different B cell subsets in the pre-RT (baseline) samples (Figure 3B and Table 1), we found that a significant increase of T-bet expression after RT was observed only in Ki67^–DNMT3a⁺ naïve B cell subset (Figure 3B). Since it has been reported that most of the age- ABCs express T-bet (55), and we previously reported the role of p90RSK activation in regulating senescence (30, 56), we hypothesized that p90RSK activation could up-regulate T-bet expression. First, we investigated the relationship between the intensities of p90RSK S380 phosphorylation (p-p90RSK) and T-bet expression. p-p90RSK expression in B cells showed no significant difference with other cell types except NK cells (Figure 4A), and p-p90RSK expression did not increase after RT in the Ki67^–DNMT3a⁺ naïve B cell subset (Figure 4B). But interestingly, we found a good correlation between p-p90RSK and T-bet in the samples from pre-RT, post-RT, and the combination of both (Figure 4C). Next, we compared T-bet expression with the intensities of total p90RSK (t-p90RSK). Although we could not find any significant correlation between T-bet and t-p90RSK in the samples from pre-RT, we found a significant correlation between T-bet and t-p90RSK in post-RT samples and a combination of both pre-and post-RT (Figure 4D). These data suggest that p90RSK activation and total p90RSK expression may play a role in T-bet expression.

TABLE 1

Table 1. Profiling subsets, which showed a significant difference in indicated protein expression level between pre- and post-RT, using the Astrolabe Platform.

Astrolabe identified all protein expression levels examined in this study in each cell type (Figure 5A) and cell subtype (Figure 5B) and found that TOP2b and CD38 expression were increased after RT in the total B cells after RT (Figure 5A), and in the subsets of CD27^– B cells (Figure 5B). A significant increase in CD38 expression was also observed in CD8⁺ T_EMRA cells (Figure 5B). These data suggest a contribution of T-bet and CD38 expression to the induction of senescent phenotype after RT. Since TOP2b can protect against DNA damage, the increase of TOP2b expression after RT may contribute to DNA damage repair after RT.

FIGURE 5

Figure 5. Only TOP2b and CD38 expression were significantly increased after RT at the cell subset level. (A) Only TOP2b and CD38 expression in B cells were increased after RT, and no significant change of studied molecules in cell types after RT compared to pre-RT level was observed. (B) CD38 and TOP2b expression in B cell (CD27^–) and CD8⁺ T_EMRA subset pre- and post-RT. Mean ± SD, *adj. P-Val < 0.05, **adj. P-Val < 0.01, ***adj. P-Val < 0.001.

Radiation-induced T-bet⁺ memory and naïve B cells and myeloid cells by 90 kDa ribosomal S6 kinase activation

We next investigated whether p90RSK activation played a role in ionizing radiation (IR)-induced increase in CD38⁺T-bet⁺ cells. PBMCs from the healthy human subject were pre-treated with vehicle or FMK-MEA (a specific inhibitor of p90RSK), then exposed to IR (2 Gy). Previously, we have confirmed the specificity of FMK-MEA by using the Ambit/DiscoverX platform to screen 443 kinases (57). Twenty-four hours following the irradiation, the numbers and frequencies of naïve and memory B cells were decreased (Figures 6A,B). However, we found a significant increase of T-bet⁺ expression in memory and naïve B cells and this effect was inhibited by the p90RSK inhibitor FMK-MEA (Figures 6C–F). Since we have reported the role of SASP-induced CD38 expression in myeloid cells (32), we also investigated the impact of IR on CD38^+/T-bet⁺ myeloid cells, and we found that pre-treatment of FMK-MEA also significantly inhibited IR-mediated CD38⁺/T-bet⁺ myeloid cell induction (Figure 7). These data suggest that IR increases CD38⁺/T-bet⁺ cells induction in both B-cells and myeloid cells even in the early phase after IR, dependent on p90RSK activation.

FIGURE 6

Figure 6. IR-induced CD38 + T-bet + CD19 + B cell via p90RSK activation. (A) Memory and naive B cell gating among CD19⁺ B cells. (B) Frequency of memory B cells (Left) and naive B cells (Right) among CD19⁺ cells 24 h after irradiation/non-irradiation with/without FMK-MEA. (C) Representative FACS plots for T-bet⁺ memory cells and their frequency among CD19⁺ B cells (D) 24 h after irradiation/non-irradiation with/without FMK. (E) Representative FACS plots for T-bet⁺ naive B-cells and their frequency among CD19⁺ B cells (F) 24 h after irradiation/non-irradiation with/without FMK. N = 3 for each group, and the experiments were done more than twice. Mean ± SD, *p < 0.05, ^**p < 0.01, and ^****p < 0.001.

FIGURE 7

Figure 7. IR-induced CD38⁺/T-bet⁺ myeloid cells via p90RSK activation. (A) Representative FACS plots of myeloid T-cell distribution are depicted. (B) Myeloid cell frequency 24 h after irradiation/non-irradiation with/without FMK-MEA (10 μM). (C,D) Representative FACS plots (C) and their frequency (D) of T-bet expression in CD14⁺ myeloid cells in each condition. N = 3 for each group, and the experiments were done more than twice. Mean ± SD, *p < 0.05 and ^**p < 0.01. (E) The ratio of cell numbers [myeloid cells (Blue) and lymphocytes (Red)] at pre-RT and post-RT. Mean ± SD, ^****p < 0.001. (F,G) PBMCs were exposed to IR (2 Gy) with pre-treatment of FMK-MEA (10 μM) or vehicle for 1 hr. Senescent B cells and monocytes from PBMCs were measured using the C12FDG, a fluorogenic substrate for SA-β-galactosidase. Data shows the Quantification of positive C12FDG⁺ cells from CD19⁺ (Left) CD14⁺ (Right) cells. (G) Mitochondrial superoxide production was measured using the Mitosox dye in CD19⁺ (left) CD14⁺ (right) cells. Mean ± SD, n = 3, ^**p < 0.01.

Immunosenescence after radiation therapy and 90 kDa ribosomal S6 kinase activation

As shown in Supplementary Table 2 and Figure 7E, we found a significant decrease in lymphocytes number after RT. In contrast, myeloid cells numbers were relatively unaffected by RT, which resulted in relative myeloid skewing in the RT patients. To determine if RT induces immunosenescent in immune cells, we stained senescence-associated β-gal (SA-β-gal) staining and mitochondrial ROS (mtROS) production using healthy human PBMCs. Interestingly, CD19⁺ B cells and CD14⁺ myeloid cells in the PB showed significant increase of SA-β-gal⁺ cells and mtROS production after IR (Figures 7F,G). The increase of SA-β-gal⁺ cells and mtROS production was clearly blocked by pre-treatment of p90RSK inhibitor, FMK-MEA, in both B and myeloid cells. These results also suggest the crucial role of IR-induced p90RSK activation in immunosenescence.

Discussion

In this study, we investigated the effects of RT on various senescence markers, DDR, efferocytosis, and clonal hematopoiesis drivers in immune cells. At 3 months after RT we observed a significant reduction in several immune cell subsets, particularly B cells.

Ki67^–DNMT3a⁺ naïve B cells were the largest B cell subset. Although the frequency of this B cell (CD27^–) subset was reduced after RT, we found a significant increase of T-bet expression that was correlated with phosphorylation of p90RSK in Ki67^–DNMT3a⁺ naïve B cells. Furthermore, we also observed an increase in CD38 expression in B cells after RT. In vitro evidence suggested that p90RSK activation played a crucial role in the induction of CD38 and T-bet expression by IR in memory, naïve B, and myeloid cells. These data link RT-induced activation of p90RSK to upregulation of CD38⁺ and T-bet⁺ which are reported to have a role in immunosenescence (31, 32, 58).

The analysis using Astrolabe Platform of the CyTOF data revealed that the Ki67^–DNMT3a⁺ naïve B cell subset was the largest subset of B cells, and had a unique response to RT. CD27 was used as a marker of memory B cells in humans based on the correlation of CD27 expression with somatic hypermutation in IgM⁺/IgD⁺ cells, and CD27 was constitutively expressed in approximately 40% of peripheral blood B cells in humans (59). CD27 is a member of the TNF-receptor family with the ligand of CD70, which is expressed on the surface of activated T cells. Activation of CD27 signaling plays a crucial role in maintaining long term immunological memory against T cell development antigen (60) by activating B cell expansion, differentiation, and antibody production (61–63). Therefore, the CD27^– B cell is recognized as a naïve B cell.

We examined various senescence and aging gene markers in various blood subsets using CyTOF, hoping to see the correlation between the irradiation and some changes of theses markers. The definition of the various subsets of T- and B- cells are listed in Supplementary Table 3. We focused B cells because only B cells showed significant change of some of these markers (CD38, TOP2β). In addition, the roles of DNMT3A in B cells have been reported. For example, DNMT3 regulates B cell lineage specification in the mouse model (64). Importantly, the depletion of DNMT3a in the B cell lineage results in the leukemic transformation of B1a cells in the CLL mouse model (65), and CD38 is a marker of human CLL B cell activation (66). Therefore, we assume that the B cell subset represented by lower expression of CD38 and higher expression of DNMT3a[i.e., (CD27^–) Ki67^lo CD38^lo DNMT3a^hi cells] are a less proliferative B cell subset in accordance with low Ki67 expression. This was the dominant B cell subset before RT in thoracic cancer patients (Figure 3A). After RT, CD38 expression was increased regardless of Dnmt3a expression, suggesting that the residual B cells after RT were activated.

ABCs are antigen-experienced B cells that are characterized by a T-bet-mediated transcriptional program (67). ABCs plays a critical role in regulating immune response to infectious agents mediated by IgG2a/c and inflammatory cytokines (68, 69). However, the excessive accumulation of ABCs can trigger auto-immune reactions in response to self-antigen and may cause detrimental effects by inducing inordinate inflammation (70). In our current report, we found that T-bet expression was increased in the Ki67^–DNMT3a⁺ naïve B cell subset, suggesting RT differentiated naïve B cells to ABC and induced inflammatory and senescence pathways. Also, the increase of T-bet expression in this largest but less active B cell subset (Ki67^–DNMT3a⁺ naïve B cells) may significantly enhance the pre-mature aging process after RT by differentiating naïve B cells to ABC. Since systemic T-bet depletion inhibited atherosclerosis formation (38), it is possible that the increase of T-bet in B cells plays a role in elevating cardiovascular events after RT. Further investigation will be necessary.

It has been reported that CD38 expression in immune cells is increased in response to infection and inflammation, which subsequently induce inflammatory response (71). In rheumatoid arthritis (RA), the inhibition of CD38 activation by anti-CD38 monoclonal antibody has been effective to reduce RA symptoms and disease progression. In systemic lupus erythematosus (SLE) patients, the upregulation of CD38 expression in marginal zone-like IgD + CD27 + B cells has been reported (72), and the depletion of CD38 in pristane-induced lupus mouse model significantly improve the symptom (73). Thus, CD38 upregulation in B cells seems to be related to inflammatory status. In our current report, we found an increase in CD38 expression in B cells after RT in cancer patients. Since the essential role of CD38 in B cells in inflammatory responses has been reported (68, 69), these data suggest that NAD⁺ depletion induced by CD38 expression in B cells may initiate SASP. The induction of SASP in B cells may contribute to CVD in cancer patients after RT.

Lastly, we found a correlation between T-bet and p-p90RSK in the Ki67^–DNMT3a⁺ naïve B cell subset. Furthermore, we found a crucial role of p90RSK activation in IR-induced CD38⁺/T-bet⁺ memory and naïve B cells and myeloid cells in vitro. Taken together, these data suggested that p90RSK activation induced by RT was involved in immunosenescence via upregulating CD38 and T-bet expression. Notably, the contribution of p90RSK activation and T-bet expression in immune cells to the process of atherosclerosis formation has been reported (38, 56), and to the best of our knowledge, this is the 1st report to show that p90RSK activation is required for T-bet induction. Interestingly, the potential role of T-bet on CD38 induction has been reported (74). Taken together, these data suggest that the activation of p90RSK induced by radiation plays a crucial role in T-bet-CD38 induction, resulting in immunosenescence and consequent atherosclerosis formation. Further investigation will be necessary.

Data availability statement

The original contributions presented in this study are included in the article/Supplementary material, further inquiries can be directed to the corresponding author/s.

Ethics statement

The Institutional Review Board of the University of Texas MD Anderson Cancer Center approved the clinical study protocol (#PA16-0971). The patients/participants provided their written informed consent to participate in this study.

Author contributions

MI, HC, and SK performed experiments, analyzed data, and drafted manuscript. SK, AD, N-TL, EC, KK, VS, and L-LL supported the experiments and interpretation of the data. JH, XX, CR-G, S-CY, KS, SY, ZL, RN, E-aA., JB, NP, and JC contributed to the interpretation of the data and edited the manuscript. SL, MK, MY, and J-iA planned and generated the study design, obtained funding, interpreted data, and wrote the manuscript. All authors contributed to the article and approved the submitted version.

Funding

This work was partially supported by the National Institutes of Health (NIH) to J-iA (grants HL-149303 and AI-156921), and JC (grants HL148338, HL133254, and HL157790), N-TL (grants HL-134740 and HL-149303), RN (grants HL141966 and HL143520), and MY (grant AI121197), the Cancer Prevention and Research Institute of Texas (CPRIT) to J-iA and KS (grant RP190256), RN (grant RP190326), and JC (grant RP200619), the Ted Nash Long Life Foundation to MY, and from Institutional Research Grant to RN. This research performed in the Flow Cytometry and Cellular Imaging Facility (NC) was supported in part by the National Institutes of Health through MD Anderson’s Support (grant CA016672).

Acknowledgments

We thank Scientific Publications, Research Medical Library at The University of Texas MD Anderson Cancer Center for editing and Carolyn J. Giancursio for her technical assistance.

Conflict of interest

Author E-aA was employed by the company Astrolabe Diagnostics, Inc. SL was an Advisory Board member of AstraZeneca, Beyond Spring Pharmaceuticals, STCube Pharmaceuticals.

The remaining authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Publisher’s note

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.

Supplementary material

The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fcvm.2022.988713/full#supplementary-material

Supplementary Figure 1 | The complete profiling subset heat map. Rows are cell profiling; columns are channels. Tile intensity is the median of the channel in each profiling.

Supplementary Figure 2 | Profiling subsets of T cell, NKT cell, and monocyte revealed a significant difference in frequency between pre-RT and post-RT. (A) CD8⁺ T cell profiling subsets showed a significant difference between pre-RT and post-RT. (B) CD4⁺ T cell profiling subsets showed a significant difference between pre-RT and post-RT. (C) NKT cell profiling subset showed a significant difference between pre-RT and post-RT. (D) The monocyte profiling subset showed a significant difference between pre-RT and post-RT. Mean ± SD, *p < 0.05 (or ^**p < 0.01) by Student’s t-test, –log(FDR) > 0.5 and log(FC) < –1 or 1 < log(FC).

Supplementary Figure 3 | T-bet expression in each cell type at pre-RT. Mean ± SD, ^***p < 0.001. Statistical significance was assessed by Dunn’s multiple comparisons test.

Supplementary Table 1 | CyTOF antibody panel and the markers with each protein function. To establish our CyTOF antibody panel, some cell surface markers were picked from the previous report (40).

Supplementary Table 2 | Patient profiles. NSCLC, non-small cell lung cancer; SCLC, small cell lung cancer.

Supplementary Table 3 | The definitions of each cell subset labeled by the Astrolabe platform.

Supplementary File 1 | Expression level of each protein in each cell type at Pre-RT and Post-RT in samples obtained from 16 patients.

Supplementary File 2 | Expression level of each protein in each cell subtype at Pre-RT and Post-RT in samples obtained from 16 patients.

Supplementary File 3 | Expression level of each protein in each profiling subset at Pre-RT and Post-RT in samples obtained from 16 patients.

Supplementary File 4 | Frequency of each profiling subset at Pre-RT and Post-RT in samples obtained from 16 patients.

Footnotes

1. ^ https://www.nextprot.org

References

1. Darby SC, Ewertz M, McGale P, Bennet AM, Blom-Goldman U, Bronnum D, et al. Risk of ischemic heart disease in women after radiotherapy for breast cancer. N Engl J Med. (2013) 368:987–98. doi: 10.1056/NEJMoa1209825

PubMed Abstract | CrossRef Full Text | Google Scholar

2. Hooning MJ, Botma A, Aleman BM, Baaijens MH, Bartelink H, Klijn JG, et al. Long-term risk of cardiovascular disease in 10-year survivors of breast cancer. J Natl Cancer Inst. (2007) 99:365–75. doi: 10.1093/jnci/djk064

PubMed Abstract | CrossRef Full Text | Google Scholar

3. Cheng YJ, Nie XY, Ji CC, Lin XX, Liu LJ, Chen XM, et al. Long-term cardiovascular risk after radiotherapy in women with breast cancer. J Am Heart Assoc. (2017) 6:e005633. doi: 10.1161/JAHA.117.005633

PubMed Abstract | CrossRef Full Text | Google Scholar

4. Zagar TM, Cardinale DM, Marks LB. Breast cancer therapy-associated cardiovascular disease. Nat Rev Clin Oncol. (2016) 13:172–84. doi: 10.1038/nrclinonc.2015.171

PubMed Abstract | CrossRef Full Text | Google Scholar

5. Swerdlow AJ, Higgins CD, Smith P, Cunningham D, Hancock BW, Horwich A, et al. Myocardial infarction mortality risk after treatment for Hodgkin disease: A collaborative British cohort study. J Natl Cancer Inst. (2007) 99:206–14. doi: 10.1093/jnci/djk029

PubMed Abstract | CrossRef Full Text | Google Scholar

6. Singal PK, Iliskovic N. Doxorubicin-induced cardiomyopathy. N Engl J Med. (1998) 339:900–5. doi: 10.1056/NEJM199809243391307

PubMed Abstract | CrossRef Full Text | Google Scholar

7. Maraldo MV, Giusti F, Vogelius IR, Lundemann M, van der Kaaij MA, Ramadan S, et al. Cardiovascular disease after treatment for Hodgkin’s lymphoma: An analysis of nine collaborative EORTC-LYSA trials. Lancet Haematol. (2015) 2:e492–502. doi: 10.1016/S2352-3026(15)00153-2

CrossRef Full Text | Google Scholar

8. Aleman BM, van den Belt-Dusebout AW, De Bruin ML, van ’t Veer MB, Baaijens MH, de Boer JP, et al. Late cardiotoxicity after treatment for Hodgkin lymphoma. Blood. (2007) 109:1878–86. doi: 10.1182/blood-2006-07-034405

PubMed Abstract | CrossRef Full Text | Google Scholar

9. Armenian SH, Xu L, Ky B, Sun C, Farol LT, Pal SK, et al. Cardiovascular disease among survivors of adult-onset cancer: A community-based retrospective cohort study. J Clin Oncol. (2016) 34:1122–30. doi: 10.1200/JCO.2015.64.0409

PubMed Abstract | CrossRef Full Text | Google Scholar

10. Lally BE, Detterbeck FC, Geiger AM, Thomas CR Jr., Machtay M, Miller AA, et al. The risk of death from heart disease in patients with nonsmall cell lung cancer who receive postoperative radiotherapy: Analysis of the Surveillance, Epidemiology, and End Results database. Cancer. (2007) 110:911–7. doi: 10.1002/cncr.22845

PubMed Abstract | CrossRef Full Text | Google Scholar

11. Stewart FA. Mechanisms and dose-response relationships for radiation-induced cardiovascular disease. Ann ICRP. (2012) 41:72–9. doi: 10.1016/j.icrp.2012.06.031

PubMed Abstract | CrossRef Full Text | Google Scholar

12. Lipshultz SE, Landy DC, Lopez-Mitnik G, Lipsitz SR, Hinkle AS, Constine LS, et al. Cardiovascular status of childhood cancer survivors exposed and unexposed to cardiotoxic therapy. J Clin Oncol. (2012) 30:1050–7. doi: 10.1200/JCO.2010.33.7907

PubMed Abstract | CrossRef Full Text | Google Scholar

13. Sanoff HK, Deal AM, Krishnamurthy J, Torrice C, Dillon P, Sorrentino J, et al. Effect of cytotoxic chemotherapy on markers of molecular age in patients with breast cancer. J Natl Cancer Inst. (2014) 106:dju057. doi: 10.1093/jnci/dju057

PubMed Abstract | CrossRef Full Text | Google Scholar

14. Hurria A, Jones L, Muss HB. Cancer treatment as an accelerated aging process: Assessment, biomarkers, and interventions. Am Soc Clin Oncol Educ Book. (2016) 35:e516–22. doi: 10.1200/EDBK_156160

CrossRef Full Text | Google Scholar

15. Parrinello S, Samper E, Krtolica A, Goldstein J, Melov S, Campisi J. Oxygen sensitivity severely limits the replicative lifespan of murine fibroblasts. Nat Cell Biol. (2003) 5:741–7. doi: 10.1038/ncb1024

PubMed Abstract | CrossRef Full Text | Google Scholar

16. Bennett G, Papamichos-Chronakis M, Peterson CL. DNA repair choice defines a common pathway for recruitment of chromatin regulators. Nat Commun. (2013) 4:2084. doi: 10.1038/ncomms3084

PubMed Abstract | CrossRef Full Text | Google Scholar

17. Fumagalli M, Rossiello F, Clerici M, Barozzi S, Cittaro D, Kaplunov JM, et al. Telomeric DNA damage is irreparable and causes persistent DNA-damage-response activation. Nat Cell Biol. (2012) 14:355–65. doi: 10.1038/ncb2466

PubMed Abstract | CrossRef Full Text | Google Scholar

18. Salminen A, Kaarniranta K, Kauppinen A. Immunosenescence: The potential role of myeloid-derived suppressor cells (MDSC) in age-related immune deficiency. Cell Mol Life Sci. (2019) 76:1901–18. doi: 10.1007/s00018-019-03048-x

PubMed Abstract | CrossRef Full Text | Google Scholar

19. Lian J, Yue Y, Yu W, Zhang Y. Immunosenescence: A key player in cancer development. J Hematol Oncol. (2020) 13:151. doi: 10.1186/s13045-020-00986-z

PubMed Abstract | CrossRef Full Text | Google Scholar

20. Del Pinto R, Ferri C. Inflammation-accelerated senescence and the cardiovascular system: Mechanisms and perspectives. Int J Mol Sci. (2018) 19:3701. doi: 10.3390/ijms19123701

PubMed Abstract | CrossRef Full Text | Google Scholar

21. Bick AG, Weinstock JS, Nandakumar SK, Fulco CP, Bao EL, Zekavat SM, et al. Inherited causes of clonal haematopoiesis in 97,691 whole genomes. Nature. (2020) 586:763–8.

Google Scholar

22. Jaiswal S, Natarajan P, Silver AJ, Gibson CJ, Bick AG, Shvartz E, et al. Clonal hematopoiesis and risk of atherosclerotic cardiovascular disease. N Engl J Med. (2017) 377:111–21. doi: 10.1056/NEJMoa1701719

PubMed Abstract | CrossRef Full Text | Google Scholar

23. Yu B, Roberts MB, Raffield LM, Zekavat SM, Nguyen NQH, Biggs ML, et al. Supplemental association of clonal hematopoiesis with incident heart failure. J Am Coll Cardiol. (2021) 78:42–52. doi: 10.1016/j.jacc.2021.04.085

PubMed Abstract | CrossRef Full Text | Google Scholar

24. Fuster JJ, MacLauchlan S, Zuriaga MA, Polackal MN, Ostriker AC, Chakraborty R, et al. Clonal hematopoiesis associated with TET2 deficiency accelerates atherosclerosis development in mice. Science. (2017) 355:842–7. doi: 10.1126/science.aag1381

PubMed Abstract | CrossRef Full Text | Google Scholar

25. Coombs CC, Zehir A, Devlin SM, Kishtagari A, Syed A, Jonsson P, et al. Therapy-related clonal hematopoiesis in patients with non-hematologic cancers is common and associated with adverse clinical outcomes. Cell Stem Cell. (2017) 21:374.e–82.e. doi: 10.1016/j.stem.2017.07.010

PubMed Abstract | CrossRef Full Text | Google Scholar

26. Bonnefond A, Skrobek B, Lobbens S, Eury E, Thuillier D, Cauchi S, et al. Association between large detectable clonal mosaicism and type 2 diabetes with vascular complications. Nat Genet. (2013) 45:1040–3. doi: 10.1038/ng.2700

PubMed Abstract | CrossRef Full Text | Google Scholar

27. Dorsheimer L, Assmus B, Rasper T, Ortmann CA, Ecke A, Abou-El-Ardat K, et al. Association of mutations contributing to clonal hematopoiesis with prognosis in chronic ischemic heart failure. JAMA Cardiol. (2019) 4:25–33. doi: 10.1001/jamacardio.2018.3965

PubMed Abstract | CrossRef Full Text | Google Scholar

28. Mas-Peiro S, Hoffmann J, Fichtlscherer S, Dorsheimer L, Rieger MA, Dimmeler S, et al. Clonal haematopoiesis in patients with degenerative aortic valve stenosis undergoing transcatheter aortic valve implantation. Eur Heart J. (2020) 41:933–9. doi: 10.1093/eurheartj/ehz591

PubMed Abstract | CrossRef Full Text | Google Scholar

29. Yoshino J, Baur JA, Imai SI. NAD(+) Intermediates: The biology and therapeutic potential of NMN and NR. Cell Metab. (2018) 27:513–28. doi: 10.1016/j.cmet.2017.11.002

PubMed Abstract | CrossRef Full Text | Google Scholar

30. Kotla S, Zhang A, Imanishi M, Ko KA, Lin SH, Gi YJ, et al. Nucleus-mitochondria positive feedback loop formed by ERK5 S496 phosphorylation-mediated poly (ADP-ribose) polymerase activation provokes persistent pro-inflammatory senescent phenotype and accelerates coronary atherosclerosis after chemo-radiation. Redox Biol. (2021) 47:102132. doi: 10.1016/j.redox.2021.102132

PubMed Abstract | CrossRef Full Text | Google Scholar

31. Hogan KA, Chini CCS, Chini EN. The multi-faceted ecto-enzyme CD38: Roles in immunomodulation, cancer, aging, and metabolic diseases. Front Immunol. (2019) 10:1187. doi: 10.3389/fimmu.2019.01187

PubMed Abstract | CrossRef Full Text | Google Scholar

32. Chini C, Hogan KA, Warner GM, Tarrago MG, Peclat TR, Tchkonia T, et al. The NADase CD38 is induced by factors secreted from senescent cells providing a potential link between senescence and age-related cellular NAD(+) decline. Biochem Biophys Res Commun. (2019) 513:486–93. doi: 10.1016/j.bbrc.2019.03.199

PubMed Abstract | CrossRef Full Text | Google Scholar

33. Lazarevic V, Glimcher LH, Lord GM. T-bet: A bridge between innate and adaptive immunity. Nat Rev Immunol. (2013) 13:777–89. doi: 10.1038/nri3536

PubMed Abstract | CrossRef Full Text | Google Scholar

34. Peng SL, Szabo SJ, Glimcher LH. T-bet regulates IgG class switching and pathogenic autoantibody production. Proc Natl Acad Sci U.S.A. (2002) 99:5545–50. doi: 10.1073/pnas.082114899

PubMed Abstract | CrossRef Full Text | Google Scholar

35. Wang NS, McHeyzer-Williams LJ, Okitsu SL, Burris TP, Reiner SL, McHeyzer-Williams MG. Divergent transcriptional programming of class-specific B cell memory by T-bet and RORalpha. Nat Immunol. (2012) 13:604–11. doi: 10.1038/ni.2294

PubMed Abstract | CrossRef Full Text | Google Scholar

36. Serre K, Cunningham AF, Coughlan RE, Lino AC, Rot A, Hub E, et al. CD8 T cells induce T-bet-dependent migration toward CXCR3 ligands by differentiated B cells produced during responses to alum-protein vaccines. Blood. (2012) 120:4552–9. doi: 10.1182/blood-2012-03-417733

PubMed Abstract | CrossRef Full Text | Google Scholar

37. Cancro MP. Age-associated B cells. Annu Rev Immunol. (2020) 38:315–40. doi: 10.1146/annurev-immunol-092419-031130

PubMed Abstract | CrossRef Full Text | Google Scholar

38. Buono C, Binder CJ, Stavrakis G, Witztum JL, Glimcher LH, Lichtman AH. T-bet deficiency reduces atherosclerosis and alters plaque antigen-specific immune responses. Proc Natl Acad Sci U.S.A. (2005) 102:1596–601. doi: 10.1073/pnas.0409015102

PubMed Abstract | CrossRef Full Text | Google Scholar

39. Banerjee P, Kotla S, Reddy Velatooru L, Abe RJ, Davis EA, Cooke JP, et al. Senescence-associated secretory phenotype as a hinge between cardiovascular diseases and cancer. Front Cardiovasc Med. (2021) 8:763930. doi: 10.3389/fcvm.2021.763930

PubMed Abstract | CrossRef Full Text | Google Scholar

40. Hartmann FJ, Babdor J, Gherardini PF, Amir ED, Jones K, Sahaf B, et al. Comprehensive immune monitoring of clinical trials to advance human immunotherapy. Cell Rep. (2019) 28:819.e–31.e. doi: 10.1016/j.celrep.2019.06.049

PubMed Abstract | CrossRef Full Text | Google Scholar

41. Van Gassen S, Callebaut B, Van Helden MJ, Lambrecht BN, Demeester P, Dhaene T, et al. FlowSOM: Using self-organizing maps for visualization and interpretation of cytometry data. Cytometry A. (2015) 87:636–45. doi: 10.1002/cyto.a.22625

PubMed Abstract | CrossRef Full Text | Google Scholar

42. Amir ED, Lee B, Badoual P, Gordon M, Guo XV, Merad M, et al. Development of a comprehensive antibody staining database using a standardized analytics pipeline. Front Immunol. (2019) 10:1315. doi: 10.3389/fimmu.2019.01315

PubMed Abstract | CrossRef Full Text | Google Scholar

43. Maecker HT, McCoy JP, Nussenblatt R. Standardizing immunophenotyping for the human immunology project. Nat Rev Immunol. (2012) 12:191–200. doi: 10.1038/nri3158

PubMed Abstract | CrossRef Full Text | Google Scholar

44. Finak G, Langweiler M, Jaimes M, Malek M, Taghiyar J, Korin Y, et al. Standardizing flow cytometry immunophenotyping analysis from the human immunophenotyping consortium. Sci Rep. (2016) 6:20686. doi: 10.1038/srep20686

PubMed Abstract | CrossRef Full Text | Google Scholar

45. Crosby EJ, Hobeika AC, Niedzwiecki D, Rushing C, Hsu D, Berglund P, et al. Long-term survival of patients with stage III colon cancer treated with VRP-CEA(6D), an alphavirus vector that increases the CD8+ effector memory T cell to Treg ratio. J Immunother Cancer. (2020) 8:e001662. doi: 10.1136/jitc-2020-001662

PubMed Abstract | CrossRef Full Text | Google Scholar

46. Robinson MD, McCarthy DJ, Smyth GK. edgeR: A Bioconductor package for differential expression analysis of digital gene expression data. Bioinformatics. (2010) 26:139–40. doi: 10.1093/bioinformatics/btp616

PubMed Abstract | CrossRef Full Text | Google Scholar

47. McCarthy DJ, Chen Y, Smyth GK. Differential expression analysis of multifactor RNA-Seq experiments with respect to biological variation. Nucleic Acids Res. (2012) 40:4288–97. doi: 10.1093/nar/gks042

PubMed Abstract | CrossRef Full Text | Google Scholar

48. Lun ATL, Richard AC, Marioni JC. Testing for differential abundance in mass cytometry data. Nat Methods. (2017) 14:707–9. doi: 10.1038/nmeth.4295

PubMed Abstract | CrossRef Full Text | Google Scholar

49. Ritchie ME, Phipson B, Wu D, Hu Y, Law CW, Shi W, et al. limma powers differential expression analyses for RNA-sequencing and microarray studies. Nucleic Acids Res. (2015) 43:e47. doi: 10.1093/nar/gkv007

PubMed Abstract | CrossRef Full Text | Google Scholar

50. Weber LM, Nowicka M, Soneson C, Robinson MD. diffcyt: Differential discovery in high-dimensional cytometry via high-resolution clustering. Commun Biol. (2019) 2:183. doi: 10.1038/s42003-019-0415-5

PubMed Abstract | CrossRef Full Text | Google Scholar

51. Weiskopf D, Bangs DJ, Sidney J, Kolla RV, De Silva AD, de Silva AM, et al. Dengue virus infection elicits highly polarized CX3CR1+ cytotoxic CD4+ T cells associated with protective immunity. Proc Natl Acad Sci U.S.A. (2015) 112:E4256–63. doi: 10.1073/pnas.1505956112

PubMed Abstract | CrossRef Full Text | Google Scholar

52. Nishimura M, Umehara H, Nakayama T, Yoneda O, Hieshima K, Kakizaki M, et al. Dual functions of fractalkine/CX3C ligand 1 in trafficking of perforin+/granzyme B+ cytotoxic effector lymphocytes that are defined by CX3CR1 expression. J Immunol. (2002) 168:6173–80. doi: 10.4049/jimmunol.168.12.6173

PubMed Abstract | CrossRef Full Text | Google Scholar

53. Bottcher JP, Beyer M, Meissner F, Abdullah Z, Sander J, Hochst B, et al. Functional classification of memory CD8(+) T cells by CX3CR1 expression. Nat Commun. (2015) 6:8306. doi: 10.1038/ncomms9306

PubMed Abstract | CrossRef Full Text | Google Scholar

54. Pachnio A, Ciaurriz M, Begum J, Lal N, Zuo J, Beggs A, et al. Cytomegalovirus infection leads to development of high frequencies of cytotoxic virus-specific CD4+ T cells targeted to vascular endothelium. PLoS Pathog. (2016) 12:e1005832. doi: 10.1371/journal.ppat.1005832

PubMed Abstract | CrossRef Full Text | Google Scholar

55. Karnell JL, Kumar V, Wang J, Wang S, Voynova E, Ettinger R. Role of CD11c(+) T-bet(+) B cells in human health and disease. Cell Immunol. (2017) 321:40–5. doi: 10.1016/j.cellimm.2017.05.008

PubMed Abstract | CrossRef Full Text | Google Scholar

56. Singh MV, Kotla S, Le NT, Ae Ko K, Heo KS, Wang Y, et al. Senescent phenotype induced by p90RSK-NRF2 signaling sensitizes monocytes and macrophages to oxidative stress in HIV-positive individuals. Circulation. (2019) 139:1199–216. doi: 10.1161/CIRCULATIONAHA.118.036232

PubMed Abstract | CrossRef Full Text | Google Scholar

57. Le NT, Heo KS, Takei Y, Lee H, Woo CH, Chang E, et al. A crucial role for p90RSK-mediated reduction of ERK5 transcriptional activity in endothelial dysfunction and atherosclerosis. Circulation. (2013) 127:486–99. doi: 10.1161/CIRCULATIONAHA.112.116988

PubMed Abstract | CrossRef Full Text | Google Scholar

58. Winslow GM, Papillion AM, Kenderes KJ, Levack RC. CD11c+ T-bet+ memory B cells: Immune maintenance during chronic infection and inflammation? Cell Immunol. (2017) 321:8–17. doi: 10.1016/j.cellimm.2017.07.006

PubMed Abstract | CrossRef Full Text | Google Scholar

59. Klein U, Rajewsky K, Kuppers R. Human immunoglobulin (Ig)M+IgD+ peripheral blood B cells expressing the CD27 cell surface antigen carry somatically mutated variable region genes: CD27 as a general marker for somatically mutated (memory) B cells. J Exp Med. (1998) 188:1679–89. doi: 10.1084/jem.188.9.1679

PubMed Abstract | CrossRef Full Text | Google Scholar

60. Agematsu K, Nagumo H, Yang FC, Nakazawa T, Fukushima K, Ito S, et al. B cell subpopulations separated by CD27 and crucial collaboration of CD27+ B cells and helper T cells in immunoglobulin production. Eur J Immunol. (1997) 27:2073–9. doi: 10.1002/eji.1830270835

PubMed Abstract | CrossRef Full Text | Google Scholar

61. Kobata T, Jacquot S, Kozlowski S, Agematsu K, Schlossman SF, Morimoto C. CD27-CD70 interactions regulate B-cell activation by T cells. Proc Natl Acad Sci U.S.A. (1995) 92:11249–53. doi: 10.1073/pnas.92.24.11249

PubMed Abstract | CrossRef Full Text | Google Scholar

62. Lens SM, Keehnen RM, van Oers MH, van Lier RA, Pals ST, Koopman G. Identification of a novel subpopulation of germinal center B cells characterized by expression of IgD and CD70. Eur J Immunol. (1996) 26:1007–11. doi: 10.1002/eji.1830260508

PubMed Abstract | CrossRef Full Text | Google Scholar

63. Arens R, Nolte MA, Tesselaar K, Heemskerk B, Reedquist KA, van Lier RA, et al. Signaling through CD70 regulates B cell activation and IgG production. J Immunol. (2004) 173:3901–8. doi: 10.4049/jimmunol.173.6.3901

PubMed Abstract | CrossRef Full Text | Google Scholar

64. Mahajan VS, Mattoo H, Sun N, Viswanadham V, Yuen GJ, Allard-Chamard H, et al. B1a and B2 cells are characterized by distinct CpG modification states at DNMT3A-maintained enhancers. Nat Commun. (2021) 12:2208. doi: 10.1038/s41467-021-22458-9

PubMed Abstract | CrossRef Full Text | Google Scholar

65. Haney SL, Upchurch GM, Opavska J, Klinkebiel D, Hlady RA, Suresh A, et al. Promoter hypomethylation and expression is conserved in mouse chronic lymphocytic leukemia induced by decreased or inactivated Dnmt3a. Cell Rep. (2016) 15:1190–201. doi: 10.1016/j.celrep.2016.04.004

PubMed Abstract | CrossRef Full Text | Google Scholar

66. Pittner BT, Shanafelt TD, Kay NE, Jelinek DF. CD38 expression levels in chronic lymphocytic leukemia B cells are associated with activation marker expression and differential responses to interferon stimulation. Leukemia. (2005) 19:2264–72. doi: 10.1038/sj.leu.2403975

PubMed Abstract | CrossRef Full Text | Google Scholar

67. Mouat IC, Goldberg E, Horwitz MS. Age-associated B cells in autoimmune diseases. Cell Mol Life Sci. (2022) 79:402. doi: 10.1007/s00018-022-04433-9

PubMed Abstract | CrossRef Full Text | Google Scholar

68. Rubtsov AV, Rubtsova K, Fischer A, Meehan RT, Gillis JZ, Kappler JW, et al. Toll-like receptor 7 (TLR7)-driven accumulation of a novel CD11c(+) B-cell population is important for the development of autoimmunity. Blood. (2011) 118:1305–15. doi: 10.1182/blood-2011-01-331462

PubMed Abstract | CrossRef Full Text | Google Scholar

69. Rubtsova K, Rubtsov AV, van Dyk LF, Kappler JW, Marrack P. T-box transcription factor T-bet, a key player in a unique type of B-cell activation essential for effective viral clearance. Proc Natl Acad Sci U.S.A. (2013) 110:E3216–24. doi: 10.1073/pnas.1312348110

PubMed Abstract | CrossRef Full Text | Google Scholar

70. Rubtsova K, Rubtsov AV, Cancro MP, Marrack P. Age-associated B cells: A T-bet-dependent effector with roles in protective and pathogenic immunity. J Immunol. (2015) 195:1933–7. doi: 10.4049/jimmunol.1501209

PubMed Abstract | CrossRef Full Text | Google Scholar

71. Piedra-Quintero ZL, Wilson Z, Nava P, Guerau-de-Arellano M. CD38: An immunomodulatory molecule in inflammation and autoimmunity. Front Immunol. (2020) 11:597959. doi: 10.3389/fimmu.2020.597959

PubMed Abstract | CrossRef Full Text | Google Scholar

72. Burns M, Ostendorf L, Biesen R, Grutzkau A, Hiepe F, Mei HE, et al. Dysregulated CD38 expression on peripheral blood immune cell subsets in SLE. Int J Mol Sci. (2021) 22:2424. doi: 10.3390/ijms22052424

PubMed Abstract | CrossRef Full Text | Google Scholar

73. Burlock B, Richardson G, Garcia-Rodriguez S, Guerrero S, Zubiaur M, Sancho J. The role of CD38 on the function of regulatory B cells in a murine model of lupus. Int J Mol Sci. (2018) 19:2906. doi: 10.3390/ijms19102906

PubMed Abstract | CrossRef Full Text | Google Scholar

74. Burgler S, Gimeno A, Parente-Ribes A, Wang D, Os A, Devereux S, et al. Chronic lymphocytic leukemia cells express CD38 in response to Th1 cell-derived IFN-gamma by a T-bet-dependent mechanism. J Immunol. (2015) 194:827–35. doi: 10.4049/jimmunol.1401350

PubMed Abstract | CrossRef Full Text | Google Scholar