Abstract
Sipuleucel-T is an autologous cellular immunotherapy, administered as three infusions, for metastatic castration-resistant prostate cancer (mCRPC). Sipuleucel-T induces T- and B-cell responses to prostatic acid phosphatase (PAP), correlating to improved survival. The long-term impact of sipuleucel-T on tumor antigen–specific immunologic memory remains unknown, in particular, B-cell responses, as measured by antigen-specific antibody responses and B-cell receptor (BCR) sequences. To evaluate whether sipuleucel-T could induce long-term immunologic memory, we examined circulating B-cell responses before and after sipuleucel-T treatment in two groups of patients with mCRPC: those who had previously received sipuleucel-T (treated; median, 8.9 years since the previous treatment) versus those who had not (naïve). Before re-treatment, previously treated patients exhibited persistent antibody responses as well as more focused and convergent BCR repertoires with distinct V(D)J gene usage compared with naïve patients. After re-treatment, previously treated patients maintained high-frequency clones and developed more convergent BCRs at earlier time points unlike naïve patients. With the first sipuleucel-T infusion specifically, previously treated patients had less shuffling within the 100 most abundant baseline clones. In contrast, naïve patients exhibited great BCR turnover with a continued influx of new B-cell clones. Social network analysis showed that previously treated patients had more highly organized B-cell repertoires, consistent with greater clonal maturation. Higher treatment-induced BCR clonality correlated with longer survival for naïve patients. These results demonstrated the capacity of sipuleucel-T to induce long-term immune memory and lasting changes to the B-cell repertoire.
Introduction
Sipuleucel-T is an autologous cellular immunotherapy FDA-approved for the treatment of asymptomatic or minimally symptomatic metastatic castration-resistant prostate cancer (mCRPC; ref. 1). Sipuleucel-T is manufactured by culturing isolated peripheral blood mononuclear cells (PBMC) with PA2024, a recombinant fusion protein composed of prostatic acid phosphatase (PAP) linked to granulocyte-macrophage colony–stimulating factor (1). Sipuleucel-T functions as a cancer therapy by priming T- and B-cell immune responses to PA2024, both of which associate improved overall survival (2–4). This treatment functions like a vaccine, inducing increases in both T-cell repertoire diversity and the number of activated T cells within the tumor microenvironment (5, 6). Most studies of cancer immunotherapies thus far have focused on their impact on the T-cell repertoire, largely ignoring the B-cell response. The durability of antigen-specific responses induced with cancer immunotherapy is unknown. Here, we assessed the capacity of sipuleucel-T to induce acute and memory B-cell responses by measuring antigen-specific antibodies and determining the treatment induced effects in the B-cell repertoire. We used social network analysis to assess the maturity of the B-cell memory response, because affinity maturation can result in highly related clonotypes.
The B-cell receptor (BCR) is an antibody with a heterodimeric transmembrane protein complex that comprises two identical heavy-chain (IGH) and two identical light-chain proteins (7, 8). BCR genes include variable (V), diversity (D), and joining (J) gene segments that contribute to receptor diversity through their rearrangements (9, 10). Complementarity-determining region 3 (CDR3) is a part of the rearranged BCR critical for antigen recognition. Antigen-specific BCRs can develop improved affinity through somatic hypermutation and clonal affinity selection resulting in increased diversity, and the sequences of these clones will be highly related (11, 12). Analysis of these antibodies could reveal clonal families of B cells that persist over time and show evidence of progressive affinity maturation and class switching (13).
Immunologic studies conducted on samples from a phase III trial that demonstrates improved survival with sipuleucel-T treatment in mCRPC shows that 68% of treated patients developed antibody responses against PA2024, the immunogen used to make the product (2, 3). To investigate whether prior treatment with sipuleucel-T induced long-term immunologic memory, we compared treatment-induced B-cell responses in patients who were previously treated with sipuleucel-T (treated; P10-1 study, NCT01338012; Supplementary Fig. S1) and patients who received their first sipuleucel-T treatment (naïve; STRIDE study, NCT01981122; Supplementary Fig. S1). Among the 8 patients in P10-1 study, the median overall survival (OS) defined as from the registration date to death of any cause or last follow up was 31.3 months. For the 52 patients in STRIDE study, the median OS was 32.7 months, which was defined as from the date of enrollment to death of any cause or last follow up. However, the two original studies, STRIDE and P10-1, were not designed to assess BCR sequencing. This current work encompasses experiments using available specimens with the goal to conduct this exploratory study focusing on B-cell memory induced with sipuleucel-T treatment in patients with mCRPC.
Materials and Methods
Clinical studies
Naïve patients in the current study previously participated in a multicenter, randomized, open-label phase II clinical trial (STRIDE trial, NCT01981122) in which patients with asymptomatic or minimally symptomatic mCRPC received their first course of sipuleucel-T with either concurrent enzalutamide (starting 2 weeks before sipuleucel-T initiation) or sequential enzalutamide (10 weeks after sipuleucel-T initiation; Supplementary Fig. S1; ref. 14). As there was no difference in the metrics of either BCR diversity or dynamics across time points between concurrent and sequential arms within naïve patients, we combined the two arms.
Treated patients in this study participated in a single-arm, open-label, phase II, clinical trial that enrolled patients previously treated with sipuleucel-T (P10-1 trial, NCT01338012). These patients previously had received sipuleucel-T in the androgen-sensitive setting (PROTECT trial, NCT00779402). The P10-1 study that explored re-treatment with sipuleucel-T only included those patients from PROTECT subsequently progressed to mCRPC (Supplementary Fig. S1) and none of them were complete responders. Patients were retreated with a second full course of sipuleucel-T after a median of 8.9 years as they developed mCRPC (15). This is actually the disease setting where sipuleucel-T confers a survival benefit in patients with mCRPC in phase III clinical trials and is FDA approved. Induction of antibody responses is associated with improved outcomes. This study was conducted to examine the durability of immune response to sipuleucel-T, as well as to determine whether additional treatment with sipuleucel-T further alters the immune response.
Blood samples evaluated in this study came from the above two previous clinical trials of sipuleucel-T, with 7 patients from P10-1 trial (NCT01338012) and 19 patients from STRIDE trial (NCT01981122). Each patient in these trials provided informed consent; and relevant institutional review boards approved the study protocols, including the collection of blood, for these previous trials.
Both studies were conducted in accordance with the ethical guideline of Declaration of Helsinki and were performed after approval by an institutional review board (IRB). All the investigators were obtained informed written consent from the subjects.
Samples and BCR sequencing
In each study, peripheral blood mononuclear cells (PBMC) were collected and cryopreserved at baseline (week 0), during (weeks 2 and 4), and after sipuleucel-T therapy (weeks 6, 26, and 52). For each PBMC sample, immunoglobulin heavy-chain CDR3 regions were amplified and sequenced using 2 μg of genomic DNA (or all available extracted DNA if the amount available was less than 2 μg) by Adaptive Biotechnologies using the immunoSEQ platform (16). BCR-sequencing datasets generated during the study are available via https://clients.adaptivebiotech.com/pub/zhang-2020-cir.
Antibody response assessment
IgG levels against PA2024 and prostate-specific antigen (PSA) were evaluated by Life Technologies Corporation using Luminex xMAP technology, which uses multiplexed antigen-coated spectrally distinguishable fluorescence dyed beads (17). The beads were conjugated to PA2024, PAP, PSA, ERAS, KRAS, LGALS3, LGALS8, and KLK2 as part of a custom kit. Serum samples (30 μL) were assessed at 1:200 dilution, and normalized signal intensities were log2-transformed for analysis. Nonspecific background signal was measured using bovine serum albumin, Glutathione S-transferase (GST) and anti-GST–conjugated beads, and data was normalization based on linear model to reduce technical variability (18). Changes in serum IgG levels were reported relative to baseline (week 0).
BCR data assessment
BCR data analysis was performed using the TCR3D software developed by our group (19). Only clones with a count of at least 5 were retained in the analysis. Diversity of the BCR repertoire at each time point was measured by clonality on a scale of 0 to 1, where 0 indicates total diversity (i.e., all BCR clonotypes equally common) and 1 indicates no diversity (i.e., the BCR repertoire is dominated by a single clone). Clonality can be considered as a normalized Shannon entropy index (H) over the number of unique clones (n): Clonality = 1 − H/loge(n), where H/loge(n) is Pielou evenness (equability; ref. 19). Convergent frequency was calculated as the aggregate frequency of clones sharing an amino acid sequence with at least one other clone (20).
Morisita distance, a distance measurement ranging from 0 (maximally dissimilar) to 1 (minimally dissimilar), was calculated to examine the dynamic change in BCR repertoire from one time point to another time point for each patient (19). In addition, clonotypes were clustered into three groups: preexisting clones that were only present in the earlier time point, newly generated clones that were only present in the later time point, and overlap clones that were present at both time points. The proportion of each of the three groups was calculated by the number of the clones belong to the corresponding group divided by the total number of clones presented at either time point.
In addition, the top 100 clonotypes were identified at every time point based on their abundances and then the ranks of those clones at other time points were also obtained. Intraclass correlation coefficients (ICC) of the ranks of those top 100 clones at each time point were calculated. The higher the ICC, the more consistent the rank order was across time points, that is, less changes across the time points and therefore lower clonal shuffling.
VJ gene usage was defined as the number of clonotypes that utilize the same combination of V and J genes normalized by the total number of unique clones. The random forest method was used to identify the genes whose usages at each time point are differential (21). Genes were clustered by unsupervised hierarchical clustering. Results were illustrated using heatmap.
Visualization of somatic hypermutation and affinity maturation
Similarly, as in ref. 22, for each patient, pairwise distance matrix of each pair of full-length of BCR clonotype sequences was calculated on Levenshtein distance (R Package: RecordLinkage; ref. 23) within each V gene family. V gene families with greater than 2,000 clones were separated into V, D, and J gene levels. Because Levenshtein Distance detects the similarities that include insertion, deletion, and transformation, we included sequences with different length. A convergent group was defined as the cluster that included the clones with the distance less than or equal to 1 (allowing maximum of one base pair mutation among clone sequences sharing the same V-gene, J-gene, and CDR3 length).
Social network visualization was performed for each V gene family using R packages: Ape (24) and igraph (25). Each node represents a single full-length BCR clonotype sequence colored by its first presenting time point or by corresponding CDR3 amino acid sequence, and node size represents the corresponding abundance of that clone. The phylogenetic tree was plotted by R package: ggtree (26) with edge color coded by CDR3 amino acid sequences. For a selected representative convergent group, a chord diagram (R package: circlize; ref. 27) was plotted to show the clonal development across time points. The outer circle color stands for the time point and the width of the arch for the abundance information (log10 transformed). Motif analysis was performed by R package: ggseglogo (28). The group size (the number of clones belong to each group), and the maximum group size (the maximized number of clones belong to each group across all groups within the same sample) within each network were used to describe the groups for each sample.
As clonal grouping can be impacted by experimental factors such as sampling and sequencing depth, we performed a sensitivity analysis. We subsampled clones from representative patients to achieve similar sequencing depth, where subsampling was weighted by the distribution of normalized abundance, and then performed the above visualization and analysis pipeline.
Statistical analysis
In general, frequency distribution and percentages were used to summarize categorical variables, and median with range (min–max) was used to summarize continuous variables. Comparison of two studies was performed using Wilcoxon rank sum test for continuous variables and Pearson χ2 test for categorical variables. Diversity or dynamic indices between time points were compared using the Wilcoxon signed rank test for within study comparisons and comparison between two studies was performed using Wilcoxon rank sum test. For the STRIDE trial, OS, which was defined as the duration between the date of trial enrollment to the date of death due to any cause, was estimated by the Kaplan–Meier method. The relationship between OS and BCR diversity or dynamics was analyzed using a log-rank test and Cox proportional hazards models. Statistical significance was declared at P < 0.05, and no multiple testing adjustment was done. All statistical analysis was done with the statistical computing software R (https://www.r-project.org/).
Results
Prior sipuleucel-T treatment resulted in long-term B-cell responses
Table 1 lists the key demographic and clinical characteristics along with BCR-sequencing characteristics of both naïve and treated patients. There were 7 treated patients and 19 naïve patients with BCR sequencing, who had similar demographic and clinical characteristics. All were Caucasians with median age of 71 years old.
Key demographic and clinical characteristics of patients from sipuleucel-T treated and naïve trials with BCR sequencing data.
. | Treated . | Naïve . |
---|---|---|
Characteristics . | (N = 7) . | (N = 19) . |
Age (years), median (range)a | 71 (67–84) | 71 (51–84) |
Caucasian race, n (%) | 7 (100) | 19 (100) |
ECOG PS 0, n (%) | 6 (86) | 16 (84) |
Gleason score, n (%) | ||
≤7 | 3 (43) | 6 (32) |
≥8 | 4 (57) | 13 (68) |
PSA (ng/mL) median (range)a | 20 (5–167) | 9 (2–59) |
LDH (U/L), median (range)a | 163 (138–205) | 189 (86–237) |
Alkaline phosphatase (U/L), median (range)a | 77 (65–83) | 87 (50–271) |
Prior chemotherapy, n (%) | 1 (14) | 0 (0) |
. | Treated . | Naïve . |
---|---|---|
Characteristics . | (N = 7) . | (N = 19) . |
Age (years), median (range)a | 71 (67–84) | 71 (51–84) |
Caucasian race, n (%) | 7 (100) | 19 (100) |
ECOG PS 0, n (%) | 6 (86) | 16 (84) |
Gleason score, n (%) | ||
≤7 | 3 (43) | 6 (32) |
≥8 | 4 (57) | 13 (68) |
PSA (ng/mL) median (range)a | 20 (5–167) | 9 (2–59) |
LDH (U/L), median (range)a | 163 (138–205) | 189 (86–237) |
Alkaline phosphatase (U/L), median (range)a | 77 (65–83) | 87 (50–271) |
Prior chemotherapy, n (%) | 1 (14) | 0 (0) |
Abbreviations: ECOG PS, Eastern Cooperative Oncology Group performance status; LDH, lactate dehydrogenase; PSA, prostate-specific antigen.
aRange = min–max.
Baseline sera from each study was examined; week 0 was defined as before initial exposure among the naïve patients, and before reexposure in the treated patients. First, we compared the baseline sera in the two groups to determine whether sipuleucel-T induced a persistent antibody response. The levels of antibodies to PA2024, the immunogen in sipuleucel-T, were significantly higher in treated patients than in naïve patients (PA2024, Wilcoxon rank sum test P < 0.001; Fig. 1A), whereas the levels of antibodies to PSA were not. Also, baseline BCR clonality was significantly higher in treated patients than in naïve patients (Wilcoxon rank sum test P = 0.003; Fig. 1B), indicating that the treated patients had a more focused repertoire before being retreated with sipuleucel-T. The baseline BCR convergent frequency was also significantly higher in treated patients than in naïve patients (Wilcoxon rank sum test P = 0.003; Fig. 1C), suggesting persistent response to previous exposure. In addition, treated patients had distinct VDJ gene usage relative to naïve patients (Fig. 1D), where the genes were selected by random forest based on the differential usages between treated and naïve patient repertoires at baseline. Genes were clustered by unsupervised hierarchical clustering.
Baseline differences in antigen and B-cell repertoire between 7 treated and 19 naïve patients. A, Log2 of fluorescence intensity of antigen PA2024 (left) and PSA (right). In each panel, left and right boxplots represent log2 of fluorescence intensity of treated and naïve patients, respectively. The asterisk shows significant difference between treated patients and naïve patients (Wilcoxon rank sum test P < 0.05). B, Clonality at baseline for treated (left boxplot) and naïve (right boxplot) patients. For each repertoire, clonality was calculated as 1 − H/loge(n), where H is Shannon entropy index and n is the number of unique clones (19). Clonality of 0 and 1 indicate total diversity and no diversity. Clonality in treated patients was significantly higher than in naïve patients (Wilcoxon rank sum test P = 0.003). The asterisk shows significant difference between treated patients and naïve patients (Wilcoxon rank sum test P < 0.05). C, BCR convergent frequency at baseline for both treated (left boxplot) and naïve (right boxplot) patients.Convergent frequency was calculated as the aggregate frequency of clones sharing an amino acid sequence with at least one other clone (20). BCR convergent frequency in treated patients was significantly higher than in naïve patients (Wilcoxon rank sum test P = 0.003). The asterisk shows significant difference between treated patients and naïve patients (Wilcoxon rank sum test P < 0.05). D, Heatmap of VJ gene usages at baseline. The genes were selected by random forest based on the differential gene usages between treated and naïve patient repertories at baseline. Genes were clustered by unsupervised hierarchical clustering. Heatmap column represents single VJ gene combination, and row represents individual patient, with blue and orange bars on left side representing treated and naïve patients, respectively. Scaled color code stands for the normalized gene usages, with darker red corresponding to higher gene usage. Social network figure at baseline for the representative treated patient (E) and for the representative naïve patient (F). For each patient, pairwise distance matrix of each pair of full-length of BCR clonotype sequences was calculated on the basis of Levenshtein distance (R package: RecordLinkage; ref. 23) within each V gene family. Social network analysis was then performed using R packages: Ape (24) and igraph (25). Each node represents a full-length of single clone colored by amino acid sequence, and the node size is attributed on the basis of the corresponding abundance. The nodes connected by lines were within the same convergent group that included the full-length of BCR clonotype sequences with the distance less than or equal to 1 (allowing maximum of one base pair mutation among clone sequences sharing the same V-gene, J-gene, and CDR3 length). One group was picked from each figure as an example of actual clone sequences (first 70 base pairs due to space limit), with the yellow highlights representing differences in the nucleotide sequence. G, Violin plots of the number of BCR convergent groups at baseline. The number of convergent groups within each network was used to describe the feature of the network. Left and right violin plots illustrate the number of BCR convergent groups at baseline for treated patients and naïve patients, respectively. The asterisk shows significant difference between treated patients and naïve patients (Wilcoxon rank sum test P < 0.05). H, Violin plots of the size of BCR convergent groups at baseline. The maximum group size within each network was used to describe the feature of the network. Left and right violin plots illustrate the maximum number of clones within each BCR convergent group at baseline for treated patients and naïve patients, respectively. The asterisk shows significant difference between treated patients and naïve patients (Wilcoxon rank sum test P < 0.05). All boxplots present 25% percentile, median, and 75% percentile of the data, with each dot representing each sample. All violin plots show distributions of the data smoothed by a kernel density estimator.
Baseline differences in antigen and B-cell repertoire between 7 treated and 19 naïve patients. A, Log2 of fluorescence intensity of antigen PA2024 (left) and PSA (right). In each panel, left and right boxplots represent log2 of fluorescence intensity of treated and naïve patients, respectively. The asterisk shows significant difference between treated patients and naïve patients (Wilcoxon rank sum test P < 0.05). B, Clonality at baseline for treated (left boxplot) and naïve (right boxplot) patients. For each repertoire, clonality was calculated as 1 − H/loge(n), where H is Shannon entropy index and n is the number of unique clones (19). Clonality of 0 and 1 indicate total diversity and no diversity. Clonality in treated patients was significantly higher than in naïve patients (Wilcoxon rank sum test P = 0.003). The asterisk shows significant difference between treated patients and naïve patients (Wilcoxon rank sum test P < 0.05). C, BCR convergent frequency at baseline for both treated (left boxplot) and naïve (right boxplot) patients.Convergent frequency was calculated as the aggregate frequency of clones sharing an amino acid sequence with at least one other clone (20). BCR convergent frequency in treated patients was significantly higher than in naïve patients (Wilcoxon rank sum test P = 0.003). The asterisk shows significant difference between treated patients and naïve patients (Wilcoxon rank sum test P < 0.05). D, Heatmap of VJ gene usages at baseline. The genes were selected by random forest based on the differential gene usages between treated and naïve patient repertories at baseline. Genes were clustered by unsupervised hierarchical clustering. Heatmap column represents single VJ gene combination, and row represents individual patient, with blue and orange bars on left side representing treated and naïve patients, respectively. Scaled color code stands for the normalized gene usages, with darker red corresponding to higher gene usage. Social network figure at baseline for the representative treated patient (E) and for the representative naïve patient (F). For each patient, pairwise distance matrix of each pair of full-length of BCR clonotype sequences was calculated on the basis of Levenshtein distance (R package: RecordLinkage; ref. 23) within each V gene family. Social network analysis was then performed using R packages: Ape (24) and igraph (25). Each node represents a full-length of single clone colored by amino acid sequence, and the node size is attributed on the basis of the corresponding abundance. The nodes connected by lines were within the same convergent group that included the full-length of BCR clonotype sequences with the distance less than or equal to 1 (allowing maximum of one base pair mutation among clone sequences sharing the same V-gene, J-gene, and CDR3 length). One group was picked from each figure as an example of actual clone sequences (first 70 base pairs due to space limit), with the yellow highlights representing differences in the nucleotide sequence. G, Violin plots of the number of BCR convergent groups at baseline. The number of convergent groups within each network was used to describe the feature of the network. Left and right violin plots illustrate the number of BCR convergent groups at baseline for treated patients and naïve patients, respectively. The asterisk shows significant difference between treated patients and naïve patients (Wilcoxon rank sum test P < 0.05). H, Violin plots of the size of BCR convergent groups at baseline. The maximum group size within each network was used to describe the feature of the network. Left and right violin plots illustrate the maximum number of clones within each BCR convergent group at baseline for treated patients and naïve patients, respectively. The asterisk shows significant difference between treated patients and naïve patients (Wilcoxon rank sum test P < 0.05). All boxplots present 25% percentile, median, and 75% percentile of the data, with each dot representing each sample. All violin plots show distributions of the data smoothed by a kernel density estimator.
Social network analysis (see the subsection of Visualization of Somatic Hypermutation and Affinity Maturation in Method Section) showed that the representative treated patient had 18 convergent groups at baseline, with each group having more than four clones (Fig. 1E), whereas the representative naïve patient had only one convergent group with more than four clones (Fig. 1F). A convergent group included all clones with maximum of one base pair mutation who share the same V-gene, J-gene, and CDR3 length. Overall, at week 0, the treated patients had more convergent groups (Fig. 1G; Wilcoxon rank sum test P = 0.003) with larger group sizes (Fig. 1H; Wilcoxon rank sum test P = 0.003). These differences were consistent with long-term immunologic memory from prior treatment.
Dynamic changes with sipuleucel-T treatment differed between naïve and treated patients
Next, we assessed immune responses induced by new sipuleucel-T treatment in both the naïve and treated patients, hypothesizing that the magnitude and kinetics of a response would differ with priming than with recall responses. Induction of antibody responses against PA2024 was significantly greater between week 0 and later time points (except week 26) for naïve patients than for treated patients (Fig. 2A; Wilcoxon rank sum test P = 0.005–0.046), a pattern that was not observed for PSA, another prostate antigen that is not contained in sipuleucel-T (Fig. 2B). Yet, increases in BCR clonality occurred faster and were larger in treated patients than in naïve patients. In addition, BCR clonality was significantly higher (Wilcoxon rank sum test P ≤ 0.001) in treated patients than in naïve patients at every post-baseline time point except week 6, suggesting a more focused antibody response (i.e., with less clonal diversity; Fig. 2C).
Changes in antigen BCR repertoire with sipuleucel-T treatment between 7 treated and 19 naïve patients. A and B, Log2 of fold difference in fluorescence intensity relative to week 0 of antigen PA2024 (A) and PSA (B) at all post-baseline time points. In each figure, left and right boxplots represent log2 of fold difference in fluorescent intensity relative to week 0 of treated and naïve patients across all post-baseline time points. The asterisk shows significant difference between treated patients and naïve patients (Wilcoxon rank sum test P < 0.05). C, Clonality at all time points for treated (left boxplot) and naïve (right boxplot) patients. For each repertoire, clonality was calculated as 1 − H/loge(n), where H is Shannon entropy index and n is the number of unique clones (19). Clonality of 0 and 1 indicate total diversity and no diversity. The asterisk shows significant increase of clonality in treated patients compared with naïve patients (Wilcoxon rank sum test P < 0.05). Changes in ranks of the top 100 BCR clones at week 0 for a representative treated patient (D) and naïve patient (E) are shown. For each patient, the top 100 clonotypes at week 0 were identified on the basis of their abundances at week 0. The ranks of those clones at other time points were also found and plotted across all six time points. Each dot connected line represents a single clone. F, ICCs of the ranks of top 100 clones at all time points. For each patient, the top 100 clonotypes at week 0 were identified based on their abundances at week 0, and then the ranks of those clones at other time points were also obtained. Intraclass correlation of the ranks of those top 100 clones at each time point was calculated (19). The higher the ICC, the more consistent the rank order was across time points, that is, less changes across the time points and therefore lower clonal shuffling. The asterisk shows that at the corresponding time point, ICC was significantly higher (Wilcoxon rank sum test P < 0.05) for treated patients compared to naïve patients. All boxplots present 25% percentile, median, and 75% percentile of the data, with each dot representing each sample.
Changes in antigen BCR repertoire with sipuleucel-T treatment between 7 treated and 19 naïve patients. A and B, Log2 of fold difference in fluorescence intensity relative to week 0 of antigen PA2024 (A) and PSA (B) at all post-baseline time points. In each figure, left and right boxplots represent log2 of fold difference in fluorescent intensity relative to week 0 of treated and naïve patients across all post-baseline time points. The asterisk shows significant difference between treated patients and naïve patients (Wilcoxon rank sum test P < 0.05). C, Clonality at all time points for treated (left boxplot) and naïve (right boxplot) patients. For each repertoire, clonality was calculated as 1 − H/loge(n), where H is Shannon entropy index and n is the number of unique clones (19). Clonality of 0 and 1 indicate total diversity and no diversity. The asterisk shows significant increase of clonality in treated patients compared with naïve patients (Wilcoxon rank sum test P < 0.05). Changes in ranks of the top 100 BCR clones at week 0 for a representative treated patient (D) and naïve patient (E) are shown. For each patient, the top 100 clonotypes at week 0 were identified on the basis of their abundances at week 0. The ranks of those clones at other time points were also found and plotted across all six time points. Each dot connected line represents a single clone. F, ICCs of the ranks of top 100 clones at all time points. For each patient, the top 100 clonotypes at week 0 were identified based on their abundances at week 0, and then the ranks of those clones at other time points were also obtained. Intraclass correlation of the ranks of those top 100 clones at each time point was calculated (19). The higher the ICC, the more consistent the rank order was across time points, that is, less changes across the time points and therefore lower clonal shuffling. The asterisk shows that at the corresponding time point, ICC was significantly higher (Wilcoxon rank sum test P < 0.05) for treated patients compared to naïve patients. All boxplots present 25% percentile, median, and 75% percentile of the data, with each dot representing each sample.
As shown in Table 2, treated patients generated more de novo clones at week 2 (Wilcoxon rank sum test P = 0.021) and week 4 (Wilcoxon rank sum test P = 0.002) relative to week 0, that is, new clones that were not evident at week 0. However, starting at week 2, more clones persisted in the repertoire in treated patients, that is, the overlapping clones, indicating that sipuleucel‐T had already induced high frequency clones even after a single infusion of sipuleucel-T re-treatment (Supplementary Table S1).
Comparison of proportions of newly generated BCR clones present only at later time points (T2) versus the earlier time point (T1) between treated and naïve patients.
. | . | Proportion of newly generated BCR clones . | . | |
---|---|---|---|---|
. | . | Median (min–max) . | . | |
. | . | [N] . | . | |
T1 . | T2 . | Treated . | Naïve . | P . |
Week 0 | Week 2 | 0.74 (0.52–0.78) | 0.42 (0.06–0.86) | 0.021 |
[n = 6] | [n = 16] | |||
Week 0 | Week 4 | 0.71 (0.54–0.83) | 0.31 (0.03–0.74) | 0.002 |
[n = 6] | [n = 14] | |||
Week 0 | Week 6 | 0.57 (0.37–0.64) | 0.36 (0.04–0.76) | 0.109 |
[n = 5] | [n = 16] | |||
Week 0 | Week 26 | 0.55 (0.50–0.69) | 0.56 (0.10–0.73) | 0.871 |
[n = 4] | [n = 13] | |||
Week 0 | Week 52 | 0.65 (0.63–0.71) | 0.49 (0.04–0.82) | 0.24 |
[n = 4] | [n = 10] | |||
Week 2 | Week 4 | 0.48 (0.25–0.63) | 0.39 (0.21–0.87) | 0.689 |
[n = 7] | [n = 14] | |||
Week 2 | Week 6 | 0.35 (0.14–0.36) | 0.51 (0.02–0.89) | 0.098 |
[n = 5] | [n = 15] | |||
Week 2 | Week 26 | 0.40 (0.22–0.87) | 0.64 (0.27–0.85) | 0.383 |
[n = 5] | [n = 12] | |||
Week 2 | Week 52 | 0.45 (0.32–0.53) | 0.62 (0.13–0.91) | 0.298 |
[n = 5] | [n = 9] | |||
Week 4 | Week 6 | 0.24 (0.14–0.41) | 0.64 (0.15–0.92) | 0.007 |
[n = 5] | [n = 14] | |||
Week 4 | Week 26 | 0.38 (0.22–0.91) | 0.74 (0.32–0.86) | 0.246 |
[n = 5] | [n = 13] | |||
Week 4 | Week 52 | 0.33 (0.24–0.63) | 0.58 (0.17–0.91) | 0.129 |
[n = 5] | [n = 10] | |||
Week 6 | Week 26 | 0.51 (0.41–0.62) | 0.53 (0.04–0.78) | 0.956 |
[n = 4] | [n = 13] | |||
Week 6 | Week 52 | 0.55 (0.51–0.62) | 0.49 (0.08–0.88) | 0.839 |
[n = 4] | [n = 10] | |||
Week 26 | Week 52 | 0.54 (0.14–0.58) | 0.60 (0.07–0.72) | 0.414 |
[n = 4] | [n = 9] |
. | . | Proportion of newly generated BCR clones . | . | |
---|---|---|---|---|
. | . | Median (min–max) . | . | |
. | . | [N] . | . | |
T1 . | T2 . | Treated . | Naïve . | P . |
Week 0 | Week 2 | 0.74 (0.52–0.78) | 0.42 (0.06–0.86) | 0.021 |
[n = 6] | [n = 16] | |||
Week 0 | Week 4 | 0.71 (0.54–0.83) | 0.31 (0.03–0.74) | 0.002 |
[n = 6] | [n = 14] | |||
Week 0 | Week 6 | 0.57 (0.37–0.64) | 0.36 (0.04–0.76) | 0.109 |
[n = 5] | [n = 16] | |||
Week 0 | Week 26 | 0.55 (0.50–0.69) | 0.56 (0.10–0.73) | 0.871 |
[n = 4] | [n = 13] | |||
Week 0 | Week 52 | 0.65 (0.63–0.71) | 0.49 (0.04–0.82) | 0.24 |
[n = 4] | [n = 10] | |||
Week 2 | Week 4 | 0.48 (0.25–0.63) | 0.39 (0.21–0.87) | 0.689 |
[n = 7] | [n = 14] | |||
Week 2 | Week 6 | 0.35 (0.14–0.36) | 0.51 (0.02–0.89) | 0.098 |
[n = 5] | [n = 15] | |||
Week 2 | Week 26 | 0.40 (0.22–0.87) | 0.64 (0.27–0.85) | 0.383 |
[n = 5] | [n = 12] | |||
Week 2 | Week 52 | 0.45 (0.32–0.53) | 0.62 (0.13–0.91) | 0.298 |
[n = 5] | [n = 9] | |||
Week 4 | Week 6 | 0.24 (0.14–0.41) | 0.64 (0.15–0.92) | 0.007 |
[n = 5] | [n = 14] | |||
Week 4 | Week 26 | 0.38 (0.22–0.91) | 0.74 (0.32–0.86) | 0.246 |
[n = 5] | [n = 13] | |||
Week 4 | Week 52 | 0.33 (0.24–0.63) | 0.58 (0.17–0.91) | 0.129 |
[n = 5] | [n = 10] | |||
Week 6 | Week 26 | 0.51 (0.41–0.62) | 0.53 (0.04–0.78) | 0.956 |
[n = 4] | [n = 13] | |||
Week 6 | Week 52 | 0.55 (0.51–0.62) | 0.49 (0.08–0.88) | 0.839 |
[n = 4] | [n = 10] | |||
Week 26 | Week 52 | 0.54 (0.14–0.58) | 0.60 (0.07–0.72) | 0.414 |
[n = 4] | [n = 9] |
Note: The proportion of newly generated BCR clones was calculated by the number of the newly generated clones divided by the total number of clones presented at either time point. Bold text highlights significant differences between treated and naïve patients (Wilcoxon rank sum test P < 0.05).
Treated patients had many high-frequency clones (i.e., top 100 B-cell clonotypes at baseline) that persisted throughout treatment (Fig. 2D); for example, for the representative patient, 57 of the top 100 B-cell clonotypes at baseline persisted throughout treatment. In contrast, naïve patients lost the high-frequency clones that were present at baseline (Fig. 2E). The response patterns in the two groups appeared numerically different but were not statistically different (Wilcoxon rank sum test P = 0.077 for ICC comparisons between the two studies). For the top 100 clones post-baseline (except week 26), the likelihood of retaining the top ranked clones at other time points was lower in naïve patients than in treated patients, as evidenced by the lower ICCs (Wilcoxon rank sum test P = 0.034 at week 2, P = 0.026 at week 4, P = 0.011 at week 6, and P = 0.028 at week 52), indicating greater clonal shuffling in naïve patients (Fig. 2F). Supplementary Figure S2A also presented the rank change of clonotypes starting from week 2 to show the gain of clonotypes for naïve patients. In addition, Morisita distances were lower in naïve patients than in treated patients, indicating a greater degree of change in the B-cell repertoire among naïve patients (Supplementary Fig. S2B).
BCR convergence was induced with treatment
Antigen-specific immune responses can drive BCRs to converge to similar amino acid sequences in the CDR3 domains that are derived from different nucleotide sequences. The BCR convergent frequency, which was calculated as the aggregate frequency of clones sharing an amino acid sequence with at least one other clone, was significantly higher for treated patients versus naïve patients at every time point (Wilcoxon rank sum test P ≤ 0.05). In addition, re-treatment increased convergent frequency between week 0 and both week 2 (Wilcoxon rank sum test P = 0.063) and week 4 (Wilcoxon rank sum test P = 0.031), a pattern that was not observed in naïve patients (Fig. 3A). On average, the number and size of BCR convergent groups based on social network analysis were larger in treated patients than in naïve ones (Wilcoxon rank sum test Fig. 3B, top, P < 0.002; and Fig. 3B, bottom, P < 0.02, respectively).
Frequency of BCR convergent families differed with prior treatment with sipuleucel-T. A, BCR convergent frequency at all time points for 7 treated (left boxplot) and 19 naïve (right boxplot) patients. For each patient, convergent frequency was calculated as the aggregate frequency of clones sharing an amino acid sequence with at least one other clone (20). The asterisk shows significant increase of clonality in treated patients compared to naïve patients (Wilcoxon rank sum test P < 0.05). B, Violin plots of the number of BCR convergent groups (top) and the size of BCR convergent groups (bottom) at each time point. For each patient at each time point, social network analysis was performed on full-length of BCR clonotype sequences using R packages: Ape (24) and igraph (25). A convergent group was defined as the cluster that included the clones with the distance less than or equal to 1 (allowing maximum of one base pair mutation among clone sequences sharing the same V-gene, J-gene, and CDR3 length). The number of convergent groups and the maximum group size within each network were used to describe the feature of the network. At each time point, left and right half violin plots plot the number of BCR convergent groups (top) or the maximum number of clones within each BCR convergent group (bottom) at the corresponding time for treated patients and naïve patients, respectively. Social network figures of V03 family across all repertories of the representative treated patient (C) and of the representative naïve patient (D). For each patient, pairwise distance matrix of each pair of full-length of BCR clonotype sequences was calculated on the basis of Levenshtein distance (R package: RecordLinkage; ref. 23) within each V gene family. Social network analysis was then performed. Each node represents a full-length of single clone coded by the first presenting time point, and the node size is attributed on the basis of the corresponding abundance. The nodes connected by lines were within the same convergent group that included the full-length of BCR clonotype sequences with the distance less than or equal to 1. All boxplots present 25% percentile, median, and 75% percentile of the data, with each dot representing each sample. All violin plots show distributions of the data smoothed by a kernel density estimator.
Frequency of BCR convergent families differed with prior treatment with sipuleucel-T. A, BCR convergent frequency at all time points for 7 treated (left boxplot) and 19 naïve (right boxplot) patients. For each patient, convergent frequency was calculated as the aggregate frequency of clones sharing an amino acid sequence with at least one other clone (20). The asterisk shows significant increase of clonality in treated patients compared to naïve patients (Wilcoxon rank sum test P < 0.05). B, Violin plots of the number of BCR convergent groups (top) and the size of BCR convergent groups (bottom) at each time point. For each patient at each time point, social network analysis was performed on full-length of BCR clonotype sequences using R packages: Ape (24) and igraph (25). A convergent group was defined as the cluster that included the clones with the distance less than or equal to 1 (allowing maximum of one base pair mutation among clone sequences sharing the same V-gene, J-gene, and CDR3 length). The number of convergent groups and the maximum group size within each network were used to describe the feature of the network. At each time point, left and right half violin plots plot the number of BCR convergent groups (top) or the maximum number of clones within each BCR convergent group (bottom) at the corresponding time for treated patients and naïve patients, respectively. Social network figures of V03 family across all repertories of the representative treated patient (C) and of the representative naïve patient (D). For each patient, pairwise distance matrix of each pair of full-length of BCR clonotype sequences was calculated on the basis of Levenshtein distance (R package: RecordLinkage; ref. 23) within each V gene family. Social network analysis was then performed. Each node represents a full-length of single clone coded by the first presenting time point, and the node size is attributed on the basis of the corresponding abundance. The nodes connected by lines were within the same convergent group that included the full-length of BCR clonotype sequences with the distance less than or equal to 1. All boxplots present 25% percentile, median, and 75% percentile of the data, with each dot representing each sample. All violin plots show distributions of the data smoothed by a kernel density estimator.
Social network analysis of the representative patients across all time points, where the node color was coded based on the initial present time point of the corresponding clone, also demonstrated that most convergent groups appeared earlier in treated patients (Fig. 3C; e.g., at week 0 and week 2) than in naïve patients (Fig. 3D; e.g., week 4 or later, i.e., after two or all three infusions). Overall, the number of the convergent groups that persisted across time points was always higher in treated versus naïve patients (Wilcoxon rank sum test P < 0.001; Fig. 4A). The number of the convergent groups that were induced from baseline to week 2 and persisted in the subsequent time points, was higher in treated versus naïve patients (Wilcoxon rank sum test P < 0.001; Fig. 4B). This was consistent with there being more affinity maturation in treated patients. For example, among 227 convergent groups in the representative treated patient, 38 of them spanned all six time points (Fig. 4C), and 48 of them were enriched after the first infusion (Fig. 4D). Figure 4E and F provide example groups from the representative treated patient, which were persistently present across all time points and induced after the first infusion, respectively. In the representative naïve patient, no convergent groups spanned over all time points and only one enriched after the first infusion (Fig. 4G). The phylogenetic trees show convergent groups with the member nucleotide sequences (Fig. 4E–G). The clonal relatedness and abundance changes across time points were also evident in Supplementary Fig. S3A–S3C. This relatedness of clones could have reflected affinity maturation through somatic hypermutation or clonal selection.
Patterns of the convergent groups of 7 treated and 19 naïve patients. A, Violin plots of the number of the persist convergent groups. If the clones within a convergent group presented at all six time points, then the convergent group was defined as a persistent convergent group. Violin plots exhibit the number of the convergent groups that spanned over all time points for the treated patients (left) and the naïve patients (right). The asterisk shows significant difference between treated patients and naïve patients (Wilcoxon rank sum test P < 0.05). B, Violin plots of the number of the induced convergent groups. If the clones within a convergent group presented at all last five time points except baseline, then the convergent group was defined as an induced convergent group. Violin plots exhibit the number of the convergent groups that were induced between baseline and week 2 and continue to exist at remaining time points for the treated patients (left) and the naïve patients (right). The asterisk shows significant difference between treated patients and naïve patients (Wilcoxon rank sum test P < 0.05). C and D, The convergent groups tracking plots. C, The convergent groups that persist across the time points for the representative treated patient. D, The induced convergent groups that persist across the time points for the representative treated patient. Each dot connected line represents a convergent group. Selected convergent groups in the representative treated patient (E and F) and in the representative naïve patient (G). These plots present the groups that persist across all time points (E) and the induced groups for the treated patient (F). The induced groups for the naive patient are presented in G. In each figure, the left panel shows the network for a representative group, with the node color coded by the first presenting time point, the node size reflecting the corresponding abundance, and the lines linking clones with only one base pair difference. The adjacent tree plot presents the presence of the clone across six time points, with filled circles representing the clone's presence at the corresponding time point, while open circles denote time points where the clone is absent. The circles or dots in the same row stand for the exact same clones, which correspond to the actual nucleotide sequence in the same row colored coded by ACTG (only first 60 base pairs were retained due to limit space). The yellow highlights present the hypermutation in the nucleotide sequence, which is also corresponding to the bottom row of the motif analysis. The right most panels are chord diagrams, where the outside circle color stands for the time point (denoted as color legend) and the width of the arch stands for the abundance information (log10 transformed). The same inside color connects the same clone within each chord diagram. All violin plots show distributions of the data smoothed by a kernel density estimator.
Patterns of the convergent groups of 7 treated and 19 naïve patients. A, Violin plots of the number of the persist convergent groups. If the clones within a convergent group presented at all six time points, then the convergent group was defined as a persistent convergent group. Violin plots exhibit the number of the convergent groups that spanned over all time points for the treated patients (left) and the naïve patients (right). The asterisk shows significant difference between treated patients and naïve patients (Wilcoxon rank sum test P < 0.05). B, Violin plots of the number of the induced convergent groups. If the clones within a convergent group presented at all last five time points except baseline, then the convergent group was defined as an induced convergent group. Violin plots exhibit the number of the convergent groups that were induced between baseline and week 2 and continue to exist at remaining time points for the treated patients (left) and the naïve patients (right). The asterisk shows significant difference between treated patients and naïve patients (Wilcoxon rank sum test P < 0.05). C and D, The convergent groups tracking plots. C, The convergent groups that persist across the time points for the representative treated patient. D, The induced convergent groups that persist across the time points for the representative treated patient. Each dot connected line represents a convergent group. Selected convergent groups in the representative treated patient (E and F) and in the representative naïve patient (G). These plots present the groups that persist across all time points (E) and the induced groups for the treated patient (F). The induced groups for the naive patient are presented in G. In each figure, the left panel shows the network for a representative group, with the node color coded by the first presenting time point, the node size reflecting the corresponding abundance, and the lines linking clones with only one base pair difference. The adjacent tree plot presents the presence of the clone across six time points, with filled circles representing the clone's presence at the corresponding time point, while open circles denote time points where the clone is absent. The circles or dots in the same row stand for the exact same clones, which correspond to the actual nucleotide sequence in the same row colored coded by ACTG (only first 60 base pairs were retained due to limit space). The yellow highlights present the hypermutation in the nucleotide sequence, which is also corresponding to the bottom row of the motif analysis. The right most panels are chord diagrams, where the outside circle color stands for the time point (denoted as color legend) and the width of the arch stands for the abundance information (log10 transformed). The same inside color connects the same clone within each chord diagram. All violin plots show distributions of the data smoothed by a kernel density estimator.
Most related groups converged to single or only few amino acid sequences, for example, the phylogenetic trees of 2 representative patients across all V gene families (Supplementary Fig. S4A; Fig. 4B), where the edge color was coded by the corresponding amino acid sequence. For the treated patient, among 227 convergent groups, 102 groups had identified amino acid sequences (Supplementary Fig. S4C), each member of the group converged to a single amino acid sequence or to two amino acid sequences. Supplementary Figure S4D presents the top six largest groups color coded by amino acid sequences, which all converged to one or two amino acid sequences. In addition, we found that there were five public clonotypes (Supplementary Table S2), for which each was shared by 2 of 6 treated patients at baseline (please note there were only 6 treated patients having baseline BCR data) and 2 of naïve patients post-baseline (i.e., week 2 to week 52). They corresponded to amino acid sequences of CAKDLDYW, CAKEDYW, CARDFDYW, CARDLDYW, CARDSSGWYYFDYW. On the basis of these public clones, five public clusters were identified that were shared by treated patients at baseline and naïve patients post-baseline. (Supplementary Fig. S4E).
We also performed a sensitivity analysis to check whether clonal grouping can be impacted by sequencing depth by subsampling (Supplementary Fig. S5A). The representative treated patient had a more complex network (Supplementary Fig. S5B) compared to the naïve patient (Supplementary Fig. S5C), and the network analysis based on subsampling of the treated patient retained a similar pattern (Supplementary Fig. S5D).
Focused BCR repertoire associated with survival
The median OS for the 19 naïve patients with BCR data was 36 months. Although at week 0, a more focused BCR repertoire was not correlated with longer survival (Fig. 5A), at week 6 (after all three infusions), a more focused BCR repertoire, that is, higher clonality (>0.1, which was the median clonality across all time points), correlated with longer survival (log rank test P = 0.045; Fig. 5B). In addition, ranks of the top 100 abundant clones at each week were compared across time points within patients with OS > 36 months and patients with OS ≤ 36 months, respectively. Higher ICCs were only seen in patients with OS ≤ 36 months at week 2 (Wilcoxon rank sum test P = 0.032) and week 4 (Wilcoxon rank sum test P = 0.063), indicating less clonal shuffling in these top 100 clones at week 2 and week 4 for patients with OS ≤36 months, but not in patients with OS > 36 months. This suggested that fewer changes in the repertoire induced by treatment were associated with shorter survival.
Association between BCR clonality and clinical outcome for 19 naïve patients. A, OS is illustrated in Kaplan–Meier curves for patients with low clonality (<0.1) and high clonality (>0.1) at week 0. For each repertoire, clonality was calculated as 1 − H/loge(n), where H is Shannon entropy index and n is the number of unique clones (19). Clonality of 0 and 1 indicate total diversity and no diversity. The median clonality was 0.1 across all time points for all naïve patients. Median OS was obtained for the patients with clonality at week 0 < 0.1 versus > 0.1, respectively, by Kaplan–Meier method. Median OS did not differ between patients with low clonality and those with high clonality. B, OS is illustrated using Kaplan–Meier curves for patients with low clonality (<0.1) and high clonality (>0.1) at week 6. Median OS was obtained for the patients with clonality at week 6 < 0.1 versus > 0.1, respectively, by Kaplan–Meier method. Patients with low clonality had a short median OS of 26.6 months versus 35.6 months for the patients with high clonality (log-rank test P = 0.045).
Association between BCR clonality and clinical outcome for 19 naïve patients. A, OS is illustrated in Kaplan–Meier curves for patients with low clonality (<0.1) and high clonality (>0.1) at week 0. For each repertoire, clonality was calculated as 1 − H/loge(n), where H is Shannon entropy index and n is the number of unique clones (19). Clonality of 0 and 1 indicate total diversity and no diversity. The median clonality was 0.1 across all time points for all naïve patients. Median OS was obtained for the patients with clonality at week 0 < 0.1 versus > 0.1, respectively, by Kaplan–Meier method. Median OS did not differ between patients with low clonality and those with high clonality. B, OS is illustrated using Kaplan–Meier curves for patients with low clonality (<0.1) and high clonality (>0.1) at week 6. Median OS was obtained for the patients with clonality at week 6 < 0.1 versus > 0.1, respectively, by Kaplan–Meier method. Patients with low clonality had a short median OS of 26.6 months versus 35.6 months for the patients with high clonality (log-rank test P = 0.045).
Discussion
Immunotherapies are being used to treat malignancies, yet the durability of the antigen-specific responses induced by these treatments and their impact on the adaptive immune system remains unclear. Most studies have investigated how the T-cell repertoire is modulated by immunotherapies and how the T-cell repertoire affects clinical response. There are emerging studies suggesting that the B-cell repertoire may in fact be important in prognosticating and predicting survival and response to immunotherapy (29, 30). In addition, the B-cell repertoire can mediate response to checkpoint inhibitors (31, 32). Radiotherapy and checkpoint inhibitor immune therapies for various cancers induce B-cell somatic hypermutation, which may lead to durable clinical response (33, 34). However, overall, the B-cell arm of the adaptive immune system remains largely unstudied in the context of cancer immunotherapies. This is to our knowledge the first study to look at the effects that sipuleucel-T, an autologous immunotherapy, has on the B-cell repertoire and the durability of antigen-specific response.
Sipuleucel-T improves OS in men with mCRPC, along with delayed opioid use and delayed use of secondary anticancer agent use (3). Antigen-specific immune responses induced by sipuleucel-T correlate with survival benefit (3). Sipuleucel-T induces B-cell immune responses to PA2024, the immunogen used to manufacture the cell product. To evaluate antigen specificity, we also assessed antibody responses to PSA, another prostate cancer–associated antigen not contained in sipuleucel-T. The lower anti-PSA responses seen in the naïve group may reflect baseline differences in immune status between the groups including the possibility that the prior sipuleucel-T–induced antigen spreading to encompass PSA. The level of antibodies to PA2024 was significantly higher in treated patients than in naïve patients. This was not seen for PSA, a nonspecific antigen, which suggested that treated patients developed antigen-specific immune response induced by sipuleucel-T and that this antigen specificity persisted over time. We found that there were five public clones shared by treated patients at baseline and naïve patients across the five post-baseline time points. Because those clones were not present at baseline, but were induced following treatment for the naïve patients, and they were already present at baseline for the treated patients, they might have defined memory B cells with potential specificity to PAP. Some of these preexisting public clonotypes in treated patients disappeared with re-treatment. This could have represented the trafficking of these B cells out of the circulation with stimulation, as one of these clonotypes did return at later time points. Alternatively, clonotypes may have undergone further somatic hypermutation so that they no longer had this phenotype.
Despite a median 8.9-year interval between the initial treatment and subsequent treatment with sipuleucel-T, treated patients exhibited a more focused BCR repertoire, higher convergent frequency, and distinct BCR gene usage before re-treatment compared with naïve patients at baseline. Both social network and phylogenetic tree analysis indicated that treated patients exhibited more convergent groups with more clones in each convergent group, and earlier established group patterns (i.e., at week 0 and week 2) compared with naïve patients; the latter formed group patterns at later time points (week 6 and week 26). These changes were consistent with clonal maturation that was driven by antigen recognition, including endogenous tumor antigens. The same amino acid sequences were derived by different new nucleotide sequences, consistent with an antigen-driven response, most likely to PAP. Some naïve patients had consistently lower clonality even at week 52 compared with week 0 of treated patients. There are many factors that could influence the complexity of an individual's B-cell repertoire including history of infections, vaccinations, and/or treatments. In this vein, this might have been explained by that BCR repertoire of P10-1 (treated) patients having been prealtered by the previous sipuleucel-T treatment; the prior treatment with sipuleucel-T could have left a permanent imprint on the repertoire with some durably expanded clonotypes and persistence of these clonotypes resulting in this higher clonality.
Our results also provide novel insight into memory B-cell recall responses induced with immunotherapy. Memory B-cell recall responses to annual influenza vaccinations do not have overlapping BCR sequences between individuals, which may be due to the low chances of different individuals generating the same antibody responses (35). However, BCR sequences may be shared between different time points within the same individual. Here, treated patients had a higher proportion of clones that persisted during and after re-treatment. In one of the treated patients, more than 16% of the convergent families (38 of 227) spanned all time points. These results were consistent with these clones being primed with initial treatment, surviving during the intervening years, and persisting with re-treatment.
Treated patients demonstrated an immunologic boost effect upon sipuleucel-T re-treatment. These patients had higher BCR clonality and more rapid changes with re-treatment that persisted until week 52 or later after treatment. Pairwise Morisita distance was consistently higher in treated patients than in naïve patients, indicating less clonal shuffling in treated patients. This finding supported the notion that naïve patients were generating de novo responses (i.e., immune priming). Consistent with this, we observed notable clonotype shuffling of the most abundant B-cell clonotypes in naïve patients. This demonstrated that re-treatment with sipuleucel-T resulted in a more rapid focusing of the BCR repertoire and rapid expansion of treatment-related antigen-specific subset of clones.
Gene expression of certain BCR gene segments correlate with overall survival and progression-free survival of patients in breast cancer and ovarian cancer (36). This finding is in-line with the changes in BCR repertoire in our data being associated with overall survival. In naïve patients, higher clonality in the BCR repertoire induced by sipuleucel-T (after all three infusions) was correlated with longer survival. Whether differences in response reflect immune-intrinsic or tumor-intrinsic determinants is unknown. If a patient does not experience a sufficient increase in clonality with standard three doses, then perhaps that patient may benefit from additional treatment with sipuleucel-T.
There are some limitations related to these analyses. First, the number of patients studied was small as only a limited number of samples were available from the two studies to conduct these ad hoc experiments. The two clinical studies were originally designed to assess immune responses with treatment in each disease context, but not for the B-cell response experiments described herein. As a result, we relied on available samples for these exploratory analyses. However, it is rare to have samples from patients who had previously received the same immunotherapy, which we think represents a critically novel element of the current work. Second, the association between overall survival and BCR clonality was weak probably due to the smaller sample size. The presented analysis was exploratory and will need further substantiation in follow-on studies with larger number of subjects. Third, because sipuleucel-T induces T- and B-cell responses to PAP that correlate with improved survival, it would be interesting to assess whether these responses are parallel and complementary in patients by combined evaluations of T- and B-cell responses, which would be examined in future studies.
In conclusion, sipuleucel-T treatment altered the BCR repertoire and resulted in long-term immunologic memory. In addition, sipuleucel-T re-treatment resulted in an immunologic boost effect. Analysis of the BCR for motifs and single-cell studies could help to further elucidate the functional states and the antigens recognized by the different B-cell clones. Nevertheless, these results demonstrated that cancer immunotherapies can induce memory responses lasting years.
Disclosure of Potential Conflicts of Interest
L. Zhang reports personal fees from Dendreon Pharmaceuticals LLC (paid consultant) during the conduct of the study, as well as personal fees from Smith-Kettlewell Eye Research Institute (paid consultant) and personal fees and other from Raydiant Oximetry, Inc. (paid consultant) outside the submitted work. D.Y. Oh reports grants from Prostate Cancer Foundation (young investigator award) during the conduct of the study, as well as other from Roche/Genentech (research support) and Merck (research support) and personal fees from Maze Therapeutics (consulting) outside the submitted work. N.A. Sheikh reports personal fees from Dendreon Pharmaceuticals LLC (salary) outside the submitted work. L. Fong reports grants and personal fees from Dendreon Pharmaceuticals LLC during the conduct of the study, as well as grants from AbbVie, Bavarian Nordic, Janssen, and Merck and grants and personal fees from Bristol-Myers Squibb and Roche/Genentech outside the submitted work.
Authors' Contributions
L. Zhang: Data curation, software, formal analysis, supervision, validation, investigation, visualization, methodology, writing–original draft, project administration, writing–review and editing. H. Kandadi: Resources, data curation, software, formal analysis, validation, investigation, visualization, methodology, writing–original draft, writing–review and editing. H. Yang: Software, visualization, methodology. J. Cham: Investigation, writing–original draft, writing–review and editing. T. He: Software, visualization, methodology, writing–review and editing. D.Y. Oh: Investigation, writing–review and editing. N.A. Sheikh: Conceptualization, resources, data curation, supervision, funding acquisition, validation, investigation, writing–original draft, project administration, writing–review and editing. L. Fong: Conceptualization, resources, supervision, funding acquisition, investigation, writing–original draft, project administration, writing–review and editing.
Acknowledgments
Editorial assistance was provided by Helen M. Wilfehrt, PhD, ISMPP CMPP of Dendreon Pharmaceuticals LLC. D.Y. Oh is supported by NIH 4T32 CA177555, 1K08 AI139375, the Harry F. Bisel, MD, Endowed Young Investigator Award from the Conquer Cancer Foundation of the American Society of Clinical Oncology, and the Prostate Cancer Foundation Young Investigator Award. L. Fong is supported by NIH R01CA223484, U01CA233100, and the Prostate Cancer Foundation. L. Zhang receives the support of the University of California San Francisco Academic Senate Committee on Research.
The costs of publication of this article were defrayed in part by the payment of page charges. This article must therefore be hereby marked advertisement in accordance with 18 U.S.C. Section 1734 solely to indicate this fact.