3.1 Patient Population
A total of 60 patients were enrolled in the study between August 2018 and May 2021. Two patients were excluded. One of them received neoadjuvant letrozole during the SarsCov2 pandemic period due to a delay in the scheduled date of surgery, and other patient only received one dose of denosumab due to withdrawal of consent after the first infusion. One patient who underwent all the procedures but did not receive denosumab due to an administrative error was transferred to the control group. Thus, a total of 37 patients were analyzed in the experimental group and 21 patients in the control group, all 58 patients completed the follow-up period (Figure 2). The analysis was done by protocol.
The clinicopathological baseline characteristics of the evaluable patients are shown in Table 1. The mean age of the study population was 56.4 years (range, 37-80 years), with a mean age of 57 years in the experimental arm and 55.4 years in the control arm. The mean time between the first administration of denosumab and surgery was 21 days. Of the total participants, 27 were pre-menopausal women, with 17 (45.9%) in the experimental arm and 10 (47.6%) in the control group, meeting the requirement of including at least 24 pre-menopausal women. The distribution of post-menopausal patients was: 20 (54.1%) in the experimental arm and 11 (52.4%) in the control group. Most patients were classified as IA clinical stage, with 26 cases (70.3%) in theT experimental arm and 17 cases (81%) in the control arm. The median tumor size was 18 mm (range 8-45 mm).
In terms of tumor characteristics, the majority of patients had invasive breast cancer without specific features, classified as no special type (NST) or not otherwise specified (NOS) (ductal), 62.07% overall; with a slightly higher percentage in the control arm (71.4% vs. 56.8% in the experimental group), with a lower number of patients with invasive lobular carcinoma in the control arm (14.3 vs. 32.4% in the experimental arm), although these differences were not significant. Most tumors had a histological grade 2 (58.62%), with 7 tumors (12.06%) classified as grade 3; notably, there was a higher percentage of cases with histological grade 1 tumors (35.1%) in the experimental group compared to the control group (19%). The data are consistent with a higher number of cases with low Ki67 (<15) in the experimental group (43.2%) compared to the control group (9.5%). Ki67 was the only parameter that showed significance (p=0.001) when comparing both groups; however, when analyzing the tumors by surrogate molecular subtype, this difference was not significant, although numerically, there were still more cases of luminal A-like tumors in the experimental arm (51.4% vs. 38.1%). The groups were well balanced as no statistically significant differences between experimental and control arms were observed, except for the percentage of patients with tumors with low Ki67 (<15%) (Table 1). Of the total of 58 evaluable cases, 48 were luminal tumors (27 luminal A-like and 21 luminal B-like), while 10 were TNBC. Initially, a higher percentage of TNBC cases was expected. However, most TNBC tumors were selected for neoadjuvant treatment and were therefore excluded from our study. Upon reviewing these 10 TNBC cases, only 5 exhibited typical aggressive characteristics (histological grade 2/3, Ki67 >30%), with 3 of these cases in the control group and 2 in the experimental arm. The remaining cases were low aggressive TNBC, including 3 cases of invasive carcinomas with apocrine differentiation, 2 in the experimental arm and 1 in the control arm: apocrine carcinomas are known to exhibit a more indolent behavior compared to typical TNBCs. Additionally, the experimental group had 1 case of lobular carcinoma and 1 case of carcinoma not otherwise specified (NOS) ductal, with Ki67 of 10% and 5%, respectively. Given these factors, caution should be exercised in drawing conclusions about the triple negative subgroup.
3.2 Safety Data
Five patients experienced localized hematomas in the breast following study biopsies, which were classified as grade 1 and did not require drainage or special measures; 2 of these cases were in the control group and 3 were in the experimental group. In the experimental arm, the most common adverse events were grade 1 or 2 bone pain occurring 24 hours after infusion in 10 of 37 patients (27.03%), grade 1 asthenia in 4 patients (10.81%), grade 1 pain at the denosumab infusion site in 3 patients (8.10%), grade 1 chills in 2 patients (5.41%), and grade 2 dental infection in 1 patient (2.70%). There were no reported cases of hypocalcemia. No grade 3 toxicities were reported (see Table S1 in the Supplementary Appendix).
Outside the study follow-up period, we observed a long-term event of osteonecrosis of the jaw in a heavy-smoker patient who received denosumab on March 10th and 17th, 2021. This patient was initially reported with a dental infection. Symptoms began one month after the last dose of denosumab (reported on April 19, 2021), with discomfort and pain in the jaw, the case was referred to the Maxillofacial Surgery Department. Computed tomography initially did not reveal signs of osteonecrosis, and the case was initially diagnosed as a tooth infection, which improved with oral antibiotic treatment. However, repeated episodes of dental infection in the same location needed specific follow-up. After 11 months, the patient underwent jaw surgery in February 2022, and the pathological report confirmed osteonecrosis of the right quadrant 44-47. This event was considered possibly related to denosumab, although other triggering risk factors such as chronic infection due to long-term smoking should be considered. The osteonecrosis of the jaw was classified as grade 3. At the time of this report, no grade 4 or 5 toxicity has been reported.
3.3 Serum Analysis
Denosumab was associated with systemic inhibition of RANKL but not with changes in bone remodeling markers
The blockade of the RANK-RANKL pathway was confirmed by the drop in serum of free RANKL (sRANKL), measured by ELISA, in the experimental group as RANKL became bound to denosumab (mean serum A 0,096 pg./L vs serum B vs. 0,000 pg./L p<0.001), while no changes were found in the control group (mean serum A 0,100 pg./L vs. serum B 0,116 pg./L; p=0.270) comparing serum A vs serum B. OPG levels tended to increase in the experimental group (p=0.071), consistent with the reduction of free RANKL. However, the serum levels of the bone resorption markers tartrate-resistant acid phosphatase 5b (TRACP5b) (n=37) and carboxy-terminal collagen crosslinks (CTX) (n=38) did not change in any group (Figure 3 and Supplementary data Figure S1). It is unclear whether the lack of changes is due to technical limitations (limited detection rate) or due to the kinetics of bone resorption. Although no significant alterations were observed in bone resorption markers, a small decrease in serum calcium -not clinically relevant- was reported in the experimental arm (Figure 3), despite the prescribed calcium and vitamin D supplementation. This finding reinforces the efficacy of denosumab in bone remodeling.
A correlation analysis was conducted between the different serum markers studied and menopausal status to better understand the biology of the pathway. Follicle-stimulating hormone (FSH) levels were higher in post-menopausal patients, consistent with menopausal physiology (p<0.0001). Levels of free sRANKL did not differ based on menopausal status, while higher levels of OPG were detected in post-menopausal women (p=0.010) at baseline (Supplementary data Figure S2A). Comparable values of the bone markers TRACP5b, and CTX were found at the time of diagnosis (baseline) between pre- and post-menopausal women, although slightly higher levels were found in the post-menopausal group. In pre-menopausal patients, no associations were found between levels of progesterone and sRANKL in serum (p=0.401). Correlation between OPG and TRACP5b values was not significant, unless one sample with high OPG and low TRACP5b was excluded. Finally, a negative correlation between OPG and sRANKL was demonstrated (p=0.0026), consistent with the known interaction between these factors (Supplementary data Figure S2B).
3.4 Tumor assessment´s results
Denosumab was not associated with a reduction in tumor cell proliferation or an increase in apoptosis
The primary endpoints of the clinical trial were a decrease in tumor cell Ki67 and an increase in apoptosis between biopsy and surgery. Denosumab did not reduce tumor cell proliferation or survival between paired biopsy and surgery samples (Figure 4A-B). The percentage of tumor cells expressing Ki67 increased in both groups (control p=0.035 and experimental p=0.012), which may be attributed to a higher quantification of fields within the surgical specimen (more fields) compared to the core biopsy. The mean Ki67 in the control arm increased from 24.52% at biopsy to 29.19% at surgery. Similarly, in the experimental arm, it increased from 20.86% to 24.81%. This nearly 5-percentage-point increase when comparing surgical samples with baseline biopsies in both groups suggests that these changes were not influenced by denosumab treatment. Indeed, a comparison between the experimental and control groups revealed identical behavior (p=0.928) (Figure 4A).
Denosumab did not induce an increase in tumor cell apoptosis, as demonstrated by the assessment of the H-score of cleaved caspase-3 between biopsy and surgery. The comparison between the experimental and control groups revealed notable inter-patient variability. Although a statistically significant difference in apoptosis was observed inter-group (p=0.042), the change in cleaved caspase-3 H-score quantification was less than 1 in both groups. Despite the contrasting trends, the alteration in apoptosis is considered clinically insignificant due to very low H-scores in all cases. The evaluating intra-patient (Paired T-test) showed no changes (control group p=0.060 and experimental p=0.238). Additionally, 3 patients in the experimental group exhibited higher levels of cleaved caspase-3 at baseline (Figure 4B) that may suggest potential deterioration or non-specific staining in some areas.
A subgroup analysis was conducted to elucidate if any specific patient group could benefit from denosumab treatment. Patients were divided according to surrogate molecular subtype (Supplementary data Table S2). Patients with luminal A-like (n=27) and luminal B-like (n=21) tumors showed similar trends to the overall population: there was no reduction in Ki67 or increase in Cleaved Caspase-3 in the experimental group and no changes were observed in OPG or TRACP5b. In the TNBC group (n=10), there was an increase in Ki67 in the experimental arm (p=0.025), which was not evident in the control group (p=0.517), but no difference was found in the inter-group comparison (p=0.197). Moreover, it is important to note the small number of TNBC cases, with 6 in the experimental group and 4 in the control group, as well as an imbalance between aggressive tumors, as previously explained (Supplementary data Table S2).
A subgroup analysis was also conducted based on menopausal status, excluding cases of triple-negative tumors to avoid biasing the information, as 9 out of 10 TNBC tumors were post-menopausal. Both pre-menopausal (n=26) and post-menopausal (n=22) tumors showed similar trends, with no reduction in proliferation or cell survival, and no other notable findings (Supplementary data Table S3).
Denosumab increased tumor infiltrating lymphocytes in early breast cancer, particularly in luminal B-Like tumors and regardless of menopausal status
Next, we interrogated the effect of denosumab on tumor immune infiltration and observed an increase in TILs in the surgery sample compared to the initial biopsy in the experimental arm (p=0.001), but not in the control arm (p=0.06). However, both groups showed a similar trend (p=0.789) (Figure 5A-B).
Performing a similar analysis based on surrogate molecular subtype, we found that while no changes in TILs were noted in Luminal A tumors (experimental (p=0.144), control (p=0.958)), there was a denosumab-induced elevation in TILs in luminal B-like tumors (p=0.012), with no significant changes in the control group (p=0.070). TILs did not change in TNBC (experimental (p=0.079), control (p=0.237)), although an increased number of representative TNBC tumors is required to reach conclusions in this subtype. Attending to menopausal stage in luminal tumors, we found that both pre-menopausal and postmenopausal patients experienced an increase in TILs in the experimental arm (premenopausal: p=0.048, postmenopausal: p=0.041), but not in the control arm (premenopausal: p=0.639, postmenopausal: p=0.062). In all comparisons the experimental and control arms showed similar trends (Table 2). Furthermore, applying a threshold of a 10% or greater increase in TILs between biopsy and surgery samples, we observed that 9 out of 37 (24.3%) patients in the experimental arm and 5 out of 21 (23.8%) patients in the control arm showed a clinically significant increase in TILs.
Tumor and stroma RANK expression was associated with highly proliferative tumors
The quantification conducted externally by NeoGenomics and in-house by the laboratory of Dr. Gonzalez-Suarez showed perfect correlation in the analysis of tumor RANK and RANKL protein expression (Supplementary data Figure S3). Given the correlation and the fact that the quantification of H-score for RANK and RANKL in the stroma was only conducted in the laboratory of Dr. Gonzalez-Suarez, all subsequent analyses were carried out using in-house quantification.
A total of 55 cases were assessable for tumor RANK and RANKL expression at baseline (Table 3, Figure 6A-B) and 19 tumors (34.54%) exhibited positive baseline expression of RANK, defined as an H-score > 0 (tumor RANK+). Of these positive cases, 14 tumors were randomized to the experimental arm (38.9%), and 5 were assigned to the control group (p=0.526, well balanced). The frequency of RANK+ tumors was comparable between tumors from premenopausal (36%) and postmenopausal patients (33%). Additionally, we compared tumor RANK expression across different molecular subtypes, 28% of luminal A-like, 35% of luminal –B-like and 50% of triple negative tumors exhibited RANK+ tumor cells, a higher frequency than that previously reported, particularly in luminal tumors (18,19). A total of 42.4% of grade 2 and 42.9% of grade 3 tumors were positive for tumor RANK, compared to only 13.3% of grade 1 tumors.
Regarding RANKL expression, 17 out of 55 tumors (30.10%) exhibited RANKL expression in tumor cells, with 12 assigned to the experimental group and 5 to the control group (p=0.819, well balanced). The frequency of tumor RANKL expression was similar between pre- and post-menopausal conditions (32% and 30%, respectively). No differences were found concerning molecular subtype or histological grade (Table 3). In only 8 samples (14.55%), both RANK+ and RANKL+ tumor cells were identified, but in these overlapping cases H-scores were low.
Based on the assessment of RANK and RANKL expression in stromal cells, 27 out of 56 evaluable samples (48,21%) and 18 out of 55 cases (32.7%) respectively exhibited an H-Score>0. Both the control and experimental groups were well balanced at baseline. There were no differences based on menopausal status; although in the premenopausal group stromal RANK expression was found in 57.7% compared to 40% in postmenopausal. Upon analysis by molecular subtype, it was noteworthy that luminal B-like tumors exhibited elevated RANK expression in the stroma, with rates of 75% compared to 34.5% in luminal A-like tumors and 30% in TNBC. This difference was statistically significant (p=0.012), representing the sole parameter where such distinction was observed (Table 3).
Observing this pattern and recognizing that RANK and RANKL expressions behave more like continuous than categorical variables, we decided to perform an analysis between Ki67, histological grade, ER expression and TILs at the level of basal biopsies (to avoid deviations related to denosumab) and correlate them with the IHC expression of RANK and RANKL in these biopsies (Figure 6C-F). A positive correlation was identified between RANK expression in tumor cells and cell proliferation (Ki67) p=0.03 and histological grade p=0.015, while a negative correlation with ER expression p=0.006 and no association with % of TILs was observed (Figure 6C). The findings indicated that cases characterized by high histological grade, high Ki67 levels, and low estrogen receptor expression showed higher expression of RANK protein in tumor cells (Figure 6C). Despite tumors with the highest RANKL scores (H>25) showed low levels of ki67, RANKL expression in tumor cells did not associate with Ki67, nor with the other parameters (Figure 6D). Strikingly, RANK expression in the stroma was associated with high Ki67 expression p=0.001, while stromal RANKL did not associate with any parameter (Figure 6E-F). The percentage of TILs at baseline did not associate with RANK or RANKL expression. The notable correlation observed in between RANK expression in both tumor and stroma underlines its association with aggressive tumors.
When restricting the analyses to luminal tumors neither RANK nor RANKL expression in tumor cells associated with any of the parameters analyzed (Supplementary data Figure S4A-B). Notably, increased stromal RANK expression remained associated with high Ki67 levels (p=0.001) and a positive correlation between RANKL in the stroma and TILs was observed (possibly as a marker of this population) (Supplementary data Figure S4 C-D). Therefore, we can infer that higher tumor proliferation correlates with higher RANK expression in the stroma in luminal tumors and the global analyses including TNBC. These results provide valuable information on the relationship between RANK, RANKL and other BC markers.
As expected, no changes in the expression of RANK or RANKL in tumor cells or in the stroma were found between biopsy and surgery, neither in the control nor in the experimental arm (Supplementary data Figure S5A-D).
Baseline RANK or RANKL expression did not predict denosumab-driven changes in TILs
Univariate and multivariate analyses were performed to identify possible factors associated with the 10% increase in TILs. Only having a high Ki67 (>30) could be related to an elevation in TILS, despite no reaching significance in the multivariate analyses (OR 7.12 (1.18-43.1) (p=0.079) (Table 4). The only significant factor identified in the multivariate analysis that correlated with an increase in TILs was RANK expression in tumor cells at baseline (p < 0.001). No other factors were found to be significantly associated with the 10% increase in TILs (Table 4).
Finally, when we analyzed the expression of RANK and RANKL as possible biomarkers of response to denosumab, we observed that the trends were similar to the overall population (Tables 5 and 6). There was no reduction in Ki67 or increase in cleaved caspase-3 after denosumab treatment when the analyses were performed only in tumors expressing RANK or RANKL protein in the tumor or stroma at baseline. Importantly, denosumab increased TILs regardless of tumor and stroma RANKL or RANK expression. This is, in tumor RANK positive samples, the experimental group showed an increase in TILs (p=0.013), similar to that observed in tumors that did not express RANK (tumor RANK-) (p=0.008). The same was observed when tumor RANKL expression was considered, the tumor RANKL positive group showed an increase in TILs (p=0.048), similar to the tumor RANKL negative group (p=0.002) (Table 6). The benefit of increased TILs after denosumab treatment was observed in tumors, irrespectively of the stromal expression of RANK or RANKL (Table 7). In conclusion, RANK and RANKL protein expression at baseline cannot be used as a biomarker capable of predicting the elevation of TILs caused by denosumab.